ATLAS: a Snakemake workflow for assembly, annotation, and genomic binning of metagenome sequence data

An orally active carbon monoxide-releasing molecule enhances beneficial gut microbial species to combat obesity in mice

Carbon monoxide (CO), a gaseous signaling molecule, has shown promise in preventing body weight gain and metabolic dysfunction induced by high fat diet (HFD), but the mechanisms underlying these effects are largely unknown. An essential component in response to HFD is the gut microbiome, which is significantly altered during obesity and represents a target for developing new therapeutic interventions to fight metabolic diseases. Here, we show that CO delivered to the gut by oral administration with a CO-releasing molecule (CORM-401) accumulates in faeces and enriches a variety of microbial species that were perturbed by a HFD regimen. Notably, Akkermansia muciniphila, which exerts salutary metabolic effects in mice and humans, was strongly depleted by HFD but was the most abundant gut species detected after CORM-401 treatment. Analysis of bacterial transcripts revealed a restoration of microbial functional activity, with partial or full recovery of the Krebs cycle, β-oxidation, respiratory chain and glycolysis. Mice treated with CORM-401 exhibited normalization of several plasma and fecal metabolites that were disrupted by HFD and are dependent on Akkermansia muciniphila's metabolic activity, including indoles and tryptophan derivatives. Finally, CORM-401 treatment led to an improvement in gut morphology as well as reduction of inflammatory markers in colon and cecum and restoration of metabolic profiles in these tissues. Our findings provide therapeutic insights on the efficacy of CO as a potential prebiotic to combat obesity, identifying the gut microbiota as a crucial target for CO-mediated pharmacological activities against metabolic disorders.

Microbial genome collection of aerobic granular sludge cultivated in sequencing batch reactors using different carbon source mixtures

Aerobic granular sludge (AGS) consists of a microbial consortium that has an important role in wastewater treatment. This study investigates AGS microorganisms cultivated in a laboratory-scale sequencing batch reactor. Metagenomic sequencing was conducted using PacBio and Illumina, resulting in 759 metagenome-assembled genomes, 331 of which remained after dereplication.

A seven-step guide to spatial, agent-based modelling of tumour evolution

Spatial agent-based models are frequently used to investigate the evolution of solid tumours subject to localized cell-cell interactions and microenvironmental heterogeneity. As spatial genomic, transcriptomic and proteomic technologies gain traction, spatial computational models are predicted to become ever more necessary for making sense of complex clinical and experimental data sets, for predicting clinical outcomes, and for optimizing treatment strategies. Here we present a non-technical step by step guide to developing such a model from first principles. Stressing the importance of tailoring the model structure to that of the biological system, we describe methods of increasing complexity, from the basic Eden growth model up to off-lattice simulations with diffusible factors. We examine choices that unavoidably arise in model design, such as implementation, parameterization, visualization and reproducibility. Each topic is illustrated with examples drawn from recent research studies and state of the art modelling platforms. We emphasize the benefits of simpler models that aim to match the complexity of the phenomena of interest, rather than that of the entire biological system. Our guide is aimed at both aspiring modellers and other biologists and oncologists who wish to understand the assumptions and limitations of the models on which major cancer studies now so often depend.

The potential of enhanced phytoremediation to clean up multi-contaminated soil - insights from metatranscriptomics

This study aimed to (i) investigate the potential for enhanced phytoremediation to remove contaminants from soil historically co-contaminated with petroleum hydrocarbons (PHs) and heavy metals (HMs) and (ii) analyze the expression of crucial bacterial genes and whole metatranscriptomics profiles for better understanding of soil processes during applied treatment. Phytoremediation was performed using Zea mays and supported by the Pseudomonas qingdaonensis ZCR6 strain and a natural biofertilizer: meat and bone meal (MBM). In previous investigations, mechanisms supporting plant growth and PH degradation were described in the ZCR6 strain. Here, ZCR6 survived in the soil throughout the experiment, but the efficacy of PH removal from all soils fertilized with MBM reached 32 % regardless of the bacterial inoculation. All experimental groups contained 2 % (w/w) MBM. The toxic effect of this amendment on plants was detected 30 days after germination, irrespective of ZCR6 inoculation. Among the 17 genes tested using the qPCR method, only expression of the acdS gene, encoding 1-aminocyclopropane-1-carboxylic acid deaminase, and the CYP153 gene, encoding cytochrome P450-type alkane hydroxylase, was detected in soils. Metatranscriptomic analysis of soils indicated increased expression of methane particulated ammonia monooxygenase subunit A (pmoA-amoA) by Nitrosomonadales bacteria in all soils enriched with MBM compared to the non-fertilized control. We suggest that the addition of 2 % (w/w) MBM caused the toxic effect on plants via the rapid release of ammonia, and this led to high pmoA-amoA expression. In parallel, due to its wide substrate specificity, enhanced bacterial hydrocarbon removal in MBM-treated soils was observed. The metatranscriptomic results indicate that MBM application should be considered to improve bioremediation of soils polluted with PHs rather than phytoremediation. However, lower concentrations of MBM could be considered for phytoremediation enhancement. From a broader perspective, these results indicated the superior capability of metatranscriptomics to investigate the microbial mechanisms driving various bioremediation techniques.

Integrating taxonomic signals from MAGs and contigs improves read annotation and taxonomic profiling of metagenomes

Metagenomic analysis typically includes read-based taxonomic profiling, assembly, and binning of metagenome-assembled genomes (MAGs). Here we integrate these steps in Read Annotation Tool (RAT), which uses robust taxonomic signals from MAGs and contigs to enhance read annotation. RAT reconstructs taxonomic profiles with high precision and sensitivity, outperforming other state-of-the-art tools. In high-diversity groundwater samples, RAT annotates a large fraction of the metagenomic reads, calling novel taxa at the appropriate, sometimes high taxonomic ranks. Thus, RAT integrative profiling provides an accurate and comprehensive view of the microbiome from shotgun metagenomics data. The package of Contig Annotation Tool (CAT), Bin Annotation Tool (BAT), and RAT is available at https://github.com/MGXlab/CAT_pack (from CAT pack v6.0). The CAT pack now also supports Genome Taxonomy Database (GTDB) annotations.

Active microbial communities facilitate carbon turnover in brine pools found in the deep Southeastern Mediterranean Sea

Discharge of gas-rich brines fuels productive chemosynthetic ecosystems in the deep sea. In these salty, methanic and sulfidic brines, microbial communities adapt to specific niches along the physicochemical gradients. However, the molecular mechanisms that underpin these adaptations are not fully known. Using metagenomics, we investigated the dense (∼106 cell ml-1) microbial communities that occupy small deep-sea brine pools found in the Southeastern Mediterranean Sea (1150 m water depth, ∼22 °C, ∼60 PSU salinity, sulfide, methane, ammonia reaching millimolar levels, and oxygen usually depleted), reaching high productivity rates of 685 μg C L-1 d-1 ex-situ. We curated 266 metagenome-assembled genomes of bacteria and archaea from the several pools and adjacent sediment-water interface, highlighting the dominance of a single Sulfurimonas, which likely fuels its autotrophy using sulfide oxidation or inorganic sulfur disproportionation. This lineage may be dominant in its niche due to genome streamlining, limiting its metabolic repertoire, particularly by using a single variant of sulfide: quinone oxidoreductase. These primary producers co-exist with ANME-2c archaea that catalyze the anaerobic oxidation of methane. Other lineages can degrade the necromass aerobically (Halomonas and Alcanivorax), or anaerobically through fermentation of macromolecules (e.g., Caldatribacteriota, Bipolaricaulia, Chloroflexota, etc). These low-abundance organisms likely support the autotrophs, providing energy-rich H2, and vital organics such as vitamin B12.

ARGprofiler-a pipeline for large-scale analysis of antimicrobial resistance genes and their flanking regions in metagenomic datasets

MOTIVATION

Analyzing metagenomic data can be highly valuable for understanding the function and distribution of antimicrobial resistance genes (ARGs). However, there is a need for standardized and reproducible workflows to ensure the comparability of studies, as the current options involve various tools and reference databases, each designed with a specific purpose in mind.

RESULTS

In this work, we have created the workflow ARGprofiler to process large amounts of raw sequencing reads for studying the composition, distribution, and function of ARGs. ARGprofiler tackles the challenge of deciding which reference database to use by providing the PanRes database of 14 078 unique ARGs that combines several existing collections into one. Our pipeline is designed to not only produce abundance tables of genes and microbes but also to reconstruct the flanking regions of ARGs with ARGextender. ARGextender is a bioinformatic approach combining KMA and SPAdes to recruit reads for a targeted de novo assembly. While our aim is on ARGs, the pipeline also creates Mash sketches for fast searching and comparisons of sequencing runs.

