Selection strength and hitchhiking around two anti-malarial resistance genes

Approximations for the hitchhiking effect caused by the evolution of antimalarial-drug resistance

An analytically feasible, deterministic model for the spread of drug resistance among human malaria parasites, which incorporates all characteristics of the complex malaria-transmission cycle was introduced by Schneider and Kim (Theor. Popul Biol, 2010). The model accounts for the fact that only a fraction of infected hosts receive drug treatment and that hosts can be co-infected by differently many parasites. Furthermore, the model also incorporates host heterogeneity. Antimalarial-drug resistance is assumed to be caused by a single locus with two alleles-a sensitive one and a resistance one. The most important result for this model is that an analytical solution for the frequencies of a linked neutral biallelic locus exists. However, the exact solution does not admit an explicit form, and cannot straightforwardly be interpreted in terms of the model parameters. Here, we establish simple approximations for the equilibrium frequency at the neutral locus. Under the assumption that the resistant allele is initially rare-the biologically most relevant assumption in this context-and that recombination is weak, the approximations become similar to the approximations in the standard hitchhiking model. However, there are crucial differences. In particular, because of the high degree of selfing among malaria parasites in their sexual phase, a genome-wide reduction of relative heterozygosity occurs if selection is sufficiently strong. It turns out that the approximations are accurate even if the recombination rates are not small and the resistant allele is initially not very rare. The main advantage of our approximations is that they are easy to interpret in terms of model parameters. Moreover, they allow to make predictions of the size of the valley of reduced heterozygosity around the selected locus for given model parameters. Reversely, for a given reduction of heterozygosity, it is possible to identify the corresponding parameters. Moreover, we will show that incorporating host heterogeneity leads to an increased hitchhiking effect.

Identification of local- and habitat-dependent selection: scanning functionally important genes in nine-spined sticklebacks (Pungitius pungitius)

Understanding the selective forces promoting adaptive population divergence is a central issue in evolutionary biology. The role of environmental salinity in driving adaptation and evolution in aquatic organisms is still poorly understood. We investigated the relative impacts of habitat type (cf. saltwater vs. freshwater) and geographic area in shaping adaptive population divergence, as well as genes responsible for adaptation to different salinities in nine-spined sticklebacks (Pungitius pungitius). To this end, we employed a hitchhiking mapping approach with 111 microsatellite loci and one insertion/deletion locus including 63 loci situated within or close to genes with important physiological functions such as osmoregulation, growth, and thermal response. Using three pairs of marine and freshwater populations from different geographic areas, we identified several loci showing consistent evidence of being under directional selection in different outlier tests. Analyses of molecular variance at the loci under selection indicated that geographic area rather than habitat type has been acting as a central force in shaping adaptive population divergence. Nevertheless, both outlier tests and a spatial analysis method indicated that two loci (growth hormone receptor 2 and DEAD box polypeptide 56) are involved in adaptation to different habitats, implying that environmental salinity has been affecting them as a selective force. These loci are promising candidates for further investigations focusing on the molecular mechanisms of adaptation to marine and freshwater environments.

An analytical model for genetic hitchhiking in the evolution of antimalarial drug resistance

We analytically study a deterministic model for the spread of drug resistance among human malaria parasites. The model incorporates all major characteristics of the complex malaria transmission cycle and accounts for the fact that only a fraction alpha of infected hosts receive drug treatment. Furthermore, the model incorporates that hosts can be co-infected. The number m of parasites co-infecting a host is either a constant or, more generally, follows a given frequency distribution. Although the model is formulated in a multilocus setup, for our results we assume that drug resistance is caused by a single locus with two alleles - a sensitive one and a resistant one. We assume that the resistant allele has a selective advantage only in treated hosts and pays metabolic costs, which causes this allele to be deleterious in untreated hosts. We provide necessary and sufficient conditions for the fixation of the resistant allele. Moreover, provided the resistant allele will sweep through the population, we derive a formula for the time until it reaches a given frequency and in particular for the time until quasi-fixation. Furthermore, we establish an analytical solution for allele frequency changes at a linked neutral biallelic locus due to the rapid increase in frequency of the resistant allele. Our solution describes a local reduction in heterozygosity among parasite chromosomes around the resistant allele, the effect commonly referred to as the hitchhiking effect, as a function of alpha and m. The result therefore allows the investigation of selective sweep patterns under specific demographic settings. We find that the hitchhiking effect is similar but different from the standard model of genetic hitchhiking that assumes random mating and homogeneous selection. In particular, the process of recombination and selection cannot be decoupled. We further explain why standard hitchhiking theory cannot be applied to drug resistance in malaria. Furthermore, we will show that a genome-wide reduction in relative heterozygosity can occur provided a fraction of hosts is infected by a single parasite haplotype. Finally, we show how to incorporate host heterogeneity, and generalize our results to this biologically more realistic case.

Origins and spread of pfdhfr mutant alleles in Plasmodium falciparum

The emergence and spread of Plasmodium falciparum parasite resistant to sulfadoxine and pyrimethamine (SP) poses a serious public health problem. Resistance is caused by point mutations in dihydrofolate reductase (pfdhfr) and dihydropteroate synthase (pfdhps), the two key enzymes in the folate biosynthetic pathway. The use of microsatellite markers flanking pfdhfr has recently shown that the invasion of limited resistant lineages may explain the widespread SP resistance in many endemic regions. In Africa, however, multiple indigenous origins of pfdhfr triple mutants have been demonstrated. More new independent lineages and routes of geographical spread of resistance may be found by further molecular evolutionary analyses using samples from various endemic regions. Here, I review recent studies about the history of SP usage and the evolution and spread of resistant lineages while addressing the technical issue of microsatellite analysis.

