Metatranscriptomics sheds light on the links between the functional traits of fungal guilds and ecological processes in forest soil ecosystems

Soil fungi belonging to different functional guilds, such as saprotrophs, pathogens, and mycorrhizal symbionts, play key roles in forest ecosystems. To date, no study has compared the actual gene expression of these guilds in different forest soils. We used metatranscriptomics to study the competition for organic resources by these fungal groups in boreal, temperate, and Mediterranean forest soils. Using a dedicated mRNA annotation pipeline combined with the JGI MycoCosm database, we compared the transcripts of these three fungal guilds, targeting enzymes involved in C- and N mobilization from plant and microbial cell walls. Genes encoding enzymes involved in the degradation of plant cell walls were expressed at a higher level in saprotrophic fungi than in ectomycorrhizal and pathogenic fungi. However, ectomycorrhizal and saprotrophic fungi showed similarly high expression levels of genes encoding enzymes involved in fungal cell wall degradation. Transcripts for N-related transporters were more highly expressed in ectomycorrhizal fungi than in other groups. We showed that ectomycorrhizal and saprotrophic fungi compete for N in soil organic matter, suggesting that their interactions could decelerate C cycling. Metatranscriptomics provides a unique tool to test controversial ecological hypotheses and to better understand the underlying ecological processes involved in soil functioning and carbon stabilization.

Multilevel analysis between Physcomitrium patens and Mortierellaceae endophytes explores potential long-standing interaction among land plants and fungi

The model moss species Physcomitrium patens has long been used for studying divergence of land plants spanning from bryophytes to angiosperms. In addition to its phylogenetic relationships, the limited number of differential tissues, and comparable morphology to the earliest embryophytes provide a system to represent basic plant architecture. Based on plant-fungal interactions today, it is hypothesized these kingdoms have a long-standing relationship, predating plant terrestrialization. Mortierellaceae have origins diverging from other land fungi paralleling bryophyte divergence, are related to arbuscular mycorrhizal fungi but are free-living, observed to interact with plants, and can be found in moss microbiomes globally. Due to their parallel origins, we assess here how two Mortierellaceae species, Linnemannia elongata and Benniella erionia, interact with P. patens in coculture. We also assess how Mollicute-related or Burkholderia-related endobacterial symbionts (MRE or BRE) of these fungi impact plant response. Coculture interactions are investigated through high-throughput phenomics, microscopy, RNA-sequencing, differential expression profiling, gene ontology enrichment, and comparisons among 99 other P. patens transcriptomic studies. Here we present new high-throughput approaches for measuring P. patens growth, identify novel expression of over 800 genes that are not expressed on traditional agar media, identify subtle interactions between P. patens and Mortierellaceae, and observe changes to plant-fungal interactions dependent on whether MRE or BRE are present. Our study provides insights into how plants and fungal partners may have interacted based on their communications observed today as well as identifying L. elongata and B. erionia as modern fungal endophytes with P. patens.

Populus MYC2 orchestrates root transcriptional reprogramming of defence pathway to impair Laccaria bicolor ectomycorrhizal development

The jasmonic acid (JA) signalling pathway plays an important role in the establishment of the ectomycorrhizal symbiosis. The Laccaria bicolor effector MiSSP7 stabilizes JA corepressor JAZ6, thereby inhibiting the activity of Populus MYC2 transcription factors. Although the role of MYC2 in orchestrating plant defences against pathogens is well established, its exact contribution to ECM symbiosis remains unclear. This information is crucial for understanding the balance between plant immunity and symbiotic relationships. Transgenic poplars overexpressing or silencing for the two paralogues of MYC2 transcription factor (MYC2s) were produced, and their ability to establish ectomycorrhiza was assessed. Transcriptomics and DNA affinity purification sequencing were performed. MYC2s overexpression led to a decrease in fungal colonization, whereas its silencing increased it. The enrichment of terpene synthase genes in the MYC2-regulated gene set suggests a complex interplay between the host monoterpenes and fungal growth. Several root monoterpenes have been identified as inhibitors of fungal growth and ECM symbiosis. Our results highlight the significance of poplar MYC2s and terpenes in mutualistic symbiosis by controlling root fungal colonization. We identified poplar genes which direct or indirect control by MYC2 is required for ECM establishment. These findings deepen our understanding of the molecular mechanisms underlying ECM symbiosis.

