Phylogeographic variation in recombination rates within a global clone of methicillin-resistant Staphylococcus aureus

Background

Next-generation sequencing (NGS) is a powerful tool for understanding both patterns of descent over time and space (phylogeography) and the molecular processes underpinning genome divergence in pathogenic bacteria. Here, we describe a synthesis between these perspectives by employing a recently developed Bayesian approach, BRATNextGen, for detecting recombination on an expanded NGS dataset of the globally disseminated methicillin-resistant Staphylococcus aureus (MRSA) clone ST239.

Results

The data confirm strong geographical clustering at continental, national and city scales and demonstrate that the rate of recombination varies significantly between phylogeographic sub-groups representing independent introductions from Europe. These differences are most striking when mobile non-core genes are included, but remain apparent even when only considering the stable core genome. The monophyletic ST239 sub-group corresponding to isolates from South America shows heightened recombination, the sub-group predominantly from Asia shows an intermediate level, and a very low level of recombination is noted in a third sub-group representing a large collection from Turkey.

Conclusions

We show that the rapid global dissemination of a single pathogenic bacterial clone results in local variation in measured recombination rates. Possible explanatory variables include the size and time since emergence of each defined sub-population (as determined by the sampling frame), variation in transmission dynamics due to host movement, and changes in the bacterial genome affecting the propensity for recombination.     

Unifying Viral Genetics and Human Transportation Data to Predict the Global Transmission Dynamics of Human Influenza H3N2

Information on global human movement patterns is central to spatial epidemiological models used to predict the behavior of influenza and other infectious diseases. Yet it remains difficult to test which modes of dispersal drive pathogen spread at various geographic scales using standard epidemiological data alone. Evolutionary analyses of pathogen genome sequences increasingly provide insights into the spatial dynamics of influenza viruses, but to date they have largely neglected the wealth of information on human mobility, mainly because no statistical framework exists within which viral gene sequences and empirical data on host movement can be combined. Here, we address this problem by applying a phylogeographic approach to elucidate the global spread of human influenza subtype H3N2 and assess its ability to predict the spatial spread of human influenza A viruses worldwide. Using a framework that estimates the migration history of human influenza while simultaneously testing and quantifying a range of potential predictive variables of spatial spread, we show that the global dynamics of influenza H3N2 are driven by air passenger flows, whereas at more local scales spread is also determined by processes that correlate with geographic distance. Our analyses further confirm a central role for mainland China and Southeast Asia in maintaining a source population for global influenza diversity. By comparing model output with the known pandemic expansion of H1N1 during 2009, we demonstrate that predictions of influenza spatial spread are most accurate when data on human mobility and viral evolution are integrated. In conclusion, the global dynamics of influenza viruses are best explained by combining human mobility data with the spatial information inherent in sampled viral genomes. The integrated approach introduced here offers great potential for epidemiological surveillance through phylogeographic reconstructions and for improving predictive models of disease control.     

Phylogeographic characteristics of hypervariable regions in the mitochondrial genome of a cosmopolitan,bloom‐forming raphidophyte,Heterosigma akashiwo

The global distribution of phytoplankton is defined by many events, including long‐term evolutionary processes and shorter time span processes (e.g., global climate change). Furthermore, human‐assisted, unintentional dispersion, including the transport of live fish and spat for aquaculture, and transfer of aquatic microorganisms contained in ship ballast water, may aid the spread of phytoplankton. To understand the phylogeographic history of a species, the development of useful molecular markers is crucial. We previously reported a hypervariable mitochondrial gene in the cosmopolitan bloom‐forming alga, Heterosigma akashiwo. In this study, we identified two additional hypervariable segments in the H. akashiwo mitochondrial genome, one a protein coding sequence, and the other an intergenic region, by comparing the whole mitogenomes of strains obtained from various geographic origins. Interestingly, the newly identified hypervariable protein coding sequence was a paralogue of the previously identified gene, and both sequences showed tendencies to correlate with latitude of geographic origin. However, the hypervariable intergenic sequence did not show a clear correlation with origin. Our results demonstrated that the protein coding sequences may serve as useful tools for understanding the phylogeographic history of H. akashiwo, and they may crucially function in adaptation of the species to the environment.     

MRSA epidemic linked to a quickly spreading colonization and virulence determinant

The molecular processes underlying epidemic waves of methicillin-resistant Staphylococcus aureus (MRSA) infection are poorly understood(1). Although a major role has been attributed to the acquisition of virulence determinants by horizontal gene transfer(2), there are insufficient epidemiological and functional data supporting that concept. We here report the spread of clones containing a previously extremely rare(3,4) mobile genetic element–encoded gene, sasX. We demonstrate that sasX has a key role in MRSA colonization and pathogenesis, substantially enhancing nasal colonization, lung disease and abscess formation and promoting mechanisms of immune evasion. Moreover, we observed the recent spread of sasX from sequence type 239 (ST239) to invasive clones belonging to other sequence types. Our study identifies sasX as a quickly spreading crucial determinant of MRSA pathogenic success and a promising target for therapeutic interference. Our results provide proof of principle that horizontal gene transfer of key virulence determinants drives MRSA epidemic waves.     

Complete genome sequence of Staphylococcus aureus T0131, an ST239-MRSA-SCCmec type III clone isolated in China

We report here the complete genome sequence of Staphylococcus aureus T0131, which is a multiresistant clinical isolate recovered in China and the first sequenced epidemic ST239-MRSA-SCCmec type III strain obtained in Asia. Comparison with two published genomes of ST239 reveals the polymorphism among strains of this type from different continents.     

Molecular analysis and antimicrobial susceptibility of methicillin resistant Staphylococcus aureus in one of the hospitals of Tehran University of Medical Sciences: high prevalence of sequence type 239 (ST239) clone

Methicillin-resistant Staphylococcus aureus (MRSA), particularly the multidrug-resistant clones, is an increasing worldwide problem. The average incidence rate of MRSA in Tehran was found to be over 40%. A total of 140 MRSA isolates obtained from patients attending a teaching hospital in Tehran, from May 2009 to December 2009, were included in this study. The antimicrobial susceptibility profile of MRSA isolates was determined by the agar disk diffusion method. Molecular analysis of MRSA strains was accomplished by Pulsed-Field Gel Electrophoresis (PFGE) and Multi-locus sequence typing (MLST). Detection of mecA gene was used to confirm resistance to methicillin among the MRSA isolates. All the MRSA isolates were susceptible to chloramphenicol, teicoplanin, tigecycline and vancomycin. All MRSAisolates were resistant to oxacillin, whilst 139 strains showed resistance against ciprofloxacin, erythromycin, gentamicin, tetracycline and trimethoprim-sulfamethoxazole. PFGE analysis of all the 140 MRSA isolates produced five distinct pulsotypes designated as pulsotypes A-E. Most of the isolates (n=132) were clustered into pulsotype A. The most prevalent sequence type (ST) was ST 239 (pulsotype A) found in 82% (37/45) of the tested isolates. The second most prevalent type was ST 1238 (pulsotypes B, C and D) found in 15% (7/45) of the isolates. The remaining type, ST 8 (pulsotype E) was found in a single isolate. The results of this study indicated that the MRSA clone ST 239 was a major clone in the selected university hospital of Tehran and that it was widely spread among the different wards as well as all the age groups of patients.     

Panton-Valentine leukocidin gene sequence variation and phage in methicillin-resistant and methicillin-susceptible Staphylococcus aureus from children in mainland China

To determine the variation in the Panton–Valentine leukocidin (PVL) gene sequences and different PVL‐encoding phages of Staphylococcus aureus strains collected from children in mainland China, fifty‐eight strains with PVL collected from 2007 to 2009 were used. Their molecular characteristics were examined. Primers were designed to sequence the PVL genes. Six PVL‐encoding phages (?PVL, ?108PVL, ?SLT, ?Sa2MW, ?Sa2958, and ?Sa2USA) were identified by PCR. Eleven sequence types (ST) were detected with ST59 (39.7%, 23/58) the most frequent ST, followed by 910 (22.4%, 13/58), and 338 (12.1%, 7/58). Single nucleotide polymorphisms (SNP) were identified at 11 locations in the PVL genes. SNP (nucleotide 1396, A→G) and SNP (nucleotide 1546, A→G) were observed in >10 sequences. Four additional SNP were non‐synonymous. Both SNP (nucleotide 16, C→A) and SNP (nucleotide 62, C→T) were present in the same ST59 strain. SNP (nucleotide 527, A→G) was present in five strains belonging to ST30, 121, 1, and 93. SNP (nucleotide 1436, A→C) was present in one ST30 strain. Fifteen strains belonging to ST910, ST217, and ST30 carried a PVL phage that had an icosahedral head morphology. Nine ST59 strains carried ?108PVL. Three ST88 strains carried a PVL phage that had an elongated head morphology. Twenty‐seven strains, including 60.9% (14/23) of ST59 and all ST338 strains, had no detectable phage. In conclusion, sequence variation in PVL genes and PVL‐encoding phages was generally related to the lineage. ST59 strains may indeed carry novel PVL phages.     

Multilocus Sequence Analysis of Xanthomonads Causing Bacterial Spot of Tomato and Pepper Plants Reveals Strains Generated by Recombination among Species and Recent Global Spread of Xanthomonas gardneri

Four Xanthomonas species are known to cause bacterial spot of tomato and pepper, but the global distribution and genetic diversity of these species are not well understood. A collection of bacterial spot-causing strains from the Americas, Africa, Southeast Asia, and New Zealand were characterized for genetic diversity and phylogenetic relationships using multilocus sequence analysis of six housekeeping genes. By examining strains from different continents, we found unexpected phylogeographic patterns, including the global distribution of a single multilocus haplotype of X. gardneri, possible regional differentiation in X. vesicatoria, and high species diversity on tomato in Africa. In addition, we found evidence of multiple recombination events between X. euvesicatoria and X. perforans. Our results indicate that there have been shifts in the species composition of bacterial spot pathogen populations due to the global spread of dominant genotypes and that recombination between species has generated genetic diversity in these populations.     

Comparative genomics and drug resistance of a geographic variant of ST239 methicillin-resistant Staphylococcus aureus emerged in Russia

Two distinct classes of methicillin-resistant Staphylococcus aureus (MRSA) are spreading in hospitals (as hospital-acquired MRSA, HA-MRSA) and in the community (as community-acquired MRSA, CA-MRSA). Multilocus sequence type (ST) 239 MRSA, one of the most worldwide-disseminated lineages, has been noted as a representative HA-MRSA. Here, we isolated ST239 MRSA (spa type 3 [t037] and staphylococcal cassette chromosome mec [SCCmec] type III.1.1.1) and its novel variant with ST239/spa351 (t030)/SCCmecIII.1.1.4 (SCCmecIII(R)) not only from hospitals but also from patients with urethritis in the community in Russia. The Russian variant (strain 16K) possessed a hybrid genome consisting of CC8 and CC30, similar to the ST239/spa3/SCCmecIII.1.1.1 HA-MRSA (TW20) genome, but with marked diversity. The 16K' CC30 section had SCCmecIII(R) carrying the dcs-carrying unit (which corresponded to the SCCmecIVc J3 joining region of ST30 CA-MRSA), lacked SCCmercury, and possessed a novel mobile element structure (MES16K) carrying the ccrC-carrying unit (with the recombinase gene ccrC1 allele 3) and drug resistance tranposons. The Russian variant included strains with a high ability to transfer its multiple drug resistance by conjugation; e.g., for strain 16K, the transfer frequency of a chloramphenicol resistance plasmid (p16K-1 with 2.9 kb in size) reached 1.4×10(-2), followed by Tn554 conjugative transfer at 3.6×l0(-4). The Russian variant, which has been increasing recently, included divergent strains with different plasmid patterns and pulsed field gel electrophoresis profiles. The data demonstrate the alternative nature of ST239 MRSA as CA-MRSA and also as a drug resistance disseminator, and its micro but dynamic evolution in Russia.     

Full genome SNP-based phylogenetic analysis reveals the origin and global spread of Brucella melitensis

Background

Brucellosis is an important zoonotic disease that affects both humans and animals. We sequenced the full genome and characterised the genetic diversity of two Brucella melitensis isolates from Malaysia and the Philippines. In addition, we performed a comparative whole-genome single nucleotide polymorphism (SNP) analysis of B. melitensis strains collected from around the world, to investigate the potential origin and the history of the global spread of B. melitensis.

Results

Single sequencing runs of each genome resulted in draft genome sequences of MY1483/09 and Phil1136/12, which covered 99.85% and 99.92% of the complete genome sequences, respectively. The B. melitensis genome sequences, and two B. abortus strains used as the outgroup strains, yielded a total of 13,728 SNP sites. Phylogenetic analysis using whole-genome SNPs and geographical distribution of the isolates revealed spatial clustering of the B. melitensis isolates into five genotypes, I, II, III, IV and V. The Mediterranean strains, identified as genotype I, occupied the basal node of the phylogenetic tree, suggesting that B. melitensis may have originated from the Mediterranean regions. All of the Asian B. melitensis strains clustered into genotype II with the SEA strains, including the two isolates sequenced in this study, forming a distinct clade denoted here as genotype IId. Genotypes III, IV and V of B. melitensis demonstrated a restricted geographical distribution, with genotype III representing the African lineage, genotype IV representing the European lineage and genotype V representing the American lineage.

Conclusion

We showed that SNPs retrieved from the B. melitensis draft full genomes were sufficient to resolve the interspecies relationships between B. melitensis strains and to discriminate between the vaccine and endemic strains. Phylogeographic reconstruction of the history of B. melitensis global spread at a finer scale by using whole-genome SNP analyses supported the origin of all B. melitensis strains from the Mediterranean region. The possible global distribution of B. melitensis following the ancient trade routes was also consistent with whole-genome SNP phylogeny. The whole genome SNP phylogenetics analysis, hence is a powerful tool for intraspecies discrimination of closely related species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1294-x) contains supplementary material, which is available to authorized users.     

Quantifying type-specific reproduction numbers for nosocomial pathogens: evidence for heightened transmission of an Asian sequence type 239 MRSA clone

An important determinant of a pathogen's success is the rate at which it is transmitted from infected to susceptible hosts. Although there are anecdotal reports that methicillin-resistant Staphylococcus aureus (MRSA) clones vary in their transmissibility in hospital settings, attempts to quantify such variation are lacking for common subtypes, as are methods for addressing this question using routinely-collected MRSA screening data in endemic settings. Here we present a method to quantify the time-varying transmissibility of different subtypes of common bacterial nosocomial pathogens using routine surveillance data. The method adapts approaches for estimating reproduction numbers based on the probabilistic reconstruction of epidemic trees, but uses relative hazards rather than serial intervals to assign probabilities to different sources for observed transmission events. The method is applied to data collected as part of a retrospective observational study of a concurrent MRSA outbreak in the United Kingdom with dominant endemic MRSA clones (ST22 and ST36) and an Asian ST239 MRSA strain (ST239-TW) in two linked adult intensive care units, and compared with an approach based on a fully parametric transmission model. The results provide support for the hypothesis that the clones responded differently to an infection control measure based on the use of topical antiseptics, which was more effective at reducing transmission of endemic clones. They also suggest that in one of the two ICUs patients colonized or infected with the ST239-TW MRSA clone had consistently higher risks of transmitting MRSA to patients free of MRSA. These findings represent some of the first quantitative evidence of enhanced transmissibility of a pandemic MRSA lineage, and highlight the potential value of tailoring hospital infection control measures to specific pathogen subtypes.     

Characteristics of the community-genotype sequence type 72 methicillin-resistant <Emphasis Type="Italic">Staphylococcus aureus</Emphasis> isolates that underlie their persistence in hospitals

Panton-Valentine leukocidin-negative methicillin-resistant Staphylococcus aureus (MRSA) clone ST72, known as a major community-associated MRSA in Korea, has emerged as an important pathogen in hospitals. To understand bacterial properties that underlie transformation of this clone into a nosocomial pathogen, we compared characteristics of the community-genotype ST72 MRSA isolates with those of ST5 and ST239 MRSA, which have been predominant nosocomial MRSA clones in Korea. Several genes associated with adhesion and virulence were absent or rarely found in ST72 isolates. Many ST72 isolates (70.1%) belonged to agr group I, but the agr group of other ST72 isolates could not be determined. As indicated by d-hemolysin production, ST72 isolates expressed fully functional agr, whereas agr dysfunction was observed in ST5 and ST239 isolates. In the biofilm formation assay, no upregulation of biofilm-forming activity of ST72 MRSA was detected. However, ST72 isolates demonstrated persistence under hypotonic and desiccating conditions (survival rates 72.3% and 33.9%, respectively), which was similar to characteristics of ST5 or ST239 isolates. ST72- MRSA isolates showed low virulence, but properties of their functional agr system could facilitate their spread in hospitals. In conclusion, tolerance to stressful environments, e.g., hypotonic and dry conditions, may also contribute to survival of the community-associated MRSA clones in healthcare facilities.     

SNP/RD typing of Mycobacterium tuberculosis Beijing strains reveals local and worldwide disseminated clonal complexes

The Beijing strain is one of the most successful genotypes of Mycobacterium tuberculosis worldwide and appears to be highly homogenous according to existing genotyping methods. To type Beijing strains reliably we developed a robust typing scheme using single nucleotide polymorphisms (SNPs) and regions of difference (RDs) derived from whole-genome sequencing data of eight Beijing strains. SNP/RD typing of 259 M. tuberculosis isolates originating from 45 countries worldwide discriminated 27 clonal complexes within the Beijing genotype family. A total of 16 Beijing clonal complexes contained more than one isolate of known origin, of which two clonal complexes were strongly associated with South African origin. The remaining 14 clonal complexes encompassed isolates from different countries. Even highly resolved clonal complexes comprised isolates from distinct geographical sites. Our results suggest that Beijing strains spread globally on multiple occasions and that the tuberculosis epidemic caused by the Beijing genotype is at least partially driven by modern migration patterns. The SNPs and RDs presented in this study will facilitate future molecular epidemiological and phylogenetic studies on Beijing strains.     

The carriage population of Staphylococcus aureus from Mali is composed of a combination of pandemic clones and the divergent Panton-Valentine leukocidin-positive genotype ST152

Staphylococcus aureus is an important human pathogen, but it appears more commonly in asymptomatic colonization of the nasopharynx than in cases of invasive disease. Evidence concerning the global population structure of S. aureus is limited by the overrepresentation in the multilocus sequence testing database of disease isolates recovered from Western Europe, the Americas, Australia, and Japan. We address this by presenting data from the S. aureus carriage population in Mali, the first detailed characterization of asymptomatic carriage from an African population. These data confirm the pandemic spread of many of the common S. aureus clones in the carriage population. We also note the high frequency (approximately 24%) of a single divergent genotype, sequence type 152 (ST152), which has not previously been recovered from nasal carriage isolates but corresponds to a sporadic Panton-Valentine leukocidin (PVL)-positive, community-acquired methicillin-resistant S. aureus clone noted mostly in Central Europe. We show that 100% of the ST152 isolates recovered from nasal carriage samples in Mali are PVL positive and discuss implications relating to the emergence and spread of this virulent genotype.     

The rise and fall of bacterial clones: Streptococcus pneumoniae

Globally spreading bacterial strains belong to clonal types that have the capacity to colonize, spread and cause disease in the community. Recent comparative genomic analyses of well-defined clinical isolates have led to the identification of bacterial properties that are required for the successful spread of bacterial clones. In this Review, we discuss the evolution of bacterial clones, the importance of recombination versus mutations for evolution of clones, common methods used to study clonal relationships among bacteria, factors that may contribute to the clonal spread of bacteria and the potential relevance of bacterial clones to clinical disease. We focus on the common pathogen Streptococcus pneumoniae, although other bacteria are also briefly discussed, such as Helicobacter pylori, Staphylococcus aureus and Mycobacterium tuberculosis.     

Evolutionary dynamics of methicillin-resistant Staphylococcus aureus within a healthcare system

BackgroundIn the past decade, several countries have seen gradual replacement of endemic multi-resistant healthcare-associated methicillin-resistant Staphylococcus aureus (MRSA) with clones that are more susceptible to antibiotic treatment. One example is Singapore, where MRSA ST239, the dominant clone since molecular profiling of MRSA began in the mid-1980s, has been replaced by ST22 isolates belonging to EMRSA-15, a recently emerged pandemic lineage originating from Europe.ResultsWe investigated the population structure of MRSA in Singaporean hospitals spanning three decades, using whole genome sequencing. Applying Bayesian phylogenetic methods we report that prior to the introduction of ST22, the ST239 MRSA population in Singapore originated from multiple introductions from the surrounding region; it was frequently transferred within the healthcare system resulting in a heterogeneous hospital population. Following the introduction of ST22 around the beginning of the millennium, this clone spread rapidly through Singaporean hospitals, supplanting the endemic ST239 population. Coalescent analysis revealed that although the genetic diversity of ST239 initially decreased as ST22 became more dominant, from 2007 onwards the genetic diversity of ST239 began to increase once more, which was not associated with the emergence of a sub-clone of ST239. Comparative genomic analysis of the accessory genome of the extant ST239 population identified that the Arginine Catabolic Mobile Element arose multiple times, thereby introducing genes associated with enhanced skin colonization into this population.ConclusionsOur results clearly demonstrate that, alongside clinical practice and antibiotic usage, competition between clones also has an important role in driving the evolution of nosocomial pathogen populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0643-z) contains supplementary material, which is available to authorized users.     

Yersinia pestis strains isolated in natural plague foci of Caucasus and Transcaucasia in the context of the global evolution of species

BackgroundPlague is a highly dangerous vector-borne infectious disease that has left a significant mark on history of humankind. There are 13 natural plague foci in the Caucasus, located on the territory of the Russian Federation, Azerbaijan, Armenia and Georgia. We performed whole-genome sequencing of Y. pestis strains, isolated in the natural foci of the Caucasus and Transcaucasia. Using the data of whole-genome SNP analysis and Bayesian phylogeny methods, we carried out an evolutionary-phylogeographic analysis of modern population of the plague pathogen in order to determine the phylogenetic relationships of Y. pestis strains from the Caucasus with the strains from other countries.ResultsWe used 345 Y. pestis genomes to construct a global evolutionary phylogenetic reconstruction of species based on whole-genome SNP analysis. The genomes of 16 isolates were sequenced in this study, the remaining 329 genomes were obtained from the GenBank database. Analysis of the core genome revealed 3315 SNPs that allow differentiation of strains. The evolutionary phylogeographic analysis showed that the studied Y. pestis strains belong to the genetic lineages 0.PE2, 2.MED0, and 2.MED1. It was shown that the Y. pestis strains isolated on the territory of the East Caucasian high-mountain, the Transcaucasian high-mountain and the Priaraksinsky low-mountain plague foci belong to the most ancient of all existing genetic lineages - 0.PE2.ConclusionsOn the basis of the whole-genome SNP analysis of 345 Y. pestis strains, we describe the modern population structure of the plague pathogen and specify the place of the strains isolated in the natural foci of the Caucasus and Transcaucasia in the structure of the global population of Y. pestis. As a result of the retrospective evolutionary-phylogeographic analysis of the current population of the pathogen, we determined the probable time frame of the divergence of the genetic lineages of Y. pestis, as well as suggested the possible paths of the historical spread of the plague pathogen.     

