Coniferous forest soils represent a globally important carbon sink, where the microbiome is essential for carbon flux between tree roots, rhizosphere, litter and soil. Soil habitats, such as roots, rhizosphere, bulk soil and litter differ in physicochemical properties and composition of highly specialized microbial communities, whose activity reflects the seasonality of temperature and tree activity of these mid- to high-latitude biomes. Here we present a multi-omic dataset encompassing 160 samples collected from four coniferous forest soil habitats in the Czech Republic and Norway, sampled in early summer, late summer, early winter and late winter that characterize the composition, genomic potential and activity of tree roots and microbiome. For each sample, we provide metabarcoding-based composition of bacterial, fungal and eukaryotic communities, results of shotgun DNA sequencing (metagenomes) and shotgun RNA sequencing (metatranscriptomes) illustrating the functional potential and activity within habitats. This dataset enables analyses of the temporal variation of taxonomic composition, functional potential and transcription across seasons in a temperate and boreal coniferous forest.

The legume family originated ca. 60-65 million years ago and soon diversified into at least six lineages (now extant subfamilies). The signal of whole genome duplications (WGD) is apparent in species sampled from all six subfamilies. The early diversification has posed difficulties for resolving the legume backbone structure and the timing of WGDs, especially in Caesalpinioideae where the diversification and WGD signals coincide. In this study, we report the genome sequences and annotations for Cercis canadensis (Cercidoideae) and Chamaecrista fasciculata (Caesalpinioideae) to help resolve the timings of WGDs relative to subfamily origins and the ancestral legume karyotype. Analyses of genome assemblies from four subfamilies within Fabaceae show that the last common ancestor of all legumes likely had seven chromosomes, with a genome structure similar to the extant Cercis genome. The retained karyotype structure, the lack of a WGD in the last 100+ Mya (Cercis and the lineage leading to it following the eudicot γ whole-genome triplication), and the unusually slow rates of nucleotide substitution and structural evolution in the Cercis genome underscore its utility as a genomic proxy for the last common ancestor of all legume species. Our analysis supports an allopolyploid origin of Caesalpinioideae, with progenitors from lineages along the backbone of the legume phylogeny. Rapid diversification and the inferred allopolyploid origin of Caesalpinioideae provide a partial explanation for the difficulty in resolving the backbone of the legume phylogeny and early Caesalpinioideae diversification.

Crassulacean acid metabolism (CAM) is a specialized photosynthetic pathway that enhances water-use efficiency by temporally separating nocturnal CO2 uptake from daytime decarboxylation and carbon fixation. To uncover the regulatory mechanisms coordinating these temporal dynamics, we generated high-resolution, 48 h time-course transcriptomes for the CAM model Kalanchoe fedtschenkoi under both 12 h/12 h light/dark (LD) cycles and continuous light (LL). A rhythmicity analysis revealed that diel light cues are the dominant driver of transcript oscillations: 16,810 genes (54.3% of annotated genes) exhibited rhythmic expression only under LD, whereas just 399 genes (1.3%) remained rhythmic under LL. A smaller set of 3009 genes (9.7%) oscillated in both conditions, indicating that the intrinsic circadian clock sustains rhythmicity for a limited subset of the transcriptome. A gene co-expression network analysis revealed extensive integration between circadian clock components, core CAM pathway enzymes, and stomatal regulators, defining regulatory modules that coordinate metabolic and physiological timing. Notably, key hub genes associated with post-translational and post-transcriptional regulation, including the E3 ubiquitin ligase HUB2 and several pentatricopeptide repeat (PPR) proteins, act as central nodes in CAM-associated networks. This discovery implicates epigenetic and organellar regulation as previously unrecognized critical tiers of control in CAM. Together, our results support a regulatory model in which CAM rhythmicity is governed by both external light/dark cues and the endogenous circadian clock through multi-level control spanning transcriptional and protein-level regulation. To support community exploration, we also provide an interactive eFP (electronic Fluorescent Pictograph) browser for visualizing time-resolved gene expression profiles.

Increasing drought threatens global agriculture, especially in water-limited regions. Sorghum bicolor , a drought-tolerant C4 grass, is a promising bioenergy crop for cultivation on marginal lands, yet its molecular drought responses remain poorly understood. To uncover single-cell-level transcriptional responses to drought, we performed single-nucleus RNA sequencing on mature sorghum leaves under well-watered and drought conditions. We identified major cell types and analyzed differential gene expression across mesophyll, bundle sheath, epidermal, vascular, and stomatal cells. Surprisingly, drought effects on transcriptomes exceeded differences due to cell identity, revealing a shared response across cell types. We leveraged this convergence to identify candidate regulators of drought-responsive gene expression. These findings advance our understanding of sorghum drought adaptation and offer new targets for engineering enhanced water use efficiency in bioenergy crops. ### Competing Interest Statement The authors have declared no competing interest.

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.

