Comparative Transcriptomics and Proteomics of Atractylodes lancea in Response to Endophytic Fungus Gilmaniella sp. AL12 Reveals Regulation in Plant Metabolism

The fungal endophyte Gilmaniella sp. AL12 can establish a beneficial association with the medicinal herb Atractylodes lancea, and improve plant growth and sesquiterpenoids accumulation, which is termed “double promotion.” Our previous studies have uncovered the underling primary mechanism based on some physiological evidences. However, a global understanding of gene or protein expression regulation in primary and secondary metabolism and related regulatory processes is still lacking. In this study, we employed transcriptomics and proteomics of Gilmaniella sp. AL12-inoculated and Gilmaniella sp. AL12-free plants to study the impact of endophyte inoculation at the transcriptional and translational levels. The results showed that plant genes involved in plant immunity and signaling were suppressed, similar to the plant response caused by some endophytic fungi and biotroph pathogen. The downregulated plant immunity may contribute to plant-endophyte beneficial interaction. Additionally, genes and proteins related to primary metabolism (carbon fixation, carbohydrate metabolism, and energy metabolism) tended to be upregulated after Gilmaniella sp. AL12 inoculation, which was consistent with our previous physiological evidences. And, Gilmaniella sp. AL12 upregulated genes involved in terpene skeleton biosynthesis, and upregulated genes annotated as β-farnesene synthase and β-caryophyllene synthase. Based on the above results, we proposed that endophyte-plant associations may improve production (biomass and sesquiterpenoids accumulation) by increasing the source (photosynthesis), expanding the sink (glycolysis and tricarboxylic acid cycle), and enhancing the metabolic flux (sesquiterpenoids biosynthesis pathway) in A. lancea. And, this study will help to further clarify plant-endophyte interactions.

Introduction

Medicinal plants are rich in active compounds such as artemisinin (Sharma and Agrawal, 2013) and ginseng saponin (Wei et al., 2018), which are important sources of modern drugs. Within the increasing population pressure, costs and side effects of drugs, demands for the uses of medicines from plants for treatment of various human ailments are increasing. Compared with synthetic medicines, compounds from plants are thought to be safe to human beings and the ecosystem (Nema et al., 2013). To meet the market demand for these medicinal materials, artificial cultivation technology has been employed for planting medicinal herbs in some Asian countries (Zhou et al., 2005; Sharma et al., 2010; Oh et al., 2014; Lv et al., 2017). However, it is difficult to guarantee the quality and quantity of artificially cultivated medicinal plants because they are more vulnerable to infections by pests and pathogens. Thus, it is imperative to seek effective methods for medicinal plant cultivation. In the past few years, endophytic fungi possessing plant growth promoting properties have been an effective tool for medicinal plant cultivation (Mandal et al., 2013; Ming et al., 2013; Sharma and Agrawal, 2013; Zheng et al., 2016; Zhai et al., 2018).

Atractylodes lancea (Thunb.) DC., belonging to the Asteraceae family, is an endangered traditional Chinese medicinal herb (Wang et al., 2008). Its bioactive component, the sesquiterpenoids, possesses various pharmacology properties such as antibacterial, antitumour, and immunomodulation abilities (Wang et al., 2008; Koonrungsesomboon et al., 2014; Na-Bangchang et al., 2017). Over the past few years, natural sources of A. lancea have been in short supply because of the excessive exploitation and slow growth rate of the herb (Zhou et al., 2016). The medicinal source of A. lancea mainly derives from artificial cultivation, but the yield and quality of this herb are relatively low (Zhou et al., 2016). At present, it is urgent to improve the quality and quantity of the herb as the market demand for A. lancea is increasing on a daily basis. The endophytic fungus Gilmaniella sp. AL12 isolated from stem of A. lancea can establish a beneficial interaction with the host plant (Wang et al., 2012) and promote plant growth and sesquiterpenoid accumulation of tissue culture seedings, which is termed the “double promotion” effect of the endophyte on A. lancea (Yuan et al., 2016b). Consistent with this phenomenon, the endophytic fungi AL12 promotes plant growth and sesquiterpenoid accumulation within two years of growth in field experiments. Therefore, a beneficial interaction of Gilmaniella sp. AL12 with A. lancea is considered suitable for cultivation of A. lancea and will provide a theoretical reference for endophytic fungi-medicinal herb interactions.

In view of the limited carbon and energy source in plants, the accumulation of secondary metabolites occurs at the cost of primary metabolism, representing a discrepancy with the “double promotion” effect of Gilmaniella sp. AL12 on A. lancea. The plant growth promotion effect of the endophyte on A. lancea has been preliminarily ascribed to nutrient assimilation, photosynthesis, and phytohormone content regulation (Yuan et al., 2016b). Moreover, the enhanced sesquiterpenoids accumulation of A. lancea has been shown to be mediated by multiple defense related signals of the host induced by the endophyte (Wang et al., 2011; Ren and Dai, 2012, 2013; Yuan et al., 2016a). Given that primary metabolism-dependent terpenoid precursor biosynthesis and secondary metabolism-related terpenoid skeleton biosynthesis and transformation are simultaneously involved in sesquiterpenoid synthesis (Dudareva et al., 2006; Chen W. et al., 2017; Sharma et al., 2017; Vattekkatte et al., 2018), the molecular and biochemical regulation of the plants relevant to primary and secondary metabolism should be considered. However, thus far, a global understanding of the -regulated expression of genes or proteins in primary and secondary metabolism and related regulatory processes is still lacking.

In this study, we employed transcriptomics and proteomics on endophyte-inoculated and endophyte free plants to better understand the impact of Gilmaniella sp. AL12 on plant metabolism and related regulatory processes of A. lancea at the transcriptional and translational level. The following four essential questions were addressed in this study: (1) Which plant metabolic or regulatory processes of A. lancea are affected by Gilmaniella sp. AL12? (2) What is the effect of the fungal endophyte on the regulation of primary metabolism-dependent terpenoid precursor biosynthesis in A. lancea? (3) What is the effect of the fungal endophyte on the regulation of secondary metabolism-related terpenoid skeleton biosynthesis and transformation in A. lancea? (4) Are other potential signals involved in the fungal endophyte-induced sesquiterpenoids biosynthesis of A. lancea?

Materials and Methods

Plantlet Material and Fungal Inoculation

Atractylodes lancea meristem cultures were established using sterilized plantlets according to our previous studies (Wang et al., 2012). Firstly, meristem cultures were established using mature A. lancea planted in Maoshan, Jiangsu Province, China (Wang et al., 2012). Sterile adventitious buds (approximately 2–3 cm long) of young stems were collected and carefully washed under running tap water. They were surface sterilized by immersing in ethanol (75%) for 30 s, followed by soaking in mercury chloride solution (1%) for 10 min and rinsing in sterile distilled water five times (Wang et al., 2012). Subsequent procedures were conducted aseptically (Wang et al., 2012). The explants were transferred in 50 mL Murashige and Skoog medium containing sucrose (30 g L⁻¹), agar (10 g L⁻¹), naphthaleneacetic acid (0.3 mg L⁻¹), and 6-benzyladenine (2.0 mg L⁻¹) in 150-mL sealed Erlenmeyer flask to emerge adventitious buds for 4 weeks (Wang et al., 2012). Then, newborn adventitious buds were separated and grown in 50 mL Murashige and Skoog medium containing sucrose (30 g L⁻¹), agar (10 g L⁻¹), naphthaleneacetic acid (0.3 mg L⁻¹) and 6-benzyladenine (2.0 mg L⁻¹) in 150-mL sealed Erlenmeyer flask for further differentiation until newborn adventitious buds were suficient for transplantation. After separation, newborn axillary buds were transplanted into 50 mL Murashige and Skoog medium containing sucrose (30 g L⁻¹), agar (10 g L⁻¹), and naphthaleneacetic acid (0.25 mg L⁻¹) in a 150-mL sealed Erlenmeyer flask to develop into rooting plantlets. Each bud was cultured in an Erlenmeyer flask. All plantlets were grown in growth chamber (PGX-600A-12H, Ningbo Lifewww Technology Co., China) with a photoperiod of 12 h, a light density of 3,400 lm m⁻², and a temperature cycle of 25/18°C day/night. Thirty-day-old rooting plantlets were chosen for the endophytic fungus inoculation experiments. The endophytic fungus Gilmaniella sp. AL12 isolated from the stem of A. lancea was grown on potato dextrose agar medium (Chen et al., 2008; Wang et al., 2012). After five days of culture at 28°C, 5-mm Gilmaniella sp. AL12 mycelial disks were placed near the plant caudexes on the medium for inoculation. Additionally, equal-sized potato dextrose agar disks were used as the control (Yuan et al., 2016a). All plantlets were grown in growth chamber and performed in triplicate.

Plant Growth and Sesquiterpenoid Content Analysis

Plantlets were harvested at 15 DPI for the shoot dry weight and sesquiterpenoid content analysis. Briefly, harvested plants were dried at 30°C until the weight was constant, and then 1 g of the ground powder was extracted with 4 mL cyclohexane at room temperature for 10 h. After 15 min of sonication (60 Hz), the mixture was centrifuged for 5 min at 5,000 g at 4°C. After filtering through 0.22-μm-diameter microporous membranes, the total cyclohexane extract was dried over anhydrous sodium sulfate and stored in a dark glass bottle at 4°C before gas chromatography (GC) analysis (Yuan et al., 2016b).

A GC system (Agilent 7890A, United States) equipped with a fame ionization detector was used for sesquiterpenoid analysis (Yuan et al., 2016b). The GC column was DB-1ms (30 m × 0.32 mm × 0.10 μm) with high-purity nitrogen as a carrier gas at a flow rate of 0.8 mL min⁻¹. For sesquiterpenoid analysis, 1 μl of cyclohexane was directly injected onto the GC column with the following temperature program: injection at 240°C, an initial temperature of 100°C (4 min hold) raised to 140°C at a rate of 10°C min⁻¹ (10-min hold), followed by an increase to 220°C at a rate of 10°C min⁻¹ (10-min hold), raised to 260°C at a rate of 10°C min⁻¹ (2-min hold). The detector temperature was then set at 350°C. As previously described, qualitative and quantitative analyses of seven sesquiterpenoids were performed according to authentic standards (Yuan et al., 2016b).

Determination of Photosynthetic Parameters

