Deep metagenomics examines the oral microbiome during dental caries, revealing novel taxa and co-occurrences with host molecules.

Dental caries, the most common chronic infectious disease worldwide, has a complex etiology involving the interplay of microbial and host factors that are not completely understood. In this study, the oral microbiome and 38 host cytokines and chemokines were analyzed across 23 children with caries and 24 children with healthy dentition. De novo assembly of metagenomic sequencing obtained 527 metagenome-assembled genomes (MAGs), representing 150 bacterial species. Forty-two of these species had no genomes in public repositories, thereby representing novel taxa. These new genomes greatly expanded the known pangenomes of many oral clades, including the enigmatic Saccharibacteria clades G3 and G6, which had distinct functional repertoires compared to other oral Saccharibacteria. Saccharibacteria are understood to be obligate epibionts, which are dependent on host bacteria. These data suggest that the various Saccharibacteria clades may rely on their hosts for highly distinct metabolic requirements, which would have significant evolutionary and ecological implications. Across the study group, Rothia, Neisseria, and Haemophilus spp. were associated with good dental health, whereas Prevotella spp., Streptococcus mutans, and Human herpesvirus 4 (Epstein-Barr virus [EBV]) were more prevalent in children with caries. Finally, 10 of the host immunological markers were significantly elevated in the caries group, and co-occurrence analysis provided an atlas of potential relationships between microbes and host immunological molecules. Overall, this study illustrated the oral microbiome at an unprecedented resolution and contributed several leads for further study that will increase the understanding of caries pathogenesis and guide therapeutic development.

A communal catalogue reveals Earth's multiscale microbial diversity.

Our growing awareness of the microbial world's importance and diversity contrasts starkly with our limited understanding of its fundamental structure. Despite recent advances in DNA sequencing, a lack of standardized protocols and common analytical frameworks impedes comparisons among studies, hindering the development of global inferences about microbial life on Earth. Here we present a meta-analysis of microbial community samples collected by hundreds of researchers for the Earth Microbiome Project. Coordinated protocols and new analytical methods, particularly the use of exact sequences instead of clustered operational taxonomic units, enable bacterial and archaeal ribosomal RNA gene sequences to be followed across multiple studies and allow us to explore patterns of diversity at an unprecedented scale. The result is both a reference database giving global context to DNA sequence data and a framework for incorporating data from future studies, fostering increasingly complete characterization of Earth's microbial diversity.

Learning representations of microbe-metabolite interactions.

Integrating multiomics datasets is critical for microbiome research; however, inferring interactions across omics datasets has multiple statistical challenges. We solve this problem by using neural networks (https://github.com/biocore/mmvec) to estimate the conditional probability that each molecule is present given the presence of a specific microorganism. We show with known environmental (desert soil biocrust wetting) and clinical (cystic fibrosis lung) examples, our ability to recover microbe-metabolite relationships, and demonstrate how the method can discover relationships between microbially produced metabolites and inflammatory bowel disease.

Impacts of the Human Gut Microbiome on Therapeutics.

The human microbiome contains a vast source of genetic and biochemical variation, and its impacts on therapeutic responses are just beginning to be understood. This expanded understanding is especially important because the human microbiome differs far more among different people than does the human genome, and it is also dramatically easier to change. Here, we describe some of the major factors driving differences in the human microbiome among individuals and populations. We then describe some of the many ways in which gut microbes modify the action of specific chemotherapeutic agents, including nonsteroidal anti-inflammatory drugs and cardiac glycosides, and outline the potential of fecal microbiota transplant as a therapeutic. Intriguingly, microbes also alter how hosts respond to therapeutic agents through various pathways acting at distal sites. Finally, we discuss some of the computational and practical issues surrounding use of the microbiome to stratify individuals for drug response, and we envision a future where the microbiome will be modified to increase everyone's potential to benefit from therapy.

Red Sea SAR11 and Prochlorococcus Single-Cell Genomes Reflect Globally Distributed Pangenomes.

Evidence suggests many marine bacteria are cosmopolitan, with widespread but sparse strains poised to seed abundant populations under conducive growth conditions. However, studies supporting this "microbial seed bank" hypothesis have analyzed taxonomic marker genes rather than whole genomes/metagenomes, leaving open the possibility that disparate ocean regions harbor endemic gene content. The Red Sea is isolated geographically from the rest of the ocean and has a combination of high irradiance, high temperature, and high salinity that is unique among the oceans; we therefore asked whether it harbors endemic gene content. We sequenced and assembled single-cell genomes of 21 SAR11 (subclades Ia, Ib, Id, and II) and 5 Prochlorococcus (ecotype HLII) samples from the Red Sea and combined them with globally sourced reference genomes to cluster genes into ortholog groups (OGs). Ordination of OG composition could distinguish clades, including phylogenetically cryptic Prochlorococcus ecotypes LLII and LLIII. Compared with reference genomes, 1% of Prochlorococcus and 17% of SAR11 OGs were unique to the Red Sea genomes (RS-OGs). Most (83%) RS-OGs had no annotated function, but 65% of RS-OGs were expressed in diel Red Sea metatranscriptomes, suggesting they are functional. Searching Tara Oceans metagenomes, RS-OGs were as likely to be found as non-RS-OGs; nevertheless, Red Sea and other warm samples could be distinguished from cooler samples using the relative abundances of OGs. The results suggest that the prevalence of OGs in these surface ocean bacteria is largely cosmopolitan, with differences in population metagenomes manifested by differences in relative abundance rather than complete presence/absence of OGs.IMPORTANCE Studies have shown that as we sequence seawater from a selected environment deeper and deeper, we approach finding every bacterial taxon known for the ocean as a whole. However, such studies have focused on taxonomic marker genes rather than on whole genomes, raising the possibility that the lack of endemism results from the method of investigation. We took a geographically isolated water body, the Red Sea, and sequenced single cells from it. We compared those single-cell genomes to available genomes from around the ocean and to ocean-spanning metagenomes. We showed that gene ortholog groups found in Red Sea genomes but not in other genomes are nevertheless common across global ocean metagenomes. These results suggest that Baas Becking's hypothesis "everything is everywhere, but the environment selects" also applies to gene ortholog groups. This widely dispersed functional diversity may give oceanic microbial communities the functional capacity to respond rapidly to changing conditions.

Immunization with a heat-killed preparation of the environmental bacterium Mycobacterium vaccae promotes stress resilience in mice.

The prevalence of inflammatory diseases is increasing in modern urban societies. Inflammation increases risk of stress-related pathology; consequently, immunoregulatory or antiinflammatory approaches may protect against negative stress-related outcomes. We show that stress disrupts the homeostatic relationship between the microbiota and the host, resulting in exaggerated inflammation. Repeated immunization with a heat-killed preparation of Mycobacterium vaccae, an immunoregulatory environmental microorganism, reduced subordinate, flight, and avoiding behavioral responses to a dominant aggressor in a murine model of chronic psychosocial stress when tested 1-2 wk following the final immunization. Furthermore, immunization with M. vaccae prevented stress-induced spontaneous colitis and, in stressed mice, induced anxiolytic or fear-reducing effects as measured on the elevated plus-maze, despite stress-induced gut microbiota changes characteristic of gut infection and colitis. Immunization with M. vaccae also prevented stress-induced aggravation of colitis in a model of inflammatory bowel disease. Depletion of regulatory T cells negated protective effects of immunization with M. vaccae on stress-induced colitis and anxiety-like or fear behaviors. These data provide a framework for developing microbiome- and immunoregulation-based strategies for prevention of stress-related pathologies.

Mass Spectrometry-Based Chemical Cartography of a Cardiac Parasitic Infection.

Trypanosoma cruzi parasites are the causative agents of Chagas disease, a leading infectious form of heart failure whose pathogenesis is still not fully characterized. In this work, we applied untargeted liquid chromatography-tandem mass spectrometry to heart sections from T. cruzi-infected and uninfected mice. We combined molecular networking and three-dimensional modeling to generate chemical cartographical heart models. This approach revealed for the first time preferential parasite localization to the base of the heart and regiospecific distributions of nucleoside derivatives and eicosanoids, which we correlated to tissue-damaging immune responses. We further detected novel cardiac chemical signatures related to the severity and ultimate outcome of the infection. These signatures included differential representation of higher- vs lower-molecular-weight carnitine and phosphatidylcholine family members in specific cardiac regions of mice infected with lethal or nonlethal T. cruzi strains and doses. Overall, this work provides new insights into Chagas disease pathogenesis and presents an analytical chemistry approach that can be broadly applied to the study of host-microbe interactions.

The Microbiome in Posttraumatic Stress Disorder and Trauma-Exposed Controls: An Exploratory Study.

OBJECTIVE: Inadequate immunoregulation and elevated inflammation may be risk factors for posttraumatic stress disorder (PTSD), and microbial inputs are important determinants of immunoregulation; however, the association between the gut microbiota and PTSD is unknown. This study investigated the gut microbiome in a South African sample of PTSD-affected individuals and trauma-exposed (TE) controls to identify potential differences in microbial diversity or microbial community structure. METHODS: The Clinician-Administered PTSD Scale for DSM-5 was used to diagnose PTSD according to Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition criteria. Microbial DNA was extracted from stool samples obtained from 18 individuals with PTSD and 12 TE control participants. Bacterial 16S ribosomal RNA gene V3/V4 amplicons were generated and sequenced. Microbial community structure, α-diversity, and ß-diversity were analyzed; random forest analysis was used to identify associations between bacterial taxa and PTSD. RESULTS: There were no differences between PTSD and TE control groups in α- or ß-diversity measures (e.g., α-diversity: Shannon index, t = 0.386, p = .70; ß-diversity, on the basis of analysis of similarities: Bray-Curtis test statistic = -0.033, p = .70); however, random forest analysis highlighted three phyla as important to distinguish PTSD status: Actinobacteria, Lentisphaerae, and Verrucomicrobia. Decreased total abundance of these taxa was associated with higher Clinician-Administered PTSD Scale scores (r = -0.387, p = .035). CONCLUSIONS: In this exploratory study, measures of overall microbial diversity were similar among individuals with PTSD and TE controls; however, decreased total abundance of Actinobacteria, Lentisphaerae, and Verrucomicrobia was associated with PTSD status.

Parkinson's disease and Parkinson's disease medications have distinct signatures of the gut microbiome.

BACKGROUND: There is mounting evidence for a connection between the gut and Parkinson's disease (PD). Dysbiosis of gut microbiota could explain several features of PD. OBJECTIVE: The objective of this study was to determine if PD involves dysbiosis of gut microbiome, disentangle effects of confounders, and identify candidate taxa and functional pathways to guide research. METHODS: A total of 197 PD cases and 130 controls were studied. Microbial composition was determined by 16S rRNA gene sequencing of DNA extracted from stool. Metadata were collected on 39 potential confounders including medications, diet, gastrointestinal symptoms, and demographics. Statistical analyses were conducted while controlling for potential confounders and correcting for multiple testing. We tested differences in the overall microbial composition, taxa abundance, and functional pathways. RESULTS: Independent microbial signatures were detected for PD (P = 4E-5), participants' region of residence within the United States (P = 3E-3), age (P = 0.03), sex (P = 1E-3), and dietary fruits/vegetables (P = 0.01). Among patients, independent signals were detected for catechol-O-methyltransferase-inhibitors (P = 4E-4), anticholinergics (P = 5E-3), and possibly carbidopa/levodopa (P = 0.05). We found significantly altered abundances of the Bifidobacteriaceae, Christensenellaceae, [Tissierellaceae], Lachnospiraceae, Lactobacillaceae, Pasteurellaceae, and Verrucomicrobiaceae families. Functional predictions revealed changes in numerous pathways, including the metabolism of plant-derived compounds and xenobiotics degradation. CONCLUSION: PD is accompanied by dysbiosis of gut microbiome. Results coalesce divergent findings of prior studies, reveal altered abundance of several taxa, nominate functional pathways, and demonstrate independent effects of PD medications on the microbiome. The findings provide new leads and testable hypotheses on the pathophysiology and treatment of PD. © 2017 International Parkinson and Movement Disorder Society.

A large scale prediction of bacteriocin gene blocks suggests a wide functional spectrum for bacteriocins.

BACKGROUND: Bacteriocins are peptide-derived molecules produced by bacteria, whose recently-discovered functions include virulence factors and signaling molecules as well as their better known roles as antibiotics. To date, close to five hundred bacteriocins have been identified and classified. Recent discoveries have shown that bacteriocins are highly diverse and widely distributed among bacterial species. Given the heterogeneity of bacteriocin compounds, many tools struggle with identifying novel bacteriocins due to their vast sequence and structural diversity. Many bacteriocins undergo post-translational processing or modifications necessary for the biosynthesis of the final mature form. Enzymatic modification of bacteriocins as well as their export is achieved by proteins whose genes are often located in a discrete gene cluster proximal to the bacteriocin precursor gene, referred to as context genes in this study. Although bacteriocins themselves are structurally diverse, context genes have been shown to be largely conserved across unrelated species. METHODS: Using this knowledge, we set out to identify new candidates for context genes which may clarify how bacteriocins are synthesized, and identify new candidates for bacteriocins that bear no sequence similarity to known toxins. To achieve these goals, we have developed a software tool, Bacteriocin Operon and gene block Associator (BOA) that can identify homologous bacteriocin associated gene blocks and predict novel ones. BOA generates profile Hidden Markov Models from the clusters of bacteriocin context genes, and uses them to identify novel bacteriocin gene blocks and operons. RESULTS AND CONCLUSIONS: We provide a novel dataset of predicted bacteriocins and context genes. We also discover that several phyla have a strong preference for bacteriocin genes, suggesting distinct functions for this group of molecules. SOFTWARE AVAILABILITY: https://github.com/idoerg/BOA.

SCOPE++: sequence classification of homoPolymer emissions.

BACKGROUND: mRNA polyadenylation, the addition of a poly(A) tail to the 3'-end of pre-mRNA, is a process critical to gene expression and regulation in eukaryotes. To understand the molecular mechanisms governing polyadenylation and other relevant biological processes, it is important to identify these poly(A) tails accurately in transcriptome sequencing data and differentiate them from artificial adapter sequences added in the sequencing process. But the annotation of these tails is complicated by the presence of sequencing errors and post-transcriptional modifications. While determining that a tail is present in a given transcript fragment is straight-forward, these obfuscations make the problem of boundary identification a challenge; conventional seed-and-extend algorithms struggle to accurately identify these poly(A) tail end-points. Further, all existing tools that we are aware of focus exclusively on the trimming of poly(A) tails, failing to provide the detailed information needed for studying the polyadenylation process. RESULTS: We have created SCOPE++, an open-source tool for finding the precise border of poly(A) tails and other homopolymers in raw mRNA sequence reads. Based on a Hidden Markov Model (HMM) approach, SCOPE++ accurately identifies specific homopolymer sequences in error-prone EST/cDNA data or RNA-Seq data at a speed appropriate for large sequence sets. CONCLUSIONS: We demonstrate that our tool can precisely identify poly(A) tails with near perfect accuracy at the speed required for high-throughput applications, providing a valuable resource for polyadenylation research.

Meta-analysis of the human gut microbiome uncovers shared and distinct microbial signatures between diseases.

Microbiome studies have revealed gut microbiota's potential impact on complex diseases. However, many studies often focus on one disease per cohort. We developed a meta-analysis workflow for gut microbiome profiles and analyzed shotgun metagenomic data covering 11 diseases. Using interpretable machine learning and differential abundance analysis, our findings reinforce the generalization of binary classifiers for Crohn's disease (CD) and colorectal cancer (CRC) to hold-out cohorts and highlight the key microbes driving these classifications. We identified high microbial similarity in disease pairs like CD vs ulcerative colitis (UC), CD vs CRC, Parkinson's disease vs type 2 diabetes (T2D), and schizophrenia vs T2D. We also found strong inverse correlations in Alzheimer's disease vs CD and UC. These findings detected by our pipeline provide valuable insights into these diseases.

Meta-analysis of the human gut microbiome uncovers shared and distinct microbial signatures between diseases.

Microbiome studies have revealed gut microbiota's potential impact on complex diseases. However, many studies often focus on one disease per cohort. We developed a meta-analysis workflow for gut microbiome profiles and analyzed shotgun metagenomic data covering 11 diseases. Using interpretable machine learning and differential abundance analysis, our findings reinforce the generalization of binary classifiers for Crohn's disease (CD) and colorectal cancer (CRC) to hold-out cohorts and highlight the key microbes driving these classifications. We identified high microbial similarity in disease pairs like CD vs ulcerative colitis (UC), CD vs CRC, Parkinson's disease vs type 2 diabetes (T2D), and schizophrenia vs T2D. We also found strong inverse correlations in Alzheimer's disease vs CD and UC. These findings, detected by our pipeline, provide valuable insights into these diseases. IMPORTANCE: Assessing disease similarity is an essential initial step preceding a disease-based approach for drug repositioning. Our study provides a modest first step in underscoring the potential of integrating microbiome insights into the disease similarity assessment. Recent microbiome research has predominantly focused on analyzing individual diseases to understand their unique characteristics, which by design excludes comorbidities in individuals. We analyzed shotgun metagenomic data from existing studies and identified previously unknown similarities between diseases. Our research represents a pioneering effort that utilizes both interpretable machine learning and differential abundance analysis to assess microbial similarity between diseases.

Learning sequence, structure, and function representations of proteins with language models.

The sequence-structure-function relationships that ultimately generate the diversity of extant observed proteins is complex, as proteins bridge the gap between multiple informational and physical scales involved in nearly all cellular processes. One limitation of existing protein annotation databases such as UniProt is that less than 1% of proteins have experimentally verified functions, and computational methods are needed to fill in the missing information. Here, we demonstrate that a multi-aspect framework based on protein language models can learn sequence-structure-function representations of amino acid sequences, and can provide the foundation for sensitive sequence-structure-function aware protein sequence search and annotation. Based on this model, we introduce a multi-aspect information retrieval system for proteins, Protein-Vec, covering sequence, structure, and function aspects, that enables computational protein annotation and function prediction at tree-of-life scales.

BIRDMAn: A Bayesian differential abundance framework that enables robust inference of host-microbe associations.

Quantifying the differential abundance (DA) of specific taxa among experimental groups in microbiome studies is challenging due to data characteristics (e.g., compositionality, sparsity) and specific study designs (e.g., repeated measures, meta-analysis, cross-over). Here we present BIRDMAn (Bayesian Inferential Regression for Differential Microbiome Analysis), a flexible DA method that can account for microbiome data characteristics and diverse experimental designs. Simulations show that BIRDMAn models are robust to uneven sequencing depth and provide a >20-fold improvement in statistical power over existing methods. We then use BIRDMAn to identify antibiotic-mediated perturbations undetected by other DA methods due to subject-level heterogeneity. Finally, we demonstrate how BIRDMAn can construct state-of-the-art cancer-type classifiers using The Cancer Genome Atlas (TCGA) dataset, with substantial accuracy improvements over random forests and existing DA tools across multiple sequencing centers. Collectively, BIRDMAn extracts more informative biological signals while accounting for study-specific experimental conditions than existing approaches.

Protein remote homology detection and structural alignment using deep learning.

Exploiting sequence-structure-function relationships in biotechnology requires improved methods for aligning proteins that have low sequence similarity to previously annotated proteins. We develop two deep learning methods to address this gap, TM-Vec and DeepBLAST. TM-Vec allows searching for structure-structure similarities in large sequence databases. It is trained to accurately predict TM-scores as a metric of structural similarity directly from sequence pairs without the need for intermediate computation or solution of structures. Once structurally similar proteins have been identified, DeepBLAST can structurally align proteins using only sequence information by identifying structurally homologous regions between proteins. It outperforms traditional sequence alignment methods and performs similarly to structure-based alignment methods. We show the merits of TM-Vec and DeepBLAST on a variety of datasets, including better identification of remotely homologous proteins compared with state-of-the-art sequence alignment and structure prediction methods.

The Local Tumor Microbiome Is Associated with Survival in Late-Stage Colorectal Cancer Patients.

The gut microbiome is associated with survival in colorectal cancer. Single organisms have been identified as markers of poor prognosis. However, in situ imaging of tumors demonstrate a polymicrobial tumor-associated community. To understand the role of these polymicrobial communities in survival, we conducted a nested case-control study in late-stage cancer patients undergoing resection for primary adenocarcinoma. The microbiome of paired tumor and adjacent normal tissue samples was profiled using 16S rRNA sequencing. We found a consistent difference in the microbiome between paired tumor and adjacent tissue, despite strong individual microbial identities. Furthermore, a larger difference between normal and tumor tissue was associated with prognosis: patients with shorter survival had a larger difference between normal and tumor tissue. Within the tumor tissue, we identified a 39-member community statistic associated with survival; for every log2-fold increase in this value, an individual's odds of survival increased by 20% (odds ratio survival 1.20; 95% confidence interval = 1.04 to 1.33). Our results suggest that a polymicrobial tumor-specific microbiome is associated with survival in late-stage colorectal cancer patients. IMPORTANCE Microbiome studies in colorectal cancer (CRC) have primarily focused on the role of single organisms in cancer progression. Recent work has identified specific organisms throughout the intestinal tract, which may affect survival; however, the results are inconsistent. We found differences between the tumor microbiome and the microbiome of the rest of the intestine in patients, and the magnitude of this difference was associated with survival, or, the more like a healthy gut a tumor looked, the better a patient's prognosis. Our results suggest that future microbiome-based interventions to affect survival in CRC will need to target the tumor community.

Greengenes2 unifies microbial data in a single reference tree.

Studies using 16S rRNA and shotgun metagenomics typically yield different results, usually attributed to PCR amplification biases. We introduce Greengenes2, a reference tree that unifies genomic and 16S rRNA databases in a consistent, integrated resource. By inserting sequences into a whole-genome phylogeny, we show that 16S rRNA and shotgun metagenomic data generated from the same samples agree in principal coordinates space, taxonomy and phenotype effect size when analyzed with the same tree.

Multi-level analysis of the gut-brain axis shows autism spectrum disorder-associated molecular and microbial profiles.

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by heterogeneous cognitive, behavioral and communication impairments. Disruption of the gut-brain axis (GBA) has been implicated in ASD although with limited reproducibility across studies. In this study, we developed a Bayesian differential ranking algorithm to identify ASD-associated molecular and taxa profiles across 10 cross-sectional microbiome datasets and 15 other datasets, including dietary patterns, metabolomics, cytokine profiles and human brain gene expression profiles. We found a functional architecture along the GBA that correlates with heterogeneity of ASD phenotypes, and it is characterized by ASD-associated amino acid, carbohydrate and lipid profiles predominantly encoded by microbial species in the genera Prevotella, Bifidobacterium, Desulfovibrio and Bacteroides and correlates with brain gene expression changes, restrictive dietary patterns and pro-inflammatory cytokine profiles. The functional architecture revealed in age-matched and sex-matched cohorts is not present in sibling-matched cohorts. We also show a strong association between temporal changes in microbiome composition and ASD phenotypes. In summary, we propose a framework to leverage multi-omic datasets from well-defined cohorts and investigate how the GBA influences ASD.