Elymus L. is the largest genus in the tribe Triticeae Dumort., encompassing approximately 150 polyploid perennial species widely distributed in the temperate regions of the world. It is considered to be an important gene pool for improving cereal crops. However, a shortage of molecular marker limits the efficiency and accuracy of genetic breeding for Elymus species. High-throughput transcriptome sequencing data is essential for gene discovery and molecular marker development. We obtained the transcriptome dataset of E. sibiricus, the type species of the genus Elymus, and identified a total of 8871 putative EST-SSRs from 6685 unigenes. Trinucleotides were the dominant repeat motif (4760, 53.66%), followed by dinucleotides (1993, 22.47%) and mononucleotides (1876, 21.15%). The most dominant trinucleotide repeat motif was CCG/CGG (1119, 23.5%). Sequencing of PCR products showed that the sequenced alleles from different Elymus species were homologous to the original SSR locus from which the primer was designed. Different types of tri-repeats as abundant SSR motifs were observed in repeat regions. Two hundred EST-SSR primer pairs were designed and selected to amplify ten DNA samples of Elymus species. Eighty-seven pairs of primer (43.5%) generated clear and reproducible bands with expected size, and showed good transferability across different Elymus species. Finally, thirty primer pairs successfully amplified ninety-five accessions of seventeen Elymus species, and detected significant amounts of polymorphism. In general, hexaploid Elymus species with genomes StStHHYY had a relatively higher level of genetic diversity (H = 0.219, I = 0.330, %P = 63.7), while tetraploid Elymus species with genomes StStYY had low level of genetic diversity (H = 0.182, I = 0.272, %P = 50.4) in the study. The cluster analysis showed that all ninety-five accessions were clustered into three major clusters. The accessions were grouped mainly according to their genomic components and origins. This study demonstrated that transcriptome sequencing is a fast and cost-effective approach to molecular marker development. These EST-SSR markers developed in this study are valuable tools for genetic diversity, evolutionary, and molecular breeding in E. sibiricus, and other Elymus species.

Elymus nutans and Elymus sibiricus are two important perennial forage grasses of the genus Elymus, widely distributed in high altitude regions of Western and Northern China, especially on the Qinghai-Tibetan Plateau. Information on phenotypic and genetic diversity is limited, but necessary for Elymus germplasm collection, conservation, and utilization. In the present study, the phenotypic and genetic differentiation of 73 accessions of the two species were evaluated using 15 phenotypic traits and 40 expressed sequence tag derived simple sequence repeat markers (EST-SSRs). The results showed that only 7.23% phenotypic differentiation (Pst) existed between the two Elymus species based on fifteen quantitative traits. Principal component analysis (PCA) revealed that leaf traits, spike traits, and some seed traits were dominant factors in phenotypic variation. Moreover, 396 (97.8%) and 331 (87.1%) polymorphic bands were generated from 40 EST-SSR primers, suggesting high levels of genetic diversity for the two species. The highest genetic diversity was found in the Northeastern Qinghai-Tibetan Plateau groups. Clustering analysis based on molecular data showed that most accessions of each Elymus species tended to group together. Similar results were described by principal coordinates analysis (PCoA) and structure analysis. The molecular variance analysis (AMOVA) revealed that 81.47% and 89.32% variation existed within the geographical groups for the two species, respectively. Pearson’s correlation analyses showed a strong positive correlation between Nei’s genetic diversity and annual mean temperature. These results could facilitate Elymus germplasm collection, conservation, and future breeding.

Elymus sibiricus, which is a perennial and self-pollinated grass, is the typical species of the genus Elymus, which plays an important role in forage production and ecological restoration. No reports have, so far, systematically described the selection of optimal reference genes for reverse transcriptase quantitative real-time polymerase chain reaction (RT-qPCR) analysis in E. sibiricus. The goals of this study were to evaluate the expression stability of 13 candidate reference genes in different experimental conditions, and to determine the appropriate reference genes for gene expression analysis in E. sibiricus. Five methods including Delta Ct (ΔCt), BestKeeper, NormFinder, geNorm, and RefFinder were used to assess the expression stability of 13 potential reference genes. The results of the RefFinder analysis showed that TBP2 and HIS3 were the most stable reference genes in different genotypes. TUA2 and PP2A had the most stable expression in different developmental stages. TBP2 and PP2A were suitable reference genes in different tissues. Under salt stress, ACT2 and TBP2 were identified as the most stable reference genes. ACT2 and TUA2 showed the most stability under heat stress. For cold stress, PP2A and ACT2 presented the highest degree of expression stability. DNAJ and U2AF were considered as the most stable reference genes under osmotic stress. The optimal reference genes were selected to investigate the expression pattern of target gene CSLE6 in different conditions. This study provides suitable reference genes for further gene expression analysis using RT-qPCR in E. sibiricus.

