Publications

We describe 35 microsatellite markers from the human parasitic nematode Ascaris lumbricoides. We found 7 sex-linked markers and demonstrate that 26 autosomal loci can be scored reliably. These markers have high genetic variability and provide the tools to address multiple questions concerning the epidemiology, fine-scale genetic structure, host specificity, and mating systems of this parasite.

Knowledge of cross-transmission and hybridization between parasites of humans and reservoir hosts is critical for understanding the evolution of the parasite and for implementing control programmes. There is now a consensus that populations of pig and human Ascaris (roundworms) show significant genetic subdivision. However, it is unclear whether this has resulted from a single or multiple host shift(s). Furthermore, previous molecular data have not been sufficient to determine whether sympatric populations of human and pig Ascaris can exchange genes. To disentangle patterns of host colonization and hybridization, we used 23 microsatellite loci to conduct Bayesian clustering analyses of individual worms collected from pigs and humans. We observed strong differentiation between populations which was primarily driven by geography, with secondary differentiation resulting from host affiliation within locations. This pattern is consistent with multiple host colonization events. However, there is low support for the short internal branches of the dendrograms. In part, the relationships among clusters may result from current hybridization among sympatric human and pig roundworms. Indeed, congruence in three Bayesian methods indicated that 4 and 7% of roundworms sampled from Guatemala and China, respectively, were hybrids. These results indicate that there is contemporary cross-transmission between populations of human and pig Ascaris.

Recent drug trials in Laos have shown high levels of Plasmodium falciparum resistance to chloroquine, but there are no published data on in vitro antimalarial drug susceptibility. We used the double-site enzyme-linked pLDH immunodetection (DELI) assay to estimate the in vitro antimalarial drug susceptibility of 108 fresh P. falciparum isolates from southern Laos. The geometric mean (95% confidence interval) 50% inhibitory concentration values (nmol/L) were 152.4 (123.8-187.6) for chloroquine, 679.8 (533.8-863.0) for quinine, 45.9 (37.9-55.7) for mefloquine, 5.0 (4.4-6.4) for artesunate, 6.3 (4.5-8.9) for dihydroartemisinin, and 59.1 (46.4-75.3) for lumefantrine. The proportion of isolates defined as resistant were 65%, 40%, and 8% for chloroquine, quinine, and mefloquine, respectively. Of 53 isolates genotyped for the pfcrt T76K chloroquine-resistance mutation, 48 (91%) were mutants. P. falciparum in Laos is multi-drug resistant; antimalarial immunity resulting from the use of ineffective chloroquine before 2005 probably contributes significantly to the therapeutic responses in clinical trials.

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.