(C) PLOS One
This story was originally published by PLOS One and is unaltered.
. . . . . . . . . .



Segregational drift hinders the evolution of antibiotic resistance on polyploid replicons [1]

['Ana Garoña', 'Institute Of General Microbiology', 'Kiel University', 'Kiel', 'Mario Santer', 'Research Group Stochastic Evolutionary Dynamics', 'Department Of Evolutionary Theory', 'Max Planck Institute For Evolutionary Biology', 'Plön', 'Nils F. Hülter']

Date: 2023-08

The emergence of antibiotic resistance under treatment depends on the availability of resistance alleles and their establishment in the population. Novel resistance alleles are encoded either in chromosomal or extrachromosomal genetic elements; both types may be present in multiple copies within the cell. However, the effect of polyploidy on the emergence of antibiotic resistance remains understudied. Here we show that the establishment of resistance alleles in microbial populations depends on the ploidy level. Evolving bacterial populations under selection for antibiotic resistance, we demonstrate that resistance alleles in polyploid elements are lost frequently in comparison to alleles in monoploid elements due to segregational drift. Integrating the experiments with a mathematical model, we find a remarkable agreement between the theoretical and empirical results, confirming our understanding of the allele segregation process. Using the mathematical model, we further show that the effect of polyploidy on the establishment probability of beneficial alleles is strongest for low replicon copy numbers and plateaus for high replicon copy numbers. Our results suggest that the distribution of fitness effects for mutations that are eventually fixed in a population depends on the replicon ploidy level. Our study indicates that the emergence of antibiotic resistance in bacterial pathogens depends on the pathogen ploidy level.

Prokaryotic replicons—chromosomes and plasmids—may be present in the cell in multiple copies, a phenomenon termed polyploidy. Here we studied the effect of replicon copy number, i.e., polyploidy level, on the emergence of beneficial traits in prokaryotes, with a focus on the evolution of antibiotic resistance. We discovered that beneficial mutations encoded in multicopy replicons are bound to be frequently lost due to stochasticity during cell division that ultimately increases genetic drift. Our study presents empirical evidence for the effect of polyploidy on the rate of adaptation in prokaryotes. A comparison to simulations of a mathematical model confirms our understanding of the allele dynamics. We conclude that drug treatment strategies against pathogenic bacteria should consider the ploidy level of the organism.

Funding: AG is supported by the International Max Planck Research School (IMPRS) for Evolutionary Biology. TD is supported by the Leibniz Science Campus EvoLUNG. MS is a member of the International Max Planck Research School for Evolutionary Biology and gratefully acknowledges the benefits provided by the program. This work was furthermore funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) - project number 418432175 (to HU) and European Union (ERC) pMolEvol, grant number 101043835 (to TD). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Here we investigate the dynamics of a novel antibiotic resistance allele in an experimentally evolving bacterial population of Acinetobacter baylyi and compare the results between populations where the allele is encoded either on a polyploid (multicopy) plasmid or on a monoploid chromosome. In our experimental setup, we manipulate the degree of stochasticity present at the level of the cell population using two measures: 1) introducing a low initial frequency of the novel resistance allele, and 2) applying serial population bottlenecks on the evolved population. This experimental design allows us to systematically study how segregational drift interacts with stochasticity in the population dynamics of cells carrying the novel allele. We furthermore combine the experimental results with a mathematical model of bacterial population genetics that is specifically tailored to the experimental setup to directly probe our mechanistic understanding of the dynamics of alleles on polyploid replicons.

The presence of multiple replicon copies in the cell has important implications for bacterial evolution. Experimental evolution of bacteria has shown that the mutational supply is positively associated with the replicon copy number, i.e., the probability of novel mutations can be higher for polyploid replicons (as demonstrated for multicopy plasmids in ref. [ 21 ]). Polyploidy furthermore allows for intracellular genetic diversity, also termed heterozygosity [ 22 – 24 ]. Alleles in monoploid chromosomes are inevitably inherited to both daughter cells during cell division. In contrast, the segregation of alleles emerging in one out of many copies of a polyploid replicon is a neutral process and it may be unbalanced, such that only one daughter cell inherits the novel allele. The effect of stochasticity on novel plasmid alleles during cell division has been termed segregational drift [ 25 ]. The dynamics of alleles encoded in polyploid chromosomes (and multicopy plasmids) are thus affected by processes at two hierarchical levels: within the cell and within the population [ 25 – 29 ]. The mechanism of replication and the mode of segregation are important determinants of the allele dynamics on polyploid replicons [ 28 , 29 ]. A recent model on the fixation process of beneficial allele dynamics shows that heterozygous cells can persist for an extended period of time during the fixation process of high copy number replicons [ 29 ]. The coexistence of multiple plasmid alleles within a cell enables the host population to survive in fluctuating environmental conditions, e.g., under alternating drug pressures selecting for multiple antibiotic resistances [ 23 ]. The replicon copy number has furthermore important implications for the dynamics of rare novel alleles. Theory predicts that segregational drift will decrease the establishment probability of beneficial alleles when encoded in multicopy plasmids, at least in the absence of gene dosage effects [ 25 , 27 ]. Nonetheless, empirical evidence for the workings of segregational drift on the dynamics of beneficial alleles remains lacking.

Prokaryotic genomes comprise chromosomes and extrachromosomal genetic elements, e.g., plasmids, collectively termed replicons. The number of chromosome copies in the cell is tightly controlled during replication and is synchronised with cell division [ 9 , 10 ]. Several prokaryotic taxa are known to harbour polyploid chromosomes. This includes cyanobacteria, e.g., Synechococcus elongatus that has a chromosome copy number of 3–4 [ 10 , 11 ]. Other examples are the human pathogens Borrelia hermsii and B. burgdorferi, which have a chromosome copy number of 4–14 and 10–20, respectively [ 12 ]. The number of chromosome copies in the cell may furthermore depend on the growth phase and nutrient conditions (e.g., refs. [ 13 , 14 ]); indeed, also monoploid bacteria, such as E. coli, can be polyploid (or rather mero-oligoploid) during exponential growth [ 15 , 16 ]. The number of plasmid copies in a host cell depends on the plasmid type and the host genetics. Low copy number plasmids are found in 1–5 copies in the cell, while high copy number plasmids may reach 200 copies in the cell [ 17 , 18 ]. Similarly to chromosomes, the plasmid copy number may vary depending on the growth conditions [ 19 , 20 ].

Bacterial adaptation to novel environmental conditions depends on the availability of beneficial alleles and the dynamics of their proliferation within the population. However, novel alleles are initially rare in the population and therefore–even if they are adaptive–prone to stochastic loss due to genetic drift [ 1 ]. The probability that beneficial alleles survive the early proliferation phase and become established is termed ’establishment probability’. The probability of establishment (and eventually fixation) of beneficial alleles is a key concept in theoretical population genetics [ 2 , 3 ], dating back to the very early days of the field [ 1 ]. Understanding the determinants of allele establishment in bacterial populations is pivotal in the context of antibiotic resistance evolution. Antibiotic treatment confers a growth advantage to bacteria carrying resistance alleles, which might establish and subsequently rise to high numbers. The emergence of novel resistant strains of diverse human and livestock pathogens due to antibiotic treatment has worrisome consequences for global human health [ 4 ]. Sustainable treatment strategies are needed [ 5 ], for the development of which a profound understanding of the dynamics of beneficial alleles–including the early stochastic phase–is necessary. Several recent studies have supplied empirical evidence for stochastic loss of adaptive antibiotic resistance genotypes and factors that may be involved in their establishment [ 6 – 8 ]. However, the effect of bacterial ploidy on the establishment of antibiotic resistance remains understudied.

Results

