Evidence of population bottleneck in Astragalus michauxii (Fabaceae), a narrow endemic of the southeastern United States

Genetic factors such as decreased genetic diversity and increased homozygosity can have detrimental effects on rare species, and may ultimately limit potential adaptation and exacerbate population declines. The Gulf and Atlantic Coastal Plain physiographic region has the second highest level of endemism in the continental USA, but habitat fragmentation and land use changes have resulted in catastrophic population declines for many species. Astragalus michauxii (Fabaceae) is an herbaceous plant endemic to the region that is considered vulnerable to extinction, with populations generally consisting of fewer than 20 individuals. We developed eight polymorphic microsatellites and genotyped 355 individuals from 24 populations. We characterized the population genetic diversity and structure, tested for evidence of past bottlenecks, and identified evidence of contemporary gene flow between populations. The mean ratios of the number of alleles to the allelic range (M ratio) across loci for A. michauxii populations were well below the threshold of 0.68 identified as indicative of a past genetic bottleneck. Genetic diversity estimates were similar across regions and populations, and comparable to other long-lived perennial species. Within-population genetic variation accounted for 92 % of the total genetic variation found in the species. Finally, there is evidence for contemporary gene flow among the populations in North Carolina. Although genetic factors can threaten rare species, maintaining habitats through prescribed burning, in concert with other interventions such as population augmentation or (re)introduction, are likely most critical to the long term survival of A. michauxii.     

Floral scent composition predicts bee pollination system in five butterfly bush (Buddleja,Scrophulariaceae) species

Traditionally, plant–pollinator interactions have been interpreted as pollination syndrome. However, the validity of pollination syndrome has been widely doubted in modern studies of pollination ecology. The pollination ecology of five Asian Buddleja species, B. asiatica, B. crispa, B. forrestii, B. macrostachya and B. myriantha, in the Sino‐Himalayan region in Asia, flowering in different local seasons, with scented inflorescences were investigated during 2011 and 2012. These five species exhibited diverse floral traits, with narrow and long corolla tubes and concealed nectar. According to their floral morphology, larger bees and Lepidoptera were expected to be the major pollinators. However, field observations showed that only larger bees (honeybee/bumblebee) were the primary pollinators, ranging from 77.95% to 97.90% of total visits. In this study, floral scents of each species were also analysed using coupled gas chromatography and mass spectrometry (GC‐MS). Although the five Buddleja species emitted differentiated floral scent compositions, our results showed that floral scents of the five species are dominated by substances that can serve as attractive signals to bees, including species‐specific scent compounds and principal compounds with larger relative amounts. This suggests that floral scent compositions are closely associated with the principal pollinator assemblages in these five species. Therefore, we conclude that floral scent compositions rather than floral morphology traits should be used to interpret plant–pollinator interactions in these Asian Buddleja species.     

Differentiation among populations of the Rhodiola iremelica Boriss. (Grassulaceae) in the Southern Urals

Isoenzyme markers and polyacrylamide gel electrophoresis have been used to study the genetic structure of populations of Rhodiola iremelica Boriss. (Grassulaceae), a southern Ural endemic protected by the state and included in the Red Data Book of Bashkortostan Republic. A relatively large genetic variation at the species level has been found. The subdivision among populations (F _ST = 0.115) is higher than in most cross-pollination angiosperms. No consistent pattern has been observed in the spatial distribution of its genetic variation. The relatively high differentiation among samples of R. iremelica characterized by small effective population sizes, may be accounted for by genetic drift, inbreeding, and a restricted gene flow. To preserve the population gene pool, in situ protection of the species in nature is insufficient. It seems advisable to create synthetic populations ex situ and reintroduce them into nature.     

Comparison of floral scent between and within Buddleja fallowiana and Buddleja officinalis (Scrophulariaceae)

Floral scents are important olfactory signals for communication between plants and pollinators. Several studies have focused on inter-specific variation of floral scents, but little is known about the intra-specific variation, especially in some polychromic species. In this study, we investigated the floral scent compositions of Buddleja fallowiana and Buddleja officinalis in situ by dynamic headspace collection and coupled GC–MS. Variations of scent compositions within and between populations as well as among species were compared. In spite of substantial intra- and inter-population variability, B. fallowiana and B. officinalis were clearly differentiated in their scent profiles. In B. fallowiana, obvious differentiation was found between studied populations, while all investigated populations in B. officinalis are part of a metapopulation. These high intra-specific variations are discussed in relation to the introgression through hybridization and founder effects from different populations.     

