Palaeoclimate explains a unique proportion of the global variation in soil bacterial communities

Nitrogen deposition caused higher increases in plant-derived organic carbon than microbial-derived organic carbon in forest soils

Plant- and microbial-derived organic carbon, two components of the soil organic carbon (SOC) pool in terrestrial ecosystems, are regulated by increased atmospheric nitrogen (N) deposition. However, the spatial patterns and driving factors of the responses of plant- and microbial-derived SOC to N deposition in forests are not clear, which hinders our understanding of SOC sequestration. In this study, we explored the spatial patterns of plant- and microbial-derived SOC, and their responses to N addition and elucidated their underlying mechanisms in forest soils receiving N addition at four sites with various soil and climate conditions. Plant- and microbial-derived SOC were quantified using lignin phenols and amino sugars, respectively. N addition increased the total microbial residues by 20.5% on average ranging from 9.4% to 34.0% in temperate forests but not in tropical forests, and the increase was mainly derived from fungal residues. Lignin phenols increased more in temperate forests (average of 63.8%) than in tropical forests (average of 15.7%) following N addition. The ratio of total amino sugars to lignin phenols was higher in temperate forests than in tropical forests and decreased with N addition in temperate forests. N addition mainly regulated soil microbial residues by affecting pH, SOC, exchangeable Ca2+, gram-negative bacteria biomass, and the C:N ratio, while it mainly had indirect effects on lignin phenols by altering SOC, soil C:N ratio, and gram-negative bacteria biomass. Overall, our findings suggested that N deposition caused a greater increase in plant-derived SOC than in microbial-derived SOC and that plant-derived SOC would have a more important role in sequestering SOC under increasing N deposition in forest ecosystems, particularly in temperate forests.

Climate and local environment co-mediate the taxonomic and functional diversity of bacteria and archaea in the Qinghai-Tibet Plateau rivers

Understanding the environmental response patterns of riverine microbiota is essential for predicting the potential impact of future environmental change on river ecosystems. Vulnerable plateau ecosystems are particularly sensitive to climate and local environmental changes, however, the environmental response patterns of the taxonomic and functional diversity of riverine microbiota remain unclear. Here, we conducted a systematic investigation of the taxonomic and functional diversity of bacteria and archaea from riparian soils, sediments, and water across the elevation of 1800- 4800 m in the Qinghai-Tibet Plateau rivers. We found that within the elevation range of 1800 to 3800 m, riparian soils and sediments exhibited similarities and stabilities in microbial taxonomic and functional diversity, and water microbiomes were more sensitive with great fluctuations in microbial diversity. Beyond the elevation of 3800 m, microbial diversity declined across all riverine matrixes. Local environmental conditions can influence the sensitivity of microbiomes to climate change. The combination of critical climate and local environmental factors, including total nitrogen, total organic carbon, as well as climate variables associated with temperature and precipitation, provided better explanations for microbial diversity than single-factor analyses. Under the extremely adverse scenario of high greenhouse gas emission concentrations (SSP585), we anticipate that by the end of this century, the bacterial, archaeal, and microbial functional diversity across the river network of the Yangtze and Yellow source basin would potentially change by -16.9- 5.2 %, -16.1- 5.7 %, and -9.3- 6.4 %, respectively. Overall, climate and local environments jointly shaped the microbial diversity in plateau river ecosystems, and water microbiomes would provide early signs of environmental changes. Our study provides effective theoretical foundations for the conservation of river biodiversity and functional stability under environmental changes.

Changing soil available substrate primarily caused by fertilization management contributed more to soil respiration temperature sensitivity than microbial community thermal adaptation

Substrate depletion and microbial community thermal adaptation are major mechanisms that regulate the temperature sensitivity (Q10) of soil microbial respiration. Traditionally, the Q10 of soil microbial respiration is measured using laboratory incubation, which has limits in the continuous input of available substrates and the time scale for microbial community thermal adaptation. How the available substrate and the soil microbial community regulate the Q10 of soil microbial respiration under natural warming conditions remains unclear. To fill this gap in knowledge, a long-term field experiment was conducted consisting of two years of soil respiration observations combined with a soil available substrate and microbial community thermal adaptation analysis under seasonal warming conditions. The Q10 of soil respiration was calculated using the square root method, and it was more affected by the available substrate than by microbial community thermal adaptation. Fertilization management has a stronger effect on soil available substrate than temperature. As the temperature increased, NH4-N proved itself to be important for the bacterial community in the process of Q10 regulation, while dissolved organic carbon and nitrogen were key factors for the fungal community. Based on the niche breadth of microbial community composition, the changing Q10 of the soil respiration was not only closely associated with the specialist community, but also the generalist and neutralist communities. Furthermore, bacterial community thermal adaptation primarily occurred through shifts in the abundances of specialists and neutralists, while changes in species richness and species replacement occurred for the fungal generalists and neutralists. This work indicates that changing available nitrogen and DOC primarily caused by fertilization management contributed more in regulating the Q10 of soil microbial respiration than microbial community thermal adaptation, and there are different mechanisms for bacterial and fungal community thermal adaptation under warming.

Microbial impact on initial soil formation in arid and semiarid environments under simulated climate change

The microbiota is attributed to be important for initial soil formation under extreme climate conditions, but experimental evidence for its relevance is scarce. To fill this gap, we investigated the impact of in situ microbial communities and their interrelationship with biocrust and plants compared to abiotic controls on soil formation in initial arid and semiarid soils. Additionally, we assessed the response of bacterial communities to climate change. Topsoil and subsoil samples from arid and semiarid sites in the Chilean Coastal Cordillera were incubated for 16 weeks under diurnal temperature and moisture variations to simulate humid climate conditions as part of a climate change scenario. Our findings indicate that microorganism-plant interaction intensified aggregate formation and stabilized soil structure, facilitating initial soil formation. Interestingly, microorganisms alone or in conjunction with biocrust showed no discernible patterns compared to abiotic controls, potentially due to water-masking effects. Arid soils displayed reduced bacterial diversity and developed a new community structure dominated by Proteobacteria, Actinobacteriota, and Planctomycetota, while semiarid soils maintained a consistently dominant community of Acidobacteriota and Proteobacteria. This highlighted a sensitive and specialized bacterial community in arid soils, while semiarid soils exhibited a more complex and stable community. We conclude that microorganism-plant interaction has measurable impacts on initial soil formation in arid and semiarid regions on short time scales under climate change. Additionally, we propose that soil and climate legacies are decisive for the present soil microbial community structure and interactions, future soil development, and microbial responses.

