Mapping the malaria parasite druggable genome by using in vitro evolution and chemogenomics.

Chemogenetic characterization through in vitro evolution combined with whole-genome analysis can identify antimalarial drug targets and drug-resistance genes. We performed a genome analysis of 262 Plasmodium falciparum parasites resistant to 37 diverse compounds. We found 159 gene amplifications and 148 nonsynonymous changes in 83 genes associated with drug-resistance acquisition, where gene amplifications contributed to one-third of resistance acquisition events. Beyond confirming previously identified multidrug-resistance mechanisms, we discovered hitherto unrecognized drug target-inhibitor pairs, including thymidylate synthase and a benzoquinazolinone, farnesyltransferase and a pyrimidinedione, and a dipeptidylpeptidase and an arylurea. This exploration of the P. falciparum resistome and druggable genome will likely guide drug discovery and structural biology efforts, while also advancing our understanding of resistance mechanisms available to the malaria parasite.

Hexahydroquinolines are antimalarial candidates with potent blood-stage and transmission-blocking activity.

Antimalarial compounds with dual therapeutic and transmission-blocking activity are desired as high-value partners for combination therapies. Here, we report the identification and characterization of hexahydroquinolines (HHQs) that show low nanomolar potency against both pathogenic and transmissible intra-erythrocytic forms of the malaria parasite Plasmodium falciparum. This activity translates into potent transmission-blocking potential, as shown by in vitro male gamete formation assays and reduced oocyst infection and prevalence in Anopheles mosquitoes. In vivo studies illustrated the ability of lead HHQs to suppress Plasmodium berghei blood-stage parasite proliferation. Resistance selection studies, confirmed by CRISPR-Cas9-based gene editing, identified the digestive vacuole membrane-spanning transporter PfMDR1 (P. falciparum multidrug resistance gene-1) as a determinant of parasite resistance to HHQs. Haemoglobin and haem fractionation assays suggest a mode of action that results in reduced haemozoin levels and might involve inhibition of host haemoglobin uptake into intra-erythrocytic parasites. Furthermore, parasites resistant to HHQs displayed increased susceptibility to several first-line antimalarial drugs, including lumefantrine, confirming that HHQs have a different mode of action to other antimalarials drugs for which PfMDR1 is known to confer resistance. This work evokes therapeutic strategies that combine opposing selective pressures on this parasite transporter as an approach to countering the emergence and transmission of multidrug-resistant P. falciparum malaria.

Esterase mutation is a mechanism of resistance to antimalarial compounds.

Pepstatin is a potent peptidyl inhibitor of various malarial aspartic proteases, and also has parasiticidal activity. Activity of pepstatin against cultured Plasmodium falciparum is highly variable depending on the commercial source. Here we identify a minor contaminant (pepstatin butyl ester) as the active anti-parasitic principle. We synthesize a series of derivatives and characterize an analogue (pepstatin hexyl ester) with low nanomolar activity. By selecting resistant parasite mutants, we find that a parasite esterase, PfPARE (P. falciparum Prodrug Activation and Resistance Esterase) is required for activation of esterified pepstatin. Parasites with esterase mutations are resistant to pepstatin esters and to an open source antimalarial compound, MMV011438. Recombinant PfPARE hydrolyses pepstatin esters and de-esterifies MMV011438. We conclude that (1) pepstatin is a potent but poorly bioavailable antimalarial; (2) PfPARE is a functional esterase that is capable of activating prodrugs; (3) Mutations in PfPARE constitute a mechanism of antimalarial resistance.

Plasmodium falciparum Cyclic Amine Resistance Locus (PfCARL), a Resistance Mechanism for Two Distinct Compound Classes.

MMV007564 is a novel antimalarial benzimidazolyl piperidine chemotype identified in cellular screens. To identify the genetic determinant of MMV007564 resistance, parasites were cultured in the presence of the compound to generate resistant lines. Whole genome sequencing revealed distinct mutations in the gene named Plasmodium falciparum cyclic amine resistance locus (pfcarl), encoding a conserved protein of unknown function. Mutations in pfcarl are strongly associated with resistance to a structurally unrelated class of compounds, the imidazolopiperazines, including KAF156, currently in clinical trials. Our data demonstrate that pfcarl mutations confer resistance to two distinct compound classes, benzimidazolyl piperidines and imidazolopiperazines. However, MMV007564 and the imidazolopiperazines, KAF156 and GNF179, have different timings of action in the asexual blood stage and different potencies against the liver and sexual blood stages. These data suggest that pfcarl is a multidrug-resistance gene rather than a common target for benzimidazolyl piperidines and imidazolopiperazines.

A broad analysis of resistance development in the malaria parasite.

Microbial resistance to chemotherapy has caused countless deaths where malaria is endemic. Chemotherapy may fail either due to pre-existing resistance or evolution of drug-resistant parasites. Here we use a diverse set of antimalarial compounds to investigate the acquisition of drug resistance and the degree of cross-resistance against common resistance alleles. We assess cross-resistance using a set of 15 parasite lines carrying resistance-conferring alleles in pfatp4, cytochrome bc1, pfcarl, pfdhod, pfcrt, pfmdr, pfdhfr, cytoplasmic prolyl t-RNA synthetase or hsp90. Subsequently, we assess whether resistant parasites can be obtained after several rounds of drug selection. Twenty-three of the 48 in vitro selections result in resistant parasites, with time to resistance onset ranging from 15 to 300 days. Our data indicate that pre-existing resistance may not be a major hurdle for novel-target antimalarial candidates, and focusing our attention on fast-killing compounds may result in a slower onset of clinical resistance.

