A genome-wide map of diversity in Plasmodium falciparum

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.

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.

Ancient diversity in host-parasite interaction genes in a model parasitic nematode

Host-parasite interactions exert strong selection pressures on the genomes of both host and parasite. These interactions can lead to negative frequency-dependent selection, a form of balancing selection that is hypothesised to explain the high levels of polymorphism seen in many host immune and parasite antigen loci. Here, we sequence the genomes of several individuals of Heligmosomoides bakeri, a model parasite of house mice, and Heligmosomoides polygyrus, a closely related parasite of wood mice. Although H. bakeri is commonly referred to as H. polygyrus in the literature, their genomes show levels of divergence that are consistent with at least a million years of independent evolution. The genomes of both species contain hyper-divergent haplotypes that are enriched for proteins that interact with the host immune response. Many of these haplotypes originated prior to the divergence between H. bakeri and H. polygyrus, suggesting that they have been maintained by long-term balancing selection. Together, our results suggest that the selection pressures exerted by the host immune response have played a key role in shaping patterns of genetic diversity in the genomes of parasitic nematodes.

Malaria Genomics, Vaccine Development, and Microbiome

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

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.

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.

Modelling the relative cost-effectiveness of the RTS,S/AS01 malaria vaccine compared to investment in vector control or chemoprophylaxis

BACKGROUND

The World Health Organization has recommended a 4-dose schedule of the RTS,S/AS01 (RTS,S) vaccine for children in regions of moderate to high P. falciparum transmission. Faced with limited supply and finite resources, global funders and domestic malaria control programs will need to examine the relative cost-effectiveness of RTS,S and identify target areas for vaccine implementation relative to scale-up of existing interventions.

METHODS

Using an individual-based mathematical model of P. falciparum, we modelled the cost-effectiveness of RTS,S across a range of settings in sub-Saharan Africa, incorporating various rainfall patterns, insecticide-treated net (ITN) use, treatment coverage, and parasite prevalence bands. We compare age-based and seasonal RTS,S administration to increasing ITN usage, switching to next generation ITNs in settings experiencing insecticide-resistance, and introduction of seasonal malaria chemoprevention (SMC) in areas of seasonal transmission.

RESULTS

For RTS,S to be the most cost-effective intervention option considered, the maximum cost per dose was less than $9.30 USD in 90.9% of scenarios. Nearly all (89.8%) values at or above $9.30 USD per dose were in settings with 60% established bed net use and / or with established SMC, and 76.3% were in the highest PfPR_2-10 band modelled (40%). Addition of RTS,S to strategies involving 60% ITN use, increased ITN usage or a switch to PBO nets, and SMC, if eligible, still led to significant marginal case reductions, with a median of 2,653 (IQR: 1,741 to 3,966) cases averted per 100,000 people annually, and 82,270 (IQR: 54,034 to 123,105) cases averted per 100,000 fully vaccinated children (receiving at least three doses).

CONCLUSIONS

Use of RTS,S results in reductions in malaria cases and deaths even when layered upon existing interventions. When comparing relative cost-effectiveness, scale up of ITNs, introduction of SMC, and switching to new technology nets should be prioritized in eligible settings.

Genetic diversity and natural selection of rif gene (PF3D7_1254800) in the Plasmodium falciparum global populations

To reveal the genetic characteristics of one member of the Plasmodium falciparum repetitive interspersed family (rif), we sequenced the rif gene (PF3D7_1254800) in 53 field isolates collected from Ghana-imported cases into China and compared them with 350 publicly available P. falciparum rif sequences from global populations. In the Ghana-imported population, the nucleotide diversities were 0.05714 and 0.06616 for the full length and variable region of rif gene, respectively. Meanwhile, 22 and 20 haplotypes were identified for the full length and variable region of rif gene (Hd = 0.843 and 0.838, respectively). Diversity of rif gene in Ghana-imported population was higher than that observed in Cambodia, Thailand, Vietnam, Myanmar, Mali, Ghana, and Senegal populations. In this analysis, we found high genetic diversity of rif gene in global P. falciparum populations and identified 158 haplotypes. Tajima's D-test shows that there are large differences in the direction of selection between the conserved and variable region of rif gene. Tajima's D value for the variable region was 0.20074, indicating that balancing selection existed in this region. We found that the variable region was the main target of selection for positive diversification, and most mutation sites were located in this region. The population structure suggested optimized cluster values of K = 6. The five groups in Ghana-imported population included a unique subpopulation. Our results reveal the dynamics of the rif gene (PF3D7_1254800) in P. falciparum populations, which can aid in the rational design of P. falciparum rif-based vaccines.

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.

QSAR Classification Models for Prediction of Hydroxamate Histone Deacetylase Inhibitor Activity against Malaria Parasites

Malaria, caused by Plasmodium parasites, results in >400,000 deaths annually. There is no effective vaccine, and new drugs with novel modes of action are needed because of increasing parasite resistance to current antimalarials. Histone deacetylases (HDACs) are epigenetic regulatory enzymes that catalyze post-translational protein deacetylation and are promising malaria drug targets. Here, we describe quantitative structure-activity relationship models to predict the antiplasmodial activity of hydroxamate-based HDAC inhibitors. The models incorporate P. falciparum in vitro activity data for 385 compounds containing a hydroxamic acid and were subject to internal and external validation. When used to screen 22 new hydroxamate-based HDAC inhibitors for antiplasmodial activity, model A7 (external accuracy 91%) identified three hits that were subsequently verified as having potent in vitro activity against P. falciparum parasites (IC₅₀ = 6, 71, and 84 nM), with 8 to 51-fold selectivity for P. falciparum versus human cells.

Advances and opportunities in malaria population genomics

Almost 20 years have passed since the first reference genome assemblies were published for Plasmodium falciparum, the deadliest malaria parasite, and Anopheles gambiae, the most important mosquito vector of malaria in sub-Saharan Africa. Reference genomes now exist for all human malaria parasites and nearly half of the ~40 important vectors around the world. As a foundation for genetic diversity studies, these reference genomes have helped advance our understanding of basic disease biology and drug and insecticide resistance, and have informed vaccine development efforts. Population genomic data are increasingly being used to guide our understanding of malaria epidemiology, for example by assessing connectivity between populations and the efficacy of parasite and vector interventions. The potential value of these applications to malaria control strategies, together with the increasing diversity of genomic data types and contexts in which data are being generated, raise both opportunities and challenges in the field. This Review discusses advances in malaria genomics and explores how population genomic data could be harnessed to further support global disease control efforts.

Variation in selective constraints along the Plasmodium life cycle

Plasmodium parasites, the cause of malaria, have a complex life cycle, infecting alternatively vertebrate hosts and female Anopheles mosquitoes and undergoing intra- and extra-cellular development in several organs of these hosts. Most of the ~5000 protein-coding genes present in Plasmodium genomes are only expressed at specific life stages, and different genes might therefore be subject to different selective pressures depending on the biological activity of the parasite and its microenvironment at this point in development. Here, we estimate the selective constraints on the protein-coding sequences of all annotated genes of rodent and primate Plasmodium parasites and, using data from scRNA-seq experiments spanning many developmental stages, analyze their variation with regard to when these genes are expressed in the parasite life cycle. Our study reveals extensive variation in selective constraints throughout the parasites' development and highlights stages that are evolving more rapidly than others. These findings provide novel insights into the biology of these parasites and could provide important information to develop better treatment strategies or vaccines against these medically-important organisms.

Genomic Approaches to Drug Resistance in Malaria

Although the last two decades have seen a substantial decline in malaria incidence and mortality due to the use of insecticide-treated bed nets and artemisinin combination therapy, the threat of drug resistance is a constant obstacle to sustainable malaria control. Given that patients can die quickly from this disease, public health officials and doctors need to understand whether drug resistance exists in the parasite population, as well as how prevalent it is so they can make informed decisions about treatment. As testing for drug efficacy before providing treatment to malaria patients is impractical, researchers need molecular markers of resistance that can be more readily tracked in parasite populations. To this end, much work has been done to unravel the genetic underpinnings of drug resistance in Plasmodium falciparum. The aim of this review is to provide a broad overview of common genomic approaches that have been used to discover the alleles that drive drug response phenotypes in the most lethal human malaria parasite.

Single-cell sequencing of the small and AT-skewed genome of malaria parasites

Single-cell genomics is a rapidly advancing field; however, most techniques are designed for mammalian cells. We present a single-cell sequencing pipeline for an intracellular parasite, Plasmodium falciparum, with a small genome of extreme base content. Through optimization of a quasi-linear amplification method, we target the parasite genome over contaminants and generate coverage levels allowing detection of minor genetic variants. This work, as well as efforts that build on these findings, will enable detection of parasite heterogeneity contributing to P. falciparum adaptation. Furthermore, this study provides a framework for optimizing single-cell amplification and variant analysis in challenging genomes.

Genomic analysis reveals independent evolution of Plasmodium falciparum populations in Ethiopia

Background

Plasmodium falciparum parasite populations in Ethiopia have been experiencing local selective pressures from drugs and immunity, leading to evolutionary adaptation. However, there was a paucity of data on genomic characterization and evolutionary adaptations of P. falciparum isolates from the central area of Ethiopia.

Methods

Whole-genome analysis of 25 P. falciparum isolates from central Ethiopia, specifically from West Arsi, were studied to determine their genetic diversity, population structures, and signatures of selection in known drug resistance alleles against global isolates from Cambodia, Thailand, DR Congo, and Malawi.

Results

A total of 18,517 high-quality single-nucleotide polymorphisms (SNPs) were identified in Ethiopian P. falciparum isolates. About 84% of the Ethiopian P. falciparum isolates had a F_WS value > 0.95 showing a dominant single genotype infection in most isolates at the time of collection with little potential for out-crossing as expected in areas with low transmission intensity. Within-host diversity of Ethiopian infections was significantly different from East African (p < 0.001), but not Southeast Asian infections (P > 0.05). A significant population structure has been observed by PCA and population differentiation between Ethiopian parasites and East African (Fst ~ 10%) and Southeast Asian populations (Fst ~ 18%), suggesting limited gene flow and the independent evolution of the Ethiopian parasite population. Moreover, a total of 125 genes under balancing selection was found that include ama1, trap, eba175, and lsa3, previously identified as targets of human host immunity. Recent directional selection analysis using integrated standardized haplotype score (IHS) did not detect any selection signatures in the Pfcrt, Pfdhfr, Pfdhps, Pfmdr1, and PfK13 genes. However, known drug resistance-conferring mutations analysis showed that at least one SNP marker was fixed in these genes, but not in Pfdhps and PfK13.

Conclusion

Plasmodium falciparum populations in the central region of Ethiopia was structurally diverged from both Southeast Asian and other East African populations. Malaria infections in Ethiopia had low within-host diversity, and parasites carry fixed chloroquine resistance markers despite the withdrawal of this drug for the treatment of P. falciparum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03660-y.

Genome-Wide Scans for Ghanaian Plasmodium falciparum Genes Under Selection From Local and Chinese Host Populations

Initial malarial infection mostly causes symptomatic illness in humans. Infection that is not fatal induces complete protection from severe illness and death, and thus complete protection from severe illness or death is granted with sufficient exposure. However, malaria parasite immunity necessitates constant exposure. Therefore, it is important to evaluate lowered immunity and recurrent susceptibility to symptomatic disease in lower transmission areas. We aimed to investigate selection pressure based on transmission levels, antimalarial drug use, and environmental factors. We whole genome sequenced (WGS) P. falciparum clinical samples from Chinese hosts working in Ghana and compared the results with the WGS data of isolates from native Ghanaians downloaded from pf3k. The P. falciparum samples were generally clustered according to their geographic origin, and Chinese imported samples showed a clear African origin with a slightly different distribution from the native Ghanaian samples. Moreover, samples collected from two host populations showed evidence of differences in the intensity of selection. Compared with native Ghanaian samples, the China-imported isolates exhibited a higher proportion of monoclonal infections, and many genes associated with RBC invasion and immune evasion were found to be under less selection pressure. There was no significant difference in the selection of drug-resistance genes due to a similar artemisinin-based combination therapy medication profile. Local selection of malarial parasites is considered to be a result of differences in the host immunity or disparity in the transmission opportunities of the host. In China, most P. falciparum infections were imported from Africa, and under these circumstances, distinct local selective pressures may be caused by varying acquired immunity and transmission intensity. This study revealed the impact of host switching on the immune system, and it may provide a better understanding of the mechanisms that enable clinical immunity to malaria.

Evolutionary analyses of the major variant surface antigen-encoding genes reveal population structure of Plasmodium falciparum within and between continents

Malaria remains a major public health problem in many countries. Unlike influenza and HIV, where diversity in immunodominant surface antigens is understood geographically to inform disease surveillance, relatively little is known about the global population structure of PfEMP1, the major variant surface antigen of the malaria parasite Plasmodium falciparum. The complexity of the var multigene family that encodes PfEMP1 and that diversifies by recombination, has so far precluded its use in malaria surveillance. Recent studies have demonstrated that cost-effective deep sequencing of the region of var genes encoding the PfEMP1 DBLα domain and subsequent classification of within host sequences at 96% identity to define unique DBLα types, can reveal structure and strain dynamics within countries. However, to date there has not been a comprehensive comparison of these DBLα types between countries. By leveraging a bioinformatic approach (jumping hidden Markov model) designed specifically for the analysis of recombination within var genes and applying it to a dataset of DBLα types from 10 countries, we are able to describe population structure of DBLα types at the global scale. The sensitivity of the approach allows for the comparison of the global dataset to ape samples of Plasmodium Laverania species. Our analyses show that the evolution of the parasite population emerging out of Africa underlies current patterns of DBLα type diversity. Most importantly, we can distinguish geographic population structure within Africa between Gabon and Ghana in West Africa and Uganda in East Africa. Our evolutionary findings have translational implications in the context of globalization. Firstly, DBLα type diversity can provide a simple diagnostic framework for geographic surveillance of the rapidly evolving transmission dynamics of P. falciparum. It can also inform efforts to understand the presence or absence of global, regional and local population immunity to major surface antigen variants. Additionally, we identify a number of highly conserved DBLα types that are present globally that may be of biological significance and warrant further characterization.

Next Generation Sequencing and Bioinformatics Methodologies for Infectious Disease Research and Public Health: Approaches, Applications, and Considerations for Development of Laboratory Capacity

Next generation sequencing (NGS) combined with bioinformatics has successfully been used in a vast array of analyses for infectious disease research of public health relevance. For instance, NGS and bioinformatics approaches have been used to identify outbreak origins, track transmissions, investigate epidemic dynamics, determine etiological agents of a disease, and discover novel human pathogens. However, implementation of high-quality NGS and bioinformatics in research and public health laboratories can be challenging. These challenges mainly include the choice of the sequencing platform and the sequencing approach, the choice of bioinformatics methodologies, access to the appropriate computation and information technology infrastructure, and recruiting and retaining personnel with the specialized skills and experience in this field. In this review, we summarize the most common NGS and bioinformatics workflows in the context of infectious disease genomic surveillance and pathogen discovery, and highlight the main challenges and considerations for setting up an NGS and bioinformatics-focused infectious disease research public health laboratory. We describe the most commonly used sequencing platforms and review their strengths and weaknesses. We review sequencing approaches that have been used for various pathogens and study questions, as well as the most common difficulties associated with these approaches that should be considered when implementing in a public health or research setting. In addition, we provide a review of some common bioinformatics tools and procedures used for pathogen discovery and genome assembly, along with the most common challenges and solutions. Finally, we summarize the bioinformatics of advanced viral, bacterial, and parasite pathogen characterization, including types of study questions that can be answered when utilizing NGS and bioinformatics.

Restricted genetic heterogeneity of the Plasmodium vivax transmission-blocking vaccine (TBV) candidate Pvs48/45 in a low transmission setting: Implications for the Plasmodium vivax malaria vaccine development

Plasmodium vivax is the most widespread malaria species parasitizing humans outside Africa, with approximately 100 million cases reported per year. Most human cases of P. vivax are asymptomatic with low parasitemia, making active case detection-based elimination programme challenging and less effective. Despite the widespread distribution of P. vivax, no effective vaccines are currently available. Transmission blocking vaccines have recently emerged as potential vaccine candidates to reduce transmission rates to below the essential levels required for the maintenance of the parasite life cycle. Here, we demonstrated that P. vivax was the predominant species found in a malaria-endemic area, although P. vivax/P. falciparum co-infections were also common. Through genomic sequence analysis and neighbor-joining algorithms, we demonstrated limited genetic heterogeneity in the P. vivax transmission-blocking vaccine candidate Pvs48/45 among clinical isolates of P. vivax. Restricted genetic polymorphism occurred at both nucleotide and amino acid levels. The most frequent mutation was A → G at nucleotide position 77 (46.7%), whereas the least frequent was C → T at nucleotide position 1230 (3.3%). The occurrence of single nucleotide polymorphisms (SNPs) distribution at 6/8 positions (75%) led to changes in amino acid sequences in the Pvs48/45 loci, whereas 2/8 (25%) of SNPs resulted in no amino acid sequence variations. Consistently, the nucleotide diversity in the Pvs48/45 locus among the P. vivax population studied was extremely low (π = 0.000525). Changes in amino acid sequences in the Pvs48/45 protein did not result in substantial conformational modifications in the tertiary structures of these proteins. Unveiling the population genetic structure and genetic heterogeneity of vaccine target antigens are necessary for rational design of transmission-blocking antibody vaccines and to monitor the vaccine efficacy in clinical trials in endemic areas for malaria.

Monitoring Plasmodium falciparum and Plasmodium vivax using microsatellite markers indicates limited changes in population structure after substantial transmission decline in Papua New Guinea

Monitoring the genetic structure of pathogen populations may be an economical and sensitive approach to quantify the impact of control on transmission dynamics, highlighting the need for a better understanding of changes in population genetic parameters as transmission declines. Here we describe the first population genetic analysis of two major human malaria parasites, Plasmodium falciparum (Pf) and Plasmodium vivax (Pv), following nationwide distribution of long-lasting insecticide-treated nets (LLINs) in Papua New Guinea (PNG). Parasite isolates from pre- (2005-2006) and post-LLIN (2010-2014) were genotyped using microsatellite markers. Despite parasite prevalence declining substantially (East Sepik Province: Pf = 54.9%-8.5%, Pv = 35.7%-5.6%, Madang Province: Pf = 38.0%-9.0%, Pv: 31.8%-19.7%), genetically diverse and intermixing parasite populations remained. Pf diversity declined modestly post-LLIN relative to pre-LLIN (East Sepik: R_s  = 7.1-6.4, H_E  = 0.77-0.71; Madang: R_s  = 8.2-6.1, H_E  = 0.79-0.71). Unexpectedly, population structure present in pre-LLIN populations was lost post-LLIN, suggesting that more frequent human movement between provinces may have contributed to higher gene flow. Pv prevalence initially declined but increased again in one province, yet diversity remained high throughout the study period (East Sepik: R_s  = 11.4-9.3, H_E  = 0.83-0.80; Madang: R_s  = 12.2-14.5, H_E  = 0.85-0.88). Although genetic differentiation values increased between provinces over time, no significant population structure was observed at any time point. For both species, a decline in multiple infections and increasing clonal transmission and significant multilocus linkage disequilibrium post-LLIN were positive indicators of impact on the parasite population using microsatellite markers. These parameters may be useful adjuncts to traditional epidemiological tools in the early stages of transmission reduction.

Naturally Acquired Humoral Immunity Against Plasmodium falciparum Malaria

Malaria remains a significant contributor to the global burden of disease, with around 40% of the world's population at risk of Plasmodium infections. The development of an effective vaccine against the malaria parasite would mark a breakthrough in the fight to eradicate the disease. Over time, natural infection elicits a robust immune response against the blood stage of the parasite, providing protection against malaria. In recent years, we have gained valuable insight into the mechanisms by which IgG acts to prevent pathology and inhibit parasite replication, as well as the potential role of immunoglobulin M (IgM) in these processes. Here, we discuss recent advances in our understanding of the mechanisms, acquisition, and maintenance of naturally acquired immunity, and the relevance of these discoveries for the development of a potential vaccine against the blood stage of Plasmodium falciparum.

Host-Malaria Parasite Interactions and Impacts on Mutual Evolution

Malaria is the most deadly parasitic disease, affecting hundreds of millions of people worldwide. Malaria parasites have been associated with their hosts for millions of years. During the long history of host-parasite co-evolution, both parasites and hosts have applied pressure on each other through complex host-parasite molecular interactions. Whereas the hosts activate various immune mechanisms to remove parasites during an infection, the parasites attempt to evade host immunity by diversifying their genome and switching expression of targets of the host immune system. Human intervention to control the disease such as antimalarial drugs and vaccination can greatly alter parasite population dynamics and evolution, particularly the massive applications of antimalarial drugs in recent human history. Vaccination is likely the best method to prevent the disease; however, a partially protective vaccine may have unwanted consequences that require further investigation. Studies of host-parasite interactions and co-evolution will provide important information for designing safe and effective vaccines and for preventing drug resistance. In this essay, we will discuss some interesting molecules involved in host-parasite interactions, including important parasite antigens. We also discuss subjects relevant to drug and vaccine development and some approaches for studying host-parasite interactions.

Genetic diversity of Plasmodium falciparum in Grande Comore Island

Background

Despite several control interventions resulting in a considerable decrease in malaria prevalence in the Union of the Comoros, the disease remains a public health problem with high transmission in Grande Comore compared to neighbouring islands. In this country, only a few studies investigating the genetic diversity of Plasmodium falciparum have been performed so far. For this reason, this study aims to examine the genetic diversity of P. falciparum by studying samples collected in Grande Comore in 2012 and 2013, using merozoite surface protein 1 (msp1), merozoite surface protein 2 (msp2) and single nucleotide polymorphism (SNP) genetic markers.

Methods

A total of 162 positive rapid diagnostic test (RDT) samples from Grande Comore were used to extract parasite DNA. Allelic families K1, Mad20 and RO33 of the msp1 gene as well as allelic families IC3D7 and FC37 of the msp2 gene were determined by using nested PCR. Additionally, 50 out of 151 samples were genotyped to study 24 SNPs by using high resolution melting (HRM).

Results

Two allelic families were predominant, the K1 family of msp1 gene (55%) and the FC27 family of msp2 gene (47.4%). Among 50 samples genotyped for 24 SNPs, 42 (84%) yielded interpretable results. Out of these isolates, 36 (85%) were genetically unique and 6 (15%) grouped into two clusters. The genetic diversity of P. falciparum calculated from msp1 and msp2 genes and SNPs was 0.82 and 0.61, respectively.

Conclusion

In summary, a large genetic diversity of P. falciparum was observed in Grande Comore. This may favour persistence of malaria and might be one of the reasons for the high malaria transmission compared to neighbouring islands. Further surveillance of P. falciparum isolates, mainly through environmental management and vector control, is warranted until complete elimination is attained.

Detection of simple and complex de novo mutations with multiple reference sequences

The characterization of de novo mutations in regions of high sequence and structural diversity from whole-genome sequencing data remains highly challenging. Complex structural variants tend to arise in regions of high repetitiveness and low complexity, challenging both de novo assembly, in which short reads do not capture the long-range context required for resolution, and mapping approaches, in which improper alignment of reads to a reference genome that is highly diverged from that of the sample can lead to false or partial calls. Long-read technologies can potentially solve such problems but are currently unfeasible to use at scale. Here we present Corticall, a graph-based method that combines the advantages of multiple technologies and prior data sources to detect arbitrary classes of genetic variant. We construct multisample, colored de Bruijn graphs from short-read data for all samples, align long-read–derived haplotypes and multiple reference data sources to restore graph connectivity information, and call variants using graph path-finding algorithms and a model for simultaneous alignment and recombination. We validate and evaluate the approach using extensive simulations and use it to characterize the rate and spectrum of de novo mutation events in 119 progeny from four Plasmodium falciparum experimental crosses, using long-read data on the parents to inform reconstructions of the progeny and to detect several known and novel nonallelic homologous recombination events.

Selective whole genome amplification of Plasmodium malariae DNA from clinical samples reveals insights into population structure

The genomic diversity of Plasmodium malariae malaria parasites is understudied, partly because infected individuals tend to present with low parasite densities, leading to difficulties in obtaining sufficient parasite DNA for genome analysis. Selective whole genome amplification (SWGA) increases the relative levels of pathogen DNA in a clinical sample, but has not been adapted for P. malariae parasites. Here we design customized SWGA primers which successfully amplify P. malariae DNA extracted directly from unprocessed clinical blood samples obtained from patients with P. malariae-mono-infections from six countries, and further test the efficacy of SWGA on mixed infections with other Plasmodium spp. SWGA enables the successful whole genome sequencing of samples with low parasite density (i.e. one sample with a parasitaemia of 0.0064% resulted in 44% of the genome covered by ≥ 5 reads), leading to an average 14-fold increase in genome coverage when compared to unamplified samples. We identify a total of 868,476 genome-wide SNPs, of which 194,709 are unique across 18 high-quality isolates. After exclusion of the hypervariable subtelomeric regions, a high-quality core subset of 29,899 unique SNPs is defined. Population genetic analysis suggests that P. malariae parasites display clear geographical separation by continent. Further, SWGA successfully amplifies genetic regions of interest such as orthologs of P. falciparum drug resistance-associated loci (Pfdhfr, Pfdhps, Pfcrt, Pfk13 and Pfmdr1), and several non-synonymous SNPs were detected in these genes. In conclusion, we have established a robust SWGA approach that can assist whole genome sequencing of P. malariae, and thereby facilitate the implementation of much-needed large-scale multi-population genomic studies of this neglected malaria parasite. As demonstrated in other Plasmodia, such genetic diversity studies can provide insights into the biology underlying the disease and inform malaria surveillance and control measures.

Population genetic structure analysis of thrombospondin-related adhesive protein (TRAP) as a vaccine candidate antigen in worldwide Plasmodium falciparum isolates

Antigenic diversity is a major concern in malaria vaccine development that requires to be considered in developing a malaria vaccine. Plasmodium falciparum thrombospondin-related adhesive protein (PfTRAP) is a leading malaria vaccine candidate antigen. In the current study, we investigated the level of genetic diversity and natural selection of pftrap sequences in P. falciparum isolates from Iran (n = 47). The gene diversity of Iranian pftrap sequences was also compared to available global pftrap sequences deposited in the GenBank or PlasmoDB databases (n = 220). Comparison of Iranian PfTRAP sequences with T9/96 reference sequence showed the presence of 35 amino acid changes in 32 positions and a limited variation in repeat sequences, leading to 13 distinct haplotypes. The overall nucleotide diversity (π) for the ectodomain of Iranian pftrap sequences was 0.00444 ± 0.00043, with the highest diversity in Domain IV. Alignment comparison of global PfTRAP sequences with T9/96 reference sequence indicated 96 amino acid replacements as well as extensive variable repeat sequences (9-23 repeats), which led to 192 haplotypes. Among the global isolates, the lowest nucleotide diversity was detected in French Guianan (0.00428 ± 0.00163) and Iranian (0.00444 ± 0.00043) pftrap sequences, and the most variation was observed in domains II and IV in all populations. The dN-dS value displayed the evidence of positive selection due to recombination and immune system pressure. The Fst analysis revealed a gene flow between African populations; however, genetic differentiation observed between Iranian and other populations probably was due to gene flow barriers. Both conserved and variable epitopes were predicted in B- and T-cell epitopes of PfTRAP antigen. The obtained results from this study could be helpful for developing a PfTRAP-based malaria vaccine.

Detecting sequence variants in clinically important protozoan parasites

Second and third generation sequencing methods are crucial for population genetic studies, and variant detection is a popular approach for exploiting this sequence data. While mini- and microsatellites are historically useful markers for studying important Protozoa such as Toxoplasma and Plasmodium spp., detecting non-repetitive variants such as those found in genes can be fundamental to investigating a pathogen's biology. These variants, namely single nucleotide polymorphisms and insertions and deletions, can help elucidate the genetic basis of an organism's pathogenicity, identify selective pressures, and resolve phylogenetic relationships. They also have the added benefit of possessing a comparatively low mutation rate, which contributes to their stability. However, there is a plethora of variant analysis tools with nuanced pipelines and conflicting recommendations for best practise, which can be confounding. This lack of standardisation means that variant analysis requires careful parameter optimisation, an understanding of its limitations, and the availability of high quality data. This review explores the value of variant detection when applied to non-model organisms such as clinically important protozoan pathogens. The limitations of current methods are discussed, including special considerations that require the end-users' attention to ensure that the results generated are reproducible, and the biological conclusions drawn are valid.

A high throughput multi-locus insecticide resistance marker panel for tracking resistance emergence and spread in Anopheles gambiae

The spread of resistance to insecticides in disease-carrying mosquitoes poses a threat to the effectiveness of control programmes, which rely largely on insecticide-based interventions. Monitoring mosquito populations is essential, but obtaining phenotypic measurements of resistance is laborious and error-prone. High-throughput genotyping offers the prospect of quick and repeatable estimates of resistance, while also allowing resistance markers to be tracked and studied. To demonstrate the potential of highly-mulitplexed genotypic screening for measuring resistance-association of mutations and tracking their spread, we developed a panel of 28 known or putative resistance markers in the major malaria vector Anopheles gambiae, which we used to screen mosquitoes from a wide swathe of Sub-Saharan Africa (Burkina Faso, Ghana, Democratic Republic of Congo (DRC) and Kenya). We found resistance association in four markers, including a novel mutation in the detoxification gene Gste2 (Gste2-119V). We also identified a duplication in Gste2 combining a resistance-associated mutation with its wild-type counterpart, potentially alleviating the costs of resistance. Finally, we describe the distribution of the multiple origins of kdr resistance, finding unprecedented diversity in the DRC. This panel represents the first step towards a quantitative genotypic model of insecticide resistance that can be used to predict resistance status in An. gambiae.

Infectivity and Screening of Anti-piperaquine Genes in Mice Infected with Piperaquine-Sensitive and Piperaquine-Resistant Plasmodium berghei

BACKGROUND

Piperaquine (PQ) is one of the major components of artemisinin-based combination therapy for malaria. However, the mechanism of PQ resistance has remained unclear.

METHODS

In this study, we infected mice with PQ-resistant Plasmodium berghei ANKA strain line (PbPQR) or PQ-sensitive P. berghei ANKA strain line (PbPQS) and their survival rates, parasitemia, and spleen sizes were compared. In addition, we constructed genomic DNA subtractive library of spleens from the infected mice, and screened the potential PQ-resistant related genes from genomic DNA of PbPQR line using the representational difference analysis (RDA) method. Clones of the subtractive library were screened by PCR, and related genes were sequenced and analyzed using BLAST software of NCBI.

RESULTS

Compared to PbPQS-infected mice, PbPQR-infected mice survived significantly longer, and had significantly lowered parasitemia rate and significantly increased splenomegaly. Among the total of 502 clones picked, 494 were sequenced and 96 unique PCR fragments were obtained; in which 24 DNA fragments were homologous to chromosomes related to immune function of mice. ORF Finder blasting showed that at the protein level, 26 encoded proteins were homologous to 18 hypothetical PbANKA proteins and 13 encoded proteins were homologous to "ferlin-like protein" family of PbANKA. In addition, there were more immune-related DNA molecules, ubiquitous PbANKA homology at the ORF fragment level, and enriched ferlin-like protein families identified from PbPQR-infected mice than those from PbPQS-infected mice.

CONCLUSION

These findings suggest that PbPQR may induce stronger protective immune response than that of PbPQS in infected mice.

Plasmodium Genomics and Genetics: New Insights into Malaria Pathogenesis, Drug Resistance, Epidemiology, and Evolution

Protozoan Plasmodium parasites are the causative agents of malaria, a deadly disease that continues to afflict hundreds of millions of people every year. Infections with malaria parasites can be asymptomatic, with mild or severe symptoms, or fatal, depending on many factors such as parasite virulence and host immune status. Malaria can be treated with various drugs, with artemisinin-based combination therapies (ACTs) being the first-line choice. Recent advances in genetics and genomics of malaria parasites have contributed greatly to our understanding of parasite population dynamics, transmission, drug responses, and pathogenesis. However, knowledge gaps in parasite biology and host-parasite interactions still remain. Parasites resistant to multiple antimalarial drugs have emerged, while advanced clinical trials have shown partial efficacy for one available vaccine. Here we discuss genetic and genomic studies of Plasmodium biology, host-parasite interactions, population structures, mosquito infectivity, antigenic variation, and targets for treatment and immunization. Knowledge from these studies will advance our understanding of malaria pathogenesis, epidemiology, and evolution and will support work to discover and develop new medicines and vaccines.

Plasticity and genetic variation in traits underpinning asexual replication of the rodent malaria parasite, Plasmodium chabaudi

BACKGROUND

The ability of malaria (Plasmodium) parasites to adjust investment into sexual transmission stages versus asexually replicating stages is well known, but plasticity in other traits underpinning the replication rate of asexual stages in the blood has received less attention. Such traits include burst size (the number of merozoites produced per schizont), the duration of the asexual cycle, and invasion preference for different ages of red blood cell (RBC).

METHODS

Here, plasticity [environment (E) effects] and genetic variation [genotype (G) effects] in traits relating to asexual replication rate are examined for 4 genotypes of the rodent malaria parasite Plasmodium chabaudi. An experiment tested whether asexual dynamics differ between parasites infecting control versus anaemic hosts, and whether variation in replication rate can be explained by differences in burst size, asexual cycle, and invasion rates.

RESULTS

The within-host environment affected each trait to different extents but generally had similar impacts across genotypes. The dynamics of asexual densities exhibited a genotype by environment effect (G×E), in which one of the genotypes increased replication rate more than the others in anaemic hosts. Burst size and cycle duration varied between the genotypes (G), while burst size increased and cycle duration became longer in anaemic hosts (E). Variation in invasion rates of differently aged RBCs was not explained by environmental or genetic effects. Plasticity in burst size and genotype are the only traits making significant contributions to the increase in asexual densities observed in anaemic hosts, together explaining 46.4% of the variation in replication rate.

CONCLUSIONS

That host anaemia induces several species of malaria parasites to alter conversion rate is well documented. Here, previously unknown plasticity in other traits underpinning asexual replication is revealed. These findings contribute to mounting evidence that malaria parasites deploy a suite of sophisticated strategies to maximize fitness by coping with, or exploiting the opportunities provided by, the variable within-host conditions experienced during infections. That genetic variation and genotype by environment interactions also shape these traits highlights their evolutionary potential. Asexual replication rate is a major determinant of virulence and so, understanding the evolution of virulence requires knowledge of the ecological (within-host environment) and genetic drivers of variation among parasites.

Genetic diversity and drug resistance surveillance of Plasmodium falciparum for malaria elimination: is there an ideal tool for resource-limited sub-Saharan Africa?

The intensification of malaria control interventions has resulted in its global decline, but it remains a significant public health burden especially in sub-Saharan Africa (sSA). Knowledge on the parasite diversity, its transmission dynamics, mechanisms of adaptation to environmental and interventional pressures could help refine or develop new control and elimination strategies. Critical to this is the accurate assessment of the parasite’s genetic diversity and monitoring of genetic markers of anti-malarial resistance across all susceptible populations. Such wide molecular surveillance will require selected tools and approaches from a variety of ever evolving advancements in technology and the changing epidemiology of malaria. The choice of an effective approach for specific endemic settings remains challenging, particularly for countries in sSA with limited access to advanced technologies. This article examines the current strategies and tools for Plasmodium falciparum genetic diversity typing and resistance monitoring and proposes how the different tools could be employed in resource-poor settings. Advanced approaches enabling targeted deep sequencing is valued as a sensitive method for assessing drug resistance and parasite diversity but remains out of the reach of most laboratories in sSA due to the high cost of development and maintenance. It is, however, feasible to equip a limited number of laboratories as Centres of Excellence in Africa (CEA), which will receive and process samples from a network of peripheral laboratories in the continent. Cheaper, sensitive and portable real-time PCR methods can be used in peripheral laboratories to pre-screen and select samples for targeted deep sequence or genome wide analyses at these CEAs.

Genomic structure and diversity of Plasmodium falciparum in Southeast Asia reveal recent parasite migration patterns

Estimates of Plasmodium falciparum migration may inform strategies for malaria elimination. Here we elucidate fine-scale parasite population structure and infer recent migration across Southeast Asia using identity-by-descent (IBD) approaches based on genome-wide single nucleotide polymorphisms called in 1722 samples from 54 districts. IBD estimates are consistent with isolation-by-distance. We observe greater sharing of larger IBD segments between artemisinin-resistant parasites versus sensitive parasites, which is consistent with the recent spread of drug resistance. Our IBD analyses reveal actionable patterns, including isolated parasite populations, which may be prioritized for malaria elimination, as well as asymmetrical migration identifying potential sources and sinks of migrating parasites.

Functional Comparison of Blood-Stage Plasmodium falciparum Malaria Vaccine Candidate Antigens

The malaria genome encodes over 5,000 proteins and many of these have also been proposed to be potential vaccine candidates, although few of these have been tested clinically. RH5 is one of the leading blood-stage Plasmodium falciparum malaria vaccine antigens and Phase I/II clinical trials of vaccines containing this antigen are currently underway. Its likely mechanism of action is to elicit antibodies that can neutralize merozoites by blocking their invasion of red blood cells (RBC). However, many other antigens could also elicit neutralizing antibodies against the merozoite, and most of these have never been compared directly to RH5. The objective of this study was to compare a range of blood-stage antigens to RH5, to identify any antigens that outperform or synergize with anti-RH5 antibodies. We selected 55 gene products, covering 15 candidate antigens that have been described in the literature and 40 genes selected on the basis of bioinformatics functional prediction. We were able to make 20 protein-in-adjuvant vaccines from the original selection. Of these, S-antigen and CyRPA robustly elicited antibodies with neutralizing properties. Anti-CyRPA IgG generally showed additive GIA with anti-RH5 IgG, although high levels of anti-CyRPA-specific rabbit polyclonal IgG were required to achieve 50% GIA. Our data suggest that further vaccine antigen screening efforts are required to identify a second merozoite target with similar antibody-susceptibility to RH5.

Comprehensive Review of Human Plasmodium falciparum-Specific CD8+ T Cell Epitopes

Control of malaria is an important global health issue and there is still an urgent need for the development of an effective prophylactic vaccine. Multiple studies have provided strong evidence that Plasmodium falciparum-specific MHC class I-restricted CD8+ T cells are important for sterile protection against Plasmodium falciparum infection. Here, we present an interactive epitope map of all P. falciparum-specific CD8+ T cell epitopes published to date, based on a comprehensive data base (IEDB), and literature search. The majority of the described P. falciparum-specific CD8+ T cells were directed against the antigens CSP, TRAP, AMA1, and LSA1. Notably, most of the epitopes were discovered in vaccine trials conducted with malaria-naïve volunteers. Only few immunological studies of P. falciparum-specific CD8+ T cell epitopes detected in patients suffering from acute malaria or in people living in malaria endemic areas have been published. Further detailed immunological mappings of P. falciparum-specific epitopes of a broader range of P. falciparum proteins in different settings and with different disease status are needed to gain a more comprehensive understanding of the role of CD8+ T cell responses for protection, and to better guide vaccine design and to study their efficacy.

Glucose 6-phosphate dehydrogenase 6-phosphogluconolactonase: characterization of the Plasmodium vivax enzyme and inhibitor studies

BACKGROUND

Since malaria parasites highly depend on ribose 5-phosphate for DNA and RNA synthesis and on NADPH as a source of reducing equivalents, the pentose phosphate pathway (PPP) is considered an excellent anti-malarial drug target. In Plasmodium, a bifunctional enzyme named glucose 6-phosphate dehydrogenase 6-phosphogluconolactonase (GluPho) catalyzes the first two steps of the PPP. PfGluPho has been shown to be essential for the growth of blood stage Plasmodium falciparum parasites.

METHODS

Plasmodium vivax glucose 6-phosphate dehydrogenase (PvG6PD) was cloned, recombinantly produced in Escherichia coli, purified, and characterized via enzyme kinetics and inhibitor studies. The effects of post-translational cysteine modifications were assessed via western blotting and enzyme activity assays. Genetically encoded probes were employed to study the effects of G6PD inhibitors on the cytosolic redox potential of Plasmodium.

RESULTS

Here the recombinant production and characterization of PvG6PD, the C-terminal and NADPH-producing part of PvGluPho, is described. A comparison with PfG6PD (the NADPH-producing part of PfGluPho) indicates that the P. vivax enzyme has higher K_M values for the substrate and cofactor. Like the P. falciparum enzyme, PvG6PD is hardly affected by S-glutathionylation and moderately by S-nitrosation. Since there are several naturally occurring variants of PfGluPho, the impact of these mutations on the kinetic properties of the enzyme was analysed. Notably, in contrast to many human G6PD variants, the mutations resulted in only minor changes in enzyme activity. Moreover, nanomolar IC₅₀ values of several compounds were determined on P. vivax G6PD (including ellagic acid, flavellagic acid, and coruleoellagic acid), inhibitors that had been previously characterized on PfGluPho. ML304, a recently developed PfGluPho inhibitor, was verified to also be active on PvG6PD. Using genetically encoded probes, ML304 was confirmed to disturb the cytosolic glutathione-dependent redox potential of P. falciparum blood stage parasites. Finally, a new series of novel small molecules with the potential to inhibit the falciparum and vivax enzymes were synthesized, resulting in two compounds with nanomolar activity.

CONCLUSION

The characterization of PvG6PD makes this enzyme accessible to further drug discovery activities. In contrast to naturally occurring G6PD variants in the human host that can alter the kinetic properties of the enzyme and thus the redox homeostasis of the cells, the naturally occurring PfGluPho variants studied here are unlikely to have a major impact on the parasites' redox homeostasis. Several classes of inhibitors have been successfully tested and are presently being followed up.

Comprehensive analysis of antibody responses to Plasmodium falciparum erythrocyte membrane protein 1 domains

Acquired antibodies directed towards antigens expressed on the surface of merozoites and infected erythrocytes play an important role in protective immunity to Plasmodium falciparum malaria. P. falciparum erythrocyte membrane protein 1 (PfEMP1), the major parasite component of the infected erythrocyte surface, has been implicated in malaria pathology, parasite sequestration and host immune evasion. However, the extent to which unique PfEMP1 domains interact with host immune response remains largely unknown. In this study, we sought to comprehensively understand the naturally acquired antibody responses targeting different Duffy binding-like (DBL), and Cysteine-rich interdomain region (CIDR) domains in a Ugandan cohort. Consequently, we created a protein library consisting of full-length DBL (n = 163) and CIDR (n = 108) domains derived from 62-var genes based on 3D7 genome. The proteins were expressed by a wheat germ cell-free system; a system that yields plasmodial proteins that are comparatively soluble, intact, biologically active and immunoreactive to human sera. Our findings suggest that all PfEMP1 DBL and CIDR domains, regardless of PfEMP1 group, are targets of naturally acquired immunity. The breadth of the immune response expands with children's age. We concurrently identified 10 DBL and 8 CIDR domains whose antibody responses were associated with reduced risk to symptomatic malaria in the Ugandan children cohort. This study highlights that only a restricted set of specific domains are essential for eliciting naturally acquired protective immunity in malaria. In light of current data, tandem domains in PfEMP1s PF3D7_0700100 and PF3D7_0425800 (DC4) are recommended for extensive evaluation in larger population cohorts to further assess their potential as alternative targets for malaria vaccine development.

Consistent signatures of selection from genomic analysis of pairs of temporal and spatial Plasmodium falciparum populations from The Gambia

Genome sequences of 247 Plasmodium falciparum isolates collected in The Gambia in 2008 and 2014 were analysed to identify changes possibly related to the scale-up of antimalarial interventions that occurred during this period. Overall, there were 15 regions across the genomes with signatures of positive selection. Five of these were sweeps around known drug resistance and antigenic loci. Signatures at antigenic loci such as thrombospodin related adhesive protein (Pftrap) were most frequent in eastern Gambia, where parasite prevalence and transmission remain high. There was a strong temporal differentiation at a non-synonymous SNP in a cysteine desulfarase (Pfnfs) involved in iron-sulphur complex biogenesis. During the 7-year period, the frequency of the lysine variant at codon 65 (Pfnfs-Q65K) increased by 22% (10% to 32%) in the Greater Banjul area. Between 2014 and 2015, the frequency of this variant increased by 6% (20% to 26%) in eastern Gambia. IC₅₀ for lumefantrine was significantly higher in Pfnfs-65K isolates. This is probably the first evidence of directional selection on Pfnfs or linked loci by lumefantrine. Given the declining malaria transmission, the consequent loss of population immunity, and sustained drug pressure, it is important to monitor Gambian P. falciparum populations for further signs of adaptation.

Identifying functional groups among the diverse, recombining antigenic var genes of the malaria parasite Plasmodium falciparum from a local community in Ghana

A challenge in studying diverse multi-copy gene families is deciphering distinct functional types within immense sequence variation. Functional changes can in some cases be tracked through the evolutionary history of a gene family; however phylogenetic approaches are not possible in cases where gene families diversify primarily by recombination. We take a network theoretical approach to functionally classify the highly recombining var antigenic gene family of the malaria parasite Plasmodium falciparum. We sample var DBLα sequence types from a local population in Ghana, and classify 9,276 of these variants into just 48 functional types. Our approach is to first decompose each sequence type into its constituent, recombining parts; we then use a stochastic block model to identify functional groups among the parts; finally, we classify the sequence types based on which functional groups they contain. This method for functional classification does not rely on an inferred phylogenetic history, nor does it rely on inferring function based on conserved sequence features. Instead, it infers functional similarity among recombining parts based on the sharing of similar co-occurrence interactions with other parts. This method can therefore group sequences that have undetectable sequence homology or even distinct origination. Describing these 48 var functional types allows us to simplify the antigenic diversity within our dataset by over two orders of magnitude. We consider how the var functional types are distributed in isolates, and find a nonrandom pattern reflecting that common var functional types are non-randomly distinct from one another in terms of their functional composition. The coarse-graining of var gene diversity into biologically meaningful functional groups has important implications for understanding the disease ecology and evolution of this system, as well as for designing effective epidemiological monitoring and intervention.

Gene copy number variation in natural populations of Plasmodium falciparum in Eastern Africa

BACKGROUND

Gene copy number variants (CNVs), which consist of deletions and amplifications of single or sets of contiguous genes, contribute to the great diversity in the Plasmodium falciparum genome. In vitro studies in the laboratory have revealed their important role in parasite fitness phenotypes such as red cell invasion, transmissibility and cytoadherence. Studies of natural parasite populations indicate that CNVs are also common in the field and thus may facilitate adaptation of the parasite to its local environment.

RESULTS

In a survey of 183 fresh field isolates from three populations in Eastern Africa with different malaria transmission intensities, we identified 94 CNV loci using microarrays. All CNVs had low population frequencies (minor allele frequency < 5%) but each parasite isolate carried an average of 8 CNVs. Nine CNVs showed high levels of population differentiation (F_ST > 0.3) and nine exhibited significant clines in population frequency across a gradient in transmission intensity. The clearest example of this was a large deletion on chromosome 9 previously reported only in laboratory-adapted isolates. This deletion was present in 33% of isolates from a population with low and highly seasonal malaria transmission, and in < 9% of isolates from populations with higher transmission. Subsets of CNVs were strongly correlated in their population frequencies, implying co-selection.

CONCLUSIONS

These results support the hypothesis that CNVs are the target of selection in natural populations of P. falciparum. Their environment-specific patterns observed here imply an important role for them in conferring adaptability to the parasite thus enabling it to persist in its highly diverse ecological environment.

Rethinking the extrinsic incubation period of malaria parasites

The time it takes for malaria parasites to develop within a mosquito, and become transmissible, is known as the extrinsic incubation period, or EIP. EIP is a key parameter influencing transmission intensity as it combines with mosquito mortality rate and competence to determine the number of mosquitoes that ultimately become infectious. In spite of its epidemiological significance, data on EIP are scant. Current approaches to estimate EIP are largely based on temperature-dependent models developed from data collected on parasite development within a single mosquito species in the 1930s. These models assume that the only factor affecting EIP is mean environmental temperature. Here, we review evidence to suggest that in addition to mean temperature, EIP is likely influenced by genetic diversity of the vector, diversity of the parasite, and variation in a range of biotic and abiotic factors that affect mosquito condition. We further demonstrate that the classic approach of measuring EIP as the time at which mosquitoes first become infectious likely misrepresents EIP for a mosquito population. We argue for a better understanding of EIP to improve models of transmission, refine predictions of the possible impacts of climate change, and determine the potential evolutionary responses of malaria parasites to current and future mosquito control tools.

Signatures of competition and strain structure within the major blood-stage antigen of Plasmodium falciparum in a local community in Ghana

The concept of niche partitioning has received considerable theoretical attention at the interface of ecology and evolution of infectious diseases. Strain theory postulates that pathogen populations can be structured into distinct nonoverlapping strains by frequency-dependent selection in response to intraspecific competition for host immune space. The malaria parasite Plasmodium falciparum presents an opportunity to investigate this phenomenon in nature, under conditions of high recombination rate and extensive antigenic diversity. The parasite's major blood-stage antigen, Pf EMP1, is encoded by the hyperdiverse var genes. With a dataset that includes thousands of var DBLα sequence types sampled from asymptomatic cases within an area of high endemicity in Ghana, we address how var diversity is distributed within isolates and compare this to the distribution of microsatellite allelic diversity within isolates to test whether antigenic and neutral regions of the genome are structured differently. With respect to var DBLα sequence types, we find that on average isolates exhibit significantly lower overlap than expected randomly, but that there also exists frequent pairs of isolates that are highly related. Furthermore, the linkage network of var DBLα sequence types reveals a pattern of nonrandom modularity unique to these antigenic genes, and we find that modules of highly linked DBLα types are not explainable by neutral forces related to var recombination constraints, microsatellite diversity, sampling location, host age, or multiplicity of infection. These findings of reduced overlap and modularity among the var antigenic genes are consistent with a role for immune selection as proposed by strain theory. Identifying the evolutionary and ecological dynamics that are responsible for the nonrandom structure in P. falciparum antigenic diversity is important for designing effective intervention in endemic areas.

Malaria Pathogenesis

In the mosquito-human life cycle, the six species of malaria parasites infecting humans (Plasmodium falciparum, Plasmodium vivax, Plasmodium ovale wallickeri, Plasmodium ovale curtisi, Plasmodium malariae, and Plasmodium knowlesi) undergo 10 or more morphological states, replicate from single to 10,000+ cells, and vary in total population from one to many more than 10⁶ organisms. In the human host, only a small number of these morphological stages lead to clinical disease and the vast majority of all malaria-infected patients in the world produce few (if any) symptoms in the human. Human clinical disease (e.g., fever, anemia, coma) is the result of the parasite preprogrammed biology in concert with the human pathophysiological response. Caveats and corollaries that add variation to this host-parasite interaction include parasite genetic diversity of key proteins, coinfections, comorbidities, delays in treatment, human polymorphisms, and environmental determinants.

Quantitative Proteomic Profiling Reveals Novel Plasmodium falciparum Surface Antigens and Possible Vaccine Candidates

Despite recent efforts toward control and elimination, malaria remains a major public health problem worldwide. Plasmodium falciparum resistance against artemisinin, used in front line combination drugs, is on the rise, and the only approved vaccine shows limited efficacy. Combinations of novel and tailored drug and vaccine interventions are required to maintain the momentum of the current malaria elimination program. Current evidence suggests that strain-transcendent protection against malaria infection can be achieved using whole organism vaccination or with a polyvalent vaccine covering multiple antigens or epitopes. These approaches have been successfully applied to the human-infective sporozoite stage. Both systemic and tissue-specific pathology during infection with the human malaria parasite P. falciparum is caused by asexual blood stages. Tissue tropism and vascular sequestration are the result of specific binding interactions between antigens on the parasite-infected red blood cell (pRBC) surface and endothelial receptors. The major surface antigen and parasite ligand binding to endothelial receptors, PfEMP1 is encoded by about 60 variants per genome and shows high sequence diversity across strains. Apart from PfEMP1 and three additional variant surface antigen families RIFIN, STEVOR, and SURFIN, systematic analysis of the infected red blood cell surface is lacking. Here we present the most comprehensive proteomic investigation of the parasitized red blood cell surface so far. Apart from the known variant surface antigens, we identified a set of putative single copy surface antigens with low sequence diversity, several of which are validated in a series of complementary experiments. Further functional and immunological investigation is underway to test these novel P. falciparum blood stage proteins as possible vaccine candidates.

Synergistic malaria vaccine combinations identified by systematic antigen screening

Malaria still kills hundreds of thousands of children each year. Malaria vaccine development is complicated by high levels of parasite genetic diversity, which makes single target vaccines vulnerable to the development of variant-specific immunity. To overcome this hurdle, we systematically screened a panel of 29 blood-stage antigens from the most deadly human malaria parasite, Plasmodium falciparum. We identified several targets that were able to inhibit erythrocyte invasion in two genetically diverse strains. Testing these targets in combination identified several pairs that blocked invasion more effectively in combination than in isolation. Video microscopy and studies of natural immune responses to malaria in patients suggest that targeting multiple steps in invasion is more likely to produce a synergistic vaccine response.

A highly effective vaccine would be a valuable weapon in the drive toward malaria elimination. No such vaccine currently exists, and only a handful of the hundreds of potential candidates in the parasite genome have been evaluated. In this study, we systematically evaluated 29 antigens likely to be involved in erythrocyte invasion, an essential developmental stage during which the malaria parasite is vulnerable to antibody-mediated inhibition. Testing antigens alone and in combination identified several strain-transcending targets that had synergistic combinatorial effects in vitro, while studies in an endemic population revealed that combinations of the same antigens were associated with protection from febrile malaria. Video microscopy established that the most effective combinations targeted multiple discrete stages of invasion, suggesting a mechanistic explanation for synergy. Overall, this study both identifies specific antigen combinations for high-priority clinical testing and establishes a generalizable approach that is more likely to produce effective vaccines.

Extreme mutation bias and high AT content in Plasmodium falciparum

For reasons that remain unknown, the Plasmodium falciparum genome has an exceptionally high AT content compared to other Plasmodium species and eukaryotes in general - nearly 80% in coding regions and approaching 90% in non-coding regions. Here, we examine how this phenomenon relates to genome-wide patterns of de novo mutation. Mutation accumulation experiments were performed by sequential cloning of six P. falciparum isolates growing in human erythrocytes in vitro for 4 years, with 279 clones sampled for whole genome sequencing at different time points. Genome sequence analysis of these samples revealed a significant excess of G:C to A:T transitions compared to other types of nucleotide substitution, which would naturally cause AT content to equilibrate close to the level seen across the P. falciparum reference genome (80.6% AT). These data also uncover an extremely high rate of small indel mutation relative to other species, primarily associated with repetitive AT-rich sequences, in addition to larger-scale structural rearrangements focused in antigen-coding var genes. In conclusion, high AT content in P. falciparum is driven by a systematic mutational bias and ultimately leads to an unusual level of microstructural plasticity, raising the question of whether this contributes to adaptive evolution.

Population Genomics of Paramecium Species

Population-genomic analyses are essential to understanding factors shaping genomic variation and lineage-specific sequence constraints. The dearth of such analyses for unicellular eukaryotes prompted us to assess genomic variation in Paramecium, one of the most well-studied ciliate genera. The Paramecium aurelia complex consists of ∼15 morphologically indistinguishable species that diverged subsequent to two rounds of whole-genome duplications (WGDs, as long as 320 MYA) and possess extremely streamlined genomes. We examine patterns of both nuclear and mitochondrial polymorphism, by sequencing whole genomes of 10-13 worldwide isolates of each of three species belonging to the P. aurelia complex: P. tetraurelia, P. biaurelia, P. sexaurelia, as well as two outgroup species that do not share the WGDs: P. caudatum and P. multimicronucleatum. An apparent absence of global geographic population structure suggests continuous or recent dispersal of Paramecium over long distances. Intergenic regions are highly constrained relative to coding sequences, especially in P. caudatum and P. multimicronucleatum that have shorter intergenic distances. Sequence diversity and divergence are reduced up to ∼100-150 bp both upstream and downstream of genes, suggesting strong constraints imposed by the presence of densely packed regulatory modules. In addition, comparison of sequence variation at non-synonymous and synonymous sites suggests similar recent selective pressures on paralogs within and orthologs across the deeply diverging species. This study presents the first genome-wide population-genomic analysis in ciliates and provides a valuable resource for future studies in evolutionary and functional genetics in Paramecium.

Malaria Genomics in the Era of Eradication

The first reference genome assembly for the Plasmodium falciparum malaria parasite was completed over a decade ago, and the impact of this and other genomic resources on malaria research has been significant. Genomic resources for other malaria parasites are being established, even as P. falciparum continues to be the focus of development of new genomic methods and applications. Here we review the impact and applications of genomic data on malaria research, and discuss future needs and directions as genomic data generation becomes less expensive and more decentralized. Specifically, we focus on how population genomic strategies can be utilized to advance the malaria eradication agenda.