A genomic portrait of the emergence, evolution, and global spread of a methicillin-resistant Staphylococcus aureus pandemic.

The widespread use of antibiotics in association with high-density clinical care has driven the emergence of drug-resistant bacteria that are adapted to thrive in hospitalized patients. Of particular concern are globally disseminated methicillin-resistant Staphylococcus aureus (MRSA) clones that cause outbreaks and epidemics associated with health care. The most rapidly spreading and tenacious health-care-associated clone in Europe currently is EMRSA-15, which was first detected in the UK in the early 1990s and subsequently spread throughout Europe and beyond. Using phylogenomic methods to analyze the genome sequences for 193 S. aureus isolates, we were able to show that the current pandemic population of EMRSA-15 descends from a health-care-associated MRSA epidemic that spread throughout England in the 1980s, which had itself previously emerged from a primarily community-associated methicillin-sensitive population. The emergence of fluoroquinolone resistance in this EMRSA-15 subclone in the English Midlands during the mid-1980s appears to have played a key role in triggering pandemic spread, and occurred shortly after the first clinical trials of this drug. Genome-based coalescence analysis estimated that the population of this subclone over the last 20 yr has grown four times faster than its progenitor. Using comparative genomic analysis we identified the molecular genetic basis of 99.8% of the antimicrobial resistance phenotypes of the isolates, highlighting the potential of pathogen genome sequencing as a diagnostic tool. We document the genetic changes associated with adaptation to the hospital environment and with increasing drug resistance over time, and how MRSA evolution likely has been influenced by country-specific drug use regimens.

Molecular tracing of the emergence, adaptation, and transmission of hospital-associated methicillin-resistant Staphylococcus aureus.

Hospital-associated infections caused by methicillin-resistant Staphylococcus aureus (MRSA) are a global health burden dominated by a small number of bacterial clones. The pandemic EMRSA-16 clone (ST36-II) has been widespread in UK hospitals for 20 y, but its evolutionary origin and the molecular basis for its hospital association are unclear. We carried out a Bayesian phylogenetic reconstruction on the basis of the genome sequences of 87 S. aureus isolates including 60 EMRSA-16 and 27 additional clonal complex 30 (CC30) isolates, collected from patients in three continents over a 53-y period. The three major pandemic clones to originate from the CC30 lineage, including phage type 80/81, Southwest Pacific, and EMRSA-16, shared a most recent common ancestor that existed over 100 y ago, whereas the hospital-associated EMRSA-16 clone is estimated to have emerged about 35 y ago. Our CC30 genome-wide analysis revealed striking molecular correlates of hospital- or community-associated pandemics represented by mobile genetic elements and nonsynonymous mutations affecting antibiotic resistance and virulence. Importantly, phylogeographic analysis indicates that EMRSA-16 spread within the United Kingdom by transmission from hospitals in large population centers in London and Glasgow to regional health-care settings, implicating patient referrals as an important cause of nationwide transmission. Taken together, the high-resolution phylogenomic approach used resulted in a unique understanding of the emergence and transmission of a major MRSA clone and provided molecular correlates of its hospital adaptation. Similar approaches for hospital-associated clones of other bacterial pathogens may inform appropriate measures for controlling their intra- and interhospital spread.

Comprehensive Profiling of Poor-Risk Paired Primary and Recurrent Triple-Negative Breast Cancers Reveals Immune Phenotype Shifts.

