Validating a 14-Drug Microtiter Plate Containing Bedaquiline and Delamanid for Large-Scale Research Susceptibility Testing of Mycobacterium tuberculosis

The UKMYC5 plate is a 96-well microtiter plate designed by the CRyPTIC Consortium (Comprehensive Resistance Prediction for Tuberculosis: an International Consortium) to enable the measurement of MICs of 14 different antituberculosis (anti-TB) compounds for >30,000 clinical Mycobacterium tuberculosis isolates. Unlike the MYCOTB plate, on which the UKMYC5 plate is based, the UKMYC5 plate includes two new (bedaquiline and delamanid) and two repurposed (clofazimine and linezolid) compounds. UKMYC5 plates were tested by seven laboratories on four continents by use of a panel of 19 external quality assessment (EQA) strains, including H37Rv. To assess the optimal combination of reading method and incubation time, MICs were measured from each plate by two readers, using three methods (mirrored box, microscope, and Vizion digital viewing system), after 7, 10, 14, and 21 days of incubation. In addition, all EQA strains were subjected to whole-genome sequencing and phenotypically characterized by the 7H10/7H11 agar proportion method (APM) and by use of MGIT960 mycobacterial growth indicator tubes. We concluded that the UKMYC5 plate is optimally read using the Vizion system after 14 days of incubation, achieving an interreader agreement of 97.9% and intra- and interlaboratory reproducibility rates of 95.6% and 93.1%, respectively. The mirrored box had a similar reproducibility. Strains classified as resistant by APM, MGIT960, or the presence of mutations known to confer resistance consistently showed elevated MICs compared to those for strains classified as susceptible. Finally, the UKMYC5 plate records intermediate MICs for one strain for which the APM measured MICs close to the applied critical concentration, providing early evidence that the UKMYC5 plate can quantitatively measure the magnitude of resistance to anti-TB compounds that is due to specific genetic variation.

Trends over time in Escherichia coli bloodstream infections, urinary tract infections, and antibiotic susceptibilities in Oxfordshire, UK, 1998-2016: a study of electronic health records

BACKGROUND

Escherichia coli bloodstream infections are increasing in the UK and internationally. The evidence base to guide interventions against this major public health concern is small. We aimed to investigate possible drivers of changes in the incidence of E coli bloodstream infection and antibiotic susceptibilities in Oxfordshire, UK, over the past two decades, while stratifying for time since hospital exposure.

METHODS

In this observational study, we used all available data on E coli bloodstream infections and E coli urinary tract infections (UTIs) from one UK region (Oxfordshire) using anonymised linked microbiological data and hospital electronic health records from the Infections in Oxfordshire Research Database (IORD). We estimated the incidence of infections across a two decade period and the annual incidence rate ratio (aIRR) in 2016. We modelled the data using negative binomial regression on the basis of microbiological, clinical, and health-care-exposure risk factors. We investigated infection severity, 30-day all-cause mortality, and community and hospital amoxicillin plus clavulanic acid (co-amoxiclav) use to estimate changes in bacterial virulence and the effect of antimicrobial resistance on incidence.

FINDINGS

From Jan 1, 1998, to Dec 31, 2016, 5706 E coli bloodstream infections occurred in 5215 patients, and 228 376 E coli UTIs occurred in 137 075 patients. 1365 (24%) E coli bloodstream infections were nosocomial (onset >48 h after hospital admission), 1132 (20%) were quasi-nosocomial (≤30 days after discharge), 1346 (24%) were quasi-community (31-365 days after discharge), and 1863 (33%) were community (>365 days after hospital discharge). The overall incidence increased year on year (aIRR 1·06, 95% CI 1·05-1·06). In 2016, 212 (41%) of 515 E coli bloodstream infections and 3921 (28%) of 13 792 E coli UTIs were co-amoxiclav resistant. Increases in E coli bloodstream infections were driven by increases in community (aIRR 1·10, 95% CI 1·07-1·13; p<0·0001) and quasi-community (aIRR 1·08, 1·07-1·10; p<0·0001) cases. 30-day mortality associated with E coli bloodstream infection decreased over time in the nosocomial (adjusted rate ratio [RR] 0·98, 95% CI 0·96-1·00; p=0·03) group, and remained stable in the quasi-nosocomial (adjusted RR 0·98, 0·95-1·00; p=0·06), quasi-community (adjusted RR 0·99, 0·96-1·01; p=0·32), and community (adjusted RR 0·99, 0·96-1·01; p=0·21) groups. Mortality was, however, substantial at 14-25% across all hospital-exposure groups. Co-amoxiclav-resistant E coli bloodstream infections increased in all groups across the study period (by 11-18% per year, significantly faster than co-amoxiclav-susceptible E coli bloodstream infections; p_{heterogeneity}<0·0001), as did co-amoxiclav-resistant E coli UTIs (by 14-29% per year; p_{heterogeneity}<0·0001). Previous year co-amoxiclav use in primary-care facilities was associated with increased subsequent year community co-amoxiclav-resistant E coli UTIs (p=0·003).