Plantlets were harvested at 10: 00 AM at 15 DPI for Chlorophyll a fluorescence measurement using a Handy-PEA instrument (Handy PEA, Hansatech, United Kingdom) (Wang et al., 2017). Prior to measurement, A. lances leaves were dark-adapted for 30 min at room temperature to relax the reaction centers. The leaves were exposed to red light of 650 nm at an excitation irradiance of 3,000 μmol m⁻² s⁻¹ for 800 ms. Then, fluorescence parameters were calculated according to Supplementary Table S1.

RNA Sequencing and Functional Annotation of Differentially Expressed Genes

Plantlets were harvested at 15 DPI. Shoot tissue of endophyte-inoculated or endophyte-free plants were harvested for RNA extraction using TRIzol reagent (Invitrogen, CA, United States) (Chen F. et al., 2017). After quality and purity checking, equal amounts of RNA samples from three biological replicates of CK (endophyte-free plants for control) and AL12 (endophyte-inoculated plants) were used for construction of six mRNA-Seq libraries. Libraries were indexed, pooled and then sequenced on an Illumina HiSeq^TM 4000 sequencing platform (Vazyme Biotech Co., Nanjing, China). After filtering, the remaining high-quality reads were assembled de novo by using the Trinity program. Certain short reads with overlapping regions were assembled into longer clusters (contigs). Paired-end reads were used to fill scaffolds gaps to obtain unigenes. Based on a BLASTX search (E value < 10⁻⁵), functional annotation of these unigenes was performed with the following public databases, including the National Center for Biotechnology Information (NCBI) non-redundant protein (NR) and nucleotide (NT) database, Swiss-Prot protein database, Gene Ontology (GO) database, Cluster of Orthologous Groups (COGs) database, and Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The fragments per kilobase of transcript per million mapped reads (FPKM) method was used to determine the unigene expression. Differentially expressed genes (DEGs) between CK and AL12 samples were identified with log₁₀ (AL12_FPKM/CK_FPKM) ≥ 1 and a false discovery rate (FDR) ≤ 0.05. Subsequently, upregulated and downregulated DEGs were conducted by GO and KEGG enrichment analysis, respectively. The RNA sequencing data for this article were submitted to the NCBI Sequence Read Archive (SRA) under accession number SRP132616.

Real-Time Quantitative PCR Validation

Real-time quantitative PCR (RT-qPCR) was performed to determine the expression levels of differentially expressed genes of A. lancea cultured in vitro. Elongation factor 1 alpha gene (EF1a) was used as an internal reference (Yuan et al., 2016b). Primers for the selected DEGs are listed in Supplementary Table S2. One microgram of total RNA was transcribed into cDNA using HisCRIPT^® II Q RT SuperMix for qPCR (+gDNA wiper) (Vazyme Biotech Co., Nanjing, China) according to the manufacturer’s instruction. Subsequently, RT-qPCR was conducted using the DNA Engine Opticon 2 Real-time PCR Detection System (Bio-Rad, Hercules, CA, United States). The reaction system consisted of a volume of 20 μL, which included 10 μL of 2 × AceQ qPCR SYBR^® Green Master Mix (High ROX Premixed) (Vazyme Biotech Co., Nanjing, China), 2 μL of the cDNA template, 0.4 μL of each primer, and 7.6 μL ddH₂O. The reaction conditions were 95°C for 5 min, followed by 40 cycles of 95°C for 30 s, 60°C for 30 s, 95°C for 15 s, 60°C for 60 s, and 95°C for 15 s. All assays were performed in triplicate. Relative expression levels for each cDNA sample were calculated by the 2^−ΔΔCt method (Wang et al., 2015; Liu et al., 2018).

Total Protein Extraction, Two-Dimensional Gel Electrophoresis (2-DE) Separation, Image Analysis, In-Gel Digestion, Protein Identification and Functional Analysis

Plant shoots were harvested at 15 DPI for protein extraction. Protein was extracted using the trichloroacetic acid/acetone precipitation method (Wang et al., 2016). For this purpose, 2 g of shoot material was frozen in liquid nitrogen, ground into a fine powder, mixed with 10 mL cold trichloroacetic acid/acetone buffer (13% (w/v) trichloroacetic acid, 0.07% (v/v) 2-mercaptoethanol in acetone, and kept overnight at −20°C. After shaking for 1 h, the proteins were separated by centrifugation (10,000 g, 15 min, 4°C), and washed twice with cold acetone and once with 80% (v/v) acetone. The proteins were then air-dried and dissolved in protein lysis buffer (7 M urea, 2 M thiourea, 4% (w/v) 3-[(3-cholamidopropyl) dimethylam-monio]-1-propanesulphonate (CHAPS), 2% (v/v) IPG buffer, 1% (v/v) phenylmethanesulphonyl fluoride (PMSF), 10 mM Na₃VO₄, 10 mM NaF, 50 mM glycerophosphate, 5 μg mL⁻¹ antiprotease, 5 μg mL⁻¹ trasylol and 5 μg mL⁻¹ leupeptin). After sonication and centrifugation, the resulting protein extraction was quantified using a commercial dye reagent (Bio-Rad Laboratories, Hercules, CA, United States).

Proteins (900 μg) were separated using immobilized pH gradient (IPG) strips (pH 3–10, nonlinear gradient, 24 cm) and rehydrated for 12 h. First-dimension isoelectric focusing was performed at 20°C for 1 h at 250 V, 2 h at 1,000 V, 5 h at 10,000 V, and 12 h at 500 V (Chen et al., 2016). The strips were then equilibrated with a buffer containing 6 M urea, 2% (w/v) sodium dodecyl sulphate (SDS), 75 mM Tris-HCl (pH 8.8), 30% (v/v) glycerol, and 1% (w/v) dithiothreitol (DTT). In the second dimension, 12% (w/v) sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) was used, and then coomassie brilliant blue stained (GE Healthcare). Gel images were scanned using an Image Scanner III (GE Healthcare) and analyzed using the Image Mater 2D Platinum 7.0 (GE Healthcare). Each spot was normalized to a relative volume. Quantitative analysis sets were created between each control and AL12-inoculated plants, and each treatment was performed in triplicate.

Differentially expressed spots were selected, manually excised, digested with trypsin (Tatli et al., 2017), and analyzed by UltrafleXtreme TOF/TOF (Bruker Daltonics, Germany). Reflector positive mode was used, with a wavelength of 355 nm and acceleration voltage of 2 kV. MS and MS/MS data were analyzed, and peak lists were generated using FlexAnalysis 3.1 (Bruker Daltonics, Germany). MS and MS/MS analyses were compared to the NCBI green plant protein database using the MASCOT 2.4 search engine (Matrix Science, London, United Kingdom) (Chen et al., 2016; Tatli et al., 2017). Search parameters were as follows: trypsin digestion with one missed cleavage, variable modifications (oxidation of methionine and carbamidomethylation of cysteine), and mass tolerance of a precursor ion and fragment ion at 0.2 Da for +1 charged ions. For all proteins successfully identified by a Peptide Mass Fingerprint and/or MS/MS, a Mascot score greater than 50 (the default MASCOT threshold for such searches) was accepted as significant (P < 0.05) (Tatli et al., 2017).

Statistical Analysis

The mean values and standard deviations were calculated using SPSS Statistics 17.0 software (SPSS Inc., Chicago, United States), and the statistical evaluation between two treatments was compared using the independent-samples t-test. Values are the means of three independent experiments. Bars represent standard deviations. Asterisks denote significant differences from the control (t-test; ^∗P < 0.05; ^∗∗P < 0.01).

Results

Gilmaniella sp. AL12 Improved Shoot Growth and Sesquiterpenoid Accumulation in A. lancea

The content of β-caryophyllene, zingiberene, β-sesquiphellandrene, caryophyllene oxide, and atractylone in shoots of A. lancea all significantly increased after fungal endophyte inoculation (Figure 1A). The shoot dry weight of endophyte-inoculated plants was higher than that of endophyte-free A. lancea (Figure 1B). And, we analyzed chlorophyll fluorescence parameter of endophyte-inoculated and endophyte-free A. lancea (Figures 1C,D). The active photosynthesis II (PSII) reaction center captures of light energy to transform into excitation energy, then converts part of excitation energy into chemical energy, thus promotes carbon assimilation reaction, and the rest of excitation energy dissipates (Li et al., 2005). Although the maximum quantum yield of primary photochemistry (F_v/F_m) of A. lancea was unchanged after endophyte inoculation (Figure 1C), the endophyte increased phenomenological fluxes per cross section including absorption flux of photons (ABS/CS_m), phenomenological fluxes for trapping (TR_o/CS_m), and potential electron transport (ET_o/CS_m) of A. lancea (Figure 1D), indicating that the fungal endophyte improved PS II reaction center performance in A. lancea.

FIGURE 1

Figure 1. The effect of the fungal endophyte Gilmaniella sp. AL12 inoculation on shoot sesquiterpenoid content (A), shoot dry weight (B), chlorophyll fluorescence relative value (C), and phenomenological fluxes per cross section (D) of A. lancea. Thirty-day-old plantlets treated with 5-mm AL12 mycelial disks were harvested at 15 DPI. Controls were established using equal sized potato dextrose agar disks. Values are the means of three independent experiments. Bars represent standard deviations. Asterisks denote significant differences from the control (t-test; ^∗P < 0.05; ^∗∗P < 0.01).

Gilmaniella sp. AL12 Induced Shoot Transcriptional Changes in A. lancea

We analyzed the transcriptome of A. lancea shoots inoculated or not inoculated with endophytic fungi AL12. A total of 1956 upregulated and 2063 downregulated genes were identified (Supplementary Figure S1). To further elaborate the function of the DEGs, we first conducted GO enrichment analysis (Supplementary Figure S2). The most upregulated DEGs belonged to oxidation-reduction-related GO terms. Among the downregulated DEGs, stress response-related GO terms were the most enriched. We further performed KEGG pathway enrichment analysis for the upregulated and downregulated DEGs (Figure 2). For the upregulated DEGs, “metabolic pathway” and “biosynthesis of secondary metabolism” were the most enriched KEGG pathways (Figure 2A), while “plant-pathogen interaction” and “plant hormone signal transduction” were significantly enriched pathways for the downregulated DEGs (Figure 2B).

FIGURE 2

Figure 2. KEGG analysis of differentially expressed genes in shoots of A. lancea after the fungal endophyte Gilmaniella sp. AL12 inoculation. (A) The top 21 enriched KEGG pathways of AL12 upregulated and (B) AL12 downregulated genes in shoots of A. lancea.

