Search

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination., Competing Interests: No competing interests were disclosed., (Copyright: © 2021 MalariaGEN et al.)

Plasmodium vivax is responsible for the majority of malaria cases outside Africa. Unlike P. falciparum, the P. vivax life-cycle includes a dormant liver stage, the hypnozoite, which can cause infection in the absence of mosquito transmission. An effective vaccine against P. vivax blood stages would limit symptoms and pathology from such recurrent infections, and therefore could play a critical role in the control of this species. Vaccine development in P. vivax, however, lags considerably behind P. falciparum, which has many identified targets with several having transitioned to Phase II testing. By contrast only one P. vivax blood-stage vaccine candidate based on the Duffy Binding Protein (PvDBP), has reached Phase Ia, in large part because the lack of a continuous in vitro culture system for P. vivax limits systematic screening of new candidates. We used the close phylogenetic relationship between P. vivax and P. knowlesi, for which an in vitro culture system in human erythrocytes exists, to test the scalability of systematic reverse vaccinology to identify and prioritise P. vivax blood-stage targets. A panel of P. vivax proteins predicted to function in erythrocyte invasion were expressed as full-length recombinant ectodomains in a mammalian expression system. Eight of these antigens were used to generate polyclonal antibodies, which were screened for their ability to recognize orthologous proteins in P. knowlesi. These antibodies were then tested for inhibition of growth and invasion of both wild type P. knowlesi and chimeric P. knowlesi lines modified using CRISPR/Cas9 to exchange P. knowlesi genes with their P. vivax orthologues. Candidates that induced antibodies that inhibited invasion to a similar level as PvDBP were identified, confirming the utility of P. knowlesi as a model for P. vivax vaccine development and prioritizing antigens for further follow up., Competing Interests: The authors have declared that no competing interests exist.

Background: The invasion of the red blood cells by Plasmodium falciparum merozoites involves the interplay of several proteins that are also targets for vaccine development. The proteins PfRh5-PfRipr-PfCyRPA-Pfp113 assemble into a complex at the apical end of the merozoite and are together essential for erythrocyte invasion. They have also been shown to induce neutralizing antibodies and appear to be less polymorphic than other invasion-associated proteins, making them high priority blood-stage vaccine candidates. Using available whole genome sequencing data (WGS) and new capillary sequencing data (CS), this study describes the genetic polymorphism in the Rh5 complex in P. falciparum isolates obtained from Kilifi, Kenya., Methods: 162 samples collected in 2013 and 2014 were genotyped by capillary sequencing (CS) and re-analysed WGS from 68 culture-adapted P. falciparum samples obtained from a drug trial conducted from 2005 to 2007. The frequency of polymorphisms in the merozoite invasion proteins, PfRh5, PfRipr, PfCyRPA and PfP113 were examined and where possible polymorphisms co-occurring in the same isolates., Results: From a total 70 variants, including 2 indels, 19 SNPs [27.1%] were identified by both CS and WGS, while an additional 15 [21.4%] and 36 [51.4%] SNPs were identified only by either CS or WGS, respectively. All the SNPs identified by CS were non-synonymous, whereas WGS identified 8 synonymous and 47 non-synonymous SNPs. CS identified indels in repeat regions in the p113 gene in codons 275 and 859 that were not identified in the WGS data. The minor allele frequencies of the SNPs ranged between 0.7 and 34.9% for WGS and 1.1-29.6% for CS. Collectively, 12 high frequency SNPs (> 5%) were identified: four in Rh5 codon 147, 148, 203 and 429, two in p113 at codons 7 and 267 and six in Ripr codons 190, 259, 524, 985, 1003 and 1039., Conclusion: This study reveals that the majority of the polymorphisms are rare variants and confirms a low level of genetic polymorphisms in all proteins within the Rh5 complex.

Genome-scale mutagenesis screens for genes essential for apicomplexan parasite survival have been completed in three species: Plasmodium falciparum, the major human malaria parasite, Plasmodium berghei, a model rodent malaria parasite, and the more distantly related Toxoplasma gondii, the causative agent of toxoplasmosis. These three species share 2606 single-copy orthologs, 1500 of which have essentiality data in all three screens. In this review, we explore the overlap between these datasets to define the core essential genes of the phylum Apicomplexa. We further discuss the implications of these groundbreaking studies for understanding apicomplexan parasite biology, and we identify promising areas of focus for developing new pan-apicomplexan parasite interventions., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Sporozoites are a motile form of malaria-causing Plasmodium falciparum parasites that migrate from the site of transmission in the dermis through the bloodstream to invade hepatocytes. Sporozoites interact with many cells within the host, but the molecular identity of these interactions and their role in the pathology of malaria is poorly understood. Parasite proteins that are secreted and embedded within membranes are known to be important for these interactions, but our understanding of how they interact with each other to form functional complexes is largely unknown. Here, we compile a library of recombinant proteins representing the repertoire of cell surface and secreted proteins from the P. falciparum sporozoite and use an assay designed to detect extracellular interactions to systematically identify complexes. We identify three protein complexes including an interaction between two components of the p24 complex that is involved in the trafficking of glycosylphosphatidylinositol-anchored proteins through the secretory pathway. Plasmodium parasites lacking either gene are strongly inhibited in the establishment of liver-stage infections. These findings reveal an important role for the p24 complex in malaria pathogenesis and show that the library of recombinant proteins represents a valuable resource to investigate P. falciparum sporozoite biology., Competing Interests: Conflict of interest The authors declare no competing interests., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Plasmodium vivax malaria is much less common in Africa than the rest of the world because the parasite relies primarily on the Duffy antigen/chemokine receptor (DARC) to invade human erythrocytes, and the majority of Africans are Duffy negative. Recently, there has been a dramatic increase in the reporting of P. vivax cases in Africa, with a high number of them being in Duffy negative individuals, potentially indicating P. vivax has evolved an alternative invasion mechanism that can overcome Duffy negativity. Here, we analyzed single nucleotide polymorphism (SNP) and copy number variation (CNV) in Whole Genome Sequence (WGS) data from 44 P. vivax samples isolated from symptomatic malaria patients in southwestern Ethiopia, where both Duffy positive and Duffy negative individuals are found. A total of 123,711 SNPs were detected, of which 22.7% were nonsynonymous and 77.3% were synonymous mutations. The largest number of SNPs were detected on chromosomes 9 (24,007 SNPs; 19.4% of total) and 10 (16,852 SNPs, 13.6% of total). There were particularly high levels of polymorphism in erythrocyte binding gene candidates including merozoite surface protein 1 (MSP1) and merozoite surface protein 3 (MSP3.5, MSP3.85 and MSP3.9). Two genes, MAEBL and MSP3.8 related to immunogenicity and erythrocyte binding function were detected with significant signals of positive selection. Variation in gene copy number was also concentrated in genes involved in host-parasite interactions, including the expansion of the Duffy binding protein gene (PvDBP) on chromosome 6 and MSP3.11 on chromosome 10. Based on the phylogeny constructed from the whole genome sequences, the expansion of these genes was an independent process among the P. vivax lineages in Ethiopia. We further inferred transmission patterns of P. vivax infections among study sites and showed various levels of gene flow at a small geographical scale. The genomic features of P. vivax provided baseline data for future comparison with those in Duffy-negative individuals and allowed us to develop a panel of informative Single Nucleotide Polymorphic markers diagnostic at a micro-geographical scale., Competing Interests: The authors have declared that no competing interests exist.