Assessment of genetic diversity and population structure of Bergenia stracheyi (Saxifragaceae) in the Western Himalaya (India)

Genetic diversity and population structure in Bergenia stracheyi, a threatened medicinal herb in the Western Himalaya of India was analysed using directed amplification of minisatellite DNA (DAMD) and inter simple sequence repeats (ISSR) markers. A total of 41 accessions of B. stracheyi representing three populations (Khillenmarg –KLM, Jalori Pass-JLP and Rohtang-RTG) were considered in the present study. The cumulative data analysis for 26 (10 DAMD + 16 ISSR) markers revealed 87.1% polymorphism. The maximum inter-population genetic distance was found between KLM and JLP, whereas the minimum genetic distance was found between RTG and JLP populations. The analysis of molecular variance (AMOVA) revealed maximum percentage of variation among individuals within populations (75%) than among the populations (25%). Clustering pattern of the three sample populations in STRUCTURE and PCoA analyses showed high genetic variation at population level. The present study revealed that distribution patterns, high altitudinal ranges, high habitat specificity, relatively high gene flow, small and isolated population size have shaped the current population structure of B. stracheyi in the Western Himalayan region. DAMD and ISSR markers have provided significant insights into characterization of B. stracheyi populations, and facilitate selection of appropriate accessions for further utilization in conservation and bioprospecting programmes.     

Landscape heterogeneity predicts gene flow in a widespread polymorphic bumble bee, Bombus bifarius (Hymenoptera: Apidae)

Bombus bifarius is a widespread bumble bee that occurs in montane regions of western North America. This species has several major color pattern polymorphisms and shows evidence of genetic structuring among regional populations, and the taxonomic status of regional populations has repeatedly been debated. We test whether observed structure is evidence for discrete gene flow barriers that might indicate isolation or instead reflects clinal variation associated with spatially limited dispersal in a complex landscape. We first consider color pattern variation and identify geographical patterns of B. bifarius color variation using cluster analysis. We then use climate data and a comprehensive set of B. bifarius natural history records with an existing genetic data set to model the distribution of environmentally suitable habitat in western North America and predict pathways of potential gene flow using circuit theory. Resistance distances among populations that incorporate environmental suitability information predict patterns of genetic structure much better than geographic distance or Bayesian clustering alone. Results suggest that there may not be barriers to gene flow warranting further taxonomic considerations, but rather that the arrangement of suitable habitat at broad scales limits dispersal sufficiently to explain observed levels of population differentiation in B. bifarius.     

Environmental Variables Explain Genetic Structure in a Beetle-Associated Nematode

The distribution of a species is a complex expression of its ecological and evolutionary history and integrating population genetic, environmental, and ecological data can provide new insights into the effects of the environment on the population structure of species. Previous work demonstrated strong patterns of genetic differentiation in natural populations of the hermaphroditic nematode Pristionchus pacificus in its La Réunion Island habitat, but gave no clear understanding of the role of the environment in structuring this variation. Here, we present what is to our knowledge the first study to statistically evaluate the role of the environment in shaping the structure and distribution of nematode populations. We test the hypothesis that genetic structure in P. pacificus is influenced by environmental variables, by combining population genetic analyses of microsatellite data from 18 populations and 370 strains, with multivariate statistics on environmental data, and species distribution modelling. We assess and quantify the relative importance of environmental factors (geographic distance, altitude, temperature, precipitation, and beetle host) on genetic variation among populations. Despite the fact that geographic populations of P. pacificus comprise vast genetic diversity sourced from multiple ancestral lineages, we find strong evidence for local associations between environment and genetic variation. Further, we show that significantly more genetic variation in P. pacificus populations is explained by environmental variation than by geographic distances. This supports a strong role for environmental heterogeneity vs. genetic drift in the divergence of populations, which we suggest may be influenced by adaptive forces.     

Regional genetic differentiation in Western Australian sandalwood (<Emphasis Type="Italic">Santalum spicatum</Emphasis>) as revealed by nuclear RFLP analysis