Rare alleles in multicopy replicons are prone to rapid loss To experimentally quantify the effect of the replicon copy number on the fate of beneficial alleles, we studied the dynamics of an allele conferring resistance to the antibiotic kanamycin (nptII). The beneficial allele was introduced either on a polyploid (multicopy) plasmid or a monoploid chromosome in the model organism A. baylyi. Performing an experimental evolution experiment with serial transfers, we examined the effects of the initial allele frequency (i.e., the number of novel allele copies), the strength of the population bottlenecks, and the selection regime on allele dynamics (Fig 1A). We compare the results for three initial frequencies of cells containing the novel allele: f 0 = 10−4 (ca. 105 cells, high), f 0 = 10−6 (ca. 103 cells, moderate), and f 0 = 10−7 (ca. 102 cells, low). The control of the initial frequency of cells containing the novel plasmid allele was achieved via natural transformation employing varying donor DNA concentrations (700, 7, or 0.7 ng/μl DNA). The novel nptII allele is introduced into an ancestral population that carries the stable model plasmid pTAD-R, which has a copy number of ca. 15. Homologous recombination between the donor DNA and the ancestral plasmid creates heterozygous cells carrying both the ancestral and the novel alleles, where the intracellular frequency of the novel allele is one plasmid copy out of 15 due to the one-hit kinetics of natural transformation (Fig 1A; ref. [28]). For the chromosomal allele, we controlled the initial allele frequency by mixing strains carrying either the novel or the ancestral allele at the abovementioned frequencies. Using this approach, we create an ancestral population that has the same initial number of novel allele copies and the same initial frequency of cells carrying the novel allele for the polyploid and the monoploid replicons. The initial novel allele frequency from the total replicons in the population is consequently lower for the plasmid experiments. The evolved populations were serially transferred using two population bottleneck sizes: 1:100 (b = 10−2; ca. 107 cells, weak bottleneck) and 1:1000 (b = 10−3; ca. 106 cells, strong bottleneck) for ca. 200 generations for the weak bottleneck and 300 for the strong bottleneck (30 serial transfers). The experiment was conducted with six replicates for each experimental condition. The fitness effect of the novel plasmid allele in the absence of antibiotics was estimated from the relative fitness of cells that are homozygous for the novel allele and cells that are homozygous for the ancestral allele; the results show a mean relative fitness of w = 0.9884 ± 0.0224 SD (or selective coefficient of s = -0.0116), which was not significantly different from w = 1 (or s = 0; p = 0.4375, using Wilcoxon test, n = 6; S2 Fig). The novel chromosomal allele had a mean relative fitness of w = 0.9892 ± 0.0086 SD, with no significant difference from w = 1 (i.e., s = 0; p = 0.0625, using Wilcoxon test, n = 6; S2 Fig). Thus, in the absence of antibiotics, changes in allele frequencies in the population over time are expected to stem from neutral processes rather than fitness differences between genotypes in the population. Changes in the allele frequency due to horizontal transfer via natural transformation are negligible in our system [28]. PPT PowerPoint slide

PNG larger image

