Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, February 2009, p. 716-727, Vol. 53, No. 2
0066-4804/09/$08.00+0 doi:10.1128/AAC.00645-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.
,
Zulfiqar A. Bhutta,5
Rumina Hasan,5
Gordon Dougan,1 and
John Wain1*
The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, United Kingdom,1 Centre for Molecular Microbiology and Infection, Imperial College, London, United Kingdom,2 Bacterial Microarray Group, Department of Medical Microbiology, St Georges Hospital Medical School, London, United Kingdom,3 Defence Science and Technology Laboratory, Porton Down, Salisbury, Wiltshire, United Kingdom,4 Aga Khan University Medical Centre, Karachi, Pakistan5
Received 16 May 2008/ Returned for modification 22 July 2008/ Accepted 13 October 2008
|
|
|---|
|
|
|---|
Since the early 1970s antibiotic resistance has become established in S. enterica serovar Typhi, the cause of human typhoid (40). Two years after chloramphenicol was first used for treating typhoid fever in 1948 (42), the first resistant serovar Typhi isolate was described (4). However, it took another 22 years before chloramphenicol resistance became a major clinical problem, with the first major outbreak of resistant typhoid fever in 1972 (21). Thereafter, the global spread of chloramphenicol-resistant serovar Typhi occurred very rapidly, and over the next two decades further antibiotic resistance genes were acquired to generate MDR serovar Typhi (40). Although occasionally IncF, IncP, and IncB/O plasmids are found (20), it is striking that the majority of resistance plasmids from serovar Typhi belong to the HI1 incompatibility group (IncHI1). Thus, IncHI1 plasmids appear to have evolved a stable relationship with the serovar Typhi bacterial host so that even in the absence of obvious antibiotic selection, a proportion of the isolates maintain the plasmid (7).
There are three fully sequenced IncHI1 plasmids publicly available: R27 (31), pHCM1 (24), and pAKU1 (13). Holt et al. (13) have reported a detailed comparison of these plasmids showing a 164.4-kb shared IncHI1-associated backbone and multiple antibiotic resistance genes carried within Tn10 and a Tn9/Tn21 composite transposon.
A study of IncHI1 plasmids isolated from outbreaks and sporadic typhoid cases in Vietnam between 1993 and 1996 revealed seven different HindIII-digested restriction fragment length polymorphism (RFLP) patterns (patterns 1 to 7), with pattern 1 being the most common (38). The trend was reversed after 1996, at which point pattern 7 became predominant, while plasmids of other patterns were no longer detected. Based on this study, it was hypothesized that a milestone of IncHI1 plasmid evolution had been captured, during which a newly acquired increase in fitness drove one plasmid type to outcompete all others circulating in the serovar Typhi population in Vietnam.
We extend our analyses here to a global collection of IncH1 plasmids isolated from serovar Typhi over the past several decades using a combination of DNA microarrays, PCR, and sequence analysis. We have devised a multilocus sequence typing (MLST) scheme (17) for the IncHI1 plasmids, which is called plasmid MLST (PMLST).
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Plasmids used for microarray analysis in this study
|
Microarray experiments. Pairs of specific oligonucleotide primers were designed for each predicted coding sequence (CDS) based on the EMBL gene sequence files of pHCM1 (AL513383) and pHCM2 (AL513384) using Microarray Design Software (K. Vass, Beatson Institute for Cancer Research, Glasgow, United Kingdom) and the primer design program Primer3 (27). Whole genomic DNA extracted from serovar Typhi CT18 was used as a template in the PCRs. PCR amplification was carried out using standard conditions in a RoboAmp4200 (MWG-Biotech, Germany). The microarray was printed by using MicroGrid II Biorobotics (MWG-Biotech) onto poly-L-lysine-coated slides, with four copies of each CDS per array. Specific PCR products corresponding to six genes located on the serovar Typhi chromosome and Cy3-dCTP- and Cy5-dCTP-labeled nucleotides were printed as controls.
Whole genomic serovar Typhi DNA was used as a template for direct incorporation of fluorescent labeled nucleotides. In each experiment the control DNA (serovar Typhi CT18 whole genomic DNA) was labeled with Cy5-dCTP (Amersham, United Kingdom) and the test DNA with Cy3-dCTP (Amersham). More than 5 µg of whole genomic DNA of control and test isolates and 3 µg of random primers (Roche, Welwyn Garden City, United Kingdom) in 41.5 µl of water were heated at 95°C for 5 min, snap cooled on ice, and centrifuged briefly. Each labeling reaction mix contained 100 µM dA/T/GTP, 40 µM dCTP, 300 µM Cy3-dCTP, or Cy5-dCTP, and 5 U of Klenow reagent in a final volume of 50 µl. Reaction mixes were incubated at 37°C in the dark for 90 min. Microarray slides were incubated in a prehybridization solution (3.5x SSC [1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate], 0.1% sodium dodecyl sulfate [SDS], and 10 mg of bovine serum albumin/ml in a final volume of 50 ml of H2O) at 65°C for 20 min. After prehybridization slides were rinsed thoroughly in H2O for 1 min and propan-2-ol for 1 min and then dried.
Labeled DNA was purified and eluted in a final volume of 21 µl, with a final concentration of 4x SSC and 0.3% SDS. The hybridization mixture was denatured at 95°C for 2 min, allowed to cool, applied to the microarray, and covered with a coverslip. Hybridization was carried out at 65°C for 16 to 20 h in a waterproof hybridization chamber. After hybridization, slides were washed at 65°C in 1x SSC and 0.06% SDS for 2 min, followed by two washes in 0.06x SSC each for 2 min. To ensure reproducibility of data, three separate hybridizations were carried out for each plasmid isolate.
Data input and analysis.
Scanning of microarrays was carried out by using an Affymetrix 428 dual laser array scanner. Scanned data were analyzed and quantified in ImaGene v4.2 (Biodiscovery, Inc.). The quantified data file was extracted and normalized using GeneSpring version 5.1 (Silicon Genetics). Genes were classed as present or absent/highly divergent as follows. Briefly, a list of "core" genes was derived. To allow for divergence in primary DNA sequence or poor hybridization signal, this was defined as genes present in at least 19 of 20 serovar Typhi plasmids hybridized to the microarray. Insertion sequences, transposons, and antibiotic resistance genes were excluded from this analysis since the plasmids were selected because of the presence of resistance and so these genes were necessarily present even though they are not part of the plasmid's survival functions. For each isolate the median of the normalized ratio (baseline) and standard deviation were calculated for the 141 "core" genes (26). Three standard deviations (P
0.01) above and below the baseline for each isolate were used to set the threshold for those genes defined as "present." Genes with normalized ratios higher than three standard deviations above the baseline were considered duplicated, whereas ratios that were less than three standard deviations below the baseline were considered absent or highly divergent.
PCR for genetic elements. Oligonucleotide primers were taken from the literature or designed (Primer3 v.0.4.0 [http://frodo.wi.mit.edu/]) to specifically amplify the elements listed in Table 2. For confirmation of microarray data, long range PCRs were designed (see Table S2 in the supplemental material). Amplification of serovar Typhi whole genomic DNA was carried out using standard PCR conditions, the Tm for each primer is given in Table 2 or in Table S2 in the supplemental material. Specific amplicons were identified by separation on 1.0% agarose gel, followed by staining with ethidium bromide and visualization under UV light.
|
View this table: [in a new window] |
TABLE 2. Oligonucleotide primers
|
|
View this table: [in a new window] |
TABLE 3. Candidate PMLST genes and their annotations
|
PMLST allele accession numbers. We have deposited the DNA sequences of PMLST alleles in the NCBI database under the following accession numbers: HCM1.043 allele 1, FJ183728; HCM1.043 allele 3, FJ183729; HCM1.064 allele 1, FJ183730; HCM1.064 allele 2, FJ183731; HCM1.099 allele 1, FJ183732; HCM1.099 allele 2, FJ183733; HCM1.099 allele 3, FJ183734; HCM1.116 allele 1, FJ183735; HCM1.116 allele 2, FJ183736; HCM1.178ac allele 01, FJ183737; HCM1.178ac allele 02, FJ183738; HCM1.259 allele 01, FJ183739; HCM1.259 allele 02, FJ183740; and HCM1.259 allele 03, FJ183741.
|
|
|---|
![]() View larger version (41K): [in a new window] |
FIG. 1. Variation in IncHI1 plasmids from serovar Typhi determined by microarray. (a) CDSs are defined in the sequence of pHCM1 in the outer ring (NC_003384). The inner concentric rings each represent the CDSs predicted as present in each different plasmid by microarray. The order, from the inside is as follows: blue, R27 (United Kingdom, 1961); magenta, 40R181 (Mexico, 1972); pink, 40R344 (India, 1972); light green, 40R311 (Thailand, 1970s); dark green, 42R917 (Vietnam, 1970s); yellow, pSTY1 (Vietnam, 1993); sand, pSTY4 (Vietnam, 1993); gray, pSTY6 (Vietnam, 1993); black, pSTY7 (Vietnam, 1996); and red and blue, pHCM1 (Vietnam, 1993) on the outside. The CDSs of pHCM1 are marked blue when present and red when absent from R27 sequence (NC_002305). (b) Details of regions A to E shown in panel a. These are the main regions of variability between plasmids. PMLST clusters are shown in Fig. 3.
|
Long-range overlapping PCR confirms microarray data. As a complement to the microarray data above, which demonstrate the presence but not the position of core genes in each plasmid, a set of oligonucleotide primers was designed to amplify overlapping fragments around the whole pHCM1 plasmid. Performing the same PCRs with other IncHI1 plasmids gives an indication of those regions similar to pHCM1, where PCR products are obtained, while suggesting which regions show variation, when PCR products are absent or a different size compared to that obtained for pHCM1.
Initially, 20 oligonucleotide pairs were designed to amplify overlapping 10-kb fragments covering the whole of pHCM1 (Fig. 2a). Thereafter, for regions in which variability was indicated, oligonucleotides for overlapping 5-kb fragments were obtained (Fig. 2b). Reactions for oligonucleotide pairs from 2 to 8 and 18 were positive in all plasmids that have been analyzed by 10-kb PCR. These PCRs cover CDSs encoding plasmid maintenance and transfer functions in pHCM1. Three adjacent reactions (1, 19, and 20 in Fig. 2a) were negative in pSTY6; further analysis confirmed a deletion of 10 to 15 kb in this region (data not shown). The regions that did not amplify in the other IncHI1 plasmids are those encoding antibiotic resistance-determining regions in pHCM1. This suggests that the IncHI1 plasmids possess an almost invariable core of genes that code for plasmid maintenance and transfer, the order of which is generally maintained between the plasmids, and variable regions incorporating the antibiotic resistance determinants.
![]() View larger version (29K): [in a new window] |
FIG. 2. PCR confirmation of microarray results. (a) Long-range PCR of 10-kb amplicons around IncHI1 plasmids of serovar Typhi. The products for R27 (inner gray ring) are predicted from the DNA sequence, and the primer binding sites are at either end of the gray bars on the inner ring. The order from the inside is R27, pSTY1, pSTY4, pSTY6, and pST7, and pHCM1 is the outside ring. The outmost ring shows positions of six PMLST loci on pHCM1. (b) PCR analysis comparing the two most successful IncHI1 plasmids of serovar Typhi using 5-kb amplicons within regions that do not give products with 10Kb PCR primers. The thin gray lines are predicted products from the pHCM1 sequence. The order from the inside is 5-kb fragments of pSTY7, pSTY1, and pHCM1 and then 10-kb fragments of pSTY7, pSTY1, and pHCM1.
|
Plasmid pHCM1 contained two different DNA regions encoding mercury resistance, annotated as HCM1.151c to HCM1.159 (similar to mer operon from Tn1696) and HCM1.230 to HCM1.235 (similar to mer operon from Tn21), and neither of these regions was present on R27. Of the five plasmids analyzed from the 1970s, only the Mexican plasmid 40R181 harbored similar mercury resistance determinants. Of the later plasmids, p61-57-85 (isolated in Pakistan in 2000) and pSTY7 (isolated in Vietnam in 1996) gave microarray signatures suggesting the presence of a single mercury resistance determinant, whereas all plasmids from Vietnam isolated from 1993 to 1994 had both operons. Mercury resistance is therefore a phenotype acquired on more than one occasion by the IncHI1 plasmids. The efflux pump encoded by these operons is very specific for mercury and may play a role in protection against mercury in the human digestive tract, from dental filings (33), or in the environment.
Antibiotic resistance is a major phenotype conferred by IncHI1 plasmids in serovar Typhi. The resistance genes identified by microarray analysis consistently match the observed resistance phenotype. It is clear that these plasmids have acquired resistance genes sequentially into the variable regions, probably by the insertion of transposons into transposons (24). In some cases a resistance phenotype is always encoded by the same gene; cat (HCM1.209), tetRA (HCM1.240 and HCM1.241c), blaTEM-1 (HCM1.216), and strAB (HCM1.223 and HCM1.224) are conserved across all plasmids conferring resistance to chloramphenicol, tetracycline, ampicillin, and streptomycin, respectively. The genetic context of these genes in different plasmids, however, is not the same. From genome sequence data we know that tetracycline resistance in R27 is encoded within a complete Tn10 but that the IS10-right, the active transposase from Tn10, and part of the tetracycline resistance operon (tetCD), are not present in pHCM1. The array data are consistent with this and also suggest that plasmids 40R181 (Mexico, 1972) and 44R315 and 44R311 (Thailand, 1970s) possess only the IS10-left transposase (HCM1.194), while plasmids 40R344 (India, 1972) and 42R917 (Vietnam, 1970s) are absent or divergent for the IS10-left transposase (HCM1.194c and HCM1.249c). Sequence data also showed that Tn10 was inserted in different locations in R27 and pHCM1. These sets of evidence suggest that Tn10 has been independently acquired by IncHI1 plasmids on at least three occasions.
Different alleles conferring trimethoprim resistance were present within the IncHI1 plasmids. In pHCM1 the dihydrofolate reductase responsible for trimethoprim resistance is a dfrA14 gene cassette within a truncated version of a class one integron. Truncation has occurred at the 59-base element within a Tn21-like element (8). However, other plasmids isolated in Vietnam from 1992 to 1994 all had an intact class one integron containing the dfrA14 trimethoprim resistance gene cassette (1-kb fragment; Table 4). Later plasmids from Vietnam (pSTY7, pCB, and pD6) had a smaller gene cassette containing dfrA7 gene (0.75-kb fragment; Table 4). It is clear that class one integrons are involved in the acquisition of trimethoprim resistance by capturing different gene cassettes.
|
View this table: [in a new window] |
TABLE 4. Presence or absence of common integrase genes as determined by PCR and product size using primers 5'CS and 3'CS to amplify gene cassettes from class 1 integrons
|
PMLST, a typing scheme with phylogenetic implications. The presence of conserved regions among these plasmids provides material for a molecular typing scheme. Here we apply the principle of MLST (17) to this group of plasmids (PMLST). Thirteen loci not associated with mobile elements were sequenced and analyzed for fourteen of the plasmids to evaluate their usefulness in the PMLST scheme (Table 3). Four loci (HCM1.054, HCM1.094, HCM1.107, and HCM1.277) were excluded from the PMLST as they showed no variation. Three more loci (HCM1.177, HCM1.280c, and HCM1.286) showed inconsistent PCR amplification among the plasmids analyzed and were therefore excluded. The final PMLST set includes 6 loci: HCM1.043, HCM1.064, HCM1.099, HCM1.116, HCM1.178ac, and HCM1.259. Upon sequencing in 36 plasmids, two to three alleles (variants) of each locus were detected, the combinations of which define eight unique plasmid sequence types (STs) (Fig. 3 and see Table S3 in the supplemental material). Figure 3 shows the relatedness of 36 IncHI1 plasmids based on their STs, determined by using eBURST. The eBURST program grouped STs sharing five out of six identical loci, resulting in two groups and a singleton ST. Group 1 contains ST1, ST2, ST3, and ST4, while group 2 includes ST6, ST7, and ST8. Plasmid R27 was assigned ST5 and stood as a singlet, sharing no more than four identical loci with any ST in either group. ST1 includes pHCM1 and four other plasmids, all isolated from sporadic typhoid cases in Vietnam in 1993. Plasmids from ST2 were isolated earlier in the 1970s from Thailand (40R311 and 40R315) and India (40R344) except for pSTY4, which was isolated in Vietnam in 1993. ST3 and ST4 have only one member each, isolated in the 1970s from Mexico (40R181) and Vietnam (42R917), respectively. In group 2, ST6 contains plasmids from Vietnam and Jordan, all isolated after 1993. ST7 consists of 13 plasmids isolated from S. enterica serovar Paratyphi A from Karachi, Pakistan (2002 to 2004). The serovar Paratyphi A plasmid AKU1, which has recently been sequenced and analyzed (13), also belongs to this ST. ST8 contains four other S. Paratyphi A plasmids from Karachi (2003 to 2004) and differs from ST7 by the deletion of locus HCM1.116.
![]() View larger version (20K): [in a new window] |
FIG. 3. Modified eBURST diagram for PMLST. The numbers inside the circles represent the STs. Lines connecting circles represent single-locus variants. The dotted circles represent groups of related plasmids. Plasmids in group 1 are all from before 1993, whereas plasmids from group 2 are from 1993 to 2004. The singlet ST is the first described IncHI1 plasmid R27.
|
Differences within the conserved regions. Of the whole plasmid set analyzed, there are three fully sequenced plasmids, each of which represents one group in our PMLST map. The sequence comparison among these three plasmids (pHCM1, R27, and pAKU1) has been described elsewhere (13). We highlight here some features within the core region that we have defined by microarray and long-range PCR (Fig. 1 and 2) and further investigated on our full set of plasmids. Sequence comparison between pHCM1 and R27 defined five variable regions, two of them (regions A and B) encode antibiotic resistance determinants and mobile elements, while the other three (regions C, D, and E) were situated among the conserved regions defined by long-range PCR and do not show any evidence of mobile elements. The genes in regions C, D, and E were present in pHCM1 and absent in both R27 and pAKU1. The sequence contexts surrounding these three regions are similar in R27 and pAKU1. Another feature that lies within the conserved region is a composite transposon containing two IS1 flanking two CDSs (SPAP0105 and SPAP0106). This transposon (hereby called Ins1056) was inserted in the middle of the transfer region of pAKU1 (right after trhC). PCR assays were designed to screen for the presence or absence of these regions in the whole set of plasmids. To avoid false negatives, these assays were designed such that we expect a large PCR product when the region is present and a smaller PCR product when the region is absent.
Figure 4 shows the presence or absence of regions C, D, E, and Ins1056 in each plasmid, along with the phylogenetic tree of plasmids of all STs. Ins1056 was specific to plasmids of group 2, which suggests the singleton plasmid R27 may be more closely related to group 1 than to group 2. Regions C, D, and E were present in all five plasmids of ST1 (group 1). The remaining plasmids of group 1 were negative for regions C and E (except for 40R181). In group 2, regions C, D, and E were generally absent. Region D failed to amplify in the four plasmids of ST8 and two from ST7, suggesting further sequence changes had occurred that affected the binding sites of PCR primers.
![]() View larger version (40K): [in a new window] |
FIG. 4. IncHI1 plasmids phylogenetic tree and insertion/deletion (Ins/Del) events. The phylogenetic tree was built based on the single nucleotide polymorphisms identified within PMLST sequences using parsimony methods for discrete character data (PHYLIP software package), with bootstrap values shown from 1,000 resamplings. The insertion/deletion events were determined by PCR, with "a" for absent, "p" for present, and "-" for negative PCR.
|
It has been suggested previously that R27 and pAKU1 were more closely related to each other than either of them to pHCM1 (13). This is supported by our data showing the absence of regions C, D, and E in 24 plasmids from group 2 (including pAKU1), as well as their absence in R27. Combining the eBURST grouping with the phylogenetic tree and PCR results, we conclude that plasmids of group 2 (post-1993) did not evolve directly from group 1 (pre-1993) but belong to two different lineages.
Replication initiation. The complete sequence of plasmid pHCM1 revealed five genes that are predicted to encode proteins involved in the initiation of plasmid replication. Two, RepHI1B and RepHI1A, were classified as "core genes" in this analysis since they were present in all 20 serovar Typhi plasmids. However, RepE, RepA, and RepC were variable between isolates and so were classified as accessory genes (Table 5). These replication (Rep) proteins relate to incompatibility, and so this is a possible way in which different plasmid cores (10), or possibly different accessory gene linkage groups, could compete for replication space in the bacterial host. Expression of different Rep proteins, as seen with Neisseria gonorrhoeae (22), may explain the occasional reports of plasmids of different incompatibility groups in serovar Typhi (20). Three possible DNA-binding proteins (HCM1.58c, HCM1.175, and HCM1.286c) clustered within the core genes, as did the plasmid partition proteins parA (HCM1.86) and a putative parB encoding CDS (HCM1.87), with a series of 54-bp repeats upstream that could act as binding sites for the partition proteins.
|
View this table: [in a new window] |
TABLE 5. Variation in pHCM1 CDSs associated with plasmid replication
|
|
|
|---|
In an attempt to elucidate phylogenetic relationships between the IncHI1 plasmids of serovar Typhi, we sequenced six core genes from each of the IncHI1 plasmids, an approach based on MLST (17). All plasmids from the 1970s, and those circulating in serovar Typhi in Vietnam in the early 1990s, were grouped together (PMLST group 1). We therefore infer that they are directly related and represent the expansion of a successful plasmid backbone in the serovar Typhi population over a period of 20 years. A recognizably distinct plasmid core, PMLST group 2, was consistently present after 1993, suggesting successful competition with those already present in the serovar Typhi population.
Although the earlier RFLP study suggested a new plasmid type arose and spread in Vietnam in the mid 1990s, the origin of the plasmid could not be determined. In contrast to our initial hypothesis of stepwise evolution from R27 to pHCM1 to pSTY7, PMLST analysis suggests that pSTY7 (PMLST group 2) and pHCM1 (PMLST group 1) belong to distinct plasmid lineages. This is also supported by whole-plasmid sequence comparisons (13). Thus, the change in the predominant plasmid type in Vietnam is best explained by the acquisition and spread of a distinct plasmid type (group 2) rather than the clonal expansion of a particular variant already present in the serovar Typhi population. Comparison to the broader collection of plasmids suggests that this replacement of plasmid types in 1993 was not unique in Vietnam but was a global phenomenon.
The typing data for both serovar Typhi host and its plasmid are needed to investigate whether the predomination of group 2 plasmids is due to the spread of the strain carrying the plasmid or the spread of plasmid itself. There is limited data on the background strains of serovar Typhi circulating in Vietnam during these periods. Some authors have shown that multiple RFLP types of serovar Typhi had become MDR between 1993 and 1997 (5), while other suggest clonal expansion of single serovar Typhi strain between 1995 and 2002 (15). Unfortunately, several of the plasmids analyzed in the present study were provided as transconjugants in Escherichia coli, and so it was not possible to type the background serovar Typhi hosts.
We also investigated the accessory gene content of different PMLST groups. The location and date of isolation were the most important factors governing the accessory gene content (Fig. 1b), while date of isolation alone had the greatest influence on core gene PMLST. This suggests that IncHI1 plasmids circulate globally but pick up different accessory genes depending on the environment they inhabit.
These plasmids conjugate with a higher frequency (3 to 4 logs) at ambient temperatures (27°C and below) than at an in vivo temperature (37°C) (R. Curtiss III, unpublished data). This, combined with the fact that contaminated water is the main route for serovar Typhi transmission (2), suggests that the plasmids are most likely to spread in contaminated water at ambient temperature. Accessory genes may be exchanged in this environment.
RFLP is a useful method for plasmid typing that can be applied to most groups of plasmids. Several plasmids in the present study gave the same PMLST and had very similar accessory gene content but demonstrated different RFLP patterns (38). The variation detected by RFLP indicates that the diversity of IncHI1 plasmids in serovar Typhi decreased after 1994 (38). Attempts to draw trees to describe phylogenetic relationships between the different RFLP types, however, failed. This is because the variation detected by RFLP is generated by several genetic mechanisms with different effects: point mutation, rearrangement, restriction protection mechanisms, and insertion/deletion. For the present study PMLST, while only applicable to IncHI1 plasmids, was used to define plasmid groups based on the accumulation of mutations in nonselected regions of the plasmid core. The accumulation of mutations is a function of time, and this typing scheme therefore represents a quantitative estimate of variation which can be related to the ancestry of the plasmids.
Since microarray data gave us no information on gene order, we screened for rearrangements and then mapped the regions of variation using specifically designed PCR primers to amplify 5- to 10-kb fragments around the entire plasmid. The results of these long-range PCRs are presented in Fig. 3. The core genes were shown to share the same arrangement in all plasmids analyzed.
It is intriguing that IncHI1 plasmids, which can be easily transferred to E. coli at high frequency in vitro, are found almost exclusively in two S. enterica serovars, Typhi and Paratyphi A. Plasmid fitness cost must play a role in maintaining IncHI1 plasmids in Salmonella. This involves the integration of plasmid and chromosomal regulatory networks, such as H-NS (6). In serovar Typhimurium there is a detectable cost of having IncHI1 plasmid (6), and so it is possible that in E. coli there is also a cost. Advantages to the bacterial host harboring IncHI1 plasmids have also been demonstrated: an increased level of survival inside monocytic cell lines in serovar Typhimurium (6) or higher levels of serovar Typhi bacteremia during typhoid fever (39) or possibly enhanced tolerance to the conditions faced during the infection cycle.
Taken together, these analyses show that the IncHI1 plasmids of serovar Typhi contain truly conserved regions, as well as truly dynamic regions. This is reminiscent of the situation in E. coli K1 F-like plasmids (19). All of the IncHI1 plasmids studied shared a conserved backbone encoding a set of 163 core genes, into which several acquisitions of different antibiotic resistance and other accessory genes have occurred. These plasmids are thus sampling the mobile gene pool and picking up different alleles of the same resistance genes, which are presumably selected for during treatment. It can therefore be hypothesized that variation in the accessory gene content is driven by selection acting on a phenotype encoded by several distinct genotypes.
The adaptation of the serovar Typhi bacterial host to acquire and retain resistance plasmids, and the subsequent evolution of this bacterium-plasmid combination, has occurred over a 20- to 30-year time frame. It is likely that pressure from chloramphenicol treatment in the 1970s drove the selection of a linkage group of genes capable of survival in the serovar Typhi host. The presence of such linkage group was reported by our results to be the core gene set of IncHI1 plasmids. This core has evolved little, but point mutations accumulated in these conserved genes can be used to discern subtly different plasmid lineages, as demonstrated by our newly developed PMLST scheme. By PMLST we showed for the first time the population structure and the dynamic evolutionary process of IncHI1 plasmids. The two lineages defined by PMLST corresponded perfectly to time of plasmid isolation, suggesting that competition between lineages resulted in a global replacement of group 1 by group 2 plasmids.
Antibiotic pressure concurrently drove the evolution of IncHI1 plasmids by acquisition of accessory genes into this core at certain regions where disruption of DNA can occur without effecting core plasmid functions. Genes conferring the same resistance phenotype appear to have been acquired on multiple occasions. Genetic components that are not involved in resistance to toxic substances have also been acquired.
In summary, the selective pressure from antimicrobial chemotherapy has forced serovar Typhi to gain resistance plasmids that compete to colonize the bacterial host population by allowing this major pathogen access to a wide variety of accessory genes.
Further investigation is required in order to identify the molecular mechanisms responsible for the success of PMLST group 2 plasmids. The spread of these plasmids is also of interest since it may shed light upon the transfer route of antibiotic resistance, whether it is the synergistic spread of a bacterium-plasmid clone or the emergence of a more competitive plasmid. Interrogation of the association between serovar Typhi strains and IncHI1 plasmid types is vital to answer this question. The findings suggest a possible approach to destabilize IncHI1 plasmids in serovar Typhi, rendering them susceptible once again to first-line antibiotics.
We thank Henry Smith from Colindale, Christiane Dolecek and Chau Tran Thuy from Oxford Clinical Research Unit in the HCMC (Vietnam), To Song Diep, Microbiology, HTD, HCMC (Vietnam), and Bianca Paglietti from Università degli Studi di Sassari (Italy) for supplying IncHI1 plasmids.
Published ahead of print on 17 November 2008. ![]()
Supplemental material for this article may be found at http://aac.asm.org/. ![]()
Present address: School of Biosciences, University of Exeter, Devon EX4 4QD, United Kingdom. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»