Building a genomic framework for prospective MRSA surveillance in the United Kingdom and the Republic of Ireland.

The correct interpretation of microbial sequencing data applied to surveillance and outbreak investigation depends on accessible genomic databases to provide vital genetic context. Our aim was to construct and describe a United Kingdom MRSA database containing over 1000 methicillin-resistant Staphylococcus aureus (MRSA) genomes drawn from England, Northern Ireland, Wales, Scotland, and the Republic of Ireland over a decade. We sequenced 1013 MRSA submitted to the British Society for Antimicrobial Chemotherapy by 46 laboratories between 2001 and 2010. Each isolate was assigned to a regional healthcare referral network in England and was otherwise grouped based on country of origin. Phylogenetic reconstructions were used to contextualize MRSA outbreak investigations and to detect the spread of resistance. The majority of isolates (n = 783, 77%) belonged to CC22, which contains the dominant United Kingdom epidemic clone (EMRSA-15). There was marked geographic structuring of EMRSA-15, consistent with widespread dissemination prior to the sampling decade followed by local diversification. The addition of MRSA genomes from two outbreaks and one pseudo-outbreak demonstrated the certainty with which outbreaks could be confirmed or refuted. We identified local and regional differences in antibiotic resistance profiles, with examples of local expansion, as well as widespread circulation of mobile genetic elements across the bacterial population. We have generated a resource for the future surveillance and outbreak investigation of MRSA in the United Kingdom and Ireland and have shown the value of this during outbreak investigation and tracking of antimicrobial resistance.

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.

Endemic melioidosis in residents of desert region after atypically intense rainfall in central Australia, 2011.

After heavy rains and flooding during early 2011 in the normally arid interior of Australia, melioidosis was diagnosed in 6 persons over a 4-month period. Although the precise global distribution of the causal bacterium Burkholderia pseudomallei remains to be determined, this organism can clearly survive in harsh and even desert environments outside the wet tropics.

Whole-genome sequencing confirms that Burkholderia pseudomallei multilocus sequence types common to both Cambodia and Australia are due to homoplasy.

Burkholderia pseudomallei isolates with shared multilocus sequence types (STs) have not been isolated from different continents. We identified two STs shared between Australia and Cambodia. Whole-genome analysis revealed substantial diversity within STs, correctly identified the Asian or Australian origin, and confirmed that these shared STs were due to homoplasy.

Patient sharing and population genetic structure of methicillin-resistant Staphylococcus aureus.

Rates of hospital-acquired infections, specifically methicillin-resistant Staphylococcus aureus (MRSA), are increasingly being used as indicators for quality of hospital hygiene. There has been much effort on understanding the transmission process at the hospital level; however, interhospital population-based transmission remains poorly defined. We evaluated whether the proportion of shared patients between hospitals was correlated with genetic similarity of MRSA strains from those hospitals. Using data collected from 30 of 32 hospitals in Orange County, California, multivariate linear regression showed that for each twofold increase in the proportion of patients shared between 2 hospitals, there was a 7.7% reduction in genetic heterogeneity between the hospitals' MRSA populations (permutation P value = 0.0356). Pairs of hospitals that both served adults had more similar MRSA populations than pairs including a pediatric hospital. These findings suggest that concerted efforts among hospitals that share large numbers of patients may be synergistic to prevent MRSA transmission.

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.

Distribution of Burkholderia pseudomallei in northern Australia, a land of diversity.

Burkholderia pseudomallei is a Gram-negative soil bacillus that is the etiological agent of melioidosis and a biothreat agent. Little is known about the biogeography of this bacterium in Australia, despite its hyperendemicity in the northern region of this continent. The population structure of 953 Australian B. pseudomallei strains representing 779 and 174 isolates of clinical and environmental origins, respectively, was analyzed using multilocus sequence typing (MLST). Bayesian population structure and network SplitsTree analyses were performed on concatenated MLST loci, and sequence type (ST) diversity and evenness were examined using Simpson's and Pielou's indices and a multivariate dissimilarity matrix. Bayesian analysis found two B. pseudomallei populations in Australia that were geographically distinct; isolates from the Northern Territory were grouped mainly into the first population, whereas the majority of isolates from Queensland were grouped in a second population. Differences in ST evenness were observed between sampling areas, confirming that B. pseudomallei is widespread and established across northern Australia, with a large number of fragmented habitats. ST analysis showed that B. pseudomallei populations diversified as the sampling area increased. This observation was in contrast to smaller sampling areas where a few STs predominated, suggesting that B. pseudomallei populations are ecologically established and not frequently dispersed. Interestingly, there was no identifiable ST bias between clinical and environmental isolates, suggesting the potential for all culturable B. pseudomallei isolates to cause disease. Our findings have important implications for understanding the ecology of B. pseudomallei in Australia and for potential source attribution of this bacterium in the event of unexpected cases of melioidosis.

