Difficulty in inferring microbial community structure based on co-occurrence network approaches

Microbiome modeling: a beginner's guide

Microbiomes, comprised of diverse microbial species and viruses, play pivotal roles in human health, environmental processes, and biotechnological applications and interact with each other, their environment, and hosts via ecological interactions. Our understanding of microbiomes is still limited and hampered by their complexity. A concept improving this understanding is systems biology, which focuses on the holistic description of biological systems utilizing experimental and computational methods. An important set of such experimental methods are metaomics methods which analyze microbiomes and output lists of molecular features. These lists of data are integrated, interpreted, and compiled into computational microbiome models, to predict, optimize, and control microbiome behavior. There exists a gap in understanding between microbiologists and modelers/bioinformaticians, stemming from a lack of interdisciplinary knowledge. This knowledge gap hinders the establishment of computational models in microbiome analysis. This review aims to bridge this gap and is tailored for microbiologists, researchers new to microbiome modeling, and bioinformaticians. To achieve this goal, it provides an interdisciplinary overview of microbiome modeling, starting with fundamental knowledge of microbiomes, metaomics methods, common modeling formalisms, and how models facilitate microbiome control. It concludes with guidelines and repositories for modeling. Each section provides entry-level information, example applications, and important references, serving as a valuable resource for comprehending and navigating the complex landscape of microbiome research and modeling.

Integration of metataxonomic data sets into microbial association networks highlights shared bacterial community dynamics in fermented vegetables

The management of food fermentation is still largely based on empirical knowledge, as the dynamics of microbial communities and the underlying metabolic networks that produce safe and nutritious products remain beyond our understanding. Although these closed ecosystems contain relatively few taxa, they have not yet been thoroughly characterized with respect to how their microbial communities interact and dynamically evolve. However, with the increased availability of metataxonomic data sets on different fermented vegetables, it is now possible to gain a comprehensive understanding of the microbial relationships that structure plant fermentation. In this study, we applied a network-based approach to the integration of public metataxonomic 16S data sets targeting different fermented vegetables throughout time. Specifically, we aimed to explore, compare, and combine public 16S data sets to identify shared associations between amplicon sequence variants (ASVs) obtained from independent studies. The workflow includes steps for searching and selecting public time-series data sets and constructing association networks of ASVs based on co-abundance metrics. Networks for individual data sets are then integrated into a core network, highlighting significant associations. Microbial communities are identified based on the comparison and clustering of ASV networks using the "stochastic block model" method. When we applied this method to 10 public data sets (including a total of 931 samples) targeting five varieties of vegetables with different sampling times, we found that it was able to shed light on the dynamics of vegetable fermentation by characterizing the processes of community succession among different bacterial assemblages.

IMPORTANCE

Within the growing body of research on the bacterial communities involved in the fermentation of vegetables, there is particular interest in discovering the species or consortia that drive different fermentation steps. This integrative analysis demonstrates that the reuse and integration of public microbiome data sets can provide new insights into a little-known biotope. Our most important finding is the recurrent but transient appearance, at the beginning of vegetable fermentation, of amplicon sequence variants (ASVs) belonging to Enterobacterales and their associations with ASVs belonging to Lactobacillales. These findings could be applied to the design of new fermented products.

Sediment DNA Records the Critical Transition of Bacterial Communities in the Arid Lake

It is necessary to predict the critical transition of lake ecosystems due to their abrupt, non-linear effects on social-economic systems. Given the promising application of paleolimnological archives to tracking the historical changes of lake ecosystems, it is speculated that they can also record the lake's critical transition. We studied Lake Dali-Nor in the arid region of Inner Mongolia because of the profound shrinking the lake experienced between the 1300 s and the 1600 s. We reconstructed the succession of bacterial communities from a 140-cm-long sediment core at 4-cm intervals and detected the critical transition. Our results showed that the historical trajectory of bacterial communities from the 1200 s to the 2010s was divided into two alternative states: state1 from 1200 to 1300 s and state2 from 1400 to 2010s. Furthermore, in the late 1300 s, the appearance of a tipping point and critical slowing down implied the existence of a critical transition. By using a multi-decadal time series from the sedimentary core, with general Lotka-Volterra model simulations, local stability analysis found that bacterial communities were the most unstable as they approached the critical transition, suggesting that the collapse of stability triggers the community shift from an equilibrium state to another state. Furthermore, the most unstable community harbored the strongest antagonistic and mutualistic interactions, which may imply the detrimental role of interaction strength on community stability. Collectively, our study showed that sediment DNA can be used to detect the critical transition of lake ecosystems.

Non-mercury methylating microbial taxa are integral to understanding links between mercury methylation and elemental cycles in marine and freshwater sediments

The goal of this study was to explore the role of non-mercury (Hg) methylating taxa in mercury methylation and to identify potential links between elemental cycles and Hg methylation. Statistical approaches were utilized to investigate the microbial community and biochemical functions in relation to methylmercury (MeHg) concentrations in marine and freshwater sediments. Sediments were collected from the methylation zone (top 15 cm) in four Hg-contaminated sites. Both abiotic (e.g., sulfate, sulfide, iron, salinity, total organic matter, etc.) and biotic factors (e.g., hgcA, abundances of methylating and non-methylating taxa) were quantified. Random forest and stepwise regression were performed to assess whether non-methylating taxa were significantly associated with MeHg concentration. Co-occurrence and functional network analyses were constructed to explore associations between taxa by examining microbial community structure, composition, and biochemical functions across sites. Regression analysis showed that approximately 80% of the variability in sediment MeHg concentration was predicted by total mercury concentration, the abundances of Hg methylating taxa, and the abundances of the non-Hg methylating taxa. The co-occurrence networks identified Paludibacteraceae and Syntrophorhabdaceae as keystone non Hg methylating taxa in multiple sites, indicating the potential for syntrophic interactions with Hg methylators. Strong associations were also observed between methanogens and sulfate-reducing bacteria, which were likely symbiotic associations. The functional network results suggested that non-Hg methylating taxa play important roles in sulfur respiration, nitrogen respiration, and the carbon metabolism-related functions methylotrophy, methanotrophy, and chemoheterotrophy. Interestingly, keystone functions varied by site and did not involve carbon- and sulfur-related functions only. Our findings highlight associations between methylating and non-methylating taxa and sulfur, carbon, and nitrogen cycles in sediment methylation zones, with implications for predicting and understanding the impact of climate and land/sea use changes on Hg methylation.

Organic farming promotes the abundance of fungi keystone taxa in bacteria-fungi interkingdom networks

Soil bacteria-fungi interactions are essential in the biogeochemical cycles of several nutrients, making these microbes major players in agroecosystems. While the impact of the farming system on microbial community composition has been extensively reported in the literature, whether sustainable farming approaches can promote associations between bacteria and fungi is still unclear. To study this, we employed 16S, ITS, and 18S DNA sequencing to uncover how microbial interactions were affected by conventional and organic farming systems on maize crops. The Bray-Curtis index revealed that bacterial, fungal, and arbuscular mycorrhizal fungi communities were significantly different between the two farming systems. Several taxa known to thrive in healthy soils, such as Nitrosophaerales, Orbiliales, and Glomus were more abundant in the organic farming system. Constrained ordination revealed that the organic farming system microbial community was significantly correlated with the β-glucosidase activity, whereas the conventional farming system microbial community significantly correlated with soil pH. Both conventional and organic co-occurrence interkingdom networks exhibited a parallel node count, however, the former had a higher number of edges, thus being denser than the latter. Despite the similar amount of fungal nodes in the co-occurrence networks, the organic farming system co-occurrence network exhibited more than 3-fold the proportion of fungal taxa as keystone nodes than the conventional co-occurrence network. The genera Bionectria, Cercophora, Geastrum, Penicillium, Preussia, Metarhizium, Myceliophthora, and Rhizophlyctis were among the fungal keystone nodes of the organic farming system network. Altogether, our results uncover that beyond differences in microbial community composition between the two farming systems, fungal keystone nodes are far more relevant in the organic farming system, thus suggesting that bacteria-fungi interactions are more frequent in organic farming systems, promoting a more functional microbial community.

