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 Analysis between Homoeologous Genome Segments of Brassica napus and Its Progenitor Species Reveals Extensive Sequence-Level Divergence

Homoeologous regions of Brassica genomes were analyzed at the sequence level. These represent segments of the Brassica A genome as found in Brassica rapa and Brassica napus and the corresponding segments of the Brassica C genome as found in Brassica oleracea and B. napus. Analysis of synonymous base substitution rates within modeled genes revealed a relatively broad range of times (0.12 to 1.37 million years ago) since the divergence of orthologous genome segments as represented in B. napus and the diploid species. Similar, and consistent, ranges were also identified for single nucleotide polymorphism and insertion-deletion variation. Genes conserved across the Brassica genomes and the homoeologous segments of the genome of Arabidopsis thaliana showed almost perfect collinearity. Numerous examples of apparent transduplication of gene fragments, as previously reported in B. oleracea, were observed in B. rapa and B. napus, indicating that this phenomenon is widespread in Brassica species. In the majority of the regions studied, the C genome segments were expanded in size relative to their A genome counterparts. The considerable variation that we observed, even between the different versions of the same Brassica genome, for gene fragments and annotated putative genes suggest that the concept of the pan-genome might be particularly appropriate when considering Brassica genomes.     

Sequence and structure of Brassica rapa chromosome A3

Background

The species Brassica rapa includes important vegetable and oil crops. It also serves as an excellent model system to study polyploidy-related genome evolution because of its paleohexaploid ancestry and its close evolutionary relationships with Arabidopsis thaliana and other Brassica species with larger genomes. Therefore, its genome sequence will be used to accelerate both basic research on genome evolution and applied research across the cultivated Brassica species.

Results

We have determined and analyzed the sequence of B. rapa chromosome A3. We obtained 31.9 Mb of sequences, organized into nine contigs, which incorporated 348 overlapping BAC clones. Annotation revealed 7,058 protein-coding genes, with an average gene density of 4.6 kb per gene. Analysis of chromosome collinearity with the A. thaliana genome identified conserved synteny blocks encompassing the whole of the B. rapa chromosome A3 and sections of four A. thaliana chromosomes. The frequency of tandem duplication of genes differed between the conserved genome segments in B. rapa and A. thaliana, indicating differential rates of occurrence/retention of such duplicate copies of genes. Analysis of 'ancestral karyotype' genome building blocks enabled the development of a hypothetical model for the derivation of the B. rapa chromosome A3.

Conclusions

We report the near-complete chromosome sequence from a dicotyledonous crop species. This provides an example of the complexity of genome evolution following polyploidy. The high degree of contiguity afforded by the clone-by-clone approach provides a benchmark for the performance of whole genome shotgun approaches presently being applied in B. rapa and other species with complex genomes.     

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 set of public SSR markers derived from genomic sequence of a rapid cycling <Emphasis Type="Italic">Brassica oleracea</Emphasis> L. genotype

The traditional development of simple sequence repeat (SSR) or microsatellite markers by probe hybridization can be time-consuming and requires the use of specialized laboratory equipment. In this study, probe hybridization was circumvented by using sequence information on 3,500 genomic clones mainly from Brassica oleracea to identify di, tri, tetra and penta-nucleotide repeats. A total of 587 primer pairs flanking SSR were developed using this approach. From these, 420 SSR markers amplified DNA in two parental lines of B. rapa (26% were polymorphic) and 523 in two parental lines of B. oleracea (32% were polymorphic). A diverse array of motif types was identified, characterized and compared with traditional SSR detection methods. The most abundant motifs found were di- (38%) and trinucleotides (33%) followed by penta- (16%) and tetranucleotide (13%) motifs. The type of motif class, motif length and repeat were not indicative of polymorphisms. The frequency of B. oleracea SSRs in genomic shotgun sequence was estimated to be 1 every 4 Kb. In general, the average motif length and repeat numbers were shorter than those obtained previously by probe hybridization, and they contained a more balanced representation of SSR motif types in the genome by identifying those that do not hybridize well to DNA probes. Brassica genomic DNA sequence information is a promising resource for developing a large number of SSR molecular markers in Brassica species. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.     

