Persistence and plasticity in bacterial gene regulation.

Organisms orchestrate cellular functions through transcription factor (TF) interactions with their target genes, although these regulatory relationships are largely unknown in most species. Here we report a high-throughput approach for characterizing TF-target gene interactions across species and its application to 354 TFs across 48 bacteria, generating 17,000 genome-wide binding maps. This dataset revealed themes of ancient conservation and rapid evolution of regulatory modules. We observed rewiring, where the TF sensing and regulatory role is maintained while the arrangement and identity of target genes diverges, in some cases encoding entirely new functions. We further integrated phenotypic information to define new functional regulatory modules and pathways. Finally, we identified 242 new TF DNA binding motifs, including a 70% increase of known Escherichia coli motifs and the first annotation in Pseudomonas simiae, revealing deep conservation in bacterial promoter architecture. Our method provides a versatile tool for functional characterization of genetic pathways in prokaryotes and eukaryotes.

Identification and recombinant expression of a cutinase from Papiliotrema laurentii that hydrolyzes natural and synthetic polyesters.

Given the multitude of extracellular enzymes at their disposal, many of which are designed to degrade nature's polymers (lignin, cutin, cellulose, etc.), fungi are adept at targeting synthetic polyesters with similar chemical composition. Microbial-influenced deterioration of xenobiotic polymeric surfaces is an area of interest for material scientists as these are important for the conservation of the underlying structural materials. Here, we describe the isolation and characterization of the Papiliotrema laurentii 5307AH (P. laurentii) cutinase, Plcut1. P. laurentii is basidiomycete yeast with the ability to disperse Impranil-DLN (Impranil), a colloidal polyester polyurethane, in agar plates. To test whether the fungal factor involved in this clearing was a secreted enzyme, we screened the ability of P. laurentii culture supernatants to disperse Impranil. Using size exclusion chromatography (SEC), we isolated fractions that contained Impranil-clearing activity. These fractions harbored a single ~22 kD band, which was excised and subjected to peptide sequencing. Homology searches using the peptide sequences identified, revealed that the protein Papla1 543643 (Plcut1) displays similarities to serine esterase and cutinase family of proteins. Biochemical assays using recombinant Plcut1 confirmed that this enzyme has the capability to hydrolyze Impranil, soluble esterase substrates, and apple cutin. Finally, we confirmed the presence of the Plcut1 in culture supernatants using a custom antibody that specifically recognizes this protein. The work shown here supports a major role for the Plcut1 in the fungal degradation of natural polyesters and xenobiotic polymer surfaces.IMPORTANCEFungi play a vital role in the execution of a broad range of biological processes that drive ecosystem function through production of a diverse arsenal of enzymes. However, the universal reactivity of these enzymes is a current problem for the built environment and the undesired degradation of polymeric materials in protective coatings. Here, we report the identification and characterization of a hydrolase from Papiliotrema laurentii 5307AH, an aircraft-derived fungal isolate found colonizing a biodeteriorated polymer-coated surface. We show that P. laurentii secretes a cutinase capable of hydrolyzing soluble esters as well as ester-based compounds forming solid surface coatings. These findings indicate that this fungus plays a significant role in biodeterioration through the production of a cutinase adept at degrading ester-based polymers, some of which form the backbone of protective surface coatings. The work shown here provides insights into the mechanisms employed by fungi to degrade xenobiotic polymers.

Genome-wide role of codon usage on transcription and identification of potential regulators.

Codon usage bias is a fundamental feature of all genomes and plays an important role in determining gene expression levels. The codon usage was thought to influence gene expression mainly due to its impact on translation. Recently, however, codon usage was shown to affect transcription of fungal and mammalian genes, indicating the existence of a gene regulatory phenomenon with unknown mechanism. In Neurospora, codon usage biases strongly correlate with mRNA levels genome-wide, and here we show that the correlation between codon usage and RNA levels is maintained in the nucleus. In addition, codon optimality is tightly correlated with both total and nuclear RNA levels, suggesting that codon usage broadly influences mRNA levels through transcription in a translation-independent manner. A large-scale RNA sequencing-based genetic screen in Neurospora identified 18 candidate factors that when deleted decreased the genome-wide correlation between codon usage and RNA levels and reduced the codon usage effect on gene expression. Most of these factors, such as the H3K36 methyltransferase, are chromatin regulators or transcription factors. Together, our results suggest that the transcriptional effect of codon usage is mediated by multiple transcriptional regulatory mechanisms.

