Strong positive selection biases identity-by-descent-based inferences of recent demography and population structure in Plasmodium falciparum.

Malaria genomic surveillance often estimates parasite genetic relatedness using metrics such as Identity-By-Decent (IBD), yet strong positive selection stemming from antimalarial drug resistance or other interventions may bias IBD-based estimates. In this study, we use simulations, a true IBD inference algorithm, and empirical data sets from different malaria transmission settings to investigate the extent of this bias and explore potential correction strategies. We analyze whole genome sequence data generated from 640 new and 3089 publicly available Plasmodium falciparum clinical isolates. We demonstrate that positive selection distorts IBD distributions, leading to underestimated effective population size and blurred population structure. Additionally, we discover that the removal of IBD peak regions partially restores the accuracy of IBD-based inferences, with this effect contingent on the population's background genetic relatedness and extent of inbreeding. Consequently, we advocate for selection correction for parasite populations undergoing strong, recent positive selection, particularly in high malaria transmission settings.

Plasmodium falciparum merozoite surface protein 1 as asexual blood stage malaria vaccine candidate.

INTRODUCTION: Malaria represents a public health challenge in tropical and subtropical regions, and currently deployed control strategies are likely insufficient to drive elimination of malaria. Development and improvement of malaria vaccines might be key to reduce disease burden. Vaccines targeting asexual blood stages of the parasite have shown limited efficacy when studied in human trials conducted over the past decades. AREAS COVERED: Vaccine candidates based on the merozoite surface protein 1 (MSP1) were initially envisioned as one of the most promising approaches to provide immune protection against asexual blood-stage malaria. Successful immunization studies in monkey involved the use of the full-length MSP1 (MSP1FL) as vaccine construct. Vaccines using MSP1FL for immunization have the potential benefit of including numerous conserved B-cell and T-cell epitopes. This could result in improved parasite strain-transcending, protective immunity in the field. We review outcomes of clinical trials that utilized a variety of MSP1 constructs and formulations, including MSP1FL, either alone or in combination with other antigens, in both animal models and humans. EXPERT OPINION: Novel approaches to analyze breadth and magnitude of effector functions of MSP1-targeting antibodies in volunteers undergoing experimental vaccination and controlled human malaria infection will help to define correlates of protective immunity.

Can incorporating genotyping data into efficacy estimators improve efficiency of early phase malaria vaccine trials?

BACKGROUND: Early phase malaria vaccine field trials typically measure malaria infection by PCR or thick blood smear microscopy performed on serially sampled blood. Vaccine efficacy (VE) is the proportion reduction in an endpoint due to vaccination and is often calculated as VEHR = 1-hazard ratio or VERR = 1-risk ratio. Genotyping information can distinguish different clones and distinguish multiple infections over time, potentially increasing statistical power. This paper investigates two alternative VE endpoints incorporating genotyping information: VEmolFOI, the vaccine-induced proportion reduction in incidence of new clones acquired over time, and VEC, the vaccine-induced proportion reduction in mean number of infecting clones per exposure. METHODS: Power of VEmolFOI and VEC was compared to that of VEHR and VERR by simulations and analytic derivations, and the four VE methods were applied to three data sets: a Phase 3 trial of RTS,S malaria vaccine in 6912 African infants, a Phase 2 trial of PfSPZ Vaccine in 80 Burkina Faso adults, and a trial comparing Plasmodium vivax incidence in 466 Papua New Guinean children after receiving chloroquine + artemether lumefantrine with or without primaquine (as these VE methods can also quantify effects of other prevention measures). By destroying hibernating liver-stage P. vivax, primaquine reduces subsequent reactivations after treatment completion. RESULTS: In the trial of RTS,S vaccine, a significantly reduced number of clones at first infection was observed, but this was not the case in trials of PfSPZ Vaccine or primaquine, although the PfSPZ trial lacked power to show a reduction. Resampling smaller data sets from the large RTS,S trial to simulate phase 2 trials showed modest power gains from VEC compared to VEHR for data like those from RTS,S, but VEC is less powerful than VEHR for trials in which the number of clones at first infection is not reduced. VEmolFOI was most powerful in model-based simulations, but only the primaquine trial collected enough serial samples to precisely estimate VEmolFOI. The primaquine VEmolFOI estimate decreased after most control arm liver-stage infections reactivated (which mathematically resembles a waning vaccine), preventing VEmolFOI from improving power. CONCLUSIONS: The power gain from the genotyping methods depends on the context. Because input parameters for early phase power calculations are often uncertain, these estimators are not recommended as primary endpoints for small trials unless supported by targeted data analysis. TRIAL REGISTRATIONS: NCT00866619, NCT02663700, NCT02143934.

The Hypervariable Tpr Multigene Family of Theileria Parasites, Defined by a Conserved, Membrane-Associated, C-Terminal Domain, Includes Several Copies with Defined Orthology Between Species.

Multigene families often play an important role in host-parasite interactions. One of the largest multigene families in Theileria parva, the causative agent of East Coast fever, is the T. parva repeat (Tpr) gene family. The function of the putative Tpr proteins remains unknown. The initial publication of the T. parva reference genome identified 39 Tpr family open reading frames (ORFs) sharing a conserved C-terminal domain. Twenty-eight of these are clustered in a central region of chromosome 3, termed the "Tpr locus", while others are dispersed throughout all four nuclear chromosomes. The Tpr locus contains three of the four assembly gaps remaining in the genome, suggesting the presence of additional, as yet uncharacterized, Tpr gene copies. Here, we describe the use of long-read sequencing to attempt to close the gaps in the reference assembly of T. parva (located among multigene families clusters), characterize the full complement of Tpr family ORFs in the T. parva reference genome, and evaluate their evolutionary relationship with Tpr homologs in other Theileria species. We identify three new Tpr family genes in the T. parva reference genome and show that sequence similarity among paralogs in the Tpr locus is significantly higher than between genes outside the Tpr locus. We also identify sequences homologous to the conserved C-terminal domain in five additional Theileria species. Using these sequences, we show that the evolution of this gene family involves conservation of a few orthologs across species, combined with gene gains/losses, and species-specific expansions.

Can incorporating genotyping data into efficacy estimators improve efficiency of early phase malaria vaccine trials?

Background: Early phase malaria vaccine field trials typically measure malaria infection by PCR or thick blood smear microscopy performed on serially sampled blood. Vaccine efficacy (VE) is the proportion reduction in an endpoint due to vaccination and is often calculated as VEHR=1 - hazard ratio or VERR=1 - risk ratio. Genotyping information can distinguish different clones and distinguish multiple infections over time, potentially increasing statistical power. This paper investigates two alternative VE endpoints incorporating genotyping information: VEmolFOI, the vaccine-induced proportion reduction in incidence of new clones acquired over time, and VEC, the vaccine-induced proportion reduction in mean number of infecting clones per exposure. Methods: We used simulations and analytic derivations to compare power of these methods to VEHR and VERR and applied them to three data sets: a Phase 3 trial of RTS,S malaria vaccine in 6912 African infants, a Phase 2 trial of PfSPZ Vaccine in 80 Burkina Faso adults, and a trial comparing Plasmodium vivax incidence in 466 Papua New Guinean children after receiving chloroquine + artemether lumefantrine with or without primaquine (as these VE methods can also quantify effects of other prevention measures). By destroying hibernating liver-stage P. vivax, primaquine reduces subsequent reactivations after treatment completion. Results: The RTS,S vaccine significantly reduced the number of clones at first infection, but PfSPZ vaccine and primaquine did not. Resampling smaller data sets from the large RTS,S trial to simulate phase 2 trials showed modest power gains from VEC compared to VEHR for data like RTS,S, but VEC is less powerful than VEHR for vaccines which do not reduce the number of clones at first infection. VEmolFOI was most powerful in model-based simulations, but only the primaquine trial collected enough serial samples to precisely estimate VEmolFOI. The primaquine VEmolFOI estimate decreased after most control arm liver-stage infections reactivated (which mathematically resembles a waning vaccine), preventing VEmolFOI from improving power. Conclusions: The power gain from the genotyping methods depends on the context. Because input parameters for early phase power calculations are often uncertain, we recommend against these estimators as primary endpoints for small trials unless supported by targeted data analysis. Trial registrations: NCT00866619, NCT02663700, NCT02143934.

Sporozoite immunization: innovative translational science to support the fight against malaria.

