Conservation of the microstructure of genome segments in Brassica napus and its diploid relatives

The cultivated Brassica species are the group of crops most closely related to Arabidopsis thaliana (Arabidopsis). They represent models for the application in crops of genomic information gained in Arabidopsis and provide an opportunity for the investigation of polyploid genome formation and evolution. The scientific literature contains contradictory evidence for the dynamics of the evolution of polyploid genomes. We aimed at overcoming the inherent complexity of Brassica genomes and clarify the effects of polyploidy on the evolution of genome microstructure in specific segments of the genome. To do this, we have constructed bacterial artificial chromosome (BAC) libraries from genomic DNA of B. rapa subspecies trilocularis (JBr) and B. napus var Tapidor (JBnB) to supplement an existing BAC library from B. oleracea. These allowed us to analyse both recent polyploidization (under 10,000 years in B. napus) and more ancient polyploidization events (ca. 20 Myr for B. rapa and B. oleracea relative to Arabidopsis), with an analysis of the events occurring on an intermediate time scale (over the ca. 4 Myr since the divergence of the B. rapa and B. oleracea lineages). Using the Arabidopsis genome sequence and clones from the JBr library, we have analysed aspects of gene conservation and microsynteny between six regions of the genome of B. rapa with the homoeologous regions of the genomes of B. oleracea and Arabidopsis. Extensive divergence of gene content was observed between the B. rapa paralogous segments and their homoeologous segments within the genome of Arabidopsis. A pattern of interspersed gene loss was identified that is similar, but not identical, to that observed in B. oleracea. The conserved genes show highly conserved collinearity with their orthologues across genomes, but a small number of species-specific rearrangements were identified. Thus the evolution of genome microstructure is an ongoing process. Brassica napus is a recently formed polyploid resulting from the hybridization of B. rapa (containing the Brassica A genome) and B. oleracea (containing the Brassica C genome). Using clones from the JBnB library, we have analysed the microstructure of the corresponding segments of the B. napus genome. The results show that there has been little or no change to the microstructure of the analysed segments of the Brassica A and C genomes as a consequence of the hybridization event forming natural B. napus. The observations indicate that, upon polyploid formation, these segments of the genome did not undergo a burst of evolution discernible at the scale of microstructure.     

Comparative physical mapping of segments of the genome of Brassica oleracea var. alboglabra that are homoeologous to sequenced regions of chromosomes 4 and 5 of Arabidopsis thaliana

Due to their relatedness to Arabidopsis thaliana (Arabidopsis), the cultivated Brassica species represent the first group of crops with which to evaluate comparative genomics approaches to understanding biological processes and manipulating traits. We have constructed a high-quality binary BAC library (JBo) from genomic DNA of Brassica oleracea var. alboglabra, in order to underpin such investigations. Using the Arabidopsis genome sequence and clones from the JBo library, we have analysed aspects of gene conservation and microsynteny between a 222 kb region of the genome of Arabidopsis and homoeologous segments of the genome of B. oleracea. All 19 predicted genes tested were found to hybridize to clones in the JBo library, indicating a high level of gene conservation. Further analyses and physical mapping with the BAC clones identified allowed us to construct clone contig maps and analyse in detail the gene content and organization in the set of paralogous segments identified in the genome of B. oleracea. Extensive divergence of gene content was observed, both between the B. oleracea paralogous segments and between them and their homoeologous segment within the genome of Arabidopsis. However, the genes present show highly conserved collinearity with their orthologues in the genome of Arabidopsis. We have identified one example of a Brassica gene in a non-collinear position and one rearrangement. Some of the genes not present in the discernible homoeologous regions appear to be located elsewhere in the B. oleracea genome. The implications of our findings for comparative map-based cloning of genes from crop species are discussed.     

Altered patterns of fractionation and exon deletions in Brassica rapa support a two-step model of paleohexaploidy

The genome sequence of the paleohexaploid Brassica rapa shows that fractionation is biased among the three subgenomes and that the least fractionated subgenome has approximately twice as many orthologs as its close (and relatively unduplicated) relative Arabidopsis than had either of the other two subgenomes. One evolutionary scenario is that the two subgenomes with heavy gene losses (I and II) were in the same nucleus for a longer period of time than the third subgenome (III) with the fewest gene losses. This "two-step" hypothesis is essentially the same as that proposed previously for the eudicot paleohexaploidy; however, the more recent nature of the B. rapa paleohexaploidy makes this model more testable. We found that subgenome II suffered recent small deletions within exons more frequently than subgenome I, as would be expected if the genes in subgenome I had already been near maximally fractionated before subgenome III was introduced. We observed that some sequences, before these deletions, were flanked by short direct repeats, a unique signature of intrachromosomal illegitimate recombination. We also found, through simulations, that short--single or two-gene--deletions appear to dominate the fractionation patterns in B. rapa. We conclude that the observed patterns of the triplicated regions in the Brassica genome are best explained by a two-step fractionation model. The triplication and subsequent mode of fractionation could influence the potential to generate morphological diversity--a hallmark of the Brassica genus.     

Biased gene fractionation and dominant gene expression among the subgenomes of Brassica rapa

Polyploidization, both ancient and recent, is frequent among plants. A "two-step theory" was proposed to explain the meso-triplication of the Brassica "A" genome: Brassica rapa. By accurately partitioning of this genome, we observed that genes in the less fractioned subgenome (LF) were dominantly expressed over the genes in more fractioned subgenomes (MFs: MF1 and MF2), while the genes in MF1 were slightly dominantly expressed over the genes in MF2. The results indicated that the dominantly expressed genes tended to be resistant against gene fractionation. By re-sequencing two B. rapa accessions: a vegetable turnip (VT117) and a Rapid Cycling line (L144), we found that genes in LF had less non-synonymous or frameshift mutations than genes in MFs; however mutation rates were not significantly different between MF1 and MF2. The differences in gene expression patterns and on-going gene death among the three subgenomes suggest that "two-step" genome triplication and differential subgenome methylation played important roles in the genome evolution of B. rapa.     

Origin and maintenance of a broad-spectrum disease resistance locus in Arabidopsis

The broad-spectrum mildew resistance genes RPW8.1 and RPW8.2 define a unique type of plant disease resistance (R) gene, and so far homologous sequences have been found in Arabidopsis thaliana only, which suggests a recent origin. In addition to RPW8.1 and RPW8.2, the RPW8 locus contains three homologs of RPW8, HR1, HR2, and HR3, which do not contribute to powdery mildew resistance. To investigate whether RPW8 has originated recently, and if so the processes involved, we have isolated and analyzed the syntenic RPW8 loci from Arabidopsis lyrata, and from Brassica rapa and B. oleracea. The A. lyrata locus contains four genes orthologous to HR1, HR2, HR3, and RPW8.2, respectively. Two syntenic loci have been characterized in Brassica; one locus contains three genes and is present in both B. oleracea and B. rapa, and the other locus contains a single gene and is detected in B. rapa only. The Brassica homologs have highest similarity to HR3. Sequence analyses suggested that the RPW8 gene family in Brassicaceae originated from an HR3-like ancestor gene through a series of duplications and that RPW8.1 and RPW8.2 evolved from functional diversification through positive selection several MYA. Examination of the sequence polymorphism of 32 A. thaliana accessions at the RPW8 locus and their disease reaction phenotypes revealed that the polymorphic RPW8 locus defines a major source of resistance to powdery mildew diseases. A possible evolutionary mechanism by which functional polymorphism at the AtRPW8 locus has been maintained in contemporary populations of A. thaliana is discussed.     

Computational tools for Brassica-Arabidopsis comparative genomics

Recent advances, such as the availability of extensive genome survey sequence (GSS) data and draft physical maps, are radically transforming the means by which we can dissect Brassica genome structure and systematically relate it to the Arabidopsis model. Hitherto, our view of the co-linearities between these closely related genomes had been largely inferred from comparative RFLP data, necessitating substantial interpolation and expert interpretation. Sequencing of the Brassica rapa genome by the Multinational Brassica Genome Project will, however, enable an entirely computational approach to this problem. Meanwhile we have been developing databases and bioinformatics tools to support our work in Brassica comparative genomics, including a recently completed draft physical map of B. rapa integrated with anchor probes derived from the Arabidopsis genome sequence. We are also exploring new ways to display the emerging Brassica-Arabidopsis sequence homology data. We have mapped all publicly available Brassica sequences in silico to the Arabidopsis TIGR v5 genome sequence and published this in the ATIDB database that uses Generic Genome Browser (GBrowse). This in silico approach potentially identifies all paralogous sequences and so we colour-code the significance of the mappings and offer an integrated, real-time multiple alignment tool to partition them into paralogous groups. The MySQL database driving GBrowse can also be directly interrogated, using the powerful API offered by the Perl BioColon, two colonsDBColon, two colonsGFF methods, facilitating a wide range of data-mining possibilities.     

Toward unraveling the structure of Brassica rapa genome

Genomic research in any organism encompasses understanding structure of the target genome and genes, their function, and evolution. Brassica rapa , which is phylogenetically related to Arabidopsis thaliana , is an important species with respect to its uses as vegetable, oil, and fodder. The availability of suitable genetic and genomic resources is a prerequisite to undertake genomic research in B. rapa . We have developed reference mapping populations of Chinese cabbage ( B. rapa ssp. pekinensis ) comprising 78 doubled haploid lines and over 250 recombinant inbred lines. Two Bacterial Artificial Chromosome (BAC) libraries, generated by restriction enzymes Hin dIII (KBrH) and Bam HI (KBrB), comprise 56 592 and 50 688 clones, respectively. We have also constructed 22 cDNA libraries from different plant tissues consisting of 104 914 clones with an average length of 575 bp. Initial BAC-end sequence analysis of 1473 clones of the KBrH library led us to understand the structure of B. rapa genome with respect to extent of genic sequences and their annotation, and relative abundance of different types of repetitive DNAs. Full-length sequence analysis of BAC clones revealed extensive triplication of B. rapa DNA segments coupled with variable gene losses within the segments. The formulation of the 'Multinational Brassica Genome Project' has laid the foundation to sequence the complete genome of B. rapa ssp. pekinensis by the international Brassica research community. It has been proposed to undertake BAC-to-BAC sequencing of genetically mapped seed BACs. In recent years, development of bioinformatics tools in Brassica has given a boost to structural genomics research in Brassica species. The research undertaken with the availability of various genomic resources in the public domain has added to our understanding of the structure of B. rapa .     

与甘蓝型油菜重要经济性状有关的DNA克隆在拟南芥遗传图谱中的整合

拟南芥与油菜同属十字花科植物芸寡族,亲缘关系很近,基因组间的同源性很高,在用拟南芥EST克隆和油菜DNA克隆作探针定位了甘蓝型油菜一系列重要性状的基础上,对25个与油菜雄性不育恢复基因,硼高效利用基因,抗菌核病QTL及油菜种间杂种营养优势相关联的克隆进行了测序,在拟南芥基因组数据库中寻找到与这25个克隆高度同源的序列,根据这些高度同源序列在拟南芥染色体上的相位位置,将油菜DNA克隆整合到了拟南芥遗传图谱上,其中油菜硼高效基因BE1两侧的标记克隆整合在拟南芥第一染色体长臂一个较小的区段内,以该目标区段内的拟南芥EST克隆PA24为探针对甘蓝型油菜基因组比较作图,将该克隆定位在油菜连锁图BE1两侧标记之间,表明了利用基因组间的相互比较作图来精细定位芸薹属作物重要基因的可能性。     

Isolation of circadian-associated genes in Brassica rapa by comparative genomics with Arabidopsis thaliana

Elucidation of the roles of circadian associated factors requires a better understanding of the molecular mechanisms of circadian rhythms, control of flowering time through photoperiodic pathways, and photosensory signal transduction. In Arabidopsis, the APRR1 quintet, APRRs 1, 3, 5, 7, and 9, are known as central oscillator genes. Other plants may share the molecular mechanism underlying the circadian rhythm. To identify and characterize these circadian response genes in Brassica crops whose genome was triplicated after divergence from Arabidopsis, we identified B. rapa BAC clones containing these genes by BLAST analysis of B. rapa BAC end sequences against the five corresponding Arabidopsis regions. Subsequent fingerprinting, Southern hybridization, and PCR allowed identification of five BAC clones, one for each of the five circadian-related genes. By draft shotgun sequencing of the BAC clones, we identified the complete gene sequences and cloned the five expressed B. rapa circadian-associated gene members, BrPRRs 1, 3, 5, 7, and 9. Phylogenetic analysis revealed that each BrPRR was orthologous to the corresponding APRR at the sequence level. Northern hybridization revealed that the five genes were transcribed at distinct points in the 24 hour period, and Southern hybridization revealed that they are present in 2, 1, 2, 2, and 1 copies, respectively in the B. rapa genome, which was triplicated and then diploidized during the last 15 million years.     

Structural and functional comparative mapping between the Brassica A genomes in allotetraploid Brassica napus and diploid Brassica rapa

Brassica napus (AACC genome) is an important oilseed crop that was formed by the fusion of the diploids B. rapa (AA) and B. oleracea (CC). The complete genomic sequence of the Brassica A genome will be available soon from the B. rapa genome sequencing project, but it is not clear how informative the A genome sequence in B. rapa (A(r)) will be for predicting the structure and function of the A subgenome in the allotetraploid Brassica species B. napus (A(n)). In this paper, we report the results of structural and functional comparative mapping between the A subgenomes of B. napus and B. rapa based on genetic maps that were anchored with bacterial artificial chromosomes (BACs)-sequence of B. rapa. We identified segmental conservation that represented by syntenic blocks in over one third of the A genome; meanwhile, comparative mapping of quantitative trait loci for seed quality traits identified a dozen homologous regions with conserved function in the A genome of the two species. However, several genomic rearrangement events, such as inversions, intra- and inter-chromosomal translocations, were also observed, covering totally at least 5% of the A genome, between allotetraploid B. napus and diploid B. rapa. Based on these results, the A genomes of B. rapa and B. napus are mostly functionally conserved, but caution will be necessary in applying the full sequence data from B. rapa to the B. napus as a result of genomic rearrangements in the A genome between the two species.     

Comparative fluorescence in situ hybridization mapping of a 431-kb Arabidopsis thaliana bacterial artificial chromosome contig reveals the role of chromosomal duplications in the expansion of the Brassica rapa genome

Comparative genome studies are important contributors to our understanding of genome evolution. Most comparative genome studies in plants have been based on genetic mapping of homologous DNA loci in different genomes. Large-scale comparative physical mapping has been hindered by the lack of efficient and affordable techniques. We report here the adaptation of fluorescence in situ hybridization (FISH) techniques for comparative physical mapping between Arabidopsis thaliana and Brassica rapa. A set of six bacterial artificial chromosomes (BACs) representing a 431-kb contiguous region of chromosome 2 of A. thaliana was mapped on both chromosomes and DNA fibers of B. rapa. This DNA fragment has a single location in the A. thaliana genome, but hybridized to four to six B. rapa chromosomes, indicating multiple duplications in the B. rapa genome. The sizes of the fiber-FISH signals from the same BACs were not longer in B. rapa than those in A. thaliana, suggesting that this genomic region is duplicated but not expanded in the B. rapa genome. The comparative fiber-FISH mapping results support that chromosomal duplications, rather than regional expansion due to accumulation of repetitive sequences in the intergenic regions, played the major role in the evolution of the B. rapa genome.     

Genome-wide comparative analysis of the Brassica rapa gene space reveals genome shrinkage and differential loss of duplicated genes after whole genome triplication

Background

Brassica rapa is one of the most economically important vegetable crops worldwide. Owing to its agronomic importance and phylogenetic position, B. rapa provides a crucial reference to understand polyploidy-related crop genome evolution. The high degree of sequence identity and remarkably conserved genome structure between Arabidopsis and Brassica genomes enables comparative tiling sequencing using Arabidopsis sequences as references to select the counterpart regions in B. rapa, which is a strong challenge of structural and comparative crop genomics.

Results

We assembled 65.8 megabase-pairs of non-redundant euchromatic sequence of B. rapa and compared this sequence to the Arabidopsis genome to investigate chromosomal relationships, macrosynteny blocks, and microsynteny within blocks. The triplicated B. rapa genome contains only approximately twice the number of genes as in Arabidopsis because of genome shrinkage. Genome comparisons suggest that B. rapa has a distinct organization of ancestral genome blocks as a result of recent whole genome triplication followed by a unique diploidization process. A lack of the most recent whole genome duplication (3R) event in the B. rapa genome, atypical of other Brassica genomes, may account for the emergence of B. rapa from the Brassica progenitor around 8 million years ago.

Conclusions

This work demonstrates the potential of using comparative tiling sequencing for genome analysis of crop species. Based on a comparative analysis of the B. rapa sequences and the Arabidopsis genome, it appears that polyploidy and chromosomal diploidization are ongoing processes that collectively stabilize the B. rapa genome and facilitate its evolution.     

Investigating ancient duplication events in the Arabidopsis genome

The complete genomic analysis of Arabidopsis thaliana has shown that a major fraction of the genome consists of paralogous genes that probably originated through one or more ancient large-scale gene or genome duplication events. However, the number and timing of these duplications still remains unclear, and several different hypotheses have been put forward recently. Here, we reanalyzed duplicated blocks found in the Arabidopsis genome described previously and determined their date of divergence based on silent substitution estimations between the paralogous genes and, where possible, by phylogenetic reconstruction. We show that methods based on averaging protein distances of heterogeneous classes of duplicated genes lead to unreliable conclusions and that a large fraction of blocks duplicated much more recently than assumed previously. We found clear evidence for one large-scale gene or even complete genome duplication event somewhere between 70 to 90 million years ago. Traces pointing to a much older (probably more than 200 million years) large-scale gene duplication event could be detected. However, for now it is impossible to conclude whether these old duplicates are the result of one or more large-scale gene duplication events. abbreviations dA, fraction of amino acid substitutions; Kn, number of nonsynonymous substitutions per nonsynonymous site; Ks, number of synonymous substitutions per synonymous site; MYA, million years ago     

The CACTA transposon Bot1 played a major role in Brassica genome divergence and gene proliferation

We isolated and characterized a Brassica C genome-specific CACTA element, which was designated Bot1 (Brassica oleracea transposon 1). After analysing phylogenetic relationships, copy numbers and sequence similarity of Bot1 and Bot1 analogues in B. oleracea (C genome) versus Brassica rapa (A genome), we concluded that Bot1 has encountered several rounds of amplification in the oleracea genome only, and has played a major role in the recent rapa and oleracea genome divergence. We performed in silico analyses of the genomic organization and internal structure of Bot1, and established which segment of Bot1 is C-genome specific. Our work reports a fully characterized Brassica repetitive sequence that can distinguish the Brassica A and C chromosomes in the allotetraploid Brassica napus, by fluorescent in situ hybridization. We demonstrated that Bot1 carries a host S locus-associated SLL3 gene copy. We speculate that Bot1 was involved in the proliferation of SLL3 around the Brassica genome. The present study reinforces the assumption that transposons are a major driver of genome and gene evolution in higher plants.     

Multiple flowering time QTLs within several Brassica species could be the result of duplicated copies of one ancestral gene.

Quantitative trait locus (QTL) analysis was used to study the evolution of genes controlling the timing of flowering in four Brassica genomes that are all extensively replicated. Comparative mapping showed that a chromosomal region from the top of Arabidopsis thaliana chromosome 5 corresponded to three homoeologous copies in each of the diploid species Brassica nigra, B. oleracea, and B. rapa and six copies in the amphidiploid B. juncea. QTLs were detected in two of the three replicated segments in each diploid genome and in three of the six replicated segments in B. juncea. These results indicate that, for the studied trait, multiple QTLs resulting from genome duplication is the rule rather than the exception. Brassica homologues to two candidate genes (CO and FLC) identified from the corresponding A. thaliana region were mapped. CO homologues mapped close to the QTL peaks in eight of nine QTLs, while FLC homologues mapped farther away in those cases where the mapping resolution allowed a comparison. Thus, our data are consistent with the hypothesis that all the major QTLs we detected in the different species of Brassica could be the result of duplicated copies of the same ancestral gene, possibly the ancestor of CO.     

Analysis of expressed sequence tags from Brassica rapa L. ssp. pekinensis

Non-redundant expressed sequence tags (ESTs) were generated from six different organs at various developmental stages of Chinese cabbage, Brassica rapa L. ssp. pekinensis. Of the 1,295 ESTs, 915 (71%) showed significantly high homology in nucleotide or deduced amino acid sequences with other sequences deposited in databases, while 380 did not show similarity to any sequences. Briefly, 598 ESTs matched with proteins of identified biological function, 177 with hypothetical proteins or non-annotated Arabidopsis genome sequences, and 140 with other ESTs. About 82% of the top-scored matching sequences were from Arabidopsis or Brassica, but overall 558 (43%) ESTs matched with Arabidopsis ESTs at the nucleotide sequence level. This observation strongly supports the idea that gene-expression profiles of Chinese cabbage differ from that of Arabidopsis, despite their genome structures being similar to each other. Moreover, sequence analyses of 21 Brassica ESTs revealed that their primary structure is different from those of corresponding annotated sequences of Arabidopsis genes. Our data suggest that direct prediction of Brassica gene expression pattern based on the information from Arabidopsis genome research has some limitations. Thus, information obtained from the Brassica EST study is useful not only for understanding of unique developmental processes of the plant, but also for the study of Arabidopsis genome structure.