Novel stem cell technologies are powerful tools to understand the impact of human factors on Plasmodium falciparum malaria

Plasmodium falciparum parasites have a complex life cycle, but the most clinically relevant stage of the disease is the invasion of erythrocytes and the proliferation of the parasite in the blood. The influence of human genetic traits on malaria has been known for a long time, however understanding the role of the proteins involved is hampered by the anuclear nature of erythrocytes that makes them inaccessible to genetic tools. Here we overcome this limitation using stem cells to generate erythroid cells with an in-vitro differentiation protocol and assess parasite invasion with an adaptation of flow cytometry to detect parasite hemozoin. We combine this strategy with reprogramming of patient cells to Induced Pluripotent Stem Cells and genome editing to understand the role of key genes and human traits in malaria infection. We show that deletion of basigin ablates invasion while deletion of ATP2B4 has a minor effect and that erythroid cells from reprogrammed patient-derived HbBart α-thalassemia samples poorly support infection. The possibility to obtain patient-secific and genetically modifed erythoid cells offers an unparalleled opportunity to study the role of human genes and polymorphisms in malaria allowing preservation of the genomic background to demonstrate their function and understand their mechanisms.

Genomic analysis of Plasmodium vivax describes patterns of connectivity and putative drivers of adaptation in Ethiopia

Ethiopia has the greatest burden of Plasmodium vivax in Africa, but little is known about the epidemiological landscape of parasites across the country. We analysed the genomic diversity of 137 P. vivax isolates collected nine Ethiopian districts from 2012 to 2016. Signatures of selection were detected by cross-country comparisons with isolates from Thailand (n = 104) and Indonesia (n = 111), representing regions with low and high chloroquine resistance respectively. 26% (35/137) of Ethiopian infections were polyclonal, and 48.5% (17/35) of these comprised highly related clones (within-host identity-by-descent > 25%), indicating frequent co-transmission and superinfection. Parasite gene flow between districts could not be explained entirely by geographic distance, with economic and cultural factors hypothesised to have an impact on connectivity. Amplification of the duffy binding protein gene (pvdbp1) was prevalent across all districts (16-75%). Cross-population haplotype homozygosity revealed positive selection in a region proximal to the putative chloroquine resistance transporter gene (pvcrt-o). An S25P variant in amino acid transporter 1 (pvaat1), whose homologue has recently been implicated in P. falciparum chloroquine resistance evolution, was prevalent in Ethiopia (96%) but not Thailand or Indonesia (35-53%). The genomic architecture in Ethiopia highlights circulating variants of potential public health concern in an endemic setting with evidence of stable transmission.

Malaria: moving beyond the search for magic bullets

Round table discussion on challenges and opportunities in malaria research with Elena Levashina, Dominique Soldati-Favre, Andrew Waters, Friedrich Frischknecht, and Julian Rayner.

A systematic analysis of the human immune response to Plasmodium vivax

BACKGROUNDThe biology of Plasmodium vivax is markedly different from that of P. falciparum; how this shapes the immune response to infection remains unclear. To address this shortfall, we inoculated human volunteers with a clonal field isolate of P. vivax and tracked their response through infection and convalescence.METHODSParticipants were injected intravenously with blood-stage parasites and infection dynamics were tracked in real time by quantitative PCR. Whole blood samples were used for high dimensional protein analysis, RNA sequencing, and cytometry by time of flight, and temporal changes in the host response to P. vivax were quantified by linear regression. Comparative analyses with P. falciparum were then undertaken using analogous data sets derived from prior controlled human malaria infection studies.RESULTSP. vivax rapidly induced a type I inflammatory response that coincided with hallmark features of clinical malaria. This acute-phase response shared remarkable overlap with that induced by P. falciparum but was significantly elevated (at RNA and protein levels), leading to an increased incidence of pyrexia. In contrast, T cell activation and terminal differentiation were significantly increased in volunteers infected with P. falciparum. Heterogeneous CD4+ T cells were found to dominate this adaptive response and phenotypic analysis revealed unexpected features normally associated with cytotoxicity and autoinflammatory disease.CONCLUSIONP. vivax triggers increased systemic interferon signaling (cf P. falciparum), which likely explains its reduced pyrogenic threshold. In contrast, P. falciparum drives T cell activation far in excess of P. vivax, which may partially explain why falciparum malaria more frequently causes severe disease.TRIAL REGISTRATIONClinicalTrials.gov NCT03797989.FUNDINGThe European Union's Horizon 2020 Research and Innovation programme, the Wellcome Trust, and the Royal Society.

Extracellular vesicles could be a putative posttranscriptional regulatory mechanism that shapes intracellular RNA levels in Plasmodium falciparum

Plasmodium falciparum secretes extracellular vesicles (PfEVs) that contain parasite-derived RNA. However, the significance of the secreted RNA remains unexplored. Here, we compare secreted and intracellular RNA from asexual cultures of six P. falciparum lines. We find that secretion of RNA via extracellular vesicles is not only periodic throughout the asexual intraerythrocytic developmental cycle but is also highly conserved across P. falciparum isolates. We further demonstrate that the phases of RNA secreted via extracellular vesicles are discernibly shifted compared to those of the intracellular RNA within the secreting whole parasite. Finally, transcripts of genes with no known function during the asexual intraerythrocytic developmental cycle are enriched in PfEVs compared to the whole parasite. We conclude that the secretion of extracellular vesicles could be a putative posttranscriptional RNA regulation mechanism that is part of or synergise the classic RNA decay processes to maintain intracellular RNA levels in P. falciparum.

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

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

The global transcriptome of Plasmodium falciparum mid-stage gametocytes (stages II-IV) appears largely conserved and gametocyte-specific gene expression patterns vary in clinical isolates

Our overall understanding of the developmental biology of malaria parasites has been greatly enhanced by recent advances in transcriptomic analysis. However, most of these investigations rely on laboratory strains (LS) that were adapted into in vitro culture many years ago, and the transcriptomes of clinical isolates (CI) circulating in human populations have not been assessed. In this study, RNA-seq was used to compare the global transcriptome of mid-stage gametocytes derived from three short-term cultured CI, with gametocytes derived from the NF54 reference laboratory strain. The core transcriptome appeared to be consistent between CI- and LS-derived gametocyte preparations, but some important differences were also observed. A majority of gametocyte-specific genes (43/53) appear to have relatively higher expression in CI-derived gametocytes than in LS-derived gametocytes, but a K-means clustering analysis showed that genes involved in flagellum- and microtubule-based processes (movement/motility) were more abundant in both groups, albeit with some differences between them. In addition, gametocytes from one CI described as CI group II gametocytes (CI:GGII) showed gene expression variation in the form of reduced gametocyte-specific gene expression compared to the other two CI-derived gametocytes (CI gametocyte group I, CI:GGI), although the mixed developmental stages used in our study is a potential confounder, only partially mitigated by the inclusion of multiple replicates for each CI. Overall, our study suggests that there may be subtle differences in the gene expression profiles of mid-stage gametocytes from CI relative to the NF54 reference strain of Plasmodium falciparum. Thus, it is necessary to deploy gametocyte-producing clinical parasite isolates to fully understand the diversity of gene expression strategies that may occur during the sequestered development of parasite sexual stages. IMPORTANCE Maturing gametocytes of Plasmodium falciparum are known to sequester away from peripheral circulation into the bone marrow until they are mature. Blocking gametocyte sequestration can prevent malaria transmission from humans to mosquitoes, but most studies aim to understand gametocyte development utilizing long-term adapted laboratory lines instead of clinical isolates. This is a particular issue for our understanding of the sexual stages, which are known to decrease rapidly during adaptation to long-term culture, meaning that many LS are unable to produce transmissible gametocytes. Using RNA-seq, we investigated the global transcriptome of mid-stage gametocytes derived from three clinical isolates and a reference strain (NF54). This identified important differences in gene expression profiles between immature gametocytes of CI and the NF54 reference strain of P. falciparum, suggesting increased investment in gametocytogenesis in clinical isolates. Our transcriptomic data highlight the use of clinical isolates in studying the morphological, cellular features and molecular biology of gametocytes.

Development of an Affordable ELISA Targeting the SARS-CoV-2 Nucleocapsid and Its Application to Samples from the Ongoing COVID-19 Epidemic in Ghana

INTRODUCTION

The true nature of the population spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in populations is often not fully known as most cases, particularly in Africa, are asymptomatic. Finding the true magnitude of SARS-CoV-2 spread is crucial to provide actionable data about the epidemiological progress of the disease for researchers and policymakers. This study developed and optimized an antibody enzyme-linked immunosorbent assay (ELISA) using recombinant nucleocapsid antigen expressed in-house using a simple bacterial expression system.

METHODS

Nucleocapsid protein from SARS-CoV-2 was expressed and purified from Escherichia coli. Plasma samples used for the assay development were obtained from Ghanaian SARS-CoV-2 seropositive individuals during the pandemic, while seronegative controls were plasma samples collected from blood donors before the coronavirus disease 2019 (COVID-19) pandemic. Another set of seronegative controls was collected during the COVID-19 pandemic. Antibody detection and levels within the samples were validated using commercial kits and Luminex. Analyses were performed using GraphPad Prism, and the sensitivity, specificity and background cut-off were calculated.

RESULTS AND DISCUSSION

This low-cost ELISA (£0.96/test) assay has a high prediction of 98.9%, and sensitivity and specificity of 97% and 99%, respectively. The assay was subsequently used to screen plasma from SARS-CoV-2 RT-PCR-positive Ghanaians. The assay showed no significant difference in nucleocapsid antibody levels between symptomatic and asymptomatic, with an increase of the levels over time. This is in line with our previous publication.

CONCLUSION

This study developed a low-cost and transferable assay that enables highly sensitive and specific detection of human anti-SARS-CoV-2 IgG antibodies. This assay can be modified to include additional antigens and used for continuous monitoring of sero-exposure to SARS-CoV-2 in West Africa.

Chemogenomic Profiling of a Plasmodium falciparum Transposon Mutant Library Reveals Shared Effects of Dihydroartemisinin and Bortezomib on Lipid Metabolism and Exported Proteins

The antimalarial activity of the frontline drug artemisinin involves generation of reactive oxygen species (ROS) leading to oxidative damage of parasite proteins. To achieve homeostasis and maintain protein quality control in the overwhelmed parasite, the ubiquitin-proteasome system kicks in. Even though molecular markers for artemisinin resistance like pfkelch13 have been identified, the intricate network of mechanisms driving resistance remains to be elucidated. Here, we report a forward genetic screening strategy that enables a broader identification of genetic factors responsible for altering sensitivity to dihydroartemisinin (DHA) and a proteasome inhibitor, bortezomib (BTZ). Using a library of isogenic piggyBac mutants in P. falciparum, we defined phenotype-genotype associations influencing drug responses and highlighted shared mechanisms between the two processes, which mainly included proteasome-mediated degradation and the lipid metabolism genes. Additional transcriptomic analysis of a DHA/BTZ-sensitive piggyBac mutant showed it is possible to find differences between the two response mechanisms on the specific components for regulation of the exportome. Our results provide further insight into the molecular mechanisms of antimalarial drug resistance. IMPORTANCE Malaria control is seriously threatened by the emergence and spread of Plasmodium falciparum resistance to the leading antimalarial, artemisinin. The potent killing activity of artemisinin results from oxidative damage unleashed by free heme activation released by hemoglobin digestion. Although the ubiquitin-proteasome system is considered critical for parasite survival of this toxicity, the diverse genetic changes linked to artemisinin resistance are complex and, so far, have not included the ubiquitin-proteasome system. In this study, we use a systematic forward genetic approach by screening a library of P. falciparum random piggyBac mutants to decipher the genetic factors driving malaria parasite responses to the oxidative stress caused by antimalarial drugs. This study compares phenotype-genotype associations influencing dihydroartemisinin responses with the proteasome inhibitor bortezomib to delineate the role of ubiquitin-proteasome system. Our study highlights shared and unique pathways from the complex array of molecular processes critical for P. falciparum survival resulting from the oxidative damage of artemisinin.

Phenotypic Screens Identify Genetic Factors Associated with Gametocyte Development in the Human Malaria Parasite Plasmodium falciparum

Transmission of the deadly malaria parasite Plasmodium falciparum from humans to mosquitoes is achieved by specialized intraerythrocytic sexual forms called gametocytes. Though the crucial regulatory mechanisms leading to gametocyte commitment have recently come to light, networks of genes that control sexual development remain to be elucidated. Here, we report a pooled-mutant screen to identify genes associated with gametocyte development in P. falciparum. Our results categorized genes that modulate gametocyte progression as hypoproducers or hyperproducers of gametocytes, and the in-depth analysis of individual clones confirmed phenotypes in sexual commitment rates and putative functions in gametocyte development. We present a new set of genes that have not been implicated in gametocytogenesis before and demonstrate the potential of forward genetic screens in isolating genes impacting parasite sexual biology, an exciting step toward the discovery of new antimalarials for a globally significant pathogen. IMPORTANCE Blocking human-to-vector transmission is an essential step toward malaria elimination. Gametocytes are solely responsible for achieving this transmission and represent an opportunity for therapeutic intervention. While these falciform-shaped parasite stages were first discovered in the 1880s, our understanding of the genetic determinants responsible for their formation and molecular mechanisms that drive their development is limited. In this work, we developed a scalable screening methodology with piggyBac mutants to identify genes that influence the development of gametocytes in the most lethal human malaria parasite, P. falciparum. By doing so, we lay the foundation for large-scale functional genomic studies specifically designed to address remaining questions about sexual commitment, maturation, and mosquito infection in P. falciparum. Such functional genetic screens will serve to expedite the identification of essential pathways and processes for the development of novel transmission-blocking agents.

The Dantu blood group prevents parasite growth in vivo: Evidence from a controlled human malaria infection study

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 × 10³ 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.

Lineage-informative microhaplotypes for spatio-temporal surveillance of Plasmodium vivax malaria parasites

Challenges in understanding the origin of recurrent Plasmodium vivax infections constrains the surveillance of antimalarial efficacy and transmission of this neglected parasite. Recurrent infections within an individual may arise from activation of dormant liver stages (relapse), blood-stage treatment failure (recrudescence) or new inoculations (reinfection). Molecular inference of familial relatedness (identity-by-descent or IBD) based on whole genome sequence data, together with analysis of the intervals between parasitaemic episodes ("time-to-event" analysis), can help resolve the probable origin of recurrences. Whole genome sequencing of predominantly low-density P. vivax infections is challenging, so an accurate and scalable genotyping method to determine the origins of recurrent parasitaemia would be of significant benefit. We have developed a P. vivax genome-wide informatics pipeline to select specific microhaplotype panels that can capture IBD within small, amplifiable segments of the genome. Using a global set of 615 P. vivax genomes, we derived a panel of 100 microhaplotypes, each comprising 3-10 high frequency SNPs within <200 bp sequence windows. This panel exhibits high diversity in regions of the Asia-Pacific, Latin America and the horn of Africa (median HE = 0.70-0.81) and it captured 89% (273/307) of the polyclonal infections detected with genome-wide datasets. Using data simulations, we demonstrate lower error in estimating pairwise IBD using microhaplotypes, relative to traditional biallelic SNP barcodes. Our panel exhibited high accuracy in predicting the country of origin (median Matthew's correlation coefficient >0.9 in 90% countries tested) and it also captured local infection outbreak and bottlenecking events. The informatics pipeline is available open-source and yields microhaplotypes that can be readily transferred to high-throughput amplicon sequencing assays for surveillance in malaria-endemic regions.

Regulators of male and female sexual development are critical for the transmission of a malaria parasite

Malaria transmission to mosquitoes requires a developmental switch in asexually dividing blood-stage parasites to sexual reproduction. In Plasmodium berghei, the transcription factor AP2-G is required and sufficient for this switch, but how a particular sex is determined in a haploid parasite remains unknown. Using a global screen of barcoded mutants, we here identify genes essential for the formation of either male or female sexual forms and validate their importance for transmission. High-resolution single-cell transcriptomics of ten mutant parasites portrays the developmental bifurcation and reveals a regulatory cascade of putative gene functions in the determination and subsequent differentiation of each sex. A male-determining gene with a LOTUS/OST-HTH domain as well as the protein interactors of a female-determining zinc-finger protein indicate that germ-granule-like ribonucleoprotein complexes complement transcriptional processes in the regulation of both male and female development of a malaria parasite.

Pf7: an open dataset of Plasmodium falciparum genome variation in 20,000 worldwide samples

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.

Characterization of a novel Plasmodium falciparum merozoite surface antigen and potential vaccine target

Introduction

Detailed analyses of genetic diversity, antigenic variability, protein localization and immunological responses are vital for the prioritization of novel malaria vaccine candidates. Comprehensive approaches to determine the most appropriate antigen variants needed to provide broad protection are challenging and consequently rarely undertaken.

Methods

Here, we characterized PF3D7_1136200, which we named Asparagine-Rich Merozoite Antigen (ARMA) based on the analysis of its sequence, localization and immunogenicity. We analyzed IgG and IgM responses against the common variants of ARMA in independent prospective cohort studies in Burkina Faso (N = 228), Kenya (N = 252) and Mali (N = 195) using a custom microarray, Div-KILCHIP.

Results

We found a marked population structure between parasites from Africa and Asia. African isolates shared 34 common haplotypes, including a dominant pair although the overall selection pressure was directional (Tajima's D = -2.57; Fu and Li's F = -9.69; P < 0.02). ARMA was localized to the merozoite surface, IgG antibodies induced Fc-mediated degranulation of natural killer cells and strongly inhibited parasite growth in vitro. We found profound serological diversity, but IgG and IgM responses were highly correlated and a hierarchical clustering analysis identified only three major serogroups. Protective IgG and IgM antibodies appeared to target both cross-reactive and distinct epitopes across variants. However, combinations of IgG and IgM antibodies against selected variants were associated with complete protection against clinical episodes of malaria.

Discussion

Our systematic strategy exploits genomic data to deduce the handful of antigen variants with the strongest potential to induce broad protection and may be broadly applicable to other complex pathogens for which effective vaccines remain elusive.

A molecular barcode and web-based data analysis tool to identify imported Plasmodium vivax malaria

Traditionally, patient travel history has been used to distinguish imported from autochthonous malaria cases, but the dormant liver stages of Plasmodium vivax confound this approach. Molecular tools offer an alternative method to identify, and map imported cases. Using machine learning approaches incorporating hierarchical fixation index and decision tree analyses applied to 799 P. vivax genomes from 21 countries, we identified 33-SNP, 50-SNP and 55-SNP barcodes (GEO33, GEO50 and GEO55), with high capacity to predict the infection's country of origin. The Matthews correlation coefficient (MCC) for an existing, commonly applied 38-SNP barcode (BR38) exceeded 0.80 in 62% countries. The GEO panels outperformed BR38, with median MCCs > 0.80 in 90% countries at GEO33, and 95% at GEO50 and GEO55. An online, open-access, likelihood-based classifier framework was established to support data analysis (vivaxGEN-geo). The SNP selection and classifier methods can be readily amended for other use cases to support malaria control programs.

Resolving drug selection and migration in an inbred South American Plasmodium falciparum population with identity-by-descent analysis

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.

Longitudinal IgG antibody responses to Plasmodium vivax blood-stage antigens during and after acute vivax malaria in individuals living in the Brazilian Amazon

BACKGROUND

To 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 FINDINGS

The 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/SIGNIFICANCE

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

Barcoding Genetically Distinct Plasmodium falciparum Strains for Comparative Assessment of Fitness and Antimalarial Drug Resistance

The repeated emergence of antimalarial drug resistance in Plasmodium falciparum, including to the current frontline antimalarial artemisinin, is a perennial problem for malaria control. Next-generation sequencing has greatly accelerated the identification of polymorphisms in resistance-associated genes but has also highlighted the need for more sensitive and accurate laboratory tools to profile current and future antimalarials and to quantify the impact of drug resistance acquisition on parasite fitness. The interplay of fitness and drug response is of fundamental importance in understanding why particular genetic backgrounds are better at driving the evolution of drug resistance in natural populations, but the impact of parasite fitness landscapes on the epidemiology of drug resistance has typically been laborious to accurately quantify in the lab, with assays being limited in accuracy and throughput. Here we present a scalable method to profile fitness and drug response of genetically distinct P. falciparum strains with well-described sensitivities to several antimalarials. We leverage CRISPR/Cas9 genome-editing and barcode sequencing to track unique barcodes integrated into a nonessential gene (pfrh3). We validate this approach in multiplex competitive growth assays of three strains with distinct geographical origins. Furthermore, we demonstrate that this method can be a powerful approach for tracking artemisinin response as it can identify an artemisinin resistant strain within a mix of multiple parasite lines, suggesting an approach for scaling the laborious ring-stage survival assay across libraries of barcoded parasite lines. Overall, we present a novel high-throughput method for multiplexed competitive growth assays to evaluate parasite fitness and drug response. IMPORTANCE The complex interplay between antimalarial resistance and parasite fitness has important implications for understanding the development and spread of drug resistance alleles and the impact of genetic background on transmission. One limitation with current methodologies to measure parasite fitness is the ability to scale this beyond simple head-to-head competition experiments between a wildtype control line and test line, with a need for a scalable approach that allows tracking of parasite growth in complex mixtures. In our study, we have used CRISPR editing to insert unique DNA barcodes into a safe-harbor genomic locus to tag multiple parasite strains and use next-generation sequencing to read out strain dynamics. We observe inherent fitness differences between the strains, as well as sensitive modulation of responses to challenge with clinically relevant antimalarials, including artemisinin.

Probing SARS-CoV-2-positive plasma to identify potential factors correlating with mild COVID-19 in Ghana, West Africa

BACKGROUND

West Africa has recorded a relatively higher proportion of asymptomatic coronavirus disease 2019 (COVID-19) cases than the rest of the world, and West Africa-specific host factors could play a role in this discrepancy. Here, we assessed the association between COVID-19 severity among Ghanaians with their immune profiles and ABO blood groups.

METHODS

Plasma samples were obtained from Ghanaians PCR-confirmed severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-positive individuals. The participants were categorized into symptomatic and asymptomatic cases. Cytokine profiling and antibody quantification were performed using Luminex™ multiplex assay whereas antigen-driven agglutination assay was used to assess the ABO blood groups. Immune profile levels between symptomatic and asymptomatic groups were compared using the two-tailed Mann-Whitney U test. Multiple comparisons of cytokine levels among and between days were tested using Kruskal-Wallis with Dunn's post hoc test. Correlations within ABO blood grouping (O's and non-O's) and between cytokines were determined using Spearman correlations. Logistic regression analysis was performed to assess the association of various cytokines with asymptomatic phenotype.

RESULTS

There was a trend linking blood group O to reduced disease severity, but this association was not statistically significant. Generally, symptomatic patients displayed significantly (p < 0.05) higher cytokine levels compared to asymptomatic cases with exception of Eotaxin, which was positively associated with asymptomatic cases. There were also significant (p < 0.05) associations between other immune markers (IL-6, IL-8 and IL-1Ra) and disease severity. Cytokines' clustering patterns differ between symptomatic and asymptomatic cases. We observed a steady decrease in the concentration of most cytokines over time, while anti-SARS-CoV-2 antibody levels were stable for at least a month, regardless of the COVID-19 status.

CONCLUSIONS

The findings suggest that genetic background and pre-existing immune response patterns may in part shape the nature of the symptomatic response against COVID-19 in a West African population. This study offers clear directions to be explored further in larger studies.

Biophysical Tools and Concepts Enable Understanding of Asexual Blood Stage Malaria

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.

Design and implementation of multiplexed amplicon sequencing panels to serve genomic epidemiology of infectious disease: a malaria case study

Multiplexed PCR amplicon sequencing (AmpSeq) is an increasingly popular application for cost-effective monitoring of threatened species and managed wildlife populations, and shows strong potential for genomic epidemiology of infectious disease. AmpSeq data from infectious microbes can inform disease control in multiple ways, including measuring drug resistance marker prevalence, distinguishing imported from local cases, and determining the effectiveness of therapeutics. We describe the design and comparative evaluation of two new AmpSeq assays for Plasmodium falciparum malaria parasites: a four-locus panel ('4CAST') composed of highly diverse antigens, and a 129-locus panel ('AMPLseq') composed of drug resistance markers, highly diverse loci for inferring relatedness, and a locus to detect Plasmodium vivax co-infection. We explore the performance of each panel in various public health use cases with in silico simulations as well as empirical experiments. The 4CAST panel appears highly suitable for evaluating the number of distinct parasite strains within samples (complexity of infection), showing strong performance across a wide range of parasitemia levels without a DNA pre-amplification step. For relatedness inference, the larger AMPLseq panel performs similarly to two existing panels of comparable size, despite differences in the data and approach used for designing each panel. Finally, we describe an R package (paneljudge) that facilitates the design and comparative evaluation of genetic panels for relatedness estimation, and we provide general guidance on the design and implementation of AmpSeq panels for genomic epidemiology of infectious disease.

An open dataset of Plasmodium vivax genome variation in 1,895 worldwide samples

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.

Controlled human malaria infection with a clone of Plasmodium vivax with high quality genome assembly

Controlled human malaria infection (CHMI) provides a highly informative means to investigate host-pathogen interactions and enable in vivo proof-of-concept efficacy testing of new drugs and vaccines. However, unlike Plasmodium falciparum, well-characterized P. vivax parasites that are safe and suitable for use in modern CHMI models are limited. Here, two healthy malaria-naïve UK adults with universal donor blood group were safely infected with a clone of P. vivax from Thailand by mosquito-bite CHMI. Parasitemia developed in both volunteers and, prior to treatment, each volunteer donated blood to produce a cryopreserved stabilate of infected red blood cells. Following stringent safety screening, the parasite stabilate from one of these donors ("PvW1") was thawed and used to inoculate six healthy malaria-naïve UK adults by blood-stage CHMI, at three different dilutions. Parasitemia developed in all volunteers, who were then successfully drug treated. PvW1 parasite DNA was isolated and sequenced to produce a high quality genome assembly by using a hybrid assembly method. We analysed leading vaccine candidate antigens and multigene families, including the Vivax interspersed repeat (VIR) genes of which we identified 1145 in the PvW1 genome. Our genomic analysis will guide future assessment of candidate vaccines and drugs, as well as experimental medicine studies.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples

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.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples

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.

Using Plasmodium knowlesi as a model for screening Plasmodium vivax blood-stage malaria vaccine targets reveals new candidates

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.

The Plasmodium falciparum Rh5 invasion protein complex reveals an excess of rare variant mutations

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.

Impact of Plasmodium falciparum small-sized extracellular vesicles on host peripheral blood mononuclear cells

Background: Exagerated immune activation has a key role in the pathogenesis of malaria. During blood-stage infection, Plasmodium falciparum can interact directly with host immune cells through infected red blood cells (PfiRBCs), or indirectly by the release of extracellular vesicles (EVs). Here, we compared the impact of PfiRBCs and P. falciparum small-sized EVs (PfsEVs, also known as exosomes) from a Kenyan clinical isolate (PfKE12) adapted to short-term laboratory culture conditions on host peripheral blood mononuclear cells (PBMC). Methods: PfsEVs were isolated from cell-free culture-conditioned media by ultracentrifugation while mature trophozoite PfiRBCs were purified by magnetic column separation. The PfsEVs and the PfiRBCs were co-cultured for 18 hours with PBMC. Cellular responses were quantified by cell surface expression of activation markers (CD25, CD69) and cytokine/chemokine levels in the supernatant. Results: Relative to negative control conditions, PfsEVs induced CD25 expression on CD4+, CD19+ and CD14+ cells, while PfiRBCs induced on CD19+ and CD14+ cells. Both PfsEVs and PfiRBCs induced CD69 on CD4+, CD8+ and CD19+ cells. In addition, PfiRBCs induced higher expression of CD69 on CD14+ cells. CD69 induced by PfiRBCs on CD4+ and CD19+ cells was significantly higher than that induced by PfsEVs. Secretion of MIP1α, MIP1β, GM-CSF, IL-6, IL-8, and TNFα were significantly induced by both PfsEVs and PfiRBCs whereas MCP-1, IL-10, IL-17α  were preferentially induced by PfsEVs and IP-10 and IFN-γ by PfiRBCs. Prior exposure to malaria (judged by antibodies to schizont extract) was associated with lower monocyte responses to PfsEVs. Conclusions: PfsEVs and PfiRBCs showed differential abilities to induce secretion of IL-17α  and IFN-γ, suggesting that the former are better at inducing Th17, whilst  the latter induce Th1 immune responses respectively. Prior exposure to malaria significantly reduces the ability of PfsEVs to activate monocytes, suggesting immune tolerance to PfsEVs may play a role in naturally acquired  anti-disease immunity.

Systematic Identification of Plasmodium Falciparum Sporozoite Membrane Protein Interactions Reveals an Essential Role for the p24 Complex in Host Infection

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.

Essential Genes of the Parasitic Apicomplexa

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.

Single cell sequencing shines a light on malaria parasite relatedness in complex infections

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.

Whole genome sequencing of Plasmodium vivax isolates reveals frequent sequence and structural polymorphisms in erythrocyte binding genes

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.

Plasmodium vivax is the most geographically widespread parasite species that causes malaria in humans. Although it occurs in Africa as a member of a mix of Plasmodium species, P. vivax is dominant in other parts of the world outside of Africa (e.g., Brazil). It was previously thought that most African populations were immune to P. vivax infections due to the absence of Duffy antigen chemokine receptor (DARC) gene expression required for erythrocyte invasion. However, several recent reports have indicated the emergence and potential spread of P. vivax across human populations in Africa. Compared to Southeast Asia and South America where P. vivax is highly endemic, data on polymorphisms in erythrocyte binding gene candidates of P. vivax from Africa is limited. Filling this knowlege gap is critical for identifying functional genes in erythrocyte invasion, biomarkers for tracking the P. vivax isolates from Africa, as well as potential gene targets for vaccine development. This paper examined the level of genetic polymorphisms in a panel of 43 potential erythrocyte binding protein genes based on whole genome sequences and described transmission patterns of P. vivax infections from different study sites in Ethiopia based on the genetic variants. Our analyses showed that chromosomes 9 and 10 of the P. vivax genomes isolated in Ethiopia had the most high-quality genetic polymorphisms. Among all erythrocyte binding protein gene candidates, the merozoite surface proteins 1 and merozoite surface protein 3 showed high levels of polymorphism. MAEBL and MSP3.8 related to immunogenicity and erythrocyte binding function were detected with significant signals of positive selection. The expansion of the Duffy binding protein and merozoite surface protein 3 gene copies was an independent process among the P. vivax lineages in Ethiopia. Various levels of gene flow were observed even at a smaller geographical scale. Our study provided baseline data for future comparison with P. vivax in Duffy negative individuals and help develop a panel of genetic markers that are informative at a micro-geographical scale.

Analysis of Plasmodium vivax schizont transcriptomes from field isolates reveals heterogeneity of expression of genes involved in host-parasite interactions

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.

Red blood cell tension protects against severe malaria in the Dantu blood group

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 regions1,2. The blood group variant Dantu provides 74% protection against all forms of severe malaria in homozygous individuals3-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.

Impact of Plasmodium falciparum small-sized extracellular vesicles on host peripheral blood mononuclear cells

Background: Exagerated immune activation has a key role in the pathogenesis of malaria. During blood-stage infection, Plasmodium falciparum can interact directly with host immune cells through infected red blood cells (PfiRBCs), or indirectly by the release of extracellular vesicles (EVs). Here, we compared the impact of PfiRBCs and P. falciparum small-sized EVs (PfsEVs, also known as exosomes) from a Kenyan clinical isolate (PfKE12) adapted to short-term laboratory culture conditions on host peripheral blood mononuclear cells (PBMC). Methods: PfsEVs were isolated from cell-free culture-conditioned media by ultracentrifugation while mature trophozoite PfiRBCs were purified by magnetic column separation. The PfsEVs and the PfiRBCs were co-cultured for 18 hours with PBMC. Cellular responses were quantified by cell surface expression of activation markers (CD25, CD69) and cytokine/chemokine levels in the supernatant. Results: Relative to negative control conditions, PfsEVs induced CD25 expression on CD4+, CD19+ and CD14+ cells, while PfiRBCs induced on CD19+ and CD14+ cells. Both PfsEVs and PfiRBCs induced CD69 on CD4+, CD8+ and CD19+ cells. In addition, PfiRBCs induced higher expression of CD69 on CD14+ cells. CD69 induced by PfiRBCs on CD4+ and CD19+ cells was significantly higher than that induced by PfsEVs. Secretion of MIP1α, MIP1β, GM-CSF, IL-6, IL-8, and TNFα were significantly induced by both PfsEVs and PfiRBCs whereas MCP-1, IL-10, IL-17α  were preferentially induced by PfsEVs and IP-10 and IFN-γ by PfiRBCs. Prior exposure to malaria (judged by antibodies to schizont extract) was associated with lower monocyte responses to PfsEVs. Conclusions: PfsEVs and PfiRBCs showed differential abilities to induce secretion of IL-17α  and IFN-γ, suggesting that the former are better at inducing Th17, whilst  the latter induce Th1 immune responses respectively. Prior exposure to malaria significantly reduces the ability of PfsEVs to activate monocytes, suggesting immune tolerance to PfsEVs may play a role in naturally acquired  anti-disease immunity.

An enhanced toolkit for the generation of knockout and marker-free fluorescent Plasmodium chabaudi

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.

An enhanced toolkit for the generation of knockout and marker-free fluorescent Plasmodium chabaudi

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.

Refining the transcriptome of the human malaria parasite Plasmodium falciparum using amplification-free RNA-seq

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.

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

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

•

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

•

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

•

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

•

Essential metabolic pathways for parasite development and mechanistic origin revealed

Molecular Characterization and Immuno-Reactivity Patterns of a Novel Plasmodium falciparum Armadillo-Type Repeat Protein, PfATRP

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.

Spatio-temporal dynamics of Plasmodium falciparum transmission within a spatial unit on the Colombian Pacific Coast

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.

Adaptation of Plasmodium falciparum to humans involved the loss of an ape-specific erythrocyte invasion ligand

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.

Here, Proto et al. show that human infective Plasmodium falciparum isolates contain an inactivating mutation in the erythrocyte invasion associated gene PfEBA165, while homologues of ape-infective Laverania species are intact, and that expression of intact PfEBA165 is incompatible with parasite growth in human erythrocytes.

Frequent expansion of Plasmodium vivax Duffy Binding Protein in Ethiopia and its epidemiological significance

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 invasion of human erythrocytes relies on interaction between the Duffy antigen and P. vivax Duffy Binding Protein (PvDBP). Whole genome sequences from P. vivax field isolates in Madagascar identified a duplication of the PvDBP gene and PvDBP duplication has also been detected in non-African P. vivax-endemic countries. Two types of PvDBP duplications have been reported, termed Cambodian and Malagasy-type duplications. Our study used a combination of PCR-based diagnostic method, a novel quantitative real-time PCR assay, and whole genome sequencing to determine the prevalence and type of PvDBP duplications, as well as PvDBP copy number on a broad number of P. vivax samples in Ethiopia. We found that over 65% of P. vivax isolated from the symptomatic infections were detected with PvDBP duplications and PvDBP varied from 1 to >4 copies. The majority of PvDBP duplications belongs to the Cambodian-type while the Malagasy-type duplications was also detected. For the asymptomatic infections, despite a small sample size, the majority of P. vivax were detected with a single-copy based on both PCR and qPCR assays. There was no significant difference in PvDBP copy number between Duffy-null heterozygote and Duffy-positive homozygote/heterozygote. Further investigation is needed with expanded Duffy-null homozygotes to examine the functional significance of PvDBP expansion.

The Malaria Cell Atlas: Single parasite transcriptomes across the complete Plasmodium life cycle

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.

Landscape of the Plasmodium Interactome Reveals Both Conserved and Species-Specific Functionality

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.

Plasmodium knowlesi as a model system for characterising Plasmodium vivax drug resistance candidate genes

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.

KILchip v1.0: A Novel Plasmodium falciparum Merozoite Protein Microarray to Facilitate Malaria Vaccine Candidate Prioritization

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.

Schizont transcriptome variation among clinical isolates and laboratory-adapted clones of the malaria parasite Plasmodium falciparum

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

A forward genetic screen reveals a primary role for Plasmodium falciparum Reticulocyte Binding Protein Homologue 2a and 2b in determining alternative erythrocyte invasion pathways

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.

Plasmodium parasites cause more than 200 million cases of malaria each year. All the symptoms of malaria are caused after Plasmodium parasites invade human red blood cells. Once inside, they grow, multiply and break open the red blood cells to release new parasites. This cycle is repeated every 48 hours, rapidly amplifying the number of parasites and causing severe anemia and other complications. Plasmodium falciparum, the parasite species responsible for almost all malaria deaths, can use multiple different pathways to invade human red blood cells, but the relative importance of each is not well understood. We tested the invasion pathways used by a collection of closely related parasites and compared their genome sequences to identify the genes responsible. This analysis revealed that expression differences in two neighboring genes of the Reticulocyte Binding Homologue family are responsible for most of the variation in two invasion pathways. P. falciparum may use variation in these genes to avoid the immune system or adapt to specific blood groups, which has important implications for vaccine development against malaria.