Characterization of terminal-repeat retrotransposon in miniature (TRIM) in Brassica relatives

We have newly identified five Terminal-repeat retrotransposon in miniature (TRIM) families, four from Brassica and one from Arabidopsis. A total of 146 elements, including three Arabidopsis families reported before, are extracted from genomics data of Brassica and Arabidopsis, and these are grouped into eight distinct lineages, Br1 to Br4 derived from Brassica and At1 to At4 derived from Arabidopsis. Based on the occurrence of TRIM elements in 434 Mb of B. oleracea shotgun sequences and 96 Mb of B. rapa BAC end sequences, total number of TRIM members of Br1, Br2, Br3, and Br4 families are roughly estimated to be present in 660 and 530 copies in B. oleracea and B. rapa genomes, respectively. Studies on insertion site polymorphisms of four elements across taxa in the tribe Brassiceae infer the taxonomic lineage and dating of the insertion time. Active roles of the TRIM elements for evolution of the duplicated genes are inferred in the highly replicated Brassica genome. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users. Tae-Jin Yang and Soo-Jin Kwon have equally contributed to this work.     

Development of public immortal mapping populations,molecular markers and linkage maps for rapid cycling <Emphasis Type="Italic">Brassica rapa</Emphasis> and <Emphasis Type="Italic">B. oleracea</Emphasis>

Publicly available genomic tools help researchers integrate information and make new discoveries. In this paper, we describe the development of immortal mapping populations of rapid cycling, self-compatible lines, molecular markers, and linkage maps for Brassica rapa and B. oleracea and make the data and germplasm available to the Brassica research community. The B. rapa population consists of 160 recombinant inbred (RI) lines derived from the cross of highly inbred lines of rapid cycling and yellow sarson B. rapa. The B. oleracea population consists of 155 double haploid (DH) lines derived from an F1 cross between two DH lines, rapid cycling and broccoli. A total of 120 RFLP probes, 146 SSR markers, and one phenotypic trait (flower color) were used to construct genetic linkage maps for both species. The B. rapa map consists of 224 molecular markers distributed along 10 linkage groups (A1–A10) with a total distance of 1125.3 cM and a marker density of 5.7 cM/marker. The B. oleracea genetic map consists of 279 molecular markers and one phenotypic marker distributed along nine linkage groups (C1–C9) with a total distance of 891.4 cM and a marker density of 3.2 cM/marker. A syntenic analysis with Arabidopsis thaliana identified collinear genomic blocks that are in agreement with previous studies, reinforcing the idea of conserved chromosomal regions across the Brassicaceae.     

Physical mapping and microsynteny of Brassica rapa ssp. pekinensis genome corresponding to a 222 kbp gene-rich region of Arabidopsis chromosome 4 and partially duplicated on chromosome 5

We constructed a bacterial artificial chromosome (BAC) library, designated as KBrH, from high molecular weight genomic DNA of Brassica rapa ssp. pekinensis (Chinese cabbage). This library, which was constructed using HindIII-cleaved genomic DNA, consists of 56,592 clones with average insert size of 115 kbp. Using a partially duplicated DNA sequence of Arabidopsis, represented by 19 and 9 predicted genes on chromosome 4 and 5, respectively, and BAC clones from the KBrH library, we studied conservation and microsynteny corresponding to the Arabidopsis regions in B. rapa ssp. pekinensis. The BAC contigs assembled according to the Arabidopsis homoeologues revealed triplication and rearrangements in the Chinese cabbage. In general, collinearity of genes in the paralogous segments was maintained, but gene contents were highly variable with interstitial losses. We also used representative BAC clones, from the assembled contigs, as probes and hybridized them on mitotic (metaphase) and/or meiotic (leptotene/pachytene/metaphase I) chromosomes of Chinese cabbage using bicolor fluorescence in situ hybridization. The hybridization pattern physically identified the paralogous segments of the Arabidopsis homoeologues on B. rapa ssp. pekinensis chromosomes. The homoeologous segments corresponding to chromosome 4 of Arabidopsis were located on chromosomes 2, 8 and 7, whereas those of chromosome 5 were present on chromosomes 6, 1 and 4 of B. rapa ssp. pekinensis.     

A rich TILLING resource for studying gene function in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Background  

The Brassicaceae family includes the model plant Arabidopsis thaliana as well as a number of agronomically important species such as oilseed crops (in particular Brassica napus, B. juncea and B. rapa) and vegetables (eg. B. rapa and B. oleracea).     

An ultradense genetic recombination map for <Emphasis Type="Italic">Brassica napus</Emphasis>, consisting of 13551 SRAP markers

Sequence related amplified polymorphism (SRAP) was used to construct an ultradense genetic recombination map for a doubled haploid (DH) population in B. napus. A total of 1,634 primer combinations including 12 fluorescently labeled primers and 442 unlabeled ones produced 13,551 mapped SRAP markers. All these SRAPs were assembled in 1,055 bins that were placed onto 19 linkage groups. Ten of the nineteen linkage groups were assigned to the A genome and the remaining nine to the C genome on the basis of the differential SRAP PCR amplification in two DH lines of B. rapa and B. oleracea. Furthermore, all 19 linkage groups were assigned to their corresponding N1–N19 groups of B. napus by comparison with 55 SSR markers used to construct previous maps in this species. In total, 1,663 crossovers were detected, resulting in a map length span of 1604.8 cM. The marker density is 8.45 SRAPs per cM, and there could be more than one marker in 100 kb physical distance. There are four linkage groups in the A genome with more than 800 SRAP markers each, and three linkage groups in the C genome with more 1,000 SRAP markers each. Our studies suggest that a single SRAP map might be applicable to the three Brassica species, B. napus, B. oleracea and B. rapa. The use of this ultra high-density genetic recombination map in marker development and map-based gene cloning is discussed. 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.     

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.     

Brassica taxonomy based on nuclear restriction fragment length polymorphisms (RFLPs)

Summary RFLPs were used to study genome evolution and phylogeny in Brassica and related genera. Thirtyeight accessions, including 10 accessions of B. rapa (syn. campestris), 9 cultivated types of B. oleracea, 13 nine-chromosome wild brassicas related to B. oleracea, and 6 other species in Brassica and allied genera, were examined with more then 30 random genomic DNA probes, which identified RFLPs mapping to nine different linkage groups of the B. rapa genome. Based on the RFLP data, phylogenetic trees were constructed using the PAUP microcomputer program. Within B. rapa, accessions of pak choi, narinosa, and Chinese cabbage from East Asia constituted a group distinct from turnip and wild European populations, consistent with the hypothesis that B. rapa had two centers of domestication. A wild B. rapa accession from India was positioned in the tree between European types and East Asian types, suggesting an evolutionary pathway from Europe to India, then to South China. Cultivated B. oleracea morphotypes showed monophyletic origin with wild B. oleracea or B. alboglabra as possible ancestors. Various kales constitute a highly diverse group, and represent the primitive morphotypes of cultivated B. oleracea from which cabbage, broccoli, cauliflower, etc. probably have evolved. Cauliflower was found to be closely related to broccoli, whereas cabbage was closely related to leafy kales. A great diversity existed among the 13 collections of nine-chromosome wild brassicas related to B. oleracea, representing various taxonomic states from subspecies to species. Results from these studies suggested that two basic evolutionary pathways exist for the diploid species examined. One pathway gave rise to B. fruticulosa, B. nigra, and Sinapis arvensis, with B. adpressa or a close relative as the initial ancestor. Another pathway gave rise to B. oleracea and B. rapa, with Diplotaxis erucoides or a close relative as the initial ancestor. Raphanus sativus and Eruca sativus represented intermediate types between the two lineages, and might have been derived from introgression or hybridization between species belonging to different lineages. Molecular evidence for an ascending order of chromosome numbers in the evolution of Brassica and allied genera was obtained on the basis of RFLP data and phylogenetic analysis.     

Transferability and genome specificity of a new set of microsatellite primers among Brassica species of the U triangle

We present a new set of 12 highly polymorphic simple sequence repeat primer sequences for use with Brassica species. These new primers, and four from A.K.S. SzewcMcFadden and colleagues, were tested in four Brassica species (B. rapa, B. napus, B. oleracea and B. nigra). Most primers successfully amplified products within all species and were polymorphic. Due to the risk of gene flow from GM oilseed rape to its wild relatives, hybrid formation in the Brassicaceae is of great interest. We identify six primer pairs as specific to the A, B or C genomes that could be used to identify such hybrids.     

Quantitative karyotyping of three diploid Brassica species by imaging methods and localization of 45s rDNA loci on the identified chromosomes

 Chromosomes of the three diploid Brassica species, B. rapa (AA), B. nigra (BB) and B. oleracea (CC), were identified based on their morphological characteristics, especially on the condensation pattern appearing at the somatic pro-metaphase stage. The morphological features of the pro-metaphase chromosomes of the three Brassica spp. were quantified by imaging methods using chromosome image analyzing system II (CHIAS 2). As a result, quantitative chromosome maps or idiograms of the three diploid Brassica spp. were developed. The fluorescence in situ hybridization (FISH) method revealed the location of 45s rDNA (the 26s-5.8s-18s ribosomal RNA gene cluster) on the chromosomes involved. The number of 45s rDNA loci in the B. rapa, B. nigra and B. oleracea are five, three and two, respectively. The loci detected were systematically mapped on the idiograms of the three Brassica spp. Received: 5 September 1997 / Accepted: 6 October 1997     

High‐throughput multiplex cpDNA resequencing clarifies the genetic diversity and genetic relationships among Brassica napus,Brassica rapa and Brassica oleracea

Brassica napus (rapeseed) is a recent allotetraploid plant and the second most important oilseed crop worldwide. The origin of B. napus and the genetic relationships with its diploid ancestor species remain largely unresolved. Here, chloroplast DNA (cpDNA) from 488 B. napus accessions of global origin, 139 B. rapa accessions and 49 B. oleracea accessions were populationally resequenced using Illumina Solexa sequencing technologies. The intraspecific cpDNA variants and their allelic frequencies were called genomewide and further validated via EcoTILLING analyses of the rpo region. The cpDNA of the current global B. napus population comprises more than 400 variants (SNPs and short InDels) and maintains one predominant haplotype (Bncp1). Whole‐genome resequencing of the cpDNA of Bncp1 haplotype eliminated its direct inheritance from any accession of the B. rapa or B. oleracea species. The distribution of the polymorphism information content (PIC) values for each variant demonstrated that B. napus has much lower cpDNA diversity than B. rapa; however, a vast majority of the wild and cultivated B. oleracea specimens appeared to share one same distinct cpDNA haplotype, in contrast to its wild C‐genome relatives. This finding suggests that the cpDNA of the three Brassica species is well differentiated. The predominant B. napus cpDNA haplotype may have originated from uninvestigated relatives or from interactions between cpDNA mutations and natural/artificial selection during speciation and evolution. These exhaustive data on variation in cpDNA would provide fundamental data for research on cpDNA and chloroplasts.     

Genome-wide identification of NBS-encoding resistance genes in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Nucleotide-binding site (NBS)-encoding resistance genes are key plant disease-resistance genes and are abundant in plant genomes, comprising up to 2% of all genes. The availability of genome sequences from several plant models enables the identification and cloning of NBS-encoding genes from closely related species based on a comparative genomics approach. In this study, we used the genome sequence of Brassica rapa to identify NBS-encoding genes in the Brassica genome. We identified 92 non-redundant NBS-encoding genes [30 CC-NBS-LRR (CNL) and 62 TIR-NBS-LRR (TNL) genes] in approximately 100 Mbp of B. rapa euchromatic genome sequence. Despite the fact that B. rapa has a significantly larger genome than Arabidopsis thaliana due to a recent whole genome triplication event after speciation, B. rapa contains relatively small number of NBS-encoding genes compared to A. thaliana, presumably because of deletion of redundant genes related to genome diploidization. Phylogenetic and evolutionary analyses suggest that relatively higher relaxation of selective constraints on the TNL group after the old duplication event resulted in greater accumulation of TNLs than CNLs in both Arabidopsis and Brassica genomes. Recent tandem duplication and ectopic deletion are likely to have played a role in the generation of novel Brassica lineage-specific resistance genes.