The DEGs were functionally related to known metabolic pathways, secondary metabolic pathways, and regulatory processes using the online software iPath2: interactive Pathways Explorer¹. And, pathways including transcription, translation, replication and repair, protein folding, and transport tend to be upregulated, while protein degradation pathway tend to be downregulated in shoots of A. lancea after the fungal endophyte inoculation (Supplementary Table S3). Compared with non-inoculated plants, in endophyte-inoculated plants, nucleotide metabolism, glycan biosynthesis and metabolism, and lipid metabolism pathways tends to be downregulated (Supplementary Table S4). Additionally, the endophytic fungus AL12 upregulated the energy, carbohydrate, amino acid, cofactor and vitamin, and secondary metabolism pathways in shoots of A. lancea (Supplementary Tables S5–S7).

Gilmaniella sp. AL12 Upregulated Genes Involved in Primary Metabolism

Our results revealed that most DEGs annotated as functioning in photosynthesis and oxidative phosphorylation were globally upregulated after AL12 inoculation (Supplementary Table S5 and Figure 3), including genes annotated as ferredoxin-NADP⁺ reductase (FNR), photosystem II oxygen-evolving enhancer protein 2 (OEE), ribulose-1,5-bisphosphate carboxylase/oxygenase activase (RCA), NADH dehydrogenase, and V-type H⁺-transporting ATPase, among others. However, three genes encoding chlorophyll A/B binding protein (LHCP) were down-regulated.

FIGURE 3

Figure 3. Differentially expression genes involved in photosynthesis, carbon/nitrogen metabolism, and oxidative phosphorylation in shoots of A. lancea after the fungal endophyte Gilmaniella sp. AL12 inoculation. LHCP, light-harvesting chlorophyll a/b-binding protein; OEE, PSII oxygen-evolving enhancer protein; FNR, Ferredoxin-NADP reductase; RCA, ribulose bisphosphate carboxylase/oxygenase activase; RPLA, Ribose 5-phosphate isomerase; FBPA, fructose-bisphosphate aldolase; SS, starch synthase; SBE, 1,4-alpha-glucan branching enzyme; amylase, PHO, starch phosphorylase; TPT, Glucose-6-phosphate/phosphate-translocator; SUS, sucrose synthase; invertase, β-fructosidases; FRK, fructokinase; GlcCerase, Beta-glucosidase-related glycosidases; GAPDH, glyceraldehyde 3-phosphate dehydrogenase; PEPCK, phosphoenolpyruvate carboxykinase (ATP); PK, pyruvate kinase; PPDK, Pyruvate orthophosphate dikinase; PDHE, pyruvate dehydrogenase; LDHD, lactate dehydrogenase; PDC, pyruvate decarboxylase; ALDH, aldehyde dehydrogenase (NAD⁺); ADH, alcohol dehydrogenase; MDH, malate dehydrogenase; NRT1, low-affinity nitrate transporter; NRT2, high-affinity nitrate transporter; NR, nitrate reductase; GS, glutamine synthetase; GOGAT, glutamate synthase; GDH, glutamate dehydrogenase; AST, asparagine transaminase; ASNS, asparagine synthetase; sugar, bidirectional sugar transporter; ALMT, aluminum-activated malate transporter; SLAH, S-type anion channel; MATE, MATE efflux family protein; CCR, cytochrome-c reductase; COX, Cytochrome c oxidase; PPases, Inorganic pyrophosphatase.

Endophyte inoculation enhanced plant carbon metabolism, with an upregulation of genes involved in carbohydrate metabolism and transport (Figure 3). For instance, most DEGs annotated as phosphate transporter (TPT), amylase, invertase, β-glucosidase (GlcCerase), fructokinase (FRK), glyceraldehyde 3-phosphate dehydrogenase (GAPDH), phosphopyruvate hydratase (enolase), and phosphoenolpyruvate carboxykinase (PEPCK) were upregulated. Interestingly, genes that function in controlling pyruvate and acetyl-CoA biosynthesis were upregulated. For example, AL12 upregulated genes encoding lactate dehydrogenase (LDHD), pyruvate decarboxylase (PDC), and aldehyde dehydrogenase (ALDH), suggesting that it promotes the conversion of lactate into pyruvate and the biosynthesis of acetic acid. Similarly, the expression of genes encoding alcohol dehydrogenase (ADH) and malate dehydrogenase (MDH) was increased after AL12 inoculation, indicating that oxaloacetate (OAA) and acetic acid biosynthesis were improved. Moreover, the endophytic fungus Gilmaniella sp. AL12 upregulated genes encoding MDH and pyruvate dehydrogenase (PDHE), thus contributing to acetyl-CoA biosynthesis and OAA regeneration. Consistent with the improved tricarboxylic acid cycle (TCA cycle), seven genes annotated as malate transporter were up-regulated.

The balance of carbon and nitrogen metabolism is very important for plant metabolism. Most DEGs annotated as functioning in nitrogen metabolism were upregulated after AL12 inoculation (Figure 3). These included genes encoding high affinity nitrate transporter (NRT2), glutamate dehydrogenase (GDH), and others. In particular, one gene encoding glutamine synthetase (GS) exhibited a 10-fold increase in expression in response AL12 inoculation. Additionally, most DEGs annotated as functioning in amino acid metabolism were upregulated (Supplementary Table S6).

Most DEGs annotated as functioning in lipid metabolism were downregulated after AL12 treatment (Supplementary Table S4). These DEGs compass fatty acid elongation, the biosynthesis of unsaturated fatty acids, wax biosynthesis, linoleic acid metabolism, alpha-linolenic acid metabolism, and others. In contrast, most genes involved in fatty acid metabolism, steroid biosynthesis, glycerophospholipid metabolism, ether lipid metabolism, and sphingolipid metabolism were upregulated. Among the upregulated genes, it is worth noting that 32 DEGs encoding phospholipase D (PLD) were upregulated after AL12 inoculation, and eight showed a greater than 10-fold increase in expression (Supplementary Table S8).

Gilmaniella sp. AL12 Upregulated Genes Involved in Secondary Metabolism

A total of 299 DEGs involved in 18 biosynthesis pathways of secondary metabolites were regulated in A. lancea after AL12 inoculation (Figure 4). Among these pathways, DEGs involved in eight biosynthesis pathways of secondary metabolites were significantly regulated in A. lancea after AL12 inoculation (Supplementary Table S9). The eight pathways include Phenylpropanoid biosynthesis, Zeatin biosynthesis, Flavonoid biosynthesis, Stilbenoid, diarylheptanoid and gingerol biosynthesis, Limonene and pinene degradation, Flavone and flavonol biosynthesis, Carotenoid biosynthesis, and Diterpenoid biosynthesis. The largest subcategory was phenylpropanoid biosynthesis, followed by zeatin biosynthesis. The upregulated DEGs involved in phenylpropanoid biosynthesis were associated with lignin biosynthesis. Additionally, most of the DEGs involved in zeatin biosynthesis were upregulated, particularly adenylate dimethylallyl transferase (IPT). Consistent with the increased sesquiterpenoids content in A. lancea shoots (Figure 1A), the expression of genes associated with terpenoid backbone biosynthesis and sesquiterpene biosynthesis was increased (Figure 4 and Supplementary Figure S3). Genes encoding farnesyl diphosphate synthase (FDPS) and β-farnesene synthase were increased after AL12 treatment (Supplementary Figures S3, S4).

FIGURE 4

Figure 4. Distribution of differentially expression genes involved in secondary metabolism in shoots of A. lancea after the fungal endophyte Gilmaniella sp. AL12 inoculation. The x axis indicates the KEGG pathway. The y axis presents the number of DEGs. Numbers in red block or green blocks represent the number of upregulated or downregulated DEGs, respectively.

Gilmaniella sp. AL12 Downregulated Genes Involved in Plant-Pathogen Interactions

Most DEGs associated with plant-pathogen interaction were downregulated after AL12 inoculation (Figure 5). Among the genes related to pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI), DEGs annotated as leucine-rich repeat (LRR) receptor-like kinase FLS2, brassinosteroid insensitive 1-associated receptor kinase 1 (BAK1), LRR receptor-like serine/threonine-protein kinase EFR, calcium-dependent protein kinase (CDPK), calmodulin/calcium-binding protein (CaM/CML), mitogen-activated protein kinase kinase kinase 1 (MEKK1), and WRKY29 tended to be downregulated. However, DEGs encoding NADPHox (Rboh) and mitogen-activated protein kinase kinase 1 (MKK1/2) tended to be upregulated in inoculated compared with control plants. Regarding effector-triggered immunity (ETI), the endophytic inoculation particularly decreased the expression of genes encoding the disease resistance proteins RPM1, RPS2, and RPS5. Additionally, defense genes such as WRKY and pathogenesis-related protein1 (PR1) also tended to be repressed after AL12 inoculation.

FIGURE 5

Figure 5. The fungal endophyte Gilmaniella sp. AL12 affects the expression of genes involved in plant-pathogen interactions in shoots of A. lancea. Red boxes indicate most upregulated genes. Green boxes indicate most downregulated genes. Gray boxes indicate equally upregulated and downregulated genes. CNGC, cyclic nucleotide gated channel; FLS2, LRR receptor-like serine/threonine-protein kinase FLS2; BAK1, brassinosteroid insensitive 1-associated receptor kinase 1; EFR, LRR receptor-like serine/threonine-protein kinase EFR; CDPK, calcium-dependent protein kinase; Rboh, respiratory burst oxidase; CaM/CML, calmodulin/calcium-binding protein CML MEKK1, mitogen-activated protein kinase kinase kinase 1; MKK1/2, mitogen-activated protein kinase kinase 1; WRKY 25, transcription factor WRKY25; WRKY 29, transcription factor WRKY29; RPM1, disease resistance protein RPM1; RPS2, disease resistance protein RPS2; RPS5, disease resistance protein RPS5; PBS1, serine/threonine-protein kinase; MIN7, guanine nucleotide-exchange factor; HSP90, heat shock protein 90; MYC2, transcription factor MYC2.

Gilmaniella sp. AL12-Regulated Genes Involved in Signaling

Our result shows that DEGs involved in auxin, cytokinin, and ethylene biosynthesis tended to be upregulated (Figure 6), and DEGs involved in cytokinin and ethylene signal transduction also tended to be upregulated (Figure 6). In total, 99 TF-coding genes were differentially expressed between CK and AL12 samples (Supplementary Figure S5). These DEGs were divided into 17 families, including WRKY, APETALA2/Ethylene-Response Factors (AP2/ERF), MYB, basic Helix-Loop-Helix (bHLH), DELLA, Heat stress, Trilelix, GATA, and NAC, among others. Of these families, DEGs encoding MYB, Trilelix, and NAC transcription factor tended to be upregulated, while many genes encoding AP2/ERF, WRKY, DELLA, Heat stress and GATA were downregulated.

FIGURE 6

Figure 6. Distribution of differentially expression genes involved in hormone biosynthesis and signal transduction in shoots of A. lancea after the fungal endophyte Gilmaniella sp. AL12 inoculation. The x axis presents number of DEGs. The y axis indicates different hormone signals.

RT-qPCR Validation of Gene Expression

To further validate the transcriptomic results, 15 transcripts related to primary metabolism, secondary metabolism, and defense were selected for RT-qPCR analysis (Supplementary Table S2 and Figure 7). As shown in Figure 7, the results of the RT-qPCR analysis were largely consistent with the RNA sequencing analysis, supporting the high quality of the RNA sequencing datasets. In particular, the expression of genes annotated as GAPDH, GS, and β-caryophyllene synthase CPS1 markedly increased after endophytic fungus inoculation (Figure 7).

FIGURE 7

Figure 7. RT-qPCR validation of differentially expressed genes involved in primary metabolism, secondary metabolism, and the defense response in shoots of A. lancea after the fungal endophyte Gilmaniella sp. AL12 inoculation. Elongation factor 1 alpha gene (EF1a) was used as an internal reference, and the 2^−ΔΔCt method was used to analyze the relative mRNA expression. Values are the means of three independent experiments. Bars represent standard deviations. RCA, ribulose-1,5-bisphosphate carboxylase/oxygenase activase; GAPDH, glyceraldehyde 3-phosphate dehydrogenase; PDHE, pyruvate dehydrogenase; GS, glutamine synthetase; GOGAT, glutamate synthase; DXS, 1-deoxy-d-xylulose 5-phosphate synthase; FS, β-farnesene synthase; CPS1, β-caryophyllene synthase 1; CPS2, β-caryophyllene synthase 2; PLD, Phospholipase D; FLS2, LRR receptor-like kinase FLS2; RPM1, disease resistance protein RPM1; RPS5, disease resistance protein RPS5; WRKY29, transcription factor WRKY29; PR1, pathogenesis-related protein1.

Shoot Proteome Changes by Gilmaniella sp. AL12 Inoculation Involved Diverse Biological Processes

To identify proteins that responded to the endophyte inoculation, we employed 2-DE to compare the shoot protein profiles of the CK and AL12 samples. A total of approximately 2136 protein spots were separated (Supplementary Figure S6). After optimization of the 2-DE gels, approximately 129 differentially expressed proteins with at least a two-fold change were identified (Supplementary Figure S6). Finally, 125 proteins were successfully identified via MALDI-TOF MS/MS (Supplementary Figure S7 and Supplementary Table S10). These proteins were classified into nine functional categories (Supplementary Figure S7), including “energy metabolism,” “carbohydrate metabolism,” “amino acid metabolism,” “lipid metabolism,” “defense/stress,” “plant secondary metabolism,” “signal transduction,” “transcription and translation,” and “cell growth/division.”

The expression of energy metabolism-related proteins, including ribulose bisphosphate carboxylase, ferredoxin-nitrite reductase, chloroplast oxygen-evolving enhancer protein, and phosphoenolpyruvate carboxylase, were upregulated (Table 1). Regarding carbohydrate metabolism, glyceraldehyde-3-phosphate dehydrogenase and NAD-dependent glyceraldehyde-3-phosphate dehydrogenase were up-regulated, whereas carbon catabolite repressor protein was down-regulated (Table 1). Four plant secondary metabolism-associated proteins were upregulated, including tropinone reductase, Type III polyketide synthase and cytochrome P450 CYP76F12 (Table 1).

TABLE 1

Table 1. Identification of differentially expressed proteins related to energy, carbohydrate, and secondary metabolism in shoots of A. lancea after endophyte inoculation.

Among the 30 proteins involved in signal transduction, 17 were upregulated. These proteins included membrane-associated protein, protein kinase, protein phosphatase, abscisic acid (ABA) responsive protein, gibberellin (GA) biosynthesis related proteins, indoleacetic acid (IAA) biosynthesis related proteins, and a few transcription factors (Supplementary Table S10). The other 13 downregulated proteins incorporated pentapeptide, ethylene biosynthesis-related protein, and ubiquitin conjugate factors, among others (Supplementary Table S10). Additionally, AL12-regulated proteins involved in defense/stress response, amino acid and lipid metabolism, transcription and translation, cell growth, and other function were differentially regulated (Supplementary Table S10).

Discussion

The pharmaceutical value of medicinal plants relies on the accumulation of active pharmaceutical ingredient, and guaranteeing the yield and quality of these herbs is the main challenge (Yang et al., 2012). Endophytes have been proven to exert multiple effects on their host plants, including plant growth promotion, secondary metabolite biosynthesis promotion, stress resistance enhancement (Ludwig-Müller, 2015; Yang et al., 2016; Dastogeer et al., 2018). The advantages provided by the endophyte fungus AL12 on the host plant A. lancea were an improved plant biomass and increased sesquiterpenoids content (Yuan et al., 2016b). In this study, transcriptomics and proteomics were employed to analyze how endophytic fungi affect the regulation of transcription and translation in A. lancea. Compared to the reported endophyte associations (Doehlemann et al., 2008; Kawahara et al., 2012; Morán-Diez et al., 2012; Perazzolli et al., 2012; De Cremer et al., 2013; Zouari et al., 2014; Adolfsson et al., 2017; Chen W. et al., 2017; Dinkins et al., 2017; Hao et al., 2017; Bajaj et al., 2018), the endophyte Gilmaniella sp. AL12 resulted in 2.7% differential gene expression, thus demonstrating a greater impact on their host than most other endophytes (Supplementary Table S11). During the plant life cycle, a dynamic trade-off between growth and defense is necessary for plant resource assignment in response to multiple developmental cues and environmental stimuli (Hou et al., 2013). In contrast to the common downregulated plant-pathogen interactions and most phytohormone signaling events, Gilmaniella sp. AL12 upregulated genes involved in primary and secondary metabolism (Figures 2–4, 8). Consistent with the improved cytokinin content of A. lancea after AL12 inoculation (Yuan et al., 2016b), DEGs involved in cytokinin biosynthesis and signal transduction both tended to be upregulated (Figure 6), which may induce cell division, accelerate chlorophyll biosynthesis, and delay leaf senescence (Cortleven and Schmülling, 2015). In addition, ethylene induced by the endophyte mediates sesquiterpenoid biosynthesis in A. lancea (Yuan et al., 2016a). Based on these results, we propose that colonization by the endophytic fungus Gilmaniella sp. AL12 will shift plant metabolism from defense to growth, allowing the host plant to utilize limited energy for carbon assimilation and, thus, to achieve an increased biomass and sesquiterpenoid content in A. lancea.

It is acknowledged that plant primary metabolism and secondary metabolism required large amounts of energy. One way to support the increased energy demands is to enhance the carbon assimilation efficiency (Rozpądek et al., 2015, 2019). We have previously observed that the net photosynthesis rate and contents of chlorophyll, rubisco, and soluble carbohydrate in A. lancea increase after AL12 inoculation (Yuan et al., 2016b). In the present study, phenomenological fluxes per cross section of the host plant was improved after AL12 inoculation (Figure 1D), indicating the improvement of PS II reaction center performance in A. lancea. As shown by the transcriptome and proteome results, endophyte inoculation would enhance NADPH production and CO₂ assimilation (Table 1, Figures 3, 8 and Supplementary Figure S8), thus explaining the positive effects of endophyte inoculation on plant photosynthesis. Similarly, Dactylis glomerate inoculated with Epichloë typhina, Medicago truncatula, and Populus alba inoculated with AMF all displayed a similar photosynthesis performance with improved plant light-driven energy production efficiency, Rubisco content, CO₂ assimilation and PS II photochemistry efficiency (Lingua et al., 2008; Aloui et al., 2009; Cicatelli et al., 2010, 2012; Rozpądek et al., 2015, 2019). Although Gilmaniella sp. AL12, Epichloë typhina and AMF belong to distinct fungal species, they all provide a benefit, such as improving similar photosynthesis apparatuses of their host plant. Plant photosynthesis and respiration processes exist with each other interdependently and share adenosine diphosphate (ADP) and NADP⁺ or NAD⁺. Additionally, O₂ produced by photosynthesis can be utilized by the respiration process, and CO₂ produced by respiration can also be assimilated by photosynthesis. As speculated based on the results of this study (Table 1 and Figure 3), NADPH and O₂ production by chloroplasts, as well as NADH dehydrogenation and triphosadenine (ATP) production by mitochondria, were improved after AL12 inoculation, supporting energy storage for plant metabolism. The TCA cycle mainly converts pyruvate to malate to produce ATP, incorporating citric acid biosynthesis, oxidation and decarboxylation, and OAA regeneration (Fernie et al., 2004). Epichloë typhina has been shown to markably strengthen NADPH-MDH enzyme activity in its host Dactylis glomerate (Rozpądek et al., 2015). Likewise, the endophytic fungus Gilmaniella sp. AL12 upregulated genes encoding MDH and pyruvate dehydrogenase (PDHE) (Figure 3), thus contributing to acetyl-CoA biosynthesis and OAA regeneration (Figure 8). Acetyl-CoA and pyruvate biosynthesis were improved in A. lancea after endophyte inoculation (Figure 8), indicating that more precursor substances and energy were available for plant secondary metabolite biosynthesis (Nema et al., 2013). As the precursor substances of terpenoid biosynthesis (Chen F. et al., 2017), pyruvate and acetyl-CoA content increased after AL12 inoculation (Yuan et al., 2016b). Additionally, in contrast to the general downregulation of fatty acid biosynthesis (Supplementary Table S4), genes associated with terpenoid biosynthesis were upregulated after AL12 inoculation (Supplementary Figure S3), indicating that a greater amount of acetyl-CoA would be diverted into terpene biosynthesis.

FIGURE 8

