Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 25 October 2022
Sec. Plant Symbiotic Interactions
This article is part of the Research Topic Arbuscular Mycorrhiza-Mediated Augmentation of Plant Secondary Metabolite Production View all 5 articles

Transcriptomics and metabolomics reveal effect of arbuscular mycorrhizal fungi on growth and development of apple plants

Shan JingShan JingYuchao LiYuchao LiLingcheng ZhuLingcheng ZhuJing SuJing SuTianyi YangTianyi YangBowen LiuBowen LiuBaiquan MaBaiquan MaFengwang MaFengwang MaMingjun Li*Mingjun Li*Manrang Zhang*Manrang Zhang*
  • State Key Laboratory of Crop Stress Biology for Arid Areas, College of Horticulture, Shaanxi Key Laboratory of Apple, Northwest A&F University, Yangling, Shaanxi, China

Arbuscular mycorrhizal fungi (AMF) and plants form a symbiotic relationship that promotes plant growth and development. However, the regulatory mechanisms through which AMF promote plant growth and development are largely unexplored. In this study, the apple rootstock M26 was assessed physiologically, transcriptionally and metabolically when grown with and without AMF inoculation. AMF significantly promoted the number of lateral root (LR) increase and shoot elongation. Root transcriptomic and metabolic data showed that AMF promoted lateral root development mainly by affecting glucose metabolism, fatty acid metabolism, and hormone metabolism. Shoot transcriptomic and metabolic data showed that AMF promoted shoot elongation mainly by affecting hormone metabolism and the expression of genes associated with cell morphogenesis. To investigate whether shoot elongation is caused by root development, we analyzed the root/shoot dry weight ratio. There was a correlation between shoot growth and root development, but analysis of root and shoot metabolites showed that the regulation of AMF on plant shoot metabolites is independent of root growth. Our study bridged the gap in the field of growth and development related to AMF.

Introduction

Arbuscular mycorrhizal fungi (AMF) are widely distributed in the ecosystem. Symbiosis with arbuscular mycorrhizal (AM) from the Order Glomales is observed in 80% to 90% of terrestrial plants worldwide (Gadkar et al., 2001). AMF can improve plant resistance to biotic stressors, such as insect pests and bacterial and fungal diseases, and abiotic stressors, such as drought, salt, and heavy metals (Karthikeyan et al., 2016; Gao et al., 2020; Huang et al., 2020; Adeyemi et al., 2021). In addition, AMF promote the growth of plants, including perennial trees like apple, trifoliate orange, and catalpa (Hosseini and Gharghani, 2015; Liu et al., 2018; Chen et al., 2020). The mechanism through which AMF improves plant stress tolerance has been extensively studied, but there is not yet intensive study on how AMF regulates plant growth and development.

AMF is a specialized, obligate biotrophic fungus that must inhabit a plant root. Once established, this symbiosis can promote root development (Hosseini and Gharghani, 2015; Liu et al., 2018; Huang et al., 2020). It remains ambiguous how AMF promotes root development in plants. Numerous researchers have proposed that AMF symbiosis promotes the uptake of mineral nutrients, such as phosphate and nitrogen, which would then promote root development (Bonfante and Genre, 2010; Andrino et al., 2021). Wang et al. (2021) concluded that auxin plays a critical role in AMF promotion of plant growth. Although it is commonly understood that AMF promote root development, the molecular mechanisms underlying this phenotype are still poorly studied.

As the same time, AMF also facilitate the development of the shoot. In apple and tomato, inoculated plants have a greater plant height, stem diameter and dry shoot weight (Gao et al., 2020; Huang et al., 2020; Wang et al., 2021). In a two-year field experiment, AMF promoted the plant height of Helianthus tuberosus in a comparison between mycorrhizal and non-mycorrhizal plants (Nacoon et al., 2021). AMF also affects the hormone content in stem, including the levels of IAA, cytokinins (CTK), and gibberellic acid (GA) (Dugassa et al., 1996; Barker and Tagu, 2000). It has long been believed that the growth potential of the above ground parts of the plant is improved through the increased root development. However, it has not been confirmed if the promotion of shoot growth is due to better root development or due to the AMF.

Apples (Malus × domestica Borkh.) are among the most valuable horticultural fruit crops cultivated and is of great interest to horticultural researchers all over the world (Wang et al., 2018). In apple, AMF improves the resistance to stress and promotes plant growth and development (Huang et al., 2020; Gao et al., 2020). Metabolomics is a widely used and highly recognized analytical method for studying plant growth and development (Fàbregas et al., 2018; Yang et al., 2022). This study is an association analysis of metabolomic and transcriptomic datasets between the roots and shoots of apple plants before and after inoculation with AMF. Based on our analysis, we provide a partial explanation of how AMF regulates root and shoot growth and development at both the metabolic and transcriptional levels. We found that the promotion of shoot growth by AMF was not merely due to better root development, as AMF also directly regulated shoot growth and development. These results may provide a new direction for the study of AMF regulation of plant growth and development, that AMF may directly affect shoot development through some pathways, such as effectors, RNAs or metabolites of AMF.

Results

Plant growth and physiological responses to inoculation with AMF

To explore the effects of AMF on plant growth and development, we inoculated young plants of the cultivar M.26 with AMF. After 30 days, the average height of the plants inoculated with the mycorrhizae was higher than that of the uninoculated plants (Figures 1A, B). Sixty days after inoculating with AMF, the average height of the inoculated plants was 65.2 cm, while that of the uninoculated plants was 49.1 cm, and the height of the plants was markedly different. Moreover, the internode length of the mycorrhizal-inoculated plants was 1.4 times that of the uninoculated plants and the number of internodes was 1.28 times that of the uninoculated plants (Figures 1C, D). Interestingly, we found that there were no significant differences in the chlorophyll contents and photosynthetic rates between the mycorrhizal plants and non-mycorrhizal plants (Figures 1E, F). At the same time, a good symbiotic relationship between the AMF and the plants was observed 60 days after inoculation (Figure 1G). Furthermore, the number of roots in mycorrhizal plants was 215, while that in the uninoculated plants was 109, and the number of inoculated roots was significantly increased 60 days of inoculation with AMF (Figure 1H). The total root length was 3379.9 cm in mycorrhizal plants and 1713.8 cm in uninoculated plants (Figures 1I, J). There was no significant difference in the root/shoot (dry weight) ratio between mycorrhizal and non-mycorrhizal plants (Figure 1K). These results showed that AMF promoted shoot growth and LR development.

FIGURE 1
www.frontiersin.org

Figure 1 Apple plant growth and physiological responses to AM symbiosis. (A) The phenotype of mycorrhizal and non-mycorrhizal apple plants inoculated with Rhizophagus irregularis for 60 days. Scale bar = 7 cm. CK, non-mycorrhizal inoculated plants. AMF, mycorrhizal-inoculated plants. (B–D) The plant height, internode length, and internode number of apple plants uninoculated or inoculated with R. irregularis after 15, 30, 45 and 60 days, respectively. (E, F) The net photosynthetic rate and chlorophyll content of R. irregularis mycorrhizae uninoculated or inoculated apple plants after 60 days. (G) Staining of mycorrhizae in apple plants uninoculated or inoculated with R. irregularis after 60 days. (H) The root phenotype of uninoculated or inoculated apple plants after 60 days. Scale bar = 5 cm. (I-K) The number of roots, total root length, and root/shoot ratio of apple plants uninoculated or inoculated with R. irregularis after 60 days. Values are the means ± SD of three biological replicates. * Indicates a significant difference at P < 0.05 level.

Overall analysis of metabolomics in root and shoot of apple plants with AM symbiosis

The metabolome is closely related to plant growth and development. To explore whether AMF affected the accumulation of metabolites, we analyzed the metabolome in the root and shoot. The quality control showed that the total ion ground curve of metabolite detection had a high overlap, indicating that the signal stability of the mass spectrometry was good (Supplemental Figures 1A–D). Supplemental Figure 1E shows that the intra-group repeatability was good. Principal component analysis (PCA) revealed that components 1 and 2 explained 77.58% and 11.44% of the variability, respectively, indicating significant metabolic diversity between roots and shoots (Supplemental Figure 1F).