The reference genome and abiotic stress responses of the model perennial grass Brachypodium sylvaticum

Perennial grasses are important forage crops and emerging biomass crops and have the potential to be more sustainable grain crops. However, most perennial grass crops are difficult experimental subjects due to their large size, difficult genetics, and/or their recalcitrance to transformation. Thus, a tractable model perennial grass could be used to rapidly make discoveries that can be translated to perennial grass crops. Brachypodium sylvaticum has the potential to serve as such a model because of its small size, rapid generation time, simple genetics, and transformability. Here, we provide a high-quality genome assembly and annotation for B. sylvaticum, an essential resource for a modern model system. In addition, we conducted transcriptomic studies under 4 abiotic stresses (water, heat, salt, and freezing). Our results indicate that crowns are more responsive to freezing than leaves which may help them overwinter. We observed extensive transcriptional responses with varying temporal dynamics to all abiotic stresses, including classic heat-responsive genes. These results can be used to form testable hypotheses about how perennial grasses respond to these stresses. Taken together, these results will allow B. sylvaticum to serve as a truly tractable perennial model system.

Convergent reductive evolution and host adaptation in Mycoavidus bacterial endosymbionts of Mortierellaceae fungi

Intimate associations between fungi and intracellular bacterial endosymbionts are becoming increasingly well understood. Phylogenetic analyses demonstrate that bacterial endosymbionts of Mucoromycota fungi are related either to free-living Burkholderia or Mollicutes species. The so-called Burkholderia-related endosymbionts or BRE comprise Mycoavidus, Mycetohabitans and Candidatus Glomeribacter gigasporarum. These endosymbionts are marked by genome contraction thought to be associated with intracellular selection. However, the conclusions drawn thus far are based on a very small subset of endosymbiont genomes, and the mechanisms leading to genome streamlining are not well understood. The purpose of this study was to better understand how intracellular existence shapes Mycoavidus and BRE functionally at the genome level. To this end we generated and analyzed 14 novel draft genomes for Mycoavidus living within the hyphae of Mortierellomycotina fungi. We found that our novel Mycoavidus genomes were significantly reduced compared to free-living Burkholderiales relatives. Using a genome-scale phylogenetic approach including the novel and available existing genomes of Mycoavidus, we show that the genus is an assemblage composed of two independently derived lineages including three well supported clades of Mycoavidus. Using a comparative genomic approach, we shed light on the functional implications of genome reduction, documenting shared and unique gene loss patterns between the three Mycoavidus clades. We found that many endosymbiont isolates demonstrate patterns of vertical transmission and host-specificity, but others are present in phylogenetically disparate hosts. We discuss how reductive evolution and host specificity reflect convergent adaptation to the intrahyphal selective landscape, and commonalities of eukaryotic endosymbiont genome evolution.

Spray-induced gene silencing to identify powdery mildew gene targets and processes for powdery mildew control

