 Methodology
 Open Access
Identifying large sets of unrelated individuals and unrelated markers
 Kuruvilla Joseph Abraham^{1, 3}Email author and
 Clara Diaz^{2}
https://doi.org/10.1186/1751047396
© Abraham and Diaz; licensee BioMed Central Ltd. 2014
 Received: 15 July 2013
 Accepted: 24 February 2014
 Published: 17 March 2014
Abstract
Background
Genetic Analyses in large sample populations are important for a better understanding of the variation between populations, for designing conservation programs, for detecting rare mutations which may be risk factors for a variety of diseases, among other reasons. However these analyses frequently assume that the participating individuals or animals are mutually unrelated which may not be the case in large samples, leading to erroneous conclusions. In order to retain as much data as possible while minimizing the risk of false positives it is useful to identify a large subset of relatively unrelated individuals in the population. This can be done using a heuristic for finding a large set of independent of nodes in an undirected graph. We describe a fast randomized heuristic for this purpose. The same methodology can also be used for identifying a suitable set of markers for analyzing population stratification, and other instances where a rapid heuristic for maximal independent sets in large graphs is needed.
Results
We present FastIndep, a fast random heuristic algorithm for finding a maximal independent set of nodes in an arbitrary undirected graph along with an efficient implementation in C++. On a 64 bit Linux or MacOS platform the execution time is a few minutes, even with a graph of several thousand nodes. The algorithm can discover multiple solutions of the same cardinality. FastIndep can be used to discover unlinked markers, and unrelated individuals in populations.
Conclusions
The methods presented here provide a quick and efficient method for identifying sets of unrelated individuals in large populations and unlinked markers in marker panels. The C++ source code and instructions along with utilities for generating the input files in the appropriate format are available at http://taurus.ansci.iastate.edu/wiki/people/jabr/Joseph_Abraham.html
Keywords
 Linkage Disequilibrium
 Undirected Graph
 Unrelated Individual
 Sparse Graph
 Unlinked Marker
