Phylogenetic diversity of 200+ isolates of the ectomycorrhizal fungus Cenococcum geophilum associated with Populus trichocarpa soils in the Pacific Northwest, USA and comparison to globally distributed representatives.

The ectomycorrhizal fungal symbiont Cenococcum geophilum is of high interest as it is globally distributed, associates with many plant species, and has resistance to multiple environmental stressors. C. geophilum is only known from asexual states but is often considered a cryptic species complex, since extreme phylogenetic divergence is often observed within nearly morphologically identical strains. Alternatively, C. geophilum may represent a highly diverse single species, which would suggest cryptic but frequent recombination. Here we describe a new isolate collection of 229 C. geophilum isolates from soils under Populus trichocarpa at 123 collection sites spanning a ~283 mile north-south transect in Western Washington and Oregon, USA (PNW). To further understanding of the phylogenetic relationships within C. geophilum, we performed maximum likelihood and Bayesian phylogenetic analyses to assess divergence within the PNW isolate collection, as well as a global phylogenetic analysis of 789 isolates with publicly available data from the United States, Japan, and European countries. Phylogenetic analyses of the PNW isolates revealed three distinct phylogenetic groups, with 15 clades that strongly resolved at >80% bootstrap support based on a GAPDH phylogeny and one clade segregating strongly in two principle component analyses. The abundance and representation of PNW isolate clades varied greatly across the North-South range, including a monophyletic group of isolates that spanned nearly the entire gradient at ~250 miles. A direct comparison between the GAPDH and ITS rRNA gene region phylogenies, combined with additional analyses revealed stark incongruence between the ITS and GAPDH gene regions, consistent with intra-species recombination between PNW isolates. In the global isolate collection phylogeny, 34 clades were strongly resolved using Maximum Likelihood and Bayesian approaches (at >80% MLBS and >0.90 BPP respectively), with some clades having intra- and intercontinental distributions. Together these data are highly suggestive of divergence within multiple cryptic species, however additional analyses such as higher resolution genotype-by-sequencing approaches are needed to distinguish potential species boundaries and the mode and tempo of recombination patterns.

One-time nitrogen fertilization shifts switchgrass soil microbiomes within a context of larger spatial and temporal variation.

Soil microbiome responses to short-term nitrogen (N) inputs remain uncertain when compared with previous research that has focused on long-term fertilization responses. Here, we examined soil bacterial/archaeal and fungal communities pre- and post-N fertilization in an 8 year-old switchgrass field, in which twenty-four plots received N fertilization at three levels (0, 100, and 200 kg N ha-1 as NH4NO3) for the first time since planting. Soils were collected at two depths, 0-5 and 5-15 cm, for DNA extraction and amplicon sequencing of 16S rRNA genes and ITS regions for assessment of microbial community composition. Baseline assessments prior to fertilization revealed no significant pre-existing divergence in either bacterial/archaeal or fungal communities across plots. The one-time N fertilizations increased switchgrass yields and tissue N content, and the added N was nearly completely removed from the soil of fertilized plots by the end of the growing season. Both bacterial/archaeal and fungal communities showed large spatial (by depth) and temporal variation (by season) within each plot, accounting for 17 and 12-22% of the variation as calculated from the Sq. root of PERMANOVA tests for bacterial/archaeal and fungal community composition, respectively. While N fertilization effects accounted for only ~4% of overall variation, some specific microbial groups, including the bacterial genus Pseudonocardia and the fungal genus Archaeorhizomyces, were notably repressed by fertilization at 200 kg N ha-1. Bacterial groups varied with both depth in the soil profile and time of sampling, while temporal variability shaped the fungal community more significantly than vertical heterogeneity in the soil. These results suggest that short-term effects of N fertilization are significant but subtle, and other sources of variation will need to be carefully accounted for study designs including multiple intra-annual sampling dates, rather than one-time "snapshot" analyses that are common in the literature. Continued analyses of these trends over time with fertilization and management are needed to understand how these effects may persist or change over time.

Rhizosphere microbiomes diverge among Populus trichocarpa plant-host genotypes and chemotypes, but it depends on soil origin.

BACKGROUND: Plants have developed defense strategies for phytopathogen and herbivore protection via coordinated metabolic mechanisms. Low-molecular weight metabolites produced within plant tissues, such as salicylic acid, represent one such mechanism which likely mediates plant - microbe interactions above and below ground. Salicylic acid is a ubiquitous phytohormone at low levels in most plants, yet are concentrated defense compounds in Populus, likely acting as a selective filter for rhizosphere microbiomes. We propagated twelve Populus trichocarpa genotypes which varied an order of magnitude in salicylic acid (SA)-related secondary metabolites, in contrasting soils from two different origins. After four months of growth, plant properties (leaf growth, chlorophyll content, and net photosynthetic rate) and plant root metabolomics specifically targeting SA metabolites were measured via GC-MS. In addition, rhizosphere microbiome composition was measured via Illumina MiSeq sequencing of 16S and ITS2 rRNA-genes. RESULTS: Soil origin was the primary filter causing divergence in bacterial/archaeal and fungal communities with plant genotype secondarily influential. Both bacterial/archaeal and fungal evenness varied between soil origins and bacterial/archaeal diversity and evenness correlated with at least one SA metabolite (diversity: populin; evenness: total phenolics). The production of individual salicylic acid derivatives that varied by host genotype resulted in compositional differences for bacteria /archaea (tremuloidin) and fungi (salicylic acid) within one soil origin (Clatskanie) whereas soils from Corvallis did not illicit microbial compositional changes due to salicylic acid derivatives. Several dominant bacterial (e.g., Betaproteobacteria, Acidobacteria, Verrucomicrobia, Chloroflexi, Gemmatimonadete, Firmicutes) and one fungal phyla (Mortierellomycota) also correlated with specific SA secondary metabolites; bacterial phyla exhibited more negative interactions (declining abundance with increasing metabolite concentration) than positive interactions. CONCLUSIONS: These results indicate microbial communities diverge most among soil origin. However, within a soil origin, bacterial/archaeal communities are responsive to plant SA production within greenhouse-based rhizosphere microbiomes. Fungal microbiomes are impacted by root SA-metabolites, but overall to a lesser degree within this experimental context. These results suggest plant defense strategies, such as SA and its secondary metabolites, may partially drive patterns of both bacterial/archaeal and fungal taxa-specific colonization and assembly.