Previous Article | Next Article ![]()
Journal of Clinical Microbiology, February 2008, p. 417-422, Vol. 46, No. 2
0095-1137/08/$08.00+0 doi:10.1128/JCM.00674-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Mucosal and Vaccine Research Center,1 Geriatric Research, Education, and Clinical Center, Veterans Affairs Medical Center,3 Departments of Medicine,2 Psychiatry, University of Minnesota, Minneapolis, Minnesota,4 E. coli Reference Center, Department of Veterinary and Biomedical Sciences, Pennsylvania State University, University Park, Pennsylvania5
Received 26 March 2007/ Returned for modification 8 June 2007/ Accepted 1 November 2007
|
|
|---|
|
|
|---|
However, these studies examined only modest numbers of isolates and accessory traits and/or relied on somewhat imprecise phylogenetic methods such as multilocus enzyme electrophoresis, outer membrane protein profiling, or random amplified polymorphic DNA analysis (3, 5, 6, 8, 9, 11, 13, 15, 16, 20, 22, 26, 30-33). Accordingly, they leave uncertainty as to the extent of commonality among human- and pet-derived E. coli isolates.
To examine more rigorously the question of human-pet commonality with respect to extraintestinal pathogenic E. coli (ExPEC) clones, we used contemporary molecular methods to characterize a large collection of archival E. coli isolates of serotypes O6:H1, O6:H7, and O6:H31 from humans, dogs, and cats, obtained from the E. coli Reference Center at Pennsylvania State University. Specifically, we assessed broad phylogenetic relationships by using PCR-based phylotyping and multilocus sequence typing (MLST), determined accessory trait profiles through PCR-based detection of 55 virulence-associated genes and variants, and defined clonal identity by pulsed-field gel electrophoresis (PFGE). We then assessed the extent to which these traits segregated by host species or were shared among host species, especially between humans and pets.
|
|
|---|
Serotyping. O typing was performed as described by Orskov and Orskov (23). H typing was performed using a fliC PCR-restriction fragment length polymorphism method (21).
Phylogenetic analysis and MLST.
The major E. coli phylogenetic group (A, B1, B2, or D) was determined by a multiplex PCR assay (4). Closer phylogenetic relationships were assessed according to maximum parsimony, using PAUP* (28), based on the partial DNA sequence (
500 bp per locus) for fumC (all isolates) or a concatenated sequence including adk, fumC, gyrB, icd, mdh, purA, and recA (several representatives of each fumC allele, selected for diversity of H types and host species) (12). Additionally, designations for alleles at each of these loci, and sequence types (STs) reflecting specific allele combinations, were assigned according to the Achtman MLST system (http://web.mpiib-berlin.mpg.de/mlst/dbs/Ecoli). Strains MG1655 (K-12; group A), 536 (O6:K15:H31; group B2), and CFT073 (O6:K2:H1; group B2) and Escherichia fergusonii were included for reference (12).
Virulence genotypes. Fifty-five ExPEC-associated VF genes and variants were detected by using established multiplex PCR-based assays (10, 14, 17) (see Table S1 in the supplemental material). A dendrogram based on extended VF profiles was inferred by using the unweighted-pair group method with arithmetic averaging (UPGMA) (27).
PFGE.
Animal and human isolates that belonged to the same ST and exhibited identical VF profiles underwent PFGE analysis using XbaI according to a standardized protocol (24a). Dendrograms based on Dice similarity coefficients were constructed from the resulting profiles according to UPGMA. Isolates with
94% similar profiles (corresponding approximately with a
3-band difference) were considered to represent the same pulsotype, or clone (29).
Statistical methods. Two multivariate ordination techniques, correspondence analysis (CA) and principal-coordinates analysis (PCA), were used to visualize the arrangement, or "ordering," of data units (either variables or cases) along gradients.
Correlations among variables (host species, ST, H type, and individual VFs) were assessed by using CA. A special case of canonical correlation, CA is an eigenanalysis ordination technique for visualizing the associations in contingency tables (7). It produces a display of two-way contingency table cells as points in a low-dimensional space. The positions of the points provide a spatial map reflecting their associations within the table and enabling a global view of the data that is useful for interpretation. The computation determines a plane defined by two principal axes of the analysis. The first axis, F1, accounts for most of the variance, and the second axis, F2, orthogonal to F1, accounts for the largest part of the variance not accounted for by F1. Here, CA was conducted from a two-way table that had 130 rows (1 per E. coli isolate) and 59 columns (1 per variable).
Similarity relationships among the individual isolates with respect to VF profiles were assessed by using PCA. Also known as metric multidimensional scaling, PCA is a distance-based ordination technique that allows one to plot similarities on underlying dimensions characterizing a multivariate data set, e.g., multiple loci and multiple samples, in a low-dimensional space where distances between points are directly related to their similarity (18). Using GenAlEx6 (24), PCA was applied to the VF data set as a way to collapse the multiple VFs for simplified among-group comparisons. Values for each isolate from the first two PCA axes, which are the axes that capture most of the variance within the data set, were used in a one-way multivariate analysis of variance to test for pairwise between-group (ST or host species) differences. The values also were plotted to depict spatially the degree of separation or overlap among groups on the axis 1-axis 2 plane.
Comparisons of proportions were tested using Fisher's exact test. In assessing the distribution of categorical traits among multiple groups, an initial overall Fisher exact test screen was used to identify the presence of a significant among-group difference. The criterion for statistical significance was a P value of <0.05.
|
|
|---|
Maximum-parsimony analysis of partial fumC sequence data supported this inference (Fig. 1). The 129 group B2 study isolates segregated into three distinct groups, largely by H type. According to the Achtman MLST system, the three groups corresponded with fumC alleles 24 (n = 39) (92% H1), 43 (n = 22) (82% H7), and 14 (n = 68) (91% H31). Three strains (a cat H1, a dog H1, and a cat H7 strain) exhibited fumC allele 24 except for a 6-bp in-frame deletion, so they were analyzed with allele 24.
![]() View larger version (28K): [in a new window] |
FIG. 1. Phylogeny of 130 Escherichia coli isolates of serogroup O6 from humans, dogs, and cats according to fumC (left) (total population) or adk, fumC, gyrB, icd, mdh, purA, and recA (right) (representatives of each fumC allele). Trees were inferred using maximum parsimony. Bootstrap values of >60% (from 1,000 iterations) are shown. In the fumC tree (left), fine tick marks indicate individual study isolates with identical sequences. In the composite tree (right), each study isolate's unique identification number and (7-locus) sequence type (ST) are given. In both trees, study isolates are labeled with the host species (H, human; D, dog; C, cat) and H type (1, 7, or 31). According to the Achtman scheme, fumC alleles 14 (91% H31), 24 (92% H1), and 43 (82% H7) correspond to STs 127, 73, and 625, respectively. Reference strains 536, J96, CFT073, RS218, and MG1655 are included for comparison.
|
To corroborate the fumC phylogeny, 4 study isolates each per major fumC allele, selected to reflect the diversity of hosts and H types encountered for that allele, underwent partial sequence analysis for six additional housekeeping genes, as did the putative group A O6 isolate. A maximum-parsimony tree based on concatenated sequences across all seven loci exhibited much the same overall conformation as did the fumC tree, but with greater group separation (Fig. 1). Within each of the three main fumC groups, all isolates tested exhibited identical sequences across the seven loci, evidence that the fumC alleles correspond with 7-locus STs. According to the Achtman MLST scheme, the three main STs observed among the O6 isolates were ST73 (H1), ST625 (H7), and ST127 (H31). Reference strains CFT073 (ST73) and 536 (ST127) clustered within STs corresponding to their respective H types, whereas other reference strains represented distinct STs. Only for one study isolate did 7-locus MLST significantly alter the phylogenetic inference based on analysis of fumC alone. This was the putative group A O6 strain, which according to MLST was a chimera that exhibited ST127-like alleles of icd, gyrB, mdh, purA, and recA, a K-12 (ST98)-like allele of fumC, and a distinct (i.e., non-K12, non-ST127) adk allele.
CA. To explore associations among bacterial characteristics in relation to the host species and ST, CA was done. The F1/F2 plane accounted for 30% of total variance, with factors F1 and F2 accounting for 16% and 14%, respectively (Fig. 2). The other factors (not shown) accounted for lower percentages of the variance than did F1 and F2 and did not further improve the interpretation of the data.
![]() View larger version (16K): [in a new window] |
FIG. 2. CA for bacterial traits and host species among 130 Escherichia coli serogroup O6 isolates. Shown is a projection of the 55 VFs (named as in Table S1 in the supplemental material) characterized for the 130 E. coli isolates plus the three H types, three STs, and three host species (human, dog, and cat), on the F1/F2 plane, as computed from the CA. The loading score of each variable for factors F1 and F2, respectively, can be inferred from its x and y coordinates on the F1/F2 plane. I, II, and III refer to the corresponding papG (P fimbrial adhesin) alleles.
|
These inferences regarding the associations of VFs with host species and STs were supported by univariate comparisons of proportions, which showed more numerous and statistically stronger associations of individual VFs with STs than with host species (see Tables S1 and S2 in the supplemental material). Notably, no VF was associated categorically with a single host species except for cvaC and iss, associated with cats (P = 0.048 and P = 0.008, respectively), or with humans versus dogs and cats combined except for the F13 papA allele (P < 0.001) (see Table S2 in the supplemental material).
PCA. PCA was used to analyze jointly the 55 VFs for among-population comparisons involving host species and STs. The first two axes of the PCA collectively accounted for 64% of total variance, with axes 1 and 2 accounting for 43% and 21%, respectively. Based on each isolate's values on these two axes, pairwise comparisons were made between different host species and, separately, between different STs, by using a one-factor multivariate analysis of variance.
In the analysis by host species, human isolates differed significantly from dog and cat isolates (P < 0.001 for each comparison), whereas dog and cat isolates did not differ significantly from each other. However, this stratification of the PCA results yielded an eta-squared value of only 0.14. Likewise, in the plot of the axis 1-axis 2 plane, human, dog, and cat isolates overlapped considerably, despite some apparent segregation (Fig. 3).
![]() View larger version (10K): [in a new window] |
FIG. 3. Distribution of 130 Escherichia coli serogroup O6 isolates from humans, dogs, and cats on the axis 1-axis 2 PCA plane. Data used were extended virulence genotypes (55 traits). Top and bottom plots contain identical data (one point per isolate), with points colored to indicate either host species (top) or ST (bottom). Axis 1 (positive values to the right and negative values to the left of the dashed vertical line) accounted for 43% of total variance. Axis 2 (positive values above and negative values below the dashed horizontal line) accounted for 22% of total variance. The axes are dimensionless, i.e., they have no units.
|
These conclusions were supported by analysis of a similarity dendrogram based on extended VF profiles, which exhibited seven prominent clusters, each representing a distinct VF profile group (see Fig. S4 in the supplemental material). The three host species were extensively intermingled within the seven VF profile clusters, whereas the three STs segregated largely by cluster.
Similarity of human and animal isolates by PFGE. To identify possible cross-species clonal commonality, human and animal isolates from the same ST that also exhibited identical VF profiles underwent comparative PFGE analysis. Multiple PFGE profiles were encountered per ST, evidence of greater discriminating power for PFGE than for MLST (Fig. 4). In the resulting PFGE-based dendrogram, isolates from ST74 (H1) and ST127 (H31) were resolved as separate clusters (Fig. 4). Within each cluster, certain animal isolates were closer to one or more human isolates than to other animal isolates, and vice versa. In one instance, a dog-human isolate pair (isolates 26 and 34) exhibited the same PFGE profile (97% similarity; PFGE type 435) and hence presumably represented the same clone.
![]() View larger version (38K): [in a new window] |
FIG. 4. PFGE profiles of Escherichia coli serogroup O6 isolates from humans (H), dogs (D), and cats (C) with identical STs and virulence profiles. The dendrogram was inferred according to UPGMA based on Dice similarity coefficients as calculated from PFGE band positions. Profiles exhibiting 94% similarity (boxed) were regarded as representing the same PFGE type. Isolates D026 and H034 (PFGE type 435; ST127) were submitted in 1985. They were isolated, respectively, from a dog in Florida and a woman in New York, each with a UTI.
|
|
|
|---|
The three main STs (as resolved similarly by fumC sequence analysis and 7-locus MLST) corresponded closely with the three H (flagellar) types studied: H1, H7, and H31. However, each ST included several isolates that exhibited alternate H types, evidence of lateral transfer of fliC alleles among lineages. This finding confirms, using a more robust comparison standard than that used by Cherifi et al., these investigators suggestion that within serogroup O6 the H antigen is a useful clonal marker (3). However, it also demonstrates that H typing has limitations in this role.
Although some segregation of host species by ST was apparent, the three main STs each included isolates from all three host species, suggesting that the STs represent broad-host-range clonal groups. Of note, ST73 and ST127 appear to correspond, respectively, with clones 1 and 2 of Cherifi et al. (3) and with multilocus enzyme electrophoretic types 1 and 21 of Whittam and colleagues (9, 31). These represent two prominent clonal groups within E. coli serogroup O6 that are common causes of extraintestinal infections in humans, dogs, and cats. Underscoring the pathogenic importance of ST73 and ST127, these two clonal groups include, respectively, reference strains CFT073 and 536, two archetypal ExPEC strains from humans with pyelonephritis that have been extensively studied to elucidate causal mechanisms in UTI and other extraintestinal E. coli infections (12).
We also found that human and animal isolates overlapped considerably according to aggregate VF profiles and individual VFs, notwithstanding some statistically significant differences. Although it is conceivable that unrecognized host group-specific VFs were omitted from our analysis, the present findings suggest that similar VFs and VF profiles are relevant for infections in multiple host species, further supporting concerns regarding the possible zoonotic potential of such strains.
To more stringently test the hypothesis of cross-species commonality of virulent E. coli clones, we used PFGE analysis to resolve finer clonal relationships among animal and human isolates from the same ST that exhibited identical VF profiles. The close similarity of PFGE profiles observed among several human and animal isolates, including one 97% similar dog-human isolate pair, further supports the possibility of pet-human exchange of virulent E. coli clones. This is consistent with the results of two recent longitudinal household surveillance studies in which multiple humans and their pets were found to be colonized with the same E. coli clones, including two cystitis-causing strains (8, 22). Notably, in one such study the most extensively shared clone, which colonized both human subjects and the family cat and caused an episode of acute cystitis in the woman, belonged to serogroup O6 and, based on the VF profile and random amplified polymorphic DNA analysis, likely derived from ST127 (22).
The present study supersedes previous similar work (3, 5, 6, 8, 9, 11, 13, 15, 16, 20, 22, 26, 30-33) in several respects. It involves more isolates and a more extensive battery of VFs than did previous studies and is the first to use sequence-based phylogenetic analysis. It also is the first to demonstrate such a close PFGE profile match between epidemiologically unrelated human and pet E. coli isolates. Limitations include the undefined clinical background of the isolates, the absence of epidemiological evidence regarding cross-species transmission, and the lack of attention to minor sequence variation within virulence genes that might confer host specificity.
In summary, we found extensive commonality among human, dog, and cat E. coli isolates of serogroup O6 with respect to phylogenetic background, VF profiles, and clonal identity. Although some segregation by host species was apparent, the findings support the hypothesis that humans and their domestic pets can be colonized and infected with similar E. coli strains and may serve as reservoirs of potential pathogens for one another. Thus, in some instances extraintestinal E. coli infections in humans, dogs, and cats may represent animal-to-human or human-to-animal zoonoses.
Dave Prentiss (VA Medical Center) prepared the figures.
Published ahead of print on 14 November 2007. ![]()
Supplemental material for this article may be found at http://jcm.asm.org/. ![]()
|
|
|---|
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»