David W. Eyre

Whole-genome sequencing to delineate Mycobacterium tuberculosis outbreaks: a retrospective observational study

Summary Background Tuberculosis incidence in the UK has risen in the past decade. Disease control depends on epidemiological data, which can be difficult to obtain. Whole-genome sequencing can detect microevolution within Mycobacterium tuberculosis strains. We aimed to estimate the genetic diversity of related M tuberculosis strains in the UK Midlands and to investigate how this measurement might be used to investigate community outbreaks. Methods In a retrospective observational study, we used Illumina technology to sequence M tuberculosis genomes from an archive of frozen cultures. We characterised isolates into four groups: cross-sectional, longitudinal, household, and community. We measured pairwise nucleotide differences within hosts and between hosts in household outbreaks and estimated the rate of change in DNA sequences. We used the findings to interpret network diagrams constructed from 11 community clusters derived from mycobacterial interspersed repetitive-unit–variable-number tandem-repeat data. Findings We sequenced 390 separate isolates from 254 patients, including representatives from all five major lineages of M tuberculosis. The estimated rate of change in DNA sequences was 0·5 single nucleotide polymorphisms (SNPs) per genome per year (95% CI 0·3–0·7) in longitudinal isolates from 30 individuals and 25 families. Divergence is rarely higher than five SNPs in 3 years. 109 (96%) of 114 paired isolates from individuals and households differed by five or fewer SNPs. More than five SNPs separated isolates from none of 69 epidemiologically linked patients, two (15%) of 13 possibly linked patients, and 13 (17%) of 75 epidemiologically unlinked patients (three-way comparison exact p<0·0001). Genetic trees and clinical and epidemiological data suggest that super-spreaders were present in two community clusters. Interpretation Whole-genome sequencing can delineate outbreaks of tuberculosis and allows inference about direction of transmission between cases. The technique could identify super-spreaders and predict the existence of undiagnosed cases, potentially leading to early treatment of infectious patients and their contacts. Funding Medical Research Council, Wellcome Trust, National Institute for Health Research, and the Health Protection Agency.

Diverse sources of C. difficile infection identified on whole-genome sequencing.

BACKGROUND It has been thought that Clostridium difficile infection is transmitted predominantly within health care settings. However, endemic spread has hampered identification of precise sources of infection and the assessment of the efficacy of interventions. METHODS From September 2007 through March 2011, we performed whole-genome sequencing on isolates obtained from all symptomatic patients with C. difficile infection identified in health care settings or in the community in Oxfordshire, United Kingdom. We compared single-nucleotide variants (SNVs) between the isolates, using C. difficile evolution rates estimated on the basis of the first and last samples obtained from each of 145 patients, with 0 to 2 SNVs expected between transmitted isolates obtained less than 124 days apart, on the basis of a 95% prediction interval. We then identified plausible epidemiologic links among genetically related cases from data on hospital admissions and community location. RESULTS Of 1250 C. difficile cases that were evaluated, 1223 (98%) were successfully sequenced. In a comparison of 957 samples obtained from April 2008 through March 2011 with those obtained from September 2007 onward, a total of 333 isolates (35%) had no more than 2 SNVs from at least 1 earlier case, and 428 isolates (45%) had more than 10 SNVs from all previous cases. Reductions in incidence over time were similar in the two groups, a finding that suggests an effect of interventions targeting the transition from exposure to disease. Of the 333 patients with no more than 2 SNVs (consistent with transmission), 126 patients (38%) had close hospital contact with another patient, and 120 patients (36%) had no hospital or community contact with another patient. Distinct subtypes of infection continued to be identified throughout the study, which suggests a considerable reservoir of C. difficile. CONCLUSIONS Over a 3-year period, 45% of C. difficile cases in Oxfordshire were genetically distinct from all previous cases. Genetically diverse sources, in addition to symptomatic patients, play a major part in C. difficile transmission. (Funded by the U.K. Clinical Research Collaboration Translational Infection Research Initiative and others.).

A pilot study of rapid benchtop sequencing of Staphylococcus aureus and Clostridium difficile for outbreak detection and surveillance

Objectives To investigate the prospects of newly available benchtop sequencers to provide rapid whole-genome data in routine clinical practice. Next-generation sequencing has the potential to resolve uncertainties surrounding the route and timing of person-to-person transmission of healthcare-associated infection, which has been a major impediment to optimal management. Design The authors used Illumina MiSeq benchtop sequencing to undertake case studies investigating potential outbreaks of methicillin-resistant Staphylococcus aureus (MRSA) and Clostridium difficile. Setting Isolates were obtained from potential outbreaks associated with three UK hospitals. Participants Isolates were sequenced from a cluster of eight MRSA carriers and an associated bacteraemia case in an intensive care unit, another MRSA cluster of six cases and two clusters of C difficile. Additionally, all C difficile isolates from cases over 6 weeks in a single hospital were rapidly sequenced and compared with local strain sequences obtained in the preceding 3 years. Main outcome measure Whole-genome genetic relatedness of the isolates within each epidemiological cluster. Results Twenty-six MRSA and 15 C difficile isolates were successfully sequenced and analysed within 5 days of culture. Both MRSA clusters were identified as outbreaks, with most sequences in each cluster indistinguishable and all within three single nucleotide variants (SNVs). Epidemiologically unrelated isolates of the same spa-type were genetically distinct (≥21 SNVs). In both C difficile clusters, closely epidemiologically linked cases (in one case sharing the same strain type) were shown to be genetically distinct (≥144 SNVs). A reconstruction applying rapid sequencing in C difficile surveillance provided early outbreak detection and identified previously undetected probable community transmission. Conclusions This benchtop sequencing technology is widely generalisable to human bacterial pathogens. The findings provide several good examples of how rapid and precise sequencing could transform identification of transmission of healthcare-associated infection and therefore improve hospital infection control and patient outcomes in routine clinical practice.

Relationship between bacterial strain type, host biomarkers, and mortality in Clostridium difficile infection.

Clostridium difficile genotype predicts 14-day mortality in 1893 enzyme immunoassay–positive/culture-positive adults. Excess mortality correlates with genotype-specific changes in biomarkers, strongly implicating inflammatory pathways as a major influence on poor outcome. Polymerase chain reaction ribotype 078/ST 11(clade 5) is associated with high mortality; ongoing surveillance remains essential.

Predictors of first recurrence of Clostridium difficile infection: implications for initial management.

Symptomatic recurrence of Clostridium difficile infection (CDI) occurs in approximately 20% of patients and is challenging to treat. Identifying those at high risk could allow targeted initial management and improve outcomes. Adult toxin enzyme immunoassay–positive CDI cases in a population of approximately 600 000 persons from September 2006 through December 2010 were combined with epidemiological/clinical data. The cumulative incidence of recurrence ≥14 days after the diagnosis and/or onset of first-ever CDI was estimated, treating death without recurrence as a competing risk, and predictors were identified from cause-specific proportional hazards regression models. A total of 1678 adults alive 14 days after their first CDI were included; median age was 77 years, and 1191 (78%) were inpatients. Of these, 363 (22%) experienced a recurrence ≥14 days after their first CDI, and 594 (35%) died without recurrence through March 2011. Recurrence risk was independently and significantly higher among patients admitted as emergencies, with previous gastrointestinal ward admission(s), last discharged 4–12 weeks before first diagnosis, and with CDI diagnosed at admission. Recurrence risk also increased with increasing age, previous total hours admitted, and C-reactive protein level at first CDI (all P < .05). The 4-month recurrence risk increased by approximately 5% (absolute) for every 1-point increase in a risk score based on these factors. Risk factors, including increasing age, initial disease severity, and hospital exposure, predict CDI recurrence and identify patients likely to benefit from enhanced initial CDI treatment.

Predicting antimicrobial susceptibilities for Escherichia coli and Klebsiella pneumoniae isolates using whole genomic sequence data

Objectives Whole-genome sequencing potentially represents a single, rapid and cost-effective approach to defining resistance mechanisms and predicting phenotype, and strain type, for both clinical and epidemiological purposes. This retrospective study aimed to determine the efficacy of whole genome-based antimicrobial resistance prediction in clinical isolates of Escherichia coli and Klebsiella pneumoniae. Methods Seventy-four E. coli and 69 K. pneumoniae bacteraemia isolates from Oxfordshire, UK, were sequenced (Illumina HiSeq 2000). Resistance phenotypes were predicted from genomic sequences using BLASTn-based comparisons of de novo-assembled contigs with a study database of >100 known resistance-associated loci, including plasmid-associated and chromosomal genes. Predictions were made for seven commonly used antimicrobials: amoxicillin, co-amoxiclav, ceftriaxone, ceftazidime, ciprofloxacin, gentamicin and meropenem. Comparisons were made with phenotypic results obtained in duplicate by broth dilution (BD Phoenix). Discrepancies, either between duplicate BD Phoenix results or between genotype and phenotype, were resolved with gradient diffusion analyses. Results A wide variety of antimicrobial resistance genes were identified, including blaCTX-M, blaLEN, blaOKP, blaOXA, blaSHV, blaTEM, aac(3′)-Ia, aac-(3′)-IId, aac-(3′)-IIe, aac(6′)-Ib-cr, aadA1a, aadA4, aadA5, aadA16, aph(6′)-Id, aph(3′)-Ia, qnrB and qnrS, as well as resistance-associated mutations in chromosomal gyrA and parC genes. The sensitivity of genome-based resistance prediction across all antibiotics for both species was 0.96 (95% CI: 0.94–0.98) and the specificity was 0.97 (95% CI: 0.95–0.98). Very major and major error rates were 1.2% and 2.1%, respectively. Conclusions Our method was as sensitive and specific as routinely deployed phenotypic methods. Validation against larger datasets and formal assessments of cost and turnaround time in a routine laboratory setting are warranted.

Microevolutionary analysis of Clostridium difficile genomes to investigate transmission

BackgroundThe control of Clostridium difficile infection is a major international healthcare priority, hindered by a limited understanding of transmission epidemiology for these bacteria. However, transmission studies of bacterial pathogens are rapidly being transformed by the advent of next generation sequencing.ResultsHere we sequence whole C. difficile genomes from 486 cases arising over four years in Oxfordshire. We show that we can estimate the times back to common ancestors of bacterial lineages with sufficient resolution to distinguish whether direct transmission is plausible or not. Time depths were inferred using a within-host evolutionary rate that we estimated at 1.4 mutations per genome per year based on serially isolated genomes. The subset of plausible transmissions was found to be highly associated with pairs of patients sharing time and space in hospital. Conversely, the large majority of pairs of genomes matched by conventional typing and isolated from patients within a month of each other were too distantly related to be direct transmissions.ConclusionsOur results confirm that nosocomial transmission between symptomatic C. difficile cases contributes far less to current rates of infection than has been widely assumed, which clarifies the importance of future research into other transmission routes, such as from asymptomatic carriers. With the costs of DNA sequencing rapidly falling and its use becoming more and more widespread, genomics will revolutionize our understanding of the transmission of bacterial pathogens.

Evolutionary History of the Clostridium difficile Pathogenicity Locus

The symptoms of Clostridium difficile infection are caused by toxins expressed from its 19 kb pathogenicity locus (PaLoc). Stable integration of the PaLoc is suggested by its single chromosomal location and the clade specificity of its different genetic variants. However, the PaLoc is variably present, even among closely related strains, and thus resembles a mobile genetic element. Our aim was to explain these apparently conflicting observations by reconstructing the evolutionary history of the PaLoc. Phylogenetic analyses and annotation of the regions spanning the PaLoc were performed using C. difficile population-representative genomes chosen from a collection of 1,693 toxigenic (PaLoc present) and nontoxigenic (PaLoc absent) isolates. Comparison of the core genome and PaLoc phylogenies demonstrated an eventful evolutionary history, with distinct PaLoc variants acquired clade specifically after divergence. In particular, our data suggest a relatively recent PaLoc acquisition in clade 4. Exchanges and losses of the PaLoc DNA have also occurred, via long homologous recombination events involving flanking chromosomal sequences. The most recent loss event occurred ∼30 years ago within a clade 1 genotype. The genetic organization of the clade 3 PaLoc was unique in containing a stably integrated novel transposon (designated Tn6218), variants of which were found at multiple chromosomal locations. Tn6218 elements were Tn916-related but nonconjugative and occasionally contained genes conferring resistance to clinically relevant antibiotics. The evolutionary histories of two contrasting but clinically important genetic elements were thus characterized: the PaLoc, mobilized rarely via homologous recombination, and Tn6218, mobilized frequently through transposition.

Effects of control interventions on Clostridium difficile infection in England: an observational study

Summary Background The control of Clostridium difficile infections is an international clinical challenge. The incidence of C difficile in England declined by roughly 80% after 2006, following the implementation of national control policies; we tested two hypotheses to investigate their role in this decline. First, if C difficile infection declines in England were driven by reductions in use of particular antibiotics, then incidence of C difficile infections caused by resistant isolates should decline faster than that caused by susceptible isolates across multiple genotypes. Second, if C difficile infection declines were driven by improvements in hospital infection control, then transmitted (secondary) cases should decline regardless of susceptibility. Methods Regional (Oxfordshire and Leeds, UK) and national data for the incidence of C difficile infections and antimicrobial prescribing data (1998–2014) were combined with whole genome sequences from 4045 national and international C difficile isolates. Genotype (multilocus sequence type) and fluoroquinolone susceptibility were determined from whole genome sequences. The incidence of C difficile infections caused by fluoroquinolone-resistant and fluoroquinolone-susceptible isolates was estimated with negative-binomial regression, overall and per genotype. Selection and transmission were investigated with phylogenetic analyses. Findings National fluoroquinolone and cephalosporin prescribing correlated highly with incidence of C difficile infections (cross-correlations >0·88), by contrast with total antibiotic prescribing (cross-correlations <0·59). Regionally, C difficile decline was driven by elimination of fluoroquinolone-resistant isolates (approximately 67% of Oxfordshire infections in September, 2006, falling to approximately 3% in February, 2013; annual incidence rate ratio 0·52, 95% CI 0·48–0·56 vs fluoroquinolone-susceptible isolates: 1·02, 0·97–1·08). C difficile infections caused by fluoroquinolone-resistant isolates declined in four distinct genotypes (p<0·01). The regions of phylogenies containing fluoroquinolone-resistant isolates were short-branched and geographically structured, consistent with selection and rapid transmission. The importance of fluoroquinolone restriction over infection control was shown by significant declines in inferred secondary (transmitted) cases caused by fluoroquinolone-resistant isolates with or without hospital contact (p<0·0001) versus no change in either group of cases caused by fluoroquinolone-susceptible isolates (p>0·2). Interpretation Restricting fluoroquinolone prescribing appears to explain the decline in incidence of C difficile infections, above other measures, in Oxfordshire and Leeds, England. Antimicrobial stewardship should be a central component of C difficile infection control programmes. Funding UK Clinical Research Collaboration (Medical Research Council, Wellcome Trust, National Institute for Health Research); NIHR Oxford Biomedical Research Centre; NIHR Health Protection Research Unit on Healthcare Associated Infection and Antimicrobial Resistance (Oxford University in partnership with Public Health England [PHE]), and on Modelling Methodology (Imperial College, London in partnership with PHE); and the Health Innovation Challenge Fund.

Detection of mixed infection from bacterial whole genome sequence data allows assessment of its role in Clostridium difficile transmission.

Bacterial whole genome sequencing offers the prospect of rapid and high precision investigation of infectious disease outbreaks. Close genetic relationships between microorganisms isolated from different infected cases suggest transmission is a strong possibility, whereas transmission between cases with genetically distinct bacterial isolates can be excluded. However, undetected mixed infections—infection with ≥2 unrelated strains of the same species where only one is sequenced—potentially impairs exclusion of transmission with certainty, and may therefore limit the utility of this technique. We investigated the problem by developing a computationally efficient method for detecting mixed infection without the need for resource-intensive independent sequencing of multiple bacterial colonies. Given the relatively low density of single nucleotide polymorphisms within bacterial sequence data, direct reconstruction of mixed infection haplotypes from current short-read sequence data is not consistently possible. We therefore use a two-step maximum likelihood-based approach, assuming each sample contains up to two infecting strains. We jointly estimate the proportion of the infection arising from the dominant and minor strains, and the sequence divergence between these strains. In cases where mixed infection is confirmed, the dominant and minor haplotypes are then matched to a database of previously sequenced local isolates. We demonstrate the performance of our algorithm with in silico and in vitro mixed infection experiments, and apply it to transmission of an important healthcare-associated pathogen, Clostridium difficile. Using hospital ward movement data in a previously described stochastic transmission model, 15 pairs of cases enriched for likely transmission events associated with mixed infection were selected. Our method identified four previously undetected mixed infections, and a previously undetected transmission event, but no direct transmission between the pairs of cases under investigation. These results demonstrate that mixed infections can be detected without additional sequencing effort, and this will be important in assessing the extent of cryptic transmission in our hospitals.