Unearthing the soil-borne microbiome of land plants.

Plant-soil biodiversity interactions are fundamental for the functioning of terrestrial ecosystems. Yet, the existence of a set of globally distributed topsoil microbial and small invertebrate organisms consistently associated with land plants (i.e., their consistent soil-borne microbiome), together with the environmental preferences and functional capabilities of these organisms, remains unknown. We conducted a standardized field survey under 150 species of land plants, including 58 species of bryophytes and 92 of vascular plants, across 124 locations from all continents. We found that, despite the immense biodiversity of soil organisms, the land plants evaluated only shared a small fraction (less than 1%) of all microbial and invertebrate taxa that were present across contrasting climatic and soil conditions and vegetation types. These consistent taxa were dominated by generalist decomposers and phagotrophs and their presence was positively correlated with the abundance of functional genes linked to mineralization. Finally, we showed that crossing environmental thresholds in aridity (aridity index of 0.65, i.e., the transition from mesic to dry ecosystems), soil pH (5.5; i.e., the transition from acidic to strongly acidic soils), and carbon (less than 2%, the lower limit of fertile soils) can result in drastic disruptions in the associations between land plants and soil organisms, with potential implications for the delivery of soil ecosystem processes under ongoing global environmental change.

Linking niche size and phylogenetic signals to predict future soil microbial relative abundances.

Bacteria provide ecosystem services (e.g., biogeochemical cycling) that regulate climate, purify water, and produce food and other commodities, yet their distribution and likely responses to change or intervention are difficult to predict. Using bacterial 16S rRNA gene surveys of 1,381 soil samples from the Biomes of Australian Soil Environment (BASE) dataset, we were able to model relative abundances of soil bacterial taxonomic groups and describe bacterial niche space and optima. Hold out sample validated hypothetical causal networks (structural equation models; SEM) were able to predict the relative abundances of bacterial taxa from environmental data and elucidate soil bacterial niche space. By using explanatory SEM properties as indicators of microbial traits, we successfully predicted soil bacterial response, and in turn potential ecosystem service response, to near-term expected changes in the Australian climate. The methods developed enable prediction of continental-scale changes in bacterial relative abundances, and demonstrate their utility in predicting changes in bacterial function and thereby ecosystem services. These capabilities will be strengthened in the future with growing genome-level data.

The global distribution and environmental drivers of the soil antibiotic resistome.

BACKGROUND: Little is known about the global distribution and environmental drivers of key microbial functional traits such as antibiotic resistance genes (ARGs). Soils are one of Earth's largest reservoirs of ARGs, which are integral for soil microbial competition, and have potential implications for plant and human health. Yet, their diversity and global patterns remain poorly described. Here, we analyzed 285 ARGs in soils from 1012 sites across all continents and created the first global atlas with the distributions of topsoil ARGs. RESULTS: We show that ARGs peaked in high latitude cold and boreal forests. Climatic seasonality and mobile genetic elements, associated with the transmission of antibiotic resistance, were also key drivers of their global distribution. Dominant ARGs were mainly related to multidrug resistance genes and efflux pump machineries. We further pinpointed the global hotspots of the diversity and proportions of soil ARGs. CONCLUSIONS: Together, our work provides the foundation for a better understanding of the ecology and global distribution of the environmental soil antibiotic resistome. Video Abstract.

Correlation between in vitro toxicity of pesticides and in vivo risk guidelines in support of complex operating site risk management: A meta-analysis.

In vitro cell systems can support hazard characterization and identify mechanisms involved in toxicity; however, using in vitro data for risk assessment still is challenging. As part of an effort to develop approaches for a complex operating site used for biocide packaging and distribution, we evaluated in vitro assays that could be used in a site management format. Across 66 studies, 108 pesticides were assessed on ten human-derived cell types at four endpoints. In vitro IC50s were compared to in vivo guidelines, NOEL/NOAELs, and ADIs using Spearman correlation and linear regression models. While human neuroblastoma cells (SH-SY5Y) were the most sensitive, HepG2 was the most used cell line in evaluating the toxicity of pesticides. Amongst the ten human cell lines, the IC50s derived from SH-SY5Y cells, using MTT-24 & 48 h (the most used assay) correlated (rho = 0.56-0.79; p < 0.05) with ADIs and NOEL/NOAELs. Although in vitro cell systems have some limitations, the correlation between in vitro data derived from SH-SY5Y cells and in vivo safety guidelines can provide site investigators with a tool to survey and prioritize areas and media of concern at complex operating sites impacted by pesticide mixtures.

Global hotspots for soil nature conservation.

Soils are the foundation of all terrestrial ecosystems1. However, unlike for plants and animals, a global assessment of hotspots for soil nature conservation is still lacking2. This hampers our ability to establish nature conservation priorities for the multiple dimensions that support the soil system: from soil biodiversity to ecosystem services. Here, to identify global hotspots for soil nature conservation, we performed a global field survey that includes observations of biodiversity (archaea, bacteria, fungi, protists and invertebrates) and functions (critical for six ecosystem services) in 615 composite samples of topsoil from a standardized survey in all continents. We found that each of the different ecological dimensions of soils-that is, species richness (alpha diversity, measured as amplicon sequence variants), community dissimilarity and ecosystem services-peaked in contrasting regions of the planet, and were associated with different environmental factors. Temperate ecosystems showed the highest species richness, whereas community dissimilarity peaked in the tropics, and colder high-latitudinal ecosystems were identified as hotspots of ecosystem services. These findings highlight the complexities that are involved in simultaneously protecting multiple ecological dimensions of soil. We further show that most of these hotspots are not adequately covered by protected areas (more than 70%), and are vulnerable in the context of several scenarios of global change. Our global estimation of priorities for soil nature conservation highlights the importance of accounting for the multidimensionality of soil biodiversity and ecosystem services to conserve soils for future generations.

Brassica napus Bacterial Assembly Processes Vary with Plant Compartment and Growth Stage but Not between Lines.

Holobiont bacterial community assembly processes are an essential element to understanding the plant microbiome. To elucidate these processes, leaf, root, and rhizosphere samples were collected from eight lines of Brassica napus in Saskatchewan over the course of 10 weeks. We then used ecological null modeling to disentangle the community assembly processes over the growing season in each plant part. The root was primarily dominated by stochastic community assembly processes, which is inconsistent with previous studies that suggest of a highly selective root environment. Leaf assembly processes were primarily stochastic as well. In contrast, the rhizosphere was a highly selective environment. The dominant rhizosphere selection process leads to more similar communities. Assembly processes in all plant compartments were dependent on plant growth stage with little line effect on community assembly. The foundations of assembly in the leaf were due to the harsh environment, leading to dominance of stochastic effects, whereas the stochastic effects in the root interior likely arise due to competitive exclusion or priority effects. Engineering canola microbiomes should occur during periods of strong selection assuming strong selection could promote beneficial bacteria. For example, engineering the microbiome to resist pathogens, which are typically aerially born, should focus on the flowering period, whereas microbiomes to enhance yield should likely be engineered postflowering as the rhizosphere is undergoing strong selection. IMPORTANCE In order to harness the microbiome for more sustainable crop production, we must first have a better understanding of microbial community assembly processes that occurring during plant development. This study examines the bacterial community assembly processes of the leaf, root, and rhizosphere of eight different lines of Brassica napus over the growing season. The influence of growth stage and B. napus line were examined in conjunction with the assembly processes. Understanding what influences the assembly processes of crops might allow for more targeted breeding efforts by working with the plant to manipulate the microbiome when it is undergoing the strongest selection pressure.

Positron-emitting radiotracers spatially resolve unexpected biogeochemical relationships linked with methane oxidation in Arctic soils.

Arctic soils are marked by cryoturbic features, which impact soil-atmosphere methane (CH4 ) dynamics vital to global climate regulation. Cryoturbic diapirism alters C/N chemistry within frost boils by introducing soluble organic carbon and nutrients, potentially influencing microbial CH4 oxidation. CH4 oxidation in soils, however, requires a spatio-temporal convergence of ecological factors to occur. Spatial delineation of microbial activity with respect to these key microbial and biogeochemical factors at relevant scales is experimentally challenging in inherently complex and heterogeneous natural soil matrices. This work aims to overcome this barrier by spatially linking microbial CH4 oxidation with C/N chemistry and metagenomic characteristics. This is achieved by using positron-emitting radiotracers to visualize millimeter-scale active CH4 uptake areas in Arctic soils with and without diapirism. X-ray absorption spectroscopic speciation of active and inactive areas shows CH4 uptake spatially associates with greater proportions of inorganic N in diapiric frost boils. Metagenomic analyses reveal Ralstonia pickettii associates with CH4 uptake across soils along with pertinent CH4 and inorganic N metabolism associated genes. This study highlights the critical relationship between CH4 and N cycles in Arctic soils, with potential implications for better understanding future climate. Furthermore, our experimental framework presents a novel, widely applicable strategy for unraveling ecological relationships underlying greenhouse gas dynamics under global change.

