Publications

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.

Schistosomiasis is one of the most prevalent parasitic diseases, affecting millions of people in developing countries. Amongst the human-infective species, Schistosoma mansoni is also the most commonly used in the laboratory and here we present the systematic improvement of its draft genome. We used Sanger capillary and deep-coverage Illumina sequencing from clonal worms to upgrade the highly fragmented draft 380 Mb genome to one with only 885 scaffolds and more than 81% of the bases organised into chromosomes. We have also used transcriptome sequencing (RNA-seq) from four time points in the parasite's life cycle to refine gene predictions and profile their expression. More than 45% of predicted genes have been extensively modified and the total number has been reduced from 11,807 to 10,852. Using the new version of the genome, we identified trans-splicing events occurring in at least 11% of genes and identified clear cases where it is used to resolve polycistronic transcripts. We have produced a high-resolution map of temporal changes in expression for 9,535 genes, covering an unprecedented dynamic range for this organism. All of these data have been consolidated into a searchable format within the GeneDB (www.genedb.org) and SchistoDB (www.schistodb.net) databases. With further transcriptional profiling and genome sequencing increasingly accessible, the upgraded genome will form a fundamental dataset to underpin further advances in schistosome research.

Malaria infections containing multiple parasite genotypes are ubiquitous in nature, and play a central role in models of recombination, intra-host dynamics, virulence, sex ratio, immunity and drug resistance evolution in Plasmodium. While these multiple infections (MIs) are often assumed to result from superinfection (bites from multiple infected mosquitoes), we know remarkably little about their composition or generation. We isolated 336 parasite clones from eight patients from Malawi (high transmission) and six from Thailand (low transmission) by dilution cloning. These were genotyped using 384 single-nucleotide polymorphisms, revealing 22 independent haplotypes in Malawi (2-6 per MI) and 15 in Thailand (2-5 per MI). Surprisingly, all six patients from Thailand and six of eight from Malawi contained related haplotypes, and haplotypes were more similar within- than between-infections. These results argue against a simple superinfection model. Instead, the observed kinship patterns may be explained by inoculation of multiple related haploid sporozoites from single mosquito bites, by immune suppression of parasite subpopulations within infections, and serial transmission of related parasites between people. That relatedness is maintained in endemic areas in the face of repeated bites from infected mosquitoes has profound implications for understanding malaria transmission, immunity and intra-host dynamics of co-infecting parasite genotypes.

The diversity in the Plasmodium falciparum genome can be used to explore parasite population dynamics, with practical applications to malaria control. The ability to identify the geographic origin and trace the migratory patterns of parasites with clinically important phenotypes such as drug resistance is particularly relevant. With increasing single-nucleotide polymorphism (SNP) discovery from ongoing Plasmodium genome sequencing projects, a demand for high SNP and sample throughput genotyping platforms for large-scale population genetic studies is required. Low parasitaemias and multiple clone infections present a number of challenges to genotyping P. falciparum. We addressed some of these issues using a custom 384-SNP Illumina GoldenGate assay on P. falciparum DNA from laboratory clones (long-term cultured adapted parasite clones), short-term cultured parasite isolates and clinical (non-cultured isolates) samples from East and West Africa, Southeast Asia and Oceania. Eighty percent of the SNPs (n = 306) produced reliable genotype calls on samples containing as little as 2 ng of total genomic DNA and on whole genome amplified DNA. Analysis of artificial mixtures of laboratory clones demonstrated high genotype calling specificity and moderate sensitivity to call minor frequency alleles. Clear resolution of geographically distinct populations was demonstrated using Principal Components Analysis (PCA), and global patterns of population genetic diversity were consistent with previous reports. These results validate the utility of the platform in performing population genetic studies of P. falciparum.