A key structural gene, AaLDOX, is involved in anthocyanin biosynthesis in all red-fleshed kiwifruit (Actinidia arguta) based on transcriptome analysis

Integrated transcriptomic and metabolomic analyses elucidate the mechanism of flavonoid biosynthesis in the regulation of mulberry seed germination under salt stress

Seed propagation is the main method of mulberry expansion in China, an important economic forest species. However, seed germination is the most sensitive stage to various abiotic stresses, especially salinity stress. To reveal the molecular regulatory mechanism of mulberry seed germination under salt stress, flavonoid metabolomics and transcriptomics analyses were performed on mulberry seeds germinated under 50 and 100 mmol/L NaCl stress. Analysis of the flavonoid metabolome revealed that a total of 145 differential flavonoid metabolites (DFMs) were classified into 9 groups, 40 flavonols, 32 flavones, 16 chalcones and 14 flavanones. Among them, 61.4% (89) of the DFMs accumulated continuously with increasing salt concentration, reaching the highest level at a 100 mmol/L salt concentration; these DFMs included quercetin-3-O-glucoside (isoquercitrin), kaempferol (3,5,7,4'-tetrahydroxyflavone), quercetin-7-O-glucoside, taxifolin (dihydroquercetin) and apigenin (4',5,7-trihydroxyflavone), indicating that these flavonoids may be key metabolites involved in the response to salt stress. Transcriptional analysis identified a total of 3055 differentially expressed genes (DEGs), most of which were enriched in flavonoid biosynthesis (ko00941), phenylpropanoid biosynthesis (ko00940) and biosynthesis of secondary metabolites (ko01110). Combined analysis of flavonoid metabolomic and transcriptomic data indicated that phenylalanine ammonia-lyase (PAL), 4-coumarate-CoA ligase (4CL), chalcone synthase (CHS), flavonol synthase (FLS), bifunctional dihydroflavonol 4-reductase/flavanone 4-reductase (DFR) and anthocyanidin reductase (ANR) were the key genes involved in flavonoid accumulation during mulberry seed germination under 50 and 100 mmol/L NaCl stress. In addition, three transcription factors, MYB, bHLH and NAC, were involved in the regulation of flavonoid accumulation under salt stress. The results of quantitative real-time PCR (qRT‒PCR) validation showed that the expression levels of 11 DEGs, including 7 genes involved in flavonoid biosynthesis, under different salt concentrations were consistent with the transcriptomic data, and parallel reaction monitoring (PRM) results showed that the expression levels of 6 key enzymes (proteins) involved in flavonoid synthesis were consistent with the accumulation of flavonoids. This study provides a new perspective for investigating the regulatory role of flavonoid biosynthesis in the regulation of mulberry seed germination under salt stress at different concentrations.

Screening and identification of photoresponse factors in kiwifruit (Actinidia arguta) development

BACKGROUND

Light is essential for kiwifruit development, in which photoresponse factors contributes greatly to the quality formation. 'Light sensitive hypocotyls, also known as light-dependent short hypocotyls' (LSH) gene family can participate in fruit development as photoresponse factor. However, the key LSH gene that determine kiwifruit development remains unclear. This study aim to screen and identify the key gene AaLSH9 in A. arguta.

MATERIALS AND METHODS

Genome-wide identification of the LSH gene family was used to analyse LSH genes in kiwifruit. Homologous cloning was used to confirm the sequence of candidate LSH genes. qRT-PCR and cluster analysis of expression pattern were used to screen the key AaLSH9 gene. Subcellular localization of AaLSH9 in tobacco leaves and overexpression of AaLSH9 in Arabidopsis thaliana hy5 mutant plants were used to define the acting place in cell and identify molecular function, respectively.

RESULTS

We identified 15 LSH genes, which were divided into two sub-families namely A and B. Domain analysis of A and B showed that they contained different domain organizations, which possibly played key roles in the evolution process. Three LSH genes, AaLSH2, AaLSH9, and AaLSH11, were successfully isolated from Actinidia arguta. The expression pattern and cluster analysis of these three AaLSH genes suggested AaLSH9 might be a key photoresponse gene participating in fruit development in A. arguta. Subcellular localization showed AaLSH9 protein was located in the nucleus. The overexpression of AaLSH9 gene in Arabidopsis thaliana hy5 mutant plants partially complemented the long hypocotyls of hy5 mutant, implying AaLSH9 played a key role as photoresponse factor in cells. In addition, the seed coat color of A. thaliana over-expressing AaLSH9 became lighter than the wide type A.thaliana. Finally, AaCOP1 was confirmed as photoresponse factor to participate in developmental process by stable transgenic A. thaliana.

CONCLUSIONS

AaLSH9 can be involved in kiwifruit (A. arguta) development as key photoresponse factor. Our results not only identified the photoresponse factors AaLSH9 and AaCOP1 but also provided insights into their key role in fruit quality improvement in the process of light response.

AcMYB10 Involved in Anthocyanin Regulation of 'Hongyang' Kiwifruit Induced via Fruit Bagging and High-Postharvest-Temperature Treatments

Light and temperature are key factors influencing the accumulation of anthocyanin in fruit crops. To assess the effects of fruit bagging during development and high post-ripening temperature on 'Hongyang' kiwifruit, we compared the pigmentation phenotypes and expression levels of anthocyanin-related genes between bagged and unbagged treatments, and between 25 °C and 37 °C postharvest storage temperatures. Both the bagging and 25 °C treatments showed better pigmentation phenotypes with higher anthocyanin concentrations. The results of the qRT-PCR analysis revealed that the gene expression levels of LDOX (leucoanthocyanidin dioxygenase), F3GT (UDP-flavonoid 3-O-glycosyltransferase ), AcMYB10, and AcbHLH42 were strongly correlated and upregulated by both the bagging treatment and 25 °C storage. The results of bimolecular fluorescence complementation and luciferase complementation imaging assays indicated an interaction between AcMYB10 and AcbHLH42 in plant cells, whereas the results of a yeast one-hybrid assay further demonstrated that AcMYB10 activated the promoters of AcLODX and AcF3GT. These results strongly suggest that enhanced anthocyanin synthesis is caused by the promoted expression of AcLODX and AcF3GT, regulated by the complex formed by AcMYB10-AcbHLH42.

Combined Widely Targeted Metabolomic, Transcriptomic, and Spatial Metabolomic Analysis Reveals the Potential Mechanism of Coloration and Fruit Quality Formation in Actinidia chinensis cv. Hongyang

