Allele frequency-based and polymorphism-versus-divergence indices of balancing selection in a new filtered set of polymorphic genes in Plasmodium falciparum

Expression of the MSPDBL2 antigen in a discrete subset of Plasmodium falciparum schizonts is regulated by GDV1 but may not be linked to sexual commitment

The Plasmodium falciparum merozoite surface protein MSPDBL2 is a polymorphic antigen targeted by acquired immune responses, and normally expressed in only a minority of mature schizonts. The potential relationship of MSPDBL2 to sexual commitment is examined, as variable mspdbl2 transcript levels and proportions of MSPDBL2-positive mature schizonts in clinical isolates have previously correlated with levels of many sexual stage parasite gene transcripts, although not with the master regulator ap2-g. It is demonstrated that conditional overexpression of the gametocyte development protein GDV1, which promotes sexual commitment, also substantially increases the proportion of MSPDBL2-positive schizonts in culture. Conversely, truncation of the gdv1 gene is shown to prevent any expression of MSPDBL2. However, across diverse P. falciparum cultured lines, the variable proportions of MSPDBL2 positivity in schizonts do not correlate significantly with variable gametocyte conversion rates, indicating it is not involved in sexual commitment. Confirming this, examining a line with endogenous hemagglutinin-tagged AP2-G showed that the individual schizonts expressing MSPDBL2 are mostly different from those expressing AP2-G. Using a selection-linked integration system, modified P. falciparum lines were engineered to express an intact or disrupted version of MSPDBL2, showing the protein is not required for sexual commitment or early gametocyte development. Asexual parasite multiplication rates were also not affected by expression of either intact or disrupted MSPDBL2 in a majority of schizonts. Occurring alongside sexual commitment, the role of the discrete MSPDBL2-positive schizont subpopulation requires further investigation in natural infections where it is under immune selection.

IMPORTANCE

Malaria parasites in the blood are remarkably variable, able to switch antigenic targets so they may survive within humans who have already developed specific immune responses. This is one of the challenges in developing vaccines against malaria. MSPDBL2 is a target of naturally acquired immunity expressed in minority proportions of schizonts, the end stages of each 2-day replication cycle in red blood cells which contain merozoites prepared to invade new red blood cells. Results show that the proportion of schizonts expressing MSPDBL2 is positively controlled by the expression of the regulatory gametocyte development protein GDV1. It was previously known that expression of GDV1 leads to increased expression of AP2-G which causes parasites to switch to sexual development, so a surprising finding here is that MSPDBL2-positive parasites are mostly distinct from those that express AP2-G. This discrete antigenic subpopulation of mostly asexual parasites is regulated alongside sexually committed parasites, potentially enabling survival under stress conditions.

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.

Complex patterns shape immune genes diversity during invasion of common raccoon in Europe - Selection in action despite genetic drift

Rapid adaptation is common in invasive populations and is crucial to their long-term success. The primary target of selection in the invasive species' new range is standing genetic variation. Therefore, genetic drift and natural selection acting on existing variation are key evolutionary processes through which invaders will evolve over a short timescale. In this study, we used the case of the raccoon Procyon lotor invasion in Europe to identify the forces shaping the diversity of immune genes during invasion. The genes involved in the defence against infection should be under intense selection pressure in the invasive range where novel pathogens are expected to occur. To disentangle the selective and demographic processes shaping the adaptive immune diversity of its invasive and expanding populations, we have developed species-specific single-nucleotide polymorphism markers located in the coding regions of targeted immune-related genes. We characterised the genetic diversity of 110 functionally important immune genes in two invasive and one native raccoon genetic clusters, each presenting a different demographic history. Despite the strong effect of demographic processes in the invasive clusters, we detected a subset of genes exhibiting the diversity pattern suggestive of selection. The most likely process shaping the variation in those genes was balancing selection. The selected genes belong to toll-like receptors and cytokine-related genes. Our results suggest that the prevalence of selection depends on the level of diversity, that is - less genetically diverse invasive population from the Czech Republic displayed fewer signs of selection. Our results highlight the role of standing genetic variation in adapting to new environment. Understanding the evolutionary mechanisms behind invasion success would enable predicting how populations may respond to environmental change.

Highly Variable Expression of Merozoite Surface Protein MSPDBL2 in Diverse Plasmodium falciparum Clinical Isolates and Transcriptome Scans for Correlating Genes

The merozoite surface protein MSPDBL2 of Plasmodium falciparum is under strong balancing selection and is a target of naturally acquired antibodies. Remarkably, MSPDBL2 is expressed in only a minority of mature schizonts of any cultured parasite line, and mspdbl2 gene transcription increases in response to overexpression of the gametocyte development inducer GDV1, so it is important to understand its natural expression. Here, MSPDBL2 in mature schizonts was analyzed in the first ex vivo culture cycle of 96 clinical isolates from 4 populations with various levels of infection endemicity in different West African countries, by immunofluorescence microscopy with antibodies against a conserved region of the protein. In most isolates, less than 1% of mature schizonts were positive for MSPDBL2, but the frequency distribution was highly skewed, as nine isolates had more than 3% schizonts positive and one had 73% positive. To investigate whether the expression of other gene loci correlated with MSPDBL2 expression, whole-transcriptome sequencing was performed on schizont-enriched material from 17 of the isolates with a wide range of proportions of schizonts positive. Transcripts of particular genes were highly significantly positively correlated with MSPDBL2 positivity in schizonts as well as with mspdbl2 gene transcript levels, showing overrepresentation of genes implicated previously as involved in gametocytogenesis but not including the gametocytogenesis master regulator ap2-g. Single-cell transcriptome analysis of a laboratory-adapted clone showed that most individual parasites expressing mspdbl2 did not express ap2-g, consistent with MSPDBL2 marking a developmental subpopulation that is distinct but likely to co-occur alongside sexual commitment. IMPORTANCE These findings contribute to understanding malaria parasite antigenic and developmental variation, focusing on the merozoite surface protein encoded by the single locus under strongest balancing selection. Analyzing the initial ex vivo generation of parasites grown from a wide sample of clinical infections, we show a unique and highly skewed pattern of natural expression frequencies of MSPDBL2, distinct from that of any other antigen. Bulk transcriptome analysis of a range of clinical isolates showed significant overrepresentation of sexual development genes among those positively correlated with MSPDBL2 protein and mspdbl2 gene expression, indicating the MSPDBL2-positive subpopulation to be often coincident with parasites developing sexually in preparation for transmission. Single-cell transcriptome data confirm the absence of a direct correlation with the ap2-g master regulator of sexual development, indicating that the MSPDBL2-positive subpopulation has a separate function in asexual survival and replication under conditions that promote terminal sexual differentiation.

Gramtools enables multiscale variation analysis with genome graphs

Genome graphs allow very general representations of genetic variation; depending on the model and implementation, variation at different length-scales (single nucleotide polymorphisms (SNPs), structural variants) and on different sequence backgrounds can be incorporated with different levels of transparency. We implement a model which handles this multiscale variation and develop a JSON extension of VCF (jVCF) allowing for variant calls on multiple references, both implemented in our software gramtools. We find gramtools outperforms existing methods for genotyping SNPs overlapping large deletions in M. tuberculosis and is able to genotype on multiple alternate backgrounds in P. falciparum, revealing previously hidden recombination.

Immune selection suppresses the emergence of drug resistance in malaria parasites but facilitates its spread

Although drug resistance in Plasmodium falciparum typically evolves in regions of low transmission, resistance spreads readily following introduction to regions with a heavier disease burden. This suggests that the origin and the spread of resistance are governed by different processes, and that high transmission intensity specifically impedes the origin. Factors associated with high transmission, such as highly immune hosts and competition within genetically diverse infections, are associated with suppression of resistant lineages within hosts. However, interactions between these factors have rarely been investigated and the specific relationship between adaptive immunity and selection for resistance has not been explored. Here, we developed a multiscale, agent-based model of Plasmodium parasites, hosts, and vectors to examine how host and parasite dynamics shape the evolution of resistance in populations with different transmission intensities. We found that selection for antigenic novelty (“immune selection”) suppressed the evolution of resistance in high transmission settings. We show that high levels of population immunity increased the strength of immune selection relative to selection for resistance. As a result, immune selection delayed the evolution of resistance in high transmission populations by allowing novel, sensitive lineages to remain in circulation at the expense of the spread of a resistant lineage.

In contrast, in low transmission settings, we observed that resistant strains were able to sweep to high population prevalence without interference. Additionally, we found that the relationship between immune selection and resistance changed when resistance was widespread. Once resistance was common enough to be found on many antigenic backgrounds, immune selection stably maintained resistant parasites in the population by allowing them to proliferate, even in untreated hosts, when resistance was linked to a novel epitope. Our results suggest that immune selection plays a role in the global pattern of resistance evolution.

Drug resistance in the malaria parasite, Plasmodium falciparum, presents an ongoing public health challenge, but aspects of its evolution are poorly understood. Although antimalarial resistance is common worldwide, it can typically be traced to just a handful of evolutionary origins. Counterintuitively, although Sub Saharan Africa bears 90% of the global malaria burden, resistance typically originates in regions where transmission intensity is low. In high transmission regions, infections are genetically diverse, and hosts have significant standing adaptive immunity, both of which are known to suppress the frequency of resistance within infections. However, interactions between immune-driven selection, transmission intensity, and resistance have not been investigated. Using a multiscale, agent-based model, we found that high transmission intensity slowed the evolution of resistance via its effect on host population immunity. High host immunity strengthened selection for antigenic novelty, interfering with selection for resistance and allowing sensitive lineages to suppress resistant lineages in untreated hosts. However, once resistance was common in the circulating parasite population, immune selection maintained it in the population at a high prevalence. Our findings provide a novel explanation for observations about the origin of resistance and suggest that adaptive immunity is a critical component of selection.

Genome-Wide Association Analysis Identifies a Genetic Basis of Infectivity in a Model Bacterial Pathogen

Knowledge of the genetic architecture of pathogen infectivity and host resistance is essential for a mechanistic understanding of coevolutionary processes, yet the genetic basis of these interacting traits remains unknown for most host-pathogen systems. We used a comparative genomic approach to explore the genetic basis of infectivity in Pasteuria ramosa, a Gram-positive bacterial pathogen of planktonic crustaceans that has been established as a model for studies of Red Queen host-pathogen coevolution. We sequenced the genomes of a geographically, phenotypically, and genetically diverse collection of P. ramosa strains and performed a genome-wide association study to identify genetic correlates of infection phenotype. We found multiple polymorphisms within a single gene, Pcl7, that correlate perfectly with one common and widespread infection phenotype. We then confirmed this perfect association via Sanger sequencing in a large and diverse sample set of P. ramosa clones. Pcl7 codes for a collagen-like protein, a class of adhesion proteins known or suspected to be involved in the infection mechanisms of a number of important bacterial pathogens. Consistent with expectations under Red Queen coevolution, sequence variation of Pcl7 shows evidence of balancing selection, including extraordinarily high diversity and absence of geographic structure. Based on structural homology with a collagen-like protein of Bacillus anthracis, we propose a hypothesis for the structure of Pcl7 and the physical location of the phenotype-associated polymorphisms. Our results offer strong evidence for a gene governing infectivity and provide a molecular basis for further study of Red Queen dynamics in this model host-pathogen system.

Plasmodium falciparum Malaria Parasites in Ghana Show Signatures of Balancing Selection at Artemisinin Resistance Predisposing Background Genes

Sub-Saharan Africa is courting the risk of artemisinin resistance (ARTr) emerging in Plasmodium falciparum malaria parasites. Current molecular surveillance efforts for ARTr have been built on the utility of P. falciparum kelch13 (pfk13) validated molecular markers. However, whether these molecular markers will serve the purpose of early detection of artemisinin-resistant parasites in Ghana is hinged on a pfk13 dependent evolution. Here, we tested the hypothesis that the background pfk13 genome may be present before the pfk13 ARTr-conferring variant(s) is selected and that signatures of balancing selection on these genomic loci may serve as an early warning signal of ARTr. We analyzed 12 198 single nucleotide polymorphisms (SNPs) in Ghanaian clinical isolates in the Pf3K MalariaGEN dataset that passed a stringent filtering regimen. We identified signatures of balancing selection in 2 genes (phosphatidylinositol 4-kinase and chloroquine resistance transporter) previously reported as background loci for ARTr. These genes showed statistically significant and high positive values for Tajima's D, Fu and Li's F, and Fu and Li's D. This indicates that the biodiversity required to establish a pfk13 background genome may have been primed in clinical isolates of P. falciparum from Ghana as of 2010. Despite the absence of ARTr in Ghana to date, our finding supports the current use of pfk13 for molecular surveillance of ARTr in Ghana and highlights the potential utility of monitoring malaria parasite populations for balancing selection in ARTr precursor background genes as early warning molecular signatures for the emergence of ARTr.

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.

Host-parasite co-evolution and its genomic signature

Studies in diverse biological systems have indicated that host-parasite co-evolution is responsible for the extraordinary genetic diversity seen in some genomic regions, such as major histocompatibility (MHC) genes in jawed vertebrates and resistance genes in plants. This diversity is believed to evolve under balancing selection on hosts by parasites. However, the mechanisms that link the genomic signatures in these regions to the underlying co-evolutionary process are only slowly emerging. We still lack a clear picture of the co-evolutionary concepts and of the genetic basis of the co-evolving phenotypic traits in the interacting antagonists. Emerging genomic tools that provide new options for identifying underlying genes will contribute to a fuller understanding of the co-evolutionary process.

Few Plasmodium falciparum merozoite ligand and erythrocyte receptor pairs show evidence of balancing selection

Erythrocyte surface proteins have been identified as receptors of Plasmodium falciparum merozoite proteins. The ligand-receptor interactions enable the parasite to invade human erythrocytes, initiating the clinical symptoms of malaria. These interactions are likely to have had an evolutionary impact on the genes that encode the ligand and receptor proteins. We used sequence data from Kilifi, Kenya to detect departures from neutrality in a paired analysis of P. falciparum merozoite ligands and their erythrocyte receptor genes from the same population. We genotyped parasite and human DNA obtained from 93 individuals with severe malaria. We examined six merozoite ligands EBA175, EBL1, EBA140, MSP1, Rh4 and Rh5, and their corresponding erythrocyte receptors, glycophorin (Gyp) A, GypB, GypC, band 3, complement receptor (CR) 1 and basigin, focusing on the regions involved in the ligand-receptor interactions. Positive Tajima's D values (>1) were observed only in the MSP1 C-terminal region and EBA175 region II, while negative values (<-1) were observed in EBL-1 region II, Rh4, basigin exons 3 and 5, CR1 exon 5, Gyp B exons 2, 3 and 4 and Gyp C exon 2. Additionally, ebl-1 region II and basigin exon 3 showed extreme negative values in all three tests, Tajima's D, Fu & Li D* and F*, ≤ - 2. A large majority of the erythrocyte receptor and merozoite genes have a negative Tajima's D even when compared with previously published whole genome data. Thus, highlighting EBA175 region II and MSP1-33, as outlier genes with a positive Tajima's D (>1). Both these genes contain multiple polymorphisms, which in the case of EBA175 may counteract receptor polymorphisms and/or evade host immune responses and in MSP1 the polymorphisms may primarily evade host immune responses.

Schizont transcriptome variation among clinical isolates and laboratory-adapted clones of the malaria parasite Plasmodium falciparum

BACKGROUND

Malaria parasites are genetically polymorphic and phenotypically plastic. In studying transcriptome variation among parasites from different infections, it is challenging to overcome potentially confounding technical and biological variation between samples. We investigate variation in the major human parasite Plasmodium falciparum, generating RNA-seq data on multiple independent replicate sample preparations of merozoite-containing intra-erythrocytic schizonts from a panel of clinical isolates and from long-term laboratory-adapted clones, with a goal of robustly identifying differentially expressed genes.

RESULTS

Analysis of biological sample replicates shows that increased numbers improve the true discovery rate of differentially expressed genes, and that six independent replicates of each parasite line allowed identification of most differences that could be detected with larger numbers. For highly expressed genes, focusing on the top quartile at schizont stages, there was more power to detect differences. Comparing cultured clinical isolates and laboratory-adapted clones, genes more highly expressed in the laboratory-adapted clones include those encoding an AP2 transcription factor (PF3D7_0420300), a ubiquitin-binding protein and two putative methyl transferases. In contrast, higher expression in clinical isolates was seen for the merozoite surface protein gene dblmsp2, proposed to be a marker of schizonts forming merozoites committed to sexual differentiation. Variable expression was extremely strongly, but not exclusively, associated with genes known to be targeted by Heterochromatin Protein 1. Clinical isolates show variable expression of several known merozoite invasion ligands, as well as other genes for which new RT-qPCR assays validate the quantitation and allow characterisation in samples with more limited material. Expression levels of these genes vary among schizont preparations of different clinical isolates in the first ex vivo cycle in patient erythrocytes, but mean levels are similar to those in continuously cultured clinical isolates.

CONCLUSIONS

Analysis of multiple biological sample replicates greatly improves identification of genes variably expressed between different cultured parasite lines. Clinical isolates recently established in culture show differences from long-term adapted clones in transcript levels of particular genes, and are suitable for analyses requiring biological replicates to understand parasite phenotypes and variable expression likely to be relevant in nature.

Molecular epidemiology and evolution of drug-resistant genes in the malaria parasite Plasmodium falciparum in southwestern Nigeria

Malaria is an age-old disease of human kind living in the tropical and sub-tropical regions of the globe, with Africa contributing the highest incidence of morbidity and mortality. Among many hurdles, evolution and spread of drug-resistant Plasmodium falciparum parasites constitute major challenges to malaria control and elimination. Information on molecular epidemiology and pattern of evolution of genes conferring resistance to different antimalarials are needed to track the route of the spread of resistant parasites and also to inform if the drug-resistant genes are adapted in the population following the Darwinian model of evolution. In the present study, we have followed molecular methods to detect both the known and emerging mutations in three genes (Pfcrt, Pfdhfr and Pfdhps) of P. falciparum conferring resistance to chloroquine and sulfadoxine-pyrimethamine from two different states (Edo: meso-endemic and Lagos: hypo-endemic) in southwestern Nigeria. High diversities in haplotypes and nucleotides in genes responsible for chloroquine (Pfcrt) and sulfadoxine (Pfdhps) resistance are recorded. About 96% of Pfdhfr and Pfdhps gene in both the meso- and hypo- endemic areas were mutant type, followed by 61% in Pfcrt gene. Many unique haplotypes of Pfdhps and Pfcrt were found to be segregated in these two populations. One particular mutant haplotype of Pfdhfr (AIRNI) was found to be in very high frequency in both Lagos and Edo. While the net haplotype diversity was highest in Pfdhps (0.81 in Lagos, 0.87 in Edo), followed by Pfcrt (0.69 in Lagos, 0.65 in Edo); highest number of haplotype was found in Pfdhps with 13 distinct haplotypes, followed by seven in Pfcrt and four in Pfdhfr gene. Moreover, detection of strong linkage among mutations of Pfcrt and Pfdhfr and feeble evidence for balancing selection in Pfdhps are indicative of evolutionary potential of mutation in genes responsible for drug resistance in Nigerian populations of P. falciparum.

An expanded global inventory of allelic variation in the most extremely polymorphic region of Plasmodium falciparum merozoite surface protein 1 provided by short read sequence data

Background

Within Plasmodium falciparum merozoite surface protein 1 (MSP1), the N-terminal block 2 region is a highly polymorphic target of naturally acquired antibody responses. The antigenic diversity is determined by complex repeat sequences as well as non-repeat sequences, grouping into three major allelic types that appear to be maintained within populations by natural selection. Within these major types, many distinct allelic sequences have been described in different studies, but the extent and significance of the diversity remains unresolved.

Methods

To survey the diversity more extensively, block 2 allelic sequences in the msp1 gene were characterized in 2400 P. falciparum infection isolates with whole genome short read sequence data available from the Pf3K project, and compared with the data from previous studies.

Results

Mapping the short read sequence data in the 2400 isolates to a reference library of msp1 block 2 allelic sequences yielded 3815 allele scores at the level of major allelic family types, with 46% of isolates containing two or more of these major types. Overall frequencies were similar to those previously reported in other samples with different methods, the K1-like allelic type being most common in Africa, MAD20-like most common in Southeast Asia, and RO33-like being the third most abundant type in each continent. The rare MR type, formed by recombination between MAD20-like and RO33-like alleles, was only seen in Africa and very rarely in the Indian subcontinent but not in Southeast Asia. A combination of mapped short read assembly approaches enabled 1522 complete msp1 block 2 sequences to be determined, among which there were 363 different allele sequences, of which 246 have not been described previously. In these data, the K1-like msp1 block 2 alleles are most diverse and encode 225 distinct amino acid sequences, compared with 123 different MAD20-like, 9 RO33-like and 6 MR type sequences. Within each of the major types, the different allelic sequences show highly skewed geographical distributions, with most of the more common sequences being detected in either Africa or Asia, but not in both.

Conclusions

Allelic sequences of this extremely polymorphic locus have been derived from whole genome short read sequence data by mapping to a reference library followed by assembly of mapped reads. The catalogue of sequence variation has been greatly expanded, so that there are now more than 500 different msp1 block 2 allelic sequences described. This provides an extensive reference for molecular epidemiological genotyping and sequencing studies, and potentially for design of a multi-allelic vaccine.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2475-2) contains supplementary material, which is available to authorized users.

Antibodies in children with malaria to PfEMP1, RIFIN and SURFIN expressed at the Plasmodium falciparum parasitized red blood cell surface

Naturally acquired antibodies to proteins expressed on the Plasmodium falciparum parasitized red blood cell (pRBC) surface steer the course of a malaria infection by reducing sequestration and stimulating phagocytosis of pRBC. Here we have studied a selection of proteins representing three different parasite gene families employing a well-characterized parasite with a severe malaria phenotype (FCR3S1.2). The presence of naturally acquired antibodies, impact on rosetting rate, surface reactivity and opsonization for phagocytosis in relation to different blood groups of the ABO system were assessed in a set of sera from children with mild or complicated malaria from an endemic area. We show that the naturally acquired immune responses, developed during malaria natural infection, have limited access to the pRBCs inside a blood group A rosette. The data also indicate that SURFIN_4.2 may have a function at the pRBC surface, particularly during rosette formation, this role however needs to be further validated. Our results also indicate epitopes differentially recognized by rosette-disrupting antibodies on a peptide array. Antibodies towards parasite-derived proteins such as PfEMP1, RIFIN and SURFIN in combination with host factors, essentially the ABO blood group of a malaria patient, are suggested to determine the outcome of a malaria infection.

Vaccine candidate discovery for the next generation of malaria vaccines

Although epidemiological observations, IgG passive transfer studies and experimental infections in humans all support the feasibility of developing highly effective malaria vaccines, the precise antigens that induce protective immunity remain uncertain. Here, we review the methodologies applied to vaccine candidate discovery for Plasmodium falciparum malaria from the pre- to post-genomic era. Probing of genomic and cDNA libraries with antibodies of defined specificities or functional activity predominated the former, whereas reverse vaccinology encompassing high throughput in silico analyses of genomic, transcriptomic or proteomic parasite data sets is the mainstay of the latter. Antibody-guided vaccine design spanned both eras but currently benefits from technological advances facilitating high-throughput screening and downstream applications. We make the case that although we have exponentially increased our ability to identify numerous potential vaccine candidates in a relatively short space of time, a significant bottleneck remains in their validation and prioritization for evaluation in clinical trials. Longitudinal cohort studies provide supportive evidence but results are often conflicting between studies. Demonstration of antigen-specific antibody function is valuable but the relative importance of one mechanism over another with regards to protection remains undetermined. Animal models offer useful insights but may not accurately reflect human disease. Challenge studies in humans are preferable but prohibitively expensive. In the absence of reliable correlates of protection, suitable animal models or a better understanding of the mechanisms underlying protective immunity in humans, vaccine candidate discovery per se may not be sufficient to provide the paradigm shift necessary to develop the next generation of highly effective subunit malaria vaccines.

Geographic-genetic analysis of Plasmodium falciparum parasite populations from surveys of primary school children in Western Kenya

Background. Malaria control, and finally malaria elimination, requires the identification and targeting of residual foci or hotspots of transmission. However, the level of parasite mixing within and between geographical locations is likely to impact the effectiveness and durability of control interventions and thus should be taken into consideration when developing control programs.

Methods. In order to determine the geographic-genetic patterns of Plasmodium falciparum parasite populations at a sub-national level in Kenya, we used the Sequenom platform to genotype 111 genome-wide distributed single nucleotide polymorphic (SNP) positions in 2486 isolates collected from children in 95 primary schools in western Kenya. We analysed these parasite genotypes for genetic structure using principal component analysis and assessed local and global clustering using statistical measures of spatial autocorrelation. We further examined the region for spatial barriers to parasite movement as well as directionality in the patterns of parasite movement.

Results. We found no evidence of population structure and little evidence of spatial autocorrelation of parasite genotypes (correlation coefficients <0.03 among parasite pairs in distance classes of 1km, 2km and 5km; p value<0.01). An analysis of the geographical distribution of allele frequencies showed weak evidence of variation in distribution of alleles, with clusters representing a higher than expected number of samples with the major allele being identified for 5 SNPs. Furthermore, we found no evidence of the existence of spatial barriers to parasite movement within the region, but observed directional movement of parasites among schools in two separate sections of the region studied.

Conclusions. Our findings illustrate a pattern of high parasite mixing within the study region. If this mixing is due to rapid gene flow, then “one-off” targeted interventions may not be currently effective at the sub-national scale in Western Kenya, due to the high parasite movement that is likely to lead to re-introduction of infection from surrounding regions. However repeated targeted interventions may reduce transmission in the surrounding regions.

Geographic-genetic analysis of Plasmodium falciparum parasite populations from surveys of primary school children in Western Kenya

Background. Malaria control, and finally malaria elimination, requires the identification and targeting of residual foci or hotspots of transmission. However, the level of parasite mixing within and between geographical locations is likely to impact the effectiveness and durability of control interventions and thus should be taken into consideration when developing control programs. Methods. In order to determine the geographic-genetic patterns of Plasmodium falciparum parasite populations at a sub-national level in Kenya, we used the Sequenom platform to genotype 111 genome-wide distributed single nucleotide polymorphic (SNP) positions in 2486 isolates collected from children in 95 primary schools in western Kenya. We analysed these parasite genotypes for genetic structure using principal component analysis and assessed local and global clustering using statistical measures of spatial autocorrelation. We further examined the region for spatial barriers to parasite movement as well as directionality in the patterns of parasite movement. Results. We found no evidence of population structure and little evidence of spatial autocorrelation of parasite genotypes (correlation coefficients <0.03 among parasite pairs in distance classes of 1km, 2km and 5km; p value<0.01). An analysis of the geographical distribution of allele frequencies showed weak evidence of variation in distribution of alleles, with clusters representing a higher than expected number of samples with the major allele being identified for 5 SNPs. Furthermore, we found no evidence of the existence of spatial barriers to parasite movement within the region, but observed directional movement of parasites among schools in two separate sections of the region studied. Conclusions. Our findings illustrate a pattern of high parasite mixing within the study region. If this mixing is due to rapid gene flow, then “one-off” targeted interventions may not be currently effective at the sub-national scale in Western Kenya, due to the high parasite movement that is likely to lead to re-introduction of infection from surrounding regions. However repeated targeted interventions may reduce transmission in the surrounding regions.

Selections, frameshift mutations, and copy number variation detected on the surf 4.1 gene in the western Kenyan Plasmodium falciparum population

BACKGROUND

Plasmodium falciparum SURFIN4.1 is a putative ligand expressed on the merozoite and likely on the infected red blood cell, whose gene was suggested to be under directional selection in the eastern Kenyan population, but under balancing selection in the Thai population. To understand this difference, surf 4.1 sequences of western Kenyan P. falciparum isolates were analysed. Frameshift mutations and copy number variation (CNV) were also examined for the parasites from western Kenya and Thailand.

RESULTS

Positively significant departures from neutral expectations were detected on the surf 4.1 region encoding C-terminus of the variable region 2 (Var2) by 3 population-based tests in the western Kenyan population as similar in the Thai population, which was not covered by the previous analysis for eastern Kenyan population. Significant excess of non-synonymous substitutions per nonsynonymous site over synonymous substitutions per synonymous site was also detected in the Var2 region. Negatively significant departures from neutral expectations was detected on the region encoding Var1 C-terminus consistent to the previous observation in the eastern Kenyan population. Parasites possessing a frameshift mutation resulting a product without intracellular Trp-rich (WR) domains were 22/23 in western Kenya and 22/36 in Thailand. More than one copy of surf 4.1 gene was detected in western Kenya (4/24), but no CNV was found in Thailand (0/36).

CONCLUSIONS

The authors infer that the high polymorphism of SURFIN4.1 Var2 C-terminus in both Kenyan and Thai populations were shaped-up by diversifying selection and maintained by balancing selection. These phenomena were most likely driven by immunological pressure. Whereas the SURFIN4.1 Var1 C-terminus is suggested to be under directional selection consistent to the previous report for the eastern Kenyan population. Most western Kenyan isolates possess a frameshift mutation that would limit the expression of SURFIN4.1 on the merozoite, but only 60% of Thai isolates possess this frameshift, which would affect the level and type of the selection pressure against this protein as seen in the two extremities of Tajima's D values for Var1 C-terminus between Kenyan and Thai populations. CNV observed in Kenyan isolates may be a consequence of this frameshift mutation to increase benefits on the merozoite surface.

Out of Africa: origins and evolution of the human malaria parasites Plasmodium falciparum and Plasmodium vivax

Plasmodium falciparum and Plasmodium vivax account for more than 95% of all human malaria infections, and thus pose a serious public health challenge. To control and potentially eliminate these pathogens, it is important to understand their origins and evolutionary history. Until recently, it was widely believed that P. falciparum had co-evolved with humans (and our ancestors) over millions of years, whilst P. vivax was assumed to have emerged in southeastern Asia following the cross-species transmission of a parasite from a macaque. However, the discovery of a multitude of Plasmodium spp. in chimpanzees and gorillas has refuted these theories and instead revealed that both P. falciparum and P. vivax evolved from parasites infecting wild-living African apes. It is now clear that P. falciparum resulted from a recent cross-species transmission of a parasite from a gorilla, whilst P. vivax emerged from an ancestral stock of parasites that infected chimpanzees, gorillas and humans in Africa, until the spread of the protective Duffy-negative mutation eliminated P. vivax from human populations there. Although many questions remain concerning the biology and zoonotic potential of the P. falciparum- and P. vivax-like parasites infecting apes, comparative genomics, coupled with functional parasite and vector studies, are likely to yield new insights into ape Plasmodium transmission and pathogenesis that are relevant to the treatment and prevention of human malaria.

Binding of Plasmodium falciparum Merozoite Surface Proteins DBLMSP and DBLMSP2 to Human Immunoglobulin M Is Conserved among Broadly Diverged Sequence Variants

Diversity at pathogen genetic loci can be driven by host adaptive immune selection pressure and may reveal proteins important for parasite biology. Population-based genome sequencing of Plasmodium falciparum, the parasite responsible for the most severe form of malaria, has highlighted two related polymorphic genes called dblmsp and dblmsp2, which encode Duffy binding-like (DBL) domain-containing proteins located on the merozoite surface but whose function remains unknown. Using recombinant proteins and transgenic parasites, we show that DBLMSP and DBLMSP2 directly and avidly bind human IgM via their DBL domains. We used whole genome sequence data from over 400 African and Asian P. falciparum isolates to show that dblmsp and dblmsp2 exhibit extreme protein polymorphism in their DBL domain, with multiple variants of two major allelic classes present in every population tested. Despite this variability, the IgM binding function was retained across diverse sequence representatives. Although this interaction did not seem to have an effect on the ability of the parasite to invade red blood cells, binding of DBLMSP and DBLMSP2 to IgM inhibited the overall immunoreactivity of these proteins to IgG from patients who had been exposed to the parasite. This suggests that IgM binding might mask these proteins from the host humoral immune system.

Comparison of allele frequencies of Plasmodium falciparum merozoite antigens in malaria infections sampled in different years in a Kenyan population

Background

Plasmodium falciparum merozoite antigens elicit antibody responses in malaria-endemic populations, some of which are clinically protective, which is one of the reasons why merozoite antigens are the focus of malaria vaccine development efforts. Polymorphisms in several merozoite antigen-encoding genes are thought to arise as a result of selection by the human immune system.

Methods

The allele frequency distribution of 15 merozoite antigens over a two-year period, 2007 and 2008, was examined in parasites obtained from children with uncomplicated malaria. In the same population, allele frequency changes pre- and post-anti-malarial treatment were also examined. Any gene which showed a significant shift in allele frequencies was also assessed longitudinally in asymptomatic and complicated malaria infections.

Results

Fluctuating allele frequencies were identified in codons 147 and 148 of reticulocyte-binding homologue (Rh) 5, with a shift from HD to YH haplotypes over the two-year period in uncomplicated malaria infections. However, in both the asymptomatic and complicated malaria infections YH was the dominant and stable haplotype over the two-year and ten-year periods, respectively. A logistic regression analysis of all three malaria infection populations between 2007 and 2009 revealed, that the chance of being infected with the HD haplotype decreased with time from 2007 to 2009 and increased in the uncomplicated and asymptomatic infections.

Conclusion

Rh5 codons 147 and 148 showed heterogeneity at both an individual and population level and may be under some degree of immune selection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-016-1304-8) contains supplementary material, which is available to authorized users.

Multiple Plasmodium falciparum Merozoite Surface Protein 1 Complexes Mediate Merozoite Binding to Human Erythrocytes

Successful invasion of human erythrocytes byPlasmodium falciparummerozoites is required for infection of the host and parasite survival. The early stages of invasion are mediated via merozoite surface proteins that interact with human erythrocytes. The nature of these interactions are currently not well understood, but it is known that merozoite surface protein 1 (MSP1) is critical for successful erythrocyte invasion. Here we show that the peripheral merozoite surface proteins MSP3, MSP6, MSPDBL1, MSPDBL2, and MSP7 bind directly to MSP1, but independently of each other, to form multiple forms of the MSP1 complex on the parasite surface. These complexes have overlapping functions that interact directly with human erythrocytes. We also show that targeting the p83 fragment of MSP1 using inhibitory antibodies inhibits all forms of MSP1 complexes and disrupts parasite growthin vitro.

Whole Genome Sequencing of Field Isolates Reveals Extensive Genetic Diversity in Plasmodium vivax from Colombia

Plasmodium vivax is the most prevalent malarial species in South America and exerts a substantial burden on the populations it affects. The control and eventual elimination of P. vivax are global health priorities. Genomic research contributes to this objective by improving our understanding of the biology of P. vivax and through the development of new genetic markers that can be used to monitor efforts to reduce malaria transmission. Here we analyze whole-genome data from eight field samples from a region in Cordóba, Colombia where malaria is endemic. We find considerable genetic diversity within this population, a result that contrasts with earlier studies suggesting that P. vivax had limited diversity in the Americas. We also identify a selective sweep around a substitution known to confer resistance to sulphadoxine-pyrimethamine (SP). This is the first observation of a selective sweep for SP resistance in this species. These results indicate that P. vivax has been exposed to SP pressure even when the drug is not in use as a first line treatment for patients afflicted by this parasite. We identify multiple non-synonymous substitutions in three other genes known to be involved with drug resistance in Plasmodium species. Finally, we found extensive microsatellite polymorphisms. Using this information we developed 18 polymorphic and easy to score microsatellite loci that can be used in epidemiological investigations in South America.

Insights into the Immunological Properties of Intrinsically Disordered Malaria Proteins Using Proteome Scale Predictions

Malaria remains a significant global health burden. The development of an effective malaria vaccine remains as a major challenge with the potential to significantly reduce morbidity and mortality. While Plasmodium spp. have been shown to contain a large number of intrinsically disordered proteins (IDPs) or disordered protein regions, the relationship of protein structure to subcellular localisation and adaptive immune responses remains unclear. In this study, we employed several computational prediction algorithms to identify IDPs at the proteome level of six Plasmodium spp. and to investigate the potential impact of protein disorder on adaptive immunity against P. falciparum parasites. IDPs were shown to be particularly enriched within nuclear proteins, apical proteins, exported proteins and proteins localised to the parasitophorous vacuole. Furthermore, several leading vaccine candidates, and proteins with known roles in host-cell invasion, have extensive regions of disorder. Presentation of peptides by MHC molecules plays an important role in adaptive immune responses, and we show that IDP regions are predicted to contain relatively few MHC class I and II binding peptides owing to inherent differences in amino acid composition compared to structured domains. In contrast, linear B-cell epitopes were predicted to be enriched in IDPs. Tandem repeat regions and non-synonymous single nucleotide polymorphisms were found to be strongly associated with regions of disorder. In summary, immune responses against IDPs appear to have characteristics distinct from those against structured protein domains, with increased antibody recognition of linear epitopes but some constraints for MHC presentation and issues of polymorphisms. These findings have major implications for vaccine design, and understanding immunity to malaria.

Fifteen years of genomewide scans for selection: trends, lessons and unaddressed genetic sources of complication

Genomewide scans for natural selection (GWSS) have become increasingly common over the last 15 years due to increased availability of genome-scale genetic data. Here, we report a representative survey of GWSS from 1999 to present and find that (i) between 1999 and 2009, 35 of 49 (71%) GWSS focused on human, while from 2010 to present, only 38 of 83 (46%) of GWSS focused on human, indicating increased focus on nonmodel organisms; (ii) the large majority of GWSS incorporate interpopulation or interspecific comparisons using, for example F(ST), cross-population extended haplotype homozygosity or the ratio of nonsynonymous to synonymous substitutions; (iii) most GWSS focus on detection of directional selection rather than other modes such as balancing selection; and (iv) in human GWSS, there is a clear shift after 2004 from microsatellite markers to dense SNP data. A survey of GWSS meant to identify loci positively selected in response to severe hypoxic conditions support an approach to GWSS in which a list of a priori candidate genes based on potential selective pressures are used to filter the list of significant hits a posteriori. We also discuss four frequently ignored determinants of genomic heterogeneity that complicate GWSS: mutation, recombination, selection and the genetic architecture of adaptive traits. We recommend that GWSS methodology should better incorporate aspects of genomewide heterogeneity using empirical estimates of relevant parameters and/or realistic, whole-chromosome simulations to improve interpretation of GWSS results. Finally, we argue that knowledge of potential selective agents improves interpretation of GWSS results and that new methods focused on correlations between environmental variables and genetic variation can help automate this approach.

Identification of large variation in pfcrt, pfmdr-1 and pfubp-1 markers in Plasmodium falciparum isolates from Ethiopia and Tanzania

Background

Plasmodium falciparum resistance to anti-malarials is a major drawback in effective malaria control and elimination globally. Artemisinin-combination therapy (ACT) is currently the key first-line treatment for uncomplicated falciparum malaria. Plasmodium falciparum genetic signatures at pfmdr-1, pfcrt, and pfubp-1 loci are known to modulate in vivo and in vitro parasite response to ACT. The objective of this study was to assess the distribution of these resistance gene markers in isolates collected from different malaria transmission intensity in Ethiopia and Tanzania.

Methods

Plasmodium falciparum clinical isolates were collected from different regions of Ethiopia and Tanzania. Genetic polymorphisms in the genes pfcrt, pfmdr-1 and pfubp-1 were analysed by PCR and sequencing. Frequencies of the different alleles in the three genes were compared within and between regions, and between the two countries.

Results

The majority of the isolates from Ethiopia were mutant for the pfcrt 76 and wild-type for pfmdr-1 86. In contrast, the majority of the Tanzanian samples were wild-type for both pfcrt and pfmdr-1 loci. Analysis of a variable linker region in pfmdr-1 showed substantial variation in isolates from Tanzania as compared to Ethiopian isolates that had minimal variation. Direct sequencing of the pfubp-1 region showed that 92.8% (26/28) of the Ethiopian isolates had identical genome sequence with the wild type reference P. falciparum strain 3D7. Of 42 isolates from Tanzania, only 13 (30.9%) had identical genome sequences with 3D7. In the Tanzanian samples, 10 variant haplotypes were identified.

Conclusion

The majority of Ethiopian isolates carried the main marker for chloroquine (CQ) resistance, while the majority of the samples from Tanzania carried markers for CQ susceptibility. Polymorphic genes showed substantially more variation in Tanzanian isolates. The low variability in the polymorphic region of pfmdr-1 in Ethiopia may be a consequence of low transmission intensity as compared to high transmission intensity and large variations in Tanzania.

Inferring natural selection signals in Plasmodium vivax-encoded proteins having a potential role in merozoite invasion

Detecting natural selection signals in Plasmodium parasites antigens might be used for identifying potential new vaccine candidates. Fifty-nine Plasmodium vivax-Sal-I genes encoding proteins having a potential role in invasion were used as query for identifying them in recent P. vivax strain genome sequences and two closely-related Plasmodium species. Several measures of DNA sequence variation were then calculated and selection signatures were detected by using different approaches. Our results may be used for determining which genes expressed during P. vivax merozoite stage could be prioritised for further population genetics or functional studies for designing a P. vivax vaccine which would avoid allele-specific immune responses.

Antibodies to the Plasmodium falciparum Proteins MSPDBL1 and MSPDBL2 Opsonize Merozoites, Inhibit Parasite Growth, and Predict Protection From Clinical Malaria

Increasing evidence suggests that antibodies against merozoite surface proteins (MSPs) play an important role in clinical immunity to malaria. Two unusual members of the MSP-3 family, merozoite surface protein duffy binding-like (MSPDBL)1 and MSPDBL2, have been shown to be extrinsically associated to MSP-1 on the parasite surface. In addition to a secreted polymorphic antigen associated with merozoite (SPAM) domain characteristic of MSP-3 family members, they also contain Duffy binding-like (DBL) domain and were found to bind to erythrocytes, suggesting that they play a role in parasite invasion. Antibody responses to these proteins were investigated in a treatment-reinfection study conducted in an endemic area of Papua New Guinea to determine their contribution to naturally acquired immunity. Antibodies to the SPAM domains of MSPDBL1 and MSPDBL2 as well as the DBL domain of MSPDBL1 were found to be associated with protection from Plasmodium falciparum clinical episodes. Moreover, affinity-purified anti-MSPDBL1 and MSPDBL2 were found to inhibit in vitro parasite growth and had strong merozoite opsonizing capacity, suggesting that protection targeting these antigens results from ≥2 distinct effector mechanisms. Together these results indicate that MSPDBL1 and MSPDBL2 are important targets of naturally acquired immunity and might constitute potential vaccine candidates.

The MSPDBL2 codon 591 polymorphism is associated with lumefantrine in vitro drug responses in Plasmodium falciparum isolates from Kilifi, Kenya

The mechanisms of drug resistance development in the Plasmodium falciparum parasite to lumefantrine (LUM), commonly used in combination with artemisinin, are still unclear. We assessed the polymorphisms of Pfmspdbl2 for associations with LUM activity in a Kenyan population. MSPDBL2 codon 591S was associated with reduced susceptibility to LUM (P = 0.04). The high frequency of Pfmspdbl2 codon 591S in Kenya may be driven by the widespread use of lumefantrine in artemisinin combination therapy (Coartem).

The merozoite surface protein 1 complex is a platform for binding to human erythrocytes by Plasmodium falciparum

Plasmodium falciparum is the causative agent of the most severe form of malaria in humans. The merozoite, an extracellular stage of the parasite lifecycle, invades erythrocytes in which they develop. The most abundant protein on the surface of merozoites is merozoite surface protein 1 (MSP1), which consists of four processed fragments. Studies indicate that MSP1 interacts with other peripheral merozoite surface proteins to form a large complex. Successful invasion of merozoites into host erythrocytes is dependent on this protein complex; however, the identity of all components and its function remain largely unknown. We have shown that the peripheral merozoite surface proteins MSPDBL1 and MSPDBL2 are part of the large MSP1 complex. Using surface plasmon resonance, we determined the binding affinities of MSPDBL1 and MSPDBL2 to MSP1 to be in the range of 2-4 × 10(-7) m. Both proteins bound to three of the four proteolytically cleaved fragments of MSP1 (p42, p38, and p83). In addition, MSPDBL1 and MSPDBL2, but not MSP1, bound directly to human erythrocytes. This demonstrates that the MSP1 complex acts as a platform for display of MSPDBL1 and MSPDBL2 on the merozoite surface for binding to receptors on the erythrocyte and invasion.

Strategies for designing and monitoring malaria vaccines targeting diverse antigens

After more than 50 years of intensive research and development, only one malaria vaccine candidate, “RTS,S,” has progressed to Phase 3 clinical trials. Despite only partial efficacy, this candidate is now forecast to become the first licensed malaria vaccine. Hence, more efficacious second-generation malaria vaccines that can significantly reduce transmission are urgently needed. This review will focus on a major obstacle hindering development of effective malaria vaccines: parasite antigenic diversity. Despite extensive genetic diversity in leading candidate antigens, vaccines have been and continue to be formulated using recombinant antigens representing only one or two strains. These vaccine strains represent only a small fraction of the diversity circulating in natural parasite populations, leading to escape of non-vaccine strains and challenging investigators’ abilities to measure strain-specific efficacy in vaccine trials. Novel strategies are needed to overcome antigenic diversity in order for vaccine development to succeed. Many studies have now cataloged the global diversity of leading Plasmodium falciparum and Plasmodium vivax vaccine antigens. In this review, we describe how population genetic approaches can be applied to this rich data source to predict the alleles that best represent antigenic diversity, polymorphisms that contribute to it, and to identify key polymorphisms associated with antigenic escape. We also suggest an approach to summarize the known global diversity of a given antigen to predict antigenic diversity, how to select variants that best represent the strains circulating in natural parasite populations and how to investigate the strain-specific efficacy of vaccine trials. Use of these strategies in the design and monitoring of vaccine trials will not only shed light on the contribution of genetic diversity to the antigenic diversity of malaria, but will also maximize the potential of future malaria vaccine candidates.

Variation in the circumsporozoite protein of Plasmodium falciparum: vaccine development implications

The malaria vaccine candidate RTS,S/AS01 is based on immunogenic regions of Plasmodium falciparum circumsporozoite protein (CSP) from the 3D7 reference strain and has shown modest efficacy against clinical disease in African children. It remains unclear what aspect(s) of the immune response elicited by this vaccine are protective. The goals of this study were to measure diversity in immunogenic regions of CSP, and to identify associations between polymorphism in CSP and the risk of P. falciparum infection and clinical disease. The present study includes data and samples from a prospective cohort study designed to measure incidence of malaria infection and disease in children in Bandiagara, Mali. A total of 769 parasite-positive blood samples corresponding to both acute clinical malaria episodes and asymptomatic infections experienced by 100 children were included in the study. Non-synonymous SNP data were generated by 454 sequencing for the T-cell epitopes, and repeat length data were generated for the B-cell epitopes of the cs gene. Cox proportional hazards models were used to determine the effect of sequence variation in consecutive infections occurring within individuals on the time to new infection and new clinical malaria episode. Diversity in the T-cell epitope-encoding regions Th2R and Th3R remained stable throughout seasons, between age groups and between clinical and asymptomatic infections with the exception of a higher proportion of 3D7 haplotypes found in the oldest age group. No associations between sequence variation and hazard of infection or clinical malaria were detected. The lack of association between sequence variation and hazard of infection or clinical malaria suggests that naturally acquired immunity to CSP may not be allele-specific.

Whole-genome scans provide evidence of adaptive evolution in Malawian Plasmodium falciparum isolates

Background

Selection by host immunity and antimalarial drugs has driven extensive adaptive evolution in Plasmodium falciparum and continues to produce ever-changing landscapes of genetic variation.

Methods

We performed whole-genome sequencing of 69 P. falciparum isolates from Malawi and used population genetics approaches to investigate genetic diversity and population structure and identify loci under selection.

Results

High genetic diversity (π = 2.4 × 10⁻⁴), moderately high multiplicity of infection (2.7), and low linkage disequilibrium (500-bp) were observed in Chikhwawa District, Malawi, an area of high malaria transmission. Allele frequency–based tests provided evidence of recent population growth in Malawi and detected potential targets of host immunity and candidate vaccine antigens. Comparison of the sequence variation between isolates from Malawi and those from 5 geographically dispersed countries (Kenya, Burkina Faso, Mali, Cambodia, and Thailand) detected population genetic differences between Africa and Asia, within Southeast Asia, and within Africa. Haplotype-based tests of selection to sequence data from all 6 populations identified signals of directional selection at known drug-resistance loci, including pfcrt, pfdhps, pfmdr1, and pfgch1.

Conclusions

The sequence variations observed at drug-resistance loci reflect differences in each country's historical use of antimalarial drugs and may be useful in formulating local malaria treatment guidelines.

Genome-wide analysis of selection on the malaria parasite Plasmodium falciparum in West African populations of differing infection endemicity

Locally varying selection on pathogens may be due to differences in drug pressure, host immunity, transmission opportunities between hosts, or the intensity of between-genotype competition within hosts. Highly recombining populations of the human malaria parasite Plasmodium falciparum throughout West Africa are closely related, as gene flow is relatively unrestricted in this endemic region, but markedly varying ecology and transmission intensity should cause distinct local selective pressures. Genome-wide analysis of sequence variation was undertaken on a sample of 100 P. falciparum clinical isolates from a highly endemic region of the Republic of Guinea where transmission occurs for most of each year and compared with data from 52 clinical isolates from a previously sampled population from The Gambia, where there is relatively limited seasonal malaria transmission. Paired-end short-read sequences were mapped against the 3D7 P. falciparum reference genome sequence, and data on 136,144 single nucleotide polymorphisms (SNPs) were obtained. Within-population analyses identifying loci showing evidence of recent positive directional selection and balancing selection confirm that antimalarial drugs and host immunity have been major selective agents. Many of the signatures of recent directional selection reflected by standardized integrated haplotype scores were population specific, including differences at drug resistance loci due to historically different antimalarial use between the countries. In contrast, both populations showed a similar set of loci likely to be under balancing selection as indicated by very high Tajima’s D values, including a significant overrepresentation of genes expressed at the merozoite stage that invades erythrocytes and several previously validated targets of acquired immunity. Between-population F_ST analysis identified exceptional differentiation of allele frequencies at a small number of loci, most markedly for five SNPs covering a 15-kb region within and flanking the gdv1 gene that regulates the early stages of gametocyte development, which is likely related to the extreme differences in mosquito vector abundance and seasonality that determine the transmission opportunities for the sexual stage of the parasite.

Multilocus patterns of genetic variation across Cryptosporidium species suggest balancing selection at the gp60 locus

Cryptosporidium is an apicomplexan protozoan that lives in most vertebrates, including humans. Its gp60 gene is functionally involved in its attachment to host cells, and its high level of genetic variation has made it the reference marker for sample typing in epidemiological studies. To understand the origin of such high diversity and to determine the extent to which this classification applies to the rest of the genome, we analysed the patterns of variation at gp60 and nine other nuclear loci in isolates of three Cryptosporidium species. Most loci showed low genetic polymorphism (πS <1%) and similar levels of between-species divergence. Contrastingly, gp60 exhibited very different characteristics: (i) it was nearly ten times more variable than the other loci; (ii) it displayed a significant excess of polymorphisms relative to between-species differences in a maximum-likelihood Hudson-Kreitman-Aguadé test; (iii) gp60 subtypes turned out to be much older than the species they were found in; and (iv) showed a significant excess of polymorphic variants shared across species from random expectations. These observations suggest that this locus evolves under balancing selection and specifically under negative frequency-dependent selection (FDS). Interestingly, genetic variation at the other loci clusters very well within the groups of isolates defined by gp60 subtypes, which may provide new tools to understand the genome-wide patterns of genetic variation of the parasite in the wild. These results suggest that gp60 plays an active and essential role in the life cycle of the parasite and that genetic variation at this locus might be essential for the parasite's long-term success.

Analysis of antibodies to newly described Plasmodium falciparum merozoite antigens supports MSPDBL2 as a predicted target of naturally acquired immunity

Prospective studies continue to identify malaria parasite genes with particular patterns of polymorphism which indicate they may be under immune selection, and the encoded proteins require investigation. Sixteen new recombinant protein reagents were designed to characterize three such polymorphic proteins expressed in Plasmodium falciparum schizonts and merozoites: MSPDBL1 (also termed MSP3.4) and MSPDBL2 (MSP3.8), which possess Duffy binding-like (DBL) domains, and SURFIN4.2, encoded by a member of the surface-associated interspersed (surf) multigene family. After testing the antigenicities of these reagents by murine immunization and parasite immunofluorescence, we analyzed naturally acquired antibody responses to the antigens in two cohorts in coastal Kenya in which the parasite was endemic (Chonyi [n = 497] and Ngerenya [n = 461]). As expected, the prevalence and levels of serum antibodies increased with age. We then investigated correlations with subsequent risk of clinical malaria among children <11 years of age during 6 months follow-up surveillance. Antibodies to the polymorphic central region of MSPDBL2 were associated with reduced risk of malaria in both cohorts, with statistical significance remaining for the 3D7 allelic type after adjustment for individuals' ages in years and antibody reactivity to whole-schizont extract (Chonyi, risk ratio, 0.51, and 95% confidence interval [CI], 0.28 to 0.93; Ngerenya, risk ratio, 0.38, and 95% CI, 0.18 to 0.82). For the MSPDBL1 Palo Alto allelic-type antigen, there was a protective association in one cohort (Ngerenya, risk ratio, 0.53, and 95% CI, 0.32 to 0.89), whereas the other antigens showed no protective associations after adjustment. These findings support the prediction that antibodies to the polymorphic region of MSPDBL2 contribute to protective immunity.

Modulation of PF10_0355 (MSPDBL2) alters Plasmodium falciparum response to antimalarial drugs

Malaria's ability to rapidly adapt to new drugs has allowed it to remain one of the most devastating infectious diseases of humans. Understanding and tracking the genetic basis of these adaptations are critical to the success of treatment and intervention strategies. The novel antimalarial resistance locus PF10_0355 (Pfmspdbl2) was previously associated with the parasite response to halofantrine, and functional validation confirmed that overexpression of this gene lowered parasite sensitivity to both halofantrine and the structurally related antimalarials mefloquine and lumefantrine, predominantly through copy number variation. Here we further characterize the role of Pfmspdbl2 in mediating the antimalarial drug response of Plasmodium falciparum. Knockout of Pfmspdbl2 increased parasite sensitivity to halofantrine, mefloquine, and lumefantrine but not to unrelated antimalarials, further suggesting that this gene mediates the parasite response to a specific class of antimalarial drugs. A single nucleotide polymorphism encoding a C591S mutation within Pfmspdbl2 had the strongest association with halofantrine sensitivity and showed a high derived allele frequency among Senegalese parasites. Transgenic parasites expressing the ancestral Pfmspdbl2 allele were more sensitive to halofantrine and structurally related antimalarials than were parasites expressing the derived allele, revealing an allele-specific effect on drug sensitivity in the absence of copy number effects. Finally, growth competition experiments showed that under drug pressure, parasites expressing the derived allele of Pfmspdbl2 outcompeted parasites expressing the ancestral allele within a few generations. Together, these experiments demonstrate that modulation of Pfmspdbl2 affects malaria parasite responses to antimalarial drugs.

Population genomic scan for candidate signatures of balancing selection to guide antigen characterization in malaria parasites

Acquired immunity in vertebrates maintains polymorphisms in endemic pathogens, leading to identifiable signatures of balancing selection. To comprehensively survey for genes under such selection in the human malaria parasite Plasmodium falciparum, we generated paired-end short-read sequences of parasites in clinical isolates from an endemic Gambian population, which were mapped to the 3D7 strain reference genome to yield high-quality genome-wide coding sequence data for 65 isolates. A minority of genes did not map reliably, including the hypervariable var, rifin, and stevor families, but 5,056 genes (90.9% of all in the genome) had >70% sequence coverage with minimum read depth of 5 for at least 50 isolates, of which 2,853 genes contained 3 or more single nucleotide polymorphisms (SNPs) for analysis of polymorphic site frequency spectra. Against an overall background of negatively skewed frequencies, as expected from historical population expansion combined with purifying selection, the outlying minority of genes with signatures indicating exceptionally intermediate frequencies were identified. Comparing genes with different stage-specificity, such signatures were most common in those with peak expression at the merozoite stage that invades erythrocytes. Members of clag, PfMC-2TM, surfin, and msp3-like gene families were highly represented, the strongest signature being in the msp3-like gene PF10_0355. Analysis of msp3-like transcripts in 45 clinical and 11 laboratory adapted isolates grown to merozoite-containing schizont stages revealed surprisingly low expression of PF10_0355. In diverse clonal parasite lines the protein product was expressed in a minority of mature schizonts (<1% in most lines and ∼10% in clone HB3), and eight sub-clones of HB3 cultured separately had an intermediate spectrum of positive frequencies (0.9 to 7.5%), indicating phase variable expression of this polymorphic antigen. This and other identified targets of balancing selection are now prioritized for functional study.

Molecular signatures identify a candidate target of balancing selection in an arcD-like gene of Staphylococcus epidermidis

A comparative population genetics study revealed high levels of nucleotide polymorphism and intermediate-frequency alleles in an arcC gene of Staphylococcus epidermidis, but not in a homologous gene of the more aggressive human pathogen, Staphylococcus aureus. Further investigation showed that the arcC genes used in the multilocus sequence typing schemes of these two species were paralogs. Phylogenetic analyses of arcC-containing loci, including the arginine catabolic mobile element, from both species, suggested that these loci had an eventful history involving gene duplications, rearrangements, deletions, and horizontal transfers. The peak signatures in the polymorphic S. epidermidis locus were traced to an arcD-like gene adjacent to arcC; these signatures consisted of unusually elevated Tajima's D and π/K ratios, which were robust to assumptions about recombination and species divergence time and among the most elevated in the S. epidermidis genome. Amino acid polymorphisms, including one that differed in polarity and hydropathy, were located in the peak signatures and defined two allelic lineages. Recombination events were detected between these allelic lineages and potential donors and recipients of S. epidermidis were identified in each case. By comparison, the orthologous gene of S. aureus showed no unusual signatures. The ArcD-like protein belonged to the unknown ion transporter 3 family and appeared to be unrelated to ArcD from the arginine deiminase pathway. These studies report the first comparative population genetics results for staphylococci and the first statistical evidence for a candidate target of balancing selection in S. epidermidis.

VarB: a variation browsing and analysis tool for variants derived from next-generation sequencing data

Summary: There is an immediate need for tools to both analyse and visualize in real-time single-nucleotide polymorphisms, insertions and deletions, and other structural variants from new sequence file formats. We have developed VarB software that can be used to visualize variant call format files in real time, as well as identify regions under balancing selection and informative markers to differentiate user-defined groups (e.g. populations). We demonstrate its utility using sequence data from 50 Plasmodium falciparum isolates comprising two different continents and confirm known signals from genomic regions that contain important antigenic and anti-malarial drug-resistance genes.

Availability and implementation: The C++-based software VarB and user manual are available from www.pathogenseq.org/varb.

Contact:taane.clark@lshtm.ac.uk

Positive diversifying selection on the Plasmodium falciparum surf4.1 gene in Thailand

Plasmodium falciparum SURFIN_4.1 is a type I transmembrane protein thought to locate on the merozoite surface and to be responsible for a reversible adherence to the erythrocyte before invasion. In this study, we evaluated surf_4.1 gene segment encoding extracellular region for polymorphism, the signature of positive selection, the degree of linkage disequilibrium, and temporal change in allele frequency distribution in P. falciparum isolates from Thailand in 1988–89, 2003, and 2005. We found that SURFIN_4.1 is highly polymorphic, particularly at the C-terminal side of the variable region located just before a predicted transmembrane region. A signature of positive diversifying selection on the variable region was detected by multiple tests and, to a lesser extent, on conserved N-terminally located cysteine-rich domain by Tajima’s D test. Linkage disequilibrium between sites over a long distance (> 1.5 kb) was detected, and multiple SURFIN_4.1 haplotype sequences detected in 1988/89 still circulated in 2003. Few of the single amino acid polymorphism allele frequency distributions were significantly different between the 1988/89 and 2003 groups, suggesting that the frequency distribution of SURFIN_4.1 extracellular region remained stable over 14 years.

Insights into Duffy binding-like domains through the crystal structure and function of the merozoite surface protein MSPDBL2 from Plasmodium falciparum

Invasion of human red blood cells by Plasmodium falciparum involves interaction of the merozoite form through proteins on the surface coat. The erythrocyte binding-like protein family functions after initial merozoite interaction by binding via the Duffy binding-like (DBL) domain to receptors on the host red blood cell. The merozoite surface proteins DBL1 and -2 (PfMSPDBL1 and PfMSPDBL2) (PF10_0348 and PF10_0355) are extrinsically associated with the merozoite, and both have a DBL domain in each protein. We expressed and refolded recombinant DBL domains for PfMSPDBL1 and -2 and show they are functional. The red cell binding characteristics of these domains were shown to be similar to full-length forms of these proteins isolated from parasite cultures. Futhermore, metal cofactors were found to enhance the binding of both the DBL domains and the parasite-derived full-length proteins to erythrocytes, which has implications for receptor binding of other DBL-containing proteins in Plasmodium spp. We solved the structure of the erythrocyte-binding DBL domain of PfMSPDBL2 to 2.09 Å resolution and modeled that of PfMSPDBL1, revealing a canonical DBL fold consisting of a boomerang shaped α-helical core formed from three subdomains. PfMSPDBL2 is highly polymorphic, and mapping of these mutations shows they are on the surface, predominantly in the first two domains. For both PfMSPDBL proteins, polymorphic variation spares the cleft separating domains 1 and 2 from domain 3, and the groove between the two major helices of domain 3 extends beyond the cleft, indicating these regions are functionally important and are likely to be associated with the binding of a receptor on the red blood cell.

Next generation sequencing to detect variation in the Plasmodium falciparum circumsporozoite protein

The malaria vaccine RTS,S/AS01, based on immunogenic regions of the Plasmodium falciparum circumsporozoite protein (CSP), has partial efficacy against clinical malaria in African children. Understanding how sequence diversity in CSP T- and B-cell epitopes relates to naturally acquired and vaccine-induced immunity may be useful in efforts to improve the efficacy of CSP-based vaccines. However, limitations in sequencing technology have precluded thorough evaluation of diversity in the immunogenic regions of this protein. In this study, 454, a next generation sequencing technology, was evaluated as a method for assessing diversity in these regions. Portions of the circumsporozoite gene (cs) were sequenced both by 454 and Sanger sequencing from samples collected in a study in Bandiagara, Mali. 454 detected more single nucleotide polymorphisms and haplotypes in the T-cell epitopes than Sanger sequencing, and it was better able to resolve genetic diversity in samples with multiple infections; however, it failed to generate sequence for the B-cell epitopes.

Harnessing genomics and genome biology to understand malaria biology

Malaria is an important human disease and is the target of a global eradication campaign. New technological and informatics advancements in population genomics are being leveraged to identify genetic loci under selection in the malaria parasite and to find variants that are associated with key clinical phenotypes, such as drug resistance. This article provides a timely Review of how population-genetics-based strategies are being applied to Plasmodium falciparum both to identify genetic loci as key targets of interventions and to develop monitoring and surveillance tools that are crucial for the successful elimination and eradication of malaria.

No recent adaptive selection on the apyrase of Mediterranean Phlebotomus: implications for using salivary peptides to vaccinate against canine leishmaniasis

Vaccine development is informed by a knowledge of genetic variation among antigen alleles, especially the distribution of positive and balancing selection in populations and species. A combined approach using population genetic and phylogenetic methods to detect selective signatures can therefore be informative for identifying vaccine candidates. Parasitic Leishmania species cause the disease leishmaniasis in humans and mammalian reservoir hosts after inoculation by female phlebotomine sandflies. Like other arthropod vectors of disease agents, sandflies use salivary peptides to counteract host haemostatic and immunomodulatory responses during bloodfeeding, and these peptides are vaccine candidates because they can protect against Leishmania infection. We detected no contemporary adaptive selection on one salivary peptide, apyrase, in 20 populations of Phlebotomus ariasi, a European vector of Leishmania infantum. Maximum likelihood branch models on a gene phylogeny showed apyrase to be a single copy in P. ariasi but an ancient duplication event associated with temporary positive selection was observed in its sister group, which contains most Mediterranean vectors of L. infantum. The absence of contemporary adaptive selection on the apyrase of P. ariasi may result from this sandfly's opportunistic feeding behaviour. Our study illustrates how the molecular population genetics of arthropods can help investigate the potential of salivary peptides for disease control and for understanding geographical variation in vector competence.

Transcriptional variation in the malaria parasite Plasmodium falciparum

Malaria genetic variation has been extensively characterized, but the level of epigenetic plasticity remains largely unexplored. Here we provide a comprehensive characterization of transcriptional variation in the most lethal malaria parasite, Plasmodium falciparum, based on highly accurate transcriptional analysis of isogenic parasite lines grown under homogeneous conditions. This analysis revealed extensive transcriptional heterogeneity within genetically homogeneous clonal parasite populations. We show that clonally variant expression controlled at the epigenetic level is an intrinsic property of specific genes and gene families, the majority of which participate in host–parasite interactions. Intrinsic transcriptional variability is not restricted to genes involved in immune evasion, but also affects genes linked to lipid metabolism, protein folding, erythrocyte remodeling, or transcriptional regulation, among others, indicating that epigenetic variation results in both antigenic and functional variation. We observed a general association between heterochromatin marks and clonally variant expression, extending previous observations for specific genes to essentially all variantly expressed gene families. These results suggest that phenotypic variation of functionally unrelated P. falciparum gene families is mediated by a common mechanism based on reversible formation of H3K9me3-based heterochromatin. In changing environments, diversity confers fitness to a population. Our results support the idea that P. falciparum uses a bet-hedging strategy, as an alternative to directed transcriptional responses, to adapt to common fluctuations in its environment. Consistent with this idea, we found that transcriptionally different isogenic parasite lines markedly differed in their survival to heat-shock mimicking febrile episodes and adapted to periodic heat-shock with a pattern consistent with natural selection of pre-existing parasites.

Stable allele frequency distribution of the polymorphic region of SURFIN(4.2) in Plasmodium falciparum isolates from Thailand

Plasmodium falciparum SURFIN₄.₂ (PFD1160w) is a polymorphic protein expressed on the surface of parasite-infected erythrocytes. Such molecules are expected to be under strong host immune pressure, thus we analyzed the nucleotide diversity of the N-terminal extracellular region of SURFIN₄.₂ using P. falciparum isolates obtained from a malaria hypoendemic area of Thailand. The extracellular region of SURFIN₄.₂ was divided into four regions based on the amino acid sequence conservation among SURFIN members and the level of polymorphism among SURFIN₄.₂ sequences; N-terminal segment (Nter), a cysteine-rich domain (CRD), a variable region 1 (Var1), and a variable region 2 (Var2). Comparison between synonymous and non-synonymous substitutions, Tajima's D test, and Fu and Li's D* and F* tests detected signatures of positive selection on Var2 and to a lesser extent Var1, suggesting that these regions were likely under host immune pressure. Strong linkage disequilibrium was detected for nucleotide pairs separated by a distance of more than 1.5 kb, and 7 alleles among 19 alleles detected in 1988-1989 still circulated 14 years later, suggesting low recombination of the analyzed surf₄.₂ sequence region in Thailand. The allele frequency distribution of polymorphic areas in Var2 did not differ between two groups collected in different time points, suggesting the allele frequency distribution of this region was stable for 14 years. The observed allele frequency distribution of SURFIN₄.₂ Var2 may be fixed in Thai P. falciparum population as similar to the observation for P. falciparum merozoite surface protein 1, for which a stable allele frequency distribution was reported.