Whole-genome analysis of Malawian Plasmodium falciparum isolates identifies possible targets of allele-specific immunity to clinical malaria

Individuals acquire immunity to clinical malaria after repeated Plasmodium falciparum infections. Immunity to disease is thought to reflect the acquisition of a repertoire of responses to multiple alleles in diverse parasite antigens. In previous studies, we identified polymorphic sites within individual antigens that are associated with parasite immune evasion by examining antigen allele dynamics in individuals followed longitudinally. Here we expand this approach by analyzing genome-wide polymorphisms using whole genome sequence data from 140 parasite isolates representing malaria cases from a longitudinal study in Malawi and identify 25 genes that encode possible targets of naturally acquired immunity that should be validated immunologically and further characterized for their potential as vaccine candidates.     

Population genetic analysis of Plasmodium falciparum parasites using a customized Illumina GoldenGate genotyping assay

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.     

Genomics meets transgenics in search of the elusive Cryptosporidium drug target

Cryptosporidium is an important pathogen of humans, and a challenging model for the laboratory. The parasite genome sequence, accessible through a comprehensive database, now provides exciting opportunities for urgently needed advances. Comparative genomics, combined with the genetic system in the related parasite Toxoplasma gondii, outlines a detailed Cryptosporidium parvum metabolic map and facilitates cell biological analyses. New targets for Cryptosporidium drug and vaccine development can be identified and validated based on this approach.     

Prediction of the P. falciparum Target Space Relevant to Malaria Drug Discovery

Malaria is still one of the most devastating infectious diseases, affecting hundreds of millions of patients worldwide. Even though there are several established drugs in clinical use for malaria treatment, there is an urgent need for new drugs acting through novel mechanisms of action due to the rapid development of resistance. Resistance emerges when the parasite manages to mutate the sequence of the drug targets to the extent that the protein can still perform its function in the parasite but can no longer be inhibited by the drug, which then becomes almost ineffective. The design of a new generation of malaria drugs targeting multiple essential proteins would make it more difficult for the parasite to develop full resistance without lethally disrupting some of its vital functions. The challenge is then to identify which set of Plasmodium falciparum proteins, among the millions of possible combinations, can be targeted at the same time by a given chemotype. To do that, we predicted first the targets of the close to 20,000 antimalarial hits identified recently in three independent phenotypic screening campaigns. All targets predicted were then projected onto the genome of P. falciparum using orthologous relationships. A total of 226 P. falciparum proteins were predicted to be hit by at least one compound, of which 39 were found to be significantly enriched by the presence and degree of affinity of phenotypically active compounds. The analysis of the chemically compatible target combinations containing at least one of those 39 targets led to the identification of a priority set of 64 multi-target profiles that can set the ground for a new generation of more robust malaria drugs.     

Cas9HF1 enhanced specificity in Ustilago maydis

The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 system is widely used as a tool to precisely manipulate genomic sequence targeted by sgRNA (single guide RNA) and is adapted in different species for genome editing. One of the major concerns of CRISPR-Cas9 is the possibility of off-target effects, which can be remedied by the deployment of high fidelity Cas9 variants. Ustilago maydis is a maize fungal pathogen, which has served as a model organism for biotrophic pathogens for decades. The successful adaption of CRISPR-Cas9 in U. maydis greatly facilitated effector biology studies. Here, we constructed an U. maydis reporter strain that allows in vivo quantification of efficiency and target specificity of three high fidelity Cas9 variants, Cas9HF1, Cas9esp1.1 and Cas9hypa. This approach identified Cas9HF1 as most specific Cas9 variant in U. maydis. Furthermore, whole genome sequencing showed absence of off-target effects in U. maydis by CRISPR-Cas9 editing.     

A Systematic In Silico Search for Target Similarity Identifies Several Approved Drugs with Potential Activity against the Plasmodium falciparum Apicoplast

Most of the drugs in use against Plasmodium falciparum share similar modes of action and, consequently, there is a need to identify alternative potential drug targets. Here, we focus on the apicoplast, a malarial plastid-like organelle of algal source which evolved through secondary endosymbiosis. We undertake a systematic in silico target-based identification approach for detecting drugs already approved for clinical use in humans that may be able to interfere with the P. falciparum apicoplast. The P. falciparum genome database GeneDB was used to compile a list of ≈600 proteins containing apicoplast signal peptides. Each of these proteins was treated as a potential drug target and its predicted sequence was used to interrogate three different freely available databases (Therapeutic Target Database, DrugBank and STITCH3.1) that provide synoptic data on drugs and their primary or putative drug targets. We were able to identify several drugs that are expected to interact with forty-seven (47) peptides predicted to be involved in the biology of the P. falciparum apicoplast. Fifteen (15) of these putative targets are predicted to have affinity to drugs that are already approved for clinical use but have never been evaluated against malaria parasites. We suggest that some of these drugs should be experimentally tested and/or serve as leads for engineering new antimalarials.     

Reduced glycerol incorporation into phospholipids contributes to impaired intra-erythrocytic growth of glycerol kinase knockout Plasmodium falciparum parasites

Background

Malaria is a devastating disease and Plasmodium falciparum is the most lethal parasite infecting humans. Understanding the biology of this parasite is vital in identifying potential novel drug targets. During every 48-hour intra-erythrocytic asexual replication cycle, a single parasite can produce up to 32 progeny. This extensive proliferation implies that parasites require substantial amounts of lipid precursors for membrane biogenesis. Glycerol kinase is a highly conserved enzyme that functions at the interface of lipid synthesis and carbohydrate metabolism. P. falciparum glycerol kinase catalyzes the ATP-dependent phosphorylation of glycerol to glycerol-3-phosphate, a major phospholipid precursor.

Methods

The P. falciparum glycerol kinase gene was disrupted using double crossover homologous DNA recombination to generate a knockout parasite line. Southern hybridization and mRNA analysis were used to verify gene disruption. Parasite growth rates were monitored by flow cytometry. Radiolabelling studies were used to assess incorporation of glycerol into parasite phospholipids.

Results

Disruption of the P. falciparum glycerol kinase gene produced viable parasites, but their growth was significantly reduced to 56.5 ± 1.8% when compared to wild type parasites. ¹⁴C-glycerol incorporation into the major phospholipids of the parasite membrane, phosphatidylcholine and phosphatidylethanolamine, was 48.4 ± 10.8% and 53.1 ± 5.7% relative to an equivalent number of wild type parasites.

Conclusions

P. falciparum glycerol kinase is required for optimal intra-erythrocytic asexual parasite development. Exogenous glycerol may be used as an alternative carbon source for P. falciparum phospholipid biogenesis, despite the lack of glycerol kinase to generate glycerol-3-phosphate.

General significance

These studies provide new insight into glycerolipid metabolism in P. falciparum.     

Multigene Editing in the Escherichia coli Genome via the CRISPR-Cas9 System

An efficient genome-scale editing tool is required for construction of industrially useful microbes. We describe a targeted, continual multigene editing strategy that was applied to the Escherichia coli genome by using the Streptococcus pyogenes type II CRISPR-Cas9 system to realize a variety of precise genome modifications, including gene deletion and insertion, with a highest efficiency of 100%, which was able to achieve simultaneous multigene editing of up to three targets. The system also demonstrated successful targeted chromosomal deletions in Tatumella citrea, another species of the Enterobacteriaceae, with highest efficiency of 100%.     

Construction and use of Plasmodium falciparum phage display libraries to identify host parasite interactions

The development of an electroporation based transfection method for Plasmodium falciparum has been very successful for the study of some genes but its efficiency remains very low. While alternative approaches have been documented, electroporation of infected red blood cells generally remains the method of choice for introducing DNA into P. falciparum. In this paper we compare four published transfection techniques in their ability to achieve stable transfections.     

An Apicoplast Localized Ubiquitylation System Is Required for the Import of Nuclear-encoded Plastid Proteins

