Recent studies have pointed to the existence of two subpopulations of Enterococcus faecium, one containing primarily commensal/community-associated (CA) strains and one that contains most clinical or hospital-associated (HA) strains, including those classified by multi-locus sequence typing (MLST) as belonging to the CC17 group. The HA subpopulation more frequently has IS16, pathogenicity island(s), and plasmids or genes associated with antibiotic resistance, colonization, and/or virulence. Supporting the two clades concept, we previously found a 3–10% difference between four genes from HA-clade strains vs. CA-clade strains, including 5% difference between pbp5-R of ampicillin-resistant, HA strains and pbp5-S of ampicillin-sensitive, CA strains. To further investigate the core genome of these subpopulations, we studied 100 genes from 21 E. faecium genome sequences; our analyses of concatenated sequences, SNPs, and individual genes all identified two distinct groups. With the concatenated sequence, HA-clade strains differed by 0–1% from one another while CA clade strains differed from each other by 0–1.1%, with 3.5–4.2% difference between the two clades. While many strains had a few genes that grouped in one clade with most of their genes in the other clade, one strain had 28% of its genes in the CA clade and 72% in the HA clade, consistent with the predicted role of recombination in the evolution of E. faecium. Using estimates for Escherichia coli, molecular clock calculations using sSNP analysis indicate that these two clades may have diverged ≥1 million years ago or, using the higher mutation rate for Bacillus anthracis, ~300,000 years ago. These data confirm the existence of two clades of E. faecium and show that the differences between the HA and CA clades occur at the core genomic level and long preceded the modern antibiotic era.
References
[1]
Galloway-Pena JR, Nallapareddy SR, Arias CA, Eliopoulos GM, Murray BE (2009) Analysis of clonality and antibiotic resistance among early clinical isolates of Enterococcus faecium in the United States. J Infect Dis 200: 1566–1573.
[2]
Hidron AI, Edwards JR, Patel J, Horan TC, Sievert DM, et al. (2008) NHSN annual update: antimicrobial-resistant pathogens associated with healthcare-associated infections: annual summary of data reported to the National Healthcare Safety Network at the Centers for Disease Control and Prevention, 2006–2007. Infect Control Hosp Epidemiol 29: 996–1011.
[3]
Hendrickx AP, van Wamel WJ, Posthuma G, Bonten MJ, Willems RJ (2007) Five genes encoding surface-exposed LPXTG proteins are enriched in hospital-adapted Enterococcus faecium clonal complex 17 isolates. J Bacteriol 189: 8321–8332.
[4]
Klare I, Konstabel C, Badstubner D, Werner G, Witte W (2003) Occurrence and spread of antibiotic resistances in Enterococcus faecium. Int J Food Microbiol 88: 269–290.
[5]
Leavis HL, Bonten MJ, Willems RJ (2006) Identification of high-risk enterococcal clonal complexes: global dispersion and antibiotic resistance. Curr Opin Microbiol 9: 454–460.
[6]
Leavis HL, Willems RJ, van Wamel WJ, Schuren FH, Caspers MP, et al. (2007) Insertion sequence-driven diversification creates a globally dispersed emerging multiresistant subspecies of E. faecium. PLoS Pathog 3: e7.
[7]
Nallapareddy SR, Weinstock GM, Murray BE (2003) Clinical isolates of Enterococcus faecium exhibit strain-specific collagen binding mediated by Acm, a new member of the MSCRAMM family. Mol Microbiol 47: 1733–1747.
[8]
Panesso D, Reyes J, Rincon S, Diaz L, Galloway-Pena J, et al. (2010) Molecular epidemiology of vancomycin-resistant Enterococcus faecium: a prospective, multicenter study in South American hospitals. J Clin Microbiol 48: 1562–1569.
[9]
Rice LB, Carias L, Rudin S, Vael C, Goossens H, et al. (2003) A potential virulence gene, hylEfm, predominates in Enterococcus faecium of clinical origin. J Infect Dis 187: 508–512.
[10]
Sillanpaa J, Prakash VP, Nallapareddy SR, Murray BE (2009) Distribution of genes encoding MSCRAMMs and Pili in clinical and natural populations of Enterococcus faecium. J Clin Microbiol 47: 896–901.
[11]
Werner G, Fleige C, Geringer U, van Schaik W, Klare I, et al. (2011) IS element IS16 as a molecular screening tool to identify hospital-associated strains of Enterococcus faecium. BMC Infect Dis 11: 80.
[12]
Top J, Willems R, Bonten M (2008) Emergence of CC17 Enterococcus faecium: from commensal to hospital-adapted pathogen. FEMS Immunol Med Microbiol 52: 297–308.
[13]
Willems RJ, Hanage WP, Bessen DE, Feil EJ (2011) Population biology of Gram-positive pathogens: high-risk clones for dissemination of antibiotic resistance. FEMS Microbiol Rev 35: 872–900.
[14]
Didelot X, Falush D (2007) Inference of bacterial microevolution using multilocus sequence data. Genetics 175: 1251–1266.
[15]
Tang J, Hanage WP, Fraser C, Corander J (2009) Identifying currents in the gene pool for bacterial populations using an integrative approach. PLoS Comput Biol 5: e1000455.
[16]
Turner KM, Hanage WP, Fraser C, Connor TR, Spratt BG (2007) Assessing the reliability of eBURST using simulated populations with known ancestry. BMC Microbiol 7: 30.
[17]
Willems RJ, van Schaik W (2009) Transition of Enterococcus faecium from commensal organism to nosocomial pathogen. Future Microbiol 4: 1125–1135.
[18]
van Schaik W, Top J, Riley DR, Boekhorst J, Vrijenhoek JE, et al. (2010) Pyrosequencing-based comparative genome analysis of the nosocomial pathogen Enterococcus faecium and identification of a large transferable pathogenicity island. BMC Genomics 11: 239.
[19]
Galloway-Pena JR, Rice LB, Murray BE (2011) Analysis of PBP5 of early U.S. isolates of Enterococcus faecium: sequence variation alone does not explain increasing ampicillin resistance over time. Antimicrob Agents Chemother 55: 3272–3277.
[20]
Foster JT, Beckstrom-Sternberg SM, Pearson T, Beckstrom-Sternberg JS, Chain PS, et al. (2009) Whole-genome-based phylogeny and divergence of the genus Brucella. J Bacteriol 191: 2864–2870.
[21]
Lenski RE, Winkworth CL, Riley MA (2003) Rates of DNA sequence evolution in experimental populations of Escherichia coli during 20,000 generations. J Mol Evol 56: 498–508.
[22]
Pearson T, Giffard P, Beckstrom-Sternberg S, Auerbach R, Hornstra H, et al. (2009) Phylogeographic reconstruction of a bacterial species with high levels of lateral gene transfer. BMC Biol 7: 78.
[23]
Van Ert MN, Easterday WR, Huynh LY, Okinaka RT, Hugh-Jones ME, et al. (2007) Global genetic population structure of Bacillus anthracis. PLoS One 2: e461.
[24]
Achtman M, Morelli G, Zhu P, Wirth T, Diehl I, et al. (2004) Microevolution and history of the plague bacillus, Yersinia pestis. Proc Natl Acad Sci U S A 101: 17837–17842.
[25]
Ochman H, Elwyn S, Moran NA (1999) Calibrating bacterial evolution. Proc Natl Acad Sci U S A 96: 12638–12643.
[26]
Ochman H, Wilson AC (1987) Evolution in bacteria: evidence for a universal substitution rate in cellular genomes. J Mol Evol 26: 74–86.
[27]
Kuo CH, Ochman H (2009) Inferring clocks when lacking rocks: the variable rates of molecular evolution in bacteria. Biol Direct 4: 35.
[28]
Zuckerkandl E, Pauling L (1965) Evolutionary divergence and convergence in proteins; Vogel VBaHJ, editor. New York: Academic Press. pp. 97–166.
[29]
Kumar S, Nei M, Dudley J, Tamura K (2008) MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform 9: 299–306.
[30]
Tamura K, Dudley J, Nei M, Kumar S (2007) MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 24: 1596–1599.
Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, et al. (2003) Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res 31: 3497–3500.
[33]
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, et al. (2007) Clustal W and Clustal X version 2.0. Bioinformatics 23: 2947–2948.
[34]
Thompson JD, Higgins DG, Gibson TJ (1994) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22: 4673–4680.
