Next Article in Journal
A Potent Antifungal Agent for Basal Stem Rot Disease Treatment in Oil Palms Based on Chitosan-Dazomet Nanoparticles
Previous Article in Journal
What Drives Symbiotic Calcium Signalling in Legumes? Insights and Challenges of Imaging
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 9
10.3390/ijms20092246
ijms-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

Illumina® Sequencing Reveals Candidate Genes of Carotenoid Metabolism in Three Pummelo Cultivars (Citrus Maxima) with Different Pulp Color

by
Cui-Cui Jiang
1,
Yan-Fang Zhang
1,
Yan-Jin Lin
1,
Yuan Chen
2 and
Xin-Kun Lu
1,*
1
Fruit Research Institute, Fujian Academy of Agricultural Sciences, Fuzhou 350013, China
2
Institute of Agricultural Engineering and Technology, Fujian Academy of Agricultural Sciences, Fuzhou 350003, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(9), 2246; https://doi.org/10.3390/ijms20092246
Submission received: 10 February 2019 / Revised: 27 April 2019 / Accepted: 5 May 2019 / Published: 7 May 2019
(This article belongs to the Section Molecular Plant Sciences)
Browse Figures
Versions Notes

Abstract

:
Pummelo (Citrus maxima) is one of important fruit trees, which belongs to Citrus species. The fruits of different pummelo cultivars have different colors and differ in the contents of carotenoid. Our results clearly showed that ‘Huangjinmiyou’ (HJMY) has the highest content of β-carotene, followed by ‘Hongroumiyou’ (HRMY) and ‘Guanximiyou’ (GXMY). Lycopene is dominantly accumulated in HRMY. However, the molecular mechanism underlying the carotenoid accumulation in pummelo flesh is not fully understood. In this study, we used the RNA-Seq technique to investigate the candidate genes of carotenoid metabolism in the flesh of pummelo cv. GXMY and its mutants HRMY and HJMY in three development periods of fruit. After data assembly and bioinformatic analysis, a total of 357 genes involved in biosynthesis of secondary metabolites were isolated, of which 12 differentially expressed genes (DEGs) are involved in carotenoid biosynthesis. Among these 12 DEGs, phytoene synthase (PSY2), lycopene β-cyclase (LYCB2), lycopene Ɛ-cyclase (LYCE), carotenoid cleavage dioxygenases (CCD4), 9-cis-epoxycarotenoid dioxygenase (NCED2), aldehyde oxidase 3 (AAO3), and ABA 8′-hydroxylases (CYP707A1) are the most distinct DEGs in three pummelo cultivars. The co-expression analysis revealed that the expression patterns of several transcription factors such as bHLH, MYB, ERF, NAC and WRKY are highly correlated with DEGs, which are involved in carotenoid biosynthesis. In addition, the expression patterns of 22 DEGs were validated by real-time quantitative PCR (RT-qPCR) and the results are highly concordant with the RNA-Seq results. Our results provide a global vision of transcriptomic profile among three pummelo cultivars with different pulp colors. These results would be beneficial to further study the molecular mechanism of carotenoid accumulation in pummelo flesh and help the breeding of citrus with high carotenoid content.

Graphical Abstract

1. Introduction

Carotenoid is the second most abundant compound in natural pigments and provides colorless, yellow, orange, and red colors for many fruits, flowers and vegetables [1]. Plant is the main source of carotenoid for human. In recent years, with the improvement of people’s living standards and health consciousness, horticultural products for high carotenoid content have been favored by many consumers. Therefore, it is critical and urgent to accelerate horticultural plant breeding for high carotenoid content through breeding technologies.
Citrus takes the leading position in fruit industry worldwide. Color of citrus flesh represents an essential fruit quality attribute. Recently, high demand for citrus with different flesh colors provides potentials for genetic improvement of citrus. It is clear that coloration of most citrus fruits is mainly caused by the accumulation of carotenoid [2,3,4]. Composition and content of carotenoid in fruits depend on cultivars and species [5]. In order to regulate coloration of citrus fruit genetically, several novel methods, such as gene cloning and transcriptomic and proteomic analysis, are used to isolate genes associated with carotenoid metabolism [2,6]. Until now, RNA-Seq has been successfully applied in fruit development of Citrus [7,8,9,10]. Previous studies mostly focused on the content and composition of carotenoid in sweet orange (C. sinensis L. Osbeck), satsuma mandarin (C. unshiu Marc.), and lemon (C. limon Burm.) [2]. However, the carotenoid regulation mechanism in pummelo flesh was not well understood. Pummelo is monoembryonic, and for this reason, its true hybrids are easily obtained by cross-breeding without nucellar seedling in progenies. It is helpful to identify genetic control of carotenoid metabolism with pummelo as materials. Moreover, pummelo, one of three basic species of citrus, shows lower heterozygosity than sweet orange, so it can overcome heterozygous loci in genetic analysis. Liu et al. [11] investigated the carotenoid profile of red-flesh pummelo cv. “Chuhong” and pale green-flesh pummelo cv. “Feicui” by high-performance liquid chromatography (HPLC) and GC/MS and concluded that there are different tissue-biased patterns of carotenoid biosynthesis in pummelos with different flesh colors. Guo et al. [10] identified the differentially expressed genes in GXMY and its orange-pericarp mutant using the RNA-Seq technique and revealed extensive changes of gene expression level between the control and its mutant. There was no report on the molecular mechanism of carotenoid accumulation in pummelo flesh.
Since 1998, we have discovered two mutants with different flesh colors from GXMY in Pinghe County, Zhangzhou City, China. These mutants can be used as materials to study the molecular mechanism of carotenoid accumulation. C. maxima cv. GXMY is traditionally famous fruit in China and has yellow flesh, in which the total carotenoid content is 0.875 μg/g in fresh weight (FW) [12]. HRMY is a bud mutant of GXMY with red flesh, of which the total carotenoid content is 96.55 μg/g FW [12]. HJMY is also a high-quality bud mutant of GXMY with golden yellow flesh, of which the total carotenoid content is 251.77 μg/g FW [12,13]. In order to unveil the molecular mechanism underlying different colorations of GXMY and its mutants HJMY and HRMY, we investigated the gene expression and transcription factors associated with carotenoid metabolism in the flesh of pummelo cv. GXMY, HJMY and HRMY in three key stages of fruit development (according to the pigment contents and the ratio of pulp and fruit) by using the RNA-Seq technique. The identified genes could be useful for studying the molecular mechanism of carotenoid accumulation in citrus flesh and help the breeding of citrus with high carotenoid content.

2. Results

2.1. Changes in β-Carotene and Lycopene Contents During Fruit Development of Three Different Pummelo Cultivars

The appearances of three pummelo cultivars of different fruit growth and development periods are presented in Figure 1. The contents of β-carotene and lycopene in three different pummelo cultivars were measured by HPLC (Figure 2). As results indicated, HJMY has the highest content of β-carotene with a value of about 681.45 mg/kg FW, almost 34 times higher than that in GXMY and HRMY at an early development stage 40 days after flowering (DAF). The β-carotene content decreases at 100 DAF and then increases at 190 DAF in HJMY and HRMY, whereas it decreases during the entire development of GXMY. Lycopene was detected only in the HRMY, showing an increasing trend with the development of the fruit and mainly accumulated in the late stage of fruit development. Based on the overall trends, HJMY has the highest content of β-carotene, followed by HRMY and GXMY, and lycopene is dominantly accumulated in HRMY.

2.2. Transcriptome Sequence and De Novo Assembly

Nine cDNA libraries were constructed from the total RNA of three pummelo cultivars at 40 DAF, 100 DAF, and 190 DAF. The original sequencing data were filtered and assembled as shown in Table 1. The results showed that 201.78 million high-quality reads were obtained with a total base number of 40.66 billion (giga or G) and a total of 17.65 million contigs were assembled. The Q30 value (base call accuracy of each read is more than 99.9%) of clean reads in each library is more than 85%. These libraries were subjected to RNA-Seq using an Illumina HiSeq2500 and the clean reads from different libraries are presented in Table 1. Using Trinity software, the clean reads from the nine libraries were assembled into 106,957 transcripts with an average length of 858 bp and 54,051 unigenes with a mean length of 858 bp (Table 1). The size distributions of these unigenes are shown in Figure 3. Among the unigenes, the lengths of 12,902 unigenes (23.87% of the total number of unigenes of 54,051) range from 300 to 500 bp, the lengths of 10,547 unigenes (19.51% of the total number of unigenes) range from 501 to 1000 bp, the lengths of 8116 unigenes (15.02% of the total number of unigenes) range from 1000 to 2000 bp, and the lengths of 5953 unigenes (11.01% of the total number of unigenes) are more than 2000 bp.

2.3. Functional Annotation and Classification of Unigenes

