Cobamides are essential nutrients for most organisms but are only biosynthesized by a limited number of taxa through aerobic and anaerobic pathways. Although the origin of these widespread shared cofactors changed ecosystems and the metabolisms of living organisms, little evolutionary information is available about the biosynthetic genes and the producers. Here, we established a timeframe for the emergence of cobamide biosynthesis genes and producers, using a series of Bayesian molecular clock analyses combined with phylogenetic reconciliation. We revealed the partial producers of tetrapyrrole precursor and corrin ring were earlier than the oldest cobamide producer, suggesting the possibility that cobamide-like compounds may have existed before the emergence of their de novo producers. We also found that the anaerobic de novo producers and corrin ring producers emerged first (Pelobacter, around 2458 Mya), and that the Great Oxidation Event postdated emergence of aerobic producers (Kribbella, around 1784 Mya). These findings reveal the chronology of cobamide biosynthesis, which greatly changed global ecological frameworks and resulted in the current biosphere, and can guide the exploration of cross-feeding and the origin of diverse organisms on the planet.

Nutrient crossfeeding critically governs microbiome-host interactions and ecosystem stability. Cobamides, synthesized only by prokaryotes, offer a powerful and tractable model for studying nutrient-mediated interdependencies in soil food webs; however, their ecological role in sustaining soil health remains unclear. Here, we construct the Soil Cobamide Producer database (SCP v.1.0) by integrating over 48,000 metagenomic and genomic datasets from 1,123 sampling sites. This database catalogs phylogenetically diverse prokaryotes (19 phyla, 302 genera) with cobamide biosynthetic potential. Using this resource, we identify host-specific colonization patterns of cobamide-producing microbes in fauna. These microbes also carry diverse functional traits that may contribute to trophic cascades and microbial community stability. In an Enchytraeid model, these colonizers support host development, modulate gene expression, and promote gut stability through transkingdom interactions, with cobamide biosynthesis serving as one representative trait among multiple microbial functions. At macroecological scales, cobamide-producing microbes occur across relatively high trophic levels, reflecting a broader principle of nutrient transfer that may also apply to other essential metabolites. This framework provides a general basis for studying nutritional microbes in soil food webs and advances One Health research.

Most studies focus on how microorganisms respond to environmental changes, yet much less is known about how microorganisms alter their surroundings. In this study, we employed Bacillus, a model genus known for their resilience and broad application, to investigate the hypothesis that microorganisms can alter environmental conditions to better suit their needs. Using a combination of global soil dataset analysis and controlled laboratory experiments, we explored the adaptive strategies and mechanisms by which microorganisms respond to varying soil pH conditions. Our results from global soil dataset analysis revealed that soil pH significantly influenced the abundance, diversity, and functional capacity of Bacillus, with the highest abundance observed at around pH 6.5. When inoculated in acidic or alkaline soil for 35 days, Bacillus shifted the pH towards neutrality. Compared with control (CK), only Bacillus inoculation significantly decreased alkaline soil pH from 8.13 to 7.36, and raised acidic soil pH by 5.86-5.91. Metabolic profiling indicated pH-dependent modulation: alkaline metabolites (e.g., laurylamine) were enriched in acidic medium of pH 5 (1529-fold), while acidic metabolites (e.g., organic acids) increased in alkaline medium of pH 8 (1.5-fold), reflecting an adaptive pH stress response. Microbial community analysis revealed that the inoculation of Bacillus not only altered the native microbial compositions but also showed close associations with specific microbial taxa, which were speculated to have a potential cooperative role in regulating soil pH. Our findings highlight that microorganisms can modify their environment in predictable ways, offering valuable insights into microbe-environment interactions.