CRISPR-Cas9-modified pfmdr1 protects Plasmodium falciparum asexual blood stages and gametocytes against a class of piperazine-containing compounds but potentiates artemisinin-based combination therapy partner drugs.

Emerging resistance to first-line antimalarial combination therapies threatens malaria treatment and the global elimination campaign. Improved therapeutic strategies are required to protect existing drugs and enhance treatment efficacy. We report that the piperazine-containing compound ACT-451840 exhibits single-digit nanomolar inhibition of the Plasmodium falciparum asexual blood stages and transmissible gametocyte forms. Genome sequence analyses of in vitro-derived ACT-451840-resistant parasites revealed single nucleotide polymorphisms in pfmdr1, which encodes a digestive vacuole membrane-bound ATP-binding cassette transporter known to alter P. falciparum susceptibility to multiple first-line antimalarials. CRISPR-Cas9 based gene editing confirmed that PfMDR1 point mutations mediated ACT-451840 resistance. Resistant parasites demonstrated increased susceptibility to the clinical drugs lumefantrine, mefloquine, quinine and amodiaquine. Stage V gametocytes harboring Cas9-introduced pfmdr1 mutations also acquired ACT-451840 resistance. These findings reveal that PfMDR1 mutations can impart resistance to compounds active against asexual blood stages and mature gametocytes. Exploiting PfMDR1 resistance mechanisms provides new opportunities for developing disease-relieving and transmission-blocking antimalarials.

Whole Genome Shotgun Sequencing Shows Selection on Leptospira Regulatory Proteins During in vitro Culture Attenuation.

Leptospirosis is the most common zoonotic disease worldwide with an estimated 500,000 severe cases reported annually, and case fatality rates of 12-25%, due primarily to acute kidney and lung injuries. Despite its prevalence, the molecular mechanisms underlying leptospirosis pathogenesis remain poorly understood. To identify virulence-related genes in Leptospira interrogans, we delineated cumulative genome changes that occurred during serial in vitro passage of a highly virulent strain of L. interrogans serovar Lai into a nearly avirulent isogenic derivative. Comparison of protein coding and computationally predicted noncoding RNA (ncRNA) genes between these two polyclonal strains identified 15 nonsynonymous single nucleotide variant (nsSNV) alleles that increased in frequency and 19 that decreased, whereas no changes in allelic frequency were observed among the ncRNA genes. Some of the nsSNV alleles were in six genes shown previously to be transcriptionally upregulated during exposure to in vivo-like conditions. Five of these nsSNVs were in evolutionarily conserved positions in genes related to signal transduction and metabolism. Frequency changes of minor nsSNV alleles identified in this study likely contributed to the loss of virulence during serial in vitro culture. The identification of new virulence-associated genes should spur additional experimental inquiry into their potential role in Leptospira pathogenesis.

Next-Generation Sequencing of Plasmodium vivax Patient Samples Shows Evidence of Direct Evolution in Drug-Resistance Genes.

Understanding the mechanisms of drug resistance in Plasmodium vivax, the parasite that causes the most widespread form of human malaria, is complicated by the lack of a suitable long-term cell culture system for this parasite. In contrast to P. falciparum, which can be more readily manipulated in the laboratory, insights about parasite biology need to be inferred from human studies. Here we analyze the genomes of parasites within 10 human P. vivax infections from the Peruvian Amazon. Using next-generation sequencing we show that some P. vivax infections analyzed from the region are likely polyclonal. Despite their polyclonality we observe limited parasite genetic diversity by showing that three or fewer haplotypes comprise 94% of the examined genomes, suggesting the recent introduction of parasites into this geographic region. In contrast we find more than three haplotypes in putative drug-resistance genes, including the gene encoding dihydrofolate reductase-thymidylate synthase and the P. vivax multidrug resistance associated transporter, suggesting that resistance mutations have arisen independently. Additionally, several drug-resistance genes are located in genomic regions with evidence of increased copy number. Our data suggest that whole genome sequencing of malaria parasites from patients may provide more insight about the evolution of drug resistance than genetic linkage or association studies, especially in geographical regions with limited parasite genetic diversity.

Identification of pathogen genomic variants through an integrated pipeline.

BACKGROUND: Whole-genome sequencing represents a powerful experimental tool for pathogen research. We present methods for the analysis of small eukaryotic genomes, including a streamlined system (called Platypus) for finding single nucleotide and copy number variants as well as recombination events. RESULTS: We have validated our pipeline using four sets of Plasmodium falciparum drug resistant data containing 26 clones from 3D7 and Dd2 background strains, identifying an average of 11 single nucleotide variants per clone. We also identify 8 copy number variants with contributions to resistance, and report for the first time that all analyzed amplification events are in tandem. CONCLUSIONS: The Platypus pipeline provides malaria researchers with a powerful tool to analyze short read sequencing data. It provides an accurate way to detect SNVs using known software packages, and a novel methodology for detection of CNVs, though it does not currently support detection of small indels. We have validated that the pipeline detects known SNVs in a variety of samples while filtering out spurious data. We bundle the methods into a freely available package.