Skip to main content

Weighted gene coexpression correlation network analysis reveals the potential molecular regulatory mechanism of citrate and anthocyanin accumulation between postharvest ‘Bingtangcheng’ and ‘Tarocco’ blood orange fruit

Abstract

Background

Organic acids and anthocyanins are the most important compounds for the flavor and nutritional quality of citrus fruit. However, there are few reports on the involvement of co-regulation of citrate and anthocyanin metabolism. Here, we performed a comparative transcriptome analysis to elucidate the genes and pathways involved in both citrate and anthocyanin accumulation in postharvest citrus fruit with ‘Tarocco’ blood orange (TBO; high accumulation) and ‘Bingtangcheng’ sweet orange (BTSO; low accumulation).

Results

A robust core set of 825 DEGs were found to be temporally associated with citrate and anthocyanin accumulation throughout the storage period through transcriptome analysis. Further according to the results of weighted gene coexpression correlation network analysis (WGCNA), the turquoise and brown module was highly positively correlated with both of the content of citrate and anthocyanin, and p-type ATPase (PH8), phosphoenolpyruvate carboxylase kinase (PEPCK), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonoid 3’-hydroxylase (F3’H) and glutathione S transferase (GST) were considered key structural genes. Moreover, MYB family transcription factor (PH4), Zinc finger PHD-type transcription factor (CHR4, HAC12), Zinc finger SWIM-type transcription factor (FAR1) and Zinc finger C3H1-type transcription factor (ATC3H64) were considered hub genes related to these structural genes. Further qRT-PCR analysis verified that these transcription factors were highly expressed in TBO fruit and their expression profiles were significantly positively correlated with the structural genes of citrate and anthocyanin metabolism as well as the content of citrate and anthocyanin content.

Conclusions

The findings suggest that the CHR4, FAR1, ATC3H64 and HAC12 may be the new transcription regulators participate in controlling the level of citrate and anthocyanin in postharvest TBO fruit in addition to PH4. These results may providing new insight into the regulation mechanism of citrate and anthocyanin accumulation in citrus fruit.

Peer Review reports

Background

Citrus is one of the most important fruit crops worldwide. China is the largest producer of citrus in the world with a current annual yield of 45 million tons, among which 90% are used for fresh consumption. With the changes in people’s consumption concept, ‘high quality’ has become the core competitiveness of fresh citrus fruit, and postharvest fresh-keeping technologies play an increasingly prominent role in citrus industry. Organic acids are important components for the flavor and quality of citrus fruit, and are closely related to the fruit storage performance [1, 2]. In addition, anthocyanin is an important component for citrus fruit color, as well as plays a critical role in citrus fruit quality formation, stress resistance and human health benefits [3,4,5]. Blood orange is the only commercial citrus variety that can accumulate anthocyanins in its fruit. The content of organic acid and anthocyanin jointly determine the quality of postharvest blood orange fruit.

The intracellular organic acid level is regulated by the synthesis, degradation, transport and vacuole storage of organic acids. A variety of enzymes involved in citrate metabolism have been reported, such as phosphoenolpyruvate carboxylase (PEPC) [6], phosphoenolpyruvate carboxylase kinase (PEPCK) [7], citrate synthase (CS) [8], cis-aconitase (Aco) [9], vacuolar proton pump V-ATPase/V-PPase and citrate/H+ cotransporter (Cit) [10, 12], GABA shunt gene and citrate lyase (ACL) [13, 14]. Moreover, transcription factors (TFs) including ERF, bHLH, NAC, WRKY, bZIP and MYB have also been reported to participate in the regulation of citrate level through the VHA [15], Aco [16, 17], ACL [18] or PH8 [12, 19, 20] pathway.

The biosynthesis of anthocyanins is a complex multi-step process involving many structural and regulatory genes. These structural genes code enzymes like phenylalanine ammonia-lyase (PAL), chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonoid 3’-hydroxylase (F3’H), dihydroflavonol 4-reductase (DFR), anthocyanidin synthase (ANS) and glutathione S transferase (GST) [21, 22]. The most extensively studied TFs for the transcriptional regulation of anthocyanin biosynthesis are MYB, bHLH and WD40 [23]. It has been demonstrated that MYB can regulate anthocyanin accumulation by directly binding to the promoters of key genes in the anthocyanin biosynthesis pathway [24, 25]. In general, bHLH needs to form a complex with the MYB protein to be involved in anthocyanin biosynthesis [4, 26,27,28]. Similarly, WD40 cannot independently regulate anthocyanin synthesis and mainly functions by forming MYB-bHLH-WD40 (MBW) complex [29]. It has been reported that anthocyanin accumulation in almost all citrus cultivars is closely associated with the activity of the Ruby1 gene (a MYB TF) [30]. The Ruby1 gene can respectively interact with bHLH and WD40 to regulate anthocyanin metabolism in citrus [4, 31].

Although numerous genes related to citrate and anthocyanin metabolism have been reported, little is known about the common regulatory pathways involved in the accumulation of both citrate and anthocyanins. In petunia, it has been reported that MYB complexes with bHLH and WD40 protein to regulate anthocyanin biosynthesis, and bHLH is actually the AN1 protein involved in regulating vacuolar acidification [32]. A recent study demonstrated that the traits of exceptionally low fruit acidity and absence of anthocyanins in leaves and flowers and proanthocyanidins in citrus seeds are due to the mutation of the Noemi gene encoding a bHLH TF (a homolog of AN1 in petunia and TT8 in Arabidopsis) [33]. Because both acidity and coloration are important fruit quality traits, increasing research has been focused on the discovery of new Noemi genes [34].

Transcriptome profiling has greatly contributed to the identification of new genes or pathways for many important traits [21, 35, 36]. Blood orange is the only commercial citrus variety that can accumulate anthocyanins in its fruit, and the acid content in its fruit is higher among sweet orange varieties. ‘Tarocco’ blood orange (TBO) is the main variety of blood orange grown in Hunan. ‘Bingtangcheng’ sweet orange (BTSO) is an excellent variety selected in Hunan province, and is characterized by a low acid content. It is also the main sweet orange variety planted in Hunan province. In this study, we comparatively analyzed the transcriptome of these two contrasting sweet orange germplasm, TBO and BTSO, throughout the postharvest storage process, aiming to identify the genes playing important roles in the postharvest process of citrus fruit, and thereby provide insights into the key pathways and regulators simultaneously associated with the accumulation of citrate and anthocyanins in citrus fruit.