After filtering by signal-to-noise ratio, 814 metabolites were detected in the root, with 48 differential metabolites between uninoculated and AM-inoculated roots. Among the differential metabolites, 15 were significantly decreased and 33 were significantly increased (Supplemental Figure 1G). In shoots, 1142 metabolites were detected, of which 173 metabolites were significantly down-regulated and 63 metabolites were significantly up-regulated (DESeq2 analysis with a |log 2(fold change)|≥1 and Variable Importance in Projection (VIP) ≥1) (Supplemental Figure 1H). The 824 identified metabolites in the roots included 133 lipids, 66 organic acids, 59 saccharides and alcohols, 43 nucleotides and derivatives, 59 amino acids and derivatives, 62 terpenoids, 46 alkaloids, 114 flavonoids, 37 lignans and coumarins, 13 tannins, 139 phenolic acids, 11 quinones, and 32 other compounds that did not fit into these main classes (Supplemental Data 1). The 1142 identified metabolites of shoot included 137 lipids, 89 organic acids, 53 nucleotides and derivatives, 85 amino acids and derivatives, 73 terpenoids, 61 alkaloids, 255 flavonoids, 31 lignans and coumarins, 14 tannins, 208 phenolic acids, 73 saccharides and alcohols, and 63 other compounds (Supplemental Data 2).

Analysis of differential metabolites in apple roots and shoots in symbiosis with AM

In the AM-inoculated roots, there were 48 differential metabolites, among which there were 30 primary metabolites (14 lipids, 8 sugars, 5 organic acids, 2 others, and 1 nucleotide) and 18 secondary metabolites (9 phenolic acids, 6 terpenoids, 2 lignans and coumarins, and 1 flavonoid). Among these differential metabolites, lipids, sugars, organic acids, terpenoids, and phenolic acids were visibly altered with AM symbiosis, accounting for 87.5% of all differential metabolites. Lipid-containing substances were the most abundant, with two (eicosenoic acid and LysoPC 22:4) down-regulated, while the other 12 were up-regulated (Figure 2A). Changes in the phenolic acids and sugars were comparatively significant, accounting for 18.8% and 16.7% of differential metabolites, respectively. Many organic acids and terpenoids were also changed (Figure 2B). We also measured the levels of hormones in the roots. The contents of abscisic acid, indole-3-acetyl-L-aspartic acid, and Indole-3-acetyl glycine were increased over 1-fold, while 3-Indoleacrylic acid, dihydrozeatin-7-glucoside, and salicylic acid 2-O-β-glucoside were decreased 61.8% (Figure 2C). These results suggested that AMF might affect the accumulation of lipids, sugars, organic acids, terpenoids, phenolic acids, and hormones in plant roots.

FIGURE 2
www.frontiersin.org

Figure 2 Heatmap of differential metabolites in roots and shoots of apple plants inoculated with R. irregularis after 60 days, with the baseline set as the level in uninoculated plants. (A) Thirty differential primary metabolites were identified in roots of apple plants. (B) Eighteen differential secondary metabolites were identified in roots of apple plants. (C) Six differential hormonal metabolites in roots of apple plants. (D) Fifty-one differential primary metabolites in shoots of apple plants. (E) One hundred eight-five differential secondary metabolites in shoots of apple plants. (F) Fourteen differential hormonal metabolites in shoots of apple plants. Primary and secondary metabolites are shown as the ratio of relative content, while hormone metabolites were the ratio of absolute content. In the case of each hormonal metabolite, the metabolites were only detected in either the uninoculated or inoculated plant, but not in both sample sets, which was considered to be significantly regulated (|log 2 FC| = 2).

In the apple shoots, there were 236 differential metabolites (Figures 2E, F), among which 51 were primary metabolites (17 lipids, 12 amino acids and derivatives, 10 organic acid, 9 others, and 3 nucleotides and derivatives) and 185 were secondary metabolites (82 flavonoids, 57 phenolic acids, 16 terpenoids, 15 alkaloids, 10 lignans and coumarins, and 5 tannins). There were three times more differential secondary metabolites than differential primary metabolites. Significant changes were seen for flavonoids, phenolic acids and lipids, accounting for 34.7%, 24.2%, and 9.1% of the differential metabolites, respectively. We also found 17 metabolites that showed specific accumulation in inoculated plants (5 flavonoids, 4 lipids, 2 phenolic acids, 2 organic acid, 1 terpenoids, 1 lignans and coumarins, 1 amino acid derivative, and 1 others), conversely, 21 metabolites that showed specific accumulation in the non-inoculated plants (7 flavonoids, 4 alkaloids, 3 phenolic acids, 2 organic acid, 2 others, 1 amino acid derivative, 1 tannin, and 1 lipids) (Figures 2D, E). These data collectively indicated that the most accumulated metabolites were flavonoids, phenolic acids, and lipids. Again, there were significant changes in the roots in hormones, such as auxin, cytokinin (CTK) and gibberellin (GA), accounting for 78.6% of the total differential hormones. Of the 14 differential hormone metabolites, four metabolites were related to GA, four to auxin, three to CTKs and three to other hormones (Figure 2F). Altogether, our results indicated that symbiosis with AMF might affect the accumulation of flavonoids, phenolic acids, lipids, and hormones in the shoot.

Analysis of transcripts by RNA-seq in apple roots and shoots in symbiosis with AM

In order to investigate whether metabolic pathways were transcriptionally regulated in apple inoculated with AMF, we subjected the RNA extracted from inoculated and uninoculated roots and shoots to RNA-sequencing. The intra-group repeatability was good, ranging from 0.906 to 0.999 (Supplemental Figure 2A). The PCA showed that components 1 and 2 were 91.7% and 5.6% of the variability, respectively, indicating significant transcriptional diversity between roots and shoots (Supplemental Figure 2B). There were 2371 differently expressed genes (DEGs) in root (1140 were up-regulated and 1231 were down-regulated; |log 2 FC| > 1; P value < 0.05; Supplemental Figure 2C; Supplemental Data 3) and 1142 DEGs in shoots (319 up-regulated and 823 down-regulated; |log 2 FC| > 1; Supplemental Figure 2D; Supplemental Data 3). To verify the accuracy of the transcriptome results, we randomly selected 4 DEGs (ARF19, PRP1, IAA17, and YUCCA6) for RT-qPCR using RNA from roots and shoots. The results showed that the trends for the relative expression of genes were consistent between the RNA-seq and the RT-qPCR (Supplemental Figures 2E, F).

Analysis of differentially expressed genes in roots and shoots of apple plants with AM symbiosis

The differential metabolites and the differentially expressed genes identified from the mycorrhizal-inoculated and uninoculated plants were combined for an enrichment analysis. The differential metabolites and the differential genes were enriched in the same pathways, thereby indicating that the transcriptome and metabolome corresponded (Figure 3A and Supplemental Data 4). And the datasets were next analyzed using MapMan. In roots, of the 120 differential genes, 99 genes were related to primary metabolites and 21 genes were related to secondary metabolites. Upon further analysis, the most changed genes were related to lipid metabolism (46) and carbohydrate metabolism (42). Of the genes related to lipid metabolism, 23 genes were related to glycerolipid metabolism, 12 genes were related to fatty acid metabolism, and 11 genes were in other pathways (Figures 3B, C). About 42 genes were related to carbohydrate metabolism, 33 genes were related to sugar metabolism, and 9 genes were related to fermentation. The 72 differential genes related to plant hormones were screened out, with genes related to auxin, CK, GA, abscisic acid (ABA) and strigolactone (SL) showing significant changes. Among them, 15 genes were related to auxin, 12 to CK, 19 to GA, 8 to ABA, and 6 to SL. There were a total of 12 genes related to brassinosteroid (BR), ethylene (ETH), and jasmonic acid (JA) metabolism (Figure 3D). Collectively, significant changes have taken place in the levels of transcripts related to lipid metabolism, carbohydrate metabolism and hormone metabolism, suggesting that AMF affected the regulation of these genes.

FIGURE 3
www.frontiersin.org

Figure 3 DEGs in root and shoot of apple plants inoculated with R. irregularis for 60 days (according to differential metabolites and MapMan). (A) Enrichment analysis of related differential genes and differential metabolites in root and shoot of transcriptome and metabolome. (B–D), Differential genes related to metabolites in root, (B–D) with ninety-nine differential genes related to primary metabolites, twenty-one differential genes related to secondary metabolites, and seventy-three differential genes related to hormonal metabolites in root, respectively. (E–G), Differential genes related to metabolites in shoot, (E) with forty-nine differential genes related to primary metabolites, (F) twenty differential genes related to secondary metabolites, and (G) forty-two differential genes related to hormonal metabolites in shoot. Expression-fold changes of DEGs (|log 2 FC| > 1, FDR < 0.05) were represented by color gradient from pink (induced) to blue (reduced).

In shoots, of the 68 differential genes, 48 genes were related to primary metabolites (32 genes were related to lipid metabolism, 11 genes to carbohydrate metabolism, and 5 genes to amino acid metabolism) and 20 genes were related to secondary metabolites. Of the 32 genes related to lipid metabolism, 29 genes were down-regulated, and only 4 were up-regulated (related to fatty acid metabolism). Genes related to glycerolipid metabolism (11), phytosterol metabolism (3), and lipid trafficking (2) were all down-regulated (Figures 3E, F). The 42 differential genes related to hormone action were screened out, with genes related to JA (20), GA (7), ABA (5), CTK (3) and BR (3) showing significant changes (Figure 3G). Of the 20 genes related to JA, the 8 genes related to JA biosynthesis were down-regulated, and the 7 genes related to JA conjugation and degradation were also down-regulated. This may explain why there was no significant difference in the JA content between the inoculated and uninoculated plants, with only a 1.5-fold down-regulation seen with inoculation (Supplemental Figure 5). These data indicated that AMF might affect the expression and/or turnover of genes related to lipid and hormone metabolism.

Carbohydrate metabolism and sugar accumulation in roots of apple plants in symbiosis with AM

The results of the joint metabolomic and transcriptomic analyses showed that sugar, hormone, and fatty acid pathways undergo significant changes when inoculated with AMF. Sugars, hormones, and fatty acids have been reported to influence root development (Trinh et al., 2019; Yamauchi et al., 2019; Valifard et al., 2021). To further explore how AMF changes of sugar, fatty acid and hormone contents, the sugar, fatty acid, and hormone metabolic pathways were mapped and integrated with the transcriptome and metabolome data (Figure 4).

FIGURE 4
www.frontiersin.org

Figure 4 Main carbohydrate metabolism and sugar accumulation pathways in roots of apple plants inoculated with R. irregularis for 60 days. CWINV, cell wall invertase; FK, fructokinase; F6P, fructose-6-phosphate; Fru, D-fructose; Glc, D-glucose; G6P, glucose-6-phosphate; Tre, trehalose; Gal, galactose; GPT, glucose-6-phosphatetransporter; SDH, sorbitol dehydrogenese; HK, hexokinase; NINV, neutral invertase; PFK, phosphofructokinase; PGI, phosphoglucoisomerase; pGLCT, plastidic Glc transporter;PGM, phosphoglucomutase; SPP, sucrose-phosphate phosphatase; SPS, sucrose-phosphate synthase; SUSY, sucrose synthase; vAINV, vacuolar acidinvertase; vGT, vacuole glucose transporter; UGP, UDPG-pyrophosphorylase; BMY, beta amylase; HGL, heteroglycan glucosidase; PHS, α-glucanphosphorylase; AMY, α- amylase; α-GWD, α-glucan water dikinase; PWD, phosphoglucan waterdikinase; SBE, starch branching enzyme; SS, starchsynthase; TST, tonoplast sugar transporter; TPS, trehalose-phosphate synthase; α-gal, α-galactosidase; GALK, galactokinase; THL, trehalase; T6P, trehalose-6-phosphate.

Among the sugars, the relative contents of glucose (Glc), fructose (Fru), sucrose (Suc), trehalose (Tre) and sorbitol (Sor) increased over 50% after inoculating AMF. The largest change in relative content was for Glc. This accumulation of Glc might be caused through modification of two pathways (Figure 4): the increased relative content of trehalose and the increased relative expression of neutral invertase (NINV: MD00G1092700). The content of Suc was increased 90% (Supplemental Figure 3), although the transcripts for sucrose phosphate synthase (SPS) were not significantly up-regulated, there was a high, up-regulated relative level of expression of three SPS1 genes (MD02G1022300, MD15G1127900 and MD15G1164900). In addition, the content of Fru increased by 90% after inoculation with AMF (Supplemental Figure 3), the transcript levels of sorbitol dehydrogenase (SDH) showed a downward trend as a whole, and another NINV (MD00G109270) was significantly up-regulated. These data indicated that Fru may accumulate through an increased cleavage of Suc caused by AMF symbiosis.

The responses of fatty acid synthesis and transport genes to AM symbiosis

The relative content of fatty acids showed an overall upward trend during AMF symbiosis. The relative expression levels of acyl-ACP thioesterase A (FatA) and acyl-ACP thioesterase B (FatB) showed no significant differences, but transcripts relevant to fatty acid synthesis, such as ketoacyl-ACP synthase (KAS: MD06G1110900, MD14G1131500, and MD09G1022600), were significantly up-regulated. The relative content of 14 diacylglycerol (DAG) was not significantly different, and the relative transcript levels of phosphatase (PAP) were little differences (Supplemental Figure 4).

Synthesis and signal transduction of hormones in roots of apple plants in symbiosis with AM

The transcriptome and metabolome both showed significant changes in auxin, CTK, GA and SL pathways, all of which are closely related to root growth (Yaxley et al., 2001; Zhou et al., 2011; Guillotin et al., 2017; Yamauchi et al., 2019). The content of indole-3-acetic acid (IAA) was increased (Supplemental Figure 5 and Figure 5A), which correlated with the increased levels of transcripts within the whole YUCCA pathway. Conversely, the genes in the indole-3-acetamide (IAM) pathway showed a downward expression trend as a whole. Almost all genes related to auxin signal transduction were down-regulated, with only ARF18 (MD07G1174000) and ARF2 (MD11G1297900) significantly up-regulated. These data suggested that symbiosis with AMF might mainly change the expression of genes in the YUCCA pathway, affecting the synthesis of auxin.

FIGURE 5
www.frontiersin.org

Figure 5 Synthesis and transduction pathways of hormones in roots of apple plants inoculated with R. irregularis for 60 days. (A) AMI1, AMIDASE 1; IAA, indole-3-acetic acid; IAAld, indole-3-acetaldehyde; IAM, indole-3-acetamide; IAN, indole-3-acetonitrile; IAOx, indole-3-acetaldoxime; IPyA, indole-3-pyruvic acid; L-Trp, L-tryptophan; NIT, nitrilase; TAR, tryptophan aminotransfearse related; TDC, tryptophan decarboxylase; CYP79B; AUX/IAA, transcriptional regulator family protein; ARF, auxin response factors; TIR1, transport inhibitor response1; (B) CTK, cytokinin; CRFs, cytokinin responsive factors; AHK3 and AHK4/CRE1, cytokinin response 1, 2; AHK2, Arabidopsis histidine kinase 2; type-B ARRs, Arabidopsis response regulators; AHPs, Arabidopsis histidine phosphotransferase proteins; IPT, isopentenyltransferase; CYP735A, cytochrome P450 monooxygenase; tRNA-specific IPT, tRNA-specific isopentenyltransferase; LOG, Lonely guy; iPRTP, iPR trphosphate; iPRDP, iPR diphosphate; iPRMP, iPR monophosphate; iPR, iP riboside; iP, N6-(2-isopentenyl)adenine; tZRTP, tZR 5′-trphosphate; tZRDP, tZR 5′-diphosphate; tZRMP, tZR 5′-monophosphate; tZR, tZ riboside; tZ, trans-zeatin; cZRMP, cZR 5′-monophosphate; cZ, cis-zeatin; DMAPP, Dimethylallylpyrophosphate; ATP, Adenosine triphosphate; ADP, Adenosine diphosphate; AMP, Adenosine monophosphate; HMBPP, 4-hydroxy-3-methyl-2-(E)-butenyl diphosphate; (C) DWARF27, D27; CCD7, carotenoid cleavage dioxygenase 7; MAX1, more axillary growth1; CCD8, carotenoid cleavage dioxygenase 8; SL, 5-deoxylstrigol; CL, carlactone; CLA, carlactonic acid; MeCLA, methy1 carlactonoate; LBO, oxidoreductase; (D) GA, gibberellin; KAO, ent-kaurenoic acid oxidase; KO, ent-kaurene oxidase; GA20ox, GA 20-oxidase; GA3ox, GA 3-oxidase; GA2ox, GA 2-oxidase; GID1, gibberellin insensitive dwarf1; DELLA, Gibberellic acid insensitive.

