ORIGINAL RESEARCH article

Front. Plant Sci., 18 June 2025

Sec. Plant Metabolism and Chemodiversity

Volume 16 - 2025 | https://doi.org/10.3389/fpls.2025.1615281

Metabolomic and transcriptomic analyses provide insights into the red pigmentation in loquat (Eriobotrya japonica) peel

Yaling ZhangYaling ZhangXiuping ChenXiuping ChenWenbing SuWenbing SuChaojun DengChaojun DengJimou Jiang*Jimou Jiang*Shaoquan Zheng*Shaoquan Zheng*
  • Fruit Research Institute, Fujian Academy of Agricultural Sciences, Fujian, Fuzhou, China

Loquat (Eriobotrya japonica Lindl.) is a subtropical evergreen tree native to China. Generally, the pigments accumulated in the fruits of cultivated loquats are carotenoids rather than anthocyanins. In this study, we showed that the peel of ‘Shanpaisanhao’ loquats can accumulate anthocyanins and turn red. We further investigated the mechanisms underlying anthocyanin accumulation in the red-pigmented peels of ‘Shanpaisanhao’ loquats. RNA-seq analysis demonstrated that anthocyanin accumulation in loquat peel is associated with the upregulation of anthocyanin biosynthetic and transport genes as well as the transcriptional factor EjMYB10. Transient overexpression and dual luciferase assays showed that EjMYB10 could induce weak anthocyanin accumulation in tobacco leaves when co-expressed with PsbHLH3, and activate the promoters of EjF3′H, EjANS and EjUFGT. These results provide information for further elucidating the molecular mechanism of anthocyanin accumulation in the peel of SP3H loquat and for breeding of new red-pigmented loquat cultivars.

1 Introduction

Loquat (Eriobotrya japonica Lindl.) is a subtropical fruit tree native to China and is cultivated commercially worldwide nowadays. Generally, fruits of cultivated loquats only accumulate carotenoids, the pigments that endow mature fruits with yellow-, orange-, or orange-red color (Zhou et al., 2007), and do not accumulate anthocyanins. Interestingly, the fruit peels of some Eriobotrya species, such as E. henryi and E. seguinii, are red or purple, which does not seem to be a result of carotenoid accumulation. Our previous study showed that the peel of E. henryi fruits appears red color due to the accumulation of anthocyanins (Su et al., 2023). It will be interesting to examine whether other Eriobotrya species accumulate anthocyanins in their fruits.

Anthocyanins are natural pigments that are responsible for the red coloration of fruits. Anthocyanins are produced by a branch of the flavonoid biosynthetic pathway that has been well documented in plants. The anthocyanin biosynthetic enzymes include phenylalanine ammonia-lyase (PAL), cinnamate-4-hydroxylase (C4H), 4-coumaroyl:CoA-ligase (4CL), chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H), flavonoid 3′-hydroxylase(F3′H), flavonoid 3′,5′-hydroxylase, dihydroflavonol 4-reductase (DFR), anthocyanidin synthase/leucoanthocyanidin dioxygenase (ANS/LDOX) and UDP-glucose:flavonoid3-O-glucosyltransferase (UFGT) (Jaakola, 2013; Tanaka et al., 2008). Finally, anthocyanins are transported to and deposited in the vacuole with the assistance of transporters such as glutathione S-transferase (GST) (Zhao and Dixon, 2010). Su et al. (2023) demonstrated that PAL, 4CL, CHS, CHI, F3H, F3′H, DFR, ANS and UFGT were highly expressed in the peel of E. henryi fruits that accumulate anthocyanins. These results implicate that anthocyanin accumulation in loquat fruits was transcriptionally regulated.

Numerous evidence demonstrates that MYB-bHLH-WD40 (MBW) complex is the key component of the transcriptional regulatory network of anthocyanin accumulation (Allan et al., 2008; Zhao et al., 2023) and many transcriptional factors were shown to modulate anthocyanin accumulation through interactions with the MBW complex (Yan et al., 2021; Yang et al., 2022). In rosaceous species, R2R3 MYB activators have been identified as key anthocyanin regulators (Espley et al., 2007; Fang et al., 2021b; Feng et al., 2015; Gu et al., 2015; Lin-Wang et al., 2010; Ravaglia et al., 2013). However, transcription factors that regulate anthocyanin accumulation in loquat fruits remain to be identified. Lin-Wang et al. (2010) isolated the genomic DNA of MYB10 from rosaceous crop species, including loquat (Eriobotrya japonica) and showed that these sequences were highly conserved. Whether EjMYB10 is expressed and functions in anthocyanin biosynthesis remains to be tested.

We found that ‘Shanpaisanhao’ (SP3H, E. japonica Lindl.) loquat fruits turn red under certain conditions. However, it is unclear whether anthocyanins are the main pigments in the red pigmented peels and the underlying mechanism remains to be elucidated. In this study, we showed that anthocyanins were the red pigments accumulated in the peel of SP3H loquat fruits. RNA-seq results showed that the transcription of anthocyanin pathway genes and transcription factor EjMYB10 were upregulated in red pigmented peels. Transient expression of EjMYB10 induced anthocyanin production in tobacco leaves. Furthermore, dual luciferase assays showed that EjMYB10 activated the promoters of EjF3′H, EjANS and EjUFGT. This study suggests that EjMYB10 acts as an anthocyanin activator and is responsible for anthocyanin accumulation in the peel of SP3H loquat.

2 Materials and methods

2.1 Plant materials

Fruits of the SP3H loquat with obvious red pigmentation were collected from the National Loquat Germplasm Bank (Fuzhou, Fujian, China) and used for further analysis. Unpigmented and red-pigmented peels were collected and sliced into small pieces. Three biological replicates were prepared, with ten fruits were included in each replicate. The peels of each replicate were pooled and rapidly frozen in liquid nitrogen and stored at -80°C until use.

2.2 Anthocyanin extraction and targeted metabolome analysis

Peel samples were subjected to anthocyanin-targeted metabolome analysis to identify anthocyanins accumulated in red-pigmented peels. All procedures of the metabolome analysis were performed by Metware Biotechnology Co., Ltd. (Wuhan, China). Loquat peels were freeze-dried using a vacuum freeze-dryer and then ground into fine powder using a mixer mill. Anthocyanins were extracted and subjected to UPLC-MS/MS analysis as described by Fang et al. (2025). Briefly, 50 mg of powder was dissolved in 0.5 mL of extraction solution (methanol: water: hydrochloric acid=500: 500: 1, V/V/V) and the mixture was vortexed for 5 min. Then the mixture was sonicated for 5 min and centrifuged at 12, 000 g at 4°C for 3 min. The supernatant was collected and the residue was re-extracted as described above. The supernatants were pooled and filtered through a 0.22 μm pore size filter prior to LC-MS/MS analysis. UPLC-MS/MS analysis was conducted using an UPLC-ESI-MS/MS system (UPLC, ExionLC™ AD; MS, Applied Biosystems 6500 Triple Quadrupole). The parameters were as follows: the column employed was a Waters ACQUITY BEH C18 (1.7 μm, 2.1 mm×100 mm); The mobile phase was a mixture of water (0.1% formic acid): methanol(0.1% formic acid); The gradient elution was 95:5 V/V for 6 min, 50:50 V/V for 6 min, 5:95 V/V for 2 min, and then returned to 95:5 V/V for 2 min; The flow rate was 0.35 mL/min, and the temperature was maintained at 40°C; The volume of injected sample was 2 μL.

2.3 Transcriptome sequencing and analysis

Extraction of total RNA from loquat peels and library construction were carried out at Beijing BioMarker Technologies (Beijing, China) as described in our previous study (Su et al., 2023). cDNA libraries were sequenced using Illumina NovaSeq 6000. RNA-seq was performed in three biological replicates. The RNA-seq reads have been deposited in the Genome Sequence Archive (GSA) (Chen et al., 2021) of the National Genomics Data Center (Members and Partners, 2024) (https://ngdc.cncb.ac.cn/gsa) and are accessible under PRJCA038065. The obtained raw 150 bp paired-end reads were processed as described previously (Su et al., 2023) and mapped to the genome of ‘Jiefangzhong’ (Su et al., 2021) using HISAT2. DESeq2 was used to calculate gene expression level and identify differentially expressed genes (DEGs). Only genes that met with the criteria of fold change >2 and a false discovery rate (FDR) <0.01 were selected as DEGs. GO and KEGG pathway enrichment analyses was performed as described previously (Su et al., 2023). Prediction of transcription factors was performed using the Transcription Factor Prediction module of PlantTFDB v5.0 (https://planttfdb.gao-lab.org/prediction.php).

2.4 qRT-PCR (quantitative real-time PCR) analysis

Total RNA was extracted from loquat peels using the EZNA Plant RNA Kit (Omega Bio-tek). qRT-PCR was performed as described previously (Zhang et al., 2021). EjActin (Ej00095133) was used as the reference gene. A LightCycler 480 real-time PCR system was employed for qRT-PCR analysis, with three biological and four technical replicates. The primer sequences used in qRT-PCR are listed in Supplementary Table S1.

2.5 Vector construction

cDNA from red pigmented SP3H loquat peels was synthesized using the HiScript III 1st Strand cDNA Synthesis Kit (Vazyme, Nanjing, China). The isolation of EjMYB10 coding sequence and insertion into pSAK277 were carried out using the 2×Phanta Max Master Mix (Vazyme, Nanjing, China) and the ClonExpress Ultra One Step Cloning Kit (Vazyme, Nanjing, China), respectively. The promoter of EjF3′H, EjANS and EjUFGT were amplified and cloned into pGreenII 0800-LUC to generate reporter construct as described previously (Zhang et al., 2021) and their sequences are provided in Supplementary Table S2. All primers used for cloning are listed in Supplementary Table S3.

2.6 Phylogenetic analysis and sequence alignment

Phylogenetic analyses were performed using the neighbor-joining method with 1000 bootstrap replicates by MEGA6. The alignment of amino acid sequences was carried out using ClustalW (https://www.genome.jp/toolsbin/clustalw). ESPript 3.0 (Robert and Gouet, 2014) was employed to shade the results of the multiple sequence alignment. The accession numbers of additional sequences from other species are provided in Supplementary Table S4.

2.7 Transient expression in tobacco leaf and determination of total anthocyanin content

Transient expression of EjMYB10 was carried out in Nicotiana tabacum young leaves as described previously (Zhang et al., 2021). Basic helix-loop-helix (bHLH) transcription factors that belong to the bHLH2 subgroup (PhAN1/AtTT8) have been suggested as indispensable partners of anthocyanin-promoting R2R3-MYBs (Espley et al., 2007; Gonzalez et al., 2008). These bHLHs are required for anthocyanin biosynthesis in tobacco leaves (Montefiori et al., 2015). Sequence analysis showed that Ej00089823 (EjbHLH3) was predicted to encode a bHLH2 subgroup bHLH protein. EjbHLH3 was highly homologous to apple MdbHLH3 (Xie et al., 2012), peach PpbHLH3 (Ravaglia et al., 2013) and plum PsbHLH3 (Fang et al., 2021a) (Supplementary Figure S1). We failed to clone EjbHLH3, so we chose the plum PsbHLH3 which was sufficient to induce anthocyanin accumulation in tobacco leaves when coinfiltrated with plum MYB10s (Fang et al., 2021a, Fang et al., 2021b) and blueberry VcMYBA (Zhang et al., 2021). Agrobacteria carrying constructs were cultivated and resuspended in the infiltration buffer (10 mM MES, 10 mM MgCl2, and 100 μM acetosyringone) to an OD600 of 0.6 and then infiltrated into tobacco leaves. Photos of transformed tobacco leaves were taken 7 d after infiltration. Anthocyanin in tobacco leaves was extracted and quantified as described by Fang et al. (2025).

2.8 Dual-luciferase assay

pGreenII 0800-LUC vectors carrying the promoter of EjF3′H, EjANS and EjUFGT were transformed into the A. tumefacien strain GV3101 with pSoup. Dual-luciferase assays were performed using 3 to 4-week-old Nicotiana benthamiana leaves. Agrobacteria were grown and resuspended in the infiltration buffer to an OD600 of 0.5 and then used for infiltration according to the protocol described above for transient expression assays. Firefly luciferase (LUC) activity detection and image capture were performed 3d after infiltration, as described by Fang et al. (2025). The detection of firefly luciferase and Renilla luciferase activities was carried out using a GloMax Discover instrument (Promega) and the Dual Luciferase Reporter Gene Assay Kit (Yeasen, 11402ES80, China).

2.9 Statistical analysis

Statistical analysis was performed using Student’s t-test (*P < 0.05, **P < 0.01, ***P<0.001 and ****P<0.0001). Data analysis was performed using GraphPad Prism 10.3.1 and TBtools-II (Chen et al., 2023).

3 Results

3.1 Anthocyanin accumulation in the peel of SP3H loquat fruits

Generally, the ripe SP3H loquat fruits are yellow. However, partial red pigmentation was found in the peel of some ‘Shanpaisanhao’ loquat fruits (Figures 1A–C). Our previous study showed that anthocyanins are the main pigments accumulated in red loquat peels (Su et al., 2023). Flavonoids metabolomics analysis was performed to investigate the composition of red pigments in the red pigmented areas. In total, 41 flavonoid metabolites, including seven cyanidin, five pelargonidin, nine delphinidin, five petunidin, two peonidin and two malvidin derivatives were identified (Figure 1D). 21 of the identified anthocyanins were significantly accumulated in the red-pigmented areas, and cyanidin-3-O-galactoside was the most predominant anthocyanin (Figure 1D). In addition, the accumulation of rutin, naringenin-7-O-glucoside and quercetin-3-O-glucoside were also detected in the red pigmented areas (Figure 1C). These results suggest that anthocyanin accumulation is responsible for the red pigmentation in the peel of SP3H loquat fruits.

Figure 1
www.frontiersin.org

Figure 1. Anthocyanin accumulation in the peel of SP3H loquat. (A) fruits of SP3H loquat. (B) normal SP3H loquat fruit. (C) red pigmented SP3H loquat fruit. (D) heatmap of flavonoids identified in the peel of SP3H loquat. Each value is the average for three biological replicates. The contents of flavonoids were normalized with the maximum content values of each flavonoid. The content of flavonoids is indicated using filled circles of different sizes and colors. The larger the circle, the higher the expression. The values on the right indicate the highest content value of each flavonoid.

3.2 RNA-seq and identification of differentially expressed genes

The yellow and red pigmented peels of SP3H loquat fruits were further subjected to RNA-seq analysis to explore the mechanism of anthocyanin biosynthesis. A total of 38.17 Gb clean data was obtained from loquat peel samples (Supplementary Table S5). Over 94% of the obtained reads could be mapped to the genome of ‘Jiefangzhong’ (Supplementary Table S5). In total, 10072 transcripts were detected in analyzed peel samples. These include 2429 putative new genes that have not been predicted in the genome of ‘Jiefangzhong’ and 1624 of them was annotated by databases (Supplementary Table S6). Transcriptome comparison showed that 2821 genes were differentially expressed (1770 upregulated and 1051 downregulated) between yellow and red pigmented areas.

3.3 Identification of differentially expressed genes involved in anthocyanin accumulation

KEGG enrichment analysis of DEGs indicated that 52 DEGs were involved in phenylpropanoid biosynthesis, flavonoid biosynthesis and anthocyanin biosynthesis pathways (Supplementary Table S7). Nine of these genes, including CHS (Ej00014264 and Ej00014465), CHI (Ej00092682, Ej00071798 and newGene15317), F3′H (Ej00025308), ANS (newGene13949) and UFGT (Ej00006885 and Ej00084941) were assigned to the anthocyanin biosynthetic pathway (Figure 2). The expression of these genes was positively correlated (r>0.65) with the cyanidin-3-O-galactoside and cyanidin-3-O-glucoside concentration (Figure 2). Other anthocyanin pathway genes, including 4CL (Ej00002975, Ej00005091, Ej00046783 and newGene1146), CHS (Ej00014720, Ej00070948 and Ej00054582), F3H (Ej00026228), DFR (Ej00081751), were not identified as differentially expressed genes, but their transcription was also enhanced in the red-pigmented peels. Additionally, the expression of a glutathione S-transferase (GST) gene Ej00043900 (annotated as GSTF12) is also positively correlated (r>0.7) with the concentrations of cyanidin-3-O-galactoside and cyanidin-3-O-glucoside (Figure 2). Phylogenetic analysis showed that Ej00043900 fell into a same clade with anthocyanin-related GSTs and was closely related to anthocyanin-related GST Raint from peach (Supplementary Figure S2). These results suggested that anthocyanin biosynthesis in peel of SP3H loquats is regulated by the transcription of these structural genes.

Figure 2
www.frontiersin.org

Figure 2. Expression of structural genes involved in anthocyanin biosynthesis, as evaluated by RNA-sequencing in the peel of SP3H loquat. The left panel indicated anthocyanin pathway including gene IDs, and the expression profile of structural genes. FPKM values of all genes were normalized with maximum FPKM values of each gene. The expression level of genes was indicated using filled circle with different size and color. The larger the circle, the higher the expression. High expression was indicated in red, while low expression was indicated in green. The values on the middle indicates the highest FPKM value of each gene. Abbreviations for pathway genes as described in the Introduction. The right panel indicated the correlation between structural genes and anthocyanin content. The correlation coefficient was indicated using filled circle with different size and color. Positive correlation was indicated in red, while negative correlation was indicated in blue.

The obvious upregulation of anthocyanin biosynthetic genes CHS, CHI, F3′H and ANS and anthocyanin transportation gene GSTF12 in red-pigmented peels suggests that anthocyanin accumulation in the SP3H loquat peel was transcriptionally regulated. To identify candidate transcriptional regulators that regulate anthocyanin accumulation, the sequences of loquat genes were submitted to PlantTFDB v5.0 to predict transcription factors. In total, 2365 genes were identified to encode transcriptional factor. We further analyzed the expression patterns of the identified transcription factors, and the results showed that 208 of them were differentially expressed (with maximum FPKM values ≥ 5, Supplementary Table S8). Among these differentially expressed transcription factors, 42 were positively correlated (r > 0.9) and 47 were negatively correlated (r|< -0.9) were positively and negatively correlated with the concentrations of cyanidin-3-O-galactoside and cyanidin-3-O-glucoside or structural genes, respectively (Figure 3).

Figure 3
www.frontiersin.org

Figure 3. Expression profiles of selected differentially expressed transcription factors. Left panel showed expression profiles of expressed transcription factors. Right panel showed the correlation coefficient between the expression of transcription factors and structural genes or anthocyanin content.

3.4 Expression analysis of candidate anthocyanin-related genes by qRT-PCR

The results of RNA-Seq were validated by analyzing the expression profiles of eight candidate anthocyanin-related genes, including six anthocyanin biosynthetic genes (CHS, CHI, F3′H, ANS and UFGT), one GST gene, and one MYB, using qRT-PCR. The results showed that all analyzed genes were upregulated in the red-pigmented peel (Figure 4).

Figure 4
www.frontiersin.org

Figure 4. Expression of candidate genes involved in anthocyanin biosynthesis. Actin was used as the internal control. The error bars represent the standard errors of three biological replicates.

3.5 EjMYB10 is a positive regulator of anthocyanin biosynthesis

MYBs serve as crucial regulators of anthocyanin biosynthesis. Our results demonstrated that four MYBs (Ej00002033, Ej00027194, Ej00029887 and newGene_15405) exhibited an expression profile that was positively correlated with anthocyanin concentration and the expression pattern of structural genes of anthocyanin biosynthesis pathway (Figure 3). Phylogenetic analysis showed that Ej00002033 (EjMYB10) encodes a MYB that fell into the same clade with MYB10s from other rosaceous species and was closely related to apple MdMYB110a (Figure 5A; Supplementary Table S). Multiple sequence alignment revealed EjMYB10 was highly homologous to MdMYB110a and a conserved R2R3 domain as well as the ‘ANDV’ and SG6 motifs, which are characteristic of anthocyanin-promoting MYBs (Figure 5B). However, the C-terminal of EjMYB10 protein is shorter than that of MdMYB10, PyMYB10, PsMYB10.1, and PsMYB10.2 (Figure 5B).

Figure 5
www.frontiersin.org

Figure 5. Sequence analysis of selected EjMYBs. (A) phylogenetic analysis of EjMYBs and R2R3 MYB activators from other plant species. (B) sequence alignment of EjMYB10 and anthocyanin MYB activators from other plant species. R2, R3, and bHLH binding domain are highlighted in gray, blue and green boxes, conserved ‘ANDV’ and SG6 motif for anthocyanin-promoting MYBs are indicated under black lines.

Transient expression in tobacco leaves was employed to validate the function of EjMYB10 in anthocyanin biosynthesis. Our results showed that faint red coloration was observed in leaves co-infiltrated with EjMYB10 and PsbHLH3 eight days after transformation. However, no red pigmentation was observed when EjMYB10 or PsbHLH3 was infiltrated alone (Figure 6A). Quantification of anthocyanins in tobacco leaves revealed that only a small amount of anthocyanin was accumulated in leaves infiltrated with EjMYB10 and PsbHLH3 (Figure 6B). In contrast, no anthocyanin was detected in leaves infiltrated with empty vector, EjMYB10 or PsbHLH3 (Figure 6B). These results suggest that EjMYB10 functions as a weak anthocyanin activator. To confirm this, the anthocyanin-promoting activity of EjMYB10 was further compared with strong activator PsMYB10.1 (Fang et al., 2021a) and weak activator PsMYB10.2 (Fang et al., 2021b) from plum. Coinfiltration of PsMYB10.1 with PsbHLH3 led to obvious red pigmentation at infiltration sites three days after transformation, and strong anthocyanin accumulation was observed eight days after transformation (Figure 6C). However, the infiltration of PsMYB10.2 and PsbHLH3 or EjMYB10 and PsbHLH3 resulted in only weak red pigmentation in the leaves (Figure 6C). Our results showed that the anthocyanin content in tobacco leaves infiltrated with EjMYB10 and PsbHLH3 was comparable to that in leaves infiltrated with PsMYB10.2 and PsbHLH3. However, it was significantly lower than the content in leaves infiltrated with PsMYB10.1 and PsbHLH3 (Figure 6D). Additionally, promoter structure analysis revealed that the promoters of the anthocyanin biosynthetic genes EjF3′H (Ej00025308), EjANS (newGene 13949) and EjUFGT (Ej00006885) contain multiple potential MYB binding sites (Figure 6E). Dual-luciferase assays were carried out to further verify whether EjMYB10 could interact with the promoters of these anthocyanin biosynthetic genes. Infiltration of EjMYB10 was able to enhance the transcriptional activity of EjF3′H, EjANS and EjUFGT (Figures 6F–M). These results suggest that EjMYB10 acts as a weak anthocyanin activator.

Figure 6
www.frontiersin.org

Figure 6. Functional analysis of the loquat EjMYB10 gene. (A, C) transient color assay of the EjMYB10 activity in tobacco leaf. (B, D) anthocyanin content in tobacco leaves. (E) schematic diagram of the EjF3′H, EjANS and EjUFGT promoter. The search of potential MYB-binding elements described by Zhou et al. (2019) was carried out using the PlantPAN4.0 (https://plantpan.itps.ncku.edu.tw/plantpan4/promoter_analysis.php). Color ellipses represent predicted MYB binding sites. Red dashed line indicated the location of sequence deletion in the promoter of EjF3′H. (F) schematic diagram showing the effector and reporter constructs for the dual-luciferase assay. (G–M) Validation of activation effect of EjMYB10 on the EjF3′H, EjANS and EjUFGT promoters using dual luciferase assay. (G, J, L) Representative images of tobacco leaves 72 h after infiltration. (H, I, K, M) relative firefly LUC-to-REN ratios from transient expression assays. Asterisks denote t-test significance: *P < 0.05, **P < 0.01.

4 Discussion

Commercially cultivated loquats are unable to accumulate anthocyanins in the peel and flesh. However, the accumulation of anthocyanins in the peel and/or flesh enhances the appearance and nutritional quality of loquat fruits. Understanding the mechanism underlying anthocyanin biosynthesis in loquat fruits will provide insights into the development of new cultivars that are rich in anthocyanins. In the present study, we showed that some SP3H loquat fruits accumulate anthocyanins in their red-pigmented peels. We further investigated the mechanism responsible for anthocyanin production in the peel of SP3H loquat fruits.

The composition of anthocyanins varies among species. Cyanidins have been reported to be the most abundant anthocyanins in fruits of rosaceous species, including apple (Feng et al., 2013; Tsao et al., 2003), pear (Pierantoni et al., 2010; Zhang et al., 2020), peach (Cheng et al., 2014; Liu et al., 2019), plum (Fang et al., 2021a, Fang et al., 2025, Fang et al., 2021b; Usenik et al., 2009) and cherry (Liu et al., 2011). In apple, cyanidin-3-galactoside was the most predominant anthocyanin (Espley et al., 2013), while cyanidin-3-glucoside represents the most abundant anthocyanin in peach (Cheng et al., 2014; Jiao et al., 2014; Zhou et al., 2015) and Chinese plum (Fang et al., 2025, Fang et al., 2021b). However, anthocyanins are dominated by cyanidin-3-rutinoside in European plum (Usenik et al., 2009), which was reported to be the dominant anthocyanin accumulated in cherry fruits (Liu et al., 2011). These results suggested the dominant anthocyanin is species-dependent. Our results indicated that the most abundant anthocyanin in the SP3H loquat peel is cyanidin-3-O-galactoside, which was identified as that dominant anthocyanin in E. henryi fruits in our previous study (Su et al., 2023). The anthocyanin composition in other red-pigmented fruits of Eriobotrya species, such as E. seguinii, remain to be investigated.

In this study, we found that SP3H loquat fruits accumulate anthocyanins in the peel. Although the exact conditions required for inducing anthocyanin accumulation in SP3H loquat fruits are unclear, we speculate that this could be a result of environmental stimulus such as high light. We observed that only peels exposed to direct sunlight were able to accumulate anthocyanins. This suggests that light is indispensable for anthocyanin accumulation in loquat peels. However, it is noteworthy that not all fruits exposed to direct sunlight accumulate anthocyanins, and this phenomenon does not occur in SP3H loquat fruits every year. These observations suggest that the conditions required for inducing anthocyanin accumulation in loquat peels are complicated. Environmental and other factors, such as light, temperature, hormone, nutrients and mechanical damage have been demonstrated to induce anthocyanin accumulation via modulating the transcription of genes that participate in anthocyanin biosynthetic and transport in fruits (An et al., 2019; Espley and Jaakola, 2023; Fang et al., 2025; Wang et al., 2023; Zhao et al., 2023). Our previous study showed that most anthocyanin pathway genes were significantly upregulated in the red-pigmented peels of E. henryi fruits (Su et al., 2023). However, only nine genes encoding CHS, CHI, F3′H, ANS and UFGT, as well as a GST gene were significantly upregulated in red-pigmented peels of SP3H loquat fruits (Figure 2). The significant correlation between anthocyanin content and the expression of structural genes suggests that anthocyanin accumulation in SP3H loquat fruits results from the activation of anthocyanin biosynthetic and transport genes.

The enhanced transcription of anthocyanin pathway genes implies that anthocyanin biosynthesis in the peel of SP3H loquat fruits was transcriptionally regulated. R2R3 MYBs have been proven to be key anthocyanin activators in Rosaceae species (Fang et al., 2021b; Gu et al., 2015; Jin et al., 2016; Lin-Wang et al., 2010; Qian et al., 2014; Rahim et al., 2014). RNA-seq and qRT-PCR results indicated that EjMYB10 was highly expressed in the red-pigmented peel of SP3H loquat fruits (Figures 3, 4). Sequence analysis indicated that EjMYB10 encodes a R2R3 MYB that is highly homologous to the apple anthocyanin activator MdMYB110a (Umemura et al., 2013) (Figure 5). Overexpression of EjMYB10 and plum PsbHLH3 resulted in anthocyanin production and red coloration in tobacco leaves (Figures 6A, B). Dual-luciferase assays also showed that EjMYB10 enhanced transcriptional activity of anthocyanin biosynthetic genes EjF3′H, EjANS and EjUFGT (Figures 6F–M). These results suggested that EjMYB10 is an anthocyanin activator. It is noteworthy that the concentration of anthocyanin pigments induced by EjMYB10 and PsbHLH3 was comparable to that induced by the weak activator PsMYB10.2 (Fang et al., 2021b) but much less than that induced by the strong activator PsMYB10.1 (Fang et al., 2021a). These results suggest that EjMYB10 is a weak anthocyanin activator.

5 Conclusions

In this study, the mechanism responsible for the red pigmentation in the peel of SP3H loquat was investigated using metabolomic and transcriptomic analyses. Our results showed that anthocyanin accumulation is responsible for the red coloration of loquat peel and the predominant anthocyanin was cyanidin-3-O-galactoside. RNA-seq revealed that anthocyanin accumulation in loquat peel was associated with the high expression of anthocyanin structural genes and the transcription factor EjMYB10. Transient overexpression and dual luciferase assays showed that EjMYB10 can induce weak anthocyanin accumulation in tobacco leaves when co-expressed with PsbHLH3 and activate the promoters of EjF3′H, EjANS, and EjUFGT. In summary, the results presented here offer novel perspectives on the molecular mechanism underlying anthocyanin biosynthesis in the peel of SP3H loquat. Further studies should investigate the conditions for inducing red coloration and identify genes that regulate the transcription of EjMYB10.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author contributions

YZ: Writing – original draft, Methodology, Formal Analysis, Funding acquisition, Data curation, Investigation, Conceptualization, Visualization, Project administration. XC: Writing – review & editing, Investigation, Resources. WS: Writing – review & editing. CD: Writing – review & editing. JJ: Supervision, Conceptualization, Writing – review & editing. SZ: Writing – review & editing, Supervision, Conceptualization.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was funded by the Basic Scientific Research Funds of Public Welfare Scientific Research Institutes of Fujian Province(2022R1028008), the Natural Science Foundation of Fujian Province (2022J01474), National Science and technology resource sharing service platform project (NHGRC2024-NH18-1), the Technology Innovation Team Program from Fujian Academy of Agricultural Science (CXTD2021004-1).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2025.1615281/full#supplementary-material

References

Allan, A. C., Hellens, R. P., and Laing, W. A. (2008). MYB transcription factors that colour our fruit. Trends Plant Sci. 13, 99–102. doi: 10.1016/j.tplants.2007.11.012

PubMed Abstract | Crossref Full Text | Google Scholar

An, J. P., Zhang, X. W., You, C. X., Bi, S. Q., Wang, X. F., and Hao, Y. J. (2019). MdWRKY40 promotes wounding-induced anthocyanin biosynthesis in association with MdMYB1 and undergoes MdBT2-mediated degradation. New Phytol. 224, 380–395. doi: 10.1111/nph.v224.1

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, C., Wu, Y., Li, J., Wang, X., Zeng, Z., Xu, J., et al. (2023). TBtools-II: A “one for all, all for one” bioinformatics platform for biological big-data mining. Mol. Plant 16, 1733–1742. doi: 10.1016/j.molp.2023.09.010

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, T., Chen, X., Zhang, S., Zhu, J., Tang, B., Wang, A., et al. (2021). The Genome Sequence Archive Family: toward explosive data growth and diverse data types. Genomics Proteomics Bioinf. 19, 578–583. doi: 10.1016/j.gpb.2021.08.001

PubMed Abstract | Crossref Full Text | Google Scholar

Cheng, J., Wei, G., Zhou, H., Gu, C., Vimolmangkang, S., Liao, L., et al. (2014). Unraveling the mechanism underlying the glycosylation and methylation of anthocyanins in peach. Plant Physiol. 166, 1044–1058. doi: 10.1104/pp.114.246876

PubMed Abstract | Crossref Full Text | Google Scholar

Espley, R. V., Bovy, A., Bava, C., Jaeger, S. R., Tomes, S., Norling, C., et al. (2013). Analysis of genetically modified red-fleshed apples reveals effects on growth and consumer attributes. Plant Biotechnol. J. 11, 408–419. doi: 10.1111/pbi.2013.11.issue-4

PubMed Abstract | Crossref Full Text | Google Scholar

Espley, R. V., Hellens, R. P., Putterill, J., Stevenson, D. E., Kutty-Amma, S., and Allan, A. C. (2007). Red colouration in apple fruit is due to the activity of the MYB transcription factor, MdMYB10. Plant J. 49, 414–427. doi: 10.1111/j.1365-313X.2006.02964.x

PubMed Abstract | Crossref Full Text | Google Scholar

Espley, R. V. and Jaakola, L. (2023). The role of environmental stress in fruit pigmentation. Plant Cell Environ. 46, 3663–3679. doi: 10.1111/pce.v46.12

PubMed Abstract | Crossref Full Text | Google Scholar

Fang, Z., Lin-Wang, K., Jiang, C., Zhou, D., Lin, Y., Pan, S., et al. (2021a). Postharvest temperature and light treatments induce anthocyanin accumulation in peel of ‘Akihime’ plum (Prunus salicina Lindl.) via transcription factor PsMYB10.1. Postharvest Biol. Technol. 179, 111592. doi: 10.1016/j.postharvbio.2021.111592

Crossref Full Text | Google Scholar

Fang, Z., Lin-Wang, K., Lin, Y., and Espley, R. V. (2025). Metabolomic and transcriptomic analyses provide insights into temperature and light regulated anthocyanin accumulation in flesh of ‘Furongli’ plum (Prunus salicina Lindl.). Postharvest Biol. Technol. 221, 113326. doi: 10.1016/j.postharvbio.2024.113326

Crossref Full Text | Google Scholar

Fang, Z., Lin-Wang, K., Zhou, D., Lin, Y., Jiang, C., Pan, S., et al. (2021b). Activation of PsMYB10.2 transcription causes anthocyanin accumulation in flesh of the red-fleshed mutant of ‘Sanyueli’ (Prunus salicina Lindl.). Front. Plant Sci. 12, 1167. doi: 10.3389/fpls.2021.680469

PubMed Abstract | Crossref Full Text | Google Scholar

Feng, F., Li, M., Ma, F., and Cheng, L. (2013). Phenylpropanoid metabolites and expression of key genes involved in anthocyanin biosynthesis in the shaded peel of apple fruit in response to sun exposure. Plant Physiol. Biochem. 69, 54–61. doi: 10.1016/j.plaphy.2013.04.020

PubMed Abstract | Crossref Full Text | Google Scholar

Feng, S., Sun, S., Chen, X., Wu, S., Wang, D., and Chen, X. (2015). PyMYB10 and PyMYB10.1 interact with bHLH to enhance anthocyanin accumulation in pears. PloS One 10, e0142112. doi: 10.1371/journal.pone.0142112

PubMed Abstract | Crossref Full Text | Google Scholar

Gonzalez, A., Zhao, M., Leavitt, J. M., and Lloyd, A. M. (2008). Regulation of the anthocyanin biosynthetic pathway by the TTG1/bHLH/Myb transcriptional complex in Arabidopsis seedlings. Plant J. 53, 814–827. doi: 10.1111/j.1365-313X.2007.03373.x

PubMed Abstract | Crossref Full Text | Google Scholar

Gu, C., Liao, L., Zhou, H., Wang, L., Deng, X., and Han, Y. (2015). Constitutive activation of an anthocyanin regulatory gene PcMYB10.6 is related to red coloration in purple-foliage plum. PloS One 10, e0135159. doi: 10.1371/journal.pone.0135159

PubMed Abstract | Crossref Full Text | Google Scholar

Jaakola, L. (2013). New insights into the regulation of anthocyanin biosynthesis in fruits. Trends Plant Sci. 18, 477–483. doi: 10.1016/j.tplants.2013.06.003

PubMed Abstract | Crossref Full Text | Google Scholar

Jiao, Y., Ma, R., Shen, Z., Yan, J., and Yu, M. (2014). Gene regulation of anthocyanin biosynthesis in two blood-flesh peach (Prunus persica (L.) Batsch) cultivars during fruit development. J. Zhejiang Univ. Sci. B 15, 809–819. doi: 10.1631/jzus.B1400086

PubMed Abstract | Crossref Full Text | Google Scholar

Jin, W., Wang, H., Li, M., Wang, J., Yang, Y., Zhang, X., et al. (2016). The R2R3 MYB transcription factor PavMYB10.1 involves in anthocyanin biosynthesis and determines fruit skin colour in sweet cherry (Prunus avium L.). Plant Biotechnol. J. 14, 2120–2133. doi: 10.1111/pbi.2016.14.issue-11

PubMed Abstract | Crossref Full Text | Google Scholar

Lin-Wang, K., Bolitho, K., Grafton, K., Kortstee, A., Karunairetnam, S., McGhie, T., et al. (2010). An R2R3 MYB transcription factor associated with regulation of the anthocyanin biosynthetic pathway in Rosaceae. BMC Plant Biol. 10, 50. doi: 10.1186/1471-2229-10-50

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, H., Cao, K., Zhu, G., Fang, W., Chen, C., Wang, X., et al. (2019). Genome-wide association analysis of red flesh character based on resequencing approach in peach. J. Am. Soc Hortic. Sci. 144, 209–216. doi: 10.21273/JASHS04622-18

Crossref Full Text | Google Scholar

Liu, Y., Liu, X., Zhong, F., Tian, R., Zhang, K., Zhang, X., et al. (2011). Comparative study of phenolic compounds and antioxidant activity in different species of cherries. J. Food Sci. 76, 633–638. doi: 10.1111/j.1750-3841.2011.02150.x

PubMed Abstract | Crossref Full Text | Google Scholar

Members, C.N., Partners (2024). Database resources of the national genomics data center, China national center for bioinformation in 2025. Nucleic Acids Res. 53, D30–D44. doi: 10.1093/nar/gkae978

PubMed Abstract | Crossref Full Text | Google Scholar

Montefiori, M., Brendolise, C., Dare, A. P., Lin-Wang, K., Davies, K. M., Hellens, R. P., et al. (2015). In the Solanaceae, a hierarchy of bHLHs confer distinct target specificity to the anthocyanin regulatory complex. J. Exp. Bot. 66, 1427–1436. doi: 10.1093/jxb/eru494

PubMed Abstract | Crossref Full Text | Google Scholar

Pierantoni, L., Dondini, L., De Franceschi, P., Musacchi, S., Winkel, B. S. J., and Sansavini, S. (2010). Mapping of an anthocyanin-regulating MYB transcription factor and its expression in red and green pear, Pyrus communis. Plant Physiol. Biochem. 48, 1020–1026. doi: 10.1016/j.plaphy.2010.09.002

PubMed Abstract | Crossref Full Text | Google Scholar

Qian, M., Sun, Y., Allan, A. C., Teng, Y., and Zhang, D. (2014). The red sport of ‘Zaosu’ pear and its red-striped pigmentation pattern are associated with demethylation of the PyMYB10 promoter. Phytochemistry 107, 16–23. doi: 10.1016/j.phytochem.2014.08.001

PubMed Abstract | Crossref Full Text | Google Scholar

Rahim, M. A., Busatto, N., and Trainotti, L. (2014). Regulation of anthocyanin biosynthesis in peach fruits. Planta 240, 913–929. doi: 10.1007/s00425-014-2078-2

PubMed Abstract | Crossref Full Text | Google Scholar

Ravaglia, D., Espley, R., Henry-Kirk, R., Andreotti, C., Ziosi, V., Hellens, R., et al. (2013). Transcriptional regulation of flavonoid biosynthesis in nectarine (Prunus persica) by a set of R2R3 MYB transcription factors. BMC Plant Biol. 13, 68. doi: 10.1186/1471-2229-13-68

PubMed Abstract | Crossref Full Text | Google Scholar

Robert, X. and Gouet, P. (2014). Deciphering key features in protein structures with the new ENDscript server. Nucleic Acids Res. 42, 320–324. doi: 10.1093/nar/gku316

PubMed Abstract | Crossref Full Text | Google Scholar

Su, W., Jing, Y., Lin, S., Yue, Z., Yang, X., Xu, J., et al. (2021). Polyploidy underlies co-option and diversification of biosynthetic triterpene pathways in the apple tribe. Proc. Natl. Acad. Sci. U.S.A. 118, e2101767118. doi: 10.1073/pnas.2101767118

PubMed Abstract | Crossref Full Text | Google Scholar

Su, W., Zhu, C., Fan, Z., Huang, M., Lin, H., Chen, X., et al. (2023). Comprehensive metabolome and transcriptome analyses demonstrate divergent anthocyanin and carotenoid accumulation in fruits of wild and cultivated loquats. Front. Plant Sci. 14, 1285456. doi: 10.3389/fpls.2023.1285456

PubMed Abstract | Crossref Full Text | Google Scholar

Tanaka, Y., Sasaki, N., and Ohmiya, A. (2008). Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant J. 54, 733–749. doi: 10.1111/j.1365-313X.2008.03447.x

PubMed Abstract | Crossref Full Text | Google Scholar

Tsao, R., Yang, R., Young, J. C., and Zhu, H. (2003). Polyphenolic profiles in eight apple cultivars using high-performance liquid chromatography(HPLC). J. Agr Food Chem. 51, 6347–6353. doi: 10.1021/jf0346298

PubMed Abstract | Crossref Full Text | Google Scholar

Umemura, H., Otagaki, S., Wada, M., Kondo, S., and Matsumoto, S. (2013). Expression and functional analysis of a novel MYB gene, MdMYB110a_JP, responsible for red flesh, not skin color in apple fruit. Planta 238, 65–76. doi: 10.1007/s00425-013-1875-3

PubMed Abstract | Crossref Full Text | Google Scholar

Usenik, V., Štampar, F., and Veberič, R. (2009). Anthocyanins and fruit colour in plums (Prunus domestica L.) during ripening. Food Chem. 114, 529–534. doi: 10.1016/j.foodchem.2008.09.083

Crossref Full Text | Google Scholar

Wang, L., Yang, S., Ni, J., Teng, Y., and Bai, S. (2023). Advances of anthocyanin synthesis regulated by plant growth regulators in fruit trees. Sci. Hortic. 307, 111476. doi: 10.1016/j.scienta.2022.111476

Crossref Full Text | Google Scholar

Xie, X., Li, S., Zhang, R., Zhao, J., Chen, Y., Zhao, Q., et al. (2012). The bHLH transcription factor MdbHLH3 promotes anthocyanin accumulation and fruit colouration in response to low temperature in apples. Plant Cell Environ. 35, 1884–1897. doi: 10.1111/j.1365-3040.2012.02523.x

PubMed Abstract | Crossref Full Text | Google Scholar

Yan, H., Pei, X., Zhang, H., Li, X., Zhang, X., Zhao, M., et al. (2021). MYB-mediated regulation of anthocyanin biosynthesis. Int. J. Mol. Sci. 22, 3103. doi: 10.3390/ijms22063103

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, J., Chen, Y., Xiao, Z., Shen, H., Li, Y., and Wang, Y. (2022). Multilevel regulation of anthocyanin-promoting R2R3-MYB transcription factors in plants. Front. Plant Sci. 13, 1008829. doi: 10.3389/fpls.2022.1008829

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Y., Lin-Wang, K., Albert, N. W., Elborough, C., Espley, R. V., Andre, C. M., et al. (2021). Identification of a strong anthocyanin activator, VbMYBA, from berries of Vaccinium bracteatum Thunb. Front. Plant Sci. 12, 697212. doi: 10.3389/fpls.2021.697212

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Z., Tian, C., Zhang, Y., Li, C., Li, X., Yu, Q., et al. (2020). Transcriptomic and metabolomic analysis provides insights into anthocyanin and procyanidin accumulation in pear. BMC Plant Biol. 20, 129. doi: 10.1186/s12870-020-02344-0

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, J. and Dixon, R. A. (2010). The ‘ins’ and ‘outs’ of flavonoid transport. Trends Plant Sci. 15, 72–80. doi: 10.1016/j.tplants.2009.11.006

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, Y., Sun, J., Cherono, S., An, J., Allan, A. C., and Han, Y. (2023). Colorful hues: insight into the mechanisms of anthocyanin pigmentation in fruit. Plant Physiol. 192, 1718–1732. doi: 10.1093/plphys/kiad160

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, H., Lin-Wang, K., Wang, F., Espley, R. V., Ren, F., Zhao, J., et al. (2019). Activator-type R2R3-MYB genes induce a repressor-type R2R3-MYB gene to balance anthocyanin and proanthocyanidin accumulation. New Phytol. 221, 1919–1934. doi: 10.1111/nph.2019.221.issue-4

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, H., Lin-Wang, K., Wang, H., Gu, C., Dare, A. P., Espley, R. V., et al. (2015). Molecular genetics of blood-fleshed peach reveals activation of anthocyanin biosynthesis by NAC transcription factors. Plant J. . 82, 105–121. doi: 10.1111/tpj.2015.82.issue-1

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, C., Xu, C., Sun, C., Li, X., and Chen, K. (2007). Carotenoids in white- and red-fleshed loquat fruits. J. Agr Food Chem. 55, 7822–7830. doi: 10.1021/jf071273h

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: loquat, RNA-seq, anthocyanin biosynthesis, EjMYB10, anthocyanin-targeted metabolome

Citation: Zhang Y, Chen X, Su W, Deng C, Jiang J and Zheng S (2025) Metabolomic and transcriptomic analyses provide insights into the red pigmentation in loquat (Eriobotrya japonica) peel. Front. Plant Sci. 16:1615281. doi: 10.3389/fpls.2025.1615281

Received: 22 April 2025; Accepted: 26 May 2025;
Published: 18 June 2025.

Edited by:

Aimin Wu, South China Agricultural University, China

Reviewed by:

Qigao Guo, Southwest University, China
Wenqiu Wang, Anhui Agricultural University, China

Copyright © 2025 Zhang, Chen, Su, Deng, Jiang and Zheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jimou Jiang, amptMjUxNkAxMjYuY29t; Shaoquan Zheng, enNxMzMzNTU1QDE2My5jb20=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.