Figure 8. Simplified schematic summary of metabolic changes in shoots of A. lancea after the fungal endophyte Gilmaniella sp. AL12 inoculation. DEGs with significant differences in endophyte inoculated plants are shown as colored arrows (upregulated DEGs in red and downregulated DEGs in blue). The yellow ellipse indicates the phosphate translocator (TPT). Letters in green blocks represent precursor substances of terpenoid biosynthesis. C3 cycle, Benson-Calvin cycle; GS/GOGAT, glutamine synthetase (GS)/glutamate synthase (GOGAT); TCA cycle, tricarboxylic acid cycle; PS I, photosystem I; PS II, photosystem II; 3-PGA, glycerate-3-phosphate; OAA, oxaloacetate; 2-OG, 2-oxoglutarate; FNR, ferredoxin-NADP reductase; OEE, photosystem II oxygen-evolving enhancer protein.

Plant terpenoids are derived from isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP) synthesized through the mevalonic acid (MVA) and the 2-C-methyl-D-erythritol-4-phosphate (MEP) pathway (Vattekkatte et al., 2018). The different molecular rates of isopentenyl diphosphate and dimethylallyl diphosphate condensation lead to geranyl diphosphate (GPP) for monoterpene, farnesyl diphosphate (FPP) for sesquiterpene and geranylgeranyl diphosphate (GGPP) for triterpene biosynthesis (Dudareva et al., 2013; Vattekkatte et al., 2018). Consistent with the increased sesquiterpenoids content in A. lancea shoots (Figure 1A), DEGs associated with terpenoid backbone biosynthesis and sesquiterpene biosynthesis were increased (Supplementary Figure S3). Unfortunately, we did not detect any spot identified associated with terpene backbone biosynthesis or sesquiterpenoids biosynthesis in proteome (Table 1). Considering that there is no database of A. lancea or related genus in NCBI, we identified proteins using NCBI green plant database. Thus, we considered that maybe some proteins related to secondary metabolism may not identified exactly. Yet despite that, comparative transcriptomics showed that genes related to terpene backbone biosynthesis and sesquiterpenoids biosynthesis tended to be upregulated in shoots of A. lancea after the endophyte inoculation (Supplementary Figure S3). In this study, high-quality reads were assembled de novo by using the Trinity program. The results showed that genes annotated as β-farnesene synthase and β-caryophyllene synthase were upregulated after AL12 inoculation, as shown by the transcriptome (Supplementary Figure S4) and RT-qPCR results (Figure 7). It is known that one terpene synthase can catalyze formation more than one terpene (Dudareva et al., 2013). For example, TPS21 and TPS11 from Arabidopsis thaliana are responsible for the biosynthesis of more than 20 sesquiterpenoids in flower fragrance (Tholl et al., 2005). Additionally, OkBCS of Ocimum can catalyze the production of β-caryophyllene and α-humulene (Jayaramaiah et al., 2016). The main sesquiterpenoids of A. lancea incorporate β-caryophyllene, zingiberene, β-sesquiphellandrene, caryophyllene oxide, hinesol, β-eudesmol and atractylone (Yuan et al., 2016b). Considering the multifunction of sesquiterpene synthase, we propose that β-farnesene synthase and β-caryophyllene synthase are crucial terpene synthases induced by AL12 to supply higher sesquiterpenoid contents in A. lancea. Further experiments will include the functional characterization of these sesquiterpene synthases, and how the fungal endophyte regulate genes coding for these sesquiterpene synthases.

The accumulation of secondary metabolites is regulated by cross-talking signaling cascades. Our previous studies have demonstrated that the endophyte AL12 activates multiple signals to induce plant sesquiterpenoids accumulation (Wang et al., 2011; Ren and Dai, 2012, 2013; Yuan et al., 2016a). In addition, mannan-binding lectin has been proposed to recognize mannan of Gilmaniella sp. AL12, thus causing signal transduction (Chen et al., 2016). It is worth noting that 32 DEGs encoding PLD tended to be upregulated after AL12 inoculation in the present study (Supplementary Table S8). PLD can hydrolyse phospholipids into phospholipid acid, which acts as an important molecular signal mediating heat stress-induced triterpene ganoderic acid biosynthesis in Ganoderma lucidum (Liu et al., 2017). These results indicate that phospholipid acid may be an essential signaling molecule in A. lancea after AL12 inoculation. Additionally, metabolomics and transcriptomics are required to study whether PLD or phospholipid signaling in A. lancea are activated and whether they mediate sesquiterpenoid biosynthesis after AL12 inoculation, which will help to understand their roles in inducing the biosynthesis of the secondary metabolites.

In contrast to the enhanced primary and secondary metabolism most genes involved in the plant-pathogen interaction were downregulated (Figures 2, 5). Plants evolve pattern recognition receptors (PRRs) to recognize evolutionarily conservative PAMPs, triggering the host MAPK cascade (Figure 5). Then, MAPK activates plant hormone signal transduction to integrate diverse aspects of PTI (Cui et al., 2015; Pandey et al., 2016). These optional signaling events will give rise to various immune responses (Akira and Hemmi, 2003). Furthermore, effectors produced by the fungus are also known to inhibit MAPK cascade signaling or block plant-resistant protein expression, thus suppressing plant basic immune response PTI (Zhang et al., 2012; Meng and Zhang, 2013). Fungi deploy these co-evolved effectors to modify host plant cell processes or to establish biotrophic or other types of symbiotic relationship with the host (Le Fevre et al., 2015). In turn, fungal effectors can be recognized by plant disease resistance proteins (R proteins) and thus trigger plant ETI, which usually causes the plant hypersensitive response (HR) (Jones and Dangl, 2006; Liu et al., 2007). With endophytic fungal AL12 inoculation, genes annotated as plant PRRs, including LRR receptor-like kinase FLS2, BAK1 and EFR, tended to be downregulated (Figure 5), indicating a weakened PAMP recognition ability of the host plant. Furthermore, genes involved in the response of PAMPs or effectors were also repressed (Figure 5), suggesting that endophyte inoculation suppressed host immune reactions. As shown in our previous study (Yang and Dai, 2013), AL12 could successfully colonize the leaves of A. lancea, resulting in an altered hyphal morphology such as smooth and unbranched hyphae. In addition, slightly pointed hyphal tips of the endophyte also appears similar as the ends of a nematode (Yang and Dai, 2013). Considering that genes related to PTI and ETI of A. lancea were downregulated after AL12 inoculation (Figure 5) and that the hyphae of the endophyte were altered morphologically (Yang and Dai, 2013), we deduced that either masked fungal PAMPs were not detected by the host plant, or that the plant defense response was suppressed by the host or the endophyte (Dupont et al., 2015), potentially contributing to the compatible association of Gilmaniella sp. AL12 with A. lancea.

Given the apparent downregulation of genes involved in the plant-pathogen interaction, endophyte inoculation would enhance the pathogen susceptibility of the host plant. However, this is not the actual fact. Although plant defense-related genes of Perennial ryegrass were down-regulated after Epichloë festucae var. lolii infection, its resistance to fungal pathogens was enhanced due to the production of phenolic compounds (Tian et al., 2008; Pańka et al., 2013; Dupont et al., 2015). Similarly, Gilmaniella sp. AL12 upregulated genes involved in the biosynthesis of phenylpropanoids and terpenes (Figure 4). Terpenoids are frequently used as phytoalexins, the accumulation of which are enhanced by biotic or abiotic stress (Schmelz et al., 2014; Vaughan et al., 2015). The sesquiterpenoid phytoalexins such as gossypol, hemigossypolone and heliocides can provide defense mechanisms against pathogens and herbivores of cotton (Yang C. Q. et al., 2013). Additionally, β-caryophyllene can defend against Helicoverpa armigera and Pseudomonas syringae pv. tomato DC3000 (Huang et al., 2012; Singh et al., 2014). Similarly, the monoterpene S-limonene and triterpene momilactones or oryzalexins function as defensive metabolites against Magnaporthe oryzae or Magnaporthe grisea (Toyomasu et al., 2008; Chen et al., 2018). Therefore, although defense-related genes of A. lancea were downregulated after Gilmaniella sp. AL12 inoculation, the enhanced synthesis of phenylpropanoids and terpenoids may function as defensive strategy against other microbial pathogens or herbivores, functioning as a possible competitive exclusion tactic of the endophyte.

In this study, shoots of A. lancea were sterile, and were appropriate for investigating the effect of the fungal endophyte on plantlets without the distraction of other biotic or abiotic factors, thus helping to understand plant-endophyte interactions. The above results showed that the fungal endophyte Gilmaniella sp. AL12 weakened the plant immune response, which might contribute to its successful and stable colonization of the host plant (Figure 9). The decreased plant immunity resulted in large amounts of energy for plant primary or secondary metabolism. Moreover, the fungal endophyte improved photosynthesis, carbon fixation, carbohydrate transformation, and the conversion from pyruvate to acetyl-CoA in the host plant, indicating that the endophyte promoted the accumulation of precursors for terpenoid biosynthesis. In addition, the fungal endophyte specifically induced sesquiterpenoid biosynthesis by regulating genes involved in sesquiterpenoid biosynthesis and related signaling events in A. lancea. Thus, we propose that the fungal endophyte-plant association increased production (biomass and sesquiterpenoid content) by increasing the source (photosynthesis), expanding the sink (glycolysis and the TCA cycle), and enhancing metabolic flux (sesquiterpenoids biosynthesis-related genes) in A. lancea (Figure 9). Consistent with the “double promotion” in this study, the endophytic fungi AL12 promotes plant growth and sesquiterpenoid accumulation within two years of growth in field experiments. Further studies will be conducted to investigate whether pre-inoculation with the fungal endophyte Gilmaniella sp. AL12 regulate metabolism and immunity of A. lancea in the field experiments, which will help to understand plant-endophyte interactions, and will contribute to the application of the fungal endophyte in cultivating of A. lancea.

FIGURE 9

Figure 9. Simplified schematic summary of the “double-promotion” effect on A. lancea bring by the fungal endophyte Gilmaniella sp. AL12. The fungal endophyte Gilmaniella sp. AL12 weakened the plant immune response of A. lancea, thus may contributing to plant-endophyte beneficial interaction. Additionally, endophyte-plant association increased production (biomass and sesquiterpenoid content) by increasing the source (photosynthesis), expanding the sink (glycolysis and the TCA cycle), and enhancing metabolic flux (sesquiterpenoids biosynthesis-related genes) in A. lancea.

