A cryptic plasmid is among the most numerous genetic elements in the human gut.

Plasmids are extrachromosomal genetic elements that often encode fitness-enhancing features. However, many bacteria carry "cryptic" plasmids that do not confer clear beneficial functions. We identified one such cryptic plasmid, pBI143, which is ubiquitous across industrialized gut microbiomes and is 14 times as numerous as crAssphage, currently established as the most abundant extrachromosomal genetic element in the human gut. The majority of mutations in pBI143 accumulate in specific positions across thousands of metagenomes, indicating strong purifying selection. pBI143 is monoclonal in most individuals, likely due to the priority effect of the version first acquired, often from one's mother. pBI143 can transfer between Bacteroidales, and although it does not appear to impact bacterial host fitness in vivo, it can transiently acquire additional genetic content. We identified important practical applications of pBI143, including its use in identifying human fecal contamination and its potential as an alternative approach to track human colonic inflammatory states.

Gene Expression Changes and Community Turnover Differentially Shape the Global Ocean Metatranscriptome.

Ocean microbial communities strongly influence the biogeochemistry, food webs, and climate of our planet. Despite recent advances in understanding their taxonomic and genomic compositions, little is known about how their transcriptomes vary globally. Here, we present a dataset of 187 metatranscriptomes and 370 metagenomes from 126 globally distributed sampling stations and establish a resource of 47 million genes to study community-level transcriptomes across depth layers from pole-to-pole. We examine gene expression changes and community turnover as the underlying mechanisms shaping community transcriptomes along these axes of environmental variation and show how their individual contributions differ for multiple biogeochemically relevant processes. Furthermore, we find the relative contribution of gene expression changes to be significantly lower in polar than in non-polar waters and hypothesize that in polar regions, alterations in community activity in response to ocean warming will be driven more strongly by changes in organismal composition than by gene regulatory mechanisms. VIDEO ABSTRACT.

Functional and evolutionary significance of unknown genes from uncultivated taxa.

Many of the Earth's microbes remain uncultured and understudied, limiting our understanding of the functional and evolutionary aspects of their genetic material, which remain largely overlooked in most metagenomic studies1. Here we analysed 149,842 environmental genomes from multiple habitats2-6 and compiled a curated catalogue of 404,085 functionally and evolutionarily significant novel (FESNov) gene families exclusive to uncultivated prokaryotic taxa. All FESNov families span multiple species, exhibit strong signals of purifying selection and qualify as new orthologous groups, thus nearly tripling the number of bacterial and archaeal gene families described to date. The FESNov catalogue is enriched in clade-specific traits, including 1,034 novel families that can distinguish entire uncultivated phyla, classes and orders, probably representing synapomorphies that facilitated their evolutionary divergence. Using genomic context analysis and structural alignments we predicted functional associations for 32.4% of FESNov families, including 4,349 high-confidence associations with important biological processes. These predictions provide a valuable hypothesis-driven framework that we used for experimental validatation of a new gene family involved in cell motility and a novel set of antimicrobial peptides. We also demonstrate that the relative abundance profiles of novel families can discriminate between environments and clinical conditions, leading to the discovery of potentially new biomarkers associated with colorectal cancer. We expect this work to enhance future metagenomics studies and expand our knowledge of the genetic repertory of uncultivated organisms.

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.

SPIRE: a Searchable, Planetary-scale mIcrobiome REsource.

Meta'omic data on microbial diversity and function accrue exponentially in public repositories, but derived information is often siloed according to data type, study or sampled microbial environment. Here we present SPIRE, a Searchable Planetary-scale mIcrobiome REsource that integrates various consistently processed metagenome-derived microbial data modalities across habitats, geography and phylogeny. SPIRE encompasses 99 146 metagenomic samples from 739 studies covering a wide array of microbial environments and augmented with manually-curated contextual data. Across a total metagenomic assembly of 16 Tbp, SPIRE comprises 35 billion predicted protein sequences and 1.16 million newly constructed metagenome-assembled genomes (MAGs) of medium or high quality. Beyond mapping to the high-quality genome reference provided by proGenomes3 (http://progenomes.embl.de), these novel MAGs form 92 134 novel species-level clusters, the majority of which are unclassified at species level using current tools. SPIRE enables taxonomic profiling of these species clusters via an updated, custom mOTUs database (https://motu-tool.org/) and includes several layers of functional annotation, as well as crosslinks to several (micro-)biological databases. The resource is accessible, searchable and browsable via http://spire.embl.de.

mBARq: a versatile and user-friendly framework for the analysis of DNA barcodes from transposon insertion libraries, knockout mutants, and isogenic strain populations.

MOTIVATION: DNA barcoding has become a powerful tool for assessing the fitness of strains in a variety of studies, including random transposon mutagenesis screens, attenuation of site-directed mutants, and population dynamics of isogenic strain pools. However, the statistical analysis, visualization, and contextualization of the data resulting from such experiments can be complex and require bioinformatic skills. RESULTS: Here, we developed mBARq, a user-friendly tool designed to simplify these steps for diverse experimental setups. The tool is seamlessly integrated with an intuitive web app for interactive data exploration via the STRING and KEGG databases to accelerate scientific discovery. AVAILABILITY AND IMPLEMENTATION: The tool is implemented in Python. The source code is freely available (https://github.com/MicrobiologyETHZ/mbarq) and the web app can be accessed at: https://microbiomics.io/tools/mbarq-app.

mTAGs: taxonomic profiling using degenerate consensus reference sequences of ribosomal RNA genes.

Profiling the taxonomic composition of microbial communities commonly involves the classification of ribosomal RNA gene fragments. As a trade-off to maintain high classification accuracy, existing tools are typically limited to the genus level. Here, we present mTAGs, a taxonomic profiling tool that implements the alignment of metagenomic sequencing reads to degenerate consensus reference sequences of small subunit ribosomal RNA genes. It uses DNA fragments, that is, paired-end sequencing reads, as count units and provides relative abundance profiles at multiple taxonomic ranks, including operational taxonomic units based on a 97% sequence identity cutoff. At the genus rank, mTAGs outperformed other tools across several metrics, such as the F1 score by >11% across data from different environments, and achieved competitive (F1 score) or better results (Bray-Curtis dissimilarity) at the sub-genus level. AVAILABILITY AND IMPLEMENTATION: The software tool mTAGs is implemented in Python. The source code and binaries are freely available (https://github.com/SushiLab/mTAGs). The data underlying this article are available in Zenodo, at https://doi.org/10.5281/zenodo.4352762. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Multi-omics data integration reveals novel drug targets in hepatocellular carcinoma.

BACKGROUND: Genetic aberrations in hepatocellular carcinoma (HCC) are well known, but the functional consequences of such aberrations remain poorly understood. RESULTS: Here, we explored the effect of defined genetic changes on the transcriptome, proteome and phosphoproteome in twelve tumors from an mTOR-driven hepatocellular carcinoma mouse model. Using Network-based Integration of multi-omiCS data (NetICS), we detected 74 'mediators' that relay via molecular interactions the effects of genetic and miRNA expression changes. The detected mediators account for the effects of oncogenic mTOR signaling on the transcriptome, proteome and phosphoproteome. We confirmed the dysregulation of the mediators YAP1, GRB2, SIRT1, HDAC4 and LIS1 in human HCC. CONCLUSIONS: This study suggests that targeting pathways such as YAP1 or GRB2 signaling and pathways regulating global histone acetylation could be beneficial in treating HCC with hyperactive mTOR signaling.

Impact of manipulation of glycerol/diol dehydratase activity on intestinal microbiota ecology and metabolism.

Glycerol/diol dehydratases (GDH) are enzymes that catalyse the production of propionate from 1,2-propanediol, and acrolein from glycerol. Acrolein reacts with dietary carcinogenic heterocyclic amines (HCA), reducing HCA mutagenicity, but is itself also an antimicrobial agent and toxicant. Gut microbial GDH activity has been suggested as an endogenous acrolein source; however, there is limited information on the potential of the intestinal microbiota to have GDH activity, and what impact it can have on the intestinal ecosystem and host health. We hypothesized that GDH activity of gut microbiota is determined by the abundance and distribution of GDH-active taxa and can be enhanced by supplementation of the GDH active Anaerobutyricum hallii, and tested this hypothesis combining quantitative profiling of gdh, model batch fermentations, microbiota manipulation, and kinetic modelling of acrolein formation. Our results suggest that GDH activity is a common trait of intestinal microbiota shared by a few taxa, which was dependent on overall gdh abundance. Anaerobutyricum hallii was identified as a key taxon in GDH metabolism, and its supplementation increased the rate of GDH activity and acrolein release, which enhanced the transformation of HCA and reduced fermentation activity. The findings of this first systematic study on acrolein release by intestinal microbiota indicate that dietary and microbial modulation might impact GDH activity, which may influence host health.

Bioinformatics for precision oncology.

Molecular profiling of tumor biopsies plays an increasingly important role not only in cancer research, but also in the clinical management of cancer patients. Multi-omics approaches hold the promise of improving diagnostics, prognostics and personalized treatment. To deliver on this promise of precision oncology, appropriate bioinformatics methods for managing, integrating and analyzing large and complex data are necessary. Here, we discuss the specific requirements of bioinformatics methods and software that arise in the setting of clinical oncology, owing to a stricter regulatory environment and the need for rapid, highly reproducible and robust procedures. We describe the workflow of a molecular tumor board and the specific bioinformatics support that it requires, from the primary analysis of raw molecular profiling data to the automatic generation of a clinical report and its delivery to decision-making clinical oncologists. Such workflows have to various degrees been implemented in many clinical trials, as well as in molecular tumor boards at specialized cancer centers and university hospitals worldwide. We review these and more recent efforts to include other high-dimensional multi-omics patient profiles into the tumor board, as well as the state of clinical decision support software to translate molecular findings into treatment recommendations.

NGS-pipe: a flexible, easily extendable and highly configurable framework for NGS analysis.

Motivation: Next-generation sequencing is now an established method in genomics, and massive amounts of sequencing data are being generated on a regular basis. Analysis of the sequencing data is typically performed by lab-specific in-house solutions, but the agreement of results from different facilities is often small. General standards for quality control, reproducibility and documentation are missing. Results: We developed NGS-pipe, a flexible, transparent and easy-to-use framework for the design of pipelines to analyze whole-exome, whole-genome and transcriptome sequencing data. NGS-pipe facilitates the harmonization of genomic data analysis by supporting quality control, documentation, reproducibility, parallelization and easy adaptation to other NGS experiments. Availability and implementation: https://github.com/cbg-ethz/NGS-pipe. Contact: niko.beerenwinkel@bsse.ethz.ch.

Gut microbial beta-glucuronidase and glycerol/diol dehydratase activity contribute to dietary heterocyclic amine biotransformation.

BACKGROUND: Consuming red and processed meat has been associated with an increased risk of colorectal cancer (CRC), which is partly attributed to exposure to carcinogens such as heterocyclic amines (HCA) formed during cooking and preservation processes. The interaction of gut microbes and HCA can result in altered bioactivities and it has been shown previously that human gut microbiota can transform mutagenic HCA to a glycerol conjugate with reduced mutagenic potential. However, the major form of HCA in the colon are glucuronides (HCA-G) and it is not known whether these metabolites, via stepwise microbial hydrolysis and acrolein conjugation, are viable precursors for glycerol conjugated metabolites. We hypothesized that such a process could be concurrently catalyzed by bacterial beta-glucuronidase (B-GUS) and glycerol/diol dehydratase (GDH) activity. We therefore investigated how the HCA-G PhIP-N2-ß-D-glucuronide (PhIP-G), a representative liver metabolite of PhIP (2-Amino-1-methyl-6-phenylimidazo [4,5-b] pyridine), which is the most abundant carcinogenic HCA in well-cooked meat, is transformed by enzymatic activity of human gut microbial representatives of the phyla Firmicutes, Bacteroidetes, and Proteobacteria. RESULTS: We employed a combination of growth and enzymatic assays, and a bioanalysis approach combined with metagenomics. B-GUS of Faecalibacterium prausnitzii converted PhIP-G to PhIP and GDH of Flavonifractor plautii, Blautia obeum, Eubacterium hallii, and Lactobacillus reuteri converted PhIP to PhIP-M1 in the presence of glycerol. In addition, B-GUS- and GDH-positive bacteria cooperatively converted PhIP-G to PhIP-M1. A screen of genes encoding B-GUS and GDH was performed for fecal microbiome data from healthy individuals (n = 103) and from CRC patients (n = 53), which revealed a decrease in abundance of taxa with confirmed GDH and HCA transformation activity in CRC patients. CONCLUSIONS: This study for the first time demonstrates that gut microbes mediate the stepwise transformation of PhIP-G to PhIP-M1 via the intermediate production of PhIP. Findings from this study suggest that targeted manipulation with gut microbes bearing specific functions, or dietary glycerol supplementation might modify gut microbial activity to reduce HCA-induced CRC risk.

Intestinal Dysbiosis, Barrier Dysfunction, and Bacterial Translocation Account for CKD-Related Systemic Inflammation.

CKD associates with systemic inflammation, but the underlying cause is unknown. Here, we investigated the involvement of intestinal microbiota. We report that collagen type 4 α3-deficient mice with Alport syndrome-related progressive CKD displayed systemic inflammation, including increased plasma levels of pentraxin-2 and activated antigen-presenting cells, CD4 and CD8 T cells, and Th17- or IFNγ-producing T cells in the spleen as well as regulatory T cell suppression. CKD-related systemic inflammation in these mice associated with intestinal dysbiosis of proteobacterial blooms, translocation of living bacteria across the intestinal barrier into the liver, and increased serum levels of bacterial endotoxin. Uremia did not affect secretory IgA release into the ileum lumen or mucosal leukocyte subsets. To test for causation between dysbiosis and systemic inflammation in CKD, we eradicated facultative anaerobic microbiota with antibiotics. This eradication prevented bacterial translocation, significantly reduced serum endotoxin levels, and fully reversed all markers of systemic inflammation to the level of nonuremic controls. Therefore, we conclude that uremia associates with intestinal dysbiosis, intestinal barrier dysfunction, and bacterial translocation, which trigger the state of persistent systemic inflammation in CKD. Uremic dysbiosis and intestinal barrier dysfunction may be novel therapeutic targets for intervention to suppress CKD-related systemic inflammation and its consequences.

MEGAN Community Edition - Interactive Exploration and Analysis of Large-Scale Microbiome Sequencing Data.

There is increasing interest in employing shotgun sequencing, rather than amplicon sequencing, to analyze microbiome samples. Typical projects may involve hundreds of samples and billions of sequencing reads. The comparison of such samples against a protein reference database generates billions of alignments and the analysis of such data is computationally challenging. To address this, we have substantially rewritten and extended our widely-used microbiome analysis tool MEGAN so as to facilitate the interactive analysis of the taxonomic and functional content of very large microbiome datasets. Other new features include a functional classifier called InterPro2GO, gene-centric read assembly, principal coordinate analysis of taxonomy and function, and support for metadata. The new program is called MEGAN Community Edition (CE) and is open source. By integrating MEGAN CE with our high-throughput DNA-to-protein alignment tool DIAMOND and by providing a new program MeganServer that allows access to metagenome analysis files hosted on a server, we provide a straightforward, yet powerful and complete pipeline for the analysis of metagenome shotgun sequences. We illustrate how to perform a full-scale computational analysis of a metagenomic sequencing project, involving 12 samples and 800 million reads, in less than three days on a single server. All source code is available here: https://github.com/danielhuson/megan-ce.

Endotoxicity of lipopolysaccharide as a determinant of T-cell-mediated colitis induction in mice.

BACKGROUND & AIMS: The intestinal microbiota is an important determinant of the mucosal response. In patients with inflammatory bowel diseases, the mucosal immune system has inappropriate interactions with the intestinal microbiota. We investigated how the composition of the intestinal microbiota affects its endotoxicity and development of colitis in mice. METHODS: Germ-free C57BL/6J-Rag(1tm1Mom) (Rag1(-/-)) mice were colonized with 2 different types of complex intestinal microbiota. Colitis was induced in Rag1(-/-) mice by transfer of CD4(+)CD62L(+) T cells from C57BL/6J mice. Colonic tissues were collected and used for histologic analysis and cell isolation. Activation of lamina propria dendritic cells and T cells was analyzed by flow cytometry. RESULTS: After transfer of CD4(+)CD62L(+) T cells, mice with intestinal Endo(lo) microbiota (a low proportion of Enterobacteriaceae, high proportion of Bacteroidetes, and low endotoxicity) maintained mucosal immune homeostasis, and mice with highly endotoxic Endo(hi) microbiota (a high proportion of Enterobacteriaceae and low proportion of Bacteroidetes) developed colitis. To determine whether the effects of Endo(hi) microbiota were related to the higher endotoxic activity of lipopolysaccharide (LPS), we compared LPS from Enterobacteriaceae with that of Bacteroidetes. Administration of Escherichia coli JM83 (wild-type LPS) to the mice exacerbated colitis, and Escherichia coli JM83 + htrBPG (mutated LPS, with lower endotoxicity, similar to that of Bacteroidetes) prevented development of colitis after transfer of the T cells to mice. CONCLUSIONS: The endotoxicity of LPS produced by the intestinal microbiota is a determinant of whether mice develop colitis after transfer of CD4(+)CD62L(+) T cells. This finding might aid the design of novel biologics or probiotics to treat inflammatory bowel disease.

Integrative analysis of environmental sequences using MEGAN4.

A major challenge in the analysis of environmental sequences is data integration. The question is how to analyze different types of data in a unified approach, addressing both the taxonomic and functional aspects. To facilitate such analyses, we have substantially extended MEGAN, a widely used taxonomic analysis program. The new program, MEGAN4, provides an integrated approach to the taxonomic and functional analysis of metagenomic, metatranscriptomic, metaproteomic, and rRNA data. While taxonomic analysis is performed based on the NCBI taxonomy, functional analysis is performed using the SEED classification of subsystems and functional roles or the KEGG classification of pathways and enzymes. A number of examples illustrate how such analyses can be performed, and show that one can also import and compare classification results obtained using others' tools. MEGAN4 is freely available for academic purposes, and installers for all three major operating systems can be downloaded from www-ab.informatik.uni-tuebingen.de/software/megan.

Niche availability and competitive loss by facilitation control proliferation of bacterial strains intended for soil microbiome interventions.

Microbiome engineering - the targeted manipulation of microbial communities - is considered a promising strategy to restore ecosystems, but experimental support and mechanistic understanding are required. Here, we show that bacterial inoculants for soil microbiome engineering may fail to establish because they inadvertently facilitate growth of native resident microbiomes. By generating soil microcosms in presence or absence of standardized soil resident communities, we show how different nutrient availabilities limit outgrowth of focal bacterial inoculants (three Pseudomonads), and how this might be improved by adding an artificial, inoculant-selective nutrient niche. Through random paired interaction assays in agarose micro-beads, we demonstrate that, in addition to direct competition, inoculants lose competitiveness by facilitating growth of resident soil bacteria. Metatranscriptomics experiments with toluene as selective nutrient niche for the inoculant Pseudomonas veronii indicate that this facilitation is due to loss and uptake of excreted metabolites by resident taxa. Generation of selective nutrient niches for inoculants may help to favor their proliferation for the duration of their intended action while limiting their competitive loss.

Phenology and ecological role of aerobic anoxygenic phototrophs in freshwaters.

BACKGROUND: Aerobic anoxygenic phototrophic (AAP) bacteria are heterotrophic bacteria that supply their metabolism with light energy harvested by bacteriochlorophyll-a-containing reaction centers. Despite their substantial contribution to bacterial biomass, microbial food webs, and carbon cycle, their phenology in freshwater lakes remains unknown. Hence, we investigated seasonal variations of AAP abundance and community composition biweekly across 3 years in a temperate, meso-oligotrophic freshwater lake. RESULTS: AAP bacteria displayed a clear seasonal trend with a spring maximum following the bloom of phytoplankton and a secondary maximum in autumn. As the AAP bacteria represent a highly diverse assemblage of species, we followed their seasonal succession using the amplicon sequencing of the pufM marker gene. To enhance the accuracy of the taxonomic assignment, we developed new pufM primers that generate longer amplicons and compiled the currently largest database of pufM genes, comprising 3633 reference sequences spanning all phyla known to contain AAP species. With this novel resource, we demonstrated that the majority of the species appeared during specific phases of the seasonal cycle, with less than 2% of AAP species detected during the whole year. AAP community presented an indigenous freshwater nature characterized by high resilience and heterogenic adaptations to varying conditions of the freshwater environment. CONCLUSIONS: Our findings highlight the substantial contribution of AAP bacteria to the carbon flow and ecological dynamics of lakes and unveil a recurrent and dynamic seasonal succession of the AAP community. By integrating this information with the indicator of primary production (Chlorophyll-a) and existing ecological models, we show that AAP bacteria play a pivotal role in the recycling of dissolved organic matter released during spring phytoplankton bloom. We suggest a potential role of AAP bacteria within the context of the PEG model and their consideration in further ecological models.

Diversity of the Pacific Ocean coral reef microbiome.

Coral reefs are among the most diverse ecosystems on Earth. They support high biodiversity of multicellular organisms that strongly rely on associated microorganisms for health and nutrition. However, the extent of the coral reef microbiome diversity and its distribution at the oceanic basin-scale remains to be explored. Here, we systematically sampled 3 coral morphotypes, 2 fish species, and planktonic communities in 99 reefs from 32 islands across the Pacific Ocean, to assess reef microbiome composition and biogeography. We show a very large richness of reef microorganisms compared to other environments, which extrapolated to all fishes and corals of the Pacific, approximates the current estimated total prokaryotic diversity for the entire Earth. Microbial communities vary among and within the 3 animal biomes (coral, fish, plankton), and geographically. For corals, the cross-ocean patterns of diversity are different from those known for other multicellular organisms. Within each coral morphotype, community composition is always determined by geographic distance first, both at the island and across ocean scale, and then by environment. Our unprecedented sampling effort of coral reef microbiomes, as part of the Tara Pacific expedition, provides new insight into the global microbial diversity, the factors driving their distribution, and the biocomplexity of reef ecosystems.

Ecology of Endozoicomonadaceae in three coral genera across the Pacific Ocean.

Health and resilience of the coral holobiont depend on diverse bacterial communities often dominated by key marine symbionts of the Endozoicomonadaceae family. The factors controlling their distribution and their functional diversity remain, however, poorly known. Here, we study the ecology of Endozoicomonadaceae at an ocean basin-scale by sampling specimens from three coral genera (Pocillopora, Porites, Millepora) on 99 reefs from 32 islands across the Pacific Ocean. The analysis of 2447 metabarcoding and 270 metagenomic samples reveals that each coral genus harbored a distinct new species of Endozoicomonadaceae. These species are composed of nine lineages that have distinct biogeographic patterns. The most common one, found in Pocillopora, appears to be a globally distributed symbiont with distinct metabolic capabilities, including the synthesis of amino acids and vitamins not produced by the host. The other lineages are structured partly by the host genetic lineage in Pocillopora and mainly by the geographic location in Porites. Millepora is more rarely associated to Endozoicomonadaceae. Our results show that different coral genera exhibit distinct strategies of host-Endozoicomonadaceae associations that are defined at the bacteria lineage level.