Soil Buffering Capacity Can Be Used To Optimize Biostimulation of Psychrotrophic Hydrocarbon Remediation.

Effective bioremediation of hydrocarbons requires innovative approaches to minimize phosphate precipitation in soils of different buffering capacities. Understanding the mechanisms underlying sustained stimulation of bacterial activity remains a key challenge for optimizing bioremediation-particularly in northern regions. Positron emission tomography (PET) can trace microbial activity within the naturally occurring soil structure of intact soils. Here, we use PET to test two hypotheses: (1) optimizing phosphate bioavailability in soil will outperform a generic biostimulatory solution in promoting hydrocarbon remediation and (2) oligotrophic biostimulation will be more effective than eutrophic approaches. In so doing, we highlight the key bacterial taxa that underlie aerobic and anaerobic hydrocarbon degradation in subarctic soils. In particular, we showed that (i) optimized phosphate bioavailability outperformed generic biostimulatory solutions in promoting hydrocarbon degradation, (ii) oligotrophic biostimulation is more effective than eutrophic approaches, and (iii) optimized biostimulatory solutions stimulated specific soil regions and bacterial consortia. The knowledge gleaned from this study will be crucial in developing field-scale biodegradation treatments for sustained stimulation of bacterial activity in northern regions.

From the Outside in: An Overview of Positron Imaging of Plant and Soil Processes.

Positron-emitting nuclides have long been used as imaging agents in medical science to spatially trace processes non-invasively, allowing for real-time molecular imaging using low tracer concentrations. This ability to non-destructively visualize processes in real time also makes positron imaging uniquely suitable for probing various processes in plants and porous environmental media, such as soils and sediments. Here, we provide an overview of historical and current applications of positron imaging in environmental research. We highlight plant physiological research, where positron imaging has been used extensively to image dynamics of macronutrients, signalling molecules, trace elements, and contaminant metals under various conditions and perturbations. We describe how positron imaging is used in porous soils and sediments to visualize transport, flow, and microbial metabolic processes. We also address the interface between positron imaging and other imaging approaches, and present accompanying chemical analysis of labelled compounds for reviewed topics, highlighting the bridge between positron imaging and complementary techniques across scales. Finally, we discuss possible future applications of positron imaging and its potential as a nexus of interdisciplinary biogeochemical research.

An intensive multilocation temporal dataset of fungal and bacterial communities in the root and rhizosphere of Brassica napus.

The plant microbiome has been recently recognized as a plant phenotype to help in the food security of the future population. However, global plant microbiome datasets are insufficient to be used effectively for breeding this new generation of crop plants. We surveyed the diversity and temporal composition of bacterial and fungal communities in the root and rhizosphere of Brassica napus, the world's second largest oilseed crop, weekly in eight diverse lines at one site and every three weeks in sixteen lines, at three sites in 2016 and 2017 in the Canadian Prairies. We sequenced the bacterial 16S ribosomal RNA gene generating a total of 127.7 million reads and the fungal internal transcribed spacer (ITS) region generating 113.4 million reads. 14,944 unique fungal amplicon sequence variants (ASV) were detected, with an average of 43 ASVs per root and 105 ASVs per rhizosphere sample. We detected 10,882 unique bacterial ASVs with an average of 249 ASVs per sample. Temporal, site-to-site, and line-driven variability were key determinants of microbial community structure. This dataset is a valuable resource to systematically extract information on the belowground microbiome of diverse B. napus lines in different environments, at different times in the growing season, in order to adapt effective varieties for sustainable crop production systems.

Is subarctic forest advance able to keep pace with climate change?

