Major subpopulations of Plasmodium falciparum in sub-Saharan Africa

Systematic bias in malaria parasite relatedness estimation

Genetic studies of Plasmodium parasites increasingly feature relatedness estimates. However, various aspects of malaria parasite relatedness estimation are not fully understood. For example, relatedness estimates based on whole-genome-sequence (WGS) data often exceed those based on sparser data types. Systematic bias in relatedness estimation is well documented in the literature geared towards diploid organisms, but largely unknown within the malaria community. We characterize systematic bias in malaria parasite relatedness estimation using three complementary approaches: theoretically, under a non-ancestral statistical model of pairwise relatedness; numerically, under a simulation model of ancestry; and empirically, using data on parasites sampled from Guyana and Colombia. We show that allele frequency estimates encode, locus-by-locus, relatedness averaged over the set of sampled parasites used to compute them. Plugging sample allele frequencies into models of pairwise relatedness can lead to systematic underestimation. However, systematic underestimation can be viewed as population-relatedness calibration, i.e., a way of generating measures of relative relatedness. Systematic underestimation is unavoidable when relatedness is estimated assuming independence between genetic markers. It is mitigated when relatedness is estimated using WGS data under a hidden Markov model (HMM) that exploits linkage between proximal markers. The extent of mitigation is unknowable when a HMM is fit to sparser data, but downstream analyses that use high relatedness thresholds are relatively robust regardless. In summary, practitioners can either resolve to use relative relatedness estimated under independence, or try to estimate absolute relatedness under a HMM. We propose various tools to help practitioners evaluate their situation on a case-by-case basis.

Potential and pitfalls of using identity-by-descent for malaria genomic surveillance

The ability to genotype malaria parasites on an epidemiological scale is crucial for genomic surveillance as it aids in understanding malaria transmission dynamics and parasite demography changes in response to antimalarial interventions. Identity-by-descent (IBD)-based methods have demonstrated potential in various aspects of malaria genomic surveillance. However, there is a need for validation of existing approaches and development of new techniques to address challenges posed by the parasites' unique evolutionary dynamics and complex biological characteristics, which differ markedly from organisms like humans. This review examines current IBD use cases, identifies limitations of IBD-based methods, and explores promising future directions to enhance malaria genomic surveillance.

Reduced Genetic Diversity of Key Fertility and Vector Competency Related Genes in Anopheles gambiae s.l. Across Sub-Saharan Africa

BACKGROUND

Insecticide resistance challenges the vector control efforts towards malaria elimination and proving the development of complementary tools. Targeting the genes that are involved in mosquito fertility and susceptibility to Plasmodium with small molecule inhibitors has been a promising alternative to curb the vector population and drive the transmission down. However, such an approach would require a comprehensive knowledge of the genetic diversity of the targeted genes to ensure the broad efficacy of new tools across the natural vector populations.

METHODS

Four fertility and parasite susceptibility genes were identified from a systematic review of the literature. The Single Nucleotide Polymorphisms (SNPs) found within the regions spanned by these four genes, genotyped across 2784 wild-caught Anopheles gambiae s.l. from 19 sub-Saharan African (SSA) countries, were extracted from the whole genome SNP data of the Ag1000G project (Ag3.0). The population genetic analysis on gene-specific data included the determination of the population structure, estimation of the differentiation level between the populations, evaluation of the linkage between the non-synonymous SNPs (nsSNPs), and a few statistical tests.

RESULTS

As potential targets for small molecule inhibitors to reduce malaria transmission, our set of four genes associated with Anopheles fertility and their susceptibility to Plasmodium comprises the mating-induced stimulator of oogenesis protein (MISO, AGAP002620), Vitellogenin (Vg, AGAP004203), Lipophorin (Lp, AGAP001826), and Haem-peroxidase 15 (HPX15, AGAP013327). The analyses performed on these potential targets of small inhibitor molecules revealed that the genes are conserved within SSA populations of An. gambiae s.l. The overall low Fst values and low clustering of principal component analysis between species indicated low genetic differentiation at all the genes (MISO, Vg, Lp and HPX15). The low nucleotide diversity (>0.10), negative Tajima's D values, and heterozygosity analysis provided ecological insights into the purifying selection that acts to remove deleterious mutations, maintaining genetic diversity at low levels within the populations. None of MISO nsSNPs were identified in linkage disequilibrium, whereas a few weakly linked nsSNPs with ambiguous haplotyping were detected at other genes.

CONCLUSIONS

This integrated finding on the genetic features of major malaria vectors' biological factors across natural populations offer new insights for developing sustainable malaria control tools. These loci were reasonably conserved, allowing for the design of effective targeting with small molecule inhibitors towards controlling vector populations and lowering global malaria transmission.

Distinct evolutionary regimes across domains of the Plasmodium falciparum CSP gene

Malaria disease caused by parasites of genus Plasmodium places an enormous disease burden across tropical regions of the world. The circumsporozoite protein (CSP) of Plasmodium has several key functions in binding and accessing host cells, with functions subdivided across multiple protein regions. While its key roles during infection make the gene a primary target for malaria vaccine development, the evolutionary dynamics that could affect the forecasting of useful strains remain poorly understood. We tested whether the gene undergoes multiple DNA substitution processes and whether these are divided across gene regions using a phylogenetic mixture model, and a global sample of CSP sequences specific to P. falciparum. These analyses reveal evolutionary processes unique to the central repeat region and the C-terminus. The central repeat region is dominated by synonymous substitutions (putatively neutral) and heavy C-T substitution bias, while the C-terminus undergoes mostly non-synonymous changes. These evolutionary processes are not strongly geographically restricted, and lineages from Africa and Asia where the parasite is most abundant appear to drive evolution across all CSP gene regions. We propose that insights about DNA substitution processes can help forecast the variants of importance to vaccine development, aided by state-of-the-art evolutionary modelling.

Identification of complex Plasmodium falciparum genetic backgrounds circulating in Africa: a multicountry genomic epidemiology analysis

BACKGROUND

The population structure of the malaria parasite Plasmodium falciparum can reveal underlying adaptive evolutionary processes. Selective pressures to maintain complex genetic backgrounds can encourage inbreeding, producing distinct parasite clusters identifiable by population structure analyses.

METHODS

We analysed population structure in 3783 P falciparum genomes from 21 countries across Africa, provided by the MalariaGEN Pf7 dataset. We used Principal Coordinate Analysis to cluster parasites, identity by descent (IBD) methods to identify genomic regions shared by cluster members, and linkage analyses to establish their co-inheritance patterns. Structural variants were reconstructed by de novo assembly and verified by long-read sequencing.

FINDINGS

We identified a strongly differentiated cluster of parasites, named AF1, comprising 47 (1·2%) of 3783 samples analysed, distributed over 13 countries across Africa, at locations over 7000 km apart. Members of this cluster share a complex genetic background, consisting of up to 23 loci harbouring many highly differentiated variants, rarely observed outside the cluster. IBD analyses revealed common ancestry at these loci, irrespective of sampling location. Outside the shared loci, however, AF1 members appear to outbreed with sympatric parasites. The AF1 differentiated variants comprise structural variations, including a gene conversion involving the dblmsp and dblmsp2 genes, and numerous single nucleotide polymorphisms. Several of the genes harbouring these mutations are functionally related, often involved in interactions with red blood cells including invasion, egress, and erythrocyte antigen export.

INTERPRETATION

We propose that AF1 parasites have adapted to some unidentified evolutionary niche, probably involving interactions with host erythrocytes. This adaptation involves a complex compendium of interacting variants that are rarely observed in Africa, which remains mostly intact despite recombination events. The term cryptotype was used to describe a common background interspersed with genomic regions of local origin.

FUNDING

Bill & Melinda Gates Foundation.

Systematic in vitro evolution in Plasmodium falciparum reveals key determinants of drug resistance

Surveillance of drug resistance and the discovery of novel targets-key objectives in the fight against malaria-rely on identifying resistance-conferring mutations in Plasmodium parasites. Current approaches, while successful, require laborious experimentation or large sample sizes. To elucidate shared determinants of antimalarial resistance that can empower in silico inference, we examined the genomes of 724 Plasmodium falciparum clones, each selected in vitro for resistance to one of 118 compounds. We identified 1448 variants in 128 recurrently mutated genes, including drivers of antimalarial multidrug resistance. In contrast to naturally occurring variants, those selected in vitro are more likely to be missense or frameshift, involve bulky substitutions, and occur in conserved, ordered protein domains. Collectively, our dataset reveals mutation features that predict drug resistance in eukaryotic pathogens.

Population genomics of Plasmodium ovale species in sub-Saharan Africa

Plasmodium ovale curtisi (Poc) and Plasmodium ovale wallikeri (Pow) are relapsing malaria parasites endemic to Africa and Asia that were previously thought to represent a single species. Amid increasing detection of ovale malaria in sub-Saharan Africa, we present a population genomic study of both species across the continent. We conducted whole-genome sequencing of 25 isolates from Central and East Africa and analyzed them alongside 20 previously published African genomes. Isolates are predominantly monoclonal (43/45), with their genetic similarity aligning with geography. Pow shows lower average nucleotide diversity (1.8×10-4) across the genome compared to Poc (3.0×10-4) (p < 0.0001). Signatures of selective sweeps involving the dihydrofolate reductase gene have been found in both species, as are signs of balancing selection at the merozoite surface protein 1 gene. Differences in the nucleotide diversity of Poc and Pow may reflect unique demographic history, even as similar selective forces facilitate their resilience to malaria control interventions.

Setting Up an NGS Sequencing Platform and Monitoring Molecular Markers of Anti-Malarial Drug Resistance in Djibouti

Djibouti is confronted with malaria resurgence, with malaria having been occurring in epidemic proportions since a decade ago. The current epidemiology of drug-resistant Plasmodium falciparum is not well known. Molecular markers were analyzed by targeted sequencing in 79 P. falciparum clinical isolates collected in Djibouti city in 2023 using the Miseq Illumina platform newly installed in the country. The objective of the study was to analyze the key codons in these molecular markers associated with antimalarial drug resistance. The prevalence of the mutant Pfcrt CVIET haplotype (92%) associated with chloroquine resistance and mutant Pfdhps-Pfdhfr haplotypes (7.4% SGEA and 53.5% IRN, respectively) associated with sulfadoxine-pyrimethamine resistance was high. By contrast, Pfmdr1 haplotypes associated with amodiaquine (YYY) or lumefantrine (NFD) resistance were not observed in any of the isolates. Although the "Asian-type" PfK13 mutations associated with artemisinin resistance were not observed, the "African-type" PfK13 substitution, R622I, was found in a single isolate (1.4%) for the first time in Djibouti. Our genotyping data suggest that most Djiboutian P. falciparum isolates are resistant to chloroquine and sulfadoxine-pyrimethamine but are sensitive to amodiaquine, lumefantrine, and artemisinin. Nonetheless, the presence of an isolate with the R622I PfK13 substitution is a warning signal that calls for a regular surveillance of molecular markers of antimalarial drug resistance.

Population genomics of Plasmodium ovale species in sub-Saharan Africa

Plasmodium ovale curtisi (Poc) and Plasmodium ovale wallikeri (Pow) are relapsing malaria parasites endemic to Africa and Asia that were previously thought to represent a single species. Amid increasing detection of ovale malaria in sub-Saharan Africa, we performed a population genomic study of both species across the continent. We conducted whole-genome sequencing of 25 isolates from Central and East Africa and analyzed them alongside 20 previously published African genomes. Isolates were predominantly monoclonal (43/45), with their genetic similarity aligning with geography. Pow showed lower average nucleotide diversity (1.8×10-4) across the genome compared to Poc (3.0×10-4) (p < 0.0001). Signatures of selective sweeps involving the dihydrofolate reductase gene were found in both species, as were signs of balancing selection at the merozoite surface protein 1 gene. Differences in the nucleotide diversity of Poc and Pow may reflect unique demographic history, even as similar selective forces facilitate their resilience to malaria control interventions.

A retrospective analysis of P. falciparum drug resistance markers detects an early (2016/17) high prevalence of the k13 C469Y mutation in asymptomatic infections in Northern Uganda

The emergence of drug-resistant Plasmodium falciparum parasites in sub-Saharan Africa will substantially challenge malaria control. Here, we evaluated the frequency of common drug resistance markers among adolescents from Northern Uganda with asymptomatic infections. We used an established amplicon deep sequencing strategy to screen dried blood spot samples collected from 2016 to 2017 during a reported malaria epidemic within the districts of Kitgum and Pader in Northern Uganda. We screened single-nucleotide polymorphisms within: kelch13 (Pfk13), dihydropteroate synthase (Pfdhps), multidrug resistance-1 (Pfmdr1), dihydrofolate reductase (Pfdhfr), and apical membrane antigen (Pfama1) genes. Within the study population, the median age was 15 years (14.3-15.0, 95% CI), and 54.9% (78/142) were Plasmodium positive by 18S rRNA qPCR, which were subsequently targeted for sequencing analysis. We observed a high frequency of resistance markers particularly for sulfadoxine-pyrimethamine (SP), with no wild-type-only parasites observed for Pfdhfr (N51I, C59R, and S108N) and Pfdhps (A437G and K540E) mutations. Within Pfmdr1, mixed infections were common for NF/NY (98.5%). While for artemisinin resistance, in kelch13, there was a high frequency of C469Y (34%). Using the pattern for Pfama1, we found a high level of polygenomic infections with all individuals presenting with complexity of infection greater than 2 with a median of 6.9. The high frequency of the quintuple SP drug-resistant parasites and the C469Y artemisinin resistance-associated mutation in asymptomatic individuals suggests an earlier high prevalence than previously reported from symptomatic malaria surveillance studies (in 2016/2017). Our data demonstrate the urgency for routine genomic surveillance programs throughout Africa and the value of deep sequencing.

National genomic profiling of Plasmodium falciparum antimalarial resistance in Zambian children participating in the 2018 Malaria Indicator Survey

The emergence of antimalarial drug resistance is a major threat to malaria control and elimination. Using whole genome sequencing of 282 P. falciparum samples collected during the 2018 Zambia National Malaria Indicator Survey, we determined the prevalence and spatial distribution of known and candidate antimalarial drug resistance mutations. High levels of genotypic resistance were found across Zambia to pyrimethamine, with over 94% (n=266) of samples having the Pfdhfr triple mutant (N51I, C59R, and S108N), and sulfadoxine, with over 84% (n=238) having the Pfdhps double mutant (A437G and K540E). In northern Zambia, 5.3% (n=15) of samples also harbored the Pfdhps A581G mutation. Although 29 mutations were identified in Pfkelch13, these mutations were present at low frequency (<2.5%), and only three were WHO-validated artemisinin partial resistance mutations: P441L (n=1, 0.35%), V568M (n=2, 0.7%) and R622T (n=1, 0.35%). Notably, 91 (32%) of samples carried the E431K mutation in the Pfatpase6 gene, which is associated with artemisinin resistance. No specimens carried any known mutations associated with chloroquine resistance in the Pfcrt gene (codons 72-76). P. falciparum strains circulating in Zambia were highly resistant to sulfadoxine and pyrimethamine but remained susceptible to chloroquine and artemisinin. Despite this encouraging finding, early genetic signs of developing artemisinin resistance highlight the urgent need for continued vigilance and expanded routine genomic surveillance to monitor these changes.

Plasmodium falciparum population dynamics in East Africa and genomic surveillance along the Kenya-Uganda border

East African countries accounted for ~ 10% of all malaria prevalence worldwide in 2022, with an estimated 23.8 million cases and > 53,000 deaths. Despite recent increases in malaria incidence, high-resolution genome-wide analyses of Plasmodium parasite populations are sparse in Kenya, Tanzania, and Uganda. The Kenyan-Ugandan border region is a particular concern, with Uganda confirming the emergence and spread of artemisinin resistant P. falciparum parasites. To establish genomic surveillance along the Kenyan-Ugandan border and analyse P. falciparum population dynamics within East Africa, we generated whole-genome sequencing (WGS) data for 38 parasites from Bungoma, Western Kenya. These sequences were integrated into a genomic analysis of available East African isolate data (n = 599) and revealed parasite subpopulations with distinct genetic structure and diverse ancestral origins. Ancestral admixture analysis of these subpopulations alongside isolates from across Africa (n = 365) suggested potential independent ancestral populations from other major African populations. Within isolates from Western Kenya, the prevalence of biomarkers associated with chloroquine resistance (e.g. Pfcrt K76T) were significantly reduced compared to wider East African populations and a single isolate contained the PfK13 V568I variant, potentially linked to reduced susceptibility to artemisinin. Overall, our work provides baseline WGS data and analysis for future malaria genomic surveillance in the region.

Transcriptomic analyses of differentially expressed human genes, micro RNAs and long-non-coding RNAs in severe, symptomatic and asymptomatic malaria infection

Malaria transmission and endemicity in Africa remains hugely disproportionate compared to the rest of the world. The complex life cycle of P. falciparum (Pf) between the vertebrate human host and the anopheline vector results in differential expression of genes within and between hosts. An in-depth understanding of Pf interaction with various human genes through regulatory elements will pave way for identification of newer tools in the arsenal for malaria control. Therefore, the regulatory elements (REs) involved in the over- or under-expression of various host immune genes hold the key to elucidating alternative control measures that can be applied for disease surveillance, prompt diagnosis and treatment. We carried out an RNAseq analysis to identify differentially expressed genes and network elucidation of non-coding RNAs and target genes associated with immune response in individuals with different clinical outcomes. Raw RNAseq datasets, retrieved for analyses include individuals with severe (Gambia-20), symptomatic (Burkina Faso-15), asymptomatic (Mali-16) malaria as well as uninfected controls (Tanzania-20; Mali-36). Of the total 107 datasets retrieved, we identified 5534 differentially expressed genes (DEGs) among disease and control groups. A peculiar pattern of DEGs was observed, with individuals presenting with severe/symptomatic malaria having the highest and most diverse upregulated genes, while a reverse phenomenon was recorded among asymptomatic and uninfected individuals. In addition, we identified 141 differentially expressed micro RNA (miRNA), of which 78 and 63 were upregulated and downregulated respectively. Interactome analysis revealed a moderate interaction between DEGs and miRNAs. Of all identified miRNA, five were unique (hsa-mir-32, hsa-mir-25, hsa-mir-221, hsa-mir-29 and hsa-mir-148) because of their connectivity to several genes, including hsa-mir-221 connected to 16 genes. Six-hundred and eight differentially expressed long non coding RNA (lncRNA) were also identified, including SLC7A11, LINC01524 among the upregulated ones. Our study provides important insight into host immune genes undergoing differential expression under different malaria conditions. It also identified unique miRNAs and lncRNAs that modify and/or regulate the expression of various immune genes. These regulatory elements we surmise, have the potential to serve a diagnostic purpose in discriminating between individuals with severe/symptomatic malaria and those with asymptomatic infection or uninfected, following further clinical validation from field isolates.

Benchmarking and Optimization of Methods for the Detection of Identity-By-Descent in High-Recombining Plasmodium falciparum Genomes

Genomic surveillance is crucial for identifying at-risk populations for targeted malaria control and elimination. Identity-by-descent (IBD) is increasingly being used in Plasmodium population genomics to estimate genetic relatedness, effective population size (N e ), population structure, and signals of positive selection. Despite its potential, a thorough evaluation of IBD segment detection tools for species with high recombination rates, such as P. falciparum, remains absent. Here, we perform comprehensive benchmarking of IBD callers - probabilistic (hmmIBD, isoRelate), identity-by-state-based (hap-IBD, phased IBD) and others (Refined IBD) - using population genetic simulations tailored for high recombination, and IBD quality metrics at both the IBD segment level and the IBD-based downstream inference level. Our results demonstrate that low marker density per genetic unit, related to high recombination relative to mutation, significantly compromises the accuracy of detected IBD segments. In genomes with high recombination rates resembling P. falciparum, most IBD callers exhibit high false negative rates for shorter IBD segments, which can be partially mitigated through optimization of IBD caller parameters, especially those related to marker density. Notably, IBD detected with optimized parameters allows for more accurate capture of selection signals and population structure; IBD-based N e inference is very sensitive to IBD detection errors, with IBD called from hmmIBD uniquely providing less biased estimates of N e in this context. Validation with empirical data from the MalariaGEN Pf 7 database, representing different transmission settings, corroborates these findings. We conclude that context-specific evaluation and parameter optimization are essential for accurate IBD detection in high-recombining species and recommend hmmIBD for quality-sensitive analysis, such as estimation of N e in these species. Our optimization and high-level benchmarking methods not only improve IBD segment detection in high-recombining genomes but also enhance overall genomic analysis, paving the way for more accurate genomic surveillance and targeted intervention strategies for malaria.

Plasmodium vivax populations in the western Greater Mekong Subregion evaluated using a genetic barcode

An improved understanding of the Plasmodium vivax populations in the Great Mekong Subregion (GMS) is needed to monitor the progress of malaria elimination. This study aimed to use a P. vivax single nucleotide polymorphism (SNP) barcode to evaluate the population dynamics and explore the gene flow among P. vivax parasite populations in the western GMS (China, Myanmar and Thailand). A total of 315 P. vivax patient samples collected in 2011 and 2018 from four regions of the western GMS were genotyped for 42 SNPs using the high-throughput MassARRAY SNP genotyping technology. Population genetic analysis was conducted to estimate the genetic diversity, effective population size, and population structure among the P. vivax populations. Overall, 291 samples were successfully genotyped at 39 SNPs. A significant difference was observed in the proportion of polyclonal infections among the five P. vivax populations (P = 0.0012, Pearson Chi-square test, χ2 = 18.1), with western Myanmar having the highest proportion (96.2%, 50/52) in 2018. Likewise, the average complexity of infection was also highest in western Myanmar (1.31) and lowest in northeast Myanmar (1.01) in 2018. The older samples from western China in 2011 had the highest pairwise nucleotide diversity (π, 0.388 ± 0.046), expected heterozygosity (He, 0.363 ± 0.02), and the largest effective population size. In comparison, in the neighboring northeast Myanmar, the more recent samples in 2018 showed the lowest values (π, 0.224 ± 0.036; He, 0.220 ± 0.026). Furthermore, the 2018 northeast Myanmar parasites showed high and moderate genetic differentiation from other populations with FST values of 0.162-0.252, whereas genetic differentiation among other populations was relatively low (FST ≤ 0.059). Principal component analysis, phylogeny, and STRUCTURE analysis showed that the P. vivax population in northeast Myanmar in 2018 substantially diverged from other populations. Although the 42 SNP barcode is a valuable tool for tracking parasite origins of worldwide parasite populations, a more extended barcode with additional SNPs is needed to distinguish the more related parasite populations in the western GMS.

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.

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.

Population genomic evidence of structured and connected Plasmodium vivax populations under host selection in Latin America

Pathogen genomic epidemiology has the potential to provide a deep understanding of population dynamics, facilitating strategic planning of interventions, monitoring their impact, and enabling timely responses, and thereby supporting control and elimination efforts of parasitic tropical diseases. Plasmodium vivax, responsible for most malaria cases outside Africa, shows high genetic diversity at the population level, driven by factors like sub-patent infections, a hidden reservoir of hypnozoites, and early transmission to mosquitoes. While Latin America has made significant progress in controlling Plasmodium falciparum, it faces challenges with residual P. vivax. To characterize genetic diversity and population structure and dynamics, we have analyzed the largest collection of P. vivax genomes to date, including 1474 high-quality genomes from 31 countries across Asia, Africa, Oceania, and America. While P. vivax shows high genetic diversity globally, Latin American isolates form a distinctive population, which is further divided into sub-populations and occasional clonal pockets. Genetic diversity within the continent was associated with the intensity of transmission. Population differentiation exists between Central America and the North Coast of South America, vs. the Amazon Basin, with significant gene flow within the Amazon Basin, but limited connectivity between the Northwest Coast and the Amazon Basin. Shared genomic regions in these parasite populations indicate adaptive evolution, particularly in genes related to DNA replication, RNA processing, invasion, and motility - crucial for the parasite's survival in diverse environments. Understanding these population-level adaptations is crucial for effective control efforts, offering insights into potential mechanisms behind drug resistance, immune evasion, and transmission dynamics.

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.

Genomics reveals heterogeneous Plasmodium falciparum transmission and population differentiation in Zambia and bordering countries

Genomic surveillance plays a critical role in monitoring malaria transmission and understanding how the parasite adapts in response to interventions. We conducted genomic surveillance of malaria by sequencing 241 Plasmodium falciparum genomes from regions with varying levels of malaria transmission across Zambia. We found genomic evidence of high levels of within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in the country. We identified country-level clustering of parasites from Zambia and neighboring countries, and distinct clustering of parasites from West Africa. Within Zambia, our identity by descent (IBD) relatedness analysis uncovered spatial clustering of closely related parasite pairs at the local level and rare cases of long-distance sharing. Genomic regions with large shared IBD segments and strong positive selection signatures identified genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Together, our findings enhance our understanding of P. falciparum transmission nationwide in Zambia and highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

Using a mobile nanopore sequencing lab for end-to-end genomic surveillance of Plasmodium falciparum: A feasibility study

Genomic epidemiology holds promise for malaria control and elimination efforts, for example by informing on Plasmodium falciparum genetic diversity and prevalence of mutations conferring anti-malarial drug resistance. Limited sequencing infrastructure in many malaria-endemic areas prevents the rapid generation of genomic data. To address these issues, we developed and validated assays for P. falciparum nanopore sequencing in endemic sites using a mobile laboratory, targeting key antimalarial drug resistance markers and microhaplotypes. Using two multiplexed PCR reactions, we amplified six highly polymorphic microhaplotypes and ten drug resistance markers. We developed a bioinformatics workflow that allows genotyping of polyclonal malaria infections, including minority clones. We validated the panels on mock dried blood spot (DBS) and rapid diagnostic test (RDT) samples and archived DBS, demonstrating even, high read coverage across amplicons (range: 580x to 3,212x median coverage), high haplotype calling accuracy, and the ability to explore within-sample diversity of polyclonal infections. We field-tested the feasibility of rapid genotyping in Zanzibar in close collaboration with the local malaria elimination program using DBS and routinely collected RDTs as sample inputs. Our assay identified haplotypes known to confer resistance to known antimalarials in the dhfr, dhps and mdr1 genes, but no evidence of artemisinin partial resistance. Most infections (60%) were polyclonal, with high microhaplotype diversity (median HE = 0.94). In conclusion, our assays generated actionable data within a few days, and we identified current challenges for implementing nanopore sequencing in endemic countries to accelerate malaria control and elimination.

Artemisinin resistance-associated gene mutations in Plasmodium falciparum: A case study of severe malaria from Mozambique

BACKGROUND

The effectiveness of artemisinin-based combination therapies (ACT) in treating Plasmodium falciparum, is vital for global malaria control efforts, particularly in sub-Saharan Africa. The examination of imported cases from endemic areas holds implications for malaria chemotherapy on a global scale.

METHOD

A 45-year-old male presented with high fever, dry cough, diarrhoea and generalized muscle pain, following a two-week trip to Mozambique. P. falciparum infection with hiperparasitemia was confirmed and the patient was treated initially with quinine and doxycycline, then intravenous artesunate. To assess drug susceptibility, ex vivo half-maximal inhibitory concentration assays were conducted, and the isolated P. falciparum genome was deep sequenced.

RESULTS

The clinical isolate exhibited elevated ex vivo half-maximal inhibitory concentration values to dihydroartemisinin, lumefantrine, mefloquine and piperaquine. Genomic analysis identified a I416V mutation in the P. falciparum Kelch13 (PF3D7_1343700) gene, and several mutations at the Kelch13 interaction candidate genes, pfkics (PF3D7_0813000, PF3D7_1138700, PF3D7_1246300), including the ubiquitin carboxyl-terminal hydrolase 1, pfubp1 (PF3D7_0104300). Mutations at the drug transporters and genes linked to next-generation antimalarial drug resistance were also present.

CONCLUSIONS

This case highlights the emergence of P. falciparum strains carrying mutations in artemisinin resistance-associated genes in Mozambique, couple with a reduction in ex vivo susceptibility to ACT drugs. Continuous surveillance of mutations linked to drug resistance and regular monitoring of drug susceptibility are imperative to anticipate the spread of potential resistant strains emerging in Mozambique and to maintain effective malaria control strategies.

Mathematical modeling of malaria transmission dynamics in humans with mobility and control states

Malaria importation is one of the hypothetical drivers of malaria transmission dynamics across the globe. Several studies on malaria importation focused on the effect of the use of conventional malaria control strategies as approved by the World Health Organization (WHO) on malaria transmission dynamics but did not capture the effect of the use of traditional malaria control strategies by vigilant humans. In order to handle the aforementioned situation, a novel system of Ordinary Differential Equations (ODEs) was developed comprising the human and the malaria vector compartments. Analysis of the system was carried out to assess its quantitative properties. The novel computational algorithm used to solve the developed system of ODEs was implemented and benchmarked with the existing Runge-Kutta numerical solution method. Furthermore, simulations of different vigilant conditions useful to control malaria were carried out. The novel system of malaria models was well-posed and epidemiologically meaningful based on its quantitative properties. The novel algorithm performed relatively better in terms of model simulation accuracy than Runge-Kutta. At the best model-fit condition of 98% vigilance to the use of conventional and traditional malaria control strategies, this study revealed that malaria importation has a persistent impact on malaria transmission dynamics. In lieu of this, this study opined that total vigilance to the use of the WHO-approved and traditional malaria management tools would be the most effective control strategy against malaria importation.

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.

Mapping the genomic landscape of multidrug resistance in Plasmodium falciparum and its impact on parasite fitness

Drug-resistant Plasmodium falciparum parasites have swept across Southeast Asia and now threaten Africa. By implementing a P. falciparum genetic cross using humanized mice, we report the identification of key determinants of resistance to artemisinin (ART) and piperaquine (PPQ) in the dominant Asian KEL1/PLA1 lineage. We mapped k13 as the central mediator of ART resistance in vitro and identified secondary markers. Applying bulk segregant analysis, quantitative trait loci mapping using 34 recombinant haplotypes, and gene editing, our data reveal an epistatic interaction between mutant PfCRT and multicopy plasmepsins 2/3 in mediating high-grade PPQ resistance. Susceptibility and parasite fitness assays implicate PPQ as a driver of selection for KEL1/PLA1 parasites. Mutant PfCRT enhanced susceptibility to lumefantrine, the first-line partner drug in Africa, highlighting a potential benefit of opposing selective pressures with this drug and PPQ. We also identified that the ABCI3 transporter can operate in concert with PfCRT and plasmepsins 2/3 in mediating multigenic resistance to antimalarial agents.

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.

Antimalarial drug discovery: progress and approaches

Recent antimalarial drug discovery has been a race to produce new medicines that overcome emerging drug resistance, whilst considering safety and improving dosing convenience. Discovery efforts have yielded a variety of new molecules, many with novel modes of action, and the most advanced are in late-stage clinical development. These discoveries have led to a deeper understanding of how antimalarial drugs act, the identification of a new generation of drug targets, and multiple structure-based chemistry initiatives. The limited pool of funding means it is vital to prioritize new drug candidates. They should exhibit high potency, a low propensity for resistance, a pharmacokinetic profile that favours infrequent dosing, low cost, preclinical results that demonstrate safety and tolerability in women and infants, and preferably the ability to block Plasmodium transmission to Anopheles mosquito vectors. In this Review, we describe the approaches that have been successful, progress in preclinical and clinical development, and existing challenges. We illustrate how antimalarial drug discovery can serve as a model for drug discovery in diseases of poverty.

A genomic platform for surveillance and antigen discovery in Plasmodium spp. using long-read amplicon sequencing

Many vaccine candidate proteins in the malaria parasite Plasmodium falciparum are under strong immunological pressure and confer antigenic diversity. We present a sequencing and data analysis platform for the genomic surveillance of the insertion or deletion (indel)-rich antigens merozoite surface protein 1 (MSP1), MSP2, glutamate-rich protein (GLURP), and CSP from P. falciparum using long-read circular consensus sequencing (CCS) in multiclonal malaria isolates. Our platform uses 40 PCR primers per gene to asymmetrically barcode and identify multiclonal infections in pools of up to 384 samples. With msp2, we validated the method using 235 mock infections combining 10 synthetic variants at different concentrations and infection complexities. We applied this strategy to P. falciparum isolates from a longitudinal cohort in Tanzania. Finally, we constructed an analysis pipeline that streamlines the processing and interpretation of epidemiological and antigenic diversity data from demultiplexed FASTQ files. This platform can be easily adapted to other polymorphic antigens of interest in Plasmodium or any other human pathogen.

AIDMAN: An AI-based object detection system for malaria diagnosis from smartphone thin-blood-smear images

Malaria is a significant public health concern, with ∼95% of cases occurring in Africa, but accurate and timely diagnosis is problematic in remote and low-income areas. Here, we developed an artificial intelligence-based object detection system for malaria diagnosis (AIDMAN). In this system, the YOLOv5 model is used to detect cells in a thin blood smear. An attentional aligner model (AAM) is then applied for cellular classification that consists of multi-scale features, a local context aligner, and multi-scale attention. Finally, a convolutional neural network classifier is applied for diagnosis using blood-smear images, reducing interference caused by false positive cells. The results demonstrate that AIDMAN handles interference well, with a diagnostic accuracy of 98.62% for cells and 97% for blood-smear images. The prospective clinical validation accuracy of 98.44% is comparable to that of microscopists. AIDMAN shows clinically acceptable detection of malaria parasites and could aid malaria diagnosis, especially in areas lacking experienced parasitologists and equipment.

Elucidation of the low-expressing erythroid CR1 phenotype by bioinformatic mining of the GATA1-driven blood-group regulome

Genetic determinants underlying most human blood groups are now clarified but variation in expression levels remains largely unexplored. By developing a bioinformatics pipeline analyzing GATA1/Chromatin immunoprecipitation followed by sequencing (ChIP-seq) datasets, we identify 193 potential regulatory sites in 33 blood-group genes. As proof-of-concept, we aimed to delineate the low-expressing complement receptor 1 (CR1) Helgeson phenotype on erythrocytes, which is correlated with several diseases and protects against severe malaria. We demonstrate that two candidate CR1 enhancer motifs in intron 4 bind GATA1 and drive transcription. Both are functionally abolished by naturally-occurring SNVs. Erythrocyte CR1-mRNA and CR1 levels correlate dose-dependently with genotype of one SNV (rs11117991) in two healthy donor cohorts. Haplotype analysis of rs11117991 with previously proposed markers for Helgeson shows high linkage disequilibrium in Europeans but explains the poor prediction reported for Africans. These data resolve the longstanding debate on the genetic basis of inherited low CR1 and form a systematic starting point to investigate the blood group regulome.

Population structure and genetic connectivity of Plasmodium falciparum in pre-elimination settings of Southern Africa

To accelerate malaria elimination in the Southern African region by 2030, it is essential to prevent cross-border malaria transmission. However, countries within the region are highly interconnected due to human migration that aids in the movement of the parasite across geographical borders. It is therefore important to better understand Plasmodium falciparum transmission dynamics in the region, and identify major parasite source and sink populations, as well as cross-border blocks of high parasite connectivity. We performed a meta-analysis using collated parasite allelic data generated by microsatellite genotyping of malaria parasites from Namibia, Eswatini, South Africa, and Mozambique (N = 5,314). The overall number of unique alleles was significantly higher (P ≤ 0.01) in Namibia (mean A = 17.3 ± 1.46) compared to South Africa (mean A = 12.2 ± 1.22) and Eswatini (mean A = 13.3 ± 1.27, P ≤ 0.05), whilst the level of heterozygosity was not significantly different between countries. The proportion of polyclonal infections was highest for Namibia (77%), and lowest for Mozambique (64%). A was significant population structure was detected between parasites from the four countries, and patterns of gene flow showed that Mozambique was the major source area and Eswatini the major sink area of parasites between the countries. This study showed strong signals of parasite population structure and genetic connectivity between malaria parasite populations across national borders. This calls for strengthening the harmonization of malaria control and elimination efforts between countries in the southern African region. This data also proves its potential utility as an additional surveillance tool for malaria surveillance on both a national and regional level for the identification of imported cases and/or outbreaks, as well as monitoring for the potential spread of anti-malarial drug resistance as countries work towards malaria elimination.

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.

Mapping the genomic landscape of multidrug resistance in Plasmodium falciparum and its impact on parasite fitness

Drug-resistant Plasmodium falciparum parasites have swept across Southeast Asia and now threaten Africa. By implementing a P. falciparum genetic cross using humanized mice, we report the identification of key determinants of resistance to artemisinin (ART) and piperaquine (PPQ) in the dominant Asian KEL1/PLA1 lineage. We mapped k13 as the central mediator of ART resistance and identified secondary markers. Applying bulk segregant analysis, quantitative trait loci mapping and gene editing, our data reveal an epistatic interaction between mutant PfCRT and multicopy plasmepsins 2/3 in mediating high-grade PPQ resistance. Susceptibility and parasite fitness assays implicate PPQ as a driver of selection for KEL1/PLA1 parasites. Mutant PfCRT enhanced susceptibility to lumefantrine, the first-line partner drug in Africa, highlighting a potential benefit of opposing selective pressures with this drug and PPQ. We also identified that the ABCI3 transporter can operate in concert with PfCRT and plasmepsins 2/3 in mediating multigenic resistance to antimalarial agents.

Performance of SNP barcodes to determine genetic diversity and population structure of Plasmodium falciparum in Africa

Panels of informative biallelic single nucleotide polymorphisms (SNPs) have been proposed to be an economical method to fast-track the population genetic analysis of Plasmodium falciparum in malaria-endemic areas. Whilst used successfully in low-transmission areas where infections are monoclonal and highly related, we present the first study to evaluate the performance of these 24- and 96-SNP molecular barcodes in African countries, characterised by moderate-to-high transmission, where multiclonal infections are prevalent. For SNP barcodes it is generally recommended that the SNPs chosen i) are biallelic, ii) have a minor allele frequency greater than 0.10, and iii) are independently segregating, to minimise bias in the analysis of genetic diversity and population structure. Further, to be standardised and used in many population genetic studies, these barcodes should maintain characteristics i) to iii) across various iv) geographies and v) time points. Using haplotypes generated from the MalariaGEN P. falciparum Community Project version six database, we investigated the ability of these two barcodes to fulfil these criteria in moderate-to-high transmission African populations in 25 sites across 10 countries. Predominantly clinical infections were analysed, with 52.3% found to be multiclonal, generating high proportions of mixed-allele calls (MACs) per isolate thereby impeding haplotype construction. Of the 24- and 96-SNPs, loci were removed if they were not biallelic and had low minor allele frequencies in all study populations, resulting in 20- and 75-SNP barcodes respectively for downstream population genetics analysis. Both SNP barcodes had low expected heterozygosity estimates in these African settings and consequently biased analyses of similarity. Both minor and major allele frequencies were temporally unstable. These SNP barcodes were also shown to identify weak genetic differentiation across large geographic distances based on Mantel Test and DAPC. These results demonstrate that these SNP barcodes are vulnerable to ascertainment bias and as such cannot be used as a standardised approach for malaria surveillance in moderate-to-high transmission areas in Africa, where the greatest genomic diversity of P. falciparum exists at local, regional and country levels.

Genetic variations of Plasmodium falciparum circumsporozoite protein and the impact on interactions with human immunoproteins and malaria vaccine efficacy

In October 2021, the world's first malaria vaccine RTS,S was endorsed by WHO for broad use in children, despite its low efficacy. This study examined polyclonal infections and the associations of parasite genetic variations with binding affinity to human leukocyte antigen (HLA). Multiplicity of infection was determined by amplicon deep sequencing of PfMSP1. Genetic variations in PfCSP were examined across 88 samples from Ghana and analyzed together with 1655 PfCSP sequences from other African and non-African isolates. Binding interactions of PfCSP peptide variants and HLA were predicted using NetChop and HADDOCK. High polyclonality was detected among infections, with each infection harboring multiple non-3D7 PfCSP variants. Twenty-seven PfCSP haplotypes were detected in the Ghanaian samples, and they broadly represented PfCSP diversity across Africa. The number of genetic differences between 3D7 and non-3D7 PfCSP variants does not influence binding to HLA. However, CSP peptide length after proteolytic degradation significantly affects its molecular weight and binding affinity to HLA. Despite the high diversity of HLA, the majority of the HLAI and II alleles interacted/bound with all Ghana CSP peptides. Multiple non-3D7 strains among P. falciparum infections could impact the effectiveness of RTS,S. Longer peptides of the Th2R/Th3R CSP regions should be considered in future versions of RTS,S.

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.

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.

Plasmodium malariae structure and genetic diversity in sub-Saharan Africa determined from microsatellite variants and linked SNPs in orthologues of antimalarial resistance genes

Plasmodium malariae, a neglected human malaria parasite, contributes up to 10% of malaria infections in sub-Saharan Africa (sSA). Though P. malariae infection is considered clinically benign, it presents mostly as coinfections with the dominant P. falciparum. Completion of its reference genome has paved the way to further understand its biology and interactions with the human host, including responses to antimalarial interventions. We characterized 75 P. malariae isolates from seven endemic countries in sSA using highly divergent microsatellites. The P. malariae infections were highly diverse and five subpopulations from three ancestries (independent of origin of isolates) were determined. Sequences of 11 orthologous antimalarial resistance genes, identified low frequency single nucleotide polymorphisms (SNPs), strong linkage disequilibrium between loci that may be due to antimalarial drug selection. At least three sub-populations were detectable from a subset of denoised SNP data from mostly the mitochondrial cytochrome b coding region. This evidence of diversity and selection calls for including P. malariae in malaria genomic surveillance towards improved tools and strategies for malaria elimination.

Low genetic diversity of Plasmodium falciparum merozoite surface protein 1 and 2 and multiplicity of infections in western Ethiopia following effective malaria interventions

BACKGROUND

Genetic diversity of malaria parasites can inform the intensity of transmission and poses a major threat to malaria control and elimination interventions. Characterization of the genetic diversity would provide essential information about the ongoing control efforts. This study aimed to explore allelic polymorphism of merozoite surface protein 1 (msp1) and merozoite surface protein 2 (msp2) to determine the genetic diversity and multiplicity of Plasmodium falciparum infections circulating in high and low transmission sites in western Ethiopia.

METHODS

Parasite genomic DNA was extracted from a total of 225 dried blood spots collected from confirmed uncomplicated P. falciparum malaria-infected patients in western Ethiopia. Of these, 72.4% (163/225) and 27.6% (62/225) of the samples were collected in high and low transmission areas, respectively. Polymorphic msp1 and msp2 genes were used to explore the genetic diversity and multiplicity of falciparum malaria infections. Genotyping of msp1 was successful in 86.5% (141/163) and 88.7% (55/62) samples collected from high and low transmission areas, respectively. Genotyping of msp2 was carried out among 85.3% (139/163) and 96.8% (60/62) of the samples collected in high and low transmission sites, respectively. Plasmodium falciparum msp1 and msp2 genes were amplified by nested PCR and the PCR products were analysed by QIAxcel ScreenGel Software. A P-value of less or equal to 0.05 was considered significant.

RESULTS

High prevalence of falciparum malaria was identified in children less than 15 years as compared with those  ≥ 15 years old (AOR = 2.438, P = 0.005). The three allelic families of msp1 (K1, MAD20, and RO33) and the two allelic families of msp2 (FC27 and 3D7), were observed in samples collected in high and low transmission areas. However, MAD 20 and FC 27 alleles were the predominant allelic families in both settings. Plasmodium falciparum isolates circulating in western Ethiopia had low genetic diversity and mean MOI. No difference in mean MOI between high transmission sites (mean MOI 1.104) compared with low transmission area (mean MOI 1.08) (p > 0.05). The expected heterozygosity of msp1 was slightly higher in isolates collected from high transmission sites (He = 0.17) than in those isolates from low transmission (He = 0.12). However, the heterozygosity of msp2 was not different in both settings (Pfmsp2: 0.04 in high transmission; pfmsp2: 0.03 in low transmission).

CONCLUSION

Plasmodium falciparum from clinical malaria cases in western Ethiopia has low genetic diversity and multiplicity of infection irrespective of the intensity of transmission at the site of sampling. These may be signaling the effectiveness of malaria control strategies in Ethiopia; although further studies are required to determine how specific intervention strategies and other parameters that drive the pattern.

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

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

Plasmodium falciparum merozoite invasion ligands, linked antimalarial resistance loci and ex vivo responses to antimalarials in The Gambia

BACKGROUND

Artemether/lumefantrine is the most commonly used artemisinin-based combination treatment (ACT) for malaria in sub-Saharan Africa. Drug resistance to ACT components is a major threat to malaria elimination efforts. Therefore, rigorous monitoring of drug efficacy is required for adequate management of malaria and to sustain the effectiveness of ACTs.

OBJECTIVES

This study identified and described genomic loci that correlate with differences in ex vivo responses of natural Plasmodium falciparum isolates from The Gambia to antimalarial drugs.

METHODS

Natural P. falciparum isolates from The Gambia were assayed for IC50 responses to four antimalarial drugs (artemether, dihydroartemisinin, amodiaquine and lumefantrine). Genome-wide SNPs from 56 of these P. falciparum isolates were applied to mixed-model regression and network analyses to determine linked loci correlating with drug responses. Genomic regions of shared haplotypes and positive selection within and between Gambian and Cambodian P. falciparum isolates were mapped by identity-by-descent (IBD) analysis of 209 genomes.

RESULTS

SNPs in 71 genes, mostly involved in stress and drug resistance mechanisms correlated with drug responses. Additionally, erythrocyte invasion and permeability loci, including merozoite surface proteins (Pfdblmsp, Pfsurfin), and high-molecular-weight rhoptry protein 2 (Pfrhops2) were correlated with responses to multiple drugs. Haplotypes of pfdblmsp2 and known drug resistance loci (pfaat1, pfcrt and pfdhfr) from The Gambia showed high IBD with those from Cambodia, indicating co-ancestry, with significant linkage disequilibrium between their alleles.

CONCLUSIONS

Multiple linked genic loci correlating with drug response phenotypes suggest a genomic backbone may be under selection by antimalarials. This calls for further analysis of molecular pathways to drug resistance in African P. falciparum.

Genomic approaches for monitoring transmission dynamics of malaria: A case for malaria molecular surveillance in Sub-Saharan Africa

Transmission dynamics is an important indicator for malaria control and elimination. As we move closer to eliminating malaria in Sub-Saharan Africa (sSA), transmission indices with higher resolution (genomic approaches) will complement our current measurements of transmission. Most of the present programmatic knowledge of malaria transmission patterns are derived from assessments of epidemiologic and clinical data, such as case counts, parasitological estimates of parasite prevalence, and Entomological Inoculation Rates (EIR). However, to eliminate malaria from endemic areas, we need to track changes in the parasite population and how they will impact transmission. This is made possible through the evolving field of genomics and genetics, as well as the development of tools for more in-depth studies on the diversity of parasites and the complexity of infections, among other topics. If malaria elimination is to be achieved globally, country-specific elimination activities should be supported by parasite genomic data from regularly collected blood samples for diagnosis, surveillance and possibly from other programmatic interventions. This presents a unique opportunity to track the spread of malaria parasites and shed additional light on intervention efficacy. In this review, various genetic techniques are highlighted along with their significance for an enhanced understanding of transmission patterns in distinct topological settings throughout Sub-Saharan Africa. The importance of these methods and their limitations in malaria surveillance to guide control and elimination strategies, are explored.

Anopheles gambiae s.l. swarms trapping as a complementary tool against residual malaria transmission in eastern Gambia

Malaria remains a major health problem and vector control is an essential approach to decrease its burden, although it is threatened by insecticide resistance. New approaches for vector control are needed. The females of Anopheles gambiae s.l. mate once in their life and in the swarms formed by males. Trapping swarms of Anopheles gambiae s.l. males is a potential new intervention for vector control, alternative to the use of insecticides, as it would disrupt mating . The proof-of-concept pilot study aiming at investigating swarm trapping as a potential vector control intervention, was carried out in 6 villages as in eastern Gambia. Swarms of Anopheles gambiae s.l. were identified and their size, height, and duration determined during the baseline year. Swarm trapping by local volunteers was implemented the following transmission season in 4 villages while the other 2 villages were taken as controls. Entomological outcomes were monitored by Human Landing Catches and Pyrethrum Spray Catches. A cross-sectional survey to determine malaria prevalence was carried out at the peak of the malaria transmission season for two consecutive years. At baseline, 23 swarming sites of Anopheles gambiae s.l. were identified. Before the intervention, mean indoor resting density per house and malaria prevalence were similar between control and intervention villages. Following the intervention, Anopheles gambiae s.l. indoor resting density was 44% lower in intervention than in control villages (adj IRR: 0.0.56; 95% CI 0.47-0.68); the odds of malaria infections were 68% lower in intervention than in control villages (OR: 0.32; 95% CI 0.11-0.97). Swarm trapping seems to be a promising, community-based vector control intervention that could reduce malaria prevalence by reducing vector density. Such results should be further investigated and confirmed by larger cluster-randomized trials.

Genome-wide SNP analysis of Plasmodium falciparum shows differentiation at drug-resistance-associated loci among malaria transmission settings in southern Mali

Plasmodium falciparum malaria cases in Africa represent over 90% of the global burden with Mali being amongst the 11 highest burden countries that account for 70% of this annual incidence. The persistence of P. falciparum despite massive global interventions is because of its genetic diversity that drives its ability to adapt to environmental changes, develop resistance to drugs, and evade the host immune system. Knowledge on P. falciparum genetic diversity across populations and intervention landscape is thus critical for the implementation of new strategies to eliminate malaria. This study assessed genetic variation with 12,177 high-quality SNPs from 830 Malian P. falciparum isolates collected between 2007 and 2017 from seven locations. The complexity of infections remained high, varied between sites, and showed a trend toward overall decreasing complexity over the decade. Though there was no significant substructure, allele frequencies varied geographically, partly driven by temporal variance in sampling, particularly for drug resistance and antigen loci. Thirty-two mutations in known drug resistance markers (pfcrt, pfdhps, pfdhfr, pfmdr1, pfmdr2, and pfk13) attained a frequency of at least 2% in the populations. SNPs within and around the major markers of resistance to quinolines (pfmdr1 and pfcrt) and antifolates (pfdhfr and pfdhps) varied temporally and geographically, with strong linkage disequilibrium and signatures of directional selection in the genome. These geo-temporal populations also differentiated at alleles in immune-related loci, including, protein E140, pfsurfin8, pfclag8, and pfceltos, as well as pftrap, which showed signatures of haplotype differentiation between populations. Several regions across the genomes, including five known drug resistance loci, showed signatures of differential positive selection. These results suggest that drugs and immune pressure are dominant selective forces against P. falciparum in Mali, but their effect on the parasite genome varies temporally and spatially. Interventions interacting with these genomic variants need to be routinely evaluated as malaria elimination strategies are implemented.