Hierarchical Analysis of Nucleotide Diversity in Geographically Structured Populations

Existing methods for analyzing nucleotide diversity require investigators to identify relevant hierarchical levels before beginning the analysis. We describe a method that partitions diversity into hierarchical components while allowing any structure present in the data to emerge naturally. We present an unbiased version of NEI's nucleotide diversity statistics and show that our modification has the same properties as WRIGHT's F(ST). We compare its statistical properties with several other F(ST) estimators, and we describe how to use these statistics to produce a rooted tree of relationships among the sampled populations in which the mean time to coalescence of haplotypes drawn from populations belonging to the same node is smaller than the mean time to coalescence of haplotypes drawn from populations belonging to different nodes. We illustrate the method by applying it to data from a recent survey of restriction site variation in the chloroplast genome of Coreopsis grandiflora.     

Genetic diversity of ten Egyptian chicken strains using 29 microsatellite markers

In this study, we assessed the genetic diversity of three Egyptian local chicken strains (Fayoumi, Dandarawi and Sinai) and six synthetic breeds derived from Fayoumi and Sinai by intercrossing with Barren Plymouth Rock, Rhode Island Red or White Cornish. Diversity measures were based on interrogation of 29 microsatellites. We identified three main clusters of chicken populations encompassing selected Fayoumi lines and Doki-4 (cluster-1), native Dandarawi (cluster-2) and Sinai, and all six synthetic breeds (cluster-3). Dandarawi and Fayoumi lines exhibited lower intra-population genetic diversity and allelic privacy than Sinai and synthetic breeds. The global inbreeding (F(IT) ) was 0.11, among-population differentiation (F(ST) ) was 0.07, and within-population differentiation (F(IS) ) was 0.04. The between-population marker-estimated kinship was lower than within-population estimates. The cluster analysis classified the Fayoumi lines, Dandarawi and Gimmizah as clearly separated populations. The other strains were configured in mosaic admixed groups.     

Effects of metapopulation processes on measures of genetic diversity

Many species persist as a metapopulation under a balance between the local extinction of subpopulations or demes and their recolonization through dispersal from occupied patches. Here we review the growing body of literature dealing with the genetic consequences of such population turnover. We focus our attention principally on theoretical studies of a classical metapopulation with a 'finite-island' model of population structure, rather than on 'continent-island' models or 'source-sink' models. In particular, we concern ourselves with the subset of geographically subdivided population models in which it is assumed that all demes are liable to extinction from time to time and that all demes receive immigrants. Early studies of the genetic effects of population turnover focused on population differentiation, such as measured by F(ST). A key advantage of F(ST) over absolute measures of diversity is its relative independence of the mutation process, so that different genes in the same species may be compared. Another advantage is that F(ST) will usually equilibrate more quickly following perturbations than will absolute levels of diversity. However, because F(ST) is a ratio of between-population differentiation to total diversity, the genetic effects of metapopulation processes may be difficult to interpret in terms of F(ST) on its own, so that the analysis of absolute measures of diversity in addition is likely to be informative. While population turnover may either increase or decrease F(ST), depending on the mode of colonization, recurrent extinction and recolonization is expected always to reduce levels of both within-population and species-wide diversity (piS and piT, respectively). One corollary of this is that piS cannot be used as an unbiased estimate of the scaled mutation rate, theta, as it can, with some assumptions about the migration process, in species whose demes do not fluctuate in size. The reduction of piT in response to population turnover reflects shortened mean coalescent times, although the distribution of coalescence times under extinction colonization equilibrium is not yet known. Finally, we review current understanding of the effect of metapopulation dynamics on the effective population size.     

Estimating the unbiased estimator theta for population genetic survey data

We consider a method of approximating Weir and Cockerham's theta, an unbiased estimator of genetic population structure, using values readily available from published studies using biased estimators (Wright's F(ST) or Nei's G(ST)). The estimation algorithm is shown to be useful for both model populations and real-world avian populations. However, the correlation between Wright's F(ST) and Weir and Cockerham's theta is strong when compared among 39 empirical avian datasets. Thus, the advantage of approximating an unbiased estimator is unclear considering the small actual effect of theta's bias-removing power on empirical datasets.     

Empirical Bayes inference of pairwise F(ST) and its distribution in the genome

Populations often have very complex hierarchical structure. Therefore, it is crucial in genetic monitoring and conservation biology to have a reliable estimate of the pattern of population subdivision. F(ST)'s for pairs of sampled localities or subpopulations are crucial statistics for the exploratory analysis of population structures, such as cluster analysis and multidimensional scaling. However, the estimation of F(ST) is not precise enough to reliably estimate the population structure and the extent of heterogeneity. This article proposes an empirical Bayes procedure to estimate locus-specific pairwise F(ST)'s. The posterior mean of the pairwise F(ST) can be interpreted as a shrinkage estimator, which reduces the variance of conventional estimators largely at the expense of a small bias. The global F(ST) of a population generally varies among loci in the genome. Our maximum-likelihood estimates of global F(ST)'s can be used as sufficient statistics to estimate the distribution of F(ST) in the genome. We demonstrate the efficacy and robustness of our model by simulation and by an analysis of the microsatellite allele frequencies of the Pacific herring. The heterogeneity of the global F(ST) in the genome is discussed on the basis of the estimated distribution of the global F(ST) for the herring and examples of human single nucleotide polymorphisms (SNPs).     

Assessment of genetic diversity of Korean native pig (Sus scrofa) using AFLP markers

In order to assess the genetic diversity and genetic relationships among the six commercial pig breeds including the Korean native pig, we performed an amplified fragment length polymorphism (AFLP) analysis. Applying the three EcoRI/TagI primer combinations to 54 individual pig samples out of six breeds, a total of 186 AFLP bands were generated, 67 (36%) of which were identified as polymorphic bands. From these polymorphic bands, the three estimates (percentage of polymorphic loci, Neis heterozygosity and Shannon index) of genetic diversity, G(ST) estimates, Neis unbiased genetic distance and two indices of genetic similarity were calculated. From all the calculations of genetic diversity, the lowest genetic diversity was exhibited in the Korean native pig, and the highest in the Chinese Yanbian pig. Given the mean G(ST) value (G(ST) = 0.390) across all pigs examined, levels of apparent breed subdivision were considerable. A UPGMA tree of individuals based on Jaccards similarity index showed that the Korean native pig formed a distinct cluster from the other five pigs. In addition, the tree displayed that all the individuals except for six individuals were grouped into their breeds. Principal component analysis based on the binary data matrix of either presence or absence confirmed the distinctness of the Korean native pig from the other pigs. Our results indicate that the Korean native pig has a low level of genetic diversity and is distinct from the five pig breeds, confirming the results from previous microsatellite data. The findings also suggest that AFLP analysis may be a valuable tool for revealing genetic relationships and genetic diversity among different pig breeds.     

Genetic structure in the coral-reef-associated Banggai cardinalfish, Pterapogon kauderni

In this study, we used 11 polymorphic microsatellite loci to show that oceanic distances as small as 2-5 km are sufficient to produce high levels of population genetic structure (multilocus F(ST) as high as 0.22) in the Banggai cardinalfish (Pterapogon kauderni), a heavily exploited reef fish lacking a pelagic larval dispersal phase. Global F(ST) among all populations, separated by a maximum distance of 203 km, was 0.18 (R(ST) = 0.35). Moreover, two lines of evidence suggest that estimates of F(ST) may actually underestimate the true level of genetic structure. First, within-locus F(ST) values were consistently close to the theoretical maximum set by the average within-population heterozygosity. Second, the allele size permutation test showed that R(ST) values were significantly larger than F(ST) values, indicating that populations have been isolated long enough for mutation to have played a role in generating allelic variation among populations. The high level of microspatial structure observed in this marine fish indicates that life history traits such as lack of pelagic larval phase and a good homing ability do indeed play a role in shaping population genetic structure in the marine realm.     

Species diversity and population density affect genetic structure and gene dispersal in a subtropical understory shrub

Aims The dispersal of pollen and seeds is spatially restricted and may vary among plant populations because of varying biotic interactions, population histories or abiotic conditions. Because gene dispersal is spatially restricted, it will eventually result in the development of spatial genetic structure (SGS), which in turn can allow insights into gene dispersal processes. Here, we assessed the effect of habitat characteristics like population density and community structure on small-scale SGS and estimate historical gene dispersal at different spatial scales.Methods In a set of 12 populations of the subtropical understory shrub Ardisia crenata, we assessed genetic variation at 7 microsatellite loci within and among populations. We investigated small-scale genetic structure with spatial genetic autocorrelation statistics and heterogeneity tests and estimated gene dispersal distances based on population differentiation and on within-population SGS. SGS was related to habitat characteristics by multiple regression.Important findings The populations showed high genetic diversity (H e = 0.64) within populations and rather strong genetic differentiation (F ′ ST = 0.208) among populations, following an isolation-by-distance pattern, which suggests that populations are in gene flow–drift equilibrium. Significant SGS was present within populations (mean Sp = 0.027). Population density and species diversity had a joint effect on SGS with low population density and high species diversity leading to stronger small-scale SGS. Estimates of historical gene dispersal from between-population differentiation and from within-population SGS resulted in similar values between 4.8 and 22.9 m. The results indicate that local-ranged pollen dispersal and inefficient long-distance seed dispersal, both affected by population density and species diversity, contributed to the genetic population structure of the species. We suggest that SGS in shrubs is more similar to that of herbs than to trees and that in communities with high species diversity gene flow is more restricted than at low species diversity. This may represent a process that retards the development of a positive species diversity–genetic diversity relationship.     

Complex pattern of genetic structuring in the Atlantic salmon (Salmo salar L.) of the River Foyle system in northwest Ireland: disentangling the evolutionary signal from population stochasticity

Little is known about the microevolutionary processes shaping within river population genetic structure of aquatic organisms characterized by high levels of homing and spawning site fidelity. Using a microsatellite panel, we observed complex and highly significant levels of intrariver population genetic substructure and Isolation-by-Distance, in the Atlantic salmon stock of a large river system. Two evolutionary models have been considered explaining mechanisms promoting genetic substructuring in Atlantic salmon, the member-vagrant and metapopulation models. We show that both models can be simultaneously used to explain patterns and levels of population structuring within the Foyle system. We show that anthropogenic factors have had a large influence on contemporary population structure observed. In an analytical development, we found that the frequently used estimator of genetic differentiation, F(ST), routinely underestimated genetic differentiation by a factor three to four compared to the equivalent statistic Jost's D(est) (Jost 2008). These statistics also showed a near-perfect correlation. Despite ongoing discussions regarding the usefulness of "adjusted"F(ST) statistics, we argue that these could be useful to identify and quantify qualitative differences between populations, which are important from management and conservation perspectives as an indicator of existence of biologically significant variation among tributary populations or a warning of critical environmental damage.     

A study of Y-chromosome microsatellite variation in sub-Saharan Africa: a comparison between F(ST) and R(ST) genetic distances

Seven Y-chromosome microsatellite loci (DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, and DYS393) were analyzed in three populations from sub-Saharan Africa: the Bamileke and Ewondo populations from Cameroon and the Hutu from Rwanda. Complete typing was obtained for 112 individuals, and a total of 53 different haplotypes was observed. The single-locus gene diversity, averaged across populations, ranges from 0.100 for the DYS392 locus to 0.610 for the DYS389I locus. The haplotype diversity ranges from 0.832 (Ewondo) to 0.965 (Hutu), with an intermediate value of 0.918 in the Bamileke. The diversity among Bamileke, Ewondo, Hutu, and other sub-Saharan populations selected from the literature was analyzed using both a classical (F(ST)) and a stepwise-based (R(ST)) genetic distance method. The pattern of interpopulational diversity based on F(ST) was congruent with anthropological knowledge, while that based on R(ST) revealed unexpected and unconvincing population affinities. From a practical point of view, our study indicates that Y-chromosome microsatellite data may provide useful information for analyses of interpopulational diversity among sub-Saharan populations if an adequate number of loci and individuals along with an appropriate genetic distance method are used. On a theoretical ground, we propose that the lesser performance of R(ST) compared to F(ST) could be explained by the important role played by genetic drift in shaping the relationships among examined populations.     

Human genetic diversity and the nonexistence of biological races

Sewall Wright's population structure statistic, F(ST), measured among samples of world populations is often 15% or less. This would indicate that 85% of genetic variation occurs within groups while only 15% can be attributed to allele frequency differences among groups. In this paper, we show that this low value reflects strong biases that result from violating hidden assumptions that define F(ST). These limitations on F(ST) are demonstrated algebraically and in the context of analyzing dinucleotide repeat allele frequencies for a set of eight loci genotyped in eight human groups and in chimpanzees. In our analyses, estimates of F(ST) fail to identify important variation. For example, when the analysis includes only humans, F(ST) = 0.119, but adding the chimpanzees increases it only a little, F(ST) = 0.183. By relaxing the underlying statistical assumptions, the results for chimpanzees become consistent with common knowledge, and we see a richer pattern of human genetic diversity. Some human groups are far more diverged than would be implied by standard computations of F(ST), while other groups are much less diverged. We discuss the relevance of these findings to the application of biological race concepts to humans. Four different race concepts are considered: typological, population, taxonomic, and lineage. Surprisingly, a great deal of genetic variation within groups is consistent with each of these concepts. However, none of the race concepts is compatible with the patterns of variation revealed by our analyses.     

Population genetic structure of Arabidopsis lyrata in Europe

Population genetic theory predicts that the self-incompatible and perennial herb, Arabidopsis lyrata, will have a genetic structure that differs from the self-fertilizing, annual Arabidopsis thaliana. We quantified the genetic structure for eight populations of A. lyrata ssp. petraea in historically nonglaciated regions of central Europe. Analysis of 20 microsatellite loci for 344 individuals demonstrated that, in accordance with predictions, diploid populations had high genome-wide heterozygosity (H(O) = 0.48; H(E) = 0.52), high within-population diversity (83% of total) compatible with mutation-drift equilibrium, and moderate differentiation among populations (F(ST) = 0.17). Within a single population, the vast majority of genetic variability (92%) was found at the smallest spatial scale (< 3 m). Although there was no evidence of biparental inbreeding or clonal propagation at this scale (F(IS) = 0.003), significant fine-scale spatial autocorrelation indicated localized gene flow presumably due to gravity dispersed seeds (Sp = 0.018). Limited gene flow between isolated population clusters (regions) separated by hundreds of kilometres has given rise to an isolation by distance pattern of diversification, with low, but significant, differentiation among regions (F(ST) = 0.05). The maintenance of geographically widespread polymorphisms and uniformly high diversity throughout central Europe is consistent with periglacial survival of A. lyrata ssp. petraea north of the Alps in steppe-tundra habitats during the last glacial maximum. As expected of northern and previously glaciated localities, A. lyrata in Iceland was genetically less diverse and highly differentiated from central Europe (H(E) = 0.37; F(ST) = 0.27).     

Effects of life history traits and sampling strategies on genetic diversity estimates obtained with RAPD markers in plants

A compilation of studies using RAPD markers for evaluating population differentiation resulted in 78 estimates of AMOVA-derived Φ_ST and 31 estimates of Nei's G_ST, as well as in 41 estimates of Nei's within-population diversity. In outcrossing taxa, estimates of between-population diversity were closely correlated with maximum geographic distance between sampled populations. A corresponding association was not found in selfing taxa. These results suggest that RAPD can be a sensitive method for detection of genetic structuring according to the isolation-by-distance model. However, it also means that sampling strategies, as applied in individual studies, can seriously influence the resulting estimates of between-population diversity. Other sampling strategies, like number of plants per population and number of scored polymorphic markers, do not seem to impart any serious artefacts. As previously verified with allozyme data, RAPD markers showed that long-lived, outcrossing, late successional taxa retain most of their genetic variability within populations. By contrast, annual, selfing and/or early successional taxa allocate most of the genetic variability among populations. Estimates for between- and within-population diversity, respectively, proved to be negatively correlated, as previously reported for allozyme data. The only major discrepancy between allozymes and RAPD markers concerns geographic range; within-population diversity was strongly affected by distributional range of the investigated species in the allozyme data but not in the RAPD data. Moreover, RAPD-based values for between-population diversity increased with increasing distributional range whereas the opposite has been reported in a large allozyme data compilation. Contrary to allozymes, RAPD marker-derived within-population diversity is probably therefore not a very good predictor of total species genetic diversity.     

Microsatellite allele sizes: a simple test to assess their significance on genetic differentiation

The mutation process at microsatellite loci typically occurs at high rates and with stepwise changes in allele sizes, features that may introduce bias when using classical measures of population differentiation based on allele identity (e.g., F(ST), Nei's Ds genetic distance). Allele size-based measures of differentiation, assuming a stepwise mutation process [e.g., Slatkin's R(ST), Goldstein et al.'s (deltamu)(2)], may better reflect differentiation at microsatellite loci, but they suffer high sampling variance. The relative efficiency of allele size- vs. allele identity-based statistics depends on the relative contributions of mutations vs. drift to population differentiation. We present a simple test based on a randomization procedure of allele sizes to determine whether stepwise-like mutations contributed to genetic differentiation. This test can be applied to any microsatellite data set designed to assess population differentiation and can be interpreted as testing whether F(ST) = R(ST). Computer simulations show that the test efficiently identifies which of F(ST) or R(ST) estimates has the lowest mean square error. A significant test, implying that R(ST) performs better than F(ST), is obtained when the mutation rate, mu, for a stepwise mutation process is (a) >/= m in an island model (m being the migration rate among populations) or (b) >/= 1/t in the case of isolated populations (t being the number of generations since population divergence). The test also informs on the efficiency of other statistics used in phylogenetical reconstruction [e.g., Ds and (deltamu)(2)], a nonsignificant test meaning that allele identity-based statistics perform better than allele size-based ones. This test can also provide insights into the evolutionary history of populations, revealing, for example, phylogeographic patterns, as illustrated by applying it on three published data sets.     

If FST does not measure neutral genetic differentiation,then comparing it with QST is misleading. Or is it?

The comparison between neutral genetic differentiation (F(ST) ) and quantitative genetic differentiation (Q(ST) ) is commonly used to test for signatures of selection in population divergence. However, there is an ongoing discussion about what F(ST) actually measures, even resulting in some alternative metrics to express neutral genetic differentiation. If there is a problem with F(ST) , this could have repercussions for its comparison with Q(ST) as well. We show that as the mutation rate of the neutral marker increases, F(ST) decreases: a higher within-population heterozygosity (He) yields a lower F(ST) value. However, the same is true for Q(ST) : a higher mutation rate for the underlying QTL also results in a lower Q(ST) estimate. The effect of mutation rate is equivalent in Q(ST) and F(ST) . Hence, the comparison between Q(ST) and F(ST) remains valid, if one uses neutral markers whose mutation rates are not too high compared to those of quantitative traits. Usage of highly variable neutral markers such as hypervariable microsatellites can lead to serious biases and the incorrect inference that divergent selection has acted on populations. Much of the discussion on F(ST) seems to stem from the misunderstanding that it measures the differentiation of populations, whereas it actually measures the fixation of alleles. In their capacity as measures of population differentiation, Hedrick's G'(ST) and Jost's D reach their maximum value of 1 when populations do not share alleles even when there remains variation within populations, which invalidates them for comparisons with Q(ST) .     

A Note on Combining Dependent Tests of Significance

In combining several tests of significance the individual test statistics are allowed to be stochastically dependent. By choosing the weighted inverse normal method for the combination, the dependency of the original test statistics is then characterized by a correlation of the transformed statistics. For this correlation a confidence region, an unbiased estimator and an unbiased estimate of its variance are derived. The combined test statistic is extended to include the case of possibly dependent original test statistics. Simulation studies show the performance of the actual significance level.     

横断山区高山绣线菊的谱系地理学研究

横断山区作为青藏高原东南部主要的一个冰期避难所,第四纪冰期气候的变化对该地区的植物地理分布和居群遗传结构都产生了重要的影响。为了揭示该地区物种分布的分子系统地理学结构,选取在该地区广泛分布的一种高山灌木-高山绣线菊的叶绿体trnL-trnF序列进行研究。采集了15个居群182个个体进行测序,共发现7个单倍型。总的遗传多样性较高(HT=0.809),但居群内遗传多样性较低(HS=0.236)。分子变异分析(AMOVA)结果表明分布区内高山绣线菊的遗传变异主要存在于居群间(84.48%),且居群间的遗传分化很高(GST=0.708,FST=0.84476,NST=0.863),有着显著的谱系地理学结构(NST>GST,P<0.01)和较低的居群间平均基因流(Nm=0.09)。单倍型的系统进化树和进化分支网络分析得到了相似的拓扑结构,7种单倍型都按照地理分布聚为三支:横断山区西部、横断山区东部以及两者的交接地带。本研究推测该物种在横断山区存在多个冰期避难所,而没有表现出大规模的种群集体扩张和迁移的现象。青藏高原隆升、第四纪气候的反复波动以及横断山区特殊的地理环境使得原来连续的居群片段化,并发生范围扩张,从而塑造了高山绣线菊的现代生物地理分布格局。     

Nuclear and chloroplast microsatellites reveal extreme population differentiation and limited gene flow in the Aegean endemic Brassica cretica (Brassicaceae)

Nuclear and chloroplast microsatellite markers were used to study population structure and gene flow among seven Cretan populations of the Aegean endemic plant species Brassica cretica (Brassicaceae). Both nuclear and chloroplast markers revealed exceptionally high levels of population differentiation (overall F(ST)=0.628 and 1.000, respectively) and relatively little within-population diversity (overall H(S)=0.211 and 0.000, respectively). Maximum-likelihood estimates of directional migration rates were low among all pairs of populations (average Nm=0.286). There was no evidence that differences in flower colour between populations had any influence on historical levels of gene flow. In addition, a haplotype network showed that all five chloroplast haplotypes found in the sample were closely related. Together, these results suggest that current patterns of diversification in B. cretica are mainly a result of genetic drift during the last half million years. The main conclusions from the present study are consistent with the prevailing hypothesis that plant diversification in the Aegean region is driven by random rather than adaptive differentiation among isolated populations.     

Genetic structure of an endangered plant, Antirrhinum microphyllum (Scrophulariaceae): allozyme and RAPD analysis

Thirteen allozyme loci and 68 random amplified polymorphic DNA (RAPD) markers were analyzed to assess the genetic diversity and population structure of threatened Antirrhinum microphyllum (Scrophulariaceae), a narrow endemic of central Spain known from only four populations. According to allozyme data, species genetic diversity (p = 46.15%, A = 2.61, and H(e) = 0.218), as well as within-population genetic diversity (p = 44.23%, A = 2.10, and H(e) = 0.204), were high when compared to average estimates for other narrowly distributed plant species. Ninety-four percent of species genetic diversity corresponded to within-population genetic diversity. Nevertheless, significant differences were found among populations in allele frequencies of four of the six polymorphic loci, and three private alleles were detected. Inbreeding coefficients (F(IS)) suggest that populations are structured in genetic neighborhoods. The RAPDs also showed high levels of genetic diversity (p = 89.71% and H(e) = 0.188 at the species level, and p = 67.65% and H(e) = 0.171 at the population level). Nei's genetic distances estimated both from allozymes and RAPDs indicated low differentiation among populations. In spite of this, the low frequencies of certain alleles and the presence of private alleles indicate that efforts should be made to conserve all four remaining populations.     

Genetic diversity, structure, gene flow and evolutionary relationships within the Sorghum bicolor wild�Cweedy�Ccrop complex in a western African region

Gene flow between domesticated plants and their wild relatives is one of the major evolutionary processes acting to shape their structure of genetic diversity. Earlier literature, in the 1970s, reported on the interfertility and the sympatry of wild, weedy and cultivated sorghum belonging to the species Sorghum bicolor in most regions of sub-Saharan Africa. However, only a few recent surveys have addressed the geographical and ecological distribution of sorghum wild relatives and their genetic structure. These features are poorly documented, especially in western Africa, a centre of diversity for this crop. We report here on an exhaustive in situ collection of wild, weedy and cultivated sorghum assembled in Mali and in Guinea. The extent and pattern of genetic diversity were assessed with 15 SSRs within the cultivated pool (455 accessions), the wild pool (91 wild and weedy forms) and between them. F (ST) and R (ST) statistics, distance-based trees, Bayesian clustering methods, as well as isolation by distance models, were used to infer evolutionary relationships within the wild-weedy-crop complex. Firstly, our analyses highlighted a strong racial structure of genetic diversity within cultivated sorghum (F (ST) = 0.40). Secondly, clustering analyses highlighted the introgressed nature of most of the wild and weedy sorghum and grouped them into two eco-geographical groups. Such closeness between wild and crop sorghum could be the result of both sorghum's domestication history and preferential post-domestication crop-to-wild gene flow enhanced by farmers' practices. Finally, isolation by distance analyses showed strong spatial genetic structure within each pool, due to spatially limited dispersal, and suggested consequent gene flow between the wild and the crop pools, also supported by R (ST) analyses. Our findings thus revealed important features for the collection, conservation and biosafety of domesticated and wild sorghum in their centre of diversity.