AVAILABILITY AND IMPLEMENTATION

The ARGprofiler pipeline is a Snakemake workflow that supports the reuse of metagenomic sequencing data and is easily installable and maintained at https://github.com/genomicepidemiology/ARGprofiler.

Anoxygenic phototroph of the Chloroflexota uses a type I reaction centre

Scientific exploration of phototrophic bacteria over nearly 200 years has revealed large phylogenetic gaps between known phototrophic groups that limit understanding of how phototrophy evolved and diversified1,2. Here, through Boreal Shield lake water incubations, we cultivated an anoxygenic phototrophic bacterium from a previously unknown order within the Chloroflexota phylum that represents a highly novel transition form in the evolution of photosynthesis. Unlike all other known phototrophs, this bacterium uses a type I reaction centre (RCI) for light energy conversion yet belongs to the same bacterial phylum as organisms that use a type II reaction centre (RCII) for phototrophy. Using physiological, phylogenomic and environmental metatranscriptomic data, we demonstrate active RCI-utilizing metabolism by the strain alongside usage of chlorosomes3 and bacteriochlorophylls4 related to those of RCII-utilizing Chloroflexota members. Despite using different reaction centres, our phylogenomic data provide strong evidence that RCI-utilizing and RCII-utilizing Chloroflexia members inherited phototrophy from a most recent common phototrophic ancestor. The Chloroflexota phylum preserves an evolutionary record of the use of contrasting phototrophic modes among genetically related bacteria, giving new context for exploring the diversification of phototrophy on Earth.

Exploring the gut DNA virome in fecal immunochemical test stool samples reveals associations with lifestyle in a large population-based study

Stool samples for fecal immunochemical tests (FIT) are collected in large numbers worldwide as part of colorectal cancer screening programs. Employing FIT samples from 1034 CRCbiome participants, recruited from a Norwegian colorectal cancer screening study, we identify, annotate and characterize more than 18000 DNA viruses, using shotgun metagenome sequencing. Only six percent of them are assigned to a known taxonomic family, with Microviridae being the most prevalent viral family. Linking individual profiles to comprehensive lifestyle and demographic data shows 17/25 of the variables to be associated with the gut virome. Physical activity, smoking, and dietary fiber consumption exhibit strong and consistent associations with both diversity and relative abundance of individual viruses, as well as with enrichment for auxiliary metabolic genes. We demonstrate the suitability of FIT samples for virome analysis, opening an opportunity for large-scale studies of this enigmatic part of the gut microbiome. The diverse viral populations and their connections to the individual lifestyle uncovered herein paves the way for further exploration of the role of the gut virome in health and disease.

Metagenomic assembly is the main bottleneck in the identification of mobile genetic elements

Antimicrobial resistance genes (ARG) are commonly found on acquired mobile genetic elements (MGEs) such as plasmids or transposons. Understanding the spread of resistance genes associated with mobile elements (mARGs) across different hosts and environments requires linking ARGs to the existing mobile reservoir within bacterial communities. However, reconstructing mARGs in metagenomic data from diverse ecosystems poses computational challenges, including genome fragment reconstruction (assembly), high-throughput annotation of MGEs, and identification of their association with ARGs. Recently, several bioinformatics tools have been developed to identify assembled fragments of plasmids, phages, and insertion sequence (IS) elements in metagenomic data. These methods can help in understanding the dissemination of mARGs. To streamline the process of identifying mARGs in multiple samples, we combined these tools in an automated high-throughput open-source pipeline, MetaMobilePicker, that identifies ARGs associated with plasmids, IS elements and phages, starting from short metagenomic sequencing reads. This pipeline was used to identify these three elements on a simplified simulated metagenome dataset, comprising whole genome sequences from seven clinically relevant bacterial species containing 55 ARGs, nine plasmids and five phages. The results demonstrated moderate precision for the identification of plasmids (0.57) and phages (0.71), and moderate sensitivity of identification of IS elements (0.58) and ARGs (0.70). In this study, we aim to assess the main causes of this moderate performance of the MGE prediction tools in a comprehensive manner. We conducted a systematic benchmark, considering metagenomic read coverage, contig length cutoffs and investigating the performance of the classification algorithms. Our analysis revealed that the metagenomic assembly process is the primary bottleneck when linking ARGs to identified MGEs in short-read metagenomics sequencing experiments rather than ARGs and MGEs identification by the different tools.

A genus in the bacterial phylum Aquificota appears to be endemic to Aotearoa-New Zealand

Allopatric speciation has been difficult to examine among microorganisms, with prior reports of endemism restricted to sub-genus level taxa. Previous microbial community analysis via 16S rRNA gene sequencing of 925 geothermal springs from the Taupō Volcanic Zone (TVZ), Aotearoa-New Zealand, revealed widespread distribution and abundance of a single bacterial genus across 686 of these ecosystems (pH 1.2-9.6 and 17.4-99.8 °C). Here, we present evidence to suggest that this genus, Venenivibrio (phylum Aquificota), is endemic to Aotearoa-New Zealand. A specific environmental niche that increases habitat isolation was identified, with maximal read abundance of Venenivibrio occurring at pH 4-6, 50-70 °C, and low oxidation-reduction potentials. This was further highlighted by genomic and culture-based analyses of the only characterised species for the genus, Venenivibrio stagnispumantis CP.B2T, which confirmed a chemolithoautotrophic metabolism dependent on hydrogen oxidation. While similarity between Venenivibrio populations illustrated that dispersal is not limited across the TVZ, extensive amplicon, metagenomic, and phylogenomic analyses of global microbial communities from DNA sequence databases indicates Venenivibrio is geographically restricted to the Aotearoa-New Zealand archipelago. We conclude that geographic isolation, complemented by physicochemical constraints, has resulted in the establishment of an endemic bacterial genus.

Evaluation of multiple displacement amplification for metagenomic analysis of low biomass samples

Combining multiple displacement amplification (MDA) with metagenomics enables the analysis of samples with extremely low DNA concentrations, making them suitable for high-throughput sequencing. Although amplification bias and nonspecific amplification have been reported from MDA-amplified samples, the impact of MDA on metagenomic datasets is not well understood. We compared three MDA methods (i.e. bulk MDA, emulsion MDA, and primase MDA) for metagenomic analysis of two DNA template concentrations (approx. 1 and 100 pg) derived from a microbial community standard "mock community" and two low biomass environmental samples (i.e. borehole fluid and groundwater). We assessed the impact of MDA on metagenome-based community composition, assembly quality, functional profiles, and binning. We found amplification bias against high GC content genomes but relatively low nonspecific amplification such as chimeras, artifacts, or contamination for all MDA methods. We observed MDA-associated representational bias for microbial community profiles, especially for low-input DNA and with the primase MDA method. Nevertheless, similar taxa were represented in MDA-amplified libraries to those of unamplified samples. The MDA libraries were highly fragmented, but similar functional profiles to the unamplified libraries were obtained for bulk MDA and emulsion MDA at higher DNA input and across these MDA libraries for the groundwater sample. Medium to low-quality bins were possible for the high input bulk MDA metagenomes for the most simple microbial communities, borehole fluid, and mock community. Although MDA-based amplification should be avoided, it can still reveal meaningful taxonomic and functional information from samples with extremely low DNA concentration where direct metagenomics is otherwise impossible.

Taxonomic distribution of metabolic functions in bacteria associated with Trichodesmium consortia

IMPORTANCE

Colonies of the cyanobacteria Trichodesmium act as a biological hotspot for the usage and recycling of key resources such as C, N, P, and Fe within an otherwise oligotrophic environment. While Trichodesmium colonies are known to interact and support a unique community of algae and particle-associated microbes, our understanding of the taxa that populate these colonies and the gene functions they encode is still limited. Characterizing the taxa and adaptive strategies that influence consortium physiology and its concomitant biogeochemistry is critical in a future ocean predicted to have increasingly resource-depleted regions.

Amino acid auxotrophies in human gut bacteria are linked to higher microbiome diversity and long-term stability

Amino acid auxotrophies are prevalent among bacteria. They can govern ecological dynamics in microbial communities and indicate metabolic cross-feeding interactions among coexisting genotypes. Despite the ecological importance of auxotrophies, their distribution and impact on the diversity and function of the human gut microbiome remain poorly understood. This study performed the first systematic analysis of the distribution of amino acid auxotrophies in the human gut microbiome using a combined metabolomic, metagenomic, and metabolic modeling approach. Results showed that amino acid auxotrophies are ubiquitous in the colon microbiome, with tryptophan auxotrophy being the most common. Auxotrophy frequencies were higher for those amino acids that are also essential to the human host. Moreover, a higher overall abundance of auxotrophies was associated with greater microbiome diversity and stability, and the distribution of auxotrophs was found to be related to the human host's metabolome, including trimethylamine oxide, small aromatic acids, and secondary bile acids. Thus, our results suggest that amino acid auxotrophies are important factors contributing to microbiome ecology and host-microbiome metabolic interactions.

Tar patties are hotspots of hydrocarbon turnover and nitrogen fixation during a nearshore pollution event in the oligotrophic southeastern Mediterranean Sea

Weathered oil, that is, tar, forms hotspots of hydrocarbon degradation by complex biota in marine environment. Here, we used marker gene sequencing and metagenomics to characterize the communities of bacteria, archaea and eukaryotes that colonized tar patties and control samples (wood, plastic), collected in the littoral following an offshore spill in the warm, oligotrophic southeastern Mediterranean Sea (SEMS). We show potential aerobic and anaerobic hydrocarbon catabolism niches on tar interior and exterior, linking carbon, sulfur and nitrogen cycles. Alongside aromatics and larger alkanes, short-chain alkanes appear to fuel dominant populations, both the aerobic clade UBA5335 (Macondimonas), anaerobic Syntropharchaeales, and facultative Mycobacteriales. Most key organisms, including the hydrocarbon degraders and cyanobacteria, have the potential to fix dinitrogen, potentially alleviating the nitrogen limitation of hydrocarbon degradation in the SEMS. We highlight the complexity of these tar-associated communities, where bacteria, archaea and eukaryotes co-exist, likely exchanging metabolites and competing for resources and space.

Desulfoferula mesophilus gen. nov. sp. nov., a mesophilic sulfate-reducing bacterium isolated from a brackish lake sediment

A novel sulfate-reducing bacterium (strain 12FAKT) was isolated from sediment sampled from a brackish lake in Japan. Respiratory growth was observed with formate and pyruvate as an electron donor. Sulfate, thiosulfate, elemental sulfur and dimethyl sulfoxide were utilized as an electron acceptor. The isolate grew over a temperature range of 18-42 °C (optimum 35-37 °C), a NaCl concentration range of 50-450 mM (optimum 150-300 mM) and a pH range of 6.6-7.5. The 12FAKT genome consists of a circular chromosome with a length of 4.5 Mbp and G + C content of 63.6%. Based on 16S rRNA gene sequence similarity, the closest cultured relative was Desulfarculus baarsii 2st14T (92.2%). Genome-based phylogenetic analysis placed strain 12FAKT within the family Desulfarculaceae but did not affiliate the strain with any existing genus. Taken together, we propose a novel species of a novel genus, Desulfoferula mesophilus gen. nov. sp. nov. with the type strain 12FAKT (= DSM 115219T = JCM 39399T).

Extending and improving metagenomic taxonomic profiling with uncharacterized species using MetaPhlAn 4

Metagenomic assembly enables new organism discovery from microbial communities, but it can only capture few abundant organisms from most metagenomes. Here we present MetaPhlAn 4, which integrates information from metagenome assemblies and microbial isolate genomes for more comprehensive metagenomic taxonomic profiling. From a curated collection of 1.01 M prokaryotic reference and metagenome-assembled genomes, we define unique marker genes for 26,970 species-level genome bins, 4,992 of them taxonomically unidentified at the species level. MetaPhlAn 4 explains ~20% more reads in most international human gut microbiomes and >40% in less-characterized environments such as the rumen microbiome and proves more accurate than available alternatives on synthetic evaluations while also reliably quantifying organisms with no cultured isolates. Application of the method to >24,500 metagenomes highlights previously undetected species to be strong biomarkers for host conditions and lifestyles in human and mouse microbiomes and shows that even previously uncharacterized species can be genetically profiled at the resolution of single microbial strains.

Shared ancestry of algal symbiosis and chloroplast sequestration in foraminifera

Foraminifera are unicellular organisms that established the most diverse algal symbioses in the marine realm. Endosymbiosis repeatedly evolved in several lineages, while some engaged in the sequestration of chloroplasts, known as kleptoplasty. So far, kleptoplasty has been documented exclusively in the rotaliid clade. Here, we report the discovery of kleptoplasty in the species Hauerina diversa that belongs to the miliolid clade. The existence of kleptoplasty in the two main clades suggests that it is more widespread than previously documented. We observed chloroplasts in clustered structures within the foraminiferal cytoplasm and confirmed their functionality. Phylogenetic analysis of 18S ribosomal RNA gene sequences showed that H. diversa branches next to symbiont-bearing Alveolinidae. This finding represents evidence of of a relationship between kleptoplastic and symbiotic foraminifera.. Analysis of ribosomal genes and metagenomics revealed that alveolinid symbionts and kleptoplasts belong to the same clade, which suggests a common ancestry.

HoCoRT: host contamination removal tool

BACKGROUND

Shotgun metagenome sequencing data obtained from a host environment will usually be contaminated with sequences from the host organism. Host sequences should be removed before further analysis to avoid biases, reduce downstream computational load, or ensure privacy in the case of a human host. The tools that we identified, as designed specifically to perform host contamination sequence removal, were either outdated, not maintained, or complicated to use. Consequently, we have developed HoCoRT, a fast and user-friendly tool that implements several methods for optimised host sequence removal. We have evaluated the speed and accuracy of these methods.

RESULTS

HoCoRT is an open-source command-line tool for host contamination removal. It is designed to be easy to install and use, offering a one-step option for genome indexing. HoCoRT employs a variety of well-known mapping, classification, and alignment methods to classify reads. The user can select the underlying classification method and its parameters, allowing adaptation to different scenarios. Based on our investigation of various methods and parameters using synthetic human gut and oral microbiomes, and on assessment of publicly available data, we provide recommendations for typical datasets with short and long reads.

CONCLUSIONS

To decontaminate a human gut microbiome with short reads using HoCoRT, we found the optimal combination of speed and accuracy with BioBloom, Bowtie2 in end-to-end mode, and HISAT2. Kraken2 consistently demonstrated the highest speed, albeit with a trade-off in accuracy. The same applies to an oral microbiome, but here Bowtie2 was notably slower than the other tools. For long reads, the detection of human host reads is more difficult. In this case, a combination of Kraken2 and Minimap2 achieved the highest accuracy and detected 59% of human reads. In comparison to the dedicated DeconSeq tool, HoCoRT using Bowtie2 in end-to-end mode proved considerably faster and slightly more accurate. HoCoRT is available as a Bioconda package, and the source code can be accessed at https://github.com/ignasrum/hocort along with the documentation. It is released under the MIT licence and is compatible with Linux and macOS (except for the BioBloom module).

Growth rate is a dominant factor predicting the rhizosphere effect

The root microbiome is shaped by plant root activity, which selects specific microbial taxa from the surrounding soil. This influence on the microorganisms and soil chemistry in the immediate vicinity of the roots has been referred to as the rhizosphere effect. Understanding the traits that make bacteria successful in the rhizosphere is critical for developing sustainable agriculture solutions. In this study, we compared the growth rate potential, a complex trait that can be predicted from bacterial genome sequences, to functional traits encoded by proteins. We analyzed 84 paired rhizosphere- and soil-derived 16S rRNA gene amplicon datasets from 18 different plants and soil types, performed differential abundance analysis, and estimated growth rates for each bacterial genus. We found that bacteria with higher growth rate potential consistently dominated the rhizosphere, and this trend was confirmed in different bacterial phyla using genome sequences of 3270 bacterial isolates and 6707 metagenome-assembled genomes (MAGs) from 1121 plant- and soil-associated metagenomes. We then identified which functional traits were enriched in MAGs according to their niche or growth rate status. We found that predicted growth rate potential was the main feature for differentiating rhizosphere and soil bacteria in machine learning models, and we then analyzed the features that were important for achieving faster growth rates, which makes bacteria more competitive in the rhizosphere. As growth rate potential can be predicted from genomic data, this work has implications for understanding bacterial community assembly in the rhizosphere, where many uncultivated bacteria reside.

Microdiversity of the vaginal microbiome is associated with preterm birth

Preterm birth (PTB) is the leading cause of neonatal morbidity and mortality. The vaginal microbiome has been associated with PTB, yet the mechanisms underlying this association are not fully understood. Understanding microbial genetic adaptations to selective pressures, especially those related to the host, may yield insights into these associations. Here, we analyze metagenomic data from 705 vaginal samples collected during pregnancy from 40 women who delivered preterm spontaneously and 135 term controls from the Multi-Omic Microbiome Study-Pregnancy Initiative. We find that the vaginal microbiome of pregnancies that ended preterm exhibited unique genetic profiles. It was more genetically diverse at the species level, a result which we validate in an additional cohort, and harbored a higher richness and diversity of antimicrobial resistance genes, likely promoted by transduction. Interestingly, we find that Gardnerella species drove this higher genetic diversity, particularly during the first half of the pregnancy. We further present evidence that Gardnerella spp. underwent more frequent recombination and stronger purifying selection in genes involved in lipid metabolism. Overall, our population genetics analyses reveal associations between the vaginal microbiome and PTB and suggest that evolutionary processes acting on vaginal microbes may play a role in adverse pregnancy outcomes such as PTB.

Characterization of Vaginal Microbial Community Dynamics in the Pathogenesis of Incident Bacterial Vaginosis, a Pilot Study

BACKGROUND

Despite more than 60 years of research, the etiology of bacterial vaginosis (BV) remains controversial. In this pilot study, we used shotgun metagenomic sequencing to characterize vaginal microbial community changes before the development of incident BV (iBV).

METHODS

A cohort of African American women with a baseline healthy vaginal microbiome (no Amsel criteria, Nugent score 0-3 with no Gardnerella vaginalis morphotypes) were followed for 90 days with daily self-collected vaginal specimens for iBV (≥2 consecutive days of a Nugent score of 7-10). Shotgun metagenomic sequencing was performed on select vaginal specimens from 4 women, every other day for 12 days before iBV diagnosis. Sequencing data were analyzed through Kraken2 and bioBakery 3 workflows, and specimens were classified into community state types. Quantitative polymerase chain reaction was performed to compare the correlation of read counts with bacterial abundance.

RESULTS

Common BV-associated bacteria such as G. vaginalis , Prevotella bivia , and Fannyhessea vaginae were increasingly identified in the participants before iBV. Linear modeling indicated significant increases in G. vaginalis and F . vaginae relative abundance before iBV, whereas the relative abundance of Lactobacillus species declined over time. The Lactobacillus species decline correlated with the presence of Lactobacillus phages. We observed enrichment in bacterial adhesion factor genes on days before iBV. There were also significant correlations between bacterial read counts and abundances measured by quantitative polymerase chain reaction.

CONCLUSIONS

This pilot study characterizes vaginal community dynamics before iBV and identifies key bacterial taxa and mechanisms potentially involved in the pathogenesis of iBV.

Aerobic methanotrophy increases the net iron reduction in methanogenic lake sediments

In methane (CH4) generating sediments, methane oxidation coupled with iron reduction was suggested to be catalyzed by archaea and bacterial methanotrophs of the order Methylococcales. However, the co-existence of these aerobic and anaerobic microbes, the link between the processes, and the oxygen requirement for the bacterial methanotrophs have remained unclear. Here, we show how stimulation of aerobic methane oxidation at an energetically low experimental environment influences net iron reduction, accompanied by distinct microbial community changes and lipid biomarker patterns. We performed incubation experiments (between 30 and 120 days long) with methane generating lake sediments amended with 13C-labeled methane, following the additions of hematite and different oxygen levels in nitrogen headspace, and monitored methane turnover by 13C-DIC measurements. Increasing oxygen exposure (up to 1%) promoted aerobic methanotrophy, considerable net iron reduction, and the increase of microbes, such as Methylomonas, Geobacter, and Desulfuromonas, with the latter two being likely candidates for iron recycling. Amendments of 13C-labeled methanol as a potential substrate for the methanotrophs under hypoxia instead of methane indicate that this substrate primarily fuels methylotrophic methanogenesis, identified by high methane concentrations, strongly positive δ13CDIC values, and archaeal lipid stable isotope data. In contrast, the inhibition of methanogenesis by 2-bromoethanesulfonate (BES) led to increased methanol turnover, as suggested by similar 13C enrichment in DIC and high amounts of newly produced bacterial fatty acids, probably derived from heterotrophic bacteria. Our experiments show a complex link between aerobic methanotrophy and iron reduction, which indicates iron recycling as a survival mechanism for microbes under hypoxia.

Gut dysbiosis in Thai intrahepatic cholangiocarcinoma and hepatocellular carcinoma

Primary liver cancer (PLC), which includes intrahepatic cholangiocarcinoma (iCCA) and hepatocellular carcinoma (HCC), has the highest incidence of all cancer types in Thailand. Known etiological factors, such as viral hepatitis and chronic liver disease do not fully account for the country's unusually high incidence. However, the gut-liver axis, which contributes to carcinogenesis and disease progression, is influenced by the gut microbiome. To investigate this relationship, fecal matter from 44 Thai PLC patients and 76 healthy controls were subjected to whole-genome metagenomic shotgun sequencing and then analyzed by marker gene-based and assembly based methods. Results revealed greater gut microbiome heterogeneity in iCCA compared to HCC and healthy controls. Two Veillonella species were found to be more abundant in iCCA samples and could distinguish iCCA from HCC and healthy controls. Conversely, Ruminococcus gnavus was depleted in iCCA patients and could distinguish HCC from iCCA samples. High Veillonella genus counts in the iCCA group were associated with enriched amino acid biosynthesis and glycolysis pathways, while enriched phospholipid and thiamine metabolism pathways characterized the HCC group with high Blautia genus counts. These findings reveal distinct landscapes of gut dysbiosis among Thai iCCA and HCC patients and warrant further investigation as potential biomarkers.

Microdiversity of the Vaginal Microbiome is Associated with Preterm Birth

Preterm birth (PTB) is the leading cause of neonatal morbidity and mortality. The vaginal microbiome has been associated with PTB, yet the mechanisms underlying this association are not fully understood. Understanding microbial genetic adaptations to selective pressures, especially those related to the host, may yield new insights into these associations. To this end, we analyzed metagenomic data from 705 vaginal samples collected longitudinally during pregnancy from 40 women who delivered preterm spontaneously and 135 term controls from the Multi-Omic Microbiome Study-Pregnancy Initiative (MOMS-PI). We find that the vaginal microbiome of pregnancies that ended preterm exhibits unique genetic profiles. It is more genetically diverse at the species level, a result which we validate in an additional cohort, and harbors a higher richness and diversity of antimicrobial resistance genes, likely promoted by transduction. Interestingly, we find that Gardnerella species, a group of central vaginal pathobionts, are driving this higher genetic diversity, particularly during the first half of the pregnancy. We further present evidence that Gardnerella spp. undergoes more frequent recombination and stronger purifying selection in genes involved in lipid metabolism. Overall, our results reveal novel associations between the vaginal microbiome and PTB using population genetics analyses, and suggest that evolutionary processes acting on the vaginal microbiome may play a vital role in adverse pregnancy outcomes such as preterm birth.

Physiological changes during torpor favor association with Endozoicomonas endosymbionts in the urochordate Botrylloides leachii

Environmental perturbations evoke down-regulation of metabolism in some multicellular organisms, leading to dormancy, or torpor. Colonies of the urochordate Botrylloides leachii enter torpor in response to changes in seawater temperature and may survive for months as small vasculature remnants that lack feeding and reproductive organs but possess torpor-specific microbiota. Upon returning to milder conditions, the colonies rapidly restore their original morphology, cytology and functionality while harboring re-occurring microbiota, a phenomenon that has not been described in detail to date. Here we investigated the stability of B. leachii microbiome and its functionality in active and dormant colonies, using microscopy, qPCR, in situ hybridization, genomics and transcriptomics. A novel lineage of Endozoicomonas, proposed here as Candidatus Endozoicomonas endoleachii, was dominant in torpor animals (53-79% read abundance), and potentially occupied specific hemocytes found only in torpid animals. Functional analysis of the metagenome-assembled genome and genome-targeted transcriptomics revealed that Endozoicomonas can use various cellular substrates, like amino acids and sugars, potentially producing biotin and thiamine, but also expressing various features involved in autocatalytic symbiosis. Our study suggests that the microbiome can be linked to the metabolic and physiological states of the host, B. leachii, introducing a model organism for the study of symbioses during drastic physiological changes, such as torpor.

Metabolic handoffs between multiple symbionts may benefit the deep-sea bathymodioline mussels

Bathymodioline mussels rely on thiotrophic and/or methanotrophic chemosynthetic symbionts for nutrition, yet, secondary heterotrophic symbionts are often present and play an unknown role in the fitness of the organism. The bathymodioline Idas mussels that thrive in gas seeps and on sunken wood in the Mediterranean Sea and the Atlantic Ocean, host at least six symbiont lineages that often co-occur. These lineages include the primary symbionts chemosynthetic methane- and sulfur-oxidizing gammaproteobacteria, and the secondary symbionts, Methylophagaceae, Nitrincolaceae and Flavobacteriaceae, whose physiology and metabolism are obscure. Little is known about if and how these symbionts interact or exchange metabolites. Here we curated metagenome-assembled genomes of Idas modiolaeformis symbionts and used genome-centered metatranscriptomics and metaproteomics to assess key symbiont functions. The Methylophagaceae symbiont is a methylotrophic autotroph, as it encoded and expressed the ribulose monophosphate and Calvin-Benson-Bassham cycle enzymes, particularly RuBisCO. The Nitrincolaceae ASP10-02a symbiont likely fuels its metabolism with nitrogen-rich macromolecules and may provide the holobiont with vitamin B12. The Urechidicola (Flavobacteriaceae) symbionts likely degrade glycans and may remove NO. Our findings indicate that these flexible associations allow for expanding the range of substrates and environmental niches, via new metabolic functions and handoffs.

Impact of High-Fiber or High-Protein Diet on the Capacity of Human Gut Microbiota To Produce Tryptophan Catabolites

This study investigated the effect of high-fiber-low-protein (HF) and high-protein-low-fiber (HP) diets on microbial catabolism of tryptophan in the proximal colon (PC) and distal colon(DC) compartments of the Simulator of the Human Intestinal Microbial Ecosystem. The microbiota in PC and DC was dominated by Bacteroidetes and Firmicutes, in which Bacteroidetes were more abundant in DC (∼60% versus 50%) and Firmicutes were more abundant in PC (∼40% versus 25%). Most of the tryptophan catabolites were determined at a higher concentration in PC samples than in DC samples, but the overall concentration of tryptophan catabolites was over 10-fold higher in DC samples than that in PC samples. Interestingly, indole-3-propionic acid and oxindole were only identified in DC samples. A two-week dietary intervention by the HF diet enriched the abundance of Firmicutes in PC, whereas the HP diet enriched the abundance of Proteobacteria. Compared to the HP diet, the HF diet favored the microbial production of indole-3-acetic acid, indole-3-lactic acid, indole-3-aldehyde, and indole-3-propionic acid in both PC and DC compartments. To conclude, these findings increase the understanding of the effect of diets on the microbial production of tryptophan catabolites in the colon.

Bacteria and Archaea Regulate Particulate Organic Matter Export in Suspended and Sinking Marine Particle Fractions

The biological carbon pump (BCP) in the Southern Ocean is driven by phytoplankton productivity and is a significant organic matter sink. However, the role of particle-attached (PA) and free-living (FL) prokaryotes (bacteria and archaea) and their diversity in influencing the efficiency of the BCP is still unclear. To investigate this, we analyzed the metagenomes linked to suspended and sinking marine particles from the Sub-Antarctic Southern Ocean Time Series (SOTS) by deploying a Marine Snow Catcher (MSC), obtaining suspended and sinking particulate material, determining organic carbon and nitrogen flux, and constructing metagenome-assembled genomes (MAGs). The suspended and sinking particle-pools were dominated by bacteria with the potential to degrade organic carbon. Bacterial communities associated with the sinking fraction had more genes related to the degradation of complex organic carbon than those in the suspended fraction. Archaea had the potential to drive nitrogen metabolism via nitrite and ammonia oxidation, altering organic nitrogen concentration. The data revealed several pathways for chemoautotrophy and the secretion of recalcitrant dissolved organic carbon (RDOC) from CO₂, with bacteria and archaea potentially sequestering particulate organic matter (POM) via the production of RDOC. These findings provide insights into the diversity and function of prokaryotes in suspended and sinking particles and their role in organic carbon/nitrogen export in the Southern Ocean. IMPORTANCE The biological carbon pump is crucial for the export of particulate organic matter in the ocean. Recent studies on marine microbes have shown the profound influence of bacteria and archaea as regulators of particulate organic matter export. Yet, despite the importance of the Southern Ocean as a carbon sink, we lack comparable insights regarding microbial contributions. This study provides the first insights regarding prokaryotic contributions to particulate organic matter export in the Southern Ocean. We reveal evidence that prokaryotic communities in suspended and sinking particle fractions harbor widespread genomic potential for mediating particulate organic matter export. The results substantially enhance our understanding of the role played by microorganisms in regulating particulate organic matter export in suspended and sinking marine fractions in the Southern Ocean.

Nanopore Sequencing Discloses Compositional Quality of Commercial Probiotic Feed Supplements

The market for the application of probiotics as a livestock health improvement supplement has increased in recent years. However, most of the available products are quality-controlled using low-resolution techniques and un-curated databases, resulting in misidentification and incorrect product labels. In this work, we deployed two workflows and compared results obtained by full-length 16S rRNA genes (16S) and metagenomic (Meta) data to investigate their reliability for the microbial composition of both liquid and solid forms of animal probiotic products using Oxford Nanopore long-read-only (without short-read). Our result revealed that 16S amplicon data permits to detect the bacterial microbiota even with the low abundance in the samples. Moreover, the 16S approach has the potential to provide species-level resolution for prokaryotes but not for assessing yeast communities. Whereas, Meta data has more power to recover of high-quality metagenome-assembled genomes that enables detailed exploration of both bacterial and yeast populations, as well as antimicrobial resistance genes, and functional genes in the population. Our findings clearly demonstrate that implementing these workflows with long-read-only monitoring could be applied to assessing the quality and safety of probiotic products for animals and evaluating the quality of probiotic products on the market. This would benefit the sustained growth of the livestock probiotic industry.

Metagenome-Assembled Genomes of Four Southern Ocean Archaea Harbor Multiple Genes Linked to Polyethylene Terephthalate and Polyhydroxybutyrate Plastic Degradation

Here, we present four archaeal metagenome-assembled genomes (MAGs) (three Thaumarchaeota MAGs and one Thermoplasmatota MAG) from a polar upwelling zone in the Southern Ocean. These archaea harbor putative genes encoding enzymes such as polyethylene terephthalate (PET) hydrolases (PETases) and polyhydroxybutyrate (PHB) depolymerases, which are associated with microbial degradation of PET and PHB plastics.

CrAssphage May Be Viable Markers of Contamination in Pristine and Contaminated River Water

Viruses are the most biologically abundant entities and may be ideal indicators of fecal pollutants in water. Anthropogenic activities have triggered drastic ecosystem changes in rivers, leading to substantial shifts in chemical and biological attributes. Here, we evaluate the viability of using the presence of crAssphage as indicators of fecal contamination in South African rivers. Shotgun analysis revealed diverse crAssphage viruses in these rivers, which are impacted by chemical and biological pollution. Overall, the diversity and relative abundances of these viruses was higher in contaminated sites compared to pristine locations. In contrast to fecal coliform counts, crAssphage sequences were detected in pristine rivers, supporting the assertion that the afore mentioned marker may be a more accurate indicator of fecal contamination. Our data demonstrate the presence of diverse putative hosts which includes members of the phyla Bacteroidota, Pseudomonadota, Verrucomicrobiota, and Bacillota. Phylogenetic analysis revealed novel subfamilies, suggesting that rivers potentially harbor distinct and uncharacterized clades of crAssphage. These data provide the first insights regarding the diversity, distribution, and functional roles of crAssphage in rivers. Taken together, the results support the potential application of crAssphage as viable markers for water quality monitoring. IMPORTANCE Rivers support substantial populations and provide important ecosystem services. Despite the application of fecal coliform tests and other markers, we lack rapid and reproducible approaches for determining fecal contamination in rivers. Waterborne viral outbreaks have been reported even after fecal indicator bacteria (FIB) were suggested to be absent or below regulated levels of coliforms. This indicates a need to develop and apply improved indicators of pollutants in aquatic ecosystems. Here, we evaluate the viability of crAssphage as indicators of fecal contamination in two South African rivers. We assess the abundance, distribution, and diversity of these viruses in sites that had been predicted pristine or contaminated by FIB analysis. We show that crAssphage are ideal and sensitive markers for fecal contamination and describe novel clades of crAss-like phages. Known crAss-like subfamilies were unrepresented in our data, suggesting that the diversity of these viruses may reflect geographic locality and dependence.

Targeting the gut-lung axis by synbiotic feeding to infants in a randomized controlled trial

BACKGROUND

Formula-fed infants are at increased risk of infections. Due to the cross-talk between the mucosal systems of the gastrointestinal and respiratory tracts, adding synbiotics (prebiotics and probiotics) to infant formula may prevent infections even at distant sites. Infants that were born full term and weaned from breast milk were randomized to prebiotic formula (fructo- and galactooligosaccharides) or the same prebiotic formula with Lactobacillus paracasei ssp. paracasei F19 (synbiotics) from 1 to 6 months of age. The objective was to examine the synbiotic effects on gut microbiota development.

RESULTS

Fecal samples collected at ages 1, 4, 6, and 12 months were analyzed using 16S rRNA gene sequencing and a combination of untargeted gas chromatography-mass spectrometry/liquid chromatography-mass spectrometry. These analyses revealed that the synbiotic group had a lower abundance of Klebsiella, a higher abundance of Bifidobacterium breve compared to the prebiotic group, and increases in the anti-microbial metabolite d-3-phenyllactic acid. We also analyzed the fecal metagenome and antibiotic resistome in the 11 infants that had been diagnosed with lower respiratory tract infection (cases) and 11 matched controls using deep metagenomic sequencing. Cases with lower respiratory tract infection had a higher abundance of Klebsiella species and antimicrobial resistance genes related to Klebsiella pneumoniae, compared to controls. The results obtained using 16S rRNA gene amplicon and metagenomic sequencing were confirmed in silico by successful recovery of the metagenome-assembled genomes of the bacteria of interest.

CONCLUSIONS

This study demonstrates the additional benefit of feeding specific synbiotics to formula-fed infants over prebiotics only. Synbiotic feeding led to the underrepresentation of Klebsiella, enrichment of bifidobacteria, and increases in microbial degradation metabolites implicated in immune signaling and in the gut-lung and gut-skin axes. Our findings support future clinical evaluation of synbiotic formula in the prevention of infections and associated antibiotic treatment as a primary outcome when breastfeeding is not feasible.

TRIAL REGISTRATION

ClinicalTrials.gov NCT01625273 . Retrospectively registered on 21 June 2012.

Limosilactobacillus reuteri administration alters the gut-brain-behavior axis in a sex-dependent manner in socially monogamous prairie voles

Research on the role of gut microbiota in behavior has grown dramatically. The probiotic L. reuteri can alter social and stress-related behaviors - yet, the underlying mechanisms remain largely unknown. Although traditional laboratory rodents provide a foundation for examining the role of L. reuteri on the gut-brain axis, they do not naturally display a wide variety of social behaviors. Using the highly-social, monogamous prairie vole (Microtus ochrogaster), we examined the effects of L. reuteri administration on behaviors, neurochemical marker expression, and gut-microbiome composition. Females, but not males, treated with live L. reuteri displayed lower levels of social affiliation compared to those treated with heat-killed L. reuteri. Overall, females displayed a lower level of anxiety-like behaviors than males. Live L. reuteri-treated females had lower expression of corticotrophin releasing factor (CRF) and CRF type-2-receptor in the nucleus accumbens, and lower vasopressin 1a-receptor in the paraventricular nucleus of the hypothalamus (PVN), but increased CRF in the PVN. There were both baseline sex differences and sex-by-treatment differences in gut microbiome composition. Live L. reuteri increased the abundance of several taxa, including Enterobacteriaceae, Lachnospiraceae NK4A136, and Treponema. Interestingly, heat-killed L. reuteri increased abundance of the beneficial taxa Bifidobacteriaceae and Blautia. There were significant correlations between changes in microbiota, brain neurochemical markers, and behaviors. Our data indicate that L. reuteri impacts gut microbiota, gut-brain axis and behaviors in a sex-specific manner in socially-monogamous prairie voles. This demonstrates the utility of the prairie vole model for further examining causal impacts of microbiome on brain and behavior.

Assessing microbiota composition in the context of aging

The gut microbiota is a complex community of different microbial species that influence many aspects of health. Consequently, shifts in the composition of gut microbiome have been proposed to exert negative effects on the host physiology, leading to the pathogenesis of various age-related disorders, including cardiovascular and neurological diseases, type 2 diabetes, obesity, non-alcoholic liver disease, and other pathological conditions. Thus, understanding how the gut microbiota influences the aging-related decline is particularly topical. Advances in next-generation sequencing techniques, together with mechanistic experiments in animal models, have provided substantial improvements in microbiome analysis. However, standardization and best practices are needed to limit experimental variation between different studies. Here, we detail a simple method for microbiota composition analysis in mouse fecal samples using 16S rRNA next-generation sequencing.

Deduplication Improves Cost-Efficiency and Yields of De Novo Assembly and Binning of Shotgun Metagenomes in Microbiome Research

In the last decade, metagenomics has greatly revolutionized the study of microbial communities. However, the presence of artificial duplicate reads raised mainly from the preparation of metagenomic DNA sequencing libraries and their impacts on metagenomic assembly and binning have never been brought to attention. Here, we explicitly investigated the effects of duplicate reads on metagenomic assemblies and binning based on analyses of five groups of representative metagenomes with distinct microbiome complexities. Our results showed that deduplication considerably increased the binning yields (by 3.5% to 80%) for most of the metagenomic data sets examined thanks to the improved contig length and coverage profiling of metagenome-assembled contigs, whereas it slightly decreased the binning yields of metagenomes with low complexity (e.g., human gut metagenomes). Specifically, 411 versus 397, 331 versus 317, 104 versus 88, and 9 versus 5 metagenome-assembled genomes (MAGs) were recovered from MEGAHIT assemblies of bioreactor sludge, surface water, lake sediment, and forest soil metagenomes, respectively. Noticeably, deduplication significantly reduced the computational costs of the metagenomic assembly, including the elapsed time (9.0% to 29.9%) and the maximum memory requirement (4.3% to 37.1%). Collectively, we recommend the removal of duplicate reads in metagenomes with high complexity before assembly and binning analyses, for example, the forest soil metagenomes examined in this study. IMPORTANCE Duplicated reads in shotgun metagenomes are usually considered technical artifacts. Their presence in metagenomes would theoretically not only introduce bias into the quantitative analysis but also result in mistakes in the coverage profile, leading to adverse effects on or even failures in metagenomic assembly and binning, as the widely used metagenome assemblers and binners all need coverage information for graph partitioning and assembly binning, respectively. However, this issue was seldom noticed, and its impacts on downstream essential bioinformatic procedures (e.g., assembly and binning) remained unclear. In this study, we comprehensively evaluated for the first time the implications of duplicate reads for the de novo assembly and binning of real metagenomic data sets by comparing the assembly qualities, binning yields, and requirements for computational resources with and without the removal of duplicate reads. It was revealed that deduplication considerably increased the binning yields of metagenomes with high complexity and significantly reduced the computational costs, including the elapsed time and the maximum memory requirement, for most of the metagenomes studied. These results provide empirical references for more cost-efficient metagenomic analyses in microbiome research.

Early human impact on lake cyanobacteria revealed by a Holocene record of sedimentary ancient DNA

Sedimentary DNA-based studies revealed the effects of human activity on lake cyanobacteria communities over the last centuries, yet we continue to lack information over longer timescales. Here, we apply high-resolution molecular analyses on sedimentary ancient DNA to reconstruct the history of cyanobacteria throughout the Holocene in a lake in north-eastern Germany. We find a substantial increase in cyanobacteria abundance coinciding with deforestation during the early Bronze Age around 4000 years ago, suggesting increased nutrient supply to the lake by local communities settling on the lakeshore. The next substantial human-driven increase in cyanobacteria abundance occurred only about a century ago due to intensified agricultural fertilisation which caused the dominance of potentially toxic taxa (e.g., Aphanizomenon). Our study provides evidence that humans began to locally impact lake ecology much earlier than previously assumed. Consequently, managing aquatic systems today requires awareness of the legacy of human influence dating back potentially several millennia.

Association between inflammation, lipopolysaccharide binding protein, and gut microbiota composition in a New Hampshire Bhutanese refugee population with a high burden of type 2 diabetes

Introduction

South Asian refugees experience a high risk of obesity and diabetes yet are often underrepresented in studies on chronic diseases and their risk factors. The gut microbiota and gut permeability, as assessed through circulating lipopolysaccharide binding protein (LBP), may underlie the link between chronic inflammation and type 2 diabetes (T2D). The composition of the gut microbiota varies according to multiple factors including demographics, migration, and dietary patterns, particularly fiber intake. However, there is no evidence on the composition of the gut microbiota and its relationship with metabolic health in refugee populations, including those migrating to the United States from Bhutan. The objective of this study was to examine glycemic status in relation to LBP, systemic inflammation fiber intake, and gut microbiota composition in Bhutanese refugee adults residing in New Hampshire (n = 50).

Methods

This cross-sectional study included a convenience sample of Bhutanese refugee adults (N = 50) in NH. Established bioinformatics pipelines for metagenomic analysis were used to determine relative genus abundance, species richness, and alpha diversity measures from shallow shotgun sequences. The relationships between inflammatory markers, gut microbiota composition, dietary fiber, and glycemic status were analyzed.

Results

We identified a substantial chronic disease burden in this study population, and observed a correlation between glycemic status, LBP, and inflammation, and a correlation between glycemic status and gut microbiome alpha diversity. Further, we identified a significant correlation between proinflammatory taxa and inflammatory cytokines. SCFA-producing taxa were found to be inversely correlated with age.

Conclusion

To date, this is the most comprehensive examination of metabolic health and the gut microbiome in a Bhutanese refugee population in NH. The findings highlight areas for future investigations of inflammation and glycemic impairment, in addition to informing potential interventions targeting this vulnerable population.

Metaphor-A workflow for streamlined assembly and binning of metagenomes

Recent advances in bioinformatics and high-throughput sequencing have enabled the large-scale recovery of genomes from metagenomes. This has the potential to bring important insights as researchers can bypass cultivation and analyze genomes sourced directly from environmental samples. There are, however, technical challenges associated with this process, most notably the complexity of computational workflows required to process metagenomic data, which include dozens of bioinformatics software tools, each with their own set of customizable parameters that affect the final output of the workflow. At the core of these workflows are the processes of assembly-combining the short-input reads into longer, contiguous fragments (contigs)-and binning, clustering these contigs into individual genome bins. The limitations of assembly and binning algorithms also pose different challenges depending on the selected strategy to execute them. Both of these processes can be done for each sample separately or by pooling together multiple samples to leverage information from a combination of samples. Here we present Metaphor, a fully automated workflow for genome-resolved metagenomics (GRM). Metaphor differs from existing GRM workflows by offering flexible approaches for the assembly and binning of the input data and by combining multiple binning algorithms with a bin refinement step to achieve high-quality genome bins. Moreover, Metaphor generates reports to evaluate the performance of the workflow. We showcase the functionality of Metaphor on different synthetic datasets and the impact of available assembly and binning strategies on the final results.

SnakeMAGs: a simple, efficient, flexible and scalable workflow to reconstruct prokaryotic genomes from metagenomes

Background: Over the last decade, we have observed in microbial ecology a transition from gene-centric to genome-centric analyses. Indeed, the advent of metagenomics combined with binning methods, single-cell genome sequencing as well as high-throughput cultivation methods have contributed to the continuing and exponential increase of available prokaryotic genomes, which in turn has favored the exploration of microbial metabolisms. In the case of metagenomics, data processing, from raw reads to genome reconstruction, involves various steps and software which can represent a major technical obstacle. Methods: To overcome this challenge, we developed SnakeMAGs, a simple workflow that can process Illumina data, from raw reads to metagenome-assembled genomes (MAGs) classification and relative abundance estimate. It integrates state-of-the-art bioinformatic tools to sequentially perform: quality control of the reads (illumina-utils, Trimmomatic), host sequence removal (optional step, using Bowtie2), assembly (MEGAHIT), binning (MetaBAT2), quality filtering of the bins (CheckM, GUNC), classification of the MAGs (GTDB-Tk) and estimate of their relative abundance (CoverM). Developed with the popular Snakemake workflow management system, it can be deployed on various architectures, from single to multicore and from workstation to computer clusters and grids. It is also flexible since users can easily change parameters and/or add new rules. Results: Using termite gut metagenomic datasets, we showed that SnakeMAGs is slower but allowed the recovery of more MAGs encompassing more diverse phyla compared to another similar workflow named ATLAS. Importantly, these additional MAGs showed no significant difference compared to the other ones in terms of completeness, contamination, genome size nor relative abundance. Conclusions: Overall, it should make the reconstruction of MAGs more accessible to microbiologists. SnakeMAGs as well as test files and an extended tutorial are available at https://github.com/Nachida08/SnakeMAGs.

SnakeMAGs: a simple, efficient, flexible and scalable workflow to reconstruct prokaryotic genomes from metagenomes

Background: Over the last decade, we have observed in microbial ecology a transition from gene-centric to genome-centric analyses. Indeed, the advent of metagenomics combined with binning methods, single-cell genome sequencing as well as high-throughput cultivation methods have contributed to the continuing and exponential increase of available prokaryotic genomes, which in turn has favored the exploration of microbial metabolisms. In the case of metagenomics, data processing, from raw reads to genome reconstruction, involves various steps and software which can represent a major technical obstacle. Methods: To overcome this challenge, we developed SnakeMAGs, a simple workflow that can process Illumina data, from raw reads to metagenome-assembled genomes (MAGs) classification and relative abundance estimate. It integrates state-of-the-art bioinformatic tools to sequentially perform: quality control of the reads (illumina-utils, Trimmomatic), host sequence removal (optional step, using Bowtie2), assembly (MEGAHIT), binning (MetaBAT2), quality filtering of the bins (CheckM), classification of the MAGs (GTDB-Tk) and estimate of their relative abundance (CoverM). Developed with the popular Snakemake workflow management system, it can be deployed on various architectures, from single to multicore and from workstation to computer clusters and grids. It is also flexible since users can easily change parameters and/or add new rules. Results: Using termite gut metagenomic datasets, we showed that SnakeMAGs is slower but allowed the recovery of more MAGs encompassing more diverse phyla compared to another similar workflow named ATLAS. Conclusions: Overall, it should make the reconstruction of MAGs more accessible to microbiologists. SnakeMAGs as well as test files and an extended tutorial are available at https://github.com/Nachida08/SnakeMAGs.

Metagenome-Assembled Genome of a Putative Chemoheterotroph from Volcanic Terrain in Hawaii

Here, we present the draft genome for a new putative species, Gaiellasilicea maunaloa, most closely related to Gaiella occulta, within the phylum Actinobacteria. This group contains Gram-negative, aerobic mesophilic species. This metagenome-assembled genome (MAG) contributes to knowledge of life in volcanic environments and may inform investigations of life beyond Earth.

Basaltic Lava Tube Hosts a Putative Novel Genus in the Family Solirubrobacteraceae

We report the draft genome sequence of a putative new genus and species, Siliceabacter maunaloa, in the family Solirubrobacteraceae. The members of this family of Actinobacteria are generally Gram positive and mesophilic. Found within a Hawaiian lava tube, this microbe illuminates the types of prokaryotes inhabiting secondary minerals in subsurface basaltic environments.

Mining the equine gut metagenome: poorly-characterized taxa associated with cardiovascular fitness in endurance athletes

Emerging evidence indicates that the gut microbiome contributes to endurance exercise performance. Still, the extent of its functional and metabolic potential remains unknown. Using elite endurance horses as a model system for exercise responsiveness, we built an integrated horse gut gene catalog comprising ~25 million unique genes and 372 metagenome-assembled genomes. This catalog represents 4179 genera spanning 95 phyla and functional capacities primed to exploit energy from dietary, microbial, and host resources. The holo-omics approach shows that gut microbiomes enriched in Lachnospiraceae taxa are negatively associated with cardiovascular capacity. Conversely, more complex and functionally diverse microbiomes are associated with higher glucose concentrations and reduced accumulation of long-chain acylcarnitines and non-esterified fatty acids in plasma, suggesting increased ß-oxidation capacity in the mitochondria. In line with this hypothesis, more fit athletes show upregulation of mitochondrial-related genes involved in energy metabolism, biogenesis, and Ca²⁺ cytosolic transport, all of which are necessary to improve aerobic work power, spare glycogen usage, and enhance cardiovascular capacity. The results identify an associative link between endurance performance and gut microbiome composition and gene function, laying the basis for nutritional interventions that could benefit horse athletes.

An integrated gene catalog of the gut microbiome in elite endurance horses is build. The holo-omics analyses identify an associative link between endurance performance and gut microbiome composition and gene function.

MAGNETO: An Automated Workflow for Genome-Resolved Metagenomics

Metagenome-assembled genomes (MAGs) represent individual genomes recovered from metagenomic data. MAGs are extremely useful to analyze uncultured microbial genomic diversity, as well as to characterize associated functional and metabolic potential in natural environments. Recent computational developments have considerably improved MAG reconstruction but also emphasized several limitations, such as the nonbinning of sequence regions with repetitions or distinct nucleotidic composition. Different assembly and binning strategies are often used; however, it still remains unclear which assembly strategy, in combination with which binning approach, offers the best performance for MAG recovery. Several workflows have been proposed in order to reconstruct MAGs, but users are usually limited to single-metagenome assembly or need to manually define sets of metagenomes to coassemble prior to genome binning. Here, we present MAGNETO, an automated workflow dedicated to MAG reconstruction, which includes a fully-automated coassembly step informed by optimal clustering of metagenomic distances, and implements complementary genome binning strategies, for improving MAG recovery. MAGNETO is implemented as a Snakemake workflow and is available at: https://gitlab.univ-nantes.fr/bird_pipeline_registry/magneto. IMPORTANCE Genome-resolved metagenomics has led to the discovery of previously untapped biodiversity within the microbial world. As the development of computational methods for the recovery of genomes from metagenomes continues, existing strategies need to be evaluated and compared to eventually lead to standardized computational workflows. In this study, we compared commonly used assembly and binning strategies and assessed their performance using both simulated and real metagenomic data sets. We propose a novel approach to automate coassembly, avoiding the requirement for a priori knowledge to combine metagenomic information. The comparison against a previous coassembly approach demonstrates a strong impact of this step on genome binning results, but also the benefits of informing coassembly for improving the quality of recovered genomes. MAGNETO integrates complementary assembly-binning strategies to optimize genome reconstruction and provides a complete reads-to-genomes workflow for the growing microbiome research community.

Exploring strain diversity of dominant human skin bacterial species using single-cell genome sequencing

To understand the role of the skin commensal bacterial community in skin health and the spread of pathogens, it is crucial to identify genetic differences in the bacterial strains corresponding to human individuals. A culture-independent genomics approach is an effective tool for obtaining massive high-quality bacterial genomes. Here we present a single-cell genome sequencing to obtain comprehensive whole-genome sequences of uncultured skin bacteria from skin swabs. We recovered 281 high-quality (HQ) and 244 medium-quality single-amplified genomes (SAGs) of multiple skin bacterial species from eight individuals, including cohabiting group. Single-cell sequencing outperformed in the genome recovery from the same skin swabs, showing 10-fold non-redundant strain genomes compared to the shotgun metagenomic sequencing and binning approach. We then focused on the abundant skin bacteria and identified intra-species diversity, especially in 47 Moraxella osloensis derived HQ SAGs, characterizing the strain-level heterogeneity at mobile genetic element profiles, including plasmids and prophages. Even between the cohabiting individual hosts, they have unique skin bacterial strains in the same species, which shows microdiversity in each host. Genetic and functional differences between skin bacterial strains are predictive of in vivo competition to adapt bacterial genome to utilize the sparse nutrients available on the skin or produce molecules that inhibit the colonization of other microbes or alter their behavior. Thus, single-cell sequencing provides a large number of genomes of higher resolution and quality than conventional metagenomic analysis and helps explore the skin commensal bacteria at the strain level, linking taxonomic and functional information.

Potential of Meta-Omics to Provide Modern Microbial Indicators for Monitoring Soil Quality and Securing Food Production

Soils are fundamental resources for agricultural production and play an essential role in food security. They represent the keystone of the food value chain because they harbor a large fraction of biodiversity—the backbone of the regulation of ecosystem services and “soil health” maintenance. In the face of the numerous causes of soil degradation such as unsustainable soil management practices, pollution, waste disposal, or the increasing number of extreme weather events, it has become clear that (i) preserving the soil biodiversity is key to food security, and (ii) biodiversity-based solutions for environmental monitoring have to be developed. Within the soil biodiversity reservoir, microbial diversity including Archaea, Bacteria, Fungi and protists is essential for ecosystem functioning and resilience. Microbial communities are also sensitive to various environmental drivers and to management practices; as a result, they are ideal candidates for monitoring soil quality assessment. The emergence of meta-omics approaches based on recent advances in high-throughput sequencing and bioinformatics has remarkably improved our ability to characterize microbial diversity and its potential functions. This revolution has substantially filled the knowledge gap about soil microbial diversity regulation and ecology, but also provided new and robust indicators of agricultural soil quality. We reviewed how meta-omics approaches replaced traditional methods and allowed developing modern microbial indicators of the soil biological quality. Each meta-omics approach is described in its general principles, methodologies, specificities, strengths and drawbacks, and illustrated with concrete applications for soil monitoring. The development of metabarcoding approaches in the last 20 years has led to a collection of microbial indicators that are now operational and available for the farming sector. Our review shows that despite the recent huge advances, some meta-omics approaches (e.g., metatranscriptomics or meta-proteomics) still need developments to be operational for environmental bio-monitoring. As regards prospects, we outline the importance of building up repositories of soil quality indicators. These are essential for objective and robust diagnosis, to help actors and stakeholders improve soil management, with a view to or to contribute to combining the food and environmental quality of next-generation farming systems in the context of the agroecological transition.

Recovery of High Quality Metagenome-Assembled Genomes From Full-Scale Activated Sludge Microbial Communities in a Tropical Climate Using Longitudinal Metagenome Sampling

The analysis of metagenome data based on the recovery of draft genomes (so called metagenome-assembled genomes, or MAG) has assumed an increasingly central role in microbiome research in recent years. Microbial communities underpinning the operation of wastewater treatment plants are particularly challenging targets for MAG analysis due to their high ecological complexity, and remain important, albeit understudied, microbial communities that play ssa key role in mediating interactions between human and natural ecosystems. Here we consider strategies for recovery of MAG sequence from time series metagenome surveys of full-scale activated sludge microbial communities. We generate MAG catalogs from this set of data using several different strategies, including the use of multiple individual sample assemblies, two variations on multi-sample co-assembly and a recently published MAG recovery workflow using deep learning. We obtain a total of just under 9,100 draft genomes, which collapse to around 3,100 non-redundant genomic clusters. We examine the strengths and weaknesses of these approaches in relation to MAG yield and quality, showing that co-assembly may offer advantages over single-sample assembly in the case of metagenome data obtained from closely sampled longitudinal study designs. Around 1,000 MAGs were candidates for being considered high quality, based on single-copy marker gene occurrence statistics, however only 58 MAG formally meet the MIMAG criteria for being high quality draft genomes. These findings carry broader broader implications for performing genome-resolved metagenomics on highly complex communities, the design and implementation of genome recoverability strategies, MAG decontamination and the search for better binning methodology.

The Distinguishing Bacterial Features From Active and Remission Stages of Ulcerative Colitis Revealed by Paired Fecal Metagenomes

Ulcerative colitis (UC) is a serious chronic intestinal inflammatory disease, with an increased incidence in recent years. The intestinal microbiota plays a key role in the pathogenesis of UC. However, there is no unified conclusion on how the intestinal microbiota changes. Most studies focus on the change between UC patients and healthy individuals, rather than the active and remission stage of the same patient. To minimize the influences of genetic differences, environmental and dietary factors, we studied the intestinal microbiota of paired fecal samples from 42 UC patients at the active and remission stages. We identified 175 species of microbes from 11 phyla and found no difference of the alpha and beta diversities between the active and remission stages. Paired t-test analysis revealed differential microbiota at levels of the phyla, class, order, family, genus, and species, including 13 species with differential abundance. For example, CAG-269 sp001916005, Eubacterium F sp003491505, Lachnospira sp000436475, et al. were downregulated in the remission, while the species of Parabacteroides distasonis, Prevotellamassilia sp900540885, CAG-495 sp001917125, et al. were upregulated in the remission. The 13 species can effectively distinguish the active and remission stages. Functional analysis showed that the sporulation and biosynthesis were downregulated, and the hydrogen peroxide catabolic process was upregulated in remission of UC. Our study suggests that the 13 species together may serve as a biomarker panel contributing to identify the active and remission stages of UC, which provides a valuable reference for the treatment of UC patients by FMT or other therapeutic methods.