Skip to main content
  • Research article
  • Open access
  • Published:

Construction of the first high-density genetic linkage map and identification of seed yield-related QTLs and candidate genes in Elymus sibiricus, an important forage grass in Qinghai-Tibet Plateau



Elymus sibiricus is an ecologically and economically important perennial, self-pollinated, and allotetraploid (StStHH) grass, widely used for forage production and animal husbandry in Western and Northern China. However, it has low seed yield mainly caused by seed shattering, which makes seed production difficult for this species. The goals of this study were to construct the high-density genetic linkage map, and to identify QTLs and candidate genes for seed-yield related traits.


An F2 mapping population of 200 individuals was developed from a cross between single genotype from “Y1005” and “ZhN06”. Specific-locus amplified fragment sequencing (SLAF-seq) was applied to construct the first genetic linkage map. The final genetic map included 1971 markers on the 14 linkage groups (LGs) and was 1866.35 cM in total. The length of each linkage group varied from 87.67 cM (LG7) to 183.45 cM (LG1), with an average distance of 1.66 cM between adjacent markers. The marker sequences of E. sibiricus were compared to two grass genomes and showed 1556 (79%) markers mapped to wheat, 1380 (70%) to barley. Phenotypic data of eight seed-related traits (2016–2018) were used for QTL identification. A total of 29 QTLs were detected for eight seed-related traits on 14 linkage groups, of which 16 QTLs could be consistently detected for two or three years. A total of 6 QTLs were associated with seed shattering. Based on annotation with wheat and barley genome and transcriptome data of abscission zone in E. sibiricus, we identified 30 candidate genes for seed shattering, of which 15, 7, 6 and 2 genes were involved in plant hormone signal transcription, transcription factor, hydrolase activity and lignin biosynthetic pathway, respectively.


This study constructed the first high-density genetic linkage map and identified QTLs and candidate genes for seed-related traits in E. sibiricus. Results of this study will not only serve as genome-wide resources for gene/QTL fine mapping, but also provide a genetic framework for anchoring sequence scaffolds on chromosomes in future genome sequence assembly of E. sibiricus.


The tribe Triticeae (Poaceae) includes several major cereal crops (wheat, barley, and rye) and many ecologically and economically important forage grasses [1]. Elymus L. is the largest genus in the Triticeae, which comprises about 150 polyploid perennial grass species widely distributed worldwide [2]. Asia is the most important center of origin where approximately 80 Elymus species were found [3]. Many Elymus species are closely related to wheat and barley, and may thus serve as potential gene pool for the improvement of stress tolerance (cold, drought and disease) and other important agronomic traits [4]. Elymus sibiricus (Siberian wild rye), which is indigenous to northern Asia, is an important perennial, cold-season and self-pollinated forage grass of the genus Elymus [5]. Based on the cytogenetic analysis, E. sibiricus is allotetraploid species, containing St and H genomes. The St genome is derived from Pseudoroegneria spicata (Pursh) A. Löve, and the H genome is derived from the genus Hordeum [6]. Elymus sibiricus is widely grown and used for forage production and grassland eco-engineering in the Qinghai-Tibet Plateau region of China, owing to its good forage quality, drought and cold tolerance, and excellent adaptability to local special environments [7, 8]. Despite E. sibiricus has various agricultural uses and economically importance, its serious seed shattering makes seed production difficult for this species. For cereal crops and forage grasses, seed yield is affected by many seed yield-related traits, such as spike length, seed width, floret number per spike, 1000-seed weight, and seed shattering, among which seed shattering is a major cause of yield loss [9]. Previous study showed that serious seed shattering may result in up to 80% seed yield losses if harvesting is delayed [10]. As a result, selection for high seed retention and genetic improvement of seed shattering are important breeding objectives for this species. Several major quantitative trait loci (QTLs) and genes for seed shattering have been reported in cereal crops like rice, wheat, barley, maize and sorghum, and a few forage grasses. For example, in rice, SH4 [11], qSH1 [12], OsCPL1 [13], SHAT1 [14], and SH5 [15] were identified as major genes for seed shattering, their functions and interactions in regulating abscission layers formation and development were also revealed. In addition, in hybrid Leymus (Triticeae) Wildryes, a major-effect QTL for seed retention was identified on linkage group (LG) 6a, which aligns to other seed shattering QTLs in American wildrice, Zea and Triticum [16]. Together, these studies indicate the presence of QTLs and genes with large effects on seed shattering, and the potential to understand which QTLs or genes play a role in regulating seed shattering.

The availability of genetic map makes feasible the identification of genes for monogenic traits or major loci for quantitative traits, it also provides an important basis for the study of genome structure and evolution [17]. It is particularly important for future positional gene cloning, marker-assisted selection, and comparative genome analysis [18]. The utility of genetic linkage map depends on the types and number of markers used [19]. High-density linkage map lays a foundation for genome assembly and fine mapping of quantitative trait loci (QTL) [20]. To date, several molecular marker systems have been used for the construction of genetic linkage map, including amplified fragment length polymorphism (AFLP) [21], restriction fragment length polymorphisms (RFLP) [22], random amplified polymorphic DNA (RAPD) [23], simple sequence repeat (SSR) [24], sequence-related amplified polymorphism (SRAP) [25], and single-nucleotide polymorphism (SNP) [26]. Among these markers, SNP marker is considered as the most promising molecular marker for high-density genetic map construction due to their abundant and wide distribution in genome. The advent of massive parallel next-generation sequencing (NGS) technologies could identify and obtain thousands of SNPs at the whole genome level, thus making it possible to construct high-density SNP genetic maps. However, whole-genome sequencing and genotyping large populations are still cost-prohibitive [27]. Reduced representation library sequencing is considered to be one efficient strategy to bring down the cost through genome reduction [28, 29]. For example, restriction site-associated sequencing (RAD-seq) sequences only the DNA fragment with restriction sites, and has been used for large-scale SNP discovery and genetic mapping in many species [30, 31]. As a modified reduced representation sequencing technique, specific-locus amplified fragment sequencing (SLAF-seq) has several distinguishing advantages such as reduced sequencing costs, deep sequencing, marker efficiency optimization through pre-designed reduced representation scheme, and double-barcode method for large populations. It is an efficient method for large-scale De Novo SNP discovery and genotyping of large population [32]. Recently, SLAF-seq has been increasingly used for high-density genetic linkage map construction in several crops [33], forage grasses [20], and animal species [34].

Toward improving the understanding of E. sibiricus genome arrangement and the genetic control of seed yield-related traits, we constructed a genetic linkage map and identified QTLs related to seed shattering as well as other seed traits. Two E. sibiricus genotypes were selected based on their variation for seed yield-related traits. We applied SLAF-seq to develop thousands of SLAF markers (SLAFs) and construct the first high-density genetic linkage map in E. sibiricus, then identified QTLs and candidate genes for eight seed yield-related traits. These results could lay a foundation for future functional genetic dissection of key genes related to seed shattering and other seed traits.


Analysis of SLAF-seq and SLAF markers

After SLAF library construction and high-throughput sequencing, 253.25 Gb of raw data containing 1267.20 M reads were generated. The average percentage of Q30 (quality scores of at least 30) bases was 93.03%. The average guanine-cytosine (GC) content was 46.69%. To estimate the validity of library construction, we used Oryza sativa ssp. japonica (genome size = 382 M) as control. A total of 901,095 reads with 92.17% Q30 bases and 45.32% GC content were generated (Table 1). The number of reads for male and female parents was 29,809,327 and 65,542,805, respectively. The average number of reads for offspring was 5,859,224.46 with 93.03% Q30 bases and 46.69% GC content. The number of SLAF markers generated for male and female were 232,429 and 326,923, respectively. The average number of SLAF marker in the progeny was 202,120 (Table 2). The average sequencing depth was 31.95-fold and 7.51-fold for parents and each progeny, respectively.

Table 1 Summary of SLAF sequencing data
Table 2 Summary of SLAF tag information

We detected 370,470 SLAF markers, among which 97,387 were polymorphic, 269,579 and 3504 were non-polymorphic (72.77%) and repetitive (0.94%), respectively. Polymorphic markers included mapped biallelic markers and unmapped biallelic markers, monomorphic markers with only one tag in parents were recognized as non-polymorphic markers, mutiallelic markers with tag number larger than 4 in parents were recognized as repetitive markers. Mutiallelic SLAFs which could not be used for recombination rate calculating were removed from further analysis. After filtering the SLAF markers lacking the parent information, 46,135 polymorphic SLAFs were successfully genotyped and further classified into eight segregation patterns (ab×cd, ef × eg, lm × ll, nn × np, aa×bb, hk × hk, cc × ab, ab×cc) (Fig. 1). The mapping population was obtained from the F1 hybrid plant of two homozygous parents, therefore, the 18,343 SLAF markers with aa×bb segregation pattern in the F2 population were used for genetic map construction.

Fig. 1
figure 1

Number of markers for eight segregation patterns

Basic characteristics of the genetic maps

We further filtered the SLAF markers using four criteria [20]. These SLAF markers that belonging to following four types were removed from mapping construction: SLAF markers from parents with sequencing depth less than 10X; SLAF markers with more than five SNPs; SLAF markers with missing in more than 10% of offspring and segregation-distorted markers (Chi-square, p < 0.01). Only the SLAF markers that passed the four-step filtering process were used for constructing a high-quality genetic map. The final map included 1971 markers with 2610 SNP on the 14 linkage groups (LGs) and was 1866.35 cM in length (Fig. 2). The length of each linkage group ranged from 87.67 cM (LG7) to 183.45 cM (LG1), with an average marker density of 1.66 cM between adjacent markers (Table 3). The maximum number of markers (565) were found on LG11, whereas LG8 possessed the minimum number of markers (29) (Additional file 4: Figure S1, Additional file 1: Table S1). The “Gap ≤ 5” value was used to reflect the degree of linkage between each marker, ranging from 73.08 to 100%, with an average of 92.09%. The largest gap on this map was 11.03 cM located in LG14. The number of SNP on each linkage group varied from 35 (LG7) to 712 (LG 11), with an average of 186.

Fig. 2
figure 2

Distribution of SLAF markers on the 14 linkage maps

Table 3 Description of basic characteristics of the 14 linkage maps

In total, only 26 markers showed a significant (p < 0.05) segregation distortion and were mapped on the final map, accounting for 1.32% of mapped markers (Table 4). Most of the linkage groups (LGs) had segregation distortion markers with the exceptions of LG1, LG3, LG4, LG13, and LG14. The frequencies of distorted markers on LG6 (19.23%) and LG12 (19.23%) were higher than those of the other linkage groups. LG11, which possessed the maximum mapped markers (565 SLAF markers), had the lowest frequency of distorted marker (3.85%).

Table 4 Distribution of segregation distortion markers on each linkage group

Quality evaluation of the genetic map

To evaluate the quality of the genetic map, haplotype mapping and heat mapping were carried out. The haplotype map reflected the double exchange of the population, which is caused by genotyping error, suggesting a possible recombination hotspot. The haplotype maps of each linkage group were developed for the parental controls and 200 offspring using 1971 SLAF markers. The results showed that most of the recombination blocks were distinctly defined. The LGs 9, 10, and 13 had no missing data, while LG 8 had the largest missing data (3.53%), with an average of 0.73%. Most of the LGs were uniformly distributed (Additional file 5: Figure S2). The heat maps were constructed based on the pair-wise recombination value from the 1971 mapped markers to reflect the recombination relationship between mapped markers on each single linkage group (Additional file 6: Figure S3). The results confirmed the order of mapped SLAF markers on each linkage group.

Phenotypic variation

Phenotypic analysis of the parents and F2 population revealed significant variations in all eight seed yield-related traits (Table 5, Additional file 2: Table S2). The coefficient of variation (CV) among all traits ranged from 7.24% (WS in 2018) to 58.08% (FN in 2016). We analyzed the correlation between years and traits (Table 6). Our results showed a correlation between phenotypic data detected in different years with exception of WS between 2016 and 2017, and SW1 between 2016 and 2017. For example the correlation for seed shattering (SSc) between 2016 and 2018, 2017 and 2018, 2016 and 2017 were 0.841, 0.783, and 0.360, respectively. Floret number per spike (FN) was significant correlated between 2016 and 2018, 2017 and 2018. Spike length (SL) was significant correlated during 3 years. We calculated the heritability of these traits, all traits had relatively high heritability. The highest heritability (0. 6718) was found for seed shattering (SSc), the lowest heritability (0.4638) was found for floret number per spike (FN). These results were consistent with the correlation analysis between different years. The correlation were found between most traits, for example, awn length (AL) was positively correlated with width of seed (WS), 1000-seed weight (SW1) and spike length (SL). Seed shattering (SS) was positively correlated with floret number per spike (FN). The absolute values of Skewness and Kurtosis for most traits with exception of FN (2017), WS (2017 and 2018), and SW1 (2017) were less than 1 (Table 5). Besides, the normal frequency distributions of eight traits were analyzed and the P-value was more than 0.05 except for SL (2017), FN, SS, WS (2017 and 2018) and SW1 (2017) (Fig. 3).

Table 5 Descriptive statistics for seed-related traits in the two parents and F2 population
Table 6 The correlation analysis between three years and eight seed-related traits among F2 population
Fig. 3
figure 3

The frequency distribution of eight seed yield-related traits in the F2 population. The x-axis shows the ranges of phenotypic traits and the y-axis represents the number of individuals in the F2 population

QTL mapping and comparative genome analysis

A total of 29 QTLs were detected for eight seed-related traits on 14 linkage groups, of which 3 for spike length (SL), 2 for floret number per spike (FN), 6 for seed shattering (SS, SSD and SSc), 7 for awn length (AL), 3 for width of seed (WS), and 8 for 1000 seed weight (SW1). The LOD and PVE (the percentage of phenotypic variation explained) for all QTLs ranged from 3 to 10.62, 2.17 to 10.85%, respectively (Fig. 4, Table 7). Six QTLs detected for seed shattering explained 2.17 to 9.48% of the phenotypic variation. Among the six QTLs, 1 QTLs were detected on LGs 6 using breaking tensile strength (BTS) data, 2 QTLs were detected on LGs 3 and 11 using seed shattering degree (SSD) data, 3 QTLs were detected on LGs 2, 3 and 11 using seed shattering rate (SSc) data. Especially, seed shattering QTLs on LG3 and LG11 could be detected using two methods and at two years (2016 and 2017), respectively. Seven QTLs for awn length (AL) were detected on five linkage groups (LG1, LG5, LG6, LG11 and LG13), among which the QTL on LG1 explained the maximum phenotypic variation of 10.37%. On LG12, a QTL for seed width (WS) was detected and explained the largest phenotypic variation of 10.85% among all QTLs. Moreover, QTLs for awn length (AL) and 1000 seed weight (SW1) were detected on more than five LGs, suggesting a complex genetic mechanism of these traits. A total of 16 QTLs could be consistently detected for two or three years, for example, two QTLs for spike length (SL) on LG14 were detected in 2017 and 2018, two QTLs for seed shattering on LG11 were detected in 2016 and 2017, three QTLs for 1000-seed weight (SW1) on LG9 and three QTLs for awn length (AL) on LG1 were detected for three years.

