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.

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.

Functional genomics gives new insights into the ectomycorrhizal degradation of chitin.

Ectomycorrhizal (EcM) fungi play a crucial role in the mineral nitrogen (N) nutrition of their host trees. While it has been proposed that several EcM species also mobilize organic N, studies reporting the EcM ability to degrade N-containing polymers, such as chitin, remain scarce. Here, we assessed the capacity of a representative collection of 16 EcM species to acquire 15 N from 15 N-chitin. In addition, we combined genomics and transcriptomics to identify pathways involved in exogenous chitin degradation between these fungal strains. Boletus edulis, Imleria badia, Suillus luteus, and Hebeloma cylindrosporum efficiently mobilized N from exogenous chitin. EcM genomes primarily contained genes encoding for the direct hydrolysis of chitin. Further, we found a significant relationship between the capacity of EcM fungi to assimilate organic N from chitin and their genomic and transcriptomic potentials for chitin degradation. These findings demonstrate that certain EcM fungal species depolymerize chitin using hydrolytic mechanisms and that endochitinases, but not exochitinases, represent the enzymatic bottleneck of chitin degradation. Finally, this study shows that the degradation of exogenous chitin by EcM fungi might be a key functional trait of nutrient cycling in forests dominated by EcM fungi.

A Remarkable Expansion of Oligopeptide Transporter Genes in Rust Fungi (Pucciniales) Suggests a Specialization in Nutrient Acquisition for Obligate Biotrophy.

Nutrient acquisition by rust fungi during their biotrophic growth has been assigned to a few transporters expressed in haustorial infection structures. We performed a comparative genomic analysis of all transporter genes (hereafter termed transportome) classified according to the Transporter Classification Database, focusing specifically on rust fungi (order Pucciniales) versus other species in the Dikarya. We also surveyed expression of transporter genes in the poplar rust fungus for which transcriptomics data are available across the whole life cycle. Despite a significant increase in gene number, rust fungi presented a reduced transportome compared with most fungi in the Dikarya. However, a few transporter families in the subclass Porters showed significant expansions. Notably, three metal transport-related families involved in the import, export, and sequestration of metals were expanded in Pucciniales and expressed at various stages of the rust life cycle, suggesting a tight regulation of metal homeostasis. The most remarkable gene expansion in the Pucciniales was observed for the oligopeptide transporter (OPT) family, with 25 genes on average compared with seven to 14 genes in the other surveyed taxonomical ranks. A phylogenetic analysis showed several specific expansion events at the root of the order Pucciniales with subsequent expansions in rust taxonomical families. The OPT genes showed dynamic expression patterns along the rust life cycle and more particularly during infection of the poplar host tree, suggesting a possible specialization for the acquisition of nitrogen and sulfur through the transport of oligopeptides from the host during biotrophic growth.

Genomic and transcriptomic characterization of the Collimonas quorum sensing genes and regulon.

Collimonads are well-adapted to nutrient-poor environments. They are known to hydrolyse chitin, produce antifungal metabolites, weather minerals, and are effective biocontrol agents protecting plants from fungal diseases. The production of N-acyl homoserine lactones (AHLs) was suggested to be a conserved trait of collimonads, but little is known about the genes that underlie this production or the genes that are controlled by AHLs. To improve our understanding of the role of AHLs in the ecology of collimonads, we carried out transcriptomic analyses, combined with chemical and functional assays, on strain Collimonas pratensis PMB3(1). The main AHLs produced by this strain were identified as 3-hydroxy-hexa- and octa-noyl-homoserine lactone. Genome analysis permitted to identify putative genes coding for the autoinducer synthase (colI) and cognate transcriptional regulator (colR). The ability to produce AHLs was lost in ΔcolI and ΔcolR mutants. Functional assays revealed that the two mutants metabolized glucose, formate, oxalate, and leucine better than the wild-type (WT) strain. Transcriptome sequencing analyses revealed an up-regulation of different metabolic pathways and of motility in the QS-mutants compared to the WT strain. Overall, our results provide insights into the role of the AHL-dependent regulation system of Collimonas in environment colonization, metabolism readjustment, and microbial interactions.

Phylogenomics and Comparative Genomics Highlight Specific Genetic Features in Ganoderma Species.

The Ganoderma species in Polyporales are ecologically and economically relevant wood decayers used in traditional medicine, but their genomic traits are still poorly documented. In the present study, we carried out a phylogenomic and comparative genomic analyses to better understand the genetic blueprint of this fungal lineage. We investigated seven Ganoderma genomes, including three new genomes, G. australe, G. leucocontextum, and G. lingzhi. The size of the newly sequenced genomes ranged from 60.34 to 84.27 Mb and they encoded 15,007 to 20,460 genes. A total of 58 species, including 40 white-rot fungi, 11 brown-rot fungi, four ectomycorrhizal fungi, one endophyte fungus, and two pathogens in Basidiomycota, were used for phylogenomic analyses based on 143 single-copy genes. It confirmed that Ganoderma species belong to the core polyporoid clade. Comparing to the other selected species, the genomes of the Ganoderma species encoded a larger set of genes involved in terpene metabolism and coding for secreted proteins (CAZymes, lipases, proteases and SSPs). Of note, G. australe has the largest genome size with no obvious genome wide duplication, but showed transposable elements (TEs) expansion and the largest set of terpene gene clusters, suggesting a high ability to produce terpenoids for medicinal treatment. G. australe also encoded the largest set of proteins containing domains for cytochrome P450s, heterokaryon incompatibility and major facilitator families. Besides, the size of G. australe secretome is the largest, including CAZymes (AA9, GH18, A01A), proteases G01, and lipases GGGX, which may enhance the catabolism of cell wall carbohydrates, proteins, and fats during hosts colonization. The current genomic resource will be used to develop further biotechnology and medicinal applications, together with ecological studies of the Ganoderma species.

The mineral weathering ability of Collimonas pratensis PMB3(1) involves a Malleobactin-mediated iron acquisition system.

Mineral weathering by microorganisms is considered to occur through a succession of mechanisms based on acidification and chelation. While the role of acidification is established, the role of siderophores is difficult to disentangle from the effect of the acidification. We took advantage of the ability of strain Collimonas pratensis PMB3(1) to weather minerals but not to acidify depending on the carbon source to address the role of siderophores in mineral weathering. We identified a single non-ribosomal peptide synthetase (NRPS) responsible for siderophore biosynthesis in the PMB3(1) genome. By combining iron-chelating assays, targeted mutagenesis and chemical analyses (HPLC and LC-ESI-HRMS), we identified the siderophore produced as malleobactin X and how its production depends on the concentration of available iron. Comparison with the genome sequences of other collimonads evidenced that malleobactin production seems to be a relatively conserved functional trait, though some collimonads harboured other siderophore synthesis systems. We also revealed by comparing the wild-type strain and its mutant impaired in the production of malleobactin that the ability to produce this siderophore is essential to allow the dissolution of hematite under non-acidifying conditions. This study represents the first characterization of the siderophore produced by collimonads and its role in mineral weathering.

Evolutionary innovations through gain and loss of genes in the ectomycorrhizal Boletales.

We aimed to identify genomic traits of transitions to ectomycorrhizal ecology within the Boletales by comparing the genomes of 21 symbiotrophic species with their saprotrophic brown-rot relatives. Gene duplication rate is constant along the backbone of Boletales phylogeny with large loss events in several lineages, while gene family expansion sharply increased in the late Miocene, mostly in the Boletaceae. Ectomycorrhizal Boletales have a reduced set of plant cell-wall-degrading enzymes (PCWDEs) compared with their brown-rot relatives. However, the various lineages retain distinct sets of PCWDEs, suggesting that, over their evolutionary history, symbiotic Boletales have become functionally diverse. A smaller PCWDE repertoire was found in Sclerodermatineae. The gene repertoire of several lignocellulose oxidoreductases (e.g. laccases) is similar in brown-rot and ectomycorrhizal species, suggesting that symbiotic Boletales are capable of mild lignocellulose decomposition. Transposable element (TE) proliferation contributed to the higher evolutionary rate of genes encoding effector-like small secreted proteins, proteases, and lipases. On the other hand, we showed that the loss of secreted CAZymes was not related to TE activity but to DNA decay. This study provides novel insights on our understanding of the mechanisms influencing the evolutionary diversification of symbiotic boletes.

Evolutionary transition to the ectomycorrhizal habit in the genomes of a hyperdiverse lineage of mushroom-forming fungi.