Results

Variations in the physical and chemical characteristics of TBO and BTSO during postharvest storage

To compare the fruit of TBO and BTSO, we examined the physiochemical characteristics of fruit at different stages of postharvest storage (Fig. 1). The two varieties exhibited similar fruit pulp colors at maturity (0 DAS, days after storage). Thereafter, the pulp of TBO fruit started to accumulate anthocyanins gradually from 15 to 90 DAS and was differentiated from BTSO fruit in color (Fig. 1A). Further determination of anthocyanin content revealed that it was gradually increased in TBO during storage and was significantly higher than that in BTSO (Fig. 1C).

Fig. 1
figure 1

BTSO and TBO fruit under postharvest storage (A), and their corresponding citrate and anthocyanin contents (B and C). Error bars represent the standard deviations of the mean in three replicates, and for each storage day, * stand for significant differences between two materials at p < 0.05

Citrate and malate are major organic acids in citrus fruit. Hence, we also analyzed their changes in TBO and BTSO fruit during postharvest storage. The results showed that mature TBO (10.49 ± 0.72 g kg−1) had a markedly higher citrate content than BTSO (1.49 ± 0.05 g kg−1), which was consistently observed throughout the whole storage period (Fig. 1B). The malate content showed a decreasing tendency in TBO fruit, while a slightly increasing tendency in BTSO fruit along with storage (Fig. S1). In general, the malate content had little influence on total organic acids.

Transcriptome analysis

To explore the regulatory mechanism of citrate and anthocyanin accumulation in the flesh of citrus fruit, the samples at five detection time points (0, 15, 30, 60 90 DAS) from the two different varieties were used for deep RNA-seq analysis. After filtering out the rRNAs and low-quality reads, a total of 133 million reads were mapped to the Citrus sinensis reference genome (Table S2). For these clean reads, the average mapped reads per sample was greater than 90% (ranging from 87.90% to 91.87%). A total of 29,655 annotated genes were obtained. Differentially expressed genes (DEGs) were identified based on their expression levels in different samples, and functional annotation and enrichment analysis were performed. A total of 14,232 genes were differentially expressed between postharvest stored TBO and BTSO fruit (Fig. 2A). To screen the candidate genes related to citrate and anthocyanin biosynthesis, our study focused on the DEGs at TBO vs. BTSO. There were 2397, 2677, 3067, 3131 and 2960 genes that were differentially expressed on 0, 15, 30, 60 and 90 DAS respectively between TBO and BTSO (Fig. S2). Among them, 839 genes were consistently differentially expressed (Fig. 2B), and 381 genes were consistently upregulated, while 444 were consistently downregulated in TBO relative to BTSO throughout the storage period. The remaining 14 genes were upregulated in some storage periods and downregulated in other storage periods (Fig. 2C).

Fig. 2
figure 2

Summary of differentially expressed genes between TBO and BTSO fruit during postharvest storage. A Number of DEGs in different DEG sets. B Venn diagram shows DEGs in TBO vs BTSO. C Number statistics of DEGs intersected of TBO vs BTSO in five storage period groups

Construction of a WGCNA and coexpression network

To obtain hub genes related to citrate and anthocyanin accumulation, the relationships of DEGs, citrate and anthocyanin content for each sample were analyzed by constructing a WGCNA (Fig. 3). Samples clustering showed that the three biological replicates of each sample were very good (Fig. 3A). Fourteen coexpression modules were identified by WGCNA (Fig. 3B), among which the turquoise module was highly positively correlated with the contents of citrate (r = 0.57, p value = 9e-04) and anthocyanin (r = 0.83, p value = 1e-08). In addition, the brown module was also highly positively correlated with the citrate content (r = 0.94, p value = 2e-14) (Fig. 3C). A total of 5154 and 1193 genes were obtained from turquoise and brown module respectively. Thinking about the content of citrate and anthocyanin were all markedly higher in TBO than that in BTSO fruit during the whole storage period (Fig. 1B and C), we focused on intersection genes between the 825 continually DEGs (Fig. 2C) and the turquoise and brown modules. Finally, 310 valuable genes was filtered out. According to GO and KEGG enrichment of the 310 candidate genes, 11 genes were mapped to the flavonoid metabolism pathway, 5 genes were mapped to the citrate cycle pathway, and 7 genes were mapped to the starch and sucrose metabolism pathway (Fig. S3 and S4).

Fig. 3
figure 3

Weighted gene coexpression network analysis of TBO vs BTSO fruit during postharvest storage period. A Sample clustering. B Hierarchical clustering showing modules of coexpression genes. C Module/trait correlations and corresponding p values. The right panel shows a colour scale for module/trait correlations from -1 to 1

Identification of candidate genes involved in citrate and anthocyanin accumulation

As a results, 3 and 5 structural genes involved in citrate metabolism and anthocyanin biosynthesis were obtained in 310 candidate genes. A heatmap of their expression profiles in the flesh of TBO and BTSO fruit was drawn based on their FPKM value (log2(FPKM)) (Figs. 4 and 5). The 3 structural genes from all major steps of the citrate metabolism pathway were distributed as follows: p-type ATPase gene (PH8), two phosphoenolpyruvate carboxylase kinase (PEPCK). Further the expression patterns of these 3 structural genes involved in citrate metabolism, PH8 (Cs1g16150), PEPCK1 (Cs3g16700) and PEPCK2 (Cs1g20920), were studied via qRT-PCR (quantitative real-time PCR), and the transcripts of PH8 and PEPCK2 genes were significantly higher in TBO fruit than in BTSO fruit during the whole storage period (Fig. 4C), which was consistent with the results of transcriptome analysis based on the correlation analysis (Fig. 4D).

Fig. 4
figure 4

Analysis of genes related to citrate metabolism in the intersection genes between the 825 continually DEGs and the turquoise and brown modules. A Citrate metabolism pathway. The red font indicates the genes obtained in the intersection genes. B Heatmap of the expression levels of differentially expressed genes (DEGs) involved in citrate metabolism. C qRT-PCR detection of citrate metabolism related structural genes, PH8 (Cs1g16150), PEPCK1 (Cs3g16700) and PEPCK2 (Cs1g20920). D Correlation analysis of the expression profiles in qRT-PCR (qPCR) and transcriptome data (FPKM) (*, ** and *** represent significant differences at p < 0.05, p < 0.01 and p < 0.001, respectively between the two sets of data)

Fig. 5
figure 5

Analysis of genes related to anthocyanin biosynthesis in the intersection genes between the 825 continually DEGs and the turquoise and brown modules. A Anthocyanin biosynthesis pathway. The red font indicates the genes obtained in the intersection genes. B Heatmap of the expression levels of differentially expressed genes (DEGs) involved in anthocyanin biosynthesis. C qRT-PCR detection of anthocyanin biosynthesis related structural genes, CHI (Cs7g29780), F3H (Cs1g25280 and Cs2g04110), F3’H (Cs5g11730) and GST (Cs6g15900). D Correlation analysis of the expression profiles in qRT-PCR (qPCR) and transcriptome data (FPKM) (*, ** and *** represent significant differences at p < 0.05, p < 0.01 and p < 0.001, respectively between the two sets of data)

The 5 structural genes from all major steps of the anthocyanin biosynthesis pathway were distributed as follows: chalcone isomerase (CHI), two flavanone 3-hydroxylase (F3H), flavonoid 3’-hydroxylase (F3’H) and glutathione S transferase (GST). Further the expression patterns of these 5 structural genes involved in anthocyanin biosynthesis, CHI (Cs7g29780), F3H (Cs1g25280 and Cs2g04110), F3’H (Cs5g11730) and GST (Cs6g15900), were studied via qRT-PCR, and the transcripts of these genes were significantly higher in TBO fruit than in BTSO fruit during the whole storage period (Fig. 5C), which was consistent with the results of transcriptome analysis based on the correlation analysis (Fig. 5D).

Identification of genes involved in carbohydrate metabolism

Carbohydrates are considered substrates for anthocyanin synthesis, and was also closely connect with organic acid metabolism. In our results, the contents of glucose, fructose and sucrose showed downward trends on the whole in TBO and BTSO fruit during postharvest storage. The glucose and fructose content in TBO fruit was significantly lower than that in BTSO throughout the storage period; while the sucrose content was markedly higher in TBO fruit than that in BTSO fruit at 15 and 60 DAS (Fig. 6A). In 310 candidate genes, two trehalose-phosphate phosphatase (TPP) genes were involved in starch and sugar metabolism which was related to citrate metabolism and anthocyanin synthesis (Fig. 6B). Correlation analysis showed that the expression patterns of the TPP genes were positively correlated with citrate and anthocyanin content, and was negatively correlated with glucose and fructose content (Fig. 6C), suggesting that higher expression levels of these genes were beneficial to carbohydrate metabolism, which contributed to citrate and anthocyanin accumulation in TBO fruit.

Fig. 6
figure 6

Soluble sugar content and expression pattern of carbohydrate metabolism related genes in TBO and BTSO fruit during postharvest storage period. A Glucose, fructose and sucrose contents in TBO fruit compared with BTSO fruit during storage. B Heatmap of the expression levels of carbohydrate metabolism related genes in the intersection genes between the 825 continually DEGs and the turquoise and brown modules. C Correlation analysis of the citrate content, anthocyanin content, soluble sugar content and expression profiles of carbohydrate metabolism related genes (*, ** and *** represent significant differences at p < 0.05, p < 0.01 and p < 0.001, respectively between the two sets of data). For each storage day, different lowercase letters stand for significant differences between two materials at p < 0.05

Screening potential transcription factors that regulate citrate and anthocyanin accumulation in postharvest TBO fruit

To explore the molecular regulatory mechanism of citrate metabolism and anthocyanin biosynthesis in TBO fruit, we constructed a coexpression network based on the genes present in the intersection of 825 continually differentially expressed genes (DEGs) (Fig. 2C) and turquoise + brown modules. In the network, we identified a total of 40 transcription factor genes, and 5 of them including PH4 (Cs9g03070), CHR4 (Cs2g11400), FAR1 (Cs2g11450), ATC3H64 (Cs4g14220), and HAC12 (Cs6g19010) were identified as hub genes (Fig. 7A and B). These TFs were found to be closely correlated with citrate and anthocyanin metabolism-related genes. Among them, PH4 was closely related to PH8, CHI, F3’H, and GST; CHR4 was closely related to F3H (Cs1g25280); FAR1 was closely related to PH8, PEPCK1, and PEPCK2; ATC3H64 was closely related to PEPCK2 and F3H (Cs1g25280 and Cs2g04110); HAC12 was closely related to PEPCK1 (Fig. 7B).

Fig. 7
figure 7

Coexpression network analysis of the potential key genes in the intersection between the 825 continually DEGs and the turquoise and brown modules. A coexpression network of the genes (weight > 0.15). B coexpression network of transcription factors and citrate/anthocyanin metabolism related structural genes. C Heatmap of the expression levels of the candidate transcription factors. D qRT-PCR detection of candidate transcription factor genes CHR4 (Cs2g11400), FAR1 (Cs2g11450), ATC3H64 (Cs4g14220), HAC12 (Cs6g19010) and PH4 (Cs9g03070). E Correlation analysis of the expression profiles in qRT-PCR (qPCR) and transcriptome data (FPKM) (*, ** and *** represent significant differences at p < 0.05, p < 0.01 and p < 0.001, respectively between the two sets of data)

The heatmap based on the transcriptome data showed that the expression levels of PH4, CHR4, HAC12, FAR1 and ATC3H64 were obviously upregulated in the TBO fruit compared with BTSO fruit throughout the storage period (Fig. 7C). Meanwhile, the qRT-PCR results indicated that the transcripts of these 5 detected genes in the TBO fruit were obviously higher than those in BTSO fruit during storage, especially the FAR1, ATC3H64 and PH4 (Fig. 7D), which was consistent with the transcriptome analysis results based on the correlation analysis (Fig. 7E). Furthermore, the expression patterns of the 5 transcription factor genes were positively correlated with the changes of structural genes expression levels and citrate/anthocyanin content (Fig. S5). It was proposed that PH4, CHR4, HAC12, FAR1 and ATC3H64 might be involved in citrate and anthocyanin accumulation by regulating the expression of structural genes associated with citrate metabolism and anthocyanin biosynthesis, especially the FAR1, ATC3H64 and PH4 genes.

Discussion

The mechanism of higher content of citrate and anthocyanin in postharvest TBO fruit

Accumulations of citrate and anthocyanins are important indicators to evaluate citrus fruit quality, as well as important traits for consumers’ concern. Blood orange is the only commercial citrus variety that can accumulate anthocyanins in its fruit, and with higher organic acid content in its fruit compared with other sweet orange, which determined its unique flavor. Some studies have suggested that PH8 might be responsible for the difference in acidity among mature or postharvest stored citrus varieties [1, 2, 7]. Here, the differentially expressed structural gene PH8 was found based on the transcriptome data, and was mapped to the citrate metabolism pathway (Fig. 4A and B). Moreover, its expression pattern was significantly positively correlated with citrate content (Fig. 4 C and D), which was similar to previous study. These results indicate that in TBO fruit, more citrate accumulated in vacuoles due to the strikingly higher expression of PH8, which increases the ability of the proton pump to uptake citrate into vacuoles.

Additionally, a series of studies showed that the accumulation of anthocyanin in blood orange fruit was attributed to the upregulation of genes across the biosynthetic pathway [37, 38]. In this study, the differentially expressed structural gene CHI, F3H, F3’H and GST were found based on the transcriptome data, and was mapped to the anthocyanin biosynthesis pathway (Fig. 5A and B). The expression pattern of these genes was significantly positively correlated with anthocyanin content, especially two F3H genes (Fig. 5C and D). These results indicate that the high expression of F3H might be the main reason for anthocyanin accumulation during postharvest storage of TBO fruit, and the significant high expression of F3’H and GST also jointly promoted the anthocyanin accumulation.

Transcription factors involved in citrate and anthocyanin accumulation in postharvest TBO fruit

The regulation of citrate metabolism has always been the concern of researchers. And so far, some transcription factors have been reported to be involved in the regulation of citrate metabolism through triggering structural genes [12, 33, 39,40,41]. MYB5a, MYB5b and MYBA1 have been proven to be involved in vacuolar acidification by inducing the expression of vacuolar acidification-related genes PH5 and CAC16.5 in grapevine [40]. The co-exprssion of PH4 and AN1 strongly induced PH1 and PH5 expression to lead to citrate accumulation in citrus fruit [12, 41]. VvMYB5 has similar functions to PhPH4 and AtMYB5, and is also involved in controlling vaculor hyper-acidification and trafficking in grapvine [39]. In the present study, both transcriptome analysis and qPCR quantitative detection that the expression of PH4 was significantly higher in TBO fruit than that in BTSO, and was coordinated with the expression patterns of PH8 (Fig. 7). Thus, the transcription factor PH4 is likely involved in regulating citrate accumulation in postharvest TBO fruit. Except for PH4, the expression patterns of FAR1 was also found to be closely related to citrate accumulation and changes in the PH8, and moreover, the expression patterns of FAR1 was similar to PH4, indicating that the transcription factors FAR1 was also involved in regulating citrate accumulation in postharvest TBO fruit. The proposed function of FAR is that it play crucial roles in controlling the growth and development of plants, and it also proved to involved in regulating light-induced myo-inositol biosynthesis and an alternative ascorbate biosynthetic pathway [42, 43]. However, the roles in regulating citrate accumulation in postharvest fruit are unclear and need to be studied further.

Previous studies have well documented that the transcription factors families such as the R2R3-MYB, bHLH and WD40-domain proteins primarily regulate anthocyanin biosynthesis via regulating the structural genes in the anthocyanin biosynthesis pathway [24, 28]. In grapevine, the MYB5b can involve in anthocyanin biosynthesis through upregulating a subset of anthocyanin structural genes [40]. Additionally, overexpression of LcMYB5 increased anthocyanin biosynthesis in tobacco and petunia either by directly activating the expression of key structural genes such as DFR or by indirectly up regulating the expression of endogenous bHLH1 [44]. Recently, research showed that PH4 could directly upregulating the expression of Noemi gene and involved in proanthocyanin biosynthesis through activating the expression of PA biosynthetic genes like DFR, ANS, ANR, LAR and UFGT2 in citrus [45]. In this study, the expression of PH4 was also coordinated with the expression patterns of CHI, F3’H and GST (Fig. 7). Thus, the transcription factor PH4 was also likely to be involved in regulating anthocyanin biosynthesis in postharvest TBO fruit. Moreover, the expression patterns of CHR4, FAR1, ATC3H64 and HAC12 were similar to PH4, indicating that these TFs were have similar functions to PH4. A FAR1 gene, 16 and 15 members from C3H and ZF subfamilies was identified and differentially expressed in purple-fleshed compared with white-fleshed sweetpotato through transcriptome analysis [46], suggesting that these transcription factors involved in anthocyanin accumulation in sweetpotato. However, the regulation mechanism of these transcription factors in anthocyanin biosynthesis in postharvest fruit are unclear and need to be studied further.

In recent years, many reports have proposed that the regulation of citrate metabolism is closely associated with that of anthocyanin biosynthesis [32, 33, 47]. The bHLH protein AN1 (also called Noemi) participates in the regulation of vacuolar acidification in petunia, and is also involved in the regulation of anthocyanin biosynthesis by complexing with MYB and WD40 protein [32]. The exceptionally low fruit acidity and the absence of anthocyanins in citrus leaves and flowers are due to mutation in the Noemi gene [33]. PH4 could directly upregulating the expression of Noemi gene and involved in both citrate accumulation and proanthocyanidin biosynthesis [45]. In this study, we found that the expression of AN1 was also upregulated in TBO fruit and was significantly positive correlated with PH4 (Fig. 7B). And so far, the transcription regulation involved in both citrate and anthocyanin metabolism are all around MYB and bHLH. Hence, with the results of this study, we hope to explore new regulatory factors from CHR4, FAR1, ATC3H64 and HAC12 that can simultaneously regulate the accumulation of citrate and anthocyanin. According to the results of previous studies [1, 4] and the present study, we established a model to explain the high citrate and anthocyanin accumulation in postharvest TBO fruits compared with BTSO fruits (Fig. 8).

Fig. 8
figure 8

Sketch of the potential regulatory mechanisms involved in high citrate and anthocyanin trait in TBO fruits compared with that in BTSO fruits. The red typeface indicates that the metabolite content or gene expression level is significantly higher in TBO fruits than BTSO fruits

Conclusion

In this study, TBO and BTSO were chosen for a transcriptome analysis to excavate the key genes and regulators involved in citrate and anthocyanin metabolism in citrus fruit during postharvest storage. Importantly, we identified the key structural genes, transporters and TFs for the genetic control of citrate and anthocyanin accumulation in citrus fruit. Since in almost all transcriptome studies, it remains a great challenge to demonstrate which candidate genes play actual regulatory roles in fruit acidity and coloration, our comprehensive analysis provides an important basis for future study to determine whether the key genes identified here play important and specific roles in citrate and anthocyanin accumulation in citrus fruit.

Materials and methods

Plant materials

Physiologically mature fruit of ‘Tarocco’ blood orange (TBO) and ‘Bingtangcheng’ sweet orange (BTSO) (Citrus sinensis) were purchased from an experimental station in Huaihua, Hunan Province, China, in December 2019, and then immediately transported to the laboratory. Fruit with uniform size and color and free of any visible damage or defects were selected as samples for further experiments. The fruit were individually packed in polyethylene bags and stored in a temperature controlled chamber at 8 ± 1 °C with a relative humidity 85%–90% for about three months.

Juice sacs were separated from five fruit in each group at 0, 15, 30, 60 and 90 days after storage (DAS) with three replicates to measure the levels of organic acids, sugars, and for transcriptome analysis and expression analysis of related genes. The samples were immediately frozen in liquid nitrogen and stored at –80 °C for further analysis.

Determination of organic acids and sugars

Organic acids (citrate and malate) and sugars (glucose, fructose and sucrose) were measured and analyzed by high-performance liquid chromatography (HPLC) with a previously described method [1]. The results were expressed on a fresh weight (FW) basis.

Determination of anthocyanins

The anthocyanin contents were determined based on a previously reported protocol [48] with minor modifications. The mixed sample was ground in liquid nitrogen, and exactly 0.1 g of the sample was homogenized with 2 mL of A buffer (50 mM KCl and 150 mM HCl, pH = 1.0). Another 0.1 g of the sample was homogenized with 2 mL of B buffer (400 mM NaAC and 240 mM HCl, pH = 4.5). The homogenate was centrifuged at 12,000 rpm for 15 min at 4 °C. The supernatant was then collected and diluted to measure the A510 values, and the anthocyanin content (g kg−1 FW) = (A510 at pH1.0 – A510 at pH4.5) × 484.8 (molecular weight of cyanidin) / 24,825 (molar absorption coefficient of cyanidin at A510) × dilution ratio.

Transcriptome sequencing and WGCNA

Total RNA was extracted from TBO and BTSO fruit samples for 0, 15, 30, 60 and 90 DAS (three biological repetitions) using the Trizol reagent (Invitrogen, CA, USA) following the manufacturer's procedure. The mRNA was purified and fragmented into small pieces. Then, the cleaved RNA fragments were reverse transcribed to create the final cDNA library in accordance with the protocol for the mRNA Seq sample preparation kit (Illumina, San Diego, USA). RNA sequencing was performed on an Illumina Hiseq 4000 platform (LC-Bio, Hangzhou, China) to produce raw reads according to the protocol.

High-quality clean reads were used to de novo assemble the transcriptomes of TBO and BTSO using the StringTie program [49]. Then, all transcriptomes from samples were merged to reconstruct a comprehensive transcriptome using the Perl scripts. After the final transcriptome was generated, StringTie and Ballgown were used to estimate the expression levels of all genes. StringTie was used to determine the expression level for mRNAs by calculating the FPKM. Differentially expressed mRNAs and genes were selected with log2 (fold change) > 1 or log2 (fold change) < -1 and with statistical significance (p value < 0.05) by R package -Ballgown [50]. Gene Ontology (GO) annotation was performed using the eggNOG-mapper (http://eggnog-mapper.embl.de/). DEGs were enriched in GO and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases so as to identify the changes in biological functions and metabolism pathways [51, 53]. Finally, the transcriptome data and phenotype data were analyzed by WGCNA using R language. The WGCNA package was employed to construct a gene coexpression network using a variant set of genes (15,318 genes, the maximum FPKM in all samples ≥ 2). The analysis was performed based on the WGCNA package in R studio software [54]. The main parameters of WGCNA program were as follows: variance data expression > 0,no mission data expression < 0.1; soft threshold = 18 (scale-free R2 = 0.9); deep split = 2; min module size = 60; merge cut height = 0.25.

Quantitative real-time PCR validation

The candidate DEGs were verified by qRT-PCR. Total RNA from samples was extracted by the CTAB method and used for quantitative real-time PCR.

The RNA samples were used for first-strand cDNA synthesis using a FastKing gDNA Dispelling RT SuperMix reverse transcriptase Kit (Tiangen, China) following the manufacturer’s instructions. qRT-PCR assays were performed according to previous reports [1], and the relative expression levels were calculated with the 2−ΔΔCt method in three biological replications. Specific primers were designed from the selected gene sequences using Primer Express 3.0 (Applied Biosystems, Foster City, CA, USA) and the primer sequences are given in Table S1.

Statistical analysis

Each experiment was performed in three replicates. Experimental results were analyzed using IBM SPSS Statistics, RStudio and Gephi software. Error bars denote standard deviations. Different lowercase letters above the bars indicate significant differences at p < 0.05, which were obtained based on one-way ANOVA using IBM SPSS Statistics software.

Availability of data and materials

The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE212015 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE212015.

Abbreviations

TBO:

‘Tarocco’ blood orange

BTSO:

‘Bingtangcheng’ sweet orange

WGCNA:

Weighted gene coexpression correlation network analysis

PEPC:

Phosphoenolpyruvate carboxylase

PEPCK:

Phosphoenolpyruvate carboxylase kinase

CS:

Citrate synthase

Aco:

Cis-aconitase

Cit:

Citrate/H+ cotransporter

ACL:

Citrate lyase

TF:

Transcription factors

PAL:

Phenylalanine ammonia-lyase

CHS:

Chalcone synthase

CHI:

Chalcone isomerase

F3H:

Flavanone 3-hydroxylase

F3’H:

Flavonoid 3’-hydroxylase

DFR:

Dihydroflavonol 4-reductase

ANS:

Anthocyanidin synthase

GST:

Glutathione S transferase

MBW:

MYB-bHLH-WD40

PH8:

P-type ATPase gene

TPP:

Trehalose-phosphate phosphatase

References

  1. Ma XQ, Li N, Guo J, Yang LQ, Hao CX, Li Y, Gentile A, Lu XP, Ma XF, Deng ZN, et al. Involvement of CsPH8 in citrate accumulation directly related to fruit storage performance of ‘Bingtang’ sweet orange mutants. Postharvest Biol Technol. 2020;170:111316.

    Article  CAS  Google Scholar 

  2. Sheng L, Shen DD, Yang W, Zhang MF, Zeng YL, Xu J, Deng XX, Cheng YJ. GABA pathway rate-limit citrate degradation in postharvest citrus fruit evidence from HB Pumelo (Citrus grandis) × Fairchild (Citrus reticulata) hybrid population. J Agric Food Chem. 2017;65:1669–76.

    Article  CAS  PubMed  Google Scholar 

  3. Alvarez-Suarez JM, Guadrado G, Redondo IB, Giampieri F, González-Paramás AM, Santos-Buelga G. Novel approaches in anthocyanin research-Plant fortification and bioavailability issues. Trends Food Sci Technol. 2021;117:92–105.

    Article  CAS  Google Scholar 

  4. Huang D, Wang X, Tang ZZ, Yuan Y, Xu YT, He JX, Jiang XL, Peng SA, Li L, Butelli E, et al. Subfunctionalization of the Ruby2-Ruby1 gene cluster during the demestication of citrus. Nat Plants. 2018;4:930–41.

    Article  CAS  PubMed  Google Scholar 

  5. Ma HY, Yang T, Li Y, Zhang J, Wu T, Song TT, Yao YC, Tian J. The long noncoding RNA MdLNC499 bridges MdWRKY1 and MdERF109 function to regulate early-stage light-induced anthocyanin accumulation in apple fruit. Plant Cell. 2021;33:3309–30.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Guillet C, Just D, Bénard N, Destrac-Irvine A, Baldet P, Hernould M, Causse M, Raymond P, Rothan C. A fruit-specific phospho enolpyruvate carboxylase is related to rapid growth of tomato fruit. Planta. 2002;214:717–26.

    Article  CAS  PubMed  Google Scholar 

  7. Guo LX, Shi CY, Liu X, Ning DY, Jing LF, Yang H. Liu YZ Citrate accumulation-related gene expression and/or enzyme activity analysis combined with metabolomics provide a novel insight for an orange mutant. Sci Rep. 2016;6:29343.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Etienne C, Moing A, Dirlewanger E, Raymond P, Monet R, Rothan C. Isolation and characterization of six peach cDNAs encoding key proteins in organic acid metabolism and solute accumulation: involvement in regulating peach fruit acidity. Physiol Plant. 2002;114:259–70.

    Article  CAS  PubMed  Google Scholar 

  9. Degu A, Hatew B, Nunes-Nesi A, Shlizerman L, Zur N, Katz E, Fernie AR, Blumwald E, Sadka A. Inhibition of aconitase in citrus fruit callus results in a metabolic shift towards amino acid biosynthesis. Planta. 2011;234:501–13.

    Article  CAS  PubMed  Google Scholar 

  10. Shimada T, Nakano R, Shulaev V, Sadka A, Blumwald E. Vacuolar citrate/H+ symporter of citrus juice cells. Planta. 2006;224:472–80.

    Article  CAS  PubMed  Google Scholar 

  11. Shi CY, Hussain SB, Yang H, Bai YX, Khan MA, Liu YZ. CsPH8, a P-type proton pump gene, plays a key role in the diversity of citric acid accumulation in citrus fruit. Plant Sci. 2019;289:110288.

    Article  CAS  PubMed  Google Scholar 

  12. Strazzer P, Spelt CE, Li SJ, Bliek M, Federici CT, Roose ML, Koes R, Quattrocchio FM. Hyperacidification of citrus fruit by a vacuolar proton-pumping P-ATPase complex. Nat Commun. 2019;10:744.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Cercós M, Soler G, Iglesias DJ, Gadea J, Forment J, Talón M. Global analysis of gene expression during development and ripening of citrus fruit flesh. A proposed mechanism for citric acid utilization. Plant Mol Biol. 2016;62:513–27.

    Article  Google Scholar 

  14. Chen M, Jiang Q, Yin XR, Lin Q, Chen JY, Allan AC, Xu CJ, Chen KS. Effect of hot air treatment on organic acid- and sugar-metabolism in Ponkan (Citrus reticulata) fruit. Sci Hortic. 2012;147:118–25.

    Article  CAS  Google Scholar 

  15. Li SJ, Yin XR, Xie XL, Allan AC, Ge H, Shen SL, Chen KS. The citrus transcription factor, CitERF13, regulates citric acid accumulation via a protein-protein interaction with the vacuolar proton pump, CitVHA-c4. Sci Rep. 2016;6:20151.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Li SJ, Yin XR, Wang WL, Liu XF, Zhang B, Chen KS. Citrus CitNAC62 cooperates with CitWRKY1 to participate in citric acid degradation via up-regulation of CitAco3. J Exp Bot. 2017;68:3419–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Lu XP, Cao XJ, Li FF, Li J, Xiong J, Long GY, Cao SY, Xie SX. Comparative transcriptome analysis reveals a global insight into molecular processes regulating citrate accumulation in sweet orange (Citrus sinensis). Physiol Plant. 2016;158:463–82.

    Article  CAS  PubMed  Google Scholar 

  18. Li SJ, Wang WL, Ma YC, Liu SC, Grierson D, Yin XR, Chen KS. The Citrus CitERF6 contributes to citric acid degradation via up-regulation of CitAclα1, encoding ATP-citrate lyase subunit α. J Agric Food Chem. 2020;68:10081–7.

    Article  CAS  PubMed  Google Scholar 

  19. He JX, Xu YT, Huang D, Fu JL, Liu Z, Wang L, Zhang Y, Xu RW, Li L, Deng XX, et al. TRIPTYCHON-LIKE regulates aspects of both fruit flavor and color in citrus. J Exp Bot. 2022;73:3610–42.

    Article  CAS  PubMed  Google Scholar 

  20. Wang L, Huang Y, Liu ZA, He JX, Jiang XL, He F, Lu ZH, Yang SZ, Chen P, Yu HW, et al. Somatic variations led to the selection of acidic and acidless orange cultivars. Nat Plants. 2021;7:954–65.

    Article  CAS  PubMed  Google Scholar 

  21. Li XL, Fan JZ, Luo SM, Yin L, Liao HY, Cui XQ, He JZ, Zeng YH, Qu JJ, Bu ZY. Comparative transcriptome analysis identified important genes and regulatory pathways for flower color variation in Paphiopedilum hirsutissimum. BMC Plant Biol. 2021;21:495.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Tanaka Y, Sasaki N, Ohmiya A. Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant J. 2010;54:733–49.

    Article  Google Scholar 

  23. Xu W, Dubos C, Lepiniec L. Transcriptional control of flavonoid biosynthesis by MYB-bHLH-WDR complexes. Trends Plant Sci. 2015;20:176–85.

    Article  CAS  PubMed  Google Scholar 

  24. Yao G, Ming M, Allan AC, Gu C, Li L, Wu X, Wang RZ, Chang YJ, Qi KJ, Zhang SL, et al. Map-based cloning of the pear gene MYB114 identifies an interaction with other transcription factors to coordinately regulate fruit anthocyanin biosynthesis. Plant J. 2017;92:437–51.

    Article  CAS  PubMed  Google Scholar 

  25. Zhu ZX, Wang HL, Wang YT, Guan S, Wang F, Tang JY, Zhang RJ, Xie LL, Lu YQ. Characterization of the cis elements in the proximal promoter regions of the anthocyanin pathway genes reveals a common regulatory logic that governs pathway regulation. J Exp Bot. 2015;66:3775–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. An JP, Zhang XW, Liu YJ, Wang XF, You CX, Hao YJ. ABI5 regulates ABA-induced anthocyanin biosynthesis by modulating the MYB1-bHLH3 complex in apple. J Exp Bot. 2021;72:1460–72.

    Article  CAS  PubMed  Google Scholar 

  27. Ni JB, Premathilake AT, Gao YH, Yu WJ, Tao RY, Teng YW, Bai SL. Ethylene-activated PpERF105 induces the expression of the repressor-type R2R3-MYB gene PpMYB140 to inhibit anthocyanin biosynthesis in red pear fruit. Plant J. 2021;105:167–81.

    Article  CAS  PubMed  Google Scholar 

  28. Qin J, Zhao CZ, Wang SW, Gao N, Wang XX, Na XF, Wang XM, Bi YR. PIF4-PAP1 interaction affects MYB-bHLH-WD40 complex formation and anthocyanin accumulation in Arabidopsis. J Plant Physiol. 2022;268:153558.

    Article  CAS  PubMed  Google Scholar 

  29. Schaart JG, Dubos C, Romero De La Fuente I, Van Houwelingen AM, De Vos RC, Jonker HH, Xu WJ, Routaboul JM, Lepiniec L, Bovy AG. Identification and characterization of MYB-bHLH-WD40 regulatory complexes controlling proanthocyanidin biosynthesis in strawberry (Fragaria × ananassa) fruit. New Phytol. 2012;197:454–67.

    Article  PubMed  Google Scholar 

  30. Butelli E, Garcia-Lor A, Licciardello C, Las Casas G, Hill L, Recupero GR, Keremane ML, Ramadugu C, Krueger R, Xu Q, et al. Changes in anthocyanin production during domestication of citrus. Plant Physiol. 2017;173:2225–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Wang JH, Liu JJ, Chen KL, Li HW, He J, Guan B, He L. Anthocyanin biosynthesis regulation in the fruit of citrus sinensis cv. Tarocco Plant Mol Biol Rep. 2016;34:1043–55.

    Article  CAS  Google Scholar 

  32. Quattrocchio F, Verweij W, Kroon A, Spelt C, Mol J, Koes R. PH4 of Petunia is an R2R3 MYB protein that activates vacuolar acidification through interactions with basic-helix-loop-helix transcription factors of the anthocyanin pathway. Plant Cell. 2006;18:1274–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Butelli E, Licciardello C, Ramadugu C, Durand-Hulak M, Celant A, Recupero GR, Froelicher Y, Martin C. Noemi controls production of flavonoid pigments and fruit acidity and illustrates the domestication routes of modern citrus varieties. Curr Biol Rep. 2019;29:158–64.

    Article  CAS  Google Scholar 

  34. Zhu F, Wen WW, Fernie AR. Finding Noemi: the transcription factor mutations underlying trait differentiation amongst citrus. Trends Plant Sci. 2019;24:384–6.

    Article  CAS  PubMed  Google Scholar 

  35. Ding RR, Che XK, Shen Z, Zhang YH. Metabolome and transcriptome profiling provide insights into green apple peel reveals light- and UV-B-responsive pathway in anthocyanins accumulation. BMC Plant Biol. 2021;21:351.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Peng Z, Wang YQ, Geng GD, Yang R, Yang ZF, Yang CM, Xu RH, Zhang QQ, Kalar KU, Li ZH, et al. Comparative analysis of physiological, enzymatic, and transcriptomic responses revealed mechanisms of salt tolerance and recovery in Tritipyrum. Front Plant Sci. 2022;12:800081.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Carmona L, Sulli M, Diretto G, Alquézar B, Alves M, Peña L. Improvement of antioxidant properties in fruit from two blood and blond orange cultivars by postharvest storage at low temperature. Antioxidants. 2022;11:547.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Habibi F, García-Pastor ME, Puente-Moreno J, Garrido-Auñón F, Serrano M, Valero D. Anthocyanin in blood oranges: a review on postharvest approaches for its enhancement and preservation. Crit Rev Food Sci Nutr. 2022;1–13.

  39. Amato A, Cavallini E, Walker AR, Pezzotti M, Bliek M, Quattrocchio F, Koes R, Ruperti B, Bertini E, Zenoni S, et al. The MYB5-driven MBW complex recruits a WRKY factor to enhance the expression of targets involved in vacuolar hyper-acidification and trafficking in grapevine. Plant J. 2019;99:1220–41.

    Article  CAS  PubMed  Google Scholar 

  40. Cavallini E, Zenoni S, Finezzo L, Guzzo F, Zamboni A, Avesani L, Tornielli GB. Functional diversification of grapevine MYB5a and MYB5b in the control of flavonoid biosynthesis in a petunia anthocyanin regulatory mutant. Plant Cell Physiol. 2014;55:517–34.

    Article  CAS  PubMed  Google Scholar 

  41. Li SJ, Liu XJ, Xie XL, Sun CD, Grierson D, Yin XR, Chen KS. CrMYB73, a PH-like gene, contributes to citric acid accumulation in citrus fruit. Sci Hortic. 2015;197:212–7.

    Article  CAS  Google Scholar 

  42. Lu Q, Liu H, Hong YB, Liang XQ, Li SX, Liu HY, Li HF, Wang RF, Deng QQ, Jiang HF, et al. Genome-wide identification and expression of FAR1 gene family provide insight into pod development in peaunt (Arachis hypogaea). Front in Plant Sci. 2022;13:893278.

    Article  Google Scholar 

  43. Sacco A, Raiola A, Calafiore R, Barone A, Rigano MM. New insights in the control of antioxidants accumulation in tomato by transcriptomic analyses of genotypes exhibiting contrasting levels of fruit metabolites. BMC Genomics. 2019;20:43.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Lai B, Du LN, Wang D, Huang XM, Zhao JT, Wang HC, Hu GB. Characterization of a novel litchi R2R3-MYB transcription factor that involves in anthocyanin biosynthesis and tissue acidification. BMC Plant Biol. 2019;19:62.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Zhang Y, Ye JL, Liu CY, Xu Q, Long LC, Deng XX. Citrus PH4-Noemi regulatory complex is involved in proanthocyanidin biosynthesis via a positive feedback loop. J Exp Bot. 2020;71:1306–21.

    Article  CAS  PubMed  Google Scholar 

  46. He LH, Liu XY, Liu SF, Zhang J, Zhang Y, Sun Y, Tang RM, Wang WB, et al. Transcriptomic and targeted metabolomic analysis identifies genes and metabolites involved in anthocyanin accumulation in tuberous roots of sweetpotato (Ipomoea batatas L.). Plant Physiol Biochem. 2020;156:323–32.

    Article  CAS  PubMed  Google Scholar 

  47. Huang DQ, Zhao YH, Cao MH, Qiao L, Zheng ZL. Integrated systems biology analysis of transcriptomes reveals candidate genes for acidity control in developing fruit of sweet orange (Citrus sinensis L. Osbeck). Front Plant Sci. 2016;7:486.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Rapisarda P, Fanella F, Maccarone E. Reliabbility of analytical methods for determining anthocyanins in blood orange juices. J Agric Food Chem. 2000;48:2249-52.

  49. Pertea M, Pertea GM, Antonescu CM, Chang TC, Mendell JT, Salzberg SL. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotech. 2015;33:290–5.

    Article  CAS  Google Scholar 

  50. Frazee AC, Pertea G, Jaffe AE, Langmead B, Salzberg SL, Jeffrey TL. Ballgown bridges the gap between transcriptome assembly and expression analysis. Nat Biotech. 2015;33:243–6.

    Article  CAS  Google Scholar 

  51. Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28:1947–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51:587–92.

    Article  Google Scholar 

  54. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinf. 2008;9(1):559.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was financially supported by General Fund Project of Hunan Natural Science Foundation (No. 2022JJ30292); the Natural Science Foundation of China (NSFC; No. 31801914); and Scientific Research Project of Education Department of Hunan Province (No. 20A251).

Author information

Authors and Affiliations

Authors

Contributions

LS planned and supervised all aspects of the study. YJ performed the statistical analysis of the transcriptome data and carried out the gene candidate validation strategy. MYL designed the primers and performed the qPCR analysis. NL collected samples for RNA isolation and transcriptome analysis. XQM and HMZ performed citrate and anthocyanin detection. JH, DZL, JFY, GYL and ZND revise the manuscript. LS and XPL wrote the manuscript. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Ling Sheng.

Ethics declarations

Ethics approval and consent to participate

The plant material ‘Tarocco’ blood orange (TBO) and ‘Bingtangcheng’ sweet orange (BTSO) (Citrus sinensis) were purchased from an experimental station under the permission of master of this experimental station in Huaihua, Hunan Province, China, in December 2019. This article did not contain any studies with human participants or animals and did not involve any endangered or protected species. And all experimental studies and experimental materials involved in this research are in full compliance with relevant institutional, national and international guidelines and legislation. All the protocols were performed in accordance with relevant guidelines and regulations of the IUCN Policy Statement on Research Involving Species at Risk of Extinction and the Convention on the Trade in Endangered Species of Wild Fauna and Flora.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Primers used for quantitative real-time PCR analysis. Table S2. Number of reads after filtering rRNA and low quality. Fig. S1. The content of malate in TBO and BTSO fruit during postharvest storage. Error bars represent the standard deviations of the mean in three replicates, and for each storage day, different lowercase letters stand for significant differences between two materials at p < 0.05. Fig. S2. Volcano plots of differentially expressed genes (DEGs) in the comparisons of (A) TBO 0d vs. BTSO 0d, (C) TBO 15d vs. BTSO 15d, (E) TBO 30d vs. BTSO 30d, (G) TBO 60d vs. BTSO 60d, and (I) TBO 90d vs. BTSO 90d. Scatter plot of the KEGG pathway enrichment of DEGs. Rich factor is the ratio of the DEG number to the background number in a certain pathway. The size of the dots represents the number of genes, and the color of the dots represents the range of the q-value. B, D, F, H, and J present the top 20 terms at 0, 15, 30, 60 and 90 DAS between TBO and BTSO, respectively. Fig. S3. KEGG enrichment of the 310 valuable genes. Fig. S4. GO enrichment of the 310 valuable genes. Fig. S5. Correlation analysis of expression profiles of citrate metabolism and anthocyanin biosynthesis related genes, citrate and anthocyanin content. *, ** and *** represent significant differences at p < 0.05, p < 0.01 and p < 0.001, respectively.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Jin, Y., Liao, M., Li, N. et al. Weighted gene coexpression correlation network analysis reveals the potential molecular regulatory mechanism of citrate and anthocyanin accumulation between postharvest ‘Bingtangcheng’ and ‘Tarocco’ blood orange fruit. BMC Plant Biol 23, 296 (2023). https://doi.org/10.1186/s12870-023-04309-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-023-04309-5

Keywords