Publications

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.

This article outlines genome-scale approaches that can be used to identify mutations in malaria (Plasmodium) parasites that underlie drug resistance and contribute to treatment failure. These approaches include genetic mapping by linkage or genome-wide association studies, drug selection and characterization of resistant mutants, and the identification of genome regions under strong recent selection. While these genomic approaches can identify candidate resistance loci, genetic manipulation is needed to demonstrate causality. We therefore also describe the growing arsenal of available transfection approaches for direct incrimination of mutations suspected to play a role in resistance. Our intention is both to review past progress and highlight promising approaches for future investigations.

We recently demonstrated that human p38 mitogen-activated protein kinase (MAPK) inhibitors reduced in vitro and in vivo replication of the protozoan parasites Toxoplasma gondii and Encephalitozoon cuniculi. In this study, we assessed the efficacy of five p38 MAPK inhibitors to block the replication of Plasmodium falciparum in human erythrocytes cultured ex vivo and demonstrate that the pyridinylimidazole RWJ67657 and the pyrrolobenzimidazole RWJ68198 reduced P. falciparum replication, yielded trophozoites that were greatly diminished in size at 24h, and that these two agents interfered with stage differentiation. Interestingly, the chloroquine-resistant strain W2 was significantly more sensitive to these drugs than was the chloroquine-sensitive strain HB3. These results suggest that pyridinylimidazoles and pyrrolobenzimidazoles designed to inhibit human p38 MAPK activation can be developed to treat malaria.

Macroparasite infections (e.g., helminths) remain a major human health concern. However, assessing transmission dynamics is problematic because the direct observation of macroparasite dispersal among hosts is not possible. We used a novel landscape genetics approach to examine transmission of the human roundworm Ascaris lumbricoides in a small human population in Jiri, Nepal. Unexpectedly, we found significant genetic structuring of parasites, indicating the presence of multiple transmission foci within a small sampling area ( approximately 14 km(2)). We analyzed several epidemiological variables, and found that transmission is spatially autocorrelated around households and that transmission foci are stable over time despite extensive human movement. These results would not have been obtainable via a traditional epidemiological study based on worm counts alone. Our data refute the assumption that a single host population corresponds to a single parasite transmission unit, an assumption implicit in many classic models of macroparasite transmission. Newer models have shown that the metapopulation-like pattern observed in our data can adversely affect targeted control strategies aimed at community-wide impacts. Furthermore, the observed metapopulation structure and local mating patterns generate an excess of homozygotes that can accelerate the spread of recessive traits such as drug resistance. Our study illustrates how molecular analyses complement traditional epidemiological information in providing a better understanding of parasite transmission. Similar landscape genetic approaches in other macroparasite systems will be warranted if an accurate depiction of the transmission process is to be used to inform effective control strategies.

Comparative genomic hybridization studies have revealed elevated copy number (CN) at the reticulocyte-binding protein 1 gene (PfRh1) in fast growing lab-adapted parasites, while genetic manipulation demonstrates a causal link between cell invasion and PfRh1 CN. We therefore examined PfRh1 copy number variation (CNV) in 202 single clone parasite isolates from four countries to quantify the extent of CNV within natural populations. Surprisingly, we found that no natural parasite infections showed elevated CN. In contrast, 4/28 independent laboratory reference strains show elevated CN. One possibility is that amplification of PfRh1 (or neighboring loci) is selected during laboratory culture. In the case of FCR3 group of parasites, clone trees show that PfRh1 amplification arose in laboratory lines following establishment in culture. These data show that CNV at PfRh1 is rare or non-existent in natural populations, but can arise during laboratory propagation. We conclude that PfRh1 CNV is not an important determinant of gene expression, cell invasion or growth rate in natural parasite populations.

Malaria parasites vary in phenotypic traits of biomedical or biological interest such as growth rate, virulence, sex ratio and drug resistance, and there is considerable interest in identifying the genes that underlie this variation. An important first step is to determine trait heritability (H(2)). We evaluate two approaches to measuring H(2) in natural parasite populations using relatedness inferred from genetic marker data. We collected single-clone Plasmodium falciparum infections from 185 patients from the Thailand-Burma border, monitored parasite clearance following treatment with artemisinin combination therapy (ACT), measured resistance to six antimalarial drugs and genotyped parasites using 335 microsatellites. We found strong relatedness structure. There were 27 groups of two to eight clonally identical (CI) parasites, and 74 per cent of parasites showed significant relatedness to one or more other parasites. Initially, we used matrices of allele sharing and variance components (VC) methods to estimate H(2). Inhibitory concentrations (IC(50)) for six drugs showed significant H(2) (0.24 to 0.79, p = 0.06 to 2.85 x 10(-9)), demonstrating that this study design has adequate power. However, a phenotype of current interest–parasite clearance following ACT–showed no detectable heritability (H(2) = 0-0.09, ns) in this population. The existence of CI parasites allows the use of a simple ANOVA approach for quantifying H(2), analogous to that used in human twin studies. This gave similar results to the VC method and requires considerably less genotyping information. We conclude (i) that H(2) can be effectively measured in malaria parasite populations using minimal genotype data, allowing rational design of genome-wide association studies; and (ii) while drug response (IC(50)) shows significant H(2), parasite clearance following ACT was not heritable in the population studied.