Oligotrophic microbes are recruited to resist multiple global change factors in agricultural subsoils

An increasing number of anthropogenic pressures can have negative effects on biodiversity and ecosystem functioning. However, our understanding of how soil microbial communities and functions in response to multiple global change factors (GCFs) is still incomplete, particularly in less frequently disturbed subsoils. In this study, we examined the impact of different levels of GCFs (0-9) on soil functions and bacterial communities in both topsoils (0-20 cm) and subsoils (20-40 cm) of an agricultural ecosystem, and characterized the intrinsic factors influencing community resistance based on microbial life history strategy. Our experimental results showed a decline in soil multifunctionality, bacterial diversity, and community resistance as the number of GCFs increased, with a more drastic reduction in community resistance of subsoils. Specifically, we observed a significantly positive relationship between the oligotroph/copiotroph ratio and community resistance in subsoils, which was also verified by the negative correlation between 16S rRNA operon (rrn) copy number and community resistance. Structural equation modeling further revealed the direct effects of community resistance in promoting the ecosystem functioning, regardless of top- and subsoils. Therefore, these results suggested that subsoils may recruit more oligotrophic microbes to enhance their originally weaker community resistance under multiple GCFs, which was essential for maintaining sustainable agroecological functions and services. Overall, our study represents a significant advance in linking microbial life history strategy to the resistance of belowground microbial community and functionality.

The neglected roles of adjacent natural ecosystems in maintaining bacterial diversity in agroecosystems

A central aim of community ecology is to understand how local species diversity is shaped. Agricultural activities are reshaping and filtering soil biodiversity and communities; however, ecological processes that structure agricultural communities have often overlooked the role of the regional species pool, mainly owing to the lack of large datasets across several regions. Here, we conducted a soil survey of 941 plots of agricultural and adjacent natural ecosystems (e.g., forest, wetland, grassland, and desert) in 38 regions across diverse climatic and soil gradients to evaluate whether the regional species pool of soil microbes from adjacent natural ecosystems is important in shaping agricultural soil microbial diversity and completeness. Using a framework of multiscales community assembly, we revealed that the regional species pool was an important predictor of agricultural bacterial diversity and explained a unique variation that cannot be predicted by historical legacy, large-scale environmental factors, and local community assembly processes. Moreover, the species pool effects were associated with microbial dormancy potential, where taxa with higher dormancy potential exhibited stronger species pool effects. Bacterial diversity in regions with higher agricultural intensity was more influenced by species pool effects than that in regions with low intensity, indicating that the maintenance of agricultural biodiversity in high-intensity regions strongly depends on species present in the surrounding landscape. Models for community completeness indicated the positive effect of regional species pool, further implying the community unsaturation and increased potential in bacterial diversity of agricultural ecosystems. Overall, our study reveals the indubitable role of regional species pool from adjacent natural ecosystems in predicting bacterial diversity, which has useful implication for biodiversity management and conservation in agricultural systems.

Enzyme adaptation to habitat thermal legacy shapes the thermal plasticity of marine microbiomes

Microbial communities respond to temperature with physiological adaptation and compositional turnover. Whether thermal selection of enzymes explains marine microbiome plasticity in response to temperature remains unresolved. By quantifying the thermal behaviour of seven functionally-independent enzyme classes (esterase, extradiol dioxygenase, phosphatase, beta-galactosidase, nuclease, transaminase, and aldo-keto reductase) in native proteomes of marine sediment microbiomes from the Irish Sea to the southern Red Sea, we record a significant effect of the mean annual temperature (MAT) on enzyme response in all cases. Activity and stability profiles of 228 esterases and 5 extradiol dioxygenases from sediment and seawater across 70 locations worldwide validate this thermal pattern. Modelling the esterase phase transition temperature as a measure of structural flexibility confirms the observed relationship with MAT. Furthermore, when considering temperature variability in sites with non-significantly different MATs, the broadest range of enzyme thermal behaviour and the highest growth plasticity of the enriched heterotrophic bacteria occur in samples with the widest annual thermal variability. These results indicate that temperature-driven enzyme selection shapes microbiome thermal plasticity and that thermal variability finely tunes such processes and should be considered alongside MAT in forecasting microbial community thermal response.

Deciphering factors driving soil microbial life-history strategies in restored grasslands

In macroecology, the concept of r- and K-strategy has been widely applied, yet, there have been limited studies on microbial life-history strategies in temperate grasslands using multiple sequencing approaches. Total phospholipid fatty acid (PLFA) analysis, high-throughput meta-genomic sequencing, and GeoChip technologies were used to examine the changes in microbial life-history traits in a chronosequence of restored grasslands (1, 5, 10, 15, 25, and 30 years since restoration). Grassland restoration increased the relative abundances of Actinobacteria, Proteobacteria, and Bacteroidetes but reduced the relative abundances of Acidobacteria, Planctomycetes, and Chloroflexi. PLFA analysis revealed that grassland restoration reduced the fungi:bacteria and Gram-positive:Gram-negative bacteria ratios. Combined with the meta-genomic data, we found that grassland restoration shifted microorganisms from oligotrophic (K-) to copiotrophic (r-) groups, consistent with the increased rRNA operon copy number of the microbial community. Structural equation modeling showed that soil properties positively (p < 0.05) while plant properties negatively (p < 0.05) affected microbial life-history traits. We built a framework to highlight the importance of plant and soil properties in driving microbial life-history traits during grassland restoration. Finally, by incorporating meta-genomic and other microbiological data, this study showed that microbial life-history traits support the idea that rRNA operon copy number is a trait that reflects resource availability to soil microorganisms.

The neglected role of micronutrients in predicting soil microbial structure

