Publications

Molecular surveillance provides a powerful approach to monitoring the resistance status of parasite populations in the field and for understanding resistance evolution. Oxamniquine was used to treat Brazilian schistosomiasis patients (mid-1970s to mid-2000s) and several cases of parasite infections resistant to treatment were recorded. The gene underlying resistance (SmSULT-OR) encodes a sulfotransferase required for intracellular drug activation. Resistance has a recessive basis and occurs when both SmSULT-OR alleles encode for defective proteins. Here we examine SmSULT-OR sequence variation in a natural schistosome population in Brazil ∼40years after the first use of this drug. We sequenced SmSULT-OR from 189 individual miracidia (1-11 per patient) recovered from 49 patients, and tested proteins expressed from putative resistance alleles for their ability to activate oxamniquine. We found nine mutations (four non-synonymous single nucleotide polymorphisms, three non-coding single nucleotide polymorphisms and two indels). Both mutations (p.E142del and p.C35R) identified previously were recovered in this field population. We also found two additional mutations (a splice site variant and 1bp coding insertion) predicted to encode non-functional truncated proteins. Two additional substitutions (p.G206V, p.N215Y) tested had no impact on oxamniquine activation. Three results are of particular interest: (i) we recovered the p.E142del mutation from the field: this same deletion is responsible for resistance in an oxamniquine selected laboratory parasite population; (ii) frequencies of resistance alleles are extremely low (0.27-0.8%), perhaps due to fitness costs associated with carriage of these alleles; (iii) that four independent resistant alleles were found is consistent with the idea that multiple mutations can generate loss-of-function alleles.

BACKGROUND: Deployment of mefloquine-artesunate (MAS3) on the Thailand-Myanmar border has led to a sustained reduction in falciparum malaria, although antimalarial efficacy has declined substantially in recent years. The role of Plasmodium falciparum K13 mutations (a marker of artemisinin resistance) in reducing treatment efficacy remains controversial.

METHODS: Between 2003 and 2013, we studied the efficacy of MAS3 in 1005 patients with uncomplicated P. falciparum malaria in relation to molecular markers of resistance.

RESULTS: Polymerase chain reaction (PCR)-adjusted cure rates declined from 100% in 2003 to 81.1% in 2013 as the proportions of isolates with multiple Pfmdr1 copies doubled from 32.4% to 64.7% and those with K13 mutations increased from 6.7% to 83.4%. K13 mutations conferring moderate artemisinin resistance (notably E252Q) predominated initially but were later overtaken by propeller mutations associated with slower parasite clearance (notably C580Y). Those infected with both multiple Pfmdr1 copy number and a K13 propeller mutation were 14 times more likely to fail treatment. The PCR-adjusted cure rate was 57.8% (95% confidence interval [CI], 45.4, 68.3) compared with 97.8% (95% CI, 93.3, 99.3) in patients with K13 wild type and Pfmdr1 single copy. K13 propeller mutation alone was a strong risk factor for recrudescence (P = .009). The combined population attributable fraction of recrudescence associated with K13 mutation and Pfmdr1 amplification was 82%.

CONCLUSIONS: The increasing prevalence of K13 mutations was the decisive factor for the recent and rapid decline in efficacy of artemisinin-based combination (MAS3) on the Thailand-Myanmar border.

The gender of cercarial larvae can only be determined using molecular methods. End point PCR methods that amplify repetitive markers on the W chromosome of the female (ZW) parasites have been developed, but sometimes results are ambiguous or incorrect. To more effectively distinguish sexes, and to determine why end point PCR can be incorrect, we quantified the W6 repeat sequence and a specific Z chromosome gene using real-time PCR. The ratio between copy number of W6 and a Z chromosome marker unambiguously identifies gender: females have higher ratios (421-4371) than males (0-21). However, some males have low numbers of W6 elements in their genome, and qPCR demonstrated significantly higher W6/Z marker ratios for male genotypes giving ambiguous end point PCR results compared with males giving clear end point results. The quantitative PCR sexing method developed will be particularly useful where reliable sexing of cercariae is critical, for example when staging genetic crosses.

BACKGROUND: Biomphalaria snails are the intermediate host of the blood fluke Schistosoma mansoni, which infect more than 67 million people in tropical areas. Phenoloxidase enzymes (POs), including tyrosinases, catecholases, and laccases, are known to play a role in the immune defenses of arthropods, but the PO activity present in Biomphalaria spp. hemolymph has not been characterized. This study was designed to characterize substrate specificity and reaction optima of PO activity in Biomphalaria spp. hemolymph as a starting point to understand the role of this important invertebrate enzyme activity in snail biology and snail-schistosome interactions.

METHODS: We used spectrophotometric assays with 3 specific substrates (L-tyrosine for tyrosinase, L-DOPA for catecholase, and PPD for laccase) and diethylthiocarbarmate (DETC) as specific PO inhibitor to characterize PO activity in the hemolymph of uninfected snails from two Biomphalaria species, and to determine the impact of the parasite Schistosoma mansoni on the PO activity of its B. glabrata vector.

RESULTS: We identified laccase activity in hemolymph from uninfected B. glabrata and B. alexandrina. For both species, the activity was optimal at 45 °C and pH 8.5, and located in the plasma. The K m and V max of PO enzymes are 1.45 mM and 0.024 OD.min(-1) for B. glabrata, and 1.19 mM and 0.025 OD.min(-1) for B. alexandrina. When the snail vector is parasitized by S. mansoni, we observed a sharp reduction in laccase activity seven weeks after snail infection.

CONCLUSIONS: We employed a highly specific spectrophotometric assay using PPD substrate which allows accurate measurement of laccase activity in Biomphalaria spp. hemolymph. We also demonstrated a strong impact of the parasite S. mansoni on laccase activity in the snail host.

Some of the most valuable antimalarial compounds, including quinine and artemisinin, originated from plants. While these drugs have served important roles over many years for the treatment of malaria, drug resistance has become a widespread problem. Therefore, a critical need exists to identify new compounds that have efficacy against drug-resistant malaria strains. In the current study, extracts prepared from plants readily obtained from local sources were screened for activity against Plasmodium falciparum. Bioassay-guided fractionation was used to identify 18 compounds from five plant species. These compounds included eight lupane triterpenes (1-8), four kaempferol 3-O-rhamnosides (10-13), four kaempferol 3-O-glucosides (14-17), and the known compounds amentoflavone and knipholone. These compounds were tested for their efficacy against multi-drug-resistant malaria parasites and counterscreened against HeLa cells to measure their antimalarial selectivity. Most notably, one of the new lupane triterpenes (3) isolated from the supercritical extract of Buxus sempervirens, the common boxwood, showed activity against both drug-sensitive and -resistant malaria strains at a concentration that was 75-fold more selective for the drug-resistant malaria parasites as compared to HeLa cells. This study demonstrates that new antimalarial compounds with efficacy against drug-resistant strains can be identified from native and introduced plant species in the United States, which traditionally have received scant investigation compared to more heavily explored tropical and semitropical botanical resources from around the world.

BACKGROUND: For over two decades, a racemic mixture of oxamniquine (OXA) was administered to patients infected by Schistosoma mansoni, but whether one or both enantiomers exert antischistosomal activity was unknown. Recently, a  30 kDa S. mansoni sulfotransferase (SmSULT) was identified as the target of OXA action.

METHODOLOGY/PRINCIPAL FINDINGS: Here, we separate the OXA enantiomers using chromatographic methods and assign their optical activities as dextrorotary [(+)-OXA] or levorotary [(-)-OXA]. Crystal structures of the parasite enzyme in complex with optically pure (+)-OXA and (-)-OXA) reveal their absolute configurations as S- and R-, respectively. When tested in vitro, S-OXA demonstrated the bulk of schistosomicidal activity, while R-OXA had antischistosomal effects when present at relatively high concentrations. Crystal structures R-OXA•SmSULT and S-OXA•SmSULT complexes reveal similarities in the modes of OXA binding, but only the S-OXA enantiomer is observed in the structure of the enzyme exposed to racemic OXA.

CONCLUSIONS/SIGNIFICANCE: Together the data suggest the higher schistosomicidal activity of S-OXA is correlated with its ability to outcompete R-OXA binding the sulfotransferase active site. These findings have important implications for the design, syntheses, and dosing of new OXA-based antischistosomal compounds.

BACKGROUND: Emergence of artemisinin resistance in southeast Asia poses a serious threat to the global control of Plasmodium falciparum malaria. Discovery of the K13 marker has transformed approaches to the monitoring of artemisinin resistance, allowing introduction of molecular surveillance in remote areas through analysis of DNA. We aimed to assess the spread of artemisinin-resistant P falciparum in Myanmar by determining the relative prevalence of P falciparum parasites carrying K13-propeller mutations.

METHODS: We did this cross-sectional survey at malaria treatment centres at 55 sites in ten administrative regions in Myanmar, and in relevant border regions in Thailand and Bangladesh, between January, 2013, and September, 2014. K13 sequences from P falciparum infections were obtained mainly by passive case detection. We entered data into two geostatistical models to produce predictive maps of the estimated prevalence of mutations of the K13 propeller region across Myanmar.

FINDINGS: Overall, 371 (39%) of 940 samples carried a K13-propeller mutation. We recorded 26 different mutations, including nine mutations not described previously in southeast Asia. In seven (70%) of the ten administrative regions of Myanmar, the combined K13-mutation prevalence was more than 20%. Geospatial mapping showed that the overall prevalence of K13 mutations exceeded 10% in much of the east and north of the country. In Homalin, Sagaing Region, 25 km from the Indian border, 21 (47%) of 45 parasite samples carried K13-propeller mutations.

INTERPRETATION: Artemisinin resistance extends across much of Myanmar. We recorded P falciparum parasites carrying K13-propeller mutations at high prevalence next to the northwestern border with India. Appropriate therapeutic regimens should be tested urgently and implemented comprehensively if spread of artemisinin resistance to other regions is to be avoided.

FUNDING: Wellcome Trust-Mahidol University-Oxford Tropical Medicine Research Programme and the Bill & Melinda Gates Foundation.

BACKGROUND: Artemisinin-resistant falciparum malaria has emerged in Southeast Asia, posing a major threat to malaria control. It is characterised by delayed asexual-stage parasite clearance, which is the reference comparator for the molecular marker 'Kelch 13' and in vitro sensitivity tests. However, current cut-off values denoting slow clearance based on the proportion of individuals remaining parasitaemic on the third day of treatment ('day-3'), or on peripheral blood parasite half-life, are not well supported. We here explore the parasite clearance distributions in an area of artemisinin resistance with the aim refining the in vivo phenotypic definitions.

METHODS AND FINDINGS: Data from 1,518 patients on the Thai-Myanmar and Thai-Cambodian borders with parasite half-life assessments after artesunate treatment were analysed. Half-lives followed a bimodal distribution. A statistical approach was developed to infer the characteristics of the component distributions and their relative contribution to the composite mixture. A model representing two parasite subpopulations with geometric mean (IQR) parasite half-lives of 3.0 (2.4-3.9) hours and 6.50 (5.7-7.4) hours was consistent with the data. For individual patients, the parasite half-life provided a predicted likelihood of an artemisinin-resistant infection which depends on the population prevalence of resistance in that area. Consequently, a half-life where the probability is 0.5 varied between 3.5 and 5.5 hours. Using this model, the current 'day-3' cut-off value of 10% predicts the potential presence of artemisinin-resistant infections in most but not all scenarios. These findings are relevant to the low-transmission setting of Southeast Asia. Generalisation to a high transmission setting as in regions of Sub-Saharan Africa will need additional evaluation.

CONCLUSIONS: Characterisation of overlapping distributions of parasite half-lives provides quantitative insight into the relationship between parasite clearance and artemisinin resistance, as well as the predictive value of the 10% cut-off in 'day-3' parasitaemia. The findings are important for the interpretation of in vitro sensitivity tests and molecular markers for artemisinin resistance and for contextualising the 'day 3' threshold to account for initial parasitaemia and sample size.

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.