Spray-induced gene silencing (SIGS) is an emerging tool for crop pest protection. It utilizes exogenously applied double-stranded RNA to specifically reduce pest target gene expression using endogenous RNA interference machinery. In this study, SIGS methods were developed and optimized for powdery mildew fungi, which are widespread obligate biotrophic fungi that infect agricultural crops, using the known azole-fungicide target cytochrome P450 51 (CYP51) in the Golovinomyces orontii-Arabidopsis thaliana pathosystem. Additional screening resulted in the identification of conserved gene targets and processes important to powdery mildew proliferation: apoptosis-antagonizing transcription factor in essential cellular metabolism and stress response; lipid catabolism genes lipase a, lipase 1, and acetyl-CoA oxidase in energy production; and genes involved in manipulation of the plant host via abscisic acid metabolism (9-cis-epoxycarotenoid dioxygenase, xanthoxin dehydrogenase, and a putative abscisic acid G-protein coupled receptor) and secretion of the effector protein, effector candidate 2. Powdery mildew is the dominant disease impacting grapes and extensive powdery mildew resistance to applied fungicides has been reported. We therefore developed SIGS for the Erysiphe necator-Vitis vinifera system and tested six successful targets identified using the G. orontii-A. thaliana system. For all targets tested, a similar reduction in powdery mildew disease was observed between systems. This indicates screening of broadly conserved targets in the G. orontii-A. thaliana pathosystem identifies targets and processes for the successful control of other powdery mildew fungi. The efficacy of SIGS on powdery mildew fungi makes SIGS an exciting prospect for commercial powdery mildew control.

High-quality RNA extraction and the regulation of genes encoding cellulosomes are correlated with growth stage in anaerobic fungi

Anaerobic fungi produce biomass-degrading enzymes and natural products that are important to harness for several biotechnology applications. Although progress has been made in the development of methods for extracting nucleic acids for genomic and transcriptomic sequencing of these fungi, most studies are limited in that they do not sample multiple fungal growth phases in batch culture. In this study, we establish a method to harvest RNA from fungal monocultures and fungal-methanogen co-cultures, and also determine an optimal time frame for high-quality RNA extraction from anaerobic fungi. Based on RNA quality and quantity targets, the optimal time frame in which to harvest anaerobic fungal monocultures and fungal-methanogen co-cultures for RNA extraction was 2-5 days of growth post-inoculation. When grown on cellulose, the fungal strain Anaeromyces robustus cocultivated with the methanogen Methanobacterium bryantii upregulated genes encoding fungal carbohydrate-active enzymes and other cellulosome components relative to fungal monocultures during this time frame, but expression patterns changed at 24-hour intervals throughout the fungal growth phase. These results demonstrate the importance of establishing methods to extract high-quality RNA from anaerobic fungi at multiple time points during batch cultivation.

DNA hypomethylation of the host tree impairs interaction with mutualistic ectomycorrhizal fungus

Ectomycorrhizas are an intrinsic component of tree nutrition and responses to environmental variations. How epigenetic mechanisms might regulate these mutualistic interactions is unknown. By manipulating the level of expression of the chromatin remodeler DECREASE IN DNA METHYLATION 1 (DDM1) and two demethylases DEMETER-LIKE (DML) in Populus tremula × Populus alba lines, we examined how host DNA methylation modulates multiple parameters of the responses to root colonization with the mutualistic fungus Laccaria bicolor. We compared the ectomycorrhizas formed between transgenic and wild-type (WT) trees and analyzed their methylomes and transcriptomes. The poplar lines displaying lower mycorrhiza formation rate corresponded to hypomethylated overexpressing DML or RNAi-ddm1 lines. We found 86 genes and 288 transposable elements (TEs) differentially methylated between WT and hypomethylated lines (common to both OX-dml and RNAi-ddm1) and 120 genes/1441 TEs in the fungal genome suggesting a host-induced remodeling of the fungal methylome. Hypomethylated poplar lines displayed 205 differentially expressed genes (cis and trans effects) in common with 17 being differentially methylated (cis). Our findings suggest a central role of host and fungal DNA methylation in the ability to form ectomycorrhizas including not only poplar genes involved in root initiation, ethylene and jasmonate-mediated pathways, and immune response but also terpenoid metabolism.

Major proliferation of transposable elements shaped the genome of the soybean rust pathogen Phakopsora pachyrhizi