INTERPRETATION

Increases in E coli bloodstream infections in Oxfordshire are primarily community associated, with substantial co-amoxiclav resistance; nevertheless, we found little or no change in mortality. Focusing interventions on primary care facilities, particularly those with high co-amoxiclav use, could be effective in reducing the incidence of co-amoxiclav-resistant E coli bloodstream infections, in this region and more generally.

FUNDING

National Institute for Health Research.

A Quantitative Evaluation of MIRU-VNTR Typing Against Whole-Genome Sequencing for Identifying Mycobacterium tuberculosis Transmission: A Prospective Observational Cohort Study

Background

Mycobacterial Interspersed Repetitive Unit-Variable Number Tandem Repeat (MIRU-VNTR) typing is widely used in high-income countries to determine Mycobacterium tuberculosis relatedness. Whole-genome sequencing (WGS) is known to deliver greater specificity, but no quantitative prospective comparison has yet been undertaken.

Methods

We studied isolates from the English Midlands, sampled consecutively between 1 January 2012 and 31 December 2015. In addition to routinely performed MIRU-VNTR typing, DNA was extracted from liquid cultures and sequenced using Illumina technology. Demographic and epidemiological data for the relevant patients were extracted from the Enhanced Tuberculosis Surveillance system run by Public Health England. Closely related samples, defined using a threshold of five single nucleotide variants (SNVs), were compared to samples with identical MIRU-VNTR profiles, to samples from individuals with shared epidemiological risk factors, and to those with both characteristics.

Findings

1999 patients were identified for whom at least one M. tuberculosis isolate had been MIRU-VNTR typed and sequenced. Comparing epidemiological risk factors with close genetic relatedness, only co-residence had a positive predictive value of over 5%. Excluding co-resident individuals, 18.6% of patients with identical MIRU-VNTR profiles were within 5 SNVs. Where patients also shared social risk factors and ethnic group, this rose to 48%. Only 8% of MIRU-VNTR linked pairs in lineage 1 were within 5 SNV, compared to 31% in lineage 4.

Interpretation

In the setting studied, this molecular epidemiological study shows MIRU-VNTR typing and epidemiological risk factors are poorly predictive of close genomic relatedness, assessed by SNV. MIRU-VNTR performance varies markedly by lineage.

Funding

Public Health England, Health Innovation Challenge Fund, NIHR Health Protection Research Unit Oxford, NIHR Oxford Biomedical Research Centre.

•

M. tuberculosis relatedness was measured by MIRU-VNTR typing and single nucleotide variants (SNV).

•

In consecutive isolates in England, only 20% of samples with identical MIRU-VNTR profiles were within 5 SNVs.

•

MIRU-VNTR is more predictive of close relatedness in lineage 4 (Euro-American) lineages than in other lineages.

Control of Artifactual Variation in Reported Intersample Relatedness during Clinical Use of a Mycobacterium tuberculosis Sequencing Pipeline

Contact tracing requires reliable identification of closely related bacterial isolates. When we noticed the reporting of artifactual variation between Mycobacterium tuberculosis isolates during routine next-generation sequencing of Mycobacterium spp., we investigated its basis in 2,018 consecutive M. tuberculosis isolates. In the routine process used, clinical samples were decontaminated and inoculated into broth cultures; from positive broth cultures DNA was extracted and sequenced, reads were mapped, and consensus sequences were determined. We investigated the process of consensus sequence determination, which selects the most common nucleotide at each position. Having determined the high-quality read depth and depth of minor variants across 8,006 M. tuberculosis genomic regions, we quantified the relationship between the minor variant depth and the amount of nonmycobacterial bacterial DNA, which originates from commensal microbes killed during sample decontamination. In the presence of nonmycobacterial bacterial DNA, we found significant increases in minor variant frequencies, of more than 1.5-fold, in 242 regions covering 5.1% of the M. tuberculosis genome. Included within these were four high-variation regions strongly influenced by the amount of nonmycobacterial bacterial DNA. Excluding these four regions from pairwise distance comparisons reduced biologically implausible variation from 5.2% to 0% in an independent validation set derived from 226 individuals. Thus, we demonstrated an approach identifying critical genomic regions contributing to clinically relevant artifactual variation in bacterial similarity searches. The approach described monitors the outputs of the complex multistep laboratory and bioinformatics process, allows periodic process adjustments, and will have application to quality control of routine bacterial genomics.

