Skip to main content

Anthocyanins accumulation analysis of correlated genes by metabolome and transcriptome in green and purple peppers (Capsicum annuum)



In order to clarify the the molecular mechanism of anthocyanin accumulation in green and purple fruits of pepper using metabolomics and transcriptomics,to identify different anthocyanin metabolites,and to analyze the differentially expressed genes involved in anthocyanin biosynthesis..


We analyzed the anthocyanin metabolome and transcriptome data of the fruits of 2 purple pepper and 1 green pepper. A total of 5 anthocyanin metabolites and 2224 differentially expressed genes were identified between the green and purple fruits of pepper. Among the 5 anthocyanin metabolites,delphin chloride was unique to purple pepper fruits,which is the mainly responsible for the purple fruit color of pepper. A total of 59 unigenes encoding 7 enzymes were identified as candidate genes involved in anthocyanin biosynthesis in pepper fruit. The six enzymes (PAL,C4H,CHI,DFR,ANS,UFGT) had higher expression levels except the F3H gene in purple compared with green fruits. In addition,seven transcription factors were also found in this study. These transcription factors may contribute to anthocyanin metabolite biosynthesis in the fruits of pepper. One of differentially expressed gene novel.2098 was founded. It was not annotated in NCBI. Though blast analysis we preliminarily considered that this gene related to MYB transcription factor and was involved in anthocyanin biosynthesis in pepper fruit.


Overall, the results of this study provide useful information for understanding anthocyanin accumulation and the molecular mechanism of anthocyanin biosynthesis in peppers.

Peer Review reports


As living standards rise,more and more people begin to pay attention to their nutrition and healthcare. The development and utilization of anthocyanin compounds have become hot spots in the fields of phytochemistry, medicine,and healthcare. Anthocyanins,an important type of water-soluble pigment in plants,belong to the flavonoids and are widely distributed in plant organs,resulting in many colors including red,blue,and purple in plants [1, 2]. Anthocyanins have biological functions that improve health including antioxidation,anticancer,and antiaging functions,protecting eyesight,preventing cardiovascular diseases,and improving memory [3,4,5]. In addition, anthocyanins can reduce the photoinhibition of photosynthesis and photobleaching of chlorophyll under strong light. The accumulation of anthocyanins can increase the photostability of the photosynthetic system in seedlings,reduce plant tissue damage caused by high levels of Ultraviolet (UV) light,and improve stress resistance in plants [6].

Pepper (Capsicum annuum L.) is a plant that belongs to the genus Capsicum of the Solanaceae. There are many types of pepper germplasm resources,and purple pepper is one of the rarer types. Purple pepper is rich in anthocyanins and has good physiological tolerance to high temperature and drought stresses. Generally,plant fruit skin has the highest anthocyanin content. However,the fruit skin of most plants such as eggplant and purple sweet potato is inedible,whereas the edible part of purple pepper is the brightly colored pericarp,which is valuable for human health and directly determines its economic value and popularity in the market.

In recent years,progress has been made in the study of the regulation of anthocyanin biosynthesis and metabolism using genetics,genetic engineering,and molecular biology approaches. The biosynthesis of anthocyanins includes a series of metabolic reactions involving 20 different organic molecules and 12 different catalytic enzymes encoded by multiple homologous genes [7]. The biosynthesis of anthocyanins is an extension of the flavonoid pathway. Phenylalanine lyase (PAL),cinnamic acid 4-hydroxylase(C4H),and coumadin CoA ligase(4CL) are the starting points for the biosynthesis of flavonoids. The successive steps catalysed by the chalcone isomerase (CHI),flavonoid 3-hydroxylase(F3H/FHT),flavonoid 3′-hydroxylase(F3’H),dihydroflavonol-4-reductase (DFR),anthocyanidin synthase (ANS) and 3-glycosyltransferase(3-GT) enzymes lead to the production of the anthocyanin pigments. The late biosynthetic genes (LBGs) F3’H,flavonoid 3′5’-hyroxylase(F3’5’H),DFR,ANS and Anthocyanidin 3-O-glucosyl -transferase (UFGT) are essential for producing anthocyanin or some specific flavonoids [8,9,10].

In higher plants,the biosynthesis of most anthocyanins is often regulated in different spatial temporal patterns by combinations of multiple regulatory factors such as R2R3-MYB,MADS-box,bHLH and WD40 [11, 12]. To date,several R2R3-MYB genes associated with either positive or negative regulation of the key genes (DFR,ANS,and UFGT) in anthocyanin biosynthesis have been identified. However,their contributions to the color of fruits and other organs vary [13,14,15]. bHLH transcription factor is also a key regulator of anthocyanin biosynthesis,usually independently regulating CHS,DFR,and UFGT [16]. WD40 proteins do not bind to the promoters of the anthocyanin biosynthetic genes. Instead,they interact with bHLH and MYB to regulate anthocyanin biosynthesis [17]. The R2R3MYB,bHLH,and WD40 proteins usually form a transcriptional activation complex (MYBbHLH-WD40,MBW) to regulate the transcription of late anthocyanin biosynthetic genes in most plants [18,19,20].

