Next Article in Journal
Eugenol β-Amino/β-Alkoxy Alcohols with Selective Anticancer Activity
Next Article in Special Issue
ABA Speeds Up the Progress of Color in Developing F. chiloensis Fruit through the Activation of PAL, CHS and ANS, Key Genes of the Phenylpropanoid/Flavonoid and Anthocyanin Pathways
Previous Article in Journal
Molecular Determinants of Neuroblastoma
Previous Article in Special Issue
A 69 kbp Deletion at the Berry Color Locus Is Responsible for Berry Color Recovery in Vitis vinifera L. Cultivar ‘Riesling Rot’
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 7
10.3390/ijms23073766
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Analysis of Apples in High-Temperature Treatments Reveals a Role of MdLBD37 in the Inhibition of Anthocyanin Accumulation

by
Yu-Feng Bu
1,†,
Shuo Wang
1,†,
Chen-Zhiyu Li
1,
Yue Fang
1,
Ya Zhang
1,
Qing-Yu Li
2,
Hai-Bo Wang
3,
Xue-Sen Chen
1,* and
Shou-Qian Feng
1,*
1
National Key Laboratory of Crop Biology, National Research Center for Apple Engineering and Technology, College of Horticulture Science and Engineering, Shandong Agricultural University, Tai’an 271018, China
2
Yantai Academy of Agricultural Science, Yantai 264000, China
3
Shandong Institute of Pomology, Tai’an 271000, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(7), 3766; https://doi.org/10.3390/ijms23073766
Submission received: 28 February 2022 / Revised: 20 March 2022 / Accepted: 25 March 2022 / Published: 29 March 2022
(This article belongs to the Special Issue Genetic Regulation of Fruit Pigmentation)
Browse Figures
Review Reports Versions Notes

Abstract

:
Coloring in apple fruit due to anthocyanin accumulation is inhibited by high temperature; however, the underlying mechanism remains unclear. In the present study, total anthocyanin and cyanidin 3-galactoside contents were determined and compared between cv. ‘Redchief Delicious’ apple fruits at 25 °C and 35 °C treatments. The high temperature (35 °C) treatment substantially decreased total anthocyanin and cyanidin 3-galactoside contents. The transcriptomes of 25 °C- and 35 °C-treated apples were analyzed by high-throughput RNA sequencing. A total of 8354 differentially expressed genes (DEGs) were detected at four time points corresponding to the two temperature treatments. The up-regulated DEGs were annotated using GO as well as KEGG databases. A network module of 528 genes (including 21 transcription factors) most associated with the total anthocyanin and cyanidin 3-galactoside contents was constructed by weighted correlation network analysis (WGCNA). In the WGCNA module, we unearthed a LOB domain-containing gene designated as MdLBD37. The expression of MdLBD37 was sharply up-regulated by high temperature and negatively correlated with the total anthocyanin and cyanidin 3-galactoside contents. Overexpression of MdLBD37 in apple fruit and calli decreased the expression of anthocyanin biosynthetic genes, such as MdCHI, MdCHS, MdF3H, MdANS, MdDFR, and MdUFGT, along with anthocyanin accumulation. Our results suggested that MdLBD37 significantly influenced the high-temperature inhibition of anthocyanin accumulation in apples. The findings shed more light on the mechanism of anthocyanin inhibition during high-temperature stress in apples.

1. Introduction

Color is an important index to measure the appearance quality of an apple (Malus domestica Borkh.), which largely determines its commodity value and market popularity [1]. Fruit color is closely related to the anthocyanin components and contents in the pericarp. The anthocyanin content in the peel directly affects fruit coloring and usually differs across the variety [2,3,4]. Generally, the anthocyanin content of a variety is mainly determined by genetic factors [5], but environmental factors, including temperature, are other important regulators [4,6,7].
HT (high temperature) inhibits anthocyanin accumulation in apples, affecting fruit coloration and economic value [8,9,10]. The accelerated urbanization and industrialization, and the corresponding release of greenhouse gases, have boosted global warming [11]. In recent years, the frequency of high-temperature weather has especially increased, thereby affecting larger areas [12]. Continuous high-temperature weather during the summer and autumn is making the phenomenon of poor apple coloration increasingly serious [8]. As a result, unraveling the molecular foundation of high-temperature inhibition of apple coloring is crucial.
Several studies reported that high temperature inhibits anthocyanin synthesis by reducing the activity of associated enzymes and the expression of anthocyanin-related genes [13,14]. Thus far, anthocyanin biosynthetic genes, including PAL (phenylalanine ammonia-lyase), CHS (chalcone synthase), CHI (chalcone isomerase), F3H (flavanone 3-hydroxylase), DFR (dihydroflavonol 4-reductase), and UFGT (uridine diphosphate-glucose: flavonoid 3-O-glucosyltransferase), have been isolated and identified from apple [4,15,16]. The expression of these structural genes is regulated by different kinds of transcription factors (TFs), such as MYB, bHLH, LBD, NAC, and WD40 [17,18,19,20,21]. MdMYB1 was identified as the key gene which regulates anthocyanin synthesis in apple peel [22]. Recent studies have found that the expression of MdMYB1, MdCHI, MdF3H, and MdUFGT was inhibited by high temperatures [8,9]. However, the regulatory mechanism of high-temperature inhibition of anthocyanin accumulation in apples has not been fully understood.
In the present study, the total anthocyanin and cyanidin 3-galactoside contents and transcriptome in the peel of ‘Redchief Delicious’ apples treated at 25 °C and 35 °C were estimated and compared to understand the mechanism of high-temperature inhibition of apple coloration.

2. Results

2.1. Decrease in Anthocyanin Accumulation in Apple Due to HT Treatment

