Identification of two loci for resistance to clubroot (<Emphasis Type="Italic">Plasmodiophora brassicae</Emphasis> Woronin) in <Emphasis Type="Italic">Brassica rapa</Emphasis> L.

In an analysis of 114 F₂ individuals from a cross between clubroot-resistant and susceptible lines of Brassica rapa L., 'G004' and 'Hakusai Chukanbohon Nou 7' (A9709), respectively, we identified two loci, Crr1 and Crr2, for clubroot (caused by Plasmodiophora brassicae Woronin) resistance. Each locus segregated independently among the F₂ population, indicating that the loci reside on a different region of chromosomes or on different chromosomes. Genetic analysis showed that each locus had little effect on clubroot resistance by itself, indicating that these two loci are complementary for clubroot resistance. The resistance to clubroot was much stronger when both loci were homozygous for resistant alleles than when they were heterozygous. These results indicate that clubroot resistance in B. rapa is under oligogenic control and at least two loci are necessary for resistance.Communicated by H.C. Becker     

Fine genetic and physical mapping of the CRb gene conferring resistance to clubroot disease in Brassica rapa

The use of clubroot resistance (CR) genes is an effective and economical approach for controlling Plasmodiophora brassicae, the causal agent of clubroot disease in Chinese cabbage (Brassica rapa) and other Brassica crops. In a previous study, we identified and mapped the CRb locus on chromosome A03 of B. rapa in the doubled-haploid (DH) line ‘CR Shinki DH line’ of Chinese cabbage. In this study, CRb, a dominant gene conferring resistance to pathotype 4 of P. brassicae, was finely mapped in combination with bulked segregant analysis and bioinformatics analysis (BIA). Using 1,486 highly susceptible individuals and 2,896 individuals from two separate F₂ populations of ‘702-5’ (B. rapa ssp. chinensis) ×  ‘CR Shinki DH line,’ the CRb locus was narrowed to a region of approximately 0.14 cM between two flanking markers, TCR79 and TCR108. The sequences of seven newly developed markers linked to CRb were landed on bacterial artificial chromosome (BAC) of the reference B. rapa ‘Chiifu-401-42’ by BIA, and a physical map consisting of three BAC clones was constructed. The CRb locus was defined as an interval of approximately 83.5 kb on a BAC clone (KBrB085J21). The target interval contained one Toll-interleukin-1 receptor/nucleotide-binding site/leucine-rich repeat (TIR–NBS–LRR) gene, one NBS–LRR gene, and several putative regulatory genes in the B. rapa genome. The CRb gene was tightly linked to two other CR genes, CRa and CRb ^Kato . These results provide useful information for isolation of the CRb gene and tightly linked molecular markers for breeding CR in B. rapa.     

Identification of Novel QTLs for Isolate-Specific Partial Resistance to Plasmodiophora brassicae in Brassica rapa

Plasmodiophora brassicae, the causal agent of clubroot disease of the Brassica crops, is widespread in the world. Quantitative trait loci (QTLs) for partial resistance to 4 different isolates of P. brassicae (Pb2, Pb4, Pb7, and Pb10) were investigated using a BC₁F₁ population from a cross between two subspecies of Brassica rapa, i.e. Chinese cabbage inbred line C59-1 as a susceptible recurrent parent and turnip inbred line ECD04 as a resistant donor parent. The BC₁F₂ families were assessed for resistance under controlled conditions. A linkage map constructed with simple sequence repeats (SSR), unigene-derived microsatellite (UGMS) markers, and specific markers linked to published clubroot resistance (CR) genes of B. rapa was used to perform QTL mapping. A total of 6 QTLs residing in 5 CR QTL regions of the B. rapa chromosomes A01, A03, and A08 were identified to account for 12.2 to 35.2% of the phenotypic variance. Two QTL regions were found to be novel except for 3 QTLs in the respective regions of previously identified Crr1, Crr2, and Crr3. QTL mapping results indicated that 1 QTL region was common for partial resistance to the 2 isolates of Pb2 and Pb7, whereas the others were specific for each isolate. Additionally, synteny analysis between B. rapa and Arabidopsis thaliana revealed that all CR QTL regions were aligned to a single conserved crucifer blocks (U, F, and R) on 3 Arabidopsis chromosomes where 2 CR QTLs were detected in A. thaliana. These results suggest that some common ancestral genomic regions were involved in the evolution of CR genes in B. rapa.     

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.     

A review of sources of resistance to turnip yellows virus (TuYV) in Brassica species

Turnip yellows virus (TuYV; previously known as beet western yellows virus) causes major diseases of Brassica species worldwide resulting in severe yield-losses in arable and vegetable crops. It has also been shown to reduce the quality of vegetables, particularly cabbage where it causes tip burn. Incidences of 100% have been recorded in commercial crops of winter oilseed rape (Brassica napus) and vegetable crops (particularly Brassica oleracea) in Europe. This review summarises the known sources of resistance to TuYV in B. napus (AACC genome), Brassica rapa (AA genome) and B. oleracea (CC genome). It also proposes names for the quantitative trait loci (QTLs) responsible for the resistances, Tu rnip Y ellows virus R esistance (TuYR), that have been mapped to at least the chromosome level in the different Brassica species. There is currently only one known source of resistance deployed commercially (TuYR1). This resistance is said to have originated in B. rapa and was introgressed into the A genome of oilseed rape via hybridisation with B. oleracea to produce allotetraploid (AACC) plants that were then backcrossed into oilseed rape. It has been utilised in the majority of known TuYV-resistant oilseed rape varieties. This has placed significant selection pressure for resistance-breaking mutations arising in TuYV. Further QTLs for resistance to TuYV (TuYR2-TuYR9) have been mapped in the genomes of B. napus, B. rapa and B. oleracea and are described here. QTLs from the latter two species have been introgressed into allotetraploid plants, providing for the first time, combined resistance from both the A and the C genomes for deployment in oilseed rape. Introgression of these new resistances into commercial oilseed rape and vegetable brassicas can be accelerated using the molecular markers that have been developed. The deployment of these resistances should lessen selection pressure for resistance-breaking isolates of TuYV and thereby prolong the effectiveness of each other and extant resistance.     

Cold pretreatment enhances microspore embryogenesis in oilseed rape (Brassica napus L.)

Stress is an essential component during embryogenesis induction in microspore culture. Cold pretreatment has been used in cereal microspore culture but very seldom attempted in Brassica microspore culture. The effect of cold pretreatment of flower buds subjected to a liquid medium on microspore embryogenesis was investigated in spring and winter Brassica napus, as well as in B. rapa and B. oleracea. Cold pretreatment significantly enhanced microspore embryogenesis (by 1–7 fold) compared to commonly used microspore culture protocol in B. napus, while it was less effective in B. rapa or even negative in B. oleracea. The appropriate duration of cold pretreatment was found to be 2–4 days, which stimulated the best microspore embryogenesis. Cold pretreatment was also able to promote embryo development including the improvement of embryo quality and acceleration of embryogenesis. When incorporating with medium refreshing, cold pretreatment could initiate the most microspore embryogenesis than any other treatment used. With further improvement cold pretreatment method may have a positive potential in Brassica breeding programmes.     

Mapping of RFLP and qualitative trait loci in Brassica rapa and comparison to the linkage maps of B. napus,B. oleracea,and Arabidopsis thaliana

A linkage map of restriction fragment length polymorphisms (RFLPs) was constructed for oilseed, Brassica rapa, using anonymous genomic DNA and cDNA clones from Brassica and cloned genes from the crucifer Arabidopsis thaliana. We also mapped genes controlling the simply inherited traits, yellow seeds, low seed erucic acid, and pubescence. The map included 139 RFLP loci organized into ten linkage groups (LGs) and one small group covering 1785 cM. Each of the three traits mapped to a single locus on three different LGs. Many of the RFLP loci were detected with the same set of probes used to construct maps in the diploid B. oleracea and the amphidiploid B. napus. Comparisons of the linkage arrangements between the diploid species B. rapa and B. oleracea revealed six LGs with at least two loci in common. Nine of the B. rapa LGs had conserved linkage arrangements with B. napus LGs. The majority of loci in common were in the same order among the three species, although the distances between loci were largest on the B. rapa map. We also compared the genome organization between B. rapa and A. thaliana using RFLP loci detected with 12 cloned genes in the two species and found some evidence for a conservation of the linkage arrangements. This B. rapa map will be used to test for associations between segregation of RFLPs, detected by cloned genes of known function, and traits of interest.     

A novel locus for clubroot resistance in<Emphasis Type="Italic"> Brassica rapa</Emphasis> and its linkage markers

An inbred turnip (Brassica rapa syn. campestris) line, N-WMR-3, which carries the trait of clubroot resistance (CR) from a European turnip, Milan White, was crossed with a clubroot-susceptible doubled haploid line, A9709. A segregating F₃ population was obtained by single-seed descent of F₂ plants and used for a genetic analysis. Segregation of CR in the F₃ population suggested that CR is controlled by a major gene. Two RAPD markers, OPC11-1 and OPC11-2, were obtained as candidates of linkage markers by bulked segregant analysis. These were converted to sequence-tagged site markers, by cloning and sequencing of the polymorphic bands, and named OPC11-1S and OPC11-2S, respectively. The specific primer pairs for OPC11-1S amplified a clear dominant band, while the primer pairs for OPC11-2S resulted in co-dominant bands. Frequency distributions and statistical analyses indicate the presence of a major dominant CR gene linked to these two markers. The present marker for CR was independent of the previously found CR loci, Crr1 andCrr2. Genotypic distribution and statistical analyses did not show any evidence of CR alleles on Crr1 andCrr2 loci in N-WMR-3. The present study clearly demonstrates that B. rapa has at least three CR loci. Therefore, the new CR locus was named Crr3. The present locus may be useful in breeding CR Chinese cabbage cultivars to overcome the decay of present CR cultivars.Communicated by C. Möllers     

Status and Perspectives of Clubroot Resistance Breeding in Crucifer Crops

Clubroot disease is a major threat to crops belonging to the Brassicaceae. It is controlled most effectively by the use of resistant cultivars. Plasmodiophora brassicae, the causal agent, shows a wide variation for pathogenicity, which can be displayed by using differential host sets. Except for Brassica juncea and B. carinata, resistant accessions can be found in all major crops. Most resistance sources are race-specific, despite some race-independent resistant accessions which can be found in B. oleracea. European field isolates from P. brassicae display great variation and show a tendency to overcome different resistance sources from either B. rapa or B. oleracea. At present, resistance genes from stubble turnips (B. rapa) are most effective and most widely used in resistance breeding of different Brassica crops. Resistance to P. brassicae from turnips was introduced into Chinese cabbage, oilseed rape, and B. oleracea. Although most turnips carry more than one resistance gene, the resistant cultivars from other crops received primarily a single, dominant resistance gene having a race-specific effect. Populations of P. brassicae that are compatible against most of the used resistance sources have been present in certain European areas for many decades. Such pathogen populations appeared in Japanese Chinese cabbage crops only a few years after the introduction of resistant cultivars. As the spread of virulent P. brassicae pathotypes seems to be slow, resistant cultivars are still a very effective method of control in many cropping areas. Mapping studies have revealed the presence of several clubroot-resistance genes in the Brassica A and C genomes; most of these genes are showing race specificity. Only in B. oleracea was one broad-spectrum locus detected. Two loci from the A genome confer resistance to more than one pathotype, but not to all isolates. Progress made in the determination of resistance loci should be used to formulate and introduce an improved differential set. Future efforts for breeding P. brassicae resistance will focus on durability by broadening the genetic basis of clubroot resistance by using either natural variation or transgenic strategies.     

Mapping loci controlling vernalization requirement in Brassica rapa

Brassica cultivars are classified as biennial or annual based on their requirement for a period of cold treatment (vernalization) to induce flowering. Genes controlling the vernalization requirement were identified in a Brassica rapa F₂ population derived from a cross between an annual and a biennial oilseed cultivar by using an RFLP linkage map and quantitative trait locus (QTL) analysis of flowering time in F₃ lines. Two genomic regions were strongly associated with variation for flowering time of unvernalized plants and alleles from the biennial parent in these regions delayed flowering. These QTLs had no significant effect on flowering time after plants were vernalized for 6 weeks, suggesting that they control flowering time through the requirement for vernalization. The two B. rapa linkage groups containing these QTLs had RFLP loci in common with two B. napus linkage groups that were shown previously to contain QTLs for flowering time. An RFLP locus detected by the cold-induced gene COR6.6 cloned from Arabidopsis thaliana mapped very near to one of the B. rapa QTLs for flowering time.     

Fine mapping of the clubroot resistance gene, Crr3, in Brassica rapa

A linkage map of Chinese cabbage (Brassica rapa) was constructed to localize the clubroot resistance (CR) gene, Crr3. Quantitative trait loci analysis using an F₃ population revealed a sharp peak in the logarithm of odds score around the sequence-tagged site (STS) marker, OPC11-2S. Therefore, this region contained Crr3. Nucleotide sequences of OPC11-2S and its proximal markers showed homology to sequences in the top arm of Arabidopsis chromosome 3, suggesting a synteny between the two species. For fine mapping of Crr3, a number of STS markers were developed based on genomic information from Arabidopsis. We obtained polymorphisms in 23 Arabidopsis-derived STS markers, 11 of which were closely linked to Crr3. The precise position of Crr3 was determined using a population of 888 F₂ plants. Eighty plants showing recombination around Crr3 locus were selected and used for the mapping. A fine map of 4.74 cM was obtained, in which two markers (BrSTS-41 and BrSTS-44) and three markers (OPC11-2S, BrSTS-54 and BrSTS-61) were cosegregated. Marker genotypes of the 21 selected F₂ families and CR tests of their progenies strongly suggested that the Crr3 gene is located in a 0.35 cM segment between the two markers, BrSTS-33 and BrSTS-78.     

In silico integration of quantitative trait loci for seed yield and yield-related traits in Brassica napus

Mapping quantitative trait loci (QTLs) is a foundation for molecular marker-assisted selection and map-based gene cloning. During the past decade, numerous QTLs for seed yield (SY) and yield-related traits in Brassica napus L. have been identified. However, integration of these results in order to compare QTLs from different mapping populations has not been undertaken, due to the lack of common molecular markers between studies. Using previously reported Brassica rapa and Brassica oleracea genome sequences, we carried out in silico integration of 1,960 QTLs associated with 13 SY and yield-related traits from 15 B. napus mapping experiments over the last decade. A total of 736 SY and yield-related QTLs were mapped onto 283 loci in the A and C genomes of B. napus. These QTLs were unevenly distributed across the 19 B. napus chromosomes, with the most on chromosome A3 and the least on chromosome C6. Our integrated QTL map identified 142 loci where the conserved QTLs were detected and 25 multifunctional loci, mostly for the traits of flowering time (FT), plant height, 1,000-seed weight, maturity time and SY. These conserved QTLs and multifunctional loci may result from pleiotropism or clustered genes. At the same time, a total of 146 genes underlying the QTLs for FT and other yield-related traits were identified by comparative mapping with the Arabidopsis genome. These results facilitate the retrieval of B. napus SY and yield-related QTLs for research communities, increase the density of targeted QTL-linked markers, validate the existence of QTLs across different populations, and advance the fine mapping of genes.     

Detection,introgression and localization of genes conferring specific resistance to <Emphasis Type="Italic">Leptosphaeria maculans</Emphasis> from <Emphasis Type="Italic">Brassica rapa</Emphasis> into <Emphasis Type="Italic">B. napus</Emphasis>

Blackleg (stem canker) caused by the fungus Leptosphaeria maculans is one of the most damaging diseases of oilseed rape (Brassica napus). Crop relatives represent a valuable source of “new” resistance genes that could be used to diversify cultivar resistance. B. rapa, one of the progenitors of B. napus, is a potential source of new resistance genes. However, most of the accessions are heterozygous so it is impossible to directly detect the plant genes conferring specific resistance due to the complex patterns of avirulence genes in L. maculans isolates. We developed a strategy to simultaneously characterize and introgress resistance genes from B. rapa, by homologous recombination, into B. napus. One B. rapa plant resistant to one L. maculans isolate was used to produce B. rapa backcross progeny and a resynthesized B. napus plant from which a population of doubled haploid lines was derived after crossing with natural B. napus. We then used molecular analyses and resistance tests on these populations to identify and map the resistance genes and to characterize their introgression from B. rapa into B. napus. Three specific genes conferring resistance to L. maculans (Rlm1, Rlm2 and Rlm7) were identified in B. rapa. Comparisons of genetic maps showed that two of these genes were located on the R7 linkage group, in a region homologous to the region on linkage group N7 in B. napus, where these genes have been reported previously. The results of our study offer new perspectives for gene introgression and cloning in Brassicas.     

Distribution of similar self-incompatibility (<Emphasis Type="Italic">S</Emphasis>) haplotypes in different genera,<Emphasis Type="Italic"> Raphanus</Emphasis> and<Emphasis Type="Italic"> Brassica</Emphasis>