Complex evolutionary history of a Neotropical lowland forest bird (Lepidothrix coronata) and its implications for historical hypotheses of the origin of Neotropical avian diversity

Here we apply a combination of phylogeographic and historical demographic analyses to the study of mtDNA sequence variation within the Blue-crowned Manakin (Lepidothrix coronata), a widespread Neotropical bird. A high degree of phylogeographic structure allowed us to demonstrate that several vicariant events, including Andean uplift, the formation of riverine barriers, and climatically induced vegetational shifts, as well as a non-vicariant process, range expansion, have all acted, at varying spatial and temporal scales, to influence genetic structure within L. coronata, suggesting that current historical hypotheses of the origin of Neotropical avian diversity that focus on single vicariant mechanisms may be overly simplistic. Our data also support an origin (>2 mybp) that is substantially older than the late Pleistocene for the genetic structure within this species and indicate that phylogeographic patterns within the species are not concordant with plumage-based subspecific taxonomy. These data add to a growing body of evidence suggesting that the origin of several Neotropical avian species may have occurred in the mid-Pliocene, thus, geological arguments surrounding putative Pleistocene vicariant events, while interesting in their own right, may have little relevance to Neotropical avian diversification at the species level.     

Comparative genomics of methicillin-resistant Staphylococcus aureus ST239: distinct geographical variants in Beijing and Hong Kong

Background

The ST239 lineage is a globally disseminated, multiply drug-resistant hospital-associated methicillin-resistant Staphylococcus aureus (HA-MRSA). We performed whole-genome sequencing of representative HA-MRSA isolates of the ST239 lineage from bacteremic patients in hospitals in Hong Kong (HK) and Beijing (BJ) and compared them with three published complete genomes of ST239, namely T0131, TW20 and JKD6008. Orthologous gene group (OGG) analyses of the Hong Kong and Beijing cluster strains were also undertaken.

Results

Homology analysis, based on highest-percentage nucleotide identity, indicated that HK isolates were closely related to TW20, whereas BJ isolates were more closely related to T0131 from Tianjin. Phylogenetic analysis, incorporating a total of 30 isolates from different continents, revealed that strains from HK clustered with TW20 into the ‘Asian clade’, whereas BJ isolates and T0131 clustered closely with strains of the ‘Turkish clade’ from Eastern Europe. HK isolates contained the typical φSPβ-like prophage with the SasX gene similar to TW20. In contrast, BJ isolates contained a unique 15 kb PT1028-like prophage but lacked φSPβ-like and φSA1 prophages. Besides distinct mobile genetic elements (MGE) in the two clusters, OGG analyses and whole-genome alignment of these clusters highlighted differences in genes located in the core genome, including the identification of single nucleotide deletions in several genes, resulting in frameshift mutations and the subsequent predicted truncation of encoded proteins involved in metabolism and antimicrobial resistance.

Conclusions

Comparative genomics, based on de novo assembly and deep sequencing of HK and BJ strains, revealed different origins of the ST239 lineage in northern and southern China and identified differences between the two clades at single nucleotide polymorphism (SNP), core gene and MGE levels. The results suggest that ST239 strains isolated in Hong Kong since the 1990s belong to the Asian clade, present mainly in southern Asia, whereas those that emerged in northern China were of a distinct origin, reflecting the complexity of dissemination and the dynamic evolution of this ST239 lineage.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-529) contains supplementary material, which is available to authorized users.