Recent climate warming and scenarios for further warming have led to expectations of rapid movement of ecological boundaries. Here we focus on the circumarctic forest-tundra ecotone (FTE), which represents an important bioclimatic zone with feedbacks from forest advance and corresponding tundra disappearance (up to 50% loss predicted this century) driving widespread ecological and climatic changes. We address FTE advance and climate history relations over the 20th century, using FTE response data from 151 sites across the circumarctic area and site-specific climate data. Specifically, we investigate spatial uniformity of FTE advance, statistical associations with 20th century climate trends, and whether advance rates match climate change velocities (CCVs). Study sites diverged into four regions (Eastern Canada; Central and Western Canada and Alaska; Siberia; and Western Eurasia) based on their climate history, although all were characterized by similar qualitative patterns of behaviour (with about half of the sites showing advancing behaviour). The main associations between climate trend variables and behaviour indicate the importance of precipitation rather than temperature for both qualitative and quantitative behaviours, and the importance of non-growing season as well as growing season months. Poleward latitudinal advance rates differed significantly among regions, being smallest in Eastern Canada (~10 m/year) and largest in Western Eurasia (~100 m/year). These rates were 1-2 orders of magnitude smaller than expected if vegetation distribution remained in equilibrium with climate. The many biotic and abiotic factors influencing FTE behaviour make poleward advance rates matching predicted 21st century CCVs (~103 -104  m/year) unlikely. The lack of empirical evidence for swift forest relocation and the discrepancy between CCV and FTE response contradict equilibrium model-based assumptions and warrant caution when assessing global-change-related biotic and abiotic implications, including land-atmosphere feedbacks and carbon sequestration.

An intensive multilocation temporal dataset of fungal communities in the root and rhizosphere of Brassica napus.

The plant microbiome has been recently recognized as a plant phenotype to help in the food security of the future population. However, global plant microbiome datasets are insufficient to be used effectively for breeding this new generation of crop plants. We surveyed the diversity and temporal composition of fungal communities in the root and rhizosphere of Brassica napus, the world's second largest oilseed crop, weekly in eight diverse lines at one site and every three weeks in sixteen lines, at three sites in 2016 and 2017 in the Canadian Prairies. 14,944 unique amplicon sequence variants (ASV) were detected based on the internal transcribed spacer region, with an average of 43 ASVs per root and 105 ASVs per rhizosphere sample. Temporal, site-to-site, and line-driven variability were key determinants of fungal community structure. This dataset is a valuable resource to systematically extract information on the belowground microbiome of diverse B. napus lines in different environments, at different times in the growing season, in order to adapt effective varieties for sustainable crop production systems.

The Charosphere Promotes Mineralization of 13C-Phenanthrene by Psychrotrophic Microorganisms in Greenland Soils.

When soil is frozen, biochar promotes petroleum hydrocarbon (PHC) degradation, yet we still do not understand why. To investigate microbial biodegradation activity under frozen conditions, we placed 60-µm mesh bags containing 6% (v/v) biochar created from fishmeal, bonemeal, bone chip, or wood into PHC-contaminated soil, which was then frozen to -5°C. This created three soil niches: biochar particles, the charosphere (biochar-contiguous soil), and bulk soil outside of the bags. After 90 d, C-phenanthrene mineralization reached 55% in bonemeal biochar and 84% in bone chip biochar charosphere soil, compared with only 43% in bulk soil and 13% in bone chip biochar particles. Soil pH remained near neutral in bone chip and bonemeal biochar treatments, unlike wood biochar, which increased alkalinity and likely made phosphate unavailable for microorganisms. Generally, charosphere soil had higher aromatic degradative gene abundances than bulk soil, but gene abundance was not directly linked to C-phenanthrene mineralization. In bone chip biochar-amended soils, phosphate successfully predicted microbial community composition, and abundances of and increased in charosphere soil. Biochar effects on charosphere soil were dependent on feedstock material and suggest that optimizing the charosphere in bone-derived biochars may increase remediation success in northern regions.

Structural equation modeling of a winnowed soil microbiome identifies how invasive plants re-structure microbial networks.