Gonorrhoea treatment failure caused by a Neisseria gonorrhoeae strain with combined ceftriaxone and high-level azithromycin resistance, England, February 2018

We describe a gonorrhoea case with combined high-level azithromycin resistance and ceftriaxone resistance. In February 2018, a heterosexual male was diagnosed with gonorrhoea in the United Kingdom following sexual intercourse with a locally resident female in Thailand and failed treatment with ceftriaxone plus doxycycline and subsequently spectinomycin. Resistance arose from two mechanisms combining for the first time in a genetic background similar to a commonly circulating strain. Urgent action is essential to prevent further spread.

The use of whole-genome sequencing in cluster investigation of a multidrug-resistant tuberculosis outbreak

We used whole-genome sequencing (WGS) to delineate transmission networks and investigate the benefits of WGS during cluster investigation.We included clustered cases of multidrug-resistant (MDR) tuberculosis (TB)/extensively drug-resistant (XDR) TB linked by mycobacterial interspersed repetitive unit variable tandem repeat (MIRU-VNTR) strain typing or epidemiological information in the national cluster B1006, notified between 2007 and 2013 in the UK. We excluded from further investigation cases whose isolates differed by greater than 12 single nucleotide polymorphisms (SNPs). Data relating to patients' social networks were collected.27 cases were investigated and 22 had WGS, eight of which (36%) were excluded as their isolates differed by more than 12 SNPs to other cases. 18 cases were ruled into the transmission network based on genomic and epidemiological information. Evidence of transmission was inconclusive in seven out of 18 cases (39%) in the transmission network following WGS and epidemiological investigation.This investigation of a drug-resistant TB cluster illustrates the opportunities and limitations of WGS in understanding transmission in a setting with a high proportion of migrant cases. The use of WGS should be combined with classical epidemiological methods. However, not every cluster will be solvable, regardless of the quality of genomic data.

mcr-1 and mcr-2 variant genes identified in Moraxella species isolated from pigs in Great Britain from 2014 to 2015

Objectives

To determine the occurrence of mcr-1 and mcr-2 genes in Gram-negative bacteria isolated from healthy pigs in Great Britain.

Methods

Gram-negative bacteria (n = 657) isolated from pigs between 2014 and 2015 were examined by WGS.

Results

Variants of mcr-1 and mcr-2 were identified in Moraxella spp. isolated from pooled caecal contents of healthy pigs at slaughter collected from six farms in Great Britain. Other bacteria, including Escherichia coli from the same farms, were not detected harbouring mcr-1 or mcr-2. A Moraxella porci-like isolate, MSG13-C03, harboured MCR-1.10 with 98.7% identity to MCR-1, and a Moraxella pluranimalium-like isolate, MSG47-C17, harboured an MCR-2.2 variant with 87.9% identity to MCR-2, from E. coli; the isolates had colistin MICs of 1-2 mg/L. No intact insertion elements were identified in either MSG13-C03 or MSG47-C17, although MSG13-C03 harboured the conserved nucleotides abutting the ISApl1 composite transposon found in E. coli plasmids and the intervening ∼2.6 kb fragment showed 97% identity. Six Moraxella osloensis isolates were positive for phosphoethanolamine transferase (EptA). They shared 62%-64.5% identity to MCR-1 and MCR-2, with colistin MICs from 2 to 4 mg/L. Phylogenetic analysis indicated that MCR and EptA have evolved from a common ancestor. In addition to mcr, the β-lactamase gene, blaBRO-1, was found in both isolates, whilst the tetracycline resistance gene, tetL, was found in MSG47-C17.

Conclusions

Our results add further evidence for the mobilization of the mcr-pap2 unit from Moraxella via composite transposons leading to its global dissemination. The presence of mcr-pap2 from recent Moraxella isolates indicates they may comprise a reservoir for mcr.

Comparison of Control of Clostridium difficile Infection in Six English Hospitals Using Whole-Genome Sequencing

Background.

Variation in Clostridium difficile infection (CDI) rates between healthcare institutions suggests overall incidence could be reduced if the lowest rates could be achieved more widely.

Methods.

We used whole-genome sequencing (WGS) of consecutive C. difficile isolates from 6 English hospitals over 1 year (2013–14) to compare infection control performance. Fecal samples with a positive initial screen for C. difficile were sequenced. Within each hospital, we estimated the proportion of cases plausibly acquired from previous cases.

Results.

Overall, 851/971 (87.6%) sequenced samples contained toxin genes, and 451 (46.4%) were fecal-toxin-positive. Of 652 potentially toxigenic isolates >90-days after the study started, 128 (20%, 95% confidence interval [CI] 17–23%) were genetically linked (within ≤2 single nucleotide polymorphisms) to a prior patient’s isolate from the previous 90 days. Hospital 2 had the fewest linked isolates, 7/105 (7%, 3–13%), hospital 1, 9/70 (13%, 6–23%), and hospitals 3–6 had similar proportions of linked isolates (22–26%) (P ≤ .002 comparing hospital-2 vs 3–6). Results were similar adjusting for locally circulating ribotypes. Adjusting for hospital, ribotype-027 had the highest proportion of linked isolates (57%, 95% CI 29–81%). Fecal-toxin-positive and toxin-negative patients were similarly likely to be a potential transmission donor, OR = 1.01 (0.68–1.49). There was no association between the estimated proportion of linked cases and testing rates.

Conclusions.

WGS can be used as a novel surveillance tool to identify varying rates of C. difficile transmission between institutions and therefore to allow targeted efforts to reduce CDI incidence.

The Hospital Water Environment as a Reservoir for Carbapenem-Resistant Organisms Causing Hospital-Acquired Infections-A Systematic Review of the Literature

Over the last 20 years there have been 32 reports of carbapenem-resistant organisms in the hospital water environment, with half of these occurring since 2010. The majority of these reports have described associated clinical outbreaks in the intensive care setting, affecting the critically ill and the immunocompromised. Drains, sinks, and faucets were most frequently colonized, and Pseudomonas aeruginosa the predominant organism. Imipenemase (IMP), Klebsiella pneumoniae carbapenemase (KPC), and Verona integron-encoded metallo-β-lactamase (VIM) were the most common carbapenemases found. Molecular typing was performed in almost all studies, with pulse field gel electrophoresis being most commonly used. Seventy-two percent of studies reported controlling outbreaks, of which just more than one-third eliminated the organism from the water environment. A combination of interventions seems to be most successful, including reinforcement of general infection control measures, alongside chemical disinfection. The most appropriate disinfection method remains unclear, however, and it is likely that replacement of colonized water reservoirs may be required for long-term clearance.

DNA extraction from primary liquid blood cultures for bloodstream infection diagnosis using whole genome sequencing

PURPOSE

Speed of bloodstream infection diagnosis is vital to reduce morbidity and mortality. Whole genome sequencing (WGS) performed directly from liquid blood culture could provide single-assay species and antibiotic susceptibility prediction; however, high inhibitor and human cell/DNA concentrations limit pathogen recovery. We develop a method for the preparation of bacterial DNA for WGS-based diagnostics direct from liquid blood culture.

METHODOLOGY

We evaluate three commercial DNA extraction kits: BiOstic Bacteraemia, Amplex Hyplex and MolYsis Plus. Differential centrifugation, filtration, selective lysis and solid-phase reversible immobilization bead clean-up are tested to improve human cells/DNA and inhibitor removal. Using WGS (Illumina/MinION), we assess human DNA removal, pathogen recovery, and predict species and antibiotic susceptibility inpositive blood cultures of 44 Gram-negative and 54 Staphylococcus species.Results/Key findings. BiOstic kit extractions yield the greatest mean DNA concentration, 94-301 ng µl^-1, versus 0-2.5 ng µl^-1 using Amplex and MolYsis kits. However, we note higher levels of inhibition (260/280 ratio 0.9-2.1) and human DNA (0.0-4.4×10⁶ copies) in BiOstic extracts. Differential centrifugation (2000 g, 1 min) prior to BiOstic extraction reduces human DNA by 63-89 % with selective lysis minimizing by a further 62 %. Post-extraction bead clean-up lowers inhibition. Overall, 67 % of sequenced samples (Illumina MiSeq) contain <10 % human DNA, with >93 % concordance between WGS-based species and susceptibility predictions and clinical diagnosis. If >60 % of sequencing reads are human (7/98 samples) susceptibility prediction becomes compromised. Novel MinION-based WGS (n=9) currently gives rapid species identification but not susceptibility prediction.

CONCLUSION

Our method for DNA preparation allows WGS-based diagnosis direct from blood culture bottles, providing species and antibiotic susceptibility prediction in a single assay.

