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

Investigation of genetic relationships within three Miscanthus species using SNP markers identified with SLAF-seq



Miscanthus, which is a leading dedicated-energy grass in Europe and in parts of Asia, is expected to play a key role in the development of the future bioeconomy. However, due to its complex genetic background, it is difficult to investigate phylogenetic relationships in this genus. Here, we investigated 50 Miscanthus germplasms: 1 female parent (M. lutarioriparius), 30 candidate male parents (M. lutarioriparius, M. sinensis, and M. sacchariflorus), and 19 offspring. We used high-throughput Specific-Locus Amplified Fragment sequencing (SLAF-seq) to identify informative single nucleotide polymorphisms (SNPs) in all germplasms.


We identified 257,889 SLAF tags, of which 87,162 were polymorphic. Each tag was 264–364 bp long. The obtained 724,773 population SNPs were used to investigate genetic relationships within three species of Miscanthus. We constructed a phylogenetic tree of the 50 germplasms using the obtained SNPs and grouped them into two clades: one clade comprised of M. sinensis alone and the other one included the offspring, M. lutarioriparius, and M. sacchariflorus. Genetic cluster analysis had revealed that M. lutarioriparius germplasm C3 was the most likely male parent of the offspring.


As a high-throughput sequencing method, SLAF-seq can be used to identify informative SNPs in Miscanthus germplasms and to rapidly characterize genetic relationships within this genus. Our results will support the development of breeding programs with the focus on utilizing Miscanthus cultivars with elite biomass- or fiber-production potential for the developing bioeconomy.


Miscanthus, a high-biomass-yielding perennial C4 grass, has great potential utility for biobased industry [1,2,3,4]. This is a heterogeneous gramineous plant that hybridizes interspecifically and hence has a complex genetic background. Recent molecular studies have revealed that the Miscanthus genome is unusually large (2.07–2.6 GB) [5,6,7,8,9]. In Miscanthus, inter- and intraspecific hybridization is common which is why this genus is characterized by rich genetic diversity and heterosis [6, 9]. The genetic diversity has been exploited to develop Miscanthus hybrids, which can deliver higher biomass yields and exhibit better adaptability to diverse climatic conditions than their parent species. The improved biomass yields of hybrids are attributed to efficient rates of photosynthesis [3]. On the other hand, this high degree of genetic diversity also increases the complexity of interspecific relationships in the Miscanthus and subsequently makes it challenging to carry out the genetic evolution analyses of this genus. Thus, it is difficult to mine functional genes in Miscanthus which can seriously limit the utilization potential of Miscanthus [5]. Molecular markers would be useful for further investigations of Miscanthus plants because such markers have been widely used in studies of genetics, molecular population genetics, species formation, evolutionary and phylogenetic relationships [10,11,12], and molecular taxonomy [6].

First-generation molecular markers cover restriction fragment length polymorphisms (RFLPs) [7, 8], random amplified polymorphic DNA (RAPD) [9, 10], and amplified fragment length polymorphisms (AFLPs) [11], while second-generation molecular markers include simple sequence repeats (SSRs) [12] and inter-simple sequence repeats (ISSRs) [13]. However, these markers have several limitations: they are low throughput, inaccurate, time-consuming, labor-intensive, and costly [3]. These drawbacks have motivated the development of third-generation polymorphic molecular markers, which are named as SNPs. These markers are generally widely distributed throughout the whole genome [14]. SNP markers are amenable to large-scale automated monitoring and have been instrumental in various crop breeding programs, such as the construction of genetic maps, the DNA fingerprinting of germplasm resources, the detection of molecular biodiversity, and the analysis of linkage disequilibrium [15]. This continuous development of molecular marker technology has accelerated functional gene identification and characterization in other crops, which has led to the development of varieties with improved functional traits [5]. Thus, these techniques might be useful for molecular genetic research in Miscanthus.

Although genotyping-by-sequencing (GBS) and restriction site-associated DNA sequencing (RAD-seq) have been used extensively in Miscanthus, there are still some difficulties and challenges associated with the application of these techniques in Miscanthus [16]. For example, one obstacle in the widespread use of GBS is the difficulty in carrying out the associated bioinformatics analysis, which is typically hampered by a large number of erroneous SNP interferences that are not easy to diagnose or correct [17].

To overcome these challenges, we aimed to develop and identify SNP markers for Miscanthus using SLAF-seq techniques. Moreover, reduction in the genomic complexity using specific digestion, develop markers via the high-throughput sequencing of representative libraries, and determine phylogenetic relationships using genotyping is also part of this study.

SLAF-seq uses bioinformatics methods to systematically analyze known genome sequences as well as of related species, bacterial artificial chromosome (BAC) sequences, or Fosmid sequences [18,19,20,21,22]. SLAF-seq techniques differ in several ways from GBS or RAD-seq techniques. The major differences cover: a) SLAF-seq identifies one tag about every 10 K; b) SLAF tags are uniformly distributed, which means that important chromosome segments are not missed; c) SLAF-seq avoids repetitive sequences, which makes it a cost effective technique. Moreover, such SLAF-seq utilizes deep sequencing to ensure genotyping accuracy, a pre-designed representation scheme to optimize marker efficiency and a double-barcode system for large populations [23].

Due to its efficient identification of SNP markers, SLAF-seq has been widely used in variety of crops [24,25,26,27,28] for the development of specific molecular markers and genetic maps [29]. For example, Sun et al. used 50,530 SLAFs with 13,291 SNPs to genotype the F1 population of the common carp [30]. SLAF-seq was also used to develop the first high-density genetic maps for several economically important species, including sesame [31], cucumber [32], the brown alga Undaria pinnatifida (Phaeophyceae) [33], wax gourd [34], watermelon (Citrullus lanatus L.) [35], tobacco [36], soybean [37], peanut [37], and Salvia miltiorrhiza [38]. In addition, an increasing number of studies for the Gramineae family have been performed using SLAF-seq [25, 27, 28, 39,40,41]. For example, SLAF-seq was used to develop the first 7E-chromosome-specific molecular markers for Thinopyrum elongatum [39], while 5142 polymorphic SLAFs were analyzed to identify a new maize inflorescence meristem mutant [40]. Zhang et al. used 69,325 high-quality SLAFs, of which 26,248 were polymorphic, to develop sufficient markers for a segregating Agropyron F1 population [28]. Furthermore, a high-density genetic linkage map for orchardgrass was developed using 2467 SLAF markers and 43 SSR markers [25], and in barley the semi-dwarf gene was fine-mapped using molecular markers developed with SLAF-seq [27]. The successful application of SLAF-seq in other species provides reference material for the development of SNPs in this study.

Genome-wide SLAF markers and SNPs for the three Miscanthus species are generated by using SLAF-seq as a part of current study. In addition, phylogenetic relationships are estimated with these species based on the generated SNPs. The genome-wide markers for Miscanthus identified in this study will lead towards the utilization of its genetic resources to develop molecular marker-assisted Miscanthus breeding programs.


Evaluation of the digestive enzymes

Enzyme digestion was performed according to set selection principles which included: a) the proportion of restriction fragment located in the repeat sequence was low; b) the fragments were evenly distributed in the control genome; c) the length of restriction fragment was appropriate for the experimental system; d) the number of SLAF tags was consistent with the expected number of tags. The pair-end digestion efficiency of EcoRV + ScaI for the control genome (Nipponbare) was 90.87%, indicating that this enzyme combination was suitable.

Analysis of the SLAF-seq and SNP data

We obtained 57.8 Mb clean sequence reads based on SLAF library construction and high-throughput sequencing. The average sequencing depth per sample was 11.76x for the female parent (sample A12), whereas 15.47x for the male parents (samples A1–A11, B1–B10, and C1–C9), and 7.85x for the progeny (samples D1–D19). The average GC content across all sequences was 41.39%. Across all sequences, the average number of bases with a quality score ≥ 30 (Q30) was 93.66%. In parallel, we obtained 7.17 Mb reads by sequencing the rice genome, which indicated that the experimental database was accurate.

Using the obtained clean sequence reads, we developed 257,889 SLAF tags, of which 87,162 were polymorphic. We also generated a map showing the distribution of SLAF tags across the Miscanthus chromosomes based on the chromosome-level of reference genome data for M. lutarioriparius [8] (Fig. 1).

Fig. 1
figure 1

Distribution of SLAF tags on Miscanthus chromosomes. Chromosome length is shown on the abscissa: each band represents a chromosome, and the genome is divided into 1 Mb segments. Each segment is colored according to the number of SLAF tags: the darker areas are those where SLAF tags are concentrated

In total, 724,773 highly consistent population SNPs with integrity of > 0.8 and minor allele frequency (MAF) of > 0.05 were identified across all samples in this study. A map showing the distribution of these SNPs across the Miscanthus chromosomes is presented in Fig. 2.

Fig. 2
figure 2

Map showing the distribution of SNPs across Miscanthus chromosomes. Chromosome length is shown on the abscissa: each band represents a chromosome, and the genome is divided into 1 Mb segments. Each segment is colored according to the number of SNPs: the darker areas are those where SNPs are concentrated

Genetic diversity of three Miscanthus species

Genetic diversity analyses can provide information about the origin and composition of individual lineages. In this study, we analyzed the genetic diversity of three Miscanthus species based on SNPs. First, population structure was analyzed under the assumption that the number of clusters (K) were ranged 1–10. At a minimum value of ΔK, there were four clusters, suggesting that all of our samples may have originated from four primitive ancestors. Cluster graphs showing K values of 1–10 of 50 samples are shown in Fig. 3. The 50 samples clustered into four groups: the A group included A1–A11; the B group included B1–B10, C1, C2, and C4–C9; the C group included C3, D8, D9, D11–D13, D15, and D17–D19; and the D group included A12, D1–D7, D10, D14, and D16.

Fig. 3
figure 3

Cluster graphs of K = 1–10 with 50 samples. Each color in the figure to the left represents a cluster, with each row representing a cluster value. In the figure to the right, ΔK is the value corresponding to each K value. ΔK was minimized when K = 4

PCA analysis indicated that the M. sinensis (A1-A11) germplasms were distinct from all other groups. The offspring (D1-D19), the female parent A12 (M. lutarioriparius), and sample C3 (M. lutarioriparius) formed a loose cluster (Fig. 4). The remaining M. lutarioriparius and most of the M. sacchariflorus formed a cluster, as did B8 and B10 (M. sacchariflorus, Fig. 4).

Fig. 4
figure 4

Principal component analysis (PCA) of the 50 Miscanthus germplasms. Each point represents a sample, and points of different color and shape correspond to different sample groups

The phylogenetic analysis indicated that the investigated 50 accessions fell into three distinct groups: 1) group containing the accessions from M. sinensis (A1-A11); 2) containing the accessions from A12, C3, and all offspring; 3) this group containing the accessions from M. lutarioriparius and M. sacchariflorus except B8 and B10 (Fig. 5).

Fig. 5
figure 5

Neighbor-joining phylogenetic analysis of all 50 Miscanthus germplasms. The M. lutarioriparius genome was used as the reference. Each node represents a sample. Scale bar corresponds to genetic distance

The results of phylogenetic analysis showed that the genetic variation within some populations was greater than between populations. Furthermore, population structure analysis, PCA, and phylogenetic analysis indicated that C3 was the male parent of the offspring.


Optimization of molecular markers based on SLAF-seq and SNPs in Miscanthus

Currently, supply of sufficient quantities of sustainably produced biomass with optimal quality characteristics is a major challenge in the development of biobased industry. Thus, genetic improvements to deliver high biomass yield with required quality traits can be a way forward. Biomass quality may be substantially improved by the development of genetic markers associated with quality traits [31]. SNPs, which are more abundant in the genome than any other molecular markers, are particularly useful for the analyses of genetic diversity and population structure [33]. In this study, we used SLAF-seq to efficiently identify SNP markers. Compared with other methods, such as GBS and RAD-seq, SLAF-seq is more accurate, faster, and less expensive. Moreover, SLAF-seq also reduces genome complexity [34]. Here, we obtained 257,889 SLAF tags and 724,773 SNPs, which is greater than the number of SNPs previously obtained in the Miscanthus genome using RAD-seq [35]. In addition, SLAF and SNP markers were evenly distributed on each chromosome based on the chromosome-level of the reference genome of M. lutarioriparius (Figs. 1, 2). These polymorphic molecular markers are highly discriminatory and can be used for genetic map construction and gene mapping in Miscanthus.

Phylogenic analysis based on Miscanthus SNPs

The heterozygosity and polyploidy that have accumulated in Miscanthus genomes over their long evolutionary history make it more difficult to sequence complete genomes in this genus [5]. Fortunately, the genomes of M. sinensis, M. lutarioriparius, M. sacchariflorus, and M. floridulus have been sequenced through the joint efforts of many researchers [6,7,8]. With these reference genomes, the accuracy of Miscanthus phylogenetic analyses can be further improved through SNPs.

In this study, the M. lutarioriparius was selected as the reference genome. The use of the M. lutarioriparius reference genome enabled us not only to cluster the 50 samples, but also to draw a rooted phylogenetic tree (Figs. 5, 6). The phylogenetic tree in combination with the morphological characteristics indicated that the offspring had been produced by the intraspecific hybridization of M. lutarioriparius, which was also consistent with previous studies [8, 11, 15]. The phylogenetic analyses revealed that M. lutarioriparius is adjacent to M. sacchariflorus and had shown that the coefficient of intraspecific genetic variation in these two species was high.

It had been believed for long that the distribution of M. lutarioriparius is extremely narrow, consisting of only the middle and lower shallows of the Yangtze River in Southern China [42]. However, in recent years, latest studies had revealed that M. lutarioriparius is extremely adaptable with diverse geographic distribution and can even thrive in marginal areas such as saline alkali soils [43] and arid conditions [44]. In addition, the recorded distribution of M. lutarioriparius now includes both shallows and wetlands [45]. M. sacchariflorus is considered highly adaptable [7], thus this diverse geographic distribution of M. lutarioriparius can be explained through the genetic similarity between M. lutarioriparius and M. sacchariflorus. The latest studies had shown that M. lutarioriparius and M. sinensis diverged very recently [8]. This argument can be supported by the interspecific genetic distance between M. lutarioriparius and M. sinensis, which was even lower than that within M. lutarioriparius.


We obtained 724,773 SNPs using SLAF-seq technology. We successfully identified the paternal parent and obtained an intraspecific hybrid polyploid population of M. lutarioriparius. Despite the high similarity between the genomes, M. lutarioriparius is morphologically distinct from M. sacchariflorus. The polyploid M. lutarioriparius stands out as an excellent strain, which produces more biomass and is highly adaptable than any other Miscanthus species, including the diploid M. sacchariflorus. Based on the SNPs obtained for this population as a part of this study, high-density molecular marker linkage maps could be constructed. Such maps would be valuable genetic resource for the development of miscanthus based bioeconomy. Our results will also support future genetic improvements in biomass yield as well as quality traits in Miscanthus species.


Plant materials

A total of 50 Miscanthus specimens were used in this study: the maternal plant Miscanthus lutarioriparius (A12; tetraploid, 2n = 4x), the paternal candidate M. sinensis (A1–A11; diploid, 2n = 2x), the paternal candidate M. saccariflorus (B1–B10; diploid, 2n = 2x), the paternal candidate M. lutarioriparius (C1–C9; diploid, 2n = 2x), and 19 progenies (D1–D19; triploid, 2n = 3x). The 19 progenies were generated via open pollination from A12. A12, which was bagged during flowering to test its self-pollination seed-setting rate. The measured self-pollinating seed-setting rate was 0, indicating that A12 was highly self-incompatible. Other Miscanthus genotypes that flowered simultaneously with A12 in the Germplasm Nursery were selected as paternal candidates. The collection of maternal and paternal candidate materials were collected from across China (Table 1). All the plants were grown in the Miscanthus Germplasm Nursery at Hunan Agricultural University, Changsha, Hunan, China (latitude 28°11′, longitude 113°04′).

Table 1 Origins of the 31 pure species Miscanthus germplasms used in this study. The female parent is sample A12; the remaining samples are the candidate male parents. All collection sites are in China

DNA extraction

Fresh leaves from 50 Miscanthus individuals were collected and frozen in liquid nitrogen before being manually ground into a fine powder. Total genomic DNA was extracted following a modified cetyltrimethyl ammonium bromide (CTAB) method [46]. The leaves of Miscanthus are rich in phenols, therefore CTAB extraction solution was supplemented with 2% poly-N-vinylpyrrolidone (PVP) and 1% β-mercaptoethanol to purify the Miscanthus DNA. The concentration and quality of the extracted DNA was detected using 0.8% agarose gel electrophoresis and an ND-1000 spectrophotometer (Nano Drop), respectively.

Enzyme digestion design

To identify the most appropriate enzymes for genomic digestion, we selected the M. lutarioriparius genome as the reference genome ( Based on the reference genome, suitable restriction enzyme combinations were predicted for the digestion by using DNassist [47]. To assess the efficiency of the predicted enzymes, the genome of Japanese rice (Oryza sativa L. subsp. japonica) was selected as a control.

SLAF library construction and high-throughput sequencing

Using the identified restriction enzyme combination (EcoRV + ScaI), the total genomic DNA of each sample was digested. After adding an A-tail to the 3′ end of each digested fragment and ligating the fragment to the Dual-index sequencing adapter [48], each fragment was amplified through PCR. The amplicons were purified, mixed, and cut. The digested fragments were chosen as target segments. The libraries were selected and sequenced using an Illumina HiSeq TM 2500 platform (Biomarker Technologies Corporation), with a read length of 2 × 100 bp. The obtained sequenced reads were mapped based on the M. lutarioriparius reference genome for subsequent mutation analysis. To assess the accuracy of SLAF library construction, the Japanese rice genome was again used as a control.

Development of SLAF tags and SNP markers

The Dual-index tags were used to classify the raw sequencing data by sample. Sequences reads from the same locus were grouped using similarity clustering [49]. In general, only high-depth fragments were selected in each cluster group, whereas low-depth segments were removed. Here, we first calculated the SLAF tags for each sample independently, and then all single-sample SLAF tags were clustered to derive population-wide SLAF tags. The positions of clean reads on the reference genome were compared, the sequencing depth of each sample was counted, and variations were detected. Each sequence was aligned with the reference sequence using bwa [47]. We used GATK [48] and SAMtools [49] to identify SNPs. The SNPs, identified by both methods were considered reliable. Of these reliable SNPs, those with integrity > 0.8 and MAF > 0.05 were considered highly consistent and were used for subsequent analyses.

Genetic relationships among samples

We used admixture software [46] to determine the population structure of the 50 Miscanthus germplasms. We also performed principal components analysis (PCA) of the germplasms using cluster software [50]. PCAs was performed for linear transformations of variables to create orthogonal axes ordered by the proportion of variance explained [51]. A rooted phylogenetic tree based on our SNP data for these 50 germplasms was constructed using the neighbor-joining (NJ) method [52] in MEGA 5.0 software [53]. The phylogenetic tree based on our samples was constructed using the Kimura 2-parameter model with 1000 bootstrap replicates [54, 55].

Availability of data and materials

The Miscanthus SLAF tags and SNPs supporting the findings of this work are available within the paper and its Supplementary Information files. The raw date of Miscanthus sequence generated and analyzed during the current study are available from NCBI SRA data PRJNA786085, and the SRA records will be accessible with the following link:



Amplified fragment length polymorphisms


Bacterial artificial chromosome


Cetyltrimethyl ammonium bromide


Genotyping by sequencing


Inter-simple sequence repeats


Minor allele frequency




Principal components analysis


Poly-N- vinylpyrrolidone


Restriction site-associated DNA sequencing


Random amplified polymorphic DNA


Restriction fragment length polymorphisms


Specific-locus amplified fragment sequencing


Single nucleotide polymorphisms


Simple sequence repeats


  1. Sang T. Toward the domestication of lignocellulosic energy crops: learning from food crop domestication. J Integr Plant Biol. 2011;53:96–104.

    PubMed  Google Scholar 

  2. Sang T. China's bioenergy potential. Glob Change Biol Bioenerg. 2011;3:79–90.

    Google Scholar 

  3. Clifton-Brown J, Harfouche A, Casler MD, Dylan Jones H, Macalpine WJ, Murphy-Bokern D, et al. Breeding progress and preparedness for mass-scale deployment of perennial lignocellulosic biomass crops switchgrass, miscanthus, willow and poplar. Glob Change Biol Bioenerg. 2019;11:118–51.

    Google Scholar 

  4. Liu W, Sang T. Potential productivity of the Miscanthus energy crop in the loess plateau of China under climate change. Environ Res Lett. 2013;8(4)4003.

  5. Kirkpatrick J. Construction and analysis of the miscanthus genespace (doctoral dissertation, University of Illinois at Urbana-- Champaign). 2014.

  6. Mitros T, Session AM, James BT, Wu GA, Belaffif MB, Clark LV, et al. Genome biology of the paleotetraploid perennial biomass crop Miscanthus. Nat Commun. 2020;11:5442.

    CAS  PubMed  PubMed Central  Google Scholar 

  7. De Vega J, Donnison I, Dyer S, Farrar K. Draft genome assembly of the biofuel grass crop Miscanthus sacchariflorus. F1000Research. 2021;10:29.

    PubMed  PubMed Central  Google Scholar 

  8. Miao J, Feng Q, Li Y, Zhao Q, Zhou C, Lu H, et al. Chromosome-scale assembly and analysis of biomass crop Miscanthus lutarioriparius genome. Nat Commun. 2021;12:2458.

    CAS  PubMed  PubMed Central  Google Scholar 

  9. Zhang G, Ge C, Xu P, Wang S, Cheng S, Han Y, et al. The reference genome of Miscanthus floridulus illuminates the evolution of Saccharinae. Nature Plants. 2021;7:608–18.

    CAS  PubMed  PubMed Central  Google Scholar 

  10. Ge CX, Liu XM, Liu SM, Xu J, Li HF, Cui TT, et al. Miscanthus sp.: Genetic Diversity and Phylogeny in China. Plant Mol Biol Report. 2017;35:600–10.

    Article  CAS  Google Scholar 

  11. Li SS, Zhou HF, Chen WL, Yan J, Cai Z, Wei RX, et al. Population genetics and evolutionary history of Miscanthus species in China. J Syst Evol. 2019;57(5):530–42.

    Google Scholar 

  12. Tang YM, Xiao L, Iqbal Y, Liao JF, Xiao LQ, Yi ZL, et al. Molecular cytogenetic characterization and phylogenetic analysis of four Miscanthus species (Poaceae). Comp Cytogenet. 2019;13(3):211–30.

    PubMed  PubMed Central  Google Scholar 

  13. Olatoye MO, Clark LV, Labonte NR, Dong H, Dwiyanti MS, Anzoua KG, et al. Training Population Optimization for Genomic Selection in Miscanthus. G3 (Bethesda). 2020;10:2465–76.

    CAS  Google Scholar 

  14. Olatoye MO, Clark LV, Wang JP, Yang XP, Yamada T, Sacks EJ, et al. Evaluation of genomic selection and marker-assisted selection in Miscanthus and energycane. Mol Breed. 2019;39:171.

  15. Graner A, Jahoor A, Schondelmaier J, Siedler H, Pillen K, Fischbeck G, et al. Construction of an RFLP map of barley. Theor Appl Genet. 1991;83:250–6.

    CAS  PubMed  Google Scholar 

  16. Paterson AH, Lander ES, Hewitt JD, Peterson S, Lincoln SE, Tanksley SD. Resolution of quantitative traits into Mendelian factors by using a complete linkage map of restriction fragment length polymorphisms. Nature. 1988;335:721–6.

    CAS  PubMed  Google Scholar 

  17. Williams JG, Kubelik AR, Livak KJ, Rafalski JA, Tingey SV. DNA polymorphisms amplified by arbitrary primers are useful as genetic markers. Nucleic Acids Res. 1990;18:6531–5.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. John W, Michael M. Fingerprinting genomes using PCR with arbitrary primers. Nucleic Acids Res. 1990;18:7213–8.

    Google Scholar 

  19. Vos P, Hogers R, Bleeker M, Reijans M, Van De Lee T, Hornes M, et al. AFLP: a new technique for DNA fingerprinting. Nucleic Acids Res. 1995;23:4407–14.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Litt M, Luty JA. A hypervariable microsatellite revealed by in vitro amplification of a dinucleotide repeat within the cardiac muscle actin gene. Am J Hum Genet. 1989;44(3):397–401 PMID: 2563634; PMCID: PMC1715430.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. Reddy MP, Sarla N, Siddiq EA. Inter simple sequence repeat (ISSR) polymorphism and its application in plant breeding. Euphytica. 2002;128:9–17.

    Google Scholar 

  22. Clark LV, Brummer JE, Glowacka K, Hall MC, Heo K, Peng JH, et al. A footprint of past climate change on the diversity and population structure of Miscanthus sinensis. Ann Bot. 2014;114:97–107.

    PubMed  PubMed Central  Google Scholar 

  23. Ge C, Ai X, Jia S, Yang Y, Che L, Yi Z, et al. Interspecific genetic maps in Miscanthus floridulus and M. sacchariflorus accelerate detection of QTLs associated with plant height and inflorescence. Mol Gen Genomics. 2019;294:35–45.

    CAS  Google Scholar 

  24. Spindel J, Wright M, Chen C, Cobb J, Gage J, Harrington S, et al. Bridging the genotyping gap: using genotyping by sequencing (GBS) to add high-density SNP markers and new value to traditional bi-parental mapping and breeding populations. Theor Appl Genet. 2013;126:2699–716.

    CAS  PubMed  Google Scholar 

  25. Zeng B, Yan HD, Liu XC, Zang WJ, Zhang AL, Zhou SF, et al. Genome-wide association study of rust traits in orchardgrass using SLAF-seq technology. Hereditas. 2017;154:5.

  26. Liu LQ, Luo QL, Li HW, Li B, Li ZS, Zheng Q. Physical mapping of the blue-grained gene from Thinopyrum ponticum chromosome 4Ag and development of blue-grain-related molecular markers and a FISH probe based on SLAF-seq technology. Theor Appl Genet. 2018;131:2359–70.

    PubMed  Google Scholar 

  27. Jia Q, Tan C, Wang J, Zhang X-Q, Zhu J, Luo H, et al. Marker development using SLAF-seq and whole-genome shotgun strategy to fine-map the semi-dwarf gene ari-e in barley. BMC Genomics. 2016;17:911.

    PubMed  PubMed Central  Google Scholar 

  28. Zhang Y, Zhang JP, Huang L, Gao AN, Zhang J, Yang XM, et al. A high-density genetic map for P genome of Agropyron Gaertn. based on specific-locus amplified fragment sequencing (SLAF-seq). Planta. 2015;242:1335–47.

    CAS  PubMed  Google Scholar 

  29. Elshire RJ, Glaubitz JC, Sun Q, Poland JA, Kawamoto K, Buckler ES, et al. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS One. 2011;6:e19379.

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Sun X, Liu D, Zhang X, Li W, Liu H, Hong W, et al. SLAF-seq: an efficient method of large-scale de novo SNP discovery and genotyping using high-throughput sequencing. Plos One. 2013;8:e58700.

    CAS  PubMed  PubMed Central  Google Scholar 

  31. Zhang Y, Wang L, Xin H, Li D, Ma C, Ding X, et al. 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 

  32. Wei QZ, Wang YZ, Qin XD, Zhang YX, Zhang ZT, Wang J, et al. 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.

    PubMed  PubMed Central  Google Scholar 

  33. Shan T, Pang S, Li J, Li X, Su L. Construction of a high-density genetic map and mapping of a sex-linked locus for the brown alga Undaria pinnatifida (Phaeophyceae) based on large scale marker development by specific length amplified fragment (SLAF) sequencing. BMC Genomics. 2015;16:902.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Jiang B, Liu WR, Xie DS, Peng QW, He XM, Lin YE, et al. High-density genetic map construction and gene mapping of pericarp color in wax gourd using specific-locus amplified fragment (SLAF) sequencing. BMC Genomics. 2015;16:1035.

  35. Shang JL, Li N, Li NN, Xu YY, Ma SW, Wang JM. Construction of a high-density genetic map for watermelon (Citrullus lanatus L.) based on large-scale SNP discovery by specific length amplified fragment sequencing (SLAF-seq). Scientia Horticulturae. 2016;203:38–46.

    CAS  Google Scholar 

  36. Gong DP, Huang L, Xu XH, Wang CY, Ren M, Wang CK, et al. Construction of a high-density SNP genetic map in fluecured tobacco based on SLAF-seq. Mol Breed. 2016;36(100)

  37. Dong ZM, Chen L, Li Z, Liu NX, Zhang SC, Liu J, et al. Identification and molecular mapping of the semi-dwarf locus (sdf-1) in soybean by SLAF-seq method. Euphytica. 2020;216:103.

  38. Liu T, Guo L, Pan Y, Zhao Q, Wang J, Song Z. Construction of the first high-density genetic linkage map of Salvia miltiorrhiza using specific length amplified fragment (SLAF) sequencing. Sci Rep. 2016;6:24070.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. Chen S, Huang Z, Dai Y, Qin S, Gao Y, Zhang L, et al. The development of 7E chromosome-specific molecular markers for Thinopyrum elongatum based on SLAF-seq technology. Plos One. 2013;8:e65122.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Xia C, Chen LL, Rong TZ, Li R, Xiang Y, Wang P, et al. Identification of a new maize inflorescence meristem mutant and association analysis using SLAF-seq method. Euphytica. 2015;202:35–44.

    Google Scholar 

  41. Zhao X, Huang L, Zhang X, Wang J, Yan D, Li JM, et al. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. Chen S & Renvoize SA. Miscanthus Andersson, Öfvers. Kongl. Vetensk.-Akad. Förh. 12: 165, 1855, (2006). Flora of China 22, 581–583.49.

  43. Zheng C, Xue S, Xiao L, Iqbal Y, Sun GR, Duan MJ, et al. “Two-steps” seed-derived plugs as an effective propagation method for the establishment of Miscanthus in saline–alkaline soil. GCB Bioenergy. 2021;3.

  44. Yan J, Zhu CY, Liu W, Luo F, Mi J, Ren YJ, et al. High photosynthetic rate and water use efficiency of Miscanthus lutarioriparius characterize an energy crop in the semiarid temperate region. GCB Bioenergy. 2015;7:207–18.

    CAS  Google Scholar 

  45. Xue S, Guo MQ, Iqbal Y, Liao JF, Yang S, Xiao L, et al. Mapping current distribution and genetic diversity of the native Miscanthus lutarioriparius across China. Renew Sustain Energy Rev. 2020;134:110386.

    Google Scholar 

  46. Xia YN, Xu J, Duan JY, Liu QB, Huang HM, Yi ZL, et al. Transgenic Miscanthus lutarioriparius that co-expresses the cry 2Aa# and Bar genes. Can J Plant Sci. 2019;99:841–51.

    CAS  Google Scholar 

  47. Parkinson J, Blaxter M. Expressed sequence tags: an overview. Methods Mol Biol. 2009;533:1–12.

    CAS  PubMed  Google Scholar 

  48. Kanazin V, Talbert H, See D, Decamp P, Nevo E, Blake T. Discovery and assay of single nucleotide polymorphisms in barley (Hordeum vulgare). Plant Mol Biol. 2002;48:529–37.

    CAS  PubMed  Google Scholar 

  49. Jones JC, Fan SH, Franchini P, Schartl M, Meyer A. The evolutionary history of Xiphophorus fish and their sexually selected sword: a genome-wide approach using restriction site-associated DNA sequencing. Mol Ecol. 2013;22:2986–3001.

    CAS  PubMed  Google Scholar 

  50. Sheng J, Zheng X, Wang J, Zeng X, Zhou F, Jin S, et al. Transcriptomics and proteomics reveal genetic and biological basis of superior biomass crop Miscanthus. Sci Rep. 2017;7:13777.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Kim C, Tang H, Paterson AH. Duplication and divergence of grass genomes: integrating the Chloridoids. Trop Plant Biol. 2009;2:51–62.

    Google Scholar 

  52. Kim C, Zhang D, Auckland SA, Rainville LK, Jakob K, Kronmiller B, et al. SSR-based genetic maps of Miscanthus sinensis and M. sacchariflorus, and their comparison to sorghum. Theor Appl Genet. 2012;124:1325–38.

    CAS  PubMed  Google Scholar 

  53. Phillips SM, Dubery IA, Van Heerden H. Molecular characterisation of two homoeologous elicitor-responsive lipin genes in cotton. Mol Gen Genomics. 2013;288:519–33.

    CAS  Google Scholar 

  54. Mahadani AK, Awasthi S, Sanyal G, Bhattacharjee P, Pippal S. Indel-K2P: a modified Kimura 2 parameters (K2P) model to incorporate insertion and deletion (Indel) information in phylogenetic analysis. Cyber Phys Syst. 2021.

  55. Kaur J, Bhambri P, Gupta OP. Distance based phylogenetic trees with bootstrapping. Int J Comp Appl. 2012;47(24):6–10.

    Google Scholar 

Download references


We would like to thank LetPub ( for providing linguistic assistance during the preparation of this manuscript.


This work was supported by The National Natural Science Foundation of China [31471557, 31871693]. The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations



Author contributions

This research was designed by ZYC, HHH and ZLY; and was performed by YCH with assistance from YLS. ZYC and YCH wrote the paper with help from the other authors. All authors read and approved the final manuscript.

Authors’ information

At the time of this work, YCH and YLS were postgraduate students in the College of Bioscience and Biotechnology at Hunan Agricultural University, Hunan, China. ZYC, HHH, and YI are Assistant Professors, and ZLY is a Professor in the College of Bioscience and Biotechnology at Hunan Agricultural University. All authors are also members of Hunan Engineering Laboratory of Miscanthus Ecological Applications, Hunan, China.

Corresponding authors

Correspondence to Zhiyong Chen, Hongmei Huang or Zili Yi.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

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.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, Z., He, Y., Iqbal, Y. et al. Investigation of genetic relationships within three Miscanthus species using SNP markers identified with SLAF-seq. BMC Genomics 23, 43 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: