Assessing the effect of natural selection in malaria parasites

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

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

Molecular epidemiology of the Old World screwworm fly (OWSF) in Iraq, and the genetic structure of various OWSF populations worldwide

Despite being endemic in Iraq, no reports have been published in the past 10 years to update the molecular epidemiology of the Old World screwworm fly (OWSF), Chrysomya bezziana, in this country. In the present study, 130 sheep from 10 Iraqi governorates were found infected with C. bezziana larvae, whose identities were PCR-confirmed based on the cytochrome b (Cytb) gene, and 23 isolates from various tested governorates were successfully sequenced. Although most isolates (n = 20) belonged to the common haplotype circulating in Iraq, two new haplotypes were detected. Significant changes in OWSF epidemiology in Iraq were also suggested, since infestations were detected, for the first time, in Nineveh governorate. Isolates of the present study were combined to those previously published from Iraq and worldwide, collected after searching the GenBank, and various genetic and population structure analyses were conducted. These isolates displayed a great statistically significant value when tested for the purifying (negative) selection, suggesting the limited occurrence of genetic variations, which was evidenced by the high sequence conservation (C = 0.937) value detected. A few isolates from Africa were revealed during our search, and clustered in a separate lineage other than that of the Asian isolates. The latter displayed different genetic variation patterns when compared. For example, isolates from geographically separate regions, e.g., the Gulf Arab countries and South-Eastern Asia had marked genetic differences. On the other hand, isolates from regions with close geographic proximity (the Gulf Arab countries and Iran) had limited genetic subdivision. This is not the case when comparing isolates from 10 islands in the Indonesian Archipelago. Populations from Sumatra and Sumba were isolated and displayed high genetic variations toward the other populations. On the contrary, populations from Sulawesi, Lombok and Sumbawa displayed limited genetic variations. This is particularly important, since it can help detecting the dynamics of establishing the sterile insect technique over various regions as an effective control strategy against the OWSFs.

The impacts of thermal heterogeneity across microhabitats on post-settlement selection of intertidal mussels

Rapid genetic selection is critical for allowing natural populations to adapt to different thermal environments such as those that occur across intertidal microhabitats with high degrees of thermal heterogeneity. To address the question of how thermal regimes influence selection and adaptation in the intertidal black mussel Mytilisepta virgata, we continuously recorded environmental temperatures in both tidal pools and emergent rock microhabitats and then assessed genetic differentiation, gene expression patterns, RNA editing level, and cardiac performance. Our results showed that the subpopulations in the tidal pool and on emergent rocks had different genetic structures and exhibited different physiological and molecular responses to high-temperature stress. These results indicate that environmental heterogeneity across microhabitats is important for driving genetic differentiation and shed light on the importance of post-settlement selection for adaptively modifying the genetic composition and thermal responses of these intertidal mussels.

Genetic diversity of merozoite surface protein-5 (MSP-5) of Plasmodium vivax isolates from Malaria patients in Iran

Malaria has not yet been eradicated in Iran, and Plasmodium vivax (P. vivax) is the main cause of malaria in the country. This study aimed to investigate and analyze the amount of genetic diversity of Plasmodium vivax merozoite surface protein-5 (PvMSP-5) exon 1 gene in the southeast of Iran.Thirty-five patients with clinical symptoms of P. vivax malaria participated. The exon 1 of PvMSP-5 was amplified by PCR, and the PCR product of all isolates was sequenced, and genetic polymorphisms were determined using various genetic software.The analysis showed that studied isolates are different from one another in the DnaSP software version. Out of the 612 sites, 477 were monomorphic and 135 were segregated. The total number of mutations was 143. The singleton variable and the parsimony informative sites were 23 and 112, respectively. There were 17 specific haplotypes with haplotype diversity equal to 0.943. Nucleotide diversity was equal to 0.06766 in the isolates. The ratio of nonsynonymous (0.06446) to synonymous (0.07909) mutations was 0.815020. Tajima's D, which expressed coding, and non-coding regions, was 0.72403, which was not deemed significant (P > 0.10).The analysis of intrapopulation diversity revealed nucleotide and haplotype diversity in the msp-5 gene of Iranian P. vivax isolates. In addition to balancing or purifying selection, intragenic recombination also contributed to the variation observed in exon 1 of PvMSP-5, according to the findings.

Genetic polymorphism and natural selection of the erythrocyte binding antigen 175 region II in Plasmodium falciparum populations from Myanmar and Vietnam

Plasmodium falciparum erythrocyte binding antigen 175 (PfEBA-175) plays essential role in erythrocyte invasion by the parasite and is a leading vaccine candidate. However, its genetic diversity in global isolates is a concern in developing an universal vaccine incorporating this protein. This study aimed to investigate genetic polymorphisms and natural selection of pfeba-175 region II (RII) in Myanmar and Vietnam P. falciparum isolates. Vietnam pfeba-175 RII displayed a low genetic polymorphism, while Myanmar pfeba-175 RII showed high levels of genetic diversity across the region. Point mutations, deletion, and recombinations were main factors contributing to genetic diversities in P. falciparum populations. Global pfeba-175 RII revealed similar, but not identical, genetic polymorphisms and natural selection profiles. Despite profiles of amino acid substitutions differed among populations, five major amino acid changes (K279E, E403K, K481I, Q584K, and R664) were commonly detected in global pfeba-175 RII populations. Haplotype network and genetic differentiation analyses of global pfeba-175 RII populations demonstrated no geographical relationships. Non-neglectable level of genetic diversity was observed in global pfeba-175 RII populations, emphasizing the need to consider this when designing an effective vaccine based on this protein. This study underscores the importance of the continuous monitoring of genetic diversity of pfeba-175 RII in the global P. falciparum populations.

Genetic Diversity of Merozoite Surface Antigens in Global Babesia bovis Populations

Cattle can be severely infected with the tick-borne protozoa Babesia bovis, giving rise to serious economic losses. Invasion of the host's RBCs by the parasite merozoite/sporozoites depends largely on the MSA (merozoite surface antigens) gene family, which comprises various fragments, e.g., MSA-1, MSA-2a1, MSA-2a2, MSA-2b and MSA-2c, highlighting the importance of these antigens as vaccine candidates. However, experimental trials documented the failure of some developed MSA-based vaccines to fully protect animals from B. bovis infection. One reason for this failure may be related to the genetic structure of the parasite. In the present study, all MSA-sequenced B. bovis isolates on the GenBank were collected and subjected to various analyses to evaluate their genetic diversity and population structure. The analyses were conducted on 199 MSA-1, 24 MSA-2a1, 193 MSA-2b and 148 MSA-2c isolates from geographically diverse regions. All these fragments displayed high nucleotide and haplotype diversities, but the MSA-1 was the most hypervariable and had the lowest inter- and intra-population gene flow values. This fragment also displayed a strong positive selection when testing its isolates for the natural selection, which suggests the potential occurrence of more genetic variations. On the contrary, the MSA-2c was the most conserved in comparison to the other fragments, and displayed the highest inter- and intra-population gene flow values, which was evidenced by a significantly negative selection and negative neutrality indices (Fu's Fs and Tajima's D). The majority of the MSA-2c tested isolates had two conserved amino acid repeats, and earlier reports have found these repeats to be highly immunogenic, which underlines the importance of this fragment in developing vaccines against B. bovis. Results of the MSA-2a1 analyses were also promising, but many more MSA-2a1 sequenced isolates are required to validating this assumption. The genetic analyses conducted for the MSA-2b fragment displayed borderline values when compared to the other fragments.

Genetic diversity and natural selection of Plasmodium vivax reticulocyte invasion genes in Ecuador

BACKGROUND

Knowledge of the diversity of invasion ligands in malaria parasites in endemic regions is essential to understand how natural selection influences genetic diversity of these ligands and their feasibility as possible targets for future vaccine development. In this study the diversity of four genes for merozoite invasion ligands was studied in Ecuadorian isolates of Plasmodium vivax.

METHODS

Eighty-eight samples from P. vivax infected individuals from the Coast and Amazon region of Ecuador were obtained between 2012 and 2015. The merozoite invasion genes pvmsp-1-19, pvdbpII, pvrbp1a-2 and pvama1 were amplified, sequenced, and compared to the Sal-1 strain. Polymorphisms were mapped and genetic relationships between haplotypes were determined.

RESULTS

Only one nonsynonymous polymorphism was detected in pvmsp-1-19, while 44 nonsynonymous polymorphisms were detected in pvdbpII, 56 in pvrbp1a-2 and 33 in pvama1. While haplotypes appeared to be more related within each area of study and there was less relationship between parasites of the coastal and Amazon regions of the country, diversification processes were observed in the two Amazon regions. The highest haplotypic diversity for most genes occurred in the East Amazon of the country. The high diversity observed in Ecuadorian samples is closer to Brazilian and Venezuelan isolates, but lower than reported in other endemic regions. In addition, departure from neutrality was observed in Ecuadorian pvama1. Polymorphisms for pvdbpII and pvama1 were associated to B-cell epitopes.

CONCLUSIONS

pvdbpII and pvama1 genetic diversity found in Ecuadorian P. vivax was very similar to that encountered in other malaria endemic countries with varying transmission levels and segregated by geographic region. The highest diversity of P. vivax invasion genes in Ecuador was found in the Amazonian region. Although selection appeared to have small effect on pvdbpII and pvrbp1a-2, pvama1 was influenced by significant balancing selection.

Genetic Diversity and Clustering of the Rhoptry Associated Protein-1 of Plasmodium knowlesi from Peninsular Malaysia and Malaysian Borneo

Human infection with simian malaria Plasmodium knowlesi is a cause for concern in Southeast Asian countries, especially in Malaysia. A previous study on Peninsular Malaysia P. knowlesi rhoptry associated protein-1 (PkRAP1) gene has discovered the existence of dimorphism. In this study, genetic analysis of PkRAP1 in a larger number of P. knowlesi samples from Malaysian Borneo was conducted. The PkRAP1 of these P. knowlesi isolates was PCR-amplified and sequenced. The newly obtained PkRAP1 gene sequences (n = 34) were combined with those from the previous study (n = 26) and analysed for polymorphism and natural selection. Sequence analysis revealed a higher genetic diversity of PkRAP1 compared to the previous study. Exon II of the gene had higher diversity (π = 0.0172) than exon I (π = 0.0128). The diversity of the total coding region (π = 0.0167) was much higher than those of RAP1 orthologues such as PfRAP-1 (π = 0.0041) and PvRAP1 (π = 0.00088). Z-test results indicated that the gene was under purifying selection. Phylogenetic tree and haplotype network showed distinct clustering of Peninsular Malaysia and Malaysian Borneo PkRAP1 haplotypes. This geographical-based clustering of PkRAP1 haplotypes provides further evidence of the dimorphism of the gene and possible existence of 2 distinct P. knowlesi lineages in Malaysia.

Why Plasmodium vivax and Plasmodium falciparum are so different? A tale of two clades and their species diversities

The global malaria burden sometimes obscures that the genus Plasmodium comprises diverse clades with lineages that independently gave origin to the extant human parasites. Indeed, the differences between the human malaria parasites were highlighted in the classical taxonomy by dividing them into two subgenera, the subgenus Plasmodium, which included all the human parasites but Plasmodium falciparum that was placed in its separate subgenus, Laverania. Here, the evolution of Plasmodium in primates will be discussed in terms of their species diversity and some of their distinct phenotypes, putative molecular adaptations, and host–parasite biocenosis. Thus, in addition to a current phylogeny using genome-level data, some specific molecular features will be discussed as examples of how these parasites have diverged. The two subgenera of malaria parasites found in primates, Plasmodium and Laverania, reflect extant monophyletic groups that originated in Africa. However, the subgenus Plasmodium involves species in Southeast Asia that were likely the result of adaptive radiation. Such events led to the Plasmodium vivax lineage. Although the Laverania species, including P. falciparum, has been considered to share “avian characteristics,” molecular traits that were likely in the common ancestor of primate and avian parasites are sometimes kept in the Plasmodium subgenus while being lost in Laverania. Assessing how molecular traits in the primate malaria clades originated is a fundamental science problem that will likely provide new targets for interventions. However, given that the genus Plasmodium is paraphyletic (some descendant groups are in other genera), understanding the evolution of malaria parasites will benefit from studying “non-Plasmodium” Haemosporida.

Genetic polymorphism and natural selection of circumsporozoite protein in Myanmar Plasmodium vivax

Background

Circumsporozoite surface protein (CSP) of malaria parasites has been recognized as one of the leading vaccine candidates. Clinical trials of vaccines for vivax malaria incorporating Plasmodium vivax CSP (PvCSP) have demonstrated their effectiveness in preventing malaria, at least in part. However, genetic diversity of pvcsp in the natural population remains a major concern.

Methods

A total of 171 blood samples collected from patients infected with Plasmodium vivax in Myanmar were analysed in this study. The pvcsp was amplified by polymerase chain reaction, followed by cloning and sequencing. Polymorphic characteristics and natural selection of pvcsp population in Myanmar were analysed using DNASTAR, MEGA6 and DnaSP programs. The polymorphic pattern and natural selection of publicly accessible global pvcsp sequences were also comparatively analysed.

Results

Myanmar pvcsp sequences were divided into two subtypes VK210 and VK247 comprising 143 and 28 sequences, respectively. The VK210 subtypes showed higher levels of genetic diversity and polymorphism than the VK247 subtypes. The N-terminal non-repeat region of pvcsp displayed limited genetic variations in the global population. Different patterns of octapeptide insertion (ANKKAEDA in VK210 and ANKKAGDA in VK247) and tetrapeptide repeat motif (GGNA) were identified in the C-terminal region of global pvcsp population. Meanwhile, the central repeat region (CRR) of Myanmar and global pvcsp, both in VK210 and VK247 variants, was highly polymorphic. The high level of genetic diversity in the CRR has been attributed to the different numbers, types and combinations of peptide repeat motifs (PRMs). Interestingly, 27 and 5 novel PRMs were found in Myanmar VK210 and VK247 variants, respectively.

Conclusion

Comparative analysis of the global pvcsp population suggests a complex genetic profile of pvcsp in the global population. These results widen understanding of the genetic make-up of pvcsp in the global P. vivax population and provide valuable information for the development of a vaccine based on PvCSP.

Early transmission of sensitive strain slows down emergence of drug resistance in Plasmodium vivax

The spread of drug resistance of Plasmodium falciparum and Plasmodium vivax parasites is a challenge towards malaria elimination. P. falciparum has shown an early and severe drug resistance in comparison to P. vivax in various countries. In fact, P. vivax differs in its life cycle and treatment in various factors: development and duration of sexual parasite forms differ, symptoms severity are unequal, relapses present only in P. vivax cases and the Artemisinin-based combination therapy (ACT) is only mandatory in P. falciparum cases. We compared the spread of drug resistance for both species through two compartmental models using ordinary differential equations. The model structure describes how sensitive and resistant parasite strains infect a human population treated with antimalarials. We found that an early transmission,i.e., before treatment and low effectiveness of drug coverage, supports the prevalence of sensitive parasites delaying the emergence of resistant P. vivax. These results imply that earlier attention of both symptomatic cases and reservoirs of P. vivax are essential in controlling transmission but also accelerate the spread of drug resistance.

Malaria Molecular Epidemiology: An Evolutionary Genetics Perspective

Malaria is a vector-borne disease that involves multiple parasite species in a variety of ecological settings. However, the parasite species causing the disease, the prevalence of subclinical infections, the emergence of drug resistance, the scale-up of interventions, and the ecological factors affecting malaria transmission, among others, are aspects that vary across areas where malaria is endemic. Such complexities have propelled the study of parasite genetic diversity patterns in the context of epidemiologic investigations. Importantly, molecular studies indicate that the time and spatial distribution of malaria cases reflect epidemiologic processes that cannot be fully understood without characterizing the evolutionary forces shaping parasite population genetic patterns. Although broad in scope, this review in the Microbiology Spectrum Curated Collection: Advances in Molecular Epidemiology highlights the need for understanding population genetic concepts when interpreting parasite molecular data. First, we discuss malaria complexity in terms of the parasite species involved. Second, we describe how molecular data are changing our understanding of malaria incidence and infectiousness. Third, we compare different approaches to generate parasite genetic information in the context of epidemiologically relevant questions related to malaria control. Finally, we describe a few Plasmodium genomic studies as evidence of how these approaches will provide new insights into the malaria disease dynamics. *This article is part of a curated collection.

Exploration of genetic diversity of Plasmodium vivax circumsporozoite protein (Pvcsp) and Plasmodium vivax sexual stage antigen (Pvs25) among North Indian isolates

BACKGROUND

Malaria is one of the important vector-borne diseases with high fatality rates in tropical countries. The pattern of emergence and spread of novel antigenic variants, leading to escape of vaccine-induced immunity might be factors responsible for severe malaria. A high level of polymorphism has been reported among malarial antigens which are under selection pressure imposed by host immunity. There are limited reports available on comparative stage-specific genetic diversity among Plasmodium vivax candidate genes in complicated vivax malaria. The present study was planned to study genetic diversity (Pvcsp and Pvs25) among complicated and uncomplicated P. vivax isolates.

METHODS

Pvcsp and Pvs2-specific PCRs and DNA sequencing were performed on P. vivax PCR positive samples. Genetic diversity was analysed using appropriate software.

RESULTS

The present study was carried out on 143 P. vivax clinical isolates, collected from Postgraduate Institute of Medical Education and Research, Chandigarh. Among the classic and variant types of Pvcsp, the VK210 (99%; 115/116) was found to be predominant in both complicated and uncomplicated group isolates. Out of the various peptide repeat motifs (PRMs) observed, GDRADGQPA (PRM1) and GDRAAGQPA (PRM2) was the most widely distributed among the P. vivax isolates. Whereas among the Pvs25 isolates, 100% of double mutants (E97Q/I130T) in both the complicated (45/45) as well as in the uncomplicated (81/81) group was observed.

CONCLUSION

An analysis of genetic variability enables an understanding of the role of genetic variants in severe vivax malaria.

Effect of low complexity regions within the PvMSP3α block II on the tertiary structure of the protein and implications to immune escape mechanisms

Background

Plasmodium vivax merozoite surface protein 3α (PvMSP3α) is a promising vaccine candidate which has shown strong association with immunogenicity and protectiveness. Its use is however complicated by evolutionary plasticity features which enhance immune evasion. Low complexity regions (LCRs) provide plasticity in surface proteins of Plasmodium species, but its implication in vaccine design remain unexplored. Here population genetic, comparative phylogenetic and structural biology analysis was performed on the gene encoding PvMSP3α.

Results

Three LCRs were found in PvMSP3α block II. Both the predicted tertiary structure of the protein and the phylogenetic trees based on this region were influenced by the presence of the LCRs. The LCRs were mainly B cell epitopes within or adjacent. In addition a repeat motif mimicking one of the B cell epitopes was found within the PvMSP3a block II low complexity region. This particular B cell epitope also featured rampant alanine substitutions which might impair antibody binding.

Conclusion

The findings indicate that PvMSP3α block II possesses LCRs which might confer a strong phenotypic plasticity. The phenomenon of phenotypic plasticity and implication of LCRs in malaria immunology in general and vaccine candidate genes in particular merits further exploration.

Electronic supplementary material

The online version of this article (10.1186/s12900-019-0104-0) contains supplementary material, which is available to authorized users.

Genetic polymorphism and natural selection in the C-terminal 42 kDa region of merozoite surface protein-1 (MSP-1) among Plasmodium knowlesi samples from Malaysia

Background

The merozoite surface protein-1 (MSP-1) gene encodes for a leading malaria vaccine candidate antigen. However, its extensive polymorphic nature represents a major obstacle to the development of a protective vaccine. Previously, a pilot study was carried out to explore the sequence variation of the C-terminal 42 kDa fragment within P. knowlesi MSP-1 gene (PkMSP-1₄₂) based on 12 clinical samples; however, further study on an adequate sample size is vital in estimating the genetic diversity of the parasite population.

Methods

In the present study, we included a larger sample size of P. knowlesi (83 samples) covering eight states of Malaysia to determine the genetic polymorphism, natural selection and haplotype groups of the gene fragment coding PkMSP-1₄₂. The region flanking PkMSP-1₄₂ was amplified by PCR and directly sequenced. Genetic diversity, haplotype diversity, population genetic differentiation and natural selection were determined in order to study the polymorphic characteristic of PkMSP-1₄₂.

Results

A high level of genetic diversity (Hd = 0.970 ± 0.007; л = 0.01079 ± 0.00033) was observed among the 83 P. knowlesi samples, confirming the extensive genetic polymorphism exhibited among the P. knowlesi population found in Malaysia. A total of 18 distinct haplotypes with 17 amino acid changes were identified, whereby 15 were new haplotypes. High population differentiation values were observed within samples from Peninsular Malaysia and Malaysian Borneo. The 42 kDa fragments of P. knowlesi from Malaysian Borneo were found to be acting on balancing selection whilst purifying selection was suggested to act on isolates from Peninsular Malaysia. The separation of PkMSP-1₄₂ haplotypes into two main groups based on geographical separation has further supported the existence of two distinct P. knowlesi lineages.

Conclusions

A high level of genetic diversity was observed among PkMSP-1₄₂ in Malaysia, whereby most of the polymorphisms were found within the 33 kDa region. Taken together, these data will be useful in order to understand the nature of P. knowlesi population in Malaysia as well as the design and development of a MSP-1₄₂ based knowlesi malaria vaccine.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3234-5) contains supplementary material, which is available to authorized users.

The role of natural selection in shaping genetic variation in a promising Chagas disease drug target: Trypanosoma cruzi trans-sialidase

Rational drug design creates innovative therapeutics based on knowledge of the biological target to provide more effective and responsible therapeutics. Chagas disease, endemic throughout Latin America, is caused by Trypanosoma cruzi, a protozoan parasite. Current therapeutics are problematic with widespread calls for new approaches. Researchers are using rational drug design for Chagas disease and one target receiving considerable attention is the T. cruzi trans-sialidase protein (TcTS). In T. cruzi, trans-sialidase catalyzes the transfer of sialic acid from a mammalian host to coat the parasite surface membrane and avoid immuno-detection. However, the role of TcTS in pathology variance among and within genetic variants of the parasite is not well understood despite numerous studies. Previous studies reported the crystalline structure of TcTS and the TS protein structure in other trypanosomes where the enzyme is often inactive. However, no study has examined the role of natural selection in genetic variation in TcTS. To understand the role of natural selection in TcTS DNA sequence and protein variation, we examined a 471 bp portion of the TcTS gene from 48 T. cruzi samples isolated from insect vectors. Because there may be multiple parasite genotypes infecting one insect and there are multiple copies of TcTS per parasite genome, all 48 sequences had multiple polymorphic bases. To resolve these polymorphisms, we examined cloned sequences from two insect vectors. The data are analyzed to understand the role of natural selection in shaping genetic variation in TcTS and interpreted in light of the possible role of TcTS as a drug target. The analysis highlights negative or purifying selection on three amino acids previously shown to be important in TcTS transfer activity. One amino acid in particular, Tyr342, is a strong candidate for a drug target because it is under negative selection and amino acid substitutions inactivate TcTS transfer activity. AUTHOR SUMMARY: Chagas disease is caused by the protozoan parasite Trypanosoma cruzi and transmitted to humans and other mammals primarily by Triatomine insects. Being endemic in many South and Central American countries and affecting millions of people the need for new more effective and safe therapies is evident. Here, we examine genetic variation and natural selection on DNA (471 bp) and amino acid (157 aa) sequence data of the T. cruzi trans-sialdiase (TcTS) protein, often suggested as a candidate for rational drug design. In our surveyed region of the protein there were five amino acid residues that have been shown to be integral to the function of TcTS. We found that three were under strong negative selection making them ideal candidates for drug design; however, one was under balancing selection and should be avoided as a drug target. Our study provides new information into identifying potential targets for a new Chagas drug.

Analysis of polymorphisms and selective pressures on ama1 gene in Plasmodium knowlesi isolates from Sabah, Malaysia

The apical membrane antigen-1 (AMA-1) of Plasmodium spp. is a merozoite surface antigen that is essential for the recognition and invasion of erythrocytes. Polymorphisms occurring in this surface antigen will cause major obstacles in developing effective malaria vaccines based on AMA-1. The objective of this study was to characterize ama1 gene in Plasmodium knowlesi isolates from Sabah. DNA was extracted from blood samples collected from Keningau, Kota Kinabalu and Kudat. The Pkama1 gene was amplified using nested PCR and subjected to bidirectional sequencing. Analysis of DNA sequence revealed that most of the nucleotide polymorphisms were synonymous and concentrated in domain I of PkAMA-1. Forteen haplotypes were identified based on amino acid variations and haplotype K5 was the most common haplotype. d_N/d_S ratios implied that purifying selection was prevalent in Pkama1 gene. Fu and Li's D and F values further provided evidence of negative selection acting on domain II of Pkama1. Lownucleotide diversitywas also detected for the Pkama1 sequences,which is similar to reports on Pkama1 from Peninsular Malaysia and Sarawak. The presence of purifying selection and low nucleotide diversity indicated that domain II of Pkama1 can be used as a target for vaccine development.

Temporal genetic changes in Plasmodium vivax apical membrane antigen 1 over 19 years of transmission in southern Mexico

Background

Mexico advanced to the pre-elimination phase in 2009 due to a significant reduction in malaria cases, and since 2000, Plasmodium vivax is the only species transmitted. During the last two decades, malaria transmission has been mostly local and isolated to a few regions. It is important to gain further insights into the impact of control measures on the parasite population structure. Hence, the aim of the current study was to determine detailed changes in P. vivax genetic diversity and population structure based on analysing the gene that encodes the apical membrane antigen 1 (pvama1). This analysis covered from control to pre-elimination (1993–2011) in a hypo-endemic region in southern Mexico.

Results

The 213 pvama1_I-II sequences presently analysed were grouped into six periods of three years each. They showed low genetic diversity, with 15 haplotypes resolved. Among the DNA sequences, there was a gradual decrease in genetic diversity, the number of mixed genotype infections and the intensity of positive selection, in agreement with the parallel decline in malaria cases. At the same time, linkage disequilibrium (R²) increased. The three-dimensional haplotype network revealed that pvama1_I-II haplotypes were separated by 1–11 mutational steps, and between one another by 0–3 unsampled haplotypes. In the temporal network, seven haplotypes were detected in at least two of the six-time layers, and only four distinct haplotypes were evidenced in the pre-elimination phase. Structure analysis indicated that three subpopulations fluctuated over time. Only 8.5% of the samples had mixed ancestry. In the pre-elimination phase, subpopulation P1 was drastically reduced, and the admixture was absent.

Conclusions

The results suggest that P. vivax in southern Mexico evolved based on local adaptation into three “pseudoclonal” subpopulations that diversified at the regional level and persisted over time, although with varying frequency. Control measures and climate events influenced the number of malaria cases and the genetic structure. The sharp decrease in parasite diversity and other related genetic parameters during the pre-elimination phase suggests that malaria elimination is possible in the near future. These results are useful for epidemiological surveillance.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2156-y) contains supplementary material, which is available to authorized users.

Genetic diversity and selection of three nuclear genes in Schistosoma japonicum populations

BACKGROUND

The blood fluke, Schistosoma japonicum still causes severe disease in China, the Philippines and Indonesia. Although there have been some studies the molecular epidemiology of this persistent and harmful parasite, few have explored the possibility and implications of selection in S. japonicum populations.

METHODS

We analyzed diversity and looked for evidence of selection at three nuclear genes (SjIpp2, SjFabp and SjT22.6) in 13 S. japonicum populations.

RESULTS

SjT22.6 was found to exhibit high nucleotide diversity and was under positive selection in the mountainous region of mainland China. As a tegumental protein, its secondary and tertiary structure differed between S. japonicum strains from the mountainous and lakes regions. In contrast, SjIpp2 and SjFabp had relatively low levels of nucleotide diversity and did not show significant departure from neutrality.

CONCLUSIONS

As a tegument-associated antigen-encoding gene of S. japonicum, SjT22.6 has high nucleotide diversity and appears to be under positive selection in the mountainous region of mainland China.

Genetic clustering and polymorphism of the merozoite surface protein-3 of Plasmodium knowlesi clinical isolates from Peninsular Malaysia

Background

The simian malaria parasite Plasmodium knowlesi has been reported to cause significant numbers of human infection in South East Asia. Its merozoite surface protein-3 (MSP3) is a protein that belongs to a multi-gene family of proteins first found in Plasmodium falciparum. Several studies have evaluated the potential of P. falciparum MSP3 as a potential vaccine candidate. However, to date no detailed studies have been carried out on P. knowlesi MSP3 gene (pkmsp3). The present study investigates the genetic diversity, and haplotypes groups of pkmsp3 in P. knowlesi clinical samples from Peninsular Malaysia.

Methods

Blood samples were collected from P. knowlesi malaria patients within a period of 4 years (2008–2012). The pkmsp3 gene of the isolates was amplified via PCR, and subsequently cloned and sequenced. The full length pkmsp3 sequence was divided into Domain A and Domain B. Natural selection, genetic diversity, and haplotypes of pkmsp3 were analysed using MEGA6 and DnaSP ver. 5.10.00 programmes.

Results

From 23 samples, 48 pkmsp3 sequences were successfully obtained. At the nucleotide level, 101 synonymous and 238 non-synonymous mutations were observed. Tests of neutrality were not significant for the full length, Domain A or Domain B sequences. However, the dN/dS ratio of Domain B indicates purifying selection for this domain. Analysis of the deduced amino acid sequences revealed 42 different haplotypes. Neighbour Joining phylogenetic tree and haplotype network analyses revealed that the haplotypes clustered into two distinct groups.

Conclusions

A moderate level of genetic diversity was observed in the pkmsp3 and only the C-terminal region (Domain B) appeared to be under purifying selection. The separation of the pkmsp3 into two haplotype groups provides further evidence of the existence of two distinct P. knowlesi types or lineages. Future studies should investigate the diversity of pkmsp3 among P. knowlesi isolates in North Borneo, where large numbers of human knowlesi malaria infection still occur.

Electronic supplementary material

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

Evolution of the Transmission-Blocking Vaccine Candidates Pvs28 and Pvs25 in Plasmodium vivax: Geographic Differentiation and Evidence of Positive Selection

Transmission-blocking (TB) vaccines are considered an important tool for malaria control and elimination. Among all the antigens characterized as TB vaccines against Plasmodium vivax, the ookinete surface proteins Pvs28 and Pvs25 are leading candidates. These proteins likely originated by a gene duplication event that took place before the radiation of the known Plasmodium species to primates. We report an evolutionary genetic analysis of a worldwide sample of pvs28 and pvs25 alleles. Our results show that both genes display low levels of genetic polymorphism when compared to the merozoite surface antigens AMA-1 and MSP-1; however, both ookinete antigens can be as polymorphic as other merozoite antigens such as MSP-8 and MSP-10. We found that parasite populations in Asia and the Americas are geographically differentiated with comparable levels of genetic diversity and specific amino acid replacements found only in the Americas. Furthermore, the observed variation was mainly accumulated in the EGF2- and EGF3-like domains for P. vivax in both proteins. This pattern was shared by other closely related non-human primate parasites such as Plasmodium cynomolgi, suggesting that it could be functionally important. In addition, examination with a suite of evolutionary genetic analyses indicated that the observed patterns are consistent with positive natural selection acting on Pvs28 and Pvs25 polymorphisms. The geographic pattern of genetic differentiation and the evidence for positive selection strongly suggest that the functional consequences of the observed polymorphism should be evaluated during development of TBVs that include Pvs25 and Pvs28.

Plasmodium vivax is the most prevalent human malarial parasite outside Africa. The fact that patients can relapse due to the parasite dormant liver stages, among other biologic and epidemiologic characteristics of vivax malaria, facilitates the persistence of the disease in many endemic areas. These challenges have fueled the search for new control tools, including transmission blocking (TB) vaccines targeting the parasite sexual stages. Here we study the genetic diversity of two major TB vaccine antigens, Pvs25 and Pvs28. We show that these genes are relatively conserved worldwide but still harbor diversity that is not evenly distributed across the genes. These patterns are shared by the same proteins in closely related parasite species suggesting their functional importance. We also identify strong geographic differentiation between the circulating variants found in Asia and the Americas. Finally, evolutionary genetic analyses indicate that the observed variation in both genes could be maintained by natural selection. Thus, these polymorphisms may confer an adaptive advantage to the parasite. These results indicate that the genetic variation found in these genes and their geographic distribution should be considered by vaccine developers.

Genetic diversity among Plasmodium vivax isolates along the Thai-Myanmar border of Thailand

Background

Knowledge of the population genetics and transmission dynamics of Plasmodium vivax is crucial in predicting the emergence of drug resistance, relapse pattern and novel parasite phenotypes, all of which are relevant to the control of vivax infections. The aim of this study was to analyse changes in the genetic diversity of P. vivax genes from field isolates collected at different times along the Thai–Myanmar border.

Methods

Two hundred and fifty-four P. vivax isolates collected during two periods 10 years apart along the Thai–Myanmar border were analysed. The parasites were genotyped by nested-PCR and PCR–RFLP targeting selected polymorphic loci of Pvmsp1, Pvmsp3α and Pvcsp genes.

Results

The total number of distinguishable allelic variants observed for Pvcsp, Pvmsp1, and Pvmsp3α was 17, 7 and 3, respectively. High genetic diversity was observed for Pvcsp (H_E = 0.846) and Pvmsp1 (H_E = 0.709). Of the 254 isolates, 4.3 and 14.6 % harboured mixed Pvmsp1 and Pvcsp genotypes with a mean multiplicity of infection (MOI) of 1.06 and 1.15, respectively. The overall frequency of multiple genotypes was 16.9 %. When the frequencies of allelic variants of each gene during the two distinct periods were analysed, significant differences were noted for Pvmsp1 (P = 0.018) and the Pvcsp (P = 0.033) allelic variants.

Conclusion

Despite the low malaria transmission levels in Thailand, P. vivax population exhibit a relatively high degree of genetic diversity along the Thai–Myanmar border of Thailand, in particular for Pvmsp1 and Pvcsp, with indication of geographic and temporal variation in frequencies for some variants. These results are of relevance to monitoring the emergence of drug resistance and to the elaboration of measures to control vivax malaria.

Genetic diversity and natural selection in the rhoptry-associated protein 1 (RAP-1) of recent Plasmodium knowlesi clinical isolates from Malaysia

BACKGROUND

The Plasmodium rhoptry-associated protein 1 (RAP-1) plays a role in the formation of the parasitophorous vacuole following the parasite's invasion of red blood cells. Although there is some evidence that the protein is recognized by the host's immune system, study of Plasmodium falciparum RAP-1 (PfRAP-1) suggests that it is not under immune pressure. A previous study on five old (1953-1962) P. knowlesi strains suggested that RAP-1 has limited genetic polymorphism and might be under negative selection. In the present study, 30 recent P. knowlesi isolates were studied to obtain a better insight into the polymorphism and natural selection of PkRAP-1.

METHODS

Blood samples from 30 knowlesi malaria patients were used. These samples were collected between 2010 and 2014. The PkRAP-1 gene, which contains two exons, was amplified by PCR, cloned into Escherichia coli and sequenced. Genetic diversity and phylogenetic analyses were performed using MEGA6 and DnaSP ver. 5.10.00 programs.

RESULTS

Thirty PkRAP-1 sequences were obtained. The nucleotide diversity (π) of exons 1, 2 and the total coding region (0.00915, 0.01353 and 0.01298, respectively) were higher than those of the old strains. Further analysis revealed a lower rate of non-synonymous (dN) than synonymous (dS) mutations, suggesting negative (purifying) selection of PkRAP-1. Tajima's D test and Fu and Li's D test values were not significant. At the amino acid level, 22 haplotypes were established with haplotype H7 having the highest frequency (7/34, 20.5 %). In the phylogenetic analysis, two distinct haplotype groups were observed. The first group contained the majority of the haplotypes, whereas the second had fewer haplotypes.

CONCLUSIONS

The present study found higher genetic polymorphism in the PkRAP-1 gene than the polymorphism level reported in a previous study. This observation may stem from the difference in sample size between the present (n = 30) and the previous (n = 5) study. Synonymous and non-synonymous mutation analysis indicated purifying (negative) selection of the gene. The separation of PkRAP-1haplotypes into two groups provides further evidence to the postulation of two distinct P. knowlesi types or lineages.

Diversity and evolutionary genetics of the three major Plasmodium vivax merozoite genes participating in reticulocyte invasion in southern Mexico

Background

Reported malaria cases in the Americas had been reduced to about one-half million by 2012. To advance towards elimination of this disease, it is necessary to gain insights into how the malaria parasite is evolving, including the emergence, spread and persistence of new haplotypes in affected regions. In here, the genetic diversity of the three major P. vivax merozoite genes was analyzed.

Methods

From P. vivax-infected blood samples obtained in southern Mexico (SMX) during 2006–2007, nucleotide sequences were achieved for: the 42 kDa carboxyl fragment of the merozoite surface protein-1 (msp1₄₂), domains I-II of the apical membrane antigen-1 (ama1_I-II), and domain II of the Duffy binding protein (dbp_II). Gene polymorphism was examined and haplotype networks were developed to depict parasite relationships in SMX. Then genetic diversity, recombination and natural selection were analyzed and the degree of differentiation was determined as F_ST values.

Results

The diversity of P. vivax merozoite genes in SMX was less than that of parasites from other geographic origins, with dbp_II < ama1_I-II < msp1₄₂. Ama1_I-II and msp1₄₂ exposed the more numerous haplotypes exclusive to SMX. While, all dbp_II haplotypes from SMX were separated from one to three mutational steps, the networks of ama1_I-II and msp1₄₂ were more complex; loops and numerous mutational steps were evidenced, likely due to recombination. Sings of local diversification were more evident for msp1₄₂. Sixteen combined haplotypes were determined; one of these haplotypes not detected in 2006 was highly frequent in 2007. The Rm value was higher for msp1₄₂than for ama1_I-II, being insignificant for dbp_II. The dN-dS value was highly significant for ama1_I-II and lesser so for dbp_II. The F_ST values were higher for dbp_II than msp1₄₂, and very low for ama1_I-II.

Conclusions

In SMX, P. vivax ama1_I-II, dbp_II and msp1₄₂ demonstrated limited diversity, and exhibited a differentiated parasite population. The results suggest that differential intensities of selective forces are operating on these gene fragments, and probably related to their timing, length of exposure and function during reticulocyte adhesion and invasion. Therefore, these finding are essential for mono and multivalent vaccine development and for epidemiological surveillance.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-1266-7) contains supplementary material, which is available to authorized users.

Genetic polymorphism in domain I of the apical membrane antigen-1 among Plasmodium knowlesi clinical isolates from Peninsular Malaysia

The simian malaria parasite Plasmodium knowlesi is now recognized as a species that can cause human malaria. The first report of large scale human knowlesi malaria was in 2004 in Malaysia Borneo. Since then, hundreds of human knowlesi malaria cases have been reported in Southeast Asia. The present study investigates the genetic polymorphism of P. knowlesi DI domain of the apical membrane antigen-1 (AMA-1), a protein considered as a promising vaccine candidate for malaria. The DI domain of AMA-1 gene of P. knowlesi clinical isolates from Peninsular Malaysia was amplified by PCR, cloned into Escherichia coli, then sequenced and analysed. Ninety-seven DI domain sequences were obtained. Comparison at the nucleotide level against P. knowlesi strain H as reference sequence showed 21 synonymous and 25 nonsynonymous mutations. Nonetheless, nucleotide sequence analysis revealed low genetic diversity of the DI domain, and it was under purifying (negative) selection. At the amino acid level, 26 different haplotypes were identified and 2 were predominant haplotypes (H1, H2) with high frequencies. Phylogenetic analysis revealed that the 26 haplotypes could be clustered into 2 distinct groups (I and II). Members of the groups were basically derived from haplotypes H1 and H2, respectively.

Population genetic structure and natural selection of apical membrane antigen-1 in Plasmodium vivax Korean isolates

Background

Plasmodium vivax apical membrane antigen-1 (PvAMA-1) is a leading candidate antigen for blood stage malaria vaccine. However, antigenic variation is a major obstacle in the development of an effective vaccine based on this antigen. In this study, the genetic structure and the effect of natural selection of PvAMA-1 among Korean P. vivax isolates were analysed.

Methods

Blood samples were collected from 66 Korean patients with vivax malaria. The entire PvAMA-1 gene was amplified by polymerase chain reaction and cloned into a TA cloning vector. The PvAMA-1 sequence of each isolate was sequenced and the polymorphic characteristics and effect of natural selection were analysed using the DNASTAR, MEGA4, and DnaSP programs.

Results

Thirty haplotypes of PvAMA-1, which were further classified into seven different clusters, were identified in the 66 Korean P. vivax isolates. Domain II was highly conserved among the sequences, but substantial nucleotide diversity was observed in domains I and III. The difference between the rates of non-synonymous and synonymous mutations suggested that the gene has evolved under natural selection. No strong evidence indicating balancing or positive selection on PvAMA-1 was identified. Recombination may also play a role in the resulting genetic diversity of PvAMA-1.

Conclusions

This study is the first comprehensive analysis of nucleotide diversity across the entire PvAMA-1 gene using a single population sample from Korea. Korean PvAMA-1 had limited genetic diversity compared to PvAMA-1 in global isolates. The overall pattern of genetic polymorphism of Korean PvAMA-1 differed from other global isolates and novel amino acid changes were also identified in Korean PvAMA-1. Evidences for natural selection and recombination event were observed, which is likely to play an important role in generating genetic diversity across the PvAMA-1. These results provide useful information for the understanding the population structure of P. vivax circulating in Korea and have important implications for the design of a vaccine incorporating PvAMA-1.

Molecular Evolution of PvMSP3α Block II in Plasmodium vivax from Diverse Geographic Origins

Block II of Plasmodium vivax merozoite surface protein 3α (PvMSP3α) is conserved and has been proposed as a potential candidate for a malaria vaccine. The present study aimed to compare sequence diversity in PvMSP3a block II at a local microgeographic scale in a village as well as from larger geographic regions (countries and worldwide). Blood samples were collected from asymptomatic carriers of P. vivax in a village at the western border of Thailand and PvMSP3α was amplified and sequenced. For population genetic analysis, 237 PvMSP3α block II sequences from eleven P. vivax endemic countries were analyzed. PvMSP3α sequences from 20 village-level samples revealed two length variant types with one type containing a large deletion in block I. In contrast, block II was relatively conserved; especially, some non-synonymous mutations were extensively shared among 11 parasite populations. However, the majority of the low-frequency synonymous variations were population specific. The conserved pattern of nucleotide diversity in block II sequences was probably due to functional/structural constraints, which were further supported by the tests of neutrality. Notably, a small region in block II that encodes a predicted B cell epitope was highly polymorphic and showed signs of balancing selection, signifying that this region might be influenced by the immune selection and may serve as a starting point for designing multi-antigen/stage epitope based vaccines against this parasite.

Genome-wide patterns of genetic polymorphism and signatures of selection in Plasmodium vivax

Plasmodium vivax is the most prevalent human malaria parasite outside of Africa. Yet, studies aimed to identify genes with signatures consistent with natural selection are rare. Here, we present a comparative analysis of the pattern of genetic variation of five sequenced isolates of P. vivax and its divergence with two closely related species, Plasmodium cynomolgi and Plasmodium knowlesi, using a set of orthologous genes. In contrast to Plasmodium falciparum, the parasite that causes the most lethal form of human malaria, we did not find significant constraints on the evolution of synonymous sites genome wide in P. vivax. The comparative analysis of polymorphism and divergence across loci allowed us to identify 87 genes with patterns consistent with positive selection, including genes involved in the “exportome” of P. vivax, which are potentially involved in evasion of the host immune system. Nevertheless, we have found a pattern of polymorphism genome wide that is consistent with a significant amount of constraint on the replacement changes and prevalent negative selection. Our analyses also show that silent polymorphism tends to be larger toward the ends of the chromosomes, where many genes involved in antigenicity are located, suggesting that natural selection acts not only by shaping the patterns of variation within the genes but it also affects genome organization.

PCR/RFLP-based analysis of genetically distinct Plasmodium vivax population of Pvmsp-3α and Pvmsp-3β genes in Pakistan

Background

Plasmodium vivax is one of the widespread human malarial parasites accounting for 75% of malaria epidemics. However, there is no baseline information about the status and nature of genetic variation of Plasmodium species circulating in various parts of Pakistan. The present study was aimed at observing the molecular epidemiology and genetic variation of Plasmodium vivax by analysing its merozoite surface protein-3α (msp-3α) and merozoite surface protein-3β (msp-3β) genes, by using suballele, species-specific, combined nested PCR/RFLP detection techniques.

Methods

A total of 230 blood samples from suspected subjects tested slide positive for vivax malaria were collected from Punjab, Sindh, Khyber Pakhtunkhwa, and Balochistan during the period May 2012 to December 2013. Combined nested PCR/RFLP technique was conducted using Pvmsp-3α and Pvmsp-3β genetic markers to detect extent of genetic variation in clinical isolates of P. vivax in the studied areas of Pakistan.

Results

By PCR, P. vivax, 202/230 (87.82%), was found to be widely distributed in the studied areas. PCR/RFLP analysis showed a high range of allelic variations for both msp-3α and msp-3β genetic markers of P. vivax, i.e., 21 alleles for msp-3α and 19 for msp-3β. Statistically a significant difference (p ≤ 0.05) was observed in the genetic diversity of the suballelic variants of msp-3α and msp-3β genes of P. vivax.

Conclusion

It is concluded that P. vivax populations are highly polymorphic and diverse allelic variants of Pvmsp-3α and Pvmsp-3β are present in Pakistan.

An unstable Th epitope of P. falciparum fosters central memory T cells and anti-CS antibody responses

Malaria is transmitted by Plasmodium-infected anopheles mosquitoes. Widespread resistance of mosquitoes to insecticides and resistance of parasites to drugs highlight the urgent need for malaria vaccines. The most advanced malaria vaccines target sporozoites, the infective form of the parasite. A major target of the antibody response to sporozoites are the repeat epitopes of the circumsporozoite (CS) protein, which span almost one half of the protein. Antibodies to these repeats can neutralize sporozoite infectivity. Generation of protective antibody responses to the CS protein (anti-CS Ab) requires help by CD4 T cells. A CD4 T cell epitope from the CS protein designated T* was previously identified by screening T cells from volunteers immunized with irradiated P. falciparum sporozoites. The T* sequence spans twenty amino acids that contains multiple T cell epitopes restricted by various HLA alleles. Subunit malaria vaccines including T* are highly immunogenic in rodents, non-human primates and humans. In this study we characterized a highly conserved HLA-DRβ1*04:01 (DR4) restricted T cell epitope (QNT-5) located at the C-terminus of T*. We found that a peptide containing QNT-5 was able to elicit long-term anti-CS Ab responses and prime CD4 T cells in HLA-DR4 transgenic mice despite forming relatively unstable MHC-peptide complexes highly susceptible to HLA-DM editing. We attempted to improve the immunogenicity of QNT-5 by replacing the P1 anchor position with an optimal tyrosine residue. The modified peptide QNT-Y formed stable MHC-peptide complexes highly resistant to HLA-DM editing. Contrary to expectations, a linear peptide containing QNT-Y elicited almost 10-fold lower long-term antibody and IFN-γ responses compared to the linear peptide containing the wild type QNT-5 sequence. Some possibilities regarding why QNT-5 is more effective than QNT-Y in inducing long-term T cell and anti-CS Ab when used as vaccine are discussed.

The origin and diversification of the merozoite surface protein 3 (msp3) multi-gene family in Plasmodium vivax and related parasites

The genus Plasmodium is a diversified group of parasites with more than 200 known species that includes those causing malaria in humans. These parasites use numerous proteins in a complex process that allows them to invade the red blood cells of their vertebrate hosts. Many of those proteins are part of multi-gene families; one of which is the merozoite surface protein-3 (msp3) family. The msp3 multi-gene family is considered important in the two main human parasites, Plasmodium vivax and Plasmodium falciparum, as its paralogs are simultaneously expressed in the blood stage (merozoite) and are immunogenic. There are large differences among Plasmodium species in the number of paralogs in this family. Such differences have been previously explained, in part, as adaptations that allow the different Plasmodium species to invade their hosts. To investigate this, we characterized the array containing msp3 genes among several Plasmodium species, including P. falciparum and P. vivax. We first found no evidence indicating that the msp3 family of P. falciparum was homologous to that of P. vivax. Subsequently, by focusing on the diverse clade of nonhuman primate parasites to which P. vivax is closely related, where homology was evident, we found no evidence indicating that the interspecies variation in the number of paralogs was an adaptation related to changes in host range or host switches. Overall, we hypothesize that the evolution of the msp3 family in P. vivax is consistent with a model of multi-allelic diversifying selection where the paralogs may have functionally redundant roles in terms of increasing antigenic diversity. Thus, we suggest that the expressed MSP3 proteins could serve as "decoys", via antigenic diversity, during the critical process of invading the host red blood cells.

Genetic diversity, haplotypes and allele groups of Duffy binding protein (PkDBPαII) of Plasmodium knowlesi clinical isolates from Peninsular Malaysia

Background

The monkey malaria parasite Plasmodium knowlesi is now recognized as the fifth species of Plasmodium that can cause human malaria. Like the region II of the Duffy binding protein of P. vivax (PvDBPII), the region II of the P. knowlesi Duffy binding protein (PkDBPαII) plays an essential role in the parasite’s invasion into the host’s erythrocyte. Numerous polymorphism studies have been carried out on PvDBPII, but none has been reported on PkDBPαII. In this study, the genetic diversity, haplotyes and allele groups of PkDBPαII of P. knowlesi clinical isolates from Peninsular Malaysia were investigated.

Methods

Blood samples from 20 knowlesi malaria patients and 2 wild monkeys (Macaca fascicularis) were used. These samples were collected between 2010 and 2012. The PkDBPαII region of the isolates was amplified by PCR, cloned into Escherichia coli, and sequenced. The genetic diversity, natural selection and haplotypes of PkDBPαII were analysed using MEGA5 and DnaSP ver. 5.10.00 programmes.

Results

Fifty-three PkDBPαII sequences from human infections and 6 from monkeys were obtained. Comparison at the nucleotide level against P. knowlesi strain H as reference sequence showed 52 synonymous and 76 nonsynonymous mutations. Analysis on the rate of these mutations indicated that PkDBPαII was under purifying (negative) selection. At the amino acid level, 36 different PkDBPαII haplotypes were identified. Twelve of the 20 human and 1 monkey blood samples had mixed haplotype infections. These haplotypes were clustered into 2 distinct allele groups. The majority of the haplotypes clustered into the large dominant group.

Conclusions

Our present study is the first to report the genetic diversity and natural selection of PkDBPαII. Hence, the haplotypes described in this report can be considered as novel. Although a high level of genetic diversity was observed, the PkDBPαII appeared to be under purifying selection. The distribution of the haplotypes was skewed, with one dominant major and one minor group. Future study should investigate PkDBPαII of P. knowlesi from Borneo, which hitherto has recorded the highest number of human knowlesi malaria.

Polymorphic patterns of the merozoite surface protein-3β in Korean isolates of Plasmodium vivax

BACKGROUND

The merozoite surface protein-3β of Plasmodium vivax (PvMSP-3β) is one of the candidate antigens for blood stage malaria vaccine development. The polymorphisms in PvMSP-3β have been reported in certain P. vivax isolates. However, the diversity of PvMSP-3β throughout its global distribution has not been well understood. In this study, the genetic diversity and the effects of natural selection in PvMSP-3β among P. vivax Korean isolates were analysed.

METHODS

Blood samples were collected from 95 patients with vivax malaria in Korea. The region flanking full-length PvMSP-3β was amplified by polymerase chain reaction and cloned into a TA cloning vector. The PvMSP-3β sequence of each isolate was determined and the polymorphic characteristics and effects of natural selection were analysed using the DNASTAR, MEGA4, and DnaSP programs.

RESULTS

Five different subtypes of PvMSP-3β were identified based on single nucleotide polymorphisms (SNPs), insertions, and deletions. Although a high level of sequence diversity was observed in the PvMSP-3β gene, the coiled-coil tertiary structure of the PvMSP-3β protein was well conserved in all of the sequences. The PvMSP-3β of Korean isolates is under natural selection. DNA polymerase slippage and intragenic recombination likely contributed to PvMSP-3β diversity in Korean P. vivax isolates.

CONCLUSIONS

The PvMSP-3β of Korean P. vivax isolates displayed polymorphisms, with SNPs, insertions and deletions scattered throughout of the gene. These results of parasite heterogeneity are relevant to the development of a PvMSP-3β based vaccine against P. vivax and the implementation of malaria control programmes in Korea.

On the road to eliminate malaria in Sri Lanka: lessons from history, challenges, gaps in knowledge and research needs

Malaria is one of the most important tropical diseases that has caused devastation throughout the history of mankind. Malaria eradication programmes in the past have had many positive effects but failed to wipe out malaria from most tropical countries, including Sri Lanka. Encouraged by the impressive levels of reduction in malaria case numbers during the past decade, Sri Lanka has launched a programme to eliminate malaria by year 2014. This article reviews the historical milestones associated with the malaria eradication programme that failed subsequently and the events that led to the launch of the ongoing malaria elimination plans at national-level and its strategies that are operational across the entire country. The existing gaps in knowledge are also discussed together with the priority areas for research to fill in these gaps that are posing as challenges to the envisaged goal of wiping out malaria from this island nation.

The evolution and diversity of a low complexity vaccine candidate, merozoite surface protein 9 (MSP-9), in Plasmodium vivax and closely related species

The merozoite surface protein-9 (MSP-9) has been considered a target for an anti-malarial vaccine since it is one of many proteins involved in the erythrocyte invasion, a critical step in the parasite life cycle. Orthologs encoding this antigen have been found in all known species of Plasmodium parasitic to primates. In order to characterize and investigate the extent and maintenance of MSP-9 genetic diversity, we analyzed DNA sequences of the following malaria parasite species: Plasmodium falciparum, Plasmodium reichenowi, Plasmodium chabaudi, Plasmodium yoelii, Plasmodium berghei, Plasmodium coatneyi, Plasmodium gonderi, Plasmodium knowlesi, Plasmodium inui, Plasmodium simiovale, Plasmodium fieldi, Plasmodium cynomolgi and Plasmodium vivax and evaluated the signature of natural selection in all MSP-9 orthologs. Our findings suggest that the gene encoding MSP-9 is under purifying selection in P. vivax and closely related species. We further explored how selection affected different regions of MSP-9 by comparing the polymorphisms in P. vivax and P. falciparum, and found contrasting patterns between these two species that suggest differences in functional constraints. This observation implies that the MSP-9 orthologs in human parasites may interact differently with the host immune response. Thus, studies carried out in one species cannot be directly translated into the other.

Merozoite surface protein-3 alpha as a genetic marker for epidemiologic studies in Plasmodium vivax: a cautionary note

Background

Plasmodium vivax is the most widespread of the human malaria parasites in terms of geography, and is thought to present unique challenges to local efforts aimed at control and elimination. Parasite molecular markers can provide much needed data on P. vivax populations, but few such markers have been critically evaluated. One marker that has seen extensive use is the gene encoding merozoite surface protein 3-alpha (MSP-3α), a blood-stage antigen known to be highly variable among P. vivax isolates. Here, a sample of complete msp-3α gene sequences is analysed in order to assess its utility as a molecular marker for epidemiologic investigations.

Methods

Amplification, cloning and sequencing of additional P. vivax isolates from different geographic locations, including a set of Venezuelan field isolates (n = 10), yielded a sample of 48 complete msp-3α coding sequences. Characterization of standard population genetic measures of diversity, phylogenetic analysis, and tests for recombination were performed. This allowed comparisons to patterns inferred from the in silico simulation of a polymerase chain reaction restriction fragment length polymorphism (PCR-RFLP) protocol used widely.

Results

The larger sample of MSP-3α diversity revealed incongruence between the observed levels of nucleotide polymorphism, which were high in all populations, and the pattern of PCR-RFLP haplotype diversity. Indeed, PCR-RFLP haplotypes were not informative of a population’s genetic diversity and identical haplotypes could be produced from analogous bands in the commonly used protocol. Evidence of frequent and variable insertion-deletion mutations and recurrent recombination between MSP-3α haplotypes complicated the inference of genetic diversity patterns and reduced the phylogenetic signal.

Conclusions

The genetic diversity of P. vivax msp-3α involves intragenic recombination events. Whereas the high genetic diversity of msp-3α makes it a promising marker for some epidemiological applications, the ability of msp-3α PCR-RFLP analysis to accurately track parasites is limited. Local studies of the circulating alleles are needed before implementing PCR-RFLP approaches. Furthermore, evidence from the global sample analysed here suggests such msp-3α PCR-RFLP methods are not suitable for broad geographic studies or tracking parasite populations for an extended period of time.

Genomics, population genetics and evolutionary history of Plasmodium vivax

Plasmodium vivax is part of a highly diverse clade that includes several Plasmodium species found in nonhuman primates from Southeast Asia. The diversity of primate malarias in Asia is staggering; nevertheless, their origin was relatively recent in the evolution of Plasmodium. We discuss how humans acquired the lineage leading to P. vivax from a nonhuman primate determined by the complex geological processes that took place in Southeast Asia during the last few million years. We conclude that widespread population genomic investigations are needed in order to understand the demographic processes involved in the expansion of P. vivax in the human populations. India represents one of the few countries with widespread vivax malaria. Earlier studies have indicated high genetic polymorphism at antigenic loci and no evidence for geographic structuring. However, new studies using genetic markers in selectively neutral genetic regions indicate that Indian P. vivax presents complex evolutionary history but possesses features consistent with being part of the ancestral distribution range of this species. Such studies are possible due to the availability of the first P. vivax genome sequences. Next generation sequencing technologies are now paving the way for the sequencing of more P. vivax genomes that will dramatically increase our understanding of the unique biology of this species.

Genetic diversity and signatures of selection of drug resistance in Plasmodium populations from both human and mosquito hosts in continental Equatorial Guinea

Background

In Plasmodium, the high level of genetic diversity and the interactions established by co-infecting parasite populations within the same host may be a source of selection on pathogen virulence and drug resistance. As different patterns have already been described in humans and mosquitoes, parasite diversity and population structure should be studied in both hosts to properly assess their effects on infection and transmission dynamics. This study aimed to characterize the circulating populations of Plasmodium spp and Plasmodium falciparum from a combined set of human blood and mosquito samples gathered in mainland Equatorial Guinea. Further, the origin and evolution of anti-malarial resistance in this area, where malaria remains a major public health problem were traced.

Methods

Plasmodium species infecting humans and mosquitoes were identified by nested-PCR of chelex-extracted DNA from dried blood spot samples and mosquitoes. Analysis of Pfmsp2 gene, anti-malarial-resistance associated genes, Pfdhps, Pfdhfr, Pfcrt and Pfmdr1, neutral microsatellites (STR) loci and Pfdhfr and Pfdhps flanking STR was undertaken to evaluate P. falciparum diversity.

Results

Prevalence of infection remains high in mainland Equatorial Guinea. No differences in parasite formula or significant genetic differentiation were seen in the parasite populations in both human and mosquito samples. Point mutations in all genes associated with anti-malarial resistance were highly prevalent. A high prevalence was observed for the Pfdhfr triple mutant in particular, associated with pyrimethamine resistance.

Analysis of Pfdhps and Pfdhfr flanking STR revealed a decrease in the genetic diversity. This finding along with multiple independent introductions of Pfdhps mutant haplotypes suggest a soft selective sweep and an increased differentiation at Pfdhfr flanking microsatellites hints a model of positive directional selection for this gene.

Conclusions

Chloroquine is no longer recommended for malaria treatment in Equatorial Guinea but sulphadoxine-pyrimethamine (SP) remains in use in combination with artesunate and is the only drug recommended in preventive chemotherapy in pregnancy. The high prevalence of point mutations in Pfdhfr and Pfdhps points to the danger of an eventual reduction in the efficacy of SP combined therapy in P. falciparum populations in Equatorial Guinea and to the essential continuous monitoring of these two genes.

Population genetic structure of the Plasmodium vivax circumsporozoite protein (Pvcsp) in Sri Lanka

Molecular methods elucidate evolutionary and ecological processes in parasites, where interaction between hosts and parasites enlighten the evolution of parasite lifestyles and host defenses. Population genetics of Plasmodium vivax parasites accurately describe transmission dynamics of the parasites and evaluation of malaria control measures. As a first generation vaccine candidate against malaria, the Circumsporozoite Protein (CSP) has demonstrated significant potential in P. falciparum. Extensive polymorphism hinders the development of a potent malaria vaccine. Hence, the genetic diversity of Pvcsp was investigated for the first time in 60 Sri Lankan clinical isolates by obtaining the nucleotide sequence of the central repeat (CR) domain and examining the polymorphism of the peptide repeat motifs (PRMs), the genetic diversity indices and phylogenetic relationships. PCR amplicons determined size polymorphism of 610, 700 and 710 bp in Pvcsp of Sri Lanka where all amino acid sequences obtained were of the VK210 variant, consisting variable repeats of 4 different PRMs. The two most abundant PRMs of the CR domain, GDRADGQPA and GDRAAGQPA consisted ~2-4 repeats, while GNRAAGQPA was unique to the island. Though, different nucleotide sequences termed repeat allotypes (RATs) were observed for each PRM, these were synonymous contributing to a less polymorphic CR domain. The genetic diversity of Pvcsp in Sri Lanka was due to the number of repetitive peptide repeat motifs, point mutations, and intragenic recombination. The 19 amino acid haplotypes defined were exclusive to Sri Lanka, whereas the 194 Pvcsp sequences of global isolates generated 57 more distinct a.a. haplotypes of the VK210 variant. Strikingly, the CR domain of both VK210 and VK247 variants was under purifying selection interpreting the scarcity of CSP non-synonymous polymorphisms. Insights to the distribution of RATs in the CR region with geographic clustering of the P. vivax VK210 variant were revealed. The cladogram reiterated this unique geographic clustering of local (VK210) and global isolates (VK210 and VK247), which was further validated by the elevated fixation index values of the VK210 variant.

Fitness components and natural selection: why are there different patterns on the emergence of drug resistance in Plasmodium falciparum and Plasmodium vivax?

BACKGROUND

Considering the distinct biological characteristics of Plasmodium species is crucial for control and elimination efforts, in particular when facing the spread of drug resistance. Whereas the evolutionary fitness of all malarial species could be approximated by the probability of being taken by a mosquito and then infecting a new host, the actual steps in the malaria life cycle leading to a successful transmission event show differences among Plasmodium species. These "steps" are called fitness components. Differences in terms of fitness components may affect how selection imposed by interventions, e.g. drug treatments, differentially acts on each Plasmodium species. Thus, a successful malaria control or elimination programme should understand how differences in fitness components among different malaria species could affect adaptive evolution (e.g. the emergence of drug resistance). In this investigation, the interactions between some fitness components and natural selection are explored.

METHODS

A population-genetic model is formulated that qualitatively explains how different fitness components (in particular gametocytogenesis and longevity of gametocytes) affect selection acting on merozoites during the erythrocytic cycle. By comparing Plasmodium falciparum and Plasmodium vivax, the interplay of parasitaemia and gametocytaemia dynamics in determining fitness is modelled under circumstances that allow contrasting solely the differences between these two parasites in terms of their fitness components.

RESULTS

By simulating fitness components, it is shown that selection acting on merozoites (e.g., on drug resistant mutations or malaria antigens) is more efficient in P. falciparum than in P. vivax. These results could explain, at least in part, why resistance against drugs, such as chloroquine (CQ) is highly prevalent in P. falciparum worldwide, while CQ is still a successful treatment for P. vivax despite its massive use. Furthermore, these analyses are used to explore the importance of understanding the dynamic of gametocytaemia to ascertain the spreading of drug resistance.

CONCLUSIONS

The strength of natural selection on mutations that express their advantage at the merozoite stage is different in P. vivax and P. falciparum. Species-specific differences in gametocytogenesis and longevity of gametocytes need to be accounted for when designing effective malaria control and elimination programmes. There is a need for reliable data on gametocytogenesis from field studies.

Analysis of immunogenicity of different protein groups from malaria parasite Plasmodium falciparum

It was observed that pressure of host immune system leads to diversifying selection (which can be measured in terms of pN/pS ratio). In this research we checked whether Plasmodium falciparum proteins containing experimentally evident epitopes from the IEDB database are subject to diversifying selection. We also investigated which life stage of this parasite and which proteins are subject to the strongest immune pressure. To answer these questions we used information about experimentally evident epitopes from P. falciparum, that interact with human immune system and sequences of different isolates of P. falciparum obtained from PlasmoDB. We confirmed the expectations that proteins containing IEDB epitopes are subject to stronger diversifying selection which is evidenced by higher pN/pS ratio. A stage characterized by the highest average pN/pS ratio is that of the sporozoite. The greatest fraction of putative antigens is also present at this stage. We also found that the sporozoite stage is particularly interesting for further analysis as it potentially contains the highest number of unidentified epitopes.

Genetic polymorphism and natural selection in the C-terminal 42 kDa region of merozoite surface protein-1 among Plasmodium vivax Korean isolates

Background

The carboxy-terminal 42 kDa region of Plasmodium vivax merozoite surface protein-1 (PvMSP-1₄₂) is a leading candidate antigen for blood stage vaccine development. However, this region has been observed to be highly polymorphic among filed isolates of P. vivax. Therefore it is important to analyse the existing diversity of this antigen in the field isolates of P. vivax. In this study, the genetic diversity and natural selection in PvMSP-1₄₂ among P. vivax Korean isolates were analysed.

Methods

A total of 149 P. vivax-infected blood samples collected from patients in Korea were used. The region flanking PvMSP-1₄₂ was amplified by PCR, cloned into Escherichia coli, and then sequenced. The polymorphic characteristic and natural selection of PvMSP-1₄₂ were analysed using the DNASTAR, MEGA4 and DnaSP programs.

Results

A total of 11 distinct haplotypes of PvMSP-1₄₂ with 40 amino acid changes, as compared to the reference Sal I sequence, were identified in the Korean P. vivax isolates. Most of the mutations were concentrated in the 33 kDa fragment (PvMSP-1₃₃), but a novel mutation was found in the 19 kDa fragment (PvMSP-1₁₉). PvMSP-1₄₂ of Korean isolates appeared to be under balancing selection. Recombination may also play a role in the resulting genetic diversity of PvMSP-1₄₂.

Conclusions

PvMSP-1₄₂ of Korean P. vivax isolates displayed allelic polymorphisms caused by mutation, recombination and balancing selection. These results will be useful for understanding the nature of the P. vivax population in Korea and for development of a PvMSP-1₄₂ based vaccine against P. vivax.

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.

Genetic diversity in the block 2 region of the merozoite surface protein-1 of Plasmodium falciparum in central India

Background

Malaria continues to be a significant health problem in India. Several of the intended Plasmodium falciparum vaccine candidate antigens are highly polymorphic. The genetic diversity of P. falciparum merozoite surface protein-1 (MSP-1) has been extensively studied from various parts of the world. However, limited data are available from India. The aim of the present study was a molecular characterization of block 2 region of MSP-1 gene from the tribal-dominated, forested region of Madhya Pradesh.

Methods

DNA sequencing analysis was carried out in 71 field isolates collected between July 2005 to November 2005 and in 98 field isolates collected from July 2009 to December 2009. Alleles identified by DNA sequencing were aligned with the strain 3D7 and polymorphism analysis was done by using Edit Sequence tool (DNASTAR).

Results

The malaria positivity was 26% in 2005, which rose to 29% in 2009 and P. falciparum prevalence was also increased from 72% in 2005 to 81% in 2009. The overall allelic prevalence was higher in K1 (51%) followed by MAD20 (28%) and RO33 (21%) in 2005 while in 2009, RO33 was highest (40%) followed by K1 (36%) and MAD20 (24%).

Conclusions

The present study reports extensive genetic variations and dynamic evolution of block 2 region of MSP-1 in central India. Characterization of antigenic diversity in vaccine candidate antigens are valuable for future vaccine trials as well as understanding the population dynamics of P. falciparum parasites in this area.

Genetic diversity and population structure of genes encoding vaccine candidate antigens of Plasmodium vivax

Background

A major concern in malaria vaccine development is genetic polymorphisms typically observed among Plasmodium isolates in different geographical areas across the world. Highly polymorphic regions have been observed in Plasmodium falciparum and Plasmodium vivax antigenic surface proteins such as Circumsporozoite protein (CSP), Duffy-binding protein (DBP), Merozoite surface protein-1 (MSP-1), Apical membrane antigen-1 (AMA-1) and Thrombospondin related anonymous protein (TRAP).

Methods

Genetic variability was assessed in important polymorphic regions of various vaccine candidate antigens in P. vivax among 106 isolates from the Amazon Region of Loreto, Peru. In addition, genetic diversity determined in Peruvian isolates was compared to population studies from various geographical locations worldwide.

Results

The structured diversity found in P. vivax populations did not show a geographic pattern and haplotypes from all gene candidates were distributed worldwide. In addition, evidence of balancing selection was found in polymorphic regions of the trap, dbp and ama-1 genes.

Conclusions

It is important to have a good representation of the haplotypes circulating worldwide when implementing a vaccine, regardless of the geographic region of deployment since selective pressure plays an important role in structuring antigen diversity.

Evidence of purifying selection on merozoite surface protein 8 (MSP8) and 10 (MSP10) in Plasmodium spp

Evidence for natural selection, positive or negative, on gene encoding antigens may indicate variation or functional constraints that are immunologically relevant. Most malaria surface antigens with high genetic diversity have been reported to be under positive-diversifying selection. However, antigens with limited genetic variation are usually ignored in terms of the role that natural selection may have in generating such patterns. We investigated orthologous genes encoding two merozoite proteins, MSP8 and MSP10, among several mammalian Plasmodium spp. These antigens, together with MSP1, are among the few MSPs that have two epidermal growth factor-like domains (EGF) at the C-terminal. Those EGF are relatively conserved (low levels of genetic polymorphism) and have been proposed to act as ligands during the invasion of RBCs. We use several evolutionary genetic methods to detect patterns consistent with natural selection acting on MSP8 and MSP10 orthologs in the human parasites Plasmodium falciparum and P. vivax, as well as closely related malarial species found in non-human primates (NHPs). Overall, these antigens have low polymorphism in the human parasites in comparison with the orthologs from other Plasmodium spp. We found that the MSP10 gene polymorphism in P. falciparum only harbor non-synonymous substitutions, a pattern consistent with a gene under positive selection. Evidence of purifying selection was found on the polymorphism observed in both orthologs from P. cynomolgi, a non-human primate parasite closely related to P. vivax, but it was not conclusive in the human parasite. Yet, using phylogenetic base approaches, we found evidence for purifying selection on both MSP8 and MSP10 in the lineage leading to P. vivax. Such antigens evolving under strong functional constraints could become valuable vaccine candidates. We discuss how comparative approaches could allow detecting patterns consistent with negative selection even when there is low polymorphism in the extant populations.

Application of genomics to field investigations of malaria by the international centers of excellence for malaria research

Success of the global research agenda toward eradication of malaria will depend on development of new tools, including drugs, vaccines, insecticides and diagnostics. Genomic information, now available for the malaria parasites, their mosquito vectors, and human host, can be leveraged to both develop these tools and monitor their effectiveness. Although knowledge of genomic sequences for the malaria parasites, Plasmodium falciparum and Plasmodium vivax, have helped advance our understanding of malaria biology, simply knowing this sequence information has not yielded a plethora of new interventions to reduce the burden of malaria. Here we review and provide specific examples of how genomic information has increased our knowledge of parasite biology, focusing on P. falciparum malaria. We then discuss how population genetics can be applied toward the epidemiological and transmission-related goals outlined by the International Centers of Excellence for Malaria Research groups recently established by the National Institutes of Health. Finally, we propose genomics is a research area that can promote coordination and collaboration between various ICEMR groups, and that working together as a community can significantly advance the value of this information toward reduction of the global malaria burden.