The evolution of plasmid-carried antibiotic resistance

  • Fabian Svara1 and

    Affiliated with

    • Daniel J Rankin1, 2Email author

      Affiliated with

      BMC Evolutionary Biology201111:130

      DOI: 10.1186/1471-2148-11-130

      Received: 25 October 2010

      Accepted: 19 May 2011

      Published: 19 May 2011

      Abstract

      Background

      Antibiotic resistance represents a significant public health problem. When resistance genes are mobile, being carried on plasmids or phages, their spread can be greatly accelerated. Plasmids in particular have been implicated in the spread of antibiotic resistance genes. However, the selective pressures which favour plasmid-carried resistance genes have not been fully established. Here we address this issue with mathematical models of plasmid dynamics in response to different antibiotic treatment regimes.

      Results

      We show that transmission of plasmids is a key factor influencing plasmid-borne antibiotic resistance, but the dosage and interval between treatments is also important. Our results also hold when plasmids carrying the resistance gene are in competition with other plasmids that do not carry the resistance gene. By altering the interval between antibiotic treatments, and the dosage of antibiotic, we show that different treatment regimes can select for either plasmid-carried, or chromosome-carried, resistance.

      Conclusions

      Our research addresses the effect of environmental variation on the evolution of plasmid-carried antibiotic resistance.

      plasmid resistance antibiotic horizontal gene transfer mobile elements

      Background

      The emergence of antibiotic resistance in pathogenic bacteria, both in hospital and community-acquired infections, represents a significant public health problem [15]. Bacterial cells are capable of transferring genes horizontally. This DNA transfer can take place in three ways, through plasmids, phages, or uptake of naked DNA [6]. Plasmids are extra-chromosomal pieces of DNA, which are capable of replicating independently of the genome, and are particularly important in the spread of antibiotic resistance genes. Plasmids have been directly implicated in the acquisition of resistance to many antibiotics [714]. This is particularly problematic since plasmids can cross many species and genus barriers, and the rate of plasmid transfer has even been shown to increase in more heterogeneous communities [15]. Plasmids thus allow resistance to spread and persist in niches that are not necessarily subject to antibiotics [16].

      In Escherichia and Shigella strains, a larger proportion of the genome of plasmids codes for antibiotic resistance than that of the chromosome (Figure 1 - see figure legend for details). However, like all mobile genetic elements, plasmids face the fundamental trade-off between the cost of their mobility and the advantage imparted by their accelerated spread [e.g. [17]]. Plasmids can be seen as genomic parasites, as they can persist despite the costs they impose on their bacterial hosts [1820]. Previous models have examined the conditions under which plasmids can be maintained in the population [18, 19, 21, 22]. It is generally accepted that parasites can persist if their R 0, their per-capita rate of increase in a susceptible population of hosts, is greater than 1 [2325]. In the case of parasitic plasmids, which do not carry beneficial traits, and inflict a net cost on their host, this means that the rate of horizontal transfer must exceed the cost that they impose on their host [19]. Despite being fundamentally parasitic by nature, plasmids often carry traits which are beneficial to a host, or a host's neighbours [18, 20, 2631]. These dynamics will likely be modulated by environmental conditions, such as antibiotic treatment. Both the dosage of antibiotics and the interval between treatments is likely to affect the evolution of plasmid-borne resistance [e.g. [18, 21]]. Here we examine the role that these parameters will have on the evolution of plasmid-borne resistance. We address this question with a simple model of antibiotic resistance in a population exposed to different treatment regimes.
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Fig1_HTML.jpg
      Figure 1

      Percentage of the genome among 37EscherichiaandShigellastrains which code for genes related to antibiotic resistance. This shows that antibiotic resistance genes are over-represented on plasmids. Bars show the mean percentage (±s.e.) in all replicons (either on the chromosome, n = 37, or on plasmids, n = 53, including 38 plasmid genomes which contained no antibiotic resistance genes). The genomes are from 37 Escherichia and Shigella strains, including all chromosomal and plasmid genes in each genome (data obtained from the NCBI at http://​www.​ncbi.​nlm.​nih.​gov/​genomes/​). There were a total of 174,862 genes, with 170,709 on all 37 chromosomes and 4,153 on all 53 plasmids. For comparison with known genes involved in antibiotic resistance, the "Antibiotic Resistance Database" [http://​ardb.​cbcb.​umd.​edu/​- [62]] was used to identify genes involved in antibiotic resistance (homologues were identified using an E-value of 10-10 and percentage of sequence identity between 60-95%). From this, there were a total of 817 genes involved in antibiotic resistance, with 765 carried by the bacterial chromosomes (and 0.48% of all chromosomal genes) and 52 carried by plasmids (1.25% of all plasmid genes). A χ2 test revealed that antibiotic resistance genes were significantly over-represented on plasmids, compared to on the bacterial chromosome (Pearson's χ2 with Yates' continuity correction: χ2 = 54.6, p < 10-12, df = 1).

      Methods and Results

      Basic model

      We start off with a population of wild-type cells ("free" cells) that do not carry plasmids or antibiotic resistance genes. Assuming logistic growth, with a growth rate of r, environmental antibiotic concentration A, the dynamics of the wild type n F cells are:
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equ1_HTML.gif
      (1)

      The density-dependent death rate is given by a. This per capita death rate a depends on the density of cells in the population. In a population of wild-type cells this is simply n F , since there are no other cell types. We assume an antibiotic-induced fitness cost Am. Antibiotics can have two actions on a cell: they can either kill the bacteria (bactericidal antibiotics such as penicillin) or they can prevent reproduction (bacteriostatic antibiotics, such as tetracycline). Our model is general and does not differentiate between the two and could thus apply to either. From equation (1), at equilibrium, dn F /dt = 0, there are two equilibrium points, n F * = 0 and n F * = (r-mA)/a. This implies that the antibiotic will drive the host population extinct if mA > r. Thus, the higher the concentration of antibiotic, the more likely it is that wild-type will be eradicated.

      To keep our model simple, we simply focus on cells which are either resistant to the antibiotic, or not, and where the gene is localized in the genome (on a plasmid or a chromosome). We therefore assume that genes for resistance against an antibiotic can be carried either by the chromosome or by a plasmid. As such, our model does not make any other assumptions about the genetics of the systems, such as compensatory mutations or the number of mutations needed to confer resistance to the antibiotic. We assume that plasmids can transfer horizontally by conjugation at a rate β. The dynamics of wild-type cells now becomes:
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equ2_HTML.gif
      (2)
      Details of parameters used in all models are given in Table 1. The total population density is given by N = n F+n P+n C+n CP, where n P is the density of wild-type cells which are infected by a plasmid carrying a gene for antibiotic resistance and n CP is the density of cells which both carry a gene for antibiotic resistance on the chromosome and are infected by a plasmid carrying a gene for antibiotic resistance. We assume that cells which carry a gene for antibiotic resistance are fully immune to the effects of the antibiotic (although we later relax this assumption).
      Table 1

      Parameters used in the model

      Parameter

      Description

      n F

      Density of wild-type/plasmid-free cells ('F' in figures)

      n C

      Density of cells with resistance on the chromosome ('C' in figures)

      n P

      Density of cells infected with a resistance-carrying plasmid ('Plasmids' in figures)

      n B

      Density of cells infected with a non-coding plasmid ('B' in figures)

      n CP

      Density of cells with resistance both on the chromosome and on a plasmid ('CP' in figures)

      n CB

      Density of cells infected with a non-coding plasmid, and with resistance on the chromosome ('CB' in figures)

      N

      Total density of cells

      N = n F +n C +n P +n B +n CP +n CB

      A

      Concentration of antibiotic

      l

      Rate of decay of antibiotic

      θ

      Treatment dosage of antibiotic

      τ

      Interval between antibiotic treatments

      r

      Intrinsic per-capita growth rate of cells

      a

      Extrinsic density-dependent death rate of cells

      m

      Death rate of non-resistant cells, due to antibiotic

      c C

      Cost of antibiotic resistance when gene carried on chromosome

      c P

      Cost of antibiotic resistance when gene carried on plasmid

      x

      Cost of plasmid carriage

      β

      Rate of horizontal transfer of plasmids

      s

      Rate of segregation

      Antibiotic resistance is likely to inflict a cost on the carrying individual [3237], which we denote c C in the case of chromosome-carried resistance and c P in the case of plasmid-carried resistance. The dynamics of chromosome-carried resistance genes are:
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equ3_HTML.gif
      (3)
      Building on previous models of the dynamics of plasmid carriage [18, 19, 38, 39], and assuming that plasmids are fundamentally genomic parasites [1822], we assume that there is a cost x to a cell from plasmid carriage. The dynamics of plasmid-infected cells are therefore:
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equ4_HTML.gif
      (4)
      Cells with resistance on the chromosome can also be infected by plasmids. In this case there will be multiple copies of the gene and the cell will pay the cost of carrying both genes (i.e. c C + c P ). The dynamics for such cells are:
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equ5_HTML.gif
      (5)

      Analytical results

      From equations (1)-(5), we can look at the conditions under which resistance (either carried by plasmids or coded on the chromosome) will invade. The condition under which a cell with resistance on the chromosome will invade, if wild-type cells are at equilibrium n F * = (r-mA)/a (and if there are no other cells in the population), is if
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equa_HTML.gif
      This occurs if the net benefit (Am) from being resistant to the antibiotic is greater than the costs c C of resistance. When the wild-type cells are at equilibrium (and there are no other cells in the population), a plasmid carrying the resistance gene will invade if
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equb_HTML.gif

      This shows that lower costs (in terms of x, s and c P ), higher rates of transfer of the plasmid given by β(r-mA)/a, or greater impact of the antibiotic on wild-type cells (given by mA) will all favour the invasion of plasmid resistance.

      Plasmids will be selected to carry antibiotic resistance genes if the fitness of cells with resistance on a plasmid is greater than the fitness of cells with resistance on a chromosome. Plasmid-bearing cells will grow faster than cells with resistance on the chromosome when both are rare (i.e. n P →0, n C →0, n CP →0) if the following condition is fulfilled:
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equ6_HTML.gif
      (7)
      From inequality 7, we find that the plasmid will be favoured if
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equ7_HTML.gif
      (8)

      This shows that plasmids cannot invade the population if there are no wild-type cells (i.e. n F = 0). As the antibiotic will kill wild-type cells, cells with resistance on the chromosome will be favoured over cells with resistance on a plasmid if the antibiotic has substantially reduced the density n F of wild-type cells. Plasmid-carried resistance will also be favoured over chromosome-carried resistance if the costs of expressing a resistance gene on the chromosome c C , is greater than the cost of expressing a resistance gene on a plasmid c P .

      If the antibiotic fails to completely eradicate the population of wild-type cells, a wild-type population with no plasmids will have an equilibrium density of n F * = (r-Am)/a, assuming that the concentration of antibiotics is not enough to drive them extinct (which requires that r > Am). From this, plasmid-borne resistance will thus be favoured over chromosomal resistance when
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equc_HTML.gif
      This implies that higher per-capita growth rate, or low doses of antibiotic, will favour plasmid-borne resistance. Consequently, plasmids will be favoured when the wild-type density is not much constrained by antibiotic-induced mortality. Plasmids will have a faster growth rate than wild-type cells (cells without any resistance gene or plasmid) if
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equd_HTML.gif
      which will occur when
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Eque_HTML.gif
      When the wild-type cells are at equilibrium n F * = (r-Am)/a then this criteria becomes:
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equf_HTML.gif

      This shows that higher transfer rates β of the plasmid, as well as a greater effect of the antibiotic, or lower overall costs from carrying the plasmid (x +s + c p ) will all result plasmid-carried resistance have a greater over all rate of spread than wild-type cells. Interestingly, greater growth rates r will also favour plasmid-carried resistance to have a larger growth rate than wild-type cells.

      Numerical analysis model

      We now explore how different treatments select for antibiotic resistance to be carried either on a plasmid, or a chromosome, using numerical simulations of our system. We built a numerical simulation to investigate how dosage strength and interval affect whether resistance will be found on the host chromosome or a plasmid. We start with a population of wild-type cells at the carrying capacity (i.e. n F = 1). Following [40], we assume exponential degradation at a rate l for the antibiotic:
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equ8_HTML.gif
      (6)
      We use the Euler method of integration with a time step of 0.001. Once the cells are at the equilibrium density, we then assume that the antibiotic is introduced, at a dosage of θ. We then assume that an amount θ of the antibiotic is re-introduced in regular treatment intervals of τ time-steps. At the beginning of the simulation, we assume that only 0.1% of cells contain the resistance gene, and these are carried either on a plasmid or on the host chromosome (i.e. n C = 0.001 and n P = 0.001 at the start of the simulation). We then run the simulation for 5,000 time-steps, or until the cells have reached a stable dynamical equilibrium. Examples of the dynamics of the model are given in Figure 2, illustrating the case where the interval is either 5 (top panel) or 50 (bottom panel).
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Fig2_HTML.jpg
      Figure 2

      Dynamics of the basic model (consisting of wild-type cells, cells with resistance on the chromosome and plasmid-carried resistance) in response to different treatment regimes. "F" denotes wild-type cells, "P" denotes cells infected with a plasmid carrying resistance, "C" denotes cells with resistance on the chromosome and "CP" stands for cells with resistance on the chromosome that also carry a resistance plasmid. Parameters used are r = 1, a = 1, β = 0.006, s = 0.001, cc = 0.002, cp = 0.002, x = 0.002, m = 0.1 and l = 0.005. The dosage is 50 in both panels, while the treatment interval is τ = 5 time units in the top panel and τ = 50 time units in the bottom panel. Time in our model is scaled, and is thus measured in arbitrary units (AU).

      Figure 3 shows the qualitative outcome for a large set of treatment intervals and dosages. The boundaries were interpolated for increased clarity. Under the conditions shown in Figure 3, chromosomal resistance does not evolve, and that lower intervals between treatments, or high dosages of antibiotics, increasingly favour a population only consisting of plasmid resistance, as the wild-type cells are killed by the antibiotic. However, as shown with our analytical model, the strain that wins depends strongly on the transmission rate (Additional File 1). In the case that the transmission rate beta is low, chromosomal resistance will be favoured if the interval is short and the treatment dosage is high, while as the transmission rate increases, plasmids increasingly dominate across the population. These confirm the results of our analytical model and suggest that transmission is the most important factor promoting plasmid-carried resistance.
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Fig3_HTML.jpg
      Figure 3

      The effect of antibiotic dosage intensity and the interval between treatments on the cell types persisting at equilibrium for the basic model (consisting of wild-type cells, cells with resistance on the chromosome and plasmid-carried resistance). Here, chromosomal resistance can be outcompeted by plasmid-carried resistance. "F" denotes wild-type cells and "P" denotes cells infected with a plasmid carrying resistance. Cell types that are present in the population at a density greater than exceeding 0.001 are shown. The plots were calculated by running the simulation for a number of parameter values for 5,000 time-steps. Lines were then smoothed by interpolation. Parameters used are r = 1, a = 1, β = 0.1, c c = 0.02, c p = 0.02, x = 0.05, m = 0.1, s = 0.001 and l = 0.5.

      Model with non-coding plasmids

      Up to now, we used the simplifying assumption that plasmids would always carry resistance genes. It has been shown previously that plasmids which do not carry the focal gene can have a substantial impact on the invasion of horizontally-transferred genes [39]. While a significant proportion of the genome of plasmids is involved in antibiotic resistance, many plasmids contain no genes involved in antibiotics at all (Figure 1 - see figure legend for details). Thus, we now additionally consider plasmids which are infectious, but do not code for resistance, which we refer to as "non-resistant plasmids". These cells pay a cost of plasmid carriage, but do not carry resistance to the antibiotic, and are therefore susceptible to the antibiotic. The dynamics of wild-type, cells with resistance on the chromosome and cells with resistance on a plasmid now become:
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equ9_HTML.gif
      (9)
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equ10_HTML.gif
      (10)
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equ11_HTML.gif
      (11)
      We assume that non-resistant plasmids are incompatible with resistance plasmids, and hence a cell can be infected by only one type of plasmid. The density of cells infected with non-resistant plasmids n B is given by:
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equ12_HTML.gif
      (12)
      The density of cells which carry resistance on both the chromosome, and on a plasmid, is given by n CP . The dynamics of these cells are:
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equ13_HTML.gif
      (13)
      Finally, we also assume that cells which have resistance on the chromosome can be infected with a non-resistant plasmid, and in this case the cell pays the cost x of bearing a plasmid and the cost c c of antibiotic resistance. The density of such cells is given by n CB and their dynamics is given by:
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Equ14_HTML.gif
      (14)

      In the basic model, plasmids were favoured due to their ability to transfer infectiously. In the absence of antibiotics, non-resistant plasmids will outcompete resistance plasmids, because resistance plasmids additionally pay the cost of antibiotic resistance gene expression c P . As in the basic model, we assume that the antibiotic degrades at an exponential rate described by equation 6. We analyse the system described by equations 9-14 using the numerical method described in the previous section, but here we start at the point where both n F and n B have reached equilibria (i.e. the population already consists of wild-type cells and non-resistant plasmids coexisting), before introducing a small number of cells with resistance on the plasmid, or on the chromosome, respectively (i.e. n P = 0.01, n C = 0.01, n CP = 0, n CB = 0 at the start of the simulation).

      Figure 4 shows examples of the dynamics that result for different treatment regimes. Figure 5 shows the conditions under which resistant plasmids persist under a wide range of different treatment regimes. This shows that plasmid-carried resistance is favoured under more intense, but less frequent, treatment regimes. Figure 5 shows that plasmids that carry resistance are outcompeted by cells with resistance on the chromosome at lower dosages and shorter treatment intervals (Figure 5). In contrast, non-resistant plasmids are favoured only when there is a low dosage of antibiotic applied with long intervals between treatments. Under weak treatments, the non-resistant plasmids outcompete plasmids carrying the resistance gene due to the costs c p of expressing the resistance gene. In general, the results of our extended model support our basic model, and confirm our predictions that plasmid-coded resistance will outcompete chromosomally-coded resistance under long-interval treatment regimes.
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Fig4_HTML.jpg
      Figure 4

      Dynamics of the full model (consisting of wild-type cells, cells with resistance on the chromosome, plasmids carrying resistance and plasmids with no resistance) in response to different treatment regimes. "F" denotes wild-type cells, "P" denotes cells infected with a plasmid carrying resistance, "C" denotes cells with resistance on the chromosome, "B" denotes cells infected plasmids which do not carry the resistance gene, and "CP" and "CB" stand for the chromosomally resistant cell types also carrying the respective plasmid. Parameters used are r = 1, a = 1, β = 0.1, s = 0.001, cc = 0.02, cp = 0.02, x = 0.05, m = 0.1 and l = 0.5. The dosage is 50 in the top two panels and 5 in the bottom panel, while the treatment interval is τ = 1 in the top panel and τ = 80 in the bottom two panels.

      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Fig5_HTML.jpg
      Figure 5

      The effect of antibiotic dosage intensity and the interval between treatments on the cell types persisting at equilibrium for the full model (consisting of wild-type cells, cells with resistance on the chromosome, plasmids carrying resistance and plasmids with no resistance). "F" denotes wild-type cells, "P" denotes cells infected with a plasmid carrying resistance and "C" denotes cells with resistance on the chromosome and "B" denotes cells carrying plasmids that do not code for resistance genes. Here, plasmid-carried resistance is favoured when dosage and interval between treatments are high. Cell types that are present in the population at a density greater than exceeding 0.001 are shown. The plots were calculated by running the simulation for a number of parameter values for 5,000 time-steps. Lines were then smoothed by interpolation. Parameters used are r = 1, a = 1, β = 0.1, cc = 0.02, cp = 0.02, x = 0.05, m = 0.1, s = 0.001 and l = 0.5.

      It is likely that cells that have resistance genes on a plasmid will be less susceptible to the antibiotic than cells which have resistance genes on a plasmid. To explore this, we looked at the case where all cell types could have differing levels of susceptibility. In this case, cells with resistance on the chromosome suffered a mortality of m = 0.002 per cell per unit of antibiotic, while cells with resistance on a plasmids suffered a mortality of m = 0.001 per cell per unit of antibiotic. If a cell carried the resistance gene both on the chromosome, or on a plasmid, we assumed the effective of having both genes would be additive, and the antibiotic-induced mortality of 0.000002 (i.e. the product of the two mortalities). As in the previous models, any cell which did not have either resistance gene suffered a per unit antibiotic mortality of m = 0.1. The results are plotted in Figure 6 (this is in comparison to Figure 5 as the comparison case where resistance is 100% effective). In this case, more intensive treatment regimes (i.e. high dosages, with short intervals between doses) favour a cells with both resistance on a plasmid and on the chromosome (Figure 6). This is in contrast to the case where resistance is 100% effective, where only cells carrying resistance on the plasmid persist at intensive treatment regimes.
      http://static-content.springer.com/image/art%3A10.1186%2F1471-2148-11-130/MediaObjects/12862_2010_1738_Fig6_HTML.jpg
      Figure 6

      The effect of antibiotic dosage intensity and the interval between treatments on the cell types persisting at equilibrium for the full model (consisting of wild-type cells, cells with resistance on the chromosome, plasmids carrying resistance and plasmids with no resistance) without segregation (i.e.s= 0). Here, plasmid-carried resistance is favoured over non-resistant plasmids at higher dosages, but shorter treatment intervals, but are favoured over chromosomal resistance at lower dosages but longer treatment intervals. This figure shows the case where all cells are susceptible to the antibiotic but where "C" cells suffer a mortality of m = 0.002, "P" cells suffer a mortality of m = 0.001 and "CP" cells have a combined mortality of 0.000002 (the product of the two mortalities). Wilt-type cells suffer a mortality m = 0.1. "F" denotes wild-type cells, "P" denotes cells infected with a plasmid carrying resistance and "C" denotes cells with resistance on the chromosome and "B" denotes cells carrying plasmids that do not code for resistance genes. Cell types that are present in the population at a density greater than exceeding 0.001 are shown. The plots were calculated by running the simulation for a number of parameter values for 5,000 time-steps. Lines were then smoothed by interpolation. Parameters used are r = 1, a = 1, β = 0.1, c c = 0.02, c p = 0.02, x = 0.05 and l = 0.5.

      We further explored our results for different parameter value combinations, such as the transmission rate (Additional Files 1 and 2), the segregation rate (Additional File 3), the respective cost of having resistance on a plasmid or a chromosome (Additional File 4) and the mortality cost m and the degradation rate l of the antibiotic (Additional File 5). All graphs can be found in the supplementary material. Our results generally held when changing these parameter values.

      Discussion

      Plasmids are favoured when their rate of horizontal transfer outweighs the costs involved in plasmid carriage to a cell [19]. However, in the presence of less costly plasmids, the benefit of transfer disappears [39]. In our basic model in the absence of non-resistant plasmids, plasmids can persist as long as the rate of transfer is sufficient (Figure 3, Additional File 1). Our full model, incorporating non-resistant plasmids, shows that resistant plasmids outcompete non-resistant plasmids under higher dosages, and outcompete cells which resistance on the chromosome under longer intervals between treatments (Figures 5 and 6). This is due to horizontal gene transfer: wild-type cells are common when treatment is weak and their presence allows plasmids to invade. Resistance genes encoded by the host genome are favoured by treatment regimes that result in a continual presence of antibioitic. We can understand this intuitively by considering that the infectious maintenance of plasmids in a population consisting only of plasmid-carrying and wild-type cells requires the availability of plasmid-free cells. This is because carrying a plasmid and expressing resistance has a cost which has to be offset by the benefit of infecting susceptible cells. When susceptible cells are killed by antibiotics, plasmids suffer indirectly as they have fewer cells to infect. It is worth noting that our model assumes relatively efficient horizontal transfer, and thus our results may change if plasmids do not transmit themselves well. Our model would therefore predict that plasmids with lower transmission rates would not carry genes for antibiotic resistance. However the generality of our model of horizontal transfer of plasmids means that our findings will apply to other vehicles of gene transfer, such as phages.

      We assumed that chromosomal- and plasmid-borne resistant types were phenotypically identical. This may not be the case [41] as changes in gene dosage can influence the resistance level [4247]. We predict that higher gene dosage from multi-copy plasmids would enlarge the treatment space beneficial to plasmids, due to the additional resistance compared to chromosomally carried resistance genes. In the absence of antibiotics, hosts carrying multiple copies of resistance genes would be at a disadvantage due to the extra cost of producing more antibiotic resistance proteins. An important result of our study is to illustrate how the ecology of the bacteria, in this case the density of wild-type cells (mediated by the concentration of antibiotic), affects the conditions under which plasmid-borne resistance is favoured.

      Our model makes qualitative predictions regarding the conditions under which antibiotic resistance genes will be selected to be carried on plasmids, as opposed to on the chromosome. Specifically, we were interested in asking whether the intensity and interval between treatments will select for resistance to be on a plasmid, or on a chromosome. However, we find that our results generally hold under a wide range of different parameter ranges (see supplementary material for details). There is a wide range of literature available with which such models could be parameterised, and it would be possible to design models to test when selection will result in resistance to be carried on plasmids. For example, it is possible to add more realism to model the pharmacodynamics of drugs and how they kill bacteria [e.g. [40, 48]], and there are some studies from which the costs of antibiotic resistance can be estimated [for example, see references in [49]]. Rates of horizontal gene transfer range widely and in the case of conjugation can range between 10-3 and 10-1 per donor in biofilms and more than 10-5 in water or soil [5052]. As plasmids have been described as the Achilles' heal of drug resistant bacteria [53, 54], understanding the conditions under which antibiotic resistant genes are carried by plasmids could help to develop strategies to minimise the spread of resistance [53]. Developing and parameterising more explicit models in combination with laboratory studies, in contrast to the simple qualitative models described here, will be essential to developing treatments than minimise the rates of transfer.

      The use of multiple antibiotics and antibiotic cycling has been proposed as a way to prevent the evolution of antibiotic resistance [e.g. [55]]. A previous study showed that plasmids could play a role in the acquisition of multiple drug resistance by repeated gene transfer [56]. Thus, plasmids may reduce the effect of antimicrobial cycling, as multiple resistance could be acquired more quickly. In the case of multidrug therapy, we would predict that the results of our model would hold, and longer intervals between treatment regimes would favour plasmid-carried resistance, which would in turn favour the evolution of multiple resistance.

      Given the biological variability in the mechanisms of resistance, horizontal transfer and persistence, some of our assumptions may not hold in specific settings. Our work rests on costly antibiotic resistance and should thus apply particularly to mechanisms that need to be strongly expressed or that require a lot of cellular energy. When resistance proteins catalyse reactions that modify the antibiotic, such as in the case of β-lactamases, chloramphenicol acetyltransferase or aminoglycoside acetyltransferases, a higher expression level is expected to translate to a higher resistance level and the cell will have to pay high protein synthesis costs in order to be resistant. The same is true for efflux pumps, such as ABC -transporters or MF-type pumps [57]. The cost of resistance plasmids, including the cost of maintenance and expression, could be significantly ameliorated in experimental evolution experiments [5860]. Based on the above arguments and the fact that even small costs are significant over evolutionary timescales, we expect there to be a large number of environments in which our findings would hold despite this effect, but it illustrates that with respect to the clinical view on this fundamental problem, more detailed case-by-case analyses will be required.

      Conclusions

      To the best of our knowledge, our study represents the first model to explicitly examine the effect of treatment regimes plasmid-borne antibiotic resistance. We show that high transmission favours plasmid-carried resistance, in the absence of competing non-resistant plasmids. When other plasmids are present, low frequency treatments favour plasmid-carried resistance over chromosomal resistance, but a high intensity of antibiotic is needed for resistance plasmids to outcompete non-resistance plasmids. Targeting plasmid spread has been proposed to manage antibiotic resistance [54] and our results suggest that paying attention to the treatment regime is an essential requirement of any such strategy. The solution, as recommend by Paul Ehrlich almost a century ago, is to "hit hard and hit quickly" [61]. Gaining a quantitative understanding of the dynamics of plasmids will allow us to understand harmful patterns of antibiotic use more effectively.

      Declarations

      Acknowledgements

      We thank Sam Brown, Jeffrey Lawrence, Sorcha Mc Ginty, Roland Regoes, jeff smith and three anonymous reviewers for helpful comments on the manuscript. We also thank Evandro Ferrada for help obtaining the genome data used in Figure 1, the Swiss National Science Foundation (grants 31003A-125457 and PZ00P3-121800) and the University of Zürich (both to DJR) for funding.

      Authors’ Affiliations

      (1)
      Institute of Evolutionary Biology and Environmental Studies, University of Zürich
      (2)
      Swiss Institute of Bioinformatics, Quartier Sorge, Bâtiment Génopode

      References

      1. WHO: Overcoming antibiotic resistance. In World Health Organisation Report in Infectious Disease. Geneva: World Health Organisation; 2000.
      2. Lipsitch M, Bergstrom CT: Modeling of antibiotic resistance in the ICU - U.S. slant. In Infection control in the ICU environment. Edited by: Weinstein RA. Bonten M: Kluwer; 2002.
      3. Bergstrom CT, Feldgarden M: The ecology and evolution of antibiotic-resistant bacteria. In Evolution in Health and Disease. 2nd edition. Edited by: Stearns S, Koella J. Oxford: Oxford University Press; 2008.
      4. Lipsitch M: The rise and fall of antimicrobial resistance. Trends in Microbiology 2001,9(9):438–444.PubMedView Article
      5. Lipsitch M, Samore MH: Antimicrobial use and antimicrobial resistance: a population perspective. Emerging Infectious Diseases 2002,8(4):347–354.PubMedView Article
      6. Thomas CM, Nielsen KM: Mechanisms of, and Barriers to, Horizontal Gene Transfer between Bacteria. Nat Rev Micro 2005,3(9):711–721.View Article
      7. Falkow S, Citarella RV, Wohlhieter JA, Watanabe T: The molecular nature of R-factors. Journal of Molecular Biology 1966,17(1):102–116.PubMedView Article
      8. Sugino Y, Hirota Y: Conjugal fertility associated with resistance factor R in Escherichia coli . Journal of Bacteriology 1962,84(5):902–910.PubMed
      9. Clewell DB, Yagi Y, Bauer B: Plasmid-determined tetracycline resistance in Streptococcus faecalis: evidence for gene amplification during growth in presence of tetracycline. Proceedings of the National Academy of Sciences 1975, 72:1720–1724.View Article
      10. Livermore DM: Beta-Lactamases in laboratory and clinical resistance. Clinical Microbiolologocal Reviews 1995,8(4):557–584.
      11. Johnsen PJ, Simonsen GS, Olsvik O, Midtvedt T, Sundsfjord A: Stability, persistence, and evolution of plasmid-encoded VanA glycopeptide resistance in enterococci in the absence of antibiotic selection in vitro and in mice. Microbial Drug Resistance 2002, 8:161–170.PubMedView Article
      12. Weigel LM, Clewell DB, Gill SR, Clark NC, McDougal LK, Flannagan SE, Kolonay JF, Shetty J, Killgore GE, Tenover FC: Genetic Analysis of a High-Level Vancomycin-Resistant Isolate of Staphylococcus aureus . Science 2003,302(5650):1569–1571.PubMedView Article
      13. Furuya EY, Lowy FD: Antimicrobial-resistant bacteria in the community setting. Nat Rev Micro 2006,4(1):36–45.View Article
      14. Robicsek A, Jacoby GA, Hooper DC: The worldwide emergence of plasmid-mediated quinolone resistance. The Lancet Infectious Diseases 2006, 6:629–640.PubMedView Article
      15. Dionisio F, Matic I, Radman M, Rodrigues OR, Taddei F: Plasmids spread very fast in heterogeneous bacterial communities. Genetics 2002,162(4):1525–1532.PubMed
      16. Hughes VM, Datta N: Conjugative plasmids in bacteria of the 'pre-antibiotic' era. Nature 1983,302(5910):725–726.PubMedView Article
      17. Turner PE, Cooper VS, Lenski RE: Tradeoff between horizontal and vertical modes of transmission in bacterial plasmids. Evolution 1998,52(2):315–329.View Article
      18. Bergstrom CT, Lipsitch M, Levin BR: Natural selection, infectious transfer and the existence conditions for bacterial plasmids. Genetics 2000, 155:1505–1519.PubMed
      19. Lili LN, Britton NF, Feil EJ: The persistence of parasitic plasmids. Genetics 2007,177(1):399–405.PubMedView Article
      20. Rankin DJ, Rocha EPC, Brown SP: What genes are carried on mobile elements, and why? Heredity 2011, 106:1–10.PubMedView Article
      21. Levin BR, Stewart FM: Probability of establishing chimeric plasmids in natural populations of bacteria. Science 1977, 196:218–220.PubMedView Article
      22. Stewart FM, Levin BR: The population biology of bacterial plasmids: a priori conditions for the existence of conjugationally transmitted factors. Genetics 1977, 87:209–228.PubMed
      23. Anderson RM, May RM: Population biology of infectious diseases: Part I. Nature 1979,280(5721):361–367.PubMedView Article
      24. Anderson RM, May RM: Infectious Diseases of Humans: Dynamics and Control. Oxford: Oxford University Press; 1992.
      25. May RM, Anderson RM: Population biology of infectious diseases: Part II. Nature 1979,280(5722):455–461.PubMedView Article
      26. Smith J: The social evolution of bacterial pathogenesis. Proceedings of the Royal Society of London Series B-Biological Sciences 2001, 268:61–69.View Article
      27. Dionisio F, Conceicao IC, Marques ACR, Fernandes L, Gordo I: The evolution of a conjugative plasmid and its ability to increase bacterial fitness. Biology Letters 2005,1(2):250–252.PubMedView Article
      28. Yates CM, Shaw DJ, Roe AJ, Woolhouse MEJ, Amyes SGB: Enhancement of bacterial competitive fitness by apramycin resistance plasmids from non-pathogenic Escherichia coli . Biology Letters 2006,2(3):463–465.PubMedView Article
      29. Ellis RJ, Lilley AK, Lacey SJ, Murrell D, Godfray HCJ: Frequency-dependent advantages of plasmid carriage by Pseudomonas in homogeneous and spatially structured environments. Isme J 2007,1(1):92–95.PubMedView Article
      30. Slater FR, Bailey MJ, Tett AJ, Turner SL: Progress towards understanding the fate of plasmids in bacterial communities. FEMS Microbiol Ecology 2008, 66:3–13.View Article
      31. Nogueira T, Rankin DJ, Touchon M, Taddei F, Brown SP, Rocha EP: Gene mobility drives the evolution of bacterial cooperation and virulence. Current Biology 2009, 19:1683–1691.PubMedView Article
      32. Lenski RE: Bacterial evolution and the cost of antibiotic resistance. International Microbiology 1998, 1:265–270.PubMed
      33. Andersson DI, Levin BR: The biological cost of antibiotic resistance. Current Opinion in Microbiology 1999,2(5):489–493.PubMedView Article
      34. Björkman J, Nagaev I, Berg OG, Hughes D, Andersson DI: Effects of Environment on Compensatory Mutations to Ameliorate Costs of Antibiotic Resistance. Science 2000,287(5457):1479–1482.PubMedView Article
      35. Schrag SJ, Perrot V, Levin BR: Adaptation to the Fitness Costs of Antibiotic Resistance in Escherichia coli. Proceedings: Biological Sciences 1997,264(1386):1287–1291.View Article
      36. Gagneux S, Long CD, Small PM, Van T, Schoolnik GK, Bohannan BJM: The Competitive Cost of Antibiotic Resistance in Mycobacterium tuberculosis. Science 2006,312(5782):1944–1946.PubMedView Article
      37. Andersson DI: The biological cost of mutational antibiotic resistance: any practical conclusions? Current Opinion in Microbiology 2006, 9:461–465.PubMedView Article
      38. Rankin DJ, Bichsel M, Wagner A: Mobile DNA can drive lineage extinction in prokaryotic populations. Journal of Evolutionary Biology 2010, 23:2422–2341.PubMedView Article
      39. Mc Ginty SE, Rankin DJ, Brown SP: Horizontal gene transfer and the evolution of bacterial cooperation. Evolution 2011, 65:21–32.PubMedView Article
      40. Regoes R, Wiuff C, Zappala RM, Garner KN, Baquero F, Levin BR: Pharmacodynamic functions: a multiparameter approach to the design of anitbiotic treatment regimens. Antimicrobial Agents and Chemotherapy 2004,48(10):3670–3676.PubMedView Article
      41. Lobner-Olesen A: Distribution of minichromosomes in individual Escherichia coli cells: implications for replication control. EMBO J 1999,18(6):1712–1721.PubMedView Article
      42. Nordstrom K, Uhlin BE: Runaway-Replication Plasmids as Tools to Produce Large Quantities of Proteins from Cloned Genes in Bacteria. Nat Biotech 1992,10(6):661–666.View Article
      43. Makrides SC: Strategies for Achieving High-Level Expression of Genes in Escherichia coli. Microbiological Reviews 1996,60(3):512–538.PubMed
      44. Morino T, Morita M, Seya K, Sukenaga Y, Kato K, Nakamura T: Construction of a runaway vector and its use for a high-level expression of a cloned human superoxide dismutase gene. Applied Microbiology and Biotechnology 1988, 28:170–175.View Article
      45. Andersson DI, Hughes D: Gene amplification and adaptive evolution in bacteria. Annual Review of Genetics 2009, 43:167–195.PubMedView Article
      46. Sandegren L, Andersson DI: Bacterial gene amplification: implications for the evolution of antibiotic resistance. Nat Rev Micro 2009,7(8):578–588.View Article
      47. Sun S, Berg OG, Roth JR, Andersson DI: Contribution of Gene Amplification to Evolution of Increased Antibiotic Resistance in Salmonella typhimurium. Genetics 2009, 182:1183–1195.PubMedView Article
      48. Vinks A, Derendorf H, Mouton J, (eds): Fundamentals of Antimicrobial Pharmacokinetics and Pharmacodynamics. New York: Springer; 2011.
      49. Andersson DI, Hughes D: Antibiotic resistance and its cost: is it possible to reverse resistance? Nature reviews Microbiology 2010,8(4):260–271.PubMed
      50. Sorensen SJ, Bailey M, Hansen LH, Kroer N, Wuertz S: Studying plasmid horizontal transfer in situ: a critical review. Nature reviews Microbiology 2005,3(9):700–710.PubMedView Article
      51. Hausner M, Wuertz S: High rates of conjugation in bacterial biofilms as determined by quantitative in situ analysis. Applied and Environmental Microbiology 1999,65(8):3710–3713.PubMed
      52. Dahlberg C, Bergstrom M, Hermansson M: In Situ Detection of High Levels of Horizontal Plasmid Transfer in Marine Bacterial Communities. Applied and Environmental Microbiology 1998,64(7):2670–2675.PubMed
      53. Williams JJ, Hergenrother PJ: Exposing plasmids as the Achilles' heel of drug-resistant bacteria. Current Opinion in Chemical Biology 2008,12(4):389–399.PubMedView Article
      54. Amabile-Cuevas CF, Heinemann JA: Shooting the messenger of antibiotic resistance: plasmid elimination as a counter-evolutionary tactic. Drug Discovery Today 2004,9(11):465–467.PubMedView Article
      55. Bergstrom CT, Lo M, Lipsitch M: Ecological theory suggests that antimicrobial cycling will not reduce antimicrobial resistance in hospitals. Proceedings of the National Academy of Sciences of the United States of America 2004, 101:13285–13290.PubMedView Article
      56. Willms AR, Roughan PD, Heinemann JA: Static recipient cells as reservoirs of antibiotic resistance during antibiotic therapy. Theoretical Population Biology 2006,70(4):436–451.PubMedView Article
      57. Webber MA, Piddock LJV: The importance of efflux pumps in bacterial antibiotic resistance. Journal of Antimicrobial Chemotherapy 2003,51(1):9–11.PubMedView Article
      58. Lenski R: Bacterial evolution and the cost of antibiotic resistance. International Microbiology 1998, 1:265–270.PubMed
      59. Modi RI, Adams J: Coevolution in bacterial-plasmid populations. Evolution 1991,43(3):656–667.View Article
      60. Dahlberg C, Chao L: Amelioration of the cost of conjugative plasmid carriage in Eshericha coli K12. Genetics 2003, 165:1641–1649.PubMed
      61. Ehrlich P: Address in Pathology on Chemotherapeutics: Scientific Principles, Methods and Results. The Lancet 1913, 182:445–451.View Article
      62. Liu B, Pop M: ARDB-Antibiotic Resistance Genes Database. Nucleic Acids Research 2009, 37:D443–447.PubMedView Article

      Copyright

      © Svara and Rankin; licensee BioMed Central Ltd. 2011

      This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

      Advertisement