Fig. 4
figure 4

Quantitative trait loci (QTLs) for eight seed yield-related traits. Each QTLs were compared with barley and wheat genomes, respectively

Table 7 Seed-related QTLs detected in F2 population of E. sibiricus and a comparative genome analysis with Barley and Wheat

The 1971 mapped SLAF markers generated from E. sibiricus were compared with the genome sequences of wheat and barley. The Circos plot and Colinear graph was constructed to show the linear relationships between E. sibiricus and wheat and barley, illustrating a corresponding relationship between the mapped markers and their genomic locations (Fig. 5). The numbers of matching markers between E. sibiricus and each species were 1556 (79%) for wheat, 1380 (70%) for barley (Fig. 5a). We further broken down alignments to each subgenome of wheat (A, B and D), the number of matching markers on ChrA, ChrB and ChrD was 311, 523 and 521, respectively (Fig. 5b and c). The largest number of matching markers was found on Chr2A (79), CHr3B (129) and Chr7D (98) for each subgenome, respectively. For barely, the number of matching marker on each chromosome ranged from 139 (Chr1) to 237(Chr2), with an average of 184.

Fig. 5
figure 5

Comparative genome analysis between E. sibiricus and other two grass species. a, circos plot showing linear relationship between E. sibiricus with barley (left) and wheat (right); b, colinear graph between E. sibiricus and three subgenomes of wheat (A, B and D); c, the number of matching markers on ChrA, ChrB and ChrD

A total of 43 markers within seed yield-related QTL regions could be identified on barley and wheat chromosomes, of which 3 markers for spike length (SL), 2 markers for floret number per spike (FN), 8 markers for seed shattering (SS, SSD and SSC), 7 markers for awn length (AL), 6 markers for seed width (WS), and 17 markers for 1000-seed weight (SW1) (Table 7). These markers were distributed across different chromosomes of each species. For example, 8 markers linked with seed shattering distributed on wheat chromosomes 1A, 3B, 6A, 6B, and 7B, and barley chromosomes 1H, 2H, 3H and 5H. We further identified 30 candidate genes for seed shattering within six QTL regions based on the functional annotation of barley and wheat genomes, of which 15, 7, 6 and 2 genes were involved in plant hormone signal transcription, transcription factor, hydrolase activity and lignin biosynthetic process, respectively (Fig. 6, Additional file 3: Table S3). In particular, among candidate genes for plant hormone, 3 genes were involved in regulation of abscisic acid-activated signaling pathway, 6 genes were involved in ethylene response pathway, 3 genes were involved in auxin-activated signaling pathway, 1 genes for gibberellin, and 1 gene for jasmonic acid. Based on our abscission zone transcriptome data of two E. sibiricus genotypes (XH, high seed shattering, ZN, low seed shattering). A total of 20 unigenes involved in plant hormone, transcription factor, and hydrolase activity were predicted from “XH-WS vs ZN-WS”, of which 14 genes were up regulated in high seed shattering genotype XH, 6 genes (2 for ethylene activity, 1 for gibberellin activity, 1 for MYB transcription factor activity, 1 for xylanase activity and 1 for glycosyl hydrolase activity) were up regulated in low seed shattering genotype (Fig. 7). Together, these results suggested these candidate genes might be associated with the regulation of seed shattering in E. sibiricus.

Fig. 6
figure 6

The mapping position and thirty identified candidate genes for seed shattering QTLs on LG3 and LG11

Fig. 7
figure 7

Heatmap diagram of the expression levels of 20 differentially expressed genes involved in seed shattering of two E. sibiricus genotypes. XH, high seed shattering genotype, ZN, low seed shattering genotype. WS, 28 days after heading


The first linkage map for E. sibiricus by next-generation sequencing

E. sibiricus breeders have been working to improve economically important traits including yield potential, stress tolerance and persistence. To understand the molecular and genetic mechanisms underlying these important traits, we need some modern methods and tools to effectively identify QTLs and candidate genes involved in these traits. Molecular marker development and genetic linkage map construction are important preliminary basic works for undertaking molecular breeding activities in any crops [35]. Especially, markers strongly related to preferred traits could be used in marker-assisted selection (MAS) to speed up the genetic improvement of desired agronomic traits. To date, a variety of molecular markers including RAPD, SRAP, AFLP and SSRs have been wildly used for linkage map construction in a number of forage grasses such as tall fescue [36], meadow fescue [18], ryegrass [37], orchardgrass [38, 39], alfalfa [40], white clover [41], red clover [42], and Elymus wheatgrass [2]. However, the number of mapped markers are limited in many previously reported maps. SLAF sequencing, as a next-generation sequencing technology for genome-wide SNP discovery, has been successfully used for high-density genetic map construction in many plant species including rice [33], cucumber [26], kiwifruit (Actinidia chinensis) [43], mei (Prunus mume) [44], sesame [45], even in animals such as chicken [46], and white shrimp [34]. However, large scale SNP mining in E. sibiricus lagged behind other species due to its large, complex, and nature polyploidy. Here, we used the SLAF-seq to develop and identify 370,470 SLAF markers, of which 97,387 were polymorphic with a polymorphism rate of 26.29%. The polymorphism rate of SLAF markers between the two parents was 26.29%, higher than previous reports in cucumber (9.57%) [26] and sesame (5.12%) [45], suggesting a considerable difference between the two parental genotypes. The linkage map contained 14 linkage groups and spanned 1866.35 cM. Based on SLAF-seq, the map quality of the present genetic maps was similar to previously reported genetic maps for other several species, although, they were the first ones reported for this species. For example, Zhang et al. [45] reported the first high-density genetic map for sesame. A total of 1233 markers were mapped on the 15 linkage groups, with an average marker density of 1.20 cM. In general, these results proved that SLAF-seq is a powerful high-throughput technology for the whole-genome wide SNP discovery and is effective for E. sibiricus linkage map construction.

According to our results, the number of SLAF markers on each linkage group varied from 27 (LG7) to 565 (LG 11). Especially, some markers tended to highly cluster in some regions on LG9, LG10 and LG 11. The similar results were found in other plants such as grape [47], sunflower [48], and tree peony [49]. This phenomenon may result from the non-random distribution of mapped markers on linkage group and the uneven recombination rates and marker polymorphism between mapping parents on some chromosomes [50]. Moreover, four gaps larger than 10 cM were located on LG1, LG7, LG8, and LG14. The lack of maker polymorphism and a shortage of marker detection in these regions may have contributed to this finding [51, 52].

Application of this map for the QTL detection of seed yield-related traits

Like most native grasses, E. sibiricus has serious seed shattering which cause large seed yield losses during harvest, making commercial seed production difficult. In this study, Y1005 and ZhN06 were selected and used as parents because they are genetically divergent and have several contrasting seed traits. In the F2 mapping population, eight seed yield-related traits showed considerable phenotypic variation. For example, seed shattering in mapping population ranged from 5.14 gf to 18.80 gf during 2016, from 5.66 gf to 20.68 gf during 2017, and from 6.53 gf to 19.62 gf during 2018. Some progenies had lower seed shattering than their parents. Thus, the use of low seed shattering genotype provides QTLs for low seed shattering that could prove robust for breeding lower seed shattering genotypes.

Genetic linkage map allows for comparative genetic studies with other species, provides important information about the genome structure and evolution of a species, and lays a foundation for studying complex and important agronomic traits [53]. The major objective of this study was to reveal the genetic mechanisms controlling seed yield-related traits in E. sibiricus. Our comparative genome analysis indicated that E. sibiricus was more closely related to wheat and barley. In this study, a total of 8 markers within seed shattering QTL (SS, SSD and SSC) regions could be identified on wheat chromosomes 1A, 3B, 6A, 6B, and 7B, and barley chromosomes 1H, 2H, 3H and 5H. In wheat, the Q gene on chromosome 5A was identified as a major domestication gene, which encodes an AP2 transcription factor largely responsible for free seed shattering [54]. TaqSH1, encoding a BEL1-like protein, was located on the homoeologous group 3 chromosome in wheat. Overexpression of TaqSH1 gene in transgenic Arabidopsis plants could down-regulate several well-known abscission-related genes such as HAESA and KNAT1/6, suggesting this gene might play as a key upstream regulator for abscission zone development [55]. In barley, the non-brittle rachis trait was controlled by btr1 and btr2 genes on chromosome 3H. Further, Pourkheirandish et al. [56] cloned and identified the Btr1 and Btr2 genes and elucidated the mechanism underlying the disarticulation of the wild type barley spike.

Based on our comparative genome analysis, marker 5164 and marker 144,585 were on wheat Chr3B, marker 124,682 was on barely Chr 3H. Together, these results highlight the possibility for seed shattering candidate gene identification in E. sibiricus based on these genetic maps and the synteny with model grasses. QTLs on other genomic regions may provide new candidate genes for seed shattering in E. sibiricus.

Candidate genes for seed shattering

Plant hormone play an important role in regulating plant growth and development processes. In this study, we identified 15 potential seed shattering candidate genes involved plant hormone including abscisic acid (ABA), ethylene, auxin, gibberellin and jasmonic acid. Abscisic acid regulates many agronomically important development processes and numerous adaptive stress responses in plants [57]. It plays a direct role in abscission of many plant organs such as seed, leaf, flower and fruit [5, 58]. Our previous studies found some ABA-related genes were up-regulated in the abscission zone and suggested these ABA-responsive genes may affect seed shattering [5, 8]. Ethylene is an important regulator of abscission of many plant organ such as seeds, fruit and leaf [58]. Ethylene receptor genes (ETR1), the ethylene insensitive mutant of Arabidopsis, plays a role in delaying the shedding of floral parts [59]. In our previous study, we reported that 5 ETR1 gene were up-regulated in abscission zone in Elymus nutans, suggesting the roles of ETR genes in regulating abscission [60]. Ethylene-responsive factors (ERF) also regulate abscission of plant organs, including floral organs, leaf and seed. In tomato (Solanum lycopersicum) Ethylene-responsive factor 52 (SIERF52) is specifically expressed in pedicel abscssion zones. When the expression of SIERF52 was suppressed in transgenic plants flower abscission will be significantly delayed compared with wild type [61]. Our transcriptome analysis showed 4 ethylene response factor genes were up-regulated in the abscission zone of high seed shattering genotype, suggesting their roles in promoting seed shattering. Jasmonic acid (JA) is an important regulator of plant growth, development and defense. It had good abscission-promoting effect and positively promoted the abscission of bean petiole expiants in the dark and light without enhancing ethylene production in bean petiole expiants [62]. Similar results were reported in our study, we found three genes involved in JA-mediated signaling pathway were up-regulated in abscission zone of high seed shattering genotype XH-WS. However, it is difficult to identify which hormone or gene is key factor for regulating abscission process, as plant organ abscission is a complex and highly coordinated process involving multiple gene expressions in plant hormone signing pathway. A balance and interaction of these genes may have contributed to seed shattering in E. sibiricus.

Transcription factors play an important role in the signal transduction pathways. In this study we identified 7 transcription factor genes, such as MYB, MADS-box and WRKY genes. According to previous study in major crops. Many well-known genes for seed shattering are transcription factors genes like SH4 [11] and qSH1 [12], STK [63]. In rice, SH4 is a major QTL for seed shattering, which encodes a MBY3 transcription factor [11]. And many MYB proteins are critical components of multiple hormone-mediated transcriptional regulatory, including ethylene, abscission acid and auxin, which act as important regulators of plant organ abscission [64]. In our previous study, we identified 14 MYB genes up-regulated in the abscission zone in E. nutans, indicating their potential roles in regulating the development of abscission zone [60]. In Arabidopsis, STK is a MADS-box transcription factor gene, which regulates the formation of seed abscission zone [63]. WRKY transcription factors are key components in abscisic acid signaling, and play roles in regulating many plant processes, including seed development, the response to stress, and seed shattering [65]. Therefore, we inferred these candidate transcription factors identified in this study may have the similar functions in the regulation of seed shattering.

Seed shattering is generally caused by the development and degradation of abscission layers that is located in the rachilla just below each seed. Our previous study showed that increased hydrolytic enzymes activity like cellulase and polygalacturonase in abscission zone are highly related to high seed shattering degree in E. sibiricus [5]. Cellulase is a key hydrolytic enzyme, which plays role in plant cell wall loosening during plant organ abscission [66]. In rice, an endo-1,4,-β- glucanase gene named as OsCel9D, is an important regulator in modifying cell wall structure and component during abscission, and mutations of this gene will reduce cell elongation and affect cellulose biosynthesis and increase the pectin content, which finally hamper the abscission process in seed shattering [67]. In addition, xyloglucan endotransglucosylase (XTHs) were suggested to loosen plant cell wall through cutting and rejoining the xyloglucans that tether adjacent cellulose miscrofibrils [68]. In this study we identified 2 cellulase genes and 2 XTHs genes, and transcriptome analysis showed that these candidate genes were differently expressed in the abscission zone of two E. sibiricus genotypes, indicating these genes may have contributed to seed shattering.


In general, seed shattering is an important agronomic trait that need improvement during wild grass domestication. Seed shattering is a complex biological process affected by environment factors, cultivation management and multiple changes in the metabolism process, abscission layer cell structure, and functional gene expression level. Previous studies have reported many genes associated with seed shattering, they are involved in hydrolytic enzymes activity, lignin biosynthesis and degradation, plant hormone signaling and response, transcription factors and protein kinase activity [5, 8, 15, 58]. In this study, we constructed the first genetic linkage map and identified seed-related QTLs and candidate genes for seed shattering. Results from this study mostly confirmed previous findings and also reported some new potential candidate genes for seed shattering. More studies are needed to explore their potential roles and functions in seed shattering in the future. In addition, a combination of multiple methods, including genetic mapping and whole-genome association analysis and transcriptome analysis could help us to identify more major loci underlying seed shattering in E. sibiricus.

Materials and methods

Plant materials and DNA extraction

The mapping population of 200 F2 individuals was obtained from self-pollinating a single F1 plant, which was developed from a cross between “Y1005” (male parent) and “ZhN06” (female parent). Y1005 and ZhN06 were collected from Sichuan and Gansu provinces, China, respectively. The parents were selected based on a previous evaluation for agronomic traits and genetic diversity [69]. The two parental genotypes are genetically divergent and have several contrasting seed traits, including seed shattering, 1000-seed weight, seed width, spike length, awn length, and floret number per spike.

DNA extraction was carried out with young healthy leaves using the Qiagen DNeasy 96-well procedure (QIAGEN, Valencia, Calif). The quantity and quality of genomic DNA samples were evaluated using the NanoDrop ND1000 spectrophotometer (NanoDrop, Wilmington, DE, USA) and by 0.8% agarose gel electrophoresis.

Phenotypic evaluation

The F2 population of 200 individuals and parents were grown and evaluated in the field of Yuzhong research farm, Lanzhou University, Gansu, China (elevation 1720 m, longitude 103°34′ E, latitude 35°34′ N) for mapping. Plants were spaced 0.5 m within rows and 0.5 m between rows. Phenotypic data for seed yield-related traits were evaluated, including spike length (SL), floret number per spike (FN), seed shattering (SS), awn length (AL), width of seed (WS), and 1000-seed weight (SW1) for three consecutive years (2016–2018). E. sibiricus plant has a spike inflorescence containing 15–30 spikelets, each spikelet consists of approximately 5 florets [5]. Three previously reported methods were used to accurately evaluate seed shattering in this study. Breaking tensile strength (BTS) method was firstly used to determine seed shattering. The BTS value was measured upon detachment of seed from the pedicels by pulling, which is negatively related with seed shattering degree. Totally, thirty randomly chosen spikelets from middle part of the spike of each plant were examined at maturity (4 weeks after heading), and their average BTS values were calculated [69]. The second method is that seed shattering degree (SSD) was measured following the described procedure by Yao [70] with a minor modification. Three spikes from each plant were released at a height of 1 m and freely fell down onto a hard surface. The seed shattering degree was expressed by a percentage (%) of the number of shattered seeds to the total number of seeds. The third method is that seed shattering rate (SSc) was evaluated at maturity. Based on the number of naturally shattered seeds, seed shattering was rated as follows: 1 (> 80% shattering), 2 (60–80% shattering), 3 (40–60% shattering), 4(20–40% shattering), 5 (< 20% shattering). Other seed traits: SL, FN, AL, WS, and SW1 were measured according to the methods described by Zhang et al. [69]. The descriptive statistics of phenotypic data and the correlation analysis between years and traits were calculated by using SPSS software (SPSS, version 19 for Windows, SPSS Inc., Chicago, IL, USA). Heritability (h2) for each trait was estimated based on previously reported method: h2 = σg2/ (σg2 + σgy2/n + σe2/nr), where σg2 is the genotypic variance, σgy2 is the variance caused by the interaction between genotype and year, σe2 is the error variance, n is the number of years, and r is the number of replications [71].

SLAF library construction and high-throughput sequencing

SLAF library construction was performed according to the methods described by Sun et al. [32]. For maximum SLAF-seq efficiency, a pilot experiment was carried out to establish and optimize the conditions required to avoid repetitive SLAFs, achieve an uniform distribution of SLAFs, and obtain optimal SLAF yield. Briefly, HaeIII restriction enzymes (New England Biolabs, NEB, USA) was used to digest the genomic DNA of the two parents and F2 population. A poly-A was added to the 3′ ends of digested fragments. These digested fragments were then ligated with Dual-index sequencing adaptors, and amplified by PCR. The PCR was carried out in reaction solutions containing the diluted restriction-ligation DNA samples, dNTPs, Q5® High-Fidelity DNA polymerase (NEB) and PCR forward primers 5′- AATGATACGGCGACCACCGA-3′ and reverse primer 5′-CAAGCAGAAGACGGCATACG-3′. The PCR products were then purified using Agencourt AMPure XP beads (Beckman Coulter, High Wycombe, UK) and pooled.

The pooled samples were separated with 2% agarose gel. Fragments of 464–494 bp (with barcodes and adaptors) were excised and purified using a QIAquick gel extraction kit (Qiagen, Hilden, Germany). The obtained SLAFs in the quality-tested library were used for paired-end sequencing on an Illumina HiSeq 2500 sequencing platform (Illumina, San Diego, CA, USA). To check the reliability of testing processes, we used the genome of Oryza sativa as a quality control to undergo the same procedures of library construction and sequencing as the E. sibiricus mapping population.

Sequence data analysis and genotyping

SLAF marker identification and genotyping were performed according to the procedures described by Sun et al. [32]. Low-quality reads (with quality score < Q30) were deleted and then the left reads were assigned to the two parents and F2 individuals according to the duplex barcodes. The clean reads were obtained after filtering the barcodes and the terminal 5-bp positions from each read. All paired-end reads (200 bp per read) generated from SLAF-seq raw reads were clustered according to their sequence similarity. Sequences with over 90% identity were grouped into one SLAF locus. As E. sibiricus is self-pollinated species, an F2 population was obtained by self-pollinating the F1 plant of a cross between two fully homozygous parents with genotype aa or bb. Therefore, we only used the SLAF markers with aa×bb segregation pattern for genetic map construction.

Map construction and QTL analysis

HighMap software was used for linkage map construction with four steps: SLAF marker grouping, SLAF marker ordering, genotyping error correction, and genetic map evaluation [72]. The single-linkage clustering algorithm was applied to assign the markers into linkage groups. The modified logarithm of odds (MLOD) score > 5 was set up to partition marker loci into linkage groups (LGs). For comparative genome analysis, we carried out the BLAST search between the mapped SLAF makers and the whole genome sequences of wheat (Triticum aestivum L.) and barley (Hordeum vulgare L.) by using an E-value cutoff of 1e-10 and 90% identity cutoff [20]. For QTL identification, the genotypic data of the mapped markers on the linkage map was integrated with the field phenotypic data of eight seed traits. Logarithm of odds (LOD) scores larger than the 5% cutoff value was used to identify significant loci associated with seed traits. The threshold value was determined through 1000 permutation test according to the composite interval mapping (CIM) method from “qtl” package of R. An interval mapping model with LOD scores of 3.0 for potential QTL was used for QTL detection. MapQTL 6.0 [73] was used to estimate the percentage of phenotypic variation and additive effect explained by a QTL for a trait.

Candidate gene identification

In order to gain an in-depth understanding of the potential functions of these seed yield-related QTLs and identify candidate genes for seed-related traits, the SLAFs within QTL regions were subsequently searched against the Hordeum vulgare and Chinese spring wheat reference genome (; by using the basic local alignment search tool (BLAST). For annotation, the assembled sequences were queried using BLASTX (E-value ≤1e-5) against 7 public databases like the NCBI non-redundant protein sequence (Nr), Gene Ontology (GO), Protein family (Pfam), Cluster of Orthologous Groups (COG), Annotated protein sequence database (Swiss-Prot), and Kyoto Encyclopedia of Genes and Genomes (KEGG), euKaryotic Orthologous Groups (KOG). The expression profiles of candidate gene were obtained from E. sibiricus abscission zone transcriptome data ( The formula log2 (FC) was used to calculate the transcript fold-change, and the false discovery rate (FDR) control method was applied for the correction for multiple tests [27]. Significant differentially expressed transcripts (DETs) between two samples were identified only when an absolute value of the log2 (FC) ≥ 2 and FDR significance score ≤ 0.01 were set as the thresholds. A heatmap was constructed for candidate genes using the Heatmap Illustrator (HemI 1.0) program (Beijing Institude of Genomics, CAS, Beijing, China) [74].

Availability of data and materials

Raw Illumina sequencing data are available in NCBI SRA: SRX6857199-SRX6857400 (, other datasets generated or analysed during this study are included in this published article and its supplementary files.



Abscisic acid


Awn length


Basic local alignment search tool


Breaking tensile strength




Composite interval mapping


Cluster of orthologous groups


Coefficient of variation


Differentially expressed transcripts


Ethylene-responsive factors


Ethylene receptor


Fold change


False discovery rate


Floret number per spike




Gene ontology


Jasmonic acid


Kyoto encyclopedia of genes and genomes


Clusters of orthologous groups for eukaryotic complete genomes


Linkage groups


Logarithm of odds


Marker-assisted selection


The modified logarithm of odds


Next-generation sequencing


RefSeq non-redundant protein


Polymerase chain reaction


Protein family


The percentage of phenotypic variation explained


Quantitative trait loci


Standard deviation


Spike length


Specific-locus amplified fragment sequencing


Single-nucleotide polymorphism


Seed shattering


Seed shattering rate evaluated by five level of classifications


Seed shattering degree assessed by dropping from a height


1000 seed weight


SNP type, Trv means transversion, Tri means transition


Width of seed


Xyloglucan endotransglucosylase


  1. Mott IW, Larson SR, Bushman BS. Simple sequence repeat (SSR) markers for Elymus, Pseudoroegneria and Pascopyrum species (Triticeae: Gramineae). Plant Genet Resour-C. 2011;9(4):489–94.

    Article  CAS  Google Scholar 

  2. Mott IW, Larson SR, Jones TA, Robins JG, Jensen KB, Peel MD. A molecular genetic linkage map identifying the St and H subgenomes of Elymus (Poaceae: Triticeae) wheatgrass. Genome. 2011;54:819–28.

    Article  CAS  PubMed  Google Scholar 

  3. McMillan E, Sun G. Genetic relationships of tetraploid Elymus species and their genomic donor species inferred from polymerase chain reaction-restriction length polymorphism analysis of chloroplast gene regions. Theor Appl Genet. 2004;108:535–42.

    Article  CAS  PubMed  Google Scholar 

  4. Xu DH, Ban T. Phylogenetic and evolutionary relationships between Elymus humidus and other Elymus species based on sequencing of non-coding regions of cpDNA and AFLP of nuclear DNA. Theor Appl Genet. 2004;108:1443–8.

    Article  CAS  PubMed  Google Scholar 

  5. Xie W, Zhang J, Zhao X, Zhang Z, Wang Y. Transcriptome profiling of Elymus sibiricus, an important forage grass in Qinghai-Tibet plateau, reveals novel insights into candidate genes that potentially connected to seed shattering. BMC Plant Biol. 2017;17:78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Dewey DR. Synthetic hybrids of Hordeum bogdanii with Elymus canadensis and Sitanion hystrix. Am J Bot. 1971;58(10):902–8.

    Article  Google Scholar 

  7. Ma X, Chen S, Zhang X, Bai S, Zhang C. Assessment of worldwide genetic diversity of Siberian wild rye (Elymus sibiricus L.) germplasm based on gliadin analysis. Molecules. 2012;17:4424–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Zhao X, Xie W, Zhang J, Zhang Z, Wang Y. Histological characteristics, cell wall hydrolytic enzymes activity and candidate genes expression associated with seed shattering of Elymus sibiricus accessions. Front Plant Sci. 2017;8:606.

    PubMed  PubMed Central  Google Scholar 

  9. Dong Y, Wang YZ. Seed shattering: from models to crops. Front Plant Sci. 2015;6:476.

    PubMed  PubMed Central  Google Scholar 

  10. You MH, Liu JP, Bai SQ, Zhang XQ, Yan JJ. Study on relationship of seed shattering, seed development and yield traits of Elymus sibiricus L. Southwest China J Agric Sci. 2011;24(4):1256–60.

    Google Scholar 

  11. Li C, Zhou A, Sang T. Rice domestication by reducing shattering. Science. 2006;311:1936–9.

    Article  CAS  PubMed  Google Scholar 

  12. Konishi S, Izawa T, Lin SY, Ebana K, Fukuta Y, Sasaki T, Yano M. An SNP caused loss of seed shattering during rice domestication. Science. 2006;312:1392–6.

    Article  CAS  PubMed  Google Scholar 

  13. Ji H, Kim SR, Kim YH, Kim H, Eun MY, Jin ID, Cha YS, Yun DW, Ahn BO, Lee MC, et al. Inactivation of the CTD phosphatase-like gene OsCPL1 enhances the development of the abscission layer and seed shattering in rice. Plant J. 2010;61:96–106.

    Article  CAS  PubMed  Google Scholar 

  14. Zhou Y, Lu D, Li C, Luo J, Zhu BF, Zhu J, Shangguan Y, Wang Z, Sang T, Zhou B, et al. Genetic control of seed shattering in rice by the APETALA2 transcription factor shattering abortion1. Plant Cell. 2012;24:1034–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Yoon J, Cho LH, Kim SL, Choi H, Koh HJ, An G. The BEL1-type homeobox gene SH5 induces seed shattering by enhancing abscission-zone development and inhibiting lignin biosynthesis. Plant J. 2014;79:717–28.

    Article  CAS  PubMed  Google Scholar 

  16. Larson SR, Kellogg EA. Genetic dissection of seed production traits and identification of a major-effect seed retention QTL in hybrid Leymus (Triticeae) wildryes. Crop Sci. 2009;49:29–40.

    Article  CAS  Google Scholar 

  17. Wang YH, Thomas CE, Dean RA. A genetic map of melon (Cucumis melo L.) based on amplified fragment length polymorphism (AFLP) markers. Theor Appl Genet. 1997;95:791–8.

    Article  CAS  Google Scholar 

  18. Alm V, Fang C, Busso CS, Devos KM, Vollan K, Grieg Z, Rognli OA. A linkage map of meadow fescue (Festuca pratensis Huds.) and comparative mapping with other Poaceae species. Theor Appl Genet. 2003;108:25–40.

    Article  CAS  PubMed  Google Scholar 

  19. Zhu WY, Huang L, Chen L, Yang JT, Wu JN, Qu ML, Yao DQ, Guo CL, Lian HL, He HL, et al. A high-density genetic linkage map for cucumber (Cucumis sativus L.): Based on specific length amplified fragment (SLAF) sequencing and QTL analysis of fruit traits in cucumber. Front Plant Sci. 2016;7:437.

    PubMed  PubMed Central  Google Scholar 

  20. Zhao X, Huang L, Zhang X, Wang J, Yan D, Li J, Tang L, Li X, Shi T. Construction of high-density genetic linkage map and identification of flowering-time QTLs in orchardgrass using SSRs and SLAF-seq. Sci Rep. 2016;6:29345.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Price AH, Steele KA, Moore BJ, Barraclough PB, Clark LJ. A combined RFLP and AFLP linkage map of upland rice (Oryza sativa L.) used to identify QTLs for root-penetration ability. Theor Appl Genet. 2000;100:49–56.

    Article  CAS  Google Scholar 

  22. Isobe S, Klimenko I, Ivashuta S, Gau M, Kozlov NN. First RFLP linkage map of red clover (Trifolium pratense L.) based on cDNA probes and its transferability to other red clover germplasm. Theor Appl Genet. 2003;108:105–12.

    Article  CAS  PubMed  Google Scholar 

  23. Hayward MD, Forster JW, Jones JG, Dolstra O, Evans C, McAdam NJ, Hossain KG, Stammers M, Will J, Humphreys MO, et al. Genetic analysis of Lolium. I. Identification of linkage groups and the establishment of a genetic map. Plant Breed. 1998;117:451–5.

    Article  Google Scholar 

  24. Cai HW, Inoue M, Yuyama N, Takahashi W, Hirata M, Sasaki T. Isolation, characterization and mapping of simple sequence repeat markers in zoysiagrass (Zoysia spp.). Theor Appl Genet. 2005;112:158–66.

    Article  CAS  PubMed  Google Scholar 

  25. Xie W, Zhang X, Cai H, Huang L, Peng Y, Ma X. Genetic maps of SSR and SRAP markers in diploid orchardgrass (Dactylis glomerata L.) using the pseudo-testcross strategy. Genome. 2011;54:212–21.

    Article  CAS  PubMed  Google Scholar 

  26. Wei Q, Wang Y, Qin X, Zhang Y, Zhang Z, Wang J, Li J, Lou Q, Chen J. An SNP-based saturated genetic map and QTL analysis of fruit-related traits in cucumber using specific-length amplified fragment (SLAF) sequencing. BMC Genomics. 2014;15:1158.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Davey JW, Hohenlohe PA, Etter PD, Boone JQ, Catchen JM, Blaxter ML. Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nat Rev Genet. 2011;12:499–510.

    Article  CAS  PubMed  Google Scholar 

  28. Lucito R, Nakimura M, West JA, Han Y, Chin K, Jensen K, McCombie R, Gray JW, Wigler M. Genetic analysis using genomic representations. P Natl Acad Sci USA. 1998;95:4487–92.

    Article  CAS  Google Scholar 

  29. Altshuler D, Pollara VJ, Cowles CR, Van Etten WJ, Baldwin J, Linton L, Lander ES. An SNP map of the human genome generated by reduced representation shotgun sequencing. Nature. 2000;407:513–6.

    Article  CAS  PubMed  Google Scholar 

  30. Miller MR, Dunham JP, Amores A, Cresko WA, Johnson EA. Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers. Genome Res. 2007;17:240–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Hegarty M, Yadav R, Lee M, Armstead I, Sanderson R, Scollan N, Powell W, Skøt L. Genotyping by RAD sequencing enables mapping of fatty acid composition traits in perennial ryegrass (Lolium perenne (L.)). Plant Biotechnol J. 2013;11:572–81.

    Article  CAS  PubMed  Google Scholar 

  32. Sun X, Liu D, Zhang X, Li W, Liu H, Hong W, Jiang C, Guan N, Ma C, Zeng H, et al. SLAF-seq: An efficient method of large-scale de novo SNP discovery and genotyping using high-throughput sequencing. PLoS One. 2013;8(3):e58700.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Xu F, Sun X, Chen Y, Huang Y, Tong C, Bao J. Rapid identification of major QTLs associated with rice grain weight and their utilization. PLoS One. 2015;10(3):e0122206.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Yu Y, Zhang X, Li F, Chen X, Zhao Y, Huang L, Zheng H, Xiang J. Development of a high density SNP-based linkage map of pacific white shrimp (Litopenaeus vannamei) using next generation sequencing. Singapore: Plant and Animal Genome Conference; 2014.

    Google Scholar 

  35. Varshney RK, Bertioli DJ, Moretzsohn MC, Vadez V, Krishnamurthy L, Aruna R, Nigam SN, Moss BJ, Seetha K, Ravi K, et al. The first SSR-based genetic linkage map for cultivated groundnut (Arachis hypogaea L.). Theor Appl Genet. 2009;118:729–39.

    Article  CAS  PubMed  Google Scholar 

  36. Saha MC, Mian R, Zwonitzer JC, Chekhovskiy K, Hopkins AA. An SSR- and AFLP-based genetic linkage map of tall fescue (Festuca arundinacea Schreb.). Theor Appl Genet. 2005;110:323–36.

    Article  CAS  PubMed  Google Scholar 

  37. Jones ES, Mahoney NL, Hayward MD, Armstead IP, Jones JG, Humphreys MO, King IP, Kishida T, Yamada T, Balfourier F, et al. An enhanced molecular marker based genetic map of perennial ryegrass (Lolium perenne) reveals comparative relationships with other Poaceae genomes. Genome. 2002;45:282–95.

    Article  CAS  PubMed  Google Scholar 

  38. Song Y, Liu F, Zhu Z, Tan L, Fu Y, Sun C, Cai H. Construction of a simple sequence repeat marker-based genetic linkage map in the autotetraploid forage grass Dactylis glomerata L. Grassl Sci. 2011;57:158–67.

    Article  CAS  Google Scholar 

  39. Xie W, Robins JG, Bushman BS. A genetic linkage map of tetraploid orchardgrass (Dactylis glomerata L.) and quantitative trait loci for heading date. Genome. 2012;55:360–9.

    Article  CAS  PubMed  Google Scholar 

  40. Julier B, Flajoulot S, Barre P, Cardinet G, Santoni S, Huguet T, Huyghe C. Construction of two genetic linkage maps in cultivated tetraploid alfalfa (Medicago sativa) using microsatellite and AFLP markers. BMC Plant Biol. 2003;3:9.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Jones ES, Hughes LJ, Drayton MC, Abberton MT, Michaelson-Yeates TPT, Bowen C, Forster JW. An SSR and AFLP molecular marker-based genetic map of white clover (Trifolium repens L.). Plant Sci. 2003;165:531–9.

    Article  CAS  Google Scholar 

  42. Isobe S, Kölliker R, Hisano H, Sasamoto S, Wada T, Klimenko I, Okumura K, Tabata S. Construction of a consensus linkage map for red clover (Trifolium pratense L.). BMC Plant Biol. 2009;9:57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Huang S, Ding J, Deng D, Tang W, Sun H, Liu D, Zhang L, Niu X, Zhang X, Meng M, et al. Draft genome of the kiwifruit Actinidia chinensis. Nat Commun. 2013;4:2640.

    Article  CAS  PubMed  Google Scholar 

  44. Zhang J, Zhang Q, Cheng T, Yang W, Pan H, Zhong J, Huang L, Liu E. High-density genetic map construction and identification of a locus controlling weeping trait in an ornamental woody plant (Prunus mume Sieb. Et Zucc). DNA Res. 2015;22(3):183–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Zhang Y, Wang L, Xin H, Li D, Ma C, Ding X, Hong W, Zhang X. Construction of a high-density genetic map for sesame based on large scale marker development by specific length amplified fragment (SLAF) sequencing. BMC Plant Biol. 2013;13:141.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Wang W, Zhang T, Zhang G, Wang J, Han K, Wang Y, Zhang Y. Genome-wide association study of antibody level response to NDV and IBV in Jinghai yellow chicken based on SLAF-seq technology. J Appl Genetics. 2015;56:365–73.

    Article  CAS  Google Scholar 

  47. Wang N, Fang L, Xin H, Wang L, Li S. Construction of a high-density genetic map for grape using next generation restriction-site associated DNA sequencing. BMC Plant Biol. 2012;12:148.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Talukder ZI, Gong L, Hulke BS, Pegadaraju V, Song Q, Schultz Q, Qi L. A high-density SNP map of sunflower derived from RAD-sequencing facilitating fine-mapping of the rust resistance gene R 12. PLoS One. 2014;9(7):e98628.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Cai C, Cheng FY, Wu J, Zhong Y, Liu G. The first high-density genetic map construction in tree peony (Paeonia Sect. Moutan) using genotyping by specific-locus amplified fragment sequencing. PLoS One. 2015;10(5):e0128584.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Sun L, Yang W, Zhang Q, Cheng T, Pan H, Xu Z, Zhang J, Chen C. Genome-wide characterization and linkage mapping of simple sequence repeats in mei (Prunus mume Sieb. et Zucc.). PLoS One. 2013;8(3):e59562.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Zhang F, Chen S, Chen F, Fang W, Chen Y, Li F. SRAP-based mapping and QTL detection for inflorescence-related traits in chrysanthemum (Dendranthema morifolium). Mol Breeding. 2011;27:11–23.

    Article  CAS  Google Scholar 

  52. Wang W, Huang S, Liu Y, Fang Z, Yang L, Hua W, Yuan S, Liu S, Sun J, Zhuang M, et al. Construction and analysis of a high-density genetic linkage map in cabbage (Brassica oleracea L. var. capitata). BMC Genomics. 2012;13:523.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Lin Z, Zhang Y, Zhang X, Guo X. A high-density integrative linkage map for Gossypium hirsutum. Euphytica. 2009;166:35–45.

    Article  Google Scholar 

  54. Simons KJ, Fellers JP, Trick HN, Zhang Z, Tai YS, Gill BS, Faris JD. Molecular characterization of the major wheat domestication gene Q. Genetics. 2006;172:547–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Zhang L, Liu D, Wang D, Zhang R, Geng S, Wu L, Li A, Mao L. Over expression of the wheat BEL1-like gene TaqSH1 affects floral organ abscission in Arabidopsis thaliana. J Plant Biol. 2013;56:98–105.

    Article  CAS  Google Scholar 

  56. Pourkheirandish M, Hensel G, Kilian B, Kumlehn J, Sato K, Komatsuda T. Evolution of the grain dispersal system in barley. Cell. 2015;162:527–39.

    Article  CAS  PubMed  Google Scholar 

  57. Taylor JE, Whitelaw CA. Signals in abscission. New Phytol. 2001;151:323–39.

    Article  CAS  Google Scholar 

  58. Sargent JA, Osborne DJ, Dunford SM. Cell separation and its hormonal control during fruit abscission in the Gramineae. J Exp Bot. 1984;35:1663–74.

    CAS  Google Scholar 

  59. Schaller GE, Bleecker AB. Ethylene-binding sites generated in yeast expressing the Arabidopsis ETR1 gene. Science. 1995;270(5243):1809–11.

    Article  CAS  PubMed  Google Scholar 

  60. Zhao Y, Zhang J, Zhang Z, Xie W. Elymus nutans genes for seed shattering and candidate gene-derived EST-SSR markers for germplasm evaluation. BMC Plant Biol. 2019;19:102.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Nakano T, Fujisawa M, Shima Y, Ito Y. The AP2/ERF transcription factor SIERF52 functions in flower pedicel abscission in tomato. J Exp Bot. 2014;65(12):3111–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Ueda J, Miyamoto K, Hashimoto M. Jasmonates promote abscission in bean petiole explants: its relationship to the metabolism of cell wall polysaccharides and cellulase activity. J Plant Growth Regul. 1996;15:189–95.

    Article  CAS  Google Scholar 

  63. Pinyopich A, Ditta GS, Savidge B, Liljegren SJ, Baumann E, Wisman E, Yanofsky MF. Assessing the redundancy of MADS-box genes during carpel and ovule development. Nature. 2003;424:85–8.

    Article  CAS  PubMed  Google Scholar 

  64. Gil-Amado JA, Gomez-Jimenez MC. Transcriptome analysis of mature fruit abscission control in olive. Plant Cell Physiol. 2013;54(2):244–69.

    Article  CAS  PubMed  Google Scholar 

  65. Rushton DL, Tripathi P, Rabara RC, Lin J, Ringler P, Boken AK, Langum TJ, Smidt L, Boomsma DD, Emme NJ, et al. WRKY transcription factors: key components in abscisic acid signalling. Plant Biotechnol J. 2012;10:2–11.

    Article  CAS  PubMed  Google Scholar 

  66. Roberts JA, Elliott KA, Gonzalez-Carranza ZH. Abscission, dehiscence, and other cell separation processes. Annu Rev Plant Biol. 2002;53:131–58.

    Article  CAS  PubMed  Google Scholar 

  67. Zhou HL, He SJ, Cao YR, Chen T, Du BX, Chu CC, Zhang JS, Chen SY. OsGLU1, a putative membrane-bound endo-1, 4-ß-d -glucanase from rice, affects plant internode elongation. Plant Mol Biol. 2006;60:137–51.

    Article  CAS  PubMed  Google Scholar 

  68. Wolf S, Hématy K, Höfte H. Growth control and cell wall signaling in plants. Annu Rev Plant Biol. 2012;63:381–407.

    Article  CAS  PubMed  Google Scholar 

  69. Zhang Z, Zhang J, Zhao X, Xie W, Wang Y. Assessing and broadening genetic diversity of Elymus sibiricus germplasm for the improvement of seed shattering. Molecules. 2016;21:869.

    Article  CAS  PubMed Central  Google Scholar 

  70. Yao N, Wang L, Yan H, Liu Y, Lu BR. Mapping quantitative trait loci (QTL) determining seed-shattering in weedy rice: evolution of seed shattering in weedy rice through de-domestication. Euphytica. 2015;204:513–22.

    Article  Google Scholar 

  71. Yu SB, Li JX, Xu CG, Tan YF, Li XH, Zhang Q. Identification of quantitative trait loci and epistatic interactions for plant height and heading date in rice. Theor Appl Genet. 2002;104:619–25.

    Article  CAS  PubMed  Google Scholar 

  72. Qi Z, Huang L, Zhu R, Xin D, Liu C, Han X, Jiang H, Hong W, Hu G, Zheng H, et al. A high-density genetic map for soybean based on specific length amplified fragment sequencing. PLoS One. 2014;9(8):e104871.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Van Ooijen JW. MapQTL®5, software for the mapping of quantitative trait loci in experimental populations. Wageningen, The Netherlands: Kyazma BV; 2004.

  74. Deng W, Wang Y, Liu Z, Cheng H, Xue Y. HemI: a toolkit for illustrating heatmaps. PLoS One. 2014;9(11):e111988.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


This work was supported by Chinese National Natural Science Foundation (No. 31971751), Chinese National Basic Research Program (2014CB138704), the Fundamental Research Funds for the Central Universities (LZUJBKY-2019-37), Major Science and Technology Projects of Gansu Province and 111 program (B12002). We also thank Mingshu Cao from AgResearch New Zealand for language editing.


This study was financially supported by grants from Chinese National Natural Science Foundation (No. 31971751), Chinese National Basic Research Program (2014CB138704), and the Fundamental Research Funds for the Central Universities (LZUJBKY-2019-37).

Author information

Authors and Affiliations



WX and YW conceived and designed the research, and WX also wrote the paper and analyzed the data. ZZ performed the research, analyzed the data, and wrote this paper, JZ, YZ and NW collected the data, SB revised the paper. All authors read and approved the manuscript.

Corresponding authors

Correspondence to Wengang Xie or Yanrong Wang.

Ethics declarations

Ethics approval and consent to participate

These plant materials used in this study were obtained from Lanzhou University. E. sibiricus are not endangered or protected species, thus, no permissions or licences were required for collecting these samples and conducting this experiment.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

SLAF marker and their locations on each linkage.

Additional file 2: Table S2.

Phenotypic data of parents and F2 population in E. sibiricus.

Additional file 3: Table S3.

Thirty identified candidate genes for seed shattering within six QTL regions based on the functional annotation of barley and wheat genomes.

Additional file 4: Figure S1.

High-density genetic linkage maps of E. sibiricus.

Additional file 5: Figure S2.

Haplotype map of linkage map.

Additional file 6: Figure S3.

Heat map of linkage map.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, Z., Xie, W., Zhang, J. et al. Construction of the first high-density genetic linkage map and identification of seed yield-related QTLs and candidate genes in Elymus sibiricus, an important forage grass in Qinghai-Tibet Plateau. BMC Genomics 20, 861 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: