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.     

Sequence-level analysis of the diploidization process in the triplicated FLOWERING LOCUS C region of Brassica rapa

Strong evidence exists for polyploidy having occurred during the evolution of the tribe Brassiceae. We show evidence for the dynamic and ongoing diploidization process by comparative analysis of the sequences of four paralogous Brassica rapa BAC clones and the homologous 124-kb segment of Arabidopsis thaliana chromosome 5. We estimated the times since divergence of the paralogous and homologous lineages. The three paralogous subgenomes of B. rapa triplicated 13 to 17 million years ago (MYA), very soon after the Arabidopsis and Brassica divergence occurred at 17 to 18 MYA. In addition, a pair of BACs represents a more recent segmental duplication, which occurred approximately 0.8 MYA, and provides an exception to the general expectation of three paralogous segments within the B. rapa genome. The Brassica genome segments show extensive interspersed gene loss relative to the inferred structure of the ancestral genome, whereas the Arabidopsis genome segment appears little changed. Representatives of all 32 genes in the Arabidopsis genome segment are represented in Brassica, but the hexaploid complement of 96 has been reduced to 54 in the three subgenomes, with compression of the genomic region lengths they occupy to between 52 and 110 kb. The gene content of the recently duplicated B. rapa genome segments is identical, but intergenic sequences differ.     

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.     

Characterization of DNA methyltransferase genes in Brassica rapa

DNA methylation is essential for normal development and plays important roles in regulating gene expression in plants. Analysis of the key enzymes catalyzing DNA methylation is important to understand epigenetic phenomena. In this study, three putative methyltransferase genes, BrMET1a, BrMET1b, and BrCMT, were isolated from a genome library of Brassica rapa. Structural conservation of the amino acid sequence between BrMET1a/BrMET1b and AtMET1 and that between BrCMT and AtCMT3 suggests that they may function as DNA methyltransferase. BrMET1a was expressed in vegetative and reproductive organs, while BrMET1b was expressed only in pistils, indicating that these two genes have different functions. BrCMT was expressed especially in stamens at the stage of 2-4 days before anthesis. We isolated three DNA methyltransferase genes in Brassica rapa and indicated differences of expression patterns of these DNA methyltransferase genes and expression levels in different tissues and developmental stages, suggesting that these genes might play important roles in epigenetic gene regulation in B. rapa.     

白菜和甘蓝基因组转座子表达及其对基因调控的潜在影响

转座子或转座元件是大多数真核生物基因组的主要组成成分。甘蓝(Brassica oleracea)基因组比白菜(B. rapa)大主要是转座子的扩增差异造成的。然而, 这两个芸薹属近缘物种转座子表达水平以及对基因的调控和功能的影响目前还不清楚。文章对白菜和甘蓝叶、根、茎3个器官的转录组数据进行了初步分析。结果显示, 转座子的表达量很低, 转录组reads中有1%来自转座子的转录本; 转座子的表达存在器官差异, 且不同类别和家族的转座子表达量相差很大, 相同类别和同一家族的转座子在白菜和甘蓝基因组中的表达活性也不相同。进一步鉴定到转录读出的LTR反转座子, 其与下游基因距离小于2 kb的有41个, 小于100 bp的有9个, 这些LTR的转录读出很可能通过正义或反义的转录本激活或干扰下游基因的表达。同时, 具有转录读出的intact LTR比solo LTR具有更强的读出活性。通过深入分析转座子的插入位点发现, 白菜基因组中转座子插入基因内部的频率比甘蓝基因组中的高; 与反转座子相比, DNA转座子更偏向于插入或保留在基因的内含子当中。这些结果为认识转座子对其他蛋白编码基因的影响提供了基础。     

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.     

通过白菜型油菜和埃塞俄比亚芥的Ar和Cc基因组导入创造甘蓝型油菜新材料

油菜是目前我国主要种植的油料作物之一,但现有的种质资源限制了产量的进一步提高。本研究采取了一种新的育种方式来增加甘蓝型油菜的种质资源,即通过远缘杂交结合分子标记辅助选择的方式将白菜型油菜的Ar基因组和埃塞俄比亚芥的Cc对现有的甘蓝型油菜品种的基因组(AnAnCnCn)进行部分替换。通过对五倍体杂交后代(ArAnBcCcCn)进行染色体选择,找到了染色体数目为38的材料。为了和现有的甘蓝型油菜进行区分,得到的新材料被认定为甘蓝型油菜新材料。实验结果表明,得到的部分甘蓝型油菜新材料具有基本正常的减数分裂过程、正常的花粉萌发以及胚囊发育过程,这说明甘蓝型油菜新材料达到了遗传平衡。分子标记分析表明:甘蓝型油菜新材料的约50%的基因组被白菜型油菜的Ar基因组和埃塞俄比亚芥的Cc替换,并且这些甘蓝型油菜新材料之间具有丰富的遗传多样性。因此,白菜型油菜的Ar基因组和埃塞俄比亚芥的Cc基因组导入对于丰富现有的甘蓝型油菜种质资源具有明显的效果。     

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.     

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.     

Genes identified by visible mutant phenotypes show increased bias toward one of two subgenomes of maize

Not all genes are created equal. Despite being supported by sequence conservation and expression data, knockout homozygotes of many genes show no visible effects, at least under laboratory conditions. We have identified a set of maize (Zea mays L.) genes which have been the subject of a disproportionate share of publications recorded at MaizeGDB. We manually anchored these "classical" maize genes to gene models in the B73 reference genome, and identified syntenic orthologs in other grass genomes. In addition to proofing the most recent version 2 maize gene models, we show that a subset of these genes, those that were identified by morphological phenotype prior to cloning, are retained at syntenic locations throughout the grasses at much higher levels than the average expressed maize gene, and are preferentially found on the maize1 subgenome even with a duplicate copy is still retained on the opposite subgenome. Maize1 is the subgenome that experienced less gene loss following the whole genome duplication in maize lineage 5-12 million years ago and genes located on this subgenome tend to be expressed at higher levels in modern maize. Links to the web based software that supported our syntenic analyses in the grasses should empower further research and support teaching involving the history of maize genetic research. Our findings exemplify the concept of "grasses as a single genetic system," where what is learned in one grass may be applied to another.     

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.     

Stable Bacillus thuringiensis (Bt) toxin content in interspecific F1 and backcross populations of wild Brassica rapa after Bt gene transfer

Stable expression of a transgene may lead to increased fitness for wild plants after acquiring the transgene via crop-weed hybridization. Here, we investigate the stability of Bt toxin content in wild Brassica rapa acquiring the Bt gene from Bt Brassica napus. The Bt toxin content in nine Bt-expressing B. napus lines was 0.80-1.70 micro g/g leaf tissue throughout the growing season. These nine lines were crossed with three accessions of wild B. rapa and the Bt gene was successfully transferred to interspecific hybrids (F1) and successive backcross generations (BC1 to BC4). The Bt toxin level in F1 and BC progenies containing the Bt gene remained at 0.90-3.10 micro g/g leaf tissue. This study indicates that the Bt gene can persist and be stably expressed in wild B. rapa.     

Genome-wide characterization of the pectate lyase-like (PLL) genes in Brassica rapa

Pectate lyases (PL) depolymerize demethylated pectin (pectate, EC 4.2.2.2) by catalyzing the eliminative cleavage of α-1,4-glycosidic linked galacturonan. Pectate Lyase-like (PLL) genes are one of the largest and most complex families in plants. However, studies on the phylogeny, gene structure, and expression of PLL genes are limited. To understand the potential functions of PLL genes in plants, we characterized their intron–exon structure, phylogenetic relationships, and protein structures, and measured their expression patterns in various tissues, specifically the reproductive tissues in Brassica rapa. Sequence alignments revealed two characteristic motifs in PLL genes. The chromosome location analysis indicated that 18 of the 46 PLL genes were located in the least fractionated sub-genome (LF) of B. rapa, while 16 were located in the medium fractionated sub-genome (MF1) and 12 in the more fractionated sub-genome (MF2). Quantitative RT-PCR analysis showed that BrPLL genes were expressed in various tissues, with most of them being expressed in flowers. Detailed qRT-PCR analysis identified 11 pollen specific PLL genes and several other genes with unique spatial expression patterns. In addition, some duplicated genes showed similar expression patterns. The phylogenetic analysis identified three PLL gene subfamilies in plants, among which subfamily II might have evolved from gene neofunctionalization or subfunctionalization. Therefore, this study opens the possibility for exploring the roles of PLL genes during plant development.     

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 .     

Expression, mapping, and genetic variability of Brassica napus disease resistance gene analogues.

Numerous sequences analogous to resistance (R) genes exist in plant genomes and could be involved in resistance traits. The aim of this study was to identify a large number of Brassica napus sequences related to R genes and also to test the adequacy of specific PCR-based tools for studying them. Different consensus primers were compared for their efficiency in amplifying resistance-gene analogues (RGAs) related to the nucleotide-binding-site subgroup of R genes. Specific primers were subsequently designed to fine-study the different RGAs and we tested their efficiency in three species related to B. napus: Brassica oleracea, Brassica rapa, and Arabidopsis thaliana. Forty-four B. napus RGAs were identified. Among 29 examined, at least one-third were expressed. Eighteen RGAs were mapped on 10 of the 19 B. napus linkage groups. The high variability within these sequences permitted discrimination of each genotype within a B. napus collection. The RGA-specific primers amplified RGAs in the B. oleracea and B. rapa genomes, but the sequences appear to be poorly conserved in A. thaliana. Specific RGA primers are a precise tool for studying known-sequence RGAs. These sequences represent interesting markers that could be correlated with resistance traits in B. napus or related Brassica genomes.     

Comparison of Flowering Time Genes in Brassica Rapa, B. Napus and Arabidopsis Thaliana

The major difference between annual and biennial cultivars of oilseed Brassica napus and B. rapa is conferred by genes controlling vernalization-responsive flowering time. These genes were compared between the species by aligning the map positions of flowering time quantitative trait loci (QTLs) detected in a segregating population of each species. The results suggest that two major QTLs identified in B. rapa correspond to two major QTLs identified in B. napus. Since B. rapa is one of the hypothesized diploid parents of the amphidiploid B. napus, the vernalization requirement of B. napus probably originated from B. rapa. Brassica genes also were compared to flowering time genes in Arabidopsis thaliana by mapping RFLP loci with the same probes in both B. napus and Arabidopsis. The region containing one pair of Brassica QTLs was collinear with the top of chromosome 5 in A. thaliana where flowering time genes FLC, FY and CO are located. The region containing the second pair of QTLs showed fractured collinearity with several regions of the Arabidopsis genome, including the top of chromosome 4 where FRI is located. Thus, these Brassica genes may correspond to two genes (FLC and FRI) that regulate flowering time in the latest flowering ecotypes of Arabidopsis.