Medicinal plants are rich in active compounds such as artemisinin (Sharma and Agrawal, 2013) and ginseng saponin (Wei et al., 2018), which are important sources of modern drugs. You-You Tu, the mother of artemisinin, was awarded Nobel Prize in medicine in 2015. The development of the medicinal plants gradually becomes a focused issue, and medicinal plants-endophytes interactions have received much attentions. A review of the medicinal plant microbiome has shown that the secondary metabolite content of medicinal herbs depends on where they are cultivated, which can be partly ascribed to different rhizospheric or endophytic microbes associated with their cultivation location (Köberl et al., 2013; Huang et al., 2018). Plantlets of A. lancea grown in the Maoshan area of southeast China represent the geo–authentic medicinal plant, with a much greater sesquiterpenoid content and diversity than A. lancea in other cultivation area (Zhou et al., 2018). The endophytic fungus Gilmaniella sp. AL12 was isolated from the stem of the geo–authentic A. lancea, and the compatible association contributes some benefits to the host plant, including the promotion of plant growth and sesquiterpenoid content, and the prevention of root rot (Wang et al., 2011; Ren and Dai, 2012, 2013; Wang et al., 2012; Yuan et al., 2016a,b; Ren et al., 2017). Whether Gilmaniella sp. AL12 is a key member of the core microbiome of geo–authentic A. lancea remain unanswered. In this regard, our previous studies have shown that prioritizing endophytic Gilmaniella sp. AL12 inoculation enhances the diversity and size of phyllospheric microbial communities of A. lancea by increasing the soluble sugar content in the rhizosphere (Yang T. et al., 2013). Given that the endophyte provided physiological and biochemical benefits to the host and affected host phyllospheric microbial communities, we more broadly propose a beneficial association of endophytic Gilmaniella sp. AL12 with A. lancea as a potential model for endophytic fungi-medicinal herb interaction. The endophyte-A. lancea association contributes to medicinal herb cultivation and helps to further clarify plant-endophyte interactions.

Conclusion

In summary, this study showed that the fungal endophyte Gilmaniella sp. AL12 weakened the plant immune response of A. lancea, thus may contributing to the beneficial plant-endophyte interaction. Additionally, endophyte-plant association increased production (biomass and sesquiterpenoid content) by increasing the source (photosynthesis), expanding the sink (glycolysis and the TCA cycle), and enhancing metabolic flux (sesquiterpenoids biosynthesis-related genes) in A. lancea. This study revealed the regulation of Gilmaniella sp. AL12 on plant metabolism and related regulatory processes in shoots of A. lancea at the transcriptional and translational level. This study provides a theoretical basis for medicinal herb cultivation and helps to further clarify plant-endophyte interactions.

Data Availability

The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.

Author Contributions

C-CD and JY designed the experiments, analyzed the data, and wrote the manuscript. P-XC and M-JT helped extract plant RNA and protein. WZ, KS, and XL helped supervise the manuscript writing. All the authors read and approved the final manuscript.

Funding

This work was supported by the National Key Research and Development Program of China (Grant No. 2017YFD0800705).

Conflict of Interest Statement

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/fmicb.2019.01208/full#supplementary-material

Abbreviations

A. lancea, Atractylodes lancea; AL12, Gilmaniella sp. AL12; AMF, arbuscular mycorrhizal fungi; ATP, triphosadenine; BAK1, brassinosteroid insensitive 1-associated receptor kinase 1; DEGs, differentially expressed genes; DPI, days post-AL12 inoculation; EFR, LRR receptor-like serine/threonine-protein kinase EFR; ETI, effector triggered immunity; FPKM, the fragments per kilobase of transcript per million mapped reads; GC, gas chromatography; GO, gene ontology database; KEGG, Kyoto Encyclopedia of Genes and Genomes database; LRR, leucine-rich repeat; MDH, malate dehydrogenase; NCBI, National Center for Biotechnology Information; OAA, oxaloacetate; PAMPs, pathogen-associated molecular patterns; Pi, orthophosphate; PLD, phospholipase D; PS II, photosystem II; PTI, PAMP triggered immunity; RT-qPCR, real-time quantitative PCR; TCA cycle, tricarboxylic acid cycle; 2-DE, two-dimensional gel electrophoresis.

Footnotes

1. ^ http://pathways.embl.de/iPath2.cgi#

References

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

Akira, S., and Hemmi, H. (2003). Recognition of pathogen-associated molecular patterns by TLR family. Immunol. Lett. 85, 85–95. doi: 10.1016/S0165-2478(02)00228-6

CrossRef Full Text | Google Scholar

Aloui, A., Recorbet, G., Gollotte, A., Robert, F., Valot, B., Gianinazzi-Pearson, V., et al. (2009). On the mechanisms of cadmium stress alleviation in Medicago truncatula by arbuscular mycorrhizal symbiosis: a root proteomic study. Proteomics 9, 420–433. doi: 10.1002/pmic.200800336

PubMed Abstract | CrossRef Full Text | Google Scholar

Bajaj, R., Huang, Y., Gebrechristos, S., Mikolajczyk, B., Brown, H., Prasad, R., et al. (2018). Transcriptional responses of soybean roots to colonization with the root endophytic fungus Piriformospora indica reveals altered phenylpropanoid and secondary metabolism. Sci. Rep. 8:10227. doi: 10.1038/s41598-018-26809-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, F., Ren, C. G., Zhou, T., Wei, Y. J., and Dai, C. C. (2016). A novel exopolysaccharide elicitor from endophytic fungus Gilmaniella sp. AL12 on volatile oils accumulation in Atractylodes lancea. Sci. Rep. 6:34735. doi: 10.1038/srep34735

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, F., Wei, Y. X., Zhang, J. M., Sang, X. M., and Dai, C. C. (2017). Transcriptomics analysis investigates sesquiterpenoids accumulation pattern in different tissues of Atractylodes lancea (Thunb.) DC. plantlet. Plant Cell Tiss. Org. 130, 73–90. doi: 10.1007/s11240-017-1205-8

CrossRef Full Text | Google Scholar

Chen, W., Li, J., Zhu, H., Xu, P., Chen, J., and Yao, Q. (2017). Arbuscular mycorrhizal fungus enhances lateral root formation in Poncirus trifoliata (L.) as revealed by RNA-seq analysis. Front. Plant Sci. 8:2039. doi: 10.3389/fpls.2017.02039

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Chen, H., Yuan, J. S., Köllner, T. G., Chen, Y., Guo, Y., et al. (2018). The rice terpene synthase gene OsTPS19 functions as an (S)-limonene synthase in planta, and its overexpression leads to enhanced resistance to the blast fungus Magnaporthe oryzae. Plant Biotechnol. J. 16, 1778–1787. doi: 10.1111/pbi.12914

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J. X., Dai, C. C., Li, X., Tian, L. S., and Xie, H. (2008). Endophytic fungi screening from Atractylodes lancea and inoculating into the host plantlet. Guihaia 28, 256–260.

Google Scholar

Cicatelli, A., Lingua, G., Todeschini, V., Biondi, S., Torrigiani, P., and Castiglione, S. (2010). Arbuscular mycorrhizal fungi restore normal growth in a white poplar clone grown on heavy metal-contaminated soil, and this is associated with upregulation of foliar metallothionein and polyamine biosynthetic gene expression. Ann. Bot. 106, 791–802. doi: 10.1093/aob/mcq170

PubMed Abstract | CrossRef Full Text | Google Scholar

Cicatelli, A., Lingua, G., Todeschini, V., Biondi, S., Torrigiani, P., and Castiglione, S. (2012). Arbuscular mycorrhizal fungi modulate the leaf transcriptome of a Populus alba L. clone grown on a zinc and copper-contaminated soil. Environ. Exp. Bot. 75, 25–35. doi: 10.1016/j.envexpbot.2011.08.012

CrossRef Full Text | Google Scholar

Cortleven, A., and Schmülling, T. (2015). Regulation of chloroplast development and function by cytokinin. J. Exp. Bot. 66, 4999–5013. doi: 10.1093/jxb/erv132

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, H. T., Tsuda, K., and Parker, J. E. (2015). Effector-triggered immunity: from pathogen perception to robust defense. Annu. Rev. Plant Biol. 66, 487–511. doi: 10.1146/annurev-arplant-050213-040012

PubMed Abstract | CrossRef Full Text | Google Scholar

Dastogeer, K. M. G., Li, H., Sivasithamparam, K., Jones, M. G. K., and Wylie, S. J. (2018). Fungal endophytes and a virus confer drought tolerance to Nicotiana benthamiana plants through modulating osmolytes, antioxidant enzymes and expression of host drought responsive genes. Env. Exp. Bot. 149, 95–108. doi: 10.1016/j.envexpbot.2018.02.009

CrossRef Full Text | Google Scholar

De Cremer, K., Mathys, J., Vos, C., Froenicke, L., Michelmore, R. W., Cammue, B. P., et al. (2013). RNAseq-based transcriptome analysis of Lactuca sativa infected by the fungal necrotroph Botrytis cinereal. Plant Cell Environ. 36, 1992–2007. doi: 10.1111/pce.12106

PubMed Abstract | CrossRef Full Text | Google Scholar

Dinkins, R. D., Nagabhyru, P., Graham, M. A., Boykin, D., and Schardl, C. L. (2017). Transcriptome response of Lolium arundinaceum to its fungal endophyte Epichloë coenophiala. New Phytol. 213, 324–337. doi: 10.1111/nph.14103

PubMed Abstract | CrossRef Full Text | Google Scholar

Doehlemann, G., Wahl, R., Horst, R. J., Voll, L. M., Usadel, B., Poree, F., et al. (2008). Reprogramming a maize plant: transcriptional and metabolic changes induced by the fungal biotroph Ustilago maydis. Plant J. 2008, 181–195. doi: 10.1111/j.1365-313X.2008.03590.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dudareva, N., Klempien, A., Muhlemann, J. K., and Kaplan, I. (2013). Biosynthesis, function and metabolic engineering of plant volatile organic compounds. New Phytol. 198, 16–32. doi: 10.1111/nph.12145

PubMed Abstract | CrossRef Full Text | Google Scholar

Dudareva, N., Negre, F., Nagegowda, D. A., and Orlova, I. (2006). Plant volatiles: recent advances and future perspectives. Crit. Rev. Plant Sci. 25, 417–440. doi: 10.1080/07352680600899973

CrossRef Full Text | Google Scholar

Dupont, P. Y., Eaton, C. J., Wargent, J. J., Fechtner, S., Solomon, P., Schmid, J., et al. (2015). Fungal endophyte infection of ryegrass reprograms host metabolism and alters development. New Phytol. 208, 1227–1240. doi: 10.1111/nph.13614

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernie, A. R., Carrari, F., and Sweetlove, L. J. (2004). Respiratory metabolism: glycolysis, the TCA cycle and mitochondrial electron transport. Curr. Opin. Plant Biol. 7, 254–261. doi: 10.1016/j.pbi.2004.03.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Hao, K., Wang, F., Nong, X., McNeill, M. R., Liu, S., Wang, G., et al. (2017). Response of peanut Arachis hypogaea roots to the presence of beneficial and pathogenic fungi by transcriptome analysis. Sci. Rep. 7:964. doi: 10.1038/s41598-017-01029-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Hou, X. L., Ding, L. H., and Yu, H. (2013). Crosstalk between GA and JA signaling mediates plant growth and defense. Plant Cell Rep. 32, 1067–1074. doi: 10.1007/s00299-013-1423-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, M., Sanchez-Moreiras, A. M., Abel, C., Sohrabi, R., Lee, S., Gershenzon, J., et al. (2012). The major volatile organic compound emitted from Arabidopsis thaliana flowers, the sesquiterpene (E)-beta-caryophyllene, is a defense against a bacterial pathogen. New Phytol. 193, 997–1008. doi: 10.1111/j.1469-8137.2011.04001.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, W., Long, C., and Lam, E. (2018). Roles of plant-associated microbiota in traditional herbal medicine. Trends Plant Sci. 23, 559–562. doi: 10.1016/j.tplants.2018.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Jayaramaiah, R. H., Anand, A., Beedkar, S. D., Dholakia, B. B., Punekar, S. A., Kalunke, R. M., et al. (2016). Functional characterization and transient expression manipulation of a new sesquiterpene synthase involved in beta-caryophyllene accumulation in Ocimum. Biochem Bioph Res. Co. 473, 265–271. doi: 10.1016/j.bbrc.2016.03.090

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, J. D., and Dangl, J. L. (2006). The plant immune system. Nature 444, 323–329. doi: 10.1038/nature05286

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawahara, Y., Oono, Y., Kanamori, H., Matsumoto, T., Itoh, T., and Minami, E. (2012). Simultaneous RNA-seq analysis of a mixed transcriptome of rice and blast fungus interaction. PLoS One 7:e49423. doi: 10.1371/journal.pone.0049423

PubMed Abstract | CrossRef Full Text | Google Scholar

Köberl, M., Schmidt, R., Ramadan, E. M., Bauer, R., and Berg, G. (2013). The microbiome of medicinal plants: diversity and importance for plant growth, quality and health. Front. Microbiol. 4:400. doi: 10.3389/fmicb.2013.00400

PubMed Abstract | CrossRef Full Text | Google Scholar

Koonrungsesomboon, N., Na-Bangchang, K., and Karbwang, J. (2014). Therapeutic potential and pharmacological activities of Atractylodes lancea (Thunb.) DC. Asian Pac. J. Trop. Med. 7, 421–428. doi: 10.1016/S1995-7645(14)60069-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Le Fevre, R., Evangelisti, E., Rey, T., and Schornack, S. (2015). Modulation of host cell biology by plant pathogenic microbes. Annu. Rev. Cell Dev. Biol. 31, 201–229. doi: 10.1146/annurev-cellbio-102314-112502

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, P. M., Gao, H. Y., and Reto, J. S. (2005). Application of the fast chlorophyll fluorescence induction dynamics analysis in photosynthesis study. J. Plant Physiol. Mol. Biol. 31, 559–566. doi: 10.3321/j.issn:1671-3877.2005.06.001

CrossRef Full Text | Google Scholar

Lingua, G., Franchin, C., Todeschini, V., Castiglione, S., Biondi, S., Burlando, B., et al. (2008). Arbuscular mycorrhizal fungi differentially affect the response to high zinc concentrations of two registered poplar clones. Environ. Pollut. 153, 137–147. doi: 10.1016/j.envpol.2007.07.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, G., Liu, J., Zhang, C., You, X., Zhao, T., Jiang, J., et al. (2018). Physiological and RNA-seq analyses provide insights into the response mechanism of the Cf-10-mediated resistance to Cladosporium fulvum infection in tomato. Plant Mol. Biol. 96, 403–416. doi: 10.1007/s11103-018-0706-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Ren, D., Pike, S., Pallardy, S., Gassmann, W., and Zhang, S. (2007). Chloroplast-generated reactive oxygen species are involved in hypersensitive response-like cell death mediated by a mitogen-activated protein kinase cascade. Plant J. 51, 941–954. doi: 10.1111/j.1365-313X.2007.03191.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y. N., Lu, X. X., Chen, D., Lu, Y. P., Ren, A., Shi, L., et al. (2017). Phospholipase D and phosphatidic acid mediate heat stress induced secondary metabolism in Ganoderma lucidum. Environ. Microbiol. 19, 4657–4669. doi: 10.1111/1462-2920.13928

PubMed Abstract | CrossRef Full Text | Google Scholar

Ludwig-Müller, J. (2015). Plants and endophytes: equal partners in secondary metabolite production? Biotech Lett. 37, 1325–1334. doi: 10.1007/s10529-015-1814-s1814

PubMed Abstract | CrossRef Full Text | Google Scholar

Lv, Z. Y., Zhang, L., and Tang, K. X. (2017). New insights into artemisinin regulation. Plant Signal Behav. 12:10. doi: 10.1080/15592324.2017.1366398

PubMed Abstract | CrossRef Full Text | Google Scholar

Mandal, S., Evelin, H., Girl, B., Singh, V. P., and Kapoor, R. (2013). Arbuscular mycorrhiza enhances the production of stevioside and rebaudioside-A in Stevia rebaudiana via nutritional and non-nutritional mechanisms. Appl. Soil Ecol. 72, 187–194. doi: 10.1016/j.apsoil.2013.07.003

CrossRef Full Text | Google Scholar

Meng, X. Z., and Zhang, S. Q. (2013). MAPK Cascades in plant disease resistance signaling. Annu. Rev. Phytopathol. 51, 245–266. doi: 10.1146/annurev-phyto-082712-102314

PubMed Abstract | CrossRef Full Text | Google Scholar

Ming, Q., Su, C., Zheng, C., Jia, M., Zhang, Q., Zhang, H., et al. (2013). Elicitors from the endophytic fungus Trichoderma atroviride promote Salvia miltiorrhiza hairy root growth and tanshinone biosynthesis. J. Exp. Bot. 64, 5687–5694. doi: 10.1093/jxb/ert342

PubMed Abstract | CrossRef Full Text | Google Scholar

Morán-Diez, E., Rubio, B., Domínguez, S., Hermosa, R., Monte, E., and Nicolás, C. (2012). Transcriptomic response of Arabidopsis thaliana after 24 h incubation with the biocontrol fungus Trichoderma harzianum. J. Plant Physiol. 169, 614–620. doi: 10.1016/j.jplph.2011.12.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Na-Bangchang, K., Plengsuriyakarn, T., and Karbwang, J. (2017). Research and Development of Atractylodes lancea (Thunb) DC. as a Promising Candidate for Cholangiocarcinoma chemotherapeutics. Evid-based Compl. Alt. 2017:5929234. doi: 10.1155/2017/5929234

PubMed Abstract | CrossRef Full Text | Google Scholar

Nema, R., Khare, S., Jain, P., Pradhan, A., Gupta, A., and Singh, D. (2013). Natural Products Potential and Scope for Modern Cancer Research. Am. J. Plant Sci. 4, 1270–1277. doi: 10.4236/ajps.2013.46157

CrossRef Full Text | Google Scholar

Oh, H. J., Park, J. E., Park, Y. G., and Jeong, B. R. (2014). Growth and quality of plug seedlings of three indigenous medicinal plants as affected by ionic strength of the nutrient solution. Hort. Environ. Biotechnol. 55, 63–69. doi: 10.1007/s13580-014-0074-x

CrossRef Full Text | Google Scholar

Pandey, D., Rajendran, S. R. C. K., Gaur, M., Sajeesh, P. K., and Kumar, A. (2016). Plant defense signaling and responses against necrotrophic fungal pathogens. J. Plant Growth Reg. 35, 1159–1174. doi: 10.1007/s00344-016-9600-7

CrossRef Full Text | Google Scholar

Pańka, D., Piesik, D., Jeske, M., and Baturo-Cieśniewska, A. (2013). Production of phenolics and the emission of volatile organic compounds by perennial ryegrass (Lolium perenne L.)/Neotyphodium lolii association as a response to infection by Fusarium poae. J. Plant Physiol. 170, 1010–1019. doi: 10.1016/j.jplph.2013.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Perazzolli, M., Moretto, M., Fontana, P., Ferrarini, A., Velasco, R., Moser, C., et al. (2012). Downy mildew resistance induced by Trichoderma harzianum T39 in susceptible grapevines partially mimics transcriptional changes of resistant genotypes. BMC Genomics 13:660. doi: 10.1186/1471-2164-13-660

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, C. G., Chen, F., and Dai, C. C. (2017). Fungal endophyte protects Atractylodes lancea from root rot caused by Fusarium oxysporum. Plant Pathol. 66, 223–229. doi: 10.1111/ppa.12567

CrossRef Full Text | Google Scholar

Ren, C. G., and Dai, C. C. (2012). Jasmonic acid is involved in the signaling pathway for fungal endophyte-induced volatile oil accumulation of Atractylodes lancea plantlets. BMC Plant Biol. 12:128. doi: 10.1186/1471-2229-12-128

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, C. G., and Dai, C. C. (2013). Nitric oxide and brassinosteroids mediated fungal endophyte-induced volatile oil production through protein phosphorylation pathways in Atractylodes lancea plantlets. J. Integr. Plant Biol. 55, 1136–1146. doi: 10.1111/jipb.12087

PubMed Abstract | CrossRef Full Text | Google Scholar

Rozpądek, P., Nosek, M., Domka, A., Ważny, R., Jȩdrzejczyk, R. J., Tokarz, K., et al. (2019). Acclimation of the photosynthetic apparatus and alterations in sugar metabolism in response to inoculation with endophytic fungi. Plant Cell Environ. 42, 1408–1423. doi: 10.1111/pce.13485

PubMed Abstract | CrossRef Full Text | Google Scholar

Rozpądek, P., Wezowicz, K., Nosek, M., Ważny, R., Tokarz, K., Lembicz, M., et al. (2015). The fungal endophyte Epichloë typhina improves photosynthesis efficiency of its host orchard grass (Dactylis glomerata). Planta 242, 1025–1035. doi: 10.1007/s00425-015-2337-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmelz, E. A., Huffaker, A., Sims, J. W., Christensen, S. A., Lu, X., Okada, K., et al. (2014). Biosynthesis, elicitation and roles of monocot terpenoid phytoalexins. Plant J. 79, 659–678. doi: 10.1111/tpj.12436

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, E., Anand, G., and Kapoor, R. (2017). Terpenoids in plant and arbuscular mycorrhiza-reinforced defence against herbivorous insects. Ann. Bot. 119, 791–801. doi: 10.1093/aob/mcw263

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, G., and Agrawal, V. (2013). Marked enhancement in the artemisinin content and biomass productivity in Artemisia annua L. shoots co-cultivated with Piriformospora indica. World J. Microb. Biot. 29, 1133–1138. doi: 10.1007/s11274-013-1263-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, S., Rathi, N., Kamal, B., Pundir, D., Kaur, B., and Arya, S. (2010). Conservation of biodiversity of highly important medicinal plants of India through tissue culture technology- a review. Agric. Biol. J. N. Am. 1, 827–833. doi: 10.5251/abjna.2010.1.5.827.833

CrossRef Full Text | Google Scholar

Singh, P., Jayaramaiah, R. H., Sarate, P., Thulasiram, H. V., Kulkarni, M. J., and Giri, A. P. (2014). Insecticidal potential of defense metabolites from Ocimum kilimandscharicum against Helicoverpa armigera. PLoS One 9:e104377. doi: 10.1016/j.compbiolchem.2016.06.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Tatli, O., Sogutmaz Ozdemir, B., and Dinler Doganay, G. (2017). Time-dependent leaf proteome alterations of Brachypodium distachyon in response to drought stress. Plant Mol. Biol. 94, 609–623. doi: 10.1007/s11103-017-0628-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Tholl, D., Chen, F., Petri, J., Gershenzon, J., and Pichersky, E. (2005). Two sesquiterpene synthases are responsible for the complex mixture of sesquiterpenes emitted from Arabidopsis flowers. Plant J. 42, 757–771. doi: 10.1111/j.1365-313X.2005.02417.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, P., Nan, Z. B., Li, C. J., and Spangenberg, G. (2008). Effect of the endophyte Neotyphodium lolii on susceptibility and host physiological response of perennial ryegrass to fungal pathogens. Eur. J. Plant Pathol. 122, 593–602. doi: 10.1007/s10658-008-9329-7

CrossRef Full Text | Google Scholar

Toyomasu, T., Kagahara, T., Okada, K., Koga, J., Hasegawa, M., Mitsuhashi, W., et al. (2008). Diterpene phytoalexins are biosynthesized in and exuded from the roots of rice seedlings. Biosci. Biotech. Bioch. 72, 562–567. doi: 10.1271/bbb.70677

PubMed Abstract | CrossRef Full Text | Google Scholar

Vattekkatte, A., Garms, S., Brandt, W., and Boland, W. (2018). Enhanced structural diversity in terpenoid biosynthesis: enzymes, substrates and cofactors. Org. Biomol. Chem. 16, 348–362. doi: 10.1039/c7ob02040f

PubMed Abstract | CrossRef Full Text | Google Scholar

Vaughan, M. M., Christensen, S., Schmelz, E. A., Huffaker, A., McAuslane, H. J., Alborn, H. T., et al. (2015). Accumulation of terpenoid phytoalexins in maize roots is associated with drought tolerance. Plant Cell Environ. 38, 2195–2207. doi: 10.1111/pce.12482

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H. X., Liu, C. M., Liu, Q., and Gao, K. (2008). Three types of sesquiterpenes from rhizomes of Atractylodes lancea. Phytochemistry 69, 2088–2094. doi: 10.1016/j.phytochem.2008.04.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, N., Wu, X. L., Ku, L. X., Chen, Y. H., and Wang, W. (2016). evaluation of three protein-extraction methods for proteome analysis of maize leaf midrib, a compound tissue rich in sclerenchyma cells. Front. Plant Sci. 7:856. doi: 10.3389/fpls.2016.00856

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X. M., Yang, B., Ren, C. G., Wang, H. W., Wang, J. Y., and Dai, C. C. (2015). Involvement of abscisic acid and salicylic acid in signal cascade regulating bacterial endophyte-induced volatile oil biosynthesis in plantlets of Atractylodes lancea. Physiol. Plantarum 153, 30–42. doi: 10.1111/ppl.12236

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Dai, C. C., Cao, J. L., and Xu, D. S. (2012). Comparison of the effects of fungal endophyte Gilmaniella sp. and its elicitor on Atractylodes lancea plantlets. World J. Microb. Biot. 28, 575–584. doi: 10.1007/s11274-011-0850-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Dai, C. C., Zhao, Y. W., and Peng, Y. (2011). Fungal endophyte-induced volatile oil accumulation in Atractylodes lancea plantlets is mediated by nitric oxide, salicylic acid and hydrogen peroxide. Process Biochem. 46, 730–735. doi: 10.1016/j.procbio.2010.11.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y. W., Xu, C., Li, K., Cai, X. J., Wu, M., and Chen, G. X. (2017). Fe deficiency induced changes in rice (Oryza sativa L.) thylakoids. Environ. Sci. Pollut. R. 24, 1380–1388. doi: 10.1007/s11356-016-7900-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, G. F., Dong, L. L., Yang, J., Zhang, L. J., Xu, J., Yang, F., et al. (2018). Integrated metabolomic and transcriptomic analyses revealed the distribution of saponins in Panax notoginseng. Acta Pharm. Sin. B 8, 458–465. doi: 10.1016/j.apsb.2017.12.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, C. Q., Fang, X., Wu, X. M., Mao, Y. B., Wang, L. J., and Chen, X. Y. (2012). Transcriptional regulation of plant secondary metabolism. J. Integr. Plant Biol. 54, 703–712. doi: 10.1111/j.1744-7909.2012.01161.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, C. Q., Wu, X. M., Ruan, J. X., Hu, W. L., Mao, Y. B., Chen, X. Y., et al. (2013). Isolation and characterization of terpene synthases in cotton (Gossypium hirsutum). Phytochemistry 96, 46–56. doi: 10.1016/j.phytochem.2013.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, M. Z., Ma, M. D., Yuan, M. Q., Huang, Z. Y., Yang, W. X., Zhang, H. B., et al. (2016). Fungal endophytes as a metabolic fine-tuning regulator for wine grape. PLoS One 11:e0163186. doi: 10.1371/journal.pone.0163186

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, T., and Dai, C. C. (2013). Interactions of two endophytic fungi colonizing Atractylodes lancea and effects on the host’s essential oils. Acta Ecol. Sinica 33, 87–93. doi: 10.1016/j.chnaes.2013.01.004

CrossRef Full Text | Google Scholar

Yang, T., Du, W., Zhou, J. Y., Wang, X. M., and Dai, C. C. (2013). Effects of the symbiosis between fungal endophytes and Atractylodes lancea on rhizosphere and phyllosphere microbial communities. Symbiosis 61, 23–36. doi: 10.1007/s13199-013-0254-y

CrossRef Full Text | Google Scholar

Yuan, J., Sun, K., Deng-Wang, M. Y., and Dai, C. C. (2016a). The Mechanism of ethylene signaling induced by endophytic fungus Gilmaniella sp. AL12 mediating sesquiterpenoids biosynthesis in Atractylodes lancea. Front. Plant Sci. 7:361. doi: 10.3389/fpls.2016.00361

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, J., Zhou, J. Y., Li, X., and Dai, C. C. (2016b). The primary mechanism of endophytic fungus Gilmaniella sp AL12 promotion of plant growth and sesquiterpenoid accumulation in Atractylodes lancea. Plant Cell Tiss. Org. 125, 571–584. doi: 10.1007/s11240-016-0971-z

CrossRef Full Text | Google Scholar

Zhai, X., Luo, D., Li, X. Q., Han, T., Jia, M., Kong, Z. Y., et al. (2018). Endophyte Chaetomium globosum D38 promotes bioactive constituents accumulation and root production in Salvia miltiorrhiza. Front. Microbiol. 8:2694. doi: 10.1080/1040841X.2016.1201041

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z. B., Wu, Y. L., Gao, M. H., Zhang, J., Kong, Q., Liu, Y. A., et al. (2012). Disruption of PAMP-induced map kinase cascade by a Pseudomonas syringae effector activates plant immunity mediated by the NB-LRR protein SUMM2. Cell Host Microbe 11, 253–263. doi: 10.1016/j.chom.2012.01.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, L. P., Tian, H., Yuan, Y. F., and Wang, J. W. (2016). The influence of endophytic Penicillium oxalicum B4 on growth and artemisinin biosynthesis of in vitro propagated plantlets of Artemisia annua L. Plant Growth Regul. 80, 93–102. doi: 10.1007/s10725-016-0162-2

CrossRef Full Text | Google Scholar

Zhou, J. Y., Li, X., Zhao, D., Deng-Wang, M. Y., and Dai, C. C. (2016). Reactive oxygen species and hormone signaling cascades in endophytic bacterium induced essential oil accumulation in Atractylodes lancea. Planta 244, 699–712. doi: 10.1007/s00425-016-2536-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, J. Y., Sun, K., Chen, F., Yuan, J., Li, X., and Dai, C. C. (2018). Endophytic Pseudomonas induces metabolic flux changes that enhance medicinal sesquiterpenoid accumulation in Atractylodes lancea. Plant Physiol. Bioch. 130, 473–481. doi: 10.1016/j.plaphy.2018.07.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, L., Zuo, Z., and Chow, M. S. (2005). Danshen: an overview of its chemistry, pharmacology, pharmacokinetics, and clinical use. J. Clin. Pharmacol. 45, 1345–1359. doi: 10.1002/jcph.96

PubMed Abstract | CrossRef Full Text | Google Scholar

Zouari, I., Salvioli, A., Chialva, M., Novero, M., Miozzi, L., Tenore, G. C., et al. (2014). From root to fruit: RNA-Seq analysis shows that arbuscular mycorrhizal symbiosis may affect tomato fruit metabolism. BMC Genomics 15:221. doi: 10.1186/1471-2164-15-221

PubMed Abstract | CrossRef Full Text | Google Scholar