PlasmoDB: the Plasmodium genome resource. A database integrating experimental and computational data

The spatiotemporal transcriptional profiling of murine brain during cerebral malaria progression and after artemisinin treatment

Cerebral malaria (CM) is a severe encephalopathy caused by Plasmodium parasite infection, resulting in thousands of annual deaths and neuro-cognitive sequelae even after anti-malarial drugs treatment. Despite efforts to dissect the mechanism, the cellular transcriptomic reprogramming within the spatial context remains elusive. Here, we constructed single-cell and spatial transcriptome atlases of experimental CM (ECM) male murine brain tissues with or without artesunate (ART) treatment. We identified activated inflammatory endothelial cells during ECM, characterized by a disrupted blood-brain barrier, increased antigen presentation, and leukocyte adhesion. We also observed that inflammatory microglia enhance antigen presentation pathway such as MHC-I to CD8+ cytotoxic T cells. The latter underwent an inflammatory state transition with up-regulated cytokine expression and cytotoxic activity. Multi-omics analysis revealed that the activated interferon-gamma response of injured neurons during ECM and persisted after ART treatment. Overall, our research provides valuable resources for understanding malaria parasite-host interaction mechanisms and adjuvant therapy development.

Profiling antimalarial drug-resistant haplotypes in Pfcrt, Pfmdr1, Pfdhps and Pfdhfr genes in Plasmodium falciparum causing malaria in the Central Region of Ghana: a multicentre cross-sectional study

Background

The proliferation of Plasmodium parasites resistant to antimalarial drugs poses a serious threat to human life and remains an obstacle to managing and eradicating Plasmodium falciparum. The surveillance of molecular markers has become necessary to monitor the spread of resistant haplotypes and discover emerging mutations.

Objective

This molecular epidemiological study aimed to evaluate the prevalence of known mutations in the drug resistance genes Pfcrt, Pfmdr1, Pfdhfr and Pfdhps in the Central Region of Ghana.

Design

A multi-centre cross-sectional study.

Methods

This prospective study utilised dried blood spots from individuals with P. falciparum-infection from five districts in the Central Region of Ghana. Selective Whole Genome Amplification (sWGA) and Single Nucleotide Polymorphisms (SNPs) in P. falciparum chloroquine transporter genes (Pfcrt), P. falciparum multidrug resistance 1 (Pfmdr1), P. falciparum dihydropteroate synthase (Pfdhps) and P. falciparum dihydrofolate reductase (Pfdhfr) were analysed.

Results

Whole genome sequencing was carried out on 522 samples. Of these, 409 (78%) samples were successfully sequenced. Six (6) of the sequenced samples were of co-infection of other parasite species with P. falciparum and excluded from the analysis. Analysis of the Pfcrt gene revealed 0.5% were CVIET (C72, V73, M74I, N75E, K76T) while the Pfcrt CVMNK (C72, V73, M74, N75, K76) wild-type haplotypes were 97% with (2.5%) (CV[M/I][N/E][K/T]) being mixed haplotypes. In the Pfmdr1 gene, monoclonal haplotypes; NFD (N86, Y184F, D1246) and YFN (N86Y, Y184F, D1246N) occurred at 44% and 9.8%, respectively, whereas mixed- haplotypes (N[Y/F]D and [N/Y][Y/F]D) were 23.5% and 0.3%, respectively. Combined Pfdhfr/Pfdhps genes yielded about 88% Pfdhfr IRNI (N51I, C59R, S108N, I164) + Pfdhps A437G haplotypes (conferring partial resistance to Sulphadoxine-Pyrimethamine (SP)) while 9% of the parasites had Pfhdfr IRNI + Pfdhps A437G + K540E haplotypes (conferring full resistance to SP). The wild-type haplotype, Pfdhfr (N51, C59, S108, I164) and Pfdhps (S436, A437, K540, A581, A613) was not observed.

Conclusion

The findings show a low prevalence of CVIET and relatively higher rates for Pfmdr1 NFD and parasites with Pfdhfr IRNI (N51I, C59R, S108N, I164) + Pfdhps A437G haplotypes. These observations advocate for enhanced surveillance which is inimical to malaria management in an endemic area.

Proteomes of plasmodium knowlesi early and late ring-stage parasites and infected host erythrocytes

The emerging malaria parasite Plasmodium knowlesi threatens the goal of worldwide malaria elimination due to its zoonotic spread in Southeast Asia. After brief ex-vivo culture we used 2D LC/MS/MS to examine the early and late ring stages of infected Macaca mulatta red blood cells harboring P. knowlesi. The M. mulatta clathrin heavy chain and T-cell and macrophage inhibitor ERMAP were overexpressed in the early ring stage; glutaredoxin 3 was overexpressed in the late ring stage; GO term differential enrichments included response to oxidative stress and the cortical cytoskeleton in the early ring stage. P. knowlesi clathrin heavy chain and 60S acidic ribosomal protein P2 were overexpressed in the late ring stage; GO term differential enrichments included vacuoles in the early ring stage, ribosomes and translation in the late ring stage, and Golgi- and COPI-coated vesicles, proteasomes, nucleosomes, vacuoles, ion-, peptide-, protein-, nucleocytoplasmic- and RNA-transport, antioxidant activity and glycolysis in both stages. SIGNIFICANCE: Due to its zoonotic spread, cases of the emerging human pathogen Plasmodium knowlesi in southeast Asia, and particularly in Malaysia, threaten regional and worldwide goals for malaria elimination. Infection by this parasite can be fatal to humans, and can be associated with significant morbidity. Due to zoonotic transmission from large macaque reservoirs that are untreatable by drugs, and outdoor biting mosquito vectors that negate use of preventive measures such as bed nets, its containment remains a challenge. Its biology remains incompletely understood. Thus we examine the expressed proteome of the early and late ex-vivo cultured ring stages, the first intraerythrocyte developmental stages after infection of host rhesus macaque erythrocytes. We used GO term enrichment strategies and differential protein expression to compare early and late ring stages. The early ring stage is characterized by the enrichment of P. knowlesi vacuoles, and overexpression of the M. mulatta clathrin heavy chain, important for clathrin-coated pits and vesicles, and clathrin-mediated endocytosis. The M. mulatta protein ERMAP was also overexpressed in the early ring stage, suggesting a potential role in early ring stage inhibition of T-cells and macrophages responding to P. knowlesi infection of reticulocytes. This could allow expansion of the host P. knowlesi cellular niche, allowing parasite adaptation to invasion of a wider age range of RBCs than the preferred young RBCs or reticulocytes, resulting in proliferation and increased pathogenesis in infected humans. Other GO terms differentially enriched in the early ring stage include the M. mulatta cortical cytoskeleton and response to oxidative stress. The late ring stage is characterized by overexpression of the P. knowlesi clathrin heavy chain. Combined with late ring stage GO term enrichment of Golgi-associated and coated vesicles, and enrichment of COPI-coated vesicles in both stages, this suggests the importance to P. knowlesi biology of clathrin-mediated endocytosis. P. knowlesi ribosomes and translation were also differentially enriched in the late ring stage. With expression of a variety of heat shock proteins, these results suggest production of folded parasite proteins is increasing by the late ring stage. M. mulatta endocytosis was differentially enriched in the late ring stage, as were clathrin-coated vesicles and endocytic vesicles. This suggests that M. mulatta clathrin-based endocytosis, perhaps in infected reticulocytes rather than mature RBC, may be an important process in the late ring stage. Additional ring stage biology from enriched GO terms includes M. mulatta proteasomes, protein folding and the chaperonin-containing T complex, actin and cortical actin cytoskeletons. P knowlesi biology also includes proteasomes, as well as nucleosomes, antioxidant activity, a variety of transport processes, glycolysis, vacuoles and protein folding. Mature RBCs have lost internal organelles, suggesting infection here may involve immature reticulocytes still retaining organelles. P. knowlesi parasite proteasomes and translational machinery may be ring stage drug targets for known selective inhibitors of these processes in other Plasmodium species. To our knowledge this is the first examination of more than one timepoint within the ring stage. Our results expand knowledge of both host and parasite proteins, pathways and organelles underlying P. knowlesi ring stage biology.

Infected erythrocytes and plasma proteomics reveal a specific protein signature of severe malaria

Cerebral malaria (CM), the most lethal complication of Plasmodium falciparum severe malaria (SM), remains fatal for 15-25% of affected children despite the availability of treatment. P. falciparum infects and multiplies in erythrocytes, contributing to anemia, parasite sequestration, and inflammation. An unbiased proteomic assessment of infected erythrocytes and plasma samples from 24 Beninese children was performed to study the complex mechanisms underlying CM. A significant down-regulation of proteins from the ubiquitin-proteasome pathway and an up-regulation of the erythroid precursor marker transferrin receptor protein 1 (TFRC) were associated with infected erythrocytes from CM patients. At the plasma level, the samples clustered according to clinical presentation. Significantly, increased levels of the 20S proteasome components were associated with SM. Targeted quantification assays confirmed these findings on a larger cohort (n = 340). These findings suggest that parasites causing CM preferentially infect reticulocytes or erythroblasts and alter their maturation. Importantly, the host plasma proteome serves as a specific signature of SM and presents a remarkable opportunity for developing innovative diagnostic and prognostic biomarkers.

Long-Read Sequencing and De Novo Genome Assembly Pipeline of Two Plasmodium falciparum Clones (Pf3D7, PfW2) Using Only the PromethION Sequencer from Oxford Nanopore Technologies without Whole-Genome Amplification

Antimalarial drug resistance has become a real public health problem despite WHO measures. New sequencing technologies make it possible to investigate genomic variations associated with resistant phenotypes at the genome-wide scale. Based on the use of hemisynthetic nanopores, the PromethION technology from Oxford Nanopore Technologies can produce long-read sequences, in contrast to previous short-read technologies used as the gold standard to sequence Plasmodium. Two clones of P. falciparum (Pf3D7 and PfW2) were sequenced in long-read using the PromethION sequencer from Oxford Nanopore Technologies without genomic amplification. This made it possible to create a processing analysis pipeline for human Plasmodium with ONT Fastq only. De novo assembly revealed N50 lengths of 18,488 kb and 17,502 kb for the Pf3D7 and PfW2, respectively. The genome size was estimated at 23,235,407 base pairs for the Pf3D7 clone and 21,712,038 base pairs for the PfW2 clone. The average genome coverage depth was estimated at 787X and 653X for the Pf3D7 and PfW2 clones, respectively. This study proposes an assembly processing pipeline for the human Plasmodium genome using software adapted to large ONT data and the high AT percentage of Plasmodium. This search provides all the parameters which were optimized for use with the software selected in the pipeline.

The RNA-binding KH-domain in the unique transcription factor of the malaria parasite is responsible for its transcriptional regulatory activity

To date, only a small number of regulatory transcription factors have been predicted from the genome of Plasmodium and Apicomplexan parasites. We previously identified a nuclear factor named Prx regulatory element-binding protein (PREBP) from Plasmodium falciparum. PREBP had been suggested to bind to the cis-element in the promoter of an antioxidant pf1-cys-prx gene, thereby promoting the expression of downstream genes. PREBP has 4 putative K homology (KH) domains, which are known to bind RNA and single-stranded DNA. In this study, to understand the detailed action of PREBP in parasite cells, we first observed that in living parasite cells, PREBP was localized in the nucleus in the trophozoite and schizont stages, in which the expression of the target pf1-cys-prx was enhanced. The interaction of PREBP and the cis-element of pf1-cys-prx in the parasite cells was also confirmed. Further, the activities of PREBP deletion mutants were analyzed, and regions with repeated KH domains in PREBP seemed to be responsible for the recognition of the cis-element. These results led us to hypothesize that Plasmodium and other Apicomplexan parasites might have a transcription factor family with KH domains. Bioinformatic analysis revealed a putative ortholog group including PREBP and several Plasmodium and Apicomplexan factors with KH domains. One of the P. falciparum-derived factors, which were included in the putative ortholog group, was found to be localized at the nucleus in the trophozoite stage, indicating that it might be a novel transcription factor. The discovery of PREBP and putative transcription factors with KH domains suggested that multi-functional proteins with KH domains possibly evolved in the Apicomplexan organisms. They might play key roles in transcriptional regulatory processes that are essential for living organisms and may even represent unique drug targets for malaria therapy.

Protein kinase PfPK2 mediated signalling is critical for host erythrocyte invasion by malaria parasite

Signalling pathways in malaria parasite remain poorly defined and major reason for this is the lack of understanding of the function of majority of parasite protein kinases and phosphatases in parasite signalling and its biology. In the present study, we have elucidated the function of Protein Kinase 2 (PfPK2), which is known to be indispensable for the survival of human malaria parasite Plasmodium falciparum. We demonstrate that it is involved in the invasion of host erythrocytes, which is critical for establishing infection. In addition, PfPK2 may also be involved in the maturation of the parasite post-invasion. PfPK2 regulates the release of microneme proteins like Apical Membrane Antigen 1 (AMA1), which facilitates the formation of Tight Junction between the merozoite and host erythrocyte- a key step in the process of invasion. Comparative phosphoproteomics studies revealed that PfPK2 may be involved in regulation of several key proteins involved in invasion and signalling. Furthermore, PfPK2 regulates the generation of cGMP and the release of calcium in the parasite, which are key second messengers for the process of invasion. These and other studies have shed light on a novel signalling pathway in which PfPK2 acts as an upstream regulator of important cGMP-calcium signalling, which plays an important role in parasite invasion.

The structure of a Plasmodium vivax Tryptophan Rich Antigen domain suggests a lipid binding function for a pan-Plasmodium multi-gene family

Tryptophan Rich Antigens (TRAgs) are encoded by a multi-gene family found in all Plasmodium species, but are significantly expanded in P. vivax and closely related parasites. We show that multiple P. vivax TRAgs are expressed on the merozoite surface and that one, PVP01_0000100 binds red blood cells with a strong preference for reticulocytes. Using X-ray crystallography, we solved the structure of the PVP01_0000100 C-terminal tryptophan rich domain, which defines the TRAg family, revealing a three-helical bundle that is conserved across Plasmodium and has structural homology with lipid-binding BAR domains involved in membrane remodelling. Biochemical assays confirm that the PVP01_0000100 C-terminal domain has lipid binding activity with preference for sulfatide, a glycosphingolipid present in the outer leaflet of plasma membranes. Deletion of the putative orthologue in P. knowlesi, PKNH_1300500, impacts invasion in reticulocytes, suggesting a role during this essential process. Together, this work defines an emerging molecular function for the Plasmodium TRAg family.

Amplicon Deep Sequencing Reveals Multiple Genetic Events Lead to Treatment Failure with Atovaquone-Proguanil in Plasmodium falciparum

Atovaquone-proguanil (AP) is used as treatment for uncomplicated malaria, and as a chemoprophylactic agent against Plasmodium falciparum. Imported malaria remains one of the top causes of fever in Canadian returning travelers. Twelve sequential whole-blood samples before and after AP treatment failure were obtained from a patient diagnosed with P. falciparum malaria upon their return from Uganda and Sudan. Ultradeep sequencing was performed on the cytb, dhfr, and dhps markers of treatment resistance before and during the episode of recrudescence. Haplotyping profiles were generated using three different approaches: msp2-3D7 agarose and capillary electrophoresis, and cpmp using amplicon deep sequencing (ADS). A complexity of infection (COI) analysis was conducted. De novo cytb Y268C mutants strains were observed during an episode of recrudescence 17 days and 16 h after the initial malaria diagnosis and AP treatment initiation. No Y268C mutant reads were observed in any of the samples prior to the recrudescence. SNPs in the dhfr and dhps genes were observed upon initial presentation. The haplotyping profiles suggest multiple clones mutating under AP selection pressure (COI > 3). Significant differences in COI were observed by capillary electrophoresis and ADS compared to the agarose gel results. ADS using cpmp revealed the lowest haplotype variation across the longitudinal analysis. Our findings highlight the value of ultra-deep sequencing methods in the understanding of P. falciparum haplotype infection dynamics. Longitudinal samples should be analyzed in genotyping studies to increase the analytical sensitivity.

Heat-shock responses: systemic and essential ways of malaria parasite survival

Fever is a part of the human innate immune response that contributes to limiting microbial growth and development in many infectious diseases. For the parasite Plasmodium falciparum, survival of febrile temperatures is crucial for its successful propagation in human populations as well as a fundamental aspect of malaria pathogenesis. This review discusses recent insights into the biological complexity of the malaria parasite's heat-shock response, which involves many cellular compartments and essential metabolic processes to alleviate oxidative stress and accumulation of damaged and unfolded proteins. We highlight the overlap between heat-shock and artemisinin resistance responses, while also explaining how the malaria parasite adapts its fever response to fight artemisinin treatment. Additionally, we discuss how this systemic and essential fight for survival can also contribute to parasite transmission to mosquitoes.

Performance of SNP barcodes to determine genetic diversity and population structure of Plasmodium falciparum in Africa

