Publications

BACKGROUND: Identification of parasite genes that underlie traits such as drug resistance and host specificity is challenging using classical linkage mapping approaches. Extreme QTL (X-QTL) methods, originally developed by rodent malaria and yeast researchers, promise to increase the power and simplify logistics of linkage mapping in experimental crosses of schistosomes (or other helminth parasites), because many 1000s of progeny can be analysed, phenotyping is not required, and progeny pools rather than individuals are genotyped. We explored the utility of this method for mapping a drug resistance gene in the human parasitic fluke Schistosoma mansoni.

RESULTS: We staged a genetic cross between oxamniquine sensitive and resistant parasites, then between two F1 progeny, to generate multiple F2 progeny. One group of F2s infecting hamsters was treated with oxamniquine, while a second group was left untreated. We used exome capture to reduce the size of the genome (from 363 Mb to 15 Mb) and exomes from pooled F2 progeny (treated males, untreated males, treated females, untreated females) and the two parent parasites were sequenced to high read depth (mean = 95-366×) and allele frequencies at 14,489 variants compared. We observed dramatic enrichment of alleles from the resistant parent in a small region of chromosome 6 in drug-treated male and female pools (combined analysis: Z = 11.07, p = 8.74 × 10(-29)). This region contains Smp_089320 a gene encoding a sulfotransferase recently implicated in oxamniquine resistance using classical linkage mapping methods.

CONCLUSIONS: These results (a) demonstrate the utility of exome capture for generating reduced representation libraries in Schistosoma mansoni, and (b) provide proof-of-principle that X-QTL methods can be successfully applied to an important human helminth. The combination of these methods will simplify linkage analysis of biomedically or biologically important traits in this parasite.

The analysis of in vitro anti-malarial drug susceptibility testing is vulnerable to the effects of different statistical approaches and selection biases. These confounding factors were assessed with respect to pfmdr1 gene mutation and amplification in 490 clinical isolates. Two statistical approaches for estimating the drug concentration associated with 50% effect (EC50 ) were compared: the commonly used standard two-stage (STS) method, and nonlinear mixed-effects modelling. The in vitro concentration-effect relationships for, chloroquine, mefloquine, lumefantrine and artesunate, were derived from clinical isolates obtained from patients on the western border of Thailand. All isolates were genotyped for polymorphisms in the pfmdr1 gene. The EC50 estimates were similar for the two statistical approaches but 15-28% of isolates in the STS method had a high coefficient of variation (>15%) for individual estimates of EC50 and these isolates had EC50 values that were 32 to 66% higher than isolates derived with more precision. In total 41% (202/490) of isolates had amplification of pfmdr1 and single nucleotide polymorphisms were found in 50 (10%). Pfmdr1 amplification was associated with an increase in EC50 for mefloquine (139% relative increase in EC50 for 2 copies, 188% for 3+ copies), lumefantrine (82% and 75% for 2 and 3+ copies respectively) and artesunate (63% and 127% for 2 and 3+ copies respectively). In contrast pfmdr1 mutation at codons 86 or 1042 were associated with an increase in chloroquine EC50 (44-48%). Sample size calculations showed that to demonstrate an EC50 shift of 50% or more with 80% power if the prevalence was 10% would require 430 isolates and 245 isolates if the prevalence was 20%. In conclusion, although nonlinear mixed-effects modelling did not demonstrate any major advantage for determining estimates of anti-malarial drug susceptibility, the method includes all isolates, thereby, potentially improving confirmation of candidate molecular markers of anti-malarial drug susceptibility.

BACKGROUND: Resistance to chloroquine and antifolate drugs has evolved independently in South America, suggesting that genotype - phenotype studies aimed at understanding the genetic basis of resistance to these and other drugs should be conducted in this continent. This research was conducted to better understand the population structure of Colombian Plasmodium falciparum in preparation for such studies.

RESULTS: A set of 384 SNPs were genotyped in blood spot DNA samples from 447 P. falciparum infected subjects collected over a ten year period from four provinces of the Colombian Pacific coast to evaluate clonality, population structure and linkage disequilibrium (LD). Most infections (81%) contained a single predominant clone. These clustered into 136 multilocus genotypes (MLGs), with 32% of MLGs recovered from multiple (2 - 28) independent subjects. We observed extremely low genotypic richness (R = 0.42) and long persistence of MLGs through time (median = 537 days, range = 1 - 2,997 days). There was a high probability (>5%) of sampling parasites from the same MLG in different subjects within 28 days, suggesting caution is needed when using genotyping methods to assess treatment success in clinical drug trials. Panmixia was rejected as four well differentiated subpopulations (FST = 0.084 - 0.279) were identified. These occurred sympatrically but varied in frequency within the four provinces. Linkage disequilibrium (LD) decayed more rapidly (r2 = 0.17 for markers <10 kb apart) than observed previously in South American samples.

CONCLUSIONS: We conclude that Colombian populations have several advantages for association studies, because multiple clone infections are uncommon and LD decays over the scale of one or a few genes. However, the extensive population structure and low genotype richness will need to be accounted for when designing and analyzing association studies.

Accurate measurement of malaria parasite clearance rates (CRs) following artemisinin (ART) treatment is critical for resistance surveillance and research, and various CR metrics are currently used. We measured 13 CR metrics in 1472 ART-treated hyperparasitemia infections for which 6-hour parasite counts and parasite genotypes (93 single nucleotide polymorphisms [SNPs]) were available. We used heritability to evaluate the performance of each metric. Heritability ranged from 0.06 ± 0.06 (SD) for 50% parasite clearance times to 0.67 ± 0.04 (SD) for clearance half-lives estimated from 6-hour parasite counts. These results identify the measures that should be avoided and show that reliable clearance measures can be obtained with abbreviated monitoring protocols.

BACKGROUND: Artemisinin-resistant Plasmodium falciparum has been reported in Pailin, western Cambodia, detected as a slow parasite clearance rate in vivo. Emergence of this phenotype in western Thailand and possibly elsewhere threatens to compromise the effectiveness of all artemisinin-based combination therapies. Parasite genetics is associated with parasite clearance rate but does not account for all variation. We investigated contributions of both parasite genetics and host factors to the artemisinin-resistance phenotype in Pursat, western Cambodia.

METHODS: Between June 19 and Nov 28, 2009, and June 26 and Dec 6, 2010, we enrolled patients aged 10 years or older with uncomplicated falciparum malaria, a density of asexual parasites of at least 10,000 per μL of whole blood, no symptoms or signs of severe malaria, no other cause of febrile illness, and no chronic illness. We gave participants 4 mg/kg artesunate at 0, 24, and 48 h, 15 mg/kg mefloquine at 72 h, and 10 mg/kg mefloquine at 96 h. We assessed parasite density on thick blood films every 6 h until undetectable. The parasite clearance half-life was calculated from the parasite clearance curve. We genotyped parasites with 18 microsatellite markers and patients for haemoglobin E, α-thalassaemia, and a mutation of G6PD, which encodes glucose-6-phosphate dehydrogenase. To account for the possible effects of acquired immunity on half-life, we used three surrogates for increased likelihood of exposure to P falciparum: age, sex, and place of residence. This study is registered with ClinicalTrials.gov, number NCT00341003.

FINDINGS: We assessed 3504 individuals from all six districts of Pursat province seeking treatment for malaria symptoms. We enrolled 168 patients with falciparum malaria who met inclusion criteria. The geometric mean half-life was 5·85 h (95% CI 5·54-6·18) in Pursat, similar to that reported in Pailin (p=0·109). We identified two genetically different parasite clone groups: parasite group 1 (PG1) and parasite group 2 (PG2). Non-significant increases in parasite clearance half-life were seen in patients with haemoglobin E (0·55 h; p=0·078), those of male sex (0·96 h; p=0·064), and in 2010 (0·68 h; p=0·068); PG1 was associated with a significant increase (0·79 h; p=0·033). The mean parasite heritability of half-life was 0·40 (SD 0·17).

INTERPRETATION: Heritable artemisinin resistance is established in a second Cambodian province. To accurately identify parasites that are intrinsically susceptible or resistant to artemisinins, future studies should explore the effect of erythrocyte polymorphisms and specific immune responses on half-life variation.

FUNDING: Division of Intramural Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health.