Background
The analysis of genotypes collected from large numbers of individuals is necessary for Genome Wide Asociation Studies (GWAS), conservation biology, and population genetics among other purposes. Most statistical analyses in these areas assume the individuals are unrelated, which may not always be the case in large sample populations. One way to avoid false positives due to the presence of related individuals is to identify the largest subset of unrelated individuals in the study population and retain only those individuals. A similar computational issue arises when identifying sets of unrelated markers for use in analyzing population stratification. Analyzing population stratification with the dense marker maps which are now currently available can be very time consuming due to the large number of markers being used, however due to linkage disequilibrium between the markers in dense maps not all the information provided by the markers can be considered independent. In this situation it makes sense to select a subset of markers which are mutually independent and which can provide sufficient information to analyze population stratification; this strategy reduces the computational burden arising from using all available markers while attempting to retain as many markers as possible. Indeed, certain well established approaches to analyzing population stratification assume that the input markers are unrelated [1], thus there are both conceptual and computational reasons for using a subset of weakly correlated markers. However, in populations with a complex history, identifying this set of markers may not be possible using just genetic map information, but may also require a more detailed analysis which takes into account the patterns of linkage disequilibrium between markers. In this manuscript we show how these two problems are related and present methods which address both problems.
In order to proceed with identifying a large subset of unrelated individuals it is important to quantify relatedness between individuals; one way to do this is to use genotype information to assess Identity by State (IBS) which in turn can be used to estimate pairwise Identity by Descent (IBD), for example by using a Hidden Markov Model as is done in PLINK [2]. This is not the only possibility, in [3] genotype information is used to find an estimator for the coefficient of kinship between individuals which is used as a measure of relatedness. Other possibilities for relating allele sharing to similarities between individuals are discussed in [4] and [5]. Once the pairwise IBD (or some other similarity measure) is known for all pairs of individuals in the population, we can assemble a symmetric matrix of distances between individuals. As has been pointed out in [6], once a suitable threshold for defining unrelatedness has been specified it becomes possible to define an undirected graph whose nodes correspond to individuals with edges connecting unrelated individuals. This is not the only possibility, it is also possible to define a graph where edges are present between individuals who are related. Both definitions lead to graphs whose structure can be analyzed to extract subsets of mutually unrelated individuals. FastIndep makes use of a graph in which related individuals are connected by edges. Regardless of how the graph is specified a threshold must be defined, and the value of this threshold will depend on the similarity measure used. For example in [3], it is shown how the estimate for the coefficient of kinship along with the information on the number of markers can be used to fix this threshold and thus the edge structure of the graph. The graph, in addition to being undirected is also unweighted, i.e. the actual extent of similarity between the individuals is not as important as whether or not the similarity falls above or below a certain user defined threshold. For addressing the problem of finding a large number of unlinked markers, the undirected graph arises from assigning a node to each marker and an edge between two nodes if the linkage disequilibrium between the corresponding markers is above some threshold for statistical significance as defined by the user.
Finding a large number of mutually unrelated individuals or mutually unrelated markers corresponds to finding a large subset of the nodes in the graph such that no two nodes in the subset are connected by one of the edges in the graph. Such subsets of nodes constitute independent sets, and for our purposes we seek maximal independent sets i.e. independent sets with the property that there exists no other node in the graph which can be added to the independent set while retaining the property that all nodes remain mutually unrelated. In graphs with sizes corresponding to real data there may be many maximal independent sets and the question which arises is, which one of these to choose. Frequently, we wish to maximize the number of markers or unrelated individuals. This requires finding the largest maximal independent set present in the graph. The maximal independent set of largest size is the maximum independent set, and in many graphs there may be more than one maximum independent set. Finding the maximum independent set is NPhard, i.e. there is no known efficient algorithm for finding the maximum independent set for an arbitrary graph. That is why exact algorithms, such as the Bron Kerbosch algorithm [7] which can be used to find all maximal independent sets (including maximum independent sets) on an arbitrary graph, are prohibitively slow for graphs beyond a certain size. This necessitates the use of heuristics for problems with realistic sizes, for reasons which will be clear later FastIndep uses a stochastic heuristic. The need for heuristics has been recognized in earlier work [6] where the use of the Bron Kerbosch algorithm is restricted to small graphs with upto a few hundred nodes at most, and for larger graphs a deterministic heuristic algorithm which outputs a single maximal independent set is used.
The FastIndep algorithm was first introduced for selecting unlinked markers for analyzing population stratification and was first discussed in [8]. FastIndep differs from the Primus software of [6] in a number of respects most notably in that the FastIndep algorithm is a stochastic greedy heuristic regardless of the size of the problem. The use of a stochastic heuristic is motivated by [9] where it is demonstrated that randomization may improve on deterministic heuristics for graph triangulations arising in genetic linkage analysis. Furthermore as pointed out in [6], when analyzing a population containing both healthy individuals and those affected with some disease it may be useful to consider maximal independent sets containing a large number of affected individuals; such sets may or may not be maximum independent sets. If the graph is sufficiently large, then it may not be possible to use the Bron Kerbosch algorithm to enumerate all maximal independent sets and pick the most suitable one. In this situation, it may be useful to work with a stochastic heuristic which outputs a number of different maximal independent sets, one or more of which can be choosen. A deterministic heuristic by contrast may not offer the possibility of easily checking for alternatives.
Another advantage to using a stochastic heuristic arises when selecting a subset of markers for analyzing population stratification. If the original set of markers is sufficiently large it is not feasible to use the exact Bron Kerbosch algorithm to select the largest subset of unrelated markers, necessitating the use of some approximate algorithm. This raises the question of the dependence of the final results on the choice of markers. One way to check the extent to which the final results are dependent on the choice of markers is to repeat the analysis of population stratification using alternative sets of markers generated by running the algorithm repeatedly. If the results from using different sets of markers are consistent then the variation of the final results with the choice of markers should be quite small. This consistency check is unique to FastIndep and is not available in other publicly available codes for finding maximal independent sets. The FastIndep code is sufficiently general that it is not restricted to a particular choice of correlation measure, all that matters is that the larger the entry in the matrix, the stronger the correlation between the corresponding individuals (or markers). Any correlation measure that satisfies this criterion may be used so long as the entries of the matrix are larger than or equal to zero. With different correlation measures different thresholds for independence may be required; it is the responsibility of the user to define these thresholds depending on the correlation measure used. For example analysing a marker panel using the Linkage Disequilibrium measure r^{2} the threshold could be choosen based on the relation between the r^{2} and the Pearson Correlation coefficient and thus the χ^{2} distribution. For finding unrelated individuals using the method of [3] the threshold would be fixed based on the sampling distributions for the coancestry described in [3].
Methods
The current C++ implementation of FastIndep is closely related to the one described in [8] for the purpose of selecting a large subset of unlinked markers for analyzing population stratification and will be briefly described below. The key idea behind the algorithm is to start with a deterministic greedy heuristic (which typically finds a reasonably large maximal independent set), and then randomize the heuristic in order to explore solutions close to the one found by the deterministic algorithm. The algorithm is run many times, and due to the stochastic nature of the algorithm the output from the different runs will contain a variety of different maximal independent sets. Some of these maximal independent sets may be larger than the one uncovered by the deterministic algorithm and may not be easy to discover without the randomization, which is the motivation behind introducing a stochastic element to the algorithm. Furthermore, if the algorithm is efficiently implemented, the additional time needed for the multiple runs is not excessive. We emphasize that the algorithm does not make use of map information, this lack of dependence on map information allows the application of the algorithm to genetic data with very complex patterns of Linkage Disequilbrium, as well as to problems where there is no analogue of map information, such as identifying large numbers of unrelated individuals in populations. A description of the algorithm follows, although the same algorithm can be used for identifying both unlinked markers and unrelated invididuals, the description will be in terms of finding unlinked markers. The rest of the section can safely be skipped without any loss of continuity by readers not interested in technical details of the FastIndep algorithm.
As before we assume we have an undirected graph in which is the set of nodes, ={V_{ i }} for i ≤ i ≤ n, and is the set of edges. is not assumed related to any genetic marker map so the algorithm is at this stage perfectly general. In the most general case, the algorithm also requires as input a positive parameter γ which is a measure of how far the algorithm deviates from a deterministic greedy heuristic. The larger the value of γ the closer the algorithm resembles a deterministic heuristic. In the implementation presented here γ is set equal to one. We define D_{ i } the set of non neighbors of V_{ i } and n_{ i } size of D_{ i }, and consider only nodes for which n_{ i } > 0, i.e. we ignore in the discussion nodes connected to all other nodes. This corresponds to ignoring markers in tight Linkage Disequilibrium with all other markers in the panel. We define sets of nodes CandSet, TempSet as well as ReturnSet, which is the output from the program. ReturnSet is initialized by the set of all nodes not connected to any other nodes. The algorithm is initialized by selecting a node disconnected to many others, with a probability which depends on how few neighbors it has.

