Although disinfection is key to infection control, the colonization patterns and resistomes of hospital-environment microbes remain underexplored. We report the first extensive genomic characterization of microbiomes, pathogens and antibiotic resistance cassettes in a tertiary-care hospital, from repeated sampling (up to 1.5 years apart) of 179 sites associated with 45 beds. Deep shotgun metagenomics unveiled distinct ecological niches of microbes and antibiotic resistance genes characterized by biofilm-forming and human-microbiome-influenced environments with corresponding patterns of spatiotemporal divergence. Quasi-metagenomics with nanopore sequencing provided thousands of high-contiguity genomes, phage and plasmid sequences (>60% novel), enabling characterization of resistome and mobilome diversity and dynamic architectures in hospital environments. Phylogenetics identified multidrug-resistant strains as being widely distributed and stably colonizing across sites. Comparisons with clinical isolates indicated that such microbes can persist in hospitals for extended periods (>8 years), to opportunistically infect patients. These findings highlight the importance of characterizing antibiotic resistance reservoirs in hospitals and establish the feasibility of systematic surveys to target resources for preventing infections.

Background The dynamics of microbial communities is driven by a range of interactions from symbiosis to predator-prey relationships, the majority of which are poorly understood. With the increasing availability of high-throughput microbiome taxonomic profiling data, it is now conceivable to directly learn the ecological models that explicitly define microbial interactions and explain community dynamics. The applicability of these approaches is severely limited by the lack of accurate absolute cell density measurements (biomass). Methods We present a new computational approach that resolves this key limitation in the inference of generalized Lotka-Volterra models (gLVMs) by coupling biomass estimation and model inference with an expectation-maximization algorithm (BEEM). Results BEEM outperforms the state-of-the-art methods for inferring gLVMs, while simultaneously eliminating the need for additional experimental biomass data as input. BEEM’s application to previously inaccessible public datasets (due to the lack of biomass data) allowed us to construct ecological models of microbial communities in the human gut on a per-individual basis, revealing personalized dynamics and keystone species. Conclusions BEEM addresses a key bottleneck in “systems analysis” of microbiomes by enabling accurate inference of ecological models from high throughput sequencing data without the need for experimental biomass measurements. Electronic supplementary material The online version of this article (10.1186/s40168-019-0729-z) contains supplementary material, which is available to authorized users.

There is growing attention surrounding hospital acquired infections (HAIs) due to high associated healthcare costs, compounded by the scourge of widespread multi-antibiotic resistance. Although hospital environment disinfection is well acknowledged to be key for infection control, an understanding of colonization patterns and resistome profiles of environment-dwelling microbes is currently lacking. We report the first extensive genomic characterization of microbiomes (355), common HAI-associated microbes (891) and transmissible drug resistance cassettes (1435) in a tertiary hospital environment based on a 2-timepoint sampling of 179 sites from 45 beds. Deep shotgun metagenomic sequencing unveiled two distinct ecological niches of microbes and antibiotic resistance genes characterized by biofilm-forming and human microbiome influenced environments that display corresponding patterns of divergence over space and time. To study common nosocomial pathogens that were typically present at low abundances, a combination of culture enrichment and long-read nanopore sequencing was used to obtain thousands of high contiguity genomes (2347) and closed plasmids (5910), a significant fraction of which (>58%) are not represented in current sequence databases. These high-quality assemblies and metadata enabled a rich characterization of resistance gene combinations, plasmid architectures, and the dynamic nature of hospital environment resistomes and their reservoirs. Phylogenetic analysis identified multidrug resistant clonal strains as being more widely disseminated and stably colonizing across hospital sites. Further genomic comparisons with clinical isolates across multiple species supports the hypothesis that multidrug resistant strains can persist in the hospital environment for extended periods (>8 years) to opportunistically infect patients. These findings highlight the importance of characterizing antibiotic resistance reservoirs in the hospital environment and establishes the feasibility of systematic genomic surveys to help target resources more efficiently for preventing HAIs.

The structure and function of diverse microbial communities is underpinned by ecological interactions that remain uncharacterized. With rapid adoption of next-generation sequencing for studying microbiomes, data-driven inference of microbial interactions based on abundance correlations is widely used, but with the drawback that ecological interpretations may not be possible. Leveraging cross-sectional microbiome datasets for unravelling ecological structure in a scalable manner thus remains an open problem. We present an expectation-maximization algorithm (BEEM-Static) that can be applied to cross-sectional datasets to infer interaction networks based on an ecological model (generalized Lotka-Volterra). The method exhibits robustness to violations in model assumptions by using statistical filters to identify and remove corresponding samples. Benchmarking against 10 state-of-the-art correlation based methods showed that BEEM-Static can infer presence and directionality of ecological interactions even with relative abundance data (AUC-ROC>0.85), a task that other methods struggle with (AUC-ROC<0.63). In addition, BEEM-Static can tolerate a high fraction of samples (up to 40%) being not at steady state or coming from an alternate model. Applying BEEM-Static to a large public dataset of human gut microbiomes (n = 4,617) identified multiple stable equilibria that better reflect ecological enterotypes with distinct carrying capacities and interactions for key species. Conclusion BEEM-Static provides new opportunities for mining ecologically interpretable interactions and systems insights from the growing corpus of microbiome data.

MotivationThe structure and function of diverse microbial communities is underpinned by ecological interactions that remain uncharacterized. With rapid adoption of metagenomic sequencing for studying microbiomes, data-driven inference of microbial interactions based on abundance correlations is widely used, but with the drawback that ecological interpretations may not be possible. Leveraging cross-sectional metagenomic datasets for unravelling ecological structure in a scalable manner thus remains an open problem.MethodsWe present an expectation-maximization algorithm (BEEM-Static) that can be applied to cross-sectional datasets to infer interaction networks based on an ecological model (generalized Lotka-Volterra). The method exhibits robustness to violations in model assumptions by using statistical filters to identify and remove corresponding samples.ResultsBenchmarking against 10 state-of-the-art correlation based methods showed that BEEM-Static can infer presence and directionality of ecological interactions even with relative abundance data (AUC-ROC>0.85), a task that other methods struggle with (AUC-ROC<0.63). In addition, BEEM-Static can tolerate a high fraction of samples (up to 40%) being not at steady state or coming from an alternate model. Applying BEEM-Static to a large public dataset of human gut microbiomes (n=4,617) identified multiple stable equilibria that better reflect ecological enterotypes with distinct carrying capacities and interactions for key species.ConclusionBEEM-Static provides new opportunities for mining ecologically interpretable interactions and systems insights from the growing corpus of metagenomic data.