To determine and annotate the transcriptome data of three pummelo cultivars, approximately 24,206 unigenes were searched against five public databases [NCBI non-redundant protein (Nr), Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups (COG), and Gene Ontology (GO)] with an E-value of ≤ 10−5 by BLASTX and BLASTN searching engines [14]. The functional annotation results showed that there were 8876, 17,764, 6093, 16,857, and 24,125 unigenes annotated to COG, GO, KEGG, Swiss-Prot, and Nr, respectively. Totally, only 44.8% of the unigenes (24,206) could be successfully annotated to these five public databases. The remaining unigenes (55.2% of the total number of unigenes) have no homologs with the known genes, which is likely due to the presence of short sequences. Pathway-based analysis can help to understand the biological functions and interactions of genes. There were 6093 unigenes assigned to 117 KEGG pathways (24,206) (Table S1, Supplementary Materials). The majority of unigenes were associated with metabolic pathways (1847). We focused on the “Biosynthesis of secondary metabolites” category in relation to fruit pigmentation. In this category, 357 unigenes were classified into nine subcategories (Figure 4). Among these nine subcategories, “Carotenoid biosynthesis” is the highest clustered, followed by “Flavonoid biosynthesis”. “Flavone and flavonol biosynthesis” and “Caffeine metabolism” appear to be the smallest groups.

2.4. Analysis of DEG During Fruit Development of Three Different Pummelo Cultivars

Using the criteria of p < 0.01 and |log2 (fold change)| ≥ 2, a differential expression analysis of unigenes was performed during fruit development of three different pummelo cultivars. The relationships between different DEG groups are displayed as Venn diagram (Figure 5), and 616, 258, and 256 unigenes were shared in all development stages of three pummelo cultivars, respectively. Through comparing the libraries, we found a number of differentially related genes, as shown in Table 2. The results suggested that the developmental periods with the most dynamic changes in the transcriptome are between 40 DAF and 190 DAF in GXMY, and between 100 DAF and 190 DAF in HRMY and HJMY.

2.5. Candidate Genes Involved in Carotenoid Biosynthesis

To further search for the genes involved in carotenoid biosynthesis, we analyzed the transcriptomic data of three pummelo cultivars at different developmental stages. Based on the KEGG annotation, 12 DEGs that encode enzymes related to carotenoid biosynthesis pathway were isolated (Table 3 and Table S2 (Supplementary Materials); Figure 6).
Transcription factors (TFs) are key regulators of gene expression associated with various aspects of fruit growth and development including carotenoid biosynthesis. In order to identify differentially expressed transcription factors from three pummelo cultivars, a BLASTX search was conducted using all the transcription factors sequences available in the PlnTFDB database. In total, 435 unigenes (Table S3, Supplementary Materials) were predicted to encode transcription factors from 53 different families (Table S4, Supplementary Materials) and 227 of them were differentially expressed (Table S5, Supplementary Materials). A transcription abundance correlation analysis was carried out between the differentially expressed transcription factors and structural genes from the carotenoid biosynthesis pathway. The expression levels of the 48 transcription factors are highly correlated with those of the candidate structural genes (Table S6, Supplementary Materials). The identified transcription factors included the transcription factors that are implicated in regulating carotenoid biosynthesis, such as bHLH, ERF, and MYB (Table 4). The 34 bHLH were differentially expressed during fruit development stages of three pummelo cultivars and nine of them [c11225.graph_c0 (homologous to CsbHLH107), c12645.graph_c0 (homologous to BIM1 isoform X1), c21405.graph_c0 (homologous to CsbHLH162), c22086.graph_c0 (homologous to CsbHLH57), c2272.graph_c0 (homologous to CsbHLH35), c22910.graph_c0 (homologous to CsbHLH61), c27501.graph_c0 (homologous to CsbHLH130), c6149.graph_c0 (homologous to ILR3), and c589.graph_c1 (homologous to HEC2)] were associated with the carotenoid biosynthesis pathway (Table 4, Tables S3 and S6 (Supplementary Materials)). Eight of the differentially expressed transcription factors annotated as ERF or AP2/ERF showed a significant positive correlation with carotenoid biosynthesis-related genes (Table 4 and Table S6 (Supplementary Materials)). Although they are correlated with none of the structural genes, the expression levels of the homolog of CsARF19 (c23948.graph_c1) and CsARF5 (c26184.graph_c0) show a differential expression during fruit development of three different cultivars. Four of the differentially expressed transcription factors annotated as MYB show a significant positive correlation with carotenoid biosynthesis-related genes (Table 4 and Table S6 (Supplementary Materials)), while only two of the MYB genes (c22501.graph_c0 and c23847.graph_c0) show a negative correlation with the expression of structural genes. A MYB gene (c21421.graph_c0), which shares a high similarity with transcription factor AS1 (XP_006429090.1) from C. clementina, significantly correlates with six carotenoid biosynthesis-related structural genes (Table S6, Supplementary Materials). In addition, one NAC (c13516.graph_c0) and two WRKYs (c21980.graph_c0 and c34311.graph_c0), which share a high sequence similarity with the NAC domain-containing protein 72 (XM_006451877.2), WRKY22 (XM_0064444887.2), and WRKY17 (XM_006450132.2) from C. clementina, respectively, were also found to be correlated to the genes involved in the carotenoid pathway.

2.6. Real-Time Quantitative PCR (RT-qPCR) Analysis of Selected DEGs During Fruit Development of Three Pummelo Cultivars

To further confirm the results of transcriptomic analysis, 11 enzyme-encoding genes and twelve transcription factors associated with carotenoid biosynthesis were analyzed by RT-qPCR in three pummelo cultivars. As RT-qPCR indicated, all selected DEGs were expressed at varying levels during the three stages and their RNA-Seq data have a close correlation (Figure 7; Figure 8).

3. Discussion

