Next Article in Journal
Long Non-Coding RNAs at the Chromosomal Risk Loci Identified by Prostate and Breast Cancer GWAS
Previous Article in Journal
Formation of the Codon Degeneracy during Interdependent Development between Metabolism and Replication
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 12
Issue 12
10.3390/genes12122027
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Carotenoid Biosynthetic Genes in Cabbage: Genome-Wide Identification, Evolution, and Expression Analysis

by
Wenxue Cao
1,2,
Peng Wang
1,
Limei Yang
1,
Zhiyuan Fang
1,
Yangyong Zhang
1,
Mu Zhuang
1,
Honghao Lv
1,
Yong Wang
1 and
Jialei Ji
1,*
1
Key Laboratory of Biology and Genetic Improvement of Horticultural Crops, Ministry of Agriculture, Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, China
2
College of Horticulture and Landscape Architecture, Hunan Agricultural University, 1 Nongda Road, Changsha 410128, China
*
Author to whom correspondence should be addressed.
Genes 2021, 12(12), 2027; https://0-doi-org.brum.beds.ac.uk/10.3390/genes12122027
Submission received: 2 December 2021 / Revised: 16 December 2021 / Accepted: 17 December 2021 / Published: 20 December 2021
(This article belongs to the Section Plant Genetics and Genomics)
Browse Figures
Review Reports Versions Notes

Abstract

:
Carotenoids are natural functional pigments produced by plants and microorganisms and play essential roles in human health. Cabbage (Brassica oleracea L. var. capitata L.) is an economically important vegetable in terms of production and consumption. It is highly nutritious and contains β-carotene, lutein, and other antioxidant carotenoids. Here, we systematically analyzed carotenoid biosynthetic genes (CBGs) on the whole genome to understand the carotenoid biosynthetic pathway in cabbage. In total, 62 CBGs were identified in the cabbage genome, which are orthologs of 47 CBGs in Arabidopsis thaliana. Out of the 62 CBGs, 46 genes in cabbage were mapped to nine chromosomes. Evolutionary analysis of carotenoid biosynthetic orthologous gene pairs among B. oleracea, B. rapa, and A. thaliana revealed that orthologous genes of B. oleracea underwent a negative selection similar to that of B. rapa. Expression analysis of the CBGs showed functional differentiation of orthologous gene copies in B. oleracea and B. rapa. Exogenous phytohormone treatment suggested that ETH, ABA, and MeJA can promote some important CBGs expression in cabbage. Phylogenetic analysis showed that BoPSYs exhibit high conservatism. Subcellular localization analysis indicated that BoPSYs are located in the chloroplast. This study is the first to study carotenoid biosynthesis genes in cabbage and provides a basis for further research on carotenoid metabolic mechanisms in cabbage.

1. Introduction

Cabbage (Brassica oleracea L. var. capitata L., 2n = 18) is a member of the family Cruciferae and one of the most economically important vegetable crops cultivated worldwide [1]. Cabbage is rich in dietary fiber, vitamin B1 (VB1), vitamin B (VB2), calcium, and iron. Most importantly, cabbage also contains β-carotene, lutein, and other antioxidant carotenoids. It is one of the best vegetables recommended by the World Health Organization [2,3].
Carotenoids are naturally occurring pigments mainly found in plants, algae, and photosynthetic bacteria. Carotenoid pigments are mostly C40 lipophilic isoprenoids, comprising eight isoprene units joined in a head-to-tail fashion, and the central unit has a reverse connection [4]. The number of conjugated double bonds determines the color of carotenoids; the more the number of conjugated double bonds, the deeper the red color [5]. So far, over 800 carotenoids have been found in nature [6]. Carotenoids are divided into two main groups based on whether they contain oxygen or not: carotenes and xanthophylls. Carotenes are composed of carbon and hydrogen atoms and include phytoene, lycopene, α-carotene, and β-carotene. On the other hand, xanthophylls are oxygenated hydrocarbons and include lutein, zeaxanthin, antheraxanthin, and violaxanthin [7,8]. Carotenoids are one of the components of the photosystem elements. They absorb light in the blue range of the spectrum and then transfer energy to chlorophylls in antenna complexes [4]. In plants, carotenoids protect proteins from excessive light incidence via thermal dissipation and free radical detoxification [9,10]. Carotenoids also act as vitamin A precursors and antioxidants, which play a vital role in human health and nutrition [11,12]. β-carotene is provitamin A, and its deficiency in a diet causes symptoms, such as night-blindness, keratomalacia, and xerophthalmia in humans [13]. Vitamin A deficiency also exacerbates afflictions like diarrhea, measles, and respiratory diseases [14,15]. α-carotene and β-carotene also serve as antioxidants in human health. For example, studies have shown that supplementing the human diet with carotenoids can reduce the risk of developing macular degeneration and neurodegenerative diseases [10,16].
Carotenoids are mainly synthesized from isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP) produced by the plastidic methylerythritol 4-phosphate (MEP) pathway [17,18] (Figure S1). The MEP pathway is localized in plastids and uses glyceraldehyde 3-phosphate and pyruvate as initial substrates to form deoxy-D-xylulose 5-phosphate (DXP). Subsequently, IPP and DMAP are generated by DXP synthase (DXS) and DXP reductoisomerase (DXR) and undergo a sequential series of condensation reactions to produce geranylgeranyl diphosphate (GGPP) [19]. Two GGPP molecules are condensed to phytoene by phytoene synthase (PSY). This is the first rate-limiting step of carotenoid biosynthesis, and PSY is the rate-limiting enzyme in this process [20,21]. Phytoene is further converted to lycopene by phytoene desaturase (PDS) and ζ-carotene desaturase (ZDS) [22,23]. Since carotenoids mainly exist in nature as trans structures, isomerization is key in carotenoid biosynthesis. Carotenoid isomerase (CRTISO) converts tetra-cis-lycopene to all-trans-lycopene [24,25,26]. The carotenoid pathway branches at the cyclization of lycopene, generating diverse carotenoids distinguished by different cyclic end groups. α-carotene and β-carotene are produced by lycopene cyclase (ε-LCY and β-LCY), catalyzing all-trans-lycopene [27,28], and then further hydroxylated to produce xanthophylls, including lutein and zeaxanthin [29]. Antheraxanthin and violaxanthin are sequentially synthesized by zeaxanthin epoxidase (ZEP). Violaxanthin is converted to neoxanthin by neoxanthin synthase (NSY). Neoxanthin can also be used to synthesize abscisic acid (ABA) [30,31].
Several transcription factors and genes participate in the regulatory network of carotenoid biosynthesis. Three PSY genes (PSY1, PSY2, and PSY3) were identified in tomato (Solanum lycopersicum). PSY1 was mapped on chromosome 3 and is up-regulated in ripening tomato fruit. PSY2 is specifically expressed in leaf tissues and is implicated in carotenoid synthesis in the leaves. PSY3 is up-regulated in roots under stress conditions [32,33,34,35]. Additionally, several transcription factors directly or indirectly participate in carotenoid synthesis and metabolism. These transcription factors include MYB family, BBX family, NAC family, MADS-Box family, and AP2/ERF family among others [36,37,38,39,40]. Environmental factors such as light, temperature, and CO2 concentration also influence carotenoid synthesis to a great extent [9,41,42].
Although more in-depth studies on carotenoids synthesis have been carried out in Arabidopsis thaliana and tomato, there are few reports on the synthesis and metabolism of carotenoids in cabbage. In order to understand the mechanism of cabbage carotenoid synthesis more deeply and use modern molecular techniques to improve the carotenoids content in cabbage, we applied comparative genomics to identify cabbage carotenoids synthesis genes in the whole genome. Further, we determined the phylogenetic characteristics, chromosome location, non-synonymous/synonymous substitution (Ka/Ks), gene expression, and subcellular localization of genes related to carotenoid biosynthesis. The present study uncovers the genetic mechanisms underlying carotenoid biosynthesis in cabbage, thus laying a basis for future studies on carotenoid biosynthesis in plants.

2. Materials and Methods

2.1. Data Resource

The genomic and annotation data of A. thaliana, B. oleracea, and Brassica rapa were retrieved from the BRAD database (http://brassicadb.cn, accessed on 20 March 2021) and TAIR database (http://www.arabidopsis.org, accessed on 20 March 2021). The carotenoid biosynthetic gene (CBG) sequences of A. thaliana were acquired from the KEGG pathway database (http://www.genome.jp/kegg/pathway.html, accessed on 18 April 2021).

2.2. Identification and Analysis of Orthologs between B. oleracea and A. thaliana

To identify the CBGs and protein sequences of cabbage, we performed BLAST search on the cabbage genome and proteome database using A. thaliana carotenoid biosynthetic genes and protein sequences, with a cutoff E-value ≤ 10−10 and coverage ≥ 0.75. Syntenic orthologous genes between A. thaliana and B. oleracea were identified based on the collinearity of flanking genes sequence similarity (E-value ≤ 10−20). The specific distribution of B. oleracea carotenoid biosynthetic genes (BoCBGs) on the chromosome was analyzed using MapChart software (http://mapinspect.software.informer.com/, accessed on 24 April 2021).

2.3. Non-Synonymous/Synonymous Substitution (Ka/Ks) Ratios of Gene Pairs among B. oleracea, A. thaliana, and B. rapa

PAML software [43] was used to calculate the ratio of Ka/Ks rates of all orthologous gene pairs to estimate the selection mode of CBGs between B. oleracea, A. thaliana, and B. rapa. Ka/Ks ratio less than 1 indicates pure selection, ratio equal to 1 represents neutral selection, and ratio greater than 1 indicates positive selection.

2.4. Expression Analysis of Carotenoid Biosynthetic Genes in B. oleracea

The expression pattern of BoCBGs was investigated using RNA-Seq data [44]. Seven tissues, including root, leaf, stem, flower, bud, callus, and silique of B. oleracea accession ‘02-12’ were used for RNA extraction. Total RNA was isolated from the cabbage tissues using TIANGEN RNAprep Pure Plant Kit and then reverse-transcribed to cDNA using PrimeScript™ RT reagent kit (TaKaRa, Kyoto, Japan), according to the manufacturer’s instructions.

2.5. Exogenous Phytohormone Spraying of Cabbage Leaves

The leaves of cabbage were sprayed with different phytohormones to explore their effect on CBGs expression. Four exogenous phytohormones: ETH, ABA, SA, and MeJA at concentrations of 100 mg/L, 50 mg/L, 100 μM, and 100 μM, respectively, were sprayed. The leaves were sampled at 2, 4, 6, 12, and 24 h after spraying. The control was sprayed with double-distilled water alone. The relative expressions of BoPSY.1, BoPDS.1, BoZDS, and BoLYC genes were determined using qRT-PCR. The primers for qRT-PCR are listed in Table S1.

2.6. Phylogenetic Analysis

PSY protein sequences from 24 plant species were analyzed. MEGA6.0 was used to construct the phylogenetic tree, and Poisson correction model was used for distance computation. Node support was assessed by 1000 bootstrap replicates.

2.7. Subcellular Localization of PSY in B. oleracea

The coding sequences of BoPSY.1 (Bol019820), BoPSY.2 (Bol034439), and BoPSY.3 (Bol021326) were amplified using primers PSY1-1, PSY1-2, PSY2, and PSY3 (Table S1), and then inserted into a pBWA(V)HS-GFP vector (BIORUN Technologies Company, Wuhan, China). The fusion constructs were introduced into tobacco leaf epidermis as previously described [45]. The fluorescence signals were observed using Olympus FV3000 confocal laser-scanning microscope.

3. Results

3.1. Identification of Carotenoid Biosynthetic Genes in B. oleracea

Genes associated with the carotenoid synthesis pathway in A. thaliana were analyzed using KEGG pathway database and TAIR database. A total of 47 genes were implicated in carotenoid biosynthesis in A. thaliana, of which 21 genes were shown to participate in the MEP pathway and 26 gene-encoding carotenoid biosynthetic enzymes (Table 1). A total of 62 CBGs were identified in cabbage, and 9 AtCBGs (GGPS4, GGPS7, GGPS6, GGPS8, GGPS9, GGPS10, GGPS11, GGPS12, and LUT1) showed no cabbage orthologs. Among the 62 BoCBGs, 58 were syntenic orthologs of the AtCBGs, and 4 BoCBGs had no syntenic relationships (Figure 1).
In this study, BoCBGs were divided into three sub-genomes (LF, MF1, and MF2). Among the 58 syntenic homologous genes, there are 27, 17, and 14 genes in LF, MF1, and MF2, respectively (Table 1). In addition, a collinearity relationship was observed between B. oleracea and B. rapa. In total, 43 pairs of carotenoid biosynthetic genes were mapped on chromosomes using Circos software (Figure 1).

3.2. Genomic Distribution on Chromosomes

Based on the genome annotation file, 62 BoCBGs were mapped to nine chromosomes. Among them, 46 were located on the chromosomes and 16 on the scaffold. The distribution of the 46 genes mapped on the chromosomes is shown in Figure 2, with 6, 3, 10, 4, 7, 2, 5, 6, and 3 genes located on chromosome C01–C09, respectively.

3.3. Comparative Evolutionary Analyses of Orthologous Gene Pairs for Carotenoid Biosynthetic Genes

Ka/Ks value is the ratio of non-synonymous mutation (Ka) to synonymous mutation (Ks) of genes. In this study, PAML software [43] was used to calculate the Ka/Ks of homologous gene pairs to explore the evolutionary selection mode of CBGs in A. thaliana, B. oleracea and B. rapa. Among the BoCBGs or BrCBGs, there were three copies of AtPSY gene, and two copies of AtPDS3, AtLUT2, and AtCHY1 genes. Eventually, 20 orthologous gene pairs were selected from A. thaliana, B. oleracea, and B. rapa (Figure 3). Relative to A. thaliana, the Ka/Ks values of carotenoid homologous genes in B. oleracea and B. rapa were roughly the same, although the Ka/Ks ratios of BoGGPS2, BoCHY1.1, BoCHY1.2, and BoCHY2.1 in B. oleracea genome were slightly higher than that of B. rapa genome. Also, BoLUT2.1 and BoLUT5.1 were lower than those in B. rapa. The Ka/Ks values shown in Figure 3 are all less than 1, indicating that these genes are subjected to purify selection in evolution.

3.4. Expression Analysis of Orthologous for Carotenoid Biosynthetic Genes among B. oleracea, B. rapa and A. thaliana

In this study, 20 CBGs in A. thaliana were studied to explore the differential expression of homologous genes in B. oleracea and B. rapa genomes. Four genes (BoPSY.1, BoPSY.2, BrPSY1, and BrPSY2) were homologous to At5G17230 gene. BoPSY.1 was expressed in all parts of B. oleracea, except roots, while BoPSY.2 was expressed only in flowers and callus (Figure 4). BrPSY1 was expressed in stems, leaves, flowers, and silique of B. rapa, while BrPSY2 was almost exclusively expressed in flowers. The BoPDS3.1 gene in B. oleracea was homologous to At4G14210 and expressed in all tissues, but the BrPDS3.1 was only expressed in silique of B. rapa. At5G57030 exist in two copies in B. oleracea and B. rapa. BoLUT2.1 expression was detected in all parts except roots, while BrLUT2.1 expression was only observed in flowers and silique. BoZEP.1, BoZEP.1, BrZEP1, and BrZEP2 are two copies of At5G67030. The expression level of BoZEP.1 was highest in silique and the lowest in leaves, while that of BrZEP1 exhibited an opposite trend. BoNCED3.2 was expressed in all tissues of B. oleracea, but BrNCED3.2 was detected in the callus of B. rapa. Analysis of differential expression between carotenoid homologous genes revealed functional variation among the homologous genes duplicated in B. oleracea and B. rapa during evolution.
The tissue expression of important genes in the carotenoid pathway was analyzed to further understand the expression of CBGs in B. oleracea (Figure 5). The results showed that BoGGPS1 was expressed in all organs, including roots, leaves, stems, flowers, buds, silique, and callus, while BoGGPS2 was only expressed in flower buds. The expression of BoPSY.1 was highest in callus and flower buds, followed by flowers, stems, leaves, and silique. Almost no expression of BoPSY.1 was detected in the roots. Meanwhile, the expression of BoPSY.2 was highest in flowers, followed by buds and callus, and almost no expression was detected in other tissues. BoPSY.3 was expressed in all tissues, with the highest expression detected in the stem. On the other hand, BoPDS3.1 and BoPDS3.2 were expressed in all tissues, however, BoPDS3.1 was mainly expressed in callus, while BoPDS3.2 was mainly expressed in the flowers. The expression of BoZDS was highest in the roots, while that of BoLYC was highest in the flowers. Gene expression analysis revealed an extensive variance between paralogs of each CBG.

3.5. BoPSY.1, BoPDS.1, BoZDS, and BoLYC Respond to Exogenous Phytohormone Treatments

Specific concentrations of ETH, ABA, SA, and MeJA were sprayed on cabbage leaves at the seedling stage to explore their effects on CBGs (BoPSY.1, BoPDS.1, BoZDS, and BoLYC). ETH at 100 mg/L promoted the expression of the four genes (Figure 6). Specifically, the expression of BoPSY.1 increased at 2 h after treatment and was the highest at 12 h. Meanwhile, the expressions of BoPDS.1, BoZDS, and BoLYC were highest at 4 and 6 h and gradually decreased at 12 and 24 h. ABA also promoted the carotenoid pathway in cabbage. After spraying 50 mg/L ABA, the expression of BoPSY.1 increased gradually and peaked after 24 h, while the expression of BoZDS increased first and then decreased. SA at 100 μmol/L promoted the expression of BoPSY.1, but had minimal effect on the other three genes. MeJA at 100 μmol/L increased the expressions of BoPSY.1 and BoPDS.1 to a certain extent, which then decreased and peaked after 24 h. Meanwhile, the expression of BoZDS gene peaked after 2 h of MeJA treatment. Collectively, these results indicate that ETH, ABA, and MeJA can promote CBGs expression in cabbage.

3.6. Evolution of PSY Genes

The protein sequences of PSY in 24 species were analyzed using MEGA 6.0 software to investigate the evolutionary relationship (Figure 7 and Table S2). As shown in Figure 7, the phylogenetic tree was consistent with the taxonomic conclusion, indicating that PSY has a high conservatism, which provides a theoretical basis for plant genetic molecular evolution. The PSY proteins of 24 species were divided into two categories. Haematococcus clustered into one category alone. The protein cluster of the other plant PSY belonged to angiosperms and was further divided into five subgroups: rice, maize, and sorghum were grouped under Commelinidae; alfalfa, carrot, strawberry, and apple belonged to Rosidae; lilium brownii and narcissus belonged to the Liliidae; Dilleniidae was divided into two groups, one group includes Arabidopsis and cabbage (Cruciferae), papaya (Carica Linn), pumpkin, and watermelon (Cucurbitaceae), the other group includes persimmon and kiwifruit clustered with Asteridae; Asteridae includes tomato, chilli, medlar, sweet potato, cape jasmine, Osmanthus, and agastache. Notably, BoPSY.1 and BoPSY.3 clustered on one branch and were more closely related to AtPSY than BoPSY.2 (Figure 7). Motif-based sequence analysis showed that all PSY proteins derived from angiosperms contain six conserved motifs, but the Haematococcus PSY only contains four motifs, lacking motif3 and motif6.

3.7. Subcellular Localization of BoPSY.1, BoPSY.2, and BoPSY.3

The results of Plant-PLoc subcellular localization prediction indicated that BoPSY proteins are localized in the chloroplast. Subcellular location of BoPSYs in tobacco (Nicotiana benthamiana L.) was examined using Confocal Laser Scanning Microscopy (FV3000). All three fluorescence signals of BoPSY-GFP fusion vector (BoPSY.1-GFP, BoPSY.2-GFP, and BoPSY.3-GFP) were observed in the chloroplast, consistent with the prediction of subcellular localization (Figure 8).

4. Discussion

4.1. Characterization of Carotenoid Biosynthetic Genes in B. oleracea

Whole genome duplication provides rich genetic material for the expansion of gene families or the evolution of new genes in plants [46,47]. A whole-genome triplication (WGT) event occurred in B. oleracea after its divergence from A. thaliana [44]. There are 47 genes associated with carotenoid synthesis in Arabidopsis. In this study, 62 CBGs were identified in cabbage, among which 58 were syntenic orthologs of Arabidopsis, and only four had no syntenic relationships, suggesting that they could have evolved from some duplicate genes of Arabidopsis. However, nine AtCBGs (GGPS4, GGPS7, GGPS6, GGPS8, GGPS9, GGPS10, GGPS11, GGPS12, and LUT1) showed no B. oleracea orthologs, indicating that they might have been lost during evolution of B. oleracea. Notably, B. oleracea and A. thaliana belong to the Brassicaceae family. Phylogenetic analysis revealed that the triplicated B. oleracea genome segments diverged from a common ancestor soon after the divergence of the A. thaliana and Brassica lineages [48,49]. The multiple copies of the BoCBGs that are syntenic to genes in A. thaliana were generated from the WGT. The fragment with the highest gene densities was in the sub-genome LF, the fragment with moderate gene densities was in the sub-genome MF1, and the fragment with the least genes was in the sub-genome MF2. There were 27, 17, and 14 genes located in LF, MF1, and MF2, respectively. The distribution of CBGs was in line with the gene separation status in the whole genomic sequences [44].
Further, we compared the non-synonymous/synonymous (Ka/Ks) ratio of orthologous gene pairs between A. thaliana, B. oleracea, and B. rapa to analyze the evolutionary relationship of orthologous gene pairs for CBGs among the three species (Figure 3). The Ka/Ks ratio is a measure of selective pressures acting on genes. There were no significant differences in orthologous gene pairs between A. thaliana-B. oleracea and A. thaliana-B. rapa lineages. Therefore, we speculated that BoCBGs and BrCBGs could have undergone similar negative selection. However, comparing the expression differences between orthologous gene pairs in CBG revealed functional variation of a few gene pairs in B. oleracea and B. rapa.

4.2. Exogenous Phytohormones Regulate the Expression of Carotenoid Biosynthetic Genes

Regulation of carotenoid biosynthesis in plants is a complex process that is regulated by multilevel factors. Exogenous phytohormones affect the accumulation of carotenoids. Ethylene (ET) plays a vital role in plant ripening and carotenoid accumulation. Thus, inhibiting ET synthesis in fruits can hinder fruit ripening and lycopene accumulation [50]. Abscisic acid (ABA) is one of the most important plant phytohormones. Its synthetic precursors include neoxanthin and violaxanthin. These precursors regulate ABA synthesis, which in turn mediate the synthesis of carotenoids in plants [51,52]. Salicylic acid (SA) and Methyl Jasmonate (MeJA) also play an essential role in plant growth and stress resistance.
In the present study, different concentrations of exogenous phytohormones (ETH, ABA, SA, and MeJA) were sprayed on cabbage leaves to explore their effects on CBGs expression. The expression of four CBGs (BoPSY.1, BoPDS.1, BoZDS, and BoLYC) was analyzed. The result showed that ETH, ABA, and MeJA could promote the expression of CBGs in cabbage at specific concentrations. SA at 100 μmol/L did not significantly affect the expression of the four genes, consistent with the report of Gao, who concluded that the carotenoid synthesis pathway is insensitive to SA [53]. These results will guide the adoption of exogenous phytohormone spraying to increase carotenoid contents in cabbage.

4.3. Phytoene Synthase (PSY) Is a Key Enzyme in the Carotenoid Biosynthetic Pathway

PSY is the first special enzyme in the carotenoid biosynthesis pathway, catalyzing the condensation of two GGPP molecules into phytoene [54]. Overexpression of PSY increased carotenoid content and substantially improved β-carotene synthesis in canola seeds, cassava roots, potato tubers, and endosperms [55,56,57]. In tomato, PSY1 is highly expressed in fruits, PSY2 is essential for carotenoid synthesis in leaf tissues, and PSY3 potentially functions in roots [32,33,34,35].
In this study, BoPSY.1 and BoPSY.3 clustered on one branch, with BoPSY.2 in a separate branch. However, the three BoPSYs clustered on a specific branch in Arabidopsis, where a single PSY gene was shown to regulate phytoene synthesis in all tissues [21]. Three genes encoding PSY enzymes were identified in B. oleracea and B. rapa after WGT event [58]. In this study, three BoPSYs genes in cabbage were all localized in the chloroplast, which is similar to maize PSY1 [59]. However, examination of tissue expression of BoPSYs showed that they are differentially expressed in cabbage tissues (Figure 5). Thus, we speculate that the BoPSYs might function in different cabbage tissues, like the PSYs expression patterns in tomato [32,33,34,35]; however, further studies are needed to verify this hypothesis.

5. Conclusions

In total, 62 carotenoid biosynthetic genes were identified in cabbage. These genes in B. oleracea species underwent similar negative selection with those in B. rapa species. However, the differential expression pattern of the duplicated CBGs occurred after polyploidization. The expression of BoCBGs can be regulated by phytohormones. ETH, ABA, and MeJA can promote the expression of BoPSY.1, BoPDS.1, BoZDS, and BoLYC, which provides a basis for spraying phytohormones to improve the carotenoid content of cabbage. As the key enzyme in the carotenoid biosynthetic pathway, three BoPSY proteins were all localized in chloroplasts and exhibited different expression patterns in various cabbage tissues. This study uncovered carotenoid metabolic mechanisms in B. oleracea, proving a basis for developing cabbage varieties with high carotenoid content via genetic engineering.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/genes12122027/s1, Figure S1: Carotenoid synthesis pathway; Table S1: The primer sequences information used in this study; Table S2: The accession No. of amino acid sequences of PSY in 23 kinds of plants.

Author Contributions

Conceptualization, J.J.; methodology, J.J. and W.C.; software, J.J. and W.C.; validation, W.C., J.J. and P.W.; formal analysis, H.L.; investigation, L.Y.; resources, Z.F. and Y.Z.; data curation, W.C. and J.J.; writing—original draft preparation, W.C.; writing—review and editing, J.J.; visualization, J.J. and W.C.; supervision, H.L.; project administration, M.Z. and Y.W.; funding acquisition, J.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (32002034), National Key Research and Development Program (2017YFD0101804), the Science and Technology Innovation Program of the Chinese Academy of Agricultural Sciences (CAAS-ASTIP-2013-IVFCAAS), and the Modern Agro-Industry Technology Research System (CARS-25-B-01).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be available on reasonable request.

Acknowledgments

We acknowledge all the members of the research group for their helpful comments and inspiration.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lv, H.; Wang, Y.; Han, F.; Ji, J.; Fang, Z.; Zhuang, M.; Li, Z.; Zhang, Y.; Yang, L. A high-quality reference genome for cabbage obtained with SMRT reveals novel genomic features and evolutionary characteristics. Sci. Rep. 2020, 10, 12394. [Google Scholar] [CrossRef] [PubMed]
  2. Cartea, M.E.; Velasco, P. Glucosinolates in Brassica foods: Bioavailability in food and significance for human health. Phytochem. Rev. 2007, 7, 213–229. [Google Scholar] [CrossRef]
  3. Wiczkowski, W.; Szawara-Nowak, D.; Topolska, J. Red cabbage anthocyanins: Profile, isolation, identification, and antioxidant activity. Food Res. Int. 2013, 51, 303–309. [Google Scholar] [CrossRef]
  4. Alvarez, R.; Vaz, B.; Gronemeyer, H.; de Lera, A.R. Functions, therapeutic applications, and synthesis of retinoids and carotenoids. Chem. Rev. 2014, 114, 1–125. [Google Scholar] [CrossRef] [PubMed]
  5. Nisar, N.; Li, L.; Lu, S.; Khin, N.C.; Pogson, B.J. Carotenoid metabolism in plants. Mol. Plant 2015, 8, 68–82. [Google Scholar] [CrossRef] [Green Version]
  6. Ding, B.Y.; Niu, J.; Shang, F.; Yang, L.; Chang, T.Y.; Wang, J.J. Characterization of the geranylgeranyl diphosphate synthase gene in Acyrthosiphon pisum (Hemiptera: Aphididae) and its association with carotenoid biosynthesis. Front. Physiol. 2019, 10, 1398. [Google Scholar] [CrossRef] [PubMed]
  7. Bhosale, P.; Bernstein, P.S. Microbial xanthophylls. Appl. Microbiol. Biotechnol. 2005, 68, 445–455. [Google Scholar] [CrossRef] [PubMed]
  8. Rodrigo-Banos, M.; Garbayo, I.; Vilchez, C.; Bonete, M.J.; Martinez-Espinosa, R.M. Carotenoids from Haloarchaea and their potential in biotechnology. Mar. Drugs 2015, 13, 5508–5532. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Llorente, B.; Martinez-Garcia, J.F.; Stange, C.; Rodriguez-Concepcion, M. Illuminating colors: Regulation of carotenoid biosynthesis and accumulation by light. Curr. Opin. Plant Biol. 2017, 37, 49–55. [Google Scholar] [CrossRef]
  10. DellaPenna, D.; Pogson, B.J. Vitamin synthesis in plants: Tocopherols and carotenoids. Annu. Rev. Plant Biol. 2006, 57, 711–738. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Fassett, R.G.; Coombes, J.S. Astaxanthin in cardiovascular health and disease. Molecules 2012, 17, 2030–2048. [Google Scholar] [CrossRef] [PubMed]
  12. Mordente, A.; Guantario, B.; Meucci, E.; Silvestrini, A.; Lombardi, E.; Martorana, G.E.; Giardina, B.; Bohm, V. Lycopene and cardiovascular diseases: An update. Curr. Med. Chem. 2011, 18, 1146–1163. [Google Scholar] [CrossRef]
  13. Sommer, A.; Vyas, K.S. A global clinical view on vitamin A and carotenoids. Am. J. Clin. Nutr. 2012, 96, 1204S–1206S. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Silva, J.S.; Chaves, G.V.; Stenzel, A.P.; Pereira, S.E.; Saboya, C.J.; Ramalho, A. Improvement of anthropometric and biochemical, but not of vitamin A, status in adolescents who undergo Roux-en-Y gastric bypass: A 1-year follow up study. Surg. Obes. Relat. Dis. 2017, 13, 227–233. [Google Scholar] [CrossRef] [PubMed]
  15. Sandmann, G. Carotenoids of biotechnological importance. In Biotechnology of Isoprenoids, Advances in Biochemical Engineering/Biotechnology; Schrader, J., Bohlmann, J., Eds.; Springer: Berlin/Heidelberg, Germany, 2014; Volume 148, pp. 449–467. [Google Scholar] [CrossRef]
  16. Obulesu, M.; Dowlathabad, M.R.; Bramhachari, P.V. Carotenoids and Alzheimer’s disease: An insight into therapeutic role of retinoids in animal models. Neurochem. Int. 2011, 59, 535–541. [Google Scholar] [CrossRef]
  17. Eisenreich, W.; Bacher, A.; Arigoni, D.; Rohdich, F. Biosynthesis of isoprenoids via the non-mevalonate pathway. Cell. Mol. Life Sci. 2004, 61, 1401–1426. [Google Scholar] [CrossRef] [PubMed]
  18. Rodriguez-Concepcion, M.; Boronat, A. Elucidation of the methylerythritol phosphate pathway for isoprenoid biosynthesis in bacteria and plastids. A metabolic milestone achieved through genomics. Plant Physiol. 2002, 130, 1079–1089. [Google Scholar] [CrossRef] [Green Version]
  19. Aharoni, A.; O’Connell, A.P. Gene expression analysis of strawberry achene and receptacle maturation using DNA microarrays. J. Exp. Bot. 2002, 53, 2073–2087. [Google Scholar] [CrossRef] [Green Version]
  20. Scolnik, P.A.; Bartley, G.E. Nucleotide sequence of an Arabidopsis cDNA for phytoene synthase. Plant Physiol. 1994, 104, 1471. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Welsch, R.; Medina, J.; Giuliano, G.; Beyer, P.; Von Lintig, J. Structural and functional characterization of the phytoene synthase promoter from Arabidopsis thaliana. Planta 2003, 216, 523–534. [Google Scholar] [CrossRef]
  22. Naing, A.H.; Kyu, S.Y.; Pe, P.P.W.; Park, K.I.; Lee, J.M.; Lim, K.B.; Kim, C.K. Silencing of the phytoene desaturase (PDS) gene affects the expression of fruit-ripening genes in tomatoes. Plant Methods 2019, 15, 110. [Google Scholar] [CrossRef] [PubMed]
  23. Scolnik, P.A.; Bartley, G.E. Phytoene desaturase from Arabidopsis. Plant Physiol. 1993, 103, 1475. [Google Scholar] [CrossRef]
  24. Cazzonelli, C.I.; Pogson, B.J. Source to sink: Regulation of carotenoid biosynthesis in plants. Trends Plant Sci. 2010, 15, 266–274. [Google Scholar] [CrossRef] [PubMed]
  25. Chen, Y.; Li, F.; Wurtzel, E.T. Isolation and characterization of the Z-ISO gene encoding a missing component of carotenoid biosynthesis in plants. Plant Physiol. 2010, 153, 66–79. [Google Scholar] [CrossRef] [Green Version]
  26. Park, H.; Kreunen, S.S.; Cuttriss, A.J.; DellaPenna, D.; Pogson, B.J. Identification of the carotenoid isomerase provides insight into carotenoid biosynthesis, prolamellar body formation, and photomorphogenesis. Plant Cell 2002, 14, 321–332. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Cunningham, F.X.; Chamovitz, D.; Misawa, N.; Gantt, E.; Hirschberg, J. Cloning and functional expression in Escherichia coli of a cyanobacterial gene for lycopene cyclase, the enzyme that catalyzes the biosynthesis of β-carotene. FEBS Lett. 1993, 328, 130–138. [Google Scholar] [CrossRef] [Green Version]
  28. Cunningham, F.X.; Pogson, B.; Sun, Z.R.; McDonald, K.A.; DellaPenna, D.; Gantt, E. Functional analysis of the β and ε lycopene cyclase enzymes of Arabidopsis reveals a mechanism for control of cyclic carotenoid formation. Plant Cell 1996, 8, 1613–1626. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Yu, B.; Lydiate, D.J.; Young, L.W.; Schafer, U.A.; Hannoufa, A. Enhancing the carotenoid content of Brassica napus seeds by downregulating lycopene ε cyclase. Transgenic Res. 2008, 17, 573–585. [Google Scholar] [CrossRef]
  30. Gomez-Garcia Mdel, R.; Ochoa-Alejo, N. Biochemistry and molecular biology of carotenoid biosynthesis in chili peppers (Capsicum spp.). Int. J. Mol. Sci. 2013, 14, 19025–19053. [Google Scholar] [CrossRef] [Green Version]
  31. Lu, S.; Li, L. Carotenoid metabolism: Biosynthesis, regulation, and beyond. J. Integr. Plant Biol. 2008, 50, 778–785. [Google Scholar] [CrossRef] [PubMed]
  32. Giuliano, G. Plant carotenoids: Genomics meets multi-gene engineering. Curr. Opin. Plant Biol. 2014, 19, 111–117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Fantini, E.; Falcone, G.; Frusciante, S.; Giliberto, L.; Giuliano, G. Dissection of tomato lycopene biosynthesis through virus-induced gene silencing. Plant Physiol. 2013, 163, 986–998. [Google Scholar] [CrossRef]
  34. Kachanovsky, D.E.; Filler, S.; Isaacson, T.; Hirschberg, J. Epistasis in tomato color mutations involves regulation of phytoene synthase 1 expression by cis-carotenoids. Proc. Natl. Acad. Sci. USA 2012, 109, 19021–19026. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Fraser, P.D.; Römer, S.; Shipton, C.A.; Mills, P.B.; Kiano, J.W.; Misawa, N.; Drake, R.G.; Schuch, W.; Bramley, P.M. Evaluation of transgenic tomato plants expressing an additional phytoene synthase in a fruit-specifific manner. Proc. Natl. Acad. Sci. USA 2002, 99, 1092–1097. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Wu, M.; Xu, X.; Hu, X.; Liu, Y.; Cao, H.; Chan, H.; Gong, Z.; Yuan, Y.; Luo, Y.; Feng, B.; et al. SlMYB72 regulates the metabolism of chlorophylls, carotenoids, and flavonoids in tomato fruit. Plant Physiol. 2020, 183, 854–868. [Google Scholar] [CrossRef] [PubMed]
  37. Xiong, C.; Luo, D.; Lin, A.; Zhang, C.; Shan, L.; He, P.; Li, B.; Zhang, Q.; Hua, B.; Yuan, Z.; et al. A tomato B-box protein SlBBX20 modulates carotenoid biosynthesis by directly activating Phytoene Synthase 1, and is targeted for 26S proteasome-mediated degradation. New Phytol. 2019, 221, 279–294. [Google Scholar] [CrossRef] [Green Version]
  38. Zhu, M.; Chen, G.; Zhou, S.; Tu, Y.; Wang, Y.; Dong, T.; Hu, Z. A new tomato NAC (NAM/ATAF1/2/CUC2) transcription factor, SlNAC4, functions as a positive regulator of fruit ripening and carotenoid accumulation. Plant Cell Physiol. 2014, 55, 119–135. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Li, L.; Wang, X.; Zhang, X.; Guo, M.; Liu, T. Unraveling the target genes of RIN transcription factor during tomato fruit ripening and softening. J. Sci. Food Agric. 2017, 97, 991–1000. [Google Scholar] [CrossRef] [PubMed]
  40. Lee, J.M.; Joung, J.G.; McQuinn, R.; Chung, M.Y.; Fei, Z.; Tieman, D.; Klee, H.; Giovannoni, J. Combined transcriptome, genetic diversity and metabolite profiling in tomato fruit reveals that the ethylene response factor SlERF6 plays an important role in ripening and carotenoid accumulation. Plant J. 2012, 70, 191–204. [Google Scholar] [CrossRef] [PubMed]
  41. Robertson, G.H.; Mahoney, N.E.; Goodman, N.; Pavlath, A.E. Regulation of lycopene formation in cell suspension culture of VFNT tomato (Lycopersicon esculentum) by CPTA, growth regulators, sucrose, and temperature. J. Exp. Bot. 1995, 46, 667–673. [Google Scholar] [CrossRef] [Green Version]
  42. Zhang, Z.; Liu, L.; Zhang, M.; Zhang, Y.; Wang, Q. Effect of carbon dioxide enrichment on health-promoting compounds and organoleptic properties of tomato fruits grown in greenhouse. Food Chem. 2014, 153, 157–163. [Google Scholar] [CrossRef]
  43. Yang, Z. PAML 4: Phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 2007, 24, 1586–1591. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Liu, S.; Liu, Y.; Yang, X.; Tong, C.; Edwards, D.; Parkin, I.A.; Zhao, M.; Ma, J.; Yu, J.; Huang, S.; et al. The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes. Nat. Commun. 2014, 5, 3930. [Google Scholar] [CrossRef] [PubMed]
  45. Sparkes, I.A.; Runions, J.; Kearns, A.; Hawes, C. Rapid, transient expression of fluorescent fusion proteins in tobacco plants and generation of stably transformed plants. Nat. Protoc. 2006, 1, 2019–2025. [Google Scholar] [CrossRef] [PubMed]
  46. Hittinger, C.T.; Carroll, S.B. Gene duplication and the adaptive evolution of a classic genetic switch. Nature 2007, 449, 677–681. [Google Scholar] [CrossRef]
  47. Spillane, C.; Schmid, K.J.; Laoueille-Duprat, S.; Pien, S.; Escobar-Restrepo, J.M.; Baroux, C.; Gagliardini, V.; Page, D.R.; Wolfe, K.H.; Grossniklaus, U. Positive darwinian selection at the imprinted MEDEA locus in plants. Nature 2007, 448, 349–352. [Google Scholar] [CrossRef]
  48. Town, C.D.; Cheung, F.; Maiti, R.; Crabtree, J.; Haas, B.J.; Wortman, J.R.; Hine, E.E.; Althoff, R.; Arbogast, T.S.; Tallon, L.J.; et al. Comparative genomics of Brassica oleracea and Arabidopsis thaliana reveal gene loss, fragmentation, and dispersal after polyploidy. Plant Cell 2006, 18, 1348–1359. [Google Scholar] [CrossRef] [Green Version]
  49. Lysak, M.A.; Koch, M.A.; Pecinka, A.; Schubert, I. Chromosome triplication found across the tribe Brassiceae. Genome Res. 2005, 15, 516–525. [Google Scholar] [CrossRef] [Green Version]
  50. Gao, H.Y.; Zhu, B.Z.; Zhu, H.L.; Zhang, Y.L.; Xie, Y.H.; Li, Y.C.; Luo, Y.B. Effect of suppression of ethylene biosynthesis on flavor products in tomato fruits. Russ. J. Plant Physiol. 2007, 54, 80–88. [Google Scholar] [CrossRef]
  51. Sun, L.; Yuan, B.; Zhang, M.; Wang, L.; Cui, M.; Wang, Q.; Leng, P. Fruit-specific RNAi-mediated suppression of SlNCED1 increases both lycopene and β-carotene contents in tomato fruit. J. Exp. Bot. 2012, 63, 3097–3108. [Google Scholar] [CrossRef]
  52. McQuinn, R.P.; Gapper, N.E.; Gray, A.G.; Zhong, S.; Tohge, T.; Fei, Z.; Fernie, A.R.; Giovannoni, J.J. Manipulation of ZDS in tomato exposes carotenoid- and ABA-specific effects on fruit development and ripening. Plant Biotechnol. J. 2020, 18, 2210–2224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Zhu, Y.N.; Yang, Q.Q.; Kou, E.F.; Song, S.W.; Chen, R.Y.; Hao, Y.W. Effects of salicylic acid treatment on growth and quality of flower Chinese cabbage. North Hortic. 2020, 22, 53–59. [Google Scholar] [CrossRef]
  54. Cunningham, F.X.; Gantt, E. Genes and enzymes of carotenoid biosynthesis in plants. Annu. Rev. Plant Physiol. Plant Mol. Biol. 1998, 49, 557–583. [Google Scholar] [CrossRef] [PubMed]
  55. Welsch, R.; Arango, J.; Bar, C.; Salazar, B.; Al-Babili, S.; Beltran, J.; Chavarriaga, P.; Ceballos, H.; Tohme, J.; Beyer, P. Provitamin A accumulation in cassava (Manihot esculenta) roots driven by a single nucleotide polymorphism in a phytoene synthase gene. Plant Cell 2010, 22, 3348–3356. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Maass, D.; Arango, J.; Wust, F.; Beyer, P.; Welsch, R. Carotenoid crystal formation in Arabidopsis and carrot roots caused by increased phytoene synthase protein levels. PLoS ONE 2009, 4, e6373. [Google Scholar] [CrossRef] [Green Version]
  57. Shewmaker, C.; Sheehy, J.; Daley, M.; Colburn, S.; Ke, D. Seed-specifific overexpression of phytoene synthase: Increase in carotenoids and other metabolic effects. Plant J. 1999, 20, 401–412. [Google Scholar] [CrossRef]
  58. Li, P.; Zhang, S.; Zhang, S.; Li, F.; Zhang, H.; Cheng, F.; Wu, J.; Wang, X.; Sun, R. Carotenoid biosynthetic genes in Brassica rapa: Comparative genomic analysis, phylogenetic analysis, and expression profiling. BMC Genom. 2015, 16, 492. [Google Scholar] [CrossRef] [Green Version]
  59. Fu, Z.; Chai, Y.; Zhou, Y.; Yang, X.; Warburton, M.L.; Xu, S.; Cai, Y.; Zhang, D.; Li, J.; Yan, J. Natural variation in the sequence of PSY1 and frequency of favorable polymorphisms among tropical and temperate maize germplasm. Theor. Appl. Genet. 2013, 126, 923–935. [Google Scholar] [CrossRef]
Genes 12 02027 g001 550
Figure 1. Circos diagram of syntenic carotenoid biosynthetic genes in B. oleracea, B. rapa, and A. thaliana. C01 to C09 indicate B. oleracea chromosomes, A01 to A10 indicate B. rapa chromosomes, and Chr1 to Chr5 indicate A. thaliana chromosomes.
Figure 1. Circos diagram of syntenic carotenoid biosynthetic genes in B. oleracea, B. rapa, and A. thaliana. C01 to C09 indicate B. oleracea chromosomes, A01 to A10 indicate B. rapa chromosomes, and Chr1 to Chr5 indicate A. thaliana chromosomes.
Genes 12 02027 g001
Genes 12 02027 g002 550
Figure 2. Genomic localization of the carotenoid biosynthetic genes on the nine chromosomes of B. oleracea. C01 to C09 indicate B. oleracea chromosomes.
Figure 2. Genomic localization of the carotenoid biosynthetic genes on the nine chromosomes of B. oleracea. C01 to C09 indicate B. oleracea chromosomes.
Genes 12 02027 g002
Genes 12 02027 g003 550
Figure 3. Ka/Ks values of orthologous gene pairs for carotenoid biosynthetic genes among A. thaliana, B. oleracea and B. rapa.
Figure 3. Ka/Ks values of orthologous gene pairs for carotenoid biosynthetic genes among A. thaliana, B. oleracea and B. rapa.
Genes 12 02027 g003
Genes 12 02027 g004 550
Figure 4. Heat map representation of orthologous gene pairs of carotenoid biosynthetic genes of B. oleracea and B. rapa. Color scale bar represents log2 transformed FPKM values.
Figure 4. Heat map representation of orthologous gene pairs of carotenoid biosynthetic genes of B. oleracea and B. rapa. Color scale bar represents log2 transformed FPKM values.
Genes 12 02027 g004
Genes 12 02027 g005 550
Figure 5. Relative expression levels of nine key carotenoid biosynthetic genes in different tissues of cabbage.
Figure 5. Relative expression levels of nine key carotenoid biosynthetic genes in different tissues of cabbage.
Genes 12 02027 g005
Genes 12 02027 g006 550
Figure 6. The expression of BoPSY.1, BoPDS.1, BoZDS, and BoLYC after treatment with exogenous phytohormone. The Y-axis and X-axis represent the relative expression level and the time course of exogenous phytohormone treatment, respectively. Leaves were sampled at 0, 2, 4, 6, 12, and 24 h after ETH, ABA, SA, and MeJA spraying. Data represent the mean ± SD of three technical repetitions.
Figure 6. The expression of BoPSY.1, BoPDS.1, BoZDS, and BoLYC after treatment with exogenous phytohormone. The Y-axis and X-axis represent the relative expression level and the time course of exogenous phytohormone treatment, respectively. Leaves were sampled at 0, 2, 4, 6, 12, and 24 h after ETH, ABA, SA, and MeJA spraying. Data represent the mean ± SD of three technical repetitions.
Genes 12 02027 g006
Genes 12 02027 g007 550
Figure 7. Phylogenetic relationship and conserved motif analysis of PSY genes. The phylogenetic tree was generated using the Muscle module and neighbor-joining (NJ) method with 1000 bootstrap replicates implemented in MEGA 6.0. Sl: Solanum lycopersicum; Ca: Capsicum annuum; Lb: Lycium barbarum; Ib: Ipomoea batatas; Gj: Gardenia jasminoides; Of: Osmanthus fragrans; Pc: Pogostemon cablin; Dk: Diospyros kaki; Ad: Actinidia deliciosa; Cl: Citrullus lanatus; Cm: Cucurbita moschata; Cp: Carica papaya; Bo: Brassica oleracea; At: Arabidopsis thaliala; Nt: Narcissus tazetta; Ll: Lilium lancifolium; Ms: Medicago sativa; Dc: Daucus carota; Fa: Fragaria ananassa; Md: Malus domestica; Os: Oryza sativa Indica; Zm: Zea mays; Sb: Sorghum bicolor; Hl: Haematococcus lacustris.
Figure 7. Phylogenetic relationship and conserved motif analysis of PSY genes. The phylogenetic tree was generated using the Muscle module and neighbor-joining (NJ) method with 1000 bootstrap replicates implemented in MEGA 6.0. Sl: Solanum lycopersicum; Ca: Capsicum annuum; Lb: Lycium barbarum; Ib: Ipomoea batatas; Gj: Gardenia jasminoides; Of: Osmanthus fragrans; Pc: Pogostemon cablin; Dk: Diospyros kaki; Ad: Actinidia deliciosa; Cl: Citrullus lanatus; Cm: Cucurbita moschata; Cp: Carica papaya; Bo: Brassica oleracea; At: Arabidopsis thaliala; Nt: Narcissus tazetta; Ll: Lilium lancifolium; Ms: Medicago sativa; Dc: Daucus carota; Fa: Fragaria ananassa; Md: Malus domestica; Os: Oryza sativa Indica; Zm: Zea mays; Sb: Sorghum bicolor; Hl: Haematococcus lacustris.
Genes 12 02027 g007
Genes 12 02027 g008 550
Figure 8. Subcellular localization of BoPSY proteins in tobacco leaf epidermis cells. (a–d) Subcellular localization of BoPSY.1. (e–h) Subcellular localization of BoPSY.2. (i–l) Subcellular localization of BoPSY.3. (a,e,i) Fluorescence signals of BoPSYs fused with GFP. (b,f,j) Fluorescence signal of chloroplast. (c,g,k) Bright field images. (d,h,l) Superposition images of bright field and fluorescence.
Figure 8. Subcellular localization of BoPSY proteins in tobacco leaf epidermis cells. (a–d) Subcellular localization of BoPSY.1. (e–h) Subcellular localization of BoPSY.2. (i–l) Subcellular localization of BoPSY.3. (a,e,i) Fluorescence signals of BoPSYs fused with GFP. (b,f,j) Fluorescence signal of chloroplast. (c,g,k) Bright field images. (d,h,l) Superposition images of bright field and fluorescence.
Genes 12 02027 g008
Table
Table 1. Carotenoid biosynthetic genes identified in cabbage.
Table 1. Carotenoid biosynthetic genes identified in cabbage.
EnzymeArabidopsis thalianaB. olerecea
Syntenic OrthologsNon-Syntenic Orthologs
LFMF1MF2
DXSAT4G15560BoDXS.1 (Bol020307)BoDXS.2 (Bol005061)-BoDXS.3 (Bol005062)
DXRAT5G62790BoDXR.1 (Bol019430)-BoDXR.2 (Bol019181)-
MCTAT2G02500-BoMCT (Bol015155)--
CMKAT2G26930BoCMK (Bol027780)---
MDSAT1G63970BoMDS.1 (Bol021045)-BoMDS.2 (Bol022634)-
HDSAT5G60600BoHDS (Bol003103)---
HDRAT4G34350BoHDR.1 (Bol013657)-BoHDR.2 (Bol001484)-
IPPI1AT5G16440-BoIPPI1 (Bol034402)--
IPPI2AT3G02780BoIPPI2.1 (Bol000950)BoIPPI2.2 (Bol001738)BoIPPI2.3 (Bol002232)-
GGPS1AT4G36810BoGGPS1 (Bol028967)---
GGPS2AT2G23800BoGGPS2 (Bol045796)---
GGPS3AT3G14550BoGGPS3.1 (Bol005099)-BoGGPS3.2 (Bol025714)-
GGPS7AT2G18620----
GGPS8AT3G14510----
GGPS11AT3G29430----
GGPS4AT2G18640----
GGPS6AT1G49530----
GGPS9AT3G14530----
GGPS10AT3G20160----
GGPS12AT3G32040----
GGRAT4G38460BoGGR (Bol000792)---
PSYAT5G17230BoPSY.1 (Bol019820)BoPSY.2 (Bol034439)BoPSY.3 (Bol021326)-
PDS3AT4G14210-BoPDS3.1 (Bol009962)BoPDS3.2 (Bol016089)-
Z-ISOAT1G10830BoZ-ISO (Bol036691)---
ZDSAT3G04870-BoZDS (Bol003131)--
CRTISOAT1G06820--BoCRTISO (Bol004146)-
CRTISO2AT1G57770--BoCRTISO2 (Bol016311)-
LYCAT3G10230BoLYC (Bol011368)---
LUT2AT5G57030BoLUT2.1 (Bol037991)BoLUT2.2 (Bol026053)--
CHY1AT4G25700BoCHY1.1 (Bol039555)BoCHY1.2 (Bol042254)--
CHY2AT5G52570BoCHY2.1 (Bol030235)-BoCHY2.2 (Bol045243)-
LUT5AT1G31800-BoLUT5.1 (Bol027064)-BoLUT5.2 (Bol027080)
BoLUT5.3 (Bol014196)
CYP97B3AT4G15110--BoCYP97B3 (Bol028312)-
LUT1AT3G53130----
ZEPAT5G67030BoZEP.1 (Bol027222)-BoZEP.2 (Bol019241)-
VDEAT1G08550BoVDE (Bol041221)---
NSYAT1G67080-BoNSY (Bol045096)--
CCD7AT2G44990-BoCCD7 (Bol021707)--
CCD8AT4G32810BoCCD8.1 (Bol017864)--BoCCD8.2 (Bol000605)
NCED2AT4G18350BoNCED2.1 (Bol009433)BoNCED2.2 (Bol037062)--
NCED3AT3G14440BoNCED3.1 (Bol005093)BoNCED3.2 (Bol011830)BoNCED3.3 (Bol035582)-
NCED4AT4G19170BoNCED4.1 (Bol009345)-BoNCED4.2 (Bol029878)-
NCED5AT1G30100BoNCED5 (Bol022516)---
NCED6AT3G24220BoNCED6 (Bol007451)---
NCED9AT1G78390BoNCED9.1 (Bol027485)BoNCED9.2 (Bol018961)--
ABA2AT1G52340BoABA2.1 (Bol000478)BoABA2.2 (Bol035060)--
AAO3AT2G27150-BoAAO3.1 (Bol026459)--
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cao, W.; Wang, P.; Yang, L.; Fang, Z.; Zhang, Y.; Zhuang, M.; Lv, H.; Wang, Y.; Ji, J. Carotenoid Biosynthetic Genes in Cabbage: Genome-Wide Identification, Evolution, and Expression Analysis. Genes 2021, 12, 2027. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12122027

AMA Style

Cao W, Wang P, Yang L, Fang Z, Zhang Y, Zhuang M, Lv H, Wang Y, Ji J. Carotenoid Biosynthetic Genes in Cabbage: Genome-Wide Identification, Evolution, and Expression Analysis. Genes. 2021; 12(12):2027. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12122027

Chicago/Turabian Style

Cao, Wenxue, Peng Wang, Limei Yang, Zhiyuan Fang, Yangyong Zhang, Mu Zhuang, Honghao Lv, Yong Wang, and Jialei Ji. 2021. "Carotenoid Biosynthetic Genes in Cabbage: Genome-Wide Identification, Evolution, and Expression Analysis" Genes 12, no. 12: 2027. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12122027

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop