FastSpar: rapid and scalable correlation estimation for compositional data

Amendment-driven soil health restoration through soil pH and microbial robustness in a Cd/Cu-combined acidic soil: A ten-year in-situ field experiment

Soil health arguably depends on biodiversity and has received wide attention in heavy-metal (HM) contaminated farmland remediation in recent years. However, long-term effects and mechanisms of soil amendment remain poorly understood with respect to soil microbal community. In this in-situ field study, four soil amendments (attapulgite-At, apatite-Ap, montmorillonite-M, lime-L) at three rates were applied once only for ten years in a cadmium (Cd)-copper (Cu) contaminated paddy soil deprecated for over five years. Results showed that after ten years and in compared with CK (no amendment), total Cd concentration and its risk in plot soils were not altered by amendments (p > 0.05), but total Cu concentration and its risk were significantly increased by both Ap and L, especially the former, rather than At and M (p < 0.05), through increased soil pH and enhanced bacterial alpha diversity as well as plant community. Soil microbial communities were more affected by amendment type (30%) than dosage (11%), microbial network characteristics were dominated by rare taxa, and soil multifunctionality was improved in Ap- and L-amended soils. A structural equation model (SEM) indicated that 57.3% of soil multifunctionality variances were accounted for by soil pH (+0.696) and microbial network robustness (-0.301). Moreover, microbial robustness could be potentially used as an indicator of soil multifunctionality, and Ap could be optimized to improve soil health in combined with biomass removal. These findings would advance the understanding of soil microbial roles, especially its network robustness, on soil multifunctionality for the remediation of metal contaminated soils and metal control management strategies in acidic soils. ENVIRONMENTAL IMPLICATION: Farmland soil contamination by heavy metals (HMs) has been becoming a serious global environmental challenge. However, most studies have been conducted over the short term, leading to a gap in the long-term remediation efficiency and ecological benefits of soil amendments. For the successful deployment of immobilization technologies, it is critical to understand the long-term stability of the immobilized HMs and soil health. Our study, to the best of our knowlege, is the first to state the long-term effects and mechanisms of soil amendments on soil health and optimize an effective and eco-friendly amendment for long-term Cd/Cu immobilization.

Good and bad dispositions between archaea and bacteria in the human gut: New insights from metagenomic survey and co-occurrence analysis

Archaea are an understudied component of the human microbiome. In this study, the gut archaeome and bacteriome of 60 healthy adults from different region were analyzed by whole-genome shotgun sequencing. Archaea were ubiquitously found in a wide range of abundances, reaching up to 7.2 %. The dominant archaeal phylum was Methanobacteriota, specifically the family Methanobacteriaceae, encompassing more than 50 % of Archaea in 50 samples. The previously underestimated Thermoplasmatota, mostly composed of Methanomassiliicoccaceae, dominated in 10 subjects (>50 %) and was present in all others except one. Halobacteriota, the sole other archaeal phylum, occurred in negligible concentration, except for two samples (4.6-4.8 %). This finding confirmed that the human gut archaeome is primarily composed of methanogenic organisms and among the known methanogenic pathway: i) hydrogenotrophic reduction of CO2 is the predominant, being the genus Methanobrevibacter and the species Methanobrevibacter smithii the most abundant in the majority of the samples; ii) the second pathway, that involved Methanomassiliicoccales, was the hydrogenotrophic reduction of methyl-compounds; iii) dismutation of acetate or methyl-compounds seemed to be absent. Co-occurrence analysis allowed to unravel correlations between Archaea and Bacteria that shapes the overall structure of the microbial community, allowing to depict a clearer picture of the human gut archaeome.

Diverse plasmid systems and their ecology across human gut metagenomes revealed by PlasX and MobMess

Plasmids alter microbial evolution and lifestyles by mobilizing genes that often confer fitness in changing environments across clades. Yet our ecological and evolutionary understanding of naturally occurring plasmids is far from complete. Here we developed a machine-learning model, PlasX, which identified 68,350 non-redundant plasmids across human gut metagenomes and organized them into 1,169 evolutionarily cohesive 'plasmid systems' using our sequence containment-aware network-partitioning algorithm, MobMess. Individual plasmids were often country specific, yet most plasmid systems spanned across geographically distinct human populations. Cargo genes in plasmid systems included well-known determinants of fitness, such as antibiotic resistance, but also many others including enzymes involved in the biosynthesis of essential nutrients and modification of transfer RNAs, revealing a wide repertoire of likely fitness determinants in complex environments. Our study introduces computational tools to recognize and organize plasmids, and uncovers the ecological and evolutionary patterns of diverse plasmids in naturally occurring habitats through plasmid systems.

The gut ileal mucosal virome is disturbed in patients with Crohn's disease and exacerbates intestinal inflammation in mice

Gut bacteriome dysbiosis is known to be implicated in the pathogenesis of inflammatory bowel disease (IBD). Crohn's disease (CD) is an IBD subtype with extensive mucosal inflammation, yet the mucosal virome, an empirical modulator of the bacteriome and mucosal immunity, remains largely unclear regarding its composition and role. Here, we exploited trans-cohort CD patients and healthy individuals to compositionally and functionally investigate the small bowel (terminal ileum) virome and bacteriome. The CD ileal virome was characterised by an under-representation of both lytic and temperate bacteriophages (especially those targeting bacterial pathogens), particularly in patients with flare-up. Meanwhile, the virome-bacteriome ecology in CD ileal mucosa was featured by a lack of Bifidobacterium- and Lachnospiraceae-led mutualistic interactions between bacteria and bacteriophages; surprisingly it was more pronounced in CD remission than flare-up, underlining the refractory and recurrent nature of mucosal inflammation in CD. Lastly, we substantiated that ileal virions from CD patients causally exacerbated intestinal inflammation in IBD mouse models, by reshaping a gut virome-bacteriome ecology preceding intestinal inflammation (microbial trigger) and augmenting microbial sensing/defence pathways in the intestine cells (host response). Altogether, our results highlight the significance of mucosal virome in CD pathogenesis and importance of mucosal virome restoration in CD therapeutics.

Loss of symbiotic and increase of virulent bacteria through microbial networks in Lynch syndrome colon carcinogenesis

Purpose

Through a pilot study, we performed whole gut metagenomic analysis in 17 Lynch syndrome (LS) families, including colorectal cancer (CRC) patients and their healthy first-degree relatives. In a second asymptomatic LS cohort (n=150) undergoing colonoscopy-screening program, individuals with early precancerous lesions were compared to those with a normal colonoscopy. Since bacteria are organized into different networks within the microbiota, we compared related network structures in patients and controls.

Experimental design

Fecal prokaryote DNA was extracted prior to colonoscopy for whole metagenome (n=34, pilot study) or 16s rRNA sequencing (validation study). We characterized bacteria taxonomy using Diamond/MEGAN6 and DADA2 pipelines and performed differential abundances using Shaman website. We constructed networks using SparCC inference tools and validated the construction's accuracy by performing qPCR on selected bacteria.

Results

Significant differences in bacterial communities in LS-CRC patients were identified, with an enrichment of virulent bacteria and a depletion of symbionts compared to their first-degree relatives. Bacteria taxa in LS asymptomatic individuals with colonic precancerous lesions (n=79) were significantly different compared to healthy individuals (n=71). The main bacterial network structures, constructed based on bacteria-bacteria correlations in CRC (pilot study) and in asymptomatic precancerous patients (validation-study), showed a different pattern than in controls. It was characterized by virulent/symbiotic co-exclusion in both studies and illustrated (validation study) by a higher Escherichia/Bifidobacterium ratio, as assessed by qPCR.

Conclusion

Enhanced fecal virulent/symbiotic bacteria ratios influence bacterial network structures. As an early event in colon carcinogenesis, these ratios can be used to identify asymptomatic LS individual with a higher risk of CRC.

Host-microbiota interactions contributing to the heterogeneous tumor microenvironment in colorectal cancer

Colorectal cancer (CRC) exhibits pronounced heterogeneity and is categorized into four widely accepted consensus molecular subtypes (CMSs) with unique tumor microenvironments (TMEs). However, the intricate landscape of the microbiota and host-microbiota interactions within these TMEs remains elusive. Using RNA-sequencing data from The Cancer Genome Atlas, we analyzed the host transcriptomes and intratumoral microbiome profiles of CRC samples. Distinct host genes and microbial genera were identified among the CMSs. Immune microenvironments were evaluated using CIBERSORTx and ESTIMATE, and microbial coabundance patterns were assessed with FastSpar. Through LASSO penalized regression, we explored host-microbiota associations for each CMS. Our analysis revealed distinct host gene signatures within the CMSs, which encompassed ferroptosis-related genes and specific immune microenvironments. Moreover, we identified 293, 153, 66, and 109 intratumoral microbial genera with differential abundance, and host-microbiota associations contributed to distinct TMEs, characterized by 829, 1,270, 634, and 1,882 robust gene-microbe associations for each CMS in CMS1-CMS4, respectively. CMS1 featured inflammation-related HSF1 activation and gene interactions within the endothelin pathway and Flammeovirga. Integrin-related genes displayed positive correlations with Sutterella in CMS2, whereas CMS3 spotlighted microbial associations with biosynthetic and metabolic pathways. In CMS4, genes involved in collagen biosynthesis showed positive associations with Sutterella, contributing to disruptions in homeostasis. Notably, immune-rich subtypes exhibited pronounced ferroptosis dysregulation, potentially linked to tissue microbial colonization. This comprehensive investigation delineates the diverse landscapes of the TME within each CMS, incorporating host genes, intratumoral microbiota, and their complex interactions. These findings shed light on previously uncharted mechanisms underpinning CRC heterogeneity and suggest potential therapeutic targets.NEW & NOTEWORTHY This study determined the following: 1) providing a comprehensive landscape of consensus molecular subtype (CMS)-specific tumor microenvironments (TMEs); 2) constructing CMS-specific networks, including host genes, intratumoral microbiota, and enriched pathways, analyzing their associations to uncover unique patterns that demonstrate the intricate interplay within the TME; and 3) revealing a connection between immune-rich subtypes and ferroptosis activation, suggesting a potential regulatory role of the microbiota in ferroptosis dysregulation of the colorectal cancer TME.

Effect of aridity on the β-diversity of alpine soil potential diazotrophs: insights into community assembly and co-occurrence patterns

Microbial diversity plays a vital role in the maintenance of ecosystem functions. However, the current understanding of mechanisms that shape microbial diversity along environmental gradients at broad spatial scales is relatively limited, especially for specific functional groups, such as potential diazotrophs. Here, we conducted an aridity-gradient transect survey from 60 sites across the Tibetan Plateau, the largest alpine ecosystem of the planet, to investigate the ecological processes (e.g., local species pools, community assembly processes, and co-occurrence patterns) that underlie the β-diversity of alpine soil potential diazotrophic communities. We found that aridity strongly and negatively affected the abundance, richness, and β-diversity of soil diazotrophs. Diazotrophs displayed a distance-decay pattern along the aridity gradient, with organisms living in lower aridity habitats having a stronger distance-decay pattern. Arid habitats had lower co-occurrence complexity, including the number of edges and vertices, the average degree, and the number of keystone taxa, as compared with humid habitats. Local species pools explained limited variations in potential diazotrophic β-diversity. In contrast, co-occurrence patterns and stochastic processes (e.g., dispersal limitation and ecological drift) played a significant role in regulating potential diazotrophic β-diversity. The relative importance of stochastic processes and co-occurrence patterns changed with increasing aridity, with stochastic processes weakening whereas that of co-occurrence patterns enhancing. The genera Geobacter and Paenibacillus were identified as keystone taxa of co-occurrence patterns that are associated with β-diversity. In summary, aridity affects the co-occurrence patterns and community assembly by regulating soil and vegetation characteristics and ultimately shapes the β-diversity of potential diazotrophs. These findings highlight the importance of co-occurrence patterns in structuring microbial diversity and advance the current understanding of mechanisms that drive belowground communities.IMPORTANCERecent studies have shown that community assembly processes and species pools are the main drivers of β-diversity in grassland microbial communities. However, co-occurrence patterns can also drive β-diversity formation by influencing the dispersal and migration of species, the importance of which has not been reported in previous studies. Assessing the impact of co-occurrence patterns on β-diversity is important for understanding the mechanisms of diversity formation. Our study highlights the influence of microbial co-occurrence patterns on β-diversity and combines the drivers of community β-diversity with drought variation, revealing that drought indirectly affects β-diversity by influencing diazotrophic co-occurrence patterns and community assembly.

Marked variations in diversity and functions of gut microbiota between wild and domestic stag beetle Dorcus Hopei Hopei

BACKGROUND

Although stag beetles are a popular saprophytic insect, their gut microbiome has been poorly studied. Here, 16 S rRNA gene sequencing was employed to reveal the gut microbiota composition and functional variations between wild and domestic Dorcus hopei hopei (Dhh) larval individuals.

RESULTS

The results indicated a significant difference between the wild and domestic Dhh gut microbiota., the domestic Dhh individuals contained more gut microbial taxa (e.g. genera Ralstonia and Methyloversatilis) with xenobiotic degrading functions. The wild Dhh possesses gut microbiota compositions (e.g. Turicibacter and Tyzzerella ) more appropriate for energy metabolism and potential growth. This study furthermore assigned all Dhh individuals by size into groups for data analysis; which indicated limited disparities between the gut microbiota of different-sized D. hopei hopei larvae.

CONCLUSION

The outcome of this study illustrated that there exists a significant discrepancy in gut microbiota composition between wild and domestic Dhh larvae. In addition, the assemblage of gut microbiome in Dhh was primarily attributed to environmental influences instead of individual differences such as developmental potential or size. These findings will provide a valuable theoretical foundation for the protection of wild saprophytic insects and the potential utilization of the insect-associated intestinal microbiome in the future.

Aridity drives the variability of desert soil microbiomes across north-western China

Dryland covers >35 % of the terrestrial surface and the global extent of dryland increases due to the forecasted increase in aridity driven by climate change. Due to the climate change-driven aridity ecosystems, deserts provide one of the most hostile environments for microbial life and survival. Therefore, a detailed study was carried out to explore the deserts with different aridity levels (exposed to severe climate change) influence on microbial (bacteria, fungi, and protist) diversity patterns, assembly processes, and co-occurrence. The results revealed that the aridity (semi-arid, arid, and hyper-arid) patterns caused distinct changes in environmental heterogeneity in desert ecosystems. Similarly, microbial diversities were also reduced with increasing the aridity pattern, and it was found that environmental heterogeneity is highly involved in affecting microbial diversities under different ecological niches. Interestingly, it was found that certain microbes, including bacterial (Firmicutes), fungal (Sordariomycetes), and protistan (Ciliophora) abundance increased with increasing aridity levels, indicating that these microbes might possess the capability to tolerate the environmental stress conditions. Moreover, microbial community turnover analysis revealed that bacterial diversities followed homogenous selection, whereas fungi and protists were mostly driven by the dispersal limitation pattern. Co-occurrence network analysis showed that hyper-arid and arid conditions tightened the bacterial and fungal communities and had more positive associations compared to protistan. In conclusion, multiple lines of evidence were provided to shed light on the habitat specialization impact on microbial (bacteria, fungi, and protists) communities and composition under different desert ecosystems.

Microbial communities living inside plant leaves or on the leaf surface are differently shaped by environmental cues

Leaf-associated microbial communities can promote plant health and resistance to biotic and abiotic stresses. However, the importance of environmental cues in the assembly of the leaf endo- and epi-microbiota remains elusive. Here, we aimed to investigate the impact of seasonal environmental variations, on the establishment of the leaf microbiome, focusing on long-term changes (five years) in bacterial, fungal, and nonfungal eukaryotic communities colonizing the surface and endosphere of six wild Arabidopsis thaliana populations. While leaf-microbial communities were found to be highly stochastic, the leaf niche had a predominant importance with endophytic microbial communities consistently exhibiting a lower diversity and variability. Among environmental factors, radiation- and humidity-related factors are the most important drivers of diversity patterns in the leaf, with stronger effects on epiphytic communities. Using linear models, we identified 30 important genera whose relative abundance in leaf compartments could be modeled from environmental variables, suggesting specific niche preferences for these taxa. With the hypothesis that environmental factors could impact interactions within microbial communities, we analyzed the seasonal patterns of microbial interaction networks across leaf compartments. We showed that epiphytic networks are more complex than endophytic and that the complexity and connectivity of these networks are partially correlated with the mentioned environmental cues. Our results indicate that humidity and solar radiation function as major environmental cues shaping the phyllosphere microbiome at both micro (leaf compartment) and macro (site) scales. These findings could have practical implications for predicting and developing field-adapted microbes in the face of global change.

Characterizing Gut Microbiota in Older Chinese Adults with Cognitive Impairment: A Cross-Sectional Study

Background

Cognitive impairment is a clinical manifestation that occurs in the course of dementia like Alzheimer's disease. The association between cognitive impairment and gut microbiota is unclear.

Objective

We aimed to identify gut microbiota characteristics and key gut microbiota biomarkers associated with cognitive impairment in a relatively large cohort of older adults in China.

Methods

A total of 229 adults aged ≥60 years from Shenzhen, China were recruited into this cross-sectional study. Participants were divided into cognitive impairment (CI) and no cognitive impairment (NCI) groups according to the results of the Mini-Mental State Examination. Diversity analysis and network analysis were used to characterize the gut microbiota between the two groups. The linear discriminant analysis effect size method and machine learning approaches were sequentially performed to identify gut microbiota biomarkers. The relationship between biomarkers and lifestyle factors was explored using Transformation-based redundancy analysis (tb-RDA).

Results

A total of 74 CI participants and 131 NCI participants were included in the analysis. The CI group demonstrated lower α-diversity compared to the NCI group (Shannon: 2.798 versus 3.152, p < 0.001). The density of the gut microbiota interaction network was lower in the CI group (0.074) compared to the NCI group (0.081). Megamonas, Blautia, Pseudomonas, Stenotrophomonas, and Veillonella were key biomarkers for CI. The tb-RDA revealed that increased fruit intake and exercise contribute to a higher abundance of Megamonas, Blautia, and Veillonella.

Conclusions