Detecting signatures of balancing selection to identify targets of anti-parasite immunity

Parasite antigen genes might evolve under frequency-dependent immune selection. The distinctive patterns of polymorphism that result can be detected using population genetic methods that test for signatures of balancing selection, allowing genes encoding important targets of immunity to be identified. Analyses can be complicated by population structures, histories and features of a parasite's genome. However, new sequencing technologies facilitate scans of polymorphism throughout parasite genomes to identify the most exceptional gene specific signatures. We focus on malaria parasites to illustrate challenges and opportunities for detecting targets of frequency-dependent immune selection to discover new potential vaccine candidates.

Darwin and genetics

Darwin's theory of natural selection lacked an adequate account of inheritance, making it logically incomplete. We review the interaction between evolution and genetics, showing how, unlike Mendel, Darwin's lack of a model of the mechanism of inheritance left him unable to interpret his own data that showed Mendelian ratios, even though he shared with Mendel a more mathematical and probabilistic outlook than most biologists of his time. Darwin's own "pangenesis" model provided a mechanism for generating ample variability on which selection could act. It involved, however, the inheritance of characters acquired during an organism's life, which Darwin himself knew could not explain some evolutionary situations. Once the particulate basis of genetics was understood, it was seen to allow variation to be passed intact to new generations, and evolution could then be understood as a process of changes in the frequencies of stable variants. Evolutionary genetics subsequently developed as a central part of biology. Darwinian principles now play a greater role in biology than ever before, which we illustrate with some examples of studies of natural selection that use DNA sequence data and with some recent advances in answering questions first asked by Darwin.

Origin and evolution of sulfadoxine resistant Plasmodium falciparum

The Thailand-Cambodia border is the epicenter for drug-resistant falciparum malaria. Previous studies have shown that chloroquine (CQ) and pyrimethamine resistance originated in this region and eventually spread to other Asian countries and Africa. However, there is a dearth in understanding the origin and evolution of dhps alleles associated with sulfadoxine resistance. The present study was designed to reveal the origin(s) of sulfadoxine resistance in Cambodia and its evolutionary relationship to African and South American dhps alleles. We sequenced 234 Cambodian Plasmodium falciparum isolates for the dhps codons S436A/F, A437G, K540E, A581G and A613S/T implicated in sulfadoxine resistance. We also genotyped 10 microsatellite loci around dhps to determine the genetic backgrounds of various alleles and compared them with the backgrounds of alleles prevalent in Africa and South America. In addition to previously known highly-resistant triple mutant dhps alleles SGEGA and AGEAA (codons 436, 437, 540, 581, 613 are sequentially indicated), a large proportion of the isolates (19.3%) contained a 540N mutation in association with 437G/581G yielding a previously unreported triple mutant allele, SGNGA. Microsatellite data strongly suggest the strength of selection was greater on triple mutant dhps alleles followed by the double and single mutants. We provide evidence for at least three independent origins for the double mutants, one each for the SGKGA, AGKAA and SGEAA alleles. Our data suggest that the triple mutant allele SGEGA and the novel allele SGNGA have common origin on the SGKGA background, whereas the AGEAA triple mutant was derived from AGKAA on multiple, albeit limited, genetic backgrounds. The SGEAA did not share haplotypes with any of the triple mutants. Comparative analysis of the microsatellite haplotypes flanking dhps alleles from Cambodia, Kenya, Cameroon and Venezuela revealed an independent origin of sulfadoxine resistant alleles in each of these regions.

Widespread resistance to chloroquine (CQ) and sulfadoxine-pyrimethamine (SP), the two least expensive and widely available antimalarial drugs, has become a major global public health challenge. It is known that point mutations in Plasmodium falciparum crt, dhfr and dhps genes contribute to resistance to CQ, pyrimethamine and sulfadoxine, respectively. CQ and pyrimethamine resistance spread to Africa and Asia from a few founding mutant lineages originating from the Thailand-Cambodia border. Here, we define the origins of dhps alleles in Cambodia and their relationships to African and South American counterparts. Three different triple mutant alleles including a novel allele comprised of 437G, 540N, and 581G mutations (S₄₃₆G₄₃₇N₅₄₀G₅₈₁A₆₁₃) were found in Cambodia as opposed to a single triple mutant allele in South America and two common double mutant alleles in Africa. Microsatellite data suggest strong selection operating on triple mutant alleles as compared to double and single mutants in Cambodia. We report three major independent origins for the double mutants and at least two independent origins for the highly resistant triple mutant dhps alleles in Cambodia. We also show that the resistant dhps alleles in Africa and South America have distinct origins from Cambodia. These results suggest that the evolution and spread of sulfadoxine resistance is different from CQ and pyrimethamine resistance.

Evidence of selective sweeps in genes conferring resistance to chloroquine and pyrimethamine in Plasmodium falciparum isolates in India

Treatment of Plasmodium falciparum is complicated by the emergence and spread of parasite resistance to many of the first-line drugs used to treat malaria. Antimalarial drug resistance has been associated with specific point mutations in several genes, suggesting that these single nucleotide polymorphisms can be useful in tracking the emergence of drug resistance. In India, P. falciparum infection can manifest itself as asymptomatic, mild, or severe malaria, with or without cerebral involvement. We tested whether chloroquine- and antifolate drug-resistant genotypes would be more commonly associated with cases of cerebral malaria than with cases of mild malaria in the province of Jabalpur, India, by genotyping the dhps, dhfr, pfmdr-1, and pfcrt genes using pyrosequencing, direct sequencing, and real-time PCR. Further, we used microsatellites surrounding the genes to determine the origins and spread of the drug-resistant genotypes in this area. Resistance to chloroquine was essentially fixed, with 95% of the isolates harboring the pfcrt K76T mutation. Resistant genotypes of dhfr, dhps, and pfmdr-1 were found in 94%, 17%, and 77% of the isolates, respectively. Drug-resistant genotypes were equally likely to be associated with cerebral malaria as with mild malaria. We found evidence of a selective sweep in pfcrt and, to a lesser degree, in dhfr, indicating high levels of resistance to chloroquine and evolving resistance to pyrimethamine. Microsatellites surrounding pfcrt indicate that the resistant genotypes (SVMNT) were most similar to those found in Papua New Guinea.

pfmdr1 amplification and fixation of pfcrt chloroquine resistance alleles in Plasmodium falciparum in Venezuela

Molecular tools are valuable for determining evolutionary history and the prevalence of drug-resistant malaria parasites. These tools have helped to predict decreased sensitivity to antimalarials and fixation of multidrug resistance genotypes in some regions. In order to assess how historical drug policies impacted Plasmodium falciparum in Venezuela, we examined molecular changes in genes associated with drug resistance. We examined pfmdr1 and pfcrt in samples from Sifontes, Venezuela, and integrated our findings with earlier work describing dhfr and dhps in these samples. We characterized pfmdr1 genotypes and copy number variation, pfcrt genotypes, and proximal microsatellites in 93 samples originating from surveillance from 2003 to 2004. Multicopy pfmdr1 was found in 12% of the samples. Two pfmdr1 alleles, Y184F/N1042D/D1246Y (37%) and Y184F/S1034C/N1042D/D1246Y (63%), were found. These alleles share ancestry, and no evidence of strong selective pressure on mutations was found. pfcrt chloroquine resistance alleles are fixed with two alleles: S(tct)VMNT (91%) and S(agt)VMNT (9%). These alleles are associated with strong selection. There was also an association between pfcrt, pfmdr1, dhfr, and dhps genotypes/haplotypes. Duplication of pfmdr1 suggests a potential shift in mefloquine sensitivity in this region, which warrants further study. A bottleneck occurred in P. falciparum in Sifontes, Venezuela, and multidrug resistance genotypes are present. This population could be targeted for malaria elimination programs to prevent the possible spread of multidrug-resistant parasites.

The dynamics of mutations associated with anti-malarial drug resistance in Plasmodium falciparum

The evolution of resistance in Plasmodium falciparum against safe and affordable drugs such as chloroquine (CQ) and sulfadoxine-pyrimethamine (SP) is a major global health threat. Investigating the dynamics of resistance against these antimalarial drugs will lead to approaches for addressing the problem of resistance in malarial parasites that are solidly based in evolutionary genetics and population biology. In this article, we discuss current developments in population biology modeling and evolutionary genetics. Despite great advancements achieved in the past decade, understanding the complex dynamics of mutations conferring drug resistance in P. falciparum requires approaches that consider the parasite population structure among other demographic processes.

Characteristics of genetic hitchhiking around dihydrofolate reductase gene associated with pyrimethamine resistance in Plasmodium falciparum isolates from India

Sulfadoxine-pyrimethamine (SP) resistance in Plasmodium falciparum has been widespread across continents, causing the major hurdle of controlling malaria. Resistance is encoded mainly by point mutations in P. falciparum dihydrofolate reductase (pfdhfr) and dihydropteroate synthase (pfdhps) target genes. To study the origin and evolution of pyrimethamine resistance on the Indian subcontinent, microsatellite markers flanking the pfdhfr gene were mapped. Here we describe the characteristics of genetic hitchhiking around the pfdhfr gene among 190 P. falciparum isolates. These isolates were collected from five different geographical regions of India (Uttar Pradesh, Madhya Pradesh, Assam, Orissa, and Andaman and Nicobar Islands) where malarial transmission rates and levels of drug resistance vary across regions. Among the isolates, we observed a significant reduction in genetic variation in the +/-20-kb vicinity of the mutant pfdhfr alleles due to hitchhiking. This reduction in genetic diversity was more prominent around quadruple pfdhfr alleles (heterozygosity [H(e)] = 0.23) than around double (H(e) = 0.365) and single (H(e) = 0.465) mutant alleles. Asymmetry in the selective sweep flanking the pfdhfr alleles was observed with regional isolates, emphasizing the drug usage with the parasite population. All the pfdhfr alleles share a single microsatellite haplotype and seem to have originated from a single progenitor similar to that of Southeast Asian (Thailand) pfdhfr mutants. Results of the present study also indicate that the emergence of drug-resistant alleles is a recent phenomenon in India compared to Southeast Asian countries.

Impaired fitness of drug-resistant malaria parasites: evidence and implication on drug-deployment policies

Malaria, a leading parasitic disease, inflicts an enormous toll on human lives and is caused by protozoal parasites belonging to the genus Plasmodium. Antimalarial drugs targeting essential biochemical processes in the parasite are the primary resources for management and control. However, the parasite has established mutations, substantially reducing the efficacy of these drugs. First-line therapy is faced the with the consistent evolution of drug-resistant genotypes carrying these mutations. However, drug-resistant genotypes are likely to be less fit than the wild-type, suggesting that they might disappear by reducing the volume of drug pressure. A substantial body of epidemiological evidence confirmed that the frequency of resistant genotypes wanes when active drug selection declines. Drug selection on the parasite genome that removes genetic variation in the vicinity of drug-resistant genes (hitch-hiking) is common among resistant parasites in the field. This can further disadvantage drug-resistant strains and limit their variability in the face of a mounting immune response. Attempts to provide unequivocal evidence for the fitness cost of drug resistance have monitored the outcomes of laboratory competition experiments of deliberate mixtures of sensitive and resistant strains, in the absence of drug pressure, using isogenic clones produced either by drug selection or gene manipulation. Some of these experiments provided inconclusive results, but they all suggested reduced fitness of drug-resistant clones in the absence of drug pressure. In addition, biochemical analyses provided clearer information demonstrating that the mutation of some antimalarial-targeted enzymes lowers their activity compared with the wild-type enzyme. Here, we review current evidences for the disadvantage of drug-resistance mutations, and discuss some strategies of drug deployment to maximize the cost of resistance and limit its spread.

Spread and evolution of Plasmodium falciparum drug resistance

Worldwide spread of Plasmodium falciparum drug resistance to conventional antimalarials, chloroquine and sulfadoxine/pyrimethamine, has been imposing a serious public health problem in many endemic regions. Recent discovery of drug resistance-associated genes, pfcrt, pfmdr1, dhfr, and dhps, and applications of microsatellite markers flanking the genes have revealed the evolution of parasite resistance to these antimalarials and the geographical spread of drug resistance. Here, we review our recent knowledge of the evolution and spread of parasite resistance to chloroquine and sulfadoxine/pyrimethamine. In both antimalarials, resistance appears to be largely explained by the invasion of limited resistant lineages to many endemic regions. However, multiple, indigenous evolutionary origins of resistant lineages have also been demonstrated. Further molecular evolutionary and population genetic approaches will greatly facilitate our understanding of the evolution and spread of parasite drug resistance, and will contribute to developing strategies for better control of malaria.

Multiple origins and regional dispersal of resistant dhps in African Plasmodium falciparum malaria

BACKGROUND

Although the molecular basis of resistance to a number of common antimalarial drugs is well known, a geographic description of the emergence and dispersal of resistance mutations across Africa has not been attempted. To that end we have characterised the evolutionary origins of antifolate resistance mutations in the dihydropteroate synthase (dhps) gene and mapped their contemporary distribution.

METHODS AND FINDINGS

We used microsatellite polymorphism flanking the dhps gene to determine which resistance alleles shared common ancestry and found five major lineages each of which had a unique geographical distribution. The extent to which allelic lineages were shared among 20 African Plasmodium falciparum populations revealed five major geographical groupings. Resistance lineages were common to all sites within these regions. The most marked differentiation was between east and west African P. falciparum, in which resistance alleles were not only of different ancestry but also carried different resistance mutations.

CONCLUSIONS

Resistant dhps has emerged independently in multiple sites in Africa during the past 10-20 years. Our data show the molecular basis of resistance differs between east and west Africa, which is likely to translate into differing antifolate sensitivity. We have also demonstrated that the dispersal patterns of resistance lineages give unique insights into recent parasite migration patterns.

Dynamics of malaria drug resistance patterns in the Amazon basin region following changes in Peruvian national treatment policy for uncomplicated malaria

Monitoring changes in the frequencies of drug-resistant and -sensitive genotypes can facilitate in vivo clinical trials to assess the efficacy of drugs before complete failure occurs. Peru changed its national treatment policy for uncomplicated malaria to artesunate (ART)-plus-mefloquine (MQ) combination therapy in the Amazon basin in 2001. We genotyped isolates collected in 1999 and isolates collected in 2006 to 2007 for mutations in the Plasmodium falciparum dihydrofolate reductase (Pfdhfr) and dihydropteroate synthase (Pfdhps) genes, multidrug resistance gene 1 (Pfmdr-1), the chloroquine (CQ) resistance transporter gene (Pfcrt), and the Ca(2+) ATPase gene (PfATP6); these have been shown to be involved in resistance to sulfadoxine-pyrimethamine (SP), MQ, CQ, and possibly ART, respectively. Microsatellite haplotypes around the Pfdhfr, Pfdhps, Pfcrt, and Pfmdr-1 loci were also determined. There was a significant decline in the highly SP resistant Pfdhfr and Pfdhps genotypes from 1999 to 2006. In contrast, a CQ-resistant Pfcrt genotype increased in frequency during the same period. Among five different Pfmdr-1 allelic forms noted in 1999, two genotypes increased in frequency while one genotype decreased by 2006. We also noted previously undescribed polymorphisms in the PfATP6 gene as well as an increase in the frequency of a deletion mutant during this period. In addition, microsatellite analysis revealed that the resistant Pfdhfr, Pfdhps, and Pfcrt genotypes have each evolved from a single founder haplotype, while Pfmdr-1 genotypes have evolved from at least two independent haplotypes. Importantly, this study demonstrates that the Peruvian triple mutant Pfdhps genotypes are very similar to those found in other parts of South America.

Extensive heterozygosity in flanking microsatellites of Plasmodium falciparum Na+/H+ exchanger (pfnhe-1) gene among Indian isolates

Plasmodium falciparum Na(+)/H(+) exchanger-1 (pfnhe-1) gene has been proposed to be a possible marker for quinine resistance. Here, we describe the sequence analysis of the flanking microsatellites of the pfnhe-1 gene among 108 Indian P.falciparum isolates. Among the parasite population, a high degree of polymorphism was observed at all the 10 microsatellite loci within +/-40kb region of the pfnhe-1 gene where the number of alleles varied from 2 to 16 with a high expected heterozygosity ranging from 0.43 to 0.91 at these loci. Also, higher levels of heterozygosity have been observed in P.falciparum isolates collected from both low and high transmission and drug resistant areas. Furthermore, there was no association between QN resistance associated DNNND repeats in PFNHE-1 and the flanking microsatellite haplotypes. In conclusion, the observed high level of microsatellite polymorphism and absence of selective sweep in the flanking +/-40kb region of the pfnhe-1 gene could be an indication that there is no strong selection pressure on this target gene.

Genome-wide SNP genotyping highlights the role of natural selection in Plasmodium falciparum population divergence

BACKGROUND

The malaria parasite Plasmodium falciparum exhibits abundant genetic diversity, and this diversity is key to its success as a pathogen. Previous efforts to study genetic diversity in P. falciparum have begun to elucidate the demographic history of the species, as well as patterns of population structure and patterns of linkage disequilibrium within its genome. Such studies will be greatly enhanced by new genomic tools and recent large-scale efforts to map genomic variation. To that end, we have developed a high throughput single nucleotide polymorphism (SNP) genotyping platform for P. falciparum.

RESULTS

Using an Affymetrix 3,000 SNP assay array, we found roughly half the assays (1,638) yielded high quality, 100% accurate genotyping calls for both major and minor SNP alleles. Genotype data from 76 global isolates confirm significant genetic differentiation among continental populations and varying levels of SNP diversity and linkage disequilibrium according to geographic location and local epidemiological factors. We further discovered that nonsynonymous and silent (synonymous or noncoding) SNPs differ with respect to within-population diversity, inter-population differentiation, and the degree to which allele frequencies are correlated between populations.

CONCLUSIONS

The distinct population profile of nonsynonymous variants indicates that natural selection has a significant influence on genomic diversity in P. falciparum, and that many of these changes may reflect functional variants deserving of follow-up study. Our analysis demonstrates the potential for new high-throughput genotyping technologies to enhance studies of population structure, natural selection, and ultimately enable genome-wide association studies in P. falciparum to find genes underlying key phenotypic traits.

Adaptive copy number evolution in malaria parasites

Copy number polymorphism (CNP) is ubiquitous in eukaryotic genomes, but the degree to which this reflects the action of positive selection is poorly understood. The first gene in the Plasmodium folate biosynthesis pathway, GTP-cyclohydrolase I (gch1), shows extensive CNP. We provide compelling evidence that gch1 CNP is an adaptive consequence of selection by antifolate drugs, which target enzymes downstream in this pathway. (1) We compared gch1 CNP in parasites from Thailand (strong historical antifolate selection) with those from neighboring Laos (weak antifolate selection). Two percent of chromosomes had amplified copy number in Laos, while 72% carried multiple (2–11) copies in Thailand, and differentiation exceeded that observed at 73 synonymous SNPs. (2) We found five amplicon types containing one to greater than six genes and spanning 1 to >11 kb, consistent with parallel evolution and strong selection for this gene amplification. gch1 was the only gene occurring in all amplicons suggesting that this locus is the target of selection. (3) We observed reduced microsatellite variation and increased linkage disequilibrium (LD) in a 900-kb region flanking gch1 in parasites from Thailand, consistent with rapid recent spread of chromosomes carrying multiple copies of gch1. (4) We found that parasites bearing dhfr-164L, which causes high-level resistance to antifolate drugs, carry significantly (p = 0.00003) higher copy numbers of gch1 than parasites bearing 164I, indicating functional association between genes located on different chromosomes but linked in the same biochemical pathway. These results demonstrate that CNP at gch1 is adaptive and the associations with dhfr-164L strongly suggest a compensatory function. More generally, these data demonstrate how selection affects multiple enzymes in a single biochemical pathway, and suggest that investigation of structural variation may provide a fast-track to locating genes underlying adaptation.

Recent comparative genomic hybridization studies have revealed extensive copy number variation in eukaryotic genomes. The first gene in the Plasmodium folate biosynthesis pathway, GTP-cyclohydrolase I (gch1), shows extensive copy number polymorphism (CNP). We provide compelling evidence that gch1 CNP is adaptive and most likely results from selection by antifolate drugs, which target enzymes downstream in this pathway. Gch1 CNP shows extreme geographical differentiation; hitchhiking reduces diversity and increases LD in flanking sequence, indicating recent rapid spread within Thailand, while amplicon structure reveals multiple origins and parallel evolution. Furthermore, strong association between elevated copy number and a critical mutation dhfr-I164L that underlies high-level antifolate resistance indicates functional linkage and fitness epistasis between genes on different chromosomes. These data reveal hidden complexity in the evolutionary response to antifolate treatment and demonstrate that analysis of structural variation can provide a fast-track to locating genes that underlie adaptation.

Hitchhiking and selective sweeps of Plasmodium falciparum sulfadoxine and pyrimethamine resistance alleles in a population from central Africa

Sulfadoxine-pyrimethamine (SP) resistance in Plasmodium falciparum is encoded by a number of mutations in the dihydrofolate reductase (dhfr) and dihydropteroate synthetase (dhps) genes. Here, we have characterized point mutations in dhfr and dhps and microsatellite loci around dhfr on chromosome 4 and dhps on chromosome 8 as well as neutral markers on chromosomes 2 and 3 in 332 samples from Yaoundé, Cameroon. The triple mutant dhfr haplotype that originated in Southeast Asia is the most predominant in this sample set, but we also find additional independent haplotypes at low frequency and an incipient process of genetic differentiation among alleles of Southeast Asian origin. As reported for other African populations, we find evidence of a selective sweep for resistant dhfr mutants in this Cameroonian population due to drug selection. Although we find evidence for a selective sweep in dhps mutants associated with SP resistance, the dynamics of dhps mutants appear different than those observed for dhfr mutants. Overall, our results yield support for the use of microsatellite markers to track resistant parasites; however, the detection of resistant dhfr alleles in low frequency, the evidence of divergence among dhfr alleles that share a common evolutionary origin, and the distinct dynamics of resistant dhps alleles emphasize the importance of comprehensive, population-based investigations to evaluate the effects of drug selection on parasite populations.

Multiple origins of resistance-conferring mutations in Plasmodium vivax dihydrofolate reductase

Background

In order to maximize the useful therapeutic life of antimalarial drugs, it is crucial to understand the mechanisms by which parasites resistant to antimalarial drugs are selected and spread in natural populations. Recent work has demonstrated that pyrimethamine-resistance conferring mutations in Plasmodium falciparum dihydrofolate reductase (dhfr) have arisen rarely de novo, but spread widely in Asia and Africa. The origin and spread of mutations in Plasmodium vivax dhfr were assessed by constructing haplotypes based on sequencing dhfr and its flanking regions.

Methods

The P. vivax dhfr coding region, 792 bp upstream and 683 bp downstream were amplified and sequenced from 137 contemporary patient isolates from Colombia, India, Indonesia, Papua New Guinea, Sri Lanka, Thailand, and Vanuatu. A repeat motif located 2.6 kb upstream of dhfr was also sequenced from 75 of 137 patient isolates, and mutational relationships among the haplotypes were visualized using the programme Network.

Results

Synonymous and non-synonymous single nucleotide polymorphisms (SNPs) within the dhfr coding region were identified, as was the well-documented in-frame insertion/deletion (indel). SNPs were also identified upstream and downstream of dhfr, with an indel and a highly polymorphic repeat region identified upstream of dhfr. The regions flanking dhfr were highly variable. The double mutant (58R/117N) dhfr allele has evolved from several origins, because the 58R is encoded by at least 3 different codons. The triple (58R/61M/117T) and quadruple (57L/61M/117T/173F, 57I/58R/61M/117T and 57L/58R/61M/117T) mutant alleles had at least three independent origins in Thailand, Indonesia, and Papua New Guinea/Vanuatu.

Conclusion

It was found that the P. vivax dhfr coding region and its flanking intergenic regions are highly polymorphic and that mutations in P. vivax dhfr that confer antifolate resistance have arisen several times in the Asian region. This contrasts sharply with the selective sweep of rare antifolate resistant alleles observed in the P. falciparum populations in Asia and Africa. The finding of multiple origins of resistance-conferring mutations has important implications for drug policy.

Discordant patterns of genetic variation at two chloroquine resistance loci in worldwide populations of the malaria parasite Plasmodium falciparum

Mutations in the chloroquine resistance (CQR) transporter gene of Plasmodium falciparum (Pfcrt; chromosome 7) play a key role in CQR, while mutations in the multidrug resistance gene (Pfmdr1; chromosome 5) play a significant role in the parasite's resistance to a variety of antimalarials and also modulate CQR. To compare patterns of genetic variation at Pfcrt and Pfmdr1 loci, we investigated 460 blood samples from P. falciparum-infected patients from four Asian, three African, and three South American countries, analyzing microsatellite (MS) loci flanking Pfcrt (five loci [approximately 40 kb]) and Pfmdr1 (either two loci [approximately 5 kb] or four loci [approximately 10 kb]). CQR Pfmdr1 allele-associated MS haplotypes showed considerably higher genetic diversity and higher levels of subdivision than CQR Pfcrt allele-associated MS haplotypes in both Asian and African parasite populations. However, both Pfcrt and Pfmdr1 MS haplotypes showed similar levels of low diversity in South American parasite populations. Median-joining network analyses showed that the Pfcrt MS haplotypes correlated well with geography and CQR Pfcrt alleles, whereas there was no distinct Pfmdr1 MS haplotype that correlated with geography and/or CQR Pfmdr1 alleles. Furthermore, multiple independent origins of CQR Pfmdr1 alleles in Asia and Africa were inferred. These results suggest that variation at Pfcrt and Pfmdr1 loci in both Asian and African parasite populations is generated and/or maintained via substantially different mechanisms. Since Pfmdr1 mutations may be associated with resistance to artemisinin combination therapies that are replacing CQ, particularly in Africa, it is important to determine if, and how, the genetic characteristics of this locus change over time.

Local adaptation and vector-mediated population structure in Plasmodium vivax malaria

Plasmodium vivax in southern Mexico exhibits different infectivities to 2 local mosquito vectors, Anopheles pseudopunctipennis and Anopheles albimanus. Previous work has tied these differences in mosquito infectivity to variation in the central repeat motif of the malaria parasite's circumsporozoite (csp) gene, but subsequent studies have questioned this view. Here we present evidence that P. vivax in southern Mexico comprised 3 genetic populations whose distributions largely mirror those of the 2 mosquito vectors. Additionally, laboratory colony feeding experiments indicate that parasite populations are most compatible with sympatric mosquito species. Our results suggest that reciprocal selection between malaria parasites and mosquito vectors has led to local adaptation of the parasite. Adaptation to local vectors may play an important role in generating population structure in Plasmodium. A better understanding of coevolutionary dynamics between sympatric mosquitoes and parasites will facilitate the identification of molecular mechanisms relevant to disease transmission in nature and provide crucial information for malaria control.

Hitchhiking mapping reveals a candidate genomic region for natural selection in three-spined stickleback chromosome VIII

Identification of genes and genomic regions under directional natural selection has become one of the major goals in evolutionary genetics, but relatively little work to this end has been done by applying hitchhiking mapping to wild populations. Hitchhiking mapping starts from a genome scan using a randomly spaced set of molecular markers followed by a fine-scale analysis in the flanking regions of the candidate regions under selection. We used the hitchhiking mapping approach to narrow down a selective sweep in the genomic region flanking a candidate locus (Stn90) in chromosome VIII in the three-spined stickleback (Gasterosteus aculeatus). Twenty-four microsatellite markers were screened in an approximately 800-kb region around the candidate locus in three marine and four freshwater populations. The patterns of genetic diversity and differentiation in the candidate region were compared to those of a putatively neutral set of markers. The Bayesian FST-test indicated an elevated genetic differentiation, deviating significantly from neutral expectations, at a continuous region of approximately 20 kb upstream from the candidate locus. Furthermore, a method developed for an array of microsatellite markers rejected neutrality in a region of approximately 90 kb flanking the candidate locus supporting the selective sweep hypothesis. Likewise, the genomewide pattern of genetic diversity differed from the candidate region in a bottleneck analysis suggesting that selection, rather than demography, explains the reduced genetic diversity at the candidate interval. The neutrality tests suggest that the selective sweep had occurred mainly in the Lake Pulmanki population, but the results from bottleneck analyses indicate that selection might have operated in other populations as well. These results suggest that the narrow interval around locus Stn90 has likely been under directional selection, but the region contains several predicted genes, each of which can be the actual targets of selection. Understanding of the functional significance of this genomic region in an ecological context will require a more detailed sequence analysis.

Molecular epidemiology of malaria

Malaria persists as an undiminished global problem, but the resources available to address it have increased. Many tools for understanding its biology and epidemiology are well developed, with a particular richness of comparative genome sequences. Targeted genetic manipulation is now effectively combined with in vitro culture assays on the most important human parasite, Plasmodium falciparum, and with in vivo analysis of rodent and monkey malaria parasites in their laboratory hosts. Studies of the epidemiology, prevention, and treatment of human malaria have already been influenced by the availability of molecular methods, and analyses of parasite polymorphisms have long had useful and highly informative applications. However, the molecular epidemiology of malaria is currently undergoing its most substantial revolution as a result of the genomic information and technologies that are available in well-resourced centers. It is a challenge for research agendas to face the real needs presented by a disease that largely exists in extremely resource-poor settings, but it is one that there appears to be an increased willingness to undertake. To this end, developments in the molecular epidemiology of malaria are reviewed here, emphasizing aspects that may be current and future priorities.

Microsatellite polymorphism within pfcrt provides evidence of continuing evolution of chloroquine-resistant alleles in Papua New Guinea

Background

Polymorphism in the pfcrt gene underlies Plasmodium falciparum chloroquine resistance (CQR), as sensitive strains consistently carry lysine (K), while CQR strains carry threonine (T) at the codon 76. Previous studies have shown that microsatellite (MS) haplotype variation can be used to study the evolution of CQR polymorphism and to characterize intra- and inter-population dispersal of CQR in Papua New Guinea (PNG).

Methods

Here, following identification of new polymorphic MS in introns 2 and 3 within the pfcrt gene (msint2 and msint3, respectively), locus-by-locus and haplotype heterozygosity (H) analyses were performed to determine the distribution of this intronic polymorphism among pfcrt chloroquine-sensitive and CQR alleles.

Results

For MS flanking the pfcrt CQR allele, H ranged from 0.07 (B5M77, -18 kb) to 0.094 (9B12, +2 kb) suggesting that CQ selection pressure was responsible for strong homogenisation of this gene locus. In a survey of 206 pfcrt-SVMNT allele-containing field samples from malaria-endemic regions of PNG, H for msint2 was 0.201. This observation suggests that pfcrt msint2 exhibits a higher level of diversity than what is expected from the analyses of pfcrt flanking MS. Further analyses showed that one of the three haplotypes present in the early 1980's samples has become the predominant haplotype (frequency = 0.901) in CQR parasite populations collected after 1995 from three PNG sites, when CQR had spread throughout malaria-endemic regions of PNG. Apparent localized diversification of pfcrt haplotypes at each site was also observed among samples collected after 1995, where minor CQR-associated haplotypes were found to be unique to each site.

Conclusion

In this study, a higher level of diversity at MS loci within the pfcrt gene was observed when compared with the level of diversity at pfcrt flanking MS. While pfcrt (K76T) and its immediate flanking region indicate homogenisation in PNG CQR parasite populations, pfcrt intronic MS variation provides evidence that the locus is still evolving. Further studies are needed to determine whether these intronic MS introduce the underlying genetic mechanisms that may generate pfcrt allelic diversity.

Common origin and fixation of Plasmodium falciparum dhfr and dhps mutations associated with sulfadoxine-pyrimethamine resistance in a low-transmission area in South America