The development of microbial networks is central to ecosystem functioning and is the hallmark of complex natural systems. Characterizing network development over time and across environmental gradients is hindered by the millions of potential interactions among community members, limiting interpretations of network evolution. We developed a feature selection approach using data winnowing that identifies the most ecologically influential microorganisms within a network undergoing change. Using a combination of graph theory, leave-one-out analysis, and statistical inference, complex microbial communities are winnowed to identify the core organisms responding to external gradients or functionality, and then network development is evaluated against these externalities. In a plant invasion case study, the winnowed microbial network became more influential as the plant invasion progressed as a result of direct plant-microbe links rather than the expected indirect plant-soil-microbe links. This represents the first use of structural equation modeling to predict microbial network evolution, which requires identification of keystone taxa and quantification of the ecological processes underpinning community structure and function patterns.

Core and Differentially Abundant Bacterial Taxa in the Rhizosphere of Field Grown Brassica napus Genotypes: Implications for Canola Breeding.

Modifying the rhizosphere microbiome through targeted plant breeding is key to harnessing positive plant-microbial interrelationships in cropping agroecosystems. Here, we examine the composition of rhizosphere bacterial communities of diverse Brassica napus genotypes to identify: (1) taxa that preferentially associate with genotypes, (2) core bacterial microbiota associated with B. napus, (3) heritable alpha diversity measures at flowering and whole growing season, and (4) correlation between microbial and plant genetic distance among canola genotypes at different growth stages. Our aim is to identify and describe signature microbiota with potential positive benefits that could be integrated in B. napus breeding and management strategies. Rhizosphere soils of 16 diverse genotypes sampled weekly over a 10-week period at single location as well as at three time points at two additional locations were analyzed using 16S rRNA gene amplicon sequencing. The B. napus rhizosphere microbiome was characterized by diverse bacterial communities with 32 named bacterial phyla. The most abundant phyla were Proteobacteria, Actinobacteria, and Acidobacteria. Overall microbial and plant genetic distances were highly correlated (R = 0.65). Alpha diversity heritability estimates were between 0.16 and 0.41 when evaluated across growth stage and between 0.24 and 0.59 at flowering. Compared with a reference B. napus genotype, a total of 81 genera were significantly more abundant and 71 were significantly less abundant in at least one B. napus genotype out of the total 558 bacterial genera. Most differentially abundant genera were Proteobacteria and Actinobacteria followed by Bacteroidetes and Firmicutes. Here, we also show that B. napus genotypes select an overall core bacterial microbiome with growth-stage-related patterns as to how taxa joined the core membership. In addition, we report that sets of B. napus core taxa were consistent across our three sites and 2 years. Both differential abundance and core analysis implicate numerous bacteria that have been reported to have beneficial effects on plant growth including disease suppression, antifungal properties, and plant growth promotion. Using a multi-site year, temporally intensive field sampling approach, we showed that small plant genetic differences cause predictable changes in canola microbiome and are potential target for direct and indirect selection within breeding programs.

Biotic and abiotic drivers of tree seedling recruitment across an alpine treeline ecotone.

Treeline responses to climate change ultimately depend on successful seedling recruitment, which requires dispersal of viable seeds and establishment of individual propagules in novel environments. In this study, we evaluated the effects of several abiotic and biotic drivers of early tree seedling recruitment across an alpine treeline ecotone. In two consecutive years, we sowed seeds of low- and high-elevation provenances of Larix decidua (European larch) and Picea abies (Norway spruce) below, at, and above the current treeline into intact vegetation and into open microsites with artificially removed surface vegetation, as well as into plots protected from seed predators and herbivores. Seedling emergence and early establishment in treatment and in control plots were monitored over two years. Tree seedling emergence occurred at and several hundred metres above the current treeline when viable seeds and suitable microsites for germination were available. However, dense vegetation cover at lower elevations and winter mortality at higher elevations particularly limited early recruitment. Post-dispersal predation, species, and provenance also affected emergence and early establishment. This study demonstrates the importance of understanding multiple abiotic and biotic drivers of early seedling recruitment that should be incorporated into predictions of treeline dynamics under climate change.

Who Is the Rock Miner and Who Is the Hunter? The Use of Heavy-Oxygen Labeled Phosphate (P18O4) to Differentiate between C and P Fluxes in a Benzene-Degrading Consortium.