Predicting the distribution patterns of soil microbial communities requires consideration of more environmental drivers. The effects of soil micronutrients on composition of microbial communities are largely unknown despite micronutrients closely relating to soil fertility and plant communities. Here we used data from 228 agricultural fields to identify the importance of micronutrients (iron, zinc, copper and manganese) in shaping structure of soil microbial communities (bacteria, fungi and protist) along latitudinal gradient over 3400 km, across diverse edaphic conditions and climatic gradients. We found that micronutrients explained more variations in the structure of microbial communities than macronutrients in maize soils. Moreover, micronutrients, particularly iron and copper, explained a unique percentage of the variation in structure of microbial communities in maize soils even after controlling for climate, soil physicochemical properties and macronutrients, but these effects were stronger for fungi and protist than for bacteria. The ability of micronutrients to predict the structure of soil microbial communities declined greatly in paddy soils. Machine learning approach showed that the addition of micronutrients substantially increased the predictive power by 9-17% in predicting the structure of soil microbial communities with up to 69-78% accuracy. These results highlighted the considerable contributions of soil micronutrients to microbial community structure, and advocated that soil micronutrients should be considered when predicting the structure of microbial communities in a changing world.

Climate legacies drive the distribution and future restoration potential of dryland forests

Knowing the extent and environmental drivers of forests is key to successfully restore degraded ecosystems, and to mitigate climate change and desertification impacts using tree planting. Water availability is the main limiting factor for the development of forests in drylands, yet the importance of groundwater resources and palaeoclimate as drivers of their current distribution has been neglected. Here we report that mid-Holocene climates and aquifer trends are key predictors of the distribution of dryland forests worldwide. We also updated the global extent of dryland forests to 1,283 million hectares and showed that failing to consider past climates and aquifers has resulted in ignoring or misplacing up to 130 million hectares of forests in drylands. Our findings highlight the importance of a wetter past and well-preserved aquifers to explain the current distribution of dryland forests, and can guide restoration actions by avoiding unsuitable areas for tree establishment in a drier world.

Continental-scale niche differentiation of dominant topsoil archaea in drylands

Archaea represent a diverse group of microorganisms often associated with extreme environments. However, an integrated understanding of biogeographical patterns of the specialist Haloarchaea and the potential generalist ammonia-oxidizing archaea (AOA) across large-scale environmental gradients remains limited. We hypothesize that niche differentiation determines their distinct distributions along environmental gradients. To test the hypothesis, we use a continental-scale research network including 173 dryland sites across northern China. Our results demonstrate that Haloarchaea and AOA dominate topsoil archaeal communities. As hypothesized, Haloarchaea and AOA show strong niche differentiation associated with two ecosystem types mainly found in China's drylands (i.e., deserts vs. grasslands), and they differ in the degree of habitat specialization. The relative abundance and richness of Haloarchaea are higher in deserts due to specialization to relatively high soil salinity and extreme climates, while those of AOA are greater in grassland soils. Our results further indicate a divergence in ecological processes underlying the segregated distributions of Haloarchaea and AOA. Haloarchaea are governed primarily by environmental-based processes while the more generalist AOA are assembled mostly via spatial-based processes. Our findings add to existing knowledge of large-scale biogeography of topsoil archaea, advancing our predictive understanding on changes in topsoil archaeal communities in a drier world. This article is protected by copyright. All rights reserved.

Microbial assemblies associated with temperature sensitivity of soil respiration along an altitudinal gradient

Identifying the drivers of the response of soil microbial respiration to warming is integral to accurately forecasting the carbon-climate feedbacks in terrestrial ecosystems. Microorganisms are the fundamental drivers of soil microbial respiration and its response to warming; however, the specific microbial communities and properties involved in the process remain largely undetermined. Here, we identified the associations between microbial community and temperature sensitivity (Q₁₀) of soil microbial respiration in alpine forests along an altitudinal gradient (from 2974 to 3558 m) from the climate-sensitive Tibetan Plateau. Our results showed that changes in microbial community composition accounted for more variations of Q₁₀ values than a wide range of other factors, including soil pH, moisture, substrate quantity and quality, microbial biomass, diversity and enzyme activities. Specifically, co-occurring microbial assemblies (i.e., ecological clusters or modules) targeting labile carbon consumption were negatively correlated with Q₁₀ of soil microbial respiration, whereas microbial assemblies associated with recalcitrant carbon decomposition were positively correlated with Q₁₀ of soil microbial respiration. Furthermore, there were progressive shifts of microbial assemblies from labile to recalcitrant carbon consumption along the altitudinal gradient, supporting relatively high Q₁₀ values in high-altitude regions. Our results provide new insights into the link between changes in major microbial assemblies with different trophic strategies and Q₁₀ of soil microbial respiration along an altitudinal gradient, highlighting that warming could have stronger effects on microbially-mediated soil organic matter decomposition in high-altitude regions than previously thought.

Relative Importance of Deterministic and Stochastic Processes on Soil Microbial Community Assembly in Temperate Grasslands

Changes in species composition across communities, i.e., β-diversity, is a central focus of ecology. Compared to macroorganisms, the β-diversity of soil microbes and its drivers are less studied. Whether the determinants of soil microbial β-diversity are consistent between soil depths and between abundant and rare microorganisms remains controversial. Here, using the 16S-rRNA of soil bacteria and archaea sampled at different soil depths (0–10 and 30–50 cm) from 32 sites along an aridity gradient of 1500 km in the temperate grasslands in northern China, we compared the effects of deterministic and stochastic processes on the taxonomic and phylogenetic β-diversity of soil microbes. Using variation partitioning and null models, we found that the taxonomic β-diversity of the overall bacterial communities was more strongly determined by deterministic processes in both soil layers (the explanatory power of environmental distance in topsoil: 25.4%; subsoil: 47.4%), while their phylogenetic counterpart was more strongly determined by stochastic processes (the explanatory power of spatial distance in topsoil: 42.1; subsoil 24.7%). However, in terms of abundance, both the taxonomic and phylogenetic β-diversity of the abundant bacteria in both soil layers was more strongly determined by deterministic processes, while those of rare bacteria were more strongly determined by stochastic processes. In comparison with bacteria, both the taxonomic and phylogenetic β-diversity of the overall abundant and rare archaea were strongly determined by deterministic processes. Among the variables representing deterministic processes, contemporary and historical climate and aboveground vegetation dominated the microbial β-diversity of the overall and abundant microbes of both domains in topsoils, but soil geochemistry dominated in subsoils. This study presents a comprehensive understanding on the β-diversity of soil microbial communities in the temperate grasslands in northern China. Our findings highlight the importance of soil depth, phylogenetic turnover, and species abundance in the assembly processes of soil microbial communities.

Soil properties and root traits jointly shape fine-scale spatial patterns of bacterial community and metabolic functions within a Korean pine forest

Spatial heterogeneity of soil bacterial community depends on scales. The fine-scale spatial heterogeneity of bacterial community composition and functions remains unknown. We analyzed the main driving factors of fine-scale spatial patterns of soil bacterial community composition and carbon metabolic functions across a 30 m × 40 m plot within a Korean pine forest by combining Illumina 16S rRNA sequencing with Biolog Ecoplates based on 53 soil samples. Clear spatial patterns in bacterial community composition and metabolic functions were observed in the forest soil. The bacterial community composition and metabolic functions both showed distance-decay of similarity within a distance of meters. Structural equation model analysis revealed that environmental variables and geographic distance together explained 37.9% and 63.1% of community and metabolic functions, respectively. Among all environmental factors, soil organic carbon (SOC) and root biomass emerged as the most important drivers of the bacterial community structure. In contrast, soil pH explained the largest variance in metabolic functions. Root biomass explained the second-largest variance in soil bacterial community composition, but root traits made no difference in metabolic functions variance. These results allow us to better understand the mechanisms controlling belowground diversity and plant-microbe interactions in forest ecosystems.

The biogeography of Streptomyces in New Zealand enabled by high-throughput sequencing of genus-specific rpoB amplicons

We evaluated Streptomyces biogeography in soils along a 1200 km latitudinal transect across New Zealand (NZ). Streptomyces diversity was examined using high-throughput sequencing of rpoB amplicons generated with a Streptomyces specific primer set. We detected 1287 Streptomyces rpoB operational taxonomic units (OTUs) with 159 ± 92 (average ± SD) rpoB OTUs per site. Only 12% (n = 149) of these OTUs matched rpoB sequences from cultured specimens (99% nucleotide identity cutoff). Streptomyces phylogenetic diversity (Faith's PD) was correlated with soil pH, mean annual temperature and plant community richness (Spearman's r: 0.77, 0.64 and -0.79, respectively; P < 0.05), but not with latitude. In addition, soil pH and plant community richness both explained significant variation in Streptomyces beta diversity. Streptomyces communities exhibited both high dissimilarity and strong dominance of one or a few species at each site. Taken together, these results suggest that dispersal limitation due to competitive interactions limits the colonization success of spores that relocate to new sites. Cultivated Streptomyces isolates represent a major source of clinically useful antibiotics, but only a small fraction of extant diversity within the genus have been identified and most species of Streptomyces have yet to be described.

Homogeneous selection drives antibiotic resistome in two adjacent sub-watersheds, China

Rivers are a significant reservoir of antibiotic resistance genes (ARGs), yet the biogeographic pattern of riverine ARGs and its underlying driving forces remain poorly understood. Here, we used metagenomic approach to investigate the spatio-temporal variation of ARGs in two adjacent sub-watersheds viz. North River (NR) and West River (WR), China. The results demonstrated that Bacitracin (22.8 % of the total ARGs), multidrug (20.7 %), sulfonamide (15.2 %) and tetracycline (10.9 %) were the dominant ARG types. SourceTracker analysis indicated that sewage treatment plants as the main source of ARGs, while animal feces mainly contributed in spreading the ARGs in the upstream of NR. Random forest and network analyses confirmed that NR was under the influence of fecal pollution. PCoA analysis demonstrated that the composition of ARGs changed along with the anthropogenic gradients, while the Raup-Crick null model showed that homogenizing selection mediated by class 1 integron intI1 resulted in stable ARG communities at whole watershed scale. Structural equation models revealed that microbial community, grassland and several non-antibiotic micropollutants may also play certain roles in influencing the distribution of ARGs. Overall, the observed deterministic formation of ARGs in riverine systems calls effective management strategies to mitigate the risks of antibiotic resistance on public health.

Precipitation controls seed bank size and its role in alpine meadow community regeneration with increasing altitude

The Tibetan Plateau has undergone significant climate warming in recent decades, and precipitation has also become increasingly variable. Much research has explored the effects of climate change on vegetation on this plateau. As potential vegetation buried in the soil, the soil seed bank is an important resource for ecosystem restoration and resilience. However, almost no studies have explored the effects of climate change on seed banks and the mechanisms of these effects. We used an altitudinal gradient to represent a decrease in temperature and collected soil seed bank samples from 27 alpine meadows (3,158-4,002 m) along this gradient. A structural equation model was used to explore the direct effects of mean annual precipitation (MAP) and mean annual temperature (MAT) on the soil seed bank and their indirect effects through aboveground vegetation and soil environmental factors. The species richness and abundance of the aboveground vegetation varied little along the altitudinal gradient, while the species richness and density of the seed bank decreased. The similarity between the seed bank and aboveground vegetation decreased with altitude; specifically, it decreased with MAP but was not related to MAT. The increase in MAP with increasing altitude directly decreased the species richness and density of the seed bank, while the increase in MAP and decrease in MAT with increasing altitude indirectly increased and decreased the species richness of the seed bank, respectively, by directly increasing and decreasing the species richness of the plant community. The size of the soil seed bank declined with increasing altitude. Increases in precipitation directly decreased the species richness and density and indirectly decreased the species richness of the seed bank with increasing elevation. The role of the seed bank in aboveground plant community regeneration decreases with increasing altitude, and this process is controlled by precipitation but not temperature.