Initialization
 1.
Evaluate normalized ${p}_{i}=({n}_{i}^{\gamma}/\mathit{\text{norm}})\phantom{\rule{2.77626pt}{0ex}}\phantom{\rule{2.77626pt}{0ex}}\forall \phantom{\rule{2.77626pt}{0ex}}i\phantom{\rule{2.77626pt}{0ex}}$ where $\mathit{\text{norm}}=\sum _{i=1}\phantom{\rule{2.77626pt}{0ex}}{n}_{i}^{\gamma}$
 2.
Pick some V _{ j } with probability p _{ j } and insert in ReturnSet.
 3.
CandSet ← D _{ j }.
 4.
TempSet ← ∅.
 1.

Main Loop
while CandSet ≠ ∅ do 1.
Select some V _{ k } ∈ CandSet with probability p _{ k } normalized over CandSet to insert in ReturnSet
 2.
TempSet ← {V _{ m } ∈ D _{ k } : V _{ m } ∉ ReturnSet where 1 ≤ m ≤ N}
 3.
CandSet←(CandSet ∩ TempSet)
 1.
If CandSet ≡ ∅ return ReturnSet.
In step 1 of the initialization, the normalized p_{ i } are computed for all nodes except nodes which are connected to all others or disconnected to all others. Nodes disconected to all others are automatically included in the final output, while nodes connected to all others do not contribute to maximal independent sets of any interest for our purposes. This leads to a substantial speed up of the program with little loss of information.
It is also possible to define a deterministic greedy heuristic in which the nodes selected in step 2 of the initialization and step 1 of the main loop are just those with the largest associated probability.
Results and discussions
The input to FastIndep consists of a carefully formatted symmetric matrix of positive correlation values (which may exceed one in value), a threshold value taking values between 0 and 1, and an integer specifying the number of times the algorithm is to be run. The input has a structure similar to a dataframe in the R programing language [10], where the column names are the names of the individuals or markers and the row names are the same as the column names. An undirected graph is constructed by assigning a node to each individual and then declaring that edges exist between nodes provided that the correlation between the individuals is greater than the user specified threshold value. As the algorithm is stochastic in nature running the algorithm repeatedly may occasionally yield the same maximal independent set more than once, thus the output of the algorithm is written to a text file which lists only the unique independent sets found. Included in the output of FastIndep is the result from the deterministic greedy heuristic discussed earlier. In all results presented, the random numbers used by FastIndep are generated using the methods of [11].
Results on simulated data
Quantile distribution of set sizes for artificial data set
Graph size  Threshold  Minimum  25%  Median  75%  Maximum 

