Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

A high-density genetic map and QTL analysis of agronomic traits in foxtail millet [Setaria italica (L.) P. Beauv.] using RAD-seq

  • Jun Wang ,

    128wan@163.com (JW); guoerhu2003@163.com (EG)

    Affiliations Millet Research Institute, Shanxi Academy of Agricultural Sciences, Changzhi, Shanxi, China, Shanxi Key Laboratory of Genetic Resources and Breeding in Minor Crops, Changzhi, Shanxi, China

    ⨯
  • Zhilan Wang,

    Affiliations Millet Research Institute, Shanxi Academy of Agricultural Sciences, Changzhi, Shanxi, China, Shanxi Key Laboratory of Genetic Resources and Breeding in Minor Crops, Changzhi, Shanxi, China

    ⨯
  • Xiaofen Du,

    Affiliations Millet Research Institute, Shanxi Academy of Agricultural Sciences, Changzhi, Shanxi, China, Shanxi Key Laboratory of Genetic Resources and Breeding in Minor Crops, Changzhi, Shanxi, China

    ⨯
  • Huiqing Yang,

    Affiliations Millet Research Institute, Shanxi Academy of Agricultural Sciences, Changzhi, Shanxi, China, Shanxi Key Laboratory of Genetic Resources and Breeding in Minor Crops, Changzhi, Shanxi, China

    ⨯
  • Fang Han,

    Affiliation Research Institute of Agriculture Sciences of Yanan, Yanan, Shaanxi, China

    ⨯
  • Yuanhuai Han,

    Affiliation Shanxi Agricultural University, Taigu, Shanxi, China

    ⨯
  • Feng Yuan,

    Affiliations Millet Research Institute, Shanxi Academy of Agricultural Sciences, Changzhi, Shanxi, China, Shanxi Key Laboratory of Genetic Resources and Breeding in Minor Crops, Changzhi, Shanxi, China

    ⨯
  • Linyi Zhang,

    Affiliations Millet Research Institute, Shanxi Academy of Agricultural Sciences, Changzhi, Shanxi, China, Shanxi Key Laboratory of Genetic Resources and Breeding in Minor Crops, Changzhi, Shanxi, China

    ⨯
  • Shuzhong Peng,

    Affiliations Millet Research Institute, Shanxi Academy of Agricultural Sciences, Changzhi, Shanxi, China, Shanxi Key Laboratory of Genetic Resources and Breeding in Minor Crops, Changzhi, Shanxi, China

    ⨯
  • Erhu Guo

    128wan@163.com (JW); guoerhu2003@163.com (EG)

    Affiliations Millet Research Institute, Shanxi Academy of Agricultural Sciences, Changzhi, Shanxi, China, Shanxi Key Laboratory of Genetic Resources and Breeding in Minor Crops, Changzhi, Shanxi, China

    ⨯

Abstract

Foxtail millet (Setaria italica), a very important grain crop in China, has become a new model plant for cereal crops and biofuel grasses. Although its reference genome sequence was released recently, quantitative trait loci (QTLs) controlling complex agronomic traits remains limited. The development of massively parallel genotyping methods and next-generation sequencing technologies provides an excellent opportunity for developing single-nucleotide polymorphisms (SNPs) for linkage map construction and QTL analysis of complex quantitative traits. In this study, a high-throughput and cost-effective RAD-seq approach was employed to generate a high-density genetic map for foxtail millet. A total of 2,668,587 SNP loci were detected according to the reference genome sequence; meanwhile, 9,968 SNP markers were used to genotype 124 F2 progenies derived from the cross between Hongmiaozhangu and Changnong35; a high-density genetic map spanning 1648.8 cM, with an average distance of 0.17 cM between adjacent markers was constructed; 11 major QTLs for eight agronomic traits were identified; five co-dominant DNA markers were developed. These findings will be of value for the identification of candidate genes and marker-assisted selection in foxtail millet.

Citation: Wang J, Wang Z, Du X, Yang H, Han F, Han Y, et al. (2017) A high-density genetic map and QTL analysis of agronomic traits in foxtail millet [Setaria italica (L.) P. Beauv.] using RAD-seq. PLoS ONE 12(6): e0179717. https://doi.org/10.1371/journal.pone.0179717

Editor: Manoj Prasad, National Institute of Plant Genome Research, INDIA

Received: February 8, 2017; Accepted: June 2, 2017; Published: June 23, 2017

Copyright: © 2017 Wang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: This research was supported by the Science and Technology Independent Innovation Ability Upgrading Project of Shanxi Province Academy of Agricultural Sciences (2015ZZCX-09, 2016ZZCX-09), Shanxi Province Youth Fund (2015021143), the Agricultural Science and Technology Innovation Research Project of Shanxi Academy of Agricultural Sciences (ZDSYS1504), and the Chinese Agricultural Research System (CRRS-07).

Competing interests: The authors have declared that no competing interests exist.

Abbreviations: EST, expressed sequence tag; FMID, first main internode diameter; ILP, intron length polymorphic; MAS, marker-assisted selection; MGWP, main grain weight per plant; MPD, main panicle diameter; MPL, main panicle length; MPWP, main panicle weight per plant; NGS, next-generation sequencing; PH, plant height; QTL, quantitative trait locus; RAD-seq, restriction site-associated DNA sequencing; RFLP, restriction fragment length polymorphism; SMID, second main internode diameter; SNP, single-nucleotide polymorphism; TMID, third main internode diameter

Introduction

Foxtail millet (Setaria italica) is one of the oldest cereals in the world, and is thought to have been domesticated from the wild species green foxtail, more than 8,000 years ago in northern China [1]. Foxtail millet has many excellent characteristics, among which C4 photosynthesis, which is the primary mode of carbon capture for some of the world’s most important food, feed, and fuel crops, such as maize, sorghum, sugarcane and switchgrass [2]. In addition, foxtail millet is known for its nutritional value: its grains have high protein, folic acid, vitamin E, carotenoids, and selenium [3–5]. In recent years, foxtail millet has become a valuable model for investigating plant architecture, drought tolerance, and C4 photosynthesis of grain and bioenergy crops because of its small genome size, self-fertilization, and short growth cycle [6–8]. Therefore, it is essential to assess its agronomic traits by developing a genetic linkage map and identifying genes or quantitative trait loci (QTLs).

Genome mapping of foxtail millet using molecular markers started in the 1990s, and a map including 160 restriction fragment length polymorphism (RFLP) markers was constructed in an intervarietal cross [9]. Later, Jia et al. [10, 11] developed 30 expressed sequence tag (EST)-derived simple sequence repeats (SSR) in foxtail millet, and constructed an integrated map with 81 SSR markers and 20 RFLP markers. By exploiting EST sequences, Gupta et al. [12] reported 98 potential intron length polymorphic (ILP) markers. However, these are low-throughput molecular markers that limit the efficiency and accuracy of QTL mapping.

The assembled reference genome of foxtail millet was released in 2012 by two independent groups [13, 14]. The availability of the foxtail millet genome sequence to the public provides an important resource for crop genetics and breeding. On the basis of the reference genome, a large number of markers (SSRs, EST-SSRs, SNPs, InDels, SVs, and TEs) have been developed and utilized [15–21]. For example, Pandey et al. [19] scanned the whole genome sequence of foxtail millet and identified a total of 28,342 microsatellite repeat motifs (SSRs); these markers showed a high percentage (~90%) of cross-genera transferability across millets, including green foxtail, cereals and bioenergy grasses. In previous studies, some traits were analyzed in foxtail millet. Doust et al. [22] identified 25 QTLs for vegetative branching. Meanwhile, Wang et al. [23] detected five QTLs related to plant height, panicle length, panicle weight, and grain weight. Sato et al. [24] analyzed and mapped the stb1 gene, which is responsible for the trait of “spikelet-tipped bristles”. Mauro-Herrera et al. [25] analyzed the differences in flowering time under varying environmental conditions and identified 18 QTLs. Qie et al. [26] detected 18 QTLs for germination and early seeding drought tolerance. Fang et al. [16] described 29 QTLs responsible for 11 agronomic and yield traits. Moreover, SiDWARF2, SiYGL1, and SiAGO1b were fine mapped and cloned using dwarf mutant, AGO1 mutant, and yellow-green leaf mutant by map-based cloning [27–29]. Recently, multiple essential agronomic and quality traits have been studied by identifying QTLs or genes using NGS technologies in foxtail millet. Jia et al. [30] sequenced 916 foxtail millet varieties and identified 512 loci associated with 47 agronomic traits using genome-wide association studies (GWAS). Bai et al. [15] re-sequenced the foxtail millet landrace ‘Shi-Li-Xiang’ (SLX) and finely mapped a waxy gene using newly developed DNA markers. Masumot et al. [31] carried out QTL-seq and rapidly mapped the NEKODE1 gene responsible for tip-branched panicle.

Indeed, most agronomic traits are determined by QTLs [32–34]. It is important to rapidly identify each locus or the major locus of QTLs for efficient crop breeding by marker-assisted selection (MAS). However, whole-genome deep re-sequencing remains cost-prohibitive for sequencing and genotyping of large populations, and is generally unnecessary. Restriction site-associated DNA sequencing (RAD-seq) reduces genome complexity by sequencing only the DNA fragments with restriction sites regardless of length, and is considered a useful tool for SNP discovery and genetic mapping [35–37].

In this study, a high-throughput and cost-effective RAD-seq approach was employed to generate a high-density genetic map for foxtail millet. The characteristics of this genetic map were analyzed and discussed in detail below. Moreover, 11 major QTLs were identified for eight distinct agronomic traits.

Materials and methods

Plant materials

Two foxtail millet cultivars, Hongmiaozhangu (P1) and Changnong35 (P2) were selected as female and male parents, respectively, for the mapping population. Hongmiaozhangu is characterized by low plant height, multiple tillers, and a small panicle. Changnong35 is characterized by high plant height, no tillering, and large panicle. The mapping parents were crossed at the Millet Research Institute (Changzhi, Shanxi) in 2013, and three real F1 hybrids were obtained. F2 seeds were obtained from a self-pollinated F1 individual in 2014. In 2015, the parents and 124 F2 plants obtained from one F1 individual were planted, and plant height (PH, cm), main panicle length (MPL, cm), main panicle diameter (MPD, cm), first main internode diameter (FMID, cm), second main internode diameter (SMID, cm), and third main internode diameter (TMID, cm) were assessed in the mature stage. Main panicle weight per plant (MPWP, g) and main grain weight per plant (MGWP, g) were measured after harvest. Data were analyzed using SPSS 17. Genomic DNA was extracted from young leaf tissues of parental and 124 F2 plants using a Plant DNA Kit (OMEGA, USA, D3485-02).

RAD-seq of the parental lines and F2 population

We employed the RAD protocol described by Baird et al. [35]. The enzymes and restriction fragment sizes were evaluated based on the reference genome sequence (https://www.ncbi.nlm.nih.gov/genome/?term=foxtail+millet). TapI was selected for RAD library construction. The library for Illumina sequencing was constructed from 200 ng of each DNA sample. All library were sequenced using Illumina HiSeq X Ten at Shanghai Major Biological Medicine Technology Co., Ltd.

SNP identification and genotyping

For SNP calling, the Burrows-Wheeler Aligner [38] was applied for sequence alignment between the individual reads and the reference genome sequence, the Genome Analysis ToolKit [39] was used to detect SNP loci, and SAMtools [40] was used to filter out SNP loci. In this study, filtering of SNP loci was based on three criteria: (i) average sequence depth is < 5-fold in parents and < 3-fold in the progeny; (ii) no polymorphism between the parents; (iii) heterozygous in parents.

Development of new DNA markers

Primers were designed according to the flanking sequences (300 bp upstream and 300 bp downstream of the selected SNPs). The primer sequences used for the new DNA markers are listed in S4 Table. PCR was carried out in a 10 μL volume containing 40 ng genomic DNA, 1.0 μL 10× reaction buffer, 0.2 μL 10 mmol L–1 dNTPs, 1.0 μL primer, 1 U rTaq DNA polymerase (TaKaRa, Dalian), and 5.7 μL ddH2O. The PCR was performed by initially denaturing the template DNA at 94°C for 5 min, followed by 35 cycles at 94°C for 30 s, 58°C for 30 s, and 72°C for 30 s, and then terminated by a final extension for 10 min at 72°C. PCR fragments were separated on 8% non-denatured polyacrylamide gel electrophoresis (PAGE) and visualized by silver staining [41].

Linkage map construction

The poorly performing markers were removed before map construction, which excessively missed with more than 30% missing data in the F2 population. Markers with significant segregation distortion (χ2 test, P < 0.05) were excluded from the subsequent linkage map construction. Construction of the linkage map was performed using MSTmap software [42]. The major parameters for loci partition were as follows: distance_function for Haldane; p_value for 0.0000001; no_map_dist for 10; missing_threshold for 0.3; objective_function for Maximum likelihood. A total of 10,016 SNP markers were considered for linkage map construction. LOD = 10 was used to partition the SNP markers into linkage groups. Linkage groups from the same chromosome were merged together, and SNP markers of the same chromosome were again reordered with MSTmap.

QTL analysis

QTL analysis was conducted using CIM (composite interval mapping) of the R/qtl package [43]. LOD thresholds for testing the significance of QTL peaks of each trait were calculated using 1,000 permutations, with a confidence interval of at least 80%. The position and effect of significant QTL was assessed for additive effects and percentage phenotypic variation explained (PVE%) by fitting a model containing all QTL identified for a given trait in R/qtl. The method of naming QTLs was as follows: q plus trait abbreviation and chromosome number; plus -1, -2, etc., when multiple QTLs for one trait were detected on the same chromosome.

Results and discussion

Phenotypic data of eight agronomic traits

Phenotypic data of the eight agronomic traits were analyzed (Table 1 and S1 Fig). A wide range of variation was observed in the eight agronomic traits in the F2 population; the absolute value of the skewness and kurtosis of most traits was less than 1, indicating that these traits had approximately normal distribution. The correlation among the eight traits were calculated. In previous study, Wang et al. [23] found significant associations of PH and MPL with PWP and GWP, respectively. Fang et al. [16] suggested that MPL, MPD, PWP, and GWP are positively correlated with one another. Similar results were reported by Wang et al. [44]. In the present study, all eight traits showed positive correlations with one another, except PH which had non-significant correlations with MPD, FMID, SMID, and TMID (Table 2), corroborating the results of previous studies.

thumbnail
Download:
Table 1. Phenotypic data analyses of eight agronomic traits for 124 F2 individuals.

https://doi.org/10.1371/journal.pone.0179717.t001

thumbnail
Download:
Table 2. Correlation coefficients among agronomic traits in 124 F2 individuals.

https://doi.org/10.1371/journal.pone.0179717.t002

RAD-seq analysis and SNP identification in parental lines and F2 individuals

High-throughput genotyping by sequencing is an option for efficient marker-assisted breeding [45]. Currently, there are many new methods using NGS for identifying genes or QTLs, including GWAS [46], QTL-seq [47], MutMap [48], RAD-seq [35], and SLAF-seq [49], and some have been applied to foxtail millet [15, 30, 31]. However, the utility of RAD-seq has yet not been reported. RAD-seq sequences short DNA fragments with restriction sites digested by restriction endonucleases, regardless of length. It has been applied for SNP identification and linkage map construction in various organisms, including barley, snail, and ryegrass [50–53].

In the present study, a total of 830,740,674 reads were obtained for the parental lines and F2 population. After removal of low-quality reads, 760,378,099 high-quality reads were obtained, representing 91.5% of all reads. The high-quality reads were subsequently mapped to the reference genomic sequence, with a mapping ratio of at least 82% (S1 Table).

A total of 2,668,587 SNP loci were detected, including 2,629,567 located on nine chromosomes and the remaining 39,020 found on scaffolds. Among the 2,629,567 SNP loci, 97.5% (2,564,566) were filtered because of low sequence depth and lack of polymorphism between parents. The remaining 65,001 SNP loci were further screened to remove markers unsuitable for genetic map construction.

Among the 65,001 SNP markers, the numbers of markers on the chromosomes ranged from 4,414 to 10,625; the markers covered at least 98.37% of the physical length of the genome (Table 3). The marker density along each chromosome ranged from 104.82 to 261.63 markers per Mb, averaging 164.64 markers per Mb. The highest marker density (261.63/Mb) was found on chromosome 8, followed by chromosome 7 (213.19/Mb); the lowest marker density (104.82/Mb) was found on chromosome 1.

thumbnail
Download:
Table 3. Number and coverage of SNP markers on the nine chromosomes.

https://doi.org/10.1371/journal.pone.0179717.t003

In this study, we genotyped the parental lines using different letters and determined the segregation patterns of the mapping population. The genotypes of the SNP loci were encoded according to the paternal and maternal genotypes instead of the reference sequence. One locus with a homozygous SNP in the paternal and maternal genotypes would be encoded as aa × bb; if a certain SNP was heterozygous for one or two parents, the genotypes of the markers for the two parents were encoded as, for example, cc × ab or ef × eg. In total, we successfully coded 65,001 polymorphic SNP loci. Furthermore, these SNPs were classified into six segregation patterns (ab × cc, cc × ab, ef × eg, nn × np, lm × ll, and aa × bb). Finally, according to the two parent genotypes, 39,299 markers, which fell into the aa × bb segregation pattern, were used in linkage analysis (Fig 1).

thumbnail
Download:
Fig 1. Number of markers for each segregation pattern.

https://doi.org/10.1371/journal.pone.0179717.g001

High-density genetic linkage map with SNPs

Genetic linkage maps play a major role in clarifying the genetic control of important traits. In particular, high-density molecular markers can be used to quickly map agronomic traits and to identify candidate genes within a region of interest [54]. In this study, after removing incomplete (26,906), significant segregation distortion (24,706), and non-aa × bb (3.373) SNP markers, 10,016 SNP markers were retained for genetic map construction. Finally, a high-density genetic map was constructed containing a total of 9,968 SNP markers. The markers were grouped into nine linkage groups and ordered (Fig 2 and S2 Table). The total genetic distance of the generated map was 1648.8 cM, with an average distance of 0.17 cM between adjacent markers. The largest linkage group was Chr. 8 with 2541 SNP markers and a length of 199.9 cM; the smallest was Chr. 6 with 841 SNP markers and a length of 144.3 cM. Two large gaps (>20 cM) were identified on Chr. 1 (34.23 cM) and Chr. 4 (23.21 cM), respectively. The highest missing data was found on Chr. 6 with 8.42%; the lowest missing data (6.52%) was found on Chr. 7 (Table 4).

thumbnail
Download:
Fig 2. Genetic linkage map and QTLs controlling agronomic traits.

https://doi.org/10.1371/journal.pone.0179717.g002

thumbnail
Download:
Table 4. Characteristics of the high-density genetic map.

https://doi.org/10.1371/journal.pone.0179717.t004

SNP markers are efficient for high-density genetic map construction since they allow high-throughput assessment, compared to RFLP and SSR markers. Currently, the most saturated intervarietal map was constructed by Fang et al. [16]. Compared with this map, the number of mapped loci (1035 SSR makers vs. 9968 SNP markers), marker density (26.25 marker/Mb vs. 164.64 marker/Mb), average distance between adjacent markers (1.27 cM vs. 0.17 cM) and total map length (1318.8 cM vs 1648.8 cM) were significantly increased in the newly constructed SNP maker genetic map. The marker number (9,968 SNPs) in this map was also higher than that of the two consensus maps constructed by Zhang et al. [14] (118 SNPs) and Bennetzen et al. [13] (992 SNPs). The current map provides not only a large number of SNP markers for foxtail millet, but also useful data for QTL analysis, gene fine mapping, and molecular breeding.

The collinearity of each chromosome with the reference genome was also analyzed (Fig 3). The average ratios of genetic-to-physical distance in low- and high-recombination chromosomes were 3.3 cM/Mb (Chr. 9) and 4.98 cM/Mb (Chr. 8), respectively (S3 Table). The 9,968 SNP markers in the genetic map covered 393.53 Mb of the physical length, spanning approximately 76.4% of the foxtail millet genome (∼515 Mb).

thumbnail
Download:
Fig 3. Genetic distance vs. physical distance for 9,968 SNPs in foxtail millet.

https://doi.org/10.1371/journal.pone.0179717.g003

In this study, high levels of collinearity for each chromosome were revealed compared with the reference genome. A relatively low collinearity was observed between Chr. 8 and the reference genome. Non-collinearity could result from multiple factors, including intra- and inter-chromosomal locus duplication, genome rearrangement, transposon-mediated marker transposition, discrepancy in recombination rate among different genomic regions, small mapping population size, compromised marker ordering in the consensus map, missing data, and genotyping errors [55–57].

QTLs for agronomic traits

Agronomic traits play an important role in the breeding of crops. The more QTLs of agronomic traits that are identified, the more they promote breeding by MAS. Recently, Chinese scientists have reported the successful development of new elite varieties in rice by pyramiding major genes that significantly contribute to grain quality and yield from three parents over 5 years, and demonstrated that rational design is a powerful strategy for meeting the challenges of future crop breeding [58]. In foxtail millet, QTLs controlling agronomic traits have been detected in previous studies [16, 22, 23, 30]. Wang et al. [23] and Fang et al. [16] reported QTLs for PH, MPL, MPD, PWP, and GWP, using F2 populations and SSR markers. Jia et al. [30] assessed a natural population of foxtail millet and SNP markers under five different environmental conditions, and identified 512 loci associated with 47 agronomic traits. In the present study, with the exception of MGWP, a total of 11 QTLs were identified for eight agronomic traits using the F2 population and SNP markers (Table 5, Fig 2). The 11 QTLs were mapped to Chr. 1, Chr. 2, Chr. 5, Chr. 7, Chr. 8, and Chr. 9.

thumbnail
Download:
Table 5. QTLs controlling agronomic traits in the Hongmiaozhangu × Changnong35 F2 population.

https://doi.org/10.1371/journal.pone.0179717.t005

MPL had three QTLs, the most prominent of which was designated qMPL8.1, that explained 12.1% of the phenotypic variance. Sixty-five SNP markers covered this interval. The other two QTLs (qMPL1.1 and qMPL1.2) were detected on Chr. 1 with LOD scores of 3.48 and 3.57, and explained 23.3% of the phenotypic variance. Of these, qMPL1.1 and qMPL1.2 have also been detected in the haplotype and SSR maps [16, 30]. Two QTLs were detected for FMID, with the largest effect displayed by qFMID9.1, which explained 15.5% of the phenotypic variance. Forty-seven SNP markers were identified within the chromosomal region of qFMID9.1. Two QTLs were detected for TMID, with the largest effect displayed by qTMID2.1, which explained 11.0% of the phenotypic variance. A total of twenty-two SNP markers were found within the chromosomal region of qTMID2.1. The other QTL, qTMID5.1 with 10.8% of the phenotypic variance was also detected in the haplotype map [30].

PH, MPD, and MPWP had only one QTL each, with LOD scores of 3.36, 3.85, and 3.92, and explained 11.1%, 12.6%, and 12.9% of the phenotypic variance, respectively. qFMID5.1 and qTMID5.1 shared the same chromosomal region, with LOD scores of 3.67 and 3.26, and explained 12.1% and 10.8% of the phenotypic variance, respectively.

With the exception of qMPL1.1 and qTMID5.1, the remaining nine QTLs were identified for the first time, indicating that differences in QTL number and position might be attributed to different mapping populations (genotypes, population sizes, etc.), type and number of markers (i.e., RFLPs, SSRs, and SNPs), and environmental effects.

The additive effects of qPH1.1, qMPL1.1, qMPL1.2, qMPL8.1, qFMID9.1, qSMID9.1, and qTMID2.1 were derived mainly from the female parent (Hongmiaozhangu), whereas those of qMPD7.1, qFMID5.1, qTMID5.1, and qMPWP9.1 stemmed mainly from the male parent (Changnong35).

Newly developed DNA markers

In order to utilize the SNPs identified between Hongmiaozhangu and Changnong35, five SNPs were randomly selected in the intervals of identified QTLs. Based on the flanking sequences of selected SNPs, we designed primers and amplified the target sequences by PCR using genomic DNA of the two parents and their progenies. Finally, five co-dominant DNA markers were developed (Fig 4 and S4 Table). These markers will be useful for gene cloning and molecular breeding of foxtail millet.

thumbnail
Download:
Fig 4. The amplification of newly developed DNA markers in the parents and F2 individuals.

https://doi.org/10.1371/journal.pone.0179717.g004

Mapping population

The construction of F2 populations is general straightforward, and F2 populations provide abundant information suitable for gene or QTL mapping and genetic analysis for many qualitative and quantitative traits. To date, there have been many attempts to use F2 populations directly for QTL analysis in crops. For example, QTLs were obtained that were associated with: crown rust susceptibility in ryegrass; drought-induced flag leaf senescence in wheat; awn, incomplete panicle exertion, and total spikelet number in rice; and aluminum-toxicity tolerance in soybean [32, 59–61]. In this study, the two parental lines have contrasting values across a wide range of agronomic traits, i.e. PH (154.0 and 176.8), MPD (1.8 and 3.1) (Table 1), which are essential for QTL identification. Finally, 11 major QTLs (PVE% > 10) were identified for eight agronomic traits. However, QTLs are sensitive to different environmental factors, such as years, regions, and to sometimes even materials [30, 62]. Therefore, different QTLs were more likely to be obtained in various studies. In particular, in the F2 population, it is important to verify the identified QTLs because of the lack of repeats. In future experiments, we will construct an recombinant inbred line (RIL) population to repeatedly verify the QTLs identified in the present study.

Conclusions

Using a high-throughput and cost-effective RAD-seq approach, we developed a total of 9,968 SNPs to construct a high-density genetic linkage map for foxtail millet, spanning 1648.8 cM, with an average distance of 0.17 cM between adjacent markers. In total, 11 major QTLs were identified for eight agronomic traits, and nine QTLs (qPH1.1, qMPL1.2, qMPL8.1, qFMID9.1, qSMID9.1, qTMID2.1, qFMID5.1, qMPD7.1, and qMPWP9.1) were newly identified. Moreover, five co-dominant markers were developed based on the SNPs between the two parents in the region of the identified QTLs. Our results lay an important foundation for candidate gene identification and MAS breeding of foxtail millet.

Accession number

Raw sequence data obtained in this study have been deposited in the NCBI Sequence Read Archive (SRA) with accession number SRP102319.

Supporting information

S1 Fig. Phenotypic frequency distribution of eight agronomic traits in Hongmiaozhangu × Changnong35 124 F2 individuals.

https://doi.org/10.1371/journal.pone.0179717.s001

(TIF)

S1 Table. The sequencing summary of the parental lines and F2 individuals.

https://doi.org/10.1371/journal.pone.0179717.s002

(XLS)

S2 Table. 9,968 SNP markers in the genetic map.

The genotypes of the 124 individuals were also included.

https://doi.org/10.1371/journal.pone.0179717.s003

(XLSX)

S3 Table. Genetic and physical distances among 9,968 SNP markers in foxtail millet.

https://doi.org/10.1371/journal.pone.0179717.s004

(XLS)

S4 Table. List of primers of newly developed DNA markers.

https://doi.org/10.1371/journal.pone.0179717.s005

(XLSX)

Author Contributions

  1. Conceptualization: JW.
  2. Data curation: ZW HY.
  3. Formal analysis: ZW HY.
  4. Investigation: FY LZ SP.
  5. Methodology: ZW HY.
  6. Project administration: XD ZW.
  7. Resources: FH.
  8. Software: XD.
  9. Supervision: JW YH EG.
  10. Validation: HY.
  11. Visualization: XD.
  12. Writing – original draft: XD.
  13. Writing – review & editing: JW XD.

References

  1. 1. Barton L, Newsome SD, Chen FH, Wang H, Guilderson TP, Bettinger RL. Agricultural origins and the isotopic identity of domestication in northern China. Proceedings of the National Academy of Sciences of the United States of America. 2009; 106(14):5523–8. pmid:19307567.
  2. 2. Brutnell TP, Wang L, Swartwood K, Goldschmidt A, Jackson D, Zhu XG, et al. Setaria viridis: a model for C4 photosynthesis. The Plant cell. 2010; 22(8):2537–44. pmid:20693355.
  3. 3. Liu S, Zhu Z, Li W, Liu F, Li Y, Huang R. Evaluation of selenium and protein content of foxtail millet landraces originated from different ecogical regions of China. Scientia Agricultura Sinica. 2009; 42(11):3812–8.
  4. 4. Liu M, Lu P. Distribution of vitamin E content and its correlation with agronomic traits and carotenoids content in foxtail millet varieties in China. Acta Agronomica Sinica. 2013; 39(3):398.
  5. 5. Shao L, Wang L, Bai W, Liu Y. Evaluation and analysis of folic acid content in millet from different ecological regions in Shanxi province. Scientia Agricultura Sinica. 2014; 47(7):1265–72.
  6. 6. Lata C, Gupta S, Prasad M. Foxtail millet: a model crop for genetic and genomic studies in bioenergy grasses. Critical reviews in biotechnology. 2012; 33(3):328–43. pmid:22985089.
  7. 7. Li P, Brutnell TP. Setaria viridis and Setaria italica, model genetic systems for the Panicoid grasses. Journal of experimental botany. 2011; 62(9):3031–7. pmid:21459768.
  8. 8. Muthamilarasan M, Prasad M. Advances in Setaria genomics for genetic improvement of cereals and bioenergy grasses. Theoretical and applied genetics. 2015; 128(1):1–14. pmid:25239219.
  9. 9. Wang ZM, Devos KM, Liu CJ, Wang RQ, Gale MD. Construction of RFLP-based maps of foxtail millet, Setaria italica (L.) P. Beauv. Theoretical and applied genetics. 1998; 96:31–6.
  10. 10. Jia X, Shi Y, Song Y, Wang G, Wang T, Li Y. Development of EST-SSR in foxtail millet (Setaria italica). Genetic Resources and Crop Evolution. 2007;54(2):233–6.
  11. 11. Jia X, Zhang Z, Liu Y, Zhang C, Shi Y, Song Y, et al. Development and genetic mapping of SSR markers in foxtail millet [Setaria italica (L.) P. Beauv.]. Theoretical and applied genetics. 2009; 118(4):821–9. pmid:19139840.
  12. 12. Gupta S, Kumari K, Das J, Lata C, Puranik S, Prasad M. Development and utilization of novel intron length polymorphic markers in foxtail millet [Setaria italica (L.) P. Beauv.]. Genome. 2011; 54:586–602. pmid:21751869
  13. 13. Bennetzen JL, Schmutz J, Wang H, Percifield R, Hawkins J, Pontaroli AC, et al. Reference genome sequence of the model plant Setaria. Nature biotechnology. 2012; 30(6):555–61. pmid:22580951.
  14. 14. Zhang G, Liu X, Quan Z, Cheng S, Xu X, Pan S, et al. Genome sequence of foxtail millet (Setaria italica) provides insights into grass evolution and biofuel potential. Nature biotechnology. 2012; 30(6):549–54. pmid:22580950.
  15. 15. Bai H, Cao Y, Quan J, Dong L, Li Z, Zhu Y, et al. Identifying the genome-wide sequence variations and developing new molecular markers for genetics research by re-sequencing a landrace cultivar of foxtail millet. PloS one. 2013; 8(9):e73514. pmid:24039970
  16. 16. Fang X, Dong K, Wang X, Liu T, He J, Ren R, et al. A high density genetic map and QTL for agronomic and yield traits in foxtail millet [Setaria italica (L.) P. Beauv]. BMC genomics. 2016; 17:336. pmid:27146360.
  17. 17. Kumari K, Muthamilarasan M, Misra G, Gupta S, Subramanian A, Parida SK, et al. Development of eSSR-markers in Setaria italica and their applicability in studying genetic diversity, cross-transferability and comparative mapping in millet and non-millet species. PloS one. 2013; 8(6):e67742. pmid:23805325.
  18. 18. Obidiegwu ON, Obidiegwu JE, Parzies H. Development of SSR for foxtail millet (Setaria italica (L.) P. Beauv.) and its utility in genetic discrimination of a core set. Genes & Genomics. 2013; 35(5):609–15.
  19. 19. Pandey G, Misra G, Kumari K, Gupta S, Parida SK, Chattopadhyay D, et al. Genome-wide development and use of microsatellite markers for large-scale genotyping applications in foxtail millet [Setaria italica (L.)]. DNA research. 2013; 20(2):197–207. pmid:23382459.
  20. 20. Yadav CB, Bonthala VS, Muthamilarasan M, Pandey G, Khan Y, Prasad M. Genome-wide development of transposable elements-based markers in foxtail millet and construction of an integrated database. DNA research. 2015; 22(1):79–90. pmid:25428892.
  21. 21. Zhang S, Tang C, Zhao Q, Li J, Yang L, Qie L, et al. Development of highly polymorphic simple sequence repeat markers using genome-wide microsatellite variant analysis in foxtail millet [Setaria italica (L.) P. Beauv.]. BMC genomics. 2014; 15(1):78.
  22. 22. Doust AN, Devos KM, Gadberry MD, Gale MD, Kellogg EA. Genetic control of branching in foxtail millet. Proceedings of the National Academy of Sciences of the United States of America. 2004; 101(24):9045–50. pmid:15184666.
  23. 23. Wang X, Diao X, Wang J, Wang C, Wang G, Hao X, et al. Construction of genetic map and QTL analysis of some main agronomic traits in millet. Journal of Plant Genetic Resoures. 2013; 14(5):871–8.
  24. 24. Sato K, Mukainari Y, Naito K, Fukunaga K. Construction of a foxtail millet linkage map and mapping of spikelet-tipped bristles 1(stb1) by using transposon display markers and simple sequence repeat markers with genome sequence information. Molecular Breeding. 2013; 31(3):675–84.
  25. 25. Mauro-Herrera M, Wang X, Barbier H, Brutnell TP, Devos KM, Doust AN. Genetic control and comparative genomic analysis of flowering time in Setaria (Poaceae). G3 Genes Genomes Genetics. 2013; 3(2):283–95. pmid:23390604
  26. 26. Qie L, Jia G, Zhang W, Schnable J, Shang Z, Li W, et al. Mapping of quantitative trait locus (QTLs) that contribute to germination and early seedling drought tolerance in the interspecific cross Setaria italica × Setaria viridis. PloS one. 2014; 9(7):e101868. pmid:25033201.
  27. 27. Xue C, Zhi H, Fang X, Liu X, Tang S, Chai Y, et al. Characterization and fine mapping of SiDWARF2 (D2) in foxtail Millet. Crop Science. 2016;56(1):95.
  28. 28. Li W, Tang S, Zhang S, Shan J, Tang C, Chen Q, et al. Gene mapping and functional analysis of the novel leaf color gene SiYGL1 in foxtail millet [Setaria italica (L.) P. Beauv]. Physiologia plantarum. 2016; 157(1):24–37. pmid:26559175.
  29. 29. Liu X, Tang S, Jia G, Schnable JC, Su H, Tang C, et al. The C-terminal motif of SiAGO1b is required for the regulation of growth, development and stress responses in foxtail millet (Setaria italica (L.)P. Beauv). Journal of experimental botany. 2016; 67(11):3237–49. pmid:27045099.
  30. 30. Jia G, Huang X, Zhi H, Zhao Y, Zhao Q, Li W, et al. A haplotype map of genomic variations and genome-wide association studies of agronomic traits in foxtail millet (Setaria italica). Nature genetics. 2013; 45(8):957–61. pmid:23793027.
  31. 31. Masumoto H, Takagi H, Mukainari Y, Terauchi R, Fukunaga K. Genetic analysis of NEKODE1 gene involved in panicle branching of foxtail millet, Setaria italica (L.) P. Beauv., and mapping by using QTL-seq. Molecular Breeding. 2016; 36(5).
  32. 32. Herlina L, Sobir S, Trijatmiko KR. Identification of quantitative trait loci (QTL) for awn, incomplete panicle exertion and total spikelet number in an F2 population derived from a backcross inbred line, Bio-148, and the recurrent parent, IR64. Makara Journal of Science. 2016; 20(1).
  33. 33. Zhai H, Feng Z, Li J, Liu X, Xiao S, Ni Z, et al. QTL analysis of spike morphological traits and plant height in winter wheat (Triticum aestivum L.) using a high-density SNP and SSR-based linkage map. Frontiers in plant science. 2016; 7:1617. pmid:27872629.
  34. 34. Li H, Yang Q, Fan N, Zhang M, Zhai H, Ni Z, et al. Quantitative trait locus analysis of heterosis for plant height and ear height in an elite maize hybrid zhengdan 958 by design III. BMC genetics. 2017; 18(1):36. pmid:28415964.
  35. 35. Baird NA, Etter PD, Atwood TS, Currey MC, Shiver AL, Lewis ZA, et al. Rapid SNP discovery and genetic mapping using sequenced RAD markers. PloS one. 2008; 3(10):e3376. pmid:18852878
  36. 36. Davey JW, Blaxter ML. RADSeq: next-generation population genetics. Briefings in functional genomics. 2010; 9(5–6):416–23. pmid:21266344.
  37. 37. 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 Research. 2007; 17(2):240–8. pmid:17189378
  38. 38. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009; 25(14):1754–60. pmid:19451168.
  39. 39. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome research. 2010; 20(9):1297–303. pmid:20644199.
  40. 40. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment / map format and SAMtools. Bioinformatics. 2009; 25(16):2078–9. pmid:19505943.
  41. 41. Marklund S, Chaudhary R, Marklund L, Sandberg K, Andersson L. Extensive mtDNA diversity in horses revealed by PCR-SSCP analysis. Animal genetics. 1995; 26(3):193–6. pmid:7793689.
  42. 42. Wu Y, Bhat PR, Close TJ, Lonardi S. Efficient and accurate construction of genetic linkage maps from the minimum spanning tree of a graph. PLoS genetics. 2008; 4(10):e1000212. pmid:18846212.
  43. 43. Broman KW, Wu H, Sen S, Churchill GA. R/qtl: QTL mapping in experimental crosses. Bioinformatics. 2003; 19(7):889–90. pmid:12724300
  44. 44. Wang H, Jia G, Zhi H, Wen Q, Dong J, Chen L, et al. Phenotypic diversity evaluations of foxtail millet core collections. Acta Agronomica Sinica. 2016; 42(1):19.
  45. 45. Cabezas JA, Ibanez J, Lijavetzky D, Velez D, Bravo G, Rodriguez V, et al. A 48 SNP set for grapevine cultivar identification. BMC plant biology. 2011; 11:153. pmid:22060012.
  46. 46. Huang X, Wei X, Sang T, Zhao Q, Feng Q, Zhao Y, et al. Genome-wide association studies of 14 agronomic traits in rice landraces. Nature genetics. 2010; 42(11):961–7. pmid:20972439.
  47. 47. Takagi H, Abe A, Yoshida K, Kosugi S, Natsume S, Mitsuoka C, et al. QTL-seq: rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations. The Plant journal: for cell and molecular biology. 2013; 74(1):174–83. pmid:23289725.
  48. 48. Abe A, Kosugi S, Yoshida K, Natsume S, Takagi H, Kanzaki H, et al. Genome sequencing reveals agronomically important loci in rice using MutMap. Nature biotechnology. 2012; 30(2):174–8. pmid:22267009.
  49. 49. Sun X, Liu D, Zhang X, Li W, Liu H, Hong W, et al. SLAF-seq: An efficient method of large scale SNP discovery and genotyping using high throughput sequencing. PloS one. 2013; 8(3):e58700. pmid:23527008
  50. 50. Chutimanitsakun Y, Nipper RW, Cuesta-Marcos A, Cistue L, Corey A, Filichkina T, et al. Construction and application for QTL analysis of a Restriction Site Associated DNA (RAD) linkage map in barley. BMC genomics. 2011; 12:4. pmid:21205322.
  51. 51. Liu MM, Davey JW, Banerjee R, Han J, Yang F, Aboobaker A, et al. Fine mapping of the pond snail left-right asymmetry (chirality) locus using RAD-Seq and fibre-FISH. PloS one. 2013; 8(8):e71067. pmid:23951082.
  52. 52. Pfender WF, Saha MC, Johnson EA, Slabaugh MB. Mapping with RAD (restriction-site associated DNA) markers to rapidly identify QTL for stem rust resistance in Lolium perenne. Theoretical and applied genetics. 2011; 122(8):1467–80. pmid:21344184.
  53. 53. Richards PM, Liu MM, Lowe N, Davey JW, Blaxter ML, Davison A. RAD-Seq derived markers flank the shell colour and banding loci of the Cepaea nemoralis supergene. Molecular ecology. 2013; 22(11):3077–89. pmid:23496771.
  54. 54. Bolger ME, Weisshaar B, Scholz U, Stein N, Usadel B, Mayer KF. Plant genome sequencing—applications for crop improvement. Current opinion in biotechnology. 2014; 26:31–7. pmid:24679255.
  55. 55. Hackett CA, Broadfoot LB. Effects of genotyping errors, missing values and segregation distortion in molecular marker data on the construction of linkage maps. Heredity. 2003; 90(1):33–8. pmid:12522423.
  56. 56. Hudson CJ, Kullan ARK, Freeman JS, Faria DA, Grattapaglia D, Kilian A, et al. High synteny and colinearity among Eucalyptus genomes revealed by high-density comparative genetic mapping. Tree Genetics & Genomes. 2011; 8(2):339–52.
  57. 57. Raman H, Raman R, Kilian A, Detering F, Long Y, Edwards D, et al. A consensus map of rapeseed (Brassica napus L.) based on diversity array technology markers: applications in genetic dissection of qualitative and quantitative traits. BMC genomics. 2013; 14:277. pmid:23617817
  58. 58. Zeng D, Tian Z, Rao Y, Dong G, Yang Y, Huang L, et al. Rational design of high-yield and superior-quality rice. Nature plants. 2017; 3:17031. pmid:28319055.
  59. 59. Tomaszewski C, Byrne SL, Foito A, Kildea S, Kopecký D, Doležel J, et al. Genetic linkage mapping in an F2 perennial ryegrass population using DArT markers. Plant Breeding. 2012; 131(2):345–9.
  60. 60. Barakat MN, Wahba LE, Milad SI. Molecular mapping of QTLs for wheat flag leaf senescence under water-stress. Biologia Plantarum. 2012; 57(1):79–84.
  61. 61. Tasma IM, Warsun A, Asadi A. Development and characterization of F2 population for molecular mapping of aluminum-toxicity tolerant QTL in soybean. Jurnal AgroBiogen. 2016; 4(1):1–8.
  62. 62. Mauro-Herrera M, Doust AN. Development and genetic control of plant architecture and biomass in the Panicoid grass, setaria. PloS one. 2016; 11(3):e0151346. pmid:26985990.