With >7000 species the order of rust fungi has a disproportionately large impact on agriculture, horticulture, forestry and foreign ecosystems. The infectious spores are typically dikaryotic, a feature unique to fungi in which two haploid nuclei reside in the same cell. A key example is Phakopsora pachyrhizi, the causal agent of Asian soybean rust disease, one of the world's most economically damaging agricultural diseases. Despite P. pachyrhizi's impact, the exceptional size and complexity of its genome prevented generation of an accurate genome assembly. Here, we sequence three independent P. pachyrhizi genomes and uncover a genome up to 1.25 Gb comprising two haplotypes with a transposable element (TE) content of ~93%. We study the incursion and dominant impact of these TEs on the genome and show how they have a key impact on various processes such as host range adaptation, stress responses and genetic plasticity.

Acquisition of host-derived carbon in biomass of 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 colonised 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.

Speciation Underpinned by Unexpected Molecular Diversity in the Mycorrhizal Fungal Genus Pisolithus

The mutualistic ectomycorrhizal (ECM) fungal genus Pisolithus comprises 19 species defined to date which colonize the roots of >50 hosts worldwide suggesting that substantial genomic and functional evolution occurred during speciation. To better understand this intra-genus variation, we undertook a comparative multi-omic study of nine Pisolithus species sampled from North America, South America, Asia, and Australasia. We found that there was a small core set of genes common to all species (13%), and that these genes were more likely to be significantly regulated during symbiosis with a host than accessory or species-specific genes. Thus, the genetic "toolbox" foundational to the symbiotic lifestyle in this genus is small. Transposable elements were located significantly closer to gene classes including effector-like small secreted proteins (SSPs). Poorly conserved SSPs were more likely to be induced by symbiosis, suggesting that they may be a class of protein that tune host specificity. The Pisolithus gene repertoire is characterized by divergent CAZyme profiles when compared with other fungi, both symbiotic and saprotrophic. This was driven by differences in enzymes associated with symbiotic sugar processing, although metabolomic analysis suggest that neither copy number nor expression of these genes is sufficient to predict sugar capture from a host plant or its metabolism in fungal hyphae. Our results demonstrate that intra-genus genomic and functional diversity within ECM fungi is greater than previously thought, underlining the importance of continued comparative studies within the fungal tree of life to refine our focus on pathways and evolutionary processes foundational to this symbiotic lifestyle.

Co‑cultivation of anaerobic fungi with Clostridium acetobutylicum bolsters butyrate and butanol production from cellulose and lignocellulose

A system for co-cultivation of anaerobic fungi with anaerobic bacteria was established based on lactate cross-feeding to produce butyrate and butanol from plant biomass. Several co-culture formulations were assembled that consisted of anaerobic fungi (Anaeromyces robustus, Neocallimastix californiae, or Caecomyces churrovis) with the bacterium Clostridium acetobutylicum. Co-cultures were grown simultaneously (e.g., 'one pot'), and compared to cultures where bacteria were cultured in fungal hydrolysate sequentially. Fungal hydrolysis of lignocellulose resulted in 7-11 mM amounts of glucose and xylose, as well as acetate, formate, ethanol, and lactate to support clostridial growth. Under these conditions, one-stage simultaneous co-culture of anaerobic fungi with C. acetobutylicum promoted the production of butyrate up to 30 mM. Alternatively, two-stage growth slightly promoted solventogenesis and elevated butanol levels (∼4-9 mM). Transcriptional regulation in the two-stage growth condition indicated that this cultivation method may decrease the time required to reach solventogenesis and induce the expression of cellulose-degrading genes in C. acetobutylicum due to relieved carbon-catabolite repression. Overall, this study demonstrates a proof of concept for biobutanol and bio-butyrate production from lignocellulose using an anaerobic fungal-bacterial co-culture system.

Genetic mapping and prediction of flowering time and plant height in a maize Stiff Stalk MAGIC population

