Publications

We explored the potential of pooled sequencing to swiftly and economically identify selective sweeps due to emerging artemisinin (ART) resistance in a South-East Asian malaria parasite population. ART resistance is defined by slow parasite clearance from the blood of ART-treated patients and mutations in the kelch gene (chr. 13) have been strongly implicated to play a role. We constructed triplicate pools of 70 slow-clearing (resistant) and 70 fast-clearing (sensitive) infections collected from the Thai-Myanmar border and sequenced these to high (∼ 150-fold) read depth. Allele frequency estimates from pools showed almost perfect correlation (Lin's concordance = 0.98) with allele frequencies at 93 single nucleotide polymorphisms measured directly from individual infections, giving us confidence in the accuracy of this approach. By mapping genome-wide divergence (FST) between pools of drug-resistant and drug-sensitive parasites, we identified two large (>150 kb) regions (on chrs. 13 and 14) and 17 smaller candidate genome regions. To identify individual genes within these genome regions, we resequenced an additional 38 parasite genomes (16 slow and 22 fast-clearing) and performed rare variant association tests. These confirmed kelch as a major molecular marker for ART resistance (P = 6.03 × 10(-6)). This two-tier approach is powerful because pooled sequencing rapidly narrows down genome regions of interest, while targeted rare variant association testing within these regions can pinpoint the genetic basis of resistance. We show that our approach is robust to recurrent mutation and the generation of soft selective sweeps, which are predicted to be common in pathogen populations with large effective population sizes, and may confound more traditional gene mapping approaches.

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.