Original Research ARTICLE
The Hydrophobin-Like OmSSP1 May Be an Effector in the Ericoid Mycorrhizal Symbiosis
- 1Department of Life Sciences and Systems Biology, University of Turin, Turin, Italy
- 2INRA (Institut National de la Recherche Agronomique), UMR 1136 INRA-Université de Lorraine Interactions Arbres/Microorganismes, Laboratoire d'Excellence ARBRE, Centre INRA-Lorraine, Champenoux, France
- 3Université de Lorraine, UMR 1136 INRA-Université de Lorraine Interactions Arbres/Microorganismes, Laboratoire d'Excellence ARBRE, Faculté des Sciences et Technologies, Vandoeuvre les Nancy, France
Mutualistic and pathogenic plant-colonizing fungi use effector molecules to manipulate the host cell metabolism to allow plant tissue invasion. Some small secreted proteins (SSPs) have been identified as fungal effectors in both ectomycorrhizal and arbuscular mycorrhizal fungi, but it is currently unknown whether SSPs also play a role as effectors in other mycorrhizal associations. Ericoid mycorrhiza is a specific endomycorrhizal type that involves symbiotic fungi mostly belonging to the Leotiomycetes (Ascomycetes) and plants in the family Ericaceae. Genomic and RNASeq data from the ericoid mycorrhizal fungus Oidiodendron maius led to the identification of several symbiosis-upregulated genes encoding putative SSPs. OmSSP1, the most highly symbiosis up-regulated SSP, was found to share some features with fungal hydrophobins, even though it lacks the Pfam hydrophobin domain. Sequence alignment with other hydrophobins and hydrophobin-like fungal proteins placed OmSSP1 within Class I hydrophobins. However, the predicted features of OmSSP1 may suggest a distinct type of hydrophobin-like proteins. The presence of a predicted signal peptide and a yeast-based signal sequence trap assay demonstrate that OmSSP1 is secreted. OmSSP1 null-mutants showed a reduced capacity to form ericoid mycorrhiza with Vaccinium myrtillus roots, suggesting a role as effectors in the ericoid mycorrhizal interaction.
Fungi secrete a wide range of enzymatic and non-enzymatic proteins that function in the break-down of complex organic molecules but also in the interaction with microbial competitors or with animal and plant partners (Stergiopoulos and de Wit, 2009; Tian et al., 2009; Scherlach et al., 2013; Talbot et al., 2013). Fungi can establish different types of interactions with plants, ranging from mutualistic to antagonistic. Whatever their lifestyle, plant-colonizing fungi are recognized by the plant immune system through invariant molecular patterns known as microbe- or pathogen-associated molecular patterns (MAMPS or PAMPs) (Jones and Dangl, 2006). To successfully colonize plant tissues, fungi must prevent the PAMP-triggered immunity (PTI) reaction (Lo Presti et al., 2015). For this purpose, fungi secrete effector molecules that may play different functions depending on the fungal lifestyles. For example, they can be toxic compounds that kill the host plant (in necrotrophs), or secreted proteins that shield the fungus and suppress the host immune response, or proteins that manipulate the host cell metabolism to allow plant tissue invasion and nutrient uptake (de Jonge et al., 2011; Giraldo et al., 2013; Selin et al., 2016). Many small secreted proteins (SSPs) have been reported to function as effectors (Lo Presti et al., 2015).
Effectors were initially considered as virulence factors secreted exclusively by pathogens (van Esse et al., 2008; Stergiopoulos and de Wit, 2009). However, it has become apparent that effectors can manipulate the plant immune system also in mutualistic associations (Kim et al., 2016). Mutualistic fungi establish intimate contacts with plants by forming specialized fungal structures involved in nutrient exchange with the host (Martin et al., 2017). Effector-like SSPs have been functionally characterized as effector-like molecules both in arbuscular (AM) and ectomycorrhizal (ECM) fungi, as well as in some endophytic fungi (Plett and Martin, 2015). For example, the ECM fungus Laccaria bicolor requires MiSSP7 (Mycorrhiza-induced Small Secreted Protein 7) to establish symbiosis. MiSSP7 suppresses the plant defense reactions by interacting with the jasmonate co-receptor JAZ6 (Plett et al., 2011, 2014). Similarly, the AM fungus Rhizophagus irregularis secretes SP7, an effector protein that counteracts the plant immune program by interacting with the pathogenesis-related transcription factor ERF19, leading to increased mycorrhization (Kloppholz et al., 2011). Tsuzuki et al. (2016) also showed that host-induced gene silencing of SlS1, a putative secreted R. irregularis SSP expressed in symbiosis, resulted in suppression of colonization and formation of stunted arbuscules. In a similar manner, the candidate effector protein PIIN_08944, secreted by the fungal endophyte Piriformospora indica during colonization of both Arabidopsis and barley plants, was demonstrated to play a crucial role in reducing the expression of PTI genes and of the salicylic acid defense pathway (Akum et al., 2015).
Protein effectors are normally secreted following the endoplasmic reticulum-Golgi apparatus pathway, and bioinformatic identification of effector candidates can thus be based on the presence of the N-terminal signal peptide (Lo Presti et al., 2015), even though alternative secretion pathways have been reported for Magnaporthe oryzae and Phytopthora infestans (Giraldo et al., 2013; Wang et al., 2017). General SSPs features are: the presence of a signal peptide and the absence of transmembrane domains or GPI-anchor sites; a small size, with a mature length smaller than 300 amino acids; a richness in cysteine residues and, sometimes, the presence of conserved motifs (Martin et al., 2008; Stergiopoulos and de Wit, 2009; Hacquard et al., 2012; Zuccaro et al., 2014; Lo Presti et al., 2015).
Bioinformatic analyses of about 50 fungal genomes has highlighted that, when compared with saprotrophic and pathogenic fungi, the ECM fungal secretome is enriched in SSPs and contains species-specific SSPs likely dedicated to the molecular cross-talk between fungal and plant partners (Pellegrin et al., 2015). In line with this finding, a comparative in silico analysis of the AM fungi Rhizophagus clarus, R. irregularis and Gigaspora rosea highlighted the presence of shared SSPs (Sedzielewska Toro and Brachmann, 2016; Kamel et al., 2017), supporting a general conserved role of SSPs in AM. These data suggest that effector SSPs may represent an important fungal “toolkit” that enables the establishment/maintenance of host plant colonization in mycorrhiza (Plett and Martin, 2015; Martin et al., 2016). However, there is increasing awareness that SSPs in mycorrhizal fungi are likely involved in additional functions unrelated to symbiosis. For example, several SSPs are secreted by the ECM fungi L. bicolor and Hebeloma cylindrosporum during the free-living phase (Vincent et al., 2012; Doré et al., 2015). Moreover, large-scale transcriptomic and genomic analyses including fungi with different lifestyles revealed a wide array of SSPs in most saprotrophic fungi (Pellegrin et al., 2015; Valette et al., 2017), suggesting a possible role for SSPs in competition and rhizospheric communication (Rovenich et al., 2014).
Whereas effector SSPs have been characterized in AM and ECM fungi, there is currently no information on the occurrence of SSPs in ericoid mycorrhizal (ERM) fungi and on their potential role in symbiosis. ERM fungi are soil-borne fungi mostly belonging to Leotiomycetes (Ascomycetes). They form a peculiar endomycorrhizal type by colonizing the root epidermal cells of plants within the family Ericaceae and promote growth of their host plant in stressful habitats (Perotto et al., 2012). A common ERM fungal species is Oidiodendron maius (Dalpé, 1986) and O. maius strain Zn, an isolate from a metal polluted soil whose genome and transcriptome have been recently sequenced (Kohler et al., 2015), has become a model system to investigate metal stress tolerance in these fungi (Perotto et al., 2012; Daghino et al., 2016; Ruytinx et al., 2016). Genomic data have revealed a large set of carbohydrate-active enzymes (CAZymes) in O. maius, with many plant cell wall degrading enzymes being expressed during symbiosis (Kohler et al., 2015). As these enzymes could potentially elicit defense reactions through oligosaccharide release, symbiosis development likely requires a tight control of the plant defense reactions and, based on our current knowledge of arbuscular and ectomycorrhizal interactions, effectors to control plant immunity. Aim of this work was to identify, through the analysis of O. maius genomic and transcriptomic data, fungal SSPs potentially involved in the molecular dialog governing the ERM symbiosis.
Materials and Methods
Fungal Strains and Growth Conditions
Oidiodendron maius strain Zn (hereafter O. maius) was isolated from the roots of V. myrtillus growing in the Niepolomice Forest (Poland), and first described by Martino et al. (2000). This O. maius strain is deposited at the Mycotheca Universitatis Taurinensis collection (MUT1381; University of Turin, Italy) and at the American Type Culture Collection (ATCC MYA-4765; Manassas, VA, US), and was maintained on Czapek-Dox solid medium (NaNO3 2 g L−1, KCl 0.5 g L−1, glycerol phosphate*H2O 0.5 g L−1, K2HPO4 0.35 g L−1, FeSO4 0.01 g L−1, sucrose 30 g L−1, agar 10 gL−1, adjusted to pH 6). O. maius and OmΔSSP1-null mutants were also grown in the presence of different stressor compounds. Czapek-Dox medium was supplemented with 0.3 mM Cd (as 38H2O), 15 mM Zn (as 7H2O), 117.6 mM H2O2, 0.75 mM menadione, 0.1% (w:v) caffeic acid, 0.5% (w:v) tannic acid, 0.5% (w:v) gallic acid and 0.5% (w:v) quercetin. Prior to fungal inoculation, sterile cellophane membranes were placed on the agar surface to provide a convenient means of removing the mycelium from the plate. The membranes were first boiled for 15 min in 10 mM EDTA (disodium salt, dihydrate, SIGMA), rinsed and then autoclaved in ddH2O. Fungal colonies were removed after 30 days, dried over-night and weighted.
In Vitro Mycorrhizal Synthesis
Axenic V. myrtillus seedlings were obtained from seeds (Les Semences du Puy, Le Puy-En-Velay, France) surface sterilized in 70% ethanol (v:v) 0.2% Tween 20 for 3 min, rinsed with sterile water, submerged in 0.25% sodium hypochlorite for 15 min and rinsed again with sterile water. Seeds were germinated on 1% water agar for 2 weeks in darkness before transfer to a growth chamber for 1 month.
Mycorrhiza was synthesized in petri plates containing Modified Melin-Norkrans (MMN) medium with the following composition: KH2PO4 0.5 gL−1, Bovine Serum Albumin (BSA) 0.1 gL−1, 2H2O 0.066 gL−1, NaCl 0.025 gL−1, 7H2O 0.15 gL−1, thiamine-HCl 0.1 gL−1, 6H2O 0.001 gL−1, agar 10 gL−1, final pH 4.7. Sterile cellophane membranes, prepared as described before, were placed on the agar surface before fungal inoculation. A suspension of O. maius conidia in sterile deionised water was distributed on the cellophane membranes in the bottom half of the MMN petri plates. Ten germinated V. myrtillus seedlings were then transferred just above the conidia suspension. Plates were sealed and placed in a growth chamber (16-h photoperiod, light at 170 μmol m−2 s−1, temperatures at 23°C day and 21°C night). Roots were collected and the percentage of mycorrhization evaluated after 45 day.
As a control for the asymbiotic condition, O. maius was grown on the same medium used for mycorrhizal synthesis. Plates covered by cellophane membranes were inoculated with 5 mm fungal plugs and fungal colonies were removed after 45 days. Three biological replicates were prepared for each sample for the RNA-Seq experiment.
RNA Extraction and RNA-Seq Data Analyses
Total RNA was extracted and quantified from 100 mg aliquots of O. maius mycelium and O. maius-inoculated V. myrtillus collected 45 days after inoculation, frozen in liquid nitrogen and mechanically ground. Total RNA was extracted from O. maius mycelium using a Tris-HCl extraction buffer and from V. myrtillus mycorrhizal roots using the CTAB method, as described by Kohler et al. (2015).
Preparation of libraries from total RNA and 2 × 100 bp Illumina HiSeq sequencing (RNA-Seq) was performed by IGA Technology Services (Udine, Italy). Raw reads were trimmed and aligned to the respective reference transcripts available at the JGI MycoCosm database (http://genome.jgi-psf.org/programs/fungi/index.jsf) using CLC Genomics Workbench v6. For mapping, the minimum length fraction was 0.9, the minimum similarity fraction 0.8 and the maximum number of hits for a read was set to 10. The unique and total mapped reads number for each transcript were determined, and then normalized to RPKM (Reads Per Kilobase of exon model per Million mapped reads). A summary of the aligned reads is given in Table S1. The complete data set has been deposited in NCBI Gene Expression Omnibus and is accessible through GEO Series accession numbers GSE63947. To identify differentially regulated transcripts in mycorrhizal tissues compared to free-living mycelium, the Baggerly's Test (Baggerly et al., 2003) implemented in CLC Genomic workbench was used. This test compares the proportions of counts in a group of samples against those of another group of samples. The samples are given different weights depending on their sizes (total counts). The weights are obtained by assuming a Beta distribution on the proportions in a group, and estimating these, along with the proportion of a binomial distribution, by the method of moments. The result is a weighted t-type test statistic. In addition, Benjamini & Hochberg multiple-hypothesis testing corrections with False Discovery Rate (FDR) were used. Transcripts with a more then 5-fold change and a FDR corrected p < 0.05 were kept for further analysis.
O. maius Small Secreted Proteins (SSPs) were identified using a custom pipeline including SignalP v4 (1), WolfPSort (2), TMHMM, TargetP (3), and PS-Scan algorithms (4) as reported in Pellegrin et al. (2015). To assess whether symbiosis-regulated transcripts were conserved or lineage-specific (i.e., orphan genes with no similarity to known sequences in DNA databases), their protein sequences were queried against the protein repertoires of 59 fungal genomes using BLASTP with e-value 1e-5. Proteins were considered as orthologs of symbiosis-regulated transcripts pending they showed 70% coverage over the regulated sequence and at least 30% amino acid identity.
cDNA Synthesis and Quantitative RT-PCR (RT-qPCR)
The expression of seven selected SSPs was evaluated by RT-qPCR. cDNA was obtained from about 1,000 ng of total RNA with a reaction mix containing 10 μM random primers, 0.5 mM dNTPs, 4 μl 5 × buffer, 2 μl 0.1 M DTT, and 1 μl Superscript II Reverse Transcriptase (Invitrogen) in a final volume of 20 μl. Temperature regime was: 65°C for 5 min, 25°C for 10 min, 42°C for 50 min, and 70°C for 15 min. Possible DNA contamination was tested with an additional PCR reaction using a specific primer pair that amplifies an intron containing region of the O. maius Elongation Factor1α (OmEF1α) (Table S2). RT-qPCR was performed with the Rotor-Gene Q (Qiagen) apparatus. The reactions were carried out in a final volume of 15 μl with 7.5 μl of iQ SYBR Green Supermix (Bio-rad), 5.5 μl of forward and reverse primers (10 μM stock concentration; Table S2) and 2 μl of cDNA (diluited 1:10). qPCR cycling program consisted of a 10 min/95°C holding step followed by 40 cycles of two steps (15 s/95°C and 1 min/60°C). The relative expression of the target transcript was measured using the 2−ΔCt method (Livak and Schmittgen, 2001). The Omβ-Tubulin (OmβTub) (Table S2) was used as reference house-keeping gene. Three to five biological replicates and two technical replicates were analyzed for each condition tested. qPCR primers were designed with Primer3Plus (http://www.bioinformatics.nl/cgi-bin/primer3plus/primer3plus.cgi) and checked for specificity and secondary structure formation with PrimerBlast (http://www.ncbi.nlm.nih.gov/tools/primer-blast/) and OligoAnalyzer (eu.idtdna.com/calc/analyzer). Primers were synthesized by Eurogentec (Belgium).
Construction of the OmSSP1-Disruption Vector and Agrobacterium-Mediated Transformation
OmSSP1-null mutants were obtained through Agrobacterium tumefaciens-mediated (ATM) homologous recombination. PCR reactions were used to produce the 5′ upstream flanking region (1502 bp) and the 3′ downstream flanking region (1,533 bp) of the OmSSP1 gene. PCR reactions were carried out in a final volume of 50 μl containing: 50 ng of genomic DNA of O. maius Zn, 1 μl dNTPs 10 mM, 2.5 μl of each primer (10 μM stock concentration; Table S2), 10 μl of 5 × Phusion HF Buffer and 0.5 units of Phusion Hot Start II High-Fidelity (Thermo Scientific). The PCR program was as follows: 30 s at 98°C for 1 cycle, 10 s at 98°C, 30 s at 60°C, 45 s at 72°C for 30 cycles, 10 min at 72°C for 1 cycle. Amplicons were then purified with Wizard® SV gel and PCR clean-up system (PROMEGA) following the manufacturer's instructions. PCRs amplicons were cut with XmaI-HindIII (for the 5′) and BglII-HpaI (for the 3′) and cloned into the pCAMBIA0380_HYG vector (Fiorilli et al., 2016) in order to obtain the pCAMBIA0380_HYG-ΔOmSSP1 vector (Figure S1). The restriction reactions were performed in 30 μl final volume containing 0.5 μg of DNA (1 μg for the plasmid), 0.5 μl of each enzyme (from PROMEGA), 0.3 μl of BSA 100X and 3 μl of buffer 10X, overnight at 37°C. The ligase reaction was carried out in 20 μl final volume containing 50 ng of vector, 18 ng of the amplicon, 2 μl of buffer 10X and 1 μl of T4 enzyme (PROMEGA), overnight at 4°C. The vector sequence was checked by PCR and DNA sequencing.
The vector was cloned into Agrobacterium tumefaciens LBA1100, that was used to transform ungerminated O. maius conidia according to the protocol described in Abbà et al. (2009).
Identification of OmSSP1-Null Mutants by PCR and Southern Blot
Fungal transformants were screened by PCR. A small portion of each fungal colony was collected and boiled for 15 min in 20 μl of 10 mM Tris HCl pH 8.2, vortexed for 1 min and centrifuged 15 min at room temperature. Then, 2 μl of the supernatant were used directly for PCR amplification without any further purification, using two sets of primers (Table S2). The first primer set was designed to amplify the OmSSP1 gene (OmSSPb1r and OmSSPb1f) whereas the second set (Hyg4f e Hyg2r) was designed to amplify the portion of the hph gene corresponding to the Hyg-probe (Figure S1). A OmSSP1-null mutant would yield an amplified product only with the second primer set. The putative OmSSP1-null mutants identified were validated by PCR using primers (Table S2) designed on the genome at the 5′ and 3′ of the homologous recombination site (respectively preOmSSPb1f3/Hyg6r and postOmSSPb1r3/Hyg3f). The positive OmSSP1-null mutants were further analyzed through Southern blot hybridization to verify single-copy integration of the disruption cassette in the genome. Fifteen micrograms of genomic DNA from the deletion mutants and from the wild-type were digested with BamHI (PROMEGA) and size-fractionated on 1% (w:v) agarose TAE 1X gel. The separated restriction fragments were blotted onto a nylon membrane following standard procedures (see Abbà et al., 2009). Hybridization with a probe designed on the Hygromycin-resistance cassette (Hyg-probe, Figure S1) was performed with a chemiluminescent detection system (ECL direct DNA labeling and detection system; GE Healthcare, U.K.) according to the manufacturer's recommendations.
Quantification of the Degree of Mycorrhization
To determine differences in root colonization between O. maius wild-type and OmSSP1-null mutants, the percentage of mycorrhization was recorded after 1.5 months. The roots of 3–6 seedlings colonized by each mutant strain were collected and the whole root system was stained overnight in a solution of lactic acid:glycerol:H2O (14:1:1) containing acid fuchsin 0.01% (w:v), destained twice with 80% lactic acid and observed using a Nikon Eclipse E400 optical microscope. The magnified intersections method (Villarreal-Ruiz et al., 2004) was adapted to quantify the percentage of fungal colonization of V. myrtillus hair roots, under the microscope, using the rectangle around the cross-hair as intersection area at 40 × magnification. A total of 60 intersections per seedling root system were scored. Counts were recorded as percentage of root cells colonized (RC) by the fungus using the formula: RC% = 100 × Σ of coils counted for all the intersections/Σ of epidermal cells counted for all the intersections.
Phylogenetic and Bioinformatic Analyses
The aminoacid sequences comprised between cysteine 1 and cysteine 8 of the fungal proteins listed in Table S3 were aligned using MUSCLE (MUltiple Sequence Comparison by Log-Expectation; Gap Open penalty−2, Edgar, 2004) tool implemented in MEGA7 (Tamura et al., 2007). Maximum likelihood analysis (Guindon and Gascuel, 2003) was conducted using www.phylogeny.fr in advanced mode (Dereeper et al., 2008). The phylogenetic tree was reconstructed using the maximum likelihood method implemented in the PhyML program (v3.1/3.0 aLRT).
Bioinformatic analyses of protein primary sequences were performed using online tools. Blastp searches on the Uniprot database (The UniProt Consortium, 2017) identified the closest protein matches. Hydropathy profiles were generated with ProtScale tools (http://web.expasy.org/protscale/), while the grand average of hydropathicity (GRAVY), the predicted amino acid number and molecular weight were calculated using the ProtParam tool (http://web.expasy.org/protparam/). The intrinsic solubility profiles were obtained with the camsolintrinsic calculator (http://www-mvsoftware.ch.cam.ac.uk/index.php/camsolintrinsic). The representation of the residues hydrophobicity of the aligned sequences of each clade was obtained using http://www.ibi.vu.nl (Simossis and Heringa, 2005).
The Yeast Signal Sequence Trap Assay
Functional validation of the predicted signal peptide of OmSSP1 was conducted with a yeast signal sequence trap assay (Plett et al., 2011). The pSUC2-GW gateway vector carries a truncated invertase (SUC2) lacking both its initiation methionine and signal peptide. cDNA encoding the predicted OmSSP1 signal peptide was cloned into pSUC2-GW plasmids using BP/LR Gateway technologies (Invitrogen). Then, yeast cells (YTK12 strain) were transformed with 200 ng of the pSUC2-GW/OmSSP1 plasmid using the lithium acetate method (Gietz and Schiestl, 2007). Transformants were grown on SD-W yeast minimal medium (6.7 gL−1 Yeast Nitrogen Base without amino acids, 0.7 gL−1 tryptophan dropout supplement, 20 gL−1 glucose, 20 gL−1 agar, pH 5.6) and on YPGA medium (10 gL−1 yeast extract, 20 gL−1 peptone, 20 gL−1 agar amended with 20 gL−1 glucose and 60 μg mL−1 antimycin A after autoclaving, pH 6.5). To assay for invertase secretion, colonies were grown overnight at 30°C with shaking (200 rpm) and diluted to an OD600 = 1, then 5 μl of serial dilution of the yeast culture were plated onto YPSA medium containing sucrose (10 gL−1 yeast extract, 20 gL−1 peptone, 20 gL−1 agar amended with 2 gL−1 sucrose and 60 μg mL−1 antimycin A after autoclaving, pH 6.5).
The significance of differences among the different treatments was statistically evaluated by ANOVA with Tukey's pairwise comparison as post-hoc test for multiple comparisons for normally distributed data. Statistical elaborations of growth and biomass data were performed using PAST statistical package, version 2.17 (Hammer et al., 2001). The differences were considered significant at a probability level of p < 0.05.
The O. maius Genome Contains Several SSPs Up-Regulated during Mycorrhizal Symbiosis with V. myrtillus
Among the 16,703 genes found in the O. maius genome (Kohler et al., 2015), 445 genes (~38% of the total O. maius predicted secretome) code for putatively secreted proteins smaller than 300 amino acids (Table S4a). The transcriptomic analysis of O. maius under free-living conditions (FLM) and in symbiosis with V. myrtillus (MYC) indicated that about 24% (278/1163) of the genes coding for putatively secreted proteins were induced in symbiosis (Fold Change > 5, p<0.05), 90 of them corresponding to SSPs. Many of these mycorrhiza-induced SSPs (MiSSP) were strongly up-regulated in symbiosis (13 with FC ≥ 400) or mycorrhiza-specific (Table S4b). Only 32 were cysteine enriched (C > 3%), a feature normally attributed to SSPs (Kim et al., 2016). About half (49/90) contained PFAM motifs specific of CAZymes (especially glycoside hydrolases, GHs), lipases, hydrophobins and peptidases, whereas the remaining 41 contained no known PFAM domain. None of these MiSSPs contained a nuclear localization signal motif but 6 of them featured a KR rich sequence, i.e., a motif characterized by basic aminoacids (like lysine and arginine) supporting entry in the plant nucleus (Table S4b).
Genes orthologous to the O. maius MiSSPs were identified by genomic comparative analyses with 59 taxonomically and ecologically distinct fungi (Table S4d), including three other ERM fungi in the Leotiomycetes (Meliniomyces bicolor, M. variabilis and Rhizoscyphus ericae) belonging to the “R. ericae” aggregate (Vrålstad et al., 2002). Ten of the 90 O. maius symbiosis-induced SSPs were specific for O. maius, whereas 2 to 834 orthologous genes were found for the other 80 SSPs (Table S4b). Many O. maius SSPs orthologs were found in the other three ERM fungal species, although no ERM specific SSPs could be identified. The highest number of O. maius SSPs orthologs was found in the genomes of pathogenic and saprotrophic fungi (64 and 73 respectively), whereas only 38 orthologous genes were found in the genomes of 12 ECM fungi (Table S4b).
We selected seven O. maius SSPs for further analyses, with a preference for those uniquely or highly expressed in symbiosis that did not contain PFAM motifs with known functions (Table 1). Three different software (PrediSi, SIgnalP 4.1, Phobius) confirmed the presence of a signal peptide (Table 1), but alignment of the primary sequences indicated very low similarities. Cysteine enrichment > 3% was only observed for OmSSP1 (8.6%), the most highly symbiosis-induced SSP (Table 1). The expression of these OmSSPs was investigated in free-living and in mycorrhizal O. maius by RT-qPCR, which confirmed a significant up-regulation of all selected SSP genes in symbiosis, and in particular a very strong induction of OmSSP1 (Figure 1).
Table 1. Expression level, biochemical and bioinformatic features of seven selected O. maius mycorrhiza-induced SSPs.
Figure 1. RT-qPCR validation of selected O. maius SSPs expression level in mycorrhizal tissues of V. myrtillus (MYC) as compared to the free living mycelium (FLM). Relative expression (mean of ΔCt values) of 7 O. maius SSPs was normalized to OmTub transcripts. Bars represent the mean ±SD. Asterisks (*) indicates a statistically significant difference (p < 0.05) (ANOVA, Tukey's post hoc test).
OmSSP1 Shows Molecular Similarities With Fungal Proteins Annotated as Hydrophobins
OmSSP1 is a single copy gene located in the scaffold 14 of the O. maius genome. The coding region contains 354 nucleotides, with 2 exons and 1 intron. In addition to the signal peptide, OmSSP1 is rich in glycine (16.1%) and leucine (12.9%) and contains 8 cysteine residues (Figure 2). The intrinsic calculated solubility and solvent accessibility highlighted at least three poorly water soluble regions. The majority of hydrophobic residues in OmSSP1 are clustered between amino acid residues 41–56 and 79–93 (Figure 2). OmSSP1 secondary structure predictions indicated a disordered folded state for 46% of the protein, partly due to the presence of a low complexity region (LCR) before the first cysteine (Figure 2), whereas 16 and 20% of the protein structure could form helix and beta-sheet structures, respectively (not shown).
Figure 2. OmSSP1 structural and biochemical features predicted through bioinformatics tools. (A) Schematic representation of OmSSP1 primary amino acid sequence. The green bar represents the signal peptide (SP), the blue bar the low complexity region (LCR) and the red bars the cysteine (C) residues. (B) Calculation of the intrinsic solubility profile (ISP) and (C) hydrophaty profile (HP). For the ISP, scores larger than 1 indicate highly soluble regions, while scores smaller than −1 indicate poorly soluble regions. For the HP, hydrophobic regions show positive peaks with values above 0 whereas hydrophilic regions show negative peaks. The position of the C3–C4 loop, characterized by a conserved hydrophobic core, is indicated. The sequences of A–C have been graphically aligned in order to show the correspondence between the elements of the primary sequence (A), their solubility (B) and hydrophobicity profile (C).
Although no known domains could be found in the predicted OmSSP1 protein sequence by PFAM database searching, nor a functional classification in the InterPro database, BlastP searches both in the UniProt and in the RefSeq databases yielded, as closest protein matches, fungal hydrophobins or hydrophobin-like proteins (Table 2). Indeed, the two OmSSP1 orthologous genes in T. terrestris (protein ID 2089872 and 209296) showed high homology with fungal hydrophobins. The hydropathy profile of OmSSP1 was analyzed and a positive GRAVY value was found (0.53), indicating overall hydrophobicity (Table 1). Therefore, we compared the protein sequence of OmSSP1 with those of the four O. maius annotated hydrophobins and of other fungal hydrophobins (Linder et al., 2005; Seidl-Seiboth et al., 2011; Grigoriev et al., 2014). Three out of the four O. maius hydrophobins showed a GRAVY score above 0.6, indicating an overall hydrophobicity higher than OmSSP1 (0.53), whereas the GRAVY score for O. maius hydrophobin 4 was only 0.32 (Table S5).
Two classes of fungal hydrophobins have been distinguished by Wessels (1994) on the basis of the hydropathy profile and of slightly different motifs between the eight positionally conserved cysteine residues: Class I (C-X5−8-CC-X17−39-C-X8−23-C-X5−6-CC-X6−18-C-X2−13), and Class II (C-X9−10-CC-X11-C-X16-C-X8−9-CC-X10-C-X6−7). OmSSP1 bears a signature motif similar to class I hydrophobins (C-X7-CC-X7-C-X8-C-X5-CC-X16-C-X2), with the exception of a shorter sequence (only 7 aminoacids) between cysteines 3 and 4, expected to host the most hydrophobic protein region. Although shorter, the C3-C4 loop of OmSSP1 showed a hydrophobic core composed by valine (V) and leucine (L), thus explaining the negative values in solubility profile and positive values in hydropathy profile (Figure 2). Hydrophobins are amphiphilic molecules (Whiteford and Spanu, 2002; Rineau et al., 2017), and OmSSP1 has a very hydrophilic stretch rich in G before the C3-C4 loop, corresponding to the LCR (Figure 2).
To better understand the phylogenetic relatedness of OmSSP1 with fungal hydrophobins, we aligned the OmSSP1 protein sequence with the four annotated O. maius hydrophobins and with Class I and Class II annotated hydrophobins from Ascomycetes (listed in Table S3). A phylogenetic tree built on the complete C1-C8 sequence alignment is shown in Figure 3. However, since the different proteins showed a highly variable sequence length of the C3-C4 loop, ranging from 4 to 39 amino acids, a more conserved phylogenetic tree was generated by Maximum Likelihood (ML) without this protein region, to avoid possible bias due to the different protein lengths (Figure S2). In both trees, class II hydrophobins and two of the four O. maius hydrophobins (Oidma2 and Oidma3) grouped in a single, well supported cluster. The two other O. maius hydrophobins grouped in a well-supported cluster together with characterized Class I hydrophobins (Figure 3 and Figure S2). Although the position of some fungal proteins differed in the two ML trees, most terminal clades (A–D) were maintained and well supported. OmSSP1 clustered in Clade C (Figure 3 and Figure S2) together with other proteins reported as hydrophobins (Table S3) and featuring a short C3-C4 loop (X4−9). The complete C1-C8 sequence alignment of proteins in Clade C is shown in Figure 4. As in OmSSP1, most amino acids in the short C3-C4 loop of other proteins in this clade were hydrophobic. Clade B (Figure 3 and Figure S2) was another well-supported clade containing proteins also featuring a very short (X8−9) C3-C4 loop (Figure 4). Clade B included a Trichoderma atroviride hydrophobin (Triat1) described by Seidl-Seiboth et al. (2011) as a member of a novel subclass in Class I hydrophobins. However, unlike Clade C, proteins in Clade B featured mainly hydrophilic amino acids in the C3-C4 loop (Figure 4). Clade A included the two O. maius hydrophobins Oidma1 and Oidma4.
Figure 3. Phylogenetic tree of OmSSP1 and O. maius hydrophobins with other annotated hydrophobins from Ascomycetes. The analysis included protein sequences annotated or described as Ascomycetes class I and class II hydrophobins (listed in Table S3). This sequence alignment considered the complete amino acid sequence comprised between C1 and C8. Muscle algorithm implemented in MEGA7 (Tamura et al., 2007) was used to generate the multiple protein sequence alignment. The phylogenetic tree was reconstructed on the Phylogeny.fr platform (Dereeper et al., 2008) using the maximum likelihood method (Guindon and Gascuel, 2003) implemented in the PhyML program (v3.1/3.0 aLRT). The WAG substitution model was selected assuming an estimated proportion of invariant sites (of 0.088) and 4 gamma-distributed rate categories to account for rate heterogeneity across sites. The gamma shape parameter was estimated directly from the data (gamma = 7.097). Reliability for internal branch was assessed using the aLRT test (SH-Like). Graphical representation and edition of the phylogenetic tree were performed with TreeDyn (v198.3; Chevenet et al., 2006) A–D, terminal clades.
Figure 4. Aminoacid hydrophobicity properties of the aligned Class I hydrophobins. The OmSSP1 and hydrophobin protein sequences belonging to the four clades (A–D) in Figure 3 were aligned with the Praline tool of the IBIVU server (http://www.ibi.vu.nl/programs/pralinewww/; Simossis et al., 2005). The hydrophobicity scale used is from Eisenberg et al. (1984). The C3-C4 loop is indicated above each clade alignment by a string of asterisks.
The Yeast Invertase Secretion Assay Indicates That OmSSP1 Is Secreted
The predicted OmSSP1 signal peptide was functionally validated in a yeast signal sequence trap assay. This test is based on the yeast requirement for a secreted invertase (SUC) to grow on sucrose amended media (Klein et al., 1996). The pSUC-GW (Jacobs et al., 1997) gateway vector carries a truncated invertase that lacks its signal peptide (SUC2). The cDNA coding for the putative OmSSP1 signal peptide (OmSSP1_SP) was fused in frame to the yeast SUC2 invertase, and the recombinant pSUC-GW vector was transformed into the invertase-deficient yeast strain YTK12. As positive control, the MiSSP7 (Plett et al., 2011) signal peptide (MiSSP7_SP) and the yeast wild-type signal peptide (SUC2_SP+) were used, whereas the empty vector was used as negative (mock) control. All transformants grew on SD-W and YPGA control media containing glucose, whereas only OmSSP1_SP and the two positive controls rescued YTK12 yeast's growth on YPSA media containing sucrose, indicating that OmSSP1 contains a secretion signal that is functional in yeast (Figure 5).
Figure 5. OmSSP1 contains a signal sequence that is functional in yeast. (A) Upon sucrose selection (YPSA medium) the OmSSP1 signal peptide (OmSSP1_SP) rescues the functionality of the yeast invertase, like the other two positive controls (The MiSSP7 signal peptide - MiSSP7_SP - and the wild-type sequence of the yeast signal peptide - SUC2_SP +). The empty vector (mock) was used as a negative control. Control media (B) SD-W and (C) YPGA, both containing glucose, restore the growth ability of all yeast transformants. Five microliter of undiluted (ud) or serially diluted (10−1 and 10−2) yeast culture were plated.
Growth of OmΔSSP1 Mutants Is Not Impaired Under Stressful Conditions, but They Have a Reduced Capability to Colonize V. myrtillus Roots
O. maius can be genetically transformed and gene disruption can be obtained by homologous recombination (Martino et al., 2007; Abbà et al., 2009). To investigate the biological function of OmSSP1 in O. maius, knock-out mutants (OmΔSSP1) were obtained through AMT transformation using a vector containing the hpd-cassette (Figure S1). Hygromycin-resistant colonies were screened by PCR, and eight putative homologous recombinants out of 742 screened fungal transformants could be identified. Southern blot hybridization with a probe to the Hygromycin cassette showed a single band of the expected size for four out of the eight candidate mutants (Figure S3). To further confirm the vector insertion site, PCR amplifications were performed with primers designed to amplify the genome regions flanking the inserted pCAMBIA0380ΔOmSSP1 disruption cassette (Figure S1, Table S2), followed by sequencing of the amplicons. Three O. maius transformants (OmΔSSP1150, OmΔSSP1377, OmΔSSP1412) were confirmed as OmSSP1 deletion mutants. These three OmSSP1-null mutants were not affected in mycelium morphology or growth rate when inoculated on Czapek-Dox solid medium (not shown). OmSSP1 deletion did not modify the wettability phenotype of the O. maius mycelium (Figure S4), but it should be noted that OmSSP1 gene expression was very low in the FLM (Table S4b).
It has been recently suggested that, in FLM, SSPs may increase fungal tolerance to toxic compounds (such as aromatic compounds or reactive oxygen species) released during substrate degradation (Valette et al., 2017). We therefore investigated whether deletion of the OmSSP1 gene reduced O. maius fitness when exposed to different stress inducers. As shown in Figure S5, growth of the three OmΔSSP1 mutants was not significantly affected by any of the stress conditions tested, namely two toxic heavy metals (Cd and Zn), molecules causing oxidative stress (H2O2 and menadione) and plant-derived organic compounds displaying toxic/antimicrobical effects (caffeine, tannic acid, gallic acid, quercetin and caffeic acid).
We then investigated the symbiotic ability of the three OmΔSSP1 mutants on seedlings of the host plant V. myrtillus as compared with the wild-type O. maius strain (Figure 6). Although there were no statistically significant differences in plant biomass after 45 days of co-culture, a statistically significant (p < 0.05) reduction in the percentage of root colonization was measured for all OmΔSSP1 mutants, when compared with the wild-type O. maius strain (Figure 6).
Figure 6. OmΔSSP1 mutants have a reduced capability to colonize V. myrtillus roots. V. myrtillus roots were observed after 1.5 months of co-culture with O. maius WT and with three OmSSP1 null-mutant strains. (A) The percentage of the root colonization was significantly lower for the OmSSP1 null-mutants as compared to the O. maius WT. (B) Quantification of fresh plant biomass (roots - gray bars - and aboveground portions - white bars) of V. myrtillus plants grown alone, in the presence of the O. maius WT strain or of the OmSSP1 null-mutants. Bars represent the mean ± SD, n = 5. Each biological replicate represents the total biomass of eight V. myrtillus seedlings grown in an individual plate. Different letters indicate statistically significant difference (p < 0.05, ANOVA, Tukey's post hoc test).
Similar to Other Mycorrhizal Fungi, the ERM Fungus O. maius Has a Wide Array of SSPs
Effector-like SSPs have been found to be secreted by ECM and AM fungi and to be instrumental for plant colonization (Plett and Martin, 2015; Martin et al., 2016). SSPs are also encoded in the genome of the model ERM fungus O. maius. Overall, the O. maius genome contains 445 SSPs, corresponding to 2.6% of the total number of O. maius genes. Similar percentages were reported for the ECM fungus L. bicolor (Pellegrin et al., 2015) and for saprotrophic fungi (Valette et al., 2017). The 90 mycorrhiza-induced SSPs (MiSSPs) correspond to about 20% of the total O. maius SSPs, a percentage similar to the ECM fungus L. bicolor (Kohler et al., 2015) and the AM fungus R. irregularis (Tisserant et al., 2013).
Bioinformatic analysis of O. maius MiSSPs revealed that 45.5% (41/90) correspond to orphans genes with no known PFAM domains. Many fungal effector SSPs are targeted to the host plant nucleus (Lo Presti et al., 2015). Although none of the O. maius MiSSPs contained a nuclear localization signal motif, six O. maius MiSSPs featured a KR reach sequence supporting a localization in the plant nucleus.
Comparative genomics revealed that 27% of the total O. maius SSPs are species-specific (Table S4a). Species-specific SSPs (SSSPs), defined as SSPs with no homology in other species, have been found in both AM (Salvioli et al., 2016; Sedzielewska Toro and Brachmann, 2016; Tang et al., 2016) and ECM fungi (Pellegrin et al., 2015) and they are considered to be likely involved in the promotion of host-specific interactions (Pellegrin et al., 2015). The number of SSSPs found in O. maius falls in the proportion predicted by Kim et al. (2016) for symbiotic organisms (25-50%), although only 10 out of the 90 O. maius SSPs up-regulated in symbiosis were species-specific (Table S4b).
Whereas lifestyle-specific SSPs have been found in ECM fungi (Pellegrin et al., 2015), no SSPs shared by all four ERM fungi (with no orthologs in the other fungi used for comparative analysis) could be found. ERM and ECM fungi also differed because comparative analyses showed that ECM fungi share most of their SSPs with saprotrophs such as brown rot, white rot, and litter decayers (Pellegrin et al., 2015), whereas the highest number of O. maius symbiosis-induced SSPs orthologous belong to pathogenic and saprotrophic fungi.
SSPs have been recently found in most fungal species regardless of their lifestyle, suggesting that they could be involved in a variety of processes, both common and lifestyle-specific (Pellegrin et al., 2015). Interestingly, three orthologous genes of OmSSP1, the most highly up-regulated (ca. 20,000-folds) O. maius MiSSP, were found not in other ERM fungi but in the two saprotrophic fungi Neurospora crassa and Thelavia terrestris. Switching among endophytic, pathogenic, and saprotrophic lifestyles, likely controlled by both environmental and host factors, were suggested for N. crassa (Kuo et al., 2014). T. terrestris, a thermophilic saprotrophic fungus phylogenetically related to N. crassa, shows important cellulases and hemicellulases activities (Berka et al., 2015). Thus, N. crassa and T. terrestris share important characteristics with O. maius, which besides being an endomycorrhizal fungus, is also commonly isolated from roots of other plants as well as from substrates rich in plant-derived organic matter (Rice and Currah, 2006). Moreover, the O. maius gene content for cell wall degrading enzymes, proteases and lipases places this fungus closer to saprotrophs and pathogens than to other mycorrhizal fungi (Martino et al., 2018), and may explain the fact that OmSSP1 orthologs are also found in SAP S/L/O (soil/litter/organic matter saprotrophs) fungi.
OmSSP1, the Most Highly Expressed O. maius MiSSP, May be a Distinctive Type of Hydrophobin
Hydrophobins are small secreted proteins less than 200 amino acids long, with a secretion signal and a pattern of eight cysteine residues recurring in the sequence (Whiteford and Spanu, 2002). A unique three-dimensional folding comes from these features, keeping exposed the hydrophobic residues and rendering them amphiphilic (Rineau et al., 2017). Two classes of hydrophobins have been recognized (Wessels, 1994): Class I, with higher sequence variability and more stable superstructures, is found in both Asco- and Basidiomycetes, while Class II has only been found in Ascomycetes (Kershaw and Talbot, 1998). Hydrophobins are abundantly expressed during fungal development, pathogenesis and symbiosis (Wösten, 2001; Whiteford and Spanu, 2002). Being amphiphilic, they could behave as biosurfactants and facilitate fungal adhesion to organic matter and its decomposition (Rineau et al., 2017). Hydrophobins are also instrumental for fungal hyphae to form aerial structures and to adhere to each other and/or to hydrophobic surfaces, such as the plant leaf surface during pathogenesis. Symbiosis-upregulated hydrophobins have been found in the ECM fungi Pisolithus tinctorius (Tagu et al., 2001) and in L. bicolor (Martin et al., 2008; Plett et al., 2012), where they could play a role in establishing hyphal aggregation in the symbiotic interfaces (Raudaskoski and Kothe, 2015).
The O. maius genome features four annotated hydrophobins containing the PFAM and InterPro hydrophobin domains. Rineau et al. (2017) suggested that all O. maius hydrophobins belong to Class I, but our phylogenetic analysis (that also included Class II hydrophobins) showed that two proteins (Oidma1 and Oidma4) belong to Class I and two (Oidma2 and Oidma3) to Class II hydrophobins. OmSSP1 shares some features with Class I hydrophobins and clusters with annotated hydrophobins in this Class, but it was not identified as a hydrophobin because it lacks the corresponding PFAM and InterPro domains, possibly because of the shorter C3–C4 region.
Amino acid features, such as charge and hydrophobicity, can influence hydrophobin structure and function. Thus, the amino acidic composition of the C3-C4 loop as well as of the N-terminal region of hydrophobins may influence the wettability and the substrate-attachment preference of the protein (Linder et al., 2005; Kwan et al., 2006). In this respect, it is interesting to note that proteins in Clade B and Clade C, both showing a C3-C4 loop unusually short for Class I hydrophobins, feature amino acid sequences with very different hydrophobicity (Figure 4), suggesting they may represent structurally and functionally diverse subclasses of Class I hydrophobins. The low complexity region found in OmSSP1 is also unusual for hydrophobins and could be considered a recently evolved trait of this protein (Toll-Riera et al., 2012). Its presence suggests for OmSSP1 a low propension to aggregate and to form alpha-helices and beta-sheets, three properties often correlated with the ability of hydrophobins to pile up in needle-like (amyloid) structures (Rineau et al., 2017).
A high number of hydrophilic residues (asparagine especially) were found in the N-amino terminal region of OmSSP1. According to Linder et al. (2005), the amino terminal region of hydrophobins could have important roles in the specific function of individual proteins. For example, hydrophobins featuring high number of exposed hydrophilic residues at the N-terminal region were found to be overexpressed in mycorrhizal tissues (Whiteford and Spanu, 2002; Rineau et al., 2017).
OmSSP1 Null-Mutants Have a Reduced Ability to Colonize V. myrtillus Roots
There is increasing awareness that SSPs may play important roles during saprotrophic fungal growth, as they have been identified in saprotrophic fungi and they can be expressed by mycorrhizal fungi during asymbiotic growth (Vincent et al., 2012; Doré et al., 2015; Valette et al., 2017). However, although a limited range of growth conditions were tested, OmSSP1 did not appear to be necessary in the FLM, as the three OmSSP1-null mutants were not affected in mycelium morphology or growth rate, even when they were exposed to toxic and oxidative chemical compounds. By contrast, when they were tested for symbiotic capabilities on V. myrtillus plants, a significant reduction in the percentage of mycorrhization (from about 37% to about 23-24%) was measured as compared with the wild-type strain, thus suggesting a specific role of OmSSP1 in the mycorrhization process. The OmSSP1 deletion did not fully prevent root mycorrhization. However, 20% of O. maius SSPs are induced in symbiosis, and although OmSSP1 was the most highly up-regulated, we cannot exclude a functional redundancy, as already reported for the effectors of pathogenic fungi (Selin et al., 2016). Thus, the absence of OmSSP1 could be partly compensated by other O. maius SSPs with similar function, thus lowering the impact of the OmSSP1 deletion.
In conclusion, the genome of the ERM fungus O. maius contains several SSPs that are up-regulated in symbiosis. Decreased colonization of V. myrtillus roots by OmSSP1-null mutants indicates that this protein, the most highly induced in the ERM symbiosis, is a hydrophobin-like effector that participates in the molecular fungal-plant interaction occurring during mycorrhizal formation. Our data demonstrate for the first time the importance of MiSSPs in ERM, although several questions remain open on the cellular localisation of OmSSP1 and its role in symbiosis.
In ECM, hydrophobins likely play an important role in hyphal aggregation during the formation of the extraradical fungal mantle and the Hartig net (Tagu et al., 2001). However, a similar role is unlikely for OmSSP1 because ERM fungi do not form any extraradical hyphal aggregate and individual ERM fungal hyphae take contact with the hydrophilic surface of root epidermal cells (the only cell type colonized in ERM) prior to cell wall penetration (Perotto et al., 2012).
Other intriguing roles have been proposed for fungal hydrophobins in plant-microbe interactions. For example, Whiteford and Spanu (2002) suggested that, by forming a layer on the hyphal surface, hydrophobins could shield potential microbe-associated molecular patterns (MAMPs) recognized by the host and thus avoid plant defense reactions. In Trichoderma, a genus comprising several soil borne plant-growth promoting fungi, adherence to the root surface and root colonization have been suggested to be mediated by hydrophobins (Hermosa et al., 2012). In fact, deletion of the Class I hydrophobin TasHyd1 decreased root attachment and intercellular growth in T. asperellum (Viterbo and Chet, 2006). Similarly we could speculate that a role of OmSSP1 in ERM may be to enhance O. maius attachment to the root surface and to protect the fungal hypha from plant defense compounds.
SC, SD, ElM, and SP planned and designed the research; SC, SD, ElM, SP, AK, EmM, H-RK, YD, and CV-F performed experiments or sequencing, collected, analyzed or interpreted the data; ElM, SP, and FM contributed reagents, materials, analysis; SC, ElM, and SP wrote the manuscript; SD, CV-F, and FM revised it critically. All authors read and approved the final manuscript.
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.
We thank N. Colombano for her help with some of the experiments and F. Rineau for helpful comments on the OmSSP1 biochemical features. SC was supported by a PhD fellowship from the Italian MIUR. The authors acknowledge financial support from local funding of the University of Turin and from the Laboratory of Excellence ARBRE (ANR-11-LABX-0002-01).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.00546/full#supplementary-material
Abbà, S., Khouja, H. R., Martino, E., Archer, D. B., and Perotto, S. (2009). SOD1-targeted gene disruption in the ericoid mycorrhizal fungus Oidiodendron maius reduces conidiation and the capacity for mycorrhization. Mol. Plant Microbe Int. 22, 1412–1421. doi: 10.1094/MPMI-22-11-1412
Akum, F. N., Steinbrenner, J., Biedenkopf, D., Imani, J., and Kogel, K.-H. (2015). The Piriformospora indica effector PIIN_08944 promotes the mutualistic Sebacinalean symbiosis. Front. Plant Sci. 6:906. doi: 10.3389/fpls.2015.00906
Baggerly, K. A., Deng, L., Morris, J. S., and Aldaz, C. M. (2003). Differential expression in SAGE: accounting for normal between-library variation. Bioinformatics 19, 1477–1483. doi: 10.1093/bioinformatics/btg173
Berka, R. M., Grigoriev, I. V., Otillar, R., Salamov, A., Grimwood, J., Reid, I., et al. (2015). Comparative genomic analysis of the thermophilic biomass-degrading fungi Myceliophthora thermophila and Thielavia terrestris. Nat. Biotechnol. 29, 922–927. doi: 10.1038/nbt.1976
Chevenet, F., Brun, C., Banuls, A. L., Jacq, B., and Chisten, R. (2006). TreeDyn: towards dynamic graphics and annotations for analyses of trees. BMC Bioinformatics 7:439. doi: 10.1186/1471-2105-7-439
Daghino, S., Martino, E., and Perotto, S. (2016). Model systems to unravel the molecular mechanisms of heavy metal tolerance in the ericoid mycorrhizal symbiosis. Mycorrhiza 26, 263–274. doi: 10.1007/s00572-015-0675-y
Dereeper, A., Guignon, V., Blanc, G., Audic, S., Buffet, S., Chevenet, F., et al. (2008). Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucl. Acids Res. 1, W465–W469. doi: 10.1093/nar/gkn180
Doré, J., Perraud, M., Dieryckx, C., Kohler, A., Morin, E., Henrissat, B., et al. (2015). Comparative genomics, proteomics and transcriptomics give new insight into the exoproteome of the basidiomycete Hebeloma cylindrosporum and its involvement in ectomycorrhizal symbiosis. New Phytol. 208, 1169–1187. doi: 10.1111/nph.13546
Eisenberg, D., Schwarz, E., Komaromy, M., and Wall, R. (1984). Analysis of membrane and surface protein sequences with the hydrophobic moment plot. J. Mol. Biol. 179, 125–142. doi: 10.1016/0022-2836(84)90309-7
Fiorilli, V., Belmondo, S., Khouja, H.-R., Abb,à, S., Faccio, A., Daghino, S., et al. (2016). RiPEIP1, a gene from the arbuscular mycorrhizal fungus Rhizophagus irregularis, is preferentially expressed in planta and may be involved in root colonization. Mycorrhiza 26, 609–621. doi: 10.1007/s00572-016-0697-0
Giraldo, M. C., Dagdas, Y. F., Gupta, Y. K., Mentlak, T. A., Yi, M., Martinez-Rocha, A. L., et al. (2013). Two distinct secretion systems facilitate tissue invasion by the rice blast fungus Magnaporthe oryzae. Nature Comm. 4:1996. doi: 10.1038/ncomms2996
Hacquard, S., Joly, D. L., Lin, Y.-C., Tisserant, E., Feau, N., Delaruelle, C., et al. (2012). A comprehensive analysis of genes encoding small secreted proteins identifies candidate effectors in Melampsora larici-populina (poplar leaf rust). Mol. Plant-Microbe Int. 25, 279–293. doi: 10.1094/MPMI-09-11-0238
Hammer, Ø., Harper, D. A. T., and Ryan, P. D. (2001). PAST - Paleontological statistics software. package for education and data analysis. Palaeont. Electr. 4:9. Available online at: http://palaeo-electronica.org/2001_1/past/issue1_01.htm
Jacobs, K. A., Collins-Racie, L. A., Colbert, M., Duckett, M., Golden-Fleet, M., Kelleher, K., et al. (1997). A genetic selection for isolating cDNAs encoding secreted proteins. Gene 198, 289–296. doi: 10.1016/S0378-1119(97)00330-2
Kamel, L., Tang, N., Malbreil, M., San Clemente, H., Le Marquer, M., Roux, C., et al. (2017). The comparison of expressed candidate secreted proteins from two arbuscular mycorrhizal fungi unravels common and specific molecular tools to invade different host plants. Front. Plant Sci. 8:124. doi: 10.3389/fpls.2017.00124
Kim, K.-T., Jeon, J., Choi, J., Cheong, K., Song, H., Choi, G., et al. (2016). Kingdom-wide analysis of fungal small secreted proteins (SSPs) reveals their potential role in host association. Front. Plant Sci. 7:186. doi: 10.3389/fpls.2016.00186
Kohler, A., Kuo, A., Nagy, L. G., Morin, E., Barry, K. W., Buscot, F., et al. (2015). Convergent losses of decay mechanisms and rapid turnover of symbiosis genes in mycorrhizal mutualists. Nat. Gen. 47, 410–415. doi: 10.1038/ng.3223
Kwan, A. H., Winefield, R. D., Sunde, M., Matthews, J. M., Haverkamp, R. G., Templeton, M. D., et al. (2006). Structural basis for rodlet assembly in fungal hydrophobins. Proc. Nat. Acad. Sci. U.S.A. 103, 3621–3626. doi: 10.1073/pnas.0505704103
Linder, M. B., Szilvay, G. R., Nakari-Setälä, T., and Penttilä, M. E. (2005). Hydrophobins: the protein-amphiphiles of filamentous fungi. FEMS Microbiol. Rev. 29, 877–896. doi: 10.1016/j.femsre.2005.01.004
Lo Presti, L., Lanver, D., Schweizer, G., Tanaka, S., Liang, L., Tollot, M., et al. (2015). Fungal effectors and plant susceptibility. Ann. Rev. Plant Biol. 66, 513–545. doi: 10.1146/annurev-arplant-043014-114623
Martin, F., Aerts, A., Ahrén, D., Brun, A., Danchin, E. G., Duchaussoy, F., et al. (2008). The genome of Laccaria bicolor provides insights into mycorrhizal symbiosis. Nature 452, 88–92. doi: 10.1038/nature06556
Martino, E., Morin, E., Grelet, G. A., Kuo, A., Kohler, A., Daghino, S., et al. (2018). Comparative genomics and transcriptomics depict ericoid mycorrhizal fungi as versatile saprotrophs and plant mutualists. New Phytol. 217, 1213–1229. doi: 10.1111/nph.14974
Martino, E., Murat, C., Vallino, M., Bena, A., Perotto, S., and Spanu, P. (2007). Imaging mycorrhizal fungal transformants that express EGFP during ericoid endosymbiosis. Curr. Genet. 52, 65–75. doi: 10.1007/s00294-007-0139-9
Martino, E., Turnau, K., Girlanda, M., Bonfante, P., and Perotto, S. (2000). Ericoid mycorrhizal fungi from heavy metal polluted soils: their identification and growth in the presence of zinc ions. Mycol. Res. 104, 338–344. doi: 10.1017/S0953756299001252
Pellegrin, C., Morin, E., Martin, F. M., and Veneault-Fourrey, C. (2015). Comparative analysis of secretomes from ectomycorrhizal fungi with an emphasis on small-secreted proteins. Front. Microbiol. 6:1278. doi: 10.3389/fmicb.2015.01278
Perotto, S., Martino, E., Abb,à, S., and Vallino, M. (2012). “14 Genetic diversity and functional aspects of ericoid mycorrhizal fungi,” in Fungal Associations, ed B. Hock (Berlin, Heidelberg: Springer), 255–285.
Plett, J. M., Daguerre, Y., Wittulsky, S., Vayssières, A., Deveau, A., Melton, S. J., et al. (2014). Effector MiSSP7 of the mutualistic fungus Laccaria bicolor stabilizes the Populus JAZ6 protein and represses jasmonic acid (JA) responsive genes. Proc. Nat. Acad. Sci. U.S.A. 111, 8299–8304. doi: 10.1073/pnas.1322671111
Plett, J. M., Gibon, J., Kohler, A., Duffy, K., Hoegger, P. J., Velagapudi, R., et al. (2012). Phylogenetic, genomic organization and expression analysis of hydrophobin genes in the ectomycorrhizal basidiomycete Laccaria bicolor. Fungal Gen. Biol. 49, 199–209. doi: 10.1016/j.fgb.2012.01.002
Plett, J. M., Kemppainen, M., Kale, S. D., Kohler, A., Legué, V., Brun, A., et al. (2011). A secreted effector protein of Laccaria bicolor is required for symbiosis development. Curr. Biol. 21, 1197–1203. doi: 10.1016/j.cub.2011.05.033
Punt, P. J., Oliver, R. P., Dingemanse, M. A., Pouwels, P. H., and van den Hondel, C. A. (1987). Transformation of Aspergillus based on the hygromycin B resistance marker from Escherichia coli. Gene 56, 117–124. doi: 10.1016/0378-1119(87)90164-8
Rice, A. V., and Currah, R. S. (2006). “Oidiodendron maius: saprobe in sphagnum peat, mutualist in ericaceous roots?” in Microbial Roots Endophytes, eds B. J. E. Shulz, C. J. C. Boyle, and T. N. Sieber (Berlin: Springer-Verlag), 227–246.
Rineau, F., Lmalem, H., Ahren, D., Shah, F., Johansson, T., Coninx, L., et al. (2017). Comparative genomics and expression levels of hydrophobins from eight mycorrhizal genomes. Mycorrhiza 27, 383–396. doi: 10.1007/s00572-016-0758-4
Rovenich, H., Boshoven, J. C., and Thomma, B. P. (2014). Filamentous pathogen effector functions: of pathogens, hosts and microbiomes. Curr. Opin. Plant Biol. 20, 96–103. doi: 10.1016/j.pbi.2014.05.001
Ruytinx, J., Martino, E., Rozpadek, P., Daghino, S., Turnau, K., Colpaert, J., et al. (2016). “Homeostasis of trace elements in mycorrhizal fungi,” in Molecular Mycorrhizal Symbiosis, ed F. Martin (Hoboken, NJ: John Wiley & Sons, Inc.), 277–298.
Salvioli, A., Ghignone, S., Novero, M., Navazio, L., Bagnaresi, P., and Bonfante, P. (2016). Symbiosis with an endobacterium increases the fitness of a mycorrhizal fungus, raising its bioenergetic potential. ISME J. 10, 130–144. doi: 10.1038/ismej.2015.91
Scherlach, K., Graupner, K., and Hertweck, C. (2013). Molecular bacteria-fungi interactions: effects on environment, food, and medicine. Ann. Rev. Microbiol. 67, 375–397. doi: 10.1146/annurev-micro-092412-155702
Seidl-Seiboth, V., Gruber, S., Sezerman, U., Schwecke, T., Albayrak, A., Neuhof, T., et al. (2011). Novel hydrophobins from Trichoderma define a new hydrophobin subclass: protein properties, evolution, regulation and processing. J. Mol. Evol. 72, 339–351. doi: 10.1007/s00239-011-9438-3
Selin, C., de Kievit, T. R., Belmonte, M. F., and Fernando, W. G. (2016). Elucidating the role of effectors in plant-fungal interactions: progress and challenges. Front. Microbiol. 7:600. doi: 10.3389/fmicb.2016.00600
Simossis, V. A., and Heringa, J. (2005). PRALINE: a multiple sequence alignment toolbox that integrates homology-extended and secondary structure information. Nucl. Acids Res. 33, W289–W294. doi: 10.1093/nar/gki390
Tagu, D., De Bellis, R., Balestrini, R., De Vries, O. M. H., Piccoli, G., Stocchi, V., et al. (2001). Immunolocalization of hydrophobin HYDPt-1 from the ectomycorrhizal basidiomycete Pisolithus tinctorius during colonization of Eucalyptus globulus roots. New Phytol. 149, 127–135. doi: 10.1046/j.1469-8137.2001.00009.x
Talbot, J. M., Bruns, T. D., Smith, D. P., Branco, S., Glassman, S. I., Erlandson, S., et al. (2013). Independent roles of ectomycorrhizal and saprotrophic communities in soil organic matter decomposition. Soil Biol. Biochem. 57, 282–291. doi: 10.1016/j.soilbio.2012.10.004
Tang, N., San Clemente, H., Roy, S., Bécard, G., Zhao, B., and Roux, C. (2016). A survey of the gene repertoire of Gigaspora rosea unravels conserved features among Glomeromycota for obligate biotrophy. Front. Microbiol. 7:233. doi: 10.3389/fmicb.2016.00233
Tian, C., Beeson, W. T., Iavarone, A. T., Sun, J., Marletta, M. A., Cate, J. H., et al. (2009). Systems analysis of plant cell wall degradation by the model filamentous fungus Neurospora crassa. Proc. Nat. Acad. Sci. U.S.A. 106, 22157–22162. doi: 10.1073/pnas.0906810106
Tisserant, E., Malbreil, M., Kuo, A., Kohler, A., Symeonidi, A., Balestrini, R., et al. (2013). Genome of an arbuscular mycorrhizal fungus provides insight into the oldest plant symbiosis. Proc. Nat. Acad. Sci. U.S.A. 110, 20117–20122. doi: 10.1073/pnas.1313452110
Toll-Riera, M., Radó-Trilla, N., Martys, F., and Albà, M. M. (2012). Role of low-complexity sequences in the formation of novel protein coding sequences. Mol. Biol. Evol. 29, 883–886. doi: 10.1093/molbev/msr263
Tsuzuki, S., Handa, Y., Takeda, N., and Kawaguchi, M. (2016). Strigolactone-induced putative secreted protein 1 is required for the establishment of symbiosis by the arbuscular mycorrhizal fungus Rhizophagus irregularis. Mol. Plant Microbe Interact. 29, 277–286. doi: 10.1094/MPMI-10-15-0234-R
Valette, N., Benoit-Gelber, I., Falco, M. D., Wiebenga, A., de Vries, R. P., Gelhaye, E., et al. (2017). Secretion of small proteins is species-specific within Aspergillus sp. Microb. Biotech. 10, 323–329. doi: 10.1111/1751-7915.12361
van Esse, H. P., van't Klooster, J. W., Bolton, M. D., Yadeta, K. A., van Baarlen, P., Boeren, S., et al. (2008). The Cladosporium fulvum virulence protein Avr2 inhibits host proteases required for basal defense. Plant Cell Online 20, 1948–1963. doi: 10.1105/tpc.108.059394
Villarreal-Ruiz, L., Anderson, I. C., and Alexander, I. J. (2004). Interaction between an isolate from the Hymenoscyphus ericae aggregate and roots of Pinus and Vaccinium. New Phytol. 164, 183–192. doi: 10.1111/j.1469-8137.2004.01167.x
Vincent, D., Kohler, A., Claverol, S., Solier, E., Joets, J., Gibon, J., et al. (2012). Secretome of the free-living mycelium from the ectomycorrhizal Basidiomycete Laccaria bicolor. J. Proteome Res. 11, 157–171. doi: 10.1021/pr200895f
Viterbo, A., and Chet, I. (2006). TasHyd1, a new hydrophobin gene from the biocontrol agent Trichoderma asperellum, is involved in plant root colonization. Mol. Plant Pathol. 7, 249–258. doi: 10.1111/j.1364-3703.2006.00335.x
Vrålstad, T., Myhre, E., and Schumacher, T. (2002). Molecular diversity and phylogenetic affinities of symbiotic root-associated ascomycetes of the Helotiales in burnt and metal polluted habitats. New Phytol. 155, 131–148. doi: 10.1046/j.1469-8137.2002.00444.x
Wang, S., Boevink, P. C., Welsh, L., Zhang, R., Whisson, S. C., and Birch, P. R. J. (2017). Delivery of cytoplasmic and apoplastic effectors from Phytophthora infestans haustoria by distinct secretion pathways. New Phytol. 216, 205–215. doi: 10.1111/nph.14696
Keywords: ericoid mycorrhiza, Oidiodendron maius, small secreted proteins, hydrophobins, homologous recombination
Citation: Casarrubia S, Daghino S, Kohler A, Morin E, Khouja H-R, Daguerre Y, Veneault-Fourrey C, Martin FM, Perotto S and Martino E (2018) The Hydrophobin-Like OmSSP1 May Be an Effector in the Ericoid Mycorrhizal Symbiosis. Front. Plant Sci. 9:546. doi: 10.3389/fpls.2018.00546
Received: 31 October 2017; Accepted: 09 April 2018;
Published: 01 May 2018.
Edited by:Uta Paszkowski, University of Cambridge, United Kingdom
Reviewed by:Yao-Cheng Lin, Academia Sinica, Taiwan
Kevin Garcia, South Dakota State University, United States
Copyright © 2018 Casarrubia, Daghino, Kohler, Morin, Khouja, Daguerre, Veneault-Fourrey, Martin, Perotto and Martino. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Elena Martino, email@example.com