Postharvest kiwifruit (Actinidia chinensis cv. Hongyang) pulp is mainly composed of outer yellow-flesh (LR) and inner red-flesh (HR). However, information about the differences in coloration and fruit quality between these two parts are limited. In this study, widely targeted metabolomic, transcriptomic, and spatial metabolomic analyses were used to reveal the potential mechanism of coloration and fruit quality formation. The results show that a total of 1001 metabolites were identified in Hongyang kiwifruit, and the accumulation of 211 metabolites were significantly higher in the HR than LR, including 69 flavonoids, 53 phenolic acids, and 38 terpenoids. There were no significant differences in the content of citric acid, quinic acid, glucose, fructose, or sucrose between the LR and HR. These results were consistent with the results from the RNA-seq profile and spatial metabolomic analysis. In addition, a total of 23 key candidate genes related to flesh color and fruit quality formation were identified and validated by qRT-PCR analysis. This study provides a theoretical basis for elucidating the underlying mechanism of the formation of kiwifruit flesh color and fruit quality.

Methylation of AcGST1 Is Associated with Anthocyanin Accumulation in the Kiwifruit Outer Pericarp

Most red-fleshed kiwifruit cultivars, such as Hongyang, only accumulate anthocyanins in the inner pericarp; the trait of full red flesh becomes the goal pursued by breeders. In this study, we identified a mutant "H-16" showing a red color in both the inner and outer pericarps, and the underlying mechanism was explored. Through transcriptome analysis, a key differentially expressed gene AcGST1 was screened out, which was positively correlated with anthocyanin accumulation in the outer pericarp. The result of McrBC-PCR and bisulfite sequencing revealed that the SG3 region (-292 to -597 bp) of AcGST1 promoter in "H-16" had a significantly lower CHH cytosine methylation level than that in Hongyang, accompanied by low expression of methyltransferase genes (MET1 and CMT2) and high expression of demethylase genes (ROS1 and DML1). Transient calli transformation confirmed that demethylase gene DML1 can activate transcription of AcGST1 to enhance its expression. Overexpression of AcGST1 enhanced the anthocyanin accumulation in the fruit flesh and leaves of the transgenic lines. These results suggested that a decrease in the methylation level of the AcGST1 promoter may contribute to accumulation of anthocyanin in the outer pericarp of "H-16".

Transcriptome and metabolome analyses of anthocyanin biosynthesis in post-harvest fruits of a full red-type kiwifruit (Actinidia arguta) 'Jinhongguan'

Anthocyanin is the main component of pigment in red-fleshed kiwifruit. 'Jinhongguan' is a new cultivar of Actinidia arguta with red peel and flesh after harvest. However, the specific types of anthocyanin in the 'Jinhongguan' fruit and its biosynthesis pathways remain largely unknown. Here, the total anthocyanin content in the fruit color conversion process was determined. The results showed that total anthocyanin content increased with the deepening color of the peel and flesh. To identify the genes related to anthocyanin biosynthesis and the types of anthocyanins in the 'Jinhongguan' fruit, a combined analysis of transcriptome and anthocyanin-targeted metabolome was carried out. A total of 5751 common differentially expressed genes (DEGs) at different stages of peel and flesh were identified, of which 2767 were common up-DEGs and 2976 were common down-DEGs. KEGG and GO enrichment analyses showed that the common up-DEGs were significantly enriched in anthocyanin synthesis-related pathways, suggesting some up-DEGs are involved in anthocyanin biosynthesis. In total, 29 metabolites were detected in the flesh by anthocyanin-targeted metabolome. Among these, nine were differential accumulation metabolites (DAMs) in comparison to red flesh vs green flesh. Six DAMs were up-regulated, with five of them were cyanidins. The content of cyanidin-3-O-galactoside was much higher than that of other DAMs, making it the main pigment in 'Jinhongguan'. Moreover, a total of 36 anthocyanin synthesis-related structural genes, 27 MYB transcription factors (TFs), 37 bHLH TFs and 9 WDR TFs were screened from the common DEGs. Correlation analysis of transcriptome and metabolome revealed that 9 structural genes, 6 MYB TFs, 6 bHLH TFs and 1 WDR TF were significantly associated with cyanidin-3-O-galactoside. Further, qRT-PCR analysis demonstrated that structural genes (AaPAL3, Aa4CL3, AaCHS2/3/8/9/11, AaDFR1/2, AaANR1, UFGT3a and UFGT6b) and TFs (MYB108, bHLH30, bHLH94-1 and WD43) play important roles in cyanidin biosynthesis. Overall, this study identified cyanidin-3-O-galactoside as the main anthocyanin type and revealed key candidate genes of red coloration of post-harvest fruit in Actinidia arguta. These findings provided new insights into the color formation mechanism of post-harvest fruit and offered a theoretical basis for color regulation in kiwifruit.

An effective method for establishing a regeneration and genetic transformation system for Actinidia arguta

The all-red A. arguta (Actinidia arguta) is an anthocyanin-rich and excellent hardy fruit. Many studies have focused on the green-fleshed A. arguta, and fewer studies have been conducted on the all-red A. arguta. Here we reported a regeneration and Agrobacterium-mediated transformation protocol by using leaves of all-red A. arguta as explants. Aseptic seedling leaves of A. arguta were used as callus-inducing materials. MS medium supplemented with 0.3 mg·L-1 2,4-D and 1.0 mg·L-1 BA was the optimal medium for callus induction of leaves, and medium supplemented with 3 mg·L-1 tZ and 0.5 mg·L-1 IAA was optimal for adventitious shoot regeneration. The best proliferation medium for adventitious buds was MS + 1.0 mg·L-1 BA + 0.3 mg·L-1 NAA. The best rooting medium was 1/2MS + 0.7 mg·L-1 IBA with a 100% rooting rate. For the red flesh hardy kiwi variety 'Purpurna Saduwa' (A. arguta var. purpurea), leaves are receptors for Agrobacterium (EHA105)-mediated transformation. The orthogonal experiment was used for the optimization of each genetic transformation parameter and the genetic transformation of the leaves was 21% under optimal conditions. Our study provides technical parameters for applying genetic resources and molecular breeding of kiwifruit with red flesh.

Transcriptomic and Metabolomic Profiles Provide Insights into the Red-Stipe Symptom of Morel Fruiting Bodies

