Analysis of Plasmodium falciparum diversity in natural infections by deep sequencing

Diversity and selection analyses identify transmission-blocking antigens as the optimal vaccine candidates in Plasmodium falciparum

Background

A highly effective vaccine for malaria remains an elusive target, at least in part due to the under-appreciated natural parasite variation. This study aimed to investigate genetic and structural variation, and immune selection of leading malaria vaccine candidates across the Plasmodium falciparum 's life cycle.

Methods

We analyzed 325 P. falciparum whole genome sequences from Zambia, in addition to 791 genomes from five other African countries available in the MalariaGEN Pf3k Rdatabase. Ten vaccine antigens spanning three life-history stages were examined for genetic and structural variations, using population genetics measures, haplotype network analysis, and 3D structure selection analysis.

Findings

Among the ten antigens analyzed, only three in the transmission-blocking vaccine category display P . falciparum 3D7 as the dominant haplotype. The antigens AMA1, CSP, MSP1 19 and CelTOS, are much more diverse than the other antigens, and their epitope regions are under moderate to strong balancing selection. In contrast, Rh5 , a blood stage antigen, displays low diversity yet slightly stronger immune selection in the merozoite-blocking epitope region. Except for CelTOS , the transmission-blocking antigens Pfs25 , Pfs48/45 , Pfs230 , Pfs47 , and Pfs28 exhibit minimal diversity and no immune selection in epitopes that induce strain-transcending antibodies, suggesting potential effectiveness of 3D7-based vaccines in blocking transmission.

Interpretations

These findings offer valuable insights into the selection of optimal vaccine candidates against P. falciparum . Based on our results, we recommend prioritizing conserved merozoite antigens and transmission-blocking antigens. Combining these antigens in multi-stage approaches may be particularly promising for malaria vaccine development initiatives.

Funding

Purdue Department of Biological Sciences; Puskas Memorial Fellowship; National Institute of Allergy and Infectious Diseases (U19AI089680).

Research in context

Evidence before this study: Decades of research on the most virulent malaria parasite, Plasmodium falciparum , have yielded multiple antigen candidates of pre-erythrocytic, blood-stage, and transmission-blocking vaccines in varying stages of development from preclinical development to more advanced clinical trials. The malaria vaccine, RTS,S/AS01, which was constructed using the C-terminal and NANP repeat region of the Circumsporozoite Protein ( CSP ) from the African reference strain 3D7, was approved and recommended for use in 2021. However, the vaccine's lower efficacy is likely a result of the genetic polymorphism of the target antigen shown by studies on natural variation in CSP . Similarly, another more recent pre-erythrocytic vaccine, R21/Matrix-M, showed great promise in clinical trials and was recommended in late 2023 by the WHO for use for prevention of malaria in children, but is also multi-dose and CSP -based. To maximize vaccine efficacy, it would be more strategic to first understand diversity and variation of antigens across the three types of vaccine classes, targeting various stages of the P. falciparum life cycle. Previous studies have reported analyses of vaccine candidate antigens but were mostly limited to pre-erythrocytic and blood-stage antigens, with less focus on transmission-blocking antigens. These studies revealed that most of the pre-erythrocytic and blood-stage antigens are of high diversity due to balancing selection, posing challenges for vaccine design to encompass the antigenic variation. A search conducted on PubMed on April 1, 2024, for relevant published research which used the terms "malaria vaccine", " Plasmodium falciparum " [not " vivax "], "selection" and "diversity" yielded 48 studies between 1996 and the present day, with only 14 published studies in the past 3 years. This emphasizes the need for more studies assessing genetic diversity and selection of potential P. falciparum vaccine candidates to aid in more effective vaccine development efforts. A similar search with the terms "transmission-blocking vaccine", "malaria", " Plasmodium falciparum ", not " vivax ", "selection" and "diversity" without any date or language restrictions revealed three relevant studies. This warrants future studies to explore transmission-blocking vaccines in this context. Added value of this study: By comparing the genetic and structural analyses of transmission-blocking antigens with pre-erythrocytic and blood-stage antigens, we identify promising P. falciparum vaccine antigens characterized by their conservation with low balancing selection and the presence of infection/transmission-blocking epitopes, which are essential for informing the development of new malaria vaccines. This comprehensive workflow can be adopted for studying the genetic and structural variation of other P. falciparum vaccine targets before developing the next generation of malaria vaccines for effectiveness against natural parasite populations. Implications of this study: Our suggested strategies for designing malaria vaccines include two possible approaches. We emphasize the development of a multi-stage vaccine that combines critical components such as anti-merozoite ( Rh5 ) and transmission-blocking antigens ( Pfs25 , Pfs28 , Pfs48/45 , Pfs230 ). Alternatively, we suggest the creation of transmission-blocking vaccines specifically targeting Pfs25 , Pfs28 and Pfs48/45 . These innovative approaches show great potential in advancing the development of more potent and effective malaria vaccines for the future.

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.

Genome Analysis of Plasmodium falciparum: A Preliminary Observation - Sierra Leone, 2022-2023

What is already known about this topic?

Sierra Leone, with a gross domestic product (GDP) per capita below $300 and significant poverty, ranks among the world's least developed countries (LDCs). Despite its modest population of 8.6 million, the nation reports approximately 2.6 million malaria cases annually. Previously, there has been no reporting on the malaria genome data from this country.

What is added by this report?

In this study, we present the first reported whole-genome sequence analysis of 19 high parasite-density Plasmodium falciparum isolates from Sierra Leone, providing insights into the genomic epidemiology of this high-prevalence area. We found a high degree of relatedness among infections and substantial genetic diversity, consistent with the gradual reduction in overall case numbers. Moreover, our whole-genome analysis revealed that, beyond drug-resistance genes, gene families related to blood cell invasion, immune evasion, and others are undergoing directional selection. This suggests that the population in Sierra Leone has developed a relatively strong acquired immunity.

What are the implications for public health practice?

The genomic data not only facilitate the creation of single nucleotide polymorphism barcodes for case tracking but also enable the analysis of evolving transmission dynamics and selection pressures. Additionally, the samples from Sierra Leone exhibited higher selective pressures on resistance genes compared to those from Asia, a trend not commonly observed in other African samples. This suggests that less stringent healthcare systems and inconsistent treatment strategies can subject parasites to increased drug pressure, thereby accelerating the development of resistant strains.

Genomic malaria surveillance of antenatal care users detects reduced transmission following elimination interventions in Mozambique

Routine sampling of pregnant women at first antenatal care (ANC) visits could make Plasmodium falciparum genomic surveillance more cost-efficient and convenient in sub-Saharan Africa. We compare the genetic structure of parasite populations sampled from 289 first ANC users and 93 children from the community in Mozambique between 2015 and 2019. Samples are amplicon sequenced targeting 165 microhaplotypes and 15 drug resistance genes. Metrics of genetic diversity and relatedness, as well as the prevalence of drug resistance markers, are consistent between the two populations. In an area targeted for elimination, intra-host genetic diversity declines in both populations (p = 0.002-0.007), while for the ANC population, population genetic diversity is also lower (p = 0.0004), and genetic relatedness between infections is higher (p = 0.002) than control areas, indicating a recent reduction in the parasite population size. These results highlight the added value of genomic surveillance at ANC clinics to inform about changes in transmission beyond epidemiological data.

Genetic complexity alters drug susceptibility of asexual and gametocyte stages of Plasmodium falciparum to antimalarial candidates

Malaria elimination requires interventions able to target both the asexual blood stage (ABS) parasites and transmissible gametocyte stages of Plasmodium falciparum. Lead antimalarial candidates are evaluated against clinical isolates to address key concerns regarding efficacy and to confirm that the current, circulating parasites from endemic regions lack resistance against these candidates. While this has largely been performed on ABS parasites, limited data are available on the transmission-blocking efficacy of compounds with multistage activity. Here, we evaluated the efficacy of lead antimalarial candidates against both ABS parasites and late-stage gametocytes side-by-side, against clinical P. falciparum isolates from southern Africa. We additionally correlated drug efficacy to the genetic diversity of the clinical isolates as determined with a panel of well-characterized, genome-spanning microsatellite markers. Our data indicate varying sensitivities of the isolates to key antimalarial candidates, both for ABS parasites and gametocyte stages. While ABS parasites were efficiently killed, irrespective of genetic complexity, antimalarial candidates lost some gametocytocidal efficacy when the gametocytes originated from genetically complex, multiple-clone infections. This suggests a fitness benefit to multiclone isolates to sustain transmission and reduce drug susceptibility. In conclusion, this is the first study to investigate the efficacy of antimalarial candidates on both ABS parasites and gametocytes from P. falciparum clinical isolates where the influence of parasite genetic complexity is highlighted, ultimately aiding the malaria elimination agenda.

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.

Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana

Plasmodium parasites, the causal agents of malaria, are eukaryotic organisms that obligately undergo sexual recombination within mosquitoes. However, in low transmission settings where most mosquitoes become infected with only a single parasite clone, parasites recombine with themselves, and the clonal lineage is propagated rather than broken up by outcrossing. We investigated whether stochastic/neutral factors drive the persistence and abundance of Plasmodium falciparum clonal lineages in Guyana, a country with relatively low malaria transmission, but the only setting in the Americas in which an important artemisinin resistance mutation (pfk13 C580Y) has been observed. To investigate whether this clonality was potentially associated with the persistence and spatial spread of the mutation, we performed whole genome sequencing on 1,727 Plasmodium falciparum samples collected from infected patients across a five-year period (2016-2021). We characterized the relatedness between each pair of monoclonal infections (n=1,409) through estimation of identity by descent (IBD) and also typed each sample for known or candidate drug resistance mutations. A total of 160 clones (mean IBD ≥ 0.90) were circulating in Guyana during the study period, comprising 13 highly related clusters (mean IBD ≥ 0.40). In the five-year study period, we observed a decrease in frequency of a mutation associated with artemisinin partner drug (piperaquine) resistance (pfcrt C350R) and limited co-occurence of pfcrt C350R with duplications of plasmepsin 2/3, an epistatic interaction associated with piperaquine resistance. We additionally report polymorphisms exhibiting evidence of selection for drug resistance or other phenotypes and reported a novel pfk13 mutation (G718S) as well as 61 nonsynonymous substitutions that increased markedly in frequency. However, P. falciparum clonal dynamics in Guyana appear to be largely driven by stochastic factors, in contrast to other geographic regions. The use of multiple artemisinin combination therapies in Guyana may have contributed to the disappearance of the pfk13 C580Y mutation.

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.

