Next Article in Journal
Impact of Coating of Urea with Bacillus-Augmented Zinc Oxide on Wheat Grown under Salinity Stress
Previous Article in Journal
The Challenge for Rapid Detection of High-Structured Circular RNA: Assay of Potato Spindle Tuber Viroid Based on Recombinase Polymerase Amplification and Lateral Flow Tests
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 9
Issue 10
10.3390/plants9101371
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Yuying Zheng
,
Zongyu Zhang
,
Yiyang Wan
,
Jiaoyang Tian
and
Wengang Xie
*
The State Key Laboratory of Grassland Agro-ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou 730020, China
*
Author to whom correspondence should be addressed.
Plants 2020, 9(10), 1371; https://doi.org/10.3390/plants9101371
Submission received: 30 August 2020 / Revised: 5 October 2020 / Accepted: 14 October 2020 / Published: 15 October 2020
Browse Figures
Versions Notes

Abstract

:
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.

1. Introduction

Elymus sibiricus L. is a perennial, predominantly self-pollinating and allotetraploid (2n = 4x = 28) forage grass of the genus Elymus. It is widely distributed in high altitude regions of western and northern China, especially on the Qinghai-Tibetan Plateau [1]. E. sibiricus is an important cold-season grass with excellent cold and drought tolerance, good palatability, and forage quality [2]. However, research on breeding and improvement for this species is limited. The quantity of E. sibiricus varieties cannot meet the growing demands of ecological restoration and animal husbandry in China [3]. In perennial forage grasses, flowering indicates the transition from vegetative growth to reproductive growth, which affects characters of economic importance such as forage quality and seed yield [4]. Therefore, accurate prediction of flowering is very important for breeding early flowering or late flowering varieties to meet different practical needs [5].
To date, molecular mechanisms for flowering in Arabidopsis thaliana [6] and many other species have been extensively studied. The main pathways of controlling flowering in Arabidopsis thaliana are the vernalization pathway, photoperiodic pathway [7], gibberellin pathway [8], and autonomic pathway [9]. Several key flowering genes such as, Flowering Locus T (FT), Suppressor of Overexpression of Constans 1(SOC1), Leafy (LFY), and VERNALIZATION 1 (VRN1) have been reported [10]. FT gene, located at the intersection of different flowering pathways, plays a key role in plant flowering regulation [11]. SOC1 encoding a MADS-box transcription factor is highly expressed in leaves and floral meristems [12]. LFY gene, an important flowering control factor, integrates external environment and endogenous information, regulates the transformation of inflorescence to flowering meristem, and controls the flowering time under the joint action of multiple flowering pathways [13]. VRN1, a MADS-box gene, plays an important role in the development of spikelet and spike, and also affects the flowering stage and plant height in wheat [14]. In addition, Fang et al. [15] investigated the characterization and function of the rice OsFTL gene. Overexpression of OsFTL10 led to early flowering of rice by 2 weeks. Liu et al. [16] reported TaZIM-A1, which is an atypical GATa-like transcription factor that negatively regulated flowering in wheat. Under long-day conditions, overexpression of Tazima-A1 led to late flowering in wheat. Kang et al. [17] identified MsFTa, an FT ortholog in alfalfa, and characterized its function in regulating flowering. MsFTa plays a promoting role in alfalfa flowering, suggesting that it may participate in the coding of florigen. Huang et al. [18] revealed the genetic network regulating flowering time, and identified four major flowering candidate genes in orchardgrass based on whole genome sequencing. However, studies on the molecular mechanisms controlling flowering have not been reported in E. sibiricus.
As a high-throughput sequencing technique, transcriptome sequencing provides almost complete transcript information of biological tissues or cells [19]. Due to the advantages of sequencing depth, RNA-seq can reveal the overall gene expression of an individual at a given moment and in a given organization more comprehensively [20]. RNA-seq is widely used in non-model plants with limited sequence information, because it focuses on the coding region [21]. With the rapid development of molecular markers, EST-SSR (expressed sequence tag simple sequence repeat) markers derived from RNA-seq can be used to analyze genetic diversity, genetic variation, population structure, and functional studies [22]. Development of EST-SSR markers based on interesting candidate genes related to target traits are more valuable for variety improvement. Singh et al. [23] developed 30 SSR markers based on salt tolerance candidate genes in wheat, and evaluated genetic diversity in salt susceptible and salt tolerant wheat genotypes. Babu et al. [24] developed 24 candidate gene-based SSR markers related to lysine and tryptophan metabolic pathways in maize. These candidate gene-based SSRs is useful in genetic improvement of crop plants. Jespersen et al. [25] identified 29 SSRs related to five summer heat and drought tolerance candidate genes. These candidate gene-based SSR markers were useful in molecular marker-assisted selection of stress tolerance in creeping bentgrass. Lai et al. [26] developed EST-SSR derived from flowering genes and found these markers were highly transferable in 17 Gossypium species. Zhao et al. [27] reported EST-SSR loci related to candidate genes for seed shattering may have potential application in identifying trait-associated markers in E. nutans.
To better understand the mechanism for flowering in E. sibiricus, we developed novel EST-SSR markers derived from putative candidate genes for flowering time based on previous transcriptome data in E. sibiricus, and used polymorphic EST-SSRs for genetic diversity analysis of E. sibiricus genotypes with different flowering time. Early-flowering and late-flowering germplasms can be distinguished by polymorphic EST-SSRs, which could shorten the breeding process in E. sibiricus. These newly developed EST-SSR markers could be a useful tool for genetic improvement and molecular markers-assisted selection of flowering traits in E. sibiricus.

2. Results

2.1. Transcriptome Sequencing, Unigenes Annotation, and Differentially Expressed Genes Analysis

To dissect the molecular mechanism and explore the putative genes related to flowering time in E. sibiricus, 18 cDNA libraries were constructed from leaves collected from early flowering genotypes and late flowering genotypes at booting stage, heading stage, and flowering stage. Then, these libraries were sequenced. These Illumina data are available in the Sequence Reads Archive (SRA) with accession number SRX9217303-SRX9217319. A total of 400,453,230 clean reads were obtained after filtering out low-quality reads that had a Q30 percentage greater than 91.56%. The average GC content of clean reads was 55% (Table S1). Then, a total of 477,220 transcripts and 184,735 unigenes were identified in all samples. The average length of transcripts and unigenes was 937.56 bp and 750.13 bp, with an N50 length of 1420 bp and 1195 bp, respectively (Table S2).
After comparing BLAST and HMMER to 184,735 unigenes with COG, GO, KEGG, KOG, Pfam, Swiss-Prot, TrEMBL, Nr, and Nt databases, 166,071 (89.90%) unigenes were successfully annotated. The number and proportion of unigenes corresponding to each database are 15,597 (8.44%), 48,508 (26.26%), 12,379 (6.70%), 37,405 (20.25%), 41,919 (22.69%), 42,321 (22.91%), 90,201 (48.83%), 83,698 (45.31%), and 161,713 (87.54%) (Table S3).
To investigate the changes in gene expression and explore the putative genes involved in flowering. The differentially expressed genes were identified between the early- and late-flowering genotypes at three development stages. A total of 3526 differentially expressed genes (DEGs) were identified. Among these DEGs, 1715 genes were up-regulated in early-flowering genotypes and 1811 genes were down-regulated in late-flowering genotypes. In the booting stage, 576 genes were up-regulated in early-flowering genotypes and 360 genes down-regulated in late flowering genotypes. In the heading stage, 584 genes up-regulated were up-regulated in early-flowering genotypes and 578 genes down-regulated in late flowering genotypes. In the flowering stage, 555 genes up-regulated were up-regulated in early-flowering genotypes and 873 genes down-regulated in late flowering genotypes. By comprehensive comparison of three DEG sets, the largest number of DEGs was found in the flowering stage (Figure 1).