The Stiff Stalk heterotic pool is a foundation of US maize seed parent germplasm and has been heavily utilized by both public and private maize breeders since its inception in the 1930's. Flowering time and plant height are critical characteristics for both inbred parents and their test crossed hybrid progeny. To study these traits, a six parent multiparent advanced generation intercross (MAGIC) population was developed including maize inbred lines B73, B84, PHB47 (B37 type), LH145 (B14 type), PHJ40 (novel early Stiff Stalk), and NKH8431 (B73/B14 type). A set of 779 doubled haploid lines were evaluated for flowering time and plant height in two field replicates in 2016 and 2017, and a subset of 689 and 561 doubled haploid lines were crossed to two testers, respectively, and evaluated as hybrids in two locations in 2018 and 2019 using an incomplete block design. Markers were derived from a Practical Haplotype Graph built from the founder whole genome assemblies and genotype-by-sequencing and exome capture-based sequencing of the population. Genetic mapping utilizing an update to R/qtl2 revealed differing profiles of significant loci for both traits between 635 of the DH lines and two sets of 570 and 471 derived hybrids. Genomic prediction was used to test the feasibility of predicting hybrid phenotypes based on the per se data. Predictive abilities were highest on direct models trained using the data they would predict (0.55 to 0.63), and indirect models trained using per se data to predict hybrid traits had slightly lower predictive abilities (0.49 to 0.55). Overall, this finding is consistent with the overlapping and non-overlapping significant QTL found within the per se and hybrid populations and suggests that selections for phenology traits can be made effectively on doubled haploid lines before hybrid data is available.

Co‑cultivation of the anaerobic fungus Caecomyces churrovis with Methanobacterium bryantii enhances transcription of carbohydrate binding modules, dockerins, and pyruvate formate lyases on specific substrates

Anaerobic fungi and methanogenic archaea are two classes of microorganisms found in the rumen microbiome that metabolically interact during lignocellulose breakdown. Here, stable synthetic co-cultures of the anaerobic fungus Caecomyces churrovis and the methanogen Methanobacterium bryantii (not native to the rumen) were formed, demonstrating that microbes from different environments can be paired based on metabolic ties. Transcriptional and metabolic changes induced by methanogen co-culture were evaluated in C. churrovis across a variety of substrates to identify mechanisms that impact biomass breakdown and sugar uptake. A high-quality genome of C. churrovis was obtained and annotated, which is the first sequenced genome of a non-rhizoid-forming anaerobic fungus. C. churrovis possess an abundance of CAZymes and carbohydrate binding modules and, in agreement with previous studies of early-diverging fungal lineages, N6-methyldeoxyadenine (6mA) was associated with transcriptionally active genes. Co-culture with the methanogen increased overall transcription of CAZymes, carbohydrate binding modules, and dockerin domains in co-cultures grown on both lignocellulose and cellulose and caused upregulation of genes coding associated enzymatic machinery including carbohydrate binding modules in family 18 and dockerin domains across multiple growth substrates relative to C. churrovis monoculture. Two other fungal strains grown on a reed canary grass substrate in co-culture with the same methanogen also exhibited high log2-fold change values for upregulation of genes encoding carbohydrate binding modules in families 1 and 18. Transcriptional upregulation indicated that co-culture of the C. churrovis strain with a methanogen may enhance pyruvate formate lyase (PFL) function for growth on xylan and fructose and production of bottleneck enzymes in sugar utilization pathways, further supporting the hypothesis that co-culture with a methanogen may enhance certain fungal metabolic functions. Upregulation of CBM18 may play a role in fungal-methanogen physical associations and fungal cell wall development and remodeling.

The Ancient Salicoid Genome Duplication Event: A Platform for Reconstruction of De Novo Gene Evolution in Populus trichocarpa