Phosphorus availability and cycling in microbial communities is a key determinant of bacterial activity. However, identifying organisms critical to P cycling in complex biodegrading consortia has proven elusive. Here we assess a new DNA stable isotope probing (SIP) technique using heavy oxygen-labeled phosphate (P18O4) and its effectiveness in pure cultures and a nitrate-reducing benzene-degrading consortium. First, we successfully labeled pure cultures of Gram-positive Micrococcus luteus and Gram-negative Bradyrhizobium elkanii and separated isotopically light and heavy DNA in pure cultures using centrifugal analyses. Second, using high-throughput amplicon sequencing of 16S rRNA genes to characterize active bacterial taxa (13C-labeled), we found taxa like Betaproteobacteria were key in denitrifying benzene degradation and that other degrading (nonhydrocarbon) inactive taxa (P18O4-labeled) like Staphylococcus and Corynebacterium may promote degradation through production of secondary metabolites (i.e., "helper" or "rock miner" bacteria). Overall, we successfully separated active and inactive taxa in contaminated soils, demonstrating the utility of P18O4-DNA SIP for identifying actively growing bacterial taxa. We also identified potential "miner" bacteria that choreograph hydrocarbon degradation by other microbes (i.e., the "hunters") without directly degrading contaminants themselves. Thus, while several taxa degrade benzene under denitrifying conditions, microbial benzene degradation may be enhanced by both direct degraders and miner bacteria.

Salix arctica changes root distribution and nutrient uptake in response to subsurface nutrients in High Arctic deserts.

Moisture is critical for plant success in polar deserts but not by the obvious pathway of reduced water stress. We hypothesized that an indirect, nutrient-linked, pathway resulting from unique water/frozen soil interactions in polar deserts creates nutrient-rich patches critical for plant growth. These nutrient-rich patches (diapirs) form deep in High Arctic polar deserts soils from water accumulating at the permafrost freezing front and ultimately rising into the upper soil horizons through cryoturbated convective landforms (frost boils). To determine if diapirs provide an enhanced source of plant-available N for Salix arctica (Arctic willow), we characterized soil, root, stem, and leaf 15 N natural abundance across 24 diapir and non-diapir frost boils in a High Arctic granitic semi-desert. When diapir horizons were available, S. arctica increased its subsurface (i.e., diapir) N uptake and plant root biomass doubled within diapir. Plant uptake of enriched 15 N injected into organic rich soil patches was 2.5-fold greater in diapir than in non-diapir frost boils. S. arctica percent cover was often higher (7.3 ± 1.0 [mean ± SE]) on diapiric frost boils, compared to frost boils without diapirs (4.4 ± 0.7), potentially reflecting the additional 20% nitrogen available in the subsurface of diapiric frost boils. Selective N acquisition from diapirs is a mechanism by which soil moisture indirectly enhances plant growth. Our work suggests that diapirs may be one mechanism contributing to Arctic greening by shrub expansion.

Archaea and bacteria mediate the effects of native species root loss on fungi during plant invasion.

Although invasive plants can drive ecosystem change, little is known about the directional nature of belowground interactions between invasive plants, native roots, bacteria, archaea and fungi. We used detailed bioinformatics and a recently developed root assay on soils collected in fescue grassland along a gradient of smooth brome (Bromus inermis Leyss) invasion to examine the links between smooth brome shoot litter and root, archaea, bacteria and fungal communities. We examined (1) aboveground versus belowground influences of smooth brome on soil microbial communities, (2) the importance of direct versus microbe-mediated impacts of plants on soil fungal communities, and (3) the web of roots, shoots, archaea, bacteria and fungi interactions across the A and B soil horizons in invaded and non-invaded sites. Archaea and bacteria influenced fungal composition, but not vice versa, as indicated by redundancy analyses. Co-inertia analyses suggested that bacterial-fungal variance was driven primarily by 12 bacterial operational taxonomic units (OTUs). Brome increased bacterial diversity via smooth brome litter in the A horizon and roots in the B horizon, which then reduced fungal diversity. Archaea increased abundance of several bacterial OTUs, and the key bacterial OTUs mediated changes in the fungi's response to invasion. Overall, native root diversity loss and bacterial mediation were more important drivers of fungal composition than were the direct effects of increases in smooth brome. Critically, native plant species displacement and root loss appeared to be the most important driver of fungal composition during invasion. This causal web likely gives rise to the plant-fungi feedbacks, which are an essential factor determining plant diversity in invaded grassland ecosystems.