Carotenoids are the main pigments in citrus fruits and different carotenoid components determined different colors in citrus fruits. Carotenoid has an important role in the sensory quality and nutritional quality of citrus fruit. Our previous studies have shown that the colorful pigments of HJMY juice are mainly due to the contents of β-carotene, phytoene, α-carotene, and β-cryptoxanthin, and the total amount of carotenoid reaches 251.77 μg·g−1·FW, which are 2.6 and 287.7 times those of HRMY and GXMY, respectively [12]. This study showed that the β-carotene and lycopene are dominantly accumulated in HJMY and HRMY during the development of pummelo fruits, respectively, which confer the golden orange color in HJMY and red color in HRMY (Figure 1 and Figure 2).
Identification and profile of the members of gene families involved in the carotenoid pathway could explain the diversity of coloration in fruit. Previous studies have revealed that the key regulatory genes, such as phytoene synthase (PSY), lycopene β-cyclase (LCYB), lycopene ε-cyclase (LCYE), CCD and NCED in carotenoid biosynthesis, might be key factors in controlling the carotenoid metabolism (Figure 6) [15,16,17,18,19,20,21]. The PSY is the rate-limiting enzyme in carotenoid biosynthesis and plays an important role in the biosynthesis of carotenoid in many plants, such as apple, carrot, and tomato [22,23,24]. PSY appears to be a small gene family in most plants including citrus. The up-regulation of PSY1 gene and accumulation of carotenoid are found in the flavedo of sweet orange [25]. However, there are few reports regarding PSY2. Recently, a new PSY member (PSY2) was reported in sweet orange and its expression is at a lower level than that of PSY1 in both flavedo and juice sacs in Valencia and Rohde Red Valencia [26], which suggests that PSY2 may also contribute to the accumulation of carotenoid in citrus fruit. In the current study, the PSY2 (c27977.graph_c0, Table 3; Figure 7A) has a high expression level at an earlier stage of GXMY and then decreases during the entire fruit development. It has relatively low expression levels in both HRMY and HJMY. According to our pigment measurements, PSY2 might not be the key regulatory gene responsible for carotenoid biosynthesis and accumulation in pummelo pulp (Figure 2 and Figure 7A). However, the PSY2 has the highest expression level in yellow apple germplasm, which is significantly positively correlated with carotenoid contents in apple [22].
Previous studies showed that the transcriptional regulation of LYCB and LYCE are the critical regulatory points in carotenoid biosynthesis [27]. In citrus, two LCYB members have also been isolated in navel orange fruit and grapefruit cultivars, “marsh” (white flesh) and “flame” (red flesh) [28]. Two LCYE genes are identified in sweet orange tissues, but only one LCYE gene is detected in leaves and fruits [29]. The expression levels of VvLCYB and VvLCYE gradually increase with the accumulation of carotenoid contents during the ripening of grapevine fruit [18]. In sweet potatoes, over-expression of IbLCYB2 significantly increases the carotenoid contents and enhances abiotic tolerance by positively regulating the biosynthesis of carotenoid and ABA [30]. The expression level of LCYB2 increases markedly during the ripening of papaya [31]. Here, our results showed that the expression level of LYCB2 increases from 40 DAF to ripening stages (190 DAF) in all the three pummelo cultivars with a slightly lower expression level in HJMY at the ripening stage (190 DAF) (Table 3; Figure 7B). The down-regulation of LYCE is only found at the development of GXMY, whereas the expression level of LYCE is firstly up-regulated and then up-regulated during fruit development in HRMY and GXMY (Table 3; Figure 7C). The expression levels of LYCE are higher in HRMY and HJMY at the color break period (100 DAF) and ripening period (190 DAF) than in GXMY, which may contribute to the high α-carotene content in HRMY and HJMY pulp [12].
In plants, the CCD family consists of nine members, which catalyze the oxidative cleavage of carotenoid [32]. Among them, five genes are related to the synthesis of ABA, named NCED (NCED2, 3, 5, 6, 9), and the remaining four are named CCD1, 4, 7, 8, which are closely related to the growth and development of lateral branches, plant hormones (such as strigolactones), the formation of volatile aromatic compounds, and colorless or red volatile apocarotenoids compounds [33]. Carotenoid accumulation in several yellow fruits and flowers has been shown to be negatively associated with CCD1 or CCD4 expression. For example, transcription inactivation of CmCCD4a causes yellow coloration of flowers in chrysanthemum [34]. In azalea petals, the expression levels of CCD4 are higher in ‘Miyamasatsuki’ (white-flowered) and the progeny than in R. japonicum f. flavum (yellow-flowered) at all stages [35]. In fruit, the negative correlation between the expression level of CCD4 and the carotenoid content is also observed in summer squash, peach, potato, and apple [36,37,38,39]. Here, we also found that the expression level of CCD4 (c20564.graph_c0) gene, which is homologous to the CCD4a in C. clementina (DQ309330.1), decreases with the fruit development in three pomelo varieties (Table 3; Figure 7E). In addition, we also found that the expression level of the CCD1 (c29184.graph_c0) increases from the early stages of fruit development (40 DAF) to the color break period (100 DAF), followed by a decrease as fruit becomes mature in HJMY and HRMY. In GXMY, the expression level of CCD1 decreases as the fruit becomes mature (Table 3; Figure 7D).
ABA is synthesized in the plastids from xanthophyll pigments and is derived from a precursor generated by the carotenoid biosynthesis pathway [40]. Previous research has demonstrated that ABA could potentially regulate carotenoid composition during ripening by increased leaf chlorophylls, carotenoid, and fruit LYCO in tomato [41]. Furthermore, foliar-applied ABA on grapes (Vitis vinifera) results in stimulatory effects on fruit color [42]. Previous researches indicated that the cleavage reaction, catalyzed by NCED, is the first and regulatory step in ABA biosynthesis [43,44]. The expression levels of NCED1 and NCED3 are down-regulated from 40 DAF to 100 DAF at first and then up-regulated until 190 DAF in three pummelo cultivars, whereas that of NCED2 is up-regulated from 40 DAF to 100 DAF and then down-regulated until 190 DAF (Table 3 Figure 7F). The expression levels of NCED1 and NCED3 are higher at both of 40 DAF and 100 DAF in HRMY and HJMY than in GXMY (Table 3), suggesting that the higher expression levels of NCED1 and NCED3 may facilitate the biosynthesis of ABA in HRMY and HJMY. Applications of exogenous ABA increase concentrations of tomato leaf carotenoid and fruits, which might be due to the increased activity of key enzymes such as BC hydroxylase and PSY3 or both carotenoid, and chlorophylls are derived from the MEP pathway [45]. Sun et al. [46] suggested that carotenoid is the precursor of ABA synthesis, which could trigger ethylene production and then induce an increase in carotenoid biosynthesis in the pre-ripening stage of fruit. Conversely, ABA production can also be achieved by the change of carotenoid content and composition [47]. Evidence showed that the regulation of ABA contents in fruit is mediated by two main early-response catabolic enzyme AAO3 and CYP707A1 [48]. Here, we showed that AAO3 and CYP707A3 have a high expression level in the color break period in HJMY and HRMY (Table 3; Figure 7G,I). Expression of CYP707A1 shows an increasing trend in GXMY and HRMY and a decreasing trend in HJMY, which may involve degrading ABA at the later stages of fruit ripening and potentially preventing the over-accumulation of ABA at this stage. These results indicate the possibility to control carotenoid composition and content during fruit ripening by regulating ABA metabolism. Similarly, ABA is also involved in fruit ripening and affects carotenoid accumulation in tomato and strawberry [49,50]. However, whether phytohormones directly regulate the expression level of key genes in carotenoid biosynthesis pathway remains unknown.
Several transcription factors have been reported to be associated with carotenoid metabolism such as RAP2.2 and PIF1 in Arabidopsis [51,52], RIN, TAGL1, ERF6 and HB-1 in watermelon [53], NAC, MYB, AP2/ERF, WRKY, bZIP, and ARF in sweet potato [54]. Furthermore, many genes such as RIN, TAG1, ERF6, AP2a, SGR, NOR, SlNAC4, and PIF have been shown to encode the respective transcriptional factors that affect the expression levels of genes involved in carotenoid biosynthesis through the regulation of ethylene biosynthesis and signaling [55,56,57]. Lu et al. [4] indicated that CsMADS6 directly regulates LCYB1 and other carotenogenic genes to coordinately and positively modulate carotenoid metabolism in citrus. bHLH transcription factors have been reported to play an important role in carotenoid biosynthesis regulation under environmental signals [57,58]. In tomatoes, a PRE-like atypical bHLH gene SlPRE2 negatively regulates carotenoid accumulation during fruit ripening by repressing the expression of PSY1, PDS, and ZDS [59]. In citrus, CubHLH1 is considered to have a similar function like Arabidopsis AIFs, which is functionally characterized as a negative regulator of the brassinolide (BR) signaling pathway and might be directly involved in carotenoid metabolism [58]. Here, our results showed that the expression patterns of two bHLH gene homolog CsbHLH107 (CsbHLH107, c11225.graph_c0) and AtFBH4 (CsbHLH130, c27501.graph_c0) negatively correlate with the genes involved in carotenoid biosynthesis (Table S6 (Supplementary Materials); Figure 8A,B). WRKY genes are known to participate in growth and developmental processes, the biosynthesis of secondary metabolites as well as biotic and abiotic stress responses through modulation of ABA signaling pathways [60]. Recently, a sweet osmanthus OfWRKY3 was shown to positively regulate the OfCCD4 gene through binding to its promoter W-box [61]. Our results indicated that two WRKY gene homolog GhWRKY17 (c34311.graph_c0) and AtWRKY22 (c21980.graph_c0) positively correlate with the expression level of CCD4 and negatively correlate with NCED2 (Table S6 (Supplementary Materials); Figure 8C,D). This result implied that WRKY gene could be involved in regulating carotenoid metabolism by modulating ABA signaling pathways, but the specific regulatory mechanisms need to be further studied. NAC proteins have also been reported to be related to carotenoid biosynthesis in tomato and papaya [62,63]. In tomato, over-expression of SlNAC1 leads to a decline in the contents of carotenoid and ABA, decreasing the synthesis of ethylene [62]. In addition, the SNAC4 (the protein structure of SNAC4 is very similar to that of ANAC072) and SNAC9 positively and negatively regulate fruit ripening and carotenoid accumulation by affecting ethylene and ABA levels [64]. Our results indicated that pummelo homologs of NAC072 (c13516.graph_c0) and NAC74 (c9680.graph_c0) are negatively and positively correlated with genes involved in carotenoid biosynthesis, respectively (Table S6 (Supplementary Materials); Figure 8E,K).
Ethylene plays a critical role in the regulation of carotenoid accumulation and carotenogenic gene expression [65,66]. In some non-climacteric fruits like papayas, CpEIN3a acts as a key transcription factor in the ethylene signal transduction pathway and interacts with CpNAC2 to activate the transcription of a subset of carotenoid biosynthesis-related genes, such as CpPDS2/4, CpZDS, CpLCYe, and CpCHYb [67]. As a non-climacteric fruit, ethylene promotes the color of citrus fruits and it increases the transcription of carotenoid biosynthesis-related genes and chlorophyll degradation-related genes [68]. Here, we showed that the ERF2 (c31821.graph_c0) is increased with the development of these three pummelo fruits, whereas the ERF23 (c8622.graph_c0) decreases in GXMY and increases in HRMY and HJMY during the development of fruit (Table S5 (Supplementary Materials); Figure 8H,I). This result indicated that ERF23-mediated biological regulation of ethylene may have a functional role in carotenoid biosynthesis in HRMY and HJMY. Previous studies showed that MYB transcription factors play important regulatory roles in anthocyanins biosynthesis, which could be responsible for fruit color formation [69,70,71]. Recently, studies have shown that MYB transcription factors might also play a role in carotenoid biosynthesis [72]. In this study, six MYB transcription factors are positively correlated with carotenoid biosynthesis-related genes, while only two pummelo MYB transcription factors (c22501.graph_c0 and c23847.graph_c0) show a negatively correlation with the expression of structural genes involved in carotenoid biosynthesis and one MYB transcription factor (c21421.graph_c0) is positively correlated with six structural genes involved in carotenoid biosynthesis (Tables S5 and S6, Supplementary Materials).

4. Material and Methods

4.1. Plant Materials

Three pummelo cultivars, named GXMY (yellow pulp) and its bud mutants HRMY (red pulp) and HJMY (golden pulp), were manually harvested from plants grown in the pummelo-producing areas in Pinghe County, Zhangzhou City, Fujian Province, China. Fruit samples were picked at 40 DAF (young fruit; the pulp radius accounts for one-fourth of the fruit radius), 100 DAF (color break period; the pulp radius accounts for half of the fruit length), 190 DAF (ripe stage, harvested for commercial sale) from three trees on a clear day, respectively (Figure 1). All the trees were grown in a square of 4 m × 5 m and received a standard horticultural practice such as irrigation, sod cultivation, bagging, and pathogen prevention. Samples were collected at sunny noon, peeled, and sliced into appropriate pieces after their weight and diameter were measured. In all cases, there were three biological replicates per cultivar for each sampling date and each biological replicate was mixed from at least five fruits. The pulp pieced from each fruit or several fruits were wrapped into one aluminum foil bag and labeled. Then, the samples were frozen in liquid nitrogen and stored at −80 °C until analysis.

4.2. β-Carotene and Lycopene Extractions and Analysis by HPLC

β-carotene and lycopene were extracted from fruits and were analyzed by HPLC according to the method previously described by Xu et al. [73]. Extracts were analyzed on Agilent 1100 HPLC system equipped with a C18 Kromasil column (250 mm × 4.6 mm, 5 µm). The acetonitrile:methanol:dichloromethane ratio was 70:15:15 as the mobile phase, with a flow rate of 1.0 mL·min−1 isocratic elution and a 10 μL injection volume with a detection wavelength of 450 nm in 35 min. β-carotene and lycopene contents were calculated using a standard curve method and expressed as mg·kg−1 FW. There were three replicates for the measurement of β-carotene and lycopene.

4.3. RNA Extraction, cDNA Library Construction and Sequencing

Total RNA was extracted from each sample using the TRIzol reagent (Invitrogen, Carlsbad, USA) and DNase I was then used to remove DNA (Takara, Dalian, China). The quality, quantity, and integrity of the total RNA were evaluated using Nanodrop (IMPLEN, CA, USA), Qubit 2.0 (Life Technologies, CA, USA), and Aglient 2100 (Agilent Technologies, CA, USA). Briefly, the first cDNA chain was synthesized using random hexamers as a six-base random primer, and the second cDNA chain was synthesized by adding a buffer solution, dNTPs, RNase H, and DNA polymerase I. The purified double-stranded cDNAs were terminal repaired, tailed, and sequenced. Then, the fragment size was selected by AMPure XP (Beckman Coulter, Brea, CA, USA) beads and the cDNA library was obtained by PCR enrichment. Finally, libraries were loaded on Illumina/Solexa HiSeq2000 platform with a sequenced read length set to PE100.

4.4. Reads Assembly and Functional Annotation

The raw reads were filtered to remove adapter sequences and reads with low-quality sequences containing more than 10% of bases which had a Q-value (quality score) of ≤ 20. The clean reads were assembled de novo using the Trinity Method [74]. Using the longest transcription as a sample of unigene, nine samples of unigene were combined. To calculate unigene expression, we used the fragments per kilobase of transcript per million mapped reads (FPKM) method [75]. Unigenes differentially expressed between two samples were screened using a false discovery rate (FDR) of < 0.01 and |log2 (fold changes)| of ≥1 as the criteria. To get annotation unigene sequences, the BLASTX alignment (E-value < 10−5) was performed between unigenes and Nr, Swiss-Prot, KEGG, COG, GO databases. The Blast2GO program [76] was used to obtain GO annotations of the unigenes based on Nr annotations.
The RNA-Seq data were deposited in NCBI database (https://www.ncbi.nlm.nih.gov/sra/) [77] with an SRA accession number PRJNA517027.

4.5. Identification of Genes Related to the Carotenoid Biosynthesis Pathways and Transcription Factor

According to the selected candidate gene information, genes responsible for the biosynthetic pathway of carotenoid were searched using the annotation database. Furthermore, the Arabidopsis Information Resource (TAIR) and sweet orange (C. sinensis L. Osbeck) annotation databases were used to analyze the homologous sequence similarly of the structural genes. The identified and searched sequences were confirmed by the NCBI BLAST program [14]. In order to understand the underlying transcriptional mechanisms of three pummelo fruits, we used the transcription factor prediction module PlantTFDB 3.0 [78] with a default E-value of 10−6.

4.6. Correlation Analysis of Structural Genes and Transcription Factors

A correlation analysis of structural genes and transcription factors was carried out to identify transcription factors that were co-expressed with the functional genes involved in carotenoid biosynthesis. In order to exclude false positives, functional genes and transcription factors with an FPKM value of ≥10.0 in at least one of the three stages during fruit development were selected, and transcription factors with correlation coefficient values of ≥0.666 by the t-test (the formula to calculate t value was t = (r√(n–2))/√(1–r2), at p < 0.05 and n = 9. |t| > t0.05,7 = 2.365 means significant correlation, so r > 0.666 means significant correlation. The co-expression analysis was performed by “CORREL” function in “EXCEL 2003” and confirmed by an in-house Perl scripts and IBM SPSS Statistics software.

4.7. RT-qPCR analysis

To validate the reliability of our transcriptome dataset, 10 putative candidate functional genes and 12 transcription factors selected from the DEGs were further analyzed using RT-qPCR. Total RNA of nine fruit samples was extracted using the TRIzol® Reagent (Invitrogen, USA). RT-qPCR was performed on a CFX96 real-time system (BIO-RAD, Hercules, CA, USA) with a TaKaRa 2×SYBR Premix ExTaqTM kit (Takara, Dalian, China). RT-qPCR was performed under the following conditions: at 95 °C for 30 s, followed by 40 cycles at 95 °C for 5 s and at 72 °C for 30 s. Relative expression analysis was performed using the 2−ΔΔCT method. The transcript levels were normalized relative to the actin gene. The RT-qPCR primer sequences are listed in Table S7. Three biological and three technical replicates were used for the RT-qPCR assays.

4.8. Experimental Design and Statistical Analysis

There were three biological replicates per cultivar for each sampling date and each biological replicate was mixed from at least five fruits. Considering that HRMY and HJMY were the bud mutation of GXMY and fruit samples at different development stages and used, one biological replicate of mixed sample was used for RNA-Seq in each time point for each cultivar. There were three replicates for the measurement of β-carotene, lycopene, and RT-qPCR assays. Differences among the nine samples were analyzed by three (cultivars) × three (sampling dates) ANOVA. Means were separated by the Duncan’s new multiple range test at p < 0.05.

5. Conclusions

Our results clearly showed that HJMY has the highest content of β-carotene, followed by HRMY and GXMY. Lycopene is dominantly accumulated in HRMY. By using RNA-Seq technique, we identified 357 genes involved in biosynthesis of secondary metabolites, of which 12 DEGs are involved in carotenoid biosynthesis. Among these 12 DEGs, PSY2, LYCB2, LYCE, CCD4, NCED2, AAO3, and CYP707A1 are the most distinct DEGs in three pummelo cultivars. The co-expression analysis revealed that the expression patterns of several transcription factors such as bHLH, MYB, ERF, NAC, and WRKY are highly correlated with DEGs, which are involved in carotenoid biosynthesis. In addition, the expression patterns of 22 DEGs are validated by real-time RT-qPCR and the results are highly concordant with the RNA-Seq results. In conclusion, our results provide a global vision of transcriptomic profile among three pummelo cultivars with different pulp colors and would be beneficial to further study the molecular mechanism of carotenoid accumulation in pummelo flesh and help the breeding of citrus with high carotenoid content.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/20/9/2246/s1. Table S1. List of KEGG pathway annotations in transcriptome of three pummelo cultivars. Table S2. The pathway enrichment analysis of DEGs during fruit development in three pummelo cultivars. Table S3. List of 435 transcription factor genes in three pummelo cultivars. Table S4. The number and percentage of transcription factors belonging to different families. Table S5. Differentially expressed transcription factor genes. Table S6. Correlation analysis of transcription factors and structural genes. Table S7. Primers sequences used for RT-qPCR analysis.

Author Contributions

C.C. conducted the data analysis and drafted the manuscript. Y.-F.Z. and Y.-J.L. participated in the RNA-Seq and the measurement of biochemical compounds. Y.C. conducted the RT-qPCR analysis. X.-K.L. designed the study and revised the manuscript. All the authors have read and approved the final manuscript.

Funding

This study was funded by the National Facilities and Information Infrastructure for Science and Technology (NICGR2018-069), Basic Scientific Research Funds of Public Welfare Scientific Research Institutes of Fujian Province (2018R1013-12 and 2016R1013-4), and the PI Project of the Innovation Team at Fujian Academy of Agricultural Sciences (STIT2017-1-4).

Acknowledgments

We sincerely thank Xin-Fu Ye (Fujian Academy of Agricultural Sciences) for the assistance in the experimental design.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. 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] [PubMed]
  2. Kato, M.; Ikoma, Y.; Matsumoto, H.; Sugiura, M.; Hyodo, H.; Yano, M. Accumulation of carotenoids and expression of carotenoid biosynthetic genes during maturation in citrus fruit. Plant Physiol. 2004, 134, 824–837. [Google Scholar] [CrossRef] [PubMed]
  3. Zhu, F.; Luo, T.; Liu, C.Y.; Wang, Y.; Yang, H.B.; Yang, W.; Zheng, L.; Xiao, X.; Zhang, M.F.; Xu, R.W.; et al. An R2R3-MYB transcription factor represses the transformation of alpha- and beta-branch carotenoids by negatively regulating expression of CrBCH2 and CrNCED5 in flavedo of Citrus reticulate. New Phytol. 2017, 216, 178–192. [Google Scholar] [CrossRef] [PubMed]
  4. Lu, S.W.; Zhang, Y.; Zhu, K.J.; Yang, W.; Ye, J.L.; Chai, L.J.; Xu, Q.; Deng, X.X. The citrus transcription factor CsMADS6 modulates carotenoid metabolism by directly regulating carotenogenic genes. Plant Physiol. 2018. [Google Scholar] [CrossRef]
  5. Ikoma, Y.; Matsumoto, H.; Kato, M. Diversity in the carotenoid profiles and the expression of genes related to carotenoid accumulation among citrus genotypes. Breed. Sci. 2016, 66, 139–147. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Rodrigo, M.J.; Alquézar, B.; Alós, E.; Lado, J.; Zacarías, L. Biochemical bases and molecular regulation of pigmentation in the peel of Citrus fruit. Sci. Hortic. 2013, 163, 46–62. [Google Scholar] [CrossRef]
  7. Yu, K.Q.; Xu, Q.; Da, X.L.; Guo, F.; Ding, Y.D.; Deng, X.X. Transcriptome changes during fruit development and ripening of sweet orange (Citrus sinensis). Bmc Genom. 2012, 13, 1–13. [Google Scholar] [CrossRef] [PubMed]
  8. Wu, J.; Xu, Z.; Zhang, Y.; Chai, L.; Yi, H.; Deng, X.X. An integrative analysis of the transcriptome and proteome of the pulp of a spontaneous late-ripening sweet orange mutant and its wild type improves our understanding of fruit ripening in citrus. J. Exp. Bot. 2014, 65, 1651–1671. [Google Scholar] [CrossRef] [Green Version]
  9. Zhang, Y.J.; Wang, X.J.; Wu, J.X.; Chen, S.Y.; Chen, H.; Chai, L.J.; Yi, H.L. Comparative transcriptome analyses between a spontaneous late ripening sweet orange mutant and its wild type suggest the functions of ABA, sucrose and JA during citrus fruit ripening. PLoS ONE 2014, 9, 116056. [Google Scholar] [CrossRef]
  10. Guo, F.; Yu, H.; Xu, Q.; Deng, X. Transcriptomic analysis of differentially expressed genes in an orange-pericarp mutant and wild type in pummelo (Citrus grandis). Bmc Plant Biol. 2015, 15, 44. [Google Scholar] [CrossRef]
  11. Liu, C.H.; Yan, F.H.; Gao, H.J.; He, M.; Wang, Z.; Cheng, Y.J.; Deng, X.X.; Xu, J. Features of citrus terpenoid production as revealed by carotenoid, limonoid and aroma profiles of two pummelos (Citrus maxima) with different flesh color. J. Sci. Food Agric. 2015, 95, 111–119. [Google Scholar] [CrossRef]
  12. Lu, X.K.; Lin, Q.H.; Lu, X.M.; Zhang, S.M.; Li, C.S.; Ye, X.F. Comparison of carotenoid compositions and contents n different sweet pomelos. Fujianj Agric. Sci. 2012, 27, 723–727. [Google Scholar]
  13. Lu, X.K.; Lin, Q.H.; Lin, Y.J.; Zhang, J.T.; Zhang, S.M.; Li, C.S. ‘Huangjinmiyou’, a new orange-yellow fleshed pomelo cultivar. J. Fruit Sci. 2013, 30, 900–902. [Google Scholar]
  14. BLASTX and BLASTN searching engines. Available online: https://blast.ncbi.nlm.nih.gov/Blast.cgi.
  15. Ampomah-Dwamena, C.; McGhie, T.; Wibisono, R.; Montefiori, M.; Hellens, R.P.; Allan, A.C. The kiwifruit lycopene beta-cyclase plays a significant role in carotenoid accumulation in fruit. J. Exp. Bot. 2009, 60, 3765–3779. [Google Scholar] [CrossRef] [Green Version]
  16. Lv, P.; Li, N.; Gu, H.H.; Zhao, W.E. Changes in carotenoid profiles and in the expression pattern of the genes in carotenoid metabolisms during fruit development and ripening in four watermelon cultivar. Food Chem. 2015, 174, 52–59. [Google Scholar] [CrossRef]
  17. Mizuno, K.; Tokiwano, T.; Yoshizawa, Y. Gene expression analysis of enzymes of the carotenoid biosynthesis pathway involved in β-cryptoxanthin accumulation in wild raspberry, Rubus palmatus. Biochem. Biophys. Commun. 2017, 484, 845–849. [Google Scholar] [CrossRef] [PubMed]
  18. Leng, X.P.; Wang, P.P.; Wang, C.; Zhu, X.D.; Li, X.P.; Li, H.Y.; Mu, Q.M.; Liu, Z.J.; Fang, J.G. Genome-wide identification and characterization of genes involved in carotenoid metabolic in three stages of grapevine fruit development. Sci. Rep. 2017, 7, 4216. [Google Scholar] [CrossRef] [Green Version]
  19. Richaud, D.; Stange, C.; Gadaleta, A.; Colasuonno, P.; Parada, R.; Schwember, A.R. Identification of Lycopene epsilon cyclase (LCYE) gene mutants to potentially increase β-carotene content in durum wheat (Triticum turgidum L.ssp. Durum) through TILLING. PLoS ONE 2018, 13, e0208948. [Google Scholar]
  20. Colasuonno, P.; Lozito, M.L.; Marcotuli, I.; Nigro, D.; Giancaspro, A.; Mangini, G.; De Vita, P.; Mastrangelo, A.M.; Pecchioni, N.; Houston, K.; et al. The carotenoid biosynthetic and catabolic genes in wheat and their association with yellow pigments. Bmc Genom. 2017, 18, 122. [Google Scholar] [CrossRef]
  21. Colasuonno, P.; Incerti, O.; Lozito, M.L.; Simeone, R.; Gadaleta, A.; Blanco, A. DHPLC technology for high-throughput detection of mutations in a durum wheat TILLING population. Bmc Genet. 2016, 17, 43. [Google Scholar] [CrossRef] [PubMed]
  22. Charles, A.D.; Driedonks, N.; Lewis, D.; Shumskaya, M.; Chen, X.Y.; Wurtzel, E.T.; Espley, R.V.; Allan, A.C. The Phytoene synthase gene family of apple (Malus x domestica) and its role in controlling fruit carotenoid content. Bmc Plant Biol. 2015, 15, 185. [Google Scholar]
  23. Clotault, J.; Peltier, D.; Berruyer, R.; Thomas, M.; Briard, M.; Geoffriau, E. Expression of carotenoid biosynthesis genes during carrot root development. J. Exp. Bot. 2008, 59, 3563–3573. [Google Scholar] [CrossRef] [Green Version]
  24. Fraser, P.D.; Eugenia, M.A.A.; Halket, J.M.; Truesdale, M.R.; Yu, D.M.; Gerrish, C.; Bramley, P.M. Manipulation of phytoene levels in tomato fruit: Effects on isoprenoids; plastids; and intermediary metabolism. Plant Cell 2007, 19, 3194–3211. [Google Scholar] [CrossRef]
  25. Rodrigo, M.J.; Marcos, J.F.; Zacarias, L. Biochemical and molecular analysis of carotenoid biosynthesis in flavedo of orange (Citrus sinensis L.) during fruit development and maturation. J. Agric. Food Chem. 2004, 52, 6724–6731. [Google Scholar] [CrossRef]
  26. Wei, X.; Chen, C.X.; Yu, Q.B.; Gady, A.; Yu, Y.; Liang, G.L.; Gmitter, F.G. Comparison of carotenoid accumulation and biosynthesis gene expression between Valencia and Rohde Red Valencia sweet oranges. Plant Sci. 2014, 227, 28–36. [Google Scholar] [CrossRef]
  27. Guo, F.; Zhou, W.; Zhang, J.; Xu, Q.; Deng, X.X. Effect of the citrus lycopene β-cyclase transgene on carotenoid metabolism in transgenic tomato fruits. PLoS ONE 2012, 7, e32221. [Google Scholar] [CrossRef] [PubMed]
  28. Mendes, A.F.S.; Chen, C.X.; Gmitter, F.G.; Moore, G.A.; Costa, M.G.C. Expression and phylogenetic analysis of two new lycopene β-cyclases from Citrus paradisi. Physiol. Plant 2011, 141, 1–10. [Google Scholar] [CrossRef]
  29. Chen, C.; Costa, M.G.C.; Yu, Q.; Moore, G.C.; Gmitter, F.G. Identification of novel members in sweet orange carotenoid biosynthesis gene families. Tree Genet. Genome 2010, 6, 905–914. [Google Scholar] [CrossRef]
  30. Kang, C.; Zhai, H.; Xue, L.Y.; Zhao, N.; He, S.Z.; Liu, Q.C. A lycopene β-cyclase gene; IbLCYB2; enhances carotenoid contents and abiotic stress tolerance in transgenic sweetpotato. Plant Sci. 2018, 272, 243–254. [Google Scholar] [CrossRef]
  31. Devitt, L.C.; Fanning, K.; Dietzgen, R.G.; Holton, T.A. Isolation and functional characterization of a lycopene β-cyclase gene that controls fruit colour of papaya (Carica papaya L.). J. Exp. Bot. 2010, 61, 33. [Google Scholar] [CrossRef] [PubMed]
  32. Auldridge, M.E.; McCarty, D.R.; Klee, H.J. Plant carotenoid cleavage oxygenases and their apocarotenoid products. Curr. Opin. Plant Biol. 2006, 9, 315–321. [Google Scholar] [CrossRef] [PubMed]
  33. Bouvier, F.; Suire, C.; Mutterer, J.; Camara, B. Oxidative remodeling of chromoplast carotenoids: Identification of the carotenoid dioxygenase CsCCD and CsZCD genes involved in Crocus secondary metabolite biogenesis. Plant Cell 2003, 15, 47–62. [Google Scholar] [CrossRef]
  34. Yoshioka, S.; Aida, R.; Yamamizo, C.; Shibata, M.; Ohmiya, A. The carotenoid cleavage dioxygenase 4 (CmCCD4a) gene family encodes a key regulator of petal color mutation in chrysanthemum. Euphytica 2012, 184, 377–387. [Google Scholar] [CrossRef]
  35. Ureshino, K.; Nakayaa, M.; Miyajima, I. Contribution made by the carotenoid cleavage dioxygenase 4 gene to yellow colour fade in azalea petals. Euphytica 2016, 207, 401–417. [Google Scholar] [CrossRef]
  36. Song, C.; Zhang, L.; Zhang, C.; Tian, Y.; Cong, P. Differential expression analysis of the fruits of yellow-fleshed apple and exploits genes involved in carotenoid pigmentation. Acta Hortic. Sin. 2017. [Google Scholar] [CrossRef]
  37. Bai, S.; Tuan, P.A.; Tatsuki, M.; Yaegaki, H.; Ohmiya, A.; Yamamizo, C.; Moriguchi, T. Knockdown of carotenoid cleavage dioxygenase 4 (CCD4) via virus-induced gene silencing confers yellow coloration in peach fruit: Evaluation of gene function related to fruit traits. Plant Mol. Biol. Rep. 2016, 34, 257–264. [Google Scholar] [CrossRef]
  38. Bruno, M.; Beyer, P.; Al-Babili, S. The potato carotenoid cleavage dioxygenase 4 catalyzes a single cleavage of biononering-containing carotenes and non-epoxidated xanthophylls. Arch. Biochem. Biophys. 2015, 572, 126–133. [Google Scholar] [CrossRef]
  39. González-Verdejo, C.; Obrero, A.; Román, B.; Gómez, P. Expression profile of carotenoid cleavage dioxygenase genes in summer squash (Cucurbita pepo L.). Plant Foods Hum. Nutr. 2015, 70, 200–206. [Google Scholar]
  40. Schwartz, S.H.; Qin, X.; Zeevaart, J.A. Elucidation of the indirect pathway of abscisic acid biosynthesis by mutants, genes, and enzymes. Plant Physiol. 2003, 131, 1591–1601. [Google Scholar] [CrossRef]
  41. Barickman, T.C.; Kopsell, D.A.; Sams, C.E. Abscisic acid increases carotenoid and chlorophyll concentrations in leaves and fruit of two tomato genotypes. J. Am. Soc. Hortic. Sci. 2014, 139, 261–266. [Google Scholar] [CrossRef]
  42. Cantin, C.M.; Fidelibus, M.W.; Crisostoc, C.H. Application of abscisic acid (ABA) at veraison advanced red color development and maintained postharvest quality of ‘Crimson Seedless’ grapes. Postharvest Biol. Technol. 2007, 46, 237–241. [Google Scholar] [CrossRef] [Green Version]
  43. Colasuonno, P.; Marcotuli, I.; Lozito, M.L.; Simeone, R.; Blanco, A.; Gadaleta, A. Characterization of aldehyde oxidase (AO) genes involved in the accumulation of carotenoid pigments in wheat grain. Bmc Genet. 2017, 8, 863. [Google Scholar] [CrossRef]
  44. Chernys, J.T.; Zeevaart, J.A. Characterization of the 9-cis-epoxycarotenoid dioxygenase gene family and the regulation of abscisic acid biosynthesis in avocado. Plant Physiol. 2000, 124, 343–353. [Google Scholar] [CrossRef]
  45. Meier, S.; Tzfadia, O.; Vallabhaneni, R.; Gehring, C.; Wurtzel, E.T. A transcriptional analysis of carotenoid; chlorophyll and plastidial isoprenoid biosynthesis genes during development and osmotic stress responses in Arabidopsis thaliana. Bmc Syst. Biol. 2011, 5, 77. [Google Scholar] [CrossRef]
  46. 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 beta-carotene contents in tomato fruit. J. Exp. Bot. 2012, 63, 3097–3108. [Google Scholar] [CrossRef]
  47. Fang, J.; Chai, C.; Qian, Q.; Li, C.; Tang, J.; Sun, L.; Huang, Z.; Guo, X.; Sun, C.; Liu, M.; et al. Mutations of genes in synthesis of the carotenoid precursors of ABA lead to preharvest sprouting and photooxidation in rice. Plant J. 2008, 54, 177–189. [Google Scholar] [CrossRef]
  48. Nambara, E.; Marion-Poll, A. Abscisic acid biosynthesis and catabolism. Annu. Rev. Plant Biol. 2005, 56, 165–185. [Google Scholar] [CrossRef]
  49. Ji, K.; Kai, W.; Zhao, B.; Sun, Y.; Yuan, B.; Dai, S.; Li, Q.; Chen, P.; Wang, Y.; Pei, Y.; et al. SlNCED1 and SlCYP707A2: Key genes involved in ABA metabolism during tomato fruit ripening. J. Exp. Bot. 2014, 65, 5243–5255. [Google Scholar] [CrossRef]
  50. Li, D.D.; Li, L.; Luo, Z.S.; Mou, W.S.; Mao, L.C.; Ying, T.J. Comparative transcriptome analysis reveals the influence of abscisic acid on the metabolism of pigments; ascorbic acid and folic acid during strawberry fruit ripening. PLoS ONE 2015, 10, e0130037. [Google Scholar] [CrossRef]
  51. Welsch, R.; Beyer, P.; Hugueney, P.; Kleinig, H.; von Lintig, J. Regulation and activation of phytoene synthase, a key enzyme in carotenoid biosynthesis, during photomorphogenesis. Planta 2000, 211, 846–854. [Google Scholar] [CrossRef]
  52. Toledo-Ortiz, G.; Huq, E.; Rodriguez-Concepcion, M. Direct regulation of phytoene synthase gene expression and carotenoid biosynthesis by phytochrome-interacting factors. Proc. Natl. Acad. Sci. USA 2010, 107, 11626–11631. [Google Scholar] [CrossRef] [Green Version]
  53. Grassi, S.; Piro, G.; Lee, J.M.; Zheng, Y.; Fei, Z.J.; Dalessandro, G.; Giovannoni, J.J.; Lenucci, M.S. Comparative genomics reveals candidate carotenoid pathway regulators of ripening watermelon fruit. Bmc Genom. 2013, 14, 781. [Google Scholar] [CrossRef]
  54. Li, R.; Zhai, H.; Kang, C.; Liu, D.; He, S.Z.; Liu, Q.C. De Novo transcriptome sequencing of the orange-fleshed sweet potato and analysis of differentially expressed genes related to carotenoid biosynthesis. Int. J. Genom. 2015, 2015, 1–10. [Google Scholar] [CrossRef] [PubMed]
  55. Lee, J.M.; Joung, J.G.; McQuinn, R.; Chung, M.Y.; Fei, Z.; Tieman, D.; Klee, H.; Giovannoni, J.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]
  56. Zhu, M.K.; Chen, G.P.; Zhou, S.; Tu, Y.; Wang, Y.; Dong, T.T.; Hu, Z.L. 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]
  57. Llorente, B.D.; Andrea, L.; Ruizsola, M.A.; Botterweg, E.; Pulido, P. Tomato fruit carotenoid biosynthesis is adjusted to actual ripening progression by a light-dependent mechanism. Plant J. 2016, 85, 107–119. [Google Scholar] [CrossRef] [PubMed]
  58. Endo, T.; Fujii, H.; Sugiyama, A.; Nakano, M.; Nakajima, N.; Ikoma, Y.; Omura, M.; Shimada, T. Overexpression of a citrus basic helix-loop-helix transcription factor (CubHLH1), which is homologous to Arabidopsis activation-tagged bri1 suppressor 1 interacting factor genes; modulates carotenoid metabolism in transgenic tomato. Plant Sci. 2016, 243, 35–48. [Google Scholar] [CrossRef] [PubMed]
  59. Zhu, Z.G.; Chen, G.P.; Guo, X.H.; Yin, W.C.; Yu, X.H.; Hu, J.T.; Hu, Z.L. Overexpression of SlPRE2, an atypical bHLH transcription factor; affects plant morphology and fruit pigment accumulation in tomato. Sci. Rep. 2017, 7, 5786. [Google Scholar] [CrossRef]
  60. Banerjee, A.; Roychoudhury, A. WRKY proteins: Signaling and regulation of expression during abiotic stress responses. Sci. World J. 2015, 807560. [Google Scholar] [CrossRef]
  61. Han, Y.J.; Wu, M.; Cao, L.Y.; Yuan, W.J.; Dong, M.; Wang, X.H.; Chen, W.C.; Shang, F.D. Characterization of Of WRKY3; a transcription factor that positively regulates the carotenoid cleavage dioxygenase gene Of CCD4 in Osmanthus fragrans. Plant Mol. Biol. 2016, 91, 485–496. [Google Scholar] [CrossRef] [PubMed]
  62. Ma, N.N.; Feng, H.L.; Meng, X.; Li, D.; Yang, D.Y.; Wu, C.G.; Meng, Q.W. Overexpression of tomato SlNAC1 transcription factor alters fruit pigmentation and softening. Bmc Plant Biol. 2014, 14, 351. [Google Scholar] [CrossRef]
  63. Fu, C.C.; Han, Y.C.; Fan, Z.Q.; Chen, J.Y.; Chen, W.X.; Lu, W.J.; Kuang, J.F. The papaya transcription factor cpnac1 modulates carotenoid biosynthesis through activating phytoene desaturase genes CpPDS2/4 during fruit ripening. J. Agric. Food Chem. 2016, 64, 5454. [Google Scholar] [CrossRef] [PubMed]
  64. Kou, X.H.; Zhao, Y.N.; Wu, C.; Jiang, B.L.; Zhang, Z.; Rathbun, J.R.; He, Y.L.; Xue, Z.H. SNAC4 and SNAC9 transcription factors show contrasting effects on tomato carotenoids biosynthesis and softening. Postharvest Biol. Tech. 2018, 144, 9–19. [Google Scholar] [CrossRef]
  65. Giovannoni, J.J. Genetic regulation of fruit development and ripening. Plant Cell 2004, 16, S170–S180. [Google Scholar] [CrossRef] [PubMed]
  66. Welsch, R.; Maass, D.; Voegel, T.; Dellapenna, D.; Beyer, P. Transcription factor RAP2.2 and its interacting partner SINAT2: Stable elements in the carotenogenesis of Arabidopsis leaves. Plant Physiol. 2007, 45, 1073–1085. [Google Scholar] [CrossRef] [PubMed]
  67. Fu, C.C.; Han, Y.C.; Kuang, J.F.; Chen, J.Y.; Lu, W.J. Papaya CpEIN3a and CpNAC2 cooperatively regulate carotenoid biosynthesis related genes CpPDS2/4, CpLCY-e, and CpCHY-b during fruit ripening. Plant Cell Physiol. 2017, 58, 2155–2165. [Google Scholar] [CrossRef] [PubMed]
  68. Fujii, H.; Shimada, T.; Sugiyama, A.; Nishikawa, F.; Endo, T.; Nakano, M.; Ikoma, Y.; Shimizu, T.; Omura, M. Profiling ethylene-responsive genes in mature mandarin fruit using a citrus 22K oligoarray. Sci. Hortic. 2007, 173, 340–348. [Google Scholar] [CrossRef]
  69. Allan, A.C.; Hellens, R.P.; Laing, W.A. MYB transcription factors that colour our fruit. Trends Plant Sci. 2008, 13, 99–102. [Google Scholar] [CrossRef] [PubMed]
  70. Deluc, L.; Bogs, J.; Walker, A.R.; Ferrier, T.; Decendit, A.; Merillon, J.M.; Robinson, S.P.; Barrieu, F. The transcription factor VvMYB5b contributes to the regulation of anthocyanin and proanthocyanidin biosynthesis in developing grape berries. Plant Physiol. 2008, 147, 2041–2053. [Google Scholar] [CrossRef] [PubMed]
  71. Martins, T.R.; Jiang, P.; Rausher, M.D. How petals change their spots: Cis-regulatory re-wiring in Clarkia (Onagraceae). New Phytol. 2016, 216, 510. [Google Scholar] [CrossRef]
  72. Sagawa, J.M.; Stanley, L.E.; LaFountain, A.M.; Frank, H.A.; Liu, C.; Yuan, Y.W. An R2R3-MYB transcription factor regulates carotenoid pigmentation in Mimulus lewisii flowers. New Phytol. 2016, 209, 1049–1057. [Google Scholar] [CrossRef] [PubMed]
  73. Xu, J.; Tao, N.G.; Liu, Q.; Deng, X.X. Presence of diverse rations of lycopene β–carotene in five pink or red-fleshed citrus cultivars. Sci. Hortic. 2006, 108, 181–184. [Google Scholar] [CrossRef]
  74. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644. [Google Scholar] [CrossRef] [PubMed]
  75. Li, B.; Dewey, C. RSEM: Accurate transcript quantification from RNASeq data with or without a reference genome. Bmc Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef]
  76. Conesa, A.; Götz, S.; García-Gómez, J.M.; Terol, J.; Talón, M.; Robles, M. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21, 3674–3676. [Google Scholar] [CrossRef] [PubMed]
  77. NCBI database. Available online: https://www.ncbi.nlm.nih.gov/sra.
  78. PlantTFDB 3.0 software. Available online: http://planttfdb.cbi.pku.edu.cn/prediction.php.
Ijms 20 02246 g001 550
Figure 1. Three pummelo cultivars at different periods of fruit growth and development.
Figure 1. Three pummelo cultivars at different periods of fruit growth and development.
Ijms 20 02246 g001
Ijms 20 02246 g002 550
Figure 2. β-carotene and lycopene contents in three pummelo cultivars of different fruit growth and development periods. Bars represent means ± SE (n = 3). Differences among the samples were analyzed by three (cultivars) × three (sampling dates) ANOVA. Different letters indicate a significant difference at p < 0.05.
Figure 2. β-carotene and lycopene contents in three pummelo cultivars of different fruit growth and development periods. Bars represent means ± SE (n = 3). Differences among the samples were analyzed by three (cultivars) × three (sampling dates) ANOVA. Different letters indicate a significant difference at p < 0.05.
Ijms 20 02246 g002
Ijms 20 02246 g003 550
Figure 3. Sequence length distribution of the unigenes in fruit transcriptomes of three pummelo cultivars. The x-axis indicates the unigene length interval from 200 bp to >2000 bp. The y-axis indicates the number of unigenes of each given sequence length.
Figure 3. Sequence length distribution of the unigenes in fruit transcriptomes of three pummelo cultivars. The x-axis indicates the unigene length interval from 200 bp to >2000 bp. The y-axis indicates the number of unigenes of each given sequence length.
Ijms 20 02246 g003
Ijms 20 02246 g004 550
Figure 4. Classification based on “Biosynthesis of secondary metabolites” categories.
Figure 4. Classification based on “Biosynthesis of secondary metabolites” categories.
Ijms 20 02246 g004
Ijms 20 02246 g005 550
Figure 5. Venn diagram illustrating the number of unigenes expressed of three pummelo cultivars at different developmental stages.
Figure 5. Venn diagram illustrating the number of unigenes expressed of three pummelo cultivars at different developmental stages.
Ijms 20 02246 g005
Ijms 20 02246 g006 550
Figure 6. Diagram showing a proposed model for carotenoid biosynthesis of three pummelo cultivars. Biosynthesis pathways are shown with solid arrows and regulatory interactions are shown with broken arrows. Up arrow (↑), down arrow (↓), and vertical bar (|) indicate significantly up-regulated, down-regulated and no obvious change of gene expression of these enzymes, respectively. The colored arrows or bars represent the change of gene expression in GXMY (gray), HRMY (red), and HJMY (golden).
Figure 6. Diagram showing a proposed model for carotenoid biosynthesis of three pummelo cultivars. Biosynthesis pathways are shown with solid arrows and regulatory interactions are shown with broken arrows. Up arrow (↑), down arrow (↓), and vertical bar (|) indicate significantly up-regulated, down-regulated and no obvious change of gene expression of these enzymes, respectively. The colored arrows or bars represent the change of gene expression in GXMY (gray), HRMY (red), and HJMY (golden).
Ijms 20 02246 g006
Ijms 20 02246 g007 550
Figure 7. Gene expression level analysis of 10 differentially expressed structural genes related to carotenoid biosynthesis of three pummelo cultivars at different developmental stages. Bars represent means ± SE (n = 3). Differences among the samples were analyzed by three (cultivars) × three (sampling dates) ANOVA. Different letters indicate a significant difference at p < 0.05.
Figure 7. Gene expression level analysis of 10 differentially expressed structural genes related to carotenoid biosynthesis of three pummelo cultivars at different developmental stages. Bars represent means ± SE (n = 3). Differences among the samples were analyzed by three (cultivars) × three (sampling dates) ANOVA. Different letters indicate a significant difference at p < 0.05.
Ijms 20 02246 g007
Ijms 20 02246 g008 550
Figure 8. Expression analysis of 12 differentially expressed transcription factors related to carotenoid biosynthesis of three pummelo cultivars at different developmental stages. Bars represent means ± SE (n = 3). Differences among the samples were analyzed by three (cultivars) × three (sampling dates) ANOVA. Different letters indicate a significant difference at p < 0.05.
Figure 8. Expression analysis of 12 differentially expressed transcription factors related to carotenoid biosynthesis of three pummelo cultivars at different developmental stages. Bars represent means ± SE (n = 3). Differences among the samples were analyzed by three (cultivars) × three (sampling dates) ANOVA. Different letters indicate a significant difference at p < 0.05.
Ijms 20 02246 g008
Table
Table 1. Summary of transcriptomic sequencing and de novo assembly.
Table 1. Summary of transcriptomic sequencing and de novo assembly.
SamplesRead numberBase numberGC content% (≥Q30)
GXMY1 (40 DAF)22,744,1514,581,699,08446.49%85.92%
GXMY2 (100 DAF)22,575,1784,548,053,12146.10%85.79%
GXMY3 (190 DAF)22,780,0994,589,605,41145.77%85.82%
HRMY1 (40 DAF)22,437,9534,515,689,21647.34%85.73%
HRMY2 (100 DAF)22,964,1754,624,970,83645.06%86.05%
HRMY3 (190 DAF)22,607,9524,557,243,04644.76%85.88%
HJMY1 (40 DAF)22,064,2214,450,273,68146.50%85.68%
HJMY2 (100 DAF)19,694,9503,977,639,96645.25%85.02%
HJMY3 (190 DAF)23,914,6704,816,251,31045.48%86.02%
Total201,783,34940,661,425,67145.86%85.76%
After assembly
Total clean reads201,783,349
Total clean nucleotides (nt)40,661,425,671
Total number of contigs17,658,287
Total length of contigs (nt)656,824,421
Mean length of contigs (nt)37
N50 of contigs37
Total number of unigenes54,051
Total length of unigenes (nt)46,385,997
Mean length of unigenes (nt)858
N50 of unigenes1,585
Table
Table 2. Differentially expressed genes at fruit developmental stages in three pummelo cultivars.
Table 2. Differentially expressed genes at fruit developmental stages in three pummelo cultivars.
Fruit development at different stages (days)GXMYHRMYHJMY
ALL DEGsUp-regulatedDown-regulatedALL DEGsUp-regulatedDown-regulatedALL DEGsUp-regulatedDown-regulated
40 DAF/100 DAF3240106821721887371151620804211,659
40 DAF/190 DAF3882140824742569753181623166391,677
100 DAF/190 DAF256213001262498818753113337117961,575
40 DAF/100 DAF/190 DAF616 258 256
Table
Table 3. Expression profiles of carotenoid biosynthesis-related genes at fruit developmental stages in three pummelo cultivars.
Table 3. Expression profiles of carotenoid biosynthesis-related genes at fruit developmental stages in three pummelo cultivars.
Gene NameUnigene IDGene LengthGXMYHRMYHJMY
40 DAF100 DAF190 DAF40 DAF100 DAF190 DAF40 DAF100 DAF190 DAF
PSY2c27977.graph_c013733.450.410.040.040.32000.280.04
LYCB2c35136.graph_c020375.5715.1519.492.315.0725.570.5815.419.32
LYCEc29745.graph_c020186.224.011.640.716.954.680.637.783.09
CCD1c29184.graph_c09036.043.810.651.45.870.471.265.310
CCD4c20564.graph_c0305760.5520.5916.3696.2818.7516.19118.9917.2412.32
NCED1c17328.graph_c1 15368.230.455.2012.651.0910.872.230.638.46
NCED2c31914.graph_c0242517.56107.1299.7110.13105.1184.134.5797.9959.51
NCED3c11942.graph_c0110111.740.467.2814.570.5110.623.261.2411.98
AAO3c33740.graph_c0455922.2718.467.992.519.436.511.2222.025.97
CYP707A1c20275.graph_c020412.112.977.525.594.1427.6485.485.5813.47
CYP707A3c27956.graph_c018113.925.032.611.364.262.131.113.720.46
CYP707A4c35953.graph_c0173417.328.771.161.229.110.890.8111.710.59
Table
Table 4. Correlation analyses of structural genes involved in carotenoid metabolism and transcription factors.
Table 4. Correlation analyses of structural genes involved in carotenoid metabolism and transcription factors.
Gene IDFPKM (max)FPKM (min)Description for the best hit in C. sinensisNumber of correlations
c11225.graph_c017.950.69bHLH1074
12645.graph_c042.681.02BIM1 isoform X13
c21405.graph_c03.450bHLH1622
c22086.graph_c09.370bHLH572
c2272.graph_c02.600bHLH352
c22910.graph_c014.150.46bHLH614
c27501.graph_c021.445.78bHLH1304
c589.graph_c13.970HEC22
c6149.graph_c02.740ILR34
c34091.graph_c0256.506.32Ethylene-responsive transcription factor ERF123
c24274.graph_c012.571.31Ethylene-responsive transcription factor-like protein3
c34297.graph_c0138.233.30Ethylene-responsive transcription factor ERF1073
c407.graph_c13.130Ethylene-responsive transcription factor 132
c8122.graph_c0111.065.19Ethylene-responsive transcription factor RAP2-34
c8622.graph_c057.791.26Ethylene-responsive transcription factor ERF232
c30160.graph_c022.310.04AP2-like ethylene-responsive transcription factor ANT1
c26695.graph_c02.830AP2-like ethylene-responsive transcription factor AIL14
c30127.graph_c0118.3913.99Auxin response factor 43
c27816.graph_c05.910Growth-regulating factor 13
c9781.graph_c05.380Growth-regulating factor 42
c22501.graph_c08.350.91Transcriptional adapter ADA2b isoform X16
c23847.graph_c01296.81106.33MYB1R1-like4
c25031.graph_c02.340MYB114
c9240.graph_c06.370.20MYB522
c1320.graph_c04.490MYB132
c13169.graph_c0104.341.68Transcription factor AS1-like6
c13009.graph_c0157.182.18REVEILLE 12
c21421.graph_c023.560.45Transcription factor AS16
c13063.graph_c03047.44242.76Mini zinc finger protein 24
c21980.graph_c06.220.19WRKY224
c34311.graph_c0114.9321.37WRKY175
c13516.graph_c026.604.39NAC723
c31393.graph_c020.633.05Transcription factor RF2a-like3
c28763.graph_c023.770.83Transcription factor TGA73
c34869.graph_c058.155.27BES1/BZR1 homolog protein 46
c35277.graph_c022.807.03Transcription factor TCP203
c35957.graph_c017.191.35Trihelix transcription factor GTL12
c31746.graph_c0330.3315.64DELLA protein GAI5
c31940.graph_c04.210FAR1-related sequence 5-like5
c28566.graph_c013.071.05Transcription factor HHO2-like6
c29740.graph_c026.661.14B3 domain-containing transcription factor VRN1-like4
c27534.graph_c0303.9913.66Effector of transcription 24
c27216.graph_c0110.233.11Transcription factor TCP46
c22856.graph_c07.470Trihelix transcription factor GT-24
c20987.graph_c030.860.03GATA transcription factor 114
c14734.graph_c063.086.70Transcription factor HHO35
c12838.graph_c049.599.91YABBY 2 isoform X15
c10930.graph_c022.791.61Trihelix transcription factor PTL6

Share and Cite

MDPI and ACS Style

Jiang, C.-C.; Zhang, Y.-F.; Lin, Y.-J.; Chen, Y.; Lu, X.-K. Illumina® Sequencing Reveals Candidate Genes of Carotenoid Metabolism in Three Pummelo Cultivars (Citrus Maxima) with Different Pulp Color. Int. J. Mol. Sci. 2019, 20, 2246. https://doi.org/10.3390/ijms20092246

AMA Style

Jiang C-C, Zhang Y-F, Lin Y-J, Chen Y, Lu X-K. Illumina® Sequencing Reveals Candidate Genes of Carotenoid Metabolism in Three Pummelo Cultivars (Citrus Maxima) with Different Pulp Color. International Journal of Molecular Sciences. 2019; 20(9):2246. https://doi.org/10.3390/ijms20092246

Chicago/Turabian Style

Jiang, Cui-Cui, Yan-Fang Zhang, Yan-Jin Lin, Yuan Chen, and Xin-Kun Lu. 2019. "Illumina® Sequencing Reveals Candidate Genes of Carotenoid Metabolism in Three Pummelo Cultivars (Citrus Maxima) with Different Pulp Color" International Journal of Molecular Sciences 20, no. 9: 2246. https://doi.org/10.3390/ijms20092246

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