Peripheral neuronal activation shapes the microbiome and alters gut physiology.

The gastrointestinal (GI) tract is innervated by intrinsic neurons of the enteric nervous system (ENS) and extrinsic neurons of the central nervous system and peripheral ganglia. The GI tract also harbors a diverse microbiome, but interactions between the ENS and the microbiome remain poorly understood. Here, we activate choline acetyltransferase (ChAT)-expressing or tyrosine hydroxylase (TH)-expressing gut-associated neurons in mice to determine effects on intestinal microbial communities and their metabolites as well as on host physiology. The resulting multi-omics datasets support broad roles for discrete peripheral neuronal subtypes in shaping microbiome structure, including modulating bile acid profiles and fungal colonization. Physiologically, activation of either ChAT+ or TH+ neurons increases fecal output, while only ChAT+ activation results in increased colonic contractility and diarrhea-like fluid secretion. These findings suggest that specific subsets of peripherally activated neurons differentially regulate the gut microbiome and GI physiology in mice without involvement of signals from the brain.

Comprehensive evaluation of shotgun metagenomics, amplicon sequencing, and harmonization of these platforms for epidemiological studies.

In a large cohort of 1,772 participants from the Hispanic Community Health Study/Study of Latinos with overlapping 16SV4 rRNA gene (bacterial amplicon), ITS1 (fungal amplicon), and shotgun sequencing data, we demonstrate that 16SV4 amplicon sequencing and shotgun metagenomics offer the same level of taxonomic accuracy for bacteria at the genus level even at shallow sequencing depths. In contrast, for fungal taxa, we did not observe meaningful agreements between shotgun and ITS1 amplicon results. Finally, we show that amplicon and shotgun data can be harmonized and pooled to yield larger microbiome datasets with excellent agreement (<1% effect size variance across three independent outcomes) using pooled amplicon/shotgun data compared to pure shotgun metagenomic analysis. Thus, there are multiple approaches to study the microbiome in epidemiological studies, and we provide a demonstration of a powerful pooling approach that will allow researchers to leverage the massive amount of amplicon sequencing data generated over the last two decades.

Standardized multi-omics of Earth's microbiomes reveals microbial and metabolite diversity.

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.

Pan-cancer analyses reveal cancer-type-specific fungal ecologies and bacteriome interactions.

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.

Postnatal exposure to ambient air pollutants is associated with the composition of the infant gut microbiota at 6-months of age.

Epidemiological studies in adults have shown that exposure to ambient air pollution (AAP) is associated with the composition of the adult gut microbiome, but these relationships have not been examined in infancy. We aimed to determine if 6-month postnatal AAP exposure was associated with the infant gut microbiota at 6 months of age in a cohort of Latino mother-infant dyads from the Southern California Mother's Milk Study (n = 103). We estimated particulate matter (PM2.5 and PM10) and nitrogen dioxide (NO2) exposure from birth to 6-months based on residential address histories. We characterized the infant gut microbiota using 16S rRNA amplicon sequencing at 6-months of age. At 6-months, the gut microbiota was dominated by the phyla Bacteroidetes, Firmicutes, Proteobacteria, and Actinobacteria. Our results show that, after adjusting for important confounders, postnatal AAP exposure was associated with the composition of the gut microbiota. As an example, PM10 exposure was positively associated with Dialister, Dorea, Acinetobacter, and Campylobacter while PM2.5 was positively associated with Actinomyces. Further, exposure to PM10 and PM2.5 was inversely associated with Alistipes and NO2 exposure was positively associated with Actinomyces, Enterococcus, Clostridium, and Eubacterium. Several of these taxa have previously been linked with systemic inflammation, including the genera Dialister and Dorea. This study provides the first evidence of significant associations between exposure to AAP and the composition of the infant gut microbiota, which may have important implications for future infant health and development.

Enhancing untargeted metabolomics using metadata-based source annotation.

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.

A comparison of six DNA extraction protocols for 16S, ITS and shotgun metagenomic sequencing of microbial communities.

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.

Phylogeny-Aware Analysis of Metagenome Community Ecology Based on Matched Reference Genomes while Bypassing Taxonomy.

We introduce the operational genomic unit (OGU) method, a metagenome analysis strategy that directly exploits sequence alignment hits to individual reference genomes as the minimum unit for assessing the diversity of microbial communities and their relevance to environmental factors. This approach is independent of taxonomic classification, granting the possibility of maximal resolution of community composition, and organizes features into an accurate hierarchy using a phylogenomic tree. The outputs are suitable for contemporary analytical protocols for community ecology, differential abundance, and supervised learning while supporting phylogenetic methods, such as UniFrac and phylofactorization, that are seldom applied to shotgun metagenomics despite being prevalent in 16S rRNA gene amplicon studies. As demonstrated in two real-world case studies, the OGU method produces biologically meaningful patterns from microbiome data sets. Such patterns further remain detectable at very low metagenomic sequencing depths. Compared with taxonomic unit-based analyses implemented in currently adopted metagenomics tools, and the analysis of 16S rRNA gene amplicon sequence variants, this method shows superiority in informing biologically relevant insights, including stronger correlation with body environment and host sex on the Human Microbiome Project data set and more accurate prediction of human age by the gut microbiomes of Finnish individuals included in the FINRISK 2002 cohort. We provide Woltka, a bioinformatics tool to implement this method, with full integration with the QIIME 2 package and the Qiita web platform, to facilitate adoption of the OGU method in future metagenomics studies. IMPORTANCE Shotgun metagenomics is a powerful, yet computationally challenging, technique compared to 16S rRNA gene amplicon sequencing for decoding the composition and structure of microbial communities. Current analyses of metagenomic data are primarily based on taxonomic classification, which is limited in feature resolution. To solve these challenges, we introduce operational genomic units (OGUs), which are the individual reference genomes derived from sequence alignment results, without further assigning them taxonomy. The OGU method advances current read-based metagenomics in two dimensions: (i) providing maximal resolution of community composition and (ii) permitting use of phylogeny-aware tools. Our analysis of real-world data sets shows that it is advantageous over currently adopted metagenomic analysis methods and the finest-grained 16S rRNA analysis methods in predicting biological traits. We thus propose the adoption of OGUs as an effective practice in metagenomic studies.

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.

Transcriptional Profiles of a Foliar Fungal Endophyte (Pestalotiopsis, Ascomycota) and Its Bacterial Symbiont (Luteibacter, Gammaproteobacteria) Reveal Sulfur Exchange and Growth Regulation during Early Phases of Symbiotic Interaction.

Symbiosis with bacteria is widespread among eukaryotes, including fungi. Bacteria that live within fungal mycelia (endohyphal bacteria) occur in many plant-associated fungi, including diverse Mucoromycota and Dikarya. Pestalotiopsis sp. strain 9143 is a filamentous ascomycete isolated originally as a foliar endophyte of Platycladus orientalis (Cupressaceae). It is infected naturally with the endohyphal bacterium Luteibacter sp. strain 9143, which influences auxin and enzyme production by its fungal host. Previous studies have used transcriptomics to examine similar symbioses between endohyphal bacteria and root-associated fungi such as arbuscular mycorrhizal fungi and plant pathogens. However, currently there are no gene expression studies of endohyphal bacteria of Ascomycota, the most species-rich fungal phylum. To begin to understand such symbioses, we developed methods for assessing gene expression by Pestalotiopsis sp. and Luteibacter sp. when grown in coculture and when each was grown axenically. Our assays showed that the density of Luteibacter sp. in coculture was greater than in axenic culture, but the opposite was true for Pestalotiopsis sp. Dual-transcriptome sequencing (RNA-seq) data demonstrate that growing in coculture modulates developmental and metabolic processes in both the fungus and bacterium, potentially through changes in the balance of organic sulfur via methionine acquisition. Our analyses also suggest an unexpected, potential role of the bacterial type VI secretion system in symbiosis establishment, expanding current understanding of the scope and dynamics of fungal-bacterial symbioses. IMPORTANCE Interactions between microbes and their hosts have important outcomes for host and environmental health. Foliar fungal endophytes that infect healthy plants can harbor facultative endosymbionts called endohyphal bacteria, which can influence the outcome of plant-fungus interactions. These bacterial-fungal interactions can be influential but are poorly understood, particularly from a transcriptome perspective. Here, we report on a comparative, dual-RNA-seq study examining the gene expression patterns of a foliar fungal endophyte and a facultative endohyphal bacterium when cultured together versus separately. Our findings support a role for the fungus in providing organic sulfur to the bacterium, potentially through methionine acquisition, and the potential involvement of a bacterial type VI secretion system in symbiosis establishment. This work adds to the growing body of literature characterizing endohyphal bacterial-fungal interactions, with a focus on a model facultative bacterial-fungal symbiosis in two species-rich lineages, the Ascomycota and Proteobacteria.

Swapping Metagenomics Preprocessing Pipeline Components Offers Speed and Sensitivity Increases.

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.

Early life gut microbiota is associated with rapid infant growth in Hispanics from Southern California.

We aimed to determine if the newborn gut microbiota is an underlying determinant of early life growth trajectories. 132 Hispanic infants were recruited at 1-month postpartum. The infant gut microbiome was characterized using 16S rRNA amplicon sequencing. Rapid infant growth was defined as a weight-for-age z-score (WAZ) change greater than 0.67 between birth and 12-months of age. Measures of infant growth included change in WAZ, weight-for-length z-score (WLZ), and body mass index (BMI) z-scores from birth to 12-months and infant anthropometrics at 12-months (weight, skinfold thickness). Of the 132 infants, 40% had rapid growth in the first year of life. Multiple metrics of alpha-diversity predicted rapid infant growth, including a higher Shannon diversity (OR = 1.83; 95% CI: 1.07-3.29; p = .03), Faith's phylogenic diversity (OR = 1.41, 95% CI: 1.05-1.94; p = .03), and richness (OR = 1.04, 95% CI: 1.01-1.08; p = .02). Many of these alpha-diversity metrics were also positively associated with increases in WAZ, WLZ, and BMI z-scores from birth to 12-months (pall<0.05). Importantly, we identified subsets of microbial consortia whose abundance were correlated with these same measures of infant growth. We also found that rapid growers were enriched in multiple taxa belonging to genera such as Acinetobacter, Collinsella, Enterococcus, Neisseria, and Parabacteroides. Moreover, measures of the newborn gut microbiota explained up to an additional 5% of the variance in rapid growth beyond known clinical predictors (R2 = 0.37 vs. 0.32, p < .01). These findings indicate that a more mature gut microbiota, characterized by increased alpha-diversity, at as early as 1-month of age, may influence infant growth trajectories in the first year of life.

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.

Microbiome Metadata Standards: Report of the National Microbiome Data Collaborative's Workshop and Follow-On Activities.

Microbiome samples are inherently defined by the environment in which they are found. Therefore, data that provide context and enable interpretation of measurements produced from biological samples, often referred to as metadata, are critical. Important contributions have been made in the development of community-driven metadata standards; however, these standards have not been uniformly embraced by the microbiome research community. To understand how these standards are being adopted, or the barriers to adoption, across research domains, institutions, and funding agencies, the National Microbiome Data Collaborative (NMDC) hosted a workshop in October 2019. This report provides a summary of discussions that took place throughout the workshop, as well as outcomes of the working groups initiated at the workshop.

A comparison of DNA/RNA extraction protocols for high-throughput sequencing of microbial communities.

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.

Feasibility of using alternative swabs and storage solutions for paired SARS-CoV-2 detection and microbiome analysis in the hospital environment.

BACKGROUND: Determining the role of fomites in the transmission of SARS-CoV-2 is essential in the hospital setting and will likely be important outside of medical facilities as governments around the world make plans to ease COVID-19 public health restrictions and attempt to safely reopen economies. Expanding COVID-19 testing to include environmental surfaces would ideally be performed with inexpensive swabs that could be transported safely without concern of being a source of new infections. However, CDC-approved clinical-grade sampling supplies and techniques using a synthetic swab are expensive, potentially expose laboratory workers to viable virus and prohibit analysis of the microbiome due to the presence of antibiotics in viral transport media (VTM). To this end, we performed a series of experiments comparing the diagnostic yield using five consumer-grade swabs (including plastic and wood shafts and various head materials including cotton, synthetic, and foam) and one clinical-grade swab for inhibition to RNA. For three of these swabs, we evaluated performance to detect SARS-CoV-2 in twenty intensive care unit (ICU) hospital rooms of patients including COVID-19+ patients. All swabs were placed in 95% ethanol and further evaluated in terms of RNase activity. SARS-CoV-2 was measured both directly from the swab and from the swab eluent. RESULTS: Compared to samples collected in VTM, 95% ethanol demonstrated significant inhibition properties against RNases. When extracting directly from the swab head as opposed to the eluent, RNA recovery was approximately 2-4× higher from all six swab types tested as compared to the clinical standard of testing the eluent from a CDC-approved synthetic (SYN) swab. The limit of detection (LoD) of SARS-CoV-2 from floor samples collected using the consumer-grade plastic (CGp) or research-grade plastic The Microsetta Initiative (TMI) swabs was similar or better than the SYN swab, further suggesting that swab type does not impact RNA recovery as measured by the abundance of SARS-CoV-2. The LoD for TMI was between 0 and 362.5 viral particles, while SYN and CGp were both between 725 and 1450 particles. Lastly microbiome analyses (16S rRNA gene sequencing) of paired samples (nasal and floor from same patient room) collected using different swab types in triplicate indicated that microbial communities were not impacted by swab type, but instead driven by the patient and sample type. CONCLUSIONS: Compared to using a clinical-grade synthetic swab, detection of SARS-CoV-2 from environmental samples collected from ICU rooms of patients with COVID was similar using consumer-grade swabs, stored in 95% ethanol. The yield was best from the swab head rather than the eluent and the low level of RNase activity and lack of antibiotics in these samples makes it possible to perform concomitant microbiome analyses. Video abstract.

A comparison of DNA/RNA extraction protocols for high-throughput sequencing of microbial communities.

One goal among 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 we previously established as field standards. To enable simultaneous SARS-CoV-2 and microbial community profiling, we compare the relative performance of two total nucleic acid extraction protocols and our previously benchmarked protocol. We included a diverse panel of environmental and host-associated sample types, including body sites commonly swabbed for COVID-19 testing. Here we 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. Accession numbers: Raw sequence data were deposited at the European Nucleotide Archive (accession#: ERP124610) and raw and processed data are available at Qiita (Study ID: 12201). All processing and analysis code is available on GitHub ( github.com/justinshaffer/Extraction_test_MagMAX ). Methods summary: To allow for downstream applications involving RNA-based organisms such as SARS-CoV-2, we compared the two extraction protocols designed to extract DNA and RNA against our previously established protocol for extracting only DNA for microbial community analyses. Across 10 diverse sample types, one of the two protocols was equivalent or better than our established DNA-based protocol. Our conclusion is based on per-sample comparisons of DNA and RNA yield, the number of quality sequences generated, microbial community alpha- and beta-diversity and taxonomic composition, the limit of detection, and extent of well-to-well contamination.

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.