Continental-scale metagenomics, BLAST searches, and herbarium specimens: The Australian Microbiome Initiative and the National Herbarium of Victoria

PREMISE

Motivated to make sensible interpretations of the massive volume of data from the Australian Microbiome Initiative (AusMic), we characterize the soil mycota of Australia. We establish operational taxonomic units (OTUs) from the data and compare these to GenBank and a data set from the National Herbarium of Victoria (MEL), Melbourne, Australia. We also provide visualizations of Agaricomycete diversity, drawn from our analyses of the AusMic sequences and taxonomy.

METHODS

The AusMic internal transcribed spacer (ITS) data were filtered to create OTUs, which were searched against the National Center for Biotechnology Information Nucleotide database and the MEL database. We further characterized a portion of our OTUs by graphing the counts of the families and orders of Agaricomycetes. We also graphed AusMic species determinations for Australian Agaricomycetes against latitude.

RESULTS

Our filtering process generated 192,325 OTUs; for Agaricomycetes, there were 27,730 OTUs. Based on the existing AusMic taxonomy at species level, we inferred the diversity of Australian Agaricomycetes against latitude to be lowest between -20 and -25 decimal degrees.

DISCUSSION

BLAST comparisons provided reciprocal insights between the three data sets, including the detection of unusual root-associated species in the AusMic data, insights into mushroom morphology from the MEL data, and points of comparison for the taxonomic determinations between AusMic, GenBank, and MEL. This study provides a tabulation of Australian fungi, different visual snapshots of a subset of those taxa, and a springboard for future studies.

Blind spots in global soil biodiversity and ecosystem function research

Soils harbor a substantial fraction of the world's biodiversity, contributing to many crucial ecosystem functions. It is thus essential to identify general macroecological patterns related to the distribution and functioning of soil organisms to support their conservation and consideration by governance. These macroecological analyses need to represent the diversity of environmental conditions that can be found worldwide. Here we identify and characterize existing environmental gaps in soil taxa and ecosystem functioning data across soil macroecological studies and 17,186 sampling sites across the globe. These data gaps include important spatial, environmental, taxonomic, and functional gaps, and an almost complete absence of temporally explicit data. We also identify the limitations of soil macroecological studies to explore general patterns in soil biodiversity-ecosystem functioning relationships, with only 0.3% of all sampling sites having both information about biodiversity and function, although with different taxonomic groups and functions at each site. Based on this information, we provide clear priorities to support and expand soil macroecological research.

Effects of historical legacies on soil nematode communities are mediated by contemporary environmental conditions

Both contemporary and historical factors are documented to be crucial in regulating species diversity and distribution. Soil fauna contribute substantially to global biodiversity and ecosystem functioning, while it is unclear whether and to what extent historical factors shape their diversity patterns. Here, we used soil nematodes as a model organism to test historical effects on soil fauna and to investigate the relative importance of climatic, soil, and historical factors. Based on nematode distribution data in 16 natural sites at a large scale (ranging from 22 to 40°N) in mainland China, we conducted elastic net regression model to test the effects of climatic (e.g., mean and seasonality of temperature/precipitation), soil (e.g., soil carbon, nitrogen, and pH), and historical (e.g., temperature/precipitation anomaly and the velocity of the change since the Last Glacial Maximum) variables on nematode genus richness and Shannon's diversity. Additionally, variation partitioning was used to determine the contribution of the three predictor sets to the explanation of both Jaccard and Bray-Curtis community dissimilarity. We found that climate generally explained more variations in both diversity and composition than soil and historical predictors in our samples. We also showed that although historical factors (e.g., temperature change velocity) were correlated with nematode diversity and composition, the pure effects of these historical factors were negligible. In other words, the historical effects were commonly represented by their interactions with current climatic and soil factors within our selected sites. Our results indicated that contemporary factors, especially climate, may outperform historical factors in regulating soil nematode diversity patterns at large scales.

Unique microbial module regulates the harmful algal bloom (Cochlodinium polykrikoides) and shifts the microbial community along the Southern Coast of Korea

Harmful algal blooms (HABs) of Cochlodinium (aka Margalefidinium) polykrikoides cause huge economic and ecological damages and thus are considered environmental problems. Previous studies uncovered that the formation and collapse of phytoplankton blooms could be closely related to their associated microbes although their roles in C. polykrikoides bloom have not been elucidated yet. To explore the potential interactions between C. polykrikoides and other microbes (archaea, bacteria, and phytoplankton), we collected water samples in the free-living (FL) (0.22 to 3 μm), nanoparticle-associated (NP) (3 to 20 μm), and microparticle-associated (MP) (>20 μm) fractions when C. polykrikoides blooms occurred from July to August in 2016, 2017, and 2018 in the South Sea of Korea. The microbial composition of the C. polykrikoides-associated microbial cluster (Module I) significantly differed from those of other modules associated with Alexandrium, Chaetoceros or Chattonella. Over half of the interspecies interactions in Module I occurred within the module. That is, specific microbial clusters were associated with the C. polykrikoides bloom. Structural equation modeling (SEM) further confirmed the stronger effects of Module I on C. polykrikoides blooms compared to environmental factors. Among the operational taxonomic units (OTUs) directly correlated with C. polykrikoides, Marine Group I was presumed to supply vitamin B₁₂, the essential element for C. polykrikoides growth, while the potential fish pathogens (Micrococcaceae and Piscirickettsiaceae) could contribute to the massive fish death together with C. polykrikoides itself. In addition, the zoospores of Syndiniales, a parasitoid to dinoflagellates, might be related to the sudden collapse of C. polykrikoides blooms. These microbial groups also contributed to significant alterations of the local microbial community structures. Collectively, network analysis and SEM revealed that the C. polykrikoides bloom is concomitant with distinct microbial communities, contributing to the rise and fall of the bloom, and finally determining the local microbial community structures.

Soil Microbial Biogeography in a Changing World: Recent Advances and Future Perspectives

