Next Article in Journal
Rapid and Specific Detection of the Poplar Black Spot Disease Caused by Marssonina brunnea Using Loop-Mediated Isothermal Amplification Assay
Next Article in Special Issue
Secondary Metabolites of Pseudomonas aeruginosa LV Strain Decrease Asian Soybean Rust Severity in Experimentally Infected Plants
Previous Article in Journal
Microbe Finder (MiFi®): Implementation of an Interactive Pathogen Detection Tool in Metagenomic Sequence Data
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 10
Issue 2
10.3390/plants10020251
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Transcriptomic Approach Provides Insights on the Mycorrhizal Symbiosis of the Mediterranean Orchid Limodorum abortivum in Nature

by
Rafael B. S. Valadares
1,†,
Fabio Marroni
2,3,†,
Fabiano Sillo
4,†,
Renato R. M. Oliveira
1,5,
Raffaella Balestrini
4,* and
Silvia Perotto
6,*
1
Instituto Tecnológico Vale, Rua Boaventura da Silva 955, 66050-000 Belém, Pará, Brazil
2
Dipartimento di Scienze Agroalimentari, Ambientali e Animali, Università di Udine, Via delle Scienze, I-33100 Udine, Italy
3
Istituto di Genomica Applicata, Via Linussio 51, I-33100 Udine, Italy
4
Consiglio Nazionale Delle Ricerche-Istituto per la Protezione Sostenibile Delle Piante, Viale P.A. Mattioli 25, I-10125 Torino, Italy
5
Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais, Av. Pres. Antônio Carlos, 6627, 31270-901 Belo Horizonte, Minas Gerais, Brazil
6
Dipartimento di Scienze della Vita e Biologia dei Sistemi, Università di Torino, Viale Mattioli 25, I-10125 Torino, Italy
*
Authors to whom correspondence should be addressed.
†
These authors equally contributed to the paper as first authors.
Plants 2021, 10(2), 251; https://0-doi-org.brum.beds.ac.uk/10.3390/plants10020251
Submission received: 5 January 2021 / Revised: 20 January 2021 / Accepted: 22 January 2021 / Published: 28 January 2021
(This article belongs to the Special Issue Genetic Control of Plant Interaction with Beneficial Microbes)
Browse Figures
Versions Notes

Abstract

:
The study of orchid mycorrhizal interactions is particularly complex because of the peculiar life cycle of these plants and their diverse trophic strategies. Here, transcriptomics has been applied to investigate gene expression in the mycorrhizal roots of Limodorum abortivum, a terrestrial mixotrophic orchid that associates with ectomycorrhizal fungi in the genus Russula. Our results provide new insights into the mechanisms underlying plant–fungus interactions in adult orchids in nature and in particular into the plant responses to the mycorrhizal symbiont(s) in the roots of mixotrophic orchids. Our results indicate that amino acids may represent the main nitrogen source in mycorrhizal roots of L. abortivum, as already suggested for orchid protocorms and other orchid species. The upregulation, in mycorrhizal L. abortivum roots, of some symbiotic molecular marker genes identified in mycorrhizal roots from other orchids as well as in arbuscular mycorrhiza, may mirror a common core of plant genes involved in endomycorrhizal symbioses. Further efforts will be required to understand whether the specificities of orchid mycorrhiza depend on fine-tuned regulation of these common components, or whether specific additional genes are involved.

1. Introduction

Mycorrhizas are among the most widespread plant–fungus symbioses on Earth. Different mycorrhizal types have been described in terrestrial plants, depending on the taxonomic nature of the partners and their interactions [1]. In ectomycorrhiza, the fungal hyphae remain outside the plant cell wall boundaries, whereas in endomycorrhiza—which includes arbuscular, ericoid, and orchid mycorrhiza—the symbiosis culminates in the formation of intracellular fungal structures. Endomycorrhiza formation is a complex multistep process where the exchange of molecular signals between mycorrhizal fungi and their host plants during the extra-radical phase is followed by close surface-to-surface interactions and accommodation of the fungus inside the root cells [2]. Due to their wide distribution and occurrence in economically important crops, the molecular bases of mycorrhizal fungus–host plant interactions have been investigated mainly in arbuscular mycorrhiza (AM). Plant metabolism is deeply influenced by the AM symbiosis, with many upregulated genes being related with nutrient transport, cell wall metabolism, and potential involvement in the accommodation of the AM fungus, as well as genes encoding for transcription factors [3,4]. As compared to the AM symbiosis, little is known on the metabolic changes occurring in orchid mycorrhiza (OM), despite the fact that orchids belong to one of the largest families of flowering plants, counting about 28,000 species distributed worldwide and including terrestrial species, epiphytes, and geophytes [5].
The study of OM interactions is particularly complex because of the peculiar life cycle of orchids and their diverse trophic strategies. In fact, germination of the minute endosperm-lacking seeds leads to the formation of a tuber-like structure named protocorm [6,7]. Protocorms are fully heterotrophic in most orchid species and are provided with carbon and other nutrients by the mycorrhizal fungal partners, a strategy named mycoheterotrophy [8]. Protocorms eventually give rise to leaflets and true seedlings that may grow into adult plants with green photosynthetic leaves. However, adult plants remain achlorophyllous in about 200 orchid species, and keep relying fully on fungal derived carbon [8]. Many forest orchids develop green leaves but nevertheless remain partially dependent on the fungal partner for carbon, a strategy termed mixotrophy [9]. Irrespective of the trophic strategies, all stages of the orchid life cycle are mycorrhizal in nature, and orchids are thought to be strictly dependent on their symbiotic fungal partners for growth and survival, starting from the germinating seeds through to adulthood [10]. With few exceptions, the taxonomic position and ecological role of OM fungi closely mirrors the trophic strategy of the adult host plant. Most photosynthetic orchids associate with fungi in the “rhizoctonia” complex, an assemblage of species taxonomically assigned to Agaricomycetes in the orders Cantharellales (Tulasnellaceae and Ceratobasidiaceae) and Sebacinales (Serendipitaceae). By contrast, most mycoheterotrophic and mixotrophic orchids associate with Basidiomycetes and Ascomycetes known to form ectomycorrhiza (ECM) on other host plants [11]. The association with ECM fungi allows these parasitic orchids to ultimately gain carbon from photoautotrophic plants through the common mycorrhizal network that connects the roots of ECM plants in forest communities [12,13].
Metabolic changes in OM tissues have been investigated mainly in the protocorms of photosynthetic orchid species, colonized in vitro by compatible rhizoctonia-like fungi. Early studies described in Smith and Read [1] already revealed transfer of organic carbon from the fungal partner to Dactylorhiza purpurella protocorms, and were confirmed by more recent studies in Spiranthes sinensis [14]. Transcriptomic [15,16,17] and proteomic [18,19,20] approaches have revealed the upregulation of several plant membrane transporters in mycorrhizal protocorms, including amino acid transporters and bidirectional sugar transporters of the SWEET (Sugars Will Eventually be Exported Transporters) family [15,21,22]. Interestingly, a recent metabolomic study revealed that the external mycelium of the mycorrhizal fungus Tulasnella calospora growing near the host protocorms showed significant metabolic changes, mainly in lipid metabolism [23].
Thus, orchid protocorms have been instrumental to reveal molecular changes in OM interactions. However, these tuber-like structures are anatomically and metabolically very different from the roots of adult orchids, especially in photosynthetic or mixotrophic species. We actually know very little on the metabolic and molecular changes occurring in the mycorrhizal roots of adult orchids. Although Cameron et al. [24,25] demonstrated a transfer of N and P from the fungal partner to the host plant Goodyera repens, actual mutualism is still debated even for photosynthetic orchid species, and a net flow of plant-derived fixed C to the mycorrhizal fungus has been so far demonstrated in vitro only for two photosynthetic orchid species, G. repens [24,26] and Serapias strictiflora [27].
Molecular studies on mycorrhizal roots of adult orchids are very few. A transcriptomic approach has been used to investigate changes in gene expression in adult OM roots of the epiphytic photosynthetic orchid Cymbidium hybridum, inoculated in vitro with the mycorrhizal fungus Epulorhiza repens [16]. To our knowledge, Oeceoclades maculata is the only orchid species where gene and protein expression in mycorrhizal and non-mycorrhizal roots was compared under natural conditions [28]. This tropical terrestrial orchid specifically associates with the wood decomposer Psathyrella candolleana [28], a fungus found as symbiont of mycoheterotrophic orchids [29]. Suetsugu et al. [30] also used RNAseq to investigate gene expression in mycorrhizal roots of adult orchids in nature, but their comparison of albino and green variants of Epipactis helleborine aimed to uncover the molecular mechanisms related to mycoheterotrophy, rather than to OM interactions.
Here, we used a transcriptomic approach to investigate changes in gene expression in the mycorrhizal roots of adult individuals of Limodorum abortivum grown in nature, as compared to non-mycorrhizal roots. We focused on L. abortivum because this terrestrial orchid form very thick roots with a large cortical parenchyma usually colonized by intracellular fungal coils, but non-mycorrhizal roots can be found in nature and were used as reference. In addition, previous studies on L. abortivum provide useful information on plant metabolism and fungal diversity. Stable isotope profiling [31] as well as biochemical and physiological analyses have in fact demonstrated that L. abortivum is a mixotrophic orchid with very inefficient photosynthesis and a specific association with ECM fungi in the genus Russula [32].
This study on L. abortivum expands our knowledge on the molecular bases of mycorrhizal interactions in adult orchids in nature, and addresses for the first time OM interactions in a species associated with an ECM fungus. Additionally, the comparison with a nonsterile non-mycorrhizal condition as reference likely facilitated the identification of mycorrhiza-specific plant responses over general plant responses to soil and root-associated microorganisms.

2. Results

2.1. Analysis of RNA-Seq Data

The number of raw and trimmed reads, together with the number of uniquely aligned reads, are shown in Table S1. Only reads uniquely aligned on the assembled transcriptome were used for differential expression analyses, corresponding to 4,220,880 reads (SD = 287,707) for mycorrhizal roots (LM), and 4,518,992 reads (SD = 501,399) for non-mycorrhizal roots (LS) samples. Reads assigned to fungi were on average 695,978 (SD = 400,597) in LM, and 190,961 (SD = 58,173) in LS.
A total of 30,591 transcripts tested for differential expression (DE) by DESeq2 [33] were assigned to Streptophyta. The 10 plant taxa identified by Taxon kit and represented by the highest number of transcripts are reported in Figure S1. The genome of L. abortivum has not been sequenced, but about 65% of plant transcripts matched with other orchid genera, namely Dendrobium, Phalaenopsis, Corallorhiza, and Neottia. Of the genes assigned to Streptophyta, 1677 were differentially expressed in LM and LS. In detail, when compared to LS samples, 1167 DE genes out of 1677 were upregulated in LM, while 510 were found to be downregulated (Table S2). The number of fungal transcripts identified (3094) was sensibly lower.

2.2. Transcriptomic Profile of L. abortivum Roots

Roots of L. abortivum were collected in nature where, in addition to mycorrhizal fungi, they interacted with a complex community of bacteria and fungi, as suggested by the range of fungal taxa identified in the metatranscriptome (Figure 1). The 10 most represented fungal genera covered only about 30% of the fungal reads and included both Basidiomycota (i.e., Russula, Trametes, Heterobasidion, Postia) and Ascomycota (i.e., Aspergillus, Fusarium, Exophiala). Russula was the most abundant fungal genus identified in mycorrhizal roots, confirming that this ECM fungus is the dominant symbiont of L. abortivum [32], whereas Aspergillus, together with other genera in the Ascomycota, was the most abundant genus in non-mycorrhizal roots (Figure 1). Besides Russula, other Basidiomycota mainly in the Polyporales were identified in mycorrhizal but not in non-mycorrhizal L. abortivum roots (i.e., Trametes, Heterobasidion, Stereum, Schizophyllum, and Postia), although their relationship with the plant is unknown.
When the plant transcriptomic profile was analyzed by PCA, mycorrhizal (LM) and non-mycorrhizal (LS) L. abortivum root samples were well separated along the first axis, representing 76% of variance (Figure 2). These results indicate a different plant gene expression profile in root tissues colonized by the mycorrhizal fungal symbiont(s), despite the fact that root samples in nature may display more variability than samples obtained in vitro because they are likely exposed to different environmental conditions and/or microbiota (e.g., see the more distant LM6 sample in Figure 2). Both LM and LS root samples were exposed to a rich microbiota that included non-mycorrhizal fungi, further suggesting that the separation in the PCA (Figure 2) mirrored mycorrhiza-specific plant responses in LM roots.

2.3. GO, KEGG Pathway, and CAZyme Enrichment Analyses of Plant Transcripts in Mycorrhizal L. abortivum Roots