INTRODUCTION: Malaria, a devastating febrile illness caused by protozoan parasites, sickened 247,000,000 people in 2021 and killed 619,000, mostly children and pregnant women in sub-Saharan Africa. A highly effective vaccine is urgently needed, especially for Plasmodium falciparum (Pf), the deadliest human malaria parasite. AREAS COVERED: Sporozoites (SPZ), the parasite stage transmitted by Anopheles mosquitoes to humans, are the only vaccine immunogen achieving >90% efficacy against Pf infection. This review describes >30 clinical trials of PfSPZ vaccines in the U.S.A., Europe, Africa, and Asia, based on first-hand knowledge of the trials and PubMed searches of 'sporozoites,' 'malaria,' and 'vaccines.' EXPERT OPINION: First generation (radiation-attenuated) PfSPZ vaccines are safe, well tolerated, 80-100% efficacious against homologous controlled human malaria infection (CHMI) and provide 18-19 months protection without boosting in Africa. Second generation chemo-attenuated PfSPZ are more potent, 100% efficacious against stringent heterologous (variant strain) CHMI, but require a co-administered drug, raising safety concerns. Third generation, late liver stage-arresting, replication competent (LARC), genetically-attenuated PfSPZ are expected to be both safe and highly efficacious. Overall, PfSPZ vaccines meet safety, tolerability, and efficacy requirements for protecting pregnant women and travelers exposed to Pf in Africa, with licensure for these populations possible within 5 years. Protecting children and mass vaccination programs to block transmission and eliminate malaria are long-term objectives.

Strong Positive Selection Biases Identity-By-Descent-Based Inferences of Recent Demography and Population Structure in Plasmodium falciparum.

Malaria genomic surveillance often estimates parasite genetic relatedness using metrics such as Identity-By-Decent (IBD). Yet, strong positive selection stemming from antimalarial drug resistance or other interventions may bias IBD-based estimates. In this study, we utilized simulations, a true IBD inference algorithm, and empirical datasets from different malaria transmission settings to investigate the extent of such bias and explore potential correction strategies. We analyzed whole genome sequence data generated from 640 new and 4,026 publicly available Plasmodium falciparum clinical isolates. Our findings demonstrated that positive selection distorts IBD distributions, leading to underestimated effective population size and blurred population structure. Additionally, we discovered that the removal of IBD peak regions partially restored the accuracy of IBD-based inferences, with this effect contingent on the population's background genetic relatedness. Consequently, we advocate for selection correction for parasite populations undergoing strong, recent positive selection, particularly in high malaria transmission settings.

Immunogenomic profile at baseline predicts host susceptibility to clinical malaria.

Introduction: Host gene and protein expression impact susceptibility to clinical malaria, but the balance of immune cell populations, cytokines and genes that contributes to protection, remains incompletely understood. Little is known about the determinants of host susceptibility to clinical malaria at a time when acquired immunity is developing. Methods: We analyzed peripheral blood mononuclear cells (PBMCs) collected from children who differed in susceptibility to clinical malaria, all from a small town in Mali. PBMCs were collected from children aged 4-6 years at the start, peak and end of the malaria season. We characterized the immune cell composition and cytokine secretion for a subset of 20 children per timepoint (10 children with no symptomatic malaria age-matched to 10 children with >2 symptomatic malarial illnesses), and gene expression patterns for six children (three per cohort) per timepoint. Results: We observed differences between the two groups of children in the expression of genes related to cell death and inflammation; in particular, inflammatory genes such as CXCL10 and STAT1 and apoptotic genes such as XAF1 were upregulated in susceptible children before the transmission season began. We also noted higher frequency of HLA-DR+ CD4 T cells in protected children during the peak of the malaria season and comparable levels cytokine secretion after stimulation with malaria schizonts across all three time points. Conclusion: This study highlights the importance of baseline immune signatures in determining disease outcome. Our data suggests that differences in apoptotic and inflammatory gene expression patterns can serve as predictive markers of susceptibility to clinical malaria.

The exception that proves the rule: Virulence gene expression at the onset of Plasmodium falciparum blood stage infections.

Controlled human malaria infections (CHMI) are a valuable tool to study parasite gene expression in vivo under defined conditions. In previous studies, virulence gene expression was analyzed in samples from volunteers infected with the Plasmodium falciparum (Pf) NF54 isolate, which is of African origin. Here, we provide an in-depth investigation of parasite virulence gene expression in malaria-naïve European volunteers undergoing CHMI with the genetically distinct Pf 7G8 clone, originating in Brazil. Differential expression of var genes, encoding major virulence factors of Pf, PfEMP1s, was assessed in ex vivo parasite samples as well as in parasites from the in vitro cell bank culture that was used to generate the sporozoites (SPZ) for CHMI (Sanaria PfSPZ Challenge (7G8)). We report broad activation of mainly B-type subtelomeric located var genes at the onset of a 7G8 blood stage infection in naïve volunteers, mirroring the NF54 expression study and suggesting that the expression of virulence-associated genes is generally reset during transmission from the mosquito to the human host. However, in 7G8 parasites, we additionally detected a continuously expressed single C-type variant, Pf7G8_040025600, that was most highly expressed in both pre-mosquito cell bank and volunteer samples, suggesting that 7G8, unlike NF54, maintains expression of some previously expressed var variants during transmission. This suggests that in a new host, the parasite may preferentially express the variants that previously allowed successful infection and transmission. Trial registration: ClinicalTrials.gov - NCT02704533; 2018-004523-36.

Role of Pfs47 in the dispersal of ancestral Plasmodium falciparum malaria through adaptation to different anopheline vectors.

Plasmodium falciparum malaria originated when Plasmodium praefalciparum, a gorilla malaria parasite transmitted by African sylvan anopheline mosquitoes, adapted to humans. Pfs47, a protein on the parasite surface mediates P. falciparum evasion of the mosquito immune system by interacting with a midgut receptor and is critical for Plasmodium adaptation to different anopheline species. Genetic analysis of 4,971 Pfs47 gene sequences from different continents revealed that Asia and Papua New Guinea harbor Pfs47 haplotypes more similar to its ortholog in P. praefalciparum at sites that determine vector compatibility, suggesting that ancestral P. falciparum readily adapted to Asian vectors. Consistent with this observation, Pfs47-receptor gene sequences from African sylvan malaria vectors, such as Anopheles moucheti and An. marshallii, were found to share greater similarity with those of Asian vectors than those of vectors of the African An. gambiae complex. Furthermore, experimental infections provide direct evidence that transformed P. falciparum parasites carrying Pfs47 orthologs of P. praefalciparum or P. reichenowi were more effective at evading the immune system of the Asian malaria vector An. dirus than An. gambiae. We propose that high compatibility of ancestral P. falciparum Pfs47 with the receptors of Asian vectors facilitated the early dispersal of human malaria to the Asian continent, without having to first adapt to sub-Saharan vectors of the An. gambiae complex.

Malaria in Angola: recent progress, challenges and future opportunities using parasite demography studies.

Over the past two decades, a considerable expansion of malaria interventions has occurred at the national level in Angola, together with cross-border initiatives and regional efforts in southern Africa. Currently, Angola aims to consolidate malaria control and to accelerate the transition from control to pre-elimination, along with other country members of the Elimination 8 initiative. However, the tremendous heterogeneity in malaria prevalence among Angolan provinces, as well as internal population movements and migration across borders, represent major challenges for the Angolan National Malaria Control Programme. This review aims to contribute to the understanding of factors underlying the complex malaria situation in Angola and to encourage future research studies on transmission dynamics and population structure of Plasmodium falciparum, important areas to complement host epidemiological information and to help reenergize the goal of malaria elimination in the country.

A genetically modified Plasmodium berghei parasite as a surrogate for whole-sporozoite vaccination against P. vivax malaria.