Temperature stress is one of the major limiting environmental factors that negatively impact global crop yields. Kalanchoë fedtschenkoi is an obligate crassulacean acid metabolism (CAM) plant species, exhibiting much higher water-use efficiency and tolerance to drought and heat stresses than C3 or C4 plant species. Previous studies on gene expression responses to low- or high-temperature stress have been focused on C3 and C4 plants. There is a lack of information about the regulation of gene expression by low and high temperatures in CAM plants. To address this knowledge gap, we performed transcriptome sequencing (RNA-Seq) of leaf and root tissues of K. fedtschenkoi under cold (8 °C), normal (25 °C), and heat (37 °C) conditions at dawn (i.e., 2 h before the light period) and dusk (i.e., 2 h before the dark period). Our analysis revealed differentially expressed genes (DEGs) under cold or heat treatment in comparison to normal conditions in leaf or root tissue at each of the two time points. In particular, DEGs exhibiting either the same or opposite direction of expression change (either up-regulated or down-regulated) under cold and heat treatments were identified. In addition, we analyzed gene co-expression modules regulated by cold or heat treatment, and we performed in-depth analyses of expression regulation by temperature stresses for selected gene categories, including CAM-related genes, genes encoding heat shock factors and heat shock proteins, circadian rhythm genes, and stomatal movement genes. Our study highlights both the common and distinct molecular strategies employed by CAM and C3/C4 plants in adapting to extreme temperatures, providing new insights into the molecular mechanisms underlying temperature stress responses in CAM species.

### Competing Interest Statement The authors have declared no competing interest.

Low-cost plant substrates, such as soybean hulls, are used for various industrial applications. Filamentous fungi are important producers of Carbohydrate Active enZymes (CAZymes) required for the degradation of these plant biomass substrates. CAZyme production is tightly regulated by several transcriptional activators and repressors. One such transcriptional activator is CLR-2/ClrB/ManR, which has been identified as a regulator of cellulase and mannanase production in several fungi. However, the regulatory network governing the expression of cellulase and mannanase encoding genes has been reported to differ between fungal species. Previous studies showed that Aspergillus niger ClrB is involved in the regulation of (hemi-)cellulose degradation, although its regulon has not yet been identified. To reveal its regulon, we cultivated an A. niger ΔclrB mutant and control strain on guar gum (a galactomannan-rich substrate) and soybean hulls (containing galactomannan, xylan, xyloglucan, pectin and cellulose) to identify the genes that are regulated by ClrB. Gene expression data and growth profiling showed that ClrB is indispensable for growth on cellulose and galactomannan and highly contributes to growth on xyloglucan in this fungus. Therefore, we show that A. niger ClrB is crucial for the utilization of guar gum and the agricultural substrate, soybean hulls. Moreover, we show that mannobiose is most likely the physiological inducer of ClrB in A. niger and not cellobiose, which is considered to be the inducer of N. crassa CLR-2 and A. nidulans ClrB.

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

Black yeasts are polyextremotolerant fungi that contain high amounts of melanin in their cell wall and maintain a primar yeast form. These fungi grow in xeric, nutrient depletes environments which implies that they require highly flexible metabolisms and have been suggested to contain the ability to form lichen-like mutualisms with nearby algae and bacteria. However, the exact ecological niche and interactions between these fungi and their surrounding community are not well understood. We have isolated 2 novel black yeasts from the genus Exophiala that were recovered from dryland biological soil crusts. Despite notable differences in colony and cellular morphology, both fungi appear to be members of the same species, which has been named Exophiala viscosa (i.e. E. viscosa JF 03-3 Goopy and E. viscosa JF 03-4F Slimy). A combination of whole genome sequencing, phenotypic experiments, and melanin regulation experiments have been performed on these isolates to fully characterize these fungi and help decipher their fundamental niche within the biological soil crust consortium. Our results reveal that E. viscosa is capable of utilizing a wide variety of carbon and nitrogen sources potentially derived from symbiotic microbes, can withstand many forms of abiotic stresses, and excretes melanin which can potentially provide ultraviolet resistance to the biological soil crust community. Besides the identification of a novel species within the genus Exophiala, our study also provides new insight into the regulation of melanin production in polyextremotolerant fungi.

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.

Currents are unique drivers of oceanic phylogeography and thus determine the distribution of marine coastal species, along with past glaciations and sea-level changes. Here we reconstruct the worldwide colonization history of eelgrass (Zostera marina L.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Two west-to-east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marina in the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present there for ~243 ky (thousand years). Mediterranean populations were founded ~44 kya, while extant distributions along western and eastern Atlantic shores were founded at the end of the Last Glacial Maximum (~19 kya), with at least one major refuge being the North Carolina region. The recent colonization and five- to sevenfold lower genomic diversity of the Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

Iron (Fe) is crucial for metabolic functions of living organisms. Plants access occluded Fe through interactions with rhizosphere microorganisms and symbionts. Yet, the interplay between Fe addition and plant-mycorrhizal interactions, especially the molecular mechanisms underlying mycorrhiza-assisted Fe processing in plants, remains largely unexplored.We conducted mesocosms in Pinus plants inoculated with different ectomycorrhizal fungi (EMF) Suillus species under conditions with and without Fe coatings. Meta-transcriptomic, biogeochemical, and X-ray fluorescence imaging analyses were applied to investigate early-stage mycorrhizal roots.While Fe addition promoted Pinus growth, it concurrently reduced mycorrhiza formation rate, symbiosis-related metabolites in plant roots, and aboveground plant carbon and macronutrient content. This suggested potential trade-offs between Fe-enhanced plant growth and symbiotic performance. However, the extent of this trade-off may depend on interactions between host plants and EMF species. Interestingly, dual EMF species were more effective at facilitating plant Fe uptake by inducing diverse Fe-related functions than single-EMF species. This subsequently triggered various Fe-dependent physiological and biochemical processes in Pinus roots, significantly contributing to Pinus growth. However, this resulted in a greater carbon allocation to roots, relatively reducing the aboveground plant carbon content.Our study offers critical insights into how EMF communities rebalance benefits of Fe-induced effects on symbiotic partners.

Fungi have evolved over millions of years and their species diversity is predicted to be the second largest on the earth. Fungi have cross-kingdom interactions with many organisms that have mutually shaped their evolutionary trajectories. Zygomycete fungi hold a pivotal position in the fungal tree of life and provide important perspectives on the early evolution of fungi from aquatic to terrestrial environments. Phylogenomic analyses have found that zygomycete fungi diversified into two separate clades, the Mucoromycota which are frequently associated with plants and Zoopagomycota that are commonly animal-associated fungi. Genetic elements that contributed to the fitness and divergence of these lineages may have been shaped by the varied interactions these fungi have had with plants, animals, bacteria, and other microbes. To investigate this, we performed comparative genomic analyses of the two clades of zygomycetes in the context of Kingdom Fungi, benefiting from our generation of a new collection of zygomycete genomes, including nine produced for this study. We identified lineage-specific genomic content that may contribute to the disparate biology observed in these zygomycetes. Our findings include the discovery of undescribed diversity in CotH, a Mucormycosis pathogenicity factor, which was found in a broad set of zygomycetes. Reconciliation analysis identified multiple duplication events and an expansion of CotH copies throughout the Mucoromycotina, Mortierellomycotina, Neocallimastigomycota, and Basidiobolus lineages. A kingdom-level phylogenomic analysis also identified new evolutionary relationships within the subphyla of Mucoromycota and Zoopagomycota, including supporting the sister-clade relationship between Glomeromycotina and Mortierellomycotina and the placement of Basidiobolus as sister to other Zoopagomycota lineages.

Peatlands are crucial sinks for atmospheric carbon but are critically threatened due to warming climates. Sphagnum (peat moss) species are keystone members of peatland communities where they actively engineer hyperacidic conditions, which improves their competitive advantage and accelerates ecosystem-level carbon sequestration. To dissect the molecular and physiological sources of this unique biology, we generated chromosome-scale genomes of two Sphagnum species: S. divinum and S. angustifolium. Sphagnum genomes show no gene colinearity with any other reference genome to date, demonstrating that Sphagnum represents an unsampled lineage of land plant evolution. The genomes also revealed an average recombination rate an order of magnitude higher than vascular land plants and short putative U/V sex chromosomes. These newly described sex chromosomes interact with autosomal loci that significantly impact growth across diverse pH conditions. This discovery demonstrates that the ability of Sphagnum to sequester carbon in acidic peat bogs is mediated by interactions between sex, autosomes and environment.

Crassulacean acid metabolism (CAM) plants exhibit elevated drought and heat tolerance compared to C3 and C4 plants through an inverted pattern of day/night stomatal closure and opening for CO2 assimilation. However, the molecular responses to water-deficit conditions remain unclear in obligate CAM species. In this study, we presented genome-wide transcription sequencing analysis using leaf samples of an obligate CAM species Kalanchoë fedtschenkoi under moderate and severe drought treatments at two-time points of dawn (2-h before the start of light period) and dusk (2-h before the dark period). Differentially expressed genes were identified in response to environmental drought stress and a whole genome wide co-expression network was created as well. We found that the expression of CAM-related genes was not regulated by drought stimuli in K. fedtschenkoi. Our comparative analysis revealed that CAM species (K. fedtschenkoi) and C3 species (Arabidopsis thaliana, Populus deltoides 'WV94') share some common transcriptional changes in genes involved in multiple biological processes in response to drought stress, including ABA signaling and biosynthesis of secondary metabolites.

The jasmonic acid (JA) signaling pathway plays an important role in the establishment of the ectomycorrhizal symbiosis (ECM) between Laccaria bicolor and poplar. We previously showed that the L. bicolor effector MiSSP7 induces the stabilization of the poplar JAZ6, a JA co-repressor protein that binds to Populus MYC2.1 and MYC2.2, orthologs of the Arabidopsis MYC2 transcription factor (TF), blocking their activity. Here we showed that both TFs play a central role in root colonization by L. bicolor mycelium, since their overexpression decreased the formation of the Hartig net, the hyphal network involved in symbiotic nutrient exchanges. By combining RNA sequencing and DNA Affinity Purification sequencing (DAP-seq) analysis, we identified a core set of JA-responsive genes directly activated by poplar MYC2.1 and MYC2.2, that need to be bypassed by the fungi to colonize root apoplastic spaces. These genes encode for other TFs, receptor-like kinases and many defense-related proteins, including terpene synthases (TPS). Monoterpenes produced by some of these TPS impact L. bicolor growth and ECM formation, suggesting a role for poplar root monoterpenes as negative regulators of in planta fungal growth and ECM symbiosis. Significance statement The ectomycorrhizal symbiosis is a predominant mutualistic plant-fungus interaction occurring in forests, sustaining tree health. Ectomycorrhizal fungi colonize the root intercellularly establishing the symbiotic interface required for bidirectional nutrients exchanges, the Hartig net. During root colonization, the fungus L. bicolor produces the effector protein MiSSP7 that binds to the jasmonate co-receptor PtJAZ6, maintaining the repression of MYC2-targeted genes. Here we showed that defensive genes are major targets of MYC2, suggesting that their strict control is required to allow fungal colonization, with special emphasis on the host root monoterpene synthesis. Future research will focus on how root terpene defenses mediate belowground mutualistic interactions and how they can be manipulated to engineer plants with enhanced disease resistance but stable mutualistic interactions.

Efficient utilization of agro-industrial waste, such as sugar beet pulp, is crucial for the bio-based economy. The fungus Aspergillus niger possesses a wide array of enzymes that degrade complex plant biomass substrates, and several regulators have been reported to play a role in their production. The role of the regulators GaaR, AraR, and RhaR in sugar beet pectin degradation has previously been reported. However, genetic regulation of the degradation of sugar beet pulp has not been assessed in detail. In this study, we generated a set of single and combinatorial deletion mutants targeting the pectinolytic regulators GaaR, AraR, RhaR, and GalX as well as the (hemi-)cellulolytic regulators XlnR and ClrB to address their relative contribution to the utilization of sugar beet pulp. We show that A. niger has a flexible regulatory network, adapting to the utilization of (hemi-)cellulose at early timepoints when pectin degradation is impaired.

Currents are unique drivers of oceanic phylogeography and so determine the distribution of marine coastal species, along with past glaciations and sea level changes. Here, we reconstruct the worldwide colonization history of eelgrass (Zostera marina L.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Multiple west to east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marina in the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present since ∼208 Kya (thousand years ago). Mediterranean populations were founded ∼53 Kya while extant distributions along western and eastern Atlantic shores coincide with the end of the Last Glacial Maximum (∼20 Kya). The recent colonization and 5-to 7-fold lower genomic diversity of Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.