500  0.15  3  3  4  4  6 
500  0.95  53  63  65  67  77 
1000  0.15  3  4  4  5  7 
1000  0.95  66  75  78  80  91 
2000  0.15  3  4  5  5  7 
2000  0.95  78  89  91  93  104 
In order to test the capacity of FastIndep to find multiple maximum independent sets in graphs of smaller sizes we generated a total of 3111 graphs of sizes ranging from 10 to 60 using a $\mathcal{U}(0,1)$ random number generator to generate symmetric random matrices which were then transformed into graphs using a threshold to define edges, 61 thresholds were used for each graph size. The number of runs in FastIndep was choosen to be 1000 times the number of nodes for each graph. In all cases the largest maximal independent sets found by FastIndep are the same size as those found using a branch and bound algorithm, [12, 13], and in 2597 cases FastIndep found multiple largest cardinality sets. This indicates the capacity of FastIndep to discover multiple maximal independent sets, a feature which is valuable when considering real data.
We next use simulated data to compare FastIndep with three other algorithms used in statistical genetics Primus, KING [14] and PLINK [2]. Both KING and PLINK are implemented in the Primus software, and it is these implmentations we will use. To make the comparison with KING we note that the algorithm in KING is simply the deterministic greedy heuristic which is the starting point for FastIndep. Since the output of FastIndep includes the result from the deterministic greedy heuristic as well as all results from the stochastic heuristic, some of which which may be larger than the result of the greedy heuristic, FastIndep either performs as well as KING or outperforms KING on any input dataset.
Comparing FastIndep with Primus
0.05  0.15  0.25  0.35  0.45  0.55  0.65  0.75  0.85  0.95  

900  4(4)  7(7)  8(8)  10(9)  12(11)  15(14)  19(17)  25(24)  39(36)  94(96) 
1000  4(4)  6(6)  8(8)  10(9)  12(12)  15(14)  19(17)  26(23)  39(34)  93(87) 
1100  5(4)  7(6)  8(8)  10(9)  13(11)  15(14)  20(16)  27(25)  39(37)  93(94) 
1200  5(4)  6(6)  8(8)  10(9)  12(11)  15 (14)  20(16)  26(23)  40(39)  95(98) 
1300  5(4)  7(6)  8(8)  10(10)  12(11)  16(14)  21(18)  26(24)  41(38)  97(96) 
Distribution of differences between Primus and FastIndep
Minimum  25%  Median  75%  Maximum 

3  0  1  2  6 
Recently, a new algorithm LEAF [15] was introduced to find maximal independent sets on graphs which arise from protein sequence comparison; as pointed out in [15] such graphs are typically sparse. In the context of finding unrelated individuals in a study population whose members are choosen to be distantly related or unrelated, the sparsity assumption may be valid; however in the context of finding subsets of unrelated markers a sparse graph would correspond to relatively little Linkage Disequilibrium between markers, which is a needlessly restrictive assumption. More particularly, the capacity of of a handfull of haplotype tagging SNPs in genetic marker maps to tag a large fraction of the markers present [16] suggests that the corresponding undirected graphs may have a very different structure from the sparse graphs arising from the analysis of protein sequences. As mentioned in [15] the authors use several BHOSLIB benchmarks to compare LEAF with other methods; we have used the four largest datsets which have 1272 nodes (frb5324), 1400 nodes (frb5625), 1534 nodes (frb5926) and 4000 nodes (frb10040). With 1000 iterations, FastIndep outperforms LEAF on all of these benchmarks based on the results shown in Figure six of [15]. On the the largest dataset used running FastIndep 1000 times yields 850 maximal independent sets larger than the one found by LEAF, and on this data set outperforms three of the five methods discussed in [15]. However none of the methods discussed in [15] permit the generation of multiple maximal independent sets, which as we have argued earlier is a desirable feature when analyzing population stratification. These considerations suggest that FastIndep has some distinct advantages for creating multiple maximal independent subsets from panels of thousands of markers for crosschecking the analysis of population stratification.
Comparison with true maximum independentset size
Connectivity  Largest set found  True size of maximum 

by FastIndep  independent set  
0.0976  33  37 
0.203  22  23 
0.304  15  17 
0.398  12  14 
0.503  10  10 
0.596  9  9 
0.700  7  7 
0.798  6  6 
0.899  4  4 
Results on real data
FastIndep has been used for finding unlinked markers in data sets with very complex Linkage Disequilibrium patterns such as [17] where the extent of Linkage Disequilbrium in the ancestral populations varies so much that there is no easy way to use map information alone to select a set of independent markers. Retaining all the markers in [17] for analyzing population stratification along the lines of STRUCTURE [18] leads to poor convergence, however the use of a maximal independent set of markers selected along the lines of [8] leads to much better results because the algorithm in [8] ignores map information in selecting markers. As an added feature, the randomized nature of the algorithm permits the generation of multiple different sets of markers which can be successively used in the analysis of population stratification to ensure that the final results are independent of the choice of the set of markers; this procedure has been followed in [17].
Quantiles of data set 1
Minimum  25%  Median  75%  Maximum 

0.0147  0.2353  0.2794  0.3235  0.7353 
Maximal independent set size distribution
Connectivity  Minimum  25%  Median  75%  Maximum 

0.899  3  5  5  6  8 
0.346  40  48  50  52  76 
0.0328  611  654  664  672  846 
0.00142  2697  2740  2751  2762  2907 
0.0000679  3854  3870  3874  3878  3895 
Maximal independent set size distribution
Minimum  25%  Median  75%  Maximum 