Panels of informative biallelic single nucleotide polymorphisms (SNPs) have been proposed to be an economical method to fast-track the population genetic analysis of Plasmodium falciparum in malaria-endemic areas. Whilst used successfully in low-transmission areas where infections are monoclonal and highly related, we present the first study to evaluate the performance of these 24- and 96-SNP molecular barcodes in African countries, characterised by moderate-to-high transmission, where multiclonal infections are prevalent. For SNP barcodes it is generally recommended that the SNPs chosen i) are biallelic, ii) have a minor allele frequency greater than 0.10, and iii) are independently segregating, to minimise bias in the analysis of genetic diversity and population structure. Further, to be standardised and used in many population genetic studies, these barcodes should maintain characteristics i) to iii) across various iv) geographies and v) time points. Using haplotypes generated from the MalariaGEN P. falciparum Community Project version six database, we investigated the ability of these two barcodes to fulfil these criteria in moderate-to-high transmission African populations in 25 sites across 10 countries. Predominantly clinical infections were analysed, with 52.3% found to be multiclonal, generating high proportions of mixed-allele calls (MACs) per isolate thereby impeding haplotype construction. Of the 24- and 96-SNPs, loci were removed if they were not biallelic and had low minor allele frequencies in all study populations, resulting in 20- and 75-SNP barcodes respectively for downstream population genetics analysis. Both SNP barcodes had low expected heterozygosity estimates in these African settings and consequently biased analyses of similarity. Both minor and major allele frequencies were temporally unstable. These SNP barcodes were also shown to identify weak genetic differentiation across large geographic distances based on Mantel Test and DAPC. These results demonstrate that these SNP barcodes are vulnerable to ascertainment bias and as such cannot be used as a standardised approach for malaria surveillance in moderate-to-high transmission areas in Africa, where the greatest genomic diversity of P. falciparum exists at local, regional and country levels.

The glycobiology of plasmodium falciparum: New approaches and recent advances

Like any other microorganism, pathogenic protozoan parasites rely heavily on glycoconjugates and glycan binding proteins to protect themselves from the environment and to interact with their diverse hosts. A thorough comprehension of how glycobiology contributes to the survival and virulence of these organisms may reveal unknown aspects of their biology and may open much needed avenues for the design of new strategies against them. In the case of Plasmodium falciparum, which causes the vast majority of malaria cases and deaths, the restricted variety and the simplicity of its glycans seemed to confer limited significance to the role played by glycoconjugates in the parasite. Nonetheless, the last 10 to 15 years of research are revealing a clearer and more defined picture. Thus, the use of new experimental techniques and the results obtained provide new avenues for understanding the biology of the parasite, as well as opportunities for the development of much required new tools against malaria.

Inhibitors of ApiAP2 protein DNA binding exhibit multistage activity against Plasmodium parasites

Plasmodium parasites are reliant on the Apicomplexan AP2 (ApiAP2) transcription factor family to regulate gene expression programs. AP2 DNA binding domains have no homologs in the human or mosquito host genomes, making them potential antimalarial drug targets. Using an in-silico screen to dock thousands of small molecules into the crystal structure of the AP2-EXP (Pf3D7_1466400) AP2 domain (PDB:3IGM), we identified putative AP2-EXP interacting compounds. Four compounds were found to block DNA binding by AP2-EXP and at least one additional ApiAP2 protein. Our top ApiAP2 competitor compound perturbs the transcriptome of P. falciparum trophozoites and results in a decrease in abundance of log2 fold change > 2 for 50% (46/93) of AP2-EXP target genes. Additionally, two ApiAP2 competitor compounds have multi-stage anti-Plasmodium activity against blood and mosquito stage parasites. In summary, we describe a novel set of antimalarial compounds that interact with AP2 DNA binding domains. These compounds may be used for future chemical genetic interrogation of ApiAP2 proteins or serve as starting points for a new class of antimalarial therapeutics.

Division and Transmission: Malaria Parasite Development in the Mosquito

The malaria parasite life cycle alternates between two hosts: a vertebrate and the female Anopheles mosquito vector. Cell division, proliferation, and invasion are essential for parasite development, transmission, and survival. Most research has focused on Plasmodium development in the vertebrate, which causes disease; however, knowledge of malaria parasite development in the mosquito (the sexual and transmission stages) is now rapidly accumulating, gathered largely through investigation of the rodent malaria model, with Plasmodium berghei. In this review, we discuss the seminal genome-wide screens that have uncovered key regulators of cell proliferation, invasion, and transmission during Plasmodium sexual development. Our focus is on the roles of transcription factors, reversible protein phosphorylation, and molecular motors. We also emphasize the still-unanswered important questions around key pathways in cell division during the vector transmission stages and how they may be targeted in future studies.

Systems biology of malaria explored with nonhuman primates

"The Primate Malarias" book has been a uniquely important resource for multiple generations of scientists, since its debut in 1971, and remains pertinent to the present day. Indeed, nonhuman primates (NHPs) have been instrumental for major breakthroughs in basic and pre-clinical research on malaria for over 50 years. Research involving NHPs have provided critical insights and data that have been essential for malaria research on many parasite species, drugs, vaccines, pathogenesis, and transmission, leading to improved clinical care and advancing research goals for malaria control, elimination, and eradication. Whilst most malaria scientists over the decades have been studying Plasmodium falciparum, with NHP infections, in clinical studies with humans, or using in vitro culture or rodent model systems, others have been dedicated to advancing research on Plasmodium vivax, as well as on phylogenetically related simian species, including Plasmodium cynomolgi, Plasmodium coatneyi, and Plasmodium knowlesi. In-depth study of these four phylogenetically related species over the years has spawned the design of NHP longitudinal infection strategies for gathering information about ongoing infections, which can be related to human infections. These Plasmodium-NHP infection model systems are reviewed here, with emphasis on modern systems biological approaches to studying longitudinal infections, pathogenesis, immunity, and vaccines. Recent discoveries capitalizing on NHP longitudinal infections include an advanced understanding of chronic infections, relapses, anaemia, and immune memory. With quickly emerging new technological advances, more in-depth research and mechanistic discoveries can be anticipated on these and additional critical topics, including hypnozoite biology, antigenic variation, gametocyte transmission, bone marrow dysfunction, and loss of uninfected RBCs. New strategies and insights published by the Malaria Host-Pathogen Interaction Center (MaHPIC) are recapped here along with a vision that stresses the importance of educating future experts well trained in utilizing NHP infection model systems for the pursuit of innovative, effective interventions against malaria.

In-Silico Functional Annotation of Plasmodium falciparum Hypothetical Proteins to Identify Novel Drug Targets

Plasmodium falciparum is one of the plasmodium species responsible for the majority of life-threatening malaria cases. The current antimalarial therapies are becoming less effective due to growing drug resistance, leading to the urgent requirement for alternative and more effective antimalarial drugs or vaccines. To facilitate the novel drug discovery or vaccine development efforts, recent advances in sequencing technologies provide valuable information about the whole genome of the parasite, yet a lot more needs to be deciphered due to its incomplete proteome annotation. Surprisingly, out of the 5,389 proteins currently annotated in the Plasmodium falciparum 3D7 strain, 1,626 proteins (∼30% data) are annotated as hypothetical proteins. In parasite genomic studies, the challenge to annotate hypothetical proteins is often ignored, which may obscure the crucial information related to the pathogenicity of the parasite. In this study, we attempt to characterize hypothetical proteins of the parasite to identify novel drug targets using a computational pipeline. The study reveals that out of the overall pool of the hypothetical proteins, 266 proteins have conserved functional signatures. Furthermore, the pathway analysis of these proteins revealed that 23 proteins have an essential role in various biochemical, signalling and metabolic pathways. Additionally, all the proteins (266) were subjected to computational structure analysis. We could successfully model 11 proteins. We validated and checked the structural stability of the models by performing molecular dynamics simulation. Interestingly, eight proteins show stable conformations, and seven proteins are specific for Plasmodium falciparum, based on homology analysis. Lastly, mapping the seven shortlisted hypothetical proteins on the Plasmodium falciparum protein-protein interaction network revealed 3,299 nodes and 2,750,692 edges. Our study revealed interesting functional details of seven hypothetical proteins of the parasite, which help learn more about the less-studied molecules and their interactions, providing valuable clues to unravel the role of these proteins via future experimental validation.

GIP: an open-source computational pipeline for mapping genomic instability from protists to cancer cells

Genome instability has been recognized as a key driver for microbial and cancer adaptation and thus plays a central role in many diseases. Genome instability encompasses different types of genomic alterations, yet most available genome analysis software are limited to just one type of mutation. To overcome this limitation and better understand the role of genetic changes in enhancing pathogenicity we established GIP, a novel, powerful bioinformatic pipeline for comparative genome analysis. Here, we show its application to whole genome sequencing datasets of Leishmania, Plasmodium, Candida and cancer. Applying GIP on available data sets validated our pipeline and demonstrated the power of our tool to drive biological discovery. Applied to Plasmodium vivax genomes, our pipeline uncovered the convergent amplification of erythrocyte binding proteins and identified a nullisomic strain. Re-analyzing genomes of drug adapted Candida albicans strains revealed correlated copy number variations of functionally related genes, strongly supporting a mechanism of epistatic adaptation through interacting gene-dosage changes. Our results illustrate how GIP can be used for the identification of aneuploidy, gene copy number variations, changes in nucleic acid sequences, and chromosomal rearrangements. Altogether, GIP can shed light on the genetic bases of cell adaptation and drive disease biomarker discovery.

Naturally acquired antibody kinetics against Plasmodium vivax antigens in people from a low malaria transmission region in western Thailand

BACKGROUND

Plasmodium vivax (P. vivax) is the dominant Plasmodium spp. causing the disease malaria in low-transmission regions outside of Africa. These regions often feature high proportions of asymptomatic patients with sub-microscopic parasitaemia and relapses. Naturally acquired antibody responses are induced after Plasmodium infection, providing partial protection against high parasitaemia and clinical episodes. However, previous work has failed to address the presence and maintenance of such antibody responses to P. vivax particularly in low-transmission regions.

METHODS

We followed 34 patients in western Thailand after symptomatic P. vivax infections to monitor antibody kinetics over 9 months, during which no recurrent infections occurred. We assessed total IgG, IgG subclass and IgM levels to up to 52 P. vivax proteins every 2-4 weeks using a multiplexed Luminex® assay and identified protein-specific variation in antibody longevity. Mathematical modelling was used to generate the estimated half-life of antibodies, long-, and short-lived antibody-secreting cells.

RESULTS

Generally, an increase in antibody level was observed within 1-week post symptomatic infection, followed by an exponential decay of different rates. We observed mostly IgG1 dominance and IgG3 sub-dominance in this population. IgM responses followed similar kinetic patterns to IgG, with some proteins unexpectedly inducing long-lived IgM responses. We also monitored antibody responses against 27 IgG-immunogenic antigens in 30 asymptomatic individuals from a similar region. Our results demonstrate that most antigens induced robust and long-lived total IgG responses following asymptomatic infections in the absence of (detected) boosting infections.

CONCLUSIONS

Our work provides new insights into the development and maintenance of naturally acquired immunity to P. vivax and will guide the potential use of serology to indicate immune status and/or identify populations at risk.

Population-level genome-wide STR discovery and validation for population structure and genetic diversity assessment of Plasmodium species

Short tandem repeats (STRs) are highly informative genetic markers that have been used extensively in population genetics analysis. They are an important source of genetic diversity and can also have functional impact. Despite the availability of bioinformatic methods that permit large-scale genome-wide genotyping of STRs from whole genome sequencing data, they have not previously been applied to sequencing data from large collections of malaria parasite field samples. Here, we have genotyped STRs using HipSTR in more than 3,000 Plasmodium falciparum and 174 Plasmodium vivax published whole-genome sequence data from samples collected across the globe. High levels of noise and variability in the resultant callset necessitated the development of a novel method for quality control of STR genotype calls. A set of high-quality STR loci (6,768 from P. falciparum and 3,496 from P. vivax) were used to study Plasmodium genetic diversity, population structures and genomic signatures of selection and these were compared to genome-wide single nucleotide polymorphism (SNP) genotyping data. In addition, the genome-wide information about genetic variation and other characteristics of STRs in P. falciparum and P. vivax have been available in an interactive web-based R Shiny application PlasmoSTR (https://github.com/bahlolab/PlasmoSTR).

The P. falciparum CSP repeat region contains three distinct epitopes required for protection by antibodies in vivo

Rare and potent monoclonal antibodies (mAbs) against the Plasmodium falciparum (Pf) circumsporozoite protein (CSP) on infective sporozoites (SPZ) preferentially bind the PfCSP junctional tetrapeptide NPDP or NVDP minor repeats while cross-reacting with NANP central repeats in vitro. The extent to which each of these epitopes is required for protection in vivo is unknown. Here, we assessed whether junction-, minor repeat- and central repeat-preferring human mAbs (CIS43, L9 and 317 respectively) bound and protected against in vivo challenge with transgenic P. berghei (Pb) SPZ expressing either PfCSP with the junction and minor repeats knocked out (KO), or PbCSP with the junction and minor repeats knocked in (KI). In vivo protection studies showed that the junction and minor repeats are necessary and sufficient for CIS43 and L9 to neutralize KO and KI SPZ, respectively. In contrast, 317 required major repeats for in vivo protection. These data establish that human mAbs can prevent malaria infection by targeting three different protective epitopes (NPDP, NVDP, NANP) in the PfCSP repeat region. This report will inform vaccine development and the use of mAbs to passively prevent malaria.

Human monoclonal antibodies are a promising approach for preventing malaria. Highly potent human antibodies show preferential binding to the junction or minor repeat regions of the circumsporozoite protein (CSP) of P. falciparum and cross-react to repetitive (NANP) repeats. The requirement for these binding sites for mediating protection in vivo remains unknown. Here, using transgenic P. berghei parasites expressing PfCSP containing deletions of these junctional or minor epitopes, or PbCSP containing additions of these PfCSP epitopes, we demonstrate these epitopes are necessary and sufficient to mediate protection and don’t require cross-reactivity to the NANP repeats. Our findings establish a mechanism of antibody mediated protection in vivo to prevent malaria and provide the rationale for multi-epitope vaccines to increase the breadth of the antibody responses.

Single-genome sequencing reveals within-host evolution of human malaria parasites

Population genomics of bulk malaria infections is unable to examine intrahost evolution; therefore, most work has focused on the role of recombination in generating genetic variation. We used single-cell sequencing protocol for low-parasitaemia infections to generate 406 near-complete single Plasmodium vivax genomes from 11 patients sampled during sequential febrile episodes. Parasite genomes contain hundreds of de novo mutations, showing strong signatures of selection, which are enriched in the ApiAP2 family of transcription factors, known targets of adaptation. Comparing 315 P. falciparum single-cell genomes from 15 patients with our P. vivax data, we find broad complementary patterns of de novo mutation at the gene and pathway level, revealing the importance of within-host evolution during malaria infections.

Allosteric pockets and dynamic residue network hubs of falcipain 2 in mutations including those linked to artemisinin resistance

Continually emerging resistant strains of malarial parasites to current drugs present challenges. Understanding the underlying resistance mechanisms, especially those linked to allostery is, thus, highly crucial for drug design. This forms the main concern of the paper through a case study of falcipain 2 (FP-2) and its mutations, some of which are linked to artemisinin (ART) drug resistance. Here, we applied a variety of in silico approaches and tools that we developed recently, together with existing computational tools. This included novel essential dynamics and dynamic residue network (DRN) analysis algorithms. We identified six pockets demonstrating dynamic differences in the presence of some mutations. We observed striking allosteric effects in two mutant proteins. In the presence of M245I, a cryptic pocket was detected via a unique mechanism in which Pocket 2 fused with Pocket 6. In the presence of the A353T mutation, which is located at Pocket 2, the pocket became the most rigid among all protein systems analyzed. Pocket 6 was also highly stable in all cases, except in the presence of M245I mutation. The effect of ART linked mutations was more subtle, and the changes were at residue level. Importantly, we identified an allosteric communication path formed by four unique averaged BC hubs going from the mutated residue to the catalytic site and passing through the interface of three identified pockets. Collectively, we established and demonstrated that we have robust tools and a pipeline that can be applicable to the analysis of mutations.

Single-Cell RNA Sequencing Reveals Cellular Heterogeneity and Stage Transition under Temperature Stress in Synchronized Plasmodium falciparum Cells

The malaria parasite has a complex life cycle exhibiting phenotypic and morphogenic variations in two different hosts by existing in heterogeneous developmental states. To investigate this cellular heterogeneity of the parasite within the human host, we performed single-cell RNA sequencing of synchronized Plasmodium cells under control and temperature treatment conditions. Using the Malaria Cell Atlas (https://www.sanger.ac.uk/science/tools/mca) as a guide, we identified 9 subtypes of the parasite distributed across known intraerythrocytic stages. Interestingly, temperature treatment results in the upregulation of the AP2-G gene, the master regulator of sexual development in a small subpopulation of the parasites. Moreover, we identified a heterogeneous stress-responsive subpopulation (clusters 5, 6, and 7 [∼10% of the total population]) that exhibits upregulation of stress response pathways under normal growth conditions. We also developed an online exploratory tool that will provide new insights into gene function under normal and temperature stress conditions. Thus, our study reveals important insights into cell-to-cell heterogeneity in the parasite population under temperature treatment that will be instrumental toward a mechanistic understanding of cellular adaptation and population dynamics in Plasmodium falciparum. IMPORTANCE The malaria parasite has a complex life cycle exhibiting phenotypic variations in two different hosts accompanied by cell-to-cell variability that is important for stress tolerance, immune evasion, and drug resistance. To investigate cellular heterogeneity determined by gene expression, we performed single-cell RNA sequencing (scRNA-seq) of about 12,000 synchronized Plasmodium cells under physiologically relevant normal (37°C) and temperature stress (40°C) conditions phenocopying the cyclic bouts of fever experienced during malarial infection. In this study, we found that parasites exhibit transcriptional heterogeneity in an otherwise morphologically synchronized culture. Also, a subset of parasites is continually committed to gametocytogenesis and stress-responsive pathways. These observations have important implications for understanding the mechanisms of drug resistance generation and vaccine development against the malaria parasite.

Advances and opportunities in malaria population genomics

Almost 20 years have passed since the first reference genome assemblies were published for Plasmodium falciparum, the deadliest malaria parasite, and Anopheles gambiae, the most important mosquito vector of malaria in sub-Saharan Africa. Reference genomes now exist for all human malaria parasites and nearly half of the ~40 important vectors around the world. As a foundation for genetic diversity studies, these reference genomes have helped advance our understanding of basic disease biology and drug and insecticide resistance, and have informed vaccine development efforts. Population genomic data are increasingly being used to guide our understanding of malaria epidemiology, for example by assessing connectivity between populations and the efficacy of parasite and vector interventions. The potential value of these applications to malaria control strategies, together with the increasing diversity of genomic data types and contexts in which data are being generated, raise both opportunities and challenges in the field. This Review discusses advances in malaria genomics and explores how population genomic data could be harnessed to further support global disease control efforts.

Atomic Resolution Homology Models and Molecular Dynamics Simulations of Plasmodium falciparum Tubulins

Microtubules are tubulin polymers present in the eukaryotic cytoskeleton essential for structural stability and cell division that are also roadways for intracellular transport of vesicles and organelles. In the human malaria parasite Plasmodium falciparum, apart from providing structural stability and cell division, microtubules also facilitate important biological activities crucial for parasite survival in hosts, such as egression and motility. Hence, parasite structures and processes involving microtubules are among the most important drug targets for discovering much-needed novel Plasmodium inhibitors. The current study aims to construct reliable and high-quality 3D models of α-, β-, and γ-tubulins using various modeling techniques. We identified a common binding pocket specific to Plasmodium α-, β-, and γ-tubulins. Molecular dynamics simulations confirmed the stability of the Plasmodium tubulin 3D structures. The models generated in the present study may be used for protein-protein and protein-drug interaction investigations targeted toward designing malaria parasite tubulin-specific inhibitors.

Composition and stage dynamics of mitochondrial complexes in Plasmodium falciparum

Our current understanding of mitochondrial functioning is largely restricted to traditional model organisms, which only represent a fraction of eukaryotic diversity. The unusual mitochondrion of malaria parasites is a validated drug target but remains poorly understood. Here, we apply complexome profiling to map the inventory of protein complexes across the pathogenic asexual blood stages and the transmissible gametocyte stages of Plasmodium falciparum. We identify remarkably divergent composition and clade-specific additions of all respiratory chain complexes. Furthermore, we show that respiratory chain complex components and linked metabolic pathways are up to 40-fold more prevalent in gametocytes, while glycolytic enzymes are substantially reduced. Underlining this functional switch, we find that cristae are exclusively present in gametocytes. Leveraging these divergent properties and stage dynamics for drug development presents an attractive opportunity to discover novel classes of antimalarials and increase our repertoire of gametocytocidal drugs.

Dual RNA Sequencing Meta-analysis in Plasmodium Infection Identifies Host-Parasite Interactions

Dual RNA sequencing (RNA-Seq) is the simultaneous transcriptomic analysis of interacting symbionts, for example, in malaria. Potential cross-species interactions identified by correlated gene expression might highlight interlinked signaling, metabolic, or gene regulatory pathways in addition to physically interacting proteins. Often, malaria studies address one of the interacting organisms-host or parasite-rendering the other "contamination." Here we perform a meta-analysis using such studies for cross-species expression analysis. We screened experiments for gene expression from host and Plasmodium. Out of 171 studies in Homo sapiens, Macaca mulatta, and Mus musculus, we identified 63 potential studies containing host and parasite data. While 16 studies (1,950 samples) explicitly performed dual RNA-Seq, 47 (1,398 samples) originally focused on one organism. We found 915 experimental replicates from 20 blood studies to be suitable for coexpression analysis and used orthologs for meta-analysis across different host-parasite systems. Centrality metrics from the derived gene expression networks correlated with gene essentiality in the parasites. We found indications of host immune response to elements of the Plasmodium protein degradation system, an antimalarial drug target. We identified well-studied immune responses in the host with our coexpression networks, as our approach recovers known broad processes interlinked between hosts and parasites in addition to individual host and parasite protein associations. The set of core interactions represents commonalities between human malaria and its model systems for prioritization in laboratory experiments. Our approach might also allow insights into the transferability of model systems for different pathways in malaria studies.IMPORTANCE Malaria still causes about 400,000 deaths a year and is one of the most studied infectious diseases. The disease is studied in mice and monkeys as lab models to derive potential therapeutic intervention in human malaria. Interactions between Plasmodium spp. and its hosts are either conserved across different host-parasite systems or idiosyncratic to those systems. Here we use correlation of gene expression from different RNA-Seq studies to infer common host-parasite interactions across human, mouse, and monkey studies. First, we find a set of very conserved interactors, worth further scrutiny in focused laboratory experiments. Second, this work might help assess to which extent experiments and knowledge on different pathways can be transferred from models to humans for potential therapy.

Roles of Ferredoxin-Dependent Proteins in the Apicoplast of Plasmodium falciparum Parasites

Ferredoxin (Fd) and ferredoxin-NADP+ reductase (FNR) form a redox system that is hypothesized to play a central role in the maintenance and function of the apicoplast organelle of malaria parasites. The Fd/FNR system provides reducing power to various iron-sulfur cluster (FeS)-dependent proteins in the apicoplast and is believed to help to maintain redox balance in the organelle. While the Fd/FNR system has been pursued as a target for antimalarial drug discovery, Fd, FNR, and the FeS proteins presumably reliant on their reducing power play an unknown role in parasite survival and apicoplast maintenance. To address these questions, we generated genetic deletions of these proteins in a parasite line containing an apicoplast bypass system. Through these deletions, we discovered that Fd, FNR, and certain FeS proteins are essential for parasite survival but found that none are required for apicoplast maintenance. Additionally, we addressed the question of how Fd and its downstream FeS proteins obtain FeS cofactors by deleting the FeS transfer proteins SufA and NfuApi. While individual deletions of these proteins revealed their dispensability, double deletion resulted in synthetic lethality, demonstrating a redundant role in providing FeS clusters to Fd and other essential FeS proteins. Our data support a model in which the reducing power from the Fd/FNR system to certain downstream FeS proteins is essential for the survival of blood-stage malaria parasites but not for organelle maintenance, while other FeS proteins are dispensable for this stage of parasite development. IMPORTANCE Ferredoxin (Fd) and ferredoxin-NADP+ reductase (FNR) form one of the few known redox systems in the apicoplast of malaria parasites and provide reducing power to iron-sulfur (FeS) cluster proteins within the organelle. While the Fd/FNR system has been explored as a drug target, the essentiality and roles of this system and the identity of its downstream FeS proteins have not been determined. To answer these questions, we generated deletions of these proteins in an apicoplast metabolic bypass line (PfMev) and determined the minimal set of proteins required for parasite survival. Moving upstream of this pathway, we also generated individual and dual deletions of the two FeS transfer proteins that deliver FeS clusters to Fd and downstream FeS proteins. We found that both transfer proteins are dispensable, but double deletion displayed a synthetic lethal phenotype, demonstrating their functional redundancy. These findings provide important insights into apicoplast biochemistry and drug development.

Identification of Plasmodium falciparum-specific protein PIESP2 as a novel virulence factor related to cerebral malaria

Cerebral malaria (CM) is the most severe complication caused by Plasmodium falciparum infection. The pathophysiological changes caused by parasite virulence factors and the human immune response to parasites contribute to CM. To date, very few parasite virulence proteins have been found to participate in CM. Here, we employed comparative genomics analysis and identified parasite-infected erythrocyte specific protein 2 (PIESP2) to be a CM-related protein. We conducted further experimental investigations and found that PIESP2 is an immunogenic protein. PIESP2 expression begins at the early trophozoite stage and progressively increases with parasite development. Although PIESP2 proteins mainly reside within infected red blood cells (IRBCs), some of them are present on the IRBC surface at the pigmented stage. Moreover, blockage of PIESP2 by antiserum apparently inhibited the adhesion of IRBCs to brain microvascular endothelial cells (BMECs). Western blot analysis detected the binding of PIESP2 to BMECs. Transcriptional analysis revealed that the binding of PIESP2 to BMECs can increase the expression of genes involved in the inflammatory response but decrease the expression of genes related to the anchoring junction. Overall, PIESP2 might be associated with CM by mediating the sequestration of IRBCs, inducing the inflammation response, and impairing the integrity of blood-brain barrier.

Structural and biophysical correlation of anti-NANP antibodies with in vivo protection against P. falciparum

The most advanced P. falciparum circumsporozoite protein-based malaria vaccine, RTS,S/AS01 (RTS,S), confers partial protection but with antibody titers that wane relatively rapidly, highlighting the need to elicit more potent and durable antibody responses. Here, we elucidate crystal structures, binding affinities and kinetics, and in vivo protection of eight anti-NANP antibodies derived from an RTS,S phase 2a trial and encoded by three different heavy-chain germline genes. The structures reinforce the importance of homotypic Fab-Fab interactions in protective antibodies and the overwhelmingly dominant preference for a germline-encoded aromatic residue for recognition of the NANP motif. In this study, antibody apparent affinity correlates best with protection in an in vivo mouse model, with the more potent antibodies also recognizing epitopes with repeating secondary structural motifs of type I β- and Asn pseudo 310 turns; such insights can be incorporated into design of more effective immunogens and antibodies for passive immunization.

Comprehensive Analysis of Transcript and Protein Relative Abundance During Blood Stages of Plasmodium falciparum Infection

Plasmodium falciparum is the main causative agent of human malaria. During the intraerythrocytic development cycle, the P. falciparum morphology changes dramatically from circulating young rings to sequestered mature trophozoites and schizonts. Sequestered forms contribute to the pathophysiology of severe malaria as the infected erythrocytes obstruct the microvascular flow in deep organs and induce local inflammation. However, the sequestration mechanism limits the access to the corresponding parasitic form in the clinical samples from patients infected with P. falciparum. To complement this deficiency, we aimed to evaluate the relevance of mRNA study as a proxy of protein expression in sequestered parasites. To do so, we conducted a proteotranscriptomic analysis using five independent P. falciparum laboratory strain samples. RNA sequencing was performed, and the mRNA expression level was assessed on circulating ring-stage parasites. The level of protein expression were measured by LC-MS/MS on the corresponding sequestered mature forms after 18-24 h of maturation. Overall, our results showed a strong transcriptome/transcriptome and a very strong proteome/proteome correlation between samples. Moreover, positive correlations of mRNA and protein expression levels were found between ring-stage transcriptomes and mature form proteomes. However, twice more transcripts were identified at the ring stage than proteins at the mature trophozoite stage. A high level of transcript expression did not guarantee the detection of the corresponding protein. Finally, we pointed out discrepancies at the individual gene level. Taken together, our results show that transcript and protein expressions are overall correlated. However, mRNA abundance is not a perfect proxy of protein expression at the individual level. Importantly, our study shows limitations of the "blind" use of RNA-seq and the importance of multiomics approaches for P. falciparum blood stage study in clinical samples.

Histone acetyltransferase PfGCN5 regulates stress responsive and artemisinin resistance related genes in Plasmodium falciparum

Plasmodium falciparum has evolved resistance to almost all front-line drugs including artemisinin, which threatens malaria control and elimination strategies. Oxidative stress and protein damage responses have emerged as key players in the generation of artemisinin resistance. In this study, we show that PfGCN5, a histone acetyltransferase, binds to the stress-responsive genes in a poised state and regulates their expression under stress conditions. Furthermore, we show that upon artemisinin exposure, genome-wide binding sites for PfGCN5 are increased and it is directly associated with the genes implicated in artemisinin resistance generation like BiP and TRiC chaperone. Interestingly, expression of genes bound by PfGCN5 was found to be upregulated during stress conditions. Moreover, inhibition of PfGCN5 in artemisinin-resistant parasites increases the sensitivity of the parasites to artemisinin treatment indicating its role in drug resistance generation. Together, these findings elucidate the role of PfGCN5 as a global chromatin regulator of stress-responses with a potential role in modulating artemisinin drug resistance and identify PfGCN5 as an important target against artemisinin-resistant parasites.

Comparative structural insight into prefoldin subunints of archaea and eukaryotes with special emphasis on unexplored prefoldin of Plasmodium falciparum

Prefoldin (PFD) is a heterohexameric molecular chaperone which bind unfolded proteins and subsequently deliver them to a group II chaperonin for correct folding. Although there is structural and functional information available for humans and archaea PFDs, their existence and functions in malaria parasite remains uncharacterized. In the present review, we have collected the available information on prefoldin family members of archaea and humans and attempted to analyze unexplored PFD subunits of Plasmodium falciparum (Pf). Our review enhances the understanding of probable functions, structure and mechanism of substrate binding of Pf prefoldin by comparing with the available information of its homologs in archaea and H. sapiens. Three PfPFD out of six and a Pf prefoldin-like protein are reported to be essential for parasite survival that signifies their importance in malaria parasite biology. Transcriptome analyses suggest that PfPFD subunits are up-regulated at the mRNA level during asexual and sexual stages of parasite life cycle. Our in silico analysis suggested several pivotal proteins like myosin E, cytoskeletal protein (tubulin), merozoite surface protein and ring exported protein 3 as their interacting partners. Based on structural information of archaeal and H. sapiens PFDs, P. falciparum counterparts have been modelled and key interface residues were identified that are critical for oligomerization of PfPFD subunits. We collated information on PFD-substrate binding and PFD-chaperonin interaction in detail to understand the mechanism of substrate delivery in archaea and humans. Overall, our review enables readers to view the PFD family comprehensively. Communicated by Ramaswamy H. SarmaAbbreviations: HSP: Heat shock proteins; CCT: Chaperonin containing TCP-1; PFD: Prefoldin; PFLP: Prefoldin like protein; PfPFD: Plasmodium falciparum prefoldin; Pf: Plasmodium falciparum; H. sapiens: Homo sapiens; M. thermoautotrophicus: Methanobacterium thermoautotrophicus; P. horikoshii: Pyrococcus horikoshii.

Prefoldin subunit 6 of Plasmodium falciparum binds merozoite surface protein-1 (MSP-1)

Malaria is a human disease caused by eukaryotic protozoan parasites of the Plasmodium genus. Plasmodium falciparum (Pf) causes the most lethal form of human malaria and is responsible for widespread mortality worldwide. Prefoldin is a heterohexameric molecular complex that binds and delivers unfolded proteins to chaperonin for correct folding. The prefoldin PFD6 is predicted to interact with merozoite surface protein‐1 (MSP‐1), a protein well known to play a pivotal role in erythrocyte binding and invasion by Plasmodium merozoites. We previously found that the P. falciparum (Pf) genome contains six prefoldin genes and a prefoldin‐like gene whose molecular functions are unidentified. Here, we analyzed the expression of PfPFD‐6 during the asexual blood stages of the parasite and investigated its interacting partners. PfPFD‐6 was found to be significantly expressed at the trophozoite and schizont stages. Pull‐down assays suggest PfPFD‐6 interacts with MSP‐1. In silico analysis suggested critical residues involved in the PfPFD‐6‐MSP‐1 interaction. Our data suggest PfPFD‐6 may play a role in stabilizing or trafficking MSP‐1.

PfKsgA1 functions as a transcription initiation factor and interacts with the N-terminal region of the mitochondrial RNA polymerase of Plasmodium falciparum

The small mitochondrial genome (mtDNA) of the malaria parasite is known to transcribe its genes polycistonically, although promoter element(s) have not yet been identified. An unusually large Plasmodium falciparum candidate mitochondrial phage-like RNA polymerase (PfmtRNAP) with an extended N-terminal region is encoded by the parasite nuclear genome. Using specific antibodies against the enzyme, we established that PfmtRNAP was targeted exclusively to the mitochondrion and interacted with mtDNA. Phylogenetic analysis showed that it is part of a separate apicomplexan clade. A search for PfmtRNAP-associated transcription initiation factors using sequence homology and in silico protein-protein interaction network analysis identified PfKsgA1. PfKsgA1 is a dual cytosol- and mitochondrion-targeted protein that functions as a small subunit rRNA dimethyltransferase in ribosome biogenesis. Chromatin immunoprecipitation showed that PfKsgA1 interacts with mtDNA, and in vivo crosslinking and pull-down experiments confirmed PfmtRNAP-PfKsgA1 interaction. The ability of PfKsgA1 to serve as a transcription initiation factor was demonstrated by complementation of yeast mitochondrial transcription factor Mtf1 function in Rpo41-driven in vitro transcription. Pull-down experiments using PfKsgA1 and PfmtRNAP domains indicated that the N-terminal region of PfmtRNAP interacts primarily with the PfKsgA1 C-terminal domain with some contacts being made with the linker and N-terminal domain of PfKsgA1. In the absence of full-length recombinant PfmtRNAP, solution structures of yeast mitochondrial RNA polymerase Rpo41 complexes with Mtf1 or PfKsgA1 were determined by small-angle X-ray scattering. Protein interaction interfaces thus identified matched with those reported earlier for Rpo41-Mtf1 interaction and overlaid with the PfmtRNAP-interfacing region identified experimentally for PfKsgA1. Our results indicate that in addition to a role in mitochondrial ribosome biogenesis, PfKsgA1 has an independent function as a transcription initiation factor for PfmtRNAP.

Genomic and transcriptomic evidence for descent from Plasmodium and loss of blood schizogony in Hepatocystis parasites from naturally infected red colobus monkeys

Hepatocystis is a genus of single-celled parasites infecting, amongst other hosts, monkeys, bats and squirrels. Although thought to have descended from malaria parasites (Plasmodium spp.), Hepatocystis spp. are thought not to undergo replication in the blood-the part of the Plasmodium life cycle which causes the symptoms of malaria. Furthermore, Hepatocystis is transmitted by biting midges, not mosquitoes. Comparative genomics of Hepatocystis and Plasmodium species therefore presents an opportunity to better understand some of the most important aspects of malaria parasite biology. We were able to generate a draft genome for Hepatocystis sp. using DNA sequencing reads from the blood of a naturally infected red colobus monkey. We provide robust phylogenetic support for Hepatocystis sp. as a sister group to Plasmodium parasites infecting rodents. We show transcriptomic support for a lack of replication in the blood and genomic support for a complete loss of a family of genes involved in red blood cell invasion. Our analyses highlight the rapid evolution of genes involved in parasite vector stages, revealing genes that may be critical for interactions between malaria parasites and mosquitoes.

CRISPR-TAPE: protein-centric CRISPR guide design for targeted proteome engineering

Rational molecular engineering of proteins with CRISPR-based approaches is challenged by the gene-centric nature of gRNA design tools. To address this, we have developed CRISPR-TAPE, a protein-centric gRNA design algorithm that allows users to target specific residues, or amino acid types within proteins. gRNA outputs can be customized to support maximal efficacy of homology-directed repair for engineering purposes, removing time-consuming post hoc curation, simplifying gRNA outputs and reducing CPU times.

Genome-Scale Identification of Essential Metabolic Processes for Targeting the Plasmodium Liver Stage

Plasmodium gene functions in mosquito and liver stages remain poorly characterized due to limitations in the throughput of phenotyping at these stages. To fill this gap, we followed more than 1,300 barcoded P. berghei mutants through the life cycle. We discover 461 genes required for efficient parasite transmission to mosquitoes through the liver stage and back into the bloodstream of mice. We analyze the screen in the context of genomic, transcriptomic, and metabolomic data by building a thermodynamic model of P. berghei liver-stage metabolism, which shows a major reprogramming of parasite metabolism to achieve rapid growth in the liver. We identify seven metabolic subsystems that become essential at the liver stages compared with asexual blood stages: type II fatty acid synthesis and elongation (FAE), tricarboxylic acid, amino sugar, heme, lipoate, and shikimate metabolism. Selected predictions from the model are individually validated in single mutants to provide future targets for drug development.

•

1,342 barcoded P. berghei knockout (KO) mutants analyzed for stage-specific phenotypes

•

Life-stage-specific metabolic models reveal reprogramming of cellular function

•

High agreement between blood/liver stage metabolic models and genetic screening data

•

Essential metabolic pathways for parasite development and mechanistic origin revealed

Dual Plasmepsin-Targeting Antimalarial Agents Disrupt Multiple Stages of the Malaria Parasite Life Cycle

Artemisin combination therapy (ACT) is the main treatment option for malaria, which is caused by the intracellular parasite Plasmodium. However, increased resistance to ACT highlights the importance of finding new drugs. Recently, the aspartic proteases Plasmepsin IX and X (PMIX and PMX) were identified as promising drug targets. In this study, we describe dual inhibitors of PMIX and PMX, including WM382, that block multiple stages of the Plasmodium life cycle. We demonstrate that PMX is a master modulator of merozoite invasion and direct maturation of proteins required for invasion, parasite development, and egress. Oral administration of WM382 cured mice of P. berghei and prevented blood infection from the liver. In addition, WM382 was efficacious against P. falciparum asexual infection in humanized mice and prevented transmission to mosquitoes. Selection of resistant P. falciparum in vitro was not achievable. Together, these show that dual PMIX and PMX inhibitors are promising candidates for malaria treatment and prevention.

ClinEpiDB: an open-access clinical epidemiology database resource encouraging online exploration of complex studies

The concept of open data has been gaining traction as a mechanism to increase data use, ensure that data are preserved over time, and accelerate discovery. While epidemiology data sets are increasingly deposited in databases and repositories, barriers to access still remain. ClinEpiDB was constructed as an open-access online resource for clinical and epidemiologic studies by leveraging the extensive web toolkit and infrastructure of the Eukaryotic Pathogen Database Resources (EuPathDB; a collection of databases covering 170+ eukaryotic pathogens, relevant related species, and select hosts) combined with a unified semantic web framework. Here we present an intuitive point-and-click website that allows users to visualize and subset data directly in the ClinEpiDB browser and immediately explore potential associations. Supporting study documentation aids contextualization, and data can be downloaded for advanced analyses. By facilitating access and interrogation of high-quality, large-scale data sets, ClinEpiDB aims to spur collaboration and discovery that improves global health.

ClinEpiDB: an open-access clinical epidemiology database resource encouraging online exploration of complex studies

The concept of open data has been gaining traction as a mechanism to increase data use, ensure that data are preserved over time, and accelerate discovery. While epidemiology data sets are increasingly deposited in databases and repositories, barriers to access still remain. ClinEpiDB was constructed as an open-access online resource for clinical and epidemiologic studies by leveraging the extensive web toolkit and infrastructure of the Eukaryotic Pathogen Database Resources (EuPathDB; a collection of databases covering 170+ eukaryotic pathogens, relevant related species, and select hosts) combined with a unified semantic web framework. Here we present an intuitive point-and-click website that allows users to visualize and subset data directly in the ClinEpiDB browser and immediately explore potential associations. Supporting study documentation aids contextualization, and data can be downloaded for advanced analyses. By facilitating access and interrogation of high-quality, large-scale data sets, ClinEpiDB aims to spur collaboration and discovery that improves global health.

Metabolic alterations in the erythrocyte during blood-stage development of the malaria parasite

Background

Human blood cells (erythrocytes) serve as hosts for the malaria parasite Plasmodium falciparum during its 48-h intraerythrocytic developmental cycle (IDC). Established in vitro protocols allow for the study of host–parasite interactions during this phase and, in particular, high-resolution metabolomics can provide a window into host–parasite interactions that support parasite development.

Methods

Uninfected and parasite-infected erythrocyte cultures were maintained at 2% haematocrit for the duration of the IDC, while parasitaemia was maintained at 7% in the infected cultures. The parasite-infected cultures were synchronized to obtain stage-dependent information of parasite development during the IDC. Samples were collected in quadruplicate at six time points from the uninfected and parasite-infected cultures and global metabolomics was used to analyse cell fractions of these cultures.

Results

In uninfected and parasite-infected cultures during the IDC, 501 intracellular metabolites, including 223 lipid metabolites, were successfully quantified. Of these, 19 distinct metabolites were present only in the parasite-infected culture, 10 of which increased to twofold in abundance during the IDC. This work quantified approximately five times the metabolites measured in previous studies of similar research scope, which allowed for more detailed analyses. Enrichment in lipid metabolism pathways exhibited a time-dependent association with different classes of lipids during the IDC. Specifically, enrichment occurred in sphingolipids at the earlier stages, and subsequently in lysophospholipid and phospholipid metabolites at the intermediate and end stages of the IDC, respectively. In addition, there was an accumulation of 18-, 20-, and 22-carbon polyunsaturated fatty acids, which produce eicosanoids and promote gametocytogenesis in infected erythrocyte cultures.

Conclusions

The current study revealed a number of heretofore unidentified metabolic components of the host–parasite system, which the parasite may exploit in a time-dependent manner to grow over the course of its development in the blood stage. Notably, the analyses identified components, such as precursors of immunomodulatory molecules, stage-dependent lipid dynamics, and metabolites, unique to parasite-infected cultures. These conclusions are reinforced by the metabolic alterations that were characterized during the IDC, which were in close agreement with those known from previous studies of blood-stage infection.

YihA GTPases localize to the apicoplast and mitochondrion of the malaria parasite and interact with LSU of organellar ribosomes

The YihA TRAFAC GTPases are critical for late-stage assembly of the ribosomal large subunit (LSU). In order to explore biogenesis of the reduced organellar ribosomes of the malaria parasite, we identified three nuclear-encoded homologs of YihA in Plasmodium falciparum. PfYihA1 targeted to the parasite apicoplast, PfYihA2 to the mitochondrion, and PfYihA3 was found in both the apicoplast and cytosol. The three PfYihA, expressed as recombinant proteins, were active GTPases and interacted with surrogate E. coli ribosomes in a nucleotide-independent manner. In vivo complexation of PfYihA with parasite organellar and/or cytosolic LSU was confirmed by co-immunoprecipitation using specific antibodies. Mitochondrial PfYihA2 carries a large C-ter extension with a strongly positively charged stretch. We hypothesise that this is important in compensating for the absence of helices of the central protuberance in the fragmented rRNA of Plasmodium mitoribosomes and may provide additional contact sites to aid in complex assembly. Combined with previous reports, our results indicate that P. falciparum mitochondria are likely to assemble ribosomes with the aid of PfEngA, PfObg1 and PfYihA2 GTPases while apicoplast ribosomes might use PfYihA1 and 3 in combination with other factors.

In-silico profiling and structural insights into the impact of nSNPs in the P. falciparum acetyl-CoA transporter gene to understand the mechanism of drug resistance in malaria

The continuous emergence of resistance to the available drugs poses major constraints in the development of effective therapeutics against malaria. Malaria drug resistance has been attributed to be the manifestation of numerous factors. For example, mutations in the parasite transporter protein acetyl-CoA transporter (Pfact) can remarkably affect its uptake affinity for a drug molecule against malaria, and hence enhance its susceptibility to resistance. To identify major contributors to its loss of functionality, we have thoroughly scrutinized eight such recently reported resistant mutants, via in-silico tools in terms of alterations in different properties. We performed molecular dynamics simulations of the selected Pfact mutants to gain deeper insights into the structural perturbation and dynamicity. Comparison of residue interaction network map of mutants with that of Wild type (WT) protein suggests structural changes as a result of the mutation(s) that translate into the weakening of intra-protein interactions, especially around the drug binding pocket. This, in turn, diminishes the affinity of drug molecules towards the binding site, which was validated by docking analysis. Finally, collating all the observations, we have delineated R108K mutant to deviate the most from WT protein, which, intriguingly suggests that replacing an amino acid with another of similar nature can even translate into greater functional effects as those with dissimilar substitutions.Communicated by Ramaswamy H. Sarma.

Influence of Plasmodium falciparum Calcium-Dependent Protein Kinase 5 (PfCDPK5) on the Late Schizont Stage Phosphoproteome

Protein kinases are important mediators of signal transduction in cellular pathways, and calcium-dependent protein kinases (CDPKs) compose a unique class of calcium-dependent kinases present in plants and apicomplexans, including Plasmodium parasites, the causative agents of malaria. During the asexual stage of infection, the human malaria parasite Plasmodium falciparum grows inside red blood cells, and P. falciparum calcium-dependent protein kinase 5 (PfCDPK5) is required for egress from the host cell. In this paper, we characterize the late-schizont-stage P. falciparum phosphoproteome by performing large-scale phosphoproteomic profiling on tightly synchronized parasites just prior to egress, identifying 2,704 phosphorylation sites on 919 proteins. Using a conditional knockdown of PfCDPK5, we identify 58 phosphorylation sites on 50 proteins with significant reduction in levels of PfCDPK5-deficient parasites. Furthermore, gene ontology analysis of the identified proteins reveals enrichment in transmembrane- and membrane-associated proteins and in proteins associated with transport activity. Among the identified proteins is PfNPT1, a member of the apicomplexan-specific novel putative transporter (NPT) family of proteins. We show that PfNPT1 is a potential substrate of PfCDPK5 and that PfNPT1 localizes to the parasite plasma membrane. Importantly, P. falciparum egress relies on many proteins unique to Apicomplexa that are therefore attractive targets for antimalarial therapeutics.IMPORTANCE The malaria parasite Plasmodium falciparum is a major cause of morbidity and mortality globally. The P. falciparum parasite proliferates inside red blood cells during the blood stage of infection, and egress from the red blood cell is critical for parasite survival. P. falciparum calcium-dependent protein kinase 5 (PfCDPK5) is essential for egress; parasites deficient in PfCDPK5 remain trapped inside their host cells. We have used a label-free quantitative mass spectrometry approach to identify the phosphoproteome of schizont-stage parasites just prior to egress and identify 50 proteins that display a significant reduction in phosphorylation in PfCDPK5-deficient parasites. We show that a member of the Apicomplexan-specific transport protein family, PfNPT1 is a potential substrate of PfCDPK5 and is localized to the parasite plasma membrane. P. falciparum egress requires several proteins not present in human cells, thus making this pathway an ideal target for new therapeutics.

Diverse Antibody Responses to Conserved Structural Motifs in Plasmodium falciparum Circumsporozoite Protein

Malaria vaccine candidate RTS,S/AS01 is based on the central and C-terminal regions of the circumsporozoite protein (CSP) of P. falciparum. mAb397 was isolated from a volunteer in an RTS,S/AS01 clinical trial, and it protects mice from infection by malaria sporozoites. However, mAb397 originates from the less commonly used VH3-15 germline gene compared to the VH3-30/33 antibodies generally elicited by RTS,S to the central NANP repeat region of CSP. The crystal structure of mAb397 with an NPNA₄ peptide shows that the central NPNA forms a type I β-turn and is the main recognition motif. In most anti-NANP antibodies studied to date, a germline-encoded Trp is used to engage the Pro in NPNA β-turns, but here the Trp interacts with the first Asn. This “conserved” Trp, however, can arise from different germline genes and be located in the heavy or the light chain. Variation in the terminal ψ angles of the NPNA β-turns results in different dispositions of the subsequent NPNA and, hence, different stoichiometries and modes of antibody binding to rsCSP. Diverse protective antibodies against NANP repeats are therefore not limited to a single germline gene response or mode of binding.

•

mAb397 from an RTS,S trial binds NANP repeats of malaria circumsporozoite protein.

•

Protective mAb397 is encoded by the VH3-15 gene that is rare for NANP antibodies.

•

Anti-NANP mAbs bind type I β-turns in three modes using germline-encoded Trp residues.

•

The Trp of mAb397 interacts with Asn instead of Pro in the NPNA type I β-turn.

•

Diversity in NANP immune responses can aid in next-generation vaccine design.

Decreased In Vitro Artemisinin Sensitivity of Plasmodium falciparum across India

Artemisinin-based combination therapy (ACT) has been used to treat uncomplicated Plasmodium falciparum infections in India since 2004. Since 2008, a decrease in artemisinin effectiveness has been seen throughout the Greater Mekong Subregion. The geographic proximity and ecological similarities of northeastern India to Southeast Asia may differentially affect the long-term management and sustainability of ACT in India.

Artemisinin-based combination therapy (ACT) has been used to treat uncomplicated Plasmodium falciparum infections in India since 2004. Since 2008, a decrease in artemisinin effectiveness has been seen throughout the Greater Mekong Subregion. The geographic proximity and ecological similarities of northeastern India to Southeast Asia may differentially affect the long-term management and sustainability of ACT in India. In order to collect baseline data on variations in ACT sensitivity in Indian parasites, 12 P. falciparum isolates from northeast India and 10 isolates from southwest India were studied in vitro. Ring-stage survival assay (RSA) showed reduced sensitivity to dihydroartemisinin in 50% of the samples collected in northeast India in 2014 and 2015. Two of the 10 assayed samples from the southwest region of India from as far back as 2012 also showed decreased sensitivity to artemisinin. In both these regions, kelch gene sequences were not predictive of reduced artemisinin sensitivity, as measured by RSA. The present data justify future investments in integrated approaches involving clinical follow-up studies, in vitro survival assays, and molecular markers for tracking potential changes in the effectiveness of artemisinin against P. falciparum throughout India.

Plasmodium helical interspersed subtelomeric family-an enigmatic piece of the Plasmodium biology puzzle

Plasmodium falciparum (Pf) refurbishes the infected erythrocytes by exporting a myriad of parasite proteins to the host cell. A novel exported protein family 'Plasmodium Helical Interspersed Subtelomeric' (PHIST) has gained attention for its significant roles in parasite biology. Here, we have collected and analysed available information on PHIST members to enhance understanding of their functions, varied localization and structure-function correlation. Functional diversity of PHIST proteins is highlighted by their involvement in PfEMP1 (Pf erythrocyte membrane protein 1) expression, trafficking and switching. This family also contributes to cytoadherence, gametocytogenesis, host cell modification and generation of extracellular vesicles. While the PHIST domain forms the hallmark of this family, existence and functions of additional domains (LyMP, TIGR01639) and the MEC motif underscores its diversity further. Since specific PHIST proteins seem to form pairs with PfEMP1 members, we have used in silico tools to predict such potential partners in Pf. This information and our analysis of structural data on a PHIST member provide important insights into their functioning. This review overall enables readers to view the PHIST family comprehensively, while highlighting key knowledge gaps in the field.

CX3CL1 binding protein-2 (CBP2) of Plasmodium falciparum binds nucleic acids

Several exported Plasmodium falciparum (Pf) proteins contribute to malaria biology through their involvement in cytoadherence, immune evasion and host cell remodelling. Many of these exported proteins and other host molecules are present in iRBC (infected red blood cell) generated extracellular vesicles (EVs), which are responsible for host cell modification and parasite development. CX3CL1 binding proteins (CBPs) present on the surface of iRBCs have been reported to contribute to cytoadhesion by binding with the chemokine 'CX3CL1' via their extracellular domains. Here, we have characterized the cytoplasmic domain of CBP2 to understand its function in parasite biology using biochemical and biophysical methods. Recombinant cytoplasmic CBP2 (cCBP2) binds nucleic acids showing interaction with DNA/RNA. cCBP2 shows dimer formation under non-reducing conditions highlighting the role of disulphide bonds in its oligomerization while ATP binding leads to structural changes in the protein. In vitro interaction studies depict its binding with a Maurer's cleft resident protein 'PfSBP1', which is influenced by ATP binding of cCBP2. Our results suggest CBP2 as a two-transmembrane (2TM) receptor responsible for targeting EVs and delivering cargo to host endothelial cells. We propose CBP2 as an important molecule having roles in cytoadherence and immune modulation through its extracellular and cytoplasmic domains respectively.