Systematic analysis of drug combinations against Gram-positive bacteria.

Drug combinations can expand options for antibacterial therapies but have not been systematically tested in Gram-positive species. We profiled ~8,000 combinations of 65 antibacterial drugs against the model species Bacillus subtilis and two prominent pathogens, Staphylococcus aureus and Streptococcus pneumoniae. Thereby, we recapitulated previously known drug interactions, but also identified ten times more novel interactions in the pathogen S. aureus, including 150 synergies. We showed that two synergies were equally effective against multidrug-resistant S. aureus clinical isolates in vitro and in vivo. Interactions were largely species-specific and synergies were distinct from those of Gram-negative species, owing to cell surface and drug uptake differences. We also tested 2,728 combinations of 44 commonly prescribed non-antibiotic drugs with 62 drugs with antibacterial activity against S. aureus and identified numerous antagonisms that might compromise the efficacy of antimicrobial therapies. We identified even more synergies and showed that the anti-aggregant ticagrelor synergized with cationic antibiotics by modifying the surface charge of S. aureus. All data can be browsed in an interactive interface ( https://apps.embl.de/combact/ ).

Location and condition based reconstruction of colon cancer microbiome from human RNA sequencing data.

BACKGROUND: The association between microbes and cancer has been reported repeatedly; however, it is not clear if molecular tumour properties are connected to specific microbial colonisation patterns. This is due mainly to the current technical and analytical strategy limitations to characterise tumour-associated bacteria. METHODS: Here, we propose an approach to detect bacterial signals in human RNA sequencing data and associate them with the clinical and molecular properties of the tumours. The method was tested on public datasets from The Cancer Genome Atlas, and its accuracy was assessed on a new cohort of colorectal cancer patients. RESULTS: Our analysis shows that intratumoural microbiome composition is correlated with survival, anatomic location, microsatellite instability, consensus molecular subtype and immune cell infiltration in colon tumours. In particular, we find Faecalibacterium prausnitzii, Coprococcus comes, Bacteroides spp., Fusobacterium spp. and Clostridium spp. to be strongly associated with tumour properties. CONCLUSIONS: We implemented an approach to concurrently analyse clinical and molecular properties of the tumour as well as the composition of the associated microbiome. Our results may improve patient stratification and pave the path for mechanistic studies on microbiota-tumour crosstalk.

Remodelling of the intestinal ecosystem during caloric restriction and fasting.

Benefits of fasting and caloric restriction on host metabolic health are well established, but less is known about the effects on the gut microbiome and how this impacts renewal of the intestinal mucosa. What has been repeatedly shown during fasting, however, is that bacteria utilising host-derived substrates proliferate at the expense of those relying on dietary substrates. Considering the increased recognition of the gut microbiome's role in maintaining host (metabolic) health, disentangling host-microbe interactions and establishing their physiological relevance in the context of fasting and caloric restriction is crucial. Such insights could aid in moving away from associations of gut bacterial signatures with metabolic diseases consistently reported in observational studies to potentially establishing causality. Therefore, this review aims to summarise what is currently known or still controversial about the interplay between fasting and caloric restriction, the gut microbiome and intestinal tissue physiology.

PyHMMER: a Python library binding to HMMER for efficient sequence analysis.

SUMMARY: PyHMMER provides Python integration of the popular profile Hidden Markov Model software HMMER via Cython bindings. This allows the annotation of protein sequences with profile HMMs and building new ones directly with Python. PyHMMER increases flexibility of use, allowing creating queries directly from Python code, launching searches, and obtaining results without I/O, or accessing previously unavailable statistics like uncorrected P-values. A new parallelization model greatly improves performance when running multithreaded searches, while producing the exact same results as HMMER. AVAILABILITY AND IMPLEMENTATION: PyHMMER supports all modern Python versions (Python 3.6+) and similar platforms as HMMER (x86 or PowerPC UNIX systems). Pre-compiled packages are released via PyPI (https://pypi.org/project/pyhmmer/) and Bioconda (https://anaconda.org/bioconda/pyhmmer). The PyHMMER source code is available under the terms of the open-source MIT licence and hosted on GitHub (https://github.com/althonos/pyhmmer); its documentation is available on ReadTheDocs (https://pyhmmer.readthedocs.io).

proGenomes3: approaching one million accurately and consistently annotated high-quality prokaryotic genomes.

The interpretation of genomic, transcriptomic and other microbial 'omics data is highly dependent on the availability of well-annotated genomes. As the number of publicly available microbial genomes continues to increase exponentially, the need for quality control and consistent annotation is becoming critical. We present proGenomes3, a database of 907 388 high-quality genomes containing 4 billion genes that passed stringent criteria and have been consistently annotated using multiple functional and taxonomic databases including mobile genetic elements and biosynthetic gene clusters. proGenomes3 encompasses 41 171 species-level clusters, defined based on universal single copy marker genes, for which pan-genomes and contextual habitat annotations are provided. The database is available at http://progenomes.embl.de/.

Cultivation-independent genomes greatly expand taxonomic-profiling capabilities of mOTUs across various environments.

BACKGROUND: Taxonomic profiling is a fundamental task in microbiome research that aims to detect and quantify the relative abundance of microorganisms in biological samples. Available methods using shotgun metagenomic data generally depend on the deposition of sequenced and taxonomically annotated genomes, usually from cultures of isolated strains, in reference databases (reference genomes). However, the majority of microorganisms have not been cultured yet. Thus, a substantial fraction of microbial community members remains unaccounted for during taxonomic profiling, particularly in samples from underexplored environments. To address this issue, we developed the mOTU profiler, a tool that enables reference genome-independent species-level profiling of metagenomes. As such, it supports the identification and quantification of both "known" and "unknown" species based on a set of select marker genes. RESULTS: We present mOTUs3, a command line tool that enables the profiling of metagenomes for >33,000 species-level operational taxonomic units. To achieve this, we leveraged the reconstruction of >600,000 draft genomes, most of which are metagenome-assembled genomes (MAGs), from diverse microbiomes, including soil, freshwater systems, and the gastrointestinal tract of ruminants and other animals, which we found to be underrepresented by reference genomes. Overall, two thirds of all species-level taxa lacked a reference genome. The cumulative relative abundance of these newly included taxa was low in well-studied microbiomes, such as the human body sites (6-11%). By contrast, they accounted for substantial proportions (ocean, freshwater, soil: 43-63%) or even the majority (pig, fish, cattle: 60-80%) of the relative abundance across diverse non-human-associated microbiomes. Using community-developed benchmarks and datasets, we found mOTUs3 to be more accurate than other methods and to be more congruent with 16S rRNA gene-based methods for taxonomic profiling. Furthermore, we demonstrate that mOTUs3 increases the resolution of well-known microbial groups into species-level taxa and helps identify new differentially abundant taxa in comparative metagenomic studies. CONCLUSIONS: We developed mOTUs3 to enable accurate species-level profiling of metagenomes. Compared to other methods, it provides a more comprehensive view of prokaryotic community diversity, in particular for currently underexplored microbiomes. To facilitate comparative analyses by the research community, it is released with >11,000 precomputed profiles for publicly available metagenomes and is freely available at: https://github.com/motu-tool/mOTUs . Video Abstract.

Depression and fatigue in active IBD from a microbiome perspective-a Bayesian approach to faecal metagenomics.

BACKGROUND: Extraintestinal symptoms are common in inflammatory bowel diseases (IBD) and include depression and fatigue. These are highly prevalent especially in active disease, potentially due to inflammation-mediated changes in the microbiota-gut-brain axis. The aim of this study was to investigate the associations between structural and functional microbiota characteristics and severity of fatigue and depressive symptoms in patients with active IBD. METHODS: We included clinical data of 62 prospectively enrolled patients with IBD in an active disease state. Patients supplied stool samples and completed the questionnaires regarding depression and fatigue symptoms. Based on taxonomic and functional metagenomic profiles of faecal gut microbiota, we used Bayesian statistics to investigate the associative networks and triangle motifs between bacterial genera, functional modules and symptom severity of self-reported fatigue and depression. RESULTS: Associations with moderate to strong evidence were found for 3 genera (Odoribacter, Anaerotruncus and Alistipes) and 3 functional modules (pectin, glycosaminoglycan and central carbohydrate metabolism) with regard to depression and for 4 genera (Intestinimonas, Anaerotruncus, Eubacterium and Clostridiales g.i.s) and 2 functional modules implicating amino acid and central carbohydrate metabolism with regard to fatigue. CONCLUSIONS: This study provides the first evidence of association triplets between microbiota composition, function and extraintestinal symptoms in active IBD. Depression and fatigue were associated with lower abundances of short-chain fatty acid producers and distinct pathways implicating glycan, carbohydrate and amino acid metabolism. Our results suggest that microbiota-directed therapeutic approaches may reduce fatigue and depression in IBD and should be investigated in future research.

Bacterial retrons encode phage-defending tripartite toxin-antitoxin systems.

Retrons are prokaryotic genetic retroelements encoding a reverse transcriptase that produces multi-copy single-stranded DNA1 (msDNA). Despite decades of research on the biosynthesis of msDNA2, the function and physiological roles of retrons have remained unknown. Here we show that Retron-Sen2 of Salmonella enterica serovar Typhimurium encodes an accessory toxin protein, STM14_4640, which we renamed as RcaT. RcaT is neutralized by the reverse transcriptase-msDNA antitoxin complex, and becomes active upon perturbation of msDNA biosynthesis. The reverse transcriptase is required for binding to RcaT, and the msDNA is required for the antitoxin activity. The highly prevalent RcaT-containing retron family constitutes a new type of tripartite DNA-containing toxin-antitoxin system. To understand the physiological roles of such toxin-antitoxin systems, we developed toxin activation-inhibition conjugation (TAC-TIC), a high-throughput reverse genetics approach that identifies the molecular triggers and blockers of toxin-antitoxin systems. By applying TAC-TIC to Retron-Sen2, we identified multiple trigger and blocker proteins of phage origin. We demonstrate that phage-related triggers directly modify the msDNA, thereby activating RcaT and inhibiting bacterial growth. By contrast, prophage proteins circumvent retrons by directly blocking RcaT. Consistently, retron toxin-antitoxin systems act as abortive infection anti-phage defence systems, in line with recent reports3,4. Thus, RcaT retrons are tripartite DNA-regulated toxin-antitoxin systems, which use the reverse transcriptase-msDNA complex both as an antitoxin and as a sensor of phage protein activities.

Biosynthetic potential of the global ocean microbiome.

Natural microbial communities are phylogenetically and metabolically diverse. In addition to underexplored organismal groups1, this diversity encompasses a rich discovery potential for ecologically and biotechnologically relevant enzymes and biochemical compounds2,3. However, studying this diversity to identify genomic pathways for the synthesis of such compounds4 and assigning them to their respective hosts remains challenging. The biosynthetic potential of microorganisms in the open ocean remains largely uncharted owing to limitations in the analysis of genome-resolved data at the global scale. Here we investigated the diversity and novelty of biosynthetic gene clusters in the ocean by integrating around 10,000 microbial genomes from cultivated and single cells with more than 25,000 newly reconstructed draft genomes from more than 1,000 seawater samples. These efforts revealed approximately 40,000 putative mostly new biosynthetic gene clusters, several of which were found in previously unsuspected phylogenetic groups. Among these groups, we identified a lineage rich in biosynthetic gene clusters ('Candidatus Eudoremicrobiaceae') that belongs to an uncultivated bacterial phylum and includes some of the most biosynthetically diverse microorganisms in this environment. From these, we characterized the phospeptin and pythonamide pathways, revealing cases of unusual bioactive compound structure and enzymology, respectively. Together, this research demonstrates how microbiomics-driven strategies can enable the investigation of previously undescribed enzymes and natural products in underexplored microbial groups and environments.

Microbiota-dependent activation of the myeloid calcineurin-NFAT pathway inhibits B7H3- and B7H4-dependent anti-tumor immunity in colorectal cancer.

Bacterial sensing by intestinal tumor cells contributes to tumor growth through cell-intrinsic activation of the calcineurin-NFAT axis, but the role of this pathway in other intestinal cells remains unclear. Here, we found that myeloid-specific deletion of calcineurin in mice activated protective CD8+ T cell responses and inhibited colorectal cancer (CRC) growth. Microbial sensing by myeloid cells promoted calcineurin- and NFAT-dependent interleukin 6 (IL-6) release, expression of the co-inhibitory molecules B7H3 and B7H4 by tumor cells, and inhibition of CD8+ T cell-dependent anti-tumor immunity. Accordingly, targeting members of this pathway activated protective CD8+ T cell responses and inhibited primary and metastatic CRC growth. B7H3 and B7H4 were expressed by the majority of human primary CRCs and metastases, which was associated with low numbers of tumor-infiltrating CD8+ T cells and poor survival. Therefore, a microbiota-, calcineurin-, and B7H3/B7H4-dependent pathway controls anti-tumor immunity, revealing additional targets for immune checkpoint inhibition in microsatellite-stable CRC.

Calorie restriction improves metabolic state independently of gut microbiome composition: a randomized dietary intervention trial.

BACKGROUND: The gut microbiota has been suggested to play a significant role in the development of overweight and obesity. However, the effects of calorie restriction on gut microbiota of overweight and obese adults, especially over longer durations, are largely unexplored. METHODS: Here, we longitudinally analyzed the effects of intermittent calorie restriction (ICR) operationalized as the 5:2 diet versus continuous calorie restriction (CCR) on fecal microbiota of 147 overweight or obese adults in a 50-week parallel-arm randomized controlled trial, the HELENA Trial. The primary outcome of the trial was the differential effects of ICR versus CCR on gene expression in subcutaneous adipose tissue. Changes in the gut microbiome, which are the focus of this publication, were defined as exploratory endpoint of the trial. The trial comprised a 12-week intervention period, a 12-week maintenance period, and a final follow-up period of 26 weeks. RESULTS: Both diets resulted in ~5% weight loss. However, except for Lactobacillales being enriched after ICR, post-intervention microbiome composition did not significantly differ between groups. Overall weight loss was associated with significant metabolic improvements, but not with changes in the gut microbiome. Nonetheless, the abundance of the Dorea genus at baseline was moderately predictive of subsequent weight loss (AUROC of 0.74 for distinguishing the highest versus lowest weight loss quartiles). Despite the lack of consistent intervention effects on microbiome composition, significant study group-independent co-variation between gut bacterial families and metabolic biomarkers, anthropometric measures, and dietary composition was detectable. Our analysis in particular revealed associations between insulin sensitivity (HOMA-IR) and Akkermansiaceae, Christensenellaceae, and Tanerellaceae. It also suggests the possibility of a beneficial modulation of the latter two intestinal taxa by a diet high in vegetables and fiber, and low in processed meat. CONCLUSIONS: Overall, our results suggest that the gut microbiome remains stable and highly individual-specific under dietary calorie restriction. TRIAL REGISTRATION: The trial, including the present microbiome component, was prospectively registered at ClinicalTrials.gov NCT02449148 on May 20, 2015.

A faecal microbiota signature with high specificity for pancreatic cancer.

BACKGROUND: Recent evidence suggests a role for the microbiome in pancreatic ductal adenocarcinoma (PDAC) aetiology and progression. OBJECTIVE: To explore the faecal and salivary microbiota as potential diagnostic biomarkers. METHODS: We applied shotgun metagenomic and 16S rRNA amplicon sequencing to samples from a Spanish case-control study (n=136), including 57 cases, 50 controls, and 29 patients with chronic pancreatitis in the discovery phase, and from a German case-control study (n=76), in the validation phase. RESULTS: Faecal metagenomic classifiers performed much better than saliva-based classifiers and identified patients with PDAC with an accuracy of up to 0.84 area under the receiver operating characteristic curve (AUROC) based on a set of 27 microbial species, with consistent accuracy across early and late disease stages. Performance further improved to up to 0.94 AUROC when we combined our microbiome-based predictions with serum levels of carbohydrate antigen (CA) 19-9, the only current non-invasive, Food and Drug Administration approved, low specificity PDAC diagnostic biomarker. Furthermore, a microbiota-based classification model confined to PDAC-enriched species was highly disease-specific when validated against 25 publicly available metagenomic study populations for various health conditions (n=5792). Both microbiome-based models had a high prediction accuracy on a German validation population (n=76). Several faecal PDAC marker species were detectable in pancreatic tumour and non-tumour tissue using 16S rRNA sequencing and fluorescence in situ hybridisation. CONCLUSION: Taken together, our results indicate that non-invasive, robust and specific faecal microbiota-based screening for the early detection of PDAC is feasible.

Identifying temporal and spatial patterns of variation from multimodal data using MEFISTO.

Factor analysis is a widely used method for dimensionality reduction in genome biology, with applications from personalized health to single-cell biology. Existing factor analysis models assume independence of the observed samples, an assumption that fails in spatio-temporal profiling studies. Here we present MEFISTO, a flexible and versatile toolbox for modeling high-dimensional data when spatial or temporal dependencies between the samples are known. MEFISTO maintains the established benefits of factor analysis for multimodal data, but enables the performance of spatio-temporally informed dimensionality reduction, interpolation, and separation of smooth from non-smooth patterns of variation. Moreover, MEFISTO can integrate multiple related datasets by simultaneously identifying and aligning the underlying patterns of variation in a data-driven manner. To illustrate MEFISTO, we apply the model to different datasets with spatial or temporal resolution, including an evolutionary atlas of organ development, a longitudinal microbiome study, a single-cell multi-omics atlas of mouse gastrulation and spatially resolved transcriptomics.

Unravelling the collateral damage of antibiotics on gut bacteria.

Antibiotics are used to fight pathogens but also target commensal bacteria, disturbing the composition of gut microbiota and causing dysbiosis and disease1. Despite this well-known collateral damage, the activity spectrum of different antibiotic classes on gut bacteria remains poorly characterized. Here we characterize further 144 antibiotics from a previous screen of more than 1,000 drugs on 38 representative human gut microbiome species2. Antibiotic classes exhibited distinct inhibition spectra, including generation dependence for quinolones and phylogeny independence for ß-lactams. Macrolides and tetracyclines, both prototypic bacteriostatic protein synthesis inhibitors, inhibited nearly all commensals tested but also killed several species. Killed bacteria were more readily eliminated from in vitro communities than those inhibited. This species-specific killing activity challenges the long-standing distinction between bactericidal and bacteriostatic antibiotic classes and provides a possible explanation for the strong effect of macrolides on animal3-5 and human6,7 gut microbiomes. To mitigate this collateral damage of macrolides and tetracyclines, we screened for drugs that specifically antagonized the antibiotic activity against abundant Bacteroides species but not against relevant pathogens. Such antidotes selectively protected Bacteroides species from erythromycin treatment in human-stool-derived communities and gnotobiotic mice. These findings illluminate the activity spectra of antibiotics in commensal bacteria and suggest strategies to circumvent their adverse effects on the gut microbiota.

mOTUs: Profiling Taxonomic Composition, Transcriptional Activity and Strain Populations of Microbial Communities.

The mOTU profiler, or mOTUs for short, is a software tool that enables the profiling of microbial communities in terms of their taxonomic composition, relative abundance of metabolically active members, and diversity of strain populations. To this end, it maintains a database of single-copy phylogenetic marker gene sequences, which are used as a reference to which short read metagenomic and metatranscriptomic reads are mapped for the identification and quantification of microbial taxa. Here, we describe the most common use cases of the mOTU profiler in two basic protocols. Additional supporting protocols provide information on its installation and in-depth guidance on adjusting its settings for increasing or decreasing the stringency with which taxa are detected and quantified, as well as for customizing the output file format. Guidelines for understanding the profiling results are provided, along with additional information on unique features, methodological details, and the development history of the tool. © 2021 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Metagenomic and metatranscriptomic mOTU profiling Basic Protocol 2: Metagenomic SNV profiling Support Protocol 1: Installing mOTUs Support Protocol 2: Profiling pipeline-step by step Support Protocol 3: The mOTUs profiling routine using advanced parameters Support Protocol 4: Metagenomic SNV calling: advanced parameters.

Commensal Clostridiales strains mediate effective anti-cancer immune response against solid tumors.

Despite overall success, T cell checkpoint inhibitors for cancer treatment are still only efficient in a minority of patients. Recently, intestinal microbiota was found to critically modulate anti-cancer immunity and therapy response. Here, we identify Clostridiales members of the gut microbiota associated with a lower tumor burden in mouse models of colorectal cancer (CRC). Interestingly, these commensal species are also significantly reduced in CRC patients compared with healthy controls. Oral application of a mix of four Clostridiales strains (CC4) in mice prevented and even successfully treated CRC as stand-alone therapy. This effect depended on intratumoral infiltration and activation of CD8+ T cells. Single application of Roseburia intestinalis or Anaerostipes caccae was even more effective than CC4. In a direct comparison, the CC4 mix supplementation outperformed anti-PD-1 therapy in mouse models of CRC and melanoma. Our findings provide a strong preclinical foundation for exploring gut bacteria as novel stand-alone therapy against solid tumors.

Microbiome meta-analysis and cross-disease comparison enabled by the SIAMCAT machine learning toolbox.

The human microbiome is increasingly mined for diagnostic and therapeutic biomarkers using machine learning (ML). However, metagenomics-specific software is scarce, and overoptimistic evaluation and limited cross-study generalization are prevailing issues. To address these, we developed SIAMCAT, a versatile R toolbox for ML-based comparative metagenomics. We demonstrate its capabilities in a meta-analysis of fecal metagenomic studies (10,803 samples). When naively transferred across studies, ML models lost accuracy and disease specificity, which could however be resolved by a novel training set augmentation strategy. This reveals some biomarkers to be disease-specific, with others shared across multiple conditions. SIAMCAT is freely available from siamcat.embl.de .

Statistical and Machine Learning Techniques in Human Microbiome Studies: Contemporary Challenges and Solutions.

The human microbiome has emerged as a central research topic in human biology and biomedicine. Current microbiome studies generate high-throughput omics data across different body sites, populations, and life stages. Many of the challenges in microbiome research are similar to other high-throughput studies, the quantitative analyses need to address the heterogeneity of data, specific statistical properties, and the remarkable variation in microbiome composition across individuals and body sites. This has led to a broad spectrum of statistical and machine learning challenges that range from study design, data processing, and standardization to analysis, modeling, cross-study comparison, prediction, data science ecosystems, and reproducible reporting. Nevertheless, although many statistics and machine learning approaches and tools have been developed, new techniques are needed to deal with emerging applications and the vast heterogeneity of microbiome data. We review and discuss emerging applications of statistical and machine learning techniques in human microbiome studies and introduce the COST Action CA18131 "ML4Microbiome" that brings together microbiome researchers and machine learning experts to address current challenges such as standardization of analysis pipelines for reproducibility of data analysis results, benchmarking, improvement, or development of existing and new tools and ontologies.

Metabolic cooperation and spatiotemporal niche partitioning in a kefir microbial community.

Microbial communities often undergo intricate compositional changes yet also maintain stable coexistence of diverse species. The mechanisms underlying long-term coexistence remain unclear as system-wide studies have been largely limited to engineered communities, ex situ adapted cultures or synthetic assemblies. Here, we show how kefir, a natural milk-fermenting community of prokaryotes (predominantly lactic and acetic acid bacteria) and yeasts (family Saccharomycetaceae), realizes stable coexistence through spatiotemporal orchestration of species and metabolite dynamics. During milk fermentation, kefir grains (a polysaccharide matrix synthesized by kefir microorganisms) grow in mass but remain unchanged in composition. In contrast, the milk is colonized in a sequential manner in which early members open the niche for the followers by making available metabolites such as amino acids and lactate. Through metabolomics, transcriptomics and large-scale mapping of inter-species interactions, we show how microorganisms poorly suited for milk survive in-and even dominate-the community, through metabolic cooperation and uneven partitioning between grain and milk. Overall, our findings reveal how inter-species interactions partitioned in space and time lead to stable coexistence.

The gut microbiome in solid organ transplantation.

Despite ground-breaking advances in allogeneic transplantation, allograft rejection and immunosuppressant-specific complications remain a major challenge in transplant medicine. Growing evidence suggests the human gut microbiome as a potential contributor to transplant outcome and patient health. After breakthrough findings in haematopoietic stem cell transplantation (HSCT), the relevance of the microbiome in solid organ transplantation (SOT) is becoming increasingly clear. Here, we review the role of the microbiome in SOT focusing on its significance for transplant-associated complications such as allograft rejection and infections, and highlight its potential impact on immunosuppressive treatment. Moreover, we shed light on the emerging role of the microbiome as a diagnostic biomarker and therapeutic target and discuss current microbial intervention strategies. In addition, this review includes some practical considerations in designing clinical microbiome trials and offers some advice for the interpretation of the resulting data. Further investigation of the gut microbiome harbours countless clinical application possibilities and holds great promise of having a lasting impact on transplant medicine.