2.2. Frequency and Distribution of EST-SSR Markers

Among the transcriptome assembly, 13,052 SSRs were identified from 40,639 unigenes and the number of SSR containing sequences was 10,286. The number of sequences containing more than 1 SSR was 2176 and the number of SSRs present in compound formation were 742. The identified SSRs comprised of mono- to hexa-nucleotide repeats with a big difference in quantity (Table 1). Mono-nucleotide repeats were the most abundant (5471, 41.92%), followed by tri-nucleotide (4551, 34.87%), di-nucleotide (2642, 20.24%), tetra-nucleotide (314, 2.41%), penta-nucleotide (48, 0.37%), and hexa-nucleotide repeats (26, 0.20%) (Table 1). Among the SSR motifs, the most dominant mono-nucleotide repeat motif was A/T (4759, 36.46%), followed by tri-nucleotide repeat motifs of CCG/CGG (1638, 12.55%) and AGG/CCT (779, 5.97%) (Figure 2).

2.3. Development of Candidate Gene-Based EST-SSR Markers

Based on gene function annotation, we selected 155 candidate genes associated with flowering time from 10,591 E. sibiricus unigenes (Table S4). These candidate genes involved in vernalization response (6), flowering (31), circadian clock (12), photoperiodic response (4), CONSTANS-Like (45), gibberellin response (17), autonomous response (16), age pathways (5), and central integrator (19). Photoperiodic pathways were the most abundant (61, 39.35%), followed by vernalization pathways (37, 23.87%), central integrator (19, 12.26%), gibberellin pathways (17, 10.97%), autonomic pathways (16, 10.32%), and age pathways (5, 3.23%) (Figure 3). A total of 125 pairs of candidate gene-based SSR primers were designed and used to amplify four E. sibiricus genotypes with different origins. The unigene sequences used for marker development and primer sequences can be found in GenBank (BankIt2383886: MW 016765–MW 016888, and MW052548).
Among the 125 designed candidate gene-based EST-SSR (expressed sequence tag simple sequence repeat) primers, 106 pairs of primers generated clear bands with expected size (Table S5). Finally, 15 polymorphic primers were used to analyze genetic diversity among 20 E. sibiricus accessions. Among the 20 E. sibiricus accessions, the flowering time varied from 150 (Z8) to 193 days (W1), with an average of 169 days. The average flowering time for early flowering accessions and late flowering accessions were 154 days and 184 days, respectively (Figure 4).
Based on our results, primer 28366 designed from HD3a can distinguish early-flowering and late-flowering accessions directly by specifically-amplified bands. The 175 bp band was found in all late-flowering genotypes, but this band was absent in most early-flowering genotypes, except for Z3 and Z4 (Figure 5a). According to the results of PCR (polymerase chain reaction) product sequencing, the Z3 and Z4 genotypes have similar sequence with late-flowering genotypes (Figure 5b) except for some SNP mutations.

2.4. Validation of Candidate Gene-Based Specific Primers Authenticity

To determine the authenticity of candidate gene-based EST-SSR primers, amplicons from 20 E. sibiricus germplasms for four EST-SSRs were sequenced. The sequenced alleles from 20 E. sibiricus germplasms were homologous to the original SSR motifs from which the primer was designed. According to the sequencing results of expected bands generated from primer 42789, six TTG repeats and two TTG repeats were found in different flowering time genotypes. For primer 33508, E. sibiricus accessions with different flowering time had 8 T repeats, 10 repeats, and 13 T repeats, respectively. Marker polymorphisms among the 20 early-flowering and late-flowering E. sibiricus accessions were found due to variation in number of repeats of SSR motifs. Alignment of sequences obtained from selected PCR bands amplified by two primers (32835 and 41988) revealed that the base mutation occurs in different E. sibiricus accessions. The homology of coding sequences was further verified by multiple sequence alignment (Figure 6).

2.5. Genetic Diversity Analysis Using Candidate Gene-Based EST-SSR Markers

All 15 polymorphic primers generated 60 alleles, ranging from 1 to 9 per SSR locus, with an average of 4 alleles. The highest number of alleles was observed in primer 34500. Polymorphic information content (PIC) values of the 15 primers ranged from 0.12 to 0.48, with an average of 0.25. The highest PIC was observed in primer 46865 (0.4767), the lowest PIC was observed in primer 33680 (0.1275). Heterozygosity (Ho) ranged from 0.05 to 1.00 with an average of 0.77. Highest Ho was observed in primer 48198, 34500, 36927, 34261, and 28366 with a mean of 1. The highest gene diversity (He) was in primer 34500 (0.8731), ranging from 0.22 to 0.88 with an average of 0.61 (Table 2).
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). The dendrogram generated through UPGMA (unweighted pair group method analysis) analysis grouped all the 20 E. sibiricus accessions into three clusters. Cluster analysis results showed that when the genetic similarity coefficient of 0.71 (Table S6), cluster one contained 8 early-flowering accessions and 2 late-flowering accessions (PI435088 and PI595149). Cluster two contained 8 late-flowering accessions and one early genotype (W630476). Cluster three contained only one early-flowering accession (W610305) (Figure 7). In general, E. sibiricus genotypes with similar flowering time tended to group together.
We employed STRUCTURE-v2.3.2 software (Stanford University, CA, USA) to analyze the population structure among 20 E. sibiricus accessions with different values for K (from 1 to 11). Generally, accessions with similar flowering time tended to group together, which was similar to cluster analysis results. When K = 3, the optimal number of groups was three based on maximum likelihood and delta K (ΔK) values. Group 1 in red contained three late-flowering accessions and one early-flowering accession (W610305). Group 2 in green contained five late-flowering accessions and one early-flowering accession (W630476). Group 3 in blue contained eight early-flowering accessions and two late-flowering accessions (PI435088 and PI595149). The flowering candidate gene-based EST-SSR markers may have the potential to distinguish the E. sibiricus genotypes with different flowering time (Figure 8).

3. Discussion

3.1. Candidate Gene-Based EST-SSR Markers with Flowering Time for Marker-Assisted Selection

Flowering time is a turning point from vegetative growth to reproductive growth, which is one of the important factors determining plant productivity and adaptability to the environment in different latitudes [28]. Early-flowering germplasms have more time and resources for reproductive growth because of earlier flowering, so they have a higher seed yield. Late-flowering germplasms, have higher biomass, which is not conducive to seed production, because they need more time and resources for nutrient growth after flowering [29]. The development of early-flowering or late-flowering varieties depends on the accurate evaluation of flowering time. It is difficult to identify flowering time accurately by phenotypic trait due to environment factors. Compared with phenotypic identification, candidate gene-based DNA markers have great potential to improve the efficiency and accuracy of traditional plant breeding through marker-assisted selection (MAS) [30].
Now, the development of transcriptome sequencing technology provides an opportunity for plant genome analysis. With the potential of high throughput, high precision, and low cost, next-generation sequencing technology (NGS) has been widely used in the qualitative and quantitative analyses of transcriptome and has been successfully applied to many plant species [31]. SSRs are classified into genomic-SSR and EST-SSR based on the original sequences used to identify simple repeats [32]. Compared with genomic SSRs, EST-SSR markers have a higher rate of transferability owing to the regions being more evolutionarily conserved than non-coding sequences [33] and have a higher probability of being functionally associated with differential gene expression [34]. The EST-SSR markers based on candidate genes are more likely to be associated with important agronomically traits. E. sibiricus is the non-model herbage of Elymus genus that does not have whole genomic sequences. Hence, the development of EST-SSR markers based on candidate genes is more suitable for the study of flowering mechanism and improvement of desired trait in E. sibiricus.
Kumar et al. [35] identified 13 EST-SSRs significantly associated with flowering time in a diverse panel of 96 accessions of lentil germplasms, and used them as functional markers in the lentil breeding program to develop short duration cultivars. In addition, candidate gene-based EST-SSR explained large phenotypic variation (2.3–21.8%) compared to genomic SSR markers (2.1–10.2%). Li et al. [36] designed 51 SSR primers based on flowering candidate genes and successfully applied them to amplify in 14 Paphiopedilum orchids genotypes. These valuable EST-SSR markers were important and useful for novel gene discovery and marker-assisted studies in Paphiopedilum. In this study, 15 EST-SSRs closely associated with flowering candidate genes were used for germplasms identification and genetic diversity. Based on our results, most genotypes with different flowering time can be distinguished by EST-SSR markers, especially primer 28366, which can identify all 10 late-flowering genotypes and 8 early-flowering genotypes. Two genotypes, Z3 and Z4, were identified as early-flowering genotypes by flowering time measurement, but they were identified as late-flowering genotypes by specific amplified bands. This indicates that the reliability of candidate gene-based EST-SSR markers was higher than that of phenotypic observation. Primer 28366 was designed from HD3a gene related to the central integrator. The photoperiod sensitive gene HD3 was initially detected as a heading date-related quantitative trait locus on chromosome 6 of rice. The HD3 region contains two tightly connected loci, HD3a and HD3b. HD3a is a rice ortholog of the Arabidopsis FT gene and the Kasalath allele of HD3a promoted heading under short-day conditions [37,38]. Armstead et al. [39] reported a genetic co-spectral region between rice and ryegrass, including the HD3 heading-date QTL and a major QTL of rice, accounted for 70% variation. In addition, newly developed markers related the other flowering candidate genes in this experiment have been reported in other species, such as ELF3, GID1, GA2OX6, XCT, and GIGANTEA [40,41,42,43,44]. In future studies, we will explore the function and regulatory mechanism of these genes in E. sibiricus by studying the effect of nucleotide variation of these genes on flowering by sequencing of candidate gene-based markers.

3.2. Genetic Diversity of Candidate Gene-Based EST-SSR Markers

Different genetic parameters are important indicators of the origin, evolution, and distribution of new SSR markers. Moreover, they are also important means to measure SSR markers, indicating the existence of moderate heterozygosity in the tested germplasm [24]. In this study, repeated nucleotide types from mono- to hexa-nucleotide were detected in 13,052 identified SSRs from 40,639 unigenes, representing 32.12% of the total detected unigenes. The significant difference was among motif types and frequencies existing in EST-SSRs [45]. In this study, mono-nucleotide repeats were the most abundant SSR motifs (41.92%), followed by tri-nucleotide (34.87%), which was the same as in Melilotus [46]. However, di-nucleotide repeats are the most abundant SSR motifs in tall fescue [34] and Lolium multiflorum [47], tri-nucleotide repeats are the most abundant SSR motifs in E. sibiricus [48] and Medicago sativa [49].
Genetic diversity is a key determinant of plant germplasm allocation and genetic improvement. The polymorphic information content (PIC) ranged from 0.12 to 0.48, with an average of 0.25 for 20 E. sibiricus accessions, which was similar to E. nutans. The PIC varied from 0.220 to 0.370, with an average of 0.318 [27]. However, the value of PIC was lower than barely (the PIC ranged from 0.08 to 0.75 with a mean of 0.46) [50], common vetch (Vicia sativa subsp. sativa) (the PIC ranged from 0.09 to 0.98) [51], and E. sibiricus (the PIC ranged from 0.39 to 0.81) [48]. The reason for this result may due to that EST-SSR primers designed for flowering candidate genes were relatively conservative. A further genome-wide set of SSR markers may result in a higher polymorphism. The dendrogram results showed that these newly SSR markers had the potential to be used for distinguishing E. sibiricus germplasms with different flowering time. In total, 20 early-flowering and late-flowering accessions were not completely distinguished possibly due to eco-geographical factors such as, altitude, longitude, latitude, and climate conditions. Climate change at different altitudes will affect flowering time [52]. Grenier et al. [53] reported flowering time were highly correlated with latitudinal and racial distributions of landraces that can be affected by day-length variation. Stuerz et al. [54] found that the key climatic factors affecting flowering were humidity and temperature. Mean air temperature explained 81% variation in duration to flowering across sites, which was furthermore significantly influenced by relative air humidity in lowland rice. Burgarella et al. [55] reported flowering time influence from climatic conditions mainly by variation at genes located upstream in the flowering pathways, close to the environmental stimuli. Variables related to annual precipitation were better than other factors (e.g., temperature, altitude, latitude, or longitude) to reflect the constraint of the flowering time gene. In the study, Z3 and Z4 accessions were identified as late-flowering accessions by primer 28366 but these accessions were still grouped into early-flowering accessions. The result may be due to that flowering time is a complex quantitative trait, which is affected by polygene effect [56]. The monogenic expression of the HD3, Z3, and Z4 accessions were characterized as late-flowering genotypes. The regulation effect of single gene expression in the whole flowering network is a micro effect. Therefore, Z3 and Z4 accessions were grouped into early-flowering accessions under co-expression of 13 flowering candidate genes. At the same time, candidate gene-based EST-SSRs development in this study were not effective in identifying all the materials, which may be due to the number of selected genes or the effectiveness of the primers. Future, we can develop recombinant inbred lines (RILs) by crossing early-flowering accessions with the late-flowering accessions to detect flowering QTLs. The development of diagnostic markers could facilitate us to understand the genetic control of flowering in E. sibiricus.
The newly designed candidate gene-based SSR markers are useful in marker-assisted selection and the genetic improvement of the flowering trait in E. sibiricus. According to different production requirements, E. sibiricus germplasms with different flowering time can be selected accurately to meet different agricultural needs.

4. Materials and Methods

4.1. Plant Materials

In total, 20 E. sibiricus accessions including 10 early-flowering accessions and 10 late-flowering accessions were selected for candidate genes-based EST-SSR markers development and diversity analysis in this study (Table 3). All plant materials were planted in the experimental field at Lanzhou University, Yuzhong, Gansu, China (latitude 35°34′ N, longitude 103°34′ E, elevation 1720 m). The flowering time of each accession was determined as the number of days from 1 January to the date when the first inflorescence fully emerged from the flag leaf according to the methods described by Xie et al. [57]. Totally, the flowering time of 20 E. sibiricus accessions were measured. According to flowering time results, early-flowering accessions were recoded as Z1 through Z10 and late-flowering accessions were recoded as W1 through W10.

4.2. Development of Candidate Gene-Based EST-SSR Markers Based on RNA-Seq in E. sibiricus

In our previous study, six E. sibiricus accessions, including three early flowering accessions (PI598781, PI655199, LQ10) and three late flowering accessions (PI531665, PI531669, PI595169) were used for transcriptome analysis. Leaf tissues collected at booting stage, heading stage, and flowering stage were used for RNA extraction, cDNA library construction, and Illumina sequencing. These plants were grown in the field plots at the Yuzhong experimental field of Lanzhou University, Gansu, China (latitude 35°34′ N, longitude 103°34’ E, elevation 1720 m). The flag leaf tissues of each accession were collected at three development stage: booting stage, heading stage, and flowering stage.
A total of 18 samples were sent to Breeding Biotechnologies (Xi’an, Shaanxi, China). RNA was extracted with the UNIQ-10 column Trizol total RNA extraction kit (Sangon, Shanghai). cDNA libraries were constructed after the sample is qualified and then all 18 libraries were sequenced using Illumina HiSeq 2500 (data not published) [58]. Raw data were filtered to remove adaptor sequences and low-quality reads to obtain high quality clean data. Unigenes were obtained by de novo transcriptome assembly by Trinity-v2.4.0. Transcripts were annotated by performing BLAST-v2.2.28 against public databases, including Nr (the NCBI non-redundant protein database), Swiss-Prot (Annotated protein sequence database), GO (Gene Ontology), CGO (Clusters of Orthologous Groups), KOG (EuKaryotic Orthologous Groups), KEGG (Kyoto Encyclopedia of Genes and Genomes), and Protein Family (Pfam). We conducted pairwise comparison at three developmental stage between the two genotypes. The differentially expressed genes (DEGs) were filtered with expression levels FDR < 0.01, log2 fold change ≥ 2. The number of unigenes was screened for analysis of SSR using the MISA-v2.1 [59] set to identify mono-, di-, tri-, tetra-, penta-, and hexa-nucleotide motifs. Results of RNA-seq were validated via qRT-PCR experiments.
In this study, candidate genes were homologous genes of the reported flowering genes, which were selected from gene function annotation related to flowering pathways (i.e., vernalization pathway, photoperiodic pathway, gibberellin pathway, autonomic pathway, and age pathway), flowering key genes, and phytohormone responses. Flowering is polygenically controlled and more than 200 predictive genes and other undiscovered flowering genes are involved. Therefore, we expanded the search scope beyond the DEGs and searched the annotation results of the whole unigenes. The candidate gene-based EST-SSR markers were more closely related to flowering regulation than SSR markers. The candidate gene-based EST-SSR primers were designed using Primer 3-v2.3.4 [60] based on the MISA result followed by primer length from 15–30 bp, GC contents from 40–60%, annealing temperature of 60 °C, and no self-complementary sequence. All specific primers were synthesized by the GenScript Company (Nanjing, China).

4.3. DNA Extraction, Genotyping and Primer Validation

Fresh leaf tissue of each accession was sampled for DNA extraction using the SDS (sodium dodecyl sulfate) method [61]. DNA quality was determined using the NanoDrop ND1000 spectrophotometer (Thermo Scientific, Waltham, MA, USA). The qualified DNA was diluted to 25 ng/μL and stored at −20 °C. Four DNA samples (Z6, Z9, W2, and W6) with large geographical differences were selected for primers screening.
PCR amplifications were conducted in a volume of 15 μL containing 7.5 μL 2× reaction mix, 5.2 μL ddH2O, 0.2 μL golden DNA polymerase (Beijing Tiangen Corporation), 0.05 μL forward primer and reverse primer, and 2 μL of genomic DNA. The PCR conditions were preheat 2 min at 94 °C; 94 °C for 15 s, 10 cycles; annealing temperature 62 °C for 30 s; 72 °C for 30 s; 94 °C for 30 s, 25 cycles; annealing temperature 52 °C for 30 s; 72 °C for 1 min; 72 °C for 7min; and hold at 4 °C. Fragments were separated on 6.0% non-denaturing polyacrylamide gels (PAGE) and visualized using silver staining. In addition, the DL500 DNA marker was used to determine the sizes of the PCR products. Unpurified PCR products were sequenced by Shanghai Sangon Biological Engineering Technology (Shanghai, China) used to the primer validation aligned by MEGA7 software [62].

4.4. Allele Scoring and Polymorphism Detection

The polymorphic EST-SSRs closely associated to flowering candidate genes were used for genetic diversity and population structure analysis among 20 E. sibiricus accessions including two contrasting panel (early-flowering and late-flowering). The number of alleles of each candidate gene-based EST-SSR locus was calculated based on presence (1) or absence (0). The genetic relationship was analyzed based on a data matrix constructed by EST-SSR markers data using the NTSYSpc-v2.10e software (Exeter Software, Setauket, NY, USA). In addition, a dendrogram was constructed for early-flowering and late-flowering E. sibiricus genotypes following the unweighted pair group method analysis (UPGMA) based on Jaccard’s similarity coefficient [63]. Heterozygosity (Ho), gene diversity (He), and the polymorphic information content (PIC) were calculated as previously reported. PIC was calculated for each polymorphic primer using the formula: PIC = 1 − p2 − q2, where p is the frequency of present band, and q is the frequency of the absent band. Heterozygosity (Ho) is the ratio of heterozygous individuals in the population. At a single locus, the formula is as follows: Ho = 1 − ∑ i = 1 k P l u u . Gene diversity (He) is often referred to as expected heterozygosity. He = 1 − Σ p i 2 , where pi is frequency of the ith allele [64]. The number of heterozygotes is calculated by the direct count method. The model-based approach implemented was used to subdivide the individuals into different subgroups in the program STRUCTURE-v2.3.2 (Stanford University, CA, USA). The flowering time heatmap of early-flowering and late-flowering E. sibiricus germplasms made by Heatmap Illustrator-v1.0 (Huazhong University of Science and Technology, Wuhan, China).

5. Conclusions

In this study, 125 candidate gene-based EST-SSRs were developed. Further, 15 EST-SSRs closely associated to 13 candidate genes were used to analyze genetic diversity and population structure among 20 E. sibiricus accessions including two contrasting panels (early-flowering and late-flowering). The UPGMA cluster and structure analysis showed that the 20 E. sibiricus accessions with similar flowering times tended to group together. These newly developed EST-SSR markers could be useful for molecular markers-assisted selection and germplasm evaluation of the flowering traits in E. sibiricus.

Supplementary Materials

The following are available online at https://www.mdpi.com/2223-7747/9/10/1371/s1. Table S1: Summary of the sequence data analysis. Table S2: Statistics of Unigene Library of E. sibiricus. Table S3: BLAST analysis of the non-redundant unigenes against public databases. Table S4: The unigene sequence of 155 candidate genes. Table S5: Information of flowering candidate genes and gene-based EST-SSRs primers: reverse and forward sequence, Tm, molecular weight of expected band. Table S6: The coefficient of genetic similarity of 20 E. sibiricus accessions.

