Incorporating phylogenetic metrics to microbial co‐occurrence networks based on amplicon sequences to discern community assembly processes

Microplastic type and concentration affect prokaryotic community structure and species coexistence in deep-sea cold seep sediments

As important methane hydrate storage sites, cold seep areas are threatened by microplastics (MPs) contamination. To assess the environmental impact of MPs on microbial communities in cold seep sediments, an incubation experiment was conducted using cold seep sediment amended with different concentration of polyamide (PA), polyethylene (PE), polyethylene terephthalate (PET), and polypropylene (PP) microplastics. The results showed that the different type and concentration of MPs significantly altered the prokaryotic community structures. The PE and PET addition highly changed the relative abundance of bacterial taxa in the bacterial community, while the proportion of archaeal species in the archaeal community was significantly altered in 0.5 % MPs treatments. All of the MPs reduced the network complexity of the bacterial and archaeal communities, such as the lower average degree and greater average path length. Furthermore, the MPs treatments also significantly decreased the network stability of prokaryotic communities. The lower network complexity led to lower network stability was observed in the archaeal community. The formation of oxidative functional groups on PE and PET MP surface based on FTIR analysis suggested that biodegradation could occur in cold seep sediment. Together, these results provide new evidence that MPs could change the structures and species coexistence of prokaryotic communities in cold seep sediments.

Microbial Biotic Associations Dominated Adaptability Differences of Dioecious Poplar Under Salt Stress

How different stress responses by male and female plants are influenced by interactions with rhizosphere microbes remains unclear. In this study, we employed poplar as a dioecious model plant and quantified biotic associations between microorganisms to explore the relationship between microbial associations and plant adaptation. We propose a health index (HI) to comprehensively characterize the physiological characteristics and adaptive capacity of plants under stress. It was found that male poplars demonstrated higher salt stress tolerance than females, and root-secreted citric acid was significantly higher in the rhizospheres of male poplars. Positive biotic association among bacteria increased poplar HI significantly under salt stress, while fungal and cross-domain biotic association (bacteria-fungi) did not. We further identified a keystone bacterial taxon regulating bacterial biotic association, ASV_22706, which was itself regulated by citric acid and significantly positively correlated with host HI. The abundance of keystone fungal taxa was positively correlated with HI of male poplars and negatively correlated with HI of female poplars. Compared with female poplars, male poplars enriched more prebiotics and probiotics under stress. This work primarily reveals the relationship between adaptation differences and microbial interactions in dioecious plants, which suggests a microbial approach to improve plant adaptability to stress conditions.

Distinct strategies of microeukaryotic generalists and specialists in Qinghai-Tibet plateau sediment driven by salinity

Unraveling how microeukaryotic generalists and specialists assemble and coexist under environmental stress is central to our understanding of the mechanisms maintaining diversity. Here, we explored the biogeographical distributions of microeukaryotic generalists and specialists in lake surface sediments along a salinity gradient on the Qinghai-Tibet Plateau. We found that relative abundances of Chlorophyta (28.6 %) and Dinophyceae (9.5 %) were higher as habitat generalists than as specialists. Conversely, relative abundances of habitat specialists were higher in the Ciliophora (22.2 %) and Cercozoa (11.6 %) than those of generalists. Environmental adaptation analysis showed a broader niche threshold for generalists than for specialists, whereas a stronger phylogenetic signal for environmental factors was observed for specialists. Thus, increases in salinity had stronger effects on specialists than on generalists through environmental selection and diversification processes. However, null model analysis indicated stochastic processes were the primary drivers of both generalists and specialists. Network analysis revealed that with increasing salinity, specialists were more important than generalists in stabilizing networks. In addition, phylogenetic relatedness indicated that microeukaryotic generalists coexisted because of niche differences, whereas specialists coexisted because of average fitness similarity. Our study will help to predict microeukaryotic responses to environmental changes in aquatic ecosystems.

Analyzing bacterial networks and interactions in skin and gills of Sparus aurata with microalgae-based additive feeding

The inclusion of microalgae in functional fish diets has a notable impact on the welfare, metabolism and physiology of the organism. The microbial communities associated with the fish are directly influenced by the host's diet, and further understanding the impact on mucosal microbiota is needed. This study aimed to analyze the microbiota associated with the skin and gills of Sparus aurata fed a diet containing 10% microalgae. Sequencing of the V3-V4 variable region of 16S rDNA molecules was employed to determine the composition of the microbial communities. The study employed bioinformatics tools to explore the taxonomic composition and interactions of the microbiota, emphasizing the integration of taxonomic analysis and abundance correlation networks as crucial for understanding microbial community dynamics and the impact of functional diets. The results indicated that there were not changes in the composition of the skin and gill microbiota. However, notable differences were observed in the bacterial interaction networks. The skin and gill networks exhibited distinct overall patterns influenced by their respective environments and functions. The gill network showed a highly connected and redundant structure, increasing its resilience, while the skin network showed a more fragmented structure, suggesting a potentially greater vulnerability to perturbations. Key taxa, such as Acinetobacter and Polaribacter, were identified as critical for maintaining the stability and functionality of these microbial ecosystems. Polaribacter could demonstrate potential protection against pathogens through negative interactions. These tentative studies open up an additional avenue to consider, such as the interactions among bacterial communities, as well as new proposals to corroborate these findings. These observations underline the importance of understanding the composition and interactions within bacterial groups to fully grasp the dynamics of microbial communities.

Community Assembly Mechanisms of nirK- and nirS-type Denitrifying Bacteria in Sediments of Eutrophic Lake Taihu, China

Denitrifying bacteria, particularly nirK- and nirS-type, are functionally equivalent and could occupy different niches, but their community assembly mechanisms and responses to environmental heterogeneity are poorly understood in eutrophic lakes. In this study, we investigated the community assembly mechanisms of nirK- and nirS-type denitrifying bacteria and clarified their responses to sediments environmental factors in Lake Taihu, China. The quantitative real-time PCR and Illumina HiSeq-based sequencing revealed that the abundance and composition of two types of denitrifying bacterial communities varied among different sites in the sediments of Lake Taihu. The functions of these two types of denitrifying bacteria were assigned to mainly nitrogen cycling along with carbon, oxygen, and sulfur cycling, indicating their diverse ecosystems functions. Neutral community model showed that majority of nirK- and nirS-type denitrifying bacteria were neutrally distributed, while dispersal and selection were the dominant drivers in shaping community assembly of nirK-type bacteria. The community assembly of nirS-type was mainly driven by homogeneous selection. We found complex network interactions between nirK- and nirS-type denitrifying bacteria with other bacterial communities, indicating the importance of other bacterial coexistence for ecosystem functions by denitrifying bacteria in lake sediments. Keystone taxa of other bacteria showed the highest interactions with denitrifying bacteria; further, a strong significant correlation between keystone taxa with environmental factors and sediment enzyme content revealed by Mantel tests. Specially, total phosphorous was the key environmental factor determining the composition and diversity of nirK and nirS-type denitrifying bacteria in lake sediments, whereas NAR, AmoA, and NIR were the key reductase enzymes directly or indirectly affected to them. Our results provide significant insights into understanding the effects of changing nirK- and nirS-type denitrifying bacterial diversities and underlying community assembly mechanisms under changing environmental conditions in eutrophic lake ecosystems.

Deep insights into the assembly mechanisms, co-occurrence patterns, and functional roles of microbial community in wastewater treatment plants

The understanding of activated sludge microbial status and roles is imperative for improving and enhancing the performance of wastewater treatment plants (WWTPs). In this study, we conducted a deep analysis of activated sludge microbial communities across five compartments (inflow, effluent, and aerobic, anoxic, anaerobic tanks) over temporal scales, employing high-throughput sequencing of 16S rRNA amplicons and metagenome data. Clearly discernible seasonal patterns, exhibiting cyclic variations, were observed in microbial diversity, assembly, co-occurrence network, and metabolic functions. Notably, summer samples exhibited higher α-diversity and were distinctly separated from winter samples. Our analysis revealed that microbial community assembly is influenced by both stochastic processes (66%) and deterministic processes (34%), with winter samples demonstrating more random assembly compared to summer. Co-occurrence patterns were predominantly mutualistic, with over 96% positive correlations, and summer networks were more organized than those in winter. These variations were significantly correlated with temperature, total phosphorus and sludge volume index. However, no significant differences were found among microbial community across five compartments in terms of β diversity. A core community of keystone taxa was identified, playing key roles in eight nitrogen and eleven phosphorus cycling pathways. Understanding the assembly mechanisms, co-occurrence patterns, and functional roles of microbial communities is essential for the design and optimization of biotechnological treatment processes in WWTPs.

Networks as tools for defining emergent properties of microbiomes and their stability

The potential promise of the microbiome to ameliorate a wide range of societal and ecological challenges, from disease prevention and treatment to the restoration of entire ecosystems, hinges not only on microbiome engineering but also on the stability of beneficial microbiomes. Yet the properties of microbiome stability remain elusive and challenging to discern due to the complexity of interactions and often intractable diversity within these communities of bacteria, archaea, fungi, and other microeukaryotes. Networks are powerful tools for the study of complex microbiomes, with the potential to elucidate structural patterns of stable communities and generate testable hypotheses for experimental validation. However, the implementation of these analyses introduces a cascade of dichotomies and decision trees due to the lack of consensus on best practices. Here, we provide a road map for network-based microbiome studies with an emphasis on discerning properties of stability. We identify important considerations for data preparation, network construction, and interpretation of network properties. We also highlight remaining limitations and outstanding needs for this field. This review also serves to clarify the varying schools of thought on the application of network theory for microbiome studies and to identify practices that enhance the reproducibility and validity of future work. Video Abstract.

Effects of Pre-Dehydration Treatments on Physicochemical Properties, Non-Volatile Flavor Characteristics, and Microbial Communities during Paocai Fermentation

The paocai industry faces challenges related to the production of large volumes of high-salinity and acidic brine by-products. Maintaining paocai quality while reducing brine production is crucial. This study utilized high-throughput sequencing technology to analyze microbial changes throughout the fermentation process, along with the non-volatile flavor compounds and physicochemical properties, to assess the impact of hot-air and salt-pressing pre-dehydration treatments on paocai quality. The findings indicate that pre-dehydration of raw material slowed the fermentation process but enhanced the concentration of non-volatile flavor substances, including free amino acids and organic acids. Hot-air pre-dehydration effectively reduced initial salinity to levels comparable to those in high-salinity fermentation of fresh vegetables. Furthermore, pre-dehydration altered microbial community structures and simplified inter-microbial relationships during fermentation. However, the key microorganisms such as Lactobacillus, Weissella, Enterobacter, Wallemia, Aspergillus, and Kazachstania remained consistent across all groups. Additionally, this study found that biomarkers influenced non-volatile flavor formation differently depending on the treatment, but these substances had minimal impact on the biomarkers and showed no clear correlation with high-abundance microorganisms. Overall, fermenting pre-dehydrated raw materials presents an environmentally friendly alternative to traditional paocai production.

Disentangling the factors defining Bacillus subtilis group species abundance in natural soils

Bacillus subtilis is ubiquitously and broadly distributed in various environments but is mostly isolated from soil. Given that B. subtilis is known as a plant growth-promoting rhizobacterium in agriculture, we aimed to describe the natural distribution of this species and uncover how biotic and abiotic factors affect its distribution. When comparing different soils, we discovered that B. subtilis group species are most abundant in grasslands but can rarely be isolated from forest soil, even if the soil sample sites are situated in proximity. Differential analysis revealed that spore-forming bacteria exhibited enrichments in the grassland, suggesting niche overlap or synergistic interactions leading to the proliferation of certain Bacillus species in grassland environments. Network analysis further revealed that Bacillus and other Bacillota established a densely interconnected hub module in the grassland, characterised by positive associations indicating co-occurrence, a pattern not observed in the forest soil. Speculating that this difference was driven by abiotic factors, we combined amplicon sequencing with physico-chemical analysis of soil samples and found multiple chemical variables, mainly pH, to affect microbial composition. Our study pinpoints the factors that influence B. subtilis abundance in natural soils and, therefore, offers insights for designing B. subtilis-based biocontrol products in agricultural settings.

Genetic hypogonadal (Gnrh1hpg) mouse model uncovers influence of reproductive axis on maturation of the gut microbiome during puberty

The gut microbiome plays a key role in human health and gut dysbiosis is linked to many sex-specific diseases including autoimmune, metabolic, and neurological disorders. Activation of the hypothalamic-pituitary-gonadal (HPG) axis during puberty leads to sexual maturation and development of sex differences through the action of gonadal sex steroids. While the gut microbiome also undergoes sex differentiation, the mechanisms involved remain poorly understood. Using a genetic hypogonadal (hpg) mouse model, we sampled the fecal microbiome of male and female wild-type and hpg mutant mice before and after puberty to determine how microbial taxonomy and function are influenced by age, sex, and the HPG axis. We showed that HPG axis activation during puberty is required for sexual maturation of the gut microbiota composition, community structure, and metabolic functions. We also demonstrated that some sex differences in taxonomic composition and amine metabolism developed independently of the HPG axis, indicating that sex chromosomes are sufficient for certain sex differences in the gut microbiome. In addition, we showed that age, independent of HPG axis activation, led to some aspects of pubertal maturation of the gut microbiota community composition and putative functions. These results have implications for microbiome-based treatments, indicating that sex, hormonal status, and age should be considered when designing microbiome-based therapeutics.

Soil microbial subcommunity assembly mechanisms are highly variable and intimately linked to their ecological and functional traits

Revealing the mechanisms underlying soil microbial community assembly is a fundamental objective in molecular ecology. However, despite increasing body of research on overall microbial community assembly mechanisms, our understanding of subcommunity assembly mechanisms for different prokaryotic and fungal taxa remains limited. Here, soils were collected from more than 100 sites across southwestern China. Based on amplicon high-throughput sequencing and iCAMP analysis, we determined the subcommunity assembly mechanisms for various microbial taxa. The results showed that dispersal limitation and homogenous selection were the primary drivers of soil microbial community assembly in this region. However, the subcommunity assembly mechanisms of different soil microbial taxa were highly variable. For instance, the contribution of homogenous selection to Crenarchaeota subcommunity assembly was 70%, but it was only around 10% for the subcommunity assembly of Actinomycetes, Gemmatimonadetes and Planctomycetes. The assembly of subcommunities including microbial taxa with higher occurrence frequencies, average relative abundance and network degrees, as well as wider niches tended to be more influenced by homogenizing dispersal and drift, but less affected by heterogeneous selection and dispersal limitation. The subcommunity assembly mechanisms also varied substantially among different functional guilds. Notably, the subcommunity assembly of diazotrophs, nitrifiers, saprotrophs and some pathogens were predominantly controlled by homogenous selection, while that of denitrifiers and fungal pathogens were mainly affected by stochastic processes such as drift. These findings provide novel insights into understanding soil microbial diversity maintenance mechanisms, and the analysis pipeline holds significant value for future research.

Spatio-temporal structuring and assembly of abundant and rare bacteria in the benthic compartment of a marginally eutrophic lagoon

