Previous Article | Next Article ![]()
Journal of Bacteriology, August 2005, p. 5709-5718, Vol. 187, No. 16
0021-9193/05/$08.00+0 doi:10.1128/JB.187.16.5709-5718.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
George M. Weinstock,3,4 and
Barbara E. Murray1,2,3*
Division of Infectious Diseases, Department of Internal Medicine,1 Center for the Study of Emerging and Re-emerging Pathogens,2 Department of Microbiology and Molecular Genetics, University of Texas Medical School,3 Human Genome Sequencing Center, Baylor College of Medicine, Houston, Texas 770304
Received 21 March 2005/ Accepted 25 May 2005
|
|
|---|
|
|
|---|
Molecular typing has shown that discrete lineages of pathogenic bacteria can arise periodically and then spread locally or globally in the presence of strong selective pressure (26). Thus, the identification of E. faecalis clones that are successful in achieving prolonged, widespread outbreaks and the unraveling of their genetic background may shed light on the question of how this opportunist adapts to clinical settings and behaves as a pathogen, causing a range of infections such as intraabdominal, genitourinary, endovascular, or meningeal infections, among others. Among typing methods for examining relatedness of bacterial genetic backgrounds, multilocus sequence typing (MLST), which is objective and less prone to human error, has gained recognition as one of the best approaches and has been used to identify pathogenic lineages of several species, including Neisseria meningitidis, Streptococcus pneumoniae, S. aureus, and Enterococcus faecium (4, 7, 10, 26), among others. Recently, we derived a four-gene MLST system using one housekeeping gene (pyrC) and three antigen-encoding genes (ace, efaA, and salA), chosen for their likely greater variation, which successfully differentiated E. faecalis at the subspecies level (23). When we subsequently used this MLST scheme to examine additional selected isolates from our 30-year collection, the results suggested an unexpected relationship among clinically important isolates disseminated in several states of the United States. We decided to further investigate these isolates using additional housekeeping genes plus the antigen-encoding genes; this combined use of different types of genes has the potential advantage of revealing both the long-term evolutionary history of the chromosome and a short-term differentiation resulting from the more variable antigen-encoding genes. After identifying a circulating E. faecalis lineage that had acquired resistances to the primary and secondary drugs of choice, ampicillin and vancomycin, we further explored the virulence gene profiles and predicted the evolution of this clonal complex based on acquired resistance genes and variations observed in the previously reported pathogenicity island (PAI) (34).
(A part of this work was presented at the 43rd Interscience Conference on Antimicrobial Agents and Chemotherapy, 2003, Chicago, Ill.)
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Enterococcus faecalis strains by source of isolation, PFGE type, and nine-gene MLST
|
|
View this table: [in a new window] |
TABLE 2. Details of amplicons and primers used for nine-gene MLST
|
1 nucleotide change was classified as a distinct allele (no weight was given to the degree of sequence divergence between alleles) and each isolate was defined by its allelic profile (a series of numbers corresponding to the alleles at the nine loci). In keeping with other studies, isolates with the same nine allelic profiles were assigned the same sequence type (ST), and isolates that shared alleles at
7 loci (single or double locus variants [SLVs or DLVs]) were called a clone. For linked STs differing by SLVs, the term clonal complex was used, with the implication that these are descendants of a common ancestor. Clonality was assessed using BURST, a clustering algorithm designed for use with MLST data sets of bacterial pathogens (9). PFGE and hybridizations. PFGE was performed with some modifications of a previously described method (19). Agarose plugs containing genomic DNA were digested with SmaI (Invitrogen), and electrophoresis was carried out using clamped homogeneous electric field (CHEF-DRII; Bio-Rad Laboratories, Richmond, Calif.), with ramped pulse times beginning with 5 s and ending with 45 s, at 200 V for 26 h. The PFGE patterns were interpreted using the criteria suggested by Tenover et al. (37), with closely and possibly related patterns being designated as belonging to a single clone. PFGE pattern names that were presented in earlier publications are used here. Southern and colony lysate hybridizations were performed under high-stringency conditions (36) with probes labeled using the RadPrime DNA labeling system (Invitrogen). Probe details and primers used for amplification are listed in Tables 3 and 4.
|
View this table: [in a new window] |
TABLE 3. Details of primers used for generating virulence-associated gene probes
|
|
View this table: [in a new window] |
TABLE 4. Details of primers used for generating PAI-associated gene probes
|
|
|
|---|
![]() View larger version (20K): [in a new window] |
FIG. 1. Chromosomal locations of the nine MLST loci in the E. faecalis V583 genome (28). The arrowheads inside the circle represent the open reading frame orientations of each locus (not drawn to scale). Antigenic gene arrowheads are shaded in black, and housekeeping gene arrowheads are shaded in gray. The arrowhead of dnaA (coding for chromosomal replication factor) positioned at nucleotide 1 is not shaded. The distance between any two loci ranged from 103 kb to 1,247 kb. The putative PAI region is boxed.
|
![]() View larger version (30K): [in a new window] |
FIG. 2. Allelic variation of the MLST loci sequenced for analysis. Allelic variations of ace, efaA, and salA were as previously published (23) and hence are not shown here. For these three genes, alleles A to I of the previous study (23) were redefined as alleles 1 to 9 for this study. The nucleotides present in each of the variable sites of allele 1 (E. faecalis OG1RF) are shown. Only those sites that differ are shown for the other alleles. The position of each variable site within the sequenced fragment is shown by the number above the nucleotide, read vertically.
|
![]() View larger version (30K): [in a new window] |
FIG. 3. Analysis of an unusual E. faecalis clonal cluster by MLST and PFGE. (A) Clustering of four clonally related STs using BURST analysis. The central circle denotes the predominant type among isolates tested using a nine-gene MLST scheme, and each surrounding circle indicates one allele difference. A dashed straight line denotes a double locus difference. (B) PFGE fingerprints. The PFGE phylogenetic tree was based on the unweighted pair group method. Tolerance of up to 5% shift in the band position was used. Isolates are generally referred to by their previously published designations. Year and place of isolation and ST type are shown. Isolates from outbreaks are marked with an asterisk. V583 is the first vancomycin-resistant enterococcus isolated in the United States, and HH-22 is the first-known Bla+ isolate of E. faecalis. E228 and E366, isolated a year apart from the same hospital and differing by three bands, are represented by E228.
|
2% of isolates, while Bla producers are even more rare; of note, however, a vancomycin-resistant E. faecalis was the donor of the vanA genes in at least one of the recent descriptions of vancomycin-resistant methicillin-resistant S. aureus (39). The BVE clone does not appear to be predominant among E. faecalis isolates in general and was found in only 2 of over 70 other independent isolates examined, analyzed by two or more methods of PFGE, multilocus enzyme electrophoresis (38), and MLST. Analysis of PFGE fingerprints of this BVE clonal complex also showed related PFGE patterns, differing by only a few fragments, which would categorize them as closely or possibly related (Fig. 3B) using criteria for analysis of potential nosocomial outbreaks (37). Hybridization of PFGE Southern blots of chromosomal digestion fragments with probes for gdh, aroE, yqiL, gki, and xpt showed hybridization to, for the majority of fragments, the same-sized band in isolates of the BVE clonal complex, with sizes as expected from the V583 genome. One of the exceptions (a difference in size of the xpt hybridizing band of V583) is known to be due to insertion of the vancomycin resistance (vanB) element in this strain (28). Although PFGE and hybridizations confirmed the relatedness of all the MLST-defined BVE clonal complex isolates, the degrees of difference inferred by the two techniques did not always strictly overlap; this is not surprising, since PFGE and MLST have different molecular bases.
Among two other new relationships identified, one clone named HV1 (Houston Vanr #1) represents the unrecognized persistence of a clone identified earlier (1, 11, 23) in a single hospital. Another unexpected observation was that two other Bla+ isolates (ST-9 and its DLV, ST-13, differing only in antigenic genes), representative of small nosocomial outbreaks in Connecticut (27) and in Argentina (17), respectively, belong to the same clone, named ACB (Argentina-Connecticut-Bla+); this clone is unrelated to the BVE clonal complex, with variations at eight of nine loci. Among other Bla+ E. faecalis isolates known to date, one isolate from Lebanon (20, 23, 38) and an outbreak strain (20, 30) from Boston (not included in this study) are not related to the BVE or ACB groups (20). This implies that the staphylococcal blaZ gene has spread only a few times into E. faecalis and suggests that certain backgrounds may be particularly receptive to blaZ acquisition.
To further assess the genetic content of the isolates in this study, we generated a profile based upon hybridization to 14 chromosomally encoded potential virulence genes (Fig. 4A and Table 3). Many, including six recently described immunoglobulin (Ig)-like fold-containing putative microbial surface adhesins (35), were present in all 21 isolates. However, four genes were variably present and the differences corresponded to different STs; the variably present genes (Fig. 4A) were fsrB (encoding part of the Fsr two-component system of E. faecalis, which regulates the virulence genes gelE and sprE [29]); ef1824 (encoding a predicted adhesin [35] with a characteristic Ig-like fold [3]); and hylA and hylB (ef3023 and ef0818, respectively), each encoding a putative hyaluronidase, an enzyme implicated in pathogenesis in other organisms. Although gelE was found in all 21 isolates, the lack of gelatinase production by isolates in ST-4, ST-8, and ST-11 is related to a previously described 23.9-kb deletion (22) which includes fsrB. Based on the observation that ace B repeats (24) are separated by recer sequences, which may promote recombination and thus variable repeat numbers (2), we also tested the number of B repeats of the E. faecalis-specific collagen adhesin Ace (Fig. 4A). While neither the gene profile nor ace B repeats were alone sufficiently discriminatory, the combination of these two profiles successfully distinguished the lineages from one another, reflective of the MLST and PFGE types (Fig. 4A). Notably, the BVE clonal complex, including the two additional members recognized by multilocus enzyme electrophoresis or PFGE (both HL-Gmr) (38) but not tested by MLST, contained all of the potential virulence genes.
![]() View larger version (30K): [in a new window] |
FIG. 4. Schematic presentation of hybridization profiles for potential virulence-related genes and PAI genes. (A) Virulence-related gene (non-PAI) profile and ace B repeat number profile of 21 MLST-analyzed isolates. The 10 genes (of ef0080, ef0089, ef0786, ef1091, ef1092, ef1093, ef1269, ef1818, ef2224, and ef3191) present in all 21 isolates are not shown. The STs with the 23.9-kb deletion involving fsrB are marked with a superscript "a" in the fsrB data. (B) Determination of the PAI insertion site by PCR using primers within and outside the PAI of V583. The double arrow denotes an expected-size PCR fragment with primers PAIout plus ef0481forward or ef0482forward. STs yielding an 1.7-kb-larger PCR product or an 1-kb-larger PCR product are marked with superscript "b" and "c," respectively. (C) Hybridization results with PAI-associated intragenic probes, representing 18 genes dispersed over the entire PAI region. The superscript numbers on ST and the +/ symbols denote the number of isolates of that type. The integrated plasmid region (ef0506 to ef0485) (28) is boxed. (D) Distribution of three PAI-associated genes among 341 clinical isolates, 58 nosocomial stool isolates, 33 community-derived stool isolates, and 22 animal isolates.
|
150-kb region as an E. faecalis PAI (ef0479 to ef0628 of V583) and indicated that there were only subtle differences in this region in two of the E. faecalis isolates described above, V583 (28, 32) and MMH594 (8). However, our recognition that these two isolates are actually members of the same ST suggests that the highly similar nature of the putative PAI of these isolates is a function of their close evolutionary relationship. To further assess whether a similar PAI was present in the other BVE clonal complex isolates and in unrelated strains, we tested the 21 isolates described above, which belong to 14 STs within nine different clonal lineages, with probes representing different genes dispersed over the entire PAI region (Fig. 4C). Colony lysate hybridization results with the 18 individually labeled intragenic probes (Table 4) showed that all 18 genes (which are all present in MMH594, in which the PAI was first identified [34]), are also present in two unrelated STs (ST-4 and ST-11, both represented by HL-Gmr isolates from Thailand isolated in 1980) (Fig. 4C). The remaining isolates representing 10 STs contain an incomplete PAI with deletions in different regions, except the ST-1 isolate (OGIRF [21]), which contains none of the 18 genes (Fig. 4C). Hybridization of BVE clonal complex members to xylA and gls24-like gene probes confirmed the presence of these PAI genes within the same-sized PFGE fragments of the BVE clonal complex. However, within this complex, there were isolate-specific PAI deletions localized to three regions, one including ef0530 and ef0534, a second in the middle (ef0571), and the third including ef0604 and ef0609 (Fig. 4C). This finding of PAI variability is not unexpected, considering the frequent occurrence of IS-like elements and integrase and recombinase genes in the PAI region of MMH594 and V583 (28, 34), including the previously described 17-kb deletion in PAI of V583 versus MMH594, both ST-14 isolates (34). Although there are many differences within the PAI of isolates of different lineages, possibly due to deletions, the finding of PAI-associated genes in eight of nine lineages containing isolates from around the world corroborates the earlier proposal that PAI is disseminated among different strains.
To investigate a possible common insertion site of apparently transferable PAI (25) in isolates from different STs, PCR was performed using one primer located outside the PAI (PAIout) and the second located within the PAI region (ef0481forward or ef0482forward [Table 4]); two isolates lacking both ef0481 and ef0482 were not tested. Products of 3.8 kb (with ef0481forward primer) and 4.2 kb (with ef0482forward primer) (sizes were as anticipated from the V583 genome sequence) were obtained with DNA from 16 of 19 isolates (Fig. 4B), suggesting the same PAI insertion site in different E. faecalis clones. Using the same sets of primers,
1.7-kb and
1-kb larger PCR products were obtained with DNA from the HV1 clone and from the ST-8 isolate, respectively, indicating further small insertions in this region. Thus, these results predict that at least seven of the nine lineages have a common PAI insertion site.
We also tested an additional 454 geographically and temporally diverse isolates for the presence of three selected PAI-associated genes, one close to each end and one in the middle. Hybridization results indicate that all three PAI-associated genes (esp, xylA, and gls24-like) are distributed worldwide and are enriched in infection-derived isolates (P < 0.0025 compared to community-derived isolates from human stools or animals), extending a previous report using 80 isolates of unknown clonal relatedness (34). Among the 341 clinical isolates, 17.6% were found to contain all three genes, 41.4% contained combinations of esp plus xylA or xylA plus gls24-like, and 27% of isolates lacked all three genes. The variability of the PAI region is consistent with the many deletions identified above for the well-characterized lineages (Fig. 4D). The frequent finding of two or more PAI genes, together with the results for a common insertion site, suggests that these PAI genes were acquired as a unit with subsequent deletions. The less frequent occurrence of PAI-associated genes in nosocomial stool isolates compared to nosocomial clinical isolates (49.9% versus 31% for esp; 62.5% versus 39.7% for xylA) is likely because fecal isolates of hospitalized patients include a mixture of both nosocomially derived and community-derived organisms. The uncommon occurrence of even one of these three PAI-associated genes in non-human-derived isolates (9.1%) plus the high frequency of occurrence of at least one of the three genes in clinical isolates (73%) support the hypothesis that this genomic region may be helpful during some stage of human infection. The results also suggested that some deletions may be favored, or may be a clonal marker, in specific clinical settings, as exemplified by the very frequent absence of the gls24-like region in endocarditis isolates.
In a further analysis of individual isolates of the unusual BVE clonal complex, we used the PAI region variability, together with locus variations and the presence of antibiotic resistances, to predict the evolutionary pathway of this distinctive lineage (Fig. 5). The most complete PAI region (like that present in two Thailand strains) was found in MMH594, which lacks blaZ and vanB, and so this isolate or some predecessor was positioned in an ancestral position. A large PAI region (including cylM, the 17-kb region previously noted as deleted from V583 [34], araC [ef0530], and ef0534) was missing from all Bla+ (ST-6 and ST-7) isolates of this clonal complex (Fig. 4C and 5,
B), and long-range PCR with primers outside this deletion (ef0521reverse and ef0539forward [Table 4]) yielded the same-sized PCR fragments (
10 kb) with all Bla+ isolates of the BVE clonal complex, confirming the same deletion in these isolates. An additional region (ef0604 to ef0628) (Fig. 5,
E) of the PAI was absent in a single ST-7 isolate, HH-22, suggesting that this isolate, although it was the first Bla+ enterococcus identified, was not the ancestor of the later, more widespread ST-6 isolates. The endocarditis isolate of the BVE clonal complex showed two independent PAI deletions, one in the middle (ef0571) and the other including ef0604 and ef0609 (Fig. 4C and 5,
C and
D). In our scheme, the acquisition of HL-Gmr is inferred to have occurred after PAI acquisition, although the reverse could also be true.
![]() View larger version (34K): [in a new window] |
FIG. 5. Predicted path of evolution of the BVE clonal complex based on MLST, deletions in the PAI region, year of isolation, and the presence of uncommon resistance genes, blaZ and vanB. Isolates of the BVE clonal complex contained all 14 tested potential virulence genes (Fig. 4). Isolates of each ST within the clonal complex are grouped by oval shading, and SLVs are denoted by overlapping shaded ovals. Arrows with continuous and dashed lines represent the predicted and alternative evolutionary paths. Putative ancestral isolates are boxed with a dashed line. Five partial deletions ( A to E) of the PAI of E. faecalis strain MMH594, including that previously described for V583 (34), are shown. VanB, vancomycin resistance encoded by the vanB gene.
|
This work is supported by NIH grant R37 AI47923 from the Division of Microbiology and Infectious Diseases to B. E. Murray.
Present address: Department of Infectious Diseases, Chongqing University of Medical Sciences, Chongqing, China 400016. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»