Differentially expressed plant transcripts showed significant enrichment for several Gene Ontology (GO) terms in the three GO domains (Figure 3a). In the molecular function domain of GO, L. abortivum transcripts showed enrichment in several terms related to hydrolases and peptidases. In the cellular component, we noticed an enrichment for several terms related to synaptic membrane, main axons and similar, which of course are not found in plants. However, these components include integral membrane proteins involved in secretion and vesicle transport in many eukaryotes. For example, two t-SNARE proteins belonging to the syntaxin SYP132 family (TRINITY_DN26107_c1_g3_i10 and TRINITY_DN26107_c1_g3_i4) were significantly upregulated in LM roots (log2 FC = 4.6 and 3.6, respectively). In the biological process, several categories related to transport, and specifically transmembrane transport, were overrepresented (Figure 3a).
Molecular changes in mycorrhizal L. abortivum roots, as compared to non-mycorrhizal roots, were further investigated by examining the enrichment of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in the plant differentially expressed genes (DEGs). Seven KEGG pathways, mostly related to glycosphingolipids synthesis, nitrogen metabolism, transport, and sugar metabolism were overrepresented (Figure 3b). Differential expression of specific genes related to these metabolic pathways is described in the next paragraphs.
Interestingly, most of the enriched pathways included upregulated genes (Figure S2), whereas only three enriched pathways included downregulated genes, namely photosynthesis, oxidative phosphorylation, and phenylpropanoid biosynthesis (Figure S3). Genes not differentially expressed in L. abortivum roots showed relatively low levels of enrichment in GO categories, with only three enriched GO terms, all in the Biological Process class and with a log2 enrichment lower than 1 (Figure S4).
Plants produce a wide range of enzymes involved in carbohydrate metabolism. Many of these Carbohydrate-Active Enzymes (CAZymes; http://www.cazy.org/) are regulated during plant development and plant responses to biotic and abiotic stresses [34]. Analysis of plant DEGs in L. abortivum mycorrhizal roots (LM; Figure 4) led to the identification of several annotated CAZyme domains, with an enrichment of domains corresponding to glycoside hydrolases (GHs) and carbohydrate esterases (CEs). The enrichment analysis showed that GHs were enriched in the upregulated genes (Figure S5), and CEs were enriched in the downregulated genes (Figure S6). In particular, most CEs corresponded to enzymes acting on pectin and were strongly downregulated in LM root samples, like a pectin acetylesterase (TRINITY_DN27496_c1_g1_i8; log2 FC = −5.35) and two pectin methylesterases (TRINITY_DN37597_c0_g1_i1; log2 FC = −3.3 and TRINITY_DN27102_c2_g1_i1; log2 FC = −2.0). Two pectin methylesterase inhibitors (TRINITY_DN27850_c3_g1_i1; log2 FC = 2.3 and TRINITY_DN27102_c2_g1_i2; log2 FC = 4.4) were by contrast upregulated in LM samples. Among the GHs, the most upregulated glycoside hydrolases were involved in the fungal cell wall modification/degradation and included acidic endochitinases, as well as a putative glucan endo-1,3-beta-glucosidases (Table S2). Chitinases hydrolyze β-1,4-glycosidic bonds and are mostly regarded as part of the plant defense mechanisms induced as a response to colonization by phytopathogenic fungi. However, plant genes mainly encoding class III chitinases are expressed in arbuscule-containing cells in AM, suggesting a potential role in fungal cell wall modification during arbuscule development, as well as in the hydrolysis of chitin-derived elicitors [35,36]. A gene coding for a putative class III endochitinase was among the most upregulated genes in L. abortivum mycorrhizal roots (Table 1), in agreement with previous observations in S. vomeracea protocorms colonized by T. calospora [15].

2.4. Most Upregulated Transcripts in Mycorrhizal Roots of L. abortivum

The 25 most upregulated plant DEGs (adjusted p < 0.05) are listed in Table 1 and their sequences were further analyzed. Their log2 fold-change (FC), when compared with expression in non-mycorrhizal (LS) roots, ranged from 8.26 to 12.02. In detail, the most upregulated gene coded for cucumisin (TRINITY_DN22095_c0_g2_i2), a subtilisin-like endoprotease [37]. Other proteinases, such as cysteine proteases, aspartic proteinase, and a serine carboxypeptidase, were also among the most upregulated L. abortivum DEGs (Table 1). In particular, the secreted aspartic proteinase Constitutive Disease Resistance (CDR1), the best Blast match for the L. abortivum transcripts (Table 1), plays an important role during the immune response to bacterial pathogens in Arabidopsis [38], whereas cysteine proteases are involved in several functions, including senescence and programmed cell death [39]. In addition to these most upregulated genes, several other strongly upregulated proteinase coding genes were identified in mycorrhizal L. abortivum roots (Table 1 and Table S2).
The list of the most upregulated plant genes in mycorrhizal L. abortivum roots included genes belonging to different functional categories previously identified in AM roots and in mycorrhizal protocorms of the terrestrial orchid S. vomeracea colonized by T. calospora [15,21] as well as genes reported to be upregulated in mycorrhizal roots of O. maculata colonized by P. candolleana [28]. Among these upregulated genes were genes encoding some membrane transporters, such as two bidirectional SWEET sugar transporters (TRINITY_DN27909_c0_g1_i5 and TRINITY_DN26416_c3_g4_i1), a peptide (TRINITY_DN28361_c1_g4_i3) and a copper transporter (TRINITY_DN63579_c0_g1_i1), a gene coding for an acidic endochitinase (TRINITY_DN28711_c0_g3_i3) and a gene coding for a thaumatin-like protein (TRINITY_DN26462_c1_g2_i2). The presence of a gene encoding an early nodulin with a plastocyanin-like domain (TRINITY_DN25136_c5_g6_i3) agrees with previous observations in S. vomeracea [15], where a gene coding for a nodulin containing a plastocyanin-like domain (SvNOD1) was expressed only in coil-containing cells [40]. In addition to the top 25 upregulated plant genes in mycorrhizal L. abortivum roots, other genes are worth mentioning as they contribute to a better understanding of the molecular and cellular interactions in OM.
Notably, many genes involved in nitrogen metabolism were significantly regulated in mycorrhizal roots of L. abortivum. In particular, genes coding for a putative amino acid permease (TRINITY_DN24555_c1_g1_i4; log2 FC = 5.8) and a lysine histidine transporter 1- (LHT1; TRINITY_DN28439_c1_g5_i2; log2 FC = 4.2) were found to be upregulated, suggesting that amino acids may be a source of nitrogen for adult L. abortivum orchids. Two additional amino acid transporter genes were found to be upregulated in mycorrhizal roots of L. abortivum (Table S2), but they are unlikely to be involved in symbiotic nitrogen exchange. The first coded for an amino acid permease BAT1 (TRINITY_DN25856_c0_g2_i2; log2 FC = 2.9), a membrane transporter likely mediating the transport of gamma-aminobutyric acid from the cytosol to mitochondria [41], whereas the second coded for a cationic amino acid transporter (TRINITY_DN19554_c0_g2_i3; log2 FC = 5.4), a permease likely located on the plastid membrane and only expressed in roots in Arabidopsis [42]. Genes coding for proteins homologous to members of the NRT1/PTR family were upregulated in mycorrhizal roots of L. abortivum (Table 1 and Table S2), and one of them (TRINITY_DN26254_c6_g1_i1; log2 FC = 10.19) was among the 25 most upregulated genes.
Some genes identified in L. abortivum roots further support the hypothesis of a supply of organic nitrogen to the host root by the mycorrhizal fungus. In particular, the expression of a high affinity plant nitrate transporter (TRINITY_DN22863_c0_g1_i2) was significantly downregulated (log2 FC = −2.5) in L. abortivum mycorrhizal roots (Table S2). Different glutamine synthetase (GS) isoforms were upregulated in mycorrhizal L. abortivum roots (Table S2), and one of them (TRINITY_DN28208_c3_g1_i5; log2 FC = 8.4) was among the 25 most upregulated genes. Cytosolic GS is particularly important for nitrogen remobilization and recycling [43] and strong upregulation of a GS coding gene was already reported in S. vomeracea mycorrhizal protocorms [21].
Early experiments by Smith [44] indicate that a transfer of carbon occurs in the form of sugars, at least in orchid protocorms. In plants, several sugar transporters, including members of the SWEET family and monosaccharide transporters (MST), are responsible of sugar transport. In L. abortivum, two monosaccharide transporters similar to MST3 (TRINITY_DN27753_c2_g4_i1 and TRINITY_DN28768_c4_g4_i7; both log2 FC = 3.8) were upregulated in LM samples, together with several transporters of the SWEET family (Table 1 and Table S2).

3. Discussion

Although the last few years have seen an increase in the use of omics approaches to investigate the molecular basis of orchid mycorrhizal symbioses [23,45], information is still patchy mainly due to the complexity of this symbiosis, that involves different developmental stages as well as orchid and fungal symbionts with a range of lifestyles. The absence of Limodorum abortivum genome sequence did not make feasible the alignment of transcript short reads directly to a reference genome, thus putatively leading to an incomplete identification of genes. However, through a de novo assembly-based method, we have contributed to increase knowledge on the orchid response in a symbiotic system so far unexplored (i.e., a terrestrial mixotrophic adult orchid colonized by an ECM fungus in nature). In addition, comparison of our results with published data indicates that many plant responses identified in L. abortivum are conserved not only in other orchid species and developmental stages [16,17,28], but also in AM symbioses [2,3,4], and may mirror a common core of plant genes involved in endomycorrhizal symbioses.

3.1. Disentangling Mycorrhiza Specific and General Responses to Microbes in L. abortivum Roots

Plant-associated microbes can influence gene expression in the host plant by inducing signals that trigger plant immune responses [46] and induce a primed status that helps plants to respond to stress factors [47]. The use of nonsterile non-mycorrhizal roots (LS samples) as reference in our study likely allowed an easier identification of mycorrhiza-specific plant responses with respect to the more general plant responses to microbes. Genes whose expression was not significantly modified in LS and LM root samples were involved in some of these general plant responses (Table S2). Among them, we observed several genes involved in the jasmonic acid (JA) signaling pathway, which is involved in the mediation of plant responses to microorganisms, including the beneficial root-associated microbiota, and may lead to induced systemic resistance [48]. Among genes expressed in both LS and LM roots were also several antioxidative enzymes, pathogenesis related proteins, including chitinases mainly belonging to class V and basic chitinases, and LysM domain receptor kinases likely involved in the recognition of chitin elicitors. The only LysM domain receptor kinase gene that was differentially expressed in LM roots (TRINITY_DN28270_c1_g7_i1; log2 FC = −1.88) was actually downregulated. Another interesting example was the expression of phosphate transporters. Transfer of labeled P to mycorrhizal orchid tissues was demonstrated by Cameron et al. [25] and three transcripts coding for putative inorganic phosphate transporters of the PTH1;4 family were identified both in LS and LM samples, but none of them were differentially upregulated in mycorrhizal L. abortivum roots. Intriguingly, Zhao et al. [16] observed that a phosphate transporter was upregulated in C. hybridum roots inoculated either with the mycorrhizal fungus Epulorhiza repens (Tulasnellaceae) or with the non-mycorrhizal root endophyte Umbelopsis nana, as compared with uninoculated roots. Thus, it is possible that beneficial components of the L. abortivum microbiota may influence expression of genes involved in nutrient acquisition. Such genes would not be identified as DE genes in our transcriptomic analysis. Overall, these observations suggest that the comparison of gene expression between LM and LS samples under natural conditions may have favored the identification of mycorrhiza-specific plant responses rather that general plant responses to the associated microbiota.

3.2. Genes Potentially Involved in Plant Cell Wall Remodeling and Interface Construction during Fungal Accommodation

Pectins are an important component of the plant cell wall and pectin remodeling is a central issue to understand the fundamental role of homogalacturonan-type pectins during plant growth and plant responses to environmental changes [49]. Pectin is synthesized in a highly methyl esterified form and is de-esterified within the wall by pectin methyl esterases (PMEs), thus exposing carboxyl residues that can be cross-linked by calcium [50]. The activity of PMEs is modulated by a family of proteinaceous inhibitors known as PME inhibitors (PMEIs), and their interplay can thus strengthen or loosen the cell wall [51]. Among the CAZymes differentially expressed in L. abortivum, some PMEs were strongly downregulated and PMEIs were upregulated in mycorrhizal roots, indicating that pectin in the peripheral cell wall and/or symbiotic interface is mostly in a highly methylated form in mycorrhizal roots. As a high degree of methyl esterification (DM) leads to a softer and more pliant gel in vitro and in planta, permissive of cell growth, we may speculate that pectin softening may be related with intracellular accommodation of fungal coils in the mycorrhizal roots and/or with cell expansion. Cell wall loosening is also supported by the upregulated expression of two expansin coding genes in mycorrhizal roots (TRINITY_DN3991_c0_g1_i1; log2 FC = 5.8 and TRINITY_DN23405_c0_g1_i3; log2 FC = 4.9). Pectin remodeling also plays an important role in biotic interactions, and plants with a higher DM in the cell walls due to PMEI overexpression showed in general increased resistance to fungal and bacterial pathogens [51]. Considering the role of oligogalacturonides as signaling molecules and the fact that oligogalacturonides with different degrees of esterification are expected to be released depending on the circumstances [52], our findings raise the question of whether pectin fragments might be involved or not in eliciting plant responses during OM symbiosis, as suggested in ECM [53].
Transcripts coding for subtilisin-like serine proteases were highly upregulated in the transcriptome of L. abortivum mycorrhizal roots. Most subtilisin-like serine proteases, or subtilases (SBTs), are targeted to the plant cell wall, where they can contribute to regulate cell wall organization and structure [54]. Intriguingly, a loss-of-function mutant for a SBT essential for normal seed mucilage release had increased PME activity in the seeds, suggesting that this SBT acts as a repressor of PME activity [55]. Subtilases can also function in mutualistic and pathogenic plant–microbe interactions, as they can be modulators of the plant immune responses. Several SBT genes are induced during formation of AM and nodule symbioses in the legume Lotus japonicus [56,57], and their expression has been found to be dependent on the common symbiotic pathway [58], a signaling pathway conserved in all land plants forming endosymbiosis [59].
Two strongly upregulated genes in mycorrhizal L. abortivum roots coded for syntaxin-132 (SYP132) proteins, also reported in mycorrhizal root of other orchid species [16,28]. Syntaxins have been extensively investigated in the AM symbiosis [60] and the SYP132α isoform, in particular, was demonstrated by knockdown mutations to be essential for arbuscule formation in Medicago truncatula [60]. Given the localization of SYP132 on the perisymbiotic plant membrane surrounding functional arbuscule branches, this t-SNARE protein, with a putative role in vesicle trafficking, has been suggested to be instrumental for the formation of a functioning plant–fungus interface. Although the localization of the SYP132 proteins in OM roots is currently unknown, it is tempting to speculate that a common exocytotic pathway is shared by AM and OM symbioses. Interestingly, another upregulated gene in mycorrhizal roots of L. abortivum coded for a putative synaptotagmin, which is involved in membrane fusion and can bind to t-SNAREs proteins as well as to v-SNAREs. Synaptotagmin was identified as an important component in the formation of the symbiotic membrane interface in root nodules [61], and an O. maculata homolog gene was reported to be upregulated in colonized adult roots [28].

3.3. Nutrient Exchanges in the Mycorrhizal Roots of an Adult Orchid

Recent findings in mycoheterotrophic S. vomeracea protocorms [17] suggest that the mycorrhizal fungus T. calospora provides the host with amino acids, likely as a source of both nitrogen and carbon. Genes coding for a putative amino acid permease and a lysine histidine transporter 1 (LHT1) were also upregulated in L. abortivum mycorrhizal roots. Despite the name, LHT1 can transport a broad range of amino acids, with preference for neutral and acidic amino acids, as observed in rice [62]. Upregulated genes coding for amino acid transporters and for a LHT1 have been also reported in mycorrhizal roots of the epiphytic orchid C. hybridum [16] and of the terrestrial orchid O. maculata [28]. Overall, these observations suggest that orchids receive from their mycorrhizal fungal partners organic nitrogen in the form of amino acids irrespective of the plant developmental stage, trophic strategy of the adult plant, and type of fungal partner. The expression pattern of the S. vomeracea LHT1 in mycorrhizal protocorm cells indicates a localized expression in cells containing viable fungal hyphae, suggesting a role in amino acid transfer across an intact plant–fungus membrane interface [21]. The LjLHT1.2 gene, encoding a LHT1-type amino acid transporter, was also consistently expressed in cortical cells of AM roots, where transcripts were localized mainly in arbusculated cells but also in the noncolonized cells of the root cortex [63].
Different NRT1/PTR family members were upregulated in L. abortivum roots. Members of this protein family display sequence and structural homologies with peptide transporters involved in the uptake of di- or tripeptides [63,64] and may be involved in N transfer from the fungus to the host plant. Members of the NRT1/PTR family were also upregulated in mycorrhizal roots of C. hybridyum [16], O. maculata [28], and in S. vomeracea protocorms [21], as well as in the AM symbiosis [3], where a PTR coding gene was expressed in arbusculated cells but not in noncolonized cells of mycorrhizal roots. However, members of this transporter family have been found to translocate several other substrates, in addition to peptides, such as nitrite, chloride, glucosinolates, auxin, abscisic acid, jasmonates, and gibberellins [65]. Their role in the mycorrhizal symbiosis remains therefore uncertain.
If amino acids seem to be the main source of nitrogen for mycorrhizal orchids, their carbon backbone may also represent an important source of carbon for the host plant, especially during the juvenile developmental stages and for mycoheterotrophic and mixotrophic species. Upregulated transcripts coding for SWEET transporters were reported in mycorrhizal roots of the putative mixotrophic O. maculata [28], in mycorrhizal mycoheterotrophic S. vomeracea protocorms [21], and in the mycorrhizal roots of E. helleborine albino variants [30]. It is worth noting that SWEET transporters were not upregulated in the mycorrhizal roots of the photosynthetic orchid C. hybridum [16]. These data indicate that SWEET transporters may play important roles in mycoheterotrophic developmental stages as well as in adult plants of orchid species that rely on the mycorrhizal fungal symbiont for carbon (i.e., mycoheterotrophic and mixotrophic orchids). The role and localization of these different sugar transporters require further investigations, but it is interesting that several MST and SWEET transporters were also regulated during the parasitic interaction of plant species that directly obtain carbon from their host plants [66].

4. Materials and Methods

4.1. Sampling and Preparation of Biological Materials

Transcriptional differences between mycorrhizal and non-mycorrhizal Limodorum abortivum roots were evaluated in samples collected in a single site in the Piedmont Region in North-West Italy (Monte San Giorgio, Torino, Italy). Plants growing on this site were shown to harbor a dominant mycorrhizal symbiont phylogenetically close to the ECM basidiomycete Russula delica (site IT1B in [32]). Collected root samples were thoroughly cleaned from soil debris and rinsed under tap water. Mycorrhizal roots (LM) are easily recognizable because they are thicker and darker in color than non-mycorrhizal roots (LS). Fresh hand-made sections were prepared to verify the mycorrhizal status at increasing distances from the root apex, and root segments (about 1 cm) were flash frozen in liquid nitrogen and stored at −80 °C until use.

4.2. RNA Extraction, Library Preparation, and Sequencing

Frozen root segments from different individuals were ground in liquid nitrogen using sterile mortar and pestles. Total RNA was extracted from samples following the protocol described by Chang et al. [67]. Genomic DNA was removed using the Turbo DNA-freeTM reagent (Ambion, Austin, TX, USA), according to the manufacturer’s instructions. Three biological replicates for each condition (mycorrhizal/non-mycorrhizal) were used. The quality of the isolated RNA was determined using a Nanodrop 1000 spectrophotometer (Thermo Scientific, Waltham, MA, USA) and a Bioanalyzer 2100 (Agilent, Palo Alto, CA, USA). Samples further processed had a RIN score ranging from 6.8 to 8.7 for mycorrhizal samples and from 7.8 to 8.0 for non-mycorrhizal samples, an OD 260 to 280 nm ratio between 1.8 and 2.1, and a 28S/18S ratio within 1.5–2.
cDNA libraries were constructed using 2 µg of total RNA and following the TruSeq RNA Sample preparation v.2 kit instructions (Low throughput protocol, Illumina, Inc., San Diego, CA, USA). Briefly, total RNA samples were polyA-enriched, reverse-transcribed, and double-stranded cDNA synthesized. TruSeq adapters were attached to double-stranded cDNA. Enrichment in fragments containing TruSeq adapters on both ends were performed using PCR. The quantity and quality of the enriched libraries were validated using Nanodrop 1000 and Bioanalyzer 2100, respectively. The libraries were normalized to 10 nM of amplicons in 10 mM Tris-HCl, pH 8.5 containing 0.1% Tween 20. Barcoded amplicons (library) were spread over one Illumina HiScan 1000 lane. The TruSeq PE Cluster Kit v3-cBot-HS (Illumina, Inc., San Diego, CA, USA) was used for cluster generation using 2 pM of pooled normalized libraries on the cBOT. Sequencing was performed on the Illumina HiScan 1000 to yield paired-end reads of 2 × 101 bases using the TruSeq SBS Kit v3-HS (Illumina, Inc., San Diego, CA, USA). Reads were quality-checked with Fastqc (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/), which computes various quality metrics for the raw reads.

4.3. Bioinformatics Analysis

Adapter removal and trimming of bases with phred quality score lower than 10 were performed using bbduk (http://sourceforge.net/projects/bbmap/). Transcriptome assembly was performed using the reads from all replicates and all conditions using Trinity [68] with default parameters. Reads from each sample were then aligned against the assembled transcriptome using Hisat2 [69], generating a sample-specific transcriptome. The assembler Stringtie was used to refine the individual transcriptome based on realignment of the reads on the transcriptome, and to merge and harmonize each sample into a unique consensus transcriptome [70], on which subsequent alignment-based analysis were performed. Taxonomic classification of the reads and of the reconstructed transcripts was performed using kraken2 [71] against the nt database. Taxonkit [72] was used to specifically select and group the reads and the transcripts belonging to plants (Streptophyta) and fungi (Ascomycota or Basidiomycota). Htseq-count [73] was used to count reads mapping on the genes classified as belonging to plants or fungi. Differential expression analysis was performed separately in plant and fungal transcripts using DEseq2 [33,74]. Taxonkit [72] was used for specifically selecting and grouping the contigs belonging to plants (Streptophyta) and fungi (Ascomycota or Basidiomycota). All the downstream analysis was performed for the two groups of organisms separately. A principal component analysis (PCA) was performed on normalized read counts, to assess differences among sample replicates and between conditions.
Functional annotation of the reconstructed transcriptome was performed relying on Gene Ontology [75] and KEGG [76]. In detail, blastx [77] against UniProt database was used to classify the reconstructed transcripts; GO annotation of the UniProt database was retrieved from http://current.geneontology.org/annotations/index.html, so that each transcript mapping to a UniProt term also inherited its GO annotation. In addition, transcripts were mapped to the corresponding KEGG term using the R package KEGGREST [78]. KEGG terms were summarized to KEGG pathways using the R package KEGGREST [78]. Classification of transcripts as Carbohydrate Active Enzymes (CAZymes) was performed using dbCAN2 [79]; dbCAN2 can use three different classification tools. As recommended by these authors, only transcripts identified as CAZymes by at least two tools were considered as CAZymes. Enrichment for GO, KEGG, and CAZyme terms was assayed using Fisher’s test. Enrichment analysis was performed to test enrichment in all differentially expressed genes. Additional analyses were performed to test enrichment in upregulated and downregulated genes separately, and on genes showing no differential expression.
Among the differentially expressed plant genes (adjusted p < 0.05), the 25 showing the highest log2 fold-change were further characterized through Blastx against the Refseq database and conserved and functional motifs were considered.

5. Conclusions

We investigated the transcriptomic profile of mycorrhizal roots of the mixotrophic orchid L. abortivum grown in nature. Comparison with available gene expression data in three other orchid species (S. vomeracea, O. maculata, and C. hybridum) indicate that many molecular markers are shared in mycorrhizal orchid tissues despite the different plant trophic strategy (photosynthetic, mycoheterotrophic or mixotrophic), developmental stage (protocorms or adult plants), habitat (tropical or temperate, terrestrial or epiphytic), and type of fungal partner (rhizoctonia-like, ECM or wood decomposing fungi). Many molecular markers (such as components involved in vesicle trafficking, amino acid transporters, and subtilases) were also reported in AM roots (and even legume nodule) symbioses. Thus, they may represent components of a common symbiotic toolkit shared by different endosymbioses. The analysis of a wider range of orchid species with diverse trophic strategies may help in understanding whether the OM specificities depend on fine-tuned regulation of common cell components, or whether specific additional genes are involved.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2223-7747/10/2/251/s1, Figure S1: Percentage of Streptophyta reads assigned to different genera, Figure S2: Log2 enrichment of GO terms and KEGG pathways in genes upregulated in Limodorum abortivum mycorrhizal roots. Left panel: log2 enrichment of the most significant terms for each GO category. MF = molecular function, CC = cellular component, BP = biological process (max 10 terms per category). Right panel: log2 enrichment of KEGG pathways, Figure S3: Log2 enrichment of GO terms and KEGG pathways in genes downregulated in Limodorum abortivum mycorrhizal roots. Left panel: log2 enrichment of the 10 most significant terms for each GO category. MF = molecular function, CC = cellular component, BP = biological process (max 10 terms per category). Right panel: log2 enrichment of KEGG pathways, Figure S4: Log2 enrichment of GO terms and KEGG pathways in genes showing no variation of expression in Limodorum abortivum mycorrhizal roots compared to the control. Left panel: log2 enrichment of the most significant terms for each GO category. MF = molecular function, CC = cellular component, BP = biological process (max 10 terms per category). Right panel: log2 enrichment of KEGG pathways, Figure S5: CAZyme families enrichment in Limodorum abortivum genes upregulated in mycorrhizal roots. GH: glycoside hydrolases; GT: glycosyl transferases; CE: carbohydrate esterases; AA: auxiliary activities. Significant values according to Fisher’s test (p < 0.05) are indicated with an asterisk (*), Figure S6: CAZyme families enrichment in Limodorum abortivum genes downregulated in mycorrhizal roots. GH: glycoside hydrolases; GT: glycosyl transferases; CE: carbohydrate esterases; AA: auxiliary activities. Significant values according to Fisher’s test (p < 0.05) are indicated with an asterisk (*), Table S1: Summary statistics of reads production and alignments, Table S2: Complete list of plant genes in the mycorrhizal root transcriptome of L. abortivum (excel file).

Author Contributions

Conceptualization, R.B.S.V., S.P.; methodology, R.B.S.V., F.M., F.S.; formal analysis, R.B.S.V., F.M., F.S., R.R.M.O.; data curation, R.B.S.V., F.M., F.S., R.R.M.O.; writing—original draft preparation, F.M., F.S., R.B., S.P.; writing—review and editing, R.B.S.V., F.M., F.S., R.R.M.O., R.B., S.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Sequencing reads and read counts per gene are available in GEO (GSE159700). Fasta file of raw transcriptome assembly, gtf of refined transcriptome assembly, and results of DE analysis are available on Open Science Foundation, with DOI 10.17605/OSF.IO/35XUP, and can be accessed at https://osf.io/35xup/. Shell scripts and functions used in this work are freely available at https://github.com/fabiomarroni/limodorum.

Acknowledgments

The authors thank Enrico Ercole and Samuele Voyron for their help with sample collection and processing. RV was supported by the São Paulo Research Foundation FAPESP (grant No. 2012/50266-4).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Smith, S.E.; Read, D.J. Mycorrhizal Symbiosis; Academic Press: Cambridge, UK, 2008; p. 800. [Google Scholar]
  2. Genre, A.; Lanfranco, L.; Perotto, S.; Bonfante, P. Unique and common traits in mycorrhizal symbioses. Nat. Rev. Microbiol. 2020, 18, 649–660. [Google Scholar] [CrossRef]
  3. Guether, M.; Balestrini, R.; Hannah, M.A.; Udvardi, M.K.; Bonfante, P. Genome-wide reprogramming of regulatory networks, transport, cell wall and membrane biogenesis during arbuscular mycorrhizal symbiosis in Lotus japonicus. New Phytol. 2009, 182, 200–212. [Google Scholar] [CrossRef]
  4. Balestrini, R.; Bonfante, P. Cell wall remodeling in mycorrhizal symbiosis: A way towards biotrophism. Front. Plant Sci. 2014, 5, 237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Christenhusz, M.J.M.; Byng, J.W. The number of known plants species in the world and its annual increase. Phytotaxa 2016, 261, 201–217. [Google Scholar] [CrossRef] [Green Version]
  6. Arditti, J. Fundamentals of Orchid Biology; John Wiley & Sons: New York, NY, USA, 1992; p. 691. [Google Scholar]
  7. Rasmussen, H.N. Recent developments in the study of orchid mycorrhiza. Plant Soil 2002, 244, 149–163. [Google Scholar] [CrossRef]
  8. Hynson, N.; Madsen, T.; Selosse, M.-A.; Adam, I.; Ogura-Tsujita, Y.; Roy, M.; Gebauer, G. The physiological ecology of mycoheterotrophy. In Mycoheterotrophy, the Biology of Plants Living on Fungi; Merckx, V.S.F.T., Ed.; Springer: New York, NY, USA, 2013; pp. 297–344. [Google Scholar]
  9. Selosse, M.-A.; Roy, M. Green plants that feed on fungi, facts and questions about mixotrophy. Trends Plant Sci. 2009, 14, 64–70. [Google Scholar] [CrossRef]
  10. Dearnaley, J.D.; Martos, F.; Selosse, M.A. 12 orchid mycorrhizas: Molecular ecology, physiology, evolution and conservation aspects. In Fungal Associations, the Mycota; Hock, B., Ed.; Springer: Berlin/Heidelberg, Germany, 2012; Volume 9, pp. 207–230. [Google Scholar]
  11. Bonnardeaux, Y.; Brundrett, M.; Batty, A.; Dixon, K.; Koch, J.; Sivasithamparam, K. Diversity of mycorrhizal fungi of terrestrial orchids: Compatibility webs, brief encounters, lasting relationships and alien invasions. Mycol. Res. 2007, 111, 51–61. [Google Scholar] [CrossRef]
  12. Simard, S.W.; Perry, D.A.; Jones, M.D.; Myrold, D.D.; Durall, D.M.; Molina, R. Net transfer of carbon between tree species with shared ectomycorrhizal fungi. Nature 1997, 388, 579–582. [Google Scholar] [CrossRef]
  13. McKendrick, S.L.; Leake, J.R.; Read, D.J. Symbiotic germination and development of myco-heterotrophic plants in nature: Transfer of carbon from ectomycorrhizal Salix repens and Betula pendula to the orchid Corallorhiza trifida through shared hyphal connections. New Phytol. 2000, 145, 539–548. [Google Scholar] [CrossRef]
  14. Kuga, Y.; Sakamoto, N.; Yurimoto, H. Stable isotope cellular imaging reveals that both live and degenerating fungal pelotons transfer carbon and nitrogen to orchid protocorms. New Phytol. 2014, 202, 594–605. [Google Scholar] [CrossRef]
  15. Perotto, S.; Benetti, A.; Sillo, F.; Ercole, E.; Rodda, M.; Girlanda, M.; Balestrini, R. Gene expression in mycorrhizal orchid protocorms suggests a friendly plant–fungus relationship. Planta 2014, 239, 1337–1349. [Google Scholar] [CrossRef] [PubMed]
  16. Zhao, X.; Zhang, J.; Chen, C.; Yang, J.; Zhu, H.; Liu, M.; Lv, F. Deep sequencing-based comparative transcriptional profiles of Cymbidium hybridum roots in response to mycorrhizal and non-mycorrhizal beneficial fungi. BMC Genom. 2014, 15, 747. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Fochi, V.; Chitarra, W.; Kohler, A.; Voyron, S.; Singan, V.R.; Lindquist, E.A.; Barry, K.W.; Girlanda, M.; Grigoriev, I.V.; Martin, F.; et al. Fungal and plant gene expression in the Tulasnella calospora–Serapias vomeracea symbiosis provides clues about nitrogen pathways in orchid mycorrhizas. New Phytol. 2017, 213, 365–379. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Valadares, R.B.S.; Perotto, S.; Santos, E.C.; Lambais, M.R. Proteome changes in Oncidium sphacelatum (Orchidaceae) at different trophic stages of symbiotic germination. Mycorrhiza 2014, 24, 349–360. [Google Scholar] [CrossRef] [PubMed]
  19. López-Chávez, M.Y.; Guillén-Navarro, K.; Bertolini, V.; Encarnación, S.; Hernández-Ortiz, M.; Sánchez-Moreno, I.; Damon, A. Proteomic and morphometric study of the in vitro interaction between Oncidium sphacelatum Lindl. (Orchidaceae) and Thanatephorus sp. RG26 (Ceratobasidiaceae). Mycorrhiza 2016, 26, 353–365. [Google Scholar] [CrossRef] [PubMed]
  20. Chen, J.; Liu, S.S.; Kohler, A.; Yan, B.; Luo, H.M.; Chen, X.M.; Guo, S.X. iTRAQ and RNA-Seq analyses provide new insights into regulation mechanism of symbiotic germination of Dendrobium officinale seeds (Orchidaceae). Proteome Res. 2017, 16, 2174–2187. [Google Scholar] [CrossRef]
  21. Fochi, V.; Falla, N.; Girlanda, M.; Perotto, S.; Balestrini, R. Cell-specific expression of plant nutrient transporter genes in orchid mycorrhizae. Plant Sci. 2017, 263, 39–45. [Google Scholar] [CrossRef]
  22. Miura, C.; Yamaguchi, K.; Miyahara, R.; Yamamoto, T.; Fuji, M.; Yagame, T.; Imaizumi-Anraku, H.; Yamato, M.; Shigenobu, S.; Kaminaka, H. The mycoheterotrophic symbiosis between orchids and mycorrhizal fungi possesses major components shared with mutualistic plant-mycorrhizal symbioses. Mol. Plant Microbe Interact. 2018, 31, 1032–1047. [Google Scholar] [CrossRef] [Green Version]
  23. Ghirardo, A.; Fochi, V.; Lange, B.; Witting, M.; Schnitzler, J.P.; Perotto, S.; Balestrini, R. Metabolomic adjustments in the orchid mycorrhizal fungus Tulasnella calospora during symbiosis with Serapias vomeracea. New Phytol. 2020, 228, 1939–1952. [Google Scholar] [CrossRef]
  24. Cameron, D.D.; Leake, J.R.; Read, D.J. Mutualistic mycorrhiza in orchids, evidence from plant-fungus carbon and nitrogen transfers in the green-leaved terrestrial orchid Goodyera repens. New Phytol. 2006, 171, 405–416. [Google Scholar] [CrossRef]
  25. Cameron, D.D.; Johnson, I.; Leake, J.R.; Read, D.J. Mycorrhizal acquisition of inorganic phosphorus by the green-leaved terrestrial orchid Goodyera repens. Ann. Bot. 2007, 99, 831–834. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Cameron, D.D.; Johnson, I.; Read, D.J.; Leake, J.R. Giving and receiving, measuring the carbon cost of mycorrhizas in the green orchid Goodyera repens. New Phytol. 2008, 180, 176–184. [Google Scholar] [CrossRef] [PubMed]
  27. Látalová, K.; Baláž, M. Carbon nutrition of mature green orchid Serapias strictiflora and its mycorrhizal fungus Epulorhiza sp. Biol. Plant 2010, 54, 97–104. [Google Scholar] [CrossRef]
  28. Valadares, R.B.S.; Perotto, S.; Lucheta, A.R.; Santos, E.C.; Oliveira, R.M.; Lambais, M.R. Proteomic and transcriptomic analyses indicate metabolic changes and reduced defense responses in mycorrhizal rootss of Oeceoclades maculata (Orchidaceae) collected in nature. J. Fungi 2020, 6, 148. [Google Scholar] [CrossRef] [PubMed]
  29. Bayman, P.; Mosquera-Espinosa, A.T.; Saladini-Aponte, C.M.; Hurtado-Guevara, N.C.; Viera-Ruiz, N.L. Age-dependent mycorrhizal specificity in an invasive orchid Oeceoclades maculata. Am. J. Bot. 2016, 103, 1880–1889. [Google Scholar] [CrossRef] [Green Version]
  30. Suetsugu, K.; Yamato, M.; Miura, C.; Yamaguchi, K.; Takahashi, K.; Ida, Y.; Shigenobu, S.; Kaminaka, H. Comparison of green and albino individuals of the partially mycoheterotrophic orchid Epipactis helleborine on molecular identities of mycorrhizal fungi; nutritional modes and gene expression in mycorrhizal roots. Mol. Ecol. 2017, 26, 1652–1669. [Google Scholar] [CrossRef]
  31. Gebauer, G.; Meyer, M. 15N and 13C natural abundance of autotrophic and myco-heterotrophic orchids provides insight into nitrogen and carbon gain from fungal association. New Phytol. 2003, 160, 209–223. [Google Scholar] [CrossRef]
  32. Girlanda, M.; Selosse, M.A.; Cafasso, D.; Brilli, F.; Delfine, S.; Fabbian, R.; Ghignone, S.; Pinelli, P.; Segreto, R.; Loreto, F.; et al. Inefficient photosynthesis in the Mediterranean orchid Limodorum abortivum is mirrored by specific association to ectomycorrhizal Russulaceae. Mol. Ecol. 2006, 15, 491–504. [Google Scholar] [CrossRef]
  33. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [Green Version]
  34. Lombard, V.; Golaconda Ramulu, H.; Drula, E.; Coutinho, P.M.; Henrissat, B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014, 42, 490–495. [Google Scholar] [CrossRef] [Green Version]
  35. Hogekamp, C.; Arndt, D.; Pereira, P.A.; Becker, J.D.; Hohnjec, N.; Küster, H. Laser microdissection unravels cell-type-specific transcription in arbuscular mycorrhizal roots, including CAAT-box transcription factor gene expression correlating with fungal contact and spread. Plant Physiol. 2011, 157, 2023–2043. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Bonanomi, A.; Wiemken, A.; Boller, T.; Salzer, P. Local induction of a mycorrhiza-specific class III chitinase gene in cortical root cells of Medicago truncatula containing developing or mature arbuscules. Plant Biol. 2001, 3, 194–200. [Google Scholar] [CrossRef]
  37. Xu, Y.; Wang, S.; Li, L.; Sahu, S.K.; Petersen, M.; Liu, X.; Melkonian, M.; Zhang, G.; Liu, H. Molecular evidence for origin, diversification and ancient gene duplication of plant subtilases (SBTs). Sci. Rep. 2019, 9, 12485. [Google Scholar] [CrossRef] [PubMed]
  38. Xia, Y.; Suzuki, H.; Borevitz, J.; Blount, J.; Guo, Z.; Patel, K.; Dixon, R.A.; Lamb, C. An extracellular aspartic protease functions in Arabidopsis disease resistance signaling. EMBO J. 2004, 23, 980–988. [Google Scholar] [CrossRef]
  39. Liu, H.; Hu, M.; Wang, Q.; Cheng, L.; Zhang, Z. Role of papain-like cysteine proteases in plant development. Front. Plant Sci. 2018, 9, 1717. [Google Scholar] [CrossRef] [Green Version]
  40. Balestrini, R.; Nerva, L.; Sillo, F.; Girlanda, M.; Perotto, S. Plant and fungal gene expression in mycorrhizal protocorms of the orchid Serapias vomeracea colonized by Tulasnella calospora. Plant Signal. Behav. 2014, 9, e977707. [Google Scholar] [CrossRef] [Green Version]
  41. Michaeli, S.; Fait, A.; Lagor, K.; Nunes-Nesi, A.; Grillich, N.; Yellin, A.; Bar, D.; Khan, M.; Fernie, A.R.; Turano, F.J.; et al. A mitochondrial GABA permease connects the GABA shunt and the TCA cycle, and is essential for normal carbon metabolism. Plant J. 2011, 67, 485–498. [Google Scholar] [CrossRef]
  42. Su, Y.H.; Frommer, W.B.; Ludewig, U. Molecular and functional characterization of a family of amino acid transporters from Arabidopsis. Plant Physiol. 2004, 136, 3104–3113. [Google Scholar] [CrossRef] [Green Version]
  43. Bernard, S.M.; Habash, D.Z. The importance of cytosolic glutamine synthetase in nitrogen assimilation and recycling. New Phytol. 2009, 182, 608–620. [Google Scholar] [CrossRef]
  44. Smith, S.E. Carbohydrate translocation in orchid mycorrhizas. New Phytol. 1967, 66, 371–378. [Google Scholar] [CrossRef]
  45. Yeh, C.-M.; Chung, K.; Liang, C.-K.; Tsai, W.-C. New insights into the symbiotic relationship between orchids and fungi. Appl. Sci. 2019, 9, 585. [Google Scholar] [CrossRef] [Green Version]
  46. Rodriguez, P.A.; Rothballer, M.; Chowdhury, S.P.; Nussbaumer, T.; Gutjahr, C.; Falter-Braun, P. Systems biology of plant-microbiome interactions. Mol. Plant 2019, 12, 804–821. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Alagna, F.; Balestrini, R.; Chitarra, W.; Marsico, A.D.; Nerva, L. Getting ready with the priming: Innovative weapons against biotic and abiotic crop enemies in a global changing scenario. In Priming-Mediated Stress and Cross-Stress Tolerance in Crop Plants; Hossain, M.A., Liu, F., Burritt, D., Fujita, M., Huang, B., Eds.; Academic Press: Cambridge, UK, 2020; pp. 35–56. [Google Scholar]
  48. Van Wees, S.C.M.; Van der Ent, S.; Pieterse, C.M.J. Plant immune responses triggered by beneficial microbes. Curr. Opin. Plant Biol. 2008, 11, 443–448. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Sénéchal, F.; Wattier, C.; Rustérucci, C.; Pelloux, J. Homogalacturonan-modifying enzymes: Structure, expression, and roles in plants. J. Exp. Bot. 2014, 65, 5125–5160. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Lionetti, V.; Cervone, F.; Bellincampi, D. Methyl esterification of pectin plays a role during plant–pathogen interactions and affects plant resistance to diseases. J. Plant Physiol. 2012, 169, 1623–1630. [Google Scholar] [CrossRef] [PubMed]
  51. Levesque-Tremblay, G.; Pelloux, J.; Braybrook, S.A.; Müller, K. Tuning of pectin methylesterification: Consequences for cell wall biomechanics and development. Planta 2015, 242, 791–811. [Google Scholar] [CrossRef]
  52. Ferrari, S.; Savatin, D.V.; Sicilia, F.; Gramegna, G.; Cervone, F.; De Lorenzo, G. Oligogalacturonides: Plant damage-associated molecular patterns and regulators of growth and development. Front. Plant Sci. 2013, 4, 49. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Sillo, F.; Fangel, J.U.; Henrissat, B.; Faccio, A.; Bonfante, P.; Martin, F.; Willats, W.G.T.; Balestrini, R. Understanding plant cell-wall remodelling during the symbiotic interaction between Tuber melanosporum and Corylus avellana using a carbohydrate microarray. Planta 2016, 244, 347–359. [Google Scholar] [CrossRef]
  54. Schaller, A.; Stintzi, A.; Rivas, S.; Serrano, I.; Chichkova, N.V.; Vartapetian, A.B.; Martínez, D.; Guiamét, J.J.; Sueldo, D.J.; van der Hoorn, R.A.L.; et al. From structure to function—A family portrait of plant subtilases. New Phytol. 2018, 218, 901–915. [Google Scholar] [CrossRef]
  55. Rautengarten, C.; Usadel, B.; Neumetzler, L.; Hartmann, J.; Büssis, D.; Altmann, T. A subtilisin-like serine protease essential for mucilage release from Arabidopsis seed coats. Plant J. 2008, 54, 466–480. [Google Scholar] [CrossRef]
  56. Kistner, C.; Winzer, T.; Pitzschke, A.; Mulder, L.; Sato, S.; Kaneko, T.; Tabata, S.; Sandal, N.; Stougaard, J.; Webb, K.J.; et al. Seven Lotus japonicus genes required for transcriptional reprogramming of the root during fungal and bacterial symbiosis. Plant Cell 2005, 17, 2217–2229. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Takeda, N.; Sato, S.; Asamizu, E.; Tabata, S.; Parniske, M. Apoplastic plant subtilases support arbuscular mycorrhiza development in Lotus japonicus. Plant J. 2009, 58, 766–777. [Google Scholar] [CrossRef] [PubMed]
  58. Takeda, N.; Haage, K.; Sato, S.; Tabata, S.; Parniske, M. Activation of a Lotus japonicus subtilase gene during arbuscular mycorrhiza is dependent on the common symbiosis genes and two cis-active promoter regions. Mol. Plant Microbe Interact. 2011, 24, 662–670. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Radhakrishnan, G.V.; Keller, J.; Rich, M.K.; Vernié, T.; Mbadinga Mbadinga, D.L.; Vigneron, N.; Cottret, L.; San Clemente, H.; Libourel, C.; Cheema, J.; et al. An ancestral signalling pathway is conserved in intracellular symbioses-forming plant lineages. Nat. Plants 2020, 6, 280–289. [Google Scholar] [CrossRef] [PubMed]
  60. Huisman, R.; Hontelez, J.; Bisseling, T.; Limpens, E. SNARE Complexity in Arbuscular Mycorrhizal Symbiosis. Front. Plant Sci. 2020, 11, 354. [Google Scholar] [CrossRef] [PubMed]
  61. Gavrin, A.; Kulikova, O.; Bisseling, T.; Fedorova, E.E. Interface symbiotic membrane formation in root nodules of Medicago truncatula: The role of synaptotagmins MtSyt1, MtSyt2 and MtSyt3. Front. Plant Sci. 2017, 8, 201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Guo, N.; Hu, J.; Yan, M.; Qu, H.; Luo, L.; Tegeder, M.; Xu, G. Oryza sativa Lysine-Histidine-type Transporter 1 functions in root uptake and root-to-shoot allocation of amino acids in rice. Plant J. 2020, 103, 395–411. [Google Scholar] [CrossRef]
  63. Guether, M.; Volpe, V.; Balestrini, R.; Requena, N.; Wipf, D.; Bonfante, P. LjLHT1.2-a mycorrhiza-inducible plant amino acid transporter from Lotus japonicus. Biol. Fert. Soils 2011, 47, 925. [Google Scholar] [CrossRef] [Green Version]
  64. Léran, S.; Varala, K.; Boyer, J.C.; Chiurazzi, M.; Crawford, N.; Daniel-Vedele, F.; David, L.; Dickstein, R.; Fernandez, E.; Forde, B.; et al. A unified nomenclature of nitrate transporter 1/peptide transporter family members in plants. Trends Plant Sci. 2014, 19, 5–9. [Google Scholar]
  65. Corratgé-Faillie, C.; Lacombe, B. Substrate (un)specificity of Arabidopsis NRT1/PTR family (NPF) proteins. J. Exp. Bot. 2017, 68, 3107–3113. [Google Scholar] [CrossRef]
  66. Misra, V.A.; Wafula, E.K.; Wang, Y.; de Pamphilis, C.W.; Timko, M.P. Genome-wide identification of MST, SUT and SWEET family sugar transporters in root parasitic angiosperms and analysis of their expression during host parasitism. BMC Plant Biol. 2019, 19, 196. [Google Scholar] [CrossRef] [PubMed]
  67. Chang, S.; Puryear, J.; Cairney, J. A simple and efficient method for isolating RNA from pine trees. Plant Mol. Biol. Rep. 1993, 11, 113–116. [Google Scholar] [CrossRef]
  68. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotech. 2011, 29, 644–652. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Daehwan, K.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar]
  70. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.C.; Mendell, J.T.; Salzberg, S.L. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotech. 2015, 33, 290–295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Wood, D.E.; Salzberg, S.L. Kraken: Ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 2014, 15, R46. [Google Scholar] [CrossRef] [Green Version]
  72. Shen, W.; Xiong, J. TaxonKit: A cross-platform and efficient NCBI taxonomy toolkit. BioRxiv 2019, 513523. [Google Scholar] [CrossRef] [Green Version]
  73. Anders, S.; Pyl, P.T.; Huber, W. HTSeq—A Python framework to work with high-throughput sequencing data. Bioinformatics 2015, 31, 166–169. [Google Scholar] [CrossRef]
  74. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol. 2010, 11, R106. [Google Scholar] [CrossRef] [Green Version]
  75. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene Ontology tool for the unification of biology. Nat. Genet. 2011, 25, 25–29. [Google Scholar] [CrossRef] [Green Version]
  76. Kanehisa, M.; Sato, Y.; Kawashima, M.; Furumichi, M.; Tanabe, M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016, 44, 457–462. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  78. Tenenbaum, D. KEGGREST: Client-Side REST Access to KEGG; R Package Version 1.28.0. 2020. Available online: https://bioconductor.org/packages/release/bioc/html/KEGGREST.html (accessed on 5 January 2021).
  79. Zhang, H.; Yohe, T.; Huang, L.; Entwistle, S.; Wu, P.; Yang, Z.; Busk, P.K.; Xu, Y.; Yin, Y. DbCAN2: A Meta server for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2018, 46, W95–W101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Plants 10 00251 g001 550
Figure 1. Percentage of fungal reads assigned to different genera in the roots of Limodorum abortivum. (a) Percentage of fungal reads in mycorrhizal (LM) root samples. (b) Percentage of fungal reads in non-mycorrhizal (LS) root samples. The error bar represents the 95% confidence interval of abundance estimate based on three replicates per condition.
Figure 1. Percentage of fungal reads assigned to different genera in the roots of Limodorum abortivum. (a) Percentage of fungal reads in mycorrhizal (LM) root samples. (b) Percentage of fungal reads in non-mycorrhizal (LS) root samples. The error bar represents the 95% confidence interval of abundance estimate based on three replicates per condition.
Plants 10 00251 g001
Plants 10 00251 g002 550
Figure 2. PCA analysis of the complete plant transcriptome in the three biological replicates of mycorrhizal (LM) and non-mycorrhizal (LS) root samples of Limodorum abortivum.
Figure 2. PCA analysis of the complete plant transcriptome in the three biological replicates of mycorrhizal (LM) and non-mycorrhizal (LS) root samples of Limodorum abortivum.
Plants 10 00251 g002
Plants 10 00251 g003 550
Figure 3. Enrichment in GO terms and KEGG pathways in the plant transcriptome of Limodorum abortivum mycorrhizal roots: (a) log2 enrichment of the 10 most significant terms for each GO category. MF = molecular function, CC = cellular component, BP = biological process; (b) enrichment of KEGG pathways in the plant differentially expressed (DE) genes.
Figure 3. Enrichment in GO terms and KEGG pathways in the plant transcriptome of Limodorum abortivum mycorrhizal roots: (a) log2 enrichment of the 10 most significant terms for each GO category. MF = molecular function, CC = cellular component, BP = biological process; (b) enrichment of KEGG pathways in the plant differentially expressed (DE) genes.
Plants 10 00251 g003
Plants 10 00251 g004 550
Figure 4. CAZyme families enrichment in Limodorum abortivum genes differentially expressed in mycorrhizal roots. GH: glycoside hydrolases; GT: glycosyl transferases; CE: carbohydrate esterases; AA: auxiliary activities. Significant values according to Fisher’s test (p < 0.05) are indicated with an asterisk (*).
Figure 4. CAZyme families enrichment in Limodorum abortivum genes differentially expressed in mycorrhizal roots. GH: glycoside hydrolases; GT: glycosyl transferases; CE: carbohydrate esterases; AA: auxiliary activities. Significant values according to Fisher’s test (p < 0.05) are indicated with an asterisk (*).
Plants 10 00251 g004
Table
Table 1. List of the 25 most upregulated plant DEGs identified in this study.
Table 1. List of the 25 most upregulated plant DEGs identified in this study.
Sequence IDlog2
Fold Change		Putative Function
(Blastx, Refseq)		Plant SpeciesQuery
Cov.
(%)		Ident. (%)E-Value
TRINITY_DN22095_c0_g2_i212.02Cucumisin-likeDendrobium catenatum8569.680.0
TRINITY_DN27918_c0_g1_i611.02Senescence-specific cysteine protease SAG39-likePopulus euphratica7279.522e-88
TRINITY_DN21619_c0_g1_i210.934,5-DOPA dioxygenase extradiolDendrobium catenatum5386.981e-98
TRINITY_DN24828_c2_g2_i310.54Putative beta-glucosidasePopulus alba4682.502e-87
TRINITY_DN26254_c6_g1_i110.19Protein NRT1/ PTR FAMILY 8.2-likeDendrobium catenatum7565.880.0
TRINITY_DN27041_c6_g2_i110.00Uncharacterized proteinDendrobium catenatum5252.403e-57
TRINITY_DN28711_c0_g3_i39.54Acidic endochitinase-likeDendrobium catenatum6280.411e-166
TRINITY_DN26462_c1_g2_i29.71Thaumatin-like protein 1bPhalaenopsis equestris5183.563e-115
TRINITY_DN25136_c5_g6_i39.53Early nodulin-like protein 2 (plastocyanin-like)Phalaenopsis equestris3276.702e-50
TRINITY_DN25901_c2_g1_i99.44Putative calcium-binding protein CML31Dendrobium catenatum3367.132e-56
TRINITY_DN22366_c0_g1_i19.22Serine carboxypeptidase II-3Dendrobium concatenatum7175.880.0
TRINITY_DN25614_c3_g3_i89.16Oryzain gamma chain-like (putative cysteine peptidase)Dendrobium catenatum4383.130.0
TRINITY_DN38600_c0_g1_i19.06Aspartic proteinase CDR1-likeDendrobium catenatum8879.730.0
TRINITY_DN28212_c0_g1_i19.02Aspartic proteinase CDR1-likePhalaenopsis equestris8880.280.0
TRINITY_DN29286_c2_g2_i78.99Alpha-humulene synthase-likePhalaenopsis equestris8267.582e-82
TRINITY_DN26369_c2_g2_i28.61Uncharacterized proteinPhalaenopsis equestris5958.793e-47
TRINITY_DN63579_c0_g1_i18.55Copper transporter 6-likeDendrobium catenatum3872.365e-47
TRINITY_DN27909_c0_g1_i58.54Bidirectional sugar transporter SWEET4-likePhalaenopsis equestris3295.243e-29
TRINITY_DN23758_c0_g1_i18.50Heparanase-like protein 3Dendrobium catenatum7068.652e-177
TRINITY_DN19642_c0_g1_i18.50Glucan endo-1,3-beta-glucosidaseDendrobium catenatum7872.677e-167
TRINITY_DN20494_c0_g1_i68.49EG45-like domain containing proteinPopulus trichocarpa3049.492e-20
TRINITY_DN28208_c3_g1_i58.39Glutamine synthetase nodule isozymePhoenix dactylifera3287.781e-172
TRINITY_DN25633_c3_g1_i108.28Beta-hexosaminidase 3Dendrobium catenatum6387.210.0
TRINITY_DN26416_c3_g4_i18.27Bidirectional sugar transporter SWEET4-likeDendrobium catenatum8186.791e-56
TRINITY_DN28361_c1_g4_i38.26Oligopeptide transporter 1-likeDendrobium catenatum6976.145e-123
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Valadares, R.B.S.; Marroni, F.; Sillo, F.; Oliveira, R.R.M.; Balestrini, R.; Perotto, S. A Transcriptomic Approach Provides Insights on the Mycorrhizal Symbiosis of the Mediterranean Orchid Limodorum abortivum in Nature. Plants 2021, 10, 251. https://0-doi-org.brum.beds.ac.uk/10.3390/plants10020251

AMA Style

Valadares RBS, Marroni F, Sillo F, Oliveira RRM, Balestrini R, Perotto S. A Transcriptomic Approach Provides Insights on the Mycorrhizal Symbiosis of the Mediterranean Orchid Limodorum abortivum in Nature. Plants. 2021; 10(2):251. https://0-doi-org.brum.beds.ac.uk/10.3390/plants10020251

Chicago/Turabian Style

Valadares, Rafael B. S., Fabio Marroni, Fabiano Sillo, Renato R. M. Oliveira, Raffaella Balestrini, and Silvia Perotto. 2021. "A Transcriptomic Approach Provides Insights on the Mycorrhizal Symbiosis of the Mediterranean Orchid Limodorum abortivum in Nature" Plants 10, no. 2: 251. https://0-doi-org.brum.beds.ac.uk/10.3390/plants10020251

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop