QTL Mapping of Leafy Heads by Genome Resequencing in the RIL Population of Brassica rapa

Leaf heads of cabbage (Brassica oleracea), Chinese cabbage (B. rapa), and lettuce (Lactuca sativa) are important vegetables that supply mineral nutrients, crude fiber and vitamins in the human diet. Head size, head shape, head weight, and heading time contribute to yield and quality. In an attempt to investigate genetic basis of leafy head in Chinese cabbage (B. rapa), we took advantage of recent technical advances of genome resequencing to perform quantitative trait locus (QTL) mapping using 150 recombinant inbred lines (RILs) derived from the cross between heading and non-heading Chinese cabbage. The resequenced genomes of the parents uncovered more than 1 million SNPs. Genotyping of RILs using the high-quality SNPs assisted by Hidden Markov Model (HMM) generated a recombination map. The raw genetic map revealed some physical assembly error and missing fragments in the reference genome that reduced the quality of SNP genotyping. By deletion of the genetic markers in which recombination rates higher than 20%, we have obtained a high-quality genetic map with 2209 markers and detected 18 QTLs for 6 head traits, from which 3 candidate genes were selected. These QTLs provide the foundation for study of genetic basis of leafy heads and the other complex traits.     

Identification of Genome-Wide Variants and Discovery of Variants Associated with Brassica rapa Clubroot Resistance Gene Rcr1 through Bulked Segregant RNA Sequencing

Clubroot, caused by Plasmodiophora brassicae, is an important disease on Brassica species worldwide. A clubroot resistance gene, Rcr1, with efficacy against pathotype 3 of P. brassicae, was previously mapped to chromosome A03 of B. rapa in pak choy cultivar “Flower Nabana”. In the current study, resistance to pathotypes 2, 5 and 6 was shown to be associated with Rcr1 region on chromosome A03. Bulked segregant RNA sequencing was performed and short read sequences were assembled into 10 chromosomes of the B. rapa reference genome v1.5. For the resistant (R) bulks, a total of 351.8 million (M) sequences, 30,836.5 million bases (Mb) in length, produced 120-fold coverage of the reference genome. For the susceptible (S) bulks, 322.9 M sequences, 28,216.6 Mb in length, produced 109-fold coverage. In total, 776.2 K single nucleotide polymorphisms (SNPs) and 122.2 K insertion / deletion (InDels) in R bulks and 762.8 K SNPs and 118.7 K InDels in S bulks were identified; each chromosome had about 87% SNPs and 13% InDels, with 78% monomorphic and 22% polymorphic variants between the R and S bulks. Polymorphic variants on each chromosome were usually below 23%, but made up 34% of the variants on chromosome A03. There were 35 genes annotated in the Rcr1 target region and variants were identified in 21 genes. The numbers of poly variants differed significantly among the genes. Four out of them encode Toll-Interleukin-1 receptor / nucleotide-binding site / leucine-rich-repeat proteins; Bra019409 and Bra019410 harbored the higher numbers of polymorphic variants, which indicates that they are more likely candidates of Rcr1. Fourteen SNP markers in the target region were genotyped using the Kompetitive Allele Specific PCR method and were confirmed to associate with Rcr1. Selected SNP markers were analyzed with 26 recombinants obtained from a segregating population consisting of 1587 plants, indicating that they were completely linked to Rcr1. Nine SNP markers were used for marker-assisted introgression of Rcr1 into B. napus canola from B. rapa, with 100% accuracy in this study.     

Genome-wide comparative analysis of NBS-encoding genes between Brassica species and Arabidopsis thaliana

Background

Plant disease resistance (R) genes with the nucleotide binding site (NBS) play an important role in offering resistance to pathogens. The availability of complete genome sequences of Brassica oleracea and Brassica rapa provides an important opportunity for researchers to identify and characterize NBS-encoding R genes in Brassica species and to compare with analogues in Arabidopsis thaliana based on a comparative genomics approach. However, little is known about the evolutionary fate of NBS-encoding genes in the Brassica lineage after split from A. thaliana.

Results

Here we present genome-wide analysis of NBS-encoding genes in B. oleracea, B. rapa and A. thaliana. Through the employment of HMM search and manual curation, we identified 157, 206 and 167 NBS-encoding genes in B. oleracea, B. rapa and A. thaliana genomes, respectively. Phylogenetic analysis among 3 species classified NBS-encoding genes into 6 subgroups. Tandem duplication and whole genome triplication (WGT) analyses revealed that after WGT of the Brassica ancestor, NBS-encoding homologous gene pairs on triplicated regions in Brassica ancestor were deleted or lost quickly, but NBS-encoding genes in Brassica species experienced species-specific gene amplification by tandem duplication after divergence of B. rapa and B. oleracea. Expression profiling of NBS-encoding orthologous gene pairs indicated the differential expression pattern of retained orthologous gene copies in B. oleracea and B. rapa. Furthermore, evolutionary analysis of CNL type NBS-encoding orthologous gene pairs among 3 species suggested that orthologous genes in B. rapa species have undergone stronger negative selection than those in B .oleracea species. But for TNL type, there are no significant differences in the orthologous gene pairs between the two species.

Conclusion

This study is first identification and characterization of NBS-encoding genes in B. rapa and B. oleracea based on whole genome sequences. Through tandem duplication and whole genome triplication analysis in B. oleracea, B. rapa and A. thaliana genomes, our study provides insight into the evolutionary history of NBS-encoding genes after divergence of A. thaliana and the Brassica lineage. These results together with expression pattern analysis of NBS-encoding orthologous genes provide useful resource for functional characterization of these genes and genetic improvement of relevant crops.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-3) contains supplementary material, which is available to authorized users.     

A Complex Recombination Pattern in the Genome of Allotetraploid Brassica napus as Revealed by a High-Density Genetic Map

Polyploidy plays a crucial role in plant evolution. Brassica napus (2n = 38, AACC), the most important oil crop in the Brassica genus, is an allotetraploid that originated through natural doubling of chromosomes after the hybridization of its progenitor species, B. rapa (2n = 20, AA) and B. oleracea (2n = 18, CC). A better understanding of the evolutionary relationship between B. napus and B. rapa, B. oleracea, as well as Arabidopsis, which has a common ancestor with these three species, will provide valuable information about the generation and evolution of allopolyploidy. Based on a high-density genetic map with single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) markers, we performed a comparative genomic analysis of B. napus with Arabidopsis and its progenitor species B. rapa and B. oleracea. Based on the collinear relationship of B. rapa and B. oleracea in the B. napus genetic map, the B. napus genome was found to consist of 70.1% of the skeleton components of the chromosomes of B. rapa and B. oleracea, with 17.7% of sequences derived from reciprocal translocation between homoeologous chromosomes between the A- and C-genome and 3.6% of sequences derived from reciprocal translocation between non-homologous chromosomes at both intra- and inter-genomic levels. The current study thus provides insights into the formation and evolution of the allotetraploid B. napus genome, which will allow for more accurate transfer of genomic information from B. rapa, B. oleracea and Arabidopsis to B. napus.     

Genome Wide Analysis of Flowering Time Trait in Multiple Environments via High-Throughput Genotyping Technique in Brassica napus L.

The prediction of the flowering time (FT) trait in Brassica napus based on genome-wide markers and the detection of underlying genetic factors is important not only for oilseed producers around the world but also for the other crop industry in the rotation system in China. In previous studies the low density and mixture of biomarkers used obstructed genomic selection in B. napus and comprehensive mapping of FT related loci. In this study, a high-density genome-wide SNP set was genotyped from a double-haploid population of B. napus. We first performed genomic prediction of FT traits in B. napus using SNPs across the genome under ten environments of three geographic regions via eight existing genomic predictive models. The results showed that all the models achieved comparably high accuracies, verifying the feasibility of genomic prediction in B. napus. Next, we performed a large-scale mapping of FT related loci among three regions, and found 437 associated SNPs, some of which represented known FT genes, such as AP1 and PHYE. The genes tagged by the associated SNPs were enriched in biological processes involved in the formation of flowers. Epistasis analysis showed that significant interactions were found between detected loci, even among some known FT related genes. All the results showed that our large scale and high-density genotype data are of great practical and scientific values for B. napus. To our best knowledge, this is the first evaluation of genomic selection models in B. napus based on a high-density SNP dataset and large-scale mapping of FT loci.     

Genome-Wide Comparative Analysis of 20 Miniature Inverted-Repeat Transposable Element Families in Brassica rapa and B. oleracea

Miniature inverted-repeat transposable elements (MITEs) are ubiquitous, non-autonomous class II transposable elements. Here, we conducted genome-wide comparative analysis of 20 MITE families in B. rapa, B. oleracea, and Arabidopsis thaliana. A total of 5894 and 6026 MITE members belonging to the 20 families were found in the whole genome pseudo-chromosome sequences of B. rapa and B. oleracea, respectively. Meanwhile, only four of the 20 families, comprising 573 members, were identified in the Arabidopsis genome, indicating that most of the families were activated in the Brassica genus after divergence from Arabidopsis. Copy numbers varied from 4 to 1459 for each MITE family, and there was up to 6-fold variation between B. rapa and B. oleracea. In particular, analysis of intact members showed that whereas eleven families were present in similar copy numbers in B. rapa and B. oleracea, nine families showed copy number variation ranging from 2- to 16-fold. Four of those families (BraSto-3, BraTo-3, 4, 5) were more abundant in B. rapa, and the other five (BraSto-1, BraSto-4, BraTo-1, 7 and BraHAT-1) were more abundant in B. oleracea. Overall, 54% and 51% of the MITEs resided in or within 2 kb of a gene in the B. rapa and B. oleracea genomes, respectively. Notably, 92 MITEs were found within the CDS of annotated genes, suggesting that MITEs might play roles in diversification of genes in the recently triplicated Brassica genome. MITE insertion polymorphism (MIP) analysis of 289 MITE members showed that 52% and 23% were polymorphic at the inter- and intra-species levels, respectively, indicating that there has been recent MITE activity in the Brassica genome. These recently activated MITE families with abundant MIP will provide useful resources for molecular breeding and identification of novel functional genes arising from MITE insertion.     

Genome-wide discovery of DNA polymorphism in Brassica rapa

Single nucleotide polymorphisms (SNPs) and/or insertion/deletions (InDels) are frequent sequence variations in the plant genome, which can be developed as molecular markers for genetic studies on crop improvement. The ongoing Brassica rapa genome sequencing project has generated vast amounts of sequence data useful in genetic research. Here, we report a genome-wide survey of DNA polymorphisms in the B. rapa genome based on the 557 bacterial artificial clone sequences of B. rapa ssp. pekinensis cv. Chiifu. We identified and characterized 21,311 SNPs and 6,753 InDels in the gene space of the B. rapa genome by re-sequencing 1,398 sequence-tagged sites (STSs) in eight genotypes. Comparison of our findings with a B. rapa genetic linkage map confirmed that STS loci were distributed randomly over the B. rapa whole genome. In the 1.4 Mb of aligned sequences, mean nucleotide polymorphism and diversity were θ = 0.00890 and π = 0.00917, respectively. Additionally, the nucleotide diversity in introns was almost three times greater than that in exons, and the frequency of observed InDel was almost 17 times higher in introns than in exons. Information regarding SNPs/InDels obtained here will provide an important resource for genetic studies and breeding programs of B. rapa.     

A Brassica rapa Linkage Map of EST-based SNP Markers for Identification of Candidate Genes Controlling Flowering Time and Leaf Morphological Traits

For identification of genes responsible for varietal differences in flowering time and leaf morphological traits, we constructed a linkage map of Brassica rapa DNA markers including 170 EST-based markers, 12 SSR markers, and 59 BAC sequence-based markers, of which 151 are single nucleotide polymorphism (SNP) markers. By BLASTN, 223 markers were shown to have homologous regions in Arabidopsis thaliana, and these homologous loci covered nearly the whole genome of A. thaliana. Synteny analysis between B. rapa and A. thaliana revealed 33 large syntenic regions. Three quantitative trait loci (QTLs) for flowering time were detected. BrFLC1 and BrFLC2 were linked to the QTLs for bolting time, budding time, and flowering time. Three SNPs in the promoter, which may be the cause of low expression of BrFLC2 in the early-flowering parental line, were identified. For leaf lobe depth and leaf hairiness, one major QTL corresponding to a syntenic region containing GIBBERELLIN 20 OXIDASE 3 and one major QTL containing BrGL1, respectively, were detected. Analysis of nucleotide sequences and expression of these genes suggested possible involvement of these genes in leaf morphological traits.     

Comparative mitochondrial genome analysis reveals the evolutionary rearrangement mechanism in Brassica

The genus Brassica has many species that are important for oil, vegetable and other food products. Three mitochondrial genome types (mitotype) originated from its common ancestor. In this paper, a B. nigra mitochondrial main circle genome with 232,407 bp was generated through de novo assembly. Synteny analysis showed that the mitochondrial genomes of B. rapa and B. oleracea had a better syntenic relationship than B. nigra. Principal components analysis and development of a phylogenetic tree indicated maternal ancestors of three allotetraploid species in Us triangle of Brassica. Diversified mitotypes were found in allotetraploid B. napus, in which napus‐type B. napus was derived from B. oleracea, while polima‐type B. napus was inherited from B. rapa. In addition, the mitochondrial genome of napus‐type B. napus was closer to botrytis‐type than capitata‐type B. oleracea. The sub‐stoichiometric shifting of several mitochondrial genes suggested that mitochondrial genome rearrangement underwent evolutionary selection during domestication and/or plant breeding. Our findings clarify the role of diploid species in the maternal origin of allotetraploid species in Brassica and suggest the possibility of breeding selection of the mitochondrial genome.     

Chloroplast and mitochondrial SSR help to distinguish allo-cytoplasmic male sterile types in cabbage (Brassica oleracea L. var. capitata)

The profiles of single sequence repeat (SSR) in six distinct allo-cytoplasmic male sterile (CMS) types of cabbage (Brassica oleracea L. var. capitata) were generated using 32 SSR primer pairs derived from the Arabidopsis thaliana chloroplast (cp) genome and another 21 SSR primers from the B. napus mitochondrial (mt) genome sequences. In total, 11 cpSSR and 4 mtSSR primers revealed polymorphism among the six cabbage CMS types, namely NigCMS, OguCMSR₁, OguCMSR₂, OguCMSR₃, OguCMSHY and PolCMS. Through cluster analysis, six cabbage CMS types could be unambiguously differentiated with just three sets of primers (ACP43, ACP47, mtSSR2). Analysis of the selected amplicon sequences showed high identity to that of the corresponding sequences in A. thaliana, B. rapa and B. napus. The aligned cluster analysis revealed that the polymorphism mainly included SSR number variation, single nucleotide polymorphism (SNP), and sequence insertion or deletion (InDel). Our results demonstrated that specific mitochondrial or chloroplast SSR analysis could be a feasible alternative means for cabbage CMS type identification.     

Genotyping-by-sequencing map permits identification of clubroot resistance QTLs and revision of the reference genome assembly in cabbage (Brassica oleracea L.)

Clubroot is a devastating disease caused by Plasmodiophora brassicae and results in severe losses of yield and quality in Brassica crops. Many clubroot resistance genes and markers are available in Brassica rapa but less is known in Brassica oleracea. Here, we applied the genotyping-by-sequencing (GBS) technique to construct a high-resolution genetic map and identify clubroot resistance (CR) genes. A total of 43,821 SNPs were identified using GBS data for two parental lines, one resistant and one susceptible lines to clubroot, and 18,187 of them showed >5× coverage in the GBS data. Among those, 4,103 were credibly genotyped for all 78 F₂ individual plants. These markers were clustered into nine linkage groups spanning 879.9 cM with an average interval of 1.15 cM. Quantitative trait loci (QTLs) survey based on three rounds of clubroot resistance tests using F_{2 : 3} progenies revealed two and single major QTLs for Race 2 and Race 9 of P. brassicae, respectively. The QTLs show similar locations to the previously reported CR loci for Race 4 in B. oleracea but are in different positions from any of the CR loci found in B. rapa. We utilized two reference genome sequences in this study. The high-resolution genetic map developed herein allowed us to reposition 37 and 2 misanchored scaffolds in the 02–12 and TO1000DH genome sequences, respectively. Our data also support additional positioning of two unanchored 3.3 Mb scaffolds into the 02–12 genome sequence.     

Genetic and physical fine mapping of a multilocular gene Bjln1 in Brassica juncea to a 208-kb region

Most of the germplasm resources in Brassica juncea produce silique with only two locules, whereas a few varieties can produce silique with three or four locules. The increase in locule number in B. juncea has been shown to cause an increase in the number of seeds per silique, resulting in an increase in the yield per plant. Thus, the development of high-locule-number varieties may be an effective way of improving the yield of B. juncea. Duoshi, a B. juncea landrace originating from the Qinghai–Tibetan plateau, produces silique with 3–4 locules. Genetic analysis has shown that the high-locule-number trait in Duoshi is determined by two recessive genes, tentatively designated as Bjln1 and Bjln2. For fine mapping of the Bjln1 gene, a BC3 population was developed from the cross between Duoshi (multilocular parent) and Xinjie (bilocular parent). Using a combination of amplified fragment length polymorphism (AFLP) and bulked segregant analysis, only two AFLP markers linked to Bjln1 were identified. Preliminary linkage analysis showed that the two AFLP markers were located on the same side of Bjln1. Blast analysis revealed that the sequences of the two AFLP markers had homologues on Scaffold000019 at the bottom of B. rapa A7. Using the results of linkage analysis and BlastN searches, simple sequence repeat (SSR) markers were subsequently developed based on the sequence information from B. rapa A7. Seven SSR markers were eventually identified, of which ln 8 was co-segregated with Bjln1. ln 7 and ln 9, the closest flanking markers, were mapped at 2.0 and 0.4 cM distant from the Bjln1 gene, respectively. The SSR markers were cloned, sequenced and mapped on A7 of B. rapa (corresponding to J7 in the A genome of B. juncea). The two closest flanking markers, ln 7 and ln 9, were mapped within a 208-kb genomic region on B. rapa A7, in which the Bjln1 gene might be included. The present study may facilitate cloning of the Bjln1 gene as well as the selection process for developing multilocular varieties in B. juncea by marker-assisted selection and genetic engineering.     

Cloning of fatty acid elongase1 gene and molecular identification of A and C genome in Brassica species

The fatty acid elongase 1 (FAE1) genes of Brassic napus were cloned from two cultivars, i.e. Zhongshuan No. 9 with low erucic acid content, and Zhongyou 821 with high erucic acid content, using the degenerate PCR primers. The sequence analysis showed that there was no intron within the FAE1 genes. The FAE1 genes from Zhongyou 821 contained a coding sequence of 1521 nucleotides, and those cloned from Zhongshuan No. 9 contained a 1517 bp coding sequence. Alignment of the FAE1 sequences from Brassica rapa, B. oleracea and B. napus detected 31 single nucleotide polymorphic sites (2.03%), which resulted in 7 amino-acid substitutions. Further analysis indicated that 19 SNPs were genome-specific, of which, 95% were synonymous mutations. The nucleotide substitution at position 1217 in the FAE1 genes led to a specific site of restricted cleavage. An AvrII cleavage site was present only in the C genome genes and absent in the A genome FAE1 genes. Digestion profile of the FAE1 sequences from B. rapa, B. oleracea and B. napus produced with AvrII confirmed that the FAE1 genes of B. oleracea origin was recognized and digested, while that of B. rapa origin could not. The results indicated that by AvrII cleavage it was possible to distinguish B. rapa from B. oleracea and between the A and C genome of B. napus. In addition, the FAE1 genes could be used as marker genes to detect the pollen flow of B. napus, thus providing an alternative method for risk assessment of gene flow. Supported by the National Natural Science Foundation of China (Grant No. 30471099), Development Plan of the State Key Fundamental Research of China (Grant No. 2006CB101600), and the National High Technology and Development Program of China (Grant No. 2006AA10A113)     

Identification of genome-wide single nucleotide polymorphisms in allopolyploid crop Brassica napus

Background

Single nucleotide polymorphisms (SNPs) are the most common type of genetic variation. Identification of large numbers of SNPs is helpful for genetic diversity analysis, map-based cloning, genome-wide association analyses and marker-assisted breeding. Recently, identifying genome-wide SNPs in allopolyploid Brassica napus (rapeseed, canola) by resequencing many accessions has become feasible, due to the availability of reference genomes of Brassica rapa (2n = AA) and Brassica oleracea (2n = CC), which are the progenitor species of B. napus (2n = AACC). Although many SNPs in B. napus have been released, the objective in the present study was to produce a larger, more informative set of SNPs for large-scale and efficient genotypic screening. Hence, short-read genome sequencing was conducted on ten elite B. napus accessions for SNP discovery. A subset of these SNPs was randomly selected for sequence validation and for genotyping efficiency testing using the Illumina GoldenGate assay.

Results

A total of 892,536 bi-allelic SNPs were discovered throughout the B. napus genome. A total of 36,458 putative amino acid variants were located in 13,552 protein-coding genes, which were predicted to have enriched binding and catalytic activity as a result. Using the GoldenGate genotyping platform, 94 of 96 SNPs sampled could effectively distinguish genotypes of 130 lines from two mapping populations, with an average call rate of 92%.

Conclusions

Despite the polyploid nature of B. napus, nearly 900,000 simple SNPs were identified by whole genome resequencing. These SNPs were predicted to be effective in high-throughput genotyping assays (51% polymorphic SNPs, 92% average call rate using the GoldenGate assay, leading to an estimated >450 000 useful SNPs). Hence, the development of a much larger genotyping array of informative SNPs is feasible. SNPs identified in this study to cause non-synonymous amino acid substitutions can also be utilized to directly identify causal genes in association studies.     

The use of species-specific DNA markers for assessing alien chromosome transfer in Brassica rapa and Brassica oleracea-monosomic additions of Raphanus sativus

Monosomic addition lines (MALs) are useful materials not only for cytogenetic and molecular genetic studies but also for plant breeding as gene sources. In our previous study, two MALs in the tribe Brassiceae were developed, one being Raphanus sativus lines with alien chromosomes of Brassica rapa (B. rapa-monosomic addition lines; BrMALs) and the second being those with alien chromosomes of Brassica oleracea (B. oleracea-monosomic addition lines; BoMALs). We developed species-specific DNA markers from the genomic sequences of B. rapa and B. oleracea comparing them with those of R. sativus, and identified chromosomes added in BrMALs and BoMALs using these markers. It was revealed that eight types of BrMALs have seven chromosomes of B. rapa and seven types of BoMALs have six chromosomes of B. oleracea. Furthermore, chromosome breakage and homoeologous recombination were suggested to have occurred in some MALs. The developed species-specific DNA markers are considered to be useful for producing MALs and also for assessing chromosome abnormality in MALs.     

Development of a core set of single-locus SSR markers for allotetraploid rapeseed (Brassica napus L.)

Brassica napus (AACC) is a recent allotetraploid species evolved through hybridization between two diploids, B. rapa (AA) and B. oleracea (CC). Due to extensive genome duplication and homoeology within and between the A and C genomes of B. napus, most SSR markers display multiple fragments or loci, which limit their application in genetics and breeding studies of this economically important crop. In this study, we collected 3,890 SSR markers from previous studies and also developed 5,968 SSR markers from genomic sequences of B. rapa, B. oleracea and B. napus. Of these, 2,701 markers that produced single amplicons were putative single-locus markers in the B. napus genome. Finally, a set of 230 high-quality single-locus SSR markers were established and assigned to the 19 linkage groups of B. napus using a segregating population with 154 DH individuals. A subset of 78 selected single-locus SSR markers was proved to be highly stable and could successfully discriminate each of the 45 inbred lines and hybrids. In addition, most of the 230 SSR markers showed the single-locus nature in at least one of the Brassica species of the U’s triangle besides B. napus. These results indicated that this set of single-locus SSR markers has a wide range of coverage with excellent stability and would be useful for gene tagging, sequence scaffold assignment, comparative mapping, diversity analysis, variety identification and association mapping in Brassica species.     

A high-throughput SNP array in the amphidiploid species Brassica napus shows diversity in resistance genes

Single-nucleotide polymorphisms (SNPs)are molecular markers based on nucleotide variation and can be used for genotyping assays across populations and to track genomic inheritance. SNPs offer a comprehensive genotyping alternative to whole-genome sequencing for both agricultural and research purposes including molecular breeding and diagnostics, genome evolution and genetic diversity analyses, genetic mapping, and trait association studies. Here genomic SNPs were discovered between four cultivars of the important amphidiploid oilseed species Brassica napus and used to develop a B. napus Infinium? array containing 5,306 SNPs randomly dispersed across the genome. Assay success was high, with >94 % of these producing a reproducible, polymorphic genotype in the 1,070 samples screened. Although the assay was designed to B. napus, successful SNP amplification was achieved in the B. napus progenitor species, Brassica rapa and Brassica oleracea, and to a lesser extent in the related species Brassica nigra. Phylogenetic analysis was consistent with the expected relationships between B. napus individuals. This study presents an efficient custom SNP assay development pipeline in the complex polyploid Brassica genome and demonstrates the utility of the array for high-throughput genotyping in a number of related Brassica species. It also demonstrates the utility of this assay in genotyping resistance genes on chromosome A7, which segregate amongst the 1,070 samples.     

Identification of a nuclear-recessive gene locus for male sterility on A2 chromosome using the Brassica 60 K SNP array in non-heading Chinese cabbage

WS24-3 is a newly bred recessive genic male sterility line of the non-heading Chinese cabbage (Brassica rapa ssp. chinensis). Here, an F₂ population was produced from the cross between WS24-3 and a male-fertile breeding line (WS135). The Illumina Brassica 60 K single nucleotide polymorphism (SNP) array was used for SNPs detecting between sterile and fertile bulks from the F₂ population, and 62 SNPs were identified. BLAST analysis of the 62 SNPs revealed that the A2 chromosome of Brassica rapa genome contained 22 SNPs, whereas the other chromosomes did not contain more than 6 SNPs each. These data indicated that the potential target gene locus, named Bra2Ms, might be located on A2. Based on 10 of the 22 SNPs, allele-specific-polymerase chain reaction (AS-PCR) primers and single sequence repeat (SSR) primers were designed, 5 AS-PCR primers and 9 SSR primers showed difference between the bulks in electrophoretic determination. Analysis of these markers in F₂ population revealed that Bra2Ms was genetically delimited to a region of 1.2 cM. We also detected two co-segregated markers SSRa2-951 and SSRa2-960 in this region. The markers identified in our study might facilitate the transfer of recessive genic male sterility alleles to other favorable genetic backgrounds. Furthermore, these markers will support a map-based clone of Bra2Ms.