Standing genetic variation and the evolution of drug resistance in HIV

Drug resistance remains a major problem for the treatment of HIV. Resistance can occur due to mutations that were present before treatment starts or due to mutations that occur during treatment. The relative importance of these two sources is unknown. We study three different situations in which HIV drug resistance may evolve: starting triple-drug therapy, treatment with a single dose of nevirapine and interruption of treatment. For each of these three cases good data are available from literature, which allows us to estimate the probability that resistance evolves from standing genetic variation. Depending on the treatment we find probabilities of the evolution of drug resistance due to standing genetic variation between 0 and 39%. For patients who start triple-drug combination therapy, we find that drug resistance evolves from standing genetic variation in approximately 6% of the patients. We use a population-dynamic and population-genetic model to understand the observations and to estimate important evolutionary parameters. We find that both, the effective population size of the virus before treatment, and the fitness of the resistant mutant during treatment, are key-parameters that determine the probability that resistance evolves from standing genetic variation. Importantly, clinical data indicate that both of these parameters can be manipulated by the kind of treatment that is used.Comment: 33 pages 6 figure

The relationship between the length of a treatment interruption and the probability that DRMs become established.

<p>Estimated probability that resistance mutations become established due to a single treatment interruption. Grey bars are data from seven clinical trials, estimated standard error (see supplementary table S3 in text S2). The number of patients (and the number of interruptions per patient) are noted at the top of the graph. The red circles are estimated from 1000 simulations, estimated standard error. The black squares are predictions using the average population size from the simulations and equation 2. Parameters as in <a href="http://www.ploscompbiol.org/article/info:doi/10.1371/journal.pcbi.1002527#pcbi-1002527-t002" target="_blank">table 2</a>.</p

Drug level and population size during and after a treatment interruption.

<p>Drug level (dashed line) and viral population size (solid line) during and after a treatment interruption. Red bars indicate when drugs are taken.</p

Abbreviations.

<p>Abbreviations.</p

Probability of the establishment of DRMs as a function of effective population size and the fitness of the resistant mutant during treatment.

<p>The predicted probability of the establishment of drug resistance mutations from standing genetic variation depending on the effective population size and the fitness of the resistant mutant during therapy. Grey scales indicate the probability of the evolution of drug resistance due to standing genetic variation. Dots indicate estimated parameter combinations for treatment with just sdNVP, with ZDV monotherapy followed by sdNVP (ZDV/sdNVP), with sdNVP followed by two additional drugs postpartum (sdNVP/PP) and with ZDV monotherapy followed by sdNVP and two additional drugs postpartum ZDV/sdNVP/PP.</p

Possible combinations of population size and fitness and the effect of population sizes on the probability that DRMs are present before treatment.

<p><b><a href="http://www.ploscompbiol.org/article/info:doi/10.1371/journal.pcbi.1002527#pcbi-1002527-g002" target="_blank">Figure 2a:</a></b> Continuous line: combinations of population size before treatment () and fitness of mutant virus during therapy () that lead to the observed probability that resistance mutations from standing genetic variation become established (). Dashed line: combinations of population size during treatment () and fitness of mutant virus during therapy () that lead to the observed probability that resistance mutations from standing genetic variation become established (). Open dot: and , closed dot: , . <a href="http://www.ploscompbiol.org/article/info:doi/10.1371/journal.pcbi.1002527#pcbi-1002527-g002" target="_blank"><b>Figure 2b</b>:</a> Probability that a patient has any pre-existing DRMs before the start of therapy for different population sizes, and . Open dot: .</p

Probability of detecting resistance per year of treatment.

<p>The probability that resistance is detected for the first time in the first, second or third year of treatment, given that it was not detected until then. Grey bars are the estimates from the Margot et al (<a href="http://www.ploscompbiol.org/article/info:doi/10.1371/journal.pcbi.1002527#pcbi.1002527-Margot1" target="_blank">[36]</a>) dataset, and the number of patients on which the estimates are based are noted at the top of the graph. The red dashed area reflects the inferred probability that resistance mutations from standing genetic variation become established. The black squares are values calculated using equations 2 and 4. The red circles are estimated from 1000 simulations. Parameters as in <a href="http://www.ploscompbiol.org/article/info:doi/10.1371/journal.pcbi.1002527#pcbi-1002527-t002" target="_blank">table 2</a>.</p