Previous Article | Next Article ![]()
Journal of Clinical Microbiology, August 2008, p. 2529-2534, Vol. 46, No. 8
0095-1137/08/$08.00+0 doi:10.1128/JCM.00813-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Departments of Microbiology,1 Infectious Diseases, Hospital Vall d'Hebron, Universitat Autònoma de Barcelona, Passeig Vall d'Hebron 119-129, 08035 Barcelona, Spain,2 Veterans Affairs Medical Center and University of Minnesota, Minneapolis, Minnesota 554173
Received 29 April 2008/ Accepted 5 May 2008
|
|
|---|
|
|
|---|
Favoring the prevalence hypothesis, in acute E. coli UTI, the urine strain tends to be the host's predominant fecal strain (4, 18, 26). Favoring the special-pathogenicity hypothesis, acute-phase UTI isolates usually exhibit more virulence factors than fecal isolates from healthy hosts and derive predominantly from phylogenetic groups B2 and D (versus groups A and B1 for fecal strains) (6, 13, 27). However, such comparisons may be confounded by between-population host differences, since women who develop UTIs tend to differ genetically or behaviorally from UTI-free women (19, 20, 23) and, conceivably, could have intrinsically different fecal E. coli populations. Consequently, the ideal comparison group for acute-phase UTI isolates may be fecal isolates from the same hosts.
Our previous comparison of same-host urine and fecal isolates from women with acute E. coli cystitis (12) suggested that the UTI strain is usually both more virulent-appearing and more abundant than coisolated fecal strains. However, the very small number of subjects (n = 11) and fecal colonies (3 per host) examined in that study greatly limited the strength of these conclusions. Furthermore, the limited sampling of the fecal population in that study and previous similar studies involving girls (15, 18), women (11), men (7), and dogs (5) with UTIs precluded valid assessment of the clonal diversity of the host's fecal population, which might be expected to vary with the nature of the colonizing clones and to influence the risk of UTI. Accordingly, in the present study we sampled a larger population of women with acute uncomplicated E. coli UTIs and analyzed 30 fecal E. coli colonies per host (i.e., 4 times as many subjects and 10 times as many colonies per host as in our pilot study). We then assessed the relative abundance, molecularly inferred virulence, and phylogenetic background of each fecal clone in comparison with those of the UTI-causing E. coli strain of the host and assessed bacterial traits in relation to the clonal complexity of the source fecal sample.
|
|
|---|
Urine microscopy and cultures. Quantitative urine culture was done using chromogenic agar, followed by conventional identification (14). E. coli isolates were stored in 5% glycerol in Trypticase soy broth at –70°C until use.
Stool cultures. Concurrent with urine collection, subjects collected a rectal swab using Amies transport medium. Rectal swabs were agitated in sterile saline, which was then diluted serially and spread onto MacConkey agar plates by using an automatic spiral spreading instrument. After 24 to 48 h of incubation (37°C), as many as 30 isolated colonies (as available) of putative E. coli were randomly picked. Colonies were identified using conventional methods (14). E. coli isolates were stored in 5% glycerol-supplemented broth at –70°C until use.
ERIC PCR.
Clonal relationships among the E. coli isolates from a given stool sample and the host's urine isolate were inferred from genomic enterobacterial repetitive intergenic consensus (ERIC) PCR profiles (25). Gel images were captured digitally (Gel Compar II, version 3.0; Applied Maths). Similarity relationships among profiles (within each subject) were estimated by using the Pearson correlation coefficient. Same-host isolates with
93% similar ERIC profiles were defined as clonally related. One representative of each putatively unique clone per fecal sample, plus each urine isolate, underwent phylotyping and virulence genotyping.
Phylotyping and virulence genotypes. The E. coli phylogenetic group (A, B1, B2, or D) was determined by using a triplex PCR method (2). Sixteen virulence-associated genes (papA, papG [alleles I, II, and III], papC, fimH, afa/draBC, sfa/focDE, hlyA, cnf1, iutA, fyuA, kpsM II, traT, ibeA, and malX) were characterized by multiplex PCR, with appropriate positive and negative controls (8). The virulence score was the number of virulence genes detected, with pap elements counting collectively as a single trait. (Such in vitro testing predicts experimental virulence in vivo [8, 17].)
Pulsed-field gel electrophoresis. E. coli isolates from a given subject that exhibited <93% similar ERIC profiles but identical results for the phylogenetic group and the virulence profile underwent pulsed-field gel electrophoretic analysis of XbaI-restricted total DNA (3). Isolates from a given subject that exhibited indistinguishable profiles (Dice similarity coefficient, 100%) were assumed to represent the same clone.
Definitions.
Each unique clone was assessed for its relative abundance in the host's fecal sample based on the proportion of fecal colonies it accounted for. The dominant clone was the sample's most abundant E. coli clone. If two clones were tied for the highest abundance in a sample, both were considered dominant. Pauciclonal fecal samples were arbitrarily defined as those with
4 E. coli clones; multiclonal samples were those with
5 E. coli clones.
Three distinct host ecological groups were defined on the basis of the relationship of each subject's fecal clone(s) with the urine clone (12). In group I, the stool culture yielded only the urine clone (designated a urine-fecal clone); in group II, it yielded the urine clone (urine-fecal clone) plus one or more nonurine (fecal-only) clones; and in group III, it yielded only one or more nonurine (fecal-only) clones, whereas the urine culture yielded a distinct (urine-only) clone.
Statistical methods. Comparisons of proportions were tested using a chi-square test or Fisher's exact test (two-tailed), as appropriate. Comparisons of deviations from a theoretically expected distribution were tested by using the binomial test. Virulence scores were compared by using the Mann-Whitney U test. Due to the nonindependence of multiple isolates acquired from the same sample, mixed-model regression modeling (21) was used. The mixed-effects model is a two-stage (hierarchical) maximum-likelihood approach. The first stage describes the distribution of the outcomes within subjects (i.e., the distribution is conditional on random effects). The second stage describes the variability across subjects and allows for the comparison of other effects (e.g., urine versus fecal source). The number of isolates per sample is treated as a random variable. The significance criterion was a P value of <0.05.
Correlations among variables were assessed by using correspondence analysis, which uses a covariance matrix based on
2 distances. 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. Correspondence analysis was conducted from a two-way table that had 109 rows, 1 per E. coli isolate, and 24 columns, 1 per variable.
Similarity relationships among the individual isolates with respect to virulence factor profiles were assessed by using principal-coordinates analysis (PCoA), a multivariate technique related to correspondence analysis that allows one to plot the major patterns within a multivariate data set, e.g., multiple loci and multiple samples (9). Using GenAlEx6 (16), PCoA was applied to the virulence factor data set as a way to collapse the multiple virulence factors for simplified among-group comparisons. Values for each clone from the first three PCoA axes, which capture most of the variance within the data set, were used in a one-way multivariate analysis of variance to test for differences between urine and fecal-only clones.
|
|
|---|
Distribution of E. coli clones by source and host. Overall, from the 42 subjects' fecal samples, 1,241 presumptive E. coli colonies were analyzed (30 colonies each for 38 subjects; 29, 28, 27, and 17 colonies each for 4 subjects, respectively). Of these, 1,117 colonies (90%) were confirmed as E. coli. According to molecular profiling, the 1,159 total E. coli colonies (1,117 fecal; 42 urine) represented 109 unique (by subject) clones. Of the 109 clones, 4 (4%) were found only in urine (urine-only clones), 38 (35%) in both urine and feces (urine-fecal clones), and 67 (61%) only in feces (fecal-only clones), for a total of 42 urine clones and 105 fecal clones.
Each subject had 1 urine clone and 1 to 11 (mean, 2.5) fecal clones. Comparison of each subject's urine clone and fecal clone(s) placed 14 subjects (33%) in ecological group I, 24 (57%) in group II, and 4 (10%) in group III (Table 1). Thus, in 90% (95% confidence interval [CI], 81% to 99%) of subjects (groups I and II), the urine clone was present in feces, with or without additional (fecal-only) clones.
|
View this table: [in a new window] |
TABLE 1. Phylogenetic distribution of 109 Escherichia coli clones from women with acute cystitis in relation to ecological group and source
|
Thirty-six (86%) of the 42 fecal samples yielded
4 clones each (i.e., were pauciclonal), collectively accounting for 66 (63%) of the 105 total fecal clones. The other six fecal samples (14%) yielded
5 clones each (i.e., were multiclonal), collectively accounting for the remaining 39 (37%) of the 105 fecal clones.
Phylogenetic group. Table 2 shows the overall distribution of phylogenetic groups among the 109 clones, with a rank order, from the most to the least abundant, of A, B2, D, and B1. Because some subjects had multiple clones, a by-subject analysis yielded higher values than these for the proportion of the 42 women with at least 1 clone from a given phylogenetic group, i.e., group A was found in 42% of women versus 38% of clones, group B1 in 21% versus 13%, group B2 in 71% versus 33%, and group D in 31% versus 16%. Group B2 showed the largest (absolute and proportional) increase in estimated prevalence in the by-subject analysis, surpassing group A to become the most prevalent (by-subject) phylogenetic group overall.
|
View this table: [in a new window] |
TABLE 2. Distribution of phylogenetic groups and virulence determinants according to source and the host's fecal clonal diversity among 109 Escherichia coli clones from women with acute cystitis
|
1 clone each, with overall prevalences ranging from 1% (afa/draBC) to 96% (fimH) (Table 2). Most individual virulence genes were significantly concentrated within groups B2 and D (data not shown). Accordingly, aggregate virulence scores were highest for group B2 (median, 7; range, 4 to 10), intermediate for group D (median, 5; range, 2 to 7), and lowest for groups A and B1 (median, 2; range, 0 to 4) (P < 0.001 for B2 versus each other group; P < 0.001 for group D versus group A or B1). Bacterial characteristics versus source and colonization behavior. Urine clones from ecological groups I, II, and III exhibited similarly high prevalences of group B2 and individual virulence genes, and similarly high aggregate virulence scores (see, e.g., Table 1). Likewise, fecal-only clones from groups II and III exhibited similarly low prevalences of group B2 (see, e.g., Table 1) and individual virulence genes, and similarly low aggregate virulence scores. Thus, all urine clones were combined for comparison with fecal-only clones (Table 2). Urine clones collectively exhibited significantly higher prevalences of group B2 and of multiple individual virulence genes, and significantly lower prevalences of groups A and B1, than the fecal-only clones. Likewise, the mean aggregate virulence scores of urine clones were significantly higher (mean scores, 6.2 [range, 1 to 10] versus 2.9 [range, 0 to 10]; P < 0.001).
Compared with nondominant fecal clones, dominant fecal clones more commonly represented the host's urine clone (i.e., were urine-fecal clones) (79% versus 24%; P < 0.001) and derived from group B2 (52% versus 15%; P < 0.001). They also exhibited higher virulence scores than other clones (5.2 versus 3.2; P < 0.001). Likewise, compared with clones from hosts with multiclonal fecal samples, clones from pauciclonal fecal samples exhibited a numerically higher prevalence of group B2, significantly higher prevalences of group D and of individual virulence genes (Table 2), a lower prevalence of group A, and higher aggregate virulence scores (mean score, 4.7, versus 3.0 for multiclonal fecal samples; P < 0.001).
Multilevel modeling. Due to the nonindependence among multiple isolates acquired from the same sample, we next employed mixed-model regression modeling (21). In the resulting univariate models for predicting urine clone status among the 109 total clones, similar associations were documented as in the comparisons that considered the clones independent. Specifically, phylogenetic group A was a significant (and group B1 a borderline significant) negative predictor of urine clone status (data not shown), whereas group B2 was a significant positive predictor of urine clone status (odds ratio, 20.75; 95% CI, 7.08 to 60.74) and was the strongest univariate predictor overall. Other significant positive predictors included 10 individual virulence traits, the aggregate virulence score (a stronger predictor than any single virulence gene [odds ratio, 1.62; 95% CI, 1.34 to 1.97]), fecal dominance (odds ratio, 11.53; 95% CI, 4.52 to 29.46), relative fecal abundance (odds ratio, 18.37; 95% CI, 5.03 to 67.11), and a pauciclonal fecal sample (odds ratio, 5.03; 95% CI, 1.84 to 13.73).
Next, multivariable mixed models for predicting urine clone status were constructed using all possible combinations of group B2 (as the phylogenetic-group variable), the virulence score (as the virulence trait variable), pauciclonal sample status, and either dominant clone status or absolute fecal prevalence (as alternate fecal-abundance variables). Those models in which all predictor variables remained statistically significant are shown in Table 3. When group B2 was included in the model, the virulence score lost significance, but each of the three colonization variables individually retained statistical significance and improved the predictive power of the model (i.e., lowered the –2 log likelihood value). Without group B2 in the model, the virulence score was statistically significant and remained so in the presence of each of the colonization variables, individually or in pairwise combinations, with the colonization variables also remaining statistically significant.
|
View this table: [in a new window] |
TABLE 3. Multivariable analysis of bacterial traits and colonization status as predictors of urine source according to multilevel modeling among 109 Escherichia coli clones from women with acute cystitis
|
![]() View larger version (11K): [in a new window] |
FIG. 1. Correspondence analysis for bacterial traits and colonization status among 109 Escherichia coli isolates. The 16 virulence traits found among the 109 E. coli isolates, plus the four phylogenetic groups (A, B1, B2, and D), the clone source (urine or fecal), the pauciclonality (Pauci) or multiclonality (Multi) of the fecal sample, dominant (Dom) clone status, and prevalence (Prev) in feces, as computed from the correspondence analysis, are projected on the F1-F2 plane. The F1 axis accounted for 26% of total variance and the F2 axis for 14.6%.
|
![]() View larger version (6K): [in a new window] |
FIG. 2. Distribution of 109 Escherichia coli isolates in relation to their source (urine or feces only) according to their values on the first axis of the PCoA.
|
|
|
|---|
This study is, to our knowledge, the first to assess a clone's quantitative fecal prevalence, or the pauciclonal versus multiclonal nature of the source fecal sample, in relation to the phylogenetic group, virulence factor profile, and UTI status, and to compare urine and fecal clones from the same hosts according to so broad an array of molecular traits by using sophisticated statistical methods such as multilevel modeling, correspondence analysis, and PCoA. Therefore, in addition to validating our previous findings (12) in a larger population and with greater methodological rigor, the present work also includes several novel features.
The study design allowed each woman to serve as her own control, thereby avoiding possible between-population differences that may confound comparisons of UTI isolates with fecal isolates from healthy hosts. This is an extremely powerful approach for removing the effects of known or unrecognized potential confounders, such as behavioral, environmental, physiological, or genetic differences between infected and uninfected hosts.
Our findings support three main conclusions. First, the urine clones, though almost always detectable in the host's fecal flora at the time of presentation, represented a highly biased subset of that flora, one enriched for traditionally recognized virulence traits and phylogenetic group B2, a result that supports the special-pathogenicity hypothesis. Second, however, the urine clone also tended to be the host's most abundant, or dominant, fecal clone, a result that supports the prevalence hypothesis. Finally, virulence traits and group B2 status were closely associated with fecal abundance, dominance, and pauciclonality. This suggests that prevalence and special pathogenicity are not alternative, mutually exclusive explanations for the occurrence of UTI. Instead, they may both be operational, perhaps contributing jointly to UTI pathogenesis, with virulence factors and other group B2-associated characteristics possibly promoting intestinal dominance (10) and thereby increasing the probability that subsequent steps in pathogenesis will occur.
In this regard, within the total clonal population, fecal abundance (and its close correlate, fecal dominance) was significantly associated with the virulence score (and its close correlate, group B2 status), and each of these domains individually was significantly associated with urine clone status. However, these domains were not precisely overlapping, since some independent predictive power for either domain was evident in the multivariable modeling. This suggests that virulence and abundance, although usually aligned, may sometimes substitute for one another in promoting UTI when either occurs alone. However, among the urine clones overall, dominance within the fecal reservoir did not appreciably reduce the requirement for virulence traits and group B2 status, evidence suggesting that the UTI-promoting effect of fecal abundance is generally insufficient of itself, in intact hosts, to allow low-virulence strains to cause UTI.
The observed association of fecal pauciclonality with group B2 status and virulence is novel and of potential biological and pathogenic significance. Why women with UTI who are colonized with many different E. coli clones tend to have primarily low-virulence, non-B2 clones, whereas those colonized with fewer E. coli clones tend to have high-virulence clones primarily from group B2, is unclear and warrants further study. It may be that certain women are particularly receptive to intestinal colonization with virulent clones, which then outcompete other clones and establish a low-diversity B2-dominated flora, whereas other women can support intestinal colonization with whatever clones they encounter, which probably are mostly low-virulence, non-B2 clones. Alternatively, this phenomenon could be a function of the strains themselves, with B2 strains tending to displace other phylogenetic groups through enhanced niche fitness or inhibition of other clones. Longitudinal studies could elucidate whether intestinal clonal complexity is a relatively stable trait specific to particular hosts or varies over time for a given host; whether expansion of minority B2 clones is typically accompanied by a reduction in overall intestinal clonal complexity; and whether antecedent intestinal clonal diversity predicts UTI risk.
Additional strengths of this study include the comparatively large number of subjects, the systematic approach to recruitment, and the clinical homogeneity of the study population. Limitations include the cross-sectional design, the use of multiple comparisons (which, however, was addressed by using multivariable analysis, correspondence analysis, and PCoA to supplement standard univariate analyses), and the absence of vaginal and periurethral sampling.
In summary, our findings suggest that the relationship between the structure of the fecal E. coli population and UTI pathogenesis is complex. Fecal pauciclonality, clonal dominance or relative abundance, molecularly inferred virulence, and group B2 status are closely intertwined. Phylogenetic group B2 status and/or virulence factors may promote fecal abundance and pauciclonality, thereby contributing to upstream steps in UTI pathogenesis. This relationship, if confirmed, suggests a possible reconciliation of the prevalence and special-pathogenicity hypotheses. Furthermore, pauciclonality of the intestinal E. coli population may represent a previously unrecognized correlate of pathogenicity deserving of further investigation.
We thank the study participants for their contributions.
Published ahead of print on 21 May 2008. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»