ABSTRACT
Environmental Enterococcus spp. were compared by BOX-PCR genotyping and 16S rRNA gene sequencing to clarify the predictive relationship of BOX-PCR fingerprints to species designation. BOX-PCR and 16S rRNA gene relationships agreed for 77% of strains. BOX-PCR provided superior intraspecies discrimination but incorrectly identified some strains to the species level and divided some species into multiple groups.
TEXT
Enterococcus species are aerotolerant fermentative, Gram-positive, catalase-negative cocci that are commonly found in the gastrointestinal (GI) tract of mammals and birds and are readily isolated from soil, surface waters, sediments, and vegetation associated with surface waters (1, 5, 12). Enterococci are used as regulatory tools to assess water quality in fresh and saline waters (36). Some Enterococcus species possess virulence factors and antibiotic resistance genes and are capable of causing disease, particularly in individuals with compromised health (29, 32, 34). Vancomycin-resistant enterococci (VRE) are important nosocomial pathogens that have been isolated from approximately 30% of the patients in intensive care units in U.S. hospitals (25, 30).
Accurate identification and classification of the different species belonging to the genus Enterococcus are important for both environmental and clinical studies (2, 3, 6, 10, 39). Phenotypic methods used to identify Enterococcus species are not very discriminatory or accurate due to the phenotypic similarity of certain species, such as E. gallinarum and E. casseliflavus, E. cecorum and E. columbae, and E. hirae and E. durans (9). Therefore, DNA-based analyses, such as genotyping, genetic sequencing, and targeting particular genes with specific probes and primers, are also employed for the identification and classification of enterococci isolated from environmental and clinical samples (8, 15, 17, 21, 31).
Sequencing of the 16S rRNA gene is considered the “gold standard” for microbial identification (7, 14, 38); however, genotyping methods, such as BOX-PCR, are high throughput and cost-effective compared to sequencing analysis for processing a large number of isolates. Furthermore, such genotyping methods are often more discriminatory within species than 16S rRNA gene sequencing. Some environmental studies have used BOX-PCR typing to determine genotypic relationships of enterococci (6, 16). These studies include dendrograms that display population similarities based solely on BOX-PCR genotyping. We are not aware of any studies demonstrating that the BOX-PCR patterns of various strains of a particular Enterococcus species are more similar than strains from different species, in spite of the fact that the BOX-PCR genotypes are assumed to reflect phylogenetic relatedness.
This study cross-validated two commonly used molecular methods for determining phylogenetic relationships of Enterococcus species isolates: BOX-PCR typing and 16S rRNA gene sequencing. Water, sediment, and vegetation samples were collected during two separate sampling events from two freshwater bodies (Lake Carroll and Hillsborough River) and an estuarine site (Ben T Davis beach) in the Tampa Bay area in Florida. Each sample was given a unique identifier (ID) consisting of the site location (LC for Lake Carroll, HR for Hillsborough River, and BD for Ben T Davis beach), isolate number, type of sample (W for water, S for sediment, and V for vegetation), and sampling event (07 for the first event and 08 for the second event). Sample processing and DNA extraction details were described previously (3).
DNA extracted from 61 individual, randomly selected isolates was amplified using the bacterial universal primers 8f and 1492r for the 16S rRNA gene (20). The PCR master mix and conditions are described elsewhere (3). Sequences were assembled using Sequencher 4.8 (Gene Codes Corp., Ann Arbor, MI) and analyzed using BLAST (http://www.ncbi.nlm.nih.gov/BLAST/) for identification. Phylogenetic dendrograms were created using the neighbor-joining method, and the bootstrap test was carried out with 500 replications (MEGA 4.0 [33a]).
Enterococcus isolates (n = 61) were fingerprinted using the horizontal, fluorophore-enhanced, repetitive extragenic palindromic (rep)-PCR (HFERP) technique (18). Working primer stock was prepared by mixing 0.09 μg of unlabeled BOX A2R primer (21) (0.68 μg μl−1) per μl and 0.03 μg of 6-FAM (fluorescein)-labeled BOX A2R primer (0.74 μg μl−1) per μl. The PCR master mix and conditions are described elsewhere (3). Enterococcus faecium C68 was used as the control strain in PCR and was run on each individual gel to determine intergel variability. A ladder plus nonmigrating loading dye mixture was prepared by mixing 5 μl of Genescan-2500 ROX internal lane standard (Applied Biosystems, Foster City, CA) and 20 μl dye (150 mg Ficoll 400 per ml and 25 mg blue dextran per ml). A total of 12.5 μl of individual PCR products was mixed with 3.3 μl of the ROX-dye mixture, loaded in a 1.5% agarose gel, and electrophoresed at 90 V for 4 h. Gel images were scanned using a Typhoon 8600 variable-mode imager (GE Healthcare), imported into Bionumerics (Applied Maths, Belgium), and analyzed using the Pearson similarity coefficient by constructing dendrograms using the unweighted-pair group method with arithmetic mean (UPGMA) (optimization 1.0%, tolerance 0.5%). A control strain (E. faecium C68) was included on all gels, and the percent similarity (89%) for the control strain patterns was considered the cutoff value during similarity analyses. BOX-PCR patterns were also compared visually to confirm results.
Known Enterococcus species, including Enterococcus faecalis ATCC 19433, Enterococcus faecalis ATCC 29212, Enterococcus faecalis ATCC 49383, Enterococcus faecalis ATCC 700802, Enterococcus faecium C68, Enterococcus faecium ATCC 49224, and Enterococcus casseliflavus ATCC 700327, were also sequenced and typed for comparison purposes.
Agreement between the typing and sequencing results was calculated by setting the DNA sequencing results as the “standard” and determining the percentage of BOX-PCR patterns that clustered with the correct clade based on their respective 16S rRNA genes. Since BOX-PCR patterns of the control strain E. faecium C68 compared from multiple gels were 89% similar, patterns with similarity values greater than 89% were considered indistinguishable here. As we hypothesized, BOX-PCR dendrograms (see Fig. 2) were in good agreement (77%) with the 16S rRNA gene phylogenetic tree (Fig. 1). Sequencing of the 16S rRNA gene was incapable of differentiating among the known strains of E. faecalis and also among the known strains of E. faecium (Fig. 1); however, it readily discriminated between the two species. The BOX-PCR patterns of three known strains of E. faecalis (ATCC 19433, ATCC 29212, and ATCC 700802) were 95% similar and identical when examined visually (Fig. 2). Minor differences (90% similarity) were observed in the BOX-PCR patterns of the two known strains of E. faecium; however, this similarity level was greater than the cutoff established by the E. faecium C68 standard. Thus, among the control Enterococcus species strains, neither 16S rRNA gene sequencing nor BOX-PCR typing could discriminate among strains of the same species.
Phylogenetic tree constructed using the neighbor-joining algorithm to evaluate the distance between 16S rRNA gene sequences of environmental enterococci and known Enterococcus species. Numbers on the left are bootstrap values (percentages). The scale bar denotes substitutions per site.
Dendrogram comparing BOX-PCR patterns of Enterococcus species isolated from environmental matrices and known Enterococcus species. The species name to the right of the strain designation represents the identification based on 16S rRNA gene sequencing. The scale at the top of the figure represents percent similarity.
In contrast, BOX-PCR typing was better able to discriminate among environmental isolates of Enterococcus species than 16S rRNA gene sequencing (Fig. 1 and 2), which can be explained in part by the difference in methodology between 16S rRNA gene sequencing and BOX-PCR typing. The 16S rRNA gene sequence (∼1,500 bp) has both highly conserved and variable regions and is widely used for determination of evolutionary relationships (40). However, the resolution capacity of the 16S rRNA gene is limited when identifying closely related organisms (11, 33). 16S rRNA gene sequencing concentrates on a much smaller portion of the genome than BOX-PCR. BOX-PCR typing targets sequences located between interspersed repetitive DNA elements, resulting in amplification products of different sizes that generate a genomic fingerprint of individual bacterial strains. Variation in genome sizes (26), as well as the location of BOX elements among different strains of a particular Enterococcus species, leads to generation of multiple strain-specific fingerprint patterns, which allows BOX-PCR typing to discriminate among different strains of the same species better than 16S rRNA gene sequencing.
The 16S rRNA gene sequences of E. faecium and E. mundtii were approximately 98% similar and formed a single cluster in the phylogenetic tree. In contrast, the BOX-PCR patterns of E. faecium and E. mundtii exhibited less than 60% similarity. This demonstrates the ability of BOX-PCR genotyping to differentiate between closely related species of enterococci. Rademaker et al. compared genotypic relationships of Xanthomonas species and strains by BOX-PCR with DNA-DNA hybridization studies and observed a high correlation between the two methods (28). This observation reflects the good agreement between the genotypic (BOX-PCR) and phylogenetic (16S rRNA gene) relationships of Enterococcus spp. observed in our study. In contrast, other studies have found poor correlation between genotyping (rep-PCR) and DNA sequencing methods used to determine species/strain relationships of other bacterial genera (4, 19).
The 16S rRNA gene sequence of an environmental strain of Lactococcus garvieae served as an effective outgroup for the 16S rRNA gene dendrogram (Fig. 1). The 16S rRNA gene sequence of L. garvieae is 11.4 to 11.8% different from that of Enterococcus species (27). Although L. garvieae clustered outside most of the Enterococcus spp. in the BOX-PCR dendrogram at approximately 50% similarity, it did not form a true outgroup, and a heterogeneous assortment of Enterococcus species clustered with it, including one E. faecalis and one E. faecium strain (Fig. 2). However, visual comparison of the patterns demonstrates that the L. garvieae BOX-PCR pattern is more dissimilar from those of the enterococci than the value the software has calculated. These data illustrate a caveat that must be considered in all types of genetic comparisons, i.e., algorithms and software are useful for comparing multiple data entries that cannot readily be compared by eye, but each comparative method has tradeoffs that must be considered when interpreting the data. In the case of genomic patterns, software-generated comparisons should always be confirmed by eye (3, 24).
Of the 61 isolates sequenced in this study, only one was identified by 16S rRNA gene sequencing as a non-Enterococcus isolate. This finding confirms the specificity of mEI agar, a medium frequently used for isolation of enterococci from surface waters and recommended for ambient water quality testing by the U.S. Environmental Protection Agency (EPA) (35). In other studies, mEI agar was found to be less specific, since organisms belonging to other genera were isolated from environmental samples (23), biosolids (37), and clinical samples (13) along with Enterococcus spp. According to the U.S. EPA, the false-positive rate for isolation of nonenterococci from environmental samples on mEI agar is 6% (22, 35). In comparison, we observed a very low false-positive rate (1.6%) for isolation of nonenterococci from environmental matrices in the subtropical waters sampled in our study.
Although BOX-PCR genotyping was found to be more discriminatory at the strain level than 16S rRNA gene sequencing (e.g., E. casseliflavus/flavescens strains, E. hirae strains), the incorrect grouping of some strains within heterogeneous groups (e.g., the L. garvieae group) or with clades representing other species raises doubts about the ability of the method to correctly project relatedness of all Enterococcus strains (Fig. 2). While studies relying solely on BOX-PCR typing should exercise caution when inferring phylogenetic relationships, the method does provide a useful approximation of phylogeny and is an excellent tool for investigating strain diversity.
ACKNOWLEDGMENTS
We thank Paul Siddall (Lake Carroll resident) for the use of his property to access Lake Carroll. In addition, we thank Miriam Brownell, Asja Korajkic, Christopher Staley, and Mythili Penugonda (USF Department of Integrative Biology) for providing help with sample collection and processing. Partial funding was provided by U.S. EPA Gulf of Mexico Alliance Regional Partnership project MX-96478707-0.
FOOTNOTES
- Received 20 January 2011.
- Accepted 18 May 2011.
- Accepted manuscript posted online 27 May 2011.
- Copyright © 2011, American Society for Microbiology
