Search

Background: The long co-evolution of Homo sapiens and Plasmodium falciparum has resulted in the selection of numerous human genetic variants that confer an advantage against severe malaria and death. One such variant is the Dantu blood group antigen, which is associated with 74% protection against severe and complicated P. falciparum malaria infections in homozygous individuals, similar to that provided by the sickle haemoglobin allele (HbS). Recent in vitro studies suggest that Dantu exerts this protection by increasing the surface tension of red blood cells, thereby impeding the ability of P. falciparum merozoites to invade them and reducing parasite multiplication. However, no studies have yet explored this hypothesis in vivo. Methods: We investigated the effect of Dantu on early phase P. falciparum (Pf) infections in a controlled human malaria infection (CHMI) study. 141 sickle-negative Kenyan adults were inoculated with 3.2 × 103 aseptic, purified, cryopreserved Pf sporozoites (PfSPZ Challenge) then monitored for blood-stage parasitaemia for 21 days by quantitative polymerase chain reaction (qPCR)analysis of the 18S ribosomal RNA P. falciparum gene. The primary endpoint was blood-stage P. falciparum parasitaemia of ≥500/μl while the secondary endpoint was the receipt of antimalarial treatment in the presence of parasitaemia of any density. On study completion, all participants were genotyped both for Dantu and for four other polymorphisms that are associated with protection against severe falciparum malaria: α+-thalassaemia, blood group O, G6PD deficiency, and the rs4951074 allele in the red cell calcium transporter ATP2B4. Results: The primary endpoint was reached in 25/111 (22.5%) non-Dantu subjects in comparison to 0/27 (0%) Dantu heterozygotes and 0/3 (0.0%) Dantu homozygotes (p=0.01). Similarly, 49/111 (44.1%) non-Dantu subjects reached the secondary endpoint in comparison to only 7/27 (25.9%) and 0/3 (0.0%) Dantu heterozygotes and homozygotes, respectively (p=0.021). No significant impacts on either outcome were seen for any of the other genetic variants under study. Conclusions: This study reveals, for the first time, that the Dantu blood group is associated with high-level protection against early, non-clinical, P. falciparum malaria infections in vivo. Learning more about the mechanisms involved could potentially lead to new approaches to the prevention or treatment of the disease. Our study illustrates the power of CHMI with PfSPZ Challenge for directly testing the protective impact of genotypes previously identified using other methods. Funding: The Kenya CHMI study was supported by an award from Wellcome (grant number 107499). SK was supported by a Training Fellowship (216444/Z/19/Z), TNW by a Senior Research Fellowship (202800/Z/16/Z), JCR by an Investigator Award (220266/Z/20/Z), and core support to the KEMRI-Wellcome Trust Research Programme in Kilifi, Kenya (203077), all from Wellcome. The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication. For the purpose of Open Access, the authors have applied a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission. Clinical trial number: NCT02739763

We describe the MalariaGEN Pf7 data resource, the seventh release of Plasmodium falciparum genome variation data from the MalariaGEN network. It comprises over 20,000 samples from 82 partner studies in 33 countries, including several malaria endemic regions that were previously underrepresented. For the first time we include dried blood spot samples that were sequenced after selective whole genome amplification, necessitating new methods to genotype copy number variations. We identify a large number of newly emerging crt mutations in parts of Southeast Asia, and show examples of heterogeneities in patterns of drug resistance within Africa and within the Indian subcontinent. We describe the profile of variations in the C-terminal of the csp gene and relate this to the sequence used in the RTS,S and R21 malaria vaccines. Pf7 provides high-quality data on genotype calls for 6 million SNPs and short indels, analysis of large deletions that cause failure of rapid diagnostic tests, and systematic characterisation of six major drug resistance loci, all of which can be freely downloaded from the MalariaGEN website.