Molecular causes of an evolutionary shift along the parasitism-mutualism continuum in a bacterial symbiont.

Symbiosis with microbes is a ubiquitous phenomenon with a massive impact on all living organisms, shaping the world around us today. Theoretical and experimental studies show that vertical transmission of symbionts leads to the evolution of mutualistic traits, whereas horizontal transmission facilitates the emergence of parasitic features. However, these studies focused on phenotypic data, and we know little about underlying molecular changes at the genomic level. Here, we combined an experimental evolution approach with infection assays, genome resequencing, and global gene expression analysis to study the effect of transmission mode on an obligate intracellular bacterial symbiont. We show that a dramatic shift in the frequency of genetic variants, coupled with major changes in gene expression, allow the symbiont to alter its position in the parasitism-mutualism continuum depending on the mode of between-host transmission. We found that increased parasitism in horizontally transmitted chlamydiae residing in amoebae was a result of processes occurring at the infectious stage of the symbiont's developmental cycle. Specifically, genes involved in energy production required for extracellular survival and the type III secretion system-the symbiont's primary virulence mechanism-were significantly up-regulated. Our results identify the genomic and transcriptional dynamics sufficient to favor parasitic or mutualistic strategies.

Cryptic genetic structure and copy-number variation in the ubiquitous forest symbiotic fungus Cenococcum geophilum.

Ectomycorrhizal (ECM) fungi associated with plants constitute one of the most successful symbiotic interactions in forest ecosystems. ECM support trophic exchanges with host plants and are important factors for the survival and stress resilience of trees. However, ECM clades often harbour morpho-species and cryptic lineages, with weak morphological differentiation. How this relates to intraspecific genome variability and ecological functioning is poorly known. Here, we analysed 16 European isolates of the ascomycete Cenococcum geophilum, an extremely ubiquitous forest symbiotic fungus with no known sexual or asexual spore-forming structures but with a massively enlarged genome. We carried out whole-genome sequencing to identify single-nucleotide polymorphisms. We found no geographic structure at the European scale but divergent lineages within sampling sites. Evidence for recombination was restricted to specific cryptic lineages. Lineage differentiation was supported by extensive copy-number variation. Finally, we confirmed heterothallism with a single MAT1 idiomorph per genome. Synteny analyses of the MAT1 locus revealed substantial rearrangements and a pseudogene of the opposite MAT1 idiomorph. Our study provides the first evidence for substantial genome-wide structural variation, lineage-specific recombination and low continent-wide genetic differentiation in C. geophilum. Our study provides a foundation for targeted analyses of intra-specific functional variation in this major symbiosis.

Gene family expansions and transcriptome signatures uncover fungal adaptations to wood decay.

Because they comprise some of the most efficient wood-decayers, Polyporales fungi impact carbon cycling in forest environment. Despite continuous discoveries on the enzymatic machinery involved in wood decomposition, the vision on their evolutionary adaptation to wood decay and genome diversity remains incomplete. We combined the genome sequence information from 50 Polyporales species, including 26 newly sequenced genomes and sought for genomic and functional adaptations to wood decay through the analysis of genome composition and transcriptome responses to different carbon sources. The genomes of Polyporales from different phylogenetic clades showed poor conservation in macrosynteny, indicative of genome rearrangements. We observed different gene family expansion/contraction histories for plant cell wall degrading enzymes in core polyporoids and phlebioids and captured expansions for genes involved in signalling and regulation in the lineages of white rotters. Furthermore, we identified conserved cupredoxins, thaumatin-like proteins and lytic polysaccharide monooxygenases with a yet uncharacterized appended module as new candidate players in wood decomposition. Given the current need for enzymatic toolkits dedicated to the transformation of renewable carbon sources, the observed genomic diversity among Polyporales strengthens the relevance of mining Polyporales biodiversity to understand the molecular mechanisms of wood decay.

Comparative genomics reveals dynamic genome evolution in host specialist ectomycorrhizal fungi.

While there has been significant progress characterizing the 'symbiotic toolkit' of ectomycorrhizal (ECM) fungi, how host specificity may be encoded into ECM fungal genomes remains poorly understood. We conducted a comparative genomic analysis of ECM fungal host specialists and generalists, focusing on the specialist genus Suillus. Global analyses of genome dynamics across 46 species were assessed, along with targeted analyses of three classes of molecules previously identified as important determinants of host specificity: small secreted proteins (SSPs), secondary metabolites (SMs) and G-protein coupled receptors (GPCRs). Relative to other ECM fungi, including other host specialists, Suillus had highly dynamic genomes including numerous rapidly evolving gene families and many domain expansions and contractions. Targeted analyses supported a role for SMs but not SSPs or GPCRs in Suillus host specificity. Phylogenomic-based ancestral state reconstruction identified Larix as the ancestral host of Suillus, with multiple independent switches between white and red pine hosts. These results suggest that like other defining characteristics of the ECM lifestyle, host specificity is a dynamic process at the genome level. In the case of Suillus, both SMs and pathways involved in the deactivation of reactive oxygen species appear to be strongly associated with enhanced host specificity.

An ectomycorrhizal fungus alters sensitivity to jasmonate, salicylate, gibberellin, and ethylene in host roots.

The phytohormones jasmonate, gibberellin, salicylate, and ethylene regulate an interconnected reprogramming network integrating root development with plant responses against microbes. The establishment of mutualistic ectomycorrhizal symbiosis requires the suppression of plant defense responses against fungi as well as the modification of root architecture and cortical cell wall properties. Here, we investigated the contribution of phytohormones and their crosstalk to the ontogenesis of ectomycorrhizae (ECM) between grey poplar (Populus tremula x alba) roots and the fungus Laccaria bicolor. To obtain the hormonal blueprint of developing ECM, we quantified the concentrations of jasmonates, gibberellins, and salicylate via liquid chromatography-tandem mass spectrometry. Subsequently, we assessed root architecture, mycorrhizal morphology, and gene expression levels (RNA sequencing) in phytohormone-treated poplar lateral roots in the presence or absence of L. bicolor. Salicylic acid accumulated in mid-stage ECM. Exogenous phytohormone treatment affected the fungal colonization rate and/or frequency of Hartig net formation. Colonized lateral roots displayed diminished responsiveness to jasmonate but regulated some genes, implicated in defense and cell wall remodelling, that were specifically differentially expressed after jasmonate treatment. Responses to salicylate, gibberellin, and ethylene were enhanced in ECM. The dynamics of phytohormone accumulation and response suggest that jasmonate, gibberellin, salicylate, and ethylene signalling play multifaceted roles in poplar L. bicolor ectomycorrhizal development.

Genome sequencing of evolved aspergilli populations reveals robust genomes, transversions in A. flavus, and sexual aberrancy in non-homologous end-joining mutants.

BACKGROUND: Aspergillus spp. comprises a very diverse group of lower eukaryotes with a high relevance for industrial applications and clinical implications. These multinucleate species are often cultured for many generations in the laboratory, which can unknowingly propagate hidden genetic mutations. To assess the likelihood of such events, we studied the genome stability of aspergilli by using a combination of mutation accumulation (MA) lines and whole genome sequencing. RESULTS: We sequenced the whole genomes of 30 asexual and 10 sexual MA lines of three Aspergillus species (A. flavus, A. fumigatus and A. nidulans) and estimated that each MA line accumulated mutations for over 4000 mitoses during asexual cycles. We estimated mutation rates of 4.2 × 10-11 (A. flavus), 1.1 × 10-11 (A. fumigatus) and 4.1 × 10-11 (A. nidulans) per site per mitosis, suggesting that the genomes are very robust. Unexpectedly, we found a very high rate of GC → TA transversions only in A. flavus. In parallel, 30 asexual lines of the non-homologous end-joining (NHEJ) mutants of the three species were also allowed to accumulate mutations for the same number of mitoses. Sequencing of these NHEJ MA lines gave an estimated mutation rate of 5.1 × 10-11 (A. flavus), 2.2 × 10-11 (A. fumigatus) and 4.5 × 10-11 (A. nidulans) per base per mitosis, which is slightly higher than in the wild-type strains and some ~ 5-6 times lower than in the yeasts. Additionally, in A. nidulans, we found a NHEJ-dependent interference of the sexual cycle that is independent of the accumulation of mutations. CONCLUSIONS: We present for the first time direct counts of the mutation rate of filamentous fungal species and find that Aspergillus genomes are very robust. Deletion of the NHEJ machinery results in a slight increase in the mutation rate, but at a rate we suggest is still safe to use for biotechnology purposes. Unexpectedly, we found GC→TA transversions predominated only in the species A. flavus, which could be generated by the hepatocarcinogen secondary metabolite aflatoxin. Lastly, a strong effect of the NHEJ mutation in self-crossing was observed and an increase in the mutations of the asexual lines was quantified.

