Transcriptome, Plant Hormone, and Metabolome Analysis Reveals the Mechanism of Purple Pericarp Formation in ‘Zihui’ Papaya (Carica papaya L.)

The color of the pericarp is a crucial characteristic that influences the marketability of papaya fruit. Prior to ripening, normal papaya exhibits a green pericarp, whereas the cultivar ‘Zihui’ displays purple ring spots on the fruit tip, which significantly affects the fruit’s visual appeal. To understand the mechanism behind the formation of purple pericarp, this study performed a thorough examination of the transcriptome, plant hormone, and metabolome. Based on the UPLC-ESI-MS/MS system, a total of 35 anthocyanins and 11 plant hormones were identified, with 27 anthocyanins and two plant hormones exhibiting higher levels of abundance in the purple pericarp. In the purple pericarp, 14 anthocyanin synthesis genes were up-regulated, including CHS, CHI, F3H, F3′5′H, F3′H, ANS, OMT, and CYP73A. Additionally, through co-expression network analysis, three MYBs were identified as potential key regulators of anthocyanin synthesis by controlling genes encoding anthocyanin biosynthesis. As a result, we have identified numerous key genes involved in anthocyanin synthesis and developed new insights into how the purple pericarp of papaya is formed.


Introduction
Papaya (Carica papaya L.) is a perennial fleshy herb and a significant fruit-bearing tree in tropical and subtropical regions.It is popular worldwide, being rich in vitamins, protein, catabolic enzymes, and minerals [1].An important breeding trait is the appearance of papaya pericarp, which is crucial to consumer evaluation.The primary objective of breeding is to select a papaya pericarp's external traits to meet public expectations [2].
The vibrant colors seen in plant flowers, fruits, and leaves are often a result of the buildup of pigment substances [3].Natural plant pigments can generally be categorized into four groups: chlorophyll, carotenoids, anthocyanins, and betalains.Chlorophyll is a critical pigment for photosynthesis and is widely distributed in plants, giving plant organs a dark or light green appearance based on chlorophyll content differences [4].Natural terpenoid lipid-soluble pigments called carotenoids show mainly as yellow, orange, and red colors [5,6].Plants are widely distributed with anthocyanins, water-soluble pigments that are known for their antioxidant properties and present in blue, purple, and red colors [7,8].Betalains are primarily found in sugar beets, distinguishing them from the aforementioned three pigments [9].The presence of these plant pigments contributes to the vibrant and colorful appearance of plants [10].
Anthocyanins are the most significant coloring compounds in flavonoids.Up to now, anthocyanins have been successfully isolated and identified in more than 600 plant types [11].The mechanism of anthocyanins that allows multiple colors is complicated, mainly involving three aspects: First, the type and concentration of pigments are the primary factors influencing color variation, anthocyanin pigments are found in a wide variety of species and present a variety of colors.The most basic types are pelargonidin (orange), cyanidin (red), and delphinidin (purple).These are further modified through methylation, glycosylation, and acylation to produce different types of anthocyanin glycosides [12].Pigments of different colors interact to produce a spectrum of colors, and the concentration of pigment directly affects the final color display.For example, pelargonidin typically appears orange, but when accumulated in large amounts, it tends to appear purple-red or black-red, while low concentrations result in a pink hue [3].Second, altering anthocyanins in different ways will dramatically alter how they render color.When anthocyanins are methylated, the color will clearly change to red [13].An important stage in the molecular alteration of anthocyanins is glycosylation; the quantity and location of this modification greatly affect how the color appears.They will have a tendency to turn bluish violet due to glycosylated anthocyanins [14].In addition to glycosylation location, acylation modification is also influenced by the number and position of acyl groups.For instance, polyaromatic acyl modification at positions 3 and 7 can cause the plant to appear blue, while acyl modification at positions 3 and 5 results in a purplish-red color [15].Furthermore, anthocyanin's color and physicochemical characteristics can be altered by processes such as fatty acid acylation and aromatic acid acylation [14].Third, the hue of anthocyanins has also been shown to be significantly affected by environmental factors including illumination, temperature, humidity, nutrition, pH, and hormonal stimulation [16,17].
'Zihui' is a newly developed papaya variety by the Fruit Tree Research Institute of the Guangdong Academy of Agricultural Sciences.It develops a well-established root system, produces high yields, and has good stress resistance, especially in high-temperature environments, where its continuous bearing ability is strong.In comparison to conventional varieties, 'Zihui' shows great potential for widespread adoption and use, but the presence of purple ring spots on the fruit tip can negatively impact consumer preference.Therefore, in our study, transcriptome, plant hormones, and metabolome analysis were used to identify key genes involved in anthocyanin synthesis and to develop new insights into how the purple pericarp of papaya is formed.

Sequencing Data Analysis of 12 RNA Libraries from Different Colored Papaya Pericarp
RNA sequencing of the YG, YP, MG, and MP generated approximately 71.74 GB of nucleotide data (Figure 1).We obtained 38,351,072,194 reads for the samples, and 85.83 to 88.60% were mapped to the papaya reference genome.Scores greater than Q30 were found in over 93.73% of the clean data, and GC content ranged from 44.28 to 45.23%.The detailed information about RNA sequencing and mapping is summarized in Table 1.

Sequencing Data Analysis of 12 RNA Libraries from Different Colored Papaya Pericarp
RNA sequencing of the YG, YP, MG, and MP generated approximately 71.74 GB of nucleotide data (Figure 1).We obtained 38,351,072,194 reads for the samples, and 85.83 to 88.60% were mapped to the papaya reference genome.Scores greater than Q30 were found in over 93.73% of the clean data, and GC content ranged from 44.28 to 45.23%.The detailed information about RNA sequencing and mapping is summarized in Table 1.Between the three biological replications, a high correlation was observed, indicating a steady expression pattern in the replicated samples (Figure S1).In total, 3382 DEGs were discovered in the comparison of YP&YG, YP&MP, YG&MG, and MP&MG.These genes included 773 (561 up-regulated and 212 down-regulated) and 650 DEGs (342 up-regulated and 308 down-regulated) between the YP&YG and MP&MG, respectively (Figure 2).
The DEGs were performed with KEGG enrichment to better understand their biological functions; DEGs between YP&YG were mainly located in plant hormone signal transduction, plant-pathogen interaction, MAPK signaling pathway-plant, starch and sucrose metabolism, pentose and glucuronate interconversions pathways, etc. DEGs between MP&MG were mainly classified into plant hormone signal transduction, plant-  Between the three biological replications, a high correlation was observed, indicating a steady expression pattern in the replicated samples (Figure S1).In total, 3382 DEGs were discovered in the comparison of YP&YG, YP&MP, YG&MG, and MP&MG.These genes included 773 (561 up-regulated and 212 down-regulated) and 650 DEGs (342 up-regulated and 308 down-regulated) between the YP&YG and MP&MG, respectively (Figure 2).
The DEGs were performed with KEGG enrichment to better understand their biological functions; DEGs between YP&YG were mainly located in plant hormone signal transduction, plant-pathogen interaction, MAPK signaling pathway-plant, starch and sucrose metabolism, pentose and glucuronate interconversions pathways, etc. DEGs between MP&MG were mainly classified into plant hormone signal transduction, plant-pathogen interaction, flavonoid biosynthesis, ABC transporters, diterpenoid biosynthesis pathways, etc.Among them, DEGs between YP&YG and MP&MG were commonly enriched in plant hormone signal transduction, plant-pathogen interaction, ABC transporters, and flavonoid biosynthesis pathways, which indicated that those pathways may contribute to the formation of purple ring spots on the pericarp of papaya fruit tips (Figure 3).pathogen interaction, flavonoid biosynthesis, ABC transporters, diterpenoid biosynthesis pathways, etc.Among them, DEGs between YP&YG and MP&MG were commonly enriched in plant hormone signal transduction, plant-pathogen interaction, ABC transporters, and flavonoid biosynthesis pathways, which indicated that those pathways may contribute to the formation of purple ring spots on the pericarp of papaya fruit tips (Figure 3).pathogen interaction, flavonoid biosynthesis, ABC transporters, diterpenoid biosynthesis pathways, etc.Among them, DEGs between YP&YG and MP&MG were commonly enriched in plant hormone signal transduction, plant-pathogen interaction, ABC transporters, and flavonoid biosynthesis pathways, which indicated that those pathways may contribute to the formation of purple ring spots on the pericarp of papaya fruit tips (Figure 3).

Different Concentrations of Plant Hormones between Purple and Green Papaya Pericarp
To determine which plant hormones are responsible for the formation of purple ring spots on the pericarp of papaya fruit tips, a UHPLC-MRM-MS method was used to detect 24 target metabolites of plant hormones between MP and MG.Among the 24 target plant hormones, 11 were detected.Only four plant hormones, abscisic acid, 1-Aminocyclopropanecarboxylic acid, N-((-)-jasmonoyl)-S-isoleucine, and trans-Zeatin, were differentially expressed between the MP and MG (Table 2).1-Aminocyclopropanecarboxylic acid and N-((-)-jasmonoyl)-S-isoleucine were up-regulated, while abscisic acid and trans-Zeatin were down-regulated in the MG pericarp.With RNA sequencing, there were 24 DEGs involved in signal transduction in plant hormones according to KEGG enrichment analysis.Among them, 15 DEGs were up-regulated, while nine DEGs were down-regulated in the MP group.The expression profiles of these GEGs are shown in Figure S2.

Metabolome Profiling of the Papaya Pericarp Related to Anthocyanin Metabolism
To study the relationship between anthocyanin metabolism and the formation of purple ring spots on papaya pericarp, UPLC-MS-MS was used to detect 108 anthocyanins based on the public metabolic database, and 35 anthocyanins were detected between the MP and MG papaya pericarp.Anthocyanin derivative accumulations at the metabolomic level were responsive to color differences in papaya pericarp based on the PCA analysis (97.64%) (Figure 4a).Biological repeat correlation was evaluated using the Spearman rank correlation coefficient r.The closer r 2 is to 1, the stronger the correlation between two duplicate samples.In this study, the intra-group correlation coefficients of the two samples were greater than 0.85, indicating the differential accumulation metabolites (DAMs) obtained were reliable (Figure 4b).

Integrated Analysis of DAMs and DEGs in Response to Different Colors of Papaya Pericarp
Analysis of combined transcriptomic/metabolomic regulatory networks was conducted to explore the differences in anthocyanin biosynthesis between MP and MG papaya pericarp.Flavonoid biosynthesis appeared to be the DEGs' active biological pathway based on KEGG enrichment analysis.Fifteen DEGs were found on the flavonoid pathway, with two CHS (evm.TU.supercontig_70.51,evm.TU.supercontig_59.102),two CHI (evm.TU.supercontig_19.60,evm.TU.supercontig_1855.1),one F3H (evm.TU.supercon-

Correlation Analysis between Transcripts and Anthocyanidins
It was still not clarified how those genes control the synthesis of anthocyanins in C papaya.A co-expression network analysis of DEGs and DAMs was conducted to investigate the regulatory network of anthocyanidin biosynthesis in the different colored papaya pericarps.Twenty genes that had a high correlation (r > 0.9) with 10 anthocyanins, including nine structural genes, were identified: three OMT, three cytochrome P450 genes (CYP), one ANS, one DFR, one F3′5′H, four transcription factors (three MYBs and one SRF), and one WAT-related protein, one glutathione S-transferase F12 (GSTF12), one alpha carbonic
anhydrase 1 (ACA1), one expansion-like B1 (EXLB1), one ribulose bisphosphate carboxylase small chain (RBCS), and two unknown proteins (Figure 6 and Table S3).MYB transcription factors play a crucial role in the synthesis of anthocyanins, and previous studies found that glutathione S-transferase was involved in flavonoid transport [29].In our study, two MYB transcription factors, evm.TU.supercontig_119.55  petunidin-3-O-sambubioside-5-O-glucoside), respectively.Evm.TU.supercontig_2742.1 (MYB113) was highly negatively related to cyanidin-3-Osophoroside.The WAT-related protein, GSTF12, ACA1, EXLB1, and RBCS in this study were reported to be involved in anthocyanin biosynthesis in C papaya for the first time.In summary, the above genes may play an important role in anthocyanin biosynthesis.To validate the transcriptome sequencing results, 10 genes were selected for qRT-PCR validation.As expected, these results coincided with transcriptome sequencing results (Figure 7).

Discussion
Papaya is a commercially popular tropical fruit tree with fruit rich in various vitamins and minerals considered beneficial for human health.Under normal conditions, the papaya fruit pericarp is green before ripening and orange-yellow when matured.In our study, the cultivar 'Zhihui' is special, with a circle of purple ring spots on the fruit tip that seriously affects its appearance.To date, there are few research reports on the color of the papaya fruit pericarp.One reason for this is the lack of diversity in research materials.This special material as the object of study can provide new insight into how the circle of purple ring spots on the fruit tip pericarp is formed [30].
Plant hormones are key intrinsic environmental factors that play a prominent role in the regulation of plant anthocyanin signaling pathways [36].Previous studies showed that low auxin levels can promote the expression of anthocyanin-related genes, thus promoting the accumulation of anthocyanins, while a high concentration of auxin inhibits it

Discussion
Papaya is a commercially popular tropical fruit tree with fruit rich in various vitamins and minerals considered beneficial for human health.Under normal conditions, the papaya fruit pericarp is green before ripening and orange-yellow when matured.In our study, the cultivar 'Zhihui' is special, with a circle of purple ring spots on the fruit tip that seriously affects its appearance.To date, there are few research reports on the color of the papaya fruit pericarp.One reason for this is the lack of diversity in research materials.This special material as the object of study can provide new insight into how the circle of purple ring spots on the fruit tip pericarp is formed [30].
Plant hormones are key intrinsic environmental factors that play a prominent role in the regulation of plant anthocyanin signaling pathways [36].Previous studies showed that low auxin levels can promote the expression of anthocyanin-related genes, thus promoting the accumulation of anthocyanins, while a high concentration of auxin inhibits it [37,38].
On the effect of ethylene treatment on photoinduced anthocyanin synthesis in sorghum, it was found that ethylene could promote anthocyanin accumulation in the early lagging stage but inhibit it in the late lagging stage [39].Those results indicated that both auxin and ethylene have positive and negative regulation roles in the anthocyanin signal pathway; auxin is presented as an enhancement at low concentrations but is inhibitory at high concentrations.Anthocyanin accumulation has been demonstrated in rice, ornamental cabbage, and wheat after ABA treatment [40][41][42].Cytokinin (CTK) can enhance the induction effect of light and can promote anthocyanin accumulation in the aleurone layer.This was also confirmed in rice and Arabidopsis thaliana [43,44].It was found that salicylic acid (SA) could induce the expression of the ubiquitin-conjugating enzyme (UBC) OgUBC1 gene, and over-expression of OgUBC1 in A. thaliana increased anthocyanin accumulation and resistance to UV radiation and pathogens [45].Through biological or abiotic stress stimulation, the level of endogenous jasmonic acid (JA) in plants was increased, and the genes related to JA and the anthocyanin synthesis signal were significantly up-regulated, resulting in anthocyanin accumulation [46].The addition of exogenous gibberellin (GA) could significantly inhibit the accumulation of anthocyanins, while inhibition of gibberellin synthesis can increase the anthocyanin content in mature embryos [47].In conclusion, plant hormones play different roles in regulating the anthocyanin signaling pathway in plants.Auxin and ethylene have both positive and negative regulatory effects; ABA, CTK, SA, and JA have positive regulatory effects, and GA has negative regulatory effects.
In our study, 11 plant hormones were detected by UHPLC-MRM-MS.Only four plant hormones, abscisic acid, 1-Aminocyclopropanecarboxylic acid, N-((-)-jasmonoyl)-S-isoleucine, and trans-Zeatin were differentially expressed between the MP and MG.Abscisic acid and trans-Zeatin were up-regulated in the MP pericarp, consistent with the result that ABA and CTK could promote the accumulation of anthocyanins.1-Aminocyclopropanecarboxylic acid and N-((-)-jasmonoyl)-S-isoleucine were down-regulated in the MP pericarp, differently from the others, indicating that the pathways of plant hormones regulating anthocyanin synthesis vary among different plants.The analysis of its regulatory pathway and regulatory network requires further study.
In previous studies, two types of genes were reportedly associated with anthocyanidin biosynthesis: one type encodes enzymes (structural genes) that catalyze the synthesis of different anthocyanin glycosides in the flavonoid synthesis pathway, and the other type gene is regulatory and mainly includes transcription factors [48].By combining transcriptome and metabolome profiling analyses, 15 structural DEGs related to the pathway of anthocyanin synthesis were identified, including two CHS, two CHI, one F3H, one F3 5 H, one F3 H, and two DFR, two ANS, two OMT, and one CYP73A.A positive correlation was found between anthocyanin accumulation and structural gene expression in fruit trees.It was determined that the structural genes SmCHS, SmCHI, SmF3H, SmF3 5 H, SmDFR, and SmANS displayed a higher abundance in purple eggplant pericarp compared with white eggplant pericarp [49].A similar conclusion was also obtained from wild-type purple eggplants and mutants with other colors treated by EMS.This indicated that structural genes play a key role in the biosynthesis of anthocyanin in purple eggplant pericarp [50].The genes responsible for anthocyanin synthesis structure, including CHS, F3H, F3 H, DER, and ANS, consistently exhibit higher levels of expression in purple plant tissue, resulting in significantly higher anthocyanin content compared to green or white tissue [51][52][53][54].In our study, 14 anthocyanin synthesis structural genes were up-regulated in the MP group, and only one DFR was down-regulated.This result was consistent with previous studies [49][50][51][52][53][54], indicating that those genes were closely related to the formation of purple ring spots.
Anthocyanin synthesis structural genes directly participate in anthocyanin synthesis.In contrast, anthocyanin biosynthesis genes were found to be regulated by transcription factors MYBs.Three MYBs played a significant role in the anthocyanin synthesis pathway through a co-expression network analysis of DEGs and DAMs.The majority of MYBs could promote the biosynthesis of anthocyanins, and overexpression of AN1 and MYBA led to the up-regulation of CHS, F3 H, and DFR.This resulted in the accumulation of anthocyanins as reported, respectively, in purple potato and purple pepper [55,56].MYB2 and Pr (MYB transcription factors) up-regulate the expression of downstream structural genes DFR, ANS, and UFGT in purple cabbage and purple cauliflower, resulting in purple leaves [57,58].In addition, the transcription factor MYB plays a tissue-specific role in anthocyanin synthesis regulation, and different organs may be regulated by different MYBs.MdMYB1, MdMYBA, MdMYB10, and MYB110a were the key positive regulatory factors regulating anthocyanin synthesis in various tissues [59].
To maintain the metabolic equilibrium of anthocyanins in plants, certain MYBs act as negative regulators in anthocyanin biosynthesis.These MYBs are characterized by the presence of a repressor domain at the carboxyl end of the protein that hinders the transcription of genes associated with anthocyanin synthesis and negatively affects the production of anthocyanins [60][61][62].In our study, MYB5 and MYB113 were up-regulated in the purple pericarp, while MYB-CC was down-regulated.MYB113 has been reported to have a positive correlation with anthocyanin accumulation in other plants [63][64][65].The highly expressed BjMYB113 regulated purple leaves by enhancing the expression of anthocyanin biosynthesis genes [66].To date, there have been no reports of MYB regulating anthocyanin synthesis in papaya.We hypothesize that MYB5 and MYB113 serve as key positive regulatory factors that enhance anthocyanin synthesis gene expression, while MYB-CC acts as a negative regulatory factor inhibiting the expression of anthocyanin biosynthesis genes.As a result of these findings, we speculate that purple ring spots on papaya fruit tips are formed when anthocyanins accumulate in the purple pericarp.

Plant Materials
The papaya variety 'Zihui' was used as the research material.Papaya fruits were sampled in two developmental stages, 60 and 120 days after anthesis (60DAA, 120DAA).The pericarps from the fruits of six individual plants were collected, pooled, and then immediately frozen and stored in liquid nitrogen for transcriptome and metabolome analysis.The biological replicates were conducted three times for each stage.The fruit pericarp of the green and purple parts for 60DAA and 120DAA were named YG, YP, MG, and MP, respectively.The fruits of 'Zihui' were planted in the Baiyun district of Guangzhou, Guangdong, China.Field management was the same as for the local papaya materials.

RNA Sequencing and Data Analyses
The total RNA of pericarp samples was extracted with three biological replicates using a FastQuant RT Kit (Takara, Dalian, China) following the manufacturer's protocol.Nan-oDrop 2000, Qubit 2.0, Agilent 2100, etc., were used to assess RNA sample purity, concentration, and integrity.The RNA samples that met the requirements (concentration greater than 100 ng•uL −1 ) were used to construct an RNA library, as previously described [67], and then subjected to an Illumina sequencing platform (Illumina HiSeq X-ten, Illumina, Foster, CA, USA) to generate 150 bp paired-end reads.Fastq format was used to store the raw data.To obtain clean reads for further analysis, reads with poly-N, low quality, and less than 150 bp in length were filtered.HISAT2 version 2.0.5 was then used to align the clean reads to the papaya reference genome (http://www.ncbi.nlm.nih.gov/gene/?term=papaya, accessed on 12 August 2022).

Differentially Expressed Gene Identification
We calculated gene expression abundance based on fragments per kilobase of transcript per million fragments (FPKM), which were mapped according to gene length.Genes that were differentially expressed (DEGs) were identified by the DESeq R package, with a standard of |log 2 (foldchange)| ≥ 1 as well as a false discovery rate (FDR) < 0.05.The enrichment analysis (p-value ≤ 0.05) was performed based on the Gene Ontology (GO) terms and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway.

Metabolites Extraction
The samples were freeze-dried by Scientz-100F and ground into powder (30 Hz, 1.5 min).A total of 50mg of each individual sample was accurately weighed and transferred to an Eppendorf tube.Then, 500 µL of extract solution (methanol/water/hydrochloric acid (500:500:1, V/V/V) was added, and after the sample was dissolved, the mixture was vortex oscillated once every 30 min for a total of 5 times to improve the extraction rate and refrigerated overnight at 4 • C. The extraction was then performed, as previously described [68].

UPLC-MS-MS Analysis
A UPLC-ESI-MS/MS system was used to analyze the sample extracts described previously [69].Anthocyanins were analyzed using scheduled multiple reaction monitoring (MRM).Those metabolites with multiple |log 2 (fold change)| greater than 1, a p-value less than 0.05, and variable importance in the project (VIP) greater than 1 were considered differentially accumulated metabolites (DAMs).

Quantitative Reverse-Transcriptase PCR Validation
Total RNA of the pericarp samples was extracted using a FastQuant RT Kit (Takara) following the manufacturer's instructions.Ten DEGs were selected for quantitative reversetranscriptase PCR (qRT-PCR) validation.Primers were designed using NCBI primer BLAST (http://www.ncbi.nlm.nih.gov/tools/primer-blast/,accessed on 7 January 2023) and acted as the internal reference control.The qRT-PCR was performed using Bio-Rad's CFX96 Real-Time System and the reactions program, following the method in the previous study [70].The relative transcriptional levels were calculated using the 2 −∆∆Ct method.The primers are listed in Table S1.

Conclusions
In summary, we conducted transcriptome, plant hormone, and metabolome profiling analyses of green and purple pericarp, aiming to identify the potential mechanism of purple pericarp formation in papaya.In total, 650 DEGs were identified in the transcriptome sequencing between MG and MP, including 15 anthocyanin synthesis structure genes.Among these, 14 DEGs were up-regulated, except for a DFR in the purple pericarp that led to the accumulation of anthocyanins.Metabolome profiling found four plant hormones (especially abscisic acid and 1-Aminocyclopropanecarboxylic acid), and 27 anthocyanins showed significant differences in the MP and MG papaya pericarp.This result was consistent with the transcriptome analysis, which found all anthocyanins to be up-regulated.Moreover, three MYBs may regulate anthocyanin biosynthesis structure genes, crucial for anthocyanin synthesis.In addition to identifying numerous genes involved in the anthocyanin synthesis pathway of 'Zihui' papaya, this study offered new insights into the mechanism for purple pericarp formation.

Figure 1 .
Figure 1.Phenotype of 'Zihui' papaya fruit in two developmental stages, 60 and 120 days after anthesis (60DAA, 120DAA).The fruit pericarp of green and purple parts for 60DAA and 120DAA were named YG, YP, MG, and MP, respectively.(a) Indicates fruit length, (b) reflects the fruit longitudinal diameter of the fruit.

Figure 1 .
Figure 1.Phenotype of 'Zihui' papaya fruit in two developmental stages, 60 and 120 days after anthesis (60DAA, 120DAA).The fruit pericarp of green and purple parts for 60DAA and 120DAA were named YG, YP, MG, and MP, respectively.(a) Indicates fruit length, (b) reflects the fruit longitudinal diameter of the fruit.

Figure 3 .
Figure 3. Enriched KEGG pathways of target genes for DEGS between YP&YG and MP&MG.Hotter colors correspond to lower p-values.The size of the circle indicates the number of DEGs.

Figure 4 .
Figure 4. Metabolome profiling of MP and MG papaya pericarp related to anthocyanin metabolism.(a) Principal component analysis (PCA) score plots for MP and MG groups.(b) Spearman rank correlation for the metabolite determination results among different samples; the darker the color, the higher the correlation coefficient, while the lighter color indicate the lower correlation coefficient.. (c) Primary classification of the DAMs identified between MP and MG groups.(d) Up-regulated and unchanged metabolite numbers between MP and MG groups; red dots indicate up-regulated metabolites and gray dots indicate unchanged metabolites.(e) Cluster heat map of the 35 anthocyanins in the MP and MG groups.Intensity values were adjusted by log transformation and then normalized.

Figure 4 .
Figure 4. Metabolome profiling of MP and MG papaya pericarp related to anthocyanin metabolism.(a) Principal component analysis (PCA) score plots for MP and MG groups.(b) Spearman rank correlation for the metabolite determination results among different samples; the darker the color, the higher the correlation coefficient, while the lighter color indicate the lower correlation coefficient.(c) Primary classification of the DAMs identified between MP and MG groups.(d) Up-regulated and unchanged metabolite numbers between MP and MG groups; red dots indicate up-regulated metabolites and gray dots indicate unchanged metabolites.(e) Cluster heat map of the 35 anthocyanins in the MP and MG groups.Intensity values were adjusted by log transformation and then normalized.

Figure 6 .
Figure 6.Correlation networks of the DEGs and DAMs involved in the anthocyanin biosynthesis pathway.The DEGs and DAMs are on the right and left of the figure, respectively.As shown by the pink line, the correlation was highly positive, while the cyan line represented a highly negative correlation.

Figure 6 .
Figure 6.Correlation networks of the DEGs and DAMs involved in the anthocyanin biosynthesis pathway.The DEGs and DAMs are on the right and left of the figure, respectively.As shown by the pink line, the correlation was highly positive, while the cyan line represented a highly negative correlation.

Figure 7 .
Figure 7. qRT-PCR verification of DEG relative expression levels.Standard deviations of three biological replicates are shown in the error bars.An evaluation of 10 DEGs using qRT-PCR and RNAseq is shown in the gray and black bars, respectively.

Figure 7 .
Figure 7. qRT-PCR verification of DEG relative expression levels.Standard deviations of three biological replicates are shown in the error bars.An evaluation of 10 DEGs using qRT-PCR and RNA-seq is shown in the gray and black bars, respectively.

Table 1 .
Summary of RNA sequencing data.

Table 1 .
Summary of RNA sequencing data.

Table 2 .
The comparison of contents of different plant hormones between MG and MP pericarp.