Soil microbial communities are fundamental to maintaining key soil processes associated with litter decomposition, nutrient cycling, and plant productivity and are thus integral to human well-being. Recent technological advances have exponentially increased our knowledge concerning the global ecological distributions of microbial communities across space and time and have provided evidence for their contribution to ecosystem functions. However, major knowledge gaps in soil biogeography remain to be addressed over the coming years as technology and research questions continue to evolve. In this minireview, we state recent advances and future directions in the study of soil microbial biogeography and discuss the need for a clearer concept of microbial species, projections of soil microbial distributions toward future global change scenarios, and the importance of embracing culture and isolation approaches to determine microbial functional profiles. This knowledge will be critical to better predict ecosystem functions in a changing world.

Global ecosystem thresholds driven by aridity

Aridity, which is increasing worldwide because of climate change, affects the structure and functioning of dryland ecosystems. Whether aridification leads to gradual (versus abrupt) and systemic (versus specific) ecosystem changes is largely unknown. We investigated how 20 structural and functional ecosystem attributes respond to aridity in global drylands. Aridification led to systemic and abrupt changes in multiple ecosystem attributes. These changes occurred sequentially in three phases characterized by abrupt decays in plant productivity, soil fertility, and plant cover and richness at aridity values of 0.54, 0.7, and 0.8, respectively. More than 20% of the terrestrial surface will cross one or several of these thresholds by 2100, which calls for immediate actions to minimize the negative impacts of aridification on essential ecosystem services for the more than 2 billion people living in drylands.

Contrasting Biogeographic Patterns of Bacterial and Archaeal Diversity in the Top- and Subsoils of Temperate Grasslands

Biogeographic patterns and drivers of soil microbial diversity have been extensively studied in the past few decades. However, most research has focused on the topsoil, while the subsoil is assumed to have microbial diversity patterns similar to those of the topsoil. Here we compared patterns and drivers of microbial alpha and beta diversity in and between topsoils (0 to 10 cm) and subsoils (30 to 50 cm) of temperate grasslands in Inner Mongolia of China, covering an ∼1,500-km transect along an aridity gradient. Counter to the conventional assumption, we find contrasting biogeographic patterns of diversity and influencing factors for different bacterial and archaeal groups and between depths. While bacterial diversity remains constant or increases with increasing aridity in topsoil and decreases in subsoil, archaeal diversity decreases in topsoil and remains constant in subsoil. Microbial diversity in the topsoil is most strongly influenced by aboveground vegetation and contemporary climate but is most strongly influenced by the factor historical temperature anomaly since the Last Glacial Maximum (LGM) and by soil pH in the subsoil. Moreover, the biogeographic patterns of topsoil-subsoil community dissimilarities vary for different microbial groups and are overall most strongly influenced by soil fertility differences between depths for bacteria and by contemporary climate for archaea. These findings suggest that diversity patterns observed in the topsoil may not be readily applied to the subsoil horizons. For the subsoil in particular, historical climate plays a vital role in the spatial variation of bacterial diversity. Overall, our study provides novel information for understanding and predicting soil microbial diversity patterns at depth.IMPORTANCE Exploring the biogeographic patterns of soil microbial diversity is critical for understanding mechanisms underlying the response of soil processes to climate change. Using top- and subsoils from an ∼1,500-km temperate grassland transect, we find divergent patterns of microbial diversity and its determinants in the topsoil versus the subsoil. Furthermore, we find important and direct legacy effects of historical climate change on the microbial diversity of subsoil yet indirect effects on topsoil. Our findings challenge the conventional assumption of similar geographic patterns of soil microbial diversity along soil profiles and help to improve our understanding of how soil microbial communities may respond to future climate change in different regions with various climate histories.

Intra-horizon differentiation of the bacterial community and its co-occurrence network in a typical Plinthic horizon

Subsurface soil bacterial community composition and the controlling factors remain largely unknown, especially the micro-zone differentiation of community composition within a horizon. We studied a plinthic horizon to determine how different micro-zones in a horizon affect the bacterial community. The plinthic horizon is a net-like horizon characterized by the segregation of iron forms as shown by contrasting red matrix and white veins, which share common macro-environmental conditions such as climate and land use but differ only in physical and chemical compositions. The studied horizon is typical of the red soils of southeastern China and is an important layer in the red soil Critical Zone. The plinthite is considered to have been formed in the Quaternary and thus is a record of the paleo-environment. We evaluated the difference in the bacterial community composition between the red matrix and white veins and explored the possible assembly mechanisms of their co-occurrence patterns. Compared to the eutrophic environments of a red matrix, higher relative abundances of Acidobacteria and Nitrospirae were observed in the white veins. Similarly, more niches led to a higher density of bacterial co-occurrence patterns in the red matrix. The differences in the bacterial community composition and association networks are due to environmental selection, including the legacy of the paleoclimate that is represented by major element contents and contemporary hydrological properties that are mainly controlled by the soil texture. Our study shows that micro-zones even within a same plinthic horizon can provide different habitats and thus select for specific bacterial communities. Furthermore, this study could improve our understanding of the differentiation of bacterial communities among microenvironments caused by both historical and contemporary processes and help to predict how these communities may respond to future environmental changes.

Multifunctionality debt in global drylands linked to past biome and climate

Past vegetation and climatic conditions are known to influence current biodiversity patterns. However, whether their legacy effects affect the provision of multiple ecosystem functions, that is, multifunctionality, remains largely unknown. Here we analyzed soil nutrient stocks and their transformation rates in 236 drylands from six continents to evaluate the associations between current levels of multifunctionality and legacy effects of the Last Glacial Maximum (LGM) desert biome distribution and climate. We found that past desert distribution and temperature legacy, defined as increasing temperature from LGM, were negatively correlated with contemporary multifunctionality even after accounting for predictors such as current climate, soil texture, plant species richness, and site topography. Ecosystems that have been deserts since the LGM had up to 30% lower contemporary multifunctionality compared with those that were nondeserts during the LGM. In addition, ecosystems that experienced higher warming rates since the LGM had lower contemporary multifunctionality than those suffering lower warming rates, with a ~9% reduction per extra degree Celsius. Past desert distribution and temperature legacies had direct negative effects, while temperature legacy also had indirect (via soil sand content) negative effects on multifunctionality. Our results indicate that past biome and climatic conditions have left a strong "functionality debt" in global drylands. They also suggest that ongoing warming and expansion of desert areas may leave a strong fingerprint in the future functioning of dryland ecosystems worldwide that needs to be considered when establishing management actions aiming to combat land degradation and desertification.

