- 1College of Horticulture and Plant Protection, Inner Mongolia Agricultural University, Hohhot, China
- 2Beijing Engineering Research Center for Edible Mushroom, Institute of Plant Protection, Beijing Academy of Agriculture and Forestry Sciences, Beijing, China
Introduction: Valsa canker, caused by Cytospora mali, is a destructive disease in apple production. However, the mechanism by which apple defend against C. mali infection remains unclear.
Methods: In this study, the integrative transcriptional and metabolic analysis were used to investigate the responses of the ‘Jin Hong’ apple branches to the invasion of C. mali.
Results and Discussion: Results showed that the differentially expressed genes were mainly enriched in the pathways of carbon metabolism, photosynthesis-antenna proteins, and biosynthesis of amino acids pathways. Additionally, the differentially accumulated metabolites were significantly enriched in aminoacyl-tRNA biosynthesis, fructose and mannose metabolism, and alanine, aspartate, and glutamate metabolism pathways. Conjoint analysis revealed that C. mali infection significantly altered 5 metabolic pathways, 8 highly relevant metabolites and 15 genes of apples. Among which the transcription factors WRKY and basic domain leucine zipper transcription family were induced, the α-linolenic acid and betaine were significantly accumulated in C. mali infected apple stems. This work presents an overview of the changes in gene expression and metabolic profiles in apple under the inoculation of C. mali, which may help to further screen out the mechanism of plant-pathogen interaction at the molecular level.
1 Introduction
Apple Valsa Canker (AVC), also known as Cytospora canker, is caused by the necrotrophic fungus Cytospora mali Grove (Yin et al., 2015; Wang D. L. et al., 2020). AVC stands out as one of the most destructive diseases affecting apples (Malus domestica) in the world, leading to significant yield losses and decreased fruit quality. Extensive surveys have revealed that the incidence of AVC averaged approximately 53% in the major apple production regions of China, and reached up to 80% in severely affected orchards (Meng et al., 2019).
AVC mainly occurs on the branches of trees, the infected tissues then gradually dry out, collapse slightly, and finally form localized cankers and even kills entire apple trees (Feng et al., 2023). C. mali can extensively penetrate into the host’s phloem and xylem, which is difficult for traditional chemical agents to access, and making it harder to prevent (Yin et al., 2015).
The molecular mechanism of apples in responses to C. mali infection has been explored in recent years. It is reported that apples could resistant against AVC by negatively regulating the production of phloridzin (Zhou et al., 2019, 2021). Mao et al. (2021) revealed that MdCN11 and MdCN19 genes could negatively regulate the AVC resistance via inducting hypersensitive response. DEGs related to hormonal, Ca2+ signaling, and phenylpropanoid biosynthesis were highly expressed to enhance the apple’s resistance to C. mali (Yin et al., 2016; Zuo et al., 2018). Wen et al. (2023) established an efficient system for screening disease-resistant genes in response to AVC. In addition, the apple activates and accumulates defense-related enzymes and metabolites to protect itself from C. mali, such as phenylalanine ammonia-lyase (PAL), β-1,3-glucanase, and chitinase, dopamine, trans-cinnamic acid, coumarin as well as chlorogenic acid and other phenylpropanoids and flavonoids (Sung and Lee, 2010; Yu et al., 2016; Moneo-Sánchez et al., 2018; Geng et al., 2020; Liu X. M. et al., 2021; Liu X. Z. et al., 2021; Zhao et al., 2021; Du et al., 2023; Li C. R. et al., 2023). In recent years, the metabolomics and transcriptomics have become a research hotspot to explore the molecular mechanisms of disease resistance in apples against pathogens. Zhao et al. (2022) revealed mechanisms of disease resistance in apples induced by Wickerhamomyces anomalus, using the integrated transcriptomic and metabonomic analysis. The transcriptomics was used to explore the molecular mechanisms of apples in responses to C. mali infection (Feng et al., 2023). However, there are no studies on combined transcriptomics and metabolomics analyses of apples by challenging with C. mali, and some highly relevant genes and metabolites remain unexplored.
This work conducted an integrative analysis of the transcriptome and metabolome of apple branches challenging with C. mali. Significant changes were observed in both the transcriptional and metabolic profiles. A conjoint analysis was then performed to identify co-enrichment pathways in both transcriptional and metabolic data. A number of altered pathways, as well as highly correlated genes and metabolites, were identified. These findings will deepen our understanding of AVC caused by C. mali, and provide valuable insights for further investigation into the molecular mechanisms of apples against AVC.
2 Materials and methods
2.1 Plant materials and fungal strains
C. mali (strain QH2), one of the main species causing AVC in northern China (Ma et al., 2020), was used in this study. The strain was isolated from the infected stems of apple tree in Hohhot, Inner Mongolia. The apple branches used in this study were obtained from the orchard in the Inner Mongolia Agricultural University’s campus, the collected branches were inoculated with C. mali QH2 following Zang’s method (Zang et al., 2007). Briefly, branches from “Jin Hong” apple tree were cut into 30 cm length, immersed in 1% sodium hypochlorite for 10 min, 75% alcohol for 3 min, and washed with sterile water 3 times. The two ends of the twigs were sealed by dipping in heated liquid wax and the middle of the twigs was scalded by a soldering iron 6 mm in diameter. C. mali QH2 disks (Φ = 6 mm) were placed on the scalded area. The fungal disks were wrapped with sterile wet degreasing cotton and preservative film. Healthy branches inoculated with agar were treated as control. Treated branches were placed in a pallet with two players of gauze, and incubated in an illumination incubator with a 16 h light/8 h dark cycle, 85% humidity at 25°C for 9 days (Figure 1). Approximately 2.0 g of bark tissue, including buds, from the borders of infected area was sampled. Each sample was divided into two parts, for transcriptomic and metabolomic analysis, respectively.
2.2 Total RNA extraction, library construction, and sequencing
The total RNA was extracted using the Takara Mini BEST plant RNA Extraction Kit (Takara, Kyoto, Japan) followed by the instructions. The RNA quality was assessed using the NanoPhotometer® spectrophotometer (IMPLEN, LA, CA, United States) and the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, Palo alto, CA, United States).
A total of 1 μg RNA per sample was used for RNA sample preparation. The NEBNext® UltraTM RNA Library Prep Kit for Illumina® (NEB, Ipswich, MA, USA) was used to generate sequencing libraries. Index codes were added to attribute sequences to each sample.
Briefly, mRNA was enriched by poly A tail selection and chemically fragmented. Subsequently, second strand cDNA synthesis was performed using DNA Polymerase I and RNase H. The synthesized cDNA was then subjected to end-repair and “A” base addition according to Illumina’s library construction protocol. To select cDNA fragments for PCR, we purified the library fragments using AMPure XP beads (Beckman Coulter, Beverly, United States) to ensure a length of 250–300 bp. We also purified the PCR products and assessed the library quality on the Agilent Bioanalyzer 2100 system. Following cluster generation, we sequenced the library preparations on an Illumina Novaseq platform, generating 150 bp paired-end reads.
2.3 Data quality control
We performed splicing and quality control on the raw data. Aligned paired-end clean reads to the reference genome using HISAT2 v2.0.5. Used Feature Counts (1.5.0-p3) to calculate the number of reads mapped to each gene. Then calculated the fragments per kilobase of transcript per million mapped reads (FPKM) for each gene based on gene length and determine the number of reads mapped to that gene.
2.4 Differential expression analysis
Differentially expressed genes (DEGs) were identified using the DESeq2 R package (version 1.16.1) (Varet et al., 2017). Differential expression analysis of two conditions was performed using the edgeR R package (version 3.18.1) (Dietz et al., 2010; McCarthy et al., 2012; Varet et al., 2017).
The resulting p were adjusted using the Benjamini and Hochberg’s approach (Yoav and Yosef, 2000) to control for false discovery rate. DEGs were determined according to the following general rules: |Log2FoldChange(FC)| > 1 and p-adjust ≤0.05. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway functional enrichment analyses were performed using the clusterProfiler R package (version 3.4.4) (Wu et al., 2021).
2.5 Metabolite extraction
The metabolite was extracted using a solution of methanol/acetonitrile/water in a 2: 2: 1 (v/v) ratio. The mixture was vortexed and subjected to low-temperature sonication for 30 min and stand at −20°C for 10 min, then centrifuged (Eppendorf, 5430R, Hamburg, Germany) at 14,000 g at 4°C for 20 min, the supernatant was collected and vacuum dried. For mass spectrometry analysis, 100 μL of acetonitrile/water solution (acetonitrile: water = 1: 1, v/v) was added for reconstitution, followed by vertexing, the reconstituted sample was then centrifuged at 14,000 g for 15 min at 4°C to obtain the supernatant.
2.6 Non-target metabolomics
Metabolites were chromatographically separated using an ultra-high performance liquid chromatography (UHPLC) system (Agilent 1290 Infinity LC, Agilent Technologies, Palo Alto, California, United States). The UHPLC system was equipped with an ACQUITY BEH C18 column (100 mm × 2.1 mm column, 1.7 μm, Waters, Milford, MA, United States). The mobile phases consisted of solvent A (Water +25 mM ammonium acetate +25 mM ammonia) and solvent B (acetonitrile). The gradient elution program was set as follows to equilibrate the systems: from 0 to 0.5 min, 95% (B); from 0.5 to 7 min, 95 to 65% (B); from 7 to 8 min, 65 to 40% (B); from 8 to 9 min, 40% (B); 9 to 9.1 min, 40 to 95% (B); from 9.1 to 12 min, 95% (B). The sample injection volume was 2 μL, the flow rate was 0.5 mL min−1, and the column temperature was maintained at 25°C throughout the chromatographic separation.
Mass spectra of the compounds were obtained using an AB Triple TOF 6600 (AB SCIEX, Framingham, United States). ESI source conditions after HILIC chromatographic separation Ion source Gas1 (Gas1): 60, Ion source Gas2 (Gas2): 60, curtain gas (CUR): 30, source temperature: 600°C Ion source voltage floating (ISVF) ± 5,500 V; TOF MS scan m/z range: 60–1,000 Da, product ion scan m/z range: 25–1,000 Da, TOF MS scan accumulation time 0.20 s/spectra, product ion scan accumulation time 0.05 s/spectra; MS2 data were acquired using information dependent acquisition (IDA) and operated in high sensitivity mode, declustering potential (DP): ± 60 V, collision energy: 35 ± 15 eV, exclude isotopes within 4 Da, candidate ions to monitor per cycle: 10. Samples were analyzed in both positive (+) and negative (−) ESI modes.
Raw data in Wiff format were converted to mzXML format using ProteoWizard (Darren et al., 2008) and then peaks were aligned and peak areas extracted using XCMS software (Domingo-Almenara and Siuzdak, 2020). The data extracted by XCMS were used for metabolite structure identification and data pre-processing, quality evaluation, and multi-variate analysis.
2.7 Differentially accumulated metabolites analysis
Orthogonal partial least squares discriminant analysis (OPLS–DA) was used to determine the differentially accumulated metabolites between pairwise groups. The p was estimated using a Hotelling’s T2 test for statistical analysis. Differentially accumulated metabolites (DAMs) (VIP ≥ 1, |Log2FC| > 1) between groups were mapped into biochemical pathways using metabolic enrichment and pathway analysis based on MetaboAnalyst.1 Those pathways with p ≤ 0.05 were considered significantly enriched.
2.8 Statistical analysis
OPLS–DA was performed using SIMCA 14.1 (Umetrics, Umea, Sweden). The data were log-transformed. The reliability of the predictive models was assessed using the analysis of variance test of cross-validated predictive residuals (CV–ANOVA), the coefficient of determination (R2) and the predicted variance (Q2) which is the proportion of the total variation of X or Y that can be predicted by a component. To test the overall associations between DEGs and DAMs, the DEGs and DAMs co-enriched in KEGG pathways were identified by KEGG analysis. The DEGs and DAMs were subjected to the integration analysis using regularized canonical correlation analysis (rCCA) (Cruz-Cano and Lee, 2014) to explore their relationships. rCCA was implemented with mixOmics (Rohart et al., 2017) using the R program. The regularization parameters (lambda1, lambda2) were optimized using the tune.rcc function. The network communities were determined using the Fast Greedy algorithms with the clusterMaker program (Morris et al., 2011; Utriainen and Morris, 2023) and visualized using Cytoscape 3.10.1 (Zhou et al., 2022). Network attributes were calculated using Cytoscape’s built-in network analyzer. Bioinformatics analysis was performed using dynamic real-time interactive online platforms including APPLIED PROTEIN TECHNOLOGY,2 Metware Cloud,3 Majorbio,4 Omicsmart,5 and bioinformatics.6
3 Results
3.1 Transcriptional profiling reveals altered pathways and genes in diseased apple branches
The transcriptome datasets generated in this work are deposited in the NCBI repository with the accession number: PRJNA1074127, following the link: https://www.ncbi.nlm.nih.gov/sra/PRJNA1074127. The PCA analysis revealed that the two principal components explained 60.9 and 15.7% of the variation and effectively separated samples from different stages (Figure 2A). The Venn analysis showed that 1,271 genes were unique to the control group, 1,107 genes were unique to AVC, and 37,511 genes were shared between the two groups (Figure 2B). The genes related to cytochrome P450, protein kinase domain, and homeobox domain-like were the most prevalent among control-specific genes. Among AVCs, the highest percentage of genes were for protein kinase domain, leucine-rich repeat, and cytochrome P450 (Figure 2C). The most predominant genes shared by both groups were for protein kinase domain, serine–threonine/tyrosine-protein kinase, and catalytic domain (Figure 2D). The results indicated that the expression of 5, 413 DEGs was up-regulated and that of 5, 298 DEGs was down-regulated in the AVC group compared to the control group (Figure 2E). The DEGs identified were annotated using KEGG enrichment analysis. The up-and down-regulated DEGs were enriched in different pathways. The up-regulated genes were mostly enriched in the proteasome, biosynthesis of amino acids and oxidative phosphorylation. On the other hand, the down-regulated genes were enriched in photosynthesis—antenna proteins, plant hormone signal transduction and photosynthesis (Figures 2F,G). The KEGG analysis of all DEGs revealed that the pathways were significantly involved in carbon metabolism, photosynthesis-antenna proteins, and biosynthesis of amino acids (Supplementary Table S1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1394447/fmicb-15-1394447-HTML/image_m/fmicb-15-1394447-g002.jpg)
Figure 2. Gene expression analysis of “Jin Hong” in response to AVC infection. (A) Principal component analysis (PCA) of transcriptome data. (B) The venn plot of differentially expressed genes. (C) The unique genes of two groups. (D) The genes shared by both groups. (E) The volcano plot showing the up-and down-regulated genes (p-adjust < 0.05, |Log2FC| > 1, red represents up-regulated genes, blue represents down-regulated genes, and gray indicates not significantly expressed genes, text highlighted for up-and down-regulated Log2FC top 5 genes). (F,G) The KEGG enriched analysis of the up-and down-regulated genes in AVC _vs_ Control (p-adjust < 0.05).
After performing GO analysis on the up-regulated genes, the most significantly enriched molecular function (MF) were threonine-type endopeptidase activity and threonine-type peptidase activity. Additionally, the most significantly enriched cellular component (CC) was organelle membrane. Moreover, the down-regulated mRNAs indicate significant enrichment in thylakoid, thylakoid part, and photosynthetic membrane (Supplementary Table S2).
3.2 Metabolic profiling reveals altered pathways and metabolites in diseased apple branches
A summary table highlighting the expressed metabolites can be found in Supplementary Table S3. The OPLS-DA analysis showed a significant difference between the control and AVC groups (Figure 3A). The identified metabolites were annotated using the HMDB database and classified into the following superclass categories: 16 phenylpropanoids and polyketides, 1 organooxygen compound, 16 organoheterocyclic compounds, 43 organic oxygen compounds, 5 organic nitrogen compounds, 40 organic acids and derivatives, 13 nucleosides, nucleotides, and analogs, 1 organonitrogen compound, 17 lipids and lipid-like molecules, 2 lignans, neolignans and related compounds, 11 benzenoids, and 20 metabolites of unknown classification. The top 30 metabolites, as shown in the loading plot (VIP > 1), were classified into 9 superclass categories (Figure 3B). The bar charts illustrated the significant regression coefficients of each compound for the different groups (Figure 3C). The number of compounds with significant correlations differed between the two groups. The metabolites of benzoic acid, L-threonate and L-gulonic gamma-lactone were positively correlated with the AVC group, while the metabolites of dulcitol, tyramine, and l-phenylalanine showed a negative correlation. Heatmap showed different superclasses of identified metabolites (Figure 3D). The volcano plot illustrated the DAMs between AVC and control. There were 2,433 up-regulated and 2,977 down-regulated DAMs (Figure 3E). The DAMs between the two groups were mapped to the KEGG database. The identified 48 DAMs were mainly enriched in pathways of aminoacyl-tRNA biosynthesis, fructose and mannose metabolism, alanine, aspartate and glutamate metabolism, galactose metabolism, and isoquinoline alkaloid biosynthesis (Figure 3F; Supplementary Table S4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1394447/fmicb-15-1394447-HTML/image_m/fmicb-15-1394447-g003.jpg)
Figure 3. Metabolic profiling of “Jin Hong” in response to AVC infection. (A) OPLS–DA analysis of the metabolome data of all samples (R2Xcum = 0.764, R2Ycum = 0.980, Q2cum = 0.969, CV–ANOVA, p = 5.577 × 10−3). (B) OPLS–DA loading showing the top 30 metabolites (VIP > 1), color indicates the superclass. (C) Correlation coefficient plots of OPLS–DA showing the relationship between the X and Y variables for the predictive components. The size of the scaled coefficient represents the change in the Y variable when the X variable varies from 0 to 1 in coded units (Selection of the top 5 substances with positive and negative correlation coefficients). (D) The heatmap shows the qualitative names of metabolites, the green bars represent the VIP values of the metabolites and the red dashed line indicates VIP = 1, different colored lines indicate the superclass classification of metabolites. (E) The volcano plot of DAMs in ‘Jin Hong’ upon exposure to AVC (p < 0.05, |Log2FC| > 1, red represents up-regulated metabolites, blue represents down-regulated metabolites, and gray indicates not significantly expressed genes, text highlighted for up-and down-regulated|Log2FC|top 5 metabolites). (F) The top 10 enriched KEGG pathways of DAMs in “Jin Hong” upon exposure to AVC infection (p < 0.05).
3.3 Comparative transcriptional and metabolic profiling revealed co-enrichment of pathways
Venn analysis results indicated that DEGs and DAMs were both enriched in 45 KEGG metabolic pathways, among which, 5 were significantly enriched (p < 0.05) (Figure 4A). Figure 4B illustrated pathways associated with plant stress tolerance, including 4 of the 45 co-enriched KEGG pathways and 8 pathways specific to DEGs. The DEGs and DAMs were significantly enriched in 5 KEGG pathways, including galactose metabolism, fructose and mannose metabolism, phenylalanine, tyrosine and tryptophan biosynthesis, glycine, serine and threonine metabolism, and isoquinoline alkaloid biosynthesis. Some metabolites are associated with plant resistance to adversity stress, and the metabolites in the 5 significantly enriched shared pathways contain sorbitol, D-mannose, tyramine, L-tryptophan, and betaine (Figures 4C,D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1394447/fmicb-15-1394447-HTML/image_m/fmicb-15-1394447-g004.jpg)
Figure 4. Transcriptome-metabolome-wide association of KEGG pathways. (A) The venn plot of DEGs and DAMs enriched all KEGG pathways and significant KEGG pathways. (B) The pathways associated with plant stress resistance among the 45 commonly enriched pathways as well as those unique to DEGs. (C,D) The heatmap of metabolites and genes among the 5 commonly significant enriched pathways, the triangles are labeled with substances associated with plant stress tolerance and the different colored columns indicate KEEG pathways.
3.4 Combined transcriptional and metabolic profiling revealed relevant genes and metabolites
Strong associations were found between genes and metabolites by procrustes analysis (Figure 5A). To identify genes and metabolites that are highly correlated with changes in apples in response to infestation with C. mali, the 45 DEGs and DAMs co-enriched pathways were subjected to a correlation analysis using the rCCA method (threshold > 0.90). The correlation network revealed 3 clear gene-metabolite regulatory communities (Figure 5B). In community 1, 94 transcripts exhibited a strong correlation coefficient with L-glutamine and 1-aminocyclopropanecarboxylic acid. Based on the magnitude of the betweenness centrality value, these metabolites were found to be highly correlated with genes MD01G1037400, MD05G1023900, MD14G1102200, MD15G1077000, and MD15G1204200. The interaction networks between the 15 metabolites and 79 transcripts were organized in control and AVC in community 2. The metabolites D-glucuronolactone and alpha-linolenic acid (ALA) had higher connectivity with genes MD00G1112500, MD02G1123600, MD04G1246300, MD10G1063600, and MD00G1061700. Community 3 was the largest gene-metabolite association, consisting of 13 metabolites and 204 genes. The genes MD16G1179500, MD17G1068900, MD07G1159300, MD12G1101800, and MD09G1078700 were found to be more correlated with the metabolites 2-oxoadipate and biotin (Supplementary Table S5). Additionally, we constructed a correlation network of DEGs and DAMs in 5 co-enriched pathways on the Pearson correlation coefficient (PCC) analysis results. Figure 5C showed that 15 DEGs and 8 DAMs in the network were highly correlated (PCC > 0.99). These genes and metabolites may play a crucial role in the defense of apples against AVC.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1394447/fmicb-15-1394447-HTML/image_m/fmicb-15-1394447-g005.jpg)
Figure 5. Transcriptome-metabolome-wide association network of all samples. (A) Procrustes analysis between all genes and metabolites from all samples. (B) Integration network of highly association genes and metabolites at threshold of 0.9. (C) Correlation analysis of genes and metabolites in 5 commonly significant pathways (threshold>0.99), blue ovals indicate metabolites, pink squares indicate genes, and red line segments indicate positive correlations and blue line segments indicate negative correlations.
4 Discussion
Almost all apple trees can be infected by AVC (Ke et al., 2014; Haxim et al., 2022). Plants have developed various mechanisms to cope with pathogenic stresses, such as forming physical barriers and activating their immune system, which mainly includes reinforcing the cell wall, producing reactive oxygen species, altering hormone levels and gene expression, activating pathogenesis-related proteins, accumulating defensive metabolites, and inducing the allergic necrosis of plant cells (Yin et al., 2015; Wang X. L. et al., 2020; Liu X. M. et al., 2021; Liu X. Z. et al., 2021; Wang et al., 2022; Xin et al., 2023). In this study, the RNA-seq and UHPLC–MS/MS techniques were utilized to investigate the changes in genes and metabolites in apple branches during C. mali infection, 5 metabolic pathways were significantly enriched: galactose metabolism, fructose and mannose metabolism, phenylalanine, tyrosine and tryptophan biosynthesis, glycine, serine and threonine metabolism, and isoquinoline alkaloid biosynthesis. Furthermore, the 15 key genes and 8 DAMs (tyramine, L-tyrosine, L-phenylalanine, sorbitol, L-tryptophan, L-threonine, melibiose, and betaine) that may be associated with plant’s reaction to pathogen stress were identified.
The cell wall provides initial protection and signal perception against pathogen attacks through sensing and regulating related components on the cell wall (Sebastian, 2022; Xin et al., 2023). Fungal pathogens have developed various combinations of plant cell wall degrading enzymes to break them down (Yin et al., 2015; Xin et al., 2023). β-galactosidase (β-Gal) is an essential cell wall glycosyl hydrolase in apple (Yin et al., 2015; Li X. L. et al., 2023). The gene MD02G1287500 encoding β-Gal was found to be down-regulated in the galactose metabolite pathway in this work, this result is consistent with previous results, that apples may slow down the degradation of the cell wall by decreasing the expression of β-Gal to resist pathogen invasion (Moneo-Sánchez et al., 2018; Zhao et al., 2021; Feng et al., 2023).
UDP-glucose catalyzes the glycosylation of phloretin to phloridzin, the pathogen could utilize phloridzin to produce toxins that facilitate necrosis in apple bark (Zhou et al., 2021). Our findings showed a 4-fold down-regulation of the UDP-glucose gene MD13G1093700 in the galactose metabolism pathway and decreased accumulation of its metabolite phenothiazine, although this change was not statistically significant. So this result is consistent with the previous finding that the apples could enhance AVC resistance by decreasing the phloridzin biosynthesis, and the down expression of UDP-glucose genes, which may indirectly modulating cell wall deposition and increasing hormone levels (Zhou et al., 2019, 2021; Feng et al., 2023).
Once the pathogen has penetrated the external cell wall, the intracellular immune system is triggered. Several reports have demonstrated that the WRKY and basic domain leucine zipper (bZIP) gene families play significant roles in various biological processes, including physiological metabolism and stress response (Alves et al., 2013; Jiang et al., 2017), and are particularly involved in resisting C. mali in apples (Tsuda and Somssich, 2015; Geng et al., 2020; Liu X. M. et al., 2021; Liu X. Z. et al., 2021; Mao et al., 2021; Wen et al., 2023). Moreover, WRKYs played a crucial role in mediating the communication between jasmonic acid and salicylic acid and are integral to the plant’s defense responses (Phukan et al., 2016; Lui et al., 2017; Han et al., 2023). Transcriptomic analysis in this study revealed that the invasion of C. mali activated the expression of WRKY and bZIP genes in apple. There were significant increase in the expression of the WRKY genes (MD13G1122100, MD09G1008800, MD17G113810, and MD03G119760) and bZIP gene (MD08G1025800) in the diseased apple branches. These high expression in WRKY and bZIP genes might associated with the apple’s defense response.
The destruction of cellular membranes in pathogen-affected plants releases a large number of polyunsaturated fatty acids, which act as signal molecules to activate plant defense responses (De Carvalho and Caramujo, 2018). For example, α-linolenic acid (ALA) serves as a precursor to jasmonic acid, a key phytohormone and signaling molecule that mediates plant stress response (Per et al., 2018; Bizuneh, 2020; Muñoz-Hoyos and Stam, 2023; Xin et al., 2023). The KEGG analyses of DEGs in this study revealed a significant up-regulation of fatty acid metabolic pathways, which is consistent with previous reports (Ke et al., 2014; Yin et al., 2016), and cluster 2 of the network showed a strong association between the up-regulated metabolite ALA and several genes. All these results underscore the crucial roles of the metabolite ALA in apple resistance to AVC.
Betaine could enhance the stability of biomacromolecular structure and function, and reduce ROS accumulation by preserving the activity of ROS scavenging enzymes (Huang et al., 2020), which leads to improved plant resistance against various abiotic stresses such as metal ions, salt, low temperature, and drought (Li C. Y. et al., 2021; Jokinen et al., 2022). Moreover, exogenous betaine was added to enhance apple resistance against anthracnose leaf blight (Liu et al., 2022). In our study, we observed a significant accumulation of the betaine in the Glycine, serine, and threonine metabolism pathway. Combining our findings with previous research, we hypothesize that betaine may play a crucial role in apple resistance to C. mali invasion.
The phenylpropanoid pathway leads to synthesis of different types of flavonoid phytoalexins and phenolic compounds, which are involved in plant defense (Li P. Q. et al., 2021). The pathway starts with phenylalanine, which can be converted into aromatic compounds and phenylpropanoids in response to disease resistance (Geng et al., 2020). In addition, phenylalanine ammonia-lyase (PAL) is a crucial enzyme in secondary phenylalanine metabolism and is widely studied in plant responses to biotic stresses (Kurihara and Yamana, 2022; Mansoor et al., 2023). We observed that many DEGs and DAMs were enriched in the phenylalanine, tyrosine, and tryptophan biosynthesis pathways in both groups. Within these pathways, the expression of PAL gene MD11G1223500 was significantly increased after being infected. These metabolites and genes in the phenylpropanoid pathway were associated with the resistance of apple to C. mali infestation, and need further exploration.
The infection of apples by C. mali disrupts the cell structure and secretes pathogenic virulence factors. In response to C. mali infection, apples have evolved multiple strategies. Our study indicates that apples primarily resist C. mali infection by mitigating cell wall degradation, activating the PAL activity, and promoting the synthesis of defensive substances.
5 Conclusion
This study conducted a comparative analysis of the metabolome and transcriptome of ‘Jinhong’ apple branches between control and AVC inoculation. Pathway analysis of the apple resistance response revealed significant changes in several metabolic pathways, including galactose metabolism, fructose and mannose metabolism, phenylalanine, tyrosine and tryptophan biosynthesis, glycine, serine and threonine metabolism, and isoquinoline alkaloid biosynthesis, during the C. mali inoculation. Moreover, the transcriptional and metabolic patterns of the affected tissues were altered. The production and accumulation of several enzymes and metabolites, including β-Gal, UGP-glucose, ALA, betaine, and PAL may play a crucial role in the presentation of apples when challenged by C. mali. These results shed light on the potential molecular mechanisms underlying the infection of AVC in ‘Jinhong’ branches, further exploration of DEGs and DEMs are in need in future work.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Author contributions
JZ: Formal analysis, Investigation, Writing – original draft, Methodology. YG: Data curation, Software, Writing – original draft. ZL: Resources, Writing – review & editing, Project administration. YS: Methodology, Writing – original draft. PS: Funding acquisition, Methodology, Project administration, Supervision, Validation, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Natural Science Foundation of China (Grant numbers: 32060594 and 31860540); Natural Science Foundation of Inner Mongolia, China (2023LHMS03020); and Start-up program of innovation and entrepreneurship for returned overseas Chinese scholars in the Inner Mongolia autonomous region (DC2100001765).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1394447/full#supplementary-material
Footnotes
References
Alves, M. S., Dadalto, S. P., Goncalves, A. B., Barros, V. A., and Fietto, L. G. (2013). Plant bZIP transcription factors responsive to pathogens: a review. Int. J. Mol. Sci. 14, 7815–7828. doi: 10.3390/ijms14047815
Bizuneh, G. K. (2020). The chemical diversity and biological activities of phytoalexins. Adv. Tradit. Med. 21, 1–13. doi: 10.1007/s13596-020-00442-w
Cruz-Cano, R., and Lee, M. L. T. (2014). Fast regularized canonical correlation analysis. Comput. Stat. Data Anal. 70, 88–100. doi: 10.1016/j.csda.2013.09.020
Darren, K., Matt, C., Robert, B., David, A., and Parag, M. (2008). ProteoWizard: open source software for rapid proteomics tools development. Bioinformatics (Oxford, England) 24, 2534–2536. doi: 10.1093/bioinformatics/btn323
De Carvalho, C. C. R., and Caramujo, M. J. (2018). The various roles of fatty acids. Molecules (Basel, Switzerland) 23:2583. doi: 10.3390/molecules23102583
Dietz, K. J., Vogel, M. O., and Viehhauser, A. (2010). AP2/EREBP transcription factors are part of gene regulatory networks and integrate metabolic, hormonal and environmental signals in stress acclimation and retrograde signalling. Protoplasma 245, 3–14. doi: 10.1007/s00709-010-0142-8
Domingo-Almenara, X., and Siuzdak, G. (2020). Metabolomics data processing using XCMS. Methods Mol. Biol. (Clifton, N.J.) 2104, 11–24. doi: 10.1007/978-1-0716-0239-3_2
Du, Y. W., Jia, H. C., Yang, Z., Wang, S. H., Liu, Y. Y., Ma, H. Y., et al. (2023). Sufficient coumarin accumulation improves apple resistance to Cytospora mali under high potassium status. Plant Physiol. 192, 1396–1419. doi: 10.1093/plphys/kiad184
Feng, H., Wang, C. L., He, Y. T., Tang, L., Han, P. L., Liang, J. H., et al. (2023). Apple Valsa canker: insights into pathogenesis and disease control. Phytopathol. Res. 5:45. doi: 10.1186/s42483-023-00200-1
Geng, D. L., Shen, X. M., Xie, Y. P., Yang, Y. S., Bian, R. L., Gao, Y. Q., et al. (2020). Regulation of phenylpropanoid biosynthesis by MdMYB88 and MdMYB124 contributes to pathogen and drought resistance in apple. Horticult. Res. 7:102. doi: 10.1038/s41438-020-0324-2
Han, P. L., Zhang, R. T., Li, R., Li, F. D., Nie, J. J., Xu, M., et al. (2023). MdVQ12 confers resistance to Valsa Mali by regulating MdHDA19 expression in apple. Mol. Plant Pathol. 25:e13411. doi: 10.1111/mpp.13411
Haxim, Y., Kahar, G., Zhang, X. C., Si, Y., Waheed, A., Liu, X., et al. (2022). Genome-wide characterization of the chitinase gene family in wild apple (Malus sieversii) and domesticated apple (Malus domestica) reveals its role in resistance to Valsa Mali. Front. Plant Sci. 13, –1007936. doi: 10.3389/fpls.2022.1007936
Huang, S., Zuo, T., and Ni, W. Z. (2020). Important roles of glycinebetaine in stabilizing the struture and function of the photosystem II complex under abiotic stersses. Planta 251:36. doi: 10.1007/s00425-019-03330-z
Jiang, J. J., Ma, S. H., Ye, N. H., Jiang, M., Cao, J. S., and Zhang, J. H. (2017). WRKY transcription factors in plant responses to stresses. J. Integr. Plant Biol. 59, 86–101. doi: 10.1111/jipb.12513
Jokinen, K., Salovaara, A. K., Wasonga, D. O., Edelmann, M., Simpura, I., and Mäkelä Pirjo, S. A. (2022). Root-applied glycinebetaine decreases intrate accumulation and improves quality in hydroponi cally grown lettuce. Food Chem. 366:130558. doi: 10.1016/j.foodchem.2021.130558
Ke, X. W., Yin, Z. Y., Song, N., Dai, Q. Q., Voegele, R. T., Liu, Y. Y., et al. (2014). Transcriptome profiling to identify genes involved in pathogenicity of Valsa Mali on apple tree. Fungal Genet. Biol. 68, 31–38. doi: 10.1016/j.fgb.2014.04.004
Kurihara, J., and Yamana, T. (2022). Detection of apple Valsa canker based on hyperspectral imaging. Remote Sens. 14:1420. doi: 10.3390/rs14061420
Li, C. R., Jing, Y. Y., Cheng, L., Si, Z. G., Mou, Z. Y., Niu, D. S., et al. (2023). The antifungal activity of trans-cinnamic acid and its priming effect in apple in response to Valsa Mali. Plant Pathol. 72, 1595–1603. doi: 10.1111/ppa.13786
Li, P. Q., Ruan, Z., Fei, Z. X., Yan, J. J., and Tang, G. H. (2021). Integrated transcriptome and metabolome analysis revealed that flavonoid biosynthesis may dominate the resistance of Zanthoxylum bungeanum against stem canker. J. Agric. Food Chem. 69, 6360–6378. doi: 10.1021/acs.jafc.1c00357
Li, X. L., Su, Q. F., Jia, R. J., Wang, Z. D., Fu, J. H., Guo, J. H., et al. (2023). Comparison of cell wall changes of two different types of apple cultivars during fruit development and ripening. J. Integr. Agric. 22, 2705–2718. doi: 10.1016/j.jia.2023.07.019
Li, C. Y., Zhang, T. P., Feng, P. W., Li, D. X., Brestic, M., Liu, Y., et al. (2021). Genetic engineering of glycinebetaine synthesis enhances cadmium tolerance in BADH-transgenic tobacco plants viareducing cadmium stress damage. Environ. Exp. Bot. 191:104602. doi: 10.1016/j.envexpbot.2021.104602
Liu, Y. T., Xu, R. X., Wang, H. T., Shi, Y. J., Li, C. Y., Ma, F. W., et al. (2022). Exogenous glycine betaine improved the resistance of apple to Glomerella leaf spot. J. Fruit Sci. 39, 1252–1261. doi: 10.13935/j.cnki.gsxb.20210566
Liu, X. M., Yuan, X., Zhang, Z. J., Wang, Y. P., Ma, F. W., and Li, C. (2021). Dopamine enhances the resistance of apple to Valsa Mali infection. Phytopathology 112, 1141–1151. doi: 10.1094/PHYTO-08-21-0328-R
Liu, X. Z., Zhang, Y. H., Zhou, T., Li, X. S., Wen, X. J., and Zhang, D. Y. (2021). Full-length transcriptome-wide characteristic and functional identification of WRKY family in Malus sieversii during the Valsa canker disease response. Forests 12:790. doi: 10.3390/f12060790
Lui, S., Luo, C. G., Zhu, L. M., Sha, R. H., Qu, S. C., Cai, B. H., et al. (2017). Identification and expression analysis of WRKY transcription factor genes in response to fungal pathogen and hormone treatments in apple (Malus domestica). J. Plant Biol. 60, 215–230. doi: 10.1007/s12374-016-0577-3
Ma, Q., Ju, M. Y., Liu, Q., Zhang, L., Li, Z. N., and Sun, P. P. (2020). Identification of the pathogenic fungi associated with apple valsa cancer in Inner Mongolia, China. J. Fruit Sci. 37, 714–722. doi: 10.13925/j.cnki.gsxb.20190474
Mansoor, S., Sakina, A., Mir, M. A., Mir, J. I., Wani, A. A., Un Nabi, S., et al. (2023). Elucidating the role of reactive oxygen species metabolism and phenylpropanoid pathway during an incompatible interaction between apple-Venturia inaequalis host-pathosystem. S. Afr. J. Bot. 160, 428–436. doi: 10.1016/j.sajb.2023.07.020
Mao, X., Wang, C., Lv, Q., Tian, Y., Wang, D., Chen, B., et al. (2021). Cyclic nucleotide gated channel genes (CNGCs) in Rosaceae: genome-wide annotation, evolution and the roles on Valsa canker resistance. Plant Cell Rep. 40, 2369–2382. doi: 10.1007/s00299-021-02778-2
McCarthy, D. J., Chen, Y., and Smyth, G. K. (2012). Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 40, 4288–4297. doi: 10.1093/nar/gks042
Meng, X. L., Qi, X. H., Han, Z. Y., Guo, Y. B., Wang, Y. N., Hu, T. L., et al. (2019). Latent infection of Valsa Mali in the seeds, seedlings and twigs of crabapple and apple trees is a potential inoculum source of Valsa canker. Sci. Rep. 9, 1–10. doi: 10.1038/s41598-019-44228-w
Moneo-Sánchez, M., Izquierdo, L., Martín, I., Hernández-Nistal, J., Albornos, L., Dopico, B., et al. (2018). Knockout mutants of Arabidopsis thaliana β-galactosidase. Modifications in the cell wall saccharides and enzymatic activities. Biol. Plant. 62, 80–88. doi: 10.1007/s10535-017-0739-2
Morris, J. H., Apeltsin, L., Newman, A. M., Baumbach, J., Wittkop, T., Su, G., et al. (2011). ClusterMaker: a multi-algorithm clustering plugin for Cytoscape. BMC Bioinf. 12:436. doi: 10.1186/1471-2105-12-436
Muñoz-Hoyos, L., and Stam, R. (2023). Metabolomics in plant pathogen-defence: from single molecules to large scale analysis. Phytopathology 113, 760–770. doi: 10.1094/PHYTO-11-22-0415-FI
Per, T. S., Khan, M. I. R., Anjum, N. A., Masood, A., Hussain, S. J., and Khan, N. A. (2018). Jasmonates in plants under abiotic stresses: crosstalk with other phytohormones matters. Environ. Exp. Bot. 145, 104–120. doi: 10.1016/j.envexpbot.2017.11.004
Phukan, U. J., Jeena, G. S., and Shukla, R. K. (2016). WRKY transcription factors: molecular regulation and stress responses in plants. Front. Plant Sci. 7:760. doi: 10.3389/fpls.2016.00760
Rohart, F., Gautier, B., Singh, A., and Lê Cao, K. A. (2017). Mix omics: an R package for 'omics feature selection and multiple data integration. PLoS Comput. Biol. 13:e1005752. doi: 10.1371/journal.pcbi.1005752
Sebastian, W. (2022). Cell Wall signaling in plant development and defense. Annu. Rev. Plant Biol. 73, 323–353. doi: 10.1146/annurev-arplant-102820-095312
Sung, W. S., and Lee, D. G. (2010). Antifungal action of chlorogenic acid against pathogenic fungi, mediated by membrane disruption. Pure Appl. Chem. 82, 219–226. doi: 10.1351/PAC-CON-09-01-08
Tsuda, K., and Somssich, I. E. (2015). Transcriptional networks in plant immunity. New Phytol. 206, 932–947. doi: 10.1111/nph.13286
Utriainen, M., and Morris, J. H. (2023). Cluster maker 2: a major update to cluster maker, a multi-algorithm clustering app for Cytoscape. BMC Bioinformatics 24:134. doi: 10.1186/s12859-023-05225-z
Varet, H., Brillet-Guéguen, L., Coppée, J., and Dillies, M. (2017). SARTools: a DESeq2-and edge R-based R pipeline for comprehensive differential analysis of RNA-Seq data. PLoS One 11:e0157022. doi: 10.1371/journal.pone.0157022
Wang, S. S., Chen, K., Zhang, J. Y., Wang, J. Q., Li, H. S., Yang, X. Y., et al. (2022). Genome-wide characterization of MATE family members in Cucumis melo L. and their expression profiles in response to abiotic and biotic stress. Horticult. Plant J. 8, 474–488. doi: 10.1016/j.hpj.2022.05.004
Wang, D. L., Liu, B., Ma, Z. Q., Feng, J. T., and Yan, H. (2020). Reticine A, a new potent natural elicitor: isolation from the fruit peel of Citrus reticulate and induction of systemic resistance against tobacco mosaic virus and other plant fungal diseases. Pest Manag. Sci. 77, 354–364. doi: 10.1002/ps.6025
Wang, X. L., Shi, C. M., Gleason, M. L., and Huang, L. L. (2020). Fungal species associated with apple Valsa canker in East Asia. Phytopathol. Res. 2:14. doi: 10.1186/s42483-020-00076-5
Wen, X. Z., Yuan, J. X., Bozorov, T. A., Waheed, A., Kahar, G., Haxim, Y., et al. (2023). An efficient screening system of disease-resistant genes from wild apple, Malus sieversii in response to Valsa Mali pathogenic fungus. Plant Methods 19:138. doi: 10.1186/s13007-023-01115-w
Wu, T. Z., Hu, E. Q., Xu, S. B., Chen, M. J., Guo, P. F., Dai, Z. H., et al. (2021). Cluster profiler 4.0: a universal enrichment tool for interpreting omics data. Innovations 2:100141. doi: 10.1016/j.xinn.2021.100141
Xin, L., Zhang, R., Wang, X. P., Liu, X. X., Wang, Y. C., Qi, P. R., et al. (2023). Extracellular and intracellular infection of Botryosphaeria dothidea and resistance mechanism in apple cells. Horticult. Plant J. 9, 209–223. doi: 10.1016/j.hpj.2022.05.001
Yin, Z. Y., Ke, X. W., Kang, Z. S., and Huang, L. L. (2016). Apple resistance responses against Valsa Mali revealed by transcriptomics analyses. Physiol. Mol. Plant Pathol. 93, 85–92. doi: 10.1016/j.pmpp.2016.01.004
Yin, Z. Y., Liu, H. Q., Li, Z. P., Ke, X. W., Dou, D. L., Gao, X. N., et al. (2015). Genome sequence of Valsa canker pathogens uncovers a potential adaptation of colonization of woody bark. New Phytol. 208, 1202–1216. doi: 10.1111/nph.13544
Yoav, B., and Yosef, H. (2000). On the adaptive control of the false discovery rate in multiple testing with independent statistics. J. Educ. Behav. Stat. 25:60. doi: 10.2307/1165312
Yu, X. Y., Bi, Y., Yan, L., Liu, X., Wang, Y., Shen, K. P., et al. (2016). Activation of phenylpropanoid pathway and PR of potato tuber against fusarium sulphureum by fungal elicitor from Trichothecium roseum. World J. Microbiol. Biotechnol. 32:142. doi: 10.1007/s11274-016-2108-2
Zang, R., Huang, L., Kang, Z., and Wang, X. (2007). Biological characteristics and pathogenicity of different isolates of Cytospora spp. isolated from apple trees in Shaanxi province. Acta Phytopathol. Sin. 37, 343–351. doi: 10.3321/j.issn:0412-0914.2007.04.002
Zhao, D., Duo, H., Lv, Q. Q., Tian, Y. Z., and Zuo, C. W. (2021). β-Galactosidase gene family responds to the signal of Valsa mail in apple. Plant Physiol. J. 57, 1319–1328. doi: 10.13592/j.cnki.ppj.2020.0515
Zhao, L. N., Shu, Y. L., Liang, L. Y., Wang, Y. J., Esa, A. G., Zhang, X. Y., et al. (2022). Integrated transcriptomic and metabonomic analysis reveal mechanisms of disease resistance in apples induced by Wickerhamomyces anomalus. Biol. Control 173:105005. doi: 10.1016/j.biocontrol.2022.105005
Zhou, K., Hu, L. Y., Li, Y. T. S., Chen, X. F., Zhang, Z. J., Liu, B. B., et al. (2019). MdUGT88F1-mediated Phloridzin biosynthesis regulates apple development and Valsa canker resistance. Plant Physiol. 180, 2290–2305. doi: 10.1104/pp.19.00494
Zhou, K., Hu, L. Y., Yue, H., Zhang, Z. J., Zhang, J. Y., Gong, X. Q., et al. (2021). MdUGT88F1-mediated phloridzin biosynthesis coordinates carbon and nitrogen accumulation in apple. J. Exp. Bot. 73, 886–902. doi: 10.1093/jxb/erab410
Zhou, G. Y., Pang, Z. Q., Lu, Y., Jessica, E., and Xia, J. G. (2022). Omics net 2.0: a web-based platform for multi-omics integration and network visual analytics. Nucleic Acids Res. 50, W527–W533. doi: 10.1093/nar/gkac376
Keywords: apple, Cytospora mali, transcriptome, metabolome, conjoint analysis
Citation: Zhao J, Guo Y, Li Z, Shi Y and Sun P (2024) Integrated transcriptome and metabolome profiling reveals mechanisms underlying the infection of Cytospora mali in “Jin Hong” branches. Front. Microbiol. 15:1394447. doi: 10.3389/fmicb.2024.1394447
Edited by:
Yonglin Wang, Beijing Forestry University, ChinaCopyright © 2024 Zhao, Guo, Li, Shi and Sun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Pingping Sun, sunpingping@imau.edu.cn