We identified a significantly reduced abundance of certain beneficial gut microbiota in older Chinese adults with cognitive impairment.

Characterizing alterations in the gut microbiota following postpartum weight change

IMPORTANCE

Previous research has reported differences in the gut microbiome associated with varying body compositions. More specifically, within populations of mothers, the focus has been on the impact of gestational weight gain. This is the first study to examine postpartum weight change and its association with changes in the gut microbiome, similarly, it is the first to use a Latina cohort to do so. The results support the idea that weight gain may be an important factor in reducing gut microbiome network connectivity, diversity, and changing abundances of specific microbial taxa, all measures thought to impact host health. These results suggest that weight gain dynamically alters mothers' gut microbial communities in the first 6 months postpartum, with comparatively little change in mothers who lost weight; further research is needed to examine the health consequences of such changes.

Gut microbiota affects the estrus return of sows by regulating the metabolism of sex steroid hormones

BACKGROUND

Sex hormones play important roles in the estrus return of post-weaning sows. Previous studies have demonstrated a complex and bi-directional regulation between sex hormones and gut microbiota. However, the extent to which the gut microbiota affects estrus return of post-weaning sows is largely unknown.

RESULTS

In this study, we first screened 207 fecal samples from well-phenotyped sows by 16S rRNA gene sequencing and identified significant associations between microbes and estrus return of post-weaning sows. Using metagenomic sequencing data from 85 fecal samples, we identified 37 bacterial species that were significantly associated with estrus return. Normally returning sows were characterized by increased abundances of L. reuteri and P. copri and decreased abundances of B. fragilis, S. suis, and B. pseudolongum. The changes in gut microbial composition significantly altered the functional capacity of steroid hormone biosynthesis in the gut microbiome. The results were confirmed in a validation cohort. Significant changes in sex steroid hormones and related compounds were found between normal and non-return sows via metabolome analysis. An integrated analysis of differential bacterial species, metagenome, and fecal metabolome provided evidence that normal return-associated bacterial species L. reuteri and Prevotella spp. participated in the degradation of pregnenolone, progesterone, and testosterone, thereby promoting estrogen biosynthesis. Furthermore, the microbial metabolites related to sow energy and nutrient supply or metabolic disorders also showed relationships with sow estrus return.

CONCLUSIONS

An integrated analysis of differentially abundant bacterial species, metagenome, and fecal metabolome revealed the involvement of L. reuteri and Prevotella spp. in sow estrus return. These findings provide deep insight into the role of gut microbiota in the estrus return of post-weaning sows and the complex cross-talk between gut microbiota and sex hormones, suggesting that the manipulation of the gut microbiota could be an effective strategy to improve sow estrus return after weaning.

Vertical connectivity of microbiome and metabolome reveals depth-dependent variations across a deep cold-seep water column

Deciphering the vertical connectivity of oceanic microbiome and metabolome is crucial for understanding the carbon sequestration and achieving the carbon neutrality. However, we lack a systematic view of the interplay among particle transport, microbial community, and metabolic trait across depths. Through integrating the biogeochemical, microbial, and metabolic characteristics of a deep cold-seep water column (∼1989 m), we find the altered connectivity of microbial community and dissolved organic matter (DOM) across depths. Both the microbial communities (bacteria and protists) and DOM show a clear compositional connectivity from surface to the depth of 1000 m, highlighting the controls of sinking particle over microbial connectivity from the epipelagic to mesopelagic zone. However, due to the biological migration and ocean mixing, the fecal-associated bacteria and protistan consumers unexpectedly emerge and the degradation index of DOM substantially alters around 1000-1200 m. Collectively, we unveil the significance of multi-faceted particle dispersion, which supports the connectivity and variability of deep ocean microbial communities.

Phylogenetic and ecophysiological novelty of subsurface mercury methylators in mangrove sediments

Mangrove sediment is a crucial component in the global mercury (Hg) cycling and acts as a hotspot for methylmercury (MeHg) production. Early evidence has documented the ubiquity of well-studied Hg methylators in mangrove superficial sediments; however, their diversity and metabolic adaptation in the more anoxic and highly reduced subsurface sediments are lacking. Through MeHg biogeochemical assay and metagenomic sequencing, we found that mangrove subsurface sediments (20-100 cm) showed a less hgcA gene abundance but higher diversity of Hg methylators than superficial sediments (0-20 cm). Regional-scale investigation of mangrove subsurface sediments spanning over 1500 km demonstrated a prevalence and family-level novelty of Hg-methylating microbial lineages (i.e., those affiliated to Anaerolineae, Phycisphaerae, and Desulfobacterales). We proposed the candidate phylum Zixibacteria lineage with sulfate-reducing capacity as a currently understudied Hg methylator across anoxic environments. Unlike other Hg methylators, the Zixibacteria lineage does not use the Wood-Ljungdahl pathway but has unique capabilities of performing methionine synthesis to donate methyl groups. The absence of cobalamin biosynthesis pathway suggests that this Hg-methylating lineage may depend on its syntrophic partners (i.e., Syntrophobacterales members) for energy in subsurface sediments. Our results expand the diversity of subsurface Hg methylators and uncover their unique ecophysiological adaptations in mangrove sediments.

Multi-omics reveal mechanisms underlying chronic kidney disease of unknown etiology (CKDu) pathogenesis using zebrafish

Chronic kidney disease of unknown etiology (CKDu) is an endemic disease in the dry zone of farming communities, Sri Lanka. The drinking water in a CKDu prevalent area contains a high concentration of F-, hardness and other environmental pollutants, including heavy metals and microcystin, which are considered possible etiology of CKDu in these areas. Here, multi-omics data with host transcriptome, metabolome and gut microbiomes were obtained using simulated local drinking water of Sri Lanka after their exposure to adult zebrafish. Based on an integrated multi-omics analysis in the context of host physiology in the kidney injury samples with different pathologic grades, two common pathways necroptosis and purine metabolism were identified as potentially important pathways that affect kidney injury. The key metabolite acetyl adenylate in the purine metabolism pathway was significantly positively correlated with Comamonas (rho = 0.72) and significantly negatively correlated with Plesiomonas (rho = -0.58). This crucial metabolite and two key gut bacteria genera may not only be potential markers but also potential therapeutic targets in the uric acid metabolic pathway, which is an important factor in the pathogenesis of acute kidney injury (AKI) in general, as well as of chronic kidney disease (CKD). Based on this, we revealed the urea metabolism pathway of kidney injury in zebrafish and provided a new avenue for the treatment of CKDu in Sri Lanka.

Prescribed fire selects for a pyrophilous soil sub-community in a northern California mixed conifer forest

Prescribed fire is a critical strategy for mitigating the effects of catastrophic wildfires. While the above-ground response to fire has been well-documented, fewer studies have addressed the effect of prescribed fire on soil microorganisms. To understand how soil microbial communities respond to prescribed fire, we sampled four plots at a high temporal resolution (two burned, two controls), for 17 months, in a mixed conifer forest in northern California, USA. Using amplicon sequencing, we found that prescribed fire significantly altered both fungal and bacterial community structure. We found that most differentially abundant fungal taxa had a positive fold-change, while differentially abundant bacterial taxa generally had a negative fold-change. We tested the null hypothesis that these communities assembled due to neutral processes (i.e., drift and/or dispersal), finding that >90% of taxa fit this neutral prediction. However, a dynamic sub-community composed of burn-associated indicator taxa that were positively differentially abundant was enriched for non-neutral amplicon sequence variants, suggesting assembly via deterministic processes. In synthesizing these results, we identified 15 pyrophilous taxa with a significant and positive response to prescribed burns. Together, these results lay the foundation for building a process-driven understanding of microbial community assembly in the context of the classical disturbance regime of fire.

A citizen-science-enabled catalogue of the vaginal microbiome and associated factors

Understanding the composition and function of the vaginal microbiome is crucial for reproductive and overall health. Here we established the Isala citizen-science project to analyse the vaginal microbiomes of 3,345 women in Belgium (18-98 years) through self-sampling, 16S amplicon sequencing and extensive questionnaires. The overall vaginal microbiome composition was strongly tied to age, childbirth and menstrual cycle phase. Lactobacillus species dominated 78% of the vaginal samples. Specific bacterial taxa also showed to co-occur in modules based on network correlation analysis. Notably, the module containing Lactobacillus crispatus, Lactobacillus jensenii and Limosilactobacillus taxa was positively linked to oestrogen levels and contraceptive use and negatively linked to childbirth and breastfeeding. Other modules, named after abundant taxa (Gardnerella, Prevotella and Bacteroides), correlated with multiple partners, menopause, menstrual hygiene and contraceptive use. With this resource-rich vaginal microbiome map and associated health, life-course, lifestyle and dietary factors, we provide unique data and insights for follow-up clinical and mechanistic research.

Quantitative Amplicon Sequencing Is Necessary to Identify Differential Taxa and Correlated Taxa Where Population Sizes Differ

