Phylogenetic analysis suggests single and multiple origins of dihydrofolate reductase mutations in Plasmodium vivax

A computational method for predicting the most likely evolutionary trajectories in the stepwise accumulation of resistance mutations

Pathogen evolution of drug resistance often occurs in a stepwise manner via the accumulation of multiple mutations that in combination have a non-additive impact on fitness, a phenomenon known as epistasis. The evolution of resistance via the accumulation of point mutations in the DHFR genes of Plasmodium falciparum (Pf) and Plasmodium vivax (Pv) has been studied extensively and multiple studies have shown epistatic interactions between these mutations determine the accessible evolutionary trajectories to highly resistant multiple mutations. Here, we simulated these evolutionary trajectories using a model of molecular evolution, parameterised using Rosetta Flex ddG predictions, where selection acts to reduce the target-drug binding affinity. We observe strong agreement with pathways determined using experimentally measured IC50 values of pyrimethamine binding, which suggests binding affinity is strongly predictive of resistance and epistasis in binding affinity strongly influences the order of fixation of resistance mutations. We also infer pathways directly from the frequency of mutations found in isolate data, and observe remarkable agreement with the most likely pathways predicted by our mechanistic model, as well as those determined experimentally. This suggests mutation frequency data can be used to intuitively infer evolutionary pathways, provided sufficient sampling of the population.

Genetic characterisation of the Theileria annulata cytochrome b locus and its impact on buparvaquone resistance in bovine

Control of tropical theileriosis, caused by the apicomplexan Theileria annulata, depends on the use of a single drug, buparvaquone, the efficacy of which is compromised by the emergence of resistance. The present study was undertaken to improve understanding of the role of mutations conferring buparvaquone resistance in T. annulata, and the effects of selection pressures on their emergence and spread. First, we investigated genetic characteristics of the cytochrome b locus associated with buparvaquone resistance in 10 susceptible and 7 resistant T. annulata isolates. The 129G (GGC) mutation was found in the Q₀₁ binding pocket and 253S (TCT) and 262S (TCA) mutations were identified within the Q₀₂ binding pocket. Next, we examined field isolates and identified cytochrome b mutations 129G (GGC), 253S (TCT) and 262S (TCA) in 21/75 buffalo-derived and 19/119 cattle-derived T. annulata isolates, providing evidence of positive selection pressure. Both hard and soft selective sweeps were identified, with striking differences between isolates. For example, 19 buffalo-derived and 7 cattle-derived isolates contained 129G (GGC) and 253S (TCT) resistance haplotypes at a high frequency, implying the emergence of resistance by a single mutation. Two buffalo-derived and 12 cattle-derived isolates contained equally high frequencies of 129G (GGC), 253S (TCT), 129G (GGC)/253S (TCT) and 262S (TCA) resistance haplotypes, implying the emergence of resistance by pre-existing or recurrent mutations. Phylogenetic analysis further revealed that 9 and 21 unique haplotypes in buffalo and cattle-derived isolates were present in a single lineage, suggesting a single origin. We propose that animal migration between farms is an important factor in the spread of buparvaquone resistance in endemic regions of Pakistan. The overall outcomes will be useful in understanding how drug resistance emerges and spreads, and this information will help design strategies to optimise the use and lifespan of the single most drug use to control tropical theileriosis.