Apicomplexan parasites are responsible for numerous important human diseases including toxoplasmosis, cryptosporidiosis, and most importantly malaria. There is a constant need for new antimalarials, and one of most keenly pursued drug targets is an ancient algal endosymbiont, the apicoplast. The apicoplast is essential for parasite survival, and several aspects of its metabolism and maintenance have been validated as targets of anti-parasitic drug treatment. Most apicoplast proteins are nuclear encoded and have to be imported into the organelle. Recently, a protein translocon typically required for endoplasmic reticulum associated protein degradation (ERAD) has been proposed to act in apicoplast protein import. Here, we show ubiquitylation to be a conserved and essential component of this process. We identify apicoplast localized ubiquitin activating, conjugating and ligating enzymes in Toxoplasma gondii and Plasmodium falciparum and observe biochemical activity by in vitro reconstitution. Using conditional gene ablation and complementation analysis we link this activity to apicoplast protein import and parasite survival. Our studies suggest ubiquitylation to be a mechanistic requirement of apicoplast protein import independent to the proteasomal degradation pathway.     

Enhancing Blockade of Plasmodium falciparum Erythrocyte Invasion: Assessing Combinations of Antibodies against PfRH5 and Other Merozoite Antigens

No vaccine has yet proven effective against the blood-stages of Plasmodium falciparum, which cause the symptoms and severe manifestations of malaria. We recently found that PfRH5, a P. falciparum-specific protein expressed in merozoites, is efficiently targeted by broadly-neutralizing, vaccine-induced antibodies. Here we show that antibodies against PfRH5 efficiently inhibit the in vitro growth of short-term-adapted parasite isolates from Cambodia, and that the EC₅₀ values of antigen-specific antibodies against PfRH5 are lower than those against PfAMA1. Since antibody responses elicited by multiple antigens are speculated to improve the efficacy of blood-stage vaccines, we conducted detailed assessments of parasite growth inhibition by antibodies against PfRH5 in combination with antibodies against seven other merozoite antigens. We found that antibodies against PfRH5 act synergistically with antibodies against certain other merozoite antigens, most notably with antibodies against other erythrocyte-binding antigens such as PfRH4, to inhibit the growth of a homologous P. falciparum clone. A combination of antibodies against PfRH4 and basigin, the erythrocyte receptor for PfRH5, also potently inhibited parasite growth. This methodology provides the first quantitative evidence that polyclonal vaccine-induced antibodies can act synergistically against P. falciparum antigens and should help to guide the rational development of future multi-antigen vaccines.     

Global diversity and balancing selection of 23 leading Plasmodium falciparum candidate vaccine antigens

Investigation of the diversity of malaria parasite antigens can help prioritize and validate them as vaccine candidates and identify the most common variants for inclusion in vaccine formulations. Studies of vaccine candidates of the most virulent human malaria parasite, Plasmodium falciparum, have focused on a handful of well-known antigens, while several others have never been studied. Here we examine the global diversity and population structure of leading vaccine candidate antigens of P. falciparum using the MalariaGEN Pf3K (version 5.1) resource, comprising more than 2600 genomes from 15 malaria endemic countries. A stringent variant calling pipeline was used to extract high quality antigen gene ‘haplotypes’ from the global dataset and a new R-package named VaxPack was used to streamline population genetic analyses. In addition, a newly developed algorithm that enables spatial averaging of selection pressure on 3D protein structures was applied to the dataset. We analysed the genes encoding 23 leading and novel candidate malaria vaccine antigens including csp, trap, eba175, ama1, rh5, and CelTOS. Our analysis shows that current malaria vaccine formulations are based on rare haplotypes and thus may have limited efficacy against natural parasite populations. High levels of diversity with evidence of balancing selection was detected for most of the erythrocytic and pre-erythrocytic antigens. Measures of natural selection were then mapped to 3D protein structures to predict targets of functional antibodies. For some antigens, geographical variation in the intensity and distribution of these signals on the 3D structure suggests adaptation to different human host or mosquito vector populations. This study provides an essential framework for the diversity of P. falciparum antigens to be considered in the design of the next generation of malaria vaccines.     

The multifunctional autophagy pathway in the human malaria parasite,Plasmodium falciparum

Autophagy is a catabolic pathway typically induced by nutrient starvation to recycle amino acids, but can also function in removing damaged organelles. In addition, this pathway plays a key role in eukaryotic development. To date, not much is known about the role of autophagy in apicomplexan parasites and more specifically in the human malaria parasite Plasmodium falciparum. Comparative genomic analysis has uncovered some, but not all, orthologs of autophagy-related (ATG) genes in the malaria parasite genome. Here, using a genome-wide in silico analysis, we confirmed that ATG genes whose products are required for vesicle expansion and completion are present, while genes involved in induction of autophagy and cargo packaging are mostly absent. We subsequently focused on the molecular and cellular function of P. falciparum ATG8 (PfATG8), an autophagosome membrane marker and key component of the autophagy pathway, throughout the parasite asexual and sexual erythrocytic stages. In this context, we showed that PfATG8 has a distinct and atypical role in parasite development. PfATG8 localized in the apicoplast and in vesicles throughout the cytosol during parasite development. Immunofluorescence assays of PfATG8 in apicoplast-minus parasites suggest that PfATG8 is involved in apicoplast biogenesis. Furthermore, treatment of parasite cultures with bafilomycin A₁ and chloroquine, both lysosomotropic agents that inhibit autophagosome and lysosome fusion, resulted in dramatic morphological changes of the apicoplast, and parasite death. Furthermore, deep proteomic analysis of components associated with PfATG8 indicated that it may possibly be involved in ribophagy and piecemeal microautophagy of the nucleus. Collectively, our data revealed the importance and specificity of the autophagy pathway in the malaria parasite and offer potential novel therapeutic strategies.     

Complete telomere-to-telomere de novo assembly of the Plasmodium falciparum genome through long-read (>11?kb), single molecule,real-time sequencing

The application of next-generation sequencing to estimate genetic diversity of Plasmodium falciparum, the most lethal malaria parasite, has proved challenging due to the skewed AT-richness [∼80.6% (A + T)] of its genome and the lack of technology to assemble highly polymorphic subtelomeric regions that contain clonally variant, multigene virulence families (Ex: var and rifin). To address this, we performed amplification-free, single molecule, real-time sequencing of P. falciparum genomic DNA and generated reads of average length 12 kb, with 50% of the reads between 15.5 and 50 kb in length. Next, using the Hierarchical Genome Assembly Process, we assembled the P. falciparum genome de novo and successfully compiled all 14 nuclear chromosomes telomere-to-telomere. We also accurately resolved centromeres [∼90–99% (A + T)] and subtelomeric regions and identified large insertions and duplications that add extra var and rifin genes to the genome, along with smaller structural variants such as homopolymer tract expansions. Overall, we show that amplification-free, long-read sequencing combined with de novo assembly overcomes major challenges inherent to studying the P. falciparum genome. Indeed, this technology may not only identify the polymorphic and repetitive subtelomeric sequences of parasite populations from endemic areas but may also evaluate structural variation linked to virulence, drug resistance and disease transmission.     

Plasmodium falciparum UvrD Helicase Translocates in 3′ to 5′ Direction,Colocalizes with MLH and Modulates Its Activity through Physical Interaction

Malaria is a global disease and a major health problem. The control of malaria is a daunting task due to the increasing drug resistance. Therefore, there is an urgent need to identify and characterize novel parasite specific drug targets. In the present study we report the biochemical characterization of parasite specific UvrD helicase from Plasmodium falciparum. The N-terminal fragment (PfUDN) containing UvrD helicase domain, which consists of helicase motifs Q, Ia–Id, II, III and most of motif IV, and the C-terminal fragment (PfUDC1) containing UvrD helicase C terminal domain, consisting of remaining part of motif IV and motifs IVa–IVc and 161 amino acids of intervening sequence between motif IV and V, possess ssDNA-dependent ATPase and DNA helicase activities in vitro. Using immunodepletion assays we show that the ATPase and helicase activities are attributable to PfUDN and PfUDC1 proteins. The helicase activity can utilize the hydrolysis of all the nucleotide and deoxynucleotide triphosphates and the direction of unwinding is 3′ to 5′. The endogenous P. falciparum UvrD contains the characteristic DNA helicase activity. PfUDN interacts with PfMLH (P. falciparum MutL homologue) and modulates the endonuclease activity of PfMLH and PfMLH positively regulates the unwinding activity of PfUDN. We show that PfUvrD is expressed in the nucleus distinctly in the schizont stages of the intraerythrocytic development of the parasite and it colocalizes with PfMLH. These studies will make an important contribution in understanding the nucleic acid transaction in the malaria parasite.