Persistent drought affects global crop production and is becoming more severe in many parts of the world in recent decades. Deciphering how plants respond to drought will facilitate the development of flexible mitigation strategies. Sorghum bicolor L. Moench (sorghum), a major cereal crop and an emerging bioenergy crop, exhibits remarkable resilience to drought. To better understand the molecular traits that underlie sorghum's remarkable drought tolerance, we undertook a large-scale sorghum gene expression profiling effort, totaling nearly 1500 transcriptome profiles, across a 3-year field study with replicated plots in California's Central Valley. This study included time-resolved gene expression data from roots and leaves of two sorghum genotypes, BTx642 and RTx430, with different pre-flowering and post-flowering drought-tolerance adaptations under control and drought conditions. Quantification of genotype-specific drought tolerance effects was enabled by de novo sequencing, assembly, and annotation of both BTx642 and RTx430 genomes. These reference-quality genomes were used to construct a pangene set for characterizing conserved and genotype-specific expression. By integrating time-resolved transcriptomic responses to drought in the field across three consecutive years, we identified a set of 726 drought-responsive genes that responded similarly in all 3 years of our field study. Functional enrichment analysis identified abiotic stress, secondary cell wall-related processes and metabolism as particularly affected under both types of drought stress. We also found that some glyoxylate cycle pathway genes, including malate synthase and isocitrate lyase, are differentially regulated particularly during post-flowering drought stress, implicating this pathway as potentially important for drought responsiveness. This expansive dataset represents a unique resource for sorghum and drought research communities and provides a methodological framework for the integration of multi-faceted time-resolved transcriptomic datasets.

Morels are iconic macrofungi known for their culinary value and spring emergence. Molecular phylogenetic studies have hitherto failed to elucidate the evolution of the ecological lifestyles of true and false morels, impeding their capacity to resolve longstanding debates regarding the extent of their saprotrophy and the potential for biotrophic associations with plants. In this study, we examined the evolutionary history and molecular innovations of Morchellaceae by generating high-quality genomes from species encompassing all the major clades within this family. We used conserved single-copy orthologs to infer the evolutionary history of true and false morels and related truffles. The genome size, transposable element content, and polysaccharide-degrading enzyme gene content are consistent with an ectomycorrhizal lifestyle for the Morchellaceae truffles Leucangium carthusianum and Kalapuya brunnea, whereas the other clades, including all true morels (Morchella spp.), exhibited genomic signatures of saprotrophic habits characterized by a highly conserved set of genes encoding plant biomass-related degrading enzymes. Morchellaceae species were found to be predominantly heterothallic, with either the MAT1-1 or MAT1-2 loci; however, the occurrence of colocalized mating-type idiomorphs, indicative of homothallism, was observed in M. rufobrunnea, M. peruviana, L. carthusianum, and the outgroup taxon Gyromitra esculenta. This study revealed a wealth of largely undiscovered genomic traits, including saprotrophic potential within Morchellaceae and a lack of biotrophic markers, and contributes to our understanding of the intricate evolutionary trajectories of the modes of nutrition in soil fungi.

Effector secretion is crucial for root endophytes to establish and protect their ecological niche. We used time-resolved transcriptomics to monitor effector gene expression dynamics in two closely related Sebacinales, Serendipita indica and Serendipita vermifera, during symbiosis with three plant species, competition with the phytopathogenic fungus Bipolaris sorokiniana, and cooperation with root-associated bacteria. We observed increased effector gene expression in response to biotic interactions, particularly with plants, indicating their importance in host colonization. Some effectors responded to both plants and microbes, suggesting dual roles in intermicrobial competition and plant-microbe interactions. A subset of putative antimicrobial effectors, including a GH18-CBM5 chitinase, was induced exclusively by microbes. Functional analyses of this chitinase revealed its antimicrobial and plant-protective properties. We conclude that dynamic effector gene expression underpins the ability of Sebacinales to thrive in diverse ecological niches with a single fungal chitinase contributing substantially to niche defense.

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.

Summary 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.

Wood-degrading white-rot fungi can efficiently degrade all plant biomass components, but the molecular mechanisms behind the degradation of plant polysaccharides remain poorly understood. For example, the gene sets and expression levels induced by the plant polysaccharide-derived monosaccharides in white-rot fungi do not reflect those induced by crude plant biomass substrates. To explore the molecular response of the white-rot fungus Dichomitus squalens to plant-derived oligo- and polysaccharides, we investigated the transcriptomes from mono- and dikaryotic strains of the fungus on 10 substrates and compared the expression of carbohydrate-active enzyme-encoding genes to that previously reported for different monosaccharides and cellobiose. Our results revealed that in D. squalens, a robust response to cellulose leads to its effective depolymerization, with an orthologue of the ascomycete Trichoderma reesei ACE3 likely acting as a central transcriptional regulator. The conserved response between cellulose and cellobiose further confirms cellobiose as the main cellulase inducer in D. squalens. Surprisingly, despite low abundance of pectin in the natural wood substrate of D. squalens, we identified polygalacturonic acid as a major inducer of a broad-targeted pectinolytic response including pectinase, pectin-related sugar transporter and catabolism genes, and four fungal specific transcription factors. This indicates that D. squalens has not only maintained its ability to degrade minor polysaccharide components in its biotope, but also a regulatory system spanning from extracellular degradation to metabolic conversion. Our study contributes to a deeper understanding of the molecular mechanisms behind white-rot fungal plant polysaccharide degradation and provides leads for functional studies of potential transcriptional regulators in basidiomycetes.

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.

