Whole genome sequencing for M/XDR tuberculosis surveillance and for resistance testing

Whole genome sequencing (WGS) can help to relate Mycobacterium tuberculosis genomes to one another to assess genetic relatedness and infer the likelihood of transmission between cases. The same sequence data are now increasingly being used to predict drug resistance and susceptibility. Controlling the spread of tuberculosis and providing patients with the correct treatment are central to the World Health Organization's target to 'End TB' by 2035, for which the global prevalence of drug-resistant tuberculosis remains one of the main obstacles to success. So far, WGS has been applied largely to drug-susceptible strains for the purposes of understanding transmission, leaving a number of analytical considerations before transferring what has been learnt from drug-susceptible disease to drug-resistant tuberculosis. We discuss these potential problems here, alongside some of the challenges to characterizing the Mycobacterium tuberculosis 'resistome'-the optimal knowledge-base required for WGS-based assays to successfully direct individualized treatment regimens through the prediction of drug resistance and susceptibility in the future.

Natural mutations in a Staphylococcus aureus virulence regulator attenuate cytotoxicity but permit bacteremia and abscess formation

Staphylococcus aureus is a major bacterial pathogen, which causes severe blood and tissue infections that frequently emerge by autoinfection with asymptomatically carried nose and skin populations. However, recent studies report that bloodstream isolates differ systematically from those found in the nose and skin, exhibiting reduced toxicity toward leukocytes. In two patients, an attenuated toxicity bloodstream infection evolved from an asymptomatically carried high-toxicity nasal strain by loss-of-function mutations in the gene encoding the transcription factor repressor of surface proteins (rsp). Here, we report that rsp knockout mutants lead to global transcriptional and proteomic reprofiling, and they exhibit the greatest signal in a genome-wide screen for genes influencing S. aureus survival in human cells. This effect is likely to be mediated in part via SSR42, a long-noncoding RNA. We show that rsp controls SSR42 expression, is induced by hydrogen peroxide, and is required for normal cytotoxicity and hemolytic activity. Rsp inactivation in laboratory- and bacteremia-derived mutants attenuates toxin production, but up-regulates other immune subversion proteins and reduces lethality during experimental infection. Crucially, inactivation of rsp preserves bacterial dissemination, because it affects neither formation of deep abscesses in mice nor survival in human blood. Thus, we have identified a spontaneously evolving, attenuated-cytotoxicity, nonhemolytic S. aureus phenotype, controlled by a pleiotropic transcriptional regulator/noncoding RNA virulence regulatory system, capable of causing S. aureus bloodstream infections. Such a phenotype could promote deep infection with limited early clinical manifestations, raising concerns that bacterial evolution within the human body may contribute to severe infection.

The Stealthy Superbug: the Role of Asymptomatic Enteric Carriage in Maintaining a Long-Term Hospital Outbreak of ST228 Methicillin-Resistant Staphylococcus aureus

Whole-genome sequencing (WGS) of 228 isolates was used to elucidate the origin and dynamics of a long-term outbreak of methicillin-resistant Staphylococcus aureus (MRSA) sequence type 228 (ST228) SCCmec I that involved 1,600 patients in a tertiary care hospital between 2008 and 2012. Combining of the sequence data with detailed metadata on patient admission and movement confirmed that the outbreak was due to the transmission of a single clonal variant of ST228, rather than repeated introductions of this clone into the hospital. We note that this clone is significantly more frequently recovered from groin and rectal swabs than other clones (P < 0.0001) and is also significantly more transmissible between roommates (P < 0.01). Unrecognized MRSA carriers, together with movements of patients within the hospital, also seem to have played a major role. These atypical colonization and transmission dynamics can help explain how the outbreak was maintained over the long term. This "stealthy" asymptomatic colonization of the gut, combined with heightened transmissibility (potentially reflecting a role for environmental reservoirs), means the dynamics of this outbreak share some properties with enteric pathogens such as vancomycin-resistant enterococci or Clostridium difficile.

IMPORTANCE

Using whole-genome sequencing, we showed that a large and prolonged outbreak of methicillin-resistant Staphylococcus aureus was due to the clonal spread of a specific strain with genetic elements adapted to the hospital environment. Unrecognized MRSA carriers, the movement of patients within the hospital, and the low detection with clinical specimens were also factors that played a role in this occurrence. The atypical colonization of the gut means the dynamics of this outbreak may share some properties with enteric pathogens.

Rapid antibiotic-resistance predictions from genome sequence data for Staphylococcus aureus and Mycobacterium tuberculosis

The rise of antibiotic-resistant bacteria has led to an urgent need for rapid detection of drug resistance in clinical samples, and improvements in global surveillance. Here we show how de Bruijn graph representation of bacterial diversity can be used to identify species and resistance profiles of clinical isolates. We implement this method for Staphylococcus aureus and Mycobacterium tuberculosis in a software package ('Mykrobe predictor') that takes raw sequence data as input, and generates a clinician-friendly report within 3 minutes on a laptop. For S. aureus, the error rates of our method are comparable to gold-standard phenotypic methods, with sensitivity/specificity of 99.1%/99.6% across 12 antibiotics (using an independent validation set, n=470). For M. tuberculosis, our method predicts resistance with sensitivity/specificity of 82.6%/98.5% (independent validation set, n=1,609); sensitivity is lower here, probably because of limited understanding of the underlying genetic mechanisms. We give evidence that minor alleles improve detection of extremely drug-resistant strains, and demonstrate feasibility of the use of emerging single-molecule nanopore sequencing techniques for these purposes.

Evolutionary Trade-Offs Underlie the Multi-faceted Virulence of Staphylococcus aureus

Bacterial virulence is a multifaceted trait where the interactions between pathogen and host factors affect the severity and outcome of the infection. Toxin secretion is central to the biology of many bacterial pathogens and is widely accepted as playing a crucial role in disease pathology. To understand the relationship between toxicity and bacterial virulence in greater depth, we studied two sequenced collections of the major human pathogen Staphylococcus aureus and found an unexpected inverse correlation between bacterial toxicity and disease severity. By applying a functional genomics approach, we identified several novel toxicity-affecting loci responsible for the wide range in toxic phenotypes observed within these collections. To understand the apparent higher propensity of low toxicity isolates to cause bacteraemia, we performed several functional assays, and our findings suggest that within-host fitness differences between high- and low-toxicity isolates in human serum is a contributing factor. As invasive infections, such as bacteraemia, limit the opportunities for onward transmission, highly toxic strains could gain an additional between-host fitness advantage, potentially contributing to the maintenance of toxicity at the population level. Our results clearly demonstrate how evolutionary trade-offs between toxicity, relative fitness, and transmissibility are critical for understanding the multifaceted nature of bacterial virulence.

This study shows that, contrary to expectation, toxin secretion inversely correlates with disease severity for the major human pathogen Staphylococcus aureus.

Global efforts to counter the growing problem of antibiotic resistance and develop alternative treatment strategies rely on a fuller understanding of when and why opportunistic pathogens cause disease. Recent advances in DNA sequencing technologies have opened up new opportunities to study infectious organisms, yet identifying the genetic variants that explain differences in disease remains challenging. Here we aimed to understand the complex relationship between toxicity—a known risk factor for disease in many bacteria—and infection severity for the major human pathogen S. aureus. Against expectations, we found that the bacteria that caused the most severe disease were the least toxic strains. We were able to determine the mutations responsible for the differences in toxicity and identified a number of novel toxicity-affecting genes. We further discovered that bacterial fitness in human serum could explain the unexpected association of low-toxicity isolates with severe, invasive disease. Invasive S. aureus infections are usually considered a dead end for these bacteria, as these infections are rarely transmitted to another person. Here we show using a simple mathematical model that this might in fact favour transmission of highly toxic bacteria on a population level and thus contribute to their global success. Our work therefore highlights the complexity of bacterial infection and should aid in devising new treatment and control strategies against this important pathogen.

Extensive Within-Host Diversity in Fecally Carried Extended-Spectrum-Beta-Lactamase-Producing Escherichia coli Isolates: Implications for Transmission Analyses

