Search

Longitudinal studies are crucial for understanding complex microbiome dynamics and their link to health. We introduce TEMPoral TEnsor Decomposition (TEMPTED), a time-informed dimensionality reduction method for high-dimensional longitudinal data that treats time as a continuous variable, effectively characterizing temporal information and handling varying temporal sampling. TEMPTED captures key microbial dynamics, facilitates beta-diversity analysis, and enhances reproducibility by transferring learned representations to new data. In simulations, it achieves 90% accuracy in phenotype classification, significantly outperforming existing methods. In real data, TEMPTED identifies vaginal microbial markers linked to term and preterm births, demonstrating robust performance across datasets and sequencing platforms.

In this study, we aimed to understand the potential role of the gut microbiome in the development of Alzheimers disease (AD). We took a multi-faceted approach to investigate this relationship. Urine metabolomics were examined in individuals with AD and controls, revealing decreased formate and fumarate concentrations in AD. Additionally, we utilised whole-genome sequencing (WGS) data obtained from a separate group of individuals with AD and controls. This information allowed us to create and investigate host-microbiome personalised whole-body metabolic models. Notably, AD individuals displayed diminished formate microbial secretion in these models. Additionally, we identified specific reactions responsible for the production of formate in the host, and interestingly, these reactions were linked to genes that have correlations with AD. This study suggests formate as a possible early AD marker and highlights genetic and microbiome contributions to its production. The reduced formate secretion and its genetic associations point to a complex connection between gut microbiota and AD. This holistic understanding might pave the way for novel diagnostic and therapeutic avenues in AD management.

Microbial breakdown of organic matter is one of the most important processes on Earth, yet the controls of decomposition are poorly understood. Here we track 36 terrestrial human cadavers in three locations and show that a phylogenetically distinct, interdomain microbial network assembles during decomposition despite selection effects of location, climate and season. We generated a metagenome-assembled genome library from cadaver-associated soils and integrated it with metabolomics data to identify links between taxonomy and function. This universal network of microbial decomposers is characterized by cross-feeding to metabolize labile decomposition products. The key bacterial and fungal decomposers are rare across non-decomposition environments and appear unique to the breakdown of terrestrial decaying flesh, including humans, swine, mice and cattle, with insects as likely important vectors for dispersal. The observed lockstep of microbial interactions further underlies a robust microbial forensic tool with the potential to aid predictions of the time since death.

INTRODUCTION: The ketogenic diet (KD) is an intriguing therapeutic candidate for Alzheimers disease (AD) given its protective effects against metabolic dysregulation and seizures. Gut microbiota are essential for KD-mediated neuroprotection against seizures as well as modulation of bile acids, which play a major role in cholesterol metabolism. These relationships motivated our analysis of gut microbiota and metabolites related to cognitive status following a controlled KD intervention compared with a low-fat-diet intervention. METHODS: Prediabetic adults, either with mild cognitive impairment (MCI) or cognitively normal (CN), were placed on either a low-fat American Heart Association diet or high-fat modified Mediterranean KD (MMKD) for 6 weeks; then, after a 6-week washout period, they crossed over to the alternate diet. We collected stool samples for shotgun metagenomics and untargeted metabolomics at five time points to investigate individuals microbiome and metabolome throughout the dietary interventions. RESULTS: Participants with MCI on the MMKD had lower levels of GABA-producing microbes Alistipes sp. CAG:514 and GABA, and higher levels of GABA-regulating microbes Akkermansia muciniphila. MCI individuals with curcumin in their diet had lower levels of bile salt hydrolase-containing microbes and an altered bile acid pool, suggesting reduced gut motility. DISCUSSION: Our results suggest that the MMKD may benefit adults with MCI through modulation of GABA levels and gut-transit time.

Colorectal cancer (CRC) is driven by genomic alterations in concert with dietary influences, with the gut microbiome implicated as an effector in disease development and progression. While meta-analyses have provided mechanistic insight into patients with CRC, study heterogeneity has limited causal associations. Using multi-omics studies on genetically controlled cohorts of mice, we identify diet as the major driver of microbial and metabolomic differences, with reductions in α diversity and widespread changes in cecal metabolites seen in high-fat diet (HFD)-fed mice. In addition, non-classic amino acid conjugation of the bile acid cholic acid (AA-CA) increased with HFD. We show that AA-CAs impact intestinal stem cell growth and demonstrate that Ileibacterium valens and Ruminococcus gnavus are able to synthesize these AA-CAs. This multi-omics dataset implicates diet-induced shifts in the microbiome and the metabolome in disease progression and has potential utility in future diagnostic and therapeutic developments.

Periodontal disease is more common in individuals with rheumatoid arthritis (RA) who have detectable anti-citrullinated protein antibodies (ACPAs), implicating oral mucosal inflammation in RA pathogenesis. Here, we performed paired analysis of human and bacterial transcriptomics in longitudinal blood samples from RA patients. We found that patients with RA and periodontal disease experienced repeated oral bacteremias associated with transcriptional signatures of ISG15+HLADRhi and CD48highS100A2pos monocytes, recently identified in inflamed RA synovia and blood of those with RA flares. The oral bacteria observed transiently in blood were broadly citrullinated in the mouth, and their in situ citrullinated epitopes were targeted by extensively somatically hypermutated ACPAs encoded by RA blood plasmablasts. Together, these results suggest that (i) periodontal disease results in repeated breaches of the oral mucosa that release citrullinated oral bacteria into circulation, which (ii) activate inflammatory monocyte subsets that are observed in inflamed RA synovia and blood of RA patients with flares and (iii) activate ACPA B cells, thereby promoting affinity maturation and epitope spreading to citrullinated human antigens.

Diabetes represents a spectrum of disease in which metabolic dysfunction damages multiple organ systems including liver, kidneys and peripheral nerves1,2. Although the onset and progression of these co-morbidities are linked with insulin resistance, hyperglycaemia and dyslipidaemia3-7, aberrant non-essential amino acid (NEAA) metabolism also contributes to the pathogenesis of diabetes8-10. Serine and glycine are closely related NEAAs whose levels are consistently reduced in patients with metabolic syndrome10-14, but the mechanistic drivers and downstream consequences of this metabotype remain unclear. Low systemic serine and glycine are also emerging as a hallmark of macular and peripheral nerve disorders, correlating with impaired visual acuity and peripheral neuropathy15,16. Here we demonstrate that aberrant serine homeostasis drives serine and glycine deficiencies in diabetic mice, which can be diagnosed with a serine tolerance test that quantifies serine uptake and disposal. Mimicking these metabolic alterations in young mice by dietary serine or glycine restriction together with high fat intake markedly accelerates the onset of small fibre neuropathy while reducing adiposity. Normalization of serine by dietary supplementation and mitigation of dyslipidaemia with myriocin both alleviate neuropathy in diabetic mice, linking serine-associated peripheral neuropathy to sphingolipid metabolism. These findings identify systemic serine deficiency and dyslipidaemia as novel risk factors for peripheral neuropathy that may be exploited therapeutically.

BackgroundMicroorganisms such as coliform-forming bacteria are commonly used to assess freshwater quality for drinking and recreational use. However, such organisms do not exist in isolation; they exist within the context of dynamic, interactive microbial communities which vary through space and time. Elucidating spatiotemporal microbial dynamics is imperative for discriminating robust community changes from ephemeral ecological trends, and for improving our overall understanding of the relationship between microbial communities and ecosystem health. We conducted a seven-year (2013-2019) microbial time-series investigation in the Chicago Area Waterways (CAWS): an urban river system which, in 2016, experienced substantial upgrades to disinfection processes at two wastewater reclamation plants (WRPs) that discharge into the CAWS and improved stormwater capture, to improve river water quality and reduce flooding. Using culture-independent and culture-dependent approaches, we compared CAWS microbial ecology before and after the intervention.ResultsExaminations of time-resolved beta distances between WRP-adjacent sites showed that community similarity measures were often consistent with the spatial orientation of site locations to one another and to the WRP outfalls. Fecal coliform results suggested that upgrades reduced coliform-associated bacteria in the effluent and the downstream river community. However, examinations of whole community changes through time suggest that the upgrades did little to affect overall riverine community dynamics, which instead were overwhelmingly driven by yearly patterns consistent with seasonality.ConclusionsThis study presents a systematic effort to combine 16S rRNA gene amplicon sequencing with traditional culture-based methods to evaluate the influence of treatment innovations and systems upgrades on the microbiome of the Chicago Area Waterway System, representing the longest and most comprehensive characterization of the microbiome of an urban waterway yet attempted. We found that the systems upgrades were successful in improving specific water quality measures immediately downstream of wastewater outflows. Additionally, we found that the implementation of the water quality improvement measures to the river system did not disrupt the overall dynamics of the downstream microbial community, which remained heavily influenced by seasonal trends. Such results emphasize the dynamic nature of microbiomes in open environmental systems such as the CAWS, but also suggest that the seasonal oscillations remain consistent even when perturbed.

Human untargeted metabolomics studies annotate only ~10% of molecular features. We introduce reference-data-driven analysis to match metabolomics tandem mass spectrometry (MS/MS) data against metadata-annotated source data as a pseudo-MS/MS reference library. Applying this approach to food source data, we show that it increases MS/MS spectral usage 5.1-fold over conventional structural MS/MS library matches and allows empirical assessment of dietary patterns from untargeted data.

Despite advances in sequencing, lack of standardization makes comparisons across studies challenging and hampers insights into the structure and function of microbial communities across multiple habitats on a planetary scale. Here we present a multi-omics analysis of a diverse set of 880 microbial community samples collected for the Earth Microbiome Project. We include amplicon (16S, 18S, ITS) and shotgun metagenomic sequence data, and untargeted metabolomics data (liquid chromatography-tandem mass spectrometry and gas chromatography mass spectrometry). We used standardized protocols and analytical methods to characterize microbial communities, focusing on relationships and co-occurrences of microbially related metabolites and microbial taxa across environments, thus allowing us to explore diversity at extraordinary scale. In addition to a reference database for metagenomic and metabolomic data, we provide a framework for incorporating additional studies, enabling the expansion of existing knowledge in the form of an evolving community resource. We demonstrate the utility of this database by testing the hypothesis that every microbe and metabolite is everywhere but the environment selects. Our results show that metabolite diversity exhibits turnover and nestedness related to both microbial communities and the environment, whereas the relative abundances of microbially related metabolites vary and co-occur with specific microbial consortia in a habitat-specific manner. We additionally show the power of certain chemistry, in particular terpenoids, in distinguishing Earth's environments (for example, terrestrial plant surfaces and soils, freshwater and marine animal stool), as well as that of certain microbes including Conexibacter woesei (terrestrial soils), Haloquadratum walsbyi (marine deposits) and Pantoea dispersa (terrestrial plant detritus). This Resource provides insight into the taxa and metabolites within microbial communities from diverse habitats across Earth, informing both microbial and chemical ecology, and provides a foundation and methods for multi-omics microbiome studies of hosts and the environment.

Cancer-microbe associations have been explored for centuries, but cancer-associated fungi have rarely been examined. Here, we comprehensively characterize the cancer mycobiome within 17,401 patient tissue, blood, and plasma samples across 35 cancer types in four independent cohorts. We report fungal DNA and cells at low abundances across many major human cancers, with differences in community compositions that differ among cancer types, even when accounting for technical background. Fungal histological staining of tissue microarrays supported intratumoral presence and frequent spatial association with cancer cells and macrophages. Comparing intratumoral fungal communities with matched bacteriomes and immunomes revealed co-occurring bi-domain ecologies, often with permissive, rather than competitive, microenvironments and distinct immune responses. Clinically focused assessments suggested prognostic and diagnostic capacities of the tissue and plasma mycobiomes, even in stage I cancers, and synergistic predictive performance with bacteriomes.

Microbial communities contain a broad phylogenetic diversity of organisms; however, the majority of methods center on describing bacteria and archaea. Fungi are important symbionts in many ecosystems and are potentially important members of the human microbiome, beyond those that can cause disease. To expand our analysis of microbial communities to include data from the fungal internal transcribed spacer (ITS) region, five candidate DNA extraction kits were compared against our standardized protocol for describing bacteria and archaea using 16S rRNA gene amplicon- and shotgun metagenomics sequencing. The results are presented considering a diverse panel of host-associated and environmental sample types and comparing the cost, processing time, well-to-well contamination, DNA yield, limit of detection and microbial community composition among protocols. Across all criteria, the MagMAX Microbiome kit was found to perform best. The PowerSoil Pro kit performed comparably but with increased cost per sample and overall processing time. The Zymo MagBead, NucleoMag Food and Norgen Stool kits were included.

Compositional oscillations of the gut microbiome are essential for normal peripheral circadian rhythms, both of which are disrupted in diet-induced obesity (DIO). Although time-restricted feeding (TRF) maintains circadian synchrony and protects against DIO, its impact on the dynamics of the cecal gut microbiome is modest. Thus, other regions of the gut, particularly the ileum, the nexus for incretin and bile acid signaling, may play an important role in entraining peripheral circadian rhythms. We demonstrate the effect of diet and feeding rhythms on the ileal microbiome composition and transcriptome in mice. The dynamic rhythms of ileal microbiome composition and transcriptome are dampened in DIO. TRF partially restores diurnal rhythms of the ileal microbiome and transcriptome, increases GLP-1 release, and alters the ileal bile acid pool and farnesoid X receptor (FXR) signaling, which could explain how TRF exerts its metabolic benefits. Finally, we provide a web resource for exploration of ileal microbiome and transcriptome circadian data.

Microbiome data have several specific characteristics (sparsity and compositionality) that introduce challenges in data analysis. The integration of prior information regarding the data structure, such as phylogenetic structure and repeated-measure study designs, into analysis, is an effective approach for revealing robust patterns in microbiome data. Past methods have addressed some but not all of these challenges and features: for example, robust principal-component analysis (RPCA) addresses sparsity and compositionality; compositional tensor factorization (CTF) addresses sparsity, compositionality, and repeated measure study designs; and UniFrac incorporates phylogenetic information. Here we introduce a strategy of incorporating phylogenetic information into RPCA and CTF. The resulting methods, phylo-RPCA, and phylo-CTF, provide substantial improvements over state-of-the-art methods in terms of discriminatory power of underlying clustering ranging from the mode of delivery to adult human lifestyle. We demonstrate quantitatively that the addition of phylogenetic information improves effect size and classification accuracy in both data-driven simulated data and real microbiome data. IMPORTANCE Microbiome data analysis can be difficult because of particular data features, some unavoidable and some due to technical limitations of DNA sequencing instruments. The first step in many analyses that ultimately reveals patterns of similarities and differences among sets of samples (e.g., separating samples from sick and healthy people or samples from seawater versus soil) is calculating the difference between each pair of samples. We introduce two new methods to calculate these differences that combine features of past methods, specifically being able to take into account the principles that most types of microbes are not in most samples (sparsity), that abundances are relative rather than absolute (compositionality), and that all microbes have a shared evolutionary history (phylogeny). We show using simulated and real data that our new methods provide improved classification accuracy of ordinal sample clusters and increased effect size between sample groups on beta-diversity distances.

Abstract While the study of the tumor microbiome and its effects on cancer biology has expanded considerably over the last few years, most of this research focused on bacteria and viruses, leaving behind the fungal kingdom. Recently, a few studies have demonstrated that specific fungi may promote tumor progression, stressing the importance of comprehensively studying the tumor mycobiome and its effects. To address this, we have characterized the mycobiome in 1,183 human tumors and their adjacent tissues, originating from eight major solid tumor types. Staining and imaging demonstrated the presence of fungi in both cancer and immune cells, with tumor-type specific distribution patterns. Quantitative PCR of the fungal 5.8s rDNA revealed the presence of fungal DNA in all tumor types. To characterize the tumor mycobiome and address potential contamination during tissue handling and processing, we subjected all samples, as well as 295 negative controls of different types, to sequencing of the ITS2 region that is situated between fungal rRNA genes. We found cancer-type specific mycobial signatures with relatively high similarity between tumors and their adjacent tissues. While the fungal mycobiome had a lower species richness as compared to the bacterial microbiome of the same tumors, fungi showed significant co-occurrences with specific bacteria, suggesting the existence of ecological niches within the tumors. We also found significant correlations with clinical parameters such as patient’s age, tumor stage, progression-free survival, overall survival, and response to immune checkpoint blockade therapy. Characterization of the tumor mycobiome may add a biologically relevant, previously overlooked, component to be considered in the study of cancer, including its effects on tumor initiation, progression, diagnosis, and response to therapy. Citation Format: Lian Narunsky Haziza, Gregory D. Sepich-Poore, Ilana Livyatan, Omer Asraf, Cameron Martino, Deborah Nejman, Nancy Gavert, Jason E. Stajich, Guy Amit, Antonio González, Stephen Wandro, Gili Perry, Ruthie Ariel, Arnon Meltser, Justin P. Shaffer, Qiyun Zhu, Nora Balint-Lahat, Iris Barshack, Maya Dadian, Einav N. Gal-Yam, Sandip P. Pate, Amir Bashan, Austin D. Swafford, Yitzhak Pilpel, Rob Knight, Ravid Straussman. Pan-cancer characterization of the tumor mycobiome and its clinical effects [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3054.

Increasing data volumes on high-throughput sequencing instruments such as the NovaSeq 6000 leads to long computational bottlenecks for common metagenomics data preprocessing tasks such as adaptor and primer trimming and host removal. Here, we test whether faster recently developed computational tools (Fastp and Minimap2) can replace widely used choices (Atropos and Bowtie2), obtaining dramatic accelerations with additional sensitivity and minimal loss of specificity for these tasks. Furthermore, the taxonomic tables resulting from downstream processing provide biologically comparable results. However, we demonstrate that for taxonomic assignment, Bowtie2's specificity is still required. We suggest that periodic reevaluation of pipeline components, together with improvements to standardized APIs to chain them together, will greatly enhance the efficiency of common bioinformatics tasks while also facilitating incorporation of further optimized steps running on GPUs, FPGAs, or other architectures. We also note that a detailed exploration of available algorithms and pipeline components is an important step that should be taken before optimization of less efficient algorithms on advanced or nonstandard hardware. IMPORTANCE In shotgun metagenomics studies that seek to relate changes in microbial DNA across samples, processing the data on a computer often takes longer than obtaining the data from the sequencing instrument. Recently developed software packages that perform individual steps in the pipeline of data processing in principle offer speed advantages, but in practice they may contain pitfalls that prevent their use, for example, they may make approximations that introduce unacceptable errors in the data. Here, we show that differences in choices of these components can speed up overall data processing by 5-fold or more on the same hardware while maintaining a high degree of correctness, greatly reducing the time taken to interpret results. This is an important step for using the data in clinical settings, where the time taken to obtain the results may be critical for guiding treatment.

Periodontitis affects up to 50% of individuals worldwide, and 8.5% are diagnosed with diabetes. The high-comorbidity rate of these diseases may suggest, at least in part, a shared etiology and pathophysiology. Changes in oral microbial communities have been documented in the context of severe periodontitis and diabetes, both independently and together. However, much less is known about the early oral microbial markers of these diseases. We used a subset of the ORIGINS project dataset, which collected detailed periodontal and cardiometabolic information from 787 healthy individuals, to identify early microbial markers of periodontitis and its association with markers of cardiometabolic health. Using state-of-the-art compositional data analysis tools, we identified the log-ratio of Treponema to Corynebacterium bacteria to be a novel Microbial Indicator of Periodontitis (MIP), and found that this MIP correlates with poor periodontal health and cardiometabolic markers early in disease pathogenesis in both subgingival plaque and saliva.

Cell interactions determine phenotypes, and intercellular communication is shaped by cellular contexts such as disease state, organismal life stage, and tissue microenvironment. Single-cell technologies measure the molecules mediating cell-cell communication, and emerging computational tools can exploit these data to decipher intercellular communication. However, current methods either disregard cellular context or rely on simple pairwise comparisons between samples, thus limiting the ability to decipher complex cell-cell communication across multiple time points, levels of disease severity, or spatial contexts. Here we present Tensor-cell2cell, an unsupervised method using tensor decomposition, which deciphers context-driven intercellular communication by simultaneously accounting for multiple stages, states, or locations of the cells. To do so, Tensor-cell2cell uncovers context-driven patterns of communication associated with different phenotypic states and determined by unique combinations of cell types and ligand-receptor pairs. As such, Tensor-cell2cell robustly improves upon and extends the analytical capabilities of existing tools. We show Tensor-cell2cell can identify multiple modules associated with distinct communication processes (e.g., participating cell-cell and ligand-receptor pairs) linked to severities of Coronavirus Disease 2019 and to Autism Spectrum Disorder. Thus, we introduce an effective and easy-to-use strategy for understanding complex communication patterns across diverse conditions.

Liver damage due to chronic alcohol use is among the most prevalent liver diseases. Alcohol consumption frequency is a strong factor of microbiota variance. Here we use isotope labeled [1-13C] ethanol, metagenomics, and metatranscriptomics in ethanol-feeding and intragastric mouse models to investigate the metabolic impacts of alcohol consumption on the gut microbiota. First, we show that although stable isotope labeled [1-13C] ethanol contributes to fatty acid pools in the liver, plasma, and cecum contents of mice, there is no evidence of ethanol metabolism by gut microbiota ex vivo under anaerobic conditions. Next, we observe through metatranscriptomics that the gut microbiota responds to ethanol-feeding by activating acetate dissimilation, not by metabolizing ethanol directly. We demonstrate that blood acetate concentrations are elevated during ethanol consumption. Finally, by increasing systemic acetate levels with glyceryl triacetate supplementation, we do not observe any impact on liver disease, but do induce similar gut microbiota alterations as chronic ethanol-feeding in mice. Our results show that ethanol is not directly metabolized by the gut microbiota, and changes in the gut microbiota linked to ethanol are a side effect of elevated acetate levels. De-trending for these acetate effects may be critical for understanding gut microbiota changes that cause alcohol-related liver disease.

BackgroundThe Spacecraft Assembly Facility (SAF) at the NASA's Jet Propulsion Laboratory is the primary cleanroom facility used in the construction of some of the planetary protection (PP)-sensitive missions developed by NASA, including the Mars 2020 Perseverance Rover that launched in July 2020. SAF floor samples (n=98) were collected, over a 6-month period in 2016 prior to the construction of the Mars rover subsystems, to better understand the temporal and spatial distribution of bacterial populations (total, viable, cultivable, and spore) in this unique cleanroom.ResultsCleanroom samples were examined for total (living and dead) and viable (living only) microbial populations using molecular approaches and cultured isolates employing the traditional NASA standard spore assay (NSA), which predominantly isolated spores. The 130 NSA isolates were represented by 16 bacterial genera, of which 97% were identified as spore-formers via Sanger sequencing. The most spatially abundant isolate was Bacillus subtilis, and the most temporally abundant spore-former was Virgibacillus panthothenticus. The 16S rRNA gene-targeted amplicon sequencing detected 51 additional genera not found in the NSA method. The amplicon sequencing of the samples treated with propidium monoazide (PMA), which would differentiate between viable and dead organisms, revealed a total of 54 genera: 46 viable non-spore forming genera and 8 viable spore forming genera in these samples. The microbial diversity generated by the amplicon sequencing corresponded to ~86% non-spore-formers and ~14% spore-formers. The most common spatially distributed genera were Sphinigobium, Geobacillus, and Bacillus whereas temporally distributed common genera were Acinetobacter, Geobacilllus, and Bacillus. Single-cell genomics detected 6 genera in the sample analyzed, with the most prominent being Acinetobacter.ConclusionThis study clearly established that detecting spores via NSA does not provide a complete assessment for the cleanliness of spacecraft-associated environments since it failed to detect several PP-relevant genera that were only recovered via molecular methods. This highlights the importance of a methodological paradigm shift to appropriately monitor bioburden in cleanrooms for not only the aeronautical industry but also for pharmaceutical, medical industries, etc., and the need to employ molecular sequencing to complement traditional culture-based assays. Video abstract.

BackgroundInfectious bacterial diseases exhibiting increasing resistance to antibiotics are a serious global health issue. Bacteriophage therapy is an anti-microbial alternative to treat patients with serious bacterial infections. However, the impacts to the host microbiome in response to clinical use of phage therapy are not well understood.ResultsOur paper demonstrates a largely unchanged microbiota profile during 4 weeks of phage therapy when added to systemic antibiotics in a single patient with Staphylococcus aureus device infection. Metabolomic analyses suggest potential indirect cascading ecological impacts to the host (skin) microbiome. We did not detect genomes of the three phages used to treat the patient in metagenomic samples taken from saliva, stool, and skin; however, phages were detected using endpoint-PCR in patient serum.ConclusionResults from our proof-of-principal study supports the use of bacteriophages as a microbiome-sparing approach to treat bacterial infections. Video abstract.

BackgroundSARS-CoV-2 is an RNA virus responsible for the coronavirus disease 2019 (COVID-19) pandemic. Viruses exist in complex microbial environments, and recent studies have revealed both synergistic and antagonistic effects of specific bacterial taxa on viral prevalence and infectivity. We set out to test whether specific bacterial communities predict SARS-CoV-2 occurrence in a hospital setting.MethodsWe collected 972 samples from hospitalized patients with COVID-19, their health care providers, and hospital surfaces before, during, and after admission. We screened for SARS-CoV-2 using RT-qPCR, characterized microbial communities using 16S rRNA gene amplicon sequencing, and used these bacterial profiles to classify SARS-CoV-2 RNA detection with a random forest model.ResultsSixteen percent of surfaces from COVID-19 patient rooms had detectable SARS-CoV-2 RNA, although infectivity was not assessed. The highest prevalence was in floor samples next to patient beds (39%) and directly outside their rooms (29%). Although bed rail samples more closely resembled the patient microbiome compared to floor samples, SARS-CoV-2 RNA was detected less often in bed rail samples (11%). SARS-CoV-2 positive samples had higher bacterial phylogenetic diversity in both human and surface samples and higher biomass in floor samples. 16S microbial community profiles enabled high classifier accuracy for SARS-CoV-2 status in not only nares, but also forehead, stool, and floor samples. Across these distinct microbial profiles, a single amplicon sequence variant from the genus Rothia strongly predicted SARS-CoV-2 presence across sample types, with greater prevalence in positive surface and human samples, even when compared to samples from patients in other intensive care units prior to the COVID-19 pandemic.ConclusionsThese results contextualize the vast diversity of microbial niches where SARS-CoV-2 RNA is detected and identify specific bacterial taxa that associate with the viral RNA prevalence both in the host and hospital environment. Video Abstract.

Microbiome data are sparse and high dimensional, so effective visualization of these data requires dimensionality reduction. To date, the most commonly used method for dimensionality reduction in the microbiome is calculation of between-sample microbial differences (beta diversity), followed by principal-coordinate analysis (PCoA). Uniform Manifold Approximation and Projection (UMAP) is an alternative method that can reduce the dimensionality of beta diversity distance matrices. Here, we demonstrate the benefits and limitations of using UMAP for dimensionality reduction on microbiome data. Using real data, we demonstrate that UMAP can improve the representation of clusters, especially when the clusters are composed of multiple subgroups. Additionally, we show that UMAP provides improved correlation of biological variation along a gradient with a reduced number of coordinates of the resulting embedding. Finally, we provide parameter recommendations that emphasize the preservation of global geometry. We therefore conclude that UMAP should be routinely used as a complementary visualization method for microbiome beta diversity studies. IMPORTANCE UMAP provides an additional method to visualize microbiome data. The method is extensible to any beta diversity metric used with PCoA, and our results demonstrate that UMAP can indeed improve visualization quality and correspondence with biological and technical variables of interest. The software to perform this analysis is available under an open-source license and can be obtained at https://github.com/knightlab-analyses/umap-microbiome-benchmarking; additionally, we have provided a QIIME 2 plugin for UMAP at https://github.com/biocore/q2-umap.

BackgroundEarly microbiota perturbations are associated with disorders that involve immunological underpinnings. Cesarean section (CS)-born babies show altered microbiota development in relation to babies born vaginally. Here we present the first statistically powered longitudinal study to determine the effect of restoring exposure to maternal vaginal fluids after CS birth.MethodsUsing 16S rRNA gene sequencing, we followed the microbial trajectories of multiple body sites in 177 babies over the first year of life; 98 were born vaginally, and 79 were born by CS, of whom 30 were swabbed with a maternal vaginal gauze right after birth.FindingsCompositional tensor factorization analysis confirmed that microbiota trajectories of exposed CS-born babies aligned more closely with that of vaginally born babies. Interestingly, the majority of amplicon sequence variants from maternal vaginal microbiomes on the day of birth were shared with other maternal sites, in contrast to non-pregnant women from the Human Microbiome Project (HMP) study.ConclusionsThe results of this observational study prompt urgent randomized clinical trials to test whether microbial restoration reduces the increased disease risk associated with CS birth and the underlying mechanisms. It also provides evidence of the pluripotential nature of maternal vaginal fluids to provide pioneer bacterial colonizers for the newborn body sites. This is the first study showing long-term naturalization of the microbiota of CS-born infants by restoring microbial exposure at birth.FundingC&D, Emch Fund, CIFAR, Chilean CONICYT and SOCHIPE, Norwegian Institute of Public Health, Emerald Foundation, NIH, National Institute of Justice, Janssen.

Obstructive sleep apnea (OSA), characterized by intermittent hypoxia and hypercapnia (IHC), affects the composition of the gut microbiome and metabolome. The gut microbiome has diurnal oscillations that play a crucial role in regulating circadian and overall metabolic homeostasis. Thus, we hypothesized that IHC adversely alters the gut luminal dynamics of key microbial families and metabolites. The objective of this study was to determine the diurnal dynamics of the fecal microbiome and metabolome of Apoe-/- mice after a week of IHC exposure. Individually housed, 10-week-old Apoe-/- mice on an atherogenic diet were split into two groups. One group was exposed to daily IHC conditions for 10 h (Zeitgeber time 2 [ZT2] to ZT12), while the other was maintained in room air. Six days after the initiation of the IHC conditions, fecal samples were collected every 4 h for 24 h (6 time points). We performed 16S rRNA gene amplicon sequencing and untargeted liquid chromatography-mass spectrometry (LC-MS) to assess changes in the microbiome and metabolome. IHC induced global changes in the cyclical dynamics of the gut microbiome and metabolome. Ruminococcaceae, Lachnospiraceae, S24-7, and Verrucomicrobiaceae had the greatest shifts in their diurnal oscillations. In the metabolome, bile acids, glycerolipids (phosphocholines and phosphoethanolamines), and acylcarnitines were greatly affected. Multi-omic analysis of these results demonstrated that Ruminococcaceae and tauro-β-muricholic acid (TβMCA) cooccur and are associated with IHC conditions and that Coriobacteriaceae and chenodeoxycholic acid (CDCA) cooccur and are associated with control conditions. IHC significantly change the diurnal dynamics of the fecal microbiome and metabolome, increasing members and metabolites that are proinflammatory and proatherogenic while decreasing protective ones. IMPORTANCE People with obstructive sleep apnea are at a higher risk of high blood pressure, type 2 diabetes, cardiac arrhythmias, stroke, and sudden cardiac death. We wanted to understand whether the gut microbiome changes induced by obstructive sleep apnea could potentially explain some of these medical problems. By collecting stool from a mouse model of this disease at multiple time points during the day, we studied how obstructive sleep apnea changed the day-night patterns of microbes and metabolites of the gut. Since the oscillations of the gut microbiome play a crucial role in regulating metabolism, changes in these oscillations can explain why these patients can develop so many metabolic problems. We found changes in microbial families and metabolites that regulate many metabolic pathways contributing to the increased risk for heart disease seen in patients with obstructive sleep apnea.

Standard workflows for analyzing microbiomes often include the creation and curation of phylogenetic trees. Here we present EMPress, an interactive web tool for visualizing trees in the context of microbiome, metabolome, and other community data scalable to trees with well over 500,000 nodes. EMPress provides novel functionality-including ordination integration and animations-alongside many standard tree visualization features and thus simplifies exploratory analyses of many forms of 'omic data.IMPORTANCE Phylogenetic trees are integral data structures for the analysis of microbial communities. Recent work has also shown the utility of trees constructed from certain metabolomic data sets, further highlighting their importance in microbiome research. The ever-growing scale of modern microbiome surveys has led to numerous challenges in visualizing these data. In this paper we used five diverse data sets to showcase the versatility and scalability of EMPress, an interactive web visualization tool. EMPress addresses the growing need for exploratory analysis tools that can accommodate large, complex multi-omic data sets.

One goal of microbial ecology researchers is to capture the maximum amount of information from all organisms in a sample. The recent COVID-19 pandemic, caused by the RNA virus SARS-CoV-2, has highlighted a gap in traditional DNA-based protocols, including the high-throughput methods the authors previously established as field standards. To enable simultaneous SARS-CoV-2 and microbial community profiling, the authors compared the relative performance of two total nucleic acid extraction protocols with the authors' previously benchmarked protocol. The authors included a diverse panel of environmental and host-associated sample types, including body sites commonly swabbed for COVID-19 testing. Here the authors present results comparing the cost, processing time, DNA and RNA yield, microbial community composition, limit of detection and well-to-well contamination between these protocols.

The translational power of human microbiome studies is limited by high interindividual variation. We describe a dimensionality reduction tool, compositional tensor factorization (CTF), that incorporates information from the same host across multiple samples to reveal patterns driving differences in microbial composition across phenotypes. CTF identifies robust patterns in sparse compositional datasets, allowing for the detection of microbial changes associated with specific phenotypes that are reproducible across datasets.

Obstructive sleep apnea (OSA), a common sleep disorder characterized by intermittent hypoxia and hypercapnia (IHC), increases atherosclerosis risk. However, the contribution of intermittent hypoxia (IH) or intermittent hypercapnia (IC) in promoting atherosclerosis remains unclear. Since gut microbiota and metabolites have been implicated in atherosclerosis, we examined whether IH or IC alters the microbiome and metabolome to induce a pro-atherosclerotic state. Apolipoprotein E deficient mice (ApoE-/- ), treated with IH or IC on a high-fat diet (HFD) for 10 weeks, were compared to Air controls. Atherosclerotic lesions were examined, gut microbiome was profiled using 16S rRNA gene amplicon sequencing and metabolome was assessed by untargeted mass spectrometry. In the aorta, IC-induced atherosclerosis was significantly greater than IH and Air controls (aorta, IC 11.1 ± 0.7% vs. IH 7.6 ± 0.4%, p < 0.05 vs. Air 8.1 ± 0.8%, p < 0.05). In the pulmonary artery (PA), however, IH, IC, and Air were significantly different from each other in atherosclerotic formation with the largest lesion observed under IH (PA, IH 40.9 ± 2.0% vs. IC 20.1 ± 2.6% vs. Air 12.2 ± 1.5%, p < 0.05). The most differentially abundant microbial families (p < 0.001) were Peptostreptococcaceae, Ruminococcaceae, and Erysipelotrichaceae. The most differentially abundant metabolites (p < 0.001) were tauro-β-muricholic acid, ursodeoxycholic acid, and lysophosphoethanolamine (18:0). We conclude that IH and IC (a) modulate atherosclerosis progression differently in distinct vascular beds with IC, unlike IH, facilitating atherosclerosis in both aorta and PA and (b) promote an atherosclerotic luminal gut environment that is more evident in IH than IC. We speculate that the resulting changes in the gut metabolome and microbiome interact differently with distinct vascular beds.

We show that SARS-CoV-2 spike protein interacts with both cellular heparan sulfate and angiotensin-converting enzyme 2 (ACE2) through its receptor-binding domain (RBD). Docking studies suggest a heparin/heparan sulfate-binding site adjacent to the ACE2-binding site. Both ACE2 and heparin can bind independently to spike protein in vitro, and a ternary complex can be generated using heparin as a scaffold. Electron micrographs of spike protein suggests that heparin enhances the open conformation of the RBD that binds ACE2. On cells, spike protein binding depends on both heparan sulfate and ACE2. Unfractionated heparin, non-anticoagulant heparin, heparin lyases, and lung heparan sulfate potently block spike protein binding and/or infection by pseudotyped virus and authentic SARS-CoV-2 virus. We suggest a model in which viral attachment and infection involves heparan sulfate-dependent enhancement of binding to ACE2. Manipulation of heparan sulfate or inhibition of viral adhesion by exogenous heparin presents new therapeutic opportunities.

In this cross-sectional study, we describe the composition and diversity of the gut microbiota among undernourished children living in urban slums of Mumbai, India, and determine how nutritional status, including anthropometric measurements, dietary intakes from complementary foods, feeding practices, and micronutrient concentrations, is associated with their gut microbiota. We collected rectal swabs from children aged 10 to 18 months living in urban slums of Mumbai participating in a randomized controlled feeding trial and conducted 16S rRNA sequencing to determine the composition of the gut microbiota. Across the study cohort, Proteobacteria dominated the gut microbiota at over 80% relative abundance, with Actinobacteria representation at

Many tools for dealing with compositional ' 'omics' data produce feature-wise values that can be ranked in order to describe features' associations with some sort of variation. These values include differentials (which describe features' associations with specified covariates) and feature loadings (which describe features' associations with variation along a given axis in a biplot). Although prior work has discussed the use of these 'rankings' as a starting point for exploring the log-ratios of particularly high- or low-ranked features, such exploratory analyses have previously been done using custom code to visualize feature rankings and the log-ratios of interest. This approach is laborious, prone to errors and raises questions about reproducibility. To address these problems we introduce Qurro, a tool that interactively visualizes a plot of feature rankings (a 'rank plot') alongside a plot of selected features' log-ratios within samples (a 'sample plot'). Qurro's interface includes various controls that allow users to select features from along the rank plot to compute a log-ratio; this action updates both the rank plot (through highlighting selected features) and the sample plot (through displaying the current log-ratios of samples). Here, we demonstrate how this unique interface helps users explore feature rankings and log-ratios simply and effectively.

QIIME 2 is a completely re-engineered microbiome bioinformatics platform based on the popular QIIME platform, which it has replaced. QIIME 2 facilitates comprehensive and fully reproducible microbiome data science, improving accessibility to diverse users by adding multiple user interfaces. QIIME 2 can be combined with Qiita, an open-source web-based platform, to re-use available data for meta-analysis. The following basic protocol describes how to install QIIME 2 on a single computer and analyze microbiome sequence data, from processing of raw DNA sequence reads through generating publishable interactive figures. These interactive figures allow readers of a study to interact with data with the same ease as its authors, advancing microbiome science transparency and reproducibility. We also show how plug-ins developed by the community to add analysis capabilities can be installed and used with QIIME 2, enhancing various aspects of microbiome analyses-e.g., improving taxonomic classification accuracy. Finally, we illustrate how users can perform meta-analyses combining different datasets using readily available public data through Qiita. In this tutorial, we analyze a subset of the Early Childhood Antibiotics and the Microbiome (ECAM) study, which tracked the microbiome composition and development of 43 infants in the United States from birth to 2 years of age, identifying microbiome associations with antibiotic exposure, delivery mode, and diet. For more information about QIIME 2, see https://qiime2.org. To troubleshoot or ask questions about QIIME 2 and microbiome analysis, join the active community at https://forum.qiime2.org. © 2020 The Authors. Basic Protocol: Using QIIME 2 with microbiome data Support Protocol: Further microbiome analyses.

Animal work indicates exposure to air pollutants may alter the composition of the gut microbiota. This study examined relationships between air pollutants and the gut microbiome in young adults residing in Southern California. Our results demonstrate significant associations between exposure to air pollutants and the composition of the gut microbiome using whole-genome sequencing. Higher exposure to 24-hour O3 was associated with lower Shannon diversity index, higher Bacteroides caecimuris, and multiple gene pathways, including L-ornithine de novo biosynthesis as well as pantothenate and coenzyme A biosynthesis I. Among other pollutants, higher NO2 exposure was associated with fewer taxa, including higher Firmicutes. The percent variation in gut bacterial composition that was explained by air pollution exposure was up to 11.2% for O3 concentrations, which is large compared to the effect size for many other covariates reported in healthy populations. This study provides the first evidence of significant associations between exposure to air pollutants and the compositional and functional profile of the human gut microbiome. These results identify O3 as an important pollutant that may alter the human gut microbiome.

Lifestyle factors, such as diet, strongly influence the structure, diversity, and composition of the microbiome. While we have witnessed over the last several years a resurgence of interest in fermented foods, no study has specifically explored the effects of their consumption on gut microbiota in large cohorts. To assess whether the consumption of fermented foods is associated with a systematic signal in the gut microbiome and metabolome, we used a multi-omic approach (16S rRNA amplicon sequencing, metagenomic sequencing, and untargeted mass spectrometry) to analyze stool samples from 6,811 individuals from the American Gut Project, including 115 individuals specifically recruited for their frequency of fermented food consumption for a targeted 4-week longitudinal study. We observed subtle but statistically significant differences between consumers and nonconsumers in beta diversity as well as differential taxa between the two groups. We found that the metabolome of fermented food consumers was enriched with conjugated linoleic acid (CLA), a putatively health-promoting molecule. Cross-omic analyses between metagenomic sequencing and mass spectrometry suggest that CLA may be driven by taxa associated with fermented food consumers. Collectively, we found modest yet persistent signatures associated with fermented food consumption that appear present in multiple -omic types which motivate further investigation of how different types of fermented food impact the gut microbiome and overall health.IMPORTANCE Public interest in the effects of fermented food on the human gut microbiome is high, but limited studies have explored the association between fermented food consumption and the gut microbiome in large cohorts. Here, we used a combination of omics-based analyses to study the relationship between the microbiome and fermented food consumption in thousands of people using both cross-sectional and longitudinal data. We found that fermented food consumers have subtle differences in their gut microbiota structure, which is enriched in conjugated linoleic acid, thought to be beneficial. The results suggest that further studies of specific kinds of fermented food and their impacts on the microbiome and health will be useful.

As metagenomic studies move to increasing numbers of samples, communities like the human gut may benefit more from the assembly of abundant microbes in many samples, rather than the exhaustive assembly of fewer samples. We term this approach leaderboard metagenome sequencing. To explore protocol optimization for leaderboard metagenomics in real samples, we introduce a benchmark of library prep and sequencing using internal references generated by synthetic long-read technology, allowing us to evaluate high-throughput library preparation methods against gold-standard reference genomes derived from the samples themselves. We introduce a low-cost protocol for high-throughput library preparation and sequencing.

Dietary habits have been associated with alterations of the human gut resident microorganisms contributing to obesity, diabetes and cancer1. In Western diets, red meat is a frequently eaten food2, but long-term consumption has been associated with increased risk of disease3,4. Red meat is enriched in N-glycolylneuraminic acid (Neu5Gc) that cannot be synthesized by humans5. However, consumption can cause Neu5Gc incorporation into cell surface glycans6, especially in carcinomas4,7. As a consequence, an inflammatory response is triggered when Neu5Gc-containing glycans encounter circulating anti-Neu5Gc antibodies8,9. Although bacteria can use free sialic acids as a nutrient source10-12, it is currently unknown if gut microorganisms contribute to releasing Neu5Gc from food. We found that a Neu5Gc-rich diet induces changes in the gut microbiota, with Bacteroidales and Clostridiales responding the most. Genome assembling of mouse and human shotgun metagenomic sequencing identified bacterial sialidases with previously unobserved substrate preference for Neu5Gc-containing glycans. X-ray crystallography revealed key amino acids potentially contributing to substrate preference. Additionally, we verified that mouse and human sialidases were able to release Neu5Gc from red meat. The release of Neu5Gc from red meat using bacterial sialidases could reduce the risk of inflammatory diseases associated with red meat consumption, including colorectal cancer4 and atherosclerosis13.

Impact statementConsidering the clear effects of microbiota on important aspects of animal biology and development (including in humans), this topic is timely and broadly appealing, as it compels us to consider the possibilities of altering the microbiome (without antibiotics) to positively affect animal health. In this review, we highlight three general approaches to manipulating the microbiome that have demonstrated success and promise for use in animal health. We also point out knowledge gaps where further inquiry would most benefit the field. Our paper not only provides a short and digestible overview of the current state of application, but also calls for further exploration of the microbial diversity at hand to expand our toolkit, while also leveraging the diversity and flexibility of animal systems to better understand mechanisms of efficacy.

The human oral cavity and upper airway serves as an early barrier and reservoir in the transmission of SARS-CoV-2. Saliva in this microenvironment may serve as a key host factor that can modulate susceptibility to infection and eventual infection of the lower respiratory tract. We sought to analyze the content and composition of heparan sulfate, a glycosaminoglycan identified as an important co-receptor for viral entry, and whether there is any correlation with SARS-CoV-2 infection. We enlisted 98 participants stratified by age, gender, race, and COVID-19 history. Notably, the concentration of heparan sulfate in saliva increased with age, and its composition showed a wide range of variability within each age group independently of age. Heparan sulfate concentration and composition did not differ significantly with gender, ethnicity or race. Compared to patients with no COVID-19 history, patients with previous infection had a similar salivary heparan sulfate concentration, but significant increases in overall sulfation were noted. Moreover, in a subset of participants, for which data was available pre- and post- infection, significant elevation in N-sulfoglucosamine in heparan sulfate was observed post- COVID-19. Examination of salivary bacterial 16S rRNA, showed a significant reduction in species predicted to possess heparan sulfate-modifying capacity among participants >60 years old, which correlates with the increase in heparan sulfate content in older individuals. These findings demonstrate a surprisingly wide variation in heparan sulfate content and composition in saliva across the sampled population and confirm other findings showing variation in content and composition of glycosaminoglycans in blood and urine. [ABSTRACT FROM AUTHOR]

The central aims of many host or environmental microbiome studies are to elucidate factors associated with microbial community compositions and to relate microbial features to outcomes. However, these aims are often complicated by difficulties stemming from high-dimensionality, non-normality, sparsity, and the compositional nature of microbiome data sets. A key tool in microbiome analysis is beta diversity, defined by the distances between microbial samples. Many different distance metrics have been proposed, all with varying discriminatory power on data with differing characteristics. Here, we propose a compositional beta diversity metric rooted in a centered log-ratio transformation and matrix completion called robust Aitchison PCA. We demonstrate the benefits of compositional transformations upstream of beta diversity calculations through simulations. Additionally, we demonstrate improved effect size, classification accuracy, and robustness to sequencing depth over the current methods on several decreased sample subsets of real microbiome data sets. Finally, we highlight the ability of this new beta diversity metric to retain the feature loadings linked to sample ordinations revealing salient intercommunity niche feature importance. IMPORTANCE By accounting for the sparse compositional nature of microbiome data sets, robust Aitchison PCA can yield high discriminatory power and salient feature ranking between microbial niches. The software to perform this analysis is available under an open-source license and can be obtained at https://github.com/biocore/DEICODE; additionally, a QIIME 2 plugin is provided to perform this analysis at https://library.qiime2.org/plugins/deicode/.

Although much work has linked the human microbiome to specific phenotypes and lifestyle variables, data from different projects have been challenging to integrate and the extent of microbial and molecular diversity in human stool remains unknown. Using standardized protocols from the Earth Microbiome Project and sample contributions from over 10,000 citizen-scientists, together with an open research network, we compare human microbiome specimens primarily from the United States, United Kingdom, and Australia to one another and to environmental samples. Our results show an unexpected range of beta-diversity in human stool microbiomes compared to environmental samples; demonstrate the utility of procedures for removing the effects of overgrowth during room-temperature shipping for revealing phenotype correlations; uncover new molecules and kinds of molecular communities in the human stool metabolome; and examine emergent associations among the microbiome, metabolome, and the diversity of plants that are consumed (rather than relying on reductive categorical variables such as veganism, which have little or no explanatory power). We also demonstrate the utility of the living data resource and cross-cohort comparison to confirm existing associations between the microbiome and psychiatric illness and to reveal the extent of microbiome change within one individual during surgery, providing a paradigm for open microbiome research and education. IMPORTANCE We show that a citizen science, self-selected cohort shipping samples through the mail at room temperature recaptures many known microbiome results from clinically collected cohorts and reveals new ones. Of particular interest is integrating n = 1 study data with the population data, showing that the extent of microbiome change after events such as surgery can exceed differences between distinct environmental biomes, and the effect of diverse plants in the diet, which we confirm with untargeted metabolomics on hundreds of samples.