Evolutionary Processes Driving the Rise and Fall of Staphylococcus aureus ST239, a Dominant Hybrid Pathogen.

Selection plays a key role in the spread of antibiotic resistance, but the evolutionary drivers of clinically important resistant strains remain poorly understood. Here, we use genomic analyses and competition experiments to study Staphylococcus aureus ST239, a prominent MRSA strain that is thought to have been formed by large-scale recombination between ST8 and ST30. Genomic analyses allowed us to refine the hybrid model for the origin of ST239 and to date the origin of ST239 to 1920 to 1945, which predates the clinical introduction of methicillin in 1959. Although purifying selection has dominated the evolution of ST239, parallel evolution has occurred in genes involved in antibiotic resistance and virulence, suggesting that ST239 has evolved toward an increasingly pathogenic lifestyle. Crucially, ST239 isolates have low competitive fitness relative to both ST8 and ST30 isolates, supporting the idea that fitness costs have driven the demise of this once-dominant pathogen strain. IMPORTANCE The rise of antibiotic resistance in most pathogenic bacteria has been driven by the spread of a small number of epidemically successful resistant strains. However, the processes that drive the rise and fall of these "superbugs" remain poorly understood. In our study, we investigated Staphylococcus aureus ST239, an important MRSA strain that has been a leading cause of serious hospital-acquired infections. We show here that ST239 was formed by the exchange of a very large fragment of DNA carrying resistance genes between strains of S. aureus sometime before 1945. The introduction of the antibiotic methicillin in 1959 provided a key advantage for Methicillin-Resistant Staphylococcus aureus (MRSA) strains. Interestingly, we found that ST239 has low competitive ability compared to other strains of S. aureus, and this evolutionary hindrance may explain why the prevalence of ST239 has declined precipitously at a global scale.

Rapid evolution and host immunity drive the rise and fall of carbapenem resistance during an acute Pseudomonas aeruginosa infection.

It is well established that antibiotic treatment selects for resistance, but the dynamics of this process during infections are poorly understood. Here we map the responses of Pseudomonas aeruginosa to treatment in high definition during a lung infection of a single ICU patient. Host immunity and antibiotic therapy with meropenem suppressed P. aeruginosa, but a second wave of infection emerged due to the growth of oprD and wbpM meropenem resistant mutants that evolved in situ. Selection then led to a loss of resistance by decreasing the prevalence of low fitness oprD mutants, increasing the frequency of high fitness mutants lacking the MexAB-OprM efflux pump, and decreasing the copy number of a multidrug resistance plasmid. Ultimately, host immunity suppressed wbpM mutants with high meropenem resistance and fitness. Our study highlights how natural selection and host immunity interact to drive both the rapid rise, and fall, of resistance during infection.

Efflux pump activity potentiates the evolution of antibiotic resistance across S. aureus isolates.

The rise of antibiotic resistance in many bacterial pathogens has been driven by the spread of a few successful strains, suggesting that some bacteria are genetically pre-disposed to evolving resistance. Here, we test this hypothesis by challenging a diverse set of 222 isolates of Staphylococcus aureus with the antibiotic ciprofloxacin in a large-scale evolution experiment. We find that a single efflux pump, norA, causes widespread variation in evolvability across isolates. Elevated norA expression potentiates evolution by increasing the fitness benefit provided by DNA topoisomerase mutations under ciprofloxacin treatment. Amplification of norA provides a further mechanism of rapid evolution in isolates from the CC398 lineage. Crucially, chemical inhibition of NorA effectively prevents the evolution of resistance in all isolates. Our study shows that pre-existing genetic diversity plays a key role in shaping resistance evolution, and it may be possible to predict which strains are likely to evolve resistance and to optimize inhibitor use to prevent this outcome.

Insidious Risk of Severe Mycobacterium chimaera Infection in Cardiac Surgery Patients.

BACKGROUND: An urgent UK investigation was launched to assess risk of invasive Mycobacterium chimaera infection in cardiothoracic surgery and a possible association with cardiopulmonary bypass heater-cooler units following alerts in Switzerland and The Netherlands. METHODS: Parallel investigations were pursued: (1) identification of cardiopulmonary bypass-associated M. chimaera infection through national laboratory and hospital admissions data linkage; (2) cohort study to assess patient risk; (3) microbiological and aerobiological investigations of heater-coolers in situ and under controlled laboratory conditions; and (4) whole-genome sequencing of clinical and environmental isolates. RESULTS: Eighteen probable cases of cardiopulmonary bypass-associated M. chimaera infection were identified; all except one occurred in adults. Patients had undergone valve replacement in 11 hospitals between 2007 and 2015, a median of 19 months prior to onset (range, 3 months to 5 years). Risk to patients increased after 2010 from <0.2 to 1.65 per 10000 person-years in 2013, a 9-fold rise for infections within 2 years of surgery (rate ratio, 9.08 [95% CI, 1.81-87.76]). Endocarditis was the most common presentation (n = 11). To date, 9 patients have died. Investigations identified aerosol release through breaches in heater-cooler tanks. Mycobacterium chimaera and other pathogens were recovered from water and air samples. Phylogenetic analysis found close clustering of strains from probable cases. CONCLUSIONS: We identified low but escalating risk of severe M. chimaera infection associated with heater-coolers with cases in a quarter of cardiothoracic centers. Our investigations strengthen etiological evidence for the role of heater-coolers in transmission and raise the possibility of an ongoing, international point-source outbreak. Active management of heater-coolers and heightened clinical awareness are imperative given the consequences of infection.

Practical Approaches for Detecting Selection in Microbial Genomes.

Microbial genome evolution is shaped by a variety of selective pressures. Understanding how these processes occur can help to address important problems in microbiology by explaining observed differences in phenotypes, including virulence and resistance to antibiotics. Greater access to whole-genome sequencing provides microbiologists with the opportunity to perform large-scale analyses of selection in novel settings, such as within individual hosts. This tutorial aims to guide researchers through the fundamentals underpinning popular methods for measuring selection in pathogens. These methods are transferable to a wide variety of organisms, and the exercises provided are designed for researchers with any level of programming experience.

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.

Bacterial phylogenetic reconstruction from whole genomes is robust to recombination but demographic inference is not.

UNLABELLED: Phylogenetic inference in bacterial genomics is fundamental to understanding problems such as population history, antimicrobial resistance, and transmission dynamics. The field has been plagued by an apparent state of contradiction since the distorting effects of recombination on phylogeny were discovered more than a decade ago. Researchers persist with detailed phylogenetic analyses while simultaneously acknowledging that recombination seriously misleads inference of population dynamics and selection. Here we resolve this paradox by showing that phylogenetic tree topologies based on whole genomes robustly reconstruct the clonal frame topology but that branch lengths are badly skewed. Surprisingly, removing recombining sites can exacerbate branch length distortion caused by recombination. IMPORTANCE: Phylogenetic tree reconstruction is a popular approach for understanding the relatedness of bacteria in a population from differences in their genome sequences. However, bacteria frequently exchange regions of their genomes by a process called homologous recombination, which violates a fundamental assumption of phylogenetic methods. Since many researchers continue to use phylogenetics for recombining bacteria, it is important to understand how recombination affects the conclusions drawn from these analyses. We find that whole-genome sequences afford great accuracy in reconstructing evolutionary relationships despite concerns surrounding the presence of recombination, but the branch lengths of the phylogenetic tree are indeed badly distorted. Surprisingly, methods to reduce the impact of recombination on branch lengths can exacerbate the problem.

Whole genome sequencing and de novo assembly identifies Sydney-like variant noroviruses and recombinants during the winter 2012/2013 outbreak in England.

BACKGROUND: Norovirus is the commonest cause of epidemic gastroenteritis among people of all ages. Outbreaks frequently occur in hospitals and the community, costing the UK an estimated £110 m per annum. An evolutionary explanation for periodic increases in norovirus cases, despite some host-specific post immunity is currently limited to the identification of obvious recombinants. Our understanding could be significantly enhanced by full length genome sequences for large numbers of intensively sampled viruses, which would also assist control and vaccine design. Our objective is to develop rapid, high-throughput, end-to-end methods yielding complete norovirus genome sequences. We apply these methods to recent English outbreaks, placing them in the wider context of the international norovirus epidemic of winter 2012. METHOD: Norovirus sequences were generated from 28 unique clinical samples by Illumina RNA sequencing (RNA-Seq) of total faecal RNA. A range of de novo sequence assemblers were attempted. The best assembler was identified by validation against three replicate samples and two norovirus qPCR negative samples, together with an additional 20 sequences determined by PCR and fractional capillary sequencing. Phylogenetic methods were used to reconstruct evolutionary relationships from the whole genome sequences. RESULTS: Full length norovirus genomes were generated from 23/28 samples. 5/28 partial norovirus genomes were associated with low viral copy numbers. The de novo assembled sequences differed from sequences determined by capillary sequencing by <0.003%. Intra-host nucleotide sequence diversity was rare, but detectable by mapping short sequence reads onto its de novo assembled consensus. Genomes similar to the Sydney 2012 strain caused 78% (18/23) of cases, consistent with its previously documented association with the winter 2012 global outbreak. Interestingly, phylogenetic analysis and recombination detection analysis of the consensus sequences identified two related viruses as recombinants, containing sequences in prior circulation to Sydney 2012 in open reading frame (ORF) 2. CONCLUSION: Our approach facilitates the rapid determination of complete norovirus genomes. This method provides high resolution of full norovirus genomes which, when coupled with detailed epidemiology, may improve the understanding of evolution and control of this important healthcare-associated pathogen.

Evolutionary dynamics of local pandemic H1N1/2009 influenza virus lineages revealed by whole-genome analysis.

Virus gene sequencing and phylogenetics can be used to study the epidemiological dynamics of rapidly evolving viruses. With complete genome data, it becomes possible to identify and trace individual transmission chains of viruses such as influenza virus during the course of an epidemic. Here we sequenced 153 pandemic influenza H1N1/09 virus genomes from United Kingdom isolates from the first (127 isolates) and second (26 isolates) waves of the 2009 pandemic and used their sequences, dates of isolation, and geographical locations to infer the genetic epidemiology of the epidemic in the United Kingdom. We demonstrate that the epidemic in the United Kingdom was composed of many cocirculating lineages, among which at least 13 were exclusively or predominantly United Kingdom clusters. The estimated divergence times of two of the clusters predate the detection of pandemic H1N1/09 virus in the United Kingdom, suggesting that the pandemic H1N1/09 virus was already circulating in the United Kingdom before the first clinical case. Crucially, three clusters contain isolates from the second wave of infections in the United Kingdom, two of which represent chains of transmission that appear to have persisted within the United Kingdom between the first and second waves. This demonstrates that whole-genome analysis can track in fine detail the behavior of individual influenza virus lineages during the course of a single epidemic or pandemic.

Evolutionary history of rabies in Ghana.

Rabies virus (RABV) is enzootic throughout Africa, with the domestic dog (Canis familiaris) being the principal vector. Dog rabies is estimated to cause 24,000 human deaths per year in Africa, however, this estimate is still considered to be conservative. Two sub-Saharan African RABV lineages have been detected in West Africa. Lineage 2 is present throughout West Africa, whereas Africa 1a dominates in northern and eastern Africa, but has been detected in Nigeria and Gabon, and Africa 1b was previously absent from West Africa. We confirmed the presence of RABV in a cohort of 76 brain samples obtained from rabid animals in Ghana collected over an eighteen-month period (2007-2009). Phylogenetic analysis of the sequences obtained confirmed all viruses to be RABV, belonging to lineages previously detected in sub-Saharan Africa. However, unlike earlier reported studies that suggested a single lineage (Africa 2) circulates in West Africa, we identified viruses belonging to the Africa 2 lineage and both Africa 1 (a and b) sub-lineages. Phylogeographic Bayesian Markov chain Monte Carlo analysis of a 405 bp fragment of the RABV nucleoprotein gene from the 76 new sequences derived from Ghanaian animals suggest that within the Africa 2 lineage three clades co-circulate with their origins in other West African countries. Africa 1a is probably a western extension of a clade circulating in central Africa and the Africa 1b virus a probable recent introduction from eastern Africa. We also developed and tested a novel reverse-transcription loop-mediated isothermal amplification (RT-LAMP) assay for the detection of RABV in African laboratories. This RT-LAMP was shown to detect both Africa 1 and 2 viruses, including its adaptation to a lateral flow device format for product visualization. These data suggest that RABV epidemiology is more complex than previously thought in West Africa and that there have been repeated introductions of RABV into Ghana. This analysis highlights the potential problems of individual developing nations implementing rabies control programmes in the absence of a regional programme.

Emergence, dominance, and possible decline of CXCR4 chemokine receptor usage during the course of HIV infection.

Binding to a chemokine receptor, either CCR5 or CXCR4, by the gp120 glycoprotein is an essential step in the pathway by which HIV enters host cells. Recently, CCR5 antagonists have been developed that obstruct binding of CCR5 by gp120, thus inhibiting host cell entry. Resistance to such CCR5 antagonists may emerge, however, through the selection of viral strains capable of utilizing CXCR4 receptors. This study explores the evolutionary context of emergence, and in many cases decline, of dominant CXCR4-usage (X4) during disease progression within a number of individuals. Of seven individuals exhibiting a switch to dominant CXCR4 usage, such dominance is transient in five of them with CCR5-usage (R5) re-emerging to dominate the viral population later in disease progression. Three individuals conform to documented X4 transience in that the re-emergence of R5 dominance is an outgrowth from the predominant R5 strain. However, in two individuals we observe a novel pathway for R5 re-emergence in that R5 strains emerge to dominate late in disease progression through continued evolution of the X4 population. This suggests that the molecular mechanism of such switches between R5 and X4-usage is strain specific and that no single mechanism is shared between individuals. These findings have implications for the understanding of the mechanisms of potential emergence of resistance to CCR5 antagonists through use of the CXCR4 receptor and support the importance to have an appropriately optimized background therapy for use with entry inhibitors and, as for all HAART, to monitor drug resistance in a comprehensive manner.