Resolving the Mortierellaceae phylogeny through synthesis of multi-gene phylogenetics and phylogenomics.

Early efforts to classify Mortierellaceae were based on macro- and micromorphology, but sequencing and phylogenetic studies with ribosomal DNA (rDNA) markers have demonstrated conflicting taxonomic groupings and polyphyletic genera. Although some taxonomic confusion in the family has been clarified, rDNA data alone is unable to resolve higher level phylogenetic relationships within Mortierellaceae. In this study, we applied two parallel approaches to resolve the Mortierellaceae phylogeny: low coverage genome (LCG) sequencing and high-throughput, multiplexed targeted amplicon sequencing to generate sequence data for multi-gene phylogenetics. We then combined our datasets to provide a well-supported genome-based phylogeny having broad sampling depth from the amplicon dataset. Resolving the Mortierellaceae phylogeny into monophyletic groups led to the definition of 14 genera, 7 of which are newly proposed. Low-coverage genome sequencing proved to be a relatively cost-effective means of generating a well-resolved phylogeny. The multi-gene phylogenetics approach enabled much greater sampling depth and breadth than the LCG approach, but was unable to resolve higher-level organization of groups. We present this work to resolve some of the taxonomic confusion and provide a genus-level framework to empower future studies on Mortierellaceae diversity, biology, and evolution.

Phylogenomics of Endogonaceae and evolution of mycorrhizas within Mucoromycota.

Endogonales (Mucoromycotina), composed of Endogonaceae and Densosporaceae, is the only known non-Dikarya order with ectomycorrhizal members. They also form mycorrhizal-like association with some nonspermatophyte plants. It has been recently proposed that Endogonales were among the earliest mycorrhizal partners with land plants. It remains unknown whether Endogonales possess genomes with mycorrhizal-lifestyle signatures and whether Endogonales originated around the same time as land plants did. We sampled sporocarp tissue from four Endogonaceae collections and performed shotgun genome sequencing. After binning the metagenome data, we assembled and annotated the Endogonaceae genomes. We performed comparative analysis on plant-cell-wall-degrading enzymes (PCWDEs) and small secreted proteins (SSPs). We inferred phylogenetic placement of Endogonaceae and estimated the ages of Endogonaceae and Endogonales with expanded taxon sampling. Endogonaceae have large genomes with high repeat content, low diversity of PCWDEs, but without elevated SSP/secretome ratios. Dating analysis estimated that Endogonaceae originated in the Permian-Triassic boundary and Endogonales originated in the mid-late Silurian. Mycoplasma-related endobacterium sequences were identified in three Endogonaceae genomes. Endogonaceae genomes possess typical signatures of mycorrhizal lifestyle. The early origin of Endogonales suggests that the mycorrhizal association between Endogonales and plants might have played an important role during the colonization of land by plants.

Draft genome sequence of the Tremellomycetes yeast Papiliotrema laurentii 5307AH, isolated from aircraft.

Papiliotrema laurentii 5307AH was isolated from an aircraft polymer-coated surface. The genome size is 19,510,785 bp with a G + C content of 56%. The genome harbors genes encoding oxygenases, cutinases, lipases, and enzymes for styrene degradation, all of which could play a critical role in survival on xenobiotic surfaces.

Genome Dynamics and Temperature Adaptation During Experimental Evolution of Obligate Intracellular Bacteria.

Evolution experiments with free-living microbes have radically improved our understanding of genome evolution and how microorganisms adapt. Yet there is a paucity of such research focusing on strictly host-associated bacteria, even though they are widespread in nature. Here, we used the Acanthamoeba symbiont Protochlamydia amoebophila, a distant relative of the human pathogen Chlamydia trachomatis and representative of a large group of protist-associated environmental chlamydiae, as a model to study how obligate intracellular symbionts evolve and adapt to elevated temperature, a prerequisite for the pivotal evolutionary leap from protist to endothermic animal hosts. We established 12 replicate populations under two temperatures (20 °C, 30 °C) for 510 bacterial generations (38 months). We then used infectivity assays and pooled whole-genome resequencing to identify any evolved phenotypes and the molecular basis of adaptation in these bacteria. We observed an overall reduction in infectivity of the symbionts evolved at 30 °C, and we identified numerous nonsynonymous mutations and small indels in these symbiont populations, with several variants persisting throughout multiple time points and reaching high frequencies. This suggests that many mutations may have been beneficial and played an adaptive role. Mutated genes within the same temperature regime were more similar than those between temperature regimes. Our results provide insights into the molecular evolution of intracellular bacteria under the constraints of strict host dependance and highly structured populations and suggest that for chlamydial symbionts of protists, temperature adaptation was facilitated through attenuation of symbiont infectivity as a tradeoff to reduce host cell burden.

Acquisition of host-derived carbon in biomass of the ectomycorrhizal fungus Pisolithus microcarpus is correlated to fungal carbon demand and plant defences.

Ectomycorrhizal (ECM) fungi are key players in forest carbon (C) sequestration, receiving a substantial proportion of photosynthetic C from their forest tree hosts in exchange for plant growth-limiting soil nutrients. However, it remains unknown whether the fungus or plant controls the quantum of C in this exchange, nor what mechanisms are involved. Here, we aimed to identify physiological and genetic properties of both partners that influence ECM C transfer. Using a microcosm system, stable isotope tracing, and transcriptomics, we quantified plant-to-fungus C transfer between the host plant Eucalyptus grandis and nine isolates of the ECM fungus Pisolithus microcarpus that range in their mycorrhization potential and investigated fungal growth characteristics and plant and fungal genes that correlated with C acquisition. We found that C acquisition by P. microcarpus correlated positively with both fungal biomass production and the expression of a subset of fungal C metabolism genes. In the plant, C transfer was not positively correlated to the number of colonized root tips, but rather to the expression of defence- and stress-related genes. These findings suggest that C acquisition by ECM fungi involves individual fungal demand for C and defence responses of the host against C drain.

Lignin deconstruction by anaerobic fungi.

Lignocellulose forms plant cell walls, and its three constituent polymers, cellulose, hemicellulose and lignin, represent the largest renewable organic carbon pool in the terrestrial biosphere. Insights into biological lignocellulose deconstruction inform understandings of global carbon sequestration dynamics and provide inspiration for biotechnologies seeking to address the current climate crisis by producing renewable chemicals from plant biomass. Organisms in diverse environments disassemble lignocellulose, and carbohydrate degradation processes are well defined, but biological lignin deconstruction is described only in aerobic systems. It is currently unclear whether anaerobic lignin deconstruction is impossible because of biochemical constraints or, alternatively, has not yet been measured. We applied whole cell-wall nuclear magnetic resonance, gel-permeation chromatography and transcriptome sequencing to interrogate the apparent paradox that anaerobic fungi (Neocallimastigomycetes), well-documented lignocellulose degradation specialists, are unable to modify lignin. We find that Neocallimastigomycetes anaerobically break chemical bonds in grass and hardwood lignins, and we further associate upregulated gene products with the observed lignocellulose deconstruction. These findings alter perceptions of lignin deconstruction by anaerobes and provide opportunities to advance decarbonization biotechnologies that depend on depolymerizing lignocellulose.

Introducing EzTaxon-e: a prokaryotic 16S rRNA gene sequence database with phylotypes that represent uncultured species.

Despite recent advances in commercially optimized identification systems, bacterial identification remains a challenging task in many routine microbiological laboratories, especially in situations where taxonomically novel isolates are involved. The 16S rRNA gene has been used extensively for this task when coupled with a well-curated database, such as EzTaxon, containing sequences of type strains of prokaryotic species with validly published names. Although the EzTaxon database has been widely used for routine identification of prokaryotic isolates, sequences from uncultured prokaryotes have not been considered. Here, the next generation database, named EzTaxon-e, is formally introduced. This new database covers not only species within the formal nomenclatural system but also phylotypes that may represent species in nature. In addition to an identification function based on Basic Local Alignment Search Tool (blast) searches and pairwise global sequence alignments, a new objective method of assessing the degree of completeness in sequencing is proposed. All sequences that are held in the EzTaxon-e database have been subjected to phylogenetic analysis and this has resulted in a complete hierarchical classification system. It is concluded that the EzTaxon-e database provides a useful taxonomic backbone for the identification of cultured and uncultured prokaryotes and offers a valuable means of communication among microbiologists who routinely encounter taxonomically novel isolates. The database and its analytical functions can be found at http://eztaxon-e.ezbiocloud.net/.

Genomic Signatures of a Major Adaptive Event in the Pathogenic Fungus Melampsora larici-populina.

The recent availability of genome-wide sequencing techniques has allowed systematic screening for molecular signatures of adaptation, including in nonmodel organisms. Host-pathogen interactions constitute good models due to the strong selective pressures that they entail. We focused on an adaptive event which affected the poplar rust fungus Melampsora larici-populina when it overcame a resistance gene borne by its host, cultivated poplar. Based on 76 virulent and avirulent isolates framing narrowly the estimated date of the adaptive event, we examined the molecular signatures of selection. Using an array of genome scan methods based on different features of nucleotide diversity, we detected a single locus exhibiting a consistent pattern suggestive of a selective sweep in virulent individuals (excess of differentiation between virulent and avirulent samples, linkage disequilibrium, genotype-phenotype statistical association, and long-range haplotypes). Our study pinpoints a single gene and further a single amino acid replacement which may have allowed the adaptive event. Although our samples are nearly contemporary to the selective sweep, it does not seem to have affected genome diversity further than the immediate vicinity of the causal locus, which can be explained by a soft selective sweep (where selection acts on standing variation) and by the impact of recombination in mitigating the impact of selection. Therefore, it seems that properties of the life cycle of M. larici-populina, which entails both high genetic diversity and outbreeding, has facilitated its adaptation.

Draft Genome Assemblies of Ionic Liquid-Resistant Yarrowia lipolytica PO1f and Its Superior Evolved Strain, YlCW001.

Adaptive laboratory evolution of Yarrowia lipolytica PO1f in the benchmark ionic liquid (IL; 1-ethyl-3-methylimidazolium acetate) produced a superior IL-tolerant microorganism, strain YlCW001. Here, we report the genome sequences of PO1f and YlCW001 to study the robustness of Y. lipolytica and its potential use as a microbial platform for producing fuels and chemicals.

Draft Genome Assemblies of Five Robust Yarrowia lipolytica Strains Exhibiting High Lipid Production, Pentose Sugar Utilization, and Sugar Alcohol Secretion from Undetoxified Lignocellulosic Biomass Hydrolysates.

Screening the genetic diversity of 45 Yarrowia lipolytica strains identified five candidates with unique metabolic capability and robustness in undetoxified switchgrass hydrolysates, including superior lipid production and efficient pentose sugar utilization. Here, we report the genome sequences of these strains to study their robustness and potential to produce fuels and chemicals.

Secretome Analysis from the Ectomycorrhizal Ascomycete Cenococcum geophilum.

Cenococcum geophilum is an ectomycorrhizal fungus with global distribution in numerous habitats and associates with a large range of host species including gymnosperm and angiosperm trees. Moreover, C. geophilum is the unique ectomycorrhizal species within the clade Dothideomycetes, the largest class of Ascomycetes containing predominantly saprotrophic and many devastating phytopathogenic fungi. Recent studies highlight that mycorrhizal fungi, as pathogenic ones, use effectors in form of Small Secreted Proteins (SSPs) as molecular keys to promote symbiosis. In order to better understand the biotic interaction of C. geophilum with its host plants, the goal of this work was to characterize mycorrhiza-induced small-secreted proteins (MiSSPs) that potentially play a role in the ectomycorrhiza formation and functioning of this ecologically very important species. We combined different approaches such as gene expression profiling, genome localization and conservation of MiSSP genes in different C. geophilum strains and closely related species as well as protein subcellular localization studies of potential targets of MiSSPs in interacting plants using in tobacco leaf cells. Gene expression analyses of C. geophilum interacting with Pinus sylvestris (pine) and Populus tremula × Populus alba (poplar) showed that similar sets of genes coding for secreted proteins were up-regulated and only few were specific to each host. Whereas pine induced more carbohydrate active enzymes (CAZymes), the interaction with poplar induced the expression of specific SSPs. We identified a set of 22 MiSSPs, which are located in both, gene-rich, repeat-poor or gene-sparse, repeat-rich regions of the C. geophilum genome, a genome showing a bipartite architecture as seen for some pathogens but not yet for an ectomycorrhizal fungus. Genome re-sequencing data of 15 C. geophilum strains and two close relatives Glonium stellatum and Lepidopterella palustris were used to study sequence conservation of MiSSP-encoding genes. The 22 MiSSPs showed a high presence-absence polymorphism among the studied C. geophilum strains suggesting an evolution through gene gain/gene loss. Finally, we showed that six CgMiSSPs target four distinct sub-cellular compartments such as endoplasmic reticulum, plasma membrane, cytosol and tonoplast. Overall, this work presents a comprehensive analysis of secreted proteins and MiSSPs in different genetic level of C. geophilum opening a valuable resource to future functional analysis.