Diversity of methicillin-resistant Staphylococcus aureus strains isolated from residents of 26 nursing homes in Orange County, California.

Nursing homes represent a unique and important methicillin-resistant Staphylococcus aureus (MRSA) reservoir. Not only are strains imported from hospitals and the community, strains can be transported back into these settings from nursing homes. Since MRSA bacteria are prevalent in nursing homes and yet relatively poorly studied in this setting, a multicenter, regional assessment of the frequency and diversity of MRSA in the nursing home reservoir was carried out and compared to that of the MRSA from hospitals in the same region. The prospective study collected MRSA from nasal swabbing of residents of 26 nursing homes in Orange County, California, and characterized each isolate by spa typing. A total of 837 MRSA isolates were collected from the nursing homes. Estimates of admission prevalence and point prevalence of MRSA were 16% and 26%, respectively. The spa type genetic diversity was heterogeneous between nursing homes and significantly higher overall (77%) than the diversity in Orange County hospitals (72%). MRSA burden in nursing homes appears largely due to importation from hospitals. As seen in Orange County hospitals, USA300 (sequence type 8 [ST8]/t008), USA100 (ST5/t002), and a USA100 variant (ST5/t242) were the dominant MRSA clones in Orange County nursing homes, representing 83% of all isolates, although the USA100 variant was predominant in nursing homes, whereas USA300 was predominant in hospitals. Control strategies tailored to the complex problem of MRSA transmission and infection in nursing homes are needed in order to minimize the impact of this unique reservoir on the overall regional MRSA burden.

Differences in methicillin-resistant Staphylococcus aureus strains isolated from pediatric and adult patients from hospitals in a large county in California.

Studies of U.S. epidemics of community- and health care-associated methicillin-resistant Staphylococcus aureus (MRSA) suggested differences in MRSA strains in adults and those in children. Comprehensive population-based studies exploring these differences are lacking. We conducted a prospective cohort study of inpatients in Orange County, CA, collecting clinical MRSA isolates from 30 of 31 Orange County hospitals, to characterize differences in MRSA strains isolated from children compared to those isolated from adults. All isolates were characterized by spa typing. We collected 1,124 MRSA isolates from adults and 159 from children. Annual Orange County population estimates of MRSA inpatient clinical cultures were 119/100,000 adults and 22/100,000 children. spa types t008, t242, and t002 accounted for 83% of all isolates. The distribution of these three spa types among adults was significantly different from that among children (χ(2) = 52.29; P < 0.001). Forty-one percent of adult isolates were of t008 (USA300), compared to 69% of pediatric isolates. In multivariate analyses, specimens from pediatric patients, wounds, non-intensive care unit (ICU) wards, and hospitals with a high proportion of Medicaid-insured patients were significantly associated with the detection of t008 strains. While community- and health care-associated MRSA reservoirs have begun to merge, significant differences remain in pediatric and adult patient populations. Community-associated MRSA spa type t008 is significantly more common in pediatric patients.

The 2011 Garrod Lecture: From penicillin-binding proteins to molecular epidemiology.

In this review, based on my Garrod Lecture to the British Society for Antimicrobial Chemotherapy, I have given a brief outline of my career over the past 40 years, starting with research in the 1970s into the properties and functions of penicillin-binding proteins (PBPs), leading to the identification of the high molecular mass PBPs as the physiological targets of penicillin, and subsequent studies showing the emergence of low-affinity PBPs in penicillin-resistant clinical isolates by inter-species recombination and the generation of mosaic PBP genes. The studies of clinical isolates of gonococci, meningococci and pneumococci with PBP-mediated resistance to penicillin led to new interests in molecular epidemiology and the population and evolutionary biology of bacterial pathogens. The development (with colleagues) of multilocus sequence typing provided a method for the unambiguous characterization of bacterial strains that has proved to be very widely used, but the recent remarkable (and ongoing) developments in DNA sequencing technologies have provided the prospect of being able routinely to use whole genome sequences to characterize pathogen isolates. These developments will soon have major implications for diagnostic microbiology, outbreak investigations and our ability to follow the spread of strains of community-acquired and nosocomial pathogens at local, national and international levels. However, there are major barriers to be overcome, particularly with respect to how the avalanche of genome sequence data will be stored so that its transformative potential for molecular epidemiology and international public health are fully realized.

Bayesian semi-supervised classification of bacterial samples using MLST databases.

BACKGROUND: Worldwide effort on sampling and characterization of molecular variation within a large number of human and animal pathogens has lead to the emergence of multi-locus sequence typing (MLST) databases as an important tool for studying the epidemiology and evolution of pathogens. Many of these databases are currently harboring several thousands of multi-locus DNA sequence types (STs) enriched with metadata over traits such as serotype, antibiotic resistance, host organism etc of the isolates. Curators of the databases have thus the possibility of dividing the pathogen populations into subsets representing different evolutionary lineages, geographically associated groups, or other subpopulations, which are defined in terms of molecular similarities and dissimilarities residing within a database. When combined with the existing metadata, such subsets may provide invaluable information for assessing the position of a new set of isolates in relation to the whole pathogen population. RESULTS: To enable users of MLST schemes to query the databases with sets of new bacterial isolates and to automatically analyze their relation to existing curated sequences, we introduce here a Bayesian model-based method for semi-supervised classification of MLST data. Our method can use an MLST database as a training set and assign simultaneously any set of query sequences into the earlier discovered lineages/populations, while also allowing some or all of these sequences to form previously undiscovered genetically distinct groups. This tool provides probabilistic quantification of the classification uncertainty and is highly efficient computationally, thus enabling rapid analyses of large databases and sets of query sequences. The latter feature is a necessary prerequisite for an automated access through the MLST web interface. We demonstrate the versatility of our approach by anayzing both real and synthesized data from MLST databases. The introduced method for semi-supervised classification of sets of query STs is freely available for Windows, Mac OS X and Linux operative systems in BAPS 5.4 software which is downloadable at http://web.abo.fi/fak/mnf/mate/jc/software/baps.html. The query functionality is also directly available for the Staphylococcus aureus database at http://www.mlst.net and shortly will be available for other species databases hosted at this web portal. CONCLUSIONS: We have introduced a model-based tool for automated semi-supervised classification of new pathogen samples that can be integrated into the web interface of the MLST databases. In particular, when combined with the existing metadata, the semi-supervised labeling may provide invaluable information for assessing the position of a new set of query strains in relation to the particular pathogen population represented by the curated database.Such information will be useful both for clinical and basic research purposes.

Burkholderia pseudomallei in unchlorinated domestic bore water, Tropical Northern Australia.

To determine whether unchlorinated bore water in northern Australia contained Burkholderia pseudomallei organisms, we sampled 55 bores; 18 (33%) were culture positive. Multilocus sequence typing identified 15 sequence types. The B. pseudomallei sequence type from 1 water sample matched a clinical isolate from a resident with melioidosis on the same property.

Analysis of emergence of quinolone-resistant gonococci in Greece by combined use of Neisseria gonorrhoeae multiantigen sequence typing and multilocus sequence typing.