Prokaryotic genomic traits offer insights into their functional roles, evolutionary processes, and ecological interactions, but global patterns in soil microbial genomes remain poorly understood. In this study, we examined 6436 metagenome-assembled genomes (MAGs) from global soil environments to explore the driving factors of prokaryotic gene content. Through random forest analysis, we found that, among numerous potential influencing factors such as climate, soil physicochemical properties, and human activities, geographic latitude was the primary factor affecting prokaryotic gene content. Our results showed a marked decrease in gene content from the tropics to the poles, with polar MAGs containing 10.4 % and 13.3 % fewer genes than those in tropical and temperate zones, respectively. This decline correlates with shifts in key metabolic processes, such as nitrogen fixation and energy conversion. Furthermore, we assessed interspecies metabolic interactions using Metabolic Resource Overlap (MRO) and Metabolic Interaction Potential (MIP) metrics. Our analysis revealed significantly lower MRO in high-latitude microbial communities, yet comparable MIP values to those in lower latitudes, indicating that reduced competition may contribute to genomic streamlining. These findings highlight the significant influence of latitude and interspecies interactions on microbial genomic characteristics, advancing our comprehension of microbial ecological adaptations.

Cobamides, a class of essential coenzymes synthesized only by a subset of prokaryotes, are model nutrients in microbial interaction studies and play significant roles in global ecosystems. Yet, their spatial patterns and functional roles remain poorly understood. Herein, we present an in-depth examination of cobamide-producing microorganisms, drawn from a comprehensive analysis of 2862 marine and 2979 soil metagenomic samples. A total of 1934 nonredundant metagenome-assembled genomes (MAGs) potentially capable of producing cobamides de novo were identified. The cobamide-producing MAGs are taxonomically diverse but habitat specific. They constituted only a fraction of all the recovered MAGs, with the majority of MAGs being potential cobamide users. By mapping the distribution of cobamide producers in marine and soil environments, distinct latitudinal gradients were observed: the marine environment showed peak abundance at the equator, whereas soil environments peaked at mid-latitudes. Importantly, significant and positive links between the abundance of cobamide producers and the diversity and functions of microbial communities were observed, as well as their promotional roles in essential biogeochemical cycles. These associations were more pronounced in marine samples than in soil samples, which suggests a heightened propensity for microorganisms to engage in cobamide sharing in fluid environments relative to the more spatially restricted soil environment. These findings shed light on the global patterns and potential ecological roles of cobamide-producing microorganisms in marine and soil ecosystems, enhancing our understanding of large-scale microbial interactions.

The distribution of microbial communities along elevational gradients has been widely reported, but whether there are linkages between soil microbial diversity, community assembly processes and microbial functions in elevational patterns, and whether the patterns are influenced by micro-topography (e.g. slope direction), remain unclear. Here, we collected soils from two slope directions (i.e. north-facing slope and east-facing slope) along the elevational gradient (similar to 2000 m elevational range) of a mountain ecosystem. Unexpectedly, for bacterial alpha diversity or microbial functions (i.e. multifunctionality), we found consistent elevational patterns for both slope directions, with alpha diversity showing a hump-shaped pattern and microbial functions exhibiting a linearly increasing trend. Similarly, for both slope directions, there were significant elevational distance-decay relationships in the bacterial beta diversity pattern. The elevational patterns of soil bacterial diversity were determined by the transitions of community assembly processes, while microbial functions were mainly influenced by bacterial community composition. These results suggest that soil bacterial diversity, community assembly processes and microbial functions are interactive along an elevational gradient, and that these interactions are independent of slope direction. Our study adds new insights into the elevational patterns of soil bacterial communities and their drivers and ecological consequences in mountain ecosystems.

Human activities have profoundly altered the Earth's phosphorus (P) cycling process and its associated microbial communities, yet their global distribution pattern and response to human influences remain unclear. Here, we estimated the abundances of P-cycling genes from 3321 global soil metagenomic samples and mapped the global distribution of five key P-cycling processes, that is, organic phosphoester hydrolysis, inorganic phosphorus solubilization, two-component system, phosphotransferase system, and transporters. Structural equation modeling and random forest analysis were employed to assess the impact of anthropogenic and environmental factors on the abundance of P-cycling genes. Our findings suggest that although less significant than the climate and soil profile, human-related factors, such as economic activities and population, are important drivers for the variations in P-cycling gene abundance. Notably, the gene abundances were increased parallel to the extent of human intervention, but generally at low and moderate levels of human activities. Furthermore, we identified critical genera, such as Pseudomonas and Lysobacter, which were sensitive to the changes in human activities. This study provides insights into the responses of P-cycling microbes to human activities at a global scale, enhancing our understanding of soil microbial P cycling and underscoring the importance of sustainable human activities in the Earth's biogeochemical cycle.