Development of circulating isolates of Plasmodium falciparum is accelerated in Anopheles vectors with reduced reproductive output

Anopheles gambiae and its sibling species Anopheles coluzzii are the most efficient vectors of the malaria parasite Plasmodium falciparum. When females of these species feed on an infected human host, oogenesis and parasite development proceed concurrently, but interactions between these processes are not fully understood. Using multiple natural P. falciparum isolates from Burkina Faso, we show that in both vectors, impairing steroid hormone signaling to disrupt oogenesis leads to accelerated oocyst growth and in a manner that appears to depend on both parasite and mosquito genotype. Consistently, we find that egg numbers are negatively linked to oocyst size, a metric for the rate of oocyst development. Oocyst growth rates are also strongly accelerated in females that are in a pre-gravid state, i.e. that fail to develop eggs after an initial blood meal. Overall, these findings advance our understanding of mosquito-parasite interactions that influence P. falciparum development in malaria-endemic regions.

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.

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

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

Genomics of Plasmodium vivax in Colombia reveals evidence of local bottle-necking and inter-country connectivity in the Americas

Colombia aims to eliminate malaria by 2030 but remains one of the highest burden countries in the Americas. Plasmodium vivax contributes half of all malaria cases, with its control challenged by relapsing parasitaemia, drug resistance and cross-border spread. Using 64 Colombian P. vivax genomes collected between 2013 and 2017, we explored diversity and selection in two major foci of transmission: Chocó and Córdoba. Open-access data from other countries were used for comparative assessment of drug resistance candidates and to assess cross-border spread. Across Colombia, polyclonal infections were infrequent (12%), and infection connectivity was relatively high (median IBD = 5%), consistent with low endemicity. Chocó exhibited a higher frequency of polyclonal infections (23%) than Córdoba (7%), although the difference was not significant (P = 0.300). Most Colombian infections carried double pvdhfr (95%) and single pvdhps (71%) mutants, but other drug resistance mutations were less prevalent (< 10%). There was no evidence of selection at the pvaat1 gene, whose P. falciparum orthologue has recently been implicated in chloroquine resistance. Global population comparisons identified other putative adaptations. Within the Americas, low-level connectivity was observed between Colombia and Peru, highlighting potential for cross-border spread. Our findings demonstrate the potential of molecular data to inform on infection spread and adaptation.

MOIRE: A software package for the estimation of allele frequencies and effective multiplicity of infection from polyallelic data

Malaria parasite genetic data can provide insight into parasite phenotypes, evolution, and transmission. However, estimating key parameters such as allele frequencies, multiplicity of infection (MOI), and within-host relatedness from genetic data has been challenging, particularly in the presence of multiple related coinfecting strains. Existing methods often rely on single nucleotide polymorphism (SNP) data and do not account for within-host relatedness. In this study, we introduce a Bayesian approach called MOIRE (Multiplicity Of Infection and allele frequency REcovery), designed to estimate allele frequencies, MOI, and within-host relatedness from genetic data subject to experimental error. Importantly, MOIRE is flexible in accommodating both polyallelic and SNP data, making it adaptable to diverse genotyping panels. We also introduce a novel metric, the effective MOI ( e M O I ) , which integrates MOI and within-host relatedness, providing a robust and interpretable measure of genetic diversity. Using extensive simulations and real-world data from a malaria study in Namibia, we demonstrate the superior performance of MOIRE over naive estimation methods, accurately estimating MOI up to 7 with moderate sized panels of diverse loci (e.g. microhaplotypes). MOIRE also revealed substantial heterogeneity in population mean MOI and mean relatedness across health districts in Namibia, suggesting detectable differences in transmission dynamics. Notably, e M O I emerges as a portable metric of within-host diversity, facilitating meaningful comparisons across settings, even when allele frequencies or genotyping panels are different. MOIRE represents an important addition to the analysis toolkit for malaria population dynamics. Compared to existing software, MOIRE enhances the accuracy of parameter estimation and enables more comprehensive insights into within-host diversity and population structure. Additionally, MOIRE's adaptability to diverse data sources and potential for future improvements make it a valuable asset for research on malaria and other organisms, such as other eukaryotic pathogens. MOIRE is available as an R package at https://eppicenter.github.io/moire/.

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.

A Plasmodium falciparum genetic cross reveals the contributions of pfcrt and plasmepsin II/III to piperaquine drug resistance

Piperaquine (PPQ) is widely used in combination with dihydroartemisinin (DHA) as a first-line treatment against malaria parasites. Multiple genetic drivers of PPQ resistance have been reported, including mutations in the Plasmodium falciparum chloroquine resistance transporter (pfcrt) and increased copies of plasmepsin II/III (pm2/3). We generated a cross between a Cambodia-derived multi-drug resistant KEL1/PLA1 lineage isolate (KH004) and a drug susceptible parasite isolated in Malawi (Mal31). Mal31 harbors a wild-type (3D7-like) pfcrt allele and a single copy of pm2/3, while KH004 has a chloroquine-resistant (Dd2-like) pfcrt allele with an additional G367C substitution and four copies of pm2/3. We recovered 104 unique recombinant progeny and examined a targeted set of progeny representing all possible combinations of variants at pfcrt and pm2/3 for detailed analysis of competitive fitness and a range of PPQ susceptibility phenotypes, including PPQ survival assay (PSA), area under the dose-response curve (AUC), and a limited point IC50 (LP-IC50). We find that inheritance of the KH004 pfcrt allele is required for PPQ resistance, whereas copy number variation in pm2/3 further enhances resistance but does not confer resistance in the absence of PPQ-R-associated mutations in pfcrt. Deeper investigation of genotype-phenotype relationships demonstrates that progeny clones from experimental crosses can be used to understand the relative contributions of pfcrt, pm2/3, and parasite genetic background, to a range of PPQ-related traits and confirm the critical role of the PfCRT G367C substitution in PPQ resistance.

SNP-Slice Resolves Mixed Infections: Simultaneously Unveiling Strain Haplotypes and Linking Them to Hosts

Multi-strain infection is a common yet under-investigated phenomenon of many pathogens. Currently, biologists analyzing SNP information have to discard mixed infection samples, because existing downstream analyses require monogenomic inputs. Such a protocol impedes our understanding of the underlying genetic diversity, co-infection patterns, and genomic relatedness of pathogens. A reliable tool to learn and resolve the SNP haplotypes from polygenomic data is an urgent need in molecular epidemiology. In this work, we develop a slice sampling Markov Chain Monte Carlo algorithm, named SNP-Slice, to learn not only the SNP haplotypes of all strains in the populations but also which strains infect which hosts. Our method reconstructs SNP haplotypes and individual heterozygosities accurately without reference panels and outperforms the state of art methods at estimating the multiplicity of infections and allele frequencies. Thus, SNP-Slice introduces a novel approach to address polygenomic data and opens a new avenue for resolving complex infection patterns in molecular surveillance. We illustrate the performance of SNP-Slice on empirical malaria and HIV datasets and provide recommendations for the practical use of the method.

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.

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.

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.

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.

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.

Population dynamics and drug resistance mutations in Plasmodium falciparum on the Bijagós Archipelago, Guinea-Bissau

Following integrated malaria control interventions, malaria burden on the Bijagós Archipelago has significantly decreased. Understanding the genomic diversity of circulating Plasmodium falciparum malaria parasites can assist infection control, through identifying drug resistance mutations and characterising the complexity of population structure. This study presents the first whole genome sequence data for P. falciparum isolates from the Bijagós Archipelago. Amplified DNA from P. falciparum isolates sourced from dried blood spot samples of 15 asymptomatic malaria cases were sequenced. Using 1.3 million SNPs characterised across 795 African P. falciparum isolates, population structure analyses revealed that isolates from the archipelago cluster with samples from mainland West Africa and appear closely related to mainland populations; without forming a separate phylogenetic cluster. This study characterises SNPs associated with antimalarial drug resistance on the archipelago. We observed fixation of the PfDHFR mutations N51I and S108N, associated with resistance to sulphadoxine-pyrimethamine, and the continued presence of PfCRT K76T, associated with chloroquine resistance. These data have relevance for infection control and drug resistance surveillance; particularly considering expected increases in antimalarial drug use following updated WHO recommendations, and the recent implementation of seasonal malaria chemoprevention and mass drug administration in the region.

Analysis of sequence diversity in Plasmodium falciparum glutamic acid-rich protein (PfGARP), an asexual blood stage vaccine candidate

Glutamic acid-rich protein of Plasmodium falciparum (PfGARP) binds to erythrocyte band 3 and may enhance cytoadherence of infected erythrocytes. Naturally acquired anti-PfGARP antibodies could confer protection against high parasitemia and severe symptoms. While whole genome sequencing analysis has suggested high conservation in this locus, little is known about repeat polymorphism in this vaccine candidate antigen. Direct sequencing was performed from the PCR-amplified complete PfGARP gene of 80 clinical isolates from four malaria endemic provinces in Thailand and an isolate from a Guinean patient. Publicly available complete coding sequences of this locus were included for comparative analysis. Six complex repeat (RI-RVI) and two homopolymeric glutamic acid repeat (E1 and E2) domains were identified in PfGARP. The erythrocyte band 3-binding ligand in domain RIV and the epitope for mAB7899 antibody eliciting in vitro parasite killing property were perfectly conserved across isolates. Repeat lengths in domains RIII and E1-RVI-E2 seemed to be correlated with parasite density of the patients. Sequence variation in PfGARP exhibited genetic differentiation across most endemic areas of Thailand. Phylogenetic tree inferred from this locus has shown that most Thai isolates formed closely related lineages, suggesting local expansion/contractions of repeat-encoding regions. Positive selection was observed in non-repeat region preceding domain RII which corresponded to a helper T cell epitope predicted to be recognized by a common HLA class II among Thai population. Predicted linear B cell epitopes were identified in both repeat and non-repeat domains. Besides length variation in some repeat domains, sequence conservation in non-repeat regions and almost all predicted immunogenic epitopes have suggested that PfGARP-derived vaccine may largely elicit strain-transcending immunity.

Patterns of Heterochromatin Transitions Linked to Changes in the Expression of Plasmodium falciparum Clonally Variant Genes

The survival of malaria parasites in the changing human blood environment largely depends on their ability to alter gene expression by epigenetic mechanisms. The active state of Plasmodium falciparum clonally variant genes (CVGs) is associated with euchromatin characterized by the histone mark H3K9ac, whereas the silenced state is characterized by H3K9me3-based heterochromatin. Expression switches are linked to euchromatin-heterochromatin transitions, but these transitions have not been characterized for the majority of CVGs. To define the heterochromatin distribution patterns associated with the alternative transcriptional states of CVGs, we compared H3K9me3 occupancy at a genome-wide level among several parasite subclones of the same genetic background that differed in the transcriptional state of many CVGs. We found that de novo heterochromatin formation or the complete disruption of a heterochromatin domain is a relatively rare event, and for the majority of CVGs, expression switches can be explained by the expansion or retraction of heterochromatin domains. We identified different modalities of heterochromatin changes linked to transcriptional differences, but despite this complexity, heterochromatin distribution patterns generally enable the prediction of the transcriptional state of specific CVGs. We also found that in some subclones, several var genes were simultaneously in an active state. Furthermore, the heterochromatin levels in the putative regulatory region of the gdv1 antisense noncoding RNA, a regulator of sexual commitment, varied between parasite lines with different sexual conversion rates. IMPORTANCE The malaria parasite P. falciparum is responsible for more than half a million deaths every year. P. falciparum clonally variant genes (CVGs) mediate fundamental host-parasite interactions and play a key role in parasite adaptation to fluctuations in the conditions of the human host. The expression of CVGs is regulated at the epigenetic level by changes in the distribution of a type of chromatin called heterochromatin. Here, we describe at a genome-wide level the changes in the heterochromatin distribution associated with the different transcriptional states of CVGs. Our results also reveal a likely role for heterochromatin at a particular locus in determining the parasite investment in transmission to mosquitoes. Additionally, this data set will enable the prediction of the transcriptional state of CVGs from epigenomic data, which is important for the study of parasite adaptation to the conditions of the host in natural malaria infections.

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.

Structural and genomic analysis of single nucleotide polymorphisms in human host factor endothelial protein C receptor (EPCR) reveals complex interplay with malaria parasites

Plasmodium parasites responsible for malaria follow a complex life cycle of which half takes place inside the human host. Parasites present diverse antigens at different stages of their life cycle and interact with many surface molecules to attach to and enter host cells. The CIDRα1 domain of Plasmodium falciparum Erythrocyte Membrane Protein 1 (PfEMP1) in infected erythrocytes adheres to one such vascular receptor endothelial protein C receptor (EPCR). EPCR is implicated in the pathogenesis of severe malaria as preferential binding of CIDRα1 to endothelium results in widespread sequestration of infected erythrocytes leading to endothelium inflammation and severe disease. A single EPCR variant S219G is clinically reported to provide protection from severe malaria. In this work, we have collated all single nucleotide polymorphisms (SNPs) in EPCR from dbSNP. We structurally mapped the SNPs on the three-dimensional complex of EPCR and PfEMP1 CIDRα1. Analysis shows that most EPCR mutations lie on the receptor surface and are non-conservative. Of the 11 mutations in the CIDRα1-interaction region of EPCR, S88P, L96V/I, and R98L/H/P/C are seen with comparably higher occurrences in diverse populations. Our structural analysis details a framework of the interactions between the parasite ligand and host factor EPCR. These structural glimpses provide a blueprint for designing both field-based variant sequencing studies and vaccine development.

Retrospective analysis of Plasmodium vivax genomes from a pre-elimination China inland population in the 2010s

Introduction

In malaria-free countries, imported cases are challenging because interconnections with neighboring countries with higher transmission rates increase the risk of parasite reintroduction. Establishing a genetic database for rapidly identifying malaria importation or reintroduction is crucial in addressing these challenges. This study aimed to examine genomic epidemiology during the pre-elimination stage by retrospectively reporting whole-genome sequence variation of 10 Plasmodium vivax isolates from inland China.

Methods

The samples were collected during the last few inland outbreaks from 2011 to 2012 when China implemented a malaria control plan. After next-generation sequencing, we completed a genetic analysis of the population, explored the geographic specificity of the samples, and examined clustering of selection pressures. We also scanned genes for signals of positive selection.

Results