184  198  202  206  234 
It is interesting to note that the connectivity of the graph arising from the LD matrix is in the regime where the discrepancy between the size of the largest maximal independent set found by FastIndep and exact answer may be small judging by the results in Table 4. This is consistent with the results in [17] where the LD matrix of a small subset of markers was used to check the discrepancy between the largest independent set found by FastIndep and the exact answer; the discrepancy was found to be small.
Conclusions
We have introduced FastIndep, a stochastic heuristic for finding maximal independent sets on large undirected graphs. The results presented in this paper indicate that FastIndep provides an efficient and competitive heuristic for generating maximal independent sets even for graphs with thousands of nodes, sizes which can easily arise in modern genetic analyses. This is a valuable feature when selecting unrelated individuals or unrelated markers for analyzing population stratification. Indeed, modern biobanks may contain many tens of thousands of samples, so efficiently identifying unrelated samples in such datasets is crucial. Furthermore the capability of FastIndep to find multiple solutions with the same cardinality may shed light on the existence of multiple maximal independent sets which deterministic heuristics may fail to uncover. We have also argued that this capability of finding multiple maximal independent subsets provides a means for cross checking results obtained in the analysis of population stratification. Since the algorithm in FastIndep is heuristic, there is no guarentee that the set of largest cardinality found is the maximum independent set, or that the algorithm has discovered all distinct maximum independent sets. If there is some doubt, the program can be allowed to run more times in order to better explore the graph structure. Comparison with exact algorithms on smaller data sets suggest the discrepancy between the size of the largest maximal independent set found by FastIndep and size of the maximum independent set is largest when connectivities are small, i.e. for sparse graphs.
Declarations
Acknowledgement
This work was funded by the program Professor Visitante Exterior of CAPES Brazil. Many of the simulations were done on computing equipment acquired through Fapesp Project No. 09/531616. We thank the Avileña NegraIberica Breed Association for providing the genotypes of the cattle population studied in this paper. CD wishes to acknowledge support from the grant S2009AGR1704GENR (NEWGAN) from Consejeria de Educacion de la Comunidad de Madrid.
Authors’ Affiliations
References
 Alexander J, Novembre J, Lange K: Fast modelbased estimation of ancestry in unrelated individuals. Genome Res. 2009, 19: 16551664. 10.1101/gr.094052.109.PubMed CentralView ArticlePubMedGoogle Scholar
 Purcell S, Neale B, ToddBrown K, Thomas L, Ferreira M, Bender D, Maller J, Sklar P, de Bakker P, Daly M, Sham P: PLINK: a tool set for wholegenome association and populationbased linkage analyses. Am J Hum Genet. 2007, 81: 559575. 10.1086/519795.PubMed CentralView ArticlePubMedGoogle Scholar
 Manichaikul A, Mychaleckyj J, Rich S, Daly K, Michèle S, WeiMin C: Robust relationship inference in genomewide association studies. Bioinformatics. 2010, 26: 28672873. 10.1093/bioinformatics/btq559.PubMed CentralView ArticlePubMedGoogle Scholar
 Bowcock A, Ruizilinares A, Tomfohrde J, Minch E, Kidd J, CavalliSforza L: Highresolution of human evolutionary trees with polymorphic microsatellites. Nature. 1994, 368: 455457. 10.1038/368455a0.View ArticlePubMedGoogle Scholar
 Stevens E, Heckenberg G, Roberson E, Baugher J, Downey T, Pevsner J: Inference of relationships in population data using identitybydescent and identitybystate. PLoS Genet. 2011, 7 (9): e100228710.1371/journal.pgen.1002287. doi:10.1371/journal.pgen.1002287PubMed CentralView ArticlePubMedGoogle Scholar
 Staples J, Nickerson D, Below J: Utilizing graph theory to select the largest set of unrelated individuals for genetic analysis. Genet Epidemiol. 2012, 99: 3354.Google Scholar
 Bron C, Kerbosch J: Algorithm 457: finding all cliques of an undirected graph. Commun Assoc Comput Mach. 1973, 16: 575577.Google Scholar
 Abraham K, Fernando R: Applications of graphical clustering algorithms in genome wide association mapping. New Frontiers in Graph Theory. Edited by: Zhang Y. 2012, Intech Open RijekaGoogle Scholar
 Fishelson M, Geiger D: Optimizing exact genetic linkage computations. J Comput Biol. 2004, 11: 263275. 10.1089/1066527041410409.View ArticlePubMedGoogle Scholar
 R Core Team: R: A Language and Environment for Statistical Computing. 2013, Vienna, Austria: R Foundation for Statistical Computing, [http://www.Rproject.org/]. [ISBN 3900051070]Google Scholar
 Matsumoto M, Nishimura T: Mersenne Twister: A 623dimensionally equidistributed uniform pseudorandom number generator. ACM Trans Model Comput Simul. 1998, 8: 330. 10.1145/272991.272995.View ArticleGoogle Scholar
 Ostergard P: A fast algorithm for the maximum clique problem. Discrete Appl Math. 2002, 120: 197207. 10.1016/S0166218X(01)002906.View ArticleGoogle Scholar
 Niskanes S, Ostergard P: Cliquer User’s Guide Version 1.0. 2003, Espoo, Finland: Helsinki University of TechnologyGoogle Scholar
 Manichaikul A, Palmas W, Rodriguez C, Peralta C, Divers J, Guo X, Chen WM, Wong Q, Williams K, Kerr K, Taylor K, Tsai M, Goodarzi M, Sale M, DiezRoux A, Rich S, Rotter J, Mychaleckyj J: Population structure of hispanics in the United States: the multiethnic study of atherosclerosis. PLoS Genet. 2012, 8 (4): e100264010.1371/journal.pgen.1002640. doi:10.1371/journal.pgen.1002640PubMed CentralView ArticlePubMedGoogle Scholar
 Bull S, Muldoon M, Doig A: Maximising the size of nonredundant protein datasets using graph theory. PLoS ONE. 2013, 8 (2): e5548410.1371/journal.pone.0055484. doi:10.1371/journal.pone.0055484PubMed CentralView ArticlePubMedGoogle Scholar
 Halldorson B, Bafna V, Lippert R, Schwartz R, de la Vega F, Clark A, Istrail S: Optimal haplotype blockfree selection of tagging SNPs for genomewide association studies. Genome Res. 2004, 14: 16331640. 10.1101/gr.2570004.View ArticleGoogle Scholar
 Hamblin M, Close T, Bhat P, Chao S, Abraham K, Blake T, Brooks W, Cooper D, Griffey C, Hayes P, Hole D, Horsley R, Obert D, Smith K, Ullrich S, Muelhlbauer G, Jannink JL: Population structure & linkage disequilibrium in US Barley germplasm, implications for association mapping. Crop Sci. 2010, 50: 556566. 10.2135/cropsci2009.04.0198.View ArticleGoogle Scholar
 Falush D, Pritchard J, Rosenberg N: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945959.Google Scholar
 PerezFigueroa A, Saura M, Fernandez J, Toro M, Caballero A: METAPOP A software for the management of and analysis of subdivided populations in conservation programs. Conserv Genet. 2009, 10: 10971099. 10.1007/s1059200897187.View ArticleGoogle Scholar
 Hill W, Weir B: Variation in actual relationship as a consequence of Mendelian sampling and linkage. Genet Res Camb. 2011, 93: 4764. 10.1017/S0016672310000480.PubMed CentralView ArticlePubMedGoogle Scholar
 Gentleman RC, Carey VJ, Bates DM: Bioconductor: Open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R8010.1186/gb2004510r80. [http://genomebiology.com/2004/5/10/R80]PubMed CentralView ArticlePubMedGoogle Scholar
 Clayton D: snpStats: SnpMatrix and XSnpMatrix classes and methods. 2012, [R package version 1.8.2]Google Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.