In recent years, Metabolomics and transcriptomics have been widely used to study the biosynthesis of metabolites and their molecular mechanisms [21, 22]. The research on the regulation mechanism of anthocyanin synthesis has been reported in capsicum annuum [23,24,25,26], but different materials of purple pepper contained different anthocyanin metabolites, so the molecular mechanism of anthocyanin biosynthesis was slightly different. Therefore, In this study,2 purple and 1 green pepper (check, CK) were analyzed by means of targeted metabolome and transcriptome technologies to discover the differential genes of enzyme and transcription factors associated with anthocyanin synthesis in peppers, To explore the molecular mechanism of purple pepper fruit formation, provide theoretical and technical support for expanding purple pepper gene resources and promoting purple pepper molecular breeding.


Analysis of pepper anthocyanin metabolome

Anthocyanidins are the most important flavonoid pigments in plants. In order to compare the composition of anthocyanin metabolites between three pepper, L66,L29,and L9(CK),the anthocyanin-like substances in pepper fruits were examined. Principal component analysis (PCA) showed that the samples between the groups scattered and the samples within the groups clustered,indicating that the anthocyanin metabolome data was reliable. On the principal component PC1 and PC2,L9(CK) was clearly separated from L29 and L66 and the difference was significant (Fig. S1). A total of 5 anthocyanin like compounds were detected in the fruits of three pepper. The contents of delphinidin 3-O-glucoside,delphinidin 3-O-rutinoside,and delphin chloride in L66 and L29 were significantly higher than L9(CK). Delphin chloride was found in the fruits of purple pepper L66 and L29,but not in L9(CK). This suggesting that it is a unique metabolite in purple pepper fruits. Delphin chloride could be mainly responsible for the purple fruit color of pepper. The cyanidin 3-O-galactoside content in the fruits of the green pepper line was higher than that in the purple pepper L29 fruits,indicating that the green pepper fruit also accumulated anthocyanins to a certain extent (Fig. 1). The KEGG database was used to annotate the differential metabolites. No annotation result was obtained for cyanidin 3-O-galactoside. The rest of the detected metabolites were found to share the metabolic pathway ko00942.It can be seen from the KEGG pathway map that three of the 4 annotated metabolites were significantly upregulated in purple pepper(L29,L66) and they were delphinidin-derived products.

Fig. 1
figure 1

Schematic diagram of pepper anthocyanin differential metabolite relative content and pathway. Note: The red circle indicates that the metabolites are significantly upregulated; The blue circle indicates that the metabolites are not changed. The histogram shows the relative contents of different metabolites of anthocyanins

Statistics of purple and green pepper transcriptome sequencing data

We further investigated the differences in gene expression among the three samples. With three biological replicates,the transcriptome sequencing of the 9 samples yielded a total of 63.26Gb clean data with 94.14% of bases scoring Q30 (Table S1). The transcriptome sequencing reads were aligned with the reference genome with efficiencies ranged from 86.14 to 95.04% (Table S2), which showed a normal rate of data utilization,suggesting that the selected reference genome was suitable for subsequent analysis.

DEGs between purple and green peppers and KEGG enrichment analysis

In order to clarify DEGs and their biological pathways between purple and green peppers,we analyzed the TrANScriptome sequencing between the fruits of 2 purple pepper (L66 and L29) and 1 green pepper (L9). The number of DEGs were 6567 (L66 vs. L9) and 5091 (L29 vs. L9). A total of 2224 DEGs were common in both of the purple pepper (Table S3). The 2224 DEGs were annotated and 111 KEGG pathways were found. Eight of the top 20 KEGG pathways were shared by the two purple pepper,of which 3 pathways (ko00360,ko00400,and ko00941) were anthocyanin-related (Fig. 2).

Fig. 2
figure 2

Venn diagram of differentially expressed genes (DEGs) and scatter plots of KEGG enrichment in 3 pepper. A Venn diagram of DEGs in 3 pepper. B KEGG enrichment scatter plot of L66 vs. L9. C KEGG enrichment scatter plot of L29 vs. L9.Red arrows indicate: the three KEGG pathways related to anthocyanin. The red arrow from top to bottom in Fig. B:Phenylalanine, tyrosine and tryptophan, Flavonoid biosynthesis, Phenylalanine metabolism; The red arrow from top to bottom in Fig. C: Flavonoid biosynthesis, Phenylalanine, tyrosine and tryptophan, Phenylalanine metabolism

Regulation of anthocyanin biosynthetic pathway genes in purple and green pepper

Anthocyanins are one of the natural products synthesized through the metabolic pathways of phenylpropanes and flavonoids. We analyzed the anthocyanin biosynthesis pathway in order to performed the key genes in the metabolism of purple and green fruits pepper. In total,59 unigenes that encoded 9 enzymes in the anthocyanidin biosynthesis pathways were studied. Fifteen unigenes of the seven enzymes had different expression levels,including 13 upregulated and 2 downregulated unigenes (Fig. 3). The six enzymes (PALC4HCHIDFRANSUFGT) had higher expression levels except the F3H gene in purple compared with green fruits. The 13 unigenes were upregulated by 1.04 to 13.13 fold (Log 2 FC) and 2 unigenes were downregulated by − 1.5 to − 3.5 fold (Table S4).

