Search

Background Whole-genome sequencing (WGS) is increasingly used to map the spread of bacterial and viral pathogens in nosocomial settings. A limiting factor for more widespread adoption of WGS for hospital infection prevention practices is the availability of standardized tools for genomic epidemiology. Methods We developed the Pathogen Sequencing Phylogenomic Outbreak Toolkit (PathoSPOT) to automate integration of genomic and medical record data for rapid detection and tracing of nosocomial outbreaks. To demonstrate its capabilities, we applied PathoSPOT to complete genome surveillance data of 197 MRSA bacteremia cases from two hospitals during a 2-year period. Results PathoSPOT identified 8 clonal clusters encompassing 33 patients (16.8% of cases), none of which had been recognized by standard practices. The largest cluster corresponded to a prolonged outbreak of a hospital-associated MRSA clone among 16 adults, spanning 9 wards over a period of 21 months. Analysis of precise timeline and location data with our toolkit suggested that an initial exposure event in a single ward led to infection and long-term colonization of multiple patients, followed by transmissions to other patients during recurrent hospitalizations. Conclusions We demonstrate that PathoSPOT genomic surveillance enables the detection of complex transmission chains that are not readily apparent from epidemiological data and that contribute significantly to morbidity and mortality, enabling more effective intervention strategies.

Background Nosocomial respiratory virus outbreaks represent serious public hgealth challenges. Rapid and precise identification of cases and tracing of transmission chains is critical to end outbreaks and to inform prevention measures. Methods We combined conventional surveillance with influenza A virus (IAV) genome sequencing to identify and contain a large IAV outbreak in a metropolitan healthcare system. A total of 381 individuals, including 91 inpatients and 290 healthcare workers (HCWs), were included in the investigation. Results During a 12-day period in early 2019, infection preventionists identified 89 HCWs and 18 inpatients as cases of influenza-like illness (ILI), using an amended definition without the requirement for fever. Sequencing of IAV genomes from available nasopharyngeal specimens identified 66 individuals infected with a nearly identical strain of influenza A H1N1pdm09 (43 HCWs, 17 inpatients, and 6 with unspecified affiliation). All HCWs infected with the outbreak strain had received the seasonal influenza virus vaccination. Characterization of 5 representative outbreak viral isolates did not show antigenic drift. In conjunction with IAV genome sequencing, mining of electronic records pinpointed the origin of the outbreak as a single patient and a few interactions in the emergency department that occurred 1 day prior to the index ILI cluster. Conclusions We used precision surveillance to delineate a large nosocomial IAV outbreak, mapping the source of the outbreak to a single patient rather than HCWs as initially assumed based on conventional epidemiology. These findings have important ramifications for more-effective prevention strategies to curb nosocomial respiratory virus outbreaks.

BACKGROUNDReported per-patient costs ofClostridium difficileinfection (CDI) vary by 2 orders of magnitude among different hospitals, implying that infection control officers need precise, local analyses to guide rational decision making between interventions.OBJECTIVEWe sought to comprehensively estimate changes in length of stay (LOS) attributable to CDI at a single urban tertiary-care facility using only data automatically extractable from the electronic medical record (EMR).METHODSWe performed a retrospective cohort study of 171,938 visits spanning a 7-year period. In total, 23,968 variables were extracted from EMR data recorded within 24 hours of admission to train elastic-net regularized logistic regression models for propensity score matching. To address time-dependent bias (reverse causation), we separately stratified comparisons by time of infection, and we fit multistate models.RESULTSThe estimated difference in median LOS for propensity-matched cohorts varied from 3.1 days (95% CI, 2.2–3.9) to 10.1 days (95% CI, 7.3–12.2) depending on the case definition; however, dependency of the estimate on time to infection was observed. Stratification by time to first positive toxin assay, excluding probable community-acquired infections, showed a minimum excess LOS of 3.1 days (95% CI, 1.7–4.4). Under the same case definition, the multistate model averaged an excess LOS of 3.3 days (95% CI, 2.6–4.0).CONCLUSIONSIn this study, 2 independent time-to-infection adjusted methods converged on similar excess LOS estimates. Changes in LOS can be extrapolated to marginal dollar costs by multiplying by average costs of an inpatient day. Infection control officers can leverage automatically extractable EMR data to estimate costs of CDI at their own institutions.Infect Control Hosp Epidemiol.2017;38:1478–1486

Whole-genome sequencing (WGS) of Staphylococcus aureus is increasingly used as part of infection prevention practices. In this study, we established a long-read technology-based WGS screening program of all first-episode methicillin-resistant Staphylococcus aureus (MRSA) blood infections at a major urban hospital., Whole-genome sequencing (WGS) of Staphylococcus aureus is increasingly used as part of infection prevention practices. In this study, we established a long-read technology-based WGS screening program of all first-episode methicillin-resistant Staphylococcus aureus (MRSA) blood infections at a major urban hospital. A survey of 132 MRSA genomes assembled from long reads enabled detailed characterization of an outbreak lasting several months of a CC5/ST105/USA100 clone among 18 infants in a neonatal intensive care unit (NICU). Available hospital-wide genome surveillance data traced the origins of the outbreak to three patients admitted to adult wards during a 4-month period preceding the NICU outbreak. The pattern of changes among complete outbreak genomes provided full spatiotemporal resolution of its progression, which was characterized by multiple subtransmissions and likely precipitated by equipment sharing between adults and infants. Compared to other hospital strains, the outbreak strain carried distinct mutations and accessory genetic elements that impacted genes with roles in metabolism, resistance, and persistence. This included a DNA recognition domain recombination in the hsdS gene of a type I restriction modification system that altered DNA methylation. Transcriptome sequencing (RNA-Seq) profiling showed that the (epi)genetic changes in the outbreak clone attenuated agr gene expression and upregulated genes involved in stress response and biofilm formation. Overall, our findings demonstrate the utility of long-read sequencing for hospital surveillance and for characterizing accessory genomic elements that may impact MRSA virulence and persistence.

Whole-genome sequencing (WGS) of Staphylococcus aureus is increasingly used as part of infection prevention practices, but most applications are focused on conserved core genomic regions due to limitations of short-read technologies. In this study we established a long-read technology-based WGS screening program of all first-episode MRSA blood infections at a major urban hospital. A survey of 132 MRSA genomes assembled from long reads revealed widespread gain/loss of accessory mobile genetic elements among established hospital- and community-associated lineages impacting >10% of each genome, and frequent megabase-scale inversions between endogenous prophages. We also characterized an outbreak of a CC5/ST105/USA100 clone among 3 adults and 18 infants in a neonatal intensive care unit (NICU) lasting 7 months. The pattern of changes among complete outbreak genomes provided full spatiotemporal resolution of its origins and progression, which was characterized by multiple sub-transmissions and likely precipitated by equipment sharing. Compared to other hospital strains, the outbreak strain carried distinct mutations and accessory genetic elements that impacted genes with roles in metabolism, resistance and persistence. This included a DNA-recognition domain recombination in the hsdS gene of a Type-I restriction-modification system that altered DNA methylation. RNA-Seq profiling showed that the (epi)genetic changes in the outbreak clone attenuated agr gene expression and upregulated genes involved in stress response and biofilm formation. Overall our findings demonstrate that long-read sequencing substantially improves our ability to characterize accessory genomic elements that impact MRSA virulence and persistence, and provides valuable information for infection control efforts.

Therapy for bacteremia caused by Staphylococcus aureus is often ineffective, even when treatment conditions are optimal according to experimental protocols. Adapted subclones, such as those bearing mutations that attenuate agr-mediated virulence activation, are associated with persistent infection and patient mortality., Therapy for bacteremia caused by Staphylococcus aureus is often ineffective, even when treatment conditions are optimal according to experimental protocols. Adapted subclones, such as those bearing mutations that attenuate agr-mediated virulence activation, are associated with persistent infection and patient mortality. To identify additional alterations in agr-defective mutants, we sequenced and assembled the complete genomes of clone pairs from colonizing and infected sites of several patients in whom S. aureus demonstrated a within-host loss of agr function. We report that events associated with agr inactivation result in agr-defective blood and nares strain pairs that are enriched in mutations compared to pairs from wild-type controls. The random distribution of mutations between colonizing and infecting strains from the same patient, and between strains from different patients, suggests that much of the genetic complexity of agr-defective strains results from prolonged infection or therapy-induced stress. However, in one of the agr-defective infecting strains, multiple genetic changes resulted in increased virulence in a murine model of bloodstream infection, bypassing the mutation of agr and raising the possibility that some changes were selected. Expression profiling correlated the elevated virulence of this agr-defective mutant to restored expression of the agr-regulated ESAT6-like type VII secretion system, a known virulence factor. Thus, additional mutations outside the agr locus can contribute to diversification and adaptation during infection by S. aureus agr mutants associated with poor patient outcomes.

Chikungunya virus (CHIKV) is a mosquito‐borne alphavirus that causes global epidemics of debilitating disease worldwide. To gain functional insight into the host cellular genes required for virus infection, we performed whole‐blood RNA‐seq, 37‐plex mass cytometry of peripheral blood mononuclear cells (PBMCs), and serum cytokine measurements of acute‐ and convalescent‐phase samples obtained from 42 children naturally infected with CHIKV. Semi‐supervised classification and clustering of single‐cell events into 57 sub‐communities of canonical leukocyte phenotypes revealed a monocyte‐driven response to acute infection, with the greatest expansions in “intermediate” CD14++ CD16+ monocytes and an activated subpopulation of CD14+ monocytes. Increases in acute‐phase CHIKV envelope protein E2 expression were highest for monocytes and dendritic cells. Serum cytokine measurements confirmed significant acute‐phase upregulation of monocyte chemoattractants. Distinct transcriptomic signatures were associated with infection timepoint, as well as convalescent‐phase anti‐CHIKV antibody titer, acute‐phase viremia, and symptom severity. We present a multiscale network that summarizes all observed modulations across cellular and transcriptomic levels and their interactions with clinical outcomes, providing a uniquely global view of the biomolecular landscape of human CHIKV infection.

Clostridioides difficileis a leading cause of health care-associated infections. Although significant progress has been made in the understanding of its genome, the epigenome ofC. difficileand its functional impact has not been systematically explored. Here, we performed the first comprehensive DNA methylome analysis ofC. difficileusing 36 human isolates and observed great epigenomic diversity. We discovered an orphan DNA methyltransferase with a well-defined specificity whose corresponding gene is highly conserved across our dataset and in all ~300 globalC. difficilegenomes examined. Inactivation of the methyltransferase gene negatively impacted sporulation, a key step inC. difficiledisease transmission, consistently supported by multi-omics data, genetic experiments, and a mouse colonization model. Further experimental and transcriptomic analysis also suggested that epigenetic regulation is associated with cell length, biofilm formation, and host colonization. These findings open up a new epigenetic dimension to characterize medically relevant biological processes in this critical pathogen. This work also provides a set of methods for comparative epigenomics and integrative analysis, which we expect to be broadly applicable to bacterial epigenomics studies.

Clostridioides (formerly Clostridium) difficile is a leading cause of healthcare-associated infections. Although considerable progress has been made in the understanding of its genome, the epigenome of C. difficile and its functional impact has not been systematically explored. Here, we perform a comprehensive DNA methylome analysis of C. difficile using 36 human isolates and observe a high level of epigenomic diversity. We discovered an orphan DNA methyltransferase with a well-defined specificity, the corresponding gene of which is highly conserved across our dataset and in all of the approximately 300 global C. difficile genomes examined. Inactivation of the methyltransferase gene negatively impacts sporulation, a key step in C. difficile disease transmission, and these results are consistently supported by multiomics data, genetic experiments and a mouse colonization model. Further experimental and transcriptomic analyses suggest that epigenetic regulation is associated with cell length, biofilm formation and host colonization. These findings provide a unique epigenetic dimension to characterize medically relevant biological processes in this important pathogen. This study also provides a set of methods for comparative epigenomics and integrative analysis, which we expect to be broadly applicable to bacterial epigenomic studies. In this work, Fang et al. analyse the epigenetic landscape of Clostridioides difficile and identify a DNA methyltransferase present across C. difficile strains that is required for optimal sporulation and in vivo colonization and disease.

Whole-genome sequencing was used to examine a persistent Enterococcus faecium bacteremia that acquired heteroresistance to three antibiotics in response to prolonged multidrug therapy. A comparison of the complete genomes before and after each change revealed the emergence of known resistance determinants for vancomycin and linezolid and suggested that a novel mutation in fabF , encoding a fatty acid synthase, was responsible for daptomycin nonsusceptibility. Plasmid recombination contributed to the progressive loss of vancomycin resistance after withdrawal of the drug.

Dengue virus (DENV) is the most prevalent mosquito-borne virus causing human disease. Of the 4 DENV serotypes, epidemiological data suggest that DENV-2 secondary infections are associated with more severe disease than DENV-4 infections. Mass cytometry by time-of-flight (CyTOF) was used to dissect immune changes induced by DENV-2 and DENV-4 in human DCs, the initial targets of primary infections that likely affect infection outcomes. Strikingly, DENV-4 replication peaked earlier and promoted stronger innate immune responses, with increased expression of DC activation and migration markers and increased cytokine production, compared with DENV-2. In addition, infected DCs produced higher levels of inflammatory cytokines compared with bystander DCs, which mainly produced IFN-induced cytokines. These high-dimensional analyses during DENV-2 and DENV-4 infections revealed distinct viral signatures marked by different replication strategies and antiviral innate immune induction in DCs, which may result in different viral fitness, transmission, and pathogenesis.

Summary: Current web-based genome browsers require repetitious user input to scroll over long distances, alter the drawing density of elements or zoom through multiple orders of magnitude. Generally, either the server or the client is responsible for the majority of data processing, resulting in either servers having to receive and handle data relevant only to one user, or clients redundantly processing widely viewed data. ChromoZoom pre-renders and caches general-use tracks into tiled images on the server and serves them in an interactive web interface with inertial scrolling and precise, fluent zooming via the mouse wheel or trackpad. Custom tracks in several formats can be rendered by client-side code alongside the pre-rendered tracks, minimizing server load because of user-specific rendering and eliminating the need to transmit private data. ChromoZoom thereby enables rapid and simultaneous exploration of curated, experimental and personal genomic datasets. Availability: Human and yeast genome researchers may browse recent assemblies within ChromoZoom at http://chromozoom.org/. Source code is available at http://github.com/rothlab/chromozoom/. Contact: fritz.roth@utoronto.ca Supplementary information: Supplementary data are available at Bioinformatics online.

In a liver transplant recipient with vancomycin-resistant Enterococcus (VRE) surgical site and bloodstream infection, a combination of pulsed-field gel electrophoresis, multilocus sequence typing, and whole genome sequencing identified that donor and recipient VRE isolates were highly similar when compared to time-matched hospital isolates. Comparison of de novo assembled isolate genomes was highly suggestive of transplant transmission rather than hospital-acquired transmission and also identified subtle internal rearrangements between donor and recipient missed by other genomic approaches. Given the improved resolution, whole-genome assembly of pathogen genomes is likely to become an essential tool for investigation of potential organ transplant transmissions.

Genotypic differences greatly influence susceptibility and resistance to disease. Understanding genotype-phenotype relationships requires that phenotypes be viewed as manifestations of network properties, rather than simply as the result of individual genomic variations. Genome sequencing efforts have identified numerous germline mutations, and large numbers of somatic genomic alterations, associated with a predisposition to cancer. However, it remains difficult to distinguish background, or 'passenger', cancer mutations from causal, or 'driver', mutations in these data sets. Human viruses intrinsically depend on their host cell during the course of infection and can elicit pathological phenotypes similar to those arising from mutations. Here we test the hypothesis that genomic variations and tumour viruses may cause cancer through related mechanisms, by systematically examining host interactome and transcriptome network perturbations caused by DNA tumour virus proteins. The resulting integrated viral perturbation data reflects rewiring of the host cell networks, and highlights pathways, such as Notch signalling and apoptosis, that go awry in cancer. We show that systematic analyses of host targets of viral proteins can identify cancer genes with a success rate on a par with their identification through functional genomics and large-scale cataloguing of tumour mutations. Together, these complementary approaches increase the specificity of cancer gene identification. Combining systems-level studies of pathogen-encoded gene products with genomic approaches will facilitate the prioritization of cancer-causing driver genes to advance the understanding of the genetic basis of human cancer.

Variable tandem repeats are frequently used for genetic mapping, genotyping, and forensics studies. Moreover, variation in some repeats underlies rapidly evolving traits or certain diseases. However, mutation rates vary greatly from repeat to repeat, and as a consequence, not all tandem repeats are suitable genetic markers or interesting unstable genetic modules. We developed a model, “SERV,” that predicts the variability of a broad range of tandem repeats in a wide range of organisms. The nonlinear model uses three basic characteristics of the repeat (number of repeated units, unit length, and purity) to produce a numeric “VARscore” that correlates with repeat variability. SERV was experimentally validated using a large set of different artificial repeats located in the Saccharomyces cerevisiae URA3 gene. Further in silico analysis shows that SERV outperforms existing models and accurately predicts repeat variability in bacteria and eukaryotes, including plants and humans. Using SERV, we demonstrate significant enrichment of variable repeats within human genes involved in transcriptional regulation, chromatin remodeling, morphogenesis, and neurogenesis. Moreover, SERV allows identification of known and candidate genes involved in repeat-based diseases. In addition, we demonstrate the use of SERV for the selection and comparison of suitable variable repeats for genotyping and forensic purposes. Our analysis indicates that tandem repeats used for genotyping should have a VARscore between 1 and 3. SERV is publicly available from http://hulsweb1.cgr.harvard.edu/SERV/.