Orphan genes are characteristic genomic features that have no detectable homology to genes in any other species and represent an important attribute of genome evolution as sources of novel genetic functions. Here, we identified 445 genes specific to Populus trichocarpa. Of these, we performed deeper reconstruction of 13 orphan genes to provide evidence of de novo gene evolution. Populus and its sister genera Salix are particularly well suited for the study of orphan gene evolution because of the Salicoid whole-genome duplication event which resulted in highly syntenic sister chromosomal segments across the Salicaceae. We leveraged this genomic feature to reconstruct de novo gene evolution from intergenera, interspecies, and intragenomic perspectives by comparing the syntenic regions within the P. trichocarpa reference, then P. deltoides, and finally Salix purpurea. Furthermore, we demonstrated that 86.5% of the putative orphan genes had evidence of transcription. Additionally, we also utilized the Populus genome-wide association mapping panel, a collection of 1,084 undomesticated P. trichocarpa genotypes to further determine putative regulatory networks of orphan genes using expression quantitative trait loci (eQTL) mapping. Functional enrichment of these eQTL subnetworks identified common biological themes associated with orphan genes such as response to stress and defense response. We also identify a putative cis-element for a de novo gene and leverage conserved synteny to describe evolution of a putative transcription factor binding site. Overall, 45% of orphan genes were captured in trans-eQTL networks.

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.

Overexpression of a Prefoldin β subunit gene reduces biomass recalcitrance in the bioenergy crop Populus

Prefoldin (PFD) is a group II chaperonin that is ubiquitously present in the eukaryotic kingdom. Six subunits (PFD1-6) form a jellyfish-like heterohexameric PFD complex and function in protein folding and cytoskeleton organization. However, little is known about its function in plant cell wall-related processes. Here, we report the functional characterization of a PFD gene from Populus deltoides, designated as PdPFD2.2. There are two copies of PFD2 in Populus, and PdPFD2.2 was ubiquitously expressed with high transcript abundance in the cambial region. PdPFD2.2 can physically interact with DELLA protein RGA1_8g, and its subcellular localization is affected by the interaction. In P. deltoides transgenic plants overexpressing PdPFD2.2, the lignin syringyl/guaiacyl ratio was increased, but cellulose content and crystallinity index were unchanged. In addition, the total released sugar (glucose and xylose) amounts were increased by 7.6% and 6.1%, respectively, in two transgenic lines. Transcriptomic and metabolomic analyses revealed that secondary metabolic pathways, including lignin and flavonoid biosynthesis, were affected by overexpressing PdPFD2.2. A total of eight hub transcription factors (TFs) were identified based on TF binding sites of differentially expressed genes in Populus transgenic plants overexpressing PdPFD2.2. In addition, several known cell wall-related TFs, such as MYB3, MYB4, MYB7, TT8 and XND1, were affected by overexpression of PdPFD2.2. These results suggest that overexpression of PdPFD2.2 can reduce biomass recalcitrance and PdPFD2.2 is a promising target for genetic engineering to improve feedstock characteristics to enhance biofuel conversion and reduce the cost of lignocellulosic biofuel production.

The White-Rot Basidiomycete Dichomitus squalens Shows Highly Specific Transcriptional Response to Lignocellulose-Related Aromatic Compounds