The ectomycorrhizal (ECM) symbiosis has independently evolved from diverse types of saprotrophic ancestors. In this study, we seek to identify genomic signatures of the transition to the ECM habit within the hyperdiverse Russulaceae. We present comparative analyses of the genomic architecture and the total and secreted gene repertoires of 18 species across the order Russulales, of which 13 are newly sequenced, including a representative of a saprotrophic member of Russulaceae, Gloeopeniophorella convolvens. The genomes of ECM Russulaceae are characterized by a loss of genes for plant cell wall-degrading enzymes (PCWDEs), an expansion of genome size through increased transposable element (TE) content, a reduction in secondary metabolism clusters, and an association of small secreted proteins (SSPs) with TE 'nests', or dense aggregations of TEs. Some PCWDEs have been retained or even expanded, mostly in a species-specific manner. The genome of G. convolvens possesses some characteristics of ECM genomes (e.g. loss of some PCWDEs, TE expansion, reduction in secondary metabolism clusters). Functional specialization in ECM decomposition may drive diversification. Accelerated gene evolution predates the evolution of the ECM habit, indicating that changes in genome architecture and gene content may be necessary to prime the evolutionary switch.

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.

Draft Genome Sequences of the Black Truffles Tuber brumale Vittad. and Tuber indicum Cook & Massee.

Tuber brumale and Tuber indicum (Pezizomycetes) are two edible black truffles establishing ectomycorrhizal symbiosis with trees and shrubs. T. brumale is ubiquitous in Europe, and T. indicum is mainly found in China. Here, we present the draft genome sequences of T. brumale and T. indicum.

Desert truffle genomes reveal their reproductive modes and new insights into plant-fungal interaction and ectendomycorrhizal lifestyle.

Desert truffles are edible hypogeous fungi forming ectendomycorrhizal symbiosis with plants of Cistaceae family. Knowledge about the reproductive modes of these fungi and the molecular mechanisms driving the ectendomycorrhizal interaction is lacking. Genomes of the highly appreciated edible desert truffles Terfezia claveryi Chatin and Tirmania nivea Trappe have been sequenced and compared with other Pezizomycetes. Transcriptomes of T. claveryi × Helianthemum almeriense mycorrhiza from well-watered and drought-stressed plants, when intracellular colonizations is promoted, were investigated. We have identified the fungal genes related to sexual reproduction in desert truffles and desert-truffles-specific genomic and secretomic features with respect to other Pezizomycetes, such as the expansion of a large set of gene families with unknown Pfam domains and a number of species or desert-truffle-specific small secreted proteins differentially regulated in symbiosis. A core set of plant genes, including carbohydrate, lipid-metabolism, and defence-related genes, differentially expressed in mycorrhiza under both conditions was found. Our results highlight the singularities of desert truffles with respect to other mycorrhizal fungi while providing a first glimpse on plant and fungal determinants involved in ecto to endo symbiotic switch that occurs in desert truffle under dry conditions.

Large-scale genome sequencing of mycorrhizal fungi provides insights into the early evolution of symbiotic traits.

Mycorrhizal fungi are mutualists that play crucial roles in nutrient acquisition in terrestrial ecosystems. Mycorrhizal symbioses arose repeatedly across multiple lineages of Mucoromycotina, Ascomycota, and Basidiomycota. Considerable variation exists in the capacity of mycorrhizal fungi to acquire carbon from soil organic matter. Here, we present a combined analysis of 135 fungal genomes from 73 saprotrophic, endophytic and pathogenic species, and 62 mycorrhizal species, including 29 new mycorrhizal genomes. This study samples ecologically dominant fungal guilds for which there were previously no symbiotic genomes available, including ectomycorrhizal Russulales, Thelephorales and Cantharellales. Our analyses show that transitions from saprotrophy to symbiosis involve (1) widespread losses of degrading enzymes acting on lignin and cellulose, (2) co-option of genes present in saprotrophic ancestors to fulfill new symbiotic functions, (3) diversification of novel, lineage-specific symbiosis-induced genes, (4) proliferation of transposable elements and (5) divergent genetic innovations underlying the convergent origins of the ectomycorrhizal guild.

Conserved white-rot enzymatic mechanism for wood decay in the Basidiomycota genus Pycnoporus.

White-rot (WR) fungi are pivotal decomposers of dead organic matter in forest ecosystems and typically use a large array of hydrolytic and oxidative enzymes to deconstruct lignocellulose. However, the extent of lignin and cellulose degradation may vary between species and wood type. Here, we combined comparative genomics, transcriptomics and secretome proteomics to identify conserved enzymatic signatures at the onset of wood-decaying activity within the Basidiomycota genus Pycnoporus. We observed a strong conservation in the genome structures and the repertoires of protein-coding genes across the four Pycnoporus species described to date, despite the species having distinct geographic distributions. We further analysed the early response of P. cinnabarinus, P. coccineus and P. sanguineus to diverse (ligno)-cellulosic substrates. We identified a conserved set of enzymes mobilized by the three species for breaking down cellulose, hemicellulose and pectin. The co-occurrence in the exo-proteomes of H2O2-producing enzymes with H2O2-consuming enzymes was a common feature of the three species, although each enzymatic partner displayed independent transcriptional regulation. Finally, cellobiose dehydrogenase-coding genes were systematically co-regulated with at least one AA9 lytic polysaccharide monooxygenase gene, indicative of enzymatic synergy in vivo. This study highlights a conserved core white-rot fungal enzymatic mechanism behind the wood-decaying process.

Modulation of Plant and Fungal Gene Expression Upon Cd Exposure and Symbiosis in Ericoid Mycorrhizal Vaccinium myrtillus.

The success of Ericaceae in stressful habitats enriched in heavy metals has been ascribed to the distinctive abilities of their mycorrhizal fungal partners to withstand heavy metal stress and to enhance metal tolerance in the host plant. Whereas heavy metal tolerance has been extensively investigated in some ericoid mycorrhizal (ERM) fungi, the molecular and cellular mechanisms that extend tolerance to the host plant are currently unknown. Here, we show a reduced Cd content in Cd-exposed mycorrhizal roots of Vaccinium myrtillus colonized by a metal tolerant isolate of the fungus Oidiodendron maius as compared to non-mycorrhizal roots. To better understand this phenotype, we applied Next Generation Sequencing technologies to analyze gene expression in V. myrtillus and O. maius Zn grown under normal and Cd-stressed conditions, in the free living and in the mycorrhizal status. The results clearly showed that Cd had a stronger impact on plant gene expression than symbiosis, whereas fungal gene expression was mainly regulated by symbiosis. The higher abundance of transcripts coding for stress related proteins in non-mycorrhizal roots may be related to the higher Cd content. Regulated plant metal transporters have been identified that may play a role in reducing Cd content in mycorrhizal roots exposed to this metal.

At the nexus of three kingdoms: the genome of the mycorrhizal fungus Gigaspora margarita provides insights into plant, endobacterial and fungal interactions.

As members of the plant microbiota, arbuscular mycorrhizal fungi (AMF, Glomeromycotina) symbiotically colonize plant roots. AMF also possess their own microbiota, hosting some uncultivable endobacteria. Ongoing research has revealed the genetics underlying plant responses to colonization by AMF, but the fungal side of the relationship remains in the dark. Here, we sequenced the genome of Gigaspora margarita, a member of the Gigasporaceae in an early diverging group of the Glomeromycotina. In contrast to other AMF, G. margarita may host distinct endobacterial populations and possesses the largest fungal genome so far annotated (773.104 Mbp), with more than 64% transposable elements. Other unique traits of the G. margarita genome include the expansion of genes for inorganic phosphate metabolism, the presence of genes for production of secondary metabolites and a considerable number of potential horizontal gene transfer events. The sequencing of G. margarita genome reveals the importance of its immune system, shedding light on the evolutionary pathways that allowed early diverging fungi to interact with both plants and bacteria.

Transcriptome Analysis of Apple Leaves Infected by the Rust Fungus Gymnosporangium yamadae at Two Sporulation Stages.

Apple rust disease caused by Gymnosporangium yamadae is one of the major threats to apple orchards. In this study, dual RNA-seq analysis was conducted to simultaneously monitor gene expression profiles of G. yamadae and infected apple leaves during the formation of rust spermogonia and aecia. The molecular mechanisms underlying this compatible interaction at 10 and 30 days postinoculation (dpi) indicate a significant reaction from the host plant and comprise detoxication pathways at the earliest stage and the induction of secondary metabolism pathways at 30 dpi. Such host reactions have been previously reported in other rust pathosystems and may represent a general reaction to rust infection. G. yamadae transcript profiling indicates a conserved genetic program in spermogonia and aecia that is shared with other rust fungi, whereas secretome prediction reveals the presence of specific secreted candidate effector proteins expressed during apple infection. Unexpectedly, the survey of fungal unigenes in the transcriptome assemblies of inoculated and mock-inoculated apple leaves reveals that G. yamadae infection may modify the fungal community composition in the apple phyllosphere at 30 dpi. Collectively, our results provide novel insights into the compatible apple-G. yamadae interaction and advance the knowledge of this heteroecious demicyclic rust fungus.

