A conserved interdomain microbial network underpins cadaver decomposition despite environmental variables.

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.

Integrated Genomic and Social Network Analyses of SARS-CoV-2 Transmission in the Healthcare Setting.

BACKGROUND: Infection prevention (IP) measures are designed to mitigate the transmission of pathogens in healthcare. Using large-scale viral genomic and social network analyses, we determined if IP measures used during the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic were adequate in protecting healthcare workers (HCWs) and patients from acquiring SARS-CoV-2. METHODS: We performed retrospective cross-sectional analyses of viral genomics from all available SARS-CoV-2 viral samples collected at UC San Diego Health and social network analysis using the electronic medical record to derive temporospatial overlap of infections among related viromes and supplemented with contact tracing data. The outcome measure was any instance of healthcare transmission, defined as cases with closely related viral genomes and epidemiological connection within the healthcare setting during the infection window. Between November 2020 through January 2022, 12 933 viral genomes were obtained from 35 666 patients and HCWs. RESULTS: Among 5112 SARS-CoV-2 viral samples sequenced from the second and third waves of SARS-CoV-2 (pre-Omicron), 291 pairs were derived from persons with a plausible healthcare overlap. Of these, 34 pairs (12%) were phylogenetically linked: 19 attributable to household and 14 to healthcare transmission. During the Omicron wave, 2106 contact pairs among 7821 sequences resulted in 120 (6%) related pairs among 32 clusters, of which 10 were consistent with healthcare transmission. Transmission was more likely to occur in shared spaces in the older hospital compared with the newer hospital (2.54 vs 0.63 transmission events per 1000 admissions, P < .001). CONCLUSIONS: IP strategies were effective at identifying and preventing healthcare SARS-CoV-2 transmission.

Effects of variation in sample storage conditions and swab order on 16S vaginal microbiome analyses.

IMPORTANCE: The composition of the human vaginal microbiome has been linked to a variety of medical conditions including yeast infection, bacterial vaginosis, and sexually transmitted infection. The vaginal microbiome is becoming increasingly acknowledged as a key factor in personal health, and it is essential to establish methods to collect and process accurate samples with self-collection techniques to allow large, population-based studies. In this study, we investigate if using AssayAssure Genelock, a nucleic acid preservative, introduces microbial biases in self-collected vaginal samples. To our knowledge, we also contribute some of the first evidence regarding the impacts of multiple swabs taken at one time point. Vaginal samples have relatively low biomass, so the ability to collect multiple swabs from a unique participant at a single time would greatly improve the replicability and data available for future studies. This will hopefully lay the groundwork to gain a more complete and accurate understanding of the vaginal microbiome.

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.

Effects of Variation in Urine Sample Storage Conditions on 16S Urogenital Microbiome Analyses.

Replicability is a well-established challenge in microbiome research with a variety of contributing factors at all stages, from sample collection to code execution. Here, we focus on voided urine sample storage conditions for urogenital microbiome analysis. Using urine samples collected from 10 adult females, we investigated the microbiome preservation efficacy of AssayAssure Genelock (Genelock), compared with no preservative, under different temperature conditions. We varied temperature over 48 h in order to examine the impact of conditions samples may experience with home voided urine collection and shipping to a central biorepository. The following common lab and shipping conditions were investigated: -20°C, ambient temperature, 4°C, freeze-thaw cycle, and heat cycle. At 48 h, all samples were stored at -80°C until processing. After generating 16S rRNA gene amplicon sequencing data using the highly sensitive KatharoSeq protocol, we observed individual variation in both alpha and beta diversity metrics below interhuman differences, corroborating reports of individual microbiome variability in other specimen types. While there was no significant difference in beta diversity when comparing Genelock versus no preservative, we did observe a higher concordance with Genelock samples shipped at colder temperatures (-20°C and 4°C) when compared with the samples shipped at -20°C without preservative. Our results indicate that Genelock does not introduce a significant amount of microbial bias when used on a range of temperatures and is most effective at colder temperatures. IMPORTANCE The urogenital microbiome is an understudied yet important human microbiome niche. Research has been stimulated by the relatively recent discovery that urine is not sterile; urinary tract microbes have been linked to health problems, including urinary infections, incontinence, and cancer. The quality of life and economic impact of UTIs and urgency incontinence alone are enormous, with $3.5 billion and $82.6 billion, respectively, spent in the United States. annually. Given the low biomass of urine, novelty of the field, and limited reproducibility evidence, it is critical to study urine sample storage conditions to optimize scientific rigor. Efficient and reliable preservation methods inform methods for home self-sample collection and shipping, increasing the potential use in larger-scale studies. Here, we examined both buffer and temperature variation effects on 16S rRNA gene amplicon sequencing results from urogenital samples, providing data on the consequences of common storage methods on urogenital microbiome results.

Compositionally Aware Phylogenetic Beta-Diversity Measures Better Resolve Microbiomes Associated with Phenotype.

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.

Efficient computation of Faith's phylogenetic diversity with applications in characterizing microbiomes.

The number of publicly available microbiome samples is continually growing. As data set size increases, bottlenecks arise in standard analytical pipelines. Faith's phylogenetic diversity (Faith's PD) is a highly utilized phylogenetic alpha diversity metric that has thus far failed to effectively scale to trees with millions of vertices. Stacked Faith's phylogenetic diversity (SFPhD) enables calculation of this widely adopted diversity metric at a much larger scale by implementing a computationally efficient algorithm. The algorithm reduces the amount of computational resources required, resulting in more accessible software with a reduced carbon footprint, as compared to previous approaches. The new algorithm produces identical results to the previous method. We further demonstrate that the phylogenetic aspect of Faith's PD provides increased power in detecting diversity differences between younger and older populations in the FINRISK study's metagenomic data.

SARS-CoV-2 detection status associates with bacterial community composition in patients and the hospital environment.

BACKGROUND: SARS-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. METHODS: We 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. RESULTS: Sixteen 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. CONCLUSIONS: These 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.

EMPress Enables Tree-Guided, Interactive, and Exploratory Analyses of Multi-omic Data Sets.

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.

Microbial context predicts SARS-CoV-2 prevalence in patients and the hospital built environment.

Synergistic effects of bacteria on viral stability and transmission are widely documented but remain unclear in the context of SARS-CoV-2. We collected 972 samples from hospitalized ICU patients with coronavirus disease 2019 (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 contextualized the massive microbial diversity in this dataset in a meta-analysis of over 20,000 samples. Sixteen percent of surfaces from COVID-19 patient rooms were positive, with the highest prevalence in floor samples next to patient beds (39%) and directly outside their rooms (29%). Although bed rail samples increasingly resembled the patient microbiome throughout their stay, SARS-CoV-2 was less frequently detected there (11%). Despite surface contamination in almost all patient rooms, no health care workers providing COVID-19 patient care contracted the disease. SARS-CoV-2 positive samples had higher bacterial phylogenetic diversity across human and surface samples, and higher biomass in floor samples. 16S microbial community profiles allowed for 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 was highly predictive of SARS-CoV-2 across sample types, and had higher prevalence in positive surface and human samples, even when comparing to samples from patients in another intensive care unit prior to the COVID-19 pandemic. These results suggest that bacterial communities contribute to viral prevalence both in the host and hospital environment.

Phylogenomics of 10,575 genomes reveals evolutionary proximity between domains Bacteria and Archaea.

Rapid growth of genome data provides opportunities for updating microbial evolutionary relationships, but this is challenged by the discordant evolution of individual genes. Here we build a reference phylogeny of 10,575 evenly-sampled bacterial and archaeal genomes, based on a comprehensive set of 381 markers, using multiple strategies. Our trees indicate remarkably closer evolutionary proximity between Archaea and Bacteria than previous estimates that were limited to fewer "core" genes, such as the ribosomal proteins. The robustness of the results was tested with respect to several variables, including taxon and site sampling, amino acid substitution heterogeneity and saturation, non-vertical evolution, and the impact of exclusion of candidate phyla radiation (CPR) taxa. Our results provide an updated view of domain-level relationships.