Chinese cabbage (Brassica rapa ssp. pekinensis) is a member of one of the most important leaf vegetables grown worldwide, which has experienced thousands of years in cultivation and artificial selection. The entire Chinese cabbage genome sequence, and more than forty thousand proteins have been obtained to date. The genome has undergone triplication events since its divergence from Arabidopsis thaliana (13 to 17 Mya), however a high degree of sequence similarity and conserved genome structure remain between the two species. Arabidopsis is therefore a viable reference species for comparative genomics studies. Variation in the number of members in gene families due to genome triplication may contribute to the broad range of phenotypic plasticity, and increased tolerance to environmental extremes observed in Brassica species. Transcription factors are important regulators involved in plant developmental and physiological processes. The AP2/ERF proteins, one of the most important families of transcriptional regulators, play a crucial role in plant growth, and in response to biotic and abiotic stressors. Our analysis will provide resources for understanding the tolerance mechanisms in Brassica rapa ssp. pekinensis.
In the present study, 291 putative AP2/ERF transcription factor proteins were identified from the Chinese cabbage genome database, and compared with proteins from 15 additional species. The Chinese cabbage AP2/ERF superfamily was classified into four families, including AP2, ERF, RAV, and Soloist. The ERF family was further divided into DREB and ERF subfamilies. The AP2/ERF superfamily was subsequently divided into 15 groups. The identification, classification, phylogenetic reconstruction, conserved motifs, chromosome distribution, functional annotation, expression patterns, and interaction networks of the AP2/ERF transcription factor superfamily were predicted and analyzed. Distribution mapping results showed AP2/ERF superfamily genes were localized on the 10 Chinese cabbage chromosomes. AP2/ERF transcription factor expression levels exhibited differences among six tissue types based on expressed sequence tags (ESTs). In the AP2/ERF superfamily, 214 orthologous genes were identified between Chinese cabbage and Arabidopsis. Orthologous gene interaction networks were constructed, and included seven CBF and four AP2 genes, primarily involved in cold regulatory pathways and ovule development, respectively.
The evolution of the AP2/ERF transcription factor superfamily in Chinese cabbage resulted from genome triplication and tandem duplications. A comprehensive analysis of the physiological functions and biological roles of AP2/ERF superfamily genes in Chinese cabbage is required to fully elucidate AP2/ERF, which provides us with rich resources and opportunities to understand crop stress tolerance mechanisms.
Chinese cabbage; AP2/ERF; Stress tolerance; Gene expression; Interaction network; Protein annotation
The western flower thrips (Frankliniella occidentalis [Pergande]) is one of the most important insect herbivores of cultivated plants. However, no pesticide provides complete control of this species, and insecticide resistance has emerged around the world. We previously reported the important role of jasmonate (JA) in the plant's immediate response to thrips feeding by using an Arabidopsis leaf disc system. In this study, as the first step toward practical use of JA in thrips control, we analyzed the effect of JA-regulated Arabidopsis defense at the whole plant level on thrips behavior and life cycle at the population level over an extended period. We also studied the effectiveness of JA-regulated plant defense on thrips damage in Chinese cabbage (Brassica rapa subsp. pekinensis).
Thrips oviposited more on Arabidopsis JA-insensitive coi1-1 mutants than on WT plants, and the population density of the following thrips generation increased on coi1-1 mutants. Moreover, thrips preferred coi1-1 mutants more than WT plants. Application of JA to WT plants before thrips attack decreased the thrips population. To analyze these important functions of JA in a brassica crop plant, we analyzed the expression of marker genes for JA response in B. rapa. Thrips feeding induced expression of these marker genes and significantly increased the JA content in B. rapa. Application of JA to B. rapa enhanced plant resistance to thrips, restricted oviposition, and reduced the population density of the following generation.
Our results indicate that the JA-regulated plant defense restricts thrips performance and preference, and plays an important role in the resistance of Arabidopsis and B. rapa to thrips damage.
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).
Separated by only 10-20 million years, Brassica species and Arabidopsis thaliana are closely related, and it is expected that knowledge obtained relating to Arabidopsis growth and development can be translated into Brassicas for crop improvement. Moreover, certain aspects of plant development are sufficiently different between Brassica and Arabidopsis to warrant studies to be carried out directly in the crop species. However, mutating individual genes in the amphidiploid Brassicas such as B. napus and B. juncea may, on the other hand, not give rise to expected phenotypes as the genomes of these species can contain up to six orthologues per single-copy Arabidopsis gene. In order to elucidate and possibly exploit the function of redundant genes for oilseed rape crop improvement, it may therefore be more efficient to study the effects in one of the diploid Brassica species such as B. rapa. Moreover, the ongoing sequencing of the B. rapa genome makes this species a highly attractive model for Brassica research and genetic resource development.
Seeds from the diploid Brassica A genome species, B. rapa were treated with ethyl methane sulfonate (EMS) to produce a TILLING (Targeting Induced Local Lesions In Genomes) population for reverse genetics studies. We used the B. rapa genotype, R-o-18, which has a similar developmental ontogeny to an oilseed rape crop. Hence this resource is expected to be well suited for studying traits with relevance to yield and quality of oilseed rape. DNA was isolated from a total of 9,216 M2 plants and pooled to form the basis of the TILLING platform. Analysis of six genes revealed a high level of mutations with a density of about one per 60 kb. This analysis also demonstrated that screening a 1 kb amplicon in just one third of the population (3072 M2 plants) will provide an average of 68 mutations and a 97% probability of obtaining a stop-codon mutation resulting in a truncated protein. We furthermore calculated that each plant contains on average ~10,000 mutations and due to the large number of plants, it is predicted that mutations in approximately half of the GC base pairs in the genome exist within this population.
We have developed the first EMS TILLING resource in the diploid Brassica species, B. rapa. The mutation density in this population is ~1 per 60 kb, which makes it the most densely mutated diploid organism for which a TILLING population has been published. This resource is publicly available through the RevGenUK reverse genetics platform http://revgenuk.jic.ac.uk.
As part of a research programme focused on flavonoid biosynthesis in the seed coat of Brassica napus L. (oilseed rape), orthologs of the BANYULS gene that encoded anthocyanidin reductase were cloned in B. napus as well as in the related species Brassica rapa and Brassica oleracea. B. napus genome contained four functional copies of BAN, two originating from each diploid progenitor. Amino acid sequences were highly conserved between the Brassicaceae including B. napus, B. rapa, B. oleracea as well as the model plant Arabidopsis thaliana. Along the 200 bp in 5′ of the ATG codon, Bna.BAN promoters (ProBna.BAN) were conserved with AtANR promoter and contained putative cis-acting elements. In addition, transgenic Arabidopsis and oilseed rape plants carrying the first 230 bp of ProBna.BAN fused to the UidA reporter gene were generated. In the two Brassicaceae backgrounds, ProBna.BAN activity was restricted to the seed coat. In B. napus seed, ProBna.BAN was activated in procyanidin-accumulating cells, namely the innermost layer of the inner integument and the micropyle-chalaza area. At the transcriptional level, the four Bna.BAN genes were expressed in the seed. Laser microdissection assays of the seed integuments showed that Bna.BAN expression was restricted to the inner integument, which was consistent with the activation profile of ProBna.BAN. Finally, Bna.BAN genes were mapped onto oilseed rape genetic maps and potential co-localisations with seed colour quantitative trait loci are discussed.
Anthocyanidin reductase; BANYULS genes; Brassica; Flavonoid metabolism; Seed coat-specific promoter
Brassica rapa, which is closely related to
Arabidopsis thaliana, is an important crop and a
model plant for studying genome evolution via
polyploidization. We report the current understanding of the
genome structure of B. rapa and efforts for the
whole-genome sequencing of the species. The tribe
Brassicaceae, which comprises ca. 240 species,
descended from a common hexaploid ancestor with a basic genome
similar to that of Arabidopsis. Chromosome
rearrangements, including fusions and/or fissions, resulted in
the present-day “diploid” Brassica
species with variation in chromosome number and phenotype.
Triplicated genomic segments of B. rapa are
collinear to those of A. thaliana with InDels.
The genome triplication has led to an approximately 1.7-fold
increase in the B. rapa gene number compared to
that of A. thaliana. Repetitive DNA of B.
rapa has also been extensively amplified and has
diverged from that of A. thaliana. For its
whole-genome sequencing, the Brassica rapa Genome
Sequencing Project (BrGSP) consortium has developed suitable
genomic resources and constructed genetic and physical maps.
Ten chromosomes of B. rapa are being allocated to
BrGSP consortium participants, and each chromosome will be
sequenced by a BAC-by-BAC approach. Genome sequencing of
B. rapa will offer a new perspective for plant
biology and evolution in the context of polyploidization.
Chromosomal synteny analysis is important in genome comparison to reveal genomic evolution of related species. Shared synteny describes genomic fragments from different species that originated from an identical ancestor. Syntenic genes are orthologs located in these syntenic fragments, so they often share similar functions. Syntenic gene analysis is very important in Brassicaceae species to share gene annotations and investigate genome evolution. Here we designed and developed a direct and efficient tool, SynOrths, to identify pairwise syntenic genes between genomes of Brassicaceae species. SynOrths determines whether two genes are a conserved syntenic pair based not only on their sequence similarity, but also by the support of homologous flanking genes. Syntenic genes between Arabidopsis thaliana and Brassica rapa, Arabidopsis lyrata and B. rapa, and Thellungiella parvula and B. rapa were then identified using SynOrths. The occurrence of genome triplication in B. rapa was clearly observed, many genes that were evenly distributed in the genomes of A. thaliana, A. lyrata, and T. parvula had three syntenic copies in B. rapa. Additionally, there were many B. rapa genes that had no syntenic orthologs in A. thaliana, but some of these had syntenic orthologs in A. lyrata or T. parvula. Only 5,851 genes in B. rapa had no syntenic counterparts in any of the other three species. These 5,851 genes could have originated after B. rapa diverged from these species. A tool for syntenic gene analysis between species of Brassicaceae was developed, SynOrths, which could be used to accurately identify syntenic genes in differentiated but closely-related genomes. With this tool, we identified syntenic gene sets between B. rapa and each of A. thaliana, A. lyrata, T. parvula. Syntenic gene analysis is important for not only the gene annotation of newly sequenced Brassicaceae genomes by bridging them to model plant A. thaliana, but also the study of genome evolution in these species.
synteny; ortholog; Brassica rapa; Arabidopsis thaliana; Arabidopsis lyrata; Thellugiella parvula; Brassicaceae
Brassica rapa (AA) contains very diverse forms which include oleiferous types and many vegetable types. Genome sequence of B. rapa line Chiifu (ssp. pekinensis), a leafy vegetable type, was published in 2011. Using this knowledge, it is important to develop genomic resources for the oleiferous types of B. rapa. This will allow more involved molecular mapping, in-depth study of molecular mechanisms underlying important agronomic traits and introgression of traits from B. rapa to major oilseed crops - B. juncea (AABB) and B. napus (AACC). The study explores the availability of SNPs in RNA-seq generated contigs of three oleiferous lines of B. rapa - Candle (ssp. oleifera, turnip rape), YSPB-24 and Tetra (ssp. trilocularis, Yellow sarson) and their use in genome-wide linkage mapping and specific-region fine mapping using a RIL population between Chiifu and Tetra.
RNA-seq was carried out on the RNA isolated from young inflorescences containing unopened floral buds, floral axis and small leaves, using Illumina paired-end sequencing technology. Sequence assembly was carried out using the Velvet de-novo programme and the assembled contigs were organised against Chiifu gene models, available in the BRAD-CDS database. RNA-seq confirmed the presence of more than 17,000 single-copy gene models described in the BRAD database. The assembled contigs and the BRAD gene models were analyzed for the presence of SSRs and SNPs. While the number of SSRs was limited, more than 0.2 million SNPs were observed between Chiifu and the three oleiferous lines. Assays for SNPs were designed using KASPar technology and tested on a F7-RIL population derived from a Chiifu x Tetra cross. The design of the SNP assays were based on three considerations - the 50 bp flanking region of the SNPs should be strictly similar, the SNP should have a read-depth of ≥7 and no exon/intron junction should be present within the 101 bp target region. Using these criteria, a total of 640 markers (580 for genome-wide mapping and 60 for specific-region mapping) marking as many genes were tested for mapping. Out of 640 markers that were tested, 594 markers could be mapped unambiguously which included 542 markers for genome-wide mapping and 42 markers for fine mapping of the tet-o locus that is involved with the trait tetralocular ovary in the line Tetra.
A large number of SNPs and PSVs are present in the transcriptome of B. rapa lines for genome-wide linkage mapping and specific-region fine mapping. Criteria used for SNP identification delivered markers, more than 93% of which could be successfully mapped to the F7–RIL population of Chiifu x Tetra cross.
Brassica rapa; RNA-seq; Next generation sequencing; Single nucleotide polymorphism (SNP); Paralog specific variation (PSV); Coding DNA Sequences (CDS); KASPar assays
FANTOM DB, the database of Functional Annotation of RIKEN Mouse cDNA Clones, is designed to store sequence information of RIKEN full-length enriched mouse cDNA clones, graphical views of sequence analysis results, curated functional annotation information and additional descriptions, including Gene Ontology terms. RIKEN’s Mouse Gene Encyclopedia Project aims to collect full-length enriched cDNA clones from various mouse tissues, determine the full-length nucleotide sequences, infer their chromosomal locations by computer and characterize gene expression patterns. FANTOM DB has been developed to facilitate this work and to facilitate functional genomic studies such as positional candidate cloning, cDNA microarrays and protein interaction analyses. FANTOM DB contains 21 076 full-length cDNA sequences with rich functional annotations and is publicly available. FANTOM DB thus provides curated functional annotation to RIKEN full-length enriched mouse clones, and has links to other public resources. FANTOM DB can be accessed at http://fantom.gsc.riken.go.jp/db/.
The species Brassica rapa includes various vegetable crops. Production of these vegetable crops is usually impaired by heat stress. Some microRNAs (miRNAs) in Arabidopsis have been considered to mediate gene silencing in plant response to abiotic stress. However, it remains unknown whether or what miRNAs play a role in heat resistance of B. rapa. To identify genomewide conserved and novel miRNAs that are responsive to heat stress in B. rapa, we defined temperature thresholds of non-heading Chinese cabbage (B. rapa ssp. chinensis) and constructed small RNA libraries from the seedlings that had been exposed to high temperature (46 °C) for 1 h. By deep sequencing and data analysis, we selected a series of conserved and novel miRNAs that responded to heat stress. In total, Chinese cabbage shares at least 35 conserved miRNA families with Arabidopsis thaliana. Among them, five miRNA families were responsive to heat stress. Northern hybridization and real-time PCR showed that the conserved miRNAs bra-miR398a and bra-miR398b were heat-inhibitive and guided heat response of their target gene, BracCSD1; and bra-miR156h and bra-miR156g were heat-induced and its putative target BracSPL2 was down-regulated. According to the criteria of miRNA and miRNA* that form a duplex, 21 novel miRNAs belonging to 19 miRNA families were predicted. Of these, four were identified to be heat-responsive by Northern blotting and/or expression analysis of the putative targets. The two novel miRNAs bra-miR1885b.3 and bra-miR5718 negatively regulated their putative target genes. 5′-Rapid amplification of cDNA ends PCR indicated that three novel miRNAs cleaved the transcripts of their target genes where their precursors may have evolved from. These results broaden our perspective on the important role of miRNA in plant responses to heat.
Brassica rapa; heat response; miRNA; small RNA
Brassica species include both vegetable and oilseed crops, which are very important to the daily life of common human beings. Meanwhile, the Brassica species represent an excellent system for studying numerous aspects of plant biology, specifically for the analysis of genome evolution following polyploidy, so it is also very important for scientific research. Now, the genome of Brassica rapa has already been assembled, it is the time to do deep mining of the genome data.
BRAD, the Brassica database, is a web-based resource focusing on genome scale genetic and genomic data for important Brassica crops. BRAD was built based on the first whole genome sequence and on further data analysis of the Brassica A genome species, Brassica rapa (Chiifu-401-42). It provides datasets, such as the complete genome sequence of B. rapa, which was de novo assembled from Illumina GA II short reads and from BAC clone sequences, predicted genes and associated annotations, non coding RNAs, transposable elements (TE), B. rapa genes' orthologous to those in A. thaliana, as well as genetic markers and linkage maps. BRAD offers useful searching and data mining tools, including search across annotation datasets, search for syntenic or non-syntenic orthologs, and to search the flanking regions of a certain target, as well as the tools of BLAST and Gbrowse. BRAD allows users to enter almost any kind of information, such as a B. rapa or A. thaliana gene ID, physical position or genetic marker.
BRAD, a new database which focuses on the genetics and genomics of the Brassica plants has been developed, it aims at helping scientists and breeders to fully and efficiently use the information of genome data of Brassica plants. BRAD will be continuously updated and can be accessed through http://brassicadb.org.
Background and Aims
The polygalacturonase (PG) gene family has been found to be enriched in pollen of several species; however, little is currently known about the function of the PG gene in pollen development. To investigate the exact role that the PG gene has played in pollen development and about this family in general, one putative PG gene, Brassica campestris Male Fertility 9 (BcMF9), was isolated from Chinese cabbage (Brassica campestris ssp. chinensis, syn. B. rapa ssp. chinensis) and characterized.
RT-PCR, northern blotting and in situ hybridization were used to analyse the expression pattern of BcMF9, and antisense RNA technology was applied to study the function of this gene.
BcMF9 is expressed in particular in the tapetum and microspore during the late stages of pollen development. Antisense RNA transgenic plants that displayed decreased expression of BcMF9 showed pollen morphological defects that resulted in reduced pollen germination efficiency. Transmission electron microscopy revealed that the homogeneous pectic exintine layer of pollen facing the exterior was over-developed and predominantly occupied the intine, reversing the normal proportional distribution of the internal endintine layer and the external exintine in transgenic pollen. Inhibition of BcMF9 also resulted in break-up of the previously formed tectum and baculae from the beginning of the binucleate stage, as a result of premature degradation of tapetum.
Several lines of evidence, including patterns of BcMF9 expression and phenotypic defects, suggest a sporophytic role in exine patterning, and a gametophytic mode of action of BcMF9 in intine formation. BcMF9 might act as a co-ordinator in the late stages of tapetum degeneration, and subsequently in the regulation of wall material secretion and, in turn, exine formation. BcMF9 might also play a role in intine formation, possibly via regulation of the dynamic metabolism of pectin.
Brassica campestris; Chinese cabbage; exine; intine; PG; pollen wall; polygalacturonase; BcMF9
To identify genes associated with genic male sterility (GMS) that could be useful for hybrid breeding in Chinese cabbage (Brassicarapa ssp. pekinensis), floral bud transcriptome analysis was carried out using a B. rapa microarray with 300,000 probes (Br300K). Among 47,548 clones deposited on a Br300K microarray with seven probes of 60 nt length within the 3' 150 bp region, a total of 10,622 genes were differentially expressed between fertile and sterile floral buds; 4,774 and 5,848 genes were up-regulated over 2-fold in fertile and sterile buds, respectively. However, the expression of 1,413 and 199 genes showed fertile and sterile bud-specific features, respectively. Genes expressed specifically in fertile buds, possibly GMS-related genes, included homologs of several Arabidopsis male sterility-related genes, genes associated with the cell wall and synthesis of its surface proteins, pollen wall and coat components, signaling components, and nutrient supplies. However, most early genes for pollen development, genes for primexine and callose formation, and genes for pollen maturation and anther dehiscence showed no difference in expression between fertile and sterile buds. Some of the known genes associated with Arabidopsis pollen development showed similar expression patterns to those seen in this study, while others did not. BrbHLH89 and BrMYP99 are putative GMS genes. Additionally, 17 novel genes identified only in B. rapa were specifically and highly expressed only in fertile buds, implying the possible involvement in male fertility. All data suggest that Chinese cabbage GMS might be controlled by genes acting in post-meiotic tapetal development that are different from those known to be associated with Arabidopsis male sterility.
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.
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.
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.
The Brassica species, related to Arabidopsis thaliana, include an important group of crops and represent an excellent system for studying the evolutionary consequences of polyploidy. Previous studies have led to a proposed structure for an ancestral karyotype and models for the evolution of the B. rapa genome by triplication and segmental rearrangement, but these have not been validated at the sequence level.
We developed computational tools to analyse the public collection of B. rapa BAC end sequence, in order to identify candidates for representing collinearity discontinuities between the genomes of B. rapa and A. thaliana. For each putative discontinuity, one of the BACs was sequenced and analysed for collinearity with the genome of A. thaliana. Additional BAC clones were identified and sequenced as part of ongoing efforts to sequence four chromosomes of B. rapa. Strikingly few of the 19 inter-chromosomal rearrangements corresponded to the set of collinearity discontinuities anticipated on the basis of previous studies. Our analyses revealed numerous instances of newly detected collinearity blocks. For B. rapa linkage group A8, we were able to develop a model for the derivation of the chromosome from the ancestral karyotype. We were also able to identify a rearrangement event in the ancestor of B. rapa that was not shared with the ancestor of A. thaliana, and is represented in triplicate in the B. rapa genome. In addition to inter-chromosomal rearrangements, we identified and analysed 32 BACs containing the end points of segmental inversion events.
Our results show that previous studies of segmental collinearity between the A. thaliana, Brassica and ancestral karyotype genomes, although very useful, represent over-simplifications of their true relationships. The presence of numerous cryptic collinear genome segments and the frequent occurrence of segmental inversions mean that inference of the positions of genes in B. rapa based on the locations of orthologues in A. thaliana can be misleading. Our results will be of relevance to a wide range of plants that have polyploid genomes, many of which are being considered according to a paradigm of comprising conserved synteny blocks with respect to sequenced, related genomes.
In Chinese cabbage (Brassica rapa), the clubroot resistance (CR) gene CRb is effective against Plasmodiophora brassicae isolate No. 14, which is classified as pathotype group 3. Although markers linked to CRb have been reported, an accurate position in the genome and the gene structure are unknown. To determine the genomic location and estimate the structure of CRb, we developed 28 markers (average distance, 20.4 kb) around CRb and constructed a high-density partial map. The precise position of CRb was determined by using a population of 2,032 F2 plants generated by selfing B. rapa ‘CR Shinki.’ We determined that CRb is located in the 140-kb genomic region between markers KB59N07 and B1005 and found candidate resistance genes. Among other CR genes on chromosome R3, a genotype of CRa closest marker clearly matched those of CRb and Crr3 did not confer resistance to isolate No. 14. Based on the genotypes of 11 markers developed near CRb and resistance to isolate No. 14, 82 of 108 cultivars showed a strong correlation between genotypes and phenotypes. The results of this study will be useful for isolating CRb and breeding cultivars with resistance to pathotype group 3 by introducing CRb into susceptible cultivars through marker-assisted selection.
clubroot; Plasmodiophora brassicae; CRb; fine mapping; marker-assisted selection; Brassica rapa
Aliphatic glucosinolates are the predominant sulfur-rich plant secondary metabolites in economically important Brassica crops. Glucosinolates and their hydrolysis products are involved in plant–microbe, plant–insect, plant–animal, and plant–human interactions. It is, therefore, important to manipulate glucosinolate profiles and contents in Brassica species. In this study, aliphatic glucosinolates were genetically manipulated through homoeologous recombination in backcross lines followed by marker assisted selection in B. rapa. A resynthesized B. napus line, from a cross between B. rapa and B. oleracea, was backcrossed with Chinese cabbage doubled haploid line, RI16. Marker assisted selection for non-functional gene was performed in each backcross generations. Advanced backcross progenies (BC3F2) were developed to identify homoeologous gene replacement and/or introgression. Reduction in 5C aliphatic glucosinolates (gluconapoleiferin, glucoalyssin, and glucobrassicanapin) was observed in BC3F2 progenies of the recurrent parent that carried the GSL-ELONG- gene. The GSL-ELONG- positive backcross progenies were also screened by the A-genome and BraGSL-ELONG gene specific marker, which linked with 5C aliphatic glucosinolates. The A-genome specific marker was absent in the plants of advanced backcross progenies which showed reduction in 5C aliphatic glucosinolates. The results suggest that the functional allele had been replaced by the non-functional GSL-ELONG- allele from B. oleracea. Some advanced backcross progenies (BC3F2) positive for the GSL-ELONG- allele and the A-genome specific SCAR marker BraMAM1-1 did not show reduction in 5C aliphatic glucosinolates, suggesting that GSL-ELONG- allele is recessive. Replacement of the functional locus in the A-genome by non-functional counterpart in the C-genome reduced the content of 5C aliphatic glucosinolates in B. rapa seeds with 20 μmol/g.
Brassica rapa; Brassica oleracea; homeologous recombination; SCAR markers; MAS; glucosinolates
Radish (Raphanus sativus L.) belongs to Brassicaceae family and is a close relative of Brassica. This species shows a wide morphological diversity, and is an important vegetable especially in Asia. However, molecular research of radish is behind compared to that of Brassica. For example, reports on SSR (simple sequence repeat) markers are limited. Here, we designed 417 radish SSR markers from SSR-enriched genomic libraries and the cDNA data. Of the 256 SSR markers succeeded in PCR, 130 showed clear polymorphisms between two radish lines; a rat-tail radish and a Japanese cultivar, ‘Harufuku’. As a test case for evaluation of the present SSRs, we conducted two studies. First, we selected 16 SSRs to calculate polymorphism information contents (PICs) using 16 radish cultivars and four other Brassicaceae species. These markers detected 3–15 alleles (average = 9.6). PIC values ranged from 0.54 to 0.92 (average = 0.78). Second, part of the present SSRs were tested for mapping using our previously-examined mapping population. The map spanned 672.7 cM with nine linkage groups (LGs). The 21 radish SSR markers were distributed throughout the LGs. The SSR markers developed here would be informative and useful for genetic analysis in radish and its related species.
Radish; Raphanus sativus; SSR (simple sequence repeat) marker
The Multinational Brassica rapa Genome Sequencing Project (BrGSP) has developed valuable genomic resources, including BAC libraries, BAC-end sequences, genetic and physical maps, and seed BAC sequences for Brassica rapa. An integrated linkage map between the amphidiploid B. napus and diploid B. rapa will facilitate the rapid transfer of these valuable resources from B. rapa to B. napus (Oilseed rape, Canola).
In this study, we identified over 23,000 simple sequence repeats (SSRs) from 536 sequenced BACs. 890 SSR markers (designated as BrGMS) were developed and used for the construction of an integrated linkage map for the A genome in B. rapa and B. napus. Two hundred and nineteen BrGMS markers were integrated to an existing B. napus linkage map (BnaNZDH). Among these mapped BrGMS markers, 168 were only distributed on the A genome linkage groups (LGs), 18 distrubuted both on the A and C genome LGs, and 33 only distributed on the C genome LGs. Most of the A genome LGs in B. napus were collinear with the homoeologous LGs in B. rapa, although minor inversions or rearrangements occurred on A2 and A9. The mapping of these BAC-specific SSR markers enabled assignment of 161 sequenced B. rapa BACs, as well as the associated BAC contigs to the A genome LGs of B. napus.
The genetic mapping of SSR markers derived from sequenced BACs in B. rapa enabled direct links to be established between the B. napus linkage map and a B. rapa physical map, and thus the assignment of B. rapa BACs and the associated BAC contigs to the B. napus linkage map. This integrated genetic linkage map will facilitate exploitation of the B. rapa annotated genomic resources for gene tagging and map-based cloning in B. napus, and for comparative analysis of the A genome within Brassica species.
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 Bio∷DB∷GFF methods,
facilitating a wide range of data-mining possibilities.
The performance of natural enemies, such as parasitoid wasps, is affected by differences in the quality of the host’s diet, frequently mediated by species or population-related differences in plant allelochemistry. Here, we compared survival, development time, and body mass in a generalist herbivore, the cabbage moth, Mamestra brassicae, and its solitary endoparasitoid, Microplitis mediator, when reared on two cultivated (CYR and STH) and three wild (KIM, OH, and WIN) populations of cabbage, Brassica oleracea. Plants either were undamaged or induced by feeding of larvae of the cabbage butterfly, Pieris rapae. Development and biomass of M. brassicae and Mi. mediator were similar on both cultivated and one wild cabbage population (KIM), intermediate on the OH population, and significantly lower on the WIN population. Moreover, development was prolonged and biomass was reduced on herbivore-induced plants. However, only the survival of parasitized hosts (and not that of healthy larvae) was affected by induction. Analysis of glucosinolates in leaves of the cabbages revealed higher levels in the wild populations than cultivars, with the highest concentrations in WIN plants. Multivariate statistics revealed a negative correlation between insect performance and total levels of glucosinolates (GS) and levels of 3-butenyl GS. However, GS chemistry could not explain the reduced performance on induced plants since only indole GS concentrations increased in response to herbivory, which did not affect insect performance based on multivariate statistics. This result suggests that, in addition to aliphatic GS, other non-GS chemicals are responsible for the decline in insect performance, and that these chemicals affect the parasitoid more strongly than the host. Remarkably, when developing on WIN plants, the survival of Mi. mediator to adult eclosion was much higher than in its host, M. brassicae. This may be due to the fact that hosts parasitized by Mi. mediator pass through fewer instars, and host growth is arrested when they are only a fraction of the size of healthy caterpillars. Certain aspects of the biology and life-history of the host and parasitoid may determine their response to chemical challenges imposed by the food plant.
Glucosinolates; Parasitoid; Induction; Cabbage; Brassica oleracea; Microplitis mediator; Mamestra brassicae; Koinobiont
Insect hosts have evolved immunity against invasion by parasitoids, and in co-evolutionary response parasitoids have also developed strategies to overcome host immune systems. The mechanisms through which parasitoid venoms disrupt the promotion of host immunity are still unclear. We report here a new mechanism evolved by parasitoid Pteromalus puparum, whose venom inhibited the promotion of immunity in its host Pieris rapae (cabbage white butterfly).
A full-length cDNA encoding a C-type lectin (Pr-CTL) was isolated from P. rapae. Quantitative PCR and immunoblotting showed that injection of bacterial and inert beads induced expression of Pr-CTL, with peaks of mRNA and Pr-CTL protein levels at 4 and 8 h post beads challenge, respectively. In contrast, parasitoid venom suppressed Pr-CTL expression when co-injected with beads, in a time and dose-dependent manner. Immunolocalization and immunoblotting results showed that Pr-CTL was first detectable in vesicles present in cytoplasm of granulocytes in host hemolymph, and was then secreted from cells into circulatory fluid. Finally, the secreted Pr-CTL bound to cellular membranes of both granulocytes and plasmatocytes. Injection of double-stranded RNA specific for target gene decreased expression of Pr-CTL, and a few other host immune-related genes. Suppression of Pr-CTL expression also down-regulated antimicrobial and phenoloxidase activities, and reducing phagocytotic and encapsulation rates in host. The inhibitory effect of parasitoid venom on host encapsulation is consistent with its effect in suppressing Pr-CTL expression. Binding assay results showed that recombinant Pr-CTL directly attached to the surface of P. puparum egges. We infer that Pr-CTL may serve as an immune signalling co-effector, first binding to parasitoid eggs, regulating expression of a set of immune-related genes and promoting host immunity.
P. puparum venom inhibits promotion of host immune responses by silencing expression of host C-type lectin gene Pr-CTL, whose expression affected transcription of other host immune-related genes.
In 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.
We 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.
We 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.
Cytokinins (CKs) have significant roles in various aspects of plant growth and development, and they are also involved in plant stress adaptations. The fine-tuning of the controlled CK levels in individual tissues, cells, and organelles is properly maintained by isopentenyl transferases (IPTs) and cytokinin oxidase/dehydrogenases (CKXs). Chinese cabbage is one of the most economically important vegetable crops worldwide. The whole genome sequencing of Brassica rapa enables us to perform the genome-wide identification and functional analysis of the IPT and CKX gene families.
In this study, a total of 13 BrIPT genes and 12 BrCKX genes were identified. The gene structures, conserved domains and phylogenetic relationships were analyzed. The isoelectric point, subcellular localization and glycosylation sites of the proteins were predicted. Segmental duplicates were found in both BrIPT and BrCKX gene families. We also analyzed evolutionary patterns and divergence of the IPT and CKX genes in the Cruciferae family. The transcription levels of BrIPT and BrCKX genes were analyzed to obtain an initial picture of the functions of these genes. Abiotic stress elements related to adverse environmental stimuli were found in the promoter regions of BrIPT and BrCKX genes and they were confirmed to respond to drought and high salinity conditions. The effects of 6-BA and ABA on the expressions of BrIPT and BrCKX genes were also investigated.
The expansion of BrIPT and BrCKX genes after speciation from Arabidopsis thaliana is mainly attributed to segmental duplication events during the whole genome triplication (WGT) and substantial duplicated genes are lost during the long evolutionary history. Genes produced by segmental duplication events have changed their expression patterns or may adopted new functions and thus are obtained. BrIPT and BrCKX genes respond well to drought and high salinity stresses, and their transcripts are affected by exogenous hormones, such as 6-BA and ABA, suggesting their potential roles in abiotic stress conditions and regulatory mechanisms of plant hormone homeostasis. The appropriate modulation of endogenous CKs levels by IPT and CKX genes is a promising approach for developing economically important high-yielding and high-quality stress-tolerant crops in agriculture.
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;
Following successful completion of the Brassica rapa sequencing project, the next step is to investigate functions of individual genes/proteins. For Arabidopsis thaliana, large amounts of protein–protein interaction (PPI) data are available from the major PPI databases (DBs). It is known that Brassica crop species are closely related to A. thaliana. This provides an opportunity to infer the B. rapa interactome using PPI data available from A. thaliana. In this paper, we present an inferred B. rapa interactome that is based on the A. thaliana PPI data from two resources: (i) A. thaliana PPI data from three major DBs, BioGRID, IntAct, and TAIR. (ii) ortholog-based A. thaliana PPI predictions. Linking between B. rapa and A. thaliana was accomplished in three complementary ways: (i) ortholog predictions, (ii) identification of gene duplication based on synteny and collinearity, and (iii) BLAST sequence similarity search. A complementary approach was also applied, which used known/predicted domain–domain interaction data. Specifically, since the two species are closely related, we used PPI data from A. thaliana to predict interacting domains that might be conserved between the two species. The predicted interactome was investigated for the component that contains known A. thaliana meiotic proteins to demonstrate its usability.
Brassica rapa; Arabidopsis thaliana; interactome; protein–protein interaction; domain–domain interaction; meiosis