Nanoparticulate-Nd2O3 (nano-Nd2O3) has been excessively utilized in agriculture, industry, and medicine. Hence, nano-Nd2O3 can have environmental implications. However, the impact of nano-Nd2O3 on alpha diversity, composition, and function of soil bacterial communities has not been thoroughly evaluated. We amended soil to achieve different concentrations of nano-Nd2O3 (0, 10, 50, and 100 mg kg−1 soil) and incubated the mesocosms for 60 days. On days 7 and 60 of the experiment, we measured the effect of nano-Nd2O3 on alpha diversity and composition of soil bacterial community. Further, the effect of nano-Nd2O3 on the function of soil bacterial community was assessed based on changes in the activities of the six potential enzymes that mediate the cycling of nutrients in the soil. Nano-Nd2O3 did not alter the alpha diversity and composition of the soil bacterial community; however, it negatively affected community function in a dose-dependent manner. Specifically, the activities of β-1,4-glucosidase and β-1,4-n-acetylglucosaminidase that mediate soil carbon and nitrogen cycling, respectively, were significantly affected on days 7 and 60 of the exposure. The effect of nano-Nd2O3 on the soil enzymes correlated with changes in relative abundances of the rare and sensitive taxa, viz., Isosphaerales, Isosphaeraceae, Ktedonobacteraceae, and Streptomyces. Overall, we provide information for the safe implementation of technological applications that use nano-Nd2O3.

Parkinson's disease (PD), a progressive and incurable neurodegenerative disease, has taken a huge economic toll and medical burden on our society. Increasing evidence has shown a strong link between PD and the gut microbiome, but studies on the relationship between the gut microbiome and the severity of PD are limited. In this study, 90 fecal samples were collected from newly diagnosed and untreated patients with PD (n = 47) and matched healthy control subjects (n = 43). The 16S rRNA amplicon and shotgun metagenomic sequencing was performed, aiming to uncover the connection between the gut microbiome and disease severity in PD. The results showed that Desulfovibrio was significantly increased in PD compared to healthy controls and positively correlated with disease severity. The increase in Desulfovibrio was mainly driven by enhanced homogeneous selection and weakened drift. Moreover, through metagenome-assembled genomes (MAGs) analysis, a Desulfovibrio MAG (MAG58) was obtained which was also positively correlated with disease severity. MAG58 possesses a complete assimilatory sulfate reduction pathway and a near-complete dissimilatory sulfate reduction pathway to produce hydrogen sulfide which may influence the development of PD. Based on these results, a potential pathogenic mechanism was presented to illustrate how the increased Desulfovibrio accelerates the development of PD by producing excessive hydrogen sulfide. The present study highlighted the vital role of Desulfovibrio in the development of PD, which may provide a new target for the diagnosis and treatment of PD. KEY POINTS: • The evidence for the link between increased Desulfovibrio and disease severity in PD • A Desulfovibrio MAG was obtained which was correlated with PD • A model was presented to illustrate how increased Desulfovibrio causes PD.

Soil bacteria, which are active in shrub encroachment, play key roles in regulating ecosystem structure and function. However, the differentiation characteristics and assembly process of bacterial communities in scrubbed grasslands remain unknown. Taking the Qinghai-Tibet Plateau, a hotspot of shrub encroachment, as the study area, we collected 192 soils near nine natural typical shrubs' roots on a trans-longitude transect (about 1800 km) and investigated the bacterial communities using 16S rRNA amplicon sequencing. We found that the bacterial communities exhibited plant-specific and geographic-specific differentiation. On the one hand, bacterial communities differed significantly across plant species, with widely distributed shrubs harboring high diversity communities but few plant-specific taxa, and narrowly distributed shrubs possessing low diversity communities but more plant-specific taxa. Besides, there was a significant negative correlation between bacterial community similarity and plant phylogenetic distance. On the other hand, bacterial communities differed across geographic sites, with a significant decay in bacterial community similarity with geographic distance. The bacterial alpha diversity varied in an inverted V-shape from west to east, peaking at 91°E, which could be largely driven by mean annual temperature, soil pH and soil total carbon content. Community differentiation increased with the heterogeneity degree of assembly processes, and the dominant assembly process in these two specific differentiations differed. Dominated by stochastic and deterministic forces, respectively, geography diverged bacterial communities primarily through increased dispersal limitation, whereas plants diverged bacterial communities primarily through increased variable selection. Our study provides new insight into the characteristics and mechanisms of root-surrounding soil bacteria differentiation in scrubbed grasslands, contributing to the scientific management of degraded grasslands and the prediction of bacterial community structure and ecosystem function in response to global change.

