Increased circulating butyrate and ursodeoxycholate during probiotic intervention in humans with type 2 diabetes

Background

An increasing body of evidence implicates the resident gut microbiota as playing a critical role in type 2 diabetes (T2D) pathogenesis. We previously reported significant improvement in postprandial glucose control in human participants with T2D following 12-week administration of a 5-strain novel probiotic formulation (‘WBF-011’) in a double-blind, randomized, placebo controlled setting (NCT03893422). While the clinical endpoints were encouraging, additional exploratory measurements were needed in order to link the motivating mechanistic hypothesis - increased short-chain fatty acids - with markers of disease.

Results

Here we report targeted and untargeted metabolomic measurements on fasting plasma (n = 104) collected at baseline and end of intervention. Butyrate and ursodeoxycholate increased among participants randomized to WBF-011, along with compelling trends between butyrate and glycated haemoglobin (HbA1c). In vitro monoculture experiments demonstrated that the formulation’s C. butyricum strain efficiently synthesizes ursodeoxycholate from the primary bile acid chenodeoxycholate during butyrogenic growth. Untargeted metabolomics also revealed coordinated decreases in intermediates of fatty acid oxidation and bilirubin, potential secondary signatures for metabolic improvement. Finally, improvement in HbA1c was limited almost entirely to participants not using sulfonylurea drugs. We show that these drugs can inhibit growth of formulation strains in vitro.

Conclusion

To our knowledge, this is the first description of an increase in circulating butyrate or ursodeoxycholate following a probiotic intervention in humans with T2D, adding support for the possibility of a targeted microbiome-based approach to assist in the management of T2D. The efficient synthesis of UDCA by C. butyricum is also likely of interest to investigators of its use as a probiotic in other disease settings. The potential for inhibitory interaction between sulfonylurea drugs and gut microbiota should be considered carefully in the design of future studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02415-8.

Comparative genomics reveals electron transfer and syntrophic mechanisms differentiating methanotrophic and methanogenic archaea

The anaerobic oxidation of methane coupled to sulfate reduction is a microbially mediated process requiring a syntrophic partnership between anaerobic methanotrophic (ANME) archaea and sulfate-reducing bacteria (SRB). Based on genome taxonomy, ANME lineages are polyphyletic within the phylum Halobacterota, none of which have been isolated in pure culture. Here, we reconstruct 28 ANME genomes from environmental metagenomes and flow sorted syntrophic consortia. Together with a reanalysis of previously published datasets, these genomes enable a comparative analysis of all marine ANME clades. We review the genomic features that separate ANME from their methanogenic relatives and identify what differentiates ANME clades. Large multiheme cytochromes and bioenergetic complexes predicted to be involved in novel electron bifurcation reactions are well distributed and conserved in the ANME archaea, while significant variations in the anabolic C1 pathways exists between clades. Our analysis raises the possibility that methylotrophic methanogenesis may have evolved from a methanotrophic ancestor.

A comparative genomics study of anaerobic methanotrophic (ANME) archaea reveals the genetic "parts list" associated with the repeated evolutionary transition between methanogenic and methanotrophic metabolism in the archaeal domain of life.

Physiological, genomic, and sulfur isotopic characterization of methanol metabolism by Desulfovibrio carbinolicus

Methanol is often considered as a non-competitive substrate for methanogenic archaea, but an increasing number of sulfate-reducing microorganisms (SRMs) have been reported to be capable of respiring with methanol as an electron donor. A better understanding of the fate of methanol in natural or artificial anaerobic systems thus requires knowledge of the methanol dissimilation by SRMs. In this study, we describe the growth kinetics and sulfur isotope effects of Desulfovibrio carbinolicus, a methanol-oxidizing sulfate-reducing deltaproteobacterium, together with its genome sequence and annotation. D. carbinolicus can grow with a series of alcohols from methanol to butanol. Compared to longer-chain alcohols, however, specific growth and respiration rates decrease by several fold with methanol as an electron donor. Larger sulfur isotope fractionation accompanies slowed growth kinetics, indicating low chemical potential at terminal reductive steps of respiration. In a medium containing both ethanol and methanol, D. carbinolicus does not consume methanol even after the cessation of growth on ethanol. Among the two known methanol dissimilatory systems, the genome of D. carbinolicus contains the genes coding for alcohol dehydrogenase but lacks enzymes analogous to methanol methyltransferase. We analyzed the genomes of 52 additional species of sulfate-reducing bacteria that have been tested for methanol oxidation. There is no apparent relationship between phylogeny and methanol metabolizing capacity, but most gram-negative methanol oxidizers grow poorly, and none carry homologs for methyltransferase (mtaB). Although the amount of available data is limited, it is notable that more than half of the known gram-positive methanol oxidizers have both enzymatic systems, showing enhanced growth relative to the SRMs containing only alcohol dehydrogenase genes. Thus, physiological, genomic, and sulfur isotopic results suggest that D. carbinolicus and close relatives have the ability to metabolize methanol but likely play a limited role in methanol degradation in most natural environments.

Comparative Genomics and Proteomic Analysis of Assimilatory Sulfate Reduction Pathways in Anaerobic Methanotrophic Archaea

Sulfate is the predominant electron acceptor for anaerobic oxidation of methane (AOM) in marine sediments. This process is carried out by a syntrophic consortium of anaerobic methanotrophic archaea (ANME) and sulfate reducing bacteria (SRB) through an energy conservation mechanism that is still poorly understood. It was previously hypothesized that ANME alone could couple methane oxidation to dissimilatory sulfate reduction, but a genetic and biochemical basis for this proposal has not been identified. Using comparative genomic and phylogenetic analyses, we found the genetic capacity in ANME and related methanogenic archaea for sulfate reduction, including sulfate adenylyltransferase, APS kinase, APS/PAPS reductase and two different sulfite reductases. Based on characterized homologs and the lack of associated energy conserving complexes, the sulfate reduction pathways in ANME are likely used for assimilation but not dissimilation of sulfate. Environmental metaproteomic analysis confirmed the expression of 6 proteins in the sulfate assimilation pathway of ANME. The highest expressed proteins related to sulfate assimilation were two sulfite reductases, namely assimilatory-type low-molecular-weight sulfite reductase (alSir) and a divergent group of coenzyme F₄₂₀-dependent sulfite reductase (Group II Fsr). In methane seep sediment microcosm experiments, however, sulfite and zero-valent sulfur amendments were inhibitory to ANME-2a/2c while growth in their syntrophic SRB partner was not observed. Combined with our genomic and metaproteomic results, the passage of sulfur species by ANME as metabolic intermediates for their SRB partners is unlikely. Instead, our findings point to a possible niche for ANME to assimilate inorganic sulfur compounds more oxidized than sulfide in anoxic marine environments.

An expanded genomic representation of the phylum cyanobacteria

Molecular surveys of aphotic habitats have indicated the presence of major uncultured lineages phylogenetically classified as members of the Cyanobacteria. One of these lineages has recently been proposed as a nonphotosynthetic sister phylum to the Cyanobacteria, the Melainabacteria, based on recovery of population genomes from human gut and groundwater samples. Here, we expand the phylogenomic representation of the Melainabacteria through sequencing of six diverse population genomes from gut and bioreactor samples supporting the inference that this lineage is nonphotosynthetic, but not the assertion that they are strictly fermentative. We propose that the Melainabacteria is a class within the phylogenetically defined Cyanobacteria based on robust monophyly and shared ancestral traits with photosynthetic representatives. Our findings are consistent with theories that photosynthesis occurred late in the Cyanobacteria and involved extensive lateral gene transfer and extends the recognized functionality of members of this phylum.

Proteomic Stable Isotope Probing Reveals Biosynthesis Dynamics of Slow Growing Methane Based Microbial Communities

Marine methane seep habitats represent an important control on the global flux of methane. Nucleotide-based meta-omics studies outline community-wide metabolic potential, but expression patterns of environmentally relevant proteins are poorly characterized. Proteomic stable isotope probing (proteomic SIP) provides additional information by characterizing phylogenetically specific, functionally relevant activity in mixed microbial communities, offering enhanced detection through system-wide product integration. Here we applied proteomic SIP to ¹⁵NH4+ and CH₄ amended seep sediment microcosms in an attempt to track protein synthesis of slow-growing, low-energy microbial systems. Across all samples, 3495 unique proteins were identified, 11% of which were ¹⁵N-labeled. Consistent with the dominant anaerobic oxidation of methane (AOM) activity commonly observed in anoxic seep sediments, proteins associated with sulfate reduction and reverse methanogenesis—including the ANME-2 associated methylenetetrahydromethanopterin reductase (Mer)—were all observed to be actively synthesized (¹⁵N-enriched). Conversely, proteins affiliated with putative aerobic sulfur-oxidizing epsilon- and gammaproteobacteria showed a marked decrease over time in our anoxic sediment incubations. The abundance and phylogenetic range of ¹⁵N-enriched methyl-coenzyme M reductase (Mcr) orthologs, many of which exhibited novel post-translational modifications, suggests that seep sediments provide niches for multiple organisms performing analogous metabolisms. In addition, 26 proteins of unknown function were consistently detected and actively expressed under conditions supporting AOM, suggesting that they play important roles in methane seep ecosystems. Stable isotope probing in environmental proteomics experiments provides a mechanism to determine protein durability and evaluate lineage-specific responses in complex microbial communities placed under environmentally relevant conditions. Our work here demonstrates the active synthesis of a metabolically specific minority of enzymes, revealing the surprising longevity of most proteins over the course of an extended incubation experiment in an established, slow-growing, methane-impacted environmental system.

Genomic Reconstruction of an Uncultured Hydrothermal Vent Gammaproteobacterial Methanotroph (Family Methylothermaceae) Indicates Multiple Adaptations to Oxygen Limitation

Hydrothermal vents are an important contributor to marine biogeochemistry, producing large volumes of reduced fluids, gasses, and metals and housing unique, productive microbial and animal communities fueled by chemosynthesis. Methane is a common constituent of hydrothermal vent fluid and is frequently consumed at vent sites by methanotrophic bacteria that serve to control escape of this greenhouse gas into the atmosphere. Despite their ecological and geochemical importance, little is known about the ecophysiology of uncultured hydrothermal vent-associated methanotrophic bacteria. Using metagenomic binning techniques, we recovered and analyzed a near-complete genome from a novel gammaproteobacterial methanotroph (B42) associated with a white smoker chimney in the Southern Lau basin. B42 was the dominant methanotroph in the community, at ∼80x coverage, with only four others detected in the metagenome, all on low coverage contigs (7x–12x). Phylogenetic placement of B42 showed it is a member of the Methylothermaceae, a family currently represented by only one sequenced genome. Metabolic inferences based on the presence of known pathways in the genome showed that B42 possesses a branched respiratory chain with A- and B-family heme copper oxidases, cytochrome bd oxidase and a partial denitrification pathway. These genes could allow B42 to respire over a wide range of oxygen concentrations within the highly dynamic vent environment. Phylogenies of the denitrification genes revealed they are the result of separate horizontal gene transfer from other Proteobacteria and suggest that denitrification is a selective advantage in conditions where extremely low oxygen concentrations require all oxygen to be used for methane activation.

Metagenomic and metaproteomic analyses of Accumulibacter phosphatis-enriched floccular and granular biofilm

Biofilms are ubiquitous in nature, forming diverse adherent microbial communities that perform a plethora of functions. Here we operated two laboratory-scale sequencing batch reactors enriched with Candidatus Accumulibacter phosphatis (Accumulibacter) performing enhanced biological phosphorus removal. Reactors formed two distinct biofilms, one floccular biofilm, consisting of small, loose, microbial aggregates, and one granular biofilm, forming larger, dense, spherical aggregates. Using metagenomic and metaproteomic methods, we investigated the proteomic differences between these two biofilm communities, identifying a total of 2022 unique proteins. To understand biofilm differences, we compared protein abundances that were statistically enriched in both biofilm states. Floccular biofilms were enriched with pathogenic secretion systems suggesting a highly competitive microbial community. Comparatively, granular biofilms revealed a high-stress environment with evidence of nutrient starvation, phage predation pressure, and increased extracellular polymeric substance and cell lysis. Granular biofilms were enriched in outer membrane transport proteins to scavenge the extracellular milieu for amino acids and other metabolites, likely released through cell lysis, to supplement metabolic pathways. This study provides the first detailed proteomic comparison between Accumulibacter-enriched floccular and granular biofilm communities, proposes a conceptual model for the granule biofilm, and offers novel insights into granule biofilm formation and stability.

The khmer software package: enabling efficient nucleotide sequence analysis

The khmer package is a freely available software library for working efficiently with fixed length DNA words, or k-mers. khmer provides implementations of a probabilistic k-mer counting data structure, a compressible De Bruijn graph representation, De Bruijn graph partitioning, and digital normalization. khmer is implemented in C++ and Python, and is freely available under the BSD license at  https://github.com/dib-lab/khmer/.

CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes

Large-scale recovery of genomes from isolates, single cells, and metagenomic data has been made possible by advances in computational methods and substantial reductions in sequencing costs. Although this increasing breadth of draft genomes is providing key information regarding the evolutionary and functional diversity of microbial life, it has become impractical to finish all available reference genomes. Making robust biological inferences from draft genomes requires accurate estimates of their completeness and contamination. Current methods for assessing genome quality are ad hoc and generally make use of a limited number of "marker" genes conserved across all bacterial or archaeal genomes. Here we introduce CheckM, an automated method for assessing the quality of a genome using a broader set of marker genes specific to the position of a genome within a reference genome tree and information about the collocation of these genes. We demonstrate the effectiveness of CheckM using synthetic data and a wide range of isolate-, single-cell-, and metagenome-derived genomes. CheckM is shown to provide accurate estimates of genome completeness and contamination and to outperform existing approaches. Using CheckM, we identify a diverse range of errors currently impacting publicly available isolate genomes and demonstrate that genomes obtained from single cells and metagenomic data vary substantially in quality. In order to facilitate the use of draft genomes, we propose an objective measure of genome quality that can be used to select genomes suitable for specific gene- and genome-centric analyses of microbial communities.

CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes

Large-scale recovery of genomes from isolates, single cells, and metagenomic data has been made possible by advances in computational methods and substantial reductions in sequencing costs. Although this increasing breadth of draft genomes is providing key information regarding the evolutionary and functional diversity of microbial life, it has become impractical to finish all available reference genomes. Making robust biological inferences from draft genomes requires accurate estimates of their completeness and contamination. Current methods for assessing genome quality are ad hoc and generally make use of a limited number of “marker” genes conserved across all bacterial or archaeal genomes. Here we introduce CheckM, an automated method for assessing the quality of a genome using a broader set of marker genes specific to the position of a genome within a reference genome tree and information about the collocation of these genes. We demonstrate the effectiveness of CheckM using synthetic data and a wide range of isolate-, single-cell-, and metagenome-derived genomes. CheckM is shown to provide accurate estimates of genome completeness and contamination and to outperform existing approaches. Using CheckM, we identify a diverse range of errors currently impacting publicly available isolate genomes and demonstrate that genomes obtained from single cells and metagenomic data vary substantially in quality. In order to facilitate the use of draft genomes, we propose an objective measure of genome quality that can be used to select genomes suitable for specific gene- and genome-centric analyses of microbial communities.

Expanding our view of genomic diversity in Candidatus Accumulibacter clades

Enhanced biological phosphorus removal (EBPR) is an important industrial wastewater treatment process mediated by polyphosphate-accumulating organisms (PAOs). Members of the genus Candidatus Accumulibacter are one of the most extensively studied PAO as they are commonly enriched in lab-scale EBPR reactors. Members of different Accumulibacter clades are often enriched through changes in reactor process conditions; however, the two currently sequenced Accumulibacter genomes show extensive metabolic similarity. Here, we expand our understanding of Accumulibacter genomic diversity through recovery of eight population genomes using deep metagenomics, including seven from phylogenetic clades with no previously sequenced representative. Comparative genomic analysis revealed a core of shared genes involved primarily in carbon and phosphorus metabolism; however, each Accumulibacter genome also encoded a substantial number of unique genes (> 700 genes). A major difference between the Accumulibacter clades was the type of nitrate reductase encoded and the capacity to perform subsequent steps in denitrification. The Accumulibacter clade IIF genomes also contained acetaldehyde dehydrogenase that may allow ethanol to be used as carbon source. These differences in metabolism between Accumulibacter genomes provide a molecular basis for niche differentiation observed in lab-scale reactors and may offer new opportunities for process optimization.

Crass: identification and reconstruction of CRISPR from unassembled metagenomic data

Clustered regularly interspaced short palindromic repeats (CRISPR) constitute a bacterial and archaeal adaptive immune system that protect against bacteriophage (phage). Analysis of CRISPR loci reveals the history of phage infections and provides a direct link between phage and their hosts. All current tools for CRISPR identification have been developed to analyse completed genomes and are not well suited to the analysis of metagenomic data sets, where CRISPR loci are difficult to assemble owing to their repetitive structure and population heterogeneity. Here, we introduce a new algorithm, Crass, which is designed to identify and reconstruct CRISPR loci from raw metagenomic data without the need for assembly or prior knowledge of CRISPR in the data set. CRISPR in assembled data are often fragmented across many contigs/scaffolds and do not fully represent the population heterogeneity of CRISPR loci. Crass identified substantially more CRISPR in metagenomes previously analysed using assembly-based approaches. Using Crass, we were able to detect CRISPR that contained spacers with sequence homology to phage in the system, which would not have been identified using other approaches. The increased sensitivity, specificity and speed of Crass will facilitate comprehensive analysis of CRISPRs in metagenomic data sets, increasing our understanding of phage-host interactions and co-evolution within microbial communities.

In-solution fluorescence in situ hybridization and fluorescence-activated cell sorting for single cell and population genome recovery

Over the past decade, technological advances in whole genome amplification, microfluidics, flow sorting, and high-throughput sequencing have led to the development of single-cell genomics. Single-cell genomic approaches are typically applied to anonymous microbial cells with only morphology providing clues to their identity. However, targeted separation of microorganisms based on phylogenetic markers, such as the 16S rRNA gene, is beginning to emerge in the single-cell genomics field. Here, we describe an in-solution fluorescence in situ hybridization (FISH) protocol which can be combined with fluorescence-activated cell sorting (FACS) for separation of single cells or populations of interest from environmental samples. Sequencing of DNA obtained from sorted cells can be used for the recovery of draft quality genomes, and when performed in parallel with deep metagenomics, can be used to validate and further scaffold metagenomic assemblies. We illustrate in this chapter the feasibility of this FISH-FACS approach by describing the targeted recovery of a novel anaerobic methanotrophic archaeon.