Anthocyanin biosynthetic genes in Brassica rapa

Background

Anthocyanins are a group of flavonoid compounds. As a group of important secondary metabolites, they perform several key biological functions in plants. Anthocyanins also play beneficial health roles as potentially protective factors against cancer and heart disease. To elucidate the anthocyanin biosynthetic pathway in Brassica rapa, we conducted comparative genomic analyses between Arabidopsis thaliana and B. rapa on a genome-wide level.

Results

In total, we identified 73 genes in B. rapa as orthologs of 41 anthocyanin biosynthetic genes in A. thaliana. In B. rapa, the anthocyanin biosynthetic genes (ABGs) have expanded and most genes exist in more than one copy. The anthocyanin biosynthetic structural genes have expanded through whole genome and tandem duplication in B. rapa. More structural genes located upstream of the anthocyanin biosynthetic pathway have been retained than downstream. More negative regulatory genes are retained in the anthocyanin biosynthesis regulatory system of B. rapa.

Conclusions

These results will promote an understanding of the genetic mechanism of anthocyanin biosynthesis, as well as help the improvement of the nutritional quality of B. rapa through the breeding of high anthocyanin content varieties.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-426) contains supplementary material, which is available to authorized users.     

Carotenoid biosynthetic genes in Brassica rapa: comparative genomic analysis,phylogenetic analysis,and expression profiling

Background

Carotenoids are isoprenoid compounds synthesized by all photosynthetic organisms. Despite much research on carotenoid biosynthesis in the model plant Arabidopsis thaliana, there is a lack of information on the carotenoid pathway in Brassica rapa. To better understand its carotenoid biosynthetic pathway, we performed a systematic analysis of carotenoid biosynthetic genes at the genome level in B. rapa.

Results

We identified 67 carotenoid biosynthetic genes in B. rapa, which were orthologs of the 47 carotenoid genes in A. thaliana. A high level of synteny was observed for carotenoid biosynthetic genes between A. thaliana and B. rapa. Out of 47 carotenoid biosynthetic genes in A. thaliana, 46 were successfully mapped to the 10 B. rapa chromosomes, and most of the genes retained more than one copy in B. rapa. The gene expansion was caused by the whole-genome triplication (WGT) event experienced by Brassica species. An expression analysis of the carotenoid biosynthetic genes suggested that their expression levels differed in root, stem, leaf, flower, callus, and silique tissues. Additionally, the paralogs of each carotenoid biosynthetic gene, which were generated from the WGT in B. rapa, showed significantly different expression levels among tissues, suggesting differentiated functions for these multi-copy genes in the carotenoid pathway.

Conclusions

This first systematic study of carotenoid biosynthetic genes in B. rapa provides insights into the carotenoid metabolic mechanisms of Brassica crops. In addition, a better understanding of carotenoid biosynthetic genes in B. rapa will contribute to the development of conventional and transgenic B. rapa cultivars with enriched carotenoid levels in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1655-5) contains supplementary material, which is available to authorized users.     

Functional innovations of three chronological mesohexaploid Brassica rapa genomes

Background

The Brassicaceae family is an exemplary model for studying plant polyploidy. The Brassicaceae knowledge-base includes the well-annotated Arabidopsis thaliana reference sequence; well-established evidence for three rounds of whole genome duplication (WGD); and the conservation of genomic structure, with 24 conserved genomic blocks (GBs). The recently released Brassica rapa draft genome provides an ideal opportunity to update our knowledge of the conserved genomic structures in Brassica, and to study evolutionary innovations of the mesohexaploid plant, B. rapa.

Results

Three chronological B. rapa genomes (recent, young, and old) were reconstructed with sequence divergences, revealing a trace of recursive WGD events. A total of 636 fast evolving genes were unevenly distributed throughout the recent and young genomes. The representative Gene Ontology (GO) terms for these genes were ‘stress response’ and ‘development’ both through a change in protein modification or signaling, rather than by enhancing signal recognition. In retention patterns analysis, 98% of B. rapa genes were retained as collinear gene pairs; 77% of those were singly-retained in recent or young genomes resulting from death of the ancestral copies, while others were multi-retained as long retention genes. GO enrichments indicated that single retention genes mainly function in the interpretation of genetic information, whereas, multi-retention genes were biased toward signal response, especially regarding development and defense. In the recent genome, 13,302, 5,790, and 20 gene pairs were multi-retained following Brassica whole genome triplication (WGT) events with 2, 3, and 4 homoeologous copies, respectively. Enriched GO-slim terms from B. rapa homomoelogues imply that a major effect of the B. rapa WGT may have been to acquire environmental adaptability or to change the course of development. These homoeologues seem to more frequently undergo subfunctionalization with spatial expression patterns compared with other possible events including nonfunctionalization and neofunctionalization.

Conclusion

We refined Brassicaceae GB information using the latest genomic resources, and distinguished three chronologically ordered B. rapa genomes. B. rapa genes were categorized into fast evolving, single- and multi-retention genes, and long retention genes by their substitution rates and retention patterns. Representative functions of the categorized genes were elucidated, providing better understanding of B. rapa evolution and the Brassica genus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-606) contains supplementary material, which is available to authorized users.     

Genic Microsatellite Markers in Brassica rapa: Development, Characterization, Mapping, and Their Utility in Other Cultivated and Wild Brassica Relatives

Genic microsatellite markers, also known as functional markers, are preferred over anonymous markers as they reveal the variation in transcribed genes among individuals. In this study, we developed a total of 707 expressed sequence tag-derived simple sequence repeat markers (EST-SSRs) and used for development of a high-density integrated map using four individual mapping populations of B. rapa. This map contains a total of 1426 markers, consisting of 306 EST-SSRs, 153 intron polymorphic markers, 395 bacterial artificial chromosome-derived SSRs (BAC-SSRs), and 572 public SSRs and other markers covering a total distance of 1245.9 cM of the B. rapa genome. Analysis of allelic diversity in 24 B. rapa germplasm using 234 mapped EST-SSR markers showed amplification of 2 alleles by majority of EST-SSRs, although amplification of alleles ranging from 2 to 8 was found. Transferability analysis of 167 EST-SSRs in 35 species belonging to cultivated and wild brassica relatives showed 42.51% (Sysimprium leteum) to 100% (B. carinata, B. juncea, and B. napus) amplification. Our newly developed EST-SSRs and high-density linkage map based on highly transferable genic markers would facilitate the molecular mapping of quantitative trait loci and the positional cloning of specific genes, in addition to marker-assisted selection and comparative genomic studies of B. rapa with other related species.     

Tobacco TTG2 and ARF8 function concomitantly to control flower colouring by regulating anthocyanin synthesis genes

• Recently we elucidated that tobacco TTG2 cooperates with ARF8 to regulate the vegetative growth and seed production.
• Here we show that TTG2 and ARF8 control flower colouring by regulating expression of ANS and DFR genes, which function in anthocyanin biosynthesis.
• Genetic modifications that substantially altered expression levels of the TTG2 gene and production quantities of TTG2 protein were correlated with flower development and colouring. Degrees of flower colour were increased by TTG2 overexpression but decreased through TTG2 silencing, in coincidence with high and low concentrations of anthocyanins in flowers. Of five genes involved in the anthocyanin biosynthesis pathway, only ANS and DFR were TTG2‐regulated and displayed enhancement and diminution of expression with TTG2 overexpression and silencing, respectively. The floral expression of ANS and DFR also needed a functional ARF8 gene, as ANS and DFR expression were attenuated by ARF8 silencing, which concomitantly diminished the role of TTG2 in anthocyanin production. While ARF8 required TTG2 to be expressed by itself and to regulate ANS and DFR expression, the concurrent presence of normally functional TTG2 and ARF8 was critical for floral production of anthocyanins and also for flower colouration.
• Our data suggest that TTG2 functions concomitantly with ARF8 to control degrees of flower colour by regulating expression of ANS and DFR, which are involved in the anthocyanin biosynthesis pathway. ARF8 depends on TTG2 to regulate floral expression of ANS and DFR with positive effects on anthocyanin production and flower colour.

     

Glucosinolate profiling of Brassica rapa cultivars after infection by Leptosphaeria maculans and Fusarium oxysporum

The glucosinolate contents of two different cultivars of Brassica rapa (Herfstraap and Oleifera) infected with Leptosphaeria maculans and Fusarium oxysporum were determined. Infection triggered the accumulation of aliphatic glucosinolates (gluconapin, progoitrin, glucobrassicanapin and gluconapoleiferin) and indole glucosinolate (4-hydroxy-glucobrassicin) in Herfstraap and of two indole glucosinolates (glucobrassicin and 4-hydroxy-glucobrassicin) in Oleifera. While total and aliphatic glucosinolates decreased significantly in Oleifera, a large increase was observed in Herfstraap after fungal infection. The indole glucosinolate glucobrassicin accumulated in Oleifera at a higher rate than Herfstraap especially after infection with F. oxysporum. Apparently the interaction between fungus and B. rapa is cultivar and fungal species specific.     

Flower colour and cytochromes P450

Cytochromes P450 play important roles in biosynthesis of flavonoids and their coloured class of compounds, anthocyanins, both of which are major floral pigments. The number of hydroxyl groups on the B-ring of anthocyanidins (the chromophores and precursors of anthocyanins) impact the anthocyanin colour, the more the bluer. The hydroxylation pattern is determined by two cytochromes P450, flavonoid 3′-hydroxylase (F3′H) and flavonoid 3′,5′-hydroxylase (F3′5′H) and thus they play a crucial role in the determination of flower colour. F3′H and F3′5′H mostly belong to CYP75B and CYP75A, respectively, except for the F3′5′Hs in Compositae that were derived from gene duplication of CYP75B and neofunctionalization. Roses and carnations lack blue/violet flower colours owing to the deficiency of F3′5′H and therefore lack the B-ring-trihydroxylated anthocyanins based upon delphinidin. Successful redirection of the anthocyanin biosynthesis pathway to delphinidin was achieved by expressing F3′5′H coding regions resulting in carnations and roses with novel blue hues that have been commercialized. Suppression of F3′5′H and F3′H in delphinidin-producing plants reduced the number of hydroxyl groups on the anthocyanidin B-ring resulting in the production of monohydroxylated anthocyanins based on pelargonidin with a shift in flower colour to orange/red. Pelargonidin biosynthesis is enhanced by additional expression of a dihydroflavonol 4-reductase that can use the monohydroxylated dihydrokaempferol (the pelargonidin precursor). Flavone synthase II (FNSII)-catalysing flavone biosynthesis from flavanones is also a P450 (CYP93B) and contributes to flower colour, because flavones act as co-pigments to anthocyanins and can cause blueing and darkening of colour. However, transgenic plants expression of a FNSII gene yielded paler flowers owing to a reduction of anthocyanins because flavanones are precursors of anthocyanins and flavones.     

Development of Full-Length cDNAs from Chinese Cabbage (Brassica rapa Subsp. pekinensis) and Identification of Marker Genes for Defence Response

Arabidopsis belongs to the Brassicaceae family and plays an important role as a model plant for which researchers have developed fine-tuned genome resources. Genome sequencing projects have been initiated for other members of the Brassicaceae family. Among these projects, research on Chinese cabbage (Brassica rapa subsp. pekinensis) started early because of strong interest in this species. Here, we report the development of a library of Chinese cabbage full-length cDNA clones, the RIKEN BRC B. rapa full-length cDNA (BBRAF) resource, to accelerate research on Brassica species. We sequenced 10 000 BBRAF clones and confirmed 5476 independent clones. Most of these cDNAs showed high homology to Arabidopsis genes, but we also obtained more than 200 cDNA clones that lacked any sequence homology to Arabidopsis genes. We also successfully identified several possible candidate marker genes for plant defence responses from our analysis of the expression of the Brassica counterparts of Arabidopsis marker genes in response to salicylic acid and jasmonic acid. We compared gene expression of these markers in several Chinese cabbage cultivars. Our BBRAF cDNA resource will be publicly available from the RIKEN Bioresource Center and will help researchers to transfer Arabidopsis-related knowledge to Brassica crops.     

A triallelic genetic male sterility locus in Brassica napus: an integrative strategy for its physical mapping and possible local chromosome evolution around it

Background and Aims

Spontaneous male sterility is an advantageous trait for both constructing efficient pollination control systems and for understanding the developmental process of the male reproductive unit in many crops. A triallelic genetic male-sterile locus (BnMs5) has been identified in Brassica napus; however, its complicated genome structure has greatly hampered the isolation of this locus. The aim of this study was to physically map BnMs5 through an integrated map-based cloning strategy and analyse the local chromosomal evolution around BnMs5.

Methods

A large F₂ population was used to integrate the existing genetic maps around BnMs5. A map-based cloning strategy in combination with comparative mapping among B. napus, Arabidopsis, Brassica rapa and Brassica oleracea was employed to facilitate the identification of a target bacterial artificial chromosome (BAC) clone covering the BnMs5 locus. The genomic sequences from the Brassica species were analysed to reveal the regional chromosomal evolution around BnMs5.

Key Results

BnMs5 was finally delimited to a 0·3-cM genetic fragment from an integrated local genetic map, and was anchored on the B. napus A8 chromosome. Screening of a B. napus BAC clone library and identification of the positive clones validated that JBnB034L06 was the target BAC clone. The closest flanking markers restrict BnMs5 to a 21-kb region on JBnB034L06 containing six predicted functional genes. Good collinearity relationship around BnMs5 between several Brassica species was observed, while violent chromosomal evolutionary events including insertions/deletions, duplications and single nucleotide mutations were also found to have extensively occurred during their divergence.

Conclusions

This work represents major progress towards the molecular cloning of BnMs5, as well as presenting a powerful, integrative method to mapping loci in plants with complex genomic architecture, such as the amphidiploid B. napus.     

Comparison of complete mitochondrial genomes of pine wilt nematode Bursaphelenchus xylophilus and Bursaphelenchus mucronatus (Nematoda: Aphelenchoidea) and development of a molecular tool for species identification

We determined the complete mitochondrial genome sequences for Bursaphelenchus mucronatus, one species of pinewood nematode. The genome is a circular-DNA molecule of 14,583 bp (195 bp smaller than its congener Bursaphelenchus xylophilus) and contains 12 protein-coding genes (lacking atp8), 22 tRNA genes, and 2 rRNA genes encoded in the same direction, consistent with most other nematodes. Based on sequence comparison of mtDNA genomes, we developed a PCR-based molecular assay to differentiate B. xylophilus (highly pathogenic) and B. mucronatus (relatively less virulent) using species-specific primers. The molecular identification system employs multiplex-PCR and is very effective and reliable for discriminating these Bursaphelenchus species, which are economically important, but difficult to distinguish based on morphology. The comparison of the mitochondrial genomes and molecular identification system of the two species of Bursaphelenchus spp. should provide a rich source of genetic information to support the effective control and management (quarantine) of the pine wilt disease caused by pinewood nematodes.