The Three-River Headwaters region is a hotspot for studying the response of soil function to climate change. To study the horizontal variation characteristics of alpine grassland soil function and vertical changes along soil genetic horizons, soil functional indicators (including respiration, nitrogen conversion rate, and enzymatic activity) of different genetic horizons in alpine grassland soil profiles and their correlations with environmental factors were analyzed. The results showed that there were no significant differences in soil functional characteristics between alpine meadows and steppes, and topsoil had higher respiration rates, nitrogen conversion rates, and enzymatic activities than those of subsoil. Total nitrogen was a key driver of soil functional characteristics in different genetic horizons, explaining 18.3%, 21.4%, and 27.5% of the horizontal variation in functional characteristics, respectively. Climate and vegetation factors mainly affected soil function indirectly by changing soil physicochemical properties in topsoil, but atmospheric nitrogen deposition still affected soil function in subsoil. These results indicate the significant nitrogen limitation of alpine grassland soil in the Three-River Headwaters region, and the findings provide a new insight into the maintenance of soil functional diversity and the response to climate change in the context of global climate change.

Soil respiration and extracellular enzyme activity are important components of the material cycle of mountain ecosystems and play key roles in maintaining ecosystem functions. To explore the coupling relationship between soil functions and environmental factors, the soil functional indicators, environmental factors, and effects of altitude on the soil function of 36 soil samples from 12 altitudes of the Meili Mountain were analyzed. The results showed that there were significant differences in soil respirations and enzyme activities among altitudes of Meili Mountain, and high-altitude areas had higher soil functions. Soil functions increased with altitudinal difference. PCA analysis showed that the first three axes explained 56.7%, 17.4%, and 8.7% of the variance in soil functional elevation change, respectively, indicating that the functional changes related to carbon and phosphorus were higher than those related to nitrogen. There were significant correlations between environmental factors and soil functional indicators; soil function indicators had stronger correlations with soil physicochemical properties than with climatic factors. Altitude mainly affected soil function indirectly by affecting soil physicochemical properties and climatic factors. These results have great scientific significance for improving the understanding of the material cycle and ecological function of the Meili Mountain ecosystem and provide an important reference for in-depth study of the altitude distribution pattern and evolution characteristics of the soil function of the mountain ecosystem.

Cobamides are required by most organisms but are only produced by specific prokaryotic taxa. These commonly shared cofactors play significant roles in shaping the microbial community and ecosystem function. Wastewater treatment plants (WWTPs) are the world's most common biotechnological systems; knowledge about sharing of cobamides among microorganisms is predicted to be important to decipher the complex microbial relationships in these systems. Herein, we explored prokaryotic potential cobamide producers in global WWTP systems based on metagenomic analyses. A set of 8253 metagenome-assembled genomes (MAGs) were recovered and 1276 (15.5%) of them were identified as cobamide producers, which could potentially be used for the practical biological manipulation of WWTP systems. Moreover, 8090 of the total recovered MAGs (98.0%) contained at least one enzyme family dependent on cobamides, indicating the sharing of cobamides among microbial members in WWTP systems. Importantly, our results showed that the relative abundance and number of cobamide producers improved the complexity of microbial co-occurrence networks and most nitrogen, sulfur, and phosphorus cycling gene abundances, indicating the significance of cobamides in microbial ecology and their potential function in WWTP systems. These findings enhance the knowledge of cobamide producers and their functions in WWTP systems, which has important implications for improving the efficiency of microbial wastewater treatment processes.

Understanding the mechanisms of biological invasion is critical to biodiversity protection. Previous studies have produced inconsistent relationships between native species richness and invasibility, referred to as the invasion paradox. Although facilitative interactions among species have been proposed to explain the non-negative diversity–invasibility relationship, little is known about the facilitation of plant-associated microbes in invasions. We established a two-year field biodiversity experiment with a native plant species richness gradient (1, 2, 4, or 8 species) and analyzed the effects of community structure and network complexity of leaf bacteria on invasion success. Our results indicated a positive relationship between invasibility and network complexity of leaf bacteria of the invader. Consistent with previous studies, we also found that native plant species richness increased the leaf bacterial diversity and network complexity. Moreover, the results of the leaf bacteria community assembly of the invader suggested that the complex bacteria community resulted from higher native diversity rather than higher invader biomass. We concluded that increased leaf bacterial network complexity along the native plant diversity gradient likely facilitated plant invasion. Our findings provided evidence of a potential mechanism by which microbes may affect the plant community invasibility, hopefully helping to explain the non-negative relationship between native diversity and invasibility.

Biodiversity is known to affect ecosystem resistance and have implications for the maintenance of ecosystem functions and services under climate change. Compared to numbers of studies focusing on above‐ground vegetation, the response of below‐ground communities to abiotic stresses along plant diversity gradients is often ignored and is considered an important knowledge gap in ecosystem ecology. Here we conducted an integrative research to evaluate the resistance of plant biomass, and soil microbial communities and associated functional profiles to drought under varying plant diversity. We carried out a 3‐year manipulation experiment by factorially controlling plant diversity gradient (1, 2, 4 and 8 species richness) and soil moisture treatment (drought and non‐drought), and investigated the responses of plant biomass, soil bacterial and fungal diversity and community composition, soil glomalin, and five key soil enzymes. We found that plant diversity significantly improved the resistance of soil fungal communities and microbial functional profiles characterized by soil glomalin and five key enzymes, which was partly driven by the availability and accessibility of soil resources (e.g. soil moisture and organic matter) mediated by plant diversity. Furthermore, our results indicated that the enhanced resistance of fungal communities was consistent with ecological insurance theory that diverse fungal communities at high plant diversity had a higher probability of containing taxa that adapt to drought. Synthesis . Our study provides novel empirical insights into the mechanism underlying the regulatory effect of plant diversity on resistance of above‐ground vegetation and below‐ground biota to drought, with implications for understanding ecosystem response to climate change and improving biodiversity conservation practices.

Under the influence of global change and overgrazing, shrub encroachment is accelerating, which threatens ecological security and economic development in pastoral areas. Soil bacteria play key roles in maintaining ecosystem structure and function. However, the assembly process of bacterial communities in scrubbed grasslands remains unknown. Taking the Qinghai-Tibet Plateau, a hotspot of global change and shrub encroachment, as the study area, we collected 192 soils near nine natural typical shrubs’ roots on the trans-longitude transect (about 1 800 km) to evaluate the differentiation characteristics of bacterial communities, quantitatively assessed the community construction mechanisms across plant species and geographic sites mediated by assembly processes. We found that the bacterial communities exhibited plant-specific and geographic-specific differentiation. Dominated by stochastic force and deterministic force respectively, geography and plant diverge bacterial communities through increased dispersal limitation and variable selection. Our study provides new insight into the characteristics and mechanisms of root-surrounding soil bacteria differentiation in scrubbed grasslands, contributing to the scientific management of degraded grasslands and the prediction of bacterial community structure and ecosystem function in response to global change.

As the second-largest neurodegenerative disease in the world, Parkinson’s disease (PD) has brought a severe economic and medical burden to our society. Growing evidence in recent years suggests that the gut microbiome may influence PD, but the exact pathogenesis of PD remains unclear. In addition, the current diagnosis of PD could be inaccurate and expensive. In this study, the largest meta-analysis currently of the gut microbiome in PD was analyzed, including 2269 samples by 16S rRNA gene and 236 samples by shotgun metagenomics, aiming to reveal the connection between PD and gut microbiome and establish a model to predict PD. The results showed that the relative abundances of potential pro-inflammatory bacteria, genes and pathways were significantly increased in PD, while potential anti-inflammatory bacteria, genes and pathways were significantly decreased. These changes may lead to a decrease in potential anti-inflammatory substances (short-chain fatty acids) and an increase in potential pro-inflammatory substances (lipopolysaccharides, hydrogen sulfide and glutamate). Notably, the results of 16S rRNA gene and shotgun metagenomic analysis have consistently identified five decreased genera ( Roseburia , Faecalibacterium , Blautia , Lachnospira, and Prevotella ) and five increased genera ( Streptococcus , Bifidobacterium , Lactobacillus , Akkermansia, and Desulfovibrio ) in PD. Furthermore, random forest models performed well for PD prediction based on 11 genera (accuracy > 80%) or 6 genes (accuracy > 90%) related to inflammation. Finally, a possible mechanism was presented to explain the pathogenesis of inflammation leading to PD. Our results provided further insights into the prediction and treatment of PD based on inflammation.

Soil microbial biogeographical patterns have been widely explored from horizontal to vertical scales. However, studies of microbial vertical distributions were still limited (e.g., how soil genetic horizons influence microbial distributions). To shed light on this question, we investigated soil bacterial communities across three soil horizons (topsoil: horizon A; midsoil: horizon B; subsoil: horizon C) of 60 soil profiles along a 3500 km transect in the Qinghai-Tibet Plateau. We found that bacterial diversity was highest in the topsoil and lowest in the subsoil, and community composition significantly differed across soil horizons. The network complexity decreased from topsoil to subsoil. There were significant geographical/environmental distance-decay relationships (DDR) in three soil horizons, with a lower slope from topsoil to subsoil due to the decreased environmental heterogeneity. Variation partitioning analysis (VPA) showed that bacterial community variations were explained more by environmental than spatial factors. Although environmental selection processes played a dominant role, null model analysis revealed that deterministic processes (mainly variable selection) decreased with deeper soil horizons, while stochastic processes (mainly dispersal limitation) increased from topsoil to subsoil. These results suggested that microbial biogeographical patterns and community assembly processes were soil horizon dependent. Our study provides new insights into the microbial vertical distributions in large-scale alpine regions and highlights the vital role of soil genetic horizons in affecting microbial community assembly, which has implications for understanding the pedogenetic process and microbial responses to extreme environment under climate change.

Denitrification by microorganisms in soil regulates ecosystem nitrogen availability and cycling. Although arbuscular mycorrhiza fungi (AMF) are best known as the key connectors between plant and microorganisms in below-ground soil, however, less knowledge is available about the interactive effects of AMF and vegetation community traits on soil denitrifiers. In this study, we manipulated experimental pots with or without AMF inoculation (AMF+ and AMF− treatments) under three plant richness levels and 27 different plant community compositions. Our results provided evidence that inoculation of AMF significantly increased the abundances of nirS-, nirK- and nosZ-type denitrifiers as revealed by quantitative PCR (qPCR), which inferred a positive role of AMF in N-cycling microorganisms. Plant community traits, including richness, community composition, biomass and species, were less important in influencing the abundances of soil denitrifiers, and no significant interactive effect was detected between AMF inoculation and plant richness or plant community compositions, indicating a weak direct relationship between plants and soil denitrifiers. Our study provided comprehensive insights into the roles of AMF-plant associate bio-system in driving the variation of soil denitrifiers.

Observational and experimental evidence has revealed the functional importance of microbial diversity. However, the effects of microbial diversity loss on ecosystem functions are not consistent across studies, which are probably tempered by microbial functional redundancy, specific taxa and functions evaluated. Here we conducted diversity manipulation experiments in two independent soils with distinct prokaryotic communities, and investigated how the initial community traits (e.g., distinct functional redundancy and taxonomic composition) modulate the contribution of prokaryotic diversity loss and composition shift to eight ecosystem functions related to soil nutrient cycling. We found that diversity loss impaired three functions (potential nitrification rate, N2 -fixation activity and phosphatase) and multifunctionality only in the communities with low functional redundancy, but all examined functions were unaffected in the communities with high functional redundancy. All significantly affected functions belonged to specialized functions, while the broad function (soil basal respiration) was unaffected. Moreover, prokaryotic composition explained more functional variation than diversity, which was ascribed to the crucial role of specific taxa that influence particular functions. Taken together, this study provides empirical evidence for identifying the mechanism underlying the ecosystem response to changes in microbial community, with implications for improving the prediction of ecosystem process models and managing microbial communities to promote ecosystem services.