Recent studies indicated that sensitive parasites could increase in frequency in a population when drugs are removed, suggesting that the life span of affordable antimalarial drugs could be expanded. We studied 97 samples from Bolivar State, Venezuela, an area where sulfadoxine-pyrimethamine (SP) has not been used for 8 years due to its ineffectiveness. We characterized point mutations in two genes that have been implicated in resistance to SP, dihydrofolate reductase (dhfr) and dihydropteroate synthase (dhps). We also assayed neutral microsatellite markers around the dhfr (chromosome 4) and dhps (chromosome 8) genes and on chromosomes 2 and 3 to track the origin and spread of resistant alleles. We found that drug-resistant SP mutants are fixed in the population. Two genotypes were present in the samples, dhfr(50R/51I/108N) dhps(437G/540E/581G) (90.7%) and dhfr(51I/108N) dhps(437G/581G) (9.3%). We show a single microsatellite haplotype for all of the dhfr and dhps alleles, and the alleles at the microsatellite loci are different from those present in Africa. Thus, in these samples from Venezuela, there is a single origin for both dhfr and dhps SP-resistant alleles, and these alleles originated independently of those characterized from Africa. Furthermore, this is the first report of a "hitchhiking effect" on the genetic variation around dhps due to selection by SP using an extensive set of microsatellite markers. Our results indicate that, in areas where there is limited gene flow, the fixation of drug-resistant parasites in the population is stable, even after drug selection is relaxed.

Characterization and application of multiple genetic markers for Plasmodium malariae

Plasmodium malariae, a protozoan parasite that causes malaria in humans, has a global distribution in tropical and subtropical regions and is commonly found in sympatry with other Plasmodium species of humans. Little is known about the genetics or population structure of P. malariae. In the present study, we describe polymorphic genetic markers for P. malariae and present the first molecular epidemiological data for this parasite. Six microsatellite or minisatellite markers were validated using 76 P. malariae samples from a diverse geographical range. The repeat unit length varied from 2 to 17 bp, and up to 10 different alleles per locus were detected. Multiple genotypes of P. malariae were detected in 33 of 70 samples from humans with naturally acquired infection. Heterozygosity was calculated to be between 0.236 and 0.811. Allelic diversity was reduced for samples from South America and, at some loci, in samples from Thailand compared with those from Malawi. The number of unique multilocus genotypes defined using the 6 markers was significantly greater in Malawi than in Thailand, even when data from single genotype infections were used. There was a significant reduction in the multiplicity of infection in symptomatic infections compared with asymptomatic ones, suggesting that clinical episodes are usually caused by the expansion of a single genotype.

Defining the origin of Plasmodium falciparum resistant dhfr isolates in Senegal

We previously reported a high baseline prevalence of mutations in the dhfr and dhps genes of Plasmodium falciparum throughout Senegal. The highest prevalence of the triple dhfr pyrimethamine associated mutations were found in isolates obtained in the western part of the country near the capital city of Dakar. In this study, we sought out to determine the relatedness of dhfr wild type and mutated strains by analyzing three microsatellite regions upstream of the dhfr locus. Twenty-six of the 31 wild type strains had a unique microsatellite pattern. In contrast, of the 17 isolates containing the triple mutation in dhfr, 11 had an identical microsatellite pattern. Diverse geographical isolates in Senegal containing the triple dhfr mutation have arisen from a limited number of ancestral strains. In addition, we demonstrate that these isolates have shared ancestry with the previously reported triple mutation haplotype found in Tanzania, South Africa, and southeast Asia. This common ancestry may have implications for the malaria control strategy for reducing the spread of sulfadoxine-pyrimethamine resistance in Senegal and elsewhere in Africa.

Geographical distribution of selected and putatively neutral SNPs in Southeast Asian malaria parasites

Loci targeted by directional selection are expected to show elevated geographical population structure relative to neutral loci, and a flurry of recent papers have used this rationale to search for genome regions involved in adaptation. Studies of functional mutations that are known to be under selection are particularly useful for assessing the utility of this approach. Antimalarial drug treatment regimes vary considerably between countries in Southeast Asia selecting for local adaptation at parasite loci underlying resistance. We compared the population structure revealed by 10 nonsynonymous mutations (nonsynonymous single-nucleotide polymorphisms [nsSNPs]) in four loci that are known to be involved in antimalarial drug resistance, with patterns revealed by 10 synonymous mutations (synonymous single-nucleotide polymorphisms [sSNPs]) in housekeeping genes or genes of unknown function in 755 Plasmodium falciparum infections collected from 13 populations in six Southeast Asian countries. Allele frequencies at known nsSNPs underlying resistance varied markedly between locations (F(ST) = 0.18-0.66), with the highest frequencies on the Thailand-Burma border and the lowest frequencies in neighboring Lao PDR. In contrast, we found weak but significant geographic structure (F(ST) = 0-0.14) for 8 of 10 sSNPs. Importantly, all 10 nsSNPs showed significantly higher F(ST) (P < 8 x 10(-5)) than simulated neutral expectations based on observed F(ST) values in the putatively neutral sSNPs. This result was unaffected by the methods used to estimate allele frequencies or the number of populations used in the simulations. Given that dense single-nucleotide polymorphism (SNP) maps and rapid SNP assay methods are now available for P. falciparum, comparing genetic differentiation across the genome may provide a valuable aid to identifying parasite loci underlying local adaptation to drug treatment regimes or other selective forces. However, the high proportion of polymorphic sites that appear to be under balancing selection (or linked to selected sites) in the P. falciparum genome violates the central assumption that selected sites are rare, which complicates identification of outlier loci, and suggests that caution is needed when using this approach.