Population genetic variability and structure of Elymus breviaristatus (Poaceae: Triticeae) endemic to Qinghai–Tibetan Plateau inferred from SSR markers

Characterization and Application of EST-SSR Markers Developed from Transcriptome Sequences in Elymus breviaristatus (Poaceae: Triticeae)

BACKGROUND

Elymus L. is the largest genus in the Triticeae tribe. Most species in this genus are highly stress resistant, with excellent forage value. Elymus breviaristatus, a rare species endemic to the Qinghai-Tibet Plateau (QTP), is declining due to habitat fragmentation. However, genetic data for E. breviaristatus are limited, with expressed sequence tag (EST) markers being particularly rare, hampering genetic studies and protection measures.

RESULTS

We obtained 9.06 Gb clean sequences from the transcriptome of E. breviaristatus, generating 171,522 unigenes, which were assembled and functionally annotated against five public databases. We identified 30,668 SSRs in the E. breviaristatus transcriptome, from which 103 EST-SSR primer pairs were randomly selected. Of these, 58 pairs of amplified products of the expected size, and 18 of the amplified products were polymorphic. Model-based Bayesian clustering, the unweighted pair group method with arithmetic average (UPGMA), and principal coordinate analysis (PCoA) of 179 wild E. breviaristatus in 12 populations using these EST-SSRs were generally consistent, grouping the 12 populations into two major clades. Analysis of molecular variance (AMOVA) found 70% of the genetic variation among the 12 populations and 30% within the populations, indicating a high level of genetic differentiation (or low gene exchange) among the 12 populations. The transferability of the 58 successful EST-SSR primers to 22 related hexaploid species was 86.2-98.3%. UPGMA analysis generally grouped species with similar genome types together.

CONCLUSIONS

Here, we developed EST-SSR markers from the transcriptome of E. breviaristatus. The transferability of these markers was evaluated, and the genetic structure and diversity of E. breviaristatus were explored. Our results provide a basis for the conservation and management of this endangered species, and the obtained molecular markers represent valuable resources for the exploration of genetic relationships among species in the Elymus genus.

Drivers of population divergence and genetic variation in Elymus breviaristatus (Keng) Keng f. (Poaceae: Triticeae), an endemic perennial herb of the Qinghai-Tibet plateau

Elymus breviaristatus, a rare grass species with excellent resistance and ecological importance, is narrowly distributed on the Qinghai-Tibet plateau. Populations of E. breviaristatus are declining due to habitat fragmentation, and thus far, characteristics of genetic differentiation and adaptive responses to climate change remain poorly understood in this species. Here, we explored the genetic structure of 18 natural populations (269 individuals) in the transition zone between Tibet and the Hengduan Mountains using 15 expressed sequence tag (EST)-SSR primer pairs and identified possible barriers to gene flow that might have caused genetic discontinuities. Additional analyses were performed to identify the environmental factors affecting genetic diversity and to test whether the patterns of genetic variation among populations were more consistent with the isolation by distance (IBD) or isolation by environment (IBE) model. Multiple measures of genetic diversity revealed that intra-population genetic variation was low, while inter-population genetic variation was high. Clustering, structure, and principal coordinate analyses identified three genetic groups: (a) Eastern Qamdo, (b) Nagqu and Western Qamdo, and (c) Lhasa and Nyingchi. A clear physical barrier to gene flow was formed by the Yarlung Zangbo Grand Canyon and the Tanggula Mountains. We found that both IBD and IBE contributed to the observed patterns of genetic variation, and the IBE model played a leading role. In addition, precipitation-related variables, soil phosphorus content and soil K:P ratio significantly affected population genetic variation. Overall, our results emphasized the genetic fragility of E. breviaristatus populations and showed that this species requires attention, as future climate changes and human activities may further threaten its survival. In addition, the genetic differences among E. breviaristatus populations should be considered when formulating conservation measures for E. breviaristatus populations in the study area.

Assessment of the Genetic Relationship and Population Structure in Oil-Tea Camellia Species Using Simple Sequence Repeat (SSR) Markers

Oil-tea camellia trees, the collective term for a class of economically valuable woody oil crops in China, have attracted extensive attention because of their rich nutritional and pharmaceutical value. This study aimed to analyze the genetic relationship and genetic diversity of oil-tea camellia species using polymorphic SSR markers. One-hundred and forty samples of five species were tested for genetic diversity using twenty-four SSR markers. In this study, a total of 385 alleles were identified using 24 SSR markers, and the average number of alleles per locus was 16.0417. The average Shannon's information index (I) was 0.1890, and the percentages of polymorphic loci (P) of oil-tea camellia trees were 7.79-79.48%, indicating that oil-tea camellia trees have low diversity. Analysis of molecular variance (AMOVA) showed that the majority of genetic variation (77%) was within populations, and a small fraction (23%) occurred among populations. Principal coordinate analysis (PCoA) results indicated that the first two principal axes explained 7.30% (PC1) and 6.68% (PC2) of the total variance, respectively. Both UPGMA and PCoA divided the 140 accessions into three groups. Camellia oleifera clustered into one class, Camellia vietnamensis and Camellia gauchowensis clustered into one class, and Camellia crapnelliana and Camellia chekiangoleosa clustered into another class. It could be speculated that the genetic relationship of C. vietnamensis and C. gauchowensis is quite close. SSR markers could reflect the genetic relationship among oil-tea camellia germplasm resources, and the results of this study could provide comprehensive information on the conservation, collection, and breeding of oil-tea camellia germplasms.

Fingerprint identification of white clover cultivars based on SSR molecular markers

White clover (Trifolium repens L.) is an important perennial legume forage with high productivity and quality. To strengthen the basic research on the genetic characteristics, fingerprint identification and adaptability of white clover germplasm resources, Simple sequence repeat (SSR) molecular markers were applied to 10 white clover cultivars to assess the genetic diversity and related lines of white clover at the molecular level in order to lay a theoretical foundation for the selection of high-quality seeds and cultivars of white clover. A total of 120 different bands were amplified by 29 pairs of SSR primers with good polymorphism, of which 103 (89.5%) were polymorphic. Meanwhile, the PIC of each primer was 0.181-0.588, with an average of 0.329. Analysis of molecular variance revealed that 57% of the genetic variation occurred within cultivars and 43% occurred among cultivars. The results of cluster analysis and the principal coordinate analysis revealed that the parental relationships of the 10 cultivars, with the 'Purple' cultivar very distantly related to the other 9 cultivars and the closest parental relationship between 'Ladino' and 'Sulky'. The fingerprints constructed by three representative primers (gtrs679, gtrs319, and gtrs678) have a strong identification ability. In summary, the SSR markers had good polymorphism and could be used for DNA fingerprint analysis of white clover cultivars.

Revelation of genetic diversity and structure of wild Elymus excelsus (Poaceae: Triticeae) collection from western China by SSR markers

Hosting unique and important plant germplasms, the Qinghai-Tibet Plateau (QTP), as the third pole of the world, and Xinjiang, located in the centre of the Eurasian continent, are major distribution areas of perennial Triticeae grasses, especially the widespread Elymus species. Elymus excelsus Turcz. ex Griseb, a perennial forage grass with strong tolerance to environmental stresses, such as drought, cold and soil impoverishment, can be appropriately used for grassland establishment due to its high seed production. To provide basic information for collection, breeding strategies and utilization of E. excelsus germplasm, microsatellite markers (SSR) were employed in the present study to determine the genetic variation and population structure of 25 wild accessions of E. excelsus from Xinjiang (XJC) and the QTP, including Sichuan (SCC) and Gansu (GSC) of western China. Based on the 159 polymorphic bands amplified by 35 primer pairs developed from three related species, the average values of the polymorphic information content (PIC), marker index (MI), resolving power (Rp), Nei's genetic diversity (H) and Shannon's diversity index (I) of each pair of primers were 0.289, 1.348, 1.897, 0.301 and 0.459, respectively, validating that these SSR markers can also be used for the evaluation of genetic diversity of E. excelsus germplasms, and demonstrating the superior versatility of EST-SSR vs. G-SSR. We found a relatively moderate differentiation (F st  = 0.151) among the XJC, SCC and GSC geo-groups, and it is worth noting that, the intra-group genetic diversity of the SCC group (H e  = 0.197) was greater than that of the GSC (H e  = 0.176) and XJC (H e  = 0.148) groups. Both the Unweighted Pair Group Method with Arithmetic (UPGMA) clustering and principal coordinates analysis (PCoA) divided the 25 accessions into three groups, whereas the Bayesian STRUCTURE analysis suggested that E. excelsus accessions fell into four main clusters. Besides, this study suggested that geographical distance and environmental variables (annual mean precipitation and average precipitation in growing seasons), especially for QTP accessions, should be combined to explain the population genetic differentiation among the divergent geographical regions. These data provided comprehensive information about these valuable E. excelsus germplasm resources for the protection and collection of germplasms and for breeding strategies in areas of Xinjiang and QTP in western China.

Genetic Diversity and Population Divergence of a Rare, Endemic Grass (Elymus breviaristatus) in the Southeastern Qinghai-Tibetan Plateau

Elymus breviaristatus is a grass species only distributed in the southeast of Qinghai-Tibetan Plateau (QTP), which has suffered from serious habitat fragmentation. Therefore, understanding patterns of genetic diversity within and among natural E. breviaristatus populations could provide insight for future conservation strategies. In this study, sequence-related amplified polymorphism markers were employed to investigate the genetic diversity and hierarchical structure of seven E. breviaristatus populations from QTP, China. Multiple measures of genetic diversity indicated that there is low to moderate genetic variation within E. breviaristatus populations, consistent with its presumed mating system. In spite of its rarity, E. breviaristatus presented high genetic diversity that was equivalent to or even higher than that of widespread species. Bayesian clustering approaches, along with clustering analysis and principal coordinate analysis partitioned the studied populations of E. breviaristatus into five genetic clusters. Differentiation coefficients (Fst, GST, etc.) and AMOVA analysis revealed considerable genetic divergence among different populations. BARRIER analyses indicated that there were two potential barriers to gene flow among the E. breviaristatus populations. Despite these patterns of differentiation, genetic distances between populations were independent of geographic distances (r = 0.2197, p = 0.2534), indicating little isolation by distance. Moreover, despite detecting a common outlier by two methods, bioclimatic factors (altitude, annual mean temperature, and annual mean precipitation) were not related to diversity parameters, indicating little evidence for isolation caused by the environment. These patterns of diversity within and between populations are used to propose a conservation strategy for E. breviaristatus.

Genetic variability evaluation and cultivar identification of tetraploid annual ryegrass using SSR markers

Annual ryegrass (Lolium multiflorum) is a widely used cool-season turf and forage grass with high productivity and ornamental characteristics. However, the abundant intra-cultivar genetic variability usually hampers the application of conventional techniques for cultivar identification. The objectives of this study were to: (1) describe an efficient strategy for identification of six tetraploid annual ryegrass cultivars and (2) investigate the genetic diversity based on SSR markers. A total of 242 reliable bands were obtained from 29 SSR primer pairs with an average of 8.3 bands for each primer pair and the average value of polymorphic information content (PIC) was 0.304. The result of analysis of molecular variance (AMOVA) revealed that 81.99% of the genetic variation occurred in within-cultivars and 18.01% among-cultivars. The principal coordinate analysis (PCoA) showed that the first two principal axes explain 8.57% (PC1) and 6.05% (PC2) of total variation, respectively. By using multi-bulk strategy based on different filtering thresholds, the results suggested that bands frequency of 40% could be used as a reliable standard for cultivar identification in annual ryegrass. Under this threshold, 12 SSR primer pairs (00-04A, 02-06G, 02-08C, 03-05A, 04-05B, 10-09E, 12-01A, 13-02H, 13-12D, 14-06F, 15-01C and 17-10D) were detected for direct identification of six tetraploid annual ryegrass cultivars, which could be incorporated into conservation schemes to protect the intellectual property of breeders, ensure purity for consumers, as well as guarantee effective use of cultivars in future.

Genetic Structure and Eco-Geographical Differentiation of Wild Sheep Fescue (Festuca ovina L.) in Xinjiang, Northwest China

Glaciation and mountain orogeny have generated new ecologic opportunities for plants, favoring an increase in the speciation rate. Moreover, they also act as corridors or barriers for plant lineages and populations. High genetic diversity ensures that species are able to survive and adapt. Gene flow is one of the most important determinants of the genetic diversity and structure of out-crossed species, and it is easily affected by biotic and abiotic factors. The aim of this study was to characterize the genetic diversity and structure of an alpine species, Festuca ovina L., in Xinjiang, China. A total of 100 individuals from 10 populations were analyzed using six amplified fragment length polymorphism (AFLP) primer pairs. A total of 583 clear bands were generated, of which 392 were polymorphic; thus, the percentage of polymorphic bands (PPB) was 67.24%. The total and average genetic diversities were 0.2722 and 0.2006 (0.1686-0.2225), respectively. The unweighted group method with arithmetic mean (UPGMA) tree, principal coordinates analysis (PCoA) and Structure analyses revealed that these populations or individuals could be clustered into two groups. The analysis of molecular variance analysis (AMOVA) suggested that most of the genetic variance existed within a population, and the genetic differentiation (Fst) among populations was 20.71%. The Shannon differentiation coefficient (G'st) among populations was 0.2350. Limited gene flow (Nm = 0.9571) was detected across all sampling sites. The Fst and Nm presented at different levels under the genetic barriers due to fragmentation. The population genetic diversity was significant relative to environmental factors such as temperature, altitude and precipitation.