High-throughput, multiplexed-amplicon sequencing has become a core tool for understanding environmental microbiomes. As researchers have widely adopted sequencing, many open-source analysis pipelines have been developed to compare microbiomes using compositional analysis frameworks. However, there is increasing evidence that compositional analyses do not provide the information necessary to accurately interpret many community assembly processes. This is especially true when there are large gradients that drive distinct community assembly processes. Recently, sequencing has been combined with Q-PCR (among other sources of total quantitation) to generate "Quantitative Sequencing" (QSeq) data. QSeq more accurately estimates the true abundance of taxa, is a more reliable basis for inferring correlation, and, ultimately, can be more reliably related to environmental data to infer community assembly processes. In this paper, we use a combination of published data sets, synthesis, and empirical modeling to offer guidance for which contexts QSeq is advantageous. As little as 5% variation in total abundance among experimental groups resulted in more accurate inference by QSeq than compositional methods. Compositional methods for differential abundance and correlation unreliably detected patterns in abundance and covariance when there was greater than 20% variation in total abundance among experimental groups. Whether QSeq performs better for beta diversity analysis depends on the question being asked, and the analytic strategy (e.g., what distance metric is being used); for many questions and methods, QSeq and compositional analysis are equivalent for beta diversity analysis. QSeq is especially useful for taxon-specific analysis; QSeq transformation and analysis should be the default for answering taxon-specific questions of amplicon sequence data. Publicly available bioinformatics pipelines should incorporate support for QSeq transformation and analysis.

Importance of mobile genetic elements for dissemination of antimicrobial resistance in metagenomic sewage samples across the world

We are facing an ever-growing threat from increasing antimicrobial resistance (AMR) in bacteria. To mitigate this, we need a better understanding of the global spread of antimicrobial resistance genes (ARGs). ARGs are often spread among bacteria by horizontal gene transfer facilitated by mobile genetic elements (MGE). Here we use a dataset consisting of 677 metagenomic sequenced sewage samples from 97 countries or regions to study how MGEs are geographically distributed and how they disseminate ARGs worldwide. The ARGs, MGEs, and bacterial abundance were calculated by reference-based read mapping. We found systematic differences in the abundance of MGEs and ARGs, where some elements were prevalent on all continents while others had higher abundance in separate geographic areas. Different MGEs tended to be localized to temperate or tropical climate zones, while different ARGs tended to separate according to continents. This suggests that the climate is an important factor influencing the local flora of MGEs. MGEs were also found to be more geographically confined than ARGs. We identified several integrated MGEs whose abundance correlated with the abundance of ARGs and bacterial genera, indicating the ability to mobilize and disseminate these genes. Some MGEs seemed to be more able to mobilize ARGs and spread to more bacterial species. The host ranges of MGEs seemed to differ between elements, where most were associated with bacteria of the same family. We believe that our method could be used to investigate the population dynamics of MGEs in complex bacterial populations.

The Bacterial Microbiome of the Coral Skeleton Algal Symbiont Ostreobium Shows Preferential Associations and Signatures of Phylosymbiosis

Ostreobium, the major algal symbiont of the coral skeleton, remains understudied despite extensive research on the coral holobiont. The enclosed nature of the coral skeleton might reduce the dispersal and exposure of residing bacteria to the outside environment, allowing stronger associations with the algae. Here, we describe the bacterial communities associated with cultured strains of 5 Ostreobium clades using 16S rRNA sequencing. We shed light on their likely physical associations by comparative analysis of three datasets generated to capture (1) all algae associated bacteria, (2) enriched tightly attached and potential intracellular bacteria, and (3) bacteria in spent media. Our data showed that while some bacteria may be loosely attached, some tend to be tightly attached or potentially intracellular. Although colonised with diverse bacteria, Ostreobium preferentially associated with 34 bacterial taxa revealing a core microbiome. These bacteria include known nitrogen cyclers, polysaccharide degraders, sulphate reducers, antimicrobial compound producers, methylotrophs, and vitamin B12 producers. By analysing co-occurrence networks of 16S rRNA datasets from Porites lutea and Paragoniastrea australensis skeleton samples, we show that the Ostreobium-bacterial associations present in the cultures are likely to also occur in their natural environment. Finally, our data show significant congruence between the Ostreobium phylogeny and the community composition of its tightly associated microbiome, largely due to the phylosymbiotic signal originating from the core bacterial taxa. This study offers insight into the Ostreobium microbiome and reveals preferential associations that warrant further testing from functional and evolutionary perspectives.

Genome-centric investigation of the potential succession pattern in gut microbiota and altered functions under high-protein diet

Excessive intake of protein has been considered as a factor leading to intestinal microecological disorder, but why and how intestinal microbes change under the high-protein diet (HPD) have yet to be fully elucidated. Here, we performed 16S rRNA gene amplicon sequencing and metagenomic sequencing on contents of cecum, colon and feces from two groups of mice with standard diet (SD) and HPD. And then the microbial alteration of composition and function were deeply analyzed by using several statistical models and bioinformatic methods. Among the three niches, the microbes in the colon are observed to show the most significant change with lower alpha-diversity and higher beta-diversity after HPD. In addition, this alteration of microbial structure may be related to the replacement process and co-occurring community. Most species are also enriched or impoverished in the colon during this process. After analyzing the functional genes related to protein and carbohydrate hydrolysis in different niches, we found that the carbon source provided by poor carbohydrates compared with the rich protein may be the potential factor driving the enrichment of mucin degraders and desulphaters in the colon under HPD. Therefore, our study provided a new insight to understand the underlying mechanism of HPD affecting intestinal health from the perspective of microbial functional ecology.

The origin of suspended particulate matter in the Great Barrier Reef

River run-off has long been regarded as the largest source of organic-rich suspended particulate matter (SPM) in the Great Barrier Reef (GBR), contributing to high turbidity, pollutant exposure and increasing vulnerability of coral reef to climate change. However, the terrestrial versus marine origin of the SPM in the GBR is uncertain. Here we provide multiple lines of evidence (13C NMR, isotopic and genetic fingerprints) to unravel that a considerable proportion of the terrestrially-derived SPM is degraded in the riverine and estuarine mixing zones before it is transported further offshore. The fingerprints of SPM in the marine environment were completely different from those of terrestrial origin but more consistent with that formed by marine phytoplankton. This result indicates that the SPM in the GBR may not have terrestrial origin but produced locally in the marine environment, which has significant implications on developing better-targeted management practices for improving water quality in the GBR.

Multidimensional specialization and generalization are pervasive in soil prokaryotes

Habitat specialization underpins biological processes from species distributions to speciation. However, organisms are often described as specialists or generalists based on a single niche axis, despite facing complex, multidimensional environments. Here, we analysed 236 environmental soil microbiomes across the United States and demonstrate that 90% of >1,200 prokaryotes followed one of two trajectories: specialization on all niche axes (multidimensional specialization) or generalization on all axes (multidimensional generalization). We then documented that this pervasive multidimensional specialization/generalization had many ecological and evolutionary consequences. First, multidimensional specialization and generalization are highly conserved with very few transitions between these two trajectories. Second, multidimensional generalists dominated communities because they were 73 times more abundant than specialists. Lastly, multidimensional specialists played important roles in community structure with ~220% more connections in microbiome networks. These results indicate that multidimensional generalization and specialization are evolutionarily stable with multidimensional generalists supporting larger populations and multidimensional specialists playing important roles within communities, probably stemming from their overrepresentation among pollutant detoxifiers and nutrient cyclers. Taken together, we demonstrate that the vast majority of soil prokaryotes are restricted to one of two multidimensional niche trajectories, multidimensional specialization or multidimensional generalization, which then has far-reaching consequences for evolutionary transitions, microbial dominance and community roles.

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.

Bacterial endophytes of sugar maple leaves vary more idiosyncratically than epiphytes across a large geographic area

Bacteria from the leaf surface and the leaf tissue have been attributed with several beneficial properties for their plant host. Though physically connected, the microbial ecology of these compartments has mostly been studied separately such that we lack an integrated understanding of the processes shaping their assembly. We sampled leaf epiphytes and endophytes from the same individuals of sugar maple across the northern portion of its range to evaluate if their community composition was driven by similar processes within and across populations differing in plant traits and overall abiotic environment. Leaf compartment explained most of the variation in community diversity and composition across samples. Leaf epiphytic communities were driven more by host and site characteristics than endophytic communities, whose community composition was more idiosyncratic across samples. Our results suggest a greater importance of priority effects and opportunistic colonization in driving community assembly of leaf endophytes. Understanding the comparative assembly of bacterial communities at the surface and inside plant leaves may be particularly useful for leveraging their respective potential for improving the health of plants in natural and anthropized ecosystems.

Vitamin interdependencies predicted by metagenomics-informed network analyses and validated in microbial community microcosms

Metagenomic or metabarcoding data are often used to predict microbial interactions in complex communities, but these predictions are rarely explored experimentally. Here, we use an organism abundance correlation network to investigate factors that control community organization in mine tailings-derived laboratory microbial consortia grown under dozens of conditions. The network is overlaid with metagenomic information about functional capacities to generate testable hypotheses. We develop a metric to predict the importance of each node within its local network environments relative to correlated vitamin auxotrophs, and predict that a Variovorax species is a hub as an important source of thiamine. Quantification of thiamine during the growth of Variovorax in minimal media show high levels of thiamine production, up to 100 mg/L. A few of the correlated thiamine auxotrophs are predicted to produce pantothenate, which we show is required for growth of Variovorax, supporting that a subset of vitamin-dependent interactions are mutualistic. A Cryptococcus yeast produces the B-vitamin pantothenate, and co-culturing with Variovorax leads to a 90-130-fold fitness increase for both organisms. Our study demonstrates the predictive power of metagenome-informed, microbial consortia-based network analyses for identifying microbial interactions that underpin the structure and functioning of microbial communities.

Preterm Prelabor Rupture of Membranes Linked to Vaginal Bacteriome of Pregnant Females in the Early Second Trimester: a Case-Cohort Design

Preterm prelabor rupture of membranes (PPROM) is a major cause of spontaneous preterm birth (sPTB), one of the greatest challenges facing obstetrics with complicated pathogenesis. This case-cohort study investigated the association between vaginal bacteriome of singleton pregnant females in the early second trimester and PPROM. The study included 35,255 and 180 pregnant females with PPROM as cases and term-birth without prelabor rupture of membranes (TWPROM) and term prelabor rupture of membranes (TPROM) pregnant females as controls, respectively. Using 16S rRNA sequencing, the vaginal microbiome traits were analyzed. Females with PPROM had higher alpha and beta diversity (P < 0.05) than TWPROM and TPROM. The presence of L. mulieris was associated with a decreased risk of PPROM (adjusted odds ratio [aOR] = 0.35; 95% confidence interval [CI]: 0.17-0.72) compared with TWPROM. Meanwhile, the presence of Megasphaera genus (aOR = 2.27; 95% CI: 1.09-4.70), Faecalibacterium genus (aOR = 3.29; 95% CI: 1.52-7.13), Bifidobacterium genus (aOR = 3.26; 95% CI: 1.47-7.24), Xanthomonadales genus (aOR = 2.76; 95% CI: 1.27-6.01), Gammaproteobacteria class (aOR = 2.36; 95% CI: 1.09-5.14), and Alphaproteobacteria class (aOR = 2.45; 95% CI: 1.14-5.26) was associated with an increased risk of PPROM compared with TWPROM. Our results indicated that the risk of PPROM can decrease with vaginal L. mulieris but increase with high alpha or beta diversity, and several vaginal bacteria in pregnant females may be involved in the occurrence of PPROM.

Understanding the role of graphene oxide in affecting PAHs biodegradation by microorganisms: An integrated analysis using 16SrRNA, metatranscriptomic, and metabolomic approaches

Graphene oxide (GO)-promoted microbial degradation technology is considered an important strategy to eliminate polycyclic aromatic hydrocarbons (PAHs) in the environment; however, the mechanism by which GO affects microbial degradation of PAHs has not been fully studied. Thus, this study aimed to analyze the effect of GO-microbial interaction on PAHs degradation at the microbial community structure, community gene expression, and metabolic levels using multi-omics combined technology. We treated PAHs-contaminated soil samples with different concentrations of GO and analyzed the soil samples for microbial diversity after 14 and 28 days. After a short exposure, GO reduced the diversity of soil microbial community but increased potential degrading microbial abundance, promoting PAHs biodegradation. This promotion effect was further influenced by the GO concentration. In a short period of time, GO upregulated the expression of genes involved in microbial movement (flagellar assembly), bacterial chemotaxis, two-component system, and phosphotransferase system in the soil microbial community and increased the probability of microbial contact with PAHs. Biosynthesis of amino acids and carbon metabolism of microorganisms were accelerated, thereby increasing the degradation of PAHs. With the extension of time, the degradation of PAHs stagnated, which may be due to the weakened stimulation of GO on microorganisms. The results showed that screening specific degrading microorganisms, increasing the contact area between microorganisms and PAHs, and prolonging the stimulation of GO on microorganisms were important means to improve the biodegradation efficiency of PAHs in soil. This study elucidates how GO affects microbial PAHs degradation and provides important insights for the application of GO-assisted microbial degradation technology.

Oral administration of Pinus koraiensis cone essential oil reduces rumen methane emission by altering the rumen microbial composition and functions in Korean native goat (Capra hircus coreanae)

This study aimed to investigate Pinus koraiensis cone essential oil (PEO) as a methane (CH₄) inhibitor and determine its impact on the taxonomic and functional characteristics of the rumen microbiota in goats. A total of 10 growing Korean native goats (Capra hircus coreanae, 29.9 ± 1.58 kg, male) were assigned to different dietary treatments: control (CON; basal diet without additive) and PEO (basal diet +1 g/d of PEO) by a 2 × 2 crossover design. Methane measurements were conducted every 4 consecutive days for 17-20 days using a laser CH₄ detector. Samples of rumen fluid and feces were collected during each experimental period to evaluate the biological effects and dry matter (DM) digestibility after PEO oral administration. The rumen microbiota was analyzed via 16S rRNA gene amplicon sequencing. The PEO oral administration resulted in reduced CH₄ emission (eructation CH₄/body weight^0.75, p = 0.079) without affecting DM intake; however, it lowered the total volatile fatty acids (p = 0.041), molar proportion of propionate (p = 0.075), and ammonia nitrogen (p = 0.087) in the rumen. Blood metabolites (i.e., albumin, alanine transaminase/serum glutamic pyruvate transaminase, creatinine, and triglyceride) were significantly affected (p < 0.05) by PEO oral administration. The absolute fungal abundance (p = 0.009) was reduced by PEO oral administration, whereas ciliate protozoa, total bacteria, and methanogen abundance were not affected. The composition of rumen prokaryotic microbiota was altered by PEO oral administration with lower evenness (p = 0.054) observed for the PEO group than the CON group. Moreover, PICRUSt2 analysis revealed that the metabolic pathways of prokaryotic bacteria, such as pyruvate metabolism, were enriched in the PEO group. We also identified the Rikenellaceae RC9 gut group as the taxa potentially contributing to the enriched KEGG modules for histidine biosynthesis and pyruvate oxidation in the rumen of the PEO group using the FishTaco analysis. The entire co-occurrence networks showed that more nodes and edges were detected in the PEO group. Overall, these findings provide an understanding of how PEO oral administration affects CH₄ emission and rumen prokaryotic microbiota composition and function. This study may help develop potential manipulation strategies to find new essential oils to mitigate enteric CH₄ emissions from ruminants.

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.

Modified Dingchuan Decoction treats cough-variant asthma by suppressing lung inflammation and regulating the lung microbiota

ETHNOPHARMACOLOGICAL RELEVANCE

Modified Dingchuan Decoction (MDD) is a Chinese medicine formula containing 11 materials with cough suppression, asthma relief, and anti-inflammatory effects.

AIM OF THE STUDY

This study aimed to evaluate the therapeutic effect of MDD on cough-variant asthma (CVA) and to investigate its mechanism of action.

MATERIALS AND METHODS

The chemical constituents of MDD were analyzed by ultra-performance liquid chromatography-quadrupole/electrostatic field orbitrap high-resolution mass spectrometry (UPLC-Q-Orbitrap HRMS). A guinea pig CVA model was established using an intramuscular injection of ovalbumin (OVA), combined with an intraperitoneal injection of aluminum hydroxide [Al(OH)₃] and nebulized OVA. At the beginning of day 18, the low, medium, and high MDD groups were gavaged with 7.23 g/kg, 14.46 g/kg, and 28.92 g/kg of MDD, respectively, and the positive group was gavaged with 5 mg/kg of prednisone acetate combined with 1 mg/kg of montelukast sodium; the normal and model groups were given an equal volume of distilled water, once a day for 21 days. The cough was induced by 10^-3 mol/L capsaicin solution 1 h after the last administration, and the number of coughs and the latency of coughs were evaluated. Hematoxylin and eosin staining (H&E) was used to observe pathological changes in the lungs and airways. The concentration of inflammatory factors in bronchoalveolar lavage fluid (BALF) was measured by enzyme-linked immunosorbent assay (ELISA). We analyzed the lung microbiota using 16 S ribosomal DNA (16 S rDNA) high-throughput sequencing.

RESULTS

The 38 chemical components were found in MDD, and MDD reduced the number of coughs in guinea pigs with CVA, prolonged cough latency, improved pathological damage to the lungs and airways, regulated inflammatory factor levels in BALF, and modulated the lung microbiota.

CONCLUSIONS

This study demonstrated that treating CVA with MDD may be related to inhibiting lung inflammation and regulating lung microbiota.

Congeneric temperate orchids recruit similar-yet differentially abundant-endophytic bacterial communities that are uncoupled from soil, but linked to host phenology and population size

PREMISE

Besides the beneficial plant-fungus symbiosis in mycorrhizal plants, bacteria also enhance plant fitness via tripartite interactions. While bacterial associations are presumably just as important for the obligate mycorrhizal family Orchidaceae, little is known about orchid associating bacteria (OAB).

METHODS

We examined the OAB communities of two, congeneric, terrestrial orchids, Platanthera cooperi and Platanthera praeclara, which represent widely disparate North American ecosystems. We tested whether they recruit distinct OAB communities, and whether variability in OAB communities can be linked to phenology, population size, or habitat soil. Genomic DNAs from roots of seedling, vegetative, and reproductive plants and from soil were subjected to Illumina sequencing of V4 and V5 regions of the 16S rRNA gene.

RESULTS

