The MdNAC72-MdABI5 module acts as an interface integrating jasmonic acid and gibberellin signals and undergoes ubiquitination-dependent degradation regulated by MdSINA2 in apple.

Jasmonic acid (JA) and gibberellin (GA) coordinately regulate plant developmental programs and environmental cue responses. However, the fine regulatory network of the cross-interaction between JA and GA remains largely elusive. In this study, we demonstrate that MdNAC72 together with MdABI5 positively regulates anthocyanin biosynthesis through an exquisite MdNAC72-MdABI5-MdbHLH3 transcriptional cascade in apple. MdNAC72 interacts with MdABI5 to promote the transcriptional activation of MdABI5 on its target gene MdbHLH3 and directly activates the transcription of MdABI5. The MdNAC72-MdABI5 module regulates the integration of JA and GA signals in anthocyanin biosynthesis by combining with JA repressor MdJAZ2 and GA repressor MdRGL2a. MdJAZ2 disrupts the MdNAC72-MdABI5 interaction and attenuates the transcriptional activation of MdABI5 by MdNAC72. MdRGL2a sequesters MdJAZ2 from the MdJAZ2-MdNAC72 protein complex, leading to the release of MdNAC72. The E3 ubiquitin ligase MdSINA2 is responsive to JA and GA signals and promotes ubiquitination-dependent degradation of MdNAC72. The MdNAC72-MdABI5 interface fine-regulates the integration of JA and GA signals at the transcriptional and posttranslational levels by combining MdJAZ2, MdRGL2a, and MdSINA2. In summary, our findings elucidate the fine regulatory network connecting JA and GA signals with MdNAC72-MdABI5 as the core in apple.

Improving plant miRNA-target prediction with self-supervised k-mer embedding and spectral graph convolutional neural network.

Deciphering the targets of microRNAs (miRNAs) in plants is crucial for comprehending their function and the variation in phenotype that they cause. As the highly cell-specific nature of miRNA regulation, recent computational approaches usually utilize expression data to identify the most physiologically relevant targets. Although these methods are effective, they typically require a large sample size and high-depth sequencing to detect potential miRNA-target pairs, thereby limiting their applicability in improving plant breeding. In this study, we propose a novel miRNA-target prediction framework named kmerPMTF (k-mer-based prediction framework for plant miRNA-target). Our framework effectively extracts the latent semantic embeddings of sequences by utilizing k-mer splitting and a deep self-supervised neural network. We construct multiple similarity networks based on k-mer embeddings and employ graph convolutional networks to derive deep representations of miRNAs and targets and calculate the probabilities of potential associations. We evaluated the performance of kmerPMTF on four typical plant datasets: Arabidopsis thaliana, Oryza sativa, Solanum lycopersicum, and Prunus persica. The results demonstrate its ability to achieve AUPRC values of 84.9%, 91.0%, 80.1%, and 82.1% in 5-fold cross-validation, respectively. Compared with several state-of-the-art existing methods, our framework achieves better performance on threshold-independent evaluation metrics. Overall, our study provides an efficient and simplified methodology for identifying plant miRNA-target associations, which will contribute to a deeper comprehension of miRNA regulatory mechanisms in plants.

Apple E3 ligase MdPUB23 mediates ubiquitin-dependent degradation of MdABI5 to delay ABA-triggered leaf senescence.

ABSCISIC ACID-INSENSITIVE5 (ABI5) is a core regulatory factor that mediates the ABA signaling response and leaf senescence. However, the molecular mechanism underlying the synergistic regulation of leaf senescence by ABI5 with interacting partners and the homeostasis of ABI5 in the ABA signaling response remain to be further investigated. In this study, we found that the accelerated effect of MdABI5 on leaf senescence is partly dependent on MdbHLH93, an activator of leaf senescence in apple. MdABI5 directly interacted with MdbHLH93 and improved the transcriptional activation of the senescence-associated gene MdSAG18 by MdbHLH93. MdPUB23, a U-box E3 ubiquitin ligase, physically interacted with MdABI5 and delayed ABA-triggered leaf senescence. Genetic and biochemical analyses suggest that MdPUB23 inhibited MdABI5-promoted leaf premature senescence by targeting MdABI5 for ubiquitin-dependent degradation. In conclusion, our results verify that MdABI5 accelerates leaf senescence through the MdABI5-MdbHLH93-MdSAG18 regulatory module, and MdPUB23 is responsible for the dynamic regulation of ABA-triggered leaf senescence by modulating the homeostasis of MdABI5.

Phenomics and transcriptomic profiling of fruit development in distinct apple varieties.

Apple is one of the most economically important and popular temperate fruit trees. The domestication of apple has resulted in substantial phenotypic differences, particularly between wild and cultivated varieties. However, the relationship between gene expression and phenotypic variations in apple remains poorly understood. Here, we present a comprehensive dataset featuring five distinct apple varieties, including two wild varieties and three representative cultivated varieties. The dataset comprises of both phenomics data, encompassing twelve fruit quality-related traits continuously measured over two years, and transcriptomic data obtained at different developmental stages with three biological replicates. We performed basic quality control process, gene expression normalization and differential gene expression analysis to demonstrate the utility and reliability of the dataset. Our findings indicate that gene expression strongly related with phenotypic variations in apple. This dataset serves as a valuable resource, encompassing phenomics and transcriptomic data in multiple formats, thereby facilitating further exploration of the relationships between gene expression and phenotypic traits in apple.

MdbHLH162 connects the gibberellin and jasmonic acid signals to regulate anthocyanin biosynthesis in apple.

Anthocyanins are secondary metabolites induced by environmental stimuli and developmental signals. The positive regulators of anthocyanin biosynthesis have been reported, whereas the anthocyanin repressors have been neglected. Although the signal transduction pathways of gibberellin (GA) and jasmonic acid (JA) and their regulation of anthocyanin biosynthesis have been investigated, the cross-talk between GA and JA and the antagonistic mechanism of regulating anthocyanin biosynthesis remain to be investigated. In this study, we identified the anthocyanin repressor MdbHLH162 in apple and revealed its molecular mechanism of regulating anthocyanin biosynthesis by integrating the GA and JA signals. MdbHLH162 exerted passive repression by interacting with MdbHLH3 and MdbHLH33, which are two recognized positive regulators of anthocyanin biosynthesis. MdbHLH162 negatively regulated anthocyanin biosynthesis by disrupting the formation of the anthocyanin-activated MdMYB1-MdbHLH3/33 complexes and weakening transcriptional activation of the anthocyanin biosynthetic genes MdDFR and MdUF3GT by MdbHLH3 and MdbHLH33. The GA repressor MdRGL2a antagonized MdbHLH162-mediated inhibition of anthocyanins by sequestering MdbHLH162 from the MdbHLH162-MdbHLH3/33 complex. The JA repressors MdJAZ1 and MdJAZ2 interfered with the antagonistic regulation of MdbHLH162 by MdRGL2a by titrating the formation of the MdRGL2a-MdbHLH162 complex. Our findings reveal that MdbHLH162 integrates the GA and JA signals to negatively regulate anthocyanin biosynthesis. This study provides new information for discovering more anthocyanin biosynthesis repressors and explores the cross-talk between hormone signals.

Two adjacent NAC transcription factors regulate fruit maturity date and flavor in peach.

Although maturity date (MD) is an essential factor affecting fresh fruit marketing and has a pleiotropic effect on fruit taste qualities, the underlying mechanisms remain largely unclear. In this study, we functionally characterized two adjacent NAM-ATAF1/2-CUC2 (NAC) transcription factors (TFs), PpNAC1 and PpNAC5, both of which were associated with fruit MD in peach. PpNAC1 and PpNAC5 were found capable of activating transcription of genes associated with cell elongation, cell wall degradation and ethylene biosynthesis, suggesting their regulatory roles in fruit enlargement and ripening. Furthermore, PpNAC1 and PpNAC5 had pleiotropic effects on fruit taste due to their ability to activate transcription of genes for sugar accumulation and organic acid degradation. Interestingly, both PpNAC1 and PpNAC5 orthologues were found in fruit-producing angiosperms and adjacently arranged in all 91 tested dicots but absent in fruitless gymnosperms, suggesting their important roles in fruit development. Our results provide insight into the regulatory roles of NAC TFs in MD and fruit taste.

E3 ubiquitin ligases SINA4 and SINA11 regulate anthocyanin biosynthesis by targeting the IAA29-ARF5-1-ERF3 module in apple.

Auxin/indole-3-acetic acid (AUX/IAA) and auxin response factor (ARF) proteins are important components of the auxin signalling pathway, but their ubiquitination modification and the mechanism of auxin-mediated anthocyanin biosynthesis remain elusive. Here, the ARF MdARF5-1 was identified as a negative regulator of anthocyanin biosynthesis in apple, and it integrates auxin and ethylene signals by inhibiting the expression of the ethylene response factor MdERF3. The auxin repressor MdIAA29 decreased the inhibitory effect of MdARF5-1 on anthocyanin biosynthesis by attenuating the transcriptional inhibition of MdERF3 by MdARF5-1. In addition, the E3 ubiquitin ligases MdSINA4 and MdSINA11 played negative and positive regulatory roles in anthocyanin biosynthesis by targeting MdIAA29 and MdARF5-1 for ubiquitination degradation, respectively. MdSINA4 destabilized MdSINA11 to regulate anthocyanin accumulation in response to auxin signalling. In sum, our data revealed the crosstalk between auxin and ethylene signals mediated by the IAA29-ARF5-1-ERF3 module and provide new insights into the ubiquitination modification of the auxin signalling pathway.

PpARF6 acts as an integrator of auxin and ethylene signaling to promote fruit ripening in peach.

Although auxin is known to induce ethylene biosynthesis in some Rosaceae fruit crops, the mechanisms underlying the auxin-ethylene interaction during fruit ripening remain largely unknown. Here, the regulatory role of an auxin response factor, PpARF6, in fruit ripening was investigated in peach. Peach fruits showed accelerated ripening after treatment with auxin and PpARF6 was found to be significantly induced. PpARF6 not only could induce ethylene synthesis by directly activating the transcription of ethylene biosynthetic genes, but also competed with EIN3-binding F-box proteins PpEBF1/2 for binding to ethylene-insensitive3-like proteins PpEIL2/3, thereby keeping PpEIL2/3 active. Moreover, PpARF6 showed an interaction with PpEIL2/3 to enhance the PpEIL2/3-activated transcription of ethylene biosynthetic genes. Additionally, ectopic overexpression of PpARF6 in tomato accelerated fruit ripening by promoting the expression of genes involved in ethylene synthesis and fruit texture. In summary, our results revealed a positive regulatory role of PpARF6 in peach fruit ripening via integrating auxin and ethylene signaling.

Brassinosteroid signaling regulator BIM1 integrates brassinolide and jasmonic acid signaling during cold tolerance in apple.

Although brassinolide (BR) and jasmonic acid (JA) play essential roles in the regulation of cold stress responses, the molecular basis of their crosstalk remains elusive. Here, we show a key component of BR signaling in apple (Malus × domestica), BR INSENSITIVE1 (BRI1)-EMS-SUPPRESSOR1 (BES1)-INTERACTING MYC-LIKE PROTEIN1 (MdBIM1), increases cold tolerance by directly activating expression of C-REPEAT BINDING FACTOR1 (MdCBF1) and forming a complex with C-REPEAT BINDING FACTOR2 (MdCBF2) to enhance MdCBF2-activated transcription of cold-responsive genes. Two repressors of JA signaling, JAZMONATE ZIM-DOMAIN1 (MdJAZ1) and JAZMONATE ZIM-DOMAIN2 (MdJAZ2), interact with MdBIM1 to integrate BR and JA signaling under cold stress. MdJAZ1 and MdJAZ2 reduce MdBIM1-promoted cold stress tolerance by attenuating transcriptional activation of MdCBF1 expression by MdBIM1 and interfering with the formation of the MdBIM1-MdCBF2 complex. Furthermore, the E3 ubiquitin ligase ARABIDOPSIS TÓXICOS en LEVADURA73 (MdATL73) decreases MdBIM1-promoted cold tolerance by targeting MdBIM1 for ubiquitination and degradation. Our results not only reveal crosstalk between BR and JA signaling mediated by a JAZ-BIM1-CBF module but also provide insights into the posttranslational regulatory mechanism of BR signaling.

The E3 ubiquitin ligase SINA1 and the protein kinase BIN2 cooperatively regulate PHR1 in apple anthocyanin biosynthesis.

PHR1 (PHOSPHATE STARVATION RESPONSE1) plays key roles in the inorganic phosphate (Pi) starvation response and in Pi deficiency-induced anthocyanin biosynthesis in plants. However, the post-translational regulation of PHR1 is unclear, and the molecular basis of PHR1-mediated anthocyanin biosynthesis remains elusive. In this study, we determined that MdPHR1 was essential for Pi deficiency-induced anthocyanin accumulation in apple (Malus × domestica). MdPHR1 interacted with MdWRKY75, a positive regulator of anthocyanin biosynthesis, to enhance the MdWRKY75-activated transcription of MdMYB1, leading to anthocyanin accumulation. In addition, the E3 ubiquitin ligase SEVEN IN ABSENTIA1 (MdSINA1) negatively regulated MdPHR1-promoted anthocyanin biosynthesis via the ubiquitination-mediated degradation of MdPHR1. Moreover, the protein kinase apple BRASSINOSTEROID INSENSITIVE2 (MdBIN2) phosphorylated MdPHR1 and positively regulated MdPHR1-mediated anthocyanin accumulation by attenuating the MdSINA1-mediated ubiquitination degradation of MdPHR1. Taken together, these findings not only demonstrate the regulatory role of MdPHR1 in Pi starvation induced anthocyanin accumulation, but also provide an insight into the post-translational regulation of PHR1.

A comprehensive analysis of copy number variations in diverse apple populations.

BACKGROUND: As an important source of genetic variation, copy number variation (CNV) can alter the dosage of DNA segments, which in turn may affect gene expression level and phenotype. However, our knowledge of CNV in apple is still limited. Here, we obtained high-confidence CNVs and investigated their functional impact based on genome resequencing data of two apple populations, cultivars and wild relatives. RESULTS: In this study, we identified 914,610 CNVs comprising 14,839 CNV regions (CNVRs) from 346 apple accessions, including 289 cultivars and 57 wild relatives. CNVRs summed to 71.19 Mb, accounting for 10.03% of the apple genome. Under the low linkage disequilibrium (LD) with nearby SNPs, they could also accurately reflect the population structure of apple independent of SNPs. Furthermore, A total of 3,621 genes were covered by CNVRs and functionally involved in biological processes such as defense response, reproduction and metabolic processes. In addition, the population differentiation index ([Formula: see text]) analysis between cultivars and wild relatives revealed 127 CN-differentiated genes, which may contribute to trait differences in these two populations. CONCLUSIONS: This study was based on identification of CNVs from 346 diverse apple accessions, which to our knowledge was the largest dataset for CNV analysis in apple. Our work presented the first comprehensive CNV map and provided valuable resources for understanding genomic variations in apple.

The E3 ubiquitin ligases SINA1 and SINA2 integrate with the protein kinase CIPK20 to regulate the stability of RGL2a, a positive regulator of anthocyanin biosynthesis.

Although DELLA protein destabilization mediated by post-translational modifications is essential for gibberellin (GA) signal transduction and GA-regulated anthocyanin biosynthesis, the related mechanisms remain largely unknown. In this study, we report the ubiquitination and phosphorylation of an apple DELLA protein MdRGL2a in response to GA signaling and its regulatory role in anthocyanin biosynthesis. MdRGL2a could interact with MdWRKY75 to enhance the MdWRKY75-activated transcription of anthocyanin activator MdMYB1 and interfere with the interaction between anthocyanin repressor MdMYB308 and MdbHLH3 or MdbHLH33, thereby promoting anthocyanin accumulation. A protein kinase MdCIPK20 was found to phosphorylate and protect MdRGL2a from degradation, and it was essential for MdRGL2a-promoting anthocyanin accumulation. However, MdRGL2a and MdCIPK20 were ubiquitinated and degraded by E3 ubiquitin ligases MdSINA1 and MdSINA2, respectively, both of which were activated in the presence of GA. Our results display the integration of SINA1/2 with CIPK20 to dynamically regulate GA signaling and will be helpful toward understanding the mechanism of GA signal transduction and GA-inhibited anthocyanin biosynthesis. The discovery of extensive interactions between DELLA and SINA and CIPK proteins in apple will provide reference for the study of ubiquitination and phosphorylation of DELLA proteins in other species.

Dexmedetomidine enables copper homeostasis in cerebral ischemia/reperfusion via ferredoxin 1.

Excessive oxygen free radicals and toxic substances are generated in cerebral ischemia-reperfusion (I/R) process. Dexmedetomidine (DEX), a common anesthetic and sedative drug, can considerably boost glutathione (GSH), which has anti-copper influx effects. Focusing on cuproptosis, the mechanism of DEX in the I/R was revealed. Using the I/R rat model, the effects of DEX and the copper chelator D-penicillamine on cerebral infarct volume, copper levels, mitochondrial respiration and membrane potential, GSH content, and enrichment of cuproptosis functional proteins were examined. The involvement of ferredoxin 1 (FDX1) in the DEX regulatory pathway was verified by overexpressing FDX1 in vitro. DEX could significantly reduce cerebral infarction in rats, reduce copper levels, maintain mitochondrial functions, increase GSH, and reduce the content of key proteins related to cuproptosis. These aspects were replicated in vitro and revealed that FDX1 overexpression partially reversed the impacts of DEX. Together, cuproptosis occurs in the brain I/R process and DEX can enhance cell survival by blocking the primary pathway mediated by FDX1.KEY MESSAGESDexmedetomidine reduces cerebral infarction in the I/R rat models.Dexmedetomidine reduces cuproptosis in the I/R rat models.FDX1, an upstream of protein fatty acylation, mediates regulation of Dexmedetomidine.

A metabolic perspective of selection for fruit quality related to apple domestication and improvement.

BACKGROUND: Apple is an economically important fruit crop. Changes in metabolism accompanying human-guided evolution can be revealed using a multiomics approach. We perform genome-wide metabolic analysis of apple fruits collected from 292 wild and cultivated accessions representing various consumption types. RESULTS: We find decreased amounts of certain metabolites, including tannins, organic acids, phenolic acids, and flavonoids as the wild accessions transition to cultivated apples, while lysolipids increase in the "Golden Delicious" to "Ralls Janet" pedigree, suggesting better storage. We identify a total of 222,877 significant single-nucleotide polymorphisms that are associated with 2205 apple metabolites. Investigation of a region from 2.84 to 5.01 Mb on chromosome 16 containing co-mapping regions for tannins, organic acids, phenolic acids, and flavonoids indicates the importance of these metabolites for fruit quality and nutrition during breeding. The tannin and acidity-related genes Myb9-like and PH4 are mapped closely to fruit weight locus fw1 from 3.41 to 3.76 Mb on chromosome 15, a region under selection during domestication. Lysophosphatidylethanolamine (LPE) 18:1, which is suppressed by fatty acid desaturase-2 (FAD2), is positively correlated to fruit firmness. We find the fruit weight is negatively correlated with salicylic acid and abscisic acid levels. Further functional assays demonstrate regulation of these hormone levels by NAC-like activated by Apetala3/Pistillata (NAP) and ATP binding cassette G25 (ABCG25), respectively. CONCLUSIONS: This study provides a metabolic perspective for selection on fruit quality during domestication and improvement, which is a valuable resource for investigating mechanisms controlling apple metabolite content and quality.

PpHY5 is involved in anthocyanin coloration in the peach flesh surrounding the stone.

Red coloration around the stone (Cs) is an important trait of canned peaches (Prunus persica). In this study, an elongated hypocotyl 5 gene in peach termed PpHY5 was identified to participate in the regulation of the Cs trait. The E3 ubiquitin ligase PpCOP1 was expressed in the flesh around the stone and could interact with PpHY5. Although HY5 is known to be degraded by COP1 in darkness, the PpHY5 gene was activated in the flesh tissue surrounding the stone at the ripening stages and its expression was consistent with anthocyanin accumulation. PpHY5 was able to promote the transcription of PpMYB10.1 through interacting with its partner PpBBX10. Silencing of PpHY5 in the flesh around the stone caused a reduction in anthocyanin pigmentation, while transient overexpression of PpHY5 and PpBBX10 resulted in anthocyanin accumulation in peach fruits. Moreover, transgenic Arabidopsis seedlings overexpressing PpHY5 showed increased anthocyanin accumulation in leaves. Our results improve our understanding of the mechanisms of anthocyanin coloration in plants.

Colorful hues: insight into the mechanisms of anthocyanin pigmentation in fruit.

Anthocyanin is a vital indicator for both fruit nutritional and commercial value. Anthocyanin accumulation is a surprisingly complicated process mediated by multiple networks associated with genetic, developmental, hormonal, and environmental factors. Transcriptional regulation along with epigenetic regulation constitutes the dominant molecular framework for anthocyanin biosynthesis. Here, we focus on current knowledge on regulatory mechanisms of anthocyanin accumulation, with emphasis on the latest progress in transcriptional and epigenetic regulation and the crosstalk between various signaling pathways. We present an emerging picture of how various internal and external stimuli control anthocyanin biosynthesis. Additionally, we discuss the synergistic or antagonistic effect of developmental, hormonal and environmental cues on anthocyanin accumulation in fruit.

Allelic variation of MdMYB123 controls malic acid content by regulating MdMa1 and MdMa11 expression in apple.

Acidity is a key determinant of fruit organoleptic quality. Here, a candidate gene for fruit acidity, designated MdMYB123, was identified from a comparative transcriptome study of two Ma1Ma1 apple (Malus domestica) varieties, "Qinguan (QG)" and "Honeycrisp (HC)" with different malic acid content. Sequence analysis identified an A→T SNP, which was located in the last exon, resulting in a truncating mutation, designated mdmyb123. This SNP was significantly associated with fruit malic acid content, accounting for 9.5% of the observed phenotypic variation in apple germplasm. Differential MdMYB123- and mdmyb123-mediated regulation of malic acid accumulation was observed in transgenic apple calli, fruits, and plantlets. Two genes, MdMa1 and MdMa11, were up- and down-regulated in transgenic apple plantlets overexpressing MdMYB123 and mdmyb123, respectively. MdMYB123 could directly bind to the promoter of MdMa1 and MdMa11, and induce their expression. In contrast, mdmyb123 could directly bind to the promoters of MdMa1 and MdMa11, but with no transcriptional activation of both genes. In addition, gene expression analysis in 20 different apple genotypes based on SNP locus from "QG" × "HC" hybrid population confirmed a correlation between A/T SNP with expression levels of MdMa1 and MdMa11. Our finding provides valuable functional validation of MdMYB123 and its role in the transcriptional regulation of both MdMa1 and MdMa11, and apple fruit malic acid accumulation.

Genome-Wide Analysis of ATP Binding Cassette (ABC) Transporters in Peach (Prunus persica) and Identification of a Gene PpABCC1 Involved in Anthocyanin Accumulation.

The ATP-binding cassette (ABC) transporter family is a large and diverse protein superfamily that plays various roles in plant growth and development. Although the ABC transporters are known to aid in the transport of a wide range of substrates across biological membranes, their role in anthocyanin transport remains elusive. In this study, we identified a total of 132 putative ABC genes in the peach genome, and they were phylogenetically classified into eight subfamilies. Variations in spatial and temporal gene expression levels resulted in differential expression patterns of PpABC family members in various tissues of peach. PpABCC1 was identified as the most likely candidate gene essential for anthocyanin accumulation in peach. Transient overexpression of PpABCC1 caused a significant increase in anthocyanin accumulation in tobacco leaves and peach fruit, whereas virus-induced gene silencing of PpABCC1 in the blood-fleshed peach resulted in a significant decrease in anthocyanin accumulation. The PpABCC1 promoter contained an MYB binding cis-element, and it could be activated by anthocyanin-activator PpMYB10.1 based on yeast one-hybrid and dual luciferase assays. Thus, it seems that PpABCC1 plays a crucial role in anthocyanin accumulation in peach. Our results provide a new insight into the vacuolar transport of anthocyanins in peach.

Comparative network analysis reveals the dynamics of organic acid diversity during fruit ripening in peach (Prunus persica L. Batsch).

BACKGROUND: Organic acids are important components that determine the fruit flavor of peach (Prunus persica L. Batsch). However, the dynamics of organic acid diversity during fruit ripening and the key genes that modulate the organic acids metabolism remain largely unknown in this kind of fruit tree which yield ranks sixth in the world. RESULTS: In this study, we used 3D transcriptome data containing three dimensions of information, namely time, phenotype and gene expression, from 5 different varieties of peach to construct gene co-expression networks throughout fruit ripening of peach. With the network inferred, the time-ordered network comparative analysis was performed to select high-acid specific gene co-expression network and then clarify the regulatory factors controlling organic acid accumulation. As a result, network modules related to organic acid synthesis and metabolism under high-acid and low-acid comparison conditions were identified for our following research. In addition, we obtained 20 candidate genes as regulatory factors related to organic acid metabolism in peach. CONCLUSIONS: The study provides new insights into the dynamics of organic acid accumulation during fruit ripening, complements the results of classical co-expression network analysis and establishes a foundation for key genes discovery from time-series multiple species transcriptome data.