Plasmodium vivax gene regulation remains difficult to study due to the lack of a robust in vitro culture method, low parasite densities in peripheral circulation and asynchronous parasite development. We adapted an RNA-seq protocol "DAFT-seq" to sequence the transcriptome of four P. vivax field isolates that were cultured for a short period ex vivo before using a density gradient for schizont enrichment. Transcription was detected from 78% of the PvP01 reference genome, despite being schizont-enriched samples. This extensive data was used to define thousands of 5' and 3' untranslated regions, some of which overlapped with neighbouring transcripts, and to improve the gene models of 352 genes, including identifying 20 novel gene transcripts. This dataset has also significantly increased the known amount of heterogeneity between P. vivax schizont transcriptomes from individual patients. The majority of genes found to be differentially expressed between the isolates lack Plasmodium falciparum homologs and are predicted to be involved in host-parasite interactions, with an enrichment in reticulocyte binding proteins, merozoite surface proteins and exported proteins with unknown function. An improved understanding of the diversity within P. vivax transcriptomes will be essential for the prioritisation of novel vaccine targets.

Malaria has had a major effect on the human genome, with many protective polymorphisms-such as the sickle-cell trait-having been selected to high frequencies in malaria-endemic regions 1,2 . The blood group variant Dantu provides 74% protection against all forms of severe malaria in homozygous individuals 3-5 , a similar degree of protection to that afforded by the sickle-cell trait and considerably greater than that offered by the best malaria vaccine. Until now, however, the protective mechanism has been unknown. Here we demonstrate the effect of Dantu on the ability of the merozoite form of the malaria parasite Plasmodium falciparum to invade red blood cells (RBCs). We find that Dantu is associated with extensive changes to the repertoire of proteins found on the RBC surface, but, unexpectedly, inhibition of invasion does not correlate with specific RBC-parasite receptor-ligand interactions. By following invasion using video microscopy, we find a strong link between RBC tension and merozoite invasion, and identify a tension threshold above which invasion rarely occurs, even in non-Dantu RBCs. Dantu RBCs have higher average tension than non-Dantu RBCs, meaning that a greater proportion resist invasion. These findings provide both an explanation for the protective effect of Dantu, and fresh insight into why the efficiency of P. falciparum invasion might vary across the heterogenous populations of RBCs found both within and between individuals.

The recurrent emergence of drug resistance in Plasmodium falciparum increases the urgency to genetically validate drug resistance mechanisms and identify new targets. Reverse genetics have facilitated genome-scale knockout screens in Plasmodium berghei and Toxoplasma gondii, in which pooled transfections of multiple vectors were critical to increasing scale and throughput. These approaches have not yet been implemented in human malaria species such as P. falciparum and P. knowlesi, in part because the extent to which pooled transfections can be performed in these species remains to be evaluated. Here we use next-generation sequencing to quantitate uptake of a pool of 94 barcoded vectors. The distribution of vector acquisition allowed us to estimate the number of barcodes and DNA molecules taken up by the parasite population. Dilution cloning of P. falciparum transfectants showed that individual clones possess as many as seven episomal barcodes, revealing that an intake of multiple vectors is a frequent event despite the inefficient transfection efficiency. Transfection of three spectrally-distinct fluorescent reporters allowed us to evaluate different transfection methods and revealed that schizont-stage transfection limited the tendency for parasites to take up multiple vectors. In contrast to P. falciparum, we observed that the higher transfection efficiency of P. knowlesi resulted in near complete representation of the library. These findings have important implications for how reverse genetics can be scaled in culturable Plasmodium species.

The rodent parasite Plasmodium chabaudi is an important in vivo model of malaria. The ability to produce chronic infections makes it particularly useful for investigating the development of anti- Plasmodium immunity, as well as features associated with parasite virulence during both the acute and chronic phases of infection. P. chabaudi also undergoes asexual maturation (schizogony) and erythrocyte invasion in culture, so offers an experimentally-amenable in vivo to in vitro model for studying gene function and drug activity during parasite replication. To extend the usefulness of this model, we have further optimised transfection protocols and plasmids for P. chabaudi and generated stable, fluorescent lines that are free from drug-selectable marker genes. These mother-lines show the same infection dynamics as wild-type parasites throughout the lifecycle in mice and mosquitoes; furthermore, their virulence can be increased by serial blood passage and reset by mosquito transmission. We have also adapted the large-insert, linear Plasmo GEM vectors that have revolutionised the scale of experimental genetics in another rodent malaria parasite and used these to generate barcoded P. chabaudi gene-deletion and -tagging vectors for transfection in our fluorescent P. chabaudi mother-lines. This produces a tool-kit of P. chabaudi lines, vectors and transfection approaches that will be of broad utility to the research community., Competing Interests: No competing interests were disclosed., (Copyright: © 2020 Marr EJ et al.)

Background: Plasmodium parasites undergo several major developmental transitions during their complex lifecycle, which are enabled by precisely ordered gene expression programs. Transcriptomes from the 48-h blood stages of the major human malaria parasite Plasmodium falciparum have been described using cDNA microarrays and RNA-seq, but these assays have not always performed well within non-coding regions, where the AT-content is often 90-95%., Results: We developed a directional, amplification-free RNA-seq protocol (DAFT-seq) to reduce bias against AT-rich cDNA, which we have applied to three strains of P. falciparum (3D7, HB3 and IT). While strain-specific differences were detected, overall there is strong conservation between the transcriptional profiles. For the 3D7 reference strain, transcription was detected from 89% of the genome, with over 78% of the genome transcribed into mRNAs. We also find that transcription from bidirectional promoters frequently results in non-coding, antisense transcripts. These datasets allowed us to refine the 5' and 3' untranslated regions (UTRs), which can be variable, long (> 1000 nt), and often overlap those of adjacent transcripts., Conclusions: The approaches applied in this study allow a refined description of the transcriptional landscape of P. falciparum and demonstrate that very little of the densely packed P. falciparum genome is inactive or redundant. By capturing the 5' and 3' ends of mRNAs, we reveal both constant and dynamic use of transcriptional start sites across the intraerythrocytic developmental cycle that will be useful in guiding the definition of regulatory regions for use in future experimental gene expression studies.

