Global emergence of a hypervirulent carbapenem-resistant Escherichia coli ST410 clone

Carbapenem-resistant Escherichia coli (CREC) ST410 has recently emerged as a major global health problem. Here, we report a shift in CREC prevalence in Chinese hospitals between 2017 and 2021 with ST410 becoming the most commonly isolated sequence type. Genomic analysis identifies a hypervirulent CREC ST410 clone, B5/H24RxC, which caused two separate outbreaks in a children's hospital. It may have emerged from the previously characterised B4/H24RxC in 2006 and has been isolated in ten other countries from 2015 to 2021. Compared with B4/H24RxC, B5/H24RxC lacks the blaOXA-181-bearing X3 plasmid, but carries a F-type plasmid containing blaNDM-5. Most of B5/H24RxC also carry a high pathogenicity island and a novel O-antigen gene cluster. We find that B5/H24RxC grew faster in vitro and is more virulent in vivo. The identification of this newly emerged but already globally disseminated hypervirulent CREC clone, highlights the ongoing evolution of ST410 towards increased resistance and virulence.

TheiaEuk: a species-agnostic bioinformatics workflow for fungal genomic characterization

Introduction

The clinical incidence of antimicrobial-resistant fungal infections has dramatically increased in recent years. Certain fungal pathogens colonize various body cavities, leading to life-threatening bloodstream infections. However, the identification and characterization of fungal isolates in laboratories remain a significant diagnostic challenge in medicine and public health. Whole-genome sequencing provides an unbiased and uniform identification pipeline for fungal pathogens but most bioinformatic analysis pipelines focus on prokaryotic species. To this end, TheiaEuk_Illumina_PE_PHB (TheiaEuk) was designed to focus on genomic analysis specialized to fungal pathogens.

Methods

TheiaEuk was designed using containerized components and written in the workflow description language (WDL) to facilitate deployment on the cloud-based open bioinformatics platform Terra. This species-agnostic workflow enables the analysis of fungal genomes without requiring coding, thereby reducing the entry barrier for laboratory scientists. To demonstrate the usefulness of this pipeline, an ongoing outbreak of C. auris in southern Nevada was investigated. We performed whole-genome sequence analysis of 752 new C. auris isolates from this outbreak. Furthermore, TheiaEuk was utilized to observe the accumulation of mutations in the FKS1 gene over the course of the outbreak, highlighting the utility of TheiaEuk as a monitor of emerging public health threats when combined with whole-genome sequencing surveillance of fungal pathogens.

Results

A primary result of this work is a curated fungal database containing 5,667 unique genomes representing 245 species. TheiaEuk also incorporates taxon-specific submodules for specific species, including clade-typing for Candida auris (C. auris). In addition, for several fungal species, it performs dynamic reference genome selection and variant calling, reporting mutations found in genes currently associated with antifungal resistance (FKS1, ERG11, FUR1). Using genome assemblies from the ATCC Mycology collection, the taxonomic identification module used by TheiaEuk correctly assigned genomes to the species level in 126/135 (93.3%) instances and to the genus level in 131/135 (97%) of instances, and provided zero false calls. Application of TheiaEuk to actual specimens obtained in the course of work at a local public health laboratory resulted in 13/15 (86.7%) correct calls at the species level, with 2/15 called at the genus level. It made zero incorrect calls. TheiaEuk accurately assessed clade type of Candida auris in 297/302 (98.3%) of instances.

Discussion

TheiaEuk demonstrated effectiveness in identifying fungal species from whole genome sequence. It further showed accuracy in both clade-typing of C. auris and in the identification of mutations known to associate with drug resistance in that organism.

Endemicity and diversification of carbapenem-resistant Acinetobacter baumannii in an intensive care unit

Background

Carbapenem-resistant Acinetobacter baumannii (CRAB) is a major public health concern globally. Often studied in the context of hospital outbreaks, little is known about the persistence and evolutionary dynamics of endemic CRAB populations.

Methods

A three-month cross-sectional observational study was conducted in a 28-bed intensive care unit (ICU) in Hangzhou, China. A total of 5068 samples were collected from the hospital environment (n = 3985), patients (n = 964) and staff (n = 119). CRAB isolates were obtained from 10.5% of these samples (n = 532). All of these isolates, plus an additional 19 from clinical infections, were characterised through whole-genome sequencing.

