Abstract
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 and of roots by a non-invasive micro-test technology (NMT), determine N assimilation enzyme activity by the ELISA, and analyze the expression level of differentially expressed genes (DEGs) of N transporters and DEGs involved in N metabolism by RT-qPCR. Combined transcriptome and metabolome analyses showed that MF23 may influence N metabolism in D. officinale. The expression of DoNAR2.1 (nitrate transporter-activating protein), DoAMT11 (ammonium transporter), DoATFs (amino acid transporters), DoOPTs (oligopeptide transporters), and DoGDHs (glutamate dehydrogenases) in symbiotic D. officinale was upregulated. NMT results showed a preference for in D. officinale and indicated that MF23 could promote the uptake of and, especially for . ELISA results showed that MF23 could increase the activity of glutamine synthetase (GS) and glutamate dehydrogenase. This study suggested that MF23 increases the production of D. officinale by affecting N uptake and assimilation capacity.
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 (). 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 (). 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 (; ; ; ). Ceratobasidium sp. AR2 inoculation greatly enhanced the fresh weight, dry weight, plant height, root number, and bud number of Anoectochilus roxburghii (). Fungal strains isolated from the roots of natural Cymbidium plants increased the fresh weight, plant height, leaf number, and root length of Cymbidium goeringii (). 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 (, ; ). 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 (). Among these, N is an essential nutrient for plant growth and the fourth most abundant element in plants (). 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 (). Mycorrhiza makes the N concentrations in orchid tissues higher than surrounding autotrophic plants (). Mycorrhizal fungi supply both inorganic and organic N to orchids for seed germination and seedling recruitment (; ; ). At present, the N transport model between arbuscular mycorrhizal fungi (AMF) and plants has been established (; ), 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 (). 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 ()]. 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; ), 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 (). 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 (; ). GOseq software was used for gene ontology (GO) annotation (). 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 (). 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 (). 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; ). 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; ). 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 (). 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
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 (). These transporters display different substrate selectivities. Nitrate transporters (NRTs) function as a transporter, and ammonium transporters (AMTs) are responsible for the uptake of (). Transporters that import organic N include amino acid, peptide, and urea transporters (). 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
Table 1
| Gene family | Subfamily | Gene ID | Description | Log2FC (T/CK) | P |
|---|---|---|---|---|---|
| NRTs: nitrate transporters | NPF (or NRT1/PTR): nitrate Transporter1/peptide transporter | 110109825 | Protein NRT1/PTR FAMILY 6.4-like | −3.31 | 9.86E-05 |
| 110107420 | Protein NRT1/ PTR FAMILY 6.4 | −2.87 | 8.10E-03 | ||
| NAR: nitrate transporter-activating protein | 110104483 | High-affinity nitrate transporter-activating protein 2.1-like isoform X1 | 1.95 | 2.90E-06 | |
| ATF (or AAAP): amino acid transporter family (or amino acid/auxin permease) | AAP: amino acid permease | 110111144 | Amino acid permease 3-like | −2.40 | 4.00E-03 |
| LHT: lysine/histidine transporter | 110100062 | Lysine histidine transporter-like 8 | −1.24 | 4.00E-03 | |
| APC: amino acid-polyamine-choline | CAT: cationic amino acid transporter | 110111677 | Cationic amino acid transporter 3, mitochondrial-like | −0.91 | 7.10E-03 |
| OPT: oligopeptide transporter | OPT: oligopeptide transporter | 110101252 | Oligopeptide transporter 5-like isoform X1 | 2.80 | 4.56E10 |
| 110098435 | Oligopeptide transporter 1-like isoform X1 | 0.82 | 1.46E-02 | ||
| 110115169 | Oligopeptide transporter 4-like | 0.65 | 4.93E-02 |
DEGs of N transporters.
N metabolism is a complicated metabolic process mediated by not just N transporters but N assimilation enzymes (). 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.
Table 2
| Description | Enzyme nomenclature | Gene ID | Log2FC (T/C) | P |
|---|---|---|---|---|
| fmdA, formamidase | [EC:3.5.1.49] | 110115292 | −2.71 | 4.80E-02 |
| NR [NAD(P)H], nitrate reductase | [EC:1.7.1.1 1.7.1.2 1.7.1.3] | 110114643 | −1.16 | 4.63E-02 |
| GDH (NADP+), glutamate dehydrogenase | [EC:1.4.1.4] | 110091953 | 1.50 | 2.58E-02 |
| 110091954 | 1.64 | 2.57E-02 | ||
| GDH[NAD(P)+] | [EC:1.4.1.3] | 110108713 | 4.07 | 4.78E-18 |
DEGs involved in N metabolism.
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 (). Stem dry weight was closely related to yield and economic benefits.
Figure 3
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
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
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
Discussion
MF23 Increases N Uptake in D. officinale
Mycorrhizal fungi are able to greatly increase and uptake in the terrestrial orchid C. goeringii (). 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 (). NPF was specifically and preferentially activated in the mycorrhizal roots of rice and was responsible for transporting N compounds between plants and fungi (). And, the plant AMT SvAMT1 and fungal AMT TcAMT2 were upregulated in the Tulasnella calospora–Serapias vomeracea symbiosis (). 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 ().
Plants mainly obtain N nutrition by absorbing and , but organic N is the major source transferred to orchid hosts from fungi (; ; ). In addition, SvAAPs, SvLHT, and SvOPTs were more abundant in mycorrhizal protocorms (). 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 (). 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 (). Therefore, extraradical mycelium prefers to absorb directly and transfer nitrogen to plants (). 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 , 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
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 (). 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 (). found that C availability to the mycorrhizal fungus by the host plant may stimulate N uptake and transfer in the fungal partner (). 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 (). The ornithine cycle of N metabolism is pivotal, linking arginine biosynthesis and C recycling in AMF symbionts (). 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 (). PPP is an important pathway for glucose metabolism in plants, and its intermediate product may be involved in the synthesis of amino acids (). 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 (). 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.
Statements
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).
Acknowledgments
We would like to thank Wei-dong Zhu from Zhejiang Shouxiangu Pharmaceutical Co., Ltd. for sample collection.
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.
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 S1GO classification of DEGs.
Supplementary Figure S2KEGG enrichment analysis of DEGs.
Supplementary Figure S3Principal component analysis (PCA) of conventional and mycorrhizal cultivation of Dendrobium officinale.
Supplementary Figure S4KEGG enrichment analysis of DEMs.
Supplementary Figure S5KEGG pathway for N metabolism from the transcriptome and metabolome.
Supplementary Figure S6Growth conditions of samples for validation.
Supplementary Table T1RT-qPCR primers of N transporter DEGs and DEGs involved in N metabolism.
Supplementary Table T2Statistical results of sequencing data.
Supplementary Table T3Differentially expressed metabolites.
References
1
AndersS.HuberW. (2010). Differential expression analysis for sequence count data. Genome Biol.11:R106. doi: 10.1186/gb-2010-11-10-r106
2
BiswasS. S.SinghD. R.DeL. C.KalaivananN. S.PalR.JanakiramT. (2021). A comprehensive scenario of orchid nutrition – a review. J. Plant Nutr.44, 905–917. doi: 10.1080/01904167.2021.1871758
3
BonfanteP.GenreA. (2010). Mechanisms underlying beneficial plant–fungus interactions in mycorrhizal symbiosis. Nat. Commun.1:48. doi: 10.1038/ncomms1046
4
CameronD. D.LeakeJ. R.ReadD. 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
5
ChenX. M.YanH. L.TianL. X.WangC. L.WangA. R.GuoS. X. (2017). Mycorrhizal cultivation of Dendrobium officinale. Chin. Pharm. J.52, 1120–1125. doi: 10.11669/cpj.2017.13.005
6
ChenX. M.YanH. L.WangC. L.TianL. X.WangA. R.GuoS. 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
7
DrechslerN.CourtyP. E.BruléD.KunzeR. (2017). Identification of arbuscular mycorrhiza-inducible nitrate transporter 1/peptide transporter family (NPF) genes in rice. Mycorrhiza28, 93–100. doi: 10.1007/s00572-017-0802-z
8
FellbaumC. R.GachomoE. W.BeesettyY.ChoudhariS.StrahanG. D.PfefferP. 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
9
FochiV.ChitarraW.KohlerA.VoyronS.SinganV. R.LindquistE. 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
10
FochiV.FallaN.GirlandaM.PerottoS.BalestriniR. (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
11
GhirardoA.FochiV.LangeB.WittingM.SchnitzlerJ.-P.PerottoS.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
12
GovindarajuluM.PfefferP. E.JinH.AbubakerJ.DoudsD. D.AllenJ. W.et al. (2005). Nitrogen transfer in the arbuscular mycorrhizal symbiosis. Nature435, 819–823. doi: 10.1038/nature03610
13
GuetherM.NeuhäuserB.BalestriniR.DynowskiM.LudewigU.BonfanteP. (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
14
HawkinsH. J.JohansenA.GeorgeE. (2000). Uptake and transport of organic and inorganic nitrogen by arbuscular mycorrhizal fungi. Plant Soil226, 275–285. doi: 10.1023/A:1026500810385
15
HuangJ.WangJ. F.ZhangH. 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
16
HynsonN. A.MadsenT. P.SelosseM. A.AdamI.GebauerG. (2013). The Physiological Ecology of Mycoheterotrophy.New York: Springer.
17
JinH. 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
18
JinH. R.PfefferP. E.DoudsD. D.PiotrowskiE.LammersP. J.Shachar-HillY. (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
19
JoseO. B.AndreaV. (2016). Nitrate transport, sensing, and responses in plants. Mol. Plant9, 837–856. doi: 10.1016/j.molp.2016.05.004
20
KimD.LangmeadB.SalzbergS. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods12, 357–360. doi: 10.1038/nmeth.3317
21
KomarovaN. Y.ThorK.GublerA.MeierS.DietrichD.WeichertA.et al. (2008). AtPTR1 and AtPTR5 transport dipeptides in planta. Plant Physiol.148, 856–869. doi: 10.1104/pp.108.123844
22
KugaY.SakamotoN.YurimotoH. (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
23
LeighJ.HodgeA.FitterA. 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
24
LiQ.DingG.LiB.GuoS. 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
25
LiY. J.LiuZ. L.HeX. Y.TianC. J. (2014). Metabolism and interaction of C and N in the arbuscular mycorrhizal symbiosis. J. Appl. Ecol.25, 903–910. PMID:
26
MinoruK.MichihiroA.SusumuG.MasahiroH.MikaH.MasumiI.et al. (2008). KEGG for linking genomes to life and the environment. Nucleic Acids Res.36, D480–D484. doi: 10.1093/nar/gkm882
27
MiyashitaY.GoodA. 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
28
NacryP.BouguyonE.GojonA. (2013). Nitrogen acquisition by roots: physiological and developmental mechanisms ensuring plant adaptation to a fluctuating resource. Plant Soil370, 1–29. doi: 10.1007/s11104-013-1645-9
29
Pharmacopoeia Commission (2015). Pharmacopoeia of People's Republic of China.Vol. 1. Beijing: China Medical Science Press.
30
RennenbergH.WildhagenH.EhltingB.PolleA.DouglasC. (2010). Nitrogen nutrition of poplar trees. Plant Biol.12, 275–291. doi: 10.1111/j.1438-8677.2009.00309.x
31
RuanL.WeiK.WangL.ChengH.ZhangF.WuL.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
32
ShanT. T.ChenX. M.GuoS. X.WangQ. Q.WangA. 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
33
SiJ. P.WangP.LiuZ. J.LiuJ. J.LuoY. 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
34
SiJ. P.ZhangY.LuoY. B.LiuJ. J.LiuZ. 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
35
StöckelM.TěšitelováT.JersákováJ.BidartondoM. I.GebauerG. (2014). Carbon and nitrogen gain during the growth of orchid seedlings in nature. New Phytol.202, 606–615. doi: 10.1111/nph.12688
36
TedersooL.PelletP.KõljalgU.SelosseM.-A. (2007). Parallel evolutionary paths to mycoheterotrophy in understorey ericaceae and orchidaceae: ecological evidence for mixotrophy in pyroleae. Oecologia151, 206–217. doi: 10.1007/s00442-006-0581-2
37
TegederM.RentschD. (2010). Uptake and partitioning of amino acids and peptides. Mol. Plant3, 997–1011. doi: 10.1093/mp/ssq047
38
TianL. X.ChenX. M.GuoS. X.ShanT. 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
39
WangQ. X.HuH. (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
40
WeiM.YuM. Y.ChaiR. 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
41
WuJ. R.MaH. C.LüM.HanS.ZhuY.JinH.et al. (2010). Rhizoctonia fungi enhance the growth of the endangered orchid Cymbidium goeringii. Botany-Botanique88, 20–29. doi: 10.1139/B09-092
42
WuJ. R.MaH. C.XuX. L.QiaoN.GuoS. T.LiuF.et al. (2013). Mycorrhizas alter nitrogen acquisition by the terrestrial orchid Cymbidium goeringii. Ann. Bot.111, 1181–1187. doi: 10.1093/aob/mct062
43
XuH.ZhangJ. (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
44
YangJ. W.LingH.ZhangY.ZengX.GuoS. X. (2018). Effects of endophytic fungi on seed germination of medicinal plants of Orchidaceae: a review. Mycosystema37, 22–34. doi: 10.13346/j.mycosystema.170198
45
YounèsD.MathieuJ.JessicaS.MaurinoV. G.MichaelH. (2016). Photorespiratory glycolate–glyoxylate metabolism. J. Exp. Bot.67, 3041–3052. doi: 10.1093/jxb/erw090
46
YoungM. D.WakefieldM. J.SmythG. K.OshlackA. (2010). Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol.11:R14. doi: 10.1186/gb-2010-11-2-r14
47
ZhangY.LiY. Y.ChenX. M.GuoS. X.LeeY. 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
48
ZhangY.LiY. Y.GuoS. X. (2020b). Effects of the mycorrhizal fungus Ceratobasidium sp. AR2 on growth and flavonoid accumulation in Anoectochilus roxburghii. PeerJ8:e8346. doi: 10.7717/peerj.8346
49
ZhangG.SongC.ZhaoM. M.LiB.GuoS. X. (2012). Characterization of an A-type cyclin-dependent kinase gene from Dendrobium candidum. Biologia67, 360–368. doi: 10.2478/s11756-012-0016-y
50
ZhangS. B.YangY. J.LiJ. W.QinJ.ZhangW.HuangW.et al. (2018). Physiological diversity of orchids. Plant Diversity40, 196–208. doi: 10.1016/j.pld.2018.06.003
Summary
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 Assimilation. Front. Plant Sci. 12:693561. doi: 10.3389/fpls.2021.693561
Received
27 April 2021
Accepted
11 June 2021
Published
15 July 2021
Volume
12 - 2021
Edited by
Shangguo Feng, Hangzhou Normal University, China
Reviewed by
Xiaoyu Ding, Nanjing Normal University, China; Shao Shicheng, Xishuangbanna Tropical Botanical Garden (CAS), China; Xiaohua Jin, Institute of Botany, Chinese Academy of Sciences, China
Updates
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, cxm_implad@163.com; Shunxing Guo, sxguo1986@163.com
†These authors have contributed equally to this work
This article was submitted to Plant Symbiotic Interactions, a section of the journal Frontiers in Plant Science
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.