The coloration of ‘Redchief Delicious’ apple fruit was markedly reduced by the HT (35 °C) treatment; hence, no obvious color changes were noted on the fruit peel throughout the treatment period. In contrast, red coloration of the 25 °C-treated apple fruit was first detected 3 days after initiating the treatment, following which the color continued to intensify until the fruit peel was completely red (Figure 1A). Consistent with the differences in coloration, corresponding contents of total anthocyanin and cy-3-gala continued to increase in the 25 °C-treated fruit over the treatment period, whereas they steadily decreased in the 35 °C-treated fruit (Figure 1B,C). Thus, high-temperature exposure repressed anthocyanin synthesis in the apple fruit.

2.2. Gene Expression Differences between RT- and HT-Treated Apples

RNA-seq analysis of the 25 °C- and 35 °C-treated apple fruits (collected at four time points) generated a total of 153.1 Gb clean data with Q30 > 95%. The clean reads were mapped into reference genome sequences of apple (Table S2). The gene transcript levels were quantified by the RPKM values. A total of 8354 DEGs were detected from the four pairwise comparisons (absolute log2 fold change > 1, adjusted p < 0.05). Specifically, 3934, 3160, 2582, and 2752 DEGs were identified at T1, T2, T3, and T4 (3, 6, 9, and 12 days after temperature treatment) between 25 °C- and 35 °C-treated fruits (Table S3). Out of these, 1505, 1806, 1090, and 1781 genes (at T1, T2, T3, and T4, respectively) had down-regulated expression levels, while 2429, 1354, 1492, and 971 genes (corresponding to T1, T2, T3, and T4, respectively) showed up-regulated expression levels (Figure 2A). A Venn diagram presenting the distribution of DEGs revealed 344 DEGs that were common to all four comparisons. Additionally, 2021, 1189, 1226, and 1234 DEGs were specific to the 35 °C vs. 25 °C comparisons at T1, T2, T3, and T4, respectively (Figure 2B).
To assess the reliability of the RNA-seq data, 10 DEGs were chosen for a qRT-PCR assay. The trends in the expression-level changes determined by the qRT-PCR analysis were consistent with those in the RNA-seq data (Figure 2C). Thus, the RNA-seq data were confirmed as highly reliable and useful for identifying the functional genes in apple.

2.3. Functional Classification of DEGs

Most of the functionally annotated (via GO analysis) DEGs were associated with ‘phosphoric ester hydrolase activity’, ‘nucleus’, and ‘membrane-bounded organelle’ in the 35 °C vs. 25 °C comparisons at T1. At T2, most of the annotated genes were related to ‘catalytic activity’ in the 35 °C vs. 25 °C comparisons. In comparison, most of the DEGs were related to the ‘photosystem’ and ‘membrane’ in the 35 °C vs. 25 °C comparisons at T3 and T4, respectively (Figure 3). The KEGG results indicated that carbon metabolism was the most significantly enriched pathway of DEGs in the 35 °C vs. 25 °C comparisons at T1, T2, T3, and T4 (Figure 4).

2.4. Anthocyanin-Related DEGs Revealed by Analysis of Co-Expression Networks

A total of eleven WGCNA modules were identified in the analysis of co-expression networks (Figure 5A,B). Interestingly, the ‘green’ module comprising of 528 genes was most highly correlated with the total anthocyanin and cy-3-gala contents (Table S4). Twelve genes in the ‘green’ module were anthocyanin structural genes (encoding six CHSs, two F3Hs, one DFR, and two LDOXs), whereas the remaining genes were novel. The structural genes showed high expression within apple fruits at 25 °C treatment compared to those at 35 °C treatment (Figure 5D).
There were 21 TF genes in the ‘green’ module (Table S4). Most of these TF genes showed a down-regulation tendency at the 35 °C-treated apples compared to the 25 °C-treated fruits. However, homeodomain-like superfamily protein (MD02G1166700) and MdLBD37 (MD15G1294700) showed obvious up-regulated expressions in the 35 °C-treated apples compared to the 25 °C-treated fruits (Figure 5D).
Furthermore, we performed correlation analyses using the above 21 TF genes and 12 anthocyanin structural genes. A total of 98 notable correlations (correlation coefficient, R2 > 0.9) were detected between 15 TF genes and 12 anthocyanin structural genes. The results indicated that these TF genes might be important candidate anthocyanin regulators. Most of these TF genes showed positive correlations with anthocyanin structural genes, while homeodomain-like superfamily protein (MD02G1166700) and MdLBD37 (MD15G1294700) showed negative correlations with anthocyanin structural genes (Figure 5C). Next, we chose MdLBD37 to test its function in regulating anthocyanin synthesis in apples.

2.5. MdLBD37-Mediated Inhibition of Anthocyanin Synthesis in Apple

We injected the plasmids MdLBD37-TRV and MdLBD37-pRI into the pericarps of young (50 DAFB) and coloring (140 DAFB) fruits of ‘Otome’, respectively. Compared with the control, the transient silencing of MdLBD37 significantly induced the accumulation of anthocyanins around the injection site (Figure 6A–C). Consistently, anthocyanin biosynthetic genes, including MdCHI, MdCHS, MdF3H, MdDFR, and MdUFGT, were up-regulated around the injection site. On the contrary, overexpression of MdLBD37 significantly inhibited the accumulation of anthocyanins and the expression of MdCHI, MdCHS, MdF3H, MdDFR, and MdUFGT around the injection site (Figure 6D). This observation suggested that MdLBD37 inhibited anthocyanin synthesis in apples.
Furthermore, we also generated transgenic red-fleshed apple calli overexpressing MdLBD37 (MdLBD37-OE). The overexpression of MdLBD37 in red-fleshed calli weakened the calli red color and lowered the expression levels of anthocyanin biosynthetic genes MdCHI, MdCHS, MdF3H, MdDFR, and MdUFGT (Figure 6E–H). These results indicate that MdLBD37 might repress anthocyanin synthesis by inhibiting anthocyanin biosynthetic genes.

3. Discussion

Several studies have reported that fruit red color formation is weakened by a relatively higher ambient temperature [8,14,23]. Similarly, we found that high temperature significantly inhibited anthocyanin accumulation, resulting in the poor coloring of apples. Nonetheless, further investigations are required to illustrate the detailed mechanism of anthocyanin inhibition in apples at high-temperature stress. Transcriptome analysis of high-temperature stressed apple fruits was conducted in the present study to identify the genes involved in regulating anthocyanin production.
Conjoint analysis of transcriptome data by WGCNA can be used for finding the key genes responsible for the phenotypic difference. Recently, genes involved in anthocyanin-deficient yellow-skin somatic mutation in apples have been analyzed using WGCNA, and MdMYB10 was successfully detected as the candidate gene [24]. The present study identified one module (green) of 528 genes that significantly correlated with the total anthocyanin and cy-3-gala contents using the WGCNA method. The spatiotemporal difference of anthocyanin accumulation is dependent on the structural genes. In the ‘green’ gene module, 12 anthocyanin structural genes, including six CHSs, two F3Hs, one DFR, and two LDOXs, were identified. It was evident that the expression of these anthocyanin structural genes was positively related to the total anthocyanin and cy-3-gala contents. Compared to the RT-treated apples, the expression of these anthocyanin structural genes in HT-treated apples was significantly lower, indicating their key role in the high-temperature inhibition of anthocyanin biosynthesis. The observations are similar to the former studies showing that high temperature decreased anthocyanin accumulation by inhibiting anthocyanin synthesis-related enzyme activity and coding gene expression [8,9,14].
Negative regulatory genes exert a vital function in inhibiting anthocyanin synthesis in plants at high temperatures. For example, MdCOL4 was found to inhibit anthocyanin synthesis in apples at high temperatures [9]. In potatoes, the reduction of anthocyanin accumulation was caused by enhancing StMYB44 expression [25]. In our study, MdLBD37 in the ‘green’ gene module was identified as an important anthocyanin-negative gene in apples. We found that MdLBD37 inhibited anthocyanin synthesis and the expression of anthocyanin structural genes, and it was sharply induced by high temperatures. This indicated that MdLBD37 might be a new negative regulatory gene in the high-temperature regulation of anthocyanin accumulation in apples. Previously, AtLBD37, a homologous gene of MdLBD37 in Arabidopsis, was identified as an anthocyanin negative regulatory gene working together with an MBW complex, and it was up-regulated in response to N/NO3–signals [26]. However, the regulatory mechanism of MdLBD37-mediated high-temperature signal inhibition of anthocyanin synthesis has not previously been described, and this will be a very interesting topic for future research.
Anthocyanin synthesis is transcriptionally regulated by a set of TFs [27,28]. In the ‘green’ gene module, we found different TFs, including ERF, LBD, GRAS, G2-like, Dof, NAC, HD-ZIP, DBB, Trihelix, NF-YC, HSF, and MYB related members, some of which had been previously shown to regulate anthocyanin synthesis. For example, MdMYB1, encoding an R2R3-MYB TF, was necessary for anthocyanin synthesis in apple skin [22]. In addition, it was shown that HT inhibited the MdMYB1 gene. MdERF1B and MdERF3, encoding ERF TFs, and MdNAC029, encoding an NAC TF, have been shown to promote anthocyanin synthesis in apples [29,30,31]. It could be possible that TFs identified in the ‘green’ gene module might participate in regulating high-temperature effects on anthocyanin accumulation in apples.
Above all, anthocyanin biosynthesis in apples is inhibited at HT, mainly at the transcriptional level [7]. Conjoint analysis of transcriptome data by WGCNA revealed a module (green) that exhibited a significant correlation with the total anthocyanin and cy-3-gala contents. MdLBD37 in the ‘green’ gene module might be an important negative regulator located at the upstream of anthocyanin structural genes and play an important inhibitory role in anthocyanin accumulation at high temperatures (Figure 7). The results, therefore, have important implications for the molecular mechanism of anthocyanin regulation at high temperatures. We will focus on the functions of LBD TFs and other regulators identified in the WGCNA network in the future.

4. Materials and Methods

4.1. Collection and Temperature Treatment of the Apple Fruits

The fruits were reaped from apple (Malus domestica) cultivar ‘Redchief Delicious’ trees growing in an orchard in Yantai, Shandong, China (37°5′ N, 122°1′ E) at 126 DAFB (days after full bloom), when the peel began coloring. The collected fruits were exposed to 20,000 lx light intensity at 35 °C and 25 °C temperatures. We harvested fruit peels at 3, 6, 9, and 12 days of post-treatment, followed by rapid freezing within liquid nitrogen and preservation under −80 °C till further analysis.
Fruits were harvested from apple (Malus domestica) cultivar ‘Otome’ trees growing in the experimental farm of the Shandong Pomology Institute (Taian, Shandong, China, 36°11′ N, 117°61′ E) at 50 and 140 DAFB, for injection assays. Transient silencing and overexpression assays using fruits at 50 and 140 DAFB, respectively. ‘Purple 3’ red-fleshed apple calli were cultured on Murashig and Skoog (MS) medium at 25 °C in the dark. The calli were cultured for 15 days before being used for genetic transformation.

4.2. Anthocyanin Assessment

Apple fruit peels were ground in liquid nitrogen. After that, 0.5 g ground sample was subjected to 12 h 1% (v/v) HCl-methanol (10 mL) extraction in the dark at 4 °C. The resultant extract for each sample was centrifuged for 3-min at 6000× g followed by absorbance measurement at 510 and 700 nm using the UV-2450 spectrophotometer (Shimadzu, Kyoto, Japan). Total anthocyanin content was determined using the pH-differential method [32]. The resultant extract was filtered through a 0.2 µm polyethersulfone (PES) filter (Krackeler Scientific, Albany, NY, USA) and analyzed using a high-performance liquid chromatography (HPLC) system (LC2010C, Shimadzu, Japan). Cy-3-gala (cyanidin 3-galactoside) was identified with the HPLC system, followed by detection of its content using an ultraviolet detector at 520 nm and calibration curves derived from a cy-3-gala standard (Sigma-Aldrich, St. Louis, MO, USA). Three independent biological replicates were used for all the analyses.

4.3. RNA-Seq and Data Analysis

Total RNA was extracted using TRIzol reagent (Invitrogen, Carlsbad, CA, USA). Total RNA was treated with the RNase-Free DNase Set (Qiagen, Nasdaq, NY, USA). RNA quality was assessed using the Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). This was followed by cDNA preparation based on high-quality RNA. Thereafter, we adopted Illumina HiSeq™ 2000 system to sequence cDNA libraries. Low-quality reads were removed from raw RNA-seq reads. Afterward, clean reads were mapped to build a reference genome sequence for apple. Fragments Per Kilobase of transcript (FPKM) approach was utilized to determine gene expression levels [33]. Differentially expressed genes (DEGs) were detected using the DEGseq R package [34], with an adjusted p-value ≤ 0.05 and a fold-change ≥ 2 as the significance thresholds.

4.4. GO and KEGG Enrichment Analysis

DEGs were subjected to Gene Ontology (GO) annotation using the Blast2GO software. The GO terms of those DEGs showing adjusted p < 0.05 were significant [35]. We also conducted a Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the DEGs in KOBAS software [36,37].

4.5. WGCNA and Correlation Analyses of Anthocyanin-Related Genes

We built gene co-expression networks based on the DEGs using the weighted correlation network analysis (WGCNA) in the R package [38]. The DEGs from the HT- and RT-treated samples were used for co-expression analysis for each of the four time points. Total anthocyanin and cy-3-gala contents were examined for correlation with the modules and all the genes in each module. Significant anthocyanin-related modules were detected based on the highest correlation values with total anthocyanin as well as cy-3-gala contents (Figure S1A,B).

4.6. qRT-PCR Analysis

We performed qRT-PCR with a 20 µL sample in the iCycler iQ5 system (Bio-Rad, Hercules, CA, USA). The reaction mixture contained 10 µL SYBR Green PCR Master Mix (TransGen, Beijing, China), 1 µL respective primer, and 1 µL cDNA. Each qRT-PCR reaction was performed with three biological replicates. The transcripts were analyzed using the (Ct) 2–ΔΔCt method [39]. MdActin was used as an internal control. Table S1 presents all the qRT-PCR primers.

4.7. Transient Expression of MdLBD37 in Apple Fruit

The transient expression assay was conducted using the apple cultivar ‘Otome’. The CDS and another 393 bp fragment of MdLBD37 were inserted into the pRI101 and pTRV2 vectors [40]. The recombinant plasmids MdLBD37-pRI and MdLBD37-pTRV were transformed into Agrobacterium tumefaciens GV3101. Fruits of ‘Otome’ collected at 50 DAFB and 140 DAFB were injected with transformed A. tumefaciens containing MdLBD37-pTRV and MdLBD37-pRI, respectively. Five days after injection, the apple peel around the injection site was sampled for phenotypic analysis, anthocyanin content determination, and qRT-PCR assays. Ten successfully infected fruits were considered to determine anthocyanin content and gene expression.

4.8. Transformation of the Red-Fleshed Apple Calli with MdLBD37

The CDS of MdLBD37 was ligated into the SalI and SmaI sites of the pRI101 vector to generate the overexpression vector MdLBD37-pRI. The recombinant plasmids were transformed into Agrobacterium tumefaciens LBA4404. ‘Purple 3’ red-fleshed apple calli were used for transformation. Transgenic apple calli were selected on MS medium containing 50 mg L−1 kanamycin and 250 mg L−1 carbenicillin. Each successfully infected apple calli line is used as a biological repeat to determine anthocyanin content and gene expression.

5. Conclusions

In this study, we compared the total anthocyanin and cyanidin 3-galactoside contents and transcriptome in the peel of ‘Redchief Delicious’ apples treated at 25 °C and 35 °C. We identified one module that exhibited a significant correlation with the total anthocyanin and cy-3-gala contents using the WGCNA method. We characterized MdLBD37 in the module as an important negative regulator in anthocyanin accumulation at high temperatures. Our findings provide useful insights into the molecular mechanism of anthocyanin inhibition at high temperatures in apples, which lays the foundation for the genetic improvement of fruit coloration in apples under high temperatures.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23073766/s1.

Author Contributions

