High rate of mobilization for [bla.sub.CTX-M]s.

We constructed a phylogenetic analysis of class A [beta]-lactamases and found that the [blac.sub.CTX-M.]s have been mobilized to plasmids [approximately equal to] 10 times more frequently than other class A [beta]-lactamases. We also found that the [bla.sub.CTX-M]s are descended from a common ancestor that was incorporated in ancient times into the chromosome of the ancestor of Kluyvera species through horizontal transfer. Considerable sequence divergence has occurred among the descendents of that ancestral gene sequence since that gene was inserted. That divergence has mainly occurred in the presence of purifying selection, which indicates a slow rate of evolution for [bla.sub.CTX-M]s in the pre-antimicrobial drug era.

**********

Antimicrobial drug--sensitive bacteria become resistant to antimicrobial drugs through a variety of mechanisms, such as chromosomal mutations that up-regulate the expression of antibiotic-resistance genes, DNA uptake through transformation, or the process of conjugation. The ability of plasmids to evolve independently of their hosts has allowed numerous resistance genes from diverse species of bacteria to assemble within single plasmids and spread into a wide variety of organisms (1). The mobilization of a chromosomal resistance gene to a plasmid is an important event because the mobilized gene is now capable of spreading widely throughout diverse species of bacteria and because the fitness advantage that a plasmid confers generally increases as it acquires more resistance genes (1).

The class A [beta]-lactamases have been the most frequently encountered plasmidic resistance genes. Class A [beta]-lactamases from the TEM group have occurred at a particularly high frequency; in many surveillance studies, they have been identified as the resistance determinants most frequently encountered (2-9). The first [bla.sub.TEM] allele, [bla.sub.TEM-1], is a plasmidic allele that was first isolated in 1963 (10,11). Currently, [approximately equal to] 160 different plasmidic allele s encode unique TEM [beta]-lactamase enzymes (www.lahey.org/Studies), and all are descended from a single plasmidic ancestor, [bla.sub.TEM-1] (12).

The SHVs are another group of class A [beta]-lactamases that have been frequently observed in surveillance studies. As with the TEMs, numerous alleles encode unique SHV enzymes ([approximately equal] 105), and the SHVs are all descended from a single ancestor (13). The first [blas.sub.SHV] allele was detected in 1974 (10,14). Unlike [bla.sub.TEM]s, the [bla.sub.SHV]s are present in the chromosome of nearly all Klebsiella pneumoniae isolates belonging to the KP1 group. Evidence suggests that [bla.sub.SHV]s have been chromosomally located since the pre--antimicrobial drug era (15), and they may have been mobilized to plasmids up to 4 times, although the sequence divergence among them is insufficient to clearly resolve the independent mobilizations of the [bla.sub.SHV]s.

The CTX-Ms are another group of class A [beta]-lactamases that are located on plasmids and that have been of particular clinical interest because they are rapidly spreading through clinical populations of bacteria. The first plasmidic [bla.sub.CTX-M] observed in human-associated clinical populations was isolated in 1989 (16,17). Unlike the usual pattern of class A [beta]-lactamase mobilizations in which the plasmidic alleles are all descended from a single common plasmidic ancestor, evidence shows that CTX-Ms have been mobilized numerous times from the chromosomes of Kluyvera (16,18-21). Because Kluyvera chromosomal genes have been found that exactly match the sequence of plasmidic CTX-Ms (18), many of the mobilizations have likely occurred recently. To investigate the mobilization of the [bla.sub.CTX-M]s to plasmids, we generated a phylogenetic analysis of the CTX-Ms that included a representative sampling of other class A [beta]-lactamases.

Methods

BLAST Search

[bla.sub.CTX-M] and [bla.sub.CTX-M] homologue DNA sequences were identified with a TBLASTN (www.ncbi.nlm.nih.gov/blast) (22,23) search of the nonredundant National Center for Biotechnology Information (NCBI) sequence database and the completed microbial genomes database by using characterized [bla.sub.CTX-MS] vs as query sequences ([bla.sub.CTX-M1] and [bla.sub.CTX-M2]). The BLAST search of the completed microbial genomes identified positive matches for organisms that contain [bla.sub.CTX-M] homologs. The BLAST search of completed genomes also showed which microbes have no close [bla.sub.CTX-M] homologue, and thus enabled horizontal transfer events to be identified.

Alignment

The protein sequences of the [bla.sub.CTX-M]s and their homologs were aligned with ClustalX 1.8 (24) by using the Gonet 250 similarity matrix with a gap-opening penalty of 35 and a gap-extension penalty of 0.75 for the pairwise alignment stage, and a gap-opening penalty of 15 and a gap extension penalty of 0.3 for the multiple alignment stage. The corresponding DNA coding sequences were aligned by introducing triplet gaps between codons corresponding to gaps in the aligned protein sequences with the CodonAlign program. (CodonAlign for Macintosh and for PC [Windows] computers and source code that can be compiled for other platforms are available at no charge from http://sinauer.com/hall.)

Estimation of Positive Selection

Estimation of the nonsynonymous ([d.sub.N]) and synonymous ([d.sub.s]) substitution rates is an important means of understanding mechanisms of molecular evolution. A [d.sub.N]/[d.sub.s] ratio >1 is taken as evidence of positive selection, whereas a [d.sub.N]/[d.sub.s] ratio <1 is taken as evidence of purifying selection (25). The Codeml program of the PAML package (available from http://envgen.nox.ac.uk/bioinformatics/docs/codeml.html) (25) was used to estimate [d.sub.N][d.sub.s] ratios in the phylogenetic analysis shown in Figure 1. The values were calculated by using model 1 in the program, and default parameters were used for the execution of the program.

Phylogenetic Reconstruction

Phylogenies were constructed by the Bayesian method, as implemented by the program MrBayes (26) (available at no charge from www.mrbayes.net). The evolutionary model used was the General Time Reversible model (27). Because evolutionary rates are not homogeneous for every site in a gene, among-site variation in evolutionary rate was estimated separately for first, second, and third positions of sites within codons. Four chains, with a "temperature" of 0.2 for the heated chains, were run for each tree. Trees were sampled every 100 generations. A total of 10 million generations were run with a buru-in of 5,000 trees. The length of each buru-in was set at a value that exceeded twice the number of trees required for convergence upon a stable likelihood value. Because the consensus trees calculated by MrBayes do not include the posterior probabilities of the clades, each entire set of trees was imported into PAUP* (28) and the same trees used by MrBayes to calculate a consensus were used to calculate a 50% majority rule consensus in PAUP* (28). The resulting tree shows the posterior probabilities of the clades, i.e., the percentage of time that those taxa are included in the clade. The consensus trees calculated by MrBayes were imported into PAUP* for the purposes of displaying and printing the tree.

Results

Ancient Horizontal Transfer of [bla.sub.CTX-M] Ancestor

The NCBI genomes database (www.ncbi.nhn.nih. gov:80/entrez/query.fcgi?db=genome) currently contains the completed genomic sequences of 139 eubacterial organisms. Because Kluyvera are members of the Enterobacteriaceae group of the gamma subdivision of Proteobacteria, the genomes of other members of the gamma subdivision, and especially the chromosomes of Enterobacteriaceae, were of the greatest interest. BLAST searches of these microbial genomes show that among the complete genomic sequences available for 9 species of Enterobacteriaceae (Escherichia coli, Salmonella typhimrium LT2, S. enterica, Shigella flexneri, Photorhabdus luminescens, Buchnera aphidicola, Candidatus Blochmannia floridanus, Wiggles-worthia glossinidia, and Yersinia pestis) none contain chromosomal homologs of the [bla.sub.CTX-M]s that are detectable through sequence comparison. BLAST searches similarly show that many non-Enterobacteriaceae members of the gamma subdivision of Proteobacteria, Shewanella oneidensis, Haemophilus ducreyi, H. influenzae, Pseudomonas aeruginosa, P. putida, P. syringae, Vibrio cholerae, V. parahemolyticus, V. vulnificus, Xanthomonas axonopodis, X. campestris, and Xylella fastidiosa, also do not contain chromosomal homologs of the [bla.sub.CTX-M]s. However, BLAST searches of the translated nonredundant nucleotide database revealed that Kluyvera species, Citrobacter sedlakii, and Klebsiella oxytoca contain close chromosomal homologs of the [bla.sub.CTX-M]s .

[FIGURE 1 OMITTED]

These results show that the [bla.sub.CTX-M] homologs originally came into the chromosomes of K. oxytoca, Kluyverra species, and C. sedlakii by horizontal transfer because most species of gamma Proteobacteria do not contain [bla.sub.CTX-M] homologs. If [bla.sub.CTX-M] holnologs were vertically transmitted into the species that contain them, numerous deletions would be required to explain absence of those homologs in the majority of gamma Proteobacteria. However, only 3 horizontal insertions are required to explain the presence of [bla.sub.CTX-M] homologs in the chromosomes of K. oxytoca, Citrobacter species. Because fewer insertions than deletions are required to explain these data, insertion of [bla.sub.CTX-M] homologs into the chromosomes of those bacteria that contain them is the most likely explanation for their current distribution.

Divergence of the [bla.sub.CTX-M]s

The GenBank DNA and protein accession numbers of the sequences included in this analysis are shown in the online Appendix Table (available from www.cdc.gov/EID/ content/14/3/423-appT.htm), along with the organism in which the gene exists and whether the gene was located on a plasmid or a chromosome. The results of our phylogenetic analysis are presented in Figure 1. The groupings of [bla.sub.CTX-M]s on our phylogenetic analysis agree with the dendrogram published by Bonnet in a recent review (16); for purposes of clarity, we will use the same group names used in that review, as shown in Figure 1.

The phylogenetic analysis shows that the [bla.sub.CTX-M]s represent a fairly divergent group of [beta]-lactamase genes descended from a common ancestor. The genes encoding the CTX-M-1 and CTX-M-2 groups are separated by over 400 mutations, which indicates considerable diversification of the [bla.sub.CTX-M]s. The average distance separating the [bla.sub.CTX-M]s from their most recent common ancestor is 226.2-nt [+ or -] 22.8-nt mutations, which indicates that the rates of evolution among the [bla.sub.CTX-M]s have been similar.

Positive selection testing within the phylogenetic analysis shows that positive selection has occurred throughout the evolutionary history of the class A [beta]-lactamases. More positive selection appears to exist at branches deep within the tree than along more recent branches. The branches during which most of the divergence of the [bla.sub.CTX-M]s occurred are characterized by purifying selection. The detection of purifying selection suggests a slow evolutionary rate and that the [bla.sub.CTX-M]s diverged in ancient times. More recent evolution of the [bla.sub.CTX-M]s likely can be characterized by intense positive selection, but the branches at the tips are still too short to obtain reliable [d.sub.N]/[d.sub.s] ratios.

Mobilization of [bla.sub.CTX-M]s to Plasmids

The [bla.sub.CTX-M]s have been mobilized from the chromosomes of various Kluyvera species to plasmids at least 8 times since they diverged from their most recent common ancestor as indicated in Figure 1. The alleles in the CTXM-2 group have been mobilized from the chromosome of Kluyvera ascorbata at least twice (29). The alleles in the CTX-M-9 group have been mobilized once from the chromosome of Kluyvera georgiana (30). The alleles from the CTX-M-8 group were mobilized once from the chromosome of K. georgiana (20). The CTX-M-25 group has been mobilized once, although the species from which it originates has not yet been determined. The alleles in the CTX-M-1 group have been mobilized at least 3 times (17,18,31), and one of those mobilizations has been from the chromosome of K. ascorbata. When compared with the [bla.sub.CTX-M]s, which have been mobilized once, and the [bla.sub.SHV]s , which have been reported to have been mobilized 2-4 times (32), the number of mobilization events that have occurred among the [bla.sub.CTX-M]s is atypically high.

To compare the number of mobilizations that have occurred in the CTX-M group with those that have occurred in the rest of the class A [beta]-lactamases, we constructed a phylogenetic analysis of class A alleles that spans the breadth of this group and that contains representatives of all groups of class A alleles known to the authors (Figure 2). Among all of the class A [beta]-lactamases, including the [bla.sub.CTX-M]s , only 22 mobilizations to plasmids were found. To quantitatively compare the numbers of times that CTX-Ms have been mobilized to plasmids with the number of times that other class A [beta]-lactamases have been mobilized to plasmids, the total number of mutations that have occurred within the [bla.sub.CTX-M] clade were summed and divided by the number of mobilizations that have occurred in that region of the phylogenetic analysis. Among the [bla.sub.CTX-M]s, the ratio of mobilizations to mutations is 1 mobilization per 191 mutations. Among the remainder of the tree when the [bla.sub.CTX-M] clade is excluded from the analysis, 14 mobilizations occur with the ratio of mobilizations to mutations being 1 mobilization per 2,471 mutations. When the complete phylogenetic analysis is considered, 1 mobilization occurs per 1,870 mutations. By that comparison, the mobilization of the [bla.sub.CTX-M] genes to plasmids has occurred at an unusually high rate. This result is unlikely to be an artifact of sampling bias or clinical interest because other class A [beta]-lactamases have been intently studied for a longer period than the [bla.sub.CTX-M]s. If any bias exists in the data, it would be the undersampling of [bla.sub.CTX-M] mobilizations relative to other class A [beta]-lactamases.

[FIGURE 2 OMITTED]

Because nearly one half of the mobilizations that have occurred in the class A phylogenetic analysis have occurred among the [bla.sub.CTX-M]s, it seemed reasonable to conclude that the circumstances associated with the mobilizations of the [bla.sub.CTX-M]s may differ from the circumstances associated with the mobilizations of other class A [beta]-lactamases. To rule out any effect that varying intensities of selection or varying evolutionary rates might have on mobilizations to plasmids, we divided the phylogenetic analysis into several monophyletic groups for further analysis. Within the class A [beta]-lactamase phylogenetic analysis (Figure 2) are several monophyletic clades that descended from a single ancestor (node A). Each of the monophyletic clades that descended from node A were considered separately during positive selection testing except for the monophyletic clade that contains the [bla.sub.CTX-M]s; it was divided into 2 separate clades so that a clade containing only the [bla.sub.CTX-M]s and their closest relatives could be considered. Monophyletic clades that diverged before the point represented at node A were also examined individually.

The [d.sub.N]/[d.sub.s] ratios were computed for each clade (Table), and a correlation coefficient for mobilizations and [d.sub.N]/[d.sub.s] ratios of 0.21 (p = 0.40) was calculated. The average distance of each clade from the root of the tree was also computed (Table), and the correlation coefficient for mobilizations and average distance from the root is 0.21 (p = 0.41). The nonsignificant p values yielded by those results mean that the unusually high number of mobilizations among the [bla.sub.CTX-M]s are probably not an artifact caused by positive selection or evolutionary rate.

Most of the mobilizations of the [bla.sub.CTX-M]s have occurred in recent years because genes that are identical ([bla.sub.CTX-M3a] [18] and [bla.sub.CTX-M-18] [19]) or nearly identical to the ancestors of plasmidic clades (Figure 1) have been found in the chromosomes of Kluyvera species, whereas many of the other plasmidic class A [beta]-lactamases have been mobilized much longer, perhaps even since ancient times. In many cases, no chromosomal ancestor is identified and the plasmidic resistance genes are not closely related to the chromosomal resistance genes of any identified groups of bacteria.

Discussion

Although the use of antimicrobial agents generally has enhanced the spread of antimicrobial drug resistance among bacteria by providing the selective pressure needed for the emergence of novel resistance determinants, selective pressure alone does not explain the increasing frequency with which [bla.sub.CTX-M] alleles have been noted in bacterial populations in recent years (16,33). Although [bla.sub.CTX-M] alleles tend to be located on transmissible plasmids and transposable elements, which clearly facilitate their dissemination, the repeated mobilization of the [bla.sub.CTX-M]s from the chromosomes principally among Kluyera species is most intriguing. The mechanistic basis underlying this repeated mobilization to plasmids remains unknown. Whether the chromosomes of Kluyvera species have some unique aspect that enhances the mobilization of the [bla.sub.CTX-M] genes remains to be determined. Other factors, such as exposure of the isolates to specific antimicrobial agents or to environmental changes that facilitate the mobilization of [bla.sub.CTX-M]s to plasmids also need to be investigated.

Two insertion elements are known to contribute to the mobilization of [bla.sub.CTX-M]s. The first, which is associated with the CTX-M-2 and CTX-M-9 groups, is ISCR1 (34) and the second, which is associated with the CTX-M-1, CTX-M-8, and CTX-M-25 groups is ISEcp1 (35). According to our phylogenetic analysis (Figure 1), 4 mobilizations can be attributed to each of these insertion type elements. Thus, both elements seem to promote equal frequencies of mobilization of [bla.sub.CTX-M]s. Notably, ISCR1 was also reported to be responsible for the mobilization of both [bla.sub.VEB] and [bla.sub.PER] alleles, but neither of these resistance determinants has been reported to have an unusually high rate of mobilizations from chromosomal locations to plasmids.

Another factor that may contribute to the rate of mobilizations of the [bla.sub.CTX-M] resistance determinants is the frequency of plasmids in bacterial populations. As the number of plasmids increases in microbial populations, so does the number of target replicons. A comparison of the percentage of bacterial strains that contained plasmids in the pre--anti-microbial drug era (36) with the percentage of contemporary strains that carry plasmids (37,38) indicates that the frequency of plasmid carriage has increased from 19% to 58%-100%, depending on the species surveyed. Although the collection methods and resistance detection assays varied in the studies used for this comparison (which may have introduced biases toward an increasing frequency of plasmids), few doubt that plasmid carriage is much more common among bacterial strains in the antimicrobial drug era (39,40). Unfortunately, specific information about plasmid carriage of Kluyvera species versus other Enterobacteriaceae is not available.

Regardless of the mechanism, the increased number of mobilizations of [bla.sub.CTX-M]s from their chromosomal locations among relatively rare human pathogens to plasmids that circulate widely among several important human and animal pathogens (particularly among E. coli) is a serious public health concern. The results of our study indicate the potential for an increase in the rate of mobilization of a variety of other resistance determinants to plasmids. Such an increase could result in more rapid mobilizations of novel resistance determinants and contribute to the accelerated spread of antimicrobial resistance determinants among a large spectrum of bacterial pathogens.

Use of trade names is for identification only and does not imply endorsement by the Public Health Service or by the U.S. Department of Health and Human Services.

This study was supported as part of Phase V of Project ICARE through unrestricted research grants to the Rollins School of Public Health of Emory University from AstraZeneca, bioMerieux, Elan, J&J PRD, Pfizer, and 3M Health Care; and by start-up support from University of California, Merced.

References

(1.) Lawrence JG. Clustering of antibiotic resistance genes: beyond the selfish operon. American Society for Microbiology. ASM News. 2000;66:281-6 [cited 26 Jan 2008]. Available from http://newsarchive.asm.org/may00/feature2.asp

(2.) Kim S, Kim J, Kang Y, Park Y, Lee B. Occurrence of extended-spectrum beta-lactamases in members of the genus Shigella in the Republic of Korea. J Clin Microbiol. 2004;42:5264-9.

(3.) Baraniak A, Fiett J, Mrowka A, Walory J, Hryniewicz W, Gniadkowski M. Evolution of TEM-type extended-spectrum beta-lactamases in clinical Enterobacteriaceae strains in Poland. Antimicrob Agents Chemother. 2005;49:1872-80.

(4.) Brinas L, Lantero M, de Diego I, Alvarez M, Zarazaga M, Tortes C. Mechanisms of resistance to expanded-spectrum cephalosporins in Escherichia coli isolates recovered in a Spanish hospital. J Antimicrob Chemother. 2005;56:1107-10.

(5.) Ho PL, Ho AY, Chow KH, Wong RC, Duan RS, Ho WL, et al. Occurrence and molecular analysis of extended-spectrum beta-lactamase-producing Proteus mirabilis in Hong Kong, 1999-2002. J Antimicrob Chemother. 2005;55:840-5.

(6.) Ho PL, Shek RH, Chow KH, Duan RS, Mak GC, Lai EL, et al. Detection and characterization of extended-spectrum beta-lactamases among bloodstream isolates of Enterobacter spp. in Hong Kong, 2000-2002. J Antimicrob Chemother. 2005;55:326-32.

(7.) Kruger T, Szabo D, Keddy KH, Deeley K, Marsh JW, Hujer AM, et al. Infections with nontyphoidal Salmonella species producing TEM-63 or a novel TEM enzyme, TEM-131, in South Africa. Antimicrob Agents Chemother. 2004;48:4263-70.

(8.) Schlesinger J, Navon-Venezia S, Chmelnitsky I, Hammer-Munz O, Leavitt A, Gold HS, et al. Extended-spectrum beta-lactamases among Enterobacter isolates obtained in Tel Aviv, Israel. Antimicrob Agents Chemother. 2005;49:1150-6.

(9.) Tasli H, Bahar IH. Molecular characterization of TEM- and SHV-derived extended-spectrum beta-lactamases in hospital-based Enterobacteriaceae in Turkey. Jpn J Infect Dis. 2005;58:162-7.

(10.) Medeiros AA. Evolution and dissemination of beta-lactamases accelerated by generations of beta-lactam antibiotics. Clin Infect Dis. 1997;24(Suppl 1):S19-45.

(11.) Datta N, Kontomichalou P. Penicillinase synthesis controlled by infectious R factors in Enterobacteriacae. Nature. 1965;208:239-41.

(12.) Barlow M, Hall BG. Predicting evolutionary potential: in vitro evolution accurately reproduces natural evolution of the TEM beta-lactamase. Genetics. 2002; 160:823-32.

(13.) Hall BG, Barlow M. Evolution of the serine beta-lactamases: past, present and future. Drug Resist Updat. 2004;7:111-23.

(14.) Roupas A, Pitton JS. R factor-mediated and chromosomal resistance to ampicillin in Escherichia coli. Antimicrob Agents Chemother. 1974;5:186-91.

(15.) Haeggman S, Lofdahl S, Paauw A, Verhoef J, Brisse S. Diversity and evolution of the class A chromosomal beta-lactamase gene in Klebsiella pneumoniae. Antimicrob Agents Chemother. 2004;48:2400-8.

(16.) Bonnet R. Growing group of extended-spectrum beta-lactamases: the CTX-M enzymes. Antimicrob Agents Chemother. 2004;48: 1-14.

(17.) Bauernfeind A, Grimm H, Schweighart S. A new plasmidic cefotaximase in a clinical isolate of Escherichia coli. Infection. 1990;18:294-8.

(18.) Rodriguez MM, Power P, Radice M, Vay C, Famiglietti A, Galleni M, et al. Chromosome-encoded CTX-M-3 from Kluyvera ascorbata: a possible origin of plasmid-borne CTX-M-1-derived cefotaximases. Antimicrob Agents Chemother. 2004;48:4895-7.

(19.) Boyd DA, Olson AB, Silverman M, McGeer A, Willey BM, Pong-Porter V, et al. Identification of a progenitor of the CTX-M-9 group of extended spectrum beta-lactamases from Kluyvera spp. isolated in Guyana. In: 44th Interscience Conference on Antimicrobial Agents and Chemotherapy. 2004 Oct 30-Nov 2. Washington: American Society for Microbiology; 2004.

(20.) Poirel L, Kampfer P, Nordmann P. Chromosome-encoded Ambler class A beta-lactamase of Kluyvera georgiana, a probable progenitor of a subgroup of CTX-M extended-spectrum beta-lactamases. Antimicrob Agents Chemother. 2002;46:4038-40.

(21.) Humeniuk C, Arlet G, Gautier V, Grimont P, Labia R, Philippon A. Beta-lactamases of Kluyvera ascorbata, probable progenitors of some plasmid-encoded CTX-M types. Antimicrob Agents Chemother. 2002;46:3045-9.

(22.) Altschul SK Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403-10.

(23.) Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389-402.

(24.) Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The ClustalX window's interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997;25:4876-82.

(25.) Yang Z. Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol Biol Evol. 1998;15:568-73.

(26.) Huelsenbeck JP, Ronquist F. MrBayes: Bayesian inference of phylogenetic analysis. Bioinformatics. 2001;17:754-5.

(27.) Tavare L. Some probabilistic and statistical problems on the analysis of DNA sequences. Lec Math Life Sci. 1986;17:57-86.

(28.) Swofford DL. PAUP*. Phylogcnetic analysis using parsimony (*and other methods). Sunderland (MA): Sinauer Associates; 2000.

(29.) Bauernfeind A, Casellas JM, Goldberg M, Holley M, Jungwirth R, Mangold P, et al. A new plasmidic cefotaximase from patients infected with Salmonella typhimurium. Infection. 1992;20:158-63.

(30.) Olson AB, Silverman M, Boyd DA, McGeer A, Willey BM, Pong-Porter V, et al. Identification of a progenitor of the CTX-M-9 group of extended-spectrum beta-lactamases from Kluyvera georgiana isolated in Guyana. Antimicrob Agents Chemother. 2005;49:2112-5.

(31.) Baraniak A, Sadowy E, Hryniewicz W, Gniadkowski M. Two different extended-spectrum beta-lactamases (ESBLs) in one of the first ESBL-producing salmonella isolates in Poland. J Clin Microbiol. 2002;40:1095-7.

(32.) Ford PJ, Avison MB. Evolutionary mapping of the SHV beta-lactamase and evidence for two separate IS26-dependent blaSHV mobilization events from the Klebsiella pneumoniae chromosome. J Antimicrob Chemother. 2004;54:69-75.

(33.) Rodriguez-Bano J, Navarro MD, Romero L, Muniain MA, de Cueto M, Rios MJ, et al. Bacteremia due to extended-spectrum beta-lactamase-producing Escherichia coli in the CTX-M era: a new clinical challenge. Clin Infect Dis. 2006;43:1407-14.

(34.) Yoleman MA, Bennett PM, Walsh TR. ISCR elements: novel gene-capturing systems of the 21st century? Microbiol Mol Biol Rev. 2006;70:296-316.

(35.) Poirel L, Lartigue MF, Decousser JW, Nordmann P. ISEcp1B-mediated transposition of blaCTX-M in Escherichia coli. Antimicrob Agents Chemother. 2005;49:447-50.

(36.) Datta N, Hughes VM. Plasmids of the same lnc groups in Enterobacteria before and after the medical use of antibiotics. Nature. 1983;306:616-7.

(37.) Souza V, Rocha M, Valera A, Eguiarte LE. Genetic structure of natural populations of Escherichia coli in wild hosts on different continents. Appl Environ Microbiol. 1999;65:3373-85.

(38.) Wan KF, Radu S, Cheah YK, Benjamin PG, Ling CM, Hon SF, et al. Antibiotic resistance, plasmid profile and RAPD-PCR analysis of enteropathogenic Escherichia coli (EPEC) clinical isolates. Southeast Asian J Trop Med Public Health. 2003;34:620-6.

(39.) Rainey PB, Cooper TF. Evolution of bacterial diversity and the origins of modularity. Res Microbiol. 2004;155:370-5.

(40.) Mirkin BG, Fenner TI, Galperin MY, Koonin EV. Algorithms for computing parsimonious evolutionary scenarios for genome evolution, the last universal common ancestor and dominance of horizontal gene transfer in the evolution of prokaryotes. BMC Evol Biol. 2003;3:2.

Address for correspondence: Miriam Barlow, University of California, Merced, PO Box 2039, Merced, CA 95344, USA; email: mbarlow@ ucmerced.edu

Miriam Barlow, * Rebecca A. Reik, ([dagger]) Stephen D. Jacobs, * Monica Medina, * Matthew P. Meyer, * John E. McGowan, Jr., ([dagger]) and Fred C. Tenover ([double dagger])

* University of California, Merced, California, USA; ([dagger]) Emory University, Atlanta, Georgia, USA; and ([double dagger]) Centers for Disease Control and Prevention, Atlanta, Georgia, USA

Dr Barlow is a founding faculty member at the University of California, Merced. Her research focuses on the evolution of plasmidic resistance determinants, with particular emphasis on [beta]-lactamases.

Table. Average distances from root and
[d.sub.N]/[d.sub.s] ratios of monophyletic clades

         Average     [d.sub.N]/
         distance     d.sub.s
Clade   from root     ratio *     Mobilizations

1        3140.865      0.0935           5
2        2534.35       0.4535           0
3        2538.90       0.2446           0
4       2819.37634     0.0929           9
5        3139.50       0.0426           0
6        2575.30       0.5024           0
7        2975.70       0.1019           1
8        2443.35       0.2618           0
9        2532.075      0.1969           1
10       3267.81       0.0338           0
11       2429.70       0.0118           1
12       2811.90       0.016            0
13       2218.125      1.1595           0
14       1747.20       0.9331           1
15        873.60       0.2056           0
16        518.70       0.1069           1
17        559.65       0.1595           1
18        778.05       0.2496           0

* [d.sub.N], nonsynonymous substitution rate,
[d.sub.s], synonoymous substitution rate.