Author Contributions

W.X. conceived and designed the experiments. Y.W. and J.T. performed field management. Y.Z., Z.Z. Y.W. and J.T. performed the DNA extraction and EST-SSR markers development. Y.Z. and Z.Z. analyzed the data and performed data visualization. Y.Z. wrote the initial draft. W.X. and Y.Z. revised the manuscript. W.X. acquired financial support for the project leading to this publication. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Gansu Provincial Science and Technology Major Projects (19ZD2NA002), Chinese National Natural Science Foundation (31971751), The State’s Key Project of Research and Development Plan (2019YFC0507702), and the open projects of Key Laboratory of Superior Forage Germplasm in the Qinghai-Tibetan Plateau (2020-ZJ-03).

Acknowledgments

The authors are grateful to the help provided by the Changjiang Scholars and Innovation Research Team in University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xie, W.; Zhao, X.; Zhang, J.; Wang, Y.; Liu, W. Assessment of genetic diversity of Siberian wild rye (Elymus sibiricus L.) germplasms with variation of seed shattering and implication for future genetic improvement. Biochem. Syst. Ecol. 2015, 58, 211–218. [Google Scholar] [CrossRef]
  2. Ma, X.; Zhang, X.Q.; Zhou, Y.H.; Bai, S.Q.; Liu, W. Assessing genetic diversity of Elymus sibiricus (Poaceae: Triticeae) populations from Qinghai-Tibet Plateau by ISSR markers. Biochem. Syst. Ecol. 2008, 36, 514–522. [Google Scholar] [CrossRef]
  3. Xie, W.; Zhang, J.; Zhao, X.; Zhang, J.; Wang, Y. Siberian wild rye (Elymus sibiricus L.): Genetic diversity of germplasm determined using DNA fingerprinting and SCoT markers. Biochem. Syst. Ecol. 2015, 60, 186–192. [Google Scholar] [CrossRef]
  4. Craufurd, P.Q.; Wheeler, T.R. Climate change and the flowering time of annual crops. J. Exp. Bot. 2009, 60, 2529–2539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Torabi, B.; Adibnya, M.; Rahimi, A.; Azari, A. Modeling flowering response to temperature and photoperiod in safflower. Ind. Crop. Prod. 2020, 151, 10. [Google Scholar] [CrossRef]
  6. Simpson, G.G.; Dean, C. Flowering—Arabidopsis, the rosetta stone of flowering time? Science 2002, 296, 285–289. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Andrés, F.; Coupland, G. The genetic basis of flowering responses to seasonal cues. Nat. Rev. Genet. 2012, 13, 627–639. [Google Scholar] [CrossRef]
  8. Osnato, M.; Castillejo, C.; Matias-Hernandez, L.; Pelaz, S. TEMPRANILLO genes link photoperiod and gibberellin pathways to control flowering in Arabidopsis. Nat. Commun. 2012, 3, 808. [Google Scholar] [CrossRef] [Green Version]
  9. Michaels, S.D.; Amasino, R.M. Loss of FLOWERING LOCUS C activity eliminates the late-flowering phenotype of FRIGIDA and autonomous pathway mutations but not responsiveness to vernalization. Plant Cell 2001, 13, 935–941. [Google Scholar] [CrossRef] [Green Version]
  10. Levy, Y.Y.; Mesnage, S.; Mylne, J.S.; Gendall, A.R.; Dean, C. Multiple roles of Arabidopsis VRN1 in vernalization and flowering time control. Science 2002, 297, 243–246. [Google Scholar] [CrossRef]
  11. Corbesier, L.; Vincent, C.; Jang, S.H.; Fornara, F.; Fan, Q.Z.; Searle, I.; Giakountis, A.; Farrona, S.; Gissot, L.; Turnbull, C.; et al. FT protein movement contributes to long-distance signaling in floral induction of Arabidopsis. Science 2007, 316, 1030–1033. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Samach, A.; Onouchi, H.; Gold, S.E.; Ditta, G.S.; Schwarz-Sommer, Z.; Yanofsky, M.F.; Coupland, G. Distinct roles of CONSTANS target genes in reproductive development of Arabidopsis. Science 2000, 288, 1613–1616. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Yu, Y.C.; Qiao, L.F.; Chen, J.C.; Rong, Y.H.; Zhao, Y.H.; Cui, X.K.; Xu, J.P.; Hou, X.M.; Dong, C.H. Arabidopsis REM16 acts as a B3 domain transcription factor to promote flowering time via directly binding to the promoters of SOC1 and FT. Plant J. 2020. [Google Scholar] [CrossRef]
  14. Li, C.X.; Lin, H.Q.; Chen, A.; Lau, M.; Jernstedt, J.; Dubcovsky, J. Wheat VRN1, FUL2 and FUL3 play critical and redundant roles in spikelet development and spike determinacy. Development 2019, 146, 11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Fang, M.C.; Zhou, Z.J.; Zho, X.S.; Yang, H.Y.; Li, M.R.; Li, H.Q. Overexpression of OsFTL10 induces early flowering and improves drought tolerance in Oryza sativa L. PeerJ 2019, 7, 21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Liu, H.; Li, T.; Wang, Y.M.; Zheng, J.; Li, H.F.; Hao, C.Y.; Zhang, X.Y. TaZIM-A1 negatively regulates flowering time in common wheat (Triticum aestivum L.). J. Integr. Plant Biol. 2019, 61, 359–376. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Kang, J.M.; Zhang, T.J.; Guo, T.; Ding, W.; Long, R.C.; Yang, Q.C.; Wang, Z. Isolation and Functional Characterization of MsFTa, a FLOWERING LOCUS T Homolog from Alfalfa (Medicago sativa). Int. J. Mol. Sci. 2019, 20, 1968. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Huang, L.K.; Feng, G.Y.; Yan, H.D.; Zhang, Z.R.; Bushman, B.S.; Wang, J.P.; Bombarely, A.; Li, M.Z.; Yang, Z.F.; Nie, G.; et al. Genome assembly provides insights into the genome evolution and flowering regulation of orchardgrass. Plant Biotechnol. J. 2020, 18, 373–388. [Google Scholar] [CrossRef]
  19. Wang, Z.; Gerstein, M.; Snyder, M. RNA-Seq: A revolutionary tool for transcriptomics. Nat. Rev. Genet. 2009, 10, 57–63. [Google Scholar] [CrossRef]
  20. Wang, W.; Wang, Y.J.; Zhang, Q.; Qi, Y.; Guo, D.J. Global characterization of Artemisia annua glandular trichome transcriptome using 454 pyrosequencing. BMC Genom. 2009, 10, 10. [Google Scholar] [CrossRef] [Green Version]
  21. Der, J.P.; Barker, M.S.; Wickett, N.J.; de Pamphilis, C.W.; Wolf, P.G. De novo characterization of the gametophyte transcriptome in bracken fern, Pteridium aquilinum. BMC Genom. 2011, 12, 14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Arghavan, A.; Shekoufeh, E.; Sahar, A.; Mahsa, H.; Behrouz, S.; Hassan, J.; Hossein, F.; Sadegh, M.-F.; Fariba, R. Parallel consideration of SSRs and differentially expressed genes under abiotic stress for targeted development of functional markers in almond and related Prunus species. Sci. Hortic. 2015, 198, 462–472. [Google Scholar] [CrossRef]
  23. Singh, A.K.; Chaurasia, S.; Kumar, S.; Singh, R.; Kumari, J.; Yadav, M.C.; Singh, N.; Gaba, S.; Jacob, S.R. Identification, analysis and development of salt responsive candidate gene based SSR markers in wheat. BMC Plant Biol. 2018, 18, 15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Babu, B.K.; Agrawal, P.K.; Gupta, H.S.; Kumar, A.; Bhatt, J.C. Identification of candidate gene-based SSR markers for lysine and tryptophan metabolic pathways in maize (Zea mays). Plant Breed. 2012, 131, 20–27. [Google Scholar] [CrossRef]
  25. Jespersen, D.; Ma, X.Q.; Bonos, S.A.; Belanger, F.C.; Raymer, P.; Huang, B.R. Association of SSR and candidate gene markers with genetic variations in summer heat and drought performance for creeping bentgrass. Crop. Sci. 2018, 58, 2644–2656. [Google Scholar] [CrossRef]
  26. Lai, D.Y.; Li, H.Z.; Fan, S.L.; Song, M.Z.; Pang, C.Y.; Wei, H.L.; Liu, J.J.; Wu, D.; Gong, W.F.; Yu, S.X. Generation of ESTs for flowering gene discovery and SSR marker development in upland cotton. PLoS ONE 2011, 6, e28676. [Google Scholar] [CrossRef] [Green Version]
  27. Zhao, Y.Q.; Zhang, J.C.; Zhang, Z.Y.; Xie, W.G. Elymus nutans genes for seed shattering and candidate gene-derived EST-SSR markers for germplasm evaluation. BMC Plant Biol. 2019, 19. [Google Scholar] [CrossRef]
  28. Danilevskaya, O.N.; Meng, X.; McGonigle, B.; Muszynski, M.G. Beyond flowering time Pleiotropic function of the maize flowering hormone florigen. Plant Signal. Behav. 2011, 6, 1267–1270. [Google Scholar] [CrossRef] [Green Version]
  29. Bushman, B.S.; Robins, J.C.; Jensen, K.B. Dry matter yield, heading date, and plant mortality of Orchardgrass subspecies in a semiarid environment. Crop. Sci. 2012, 52, 745–751. [Google Scholar] [CrossRef]
  30. Collard, B.C.Y.; Mackill, D.J. Marker-assisted selection: An approach for precision plant breeding in the twenty-first century. Philos. Trans. R. Soc. B-Biol. Sci. 2008, 363, 557–572. [Google Scholar] [CrossRef] [Green Version]
  31. Unamba, C.I.; Nag, A.; Sharma, R.K. Next Generation Sequencing Technologies: The doorway to the unexplored genomics of non-model plants. Front. Plant Sci. 2015, 6, 1074. [Google Scholar] [CrossRef] [Green Version]
  32. Powell, W.; Machray, G.C.; Provan, J. Polymorphism revealed by simple sequence repeats. Trends Plant Sci. 1996, 1, 215–222. [Google Scholar] [CrossRef]
  33. Gupta, P.K.; Rustgi, S.; Sharma, S.; Singh, R.; Kumar, N.; Balyan, H.S. Transferable EST-SSR markers for the study of polymorphism and genetic diversity in bread wheat. Mol. Genet. Genom. 2003, 270, 315–323. [Google Scholar] [CrossRef]
  34. Saha, M.C.; Mian, M.A.R.; Eujayl, I.; Zwonitzer, J.C.; Wang, L.J.; May, G.D. Tall fescue EST-SSR markers with transferability across several grass species. Theor. Appl. Genet. 2004, 109, 783–791. [Google Scholar] [CrossRef] [PubMed]
  35. Kumar, J.; Gupta, S.; Biradar, R.S.; Gupta, P.; Dubey, S.; Singh, N.P. Association of functional markers with flowering time in lentil. J. Appl. Genet. 2018, 59, 9–21. [Google Scholar] [CrossRef]
  36. Li, D.M.; Wu, W.; Zhang, D.; Liu, X.R.; Liu, X.F.; Lin, Y.J. Floral transcriptome analyses of four Paphiopedilum Orchids with distinct flowering behaviors and development of simple sequence repeat markers. Plant Mol. Biol. Rep. 2015, 33, 1928–1952. [Google Scholar] [CrossRef]
  37. Monna, L.; Lin, H.X.; Kojima, S.; Sasaki, T.; Yano, M. Genetic dissection of a genomic region for a quantitative trait locus, Hd3, into two loci, Hd3a and Hd3b, controlling heading date in rice. Theor. Appl. Genet. 2002, 104, 772–778. [Google Scholar] [CrossRef]
  38. Kojima, S.; Takahashi, Y.; Kobayashi, Y.; Monna, L.; Sasaki, T.; Araki, T.; Yano, M. Hd3a, a rice ortholog of the Arabidopsis FT gene, promotes transition to flowering downstream of Hd1 under short-day conditions. Plant Cell Physiol. 2002, 43, 1096–1105. [Google Scholar] [CrossRef] [Green Version]
  39. Armstead, I.P.; Turner, L.B.; Farrell, M.; Skot, L.; Gomez, P.; Montoya, T.; Donnison, I.S.; King, I.P.; Humphreys, M.O. Synteny between a major heading-date QTL in perennial ryegrass (Lolium perenne L.) and the Hd3 heading-date locus in rice. Theor. Appl. Genet. 2004, 108, 822–828. [Google Scholar] [CrossRef] [PubMed]
  40. Anwer, M.U.; Davis, A.; Davis, S.J.; Quint, M. Photoperiod sensing of the circadian clock is controlled by EARLY FLOWERING 3 and GIGANTEA. Plant J. 2020, 101, 1397–1410. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Liu, Y.; Hao, X.Y.; Lu, Q.H.; Zhang, W.F.; Zhang, H.J.; Wang, L.; Yang, Y.J.; Xiao, B.; Wang, X.C. Genome-wide identification and expression analysis of flowering-related genes reveal putative floral induction and differentiation mechanisms in tea plant (Camellia sinensis). Genomics 2020, 112, 2318–2326. [Google Scholar] [CrossRef]
  42. Li, C.; Zheng, L.L.; Wang, X.N.; Hu, Z.B.; Zheng, Y.; Chen, Q.H.; Hao, X.C.; Xiao, X.; Wang, X.B.; Wang, G.D.; et al. Comprehensive expression analysis of Arabidopsis GA2-oxidase genes and their functional insights. Plant Sci. 2019, 285, 1–13. [Google Scholar] [CrossRef]
  43. Martin-Tryon, E.L.; Harmer, S.L. XAP5 CIRCADIAN TIMEKEEPER coordinates light signals for proper timing of photomorphogenesis and the circadian clock in Arabidopsis. Plant Cell 2008, 20, 1244–1259. [Google Scholar] [CrossRef] [Green Version]
  44. Zeng, J.; Mo, Y.L.; Chen, J.J.; Li, C.M.; Zhao, L.; Liu, Y.H. Expression and interaction proteins analysis of BjuFKF1 in stem mustard. Sci. Hortic. 2020, 269, 6. [Google Scholar] [CrossRef]
  45. Asp, T.; Frei, U.K.; Didion, T.; Nielsen, K.K.; Lubberstedt, T. Frequency, type, and distribution of EST-SSRs from three genotypes of Lolium perenne, and their conservation across orthologous sequences of Festuca arundinacea, Brachypodium distachyon, and Oryza sativa. BMC Plant Biol. 2007, 7. [Google Scholar] [CrossRef] [Green Version]
  46. Yan, Z.; Wu, F.; Luo, K.; Zhao, Y.; Yan, Q.; Zhang, Y.; Wang, Y.; Zhang, J. Cross-species transferability of EST-SSR markers developed from the transcriptome of Melilotus and their application to population genetics research. Sci. Rep. 2017, 7, 17959. [Google Scholar] [CrossRef]
  47. Pan, L.; Huang, T.; Yang, Z.; Tang, L.; Cheng, Y.; Wang, J.; Ma, X.; Zhang, X. EST-SSR marker characterization based on RNA-sequencing of Lolium multiflorum and cross transferability to related species. Mol. Breed. 2018, 38. [Google Scholar] [CrossRef]
  48. Zhang, Z.; Xie, W.; Zhao, Y.; Zhang, J.; Wang, N.; Ntakirutimana, F.; Yan, J.; Wang, Y. EST-SSR marker development based on RNA-sequencing of E. sibiricus and its application for phylogenetic relationships analysis of seventeen Elymus species. BMC Plant Biol. 2019, 19, 235. [Google Scholar] [CrossRef]
  49. Wang, Z.; Yu, G.; Shi, B.; Wang, X.; Qiang, H.; Gao, H. Development and characterization of simple sequence repeat (SSR) markers based on RNA-sequencing of Medicago sativa and in silico mapping onto the M. truncatula genome. PLoS ONE 2014, 9, e92029. [Google Scholar] [CrossRef]
  50. Zhang, M.A.; Mao, W.H.; Zhang, G.P.; Wu, F.B. Development and characterization of polymorphic EST-SSR and genomic SSR markers for Tibetan annual wild barley. PLoS ONE 2014, 9, e94881. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Liu, Z.P.; Liu, P.; Luo, D.; Liu, W.X.; Wang, Y.R. Exploiting illumina sequencing for the development of 95 novel polymorphic EST-SSR markers in common vetch (Vicia sativa subsp sativa). Molecules 2014, 19, 5777–5789. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Hajkova, L.; Koznarova, V.; Mozny, M.; Bartosova, L. Influence of climate change on flowering season of birch in the Czech Republic. Int. J. Biometeorol. 2020, 64, 791–801. [Google Scholar] [CrossRef] [PubMed]
  53. Grenier, C.; Bramel-Cox, P.J.; Hamon, P. Core collection of sorghum: I. Stratification based on eco-geographical data. Crop. Sci. 2001, 41, 234–240. [Google Scholar] [CrossRef]
  54. Stuerz, S.; Shrestha, S.P.; Schmierer, M.; Vu, D.H.; Hartmann, J.; Sow, A.; Razafindrazaka, A.; Abera, B.B.; Chuma, B.A.; Asch, F. Climatic determinants of lowland rice development. J. Agron. Crop. Sci. 2020, 206, 466–477. [Google Scholar] [CrossRef]
  55. Burgarella, C.; Chantret, N.; Gay, L.; Prosperi, J.M.; Bonhomme, M.; Tiffin, P.; Young, N.D.; Ronfort, J. Adaptation to climate through flowering phenology: A case study in Medicago truncatula. Mol. Ecol. 2016, 25, 3397–3415. [Google Scholar] [CrossRef]
  56. Wang, Q.B.; Wang, Y.P.; Zhang, L. Inheritance and molecular marker for flowering time in radish (Raphanus sativus L.). Plant Mol. Biol. Rep. 2018, 36, 878–887. [Google Scholar] [CrossRef]
  57. Xie, W.G.; Robins, J.G.; Bushman, B.S. A genetic linkage map of tetraploid orchardgrass (Dactylis glomerata L.) and quantitative trait loci for heading date. Genome 2012, 55, 360–369. [Google Scholar] [CrossRef]
  58. Wang, N. Evaluation of Flowering Time and Mining of Candidate Genes in Elymus Sibiricus. Master’s Thesis, Lanzhou University, Lanzhou, China, 1 May 2020. [Google Scholar] [CrossRef]
  59. Beier, S.; Thiel, T.; Munch, T.; Scholz, U.; Mascher, M. MISA-web: A web server for microsatellite prediction. Bioinformatics 2017, 33, 2583–2585. [Google Scholar] [CrossRef] [Green Version]
  60. Koressaar, T.; Remm, M. Enhancements and modifications of primer design program Primer3. Bioinformatics 2007, 23, 1289–1291. [Google Scholar] [CrossRef] [Green Version]
  61. Guerra, V.; Beule, L.; Lehtsaar, E.; Liao, H.L.; Karlovsky, P. Improved protocol for DNA extraction from subsoils using phosphate lysis buffer. Microorganisms 2020, 8, 532. [Google Scholar] [CrossRef] [Green Version]
  62. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  63. Paul, J. Nouvelles recherches sur la distribution florale. Bull. Soc. Vaud. Sci. Nat. 1908, 44, 223–270. [Google Scholar] [CrossRef]
  64. Hudson, R.R. Gennetic data-analysis-methods for discrete population genetic data-Weir, BS. Science 1990, 250, 575. [Google Scholar] [CrossRef]
Plants 09 01371 g001 550
Figure 1. Volcano plot of DEGs (differentially expressed genes) in booting stage, heading stage, and flowering stage of early-flowering and late-flowering accessions. (a) Early booting vs. Late booting, (b) Early heading vs. Late heading, and (c) Early flowering vs. Late flowering.
Figure 1. Volcano plot of DEGs (differentially expressed genes) in booting stage, heading stage, and flowering stage of early-flowering and late-flowering accessions. (a) Early booting vs. Late booting, (b) Early heading vs. Late heading, and (c) Early flowering vs. Late flowering.
Plants 09 01371 g001
Plants 09 01371 g002 550
Figure 2. Different types of SSR repeats and their percentage based on RNA-seq.
Figure 2. Different types of SSR repeats and their percentage based on RNA-seq.
Plants 09 01371 g002
Plants 09 01371 g003 550
Figure 3. The number of candidate genes involved in six flowering regulatory pathways.
Figure 3. The number of candidate genes involved in six flowering regulatory pathways.
Plants 09 01371 g003
Plants 09 01371 g004 550
Figure 4. The flowering days of 20 E. sibiricus germplasms from 1 January 2018.
Figure 4. The flowering days of 20 E. sibiricus germplasms from 1 January 2018.
Plants 09 01371 g004
Plants 09 01371 g005 550
Figure 5. Amplified results of primer 28366 in 20 E. sibiricus germplasms. (a) The result of PAGE (polyacrylamide gel electrophoresis) of primer 28366 in 20 E. sibiricus germplasms. (b) Multiple sequence alignment of PCR (polymerase chain reaction) products of primer 28366 in 20 E. sibiricus germplasms.
Figure 5. Amplified results of primer 28366 in 20 E. sibiricus germplasms. (a) The result of PAGE (polyacrylamide gel electrophoresis) of primer 28366 in 20 E. sibiricus germplasms. (b) Multiple sequence alignment of PCR (polymerase chain reaction) products of primer 28366 in 20 E. sibiricus germplasms.
Plants 09 01371 g005
Plants 09 01371 g006 550
Figure 6. Multiple sequence alignment of RNA-seq and in the late-flowering and early-flowering E. sibiricus accessions. (a) Comparative electropherogram analysis of two EST-SSR loci (42789 and 33508) in early-flowering and late-flowering E. sibiricus accessions. (b) Alignment of sequences obtained from selected PCR bands amplified by two primers (32835 and 41988).
Figure 6. Multiple sequence alignment of RNA-seq and in the late-flowering and early-flowering E. sibiricus accessions. (a) Comparative electropherogram analysis of two EST-SSR loci (42789 and 33508) in early-flowering and late-flowering E. sibiricus accessions. (b) Alignment of sequences obtained from selected PCR bands amplified by two primers (32835 and 41988).
Plants 09 01371 g006
Plants 09 01371 g007 550
Figure 7. Cluster analysis of 20 E. sibiricus accessions with 15 candidate gene-based EST-SSR markers and the heatmap analysis of flowering days of the 20 E. sibiricus accessions.
Figure 7. Cluster analysis of 20 E. sibiricus accessions with 15 candidate gene-based EST-SSR markers and the heatmap analysis of flowering days of the 20 E. sibiricus accessions.
Plants 09 01371 g007
Plants 09 01371 g008 550
Figure 8. Population structure of 20 E. sibiricus accessions inferred from the STRUCTURE program with 15 candidate gene-based EST-SSR markers data set. (a) Mean Ln P(D) over 20 runs for each K value, and maximum delta K (ΔK) values were used to determine the uppermost level of structure for K ranging from 1 to 11, where K is three and three clusters; (b) three major groups of 20 E. sibiricus genotypes. The vertical coordinate of each group indicates the membership coefficients for each genotype. Different codes and corresponding vertical lines represent individual genotype and different colors represent stock.
Figure 8. Population structure of 20 E. sibiricus accessions inferred from the STRUCTURE program with 15 candidate gene-based EST-SSR markers data set. (a) Mean Ln P(D) over 20 runs for each K value, and maximum delta K (ΔK) values were used to determine the uppermost level of structure for K ranging from 1 to 11, where K is three and three clusters; (b) three major groups of 20 E. sibiricus genotypes. The vertical coordinate of each group indicates the membership coefficients for each genotype. Different codes and corresponding vertical lines represent individual genotype and different colors represent stock.
Plants 09 01371 g008
Table
Table 1. Statistical table of SSR (simple sequence repeat) analysis results.
Table 1. Statistical table of SSR (simple sequence repeat) analysis results.
Searching ItemNumber
Total number of sequences examined40,639
Total size of examined sequences (bp)78,157,064
Total number of identified SSRs13,052
Number of SSR containing sequences10,286
Number of sequences containing more than 1 SSR2176
Number of SSRs present in compound formation742
Mono-nucleotide5471
Di-nucleotide2642
Tri-nucleotide4551
Tetra-nucleotide314
Penta-nucleotide48
Hexa-nucleotide26
Table
Table 2. He (gene diversity), Ho (heterozygosity), and PIC (polymorphic information content) of 15 polymorphic primers amplified in 20 E. sibiricus germplasms.
Table 2. He (gene diversity), Ho (heterozygosity), and PIC (polymorphic information content) of 15 polymorphic primers amplified in 20 E. sibiricus germplasms.
Primer CodeFlowering GenesFlowering Pathways/FunctionsHeHoPIC
Primer 46865FPF1Flowering0.65920.92310.4767
Primer 33680FPF1Flowering0.22680.15000.1275
Primer 43088ELF3Flowering0.79660.94120.4210
Primer 41496GID1Gibberellin0.49110.50000.2367
Primer 40505GA2OX6Gibberellin0.66550.70000.2513
Primer 37852CIGRGibberellin0.19750.05880.2500
Primer 34975XCTCircadian Clock0.51700.80000.1567
Primer 48198GIGANTEACircadian Clock0.65841.00000.1400
Primer 34500MBD9CONSTANS-Like0.87311.00000.2606
Primer 31670NFYC4CONSTANS-Like0.32000.59460.1875
Primer 43287COL4CONSTANS-Like0.71380.97300.3325
Primer 36067COL13CONSTANS-Like0.81170.97300.3025
Primer 36927FLCCentral Integrator0.81471.00000.1500
Primer 34261HD3Central Integrator0.77941.00000.2590
Primer 28366HD3aCentral Integrator0.65091.00000.1600
Mean 0.61170.77420.2475
Table
Table 3. The origins, longitude, latitude, and altitude of 20 E. sibiricus accessions.
Table 3. The origins, longitude, latitude, and altitude of 20 E. sibiricus accessions.
CodeAccession NumberSample IDOriginsLongitude/LatitudeAltitude (m)
1PI598781Z1Baikal, Buryat, Russia--
2PI610876Z2---
3PI610994Z3Siberia, Russia-1250
4PI636676Z4Xiahe, Gansu, China35.19 N, 102.65 E2720
5PI655199Z5Hongyuan, Sichuan, China32.08 N, 102.57 E3280
6W610305Z6Siberia, Russia--
7W630476Z7Kyrgyzstan--
8HZ03Z8Hezuo, Gansu, China--
9LT05Z9Lintan, Gansu, China--
10LQ02Z10Luqu, Gansu, China--
11PI435088W1Ninshan, Xinjiang, China--
12PI504462W2Xining, Qinghai, China--
13PI531665W3Beijing, China--
14PI531669W4Xining, Qinghai, China-2400
15PI595149W5Xinjiang, China44.17 N, 84.57 E1500
16PI595156W6Xinjiang, China45.03 N, 81.12 E1300
17PI595162W7Xinjiang, China44.12 N, 87.97 E1680
18PI595169W8Xinjiang, China43.77 N, 89.45 E1300
19PI595174W9Xinjiang, China43.68 N, 89.30 E1870
20PI655092W10Mongolia49.80 N, 94.96 E1370
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zheng, Y.; Zhang, Z.; Wan, Y.; Tian, J.; Xie, W. Development of EST-SSR Markers Linked to Flowering Candidate Genes in Elymus sibiricus L. Based on RNA Sequencing. Plants 2020, 9, 1371. https://doi.org/10.3390/plants9101371

AMA Style

Zheng Y, Zhang Z, Wan Y, Tian J, Xie W. Development of EST-SSR Markers Linked to Flowering Candidate Genes in Elymus sibiricus L. Based on RNA Sequencing. Plants. 2020; 9(10):1371. https://doi.org/10.3390/plants9101371

Chicago/Turabian Style

Zheng, Yuying, Zongyu Zhang, Yiyang Wan, Jiaoyang Tian, and Wengang Xie. 2020. "Development of EST-SSR Markers Linked to Flowering Candidate Genes in Elymus sibiricus L. Based on RNA Sequencing" Plants 9, no. 10: 1371. https://doi.org/10.3390/plants9101371

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop