- Research article
- Open Access
Investigation of genetic relationships within three Miscanthus species using SNP markers identified with SLAF-seq
BMC Genomics volume 23, Article number: 43 (2022)
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 . 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 . 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 .
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) , while second-generation molecular markers include simple sequence repeats (SSRs)  and inter-simple sequence repeats (ISSRs) . However, these markers have several limitations: they are low throughput, inaccurate, time-consuming, labor-intensive, and costly . 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 . 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 . 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 . 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 . 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 .
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 .
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 . For example, Sun et al. used 50,530 SLAFs with 13,291 SNPs to genotype the F1 population of the common carp . SLAF-seq was also used to develop the first high-density genetic maps for several economically important species, including sesame , cucumber , the brown alga Undaria pinnatifida (Phaeophyceae) , wax gourd , watermelon (Citrullus lanatus L.) , tobacco , soybean , peanut , and Salvia miltiorrhiza . 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 , while 5142 polymorphic SLAFs were analyzed to identify a new maize inflorescence meristem mutant . 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 . Furthermore, a high-density genetic linkage map for orchardgrass was developed using 2467 SLAF markers and 43 SSR markers , and in barley the semi-dwarf gene was fine-mapped using molecular markers developed with SLAF-seq . 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  (Fig. 1).
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.
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.
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).
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).
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 . 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 . 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 . 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 . 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 . 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 . 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  and arid conditions . In addition, the recorded distribution of M. lutarioriparius now includes both shallows and wetlands . M. sacchariflorus is considered highly adaptable , 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 . 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.
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′).
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 . 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 (https://www.nature.com/articles/s41467-020-18923-6). Based on the reference genome, suitable restriction enzyme combinations were predicted for the digestion by using DNassist . 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 , 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 . 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 . We used GATK  and SAMtools  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  to determine the population structure of the 50 Miscanthus germplasms. We also performed principal components analysis (PCA) of the germplasms using cluster software . PCAs was performed for linear transformations of variables to create orthogonal axes ordered by the proportion of variance explained . A rooted phylogenetic tree based on our SNP data for these 50 germplasms was constructed using the neighbor-joining (NJ) method  in MEGA 5.0 software . 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: https://www.ncbi.nlm.nih.gov/sra/PRJNA786085.
Amplified fragment length polymorphisms
Bacterial artificial chromosome
Cetyltrimethyl ammonium bromide
Genotyping by sequencing
Inter-simple sequence repeats
Minor allele frequency
Principal components analysis
Restriction site-associated DNA sequencing
Random amplified polymorphic DNA
Restriction fragment length polymorphisms
Specific-locus amplified fragment sequencing
Single nucleotide polymorphisms
Simple sequence repeats
Sang T. Toward the domestication of lignocellulosic energy crops: learning from food crop domestication. J Integr Plant Biol. 2011;53:96–104.
Sang T. China's bioenergy potential. Glob Change Biol Bioenerg. 2011;3:79–90.
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.
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. https://doi.org/10.1088/1748-9326/8/4/044003.
Kirkpatrick J. Construction and analysis of the miscanthus genespace (doctoral dissertation, University of Illinois at Urbana-- Champaign). 2014. https://hdl.handle.net/2142/46925.
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.
De Vega J, Donnison I, Dyer S, Farrar K. Draft genome assembly of the biofuel grass crop Miscanthus sacchariflorus. F1000Research. 2021;10:29.
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.
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.
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. https://doi.org/10.1007/s11105-017-1048-9.
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.
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.
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.
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. https://doi.org/10.1007/s11032-019-1081-5.
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.
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.
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.
John W, Michael M. Fingerprinting genomes using PCR with arbitrary primers. Nucleic Acids Res. 1990;18:7213–8.
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.
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.
Reddy MP, Sarla N, Siddiq EA. Inter simple sequence repeat (ISSR) polymorphism and its application in plant breeding. Euphytica. 2002;128:9–17.
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.
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.
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.
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. https://doi.org/10.1186/s41065-017-0027-3.
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.
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.
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.
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.
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.
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. https://doi.org/10.1186/1471-2229-13-141.
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.
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. https://doi.org/10.1186/s12864-015-2184-y.
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. https://doi.org/10.1186/s12864-015-2220-y.
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.
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) https://doi.org/10.1007/s11032-016-0514-7.
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. https://doi.org/10.1007/s10681-020-02633-7.
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.
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.
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.
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.
Chen S & Renvoize SA. Miscanthus Andersson, Öfvers. Kongl. Vetensk.-Akad. Förh. 12: 165, 1855, (2006). Flora of China 22, 581–583.49.
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. https://doi.org/10.1111/gcbb.12820.
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.
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.
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.
Parkinson J, Blaxter M. Expressed sequence tags: an overview. Methods Mol Biol. 2009;533:1–12.
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.
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.
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. https://doi.org/10.1038/s41598-017-14151-z.
Kim C, Tang H, Paterson AH. Duplication and divergence of grass genomes: integrating the Chloridoids. Trop Plant Biol. 2009;2:51–62.
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.
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.
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. https://doi.org/10.1080/23335777.2021.1879274.
Kaur J, Bhambri P, Gupta OP. Distance based phylogenetic trees with bootstrapping. Int J Comp Appl. 2012;47(24):6–10.
We would like to thank LetPub (www.letpub.com) 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.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
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). https://doi.org/10.1186/s12864-021-08277-8
- SLAF tags
- High-throughput sequencing
- Identification of genetic relationship