The cultivation of true morels (Morchella spp., Morchellaceae, Ascomycota) has rapidly expanded in recent years, especially in China. Red stipe is a symptom wherein the stipe of morel fruiting bodies becomes red-gray, resulting in the gradual death of the affected fruiting bodies. The impact of red-stipe symptom occurrence on the development and nutritional quality of morel fruiting bodies remains unclear. Herein, morel ascocarps with the red-stipe symptom (R) and normal (N), artificially cultivated in the Fujian Province of China, were selected for the transcriptome and metabolome analysis to study the physiological and biochemical responses of morel fruiting bodies to the red-stipe symptom. Transcriptome data revealed several differentially expressed genes between the R and N groups significantly enriched in the tyrosine, riboflavin, and glycerophospholipid metabolism pathways. Similarly, the differentially accumulated metabolites were mainly assigned to metabolic pathways, including tyrosine, the biosynthesis of plant secondary metabolites, and the biosynthesis of amino acids. Moreover, the transcriptome and metabolome data combination revealed that tyrosine metabolism was the most enriched pathway, which was followed by ATP-binding cassette (ABC) transport, alanine, aspartate, and glutamate metabolism. Overall, the integration of transcriptomic and metabolomic data of M. sextelata affected by red-stipe symptoms identified several important genes, metabolites, and pathways. These findings further improve our understanding of the mechanisms underlying the red-stipe symptom development of M. sextelata and provide new insights into how to optimize its cultivation methods.

Interaction of AcMADS68 with transcription factors regulates anthocyanin biosynthesis in red-fleshed kiwifruit

In red-fleshed kiwifruit, anthocyanin pigmentation is a crucial commercial trait. The MYB-bHLH-WD40 (MBW) complex and other transcription factors regulate its accumulation. Herein, a new SEP gene, AcMADS68, was identified as a regulatory candidate for anthocyanin biosynthesis in the kiwifruit by transcriptome data and bioinformatic analyses. AcMADS68 alone could not induce the accumulation of anthocyanin both in Actinidia arguta fruit and tobacco leaves. However, in combination with AcMYBF110, AcMYB123, and AcbHLH1, AcMADS68 co-overexpression increased anthocyanin biosynthesis, whereas its silencing reduced anthocyanin accumulation. The results of the dual-luciferase reporter, firefly luciferase complementation, yeast two-hybrid and co-immunoprecipitation assays showed that AcMADS68 could interact with both AcMYBF110 and AcMYB123 but not with AcbHLH1, thereby co-regulating anthocyanin biosynthesis by promoting the activation of the target genes, including AcANS, AcF3GT1, and AcGST1. Moreover, AcMADS68 also could activate the promoter of AcbHLH1 surported by dual-luciferase reporter and yeast one-hybrid assays, thereby further amplifying the regulation signals from the MBW complex, thus resulting in enhanced anthocyanin accumulation in the kiwifruit. These findings may facilitate better elucidation of various regulatory mechanisms underlying anthocyanin accumulation and contribute to the quality enhancement of red-fleshed kiwifruit.

Comprehensive Analysis of Metabolome and Transcriptome in Fruits and Roots of Kiwifruit

Kiwifruit (Actinidia chinensis) roots instead of fruits are widely used as Chinese medicine, but the functional metabolites remain unclear. In this study, we conducted comparative metabolome analysis between root and fruit in kiwifruit. A total of 410 metabolites were identified in the fruit and root tissues, and of them, 135 metabolites were annotated according to the Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway. Moreover, 54 differentially expressed metabolites (DEMs) were shared in root and fruit, with 17 DEMs involved in the flavonoid pathway. Of the 17 DEMs, three flavonols (kaempferol-3-rhamnoside, L-Epicatechin and trifolin) and one dihydrochalcone (phloretin) showed the highest differences in the content level, suggesting that flavonols and dihydrochalcones may act as functional components in kiwifruit root. Transcriptome analysis revealed that genes related to flavonols and dihydrochalcones were highly expressed in root. Moreover, two AP2 transcription factors (TFs), AcRAP2-4 and AcAP2-4, were highly expressed in root, while one bHLH TF AcbHLH62 showed extremely low expression in root. The expression profiles of these TFs were similar to those of the genes related to flavonols and dihydrochalcones, suggesting they are key candidate genes controlling the flavonoid accumulation in kiwifruit. Our results provided an insight into the functional metabolites and their regulatory mechanism in kiwifruit root.

Metabolomic and transcriptomic exploration of the uric acid-reducing flavonoids biosynthetic pathways in the fruit of Actinidia arguta Sieb. Zucc

Flavonoids from Actinidia arguta Sieb. Zucc. can reduce uric acid in mice. However, the molecular basis of its biosynthesis is still unclear. In this paper, we used a combination of extensively targeted metabolomics and transcriptomics analysis to determine the types and differences of flavonoids in the fruit ripening period (August to September) of two main cultivated varieties in northern China. The ethanol extract was prepared, and the potential flavonoids of Chrysin (Flavone1), Rutin (Flavone2), and Daidzein (Flavone3) in Actinidia arguta Sieb. Zucc. were separated and purified by HPD600 macroporous adsorption resin and preparative liquid chromatography. The structure was identified by MS-HPLC, and the serum uric acid index of male Kunming mice was determined by an animal model test.125 flavonoids and 50 differentially regulated genes were identified. The contents of UA (uric acid), BUN (urea nitrogen), Cr (creatinine), and GAPDH in mouse serum and mouse liver glycogen decreased or increased in varying degrees. This paper reveals the biosynthetic pathway of uric acid-reducing flavonoids in the fruit of Actinidia arguta Sieb. Zucc., a major cultivar in northern China, provides valuable information for the development of food and drug homologous functional foods.

Integrated metabolomics and transcriptomic analysis of the flavonoid regulatory networks in Sorghum bicolor seeds

Background

The objective of this study was to reveal the flavonoid biosynthesis pathway in white (Z6), red (Z27) and black (HC4) seeds of the sweet sorghum (Sorghum bicolor) using metabolomics and transcriptomics, to identify different flavonoid metabolites, and to analyze the differentially expressed genes involved in flavonoid biosynthesis.

Results

We analyzed the metabolomics and transcriptomics data of sweet sorghum seeds. Six hundred and fifty-one metabolites including 171 flavonoids were identified in three samples. Integrated analysis of transcriptomics and metabolomics showed that 8 chalcone synthase genes (gene19114, gene19115, gene19116, gene19117, gene19118, gene19120, gene19122 and gene19123) involved in flavonoid biosynthesis, were identified and play central role in change of color. Six flavanone including homoeriodictyol, naringin, prunin, naringenin, hesperetin and pinocembrin were main reason for the color difference.

Conclusions

Our results provide valuable information on the flavonoid metabolites and the candidate genes involved in the flavonoid biosynthesis pathway in sweet sorghum seeds.

Metabolomic and Transcriptomic Analysis of Flavonoid Biosynthesis in Two Main Cultivars of Actinidia arguta Sieb.Zucc. Grown in Northern China

Actinidia arguta Sieb.Zucc. is a fruit that is rich in flavonoids. Nevertheless, details of flavonoid formation and the potential mechanism behind flavonoid biosynthesis have not previously been reported. In order to explore the biosynthetic regulation mechanism of flavonoids in A. arguta Sieb.Zucc., we conducted a combination of extensive targeted metabolite analysis and analyzed transcriptomes to determine the flavonoids present and the genes bound up with flavonoid biosynthesis in the two main cultivated varieties of A. arguta Sieb.Zucc. in Northern China. The maturity period is from August to September. A total of 118 flavonoids were found in fruits. Among them, 39 flavonoids were accumulated at significant levels after fruit ripening. Transcriptome analysis indicated that most flavonoid biosynthesis structural genes and certain regulatory genes exhibited differential expression between the two varieties. Correlation analysis of transcriptome and metabolite profiles showed that the ways of expression of 21 differentially expressed genes related to structure and regulation between the 2 varieties were more highly correlated with 7 flavonoids after fruit ripening. These results contribute to the development of A. arguta Sieb.Zucc. as a food and drug homologous functional food.

Metabolome and Transcriptome Analyses of Anthocyanin Accumulation Mechanisms Reveal Metabolite Variations and Key Candidate Genes Involved in the Pigmentation of Prunus tomentosa Thunb. Cherry Fruit

Prunus tomentosa Thunb. has excellent nutritional, economic, and ornamental values with different fruit color. The red coloration of fruit is determined by anthocyanin pigmentation, which is an attractive trait for consumers. However, the mechanisms underlying fruit color formation in the P. tomentosa cherry are not well understood. In this research, the pigmentation patterns in red-color P. tomentosa (RP) fruit and white-color P. tomentosa (WP) were evaluated. Anthocyanin content in matured RP fruit was significantly abundant compared with WP fruit. Metabolomic profiling revealed that pelargonidin 3-O-glucoside, cyanidin 3-O-rutinoside, and pelargonidin 3-O-rutinoside were the predominant anthocyanin compounds in the RP fruit, while, WP fruit had less anthocyanin compositions and lower level. Then, integrative analyses of transcriptome and metabolome identified 285 significant differentially expressed genes (DEGs) closely related to anthocyanin differentially expressed metabolites (DEMs). Among them, nine genes were involved in anthocyanin biosynthesis, transport and degradation pathway, including four biosynthesis genes (PtPAL1, PtDFR, PtANS, and PtUFGT), two transport genes (PtGST11, PtABC10), and three degradation genes (PtPOD1, PtPOD16, PtPOD73). Transcriptome data and real-time PCR showed that the transcript levels of biosynthesis and transport genes were significantly higher in RP than in WP, especially PtANS, PtUFGT, and PtGST11, suggesting they may play key roles in red-colored fruit formation. Meanwhile, the degradation-related genes PtPOD1/16/73 took on exactly opposite trend, suggesting their potential effects on anthocyanin degradation. These results provide novel insights into color patterns formation mechanisms of cherries fruit, and the candidate key genes identified in anthocyanin biosynthesis, transport and degradation may provide a valuable resource for cherry breeding research in future.

Novel Role of AaMYBC1 in Regulating Actinidia arguta Vine Architecture by Elongating Internode Based on Multi-Omics Analysis of Transgenic Tobacco

The internode length affects the status of fruiting branches and shapes the vine architecture. MYB TFs (transcription factors) have been widely studied and reported to control many biological processes including secondary metabolism, abiotic stresses, growth and development, etc. However, the roles of MYB TFs in regulating internode length remain poorly understood. Here, we demonstrated that a secondary metabolism-related R2R3-MYB TF AaMYBC1 from Actinidia arguta was involved in the regulation of internode length by combined analysis of transcriptome and metabolome of transgenic tobacco plants. The metabolome analysis of OE (over-expressed tobacco) and WT (wild-typed tobacco) showed that there were a total of 1000 metabolites, 176 of which had significant differences. A key metabolite pme1651 annotated as indole 3-acetic acid belonged to phytohormone that was involved in internode length regulation. The RNA-seq analysis presented 446 differentially expressed genes (DEGs) between OE and WT, 14 of which were common DEGs in KEGG and GO enrichment. Through the combined analysis of metabolome and transcriptome in transgenic and wild-type tobacco, three key genes including two SAUR and a GH3 gene were possibly involved in internode elongation. Finally, a regulatory module was deduced to show the role of AaMYBC1 in internode elongation. Our results proposed a molecular mechanism of AaMYBC1 regulating internode length by mediated auxin signaling, implying the potential role in regulating the vine architecture.

Integrated LC-MS/MS and Transcriptome Sequencing Analysis Reveals the Mechanism of Color Formation During Prickly Ash Fruit Ripening

Prickly ash peel is one of the eight major condiments in China and is widely used in cooking because of its unique fragrance and numbing taste. The color of prickly ash fruit is the most intuitive quality that affects consumer choice. However, the main components and key biosynthetic genes responsible for prickly ash fruit color have not yet been determined. To better understand the biosynthetic mechanisms and accumulation of prickly ash fruit color components, we performed an integrated transcriptomic and metabolomic analysis of red and green prickly ash fruit at different growth periods. The transcriptome analysis identified 17,269 differentially expressed genes (DEGs) between fruit of red and green prickly ash: 7,236 upregulated in green fruit and 10,033 downregulated. Liquid chromatography tandem mass spectrometry (LC-MS/MS) identified 214 flavonoids of 10 types. Flavonoids and flavonols are the main flavonoids in prickly ash, and the total flavonoid content of red prickly ash is higher than that of green prickly ash. Comprehensive analysis showed that the main colored metabolites that differed between green and red prickly ash were cyanidin-3-O-galactoside and cyanidin-3-O-glucoside, and differences in the contents of these metabolites were due mainly to differences in the expression of ANS and UFGT. Our results provide insight into the mechanisms underlying color differences in red and green prickly ash and will be useful for improving the quality of prickly ash fruit.

Integrative analyses of metabolome and transcriptome reveals metabolomic variations and candidate genes involved in sweet cherry (Prunus avium L.) fruit quality during development and ripening

Sweet cherry (Prunus avium L.), one of the most appreciated and most important commercial temperate fruits, has high sensory quality and nutritional value. Investigating its metabolic variations provides valuable information on the formation of fruit quality. In this study, widely targeted LC-MS/MS based metabolomics was used to identify and quantify metabolic changes during 'Black Pearl' sweet cherry development and ripening. A total of 263 significant differentially expressed metabolites (DEMs) were detected during the four fruit-development stages. Significant differences were observed in the composition and content of compounds in the four stages of cherry development, especially sugars, organic acids, and flavonoids. Moreover, transcriptome analysis provided a molecular basis for metabolic variations during fruit development. A total of 6724 significant differentially expressed genes (DEGs) were identified. Further correlation analysis of major DEMs and DEGs showed that 19 key DEGs were involved in sugar metabolism, 23 key DEGs in organic acid metabolism, and 13 key DEGs in flavonoid metabolism. The upregulated genes involved in the flavonoid pathway probably play an important role in regulating the rapid increase of anthocyanin content during fruit development. These comprehensive analysis data provide a better understanding to improve fruit quality traits based on molecular and metabolic levels.

AvNAC030, a NAC Domain Transcription Factor, Enhances Salt Stress Tolerance in Kiwifruit

Kiwifruit (Actinidia chinensis Planch) is suitable for neutral acid soil. However, soil salinization is increasing in kiwifruit production areas, which has adverse effects on the growth and development of plants, leading to declining yields and quality. Therefore, analyzing the salt tolerance regulation mechanism can provide a theoretical basis for the industrial application and germplasm improvement of kiwifruit. We identified 120 NAC members and divided them into 13 subfamilies according to phylogenetic analysis. Subsequently, we conducted a comprehensive and systematic analysis based on the conserved motifs, key amino acid residues in the NAC domain, expression patterns, and protein interaction network predictions and screened the candidate gene AvNAC030. In order to study its function, we adopted the method of heterologous expression in Arabidopsis. Compared with the control, the overexpression plants had higher osmotic adjustment ability and improved antioxidant defense mechanism. These results suggest that AvNAC030 plays a positive role in the salt tolerance regulation mechanism in kiwifruit.

Effects of Anthocyanins on Vascular Health

Cardiovascular disorders are leading mortality causes worldwide, often with a latent evolution. Vascular health depends on endothelial function, arterial stiffness, and the presence of atherosclerotic plaques. Preventive medicine deserves special attention, focusing on modifiable cardiovascular risk factors, including diet. A diet rich in fruits and vegetables has well-known health benefits, especially due to its polyphenolic components. Anthocyanins, water-soluble flavonoid species, responsible for the red-blue color in plants and commonly found in berries, exert favorable effects on the endothelial function, oxidative stress, inhibit COX-1, and COX-2 enzymes, exert antiatherogenic, antihypertensive, antiglycation, antithrombotic, and anti-inflammatory activity, ameliorate dyslipidemia and arterial stiffness. The present review aims to give a current overview of the mechanisms involved in the vascular protective effect of anthocyanins from the human diet, considering epidemiological data, in vitro and in vivo preclinical research, clinical observational, retrospective, intervention and randomized studies, dietary and biomarker studies, and discussing preventive benefits of anthocyanins and future research directions.

Transcriptome Co-Expression Network Analysis Identifies Key Genes and Regulators of Sweet Cherry Anthocyanin Biosynthesis

Anthocyanin is the key factor that results in the attractive color of sweet cherry fruits. However, information regarding sweet cherry coloration and the potential mechanisms underlying anthocyanin biosynthesis is limited. In this study, we found that the anthocyanin accumulation varied in sweet cherry flesh and peel, while the anthocyanin content increased sharply in the dark red (DR) stage. Correlations between anthocyanin concentrations and RNA sequencing (RNA-seq), constructed with Weighted Gene Co-Expression Network Analysis (WGCNA), indicated that two structural genes (Pac4CL2, PacANS) and 11 transcription factors (PacbHLH13/74, PacDIV, PacERF109/115, PacGATA8, PacGT2, PacGTE10, PacMYB308, PacPosF21, and PacWRKY7) had similar expression patterns with the changes in anthocyanin content. Additionally, real-time PCR verified all of these gene expression patterns and revealed that PacANS exhibited the highest transcription level. In order to search for potential regulators for anthocyanin biosynthesis, a dual-luciferase assay was performed to investigate the regulatory activities of 11 transcription factors on the PacANS promoter. The results revealed that two novelty bHLHs, PacbHLH13 and PacbHLH74, can trans-activate the PacANS promoter and they might be the candidate genes for regulating anthocyanin synthesis in sweet cherry fruits. The present findings provide a novel viewpoint with regard to anthocyanin biosynthesis mechanisms and the regulatory transcriptional network of fruit coloration in sweet cherries.

Transcriptional Regulation of Anthocyanin Synthesis by MYB-bHLH-WDR Complexes in Kiwifruit (Actinidia chinensis)

The anthocyanin synthetic pathway is regulated centrally by an MYB-bHLH-WD40 (MBW) complex. Anthocyanin pigmentation is an important fruit quality trait in red-fleshed kiwifruit; however, the underlying regulatory mechanisms involving the MBW complex are not well understood. In this study, one R2R3MYB (AcMYBF110 expressed in fruit characteristically), one bHLH (AcbHLH1), two upstream regulators of AcbHLH1 (AcbHLH4 and AcbHLH5), and one WDR (AcWDR1) are characterized as being involved in the regulation of anthocyanin synthesis in kiwifruit. AcMYBF110 plays an important role in the regulation of anthocyanin accumulation by specifically activating the promoters of several anthocyanin pathway genes including AcCHS, AcF3'H, AcANS, AcUFGT3a, AcUFGT6b, and AcGST1. Coexpression of AcbHLH1, AcbHLH4, or AcbHLH5 together with AcMYBF110 induces much greater anthocyanin accumulation in both tobacco leaves and in Actinidia arguta fruit compared with AcMYBF110 alone. Moreover, this activation is further enhanced by adding AcWDR1. We found that both AcMYBF110 and AcWDR1 interact with all three AcbHLH factors, while AcMYBF110 also interacts with AcWDR1 to form three different MBW complexes that have different regulatory roles in anthocyanin accumulation of kiwifruit. The AcMYBF110-AcbHLH1-AcWDR1 complex directly targets the promoters of anthocyanin synthetic genes. Other features of the regulatory pathways identified include promotion of AcMYBF110, AcbHLH1,and AcWDR1 activities by this MBW complex, providing for both reinforcement and feedback regulation, whereas the AcMYBF110-AcbHLH4/5-AcWDR1 complex is indirectly involved in the regulation of anthocyanin synthesis by activating the promoters of AcbHLH1 and AcWDR1 to amplify the regulation signals of the first MBW complex.

Study of the Relationship between Leaf Color Formation and Anthocyanin Metabolism among Different Purple Pakchoi Lines

Purple pakchoi (Brassica rapa ssp. Chinensis) is particularly appreciated due to its high edible quality and ornamental value, but there are few studies on the underlying mechanisms of leaf color formation. To comprehensively assess the differences in purple formation in pakchoi, four lines of pakchoi with different purple leaves were used in this experiment to determine the pigment content and to investigate the distribution and components of anthocyanin using LCMS (Liquid Chromatography Mass Spectrometry) and leaf cross-sections. Moreover, the expression levels of anthocyanin synthesis-related genes in four lines were calculated by qRT-PCR. The results showed that three new purple lines rich in anthocyanin and of high-quality were bred, and the anthocyanin were mainly distributed in both the upper epidermis and lower epidermis of leaves. Thirteen anthocyanin components were separated and identified, all the anthocyanins were acylated and glycosylated cyanidins; the main anthocyanins in purple pakchoi were a diacylated form of cyanidin 3-trans-(feruloyl)diglucoside-5-(malonyl)glucoside. Both the ratio of non-aromatic acylated cyanidin to aromatic acylated cyanidin and the ratio of anthocyanin content to chlorophyll content were responsible for the color formation in different purple pakchoi lines. When the ratio was high, the leaf appeared reddish purple, and when the ratio was low, the leaf appeared deep purple, even blackish purple. The expression level of BrF3H was significantly correlated with the content of anthocyanin through the correlation coefficient, which was speculated to be the main anthocyanin synthesis-related gene resulting in color differences among the four purple pakchoi lines. These results will enhance our understanding for the cultivation of new purple pakchoi varieties.

Dynamic Changes of Phenolic Compounds and Their Associated Gene Expression Profiles Occurring during Fruit Development and Ripening of the Donghong Kiwifruit

The newly released Donghong kiwifruit is a promising commercial cultivar. The dynamic changes of major phenolic compounds (flavonols, flavanols, phenolic acids, and anthocyanins) during the representative stages of fruit development and ripening of the Donghong kiwifruit were determined by high-performance liquid chromatography. The corresponding time-course transcriptional changes were evaluated using the combined analysis of RNA-Seq and qRT-PCR. The most predominant phenolic compound in the Donghong kiwifruit was epicatechin. Cyanidin 3-O-[2-O-(β-xylosyl)-β-galactoside] and cyanidin 3-O-β-galactoside were two essential anthocyanins detected. Candidate genes and pathways involved in phenolic compounds biosynthesis were highlighted. The structural genes (AcLDOX2, Ac5GGT1, and Ac5AT2) and the transcription factor (bHLH74-2) were strongly associated with anthocyanin biosynthesis. AcMYB4-1 may be a novel transcription factor that reduces anthocyanin accumulation. Results from the study may be a very useful supplement to current knowledge of molecular mechanisms to elucidate coloration in the red-fleshed kiwifruit and could help breeders modify the kiwifruit germplasm.

The proanthocyanin-related transcription factors MYBC1 and WRKY44 regulate branch points in the kiwifruit anthocyanin pathway

The groups of plant flavonoid metabolites termed anthocyanins and proanthocyanins (PA) are responsible for pigmentation in seeds, flowers and fruits. Anthocyanins and PAs are produced by a pathway of enzymes which are transcriptionally regulated by transcription factors (TFs) that form the MYB-bHLH-WD40 (MBW) complex. In this study, transcriptomic analysis of purple-pigmented kiwifruit skin and flesh tissues identified MYBC1, from subgroup 5 of the R2R3 MYB family, and WRKY44 (highly similar to Arabidopsis TTG2) as candidate activators of the anthocyanin pathway. Transient over-expression of MYBC1 and WRKY44 induced anthocyanin accumulation in tobacco leaves. Dual luciferase promoter activation assays revealed that both MYBC1 and WRKY44 were able to strongly activate the promoters of the kiwifruit F3'H and F3'5'H genes. These enzymes are branch points of the pathway which specifies the type of anthocyanin accumulated. Stable over-expression of MYBC1 and WRKY44 in kiwifruit calli activated the expression of F3'5'H and PA-related biosynthetic genes as well as increasing levels of PAs. These results suggest that while previously characterised anthocyanin activator MYBs regulate the overall anthocyanin biosynthesis pathway, the PA-related TFs, MYBC1 and WRKY44, more specifically regulate key branch points. This adds a layer of regulatory control that potentially balances anthocyanin and PA levels.

Freezing Tolerance and Expression of β-amylase Gene in Two Actinidia arguta Cultivars with Seasonal Changes

Low temperature causes injuries to plants during winter, thereby it affects kiwi fruit quality and yield. However, the changes in metabolites and gene expression during cold acclimation (CA) and deacclimation (DA) in kiwi fruit remain largely unknown. In this study, freezing tolerance, carbohydrate metabolism, and β-amylase gene expression in two Actinidia arguta cv. "CJ-1" and "RB-3" were detected from CA to DA stages. In all acclimation stages, the "CJ-1" was hardier than "RB-3" and possessed lower semi-lethal temperature (LT50). Furthermore, "CJ-1" had a more rapid acclimation speed than "RB-3". Changes of starch, β-amylase, and soluble sugars were associated with freezing tolerance in both cultivars. Starch contents continued to follow a declining trend, while soluble sugars contents continuously accumulated in both cultivars during CA stages (from October to January). To investigate the possible molecular mechanism underlying cold response in A. arguta, in total, 16 AcBAMs genes for β-amylase were identified in the kiwi fruit genome. We carried out localization of chromosome, gene structure, the conserved motif, and the analysis of events in the duplication of genes from AcBAMs. Finally, a strong candidate gene named AaBAM3 from AcBAMs was cloned in Actinidia arguta (A. arguta), The real-time qPCR showed that AaBAM3 gene expression in seasonal changes was consistent with changes of soluble sugars. These results reveal that AaBAM3 may enhance the freezing tolerance of A. arguta through increasing soluble sugar content.

Chalcone Synthase-Encoding AeCHS is Involved in Normal Petal Coloration in Actinidia eriantha

Studies on anthocyanin biosynthesis have been mainly concentrated on the fruit, whereas few have focused the mechanism of flower coloration in kiwifruit. Here, we report that the structural gene, AeCHS, is involved in anthocyanin accumulation and indispensable for normal petal coloration in Actinidia eriantha. Petals from three different species including Actinidia eriantha (red petals), Actinidia hemsleyana (light pink petals) and Actinidia arguta (white petals) were selected for anthocyanin determination and gene expression analysis. The anthocyanin components in A. eriantha were significantly higher than in A. hemsleyana or A. arguta. Consistently, gene expression profiles suggested that AeCHS expression in A. eriantha was higher than in A. hemsleyana or A. arguta. Cluster analysis showed that AeCHS was clustered into a single group and distinctly separated from other genes, indicating the expression pattern of AeCHS gene was different from any other. Additionally, correlation analysis revealed AeCHS expression significantly correlated with anthocyanin content. The complete coding sequence of AeCHS was cloned from petals of A. eriantha 'Zaoxu', showing the length of AeCHS was 1170 bp encoding a protein of 389 amino acids. AeCHS was located in the cytoplasm, indicating it is indeed a structural gene involved in anthocyanin biosynthesis. AeCHS silencing performed by infiltration grafting-mediated virus-induced gene silencing (VIGS) reduced petal anthocyanin content and bleached red petals in A. eriantha. Our results confirm a crucial role of AeCHS in anthocyanin biosynthesis and accumulation in A. eriantha petals; furthermore, they offer important basic information and constitute a reference point for further research.

MicroRNA858-mediated regulation of anthocyanin biosynthesis in kiwifruit (Actinidia arguta) based on small RNA sequencing

As important regulators, miRNAs could play pivotal roles in regulation of fruit coloring. Actinidia arguta is a newly emerged fruit tree with extensively application prospects. However, miRNAs involved in A. arguta fruit coloring are unknown. In this study, A. arguta fruit were investigated at three developmental stages by small RNAs high-throughput sequencing. A total of 482 conserved miRNAs corresponding to 526 pre-miRNAs and 581 novel miRNAs corresponding to 619 pre-miRNAs were grouped into 46 miRNA families. Target gene prediction and analysis revealed that miR858, a strongly candidate miRNA, was involved in anthocyanin biosynthesis in which contributes to fruit coloring. The anthocyanin level was determined in three A. arguta cultivars by UPLC-MS/MS (ultra-performance liquid chromatography coupled with tandem mass spectrometry). In addition, qPCR (quantitative real-time PCR), cluster analysis were conducted as well as correlation analysis. All results were combined to propose a model in which describes an association of miRNA and anthocyanin biosynthesis in A. arguta. The data presented herein is the first report on miRNA profile analysis in A. arguta, which can provide valuable information for further research into the regulation of the miRNAs in anthocyanin biosynthesis and fruit coloring.

Overexpression of SmANS Enhances Anthocyanin Accumulation and Alters Phenolic Acids Content in Salvia miltiorrhiza and Salvia miltiorrhiza Bge f. alba Plantlets

Flavonoids play multiple roles in plant coloration and stress resistance and are closely associated with human health. Flavonoids and non-flavonoids (such as phenolic acids) are produced via the phenylpropanoid-derived pathway. Anthocyanidin synthase (ANS) catalyzes the synthesis of anthocyanins from leucoanthocyanidin in the flavonoids branched pathway. In this study, SmANS from Salvia miltiorrhiza was cloned and mainly localized in the endoplasmic reticulum (ER), plastids, Golgi, plasma membrane, and nucleus of tobacco epidermal cells, and was most highly expressed in purple petals in S. miltiorrhiza, whereas it showed almost no expression in white petals, green calyxes, and pistils in S. miltiorrhiza Bge f. alba. Overexpressed SmANS enhanced anthocyanin accumulation but reduced salvianolic acid B (SAB) and rosmarinic acid (RA) biosynthesis in S. miltiorrhiza and S. miltiorrhiza Bge f. alba plantlets, meanwhile, it restored the purple-red phenotype in S. miltiorrhiza Bge f. alba. These changes were due to reallocation of the metabolic flow, which was influenced by the SmANS gene. These findings indicate that SmANS not only plays a key role in anthocyanin accumulation in S. miltiorrhiza, but also acts as a “switch” for the coloration of S. miltiorrhiza Bge f. alba. This study provides baseline information for further research on flavonoids metabolism and improvement of anthocyanin or phenolic acid production by genetic engineering.

De Novo Assembled Transcriptome Analysis and Identification of Genic SSR Markers in Red-Flowered Strawberry

Red-flowered strawberry is a new ornamental flower derived from intergeneric hybridization (Fragaria × Potentilla). To date, few molecular markers have been reported for this plant. RNA sequencing provides a relatively fast and low-cost approach for large-scale detection of simple sequence repeats (SSRs). In the present study, we profiled the transcriptome of red-flowered strawberry by Illumina HiSeq 2500 to identify SSRs related to petal color. Based on 2 million clean reads of red and white flowers from red-flowered strawberry hybrids, we assembled 91,835 unigenes with an average length of 717 bp. After functional annotation and prediction, there were 47,058 unigenes; of these, 26,861 had a gene ontology annotation, with 14,264 SSR loci. Mononucleotide SSRs were the predominant repeat type (47.20%, n = 6724), followed by di- (32.50%, n = 4641), tri- (19.10%, n = 2729), tetra- (0.90%, n = 132), hexa- (0.2%, n = 21), and penta- (0.10%, n = 16) nucleotide repeats. The most frequent di-, tri-, and tetra-nucleotide repeats were AG/CT, AAG/CTT, and AAAG/CTTT, respectively. PCR amplification with 105 SSR primer pairs yielded four bands specific to red flowers, namely UgRFsr57622, UgRFsr94149, UgRFsr40142, and UgRFsr54608; corresponding 4 trait-specific markers were found to co-segregate with white and red flower color in hybrid population, demonstrating that the genic SSR marker is useful to discriminate between white and red flowers in strawberry. Markers to discriminate flower color in red-flowered strawberry will be useful for early selection of progeny and for breeding management.

Differential regulation of the anthocyanin profile in purple kiwifruit (Actinidia species)

Anthocyanins are a group of secondary metabolites that colour fruit and flowers orange, red, purple or blue depending on a number of factors, such as the basic structure, co-pigmentation, metal ion complexation and vacuolar pH. The biosynthesis of anthocyanin is regulated at the transcriptional level by a group of transcription factors, the MYB-bHLH-WD40 (MBW) complex. In this study, the purple colouration in several kiwifruit (Actinidia) species was identified and characterised as red cyanidin-based and blue delphinidin-based anthocyanins. The differential pigmentation in the skin and flesh can be attributed to the differential ratio of cyanidin and delphinidin derivatives accumulated in the total anthocyanin profile. The expression of anthocyanin biosynthetic genes chalcone synthase (CHS), flavonoid 3-O-glucosyltransferase (F3GT), flavonoid 3'-hydroxylase (F3'H) and flavonoid 3'5'-hydroxylase (F3'5'H) is crucial for anthocyanin accumulation. However, the balance of expression of the F3'H and F3'5'H genes appears responsible for the ratio of cyanidin and delphinidin derivatives, while a lack of CHS, F3GT and MYB110 expression is responsible for a lack of total anthocyanins. The transcriptional regulation of the F3'H and F3'5'H promoters by the R2R3 MYB transcription factor MYB110 is markedly different in tobacco transient assays. When kiwifruit MYB10 or MYB110 are over-expressed in Actinidia chinensis both cyanidin-based and delphinidin-based anthocyanins are elevated, but F3'H and F3'5'H genes are not strongly correlated with MYB expression. These results suggest that the core kiwifruit anthocyanin pathway genes are dependent on characterised MYB transcription factors, while other regulatory proteins are more directly responsible for the expression of the F3'H and F3'5'H genes.

Biochemical and Functional Characterization of Anthocyanidin Reductase (ANR) from Mangifera indica L

Mango (Mangifera indica L.) is abundant in proanthocyanidins (PAs) that are important for human health and plant response to abiotic stresses. However, the molecular mechanisms involved in PA biosynthesis still need to be elucidated. Anthocyanidin reductase (ANR) catalyzes a key step in PA biosynthesis. In this study, three ANR cDNAs (MiANR1-1,1-2,1-3) were isolated from mango, and expressed in Escherichia coli. In vitro enzyme assay showed MiANR proteins convert cyanidin to their corresponding flavan-3-ols, such as (−)-catechin and (−)-epicatechin. Despite high amino acid similarity, the recombinant ANR proteins exhibited differences in enzyme kinetics and cosubstrate preference. MiANR1-2 and MiANR1-3 have the same optimum pH of 4.0 in citrate buffer, while the optimum pH for MiANR1-1 is pH 3.0 in phosphate buffer. MiANR1-1 does not use either NADPH or NADH as co-substrate while MiANR1-2/1-3 use only NADPH as co-substrate. MiANR1-2 has the highest Km and Vmax for cyanidin, followed by MiANR1-3 and MiANR1-1. The overexpression of MiANRs in ban mutant reconstructed the biosynthetic pathway of PAs in the seed coat. These data demonstrate MiANRs can form the ANR pathway, leading to the formation of two types of isomeric flavan-3-ols and PAs in mango.

Characterization of Gene Expression Profile, Phenolic Composition, and Antioxidant Capacity in Red-Fleshed Grape Berries and Their Wines

Gene expression profile, phenolic composition, and antioxidant capacity were evaluated in red-fleshed berries and their wines (RF berries and wines) from new grape genotypes. Transcriptomic analysis revealed that ten metabolic pathways involved in polyphenol synthesis and catabolism were significantly altered, and 13 genes related to the biosynthesis and transport of phenolics were largely upregulated in RF berries compared to that of Cabernet Sauvignon (CS). Expression of MybA1 was associated with anthocyanin accumulation in red flesh. Additionally, RF berries and wines contained higher concentrations of total anthocyanins, phenols, flavonoids, and proanthocyanidins than those in CS berries and wine. Particularly, diglucosides of malvidin, peonidin, delphinidin, and cyanidin were present in red flesh and RF wines, but they were undetectable or present at very low concentrations in CS flesh and wine. Cinnamic acid and ferulic acid were clearly increased in the RF wines compared to those in the CS wine. Additionally, the RF wines had higher antioxidant capacity than that in the CS wine, and total anthocyanin content was significantly correlated to antioxidant capacity. This research provides insight into the mechanisms underlying grape flesh coloration and the composition of phenolic compounds in RF berries and wines.

Combined Analysis of the Fruit Metabolome and Transcriptome Reveals Candidate Genes Involved in Flavonoid Biosynthesis in Actinidia arguta

To assess the interrelation between the change of metabolites and the change of fruit color, we performed a combined metabolome and transcriptome analysis of the flesh in two different Actinidia arguta cultivars: "HB" ("Hongbaoshixing") and "YF" ("Yongfengyihao") at two different fruit developmental stages: 70d (days after full bloom) and 100d (days after full bloom). Metabolite and transcript profiling was obtained by ultra-performance liquid chromatography quadrupole time-of-flight tandem mass spectrometer and high-throughput RNA sequencing, respectively. The identification and quantification results of metabolites showed that a total of 28,837 metabolites had been obtained, of which 13,715 were annotated. In comparison of HB100 vs. HB70, 41 metabolites were identified as being flavonoids, 7 of which, with significant difference, were identified as bracteatin, luteolin, dihydromyricetin, cyanidin, pelargonidin, delphinidin and (-)-epigallocatechin. Association analysis between metabolome and transcriptome revealed that there were two metabolic pathways presenting significant differences during fruit development, one of which was flavonoid biosynthesis, in which 14 structural genes were selected to conduct expression analysis, as well as 5 transcription factor genes obtained by transcriptome analysis. RT-qPCR results and cluster analysis revealed that AaF3H, AaLDOX, AaUFGT, AaMYB, AabHLH, and AaHB2 showed the best possibility of being candidate genes. A regulatory network of flavonoid biosynthesis was established to illustrate differentially expressed candidate genes involved in accumulation of metabolites with significant differences, inducing red coloring during fruit development. Such a regulatory network linking genes and flavonoids revealed a system involved in the pigmentation of all-red-fleshed and all-green-fleshed A. arguta, suggesting this conjunct analysis approach is not only useful in understanding the relationship between genotype and phenotype, but is also a powerful tool for providing more valuable information for breeding.