Impact Factor 3.678

The world's most-cited Plant Sciences journal

Original Research ARTICLE

Front. Plant Sci., 21 March 2016 | https://doi.org/10.3389/fpls.2016.00334

Genome-Wide Single-Nucleotide Polymorphisms Discovery and High-Density Genetic Map Construction in Cauliflower Using Specific-Locus Amplified Fragment Sequencing

Zhenqing Zhao1, Honghui Gu1*, Xiaoguang Sheng1, Huifang Yu1, Jiansheng Wang1, Long Huang2 and Dan Wang2
  • 1Institute of Vegetables, Zhejiang Academy of Agricultural Sciences, Hangzhou, China
  • 2Biomarker Technologies Corporation, Beijing, China

Molecular markers and genetic maps play an important role in plant genomics and breeding studies. Cauliflower is an important and distinctive vegetable; however, very few molecular resources have been reported for this species. In this study, a novel, specific-locus amplified fragment (SLAF) sequencing strategy was employed for large-scale single nucleotide polymorphism (SNP) discovery and high-density genetic map construction in a double-haploid, segregating population of cauliflower. A total of 12.47 Gb raw data containing 77.92 M pair-end reads were obtained after processing and 6815 polymorphic SLAFs between the two parents were detected. The average sequencing depths reached 52.66-fold for the female parent and 49.35-fold for the male parent. Subsequently, these polymorphic SLAFs were used to genotype the population and further filtered based on several criteria to construct a genetic linkage map of cauliflower. Finally, 1776 high-quality SLAF markers, including 2741 SNPs, constituted the linkage map with average data integrity of 95.68%. The final map spanned a total genetic length of 890.01 cM with an average marker interval of 0.50 cM, and covered 364.9 Mb of the reference genome. The markers and genetic map developed in this study could provide an important foundation not only for comparative genomics studies within Brassica oleracea species but also for quantitative trait loci identification and molecular breeding of cauliflower.

Introduction

Cauliflower (Brassica oleracea var. botrytis, 2n = 2x = 18) is an important vegetable crop worldwide. It is considered a vital source of vitamins, dietary fiber, antioxidants, and anti-carcinogenic compounds (Volden et al., 2009; Picchi et al., 2012). In 2013, global cauliflower cultivation had spread across ~1.2 million hectares, with total production of ~20.9 million metric tons (http://faostat.fao.org/). Due to the economic importance and nutritional value of cauliflower, great efforts have been taken to improve its yield and quality. Nevertheless, breeding methods used are essentially conventional and less-effective, resulting in relative slow progress of the cauliflower breeding program (Gu et al., 2014). Modern strategies such as marker-assisted selection (MAS) are therefore necessary to accelerate cauliflower genetic improvement.

Molecular markers and genetic maps are considered an important foundation for quantitative trait loci (QTL) mapping and MAS (Kato et al., 2014; Li et al., 2014; Wu et al., 2014; Zhang et al., 2014). As one of the subspecies of B. oleracea, cauliflower has been employed to develop a series of segregating populations to construct several C-genome genetic maps, by crossing with other subspecies of B. oleracea, including kale (B. oleracea var. acephala) (Kianian and Quiros, 1992), broccoli (B. oleracea var. italica) (Li et al., 2003; Gao et al., 2007), collard (B. oleracea var. acephala) (Hu et al., 1998), and brussel sprouts (B. oleracea var. gemmifera) (Sebastian et al., 2000, 2002), and to identify several QTLs involved in common traits such as flowering time (Kianian and Quiros, 1992), glucosinolate profile (Gao et al., 2007), and leaf traits (Sebastian et al., 2002). However, because of the organ specificity of cauliflower curd, inter-subspecies populations are unsuitable for QTL analysis of many curd-specific traits, which are most important for cauliflower breeding. Hence, a cauliflower × cauliflower based population has extensive potential for marker discovery, genetic mapping, and QTL analysis.

Several molecular markers such as restriction fragment length polymorphism (RFLP), randomly-amplified polymorphic DNA (RAPD), amplified fragment length polymorphisms (AFLP), sequence-related amplified polymorphism (SRAP), and simple sequence repeats (SSR) have been widely used in plant genetic research. However, the utilization of these markers has always been limited by the higher time and cost requirements as well as limited marker resource (Zhao et al., 2014a). The advent of massive, parallel, next-generation sequencing (NGS) technologies have accelerated and simplified the identification of sequence variants, enabling large-scale single-nucleotide polymorphism (SNP) discovery throughout the genome (Zhou et al., 2014). Considering that whole-genome deep re-sequencing is still expensive and usually not necessary (Davey et al., 2011; Wei et al., 2014), several simplified and cost-effective methods for SNP discovery and high-throughput genotyping have been developed, such as reduced representation library (RRL) sequencing (Van Tassell et al., 2008), restriction-site associated DNA sequencing (RAD-seq; Miller et al., 2007), two-enzyme genotyping by sequencing (GBS) (Poland et al., 2012), and sequence-based genotyping (SBG) (Truong et al., 2012). More recently, specific-locus amplified fragments sequencing (SLAF-seq) was developed as a streamlined RRL sequencing approach for high-resolution de novo SNP discovery and genotyping (Sun X. et al., 2013). A high density genetic map including 1233 high-quality markers has been developed using this strategy on a sesame F2 population (Zhang et al., 2013). This study showed that SLAF sequencing was a powerful high-throughput technique for plant genome research. To date, this strategy has also been successfully applied in several other species including soybean (Li et al., 2014; Qi et al., 2014), cucumber (Wei et al., 2014; Xu et al., 2014), tea plant (Ma et al., 2015), and grape (Guo et al., 2015).

In this study, we generated a double-haploid (DH) population derived from a cross between two different types of cauliflower common in production, including an advanced inbred line of traditional compact-curd cauliflower and a DH line of loose-curd cauliflower. Based on this cauliflower × cauliflower population, SLAF-seq was then employed to detect large-scale SNPs and construct a high-density genetic map that could be used to provide a platform for future QTL mapping and MAS.

Materials and Methods

Mapping Population Development and Genomic DNA Isolation

An advanced inbred line of compact-curd cauliflower “4305” (F8) and a DH line of loose-curd cauliflower “ZN198” were used to develop the DH mapping population (Figure S1). There are significant morphological differences between the two homozygous lines, especially for several agronomically important traits including curd size, curd weight, curd shape, and curding-time. A microspore culture protocol as previously described by Gu et al. (2014) was used to produce regenerated plants from a single F1 plant of the cross “4305” × “ZN198.” The ploidy level of all the regenerated plants achieved was estimated using an FCM Ploidy Analyzer (Partec GmbH, Germany) and only diploids were selected to construct the mapping population. Parents and DH lines were planted in the experiment field of Zhejiang Academy of Agricultural Sciences in Hangzhou, China and were preserved for long-term utilization by artificial selfing.

Young leaves were collected and genomic DNA was isolated according to a modified version of the cetyltrimethyl ammonium bromide (CTAB) procedure (Doyle and Doyle, 1986). DNA concentration and quality were detected using an ND-1000 spectrophotometer (NanoDrop, Wilmington, DE, USA) and electrophoresis on 1.0% agarose gel with a standard lambda DNA.

SLAF Library Preparation and High Throughput Sequencing

A SLAF-seq strategy was used, as previously described by Sun X. et al. (2013), with modifications. First, reference genome of B. oleracea (cabbage, B. oleracea var. capitata, http://www.ocri-genomics.org/bolbase/, Liu S. et al., 2014) was used to design marker discovery experiments by simulating in silico, the number of markers produced by different enzymes. A SLAF pilot experiment was performed to determine the optimized enzymes and restriction fragment size, while the SLAF library was conducted based on the pre-designed scheme. Subsequently, the genomic DNA (2 μg) was digested with 3.6 units RsaI (New England Biolabs, NEB, USA) in 20 μl volume containing 1 × NEB buffer at 37°C for 2 h, and then a single nucleotide (A) was added by using 6 units Klenow Fragment (3′ → 5′ exo) (NEB) and 10 nmol dATP at 37°C for 1 h. Duplex tag-labeled sequencing adapters (PAGE-purified, Life Technologies, USA) were then ligated to the A-tailed fragments using T4 DNA ligase (NEB) by incubating overnight at 16°C, then 65°C for 20 min to heat deactivate the T4 ligase. Polymerase chain reaction (PCR) was performed in a 100 μl final volume reaction mixture, which contained diluted restriction-ligation DNA samples, PCR primers (forward sequence: 5′-AATGATACGGCGACCACCGA-3′, reverse sequence: 5′-CAAGCAGAAGACGGCATACG-3′), dNTP, MgCl2 and Q5® High-Fidelity DNA Polymerase (PAGE-purified, Life Technologies). The PCR cycles were 98°C for 3 min, 18 cycles of 98°C for 10 s, 65°C for 30 s, 72°C for 30 s, followed by an extension step of 5 min at 72°C before storage at 4°C. The amplification products were purified using Agencourt AMPure XP beads (Beckman Coulter, High Wycombe, UK) and pooled, followed by separation on 2% agarose gel electrophoresis. Fragments ranging from 244 to 314 bp (with indices and adaptors) in size were gel-purified using a QIAquick gel extraction kit (Qiagen, Hilden, Germany) and diluted for pair-end sequencing (125 bp at each end) using an Illumina HiSeq 2500 system (Illumina, Inc.; San Diego, CA, USA) at Beijing Biomarker Technologies Corporation, according to the manufacturer recommendations.

Sequence Data Grouping and Genotyping

SLAF marker identification and genotyping were performed as previously described by Sun X. et al. (2013) and Zhang et al. (2015). First, raw reads were demultiplexed to individuals according to the barcode sequences. The reads with quality scores < Q30 (a quality score of 30; indicating 0.1% chance of an error, and thus 99.9% confidence) were filtered out. After the barcodes and the terminal 5-bp positions were trimmed, high-quality reads from the same sample were mapped onto the reference genome sequence using SOAP software (Li et al., 2008). Subsequently, sequences locating at the same position with over 95% identity were grouped into one SLAF locus. SNP of each locus were firstly detected between parents, and SLAFs with less than three SNPs were used to define alleles. As cauliflower is diploid, only SLAFs with two to four alleles were identified as polymorphic and considered potential markers. Each polymorphic SLAF marker was then classified into eight segregation patterns in population (ab × cd, ef × eg, hk × hk, lm × ll, nn × np, aa × bb, ab × cc, and cc × ab) as previously described by Zhang et al. (2013). Since the DH mapping population used here was derived from two homozygous lines, only the SLAF markers showing aa × bb segregation pattern were used for map construction. Genotype scoring was then performed using a Bayesian approach to further ensure the genotyping quality, and high-quality SLAF markers for the genetic mapping were filtered by criteria as previously described by Sun X. et al. (2013).

Genetic Map Construction

The genetic map was constructed as previously described by Zhang et al. (2015). Marker loci were allocated primarily into nine linkage groups (LGs) based on their locations on the reference genome. Markers with the modified logarithm of odds (MLOD) score < 5 were filtered to further confirm the marker robustness. A newly developed High Map strategy (Liu D. et al., 2014) was applied to order the SLAF markers and correct genotyping errors within LGs. MSTmap algorithm was used to order SLAFs markers (Wu et al., 2008) and the SMOOTH algorithm (Hans et al., 2005) was used to correct genotyping errors following marker ordering. Map distances were estimated using the Kosambi mapping function (Kosambi, 1943). Since all the mapped markers have been mapped onto the reference genome of B. oleracea by SOAP software (Li et al., 2008), the collinearity of physical map and genetic map were visualized by RSCRIPT language following the tutorial introduction.

Results

The Mapping Population

In total, 136 regenerated plants with different genotypes were obtained from “4305” × “ZN198” F1 by microspore culture. Several ploidy levels were detected in these plants, including haploids, diploids, polyploids, and chimeras (Figure S2). Finally, 79 diploid plants showing a spontaneously doubling ratio of 58.1% were obtained and used as the mapping population.

SLAF Sequencing and Genotyping

A total of 12.47 Gb raw data containing 77.92 M pair-end reads were generated from the high-throughput sequencing (Data has been submitted to National Center of Biotechnology Information, the BioProject ID was PRJNA307521), with a GC (guanine-cytosine) content of 38.12%. The rate of high-quality reads with quality scores >30 reached 91.07%. After read clustering and filtering, 81,311 high-quality SLAFs with even distribution throughout the genome were identified (Figure 1A; Table 1; Table S1). The average sequencing depths of these SLAFs were 52.66-fold for female inbred line (4305), 49.35-fold for male DH line (ZN198), and 4.71-fold for each progeny of the DH population (Figure 2). Based on the allele number and sequence difference, the 81,311 SLAFs were grouped into three types including polymorphic, non-polymorphic, and repetitive. 6815 polymorphic SLAFs out of all SLAFs were obtained, showing a polymorphism rate of 8.38% (Figure 1B; Table S2). Interestingly, the polymorphism rate of SLAFs developed on each chromosome was quite different. In the current study, SLAFs on chromosome C01, C07, and C08 showed significantly higher polymorphism rate than the others (Table 2).

FIGURE 1
www.frontiersin.org

Figure 1. Distribution of the total SLAFs (A) polymorphic SLAFs (B) and mapped SLAF markers (C) on the reference genome. The x-axis represents chromosome length and the y-axis indicates chromosome code. Each yellow bar stands for a chromosome, and deeper color from yellow to black means more SLAFs on the corresponding location.

TABLE 1
www.frontiersin.org

Table 1. SLAF-seq data summary for the mapping population.

FIGURE 2
www.frontiersin.org

Figure 2. Sequencing depths and number of SLAFs of all the individuals of mapping population. The x-axis indicates the average depths (A) and number of SLAFs (B) the y-axis represents individual accessions.

TABLE 2
www.frontiersin.org

Table 2. Polymorphism rate of SLAFs developed on each chromosome.

Subsequently, all these polymorphic SLAFs were genotyped separately for parents and population individuals. A total of 6568 polymorphic SLAFs from total 6815 were successfully encoded, of which 4736 SLAFs were classified as the expected aa × bb segregation pattern, following the genotype encoding rule (Figure 3). These 4736 SLAFs were further filtered based on the criteria considering segregation distortion, sequencing depth, and data integrity. Finally, 1776 high-quality markers, with the average sequencing depths of 63.05-fold in the female parent, 59.75-fold in the male parent, and 4.50-fold in each DH individual, were employed to construct the genetic map (Figure 1C; Table 1).

FIGURE 3
www.frontiersin.org

Figure 3. Number of SLAF markers in each segregation patterns.

Main Characteristics of the Genetic Map

All 1776 high-quality markers were successfully assigned onto 9 LGs according to their locations on the reference genome and the MLOD scores with other markers (at least one MLOD score >5). The average data integrity of these 1776 SLAF markers reached 95.68%. The final genetic map was constructed, following linkage analysis for each of the 9 LGs, which was designated according to the corresponding chromosome number of the reference genome (Figure 4; Table 3). A total of 1776 mapped SLAF markers containing 2741 SNPs spanned a total genetic length of 890.01 cM, with an average marker interval of 0.50 cM and covered 364.9 Mb of the reference genome. The genetic length, marker number, and average marker interval of single LG ranged from 41.90 (C01) to 200.92 (C06), 144 (C01) to 277 (C03), and 0.23 (C02) to 0.87 cM (C06), respectively. The max gap was 12.94 cM, located on C07. Detailed data of the genetic map and markers are presented in Table S3.

TABLE 3
www.frontiersin.org

Table 3. Basic characteristics of the nine LGs.

FIGURE 4
www.frontiersin.org

Figure 4. Distribution of SLAF markers on eight linkage groups of cauliflower. A black bar indicates a SLAF marker. The x-axis represents linkage group number and the y-axis indicates genetic distance.

Visualization and Evaluation of the Genetic Map

The mapped markers were anchored on the reference genome and the correlation of genetic and physical position is shown in Figure S4. Generally, a sufficient genome coverage and accurate genetic location of markers were revealed by the consecutive curves. The genetic arrangements of most markers were also considered to coincide with their physical direction based on the falling trend of the curve. However, a significant inversion can be observed intuitively on pseudo-chromosome 4. The detailed collinearity and marker location are described in Figure S3; Table S3.

Haplotype maps were generated for each DH individual with two parents as control using all the mapped markers in order to detect double recombination and deletion, thus to reflect potential genotyping and marker-order errors (West et al., 2006). In this study, there was no double recombination and deletion found in any linkage group (Figure S5). In addition, the map quality was also evaluated by heat maps, which could intuitively display recombination relationships among markers within each single LG. Pair-wise recombination rates could be visualized by different color levels from yellow to purple. Here, yellow color generally showed diagonal distribution in the heat map for each of the nine LGs (Figure S6), indicating that the mapped marker had been correctly ordered.

Discussion

Due to the genome-wide abundance and high-throughput nature, SNP markers are playing an increasingly important role in plant research, such as genetic map construction, novel gene discovery, evolutionary analysis, and MAS within breeding programs (Liu et al., 2012). SLAF-seq is a newly reported enhanced RRL sequencing solution for large-scale SNP discovery and genotyping (Sun X. et al., 2013). It can produce large amounts of sequence-based information and handle any density distribution throughout the whole genome (Chen et al., 2013). Moreover, compared with another widely used, NGS-based method, RAD-seq, SLAF-seq shows higher reproducibility (Qi et al., 2014; Xu et al., 2014), and locus discrimination efficiency (Wei et al., 2014) due to its paired-ends sequencing strategy and longer read length (30–50 bp), respectively. In this study, we used this cost-effective strategy to develop and sequence a total of 81,311 SLAFs. Among these, 6815 SLAFs showing polymorphism between two parents (“4305” and “ZN198”) were identified with a polymorphic rate of 8.38%. Although the rate was much lower than that previously obtained by whole-genome re-sequencing in cabbage (Wang et al., 2012), the markers developed herein covered all pseudo-chromosomes and were evenly distributed throughout the reference genome (Figure 1B). In addition, the reads quality score of > 30 and an average sequencing depth of 57-fold for parents (Table 1) ensured high genotyping accuracy. Therefore, this set of SLAFs has great potential for use in both genomic study and breeding application of cauliflower. Additionally, our study further demonstrated that SLAF-seq is an efficient strategy for genome-wide SNP identification by sampling and sequencing a reduced set of representative genome regions instead of the whole genome.

Long-term domestication and selective breeding have resulted in abundant variations within B. oleracea species, hence several morphologically different subspecies currently available worldwide, such as cabbage, brussel sprouts, cauliflower, and broccoli (Kennard et al., 1994). In order to make the most of genetic variation within B. oleracea, traditional genomic researches usually employ different subspecies as parental materials to develop markers, construct genetic maps, and identify QTLs (Wang et al., 2012). However, these advances have been not comprehensive enough to support the genetic improvement of some specific organ, like head of cabbage or curd of cauliflower. For this reason, increasing attention has been given to genomic study based on intra-subspecies populations. For instance, Wang et al. (2012) developed more than 5000 SSR and SNP markers using two cabbage inbred lines and constructed a high-density genetic map including 1227 markers. Based on another cabbage × cabbage population, Lv et al. (2014) identified 707 InDels (insertion–deletions) and detected 13 reliable QTLs associated with five important heading traits. Among the B. oleracea varieties, cauliflower is the only plant where the immature inflorescence shows hypertrophic structure (curd). Therefore, the improvement of curd-specific traits like weight, size, shape, color, and content of nutritional compounds are important goals for the cauliflower breeding program. In this study, we used two different types of cauliflower, commonly produced in China, including an advanced inbred line of traditional compact-curd cauliflower and a DH line of loose-curd cauliflower, to identify SNPs and construct a genetic map. Firstly, the genetic base difference between loose-curd cauliflower and compact-curd cauliflower (Zhao et al., 2014b) would ensure acceptable marker polymorphism. Secondly, loose-curd cauliflower is now rapidly widespread and becoming a main cultivated type of cauliflower in China (Zhao et al., 2012). It shows varied agronomic characteristics (Zhao et al., 2013) and nutritional compounds (Gu et al., 2015) with a compact curd. The cauliflower-based SNP information should be helpful to enhance the future breeding program. In our lab, QTL mapping that aims to uncover the genetic factors controlling a series of curd-specific traits and to identify related markers based on the current study is now underway. In addition, loose-curd cauliflower originating from a portion of the genetic variation of domesticated cauliflower has undergone founder effects and intense selection for locally favorable curd characteristics (Zhao et al., 2014b). The significantly higher polymorphism percentage of SLAFs developed on C01, C07, and C08 between two parents suggest that these chromosomes may carry more genes/QTLs related to curd traits that have been strongly selected by breeders and cultivators. An interesting future task is to more thoroughly dissect the domestication of loose-curd cauliflower by using more abundant germplasm of both type of this crop.

Due to the limited genetic base within cauliflower (Zhao et al., 2014b), there has been very few genetic mapping studies using cauliflower based on intra-subspecies cross. The only case to our knowledge is a cauliflower × cauliflower-based map spanning a genetic length of 668.4 cM with 234 AFLP and 21 NBS (nucleotide binding site) markers (Gu et al., 2007). However, the limited marker numbers and density confined its further application. Besides, the markers in this map were difficult to be anchored onto the reference genome or transferred to other genetic maps. For all B. oleracea sub-species, the genetic map with highest marker density developed to date was based on B. oleracea Genome Sequencing Project (BrGSP) and whole-genome re-sequencing in cabbage. In this case, the map comprised 1227 markers, with an average interval of 0.98 cM, and was applied to anchor assembled scaffolds onto pseudo-chromosomes in BrGSP (Wang et al., 2012). In contrast, the SLAF-seq-based map in the present study contained 1776 SLAF markers including 2741 SNPs and covered 364.9 Mb of the cabbage genome. The marker quantity and resolution herein were significantly improved. However, many markers on the current map were noted to be highly clustered, some of them even located on the same locus, although their corresponding physical positions were quite different (Table S3). Therefore, despite the average marker interval being as short as 0.50 cM, there were still several obvious gaps on some regions. The similar phenomenon has been reported previously as well (Wang et al., 2012; Qi et al., 2014; Cai et al., 2015), indicating that the genetic differences or recombination between mapping parents on some genome regions were inadequate (Sun L. et al., 2013). In any case, the cluster markers may be separated distinctly if they are used for a different or broader population. For single LG, we also noted the discrepancy between the genetic length and corresponding pseudo-chromosome length of B. oleracea draft genome. For example, C06 is the longest LG in our map with a genetic length of 200.92 cM, but it is a relative short pseudo-chromosome in reference genome draft (http://www.ocri-genomics.org/bolbase/). This case, though common, could represent the extensive regions of increased/decreased recombination or could potentially indicate areas of the draft sequence that might require further refinement.

Comparative analyses of Brassica accessions have demonstrated that both genes relocation and sequence polymorphisms between species are common in the Brassica genome (Hu et al., 1998; Sebastian et al., 2000). Here, we show that the genome of cabbage and cauliflower are highly collinear, with macro collinearity punctuated by rearrangements mainly involving translocations and inversions. The majority of the observed rearrangements involved very short distances, but an apparent inversion spanning distance of more than 20 cM was noted at the central section of linkage group 04 between Marker13078 and Marker17566 (Figure S3; Table S3). Similar violation of genetic collinearity have also been identified within both intraspecific and interspecific level in other Brassica species, which may be responsible for the high degree of morphological polymorphism among these species (Sebastian et al., 2000; Bancroft et al., 2011). Since the cabbage and cauliflower subspecies shared a polyphyletic origin in primitive B. oleracea populations (Truco et al., 1996), our results imply that the inverted region seem to have been involved in the divergence of cauliflower and cabbage. Although additional genomic data from other varieties/lines are necessary to trace the genome evolution, comparative information about the collinearity between these two closely related subspecies has important implications for further marker and gene identification in cauliflower through the use of sequence data from cabbage as a genomic model.

In conclusion, we demonstrated the application of SLAF-seq strategy in cauliflower for large-scale SNP discovery and high-density genetic map construction. The parents used in this study are not only representatives of cultivated cauliflower, but also elite lines in our breeding program. Therefore, the genome-wide SNPs identified and the linkage map developed in this study will provide an important foundation not only for QTL identification and map-based cloning, especially of curd development-related traits, but also for MAS within cauliflower breeding programs.

Author Contributions

ZZ and HG conceived and performed experiments and wrote the manuscript; XS, HY, and JW were involved in data analysis; LH and DW performed the SLAF library construction and sequencing. All the authors have commented, read and approved the final manuscript.

Funding

This project is supported by the National Natural Science Foundation of China (grant number, 31501768); the Research project in Zhejiang Province Science and Technology department (grant number, 2012C12903-3-4, 2016C32102); the Research project in Zhejiang Academy of Agricultural Sciences (grant number, 2014CX031, 2015R23R08E05).

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.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.00334

References

Bancroft, I., Morgan, C., Fraser, F., Higgins, J., Wells, R., Clissold, L., et al. (2011). Dissecting the genome of the polyploid crop oilseed rape by transcriptome sequencing. Nat. Biotechnol. 29, 762–766. doi: 10.1038/nbt.1926

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, C., Cheng, F., Wu, J., Zhong, Y., and Liu, G. (2015). The first high-density genetic map construction in tree peony (Paeonia Sect. Moutan) using genotyping by specific-locus amplified fragment sequencing. PLoS ONE 10:e0128584. doi: 10.1371/journal.pone.0128584

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, S., Huang, Z., Dai, Y., Qin, S., Gao, Y., Zhang, L., et al. (2013). The development of 7E Chromosome-specific molecular markers for Thinopyrum elongatum based on SLAF-seq technology. PLoS ONE 8:e65122. doi: 10.1371/journal.pone.0065122

PubMed Abstract | CrossRef Full Text | Google Scholar

Davey, J. W., Hohenlohe, P. A., Etter, P. D., Boone, J. Q., Catchen, J. M., and Blaxter, M. L. (2011). Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nat. Rev. Genet. 12, 499–510. doi: 10.1038/nrg3012

PubMed Abstract | CrossRef Full Text | Google Scholar

Doyle, J., and Doyle, J. (1986). A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem. Bull. 19, 11–15.

Google Scholar

Gao, M., Li, G., Yang, B., Qiu, D., Farnham, M., and Quiros, C. (2007). High-density Brassica oleracea linkage map: identification of useful new linkages. Theor. Appl. Genet. 115, 277–287. doi: 10.1007/s00122-007-0568-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu, H., Wang, J., Zhao, Z., Sheng, X., Yu, H., and Huang, W. (2015). Characterization of the appearance, health-promoting compounds, and antioxidant capacity of the florets of the loose-curd cauliflower. Int. J. Food Properties 18, 392–402. doi: 10.1080/10942912.2013.831445

CrossRef Full Text | Google Scholar

Gu, H., Zhao, Z., Sheng, X., Yu, H., and Wang, J. (2014). Efficient doubled haploid production in microspore culture of loose-curd cauliflower (Brassica oleracea var. botrytis). Euphytica 195, 467–475. doi: 10.1007/s10681-013-1008-x

CrossRef Full Text | Google Scholar

Gu, Y., Zhao, Q. C., and Sun, D. L. (2007). Construction of genetic linkage map and localization of NBS-LRR like resistance gene analogues in cauliflower (Brassica oleracea var. botrytis). Hereditas 29, 751–757. doi: 10.1360/yc-007-0751

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

Hans, V. O., Stam, P., Visser, R. G. F., and Eck, H. J. V. (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

Hu, J., Sadowski, J., Osborn, T. C., Landry, B. S., and Quiros, C. F. (1998). Linkage group alignment from four independent (Brassica oleracea) RFLP maps. Genome 41, 226–235. doi: 10.1139/g98-007

CrossRef Full Text | Google Scholar

Kato, S., Sayama, T., Fujii, K., Yumoto, S., Kono, Y., Hwang, T. Y., et al. (2014). A major and stable QTL associated with seed weight in soybean across multiple environments and genetic backgrounds. Theor. Appl. Genet. 127, 1365–1374. doi: 10.1007/s00122-014-2304-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kennard, W. C., Slocum, M. K., Figdore, S. S., and Osborn, T. C. (1994). Genetic analysis of morphological variation in Brassica Oleracea using molecular markers. Theor. Appl. Genet. 87, 721–732. doi: 10.1007/BF00222898

PubMed Abstract | CrossRef Full Text | Google Scholar

Kianian, S. F., and Quiros, C. F. (1992). Generation of a Brassica oleracea composite RFLP map: linkage arrangements among various populations and evolutionary implications. Theor. Appl. Genet. 84, 544–554. doi: 10.1007/bf00224150

PubMed Abstract | CrossRef Full Text | Google Scholar

Kosambi, D. D. (1943). The estimation of map distances from recombination values. Ann. Eugen. 12, 172–175. doi: 10.1111/j.1469-1809.1943.tb02321.x

CrossRef Full Text | Google Scholar

Li, B., Ling, T., Zhang, J., Long, H., 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, G., Gao, M., Yang, B., and Quiros, C. F. (2003). Gene to gene alignment between the Brassica and Arabidopsis genomes by transcriptional mapping. Theor. Appl. Genet. 107, 168–180. doi: 10.1007/s00122-003-1236-x

CrossRef Full Text | Google Scholar

Li, R., Li, Y., Kristiansen, K., and Wang, J. (2008). SOAP: short oligonucleotide alignment program. Bioinformatics 24, 713–714. doi: 10.1093/bioinformatics/btn025

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, J., Huang, S., Sun, M., Liu, S., Liu, Y., Wang, W., et al. (2012). An improved allele-specific PCR primer design method for SNP marker analysis and its application. Plant Methods 8:34. doi: 10.1186/1746-4811-8-34

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, S., Liu, Y., Yang, X., Tong, C., Edwards, D., Parkin, I. A., et al. (2014). The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes. Nat. Commun. 5, 3930. doi: 10.1038/ncomms4930

PubMed Abstract | CrossRef Full Text | Google Scholar

Lv, H., Wang, Q., Zhang, Y., Yang, L., Fang, Z., Wang, X., et al. (2014). Linkage map construction using InDel and SSR markers and QTL analysis of heading traits in Brassica oleracea var. capitata L. Mol. Breed. 34, 1–12. doi: 10.1007/s11032-014-0019-1

CrossRef Full Text | Google Scholar

Ma, J. Q., Huang, L., Ma, C. L., Jin, J. Q., Li, C. F., Wang, R. K., 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

Miller, M. R., Dunham, J. P., Amores, A., Cresko, W. A., and Johnson, E. A. (2007). Rapid and cost effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers. Genome Res. 17, 240–248. doi: 10.1101/gr.5681207

PubMed Abstract | CrossRef Full Text | Google Scholar

Picchi, V., Migliori, C., Scalzo, R. L., Campanelli, G., Ferrari, V., and Cesare, L. F. D. (2012). Phytochemical content in organic and conventionally grown Italian cauliflower. Food Chem. 130, 501–509. doi: 10.1016/j.foodchem.2011.07.036

CrossRef Full Text | Google Scholar

Poland, J. A., Brown, P. J., Sorrells, M. E., and Jannink, J. L. (2012). Development of high-density genetic maps for barley and wheat using a novel two-enzyme genotyping-by-sequencing approach. PLoS ONE 7:e32253. doi: 10.1371/journal.pone.0032253

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

Sebastian, R. L., Howell, E. C., King, G. J., Marshall, D. F., and Kearsey, M. J. (2000). An integrated AFLP and RFLP Brassica oleracea linkage map from two morphologically distinct doubled-haploid mapping populations. Theor. Appl. Genet. 100, 75–81. doi: 10.1007/s001220050011

CrossRef Full Text | Google Scholar

Sebastian, R. L., Kearsey, M. J., and King, G. J. (2002). Identification of quantitaive trait loci controlling developmental characteristics of Brassica oleracea L. Theor. Appl. Genet. 104, 601–609. doi: 10.1007/s001220100743

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, L., Yang, W., Zhang, Q., Cheng, T., Pan, H., Xu, Z., et al. (2013). Genome-wide characterization and linkage mapping of simple sequence repeats in mei (Prunus mume Sieb. et Zucc.). PLoS ONE 8:e59562. doi: 10.1371/journal.pone.0059562

PubMed Abstract | CrossRef Full Text | 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

Truco, M. J., Hu, J., Sadowski, J., and Quiros, C. F. (1996). Inter-and intra-genomic homology of the Brassica genomes: implications for their origin and evolution. Theor. Appl. Genet. 93, 1225–3123. doi: 10.1007/BF00223454

PubMed Abstract | CrossRef Full Text | Google Scholar

Truong, H. T., Ramos, A. M., Yalcin, F., De, R. M., Hj, V. D. P., Huvenaars, K. H., et al. (2012). Sequence-based genotyping for marker discovery and co-dominant scoring in germplasm and populations. PLoS ONE 7:e37565. doi: 10.1371/journal.pone.0037565

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Tassell, C., Smith, T., Matukumalli, L., Taylor, J., Schnabel, R., Lawley, C., et al. (2008). SNP discovery and allele frequency estimation by deep sequencing of reduced representation libraries. Nat. Methods 5, 247–252. doi: 10.1038/nmeth.1185

PubMed Abstract | CrossRef Full Text | Google Scholar

Volden, J., Bengtsson, G. B., and Wicklund, T. (2009). Glucosinolates; L-ascorbic acid, total phenols, anthocyanins, antioxidant capacities and colour in cauliflower (Brassica oleracea L. ssp. botrytis); effects of long-term freezer storage. Food Chem. 112, 967–976. doi: 10.1016/j.foodchem.2008.07.018

CrossRef Full Text | Google Scholar

Wang, W., Huang, S., Liu, Y., Fang, Z., Yang, L., 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

Wei, Q., Wang, Y., Qin, X., Zhang, Y., Zhang, Z., Jing, W., 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

West, M. A., Van, L. H., Kozik, A., Kliebenstein, D. J., Doerge, R. W., St Clair, D. A., et al. (2006). High-density haplotyping with microarray-based expression and single feature polymorphism markers in Arabidopsis. Genome Res. 16, 787–795. doi: 10.1101/gr.5011206

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, J., Li, L. T., Li, M., Khan, M. A., Li, X. G., Chen, H., et al. (2014). High-density genetic linkage map construction and identification of fruit-related QTLs in pear using SNP and SSR markers. J. Exp. Bot. 65, 5771–5781. doi: 10.1093/jxb/eru311

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Y., Bhat, P. R., Close, T. J., and Lonardi, S. (2008). Efficient and accurate construction of genetic linkage maps from the minimum spanning tree of a graph. PLoS Genet. 4:e1000212. doi: 10.1371/journal.pgen.1000212

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, X., Xu, R., Zhu, B., Yu, T., Qu, W., Lu, L., et al. (2014). 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

Zhang, D., Hua, Y., Wang, X., Zhao, H., Shi, L., and Xu, F. (2014). A high-density genetic map identifies a novel major QTL for boron efficiency in oilseed rape (Brassica napus L.). PLoS ONE 9:e112089. doi: 10.1371/journal.pone.0112089

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Zhang, Q., Cheng, T., Yang, W., Pan, H., Zhong, J., et al. (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, 1–9. 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

Zhao, Z., Gu, H., Sheng, X., Yu, H., Wang, J., Cao, J., et al. (2012). Genetic purity testing of loose-curd cauliflower hybrids using SSR markers and grow out test. Seed Sci. Technol. 40, 209–214. doi: 10.15258/sst.2012.40.2.06

CrossRef Full Text | Google Scholar

Zhao, Z., Gu, H., Sheng, X., Yu, H., Wang, J., and Cao, J. (2014a). Advances and applications in crop quantitative trait loci. J. Nucl. Agric. Sci. 28, 1615–1624. doi: 10.11869/j.issn.100-8551.2014.09.1615

CrossRef Full Text

Zhao, Z., Gu, H., Sheng, X., Yu, H., Wang, J., Zhao, J., et al. (2014b). Genetic diversity and relationships among loose-curd cauliflower and related varieties as revealed by microsatellite markers. Sci. Hortic. 166, 105–110. doi: 10.1016/j.scienta.2013.12.024

CrossRef Full Text | Google Scholar

Zhao, Z., Gu, H., Wang, J., Sheng, X., and Yu, H. (2013). Development and comparison of quantitative methods to evaluate the curd solidity of cauliflower. J. Food Eng. 119, 477–482. doi: 10.1016/j.jfoodeng.2013.06.025

CrossRef Full Text | Google Scholar

Zhou, X., Xia, Y., Ren, X., Chen, Y., Huang, L., Huang, S., et al. (2014). Construction of a SNP-based genetic linkage map in cultivated peanut based on large scale marker development using next-generation double-digest restriction-site-associated DNA sequencing (ddRADseq). BMC Genomics 15:351. doi: 10.1186/1471-2164-15-351

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: cauliflower, SLAF, SNP, sequencing, genetic map

Citation: Zhao Z, Gu H, Sheng X, Yu H, Wang J, Huang L and Wang D (2016) Genome-Wide Single-Nucleotide Polymorphisms Discovery and High-Density Genetic Map Construction in Cauliflower Using Specific-Locus Amplified Fragment Sequencing. Front. Plant Sci. 7:334. doi: 10.3389/fpls.2016.00334

Received: 18 October 2015; Accepted: 04 March 2016;
Published: 21 March 2016.

Edited by:

Diego Rubiales, Consejo Superior de Investigaciones Científicas, Spain

Reviewed by:

Daniela Marone, Centre of Cereal Research - CREA-CER - Foggia, Italy
Michael Benjamin Kantar, University of Minneosta, USA

Copyright © 2016 Zhao, Gu, Sheng, Yu, Wang, Huang and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Honghui Gu, guhh@mail.zaas.ac.cn