Y.-F.B. and S.-Q.F. designed the assays and wrote the manuscript. S.W., C.-Z.L., Y.-F.B., Y.F., Y.Z., Q.-Y.L., X.-S.C. and H.-B.W. performed assays and data analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (Grant No. 2018YFD1000105), the National Natural Science Foundation of China (Grant No. 31872940), and the Agricultural Improved Seed Project of Shandong Province (Grant No. 2019LZGC008, 2021LZGC024).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, N.; Jiang, S.H.; Zhang, Z.Y.; Fang, H.C.; Xu, H.F.; Wang, Y.C.; Chen, X. Malus sieversii: The origin, flavonoid synthesis mechanism, and breeding of red-skinned and red-fleshed apples. Hortic. Res. 2018, 5, 70. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Liao, J.J.; Zang, J.; Yuan, F.; Liu, S.; Zhang, Y.B.; Li, H.Y.; Piao, Z.; Li, H. Identification and analysis of anthocyanin components in fruit color variation in Schisandra chinensis. J. Sci. Food Agric. 2016, 96, 3213–3219. [Google Scholar] [CrossRef] [PubMed]
  3. Kayesh, E.; Shangguan, L.F.; Korir, N.K.; Sun, X.; Bilkish, N.; Zhang, Y.P.; Han, J.; Song, C.; Cheng, Z.-M.; Fang, J. Fruit skin color and the role of anthocyanin. Acta Physiol. Plant. 2013, 35, 2879–2890. [Google Scholar] [CrossRef]
  4. Liu, H.M.; Liu, Z.J.; Wu, Y.; Zhang, L.M.; Zhang, G.F. Regulatory Mechanisms of Anthocyanin Biosynthesis in Apple and Pear. Int. J. Mol. Sci. 2021, 22, 8441. [Google Scholar] [CrossRef]
  5. Holton, T.A.; Cornish, E.C.J.P.C. Genetics and biochemistry of anthocyanin biosynthesis. Plant Cell 1995, 7, 1071–1083. [Google Scholar] [CrossRef]
  6. Mori, K.; Sugaya, S.; Gemma, H.J.S.H. Decreased anthocyanin biosynthesis in grape berries grown under elevated night temperature condition. Sci. Hortic. 2005, 105, 319–330. [Google Scholar] [CrossRef]
  7. Xie, X.B.; Li, S.; Zhang, R.F.; Zhao, J.; Chen, Y.C.; Zhao, Q.; Yao, Y.X.; You, C.X.; Zhang, X.S.; Hao, Y.J. The bHLH transcription factor MdbHLH3 promotes anthocyanin accumulation and fruit colouration in response to low temperature in apples. Plant Cell Environ. 2012, 35, 1884. [Google Scholar] [CrossRef] [PubMed]
  8. Lin-Wang, K.; Micheletti, D.; Palmer, J.; Volz, R.; Lozano, L.; Espley, R.; Hellens, R.P.; Chagné, D.; Rowan, D.D.; Troggio, M.; et al. High temperature reduces apple fruit colour via modulation of the anthocyanin regulatory complex. Plant Cell Environ. 2011, 34, 1176–1190. [Google Scholar] [CrossRef]
  9. Fang, H.C.; Dong, Y.H.; Yue, X.X.; Chen, X.L.; He, N.B.; Hu, J.F.; Jiang, S.; Xu, H.; Wang, Y.; Su, M.; et al. MdCOL4 interaction mediates crosstalk between UV-B and high temperature to control fruit coloration in apple. Plant Cell Physiol. 2019, 60, 1055–1066. [Google Scholar] [CrossRef] [PubMed]
  10. Ban, Y.; Kondo, S.; Ubi, B.E.; Honda, C.; Bessho, H.; Moriguchi, T. UDP-sugar biosynthetic pathway: Contribution to cyanidin 3-galactoside biosynthesis in apple skin. Planta 2009, 230, 871–881. [Google Scholar] [CrossRef] [PubMed]
  11. Battisti, D.S.; Naylor, R.L. Historical warnings of future food insecurity with unprecedented seasonal heat. Science 2009, 323, 240–244. [Google Scholar]
  12. Lobell, D.B.; Schlenker, W.; Costa-Roberts, J. Climate trends and global crop production since 1980. Science 2011, 333, 616–620. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Shaked-Sachray, L.; Weiss, D.; Reuveni, M.; Nissim-Levi, A.; Oren-Shamir, M. Increased anthocyanin accumulation in aster flowers at elevated temperatures due to magnesium treatment. Physiol. Plant 2002, 114, 559–565. [Google Scholar] [CrossRef]
  14. Mori, K.; Goto-Yamamoto, N.; Kitayama, M.; Hashizume, K. Loss of anthocyanins in red-wine grape under high temperature. J. Exp. Bot. 2007, 58, 1935–1945. [Google Scholar] [CrossRef] [PubMed]
  15. Saure, M.C. External control of anthocyanin formation in apple. Sci. Hortic. 1990, 42, 182–218. [Google Scholar] [CrossRef]
  16. Lancaster, J.E.; Dougall, D.K. Regulation of skin color in apples. Plant Sci. 1992, 10, 487–502. [Google Scholar] [CrossRef]
  17. Borevitz, J.O.; Xia, Y.; Blount, J.; Dixon, R.A.; Lamb, C. Activation tagging identifies a conserved MYB regulator of phenylpropanoid biosynthesis. Plant Cell 2000, 12, 2383–2394. [Google Scholar]
  18. Ramsay, N.A.; Walker, A.R.; Mooney, M.; Gray, J.C. Two basic-helix-loop-helix genes (MYC-146 and GL3) from Arabidopsis can activate anthocyanin biosynthesis in a white-flowered Matthiola incana mutant. Plant Mol. Biol. 2003, 52, 679–688. [Google Scholar] [CrossRef] [PubMed]
  19. Morishita, T.; Kojima, Y.; Maruta, T.; Nishizawa-Yokoi, A.; Yabuta, Y.; Shigeoka, S. Arabidopsis NAC transcription factor, ANAC078, regulates flavonoid biosynthesis under high-light. Plant Cell Physiol. 2009, 50, 2210–2222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Li, H.H.; Liu, X.; An, J.P.; Hao, Y.J.; Wang, X.F.; You, C.X. Cloning and elucidation of the functional role of apple MdLBD13 in anthocyanin biosynthesis and nitrate assimilation. Plant. Cell Tissue Organ. Cult. 2017, 130, 47–59. [Google Scholar] [CrossRef]
  21. Walker, A.R.; Davison, P.A.; Bolognesi-Winfield, A.C.; James, C.M.; Srini-vasan, N.; Blundell, T.L.; Esch, J.J.; Marks, M.D.; Gray, J.C. The TRANSPARENT TESTA GLABRA1 locus, which regulates trichome differentiation and anthocyanin biosynthesis in Arabidopsis, encodes a WD40 repeat protein. Plant Cell 1999, 11, 1337–1350. [Google Scholar] [CrossRef] [Green Version]
  22. Takos, A.M.; Jaffé, F.W.; Jacob, S.R.; Bogs, J.; Robinson, S.P.; Walker, A.R. Light-induced expression of a MYB gene regulates anthocyanin biosynthesis in red apples. Plant Physiol. 2006, 142, 1216–1232. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Niu, J.P.; Zhang, G.J.; Zhang, W.T.; Goltsev, V.; Sun, S.; Wang, J.Z.; Li, P.M.; Ma, F.W. Anthocyanin concentration depends on the counterbalance between its synthesis and degradation in plum fruit at high temperature. Sci. Rep. 2017, 7, 7684. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. El-Sharkawy, I.; Liang, D.; Xu, K.N. Transcriptome analysis of an apple (Malus × Domestica) yellow fruit somatic mutation identifies a gene network module highly associated with anthocyanin and epigenetic regulation. J. Exp. Bot. 2015, 66, 7359–7376. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Liu, Y.H.; Lin-Wang, K.; Espley, R.V.; Wang, L.; Li, Y.M.; Liu, Z.; Zhou, P.; Zeng, L.; Zhang, X.; Zhang, J.; et al. StMYB44 negatively regulates anthocyanin biosynthesis at high temperatures in tuber flesh of potato. J. Exp. Bot. 2019, 70, 3809–3824. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Rubin, G.; Tohge, T.; Matsuda, F.; Saito, K.; Scheible, W.R. Members of the LBD family of transcription factors repress anthocyanin synthesis and affect additional nitrogen responses in Arabidopsis. Plant Cell 2009, 21, 3567–3584. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Koes, R.; Verweij, W.; Quattrocchio, F. Flavonoids: A colorful model for the regulation and evolution of biochemical pathways. Trends Plant Sci. 2005, 10, 236–242. [Google Scholar]
  28. Albert, N.W.; Davies, K.M.; Lewis, D.H.; Zhang, H.; Montefiori, M.; Brendolise, C.; Boase, M.R.; Ngo, H.; Jameson, P.E.; Schwinn, K.E. A conserved network of transcriptional activators and repressors regulates anthocyanin pigmentation in eudicots. Plant Cell 2014, 26, 962–980. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Zhang, J.; Xu, H.F.; Wang, N.; Jiang, S.H.; Fang, H.C.; Zhang, Z.; Yang, G.; Wang, Y.; Su, M.; Xu, L.; et al. The ethylene response factor MdERF1B regulates anthocyanin and proanthocyanidin biosynthesis in apple. Plant Mol. Biol. 2018, 98, 205–218. [Google Scholar] [CrossRef] [PubMed]
  30. Bi, S.Q.; An, J.P.; Wang, X.F.; Hao, Y.J.; You, C.X. Ethylene response factor MdERF3 promotes anthocyanin and proanthocyanidin accumulation in apple. Acta Hortic. Sin. 2019, 46, 2277–2285. [Google Scholar]
  31. An, J.P.; Song, L.Q.; Zhao, L.L.; You, C.X.; Wang, X.F.; Hao, Y.J. Overexpression of MdNAC029 promotes anthocyanin accumulation in apple calli. Acta Hortic. Sin. 2019, 46, 2277–2285. [Google Scholar]
  32. Giusti, M.M.; Wrolstad, R.E. Characterization and measurement of anthocyanins by UV-visible spectroscopy. In Current Protocols in Food Analytical Chemistry; Wiley: Hoboken, NJ, USA, 2001; pp. F1.2.1–F1.2.13. [Google Scholar]
  33. Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; Baren, M.J.V.; Salzberg, S.L.; Wold, B.J.; Pachter, L. Transcript assembly and quantification by RNA-seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 2010, 28, 511–515. [Google Scholar] [CrossRef] [Green Version]
  34. Wang, L.; Feng, Z.; Wang, X.; Wang, X.; Zhang, X. DEGseq: An R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics 2010, 26, 136–138. [Google Scholar] [CrossRef] [PubMed]
  35. Conesa, A.; Gotz, S.; Garcia-Gomez, J.M.; Terol, J.; Talon, M.; Robles, M. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21, 3674–3676. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Mao, X.Z.; Cai, T.; Olyarchuk, J.G.; Wei, L.P. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics 2005, 21, 3787–3793. [Google Scholar] [CrossRef] [PubMed]
  37. Ogata, H.; Goto, S.; Sato, K.; Fujibuchi, W.; Bono, H.; Kanehisa, M. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 1999, 27, 29–34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Langfelder, P.; Horvath, S. WGCNA: An R package for weighted correlation network analysis. BMC Bioinform. 2008, 9, 559. [Google Scholar] [CrossRef] [Green Version]
  39. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  40. Dinesh-Kumar, S.P.; Anandalakshmi, R.; Marathe, R.; Schiff, M.; Liu, Y.L. Virus-induced gene silencing. Methods Mol. Biol. 2003, 236, 287–294. [Google Scholar] [PubMed]
Ijms 23 03766 g001 550
Figure 1. Phenotypes and total anthocyanin and cy-3-gala contents of apples at high-temperature treatment. (A) The phenotypes of apples at T1, T2, T3, and T4. (B) Cy-3-gala contents of apples at T1, T2, T3, and T4. (C) The total anthocyanin contents of apples at T1, T2, T3, and T4. T1, T2, T3, and T4 represent apples after 3, 6, 9, and 12 days of 25 °C and 35 °C treatment, respectively. Values stand for means ± SD of three independent biological replicates. Statistical significance was measured using Student’s t-test (** p < 0.01).
Figure 1. Phenotypes and total anthocyanin and cy-3-gala contents of apples at high-temperature treatment. (A) The phenotypes of apples at T1, T2, T3, and T4. (B) Cy-3-gala contents of apples at T1, T2, T3, and T4. (C) The total anthocyanin contents of apples at T1, T2, T3, and T4. T1, T2, T3, and T4 represent apples after 3, 6, 9, and 12 days of 25 °C and 35 °C treatment, respectively. Values stand for means ± SD of three independent biological replicates. Statistical significance was measured using Student’s t-test (** p < 0.01).
Ijms 23 03766 g001
Ijms 23 03766 g002 550
Figure 2. Histogram and Venn diagram showing numbers of differentially expressed genes (DEGs) in the comparisons 35 °C vs. 25 °C at T1, T2, T3, and T4. (A) Histogram showing DEGs in the comparisons 35 °C vs. 25 °C at T1, T2, T3, and T4. (B) Venn diagram presenting the overlap of DEGs in the comparisons 35 °C vs. 25 °C at T1, T2, T3, and T4. (C) Expression levels analysis of 10 selected DEGs by qRT-PCR. T1, T2, T3, and T4 represent apples after 3, 6, 9, and 12 days of 25 °C and 35 °C treatment, respectively. Values stand for means ± SD of three independent biological replicates. Statistical significance was measured using Student’s t-test (** p < 0.01).
Figure 2. Histogram and Venn diagram showing numbers of differentially expressed genes (DEGs) in the comparisons 35 °C vs. 25 °C at T1, T2, T3, and T4. (A) Histogram showing DEGs in the comparisons 35 °C vs. 25 °C at T1, T2, T3, and T4. (B) Venn diagram presenting the overlap of DEGs in the comparisons 35 °C vs. 25 °C at T1, T2, T3, and T4. (C) Expression levels analysis of 10 selected DEGs by qRT-PCR. T1, T2, T3, and T4 represent apples after 3, 6, 9, and 12 days of 25 °C and 35 °C treatment, respectively. Values stand for means ± SD of three independent biological replicates. Statistical significance was measured using Student’s t-test (** p < 0.01).
Ijms 23 03766 g002
Ijms 23 03766 g003 550
Figure 3. Gene Ontology (GO) enrichment analysis of differentially expressed genes (DEGs) in the apple in the comparisons of 35 °C vs. 25 °C at T1, T2, T3, and T4. The DEGs were presented in three main types, respectively, ‘biological process’, ‘cellular component’, and ‘molecular function’. T1, T2, T3, and T4 represent apples after 3, 6, 9, and 12 days of 25 °C and 35 °C treatment, respectively.
Figure 3. Gene Ontology (GO) enrichment analysis of differentially expressed genes (DEGs) in the apple in the comparisons of 35 °C vs. 25 °C at T1, T2, T3, and T4. The DEGs were presented in three main types, respectively, ‘biological process’, ‘cellular component’, and ‘molecular function’. T1, T2, T3, and T4 represent apples after 3, 6, 9, and 12 days of 25 °C and 35 °C treatment, respectively.
Ijms 23 03766 g003
Ijms 23 03766 g004 550
Figure 4. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of differentially expressed genes (DEGs) in apple in the comparisons of 35 °C vs. 25 °C at T1, T2, T3, and T4. T1, T2, T3, and T4 represent apples after 3, 6, 9, and 12 days of 25 °C and 35 °C treatment, respectively.
Figure 4. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of differentially expressed genes (DEGs) in apple in the comparisons of 35 °C vs. 25 °C at T1, T2, T3, and T4. T1, T2, T3, and T4 represent apples after 3, 6, 9, and 12 days of 25 °C and 35 °C treatment, respectively.
Ijms 23 03766 g004
Ijms 23 03766 g005 550
Figure 5. Weighted gene co-expression network analysis (WGCNA) of DEGs in apple at 25 °C and 35 °C treatment. (A) Hierarchical cluster tree presenting modules of co-expressed genes with color annotation. (B) Module–anthocyanin correlations and corresponding p-values. The left panel shows eleven modules. The color scale on the right displays module-trait correlation from −1 (blue) to 1 (red). (C) Connection network between 12 anthocyanin structural genes and 21 TFs (transcript factors) in the ‘green’ module. The coefficients were calculated from log2 (fold change) of each transcript using the EXCEL program. This study selected correlations with a coefficient of R2 > 0.9. Positive correlations are presented with red lines, while negative correlations are presented with green lines. The color scale on the right shows anthocyanin structural gene-TF correlation numbers from 2 (green) to 14 (red). The correlation relationships were visualized using Cytoscape (version 2.8.2). (D) The heat map for the expression of 12 anthocyanin structural genes and 21 TFs (transcript factors) in the ‘green’ module. The numerical values for the green-to-red gradient bar indicate the log2-fold change relative to the control sample. The red pentagram indicates MdLBD37.
Figure 5. Weighted gene co-expression network analysis (WGCNA) of DEGs in apple at 25 °C and 35 °C treatment. (A) Hierarchical cluster tree presenting modules of co-expressed genes with color annotation. (B) Module–anthocyanin correlations and corresponding p-values. The left panel shows eleven modules. The color scale on the right displays module-trait correlation from −1 (blue) to 1 (red). (C) Connection network between 12 anthocyanin structural genes and 21 TFs (transcript factors) in the ‘green’ module. The coefficients were calculated from log2 (fold change) of each transcript using the EXCEL program. This study selected correlations with a coefficient of R2 > 0.9. Positive correlations are presented with red lines, while negative correlations are presented with green lines. The color scale on the right shows anthocyanin structural gene-TF correlation numbers from 2 (green) to 14 (red). The correlation relationships were visualized using Cytoscape (version 2.8.2). (D) The heat map for the expression of 12 anthocyanin structural genes and 21 TFs (transcript factors) in the ‘green’ module. The numerical values for the green-to-red gradient bar indicate the log2-fold change relative to the control sample. The red pentagram indicates MdLBD37.
Ijms 23 03766 g005
Ijms 23 03766 g006 550
Figure 6. MdLBD37-mediated inhibition of anthocyanin synthesis in apple. (A) The phenotypes of WT, MdLBD37-TRV, and MdLBD37-pRI apple fruit. (B) The expression levels of MdLBD37 in WT, MdLBD37-TRV, and MdLBD37-pRI apple fruit. (C) Anthocyanin contents in WT, MdLBD37-TRV, and MdLBD37-pRI apple fruit. FW, fresh weight. (D) The expression levels of anthocyanin-related genes in WT, MdLBD37-TRV, and MdLBD37-pRI apple fruit. (E) Phenotypes of MdLBD37-overexpressing (MdLBD37-OE) apple calli. (F) The expression levels of MdLBD37 in WT and MdLBD37-OE apple calli. (G) Anthocyanin contents in WT and MdLBD37-OE apple calli. (H) The expression levels of anthocyanin-related genes in WT and MdLBD37-OE apple calli. FW, fresh weight. Values suggest means ± SD of three independent biological replicates. Statistical significance was computed with the use of Student’s t-test (** p < 0.01).
Figure 6. MdLBD37-mediated inhibition of anthocyanin synthesis in apple. (A) The phenotypes of WT, MdLBD37-TRV, and MdLBD37-pRI apple fruit. (B) The expression levels of MdLBD37 in WT, MdLBD37-TRV, and MdLBD37-pRI apple fruit. (C) Anthocyanin contents in WT, MdLBD37-TRV, and MdLBD37-pRI apple fruit. FW, fresh weight. (D) The expression levels of anthocyanin-related genes in WT, MdLBD37-TRV, and MdLBD37-pRI apple fruit. (E) Phenotypes of MdLBD37-overexpressing (MdLBD37-OE) apple calli. (F) The expression levels of MdLBD37 in WT and MdLBD37-OE apple calli. (G) Anthocyanin contents in WT and MdLBD37-OE apple calli. (H) The expression levels of anthocyanin-related genes in WT and MdLBD37-OE apple calli. FW, fresh weight. Values suggest means ± SD of three independent biological replicates. Statistical significance was computed with the use of Student’s t-test (** p < 0.01).
Ijms 23 03766 g006
Ijms 23 03766 g007 550
Figure 7. A hypothetical model depicting MdLBD37 and other TFs from the WGCNA module ‘green’ might be involved into the inhibition of high temperature on anthocyanin synthesis. The dotted line indicated possible links.
Figure 7. A hypothetical model depicting MdLBD37 and other TFs from the WGCNA module ‘green’ might be involved into the inhibition of high temperature on anthocyanin synthesis. The dotted line indicated possible links.
Ijms 23 03766 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bu, Y.-F.; Wang, S.; Li, C.-Z.; Fang, Y.; Zhang, Y.; Li, Q.-Y.; Wang, H.-B.; Chen, X.-S.; Feng, S.-Q. Transcriptome Analysis of Apples in High-Temperature Treatments Reveals a Role of MdLBD37 in the Inhibition of Anthocyanin Accumulation. Int. J. Mol. Sci. 2022, 23, 3766. https://doi.org/10.3390/ijms23073766

AMA Style

Bu Y-F, Wang S, Li C-Z, Fang Y, Zhang Y, Li Q-Y, Wang H-B, Chen X-S, Feng S-Q. Transcriptome Analysis of Apples in High-Temperature Treatments Reveals a Role of MdLBD37 in the Inhibition of Anthocyanin Accumulation. International Journal of Molecular Sciences. 2022; 23(7):3766. https://doi.org/10.3390/ijms23073766

Chicago/Turabian Style

Bu, Yu-Feng, Shuo Wang, Chen-Zhiyu Li, Yue Fang, Ya Zhang, Qing-Yu Li, Hai-Bo Wang, Xue-Sen Chen, and Shou-Qian Feng. 2022. "Transcriptome Analysis of Apples in High-Temperature Treatments Reveals a Role of MdLBD37 in the Inhibition of Anthocyanin Accumulation" International Journal of Molecular Sciences 23, no. 7: 3766. https://doi.org/10.3390/ijms23073766

APA Style

Bu, Y. -F., Wang, S., Li, C. -Z., Fang, Y., Zhang, Y., Li, Q. -Y., Wang, H. -B., Chen, X. -S., & Feng, S. -Q. (2022). Transcriptome Analysis of Apples in High-Temperature Treatments Reveals a Role of MdLBD37 in the Inhibition of Anthocyanin Accumulation. International Journal of Molecular Sciences, 23(7), 3766. https://doi.org/10.3390/ijms23073766

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop