Publications

BACKGROUND: Relapses originating from hypnozoites are characteristic of Plasmodium vivax infections. Thus, reappearance of parasitemia after treatment can result from relapse, recrudescence, or reinfection. It has been assumed that parasites causing relapse would be a subset of the parasites that caused the primary infection.

METHODS: Paired samples were collected before initiation of antimalarial treatment and at recurrence of parasitemia from 149 patients with vivax malaria in Thailand (n=36), where reinfection could be excluded, and during field studies in Myanmar (n=75) and India (n=38).

RESULTS: Combined genetic data from 2 genotyping approaches showed that novel P. vivax populations were present in the majority of patients with recurrent infection (107 [72%] of 149 patients overall [78% of patients in Thailand, 75% of patients in Myanmar {Burma}, and 63% of patients in India]). In 61% of the Thai and Burmese patients and in 55% of the Indian patients, the recurrent infections contained none of the parasite genotypes that caused the acute infection.

CONCLUSIONS: The P. vivax populations emerging from hypnozoites commonly differ from the populations that caused the acute episode. Activation of heterologous hypnozoite populations is the most common cause of first relapse in patients with vivax malaria.

Molecular markers provide a rapid and relatively inexpensive approach for assessing antimalarial drug susceptibility. We collected 884 Plasmodium falciparum-infected blood samples from 17 Lao provinces. Each sample was genotyped for 11 codons in the chloroquine resistance transporter (pfcrt), dihydrofolate reductase (pfdhfr), and dihydropteroate synthase (pfdhps) genes. The samples included 227 collected from patients recruited to clinical trials. The pfcrt K76T mutation was an excellent predictor of treatment failure for both chloroquine and chloroquine plus sulfadoxine-pyrimethamine, and mutations in both pfdhfr and pfdhps were predictive of sulfadoxine-pyrimethamine treatment failure. In multivariate analysis, the presence of the pfdhfr triple mutation (51 + 59 + 108) was strongly and independently correlated with sulfadoxine-pyrimethamine failure (odds ratio = 9.1, 95% confidence interval = 1.4-60.2, P = 0.017). Considerable geographic heterogeneity in allele frequencies occurred at all three loci with lower frequencies of mutant alleles in southern than in northern Laos. These findings suggest that chloroquine and sulfadoxine-pyrimethamine are no longer viable therapy in this country.

BACKGROUND: Our study examined the relative contributions of host, pharmacokinetic, and parasitological factors in determining the therapeutic response to artemether-lumefantrine (AL).

METHODS: On the northwest border of Thailand, patients with uncomplicated Plasmodium falciparum malaria were enrolled in prospective studies of AL treatment (4- or 6-dose regimens) and followed up for 42 days. Plasma lumefantrine concentrations were measured by high performance liquid chromatography; malaria parasite pfmdr1 copy number was quantified using a real-time polymerase chain reaction assay (PCR), and in vitro drug susceptibility was tested.

RESULTS: All treatments resulted in a rapid clinical response and were well tolerated. PCR-corrected failure rates at day 42 were 13% (95% confidence interval [CI], 9.6%-17%) for the 4-dose regimen and 3.2% (95% CI, 1.8%-4.6%) for the 6-dose regimen. Increased pfmdr1 copy number was associated with a 2-fold (95% CI, 1.8-2.4-fold) increase in lumefantrine inhibitory concentration(50) (P=.001) and an adjusted hazard ratio for risk of treatment failure following completion of a 4-dose regimen, but not a 6-dose regimen, of 4.0 (95% CI, 1.4-11; P=.008). Patients who had lumefantrine levels below 175 ng/mL on day 7 were more likely to experience recrudescence by day 42 (adjusted hazard ratio, 17; 95% CI, 5.5-53), allowing prediction of treatment failure with 75% sensitivity and 84% specificity. The 6-dose regimen ensured that therapeutic levels were achieved in 91% of treated patients.

CONCLUSIONS: The lumefantrine plasma concentration profile is the main determinant of efficacy of artemether-lumefantrine. Amplification in pfmdr1 determines lumefantrine susceptibility and, therefore, treatment responses when plasma lumefantrine levels are subtherapeutic.

Mu et al. (Mu, J., M. T. Ferdig, X. Feng, D. A. Joy, J. Duan, T. Furuya, G. Subramanian, L. Aravind, R. A. Cooper, J. C. Wootton, M. Xiong, and X. Z. Su, Mol. Microbiol. 49:977-989, 2003) recently reported exciting associations between nine new candidate transporter genes and in vitro resistance to chloroquine (CQ) and quinine (QN), with six of these loci showing association with CQ or QN in a southeast Asian population sample. We replicated and extended this work by examining polymorphisms in these genes and in vitro resistance to eight drugs in parasites collected from the Thailand-Burma border. To minimize problems of multiple testing, we used a two-phase study design, while to minimize problems caused by population structure, we analyzed parasite isolates collected from a single clinic. We first examined associations between genotype and drug response in 108 unique single-clone parasite isolates. We found strong associations between single nucleotide polymorphisms in pfmdr and mefloquine (MFQ), artesunate (AS), and lumefantrine (LUM) response. We also observed associations between an ABC transporter (G7) and response to QN and AS and between another ABC transporter (G49) and response to dihydro-artemisinin (DHA). We reexamined significant associations in an independent sample of 199 unique single-clone infections from the same location. The significant associations with pfmdr-1042 detected in the first survey remained. However, with the exception of the G7-artesunate association, all other associations observed with the nine new candidate transporters disappeared. We also examined linkage disequilibrium (LD) between markers and phenotypic correlations between drug responses. We found minimal LD between genes. Furthermore, we found no correlation between chloroquine and quinine responses, although we did find expected strong correlations between MFQ, QN, AS, DHA, and LUM. To conclude, we found no evidence for an association between 8/9 candidate genes and response to eight different antimalarial drugs. However, the consistent association observed between a 3-bp indel in G7 and AS response merits further investigation.

Neutral mutations may hitchhike to high frequency when they are situated close to sites under positive selection, generating local reductions in genetic diversity. This process is thought to be an important determinant of levels of genomic variation in natural populations. The size of genome regions affected by genetic hitchhiking is expected to be dependent on the strength of selection, but there is little empirical data supporting this prediction. Here, we compare microsatellite variation around two drug resistance genes (chloroquine resistance transporter (pfcrt), chromosome 7, and dihydrofolate reductase (dhfr), chromosome 4) in malaria parasite populations exposed to strong (Thailand) or weak selection (Laos) by anti-malarial drugs. In each population, we examined the point mutations underlying resistance and length variation at 22 (chromosome 4) or 25 (chromosome 7) microsatellite markers across these chromosomes. All parasites from Thailand carried the K76T mutation in pfcrt conferring resistance to chloroquine (CQ) and 2-4 mutations in dhfr conferring resistance to pyrimethamine. By contrast, we found both wild-type and resistant alleles at both genes in Laos. There were dramatic differences in the extent of hitchhiking in the two countries. The size of genome regions affected was smaller in Laos than in Thailand. We observed significant reduction in variation relative to sensitive parasites for 34-64 kb (2-4 cM) in Laos on chromosome 4, compared with 98-137 kb (6-8 cM) in Thailand. Similarly, on chromosome 7, we observed reduced variation for 34-69 kb (2-4 cM) around pfcrt in Laos, but for 195-268 kb (11-16 cM) in Thailand. Reduction in genetic variation was also less extreme in Laos than in Thailand. Most loci were monomorphic in a 12 kb region surrounding both genes on resistant chromosomes from Thailand, whereas in Laos, even loci immediately proximal to selective sites showed some variation on resistant chromosomes. Finally, linkage disequilibrium (LD) decayed more rapidly around resistant pfcrt and dhfr alleles from Laos than from Thailand. These results demonstrate that different realizations of the same selective sweeps may vary considerably in size and shape, in a manner broadly consistent with selection history. From a practical perspective, genomic regions containing resistance genes may be most effectively located by genome-wide association in populations exposed to strong drug selection. However, the lower levels of LD surrounding resistance alleles in populations under weak selection may simplify identification of functional mutations.

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.

Understanding the frequency with which new resistance alleles arise and their subsequent patterns of spread is critical to our attempts to manage drug resistance in parasite populations. We review recent molecular evolutionary studies utilizing marker loci situated close to resistance loci on the Plasmodium falciparum genome that have given surprising insights into the origins and spread of drug resistance loci. We discuss possible reasons for the patterns observed, and highlight the implications of these results for resistance management. In particular, we show that many resistance mutations have rather few independent origins. De novo mutation appears to be less important than migration for introducing resistance alleles into parasite populations. Attempts to manage drug resistance will be of limited effectiveness unless this is taken into account.

Here we present molecular evidence demonstrating that malaria parasites bearing high-level pyrimethamine resistance originally arrived in Africa from southeast Asia. The resistance alleles carried by these migrants are now spreading across Africa at an alarming rate, signaling the end of affordable malaria treatment and presenting sub-Saharan Africa with a public health crisis.

When alleles conferring drug resistance spread through a population of malaria parasites, they leave characteristic "scars" in the parasite genome. Flanking neutral polymorphisms "hitchhike" to high frequency with the resistance mutation, generating deep valleys of reduced variation and broad swathes of elevated linkage disequilibrium around the resistance locus. We can systematically search the genome for these scars by genotyping polymorphic marker loci at intervals throughout the genome of P. falciparum, and use them as signposts for locating drug resistance genes. In this review I outline the rational behind this approach to genetic mapping. I describe key features of P. falciparum population biology, such as recombination rate, inbreeding, and selection intensity that influence the size of genomic regions affected by selection and the choice of study population. I discuss suitable genetic markers, study designs, and statistical approaches to data analysis. Finally, to demonstrate the utility of the approach I describe two proof-of-principle studies documenting patterns of genetic variability around known drug resistance genes.

Sexual transmission occurs commonly in microparasites such as viruses and bacteria, but this is an unusual transmission route for macroparasites. Here we present evidence which suggests that a nematode parasite of Wood Mice (Apodemus sylvaticus) may be sexually transmitted and we have classified the nematode using molecular data. Wood Mice were collected annually in the course of work on their reproductive physiology. Larval nematodes were found in the epididymides of 19.6% of males. It seems likely that they would be transmitted to females at ejaculation. To identify these larval nematodes, which we were unable to do using morphological features, we sequenced the 18S rDNA. Sequence comparisons with the molecular phylogeny of Blaxter et al. (1998) demonstrated that they were bursate nematodes (Order Strongylida). The relationships between strongylid taxa were poorly resolved by 18S rDNA. However, both distance and parsimony analyses grouped the nematode with the superfamily Metastrongylidea in a clade containing Filaroides and Angiostrongylus sp. Importantly, the sequences were distinct from those of Heligmosomoides polygyrus and Angiostrongylus dujardini, two common strongylid nematodes of Apodemus. We were therefore unable positively to identify these worms by matching their sequences with those from morphologically identifiable adult strongylid nematodes infecting Apodemus. These results demonstrate that an as yet unidentified strongylid is quite commonly found in large numbers in the male reproductive tract of Wood Mice. Further work is required to understand the biology and transmission dynamics of this interesting system.