The investigations on ecological processes that structure abundant and rare sub-communities are limited from the benthic compartments of tropical brackish lagoons. We examined the spatial and temporal patterns in benthic bacterial communities of a brackish lagoon; Chilika. Abundant and rare bacteria showed differences in niche specialization but exhibited similar distance-decay patterns. Abundant bacteria were mostly habitat generalists due to their broader niche breadth, environmental response thresholds, and greater functional redundancy. In contrast, rare bacteria were mostly habitat specialists due to their narrow niche breadth, lower environmental response thresholds, and functional redundancy. The spatial patterns in abundant bacteria were largely shaped by stochastic processes (88.7 %, mostly dispersal limitation). In contrast, rare bacteria were mostly structured by deterministic processes (56.4 %, mostly heterogeneous selection). These findings provided a quantitative assessment of the different forces namely spatial, environmental, and biotic that together structured bacterial communities in the benthic compartment of a marginally eutrophic lagoon.

Ecology and risks of the global plastisphere as a newly expanding microbial habitat

Plastic offers a new niche for microorganisms, the plastisphere. The ever-increasing emission of plastic waste makes it critical to understand the microbial ecology of the plastisphere and associated effects. Here, we present a global fingerprint of the plastisphere, analyzing samples collected from freshwater, seawater, and terrestrial ecosystems. The plastisphere assembles a distinct microbial community that has a clearly higher heterogeneity and a more deterministically dominated assembly compared to natural habitats. New coexistence patterns-loose and fragile networks with mostly specialist linkages among microorganisms that are rarely found in natural habitats-are seen in the plastisphere. Plastisphere microbiomes generally have a great potential to metabolize organic compounds, which could accelerate carbon turnover. Microorganisms involved in the nitrogen cycle are also altered in the plastisphere, especially in freshwater plastispheres, where a high abundance of denitrifiers may increase the release of nitrite (aquatic toxicant) and nitrous oxide (greenhouse gas). Enrichment of animal, plant, and human pathogens means that the plastisphere could become an increasingly mobile reservoir of harmful microorganisms. Our findings highlight that if the trajectory of plastic emissions is not reversed, the expanding plastisphere could pose critical planetary health challenges.

Natural restoration exhibits better soil bacterial network complexity and stability than artificial restoration on the Loess Plateau, China

Natural restoration (NR, e.g., secondary succession) and artificial restoration (AR, e.g., afforestation) are key approaches for rehabilitating degraded land; however, a comparative assessment of microbial network between these approaches is lacking. We compared bacterial networks under NR and AR in two different watersheds on the Loess Plateau. Our findings revealed significantly heightened network complexity under NR compared to AR, including metrics such as node, edge, modularity, degree, centrality, and keystone nodes. NR's network robustness exceeded AR by 19.45-35.9% and 7.79-17.74% in the two watersheds, aligning with the ecological principle that complexity begets stability. The significantly higher negative/positive cohesion and natural connectivity under NR also support its better network stability than AR. Integrated analysis of paired sequencing data from five Loess Plateau studies conducted on the Loess Plateau further confirmed the higher complexity and stability of bacterial networks under NR. Further analysis unveiled "biological interactions" as primary drivers of bacterial co-occurrence (on average 84.21% of links), surpassing the influence of environmental filtering (5.17%) or dispersal limitation (4.2%). Importantly, networked communities under NR exhibited generally stronger linkages with various ecosystem function than AR. Overall, our study provides insights into vegetation restoration strategies from the perspective of microbial network, underscoring natural regeneration's potential as a superior remedy for degraded-land restoration.

Distinct community assembly processes and habitat specialization driving the biogeographic patterns of abundant and rare bacterioplankton in a brackish coastal lagoon

The ecological diversity patterns and community assembly processes along spatio-temporal scales are least studied in the bacterioplankton sub-communities of brackish coastal lagoons. We examined the biogeographic patterns and relative influences of different assembly processes in structuring the abundant and rare bacterioplankton sub-communities of Chilika, the largest brackish water coastal lagoon of India. Rare taxa demonstrated significantly higher α- and β-diversity and biogeochemical functions than abundant taxa in the high-throughput 16S rRNA gene sequence dataset. The majority of the abundant taxa (91.4 %) were habitat generalists with a wider niche breadth (niche breadth index, B = 11.5), whereas most of the rare taxa (95.2 %) were habitat specialists with a narrow niche breadth (B = 8.9). Abundant taxa exhibited a stronger distance-decay relationship and higher spatial turnover rate than rare taxa. β-diversity partitioning revealed that the contribution of species turnover (72.2-97.8 %) was greater than nestedness (2.2-27.8 %) in causing the spatial variation in both abundant and rare taxa. Null model analyses revealed that the distribution of abundant taxa was mostly structured by stochastic processes (62.8 %), whereas deterministic processes (54.1 %) played a greater role in the rare taxa. However, the balance of these two processes varied across spatio-temporal scales in the lagoon. Salinity was the key deterministic factor controlling the variation of both abundant and rare taxa. Potential interaction networks showed a higher influence of negative interactions, indicating that species exclusion and top-down processes played a greater role in the community assembly. Notably, abundant taxa emerged as keystone taxa across spatio-temporal scales, suggesting their greater influences on other bacterial co-occurrences and network stability. Overall, this study provided detailed mechanistic insights into biogeographic patterns and underlying community assembly processes of the abundant and rare bacterioplankton over spatio-temporal scales in a brackish lagoon.

QCMI: A method for quantifying putative biotic associations of microbes at the community level

A workflow has been compiled as "qcmi" R package-the quantifying community-level microbial interactions-to identify and quantify the putative biotic associations of microbes at the community level from ecological networks.

Determinants of synergistic cell-cell interactions in bacteria

Bacteria are ubiquitous and colonize virtually every conceivable habitat on earth. To achieve this, bacteria require different metabolites and biochemical capabilities. Rather than trying to produce all of the needed materials by themselves, bacteria have evolved a range of synergistic interactions, in which they exchange different commodities with other members of their local community. While it is widely acknowledged that synergistic interactions are key to the ecology of both individual bacteria and entire microbial communities, the factors determining their establishment remain poorly understood. Here we provide a comprehensive overview over our current knowledge on the determinants of positive cell-cell interactions among bacteria. Taking a holistic approach, we review the literature on the molecular mechanisms bacteria use to transfer commodities between bacterial cells and discuss to which extent these mechanisms favour or constrain the successful establishment of synergistic cell-cell interactions. In addition, we analyse how these different processes affect the specificity among interaction partners. By drawing together evidence from different disciplines that study the focal question on different levels of organisation, this work not only summarizes the state of the art in this exciting field of research, but also identifies new avenues for future research.

Spatial scales of competition and a growth-motility trade-off interact to determine bacterial coexistence

The coexistence of competing species is a long-lasting puzzle in evolutionary ecology research. Despite abundant experimental evidence showing that the opportunity for coexistence decreases as niche overlap increases between species, bacterial species and strains competing for the same resources are commonly found across diverse spatially heterogeneous habitats. We thus hypothesized that the spatial scale of competition may play a key role in determining bacterial coexistence, and interact with other mechanisms that promote coexistence, including a growth-motility trade-off. To test this hypothesis, we let two Pseudomonas putida strains compete at local and regional scales by inoculating them either in a mixed droplet or in separate droplets in the same Petri dish, respectively. We also created conditions that allow the bacterial strains to disperse across abiotic or fungal hyphae networks. We found that competition at the local scale led to competitive exclusion while regional competition promoted coexistence. When competing in the presence of dispersal networks, the growth-motility trade-off promoted coexistence only when the strains were inoculated in separate droplets. Our results provide a mechanism by which existing laboratory data suggesting competitive exclusion at a local scale is reconciled with the widespread coexistence of competing bacterial strains in complex natural environments with dispersal.

Storm promotes the dissemination of antibiotic resistome in an urban lagoon through enhancing bio-interactions

Antibiotic-resistance genes (ARGs) and resistant bacteria (ARB) are abundant in stormwater that could cause serious infections, posing a potential threat to public health. However, there is no inference about how stormwater contributes to ARG profiles as well as the dynamic interplay between ARGs and bacteria via vertical gene transfer (VGT) or horizontal gene transfer (HGT) in urban water ecosystems. In this study, the distribution of ARGs, their host communities, and the source and community assembly process of ARGs were investigated in Yundang Lagoon (China) via high-throughput quantitative PCR, 16S rRNA gene amplicon sequencing, and application of SourceTracker before, after and recovering from an extreme precipitation event (132.1 mm). The abundance of ARGs and mobile genetic elements (MGEs) was the highest one day after precipitation and then decreased 2 days after precipitation and so on. Based on SourceTracker and NMDS analysis, the ARG and bacterial communities in lagoon surface water from one day after precipitation were mainly contributed by the wastewater treatment plant (WWTP) influent and effluent. However, the contribution of WWTP to ARG communities was minor 11 days after the precipitation, suggesting that the storm promoted the ARG levels by introducing the input of ARGs, MGEs, and ARB from point and non-point sources, such as sewer overflow and land-applied manure. Based on a novel microbial network analysis framework, the contribution of positive biological interactions between ARGs and MGEs or bacteria was the highest one day after precipitation, indicating a promoted VGT and HGT for ARG dissemination. The microbial networks deconstructed 11 days after precipitation, suggesting the stormwater practices (e.g., tide gate opening, diversion channels, and pumping) alleviated the spread of ARGs. These results advanced our understanding of the distribution and transport of ARGs associated with their source in urban stormwater runoff.

Absence of stress‐promoted facilitation coupled with a competition decrease in the microbiome of ephemeral saline lakes

Salinity fluctuations constitute a well-known high stress factor strongly shaping global biological distributions and abundances. However, there is a knowledge gap regarding how increasing saline stress affects microbial biological interactions. We applied the combination of a probabilistic method for estimating significant co-occurrences/exclusions and a conceptual framework for filtering out associations potentially linked to environmental and/or spatial factors, in a series of connected ephemeral (hyper) saline lakes. We carried out a network analysis over the full aquatic microbiome-bacteria, eukarya, and archaea-under severe salinity fluctuations. Most of the observed co-occurrences/exclusions were potentially explained by environmental niche and/or dispersal limitation. Co-occurrences assigned to potential biological interactions remained stable, suggesting that the salt gradient was not promoting interspecific facilitation processes. Conversely, co-exclusions assigned to potential biological interactions decreased along the gradient both in number and network complexity, pointing to a decrease of interspecies competition as salinity increased. Overall, higher saline stress reduced microbial co-exclusions while co-occurrences remained stable suggesting decreasing competition coupled with lack of stress-gradient promoted facilitation in the microbiome of ephemeral saline lakes.

Co-occurrence networks reveal more complexity than community composition in resistance and resilience of microbial communities

Plant response to drought stress involves fungi and bacteria that live on and in plants and in the rhizosphere, yet the stability of these myco- and micro-biomes remains poorly understood. We investigate the resistance and resilience of fungi and bacteria to drought in an agricultural system using both community composition and microbial associations. Here we show that tests of the fundamental hypotheses that fungi, as compared to bacteria, are (i) more resistant to drought stress but (ii) less resilient when rewetting relieves the stress, found robust support at the level of community composition. Results were more complex using all-correlations and co-occurrence networks. In general, drought disrupts microbial networks based on significant positive correlations among bacteria, among fungi, and between bacteria and fungi. Surprisingly, co-occurrence networks among functional guilds of rhizosphere fungi and leaf bacteria were strengthened by drought, and the same was seen for networks involving arbuscular mycorrhizal fungi in the rhizosphere. We also found support for the stress gradient hypothesis because drought increased the relative frequency of positive correlations.

Scale dependency of ectomycorrhizal fungal community assembly processes in Mediterranean mixed forests

The assembly of biological communities depends on deterministic and stochastic processes whose influence varies across spatial and temporal scales. Although ectomycorrhizal (ECM) fungi play a key role in forest ecosystems, our knowledge on ECM community assembly processes and their dependency on spatial scales is still scarce. We analysed the assembly processes operating on ECM fungal communities associated with Cistus albidus L. and Quercus spp. in Mediterranean mixed forests (Southern Spain), for which root tip ECM fungi were characterized by high-throughput sequencing. The relative contribution of deterministic and stochastic processes that govern the ECM fungal community assembly was inferred by using phylogenetic and compositional turnover descriptors across spatial scales. Our results revealed that stochastic processes had a significantly higher contribution than selection on root tip ECM fungal community assembly. The strength of selection decreased at the smallest scale and it was linked to the plant host identity and the environment. Dispersal limitation increased at finer scales, whilst drift showed the opposite pattern likely suggesting a main influence of priority effects on ECM fungal community assembly. This study highlights the potential of phylogeny to infer ECM fungal community responses and brings new insights into the ecological processes affecting the structure and dynamics of Mediterranean forests.

Biological Microbial Interactions from Cooccurrence Networks in a High Mountain Lacustrine District

A fundamental question in biology is why some species tend to occur together in the same locations, while others are never observed coexisting. This question becomes particularly relevant for microorganisms thriving in the highly diluted waters of high mountain lakes, where biotic interactions might be required to make the most of an extreme environment. We studied a high-throughput gene data set of alpine lakes (>220 Pyrenean lakes) with cooccurrence network analysis to infer potential biotic interactions, using the combination of a probabilistic method for determining significant cooccurrences and coexclusions between pairs of species and a conceptual framework for classifying the nature of the observed cooccurrences and coexclusions. This computational approach (i) determined and quantified the importance of environmental variables and spatial distribution and (ii) defined potential interacting microbial assemblages. We determined the properties and relationships between these assemblages by examining node properties at the taxonomic level, indicating associations with their potential habitat sources (i.e., aquatic versus terrestrial) and their functional strategies (i.e., parasitic versus mixotrophic). Environmental variables explained fewer pairs in bacteria than in microbial eukaryotes for the alpine data set, with pH alone explaining the highest proportion of bacterial pairs. Nutrient composition was also relevant for explaining association pairs, particularly in microeukaryotes. We identified a reduced subset of pairs with the highest probability of species interactions (“interacting guilds”) that significantly reached higher occupancies and lower mean relative abundances in agreement with the carrying capacity hypothesis. The interacting bacterial guilds could be more related to habitat and microdispersal processes (i.e., aquatic versus soil microbes), whereas for microeukaryotes trophic roles (osmotrophs, mixotrophs, and parasitics) could potentially play a major role. Overall, our approach may add helpful information to guide further efforts for a mechanistic understanding of microbial interactions in situ.

IMPORTANCE A fundamental question in biology is why some species tend to occur together in the same locations, while others are never observed to coexist. This question becomes particularly relevant for microorganisms thriving in the highly diluted waters of high mountain lakes, in which biotic interactions might be required to make the most of an extreme environment. Microbial metacommunities are too often only studied in terms of their environmental niches and geographic barriers since they show inherent difficulties to quantify biological interactions and their role as drivers of ecosystem functioning. Our study highlights that telling apart potential interactions from both environmental and geographic niches may help for the initial characterization of organisms with similar ecologies in a large scope of ecosystems, even when information about actual interactions is partial and limited. The multilayered statistical approach carried out here offers the possibility of going beyond taxonomy to understand microbiological behavior in situ.

From diversity to complexity: Microbial networks in soils

Network analysis has been used for many years in ecological research to analyze organismal associations, for example in food webs, plant-plant or plant-animal interactions. Although network analysis is widely applied in microbial ecology, only recently has it entered the realms of soil microbial ecology, shown by a rapid rise in studies applying co-occurrence analysis to soil microbial communities. While this application offers great potential for deeper insights into the ecological structure of soil microbial ecosystems, it also brings new challenges related to the specific characteristics of soil datasets and the type of ecological questions that can be addressed. In this Perspectives Paper we assess the challenges of applying network analysis to soil microbial ecology due to the small-scale heterogeneity of the soil environment and the nature of soil microbial datasets. We review the different approaches of network construction that are commonly applied to soil microbial datasets and discuss their features and limitations. Using a test dataset of microbial communities from two depths of a forest soil, we demonstrate how different experimental designs and network constructing algorithms affect the structure of the resulting networks, and how this in turn may influence ecological conclusions. We will also reveal how assumptions of the construction method, methods of preparing the dataset, and definitions of thresholds affect the network structure. Finally, we discuss the particular questions in soil microbial ecology that can be approached by analyzing and interpreting specific network properties. Targeting these network properties in a meaningful way will allow applying this technique not in merely descriptive, but in hypothesis-driven research. Analysing microbial networks in soils opens a window to a better understanding of the complexity of microbial communities. However, this approach is unfortunately often used to draw conclusions which are far beyond the scientific evidence it can provide, which has damaged its reputation for soil microbial analysis. In this Perspectives Paper, we would like to sharpen the view for the real potential of microbial co-occurrence analysis in soils, and at the same time raise awareness regarding its limitations and the many ways how it can be misused or misinterpreted.

Distinct Co-occurrence Relationships and Assembly Processes of Active Methane-Oxidizing Bacterial Communities Between Paddy and Natural Wetlands of Northeast China

Studies of methane-oxidizing bacteria are updating our views of their composition and function in paddy and natural wetlands. However, few studies have characterized differences in the methane-oxidizing bacterial communities between paddy and natural wetlands. Here, we conducted a ¹³C stable isotope-probing experiment and high-throughput sequencing to determine the structure profiling, co-occurrence relationships, and assembly processes of methanotrophic communities in four wetlands of Northeast China. There was a clear difference in community structure between paddy and natural wetlands. LEfSe analyses revealed that Methylobacter, FWs, and Methylosinus were enriched in natural wetlands, while Methylosarcina were prevailing in paddy, all identified as indicative methanotrophs. We observed distinct co-occurrence relationships between paddy and natural wetlands: more robust and complex connections in natural wetlands than paddy wetlands. Furthermore, the relative importance of stochastic processes was greater than that of deterministic processes, as stochastic processes explained >50% of the variation in communities. These results demonstrated that the co-occurrence relationships and assembly processes of active methanotrophic communities in paddy and natural wetlands were distinct. Overall, the results of this study enhance our understanding of the communities of methane-oxidizing bacteria in paddy and natural wetlands of Northeast China.

Phylogenetic Correlation and Symbiotic Network Explain the Interdependence Between Plants and Arbuscular Mycorrhizal Fungi in a Tibetan Alpine Meadow

Plants and arbuscular mycorrhizal fungi (AMF) can form complex symbiotic networks based on functional trait selection, contributing to the maintenance of ecosystem biodiversity and stability. However, the selectivity of host plants on AMF and the characteristics of plant-AMF networks remain unclear in Tibetan alpine meadows. In this study, we studied the AMF communities in 69 root samples from 23 plant species in a Tibetan alpine meadow using Illumina-MiSeq sequencing of the 18S rRNA gene. The results showed a significant positive correlation between the phylogenetic distances of plant species and the taxonomic dissimilarity of their AMF community. The plant-AMF network was characterized by high connectance, high nestedness, anti-modularity, and anti-specialization, and the phylogenetic signal from plants was stronger than that from AMF. The high connected and nested plant-AMF network potentially promoted the interdependence and stability of the plant-AMF symbioses in Tibetan alpine meadows. This study emphasizes that plant phylogeny and plant-AMF networks play an important role in the coevolution of host plants and their mycorrhizal partners and enhance our understanding of the interactions between aboveground and belowground communities.

Phagotrophic Protists and Their Associates: Evidence for Preferential Grazing in an Abiotically Driven Soil Ecosystem

The complex relationship between ecosystem function and soil food web structure is governed by species interactions, many of which remain unmapped. Phagotrophic protists structure soil food webs by grazing the microbiome, yet their involvement in intraguild competition, susceptibility to predator diversity, and grazing preferences are only vaguely known. These species-dependent interactions are contextualized by adjacent biotic and abiotic processes, and thus obfuscated by typically high soil biodiversity. Such questions may be investigated in the McMurdo Dry Valleys (MDV) of Antarctica because the physical environment strongly filters biodiversity and simplifies the influence of abiotic factors. To detect the potential interactions in the MDV, we analyzed the co-occurrence among shotgun metagenome sequences for associations suggestive of intraguild competition, predation, and preferential grazing. In order to control for confounding abiotic drivers, we tested co-occurrence patterns against various climatic and edaphic factors. Non-random co-occurrence between phagotrophic protists and other soil fauna was biotically driven, but we found no support for competition or predation. However, protists predominately associated with Proteobacteria and avoided Actinobacteria, suggesting grazing preferences were modulated by bacterial cell-wall structure and growth rate. Our study provides a critical starting-point for mapping protist interactions in native soils and highlights key trends for future targeted molecular and culture-based approaches.

Network Properties of Local Fungal Communities Reveal the Anthropogenic Disturbance Consequences of Farming Practices in Vineyard Soils

Agroecosystems are human-managed ecosystems subject to generalized ecological rules. Understanding the ecology behind the assembly and dynamics of soil fungal communities is a fruitful way to improve management practices and plant productivity. Thus, monitoring soil health would benefit from the use of metrics that arise from ecological explanations that can also be informative for agricultural management. Beyond traditional biodiversity descriptors, community-level properties have the potential of informing about particular ecological situations. Here we assess the impact of different farming practices in a survey of 350 vineyard soils from the United States and Spain by estimating network properties based on spatial associations. Our observations using traditional approaches show results concurring with previous literature: the influence of geographic and climatic factors on sample distributions, or different operational taxonomic unit (OTU) compositions depending on agricultural managements. Furthermore, using network properties, we observe that fungal communities ranged from dense arrangements of associations to a sparser structure of associations, indicating differential levels of niche specialization. We detect fungal arrangements capable of thriving in wider or smaller ranges of temperature, revealing that niche specialization may be a critical soil process impacting soil health. Low-intervention practices (organic and biodynamic managements) promoted densely clustered networks, describing an equilibrium state based on mixed collaborative communities. In contrast, conventionally managed vineyards had highly modular sparser communities, supported by a higher coexclusion proportion. Thus, we hypothesize that network properties at the community level may help to understand how the environment and land use can affect community structure and ecological processes in agroecosystems.

IMPORTANCE Soil fungal communities play a key role in agroecosystem sustainability. The complexity of fungal communities, at both taxonomic and functional levels, makes it difficult to find clear patterns connecting community composition with ecosystem function and to understand the impact of biotic (interspecies interactions) and abiotic (e.g., climate or anthropogenic disturbances) factors on it. Here we combine network analysis methods and properties, proposing a novel analytical approach: to infer ecological properties from local networks, which we apply to the study of fungal communities in vineyard soils. We conclude that different levels of farming intensification may lead to different ecological strategies in soil fungal communities settled by particular association arrangements.

Author Video: An author video summary of this article is available.

Distinct assembly processes shape bacterial communities along unsaturated, groundwater fluctuated, and saturated zones

The subsurface soil environment through the unsaturated (vadose) zone and saturated (below groundwater table) zone is one of the most active layers in the Earth's surface with biogeochemical interactions. Geochemical variables and geographic distance are key driving forces shaping the distribution of soil microbial communities, but our understandings are mainly limited to surface soil or shallow unsaturated zone (1-3 m beneath the ground). In this study, soil and sediment samples were collected from the unsaturated zone, through groundwater fluctuated zone, to saturated zone (up to 20 m) to unravel the assembly processes mediating vertical bacterial community succession across these three zones. Our results suggested both geochemical niches and bacterial diversity had different vertical patterns in each zone. With increased depth, pH increased and nutrient levels (C, N, P, K) and bacterial diversity declined in the unsaturated zone, and nutrients and bacterial diversity remained low levels after reaching the fluctuated and saturated zones. Nutrients were the key drivers shaping bacterial variation in the unsaturated zone, but limited nutrients and only 'depth' significantly explained the variations in the fluctuated zone and saturated zone, respectively. The co-occurrence network supported a more species co-existence pattern in the unsaturated zone than that in the other two zones. Due to the geochemical variations across three zones, the assembly of phylogenetically more clustered communities was observed through deterministic processes (e.g., 55% homogenizing selection) in the unsaturated zone, but the stochastic process (e.g., 50%-70% dispersal limitation) was more important in the fluctuated and saturated zones. These findings together suggested that the vertical distribution of soil bacterial community assembly was zone-specific and shaped by the degree of deterministic vs. stochastic processes. Our results provide a novel insight into the microbial community assembly across three different ecosystems in the Earth's critical zone and shed a light on subsurface biogeochemical processes.

nir gene-based co-occurrence patterns reveal assembly mechanisms of soil denitrifiers in response to fire

Denitrification causes nitrogen losses from terrestrial ecosystems. The magnitude of nitrogen loss depends on the prevalence of denitrifiers, which show ecological differences if they harbour nirS or nirK genes encoding nitrite reductases with the same biological function. Thus, it is relevant to understand the mechanisms of co-existence of denitrifiers, including their response to environmental filters and competition due to niche similarities. We propose a framework to analyse the co-existence of denitrifiers across multiple assemblages by using nir gene-based co-occurrence networks. We applied it in Mediterranean soils before and during 1 year after an experimental fire. Burning did not modify nir community structure, but significantly impacted co-occurrence patterns. Bacteria with the same nir co-occurred in space, and those with different nir excluded each other, reflecting niche requirements: nirS abundance responded to nitrate and salinity, whereas nirK to iron content. Prior to fire, mutual exclusion between bacteria with the same nir suggested competition due to niche similarities. Burning provoked an immediate rise in mineral nitrogen and erased the signals of competition, which emerged again within days as nir abundances peaked. nir co-occurrence patterns can help infer the assembly mechanisms of denitrifying communities, which control nitrogen losses in the face of ecological disturbance.