Associations between plants and beneficial root-endophytic fungi enhance plant performance by improving nutrient uptake, abiotic stress tolerance and disease resistance. To successfully colonize different host plants and defend their host niche against competing microbes, but also to cooperate with beneficial bacterial members of the microbiota, root endophytes such as Sebacinales secrete a multitude of tightly regulated effector-proteins and carbohydrate-active enzymes. However, the functions, specificity, and regulation of these proteins remain poorly understood. In this study, we employ time-resolved transcriptomics to analyse the gene expression profiles of two Sebacinales members interacting with organisms from different kingdoms of life. We identified crucial genes for plant colonization and intermicrobial competition, including a fungal GH18-CBM5 chitinase specifically upregulated in response to the phytopathogenic fungus Bipolaris sorokiniana. This chitinase protects the plant hosts against the pathogen, reducing fungal biomass and disease symptoms in barley and Arabidopsis thaliana. Our findings shed light on interaction partner specific gene expression in Sebacinales endophytes, with potential applications in enhancing plant health and resilience.

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.

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.

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.

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.

Abstract 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.

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.

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.

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.

SUMMARY 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. Powdery mildews, widespread obligate biotrophic fungi infect agricultural crops including wheat, barley, cucurbits, grapevine, and ornamentals such as roses. In this study, SIGS methods were developed and optimized for powdery mildews using the known azole-fungicide target CYP51 and 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; lipase a, lipase 1 , and acetyl-CoA oxidase in lipid catabolism; 9-cis-epoxycarotenoid dioxygenase, xanthoxin dehydrogenase , and a putative abscisic acid G-protein coupled receptor predicted to function in manipulation of the plant hormone abscisic acid; and the secreted effector EC2 . Powdery mildew is the dominant disease impacting grapes and extensive powdery mildew resistance to applied fungicides has been reported. Therefore, we 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 mildews. The flexibility, specificity, reduced environmental and health risks, and rapid transition from the bench to the field make SIGS an exciting prospect for commercial powdery mildew control.

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 (similar to 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.

Microalgae that are of interest for biofuel production must be able to tolerate environmental changes that occur in outdoor cultivation systems. While algal cultures may experience daily temperature fluctuations and seasonal environmental changes, the underlying mechanisms that control and regulate physiological responses and adaptation to environmental pressures are largely unknown. Systems-level characterization enabled by func-tional genomics can help identify biochemical pathways that promote stability and productivity of algae in various environmental conditions. Monoraphidium minutum 26B-AM, a freshwater green microalga, was identi-fied as a top performer in biomass production in winter season screens. We sequenced the genome of M. minutum 26B-AM and applied our multi-omics pipeline to profile this high potential strain under high salt and cold temperature perturbations. Through comparative analysis, including other green algae in the class Chlor-ophyceae, we identified gene families unique to the genus Monoraphidium, including a desaturase that has been linked to cold tolerance in plants. We observed that osmolytes, such as trehalose, proline and betaine, accu-mulate under salt stress, coinciding with upregulation of genes involved in biosynthesis of these metabolites. From the genome annotation, we reconstructed a metabolic model to provide a detailed map of the metabolic pathways and can be used to simulate growth and reaction fluxes. This multi-omics analysis provides a foun-dation to explore algal strain potential for biofuel applications, guides strain engineering, and expands our understanding of metabolic and regulatory mechanisms of algae in applied systems.

Small RNAs (sRNAs) are known to regulate pathogenic plant-microbe interactions. Emerging evidence from the study of these model systems suggests that microRNAs (miRNAs) can be translocated between microbes and plants to facilitate symbiosis. The roles of sRNAs in mutualistic mycorrhizal fungal interactions, however, are largely unknown. In this study, we characterized miRNAs encoded by the ectomycorrhizal fungus Pisolithus microcarpus and investigated their expression during mutualistic interaction with Eucalyptus grandis. Using sRNA sequencing data and in situ miRNA detection, a novel fungal miRNA, Pmic_miR-8, was found to be transported into E. grandis roots after interaction with P. microcarpus Further characterization experiments demonstrate that inhibition of Pmic_miR-8 negatively impacts the maintenance of mycorrhizal roots in E. grandis, while supplementation of Pmic_miR-8 led to deeper integration of the fungus into plant tissues. Target prediction and experimental testing suggest that Pmic_miR-8 may target the host NB-ARC domain containing transcripts, suggesting a potential role for this miRNA in subverting host signaling to stabilize the symbiotic interaction. Altogether, we provide evidence of previously undescribed cross-kingdom sRNA transfer from ectomycorrhizal fungi to plant roots, shedding light onto the involvement of miRNAs during the developmental process of mutualistic symbioses.