Phylogeographic analysis reveals extensive genetic variation of native grass Elymus nutans (Poaceae) on the Qinghai-Tibetan plateau.

Introduction: Elymus nutans holds ecological and pastoral significance due to its adaptability and nutritional value, the Qinghai-Tibet Plateau (QTP) is a key hub for its genetic diversity. To conserve and harness its genetic resources in highland ecosystems, a thorough assessment is vital. However, a comprehensive phylogeographic exploration of E. nutans is lacking. The objective of this study was to unravel the genetic diversity, adaptation, and phylogenetics of E. nutans populations. Methods: Encompassing 361 individuals across 35 populations, the species' genetic landscape and dynamic responses to diverse environments were decoded by using four chloroplast DNA (cpDNA) sequences and nine microsatellite markers derived from the transcriptome. Results and discussion: This study unveiled a notable degree of genetic diversity in E. nutans populations at nuclear (I = 0.46, He = 0.32) and plastid DNA levels (Hd = 0.805, π = 0.67). Analysis via AMOVA highlighted genetic variation predominantly within populations. Despite limited isolation by distance (IBD), the Mekong-Salween Divide (MSD) emerged as a significant factor influencing genetic differentiation and conserving diversity. Furthermore, correlations were established between external environmental factors and effective alleles of three EST-SSRs (EN5, EN57 and EN80), potentially linked to glutathione S-transferases T1 or hypothetical proteins, affecting adaptation. This study deepens the understanding of the intricate relationship between genetic diversity, adaptation, and environmental factors within E. nutans populations on the QTP. The findings shed light on the species' evolutionary responses to diverse ecological conditions and contribute to a broader comprehension of plant adaptation mechanisms.

Genome constitution and evolution of Elymus atratus (Poaceae: Triticeae) inferred from cytogenetic and phylogenetic analysis.

BACKGROUND: Elymus atratus (Nevski) Hand.-Mazz. is perennial hexaploid wheatgrass. It was assigned to the genus Elymus L. sensu stricto based on morphological characters. Its genome constitution has not been disentangled yet. OBJECTIVE: To identify the genome constitution and origin of E. atratus. METHODS: In this study, genomic in situ hybridization and fluorescence in situ hybridization, and phylogenetic analysis based on the Acc1, DMC1 and matK sequences were performed. RESULTS: Genomic in situ hybridization and fluorescence in situ hybridization results reveal that E. atratus 2n = 6x = 42 is composed of 14 St genome chromosomes, 14 H genome chromosomes, and 14 Y genome chromosomes including two H-Y type translocation chromosomes, suggesting that the genome formula of E. atratus is StStYYHH. The phylogenetic analysis based on Acc1 and DMC1 sequences not only shows that the Y genome originated in a separate diploid, but also suggests that Pseudoroegneria (St), Hordeum (H), and a diploid species with Y genome were the potential donors of E. atratus. Data from chloroplast DNA showed that the maternal donor of E. atratus contains the St genome. CONCLUSION: Elymus atratus is an allohexaploid species with StYH genome, which may have originated through the hybridization between an allotetraploid Roegneria (StY) species as the maternal donor and a diploid Hordeum (H) species as the paternal donor.

Effects of Climatic Change on Phylogeography and Ecological Niche of the Endemic Herb Elymus breviaristatus on the Qinghai-Tibet Plateau.