Fig. 3
figure 3

Heat map of anthocyanin biosynthesis pathway and differential expression of enzyme gene in Pepper. Red color denotes a significant up-regulation. Green color denotes a significant down-regulation

qRT-PCR validation

We examined differential expression of unigenes between purple and green peppers to identify which unigenes had expression correlating with anthocyanin. The expression levels of anthocyanin biosynthetic genes PAL (LOC107871487,LOC107843099,LOC107854291),C4H(LOC107875407,LOC107875406),CHI (LOC107852750), F3H(LOC107843086,LOC107875955),DFR (LOC107850726,LOC107860031,LOC107842348),ANS (LOC107866341),and UFGT (LOC107843659,LOC107861697,LOC107860695) were determined by RT-qPCR using leaves of both purple pepper(L66,L29) and green pepper(L9). The results showed that the expression levels of the 6 genes were significantly higher in fruits of purple pepper than that in fruits of the green pepper line,but The expression of two genes in F3H was higher in green pepper. The qRT-PCR results of the expression of 15 genes were highly consistent with the results obtained based on transcriptome sequencing (Fig. 3), indicating that DEG analysis was highly reliable.

Transcription factors regulating anthocyanin biosynthesis in pepper fruits

Transcription factors (TF) are special genes that modulate the expression levels of other genes and represent the main players for the determination of spatiotemporal transcriptional activity [27]. Among the annotated unigenes,transcription factors related to anthocyanin synthesis were MYB,bHLH,WD40,Bzip,WRKY,MA-DS, Zinc-finger and NAC. In the same transcription factor, DEGs in the two kinds of purple pepper are not identical. All DEGs encoding eight TFs were retrieved and the results revealed that bHLH,MYB and NAC were the more DEGs which may differentially regulate structural genes participating in the Anthocyanin formation of pepper (Table S5). Most of these TFs were found up-regulated in purpel peppers(L66 and L29) when compared to the green pepper(L9),suggesting that a high transcriptional activity is required to form the anthocyanin (Table 1). Eight genes were randomly selected from each transcription factor for QRT-PCR verification,and their expression trends were highly consistent with RNA-SEQ results,indicating that DEG analysis was highly reliable (Fig. 4, Table S6).

Table 1 Transcription factors of anthocyanin biosynthesis in fruit of Capsicum annuum
Fig. 4
figure 4

Expression analysis of 8 transcription factor genes during anthocyanin biosynthesis in purple pepper and green pepper fruits

The gene involved in anthocyanin biosynthesis in fruits of Capsicum annuum

During the biosynthesis of anthocyanin in capsicum annuum,we found a differentially expressed gene novel.2098 that was not annotated in NCBI. Blast analysis indicated that the novel.2098 gene exhibited 86.48,84.75 and 84.24% amino acid identity to NM_001288387.1(Solanum tuberosum),XM_016701512.1 (Capsicum annuum) and XM_015227322.2 (Solanum pennellii). Phylogenetic analysis revealed that the novel.2098 gene had a close relationship with the MYB protein of Solanum tuberosum. Though blast analysis we preliminarily considered that this gene related to MYB transcription factor and was involved in anthocyanin biosynthesis in pepper fruit (Fig. 5).

Fig. 5
figure 5

The homology alignment analysis of novel.2098 in different species. A Blastn analysis; B Phylogenetic tree of novel.2098 in different species. Unknown: represents annotated with novel.2098


Anthocyanin identification in the fruits of Capsicum annuum

Purple pepper had been gaining more popularity because of their antioxidant function,antiaging effects and attractive colors. Anthocyanin as one of the important pigment play a crucial role in the formation of pepper fruit color [28]. The anthocyanin accumulation causes pepper fruit to appear purple and black [29]. According to the studies reported previously,delphinidin is the only anthocyanidin found in peppers [30]. In our study, The accumulation levels of these anthocyanins in particular,delphinidin 3-O-glucoside,delphinidin 3-O-rutinoside and delphin chloride positively correlated with the purple flesh coloration,implying that these compounds partly contribute to the purple color. Among the above three delphinidins,delphin chloride was unique in purple pepper. This result that delphin chloride could be mainly responsible for the purple fruit color of pepper.

Structure genes associated with anthocyanin biosynthesis in fruits of Capsicum annuum

Anthocyanidin biosynthesis belongs to a branch pathway of flavonoids that starts from phenylalanine. After carboxylation,glycosylation,methylation,and acylation,the anthocyanins are transported to and accumulated in vacuoles. A variety of enzymes participate in the biosynthesis of anthocyanins in higher plants [31, 32]. It has been reported that there was no significant difference in the expression levels of PAL,C4H,CHS and 4CL between purple and green peppers,while the expression levels of F3H, F3’5’H,DFR,ANS,and UFGT were all significantly higher in purple pepper than in green peppers [33,34,35].However,it was also suggested that the high expression of PAL, C4H and 4CL in purple pepper might be one of the reasons for the high anthocyanin content [29]. In our study,we found that the expressions of CHS and 4CL did not change significantly in purple pepper fruits, while the expressions of DFR,ANS,CHI,PAL,C4H and UFGT were significantly increased during anthocyanin biosynthesis and the expression level of F3H gene was low in purple pepper fruit, which was slightly different with the results of previous studies. Therefore, we think preliminarily DFR,ANS,CHI,PAL,C4H and UFGT play a positive regulatory role in anthocyanin biosynthesis,and their high expression is conducive to the production of anthocyanin derivatives. However,F3H gene was low expressed in purple pepper,and its decrease did not affect the increase of delphinidin products,and played a role in reverse regulation of anthocyanin biosynthesis. The structural genes of anthocyanin biosynthesis have been studied a lot,but the results are slightly different,which may be caused by different experimental materials and environmental conditions.

Transcription factor associated with anthocyanin biosynthesis in fruits of Capsicum annuum

In addition to structural genes,transcription factors are also the key factors in the synthesis of anthocyanins in plants. To date,transcription factors of MYB,bHLH,WD40,NAC,zinc finger, MADs and WRKY proteins have been found to regulate anthocyanin biosynthesis [29, 36,37,38,39]. In peppers,through qRT-PCR and RNA-Seq analysis, Tang et al. [40] found that two MYBs (CaANT1 and CaANT2),one bHLH (CaAN1),and one WD40(CaTTG1) transcription factor participated in the accumulation of anthocyanins in pepper flowers,and might activate anthocyanin accumulation by forming a new MBW complex. In recent studies,it reported that the WRKY transcription factors interacts with the MBW complex regulate the anthocyanin biosynthetic pathway [34] and some MADS-box proteins have involved in the synthesis of flavonoids [41]. It was speculated that NAC and H2C2 transcription factors may regulate the expression levels of these MYB/HD-Zip genes in the ripening periods [42]. In this study,we found that transcription factors related to anthocyanin synthesis of pepper were MYB,bHLH,WD40,Bzip,WRKY,MADS,Zinc-finger and NAC. Among them,transcription factors bHLH NAC and MYB were the more DEGs which may differentially regulate structural genes participating in the anthocyanin formation of pepper. At present,many transcription factors have been found in the process of anthocyanin biosynthesis in pepper,but the relationship between transcription factors is not clear,which needs further research and verification.

Furthermore,In this study,we detected one differentially expressed MYB genes (novel.2098) that were upregulated in green vs purple pepper. Blast analysis indicated that the novel.2098 gene is a new MYB gene in Capsicum annuum. According to the results of the forefathers [43, 44], we speculated that the high expression of the MYB gene stimulated the upregulation of structural gene expression, which led to the accumulation of anthocyanin derivatives. Therefore,the high expression level of the MYB gene may have stimulated the formation of anthocyanin, leading to the formation of purple fuirts.

By analyzing the transcriptome and metabolome data of purple and green pepper,we found that the biological metabolism of anthocyanin is complex and affected by multiple transcription factors and enzyme genes. Different gene expression levels lead to changes in the accumulation of anthocyanin metabolites,resulting in different degrees of purple in pepper. This experiment studied the metabolic pathways and key genes for the color of purple pepper fruits,which provided new views into the synthesis and accumulation of anthocyanins in pepper fruits and laid a foundation for more effective cultivation of purple pepper. In addition,the biosynthesis of pepper anthocyanins is not only affected by the key enzyme genes and transcription factors,but also affected by various biochemical reactions and environmental conditions. More in-depth research is needed in the future.


Metabolome and transcriptome data were used to reveal the anthocyanin biosynthesis metabolic pathway in pepper fruit. Changes in anthocyanin metabolites,delphinidin 3-O-glucoside and delphinidin 3-O-rutinoside contribute to the purple color. Delphin chloride could be mainly responsible for the purple fruit color of pepper. In an analysis of the anthocyanin biosynthesis pathway,we identified 59 unigenes encoding 7 enzymes that are involved in anthocyanin biosynthesis in fruits of pepper. In addition,seven transcription factors were also found that could contribute to anthocyanin metabolite biosynthesis in the fruits of pepper. Our results provide valuable information on the anthocyanin metabolites and the candidate genes involved in the anthocyanin biosynthesis pathways in pepper.


Plant materials

Two purple pepper L66 and L29 and a green pepper line L9 were provided by the Economic Crop Research Institute of Hebei Academy of Agricultural and Forestry Sciences. Purple pepper line L66 fruits have deeper purple color than L29 fruits (Fig. 6). The pepper were sown in seedling trays in early January and transplanted. Pericarp tissues of the matured fruits collected from different plants were sampled (each line had three repeats),frozen in liquid nitrogen,and then stored at − 80 °C for subsequent metabolite extraction,transcriptome sequencing (RNA-seq),and real-time quantitative PCR (qRT-PCR) analysis.

