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.     

Genetic analysis of hybrid seed formation ability of Brassica rapa in intergeneric crossings with Raphanus sativus

A hybridization barrier leads to the inability of seed formation after intergeneric crossings between Brassica rapa and Raphanus sativus. Most B. rapa lines cannot set intergeneric hybrid seeds because of embryo breakdown, but a B. rapa line obtained from turnip cultivar ‘Shogoin-kabu’ is able to produce a large number of hybrid seeds as a maternal parent by crossings with R. sativus. In ‘Shogoin-kabu’ crossed with R. sativus, developments of embryos and endosperms were slower than those in intraspecific crossings, but some of them grew to mature seeds without embryo breakdown. Intergeneric hybrid seeds were obtained in a ‘Shogoin-kabu’ line at a rate of 0.13 per pollinated flower, while no hybrid seeds were obtained in a line developed from Chinese cabbage cultivar ‘Chiifu’. F₁ hybrid plants between the lines of ‘Shogoin-kabu’ and ‘Chiifu’ set a larger number of hybrid seeds per flower, 0.68, than both the parental lines. Quantitative trait loci (QTLs) for hybrid seed formation were analyzed after intergeneric crossings using two different F₂ populations derived from the F₁ hybrids, and three QTLs with significant logarithm of odds scores were detected. Among them, two QTLs, i.e., one in linkage group A10 and the other in linkage group A01, were detected in both the F₂ populations. These two QTLs had contrary effects on the number of hybrid seeds. Epistatic interaction between these two QTLs was revealed. Possible candidate genes controlling hybrid seed formation ability in QTL regions were inferred using the published B. rapa genome sequences.     

Functional Divergence and Evolutionary Dynamics of the Putative AAAP Gene Family in Brassica rapa

The amino acid/auxin permease (AAAP) protein family is ubiquitously present in almost all eukaryotic species and functions in various aspects of growth and development. To investigate the evolution of AAAP proteins, here 83 AAAP genes in Brassica rapa were identified, and their sequence features, and evolutionary relationships were analyzed using in silico methods. According to the phylogenetic analysis, the AAAP genes of B. rapa are divided into six clades, and these clades share relatively similar sequence features, including gene structures, conserved motifs, and domain organizations. Synteny mapping strongly suggested that segmental duplications could be responsible for the expansion of this family. Adaptive evolution analysis demonstrated that most of AAAP proteins were subject to purifying selection. However, the site Tyr²⁵⁷ on eight AAAP proteins from clade 2b underwent significant positive selection. Functional divergent analysis showed that type I functional divergence coefficients (θ _I ) were significantly greater than zero in six pair-wise comparisons. However, functional divergence sites (Q _k ?>?0.95) found only in the AAAP I/II and AAAP I/III comparisons were localized mainly to the trans-membrane (TM) regions, suggesting highly divergent TM structures between these groups might be associated with group-specific functions. Our results could provide a valuable clue for further investigations of the evolutionary history and biological functions of the AAAP genes in B. rapa.     

Glucosinolate biosynthetic genes in Brassica rapa

Glucosinolates (GS) are a group of amino acid-derived secondary metabolites found throughout the Cruciferae family. Glucosinolates and their degradation products play important roles in pathogen and insect interactions, as well as in human health. In order to elucidate the glucosinolate biosynthetic pathway in Brassica rapa, we conducted comparative genomic analyses of Arabidopsis thaliana and B. rapa on a genome-wide level. We identified 102 putative genes in B. rapa as the orthologs of 52 GS genes in A. thaliana. All but one gene was successfully mapped on 10 chromosomes. Most GS genes exist in more than one copy in B. rapa. A high co-linearity in the glucosinolate biosynthetic pathway between A. thaliana and B. rapa was also established. The homologous GS genes in B. rapa and A. thaliana share 59-91% nucleotide sequence identity and 93% of the GS genes exhibit synteny between B. rapa and A. thaliana. Moreover, the structure and arrangement of the B. rapa GS (BrGS) genes correspond with the known evolutionary divergence of B. rapa, and may help explain the profiles and accumulation of GS in B. rapa.     

The use of species-specific DNA markers for assessing alien chromosome transfer in Brassica rapa and Brassica oleracea-monosomic additions of Raphanus sativus

Monosomic addition lines (MALs) are useful materials not only for cytogenetic and molecular genetic studies but also for plant breeding as gene sources. In our previous study, two MALs in the tribe Brassiceae were developed, one being Raphanus sativus lines with alien chromosomes of Brassica rapa (B. rapa-monosomic addition lines; BrMALs) and the second being those with alien chromosomes of Brassica oleracea (B. oleracea-monosomic addition lines; BoMALs). We developed species-specific DNA markers from the genomic sequences of B. rapa and B. oleracea comparing them with those of R. sativus, and identified chromosomes added in BrMALs and BoMALs using these markers. It was revealed that eight types of BrMALs have seven chromosomes of B. rapa and seven types of BoMALs have six chromosomes of B. oleracea. Furthermore, chromosome breakage and homoeologous recombination were suggested to have occurred in some MALs. The developed species-specific DNA markers are considered to be useful for producing MALs and also for assessing chromosome abnormality in MALs.     

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.     

Genome-wide analysis and stress-responsive expression of CCCH zinc finger family genes in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Background

Ubiquitous CCCH nucleic acid-binding motif is found in a wide-variety of organisms. CCCH genes are involved in plant developmental processes and biotic and abiotic stress responses. Brassica rapa is a vital economic crop and classical model plant of polyploidy evolution, but the functions of CCCH genes in B. rapa are unclear.

Results

In this study, 103 CCCH genes in B. rapa were identified. A comparative analysis of the chromosomal position, gene structure, domain organization and duplication event between B. rapa and Arabidopsis thaliana were performed. Results showed that CCCH genes could be divided into 18 subfamilies, and segmental duplication might mainly contribute to this family expansion. C-X_7/8-C-X₅-C₃-H was the most commonly found motif, but some novel CCCH motifs were also found, along with some loses of typical CCCH motifs widespread in other plant species. The multifarious gene structures and domain organizations implicated functional diversity of CCCH genes in B. rapa. Evidence also suggested functional redundancy in at least one subfamily due to high conservation between members. Finally, the expression profiles of subfamily-IX genes indicated that they are likely involved in various stress responses.

Conclusion

This study provides the first genome-wide characterization of the CCCH genes in B. rapa. The results suggest that B. rapa CCCH genes are likely functionally divergent, but mostly involved in plant development and stress response. These results are expected to facilitate future functional characterization of this potential RNA-binding protein family in Brassica crops.

     

Preferential retention,expression profile and potential functional diversity analysis of HD-Zip gene family in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Homeodomain-Leu zipper (HD-Zip) gene family performs important biological functions related to organ development, photomorphogenesis and abiotic stress response in higher plants. However, systematic analysis of HD-Zip genes in Brassica rapa has not been performed. In the present study, a bioinformatics approach was used to identify and characterize the BraHD-Zip gene family in B. rapa. A total of 88 members were identified. All putative BraHD-Zip proteins contained a clear HD and LZ combined domain. Eighty-seven BraHD-Zips were non-randomly located on ten chromosomes. This gene family was mainly expanded following the whole genome triplication event and was preferentially over-retained relative to its neighboring genes in B. rapa. On phylogenetic analysis, the BraHD-Zips could be categorized into four distinct major groups (I–IV). Each group exhibited variant gene structures and motif distributions. Some syntenic orthologous gene pairs presented diverse expression profiles, which indicate that these gene pairs may be involved in the development of new functions during evolution. In summary, our analysis provided genome-wide insights into the expansion, preferential retention, expression profiles and functional diversity of BraHD-Zip genes following whole genome triplication in B. rapa.     

Genome-wide comparative analysis of NBS-encoding genes between Brassica species and Arabidopsis thaliana

Background

Plant disease resistance (R) genes with the nucleotide binding site (NBS) play an important role in offering resistance to pathogens. The availability of complete genome sequences of Brassica oleracea and Brassica rapa provides an important opportunity for researchers to identify and characterize NBS-encoding R genes in Brassica species and to compare with analogues in Arabidopsis thaliana based on a comparative genomics approach. However, little is known about the evolutionary fate of NBS-encoding genes in the Brassica lineage after split from A. thaliana.

Results

Here we present genome-wide analysis of NBS-encoding genes in B. oleracea, B. rapa and A. thaliana. Through the employment of HMM search and manual curation, we identified 157, 206 and 167 NBS-encoding genes in B. oleracea, B. rapa and A. thaliana genomes, respectively. Phylogenetic analysis among 3 species classified NBS-encoding genes into 6 subgroups. Tandem duplication and whole genome triplication (WGT) analyses revealed that after WGT of the Brassica ancestor, NBS-encoding homologous gene pairs on triplicated regions in Brassica ancestor were deleted or lost quickly, but NBS-encoding genes in Brassica species experienced species-specific gene amplification by tandem duplication after divergence of B. rapa and B. oleracea. Expression profiling of NBS-encoding orthologous gene pairs indicated the differential expression pattern of retained orthologous gene copies in B. oleracea and B. rapa. Furthermore, evolutionary analysis of CNL type NBS-encoding orthologous gene pairs among 3 species suggested that orthologous genes in B. rapa species have undergone stronger negative selection than those in B .oleracea species. But for TNL type, there are no significant differences in the orthologous gene pairs between the two species.

Conclusion

This study is first identification and characterization of NBS-encoding genes in B. rapa and B. oleracea based on whole genome sequences. Through tandem duplication and whole genome triplication analysis in B. oleracea, B. rapa and A. thaliana genomes, our study provides insight into the evolutionary history of NBS-encoding genes after divergence of A. thaliana and the Brassica lineage. These results together with expression pattern analysis of NBS-encoding orthologous genes provide useful resource for functional characterization of these genes and genetic improvement of relevant crops.

Electronic supplementary material

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

Characterization of self-incompatibility genes in the intergeneric hybrid xBrassicoraphanus

Brassica rapa and Raphanus sativus are strictly self-incompatible (SI) plants; however, xBrassicoraphanus (baemoochae) is an intergeneric hybrid synthesized following hybridization of B. rapa and R. sativus that is self-compatible (SC). Self-incompatibility in Brassicaceae is controlled by multiple alleles of the S-locus. Two S-locus genes, SRK (S-locus receptor kinase) and SP11/SCR (S-locus protein 11 or S-locus cysteine-rich), have been reported to date, both of which are classified into class I and II. In this study, we investigate if there is an alteration in the structure or the expression in SRK or SP11 genes behind the alteration of SI to SC in baemoochae. Class I and II SRK were identified by PCR of the genomic DNA of baemoochae using SRK-specific universal primers. Cloning and sequencing of both classes of SRK was conducted and compared with SRK genes of parental species. Comparison analysis showed no genomic alterations and both of them showed expression in the stigma. Similarly, SP11 genes also showed no genomic alterations and normally are expressed in the anther. Other SI-related genes (MLPK, ARC1, THL) also showed normal expression in the stigma and anther. Taken together, these results revealed that no structural/gene expression change in these genes was responsible for the breakdown of SI in baemoochae. Rather, the transformation from SI parents to SC descendants could be responding to epigenetic mechanisms.