Discovery of Pod Shatter-Resistant Associated SNPs by Deep Sequencing of a Representative Library Followed by Bulk Segregant Analysis in Rapeseed

Abstract

Background

Single nucleotide polymorphisms (SNPs) are an important class of genetic marker for target gene mapping. As of yet, there is no rapid and effective method to identify SNPs linked with agronomic traits in rapeseed and other crop species.

Methodology/Principal Findings

We demonstrate a novel method for identifying SNP markers in rapeseed by deep sequencing a representative library and performing bulk segregant analysis. With this method, SNPs associated with rapeseed pod shatter-resistance were discovered. Firstly, a reduced representation of the rapeseed genome was used. Genomic fragments ranging from 450–550 bp were prepared from the susceptible bulk (ten F2 plants with the silique shattering resistance index, SSRI <0.10) and the resistance bulk (ten F2 plants with SSRI >0.90), and also Solexa sequencing-produced 90 bp reads. Approximately 50 million of these sequence reads were assembled into contigs to a depth of 20-fold coverage. Secondly, 60,396 ‘simple SNPs’ were identified, and the statistical significance was evaluated using Fisher's exact test. There were 70 associated SNPs whose –log₁₀p value over 16 were selected to be further analyzed. The distribution of these SNPs appeared a tight cluster, which consisted of 14 associated SNPs within a 396 kb region on chromosome A09. Our evidence indicates that this region contains a major quantitative trait locus (QTL). Finally, two associated SNPs from this region were mapped on a major QTL region.

Conclusions/Significance

70 associated SNPs were discovered and a major QTL for rapeseed pod shatter-resistance was found on chromosome A09 using our novel method. The associated SNP markers were used for mapping of the QTL, and may be useful for improving pod shatter-resistance in rapeseed through marker-assisted selection and map-based cloning. This approach will accelerate the discovery of major QTLs and the cloning of functional genes for important agronomic traits in rapeseed and other crop species.

About the Authors:

Zhiyong Hu

Contributed equally to this work with: Zhiyong Hu, Wei Hua

Affiliation: Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture, Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Wuhan, Hubei, People's Republic of China

Wei Hua

Contributed equally to this work with: Zhiyong Hu, Wei Hua

Affiliation: Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture, Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Wuhan, Hubei, People's Republic of China

Shunmou Huang

Affiliation: Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture, Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Wuhan, Hubei, People's Republic of China

Hongli Yang

Affiliation: Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture, Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Wuhan, Hubei, People's Republic of China

Gaomiao Zhan

Affiliation: Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture, Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Wuhan, Hubei, People's Republic of China

Xinfa Wang

Affiliation: Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture, Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Wuhan, Hubei, People's Republic of China

Guihua Liu

Affiliation: Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture, Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Wuhan, Hubei, People's Republic of China

Hanzhong Wang

* E-mail: [email protected]

Affiliation: Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture, Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences, Wuhan, Hubei, People's Republic of China

Introduction

Single nucleotide polymorphisms (SNPs) are the most abundant variations and currently an important class of genetic marker distributed throughout the genome [1]. SNPs are imperative to gene fine mapping and association studies aimed at identifying alleles that affect agronomic traits [2].

In the past decade, SNPs have been used as genetic markers in many cultivated plants such as rice [3]–[5], maize [6]–[8] and rapeseed [9]–[10]. Genotyping of SNPs has improved, making genome-wide linkage analysis and molecular breeding more rapid and efficient [5].

Bulk segregant analysis (BSA) is efficient for selective genotyping wherein DNA samples of extreme individuals from each tail of a phenotypic distribution are pooled and the two resultant bulks need to be genotyped [11]. Markers linked to a QTL affecting a particular trait are expected to be present at different frequencies in the contrasting tails, resulting in polymorphic expression of genotype signals between the two bulks [12]. BSA saves cost by only genotyping the pooled DNA from individuals with similar phenotypes. It is an effective technique in that it detects large effect QTL alleles in a large sample of progenies at a relatively low cost. Because it saves time and cost, it has been widely used for the genetic analysis of qualitative traits [12]–[14].

Reduced representation libraries (RRLs) and pyrosequencing technologies have facilitated the high throughput discovery of SNPs [15]–[19]. RRLs can be used to identify SNPs that are useful in mapping the genome for haplotypes associated with disease [15]. Genome complexity is reduced by using restriction digests and size selection to build an RRL [20]. The use of fragments from a size-selected digestion permits a similar subset of fragments to be obtained from different genotypes that can be deep-sequenced for accurate SNP discovery.

Rapeseed (Brassica napus; genome AACC, 2n = 38) is an allopolyploid species that originated in a limited geographic region via spontaneous hybridizations between turnip rape (B. rapa; AA, 2n = 20) and cabbage (B. oleracea; CC, 2n = 18) genotypes [21]. Rapeseed is the most important oilseed crop grown in temperate agricultural regions and has the second highest output of all oil crops in the world, providing approximately 13% of the world's supply of vegetable oil [22], [23]. Rapeseed, as a commercial oil crop in Asia, Europe, North America and Australia, has attracted interest due to its utilization in food and feed production and its growing economic importance as a novel source of renewable energy, mainly as biodiesel [24]. Rapeseed, however, is subject to significant losses in seed production as a result of premature dehiscence before and during harvest [25]. Seed shedding from pods is commonly referred to as pod shatter, which is a major production risk in rapeseed worldwide [26], [27]. Mapping quantitative trait loci (QTLs) linked to pod shatter-resistant traits will assist in developing rapeseed cultivars suitable for harvesting machinery through a molecular marker-assisted selection (MAS) strategy. Developing a large set of SNP markers for genome analyses in rapeseed will facilitate fine mapping of QTL regions and will improve the identification and exploitation of genes affecting important traits. Additionally, these makers will also allow for selective breeding through genomic selection.