Past climatic and topographic variations have created strong biogeographic barriers for alpine species and are key drivers of the distribution of genetic variation and population dynamics of species on the Qinghai-Tibet Plateau (QTP). Therefore, to better conserve and use germplasm resources, it is crucial to understand the distribution and differentiation of genetic variation within species. Elymus breviaristatus, an ecologically important rare grass species with strong resistance, is restricted to a limited area of the QTP. In this study, we investigated the phylogeography of E. breviaristatus using five chloroplast genes and spacer regions in natural populations distributed along the eastern QTP. We identified a total of 25 haplotypes among 216 individuals from 18 E. breviaristatus populations, which were further classified into four haplogroups based on geographical distribution and haplotype network analysis. Notably, we did not observe any signs of population expansion. High genetic diversity was exhibited at both species and population levels, with precipitation being the main limiting factor for population genetic diversity levels. Higher genetic diversity was exhibited by populations located near the Mekong-Salween Divide genetic barrier, suggesting that they may have served as a glacial refuge. The significant pattern of genetic differentiation by environmental isolation highlights the influence of heterogeneous environments on the genetic structure of E. breviaristatus populations. Additionally, the results of ecological niche models indicated that the geographic distribution of E. breviaristatus populations has decreased rapidly since the Last Glacial Maximum but is not threatened by future global warming.

Contribution of the bacterial community of poorly fermented oat silage to biogas emissions on the Qinghai Tibetan Plateau.

To better utilize poorly fermented oat silage on the Qinghai Tibetan Plateau, 239 samples of this biomass were collected from the plateau temperate zone (PTZ), plateau subboreal zone (PSBZ), and nonplateau climatic zone (NPCZ) in the region and analyzed for microbial community, chemical composition and in vitro gas production. Climatic factors affect the bacterial α-diversity and ß-diversity of poorly fermented oat silage, which led to the NPCZ having the highest relative abundance of Lactiplantibacillus plantarum. Furthermore, the gas production analysis showed that the NPCZ had the highest maximum cumulative gas emissions of methane. Through structural equation modeling analysis, environmental factors (solar radiation) affected methane emissions via the regulation of lactate production by L. plantarum. The enrichment of L. plantarum contributes to lactic acid production and thereby enhances methane emission from poorly fermented oat silage. Notably, there are many lactic acid bacteria detrimental to methane production in the PTZ. This knowledge will be helpful in revealing the mechanisms of environmental factors and microbial relationships influencing the metabolic processes of methane production, thereby providing a reference for the clean utilization of other poorly fermented silage.

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.

Phylogeographical Pattern and Population Evolution History of Indigenous Elymus sibiricus L. on Qinghai-Tibetan Plateau.

Elymus sibiricus L. is a perennial allotetraploid belonging to Triticeae of Poaceae, Elymus L., as the type species of genus Elymus L. The existing geographical distribution pattern and genetic spatial structure of E. sibiricus on Qinghai-Tibetan Plateau (QTP) are not yet clear. In this study, population genetic structure and demography history of 216 individuals from 44 E. sibiricus populations on QTP were studied used specific-locus amplified fragment sequencing (SLAF-seq). The result of genetic diversity showed that there was no single genetic diversity center was observed across all E. sibiricus populations. The results of genetic variation showed that 44 populations were clearly divided into the following three groups: Qinghai Plateau (Group I), South Tibet (Group II), and Hengduan Mountains (Group III). From the three analyses of AMOVA, Mantel test and Treemix, strong genetic differentiation across all populations and low genetic differentiation among populations within three groups. Molecular dating indicated that E. sibiricus diverged at 16.08 Ma (during the early Miocene) can be linked to the Himalayan Motion stage of QTP uplift. It is speculated that the reasons affecting the current phylogeographical pattern are as follows: (1) The environmental changes due to the uplift of the QTP; (2) The geographic distance between the populations (Groups I and III are close in geographic distance, and gene flow are frequent); (3) Geographical barriers (the Tanggula and Bayangela Mountains between Groups I and II). This study provides new evidence and historical perspective to the future exploration of the evolution and geographic distribution pattern of Elymus L.

Physiological and molecular responses of Phalaris arundinacea under salt stress on the Tibet plateau.

Phalaris arundinacea, with its characteristics of rapid growth and high biological yield, is regarded as an excellent forage grass in the Qinghai-Tibetan Plateau region of China. To explore the physiological and molecular response mechanism of Phalaris arundinacea under salt stress, we monitored the biomass and physiological indexes of two locally grown strains under conditions of exposure to 150 and 300 mM NaCl solution. Z0611 exhibited better salt stress tolerance than YS. Transcriptome sequencing analysis showed that YS and Z0611 had 1713 and 4290 differentially expressed genes (DEGs), respectively, including on metabolic processes, single-organism process, catalytic activity, and plant hormone signal transduction in the GO and KEGG databases. We also identified a large number of genes involved in hormone signaling, antioxidant systems, ion homeostasis, and photosynthetic systems. Our study provides physiological and molecular insight for establishing a salt resistance database and mining salt tolerance genes in Phalaris arundinacea, and also provides theoretical guidance for the restoration of saline-alkali land on the Qinghai-Tibet Plateau.

Ambient ultraviolet radiation: A new factor affecting anaerobic fermentation of oat and subsequent methane emissions.

To investigate the effect of ambient ultraviolet (UV) radiation on the Qinghai-Tibetan Plateau on anaerobic fermentation and subsequent methane emissions, oats grown under different ambient UV conditions (UV1, 100% ambient UV radiation; UV2, 50% ambient UV radiation) were fermented with molasses and Lactobacillus plantarum (LP) inoculant treatments. The UV2 treatment increased (P < 0.05) epiphytic bacterial counts and reduced (P < 0.05) the water-soluble carbohydrate contents of oat. Both the UV2 treatment and the LP inoculant substantially increased the abundance of Lactobacillus and inhibited detrimental microorganisms (enterobacteria and yeasts) during anaerobic fermentation, resulting in the absence of butyric acid. The UV2 treatment increased (P < 0.05) the acid detergent fiber contents after anaerobic fermentation and indirectly increased (P < 0.05) methane emissions by 8.52-14.69% during in vitro ruminal digestion. This study demonstrated that low ambient UV radiation during cultivation facilitated anaerobic fermentation and consequently enhanced subsequent methane emissions from the rumen.

Morphological, cytological, and molecular evidences for natural hybridization between Roegneria stricta and Roegneria turczaninovii (Triticeae: Poaceae).

Some plants with low fertility are morphologically intermediate between Roegneria stricta and Roegneria turczaninovii, and were suspected to be natural hybrids between these species. In this study, karyotype analysis showed that natural hybrids and their putative parents were tetraploids (2n = 4x = 28). Meiotic pairing in natural hybrids is more irregular than its putative parents. Results of genomic in situ hybridization and fluorescence in situ hybridization indicate that natural hybrids contain the same genome as their putative parents. The nuclear gene DNA meiotic recombinase 1 (DMC1) and the chloroplast gene rps16 of natural hybrids and their putative parents were analyzed for evidence of hybridization. The results from molecular data supported by morphology and cytology demonstrated that the plants represent natural hybrids between R. stricta and R. turczaninovii. The study is important for understanding species evolution in the genus since it demonstrates for the first time the existence of populations of natural homoploid hybrids in Roegneria. The study also reports for the first time that the composition of the genomic formula of R. turczaninovii is StY, confirming that the current taxonomic status is correct.

Confirmation of Natural Hybridization between Kengyilia (StStYYPP) and Campeiostachys (StStYYHH) (Triticeae: Poaceae) Based on Morphological and Molecular Cytogenetic Analyses.

Natural hybridization has been frequently observed in Triticeae; however, few studies have investigated the origin of natural intergeneric Triticeae hybrids. In the present study, we discovered three putative hybrid Triticeae plants in the Western Sichuan Plateau of China. Morphologically, the putative hybrids were intermediate between Kengyilia melanthera (2n = 6x = 42; StStYYPP) and Campeiostachys dahurica var. tangutorum (2n = 6x = 42; StStYYHH) with greater plant height and tiller number. Cytological analyses demonstrated that the hybrids were hexaploid with 42 chromosomes (2n = 6x = 42). At metaphase I, 12.10-12.58 bivalents and 13.81-14.18 univalents per cell were observed in the hybrid plants. Genomic in situ hybridization demonstrated that the hybrids had StStYYHP genomes. Phylogenetic analysis of Acc1 sequences indicated that the hybrids were closely related to K. melanthera and C. dahurica var. tangutorum. Our morphological, cytological, and molecular analyses indicate that these hexaploid natural hybrid plants may be hybrids of K. melanthera and C. dahurica var. tangutorum.

Physiological and transcriptional responses of Phalaris arundinacea under waterlogging conditions.

As a high-yielding forage grass, Phalaris arundinacea widely distributed in the Qinghai-Tibet Plateau region of China. To explore physiological and molecular response mechanism of Phalaris arundinacea under waterlogging, we analyzed the biomass and physiological indexes of three locally grown strains under the submerged condition of 10 cm. The material Z0611 showed the strongest waterlogging resistance while the YS showed the weakest performance. Transcriptome sequencing analysis demonstrated that the YS and Z0611 had 17010 and 7566 differently expression genes (DEGs), respectively, which were mainly concentrated in the metabolic process, cell, ribosome, phenylpropanoid biosynthesis pathway in GO and KEGG databases. We also identified a large number of genes involved in carbohydrate metabolism, hormone signaling regulation, transcription factors, antioxidant system, and ethylene signaling. Our research may provide a scientific basis for the restoration of wetland environment on the Qinghai-Tibet Plateau, and lay a foundation for further exploration of the waterlogging resistance genes of Phalaris arundinacea and breeding of new strains resistant with waterlogging stress.

Transcriptome profiles identify the common responsive genes to drought stress in two Elymus species.

Elymus, the largest genus of the Triticeae Dumort, is a forage grass in the Qinghai-Tibetan Plateau, where the climate has gradually become increasingly dry in recent years. To understand the mechanisms of the response to drought stress in Elymus species, we first investigated physiological and biochemical responses to polyethylene glycol (PEG-6000) simulated drought stress in two Elymus species, Elymus nutans and Elymus sibiricus, and found that E. nutans was more tolerant to drought stress than E. sibiricus. De novo transcriptome analysis of these two Elymus species treated with or without 10 % PEG-6000 revealed that a total of 1695 unigenes were commonly regulated by drought treatment in these two Elymus species, with 1614 unigenes up-regulated and 81 unigenes down-regulated. The coexpressed differentially expressed genes (DEGs) were enriched in regulation of transcription and gene expression in the GO database. KEGG pathway analysis indicated plant hormone signaling transduction were mostly enriched in co-expressed DEGs. Furthermore, genes annotated in the plant hormone signaling transduction were screened from co-expressed DEGs, and found that abscisic acid plays the major role in the drought stress tolerance of Elymus. Meanwhile, transcription factors screened from co-expressed DEGs were mainly classified into the ERF subfamily and WRKY, DREB, and HSF family members. Our results provide further reference for studying the response mechanism and culturing highly tolerant grasses of the Elymus species under drought stress.

Eutrema giganteum (Brassicaceae), a new species from Sichuan, southwest China.

Eutrema giganteum (Brassicaceae), a new species from Hengduan Mountains in Sichuan Province, southwest China, is described, and its relationships to the closely related E. yunnanense is discussed based on morphological, cytological, and molecular data. It is similar morphologically to E. yunnanense but is readily distinguished by having robust (vs. slender), erect (vs. decumbent), and branched (vs. mostly simple), and rather tall stems (60-110 cm vs. 20-60 cm); curved (vs. straight), smooth (vs. torulose), and shorter fruit (5-8 mm vs. 8-15 mm); and fewer ovules per ovary (1-4 vs. 6-10). All examined individuals from different populations of E. giganteum clustered into a single clade sister to E. yunnanense in phylogenetic analyses using the combined nuclear ITS and plastid DNA datasets. Our cytological studies revealed that the chromosome number of E. giganteum is 2n = 44, with a genome size of 1160 (±8) Mb, while that of E. yunnanense is 2n = 28, with a genome size of 718 (±15) Mb. Multiple lines of evidence support the recognition of E. giganteum as a distinct species well differentiated from E. yunnanense.

[NIRS for Predicting Nutritive Value of Elymus sibiricus].