Engineering agricultural soil microbiomes and predicting plant phenotypes

Plant growth-promoting rhizobacteria (PGPR) can improve crop yields, nutrient use efficiency, plant tolerance to stressors, and confer benefits to future generations of crops grown in the same soil. Unlocking the potential of microbial communities in the rhizosphere and endosphere is therefore of great interest for sustainable agriculture advancements. Before plant microbiomes can be engineered to confer desirable phenotypic effects on their plant hosts, a deeper understanding of the interacting factors influencing rhizosphere community structure and function is needed. Dealing with this complexity is becoming more feasible using computational approaches. In this review, we discuss recent advances at the intersection of experimental and computational strategies for the investigation of plant-microbiome interactions and the engineering of desirable soil microbiomes.

Stochastic Assembly Increases the Complexity and Stability of Shrimp Gut Microbiota During Aquaculture Progression

The gut microbiota of aquaculture species contributes to their food metabolism and regulates their health, which has been shown to vary during aquaculture progression of their hosts. However, limited research has examined the outcomes and mechanisms of these changes in the gut microbiota of hosts. Here, Kuruma shrimps from the beginning, middle, and late stages of aquaculture progression (about a time duration of 2 months between each stage) were collected and variations in the gut microbiota of Kuruma shrimp during the whole aquaculture process were examined. High-throughput sequencing demonstrated increases in the diversity and richness of the shrimp gut microbiota with aquaculture progression. In addition, the gut microbiota composition differed among cultural stages, with enrichment of Firmicutes, RF39, and Megamonas and a reduction in Proteobacteria in the mid-stage. Notably, only very few taxa were persistent in the shrimp gut microbiota during the whole aquaculture progression, while the number of taxa that specific to the end of aquaculture was high. Network analysis revealed increasing complexity of the shrimp gut microbiota during aquaculture progression. Moreover, the shrimp gut microbiota became significantly more stable towards the end of aquaculture. According to the results of neutral community model, contribution of stochastic processes for shaping the shrimp gut microbiota was elevated along the aquaculture progression. This study showed substantial variations in shrimp gut microbiota during aquaculture progression and explored the underlying mechanisms regulating these changes.

Viroid-like colonists of human microbiomes

Here, we describe the "Obelisks," a previously unrecognised class of viroid-like elements that we first identified in human gut metatranscriptomic data. "Obelisks" share several properties: (i) apparently circular RNA ~1kb genome assemblies, (ii) predicted rod-like secondary structures encompassing the entire genome, and (iii) open reading frames coding for a novel protein superfamily, which we call the "Oblins". We find that Obelisks form their own distinct phylogenetic group with no detectable sequence or structural similarity to known biological agents. Further, Obelisks are prevalent in tested human microbiome metatranscriptomes with representatives detected in ~7% of analysed stool metatranscriptomes (29/440) and in ~50% of analysed oral metatranscriptomes (17/32). Obelisk compositions appear to differ between the anatomic sites and are capable of persisting in individuals, with continued presence over >300 days observed in one case. Large scale searches identified 29,959 Obelisks (clustered at 90% nucleotide identity), with examples from all seven continents and in diverse ecological niches. From this search, a subset of Obelisks are identified to code for Obelisk-specific variants of the hammerhead type-III self-cleaving ribozyme. Lastly, we identified one case of a bacterial species (Streptococcus sanguinis) in which a subset of defined laboratory strains harboured a specific Obelisk RNA population. As such, Obelisks comprise a class of diverse RNAs that have colonised, and gone unnoticed in, human, and global microbiomes.

Anthropic disturbances impact the soil microbial network structure and stability to a greater extent than natural disturbances in an arid ecosystem

Growing pressure from climate change and agricultural land use is destabilizing soil microbial community interactions. Yet little is known about microbial community resistance and adaptation to disturbances over time. This hampers our ability to determine the recovery latency of microbial interactions after disturbances, with fundamental implications for ecosystem functioning and conservation measures. Here we examined the response of bacterial and fungal community networks in the rhizosphere of Haloxylon salicornicum (Moq.) Bunge ex Boiss. over the course of soil disturbances resulting from a history of different hydric constraints involving flooding-drought successions. An anthropic disturbance related to past agricultural use, with frequent successions of flooding and drought, was compared to a natural disturbance, i.e., an evaporation basin, with yearly flooding-drought successions. The anthropic disturbance resulted in a specific microbial network topology characterized by lower modularity and stability, reflecting the legacy of past agricultural use on soil microbiome. In contrast, the natural disturbance resulted in a network topology and stability close to those of natural environments despite the lower alpha diversity, and a different community composition compared to that of the other sites. These results highlighted the temporality in the response of the microbial community structure to disturbance, where long-term adaptation to flooding-drought successions lead to a higher stability than disturbances occurring over a shorter timescale.

Succession and assembly mechanisms of seawater prokaryotic communities along an extremely wide salinity gradient

Salinity is an important environmental factor in microbial ecology for affecting the microbial communities in diverse environments. Understanding the salinity adaptation mechanisms of a microbial community is a significant issue, while most previous studies only covered a narrow salinity range. Here, variations in seawater prokaryotic communities during the whole salt drying progression (salinity from 3% to 25%) were investigated. According to high-throughput sequencing results, the diversity, composition, and function of seawater prokaryotic communities varied significantly along the salinity gradient, expressing as decreased diversity, enrichment of some halophilic archaea, and powerful nitrate reduction in samples with high salt concentrations. More importantly, a sudden and dramatic alteration of prokaryotic communities was observed when salinity reached 16%, which was recognized as the change point. Combined with the results of network analysis, we found the increasing of complexity but decreasing of stability in prokaryotic communities when salinity exceeded the change point. Moreover, prokaryotic communities became more deterministic when salinity exceeded the change point due to the niche adaptation of halophilic species. Our study showed that substantial variations in seawater prokaryotic communities along an extremely wide salinity gradient, and also explored the underlying mechanisms regulating these changes.

Community assembly of the human piercing microbiome

Predicting how biological communities respond to disturbance requires understanding the forces that govern their assembly. We propose using human skin piercings as a model system for studying community assembly after rapid environmental change. Local skin sterilization provides a 'clean slate' within the novel ecological niche created by the piercing. Stochastic assembly processes can dominate skin microbiomes due to the influence of environmental exposure on local dispersal, but deterministic processes might play a greater role within occluded skin piercings if piercing habitats impose strong selection pressures on colonizing species. Here we explore the human ear-piercing microbiome and demonstrate that community assembly is predominantly stochastic but becomes significantly more deterministic with time, producing increasingly diverse and ecologically complex communities. We also observed changes in two dominant and medically relevant antagonists (Cutibacterium acnes and Staphylococcus epidermidis), consistent with competitive exclusion induced by a transition from sebaceous to moist environments. By exploiting this common yet uniquely human practice, we show that skin piercings are not just culturally significant but also represent ecosystem engineering on the human body. The novel habitats and communities that skin piercings produce may provide general insights into biological responses to environmental disturbances with implications for both ecosystem and human health.

Network analysis of patterns and relevance of enteric pathogen co-infections among infants in a diarrhea-endemic setting

Despite significant progress in recent decades toward ameliorating the excess burden of diarrheal disease globally, childhood diarrhea remains a leading cause of morbidity and mortality in low-and-middle-income countries (LMICs). Recent large-scale studies of diarrhea etiology in these populations have revealed widespread co-infection with multiple enteric pathogens, in both acute and asymptomatic stool specimens. We applied methods from network science and ecology to better understand the underlying structure of enteric co-infection among infants in two large longitudinal birth cohorts in Bangladesh. We used a configuration model to establish distributions of expected random co-occurrence, based on individual pathogen prevalence alone, for every pathogen pair among 30 enteropathogens detected by qRT-PCR in both diarrheal and asymptomatic stool specimens. We found two pairs, Enterotoxigenic E. coli (ETEC) with Enteropathogenic E. coli (EPEC), and ETEC with Campylobacter spp., co-infected significantly more than expected at random (both pairs co-occurring almost 4 standard deviations above what one could expect due to chance alone). Furthermore, we found a general pattern that bacteria-bacteria pairs appear together more frequently than expected at random, while virus-bacteria pairs tend to appear less frequently than expected based on model predictions. Finally, infants co-infected with leading bacteria-bacteria pairs had more days of diarrhea in the first year of life compared to infants without co-infection (p-value <0.0001). Our methods and results help us understand the structure of enteric co-infection which can guide further work to identify and eliminate common sources of infection or determine biologic mechanisms that promote co-infection.

More Than the Sum of Its Parts: Unlocking the Power of Network Structure for Understanding Organization and Function in Microbiomes

Plant and soil microbiomes are integral to the health and productivity of plants and ecosystems, yet researchers struggle to identify microbiome characteristics important for providing beneficial outcomes. Network analysis offers a shift in analytical framework beyond "who is present" to the organization or patterns of coexistence between microbes within the microbiome. Because microbial phenotypes are often significantly impacted by coexisting populations, patterns of coexistence within microbiomes are likely to be especially important in predicting functional outcomes. Here, we provide an overview of the how and why of network analysis in microbiome research, highlighting the ways in which network analyses have provided novel insights into microbiome organization and functional capacities, the diverse network roles of different microbial populations, and the eco-evolutionary dynamics of plant and soil microbiomes.

Inferring microbial co-occurrence networks from amplicon data: a systematic evaluation

Microbes commonly organize into communities consisting of hundreds of species involved in complex interactions with each other. 16S ribosomal RNA (16S rRNA) amplicon profiling provides snapshots that reveal the phylogenies and abundance profiles of these microbial communities. These snapshots, when collected from multiple samples, can reveal the co-occurrence of microbes, providing a glimpse into the network of associations in these communities. However, the inference of networks from 16S data involves numerous steps, each requiring specific tools and parameter choices. Moreover, the extent to which these steps affect the final network is still unclear. In this study, we perform a meticulous analysis of each step of a pipeline that can convert 16S sequencing data into a network of microbial associations. Through this process, we map how different choices of algorithms and parameters affect the co-occurrence network and identify the steps that contribute substantially to the variance. We further determine the tools and parameters that generate robust co-occurrence networks and develop consensus network algorithms based on benchmarks with mock and synthetic data sets. The Microbial Co-occurrence Network Explorer, or MiCoNE (available at https://github.com/segrelab/MiCoNE) follows these default tools and parameters and can help explore the outcome of these combinations of choices on the inferred networks. We envisage that this pipeline could be used for integrating multiple data sets and generating comparative analyses and consensus networks that can guide our understanding of microbial community assembly in different biomes. IMPORTANCE Mapping the interrelationships between different species in a microbial community is important for understanding and controlling their structure and function. The surge in the high-throughput sequencing of microbial communities has led to the creation of thousands of data sets containing information about microbial abundances. These abundances can be transformed into co-occurrence networks, providing a glimpse into the associations within microbiomes. However, processing these data sets to obtain co-occurrence information relies on several complex steps, each of which involves numerous choices of tools and corresponding parameters. These multiple options pose questions about the robustness and uniqueness of the inferred networks. In this study, we address this workflow and provide a systematic analysis of how these choices of tools affect the final network and guidelines on appropriate tool selection for a particular data set. We also develop a consensus network algorithm that helps generate more robust co-occurrence networks based on benchmark synthetic data sets.

Exploring the Frequency and Distribution of Ecological Non-monotonicity in Associations among Ecosystem Constituents

Complex links between biotic and abiotic constituents are fundamental for the functioning of ecosystems. Although non-monotonic interactions and associations are known to increase the stability, diversity, and productivity of ecosystems, they are frequently ignored by community-level standard statistical approaches. Using the copula-based dependence measure qad, capable of quantifying the directed and asymmetric dependence between variables for all forms of (functional) relationships, we determined the proportion of non-monotonic associations between different constituents of an ecosystem (plants, bacteria, fungi, and environmental parameters). Here, we show that up to 59% of all statistically significant associations are non-monotonic. Further, we show that pairwise associations between plants, bacteria, fungi, and environmental parameters are specifically characterized by their strength and degree of monotonicity, for example, microbe-microbe associations are on average stronger than and differ in degree of non-monotonicity from plant-microbe associations. Considering directed and non-monotonic associations, we extended the concept of ecosystem coupling providing more complete insights into the internal order of ecosystems. Our results emphasize the importance of ecological non-monotonicity in characterizing and understanding ecosystem patterns and processes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10021-023-00867-9.

AM fungal-bacterial relationships: what can they tell us about ecosystem sustainability and soil functioning?

Considering our growing population and our continuous degradation of soil environments, understanding the fundamental ecology of soil biota and plant microbiomes will be imperative to sustaining soil systems. Arbuscular mycorrhizal (AM) fungi extend their hyphae beyond plant root zones, creating microhabitats with bacterial symbionts for nutrient acquisition through a tripartite symbiotic relationship along with plants. Nonetheless, it is unclear what drives these AM fungal-bacterial relationships and how AM fungal functional traits contribute to these relationships. By delving into the literature, we look at the drivers and complexity behind AM fungal-bacterial relationships, describe the shift needed in AM fungal research towards the inclusion of interdisciplinary tools, and discuss the utilization of bacterial datasets to provide contextual evidence behind these complex relationships, bringing insights and new hypotheses to AM fungal functional traits. From this synthesis, we gather that interdependent microbial relationships are at the foundation of understanding microbiome functionality and deciphering microbial functional traits. We suggest using pattern-based inference tools along with machine learning to elucidate AM fungal-bacterial relationship trends, along with the utilization of synthetic communities, functional gene analyses, and metabolomics to understand how AM fungal and bacterial communities facilitate communication for the survival of host plant communities. These suggestions could result in improving microbial inocula and products, as well as a better understanding of complex relationships in terrestrial ecosystems that contribute to plant-soil feedbacks.

Ecosystem multifunctionality and soil microbial communities in response to ecological restoration in an alpine degraded grassland

Linkages between microbial communities and multiple ecosystem functions are context-dependent. However, the impacts of different restoration measures on microbial communities and ecosystem functioning remain unclear. Here, a 14-year long-term experiment was conducted using three restoration modes: planting mixed grasses (MG), planting shrub with Salix cupularis alone (SA), and planting shrub with Salix cupularis plus planting mixed grasses (SG), with an extremely degraded grassland serving as the control (CK). Our objective was to investigate how ecosystem multifunctionality and microbial communities (diversity, composition, and co-occurrence networks) respond to different restoration modes. Our results indicated that most of individual functions (i.e., soil nutrient contents, enzyme activities, and microbial biomass) in the SG treatment were significantly higher than in the CK treatment, and even higher than MG and SA treatments. Compared with the CK treatment, treatments MG, SA, and SG significantly increased the multifunctionality index on average by 0.57, 0.23 and 0.76, respectively. Random forest modeling showed that the alpha-diversity and composition of bacterial communities, rather than fungal communities, drove the ecosystem multifunctionality. Moreover, we found that both the MG and SG treatments significantly improved bacterial network stability, which exhabited stronger correlations with ecosystem multifunctionality compared to fungal network stability. In summary, this study demonstrates that planting shrub and grasses altogether is a promising restoration mode that can enhance ecosystem multifunctionality and improve microbial diversity and stability in the alpine degraded grassland.

Bakdrive: identifying a minimum set of bacterial species driving interactions across multiple microbial communities

MOTIVATION

Interactions among microbes within microbial communities have been shown to play crucial roles in human health. In spite of recent progress, low-level knowledge of bacteria driving microbial interactions within microbiomes remains unknown, limiting our ability to fully decipher and control microbial communities.

RESULTS

We present a novel approach for identifying species driving interactions within microbiomes. Bakdrive infers ecological networks of given metagenomic sequencing samples and identifies minimum sets of driver species (MDS) using control theory. Bakdrive has three key innovations in this space: (i) it leverages inherent information from metagenomic sequencing samples to identify driver species, (ii) it explicitly takes host-specific variation into consideration, and (iii) it does not require a known ecological network. In extensive simulated data, we demonstrate identifying driver species identified from healthy donor samples and introducing them to the disease samples, we can restore the gut microbiome in recurrent Clostridioides difficile (rCDI) infection patients to a healthy state. We also applied Bakdrive to two real datasets, rCDI and Crohn's disease patients, uncovering driver species consistent with previous work. Bakdrive represents a novel approach for capturing microbial interactions.

AVAILABILITY AND IMPLEMENTATION

Bakdrive is open-source and available at: https://gitlab.com/treangenlab/bakdrive.

Deciphering Interactions Within a 4-Strain Riverine Bacterial Community

The dynamics of a community of four planktonic bacterial strains isolated from river water was followed in R2 broth for 72 h in batch experiments. These strains were identified as Janthinobacterium sp., Brevundimonas sp., Flavobacterium sp. and Variovorax sp. 16S rRNA gene sequencing and flow cytometry analyses were combined to monitor the change in abundance of each individual strain in bi-cultures and quadri-culture. Two interaction networks were constructed that summarize the impact of the strains on each other's growth rate in exponential phase and carrying capacity in stationary phase. The networks agree on the absence of positive interactions but also show differences, implying that ecological interactions can be specific to particular growth phases. Janthinobacterium sp. was the fastest growing strain and dominated the co-cultures. However, its growth rate was negatively affected by the presence of other strains 10 to 100 times less abundant than Janthinobacterium sp. In general, we saw a positive correlation between growth rate and carrying capacity in this system. In addition, growth rate in monoculture was predictive of carrying capacity in co-culture. Taken together, our results highlight the necessity to take growth phases into account when measuring interactions within a microbial community. In addition, evidence that a minor strain can greatly influence the dynamics of a dominant one underlines the necessity to choose population models that do not assume a linear dependency of interaction strength to abundance of other species for accurate parameterization from such empirical data.

Wavelet clustering analysis as a tool for characterizing community structure in the human microbiome

Human microbiome research is helped by the characterization of microbial networks, as these may reveal key microbes that can be targeted for beneficial health effects. Prevailing methods of microbial network characterization are based on measures of association, often applied to limited sampling points in time. Here, we demonstrate the potential of wavelet clustering, a technique that clusters time series based on similarities in their spectral characteristics. We illustrate this technique with synthetic time series and apply wavelet clustering to densely sampled human gut microbiome time series. We compare our results with hierarchical clustering based on temporal correlations in abundance, within and across individuals, and show that the cluster trees obtained by using either method are significantly different in terms of elements clustered together, branching structure and total branch length. By capitalizing on the dynamic nature of the human microbiome, wavelet clustering reveals community structures that remain obscured in correlation-based methods.

Capturing the dynamics of microbial interactions through individual-specific networks

Longitudinal analysis of multivariate individual-specific microbiome profiles over time or across conditions remains dauntin. Most statistical tools and methods that are available to study microbiomes are based on cross-sectional data. Over the past few years, several attempts have been made to model the dynamics of bacterial species over time or across conditions. However, the field needs novel views on handling microbial interactions in temporal analyses. This study proposes a novel data analysis framework, MNDA, that combines representation learning and individual-specific microbial co-occurrence networks to uncover taxon neighborhood dynamics. As a use case, we consider a cohort of newborns with microbiomes available at 6 and 9 months after birth, and extraneous data available on the mode of delivery and diet changes between the considered time points. Our results show that prediction models for these extraneous outcomes based on an MNDA measure of local neighborhood dynamics for each taxon outperform traditional prediction models solely based on individual-specific microbial abundances. Furthermore, our results show that unsupervised similarity analysis of newborns in the study, again using the notion of a taxon's dynamic neighborhood derived from time-matched individual-specific microbial networks, can reveal different subpopulations of individuals, compared to standard microbiome-based clustering, with potential relevance to clinical practice. This study highlights the complementarity of microbial interactions and abundances in downstream analyses and opens new avenues to personalized prediction or stratified medicine with temporal microbiome data.

Disentangling temporal associations in marine microbial networks

BACKGROUND

Microbial interactions are fundamental for Earth's ecosystem functioning and biogeochemical cycling. Nevertheless, they are challenging to identify and remain barely known. Omics-based censuses are helpful in predicting microbial interactions through the statistical inference of single (static) association networks. Yet, microbial interactions are dynamic and we have limited knowledge of how they change over time. Here, we investigate the dynamics of microbial associations in a 10-year marine time series in the Mediterranean Sea using an approach inferring a time-resolved (temporal) network from a single static network.

RESULTS

A single static network including microbial eukaryotes and bacteria was built using metabarcoding data derived from 120 monthly samples. For the decade, we aimed to identify persistent, seasonal, and temporary microbial associations by determining a temporal network that captures the interactome of each individual sample. We found that the temporal network appears to follow an annual cycle, collapsing, and reassembling when transiting between colder and warmer waters. We observed higher association repeatability in colder than in warmer months. Only 16 associations could be validated using observations reported in literature, underlining our knowledge gap in marine microbial ecological interactions.

CONCLUSIONS

Our results indicate that marine microbial associations follow recurrent temporal dynamics in temperate zones, which need to be accounted for to better understand the functioning of the ocean microbiome. The constructed marine temporal network may serve as a resource for testing season-specific microbial interaction hypotheses. The applied approach can be transferred to microbiome studies in other ecosystems. Video Abstract.

A Novel E3 Probiotics Formula Restored Gut Dysbiosis and Remodelled Gut Microbial Network and Microbiome Dysbiosis Index (MDI) in Southern Chinese Adult Psoriasis Patients

Psoriasis is a common chronic immune-mediated inflammatory skin disease with the association of various comorbidities. Despite the introduction of highly effective biologic therapies over the past few decades, the exact trigger for an immune reaction in psoriasis is unclear. With the majority of immune cells residing in the gut, the effect of gut microbiome dysbiosis goes beyond the gastrointestinal site and may exacerbate inflammation and regulate the immune system elsewhere, including but not limited to the skin via the gut-skin axis. In order to delineate the role of the gut microbiome in Southern Chinese psoriasis patients, we performed targeted 16S rRNA sequencing and comprehensive bioinformatic analysis to compare the gut microbiome profile of 58 psoriasis patients against 49 healthy local subjects presumably with similar lifestyles. Blautia wexlerae and Parabacteroides distasonis were found to be enriched in psoriasis patients and in some of the healthy subjects, respectively. Metabolic functional pathways were predicted to be differentially abundant, with a clear shift toward SCFA synthesis in healthy subjects. The alteration of the co-occurrence network was also evident in the psoriasis group. In addition, we also profiled the gut microbiome in 52 of the 58 recruited psoriasis patients after taking 8 weeks of an orally administrated novel E3 probiotics formula (with prebiotics, probiotics and postbiotics). The Dermatological Life Quality Index (p = 0.009) and Psoriasis Area and Severity Index (p < 0.001) were significantly improved after taking 8 weeks of probiotics with no adverse effect observed. We showed that probiotics could at least partly restore gut dysbiosis via the modulation of the gut microbiome. Here, we also report the potential application of a machine learning-derived gut dysbiosis index based on a quantitative PCR panel (AUC = 0.88) to monitor gut dysbiosis in psoriasis patients. To sum up, our study suggests the gut microbial landscape differed in psoriasis patients at the genera, species, functional and network levels. Additionally, the dysbiosis index could be a cost-effective and rapid tool to monitor probiotics use in psoriasis patients.

Interactions between Culturable Bacteria Are Predicted by Individual Species' Growth

Predicting interspecies interactions is a key challenge in microbial ecology given that interactions shape the composition and functioning of microbial communities. However, predicting microbial interactions is challenging because they can vary considerably depending on species' metabolic capabilities and environmental conditions. Here, we employ machine learning models to predict pairwise interactions between culturable bacteria based on their phylogeny, monoculture growth capabilities, and interactions with other species. We trained our models on one of the largest available pairwise interactions data set containing over 7,500 interactions between 20 species from two taxonomic groups that were cocultured in 40 different carbon environments. Our models accurately predicted both the sign (accuracy of 88%) and the strength of effects (R² of 0.87) species had on each other's growth. Encouragingly, predictions with comparable accuracy could be made even when not relying on information about interactions with other species, which are often hard to measure. However, species' monoculture growth was essential to the model, as predictions based solely on species' phylogeny and inferred metabolic capabilities were significantly less accurate. These results bring us one step closer to a predictive understanding of microbial communities, which is essential for engineering beneficial microbial consortia. IMPORTANCE In order to understand the function and structure of microbial communities, one must know all pairwise interactions that occur between the different species within the community, as these interactions shape the community's structure and functioning. However, measuring all pairwise interactions can be an extremely difficult task especially when dealing with big complex communities. Because of that, predicting interspecies interactions is a key challenge in microbial ecology. Here, we use machine learning models in order to accurately predict the type and strength of interactions. We trained our models on one of the largest available pairwise interactions data set, containing over 7,500 interactions between 20 different species that were cocultured in 40 different environments. Our results show that, in general, accurate predictions can be made, and that the ability of each species to grow on its own in the given environment contributes the most to predictions. Being able to predict microbial interactions would put us one step closer to predicting the functionality of microbial communities and to rationally microbiome engineering.

A guide for comparing microbial co-occurrence networks

The article provides a pipeline for comparing microbial co-occurrence networks based on the R microeco package and meconetcomp package. It has high flexibility and expansibility and can help users efficiently compare networks built from different groups of samples or different construction approaches.

Haptophyte communities along the Kuroshio current reveal their geographical sources and ecological traits

Haptophytes are one of the most ecologically successful phytoplankton groups in the modern ocean and tend to maintain balanced and stable communities across various environments. However, little is known about the mechanisms that enable community stability and ecological success. To reveal the community characteristics and interactions among haptophytes, we conducted comprehensive observations from the upstream to downstream regions of the Kuroshio Current. Haptophyte abundance and taxonomy were assessed using quantitative polymerase chain reaction and metabarcoding of 18S rRNA sequences, respectively. The haptophyte community structure changed abruptly at sites on the shelf-slope of the East China Sea, indicating the strong influence of shelf waters with high phytoplankton biomass on downstream communities. Correlation network analysis combined with the phylogeny suggested that haptophytes can coexist with their close relatives, possibly owing to their nutritional flexibility, thereby escaping from resource competition. Consistently, some noncalcifying haptophyte genera with high mixotrophic capacities such as Chrysochromulina constituted a major component of the co-occurrence network, whereas coccolithophores such as Emiliania/Gephyrocapsa were rarely observed. Our study findings suggest that noncalcifying haptophytes play crucial roles in community diversity and stability, and in sustaining the food web structure in the Kuroshio ecosystems.

Harnessing the Power of Model Organisms To Unravel Microbial Functions in the Coral Holobiont

Stony corals build the framework of coral reefs, ecosystems of immense ecological and economic importance. The existence of these ecosystems is threatened by climate change and other anthropogenic stressors that manifest in microbial dysbiosis such as coral bleaching and disease, often leading to coral mortality. Despite a significant amount of research, the mechanisms ultimately underlying these destructive phenomena, and what could prevent or mitigate them, remain to be resolved. This is mostly due to practical challenges in experimentation on corals and the highly complex nature of the coral holobiont that also includes bacteria, archaea, protists, and viruses. While the overall importance of these partners is well recognized, their specific contributions to holobiont functioning and their interspecific dynamics remain largely unexplored. Here, we review the potential of adopting model organisms as more tractable systems to address these knowledge gaps. We draw on parallels from the broader biological and biomedical fields to guide the establishment, implementation, and integration of new and emerging model organisms with the aim of addressing the specific needs of coral research. We evaluate the cnidarian models Hydra, Aiptasia, Cassiopea, and Astrangia poculata; review the fast-evolving field of coral tissue and cell cultures; and propose a framework for the establishment of "true" tropical reef-building coral models. Based on this assessment, we also suggest future research to address key aspects limiting our ability to understand and hence improve the response of reef-building corals to future ocean conditions.

Environmental and Anthropogenic Factors Shape the Skin Bacterial Communities of a Semi-Arid Amphibian Species

The amphibian skin microbiome is important in maintaining host health, but is vulnerable to perturbation from changes in biotic and abiotic conditions. Anthropogenic habitat disturbance and emerging infectious diseases are both potential disrupters of the skin microbiome, in addition to being major drivers of amphibian decline globally. We investigated how host environment (hydrology, habitat disturbance), pathogen presence, and host biology (life stage) impact the skin microbiome of wild Dhofar toads (Duttaphrynus dhufarensis) in Oman. We detected ranavirus (but not Batrachochytrium dendrobatidis) across all sampling sites, constituting the first report of this pathogen in Oman, with reduced prevalence in disturbed sites. We show that skin microbiome beta diversity is driven by host life stage, water source, and habitat disturbance, but not ranavirus infection. Finally, although trends in bacterial diversity and differential abundance were evident in disturbed versus undisturbed sites, bacterial co-occurrence patterns determined through network analyses revealed high site specificity. Our results therefore provide support for amphibian skin microbiome diversity and taxa abundance being associated with habitat disturbance, with bacterial co-occurrence (and likely broader aspects of microbial community ecology) being largely site specific.

Tumour microbiomes and Fusobacterium genomics in Vietnamese colorectal cancer patients

Perturbations in the gut microbiome have been associated with colorectal cancer (CRC), with the colonic overabundance of Fusobacterium nucleatum shown as the most consistent marker. Despite its significance in the promotion of CRC, genomic studies of Fusobacterium is limited. We enrolled 43 Vietnamese CRC patients and 25 participants with non-cancerous colorectal polyps to study the colonic microbiomes and genomic diversity of Fusobacterium in this population, using a combination of 16S rRNA gene profiling, anaerobic microbiology, and whole genome analysis. Oral bacteria, including F. nucleatum and Leptotrichia, were significantly more abundant in the tumour microbiomes. We obtained 53 Fusobacterium genomes, representing 26 strains, from the saliva, tumour and non-tumour tissues of six CRC patients. Isolates from the gut belonged to diverse F. nucleatum subspecies (nucleatum, animalis, vincentii, polymorphum) and a potential new subspecies of Fusobacterium periodonticum. The Fusobacterium population within each individual was distinct and in some cases diverse, with minimal intra-clonal variation. Phylogenetic analyses showed that within four individuals, tumour-associated Fusobacterium were clonal to those isolated from non-tumour tissues. Genes encoding major virulence factors (Fap2 and RadD) showed evidence of horizontal gene transfer. Our work provides a framework to understand the genomic diversity of Fusobacterium within the CRC patients, which can be exploited for the development of CRC diagnostic and therapeutic options targeting this oncobacterium.

Toward FAIR Representations of Microbial Interactions

Despite an ever-growing number of data sets that catalog and characterize interactions between microbes in different environments and conditions, many of these data are neither easily accessible nor intercompatible. These limitations present a major challenge to microbiome research by hindering the streamlined drawing of inferences across studies. Here, we propose guiding principles to make microbial interaction data more findable, accessible, interoperable, and reusable (FAIR). We outline specific use cases for interaction data that span the diverse space of microbiome research, and discuss the untapped potential for new insights that can be fulfilled through broader integration of microbial interaction data. These include, among others, the design of intercompatible synthetic communities for environmental, industrial, or medical applications, and the inference of novel interactions from disparate studies. Lastly, we envision potential trajectories for the deployment of FAIR microbial interaction data based on existing resources, reporting standards, and current momentum within the community.

Reconstruction of a Soil Microbial Network Induced by Stress Temperature

The microbial community is viewed as a network of diverse microorganisms connected by various interspecific interactions. While the stress gradient hypothesis (SGH) predicts that positive interactions are favored in more stressful environments, the prediction has been less explored in complex microbial communities due to the challenges of identifying interactions. Here, by applying a nonlinear time series analysis to the amplicon-based diversity time series data of the soil microbiota cultured under less stressful (30°C) or more stressful (37°C) temperature conditions, we show how the microbial network responds to temperature stress. While the genera that persisted only under the less stressful condition showed fewer positive effects, the genera that appeared only under the more stressful condition received more positive effects, in agreement with SGH. However, temperature difference also induced reconstruction of the community network, leading to an increased proportion of negative interactions at the whole-community level. The anti-SGH pattern can be explained by the stronger competition caused by increased metabolic rate and population densities. IMPORTANCE By combining amplicon-based diversity survey with recently developed nonlinear analytical tools, we successfully determined the interaction networks of more than 150 natural soil microbial genera under less or more temperature stress and explored the applicability of the stress gradient hypothesis to soil microbiota, shedding new light on the well-known hypothesis.

Meta-analysis of microbiome association networks reveal patterns of dysbiosis in diseased microbiomes

The human gut microbiome is composed of a diverse and dynamic population of microbial species which play key roles in modulating host health and physiology. While individual microbial species have been found to be associated with certain disease states, increasing evidence suggests that higher-order microbial interactions may have an equal or greater contribution to host fitness. To better understand microbial community dynamics, we utilize networks to study interactions through a meta-analysis of microbial association networks between healthy and disease gut microbiomes. Taking advantage of the large number of metagenomes derived from healthy individuals and patients with various diseases, together with recent advances in network inference that can deal with sparse compositional data, we inferred microbial association networks based on co-occurrence of gut microbial species and made the networks publicly available as a resource (GitHub repository named GutNet). Through our meta-analysis of inferred networks, we were able to identify network-associated features that help stratify between healthy and disease states such as the differentiation of various bacterial phyla and enrichment of Proteobacteria interactions in diseased networks. Additionally, our findings show that the contributions of taxa in microbial associations are disproportionate to their abundances and that rarer taxa of microbial species play an integral part in shaping dynamics of microbial community interactions. Network-based meta-analysis revealed valuable insights into microbial community dynamics between healthy and disease phenotypes. We anticipate that the healthy and diseased microbiome association networks we inferred will become an important resource for human-related microbiome research.

Simultaneous methanogenesis and denitrification coupled with nitrifying biofilm for high-strength wastewater treatment: Performance and microbial mechanisms

A combined system consisting of an upflow blanket filter (UBF) and a moving-bed biofilm reactor (MBBR) was developed for the simultaneous removal of organic matters and ammonia from high-strength wastewater. With a constant COD of approximately 2000 mg/L and ammonium nitrogen in a series of concentrations (e.g., 50, 200 and 400 mg/L in stages I to III) of the influent wastewater, the removal efficiencies of COD, ammonium nitrogen and total nitrogen reached 96.10%-98.19%, 100%, and 79.12%-82.15%, respectively. With the increase of influent ammonia nitrogen concentration, the specific methanogenic activity of the UBF granules decreased significantly, while the specific denitrification rates of the UBF granules and specific nitrification rates of the MBBR biofilms increased significantly. Microbial community analysis showed that Methanobacterium and Methanosaeta were the dominant methanogens in the UBF granules, while Candidatus Competibacter, Thauera and Acinetobacter were identified as dominant denitrifiers. In addition, nitrifiers were enriched in MBBR biofilms at 11.33% and 13.87% of the average abundance of Nitrosomonas and Nitrospira, respectively, at stage III (influent ammonium at 400 mg/L, COD/NH₄⁺-N =  5). The ecological network analysis, including full-networks and sub-networks, indicated that the interactions between methanogens and denitrifiers in the UBF granules were strong when the influent ammonium concentration reached 400 mg/L. No intensive interactions were observed among the functional bacteria in the MBBR biofilms over the entire operation. Overall, this study provides a new strategy for the application and construction of efficient biological processes to achieve simultaneous removal of organic matter and nitrogen for high-strength wastewater treatment.

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.

Species abundance correlations carry limited information about microbial network interactions

Unraveling the network of interactions in ecological communities is a daunting task. Common methods to infer interspecific interactions from cross-sectional data are based on co-occurrence measures. For instance, interactions in the human microbiome are often inferred from correlations between the abundances of bacterial phylogenetic groups across subjects. We tested whether such correlation-based methods are indeed reliable for inferring interaction networks. For this purpose, we simulated bacterial communities by means of the generalized Lotka-Volterra model, with variation in model parameters representing variability among hosts. Our results show that correlations can be indicative for presence of bacterial interactions, but only when measurement noise is low relative to the variation in interaction strengths between hosts. Indication of interaction was affected by type of interaction network, process noise and sampling under non-equilibrium conditions. The sign of a correlation mostly coincided with the nature of the strongest pairwise interaction, but this is not necessarily the case. For instance, under rare conditions of identical interaction strength, we found that competitive and exploitative interactions can result in positive as well as negative correlations. Thus, cross-sectional abundance data carry limited information on specific interaction types. Correlations in abundance may hint at interactions but require independent validation.

The bacteria in and on our body (the human microbiome) largely determine how our body functions, and whether we stay healthy or get sick. These bacteria do not live on their own, but interact among each other and with their human host. Finding out which bacteria interact with each other is cumbersome, but patterns of joint occurrence between species might provide a clue to their ecological dependencies. We investigated whether correlations in species abundance can be used for the purpose of ecological network reconstruction. We simulated different bacterial communities with known interactions according to a theoretical population model. After having collected virtual samples from our simulated data, we performed a correlation analysis and then compared the correlation network with our known interaction network. We found that correlations can be informative for underlying interactions, but ecological conclusions should be drawn carefully. An obvious limitation of correlation analysis is that direction of interaction cannot be recovered from co-occurrence data, making correlations insensitive for detection of asymmetric interactions. In addition, we found that competitive and exploitative interactions can induce positive as well as negative correlations. We recommend careful interpretation and validation when inferring networks from cross-sectional abundance data.

Identification of Functional Microbial Modules Through Network-Based Analysis of Meta-Microbial Features Using Matrix Factorization

As the microbiome is composed of a variety of microbial interactions, it is imperative in microbiome research to identify a microbial sub-community that collectively conducts a specific function. However, current methodologies have been highly limited to analyzing conditional abundance changes of individual microorganisms without considering group-wise collective microbial features. To overcome this limitation, we developed a network-based method using nonnegative matrix factorization (NMF) to identify functional meta-microbial features (MMFs) that, as a group, better discriminate specific environmental conditions of samples using microbiome data. As proof of concept, large-scale human microbiome data collected from different body sites were used to identify body site-specific MMFs by applying NMF. The statistical test for MMFs led us to identify highly discriminative MMFs on sample classes, called synergistic MMFs (SYMMFs). Finally, we constructed a SYMMF-based microbial interaction network (SYMMF-net) by integrating all of the SYMMF information. Network analysis revealed core microbial modules closely related to critical sample properties. Similar results were also found when the method was applied to various disease-associated microbiome data. The developed method interprets high-dimensional microbiome data by identifying functional microbial modules on sample properties and intuitively representing their systematic relationships via a microbial network.

Bacterial species metabolic interaction network for deciphering the lignocellulolytic system in fungal cultivating termite gut microbiota

Fungus-cultivating termite Odontotermes badius developed a mutualistic association with Termitomyces fungi for the plant material decomposition and providing a food source for the host survival. The mutualistic relationship sifted the microbiome composition of the termite gut and Termitomyces fungal comb. Symbiotic bacterial communities in the O. badius gut and fungal comb have been studied extensively to identify abundant bacteria and their lignocellulose degradation capabilities. Despite several metagenomic studies, the species-wide metabolic interaction pattern of bacterial communities in termite gut and fungal comb remains unclear. The bacterial species metabolic interaction network (BSMIN) has been constructed with 230 bacteria identified from the O. badius gut and fungal comb microbiota. The network portrayed the metabolic map of the entire microbiota and highlighted several inter-species biochemical interactions like cross-feeding, metabolic interdependency, and competition. Further, the reconstruction and analysis of the bacterial influence network (BIN) quantified the positive and negative pairwise influences in the termite gut and fungal comb microbial communities. Several key macromolecule degraders and fermentative microbial entities have been identified by analyzing the BIN. The mechanistic interplay between these influential microbial groups and the crucial glycoside hydrolases (GH) enzymes produced by the macromolecule degraders execute the community-wide functionality of lignocellulose degradation and subsequent fermentation. The metabolic interaction pattern between the nine influential microbial species has been determined by considering them growing in a synthetic microbial community. Competition (30%), parasitism (47%), and mutualism (17%) were predicted to be the major mode of metabolic interaction in this synthetic microbial community. Further, the antagonistic metabolic effect was found to be very high in the metabolic-deprived condition, which may disrupt the community functionality. Thus, metabolic interactions of the crucial bacterial species and their GH enzyme cocktail identified from the O. badius gut and fungal comb microbiota may provide essential knowledge for developing a synthetic microcosm with efficient lignocellulolytic machinery.

Achieving simultaneous nitrification, denitrification, and phosphorus removal in pilot-scale flow-through biofilm reactor with low dissolved oxygen concentrations: Performance and mechanisms

In this pilot-scale study, a flow-through biofilm reactor (FTBR) was investigated for municipal wastewater treatment. The removal efficiencies for ammonium, total nitrogen, total phosphorus, and chemical oxygen demand were 87.2 ± 17.9%, 61.1 ± 13.9%, 83.5 ± 11.9%, and 92.6 ± 1.7%, respectively, at low dissolved oxygen concentrations (averaged at 0.59 mg/L), indicating the feasibility and robustness of the FTBR for a simultaneous nitrification, denitrification, and phosphorous removal (SNDPR) process. The co-occurrence network of bacteria in the dynamic biofilm was complex, with equivalent bacterial cooperation and competition. Nevertheless, the bacterial interactions in the suspended sludge were mainly cooperative. The presence of dynamic biofilms increased bacterial diversity by creating niche differentiation, which enriched keystone species closely related to nutrient removal. Overall, this study provides a novel FTBR-based SNDPR process and reveals the ecological mechanisms responsible for nutrient removal.

Co-occurrence network analysis reveals loss of microbial interactions in anaerobic digester subjected to repeated organic load shocks

Fluctuations in the anaerobic digestion (AD) organic loading rate (OLR) cause shocks to the AD microbiome, which lead to unstable methane productivity. Managing these fluctuations requires a larger digester, which is impractical for community-scale applications, limiting the potential of AD in advancing a circular economy. To allow operation of small-scale AD while managing OLR fluctuations, we need to tackle the issue through elucidation of the microbial community dynamics via 16S rRNA gene sequencing. This study elucidated the interrelation of the AD performance and the dynamics of the microbial interactions within its microbiome in response to repeated high OLR shocks at different frequencies. The OLR shocks were equivalent to 4 times the baseline OLR of 2 g VS/L/d. We found that less frequent organic load shocks result to deterioration of methane productivity. Co-occurrence network analysis shows that this coincides with the breakdown of the microbiome network structure. This suggests loss of microbial interactions necessary in maintaining stable AD. Identification of species influencing the network structure revealed that a species under the genus Anaerovorax has the greatest influence, while orders Spirochaetales and Synergistales represent the greatest number of the influential species. We inferred that the impact imposed by the OLR shocks shifted the microbiome activity towards biochemical pathways that are not contributing to methane production. Establishing a small-scale AD system that permits OLR fluctuations would require developing an AD microbiome resilient to infrequent organic loading shocks.

Cross-kingdom co-occurrence networks in the plant microbiome: Importance and ecological interpretations

Microbial co-occurrence network analysis is being widely used for data exploration in plant microbiome research. Still, challenges lie in how well these microbial networks represent natural microbial communities and how well we can interpret and extract eco-evolutionary insights from the networks. Although many technical solutions have been proposed, in this perspective, we touch on the grave problem of kingdom-level bias in network representation and interpretation. We underscore the eco-evolutionary significance of using cross-kingdom (bacterial-fungal) co-occurrence networks to increase the network’s representability of natural communities. To do so, we demonstrate how ecosystem-level interpretation of plant microbiome evolution changes with and without multi-kingdom analysis. Then, to overcome oversimplified interpretation of the networks stemming from the stereotypical dichotomy between bacteria and fungi, we recommend three avenues for ecological interpretation: (1) understanding dynamics and mechanisms of co-occurrence networks through generalized Lotka-Volterra and consumer-resource models, (2) finding alternative ecological explanations for individual negative and positive fungal-bacterial edges, and (3) connecting cross-kingdom networks to abiotic and biotic (host) environments.

Adversarial attacks on voter model dynamics in complex networks

This paper investigates adversarial attacks conducted to distort voter model dynamics in complex networks. Specifically, a simple adversarial attack method is proposed to hold the state of opinions of an individual closer to the target state in the voter model dynamics. This indicates that even when one opinion is the majority the vote outcome can be inverted (i.e., the outcome can lean toward the other opinion) by adding extremely small (hard-to-detect) perturbations strategically generated in social networks. Adversarial attacks are relatively more effective in complex (large and dense) networks. These results indicate that opinion dynamics can be unknowingly distorted.

MicNet toolbox: Visualizing and unraveling a microbial network

Applications of network theory to microbial ecology are an emerging and promising approach to understanding both global and local patterns in the structure and interplay of these microbial communities. In this paper, we present an open-source python toolbox which consists of two modules: on one hand, we introduce a visualization module that incorporates the use of UMAP, a dimensionality reduction technique that focuses on local patterns, and HDBSCAN, a clustering technique based on density; on the other hand, we have included a module that runs an enhanced version of the SparCC code, sustaining larger datasets than before, and we couple the resulting networks with network theory analyses to describe the resulting co-occurrence networks, including several novel analyses, such as structural balance metrics and a proposal to discover the underlying topology of a co-occurrence network. We validated the proposed toolbox on 1) a simple and well described biological network of kombucha, consisting of 48 ASVs, and 2) we validate the improvements of our new version of SparCC. Finally, we showcase the use of the MicNet toolbox on a large dataset from Archean Domes, consisting of more than 2,000 ASVs. Our toolbox is freely available as a github repository (https://github.com/Labevo/MicNetToolbox), and it is accompanied by a web dashboard (http://micnetapplb-1212130533.us-east-1.elb.amazonaws.com) that can be used in a simple and straightforward manner with relative abundance data. This easy-to-use implementation is aimed to microbial ecologists with little to no experience in programming, while the most experienced bioinformatics will also be able to manipulate the source code's functions with ease.

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.

Diurnal and eating-associated microbial patterns revealed via high-frequency saliva sampling

The oral microbiome is linked to oral and systemic health, but its fluctuation under frequent daily activities remains elusive. Here, we sampled saliva at 10- to 60-min intervals to track the high-resolution microbiome dynamics during the course of human activities. This dense time series data showed that eating activity markedly perturbed the salivary microbiota, with tongue-specific Campylobacter concisus and Oribacterium sinus and dental plaque-specific Lautropia mirabilis, Rothia aeria, and Neisseria oralis increased after every meal in a temporal order. The observation was reproducible in multiple subjects and across an 11-mo period. The microbiome composition showed significant diurnal oscillation patterns at different taxonomy levels with Prevotella/Alloprevotella increased at night and Bergeyella HMT 206/Haemophilus slowly increased during the daytime. We also identified microbial co-occurring patterns in saliva that are associated with the intricate biogeography of the oral microbiome. Microbial source tracking analysis showed that the contributions of distinct oral niches to the salivary microbiome were dynamically affected by daily activities, reflecting the role of saliva in exchanging microbes with other oral sites. Collectively, our study provides insights into the temporal microbiome variation in saliva and highlights the need to consider daily activities and diurnal factors in design of oral microbiome studies.

Gut Microbial Shifts Indicate Melanoma Presence and Bacterial Interactions in a Murine Model

Through a multitude of studies, the gut microbiota has been recognized as a significant influencer of both homeostasis and pathophysiology. Certain microbial taxa can even affect treatments such as cancer immunotherapies, including the immune checkpoint blockade. These taxa can impact such processes both individually as well as collectively through mechanisms from quorum sensing to metabolite production. Due to this overarching presence of the gut microbiota in many physiological processes distal to the GI tract, we hypothesized that mice bearing tumors at extraintestinal sites would display a distinct intestinal microbial signature from non-tumor-bearing mice, and that such a signature would involve taxa that collectively shift with tumor presence. Microbial OTUs were determined from 16S rRNA genes isolated from the fecal samples of C57BL/6 mice challenged with either B16-F10 melanoma cells or PBS control and analyzed using QIIME. Relative proportions of bacteria were determined for each mouse and, using machine-learning approaches, significantly altered taxa and co-occurrence patterns between tumor- and non-tumor-bearing mice were found. Mice with a tumor had elevated proportions of Ruminococcaceae, Peptococcaceae.g_rc4.4, and Christensenellaceae, as well as significant information gains and ReliefF weights for Bacteroidales.f__S24.7, Ruminococcaceae, Clostridiales, and Erysipelotrichaceae. Bacteroidales.f__S24.7, Ruminococcaceae, and Clostridiales were also implicated through shifting co-occurrences and PCA values. Using these seven taxa as a melanoma signature, a neural network reached an 80% tumor detection accuracy in a 10-fold stratified random sampling validation. These results indicated gut microbial proportions as a biosensor for tumor detection, and that shifting co-occurrences could be used to reveal relevant taxa.

Statistical Methods for Microbiome Compositional Data Network Inference: A Survey

Microbes can be found almost everywhere in the world. They are not isolated, but rather interact with each other and establish connections with their living environments. Studying these interactions is essential to an understanding of the organization and complex interplay of microbial communities, as well as the structure and dynamics of various ecosystems. A widely used approach toward this objective involves the inference of microbiome interaction networks. However, owing to the compositional, high-dimensional, sparse, and heterogeneous nature of observed microbial data, applying network inference methods to estimate their associations is challenging. In addition, external environmental interference and biological concerns also make it more difficult to deal with the network inference. In this article, we provide a comprehensive review of emerging microbiome interaction network inference methods. According to various research targets, estimated networks are divided into four main categories: correlation networks, conditional correlation networks, mixture networks, and differential networks. Their assumptions, high-level ideas, advantages, as well as limitations, are presented in this review. Since real microbial interactions can be complex and dynamic, no unifying method has, to date, captured all the aspects of interest. In addition, we discuss the challenges now confronting current microbial interaction study and future prospects. Finally, we point out several feasible directions of microbial network inference analysis and highlight that future research requires the joint promotion of statistical computation methods and experimental techniques.

Where environmental microbiome meets its host: subway and passenger microbiome relationships

Subways are urban transport systems with high capacity. Every day around the world, there are more than 150 million subway passengers. Since 2013, thousands of microbiome samples from various subways worldwide have been sequenced. Skin bacteria and environmental organisms dominate the subway microbiomes. The literature has revealed common bacterial groups in subway systems; even so, it is possible to identify cities by their microbiome. Low-frequency bacteria are responsible for specific bacterial fingerprints of each subway system. Furthermore, daily subway commuters leave their microbial clouds and interact with other passengers. Microbial exchange is quite fast; the hand microbiome changes within minutes, and after cleaning the handrails, the bacteria are re-established within minutes. To investigate new taxa and metabolic pathways of subway microbial communities, several high-quality metagenomic-assembled genomes (MAG) have been described. Subways are harsh environments unfavorable for microorganism growth. However, recent studies have observed a wide diversity of viable and metabolically active bacteria. Understanding which bacteria are living, dormant, or dead allows us to propose realistic ecological interactions. Questions regarding the relationship between humans and the subway microbiome, particularly the microbiome effects on personal and public health, remain unanswered. This review summarizes our knowledge of subway microbiomes and their relationship with passenger microbiomes.