TIFF original image Download: Fig 1. Integrating multicopy replicon allele dynamics in vivo and in silico. (A) Our experimental system enables us to create and follow bacterial populations encoding a novel allele in a plasmid or chromosomal replicon over the course of an evolution experiment. The novel allele employed is a kanamycin resistance gene (nptII), which encodes aminoglycoside 3’-phosphotransferase II [49]. The ancestral plasmid allele is a green fluorescence protein (gfp) whose transcription is under the control of a LacI repressor [28]. The chromosomal ancestral allele is a tetracycline resistance gene (tetA). The bacterial populations were serially passaged under different experimental conditions for approximately 200 generations. The frequency of cells carrying the novel allele was followed by selective plating on kanamycin-containing medium, where each viable cell will form a colony. For the plasmid replicon, the presence of both ancestral and novel allele, that is homozygotes and heterozygotes, was distinguished due to fluorescence conferred by the ancestral allele. Bacterial colonies were visually inspected employing blue light. Note that this is an assessment between homozygotic and heterozygotic colonies rather than an exact quantification of the number of alleles present. (B) Allele dynamics of representative populations evolved under non-selective conditions. The white arrow marks the point at which the frequencies of heterozygous and novel homozygous cells are equal. The plots start at the end of the first day of growth, where the final transformant frequencies, according to the amount of donor DNA added, have been achieved. The remaining experimental populations are shown in S3 and S4 Figs. (C) Schematic diagram of the mathematical model of plasmid allele dynamics under periodic population bottlenecks (transfers). Bacterial populations grow until they reach a carrying capacity N c before each transfer. A proportion b of cells is randomly chosen for the next growth phase. During each growth phase, cells carrying only the ancestral allele (ancestral homozygotes) replicate at a rate that is 1-s times smaller than the rate of cells carrying the novel allele (heterozygotes and novel homozygotes). (D) Comparison between the allele dynamics in evolution experiments using the pTAD-R plasmid and in model simulations with replicon copy number n rc = 15 under non-selective conditions (s = 0) for high initial frequencies, f 0 = 10−4, and weak bottlenecks, b = 10−2. Bold green (red) lines show the median of the heterozygote (homozygote) frequencies before each transfer from 30 simulations for both modes of replication; shaded areas show the corresponding interquartile ranges. Thin lines show the heterozygote frequency from the six replicates of the evolution experiment. Simulation results were obtained by integrating the differential equations for bacterial growth and sampling new populations at each transfer from a multinomial distribution defined by the genotype frequencies at carrying capacity. https://doi.org/10.1371/journal.pgen.1010829.g001 To demonstrate the effects of the initial frequency of cells carrying the novel allele and the bottleneck strength, we first present the observed allele dynamics under non-selective conditions for a set of representative populations (Fig 1B). Note that the novel allele has no significant fitness effect under neutral conditions (S2 Fig). Our results show that with a high initial frequency (f 0 = 10−4) and a weak bottleneck (b = 10−2), the dynamics of the novel allele in both plasmid and chromosome are characterized by a stable persistence during the first ca. 20 transfers, followed by a decrease in the fraction of kanamycin-resistant cells (i.e., mutant cells). An eventual decrease in the frequency of the rare novel allele in the population is to be expected under genetic drift; the fixation probability of a neutral allele corresponds to the initial allele frequency and is thus very low [30]. Segregation of polyploid replicon copies reduces the number of mutant cells (through flux into the two homozygous cell types) and thus increases the effect of stochasticity. Indeed, the frequency of kanamycin-resistant cells decreases faster if the allele is on the multicopy plasmid than if it is on the chromosome (Fig 1B). At the end of the experiment, the novel plasmid allele was lost in most replicates, while the chromosomal allele was present in the population at a frequency of ca. 5×10-7. Our system enables us to distinguish between heterozygous and homozygous cells for the novel plasmid allele. The subpopulation of cells that carries the novel allele initially consists of heterozygotes in which the novel allele is present in only one plasmid copy per cell. Over the course of the experiment, heterozygotes decay due to the random distribution of plasmid copies at cell divisions, which eventually causes the loss of heterozygous cells. For the highest initial frequency, f 0 = 10−4, kanamycin-resistant homozygous cells arise due to segregational drift and are maintained over several transfers if the applied bottlenecks are weak (b = 10−2). Our results show that when the initial frequency of the cells carrying the novel allele is lower, the novel allele is maintained in the population for shorter times (Fig 1B; f 0 = 10−6 and f 0 = 10−7). The bottleneck size applied during the experiment has a significant effect on the allele dynamics, which are characterized by a rapid loss in the populations transferred with a strong bottleneck (b = 10−3) compared to those populations evolved with a weak bottleneck. Note that with a strong bottleneck, fewer cells are transferred than with a weak bottleneck, increasing the strength of drift. At the same time, the bottleneck applied also has an effect on the number of generations elapsed, as with the stronger bottleneck more cell divisions will occur until the carrying capacity is reached. Therefore, the loss of the novel allele is more rapid (in terms of number of transfers) when a strong bottleneck is applied. To sum up, under non-selective conditions, both plasmid and chromosomal alleles decrease over time and eventually are lost in most populations by the end of the experiment. Nonetheless, the loss of the plasmid allele is faster in all tested combinations of initial frequencies and population bottlenecks. Our results are in agreement with previous observations on the effect of segregational drift on the dynamics of alleles encoded in polyploid replicons [25,27,28].

Mathematical modelling generates predictions on the plasmid replication mode based on the empirical plasmid allele dynamics To formulate the theoretical expectation for the allele dynamics in the experiments, we adapted a previous mathematical model of vertical replicon inheritance in bacterial populations to match the experimental design of the evolution experiment (Fig 1C; refs. [27,29]). The genotype of a cell is defined by the number of novel and ancestral replicon copies at cell birth, with the total number of novel replicon copies ranging between 0 and the replicon copy number n rc . Cells in the simulations of monoploid replicons (n rc = 1), such as the chromosome in A. baylyi, harbor either the ancestral or the novel allele. Cells in the simulations of multicopy replicons (n rc > 1), e.g., model plasmid pTAD-R, can be homozygous for the ancestral or for the novel allele, or heterozygous, carrying both allele variants in various proportions. We consider thus all n rc +1 possible cell types. Our model includes two possibilities for the mode of replication: under the regular replication mode, all n rc replicon copies are duplicated before a cell divides. Under the random replication mode, single copies are selected randomly for replication in a repetitive process until the cell carries 2 × n rc copies (refs. [27,29]; see Methods for details). For either mode of replication, the replicon copies are distributed randomly into the two daughter cells at cell division (similarly to pTAD-R), where each cell inherits n rc copies. Population growth between bottlenecks is modelled by deterministic logistic growth, where we only account for cell replication but not for cell death. At bottlenecks, cells are randomly sampled, introducing stochasticity (see detailed model description in the methods). The replication mode of the pTAD-R backbone has not been examined experimentally so far. To determine which modelled mode of replication better matches the experimental system, we simulated the segregation dynamics of a novel allele carried by a polyploid replicon and compared the simulated dynamics to the experimental results. The simulation was performed for both replication modes under neutrality (s = 0) with a replicon having n rc = 15, as model plasmid pTAD-R. In this simulation, we opted for the highest initial frequency, f 0 = 10−4, and the weak bottleneck, b = 10−2, where stochastic effects other than segregational drift are weakest. The model results show that the frequency of heterozygous cells rapidly decreases at a constant rate, while the frequency of homozygotes first increases and then reaches an equilibrium around transfer 15 (Fig 1D). The comparison between the frequency of heterozygotes and homozygotes from simulations under the mode of regular replication and the experimental results shows a good qualitative agreement for the first 11 transfers (Fig 1D). From transfer 12 onwards, the model underestimates the frequency of heterozygous cells by up to ten-fold. The deviation between the model prediction of the homozygotes frequency and the experimental results increases from transfer 15 when the homozygotes frequency in the experiment decayed. Simulations using the random replication mode overestimate the decay of heterozygotes’ frequency as well as the increase in homozygotes’ frequency already during the very first transfers (Fig 1D). To test for possible reasons for the observed deviations, we performed further simulations under the regular replication mode, assuming a different plasmid copy number (n rc between 10 and 30) or a slight negative fitness effect of 1–5% of the novel allele. Increasing the copy number in the simulations yields a better fit to the heterozygote frequencies (S1A Fig); nonetheless, the experimental quantification of pTAD-R copy number does not suggest a large deviation from n rc = 15 [28]. Simulations with a negative selection parameter (approximately between s = -0.01 and s = -0.03) yield a better fit to the homozygote frequency observed in the experiment (S1B Fig). Indeed, while the experimental results suggested no significant fitness effect of the novel plasmid allele, the average relative fitness suggested a slight deviation from a neutral regime (S2 Fig). Note that random replication always entails faster segregation dynamics than regular replication [29]. Hence achieving a good fit between the simulation with random replication and the experimental results would require a substantial increase in the plasmid copy number beyond the experimental evidence for pTAD-R. Consequently, we opted to use the regular replication mode for further simulations of pTAD-R allele dynamics. The correspondence between the simulation results under the regular replication mode and the allele dynamics in the evolution experiment suggests that the replication of plasmid pTAD-R is not random and may follow a regular mode of replication.