Siberian wildrye (Elymus sibiricus L.) is one of the predominant pasture species in Qinghai-Tibet plateau of China. It supplies a large number of fodders to domestic animals in spring and winter, and provides a large proportion of the summer and autumn grazing in these alpine regions. Crude protein (CP), acid detergent fiber (ADF), neutral detergent fiber (NDF) and in vitro dry matter digestibility (IVDMD) are the most important aspects of nutritive value of forages. A successful application of near infrared spectroscopy (NIRS) in combination with partial least square regression (PLSR) for the determination of four parameters (CP, ADF, NDF and IVDMD) of Siberian wildrye was developed. The standard errors of calibration (SEC, %DM) and Rcal(2) values (in parentheses) were 0.3299(0.9945), 0.7799(0.9499), 1.3430(0.9133), and 1.3762(0.9822) for CP, ADF, NDF and IVDMD equations, respectively. The standard errors of prediction (SEP, %DM) and Rval(2) values (in parentheses) were 0.3621(0.9938), 0.7878(0.9449), 1.3852(0.8907), and 1.4303(0.9790) for CP, ADF, NDF and IVDMD, respectively. A good correlation (r>0.9438) was found between results from NIRS and the traditional chemical method, and residual predictive deviation (RPD) ranged from 3.02 to 12.63. It was concluded that NIR spectroscopic technique associated with chemometrics is sufficiently sensitive to allow the accurate prediction of the concentrations of compo nents (CP, ADF and NDF) and IVDMD of Siberian wildrye.

[Potential Applicability of Fecal NIRs: A Review].

Near-infrared reflectance spectroscopy (NIRS) is an inexpensive, rapid, environment-friendly and non-invasive analytical technique that has been extensively applied in the analysis of the dietary attributes and the animal products. Acquisition of dietary attributes is essential for nutritional diagnoses to provide animals with reasonable diet. Traditionally, the calibration equations for the prediction of dietary attributes (e. g. crude protein) are developed from feed NIR spectra and the results of conventional chemical analysis (i. e. reference data). It is difficult to obtain the NIR spectra of forages consumed by grazing animals, so the method of this calibration is inappropriate for free-grazing herbivores. Feces, as the animal's metabolites, contain the information about both the animal's diet and the animal itself. Recently, Fecal-NIRS (F. NIRS) has been directly used to monitor diet information (botanical composition, chemical composition and digestibility), based on correlation between reference data and fecal NIR profile. Subsequently, some additional application (such as sex and species discrimination, reproductive and parasite status) of F. NIRS also is outlined. In the last, application of NIRS in animal manure is summarized. NIRS was shown to be an alternative to conventional wet chemical methods for analyzing some nutrient concentrations in animal manure rapidly. Overall, this paper proves that F. NIRS is a rapid and valid tool for the determination of the dietary attributes and of the physiological status of animal, although more efforts need to be done to improve the accuracy of the F. NIRS technique. Several researchers in English have reviewed the applications of F. NIRS. In China, however, there is a paucity of research and application regarding F. NIRS. We expect that this paper in Chinese will be helpful to the development of F. NIRS in China. At the same time, we propose NIRS as a simple and rapid analytical method for predicting the main chemical composition (dry matter, organ matter, total solid, volatile solid, total nitrogen, total Kjeldahl nitrogen and ammonium nitrogen) in animal manure.

Artificial warming facilitates growth but not survival of plateau frog (Rana kukunoris) tadpoles in presence of gape-limited predatory beetles.

BACKGROUND: Global warming has been frequently demonstrated to increase growth rate in larval amphibians that have considerable phenotypic plasticity; this may lead to an increase in larval survival because large larvae are less likely to be captured by gape-limited predators. This study is to test whether warming could improve tadpole growth and thereby enhance the tadpole survival in plateau frog Rana kukunoris. METHODOLOGY: We conducted an experiment involving growing tadpoles under two contrasting temperatures, i.e. ambient temperature vs. warming by 3.8°C, with and without their major predators--the gape-limited predaceous diving beetles Agabus sp. in eastern Tibetan Plateau, in a factorial arrangement. We recorded the survival and measured body fresh weight and morphological characteristics of the tadpoles. PRINCIPAL FINDINGS: Warming significantly increased body fresh weight in the presence of predators after three weeks of treatments. However, the predators imposed significant and similar effects on the survival of tadpoles under both ambient and elevated temperatures, with the effects mostly occurring in the first three weeks of the experiment. Changes in the body form, i.e. the greater whole length at a given fresh weight and the longer tail at a given body length, could have acted as mechanisms of defense and escape for the tadpoles. CONCLUSIONS/SIGNIFICANCE: Warming did not increase tadpole survival with or without presence of predators. Moreover, an increased growth rate (due to warming in the presence of predators) was not a major factor contributing to the tadpole survival. We postulate that even if warming increases the tadpole growth rate in the plateau frog, it does not necessarily improve their survival in the presence of gape-limited predators.

Molecular insights of genetic variation in Erianthus arundinaceus populations native to China.

BACKGROUND: E. arundinaceus (Retz.) Jeswiet is a warm-season, tall-growing perennial species native to much southern portion in China. The grass has been extensively used in sugarcane breeding and is recently targeted as a bioenergy feedstock crop. However, information on the genetic structure of the Chinese wild germplasm is limited. Knowledge of genetic variation within and among populations is essential for breeding new cultivars in the species. The major objective of this study was to quantify the magnitude of genetic variation among and within natural populations in China. METHODOLOGY/PRINCIPAL FINDINGS: In this experiment, we analyzed genetic variation of 164 individuals of 18 populations collected from natural habitats in six Chinese provinces using 20 sequence-related amplified polymorphism (SRAP) primer pairs generating 277 polymorphic bands. Among and within the populations, the percentage of polymorphic bands (PPB) was 80.00% and 27.07%, genetic diversity (HE ) was 0.245 and 0.099, effective number of alleles (NE ) was 1.350 and 1.170, and Shannon's information index (I) was 0.340 and 0.147, respectively. The populations were clustered into six groups exhibiting a high level of genetic differentiation, which was highly associated with geographic origins of respective germplasm populations, but was not significantly associated with geographic distances between the populations. CONCLUSIONS/SIGNIFICANCE: This is the first report indicating that large genetic variation exists in the Chinese E. arundinaceus germplasm based on the SRAP molecular marker analysis of native populations. The genetic structure of populations in the species has been substantially affected by geographic landforms and environments. The diverse collection will be highly valuable in genetic improvement in the species per se and likely in sugarcane.

Assessment of worldwide genetic diversity of Siberian wild rye (Elymus sibiricus L.) germplasm based on gliadin analysis.

E. sibiricus L., the type species of the genus Elymus, is a perennial, self-pollinating and allotetraploid grass indigenous to Northern Asia, which in some countries can be cultivated as an important forage grass. In the present study, eighty-six Elymus sibiricus accessions, mostly from different parts of Asia, were assayed by gliadin markers based on Acid Polyacrylamide Gel Electrophoresis to differentiate and explore their genetic relationships. The genetic similarity matrix was calculated by 47 polymorphic bands, which ranged from 0.108 to 0.952 with an average of 0.373. The total Shannon diversity index (H(o)) and the Simpson index (H(e)) was 0.460 and 0.302, respectively. Cluster analysis showed a clear demarcation between accessions from Qinghai-Tibetan Plateau, China and the others as separate groups. The clustering pattern was probably dependent on geographic origin and ecological adaptability of the accessions. The population structure analysis based on Shannon indices showed that the proportion of variance within and among the five geographic regions of the Northern Hemisphere was 55.9 and 44.1%, respectively, or 63.4 and 36.6% within and among six Chinese provinces. This distinct geographical divergence was perhaps depended on ecogeographical conditions such as climate difference and mountain distribution. The results of gladin analysis in this study are useful for the collection and preservation of E. sibiricus germplasm resources.