Population genetic structure of wild-growing ginseng (Planax ginseng C. A. Meyer) assessed using AFLP markers

Genetic variability in ten populations of wild-growing ginseng was assessed using AFLP markers with the application of fragment analysis on a genetic analyzer. The variation indices were high in the populations (P = 55.68%; H(S) = 0.1891) and for the species (P = 99.65%; H(S) = 0.2857). Considerable and statistically significant population differentiation was demonstrated (theta = 0.363; Bayesian approach, "full model"; F(ST) = 0.36, AMOVA). The results of AMOVA and Bayesian analysis indicate that 64.46% of variability is found within the populations. Mantel test showed no correlation between the genetic and geographic distances among the populations (r = -0.174; P = 0.817). Hierarchical AMOVA and analysis of genetic relationships based on Euclidean distances (NJ, PCoA, and MST) identified two divergent population groups of ginseng. Low gene flow between these groups (N(m) = 0.4) suggests their demographic independence. In accordance to the concept of evolutionary significant units (ESU), these population groups, in terms of the strategy and tactics for conservation and management of natural resources, should be treated as management units (MUs). The MS tree topology suggests recolonization of southern Sikhote-Alin by ginseng along two directions, from south and west.     

Polymorphism of RAPD,ISSR and AFLP markers of the <Emphasis Type="Italic">Panax ginseng</Emphasis> C. A. Meyer (Araliaceae) genome

The genus Panax (Araliaceae) is world-famous because many its members have important medicinal properties. Panax ginseng C. A. Meyer is more popular than other species of the genus because remedies prepared from this plant stimulate immunity, help to prevent diseases, and have antistress effects. In addition, the ginseng root extract is traditionally used as a means against aging. At present, this species is found in the wild only in Primorsky krai, Russia, but its populations are extremely exhausted and need to be restored. In this study, effectiveness of molecular DNA markers in detecting genetic variation and differentiation of the ginseng populations was tested. Genetic variation of ginseng, identified using RAPD (P = 4%; H _pop = 0.0130) and ISSR (P = 9.3%; H _pop = 0.0139) markers was low. The AFLP* approach, according to which amplicons are separated in polyacrylamide gel and visualized by means of silver staining, showed somewhat higher variability (P = 21.8%; H _pop = 0.0509), while its effectiveness in population differentiation was as low as that of RAPD and ISSR. The AFLP** technique, which included analysis of the fragments using genetic analyzer, revealed high genetic diversity of ginseng (P = 94.4%; H _pop = 0.3246). All populations examined using the AFLP** markers were statistically significantly differentiated based on the AMOVA results. Our result suggest effectiveness of AFLP** markers for characterization of the genetic structure and genetic relationships of the ginseng populations. These markers are recommended for use in large-scale population genetic studies of this species to develop measures of its conservation.     

Genetic variability of Siberian fir <Emphasis Type="Italic">Abies sibirica</Emphasis> Ledeb. inferred from AFLP markers

Genetic variability of AFLP markers was studied in 20 populations of Siberian fir (Abies sibirica, Pinaceae) and in two populations of Far-Eastern species Manchurian fir A. nephrolepis and Sakhalin fir A. sachalinensis each. Four pairs of selective primers were used. In total, 168 samples from three fir species were genotyped for 117 polymorphic loci. According to the AMOVA results, the variability proportion characterizing the differences among three Abies species was several times higher (F _CT = 0.53) than that acounting for population differences within the species (F _SC = 0.125). Differentiation of the A. sibirica populations based on AFLP markers exceeded 14% (F _ST = 0.141). Significant correlation between the genetic distances calculated from the AFLP data and the geographic distances between populations was found. The results of AFLP variability analysis supported and supplemented the conclusions inferred previously from allozyme and cpSSR data: several genetically similar geographic groups of Siberian fir were identified. These groups differ both in allele frequencies and in the levels of genetic variation.     

Protecting evolutionary significant units for the remnant populations of <Emphasis Type="Italic">Berchemiella wilsonii</Emphasis> var. <Emphasis Type="Italic">pubipetiolata</Emphasis> (Rhamnaceae)

