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.     

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.     

Identifying the chromosomes of the A- and C-genome diploid Brassica species B. rapa (syn. campestris) and B. oleracea in their amphidiploid B. napus

Oilseed rape (Brassica napus L.) is an amphidiploid species that originated from a spontaneous hybridisation of Brassica rapa L. (syn. campestris) and Brassica oleracea L., and contains the complete diploid chromosome sets of both parental genomes. The metaphase chromosomes of the highly homoeologous A genome of B. rapa and the C genome of B. oleracea cannot be reliably distinguished in B. napus because of their morphological similarity. Fluorescence in situ hybridisation (FISH) with 5S and 25S ribosomal DNA probes to prometaphase chromosomes, in combination with DAPI staining, allows more dependable identification of Brassica chromosomes. By comparing rDNA hybridisation and DAPI staining patterns from B. rapa and B. oleracea prometaphase chromosomes with those from B. napus, we were able to identify the putative homologues of B. napus chromosomes in the diploid chromosome sets of B. rapa and B. oleracea, respectively. In some cases, differences were observed between the rDNA hybridisation patterns of chromosomes in the diploid species and their putative homologue in B. napus, indicating locus losses or alterations in rDNA copy number. The ability to reliably identify A and C genome chromosomes in B. napus is discussed with respect to evolutionary and breeding aspects. Received: 13 July 2001 / Accepted: 23 August 2001     

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.     

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 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.     

Molecular and cytological analyses of A and C genomes at meiosis in synthetic allotriploid <Emphasis Type="Italic">Brassica</Emphasis> hybrids (ACC) between <Emphasis Type="Italic">B.napus</Emphasis> (AACC) and <Emphasis Type="Italic">B.oleracea</Emphasis> (CC)

Success of interspecific hybridization relies mostly on the adequate similarity between the implicated genomes to ensure synapsis, pairing and recombination between appropriate chromosomes during meiosis in allopolyploid species. Allotetraploid Brassica napus (AACC) is a model of natural hybridization between Brassica rapa (AA) and Brassica oleracea (CC), which are originally derived from a common ancestor, but genomic constitution of the same chromosomes probably varied among these species through time after establishment, giving rise to cytogenetic difference in the synthetic hybrids. Herein we investigated meiotic behaviors of A and C chromosomes of synthetic allotriploid Brassica hybrids (ACC) at molecular and cytological levels, which result from the interspecific cross between natural B. napus (AACC) and B.oleracea (CC), and the results showed that meiosis course was significantly aberrant in allotriploid Brassica hybrids, and chromosomes aligned chaotically at metaphase I, chromosome bridges and lags were frequently observed from later metaphase I to anaphase II during meiosis. Simultaneously, we also noticed that meiosis-related genes were abruptly down-regulated in allotriploid Brassica hybrids, which likely accounted for irregular scenario of meiosis observed in these synthetic hybrids. Therefore, these results indicated that inter-genomic exchanges of A and C chromosomes could occur frequently in synthetic Brassica hybrids, and provided an efficient approach for genetic changes of homeologous chromosomes during meiosis in polyploid B.napus breeding program.     

Pairing and recombination at meiosis of Brassica rapa (AA) × Brassica napus (AACC) hybrids

Interspecific crosses contribute significantly to plant evolution enabling gene exchanges between species. The efficiency of interspecific crosses depends on the similarity between the implicated genomes as high levels of genome similarity are required to ensure appropriate chromosome pairing and genetic recombination. Brassica napus (AACC) is an allopolyploid, resulting from natural hybridization between Brassica rapa (AA) and Brassica oleracea (CC), both being diploid species derived from a common ancestor. To study the relationships between genomes of these Brassica species, we have determined simultaneously the pairing and recombination pattern of A and C chromosomes during meiosis of AAC triploid hybrids, which result from the interspecific cross between natural B. napus and B. rapa. Different AAC triploid hybrids and their progenies have been analysed using cytogenetic, BAC-FISH, and molecular techniques. In 71% of the pollen mother cells, homologous A chromosomes paired regularly, and usually one chromosome of each pair was transmitted to the progeny. C chromosomes remained mainly univalent, but were involved in homoeologous pairing in 21.5% of the cells, and 13% of the transmitted C chromosomes were either recombined or broken. The rate of transmission of C chromosomes depended on the identity of the particular chromosome and on the way the hybrid was crossed, as the male or as the female parent, to B. napus or to B. rapa. Gene transfers in triploid hybrids are favoured between A genomes of B. rapa and B. napus, but also occur between A and C genomes though at lower rates.     

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.     

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.     

Extensive tRNA Gene Changes in Synthetic Brassica napus

Allopolyploidization, where two species come together to form a new species, plays a major role in speciation and genome evolution. Transfer RNAs (abbreviated tRNA) are typically 73–94 nucleotides in length, and are indispensable in protein synthesis, transferring amino acids to the cell protein synthesis machinery (ribosome). To date, the regularity and function of tRNA gene sequence variation during the process of allopolyploidization have not been well understood. In this study, the inter-tRNA gene corresponding to tRNA amplification polymorphism method was used to detect changes in tRNA gene sequences in the progeny of interspecific hybrids between Brassica rapa and B. oleracea, mimicking the original B. napus (canola) species formation event. Cluster analysis showed that tRNA gene variation during allopolyploidization did not appear to have a genotypic basis. Significant variation occurred in the early generations of synthetic B. napus (F₁ and F₂ generations), but fewer alterations were observed in the later generation (F₃). The variation-prone tRNA genes tended to be located in AT-rich regions. BlastN analysis of novel tRNA gene variants against a Brassica genome sequence database showed that the variation of these tRNA-gene-associated sequences in allopolyploidization might result in variation of gene structure and function, e.g., metabolic process and transport.     

Molecular characterization of the flowering time gene FRIGIDA in Brassica genomes A and C

An important determinant of flowering time variation in Arabidopsis, the FRIGIDA (FRI) gene has not been until recently investigated in economically important Brassica species. In diploid Brassica species, this gene exists as two paralogous loci on chromosomes A3 and A4 (B. rapa; A genome), and C3 and C9 (B. oleracea; C genome). Each locus is represented by several genome-specific alleles, which are discerned primarily by polymorphisms in C- and especially N-terminal regions. Locus- and genome-specific sequences of two FRI paralogues are conserved almost completely in the subgenomes A and C of tetraploid B. napus. The phylogenetic analysis of available FRI sequences presumes that the duplication of FRI loci preceded speciation in the genus Brassica.     

Inter- and intra-genomic homology of the Brassica genomes: implications for their origin and evolution

In order to determine the homologous regions shared by the cultivated Brassica genomes, linkage maps of the diploid cultivated B. rapa (A genome, n = 10), B. nigra (B genome, n = 8) and B. oleracea (C genome, n = 9), were compared. We found intergenomic conserved regions but with extensitve reordering among the genomes. Eighteen linkage groups from all three species could be associated on the basis of homologous segments based on at least three common markers. Intragenomic homologous conservation was also observed for some of the chromosomes of the A, B and C genomes. A possible chromosome phylogenetic pathway based on an ancestral genome of at least five, and no more than seven chromosomes, was drawn from the chromosomal inter-relationships observed. These results demonstrate that extensive duplication and rearrangement have been involved in the formation of the Brassica genomes from a smaller ancestral genome.     

Deciphering the Diploid Ancestral Genome of the Mesohexaploid Brassica rapa

The genus Brassica includes several important agricultural and horticultural crops. Their current genome structures were shaped by whole-genome triplication followed by extensive diploidization. The availability of several crucifer genome sequences, especially that of Chinese cabbage (Brassica rapa), enables study of the evolution of the mesohexaploid Brassica genomes from their diploid progenitors. We reconstructed three ancestral subgenomes of B. rapa (n = 10) by comparing its whole-genome sequence to ancestral and extant Brassicaceae genomes. All three B. rapa paleogenomes apparently consisted of seven chromosomes, similar to the ancestral translocation Proto-Calepineae Karyotype (tPCK; n = 7), which is the evolutionarily younger variant of the Proto-Calepineae Karyotype (n = 7). Based on comparative analysis of genome sequences or linkage maps of Brassica oleracea, Brassica nigra, radish (Raphanus sativus), and other closely related species, we propose a two-step merging of three tPCK-like genomes to form the hexaploid ancestor of the tribe Brassiceae with 42 chromosomes. Subsequent diversification of the Brassiceae was marked by extensive genome reshuffling and chromosome number reduction mediated by translocation events and followed by loss and/or inactivation of centromeres. Furthermore, via interspecies genome comparison, we refined intervals for seven of the genomic blocks of the Ancestral Crucifer Karyotype (n = 8), thus revising the key reference genome for evolutionary genomics of crucifers.     

Comparative sequence analysis for Brassica oleracea with similar sequences in B. rapa and Arabidopsis thaliana

We sequenced five BAC clones of Brassica oleracea doubled haploid ‘Early Big' broccoli containing major genes in the aliphatic glucosinolate pathway, and comparatively analyzed them with similar sequences in A. thaliana and B. rapa. Additionally, we included in the analysis published sequences from three other B. oleracea BAC clones and a contig of this species corresponding to segments in A. thaliana chromosomes IV and V. A total of 2,946 kb of B. oleracea, 1,069 kb of B. rapa sequence and 2,607 kb of A. thaliana sequence were compared and analyzed. We found conserved collinearity for gene order and content restricted to specific chromosomal segments, but breaks in collinearity were frequent resulting in gene absence likely not due to gene loss but rearrangements. B. oleracea has the lowest gene density of the three species, followed by B. rapa. The genome expansion of the Brassica species, B. oleracea in particular, is due to larger introns and gene spacers resulting from frequent insertion of DNA transposons and retrotransposons. These findings are discussed in relation to the possible origin and evolution of the Brassica genomes. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.     

Cross-transferability and Polymorphic Potential of Genomic STMS Markers of Brassica Species

The present study was carried out with the objective of evaluating genomic STMS markers developed earlier in Brassica napus, B. oleracea, B. rapa and B. nigra for their use in Brassica juncea and B. carinata. Ninety-six of the 100 STMS markers used under standardized annealing temperatures and gel concentrations produced clear reproducible amplification pattern. For majority of the markers 60 °C annealing temperature and 3.5% metaphor agarose gel were found suitable. High cross-transferability of STMS markers to related Brassica species including B. carinata (91.6%) and B. juncea (87.5%) suggested the possibility of utilizing these markers for genome analysis in the species where no such markers are available. The ‘B’ genome derived markers showed lower level of transferability to the ‘A’ and ‘C’ genome Brassica species. The potential of STMS markers to detect polymorphism among Brassica species and genera was 98.9%. The level of inter-specific polymorphism was much higher than the intea-specific polymorphism. The markers capable of revealing polymorphism among Brassica species and genera would be useful in Brassica introgression breeding programme. The polymorphic markers were found efficient in establishing the expected evolutionary relationships among the six different Brassica species and two related genera. Low level of intra-specific polymorphism revealed by these markers suggested use of a large set of such markers for various applications in Brassica genetics, genomics and breeding.