Previous Article | Next Article ![]()
Journal of Clinical Microbiology, December 2005, p. 5860-5864, Vol. 43, No. 12
0095-1137/05/$08.00+0 doi:10.1128/JCM.43.12.5860-5864.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Division of Infectious Diseases and Immunity, University of California at Berkeley, Berkeley, California,1 Department of Integrative Biology, University of California at Berkeley, Berkeley, California,2 Department of Epidemiology, Biostatistics, and Occupational Health, McGill University, Montreal, Quebec, Canada,3 Division of Epidemiology, School of Public Health, University of California at Berkeley, Berkeley, California4
Received 11 May 2005/ Returned for modification 8 August 2005/ Accepted 26 September 2005
|
|
|---|
|
|
|---|
The genotypic characterization of pathogens has become an important objective in epidemiologic investigations of infectious agents. Genotyping tools are amenable to validation when the diseases in question occur as recognizable outbreaks, such as those caused by E. coli O157:H7 or Salmonella spp. The problem arises in characterizing organisms that cause diseases that are not usually thought of as causing outbreaks, such as E. coli associated with community-acquired UTI. Recent studies that applied genotyping methods to differentiate uropathogenic E. coli suggest that community-acquired UTIs could occur as outbreaks.
Manges et al. reported in 2001 that a single clonal group of E. coli, designated clonal group A (CgA), based on a characteristic enterobacterial repeat intergenic consensus (ERIC2) PCR fingerprint pattern, was responsible for nearly half of trimethoprim-sulfamethoxazole-resistant E. coli isolates from 47 women with UTIs diagnosed at a California university health clinic between October 1999 and January 2000 (13). Human UTI isolates belonging to CgA were also found to share a random amplified polymorphic DNA PCR fingerprint, virulence factor profile, similar antimicrobial resistance phenotype, O antigen groups (O11, O17, O73, and O77), and indistinguishable or closely related pulsed-field gel electrophoresis (PFGE) patterns (9), suggesting that they caused an outbreak in this California community. Most recently, a gene-specific PCR was devised to identify CgA based on a single-nucleotide polymorphism (SNP) (C288T) identified within the fumC housekeeping gene (11).
Certain genotyping methods (PCR-based methods and PFGE) rely on comparison of banding patterns generated by gel electrophoresis. Such techniques are condition dependent and therefore can suffer from interlaboratory variability (20). Furthermore, visual inspection remains integral to their interpretation, introducing subjectivity and greater potential for error (22). While PFGE does offer superior reproducibility and discriminating power and remains the method most commonly used to identify food-borne outbreaks of E. coli O157:H7 at reference laboratories, it is labor intensive and technically demanding (8, 13, 21). Although gene-specific PCR is highly reproducible, a recent study evaluating 45 isolates exhibiting an ERIC2 PCR CgA pattern found that only 16 isolates contained the C288T SNP (4). Thus, a more reliable, relevant, and consistent typing method for investigating the clonal distribution of uropathogenic E. coli strains, such as CgA isolates, is needed for epidemiologic investigations. The availability of a large, well-characterized, and population-based collection of uropathogenic E. coli isolates provided us with an opportunity to evaluate another genotyping method.
Multilocus sequence typing (MLST) uses nucleotide sequences of internal fragments of selected genes as the unit of comparison and, therefore, does not suffer from the drawbacks of gel-based fingerprinting methods. Sequence data are easily comparable and transferable between laboratories and are highly reproducible (2, 12). Furthermore, the digital format of MLST data has facilitated the establishment of global, web-accessible databases for a variety of organisms and is rapidly contributing to our understanding of the clonal distribution of infectious disease agents. The most commonly used MLST schemes index the neutral genetic variation in housekeeping genes, which are believed to evolve slowly because they are under stabilizing, and not directional, selective pressure. The SNPs observed across several different MLST loci represent neutral genetic variation and presumably exhibit minimal autocorrelation, which could confound epidemiologic conclusions. MLST is thus a powerful tool for global and long-term surveillance.
In this study, we evaluated an MLST protocol standardized for E. coli maintained at the Max-Planck Institut fuer Infektionsbiologie website (http://web.mpiib-berlin.mpg.de) using a well-characterized population-based collection of CgA and other E. coli strains and compared its discriminatory power with that based on PFGE, ERIC2, and serogroup analyses.
|
|
|---|
Isolates included strains from the following sources: (i) UTI isolates from California (n = 15) women with symptoms of UTI who were seen at a university health service (they were consecutively enrolled into a study between 11 October 1999 and 31 January 2000) and UTI-causing, trimethoprim-sulfamethoxazole-resistant E. coli isolates (n = 6) from Minnesota (obtained from students who were seen at university health services with uncomplicated cystitis and were enrolled in a study between June 1998 and August 1999); (ii) animal and environmental isolates (n = 6) provided by the Gastroenteric Disease Center at Pennsylvania State University (University Park, PA); (iii) human bacteremia isolates (n = 12) from San Francisco General Hospital provided by Francoise Perdreau-Remington; and (iv) other geographic comparison isolates recovered from humans (n = 6) provided by James R. Johnson of the Minneapolis VA Medical Center, University of Minnesota.
Serotyping. Serotyping was performed on E. coli isolates at the E. coli Reference Center in University Park, Pennsylvania.
Molecular strain typing analyses. All E. coli isolates were initially typed by the ERIC2 PCR fingerprinting assay, as described elsewhere (10, 13). The CgA ERIC2 PCR electrophoretic pattern was defined by 4 predominant bands of approximately 1,145, 1,029, 908, and 720 bp; isolates exhibiting this pattern were considered to be members of CgA. A pyelonephritogenic isolate, CFT073 (provided by Harry Mobley, University of Maryland, Baltimore), was used as a reference strain for each ERIC2 PCR run.
PFGE. The standardized protocol for subtyping E. coli O157:H7 by PFGE, as established by the Centers for Disease Control and Prevention (Atlanta, GA), was used (1). XbaI-digested DNA was electrophoresed in the CHEF DR-II apparatus (Bio-Rad, Hercules, CA). Images of PFGE and ERIC2 PCR electrophoretic patterns were imported and analyzed with GelCompar II, version 2.0 (Applied Maths, Kortrijk, Belgium), with the pattern of the prototype CgA strain (SEQ102, deposited in the American Type Culture Collection as ATCC BAA-457) used as a reference, as previously described (19). To minimize gel-to-gel banding pattern variation, we used the autosearch band calling function of GelCompar with the following parameters: minimum profiling, 5%; gray zone, 0%; minimum area, 1%; shoulder sensitivity, 1. A distance matrix was calculated by GelCompar's Dice algorithm with a band position tolerance of 1%. A dendrogram was generated from the distance matrix by the neighbor-joining method.
MLST. (i) Allele templates selected. Housekeeping genes for typing were selected from the E. coli database at the MLST website maintained at the Max-Planck Institut fuer Infektionsbiologie (http://web.mpiib-berlin.mpg.de) (Table 1). The genes were shown to be unlinked on an E. coli K-12 genome map. Product lengths varied from 583 to 932 bp (Table 1). Allele templates, included on the website, were based on the genome sequence of E. coli strain MG1655.
|
View this table: [in a new window] |
TABLE 1. Forward and reverse sequences of primersc
|
PCR. Amplifications were carried out in a total volume of 50 µl with 2 µl of template DNA, 4 µl of each 10 mM primer (Sigma, Genosys, The Woodlands, TX), 5 µl 10x buffer, 15 mM MgCl2 (Applied Biosystems, Foster City, CA), 5 µl 2 mM deoxynucleoside triphosphate mix (Invitrogen, Carlsbad, CA), and 0.5 U AmpliTaq Gold (Applied Biosystems). Reaction conditions used were: 2 min of denaturation at 95°C; 30 cycles of 1 min of denaturation at 95°C, 1 min of annealing at specific temperature (Table 1), 2 min of extension at 72°C; and a final additional 5 min at 72°C in an Applied Biosystems GeneAmp PCR System 2400 thermocycler.
Sequencing. PCR products were purified for sequencing with the QIAGEN QIAquick PCR purification kit. Both the forward and reverse strands were sequenced with the PCR primer set, with the exception of mdh, which was sequenced by the forward primer 5'-TCTGGTGAAGATGCGACTCC-3' and reverse primer 5'-CCCAGGGCGATATCTTTCTT-3'. Sequencing was performed at the University of California at Berkeley Sequencing Facility, which is equipped with an Applied Biosystems 48 capillary 3730 and two Applied Biosystems 16 capillary 3100 genetic analyzer systems. A Perkin Elmer (Wellesley, MA) 9600, an Eppendorf (Hamburg, Germany) Mastercycler, and three MJ Research (Bio-Rad, Hercules, CA) PTC-200s were used for cycle sequencing. The facility runs a 25-cycle sequencing reaction with the following program: 96°C for 10 s plus 50°C for 5 s plus 60°C for 4 min.
Sequence analysis. Raw sequences were reviewed by visual inspection with Chromas, version 1.45 (32-bit) (Technelysium Pty. Ltd., Tewantin Qld 4565, Australia). DNA sequences were aligned by the neighbor-joining method with 1,000 bootstrap iterations in ClustalW (Hinxton, England). Forward and reverse sequences were aligned for each strain for comparison and editing purposes. Alignments were edited to equally sized fragments and aligned against allele templates based on E. coli strain MG1655, provided at http://www.mlst.net. Sequence editing was conducted with BioEdit (Carlsbad, CA), version 7.0.1. Sequences for the seven genes of each strain were concatenated to produce an alignment sequence of 3,423 bp. A dendrogram of SNP groupings based on these alignments was then constructed by the Phylip programs (Seattle, WA) DNADIST and NEIGHBOR, and the final dendrogram was visualized with the software TreeView (Win32) (Glasgow, Scotland), version 1.6.6.
ST designation. Gene sequences were submitted to the curator of the MLST E. coli submission page maintained at http://web.mpiib-berlin.mpg.de. Sequence types (STs) and sequence complexes were designated by the curator of the MLST database for sequence submissions that contained tracings found to be acceptable by the curator (Table 2). ST complexes were defined as STs that are linked by distances of one to two allelic differences.
|
View this table: [in a new window] |
TABLE 2. STs and ST complexes
|
|
|
|---|
![]() View larger version (73K): [in a new window] |
FIG. 1. Neighbor-joining tree constructed from the concatenated sequences of the 7 MLST genes described in the text. The 29 previously identified CgA isolates are shown with shaded backgrounds. Serogroups are indicated for the CgA isolates. The inset shows a dendrogram based on Dice distance coefficient measurements of PFGE banding patterns among the 21 isolates belonging to the ST69 complex. *, no PFGE data available; **, ST69 complex members, not ST69 (see Table 2); , prototype human uropathogenic CgA strain ATCC BAA-457.
|
Of 22 ST69 complex isolates, 19 had identical sequences in all 7 gene templates. These 19 ST69 strains were further separated into 14 subtypes by PFGE (Fig. 1, inset).
|
|
|---|
Interestingly, strain AN559, isolated from a cow in 1988 that was previously shown to be 94% similar by PFGE to a human UTI-causing CgA E. coli isolate (19), was found to occur within the ST69 complex. This was the only animal strain to cluster within this group, differing by a single nucleotide in the adk gene.
Recently, CgA isolates shown to be highly similar by multiple techniques, including C288T SNP analysis, drug susceptibility testing, presence of the 1.8-kb class I integron bearing the dfrA17 gene, and virulence gene profile analysis, were each shown to be distinct by PFGE profile (4). Our investigation of the 21 CgA strains evaluated by PFGE showed that there was no difference in clustering by MLST among the strains that shared PFGE types and the strains that were unique by PFGE. Thus, PFGE still has superior discriminatory ability for uropathogenic E. coli than any of the other methods, including this MLST scheme. A similar observation was made by Noller et al., who compared PFGE and MLST for their discriminatory ability using a collection of diarrheagenic E. coli serotype O157:H7 (16). In one study, MLST was found to be more discriminatory than PFGE only if gene templates in addition to the 7 housekeeping genes were included (15).
ST69 belongs to ECOR D group. A previous study of CgA strains by virulence factor profile analysis found them to have a profile similar to that of the O15:K52:H1 clonal group that was reported to be responsible for community outbreaks of UTI and bacteremia in Europe (8, 17, 18).
Ultimately, the appropriate level of discriminatory power of any strain typing test is determined by the epidemiologic question of interest. Just as highly discriminatory typing methods can obscure important epidemiologic relationships between related strains, typing schemes with low discriminatory power can obscure important epidemiologic associations. Uropathogenic CgA E. coli strains clearly include multiple subtypes, as was shown by the present MLST technique, PFGE analysis, and analysis by France et al. (4). Here, a housekeeping gene-based MLST applied to a well-characterized CgA collection demonstrated discriminatory power that fell between ERIC2 PCR and PFGE. Therefore, for an MLST scheme to be applicable for community-based microepidemiologic studies of uropathogenic E. coli, it will have to be based on allelic templates that show a greater degree of variation.
This work was supported by NIH grant RO1 AI059523.
|
|
|---|
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»