Studies of the transmission epidemiology of antimicrobial-resistant Escherichia coli, such as strains harboring extended-spectrum beta-lactamase (ESBL) genes, frequently use selective culture of rectal surveillance swabs to identify isolates for molecular epidemiological investigation. Typically, only single colonies are evaluated, which risks underestimating species diversity and transmission events. We sequenced the genomes of 16 E. coli colonies from each of eight fecal samples (n = 127 genomes; one failure), taken from different individuals in Cambodia, a region of high ESBL-producing E. coli prevalence. Sequence data were used to characterize both the core chromosomal diversity of E. coli isolates and their resistance/virulence gene content as a proxy measure of accessory genome diversity. The 127 E. coli genomes represented 31 distinct sequence types (STs). Seven (88%) of eight subjects carried ESBL-positive isolates, all containing blaCTX-M variants. Diversity was substantial, with a median of four STs/individual (range, 1 to 10) and wide genetic divergence at the nucleotide level within some STs. In 2/8 (25%) individuals, the same blaCTX-M variant occurred in different clones, and/or different blaCTX-M variants occurred in the same clone. Patterns of other resistance genes and common virulence factors, representing differences in the accessory genome, were also diverse within and between clones. The substantial diversity among intestinally carried ESBL-positive E. coli bacteria suggests that fecal surveillance, particularly if based on single-colony subcultures, will likely underestimate transmission events, especially in high-prevalence settings.

Whole-genome sequencing for prediction of Mycobacterium tuberculosis drug susceptibility and resistance: a retrospective cohort study

BACKGROUND

Diagnosing drug-resistance remains an obstacle to the elimination of tuberculosis. Phenotypic drug-susceptibility testing is slow and expensive, and commercial genotypic assays screen only common resistance-determining mutations. We used whole-genome sequencing to characterise common and rare mutations predicting drug resistance, or consistency with susceptibility, for all first-line and second-line drugs for tuberculosis.

METHODS

Between Sept 1, 2010, and Dec 1, 2013, we sequenced a training set of 2099 Mycobacterium tuberculosis genomes. For 23 candidate genes identified from the drug-resistance scientific literature, we algorithmically characterised genetic mutations as not conferring resistance (benign), resistance determinants, or uncharacterised. We then assessed the ability of these characterisations to predict phenotypic drug-susceptibility testing for an independent validation set of 1552 genomes. We sought mutations under similar selection pressure to those characterised as resistance determinants outside candidate genes to account for residual phenotypic resistance.

FINDINGS

We characterised 120 training-set mutations as resistance determining, and 772 as benign. With these mutations, we could predict 89·2% of the validation-set phenotypes with a mean 92·3% sensitivity (95% CI 90·7-93·7) and 98·4% specificity (98·1-98·7). 10·8% of validation-set phenotypes could not be predicted because uncharacterised mutations were present. With an in-silico comparison, characterised resistance determinants had higher sensitivity than the mutations from three line-probe assays (85·1% vs 81·6%). No additional resistance determinants were identified among mutations under selection pressure in non-candidate genes.

INTERPRETATION

A broad catalogue of genetic mutations enable data from whole-genome sequencing to be used clinically to predict drug resistance, drug susceptibility, or to identify drug phenotypes that cannot yet be genetically predicted. This approach could be integrated into routine diagnostic workflows, phasing out phenotypic drug-susceptibility testing while reporting drug resistance early.

FUNDING

Wellcome Trust, National Institute of Health Research, Medical Research Council, and the European Union.

Can rapid integrated polymerase chain reaction-based diagnostics for gastrointestinal pathogens improve routine hospital infection control practice? A diagnostic study

BACKGROUND

Every year approximately 5000-9000 patients are admitted to a hospital with diarrhoea, which in up to 90% of cases has a non-infectious cause. As a result, single rooms are 'blocked' by patients with non-infectious diarrhoea, while patients with infectious diarrhoea are still in open bays because of a lack of free side rooms. A rapid test for differentiating infectious from non-infectious diarrhoea could be very beneficial for patients.

OBJECTIVE

To evaluate MassCode multiplex polymerase chain reaction (PCR) for the simultaneous diagnosis of multiple enteropathogens directly from stool, in terms of sensitivity/specificity to detect four common important enteropathogens: Clostridium difficile, Campylobacter spp., Salmonella spp. and norovirus.

DESIGN

A retrospective study of fixed numbers of samples positive for C. difficile (n = 200), Campylobacter spp. (n = 200), Salmonella spp. (n = 100) and norovirus (n = 200) plus samples negative for all these pathogens (n = 300). Samples were sourced from NHS microbiology laboratories in Oxford and Leeds where initial diagnostic testing was performed according to Public Health England methodology. Researchers carrying out MassCode assays were blind to this information. A questionnaire survey, examining current practice for infection control teams and microbiology laboratories managing infectious diarrhoea, was also carried out.

