A hierarchy of environmental covariates control the global biogeography of soil bacterial richness

Fragmented micro-growth habitats present opportunities for alternative competitive outcomes

Bacteria in nature often thrive in fragmented environments, like soil pores, plant roots or plant leaves, leading to smaller isolated habitats, shared with fewer species. This spatial fragmentation can significantly influence bacterial interactions, affecting overall community diversity. To investigate this, we contrast paired bacterial growth in tiny picoliter droplets (1-3 cells per 35 pL up to 3-8 cells per species in 268 pL) with larger, uniform liquid cultures (about 2 million cells per 140 µl). We test four interaction scenarios using different bacterial strains: substrate competition, substrate independence, growth inhibition, and cell killing. In fragmented environments, interaction outcomes are more variable and sometimes even reverse compared to larger uniform cultures. Both experiments and simulations show that these differences stem mostly from variation in initial cell population growth phenotypes and their sizes. These effects are most significant with the smallest starting cell populations and lessen as population size increases. Simulations suggest that slower-growing species might survive competition by increasing growth variability. Our findings reveal how microhabitat fragmentation promotes diverse bacterial interaction outcomes, contributing to greater species diversity under competitive conditions.

Groundwater is a hidden global keystone ecosystem

Groundwater is a vital ecosystem of the global water cycle, hosting unique biodiversity and providing essential services to societies. Despite being the largest unfrozen freshwater resource, in a period of depletion by extraction and pollution, groundwater environments have been repeatedly overlooked in global biodiversity conservation agendas. Disregarding the importance of groundwater as an ecosystem ignores its critical role in preserving surface biomes. To foster timely global conservation of groundwater, we propose elevating the concept of keystone species into the realm of ecosystems, claiming groundwater as a keystone ecosystem that influences the integrity of many dependent ecosystems. Our global analysis shows that over half of land surface areas (52.6%) has a medium-to-high interaction with groundwater, reaching up to 74.9% when deserts and high mountains are excluded. We postulate that the intrinsic transboundary features of groundwater are critical for shifting perspectives towards more holistic approaches in aquatic ecology and beyond. Furthermore, we propose eight key themes to develop a science-policy integrated groundwater conservation agenda. Given ecosystems above and below the ground intersect at many levels, considering groundwater as an essential component of planetary health is pivotal to reduce biodiversity loss and buffer against climate change.

Improved description of terrestrial habitat types by including microbial communities as indicators

Soils host diverse communities of microorganisms essential for ecosystem functions and soil health. Despite their importance, microorganisms are not covered by legislation protecting biodiversity or habitats, such as the Habitats Directive. Advances in molecular methods have caused breakthroughs in microbial community analysis, and recent studies have shown that parts of the communities are habitat-specific. If distinct microbial communities are present in the habitat types defined in the Habitats Directive, the Directive may be improved by including these communities. Thus, monitoring and reporting of biodiversity and conservation status of habitat types could be based not only on plant communities but also on microbial communities. In the present study, bacterial and plant communities were examined in six habitat types defined in the Habitats Directive by conducting botanical surveys and collecting soil samples for amplicon sequencing across 19 sites in Denmark. Furthermore, selected physico-chemical properties expected to differ between habitat types and explain variations in community composition of bacteria and vegetation were analysed (pH, electrical conductivity (EC), soil texture, soil water repellency, soil organic carbon content (OC), inorganic nitrogen, and in-situ water content (SWC)). Despite some variations within the same habitat type and overlaps between habitat types, habitat-specific communities were observed for both bacterial and plant communities, but no correlation was observed between the alpha diversity of vegetation and bacteria. PERMANOVA analysis was used to evaluate the variables best able to explain variation in the community composition of vegetation and bacteria. Habitat type alone could explain 46% and 47% of the variation in bacterial and plant communities, respectively. Excluding habitat type as a variable, the best model (pH, SWC, OC, fine silt, and Shannon's diversity index for vegetation) could explain 37% of the variation for bacteria. For vegetation, the best model (pH, EC, ammonium content and Shannon's diversity index for bacteria) could explain 25% of the variation. Based on these results, bacterial communities could be included in the Habitats Directive to improve the monitoring, as microorganisms are more sensitive to changes in the environment compared to vegetation, which the current monitoring is based on.

The chosen few-variations in common and rare soil bacteria across biomes

Soil bacterial communities are dominated by a few abundant species, while their richness is associated with rare species with largely unknown ecological roles and biogeography. Analyses of previously published soil bacterial community data using a novel classification of common and rare bacteria indicate that only 0.4% of bacterial species can be considered common and are prevalent across biomes. The remaining bacterial species designated as rare are endemic with low relative abundances. Observations coupled with mechanistic models highlight the central role of soil wetness in shaping bacterial rarity. An individual-based model reveals systematic shifts in community composition induced by low carbon inputs in drier soils that deprive common species of exhibiting physiological advantages relative to other species. We find that only a "chosen few" common species shape bacterial communities across biomes; however, their contributions are curtailed in resource-limited environments where a larger number of rare species constitutes the soil microbiome.

Beyond Taxonomic Identification: Integration of Ecological Responses to a Soil Bacterial 16S rRNA Gene Database

High-throughput sequencing 16S rRNA gene surveys have enabled new insights into the diversity of soil bacteria, and furthered understanding of the ecological drivers of abundances across landscapes. However, current analytical approaches are of limited use in formalizing syntheses of the ecological attributes of taxa discovered, because derived taxonomic units are typically unique to individual studies and sequence identification databases only characterize taxonomy. To address this, we used sequences obtained from a large nationwide soil survey (GB Countryside Survey, henceforth CS) to create a comprehensive soil specific 16S reference database, with coupled ecological information derived from survey metadata. Specifically, we modeled taxon responses to soil pH at the OTU level using hierarchical logistic regression (HOF) models, to provide information on both the shape of landscape scale pH-abundance responses, and pH optima (pH at which OTU abundance is maximal). We identify that most of the soil OTUs examined exhibited a non-flat relationship with soil pH. Further, the pH optima could not be generalized by broad taxonomy, highlighting the need for tools and databases synthesizing ecological traits at finer taxonomic resolution. We further demonstrate the utility of the database by testing against geographically dispersed query 16S datasets; evaluating efficacy by quantifying matches, and accuracy in predicting pH responses of query sequences from a separate large soil survey. We found that the CS database provided good coverage of dominant taxa; and that the taxa indicating soil pH in a query dataset corresponded with the pH classifications of top matches in the CS database. Furthermore we were able to predict query dataset community structure, using predicted abundances of dominant taxa based on query soil pH data and the HOF models of matched CS database taxa. The database with associated HOF model outputs is released as an online portal for querying single sequences of interest (https://shiny-apps.ceh.ac.uk/ID-TaxER/), and flat files are made available for use in bioinformatic pipelines. The further development of advanced informatics infrastructures incorporating modeled ecological attributes along with new functional genomic information will likely facilitate large scale exploration and prediction of soil microbial functional biodiversity under current and future environmental change scenarios.

Opening a next-generation black box: Ecological trends for hundreds of species-like taxa uncovered within a single bacterial >99% 16S rRNA operational taxonomic unit

Current knowledge on environmental distribution and taxon richness of free-living bacteria is mainly based on cultivation-independent investigations employing 16S rRNA gene sequencing methods. Yet, 16S rRNA genes are evolutionarily rather conserved, resulting in limited taxonomic and ecological resolutions provided by this marker. The faster evolving protein-encoding gene priB was used to reveal ecological patterns hidden within a single operational taxonomic unit (OTU) defined by >99% 16S rRNA sequence similarity. The studied subcluster PnecC of the genus Polynucleobacter represents a ubiquitous group of abundant freshwater bacteria with cosmopolitan distribution, which is very frequently detected by diversity surveys of freshwater systems. Based on genome taxonomy and a large set of genome sequences, a sequence similarity threshold for delineation of species-like taxa could be established. In total, 600 species-like taxa were detected in 99 freshwater habitats scattered across three regions representing a latitudinal range of 3,400 km (42°N to 71°N) and a pH gradient of 4.2 to 8.6. In addition to the unexpectedly high richness, the increased taxonomic resolution revealed structuring of Polynucleobacter communities by a couple of macroecological trends, which was previously only demonstrated for phylogenetically much broader groups of bacteria. An unexpected pattern was the almost complete compositional separation of Polynucleobacter communities of Ca²⁺ -rich and Ca²⁺ -poor habitats. This compositional pattern strongly resembled the vicariance of plant species on silicate and limestone soils. The new cultivation-independent approach presented opened a window to an incredible, previously unseen diversity, and enables investigations aiming on deeper understanding of how environmental conditions shape bacterial communities and drive evolution of free-living bacteria.

The Utility of Macroecological Rules for Microbial Biogeography

Macroecological rules have been developed for plants and animals that describe large-scale distributional patterns and attempt to explain the underlying physiological and ecological processes behind them. Similarly, microorganisms exhibit patterns in relative abundance, distribution, diversity, and traits across space and time, yet it remains unclear the extent to which microorganisms follow macroecological rules initially developed for macroorganisms. Additionally, the usefulness of these rules as a null hypothesis when surveying microorganisms has yet to be fully evaluated. With rapid advancements in sequencing technology, we have seen a recent increase in microbial studies that utilize macroecological frameworks. Here, we review and synthesize these macroecological microbial studies with two main objectives: (1) to determine to what extent macroecological rules explain the distribution of host-associated and free-living microorganisms, and (2) to understand which environmental factors and stochastic processes may explain these patterns among microbial clades (archaea, bacteria, fungi, and protists) and habitats (host-associated and free living; terrestrial and aquatic). Overall, 78% of microbial macroecology studies focused on free living, aquatic organisms. In addition, most studies examined macroecological rules at the community level with only 35% of studies surveying organismal patterns across space. At the community level microorganisms often tracked patterns of macroorganisms for island biogeography (74% confirm) but rarely followed Latitudinal Diversity Gradients (LDGs) of macroorganisms (only 32% confirm). However, when microorganisms and macroorganisms shared the same macroecological patterns, underlying environmental drivers (e.g., temperature) were the same. Because we found a lack of studies for many microbial groups and habitats, we conclude our review by outlining several outstanding questions and creating recommendations for future studies in microbial ecology.

Diversity of Dominant Soil Bacteria Increases with Warming Velocity at the Global Scale

Understanding global soil bacterial diversity is important because of its role in maintaining a healthy global ecosystem. Given the effects of environmental changes (e.g., warming and human impact) on the diversity of animals and plants, effects on soil bacterial diversity are expected; however, they have been poorly evaluated at the global scale to date. Thus, in this study, we focused on the dominant soil bacteria, which are likely critical drivers of key soil processes worldwide, and investigated the effects of warming velocity and human activities on their diversity. Using a global dataset of bacteria, we performed spatial analysis to evaluate the effects of warming velocity and human activities, while statistically controlling for the potentially confounding effects of current climate and geographic parameters with global climate and geographic data. We demonstrated that the diversity of the dominant soil bacteria was influenced globally, not only by the aridity index (dryness) and pH but also by warming velocity from the Last Glacial Maximum (21,000 years ago) to the present, showing significant increases. The increase in bacterial diversity with warming velocity was particularly significant in forests and grasslands. An effect of human activity was also observed, but it was secondary to warming velocity. These findings provide robust evidence and advance our understanding of the effects of environmental changes (particularly global warming) on soil bacterial diversity at the global scale.

Alpha-, beta-, and gamma-diversity of bacteria varies across habitats

Bacteria are essential parts of ecosystems and are the most diverse organisms on the planet. Yet, we still do not know which habitats support the highest diversity of bacteria across multiple scales. We analyzed alpha-, beta-, and gamma-diversity of bacterial assemblages using 11,680 samples compiled by the Earth Microbiome Project. We found that soils contained the highest bacterial richness within a single sample (alpha-diversity), but sediment assemblages displayed the highest gamma-diversity. Sediment, biofilms/mats, and inland water exhibited the most variation in community composition among geographic locations (beta-diversity). Within soils, agricultural lands, hot deserts, grasslands, and shrublands contained the highest richness, while forests, cold deserts, and tundra biomes consistently harbored fewer bacterial species. Surprisingly, agricultural soils encompassed similar levels of beta-diversity as other soil biomes. These patterns were robust to the alpha- and beta- diversity metrics used and the taxonomic binning approach. Overall, the results support the idea that spatial environmental heterogeneity is an important driver of bacterial diversity.

Image-based quantification of soil microbial dead zones induced by nitrogen fertilization

Microbial communities in agricultural soils underpin many ecosystem services including the maintenance of soil structure, food production, water purification and carbon storage. However, the impact of fertilization on the health of microbial communities is not well understood. This study investigates the spatial and temporal dynamics of nitrogen (N) transport away from a fertilizer granule with pore scale resolution. Specifically, we examined how soil structure and moisture content influence fertilizer derived N movement through the soil pore network and the subsequent impact of on soil microbial communities. We develop a mathematical model to describe N transport and reactions in soil at the pore-scale. Using X-ray Computed Tomography scans, we reconstructed a microscale description of a soil-pore geometry as a computational mesh. Solving two-phase water/air model produced pore-scale water distributions at 15, 30 and 70% water-filled pore volume. The N-speciation model considered ammonium (NH₄⁺), nitrate (NO₃^-) and dissolved organic N (DON), and included N immobilization, ammonification and nitrification processes, as well as diffusion in soil solution. We simulated the dissolution of a fertilizer pellet and a pore scale N cycle at three different water saturations. To aid interpretation of the model results, microbial activity at a range of N concentrations was measured. The model showed that the diffusion and concentration of N in water films is critically dependent upon soil moisture and N species. We predict that the maximum NH₄⁺ and NO₃^- concentrations in soil solution around the pellet under dry conditions are in the order of 1 × 10³ and 1 × 10⁴ mol m^-3 respectively, and under wet conditions 2 × 10² and 1 × 10³ mol m^-3, respectively. Supporting experimental evidence suggests that these concentrations would be sufficient to reduce microbial activity in the short-term in the zone immediately around the fertilizer pellet (ranging from 0.9 to 3.8 mm), causing a major loss of soil biological functioning. This model demonstrates the importance of pore-scale processes in regulating N movement and their interactions with the soil microbiome.

Soil bacterial diversity mediated by microscale aqueous-phase processes across biomes

Soil bacterial diversity varies across biomes with potential impacts on soil ecological functioning. Here, we incorporate key factors that affect soil bacterial abundance and diversity across spatial scales into a mechanistic modeling framework considering soil type, carbon inputs and climate towards predicting soil bacterial diversity. The soil aqueous-phase content and connectivity exert strong influence on bacterial diversity for each soil type and rainfall pattern. Biome-specific carbon inputs deduced from net primary productivity provide constraints on soil bacterial abundance independent from diversity. The proposed heuristic model captures observed global trends of bacterial diversity in good agreement with predictions by an individual-based mechanistic model. Bacterial diversity is highest at intermediate water contents where the aqueous phase forms numerous disconnected habitats and soil carrying capacity determines level of occupancy. The framework delineates global soil bacterial diversity hotspots; located mainly in climatic transition zones that are sensitive to potential climate and land use changes.