The human malaria parasite Plasmodium falciparum is globally widespread, but its prevalence varies significantly between and even within countries. Most population genetic studies in P. falciparum focus on regions of high transmission where parasite populations are large and genetically diverse, such as sub-Saharan Africa. Understanding population dynamics in low transmission settings, however, is of particular importance as these are often where drug resistance first evolves. Here, we use the Pacific Coast of Colombia and Ecuador as a model for understanding the population structure and evolution of Plasmodium parasites in small populations harboring less genetic diversity. The combination of low transmission and a high proportion of monoclonal infections means there are few outcrossing events and clonal lineages persist for long periods of time. Yet despite this, the population is evolutionarily labile and has successfully adapted to changes in drug regime. Using newly sequenced whole genomes, we measure relatedness between 166 parasites, calculated as identity by descent (IBD), and find 17 distinct but highly related clonal lineages, six of which have persisted in the region for at least a decade. This inbred population structure is captured in more detail with IBD than with other common population structure analyses like PCA, ADMIXTURE, and distance-based trees. We additionally use patterns of intra-chromosomal IBD and an analysis of haplotypic variation to explore past selection events in the region. Two genes associated with chloroquine resistance, crt and aat1, show evidence of hard selective sweeps, while selection appears soft and/or incomplete at three other key resistance loci (dhps, mdr1, and dhfr). Overall, this work highlights the strength of IBD analyses for studying parasite population structure and resistance evolution in regions of low transmission, and emphasizes that drug resistance can evolve and spread in small populations, as will occur in any region nearing malaria elimination.

BackgroundTo make progress towards malaria elimination, a highly effective vaccine targeting Plasmodium vivax is urgently needed. Evaluating the kinetics of natural antibody responses to vaccine candidate antigens after acute vivax malaria can inform the design of serological markers of exposure and vaccines.Methodology/principal findingsThe responses of IgG antibodies to 9 P. vivax vaccine candidate antigens were evaluated in longitudinal serum samples from Brazilian individuals collected at the time of acute vivax malaria and 30, 60, and 180 days afterwards. Antigen-specific IgG correlations, seroprevalence, and half-lives were determined for each antigen using the longitudinal data. Antibody reactivities against Pv41 and PVX_081550 strongly correlated with each other at each of the four time points. The analysis identified robust responses in terms of magnitude and seroprevalence against Pv41 and PvGAMA at 30 and 60 days. Among the 8 P. vivax antigens demonstrating >50% seropositivity across all individuals, antibodies specific to PVX_081550 had the longest half-life (100 days; 95% CI, 83-130 days), followed by PvRBP2b (91 days; 95% CI, 76-110 days) and Pv12 (82 days; 95% CI, 64-110 days).Conclusion/significanceThis study provides an in-depth assessment of the kinetics of antibody responses to key vaccine candidate antigens in Brazilians with acute vivax malaria. Follow-up studies are needed to determine whether the longer-lived antibody responses induced by natural infection are effective in controlling blood-stage infection and mediating clinical protection.

This report describes the MalariaGEN Pv4 dataset, a new release of curated genome variation data on 1,895 samples of Plasmodium vivax collected at 88 worldwide locations between 2001 and 2017. It includes 1,370 new samples contributed by MalariaGEN and VivaxGEN partner studies in addition to previously published samples from these and other sources. We provide genotype calls at over 4.5 million variable positions including over 3 million single nucleotide polymorphisms (SNPs), as well as short indels and tandem duplications. This enlarged dataset highlights major compartments of parasite population structure, with clear differentiation between Africa, Latin America, Oceania, Western Asia and different parts of Southeast Asia. Each sample has been classified for drug resistance to sulfadoxine, pyrimethamine and mefloquine based on known markers at the dhfr, dhps and mdr1 loci. The prevalence of all of these resistance markers was much higher in Southeast Asia and Oceania than elsewhere. This open resource of analysis-ready genome variation data from the MalariaGEN and VivaxGEN networks is driven by our collective goal to advance research into the complex biology of P. vivax and to accelerate genomic surveillance for malaria control and elimination.

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.

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.

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 PlasmoGEM 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.

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 PlasmoGEM 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.

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.

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.

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.

Abstract Challenges in classifying recurrent Plasmodium vivax infections constrain surveillance of antimalarial efficacy and transmission. Recurrent infections may arise from activation of dormant liver stages (relapse), blood-stage treatment failure (recrudescence) or reinfection. Molecular inference of familial relatedness (identity-by-descent or IBD) can help resolve the probable origin of recurrences. As whole genome sequencing of P. vivax remains challenging, targeted genotyping methods are needed for scalability. We describe a P. vivax marker discovery framework to identify and select panels of microhaplotypes (multi-allelic markers within small, amplifiable segments of the genome) that can accurately capture IBD. We evaluate panels of 50–250 microhaplotypes discovered in a global set of 615 P. vivax genomes. A candidate global 100-microhaplotype panel exhibits high marker diversity in the Asia-Pacific, Latin America and horn of Africa (median H E = 0.70–0.81) and identifies 89% of the polyclonal infections detected with genome-wide datasets. Data simulations reveal lower error in estimating pairwise IBD using microhaplotypes relative to traditional biallelic SNP barcodes. The candidate global panel also exhibits high accuracy in predicting geographic origin and captures local infection outbreak and bottlenecking events. Our framework is open-source enabling customised microhaplotype discovery and selection, with potential for porting to other species or data resources.

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.

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, p90%), EBP, DBPII, RBP1a, CyRPA, and PVX_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.

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.

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.

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

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:A vaccine targeting Plasmodium vivax will be an essential component of any comprehensive malaria elimination program, but major gaps in our understanding of P. vivax biology, including the protein-protein interactions that mediate merozoite invasion of reticulocytes, hinder the search for candidate antigens. Only one ligand-receptor interaction has been identified, that between P. vivax Duffy Binding Protein (PvDBP) and the erythrocyte Duffy Antigen Receptor for Chemokines (DARC), and strain-specific immune responses to PvDBP make it a complex vaccine target. To broaden the repertoire of potential P. vivax merozoite-stage vaccine targets, we exploited a recent breakthrough in expressing full-length ectodomains of Plasmodium proteins in a functionally-active form in mammalian cells and initiated a large-scale study of P. vivax merozoite proteins that are potentially involved in reticulocyte binding and invasion. METHODOLOGY/PRINCIPAL FINDINGS:We selected 39 P. vivax proteins that are predicted to localize to the merozoite surface or invasive secretory organelles, some of which show homology to P. falciparum vaccine candidates. Of these, we were able to express 37 full-length protein ectodomains in a mammalian expression system, which has been previously used to express P. falciparum invasion ligands such as PfRH5. To establish whether the expressed proteins were correctly folded, we assessed whether they were recognized by antibodies from Cambodian patients with acute vivax malaria. IgG from these samples showed at least a two-fold change in reactivity over naïve controls in 27 of 34 antigens tested, and the majority showed heat-labile IgG immunoreactivity, suggesting the presence of conformation-sensitive epitopes and native tertiary protein structures. Using a method specifically designed to detect low-affinity, extracellular protein-protein interactions, we confirmed a predicted interaction between P. vivax 6-cysteine proteins P12 and P41, further suggesting that the proteins are natively folded and functional. This screen also identified two novel protein-protein interactions, between P12 and PVX_110945, and between MSP3.10 and MSP7.1, the latter of which was confirmed by surface plasmon resonance. CONCLUSIONS/SIGNIFICANCE:We produced a new library of recombinant full-length P. vivax ectodomains, established that the majority of them contain tertiary structure, and used them to identify predicted and novel protein-protein interactions. As well as identifying new interactions for further biological studies, this library will be useful in identifying P. vivax proteins with vaccine potential, and studying P. vivax malaria pathogenesis and immunity. TRIAL REGISTRATION:ClinicalTrials.gov NCT00663546.

During blood stage Plasmodium falciparum infection, merozoites invade uninfected erythrocytes via a complex, multistep process involving a series of distinct receptor-ligand binding events. Understanding each element in this process increases the potential to block the parasite's life cycle via drugs or vaccines. To investigate specific receptor-ligand interactions, they were systematically blocked using a combination of genetic deletion, enzymatic receptor cleavage and inhibition of binding via antibodies, peptides and small molecules, and the resulting temporal changes in invasion and morphological effects on erythrocytes were filmed using live cell imaging. Analysis of the videos have shown receptor-ligand interactions occur in the following sequence with the following cellular morphologies; 1) an early heparin-blockable interaction which weakly deforms the erythrocyte, 2) EBA and PfRh ligands which strongly deform the erythrocyte, a process dependant on the merozoite's actin-myosin motor, 3) a PfRh5-basigin binding step which results in a pore or opening between parasite and host through which it appears small molecules and possibly invasion components can flow and 4) an AMA1-RON2 interaction that mediates tight junction formation, which acts as an anchor point for internalization. In addition to enhancing general knowledge of apicomplexan biology, this work provides a rational basis to combine sequentially acting merozoite vaccine candidates in a single multi-receptor-blocking vaccine.

Plasmodium knowlesi is a newly described zoonosis that causes malaria in the human population that can be severe and fatal. The study of P. knowlesi parasites from human clinical isolates is relatively new and, in order to obtain maximum information from patient sample collections, we explored the possibility of generating P. knowlesi genome sequences from archived clinical isolates. Our patient sample collection consisted of frozen whole blood samples that contained excessive human DNA contamination and, in that form, were not suitable for parasite genome sequencing. We developed a method to reduce the amount of human DNA in the thawed blood samples in preparation for high throughput parasite genome sequencing using Illumina HiSeq and MiSeq sequencing platforms. Seven of fifteen samples processed had sufficiently pure P. knowlesi DNA for whole genome sequencing. The reads were mapped to the P. knowlesi H strain reference genome and an average mapping of 90% was obtained. Genes with low coverage were removed leaving 4623 genes for subsequent analyses. Previously we identified a DNA sequence dimorphism on a small fragment of the P. knowlesi normocyte binding protein xa gene on chromosome 14. We used the genome data to assemble full-length Pknbpxa sequences and discovered that the dimorphism extended along the gene. An in-house algorithm was developed to detect SNP sites co-associating with the dimorphism. More than half of the P. knowlesi genome was dimorphic, involving genes on all chromosomes and suggesting that two distinct types of P. knowlesi infect the human population in Sarawak, Malaysian Borneo. We use P. knowlesi clinical samples to demonstrate that Plasmodium DNA from archived patient samples can produce high quality genome data. We show that analyses, of even small numbers of difficult clinical malaria isolates, can generate comprehensive genomic information that will improve our understanding of malaria parasite diversity and pathobiology.

The most polymorphic gene family in P. falciparum is the ∼60 var genes distributed across parasite chromosomes, both in the subtelomeres and in internal regions. They encode hypervariable surface proteins known as P. falciparum erythrocyte membrane protein 1 (PfEMP1) that are critical for pathogenesis and immune evasion in Plasmodium falciparum. How var gene sequence diversity is generated is not currently completely understood. To address this, we constructed large clone trees and performed whole genome sequence analysis to study the generation of novel var gene sequences in asexually replicating parasites. While single nucleotide polymorphisms (SNPs) were scattered across the genome, structural variants (deletions, duplications, translocations) were focused in and around var genes, with considerable variation in frequency between strains. Analysis of more than 100 recombination events involving var exon 1 revealed that the average nucleotide sequence identity of two recombining exons was only 63% (range: 52.7-72.4%) yet the crossovers were error-free and occurred in such a way that the resulting sequence was in frame and domain architecture was preserved. Var exon 1, which encodes the immunologically exposed part of the protein, recombined in up to 0.2% of infected erythrocytes in vitro per life cycle. The high rate of var exon 1 recombination indicates that millions of new antigenic structures could potentially be generated each day in a single infected individual. We propose a model whereby var gene sequence polymorphism is mainly generated during the asexual part of the life cycle.

Many critical events in the Plasmodium life cycle rely on the controlled release of Ca²⁺ from intracellular stores to activate stage-specific Ca²⁺-dependent protein kinases. Using the motility of Plasmodium berghei ookinetes as a signalling paradigm, we show that the cyclic guanosine monophosphate (cGMP)-dependent protein kinase, PKG, maintains the elevated level of cytosolic Ca²⁺ required for gliding motility. We find that the same PKG-dependent pathway operates upstream of the Ca²⁺ signals that mediate activation of P. berghei gametocytes in the mosquito and egress of Plasmodium falciparum merozoites from infected human erythrocytes. Perturbations of PKG signalling in gliding ookinetes have a marked impact on the phosphoproteome, with a significant enrichment of in vivo regulated sites in multiple pathways including vesicular trafficking and phosphoinositide metabolism. A global analysis of cellular phospholipids demonstrates that in gliding ookinetes PKG controls phosphoinositide biosynthesis, possibly through the subcellular localisation or activity of lipid kinases. Similarly, phosphoinositide metabolism links PKG to egress of P. falciparum merozoites, where inhibition of PKG blocks hydrolysis of phosphatidylinostitol (4,5)-bisphosphate. In the face of an increasing complexity of signalling through multiple Ca²⁺ effectors, PKG emerges as a unifying factor to control multiple cellular Ca²⁺ signals essential for malaria parasite development and transmission.

All the symptoms and pathology of malaria are caused by the intraerythrocytic stages of the Plasmodium parasite life cycle. Because Plasmodium parasites cannot replicate outside a host cell, their ability to recognize and invade erythrocytes is an essential step for both parasite survival and malaria pathogenesis. This makes invasion a conceptually attractive vaccine target, especially because it is one of the few stages when the parasite is directly exposed to the host humoral immune system. This apparent vulnerability, however, has been countered by the parasite, which has evolved sophisticated molecular mechanisms to evade the host immune response so that parasites asymptomatically replicate within immune individuals. These mechanisms include the expansion of parasite invasion ligands, resulting in multiple and apparently redundant invasion "pathways", highly polymorphic parasite surface proteins that are immunologically distinct, and parasite proteins which are poorly immunogenic. These formidable defences have so far thwarted attempts to develop an effective blood-stage vaccine, leading many to question whether there really is an exploitable chink in the parasite's immune evasion defences. Here, we review recent advances in the molecular understanding of the P. falciparum erythrocyte invasion field, discuss some of the challenges that have so far prevented the development of blood-stage vaccines, and conclude that the parasite invasion ligand RH5 represents an essential pinch point that might be vulnerable to vaccination.

The motility and invasion of Plasmodium parasites is believed to require a cytoplasmic actin-myosin motor associated with a cell surface ligand belonging to the TRAP (thrombospondin-related anonymous protein) family. Current models of invasion usually invoke the existence of specific receptors for the TRAP-family ligands on the surface of the host cell; however, the identities of these receptors remain largely unknown. Here, we identify the GPI-linked protein Semaphorin-7A (CD108) as an erythrocyte receptor for the P. falciparum merozoite-specific TRAP homolog (MTRAP) by using a systematic screening approach designed to detect extracellular protein interactions. The specificity of the interaction was demonstrated by showing that binding was saturable and by quantifying the equilibrium and kinetic biophysical binding parameters using surface plasmon resonance. We found that two MTRAP monomers interact via their tandem TSR domains with the Sema domains of a Semaphorin-7A homodimer. Known naturally-occurring polymorphisms in Semaphorin-7A did not quantitatively affect MTRAP binding nor did the presence of glycans on the receptor. Attempts to block the interaction during in vitro erythrocyte invasion assays using recombinant proteins and antibodies showed no significant inhibitory effect, suggesting the inaccessibility of the complex to proteinaceous blocking agents. These findings now provide important experimental evidence to support the model that parasite TRAP-family ligands interact with specific host receptors during cellular invasion.