Fig. 6
figure 6

The phenotype of purple and green pepper fruits

Metabolite extraction and high-performance liquid chromatography (HPLC) conditions

Multiple reaction monitoring (MRM) was performed by Metware Biotechnology (Wuhan,China).Freeze-dried leaves were crushed using a mixer mill (MM 400,Retsch) with a zirconia bead for 1.5 min at 30 Hz.Then,100 mg of the powder was extracted overnight at 4 °C with 1.0 mL 70% aqueous methanol. Following centrifugation at 10,000 g for 10 min,the extracts were absorbed and filtrated (SCAA-104,0.22 μm). The sample extracts were analyzed using an LC-ESI-MS/MS system (HPLC,Shim-pack UFLC Shimadzu CBM30A systemMS,Applied Biosystems 4500QTRAP). The analytical conditions were as follows:high-performance liquid chromatography (HPLC): column, waters ACQUITY UPLC HSS T3 C18 (1.8 μm,2.1 mm*100 mm);solvent system,water (0.04% acetic acid): acetonitrile(0.04% acetic acid);gradient program,100:0 V/V at 0 min,5:95 V/V at 11.0 min,5:95 V/V at 12.0 min,95:5 V/V at 12.1 min,and 95:5 V/V at 15.0 min;flow rate,0.40 mL/min; temperature,40 °C;injection volume:5 μL. The effluent was alternatively connected to an ESI-triple quadrupole-linear ion trap (Q TRAP)-MS [45, 46].

ESI-q trap-MS/MS

Linear ion trap (LIT) and triple quadrupole (QQQ) scANS were acquired on a triple quadrupole-linear ion trap mass spectrometer (Q TRAP), API 6500 Q TRAP LC/MS/MS System,equipped with an ESI Turbo Ion-Spray interface,operating in a positive ion mode and controlled by Analyst 1.6 software (AB Sciex). The ESI source operation parameters were as follows:ion source,turbo spray,with source temperature at 500 °C;ion spray voltage (IS),5500 V; ion source gas I (GSI),gas II (GSII),and curtain gas (CUR) at 55,60,and 25.0 psi,respectively; and the collision gas (CAD) was 12 psi.Instrument tuning and mass calibration were performed with 10 and 100 μmol/L polypropylene glycol solutions in QQQ and LIT modes,respectively. QQQ scans were acquired as MRM experiments with the collision gas (nitrogen) set to 5 psi. Declustering potential (DP) and collision energy (CE) for individual MRM transitions were done with further DP and CE optimization. A specific set of MRM transitions were monitored for each period according to the metabolites eluted within that period.

RNA extraction and Illumina sequencing

Total RNAs for constructing nine cDNA libraries of the three pepper,each with three replicates,were isolated from frozen pericarp tissues using the modified CTAB method [47]. mRNAs were used as templates for reverse transcription into cDNAs that were then used to construct cDNA libraries after purification. Agilent 2100 bioanalyzer was used to detect the insert size to control the quality of the libraries. The qualified libraries were sequenced on the Illumina HiSeq 4000 platform. The reads of adapters, N-containing sequences (N indicates ambiguous bases),and low-quality sequences (reads with the bases that have Qphred 20 or less accounting for over 50% of the entire read length) were removed. HISAT2 software was used to quickly and accurately align the clean reads with the reference genome ( to locate the reads on it. FPKM (Fragments Per Kilobase of transcript per Million fragments mapped) was used to correct the sequencing depth and gene length [48]. DESeq software [49] was employed to analyze and screen the genes differentially expressed between sample groups,adopting |log2(Fold Change)| > 1 and padj< 0.05 as the screening standard. The differentially expressed genes (DEGs) were then analyzed through GO enrichment and KEGG [50] functional annotation of metabolic pathways.

Verification of DEGs by qRT-PCR

Real-time fluorescence qRT-PCR was used to verify the results obtained based on transcriptome sequencing. Primer Premier 5.0 software (Primer,Canada) was used to design specific primers for DEGs. The PCR program was as follows: 95 °C for 5 min followed by 40 cycles of 95 °C for 15 s and 60 °C for 30 s. The PCR was repeated 3 times for each sample. Actin (serial number: GQ339766) was used as the internal reference gene to determine the relative expression level of DEGs.

Availability of data and materials

All data generated or analysed during this study are included in this published article and its Supplementary information files. The raw RNA-seq data are freely available at: PRJNA827972.





Dihydroflavonol 4-reductase


Anthocyanidin synthase


Anthocyanidin 3-O-glucosyltransferase


Chalcone synthase


Chalcone isomerase


Flavonone 3-hydroxylase


Flavonoid 3′-monooxygenase


Flavonoid 3′,5′-hydroxylase


Anthocyanidin 3-O-glucoside 5-O-glucosyltransferase


Basic helix-loop-helix


Differentially expressed genes


Ultra- performance liquid chromatography


Linearion hydrazine-flight time


Triple quadrupole


Gas I


Gas II


Curtain gas


Collision gas


Multiple reaction monitoring


Declustering potential


Collision energy


Principal component analysis


Kyoto Encyclopedia of Gene and Genome


National Center for Biotechnology Informatio


Gene Ontology


Quantitative real-time polymerase chain reaction


Unit variance scaling


  1. Liu Y, Tikunov Y, Schouten R, et al. Anthocyanin biosynthesis and degradation mechanisms in solanaceous vegetables: a review. Front Chem. 2018;6:52.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. Naing A, Kim C. Roles of R2R3-MYB transcription factors in transcriptional regulation of anthocyanin biosynthesis in horticultural plants. Plant Mol Biol. 2018;98:1–18.

    Article  CAS  PubMed  Google Scholar 

  3. Marszalek K, Wozniak L, Kruszewski B, Skapaka S. The effect of high pressure techniques on the stability of anthocyanins in fruit and vegetables. Int J Mol Sci. 2017;18(2):277.

    Article  PubMed Central  CAS  Google Scholar 

  4. Yousuf B, Gul K, Wani A, Singh P. Health benefits of anthocyanins and their encapsulation for potential use in food systems. Crit Rev Food Sci Nutr. 2016;56(13):2223–30.

    Article  CAS  PubMed  Google Scholar 

  5. Wong MK, Takei Y. Lack of plasma kallikrein-kinin system cascade in teleosts. PLoS One. 2013;8(11):81057.

    Article  Google Scholar 

  6. Nakabayashi R, Yonekura Sakakibara K, Urano K, Suzuki M, Yamada Y, Nishizawa T, et al. Enhancement of oxidative and drought tolerance in Arabidopsis by overaccumulation of antioxidant flavonoids. Plant J. 2014;77:367–79.

    Article  CAS  PubMed  Google Scholar 

  7. Nakatsuka T, Nishihara M, Mishiba K, et al. Two different mutations are involved in the formation of white flowered gentian plants. Plant Sci. 2005;16(9):949–58.

    Article  CAS  Google Scholar 

  8. Aza-Gonzalez C, Herrera-Isidron L, Núñez-Palenius HG, Martínez De La Vega O, Ochoa-Alejo N. Anthocyanin accumulation and expression analysis of biosynthesis-related genes during chili pepper fruit development. Biol Plant. 2013;57:49–55.

    Article  CAS  Google Scholar 

  9. Zhang Z, Lid W, Jin H, Yin YX, Zhang HX, Chai WG, et al. VIGS approach reveals the modulation of anthocyanin biosynthetic genes by CaMYB in chili pepper leaves. Front Plant Sci. 2015;6:500–10.

    PubMed  PubMed Central  Google Scholar 

  10. Passeri V, Koes R, Quattrocchilio FM. New challenges for the Design of High Value Plant Products: stabilization of Anthocyanins in plant vacuoles. Front Plant Sci. 2016;7:153–62.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Hichri I, Barrieu F, Bogs J, Kappel C, Delrot S, Lauvergeat V. Recent advances in the transcriptional regulationof the flavonoid biosynthetic pathway. J Exp Bot. 2011;62(8):2465–83.

    Article  CAS  PubMed  Google Scholar 

  12. Brounp, et al. Transcriptional control of flavonoid biosynthesis: a complex network of conserved regulators involved in multiple aspects of differentiation in Arabidopsis. Curr Opin Plant Biol. 2005;8(3):272–9.

    Article  CAS  Google Scholar 

  13. Kobayashi S, Yamamoto NG, Hirochika H. Association of VvmybA1 gene expression with anthocyanin production in grape (Vitis vinifera) skincolor mutants. Hortic Sci. 2016;74:196–203.

    Google Scholar 

  14. Strygina KV, Kochetov AV, Khlestkina EK. Genetic control of anthocyanin pigmentation of potato tissues. BMC Genet. 2019;20(1):27–31.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Niu SS, Xu CJ, Zhang WS, Bo Z, Xian L, Linwang K, et al. Coordinated regulation of anthocyanin biosynthesis in chinese bayberry (Myricarubra) fruit by a R2R3-MYB transcription factor. Planta. 2010;231(4):887–99.

    Article  CAS  PubMed  Google Scholar 

  16. Stommel JR, Lightbourn GJ, Winkel BS, Griesbach RJ. Transcription factor families regulate the anthocyanin biosynthetic pathway in Capsicum annuum. J Hortic Sci. 2009;134:244–51.

    Google Scholar 

  17. Dumm JM, Stommel JR. Coordinated regulation of biosynthetic and regulatory genes coincides with AnthocyaninAccumulation in Developing Eggplant Fruit. Hortic Sci. 2015;140:129–35.

    Google Scholar 

  18. Zhao L, Gao L, Wang H, Chen X, Wang Y, Yang H, et al. The R2R3-MYB,bHLH,WD40,and related transcription factors in flavonoid biosynthesis. Funct Integr Genomics. 2013;13(1):75–82.

    Article  CAS  PubMed  Google Scholar 

  19. Docimo T, Francese G, Ruggiero A, Batelli G, Depalma M, Bassolino L, et al. Phenylpropanoids accumulation in eggplant fruit:characterization of biosynthetic genes and regulation by a MYB transcription factor. Front.Plant. 2016;6:1233–51.

    Google Scholar 

  20. Montefiori M, Brendolise C, Dare AP, Lin-Wang K, Davies KM, Hellens RP, et al. In the Solanaceae,a hierarchy of bHLHs confer distinct target specificity to the anthocyanin regulatorycomplex. J Exp Bot. 2015;66:1427–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Jiang T, Zhang M, Wen C, Xie X, Tian W, Wen S, et al. Integrated metabolomic and transcriptomic analysis of the anthocyanin regulatory networks in Salvia miltiorrhiza BMC Plant Biol. 2020;20:349.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Cai X, Lin L, Wang X, Xu C, Wang Q. Higher anthocyanin accumulation associated with higher transcription levels of anthocyanin biosynthesis genes in spinach. Genome. 2018;61(7):487–96.

    Article  CAS  PubMed  Google Scholar 

  23. Jung S, Venkatesh J, Kang MY, Kwon JK, Kang BC. A non-LTR retrotransposon activates anthocyanin biosynthesis by regulating a MYB transcription factor in Capsicum annuum. Plant Sci. 2019;287:110181.

    Article  CAS  PubMed  Google Scholar 

  24. Liu J, Ai X, Wang Y, Lu Q, Li T, Wu L, et al. Fine mapping of the Ca3GT gene controlling anthocyanin biosynthesis in mature unripe fruit of Capsicum annuum L. Theor Appl Genet. 2020;133(9):2729–42.

    Article  CAS  PubMed  Google Scholar 

  25. Guoab Y, Baiab J, Duanab X, Wangab J. Accumulation characteristics of carotenoids and adaptive fruit color variation in ornamental pepper. Sci Hortic. 2021;275:109699.

    Article  CAS  Google Scholar 

  26. Villa-Rivera MG, Ochoa-Alejo N. Transcriptional Regulation of Ripening in Chili Pepper Fruits ( Capsicum spp.). Int J Mol Sci. 2021;22(22):12151.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Banerjee N, Zhang MQ. Identifying cooperativity among transcription factors controlling the cell cycle in yeast. Nucleic Acids Res. 2003;31:7024–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Lightbourn GJ, Griesbach RJ, Novotny JA, Clevidence BA, Rao DD, Stommel JR. Effects of anthocyanin and carotenoid combinations on foliage and immature fruit color of Capsicum annuum L. J Hered. 2008;99(2):105–11.

    Article  CAS  PubMed  Google Scholar 

  29. Liu Y, Lv J, Liu Z, Wang J, Yang B, Chen W, et al. Integrative analysis of metabolome and transcriptome reveals the mechanism of color formation in pepper fruit (Capsicum annuum L.). Food Chem. 2020;306:125629.

    Article  CAS  PubMed  Google Scholar 

  30. Aza-Gonzalez C, Nuez-Palenius HG, Ochoa-Alejo N. Molecular biology of chili pepper anthocyanin biosynthesis. Chem Sci. 2012;56:93–8.

    CAS  Google Scholar 

  31. Zhang Q, Su LJ, Chen JW, et al. The antioxidative role of anthocyanins in Arabidopsis under high irradiance. Biol Plant. 2012;56(1):97–104.

    Article  CAS  Google Scholar 

  32. Tingrui MA, Wen J, Guang ZH, et al. Anthocyanin synthesis and gene regulation. Agricult Sci Technol. 2012;13(3):507–11.

    Google Scholar 

  33. Zhang H, Zhao X, Zhang J, Yang B, Yu Y, Liu T, et al. Functional analysis of an anthocyanin synthase gene StANS in potato. Sci Hortic. 2020;272:109569.

    Article  CAS  Google Scholar 

  34. Borovsky M, Oren-Shamir R, Ovadia WD, Jong I, et al. The a locus thatcontrols anthocyanin accumulation in pepper encodes a MYB transcription factorhomologous to Anthocyanin2 of Petunia. Theor Appl Genet. 2004;109:23–9.

    Article  CAS  PubMed  Google Scholar 

  35. Jung S, Venkatesh J, Kang M-Y, Kwon J-K, Kang B-C. A non-LTR retrotransposon activates anthocyanin biosynthesis by regulatinga MYB transcription factor in Capsicum annuum. Plant Sci. 2019;287:110181.

    Article  CAS  PubMed  Google Scholar 

  36. Gonzalez CA, Isidron LH, Palenius HG, et al. Anthocyanin accumulation and expression analysis of biosynthesis related genes during chill pepper fruit development. Biol Plant. 2012;99:132–8.

    Google Scholar 

  37. Schaart JG, Dubos C, Romero I, et al. Bovy. Identification and characterization of MYB-bHLH-WD40 regulatory complexes controlling proanthocyanidin biosynthesis in strawberry (Fragaria 9 ananassa) fruits. New Phytol. 2013;197:454–67.

    Article  CAS  PubMed  Google Scholar 

  38. Lloyd A, Brockman A, Aguirre L, Campbell A, Bean A, Cantero A, et al. Advances in the MYB-bHLH-WD repeat (MBW) pigment regulatory model: addition of a WRKY factor and cooption of an anthocyanin MYB for Betalain regulation. Plant Cell Physiol. 2017;58:1431–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Xin N, Yang Z, Li Z, Yi T, Wang J, Zhiping Z. Mechanisms of MYB-bHLH-WD40 complex in the Regu-lation of anthocyanin Biosynthes. Agricult Biotechnol. 2015;4(3):5–8,13.

    CAS  Google Scholar 

  40. Tang B, Li L, Hu Z, Chen Y, Tan T, Jia Y, et al. Anthocyanin accumulation and transcriptional regulation of anthocyanin biosynthesis in purple pepper. J Agric Food Chem. 2020;68(43):12152–63.

    Article  CAS  PubMed  Google Scholar 

  41. Nesi N, Debeaujon I, Jond C, Stewart AJ, Jenkins GI, Caboche M, et al. The transparent testa16 locus encodes the arabidopsis bsister MADS domain protein and is required for proper development and pigmentation of the seed coat. Plant Cell. 2002;14(10):2463–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Zhang Q, Wang L, Liu Z, Zhao Z, Zhao J, Wang Z, et al. Transcriptome and metabolome profiling unveil the mechanisms of Ziziphus jujuba Mill.peel coloration. Food Chem. 2020;312:125903.

    Article  CAS  PubMed  Google Scholar 

  43. Gomez-Maldonado J, Avila C, et al. Functional interactions between a glutamine synthetase promoter and MYB proteins. Plant J. 2004;39:513–26.

    Article  CAS  PubMed  Google Scholar 

  44. Duthie G, Crozier A. Plant-derived phenolic antioxidants. Curr Opin Lipidol. 2000;11:43–7.

    Article  CAS  PubMed  Google Scholar 

  45. Zhang S, Ying H, Pingcuo G, Wang S, Zhao F, Cui Y, Zeng X. Identification of potential metabolites mediating bird’s selective feeding on Prunus mira flowers. Biomed Res Int. 2019;(2019):1–8.

  46. Zhang X, Jiang X, He Y, Li L, Xu P, Sun Z, Li J, Xu J, Xia T, Hong G. AtHB2 a class II HD-ZIP protein, negatively regulates the expression of Csans which encodes a key enzyme in Camellia sinensis catechin biosynthesis. Physiol Plant. 2019;166(4):936–45.

    Article  CAS  PubMed  Google Scholar 

  47. Chai GH, Qi G, Cao YP, Wang ZG, Li Y, Tang XF, et al. Poplar PdC3H17 and PdC3H18 are direct targets of PdMYB3 and PdMYB21 and positively regulate secondary wall formation in Arabidopsis and poplar. New Phytol. 2014;203:520–34.

    Article  CAS  PubMed  Google Scholar 

  48. Bray N, Pimentel H, Melsted P, et al. Near-optimal RNA-Seq quantification. Xiv preprint. 2015;1505:02710.

    Google Scholar 

  49. Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51.

    Article  CAS  PubMed  Google Scholar 

Download references


We thank LetPub ( for its linguistic assistance during the preparation of this manuscript; and thank the staff of Wuhan Metware Biotechnology Co., Ltd. (Wuhan, China) for their support during the metabolite data analysis.


This work was supported by Hebei Key Research and Development Projects (No.20326326D) and the Modern Agricultural Industrial Technology System Project of Hebei Province (HBCT2018030204).

Author information

Authors and Affiliations



Libin Yan and Yanqin Fan designed the project. Yaning Meng analyzed the data, carried out the experiments, and wrote the manuscript. Hongxiao Zhang participated in the design of experiments. All authors contributed to the revision of this manuscript and approved the final manuscript.

Corresponding authors

Correspondence to Yanqin Fan or Libin Yan.

Ethics declarations

Ethics approval and consent to participate

All plant materials (not endangered materials or species) comply with local institutional guidelines and legislation.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have 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: Fig S1.

Principal component analysis (PCA) of two purple (L66,L29) and one green (L9) pepper fruits.

Additional file 2: Table s1.

Summary of sample sequencing data quality.

Additional file 3: Table s2.

Sample and reference genome comparison statistics.

Additional file 4: Table S3.

Results of gene expression in purple vs green peppers.

Additional file 5: Table S4.

Differentially expressed 15 structural genes of anthocyanin biosynthesis in purple vs green peppers.

Additional file 6: Table S5.

Differentially expressed 8 transcription factors of anthocyanin biosynthesis in purple vs green peppers.

Additional file 7: Table S6.

Differentially expressed genes of transcription factors in purple vs green peppers.

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 The Creative Commons Public Domain Dedication waiver ( 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

Meng, Y., Zhang, H., Fan, Y. et al. Anthocyanins accumulation analysis of correlated genes by metabolome and transcriptome in green and purple peppers (Capsicum annuum). BMC Plant Biol 22, 358 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: