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

Ancient Plasmodium genomes shed light on the history of human malaria

Malaria-causing protozoa of the genus Plasmodium have exerted one of the strongest selective pressures on the human genome, and resistance alleles provide biomolecular footprints that outline the historical reach of these species1. Nevertheless, debate persists over when and how malaria parasites emerged as human pathogens and spread around the globe1,2. To address these questions, we generated high-coverage ancient mitochondrial and nuclear genome-wide data from P. falciparum, P. vivax and P. malariae from 16 countries spanning around 5,500 years of human history. We identified P. vivax and P. falciparum across geographically disparate regions of Eurasia from as early as the fourth and first millennia BCE, respectively; for P. vivax, this evidence pre-dates textual references by several millennia3. Genomic analysis supports distinct disease histories for P. falciparum and P. vivax in the Americas: similarities between now-eliminated European and peri-contact South American strains indicate that European colonizers were the source of American P. vivax, whereas the trans-Atlantic slave trade probably introduced P. falciparum into the Americas. Our data underscore the role of cross-cultural contacts in the dissemination of malaria, laying the biomolecular foundation for future palaeo-epidemiological research into the impact of Plasmodium parasites on human history. Finally, our unexpected discovery of P. falciparum in the high-altitude Himalayas provides a rare case study in which individual mobility can be inferred from infection status, adding to our knowledge of cross-cultural connectivity in the region nearly three millennia ago.

Genetic polymorphism and evidence of signatures of selection in the Plasmodium falciparum circumsporozoite protein gene in Tanzanian regions with different malaria endemicity

BACKGROUND

In 2021 and 2023, the World Health Organization approved RTS,S/AS01 and R21/Matrix M malaria vaccines, respectively, for routine immunization of children in African countries with moderate to high transmission. These vaccines are made of Plasmodium falciparum circumsporozoite protein (PfCSP), but polymorphisms in the gene raise concerns regarding strain-specific responses and the long-term efficacy of these vaccines. This study assessed the Pfcsp genetic diversity, population structure and signatures of selection among parasites from areas of different malaria transmission intensities in Mainland Tanzania, to generate baseline data before the introduction of the malaria vaccines in the country.

METHODS

The analysis involved 589 whole genome sequences generated by and as part of the MalariaGEN Community Project. The samples were collected between 2013 and January 2015 from five regions of Mainland Tanzania: Morogoro and Tanga (Muheza) (moderate transmission areas), and Kagera (Muleba), Lindi (Nachingwea), and Kigoma (Ujiji) (high transmission areas). Wright's inbreeding coefficient (Fws), Wright's fixation index (FST), principal component analysis, nucleotide diversity, and Tajima's D were used to assess within-host parasite diversity, population structure and natural selection.

RESULTS

Based on Fws (< 0.95), there was high polyclonality (ranging from 69.23% in Nachingwea to 56.9% in Muheza). No population structure was detected in the Pfcsp gene in the five regions (mean FST = 0.0068). The average nucleotide diversity (π), nucleotide differentiation (K) and haplotype diversity (Hd) in the five regions were 4.19, 0.973 and 0.0035, respectively. The C-terminal region of Pfcsp showed high nucleotide diversity at Th2R and Th3R regions. Positive values for the Tajima's D were observed in the Th2R and Th3R regions consistent with balancing selection. The Pfcsp C-terminal sequences revealed 50 different haplotypes (H_1 to H_50), with only 2% of sequences matching the 3D7 strain haplotype (H_50). Conversely, with the NF54 strain, the Pfcsp C-terminal sequences revealed 49 different haplotypes (H_1 to H_49), with only 0.4% of the sequences matching the NF54 strain (Hap_49).

CONCLUSIONS

The findings demonstrate high diversity of the Pfcsp gene with limited population differentiation. The Pfcsp gene showed positive Tajima's D values, consistent with balancing selection for variants within Th2R and Th3R regions. The study observed differences between the intended haplotypes incorporated into the design of RTS,S and R21 vaccines and those present in natural parasite populations. Therefore, additional research is warranted, incorporating other regions and more recent data to comprehensively assess trends in genetic diversity within this important gene. Such insights will inform the choice of alleles to be included in the future vaccines.

Genomics reveals heterogeneous Plasmodium falciparum transmission and selection signals in Zambia

BACKGROUND

Genomic surveillance is crucial for monitoring malaria transmission and understanding parasite adaptation to interventions. Zambia lacks prior nationwide efforts in malaria genomic surveillance among African countries.

METHODS

We conducted genomic surveillance of Plasmodium falciparum parasites from the 2018 Malaria Indicator Survey in Zambia, a nationally representative household survey of children under five years of age. We whole-genome sequenced and analyzed 241 P. falciparum genomes from regions with varying levels of malaria transmission across Zambia and estimated genetic metrics that are informative about transmission intensity, genetic relatedness between parasites, and selection.

RESULTS

We provide genomic evidence of widespread within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in Zambia. Our analysis reveals country-level clustering of parasites from Zambia and neighboring regions, with distinct separation in West Africa. Within Zambia, identity by descent (IBD) relatedness analysis uncovers local spatial clustering and rare cases of long-distance sharing of closely related parasite pairs. Genomic regions with large shared IBD segments and strong positive selection signatures implicate genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Furthermore, differences in selection signatures, including drug resistance loci, are observed between eastern and western Zambian parasite populations, suggesting variable transmission intensity and ongoing drug pressure.

CONCLUSIONS

Our findings enhance our understanding of nationwide P. falciparum transmission in Zambia, establishing a baseline for analyzing parasite genetic metrics as they vary over time and space. These insights highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

Revisiting the determinants of malaria transmission

Malaria parasites have coevolved with humans over thousands of years, mirroring their migration out of Africa. They persist to this day, despite continuous elimination efforts worldwide. These parasites can adapt to changing environments during infection of human and mosquito, and when expanding the geographical range by switching vector species. Recent studies in the human malaria parasite, Plasmodium falciparum, identified determinants governing the plasticity of sexual conversion rates, sex ratio, and vector competence. Here we summarize the latest literature revealing environmental, epigenetic, and genetic determinants of malaria transmission.

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 transmission in the highlands of Ethiopia is driven by closely related and clonal parasites

Malaria cases are frequently recorded in the Ethiopian highlands even at altitudes above 2000 m. The epidemiology of malaria in the Ethiopian highlands, and, in particular, the role of importation by human migration from the highly endemic lowlands is not well understood. We sequenced 187 Plasmodium falciparum samples from two sites in the Ethiopian highlands, Gondar (n = 159) and Ziway (n = 28), using a multiplexed droplet digital PCR (ddPCR)-based amplicon sequencing method targeting 35 microhaplotypes and drug resistance loci. Here, we characterize the parasite population structure and genetic relatedness. We identify moderate parasite diversity (mean HE : 0.54) and low infection complexity (74.9% monoclonal). A significant percentage of infections share microhaplotypes, even across transmission seasons and sites, indicating persistent local transmission. We identify multiple clusters of clonal or near-clonal infections, highlighting high genetic relatedness. Only 6.3% of individuals diagnosed with P. falciparum reported recent travel. Yet, in clonal or near-clonal clusters, infections of travellers were frequently observed first in time, suggesting that parasites may have been imported and then transmitted locally. 31.1% of infections are pfhrp2-deleted and 84.4% pfhrp3-deleted, and 28.7% have pfhrp2/3 double deletions. Parasites with pfhrp2/3 deletions and wild-type parasites are genetically distinct. Mutations associated with resistance to sulphadoxine-pyrimethamine or suggested to reduce sensitivity to lumefantrine are observed at near-fixation. In conclusion, genomic data corroborate local transmission and the importance of intensified control in the Ethiopian highlands.

Role for gene conversion in the evolution of cell-surface antigens of the malaria parasite Plasmodium falciparum

While the malaria parasite Plasmodium falciparum has low average genome-wide diversity levels, likely due to its recent introduction from a gorilla-infecting ancestor (approximately 10,000 to 50,000 years ago), some genes display extremely high diversity levels. In particular, certain proteins expressed on the surface of human red blood cell-infecting merozoites (merozoite surface proteins (MSPs)) possess exactly 2 deeply diverged lineages that have seemingly not recombined. While of considerable interest, the evolutionary origin of this phenomenon remains unknown. In this study, we analysed the genetic diversity of 2 of the most variable MSPs, DBLMSP and DBLMSP2, which are paralogs (descended from an ancestral duplication). Despite thousands of available Illumina WGS datasets from malaria-endemic countries, diversity in these genes has been hard to characterise as reads containing highly diverged alleles completely fail to align to the reference genome. To solve this, we developed a pipeline leveraging genome graphs, enabling us to genotype them at high accuracy and completeness. Using our newly- resolved sequences, we found that both genes exhibit 2 deeply diverged lineages in a specific protein domain (DBL) and that one of the 2 lineages is shared across the genes. We identified clear evidence of nonallelic gene conversion between the 2 genes as the likely mechanism behind sharing, leading us to propose that gene conversion between diverged paralogs, and not recombination suppression, can generate this surprising genealogy; a model that is furthermore consistent with high diversity levels in these 2 genes despite the strong historical P. falciparum transmission bottleneck.

New reference genomes to distinguish the sympatric malaria parasites, Plasmodium ovale curtisi and Plasmodium ovale wallikeri

Despite Plasmodium ovale curtisi (Poc) and wallikeri (Pow) being important human-infecting malaria parasites that are widespread across Africa and Asia, little is known about their genome diversity. Morphologically identical, Poc and Pow are indistinguishable and commonly misidentified. Recent rises in the incidence of Poc/Pow infections have renewed efforts to address fundamental knowledge gaps in their biology, and to develop diagnostic tools to understand their epidemiological dynamics and malaria burden. A major roadblock has been the incompleteness of available reference assemblies (PocGH01, PowCR01; ~ 33.5 Mbp). Here, we applied multiple sequencing platforms and advanced bioinformatics tools to generate new reference genomes, Poc221 (South Sudan; 36.0 Mbp) and Pow222 (Nigeria; 34.3 Mbp), with improved nuclear genome contiguity (> 4.2 Mbp), annotation and completeness (> 99% Plasmodium spp., single copy orthologs). Subsequent sequencing of 6 Poc and 15 Pow isolates from Africa revealed a total of 22,517 and 43,855 high-quality core genome SNPs, respectively. Genome-wide levels of nucleotide diversity were determined to be 2.98 × 10-4 (Poc) and 3.43 × 10-4 (Pow), comparable to estimates for other Plasmodium species. Overall, the new reference genomes provide a robust foundation for dissecting the biology of Poc/Pow, their population structure and evolution, and will contribute to uncovering the recombination barrier separating these species.

Molecular markers for malaria genetic epidemiology: progress and pitfalls

Over recent years, progress in molecular markers for genotyping malaria parasites has enabled informative studies of epidemiology and transmission dynamics. Results have highlighted the value of these tools for surveillance to support malaria control and elimination strategies. There are many different types and panels of markers available for malaria parasite genotyping, and for end users, the nuances of these markers with respect to 'use case', resolution, and accuracy, are not well defined. This review clarifies issues surrounding different molecular markers and their application to malaria control and elimination. We describe available marker panels, use cases, implications for different transmission settings, limitations, access, cost, and data accuracy. The information provided can be used as a guide for molecular epidemiology and surveillance of malaria.

Genetic polymorphism and evidence of signatures of selection in the Plasmodium falciparum circumsporozoite protein gene in Tanzanian regions with different malaria endemicity

Background

In 2021 and 2023, the World Health Organization approved RTS, S/AS01 and R21/Matrix M malaria vaccines, respectively, for routine immunization of children in African countries with moderate to high transmission. These vaccines are made of Plasmodium falciparum circumsporozoite protein (Pfcsp) but polymorphisms in this gene raises concerns regarding strain-specific responses and the long-term efficacy of these vaccines. This study assessed the Pfcsp genetic diversity, population structure and signatures of selection among parasites from areas of different malaria transmission in mainland Tanzania, to generate baseline data before the introduction of the malaria vaccines in the country.

Methods

The analysis involved 589 whole genome sequences generated by and as part of the MalariaGEN Community Project. The samples were collected between 2013 and January 2015 from five regions of mainland Tanzania: Morogoro and Tanga (Muheza) (moderate transmission areas), and Kagera (Muleba), Lindi (Nachingwea), and Kigoma (Ujiji) (high transmission areas). Wright's inbreeding coefficient (Fws), Wright's fixation index (FST), principal component analysis, nucleotide diversity, and Tajima's D were used to assess within-host parasite diversity, population structure and natural selection.

Results

Based on Fws (< 0.95), there was high polyclonality (ranged from 69.23% in Nachingwea to 56.9% in Muheza). No population structure was detected in the Pfcsp gene in the five regions (mean FST= 0.0068). The average nucleotide diversity (π), nucleotide differentiation (K) and haplotype diversity (Hd) in the five regions were 4.19, 0.973 and 0.0035, respectively. The C-terminal region of Pfcsp showed high nucleotide diversity at Th2R and Th3R regions. Positive values for the Tajima's D were observed in the Th2R and Th3R regions consistent with balancing selection. The Pfcsp C-terminal sequences had 50 different haplotypes (H_1 to H_50) and only 2% of sequences matched the 3D7 strain haplotype (H_50).

Conclusions

The findings demonstrate high diversity of the Pfcsp gene with limited population differentiation. The Pfcsp gene showed positive Tajima's D values for parasite populations, consistent with balancing selection for variants within Th2R and Th3R regions. This data is consistent with other studies conducted across Africa and worldwide, which demonstrate low 3D7 haplotypes and little population structure. Therefore, additional research is warranted, incorporating other regions and more recent data to comprehensively assess trends in genetic diversity within this important gene. Such insights will inform the choice of alleles to be included in the future vaccines.

Evolution and spread of Plasmodium falciparum mutations associated with resistance to sulfadoxine-pyrimethamine in central Africa: a cross-sectional study

BACKGROUND

Efficacy of sulfadoxine-pyrimethamine, the malaria chemoprophylaxis used in pregnant women, and in children when combined with amodiaquine, is threatened by the accumulation of mutations in the Plasmodium falciparum dihydropteroate synthase (pfdhps) and dihydrofolate reductase (pfdhfr) genes. Data on the prevalence of resistant alleles in central Africa and the new pfdhps I431V mutation, particularly associated with other mutations to form the pfdhps vagKgs allele, are scarce. We explored the frequency and geographical distribution of pfdhps and pfdhfr mutations in central Africa in 2014-18, and assessed the evolutionary origin of the vagKgs allele.

METHODS

Samples were collected at 18 health-care centres in seven countries (Angola, Cameroon, Central African Republic, Democratic Republic of the Congo, Gabon, Nigeria, and Republic of the Congo) from patients who showed possible symptoms of malaria between March 1, 2014, and Oct 31, 2018. Samples that were positive for P falciparum were transported to a laboratory in Toulouse, France, and genotyped. The frequency of pfdhfr and pfdhps mutations was studied in 1749 samples. Microsatellites in pfdhps flanking regions and whole-genome analysis compared with parasite genomes from the data-sharing network MalariaGEN were performed on samples carrying the vagKgs allele.

FINDINGS

Mapping of the prevalence of single nucleotide polymorphisms and corresponding alleles of pfdhfr and pfdhps showed a substantial spread of alleles associated with sulfadoxine-pyrimethamine resistance in central Africa during the 2014-18 period, especially an increase going west to east in pfdhps alleles carrying the K540E and A581G mutations. A high prevalence of the pfdhps I431V mutation was observed in Cameroon (exceeding 50% in the northern region) and Nigeria. Genomic analysis showed a recent African emergence and a clonal expansion of the most frequent pfdhps vagKgs allele.

INTERPRETATION

Reduced sulfadoxine-pyrimethamine efficacy due to increased resistance is a worrying situation, especially because the malaria transmission level is high in central Africa. Although the resistance phenotype remains to be confirmed, the emergence and spread of the vagKgs allele in west and central Africa could challenge the use of sulfadoxine-pyrimethamine.

FUNDING

Toulouse Institute for Infectious and Inflammatory Diseases.

The Kelch13 compartment contains highly divergent vesicle trafficking proteins in malaria parasites

Single amino acid changes in the parasite protein Kelch13 (K13) result in reduced susceptibility of P. falciparum parasites to artemisinin and its derivatives (ART). Recent work indicated that K13 and other proteins co-localising with K13 (K13 compartment proteins) are involved in the endocytic uptake of host cell cytosol (HCCU) and that a reduction in HCCU results in reduced susceptibility to ART. HCCU is critical for parasite survival but is poorly understood, with the K13 compartment proteins among the few proteins so far functionally linked to this process. Here we further defined the composition of the K13 compartment by analysing more hits from a previous BioID, showing that MyoF and MCA2 as well as Kelch13 interaction candidate (KIC) 11 and 12 are found at this site. Functional analyses, tests for ART susceptibility as well as comparisons of structural similarities using AlphaFold2 predictions of these and previously identified proteins showed that vesicle trafficking and endocytosis domains were frequent in proteins involved in resistance or endocytosis (or both), comprising one group of K13 compartment proteins. While this strengthened the link of the K13 compartment to endocytosis, many proteins of this group showed unusual domain combinations and large parasite-specific regions, indicating a high level of taxon-specific adaptation of this process. Another group of K13 compartment proteins did not influence endocytosis or ART susceptibility and lacked detectable vesicle trafficking domains. We here identified the first protein of this group that is important for asexual blood stage development and showed that it likely is involved in invasion. Overall, this work identified novel proteins functioning in endocytosis and at the K13 compartment. Together with comparisons of structural predictions it provides a repertoire of functional domains at the K13 compartment that indicate a high level of adaption of endocytosis in malaria parasites.

Drug resistance and vaccine target surveillance of Plasmodium falciparum using nanopore sequencing in Ghana

Malaria results in over 600,000 deaths annually, with the highest burden of deaths in young children living in sub-Saharan Africa. Molecular surveillance can provide important information for malaria control policies, including detection of antimalarial drug resistance. However, genome sequencing capacity in malaria-endemic countries is limited. We designed and implemented an end-to-end workflow to detect Plasmodium falciparum antimalarial resistance markers and diversity in the vaccine target circumsporozoite protein (csp) using nanopore sequencing in Ghana. We analysed 196 clinical samples and showed that our method is rapid, robust, accurate and straightforward to implement. Importantly, our method could be applied to dried blood spot samples, which are readily collected in endemic settings. We report that P. falciparum parasites in Ghana are mostly susceptible to chloroquine, with persistent sulfadoxine-pyrimethamine resistance and no evidence of artemisinin resistance. Multiple single nucleotide polymorphisms were identified in csp, but their significance is uncertain. Our study demonstrates the feasibility of nanopore sequencing for malaria genomic surveillance in endemic countries.

Rapid profiling of Plasmodium parasites from genome sequences to assist malaria control

BACKGROUND

Malaria continues to be a major threat to global public health. Whole genome sequencing (WGS) of the underlying Plasmodium parasites has provided insights into the genomic epidemiology of malaria. Genome sequencing is rapidly gaining traction as a diagnostic and surveillance tool for clinical settings, where the profiling of co-infections, identification of imported malaria parasites, and detection of drug resistance are crucial for infection control and disease elimination. To support this informatically, we have developed the Malaria-Profiler tool, which rapidly (within minutes) predicts Plasmodium species, geographical source, and resistance to antimalarial drugs directly from WGS data.

RESULTS

The online and command line versions of Malaria-Profiler detect ~ 250 markers from genome sequences covering Plasmodium speciation, likely geographical source, and resistance to chloroquine, sulfadoxine-pyrimethamine (SP), and other anti-malarial drugs for P. falciparum, but also providing mutations for orthologous resistance genes in other species. The predictive performance of the mutation library was assessed using 9321 clinical isolates with WGS and geographical data, with most being single-species infections (P. falciparum 7152/7462, P. vivax 1502/1661, P. knowlesi 143/151, P. malariae 18/18, P. ovale ssp. 5/5), but co-infections were identified (456/9321; 4.8%). The accuracy of the predicted geographical profiles was high to both continental (96.1%) and regional levels (94.6%). For P. falciparum, markers were identified for resistance to chloroquine (49.2%; regional range: 24.5% to 100%), sulfadoxine (83.3%; 35.4- 90.5%), pyrimethamine (85.4%; 80.0-100%) and combined SP (77.4%). Markers associated with the partial resistance of artemisinin were found in WGS from isolates sourced from Southeast Asia (30.6%).

CONCLUSIONS

Malaria-Profiler is a user-friendly tool that can rapidly and accurately predict the geographical regional source and anti-malarial drug resistance profiles across large numbers of samples with WGS data. The software is flexible with modifiable bioinformatic pipelines. For example, it is possible to select the sequencing platform, display specific variants, and customise the format of outputs. With the increasing application of next-generation sequencing platforms on Plasmodium DNA, Malaria-Profiler has the potential to be integrated into point-of-care and surveillance settings, thereby assisting malaria control. Malaria-Profiler is available online (bioinformatics.lshtm.ac.uk/malaria-profiler) and as standalone software ( https://github.com/jodyphelan/malaria-profiler ).

A paradoxical population structure of var DBLα types in Africa

The var multigene family encodes the P. falciparum erythrocyte membrane protein 1 (PfEMP1), which is important in host-parasite interaction as a virulence factor and major surface antigen of the blood stages of the parasite, responsible for maintaining chronic infection. Whilst important in the biology of P. falciparum, these genes (50 to 60 genes per parasite genome) are routinely excluded from whole genome analyses due to their hyper-diversity, achieved primarily through recombination. The PfEMP1 head structure almost always consists of a DBLα-CIDR tandem. Categorised into different groups (upsA, upsB, upsC), different head structures have been associated with different ligand-binding affinities and disease severities. We study how conserved individual DBLα types are at the country, regional, and local scales in Sub-Saharan Africa. Using publicly-available sequence datasets and a novel ups classification algorithm, cUps, we performed an in silico exploration of DBLα conservation through time and space in Africa. In all three ups groups, the population structure of DBLα types in Africa consists of variants occurring at rare, low, moderate, and high frequencies. Non-rare variants were found to be temporally stable in a local area in endemic Ghana. When inspected across different geographical scales, we report different levels of conservation; while some DBLα types were consistently found in high frequencies in multiple African countries, others were conserved only locally, signifying local preservation of specific types. Underlying this population pattern is the composition of DBLα types within each isolate DBLα repertoire, revealed to also consist of a mix of types found at rare, low, moderate, and high frequencies in the population. We further discuss the adaptive forces and balancing selection, including host genetic factors, potentially shaping the evolution and diversity of DBLα types in Africa.

Whole-genome surveillance identifies markers of Plasmodium falciparum drug resistance and novel genomic regions under selection in Mozambique

IMPORTANCE

Malaria is a devastating disease caused by Plasmodium parasites. The evolution of parasite drug resistance continues to hamper progress toward malaria elimination, and despite extensive efforts to control malaria, it remains a leading cause of death in Mozambique and other countries in the region. The development of successful vaccines and identification of molecular markers to track drug efficacy are essential for managing the disease burden. We present an analysis of the parasite genome in Mozambique, a country with one of the highest malaria burdens globally and limited available genomic data, revealing current selection pressure. We contribute additional evidence to limited prior studies supporting the effectiveness of SWGA in producing reliable genomic data from complex clinical samples. Our results provide the identity of genomic loci that may be associated with current antimalarial drug use, including artemisinin and lumefantrine, and reveal selection pressure predicted to compromise the efficacy of current vaccine candidates.

MalariaSED: a deep learning framework to decipher the regulatory contributions of noncoding variants in malaria parasites

Malaria remains one of the deadliest infectious diseases. Transcriptional regulation effects of noncoding variants in this unusual genome of malaria parasites remain elusive. We developed a sequence-based, ab initio deep learning framework, MalariaSED, for predicting chromatin profiles in malaria parasites. The MalariaSED performance was validated by published ChIP-qPCR and TF motifs results. Applying MalariaSED to ~ 1.3 million variants shows that geographically differentiated noncoding variants are associated with parasite invasion and drug resistance. Further analysis reveals chromatin accessibility changes at Plasmodium falciparum rings are partly associated with artemisinin resistance. MalariaSED illuminates the potential functional roles of noncoding variants in malaria parasites.

Evolution of Partial Resistance to Artemisinins in Malaria Parasites in Uganda

BACKGROUND

Partial resistance of Plasmodium falciparum to the artemisinin component of artemisinin-based combination therapies, the most important malaria drugs, emerged in Southeast Asia and now threatens East Africa. Partial resistance, which manifests as delayed clearance after therapy, is mediated principally by mutations in the kelch protein K13 (PfK13). Limited longitudinal data are available on the emergence and spread of artemisinin resistance in Africa.

METHODS

We performed annual surveillance among patients who presented with uncomplicated malaria at 10 to 16 sites across Uganda from 2016 through 2022. We sequenced the gene encoding kelch 13 (pfk13) and analyzed relatedness using molecular methods. We assessed malaria metrics longitudinally in eight Ugandan districts from 2014 through 2021.

RESULTS

By 2021-2022, the prevalence of parasites with validated or candidate resistance markers reached more than 20% in 11 of the 16 districts where surveillance was conducted. The PfK13 469Y and 675V mutations were seen in far northern Uganda in 2016-2017 and increased and spread thereafter, reaching a combined prevalence of 10 to 54% across much of northern Uganda, with spread to other regions. The 469F mutation reached a prevalence of 38 to 40% in one district in southwestern Uganda in 2021-2022. The 561H mutation, previously described in Rwanda, was first seen in southwestern Uganda in 2021, reaching a prevalence of 23% by 2022. The 441L mutation reached a prevalence of 12 to 23% in three districts in western Uganda in 2022. Genetic analysis indicated local emergence of mutant parasites independent of those in Southeast Asia. The emergence of resistance was observed predominantly in areas where effective malaria control had been discontinued or transmission was unstable.

CONCLUSIONS

Data from Uganda showed the emergence of partial resistance to artemisinins in multiple geographic locations, with increasing prevalence and regional spread over time. (Funded by the National Institutes of Health.).

Malaria Genomics, Vaccine Development, and Microbiome

Recent advances in malaria genetics and genomics have transformed many aspects of malaria research in areas of molecular evolution, epidemiology, transmission, host-parasite interaction, drug resistance, pathogenicity, and vaccine development. Here, in addition to introducing some background information on malaria parasite biology, parasite genetics/genomics, and genotyping methods, we discuss some applications of genetic and genomic approaches in vaccine development and in studying interactions with microbiota. Genetic and genomic data can be used to search for novel vaccine targets, design an effective vaccine strategy, identify protective antigens in a whole-organism vaccine, and evaluate the efficacy of a vaccine. Microbiota has been shown to influence disease outcomes and vaccine efficacy; studying the effects of microbiota in pathogenicity and immunity may provide information for disease control. Malaria genetics and genomics will continue to contribute greatly to many fields of malaria research.

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.

An optimized GATK4 pipeline for Plasmodium falciparum whole genome sequencing variant calling and analysis

BACKGROUND

Accurate variant calls from whole genome sequencing (WGS) of Plasmodium falciparum infections are crucial in malaria population genomics. Here a falciparum variant calling pipeline based on GATK version 4 (GATK4) was optimized and applied to 6626 public Illumina WGS samples.

METHODS

Control WGS and accurate PacBio assemblies of 10 laboratory strains were leveraged to optimize parameters that control the heterozygosity, local assembly region size, ploidy, mapping and base quality in both GATK HaplotypeCaller and GenotypeGVCFs. From these controls, a high-quality training dataset was generated to recalibrate the raw variant data.

RESULTS

On current high-quality samples (read length = 250 bp, insert size = 405-524 bp), the optimized pipeline shows improved sensitivity (86.6 ± 1.7% for SNPs and 82.2 ± 5.9% for indels) compared to the default GATK4 pipeline (77.7 ± 1.3% for SNPs; and 73.1 ± 5.1% for indels, adjusted P < 0.001) and previous variant calling with GATK version 3 (GATK3, 70.3 ± 3.0% for SNPs and 59.7 ± 5.8% for indels, adjusted P < 0.001). Its sensitivity on simulated mixed infection samples (80.8 ± 6.1% for SNPs and 78.3 ± 5.1% for indels) was again improved relative to default GATK4 (68.8 ± 6.0% for SNPs and 38.9 ± 0.7% for indels, adjusted, adjusted P < 0.001). Precision was high and comparable across all pipelines on each type of data tested. The resulting combination of high-quality SNPs and indels increases the resolution of local population population structure detection in sub-Saharan Africa. Finally, increasing ploidy improves the detection of drug resistance mutations and estimation of complexity of infection.

CONCLUSIONS

Overall, this study provides an optimized falciparum GATK4 pipeline resource for variant calling which should help improve genomic studies of malaria.

Chloroquine resistance evolution in Plasmodium falciparum is mediated by the putative amino acid transporter AAT1

Malaria parasites break down host haemoglobin into peptides and amino acids in the digestive vacuole for export to the parasite cytoplasm for growth: interrupting this process is central to the mode of action of several antimalarial drugs. Mutations in the chloroquine (CQ) resistance transporter, pfcrt, located in the digestive vacuole membrane, confer CQ resistance in Plasmodium falciparum, and typically also affect parasite fitness. However, the role of other parasite loci in the evolution of CQ resistance is unclear. Here we use a combination of population genomics, genetic crosses and gene editing to demonstrate that a second vacuolar transporter plays a key role in both resistance and compensatory evolution. Longitudinal genomic analyses of the Gambian parasites revealed temporal signatures of selection on a putative amino acid transporter (pfaat1) variant S258L, which increased from 0% to 97% in frequency between 1984 and 2014 in parallel with the pfcrt1 K76T variant. Parasite genetic crosses then identified a chromosome 6 quantitative trait locus containing pfaat1 that is selected by CQ treatment. Gene editing demonstrated that pfaat1 S258L potentiates CQ resistance but at a cost of reduced fitness, while pfaat1 F313S, a common southeast Asian polymorphism, reduces CQ resistance while restoring fitness. Our analyses reveal hidden complexity in CQ resistance evolution, suggesting that pfaat1 may underlie regional differences in the dynamics of resistance evolution, and modulate parasite resistance or fitness by manipulating the balance between both amino acid and drug transport.

Multiplexed ddPCR-amplicon sequencing reveals isolated Plasmodium falciparum populations amenable to local elimination in Zanzibar, Tanzania

Zanzibar has made significant progress toward malaria elimination, but recent stagnation requires novel approaches. We developed a highly multiplexed droplet digital PCR (ddPCR)-based amplicon sequencing method targeting 35 microhaplotypes and drug-resistance loci, and successfully sequenced 290 samples from five districts covering both main islands. Here, we elucidate fine-scale Plasmodium falciparum population structure and infer relatedness and connectivity of infections using an identity-by-descent (IBD) approach. Despite high genetic diversity, we observe pronounced fine-scale spatial and temporal parasite genetic structure. Clusters of near-clonal infections on Pemba indicate persistent local transmission with limited parasite importation, presenting an opportunity for local elimination efforts. Furthermore, we observe an admixed parasite population on Unguja and detect a substantial fraction (2.9%) of significantly related infection pairs between Zanzibar and the mainland, suggesting recent importation. Our study provides a high-resolution view of parasite genetic structure across the Zanzibar archipelago and provides actionable insights for prioritizing malaria elimination efforts.

Targeted and whole-genome sequencing reveal a north-south divide in P. falciparum drug resistance markers and genetic structure in Mozambique

Mozambique is one of the four African countries which account for over half of all malaria deaths worldwide, yet little is known about the parasite genetic structure in that country. We performed P. falciparum amplicon and whole genome sequencing on 2251 malaria-infected blood samples collected in 2015 and 2018 in seven provinces of Mozambique to genotype antimalarial resistance markers and interrogate parasite population structure using genome-wide microhaplotyes. Here we show that the only resistance-associated markers observed at frequencies above 5% were pfmdr1-184F (59%), pfdhfr-51I/59 R/108 N (99%) and pfdhps-437G/540E (89%). The frequency of pfdhfr/pfdhps quintuple mutants associated with sulfadoxine-pyrimethamine resistance increased from 80% in 2015 to 89% in 2018 (p < 0.001), with a lower expected heterozygosity and higher relatedness of microhaplotypes surrounding pfdhps mutants than wild-type parasites suggestive of recent selection. pfdhfr/pfdhps quintuple mutants also increased from 72% in the north to 95% in the south (2018; p < 0.001). This resistance gradient was accompanied by a concentration of mutations at pfdhps-436 (17%) in the north, a south-to-north increase in the genetic complexity of P. falciparum infections (p = 0.001) and a microhaplotype signature of regional differentiation. The parasite population structure identified here offers insights to guide antimalarial interventions and epidemiological surveys.

coiaf: Directly estimating complexity of infection with allele frequencies

In malaria, individuals are often infected with different parasite strains. The complexity of infection (COI) is defined as the number of genetically distinct parasite strains in an individual. Changes in the mean COI in a population have been shown to be informative of changes in transmission intensity with a number of probabilistic likelihood and Bayesian models now developed to estimate the COI. However, rapid, direct measures based on heterozygosity or FwS do not properly represent the COI. In this work, we present two new methods that use easily calculated measures to directly estimate the COI from allele frequency data. Using a simulation framework, we show that our methods are computationally efficient and comparably accurate to current approaches in the literature. Through a sensitivity analysis, we characterize how the distribution of parasite densities, the assumed sequencing depth, and the number of sampled loci impact the bias and accuracy of our two methods. Using our developed methods, we further estimate the COI globally from Plasmodium falciparum sequencing data and compare the results against the literature. We show significant differences in the estimated COI globally between continents and a weak relationship between malaria prevalence and COI.

Generation of a mutator parasite to drive resistome discovery in Plasmodium falciparum

In vitro evolution of drug resistance is a powerful approach for identifying antimalarial targets, however, key obstacles to eliciting resistance are the parasite inoculum size and mutation rate. Here we sought to increase parasite genetic diversity to potentiate resistance selections by editing catalytic residues of Plasmodium falciparum DNA polymerase δ. Mutation accumulation assays reveal a ~5-8 fold elevation in the mutation rate, with an increase of 13-28 fold in drug-pressured lines. Upon challenge with the spiroindolone PfATP4-inhibitor KAE609, high-level resistance is obtained more rapidly and at lower inocula than wild-type parasites. Selections also yield mutants with resistance to an "irresistible" compound, MMV665794 that failed to yield resistance with other strains. We validate mutations in a previously uncharacterised gene, PF3D7_1359900, which we term quinoxaline resistance protein (QRP1), as causal for resistance to MMV665794 and a panel of quinoxaline analogues. The increased genetic repertoire available to this "mutator" parasite can be leveraged to drive P. falciparum resistome discovery.

Malaria outbreak in Laos driven by a selective sweep for Plasmodium falciparum kelch13 R539T mutants: a genetic epidemiology analysis

BACKGROUND

Malaria outbreaks are important public health concerns that can cause resurgence in endemic regions approaching elimination. We investigated a Plasmodium falciparum outbreak in Attapeu Province, Laos, during the 2020-21 malaria season, using genomic epidemiology methods to elucidate parasite population dynamics and identify its causes.

METHODS

In this genetic analysis, 2164 P falciparum dried blood spot samples were collected from southern Laos between Jan 1, 2017, and April 1, 2021, which included 249 collected during the Attapeu outbreak between April 1, 2020, and April 1, 2021, by routine surveillance. Genetic barcodes obtained from these samples were used to investigate epidemiological changes underpinning the outbreak, estimate population diversity, and analyse population structure. Whole-genome sequencing data from additional historical samples were used to reconstruct the ancestry of outbreak strains using identity-by-descent analyses.

FINDINGS

The outbreak parasite populations were characterised by unprecedented loss of genetic diversity, primarily caused by rapid clonal expansion of a multidrug-resistant strain (LAA1) carrying the kelch13 Arg539Thr (R539T) mutation. LAA1 replaced kelch13 Cys580Tyr (C580Y) mutants resistant to dihydroartemisinin-piperaquine (KEL1/PLA1) as the dominant strain. LAA1 inherited 58·8% of its genome from a strain circulating in Cambodia in 2008. A secondary outbreak strain (LAA2) carried the kelch13 C580Y allele, and a genome that is essentially identical to a Cambodian parasite from 2009. A third, low-frequency strain (LAA7) was a recombinant of KEL1/PLA1 with a kelch13 R539T mutant.

INTERPRETATION

These results strongly suggest that the outbreak was driven by a selective sweep, possibly associated with multidrug-resistant phenotypes of the outbreak strains. Established resistant populations can circulate at low frequencies for years before suddenly overwhelming dominant strains when the conditions for selection become favourable-eg, when front-line therapies change. Genetic surveillance can support elimination by characterising key properties of outbreaks such as population diversity, drug resistance marker prevalence, and the origins of outbreak strains.

FUNDING

Bill & Melinda Gates Foundation; The Global Fund to Fight AIDS, Tuberculosis and Malaria; Wellcome Trust.

TRANSLATION

For the Lao translation of the abstract see Supplementary Materials section.

Genomic variation during culture adaptation of genetically complex Plasmodium falciparum clinical isolates

Experimental studies on the biology of malaria parasites have mostly been based on laboratory-adapted lines, but there is limited understanding of how these may differ from parasites in natural infections. Loss-of-function mutants have previously been shown to emerge during culture of some Plasmodium falciparum clinical isolates in analyses focusing on single-genotype infections. The present study included a broader array of isolates, mostly representing multiple-genotype infections, which are more typical in areas where malaria is highly endemic. Genome sequence data from multiple time points over several months of culture adaptation of 28 West African isolates were analysed, including previously available sequences along with new genome sequences from additional isolates and time points. Some genetically complex isolates eventually became fixed over time to single surviving genotypes in culture, whereas others retained diversity, although proportions of genotypes varied over time. Drug resistance allele frequencies did not show overall directional changes, suggesting that resistance-associated costs are not the main causes of fitness differences among parasites in culture. Loss-of-function mutants emerged during culture in several of the multiple-genotype isolates, affecting genes (including AP2-HS, EPAC and SRPK1) for which loss-of-function mutants were previously seen to emerge in single-genotype isolates. Parasite clones were derived by limiting dilution from six of the isolates, and sequencing identified de novo variants not detected in the bulk isolate sequences. Interestingly, several of these were nonsense mutants and frameshifts disrupting the coding sequence of EPAC, the gene with the largest number of independent nonsense mutants previously identified in laboratory-adapted lines. Analysis of genomic identity by descent to explore relatedness among clones revealed co-occurring non-identical sibling parasites, illustrative of the natural genetic structure within endemic populations.

A snapshot of the prevalence of dihydropteroate synthase-431V mutation and other sulfadoxine-pyrimethamine resistance markers in Plasmodium falciparum isolates in Nigeria

BACKGROUND

Malaria is a major public health issue with substantial risks among vulnerable populations. Currently, the World Health Organization (WHO) recommends SP-IPTp in the second and third trimesters. However, the efficacy of SP-IPTp is threatened by the emergence of sulfadoxine-pyrimethamine resistant malaria parasites due to single nucleotide polymorphisms in the Plasmodium falciparum dihydrofolate reductase and dihydropteroate synthetase genes. This study aimed to assess the current prevalence of Pfdhfr/Pfdhps mutations in P. falciparum isolates collected from individuals residing in Ile-Ife, Nigeria, and also present maps of the prevalence of Pfdhps 431V and 581G within Nigeria and surrounding countries.

METHODS

Between October 2020 and April 2021, samples were collected as dried blood spots among 188 participants who showed malaria positivity with a histidine-rich-protein-based rapid diagnostic test (RDT). Nested PCR assays were used to confirm falciparum in the samples with RDT positivity, and to amplify fragments of the Pfdhfr/Pfdhps genes followed by targeted amplicon sequencing. Published data since 2007 on the prevalence of the Pfdhps genotypes in Nigeria and the neighbouring countries were used to produce maps to show the distribution of the mutant genotypes.

RESULTS

Only 74 and 61 samples were successfully amplified for the Pfdhfr and Pfdhps genes, respectively. At codons resulting in N51I, C59R, and S108N, Pfdhfr carried mutant alleles of 97.3% (72/74), 97.3% (72/74) and 98.6% (73/74), respectively. The Pfdhps gene carried mutations at codons resulting in amino acid changes at 431-436-437-540-581-613; I431V [45.9%, (28/61)], A581G [31.1% (19/61)] and A613S [49.2% (30/61)]. Constructed haplotypes were mainly the triple Pfdhfr mutant 51I-59R-108N (95.9%), and the most common haplotypes observed for the Pfdhps gene were the ISGKAA (32.8%), ISGKGS (8.2%), VAGKAA (14.8%), VAGKAS (9.8%) and VAGKGS (14.8%). In the context of the previously published data, a high prevalence of 431V/581G mutations was found in the study population. It seems quite evident that the Pfdhps 431V, 581G and 613S often co-occur as Pfdhps-VAGKGS haplotype.

CONCLUSION

This study showed that the prevalence of VAGKGS haplotype seems to be increasing in prevalence. If this is similar in effect to the emergence of 581G in East Africa, the efficacy of SP-IPTp in the presence of these novel Pfdhps mutants should be re-assessed.

Nosocomial Malaria Transmissions Resolved by Genomic Analyses-A Retrospective Case Report Study in France: 2007-2021

BACKGROUND

Exposure of blood to malaria parasites can lead to infection even in the absence of the mosquito vector. During a stay in a healthcare facility, accidental inoculation of the skin with blood from a malaria patient might occur, referred to as nosocomial malaria.

METHODS

Between 2007 and 2021, we identified 6 autochthonous malaria cases that occurred in different French hospitals, originating from nosocomial transmission and imported malaria cases being the infection source. Four cases were observed during the coronavirus disease 2019 pandemic. The genetic relatedness between source and nosocomial infections was evaluated by genome-wide short tandem repeats (STRs) and single-nucleotide polymorphisms (SNPs).

RESULTS

None of the patients with autochthonous malaria had travel history to an endemic area nor had been transfused. For each case, both the source and recipient patients stayed a few hours in the same ward. After diagnosis, autochthonous cases were treated with antimalarials and all recovered except 1. Genetically, each pair of matched source/nosocomial parasite infections showed <1% of different STRs and <6.9% (<1.5% for monoclonal infections) of different SNPs. Similar levels of genetic differences were obtained for parasite DNA samples that were independently sequenced twice as references of identical infections. Parasite phylogenomics were consistent with travel information reported by the source patients.

CONCLUSIONS

Our study demonstrates that genomics analyses may resolve nosocomial malaria transmissions, despite the uncertainty regarding the modes of contamination. Nosocomial transmission of potentially life-threatening parasites should be taken into consideration in settings or occasions where compliance with universal precautions is not rigorous.

Potent transmission-blocking monoclonal antibodies from naturally exposed individuals target a conserved epitope on Plasmodium falciparum Pfs230

Pfs230 is essential for Plasmodium falciparum transmission to mosquitoes and is the protein targeted by the most advanced malaria-transmission-blocking vaccine candidate. Prior understanding of functional epitopes on Pfs230 is based on two monoclonal antibodies (mAbs) with moderate transmission-reducing activity (TRA), elicited from subunit immunization. Here, we screened the B cell repertoire of two naturally exposed individuals possessing serum TRA and identified five potent mAbs from sixteen Pfs230 domain-1-specific mAbs. Structures of three potent and three low-activity antibodies bound to Pfs230 domain 1 revealed four distinct epitopes. Highly potent mAbs from natural infection recognized a common conformational epitope that is highly conserved across P. falciparum field isolates, while antibodies with negligible TRA derived from natural infection or immunization recognized three distinct sites. Our study provides molecular blueprints describing P. falciparum TRA, informed by contrasting potent and non-functional epitopes elicited by natural exposure and vaccination.

Highly potent, naturally acquired human monoclonal antibodies against Pfs48/45 block Plasmodium falciparum transmission to mosquitoes

Malaria transmission-blocking vaccines (TBVs) aim to induce antibodies that interrupt malaria parasite development in the mosquito, thereby blocking onward transmission, and provide a much-needed tool for malaria control and elimination. The parasite surface protein Pfs48/45 is a leading TBV candidate. Here, we isolated and characterized a panel of 81 human Pfs48/45-specific monoclonal antibodies (mAbs) from donors naturally exposed to Plasmodium parasites. Genetically diverse mAbs against each of the three domains (D1-D3) of Pfs48/45 were identified. The most potent mAbs targeted D1 and D3 and achieved >80% transmission-reducing activity in standard membrane-feeding assays, at 10 and 2 μg/mL, respectively. Co-crystal structures of D3 in complex with four different mAbs delineated two conserved protective epitopes. Altogether, these Pfs48/45-specific human mAbs provide important insight into protective and non-protective epitopes that can further our understanding of transmission and inform the design of refined malaria transmission-blocking vaccine candidates.

An Optimized GATK4 Pipeline for Plasmodium falciparum Whole Genome Sequencing Variant Calling and Analysis

Background Accurate variant calls from whole genome sequencing (WGS) of Plasmodium falciparum infections are crucial in malaria population genomics. Here we optimized a falciparum variant calling pipeline based on GATK version 4 (GATK4) and applied it to 6,626 public Illumina WGS samples. Methods We optimized parameters that control the heterozygosity, local assembly region size, ploidy, mapping and base quality in both GATK HaplotypeCaller and GenotypeGVCFs leveraging control WGS and accurate PacBio assemblies of 10 laboratory strains. From these controls we generated a high-quality training dataset to recalibrate the raw variant data. Results On current high-quality samples (read length = 250bp, insert size = 405 - 524 bp ), we show improved sensitivity (86.6 ± 1.7% for SNPs and 82.2 ± 5.9% for indels) compared to the default GATK4 pipeline (77.7 ± 1.3% for SNPs; and 73.1 ± 5.1% for indels, adjusted P < 0.001) and previous variant calling with GATK version 3 (GATK3, 70.3 ± 3.0% for SNPs and 59.7 ± 5.8% for indels, adjusted P < 0.001). The sensitivity of our pipeline on simulated mixed infection samples (80.8 ± 6.1% for SNPs and 78.3 ± 5.1% for indels) was again improved relative to default GATK4 (68.8 ± 6.0% for SNPs and 38.9 ± 0.7% for indels, adjusted P < 0.001). Precision was high and comparable across all pipelines on each type of data tested. We further show that using the combination of high-quality SNPs and indels increases the resolution of local population population structure detection in sub-Saharan Africa. We finally demonstrate that increasing ploidy improves the detection of drug resistance mutations and estimation of complexity of infection. Conclusions Overall, we provide an optimized GATK4 pipeline and resource for falciparum variant calling which should help improve genomic studies of malaria.

Unravelling var complexity: Relationship between DBLα types and var genes in Plasmodium falciparum

The enormous diversity and complexity of var genes that diversify rapidly by recombination has led to the exclusion of assembly of these genes from major genome initiatives (e.g., Pf6). A scalable solution in epidemiological surveillance of var genes is to use a small ‘tag’ region encoding the immunogenic DBLα domain as a marker to estimate var diversity. As var genes diversify by recombination, it is not clear the extent to which the same tag can appear in multiple var genes. This relationship between marker and gene has not been investigated in natural populations. Analyses of in vitro recombination within and between var genes have suggested that this relationship would not be exclusive. Using a dataset of publicly-available assembled var sequences, we test this hypothesis by studying DBLα-var relationships for four study sites in four countries: Pursat (Cambodia) and Mae Sot (Thailand), representing low malaria transmission, and Navrongo (Ghana) and Chikwawa (Malawi), representing high malaria transmission. In all study sites, DBLα-var relationships were shown to be predominantly 1-to-1, followed by a second largest proportion of 1-to-2 DBLα-var relationships. This finding indicates that DBLα tags can be used to estimate not just DBLα diversity but var gene diversity when applied in a local endemic area. Epidemiological applications of this result are discussed.

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.

Global Analysis of Plasmodium falciparum Dihydropteroate Synthase Variants Associated with Sulfadoxine Resistance Reveals Variant Distribution and Mechanisms of Resistance: A Computational-Based Study

The continual rise in sulfadoxine (SDX) resistance affects the therapeutic efficacy of sulfadoxine-pyrimethamine; therefore, careful monitoring will help guide its prolonged usage. Mutations in Plasmodium falciparum dihydropteroate synthase (Pfdhps) are being surveilled, based on their link with SDX resistance. However, there is a lack of continuous analyses and data on the potential effect of molecular markers on the Pfdhps structure and function. This study explored single-nucleotide polymorphisms (SNPs) in Pfdhps that were isolated in Africa and other countries, highlighting the regional distribution and its link with structure. In total, 6336 genomic sequences from 13 countries were subjected to SNPs, haplotypes, and structure-based analyses. The SNP analysis revealed that the key SDX resistance marker, A437G, was nearing fixation in all countries, peaking in Malawi. The mutation A613S was rare except in isolates from the Democratic Republic of Congo and Malawi. Molecular docking revealed a general loss of interactions when comparing mutant proteins to the wild-type protein. During MD simulations, SDX was released from the active site in mutants A581G and A613S before the end of run-time, whereas an unstable binding of SDX to mutant A613S and haplotype A437A/A581G/A613S was observed. Conformational changes in mutant A581G and the haplotypes A581G/A613S, A437G/A581G, and A437G/A581G/A613S were seen. The radius of gyration revealed an unfolding behavior for the A613S, K540E/A581G, and A437G/A581G systems. Overall, tracking such mutations by the continuous analysis of Pfdhps SNPs is encouraged. SNPs on the Pfdhps structure may cause protein-drug function loss, which could affect the applicability of SDX in preventing malaria in pregnant women and children.

Plasmodium falciparum Sexual Commitment Rate Variation among Clinical Isolates and Diverse Laboratory-Adapted Lines

Asexual blood-stage malaria parasites must produce sexual progeny to infect mosquitoes. It is important to understand the scope and causes of intraspecific variation in sexual commitment rates, particularly for the major human parasite P. falciparum. First, two alternative assay methods of measuring sexual commitment were compared to test a genetically modified P. falciparum line with elevated commitment rates inducible by overexpression of GDV1. The methods yielded correlated measurements with higher sensitivity and precision being achieved by one employing detection of the early gametocyte differentiation marker Pfs16. Thus, this was used to survey a diverse range of parasite lines and test each in multiple biological replicate assays in a serum-free medium supplemented with Albumax. There were differences among six recent clinical isolates from Ghana in their mean rates of sexual commitment per cycle, ranging from 3.3% to 12.2%. Among 13 diverse long-term laboratory-adapted lines, mean sexual commitment rates for most ranged from 4.7% to 13.4%, a few had lower rates with means from 0.3 to 1.6%, and one with a nonfunctional ap2-g gene always showed zero commitment. Among a subset of lines tested for the effects of exogenous choline to suppress commitment, there were significant differences. As expected, there was no effect in a line that had lost the gdv1 gene and that had generally low commitment, whereas the others showed quantitatively variable but significant responses to choline, suggesting potential trait variation. The results indicated the value of performing multiple replicate assays for understanding the variation of this key reproductive trait that likely affects transmission. IMPORTANCE Only sexual-stage malaria parasites are transmitted from human blood to mosquitoes. Thus, it is vital to understand variations in sexual commitment rates because these may be modifiable or susceptible to blocking. Two different methods of commitment rate measurement were first compared, demonstrating higher sensitivity and precision by the detection of an early differentiation marker, which was subsequently used to survey diverse lines. Clinical isolates from Ghana showed significant variation in mean per-cycle commitment rates and variation among biological replicates. Laboratory-adapted lines of diverse origins had a wider range with most being within the range observed for the clinical isolates, while a minority consistently had lower or zero rates. There was quantitative variation in the effects when adding choline to suppress commitment, indicating differing responsiveness of parasites to this environmental modification. Performing multiple assay replicates and comparisons of diverse isolates was important to understand this trait and its potential effects on transmission.

Mutant PfCRT Can Mediate Piperaquine Resistance in African Plasmodium falciparum With Reduced Fitness and Increased Susceptibility to Other Antimalarials

BACKGROUND

Additional therapeutic strategies could benefit efforts to reverse the recent increase in malaria cases in sub-Saharan Africa, which mostly affects young children. A primary candidate is dihydroartemisinin + piperaquine (DHA + PPQ), which is effective for uncomplicated malaria treatment, seasonal malaria chemoprevention, and intermittent preventive treatment. In Southeast Asia, Plasmodium falciparum parasites acquired PPQ resistance, mediated primarily by mutations in the P falciparum chloroquine resistance transporter PfCRT. The recent emergence in Africa of DHA-resistant parasites creates an imperative to assess whether PPQ resistance could emerge in African parasites with distinct PfCRT isoforms.

METHODS

We edited 2 PfCRT mutations known to mediate high-grade PPQ resistance in Southeast Asia into GB4 parasites from Gabon. Gene-edited clones were profiled in antimalarial concentration-response and fitness assays.

RESULTS

The PfCRT F145I mutation mediated moderate PPQ resistance in GB4 parasites but with a substantial fitness cost. No resistance was observed with the PfCRT G353V mutant. Both edited clones became significantly more susceptible to amodiaquine, chloroquine, and quinine.

CONCLUSIONS

A single PfCRT mutation can mediate PPQ resistance in GB4 parasites, but with a growth defect that may preclude its spread without further genetic adaptations. Our findings support regional use of drug combinations that exert opposing selective pressures on PfCRT.

R H: a genetic metric for measuring intrahost Plasmodium falciparum relatedness and distinguishing cotransmission from superinfection

Multiple-strain (polygenomic) infections are a ubiquitous feature of Plasmodium falciparum parasite population genetics. Under simple assumptions of superinfection, polygenomic infections are hypothesized to be the result of multiple infectious bites. As a result, polygenomic infections have been used as evidence of repeat exposure and used to derive genetic metrics associated with high transmission intensity. However, not all polygenomic infections are the result of multiple infectious bites. Some result from the transmission of multiple, genetically related strains during a single infectious bite (cotransmission). Superinfection and cotransmission represent two distinct transmission processes, and distinguishing between the two could improve inferences regarding parasite transmission intensity. Here, we describe a new metric, R H, that utilizes the correlation in allelic state (heterozygosity) within polygenomic infections to estimate the likelihood that the observed complexity resulted from either superinfection or cotransmission. R H is flexible and can be applied to any type of genetic data. As a proof of concept, we used R H to quantify polygenomic relatedness and estimate cotransmission and superinfection rates from a set of 1,758 malaria infections genotyped with a 24 single nucleotide polymorphism (SNP) molecular barcode. Contrary to expectation, we found that cotransmission was responsible for a significant fraction of 43% to 53% of the polygenomic infections collected in three distinct epidemiological regions in Senegal. The prediction that polygenomic infections frequently result from cotransmission stresses the need to incorporate estimates of relatedness within polygenomic infections to ensure the accuracy of genomic epidemiology surveillance data for informing public health activities.

Genome-wide analysis of Schistosoma mansoni reveals limited population structure and possible praziquantel drug selection pressure within Ugandan hot-spot communities

Populations within schistosomiasis control areas, especially those in Africa, are recommended to receive regular mass drug administration (MDA) with praziquantel (PZQ) as the main strategy for controlling the disease. The impact of PZQ treatment on schistosome genetics remains poorly understood, and is limited by a lack of high-resolution genetic data on the population structure of parasites within these control areas. We generated whole-genome sequence data from 174 individual miracidia collected from both children and adults from fishing communities on islands in Lake Victoria in Uganda that had received either annual or quarterly MDA with PZQ over four years, including samples collected immediately before and four weeks after treatment. Genome variation within and between samples was characterised and we investigated genomic signatures of natural selection acting on these populations that could be due to PZQ treatment. The parasite population on these islands was more diverse than found in nearby villages on the lake shore. We saw little or no genetic differentiation between villages, or between the groups of villages with different treatment intensity, but slightly higher genetic diversity within the pre-treatment compared to post-treatment parasite populations. We identified classes of genes significantly enriched within regions of the genome with evidence of recent positive selection among post-treatment and intensively treated parasite populations. The differential selection observed in post-treatment and pre-treatment parasite populations could be linked to any reduced susceptibility of parasites to praziquantel treatment.

Declines in prevalence alter the optimal level of sexual investment for the malaria parasite Plasmodium falciparum

Like most human pathogens, the malaria parasite Plasmodium falciparum experiences strong selection pressure from public health interventions such as drug treatment. While most commonly studied in the context of drug targets and related pathways, parasite adaptation to control measures likely extends to phenotypes beyond drug resistance. Here, we use modeling to explore how control measures can reduce levels of within-host competition between P. falciparum genotypes and favor higher rates of sexual investment. We validate these predictions with longitudinally sampled genomic data from French Guiana during a period of malaria decline and find that the most strongly selected genes are enriched for transcription factors involved in commitment to and development of the parasite’s sexual gametocyte form.

Successful infectious disease interventions can result in large reductions in parasite prevalence. Such demographic change has fitness implications for individual parasites and may shift the parasite’s optimal life history strategy. Here, we explore whether declining infection rates can alter Plasmodium falciparum’s investment in sexual versus asexual growth. Using a multiscale mathematical model, we demonstrate how the proportion of polyclonal infections, which decreases as parasite prevalence declines, affects the optimal sexual development strategy: Within-host competition in multiclone infections favors a greater investment in asexual growth whereas single-clone infections benefit from higher conversion to sexual forms. At the same time, drug treatment also imposes selection pressure on sexual development by shortening infection length and reducing within-host competition. We assess these models using 148 P. falciparum parasite genomes sampled in French Guiana over an 18-y period of intensive intervention (1998 to 2015). During this time frame, multiple public health measures, including the introduction of new drugs and expanded rapid diagnostic testing, were implemented, reducing P. falciparum malaria cases by an order of magnitude. Consistent with this prevalence decline, we see an increase in the relatedness among parasites, but no single clonal background grew to dominate the population. Analyzing individual allele frequency trajectories, we identify genes that likely experienced selective sweeps. Supporting our model predictions, genes showing the strongest signatures of selection include transcription factors involved in the development of P. falciparum’s sexual gametocyte form. These results highlight how public health interventions impose wide-ranging selection pressures that affect basic parasite life history traits.

Screening strategies and laboratory assays to support Plasmodium falciparum histidine-rich protein deletion surveillance: where we are and what is needed

Rapid diagnostic tests (RDTs) detecting Plasmodium falciparum histidine-rich protein 2 (HRP2) have been an important tool for malaria diagnosis, especially in resource-limited settings lacking quality microscopy. Plasmodium falciparum parasites with deletion of the pfhrp2 gene encoding this antigen have now been identified in dozens of countries across Asia, Africa, and South America, with new reports revealing a high prevalence of deletions in some selected regions. To determine whether HRP2-based RDTs are appropriate for continued use in a locality, focused surveys and/or surveillance activities of the endemic P. falciparum population are needed. Various survey and laboratory methods have been used to determine parasite HRP2 phenotype and pfhrp2 genotype, and the data collected by these different methods need to be interpreted in the appropriate context of survey and assay utilized. Expression of the HRP2 antigen can be evaluated using point-of-care RDTs or laboratory-based immunoassays, but confirmation of a deletion (or mutation) of pfhrp2 requires more intensive laboratory molecular assays, and new tools and strategies for rigorous but practical data collection are particularly needed for large surveys. Because malaria diagnostic strategies are typically developed at the national level, nationally representative surveys and/or surveillance that encompass broad geographical areas and large populations may be required. Here is discussed contemporary assays for the phenotypic and genotypic evaluation of P. falciparum HRP2 status, consider their strengths and weaknesses, and highlight key concepts relevant to timely and resource-conscious workflows required for efficient diagnostic policy decision making.

How has mass drug administration with dihydroartemisinin-piperaquine impacted molecular markers of drug resistance? A systematic review

The World Health Organization (WHO) recommends surveillance of molecular markers of resistance to anti-malarial drugs. This is particularly important in the case of mass drug administration (MDA), which is endorsed by the WHO in some settings to combat malaria. Dihydroartemisinin-piperaquine (DHA-PPQ) is an artemisinin-based combination therapy which has been used in MDA. This review analyses the impact of MDA with DHA-PPQ on the evolution of molecular markers of drug resistance. The review is split into two parts. Section I reviews the current evidence for different molecular markers of resistance to DHA-PPQ. This includes an overview of the prevalence of these molecular markers in Plasmodium falciparum Whole Genome Sequence data from the MalariaGEN Pf3k project. Section II is a systematic literature review of the impact that MDA with DHA-PPQ has had on the evolution of molecular markers of resistance. This systematic review followed PRISMA guidelines. This review found that despite being a recognised surveillance tool by the WHO, the surveillance of molecular markers of resistance following MDA with DHA-PPQ was not commonly performed. Of the total 96 papers screened for eligibility in this review, only 20 analysed molecular markers of drug resistance. The molecular markers published were also not standardized. Overall, this warrants greater reporting of molecular marker prevalence following MDA implementation. This should include putative pfcrt mutations which have been found to convey resistance to DHA-PPQ in vitro.

Two fatal autochthonous cases of airport malaria, Belgium, 2020

We report an outbreak investigation of two fatal cases of autochthonous Plasmodium falciparum malaria that occurred in Belgium in September 2020. Various hypotheses of the potential source of infection were investigated. The most likely route of transmission was through an infectious exotic Anopheles mosquito that was imported via the international airport of Brussels or the military airport Melsbroek and infected the cases who lived at 5 km from the airports. Based on genomic analysis of the parasites collected from the two cases, the most likely origin of the Plasmodium was Gabon or Cameroon. Further, the parasites collected from the two Belgian patients were identical by descent, which supports the assumption that the two infections originated from the bite of the same mosquito, during interrupted feeding. Although airport malaria remains a rare event, it has significant implications, particularly for the patient, as delayed or missed diagnosis of the cause of illness often results in complications and mortality. Therefore, to prevent such severe or fatal outcomes, we suggest a number of public health actions including increased awareness among health practitioners, especially those working in the vicinity of airports, and increased surveillance of exotic mosquito species at airports.

Circulation of an artemisinin-resistant malaria lineage in a traveller returning from East Africa to France

A returned traveller to Uganda presented with a Plasmodium falciparum kelch13 A675V mutant infection that exhibited delayed clearance under artesunate therapy. Parasites were genetically related to recently reported Ugandan artemisinin-resistant A675V parasites. Adequate malaria prevention measures and clinical and genotypic surveillance are important tools to avoid and track artemisinin resistance.

Malaria protection due to sickle haemoglobin depends on parasite genotype

Host genetic factors can confer resistance against malaria¹, raising the question of whether this has led to evolutionary adaptation of parasite populations. Here we searched for association between candidate host and parasite genetic variants in 3,346 Gambian and Kenyan children with severe malaria caused by Plasmodium falciparum. We identified a strong association between sickle haemoglobin (HbS) in the host and three regions of the parasite genome, which is not explained by population structure or other covariates, and which is replicated in additional samples. The HbS-associated alleles include nonsynonymous variants in the gene for the acyl-CoA synthetase family member^2–4PfACS8 on chromosome 2, in a second region of chromosome 2, and in a region containing structural variation on chromosome 11. The alleles are in strong linkage disequilibrium and have frequencies that covary with the frequency of HbS across populations, in particular being much more common in Africa than other parts of the world. The estimated protective effect of HbS against severe malaria, as determined by comparison of cases with population controls, varies greatly according to the parasite genotype at these three loci. These findings open up a new avenue of enquiry into the biological and epidemiological significance of the HbS-associated polymorphisms in the parasite genome and the evolutionary forces that have led to their high frequency and strong linkage disequilibrium in African P. falciparum populations.

A strong association has been found between three regions of the Plasmodium falciparum genome and sickle haemoglobin in children with severe malaria, suggesting parasites have adapted to overcome natural host immunity.

Global diversity of the gene encoding the Pfs25 protein-a Plasmodium falciparum transmission-blocking vaccine candidate

Background

Vaccines against the sexual stages of the malarial parasite Plasmodium falciparum are indispensable for controlling malaria and abrogating the spread of drug-resistant parasites. Pfs25, a surface antigen of the sexual stage of P. falciparum, is a leading candidate for transmission-blocking vaccine development. While clinical trials have reported that Pfs25-based vaccines are safe and effective in inducing transmission-blocking antibodies, the extent of the genetic diversity of Pfs25 in malaria endemic populations has rarely been studied. Thus, this study aimed to investigate the global diversity of Pfs25 in P. falciparum populations.

Methods

A database of 307 Pfs25 sequences of P. falciparum was established. Population genetic analyses were performed to evaluate haplotype and nucleotide diversity, analyze haplotypic distribution patterns of Pfs25 in different geographical populations, and construct a haplotype network. Neutrality tests were conducted to determine evidence of natural selection. Homology models of the Pfs25 haplotypes were constructed, subjected to molecular dynamics (MD), and analyzed in terms of flexibility and percentages of secondary structures.

Results

The Pfs25 gene of P. falciparum was found to have 11 unique haplotypes. Of these, haplotype 1 (H1) and H2, the major haplotypes, represented 70% and 22% of the population, respectively, and were dominant in Asia, whereas only H1 was dominant in Africa, Central America, and South America. Other haplotypes were rare and region-specific, resulting in unique distribution patterns in different geographical populations. The diversity in Pfs25 originated from ten single-nucleotide polymorphism (SNP) loci located in the epidermal growth factor (EGF)-like domains and anchor domain. Of these, an SNP at position 392 (GGA/GCA), resulting in amino acid substitution 131 (Gly/Ala), defined the two major haplotypes. The MD results showed that the structures of H1 and H2 variants were relatively similar. Limited polymorphism in Pfs25 could likely be due to negative selection.

Conclusions

The study successfully established a Pfs25 sequence database that can become an essential tool for monitoring vaccine efficacy, designing assays for detecting malaria carriers, and conducting epidemiological studies of P. falciparum. The discovery of the two major haplotypes, H1 and H2, and their conserved structures suggests that the current Pfs25-based vaccines could be used globally for malaria control.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-05078-6.

A population genetic perspective on the origin, spread and adaptation of the human malaria agents Plasmodium falciparum and Plasmodium vivax

Malaria is considered one of the most important scourges that humanity has faced during its history, being responsible every year for numerous deaths worldwide. The disease is caused by protozoan parasites, among which two species are responsible of the majority of the burden, Plasmodium falciparum and Plasmodium vivax. For these two parasite species, the questions of their origin (how and when they appeared in humans), of their spread throughout the world, as well as how they have adapted to humans have long been of interest to the scientific community. Here, we review the current knowledge that has accumulated on these different questions, thanks in particular to the analysis of the genetic and genomic variability of these parasites and comparison with related Plasmodium species infecting other host species (like non-human primates). In this paper we review the existing body of knowledge, including current research dealing with these questions, focusing particularly on genetic analysis and genomic variability of these parasites and comparison with related Plasmodium species infecting other species of host (such as non-human primates).

A novel CRISPR-based malaria diagnostic capable of Plasmodium detection, species differentiation, and drug-resistance genotyping

Background

CRISPR-based diagnostics are a new class of highly sensitive and specific assays with multiple applications in infectious disease diagnosis. SHERLOCK, or Specific High-Sensitivity Enzymatic Reporter UnLOCKing, is one such CRISPR-based diagnostic that combines recombinase polymerase pre-amplification, CRISPR-RNA base-pairing, and LwCas13a activity for nucleic acid detection.

Methods

We developed SHERLOCK assays capable of detecting all Plasmodium species known to cause human malaria and species-specific detection of P. vivax and P. falciparum, the species responsible for the majority of malaria cases worldwide. We further tested these assays using a diverse panel of clinical samples from the Democratic Republic of the Congo, Uganda, and Thailand and pools of Anopheles mosquitoes from Thailand. In addition, we developed a prototype SHERLOCK assay capable of detecting the dihydropteroate synthetase (dhps) single nucleotide variant A581G associated with P. falciparum sulfadoxine resistance.

Findings

The suite of Plasmodium assays achieved analytical sensitivities ranging from 2•5-18•8 parasites per reaction when tested against laboratory strain genomic DNA. When compared to real-time PCR, the P. falciparum assay achieved 94% sensitivity and 94% specificity during testing of 123 clinical samples. Compared to amplicon-based deep sequencing, the dhps SHERLOCK assay achieved 73% sensitivity and 100% specificity when applied to a panel of 43 clinical samples, with false-negative calls only at lower parasite densities.

Interpretation

These novel SHERLOCK assays demonstrate the versatility of CRISPR-based diagnostics and their potential as a new generation of molecular tools for malaria diagnosis and surveillance.

Funding

National Institutes of Health (T32GM007092, R21AI148579, K24AI134990, R01AI121558, UL1TR002489, P30CA016086)