Re-emergence of methicillin susceptibility in a resistant lineage of Staphylococcus aureus

Objectives

MRSA is a leading cause of hospital-associated infection. Acquired resistance is encoded by the mecA gene or its homologue mecC , but little is known about the evolutionary dynamics involved in gain and loss of resistance. The objective of this study was to obtain an expanded understanding of Staphylococcus aureus methicillin resistance microevolution in vivo , by focusing on a single lineage.

Methods

We compared the whole-genome sequences of 231 isolates from a single epidemic lineage [clonal complex 30 (CC30) and spa -type t018] of S. aureus that caused an epidemic in the UK.

Results

We show that resistance to methicillin in this single lineage was gained on at least two separate occasions, one of which led to a clonal expansion around 1995 presumably caused by a selective advantage. Resistance was, however, subsequently lost in vivo by nine strains isolated between 2008 and 2012. We describe the genetic mechanisms involved in this loss of resistance and the imperfect relationship between genotypic and phenotypic resistance.

Conclusions

The recent re-emergence of methicillin susceptibility in this epidemic lineage suggests a significant fitness cost of resistance and reduced selective advantage following the introduction in the mid-2000s of MRSA hospital control measures throughout the UK.

Robust Prediction of Resistance to Trimethoprim in Staphylococcus aureus

The rise of antibiotic resistance threatens modern medicine; to combat it new diagnostic methods are required. Sequencing the whole genome of a pathogen offers the potential to accurately determine which antibiotics will be effective to treat a patient. A key limitation of this approach is that it cannot classify rare or previously unseen mutations. Here we demonstrate that alchemical free energy methods, a well-established class of methods from computational chemistry, can successfully predict whether mutations in Staphylococcus aureus dihydrofolate reductase confer resistance to trimethoprim. We also show that the method is quantitatively accurate by calculating how much the most common resistance-conferring mutation, F99Y, reduces the binding free energy of trimethoprim and comparing predicted and experimentally measured minimum inhibitory concentrations for seven different mutations. Finally, by considering up to 32 free energy calculations for each mutation, we estimate its specificity and sensitivity.

Severe infections emerge from commensal bacteria by adaptive evolution

Bacteria responsible for the greatest global mortality colonize the human microbiota far more frequently than they cause severe infections. Whether mutation and selection among commensal bacteria are associated with infection is unknown. We investigated de novo mutation in 1163 Staphylococcus aureus genomes from 105 infected patients with nose colonization. We report that 72% of infections emerged from the nose, with infecting and nose-colonizing bacteria showing parallel adaptive differences. We found 2.8-to-3.6-fold adaptive enrichments of protein-altering variants in genes responding to rsp, which regulates surface antigens and toxin production; agr, which regulates quorum-sensing, toxin production and abscess formation; and host-derived antimicrobial peptides. Adaptive mutations in pathogenesis-associated genes were 3.1-fold enriched in infecting but not nose-colonizing bacteria. None of these signatures were observed in healthy carriers nor at the species-level, suggesting infection-associated, short-term, within-host selection pressures. Our results show that signatures of spontaneous adaptive evolution are specifically associated with infection, raising new possibilities for diagnosis and treatment.

Covert dissemination of carbapenemase-producing Klebsiella pneumoniae (KPC) in a successfully controlled outbreak: long- and short-read whole-genome sequencing demonstrate multiple genetic modes of transmission

BACKGROUND

Carbapenemase-producing Enterobacteriaceae (CPE), including KPC-producing Klebsiella pneumoniae (KPC-Kpn), are an increasing threat to patient safety.

OBJECTIVES

To use WGS to investigate the extent and complexity of carbapenemase gene dissemination in a controlled KPC outbreak.

MATERIALS AND METHODS

Enterobacteriaceae with reduced ertapenem susceptibility recovered from rectal screening swabs/clinical samples, during a 3 month KPC outbreak (2013-14), were investigated for carbapenemase production, antimicrobial susceptibility, variable-number-tandem-repeat profile and WGS [short-read (Illumina), long-read (MinION)]. Short-read sequences were used for MLST and plasmid/Tn4401 fingerprinting, and long-read sequence assemblies for plasmid identification. Phylogenetic analysis used IQTree followed by ClonalFrameML, and outbreak transmission dynamics were inferred using SCOTTI.

RESULTS

Twenty patients harboured KPC-positive isolates (6 infected, 14 colonized), and 23 distinct KPC-producing Enterobacteriaceae were identified. Four distinct KPC plasmids were characterized but of 20 KPC-Kpn (from six STs), 17 isolates shared a single pKpQIL-D2 KPC plasmid. All isolates had an identical transposon (Tn4401a), except one KPC-Kpn (ST661) with a single nucleotide variant. A sporadic case of KPC-Kpn (ST491) with Tn4401a-carrying pKpQIL-D2 plasmid was identified 10 months before the outbreak. This plasmid was later seen in two other species and other KPC-Kpn (ST14,ST661) including clonal spread of KPC-Kpn (ST661) from a symptomatic case to nine ward contacts.

CONCLUSIONS

WGS of outbreak KPC isolates demonstrated blaKPC dissemination via horizontal transposition (Tn4401a), plasmid spread (pKpQIL-D2) and clonal spread (K. pneumoniae ST661). Despite rapid outbreak control, considerable dissemination of blaKPC still occurred among K. pneumoniae and other Enterobacteriaceae, emphasizing its high transmission potential and the need for enhanced control efforts.

Epidemiology of Clostridium difficile in infants in Oxfordshire, UK: Risk factors for colonization and carriage, and genetic overlap with regional C. difficile infection strains

BACKGROUND

Approximately 30-40% of children <1 year of age are Clostridium difficile colonized, and may represent a reservoir for adult C. difficile infections (CDI). Risk factors for colonization with toxigenic versus non-toxigenic C. difficile strains and longitudinal acquisition dynamics in infants remain incompletely characterized.

METHODS

Predominantly healthy infants (≤2 years) were recruited in Oxfordshire, UK, and provided ≥1 fecal samples. Independent risk factors for toxigenic/non-toxigenic C. difficile colonization and acquisition were identified using multivariable regression. Infant C. difficile isolates were whole-genome sequenced to assay genetic diversity and prevalence of toxin-associated genes, and compared with sequenced strains from Oxfordshire CDI cases.

RESULTS

338/365 enrolled infants provided 1332 fecal samples, representing 158 C. difficile colonization or carriage episodes (107[68%] toxigenic). Initial colonization was associated with age, and reduced with breastfeeding but increased with pet dogs. Acquisition was associated with older age, Caesarean delivery, and diarrhea. Breastfeeding and pre-existing C. difficile colonization reduced acquisition risk. Overall 13% of CDI C. difficile strains were genetically related to infant strains. 29(18%) infant C. difficile sequences were consistent with recent direct/indirect transmission to/from Oxfordshire CDI cases (≤2 single nucleotide variants [SNVs]); 79(50%) shared a common origin with an Oxfordshire CDI case within the last ~5 years (0-10 SNVs). The hypervirulent, epidemic ST1/ribotype 027 remained notably absent in infants in this large study, as did other lineages such as STs 10/44 (ribotype 015); the most common strain in infants was ST2 (ribotype 020/014)(22%).

CONCLUSIONS

In predominantly healthy infants without significant healthcare exposure C. difficile colonization and acquisition reflect environmental exposures, with pet dogs identified as a novel risk factor. Genetic overlap between some infant strains and those isolated from CDI cases suggest common community reservoirs of these C. difficile lineages, contrasting with those lineages found only in CDI cases, and therefore more consistent with healthcare-associated spread.

Genomic epidemiology of global Klebsiella pneumoniae carbapenemase (KPC)-producing Escherichia coli

The dissemination of carbapenem resistance in Escherichia coli has major implications for the management of common infections. bla _KPC, encoding a transmissible carbapenemase (KPC), has historically largely been associated with Klebsiella pneumoniae, a predominant plasmid (pKpQIL), and a specific transposable element (Tn4401, ~10 kb). Here we characterize the genetic features of bla _KPC emergence in global E. coli, 2008-2013, using both long- and short-read whole-genome sequencing. Amongst 43/45 successfully sequenced bla _KPC-E. coli strains, we identified substantial strain diversity (n = 21 sequence types, 18% of annotated genes in the core genome); substantial plasmid diversity (≥9 replicon types); and substantial bla _KPC-associated, mobile genetic element (MGE) diversity (50% not within complete Tn4401 elements). We also found evidence of inter-species, regional and international plasmid spread. In several cases bla _KPC was found on high copy number, small Col-like plasmids, previously associated with horizontal transmission of resistance genes in the absence of antimicrobial selection pressures. E. coli is a common human pathogen, but also a commensal in multiple environmental and animal reservoirs, and easily transmissible. The association of bla _KPC with a range of MGEs previously linked to the successful spread of widely endemic resistance mechanisms (e.g. bla _TEM, bla _CTX-M) suggests that it may become similarly prevalent.

Mortality risks associated with emergency admissions during weekends and public holidays: an analysis of electronic health records

BACKGROUND

Weekend hospital admission is associated with increased mortality, but the contributions of varying illness severity and admission time to this weekend effect remain unexplored.

METHODS

We analysed unselected emergency admissions to four Oxford University National Health Service hospitals in the UK from Jan 1, 2006, to Dec 31, 2014. The primary outcome was death within 30 days of admission (in or out of hospital), analysed using Cox models measuring time from admission. The primary exposure was day of the week of admission. We adjusted for multiple confounders including demographics, comorbidities, and admission characteristics, incorporating non-linearity and interactions. Models then considered the effect of adjusting for 15 common haematology and biochemistry test results or proxies for hospital workload.

FINDINGS

257 596 individuals underwent 503 938 emergency admissions. 18 313 (4·7%) patients admitted as weekday energency admissions and 6070 (5·1%) patients admitted as weekend emergency admissions died within 30 days (p<0·0001). 9347 individuals underwent 9707 emergency admissions on public holidays. 559 (5·8%) died within 30 days (p<0·0001 vs weekday). 15 routine haematology and biochemistry test results were highly prognostic for mortality. In 271 465 (53·9%) admissions with complete data, adjustment for test results explained 33% (95% CI 21 to 70) of the excess mortality associated with emergency admission on Saturdays compared with Wednesdays, 52% (lower 95% CI 34) on Sundays, and 87% (lower 95% CI 45) on public holidays after adjustment for standard patient characteristics. Excess mortality was predominantly restricted to admissions between 1100 h and 1500 h (p_interaction=0·04). No hospital workload measure was independently associated with mortality (all p values >0·06).

INTERPRETATION

Adjustment for routine test results substantially reduced excess mortality associated with emergency admission at weekends and public holidays. Adjustment for patient-level factors not available in our study might further reduce the residual excess mortality, particularly as this clustered around midday at weekends. Hospital workload was not associated with mortality. Together, these findings suggest that the weekend effect arises from patient-level differences at admission rather than reduced hospital staffing or services.

FUNDING

NIHR Oxford Biomedical Research Centre.

WGS to predict antibiotic MICs for Neisseria gonorrhoeae

Background

Tracking the spread of antimicrobial-resistant Neisseria gonorrhoeae is a major priority for national surveillance programmes.

Objectives

We investigate whether WGS and simultaneous analysis of multiple resistance determinants can be used to predict antimicrobial susceptibilities to the level of MICs in N. gonorrhoeae.

Methods

WGS was used to identify previously reported potential resistance determinants in 681 N. gonorrhoeae isolates, from England, the USA and Canada, with phenotypes for cefixime, penicillin, azithromycin, ciprofloxacin and tetracycline determined as part of national surveillance programmes. Multivariate linear regression models were used to identify genetic predictors of MIC. Model performance was assessed using leave-one-out cross-validation.

Results

Overall 1785/3380 (53%) MIC values were predicted to the nearest doubling dilution and 3147 (93%) within ±1 doubling dilution and 3314 (98%) within ±2 doubling dilutions. MIC prediction performance was similar across the five antimicrobials tested. Prediction models included the majority of previously reported resistance determinants. Applying EUCAST breakpoints to MIC predictions, the overall very major error (VME; phenotypically resistant, WGS-prediction susceptible) rate was 21/1577 (1.3%, 95% CI 0.8%-2.0%) and the major error (ME; phenotypically susceptible, WGS-prediction resistant) rate was 20/1186 (1.7%, 1.0%-2.6%). VME rates met regulatory thresholds for all antimicrobials except cefixime and ME rates for all antimicrobials except tetracycline. Country of testing was a strongly significant predictor of MIC for all five antimicrobials.

Conclusions

We demonstrate a WGS-based MIC prediction approach that allows reliable MIC prediction for five gonorrhoea antimicrobials. Our approach should allow reasonably precise prediction of MICs for a range of bacterial species.

Resolving plasmid structures in Enterobacteriaceae using the MinION nanopore sequencer: assessment of MinION and MinION/Illumina hybrid data assembly approaches

This study aimed to assess the feasibility of using the Oxford Nanopore Technologies (ONT) MinION long-read sequencer in reconstructing fully closed plasmid sequences from eight Enterobacteriaceae isolates of six different species with plasmid populations of varying complexity. Species represented were Escherichia coli, Klebsiella pneumoniae, Citrobacter freundii, Enterobacter cloacae, Serratia marcescens and Klebsiella oxytoca, with plasmid populations ranging from 1–11 plasmids with sizes of 2–330 kb. Isolates were sequenced using Illumina (short-read) and ONT’s MinION (long-read) platforms, and compared with fully resolved PacBio (long-read) sequence assemblies for the same isolates. We compared the performance of different assembly approaches including SPAdes, plasmidSPAdes, hybridSPAdes, Canu, Canu+Pilon (canuPilon) and npScarf in recovering the plasmid structures of these isolates by comparing with the gold-standard PacBio reference sequences. Overall, canuPilon provided consistently good quality assemblies both in terms of assembly statistics (N50, number of contigs) and assembly accuracy [presence of single nucleotide polymorphisms (SNPs)/indels with respect to the reference sequence]. For plasmid reconstruction, Canu recovered 70 % of the plasmids in complete contigs, and combining three assembly approaches (Canu or canuPilon, hybridSPAdes and plasmidSPAdes) resulted in a total 78 % recovery rate for all the plasmids. The analysis demonstrated the potential of using MinION sequencing technology to resolve important plasmid structures in Enterobacteriaceae species independent of and in conjunction with Illumina sequencing data. A consensus assembly derived from several assembly approaches could present significant benefit in accurately resolving the greatest number of plasmid structures.

Enhanced Klebsiella pneumoniae Carbapenemase Expression from a Novel Tn4401 Deletion

The Klebsiella pneumoniae carbapenemase gene (blaKPC) is typically located within mobile transposon Tn4401 Enhanced KPC expression has been associated with deletions in the putative promoter region upstream of blaKPC Illumina sequences from blaKPC-positive clinical isolates from a single institution were mapped to a Tn4401b reference sequence, which carries no deletions. The novel isoform Tn4401h (188-bp deletion [between istB and blaKPC]) was present in 14% (39/281) of clinical isolates. MICs showed that Escherichia coli strains containing plasmids with Tn4401a and Tn4401h were more resistant to meropenem (≥16 and ≥16, respectively), ertapenem (≥8 and 4, respectively), and cefepime (≥64 and 4, respectively) than E. coli strains with Tn4401b (0.5, ≤0.5, and ≤1, respectively). Quantitative real-time PCR (qRT-PCR) demonstrated that Tn4401a had a 16-fold increase and Tn4401h a 4-fold increase in blaKPC mRNA levels compared to the reference Tn4401b. A lacZ reporter plasmid was used to test the activity of the promoter regions from the different variants, and the results showed that the Tn4401a and Tn4401h promoter sequences generated higher β-galactosidase activity than the corresponding Tn4401b sequence. Further dissection of the promoter region demonstrated that putative promoter P1 was not functional. The activity of the isolated P2 promoter was greatly enhanced by inclusion of the P1-P2 intervening sequence. These studies indicated that gene expression could be an important consideration in understanding resistance phenotypes predicted by genetic signatures in the context of sequencing-based rapid diagnostics.

Same-Day Diagnostic and Surveillance Data for Tuberculosis via Whole-Genome Sequencing of Direct Respiratory Samples

Routine full characterization of Mycobacterium tuberculosis is culture based, taking many weeks. Whole-genome sequencing (WGS) can generate antibiotic susceptibility profiles to inform treatment, augmented with strain information for global surveillance; such data could be transformative if provided at or near the point of care. We demonstrate a low-cost method of DNA extraction directly from patient samples for M. tuberculosis WGS. We initially evaluated the method by using the Illumina MiSeq sequencer (40 smear-positive respiratory samples obtained after routine clinical testing and 27 matched liquid cultures). M. tuberculosis was identified in all 39 samples from which DNA was successfully extracted. Sufficient data for antibiotic susceptibility prediction were obtained from 24 (62%) samples; all results were concordant with reference laboratory phenotypes. Phylogenetic placement was concordant between direct and cultured samples. With Illumina MiSeq/MiniSeq, the workflow from patient sample to results can be completed in 44/16 h at a reagent cost of £96/£198 per sample. We then employed a nonspecific PCR-based library preparation method for sequencing on an Oxford Nanopore Technologies MinION sequencer. We applied this to cultured Mycobacterium bovis strain BCG DNA and to combined culture-negative sputum DNA and BCG DNA. For flow cell version R9.4, the estimated turnaround time from patient to identification of BCG, detection of pyrazinamide resistance, and phylogenetic placement was 7.5 h, with full susceptibility results 5 h later. Antibiotic susceptibility predictions were fully concordant. A critical advantage of MinION is the ability to continue sequencing until sufficient coverage is obtained, providing a potential solution to the problem of variable amounts of M. tuberculosis DNA in direct samples.