Whole-genome surveillance identifies markers of Plasmodium falciparum drug resistance and novel genomic regions under selection in Mozambique.

IMPORTANCE: Malaria is a devastating disease caused by Plasmodium parasites. The evolution of parasite drug resistance continues to hamper progress toward malaria elimination, and despite extensive efforts to control malaria, it remains a leading cause of death in Mozambique and other countries in the region. The development of successful vaccines and identification of molecular markers to track drug efficacy are essential for managing the disease burden. We present an analysis of the parasite genome in Mozambique, a country with one of the highest malaria burdens globally and limited available genomic data, revealing current selection pressure. We contribute additional evidence to limited prior studies supporting the effectiveness of SWGA in producing reliable genomic data from complex clinical samples. Our results provide the identity of genomic loci that may be associated with current antimalarial drug use, including artemisinin and lumefantrine, and reveal selection pressure predicted to compromise the efficacy of current vaccine candidates.

The Plasmodium falciparum ABC transporter ABCI3 confers parasite strain-dependent pleiotropic antimalarial drug resistance.

Widespread Plasmodium falciparum resistance to first-line antimalarials underscores the vital need to develop compounds with novel modes of action and identify new druggable targets. Here, we profile five compounds that potently inhibit P. falciparum asexual blood stages. Resistance selection studies with three carboxamide-containing compounds, confirmed by gene editing and conditional knockdowns, identify point mutations in the parasite transporter ABCI3 as the primary mediator of resistance. Selection studies with imidazopyridine or quinoline-carboxamide compounds also yield changes in ABCI3, this time through gene amplification. Imidazopyridine mode of action is attributed to inhibition of heme detoxification, as evidenced by cellular accumulation and heme fractionation assays. For the copy-number variation-selecting imidazopyridine and quinoline-carboxamide compounds, we find that resistance, manifesting as a biphasic concentration-response curve, can independently be mediated by mutations in the chloroquine resistance transporter PfCRT. These studies reveal the interconnectedness of P. falciparum transporters in overcoming drug pressure in different parasite strains.

The Key Glycolytic Enzyme Phosphofructokinase Is Involved in Resistance to Antiplasmodial Glycosides.

Plasmodium parasites rely heavily on glycolysis for ATP production and for precursors for essential anabolic pathways, such as the methylerythritol phosphate (MEP) pathway. Here, we show that mutations in the Plasmodium falciparum glycolytic enzyme, phosphofructokinase (PfPFK9), are associated with in vitro resistance to a primary sulfonamide glycoside (PS-3). Flux through the upper glycolysis pathway was significantly reduced in PS-3-resistant parasites, which was associated with reduced ATP levels but increased flux into the pentose phosphate pathway. PS-3 may directly or indirectly target enzymes in these pathways, as PS-3-treated parasites had elevated levels of glycolytic and tricarboxylic acid (TCA) cycle intermediates. PS-3 resistance also led to reduced MEP pathway intermediates, and PS-3-resistant parasites were hypersensitive to the MEP pathway inhibitor, fosmidomycin. Overall, this study suggests that PS-3 disrupts core pathways in central carbon metabolism, which is compensated for by mutations in PfPFK9, highlighting a novel metabolic drug resistance mechanism in P. falciparumIMPORTANCE Malaria, caused by Plasmodium parasites, continues to be a devastating global health issue, causing 405,000 deaths and 228 million cases in 2018. Understanding key metabolic processes in malaria parasites is critical to the development of new drugs to combat this major infectious disease. The Plasmodium glycolytic pathway is essential to the malaria parasite, providing energy for growth and replication and supplying important biomolecules for other essential Plasmodium anabolic pathways. Despite this overreliance on glycolysis, no current drugs target glycolysis, and there is a paucity of information on critical glycolysis targets. Our work addresses this unmet need, providing new mechanistic insights into this key pathway.

Bartonella quintana Endocarditis in a Homeless Man with Cat Exposure in San Diego, California.

We report a case of Bartonella quintana endocarditis in a homeless man with congenital bicuspid aortic valve and significant cat exposure living in downtown San Diego, California.

Advances in omics-based methods to identify novel targets for malaria and other parasitic protozoan infections.

A major advance in antimalarial drug discovery has been the shift towards cell-based phenotypic screening, with notable progress in the screening of compounds against the asexual blood stage, liver stage, and gametocytes. A primary method for drug target deconvolution in Plasmodium falciparum is in vitro evolution of compound-resistant parasites followed by whole-genome scans. Several of the most promising antimalarial drug targets, such as translation elongation factor 2 (eEF2) and phenylalanine tRNA synthetase (PheRS), have been identified or confirmed using this method. One drawback of this method is that if a mutated gene is uncharacterized, a substantial effort may be required to determine whether it is a drug target, a drug resistance gene, or if the mutation is merely a background mutation. Thus, the availability of high-throughput, functional genomic datasets can greatly assist with target deconvolution. Studies mapping genome-wide essentiality in P. falciparum or performing transcriptional profiling of the host and parasite during liver-stage infection with P. berghei have identified potentially druggable pathways. Advances in mapping the epigenomic regulation of the malaria parasite genome have also enabled the identification of key processes involved in parasite development. In addition, the examination of the host genome during infection has identified novel gene candidates associated with susceptibility to severe malaria. Here, we review recent studies that have used omics-based methods to identify novel targets for interventions against protozoan parasites, focusing on malaria, and we highlight the advantages and limitations of the approaches used. These approaches have also been extended to other protozoan pathogens, including Toxoplasma, Trypanosoma, and Leishmania spp., and these studies highlight how drug discovery efforts against these pathogens benefit from the utilization of diverse omics-based methods to identify promising drug targets.

The genomic architecture of antimalarial drug resistance.

Plasmodium falciparum and Plasmodium vivax, the two protozoan parasite species that cause the majority of cases of human malaria, have developed resistance to nearly all known antimalarials. The ability of malaria parasites to develop resistance is primarily due to the high numbers of parasites in the infected person's bloodstream during the asexual blood stage of infection in conjunction with the mutability of their genomes. Identifying the genetic mutations that mediate antimalarial resistance has deepened our understanding of how the parasites evade our treatments and reveals molecular markers that can be used to track the emergence of resistance in clinical samples. In this review, we examine known genetic mutations that lead to resistance to the major classes of antimalarial medications: the 4-aminoquinolines (chloroquine, amodiaquine and piperaquine), antifolate drugs, aryl amino-alcohols (quinine, lumefantrine and mefloquine), artemisinin compounds, antibiotics (clindamycin and doxycycline) and a napthoquinone (atovaquone). We discuss how the evolution of antimalarial resistance informs strategies to design the next generation of antimalarial therapies.

Dual RNA-seq identifies human mucosal immunity protein Mucin-13 as a hallmark of Plasmodium exoerythrocytic infection.

The exoerythrocytic stage of Plasmodium infection is a critical window for prophylactic intervention. Using genome-wide dual RNA sequencing of flow-sorted infected and uninfected hepatoma cells we show that the human mucosal immunity gene, mucin-13 (MUC13), is strongly upregulated during Plasmodium exoerythrocytic hepatic-stage infection. We confirm MUC13 transcript increases in hepatoma cell lines and primary hepatocytes. In immunofluorescence assays, host MUC13 protein expression distinguishes infected cells from adjacent uninfected cells and shows similar colocalization with parasite biomarkers such as UIS4 and HSP70. We further show that localization patterns are species independent, marking both P. berghei and P. vivax infected cells, and that MUC13 can be used to identify compounds that inhibit parasite replication in hepatocytes. This data provides insights into host-parasite interactions in Plasmodium infection, and demonstrates that a component of host mucosal immunity is reprogrammed during the progression of infection.

Exploration of Plasmodium vivax transmission dynamics and recurrent infections in the Peruvian Amazon using whole genome sequencing.

BACKGROUND: Plasmodium vivax poses a significant challenge to malaria elimination due to its ability to cause relapsed infections from reactivation of dormant liver parasites called hypnozoites. We analyzed 69 P. vivax whole genome sequences obtained from subjects residing in three different villages along the Peruvian Amazon. This included 23 paired P. vivax samples from subjects who experienced recurrent P. vivax parasitemia following observed treatment with chloroquine and primaquine. METHODS: Genomic DNA was extracted from whole blood samples collected from subjects. P. vivax DNA was enriched using selective whole genome amplification and whole genome sequencing. We used single nucleotide polymorphisms (SNPs) from the core P. vivax genome to determine characteristics of the parasite population using discriminant analysis of principal components, maximum likelihood estimation of individual ancestries, and phylogenetic analysis. We estimated the relatedness of the paired samples by calculating the number of segregating sites and using a hidden Markov model approach to estimate identity by descent. RESULTS: We present a comprehensive dataset of population genetics of Plasmodium vivax in the Peruvian Amazonian. We define the parasite population structure in this region and demonstrate a novel method for distinguishing homologous relapses from reinfections or heterologous relapses with improved accuracy. The parasite population in this area was quite diverse with an estimated five subpopulations and evidence of a highly heterogeneous ancestry of some of the isolates, similar to previous analyses of P. vivax in this region. Pairwise comparison of recurrent infections determined that there were 12 homologous relapses and 3 likely heterologous relapses with highly related parasites. To the best of our knowledge, this is the first large-scale study to evaluate recurrent P. vivax infections using whole genome sequencing. CONCLUSIONS: Whole genome sequencing is a high-resolution tool that can identify P. vivax homologous relapses with increased sensitivity, while also providing data about drug resistance and parasite population genetics. This information is important for evaluating the efficacy of known and novel antirelapse medications in endemic areas and thus advancing the campaign to eliminate malaria.