The prevalence of quinolone-resistant Neisseria gonorrhoeae (QRNG) in Greece remained low from 1997 to 2003 but increased dramatically from 11% to 56% between 2004 and 2007. N. gonorrhoeae multiantigen sequence typing (NG-MAST) and multilocus sequence typing (MLST) were used to investigate trends in quinolone resistance from 1997 to 2007 and explore the origins of the recent increase in QRNG. We characterized 295 QRNG isolates from the study period and 233 quinolone-susceptible (QS) gonococci from 2004 and 2005, when the rapid increase in QRNG occurred. From 1997 to 1999, an outbreak of QRNG was due to the dissemination of isolates of serovar Arst that belonged to two closely related genotypes. Few QRNG isolates, of diverse genotypes, were present between 2001 and 2003, whereas the sharp increase in QRNG from 2004 onwards was due to the appearance of serovar Bropyst isolates of several major NG-MAST sequence type (STs) that previously had not been identified in Greece. These isolates were shown by MLST to be variants of a single multiply antibiotic-resistant QRNG strain (ST1901) that appeared in Greece and rapidly diversified into 31 NG-MAST STs. There were no isolates of MLST ST1901 or any of the 31 NG-MAST STs among QS isolates from 2004 and 2005 or among 8 representatives of multiresistant but quinolone-susceptible serovar Bropyst isolates circulating in Greece during the 1990 s, supporting the view that the recent increase in QRNG was due to importation of a QRNG strain(s) of MLST ST1901 into Greece. Recently, multiresistant QRNG isolates of ST1901 with reduced susceptibility to the newer cephalosporins have appeared in Greece.

Genomic evidence for the evolution of Streptococcus equi: host restriction, increased virulence, and genetic exchange with human pathogens.

The continued evolution of bacterial pathogens has major implications for both human and animal disease, but the exchange of genetic material between host-restricted pathogens is rarely considered. Streptococcus equi subspecies equi (S. equi) is a host-restricted pathogen of horses that has evolved from the zoonotic pathogen Streptococcus equi subspecies zooepidemicus (S. zooepidemicus). These pathogens share approximately 80% genome sequence identity with the important human pathogen Streptococcus pyogenes. We sequenced and compared the genomes of S. equi 4047 and S. zooepidemicus H70 and screened S. equi and S. zooepidemicus strains from around the world to uncover evidence of the genetic events that have shaped the evolution of the S. equi genome and led to its emergence as a host-restricted pathogen. Our analysis provides evidence of functional loss due to mutation and deletion, coupled with pathogenic specialization through the acquisition of bacteriophage encoding a phospholipase A(2) toxin, and four superantigens, and an integrative conjugative element carrying a novel iron acquisition system with similarity to the high pathogenicity island of Yersinia pestis. We also highlight that S. equi, S. zooepidemicus, and S. pyogenes share a common phage pool that enhances cross-species pathogen evolution. We conclude that the complex interplay of functional loss, pathogenic specialization, and genetic exchange between S. equi, S. zooepidemicus, and S. pyogenes continues to influence the evolution of these important streptococci.

Poor clinical outcome for meningitis caused by Haemophilus influenzae serotype A strains containing the IS1016-bexA deletion.

Since the introduction of Haemophilus influenzae type b (Hib) conjugate vaccines, meningitis caused by serotypes other than Hib has gained in importance. We conducted active hospital-based surveillance for meningitis over an 11-year period in Salvador, Brazil. H. influenzae isolates were serotyped and analyzed by polymerase chain reaction, pulsed-field gel electrophoresis, and DNA sequencing to identify strains with a specific deletion (IS1016) in the bexA gene (IS1016-bexA). We identified 43 meningitis cases caused by non-type b H. influenzae: 28 (65%) were caused by type a (Hia), 9 (21%) were caused by noncapsulated strains, and 3 (7%) each were caused by types e and f. Hia isolates clustered in 2 clonal groups; clonal group A strains (n = 9) had the IS1016-bexA deletion. Among children <5 years of age, meningitis caused by Hia from clonal group A had higher case-fatality than meningitis caused by clonal group B. Despite small numbers, these results indicate that the presence of the IS1016-bexA deletion is associated with enhanced virulence in non-type b H. influenzae.

Geographic distribution of Staphylococcus aureus causing invasive infections in Europe: a molecular-epidemiological analysis.