Western Australian sandalwood, Santalum spicatum, is widespread in the semi-arid and arid regions of Western Australia, and there is some morphological variation suggestive of two ecotypes. The level and structuring of genetic diversity within the species was investigated using anonymous nuclear RFLP loci. Santalum spicatum showed moderate levels of genetic diversity compared to other Australian tree species. The northern populations in the arid region showed greater levels of diversity and less population differentiation than the southern populations in the semi-arid region due to differences in the distribution of rare alleles. Equilibrium between drift and gene flow in the northern populations indicated that they have been established for a long period of time with stable conditions conducive to gene flow. In contrast, the southern populations showed a relationship between drift and gene flow indicative of a pattern of fragmentation and isolation where drift has greater effect than gene flow. The different patterns of diversity suggest that the ecotypes in the two regions have been subject to differences in the relative influences of drift and gene flow during their evolutionary history.Communicated by D.B. Neale     

Progeny-testing of full-sibs IBD in a SSC2 QTL region highlights epistatic interactions for fatness traits in pigs

Background

Tetraena mongolica (Zygophyllaceae), an endangered endemic species in western Inner Mongolia, China. For endemic species with a limited geographical range and declining populations, historical patterns of demography and hierarchical genetic structure are important for determining population structure, and also provide information for developing effective and sustainable management plans. In this study, we assess genetic variation, population structure, and phylogeography of T. mongolica from eight populations. Furthermore, we evaluate the conservation and management units to provide the information for conservation.

Results

Sequence variation and spatial apportionment of the atp B- rbc L noncoding spacer region of the chloroplast DNA were used to reconstruct the phylogeography of T. mongolica. A total of 880 bp was sequenced from eight extant populations throughout the whole range of its distribution. At the cpDNA locus, high levels of genetic differentiation among populations and low levels of genetic variation within populations were detected, indicating that most seed dispersal was restricted within populations.

Conclusions

Demographic fluctuations, which led to random losses of genetic polymorphisms from populations, due to frequent flooding of the Yellow River and human disturbance were indicated by the analysis of BEAST skyline plot. Nested clade analysis revealed that restricted gene flow with isolation by distance plus occasional long distance dispersal is the main evolutionary factor affecting the phylogeography and population structure of T. mongolica. For setting a conservation management plan, each population of T. mongolica should be recognized as a conservation unit.     

Patterns of population genetic structure and diversity across bumble bee communities in the Pacific Northwest

Patterns of genetic structure and diversity are largely mediated by a species’ ecological niche and sensitivity to climate variation. Some species with narrow ecological niches have been found to exhibit increased population differentiation, limited gene flow across populations, and reduced population genetic diversity. In this study, we examine patterns of population genetic structure and diversity of four bumble bee species that are broadly sympatric, but do not necessarily inhabit the same ecological niche in the Pacific Northwest of the United States. Testing for the effect of isolation by geographic distance (IBD) with linearized F _st and D _est found that Bombus sylvicola and B. mixtus exhibited significant IBD across populations. In contrast, both B. melanopygus and B. flavifrons, two species that are distributed across a broad elevation gradient, exhibited no IBD, a result further corroborated by Bayesian a priori population assignment tests. Furthermore, we discovered that B. sylvicola populations distributed on the Olympic Peninsula have significantly less average allelic diversity than populations distributed in the Cascade Mountains. Our results suggest that populations distributed in the Olympic Mountains represent a distinct genetic cluster relative to the Cascade Mountains, with B. sylvicola and B. mixtus likely experiencing the greatest degree of population genetic differentiation relative to B. flavifrons and B. melanopygus. While bumble bees are known to co-exist across a diversity of habitats, our results demonstrate that underlying population genetic structure and diversity may not necessarily be similar across species, and are largely governed by their respective niches.     

Conservation genetics of the endangered terrestrial orchid Liparis japonica in Northeast China based on AFLP markers

Levels of genetic diversity and population genetic structure of the rare, endangered terrestrial orchid Liparis japonica were examined for eight natural populations (n = 185) in Northeast China using six AFLP primer pairs, where this species has experienced severe habitat loss and fragmentation. Based on 406 DNA bands, a high level of genetic diversity was found at the species level with the PPB of 85.47 %, while the genetic diversity at the population level was low (PPB = 47.48 %). A significantly high degree of population differentiation was found with 42.69 % variation existed among populations as measured by AMOVA, indicating potential restricted gene flow. The genetic distances between populations were independent of the corresponding geographic distances, and the genetic relationship of individuals had no significant correlation with their spatial distribution. The restricted gene flow might be impacted by reduced population size, habitat destruction and fragmentation. The results in this study suggested that habitat protection and keeping a stable environment are critical for the conservation of L. japonica species.     

Geoclimatic factors influence the population genetic connectivity of Incarvillea arguta (Bignoniaceae) in the Himalaya–Hengduan Mountains biodiversity hotspot

Geoclimatic factors related to the uplift of the Himalaya and the Quaternary climatic oscillations influence the population genetic connectivity in the Himalaya–Hengduan Mountains (HHM) biodiversity hotspot. Therefore, to explore the relative roles played by these two factors, we examined the population dynamics and dispersal corridors of Incarvillea arguta (Royle) Royle incorporating ensemble species distribution modelling (SDM). Thirty‐seven populations were genotyped using plastid chloroplast DNA and low copy nuclear gene (ncpGS) sequences. Phylogeographic analysis was carried out to reveal the genetic structure and lineage differentiation. Ensemble SDMs were carried out for distributional change in the last glacial maximum, present, and future. Finally, the least cost path method was used to trace out possible dispersal corridors. The haplotypes were divided into four clades with strong geographical structure. The late Miocene origin of I. arguta in the western Himalaya ca. 7.92 Ma indicates lineage diversification related to the uplift of the HHM. The variability in habitat connectivity revealed by SDM is due to change in suitability since the Pleistocene. A putative dispersal corridor was detected along the drainage systems and river valleys, with strong support in the eastern Hengduan Mountains group. Our results support the signature of geoclimatic influence on population genetic connectivity of I. arguta in the HHM. We proposed that the major drainage systems might have assisted the rapid dispersal of isolated riverine plant species I. arguta in the HHM. The population genetic connectivity, using the fine‐tuned ensemble SDMs, enables scientists and policymakers to develop conservation strategies for the species gene pool in the HHM biodiversity hotspots.     

RAD sequencing of common whelk,Buccinum undatum,reveals fine‐scale population structuring in Europe and cryptic speciation within the North Atlantic

Buccinum undatum is a subtidal gastropod that exhibits clear spatial variation in several phenotypic shell traits (color, shape, and thickness) across its North Atlantic distribution. Studies of spatial phenotypic variation exist for the species; however, population genetic studies have thus far relied on a limited set of mitochondrial and microsatellite markers. Here, we greatly expand on previous work by characterizing population genetic structure in B. undatum across the North Atlantic from SNP variation obtained by RAD sequencing. There was a high degree of genetic differentiation between Canadian and European populations (Iceland, Faroe Islands, and England) consistent with the divergence of populations in allopatry (F _ST > 0.57 for all pairwise comparisons). In addition, B. undatum populations within Iceland, the Faroe Islands, and England are typified by weak but significant genetic structuring following an isolation‐by‐distance model. Finally, we established a significant correlation between genetic structuring in Iceland and two phenotypic traits: shell shape and color frequency. The works detailed here enhance our understanding of genetic structuring in B. undatum and establish the species as an intriguing model for future genome‐wide association studies.     

Population genetic differentiation of Schisandra chinensis and Schisandra sphenanthera as revealed by ISSR analysis

Schisandra chinensis (Turcz.) Baill. and Schisandra sphenanthera Rehd. et Wils. are well-known Chinese medicinal plants. The population genetic variation of the two species was studied using inter simple sequence repeat (ISSR) molecular markers. High levels of genetic diversity are revealed in both S. chinensis (P = 88.36%, h = 0.2894, I = 0.4396) and S. sphenanthera (P = 84.09%, H = 0.2782, I = 0.4280). However, the population genetic differentiation is significantly different between the two species. The S. sphenanthera harbors as high as 27% of the genetic variation among populations but 73% within populations, whereas in S. chinensis 17% of the genetic variation occurs among populations and 83% within populations. Both significant (P < 0.05) heterozygosity excess and shifted mode of allele frequency distribution are detected in four out of six populations of S. chinensis and one out of five populations of S. sphenanthera, suggesting the occurrence of recent population bottlenecks in the two species. The different patterns of genetic variation in S. chinensis and S. sphenanthera are discussed in relation to their differences in pollination mechanism, geographic distribution and historical events, and the level of gene flow and genetic drift.     

陕西省林麝mtDNA D-loop区序列结构和种群遗传多样性

林麝(Moschus berezovskii)曾广泛分布于中国,由于盗猎和栖息地缩小,秦岭地区野生种群数量迅速下降,圈养繁殖种群已成立了几十年,但大多数圈养种群的遗传背景不清,种群规模增长非常缓慢。为了给这一物种的保护和管理提供有用的信息,调查了陕西省林麝1个圈养种群3个野生种群线粒体DNA(mt DNA)D-Loop 632 bp片段的遗传多样性和种群结构。在69个个体中其碱基组成为A+T的平均含量63.2%高于G+C含量36.8%,共检测到变异位点171个(约占总位点数的27.05%)。核苷酸多样性(Pi)为0.04424,平均核苷酸差异数(K)为19.908。69个个体分属32个单倍型,单倍型间的平均遗传距离(P)为0.070。32个单倍型构建的NJ系统树聚为3个分支,4个林麝群体中的单倍型是随机分布的。4个群体的平均遗传距离为0.043(标准误SE为0.005),凤县养殖场群体与留坝和陇县群体的亲缘关系较远。单倍型间的平均遗传距离为0.043,可见其遗传分化尚未达到种群分化的水平。结果表明,陕西省林麝群体mt DNA D-loop区序列存在着较丰富的变异和遗传多样性,凤县野生群体和凤县养殖场群体的核苷酸多样性和单倍型多样较高,养殖场种群没有出现近亲繁殖及遗传多样性下降的情况。凤县野生群体和凤县养殖场群体两者遗传分化较小,存在着较高的基因流水平。     

Divergent selection along climatic gradients in a rare central European endemic species,Saxifraga sponhemica

Background and Aims The effects of habitat fragmentation on quantitative genetic variation in plant populations are still poorly known. Saxifraga sponhemica is a rare endemic of Central Europe with a disjunct distribution, and a stable and specialized habitat of treeless screes and cliffs. This study therefore used S. sponhemica as a model species to compare quantitative and molecular variation in order to explore (1) the relative importance of drift and selection in shaping the distribution of quantitative genetic variation along climatic gradients; (2) the relationship between plant fitness, quantitative genetic variation, molecular genetic variation and population size; and (3) the relationship between the differentiation of a trait among populations and its evolvability.Methods Genetic variation within and among 22 populations from the whole distribution area of S. sponhemica was studied using RAPD (random amplified polymorphic DNA) markers, and climatic variables were obtained for each site. Seeds were collected from each population and germinated, and seedlings were transplanted into a common garden for determination of variation in plant traits.Key Results In contrast to previous results from rare plant species, strong evidence was found for divergent selection. Most population trait means of S. sponhemica were significantly related to climate gradients, indicating adaptation. Quantitative genetic differentiation increased with geographical distance, even when neutral molecular divergence was controlled for, and Q_ST exceeded F_ST for some traits. The evolvability of traits was negatively correlated with the degree of differentiation among populations (Q_ST), i.e. traits under strong selection showed little genetic variation within populations. The evolutionary potential of a population was not related to its size, the performance of the population or its neutral genetic diversity. However, performance in the common garden was lower for plants from populations with reduced molecular genetic variation, suggesting inbreeding depression due to genetic erosion.Conclusions The findings suggest that studies of molecular and quantitative genetic variation may provide complementary insights important for the conservation of rare species. The strong differentiation of quantitative traits among populations shows that selection can be an important force for structuring variation in evolutionarily important traits even for rare endemic species restricted to very specific habitats.     

Molecular Analysis of the Freshwater Prawn Macrobrachium olfersii (Decapoda,Palaemonidae) Supports the Existence of a Single Species throughout Its Distribution

Macrobrachium olfersii is an amphidromous freshwater prawn, widespread along the eastern coasts of the Americas. This species shows great morphological modifications during ontogenesis, and several studies have verified the existence of a wide intraspecific variation. Because of this condition, the species is often misidentified, and several synonyms have been documented. To elucidate these aspects, individuals of M. olfersii from different populations along its range of distribution were investigated. The taxonomic limit was established, and the degree of genetic variability of this species was described. We extracted DNA from 53 specimens of M. olfersii, M. americanum, M. digueti and M. faustinum, which resulted in 84 new sequences (22 of 16S mtDNA, 45 of Cythocrome Oxidase I (COI) mtDNA, and 17 of Histone (H3) nDNA). Sequences of three genes (single and concatenated) from these species were used in the Maximum Likelihood and Bayesian Inference phylogenetic analyses and COI sequences from M. olfersii were used in population analysis. The genetic variation was evaluated through the alignment of 554 bp from the 16S, 638 bp from the COI, and 338 bp from the H3. The rates of genetic divergence among populations were lower at the intraspecific level. This was confirmed by the haplotype net, which showed a continuous gene flow among populations. Although a wide distribution and high morphological intraspecific variation often suggest the existence of more than one species, genetic similarity of Caribbean and Brazilian populations of M. olfersii supported them as a single species.     

Population Genetic Structure of the Invasive Red Swamp Crayfish in China Revealed by ITS1 Variation

The invasive red swamp crayfish (Procambarus clarkii) provides a valuable opportunity for studying the population genetics of invasive species that disperse rapidly. We analyzed the population genetic structure among 12 populations of the crayfish in China based on the internal transcribed spacer 1 (ITS1) region. The ITS1 of 815 bp aligned across 34 haplotypes; the average GC content was 53.9%. AMOVA showed that intrapopulation variation (95.26%) was much higher than interpopulation variation (4.74%). Genetic differentiation between the Taiwan and mainland populations (F _st = 0.160) was moderate, but the Chinese population (Taiwan and the mainland combined) and an American population were highly differentiated (0.682 and 0.977, respectively). Gene flow between the Chinese and American populations (N _m = 0.006 and 0.117, respectively) was lower than that between Taiwan and the mainland (1.536). Phylogenetic trees showed that three major genealogical clusters matched the sample locations well, suggesting that genetic differentiation is created largely by geographic isolation.     

Recent low levels of differentiation in the native <Emphasis Type="Italic">Bombus ephippiatus</Emphasis> (Hymenoptera: Apidae) along two Neotropical mountain-ranges in Guatemala

Recent anthropogenic fragmentation has led to population differentiation threatening viability of many species, including species specialized on mountainous ecosystems. Bombus ephippiatus, a widespread species mostly found in mountains in the Neotropics, seems to use the highlands as island, and deforested lowland areas may represent barriers to their dispersal, leading to isolation and potentially loss of genetic diversity. Yet, lack of knowledge of its population structure does not allow adequate management and conservation. To fill this knowledge gap, we assessed the population structure and inferred dispersion of B. ephippiatus in two mountain-ranges in Guatemala (Volcanic Chain and Sierra de las Minas). This region is characterized by high topographic variation and considerable deforestation strain. We analyzed the effects of elevation and land-use on genetic differentiation of B. ephippiatus populations and inferred its demography in the region. Our results suggest that B. ephippiatus is able to disperse long distances across most landscape types, reflected by its high genetic diversity, high effective population size, considerable gene flow, low population differentiation, as well as the lack of isolation by distance. Hence, B. ephippiatus may be a resilient species for the provision of pollination services. However, we detected a subtle divergence of B. ephippiatus into two clusters, of which Sierra de las Minas has been identified as a regional hotspot of genetic and species endemism. Yet, differentiation is very recent and hence likely caused by lowland deforestation. The combined effects of current forest cover and elevation partially explain the observed subtle patterns of differentiation suggesting that the maintenance of suitable habitat is crucial to ensure population connectivity of this keystone pollinator.     

Population structure affected by excess gene flow in self-pollinating Elymus nutans and E. burchan-buddae (Triticeae: Poaceae)

Seed-mediated gene flow can considerably affect population genetic structure of strictly self-pollinating species, but little is known on the extent and nature of such gene flow among pastoral plant populations. Molecular fingerprints provide a powerful tool to address the relevant issues. Genetic structure of 22 populations of two self-pollinating pasture species, Elymus nutans and E. burchan-buddae, collected from various altitudes of the Qinghai–Tibetan Plateau was studied using fluorescence-based amplified fragment length polymorphism technique. Analysis of molecular variance revealed 42.97% and 37.63% among-population variation for the two Elymus species, respectively, indicating that the majority of the total variation presented within populations. This result contradicts the common genetic variation pattern for a selfing plant species: lower genetic variation within populations. Further analysis suggested higher level of gene flow among populations within the same region than among different regions across the sampled area for the two Elymus species. STRUCTURE analyses of the Elymus populations indicated an evident admixture genetic structure, particularly among neighboring populations from the same region, supporting the hypothesis of considerable seed dispersal among populations. The excess within-region gene flow of E. nutans and E. burchan-buddae might be caused by grazing animals that promote seed dispersal when moved around the pastoral lands during foraging. The among-population gene flow promulgated by grazing animals may promote the maintenance of genetic diversity in the pasture species, particularly in small and fragmented populations within a given region.