Antimalarial activity of single-dose DSM265, a novel plasmodium dihydroorotate dehydrogenase inhibitor, in patients with uncomplicated Plasmodium falciparum or Plasmodium vivax malaria infection: a proof-of-concept, open-label, phase 2a study.

BACKGROUND: DSM265 is a novel, long-duration inhibitor of plasmodium dihydroorotate dehydrogenase (DHODH) with excellent selectivity over human DHODH and activity against blood and liver stages of Plasmodium falciparum. This study aimed to assess the efficacy of DSM265 in patients with P falciparum or Plasmodium vivax malaria infection. METHODS: This proof-of-concept, open-label, phase 2a study was conducted at the Asociación Civil Selva Amazónica in Iquitos, Peru. Patients aged 18-70 years, weighing 45-90 kg, who had clinical malaria (P falciparum or P vivax monoinfection) and fever within the previous 24 h were eligible. Exclusion criteria were clinical or laboratory signs of severe malaria, inability to take oral medicine, and use of other antimalarial treatment in the preceding 14 days. Patients were divided into cohorts of those with P falciparum (cohort a) or P vivax (cohort b) infection. Two initial cohorts received single oral doses of 400 mg DSM265. Patients were followed up for efficacy for 28 days and safety for 35 days. Further cohorts received escalated or de-escalated doses of DSM265, after safety and efficacy assessment of the initial dose. The primary endpoints were the proportion of patients achieving PCR-adjusted adequate clinical and parasitological response (ACPR) by day 14 for patients infected with P falciparum and the proportion of patients achieving a crude cure by day 14 for those infected with P vivax. Cohort success, the criteria for dose escalation, was defined as ACPR (P falciparum) or crude cure (P vivax) in at least 80% of patients in the cohort. The primary analysis was done in the intention-to-treat population (ITT) and the per-protocol population, and safety analyses were done in all patients who received the study drug. This study is registered at ClinicalTrials.gov (NCT02123290). FINDINGS: Between Jan 12, 2015, and Dec 2, 2015, 45 Peruvian patients (24 with P falciparum [cohort a] and 21 with P vivax [cohort b] infection) were sequentially enrolled. For patients with P falciparum malaria in the per-protocol population, all 11 (100%) in the 400 mg group and eight (80%) of ten in the 250 mg group achieved ACPR on day 14. In the ITT analysis, 11 (85%) of 13 in the 400 mg group and eight (73%) of 11 in the 250 mg group achieved ACPR at day 14. For the patients with P vivax malaria, the primary endpoint was not met. In the per-protocol analysis, none of four patients who had 400 mg, three (50%) of six who had 600 mg, and one (25%) of four who had 800 mg DSM265 achieved crude cure at day 14. In the ITT analysis, none of five in the 400 mg group, three (33%) of nine in the 600 mg group, and one (14%) of seven in the 800 mg group achieved crude cure at day 14. During the 28-day extended observation of P falciparum patients, a resistance-associated mutation in the gene encoding the DSM265 target DHODH was observed in two of four recurring patients. DSM265 was well tolerated. The most common adverse events were pyrexia (20 [44%] of 45) and headache (18 [40%] of 45), which are both common symptoms of malaria, and no patients had any treatment-related serious adverse events or adverse events leading to study discontinuation. INTERPRETATION: After a single dose of DSM265, P falciparum parasitaemia was rapidly cleared, whereas against P vivax, DSM265 showed less effective clearance kinetics. Its long duration of action provides the potential to prevent recurrence of P falciparum after treatment with a single dose, which should be further assessed in future combination studies. FUNDING: The Global Health Innovative Technology Fund, the Bill & Melinda Gates Foundation, the National Institutes of Health (R01 AI103058), the Wellcome Trust, and the UK Department of International Development.

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.

A Variant PfCRT Isoform Can Contribute to Plasmodium falciparum Resistance to the First-Line Partner Drug Piperaquine.

Current efforts to reduce the global burden of malaria are threatened by the rapid spread throughout Asia of Plasmodium falciparum resistance to artemisinin-based combination therapies, which includes increasing rates of clinical failure with dihydroartemisinin plus piperaquine (PPQ) in Cambodia. Using zinc finger nuclease-based gene editing, we report that addition of the C101F mutation to the chloroquine (CQ) resistance-conferring PfCRT Dd2 isoform common to Asia can confer PPQ resistance to cultured parasites. Resistance was demonstrated as significantly higher PPQ concentrations causing 90% inhibition of parasite growth (IC90) or 50% parasite killing (50% lethal dose [LD50]). This mutation also reversed Dd2-mediated CQ resistance, sensitized parasites to amodiaquine, quinine, and artemisinin, and conferred amantadine and blasticidin resistance. Using heme fractionation assays, we demonstrate that PPQ causes a buildup of reactive free heme and inhibits the formation of chemically inert hemozoin crystals. Our data evoke inhibition of heme detoxification in the parasite's acidic digestive vacuole as the primary mode of both the bis-aminoquinoline PPQ and the related 4-aminoquinoline CQ. Both drugs also inhibit hemoglobin proteolysis at elevated concentrations, suggesting an additional mode of action. Isogenic lines differing in their pfmdr1 copy number showed equivalent PPQ susceptibilities. We propose that mutations in PfCRT could contribute to a multifactorial basis of PPQ resistance in field isolates.IMPORTANCE The global agenda to eliminate malaria depends on the continued success of artemisinin-based combination therapies (ACTs), which target the asexual blood stages of the intracellular parasite Plasmodium Partial resistance to artemisinin, however, is now established in Southeast Asia, exposing the partner drugs to increased selective pressure. Plasmodium falciparum resistance to the first-line partner piperaquine (PPQ) is now spreading rapidly in Cambodia, resulting in clinical treatment failures. Here, we report that a variant form of the Plasmodium falciparum chloroquine resistance transporter, harboring a C101F mutation edited into the chloroquine (CQ)-resistant Dd2 isoform prevalent in Asia, can confer PPQ resistance in cultured parasites. This was accompanied by a loss of CQ resistance. Biochemical assays showed that PPQ, like CQ, inhibits the detoxification of reactive heme that is formed by parasite-mediated catabolism of host hemoglobin. We propose that novel PfCRT variants emerging in the field could contribute to a multigenic basis of PPQ resistance.

Selective Whole-Genome Amplification Is a Robust Method That Enables Scalable Whole-Genome Sequencing of Plasmodium vivax from Unprocessed Clinical Samples.

Whole-genome sequencing (WGS) of microbial pathogens from clinical samples is a highly sensitive tool used to gain a deeper understanding of the biology, epidemiology, and drug resistance mechanisms of many infections. However, WGS of organisms which exhibit low densities in their hosts is challenging due to high levels of host genomic DNA (gDNA), which leads to very low coverage of the microbial genome. WGS of Plasmodium vivax, the most widely distributed form of malaria, is especially difficult because of low parasite densities and the lack of an ex vivo culture system. Current techniques used to enrich P. vivax DNA from clinical samples require significant resources or are not consistently effective. Here, we demonstrate that selective whole-genome amplification (SWGA) can enrich P. vivax gDNA from unprocessed human blood samples and dried blood spots for high-quality WGS, allowing genetic characterization of isolates that would otherwise have been prohibitively expensive or impossible to sequence. We achieved an average genome coverage of 24×, with up to 95% of the P. vivax core genome covered by ≥5 reads. The single-nucleotide polymorphism (SNP) characteristics and drug resistance mutations seen were consistent with those of other P. vivax sequences from a similar region in Peru, demonstrating that SWGA produces high-quality sequences for downstream analysis. SWGA is a robust tool that will enable efficient, cost-effective WGS of P. vivax isolates from clinical samples that can be applied to other neglected microbial pathogens. IMPORTANCE: Malaria is a disease caused by Plasmodium parasites that caused 214 million symptomatic cases and 438,000 deaths in 2015. Plasmodium vivax is the most widely distributed species, causing the majority of malaria infections outside sub-Saharan Africa. Whole-genome sequencing (WGS) of Plasmodium parasites from clinical samples has revealed important insights into the epidemiology and mechanisms of drug resistance of malaria. However, WGS of P. vivax is challenging due to low parasite levels in humans and the lack of a routine system to culture the parasites. Selective whole-genome amplification (SWGA) preferentially amplifies the genomes of pathogens from mixtures of target and host gDNA. Here, we demonstrate that SWGA is a simple, robust method that can be used to enrich P. vivax genomic DNA (gDNA) from unprocessed human blood samples and dried blood spots for cost-effective, high-quality WGS.