BACKGROUND: Staphylococcus aureus is one of the most important human pathogens and methicillin-resistant variants (MRSAs) are a major cause of hospital and community-acquired infection. We aimed to map the geographic distribution of the dominant clones that cause invasive infections in Europe. METHODS AND FINDINGS: In each country, staphylococcal reference laboratories secured the participation of a sufficient number of hospital laboratories to achieve national geo-demographic representation. Participating laboratories collected successive methicillin-susceptible (MSSA) and MRSA isolates from patients with invasive S. aureus infection using an agreed protocol. All isolates were sent to the respective national reference laboratories and characterised by quality-controlled sequence typing of the variable region of the staphylococcal spa gene (spa typing), and data were uploaded to a central database. Relevant genetic and phenotypic information was assembled for interactive interrogation by a purpose-built Web-based mapping application. Between September 2006 and February 2007, 357 laboratories serving 450 hospitals in 26 countries collected 2,890 MSSA and MRSA isolates from patients with invasive S. aureus infection. A wide geographical distribution of spa types was found with some prevalent in all European countries. MSSA were more diverse than MRSA. Genetic diversity of MRSA differed considerably between countries with dominant MRSA spa types forming distinctive geographical clusters. We provide evidence that a network approach consisting of decentralised typing and visualisation of aggregated data using an interactive mapping tool can provide important information on the dynamics of MRSA populations such as early signalling of emerging strains, cross border spread, and importation by travel. CONCLUSIONS: In contrast to MSSA, MRSA spa types have a predominantly regional distribution in Europe. This finding is indicative of the selection and spread of a limited number of clones within health care networks, suggesting that control efforts aimed at interrupting the spread within and between health care institutions may not only be feasible but ultimately successful and should therefore be strongly encouraged.

Assigning strains to bacterial species via the internet.

BACKGROUND: Methods for assigning strains to bacterial species are cumbersome and no longer fit for purpose. The concatenated sequences of multiple house-keeping genes have been shown to be able to define and circumscribe bacterial species as sequence clusters. The advantage of this approach (multilocus sequence analysis; MLSA) is that, for any group of related species, a strain database can be produced and combined with software that allows query strains to be assigned to species via the internet. As an exemplar of this approach, we have studied a group of species, the viridans streptococci, which are very difficult to assign to species using standard taxonomic procedures, and have developed a website that allows species assignment via the internet. RESULTS: Seven house-keeping gene sequences were obtained from 420 streptococcal strains to produce a viridans group database. The reference tree produced using the concatenated sequences identified sequence clusters which, by examining the position on the tree of the type strain of each viridans group species, could be equated with species clusters. MLSA also identified clusters that may correspond to new species, and previously described species whose status needs to be re-examined. A generic website and software for electronic taxonomy was developed. This site http://www.eMLSA.net allows the sequences of the seven gene fragments of a query strain to be entered and for the species assignment to be returned, according to its position within an assigned species cluster on the reference tree. CONCLUSION: The MLSA approach resulted in the identification of well-resolved species clusters within this taxonomically challenging group and, using the software we have developed, allows unknown strains to be assigned to viridans species via the internet. Submission of new strains will provide a growing resource for the taxonomy of viridans group streptococci, allowing the recognition of potential new species and taxonomic anomalies. More generally, as the software at the MLSA website is generic, MLSA schemes and strain databases for other groups of related species can be hosted at this website, providing a portal for microbial electronic taxonomy.

Relative contributions of recombination and mutation to the diversification of the opa gene repertoire of Neisseria gonorrhoeae.

To understand the rates and mechanisms of Neisseria gonorrhoeae opa gene variation, the 11 opa genes were amplified independently so that an opa allelic profile could be defined for any isolate from the sequences at each locus. The opa allelic profiles from 14 unrelated isolates were all different, with no opa alleles shared between isolates. Examination of very closely related isolates from sexual contacts and sexual networks showed that these typically shared most opa alleles, and the mechanisms by which recent changes occurred at individual opa loci could be determined. The great majority of changes were due to recombination among existing alleles that duplicated an opa allele present at another locus or resulted in a mosaic of existing opa alleles. Single nucleotide changes or insertion/deletion of a single codon also occurred, but few of these events were assigned to mutation, the majority being assigned to localized recombination. Introduction of novel opa genes from coinfecting strains was rare, and all but one were observed in the same sexual network. Changes at opa loci occurred at a greater rate than those at the porin locus, and the opa11 locus changed more rapidly than other opa loci, almost always differing even between recent sexual contacts. Examination of the neighboring pilE gene showed that changes at opa11 and pilE often occurred together, although this linkage may not be a causal one.