In this study, our objective was to use the RRL sequencing method, via BSA, to identify SNP markers linked to rapeseed pod shatter-resistance and utilize these markers to map the QTL regions. Our study identified 70 associated SNPs, which were selected from 60,396 ‘simple SNPs’, and a major QTL was found. The QTL may be a useful tool for promoting pod shatter-resistant rapeseed through MAS and map-based cloning. Our approach will accelerate the discovery of major QTLs and the cloning of functional genes for important agronomic traits in rapeseed and many other crop species.

Results

Construction of RRLs

The silique shattering resistance index (SSRI) of two rapeseed lines (B. napus cv. zy72360 and B. napus cv. R1) and their F2 population progenies (including 276 individuals) were identified with a random impact test designed for screening pod shatter-resistance (Table S1). The SSRI values of zy72360 (resistance to pod shattering) and R1 (very susceptible germplasm to pod shattering) were 0.94 and 0.06, respectively. The SSRI of their F2 population progenies was continuous distribution ranging from 0.00 to 1.00, but didn't display a normal distribution. The frequency distribution of SSRI of the F2 population including 225 randomly selected individuals showed three peaks which followed 1∶2∶1 segregation ratio by chi-squared test (Fig. S1). We estimated there maybe exited one major gene controlling the pod shatter-resistance. Ten F2 individuals, which are susceptible to pod shattering (SSRI <0.10), and ten F2 individuals, which are resistance to pod shattering (SSRI >0.90), were pooled to form the susceptible bulk (SK) and the resistance bulk (RK), respectively.

To identify the most suitable enzymes for RRL construction, we analyzed all of our blunt end enzymes. Rapeseed (B.napus; genome AACC, 2n = 38) is a recent allopolyploid species that originated in a limited geographic region through spontaneous hybridizations between turnip rape (B. rapa; AA, 2n = 20) and cabbage (B. oleracea; CC, 2n = 18) [21]. For this reason, the combined genomes, including all sequences (790,587,732 bp) of B. rapa [28] and B. oleracea (unpublished), were used as reference genomes. Firstly, the genome sequences of B. rapa and B. oleracea were combined and the combined sequences were then simulated to be digested with all of our blunt end enzymes. In the end, four blunt end enzymes, AluI, DpnI, RsaI and RseI, were selected with the goal of minimizing repetitive content in the target size range of 450–550 bp and interpreted to be suitable for RRL construction. By digesting with these enzymes, we were assured that the size of the DNA fragments would range from 450 to 550 and account for about 10% of the size of the whole sequences (Figure 1). After validation by using megabase-size rapeseed genomic DNA, we found that the restriction enzymes DpnI, RsaI and RseI did not entirely digest the rapeseed DNA (Figure 2) and were therefore not considered for further procedures. AluI, an enzyme that completely digests rapeseed DNA (Figure 2), was ultimately selected for construction of the library with the target size range of 450–550 bp from the genome and to obtain a suitable quantity for solexa sequencing.

[Figure omitted. See PDF.]

Figure 1. Frequency distribution of the rapeseed genome sequences.

Here we show the rapeseed genome sequence after digestion with fourblunt end enzymes, AluI, DpnI, RsaI and RseI, respectively (computer simulated).

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

[Figure omitted. See PDF.]

Figure 2. Gel image of digested rapeseed DNA.

Rapeseed DNA was digested with fourblunt end enzymes, AluI, DpnI, RsaI and RseI, respectively.

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

Sequencing Results

We independently sequenced the two RRLs (SK and RK), and a total of 48,888,892 reads (90 bp) were obtained (Table 1). Because we prepared the libraries using AluI, properly ligated fragments included a two-nucleotide CT tag at each restriction site. We discarded the 1,031,506 sequences that did not begin with this tag, leaving ∼ 48 million reads for analysis. Because rapeseed genome has not been pieced together, the combined genomes, including all sequences (790,587,732 bp) of B. rapa and B. oleracea, were used as reference genomes. Subsequently, the sequences were mapped to the reference genomes using SOAP software. We discarded reads that matched exactly with more than one position in the reference genomes, and the unique sequences were chosen for SNP discovery. The most frequently occurring sequence (798,774 occurrences) was a perfect match to more than 100 locations in the reference genomes and represented a contaminating repetitive element. After removing 1,102,886 sequences with low-quality scores, 47,786,006 reads (4.3 billion bp) corresponding to 25,048,721 (12,633,495 and 12,415,226 from RK and SK libraries, respectively) unique sequences remained for SNP discovery.

[Figure omitted. See PDF.]

Table 1. Summary of DNA sequence production and filtering.

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

The sequences obtained from SK and RK libraries covered 14% (114,308,963/790,587,732) and 11% (84,456,839/790,587,732) of the reference genomes, respectively. The overlap of the two libraries covered 8% (66,622,234/790,587,732) of the reference genomes and achieved a 20-fold average coverage depth, and that coverage depth was sufficient for SNP identification.

Associated SNP Discovery

The unique sequences produced ‘simple SNPs’ between two bulks by custom perl scripts [9] and 60,396 ‘simple SNPs’ were obtained in total. The statistical significance level was evaluated using Fisher's exact test [29]. Our results indicated that there were 25,226, 3,573, 450 and 70 ‘simple SNPs’remained when the value of –log₁₀pwas over 2, 5, 10 and 16, respectively. In order to reduce the number of associated SNPs in our study, only those putative associated SNPs with the value of –log₁₀pover 16 were selected for further analysis. Thus, these 70 SNPs associated with rapeseed pod shatter-resistance were obtained and listed in Table S2.

Mapping and Distributionof Associated SNPs

Using the genome sequences from B. rapa and B. oleracea, 70 associated SNPs were mapped to the reference genomes using SOAP software (Table S2) [30]. The map of the SNP distribution within the rapeseed genome was drawn using Microsoft Excel software (Figure S2). 70 associated SNPs were distributed on 17 rapeseed chromosomes. The map reveals the appearance of a tightly clustered SNP region which included 37% (26 SNPs) of the associated SNPs on chromosome A09 (Figure S2, Table S2). Furthermore, the location of 20% (14 SNPs) associated SNPs were located within a 396 kb sequence on chromosome A09 (Figure S2, Table S2). All the other SNPs did not appear in an obvious cluster. Because these SNPs were discovered using the two libraries, RK and SK, with large phenotype differences in rapeseed pod shatter-resistant traits, the tightly clustered region of associated SNPs implies that there likely exists a QTL for rapeseed pod shatter-resistance on chromosome A09.

QTL Location and Associated SNP Validation

In this study, we simultaneously constructed a genetic map by using the F2 population derived from the cross of zy72360 and R1. The SSRI of 276 individual pods from the F2 population was measured by random impact test according to the previously described method [31] with minor modifications. The QTL was scanned by WinQTLCart software [32] and one QTL, psr1, which influenced pod shatter-resistance was detected (Fig. 3, Fig. S3). The QTL allele, psr1, is located just on chromosome A09 and its genetic contribution rate is 47%. This result is consistent with that from associated SNPs and indicates that the associated SNP clustered region was indeed convenient for marking the location of the QTL.

[Figure omitted. See PDF.]

Figure 3. Genetic linkage map of chromosome A09 from rapeseed.

The numerals in left column indicate the genetic distance [centiMorgan (cM)]. SNP or SSR marker names are indicated in the right column. Psr1 indicates the QTL region of pod shatter-resistance.

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

In order to further validate if the location of the tightly clustered SNP region was identical to the QTL allele, psr1, 14 primer pairs were designed and synthesized from the 14 SNPs within the tightly clustered SNP region on chromosome A09 (Table S3). After polymorphic SNP primers were designed from the SNP flanking sequences of the reference genomes, amplicons were generated from the F2 population (derived from the cross of zy72360 and R1), genotyped using non-denaturing polyacrylamide gels and validated for sequence variations. Surprisingly, two associated SNP markers, NS380 and NS381 were mapped and traced to the region of the QTL allele, psr1 (Fig. 3).

Discussion

In this study we chose rapeseed pod shatter-resistance as the target trait, and succeeded in isolating the associated SNP markers and defining a major QTL region by using a method that combines BSA, RRL and sequencing. Our results are proof of its easy applicability and of its validity, indicating that it is a rapid and low-cost method for the discovery of associated SNP markers. BSA was first developed by Michelmore et al. [11] as a method for rapidly identifying markers linked to any specific gene or genomic region. Feng et al. [33] identified three markers linked to the or gene in B. rapa by BSA. Xu et al. [34] and Huang et al. [35] employed AFLP technology combined with BSA to identify markers linked to the seed color gene in B. juncea and the BnMs3 gene in B. napus, respectively. By applying high-throughput sequencing to BSA, Wenger et al. [36] revealed a novel xylose utilization gene from Saccharomyces cerevisiae. BSA has also been used for the fine mapping and for shortening the confidence intervals of QTL regions [13], [37]. RRL coupled with next generation sequencing technology allows for the detection of a large number of SNPs at a much lower time investment and expense than by conventional methods, and is used widely in animal and plant species such as cattle [18], swine [19], rainbow trout [2] and soybean [20].

The key technological aspects of our study include the selection of a versatile enzyme for RRL construction and the preferences for associated SNP selection. The selection of an optimal blunt end enzyme is the most important for RRL construction. In animal studies, the restriction enzyme HaeIII was usually selected for RRL construction in cattle [18], swine [19] and rainbow trout [2]. However, in plant studies, Hyten et al [20] reported that four different blunt end restriction enzyme combinations, each containing five different restriction enzymes, were tested as an effort to reduce the likelihood of repetitive sequences. The authors created a soybean RRL that consisted of DNA digested with a combination of HaeIII. PsiI, SspI, RsaI and MslI. We simulated the digestion of the masked rapeseed sequences on a computer using all of our blunt end enzymes. From this analysis, we decided to use AluI due to its excellent in silico digestion characteristics. Our results indicate that AluI is optimal for the construction of B. napus RRLs, and this was evident because almost 98% of our sequence reads began with the ‘CT’ sequence (Table 1).

In this study, the use of RRL, coupled with next generation sequencing and BSA, provided a powerful method for the discovery of SNP markers linked to agronomic traits.The advantage of our method is that it is rapid and requires a low cost to obtain associated SNPs and find major QTLs. Schneeberger et al. also report a similar method, SHOREmap, in which they identify a causative mutation in only 8 working days of hands-on time, and can also be applied to mapping of a QTL. Their method, however, is very expensive as it used a single genomic DNA sample to prepare an Illumina library, which was sequenced to 22-fold genome coverage [38]. In our method, the sequences used were only about 10% of the size of the whole sequences of the RRLs, and the costs were greatly reduced. In this study, we only analyzed sequences of about 500bp, the target size range for RRL regions. The larger the target size of the RRL, the less sequences were needed for adequate coverage. If the target size range for the RRL was increased from 500 bp to 2 kb, then only about 25% of the sequencing power was needed to reach the same average coverage depth, and the costs were reduced by at least a half. For later studies in which we look at other species or other traits, we will consider different target size combinations depending on the enzyme that is used for RRL construction and on the genome size of the particular species in question. One major disadvantage of our method is that it cannot be precisely used to evaluate the genetic contribution of each QTL region. In our study, we only identified one major QTL for pod shatter-resistance and therefore we only discovered one tightly clustered SNP region. In theory, for traits containing several QTLs, the degree of QTL contribution is suggested by the frequency of the associated markers.

Fully mature rapeseed pods are extremely sensitive to opening, resulting in seed loss. This loss is typically 8–12% of the total seed yield [39], and can even exceed 20% if harvesting is delayed beyond the optimum time or under adverse conditions [26]. Thus, rapeseed pod shatter-resistance is a very important agronomic trait. In this study, a total of 70 SNPs associated with rapeseed pod shatter-resistance were identified. The most tightly clustered SNP region, which contains 14 associated SNPs within a 396 kb stretch on chromosome A09, indicates and validates the existence of the large effect QTL allele, psr1. The SNP markers closely linked to pod shatter-resistance were applied to fine mapping of the QTL and are useful for pod shatter-resistant improvement in rapeseed through MAS and map-based cloning. Based on its quick application and low cost, our approach can be extensively applied in discovering major QTLs and will accelerate the cloning of functional genes for important agronomic traits in rapeseed and other crop species.

Materials and Methods

Ethics Statement

No specific permits were required for the described field studies. No specific permissions were required for these locations/activities. The location is not privately-owned or protected in any way. The field studies did not involve endangered or protected species.

Plant Materials

Two rapeseed lines, B. napus cv. zy72360 (resistance to pod shattering) and B. napus cv. R1 (very susceptible germplasm to pod shattering) were selected in this study. An F2 population that included 276 individuals derived from the cross of zy72360 and R1 was used. All plant materials were planted in a field located in Wuhan, China. All lines were planted in fields and each condition duplicated.

Random Impact Test

Random impact test was conducted according to the method in [31] with minor modifications. Ten replicates of 20 intact mature pods were placed in a closed polythene drum of 10 cm diameter, 14 cm height with fifty steel balls of 8 mm diameter. The drum was then shaken mechanically in a rocking bed (HQ45Z, Wuhan Scientific Instrument Factory, China) at 300 rpm for 1 min. The drum was opened after each replicate and the number of broken and damaged pods was counted. The SSRI was counted according to the following expressions:Where x_i is the number of broken and damaged silique at each replicate, i. Each accession was tested twice and an average value was calculated.

DNA Isolation and Pooling of Samples for BSA

Genomic DNA was isolated from fresh leaves using cetyl trimethyl ammonium bromide. The concentration of DNA was adjusted to 50 ng/µL for PCR amplification. Megabase-size genomic DNA was isolated from bulk leaves of F2 individuals following the procedure described by Zhang et al. [39] with minor modifications. Ten F2 plants, which are susceptible to pod shattering (SSRI <0.10), and ten F2 plants, which are resistance to pod shattering (SSRI >0.90), were pooled to form the SK and RK bulks, respectively.

RRL Construction and Sequencing

Four different restriction enzymes (AluI, DpnI, RsaI and RseI), producing fragments with blunt ends, were tested for library construction. Digestion of DNA was performed for each enzyme as suggested by the manufacturer (New England Bio Labs) and was digested overnight to ensure complete digestion. The digested DNA was then run on a 2% agarose gel and the digestion products were excised from the agarose gel at the 450 to 550 bp region. The QIAquick Gel Extraction Kit (Qiagen, Hilden, Germany) was used following the manufacturer's protocol to obtain a total of 5,000 ng of size-selected DNA. The blunt-end DNA fragments were given to BGI (Beijin Genome Institute, Shenzhen, China) for sequencing on the HiSeq2000 (Illumina, Inc; San Diego, CA, USA) after adding a 10 bp index. The sequence data obtained from BGI contained 90 bp sequence tags in which every base was given a quality score that is similar to a Phred score. We screened the data for overall quality by looking for an average base quality score of at least 25 across the entire tag for all reads according to a previously described method [18]. The Illumina sequence data have been deposited in the NCBI, Sequence Read Archive [GenBank: SRA045576].

Associated SNP Marker Discovery

Rapeseed is a recent allopolyploid species that originated in a limited geographic region through spontaneous hybridizations between B. rapa and B. oleracea genotypes [21]. In this study, a total of 790,587,732 bp sequences which include the combined genomic sequence of B. rapa [28] and B. oleracea (unpublished) were used as reference genomes. The sequences from the two bulks were mapped to the references using SOAP software [30]. The results produced ‘simple SNPs’ between the two bulks by custom perl scripts [9]. The statistical significance level was evaluated using Fisher's exact test [29].

Linkage Analysis

The F2 population, consisting of 276 individuals, was utilized to construct a genetic map and validate the associated SNPs. The primers sequences which were used for the genetic map are listed in Table S4. PCR amplification was carried out in a solution of 20 μL containing 100 ng of genomic DNA, 6 pmol of each primer, 2.5 mmol/L of MgCl₂, 0.25 mmol/L of each dNTP, ×reaction buffer and 0.75 U Taq polymerase. Thermocycling was performed using a DNA Engine Peltier Thermal Cycler, with 30 cycles of 94°C for 30 s, 55°C for 30 s, 72°C for 1 min, and a final extension at 72°C for 5 min before holding at 10°C. The PCR products were separated on an 8% non-denaturing polyacrylamide gel, and were stained with AgNO₃ [40].

Linkage analysis and genetic map construction was done using JoinMap 3.0 [41]. The QTL was scanned using WinQTLCart software [32].

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

The frequency distribution of the silique shattering resistance index. The three dashed frames show three peaks.

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

(TIF)

Figure S2.

Distribution of associated SNPs detected among rapeseed SK and RK bulks. The star indicates the tightly clustered SNPs region.

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

(TIF)

Figure S3.

Plot of the major QTL for pod shatter-resistance on linkage group A9 in B. napus .

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

(TIF)

Table S1.

The silique shattering resistance index of two rapeseed lines ( B. napus cv. zy72360 and B. napus cv. R1) and their F2 population progenies (including 276 individuals).

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

(XLS)

Table S2.

The details of 70 SNPs associated with rapeseed pod shatter resistance from B. napus .

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

(XLSX)

Table S3.

Primers of 14 SNPs of the tightly clustered SNPs region on A09 chromosome.

https://doi.org/10.1371/journal.pone.0034253.s006

(XLSX)

Table S4.

Primers sequences of the genetic map on A09 chromosome.

https://doi.org/10.1371/journal.pone.0034253.s007

(XLSX)

Author Contributions

Conceived and designed the experiments: HW WH. Performed the experiments: ZH SH HY GZ XW GL. Analyzed the data: ZH SH. Wrote the paper: ZH WH.

Citation: Hu Z, Hua W, Huang S, Yang H, Zhan G, Wang X, et al. (2012) Discovery of Pod Shatter-Resistant Associated SNPs by Deep Sequencing of a Representative Library Followed by Bulk Segregant Analysis in Rapeseed. PLoS ONE 7(4): e34253. https://doi.org/10.1371/journal.pone.0034253

References

1. Liu N, Chen L, Wang S, Oh C, Zhao H (2005) Comparison of single nucleotide polymorphisms and microsatellites in inference of population structure. BMC Genetics 6: S26.N. LiuL. ChenS. WangC. OhH. Zhao2005Comparison of single nucleotide polymorphisms and microsatellites in inference of population structure.BMC Genetics6S26

2. Sánchez CC, Smith TP, Wiedmann RT, Vallejo RL, Salem M, et al. (2009) Single nucleotide polymorphism discovery in rainbow trout by deep sequencing of a reduced representation library. BMC Genomics 10: 559.CC SánchezTP SmithRT WiedmannRL VallejoM. Salem2009Single nucleotide polymorphism discovery in rainbow trout by deep sequencing of a reduced representation library.BMC Genomics10559

3. Feltus FA, Wan J, Schulze SR, Estill JC, Jiang N, et al. (2004) An SNP resource for rice genetics and breeding based on subspecies indica and japonica genome alignments. Genome Res 14: 1812–1819.FA FeltusJ. WanSR SchulzeJC EstillN. Jiang2004An SNP resource for rice genetics and breeding based on subspecies indica and japonica genome alignments.Genome Res1418121819

4. Monna L, Ohta R, Masuda H, Koike A, Minobe Y (2006) Genomewide searching of single-nucleotide polymorphisms among eight distantly and closely related rice cultivars (Oryza sativa L.) and a wild accession (Oryza rufipogon Griff.). DNA Res 13: 43–51.L. MonnaR. OhtaH. MasudaA. KoikeY. Minobe2006Genomewide searching of single-nucleotide polymorphisms among eight distantly and closely related rice cultivars (Oryza sativa L.) and a wild accession (Oryza rufipogon Griff.).DNA Res134351

5. Nasu S, Suzuki J, Ohta R, Hasegawa K, Yui R, et al. (2002) Search for and analysis of single nucleotide polymorphisms (SNPs) in rice (Oryza sativa, Oryza rufipogon) and establishment of SNP markers. DNA Res 9: 163–171.S. NasuJ. SuzukiR. OhtaK. HasegawaR. Yui2002Search for and analysis of single nucleotide polymorphisms (SNPs) in rice (Oryza sativa, Oryza rufipogon) and establishment of SNP markers.DNA Res9163171

6. Batley J, Barker G, O'Sullivan H, Edwards KJ, Edwards D (2003) Mining for single nucleotide polymorphisms and insertions/deletions in maize expressed sequence tag data. Plant Physiol 132: 84–91.J. BatleyG. BarkerH. O'SullivanKJ EdwardsD. Edwards2003Mining for single nucleotide polymorphisms and insertions/deletions in maize expressed sequence tag data.Plant Physiol1328491

7. Bhattramakki D, Dolan M, Hanafey M, Wineland R, Vaske D, et al. (2002) Insertion-deletion polymorphisms in 3′ regions of maize genes occur frequently and can be used as highly informative genetic markers. Plant Mol Biol 48: 539–547.D. BhattramakkiM. DolanM. HanafeyR. WinelandD. Vaske2002Insertion-deletion polymorphisms in 3′ regions of maize genes occur frequently and can be used as highly informative genetic markers.Plant Mol Biol48539547

8. Bi IV, McMullen MD, Sanchez-Villeda H, Schroeder S, Gardiner J, et al. (2006) Single nucleotide polymorphisms and insertion-deletions for genetic markers and anchoring the maize fingerprint contig physical map. Crop Sci 46: 12–21.IV BiMD McMullenH. Sanchez-VilledaS. SchroederJ. Gardiner2006Single nucleotide polymorphisms and insertion-deletions for genetic markers and anchoring the maize fingerprint contig physical map.Crop Sci461221

9. Trick M, Long Y, Meng J, Bancroft L (2009) Single nucleotide polymorphism (SNP) discovery in the polyploidy Brassica napus using Solexa transcriptome sequencing. Plant Biotechnol J 7: 334–346.M. TrickY. LongJ. MengL. Bancroft2009Single nucleotide polymorphism (SNP) discovery in the polyploidy Brassica napus using Solexa transcriptome sequencing.Plant Biotechnol J7334346

10. Westermeier P, Wenzel G, Mohler V (2009) Development and evaluation of single-nucleotide polymorphism markers in allotetraploid rapeseed (Brassica napus L.) Theor Appl Genet 119: 1301–1311.P. WestermeierG. WenzelV. Mohler2009Development and evaluation of single-nucleotide polymorphism markers in allotetraploid rapeseed (Brassica napus L.) Theor Appl Genet11913011311

11. Michelmore RW, Paranand I, Kessele RV (1991) Identification of markers linked to disease resistance genes bybulk segregant analysis: A rapid method to detect markers inspecific genome using segregant population. Proc Natl Acad Sci USA 88: 9828–9832.RW MichelmoreI. ParanandRV Kessele1991Identification of markers linked to disease resistance genes bybulk segregant analysis: A rapid method to detect markers inspecific genome using segregant population.Proc Natl Acad Sci USA8898289832

12. Venuprasad R, Dalid CO, Del Valle M, Zhao D, EspirituM , et al. (2009) Identification and characterizationof large-effect quantitative trait loci for grain yield underlowland drought stress in rice using bulk-segregant analysis.Theor Appl Genet 120: 177–190.R. VenuprasadCO DalidM. Del ValleD. ZhaoEspirituM2009Identification and characterizationof large-effect quantitative trait loci for grain yield underlowland drought stress in rice using bulk-segregant analysis.Theor Appl Genet120177190

13. Salunkhe AS, Poornima R, Prince KSJ, Kanagaraj P, Sheeba JA, et al. (2011) Fine mapping QTL for drought resistance traits in rice(Oryza sativa L.) using bulk segregant analysis. Mol Biotechnol. AS SalunkheR. PoornimaKSJ PrinceP. KanagarajJA Sheeba2011Fine mapping QTL for drought resistance traits in rice(Oryza sativa L.) using bulk segregant analysis.Mol BiotechnolDOI 10.1007/s12033-011-9382-x. DOI 10.1007/s12033-011-9382-x.

14. Sun Y, Wang J, Crouch JH, XuY (2010) Efficiency ofselective genotyping for genetic analysis of complex traits andpotential applications in crop improvement. Mol Breeding 26: 493–511.Y. SunJ. WangJH CrouchXuY2010Efficiency ofselective genotyping for genetic analysis of complex traits andpotential applications in crop improvement.MolBreeding 26493511

15. Altshuler D, Pollara VJ, Cowles CR, Van Etten WJ, Baldwin J, et al. (2000) An SNP map of the human genome generated by reduced representation shotgun sequencing. Nature 407: 513–516.D. AltshulerVJ PollaraCR CowlesWJ Van EttenJ. Baldwin2000An SNP map of the human genome generated by reduced representation shotgun sequencing.Nature407513516

16. Barbazuk WB, Emrich SJ, Chen HD, Li L, Schnable PS (2007) SNP discovery via 454 transcriptome sequencing. Plant J 51: 910–918.WB BarbazukSJ EmrichHD ChenL. LiPS Schnable2007SNP discovery via 454 transcriptome sequencing.Plant J51910918

17. Novaes E, Drost DR, Farmerie WG, Pappas GJ , Grattapaglia D, et al. (2008) High-throughput gene and SNP discovery in Eucalyptus grandis, an uncharacterized genome. BMC Genomics 9: 312.E. NovaesDR DrostWG FarmeriePappas GJ JrD. Grattapaglia2008High-throughput gene and SNP discovery in Eucalyptus grandis, an uncharacterized genome.BMC Genomics9312

18. Van Tassell CP, Smith TP, Matukumalli LK, Taylor JF, Schnabel RD, et al. (2008) SNP discovery and allele frequency estimation by deep sequencing of reduced representation libraries. Nat Methods 5: 247–252.CP Van TassellTP SmithLK MatukumalliJF TaylorRD Schnabel2008SNP discovery and allele frequency estimation by deep sequencing of reduced representation libraries.Nat Methods5247252

19. Wiedmann RT, Smith TP, Nonneman DJ (2008) SNP discovery in swine by reduced representation and high throughput pyrosequencing. BMC Genet 9: 81.RT WiedmannTP SmithDJ Nonneman2008SNP discovery in swine by reduced representation and high throughput pyrosequencing.BMC Genet981

20. Hyten DL, Cannon SB, Song Q, Weeks N, Fickus EW, et al. (2010) High-throughput SNP discovery through deep resequencing of a reduced representation library to anchor and orient scaffolds in the soybean whole genome sequence. BMC Genomics 11: 38.DL HytenSB CannonQ. SongN. WeeksEW Fickus2010High-throughput SNP discovery through deep resequencing of a reduced representation library to anchor and orient scaffolds in the soybean whole genome sequence.BMC Genomics1138

21. Hasan M, FriedtW , Pons-KühnemannJ , FreitagNM , LinkK , et al. (2008) Association of gene-linked SSR markers to seed glucosinolate content in oilseed rape (Brassica napus ssp. napus). Theor Appl Genet 116: 1035–1049.M. HasanFriedtWPons-KühnemannJFreitagNMLinkK2008Association of gene-linked SSR markers to seed glucosinolate content in oilseed rape (Brassica napus ssp. napus).Theor Appl Genet11610351049

22. Amar S, EckeW , Becker HC, MÖllers C (2008) QTL for phytosterol and sinapate ester content in Brassica napus L. collocate with the two erucic acid genes. Theor Appl Genet 116: 1051–1061.S. AmarEckeWHC BeckerC. MÖllers2008QTL for phytosterol and sinapate ester content in Brassica napus L. collocate with the two erucic acid genes.Theor Appl Genet11610511061

23. Hajduch M, Casteel JE, Hurrelmeyer KE, Song Z, Agrawal GK, et al. (2006) Proteomic analysis of seed filling in Brassica napus developmental characterization of metabolic isozymes using high-resolution two-dimensional gel ecectrophoresis. Plant Physiol 141: 32–46.M. HajduchJE CasteelKE HurrelmeyerZ. SongGK Agrawal2006Proteomic analysis of seed filling in Brassica napus developmental characterization of metabolic isozymes using high-resolution two-dimensional gel ecectrophoresis.Plant Physiol1413246

24. Mittelbach M, Gangl S (2001) Long storage stability of biodiesel made from rapeseed and used frying oil. J Am Oil Chem Soc 78: 573–577.M. MittelbachS. Gangl2001Long storage stability of biodiesel made from rapeseed and used frying oil.J Am Oil Chem Soc78573577

25. Morgan CL, Bruce DM, Child R, Ladbrooke ZL, Arthur AE (1998) Genetic variation for pod shatter resistance among lines of oilseed rape developed from synthetic B. napus. Field Crop Res 58: 153–165.CL MorganDM BruceR. ChildZL LadbrookeAE Arthur1998Genetic variation for pod shatter resistance among lines of oilseed rape developed from synthetic B. napus.Field Crop Res58153165

26. Price JS, Hobson RN, Neale MA, Bruce DM (1996) Seed losses on commercial harvesting of oilseed rape. J Agric Engng Res 65: 183–191.JS PriceRN HobsonMA NealeDM Bruce1996Seed losses on commercial harvesting of oilseed rape.J Agric Engng Res65183191

27. Wang R, Ripley VL, Rakow G (2007) Pod shatter resistance evaluation in cultivars and breeding lines of Brassica napus, B. juncea and Sinapis alba. Plant Breeding 126: 588–595.R. WangVL RipleyG. Rakow2007Pod shatter resistance evaluation in cultivars and breeding lines of Brassica napus, B. juncea and Sinapis alba.Plant Breeding126588595

28. Wang X, Wang H, Wang J, Sun R, Wu J, et al. (2011) The genome of the mesopolyploid crop species Brassica rapa. Nature genetics 43: 1035–1039.X. WangH. WangJ. WangR. SunJ. Wu2011The genome of the mesopolyploid crop species Brassica rapa.Nature genetics4310351039

29. Need AC, Ge D, Weale ME, Maia J, Feng S, et al. (2009) A genome-wide investigation of SNPs and CNVs in Schizophrenia. Plos Genet 5: e1000373.AC NeedD. GeME WealeJ. MaiaS. Feng2009A genome-wide investigation of SNPs and CNVs in Schizophrenia.Plos Genet5e1000373

30. Li R, Yu C, Li Y, Lam T, Yiu S, et al. (2009) SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics 25: 1966–1967.R. LiC. YuY. LiT. LamS. Yiu2009SOAP2: an improved ultrafast tool for short read alignment.Bioinformatics2519661967

31. WenYC , Fu TD, Tu JX, Ma CZ, Shen JX, et al. (2008) Screening and analysis of resistance to silique shattering in rape (Brassica napus L.). Acta Agronomica Sinica 34: pp. 163–166. .WenYCTD FuJX TuCZ MaJX Shen2008Screening and analysis of resistance to silique shattering in rape (Brassica napus L.).Acta Agronomica Sinica 34163166. In Chinese

32. Wang S, Basten CJ, Zeng ZB (2011) Windows QTL Cartographer 2.5. Department of Statistics, North Carolina State University, Raleigh, NC. Statistival genetics and bioinformation website. 1: S. WangCJ BastenZB Zeng2011Windows QTL Cartographer 2.5. Department of Statistics, North Carolina State University, Raleigh, NC. Statistival genetics and bioinformation website.1Available: http://statgen.ncsu.edu/qtlcart/WQTLCart.htm. Accessed 2012 Mar. Available: http://statgen.ncsu.edu/qtlcart/WQTLCart.htm. Accessed 2012 Mar.

33. Feng H, Li Y, Liu Z, Liu J (2010) Mapping of or, a gene conferring orange color on the innerleaf of the Chinese cabbage (Brassica rapa L. ssp. pekinensis). Mol Breeding. H. FengY. LiZ. LiuJ. Liu2010Mapping of or, a gene conferring orange color on the innerleaf of the Chinese cabbage (Brassica rapa L. ssp.pekinensis). Mol BreedingDOI 10.1007/s11032-010-9542-x. DOI 10.1007/s11032-010-9542-x.

34. Xu A, Huang Z, Ma C, Xiao E, Tian G, et al. (2010) Inheritance of seed color and molecular markers linkedto the seed color gene in Brassica juncea. Mol Breeding 25: 57–65.A. XuZ. HuangC. MaE. XiaoG. Tian2010Inheritance of seed color and molecular markers linkedto the seed color gene in Brassica juncea.Mol Breeding255765

35. Huang Z, Xiao L, Dun X, Xia S, Yi B, et al. (2010) Improvement of the recessive genic male sterile lineswith a subgenomic background in Brassica napusby molecular marker-assisted selection. Mol Breeding. Z. HuangL. XiaoX. DunS. XiaB. Yi2010Improvement of the recessive genic male sterile lineswith a subgenomic background in Brassica napusby molecular marker-assisted selection.Mol BreedingDOI:10.1007/s11032-010-9536-8. DOI:10.1007/s11032-010-9536-8.

36. Wenger JW, Schwartz K, Sherlock G (2010) Bulk segregant analysis by high-throughput sequencing reveals a novel xylose utilization gene fromSaccharomyces cerevisiae. Plos Genet, 6: e1000942.JW WengerK. SchwartzG. Sherlock2010Bulk segregant analysis by high-throughput sequencing reveals a novel xylose utilization gene fromSaccharomyces cerevisiae.Plos Genet,6e1000942

37. Li W, Wu J, Weng S, Zhang D, Zhang Y, et al. (2010) Characterization and fine mapping of the glabrous leafand hull mutants (gl1) in rice (Oryza sativa L.) Plant Cell Rep 29: 617–627.W. LiJ. WuS. WengD. ZhangY. Zhang2010Characterization and fine mapping of the glabrous leafand hull mutants (gl1) in rice (Oryza sativa L.) Plant Cell Rep29617627

38. Kadkol GP, Macmillan RH, Burrow RP, Halloran GM (1984) Evaluation of Brassica genotypes for resistance to shatter I. Development of a laboratory test. Euphytica 33: 63–73.GP KadkolRH MacmillanRP BurrowGM Halloran1984Evaluation of Brassica genotypes for resistance to shatter I. Development of a laboratory test.Euphytica336373

39. Zhang H, Zhao X, Ding X, Paterson AH, Wing R (1995) Preparation of Megabase-size DNA from Plant nuclei. Plant J 7: 175–184.H. ZhangX. ZhaoX. DingAH PatersonR. Wing1995Preparation of Megabase-size DNA from Plant nuclei.Plant J7175184

40. Martin-Lopes P, Zhang H, Koebner R (2001) Detection of single nucleotide mutations in wheat using single strand conformation polymorphism gels. Plant Mol Biol Rep 19: 159–162.P. Martin-LopesH. ZhangR. Koebner2001Detection of single nucleotide mutations in wheat using single strand conformation polymorphism gels.Plant Mol Biol Rep19159162

41. Van Ooijen JW, Voorrips RE (2001) JoinMap3.0, Softwarefor the calculation of genetic linkage maps. PlantResearch International, Wageningen, the Netherlands. JW Van OoijenRE Voorrips2001JoinMap3.0, Softwarefor the calculation of genetic linkage maps.PlantResearch International, Wageningen, the Netherlands