- Institute of Medicinal Plant Development, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
Dendrobium officinale Kimura et Migo is a traditional and scarce medicinal orchid in China. Mycorrhizal fungi could supply nitrogen (N) to orchids for seed germination and seedling recruitment. However, the N transport mechanism between orchids and the fungus is poorly understand. Early studies found that the fungus MF23 (Mycena sp.) could promote the growth of D. officinale. To better dissect the molecular interactions involved in N transport between D. officinale and MF23, transcriptome and metabolome analyses were conducted on conventional and mycorrhizal cultivations of D. officinale. Moreover, validation tests were carried out in the greenhouse to measure net fluxes of
Introduction
Dendrobium officinale Kimura et Migo (Orchidaceae) is the only plant in Tiepishihu, which has the effects of nourishing Yin, clearing heat, benefiting the stomach, and increasing saliva and has been used as a Chinese traditional medicine and health food for more than 2,000 years (Pharmacopoeia Commission, 2015). The wild resources of D. officinale are on the verge of exhaustion because of its low seed germination rate, low growth tendency, and the excessive demand. More than 90% of the medicinal materials in China depend on facility cultivation, with a planting area of nearly 10,000 ha, mainly distributed in Yunnan and Zhejiang (Si et al., 2017a). The industrial scale of D. officinale has reached 10 billion RMB, and it has been in constant expansion. However, there are some problems in planting, such as the long growth cycle, high cost, and low nutrient conversion rate. Orchid plants generally have a symbiotic relationship with fungi. Numerous studies have shown that mycorrhiza can have positive effects on orchid plant growth. Mycorrhizal fungi can promote tiller growth and increase the total biomass of D. officinale (Tian et al., 2018; Wei et al., 2018; Wang and Hu, 2019; Zhang et al., 2020a). Ceratobasidium sp. AR2 inoculation greatly enhanced the fresh weight, dry weight, plant height, root number, and bud number of Anoectochilus roxburghii (Zhang et al., 2020b). Fungal strains isolated from the roots of natural Cymbidium plants increased the fresh weight, plant height, leaf number, and root length of Cymbidium goeringii (Wu et al., 2010). Therefore, the fungus has substantial potential to solve the planting problem of D. officinale. In our previous studies, the fungus MF23 (Mycena sp.) was shown to promote the growth of D. officinale in different regions, improving the root, stem, and leaf weight of D. officinale (Chen et al., 2016, 2017; Tian et al., 2018). However, the promoting mechanism is still unclear.
Mycorrhizal fungi supply plants with water and nutrients, such as carbon (C), nitrogen (N), and phosphorus (P), or produce signal molecules, such as rhizosphere secretions, fungal secretions, and fungal elicitors, to promote plant growth (Xu and Zhang, 2007). Among these, N is an essential nutrient for plant growth and the fourth most abundant element in plants (Biswas et al., 2021). Limited N availability adversely affects plant growth and causes serious effects on crop production. N is a major nutrient transferred by mycorrhiza, especially in orchids (Stöckel et al., 2014). Mycorrhiza makes the N concentrations in orchid tissues higher than surrounding autotrophic plants (Hynson et al., 2013). Mycorrhizal fungi supply both inorganic and organic N to orchids for seed germination and seedling recruitment (Cameron et al., 2006; Tedersoo et al., 2007; Yang et al., 2018). At present, the N transport model between arbuscular mycorrhizal fungi (AMF) and plants has been established (Govindarajulu et al., 2005; Jin et al., 2005), in this respect, research on N nutrition transport between orchids and mycorrhizal is not as well studied.
To this end, transcriptome and metabolome analyses were conducted to decipher the mechanism of the interaction between D. officinale and MF23. The effects of MF23 on growth, N uptake, and metabolism of D. officinale were systematically studied, which revealed the promoting mechanism of MF23 from the perspective of N nutrition for the first time. Studies of symbiotic D. officinale provide a reliable molecular basis for further investigations of N transport between fungi and orchids. These lay the theoretical foundation for the usage of mycorrhizal fungi in promoting the production of D. officinale and protection of endangered medicinal plants.
Materials and Methods
Fungal Materials
Mycorrhizal fungus MF23 (Mycena sp.) was isolated from the roots of D. officinale, it was preserved by the Mycorrhizal Laboratory of Biotechnology Center, Institute of Medicinal Plant Development (Chen et al., 2016). MF23 was activated and amplified on potato dextrose agar in the dark at 25°C for 20 days. Then, fungal hyphae were transferred to solid medium containing wheat bran and saw dust (1:1, v/v) and were incubated in the dark at 25°C for 30 days to generate a solid inoculum.
Plant Materials
The D. officinale orchids used for transcriptome and metabolome analyses were collected from Jiangsu Yicaotang Zhihu Co. Ltd. (Jiangsu, China). Samples were stems of D. officinale grown in a plastic shed for 15 months after acclimatizing for 1 year [collection time and organ specified in the Pharmacopeia (Pharmacopoeia Commission, 2015)]. The substrate in the seedbed was bark mixed with a small amount of stone. About 2 g of fresh, solid form inocula of MF23 were placed near the seedlings roots when seedlings were transplanted to the seedbed (mycorrhizal cultivation of D. officinale, T; Chen et al., 2016), and seedlings without fungus were conventionally cultivated (CK). For initial inoculum, 1.5 g clump−1 of MF23 inocula was used, and 1.5 g clump−1 of supplementary inocula was added every 6 months.
The D. officinale orchids used for validation tests were collected from a conventional greenhouse (Institute of Medicinal Plant Development, Beijing). Tissue culture seedlings of D. officinale (3–5 cm in height) were cultured on solid medium containing pine bark and sawdust (4:1, v/v). Ten seedlings were transplanted to a solid medium, and 0.5 g of fresh, solid form inocula of MF23 were placed near the seedling roots (dual culture, DC). These conditions were the same for the axenic culture group but without the fungus (AC). The cultures were kept in the greenhouse with a 10 h light/14 h dark photoperiod at 24 ± 1°C and an illumination intensity of 1,500 Lx. Roots of AC and DC were collected at different stages, including 4, 10, and 16 weeks after transplanting.
Growth Index Measurement
For the Omics samples, one clump of D. officinale was used as a repetition. Eighteen repetitions were used for each group. Roots, stems, and leave were collected, dried at 60°C to constant weight, and recorded the dry weight of each organ.
For the validated sample, one culture bottle of D. officinale was used as a repetition. Six repetitions were used for each group. At each time point, the dry weight of roots, stems, and leaves was measured, as well as the highest root length, stem diameter, and seedling height.
Transcriptomic Analysis
Total RNA of the stem samples was extracted using the Improved HUAYUEYANG Quick RNA Isolation Kit (Shan et al., 2018). RNA quality analysis, library construction, sequencing, data filtering and mapping, and differential expression analysis were performed by the Novogene Bioinformatics Technology Co. Ltd. (Beijing, China). Briefly, total RNA degradation and contamination were verified by electrophoresis in a 1.0% agarose gel. RNA concentration and purity were assessed using a NanoDrop™ 2000 Spectrophotometers (Thermo Fisher, United States). Sequencing libraries were generated using the NEBNext® Ultra™ RNA Library Prep Kit for Illumina® (NEB, United States) using 3 μg of total RNA. Clean reads were obtained from the raw data by removing adapter-containing reads, poly-N-containing reads, and low-quality reads. Then, clean reads were mapped to the Dendrobium catenatum reference genome1 using the HISAT2 v.2.0.4 software (Kim et al., 2015; Si et al., 2017b). GOseq software was used for gene ontology (GO) annotation (Young et al., 2010). Differential expression analysis was performed using the DESeq R package (1.18.0) with a cutoff probability of ≥0.8 and │log2Ratio (R2P2CMS/R2P2)│≥1 (Anders and Huber, 2010). Differentially expressed genes (DEGs) were further assigned to the KEGG (Kyoto Encyclopedia of Genes and Genomes) database. KEGG was used to perform pathway analysis. Pathways with a Q-value ≤0.05 were considered significantly enriched pathways (Minoru et al., 2008). Three biological replicates were used for each sample.
Metabolomic Analysis
Chromatography was performed on an UHPLC-LTQ Orbitrap (Thermo Fisher Scientific, Bremen, Germany), using a BEH C18 column (100 mm × 2.1 mm i.d., 1.7 μm; Waters, Milford, United States). The mobile phase consisted of (A) 0.1% formic acid in water and (B) 0.1% formic acid in acetonitrile. They were eluted with the following gradient: 0–1.5 min: 5–25% B, 1.5–10.0 min: 25–100% B, 10.0–13.0 min: 100% B, and 13.0–13.5 min: 100–5, 5% B. The final solvent composition was maintained for 10min (5% B). The flow rate was 0.40 ml/min, and the injection volume was 3 μl; the column oven was set to 45°C. MS acquisition was performed with positive and negative ionization. The electrospray capillary, injection, and collision voltage were 3.0 kV, 40 V, and 30 eV, respectively. The capillary temperature was 350°C, and carrier gas flow was 45 l/h. The mass range was 50–1,000 m/z, and resolution was 30,000.
All plasma samples were detected by LC-MS to obtain plasma metabolic profiles between the groups, and Chroma TOF software was used to analyze peak figures and obtain raw data. The standardized LC-MS data matrix was imported into SIMCA-P+14.0 (Umetrics, Umeå, Sweden) to conduct multivariate statistical analysis, including principal component analysis (PCA) to observe the overall distribution between the samples and the stability of the whole analysis process and (orthogonal) signal correction partial least squares-discriminant analysis (O)PLS-DA to distinguish the overall differences in metabolic profiles between groups and find differential metabolites between groups. To prevent the model from over-fitting, the quality of the model was examined by seven cycles of reciprocal verification and 200 response sequencing tests. Differential metabolites between groups were screened by (O)PLS-DA and a t-test. The metabolites with VIP >1, p < 0.01, and fold change >1.5 were screened as differentially expressed metabolites (DEMs). Finally, plasma differential metabolites between two groups were analyzed by the KEGG pathway database and included channel enrichment and interaction network construction.
Transcriptome and Metabolome Data Set Integration
To assess the potential relevance of quantitative information between mRNA and metabolites, the D. catenatum genome was used to link DEGs and DEMs. R language was used to calculate Pearson’s correlation coefficient. KEGG was used to perform pathway analysis of all correlated DEGs/DEMs.
Measurements of Net Fluxes of and
Plants were N starved for 7 days using 1/2 Murashige and Skoog (MS) medium without NH4NO3 and KNO3 prior to flux analysis using the non-invasive micro-test technology (NMT Physiolyzer® Younger USA LLC, Amherst, MA, United States; Li et al., 2017). Six healthy roots were used for nitrate () and ammonium () flux analyses in each group. The net and fluxes were measured using NMT at the Xuyue (Beijing) Sci. & Tech. Co., Ltd. (Beijing, China).
Briefly, glass micropipettes (Φ4.5 ± 1 μm, XY-CGQ-01) were pre-pulled and salinized in an electrode controller. After finishing preparation, the end of the micropipettes was filled with backfilling solution (, 100 mm NH4Cl; , 10 mm KNO3; Ruan et al., 2016). The remaining portion was filled with 50–80 μm columns of selective liquid ion-exchange cocktails ( LIX, XY-SJ-NH4; LIX, XY-SJ-NO3). An Ag/AgCl wire microsensor holder (YG003-Y11, Younger USA) was inserted in the back of the electrode to contact the electrolyte solution. YG003-Y11 was used as the reference microsensor. The flux microsensor was calibrated at the start and end time of each test with the calibration liquid ( and : 0.05 and 0.5 mm). The flux microsensor with Nernstian slopes greater than 50 mV for and less than −50 mV for per 10-fold concentration difference was used for the following measurements. Before analysis, roots were immersed in a Petri dish that contained 10 ml measuring solution (0.1 mm NH4NO3, 0.1 mm CaCl2, and pH 6.0) and equilibrated for 30 min to lower the impact due to environmental changes.
To determine the area along the root axis corresponding with the maximal net and influx, an initial measurement was made along the root, followed by around 300 μm walk steps, ranging from 300 to 20,000 μm away from the root apex. Then, nine regions were measured as the distance from the root tip: 300, 500, 800, 1,500, 3,000, 5,000, 7,500, 10,000, and 20,000 μm. and fluxes (ca. 5 μm above the root surface) were measured by a reciprocating motion of the ion-selective microelectrode between two positions, in which the vertical movement distance was 30 μm. The ion flux was read every 6 s. The ion flux was recorded at each measurement point for 5 min. Flux data were acquired and analyzed by the imFluxes V2.0 software [imfluxes.com, Xuyue (Beijing) Sci. & Tech. Co., Ltd., Beijing, China].
Validation of RNA-seq Data by RT-qPCR
Total RNA of the root samples was extracted by the method described above. The primers designed with Primer Premier 6.0 are shown in Supplementary Table S1. A PrimeScript™ RT reagent Kit (TaKaRa, Japan) was used for reverse transcription. First, the RT product (1 μl) was diluted with 20 μl ddH2O and used as a template. Then, qPCR was performed in a 15 μl reaction mixture containing 7.5 μl of 2 × SYBR® Premix Ex Taq™ II (TaKaRa, Japan), 1.5 μl of cDNA template, and 0.3 μl of each gene specific primers. Three biological replicates and three technical replicates were conducted using the LightCycler® 480II RT-PCR System (Roche, Switzerland). The parameters for the reactions were 95°C for 30 s, 40 cycles of 95°C for 5 s, and 60°C for 30s. The cDNA libraries were standardized to housekeeping gene 18S rRNA (Zhang et al., 2012). The 2−ΔΔCt method was used for evaluating gene expression.
Determination of N Assimilation Enzyme Activities
Nitrate reductase (NR), nitrite reductase (NiR), glutamine synthetase (GS), glutamate synthase (GOGAT), and glutamate dehydrogenase (GDH) activity in the roots was determined by ELISA, according to the manufacturer’s instructions (ELISA kit, Catalog no. CGM-PMCA1-96, Wuhan Colorful Gene Biological Technology, Wuhan, China).
Statistical Analyses
Data are expressed as the mean ± standard deviation (SD) from at least three independent biological replicates. Significant differences were determined via independent sample t-test and one-way ANOVA (p < 0.05) by the software package SPSS Statistics 17.0 software (SPSS Inc., Chicago, IL, United States).
Results
Sample Information for the Transcriptome and Metabolome
Transcriptome and metabolome analyses were conducted on conventional and mycorrhizal cultivated D. officinale at harvesting time (Figures 1A,B). Mycorrhizal cultivation increased root and leaf dry weight (Figure 1C). The root and leaf dry weight of T was 1.51 ± 0.40 and 0.71 ± 0.18 g, respectively, and was 27.4 and 92.7% higher than CK'.
Figure 1. Sample information for transcriptome and metabolome analyses. (A) Growth conditions. (B) Growth morphology. (C) Growth index. CK: conventional cultivation of Dendrobium officinale and T: mycorrhizal cultivation of D. officinale. *p < 0.05; **p < 0.01.
Transcriptomic Profiles
In this study, transcriptomic analysis was carried out between T and CK to understand the molecular interactions between D. officinale and MF23. In this analysis, 299,481,298 raw reads and 289,048,662 clean reads were obtained (Supplementary Table S2). The sequence data are available in BioProject (Accession number: PRJNA524260) of NCBI. Of the 263 identified DEGs, 94 were upregulated and 169 were downregulated. According to GO gene functional classification system, DEGs were classified to biological process (BP), cellular component (CC), and molecular function (MF) categories and their subcategories. Most DEGs were enriched in BP (e.g., DNA metabolic process, DNA integration, and small GTPase-mediated signal transduction) and MF (e.g., GTP binding, guanyl ribonucleotide binding, and guanyl nucleotide binding; Supplementary Figure S1). Additionally, KEGG pathway analysis revealed 60 enriched pathways, of which, seven pathways exhibited significant differences (p < 0.01; Supplementary Figure S2). The differential pathways were ribosome (dct03010), ribosome biogenesis in eukaryotes (dct03008), 2-oxocarboxylic acid metabolism (dct01210), alanine, aspartate, and glutamate metabolism (dct00250), glycerophospholipid metabolism (dct00564), glyoxylate and dicarboxylate metabolism (dct00630), and pentose phosphate pathway (PPP; dct00030).
Metabolomic Profiles
Principal component analysis was performed and metabolites of T and CK were located on the left and right part, indicating there were differences between groups (Supplementary Figure S3). Differential metabolites were screened by (O)PLS-DA and a student’s t-test. Of the 144 metabolites identified as DEMs, 84 were upregulated and 60 were downregulated (Supplementary Table S3). DEMs were annotated in the KEGG databases. Additionally, KEGG pathway analysis revealed 41 enriched pathways, of which, four pathways exhibited significant differences (p < 0.01; Supplementary Figure S4). The differential pathways were sphingolipid metabolism (dct00600), biosynthesis of amino acids (dct01230), glyoxylate and dicarboxylate metabolism (dct00630), and glutamate metabolism (dct00250).
Association Analysis of the Transcriptome and Metabolome Point to N Metabolism
Both DEGs and DEMs were annotated in KEGG databases. KEGG pathway analysis revealed 39 enriched pathways, of which, nitrogen metabolism (dct00910) and several amino acid metabolisms were involved (Figure 2). These results suggested mycorrhizal cultivation could affect N metabolism of D. officinale. N uptake by plants is actively regulated by numerous N transporters (Nacry et al., 2013). These transporters display different substrate selectivities. Nitrate transporters (NRTs) function as a transporter, and ammonium transporters (AMTs) are responsible for the uptake of (Rennenberg et al., 2010). Transporters that import organic N include amino acid, peptide, and urea transporters (Tegeder and Rentsch, 2010). In the RNA-seq data, nine DEGs of N transporters were found (Table 1). DEGs included two nitrate transporters (NPFs), one nitrate transporter-activating protein (NAR), one amino acid permease (AAP), one lysine/histidine transporter (LHT), one cationic amino acid transporter (CAT), and three oligopeptide transporters (OPTs).
Figure 2. KEGG enrichment analysis of differentially expressed gene (DEGs) and differentially expressed metabolites.
N metabolism is a complicated metabolic process mediated by not just N transporters but N assimilation enzymes (Jose and Andrea, 2016). NR and NiR catalyze the further reduction of to . GS, GOGAT, and GDH allow to be assimilated into glutamine (Gln). In this study, genes encoding NR were significantly downregulated in mycorrhizal cultivation of D. officinale, while genes encoding GDH were significantly upregulated (Table 2). DEMs included Gln (C00064; Supplementary Figure S5), which is the central amino acid linked to assimilation with amino acid biosynthesis. These results indicated that mycorrhizal cultivation could promote assimilation into Gln in D. officinale.
Sample Information for Validation Tests
Transcriptome and metabolome results indicated that mycorrhizal cultivation could influence N and amino acid metabolism in D. officinale. To further study the mechanism of MF23 affecting N metabolism, validation tests were conducted under controllable greenhouse conditions (Supplementary Figures S6, S7). Validation tests showed that MF23 could also increase root and stem dry weight, root length, and seedling height at different growth stages (Figure 3). At week 16, the root and stem dry weight, root length, and seedling height of DC were 0.13 ± 0.01 g, 0.30 ± 0.04 g, 11.12 ± 1.81 cm, and 5.68 ± 0.46 cm, respectively, and were 79.2, 183.5, 71.0, and 27.8% higher than AC', respectively (p < 0.01). With the increase of culture time, the effect of MF23 in increasing root and stem dry weight was more significant. Roots are an important plant organ that have multiple functions, including acquisition and fixation of water and nutrients. The stem is the medicinal part of D. officinale, as specified in the Pharmacopeia (Pharmacopoeia Commission, 2015). Stem dry weight was closely related to yield and economic benefits.
Figure 3. Growth index of samples for validation. (A) Root dry weight. (B) Stem dry weight. (C) Leaf dry weight. (D) Root length. (E) Stem diameter. (F) Seedling height. AC, axenic culture and DC, dual culture. *p < 0.05; **p < 0.01.
Effects of MF23 on and Net Fluxes of D. officinale
NMT is an ion-selective microelectrode technology for the dynamic detection of various ion fluxes of living plants. To determine the effects of MF23 on N uptake in D. officinale, NMT was used to analyze the net fluxes of and in the roots. As shown in Figure 4, at week 4, influxes of AC were higher than , but AC preferred to after that assessment time. The results of were consistent with influxes between AC and DC. and influxes of DC were lower than AC' at week 4, but higher at week 10 and 16. At week 16, the net influxes of and of DC were 35.60 ± 1.49 pmol·cm−2·s−1 and 128.60 ± 9.83 pmol·cm−2·s−1, respectively, which were 44.8 and 154.7% higher, respectively, than that in AC. The amount of was three times that of . In conclusion, with the increase in culture time, the effect of MF23 in promoting and was more significant, especially in promoting the absorption of .
Figure 4. Net influxes of and in D. officinale (n = 6). (A) Detected image. (B) Effects of fungus MF23 on net and influxes in D. officinale (n = 6). AC, axenic culture and DC, dual culture. The Y-axis represents the maximal net and influx. Bars labeled with different letters indicate significant differences between treatments.
Analysis of Transcript Levels of N Transporter DEGs and DEGs Involved in N Metabolism
Of the nine N transporter DEGs, five DEGs involved in N metabolism, and AMT with minimum p value was performed by RT-qPCR analysis (Figure 5). Overall, the expression level of most genes was consistent with the RNA-seq results. Many genes showed high expression in DC, and only DoNPF6.4-1, DoNPF6.4-2, DofmdA, and DoNR had low expression at some stages. These downregulated genes were related to uptake and metabolism. Genes responsible for uptake and assimilation, such as DoAMT11 and DoGDHs, and genes responsible for amino acid and oligopeptide uptake, such as DoAAP3, DoLHT8, DoCAT3, and DoOPTs, were upregulated in DC. DoOPT5 had the highest expression among N transporter genes in DC; the fold change between DC and AC was highest at week 16, which was 10.18 ± 0.35. In contrast to AC, DoGDHs’ expression in DC had higher transcript levels. These results indicated that MF23 promotes and organic N uptake and metabolism.
Figure 5. qRT-PCR analysis of DEGs related to N metabolism in D. officinale (n = 3). (A–J): DEGs of N transporters. (J–O): DEGs involved in N metabolism. The Y-axis represents the mean fold changes in expression values. AC, axenic culture and DC, dual culture. *p < 0.05; **p < 0.01.
Effects of MF23 on Assimilation Enzymes in the N Metabolism of D. officinale
NR, NiR, GS, GOGAT, and GDH were critical rate-limiting enzymes in N metabolism. MF23 increased the activity of GS and GDH in D. officinale. At week 10, GS and GDH activity in DC was 137.92 ± 13.61 and 59.74 ± 5.08 U·L−1, respectively, which was 58.57 and 35.14% higher than that in AC. At week 16, GS and GDH activity in DC was 131.72 ± 9.84 and (58.17 ± 5.15) U·L−1, respectively, which was 23.67 and 31.21% higher than that in AC (p < 0.05; Figure 6). These results showed that MF23 promotes assimilation into Gln and improves the ability to convert inorganic N to organic N by increasing GS and GDH activity.
Figure 6. Effects of fungus MF23 on key enzyme activity of nitrogen assimilation in D. officinale (n = 3). (A) Nitrate reductase (NR) activity. (B) Nitrite reductase (NiR) activity. (C) Glutamine synthetase (GS) activity. (D) Glutamate synthase (GOGAT) activity. (E) Glutamate dehydrogenase (GDH) activity. AC, axenic culture and DC, dual culture. *p < 0.05; **p < 0.01.
Discussion
MF23 Increases N Uptake in D. officinale
Mycorrhizal fungi are able to greatly increase and uptake in the terrestrial orchid C. goeringii (Wu et al., 2013). Recently numerous studies have demonstrated that some NRTs and AMTs are exclusively expressed in the symbiotic system. LjAMT2;2 was exclusively expressed in the mycorrhizal roots of Cynodon dactylon, which is responsible for recruiting and releasing uncharged NH3 to the plant cytoplasm (Guether et al., 2009). NPF was specifically and preferentially activated in the mycorrhizal roots of rice and was responsible for transporting N compounds between plants and fungi (Drechsler et al., 2017). And, the plant AMT SvAMT1 and fungal AMT TcAMT2 were upregulated in the Tulasnella calospora–Serapias vomeracea symbiosis (Fochi et al., 2017a). These studies emphasized the crucial role of fungi for N transporter activation and N transport in plants. Our observations showed that DoNAR2.1 and DoAMT11 were upregulated in the roots of D. officinale under a symbiotic system, and more direct evidence showed a high increase in the net influxes of and in DC. Moreover, MF23 promoted the uptake of and in D. officinale, and similar results were also obtained for the orchid C. goeringii (Wu et al., 2013).
Plants mainly obtain N nutrition by absorbing and , but organic N is the major source transferred to orchid hosts from fungi (Cameron et al., 2006; Kuga et al., 2014; Ghirardo et al., 2020). In addition, SvAAPs, SvLHT, and SvOPTs were more abundant in mycorrhizal protocorms (Fochi et al., 2017b). Our studies also showed that DoAAP3, DoLHT8, DoCAT3, and DoOPTs were upregulated in the roots of D. officinale under a symbiotic system. Among these, DoOPT5 showed the highest expression, which suggested that it may play an important role in N acquisition in the MF23-D. officinale symbiosis. Oligopeptide transport was thought to be the most effective and fastest way for plants to utilize N (Komarova et al., 2008). Therefore, the mechanisms of how MF23 activates the expression of DoOPT5 should be further investigated.
MF23 Increases the Assimilation Capacity in D. officinale
Arbuscular mycorrhizal fungi can absorb and from the soil, but needs to be catalyzed to first to make it more readily available to plants; this conversion consumes energy (Hawkins et al., 2000). Therefore, extraradical mycelium prefers to absorb directly and transfer nitrogen to plants (Leigh et al., 2009). Our studies found that D. officinale preferred to absorb . The influxes of were more than twice the net influxes of at week 16, which was consistent with Zhang et al. (2018), who found that the absorption rate for was higher than that for in epiphytic orchids. More interestingly, the function of MF23 to increase influxes of in D. officinale was significantly higher than the increased influxes of . In this study, after dual culture with MF23 for 16 weeks, the increased influxes of D. officinale was three times that of . Moreover, the DoGDHs were upregulated in RNA-seq data, and the results of RT-qPCR and enzyme activity under symbiotic system were consistent. All results indicated that MF23 increased the assimilation capacity in D. officinale. Based on the results obtained in this study, a schematic map was summarized (Figure 7). MF23 promoted the growth of D. officinale by enhancing N uptake and the assimilation capacity.
Figure 7. Putative nitrogen uptake and metabolism in the D. officinale-MF23 symbiosis. Black font represents the background subjects; red, blue, and gray font represent up-, down-, and undifferentiated-regulated genes or metabolites, respectively. White boxes represent enzymes.
MF23 Influences N Metabolism in D. officinale by Regulating C Metabolism
Symbioses between fungi and plants abound in nature. Host plants transfer excess photosynthates to fungi, and in return, fungal partners uptake and transfer soil-derived nutrients and water to the plant. Nutrient exchanges between AMF and host plants are completed by a mycorrhizal symbiosis (Bonfante and Genre, 2010). Li et al. (2014) suggested the transfer and metabolism of C and N in symbioses may play a critical role in maintaining nutrient balance and resource reallocation between the host plant and fungus. C skeletons and energy provided by host plants are used by AMF to transform and transport N; eventually, these N compounds are returned as C donors to host plants (Li et al., 2014). Fellbaum et al. (2012) found that C availability to the mycorrhizal fungus by the host plant may stimulate N uptake and transfer in the fungal partner (Fellbaum et al., 2012). In non-photosynthetic orchid stages, fungal partners provide amino acids to protocorms for supplying C and N sources, and in return, the plant exports to the fungus (Fochi et al., 2017a). The ornithine cycle of N metabolism is pivotal, linking arginine biosynthesis and C recycling in AMF symbionts (Jin, 2009). These studies suggest that C metabolism is closely associated with N metabolism in mycorrhizal symbioses. In this study, many KEGG significant enrichment pathways were linked to C and N metabolism; for example, alanine, aspartate, and glutamate metabolism belongs to amino acid metabolism. Glyoxylate and dicarboxylate metabolism and the PPP belong to carbohydrate metabolism. Glyoxylate is an intermediate product of photorespiration, which is accompanied by C and N (Younès et al., 2016). PPP is an important pathway for glucose metabolism in plants, and its intermediate product may be involved in the synthesis of amino acids (Huang et al., 2004). Moreover, the results of transcriptome and validation tests showed that DoGDHs were upregulated in the roots of D. officinale under a symbiotic system. GDH plays a pivotal role in plant C and N metabolism (Miyashita and Good, 2008). MF23 promoted the growth of D. officinale by influencing N metabolism associated with C metabolism.
Conclusion
Transcriptome and metabolome analyses were conducted using conventional and mycorrhizal cultivation of D. officinale, which revealed the molecular interactions involved in N metabolism from D. officinale’s response to MF23 infection. The results of validation tests indicated that MF23 promotes the production of D. officinale by increasing N uptake and assimilation capacity. These results also demonstrated the importance of organic N transport in the D. officinale-MF23 symbiosis. The results of this study laid the foundation for further studies of the molecular mechanisms in orchid-fungus symbioses related to N uptake and metabolism.
Data Availability Statement
The original contributions presented in the study are publicly available. This data can be found here: NCBI or DDBJ repositories, accession number PRJNA524260 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA524260).
Author Contributions
SG and XC designed and directed the entire study. TS performed most of the experiments, analyzed the data, and drafted the manuscript. LZ and JL performed the RT-qPCR analysis. XC, LZ, and LT assisted with sampling and data analysis. XC, LZ, and BL edited the manuscript. All authors have read and approved the final manuscript.
Funding
This research was supported by the CAMS Innovation Fund for Medical Sciences (CIFMS; no. 2017-I2M-3-013) and the National Natural Science Foundation of China (nos. 81573526 and 81903749).
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.
Acknowledgments
We would like to thank Wei-dong Zhu from Zhejiang Shouxiangu Pharmaceutical Co., Ltd. for sample collection.
Supplementary Material
The Supplementary Material for this article can be found online at https://www.frontiersin.org/articles/10.3389/fpls.2021.693561/full#supplementary-material
Supplementary Figure S1 | GO classification of DEGs.
Supplementary Figure S2 | KEGG enrichment analysis of DEGs.
Supplementary Figure S3 | Principal component analysis (PCA) of conventional and mycorrhizal cultivation of Dendrobium officinale.
Supplementary Figure S4 | KEGG enrichment analysis of DEMs.
Supplementary Figure S5 | KEGG pathway for N metabolism from the transcriptome and metabolome.
Supplementary Figure S6 | Growth conditions of samples for validation.
Supplementary Table T1 | RT-qPCR primers of N transporter DEGs and DEGs involved in N metabolism.
Supplementary Table T2 | Statistical results of sequencing data.
Supplementary Table T3 | Differentially expressed metabolites.
Footnotes
References
Anders, S., and Huber, W. (2010). Differential expression analysis for sequence count data. Genome Biol. 11:R106. doi: 10.1186/gb-2010-11-10-r106
Biswas, S. S., Singh, D. R., De, L. C., Kalaivanan, N. S., Pal, R., and Janakiram, T. (2021). A comprehensive scenario of orchid nutrition – a review. J. Plant Nutr. 44, 905–917. doi: 10.1080/01904167.2021.1871758
Bonfante, P., and Genre, A. (2010). Mechanisms underlying beneficial plant–fungus interactions in mycorrhizal symbiosis. Nat. Commun. 1:48. doi: 10.1038/ncomms1046
Cameron, D. D., Leake, J. R., and Read, D. J. (2006). Mutualistic mycorrhiza in orchids: evidence from plant-fungus carbon and nitrogen transfers in the green-leaved terrestrial orchid Goodyera repens. New Phytol. 171, 405–416. doi: 10.1111/j.1469-8137.2006.01767.x
Chen, X. M., Yan, H. L., Tian, L. X., Wang, C. L., Wang, A. R., and Guo, S. X. (2017). Mycorrhizal cultivation of Dendrobium officinale. Chin. Pharm. J. 52, 1120–1125. doi: 10.11669/cpj.2017.13.005
Chen, X. M., Yan, H. L., Wang, C. L., Tian, L. X., Wang, A. R., and Guo, S. X. (2016). Effects of mycorrhizal fungus Mycena sp. on the growth and polysaccharide properties of Dendrobium officinale. Sci. China Life Sci. 59, 974–976. doi: 10.1007/s11427-015-0367-0
Drechsler, N., Courty, P. E., Brulé, D., and Kunze, R. (2017). Identification of arbuscular mycorrhiza-inducible nitrate transporter 1/peptide transporter family (NPF) genes in rice. Mycorrhiza 28, 93–100. doi: 10.1007/s00572-017-0802-z
Fellbaum, C. R., Gachomo, E. W., Beesetty, Y., Choudhari, S., Strahan, G. D., Pfeffer, P. E., et al. (2012). Carbon availability triggers fungal nitrogen uptake and transport in arbuscular mycorrhizal symbiosis. Proc. Natl. Acad. Sci. U. S. A. 109, 2666–2671. doi: 10.1073/pnas.1118650109
Fochi, V., Chitarra, W., Kohler, A., Voyron, S., Singan, V. R., Lindquist, E. A., et al. (2017a). Fungal and plant gene expression in the Tulasnella calospora-Serapias vomeracea symbiosis provides clues about nitrogen pathways in orchid mycorrhizas. New Phytol. 213, 365–379. doi: 10.1111/nph.14279
Fochi, V., Falla, N., Girlanda, M., Perotto, S., and Balestrini, R. (2017b). Cell-specific expression of plant nutrient transporter genes in orchid mycorrhizae. Plant Sci. 263, 39–45. doi: 10.1016/j.plantsci.2017.06.015
Ghirardo, A., Fochi, V., Lange, B., Witting, M., Schnitzler, J.-P., Perotto, S., et al. (2020). Metabolomic adjustments in the orchid mycorrhizal fungus Tulasnella calospora during symbiosis with Serapias vomeracea. New Phytol. 228, 1939–1952. doi: 10.1111/nph.16812
Govindarajulu, M., Pfeffer, P. E., Jin, H., Abubaker, J., Douds, D. D., Allen, J. W., et al. (2005). Nitrogen transfer in the arbuscular mycorrhizal symbiosis. Nature 435, 819–823. doi: 10.1038/nature03610
Guether, M., Neuhäuser, B., Balestrini, R., Dynowski, M., Ludewig, U., and Bonfante, P. (2009). A mycorrhizal-specific ammonium transporter from Lotus japonicus acquires nitrogen released by arbuscular mycorrhizal fungi. Plant Physiol. 150, 73–83. doi: 10.1104/pp.109.136390
Hawkins, H. J., Johansen, A., and George, E. (2000). Uptake and transport of organic and inorganic nitrogen by arbuscular mycorrhizal fungi. Plant Soil 226, 275–285. doi: 10.1023/A:1026500810385
Huang, J., Wang, J. F., and Zhang, H. S. (2004). Advances on plant pentose phosphate pathway and its key enzymes. Chin. Bull. Bot. 21, 139–145. doi: 10.3969/j.issn.1674-3466.2004.02.002
Hynson, N. A., Madsen, T. P., Selosse, M. A., Adam, I., and Gebauer, G. (2013). The Physiological Ecology of Mycoheterotrophy. New York: Springer.
Jin, H. R. (2009). Arginine bi-directional translocation and breakdown into ornithine along the arbuscular mycorrhizal mycelium. Sci. China C Life Sci. 52, 381–389. doi: 10.1007/s11427-009-0044-1
Jin, H. R., Pfeffer, P. E., Douds, D. D., Piotrowski, E., Lammers, P. J., and Shachar-Hill, Y. (2005). The uptake, metabolism, transport and transfer of nitrogen in an arbuscular mycorrhizal symbiosis. New Phytol. 168, 687–696. doi: 10.1111/j.1469-8137.2005.01536.x
Jose, O. B., and Andrea, V. (2016). Nitrate transport, sensing, and responses in plants. Mol. Plant 9, 837–856. doi: 10.1016/j.molp.2016.05.004
Kim, D., Langmead, B., and Salzberg, S. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357–360. doi: 10.1038/nmeth.3317
Komarova, N. Y., Thor, K., Gubler, A., Meier, S., Dietrich, D., Weichert, A., et al. (2008). AtPTR1 and AtPTR5 transport dipeptides in planta. Plant Physiol. 148, 856–869. doi: 10.1104/pp.108.123844
Kuga, Y., Sakamoto, N., and Yurimoto, H. (2014). Stable isotope cellular imaging reveals that both live and degenerating fungal pelotons transfer carbon and nitrogen to orchid protocorms. New Phytol. 202, 594–605. doi: 10.1111/nph.12700
Leigh, J., Hodge, A., and Fitter, A. H. (2009). Arbuscular mycorrhizal fungi can transfer substantial amounts of nitrogen to their host plant from organic material. New Phytol. 181, 199–207. doi: 10.1111/j.1469-8137.2008.02630.x
Li, Q., Ding, G., Li, B., and Guo, S. X. (2017). Transcriptome analysis of genes involved in dendrobine biosynthesis in Dendrobium nobile Lindl. Infected with mycorrhizal fungus MF23 (Mycena sp.). Sci. Rep. 7:316. doi: 10.1038/s41598-017-00445-9
Li, Y. J., Liu, Z. L., He, X. Y., and Tian, C. J. (2014). Metabolism and interaction of C and N in the arbuscular mycorrhizal symbiosis. J. Appl. Ecol. 25, 903–910.
Minoru, K., Michihiro, A., Susumu, G., Masahiro, H., Mika, H., Masumi, I., et al. (2008). KEGG for linking genomes to life and the environment. Nucleic Acids Res. 36, D480–D484. doi: 10.1093/nar/gkm882
Miyashita, Y., and Good, A. G. (2008). NAD(H)-dependent glutamate dehydrogenase is essential for the survival of Arabidopsis thaliana during dark-induced carbon starvation. J. Exp. Bot. 59, 667–680. doi: 10.1093/jxb/erm340
Nacry, P., Bouguyon, E., and Gojon, A. (2013). Nitrogen acquisition by roots: physiological and developmental mechanisms ensuring plant adaptation to a fluctuating resource. Plant Soil 370, 1–29. doi: 10.1007/s11104-013-1645-9
Pharmacopoeia Commission (2015). Pharmacopoeia of People's Republic of China. Vol. 1. Beijing: China Medical Science Press.
Rennenberg, H., Wildhagen, H., Ehlting, B., Polle, A., and Douglas, C. (2010). Nitrogen nutrition of poplar trees. Plant Biol. 12, 275–291. doi: 10.1111/j.1438-8677.2009.00309.x
Ruan, L., Wei, K., Wang, L., Cheng, H., Zhang, F., Wu, L., et al. (2016). Characteristics of and fluxes in tea (Camellia sinensis) roots measured by scanning ion-selective electrode technique. Sci. Rep. 6:38370. doi: 10.1038/srep38370
Shan, T. T., Chen, X. M., Guo, S. X., Wang, Q. Q., and Wang, A. R. (2018). Studies on total RNA extraction methods from the stems of Dendrobium species. Biotechnol. Bull. 34, 54–58. doi: 10.13560/j.cnki.biotech.bull.1985.2017-0920
Si, J. P., Wang, P., Liu, Z. J., Liu, J. J., and Luo, Y. B. (2017a). Breakthrough in key science and technologies in Dendrobium catenatum industry. China J. Chin. Mater. Med. 42, 2223–2227. doi: 10.19540/j.cnki.cjcmm.2017.0102
Si, J. P., Zhang, Y., Luo, Y. B., Liu, J. J., and Liu, Z. J. (2017b). Herbal textual research on relationship between Chinese medicine “Shihu” (Dendrobium spp.) and “Tiepi Shihu” (D. catenatum). China J. Chin. Mater. Med. 42, 2001–2005. doi: 10.19540/j.cnki.cjcmm.20170415.002
Stöckel, M., Těšitelová, T., Jersáková, J., Bidartondo, M. I., and Gebauer, G. (2014). Carbon and nitrogen gain during the growth of orchid seedlings in nature. New Phytol. 202, 606–615. doi: 10.1111/nph.12688
Tedersoo, L., Pellet, P., Kõljalg, U., and Selosse, M.-A. (2007). Parallel evolutionary paths to mycoheterotrophy in understorey ericaceae and orchidaceae: ecological evidence for mixotrophy in pyroleae. Oecologia 151, 206–217. doi: 10.1007/s00442-006-0581-2
Tegeder, M., and Rentsch, D. (2010). Uptake and partitioning of amino acids and peptides. Mol. Plant 3, 997–1011. doi: 10.1093/mp/ssq047
Tian, L. X., Chen, X. M., Guo, S. X., and Shan, T. T. (2018). Agronomic traits, chemical compositions and response ability to mycorrhizal fungi of two Dendrobium officinale phenotypes. Chin. Pharm. J. 53, 98–103. doi: 10.11669/cpj.2018.02.005
Wang, Q. X., and Hu, H. (2019). Effects of mycorrhizal symbiosis on growth and development of Dendrobium officinale. Shandong Agri. Sci. 51, 61–66. doi: 10.14083/j.issn.1001-4942.2019.01.012
Wei, M., Yu, M. Y., and Chai, R. J. (2018). Effect of orchid mycorrhizal fungi on photosynthetic characteristics and pepc gene expression of Dendrobium officinale. Acta Bot. Boreal-Occid. Sin. 38, 2276–2283. doi: 10.7606/j.issn.1000-4025.2018.12.2276
Wu, J. R., Ma, H. C., Lü, M., Han, S., Zhu, Y., Jin, H., et al. (2010). Rhizoctonia fungi enhance the growth of the endangered orchid Cymbidium goeringii. Botany-Botanique 88, 20–29. doi: 10.1139/B09-092
Wu, J. R., Ma, H. C., Xu, X. L., Qiao, N., Guo, S. T., Liu, F., et al. (2013). Mycorrhizas alter nitrogen acquisition by the terrestrial orchid Cymbidium goeringii. Ann. Bot. 111, 1181–1187. doi: 10.1093/aob/mct062
Xu, H., and Zhang, J. (2007). Effects of arbuscular mycorrhizal fungi on plant growth. Bull. Bot. Res. 27, 636–640. doi: 10.7525/j.issn.1673-5102.2007.05.029
Yang, J. W., Ling, H., Zhang, Y., Zeng, X., and Guo, S. X. (2018). Effects of endophytic fungi on seed germination of medicinal plants of Orchidaceae: a review. Mycosystema 37, 22–34. doi: 10.13346/j.mycosystema.170198
Younès, D., Mathieu, J., Jessica, S., Maurino, V. G., and Michael, H. (2016). Photorespiratory glycolate–glyoxylate metabolism. J. Exp. Bot. 67, 3041–3052. doi: 10.1093/jxb/erw090
Young, M. D., Wakefield, M. J., Smyth, G. K., and Oshlack, A. (2010). Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 11:R14. doi: 10.1186/gb-2010-11-2-r14
Zhang, Y., Li, Y. Y., Chen, X. M., Guo, S. X., and Lee, Y. I. (2020a). Effect of different mycobionts on symbiotic germination and seedling growth of Dendrobium officinale, an important medicinal orchid. Bot. Stud. 61:2. doi: 10.1186/s40529-019-0278-6
Zhang, Y., Li, Y. Y., and Guo, S. X. (2020b). Effects of the mycorrhizal fungus Ceratobasidium sp. AR2 on growth and flavonoid accumulation in Anoectochilus roxburghii. PeerJ 8:e8346. doi: 10.7717/peerj.8346
Zhang, G., Song, C., Zhao, M. M., Li, B., and Guo, S. X. (2012). Characterization of an A-type cyclin-dependent kinase gene from Dendrobium candidum. Biologia 67, 360–368. doi: 10.2478/s11756-012-0016-y
Keywords: Dendrobium officinale, Mycena, nitrogen uptake, metabolome, glutamate dehydrogenase
Citation: Shan T, Zhou L, Li B, Chen X, Guo S, Wang A, Tian L and Liu J (2021) The Plant Growth-Promoting Fungus MF23 (Mycena sp.) Increases Production of Dendrobium officinale (Orchidaceae) by Affecting Nitrogen Uptake and
Edited by:
Shangguo Feng, Hangzhou Normal University, ChinaReviewed by:
Xiaoyu Ding, Nanjing Normal University, ChinaShao Shicheng, Xishuangbanna Tropical Botanical Garden (CAS), China
Xiaohua Jin, Institute of Botany, Chinese Academy of Sciences, China
Copyright © 2021 Shan, Zhou, Li, Chen, Guo, Wang, Tian and Liu. 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: Xiaomei Chen, Y3htX2ltcGxhZEAxNjMuY29t; Shunxing Guo, c3hndW8xOTg2QDE2My5jb20=
†These authors have contributed equally to this work