Berchemiella wilsonii var. pubipetiolata (Rhamnaceae) is an endangered tree in eastern China. Habitat destruction has resulted in fragmentation of remnant populations and extinction of local populations. AFLP and cpDNA markers were used to determine the population structure of remnant populations of B. wilsonii var. pubipetiolata. Moderate nuclear genomic diversity was found within each of the four remnant populations (H _S = 0.141–0.172), while the cpDNA haplotype diversity in each population ranged from 0.356 to 0.681. Six haplotypes were identified by a combined cpRFLP and cpSSR analysis in a total of 89 individuals. AMOVA revealed significantly AFLP genetic differentiation within and between regions (Φ_SC = 0.196, Φ_CT = 0.396, respectively), and a high cpDNA haplotype differentiation between regions (Φ_CT = 0.849). The results suggest low gene flow between populations of B. wilsonii var. pubipetiolata. Strong genetic divergence between two regional populations as revealed by both AFLP and cpDNA markers provided convincing evidence that two distinct evolutionary lineages existed, and should be recognized as ‘evolutionary significant units’ (ESUs) for conservation concerns.     

Genetic Diversity and Geographical Differentiation of <Emphasis Type="Italic">Desmodium triflorum</Emphasis> (L.) DC. in South China Revealed by AFLP Markers

High levels of genetic variation enable species to adapt to changing environments and provide plant breeders with the raw materials necessary for artificial selection. In the present study, six AFLP primer pairs were used to assess the genetic diversity of Desmodium triflorum (L.) DC. from 12 populations in South China. A high percentage of polymorphic loci (P = 76.16%) and high total gene diversity (H _T = 0.310) were found, indicating that the genetic diversity of D. triflorum is high at the species level. Genetic diversity was also relatively high at the population level (P = 55.85%, H _e = 0.230). The coefficient of gene differentiation among populations (G _ST) was 0.255, indicating that while most genetic diversity resided within populations, there was also considerable differentiation among populations. AMOVA also indicated 24.29% of the total variation to be partitioned among populations (Φ_ST = 0.243). UPGMA clustering analysis based on genetic distances showed that the 12 populations could be separated into three subgroups: an eastern, a western, and a central-southern subgroup. However, a Mantel test revealed no significant correlation (r = 0.286, p = 0.983) between the geographical distances and genetic distances separating these populations; mountain barriers to gene flow and human disturbance may have confounded these correlations. The present study has provided some fundamental genetic data that will be of use in the exploitation of D. triflorum.     

High genetic diversity detected in the endemic Primula apennina Widmer (Primulaceae) using ISSR fingerprinting

Primula apennina Widmer is endemic to the North Apennines (Italy). ISSR were used to detect the genetic diversity within and among six populations representative of the species distribution range. High levels of genetic diversity were revealed both at population percentage of polymorphic band (PPB = 75.92%, H _S = 0.204, H _pop = 0.319) and at species level (PPB = 96.95%, H _T = 0.242, H _sp = 0.381). Nei gene diversity statistics (15.7%), Shannon diversity index (16.3%) and AMOVA (14%) detected a moderate level of interpopulation diversity. Principal coordinate and Bayesian analyses clustered the populations in three major groups along a geographic gradient. The correlation between genetic and geographic distances was positive (Mantel test, r = 0.232). All together, these analyses revealed a weak but significant spatial genetic structure in P. apennina, with gene flow acting as a homogenizing force that prevents a stronger differentiation of populations. Conservation measures are suggested based on the observed pattern of genetic variability.     

Determination of the population genetic structure of the invasive weed <Emphasis Type="Italic">Ageratina adenophora</Emphasis> using ISSR-PCR markers

A genetic analysis of nuclear DNA was performed by inter-simple sequence repeat-polymerase chain reaction (ISSR-PCR) technique on 32 populations of Ageratina adenophora, an invasive triploid weed in China. Among the 100 ISSR markers detected, 12 showed genetic variation both within and among the populations. Among the 446 amplified bands, 93.5% were found polymorphic. Most individuals (99%) displayed a unique ISSR fingerprint pattern, which yielded a high level of polymorphism (P _o = 93.5%) and genetic diversity (Nei’s H _T = 0.2354). The estimates of population variation, based on ISSR-PCR, were high, as measured by the analysis of molecular variance (AMOVA, F _ST = 0.3140), the Wright’s F-statistics (G _ST = 0.3453), and the Shannon’s information index (H _sp = 0.3716). AMOVA revealed 68.6% genetic variation within the populations and 91.2% within the provinces. The Mantel test showed that genetic distance was significantly correlated with geographic distance. Published in Russian in Fiziologiya Rastenii, 2009, Vol. 56, No. 3, pp. 453–459. This text was submitted by the authors in English.     

Genetic analysis of two Portuguese populations of <Emphasis Type="Italic">Ruditapes decussatus</Emphasis> by RAPD profiling

The clam Ruditapes decussatus is commercially important in the south of Portugal. The random amplified polymorphic DNA (RAPD) technique was applied to assess the genetic diversity and population structure of two Portuguese populations occurring in the Ria Formosa (Faro) and the Ria de Alvor, respectively. Twenty-five individuals of each population were investigated by RAPD profiles. Genetic diversity within populations, measured by the percentage of polymorphic loci (%P), varied between 68.57% (Alvor) and 73.88% (Faro). Shannon’s information index (H) and Nei’s gene diversity (h) were 0.281 and 0.176, respectively, for the Alvor population and 0.356 and 0.234 for the Faro population. Overall, genetic variation within R. decussatus populations was high. The total genetic diversity (H _T) was explained by a low variation between populations (G _ST = 0.145), which is consistent with high gene flow (N _m = 2.9). The analysis of molecular variance (AMOVA) showed that 65% of variability is within populations and 35% between populations (Φ_PT = 0.345; P ≥ 0.001). The value of Nei’s genetic distance was 0.0881, showing a low degree of population genetic distance, despite the different geographic origin. This is the first study on the population genetics of R. decussatus by RAPD technique. The results may be useful for restocking programs and aquaculture.     

Patterns of nuclear and mitochondrial DNA variation in Iberian populations of <Emphasis Type="Italic">Emys orbicularis</Emphasis> (Emydidae): conservation implications

The European pond turtle (Emys orbicularis) is threatened and in decline in several regions of its natural range, due to habitat loss combined with population fragmentation. In this work, we have focused our efforts on studying the genetic diversity and structure of Iberian populations with a fine-scale sampling (254 turtles in 10 populations) and a representation from North Africa and Balearic island populations. Using both nuclear and mitochondrial markers (seven microsatellites, ∼1048 bp nDNA and ∼1500 bp mtDNA) we have carried out phylogenetic and demographic analyses. Our results show low values of genetic diversity at the mitochondrial level although our microsatellite dataset revealed relatively high levels of genetic variability with a latitudinal genetic trend decreasing from southern to northern populations. A moderate degree of genetic differentiation was estimated for Iberian populations (genetic distances, F _ST values and clusters in the Bayesian analysis). The results in this study combining mtDNA and nDNA, provide the most comprehensive population genetic data for E. orbicularis in the Iberian Peninsula. Our results suggest that Iberian populations within the Iberian–Moroccan lineage should be considered as a single subspecies with five management units, and emphasize the importance of habitat management rather than population reinforcement (i.e. captive breeding and reintroduction) in this long-lived species.     

Population Genetic Structure of a Widespread Bat-Pollinated Columnar Cactus

Bats are the main pollinators and seed dispersers of Stenocereus thurberi, a xenogamous columnar cactus of northwestern Mexico and a good model to illustrate spatial dynamics of gene flow in long-lived species. Previous studies in this cactus showed differences among populations in the type and abundance of pollinators, and in the timing of flowering and fruiting. In this study we analyzed genetic variability and population differentiation among populations. We used three primers of ISSR to analyze within and among populations genetic variation from eight widely separated populations of S. thurberi in Sonora, Mexico. Sixty-six out of 99 of the ISSR bands (P = 66.7%) were polymorphic. Total heterozygosity for all populations sampled revealed high genetic diversity (H_sp = 0.207, H_BT = 0.224). The AMOVA showed that most of the genetic variation was within populations (80.5%). At the species level, estimates of population differentiation, θ = 0.175 and θ^B = 0.194, indicated moderate gene flow among populations. The absence of a significant correlation between genetic and geographic distances indicated little isolation by geographic distance. The large genetic variation and diversity found in S. thurberi is consistent with its open reproductive system and the high mobility of bats, a major pollinator. However, small changes in number or kind of pollinators and seed dispersal agents, in the directionality of migratory routes, and/or in the timing of flowering and fruiting among populations, can critically affect gene flow dynamics.     

RAPD variation and population genetic structure of Physalaemus cuvieri (Anura: Leptodactylidae) in Central Brazil

Studies about the organization of the genetic variability and population structure in natural populations are used either to understand microevolutionary processes or the effects of isolation by human-inducted landscape modifications. In this paper, we analyzed patterns of genetic population structure using 126 RAPD loci scored for 214 individuals of Physalaemus cuvieri, sampled from 18 local populations. Around 97% of these loci were polymorphic. The among-population variation component (Φ_ST) obtained by AMOVA was equal to 0.101 and θ ^B obtained using a Bayesian approach for dominant markers was 0.103. Genetic divergence, analyzed by Mantel spatial correlogram, revealed only a short-distance significant correlation between genetic and geographic distances. This is expected if low levels of population differentiation, due to high abundance buffering the effect of stochastic processes, are combined with low spatially restricted gene flow. Although this may be consistent with the current knowledge of species’ biology, the spatial distribution of local populations observed in this study also suggest that, at least in part, recent human occupation and habitat fragmentation may also explain part of the interpopulational component of the genetic variation.     

Genetic diversity of the endangered Chinese endemic herb Primulina tabacum (Gesneriaceae) revealed by amplified fragment length polymorphism (AFLP)

Primulina tabacum Hance, is a critically endangered perennial endemic to limestone area in South China. Genetic variability within and among four extant populations of this species was assessed using AFLP markers. We expected a low genetic diversity level of this narrowly distributed species, but our results revealed that a high level of genetic diversity remains, both at population level (55.5% of markers polymorphic, H _E = 0.220, I _S = 0.321), and at species level (P = 85.6% of markers polymorphic, H _E = 0.339, I _S = 0.495), probably resulting from its refugial history and/or breeding system. High levels of genetic differentiation among populations was apparent based on Nei’s genetic diversity analysis (G _st=0.350). The restricted gene flow between populations is a potential reason for the high genetic differentiation. The population genetic diversity of P. tabacum revealed here has clear implications for conservation and management. To maintain present levels of genetic diversity, in situ conservation of all populations is necessary.     

Comparative analysis of genetic structure in natural populations of two <Emphasis Type="Italic">Arabidopsis</Emphasis> species with different degree of panmixia

Comparative analysis of genetic structure of northern natural populations of two Arabidopsis species with different degrees of panmixia was performed. The variability of 121 RAPD loci in seven populations of model plant A. thaliana possessing high degree of self fertility was studied together with 93 RAPD loci in population of cross-pollinating species A. lyrata ssp. petraea. The population of A. l. petraea demonstrated higher level of genetic variability (P _99% = 62.50%; H _exp = 0.169) than the populations of A. thaliana, which is obviously connected with biological features of reproduction of the species. A significant level of genetic variability (P _99% = 42.27%; H _exp = 0.126) was revealed in populations of A. thaliana, which is not typical for self-pollinating plant species. The high population polymorphism of A. thaliana in the northern part of its range may be connected with adverse environmental conditions. The genetic distances between populations of the species studied (average D _N = 0.494) confirm close relatedness between A. thaliana and A. l. petraea.     

Spatial Genetic Structure of Campanula sabatia, a Threatened Narrow Endemic Species of the Mediterranean Basin

Campanula sabatia is endemic to NW Italy (Liguria) and it was included in the European Red list of endangered species due to the heavy human pressure on its habitat. AFLP markers were used to detect the genetic diversity within and among ten populations (totalling 83 individuals) representing the range of the species. In spite of the limited distribution of this endemic taxon, high levels in percentage of polymorphic bands (PPB), gene diversity (H _S and H _T) and Shannon??s information index (I) were detected both at population (PPB?=?60?%, H _S?=?0.1853, I?=?0.2836) and at species level (PPB?=?100?%, H _T?=?0.2415, I?=?0.3871). The coefficient of genetic differentiation among population (G??_ST) was 0.1935, while the level of gene flow (N?? _m) was estimated to be 2.0832. AMOVA analysis identified a genetic variation within populations of 83?% of the total. Bayesian clustering methods assigned individuals to two geographical groups partly found within the same population, probably due to a high rate of genetic exchange among its populations. Conservation measures are suggested on the basis of the genetic diversity detected to ensure an effective protection for this endemic species.     

RAPD variation within and among natural populations of outcrossing willow-leaved foxglove (Digitalis obscura L.)

 Random amplified polymorphic DNA (RAPD) markers were used to assess levels and patterns of genetic diversity in Digitalis obscura L. (Scrophulariaceae), an outcrossing cardenolide-producing medicinal plant species. A total of 50 plants from six natural populations on the Iberian Peninsula were analysed by six arbitrarily chosen decamer primers resulting in 96 highly reproducible RAPD bands. To avoid bias in parameter estimation, analyses of population genetic structure were restricted to bands (35 of 96) whose observed frequencies were less than 1–3/n in each population. The analysis of molecular variance (AMOVA) with distances among individuals corrected for the dominant nature of RAPDs (genotypic analysis) showed that most of the variation (84.8%) occurred among individuals within populations, which is expected for an outcrossing organism. Of the remaining variance, 9.7% was attributed to differences between regions, and 5.5% for differences among populations within regions. Estimates of the Wright, Weir and Cockerham and Lynch and Milligan F_ST from null-allele frequencies corroborated AMOVA partitioning and provided significant evidence for population differentiation in D. obscura. A non-parametric test for the homogeneity of molecular variance (HOMOVA) also showed significant differences in the amount of genetic variability present in the six populations. UPGMA cluster analyses, based on Apostol genetic distance, revealed grouping of some geographically proximate populations. Nevertheless, a Mantel test did not give a significant correlation between geographic and genetic distances. This is the first report of the partitioning of genetic variability within and between populations of D. obscura and provides important baseline data for optimising sampling strategies and for conserving the genetic resources of this medicinal species. Received: 7 September 1998 / Accepted: 28 November 1998     

Genetic variation and phylogeographic history of <Emphasis Type="Italic">Picea likiangensis</Emphasis> revealed by RAPD markers

Repeated cycles of retreat and recolonization during the Quaternary ice ages are thought to have greatly influenced current species distributions and their genetic diversity. It remains unclear how this climatic oscillation has affected the distribution of genetic diversity between populations of wind-pollinated conifers in the Qinghai-Tibetan region. In this study, we investigated the within-species genetic diversity and phylogenetic relationships of Picea likiangensis, a dominant forest species in this region using polymorphic DNA (RAPD) markers. Our results suggest that this species has high overall genetic diversity, with 85.42% of loci being polymorphic and an average expected heterozygosity (H _E) of 0.239. However, there were relatively low levels of polymorphism at population levels and the differences between populations were not significant, with percentages of polymorphic bands (PPB) ranging from 46.88 to 69.76%, Nei’s gene diversity (H _E) from 0.179 to 0.289 and Shannon’s indices (Hpop) from 0.267 to 0.421. In accordance with our proposed hypothesis, a high level of genetic differentiation among populations was detected based on Nei’s genetic diversity (G _ST = 0.256) and AMOVA analysis (Phi _st = 0.236). Gene flow between populations was found to be limited (Nm = 1.4532) and far lower than reported for other conifer species with wide distribution ranges from other regions. No clusters corresponding to three morphological varieties found in the south, north and west, respectively, were detected in either UPGMA or PCO analyses. Our results suggest that this species may have had different refugia during the glacial stages in the southern region and that the northern variety may have multiple origins from these different refugia.     

Genetic Variation in the Endangered Endemic Species Cycas fairylakea (Cycadaceae) in China and Implications for Conservation

Cycas fairylakea is an endangered endemic species in China. Genetic diversity within and among four natural populations of this species in China was investigated using amplified fragments length polymorphism (AFLP). A moderate to low level of intraspecific genetic diversity was detected in this species (at population level: P = 39.57 %, H₀ = 0.244; at species level: P = 60.22%, H₀ = 0.356). The among-population component accounted for, respectively, 25.7 and 31.5% of the genetic variation, according to AMOVA and Shannon’s index, indicating most of the genetic variation was found between individuals within populations. All four populations have opposite pyramid age structure, and few coning individuals, which is still decreasing. Possibly because of habitat degradation and environmental pollution, plant diseases and insect pests in the populations were extremely serious, suggesting that the main factors threatening the survival of C. fairylakea populations were not genetic variation, but human activities and the breeding system of this species.     

Genetic structure in peripheral Western European populations of the endangered species Cochlearia pyrenaica (Brassicaceae)