Two malaria parasite species, Plasmodium falciparum (Pf) and P. vivax (Pv) are responsible for most of the disease burden caused by malaria. Vaccine development against this disease has focused mainly on Pf. Whole-sporozoite (WSp) vaccination, targeting pre-erythrocytic (PE) parasite stages, is a promising strategy for immunization against malaria and several PfWSp-based vaccine candidates are currently undergoing clinical evaluation. In contrast, no WSp candidates have been developed for Pv, mainly due to constraints in the production of Pv sporozoites in the laboratory. Recently, we developed a novel approach for WSp vaccination against Pf based on the use of transgenic rodent P. berghei (Pb) sporozoites expressing immunogens of this human-infective parasite. We showed that this platform can be used to deliver PE Pf antigens, eliciting both targeted humoral responses and cross-species cellular immune responses against Pf. Here we explored this WSp platform for the delivery of Pv antigens. As the Pv circumsporozoite protein (CSP) is a leading vaccine candidate antigen, we generated a transgenic Pb parasite, PbviVac, that, in addition to its endogenous PbCSP, expresses PvCSP under the control of a strictly PE promoter. Immunofluorescence microscopy analyses confirmed that both the PbCSP and the PvCSP antigens are expressed in PbviVac sporozoites and liver stages and that PbviVac sporozoite infectivity of hepatic cells is similar to that of its wild-type Pb counterpart. Immunization of mice with PbviVac sporozoites elicits the production of anti-PvCSP antibodies that efficiently recognize and bind to Pv sporozoites. Our results warrant further development and evaluation of PbviVac as a surrogate for WSp vaccination against Pv malaria.

Gene Coverage Count and Classification (GC3): a locus sequence coverage assessment tool using short-read whole genome sequencing data, and its application to identify and classify histidine-rich protein 2 and 3 deletions in Plasmodium falciparum.

BACKGROUND: The ability of malaria rapid diagnostic tests (RDTs) to effectively detect active infections is being compromised by the presence of malaria strains with genomic deletions at the hrp2 and hrp3 loci, encoding the antigens most commonly targeted in diagnostics for Plasmodium falciparum detection. The presence of such deletions can be determined in publically available P. falciparum whole genome sequencing (WGS) datasets. A computational approach was developed and validated, termed Gene Coverage Count and Classification (GC3), to analyse genome-wide sequence coverage data and provide informative outputs to assess presence and coverage profile of a target locus in WGS data. GC3 was applied to detect deletions at hrp2 and hrp3 (hrp2/3) and flanking genes in different geographic regions and across time points. METHODS: GC3 uses Python and R scripts to extract locus read coverage metrics from mapped WGS data according to user-defined parameters and generates relevant tables and figures. GC3 was tested using WGS data for laboratory reference strains with known hrp2/3 genotypes, and its results compared to those of a hrp2/3-specific qPCR assay. Samples with at least 25% of coding region positions with zero coverage were classified as having a deletion. Publicly available sequence data was analysed and compared with published deletion frequency estimates. RESULTS: GC3 results matched the expected coverage of known laboratory reference strains. Agreement between GC3 and a hrp2/3-specific qPCR assay reported for 19/19 (100%) hrp2 deletions and 18/19 (94.7%) hrp3 deletions. Among Cambodian (n = 127) and Brazilian (n = 20) WGS datasets, which had not been previously analysed for hrp2/3 deletions, GC3 identified hrp2 deletions in three and four samples, and hrp3 deletions in 10 and 15 samples, respectively. Plots of hrp2/3 coding regions, grouped by year of sample collection, showed a decrease in median standardized coverage among Malawian samples (n = 150) suggesting the importance of a careful, properly controlled follow up to determine if an increase in frequency of deletions has occurred between 2007-2008 and 2014-2015. Among Malian (n = 90) samples, median standardized coverage was lower in 2002 than 2010, indicating widespread deletions present at the gene locus in 2002. CONCLUSIONS: The GC3 tool accurately classified hrp2/3 deletions and provided informative tables and figures to analyse targeted gene coverage. GC3 is an appropriate tool when performing preliminary and exploratory assessment of locus coverage data.

An In Silico Analysis of Malaria Pre-Erythrocytic-Stage Antigens Interpreting Worldwide Genetic Data to Suggest Vaccine Candidate Variants and Epitopes.

Failure to account for genetic diversity of antigens during vaccine design may lead to vaccine escape. To evaluate the vaccine escape potential of antigens used in vaccines currently in development or clinical testing, we surveyed the genetic diversity, measured population differentiation, and performed in silico prediction and analysis of T-cell epitopes of ten such Plasmodium falciparum pre-erythrocytic-stage antigens using whole-genome sequence data from 1010 field isolates. Of these, 699 were collected in Africa (Burkina Faso, Cameroon, Guinea, Kenya, Malawi, Mali, and Tanzania), 69 in South America (Brazil, Colombia, French Guiana, and Peru), 59 in Oceania (Papua New Guinea), and 183 in Asia (Cambodia, Myanmar, and Thailand). Antigens surveyed include cell-traversal protein for ookinetes and sporozoites, circumsporozoite protein, liver-stage antigens 1 and 3, sporozoite surface proteins P36 and P52, sporozoite asparagine-rich protein-1, sporozoite microneme protein essential for cell traversal-2, and upregulated-in-infectious-sporozoite 3 and 4 proteins. The analyses showed that a limited number of these protein variants, when combined, would be representative of worldwide parasite populations. Moreover, predicted T-cell epitopes were identified that could be further explored for immunogenicity and protective efficacy. Findings can inform the rational design of a multivalent malaria vaccine.

Plasmodium falciparum 7G8 challenge provides conservative prediction of efficacy of PfNF54-based PfSPZ Vaccine in Africa.

Controlled human malaria infection (CHMI) has supported Plasmodium falciparum (Pf) malaria vaccine development by providing preliminary estimates of vaccine efficacy (VE). Because CHMIs generally use Pf strains similar to vaccine strains, VE against antigenically heterogeneous Pf in the field has been required to establish VE. We increased the stringency of CHMI by selecting a Brazilian isolate, Pf7G8, which is genetically distant from the West African parasite (PfNF54) in our PfSPZ vaccines. Using two regimens to identically immunize US and Malian adults, VE over 24 weeks in the field was as good as or better than VE against CHMI at 24 weeks in the US. To explain this finding, here we quantify differences in the genome, proteome, and predicted CD8 T cell epitopes of PfNF54 relative to 704 Pf isolates from Africa and Pf7G8. We show that Pf7G8 is more distant from PfNF54 than any African isolates tested. We propose VE against Pf7G8 CHMI for providing pivotal data for malaria vaccine licensure for travelers to Africa, and potentially for endemic populations, because the genetic distance of Pf7G8 from the Pf vaccine strain makes it a stringent surrogate for Pf parasites in Africa.

CMTR1 is recruited to transcription start sites and promotes ribosomal protein and histone gene expression in embryonic stem cells.

CMTR1 (cap methyltransferase 1) catalyses methylation of the first transcribed nucleotide of RNAPII transcripts (N1 2'-O-Me), creating part of the mammalian RNA cap structure. In addition to marking RNA as self, N1 2'-O-Me has ill-defined roles in RNA expression and translation. Here, we investigated the gene specificity of CMTR1 and its impact on RNA expression in embryonic stem cells. Using chromatin immunoprecipitation, CMTR1 was found to bind to transcription start sites (TSS) correlating with RNAPII levels, predominantly binding at histone genes and ribosomal protein (RP) genes. Repression of CMTR1 expression resulted in repression of RNAPII binding at the TSS and repression of RNA expression, particularly of histone and RP genes. In correlation with regulation of histones and RP genes, CMTR1 repression resulted in repression of translation and induction of DNA replication stress and damage. Indicating a direct role for CMTR1 in transcription, addition of recombinant CMTR1 to purified nuclei increased transcription of the histone and RP genes. CMTR1 was found to be upregulated during neural differentiation and there was an enhanced requirement for CMTR1 for gene expression and proliferation during this process. We highlight the distinct roles of the cap methyltransferases RNMT and CMTR1 in target gene expression and differentiation.

STRIDE: a command-line HMM-based identifier and sub-classifier of Plasmodium falciparum RIFIN and STEVOR variant surface antigen families.

BACKGROUND: RIFINs and STEVORs are variant surface antigens expressed by P. falciparum that play roles in severe malaria pathogenesis and immune evasion. These two highly diverse multigene families feature multiple paralogs, making their classification challenging using traditional bioinformatic methods. RESULTS: STRIDE (STevor and RIfin iDEntifier) is an HMM-based, command-line program that automates the identification and classification of RIFIN and STEVOR protein sequences in the malaria parasite Plasmodium falciparum. STRIDE is more sensitive in detecting RIFINs and STEVORs than available PFAM and TIGRFAM tools and reports RIFIN subtypes and the number of sequences with a FHEYDER amino acid motif, which has been associated with severe malaria pathogenesis. CONCLUSIONS: STRIDE will be beneficial to malaria research groups analyzing genome sequences and transcripts of clinical field isolates, providing insight into parasite biology and virulence.