SETTING

MassCode assays were carried out at Oxford University Hospitals NHS Trust. Further multiplex assays, carried out using Luminex, were run on the same set of samples at Leeds Teaching Hospitals NHS Trust. The questionnaire was completed by various NHS trusts.

MAIN OUTCOME MEASURES

Sensitivity and specificity to detect C. difficile, Campylobacter spp., Salmonella spp., and norovirus.

RESULTS

Nucleic acids were extracted from 948 clinical samples using an optimised protocol (200 Campylobacter spp., 199 C. difficile, 60 S. enterica, 199 norovirus and 295 negative samples; some samples contained more than one pathogen). Using the MassCode assay, sensitivities for each organism compared with standard microbiological testing ranged from 43% to 94% and specificities from 95% to 98%, with particularly poor performance for S. enterica. Relatively large numbers of unexpected positives not confirmed with quantitative PCR were also observed, particularly for S. enterica, Giardia lamblia and Cryptosporidium spp. As the results indicated that S. enterica detection might provide generic challenges to other multiplex assays for gastrointestinal pathogens, the Luminex xTag(®) gastrointestinal assay was also run blinded on the same extracts (937/948 remaining) and on re-extracted samples (839/948 with sufficient material). For Campylobacter spp., C. difficile and norovirus, high sensitivities (> 92%) and specificities (> 96%) were observed. For S. enterica, on the original MassCode/Oxford extracts, Luminex sensitivity compared with standard microbiological testing was 84% [95% confidence interval (CI) 73% to 93%], but this dropped to 46% on a fresh extract, very similar to MassCode, with a corresponding increase in specificity from 92% to 99%. Overall agreement on the per-sample diagnosis compared with combined microbiology plus PCR for the main four/all pathogens was 85.6%/64.7%, 87.0%/82.9% and 89.8%/86.8% for the MassCode assay, Luminex assay/MassCode extract and Luminex assay/fresh extract, respectively. Luminex assay results from fresh extracts implied that 5% of samples did not represent infectious diarrhoea, even though enteropathogens were genuinely present. Managing infectious diarrhoea was a significant burden for infection control teams (taking 21% of their time) and better diagnostics were identified as having major potential benefits for patients.

CONCLUSIONS

Overall, the Luminex xTag gastrointestinal panel showed similar or superior sensitivity and specificity to the MassCode assay. However, on fresh extracts, this test had low sensitivity to detect a key enteric pathogen, S. enterica; making it an unrealistic option for most microbiology laboratories. Extraction efficiency appears to be a major obstacle for nucleic acid-based tests for this organism, and possibly the whole Enterobacteriaceae family. To improve workflows in service microbiology laboratories, to reduce workload for infection control practitioners, and to improve outcomes for NHS patients, further research on deoxyribonucleic acid-based multiplex gastrointestinal diagnostics is urgently needed.

FUNDING

The Health Technology Assessment programme of the National Institute for Health Research.

Monitoring meticillin resistant Staphylococcus aureus and its spread in Copenhagen, Denmark, 2013, through routine whole genome sequencing

Typing of meticillin resistant Staphylococcus aureus (MRSA) by whole genome sequencing (WGS) is performed routinely in Copenhagen since January 2013. We describe the relatedness, based on WGS data and epidemiological data, of 341 MRSA isolates. These comprised all MRSA (n = 300) identified in Copenhagen in the first five months of 2013. Moreover, because MRSA of staphylococcal protein A (spa)-type 304 (t304), sequence type (ST) 6 had been associated with a continuous neonatal ward outbreak in Copenhagen starting in 2011, 41 t304 isolates collected in the city between 2010 and 2012 were also included. Isolates from 2013 found to be of t304, ST6 (n=14) were compared to the 41 earlier isolates. In the study, isolates of clonal complex (CC) 22 were examined in detail, as this CC has been shown to include the hospital-acquired epidemic MRSA (EMRSA-15) clone. Finally, all MRSA ST80 were also further analysed, as representatives of an important community-acquired MRSA in Europe. Overall the analysis identified 85 spa-types and 35 STs from 17 CCs. WGS confirmed the relatedness of epidemiologically linked t304 neonatal outbreak isolates. Several non-outbreak related patients had isolates closely related to the neonatal isolates suggesting unrecognised community chains of transmission and insufficient epidemiological data. Only four CC22 isolates were related to EMRSA-15. No community spread was observed among the 13 ST80 isolates. WGS successfully replaced conventional typing and added information to epidemiological surveillance. Creation of a MRSA database allows clustering of isolates based on single nucleotide polymorphism (SNP) calling and has improved our understanding of MRSA transmission.

Rare variants in MYD88, IRAK4 and IKBKG and susceptibility to invasive pneumococcal disease: a population-based case-control study

Although rare variants within the Toll-like receptor signalling pathway genes have been found to underlie human primary immunodeficiencies associated with selective predisposition to invasive pneumococcal disease (IPD), the contribution of variants in these genes to IPD susceptibility at the population level remains unknown. Complete re-sequencing of IRAK4, MYD88 and IKBKG genes was undertaken in 164 IPD cases from the UK and 164 geographically-matched population-based controls. 233 single-nucleotide variants (SNVs) were identified, of which ten were in coding regions. Four rare coding variants were predicted to be deleterious, two variants in MYD88 and two in IRAK4. The predicted deleterious variants in MYD88 were observed as two heterozygote cases but not seen in controls. Frequencies of predicted deleterious IRAK4 SNVs were the same in cases and controls. Our findings suggest that rare, functional variants in MYD88, IRAK4 or IKBKG do not significantly contribute to IPD susceptibility in adults at the population level.

Dynamics of MDR Enterobacter cloacae outbreaks in a neonatal unit in Nepal: insights using wider sampling frames and next-generation sequencing

OBJECTIVES

There are limited data on Enterobacter cloacae outbreaks and fewer describing these in association with NDM-1. With whole-genome sequencing, we tested the hypothesis that a cluster of 16 E. cloacae bacteraemia cases in a Nepali neonatal unit represented a single clonal outbreak, using a wider set of epidemiologically unrelated clinical E. cloacae isolates for comparison.

METHODS

Forty-three isolates were analysed, including 23 E. cloacae and 3 Citrobacter sp. isolates obtained from blood cultures from 16 neonates over a 3 month period. These were compared with two contemporaneous community-associated drug-resistant isolates from adults, a unit soap dispenser isolate and a set of historical invasive isolates (n=14) from the same geographical locality.

RESULTS

There were two clear neonatal outbreaks and one isolated case in the unit. One outbreak was associated with an NDM-1 plasmid also identified in a historical community-associated strain. The smaller, second outbreak was likely associated with a contaminated soap dispenser. The two community-acquired adult cases and three sets of historical hospital-associated neonatal isolates represented four additional genetic clusters.

CONCLUSIONS

E. cloacae infections in this context represent several different transmission networks, operating at the community/hospital and host strain/plasmid levels. Wide sampling frames and high-resolution typing methods are needed to describe the complex molecular epidemiology of E. cloacae outbreaks, which is not appropriately reflected by routine susceptibility phenotypes. Soap dispensers may represent a reservoir for E. cloacae and bacterial strains and plasmids may persist in hospitals and in the community for long periods, sporadically being involved in outbreaks of disease.

Emergence and spread of predominantly community-onset Clostridium difficile PCR ribotype 244 infection in Australia, 2010 to 2012

We describe an Australia-wide Clostridium difficile outbreak in 2011 and 2012 involving the previously uncommon ribotype 244. In Western Australia, 14 of 25 cases were community-associated, 11 were detected in patients younger than 65 years, 14 presented to emergency/outpatient departments, and 14 to non-tertiary/community hospitals. Using whole genome sequencing, we confirm ribotype 244 is from the same C. difficile clade as the epidemic ribotype 027. Like ribotype 027, it produces toxins A, B, and binary toxin, however it is fluoroquinolone-susceptible and thousands of single nucleotide variants distinct from ribotype 027. Fifteen outbreak isolates from across Australia were sequenced. Despite their geographic separation, all were genetically highly related without evidence of geographic clustering, consistent with a point source, for example affecting the national food chain. Comparison with reference laboratory strains revealed the outbreak clone shared a common ancestor with isolates from the United States and United Kingdom (UK). A strain obtained in the UK was phylogenetically related to our outbreak. Follow-up of that case revealed the patient had recently returned from Australia. Our data demonstrate new C. difficile strains are an on-going threat, with potential for rapid spread. Active surveillance is needed to identify and control emerging lineages.

Effects of proton pump inhibitors and histamine-2 receptor antagonists on response to fidaxomicin or vancomycin in patients with Clostridium difficile-associated diarrhoea

OBJECTIVE

It has been established that use of proton pump inhibitors (PPIs) is associated with an increased risk of acquiring Clostridium difficile-associated diarrhoea (CDAD). However, it is not known whether the use of PPIs or histamine-2 receptor antagonists (H2RAs) concurrently with CDAD-targeted antibiotic treatment affects clinical response or recurrence rates.

DESIGN

In two phase 3 trials, patients with toxin-positive CDAD were randomised to receive fidaxomicin 200 mg twice daily or vancomycin 125 mg four times daily for 10 days. Only inpatients with CDAD (due to complete medication record availability) were included in this post hoc analysis: 701 patients, of whom 446 (64%) used PPIs or H2RAs during study drug treatment or follow-up. Baseline factors that were statistically significant in univariate analyses were analysed in multivariate analyses of effects on clinical response and recurrence.

RESULTS

Multivariate analysis showed that leukocytosis, elevated creatinine and hypoalbuminemia, but not PPI or H2RA use, were significant factors associated with poor clinical responses. Treatment group was the single significant predictor of recurrence; the probability of recurrence after fidaxomicin therapy was half that following vancomycin therapy.

CONCLUSIONS

Acid-suppressing drugs, used by nearly two-thirds of inpatients with CDAD, did not worsen clinical response or recurrence when used concurrently with fidaxomicin or vancomycin. Therefore, development of CDAD does not require discontinuation of anti-acid treatment in patients who have an indication for continuing PPI or H2RA therapy, such as gastro-oesophageal reflux disease and risk of gastrointestinal bleed.

Colonization with Enterobacteriaceae producing ESBLs in children attending pre-school childcare facilities in the Lao People's Democratic Republic

Objectives

Intestinal carriage constitutes an important reservoir of antimicrobial-resistant bacteria, with some of the highest rates reported from Asia. Antibiotic resistance has been little studied in Laos, where some antibiotics are available without restriction, but others such as carbapenems are not available.

Patients and methods

We collected stools from 397 healthy children in 12 randomly selected pre-school childcare facilities in and around Vientiane. Colonization with ESBL-producing Enterobacteriaceae (ESBLE) and carbapenemase-producing Enterobacteriaceae (CPE) was detected using a disc diffusion screening test and ESBLE were characterized using WGS. Risk factor data were collected by questionnaire.

Results

Ninety-two children (23%) were colonized with ESBLE, mainly Escherichia coli carrying bla_CTX-M and Klebsiella pneumoniae carrying bla_SHV or bla_CTX-M, which were frequently resistant to multiple antibiotic classes. Although residence in Vientiane Capital, foreign travel, higher maternal level of education, antibiotic use in the preceding 3 months and attending a childcare facility with a ‘good’ level of hygiene were all associated with ESBLE colonization on univariable analysis, a significant association remained only for antibiotic use when a stepwise approach was used with a multivariate random-effects model. WGS analysis suggested transmission in both childcare facilities and community settings.

Conclusions

The high prevalence of paediatric colonization with ESBLE in Laos, one of the highest reported in Asia, is probably the result of inappropriate antibiotic use. Paediatric colonization with CPE was not identified in this study, but it is important to continue to monitor the spread of antibiotic-resistant Enterobacteriaceae in Laos.

The hospital microbiome project: meeting report for the UK science and innovation network UK-USA workshop ‘beating the superbugs: hospital microbiome studies for tackling antimicrobial resistance’, October 14th 2013

The UK Science and Innovation Network UK-USA workshop ‘Beating the Superbugs: Hospital Microbiome Studies for tackling Antimicrobial Resistance’ was held on October 14th 2013 at the UK Department of Health, London. The workshop was designed to promote US-UK collaboration on hospital microbiome studies to add a new facet to our collective understanding of antimicrobial resistance. The assembled researchers debated the importance of the hospital microbial community in transmission of disease and as a reservoir for antimicrobial resistance genes, and discussed methodologies, hypotheses, and priorities. A number of complementary approaches were explored, although the importance of the built environment microbiome in disease transmission was not universally accepted. Current whole genome epidemiological methods are being pioneered in the UK and the benefits of moving to community analysis are not necessarily obvious to the pioneers; however, rapid progress in other areas of microbiology suggest to some researchers that hospital microbiome studies will be exceptionally fruitful even in the short term. Collaborative studies will recombine different strengths to tackle the international problems of antimicrobial resistance and hospital and healthcare associated infections.

Genome sequencing of an extended series of NDM-producing Klebsiella pneumoniae isolates from neonatal infections in a Nepali hospital characterizes the extent of community- versus hospital-associated transmission in an endemic setting

NDM-producing Klebsiella pneumoniae strains represent major clinical and infection control challenges, particularly in resource-limited settings with high rates of antimicrobial resistance. Determining whether transmission occurs at a gene, plasmid, or bacterial strain level and within hospital and/or the community has implications for monitoring and controlling spread. Whole-genome sequencing (WGS) is the highest-resolution typing method available for transmission epidemiology. We sequenced carbapenem-resistant K. pneumoniae isolates from 26 individuals involved in several infection case clusters in a Nepali neonatal unit and 68 other clinical Gram-negative isolates from a similar time frame, using Illumina and PacBio technologies. Within-outbreak chromosomal and closed-plasmid structures were generated and used as data set-specific references. Three temporally separated case clusters were caused by a single NDM K. pneumoniae strain with a conserved set of four plasmids, one being a 304,526-bp plasmid carrying bla(NDM-1). The plasmids contained a large number of antimicrobial/heavy metal resistance and plasmid maintenance genes, which may have explained their persistence. No obvious environmental/human reservoir was found. There was no evidence of transmission of outbreak plasmids to other Gram-negative clinical isolates, although bla(NDM) variants were present in other isolates in different genetic contexts. WGS can effectively define complex antimicrobial resistance epidemiology. Wider sampling frames are required to contextualize outbreaks. Infection control may be effective in terminating outbreaks caused by particular strains, even in areas with widespread resistance, although this study could not demonstrate evidence supporting specific interventions. Larger, detailed studies are needed to characterize resistance genes, vectors, and host strains involved in disease, to enable effective intervention.

The complexity and diversity of the Pathogenicity Locus in Clostridium difficile clade 5

The symptoms of Clostridium difficile infection are caused by two closely related toxins, TcdA and TcdB, which are encoded by the 19.6 kb Pathogenicity Locus (PaLoc). The PaLoc is variably present among strains, and in this respect it resembles a mobile genetic element. The C. difficile population structure consists mainly of five phylogenetic clades designated 1–5. Certain genotypes of clade 5 are associated with recently emergent highly pathogenic strains causing human disease and animal infections. The aim of this study was to explore the evolutionary history of the PaLoc in C. difficile clade 5. Phylogenetic analyses and annotation of clade 5 PaLoc variants and adjoining genomic regions were undertaken using a representative collection of toxigenic and nontoxigenic strains. Comparison of the core genome and PaLoc phylogenies obtained for clade 5 and representatives of the other clades identified two distinct PaLoc acquisition events, one involving a toxin A⁺B⁺ PaLoc variant and the other an A⁻B⁺ variant. Although the exact mechanism of each PaLoc acquisition is unclear, evidence of possible homologous recombination with other clades and between clade 5 lineages was found within the PaLoc and adjacent regions. The generation of nontoxigenic variants by PaLoc loss via homologous recombination with PaLoc-negative members of other clades was suggested by analysis of cdu2, although none is likely to have occurred recently. A variant of the putative holin gene present in the clade 5 A⁻B⁺ PaLoc was likely acquired via allelic exchange with an unknown element. Fine-scale phylogenetic analysis of C. difficile clade 5 revealed the extent of its genetic diversity, consistent with ancient evolutionary origins and a complex evolutionary history for the PaLoc.

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.

Mobile elements drive recombination hotspots in the core genome of Staphylococcus aureus

Horizontal gene transfer is an important driver of bacterial evolution, but genetic exchange in the core genome of clonal species, including the major pathogen Staphylococcus aureus, is incompletely understood. Here we reveal widespread homologous recombination in S. aureus at the species level, in contrast to its near-complete absence between closely related strains. We discover a patchwork of hotspots and coldspots at fine scales falling against a backdrop of broad-scale trends in rate variation. Over megabases, homoplasy rates fluctuate 1.9-fold, peaking towards the origin-of-replication. Over kilobases, we find core recombination hotspots of up to 2.5-fold enrichment situated near fault lines in the genome associated with mobile elements. The strongest hotspots include regions flanking conjugative transposon ICE6013, the staphylococcal cassette chromosome (SCC) and genomic island νSaα. Mobile element-driven core genome transfer represents an opportunity for adaptation and challenges our understanding of the recombination landscape in predominantly clonal pathogens, with important implications for genotype–phenotype mapping.

Horizontal gene transfer occurs in most bacteria, yet it is unclear whether it happens in clonal species. Here, Everitt et al. show widespread within-species recombination, driven by mobile elements, in the genome of the pathogen Staphylococcus aureus, but no recombination between closely related strains.

Reduction of invasive pneumococcal disease 3 years after the introduction of the 13-valent conjugate vaccine in the Oxfordshire region of England

BACKGROUND

The 7-valent pneumococcal conjugate (PCV7) vaccine's impact on invasive pneumococcal disease (IPD) is well described, but few reports exist on the additional impact of the 13-valent vaccine (PCV13).

METHODS

We calculated the IPD incidence across all ages in a surveillance project following implementation of PCV7 (in September 2006) and PCV13 (in April 2010) in children aged <2 years (11 hospitals; 4935 cases).

RESULTS

The overall incidence decreased from 10 cases/100 000 persons per year in 1996-1997 to 8 cases/100 000 persons per year in 2007-2008 and 7 cases/100 000 in 2012-2013. Declines were greater in children aged <2 years (from 37 cases/100 000 in 1996-1997 to 29 and 14 cases/100 000 in 2007-2008 and 2012-2013, respectively). The incidence of IPD due to PCV7 serotypes decreased in all ages after PCV7 introduction (P < .001), whereas the incidence of IPD due to the additional 6 serotypes in PCV13 and to nonvaccine types (NVTs) increased in children aged ≥2 years (P < .001 for both comparisons). The incidence of IPD due to the 6 additional serotypes in PCV13 declined significantly after PCV13 introduction in all ages (P ≤ .01), and the incidence of IPD due to NVTs declined significantly in children aged ≥2 years (P = .003). In 2011-2013, the overall incidences of IPD due to PCV7 serotypes, the 6 additional serotypes in PCV13, and NVTs were 0.3, 2.8, and 4.4 cases/100 000; the incidences among children aged <2 years were 0.9, 2.4, and 10.8 cases/100 000, respectively.

CONCLUSIONS

The annual incidence of IPD due to vaccine serotypes (1-3 cases/100 000) among children aged <2 years and nontarget groups demonstrates the success of PCV7 and PCV13. A substantially higher incidence of IPD due to NVTs indicates the importance of ongoing surveillance and extension of vaccine polyvalency.

Prevalence of Staphylococcus aureus protein A (spa) mutants in the community and hospitals in Oxfordshire

Background

Staphylococcal protein A (spa) is an important virulence factor which enables Staphylococcus aureus to evade host immune responses. Genotypes known as “spa-types”, based on highly variable Xr region sequences of the spa-gene, are frequently used to classify strains. A weakness of current spa-typing primers is that rearrangements in the IgG-binding region of the gene cause 1-2% of strains to be designated as “non-typeable”.

Results

We developed an improved primer which enabled sequencing of all strains, containing any type of genetic rearrangement, in a large study among community carriers and hospital inpatients in Oxfordshire, UK (6110 isolates). We identified eight novel spa-gene variants, plus one previously described. Three of these rearrangements would be designated “non-typeable” using current spa-typing methods; they occurred in 1.8% (72/3905) asymptomatically carried and 0.6% (14/2205) inpatient S. aureus strains. Some individuals were simultaneously colonized by both formerly non-typeable and typeable strains; previously such patients would have been identified as carrying only currently typeable strains, underestimating mixed carriage prevalence and diversity. Formerly non-typeable strains were found in more spa-types associated with multilocus sequence type ST398 (35%), common among livestock, compared to other groups with any non-typeable strains (1-4%), suggesting particular spa-types may have been under-represented in previous human studies.

Conclusions

This improved method allows us to spa-type previously non-typeable strains with rearrangements in the spa-gene and to resolve cases of mixed colonization with deletions in one or more strains, thus accounting for hidden diversity of S. aureus in both community and hospital environments.