China's inland populations were highly structured compared to the surrounding area, with a single potential ancestor. Additionally, we identified genes under selection and evaluated the selection pressure on drug-resistance genes. In the inland population, positive selection was detected in some critical gene families, including sera, msp3, and vir. Meanwhile, we identified selection signatures in drug resistance, such as ugt, krs1, and crt, and noticed that the ratio of wild-type dhps and dhfr-ts increased after China banned sulfadoxine-pyrimethamine (SP) for decades.

Discussion

Our data provides an opportunity to investigate the molecular epidemiology of pre-elimination inland malaria populations, which exhibited lower selection pressure on invasion and immune evasion genes than neighbouring areas, but increased drug resistance in low transmission settings. Our results revealed that the inland population was severely fragmented with low relatedness among infections, despite a higher incidence of multiclonal infections, suggesting that superinfection or co-transmission events are rare in low-endemic circumstances. We identified selective signatures of resistance and found that the proportion of susceptible isolates fluctuated in response to the prohibition of specific drugs. This finding is consistent with the alterations in medication strategies during the malaria elimination campaign in inland China. Such findings could provide a genetic basis for future population studies, assessing changes in other pre-elimination countries.

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

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

Short tandem repeat polymorphism in the promoter region of cyclophilin 19B drives its transcriptional upregulation and contributes to drug resistance in the malaria parasite Plasmodium falciparum

Resistance of the human malaria parasites, Plasmodium falciparum, to artemisinins is now fully established in Southeast Asia and is gradually emerging in Sub-Saharan Africa. Although nonsynonymous SNPs in the pfk13 Kelch-repeat propeller (KREP) domain are clearly associated with artemisinin resistance, their functional relevance requires cooperation with other genetic factors/alterations of the P. falciparum genome, collectively referred to as genetic background. Here we provide experimental evidence that P. falciparum cyclophilin 19B (PfCYP19B) may represent one putative factor in this genetic background, contributing to artemisinin resistance via its increased expression. We show that overexpression of PfCYP19B in vitro drives limited but significant resistance to not only artemisinin but also piperaquine, an important partner drug in artemisinin-based combination therapies. We showed that PfCYP19B acts as a negative regulator of the integrated stress response (ISR) pathway by modulating levels of phosphorylated eIF2α (eIF2α-P). Curiously, artemisinin and piperaquine affect eIF2α-P in an inverse direction that in both cases can be modulated by PfCYP19B towards resistance. Here we also provide evidence that the upregulation of PfCYP19B in the drug-resistant parasites appears to be maintained by a short tandem repeat (SRT) sequence polymorphism in the gene's promoter region. These results support a model that artemisinin (and other drugs) resistance mechanisms are complex genetic traits being contributed to by altered expression of multiple genes driven by genetic polymorphism at their promoter regions.

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

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

Measurably recombining malaria parasites

Genomic epidemiology has guided research and policy for various viral pathogens and there has been a parallel effort towards using genomic epidemiology to combat diseases that are caused by eukaryotic pathogens, such as the malaria parasite. However, the central concept of viral genomic epidemiology, namely that of measurably mutating pathogens, does not apply easily to sexually recombining parasites. Here we introduce the related but different concept of measurably recombining malaria parasites to promote convergence around a unifying theoretical framework for malaria genomic epidemiology. Akin to viral phylodynamics, we anticipate that an inferential framework developed around recombination will help guide practical research and thus realize the full public health potential of genomic epidemiology for malaria parasites and other sexually recombining pathogens.

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.

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

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

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.

Structure-guided insights into potential function of novel genetic variants in the malaria vaccine candidate PfRh5

The recent stall in the global reduction of malaria deaths has made the development of a highly effective vaccine essential. A major challenge to developing an efficacious vaccine is the extensive diversity of Plasmodium falciparum antigens. While genetic diversity plays a major role in immune evasion and is a barrier to the development of both natural and vaccine-induced protective immunity, it has been under-prioritized in the evaluation of malaria vaccine candidates. This study uses genomic approaches to evaluate genetic diversity in next generation malaria vaccine candidate PfRh5. We used targeted deep amplicon sequencing to identify non-synonymous Single Nucleotide Polymorphisms (SNPs) in PfRh5 (Reticulocyte-Binding Protein Homologue 5) in 189 P. falciparum positive samples from Southern Senegal and identified 74 novel SNPs. We evaluated the population prevalence of these SNPs as well as the frequency in individual samples and found that only a single SNP, C203Y, was present at every site. Many SNPs were unique to the individual sampled, with over 90% of SNPs being found in just one infected individual. In addition to population prevalence, we assessed individual level SNP frequencies which revealed that some SNPs were dominant (frequency of greater than 25% in a polygenomic sample) whereas most were rare, present at 2% or less of total reads mapped to the reference at the given position. Structural modeling uncovered 3 novel SNPs occurring under epitopes bound by inhibitory monoclonal antibodies, potentially impacting immune evasion, while other SNPs were predicted to impact PfRh5 structure or interactions with the receptor or binding partners. Our data demonstrate that PfRh5 exhibits greater genetic diversity than previously described, with the caveat that most of the uncovered SNPs are at a low overall frequency in the individual and prevalence in the population. The structural studies reveal that novel SNPs could have functional implications on PfRh5 receptor binding, complex formation, or immune evasion, supporting continued efforts to validate PfRh5 as an effective malaria vaccine target and development of a PfRh5 vaccine.

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.