For CTK, the contents of N6-(2-isopentenyl) adenine (iP) and cZR (CZ) were increased, and trans-zeatin riboside (tZR) and iP riboside (IPR) were decreased in AMF-inoculated apple roots. The form of CTK studied in the root was tZR (Li et al., 2018). In the CTK synthetic pathway, a cytochrome P450 monooxygenase (CYP735A: MD09G1087700) and a cytokinin phosphoribohydrolase (Lonely guy, LOG: MD05G1130400) were significantly down-regulated. Within the CTK signal transduction pathway, almost all genes showed no significant differences, only two Arabidopsis histidine phosphotransferase proteins (AHP: MD04G1212100, MD11G1293900) were significantly changed (Supplemental Figure 5 and Figure 5B). These data indicated that AMF may affect the synthesis of tZR through changing the expression of CTK biosynthesis genes (MD09G1087700, MD05G1130400) in root.

Strigolactone (SL) was also up-regulated. The SL synthesis pathway showed an overall upward trend in expression. Conversely, the genes related to the SL transduction pathway showed downward trends (Supplemental Figure 5 and Figure 5C). These data indicated that AMF may affect the expression of genes related to both synthesis and signal transduction of SL in root.

The content of GA3 was also significantly increased with AMF symbiosis. Genes related to both GA synthesis and signal transduction pathways showed overall up-regulation. In particular, the relative transcript levels of DELLA (MD09G1007200 and MD17G1011100) were up-regulated over 38 times after inoculation with AMF (Supplemental Figure 5 and Figure 5D), suggesting that AMF has a great influence on the transcription level of DELLA.

Combined analysis of root and shoot metabolomes in symbiosis with AM

The environmental effects on plant growth and development have been studied extensively in both the below and above-ground biomass (Niklas, 2005). We found that AMF symbiosis significantly increased plant height (Figure 1). Next, we explored whether this increased shoot growth was caused by root development by analyzing the root to shoot ratio of metabolites, as follows: (AMF R/S - CK R/S)/CK R/S [Changes in root to shoot ratio before and after treatment], (AMF -CK)/CK R [Changes in root before and after treatment], and (AMF -CK)/CK S [Changes in shoot before and after treatment] (Table 1). This analysis revealed that only 26.4% of metabolites has unaltered R/S ratios (-20% < (AMF R/S-CK R/S)/CK R/S < 20%), while 73.6% of metabolites showed large differences in their R/S ratio (AMF R/S-CK R/S)/CK R/S >20%; (AMF R/S-CK R/S)/CK R/S <-20%). Among the metabolites with significantly different R/S ratios, about 40.3% of metabolites showed corresponding opposite trends (R+ S- or R- S+). These data showed that the trends in roots and shoots were very different. Therefore, we speculated that the shoot growth was not only caused by root development, but may also be closely related to AMF. How AMF mediated shoot development was analyzed next (Table 1 and Supplemental Table 2).

TABLE 1
www.frontiersin.org

Table 1 Variation of the root to shoot ratios of various metabolites in apple plants inoculated with R. irregularis for 60 days.

Phenotypic analysis of shoot cells in apple plants in symbiosis with AM for 60 days

Cell development is closely related to plant height (Fu et al., 2002), so we wanted to determine if AMF affected cell development by observing the anatomical structures in the shoots of mycorrhizal and non-mycorrhizal plants. The cross-sections of mycorrhizal plants were not significantly different from non-mycorrhizal plants. In longitudinal sections, the cortex cells were longer and the node spacing increased in shoots with AMF (Supplemental Figure 6A). To confirm this visual observation, we measured the diameter and length of parenchyma cells. There was not a significant difference in the diameter of parenchyma cells between mycorrhizal and non-mycorrhizal plants, while the length of parenchyma cells was increased by 90% with AM symbiosis, which was significantly longer than non-mycorrhizal plants (Supplemental Figure 6B). Statistics of the relative expression levels of genes related to cell development (47 related to cell wall, 4 cell respiration, 4 cytoskeleton organization, and 5 cell division). The genes were related to cell wall, accounting for 78.3%.

Polar cell expansion in different tissues is critical for the development and morphogenesis of plant organs (Fu et al., 2002). Statistics of the relative expression levels of EXPA genes showed the EXPA8 (MD01G1166700 and MD11G1054500), EXPA15 (MD00G1125400) and EXPA2 (MD03G1052400) were significantly up-regulated with AMF symbiosis. Four genes related to fatty acyl omega-hydroxylase (CYP) were significantly down-regulated. These data indicated that the change in cell morphology may be caused by changes in cell wall synthesis (Supplemental Figure 6C and Supplemental Data 3).

Synthesis and signal transduction of hormones in shoots of apple plants in symbiosis with R. irregularis

Combined analysis of the transcriptome and the metabolome showed that there were significant changes related to auxin, CTK and GA, all of which influence shoot development (Ma et al., 2016; Zheng et al., 2018; Zheng et al., 2022). Significant changes related to auxin occurred at the metabolic level. For instance, the content of L-tryptophan (TRP) decreased 63.5%, while indole-3-acetic acid (IAA) increased 34.5% and N-(3-indolylacetyl)-L-valine (IAA-Val), 3-indoleacetonitrile (IAN), and indole-3-acetyl glutamic acid (IAA-Glu) increased over 2.5-fold (Supplemental Figure 7, Figure 2 and Figure 6A). However, the transcriptome showed no significant changes in genes related to synthesis or signal transduction of auxin.

FIGURE 6
www.frontiersin.org

Figure 6 Synthesis, degradation, and transduction pathways of hormones in shoots of apple plants inoculated with R. irregularis for 60 days. (A) AMI1, amidase1; IAA, indole-3-acetic acid; IAAld, indole-3-acetaldehyde; IAM, indole-3-acetamide; IAN, indole-3-acetonitrile; IAOx, indole-3-acetaldoxime; IPyA, indole-3-pyruvic acid; L-Trp, L-tryptophan; NIT, nitrilase; TAR, tryptophan aminotransferase-related; TDC, tryptophan decarboxylase; TIR1, transport inhibitor response1; AUX/IAA, transcriptional regulator family protein; ARF, auxin response factors; (B) CTK, Cytokinin; ABCG, ABC transporter G family member; type-B ARRs, Arabidopsis response regulators; AHPs, Arabidopsis histidine phosphotransferase proteins; (C) GA, gibberellin; KAO, ent-kaurenoic acid oxidase; KO, ent-kaurene oxidase; GA20ox, GA 20-oxidase; GA3ox, GA 3-oxidase; GA2ox, GA 2-oxidase; gibberellin insensitive dwarf1, GID1; Gibberellic acid insensitive, DELLA.

The content of the CTK trans-Zeatin riboside (tZR) increased 56.9% in the shoot of AMF-inoculated apple plants compared to the uninoculated mycorrhizal control plants (Supplemental Figure 7). The relative expression of the ABC transporter G family member (ABCG: MD13G1049400) was significantly up-regulated. In the CTK signal transduction pathway, the AHP was down-regulated, but the B-type ARR (B-ARR) was up-regulated (Figure 6B). Collectively, the AMF might affect expression of genes related to the transport and transduction of CTK, thereby causing a change in the CTK content.

The content of GA3 increased 87.7% in shoot of apple plants after inoculation with AMF (Supplemental Figure 7). The relative expression levels of almost all genes within the GA synthetic pathway, were decreased, especially GA2-oxidase (GA2ox1: MD05G1207000), which was significantly down-regulated. Genes in the GA transduction pathway did not show any significant differences in their transcript levels (Figure 6C). These data revealed that AMF may influence the content of GA3 by affecting its synthesis and degradation.

In short, the combined analysis of the differential metabolites and differentially expressed genes revealed a significant link between the transcriptomic and metabolomic signatures in apple plants in symbiosis with the AM R. irregularis. These results suggested that the metabolic signature of apple plants was transcriptionally controlled by inoculation with AMF. We used the metabolic and transcriptomic signatures to identify deregulated metabolic pathways. This analysis suggested that constitutive deregulation of sugar metabolism, hormone metabolism, and fatty acid metabolism in the roots of inoculated apple plants affected root growth and that the constitutive deregulation of hormone metabolism in the shoots affected shoot growth after inoculation with AMF. Cell development may also play an important role in response of shoot growth to symbiosis. According to the ratio (root/shoot) of metabolites, the increased shoot growth was a result of not only root development, but also the AMF symbiosis.

Discussion

AMF promotes plant development and growth and can increase the tolerance to stress (Huang et al., 2020; Gao et al., 2020). However, it remains unknown how symbiosis with AMF precisely activates growth. In apple, stress resistance was improved after inoculation with AMF, which also promoted root and shoot development (Huang et al., 2020; Gao et al., 2020). The promotion of shoot and root development after AMF infection was also confirmed in our study (Figure 1). Exploring the molecular mechanism underlying the influence of AMF on plant growth and development will answer many questions about the symbiotic relationship. To determine how AMF promotes both shoot growth and root development, a combined analysis of both the root and shoot of uninoculated and AMF-inoculated apple plants was conducted through transcriptomics and metabolomics.

The trends of the metabolome and the transcriptome were basically compatible with each other, showing that AMF promotes plant development at both levels (Figure 2 and Figure 3). We speculated that AMF mediates changes through the transcriptional control of metabolic pathways. In apple roots, we found that AMF regulation of development may involve multiple pathways, such as sugar, fatty acid, organic acid, and hormone metabolism (Figure 2), all of which may induce root development.

After AMF infection, a large amount of sugars usually accumulate in the root (Beltrano et al., 2013) and Figure 4 and Supplemental Figure 3. At first, most researchers thought that sugars might be the carbon source for the AMF. The fungal symbiotic can take up and use hexose within the root, so the fungus would regulate the accumulation of sugar in the root system to meet its own needs (Solaiman and Saito, 1997). However, Jiang et al. (2017) found that fatty acids from the plant are transferred to the AMF, then converted to Glc by gluconeogenesis in the fungus. Therefore, it remains controversial why AMF regulates the accumulation of sugars in roots. However, accumulation of sugars in roots can affect root development. The accumulation of Fru can repress primary root growth, but induce the development of primary lateral roots, while the accumulation of Suc and Glc promotes the growth of primary and lateral root (Valifard et al., 2021). In our study, Glc, Fru, and Suc were all increased after inoculation with AMF (Figure 4 and Supplemental Figure 3). However, the number of roots increased but the average root length did not after inoculation with AMF (Liu et al., 2018), which is the same as our phenotype (Figure 1). This may be a result of coordination between Glc, Fru, and Suc. Therefore, the accumulation of carbohydrates in roots may be another way for AMF to affect plant root development. Our analysis indicated that AMF increased the accumulation of carbohydrates in roots by regulating carbohydrate synthesis and transport. In our results, AMF promoted the relative expression levels of the alkaline/neutral invertases (NINVs: MD00G1092700) and Trehalase (THL: MD06G1113500), which may increase the content of Glc. The relative expression level of the sugar transporter SWEET (MD10G1269400) was significantly increased (Figure 4), which may affect the accumulation of sugar (Valifard et al., 2021). Fatty acids, as the form of carbon transport from host plants to AMF, also over accumulated in the root system (Jiang et al., 2017; Wu et al., 2019). The content of total fatty acids in the root of mycorrhizae-associated jujube plants was significantly higher non-mycorrhizal jujube plants (Ma et al., 2022). In our results, the relative content of fatty acids (Cis-5,8,11,14,17-Eicosapentaenoic Acid, Arachidonic Acid, and 9-Hydroxy-13-oxo-10-octadecenoic Acid) was also significantly increased (Figure 2. and Supplemental Figure 4). Interestingly, fatty acids also significantly promote root development (Trinh et al., 2019). It seems that AMF regulates the accumulation of fatty acids in roots, which may be another way that AMF affects plant root development.

The relative content of all the differentially regulated organic acids were increased in roots inoculated with mycorrhizae (Figure 2A). Accumulated organic acids can increase the acidity of the rhizosphere (Harker et al., 2002). Acidic conditions are more favorable for apple root growth (Kochian et al., 2015). Although it is common that AMF promotes the accumulation of organic acids in plant roots, not all plant roots are adapted to acidic conditions, such as tomato, rice and barley (Liu et al., 2016; Cristina et al., 2020; Szurman-Zubrzycka et al., 2021). Hage-Ahmed et al. (2013) found that tomato inoculated Glomus mosseae (AMF) had increased levels of protocatechuic acid and salicylic acid, which is in line with our findings. Therefore, we think that the accumulation of organic acids in roots may be a factor of AMF promoting root development of plants adapted to acidic growth conditions, but it may not be applicable to all species.

An enrichment of genes related to the response to hormones indicated the altered hormonal responses in apple plants under inoculation (Figure 3D; Supplemental Figure 5, and Supplemental Data 3). Further analyses revealed that the IAA, CTK, GA and SL responses were altered. Symbiosis between AMF and plants is often accompanied by changes in hormone levels. Auxin is essential for root growth (Yamauchi et al., 2019). The biosynthesis of auxin was promoted following AMF inoculation (Kaldorf and Ludwig-Müller, 2000). In our results, the content of IAA was increased (Figure 5A and Supplemental Figure 5). A previous study showed that YUCCA can convert tryptamine (TAM) to N-hydroxytryptamine (NHT) in vitro (Zhao et al., 2001), suggesting that YUCCA catalyzes the N-oxygenation of TAM, a rate-limiting step in auxin biosynthesis in many plants. The relative expression levels of three YUCCA genes (MD15G1098700, MD00G1056000 and MD05G1184800), ARF18 (MD07G1174000) and ARF2 (MD11G1297900) were upregulated (Figure 5A). This may induce the increased IAA content, which would promote root development.

Likewise, CTK also affects root development, as a higher ratio of AUX/CTK promotes DEFINE (LR) formation (Zhou et al., 2011). In our result, the relative expression of the CTK synthetic gene CYP735A (MD09G1087700) was significantly reduced, which may reduce the content of CTK (tZR) (Figure 5B). The increased IAA and decreased tZR would alter the AUX/CTK ratio of mycorrhizae-associated apple plants, which was higher than the uninoculated apple plants. In summary, the changes in the AUX and CTK contents due to AMF symbiosis are beneficial to root development.

A previous study showed that the content of SL was increased after inoculation with AMF (Guillotin et al., 2017), which is consistent with our results. The relative expression of SL synthesis genes, such as D27 (MD03G1256200, MD11G1277000), CCD8 (MD15G1093500) and LBO (MD16G11720000), were significantly up-regulated, but the relative expression levels of SL signal transduction genes D14 (MD08G1084500) and SMXL6 (MD03G1068500) were significantly down-regulated (Figure 5C). SL promotes primary root development, while the analysis of mutants in SL synthesis or signaling suggest that the absence of SL enhances lateral root formation (Kapulnik et al., 2011; Kumar et al., 2015). These results indicated the AMF may inhibit SL signal transduction to promote root development.

The growth of the aerial parts of the plant is closely related to the growth condition of the underground parts (Ko and Helariutta, 2017). In our study, the ratio of the aerial to underground parts of apple plants did not change significantly before and after AMF infection (Figure 1K). Therefore, the promotion of root development by AMF is likely a key factor affecting shoot growth. On the other hand, a comprehensive analysis of the transcriptomes and metabolomes of roots and shoots showed that both primary and secondary metabolites were quite different between the shoots and roots. The most drastic change was the 38 metabolites that specifically accumulated in the shoot. In the shoots, there were 17 metabolites that specifically accumulated with AMF inoculation, and there were 21 metabolites specific to the uninoculated plants, but no metabolites were specific to the inoculated roots (Table 1 and Supplemental Table 1). Comparing the anatomical structure of the shoot of inoculated and uninoculated plants showed that the shoot cells were significantly longer after inoculation (Supplemental Figure 6A). Therefore, we have full reason to speculate that the promotion of shoot growth by AMF was not merely due to better root development, as AMF also directly regulates shoot growth and development.

AMF clearly regulates aerial growth and development, as the plants grew significantly taller after AMF inoculation (Figure 1B) and AMF regulates some metabolites and genes related to plant growth, such as hormone (IAA, CTK, and GA) synthesis and signaling and cell morphogenesis. IAA, CTK, and GA promote shoot development (Ma et al., 2016; Zheng et al., 2018; Zheng et al., 2022), and AMF promotes stem growth by involving IAA, CTK, GA and other hormones (Dugassa et al., 1996; Barker and Tagu, 2000). Our combined shoot transcriptome and metabolome indicates that AMF regulates hormone synthesis, by affecting transcription related to hormone synthesis in the shoot, such as the genes YUCCA10 (MD00G1056000) and GA2ox (MD05G1207000) (Figure 6A, C), and hormone signal transduction, such as the genes ARF19 (MD15G1221400) and B-ARR (MD16G1159400) (Figures 6A, B). The cell morphology of shoot tissue was significantly changed after inoculation with AMF (Supplemental Figures 6A, B) and the transcriptome showed that genes related to cell morphogenesis were significantly changed (Supplemental Figures 6C). We hypothesized that AMF regulates plant growth by affecting the transcription of morphogenesis-related genes.

Sugar did not seem to play an important role in shoot growth, as neither the metabolites nor the genes changed a lot (Figure 2 and Figure 3). Eighty-two flavonoids showed large changes in response to AMF inoculation (Figure 2, Figure 3; Supplemental Data 2, and Supplemental Data 3), possibly triggered by MYC2 in a JA- and ABA-dependent manner (Adolfsson et al., 2017). It is unknown what roles this change might play in shoot growth in plants in symbiosis with AMF; this might be a meaningful scientific question to study.

Although AMF promoted the growth and development of plant roots and shoots, the colonization rate of AMF varies greatly among species under natural conditions. For instance, the natural colonization rate of soybean was around 40%, and that of wheat was lower than 20% (Gao et al., 2012; Brito et al., 2012). The colonization rate of AMF also was influenced by external environmental factors such as climate, soil, and management patterns (Gosling et al., 2006; Gosling et al., 2010). A high colonization rate is essential for AMF to promote the plants growth and development, and two years of organic management significantly increased the colonization rate of AMF on onions (Gosling et al., 2010), but less research has been studied on other species, so the widespread use of AMF in production still requires extensive field experiments to explore.

Materials and methods

Plant material and experimental treatments

This experiment was carried out in the Northwest A&F University Greenhouse of Shaanxi Province. Plants were derived from cuttings and were maintained in the greenhouse in plastic pots (21 cm ×15 cm) filled with sand (sterilized by steaming at 121°C for 2 h). There were a total of 200 plants of the rootstock M.26. Plants with a height of about 10 cm were selected, of which 100 plants were in the treatment group and 100 plants in the control group. The AMF was Rhizophagus irregularis (strain number: BGC BJ09). The inoculant contained about 60 spores per gram. The treatment group was inoculated with 20 g AMF inoculant, while the control group was inoculated with 20 g sterilized AMF. Place the plants in the greenhouse after treatment (the relative humidity was 55-75%, and the temperatures were 12-30 °C). The plants were watered normally, and all plants were watered with 100 mL ½-strength low-phosphorus Hoagland’s weekly. Sixty days after inoculation with AMF, the plants were measured, and samples were collected.

Mycorrhizal colonization and morphological indices assays

The method used for mycorrhizal colonization observation was described by He et al. (2016), with minor modifications. As follows, the roots were cleaned and then soaked in 10% KOH until transparent and softened. Stained with 0.05% trypan blue, and decolorized with 1:1 actic acid-glycerin for 3 times. Final compression and microscopic examination are performed. At the end of the 60-day co-culture period, plant height, internode length, internode number, root number, and root length were measured. We used an EPSON EXPRESSION 10000XL scanner (LAl600 scanner, Canada) to scan the samples and obtain the total root length. As for the number of roots, we performed a manual count. The net photosynthetic rate was measured by Li6800 (LICOR, Huntington Beach, CA, USA). Selected fully developed mature leaves at the same node to measure net photosynthetic rate, and at least 9 plants were as one biological replicate. At the same time, the chlorophyll content was measured by the acetone method (Chazaux et al., 2022), with the leaves in the same position were selected for determination. The rhizome ratio was calculated using the ratio of root to shoot.

Transcriptome analysis of apple plants in root and shoot

The root tips and the shoot tips were used to extract the RNA, with the protocol of *kit (TIANGEN). All primers used are shown in Supplemental Table 1. A high-quality cDNA library was constructed using. Transcriptome sequence was performed on an Illumina HiSeq platform (Illumina, San Diego, CA, USA). Raw data onto the fast format (original reads) were processed by internal perl script. After program processing, high quality reads were obtained. High-quality reads were mapped to the reference genome (https://www.rosaceae.org/species/malus/all), and the FPKM value of each gene was calculated. Using PCA to assess the relationships of samples. Differential genes were screened on P value <0.05 and fold change ≥2.0. The KO (KEGG orthogonal) and MapMan databases were used to annotate the functions of DEGs.

Metabolome analysis of apple plants in root and shoot

Untargeted metabolomics were used to differentially identify the primary metabolites and secondary metabolites by UPLC - MS/MS analysis. Using PCA to assess the relationships of samples. Differential metabolites were screened on VIP ≥1.0 and fold change ≥2.0. Add 0.1g of sample to 1.2ml of 70% methanol, vortex and shake to mix well. The UPLC system was Thermo Vanquish and column was Accucore Vanquish (Agilent SB-C18 1.8 µm, 2.1 mm * 100 mm). The program was set as follows, with a flow rate of 0.35 mL/min, a column temperature of 40°C, and an injection volume of 4 µL each time. Solvent A was water/0.1% formic acid and solvent B was ACN/0.1% formic acid. For mass spectrometry (MS) analysis, the setup procedure was as follows: the electrospray ionization (ESI) temperature was 550°C, and the GSI, GSIII and curtain gas (CUR) pressures were 50 psi, 60 psi and 25 psi, respectively.

Hormone metabolomics were used to identify differential hormones by LC-MS/MS analysis. Using PCA to assess the relationships of samples. Differential hormone metabolites were screened on P value <0.05 and fold change ≥2.0. The 0.05 g of samples, 10 µl of a mixed internal standard solution (10 ng/ml) and 1 ml methanol/water/formic acid (15:4:1, V/V/V) by vortex to mix well. The UPLC system was Thermo Vanquish and column was Waters ACQUITY UPLC HSS T3 C18 (1.8 µm, 100 mm×2.1 mm i.d.). The program was set as follows, with a flow rate of 0.35 mL/min, a column temperature of 40°C, and an injection volume of 4 µL each time. Solvents A and B were water/0.04% acetic acid and ACN/0.04% acetic acid, respectively. For mass spectrometry (MS) analysis, the Electrospray Ionization (ESI) temperature was 550°C, the Curtain Gas (CUR) pressure was 35 psi.

Statistical analysis of data

Data analysis with SPSS 22.0, the significant differences were determined using Duncan’s test in a one-way ANOVA (p < 0.05).

Data availability statement

The data presented in the study are deposited in the NCBI repository, accession number PRJNA889983.

Author contributions

MZ, ML and FM: conceived and supervised this study; SJ, TY and BL: performed the experiments; YL, LZ and JS: performed the bioinformatics analysis; SJ: wrote the manuscript; ML, BM, MZ and SJ: discussed the study and revised the manuscript. The authors responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors are MZ (mrz@nwsuaf.edu.cn), ML (limingjun@nwsuaf.edu.cn) and FM (fwm64@nwsuaf.edu.cn). All authors contributed to the article and approved the submitted version.

Funding

This research was funded by the Apple Research System (CARS-27).

Acknowledgments

We thank Dr Jing Zhang, Miss Jing Zhao, and Miss Beibei He (Horticulture Science Research Center, Northwest A&F University, Yangling, China) for providing professional technical assistance with UPLC–MS/MS analysis and Microscopy technique.

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/fpls.2022.1052464/full#supplementary-material

References

Adeyemi, N. O., Atayese, M. O., Sakariyawo, O. S., Azeez, J. O., Abayomi Sobowale, S. P., Olubode, A., et al. (2021). Alleviation of heavy metal stress by arbuscular mycorrhizal symbiosis in glycine max (L.) grown in copper, lead and zinc contaminated soils. Rhizosphere 18, 100325. doi: 10.1016/j.rhisph.2021.100325

CrossRef Full Text | Google Scholar

Adolfsson, L., Nziengui, H., Abreu, I. N., Šimura, J., Beebo, A., Herdean, A., et al. (2017). Enhanced secondary- and hormone metabolism in leaves of arbuscular mycorrhizal medicago truncatula. Plant Physiol. 175, 392–411. doi: 10.1104/pp.16.01509

PubMed Abstract | CrossRef Full Text | Google Scholar

Andrino, A., Guggenberger, G., Kernchen, S., Mikutta, R., Sauheitl, L., Boy, J. (2021). Production of organic acids by arbuscular mycorrhizal fungi and their contribution in the mobilization of phosphorus bound to iron oxides. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.661842

PubMed Abstract | CrossRef Full Text | Google Scholar

Barker, S. J., Tagu, D. (2000). The roles of auxins and cytokinins in mycorrhizal symbioses. J. Plant Growth Regul. 19, 144–154. doi: 10.1007/s003440000021

PubMed Abstract | CrossRef Full Text | Google Scholar

Beltrano, J., Ruscitti, M., Arango, M. C., Ronco, M. (2013). Effects of arbuscular mycorrhiza inoculation on plant growth, biological and physiological parameters and mineral nutrition in pepper grown under different salinity and p levels. J. Soil Sci. Plant Nutr. 13 (1), 123–14. doi: 10.4067/S0718-95162013005000012

CrossRef Full Text | Google Scholar

Bonfante, P., Genre, A. (2010). Mechanisms underlying beneficial plant–fungus interactions in mycorrhizal symbiosis. Nat. Commun. 1, 48. doi: 10.1038/ncomms1046

PubMed Abstract | CrossRef Full Text | Google Scholar

Brito, I., Goss, M. J., De Carvalho, M. (2012). Effect of tillage and crop on arbuscular mycorrhiza colonization of winter wheat and triticale under Mediterranean conditions: Effect of tillage and crop on AM. Soil Use Manage. 28, 202–208. doi: 10.1111/j.1475-2743.2012.00404.x

CrossRef Full Text | Google Scholar

Chazaux, M., Schiphorst, C., Lazzari, G., Caffarri, S. (2022). Precise estimation of chlorophyll a , b and carotenoid content by deconvolution of the absorption spectrum and new simultaneous equations for chl determination. Plant J. 109, 1630–1648. doi: 10.1111/tpj.15643

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, W., Meng, P., Feng, H., Wang, C. (2020). Effects of arbuscular mycorrhizal fungi on growth and physiological performance of catalpa bungei C.A.Mey. under drought stress. Forests 11, 1117. doi: 10.3390/f11101117

CrossRef Full Text | Google Scholar

Cristina, G., Camelin, E., Tommasi, T., Fino, D., Pugliese, M. (2020). Anaerobic digestates from sewage sludge used as fertilizer on a poor alkaline sandy soil and on a peat substrate: Effects on tomato plants growth and on soil properties. J. Environ. Manage. 269, 110767. doi: 10.1016/j.jenvman.2020.110767

PubMed Abstract | CrossRef Full Text | Google Scholar

Dugassa, G. D., von Alten, H., Schönbeck, F. (1996). Effects of arbuscular mycorrhiza (AM) on health ofLinum usitatissimum l. infected by fungal pathogens. Plant Soil 185, 173–182. doi: 10.1007/BF02257522

CrossRef Full Text | Google Scholar

Fàbregas, N., Lozano-Elena, F., Blasco-Escámez, D., Tohge, T., Martínez-Andújar, C., Albacete, A., et al. (2018). Overexpression of the vascular brassinosteroid receptor BRL3 confers drought resistance without penalizing plant growth. Nat. Commun. 9, 4680. doi: 10.1038/s41467-018-06861-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, Y., Li, H., Yang, Z. (2002). The ROP2 GTPase controls the formation of cortical fine f-actin and the early phase of directional cell expansion during arabidopsis organogenesis. Plant Cell 14, 777–794. doi: 10.1105/tpc.001537

PubMed Abstract | CrossRef Full Text | Google Scholar

Gadkar, V., David-Schwartz, R., Kunik, T., Kapulnik, Y. (2001). Arbuscular mycorrhizal fungal colonization. factors involved in host recognition. Plant Physiol. 127, 1493–1499. doi: 10.1104/pp.010783

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, T., Liu, X., Shan, L., Wu, Q., Liu, Y., Zhang, Z., et al. (2020). Dopamine and arbuscular mycorrhizal fungi act synergistically to promote apple growth under salt stress. Environ. Exp. Bot. 178, 104159. doi: 10.1016/j.envexpbot.2020.104159

CrossRef Full Text | Google Scholar

Gao, X., Lu, X., Wu, M., Zhang, H., Pan, R., Tian, J., et al. (2012). Co-Inoculation with rhizobia and AMF inhibited soybean red crown rot: From field study to plant defense-related gene expression analysis. PloS One 7, e33977. doi: 10.1371/journal.pone.0033977

PubMed Abstract | CrossRef Full Text | Google Scholar

Gosling, P., Hodge, A., Goodlass, G., Bending, G. D. (2006). Arbuscular mycorrhizal fungi and organic farming. Agriculture Ecosyst. Environ. 113, 17–35. doi: 10.1016/j.agee.2005.09.009

CrossRef Full Text | Google Scholar

Gosling, P., Ozaki, A., Jones, J., Turner, M., Rayns, F., Bending, G. D. (2010). Organic management of tilled agricultural soils results in a rapid increase in colonisation potential and spore populations of arbuscular mycorrhizal fungi. Agriculture Ecosyst. Environ. 139, 273–279. doi: 10.1016/j.agee.2010.08.013

CrossRef Full Text | Google Scholar

Guillotin, B., Etemadi, M., Audran, C., Bouzayen, M., Bécard, G., Combier, J. (2017). Sl-IAA27 regulates strigolactone biosynthesis and mycorrhization in tomato (var. MicroTom ). New Phytol. 213, 1124–1132. doi: 10.1111/nph.14246

PubMed Abstract | CrossRef Full Text | Google Scholar

Hage-Ahmed, K., Moyses, A., Voglgruber, A., Hadacek, F., Steinkellner, S. (2013). Alterations in root exudation of intercropped tomato mediated by the arbuscular mycorrhizal fungus glomus mosseae and the soilborne pathogen fusarium oxysporum f.sp. lycopersici. J. Phytopathol. 161, 763–773. doi: 10.1111/jph.12130

CrossRef Full Text | Google Scholar

Harker, F. R., Marsh, K. B., Young, H., Murray, S. H., Gunson, F. A., Walker, S. B. (2002). Sensory interpretation of instrumental measurements 2: sweet and acid taste of apple fruit. Postharvest Biol. Tec. 24, 241–250. doi: 10.1016/S0925-5214(01)00157-0

CrossRef Full Text | Google Scholar

He, F., Zhang, H., Tang, M. (2016). Aquaporin gene expression and physiological responses of robinia pseudoacacia l. to the mycorrhizal fungus rhizophagus irregularis and drought stress. Mycorrhiza 26, 311–323. doi: 10.1007/s00572-015-0670-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Hosseini, A., Gharghani, A. (2015). Effects of arbuscular mycorrhizal fungi on growth and nutrient uptake of apple rootstocks in calcareous soil. Int. J. Hortic. Sci. Technol. 2, 173–185. doi: 10.22059/ijhst.2015.56434

CrossRef Full Text | Google Scholar

Huang, D., Ma, M., Wang, Q., Zhang, M., Jing, G., Li, C., et al. (2020). Arbuscular mycorrhizal fungi enhanced drought resistance in apple by regulating genes in the MAPK pathway. Plant Physiol. Biochem. 149, 245–255. doi: 10.1016/j.plaphy.2020.02.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Y., Wang, W., Xie, Q., Liu, N., Liu, L., Wang, D., et al. (2017). Plants transfer lipids to sustain colonization by mutualistic mycorrhizal and parasitic fungi. Science 356, 1172–1175. doi: 10.1126/science.aam9970

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaldorf, M., Ludwig-Müller, J. (2000). AM fungi might affect the root morphology of maize by increasing indole-3-butyric acid biosynthesis. Physiologia Plantarum 109 (1), 58-–67. doi: 10.1034/j.1399-3054.2000.100109.x

CrossRef Full Text | Google Scholar

Kapulnik, Y., Delaux, P.-M., Resnick, N., Mayzlish-Gati, E., Wininger, S., Bhattacharya, C., et al. (2011). Strigolactones affect lateral root formation and root-hair elongation in arabidopsis. Planta 233, 209–216. doi: 10.1007/s00425-010-1310-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Karthikeyan, B., Abitha, B., Henry, A. J., Sa, T., Joe, M. M. (2016). “Interaction of rhizobacteria with arbuscular mycorrhizal fungi (AMF) and their role in stress abetment in agriculture,” in Recent advances on mycorrhizal fungi, fungal biology. Ed. Pagano, M. C. (Cham: Springer International Publishing), 117–142. doi: 10.1007/978-3-319-24355-9_11

CrossRef Full Text | Google Scholar

Kochian, L. V., Piñeros, M. A., Liu, J., Magalhaes, J. V. (2015). Plant adaptation to acid soils: The molecular basis for crop aluminum resistance. Annu. Rev. Plant Biol. 66, 571–598. doi: 10.1146/annurev-arplant-043014-114822

PubMed Abstract | CrossRef Full Text | Google Scholar

Ko, D., Helariutta, Y. (2017). Shoot–root communication in flowering plants. Curr. Biol. 27, R973–R978. doi: 10.1016/j.cub.2017.06.054

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, M., Pandya-Kumar, N., Kapulnik, Y., Koltai, H. (2015). Strigolactone signaling in root development and phosphate starvation. Plant Signal Behav. 10, e1045174. doi: 10.1080/15592324.2015.1045174

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, K., Liang, Y., Xing, L., Mao, J., Liu, Z., Dong, F., et al. (2018). Transcriptome analysis reveals multiple hormones, wounding and sugar signaling pathways mediate adventitious root formation in apple rootstock. Int. J. Mol. Sci. 19. doi: 10.3390/ijms19082201

CrossRef Full Text | Google Scholar

Liu, Z., Wang, H., Xu, R. (2016). The effects of root surface charge and nitrogen forms on the adsorption of aluminum ions by the roots of rice with different aluminum tolerances. Plant Soil 408, 43–53. doi: 10.1007/s11104-016-2909-y

CrossRef Full Text | Google Scholar

Liu, C. Y., Wang, P., Zhang, D. J., Zou, Y. N., Kuča, K., Wu, Q. S. (2018). Mycorrhiza-induced change in root hair growth is associated with IAA accumulation and expression of EXPs in trifoliate orange under two p levels. Sci. Hortic. 234, 227–235. doi: 10.1016/j.scienta.2018.02.052

CrossRef Full Text | Google Scholar

Ma, Y., Xue, H., Zhang, L., Zhang, F., Ou, C., Wang, F., et al. (2016). Involvement of auxin and brassinosteroid in dwarfism of autotetraploid apple (Malus × domestica). Sci. Rep. 6, 26719. doi: 10.1038/srep26719

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, Z., Zhao, X., He, A., Cao, Y., Han, Q., Lu, Y., et al. (2022). Mycorrhizal symbiosis reprograms ion fluxes and fatty acid metabolism in wild jujube during salt stress. Plant Physio. 189, 2481–2499. doi: 10.1093/plphys/kiac239

CrossRef Full Text | Google Scholar

Nacoon, S., Jogloy, S., Riddech, N., Mongkolthanaruk, W., Ekprasert, J., Cooper, J., et al. (2021). Combination of arbuscular mycorrhizal fungi and phosphate solubilizing bacteria on growth and production of helianthus tuberosus under field condition. Sci. Rep. 11, 6501. doi: 10.1038/s41598-021-86042-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Niklas, K. J. (2005). Modelling below- and above-ground biomass for non-woody and woody plants. Ann. Botany. 95, 315–321. doi: 10.1093/aob/mci028

CrossRef Full Text | Google Scholar

Solaiman, M., Saito, M. (1997). Use of sugars by intraradical hyphae of arbuscular mycorrhizal fungi revealed by radiorespirometry. New Phytol. 136, 533–538. doi: 10.1046/j.1469-8137.1997.00757.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Szurman-Zubrzycka, M., Chwiałkowska, K., Niemira, M., Kwaśniewski, M., Nawrot, M., Gajecka, M., et al. (2021). Aluminum or low pH – which is the bigger enemy of barley? transcriptome analysis of barley root meristem under Al and low pH stress. Front. Genet. 12. doi: 10.3389/fgene.2021.675260

PubMed Abstract | CrossRef Full Text | Google Scholar

Trinh, D.-C., Lavenus, J., Goh, T., Boutté, Y., Drogue, Q., Vaissayre, V., et al. (2019). PUCHI regulates very long chain fatty acid biosynthesis during lateral root and callus formation. Proc. Natl. Acad. Sci. U.S.A. 116, 14325–14330. doi: 10.1073/pnas.1906300116

PubMed Abstract | CrossRef Full Text | Google Scholar

Valifard, M., Le Hir, R., Müller, J., Scheuring, D., Neuhaus, H. E., Pommerrenig, B. (2021). Vacuolar fructose transporter SWEET17 is critical for root development and drought tolerance. Plant Physiol. 187, 2716–2730. doi: 10.1093/plphys/kiab436

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, N., Jiang, S., Zhang, Z., Fang, H., Xu, H., Wang, Y., et al. (2018). Malus sieversii: the origin, flavonoid synthesis mechanism, and breeding of red-skinned and red-fleshed apples. Hortic. Res. 5, 70. doi: 10.1038/s41438-018-0084-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Zhang, M., Dong, R., Liu, C., Guan, H., Liu, Q., et al. (2021). Heterologous expression of ZmGS5 enhances organ size and seed weight by regulating cell expansion in arabidopsis thaliana. Gene 793, 145749. doi: 10.1016/j.gene.2021.145749

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Q. S., He, J. D., Srivastava, A. K., Zou, Y. N., Kuča, K. (2019). Mycorrhizas enhance drought tolerance of citrus by altering root fatty acid compositions and their saturation levels. Tree Physiol. 39, 1149–1158. doi: 10.1093/treephys/tpz039

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamauchi, T., Tanaka, A., Inahashi, H., Nishizawa, N. K., Tsutsumi, N., Inukai, Y., et al. (2019). Fine control of aerenchyma and lateral root development through AUX/IAA- and ARF-dependent auxin signaling. Proc. Natl. Acad. Sci. U.S.A. 116, 20770–20775. doi: 10.1073/pnas.1907181116

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, C., Shen, S., Zhou, S., Li, Y., Mao, Y., Zhou, J., et al. (2022). Rice metabolic regulatory network spanning the entire life cycle. Mol. Plant 15, 258–275. doi: 10.1016/j.molp.2021.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Yaxley, J. R., Ross, J. J., Sherriff, L. J., Reid, J. B. (2001). Gibberellin biosynthesis mutations and root development in pea. Plant Physiol. 125, 627–633. doi: 10.1104/pp.125.2.627

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y., Christensen, S. K., Fankhauser, C., Cashman, J. R., Cohen, J. D., Weigel, D., et al. (2001). A role for flavin monooxygenase-like enzymes in auxin biosynthesis. Science 291, 306–309. doi: 10.1126/science.291.5502.306

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, X., Li, Y., Ma, C., Chen, B., Sun, Z., Tian, Y., et al. (2022). A mutation in the promoter of the arabinogalactan protein 7-like gene PcAGP7-1 affects cell morphogenesis and brassinolide content in pear ( pyrus communis l.) stems. Plant J. 109, 47–63. doi: 10.1111/tpj.15548

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, X., Zhao, Y., Shan, D., Shi, K., Wang, L., Li, Q., et al. (2018). MdWRKY9 overexpression confers intensive dwarfing in the M26 rootstock of apple by directly inhibiting brassinosteroid synthetase MdDWF4 expression. New Phytol. 217, 1086–1098. doi: 10.1111/nph.14891

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Z.-Y., Zhang, C.-G., Wu, L., Zhang, C.-G., Chai, J., Wang, M., et al. (2011). Functional characterization of the CKRC1/TAA1 gene and dissection of hormonal actions in the arabidopsis root: Functional characterization of the CKRC1/TAA1 gene. Plant J. 66, 516–527. doi: 10.1111/j.1365-313X.2011.04509.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Malus × domestica, mycorrhizal symbiosis, transcriptome, metabolome, root-associated fungus

Citation: Jing S, Li Y, Zhu L, Su J, Yang T, Liu B, Ma B, Ma F, Li M and Zhang M (2022) Transcriptomics and metabolomics reveal effect of arbuscular mycorrhizal fungi on growth and development of apple plants. Front. Plant Sci. 13:1052464. doi: 10.3389/fpls.2022.1052464

Received: 24 September 2022; Accepted: 10 October 2022;
Published: 25 October 2022.

Edited by:

Qiang-Sheng Wu, Yangtze University, China

Reviewed by:

Shunfeng Ge, Shandong Agricultural University, China
Yi Wang, China Agricultural University, China

Copyright © 2022 Jing, Li, Zhu, Su, Yang, Liu, Ma, Ma, Li and Zhang. 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: Mingjun Li, bGltaW5nanVuQG53c3VhZi5lZHUuY24=; Manrang Zhang, bXJ6QG53c3VhZi5lZHUuY24=

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