PURPOSE: Emerging data suggest immune checkpoint inhibitors have reduced efficacy in heavily pretreated triple-negative breast cancers (TNBC), but underlying mechanisms are poorly understood. To better understand the phenotypic evolution of TNBCs, we studied the genomic and transcriptomic profiles of paired tumors from patients with TNBC. EXPERIMENTAL DESIGN: We collected paired primary and metastatic TNBC specimens from 43 patients and performed targeted exome sequencing and whole-transcriptome sequencing. From these efforts, we ascertained somatic mutation profiles, tumor mutational burden (TMB), TNBC molecular subtypes, and immune-related gene expression patterns. Stromal tumor-infiltrating lymphocytes (stromal TIL), recurrence-free survival, and overall survival were also analyzed. RESULTS: We observed a typical TNBC mutational landscape with minimal shifts in copy number or TMB over time. However, there were notable TNBC molecular subtype shifts, including increases in the Lehmann/Pietenpol-defined basal-like 1 (BL1, 11.4%-22.6%) and mesenchymal (M, 11.4%-22.6%) phenotypes, and a decrease in the immunomodulatory phenotype (IM, 31.4%-3.2%). The Burstein-defined basal-like immune-activated phenotype was also decreased (BLIA, 42.2%-17.2%). Among downregulated genes from metastases, we saw enrichment of immune-related Kyoto Encyclopedia of Genes and Genomes pathways and gene ontology (GO) terms, and decreased expression of immunomodulatory gene signatures (P < 0.03) and percent stromal TILs (P = 0.03). There was no clear association between stromal TILs and survival. CONCLUSIONS: We observed few mutational shifts, but largely consistent transcriptomic shifts in longitudinally paired TNBCs. Transcriptomic and IHC analyses revealed significantly reduced immune-activating gene expression signatures and TILs in recurrent TNBCs. These data may explain the observed lack of efficacy of immunotherapeutic agents in heavily pretreated TNBCs. Further studies are ongoing to better understand these initial observations.See related commentary by Savas and Loi, p. 526.

Tumor Fusion Burden as a Hallmark of Immune Infiltration in Prostate Cancer.

Prostate cancer is the second leading cause of cancer-related death in men. Despite having a relatively lower tumor mutational burden than most tumor types, multiple gene fusions such as TMPRSS2:ERG have been characterized and linked to more aggressive disease. Individual tumor samples have been found to contain multiple fusions, and it remains unknown whether these fusions increase tumor immunogenicity. Here, we investigated the role of fusion burden on the prevalence and expression of key molecular and immune effectors in prostate cancer tissue specimens that represented the different stages of disease progression and androgen sensitivity, including hormone-sensitive and castration-resistant prostate cancer. We found that tumor fusion burden was inversely correlated with tumor mutational burden and not associated with disease stage. High fusion burden correlated with high immune infiltration, PD-L1 expression on immune cells, and immune signatures, representing activation of T cells and M1 macrophages. High fusion burden inversely correlated with immune-suppressive signatures. Our findings suggest that high tumor fusion burden may be a more appropriate biomarker than tumor mutational burden in prostate cancer, as it more closely associates with immunogenicity, and suggests that tumors with high fusion burden could be potential candidates for immunotherapeutic agents.

Evolutionary Dynamics of Pandemic Methicillin-Sensitive Staphylococcus aureus ST398 and Its International Spread via Routes of Human Migration.

Methicillin-susceptible Staphylococcus aureus (MSSA) accounts for the majority of S. aureus infections globally, and yet surprisingly little is known about its clonal evolution. We applied comparative whole-genome sequencing (WGS) analyses to epidemiologically and geographically diverse ST398-MSSA, a pandemic lineage affecting both humans and livestock. Bayesian phylogenetic analysis predicted divergence of human-associated ST398-MSSA ~40 years ago. Isolates from Midwestern pigs and veterinarians differed substantially from those in New York City (NYC). Pig ST398 strains contained a large region of recombination representing imports from multiple sequence types (STs). Phylogeographic analyses supported the spread of ST398-MSSA along local cultural and migratory links between parts of the Caribbean, North America, and France, respectively. Applying pairwise single-nucleotide polymorphism (SNP) distances as a measure of genetic relatedness between isolates, we observed that ST398 not only clustered in households but also frequently extended across local social networks. Isolates collected from environmental surfaces reflected the full diversity of colonizing individuals, highlighting their potentially critical role as reservoirs for transmission and diversification. Strikingly, we observed high within-host SNP variability compared to our previous studies on the dominant methicillin-resistant Staphylococcus aureus (MRSA) clone USA300. Our data indicate that the dynamics of colonization, persistence, and transmission differ substantially between USA300-MRSA and ST398-MSSA. Taken together, our study reveals local and international routes of transmission for a major MSSA clone, indicating key impacts of recombination and mutation on genetic diversification and highlighting important ecological differences from epidemic USA300. Our study demonstrates extensive local and international routes of transmission for a major MSSA clone despite the lack of substantial antibiotic resistance. IMPORTANCE: Unlike methicillin-resistant Staphylococcus aureus (MRSA), surprisingly little is known about the clonal evolution of methicillin-susceptible S. aureus (MSSA), although these strains account for the majority of S. aureus infections. To better understand how MSSA spreads and becomes established in communities, we applied comparative bacterial whole-genome sequencing to pandemic ST398-MSSA, a clone of clinical importance affecting humans and livestock in different geographic regions. Phylogeographic analyses identified that ST398-MSSA spread along local cultural and migratory links between parts of the Caribbean, North America, and France, respectively. We observed high within-host SNP variability compared to our previous studies on the dominant MRSA clone USA300. Our data indicate that the dynamics of colonization, persistence, and transmission differ substantially between USA300 MRSA and ST398 MSSA.

Identification of source and sink populations for the emergence and global spread of the East-Asia clone of community-associated MRSA.

BACKGROUND: Our understanding of the factors influencing the emergence, dissemination and global distribution of epidemic clones of bacteria is limited. ST59 is a major epidemic clone of community-associated MRSA in East Asia, responsible for extensive morbidity and mortality, but has a much lower prevalence in other parts of the world. The geographic origin of ST59 and its international routes of dissemination are unclear and disputed in the literature. RESULTS: To investigate the origin and spread of the ST59 clone, we obtained whole genome sequences of isolates from four continents, sampled over more than a decade, and carried out a time-scaled phylogeographic analysis. We discover that two distinct ST59 clades emerged concurrently, in East Asia and the USA, but underwent clonal expansion at different times. The East Asia clade was strongly enriched for gene determinants associated with antibiotic resistance, consistent with regional differences in antibiotic usage. Both clones spread independently to Australia and Europe, and we found evidence of the persistence of multi-drug resistance following export from East Asia. Direct transfer of strains between Taiwan and the USA was not observed in either direction, consistent with geographic niche exclusion. CONCLUSIONS: Our results resolve a longstanding controversy regarding the origin of the ST59 clone, revealing the major global source and sink populations and routes for the spread of multi-drug resistant clones. Additionally, our findings indicate that diversification of the accessory genome of epidemic clones partly reflects region-specific patterns of antibiotic usage, which may influence bacterial fitness after transmission to different geographic locations.

A single natural nucleotide mutation alters bacterial pathogen host tropism.

The capacity of microbial pathogens to alter their host tropism leading to epidemics in distinct host species populations is a global public and veterinary health concern. To investigate the molecular basis of a bacterial host-switching event in a tractable host species, we traced the evolutionary trajectory of the common rabbit clone of Staphylococcus aureus. We report that it evolved through a likely human-to-rabbit host jump over 40 years ago and that only a single naturally occurring nucleotide mutation was required and sufficient to convert a human-specific S. aureus strain into one that could infect rabbits. Related mutations were identified at the same locus in other rabbit strains of distinct clonal origin, consistent with convergent evolution. This first report of a single mutation that was sufficient to alter the host tropism of a microorganism during its evolution highlights the capacity of some pathogens to readily expand into new host species populations.

Recombination-mediated remodelling of host-pathogen interactions during Staphylococcus aureus niche adaptation.

Large-scale recombination events have led to the emergence of epidemic clones of several major bacterial pathogens. However, the functional impact of the recombination on clonal success is not understood. Here, we identified a novel widespread hybrid clone (ST71) of livestock-associated Staphylococcus aureus that evolved from an ancestor belonging to the major bovine lineage CC97, through multiple large-scale recombination events with other S. aureus lineages occupying the same ruminant niche. The recombination events, affecting a 329 kb region of the chromosome spanning the origin of replication, resulted in allele replacement and loss or gain of an array of genes influencing host-pathogen interactions. Of note, molecular functional analyses revealed that the ST71 hybrid clone has acquired multiple novel pathogenic traits associated with acquired and innate immune evasion and bovine extracellular matrix adherence. These findings provide a paradigm for the impact of large-scale recombination events on the rapid evolution of bacterial pathogens within defined ecological niches.