We obtained 809 OAB Zero-radius Operational Taxonomic Units (ZOTUs). Despite an overlap of 209 ZOTUs that accounted for >75% relative abundances of their respective OAB communities, the overall community structures of the two orchids were distinct. Within each orchid, distinctions were detected in the OAB communities of large and small populations and the three phenological stages. The OAB ZOTUs were either absent or present with low abundances in soil associated with both orchids.

CONCLUSIONS

The two orchids exhibited preferential recruitment of known growth-promoting OAB communities from soil. Their OAB communities also showed considerable overlap despite the large environmental and geographical separation of the two host taxa. Our results lend further support to the emerging evidence that not only the fungi, but root-associated bacteria also have functional importance for orchid ecology.

Plant domestication shapes rhizosphere microbiome assembly and metabolic functions

BACKGROUND

The rhizosphere microbiome, which is shaped by host genotypes, root exudates, and plant domestication, is crucial for sustaining agricultural plant growth. Despite its importance, how plant domestication builds up specific rhizosphere microbiomes and metabolic functions, as well as the importance of these affected rhizobiomes and relevant root exudates in maintaining plant growth, is not well understood. Here, we firstly investigated the rhizosphere bacterial and fungal communities of domestication and wild accessions of tetraploid wheat using amplicon sequencing (16S and ITS) after 9 years of domestication process at the main production sites in China. We then explored the ecological roles of root exudation in shaping rhizosphere microbiome functions by integrating metagenomics and metabolic genomics approaches. Furthermore, we established evident linkages between root morphology traits and keystone taxa based on microbial culture and plant inoculation experiments.

RESULTS

Our results suggested that plant rhizosphere microbiomes were co-shaped by both host genotypes and domestication status. The wheat genomes contributed more variation in the microbial diversity and composition of rhizosphere bacterial communities than fungal communities, whereas plant domestication status exerted much stronger influences on the fungal communities. In terms of microbial interkingdom association networks, domestication destabilized microbial network and depleted the abundance of keystone fungal taxa. Moreover, we found that domestication shifted the rhizosphere microbiome from slow growing and fungi dominated to fast growing and bacteria dominated, thereby resulting in a shift from fungi-dominated membership with enrichment of carbon fixation genes to bacteria-dominated membership with enrichment of carbon degradation genes. Metagenomics analyses further indicated that wild cultivars of wheat possess higher microbial function diversity than domesticated cultivars. Notably, we found that wild cultivar is able to harness rhizosphere microorganism carrying N transformation (i.e., nitrification, denitrification) and P mineralization pathway, whereas rhizobiomes carrying inorganic N fixation, organic N ammonification, and inorganic P solubilization genes are recruited by the releasing of root exudates from domesticated wheat. More importantly, our metabolite-wide association study indicated that the contrasting functional roles of root exudates and the harnessed keystone microbial taxa with different nutrient acquisition strategies jointly determined the aboveground plant phenotypes. Furthermore, we observed that although domesticated and wild wheats recruited distinct microbial taxa and relevant functions, domestication-induced recruitment of keystone taxa led to a consistent growth regulation of root regardless of wheat domestication status.

CONCLUSIONS

Our results indicate that plant domestication profoundly influences rhizosphere microbiome assembly and metabolic functions and provide evidence that host plants are able to harness a differentiated ecological role of root-associated keystone microbiomes through the release of root exudates to sustain belowground multi-nutrient cycles and plant growth. These findings provide valuable insights into the mechanisms underlying plant-microbiome interactions and how to harness the rhizosphere microbiome for crop improvement in sustainable agriculture. Video Abstract.

Univariate and Multivariate Statistical Analysis of Microbiome Data: An Overview

Microbiome data is high dimensional, sparse, compositional, and over-dispersed. Therefore, modeling microbiome data is very challenging and it is an active research area. Microbiome analysis has become a progressing area of research as microorganisms constitute a large part of life. Since many methods of microbiome data analysis have been presented, this review summarizes the challenges, methods used, and the advantages and disadvantages of those methods, to serve as an updated guide for those in the field. This review also compared different methods of analysis to progress the development of newer methods.

Exploration of lung mycobiome in the patients with non-small-cell lung cancer

As the Human Microbiome Project (HMP) progresses, the relationship between microbes and human health has been receiving increasing attention. A growing number of reports support the correlation between cancer and microbes. However, most studies have focused on bacteria, rather than fungal communities. In this study, we studied the alteration in lung mycobiome in patients with non-small-cell lung cancer (NSCLC) using metagenomic sequencing and qPCR. The higher fungal diversity and more complex network were observed in the patients with NSCLC. In addition, Alternaria arborescens was found as the most relevant fungus to NSCLC, and the enrichment of it in cancerous tissue was also detected. This study proposes that the changes in fungal communities may be closely related to lung cancer, and provides insights into further exploration the relationship between lung cancer and fungi.

Clinical NEC prevention practices drive different microbiome profiles and functional responses in the preterm intestine

Preterm infants with very low birthweight are at serious risk for necrotizing enterocolitis. To functionally analyse the principles of three successful preventive NEC regimens, we characterize fecal samples of 55 infants (<1500 g, n = 383, female = 22) longitudinally (two weeks) with respect to gut microbiome profiles (bacteria, archaea, fungi, viruses; targeted 16S rRNA gene sequencing and shotgun metagenomics), microbial function, virulence factors, antibiotic resistances and metabolic profiles, including human milk oligosaccharides (HMOs) and short-chain fatty acids (German Registry of Clinical Trials, No.: DRKS00009290). Regimens including probiotic Bifidobacterium longum subsp. infantis NCDO 2203 supplementation affect microbiome development globally, pointing toward the genomic potential to convert HMOs. Engraftment of NCDO 2203 is associated with a substantial reduction of microbiome-associated antibiotic resistance as compared to regimens using probiotic Lactobacillus rhamnosus LCR 35 or no supplementation. Crucially, the beneficial effects of Bifidobacterium longum subsp. infantis NCDO 2203 supplementation depends on simultaneous feeding with HMOs. We demonstrate that preventive regimens have the highest impact on development and maturation of the gastrointestinal microbiome, enabling the establishment of a resilient microbial ecosystem that reduces pathogenic threats in at-risk preterm infants.

Guild-Level Microbiome Signature Associated with COVID-19 Severity and Prognosis

Coronavirus disease 2019 (COVID-19) severity has been associated with alterations of the gut microbiota. However, the relationship between gut microbiome alterations and COVID-19 prognosis remains elusive. Here, we performed a genome-resolved metagenomic analysis on fecal samples from 300 in-hospital COVID-19 patients, collected at the time of admission. Among the 2,568 high quality metagenome-assembled genomes (HQMAGs), redundancy analysis identified 33 HQMAGs which showed differential distribution among mild, moderate, and severe/critical severity groups. Co-abundance network analysis determined that the 33 HQMAGs were organized as two competing guilds. Guild 1 harbored more genes for short-chain fatty acid biosynthesis, and fewer genes for virulence and antibiotic resistance, compared with Guild 2. Based on average abundance difference between the two guilds, the guild-level microbiome index (GMI) classified patients from different severity groups (average AUROC [area under the receiver operating curve] = 0.83). Moreover, age-adjusted partial Spearman's correlation showed that GMIs at admission were correlated with 8 clinical parameters, which are predictors for COVID-19 prognosis, on day 7 in hospital. In addition, GMI at admission was associated with death/discharge outcome of the critical patients. We further validated that GMI was able to consistently classify patients with different COVID-19 symptom severities in different countries and differentiated COVID-19 patients from healthy subjects and pneumonia controls in four independent data sets. Thus, this genome-based guild-level signature may facilitate early identification of hospitalized COVID-19 patients with high risk of more severe outcomes at time of admission. IMPORTANCE Previous reports on the associations between COVID-19 and gut microbiome have been constrained by taxonomic-level analysis and overlook the interaction between microbes. By applying a genome-resolved, reference-free, guild-based metagenomic analysis, we demonstrated that the relationship between gut microbiota and COVID-19 is genome-specific instead of taxon-specific or even species-specific. Moreover, the COVID-19-associated genomes were not independent but formed two competing guilds, with Guild 1 potentially beneficial and Guild 2 potentially more detrimental to the host based on comparative genomic analysis. The dominance of Guild 2 over Guild 1 at time of admission was associated with hospitalized COVID-19 patients at high risk for more severe outcomes. Moreover, the guild-level microbiome signature is not only correlated with the symptom severity of COVID-19 patients, but also differentiates COVID-19 patients from pneumonia controls and healthy subjects across different studies. Here, we showed the possibility of using genome-resolved and guild-level microbiome signatures to identify hospitalized COVID-19 patients with a high risk of more severe outcomes at the time of admission.

Unravelling microalgal-bacterial interactions in aquatic ecosystems through 16S rRNA gene-based co-occurrence networks

Interactions between microalgae and bacteria can directly influence the global biogeochemical cycles but the majority of such interactions remain unknown. 16S rRNA gene-based co-occurrence networks have potential to help identify microalgal-bacterial interactions. Here, we used data from 10 Earth microbiome projects to identify potential microalgal-bacterial associations in aquatic ecosystems. A high degree of clustering was observed in microalgal-bacterial modules, indicating densely connected neighbourhoods. Proteobacteria and Bacteroidetes predominantly co-occurred with microalgae and represented hubs of most modules. Our results also indicated that species-specificity may be a global characteristic of microalgal associated microbiomes. Several previously known associations were recovered from our network modules, validating that biologically meaningful results can be inferred using this approach. A range of previously unknown associations were recognised such as co-occurrences of Bacillariophyta with uncultured Planctomycetes OM190 and Deltaproteobacteria order NB1-j. Planctomycetes and Verrucomicrobia were identified as key associates of microalgae due to their frequent co-occurrences with several microalgal taxa. Despite no clear taxonomic pattern, bacterial associates appeared functionally similar across different environments. To summarise, we demonstrated the potential of 16S rRNA gene-based co-occurrence networks as a hypothesis-generating framework to guide more focused research on microalgal-bacterial associations.

Impact of diet and host genetics on the murine intestinal mycobiome

The mammalian gut is home to a diverse microbial ecosystem, whose composition affects various physiological traits of the host. Next-generation sequencing-based metagenomic approaches demonstrated how the interplay of host genetics, bacteria, and environmental factors shape complex traits and clinical outcomes. However, the role of fungi in these complex interactions remains understudied. Here, using 228 males and 363 females from an advanced-intercross mouse line, we provide evidence that fungi are regulated by host genetics. In addition, we map quantitative trait loci associated with various fungal species to single genes in mice using whole genome sequencing and genotyping. Moreover, we show that diet and its' interaction with host genetics alter the composition of fungi in outbred mice, and identify fungal indicator species associated with different dietary regimes. Collectively, in this work, we uncover an association of the intestinal fungal community with host genetics and a regulatory role of diet in this ecological niche.

Network-based quantitative trait linkage analysis of microbiome composition in inflammatory bowel disease families

Introduction: Inflammatory bowel disease (IBD) is characterized by a dysbiosis of the gut microbiome that results from the interaction of the constituting taxa with one another, and with the host. At the same time, host genetic variation is associated with both IBD risk and microbiome composition. Methods: In the present study, we defined quantitative traits (QTs) from modules identified in microbial co-occurrence networks to measure the inter-individual consistency of microbial abundance and subjected these QTs to a genome-wide quantitative trait locus (QTL) linkage analysis. Results: Four microbial network modules were consistently identified in two cohorts of healthy individuals, but three of the corresponding QTs differed significantly between IBD patients and unaffected individuals. The QTL linkage analysis was performed in a sub-sample of the Kiel IBD family cohort (IBD-KC), an ongoing study of 256 German families comprising 455 IBD patients and 575 first- and second-degree, non-affected relatives. The analysis revealed five chromosomal regions linked to one of three microbial module QTs, namely on chromosomes 3 (spanning 10.79 cM) and 11 (6.69 cM) for the first module, chr9 (0.13 cM) and chr16 (1.20 cM) for the second module, and chr13 (19.98 cM) for the third module. None of these loci have been implicated in a microbial phenotype before. Discussion: Our study illustrates the benefit of combining network and family-based linkage analysis to identify novel genetic drivers of microbiome composition in a specific disease context.

Obese dogs exhibit different fecal microbiome and specific microbial networks compared with normal weight dogs

Canine obesity is a major health concern that predisposes dogs to various disorders. The fecal microbiome has been attracting attention because of their impact on energy efficiency and metabolic disorders of host. However, little is known about specific microbial interactions, and how these may be affected by obesity in dogs. The objective of this study was to investigate the differences in fecal microbiome and specific microbial networks between obese and normal dogs. A total of 20 beagle dogs (males = 12, body weight [BW]: 10.5 ± 1.08 kg; females = 8, BW: 11.3 ± 1.71 kg; all 2-year-old) were fed to meet the maintenance energy requirements for 18 weeks. Then, 12 beagle dogs were selected based on body condition score (BCS) and divided into two groups: high BCS group (HBCS; BCS range: 7-9, males = 4, females = 2) and normal BCS group (NBCS; BCS range: 4-6, males = 4, females = 2). In the final week of the experiment, fecal samples were collected directly from the rectum, before breakfast, for analyzing the fecal microbiome using 16S rRNA gene amplicon sequencing. The HBCS group had a significantly higher final BW than the NBCS group (P < 0.01). The relative abundances of Faecalibacterium, Phascolarctobacterium, Megamonas, Bacteroides, Mucispirillum, and an unclassified genus within Ruminococcaceae were significantly higher in the HBCS group than those in the NBCS group (P < 0.05). Furthermore, some Kyoto Encyclopedia of Genes and Genomes (KEGG) modules related to amino acid biosynthesis and B vitamins biosynthesis were enriched in the HBCS group (P < 0.10), whereas those related to carbohydrate metabolism were enriched in the NBCS group (P < 0.10). Microbial network analysis revealed distinct co-occurrence and mutually exclusive interactions between the HBCS and NBCS groups. In conclusion, several genera related to short-chain fatty acid production were enriched in the HBCS group. The enriched KEGG modules in the HBCS group enhanced energy efficiency through cross-feeding between auxotrophs and prototrophs. However, further studies are needed to investigate how specific networks can be interpreted in the context of fermentation characteristics in the lower gut and obesity in dogs.

Factors driving the halophyte rhizosphere bacterial communities in coastal salt marshes

Introduction

Root-associated microorganisms promote plant growth and provide protection from stresses. Halophytes are the fundamental components maintaining ecosystem functions of coastal salt marshes; however, it is not clear how their microbiome are structured across large spatial scales. Here, we investigated the rhizosphere bacterial communities of typical coastal halophyte species (Phragmites australis and Suaeda salsa) in temperate and subtropical salt marshes across 1,100 km in eastern China.

Methods

The sampling sites were located from 30.33 to 40.90°N and 119.24 to 121.79°E across east China. A total of 36 plots were investigated in the Liaohe River Estuary, the Yellow River Estuary, Yancheng, and Hangzhou Bay in August 2020. We collected shoot, root, and rhizosphere soil samples. the number of pakchoi leaves, total fresh and dry weight of the seedlings was counted. The soil properties, plant functional traits, the genome sequencing, and metabolomics assay were detected.

Results

The results showed that soil nutrients (total organic carbon, dissolved organic carbon, total nitrogen, soluble sugars, and organic acids) are high in the temperate marsh, while root exudates (measured by metabolite expressions) are significantly higher in the subtropical marsh. We observed higher bacterial alpha diversity, more complex network structure, and more negative connections in the temperate salt marsh, which suggested intense competition among bacterial groups. Variation partitioning analysis showed that climatic, edaphic, and root exudates had the greatest effects on the bacteria in the salt marsh, especially for abundant and moderate subcommunities. Random forest modeling further confirmed this but showed that plant species had a limited effect.

Conclutions

Taken together, the results of this study revealed soil properties (chemical properties) and root exudates (metabolites) had the greatest influence on the bacterial community of salt marsh, especially for abundant and moderate taxa. Our results provided novel insights into the biogeography of halophyte microbiome in coastal wetlands and can be beneficial for policymakers in decision-making on the management of coastal wetlands.

SCNIC: Sparse correlation network investigation for compositional data

Microbiome studies are often limited by a lack of statistical power due to small sample sizes and a large number of features. This problem is exacerbated in correlative studies of multi-omic datasets. Statistical power can be increased by finding and summarizing modules of correlated observations, which is one dimensionality reduction method. Additionally, modules provide biological insight as correlated groups of microbes can have relationships among themselves. To address these challenges, we developed SCNIC: Sparse Cooccurrence Network Investigation for compositional data. SCNIC is open-source software that can generate correlation networks and detect and summarize modules of highly correlated features. Modules can be formed using either the Louvain Modularity Maximization (LMM) algorithm or a Shared Minimum Distance algorithm (SMD) that we newly describe here and relate to LMM using simulated data. We applied SCNIC to two published datasets and we achieved increased statistical power and identified microbes that not only differed across groups, but also correlated strongly with each other, suggesting shared environmental drivers or cooperative relationships among them. SCNIC provides an easy way to generate correlation networks, identify modules of correlated features and summarize them for downstream statistical analysis. Although SCNIC was designed considering properties of microbiome data, such as compositionality and sparsity, it can be applied to a variety of data types including metabolomics data and used to integrate multiple data types. SCNIC allows for the identification of functional microbial relationships at scale while increasing statistical power through feature reduction.

Comparative metagenomic analysis reveals rhizosphere microbial community composition and functions help protect grapevines against salt stress

Introduction

Soil salinization is a serious abiotic stress for grapevines. The rhizosphere microbiota of plants can help counter the negative effects caused by salt stress, but the distinction between rhizosphere microbes of salt-tolerant and salt-sensitive varieties remains unclear.

Methods

This study employed metagenomic sequencing to explore the rhizosphere microbial community of grapevine rootstocks 101-14 (salt tolerant) and 5BB (salt sensitive) with or without salt stress.

Results and Discussion

Compared to the control (treated with ddH₂O), salt stress induced greater changes in the rhizosphere microbiota of 101-14 than in that of 5BB. The relative abundances of more plant growth-promoting bacteria, including Planctomycetes, Bacteroidetes, Verrucomicrobia, Cyanobacteria, Gemmatimonadetes, Chloroflexi, and Firmicutes, were increased in 101-14 under salt stress, whereas only the relative abundances of four phyla (Actinobacteria, Gemmatimonadetes, Chloroflexi, and Cyanobacteria) were increased in 5BB under salt stress while those of three phyla (Acidobacteria, Verrucomicrobia, and Firmicutes) were depleted. The differentially enriched functions (KEGG level 2) in 101-14 were mainly associated with pathways related to cell motility; folding, sorting, and degradation functions; glycan biosynthesis and metabolism; xenobiotics biodegradation and metabolism; and metabolism of cofactors and vitamins, whereas only the translation function was differentially enriched in 5BB. Under salt stress, the rhizosphere microbiota functions of 101-14 and 5BB differed greatly, especially pathways related to metabolism. Further analysis revealed that pathways associated with sulfur and glutathione metabolism as well as bacterial chemotaxis were uniquely enriched in 101-14 under salt stress and therefore might play vital roles in the mitigation of salt stress on grapevines. In addition, the abundance of various sulfur cycle-related genes, including genes involved in assimilatory sulfate reduction (cysNC, cysQ, sat, and sir), sulfur reduction (fsr), SOX systems (soxB), sulfur oxidation (sqr), organic sulfur transformation (tpa, mdh, gdh, and betC), increased significantly in 101-14 after treatment with NaCl; these genes might mitigate the harmful effects of salt on grapevine. In short, the study findings indicate that both the composition and functions of the rhizosphere microbial community contribute to the enhanced tolerance of some grapevines to salt stress.

Differential Effect of Vaginal Microbiota on Spontaneous Preterm Birth among Chinese Pregnant Women

Objective

The effect of vaginal microbiota on spontaneous preterm birth (sPTB) has not been fully addressed, and few studies have explored the associations between vaginal taxa and sPTB in the gestational diabetes mellitus (GDM) and non-GDM groups, respectively. Study Design. To minimize external interference, a total of 41 pregnant women with sPTB and 308 controls (pregnant women without sPTB) from same regain were enrolled in this case-cohort study. Controls were randomly selected at baseline. With the exception of GDM, other characteristics were not significantly different between the two groups. Vaginal swabs were collected at early second trimester. Using 16S amplicon sequencing, the main bioinformatics analysis was performed on the platform of QIIME 2. Vaginal microbiota traits of the sPTB group were compared with controls. Finally, the effects of binary taxa on sPTB in the GDM group and the non-GDM group were analyzed, respectively.

Results

The proportion of GDM in the sPTB (19.51%) was higher than the controls (7.47%, P = 0.018). The vaginal microbiota of pregnant women with sPTB exhibited higher alpha diversity metrics (observed features, P = 0.001; Faith's phylogenetic diversity, P = 0.013) and different beta diversity metrics (unweighted UniFrac, P = 0.006; Jaccard's distance, P = 0.004), compared with controls. The presence of Lactobacillus paragasseri/gasseri (aOR: 3.12, 95% CI: 1.24-7.84), Streptococcus (aOR: 3.58, 95% CI: 1.68-7.65), or Proteobacteria (aOR: 3.39, 95% CI: 1.55-7.39) was associated with an increased risk of sPTB in the non-GDM group (P < 0.05). However, the relative abundance of novel L. mulieris (a new species of the L. delbrueckii group) was associated with a decreased risk of sPTB (false discovery rate, 0.10) in all pregnant women.

Conclusion

GDM may modify the association of vaginal taxa with sPTB, suggesting that maternal GDM should be considered when using vaginal taxa to identify pregnant women at high risk of sPTB.

Strong pathogen competition in neonatal gut colonisation

Opportunistic bacterial pathogen species and their strains that colonise the human gut are generally understood to compete against both each other and the commensal species colonising this ecosystem. Currently we are lacking a population-wide quantification of strain-level colonisation dynamics and the relationship of colonisation potential to prevalence in disease, and how ecological factors might be modulating these. Here, using a combination of latest high-resolution metagenomics and strain-level genomic epidemiology methods we performed a characterisation of the competition and colonisation dynamics for a longitudinal cohort of neonatal gut microbiomes. We found strong inter- and intra-species competition dynamics in the gut colonisation process, but also a number of synergistic relationships among several species belonging to genus Klebsiella, which includes the prominent human pathogen Klebsiella pneumoniae. No evidence of preferential colonisation by hospital-adapted pathogen lineages in either vaginal or caesarean section birth groups was detected. Our analysis further enabled unbiased assessment of strain-level colonisation potential of extra-intestinal pathogenic Escherichia coli (ExPEC) in comparison with their propensity to cause bloodstream infections. Our study highlights the importance of systematic surveillance of bacterial gut pathogens, not only from disease but also from carriage state, to better inform therapies and preventive medicine in the future.

Red seaweed extracts reduce methane production by altering rumen fermentation and microbial composition in vitro

A series of in vitro batch culture incubations were carried out to investigate changes in rumen fermentation characteristics, methane (CH4) production, and microbial composition in response to supplementation with five different red seaweed species (Amphiroa anceps, AANC; Asparagopsis taxiformis, ATAX; Chondracanthus tenellus, CTEN; Grateloupia elliptica, GELL; and Gracilaria parvispora, GPAR). Prior to the incubations, the total flavonoid and polyphenol content of the red seaweed extracts was quantified. The incubated substrate consisted of timothy hay and corn grain [60:40 dry matter (DM) basis]. Treatments were substrate mixtures without seaweed extract (CON) or substrate mixtures supplemented with 0.25 mg/mL of red seaweed extract. Samples were incubated for 6, 12, 24, 36, and 48 h. Each sample was incubated in triplicates in three separate runs. In vitro DM degradability, fermentation parameters (i.e., pH, volatile fatty acids, and ammonia nitrogen), total gas production, and CH4 production were analyzed for all time points. Microbial composition was analyzed using 16S rRNA amplicon sequencing after 24 h of incubation. The highest CH4 reduction (mL/g DM, mL/g digested DM, and % of total gas production) was observed in ATAX (51.3, 50.1, and 51.5%, respectively, compared to CON; P < 0.001) after 12 h of incubation. The other red seaweed extracts reduced the CH4 production (mL/g DM; P < 0.001) in the range of 4.6-35.0% compared to CON after 24 h of incubation. After 24 h of incubation, supplementation with red seaweed extracts tended to increase the molar proportion of propionate (P = 0.057) and decreased the acetate to propionate ratio (P = 0.033) compared to the CON. Abundances of the genus Methanobrevibacter and total methanogens were reduced (P = 0.050 and P = 0.016) by red seaweed extract supplementation. The linear discriminant analysis effect size (P < 0.05, LDA ≥ 2.0) showed that UG Succinivibrionaceae, Anaeroplasma, and UG Ruminococcaceae, which are associated with higher propionate production, starch degradation, and amylase activity were relatively more abundant in red seaweed extracts than in the CON. Our results suggest that supplementation with red seaweed extracts altered the microbiota, leading to the acceleration of propionate production and reduction in CH4 production.

A distinct clade of Bifidobacterium longum in the gut of Bangladeshi children thrives during weaning

The gut microbiome has an important role in infant health and development. We characterized the fecal microbiome and metabolome of 222 young children in Dhaka, Bangladesh during the first two years of life. A distinct Bifidobacterium longum clade expanded with introduction of solid foods and harbored enzymes for utilizing both breast milk and solid food substrates. The clade was highly prevalent in Bangladesh, present globally (at lower prevalence), and correlated with many other gut taxa and metabolites, indicating an important role in gut ecology. We also found that the B. longum clades and associated metabolites were implicated in childhood diarrhea and early growth, including positive associations between growth measures and B. longum subsp. infantis, indolelactate and N-acetylglutamate. Our data demonstrate geographic, cultural, seasonal, and ecological heterogeneity that should be accounted for when identifying microbiome factors implicated in and potentially benefiting infant development.

Genome-centric analysis of short and long read metagenomes reveals uncharacterized microbiome diversity in Southeast Asians

Despite extensive efforts to address it, the vastness of uncharacterized 'dark matter' microbial genetic diversity can impact short-read sequencing based metagenomic studies. Population-specific biases in genomic reference databases can further compound this problem. Leveraging advances in hybrid assembly (using short and long reads) and Hi-C technologies in a cross-sectional survey, we deeply characterized 109 gut microbiomes from three ethnicities in Singapore to comprehensively reconstruct 4497 medium and high-quality metagenome assembled genomes, 1708 of which were missing in short-read only analysis and with >28× N50 improvement. Species-level clustering identified 70 (>10% of total) novel gut species out of 685, improved reference genomes for 363 species (53% of total), and discovered 3413 strains unique to these populations. Among the top 10 most abundant gut bacteria in our study, one of the species and >80% of strains were unrepresented in existing databases. Annotation of biosynthetic gene clusters (BGCs) uncovered more than 27,000 BGCs with a large fraction (36-88%) unrepresented in current databases, and with several unique clusters predicted to produce bacteriocins that could significantly alter microbiome community structure. These results reveal significant uncharacterized gut microbial diversity in Southeast Asian populations and highlight the utility of hybrid metagenomic references for bioprospecting and disease-focused studies.