The nucleotide sequences of ten SP11 and nine SRK alleles in Raphanus sativus were determined, and deduced amino acid sequences were compared with those of Brassica SP11 and SRK. The amino acid sequence identity of class-I SP11s in R. sativus was about 30% on average, the highest being 52.2%, while that of the S domain of class-I SRK was 77.0% on average and ranged from 70.8% to 83.9%. These values were comparable to those of SP11 and SRK in Brassica oleracea and B. rapa. SP11 of R. sativus S-21 was found to be highly similar to SP11 of B. rapa S-9 (89.5% amino acid identity), and SRK of R. sativus S-21 was similar to SRK of B. rapa S-9 (91.0%). SP11 and SRK of R. sativus S-19 were also similar to SP11 and SRK of B. oleracea S-20, respectively. These similarities of both SP11 and SRK alleles between R. sativus and Brassica suggest that these S haplotype pairs originated from the same ancestral S haplotypes.     

Sequence‐level comparative analysis of the Brassica napus genome around two stearoyl‐ACP desaturase loci

We conducted a sequence‐level comparative analyses, at the scale of complete bacterial artificial chromosome (BAC) clones, between the genome of the most economically important Brassica species, Brassica napus (oilseed rape), and those of Brassica rapa, the genome of which is currently being sequenced, and Arabidopsis thaliana. We constructed a new B. napus BAC library and identified and sequenced clones that contain homoeologous regions of the genome including stearoyl‐ACP desaturase‐encoding genes. We sequenced the orthologous region of the genome of B. rapa and conducted comparative analyses between the Brassica sequences and those of the orthologous region of the genome of A. thaliana. The proportion of genes conserved (～56%) is lower than has been reported previously between A. thaliana and Brassica (～66%). The gene models for sets of conserved genes were used to determine the extent of nucleotide conservation of coding regions. This was found to be 84.2 ± 3.9% and 85.8 ± 3.7% between the B. napus A and C genomes, respectively, and that of A. thaliana, which is consistent with previous results for other Brassica species, and 97.5 ± 3.1% between the B. napus A genome and B. rapa, and 93.1 ± 4.9% between the B. napus C genome and B. rapa. The divergence of the B. napus genes from the A genome and the B. rapa genes was greater than anticipated and indicates that the A genome ancestor of the B. napus cultivar studied was relatively distantly related to the cultivar of B. rapa selected for genome sequencing.     

Comparative mapping of loci controlling winter survival and related traits in oilseed Brassica rapa and B. napus

Winter survival is an important characteristic of oilseedBrassica that is seeded in the fall in northern climates,and it may be affected by genetic variation for other cold-regulated traits,such as freezing tolerance and vernalization responsive flowering time. Weanalyzed immortalized populations of oilseed Brassica rapa(recombinant inbred lines) and B. napus (double haploidlines) derived from crosses of annual and biennial types in order to comparethe map positions and effects of quantitative trait loci controlling wintersurvival, nonacclimated and acclimated freezing tolerances, and flowering time.The B. napus population was evaluated in multiple winters,and six of the 16 total significant QTL for winter survival were detected inmore than one winter. Correspondence in the map positions of QTL controllingdifferent traits within species provided evidence that some alleles causinggreater acclimated freezing tolerance and later flowering time also contributedto increased winter survival. Correspondence in the map positions of QTLbetween species provided evidence for allelic variation at homologous loci inB. rapa and B. napus. The potentialrole of some candidate genes in regulating these traits is discussed.     

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.     

Sequence and expression of endogenous S-locus glycoprotein genes in self-compatible Brassica napus

Brassica napus is an amphidiploid plant which is self-compatible even though it is derived from hybridisation of the self-incompatible species B. oleracea and B. campestris. Experiments were undertaken to establish if S-locus glycoprotein (SLG) genes exist in B. napus and whether these are expressed as in self-incompatible Brassica species. Two different stigma-specific cDNA sequences homologous to SLG genes were obtained from the B. napus cultivar Westar. One of these sequences, SLG _WS1, displayed highest homology to class I SLG alleles, whereas the other, SLG _WS2, showed greatest homology to class II SLG genes. Both were expressed at high levels in Westar stigmas following a developmental pattern typical of SLG genes in the self-incompatible diploids. We infer that they represent the endogenous SLG genes at the two homoeologous S-loci. The occurrence of normally expressed SLG genes and its relevance to the self-compatible phenotype of B. napus is discussed.     

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.