High-throughput sequencing for the study of bacterial pathogen biology.

A revolution in sequencing technologies in recent years has led to dramatically increased throughput and reduced cost of bacterial genome sequencing. An increasing number of applications of the new technologies are providing broad insights into bacterial evolution, epidemiology, and pathogenesis. For example, the capacity to sequence large numbers of bacterial isolates is enabling high resolution phylogenetic analyses of bacterial populations leading to greatly enhanced understanding of the emergence, adaptation, and transmission of pathogenic clones. In addition, RNA-seq offers improved quantification and resolution for transcriptomic analysis, and the combination of high-throughput sequencing with transposon mutagenesis is a powerful approach for the identification of bacterial determinants required for survival in vivo. In this concise review we provide selected examples of how high throughput sequencing is being applied to understand the biology of bacterial pathogens, and discuss future technological advances likely to have a profound impact on the field.

Gene flow in environmental Legionella pneumophila leads to genetic and pathogenic heterogeneity within a Legionnaires' disease outbreak.

BACKGROUND: Legionnaires' disease is a severe form of pneumonia caused by the environmental bacterium Legionella pneumophila. Outbreaks commonly affect people with known risk factors, but the genetic and pathogenic complexity of L. pneumophila within an outbreak is not well understood. Here, we investigate the etiology of the major Legionnaires' disease outbreak that occurred in Edinburgh, UK, in 2012, by examining the evolutionary history, genome content, and virulence of L. pneumophila clinical isolates. RESULTS: Our high resolution genomic approach reveals that the outbreak was caused by multiple genetic subtypes of L. pneumophila, the majority of which had diversified from a single progenitor through mutation, recombination, and horizontal gene transfer within an environmental reservoir prior to release. In addition, we discover that some patients were infected with multiple L. pneumophila subtypes, a finding which can affect the certainty of source attribution. Importantly, variation in the complement of type IV secretion systems encoded by different genetic subtypes correlates with virulence in a Galleria mellonella model of infection, revealing variation in pathogenic potential among the outbreak source population of L. pneumophila. CONCLUSIONS: Taken together, our study indicates previously cryptic levels of pathogen heterogeneity within a Legionnaires' disease outbreak, a discovery that impacts on source attribution for future outbreak investigations. Furthermore, our data suggest that in addition to host immune status, pathogen diversity may be an important influence on the clinical outcome of individual outbreak infections.

Molecular longitudinal tracking of Mycobacterium abscessus spp. during chronic infection of the human lung.

The Mycobacterium abscessus complex is an emerging cause of chronic pulmonary infection in patients with underlying lung disease. The M. abscessus complex is regarded as an environmental pathogen but its molecular adaptation to the human lung during long-term infection is poorly understood. Here we carried out a longitudinal molecular epidemiological analysis of 178 M. abscessus spp. isolates obtained from 10 cystic fibrosis (CF) and 2 non CF patients over a 13 year period. Multi-locus sequence and molecular typing analysis revealed that 11 of 12 patients were persistently colonized with the same genotype during the course of the infection while replacement of a M. abscessus sensu stricto strain with a Mycobacterium massiliense strain was observed for a single patient. Of note, several patients including a pair of siblings were colonized with closely-related strains consistent with intra-familial transmission or a common infection reservoir. In general, a switch from smooth to rough colony morphology was observed during the course of long-term infection, which in some cases correlated with an increasing severity of clinical symptoms. To examine evolution during long-term infection of the CF lung we compared the genome sequences of 6 sequential isolates of Mycobacterium bolletii obtained from a single patient over an 11 year period, revealing a heterogeneous clonal infecting population with mutations in regulators controlling the expression of virulence factors and complex lipids. Taken together, these data provide new insights into the epidemiology of M. abscessus spp. during long-term infection of the CF lung, and the molecular transition from saprophytic organism to human pathogen.