Alveolins, or inner membrane complex (IMC) proteins, are components of the subpellicular network that forms a structural part of the pellicle of malaria parasites. In Plasmodium berghei, deletions of three alveolins, IMC1a, b, and h, each resulted in reduced mechanical strength and gliding velocity of ookinetes or sporozoites. Using time lapse imaging, we show here that deletion of IMC1h (PBANKA_143660) also has an impact on the directionality and motility behaviour of both ookinetes and sporozoites. Despite their marked motility defects, sporozoites lacking IMC1h were able to invade mosquito salivary glands, allowing us to investigate the role of IMC1h in colonisation of the mammalian host. We show that IMC1h is essential for sporozoites to progress through the dermis in vivo but does not play a significant role in hepatoma cell transmigration and invasion in vitro. Colocalisation of IMC1h with the residual IMC in liver stages was detected up to 30 hours after infection and parasites lacking IMC1h showed developmental defects in vitro and a delayed onset of blood stage infection in vivo. Together, these results suggest that IMC1h is involved in maintaining the cellular architecture which supports normal motility behaviour, access of the sporozoites to the blood stream, and further colonisation of the mammalian host.

The genomes of Plasmodium parasites that cause malaria in humans, other primates, birds, and rodents all encode multiple 6-cys proteins. Distinct 6-cys protein family members reside on the surface at each extracellular life cycle stage and those on the surface of liver infective and sexual stages have been shown to play important roles in hepatocyte growth and fertilization respectively. However, 6-cys proteins associated with the blood-stage forms of the parasite have no known function. Here we investigate the biochemical nature and function of two blood-stage 6-cys proteins in Plasmodium falciparum, the most pathogenic species to afflict humans. We show that native P12 and P41 form a stable heterodimer on the infective merozoite surface and are secreted following invasion, but could find no evidence that this complex mediates erythrocyte-receptor binding. That P12 and P41 do not appear to have a major role as adhesins to erythrocyte receptors was supported by the observation that antisera to these proteins did not substantially inhibit erythrocyte invasion. To investigate other functional roles for these proteins their genes were successfully disrupted in P. falciparum, however P12 and P41 knockout parasites grew at normal rates in vitro and displayed no other obvious phenotypic changes. It now appears likely that these blood-stage 6-cys proteins operate as a pair and play redundant roles either in erythrocyte invasion or in host-immune interactions.

Forces and mechanical properties of cells and tissues set constraints on biological functions, and are key determinants of human physiology. Changes in cell mechanics may arise from disease, or directly contribute to pathogenesis. Malaria gives many striking examples. Plasmodium parasites, the causative agents of malaria, are single-celled organisms that cannot survive outside their hosts; thus, thost-pathogen interactions are fundamental for parasite’s biological success and to the host response to infection. These interactions are often combinations of biochemical and mechanical factors, but most research focuses on the molecular side. However, Plasmodium infection of human red blood cells leads to changes in their mechanical properties, which has a crucial impact on disease pathogenesis because of the interaction of infected red blood cells with other human tissues through various adhesion mechanisms, which can be probed and modelled with biophysical techniques. Recently, natural polymorphisms affecting red blood cell biomechanics have also been shown to protect human populations, highlighting the potential of understanding biomechanical factors to inform future vaccines and drug development. Here we review biophysical techniques that have revealed new aspects of Plasmodium falciparum invasion of red blood cells and cytoadhesion of infected cells to the host vasculature. These mechanisms occur differently across Plasmodium species and are linked to malaria pathogenesis. We highlight promising techniques from the fields of bioengineering, immunomechanics, and soft matter physics that could be beneficial for studying malaria. Some approaches might also be applied to other phases of the malaria lifecycle and to apicomplexan infections with complex host-pathogen interactions.