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.     

Identification and characterization of LIM gene family in Brassica rapa

Background

LIM (Lin-11, Isl-1 and Mec-3 domains) genes have been reported to trigger the formation of actin bundles, a major higher-order cytoskeletal assembly, in higher plants; however, the stress resistance related functions of these genes are still not well known. In this study, we collected 22 LIM genes designated as Brassica rapa LIM (BrLIM) from the Brassica database, analyzed the sequences, compared them with LIM genes of other plants and analyzed their expression after applying biotic and abiotic stresses in Chinese cabbage.

Results

Upon sequence analysis these genes were confirmed as LIM genes and found to have a high degree of homology with LIM genes of other species. These genes showed distinct clusters when compared to other recognized LIM proteins upon phylogenetic analysis. Additionally, organ specific expression of these genes was observed in Chinese cabbage plants, with BrPLIM2a, b, c, BrDAR1, BrPLIM2e, f and g only being expressed in flower buds. Furthermore, the expression of these genes (except for BrDAR1 and BrPLIM2e) was high in the early flowering stages. The remaining genes were expressed in almost all organs tested. All BrDAR genes showed higher expression in flower buds compared to other organs. These organ specific expressions were clearly correlated with the phylogenetic grouping. In addition, BrWLIM2c and BrDAR4 responded to Fusarium oxysporum f. sp. conglutinans infection, while commonly two BrDARs and eight BrLIMs responded to cold, ABA and pH (pH5, pH7 and pH9) stress treatments in Chinese cabbage plants.

Conclusion

Taken together, the results of this study indicate that BrLIM and BrDAR genes may be involved in resistance against biotic and abiotic stresses in Brassica.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-641) 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.     

The quantitative-genetic and QTL architecture of trait integration and modularity in Brassica rapa across simulated seasonal settings

Within organisms, groups of traits with different functions are frequently modular, such that variation among modules is independent and variation within modules is tightly integrated, or correlated. Here, we investigated patterns of trait integration and modularity in Brassica rapa in response to three simulated seasonal temperature/photoperiod conditions. The goals of this research were to use trait correlations to understand patterns of trait integration and modularity within and among floral, vegetative and phenological traits of B. rapa in each of three treatments, to examine the QTL architecture underlying patterns of trait integration and modularity, and to quantify how variation in temperature and photoperiod affects the correlation structure and QTL architecture of traits. All floral organs of B. rapa were strongly correlated, and contrary to expectations, floral and vegetative traits were also correlated. Extensive QTL co-localization suggests that covariation of these traits is likely due to pleiotropy, although physically linked loci that independently affect individual traits cannot be ruled out. Across treatments, the structure of genotypic and QTL correlations was generally conserved. Any observed variation in genetic architecture arose from genotype × environment interactions (GEIs) and attendant QTL × E in response to temperature but not photoperiod.     

Time-lapse imaging of self- and cross-pollinations in Brassica rapa

Background and Aims

Pollination is an important process in the life cycle of plants and is the first step in bringing together the male and female gametophytes for plant reproduction. While pollination has been studied for many years, accurate knowledge of the morphological aspects of this process is still far from complete. This study therefore focuses on a morphological characterization of pollination, using time-series image analysis of self- and cross-pollinations in Brassica rapa.

Methods

Time-lapse imaging of pollen behaviour during self- and cross-pollinations was recorded for 90 min, at 1 min intervals, using a stereoscopic microscope. Using time-series digital images of pollination, characteristic features of pollen behaviours during self- and cross-pollinations were studied.

Key Results

Pollen exhibited various behaviours in both self- and cross-pollinations, and these were classified into six representative patterns: germination, expansion, contraction, sudden contraction, pulsation and no change. It is noteworthy that in ‘contraction’ pollen grains shrunk within a short period of 30–50 min, and in ‘pulsation’ repeated expansion and contraction occurred with an interval of 10 min, suggesting that a dehydration system is operating in pollination. All of the six patterns were observed on an individual stigma with both self- and cross-pollinations, and the difference between self- and cross-pollinations was in the ratios of the different behaviours. With regard to water transport to and from pollen grains, this occurred in multiple steps, before, during and after hydration. Thus, pollination is regulated by a combination of multiple components of hydration, rehydration and dehydration systems.

Conclusions

Regulated hydration of pollen is a key process for both pollination and self-incompatibility, and this is achieved by a balanced complex of hydration, dehydration and nutrient supply to pollen grains from stigmatic papilla cells.     

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.     

Seed predation,pathogen infection and life-history traits in Brassica rapa

Herbivory and disease can shape the evolution of plant populations, but their joint effects are rarely investigated. Families of plants of Brassica rapa (Brassicaceae) were grown from seeds collected in two naturalized populations in an experimental garden. We examined leaf infection by the fungus Alternaria, seed predation by a gall midge (Cecidomyiidae) and plant life-history traits. Plants from one population had heavier seeds, were more likely to flower, had less fungal infection, had more seed predation and were more fecund. Fungal infection score and seed predation rate increased with plant size, but large plants still had the greatest number of undamaged fruits. Spatial heterogeneity in the experimental garden was significant; seed predation rate and fecundity varied among blocks. An apparent tradeoff existed between susceptibility to disease and seed predation: plants with the highest fungal infection score had the lowest seed predation rate. Alternaria infection varied between populations, but the disease had no effect on fecundity. Seed predation did reduce fecundity. Damaged fruits had 31.4% fewer intact seeds. However, evidence for additive genetic variation in resistance to seed predation was weak. Therefore, neither disease nor seed predation was likely to be a strong agent of genetically based fecundity selection.     

Mapping of genes affecting linolenic acid content in Brassica rapa ssp. oleifera

F₂ progeny segregating for linolenic acid content were used to identify genes and develop markers for linolenic acid in spring turnip rape (Brassica rapa ssp. oleifera). A candidate gene approach applying the rapeseed fad3 gene and bulked segregant analysis with RAPD markers was used. A total of 27 markers were distributed in three linkage groups which each exhibited a QTL for linolenic acid. Jointly the three QTLs accounted for 73.5% of the variation in linolenic acid level in this population. The fad3 gene was mapped near one QTL controlling 23.5% of the variation. Allele-specific markers were developed for fad3 and can be used for marker-assisted selection in future spring turnip rape breeding programmes.     

Distribution of calcium (Ca) and magnesium (Mg) in the leaves of Brassica rapa under varying exogenous Ca and Mg supply

Background and Aims

Leafy vegetable Brassica crops are an important source of dietary calcium (Ca) and magnesium (Mg) and represent potential targets for increasing leaf Ca and Mg concentrations through agronomy or breeding. Although the internal distribution of Ca and Mg within leaves affects the accumulation of these elements, such data are not available for Brassica. The aim of this study was to characterize the internal distribution of Ca and Mg in the leaves of a vegetable Brassica and to determine the effects of altered exogenous Ca and Mg supply on this distribution.

Methods

Brassica rapa ssp. trilocularis ‘R-o-18’ was grown at four different Ca:Mg treatments for 21 d in a controlled environment. Concentrations of Ca and Mg were determined in fully expanded leaves using inductively coupled plasma-mass spectrometry (ICP-MS). Internal distributions of Ca and Mg were determined in transverse leaf sections at the base and apex of leaves using energy-dispersive X-ray spectroscopy (EDS) with cryo-scanning electron microscopy (cryo-SEM).

Key Results

Leaf Ca and Mg concentrations were greatest in palisade and spongy mesophyll cells, respectively, although this was dependent on exogenous supply. Calcium accumulation in palisade mesophyll cells was enhanced slightly under high Mg supply; in contrast, Mg accumulation in spongy mesophyll cells was not affected by Ca supply.

Conclusions

The results are consistent with Arabidopsis thaliana and other Brassicaceae, providing phenotypic evidence that conserved mechanisms regulate leaf Ca and Mg distribution at a cellular scale. The future study of Arabidopsis gene orthologues in mutants of this reference B. rapa genotype will improve our understanding of Ca and Mg homeostasis in plants and may provide a model-to-crop translation pathway for targeted breeding.     

Mapping of a QTL for oleic acid concentration in spring turnip rape (Brassica rapa ssp. oleifera)

Bulk segregant analysis was used to search for RAPD (random amplified polymorphic DNA) markers linked to gene(s) affecting oleic acid concentration in an F₂ population from the Brassica rapa ssp. oleifera cross Jo4002 x a high oleic acid individual from line Jo4072. Eight primers (=8 markers) out of 104 discriminated the high and low bulks consisting of extreme individuals from the oleic acid distribution. These markers were analysed throughout the entire F₂ population, and their association with oleic acid was studied using both interval mapping and ANOVA analysis. Six of the markers mapped to one linkage group. A quantitative trait locus (QTL) affecting oleic acid concentration was found to reside within this linkage group with a LOD score >15. The most suitable marker for oleic acid content is OPH-17, a codominant marker close (<4cM) to the QTL. The mean seed oleic acid content in the F₂ individuals carrying the larger allele of this marker was 80.14±9.76%; in individuals with the smaller allele, 54.53±6.83%; in the heterozygotes, 65.47±8.15%. To increase reproducibility, the RAPD marker was converted into a SCAR (sequence characterized amplied region) marker with specific primers. Marker OPH-17 can be used to select spring turnip rape individuals with the desired oleic acid content.     

Seed bank dynamics govern persistence of Brassica hybrids in crop and natural habitats

Background and Aims Gene flow from crops to their wild relatives has the potential to alter population growth rates and demography of hybrid populations, especially when a new crop has been genetically modified (GM). This study introduces a comprehensive approach to assess this potential for altered population fitness, and uses a combination of demographic data in two habitat types and mathematical (matrix) models that include crop rotations and outcrossing between parental species.Methods Full life-cycle demographic rates, including seed bank survival, of non-GM Brassica rapa × B. napus F₁ hybrids and their parent species were estimated from experiments in both agricultural and semi-natural habitats. Altered fitness potential was modelled using periodic matrices including crop rotations and outcrossing between parent species.Key Results The demographic vital rates (i.e. for major stage transitions) of the hybrid population were intermediate between or lower than both parental species. The population growth rate (λ) of hybrids indicated decreases in both habitat types, and in a semi-natural habitat hybrids became extinct at two sites. Elasticity analyses indicated that seed bank survival was the greatest contributor to λ. In agricultural habitats, hybrid populations were projected to decline, but with persistence times up to 20 years. The seed bank survival rate was the main driver determining persistence. It was found that λ of the hybrids was largely determined by parental seed bank survival and subsequent replenishment of the hybrid population through outcrossing of B. rapa with B. napus.Conclusions Hybrid persistence was found to be highly dependent on the seed bank, suggesting that targeting hybrid seed survival could be an important management option in controlling hybrid persistence. For local risk mitigation, an increased focus on the wild parent is suggested. Management actions, such as control of B. rapa, could indirectly reduce hybrid populations by blocking hybrid replenishment.     

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.     

Mapping and analysis of a novel candidate Fusarium wilt resistance gene FOC1 in Brassica oleracea

Background

Cabbage Fusarium wilt is a major disease worldwide that can cause severe yield loss in cabbage (Brassica olerecea). Although markers linked to the resistance gene FOC1 have been identified, no candidate gene for it has been determined so far. In this study, we report the fine mapping and analysis of a candidate gene for FOC1 using a double haploid (DH) population with 160 lines and a F₂ population of 4000 individuals derived from the same parental lines.

Results

We confirmed that the resistance to Fusarium wilt was controlled by a single dominant gene based on the resistance segregation ratio of the two populations. Using InDel primers designed from whole-genome re-sequencing data for the two parental lines (the resistant inbred-line 99–77 and the highly susceptible line 99–91) and the DH population, we mapped the resistance gene to a 382-kb genomic region on chromosome C06. Using the F₂ population, we narrowed the region to an 84-kb interval that harbored ten genes, including four probable resistance genes (R genes): Bol037156, Bol037157, Bol037158 and Bol037161 according to the gene annotations from BRAD, the genomic database for B. oleracea. After correcting the model of the these genes, we re-predicted two R genes in the target region: re-Bol037156 and re-Bol0371578. The latter was excluded after we compared the two genes’ sequences between ten resistant materials and ten susceptible materials. For re-Bol037156, we found high identity among the sequences of the resistant lines, while among the susceptible lines, there were two types of InDels (a 1-bp insertion and a 10-bp deletion), each of which caused a frameshift and terminating mutation in the cDNA sequences. Further sequence analysis of the two InDel loci from 80 lines (40 resistant and 40 susceptible) also showed that all 40 R lines had no InDel mutation while 39 out of 40 S lines matched the two types of loci. Thus re-Bol037156 was identified as a likely candidate gene for FOC1 in cabbage.

Conclusions

This work may lay the foundation for marker-assisted selection as well as for further function analysis of the FOC1 gene.     

Mapping loci controlling vernalization requirement and flowering time in Brassica napus

Rapeseed cultivars (Brassica napus L.) can be classified into annual and biennial groups according to their requirement for vernalization in order to induce flowering. The genetic control of these phenotypic differences is not well understood, but this information could be valuable for the design of breeding approaches to accelerate rapeseed improvement. In order to map loci controlling this variation, a doubled haploid population, derived from a cross between annual and biennial cultivars, was evaluated for vernalization requirement and days-to-flowering in a replicated field experiment using three treatments: no vernalization, 4 weeks of vernalization and 8 weeks of vernalization. A linkage map of 132 RFLP loci was used to locate loci controlling these traits. Marker segregation in one region of linkage group 9 was strongly associated with the annual/biennial growth habit in the unvernalized treatment and with days-to-flowering in all three treatments. Two other regions with smaller effects on days-to-flowering were also identified.     

Development of 5123 Intron-Length Polymorphic Markers for Large-Scale Genotyping Applications in Foxtail Millet

Generating genomic resources in terms of molecular markers is imperative in molecular breeding for crop improvement. Though development and application of microsatellite markers in large-scale was reported in the model crop foxtail millet, no such large-scale study was conducted for intron-length polymorphic (ILP) markers. Considering this, we developed 5123 ILP markers, of which 4049 were physically mapped onto 9 chromosomes of foxtail millet. BLAST analysis of 5123 expressed sequence tags (ESTs) suggested the function for ∼71.5% ESTs and grouped them into 5 different functional categories. About 440 selected primer pairs representing the foxtail millet genome and the different functional groups showed high-level of cross-genera amplification at an average of ∼85% in eight millets and five non-millet species. The efficacy of the ILP markers for distinguishing the foxtail millet is demonstrated by observed heterozygosity (0.20) and Nei''s average gene diversity (0.22). In silico comparative mapping of physically mapped ILP markers demonstrated substantial percentage of sequence-based orthology and syntenic relationship between foxtail millet chromosomes and sorghum (∼50%), maize (∼46%), rice (∼21%) and Brachypodium (∼21%) chromosomes. Hence, for the first time, we developed large-scale ILP markers in foxtail millet and demonstrated their utility in germplasm characterization, transferability, phylogenetics and comparative mapping studies in millets and bioenergy grass species.     

SNP Discovery and Linkage Map Construction in Cultivated Tomato

Few intraspecific genetic linkage maps have been reported for cultivated tomato, mainly because genetic diversity within Solanum lycopersicum is much less than that between tomato species. Single nucleotide polymorphisms (SNPs), the most abundant source of genomic variation, are the most promising source of polymorphisms for the construction of linkage maps for closely related intraspecific lines. In this study, we developed SNP markers based on expressed sequence tags for the construction of intraspecific linkage maps in tomato. Out of the 5607 SNP positions detected through in silico analysis, 1536 were selected for high-throughput genotyping of two mapping populations derived from crosses between ‘Micro-Tom’ and either ‘Ailsa Craig’ or ‘M82’. A total of 1137 markers, including 793 out of the 1338 successfully genotyped SNPs, along with 344 simple sequence repeat and intronic polymorphism markers, were mapped onto two linkage maps, which covered 1467.8 and 1422.7 cM, respectively. The SNP markers developed were then screened against cultivated tomato lines in order to estimate the transferability of these SNPs to other breeding materials. The molecular markers and linkage maps represent a milestone in the genomics and genetics, and are the first step toward molecular breeding of cultivated tomato. Information on the DNA markers, linkage maps, and SNP genotypes for these tomato lines is available at http://www.kazusa.or.jp/tomato/.