Beneficial alleles on polyploid replicons are prone to loss due to segregational drift To examine the effect of replicon polyploidy on the dynamics of beneficial alleles, we experimentally evolved the bacterial populations under conditions where the novel allele is beneficial, i.e., under selection for kanamycin resistance. In order to identify a comparable selection regime for both replicon types, we quantified the fitness effect of the nptII allele when encoded in the multicopy plasmid or single copy chromosome, employing increasing kanamycin concentrations (S5 and S6 Figs). During the evolution experiment, we applied a constant antibiotic pressure, which confers cells that are homozygous for the novel allele a fitness benefit of ca. 10% (i.e., selection coefficient s = 0.1185 ± 0.0415 SD for the novel resistant strain). The concentration of antibiotics employed to achieve a fitness benefit of ca. 10% in multicopy and single-copy replicon in novel homozygotes is similar. Therefore, we assume that the effect of gene dosage in our system is negligible; that is having more than one copy (1 to 15 copies) does not produce a higher fitness benefit to the bacteria at the concentrations employed. Note that the effect of selection is restricted to processes at the population level; antibiotic selection has no known effect on plasmid replication and segregation. We examined the effect of selection on the dynamics of the novel allele in combination with the effect of initial frequency and population bottleneck size. For each replicon type, we thus tested a total of 12 parameter combinations, each with 6 replicates (Figs 2, S3 and S4). PPT PowerPoint slide

PNG larger image

TIFF original image Download: Fig 2. Plasmid and chromosomal allele dynamics under selective and non-selective conditions. (A) The dynamics of a plasmid allele over approximately 200 generations. Points correspond to the proportion of kanamycin-resistant cells in each population and sampling time. Stacked bar plots present the state of the kanamycin-resistant populations at the end of the experiment. The fixation threshold is set at a frequency of 0.9. (B) Plasmid allele dynamics of representative populations from experimental conditions 2, 6 and 10 evolved under selective conditions and weak bottlenecks. The white arrow marks the point at which the frequencies of heterozygous and novel homozygous cells are equal. The remaining experimental populations are shown in S3 and S4 Figs). (C) The dynamics of a chromosomal allele over approximately 200 generations. https://doi.org/10.1371/journal.pgen.1010829.g002 We first focus on the results for the plasmid-encoded allele. As expected, the novel allele is maintained for longer and lost less often under selection than under neutral conditions. With high initial frequencies (f 0 = 10−4), the fraction of kanamycin-resistant cells increased in the population over time, regardless of the population bottleneck (Fig 2A, plots 1,3). The frequency of cells that carry the novel allele at the end of the experiment is close to the fixation threshold (that was set to 90%) in most replicates. In contrast, at low or moderate initial frequencies, the novel plasmid allele was quickly lost in most populations (Fig 2A). In populations that evolved under weak bottlenecks, the novel plasmid allele is lost in 1 out of 6 replicates for moderate initial frequencies and in 4 out of 6 replicates for low initial frequencies (Fig 2A, plots 5,9). In some replicate populations where the allele is maintained, we observed a slow decrease in the fraction of kanamycin-resistant cells from the 20th transfer on. Under the strong bottleneck, none of the replicate populations maintained the novel allele for moderate or low initial frequencies (Fig 2A, plots 7,11). The combination of selection and initial frequency has an effect also on the proportion of of homozygous and heterozygous cells within the resistant sub-population. For high initial frequency (f 0 = 10−4) under selective conditions, the frequency of cells that are homozygous for the novel allele increases rapidly and reaches a similar frequency to that of the heterozygotes within ca. 5–10 transfers (Fig 2B, plots 2.2 and 2.6). This initial rise in frequency is faster in comparison to that observed under non-selective conditions (ca. 10–15 transfers; Fig 1B). Under selection, both resistant genotypes, heterozygotes and homozygotes, are maintained in the population for a longer period of time and the frequency of heterozygotes decreases more slowly than under non-selective conditions. Nonetheless, heterozygotes are eventually lost from the population in several replicates (from ca. 25th onwards) (Figs 2B, S3 and S4). With medium and low initial frequencies, the homozygotes’ frequency increases more slowly, such that similar frequencies of homozygotes and heterozygotes are observed later in the experiment (Fig 2B). In summary, under selective conditions, both resistant genotypes are maintained for a longer time, where the relative frequencies of homozygotes and heterozygotes change as a consequence of segregational drift. To further examine the effect of polyploidy on the dynamics of novel alleles, we compared the dynamics observed for the plasmid-encoded allele to those observed for the chromosome-encoded allele. When the allele is encoded on the plasmid, the fraction of kanamycin-resistant cells increases more slowly and the time to fixation is longer in comparison to the chromosome-encoded allele (Fig 2C). The slower increase of the frequency of kanamycin-resistant cells in the experiments with the plasmid-encoded allele is best observed at mid-experiment (ca. transfer 15) for the high initial frequency (f 0 = 10−4), where the frequency of kanamycin-resistant cells is 100 times lower in comparison to the chromosome experiment. At the end of the experiment, the kanamycin resistance phenotype was fixed in a single population within the plasmid experiment, in contrast to the chromosome experiment where the resistance phenotype was fixed in the majority of populations (Fig 2A and 2C, plots 1,3). At lower initial frequencies (f 0 = 10−6, f 0 = 10−7), the differences between the frequency of kanamycin-resistant cells observed in the plasmid and chromosome experiments are even more pronounced, with up to ca. 1000-fold in favour of the chromosome populations (Fig 2A and 2C, plot 5). The most striking difference between plasmid and chromosome experiments is observed in populations that evolved under a strong bottleneck (b = 10−3). In those populations, the plasmid-encoded allele is lost in all replicate populations (Fig 2A and 2C, plots 7,11), in contrast to the chromosome experiment where the allele is maintained. To conclude, the opposing effects of segregational drift and selection lead to longer fixation time and loss of plasmid alleles in comparison to chromosomal alleles. Our results demonstrate that segregational drift has a pivotal role in the evolution of polyploid replicons. An alternative explanation for the observed loss of the novel beneficial plasmid allele could be the adaptation of the originally kanamycin-sensitive cells to the antibiotic environment, such that the novel allele is no longer beneficial. To rule out this alternative, we performed a competition experiment where the ancestral kanamycin-resistant populations competed against the evolved kanamycin-sensitive populations (S7 Fig). Our results show that the evolved kanamycin-sensitive populations still had a lower fitness under selective conditions compared to the ancestral kanamycin-resistant populations. Consequently, we conclude that the dynamics of the novel plasmid allele observed during the evolution experiment are the result of plasmid allele segregation rather than adaptation of the kanamycin-sensitive population to kanamycin.