Findings

The ICU CRAB population was dominated by OXA-23-producing global clone 2 isolates (99.3% of all isolates) that could be divided into 20 distinct clusters, defined through genome sequencing. CRAB was persistently present in the ICU, driven by regular introductions of distinct clusters. The hospital environment was heavily contaminated, with CRAB isolated from bed units on 183/335 (54.6%) sampling occasions but from patients on only 72/299 (24.1%) occasions. CRAB was spread to adjacent bed units and rooms, and following re-location of patients within the ICU. We also observed three horizontal gene transfer events between CRAB strains in the ICU, involving three different plasmids.

Interpretation

The epidemiology of CRAB in this setting contrasted with previously described clonal outbreaks in high-income countries, highlighting the importance of environmental CRAB reservoirs in ICU epidemiology and the unique challenges in containing the spread of CRAB in ICUs where this important multidrug-resistant pathogen is endemic.

Funding

This work was undertaken as part of the DETECTIVE research project funded by the Medical Research Council (MR/S013660/1), National Natural Science Foundation of China (81861138054, 32011530116, 31970128, 31770142), Zhejiang Province Medical Platform Backbone Talent Plan (2020RC075), and the National Key Research and Development Program of China grant (2018YFE0102100). W.v.S was also supported by a Wolfson Research Merit Award (WM160092).

Accelerating bioinformatics implementation in public health

We have adopted an open bioinformatics ecosystem to address the challenges of bioinformatics implementation in public health laboratories (PHLs). Bioinformatics implementation for public health requires practitioners to undertake standardized bioinformatic analyses and generate reproducible, validated and auditable results. It is essential that data storage and analysis are scalable, portable and secure, and that implementation of bioinformatics fits within the operational constraints of the laboratory. We address these requirements using Terra, a web-based data analysis platform with a graphical user interface connecting users to bioinformatics analyses without the use of code. We have developed bioinformatics workflows for use with Terra that specifically meet the needs of public health practitioners. These Theiagen workflows perform genome assembly, quality control, and characterization, as well as construction of phylogeny for insights into genomic epidemiology. Additonally, these workflows use open-source containerized software and the WDL workflow language to ensure standardization and interoperability with other bioinformatics solutions, whilst being adaptable by the user. They are all open source and publicly available in Dockstore with the version-controlled code available in public GitHub repositories. They have been written to generate outputs in standardized file formats to allow for further downstream analysis and visualization with separate genomic epidemiology software. Testament to this solution meeting the requirements for bioinformatic implementation in public health, Theiagen workflows have collectively been used for over 5 million sample analyses in the last 2 years by over 90 public health laboratories in at least 40 different countries. Continued adoption of technological innovations and development of further workflows will ensure that this ecosystem continues to benefit PHLs.

Extended-Spectrum β-Lactamase Genes Traverse the Escherichia coli Populations of Intensive Care Unit Patients, Staff, and Environment

Over a 3-month period, we monitored the population of extended-spectrum β-lactam-resistant Escherichia coli (ESBL-EC) associated with the patients, staff, and environment of an intensive care unit (ICU) in Guangzhou, China. Thirty-four clinical isolates were obtained from the same hospital 12 months later. A total of 165 isolates were characterized and whole-genome sequenced, with 24 isolates subjected to long-read sequencing. The diverse population included representatives of 59 different sequence types (STs). ICU patient and environmental isolates were largely distinct from staff isolates and clinical isolates. We observed five instances of highly similar isolates (0 to 13 single nucleotide polymorphisms [SNPs]) being obtained from different patients or bed unit environments. ESBL resistance in this collection was largely conferred by blaCTX-M genes, which were found in 96.4% of all isolates. The contexts of blaCTX-M genes were diverse, situated in multiple chromosomal positions and in various plasmids. We identified blaCTX-M-bearing plasmid lineages that were present in multiple STs across the surveillance, staff, and clinical collections. Closer examination of ISEcp1-blaCTX-M transposition units shed light on the dynamics of their transmission, with evidence for the acquisition of chromosomal copies of blaCTX-M genes from specific plasmid lineages and for the movement of blaCTX-M-55 from a ST1193 chromosome to a small mobilizable plasmid. A carbapenem-resistant ST167 strain isolated from a patient that had been treated with meropenem and piperacillin-tazobactam contained seven copies of blaCMY-146, which appears to have been amplified by IS1. Our data revealed limited persistence and movement of ESBL-EC strains in the ICU environment, but we observed circulating plasmid lineages playing an essential and ongoing role in shaping the cephalosporin-resistance landscape in the population examined. IMPORTANCE ESBL resistance significantly impacts clinical management of E. coli infections in hospitals globally. It is important to understand the structures of ESBL-EC populations carried by hospital patients and staff, their capacity to persist in hospital environments, and the dynamics of mobile genes that drive the spread of ESBL resistance. In our 3-month study, ESBL-EC strains found in the ICU environment were strongly associated with patient carriage but distinct from strains found in staff. However, plasmid lineages carrying blaCTX-M genes were found across the ICU populations and in a collection of clinical isolates obtained 1 year later. By examining their content and contexts, we have traced the recent histories of chromosomal and plasmid-borne ISEcp1-blaCTX-M transposition units in the ICU population. This information allowed us to implicate specific plasmid lineages in the acquisition of chromosomal blaCTX-M genes, even when the plasmids were no longer present, and to detect recent transposition of blaCTX-M-55 from a chromosome to a mobilizable plasmid. Similar high-resolution approaches to the study of mobile genetic elements will be essential if the transmission routes associated with the spread of ESBL resistance are to be understood and subjected to interventions.

GR13-type plasmids in Acinetobacter potentiate the accumulation and horizontal transfer of diverse accessory genes

Carbapenem and other antibiotic resistance genes (ARGs) can be found in plasmids in Acinetobacter, but many plasmid types in this genus have not been well-characterized. Here we describe the distribution, diversity and evolutionary capacity of rep group 13 (GR13) plasmids that are found in Acinetobacter species from diverse environments. Our investigation was prompted by the discovery of two GR13 plasmids in A. baumannii isolated in an intensive care unit (ICU). The plasmids harbour distinct accessory genes: pDETAB5 contains bla_NDM-1 and genes that confer resistance to four further antibiotic classes, while pDETAB13 carries putative alcohol tolerance determinants. Both plasmids contain multiple dif modules, which are flanked by pdif sites recognized by XerC/XerD tyrosine recombinases. The ARG-containing dif modules in pDETAB5 are almost identical to those found in pDETAB2, a GR34 plasmid from an unrelated A. baumannii isolated in the same ICU a month prior. Examination of a further 41 complete, publicly available plasmid sequences revealed that the GR13 pangenome consists of just four core but 1186 accessory genes, 123 in the shell and 1063 in the cloud, reflecting substantial capacity for diversification. The GR13 core genome includes genes for replication and partitioning, and for a putative tyrosine recombinase. Accessory segments encode proteins with diverse putative functions, including for metabolism, antibiotic/heavy metal/alcohol tolerance, restriction-modification, an anti-phage system and multiple toxin–antitoxin systems. The movement of dif modules and actions of insertion sequences play an important role in generating diversity in GR13 plasmids. Discrete GR13 plasmid lineages are internationally disseminated and found in multiple Acinetobacter species, which suggests they are important platforms for the accumulation, horizontal transmission and persistence of accessory genes in this genus.

Transferable Acinetobacter baumannii plasmid pDETAB2 encodes OXA-58 and NDM-1 and represents a new class of antibiotic resistance plasmids

OBJECTIVES

To characterize a blaOXA-58- and blaNDM-1-containing MDR plasmid from a rare Acinetobacter baumannii lineage and compare it with related plasmids to explore the distribution and evolution of a new plasmid group.

METHODS

A. baumannii DETAB-P2 was isolated from a rectal swab of an intensive care patient. Antibiotic susceptibility was determined using broth microdilution. DETAB-P2 was mated with A. baumannii ATCC 17978 and putative transconjugants were characterized by S1/PFGE and Southern hybridization. WGS was performed on the Illumina and Oxford Nanopore platforms. MLST was performed with the Pasteur and Oxford schemes. Antibiotic resistance genes were identified with ABRicate. Plasmid sequence annotation was performed manually. Complete plasmids in GenBank with the same rep gene were used for comparative analyses.

RESULTS

A. baumannii DETAB-P2 was ST138 by the Pasteur scheme and a novel Oxford type, ST2209. It transferred blaOXA-58 and blaNDM-1 to ATCC 17978 in the 100 072 bp plasmid pDETAB2 that also carried bleMBL, sul2, aacC2d, tet(39), msr(E)-mph(E) and putative mercury resistance and RND efflux system determinants. pDETAB2 represents a new plasmid type, GR34, and contained 16 pdif sites and several novel dif modules. Only a 10 kbp core sequence is shared amongst pDETAB2 and 18 further GR34 plasmids in GenBank, with diverse accessory regions comprised of various dif modules.

CONCLUSIONS

GR34 plasmids are found in several Acinetobacter species from diverse environments. They display considerable variation in accessory content owing to the presence of pdif sites and an array of dif modules, some of which contain antibiotic resistance genes.

Estimating the time-varying reproduction number of SARS-CoV-2 using national and subnational case counts

Background: Assessing temporal variations in transmission in different countries is essential for monitoring the epidemic, evaluating the effectiveness of public health interventions and estimating the impact of changes in policy. Methods: We use case and death notification data to generate daily estimates of the time-varying reproduction number globally, regionally, nationally, and subnationally over a 12-week rolling window. Our modelling framework, based on open source tooling, accounts for uncertainty in reporting delays, so that the reproduction number is estimated based on underlying latent infections. Results: Estimates of the reproduction number, trajectories of infections, and forecasts are displayed on a dedicated website as both maps and time series, and made available to download in tabular form. Conclusions:  This decision-support tool can be used to assess changes in virus transmission both globally, regionally, nationally, and subnationally. This allows public health officials and policymakers to track the progress of the outbreak in near real-time using an epidemiologically valid measure. As well as providing regular updates on our website, we also provide an open source tool-set so that our approach can be used directly by researchers and policymakers on confidential data-sets. We hope that our tool will be used to support decisions in countries worldwide throughout the ongoing COVID-19 pandemic.

Evolution and dissemination of L and M plasmid lineages carrying antibiotic resistance genes in diverse Gram-negative bacteria

Conjugative, broad host-range plasmids of the L/M complex have been associated with antibiotic resistance since the 1970s. They are found in Gram-negative bacterial genera that cause human infections and persist in hospital environments. It is crucial that these plasmids are typed accurately so that their clinical and global dissemination can be traced in epidemiological studies. The L/M complex has previously been divided into L, M1 and M2 subtypes. However, those types do not encompass all diversity seen in the group. Here, we have examined 148 complete L/M plasmid sequences in order to understand the diversity of the complex and trace the evolution of distinct lineages. The backbone sequence of each plasmid was determined by removing translocatable genetic elements and reversing their effects in silico. The sequence identities of replication regions and complete backbones were then considered for typing. This supported the distinction of L and M plasmids and revealed that there are five L and eight M types, where each type is comprised of further sub-lineages that are distinguished by variation in their backbone and translocatable element content. Regions containing antibiotic resistance genes in L and M sub-lineages have often formed by initial rare insertion events, followed by insertion of other translocatable elements within the inceptive element. As such, islands evolve in situ to contain genes conferring resistance to multiple antibiotics. In some cases, different plasmid sub-lineages have acquired the same or related resistance genes independently. This highlights the importance of these plasmids in acting as vehicles for the dissemination of emerging resistance genes. Materials are provided here for typing plasmids of the L/M complex from complete sequences or draft genomes. This should enable rapid identification of novel types and facilitate tracking the evolution of existing lineages.

Estimating the time-varying reproduction number of SARS-CoV-2 using national and subnational case counts

Background: Interventions are now in place worldwide to reduce transmission of the novel coronavirus. Assessing temporal variations in transmission in different countries is essential for evaluating the effectiveness of public health interventions and the impact of changes in policy. Methods: We use case notification data to generate daily estimates of the time-dependent reproduction number in different regions and countries. Our modelling framework, based on open source tooling, accounts for reporting delays, so that temporal variations in reproduction number estimates can be compared directly with the times at which interventions are implemented. Results: We provide three example uses of our framework. First, we demonstrate how the toolset displays temporal changes in the reproduction number. Second, we show how the framework can be used to reconstruct case counts by date of infection from case counts by date of notification, as well as to estimate the reproduction number. Third, we show how maps can be generated to clearly show if case numbers are likely to decrease or increase in different regions. Results are shown for regions and countries worldwide on our website (https://epiforecasts.io/covid/) and are updated daily. Our tooling is provided as an open-source R package to allow replication by others. Conclusions: This decision-support tool can be used to assess changes in virus transmission in different regions and countries worldwide. This allows policymakers to assess the effectiveness of current interventions, and will be useful for inferring whether or not transmission will increase when interventions are lifted. As well as providing daily updates on our website, we also provide adaptable computing code so that our approach can be used directly by researchers and policymakers on confidential datasets. We hope that our tool will be used to support decisions in countries worldwide throughout the ongoing COVID-19 pandemic.

Whole-genome sequencing illuminates the evolution and spread of multidrug-resistant tuberculosis in Southwest Nigeria

Nigeria has an emerging problem with multidrug-resistant tuberculosis (MDR-TB). Whole-genome sequencing was used to understand the epidemiology of tuberculosis and genetics of multi-drug resistance among patients from two tertiary referral centers in Southwest Nigeria. In line with previous molecular epidemiology studies, most isolates of Mycobacterium tuberculosis from this dataset belonged to the Cameroon clade within the Euro-American lineage. Phylogenetic analysis showed this clade was undergoing clonal expansion in this region, and suggests that it was involved in community transmission of sensitive and multidrug-resistant tuberculosis. Five patients enrolled for retreatment were infected with pre-extensively drug resistant (pre-XDR) due to fluoroquinolone resistance in isolates from the Cameroon clade. In all five cases resistance was conferred through a mutation in the gyrA gene. In some patients, genomic changes occurred in bacterial isolates during the course of treatment that potentially led to decreased drug susceptibility. We conclude that inter-patient transmission of resistant isolates, principally from the Cameroon clade, contributes to the spread of MDR-TB in this setting, underscoring the urgent need to curb the spread of multi-drug resistance in this region.

Culture-independent detection and characterisation of Mycobacterium tuberculosis and M. africanum in sputum samples using shotgun metagenomics on a benchtop sequencer

Tuberculosis remains a major global health problem. Laboratory diagnostic methods that allow effective, early detection of cases are central to management of tuberculosis in the individual patient and in the community. Since the 1880s, laboratory diagnosis of tuberculosis has relied primarily on microscopy and culture. However, microscopy fails to provide species- or lineage-level identification and culture-based workflows for diagnosis of tuberculosis remain complex, expensive, slow, technically demanding and poorly able to handle mixed infections. We therefore explored the potential of shotgun metagenomics, sequencing of DNA from samples without culture or target-specific amplification or capture, to detect and characterise strains from the Mycobacterium tuberculosis complex in smear-positive sputum samples obtained from The Gambia in West Africa. Eight smear- and culture-positive sputum samples were investigated using a differential-lysis protocol followed by a kit-based DNA extraction method, with sequencing performed on a benchtop sequencing instrument, the Illumina MiSeq. The number of sequence reads in each sputum-derived metagenome ranged from 989,442 to 2,818,238. The proportion of reads in each metagenome mapping against the human genome ranged from 20% to 99%. We were able to detect sequences from the M. tuberculosis complex in all eight samples, with coverage of the H37Rv reference genome ranging from 0.002X to 0.7X. By analysing the distribution of large sequence polymorphisms (deletions and the locations of the insertion element IS6110) and single nucleotide polymorphisms (SNPs), we were able to assign seven of eight metagenome-derived genomes to a species and lineage within the M. tuberculosis complex. Two metagenome-derived mycobacterial genomes were assigned to M. africanum, a species largely confined to West Africa; the others that could be assigned belonged to lineages T, H or LAM within the clade of "modern" M. tuberculosis strains. We have provided proof of principle that shotgun metagenomics can be used to detect and characterise M. tuberculosis sequences from sputum samples without culture or target-specific amplification or capture, using an accessible benchtop-sequencing platform, the Illumina MiSeq, and relatively simple DNA extraction, sequencing and bioinformatics protocols. In our hands, sputum metagenomics does not yet deliver sufficient depth of coverage to allow sequence-based sensitivity testing; it remains to be determined whether improvements in DNA extraction protocols alone can deliver this or whether culture, capture or amplification steps will be required. Nonetheless, we can foresee a tipping point when a unified automated metagenomics-based workflow might start to compete with the plethora of methods currently in use in the diagnostic microbiology laboratory.