QTL Mapping of Fiber-Related Traits Based on a High-Density Genetic Map in Flax (Linum usitatissimum L.)

A genetic map is an important and valuable tool for quantitative trait locus (QTL) mapping, marker-assisted selection (MAS)-based breeding, and reference-assisted chromosome assembly. In this study, 112 F₂ plants from a cross between Linum usitatissimum L. “DIANE” and “NY17” and parent plants were subjected to high-throughput sequencing and specific-locus amplified fragment (SLAF) library construction. After preprocessing, 61.64 Gb of raw data containing 253.71 Mb paired-end reads, each 101 bp in length, were obtained. A total of 192,797 SLAFs were identified, of which 23,115 were polymorphic, with a polymorphism rate of 11.99%. Finally, 2,339 SLAFs were organized into a linkage map consisting of 15 linkage groups (LGs). The total length of the genetic map was 1483.25 centimorgans (cM) and the average distance between adjacent markers was 0.63 cM. Combined with flax chromosome-scale pseudomolecules, 12 QTLs associating with 6 flax fiber-related traits were mapped on the chromosomal scaffolds. This high-density genetic map of flax should serve as a foundation for flax fine QTL mapping, draft genome assembly, and MAS-guided breeding. Ultimately, the genomic regions identified in this research could potentially be valuable for improving flax fiber cultivars, as well as for identification of candidate genes involved in flax fiber formation processes.

SIGNIFICANCE STATEMENT

A high-density genetic map of flax was constructed, and QTLs were identified on the sequence scaffolds to be interrelated with fiber-related traits. The results of this study will not only provide a platform for gene/QTL fine mapping, map-based gene isolation, and molecular breeding for flax, but also provide a reference to help position sequence scaffolds on the physical map and assist in the process of assembling the flax genome sequence.

Introduction

Flax (Linum usitatissimum L., 2n = 2x = 30) is a dual-purpose crop that is grown for its stem fiber (for making linen) and oilseed (for making linseed oil). Although fiber and oilseed flax belong to the same species, they exhibit morphologically distinct characteristics, including height and fiber composition, among others. Linen flax fiber is an important raw material for manufacture of textile products. Flax fiber is needed in abundance in China, the top linen exporter worldwide, with Heilongjiang contributing 85% of China's total linen production (FAO, 2014¹). Meanwhile, with regard to linseed flax, China is the third largest source after India and Canada, with the United States coming in fourth. Gansu, Inner Mongolia, and Xinjiang are the three provinces/regions in China where the majority of linseed flax is grown.

Optimal flax cultivars for textile production can yield flax of high quality, high yield, high fiber content, and other desirable characteristics. However, compared to other crops, little systematic in-depth research has been focused on improving flax cultivars for use with traditional breeding methods. Marker-assisted selection (MAS) would enhance the effectiveness of existing breeding methods to improve flax stem fiber, but currently insufficient molecular resources are available for achieving that goal. Indeed, in spite of the economic importance of flax, until now breeding programs utilizing MAS resources have only been used to construct genetic linkage maps for various crops such as rice (Orjuela et al., 2010; Chen et al., 2012), wheat (Gardner et al., 2016; Holtz et al., 2016), corn (Foiada et al., 2015; Chen et al., 2016), and soybean (Li et al., 2014, 2016; Qi et al., 2014; Song et al., 2015). Fortunately, MAS resources of other crops can provide a valuable foundation of knowledge for future molecular biological studies of flax.

Until recently, flax improvement has relied mainly on conventional breeding methods based upon limited germplasm resources. However, most important agronomic characteristics, such as stem fiber content, are quantitative innate traits that are regulated by micro-effect polygenes and environmental factors. Because increasing the fiber content is an ultimate goal in flax genetics and breeding research, identification of quantitative trait loci (QTL) would be advantageous for MAS and map-based cloning. Unfortunately, a linkage map with suitable marker density as a prerequisite for QTL detection does not yet exist (Thoday, 1961; Paterson et al., 1998). Indeed, only four limited flax linkage maps and one integrated low-density map have been published (Spielmeyer et al., 1998; Oh et al., 2000; Cloutier et al., 2011, 2012; Yi et al., 2017). Moreover, integration of these maps to generate a draft flax genome would require additional information before it could be used to effectively support MAS-based breeding. Therefore, the goal of this study was to construct a suitably dense flax genetic linkage map composed of previously published flax genome assemblies (Wang Z. et al., 2012) in combination with mapped flax fiber trait QTLs discovered in this study. These genetic resources should lay a foundation for further mining of flax genes important for optimal fiber generation for use in MAS-based breeding of flax.

Results and Analysis

Analysis of SLAF Sequencing Data and Genotyping

A total of 61.64 Gb of raw sequence data was generated using Illumina sequencing after SLAF library construction, which produced 253.71 Mb of paired-end reads, each 101 bp in length. The GC (guanine-cytosine) content was 39.18% and the Q30 ratio (quality scores of at least 30, indicating a 1% chance of error) was 86.64%. For the paternal inbred line (“DIANA”), 9,868,940 reads and 99,043 SLAFs were generated, with an average coverage of 99.64-fold for each SLAF. For the maternal line (“NY17”), the number of reads producing 97,769 SLAFs was 10,899,810 and the average coverage for each SLAF marker was 111.49-fold (Table 1). For the analysis of the F₂ mapping population, 1,061,346–3,035,345 reads were generated for the development of 67,607–90,504 SLAF markers for each plant, with marker coverage ranging from 3.64- to 9.16-fold for an average coverage of 5.48-fold (Figure 1A). The average count of SLAFs per individual plant was 80,998 (Figure 1B).

TABLE 1

Table 1. Statistical tables of sequencing data.

FIGURE 1

Figure 1. Coverage and number of markers for each of the F2 individuals. The x-axes in both (A,B) indicate each of the F₂ individuals; the y-axes indicate the marker coverage (A) and the number of markers developed for each F₂ plant (B).

After correcting or discarding low-depth SLAF tags, 192,797 high-quality SLAFs were identified, of which 23,115 were polymorphic, for a polymorphism rate of 11.99% (Table 2). The remaining 169,682 SLAFs were non-polymorphic or repetitive. After the parental lines were assigned different alphabetic letter designations as genotypes to determine segregation patterns, 10,434 polymorphic SLAFs were successfully encoded and grouped into eight segregation patterns (ab × cd, ef × eg, ab × cc, cc × ab, hk × hk, lm × ll, nn × np, and aa × bb) following a genotype encoding rule (Figure 2). Since the two parents (“NY17” and “DIANA”) are homozygous inbred lines with genotypes of aa and bb, respectively, 7,417 markers that fell within the aa × bb segregation pattern were used for performance of linkage analysis.

TABLE 2

Table 2. Discovery of SLAF markers.

FIGURE 2

Figure 2. Number of polymorphic SLAF markers for eight segregation patterns. The x-axis indicates eight segregation patterns of polymorphic SLAF markers; the y-axis indicates the number of markers.

Basic Characteristics of the Genetic Map

After removing incomplete and significant segregation distortion markers, 7,417 SLAFs were retained for the final genetic linkage map construction. The average integrity of the mapped markers was 100%, predicting potential generation of a relatively high quality genetic map. The mLOD values were calculated between pairs of tags following a standard procedure as previously described (Vision et al., 2000). After removal of each tag with a LOD score < 10 when paired with all other tags, 2,339 markers were localized to the map for an overall genome coverage percentage estimate of 31.54% (Table S1). HighMap software (Liu et al., 2014) was used to arrange markers into linkage groups and genetic distances were estimated between adjacent markers. Subsequently, a genetic linkage map of overall length 1,483.25 cM was obtained (Figure 3), with an average distance between adjacent markers of 0.63 cM. The basic characteristics of all linkage groups (LGs) obtained are shown in Table 3. In the overall linkage map, the largest linkage group (LG15) contained 196 SLAFs, while the smallest linkage group (LG12, 34.05 cM) contained 39 SLAFs. LG11 contained the maximum number of markers (765), while LG12 contained the minimum marker number (39). On average, each LG contained 156 SLAF markers. The genetic length of the 15 linkage groups ranged from 34.05 cM (LG12) to 157.28 cM (LG15), with an average interval between adjacent markers ranging from 0.20 cM (LG11) to 1.42 cM (LG2). The “Max Gap,” which reflects the largest degree of linkage between markers, was 16.52 cM and was located in LG3. A total of 5,106 SNP loci were identified among the 2,339 mapped SLAF markers (Table 4) with Tri/Trv values ranging from 1.40 to 2.67. LG11 contained the highest number of SNP loci (1,737), while LG12 contained the lowest number of SNPs (66).

FIGURE 3

Figure 3. Distribution of SLAF markers on the 15 linkage groups of flax. A black bar indicates a SLAF marker. The x-axis represents the linkage group number and the y-axis indicates the genetic distance (cM) within each linkage group.

TABLE 3

Table 3. Basic characteristics of the 15 linkage groups of flax.

TABLE 4

Table 4. Distribution of SNP loci on the 15 linkage groups of flax.

Quality Evaluation of the Genetic Map

The quality of the flax genetic map was evaluated by analyzing the integrity of mapped markers. The average integrity of each individual marker was 99.33% (Figure 4), the average depth of offspring was 7.72, and the parental depth was five times greater than the depth of offspring (Table 5), all of which suggest genotyping accuracy (Table S2).

FIGURE 4

Figure 4. The integrity distribution map of all individuals. The x-axis represents the 112 individuals and y-axis represents the complete degree of mapped markers.

TABLE 5

Table 5. Statistics of the mapped marker depth.

Haplotype maps, which reflect chromosomal-exchange events among genomes within a population, were developed for the 112 offspring and parental controls using the 2,339 SLAF markers (Presentation S1). Dual-exchange sites might result from two scenarios: (1) a recombinant hotspot region within the genome and (2) genotyping error caused by sequencing within a linkage group. Consequently, a higher proportion of dual-exchange events correlated with greater success of genotyping and more effective ordering of markers. In this study, the percentage of dual-exchanges ranged from 0.01% (LG1) to 0.25% (LG12), while the percentage of missing data for each LG ranged from 0.01% (LG15) to 2.92% (LG9) (Table 6). Ultimately, most recombination blocks were clearly defined and the LGs were uniformly distributed, suggesting that genetic mapping was of high quality.

TABLE 6

Table 6. Double exchange and missing ratio of the mapped markers.

A genetic map essentially reflects multipoint recombination analysis, with closer distances between adjacent markers reflecting smaller observed recombination rates. To analyze recombination relationships between markers, we determined the potential layout of mapped markers. The quality of the genetic map was also evaluated using heat maps which directly depicted recombination relationships among markers for all fifteen linkage groups (Presentation S2). Each cell of the heat map represents a recombination rate between two adjacent markers whereby the rate level was depicted using different colors ranging from yellow to purple (yellow indicating a lower recombination rate; purple indicating a higher rate). In this way, heat maps were generated for each LG using recombination scores for pair-wise comparative analyses of the 2,339 markers (Presentation S2). The resulting heat maps indicate that SLAF markers in most LGs were well ordered.

QTL Mapping for Related Traits

Phenotypic data for F₂ families are presented in Table 7 and the frequency distribution of all measured traits is shown in Figure 5. Plant height, stem length, and stem yield all strictly followed normal distributions, with kurtosis and skewing values close to zero. Meanwhile, the frequency distributions of seed yield and fiber content exhibited approximately normal distributions, but with higher kurtoses (although only the left peaks were characterized). The results indicate that quantitative traits are controlled by multiple genes, with overall trait values of offspring biased approaching value of the pure-bred parents. These results therefore suggest a heterosis degradative phenomenon may be operating for some traits in flax.

TABLE 7

Table 7. General statistics of offspring.

FIGURE 5

Figure 5. Frequency distribution of flax fiber related traits. (A: plant height; B: stem length; C: seed yield; D: stem yield; E: fiber yield; F: fiber content).

Based on the assembly information of flax genome, which chromosome-scale pseudomolecules were published (You et al., 2018), it is now possible to realize the positions of QTLs for various traits once all traits are mapped to the scaffold fragments over the entire genome. A total of 12 QTLs were detected within the linkage groups (Table 8), including a QTL detected for plant height on chrLG1 linked with scaffold156, with a LOD score of 3.85 that accounted for 18.77% of the total phenotypic variance for that trait., therefore, the candidate gene of the target trait may actually reside nearby. For stem length, a QTL with LOD score of 3.17 was detected within chrLG8 linked with scaffold2786 that accounted for 11.17% of the total phenotypic variance. Meanwhile, three QTLs were detected for seed yield on chrLG10 and chrLG12 (linked with the scaffold319, scaffold117, and scaffold27), with LOD scores of 2.73, 2.13, and 2.11 that accounted for 10.11–19.33% of the total phenotypic variance. In addition, three QTLs were detected for stem yield on chrLG5 and chrLG15 (linked with the scaffold33 and scaffold355), with LOD scores of 2.97, 2.62, and 2.70 that accounted for 10.91–15.81% of the total phenotypic variance. Two QTLs were detected for fiber yield on chrLG1 and chrLG11 (linked with the scaffold156 and scaffold127), with LOD scores of 3.31 and 3.68, respectively, accounting for 19.09–25.98% of the total phenotypic variance. Finally, two QTLs were detected for fiber content on chrLG5 and chrLG11 (linked with the scaffold680 and scaffold376), with LOD scores of 3.51 and 3.91 that accounted for 15.14 and 13.27% of the total phenotypic variance, respectively.

TABLE 8

Table 8. QTL analysis of fiber related traits in F2 populations.

Interestingly, main effect QTLs for plant height and fiber yield were both detected to the same extent on chrLG1 which spanned the genetic distance from 87.85 to 84.19 cM and contained the same scaffold156. A similar result was observed for stem length, although the different genetic distance spanned, whereby the main effect QTLs were detected to the same scaffold within chrLG5. Therefore, a bold idea emerged that a pleiotropic gene or neighboring genes involved in plant height and/or fiber yield is likely to exist within shared DNA fragments. All of these anchored QTLs should serve as a foundation for later accurate identification of related genes.

Previously, flax QTL mapping based on the genetic linkage map alone could only successfully map a few markers or segments linked to a particular target trait gene, but could not localize target trait genes to physical chromosomal regions. In this study using genomic assembly information in addition to QTL mapping, target trait QTLs could be identified on chromosomes and could be localized to smaller chromosomal segments. However, because each segment may contain several or even dozens of genes, further localization and identification of target genes will still require great effort. Nevertheless, the results presented here provide a foundation for further flax gene mining and characterization.

Discussion

Large-scale genotyping methodologies play an important role in genetic association studies. One such technology, SLAF-seq, has been applied in many plant studies and has produced remarkable results. Due to its relatively higher density, excellent consistency, effectiveness, and lower cost when compared with traditional methods, this method has become a popular genotyping method. Indeed, in our previous study SLAF-seq technology was used to develop numerous novel SSRs in flax (Wu et al., 2017). Subsequently, in this work sequencing results provided a large quantity of SLAF markers to further drive flax genomic research. Moreover, the SLAFs developed here may also be valuable in studies of other flax cultivars and offspring for identification of germplasm or hybrids and for analysis of genetic diversity among cultivars.

SLAF-seq has been widely successful due to several distinguishing characteristics: (i) deep sequencing ensures genotyping accuracy; (ii) reduced representation strategy reduces sequencing costs; (iii) a pre-designed reduced representation scheme optimizes marker efficiency; and (iv) a double barcode system facilitates sequencing of large populations (Sun et al., 2013). Consequently, this technology has been used in numerous genetic linkage map construction studies (Zhang et al., 2013, 2015; Li et al., 2014; Liu et al., 2014, 2016; Wei et al., 2014; Guo et al., 2015; Jiang et al., 2015; Ma et al., 2015; Xu et al., 2015; Zhu et al., 2015, 2016; Yi et al., 2017). In the present study, the total length of the new genetic linkage map, which was constructed using 2,339 SLAF markers, was 1,483.25 cM and spanned 15 LGs with an average distance between adjacent markers of 0.63 cM. To our knowledge, this map is the highest density genetic map of flax currently available and is of high quality.

The general steps of genome assembly mainly include: determination, arrangement, and directional orientation of scaffolds/contigs onto chromosomes (Burton et al., 2013). Although sequencing completed to date toward obtaining the entire flax genome sequence has produced high-quality genomic information, assembly has been difficult due to the high complexity of the genome. However, much work is still needed before the entire flax genome can be assembled, due to genomic complexity and a large number of repeated regions within the genome (Wang Z. et al., 2012). As the carrier of flax genetic information, the genome is the basis for research on genetic mechanisms of flax. Fortunately, chromosome-scale pseudomolecules were refined by optical, physical, and genetic maps in flax (You et al., 2018), this information provides a genome-wide scaffold-based chromosome assembly to guide future physically mapped of target genes. We employed the genome assembly information, with the high-density map of the flax genome constructed here, 12 QTLs were detected and linked with the chromosome scaffolds. Although the order is different in chrLGs and chromosome-scale pseudomolecules, it still in conjunction with sequence information to generate 15 chromosomes. Notably, this method has generated ideas and strategies with wider applicability for use in whole genome assembly of related species.

Interestingly, characteristic fiber-related traits were observed to segregate in flax F₂ offspring. Since “DIANE” is a cultivar of fiber flax and “NY17” is a cultivar of oilseed flax, crossing these cultivars may therefore result in generation of new varieties with high fiber content and would provide information regarding inheritance of fiber-based traits. Moreover, because plant phenotypic traits of height, stem length, stem yield, seed yield, fiber yield, and fiber content exhibit characteristics of quantitative traits, the mapped QTLs of these fiber-related traits to physical locations on chromosome scaffolds. Moreover, the results have helped us make provisional inferences of pleiotropic genes or neighboring genes that influence plant height and fiber yield to explain observed morphological correlations. Selection of new cultivars with high fiber content should thus allow for economically feasible flax production. In addition, our findings suggest that reduced representation genome sequencing (RRGS) might be effectively utilized for QTL mapping of fiber related traits in flax. Using this strategy, the major QTL identified here would serve as a promising starting point for the study of genetic mechanisms influencing fiber yield in flax. Collectively, the results of this study should greatly facilitate subsequent fine mapping of target genes and MAS-based breeding toward the goal of determining the whole flax genome.

Experimental Procedures

Plant Materials

The F₂ mapping population consisted of 112 individuals from the intraspecific crossing of L. usitatissimum L. “DIANE” (male parent, fiber flax) and “NY17” (female parent, oilseed flax). All individual progeny and parent cultivars were grown at the Institute of Industrial Crops of Heilongjiang Academy of Agricultural Sciences, Harbin, China. “DIANE” is a typical fiber flax variety that attains a maximum plant height of 98.2 ~ 110.5 cm and exhibits a thousand seed weight (TSW) of 4.3 ~ 5.0 g. An important trait of “DIANE” is its resistance to saline-alkaline soil and a high fiber content yield. “NY17” is an oil flax variety that has been used for dual production of oil and fiber, with a maximum plant height of 52.9 ~ 62.6 cm and a TSW of 7.3 ~ 8.4 g. An important trait of “NY17” is its high oil content ratio. Field management adhered to essentially standard breeding practices and all flax materials were harvested on time when seeds were mature. All traits were investigated or measured as described at http://www.cgris.net/.

Phenotypic Characterization

In this study, phenotypic descriptions of flax agronomic traits were performed as follows:

Plant height: the height from the junction of cotyledons with the stem to the top of the plant in centimeters (cm); Stem length: the distance between the junction of cotyledons and stem to the base of the first branch (cm); Seed yield: seed weight per plant, with grams (g) as the unit; Stem yield: after removal of soil, leaves, and capsules after plants were threshed, the original stem weight of each plant is called the stem yield (g); Fiber yield: the yield of the stem fiber stripped from the dried stem (g). After measuring plant height, stem length, threshing yield, and seed yield, stem yield was determined after the retting process, fiber yield was determined after removing the hards, and fiber content was calculated as the proportion of the total fiber yield contributed by post-retting stem yield.

DNA Extraction

Young healthy leaves from both parents and 112 F₂ individuals were collected, frozen in liquid nitrogen, then used for DNA extraction. Total genomic DNA was prepared from each plant according to the manufacturer's instructions using a Qiagen DNeasy 96 Plant Kit (Qiagen, CA, USA). DNA concentration and quality were estimated using an ND-1000 spectrophotometer (NanoDrop, DE, USA). Genomic DNA was visualized by electrophoresis on 0.8% agarose gels and was quantified using a NanoDrop 2000 Spectrophotometer (Thermo Scientific, MA, USA). DNA concentration was measured and adjusted to the same level for all samples for use in construction of the SLAF library.

SLAF Library Construction and High-Throughput Sequencing

Here, the L. usitatissimum L. assembled genome size was only 316 Mb with a GC content of 40%. A reference genome of estimated size of 373 Mb was downloaded from GenBank (NCBI GenBank under GenomeProject ID #68161 and Sequence Read Archive accession SRA038451, https://www.ncbi.nlm.nih.gov/sra/?term=SRA038451). Optimal enzymes used to digest DNA were chosen according to the reference genome sequence in combination with the GC content of flax. The SLAF library was constructed based on a pilot SLAF experiment. After samples were purified, genomic DNA was digested using restriction enzymes RsaI and HaeIII. Next, using a dual-indexing strategy (The Arabidopsis Genome Initiative, 2000) dATP was used to add a single nucleotide (A) overhang to the digested fragments at 37°C. Next, T4 DNA ligase was used to ligate duplex tag-labeled sequencing adapters (polyacrylamide gel electrophoresis (PAGE)-purified; Life Technologies, USA) to the A-tailed fragments. DNA fragments (SLAFs) of 314–414 bp (with indexes and adaptors) were excised and diluted for paired-end sequencing on an Illumina HiSeq 2000 sequencing platform (Illumina, Inc; San Diego, CA, USA) conducted by Biomarker Technologies Corporation in Beijing (http://www.biomarker.com.cn/english/) (Figure 6). Real-time monitoring was performed for each cycle during sequencing and ultimately the ratio of high quality raw reads with quality scores greater than Q30 (quality score >30 indicates a 1% chance of error and thus 99% confidence) were used for subsequent analyses. Guanine-cytosine (GC) content was calculated for quality control. Gel-purified products were then diluted and high-throughput paired-end sequencing (101 bp from both ends) was performed using an Illumina HiSeq 2500 System (Illumina, Inc., San Diego, CA, USA) according to the manufacturer's instructions.

FIGURE 6

Figure 6. Distribution map of the insert fragment. The x-axis represents the length of the inserted fragment of reads, and y-axis represents the number of reads.

SLAF-Seq Data Grouping and Genotype Definition

SLAF marker identification and genotyping were performed using procedures described by Sun et al. (2013). Low-quality reads (quality score < 30) were deleted and the remaining raw reads were assigned to 112 individual plant samples according to duplex barcodes of sequences. After trimming the barcodes and the terminal 5-bp positions from each high-quality read, clean reads were clustered together according to their sequence identities. Sequences mapping to the same locus with over 90% identity were assigned to a single SLAF locus (Zhang et al., 2015).

SLAFs with 2, 3, or 4 tags were identified as polymorphic SLAFs and considered to be potential markers. Polymorphic markers were classified into eight segregation patterns (ab × cd, ef × eg, hk × hk, lm × ll, nn × np, aa × bb, ab × cc, and cc × ab). The genetic map was constructed using the F₂ population from a cross between two heterozygous parents after the markers from the segregation pattern of aa × bb were filtered out. Average sequence depths of SLAF markers were >20-fold for parents and >14-fold for progeny. All progeny contained more than 80% of the SLAF markers of the parents, i.e., 80% integrity of SLAF markers were observed in individual plants.

Genetic Map Construction and Quality Evaluation

After removing incomplete and significant segregation distortion markers, 7,417 SLAFs were retained for genetic map construction using a logarithm of odds (LOD) threshold ≥3.0 and a maximum recombination fraction of 0.4, as described by van Ooijen (2011). A HighMap strategy was employed to order the SLAF markers according to the protocol detailed by Liu et al. (2014). Genotyping errors were corrected by the SMOOTH algorithm (van Os et al., 2005). Missing genotypes were imputed using a k-nearest neighbor algorithm (Huang et al., 2011). Map distances in cM were calculated using the Kosambi mapping function (Kosambi, 1944).

Genetic map quality was evaluated using three criteria. The integrity of mapped markers was analyzed in order to correct marker ordering to produce better results. Haplotype maps and heat maps were constructed that directly reflected recombination relationships for determination of the maximum possible double exchange sites among markers within the fifteen linkage groups. Haplotype maps were used to count double exchange and deletion events for all linkage groups such that positions indicated by sites marked by color changes were where reorganization events occurred. Each cell represents a recombination rate between two adjacent markers for each heat map, with recombination rates depicted using different colors ranging from yellow to purple (yellow indicating a lower recombination rate; purple indicating a higher rate). Heat maps also detected potential marker ordering conflicts.

QTL Analysis

QTL analysis was carried out using an R/QTL procedure employing an interval mapping method as previously described (Peichel et al., 2001). Composite interval mapping (CIM) was adopted using a walking speed of 1 cM (Wang W. X. et al., 2012). The significance of each QTL interval was tested by determining a statistical likelihood-ratio (logarithm of the odds or LOD score). The threshold of significance of the LOD score (P = 0.05) was determined using a permutation test (PT) of 1,000 permutations. First, the threshold value was set using the 1,000-permutation test to determine the LOD significance threshold value corresponding to a 0.99 confidence level. If there was no mapping interval, 2-LOD support intervals were constructed using 0.95 confidence intervals (van Ooijen, 1992). If the results of the PT test were not considered, the threshold value was reduced to 3.0 manually. If no interval was observed, then the 3.0 threshold was reduced to 2.0. Calculation of the percentage of phenotypic variance explained (PVE%) by each QTL was obtained from the QTL peak area based on the population variance found within the segregating population.

Author Contributions

JW and QZ conceived and designed the research experiments. JW, GW, HY, and DS implemented the filed experiments. JW, LZ, YY, and TJ performed the phenotyping and genotyping. JW and QZ analyzed the data. JW and TJ designed the overall project. JW wrote the manuscript. JW, ZL, XY, and XW revised the final version of the manuscript. All authors reviewed and approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors gratefully acknowledge the financial support from the National Natural Science Foundation of China (31401451, 31471546), National Bast Fiber Research System of China (CARS-16-S03), Construction of China-CIS-CEE joint laboratory for agricultural science and technology cooperation under the framework of The Belt and Road (2016AE6AE001). The funders had no role in study design, the collection, analysis, interpretation of data, writing of the manuscript, the preparation or decision of the manuscript to publication.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.00885/full#supplementary-material

Table S1. Summary of the sequencing data.

Table S2. SLAF marker numbers and depth of the sequencing data.

Presentation S1. Haplotype map of the genetic maps. Each row represents a marker. Markers are ranked in accordance with the map order. Each of the two columns represents an individual plant allelic set, with green and blue shading of columns representing the parental alleles. The white column represents that the source of alleles cannot be confirmed. The gray columns represent the deleted alleles.

Presentation S2 Heat map of the genetic maps. Markers of each row and column are ranked according to the map order; each small square represents the rate of recombination (r) between the two markers.

Footnotes

1. ^Food and Agriculture Organization. http://www.fao.org/home/en/

References

Burton, J. N., Adey, A., Patwardhan, R. P., Qiu, R., Kitzman, J. O., and Shendure, J. (2013). Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nat. Biotechnol. 31, 1119–1125. doi: 10.1038/nbt.2727

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, G., Wang, X., Long, S., Jaqueth, J., Li, B., Yan, J., et al. (2016). Mapping of QTL conferring resistance to northern corn leaf blight using high-density SNPs in maize. Mol. Breed. 36:4. doi: 10.1007/s11032-015-0421-3

CrossRef Full Text | Google Scholar

Chen, M., Presting, G., Barbazuk, W. B., Goicoechea, J. L., Blackmon, B., Fang, G., et al. (2012). An integrated physical and genetic map of the rice genome. Plant Cell 14, 537–545. doi: 10.1105/tpc.010485

PubMed Abstract | CrossRef Full Text

Cloutier, S., Ragupathy, R., Miranda, E., Radovanovic, N., Reimer, E., Walichnowski, A., et al. (2012). Integrated consensus genetic and physical maps of flax (Linum usitatissimum L.). Theor. Appl. Genet. 125, 1783–1795. doi: 10.1007/s00122-012-1953-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Cloutier, S., Ragupathy, R., Niu, Z., and Duguid, S. (2011). SSR-based linkage map of flax (Linum usitatissimum L.) and mapping of QTLs underlying fatty acid composition traits. Mol. Breed. 28, 437–451. doi: 10.1007/s11032-010-9494-1

CrossRef Full Text | Google Scholar

Foiada, F., Westermeier, P., Kessel, B., Ouzunova, M., Wimmer, V., Mayerhofer, W., et al. (2015). Improving resistance to the European corn borer: a comprehensive study in elite maize using QTL mapping and genome-wide prediction. Theor. Appl. Genet. 128, 875–891. doi: 10.1007/s00122-015-2477-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Gardner, K. A., Wittern, L. M., and Mackay, I. J. (2016). A highly recombined, high-density, eight-founder wheat MAGIC map reveals extensive segregation distortion and genomic locations of introgression segments. Plant Biotechnol. J. 14, 1406–1417. doi: 10.1111/pbi.12504

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, Y., Shi, G., Liu, Z., Zhao, Y., Yang, X., Zhu, J., et al. (2015). Using specific length amplified fragment sequencing to construct the high-density genetic map for Vitis (Vitis vinifera L. × Vitis amurensis Rupr.). Front. Plant Sci. 6:393. doi: 10.3389/fpls.2015.00393

PubMed Abstract | CrossRef Full Text | Google Scholar

Holtz, Y., Ardisson, M., Ranwez, V., Besnard, A., Leroy, P., Poux, G., et al. (2016). Genotyping by sequencing using specific allelic capture to build a high-density genetic map of Durum wheat. PLoS ONE 11:e0154609. doi: 10.1371/journal.pone.0154609

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, X., Zhao, Y., Wei, X., Li, C., Wang, A., Zhao, Q., et al. (2011). Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat. Genet. 44, 32–39. doi: 10.1038/ng.1018

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, B., Liu, W., Xie, D., Peng, Q., He, X., Lin, Y., et al. (2015). High-density genetic map construction and gene mapping of pericarp color in wax gourd using specific-locus amplified fragment (SLAF) sequencing. BMC Genomics 16:1035. doi: 10.1186/s12864-015-2220-y

CrossRef Full Text | Google Scholar

Kosambi, D. D. (1944). The estimation of map distance from recombination values. Ann. Eugen. 12, 172–175.

Google Scholar

Li, B., Tian, L., Zhang, J., Huang, L., Han, F., Yan, S., et al. (2014). Construction of a high-density genetic map based on large-scale markers developed by specific length amplified fragment sequencing (SLAF-seq) and its application to QTL analysis for isoflavone content in Glycine max. BMC Genomics 15:1086. doi: 10.1186/1471-2164-15-1086

CrossRef Full Text | Google Scholar

Li, H., Yang, Y., Zhang, H., Chu, S., Zhang, X., Yin, D., et al. (2016). A genetic relationship between phosphorus efficiency and photosynthetic traits in soybean as revealed by QTL analysis using a high-density genetic map. Front. Plant Sci. 7:924. doi: 10.3389/fpls.2016.00924

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, D., Ma, C., Hong, W., Huang, L., Liu, M., Liu, H., et al. (2014). Construction and analysis of high-density linkage map using high-throughput sequencing data. PLoS ONE 9:e98855. doi: 10.1371/journal.pone.0098855

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, J., Huang, L., Ma, C., Jin, J., Li, C., Wang, R., et al. (2015). Large-scale SNP discovery and genotyping for constructing a high-density genetic map of tea plant using specific-locus amplified fragment sequencing (SLAF-seq). PLoS ONE 10:e0128798. doi: 10.1371/journal.pone.0128798

PubMed Abstract | CrossRef Full Text | Google Scholar

Oh, T. J., Gorman, M., and Cullis, C. A. (2000). RFLP and RAPD mapping in flax (Linum usitatissimum). Theor. Appl. Genet. 101, 590–593. doi: 10.1007/s001220051520

CrossRef Full Text | Google Scholar

Orjuela, J., Garavito, A., Bouniol, M., Arbelaez, J. D., Moreno, L., Kimball, J., et al. (2010). A universal core genetic map for rice. Theor. Appl. Genet. 120, 563–572. doi: 10.1007/s00122-009-1176-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Paterson, S. R., Fowler, T. K., Schmidt, K. L., Yoshinobu, A. S., Yuan, E. S., and Miller, R. B. (1998). Interpreting magmatic fabric patterns in plutons. Lithos 44, 53–82.

Google Scholar

Peichel, C. L., Nereng, K. S., Ohgi, K. A., Cole, B. L., Colosimo, P. F., Buerkle, C. A., et al. (2001). The genetic architecture of divergence between three spine stickleback species. Nature 414, 901–905. doi: 10.1038/414901a

PubMed Abstract | CrossRef Full Text | Google Scholar

Qi, Z., Huang, L., Zhu, R., Xin, D., Liu, C., Han, X., et al. (2014). A high-density genetic map for soybean based on specific length amplified fragment sequencing. PLoS ONE 9:e104871. doi: 10.1371/journal.pone.0104871

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, Q., Jia, G., Hyten, D. L., Jerry, J., Eun-Young, H., Schroeder, S. G., et al. (2015). SNP assay development for linkage map construction, anchoring whole-genome sequence, and other genetic and genomic applications in common bean. G3 Genes Genet. 5, 2285–2290. doi: 10.1534/g3.115.020594

PubMed Abstract | CrossRef Full Text | Google Scholar

Spielmeyer, W., Green, A. G., Bittisnich, D., Mendham, N., and Lagudah, E. S. (1998). Identification of quantitative trait loci contributing to Fusarium wiltresistance on an AFLP linkage map of flax (Linum usitatissimum). Theor. Appl. Genet. 97, 633–641.

Google Scholar

Sun, X., Liu, D., Zhang, X., Li, W., Liu, H., Hong, W., et al. (2013). SLAF-seq: an efficient method of large-scale De novo SNP discovery and genotyping using high-throughput sequencing. PLoS ONE 8:e58700. doi: 10.1371/journal.pone.0058700

PubMed Abstract | CrossRef Full Text | Google Scholar

The Arabidopsis Genome Initiative (2000). >Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature 408(6814): 796–815. doi: 10.1038/35048692

CrossRef Full Text

Thoday, J. M. (1961). Location of polygenes. Nature 191, 368–370. doi: 10.1038/191368a0

CrossRef Full Text | Google Scholar

van Ooijen, J. W. (2011). Multipoint maximum likelihood mapping in a full-sib family of an outbreeding species. Genet. Res. 93, 343–349. doi: 10.1017/S0016672311000279

PubMed Abstract | CrossRef Full Text | Google Scholar

van Ooijen, J. W. (1992). Accuracy of mapping quantitative trait loci in autogamous species. Theor. Appl. Genet. 84, 803–811.

PubMed Abstract | Google Scholar

van Os, H., Stam, P., Visser, R. G., and van Eck, H. J. (2005). SMOOTH: a statistical method for successful removal of genotyping errors from high-density genetic linkage data. Theor. Appl. Genet. 112, 187–194. doi: 10.1007/s00122-005-0124-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Vision, T.J., Brown, D. G., Shmoys, D. B., Durrett, R. T., and Tanksley, S. D. (2000). Selective mapping: a strategy for optimizing the construction of high-density linkage maps. Genetics 155, 407–420.

PubMed Abstract | Google Scholar

Wang, W. X., Huang, S. M., Liu, Y. M., Fang, Z. Y., Yang, L. M., Hua, W., et al. (2012). Construction and analysis of a high-density genetic linkage map in cabbage (Brassica oleracea L. var. capitata). BMC Genomics 13:523. doi: 10.1186/1471-2164-13-523

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z., Hobson, N., Galindo, L., Zhu, S., Shi, D., McDill, J., et al. (2012). The genome of flax (Linum usitatissimum) assembled de novo from short shotgun sequence reads. Plant J. 72, 461–473. doi: 10.1111/j.1365-313X.2012.05093.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, Q., Wang, Y., Qin, X., Zhang, Y., Zhang, Z., Wang, J., et al. (2014). An SNP-based saturated genetic map and QTL analysis of fruit-related traits in cucumber using specific-length amplified fragment (SLAF) sequencing. BMC Genomics 15:1158. doi: 10.1186/1471-2164-15-1158

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, J., Zhao, Q., Wu, G., Zhang, S., and Jiang, T. (2017). Development of novel SSR markers for flax (Linum usitatissimum L.) using reduced-representation genome sequencing. Front. Plant Sci. 7:2018. doi: 10.3389/fpls.2016.02018

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, X., Xu, R., Zhu, B., Yu, T., Qu, W., Lu, L., et al. (2015). A high-density genetic map of cucumber derived from specific length amplified fragment sequencing (SLAF-seq). Front. Plant Sci. 5:768. doi: 10.3389/fpls.2014.00768

PubMed Abstract | CrossRef Full Text | Google Scholar

Yi, L., Gao, F., Siqin, B., Zhou, Y., Li, Q., Zhao, X., et al. (2017). Construction of an SNP-based high-density linkage map for flax (Linum usitatissimum L.) using specific length amplified fragment sequencing (SLAF-seq) technology. PLoS ONE 12:e0189785. doi: 10.1371/journal.pone.0189785

PubMed Abstract | CrossRef Full Text | Google Scholar

You, F. M., Xiao, J., Li, P., Yao, Z., Jia, G., He, L., et al. (2018). Chromosome-scale pseudomolecules refined by optical, physical, and genetic maps in flax. Plant J. Cell Mol. Biol. doi: 10.1111/tpj.13944. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Zhang, Q., Cheng, T., Yang, W., Pan, H., and Zhong, J. (2015). High-density genetic map construction and identification of a locus controlling weeping trait in an ornamental woody plant (Prunus mume Sieb. et Zucc). DNA Res. 22:183–191. doi: 10.1093/dnares/dsv003

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Wang, L., Xin, H., Li, D., Ma, C., Ding, X., et al. (2013). 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. 13:141. doi: 10.1186/1471-2229-13-141

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Y., Huang, L., Chen, L., Yang, J., Wu, J., Qu, M., et al. (2016). A high-density genetic linkage map for cucumber (Cucumis sativus L.): based on specific length amplified fragment (SLAF) sequencing and QTL analysis of fruit traits in cucumber. Front. Plant Sci. 7:437. doi: 10.3389/fpls.2016.00437

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Y., Yin, Y., Yang, K., Li, J., Sang, Y., Huang, L., et al. (2015). Construction of a high-density genetic map using specific length amplified fragment markers and identification of a quantitative trait locus for anthracnose resistance in walnut (Juglans regia L.). BMC Genomics 16:614. doi: 10.1186/s12864-015-1822-8

PubMed Abstract | CrossRef Full Text | Google Scholar