Background: Elymus sibiricus is an ecologically and economically important perennial, self-pollinated, and allotetraploid (StStHH) grass, widely used for forage production and animal husbandry in Western and Northern China. However, it has low seed yield mainly caused by seed shattering, which makes seed production difficult for this species. The goals of this study were to construct the high-density genetic linkage map, and to identify QTLs and candidate genes for seed-yield related traits. Results: An F₂ mapping population of 200 individuals was developed from a cross between single genotype from “Y1005” and “ZhN06”. Specific-locus amplified fragment sequencing (SLAF-seq) was applied to construct the first genetic linkage map. The final genetic map included 1971 markers on the 14 linkage groups (LGs) and was 1866.35 cM in total. The length of each linkage group varied from 87.67 cM (LG7) to 183.45 cM (LG1), with an average distance of 1.66 cM between adjacent markers. The marker sequences of E. sibiricus were compared to two grass genomes and showed 1556 (79%) markers mapped to wheat, 1380 (70%) to barley. Phenotypic data of eight seed-related traits (2016–2018) were used for QTL identification. A total of 29 QTLs were detected for eight seed-related traits on 14 linkage groups, of which 16 QTLs could be consistently detected for two or three years. A total of 6 QTLs were associated with seed shattering. Based on annotation with wheat and barley genome and transcriptome data of abscission zone in E. sibiricus, we identified 30 candidate genes for seed shattering, of which 15, 7, 6 and 2 genes were involved in plant hormone signal transcription, transcription factor, hydrolase activity and lignin biosynthetic pathway, respectively. Conclusion: This study constructed the first high-density genetic linkage map and identified QTLs and candidate genes for seed-related traits in E. sibiricus. Results of this study will not only serve as genome-wide resources for gene/QTL fine mapping, but also provide a genetic framework for anchoring sequence scaffolds on chromosomes in future genome sequence assembly of E. sibiricus.

Awns, needle-like structures formed on the distal of the lemmas in the florets, are of interest because of their essential roles in seed dispersal, germination and photosynthesis. Previous research has reported the potential benefits of awns in major cereal grasses, yet reports on the agronomic and economic implications of awn length variation in forage grasses remain scarce. This study investigated the variation of awn length among 20 Siberian wildrye populations and the effect of awn length on seed yield and yield components. This work then studied the impact of awn length on seed dispersal and germination. The analyses indicated a high level of awn length variation among populations. Awn length showed a significant influence on harvested seed yield per plant (p < 0.05) mostly driven by interactions between awn length and the majority of seed yield components. Principal component analysis clearly revealed that the final impact of awn length on seed yield depends on the balance of its positive and negative effects on traits determining seed yield. Furthermore, awn length tended to increase seed dispersal distance, although little diversity in the nature of this progression was observed in some populations. Awn length exhibited a significant relationship (p < 0.05) with germination percentage. It also tended to shorten germination duration, although this interaction was not statistically significant. Collectively, these results provide vital information for breeding and agronomic programs aiming to maintain yield in grasses. This is the first report to demonstrate in Siberian wildrye the agronomic impacts of awn length variation.

Elymus nutans and E. sibiricus are two important forage grasses of the genus Elymus. But they are difficult to grow for commercial seed production due to serious seed shattering. We conducted a comparative transcriptome analysis of abscission zone to find possible transcription changes associated with seed shattering, explore candidate genes involved in seed shattering and identify candidate gene-based EST-SSR markers for germplasm evaluation. cDNA libraries from abscission zone (AZ) and non-abscission zone (NAZ) tissues of E. nutans were constructed and sequenced. A total of 111,667 unigenes were annotated and 7644 differentially expressed transcripts (DETs) were predicted, corresponding to 6936 up-regulated in AZ and 708 down-regulated in NAZ. We identified 489 candidate genes related to transcription factor, cell wall hydrolysis or modification, hydrolase activity, phytohormone signaling and response, lignin biosynthesis, and signal transduction or protein turnover. Eleven similar candidate genes involved in polygalacturonase activity, hydrolase activity, and mitogen-activated protein kinase were up-regulated in the abscission zone of the two Elymus species, suggesting these genes may have specific function for abscission zone development and seed shattering. A total of 67 polymorphic EST-SSR markers were developed and characterized based on the sequences of these candidate genes. Fourteen polymorphic EST-SSR primers were finally used to study genetic diversity in 48 E. nutans genotypes with contrasting seed shattering habit. The dendrogram based on molecular data showed that most accessions with similar seed shattering degree tended to group together. The expression data generated from this study provides an important resource for future molecular biological research. Many DETs were associated with abscission zone development, and EST-SSR loci related to candidate genes may have potential application in identifying trait-associated markers in E. nutans in the future.