The negative effect of segregational drift on the fixation of the beneficial allele is predicted by model simulations The experimental results in the previous sections are in line with previous theoretical predictions on the effects of replicon polyploidy on the dynamics of beneficial alleles [27]. For a direct comparison of the experimental results to theoretical predictions, we performed simulations with the model that is adapted to the experimental set-up. In our simulations, we assume that cells carrying only the ancestral allele replicate at a rate 1–s, and cells carrying at least one copy of the novel allele replicate at a rate 1. The model assumes no effect of gene dosage. To obtain a selection parameter that optimally reflects the allele dynamics over the entire duration of the experiment, we first estimated the parameter s from the mean frequencies of the kanamycin-resistant subpopulation in the evolution experiment for large initial frequencies (f 0 = 10−4) and weak bottlenecks (b = 10−2). The s parameters yielding an optimal fit between the simulated dynamics of the novel allele and the corresponding dynamics were selected for further simulations (S8 Fig). Those deviated from the experimental assessment of the fitness benefit of the novel allele: s = 0.087 ± 0.016 for the polyploid replicon (n rc = 15) and s = 0.136 ± 0.010 for the monoploid replicon (n rc = 1) (see S8 Fig). To examine how well our model captures the experimental results, we first performed simulations that correspond to the experimental conditions (Fig 3). The dynamics of the novel phenotype (kanamycin resistance) were overall in agreement with the results of the evolution experiment for both replicon types, albeit the frequency of the novel phenotype is slightly overestimated by the model given low or moderate initial frequencies (f 0 = 10−6 or 10−7) or strong bottlenecks, b = 10−3. The simulation results for the polyploid replicon show that the rise in frequency of heterozygotes in the experiment is well captured by the model simulations for high and moderate initial frequencies (f 0 = 10−4 and 10−6) and bottleneck size combinations. Nonetheless, the model underestimates the frequency of heterozygotes observed in the evolution experiment, which is consistent with the overestimation of the loss rate of heterozygotes under neutral conditions. The frequency of heterozygotes in the simulation starts decaying at nearly the same transfer as in the experiment, yet heterozygotes in the evolution experiment reach frequencies that are 10- to 100-fold higher in comparison to the simulations. Note that a comparison of the simulated allele dynamics to the plasmid experiment for the low initial frequency (f 0 = 10−7) is not possible, as most replicates in those experiments went extinct. Using the model simulation for the low initial frequency thus supplies a theoretical expectation for the beneficial allele dynamics in rare cases where the novel allele may persist in the population. The simulation results for moderate and low initial frequencies (f 0 = 10−6 and 10−7) predict a higher variability in the kanamycin-resistant cell frequencies per transfer for the polyploid replicon in comparison to the monoploid replicon (Fig 3A, plots 5–12). This prediction conforms to the results of the evolution experiments (S9 Fig). To further test the robustness of the comparison between theory and simulations regarding the choice of the parameter s in the model, we further performed computer simulations with s = 0.1. The resulting dynamics were qualitatively in agreement with our previous results as well (S10 Fig). We further compared the survival probability of the novel allele over time with the experimental results for the fraction of populations in which the allele is still present after a given transfer (Fig 3B). The model showed a remarkable agreement with the experimental results, capturing the main trends and often even making quantitatively accurate predictions. PPT PowerPoint slide

PNG larger image

TIFF original image Download: Fig 3. Model simulations of beneficial allele dynamics for polyploid (n rc = 15) and monoploid replicons. Predicted frequencies of the novel phenotype, given has not (yet) gone extinct, are shown in (A) and the corresponding survival of the novel allele are shown in (B). (A) Cell frequencies of the novel phenotype (parallel to kanamycin-resistant cells) (black) and heterozygous cells (green) at the end of each population transfer obtained from 100 simulations for a replicon with copy number n rc = 15 (1,3,5,7,9,11, for the multicopy plasmid) and with n rc = 1 (2,4,6,8,10,12, reflecting the chromosome). Lines show the mean of positive frequencies over the number of transfers, i.e., the average proportion of cells having the novel phenotype from all simulation trajectories, where the novel allele has not (yet) gone extinct (cf. (B) for the survival of the allele). The initial population at the start of the first transfer consists of ancestral cells and a small proportion, f, of cells with one novel replicon copy, where f 0 = 10−4 in (1–4), f 0 = 10−6 in (5–8), and f 0 = 10−7 in (8–12). The grey (dark grey) areas show the 99% (67%) confidence intervals of positive frequencies of the novel phenotype. The green lines show the mean frequency of heterozygous cells carrying both the ancestral and novel plasmid allele (conditioned on positive frequencies). Grey and green markers show the frequencies of kanamycin-resistant cells and heterozygous cells, respectively, of the corresponding experimental replicates. (B) Survival of the novel allele over number of transfers in model simulations (lines) and in the evolution experiment (markers) f 0 = 10−7. The survival of the novel allele initially present at high initial frequencies, f 0 = 10−4, are not shown since it is always 100%. The fraction of kanamycin-resistant cells of 100 simulation trajectories and 6 experimental replicates, respectively, is shown for each parameter combination. Simulation results were obtained by integrating Eq (1) with s = 0.087 for the polyploid replicon and s = 0.136 for the monoploid replicon until the population size has reached N c (t) cells in each growth phase (transfer), where N c (t) is taken from the mean experimental carrying capacity of each transfer t. A population bottleneck is applied after each growth phase, where a proportion of b = 10−2 or b = 10−3 cells are randomly selected from a multinomial distribution defined by the genotype frequencies at the end of the growth phase. https://doi.org/10.1371/journal.pgen.1010829.g003

[END]
---
[1] Url: https://journals.plos.org/plosgenetics/article?id=10.1371/journal.pgen.1010829

Published and (C) by PLOS One
Content appears here under this condition or license: Creative Commons - Attribution BY 4.0.

via Magical.Fish Gopher News Feeds:
gopher://magical.fish/1/feeds/news/plosone/

You are viewing proxied material from magical.fish. The copyright of proxied material belongs to its original authors. Any comments or complaints in relation to proxied material should be directed to the original authors of the content concerned. Please see the disclaimer for more details.