Population genomics in neglected malaria parasites

Malaria elimination includes neglected human malaria parasites Plasmodium vivax, Plasmodium ovale spp., and Plasmodium malariae. Biological features such as association with low-density infection and the formation of hypnozoites responsible for relapse make their elimination challenging. Studies on these parasites rely primarily on clinical samples due to the lack of long-term culture techniques. With improved methods to enrich parasite DNA from clinical samples, whole-genome sequencing of the neglected malaria parasites has gained increasing popularity. Population genomics of more than 2200 P. vivax global isolates has improved our knowledge of parasite biology and host-parasite interactions, identified vaccine targets and potential drug resistance markers, and provided a new way to track parasite migration and introduction and monitor the evolutionary response of local populations to elimination efforts. Here, we review advances in population genomics for neglected malaria parasites, discuss how the rich genomic information is being used to understand parasite biology and epidemiology, and explore opportunities for the applications of malaria genomic data in malaria elimination practice.

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.

Genomic analysis of single nucleotide polymorphisms in malaria parasite drug targets

Malaria is a life-threatening parasitic disease caused by members of the genus Plasmodium. The development and spread of drug-resistant strains of Plasmodium parasites represent a major challenge to malaria control and elimination programmes. Evaluating genetic polymorphism in a drug target improves our understanding of drug resistance and facilitates drug design. Approximately 450 and 19 whole-genome assemblies of Plasmodium falciparum and Plasmodium vivax, respectively, are currently available, and numerous sequence variations have been found due to the presence of single nucleotide polymorphism (SNP). In the study reported here, we analysed global SNPs in the malaria parasite aminoacyl-tRNA synthetases (aaRSs). Our analysis revealed 3182 unique SNPs in the 20 cytoplasmic P. falciparum aaRSs. Structural mapping of SNPs onto the three-dimensional inhibitor-bound complexes of the three advanced drug targets within aaRSs revealed a remarkably low mutation frequency in the crucial aminoacylation domains, low overall occurrence of mutations across samples and high conservation in drug/substrate binding regions. In contrast to aaRSs, dihydropteroate synthase (DHPS), also a malaria drug target, showed high occurrences of drug resistance-causing mutations. Our results show that it is pivotal to screen potent malaria drug targets against global SNP profiles to assess genetic variances to ensure success in designing drugs against validated targets and tackle drug resistance early on.

A GPU-accelerated compute framework for pathogen genomic variant identification to aid genomic epidemiology of infectious disease: a malaria case study

As recently demonstrated by the COVID-19 pandemic, large-scale pathogen genomic data are crucial to characterize transmission patterns of human infectious diseases. Yet, current methods to process raw sequence data into analysis-ready variants remain slow to scale, hampering rapid surveillance efforts and epidemiological investigations for disease control. Here, we introduce an accelerated, scalable, reproducible, and cost-effective framework for pathogen genomic variant identification and present an evaluation of its performance and accuracy across benchmark datasets of Plasmodium falciparum malaria genomes. We demonstrate superior performance of the GPU framework relative to standard pipelines with mean execution time and computational costs reduced by 27× and 4.6×, respectively, while delivering 99.9% accuracy at enhanced reproducibility.

Potential Opportunities and Challenges of Deploying Next Generation Sequencing and CRISPR-Cas Systems to Support Diagnostics and Surveillance Towards Malaria Control and Elimination in Africa

Recent developments in molecular biology and genomics have revolutionized biology and medicine mainly in the developed world. The application of next generation sequencing (NGS) and CRISPR-Cas tools is now poised to support endemic countries in the detection, monitoring and control of endemic diseases and future epidemics, as well as with emerging and re-emerging pathogens. Most low and middle income countries (LMICs) with the highest burden of infectious diseases still largely lack the capacity to generate and perform bioinformatic analysis of genomic data. These countries have also not deployed tools based on CRISPR-Cas technologies. For LMICs including Tanzania, it is critical to focus not only on the process of generation and analysis of data generated using such tools, but also on the utilization of the findings for policy and decision making. Here we discuss the promise and challenges of NGS and CRISPR-Cas in the context of malaria as Africa moves towards malaria elimination. These innovative tools are urgently needed to strengthen the current diagnostic and surveillance systems. We discuss ongoing efforts to deploy these tools for malaria detection and molecular surveillance highlighting potential opportunities presented by these innovative technologies as well as challenges in adopting them. Their deployment will also offer an opportunity to broadly build in-country capacity in pathogen genomics and bioinformatics, and to effectively engage with multiple stakeholders as well as policy makers, overcoming current workforce and infrastructure challenges. Overall, these ongoing initiatives will build the malaria molecular surveillance capacity of African researchers and their institutions, and allow them to generate genomics data and perform bioinformatics analysis in-country in order to provide critical information that will be used for real-time policy and decision-making to support malaria elimination on the continent.

Pfs47 as a Malaria Transmission-Blocking Vaccine Target

Transmission-blocking vaccines (TBVs), pioneered by Richard Carter and others, aim to prevent parasite development in the mosquito vector and are a promising new tool for malaria elimination. Pfs47, recently identified as a TBV target, is a three-domain 6-cysteine protein on the surface of Plasmodium falciparum sexual stages. Pfs47 allows the parasite to evade mosquito immunity and is key for P. falciparum infection of the dominant malaria vectors Anopheles gambiae, Anopheles dirus, and Anopheles albimanus. Antibodies against Pfs47 domain 2 (D2) have significant transmission-blocking activity that prevents Plasmodium ookinete development and is independent of human complement. Strong transmission-blocking activity has been mapped to a region of 52 amino acids in Pfs47 D2. Efforts to optimize the immunogenicity of the Pfs47 D2 antigen with a viral-like particle have been successful, and the efficacy of a P47-based TBV was confirmed in vivo with Pbs47, the orthologue of Pfs47 in the mouse malaria parasite Plasmodium berghei. The current evidence warrants further development and clinical testing of a Pfs47-based TBV.

Evaluation of the Pfs25-IMX313/Matrix-M malaria transmission-blocking candidate vaccine in endemic settings

Malaria control relies heavily on the use of anti-malarial drugs and insecticides against malaria parasites and mosquito vectors. Drug and insecticide resistance threatens the effectiveness of conventional malarial interventions; alternative control approaches are, therefore, needed. The development of malaria transmission-blocking vaccines that target the sexual stages in humans or mosquito vectors is among new approaches being pursued. Here, the immunological mechanisms underlying malaria transmission blocking, status of Pfs25-based vaccines are viewed, as well as approaches and capacity for first in-human evaluation of a transmission-blocking candidate vaccine Pfs25-IMX313/Matrix-M administered to semi-immune healthy individuals in endemic settings. It is concluded that institutions in low and middle income settings should be supported to conduct first-in human vaccine trials in order to stimulate innovative research and reduce the overdependence on developed countries for research and local interventions against many diseases of public health importance.

Recent genetic exchanges and admixture shape the genome and population structure of the zoonotic pathogen Cryptosporidium parvum

Cryptosporidium parvum is a globally distributed zoonotic pathogen and a major cause of diarrhoeal disease in humans and ruminants. The parasite's life cycle comprises an obligatory sexual phase, during which genetic exchanges can occur between previously isolated lineages. Here, we compare 32 whole genome sequences from human- and ruminant-derived parasite isolates collected across Europe, Egypt and China. We identify three strongly supported clusters that comprise a mix of isolates from different host species, geographic origins, and subtypes. We show that: (1) recombination occurs between ruminant isolates into human isolates; (2) these recombinant regions can be passed on to other human subtypes through gene flow and population admixture; (3) there have been multiple genetic exchanges, and most are likely recent; (4) putative virulence genes are significantly enriched within these genetic exchanges, and (5) this results in an increase in their nucleotide diversity. We carefully dissect the phylogenetic sequence of two genetic exchanges, illustrating the long-term evolutionary consequences of these events. Our results suggest that increased globalisation and close human-animal contacts increase the opportunity for genetic exchanges between previously isolated parasite lineages, resulting in spillover and spillback events. We discuss how this can provide a novel substrate for natural selection at genes involved in host-parasite interactions, thereby potentially altering the dynamic coevolutionary equilibrium in the Red Queens arms race.

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

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

Amplicon sequencing as a potential surveillance tool for complexity of infection and drug resistance markers in Plasmodium falciparum asymptomatic infections

Background

Genotyping Plasmodium falciparum subpopulations in malaria infections is an important aspect of malaria molecular epidemiology to understand within-host diversity and the frequency of drug resistance markers.

Methods

We characterized P. falciparum genetic diversity in asymptomatic infections and subsequent first febrile infections using amplicon sequencing (AmpSeq) of ama1 in Coastal Kenya. We also examined temporal changes in haplotype frequencies of mdr1, a drug-resistant marker.

Results

We found >60% of the infections were polyclonal (complexity of infection [COI] >1) and there was a reduction in COI over time. Asymptomatic infections had a significantly higher mean COI than febrile infections based on ama1 sequences (2.7 [95% confidence interval {CI}, 2.65–2.77] vs 2.22 [95% CI, 2.17–2.29], respectively). Moreover, an analysis of 30 paired asymptomatic and first febrile infections revealed that many first febrile infections (91%) were due to the presence of new ama1 haplotypes. The mdr1-YY haplotype, associated with chloroquine and amodiaquine resistance, decreased over time, while the NY (wild type) and the NF (modulates response to lumefantrine) haplotypes increased.

Conclusions

This study emphasizes the utility of AmpSeq in characterizing parasite diversity as it can determine relative proportions of clones and detect minority clones. The usefulness of AmpSeq in antimalarial drug resistance surveillance is also highlighted.