Livestock origin for a human pandemic clone of community-associated methicillin-resistant Staphylococcus aureus.

UNLABELLED: The importance of livestock as a source of bacterial pathogens with the potential for epidemic spread in human populations is unclear. In recent years, there has been a global increase in community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) infections of healthy humans, but an understanding of the different evolutionary origins of CA-MRSA clones and the basis for their recent expansion is lacking. Here, using a high-resolution phylogenetic approach, we report the discovery of two emergent clones of human epidemic CA-MRSA which resulted from independent livestock-to-human host jumps by the major bovine S. aureus complex, CC97. Of note, one of the new clones was isolated from human infections on four continents, demonstrating its global dissemination since the host jump occurred over 40 years ago. The emergence of both human S. aureus clones coincided with the independent acquisition of mobile genetic elements encoding antimicrobial resistance and human-specific mediators of immune evasion, consistent with an important role for these genetic events in the capacity to survive and transmit among human populations. In conclusion, we provide evidence that livestock represent a reservoir for the emergence of new human-pathogenic S. aureus clones with the capacity for pandemic spread. These findings have major public health implications highlighting the importance of surveillance for early identification of emergent clones and improved transmission control measures at the human-livestock interface. IMPORTANCE: Animals are the major source of new pathogens affecting humans. However, the potential for pathogenic bacteria that originally were found in animals to switch hosts and become widely established in human populations is not clear. Here, we report the discovery of emergent clones of methicillin-resistant Staphylococcus aureus (MRSA) that originated in livestock and switched to humans, followed by host-adaptive evolution and epidemic spread in global human populations. Our findings demonstrate that livestock can act as a reservoir for the emergence of new human bacterial clones with potential for pandemic spread, highlighting the potential role of surveillance and biosecurity measures in the agricultural setting for preventing the emergence of new human pathogens.

Adaptive evolution of Staphylococcus aureus during chronic endobronchial infection of a cystic fibrosis patient.

The molecular adaptation of Staphylococcus aureus to its host during chronic infection is not well understood. Comparative genome sequencing of 3 S. aureus isolates obtained sequentially over 26 months from the airways of a cystic fibrosis patient, revealed variation in phage content, and genetic polymorphisms in genes which influence antibiotic resistance, and global regulation of virulence. The majority of polymorphisms were isolate-specific suggesting the existence of an heterogeneous infecting population that evolved from a single infecting strain of S. aureus. The genetic variation identified correlated with differences in growth rate, hemolytic activity, and antibiotic sensitivity, implying a profound effect on the ecology of S. aureus. In particular, a high frequency of mutations in loci associated with the alternate transcription factor SigB, were observed. The identification of genes under diversifying selection during long-term infection may inform the design of novel therapeutics for the control of refractory chronic infections.

Improving reverse vaccinology with a machine learning approach.

Reverse vaccinology aims to accelerate subunit vaccine design by rapidly predicting which proteins in a pathogenic bacterial proteome are putative protective antigens. Support vector machine classification is a machine learning approach that has been applied to solve numerous classification problems in biological sciences but has not previously been incorporated into a reverse vaccinology approach. A training data set of 136 bacterial protective antigens paired with 136 non-antigens was constructed and bioinformatic tools were used to annotate this data for predicted protein features, many of which are associated with antigenicity (i.e. extracellular localization, signal peptides and B-cell epitopes). Annotation was used to train support vector machine classifiers that exhibited a maximum accuracy of 92% for discriminating protective antigens from non-antigens as assessed by a leave-tenth-out cross-validation approach. These accuracies were superior to those achieved when annotating training data with auto and cross covariance transformations of z-descriptors for hydrophobicity, molecular size and polarity, or when classification was performed using regression methods. To further validate support vector machine classifiers, they were used to rank all the proteins in six bacterial proteomes for their antigenicity. Protective antigens from the training data were significantly recalled (enriched) in the top 75 ranked proteins for all six proteomes as assessed by a Fisher's exact test (p<0.05). This paper describes a superior workflow for performing reverse vaccinology studies and provides a benchmark training data set that can be used to evaluate future methodological improvements.