Late Quaternary climate change explains soil fungal community composition rather than fungal richness in forest ecosystems

The dramatic climate fluctuations of the late Quaternary have influenced the diversity and composition of macroorganism communities, but how they structure belowground microbial communities is less well known. Fungi constitute an important component of soil microorganism communities. They play an important role in biodiversity maintenance, community assembly, and ecosystem functioning, and differ from many macroorganisms in many traits. Here, we examined soil fungal communities in Chinese temperate, subtropical, and tropic forests using Illumina MiSeq sequencing of the fungal ITS1 region. The relative effect of late Quaternary climate change and contemporary environment (plant, soil, current climate, and geographic distance) on the soil fungal community was analyzed. The richness of the total fungal community, along with saprotrophic, ectomycorrhizal (EM), and pathogenic fungal communities, was influenced primarily by the contemporary environment (plant and/or soil) but not by late Quaternary climate change. Late Quaternary climate change acted in concert with the contemporary environment to shape total, saprotrophic, EM, and pathogenic fungal community compositions and with a stronger effect in temperate forest than in tropic-subtropical forest ecosystems. Some contemporary environmental factors influencing total, saprotrophic, EM, and pathogenic fungal communities in temperate and tropic-subtropical forests were different. We demonstrate that late Quaternary climate change can help to explain current soil fungal community composition and argue that climatic legacies can help to predict soil fungal responses to climate change.

Drought and its legacy modulate the post-fire recovery of soil functionality and microbial community structure in a Mediterranean shrubland

The effects of drought on soil dynamics after fire are poorly known, particularly its long-term (i.e., years) legacy effects once rainfall returns to normal. Understanding this is particularly important for nutrient-poor soils in semi-arid regions affected by fire, in which rainfall is projected to decrease with climate change. Here, we studied the effects of post-fire drought and its legacy on soil microbial community structure and functionality in a Cistus-Erica shrubland (Spain). Rainfall total and patterns were experimentally modified to produce an unburned control (natural rainfall) and four burned treatments: control (natural rainfall), historical control (long-term average rainfall), moderate drought (percentile 8 historical rainfall, 5 months of drought per year), and severe drought (percentile 2, 7 months of drought). Soil nutrients and microbial community composition (ester-linked fatty acid approach) and functionality (enzyme activities and C mineralization rate) were monitored during the first 4 years after fire under rainfall treatments, plus two additional ones without them (six post-fire years). We found that the recovery of burned soils was lower under drought. Post-fire drought increased nitrate in the short term and reduced available phosphorus, exchangeable potassium, soil organic matter, enzyme activities, and carbon mineralization rate. Moreover, drought decreased soil total microbial biomass and fungi, with bacteria becoming relatively more abundant. Two years after discontinuing the drought treatments, the drought legacy was significant for available phosphorus and enzyme activities. Although microbial biomass did not show any drought legacy effect, the proportion of fungi and bacteria (mainly gram-positive) did, being lower and higher, respectively, in former drought-treated plots. We show that drought has an important impact on soil processes, and that some of its effects persist for at least 2 years after the drought ended. Therefore, drought and its legacy effects can be important for modeling biogeochemical processes in burned soils under future climate change.

The influence of climatic legacies on the distribution of dryland biocrust communities

Predicting the distribution of biocrust species, mosses, lichens and liverworts associated with surface soils is difficult, but climatic legacies (changes in climate over the last 20 k years) can improve our prediction of the distribution of biocrust species. To provide empirical support for this hypothesis, we used a combination of network analyses and structural equation modelling to identify the role of climatic legacies in predicting the distribution of ecological clusters formed by species of mosses, lichens and liverworts using data from 282 large sites distributed across 0.6 million km² of eastern Australia. Two ecological clusters contained 87% of the 120 moss, lichen and liverwort species. Both clusters contained lichen, moss and liverwort species, but were dominated by different families. Sites where the air temperature increased the most over 20k years (positive temperature legacies) were associated with reductions in the relative abundance of species from the lichen (Peltulaceae and Teloschistaceae) and moss (Bryaceae) families (Cluster A species), greater groundstorey plant cover and lower soil pH. Sites where precipitation has increased over the past 20k years (positive precipitation legacy) were associated with increases in the relative abundance of lichen (Cladoniaceae, Lecideaceae and Thelotremataceae) and moss (Pottiaceae) families (Cluster B species) and lower levels of soil pH. Sites where temperatures have increased the most in the past 20k years suppressed the negative effects of plant cover on Cluster B by reducing plant cover. Increased intensity of grazing suppressed the negative effect of soil pH and the positive effect of soil carbon, on the relative abundance of Cluster B taxa. Finally, increasing temperature and precipitation legacies reduced the negative effect of soil pH on Cluster B. Understanding of the importance of climatic legacies improves our ability to predict how biocrust assemblies might respond to ongoing global environmental change associated with increasing land use intensification, increasing temperature and reduced rainfall.

Existing Climate Change Will Lead to Pronounced Shifts in the Diversity of Soil Prokaryotes

Soil bacteria are key to ecosystem function and maintenance of soil fertility. Leveraging associations of current geographic distributions of bacteria with historic climate, we predict that soil bacterial diversity will increase across the majority (∼75%) of the Tibetan Plateau and northern North America if bacterial communities equilibrate with existing climatic conditions. This prediction is possible because the current distributions of soil bacteria have stronger correlations with climate from ∼50 years ago than with current climate. This lag is likely associated with the time it takes for soil properties to adjust to changes in climate. The predicted changes are location specific and differ across bacterial taxa, including some bacteria that are predicted to have reductions in their distributions. These findings illuminate the widespread potential of climate change to influence belowground diversity and the importance of considering bacterial communities when assessing climate impacts on terrestrial ecosystems. IMPORTANCE There have been many studies highlighting how plant and animal communities lag behind climate change, causing extinction and diversity debts that will slowly be paid as communities equilibrate. By virtue of their short generation times and dispersal abilities, soil bacteria might be expected to respond to climate change quickly and to be effectively in equilibrium with current climatic conditions. We found strong evidence to the contrary in Tibet and North America. These findings could significantly improve understanding of climate impacts on soil microbial communities.

Effects of climate legacies on above- and belowground community assembly

The role of climatic legacies in regulating community assembly of above- and belowground species in terrestrial ecosystems remains largely unexplored and poorly understood. Here, we report on two separate regional and continental empirical studies, including >500 locations, aiming to identify the relative importance of climatic legacies (climatic anomaly over the last 20,000 years) compared to current climates in predicting the relative abundance of ecological clusters formed by species strongly co-occurring within two independent above- and belowground networks. Climatic legacies explained a significant portion of the variation in the current community assembly of terrestrial ecosystems (up to 15.4%) that could not be accounted for by current climate, soil properties, and management. Changes in the relative abundance of ecological clusters linked to climatic legacies (e.g., past temperature) showed the potential to indirectly alter other clusters, suggesting cascading effects. Our work illustrates the role of climatic legacies in regulating ecosystem community assembly and provides further insights into possible winner and loser community assemblies under global change scenarios.

Phylogenetic conservatism of thermal traits explains dispersal limitation and genomic differentiation of Streptomyces sister-taxa

The latitudinal diversity gradient is a pattern of biogeography observed broadly in plants and animals but largely undocumented in terrestrial microbial systems. Although patterns of microbial biogeography across broad taxonomic scales have been described in a range of contexts, the mechanisms that generate biogeographic patterns between closely related taxa remain incompletely characterized. Adaptive processes are a major driver of microbial biogeography, but there is less understanding of how microbial biogeography and diversification are shaped by dispersal limitation and drift. We recently described a latitudinal diversity gradient of species richness and intraspecific genetic diversity in Streptomyces by using a geographically explicit culture collection. Within this geographically explicit culture collection, we have identified Streptomyces sister-taxa whose geographic distribution is delimited by latitude. These sister-taxa differ in geographic distribution, genomic diversity, and ecological traits despite having nearly identical SSU rRNA gene sequences. Comparative genomic analysis reveals genomic differentiation of these sister-taxa consistent with restricted gene flow across latitude. Furthermore, we show phylogenetic conservatism of thermal traits between the sister-taxa suggesting that thermal trait adaptation limits dispersal and gene flow across climate regimes as defined by latitude. Such phylogenetic conservatism of thermal traits is commonly associated with latitudinal diversity gradients for plants and animals. These data provide further support for the hypothesis that the Streptomyces latitudinal diversity gradient was formed as a result of historical demographic processes defined by dispersal limitation and driven by paleoclimate dynamics.

Comment on "Climate legacies drive global soil carbon stocks in terrestrial ecosystems"

Delgado-Baquerizo et al. (Science Advances, 12 April 2017, e1602008) use statistical correlations to infer that paleoclimate (6000 to 22,000 years ago) is a more important driver of current soil organic carbon stocks than the current-day climate. On the other hand, a wealth of radiocarbon measurements indicates that the organic carbon in most topsoils is only a few decades to perhaps a few centuries old. These seemingly incongruous results can perhaps be reconciled by considering that the long-term pedogenic development of a soil strongly influences the physiochemical properties, which lead to stabilization of new carbon entering that soil regardless of current climate.

Response to comment on "Climate legacies drive global soil carbon stocks in terrestrial ecosystem"

The technical comment from Sanderman provides a unique opportunity to deepen our understanding of the mechanisms explaining the role of paleoclimate in the contemporary distribution of global soil C content, as reported in our article. Sanderman argues that the role of paleoclimate in predicting soil C content might be accounted for by using slowly changing soil properties as predictors. This is a key point that we highlighted in the supplementary materials of our article, which demonstrated, to the degree possible given available data, that soil properties alone cannot account for the unique portion of the variation in soil C explained by paleoclimate. Sanderman also raised an interesting question about how paleoclimate might explain the contemporary amount of C in our soils if such a C is relatively new, particularly in the topsoil layer. There is one relatively simple, yet plausible, reason. A soil with a higher amount of C, a consequence of accumulation over millennia, might promote higher contemporary C fixation rates, leading to a higher amount of new C in our soils. Thus, paleoclimate can be a good predictor of the amount of soil C in soil, but not necessarily of its age. In summary, Sanderman did not question the validity of our results but rather provides an alternative potential mechanistic explanation for the conclusion of our original article, that is, that paleoclimate explains a unique portion of the global variation of soil C content that cannot be accounted for by current climate, vegetation attributes, or soil properties.

Why decadal to century timescale palaeoclimate data are needed to explain present-day patterns of biological diversity and change

The current distribution of species, environmental conditions and their interactions represent only one snapshot of a planet that is continuously changing, in part due to human influences. To distinguish human impacts from natural factors, the magnitude and pace of climate shifts, since the Last Glacial Maximum, are often used to determine whether patterns of diversity today are artefacts of past climate change. In the absence of high-temporal resolution palaeoclimate reconstructions, this is generally done by assuming that past climate change occurred at a linear pace between widely spaced (usually, ≥1,000 years) climate snapshots. We show here that this is a flawed assumption because regional climates have changed significantly across decades and centuries during glacial-interglacial cycles, likely causing rapid regional replacement of biota. We demonstrate how recent atmosphere-ocean general circulation model (AOGCM) simulations of the climate of the past 21,000 years can provide credible estimates of the details of climate change on decadal to centennial timescales, showing that these details differ radically from what might be inferred from longer timescale information. High-temporal resolution information can provide more meaningful estimates of the magnitude and pace of climate shifts, the location and timing of drivers of physiological stress, and the extent of novel climates. They also produce new opportunities to directly investigate whether short-term climate variability is more important in shaping biodiversity patterns rather than gradual changes in long-term climatic means. Together, these more accurate measures of past climate instability are likely to bring about a better understanding of the role of palaeoclimatic change and variability in shaping current macroecological patterns in many regions of the world.