Lignocellulosic plant biomass is an important feedstock for bio-based economy. In particular, it is an abundant renewable source of aromatic compounds, which are present as part of lignin, as side-groups of xylan and pectin, and in other forms, such as tannins. As filamentous fungi are the main organisms that modify and degrade lignocellulose, they have developed a versatile metabolism to convert the aromatic compounds that are toxic at relatively low concentrations to less toxic ones. During this process, fungi form metabolites some of which represent high-value platform chemicals or important chemical building blocks, such as benzoic, vanillic, and protocatechuic acid. Especially basidiomycete white-rot fungi with unique ability to degrade the recalcitrant lignin polymer are expected to perform highly efficient enzymatic conversions of aromatic compounds, thus having huge potential for biotechnological exploitation. However, the aromatic metabolism of basidiomycete fungi is poorly studied and knowledge on them is based on the combined results of studies in variety of species, leaving the overall picture in each organism unclear. Dichomitus squalens is an efficiently wood-degrading white-rot basidiomycete that produces a diverse set of extracellular enzymes targeted for lignocellulose degradation, including oxidative enzymes that act on lignin. Our recent study showed that several intra- and extracellular aromatic compounds were produced when D. squalens was cultivated on spruce wood, indicating also versatile aromatic metabolic abilities for this species. In order to provide the first molecular level systematic insight into the conversion of plant biomass derived aromatic compounds by basidiomycete fungi, we analyzed the transcriptomes of D. squalens when grown with 10 different lignocellulose-related aromatic monomers. Significant differences for example with respect to the expression of lignocellulose degradation related genes, but also putative genes encoding transporters and catabolic pathway genes were observed between the cultivations supplemented with the different aromatic compounds. The results demonstrate that the transcriptional response of D. squalens is highly dependent on the specific aromatic compounds present suggesting that instead of a common regulatory system, fine-tuned regulation is needed for aromatic metabolism.

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.

Comparative genomics of Coniophora olivacea reveals different patterns of genome expansion in Boletales

Background

Coniophora olivacea is a basidiomycete fungus belonging to the order Boletales that produces brown-rot decay on dead wood of conifers. The Boletales order comprises a diverse group of species including saprotrophs and ectomycorrhizal fungi that show important differences in genome size.

Results

In this study we report the 39.07-megabase (Mb) draft genome assembly and annotation of C. olivacea. A total of 14,928 genes were annotated, including 470 putatively secreted proteins enriched in functions involved in lignocellulose degradation. Using similarity clustering and protein structure prediction we identified a new family of 10 putative lytic polysaccharide monooxygenase genes. This family is conserved in basidiomycota and lacks of previous functional annotation. Further analyses showed that C. olivacea has a low repetitive genome, with 2.91% of repeats and a restrained content of transposable elements (TEs). The annotation of TEs in four related Boletales yielded important differences in repeat content, ranging from 3.94 to 41.17% of the genome size. The distribution of insertion ages of LTR-retrotransposons showed that differential expansions of these repetitive elements have shaped the genome architecture of Boletales over the last 60 million years.

Conclusions

Coniophora olivacea has a small, compact genome that shows macrosynteny with Coniophora puteana. The functional annotation revealed the enzymatic signature of a canonical brown-rot. The annotation and comparative genomics of transposable elements uncovered their particular contraction in the Coniophora genera, highlighting their role in the differential genome expansions found in Boletales species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4243-z) contains supplementary material, which is available to authorized users.

Genome sequencing reveals complex secondary metabolome in the marine actinomycete Salinispora tropica

Recent fermentation studies have identified actinomycetes of the marine-dwelling genus Salinispora as prolific natural product producers. To further evaluate their biosynthetic potential, we sequenced the 5,183,331-bp S. tropica CNB-440 circular genome and analyzed all identifiable secondary natural product gene clusters. Our analysis shows that S. tropica dedicates a large percentage of its genome ( approximately 9.9%) to natural product assembly, which is greater than previous Streptomyces genome sequences as well as other natural product-producing actinomycetes. The S. tropica genome features polyketide synthase systems of every known formally classified family, nonribosomal peptide synthetases, and several hybrid clusters. Although a few clusters appear to encode molecules previously identified in Streptomyces species, the majority of the 17 biosynthetic loci are novel. Specific chemical information about putative and observed natural product molecules is presented and discussed. In addition, our bioinformatic analysis not only was critical for the structure elucidation of the polyene macrolactam salinilactam A, but its structural analysis aided the genome assembly of the highly repetitive slm loci. This study firmly establishes the genus Salinispora as a rich source of drug-like molecules and importantly reveals the powerful interplay between genomic analysis and traditional natural product isolation studies.