Mina Jin

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

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

BackgroundBrassica 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.ResultsWe 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.ConclusionsThis 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.

A Sequence-Tagged Linkage Map of Brassica rapa

A detailed genetic linkage map of Brassica rapa has been constructed containing 545 sequence-tagged loci covering 1287 cM, with an average mapping interval of 2.4 cM. The loci were identified using a combination of 520 RFLP and 25 PCR-based markers. RFLP probes were derived from 359 B. rapa EST clones and amplification products of 11 B. rapa and 26 Arabidopsis. Including 21 SSR markers provided anchors to previously published linkage maps for B. rapa and B. napus and is followed as the referenced mapping of R1–R10. The sequence-tagged markers allowed interpretation of the pattern of chromosome duplications within the B. rapa genome and comparison with Arabidopsis. A total of 62 EST markers showing a single RFLP band were mapped through 10 linkage groups, indicating that these can be valuable anchoring markers for chromosome-based genome sequencing of B. rapa. Other RFLP probes gave rise to 2–5 loci, inferring that B. rapa genome duplication is a general phenomenon through 10 chromosomes. The map includes five loci of FLC paralogues, which represent the previously reported BrFLC-1, -2, -3, and -5 and additionally identified BrFLC3 paralogues derived from local segmental duplication on R3.

Sequenced BAC anchored reference genetic map that reconciles the ten individual chromosomes of Brassica rapa.

BackgroundIn view of the immense value of Brassica rapa in the fields of agriculture and molecular biology, the multinational Brassica rapa Genome Sequencing Project (BrGSP) was launched in 2003 by five countries. The developing BrGSP has valuable resources for the community, including a reference genetic map and seed BAC sequences. Although the initial B. rapa linkage map served as a reference for the BrGSP, there was ambiguity in reconciling the linkage groups with the ten chromosomes of B. rapa. Consequently, the BrGSP assigned each of the linkage groups to the project members as chromosome substitutes for sequencing.ResultsWe identified simple sequence repeat (SSR) motifs in the B. rapa genome with the sequences of seed BACs used for the BrGSP. By testing 749 amplicons containing SSR motifs, we identified polymorphisms that enabled the anchoring of 188 BACs onto the B. rapa reference linkage map consisting of 719 loci in the 10 linkage groups with an average distance of 1.6 cM between adjacent loci. The anchored BAC sequences enabled the identification of 30 blocks of conserved synteny, totaling 534.9 cM in length, between the genomes of B. rapa and Arabidopsis thaliana. Most of these were consistent with previously reported duplication and rearrangement events that differentiate these genomes. However, we were able to identify the collinear regions for seven additional previously uncharacterized sections of the A genome. Integration of the linkage map with the B. rapa cytogenetic map was accomplished by FISH with probes representing 20 BAC clones, along with probes for rDNA and centromeric repeat sequences. This integration enabled unambiguous alignment and orientation of the maps representing the 10 B. rapa chromosomes.ConclusionWe developed a second generation reference linkage map for B. rapa, which was aligned unambiguously to the B. rapa cytogenetic map. Furthermore, using our data, we confirmed and extended the comparative genome analysis between B. rapa and A. thaliana. This work will serve as a basis for integrating the genetic, physical, and chromosome maps of the BrGSP, as well as for studies on polyploidization, speciation, and genome duplication in the genus Brassica.

Genome-wide identification of glucosinolate synthesis genes in Brassica rapa.

Glucosinolates play important roles in plant defense against herbivores and microbes, as well as in human nutrition. Some glucosinolate‐derived isothiocyanate and nitrile compounds have been clinically proven for their anticarcinogenic activity. To better understand glucosinolate biosynthesis in Brassica rapa, we conducted a comparative genomics study with Arabidopsis thaliana and identified total 56 putative biosynthetic and regulator genes. This established a high colinearity in the glucosinolate biosynthesis pathway between Arabidopsis and B. rapa. Glucosinolate genes in B. rapa share 72–94% nucleotide sequence identity with the Arabidopsis orthologs and exist in different copy numbers. The exon/intron split pattern of B. rapa is almost identical to that of Arabidopsis, although inversion, insertion, deletion and intron size variations commonly occur. Four genes appear to be nonfunctional as a result of the presence of a frame shift mutation and retrotransposon insertion. At least 12 paralogs of desulfoglucosinolate sulfotransferase were found in B. rapa, whereas only three were found in Arabidopsis. The expression of those paralogs was not tissue‐specific but varied greatly depending on B. rapa tissue types. Expression was also developmentally regulated in some paralogs but not in other paralogs. Most of the regulator genes are present as triple copies. Accordingly, glucosinolate synthesis and regulation in B. rapa appears to be more complex than that of Arabidopsis. With the isolation and further characterization of the endogenous genes, health‐beneficial vegetables or desirable animal feed crops could be developed by metabolically engineering the glucosinolate pathway.

The Korea Brassica Genome Project: A glimpse of the Brassica genome based on comparative genome analysis with Arabidopsis

A complete genome sequence provides unlimited information in the sequenced organism as well as in related taxa. According to the guidance of the Multinational Brassica Genome Project (MBGP), the Korea Brassica Genome Project (KBGP) is sequencing chromosome 1 (cytogenetically oriented chromosome #1) of Brassica rapa. We have selected 48 seed BACs on chromosome 1 using EST genetic markers and FISH analyses. Among them, 30 BAC clones have been sequenced and 18 are on the way. Comparative genome analyses of the EST sequences and sequenced BAC clones from Brassica chromosome 1 revealed their homeologous partner regions on the Arabidopsis genome and a syntenic comparative map between Brassica chromosome 1 and Arabidopsis chromosomes. In silico chromosome walking and clone validation have been successfully applied to extending sequence contigs based on the comparative map and BAC end sequences. In addition, we have defined the (peri)centromeric heterochromatin blocks with centromeric tandem repeats, rDNA and centromeric retrotransposons. In-depth sequence analyses of five homeologous BAC clones and an Arabidopsis chromosomal region reveal overall co-linearity, with 82% sequence similarity. The data indicate that the Brassica genome has undergone triplication and subsequent gene losses after the divergence of Arabidopsis and Brassica. Based on in-depth comparative genome analyses, we propose a comparative genomics approach for conquering the Brassica genome. In 2005 we intend to construct an integrated physical map, including sequence information from 500 BAC clones and integration of fingerprinting data and end sequence data of more than 100 000 BAC clones. The sequences have been submitted to GenBank with accession numbers: 10 204 BAC ends of the KBrH library (CW978640–CW988843); KBrH138P04, AC155338; KBrH117N09, AC155337; KBrH097M21, AC155348; KBrH093K03, AC155347; KBrH081N08, AC155346; KBrH080L24, AC155345; KBrH077A05, AC155343; KBrH020D15, AC155340; KBrH015H17, AC155339; KBrH001H24, AC155335; KBrH080A08, AC155344; KBrH004D11, AC155341; KBrH117M18, AC146875; KBrH052O08, AC155342.

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 (Ar) will be for predicting the structure and function of the A subgenome in the allotetraploid Brassica species B. napus (An). 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.

A Systems Genetics Approach Identifies Gene Regulatory Networks Associated with Fatty Acid Composition in Brassica rapa Seed

A systems genetics approach identifies gene regulatory networks associated with fatty acid composition in Brassica rapa seed. Fatty acids in seeds affect seed germination and seedling vigor, and fatty acid composition determines the quality of seed oil. In this study, quantitative trait locus (QTL) mapping of fatty acid and transcript abundance was integrated with gene network analysis to unravel the genetic regulation of seed fatty acid composition in a Brassica rapa doubled haploid population from a cross between a yellow sarson oil type and a black-seeded pak choi. The distribution of major QTLs for fatty acids showed a relationship with the fatty acid types: linkage group A03 for monounsaturated fatty acids, A04 for saturated fatty acids, and A05 for polyunsaturated fatty acids. Using a genetical genomics approach, expression quantitative trait locus (eQTL) hotspots were found at major fatty acid QTLs on linkage groups A03, A04, A05, and A09. An eQTL-guided gene coexpression network of lipid metabolism-related genes showed major hubs at the genes BrPLA2-ALPHA, BrWD-40, a number of seed storage protein genes, and the transcription factor BrMD-2, suggesting essential roles for these genes in lipid metabolism. Three subnetworks were extracted for the economically important and most abundant fatty acids erucic, oleic, linoleic, and linolenic acids. Network analysis, combined with comparison of the genome positions of cis- or trans-eQTLs with fatty acid QTLs, allowed the identification of candidate genes for genetic regulation of these fatty acids. The generated insights in the genetic architecture of fatty acid composition and the underlying complex gene regulatory networks in B. rapa seeds are discussed.

Construction of random sheared fosmid library from Chinese cabbage and its use for Brassica rapa genome sequencing project.

As a part of the Multinational Genome Sequencing Project of Brassica rapa, linkage group R9 and R3 were sequenced using a bacterial artificial chromosome (BAC) by BAC strategy. The current physical contigs are expected to cover approximately 90% euchromatins of both chromosomes. As the project progresses, BAC selection for sequence extension becomes more limited because BAC libraries are restriction enzyme-specific. To support the project, a random sheared fosmid library was constructed. The library consists of 97536 clones with average insert size of approximately 40 kb corresponding to seven genome equivalents, assuming a Chinese cabbage genome size of 550 Mb. The library was screened with primers designed at the end of sequences of nine points of scaffold gaps where BAC clones cannot be selected to extend the physical contigs. The selected positive clones were end-sequenced to check the overlap between the fosmid clones and the adjacent BAC clones. Nine fosmid clones were selected and fully sequenced. The sequences revealed two completed gap filling and seven sequence extensions, which can be used for further selection of BAC clones confirming that the fosmid library will facilitate the sequence completion of B. rapa.

Molecular characterization of glucosinolates and carotenoid biosynthetic genes in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

The present study aimed to investigate the contents of glucosinolates (GSLs) and carotenoids in eleven varieties of Chinese cabbage in relation to the expression level of the important transcription factors. MS and HPLC analysis identified the presence of 13 GSLs (progoitrin, sinigrin, glucoalyssin, gluconapoleiferin, gluconapin, glucocochlearin, glucobrassicanapin, glucoerucin, 4-hydroxyglucobrassicin, glucobrassicin, 4-methoxyglucobrassicin, neoglucobrassicin and gluconasturtiin) and four carotenoids (lutein, zeaxanthin, α-carotene and β-carotene). GSL contents were varied among the different cabbage varieties. The total GSL content ranged from 2.7 to 57.88 μmol/g DW. The proportion of gluconapin (54%) and glucobrassicanapin (22%) was higher in all the varieties, respectively. Results documented the variation in total and individual carotenoid contents that have also been observed among different varieties; however, the total carotenoid contents ranged from 289.12 to 1001.41 mg kg−1 DW (mean 467.66). Interestingly, the proportion of lutein (66.5) and β-carotene (25.9) were higher than α-carotene (5.1) and zeaxanthin (2.5%). Consequently, the expression level of the regulatory gene, MYB28 was higher in ‘K0648’ and was directly proportional to GSL content. Similarly, the expression levels of 1-PSY were higher in ‘K0112’; however, the expression levels of 2-ZDS, 3-LCYB, 4-LCYE, 5-CHXB and 7-NCED genes showed no significant difference. In addition, the correlation between GSL and carotenoid contents and gene expression level showed moderate significant difference in each Chinese cabbage.