Cochlearia pyrenaica is one of the most endangered plant species in Europe, listed in many European and regional conservation policy documents (e.g. Spain, France, Belgium, Switzerland). To study its genetic structure, define its conservation units and propose a management strategy for this species, amplified fragment length polymorphism markers were used to analyse the genetic diversity within and between five representative populations of the species distribution in Western Europe (Cantabrian Range, North of Spain; Pyrenees, France; Wallonia, Belgium). Low levels of genetic diversity were revealed by the population percentage of polymorphic bands (PPB?=?36.56%), average within-population diversity (H _S?=?0.0990) and genetic diversity within populations (H _pop?=?0.1541), although high levels were reported at species level (PPB?=?81.16%; total genetic diversity for the species, H _T?=?0.0990; and genetic diversity within whole species, H _sp?=?0.2515). The coefficient of genetic differentiation among populations (G _ST) was 0.3869. The analysis of Shannon diversity index in population and for the total data set partitioned (38.72%) and AMOVA (53%) detected a high level of interpopulation diversity, in broad agreement with the result of genetic differentiation analysis. NeighborNet network and principal coordinate analyses clustered the populations in three major groups congruent with geographical regions. Bayesian clustering also confirmed these three distinct genetic clusters. The level of gene flow (Nm) was estimated as 0.3961 individuals per generation among populations, with the genetic identity (I) and genetic distance (D) among populations ranging from 0.8679 to 0.9651 and from 0.0355 to 0.1417, respectively. Therefore, the low levels of genetic variation and high divergence of regional gene pools indicate that there is a need to protect each disjunct region of Western Europe.     

Genetic structure and differentiation of four populations of Afghan Pika (Ochotona rufescens) in Iran based on mitochondrial cytochrome b gene

The population genetics of the Afghan Pika (Ochotona rufescens) was studied in Northern Khorasan Province, Iran. For prediction of the genetic differentiation of four populations, the DNA of mitochondrial cytochrome b of 32 individuals from four areas was sequenced and a Bayesian analysis based on the HKY model was constructed. In total, 15 polymorphic sites, 1125 conserved sites (98.7%) and 14 different haplotypes were found. The phylogenetic tree resulting from the Bayesian analysis and network analysis showed that all samples were clustered in two major groups and the haplotypes of the four populations did not separate geographically. An analysis of molecular variance (AMOVA) indicated that a large majority of the genetic variance was due to the variance within populations. The results of fixation indices showed significant genetic structure among populations in both methods. The pairwise F_st revealed that two northern populations have a significant genetic differentiation from two southern populations, but no significance pairwise F_st value was demonstrated between the closed populations. Nei's genetic distances between closed populations were not significant, while significant values occurred between distant populations. It seems that there is not a major discontinuity between populations of Afghan Pika based on cyt-b mitochondrial gene. However, phylogenetic analysis did not separate populations and a large majority of the genetic variance was found within populations. However, AMOVA analysis showed a significant level of genetic structure among populations (p?<?0.001) and between groups (p?<?0.5). It seems that these results suggest shallow genetic differentiation between populations of different geographic groups.     

Local differentiation amidst extensive allele sharing in Oryza nivara and O. rufipogon

Genetic variation patterns within and between species may change along geographic gradients and at different spatial scales. This was revealed by microsatellite data at 29 loci obtained from 119 accessions of three Oryza series Sativae species in Asia Pacific: Oryza nivara Sharma and Shastry, O. rufipogon Griff., and O. meridionalis Ng. Genetic similarities between O. nivara and O. rufipogon across their distribution are evident in the clustering and ordination results and in the large proportion of shared alleles between these taxa. However, local‐level species separation is recognized by Bayesian clustering and neighbor‐joining analyses. At the regional scale, the two species seem more differentiated in South Asia than in Southeast Asia as revealed by F_ST analysis. The presence of strong gene flow barriers in smaller spatial units is also suggested in the analysis of molecular variance (AMOVA) results where 64% of the genetic variation is contained among populations (as compared to 26% within populations and 10% among species). Oryza nivara (H_E = 0.67) exhibits slightly lower diversity and greater population differentiation than O. rufipogon (H_E = 0.70). Bayesian inference identified four, and at a finer structural level eight, genetically distinct population groups that correspond to geographic populations within the three taxa. Oryza meridionalis and the Nepalese O. nivara seemed diverged from all the population groups of the series, whereas the Australasian O. rufipogon appeared distinct from the rest of the species.