Siberian wild rye (Elymus sibiricus L.): Genetic diversity of germplasm determined using DNA fingerprinting and SCoT markers

High-Altitude Genetic Selection and Genome-Wide Association Analysis of Yield-Related Traits in Elymus sibiricus L. Using SLAF Sequencing

The genetic adaptations to harsh climatic conditions in high altitudes and genetic basis of important agronomic traits are poorly understood in Elymus sibiricus L. In this study, an association population of 210 genotypes was used for population structure, selective sweep analysis, and genome-wide association study (GWAS) based on 88,506 single nucleotide polymorphisms (SNPs). We found 965 alleles under the natural selection of high altitude, which included 7 hub genes involved in the response to UV, and flavonoid and anthocyanin biosynthetic process based on the protein-protein interaction (PPI) analysis. Using a mixed linear model (MLM), the GWAS test identified a total of 1,825 significant loci associated with 12 agronomic traits. Based on the gene expression data of two wheat cultivars and the PPI analysis, we finally identified 12 hub genes. Especially, in plant height traits, the top hub gene (TOPLESS protein) encoding auxins and jasmonic acid signaling pathway, shoot apical meristem specification, and xylem and phloem pattern formation was highly overexpressed. These genes might play essential roles in controlling the growth and development of E. sibiricus. Therefore, this study provides fundamental insights relevant to hub genes and will benefit molecular breeding and improvement in E. sibiricus and other Elymus species.

Identification of Flowering Regulatory Networks and Hub Genes Expressed in the Leaves of Elymus sibiricus L. Using Comparative Transcriptome Analysis

Flowering is a significant stage from vegetative growth to reproductive growth in higher plants, which impacts the biomass and seed yield. To reveal the flowering time variations and identify the flowering regulatory networks and hub genes in Elymus sibiricus, we measured the booting, heading, and flowering times of 66 E. sibiricus accessions. The booting, heading, and flowering times varied from 136 to 188, 142 to 194, and 148 to 201 days, respectively. The difference in flowering time between the earliest- and the last-flowering accessions was 53 days. Furthermore, transcriptome analyses were performed at the three developmental stages of six accessions with contrasting flowering times. A total of 3,526 differentially expressed genes (DEGs) were predicted and 72 candidate genes were identified, including transcription factors, known flowering genes, and plant hormone-related genes. Among them, four candidate genes (LATE, GA2OX6, FAR3, and MFT1) were significantly upregulated in late-flowering accessions. LIMYB, PEX19, GWD3, BOR7, PMEI28, LRR, and AIRP2 were identified as hub genes in the turquoise and blue modules which were related to the development time of flowering by weighted gene co-expression network analysis (WGCNA). A single-nucleotide polymorphism (SNP) of LIMYB found by multiple sequence alignment may cause late flowering. The expression pattern of flowering candidate genes was verified in eight flowering promoters (CRY, COL, FPF1, Hd3, GID1, FLK, VIN3, and FPA) and four flowering suppressors (CCA1, ELF3, Ghd7, and COL4) under drought and salt stress by qRT-PCR. The results suggested that drought and salt stress activated the flowering regulation pathways to some extent. The findings of the present study lay a foundation for the functional verification of flowering genes and breeding of new varieties of early- and late-flowering E. sibiricus.

Development of EST-SSR Markers Linked to Flowering Candidate Genes in Elymus sibiricus L. Based on RNA Sequencing

Elymus sibiricus L. is an important cold-season grass with excellent cold and drought tolerance, good palatability, and nutrition. Flowering time is a key trait that affects forage and seed yield. Development of EST-SSR (expressed sequence tag simple sequence repeat) markers based on flowering genes contributes to the improvement of flowering traits. In the study, we detected 155 candidate genes related to flowering traits from 10,591 unigenes via transcriptome sequencing in early- and late-flowering genotypes. These candidate genes were mainly involved in the photoperiodic pathway, vernalization pathway, central integrator, and gibberellin pathway. A total of 125 candidate gene-based EST-SSRs were developed. Further, 15 polymorphic EST-SSRs closely associated to 13 candidate genes were used for genetic diversity and population structure analysis among 20 E. sibiricus accessions, including two contrasting panels (early-flowering and late-flowering). Among them, primer 28366, designed from heading date 3a (HD3a), effectively distinguished early- and late-flowering genotypes using a specifically amplified band of 175 bp. The polymorphic information content (PIC) value ranged from 0.12 to 0.48, with an average of 0.25. The unweighted pair group method analysis (UPGMA) cluster and structure analysis showed that the 20 E. sibiricus genotypes with similar flowering times tended to group together. These newly developed EST-SSR markers have the potential to be used for molecular markers assisted selection and germplasm evaluation of flowering traits in E. sibiricus.

Selection of Suitable Reference Genes for RT-qPCR Gene Expression Analysis in Siberian Wild Rye (Elymus sibiricus) under Different Experimental Conditions

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.

Genetic diversity and structure of Elymus tangutorum accessions from western China as unraveled by AFLP markers

BACKGROUND

Understanding genetic diversity of wild plant germplasm and the relationships between ecogeographic and genetic characteristics may provide insights for better utilizing and conserving genetic resources. Elymus tangutorum (Nevski) Hand.-Mazz, a cool-season hexaploid perennial, is an important pasture bunchgrass species used for forages and grassland restoration in Qinghai-Tibet Plateau and northwest China. In this study, 27 E. tangutorum accessions from diverse origins of western China were evaluated using AFLP markers in an effort to delve into the genetic relationships among them. The effects of eco-environmental factors and geographical isolation on the genetic diversity and population structure were also elucidated.

RESULTS

On account of 554 polymorphic fragments amplified with 14 primer combinations, the mean values of some marker parameters including polymorphic information content, resolving power and marker index were 0.2504, 14.10 and 23.07, respectively, validating the high efficiency and reliability of the markers selected. Genetic dissimilarity index values among accessions ranged from 0.1024 to 0.7137 with a mean of 0.2773. STRUCTURE, UPGMA clustering and PCoA analyses showed that all accessions could be divided into the three main clusters; however, this results do not exactly coincide with geographic groups. We found medium differentiation (F ST  = 0.162) between Qinghai-Tibet Plateau (QTP) and Xinjiang (XJC), and high differentiation (F ST  = 0.188) among three Bayesian subgroups. A significant correlation (r = 0.312) between genetic and geographical distance was observed by Mantel test in the species level, while the weak correlation was detected between genetic and environmental distance for all accessions and most of geographical groups. In addition, a significant ecological influence of average annual precipitation on genetic distance was revealed in XJC group and the Bayesian subgroup A.

CONCLUSION

This study indicates that AFLP technique are a useful tool to measure genetic diversity in E. tangutorum, showing that geographical and environmental factors (especially precipitation) together, play a crucial role in genetic differentiation patterns. These findings underline the importance of local adaptation in shaping patterns of genetic variability and population structure in E. tangutorum germplasm collected in Western China.

Hybrid identification and genetic variation of Elymus sibiricus hybrid populations using EST-SSR markers

Background

Elymus sibiricus is an important native grass in Qinghai-Tibetan Plateau. Seed shattering is a serious problem for E. sibiricus, especially at harvest time. Cross breeding is an effective way to create new varieties with beneficial characteristic or improved traits, and to broaden genetic base.

Results

In this study, we created five hybrid populations by crossing seven E. sibiricus genotypes that have seed shattering variation. Then, nine EST-SSR primers were used for hybrid identification based on DNA fingerprinting, and genetic diversity analysis of hybrid populations and their respective parents. A total of 15 hybrids were identified. An analysis of amplified polymorphic bands among genuine hybrids and their respective parents revealed mainly two types of markers: 1) hybrids shared bands exclusively amplified in both parents; 2)hybrids shared bands exclusively amplified in male parents. For each hybrid population, the total number of amplified bands ranged from 37 to 57, the percentage of polymorphism varied from 65.12% to 75.68%, with an average of 70.51%. Novel bands found in each hybrid population varied from 0 to 9 bands, suggesting an occurrence of rearrangements in the hybrid population. The structure analysis revealed that all hybrid populations and parents were assigned to eight groups. The principal coordinate analysis (PCoA) showed similar results.

Conclusions

In general, this study proved EST-SSR markers are efficient for hybrid identification, and suggested more genetic variation could be captured in hybrid populations by crossing breeding.