Fungal ecological strategies reflected in gene transcription - a case study of two litter decomposers.

Microbial communities interplay with their environment through their functional traits that can be a response or an effect on the environment. Here, we explore how a functional trait-the decomposition of organic matter, can be addressed based on genetic markers and how the expression of these markers reflect ecological strategies of two fungal litter decomposer Gymnopus androsaceus and Chalara longipes. We sequenced the genomes of these two fungi, as well as their transcriptomes at different steps of Pinus sylvestris needles decomposition in microcosms. Our results highlighted that if the gene content of the two species could indicate similar potential decomposition abilities, the expression levels of specific gene families belonging to the glycoside hydrolase category reflected contrasting ecological strategies. Actually, C. longipes, the weaker decomposer in this experiment, turned out to have a high content of genes involved in cell wall polysaccharides decomposition but low expression levels, reflecting a versatile ecology compare to the more competitive G. androsaceus with high expression levels of keystone functional genes. Thus, we established that sequential expression of genes coding for different components of the decomposer machinery indicated adaptation to chemical changes in the substrate as decomposition progressed.

Comparative transcriptomics of Gymnosporangium spp. teliospores reveals a conserved genetic program at this specific stage of the rust fungal life cycle.

BACKGROUND: Gymnosporangium spp. are fungal plant pathogens causing rust disease and most of them are known to infect two different host plants (heteroecious) with four spore stages (demicyclic). In the present study, we sequenced the transcriptome of G. japonicum teliospores on its host plant Juniperus chinensis and we performed comparison to the transcriptomes of G. yamadae and G. asiaticum at the same life stage, that happens in the same host but on different organs. RESULTS: Functional annotation for the three Gymnosporangium species showed the expression of a conserved genetic program with the top abundant cellular categories corresponding to energy, translation and signal transduction processes, indicating that this life stage is particularly active. Moreover, the survey of predicted secretomes in the three Gymnosporangium transcriptomes revealed shared and specific genes encoding carbohydrate active enzymes and secreted proteins of unknown function that could represent candidate pathogenesis effectors. A transcript encoding a hemicellulase of the glycoside hydrolase 26 family, previously identified in other rust fungi, was particularly highly expressed suggesting a general role in rust fungi. The comparison between the transcriptomes of the three Gymnosporangium spp. and selected Pucciniales species in different taxonomical families allowed to identify lineage-specific protein families that may relate to the biology of teliospores in rust fungi. Among clustered gene families, 205, 200 and 152 proteins were specifically identified in G. japonicum, G. yamadae and G. asiaticum, respectively, including candidate effectors expressed in teliospores. CONCLUSIONS: This comprehensive comparative transcriptomics study of three Gymnosporangium spp. identified gene functions and metabolic pathways particularly expressed in teliospores, a stage of the life cycle that is mostly overlooked in rust fungi. Secreted protein encoding transcripts expressed in teliospores may reveal new candidate effectors related to pathogenesis. Although this spore stage is not involved in host plant infection but in the production of basidiospores infecting plants in the Amygdaloideae, we speculate that candidate effectors may be expressed as early as the teliospore stage for preparing further infection by basidiospores.

Multi-omic analyses of exogenous nutrient bag decomposition by the black morel Morchella importuna reveal sustained carbon acquisition and transferring.

The black morel (Morchella importuna Kuo, O'Donnell and Volk) was once an uncultivable wild mushroom, until the development of exogenous nutrient bag (ENB), making its agricultural production quite feasible and stable. To date, how the nutritional acquisition of the morel mycelium is fulfilled to trigger its fruiting remains unknown. To investigate the mechanisms involved in ENB decomposition, the genome of a cultivable morel strain (M. importuna SCYDJ1-A1) was sequenced and the genes coding for the decay apparatus were identified. Expression of the encoded carbohydrate-active enzymes (CAZymes) was then analyzed by metatranscriptomics and metaproteomics in combination with biochemical assays. The results show that a diverse set of hydrolytic and redox CAZymes secreted by the morel mycelium is the main force driving the substrate decomposition. Plant polysaccharides such as starch and cellulose present in ENB substrate (wheat grains plus rice husks) were rapidly degraded, whereas triglycerides were accumulated initially and consumed later. ENB decomposition led to a rapid increase in the organic carbon content in the surface soil of the mushroom bed, which was thereafter consumed during morel fruiting. In contrast to the high carbon consumption, no significant acquisition of nitrogen was observed. Our findings contribute to an increasingly detailed portrait of molecular features triggering morel fruiting.