Selection strength and hitchhiking around two anti-malarial resistance genes

Molecular surveillance for anti-malarial drug resistance and genetic diversity of Plasmodium falciparum after chloroquine and sulfadoxine-pyrimethamine withdrawal in Quibdo, Colombia, 2018

Background

Resistance to anti-malarial drugs is associated with polymorphisms in target genes and surveillance for these molecular markers is important to detect the emergence of mutations associated with drug resistance and signal recovering sensitivity to anti-malarials previously used.

Methods

The presence of polymorphisms in genes associated with Plasmodium falciparum resistance to chloroquine and sulfadoxine-pyrimethamine was evaluated by Sanger sequencing, in 85 P. falciparum day of enrollment samples from a therapeutic efficacy study of artemether–lumefantrine conducted in 2018–2019 in Quibdo, Colombia. Samples were genotyped to assess mutations in pfcrt (codons 72–76), pfdhfr (codons 51, 59, 108, and 164), and pfdhps genes (codons 436, 437, 540, and 581). Further, the genetic diversity of infections using seven neutral microsatellites (NMSs) (C2M34, C3M69, Poly α, TA1, TA109, 2490, and PfPK2) was assessed.

Results

All isolates carried mutant alleles for pfcrt (K76T and N75E), and for pfdhfr (N51I and S108N), while for pfdhps, mutations were observed only for codon A437G (32/73, 43.8%). Fifty samples (58.8%) showed a complete neutral microsatellites (NMS) profile. The low mean number of alleles (2 ± 0.57) per locus and mean expected heterozygosity (0.17 ± 0.03) showed a reduced genetic diversity. NMS multilocus genotypes (MMG) were built and nine MMG were identified.

Conclusions

Overall, these findings confirm the fixation of chloroquine and pyrimethamine-resistant alleles already described in the literature, implying that these drugs are not currently appropriate for use in Colombia. In contrast, mutations in the pfdhps gene were only observed at codon 437, an indication that full resistance to sulfadoxine has not been achieved in Choco. MMGs found matched the clonal lineage E variant 1 previously reported in northwestern Colombia.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04328-x.

The many definitions of multiplicity of infection

The presence of multiple genetically different pathogenic variants within the same individual host is common in infectious diseases. Although this is neglected in some diseases, it is well recognized in others like malaria, where it is typically referred to as multiplicity of infection (MOI) or complexity of infection (COI). In malaria, with the advent of molecular surveillance, data is increasingly being available with enough resolution to capture MOI and integrate it into molecular surveillance strategies. The distribution of MOI on the population level scales with transmission intensities, while MOI on the individual level is a confounding factor when monitoring haplotypes of particular interests, e.g., those associated with drug-resistance. Particularly, in high-transmission areas, MOI leads to a discrepancy between the likelihood of a haplotype being observed in an infection (prevalence) and its abundance in the pathogen population (frequency). Despite its importance, MOI is not universally defined. Competing definitions vary from verbal ones to those based on concise statistical frameworks. Heuristic approaches to MOI are popular, although they do not mine the full potential of available data and are typically biased, potentially leading to misinferences. We introduce a formal statistical framework and suggest a concise definition of MOI and its distribution on the host-population level. We show how it relates to alternative definitions such as the number of distinct haplotypes within an infection or the maximum number of alleles detectable across a set of genetic markers. It is shown how alternatives can be derived from the general framework. Different statistical methods to estimate the distribution of MOI and pathogenic variants at the population level are discussed. The estimates can be used as plug-ins to reconstruct the most probable MOI of an infection and set of infecting haplotypes in individual infections. Furthermore, the relation between prevalence of pathogenic variants and their frequency (relative abundance) in the pathogen population in the context of MOI is clarified, with particular regard to seasonality in transmission intensities. The framework introduced here helps to guide the correct interpretation of results emerging from different definitions of MOI. Especially, it excels comparisons between studies based on different analytical methods.

Genomic Approaches to Drug Resistance in Malaria

Although the last two decades have seen a substantial decline in malaria incidence and mortality due to the use of insecticide-treated bed nets and artemisinin combination therapy, the threat of drug resistance is a constant obstacle to sustainable malaria control. Given that patients can die quickly from this disease, public health officials and doctors need to understand whether drug resistance exists in the parasite population, as well as how prevalent it is so they can make informed decisions about treatment. As testing for drug efficacy before providing treatment to malaria patients is impractical, researchers need molecular markers of resistance that can be more readily tracked in parasite populations. To this end, much work has been done to unravel the genetic underpinnings of drug resistance in Plasmodium falciparum. The aim of this review is to provide a broad overview of common genomic approaches that have been used to discover the alleles that drive drug response phenotypes in the most lethal human malaria parasite.

Study of the epidemiological behavior of malaria in the Darien Region, Panama. 2015-2017

BACKGROUND

Malaria is endemic in Darién and an assessment of the different factors affecting its epidemiology is crucial for the development of adequate strategies of surveillance, prevention, and disease control. The objective of this study was to determine the main characteristics of the epidemiological behavior of malaria in the Darien region.

METHODS

This research was comprised of a retrospective analysis to determine the incidence and malaria distribution in the Darien region from 2015 to 2017. We evaluated malaria indicators, disease distribution, incidence (by age group and sex), diagnostic methods, treatment, and control measures. In addition, we examined the cross-border migration activity and its possible contribution to the maintenance and distribution of malaria.

RESULTS

During the period of 2015-2017, we examined 41,141 thick blood smear samples, out of which 501 tested positive for malaria. Plasmodium vivax was responsible for 92.2% of those infections. Males comprised 62.7% of the total diagnosed cases. Meanwhile, a similar percentage, 62.7%, of the total cases were registered in economically active ages. The more frequent symptoms included fever (99.4%) and chills (97.4%), with 53.1% of cases registering between 2,000 and 6,000 parasites/μl of blood. The annual parasitic incidence (API) average was 3.0/1,000 inhabitants, while the slide positivity rate (SPR) was 1.2% and the annual blood examination rate (ABER) 22.5%. In Darién there is a constant internal and cross-border migration movement between Panama and Colombia. Malaria control measures consisted of the active and passive search of suspected cases and of the application of vector control measures.

CONCLUSION

This study provides an additional perspective on malaria epidemiology in Darién. Additional efforts are required to intensify malaria surveillance and to achieve an effective control, eventually moving closer to the objective of malaria elimination. At the same time, there is a need for more eco-epidemiological, entomological and migratory studies to determine how these factors contribute to the patterns of maintenance and dissemination of malaria.

Beyond drugs: the evolution of genes involved in human response to medications

The genetic variation of our species reflects human demographic history and adaptation to diverse local environments. Part of this genetic variation affects individual responses to exogenous substances, such as food, pollutants and drugs, and plays an important role in drug efficacy and safety. This review provides a synthesis of the evolution of loci implicated in human pharmacological response and metabolism, interpreted within the theoretical framework of population genetics and molecular evolution. In particular, I review and discuss key evolutionary aspects of different pharmacogenes in humans and other species, such as the relationship between the type of substrates and rate of evolution; the selective pressure exerted by landscape variables or dietary habits; expected and observed patterns of rare genetic variation. Finally, I discuss how this knowledge can be translated directly or after the implementation of specific studies, into practical guidelines.

Multiple Novel Mutations in Plasmodium falciparum Chloroquine Resistance Transporter Gene during Implementation of Artemisinin Combination Therapy in Thailand

Mutations in the chloroquine resistance transporter gene of Plasmodium falciparum (Pfcrt) are associated with drug susceptibility status of chloroquine and other antimalarials that interfere with heme detoxification process including artemisinin. We aim to investigate whether an increase in duration of artemisinin combination therapy (ACT) in Thailand could affect mutations in Pfcrt. The complete coding sequences of Pfcrt and dihydrofolate reductase (Pfdhfr), and size polymorphisms of the merozoite surface proteins-1 and 2 (Pfmsp-1 and Pfmsp-2) of 189 P. falciparum isolates collected during 1991 and 2016 were analyzed. In total, 12 novel amino acid substitutions and 13 novel PfCRT haplotypes were identified. The most prevalent haplotype belonged to the Dd2 sequence and no wild type was found. A significant positive correlation between the frequency of Pfcrt mutants and the year of sample collection was observed during nationwide ACT implementation (r = 0.780; P = 0.038). The number of haplotypes and nucleotide diversity of isolates collected during 3-day ACT (2009-2016) significantly outnumbered those collected before this treatment regimen. Positive Darwinian selection occurred in the transmembrane domains only among isolates collected during 3-day ACT but not among those collected before this period. No remarkable change was observed in the molecular indices for other loci analyzed when similar comparisons were performed. An increase in the duration of artesunate in combination therapy in Thailand could exert selective pressure on the Pfcrt locus, resulting in emergence of novel variants. The impact of these novel haplotypes on antimalarial susceptibilities requires further study.

Variation in intronic microsatellites and exon 2 of the Plasmodium falciparum chloroquine resistance transporter gene during modification of artemisinin combination therapy in Thailand

The amino acid substitution at residue 76 of the food vacuolar transmembrane protein encoded by the chloroquine resistance transporter gene of Plasmodium falciparum (Pfcrt) is an important, albeit imperfect, determinant of chloroquine susceptibility status of the parasite. Other mutations in Pfcrt can modulate susceptibility of P. falciparum to other antimalarials capable of interfering with heme detoxification process, and may exert compensatory effect on parasite growth rate. To address whether nationwide implementation of artemisinin combination therapy (ACT) in Thailand could affect sequence variation in exon 2 and introns of Pfcrt, we analyzed 136 P. falciparum isolates collected during 1997 and 2016 from endemic areas bordering Myanmar, Cambodia and Malaysia. Results revealed 6 haplotypes in exon 2 of Pfcrt with 2 novel substitutions at c.243A > G (p.R81) and c.251A > T (p.N84I). Positive selection was observed at amino acid residues 75, 76 and 97. Four, 3, and 2 alleles of microsatellite (AT/TA) repeats occurred in introns 1, 2 and 4, respectively, resulting in 7 different 3-locus haplotypes. The number of haplotypes and haplotype diversity of exon 2, and introns 1, 2 and 4 were significantly greater among isolates collected during 2009 and 2016 than those collected during 1997 and 2008 when 3-day ACT and 2-day ACT regimens were implemented nationwide, respectively (p < 0.05). By contrast, the number of haplotypes and haplotype diversity of the merozoite surface proteins 1 and 2 of these parasite populations did not differ significantly between these periods. Therefore, the Pfcrt locus of P. falciparum in Thailand continues to evolve and could have been affected by selective pressure from modification of ACT regimen.

Pooled-DNA sequencing identifies genomic regions of selection in Nigerian isolates of Plasmodium falciparum

Background

The burden of falciparum malaria is especially high in sub-Saharan Africa. Differences in pressure from host immunity and antimalarial drugs lead to adaptive changes responsible for high level of genetic variations within and between the parasite populations. Population-specific genetic studies to survey for genes under positive or balancing selection resulting from drug pressure or host immunity will allow for refinement of interventions.

Methods

We performed a pooled sequencing (pool-seq) of the genomes of 100 Plasmodium falciparum isolates from Nigeria. We explored allele-frequency based neutrality test (Tajima’s D) and integrated haplotype score (iHS) to identify genes under selection.

Results

Fourteen shared iHS regions that had at least 2 SNPs with a score > 2.5 were identified. These regions code for genes that were likely to have been under strong directional selection. Two of these genes were the chloroquine resistance transporter (CRT) on chromosome 7 and the multidrug resistance 1 (MDR1) on chromosome 5. There was a weak signature of selection in the dihydrofolate reductase (DHFR) gene on chromosome 4 and MDR5 genes on chromosome 13, with only 2 and 3 SNPs respectively identified within the iHS window. We observed strong selection pressure attributable to continued chloroquine and sulfadoxine-pyrimethamine use despite their official proscription for the treatment of uncomplicated malaria. There was also a major selective sweep on chromosome 6 which had 32 SNPs within the shared iHS region. Tajima’s D of circumsporozoite protein (CSP), erythrocyte-binding antigen (EBA-175), merozoite surface proteins - MSP3 and MSP7, merozoite surface protein duffy binding-like (MSPDBL2) and serine repeat antigen (SERA-5) were 1.38, 1.29, 0.73, 0.84 and 0.21, respectively.

Conclusion

We have demonstrated the use of pool-seq to understand genomic patterns of selection and variability in P. falciparum from Nigeria, which bears the highest burden of infections. This investigation identified known genomic signatures of selection from drug pressure and host immunity. This is evidence that P. falciparum populations explore common adaptive strategies that can be targeted for the development of new interventions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2260-z) contains supplementary material, which is available to authorized users.

Prospective Study of Plasmodium vivax Malaria Recurrence after Radical Treatment with a Chloroquine-Primaquine Standard Regimen in Turbo, Colombia

Plasmodium vivax recurrences help maintain malaria transmission. They are caused by recrudescence, reinfection, or relapse, which are not easily differentiated. A longitudinal observational study took place in Turbo municipality, Colombia. Participants with uncomplicated P. vivax infection received supervised treatment concomitantly with 25 mg/kg chloroquine and 0.25 mg/kg/day primaquine for 14 days. Incidence of recurrence was assessed over 180 days. Samples were genotyped, and origins of recurrences were established. A total of 134 participants were enrolled between February 2012 and July 2013, and 87 were followed for 180 days, during which 29 recurrences were detected. The cumulative incidence of first recurrence was 24.1% (21/87) (95% confidence interval [CI], 14.6 to 33.7%), and 86% (18/21) of these events occurred between days 51 and 110. High genetic diversity of P. vivax strains was found, and 12.5% (16/128) of the infections were polyclonal. Among detected recurrences, 93.1% (27/29) of strains were genotyped as genetically identical to the strain from the previous infection episode, and 65.5% (19/29) of infections were classified as relapses. Our results indicate that there is a high incidence of P. vivax malaria recurrence after treatment in Turbo municipality, Colombia, and that a large majority of these episodes are likely relapses from the previous infection. We attribute this to the primaquine regimen currently used in Colombia, which may be insufficient to eliminate hypnozoites.

Population Structure Shapes Copy Number Variation in Malaria Parasites

If copy number variants (CNVs) are predominantly deleterious, we would expect them to be more efficiently purged from populations with a large effective population size (Ne) than from populations with a small Ne. Malaria parasites (Plasmodium falciparum) provide an excellent organism to examine this prediction, because this protozoan shows a broad spectrum of population structures within a single species, with large, stable, outbred populations in Africa, small unstable inbred populations in South America and with intermediate population characteristics in South East Asia. We characterized 122 single-clone parasites, without prior laboratory culture, from malaria-infected patients in seven countries in Africa, South East Asia and South America using a high-density single-nucleotide polymorphism/CNV microarray. We scored 134 high-confidence CNVs across the parasite exome, including 33 deletions and 102 amplifications, which ranged in size from <500 bp to 59 kb, as well as 10,107 flanking, biallelic single-nucleotide polymorphisms. Overall, CNVs were rare, small, and skewed toward low frequency variants, consistent with the deleterious model. Relative to African and South East Asian populations, CNVs were significantly more common in South America, showed significantly less skew in allele frequencies, and were significantly larger. On this background of low frequency CNV, we also identified several high-frequency CNVs under putative positive selection using an FST outlier analysis. These included known adaptive CNVs containing rh2b and pfmdr1, and several other CNVs (e.g., DNA helicase and three conserved proteins) that require further investigation. Our data are consistent with a significant impact of genetic structure on CNV burden in an important human pathogen.

The effect of local variation in malaria transmission on the prevalence of sulfadoxine-pyrimethamine resistant haplotypes and selective sweep characteristics in Malawi

BACKGROUND

Persistence of sulfadoxine-pyrimethamine (SP) resistance has been described in an urban setting in Malawi where malaria transmission is relatively low. Higher malaria transmission is associated with greater genetic diversity and more frequent genetic recombination, which could lead to a more rapid re-emergence of SP-sensitive parasites, as well as more rapid degradation of selective sweeps. In this study, the impact of local variation in malaria transmission on the prevalence of SP-resistant haplotypes and selective sweep characteristics was investigated at an urban site with low parasite prevalence and two rural sites with moderate and high parasite prevalence.

METHODS

Samples from three sites with different parasite prevalence were genotyped for resistance markers within pfdhfr-ts and pfdhps and at microsatellites flanking these genes. Expected heterozygosity (He) was estimated to evaluate genetic diversity.

RESULTS

No difference in the prevalence of highly resistant DHFR 51I/59R/108N and DHPS 437G/540E was found between sites. Small differences in He flanking pfdhfr-ts and pfdhps were seen between rural-moderate and the other sites, as well as some shared haplotypes between the rural-high and urban-low sites.

CONCLUSIONS

The results do not show an effect of local variation in malaria transmission, as inferred from parasite prevalence, on SP-resistant haplotype prevalence.

Longitudinal analysis of Plasmodium falciparum genetic variation in Turbo, Colombia: implications for malaria control and elimination

Background

Malaria programmes estimate changes in prevalence to evaluate their efficacy. In this study, parasite genetic data was used to explore how the demography of the parasite population can inform about the processes driving variation in prevalence. In particular, how changes in treatment and population movement have affected malaria prevalence in an area with seasonal malaria.

Methods

Samples of Plasmodium falciparum collected over 8 years from a population in Turbo, Colombia were genotyped at nine microsatellite loci and three drug-resistance loci. These data were analysed using several population genetic methods to detect changes in parasite genetic diversity and population structure. In addition, a coalescent-based method was used to estimate substitution rates at the microsatellite loci.

Results

The estimated mean microsatellite substitution rates varied between 5.35 × 10⁻³ and 3.77 × 10⁻² substitutions/locus/month. Cluster analysis identified six distinct parasite clusters, five of which persisted for the full duration of the study. However, the frequencies of the clusters varied significantly between years, consistent with a small effective population size.

Conclusions

Malaria control programmes can detect re-introductions and changes in transmission using rapidly evolving microsatellite loci. In this population, the steadily decreasing diversity and the relatively constant effective population size suggest that an increase in malaria prevalence from 2004 to 2007 was primarily driven by local rather than imported cases.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-015-0887-9) contains supplementary material, which is available to authorized users.

How genomics is contributing to the fight against artemisinin-resistant malaria parasites

Plasmodium falciparum, the malignant malaria parasite, has developed resistance to artemisinin, the most important and widely used antimalarial drug at present. Currently confined to Southeast Asia, the spread of resistant parasites to Africa would constitute a public health catastrophe. In this review we highlight the recent contributions of genomics to our understanding how the parasite develops resistance to artemisinin and its derivatives, and how resistant parasites may be monitored and tracked in real-time, using molecular approaches.

Balancing drug resistance and growth rates via compensatory mutations in the Plasmodium falciparum chloroquine resistance transporter

The widespread use of chloroquine to treat Plasmodium falciparum infections has resulted in the selection and dissemination of variant haplotypes of the primary resistance determinant PfCRT. These haplotypes have encountered drug pressure and within-host competition with wild-type drug-sensitive parasites. To examine these selective forces in vitro, we genetically engineered P. falciparum to express geographically diverse PfCRT haplotypes. Variant alleles from the Philippines (PH1 and PH2, which differ solely by the C72S mutation) both conferred a moderate gain of chloroquine resistance and a reduction in growth rates in vitro. Of the two, PH2 showed higher IC50 values, contrasting with reduced growth. Furthermore, a highly mutated pfcrt allele from Cambodia (Cam734) conferred moderate chloroquine resistance and enhanced growth rates, when tested against wild-type pfcrt in co-culture competition assays. These three alleles mediated cross-resistance to amodiaquine, an antimalarial drug widely used in Africa. Each allele, along with the globally prevalent Dd2 and 7G8 alleles, rendered parasites more susceptible to lumefantrine, the partner drug used in the leading first-line artemisinin-based combination therapy. These data reveal ongoing region-specific evolution of PfCRT that impacts drug susceptibility and relative fitness in settings of mixed infections, and raise important considerations about optimal agents to treat chloroquine-resistant malaria.

Increasing prevalence of a novel triple-mutant dihydropteroate synthase genotype in Plasmodium falciparum in western Kenya

The molecular basis of sulfadoxine-pyrimethamine (SP) resistance lies in a combination of single-nucleotide polymorphisms (SNPs) in two genes coding for Plasmodium falciparum dihydrofolate reductase (Pfdhfr) and P. falciparum dihydropteroate synthase (Pfdhps), targeted by pyrimethamine and sulfadoxine, respectively. The continued use of SP for intermittent preventive treatment in pregnant women in many African countries, despite SP's discontinuation as a first-line antimalarial treatment option due to high levels of drug resistance, may further increase the prevalence of SP-resistant parasites and/or lead to the selection of new mutations. An antimalarial drug resistance surveillance study was conducted in western Kenya between 2010 and 2013. A total of 203 clinical samples from children with uncomplicated malaria were genotyped for SNPs associated with SP resistance. The prevalence of the triple-mutant Pfdhfr C50 I51R59N108: I164 genotype and the double-mutant Pfdhps S436 G437E540: A581A613 genotype was high. Two triple-mutant Pfdhps genotypes, S436 G437E540G581: A613 and H436G437E540: A581A613, were found, with the latter thus far being uniquely found in western Kenya. The prevalence of the S436 G437E540G581: A613 genotype was low. However, a steady increase in the prevalence of the Pfdhps triple-mutant H436G437E540: A581A613 genotype has been observed since its appearance in early 2000. Isolates with these genotypes shared substantial microsatellite haplotypes with the most common double-mutant allele, suggesting that this triple-mutant allele may have evolved locally. Overall, these findings show that the prevalence of the H436G437E540: A581A613 triple mutant may be increasing in this population and could compromise the efficacy of SP for intermittent preventive treatment in pregnant women if it increases the resistance threshold further.

Persistence of Sulfadoxine-Pyrimethamine Resistance Despite Reduction of Drug Pressure in Malawi

BACKGROUND

In 2007, Malawi replaced sulfadoxine-pyrimethamine (SP) with an artemisinin-based combination therapy as the first-line treatment for uncomplicated Plasmodium falciparum malaria in response to failing SP efficacy. Here we estimate the effect of reduced SP pressure on the prevalence of SP-resistant parasites and the characteristics of the associated selective sweeps flanking the resistance loci.

METHODS

Samples obtained from individuals with clinical malaria during a period of high SP use (1999-2001), a transitional period (2007-2008), and a period of low SP use (2012) were genotyped for resistance markers at pfdhfr-ts codons 51, 59, and 108 and pfdhps codons 437, 540, and 581. Expected heterozygosity was estimated to evaluate the genetic diversity flanking pfdhfr-ts and pfdhps.

RESULTS

An increase in the prevalence of the resistance haplotypes DHFR 51I/59R/108N and DHPS 437G/540E occurred under sustained drug pressure, with no change in haplotype prevalence 5 years after reduction in SP pressure. The DHPS 437G/540E/581G haplotype was observed in 2007 and increased in prevalence during a period of reduced SP pressure. Changes to the sweep characteristics flanking pfdhfr-ts and pfdhps were minimal.

CONCLUSIONS

In contrast to the rapid and complete return of chloroquine-susceptible falciparum malaria after chloroquine was withdrawn from Malawi, a reemergence of SP efficacy is unlikely in the near future.

Genetic architecture of resistance in Daphnia hosts against two species of host-specific parasites

Understanding the genetic architecture of host resistance is key for understanding the evolution of host-parasite interactions. Evolutionary models often assume simple genetics based on few loci and strong epistasis. It is unknown, however, whether these assumptions apply to natural populations. Using a quantitative trait loci (QTL) approach, we explore the genetic architecture of resistance in the crustacean Daphnia magna to two of its natural parasites: the horizontally transmitted bacterium Pasteuria ramosa and the horizontally and vertically transmitted microsporidium Hamiltosporidium tvaerminnensis. These two systems have become models for studies on the evolution of host-parasite interactions. In the QTL panel used here, Daphnia's resistance to P. ramosa is controlled by a single major QTL (which explains 50% of the observed variation). Resistance to H. tvaerminnensis horizontal infections shows a signature of a quantitative trait based in multiple loci with weak epistatic interactions (together explaining 38% variation). Resistance to H. tvaerminnensis vertical infections, however, shows only one QTL (explaining 13.5% variance) that colocalizes with one of the QTLs for horizontal infections. QTLs for resistance to Pasteuria and Hamiltosporidium do not colocalize. We conclude that the genetics of resistance in D. magna are drastically different for these two parasites. Furthermore, we infer that based on these and earlier results, the mechanisms of coevolution differ strongly for the two host-parasite systems. Only the Pasteuria-Daphnia system is expected to follow the negative frequency-dependent selection (Red Queen) model. How coevolution works in the Hamiltosporidium-Daphnia system remains unclear.

Population genomics of the killer whale indicates ecotype evolution in sympatry involving both selection and drift

The evolution of diversity in the marine ecosystem is poorly understood, given the relatively high potential for connectivity, especially for highly mobile species such as whales and dolphins. The killer whale (Orcinus orca) has a worldwide distribution, and individual social groups travel over a wide geographic range. Even so, regional populations have been shown to be genetically differentiated, including among different foraging specialists (ecotypes) in sympatry. Given the strong matrifocal social structure of this species together with strong resource specializations, understanding the process of differentiation will require an understanding of the relative importance of both genetic drift and local adaptation. Here we provide a high-resolution analysis based on nuclear single-nucleotide polymorphic markers and inference about differentiation at both neutral loci and those potentially under selection. We find that all population comparisons, within or among foraging ecotypes, show significant differentiation, including populations in parapatry and sympatry. Loci putatively under selection show a different pattern of structure compared to neutral loci and are associated with gene ontology terms reflecting physiologically relevant functions (e.g. related to digestion). The pattern of differentiation for one ecotype in the North Pacific suggests local adaptation and shows some fixed differences among sympatric ecotypes. We suggest that differential habitat use and resource specializations have promoted sufficient isolation to allow differential evolution at neutral and functional loci, but that the process is recent and dependent on both selection and drift.

High prevalence of pfcrt-CVIET haplotype in isolates from asymptomatic and symptomatic patients in south-central Oromia, Ethiopia

Background

As a result of extensive chloroquine resistance (CQR) in Plasmodium falciparum in late 1990s, Ethiopia replaced CQ with sulphadoxine-pyrimethamine (SP) as first-line drug, which in turn was replaced by artemisinin combination therapy in 2004. Plasmodium falciparum resistance to CQ is determined by the mutation at K76T of the P. falciparum chloroquine resistance transporter (pfcrt) gene. Understanding diversity in the P. falciparum genome is crucial since it has the potential to influence important phenotypes of the parasite such as drug resistance. Limited data is available regarding the type of pfcrt mutant allelic type, the effect of CQ withdrawal and diversity of the parasite population in south-central Oromia, Ethiopia.

Methods

Finger-pricked blood spotted on Whatman 3MM filter papers were collected from falciparum malaria patients. Parasite DNA was extracted from individual blood spots on the filter papers. The presence of K76T mutations was determined using nested PCR for all isolates. Complete sequencing of mutations in pfcrt 72-76 was done for a set of randomly selected resistant isolates. Four microsatellite (MS) markers were analysed to determine the heterozygosity.

Results

Although CQ was withdrawn for more than a decade, 100% of the parasites still carried the pfcrt K76T mutation. All isolates were mutant at the K76T polymorphism. Based on combinations of MS markers, seven different Ethiopian CQR variants (E1-E7) were identified. Heterozygosity (H_e) for MS flanking the pfcrt chloroquine resistance allele ranged from 0.00 (mscrt -29, -29.268 kb) to 0.21 (mscrt -2, -2.814 kb). H_e ranged from 0.00 (msint 3, 0 kb) to 0.19 (msint 2, 0 kb) for MS within the pfcrt gene. Both intronic and MS flanking the pfcrt gene showed low levels of diversity.

Conclusion

pfcrt CQR allele seems to be fixed in the study area. Of the different haplotypes associated with CQR, only the CVIET genotype was identified. No reversal to the wild-type has occurred in Ethiopia unlike in many Africa countries where CQR parasites declined after cessation of CQ use. Decreased diversity in CQR isolates surrounding pfcrt suggests CQ selection and homogenization among CQR parasite population. While mutation in msint 3 and mscrt -29 of the mutant pfcrt allele is being fixed, it seems that mutations in msint 2 and mscrt -2 are still evolving and may indicate the start of re-diversification of the population from a fixed 76 T population.

Genome-wide analysis of selection on the malaria parasite Plasmodium falciparum in West African populations of differing infection endemicity

Locally varying selection on pathogens may be due to differences in drug pressure, host immunity, transmission opportunities between hosts, or the intensity of between-genotype competition within hosts. Highly recombining populations of the human malaria parasite Plasmodium falciparum throughout West Africa are closely related, as gene flow is relatively unrestricted in this endemic region, but markedly varying ecology and transmission intensity should cause distinct local selective pressures. Genome-wide analysis of sequence variation was undertaken on a sample of 100 P. falciparum clinical isolates from a highly endemic region of the Republic of Guinea where transmission occurs for most of each year and compared with data from 52 clinical isolates from a previously sampled population from The Gambia, where there is relatively limited seasonal malaria transmission. Paired-end short-read sequences were mapped against the 3D7 P. falciparum reference genome sequence, and data on 136,144 single nucleotide polymorphisms (SNPs) were obtained. Within-population analyses identifying loci showing evidence of recent positive directional selection and balancing selection confirm that antimalarial drugs and host immunity have been major selective agents. Many of the signatures of recent directional selection reflected by standardized integrated haplotype scores were population specific, including differences at drug resistance loci due to historically different antimalarial use between the countries. In contrast, both populations showed a similar set of loci likely to be under balancing selection as indicated by very high Tajima’s D values, including a significant overrepresentation of genes expressed at the merozoite stage that invades erythrocytes and several previously validated targets of acquired immunity. Between-population F_ST analysis identified exceptional differentiation of allele frequencies at a small number of loci, most markedly for five SNPs covering a 15-kb region within and flanking the gdv1 gene that regulates the early stages of gametocyte development, which is likely related to the extreme differences in mosquito vector abundance and seasonality that determine the transmission opportunities for the sexual stage of the parasite.

Genetic variation and recurrent parasitaemia in Peruvian Plasmodium vivax populations

Background

Plasmodium vivax is a predominant species of malaria in parts of South America and there is increasing resistance to drugs to treat infections by P. vivax. The existence of latent hypnozoites further complicates the ability to classify recurrent infections as treatment failures due to relapse, recrudescence of hyponozoites or re-infections. Antigen loci are putatively under natural selection and may not be an optimal molecular marker to define parasite haplotypes in paired samples. Putatively neutral microsatellite loci, however, offer an assessment of neutral haplotypes. The objective here was to assess the utility of neutral microsatellite loci to reconcile cases of recurrent parasitaemia in Amazonian P. vivax populations in Peru.

Methods

Patient blood samples were collected from three locations in or around Iquitos in the Peruvian Amazon. Five putatively neutral microsatellite loci were characterized from 445 samples to ascertain the within and amongst population variation. A total of 30 day 0 and day of recurrent parasitaemia samples were characterized at microsatellite loci and five polymorphic antigen loci for haplotype classification.

Results

The genetic diversity at microsatellite loci was consistent with neutral levels of variation measured in other South American P. vivax populations. Results between antigen and microsatellite loci for the 30 day 0 and day of recurrent parasitaemia samples were the same for 80% of the pairs. The majority of non-concordant results were the result of differing alleles at microsatellite loci. This analysis estimates that 90% of the paired samples with the same microsatellite haplotype are unlikely to be due to a new infection.

Conclusions

A population-level approach was used to yield a better estimate of the probability of a new infection versus relapse or recrudescence of homologous hypnozoites; hypnozoite activation was common for this cohort. Population studies are critical with the evaluation of genetic markers to assess P. vivax biology and epidemiology. The additional demonstration of microsatellite loci as neutral markers capable of distinguishing the origin of the recurrent parasites (new infection or originating from the patient) lends support to their use in assessment of treatment outcomes.

Changes in malaria parasite drug resistance in an endemic population over a 25-year period with resulting genomic evidence of selection

Background. Analysis of genome-wide polymorphism in many organisms has potential to identify genes under recent selection. However, data on historical allele frequency changes are rarely available for direct confirmation.

Methods. We genotyped single nucleotide polymorphisms (SNPs) in 4 Plasmodium falciparum drug resistance genes in 668 archived parasite-positive blood samples of a Gambian population between 1984 and 2008. This covered a period before antimalarial resistance was detected locally, through subsequent failure of multiple drugs until introduction of artemisinin combination therapy. We separately performed genome-wide sequence analysis of 52 clinical isolates from 2008 to prospect for loci under recent directional selection.

Results. Resistance alleles increased from very low frequencies, peaking in 2000 for chloroquine resistance-associated crt and mdr1 genes and at the end of the survey period for dhfr and dhps genes respectively associated with pyrimethamine and sulfadoxine resistance. Temporal changes fit a model incorporating likely selection coefficients over the period. Three of the drug resistance loci were in the top 4 regions under strong selection implicated by the genome-wide analysis.

Conclusions. Genome-wide polymorphism analysis of an endemic population sample robustly identifies loci with detailed documentation of recent selection, demonstrating power to prospectively detect emerging drug resistance genes.

RTK/ERK pathway under natural selection associated with prostate cancer

Prostate cancer (PCa) is a global disease causing large numbers of deaths every year. Recent studies have indicated the RTK/ERK pathway might be a key pathway in the development of PCa. However, the exact association and evolution-based mechanism remain unclear. This study was conducted by combining genotypic and phenotypic data from the Chinese Consortium for Prostate Cancer Genetics (ChinaPCa) with related databases such as the HapMap Project and Genevar. In this analysis, expression of quantitative trait loci (eQTLs) analysis, natural selection and gene-based pathway analysis were involved. The pathway analysis confirmed the positive relationship between PCa risk and several key genes. In addition, combined with the natural selection, it seems that 4 genes (EGFR, ERBB2, PTK2, and RAF1) with five SNPs (rs11238349, rs17172438, rs984654, rs11773818, and rs17172432) especially rs17172432, might be pivotal factors in the development of PCa. The results indicate that the RTK/ERK pathway under natural selection is a key link in PCa risk. The joint effect of the genes and loci with positive selection might be one reason for the development of PCa. Dealing with all the factors simultaneously might give insight into prevention and aid in predicting the success of potential therapies for PCa.

Analyses of genetic variations at microsatellite loci present in-and-around the Pfcrt gene in Indian Plasmodium falciparum

Evolution and spread of chloroquine resistant (CQR) malaria parasite Plasmodium falciparum have posed great threat in malaria intervention across the globe. The occurrence of K76T mutation in the P. falciparum chloroquine resistance transporter (pfcrt) gene has been widely attributed to CQR with four neighboring mutations providing compensatory fitness benefit to the parasite survival. Understanding evolutionary patterns of the pfcrt gene is of great relevance not only for devising new malaria control measures but also could serve as a model to understand evolution and spread of other human drug-resistant pathogens. Several studies, mainly based on differential patterns of diversities of the microsatellite loci placed in-and-around the pfcrt gene have indicated the role of positive natural selection under the 'hitchhiking' model of molecular evolution. However, the studies were restricted to limited number of microsatellite loci present inside the pfcrt gene. Moreover, comparatively higher level of diversities in microsatellite loci present inside the pfcrt gene than the loci flanking the pfcrt gene are hallmarks of Indian P. falciparum, presenting contrasting evolutionary models to global isolates. With a view to infer evolutionary patterns of the pfcrt gene in Indian P. falciparum, we have adopted a unique sampling scheme of two types of populations (cultured and field collected) and utilized 20 polymorphic microsatellite loci (16 located inside the pfcrt gene and four in the two flanking regions) to disentangle between genetic drift (inbred cultured isolates) and natural selection (field isolates). Data analyses employing different population genetic tests could not straightforwardly explain either the model invoking 'genetic hitchhiking' or 'genetic drift'. However, complex evolutionary models influenced by both demography and natural selection or an alternative model of natural selection (e.g. diversifying/balancing selection) might better explain the observed microsatellite variation in-and-around the pfcrt gene in Indian P. falciparum.

A comprehensive survey of polymorphisms conferring anti-malarial resistance in Plasmodium falciparum across Pakistan

Background

Few studies have been conducted in Pakistan to determine the efficacy of chloroquine and sulphadoxine-pyrimethamine (SP), which remain in use as treatment for Plasmodium vivax and in combination with artesunate to treat Plasmodium falciparum, respectively. In this study, samples from several sites across Pakistan were characterized to determine prevalence of molecular resistance markers in the P. falciparum chloroquine resistance transporter (pfcrt), multidrug resistance (pfmdr1), dihydrofolate reductase (pfdhfr) and dihydropteroate synthase (pfdhps) genes and the origin of chloroquine-resistant P. falciparum parasites.

Methods

Microscopy-confirmed malaria parasite-positive blood samples from 801 patients across the country were collected in 2011. Of these, 171 infections were identified by polymerase chain reaction (PCR) as P. falciparum and analysed by pyrosequencing for mutations conferring chloroquine resistance (pfcrt codons 72–76), multidrug resistance (pfmdr1 N86Y, Y184F, S1034C, N1042D and D1246Y), pyrimethamine resistance (pfdhfr, C50R, N51I, C59R, S108N and I164L) and sulphadoxine resistance (pfdhps, S436A, A437G, K540E, A581G and A613T/S). pfmdr1 gene copy number variation was determined by real-time PCR, and microsatellites flanking the pfcrt locus were typed to determine the origin of the chloroquine-resistant haplotype.

Results

The pfcrt K76T mutation was found in all samples as part of the S72/V73/M74/N75/T76 (SVMNT) haplotype. Microsatellites flanking pfcrt showed high similarity to the signature found in India and Papua New Guinea. pfmdr1 N86Y was found in 20% of samples and all samples harboured a single copy of the pfmdr1 gene. The pfdhfr double mutation C59R + S108N was present in 87% of samples while the pfdhfr triple mutant (N51I + C59R + S108N) was not detected. Pfdhps A437G was found in 60% of samples. Pure pfdhps K540E was rare, at 4%, but mixed genotype 540 K/E was found in 77% of samples. Similarly, pure pfdhps A581G was found in 4% of the isolates while mixed 581A/G was found in 39% of samples.

Conclusions

These results suggest an emerging problem with multidrug resistant P. falciparum in Pakistan. The chloroquine resistance genotype has reached complete fixation in the population, with a microsatellite pattern indicative of a selective sweep. Moreover, the prevalence of mutations in both pfdhfr and pfdhps, albeit without the presence of the pfdhfr triple mutant, indicates that continued monitoring is warranted to assess whether SP remains efficacious as a partner drug for artesunate for the treatment of P. falciparum.

Genetic and genomic approaches for the discovery of parasite genes involved in antimalarial drug resistance

The biggest threat to the war on malaria is the continued evolution of drug resistance by the parasite. Resistance to almost all currently available antimalarials now exists in Plasmodium falciparum which causes the most suffering among all human malaria parasites. Monitoring of antimalarial efficacy and the development and subsequent spread of resistance has become an important part in the treatment and control of malaria. With recent reports of reduced efficacy of artemisinin, the current recommended treatment for uncomplicated malaria, there is urgent need for better methods to recognize and monitor drug resistance for effective treatment. Molecular markers have become a welcome addition to complement the more laborious and costly in vitro and in vivo methods that have traditionally been used to monitor drug resistance. However, there are currently no molecular markers for resistance to some antimalarials. This review highlights the role of the various genetic and genomic approaches that have been used in identifying the molecular markers that underlie drug resistance in P. falciparum. These approaches include; candidate genes, genetic linkage and genome-wide association studies. We discuss the requirements and limitations of each approach and use various examples to illustrate their contributions in identifying genomic regions of the parasite associated with antimalarial drug responses.

Microsatellite analysis of chloroquine resistance associated alleles and neutral loci reveal genetic structure of Indian Plasmodium falciparum

Efforts to control malignant malaria caused by Plasmodium falciparum are hampered by the parasite's acquisition of resistance to antimalarial drugs, e.g., chloroquine. This necessitates evaluating the spread of chloroquine resistance in any malaria-endemic area. India displays highly variable malaria epidemiology and also shares porous international borders with malaria-endemic Southeast Asian countries having multi-drug resistant malaria. Malaria epidemiology in India is believed to be affected by two major factors: high genetic diversity and evolving drug resistance in P. falciparum. How transmission intensity of malaria can influence the genetic structure of chloroquine-resistant P. falciparum population in India is unknown. Here, genetic diversity within and among P. falciparum populations is analyzed with respect to their prevalence and chloroquine resistance observed in 13 different locations in India. Microsatellites developed for P. falciparum, including three putatively neutral and seven microsatellites thought to be under a hitchhiking effect due to chloroquine selection were used. Genetic hitchhiking is observed in five of seven microsatellites flanking the gene responsible for chloroquine resistance. Genetic admixture analysis and F-statistics detected genetically distinct groups in accordance with transmission intensity of different locations and the probable use of chloroquine. A large genetic break between the chloroquine-resistant parasite of the Northeast-East-Island group and Southwest group (FST=0.253, P<0.001) suggests a long period of isolation or a possibility of different origin between them. A pattern of significant isolation by distance was observed in low transmission areas (r=0.49, P=0.003, N=83, Mantel test). An unanticipated pattern of spread of hitchhiking suggests genetic structure for Indian P. falciparum population. Overall, the study suggests that transmission intensity can be an efficient driver for genetic differentiation at both neutral and adaptive loci across India.

In vivo efficacy of artemether-lumefantrine and chloroquine against Plasmodium vivax: a randomized open label trial in central Ethiopia

BACKGROUND

In vivo efficacy assessments of antimalarials are essential for ensuring effective case management. In Ethiopia, chloroquine (CQ) without primaquine is the first-line treatment for Plasmodium vivax in malarious areas, but artemether-lumefantrine (AL) is also commonly used.

METHODS AND FINDINGS

In 2009, we conducted a 42-day efficacy study of AL or CQ for P. vivax in Oromia Regional State, Ethiopia. Individuals with P. vivax monoinfection were enrolled. Primary endpoint was day 28 cure rate. In patients with recurrent parasitemia, drug level and genotyping using microsatellite markers were assessed. Using survival analysis, uncorrected patient cure rates at day 28 were 75.7% (95% confidence interval (CI) 66.8-82.5) for AL and 90.8% (95% CI 83.6-94.9) for CQ. During the 42 days of follow-up, 41.6% (47/113) of patients in the AL arm and 31.8% (34/107) in the CQ arm presented with recurrent P. vivax infection, with the median number of days to recurrence of 28 compared to 35 days in the AL and CQ arm, respectively. Using microsatellite markers to reclassify recurrent parasitemias with a different genotype as non-treatment failures, day 28 cure rates were genotype adjusted to 91.1% (95% CI 84.1-95.1) for AL and to 97.2% (91.6-99.1) for CQ. Three patients (2.8%) with recurrent parasitemia by day 28 in the CQ arm were noted to have drug levels above 100 ng/ml.

CONCLUSIONS

In the short term, both AL and CQ were effective and well-tolerated for P. vivax malaria, but high rates of recurrent parasitemia were noted with both drugs. CQ provided longer post-treatment prophylaxis than AL, resulting in delayed recurrence of parasitemia. Although the current policy of species-specific treatment can be maintained for Ethiopia, the co-administration of primaquine for treatment of P. vivax malaria needs to be urgently considered to prevent relapse infections.

TRIAL REGISTRATION

ClinicalTrials.gov NCT01052584.

Genetic hitchhiking under heterogeneous spatial selection pressures

During adaptive evolutionary processes substantial heterogeneity in selective pressure might act across local habitats in sympatry. Examples are selection for drug resistance in malaria or herbicide resistance in weeds. In such setups standard population-genetic assumptions (homogeneous constant selection pressures, random mating etc.) are likely to be violated. To avoid misinferences on the strength and pattern of natural selection it is therefore necessary to adjust population-genetic theory to meet the specifics driving adaptive processes in particular organisms. We introduce a deterministic model in which selection acts heterogeneously on a population of haploid individuals across different patches over which the population randomly disperses every generation. A fixed proportion of individuals mates exclusively within patches, whereas the rest mates randomly across all patches. We study how the allele frequencies at neutral markers are affected by the spread of a beneficial mutation at a closely linked locus (genetic hitchhiking). We provide an analytical solution for the frequency change and the expected heterozygosity at the neutral locus after a single copy of a beneficial mutation became fixed. We furthermore provide approximations of these solutions which allow for more obvious interpretations. In addition, we validate the results by stochastic simulations. Our results show that the application of standard population-genetic theory is accurate as long as differences across selective environments are moderate. However, if selective differences are substantial, as for drug resistance in malaria, herbicide resistance in weeds, or insecticide resistance in agriculture, it is necessary to adapt available theory to the specifics of particular organisms.

Historical shifts in Brazilian P. falciparum population structure and drug resistance alleles

Previous work suggests that Brazilian Plasmodium falciparum has limited genetic diversity and a history of bottlenecks, multiple reintroductions due to human migration, and clonal expansions. We hypothesized that Brazilian P. falciparum would exhibit clonal structure. We examined isolates collected across two decades from Amapá, Rondônia, and Pará state (n = 190). By examining more microsatellites markers on more chromosomes than previous studies, we hoped to define the extent of low diversity, linkage disequilibrium, bottlenecks, population structure, and parasite migration within Brazil. We used retrospective genotyping of samples from the 1980s and 1990s to explore the population genetics of SP resistant dhfr and dhps alleles. We tested an existing hypothesis that the triple mutant dhfr mutations 50R/51I/108N and 51I/108N/164L developed in southern Amazon from a single origin of common or similar parasites. We found that Brazilian P. falciparum had limited genetic diversity and isolation by distance was rejected, which suggests it underwent bottlenecks followed by migration between sites. Unlike Peru, there appeared to be gene flow across the Brazilian Amazon basin. We were unable to divide parasite populations by clonal lineages and pairwise FST were common. Most parasite diversity was found within sites in the Brazilian Amazon, according to AMOVA. Our results challenge the hypothesis that triple mutant alleles arose from a single lineage in the Southern Amazon. SP resistance, at both the double and triple mutant stages, developed twice and potentially in different regions of the Brazilian Amazon. We would have required samples from before the 1980s to describe how SP resistance spread across the basin or describe the complex internal migration of Brazilian parasites after the colonization efforts of past decades. The Brazilian Amazon basin may have sufficient internal migration for drug resistance reported in any particular region to rapidly spread to other parts of basin under similar drug pressure.

Sulfadoxine-pyrimethamine resistance in Plasmodium falciparum: a zoomed image at the molecular level within a geographic context

Antimalarial chemotherapy is one of the main pillars in the prevention and control of malaria. Following widespread resistance of Plasmodium falciparum to chloroquine, sulfadoxine-pyrimethamine came to the scene as an alternative to the cheap and well-tolerated chloroquine. However, widespread resistance to sulfadoxine-pyrimethamine has been documented. In vivo efficacy tests are the gold standard for assessing drug resistance and treatment failure. However, they have many disadvantages, such as influence of host immunity and drug pharmacokinetics. In vitro tests of antimalarial drug efficacy also have many technical difficulties. Molecular markers of resistance have emerged as epidemiologic tools to investigate antimalarial drug resistance even before becoming clinically evident. Mutations in P. falciparum dihydrofolate reductase and dihydrofolate synthase have been extensively studied as molecular markers for resistance to pyrimethamine and sulfadoxine, respectively. This review highlights the resistance of P. falciparum at the molecular level presenting both supporting and opposing studies on the utility of molecular markers.

Local population structure of Plasmodium: impact on malaria control and elimination

BACKGROUND

Regardless of the growing interest in detecting population structures in malarial parasites, there have been limited discussions on how to use this concept in control programmes. In such context, the effects of the parasite population structures will depend on interventions' spatial or temporal scales. This investigation explores the problem of identifying genetic markers, in this case microsatellites, to unveil Plasmodium genetic structures that could affect decisions in the context of elimination. The study was performed in a low-transmission area, which offers a good proxy to better understand problems associated with surveillance at the final stages of malaria elimination.

METHODS

Plasmodium vivax samples collected in Tumeremo, Venezuela, between March 2003 and November 2004 were analysed. Since Plasmodium falciparum also circulates in many low endemic areas, P. falciparum samples from the same locality and time period were included for comparison. Plasmodium vivax samples were assayed for an original set of 25 microsatellites and P. falciparum samples were assayed for 12 microsatellites.

RESULTS

Not all microsatellite loci assayed offered reliable local data. A complex temporal-cluster dynamics is found in both P. vivax and P. falciparum. Such dynamics affect the numbers and the type of microsatellites required for identifying individual parasites or parasite clusters when performing cross-sectional studies. The minimum number of microsatellites required to differentiate circulating P. vivax clusters differs from the minimum number of hyper-variable microsatellites required to distinguish individuals within these clusters. Regardless the extended number of microsatellites used in P. vivax, it was not possible to separate all individual infections.

CONCLUSIONS

Molecular surveillance has great potential; however, it requires preliminary local studies in order to properly interpret the emerging patterns in the context of elimination. Clonal expansions and clusters turnovers need to be taken into account when using molecular markers. Those affect the number and type of microsatellite markers, as well as, the expected genetic patterns in the context of operational investigations. By considering the local dynamics, elimination programmes could cost-effectively use molecular markers. However, population level studies need to consider the local limitations of a given set of loci in terms of providing epidemiologically relevant information.

Large-scale survey for novel genotypes of Plasmodium falciparum chloroquine-resistance gene pfcrt

Background

In Plasmodium falciparum, resistance to chloroquine (CQ) is conferred by a K to T mutation at amino acid position 76 (K76T) in the P. falciparum CQ transporter (PfCRT). To date, at least 15 pfcrt genotypes, which are represented by combinations of five amino acids at positions 72-76, have been described in field isolates from various endemic regions. To identify novel mutant pfcrt genotypes and to reveal the genetic relatedness of pfcrt genotypes, a large-scale survey over a wide geographic area was performed.

Methods

Sequences for exon 2 in pfcrt, including known polymorphic sites at amino acid positions 72, 74, 75 and 76, were obtained from 256 P. falciparum isolates collected from eight endemic countries in Asia (Bangladesh, Cambodia, Lao P.D.R., the Philippines and Thailand), Melanesia (Papua New Guinea and Vanuatu) and Africa (Ghana). A haplotype network was constructed based on six microsatellite markers located -29 kb to 24 kb from pfcrt in order to examine the genetic relatedness among mutant pfcrt genotypes.

Results

In addition to wild type (CVMNK at positions 72-76), four mutant pfcrt were identified; CVIET, CVIDT, SVMNT and CVMNT (mutated amino acids underlined). Haplotype network revealed that there were only three mutant pfcrt lineages, originating in Indochina, Philippines and Melanesia. Importantly, the Indochina lineage contained two mutant pfcrt genotypes, CVIET (n = 95) and CVIDT (n = 14), indicating that CVIDT shares a common origin with CVIET. Similarly, one major haplotype in the Melanesian lineage contained two pfcrt genotypes; SVMNT (n = 71) and CVMNT (n = 3). In Africa, all mutant pfcrt genotypes were the CVIET of the Indochina lineage, probably resulting from the intercontinental migration of CQ resistance from Southeast Asia.

Conclusions

The number of CQ-mutant lineages observed in this study was identical to that found in previous studies. This supports the hypothesis that the emergence of novel CQ resistance is rare. However, in the mutant pfcrt genotypes, amino acid changes at positions 72, 74 and 75 appear to have recently been generated at least several times, producing distinct pfcrt mutant genotypes. The occurrence of new mutations flanking K76T may yield stronger resistance to CQ and/or a higher fitness than the original pfcrt mutant.

Differences in selective pressure on dhps and dhfr drug resistant mutations in western Kenya

Background

Understanding the origin and spread of mutations associated with drug resistance, especially in the context of combination therapy, will help guide strategies to halt and prevent the emergence of resistance. Unfortunately, studies have assessed these complex processes when resistance is already highly prevalent. Even further, information on the evolutionary dynamics leading to multidrug-resistant parasites is scattered and limited to areas with low or seasonal malaria transmission. This study describes the dynamics of strong selection for mutations conferring resistance against sulphadoxine-pyrimethamine (SP), a combination therapy, in western Kenya between 1992 and 1999, just before SP became first-line therapy (1999). Importantly, the study is based on longitudinal data, which allows for a comprehensive analysis that contrasts with previous cross-sectional studies carried out in other endemic regions.

Methods

This study used 236 blood samples collected between 1992 and 1999 in the Asembo Bay area of Kenya. Pyrosequencing was used to determine the alleles of dihydrofolate reductase (dhfr) and dihydropterote synthase (dhps) genes. Microsatellite alleles spanning 138 kb around dhfr and dhps, as well as, neutral markers spanning approximately 100 kb on chromosomes 2 and 3 were characterized.

Results

By 1992, the South-Asian dhfr triple mutant was already spreading, albeit in low frequency, in this holoendemic Kenyan population, prior to the use of SP as a first-line therapy. Additionally, dhfr triple mutant alleles that originated independently from the predominant Southeast Asian lineage were present in the sample set. Likewise, dhps double mutants were already present as early as 1992. There is evidence for soft selective sweeps of two dhfr mutant alleles and the possible emergence of a selective sweep of double mutant dhps alleles between 1992 and 1997. The longitudinal structure of the dataset allowed estimation of selection pressures on various dhfr and dhps mutants relative to each other based on a theoretical model tailored to P. falciparum. The data indicate that drug selection acted differently on the resistant alleles of dhfr and dhps, as evidenced by fitness differences. Thus a combination drug therapy such as SP, by itself, does not appear to select for "multidrug"-resistant parasites in areas with high recombination rate.

Conclusions

The complexity of these observations emphasizes the importance of population-based studies to evaluate the effects of strong drug selection on Plasmodium falciparum populations.

Malaria control in Malawi: current status and directions for the future

The last decade has seen an increase in investment and concerted efforts by the Malawi Ministry of Health and partners to control malaria disease. This report summarizes what is known about the burden of malaria and the strategies being implemented to control it in Malawi. Over the past 5 years, roll out of treatment and prevention efforts have been successful in the country, as demonstrated by increased use of insecticide treated nets, improved access to prompt and effective treatment and the initiation of pilot studies of indoor residual spraying. However, unlike other countries in the region, the recent data have not suggested a decrease in the burden of disease. We describe the environment in which the activities of Malawi's International Center for Excellence in Malaria Research (ICEMR) will be carried out and provide the rationale for the clinical, entomological and molecular studies. Our approach is to establish consistent, stainable data collection systems that are embedded within the public health sector. Through standardized and long-term studies of hosts, parasites and vectors, we hope to contribute to assessment of malaria disease burden, the appropriate application of interventions and policies and provide both the data collection and the health care infrastructure to ultimately eliminate the disease.

Differential genetic hitchhiking around mutant pfcrt alleles in the Indian Plasmodium falciparum population

OBJECTIVES

To study the origin and spread of the chloroquine-resistant Plasmodium falciparum population in the Indian subcontinent.

METHODS

Fourteen microsatellites spanning a ∼120 kb region, flanking the P. falciparum chloroquine resistance transporter (pfcrt) gene, were analysed in 185 parasite isolates.

RESULTS

The Indian P. falciparum population exhibited a selective valley of reduced genetic variation in the flanking microsatellites of the mutant pfcrt alleles (up to ±29 kb) as compared with the wild-type allele. This valley is much narrower than the ±200 kb valley reported from African and South-East Asian countries. The majority of the isolates showed asymmetry in the selective valley, where upstream microsatellites showed less genetic variation than the downstream microsatellites. Regional variation in the width and symmetry of the selective valley was noticed, which seems to be related to the number of pfcrt alleles present in the parasite population of a region. Forty-six different microsatellite haplotypes were observed among the P. falciparum isolates containing mutant pfcrt alleles. Parasite populations from different regions of mainland India shared microsatellite haplotypes between them, but they shared none with the isolates from the Andaman and Nicobar Islands, and vice versa. Indian isolates shared microsatellite haplotypes with the isolates from Papua New Guinea and Thailand.

CONCLUSIONS

With regard to chloroquine there is regional variation in the selection pressure on the P. falciparum population in India. These findings will help the regional implementation of drug policy in India's malaria control programme.

Genetic backgrounds of the Plasmodium falciparum chloroquine resistant transporter (pfcrt) alleles in Pakistan

Chloroquine (CQ) resistance in Plasmodium falciparum has been associated with point mutations in the P. falciparum CQ resistance transporter gene (pfcrt). Previous studies have shown 4-5 independent origins for CQ resistant pfcrt alleles globally, two in South America, one each in Southeast Asia, Papua New Guinea (PNG) and Philippines. In Asia, at least two different alleles corresponding to amino acids 72-76 (CVIET and SVMNT) have been found. The CVIET allele originated in Southeast Asia and then spread to Asia and Africa as well. The SVMNT allele, originating from PNG, has been found in India. This study was undertaken to investigate the genetic background of the CQ resistant pfcrt haplotypes in Pakistan. We genotyped microsatellite markers surrounding the pfcrt gene (six different markers at -12.3, -4.8, -1, 1.5, 3.9, 18.8 kb) in 114 clinical isolates of P. falciparum collected from different regions in Pakistan. Microsatellite analysis showed a significant reduction in genetic variation among the mutant SVMNT pfcrt alleles when compared to wild type alleles. The predominant SVMNT haplotype found in this study shared the same microsatellite haplotype found in both PNG and India. Two isolates with CVIET haplotypes showed similar microsatellite background to those found in Africa and Asia. In conclusion, this study suggests that CQ resistant SVMNT haplotypes in India and Pakistan have a common ancestral origin similar to that of Papua New Guinean isolates.

South American Plasmodium falciparum after the malaria eradication era: clonal population expansion and survival of the fittest hybrids

Malaria has reemerged in many regions where once it was nearly eliminated. Yet the source of these parasites, the process of repopulation, their population structure, and dynamics are ill defined. Peru was one of malaria eradication's successes, where Plasmodium falciparum was nearly eliminated for two decades. It reemerged in the 1990s. In the new era of malaria elimination, Peruvian P. falciparum is a model of malaria reinvasion. We investigated its population structure and drug resistance profiles. We hypothesized that only populations adapted to local ecological niches could expand and repopulate and originated as vestigial populations or recent introductions. We investigated the genetic structure (using microsatellites) and drug resistant genotypes of 220 parasites collected from patients immediately after peak epidemic expansion (1999-2000) from seven sites across the country. The majority of parasites could be grouped into five clonal lineages by networks and AMOVA. The distribution of clonal lineages and their drug sensitivity profiles suggested geographic structure. In 2001, artesunate combination therapy was introduced in Peru. We tested 62 parasites collected in 2006-2007 for changes in genetic structure. Clonal lineages had recombined under selection for the fittest parasites. Our findings illustrate that local adaptations in the post-eradication era have contributed to clonal lineage expansion. Within the shifting confluence of drug policy and malaria incidence, populations continue to evolve through genetic outcrossing influenced by antimalarial selection pressure. Understanding the population substructure of P. falciparum has implications for vaccine, drug, and epidemiologic studies, including monitoring malaria during and after the elimination phase.

Recent Progress in Functional Genomic Research in Plasmodium falciparum

With the completion and near completion of many malaria parasite genome-sequencing projects, efforts are now being directed to a better understanding of gene functions and to the discovery of vaccine and drug targets. Inter- and intraspecies comparisons of the parasite genomes will provide invaluable insights into parasite evolution, virulence, drug resistance, and immune invasion. Genome-wide searches for loci under various selection pressures may lead to discovery of genes conferring drug resistance or encoding for protective antigens. In addition, the Plasmodium falciparum genome sequence provides the basis for the development of various microarrays to monitor gene expression and to detect nucleotide substitution and deletion/amplification. Genome-wide profiling of the parasite proteome, chromatin modification, and nucleosome position also depend on availability of the parasite genome. In this brief review, we will highlight some recent advances and studies in characterizing gene function and related phenotype in P. falciparum that were made possible by the genome sequence, particularly the development of a genome-wide diversity map and various high-throughput genotyping methods for genome-wide association studies (GWAS).

Multiple origins of Plasmodium falciparum dihydropteroate synthetase mutant alleles associated with sulfadoxine resistance in India

With the spread of chloroquine (CQ)-resistant malaria in India, sulfadoxine-pyrimethamine (SP) alone or in combination with artesunate is used as an alternative antimalarial drug. Due to continuous drug pressure, the Plasmodium falciparum parasite is exhibiting resistance to antifolates because of mutations in candidate genes dihydrofolate reductase (dhfr) and dihydropteroate synthetase (dhps). Our earlier study on flanking microsatellite markers of dhfr mutant alleles from India had shown a single origin of the pyrimethamine resistance and some minor haplotypes which shared haplotypes with Southeast Asian (Thailand) strains. In the present study, we have analyzed 193 of these Indian P. falciparum isolates for 15 microsatellite loci around dhps to investigate the genetic lineages of the mutant dhps alleles in different parts of the country. Eighty-one of these samples had mutant dhps alleles, of which 62 were from Andaman and Nicobar Islands and the remaining 19 were from mainland India. Of 112 isolates with a wild-type dhps allele, 109 were from mainland India and only 3 were from Andaman and Nicobar Islands. Consistent with the model of selection, the mean expected heterozygosity (H(e)) around mutant dhps alleles (H(e) = 0.55; n = 81) associated with sulfadoxine resistance was lower (P ≤ 0.05) than the mean H(e) around the wild-type dhps allele (H(e) = 0.80; n = 112). There was more genetic diversity in flanking microsatellites of dhps than dhfr among these isolates, which confirms the assertion that dhps mutations are at a very early stage of fixation in the parasite population. Microsatellite haplotypes around various mutant dhps alleles suggest that the resistant dhps alleles have multiple independent origins in India, especially in Andaman and Nicobar Islands. Determining the genetic lineages of the resistant dhps alleles on Andaman and Nicobar Islands and mainland India is significant, given the role of Asia in the intercontinental spread of chloroquine- and pyrimethamine-resistant parasites in the past.

Patterns of chloroquine use and resistance in sub-Saharan Africa: a systematic review of household survey and molecular data

BACKGROUND

As a result of widespread chloroquine and sulphadoxine-pyrimethamine (SP) resistance, 90% of sub-Saharan African countries had adopted policies of artemisinin-based combination therapy (ACT) for treatment of uncomplicated malaria by 2007. In Malawi, cessation of chloroquine use was followed by the re-emergence of chloroquine-susceptible malaria. It was expected that introduction of ACT would lead to a return in chloroquine susceptibility throughout Africa, but this has not yet widely occurred. This observation suggests that there is continuing use of ineffective anti-malarials in Africa and that persistent chloroquine-resistant malaria is due to ongoing drug pressure despite national policy changes.

METHODS

To estimate drug use on a national level, 2006-2007 Demographic Health Survey and Multiple Indicator Cluster Survey data from 21 African countries were analysed. Resistance data were compiled by systematic review of the published literature on the prevalence of the Plasmodium falciparum chloroquine resistance transporter polymorphism at codon 76, which causes chloroquine resistance.

RESULTS

Chloroquine was the most common anti-malarial used according to surveys from 14 of 21 countries analysed, predominantly in West Africa. SP was most commonly reported in two of 21 countries. Among eight countries with longitudinal molecular resistance data, the four countries where the highest proportion of children treated for fever received chloroquine (Uganda, Burkina Faso, Guinea Bissau, and Mali) also showed no significant declines in the prevalence of chloroquine-resistant infections. The three countries with low or decreasing chloroquine use among children who reported fever treatment (Malawi, Kenya, and Tanzania) had statistically significant declines in the prevalence of chloroquine resistance.

CONCLUSIONS

This study demonstrates that in 2006-2007, chloroquine and SP continued to be used at high rates in many African countries. In countries reporting sustained chloroquine use, chloroquine-resistant malaria persists. In contrast, a low level of estimated chloroquine use is associated with a declining prevalence of chloroquine resistance.

Selective sweeps and genetic lineages of Plasmodium falciparum drug -resistant alleles in Ghana

BACKGROUND

In 2005, Ghana adopted artemisinin-based combination therapy (ACT) for primary treatment of falciparum malaria. A comprehensive study of the drug-resistance-associated mutations and their genetic lineages will lead to a better understanding of the evolution of antimalarial drug resistance in this region.

METHODS

The pfcrt, pfmdr1, dhps, and dhfr mutations associated with chloroquine (CQ) and sulfadoxine-pyrimethamine (SP) resistance and the microsatellite loci flanking these genes were genotyped in Plasmodium falciparum isolates from Ghana.

RESULTS

The prevalence of mutations associated with both CQ and SP resistance was high in Ghana. However, we observed a decrease in prevalence of the pfcrt K76T mutation in northern Ghana after the change in drug policy from CQ to ACT. Analysis of genetic diversity and differentiation at microsatellite loci flanking all 4 genes indicated that they have been under strong selection, because of CQ and SP use. The triple-mutant pfcrt and dhfr alleles in Ghana were derived from Southeast Asia, whereas the double-mutant dhfr, dhps, and pfmdr1 alleles were of African lineage.

CONCLUSION

Because of the possible role of pfmdr1 in amodiaquine and mefloquine resistance, demonstrating selection on pfmdr1 and defining lineages of resistant alleles in an African population holds great importance.

Tracking origins and spread of sulfadoxine-resistant Plasmodium falciparum dhps alleles in Thailand

The emergence and spread of drug-resistant Plasmodium falciparum have been a major impediment for the control of malaria worldwide. Earlier studies have shown that similar to chloroquine (CQ) resistance, high levels of pyrimethamine resistance in P. falciparum originated independently 4 to 5 times globally, including one origin at the Thailand-Cambodia border. In this study we describe the origins and spread of sulfadoxine-resistance-conferring dihydropteroate synthase (dhps) alleles in Thailand. The dhps mutations and flanking microsatellite loci were genotyped for P. falciparum isolates collected from 11 Thai provinces along the Burma, Cambodia, and Malaysia borders. Results indicated that resistant dhps alleles were fixed in Thailand, predominantly being the SGEGA, AGEAA, and SGNGA triple mutants and the AGKAA double mutant (mutated codons are underlined). These alleles had different geographical distributions. The SGEGA alleles were found mostly at the Burma border, while the SGNGA alleles occurred mainly at the Cambodia border and nearby provinces. Microsatellite data suggested that there were two major genetic lineages of the triple mutants in Thailand, one common for SGEGA/SGNGA alleles and another one independent for AGEAA. Importantly, the newly reported SGNGA alleles possibly originated at the Thailand-Cambodia border. All parasites in the Yala province (Malaysia border) had AGKAA alleles with almost identical flanking microsatellites haplotypes. They were also identical at putatively neutral loci on chromosomes 2 and 3, suggesting a clonal nature of the parasite population in Yala. In summary, this study suggests multiple and independent origins of resistant dhps alleles in Thailand.

Whole genome re-sequencing identifies a mutation in an ABC transporter (mdr2) in a Plasmodium chabaudi clone with altered susceptibility to antifolate drugs

In malaria parasites, mutations in two genes of folate biosynthesis encoding dihydrofolate reductase (dhfr) and dihydropteroate synthase (dhps) modify responses to antifolate therapies which target these enzymes. However, the involvement of other genes which modify the availability of exogenous folate, for example, has been proposed. Here, we used short-read whole-genome re-sequencing to determine the mutations in a clone of the rodent malaria parasite, Plasmodium chabaudi, which has altered susceptibility to both sulphadoxine and pyrimethamine. This clone bears a previously identified S106N mutation in dhfr and no mutation in dhps. Instead, three additional point mutations in genes on chromosomes 2, 13 and 14 were identified. The mutated gene on chromosome 13 (mdr2 K392Q) encodes an ABC transporter. Because Quantitative Trait Locus analysis previously indicated an association of genetic markers on chromosome 13 with responses to individual and combined antifolates, MDR2 is proposed to modulate antifolate responses, possibly mediated by the transport of folate intermediates.

Return of chloroquine-susceptible falciparum malaria in Malawi was a reexpansion of diverse susceptible parasites

The spread of drug-resistant Plasmodium falciparum malaria has been a major impediment to malaria control and threatens prospects for elimination. We recently demonstrated the return of chloroquine-susceptible malaria in Malawi after chloroquine use was abandoned. In this study, we trace the origins of chloroquine-resistant and chloroquine-susceptible parasites in Malawi by sequencing the P. falciparum chloroquine resistance transporter gene (pfcrt) and by genotyping microsatellites flanking this gene in isolates from infections that occurred in Malawi from 1992 through 2005. Malaria parasites from 2005 harbored the expected wild-type pfcrt haplotype associated with chloroquine susceptibility and have maintained high levels of diversity without linkage disequilibrium, which suggests that the return of chloroquine susceptibility is not the result of a back mutation in a formerly resistant parasite or a new selective sweep. Chloroquine-susceptible parasites that predominate in Malawi likely represent a reexpansion of the susceptible parasites that survived in the population despite widespread drug pressure in the region.

A microsatellite-based analysis for the detection of selection on BTA1 and BTA20 in northern Eurasian cattle (Bos taurus) populations

BACKGROUND

Microsatellites surrounding functionally important candidate genes or quantitative trait loci have received attention as proxy measures of polymorphism level at the candidate loci themselves. In cattle, selection for economically important traits is a long-term strategy and it has been reported that microsatellites are linked to these important loci.

METHODS

We have investigated the variation of seven microsatellites on BTA1 (Bos taurus autosome 1) and 16 on BTA20, using bovine populations of typical production types and horn status in northern Eurasia. Genetic variability of these loci and linkage disequilibrium among these loci were compared with those of 28 microsatellites on other bovine chromosomes. Four different tests were applied to detect molecular signatures of selection.

RESULTS

No marked difference in locus variability was found between microsatellites on BTA1, BTA20 and the other chromosomes in terms of different diversity indices. Average D' values of pairwise syntenic markers (0.32 and 0.28 across BTA 1 and BTA20 respectively) were significantly (P < 0.05) higher than for non-syntenic markers (0.15). The Ewens-Watterson test, the Beaumont and Nichol's modified frequentist test and the Bayesian FST-test indicated elevated or decreased genetic differentiation, at SOD1 and AGLA17 markers respectively, deviating significantly (P < 0.05) from neutral expectations. Furthermore, lnRV, lnRH and lnRtheta' statistics were used for the pairwise population comparison tests and were significantly less variable in one population relative to the other, providing additional evidence of selection signatures for two of the 51 loci. Moreover, the three Finnish native populations showed evidence of subpopulation divergence at SOD1 and AGLA17. Our data also indicate significant intergenic linkage disequilibrium around the candidate loci and suggest that hitchhiking selection has played a role in shaping the pattern of observed linkage disequilibrium.

CONCLUSION

Hitchhiking due to tight linkage with alleles at candidate genes, e.g. the POLL gene, is a possible explanation for this pattern. The potential impact of selective breeding by man on cattle populations is discussed in the context of selection effects. Our results also suggest that a practical approach to detect loci under selection is to simultaneously apply multiple neutrality tests based on different assumptions and estimations.

A high incidence of selection on physiologically important genes in the three-spined stickleback, Gasterosteus aculeatus

Genome scan approaches to detect footprints of directional selection in the genomes of wild animal and plant populations have become popular tools to study local adaptation and speciation at the molecular level. Most studies thus far have used random molecular markers and found footprints of directional selection at, on average, 5% (range: 1-15%) of the examined loci. We focused on physiologically important genes that exhibit transcriptional responses to specific environmental or developmental conditions and assessed if these genes have been subject to directional selection and are responsible for local adaptation in the three-spined stickleback (Gasterosteus aculeatus). Using microsatellite markers located within or closely linked to (<6 kb) target genes, we investigated footprints of directional selection for 157 genes with known physiological functions in three marine and six freshwater populations. A high incidence (16.6%) of footprints of directional selection for these genes was revealed by four different outlier tests. In a subset of four populations screened with both physiologically important and random genes, footprints of directional selection were more frequent in physiologically important genes (13.4%) as compared with random genes (2.4%). In general, our findings indicate strong selective pressures on physiologically important genes, suggesting that these genes have significant functions in evolutionary adaptation to environmental heterogeneity.