Nearly half of the genes in the Plasmodium falciparum genome have not yet been functionally investigated. We used homology-based structural modeling to identify multiple copies of Armadillo repeats within one uncharacterized gene expressed during the intraerythrocytic stages, PF3D7_0410600, subsequently referred to as P. falciparum Armadillo-Type Repeat Protein (PfATRP). Soluble recombinant PfATRP was expressed in a bacterial expression system, purified to apparent homogeneity and the identity of the recombinant PfATRP was confirmed by mass spectrometry. Affinity-purified α-PfATRP rabbit antibodies specifically recognized the recombinant protein. Immunofluorescence assays revealed that α-PfATRP rabbit antibodies reacted with P. falciparum schizonts. Anti-PfATRP antibody exhibited peripheral staining patterns around the merozoites. Given the localization of PfATRP in merozoites, we tested for an egress phenotype during schizont arrest assays and demonstrated that native PfATRP is inaccessible on the surface of merozoites in intact schizonts. Dual immunofluorescence assays with markers for the inner membrane complex (IMC) and microtubules suggest partial colocalization in both asexual and sexual stage parasites. Using the soluble recombinant PfATRP in a screen of plasma samples revealed that malaria-infected children have naturally acquired PfATRP-specific antibodies., (Copyright © 2020 Amlabu, Ilani, Opoku, Nyarko, Quansah, Thiam, Anim, Ayivor-Djanie, Akuh, Mensah-Brown, Rayner and Awandare.)

As malaria control programmes concentrate their efforts towards malaria elimination a better understanding of malaria transmission patterns at fine spatial resolution units becomes necessary. Defining spatial units that consider transmission heterogeneity, human movement and migration will help to set up achievable malaria elimination milestones and guide the creation of efficient operational administrative control units. Using a combination of genetic and epidemiological data we defined a malaria transmission unit as the area contributing 95% of malaria cases diagnosed at the catchment facility located in the town of Guapi in the South Pacific Coast of Colombia. We provide data showing that P. falciparum malaria transmission is heterogeneous in time and space and analysed, using topological data analysis, the spatial connectivity, at the micro epidemiological level, between parasite populations circulating within the unit. To illustrate the necessity to evaluate the efficacy of malaria control measures within the transmission unit in order to increase the efficiency of the malaria control effort, we provide information on the size of the asymptomatic reservoir, the nature of parasite genotypes associated with drug resistance as well as the frequency of the Pfhrp2/3 deletion associated with false negatives when using Rapid Diagnostic Tests.

Recent genomic studies are investigating the wide-ranging implications of malaria complex infections on parasite diversity, transmission, and downstream disease eradication efforts. Using single cell sequencing, new work by Nkhoma et al. provides evidence that there is unexpectedly frequent co-transmission of related parasites in the intense transmission setting in Malawi., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

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., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

Plasmodium species are frequently host-specific, but little is currently known about the molecular factors restricting host switching. This is particularly relevant for P. falciparum, the only known human-infective species of the Laverania sub-genus, all other members of which infect African apes. Here we show that all tested P. falciparum isolates contain an inactivating mutation in an erythrocyte invasion associated gene, PfEBA165, the homologues of which are intact in all ape-infective Laverania species. Recombinant EBA165 proteins only bind ape, not human, erythrocytes, and this specificity is due to differences in erythrocyte surface sialic acids. Correction of PfEBA165 inactivating mutations by genome editing yields viable parasites, but is associated with down regulation of both PfEBA165 and an adjacent invasion ligand, which suggests that PfEBA165 expression is incompatible with parasite growth in human erythrocytes. Pseudogenization of PfEBA165 may represent a key step in the emergence and evolution of P. falciparum.

Plasmodium vivax invasion of human erythrocytes depends on the Duffy Binding Protein (PvDBP) which interacts with the Duffy antigen. PvDBP copy number has been recently shown to vary between P. vivax isolates in Sub-Saharan Africa. However, the extent of PvDBP copy number variation, the type of PvDBP multiplications, as well as its significance across broad samples are still unclear. We determined the prevalence and type of PvDBP duplications, as well as PvDBP copy number variation among 178 Ethiopian P. vivax isolates using a PCR-based diagnostic method, a novel quantitative real-time PCR assay and whole genome sequencing. For the 145 symptomatic samples, PvDBP duplications were detected in 95 isolates, of which 81 had the Cambodian and 14 Malagasy-type PvDBP duplications. PvDBP varied from 1 to >4 copies. Isolates with multiple PvDBP copies were found to be higher in symptomatic than asymptomatic infections. For the 33 asymptomatic samples, PvDBP was detected with two copies in two of the isolates, and both were the Cambodian-type PvDBP duplication. PvDBP copy number in Duffy-negative heterozygotes was not significantly different from that in Duffy-positives, providing no support for the hypothesis that increased copy number is a specific association with Duffy-negativity, although the number of Duffy-negatives was small and further sampling is required to test this association thoroughly., Competing Interests: The authors have declared that no competing interests exist.

Malaria parasites adopt a remarkable variety of morphological life stages as they transition through multiple mammalian host and mosquito vector environments. We profiled the single-cell transcriptomes of thousands of individual parasites, deriving the first high-resolution transcriptional atlas of the entire Plasmodium berghei life cycle. We then used our atlas to precisely define developmental stages of single cells from three different human malaria parasite species, including parasites isolated directly from infected individuals. The Malaria Cell Atlas provides both a comprehensive view of gene usage in a eukaryotic parasite and an open-access reference dataset for the study of malaria parasites., (Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Malaria represents a major global health issue, and the identification of new intervention targets remains an urgent priority. This search is hampered by more than one-third of the genes of malaria-causing Plasmodium parasites being uncharacterized. We report a large-scale protein interaction network in Plasmodium schizonts, generated by combining blue native-polyacrylamide electrophoresis with quantitative mass spectrometry and machine learning. This integrative approach, spanning 3 species, identifies >20,000 putative protein interactions, organized into 600 protein clusters. We validate selected interactions, assigning functions in chromatin regulation to previously unannotated proteins and suggesting a role for an EELM2 domain-containing protein and a putative microrchidia protein as mechanistic links between AP2-domain transcription factors and epigenetic regulation. Our interactome represents a high-confidence map of the native organization of core cellular processes in Plasmodium parasites. The network reveals putative functions for uncharacterized proteins, provides mechanistic and structural insight, and uncovers potential alternative therapeutic targets., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

Plasmodium vivax causes the majority of malaria outside Africa, but is poorly understood at a cellular level partly due to technical difficulties in maintaining it in in vitro culture conditions. In the past decades, drug resistant P. vivax parasites have emerged, mainly in Southeast Asia, but while some molecular markers of resistance have been identified, none have so far been confirmed experimentally, which limits interpretation of the markers, and hence our ability to monitor and control the spread of resistance. Some of these potential markers have been identified through P. vivax genome-wide population genetic analyses, which highlighted genes under recent evolutionary selection in Southeast Asia, where chloroquine resistance is most prevalent. These genes could be involved in drug resistance, but no experimental proof currently exists to support this hypothesis. In this study, we used Plasmodium knowlesi, the most closely related species to P. vivax that can be cultured in human erythrocytes, as a model system to express P. vivax genes and test for their role in drug resistance. We adopted a strategy of episomal expression, and were able to express fourteen P. vivax genes, including two allelic variants of several hypothetical resistance genes. Their expression level and localisation were assessed, confirming cellular locations conjectured from orthologous species, and suggesting locations for several previously unlocalised proteins, including an apical location for PVX_101445. These findings establish P. knowlesi as a suitable model for P. vivax protein expression. We performed chloroquine and mefloquine drug assays, finding no significant differences in drug sensitivity: these results could be due to technical issues, or could indicate that these genes are not actually involved in drug resistance, despite being under positive selection pressure in Southeast Asia. These data confirm that in vitro P. knowlesi is a useful tool for studying P. vivax biology. Its close evolutionary relationship to P. vivax, high transfection efficiency, and the availability of markers for colocalisation, all make it a powerful model system. Our study is the first of its kind using P. knowlesi to study unknown P. vivax proteins and investigate drug resistance mechanisms., Competing Interests: The authors have declared that no competing interests exist.

Passive transfer studies in humans clearly demonstrated the protective role of IgG antibodies against malaria. Identifying the precise parasite antigens that mediate immunity is essential for vaccine design, but has proved difficult. Completion of the Plasmodium falciparum genome revealed thousands of potential vaccine candidates, but a significant bottleneck remains in their validation and prioritization for further evaluation in clinical trials. Focusing initially on the Plasmodium falciparum merozoite proteome, we used peer-reviewed publications, multiple proteomic and bioinformatic approaches, to select and prioritize potential immune targets. We expressed 109 P. falciparum recombinant proteins, the majority of which were obtained using a mammalian expression system that has been shown to produce biologically functional extracellular proteins, and used them to create KILchip v1.0: a novel protein microarray to facilitate high-throughput multiplexed antibody detection from individual samples. The microarray assay was highly specific; antibodies against P. falciparum proteins were detected exclusively in sera from malaria-exposed but not malaria-naïve individuals. The intensity of antibody reactivity varied as expected from strong to weak across well-studied antigens such as AMA1 and RH5 (Kruskal-Wallis H test for trend: p < 0.0001). The inter-assay and intra-assay variability was minimal, with reproducible results obtained in re-assays using the same chip over a duration of 3 months. Antibodies quantified using the multiplexed format in KILchip v1.0 were highly correlated with those measured in the gold-standard monoplex ELISA [median (range) Spearman's R of 0.84 (0.65-0.95)]. KILchip v1.0 is a robust, scalable and adaptable protein microarray that has broad applicability to studies of naturally acquired immunity against malaria by providing a standardized tool for the detection of antibody correlates of protection. It will facilitate rapid high-throughput validation and prioritization of potential Plasmodium falciparum merozoite-stage antigens paving the way for urgently needed clinical trials for the next generation of malaria vaccines.

Background: Malaria parasites are genetically polymorphic and phenotypically plastic. In studying transcriptome variation among parasites from different infections, it is challenging to overcome potentially confounding technical and biological variation between samples. We investigate variation in the major human parasite Plasmodium falciparum, generating RNA-seq data on multiple independent replicate sample preparations of merozoite-containing intra-erythrocytic schizonts from a panel of clinical isolates and from long-term laboratory-adapted clones, with a goal of robustly identifying differentially expressed genes., Results: Analysis of biological sample replicates shows that increased numbers improve the true discovery rate of differentially expressed genes, and that six independent replicates of each parasite line allowed identification of most differences that could be detected with larger numbers. For highly expressed genes, focusing on the top quartile at schizont stages, there was more power to detect differences. Comparing cultured clinical isolates and laboratory-adapted clones, genes more highly expressed in the laboratory-adapted clones include those encoding an AP2 transcription factor (PF3D7&#95;0420300), a ubiquitin-binding protein and two putative methyl transferases. In contrast, higher expression in clinical isolates was seen for the merozoite surface protein gene dblmsp2, proposed to be a marker of schizonts forming merozoites committed to sexual differentiation. Variable expression was extremely strongly, but not exclusively, associated with genes known to be targeted by Heterochromatin Protein 1. Clinical isolates show variable expression of several known merozoite invasion ligands, as well as other genes for which new RT-qPCR assays validate the quantitation and allow characterisation in samples with more limited material. Expression levels of these genes vary among schizont preparations of different clinical isolates in the first ex vivo cycle in patient erythrocytes, but mean levels are similar to those in continuously cultured clinical isolates., Conclusions: Analysis of multiple biological sample replicates greatly improves identification of genes variably expressed between different cultured parasite lines. Clinical isolates recently established in culture show differences from long-term adapted clones in transcript levels of particular genes, and are suitable for analyses requiring biological replicates to understand parasite phenotypes and variable expression likely to be relevant in nature.

Invasion of human erythrocytes is essential for Plasmodium falciparum parasite survival and pathogenesis, and is also a complex phenotype. While some later steps in invasion appear to be invariant and essential, the earlier steps of recognition are controlled by a series of redundant, and only partially understood, receptor-ligand interactions. Reverse genetic analysis of laboratory adapted strains has identified multiple genes that when deleted can alter invasion, but how the relative contributions of each gene translate to the phenotypes of clinical isolates is far from clear. We used a forward genetic approach to identify genes responsible for variable erythrocyte invasion by phenotyping the parents and progeny of previously generated experimental genetic crosses. Linkage analysis using whole genome sequencing data revealed a single major locus was responsible for the majority of phenotypic variation in two invasion pathways. This locus contained the PfRh2a and PfRh2b genes, members of one of the major invasion ligand gene families, but not widely thought to play such a prominent role in specifying invasion phenotypes. Variation in invasion pathways was linked to significant differences in PfRh2a and PfRh2b expression between parasite lines, and their role in specifying alternative invasion was confirmed by CRISPR-Cas9-mediated genome editing. Expansion of the analysis to a large set of clinical P. falciparum isolates revealed common deletions, suggesting that variation at this locus is a major cause of invasion phenotypic variation in the endemic setting. This work has implications for blood-stage vaccine development and will help inform the design and location of future large-scale studies of invasion in clinical isolates., Competing Interests: The authors have declared that no competing interests exist.

Wild-living African apes are endemically infected with parasites that are closely related to human Plasmodium vivax , a leading cause of malaria outside Africa. This finding suggests that the origin of P. vivax was in Africa, even though the parasite is now rare in humans there. To elucidate the emergence of human P. vivax and its relationship to the ape parasites, we analyzed genome sequence data of P. vivax strains infecting six chimpanzees and one gorilla from Cameroon, Gabon, and Côte d'Ivoire. We found that ape and human parasites share nearly identical core genomes, differing by only 2% of coding sequences. However, compared with the ape parasites, human strains of P. vivax exhibit about 10-fold less diversity and have a relative excess of nonsynonymous nucleotide polymorphisms, with site-frequency spectra suggesting they are subject to greatly relaxed purifying selection. These data suggest that human P. vivax has undergone an extreme bottleneck, followed by rapid population expansion. Investigating potential host-specificity determinants, we found that ape P. vivax parasites encode intact orthologs of three reticulocyte-binding protein genes ( rbp2d , rbp2e , and rbp3 ), which are pseudogenes in all human P. vivax strains. However, binding studies of recombinant RBP2e and RBP3 proteins to human, chimpanzee, and gorilla erythrocytes revealed no evidence of host-specific barriers to red blood cell invasion. These data suggest that, from an ancient stock of P. vivax parasites capable of infecting both humans and apes, a severely bottlenecked lineage emerged out of Africa and underwent rapid population growth as it spread globally., Competing Interests: The authors declare no conflict of interest., (Copyright © 2018 the Author(s). Published by PNAS.)

The development of highly effective and durable vaccines against the human malaria parasites Plasmodium falciparum and P. vivax remains a key priority. Decades of endeavor have taught that achieving this goal will be challenging; however, recent innovation in malaria vaccine research and a diverse pipeline of novel vaccine candidates for clinical assessment provides optimism. With first-generation pre-erythrocytic vaccines aiming for licensure in the coming years, it is important to reflect on how next-generation approaches can improve on their success. Here we review the latest vaccine approaches that seek to prevent malaria infection, disease, and transmission and highlight some of the major underlying immunological and molecular mechanisms of protection. The synthesis of rational antigen selection, immunogen design, and immunization strategies to induce quantitatively and qualitatively improved immune effector mechanisms offers promise for achieving sustained high-level protection., (Copyright © 2018 The Author(s). Published by Elsevier Inc. All rights reserved.)

The malaria parasite Plasmodium falciparum has evolved an unusual genome structure. The majority of the genome is relatively stable, with mutation rates similar to most eukaryotic species. However, some regions are very unstable with high recombination rates, driving the generation of new immune evasion-associated var genes. The molecular factors controlling the inconsistent stability of this genome are not known. Here we studied the roles of the two putative RecQ helicases in P. falciparum, PfBLM and PfWRN. When PfWRN was knocked down, recombination rates increased four-fold, generating chromosomal abnormalities, a high rate of chimeric var genes and many microindels, particularly in known 'fragile sites'. This is the first identification of a gene involved in suppressing recombination and maintaining genome stability in Plasmodium. By contrast, no change in mutation rate appeared when the second RecQ helicase, PfBLM, was mutated. At the transcriptional level, however, both helicases evidently modulate the transcription of large cohorts of genes, with several hundred genes-including a large proportion of vars-showing deregulated expression in each RecQ mutant. Aberrant processing of stalled replication forks is a possible mechanism underlying elevated mutation rates and this was assessed by measuring DNA replication dynamics in the RecQ mutant lines. Replication forks moved slowly and stalled at elevated rates in both mutants, confirming that RecQ helicases are required for efficient DNA replication. Overall, this work identifies the Plasmodium RecQ helicases as major players in DNA replication, antigenic diversification and genome stability in the most lethal human malaria parasite, with important implications for genome evolution in this pathogen., Competing Interests: The authors have declared that no competing interests exist.

Malaria has been one of the strongest selective forces on the human genome. The increased frequency of haemoglobinopathies, as well as numerous other blood groups, in malaria endemic regions is commonly attributed to a protective effect of these alleles against malaria. In the majority of these cases however there have been no systematic functional studies to test protective mechanisms, in large part because most host-parasite interaction assays are not quantitative or scalable. We describe the development of an erythrocyte preference assay which uses differential labelling with fluorescent dyes to distinguish invasion into four different erythrocyte populations which are all co-incubated with a single Plasmodium falciparum parasite culture. Testing this assay on erythrocytes across the ABO blood system from forty independent donors reveals for the first time that P. falciparum parasites preferentially invade group O over Group A erythrocytes. This runs counter to the known protective effect of group O against severe malaria, but emphasises the complexities of host-pathogen interactions, and the need for highly quantitative and scalable assays to systematically explore them.

Severe malaria is caused by the apicomplexan parasite Plasmodium falciparum. Despite decades of research, the distinct biology of these parasites has made it challenging to establish high-throughput genetic approaches to identify and prioritize therapeutic targets. Using transposon mutagenesis of P. falciparum in an approach that exploited its AT-rich genome, we generated more than 38,000 mutants, saturating the genome and defining mutability and fitness costs for over 87% of genes. Of 5399 genes, our study defined 2680 genes as essential for optimal growth of asexual blood stages in vitro. These essential genes are associated with drug resistance, represent leading vaccine candidates, and include approximately 1000 Plasmodium -conserved genes of unknown function. We validated this approach by testing proteasome pathways for individual mutants associated with artemisinin sensitivity., (Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

The symptoms of malaria are brought about by blood-stage parasites, which are established when merozoites invade human erythrocytes. Our understanding of the molecular events that underpin erythrocyte invasion remains hampered by the short-period of time that merozoites are invasive. To address this challenge, a Plasmodium falciparum gamma-irradiated long-lived merozoite (LLM) line was developed and investigated. Purified LLMs invaded erythrocytes by an increase of 10-300 fold compared to wild-type (WT) merozoites. Using an integrated omics approach, we investigated the basis for the phenotypic difference. Only a few single nucleotide polymorphisms within the P. falciparum genome were identified and only marginal differences were observed in the merozoite transcriptomes. By contrast, using label-free quantitative mass-spectrometry, a significant change in protein abundance was noted, of which 200 were proteins of unknown function. We determined the relative molar abundance of over 1100 proteins in LLMs and further characterized the major merozoite surface protein complex. A unique processed MSP1 intermediate was identified in LLM but not observed in WT suggesting that delayed processing may be important for the observed phenotype. This integrated approach has demonstrated the significant role of the merozoite proteome during erythrocyte invasion, while identifying numerous unknown proteins likely to be involved in invasion.

A highly effective vaccine would be a valuable weapon in the drive toward malaria elimination. No such vaccine currently exists, and only a handful of the hundreds of potential candidates in the parasite genome have been evaluated. In this study, we systematically evaluated 29 antigens likely to be involved in erythrocyte invasion, an essential developmental stage during which the malaria parasite is vulnerable to antibody-mediated inhibition. Testing antigens alone and in combination identified several strain-transcending targets that had synergistic combinatorial effects in vitro, while studies in an endemic population revealed that combinations of the same antigens were associated with protection from febrile malaria. Video microscopy established that the most effective combinations targeted multiple discrete stages of invasion, suggesting a mechanistic explanation for synergy. Overall, this study both identifies specific antigen combinations for high-priority clinical testing and establishes a generalizable approach that is more likely to produce effective vaccines., Competing Interests: Conflict of interest statement: Three of the authors (J.C.R., G.J.W., and L.Y.B.) are named on patent applications relating to some of the vaccine candidates studied in this work., (Copyright © 2017 the Author(s). Published by PNAS.)

The study of antigenic targets of naturally-acquired immunity is essential to identify and prioritize antigens for further functional characterization. We measured total IgG antibodies to 38 P. vivax antigens, investigating their relationship with prospective risk of malaria in a cohort of 1-3 years old Papua New Guinean children. Using simulated annealing algorithms, the potential protective efficacy of antibodies to multiple antigen-combinations, and the antibody thresholds associated with protection were investigated for the first time. High antibody levels to multiple known and newly identified proteins were strongly associated with protection (IRR 0.44-0.74, p<0.001-0.041). Among five-antigen combinations with the strongest protective effect (>90%), EBP, DBPII, RBP1a, CyRPA, and PVX&#95;081550 were most frequently identified; several of them requiring very low antibody levels to show a protective association. These data identify individual antigens that should be prioritized for further functional testing and establish a clear path to testing a multicomponent P. vivax vaccine.

Two decades after the first Plasmodium transfection, attempts have been made to disrupt more than 3,151 genes in malaria parasites, across five Plasmodium species. While results from rodent malaria transfections have been curated and systematised, empowering large-scale analysis, phenotypic data from human malaria parasite transfections currently exists as individual reports scattered across a the literature. To facilitate systematic analysis of published experimental genetic data across Plasmodium species, we have built PhenoPlasm ( http://www.phenoplasm.org), a database of phenotypes generated by transfection experiments in all Plasmodium parasites. The site provides a simple interface linking citation-backed Plasmodium reverse-genetic phenotypes to gene IDs. The database has been populated with phenotypic data on 367 P. falciparum genes, curated from 176 individual publications, as well as existing data on rodent Plasmodium species from RMgmDB and PlasmoGEM. This is the first time that all available data on P. falciparum transfection experiments has been brought together in a single place. These data are presented using ortholog mapping to allow a researcher interested in a gene in one species to see results across other Plasmodium species. The collaborative nature of the database enables any researcher to add new phenotypes as they are discovered. As an example of database utility, we use the currently available datasets to identify RAP (RNA-binding domain abundant in Apicomplexa)-domain containing proteins as crucial to parasite survival., Competing Interests: Competing interests: No competing interests were disclosed.

Malaria is caused by five different Plasmodium spp. in humans each of which modifies the host erythrocyte to survive and replicate. The two main causes of malaria, P. falciparum and P. vivax, differ in their ability to cause severe disease, mainly due to differences in the cytoadhesion of infected erythrocytes (IE) in the microvasculature. Cytoadhesion of P. falciparum in the brain leads to a large number of deaths each year and is a consequence of exported parasite proteins, some of which modify the erythrocyte cytoskeleton while others such as PfEMP1 project onto the erythrocyte surface where they bind to endothelial cells. Here we investigate the effects of knocking out an exported Hsp70-type chaperone termed Hsp70-x that is present in P. falciparum but not P. vivax. Although the growth of Δhsp70-x parasites was unaffected, the export of PfEMP1 cytoadherence proteins was delayed and Δhsp70-x IE had reduced adhesion. The Δhsp70-x IE were also more rigid than wild-type controls indicating changes in the way the parasites modified their host erythrocyte. To investigate the cause of this, transcriptional and translational changes in exported and chaperone proteins were monitored and some changes were observed. We propose that PfHsp70-x is not essential for survival in vitro, but may be required for the efficient export and functioning of some P. falciparum exported proteins.

The genomes of malaria parasites contain many genes of unknown function. To assist drug development through the identification of essential genes and pathways, we have measured competitive growth rates in mice of 2,578 barcoded Plasmodium berghei knockout mutants, representing >50% of the genome, and created a phenotype database. At a single stage of its complex life cycle, P. berghei requires two-thirds of genes for optimal growth, the highest proportion reported from any organism and a probable consequence of functional optimization necessitated by genomic reductions during the evolution of parasitism. In contrast, extreme functional redundancy has evolved among expanded gene families operating at the parasite-host interface. The level of genetic redundancy in a single-celled organism may thus reflect the degree of environmental variation it experiences. In the case of Plasmodium parasites, this helps rationalize both the relative successes of drugs and the greater difficulty of making an effective vaccine., (Copyright © 2017 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: Thailand is aiming to eliminate malaria by the year 2024. Plasmodium vivax has now become the dominant species causing malaria within the country, and a high proportion of infections are asymptomatic. A better understanding of antibody dynamics to P. vivax antigens in a low-transmission setting, where acquired immune responses are poorly characterized, will be pivotal for developing new strategies for elimination, such as improved surveillance methods and vaccines. The objective of this study was to characterize total IgG antibody levels to 11 key P. vivax proteins in a village of western Thailand., Methods: Plasma samples from 546 volunteers enrolled in a cross-sectional survey conducted in 2012 in Kanchanaburi Province were utilized. Total IgG levels to 11 different proteins known or predicted to be involved in reticulocyte binding or invasion (ARP, GAMA, P41, P12, PVX&#95;081550, and five members of the PvRBP family), as well as the leading pre-erythrocytic vaccine candidate (CSP) were measured using a multiplexed bead-based assay. Associations between IgG levels and infection status, age, and spatial location were explored., Results: Individuals from a low-transmission region of western Thailand reacted to all 11 P. vivax recombinant proteins. Significantly greater IgG levels were observed in the presence of a current P. vivax infection, despite all infected individuals being asymptomatic. IgG levels were also higher in adults (18 years and older) than in children. For most of the proteins, higher IgG levels were observed in individuals living closer to the Myanmar border and further away from local health services., Conclusions: Robust IgG responses were observed to most proteins and IgG levels correlated with surrogates of exposure, suggesting these antigens may serve as potential biomarkers of exposure, immunity, or both.

For reasons that remain unknown, the Plasmodium falciparum genome has an exceptionally high AT content compared to other Plasmodium species and eukaryotes in general - nearly 80% in coding regions and approaching 90% in non-coding regions. Here, we examine how this phenomenon relates to genome-wide patterns of de novo mutation. Mutation accumulation experiments were performed by sequential cloning of six P. falciparum isolates growing in human erythrocytes in vitro for 4 years, with 279 clones sampled for whole genome sequencing at different time points. Genome sequence analysis of these samples revealed a significant excess of G:C to A:T transitions compared to other types of nucleotide substitution, which would naturally cause AT content to equilibrate close to the level seen across the P. falciparum reference genome (80.6% AT). These data also uncover an extremely high rate of small indel mutation relative to other species, primarily associated with repetitive AT-rich sequences, in addition to larger-scale structural rearrangements focused in antigen-coding var genes. In conclusion, high AT content in P. falciparum is driven by a systematic mutational bias and ultimately leads to an unusual level of microstructural plasticity, raising the question of whether this contributes to adaptive evolution., (© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Invasion of erythrocytes by Plasmodium falciparum merozoites is necessary for malaria pathogenesis and is therefore a primary target for vaccine development. RH5 is a leading subunit vaccine candidate because anti-RH5 antibodies inhibit parasite growth and the interaction with its erythrocyte receptor basigin is essential for invasion. RH5 is secreted, complexes with other parasite proteins including CyRPA and RIPR, and contains a conserved N-terminal region (RH5Nt) of unknown function that is cleaved from the native protein. Here, we identify P113 as a merozoite surface protein that directly interacts with RH5Nt. Using recombinant proteins and a sensitive protein interaction assay, we establish the binding interdependencies of all the other known RH5 complex components and conclude that the RH5Nt-P113 interaction provides a releasable mechanism for anchoring RH5 to the merozoite surface. We exploit these findings to design a chemically synthesized peptide corresponding to RH5Nt, which could contribute to a cost-effective malaria vaccine.

Plasmodium vivax parasites cause the majority of malaria cases outside Africa, and are increasingly being acknowledged as a cause of severe disease. The unique attributes of P. vivax biology, particularly the capacity of the dormant liver stage, the hypnozoite, to maintain blood-stage infections even in the absence of active transmission, make blood-stage vaccines particularly attractive for this species. However, P. vivax vaccine development remains resolutely in first gear, with only a single blood-stage candidate having been evaluated in any depth. Experience with Plasmodium falciparum suggests that a much broader search for new candidates and a deeper understanding of high priority targets will be required to make significant advances. This review discusses some of the particular challenges of P. vivax blood-stage vaccine development, highlighting both recent advances and key remaining barriers to overcome in order to move development forward., (Copyright © 2016. Published by Elsevier Ltd.)

A family of apicomplexa-specific proteins containing AP2 DNA-binding domains (ApiAP2s) was identified in malaria parasites. This family includes sequence-specific transcription factors that are key regulators of development. However, functions for the majority of ApiAP2 genes remain unknown. Here, a systematic knockout screen in Plasmodium berghei identified ten ApiAP2 genes that were essential for mosquito transmission: four were critical for the formation of infectious ookinetes, and three were required for sporogony. We describe non-essential functions for AP2-O and AP2-SP proteins in blood stages, and identify AP2-G2 as a repressor active in both asexual and sexual stages. Comparative transcriptomics across mutants and developmental stages revealed clusters of co-regulated genes with shared cis promoter elements, whose expression can be controlled positively or negatively by different ApiAP2 factors. We propose that stage-specific interactions between ApiAP2 proteins on partly overlapping sets of target genes generate the complex transcriptional network that controls the Plasmodium life cycle., (Copyright © 2017 The Author(s). Published by Elsevier Inc. All rights reserved.)

As the Plasmodium parasite transitions between mammalian and mosquito host, it has to adjust quickly to new environments. Palmitoylation, a reversible and dynamic lipid post-translational modification, plays a central role in regulating this process and has been implicated with functions for parasite morphology, motility and host cell invasion. While proteins associated with the gliding motility machinery have been described to be palmitoylated, no palmitoyl transferase responsible for regulating gliding motility has previously been identified. Here, we characterize two palmityol transferases with gene tagging and gene deletion approaches. We identify DHHC3, a palmitoyl transferase, as a mediator of ookinete development, with a crucial role for gliding motility in ookinetes and sporozoites, and we co-localize the protein with a marker for the inner membrane complex in the ookinete stage. Ookinetes and sporozoites lacking DHHC3 are impaired in gliding motility and exhibit a strong phenotype in vivo; with ookinetes being significantly less infectious to their mosquito host and sporozoites being non-infectious to mice. Importantly, genetic complementation of the DHHC3-ko parasite completely restored virulence. We generated parasites lacking both DHHC3, as well as the palmitoyl transferase DHHC9, and found an enhanced phenotype for these double knockout parasites, allowing insights into the functional overlap and compensational nature of the large family of PbDHHCs. These findings contribute to our understanding of the organization and mechanism of the gliding motility machinery, which as is becoming increasingly clear, is mediated by palmitoylation., (© 2016 John Wiley & Sons Ltd.)

Palmitoylation is the post-translational reversible addition of the acyl moiety, palmitate, to cysteine residues of proteins and is involved in regulating protein trafficking, localization, stability and function. The Aspartate-Histidine-Histidine-Cysteine (DHHC) protein family, named for their highly conserved DHHC signature motif, is thought to be responsible for catalysing protein palmitoylation. Palmitoylation is widespread in all eukaryotes, including the malaria parasite, Plasmodium falciparum, where over 400 palmitoylated proteins are present in the asexual intraerythrocytic schizont stage parasites, including proteins involved in key aspects of parasite maturation and development. The P. falciparum genome includes 12 proteins containing the conserved DHHC motif. In this study, we adapted a palmitoyl-transferase activity assay for use with P. falciparum proteins and demonstrated for the first time that P. falciparum DHHC proteins are responsible for the palmitoylation of P. falciparum substrates. This assay also reveals that multiple DHHCs are capable of palmitoylating the same substrate, indicating functional redundancy at least in vitro. To test whether functional redundancy also exists in vivo, we investigated the endogenous localization and essentiality of a subset of schizont-expressed PfDHHC proteins. Individual PfDHHC proteins localized to distinct organelles, including parasite-specific organelles such as the rhoptries and inner membrane complex. Knock-out studies identified individual DHHCs that may be essential for blood-stage growth and others that were functionally redundant in the blood stages but may have functions in other stages of parasite development. Supporting this hypothesis, disruption of PfDHHC9 had no effect on blood-stage growth but reduced the formation of gametocytes, suggesting that this protein could be exploited as a transmission-blocking target. The localization and stage-specific expression of the DHHC proteins may be important for regulating their substrate specificity and thus may provide a path for inhibitor development., (© 2016 The Authors Cellular Microbiology Published by John Wiley & Sons Ltd.)

Background: Plasmodium vivax causes the majority of malaria episodes outside Africa, but remains a relatively understudied pathogen. The pathology of P. vivax infection depends critically on the parasite's ability to recognize and invade human erythrocytes. This invasion process involves an interaction between P. vivax Duffy Binding Protein (PvDBP) in merozoites and the Duffy antigen receptor for chemokines (DARC) on the erythrocyte surface. Whole-genome sequencing of clinical isolates recently established that some P. vivax genomes contain two copies of the PvDBP gene. The frequency of this duplication is particularly high in Madagascar, where there is also evidence for P. vivax infection in DARC-negative individuals. The functional significance and global prevalence of this duplication, and whether there are other copy number variations at the PvDBP locus, is unknown., Methodology/principal Findings: Using whole-genome sequencing and PCR to study the PvDBP locus in P. vivax clinical isolates, we found that PvDBP duplication is widespread in Cambodia. The boundaries of the Cambodian PvDBP duplication differ from those previously identified in Madagascar, meaning that current molecular assays were unable to detect it. The Cambodian PvDBP duplication did not associate with parasite density or DARC genotype, and ranged in prevalence from 20% to 38% over four annual transmission seasons in Cambodia. This duplication was also present in P. vivax isolates from Brazil and Ethiopia, but not India., Conclusions/significance: PvDBP duplications are much more widespread and complex than previously thought, and at least two distinct duplications are circulating globally. The same duplication boundaries were identified in parasites from three continents, and were found at high prevalence in human populations where DARC-negativity is essentially absent. It is therefore unlikely that PvDBP duplication is associated with infection of DARC-negative individuals, but functional tests will be required to confirm this hypothesis., Competing Interests: The authors have declared that no competing interests exist.

The widespread distribution and relapsing nature of Plasmodium vivax infection present major challenges for the elimination of malaria. To characterize the genetic diversity of this parasite in individual infections and across the population, we performed deep genome sequencing of >200 clinical samples collected across the Asia-Pacific region and analyzed data on >300,000 SNPs and nine regions of the genome with large copy number variations. Individual infections showed complex patterns of genetic structure, with variation not only in the number of dominant clones but also in their level of relatedness and inbreeding. At the population level, we observed strong signals of recent evolutionary selection both in known drug resistance genes and at new loci, and these varied markedly between geographical locations. These findings demonstrate a dynamic landscape of local evolutionary adaptation in the parasite population and provide a foundation for genomic surveillance to guide effective strategies for control and elimination of P. vivax., Competing Interests: The authors declare no competing financial interests.

Plasmodium falciparum, the major cause of malaria morbidity and mortality worldwide, is only distantly related to other human malaria parasites and has thus been placed in a separate subgenus, termed Laverania Parasites morphologically similar to P. falciparum have been identified in African apes, but only one other Laverania species, Plasmodium reichenowi from chimpanzees, has been formally described. Although recent studies have pointed to the existence of additional Laverania species, their precise number and host associations remain uncertain, primarily because of limited sampling and a paucity of parasite sequences other than from mitochondrial DNA. To address this, we used limiting dilution polymerase chain reaction to amplify additional parasite sequences from a large number of chimpanzee and gorilla blood and fecal samples collected at two sanctuaries and 30 field sites across equatorial Africa. Phylogenetic analyses of more than 2,000 new sequences derived from the mitochondrial, nuclear, and apicoplast genomes revealed six divergent and well-supported clades within the Laverania parasite group. Although two of these clades exhibited deep subdivisions in phylogenies estimated from organelle gene sequences, these sublineages were geographically defined and not present in trees from four unlinked nuclear loci. This greatly expanded sequence data set thus confirms six, and not seven or more, ape Laverania species, of which P. reichenowi, Plasmodium gaboni, and Plasmodium billcollinsi only infect chimpanzees, whereas Plasmodium praefalciparum, Plasmodium adleri, and Pladmodium blacklocki only infect gorillas. The new sequence data also confirm the P. praefalciparum origin of human P. falciparum., (© The Author(s) 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Genetic screening using random transposon insertions has been a powerful tool for uncovering biology in prokaryotes, where whole-genome saturating screens have been performed in multiple organisms. In eukaryotes, such screens have proven more problematic, in part because of the lack of a sensitive and robust system for identifying transposon insertion sites. We here describe quantitative insertion-site sequencing, or QIseq, which uses custom library preparation and Illumina sequencing technology and is able to identify insertion sites from both the 5' and 3' ends of the transposon, providing an inbuilt level of validation. The approach was developed using piggyBac mutants in the human malaria parasite Plasmodium falciparum but should be applicable to many other eukaryotic genomes. QIseq proved accurate, confirming known sites in >100 mutants, and sensitive, identifying and monitoring sites over a >10,000-fold dynamic range of sequence counts. Applying QIseq to uncloned parasites shortly after transfections revealed multiple insertions in mixed populations and suggests that >4000 independent mutants could be generated from relatively modest scales of transfection, providing a clear pathway to genome-scale screens in P. falciparum QIseq was also used to monitor the growth of pools of previously cloned mutants and reproducibly differentiated between deleterious and neutral mutations in competitive growth. Among the mutants with fitness defects was a mutant with a piggyBac insertion immediately upstream of the kelch protein K13 gene associated with artemisinin resistance, implying mutants in this gene may have competitive fitness costs. QIseq has the potential to enable the scale-up of piggyBac-mediated genetics across multiple eukaryotic systems., (© 2016 Bronner et al.; Published by Cold Spring Harbor Laboratory Press.)

Diversity at pathogen genetic loci can be driven by host adaptive immune selection pressure and may reveal proteins important for parasite biology. Population-based genome sequencing of Plasmodium falciparum, the parasite responsible for the most severe form of malaria, has highlighted two related polymorphic genes called dblmsp and dblmsp2, which encode Duffy binding-like (DBL) domain-containing proteins located on the merozoite surface but whose function remains unknown. Using recombinant proteins and transgenic parasites, we show that DBLMSP and DBLMSP2 directly and avidly bind human IgM via their DBL domains. We used whole genome sequence data from over 400 African and Asian P. falciparum isolates to show that dblmsp and dblmsp2 exhibit extreme protein polymorphism in their DBL domain, with multiple variants of two major allelic classes present in every population tested. Despite this variability, the IgM binding function was retained across diverse sequence representatives. Although this interaction did not seem to have an effect on the ability of the parasite to invade red blood cells, binding of DBLMSP and DBLMSP2 to IgM inhibited the overall immunoreactivity of these proteins to IgG from patients who had been exposed to the parasite. This suggests that IgM binding might mask these proteins from the host humoral immune system., (© 2016 by The American Society for Biochemistry and Molecular Biology, Inc.)

The clinical complications of malaria are caused by the parasite expansion in the blood. Invasion of erythrocytes is a complex process that depends on multiple receptor-ligand interactions. Identification of host receptors is paramount for fighting the disease as it could reveal new intervention targets, but the enucleated nature of erythrocytes makes genetic approaches impossible and many receptors remain unknown. Host-parasite interactions evolve rapidly and are therefore likely to be species-specific. As a results, understanding of invasion receptors outside the major human pathogen Plasmodium falciparum is very limited. Here we use mouse embryonic stem cells (mESCs) that can be genetically engineered and differentiated into erythrocytes to identify receptors for the rodent malaria parasite Plasmodium berghei. Two proteins previously implicated in human malaria infection: glycophorin C (GYPC) and Band-3 (Slc4a1) were deleted in mESCs to generate stable cell lines, which were differentiated towards erythropoiesis. In vitro infection assays revealed that while deletion of Band-3 has no effect, absence of GYPC results in a dramatic decrease in invasion, demonstrating the crucial role of this protein for P. berghei infection. This stem cell approach offers the possibility of targeting genes that may be essential and therefore difficult to disrupt in whole organisms and has the potential to be applied to a variety of parasites in diverse host cell types.

Background: Elimination of Plasmodium vivax malaria would be greatly facilitated by the development of an effective vaccine. A comprehensive and systematic characterization of antibodies to P. vivax antigens in exposed populations is useful in guiding rational vaccine design., Methodology/principal Findings: In this study, we investigated antibodies to a large library of P. vivax entire ectodomain merozoite proteins in 2 Asia-Pacific populations, analysing the relationship of antibody levels with markers of current and cumulative malaria exposure, and socioeconomic and clinical indicators. 29 antigenic targets of natural immunity were identified. Of these, 12 highly-immunogenic proteins were strongly associated with age and thus cumulative lifetime exposure in Solomon Islanders (P<0.001-0.027). A subset of 6 proteins, selected on the basis of immunogenicity and expression levels, were used to examine antibody levels in plasma samples from a population of young Papua New Guinean children with well-characterized individual differences in exposure. This analysis identified a strong association between reduced risk of clinical disease and antibody levels to P12, P41, and a novel hypothetical protein that has not previously been studied, PVX&#95;081550 (IRR 0.46-0.74; P<0.001-0.041)., Conclusion/significance: These data emphasize the benefits of an unbiased screening approach in identifying novel vaccine candidate antigens. Functional studies are now required to establish whether PVX&#95;081550 is a key component of the naturally-acquired protective immune response, a biomarker of immune status, or both.