Molecular evolution and intragenic recombination of the merozoite surface protein MSP-3alpha from the malaria parasite Plasmodium vivax in Thailand

Exploring the genetic diversity pattern of PvEBP/DBP2: A promising candidate for an effective Plasmodium vivax vaccine

Malaria remains a public health challenge. Since many control strategies have proven ineffective in eradicating this disease, new strategies are required, among which the design of a multivalent vaccine stands out. However, the effectiveness of this strategy has been hindered, among other reasons, by the genetic diversity observed in parasite antigens. In Plasmodium vivax, the Erythrocyte Binding Protein (PvEBP, also known as DBP2) is an alternate ligand to Duffy Binding Protein (DBP); given its structural resemblance to DBP, EBP/DBP2 is proposed as a promising antigen for inclusion in vaccine design. However, the extent of genetic diversity within the locus encoding this protein has not been comprehensively assessed. Thus, this study aimed to characterize the genetic diversity of the locus encoding the P. vivax EBP/DBP2 protein and to determine the evolutionary mechanisms modulating this diversity. Several intrapopulation genetic variation parameters were estimated using 36 gene sequences of PvEBP/DBP2 from Colombian P. vivax clinical isolates and 186 sequences available in databases. The study then evaluated the worldwide genetic structure and the evolutionary forces that may influence the observed patterns of genetic variation. It was found that the PvEBP/DBP2 gene exhibits one of the lowest levels of genetic diversity compared to other vaccine-candidate antigens. Four major haplotypes were shared worldwide. Analysis of the protein's 3D structure and epitope prediction identified five regions with potential antigenic properties. The results suggest that the PvEBP/DBP2 protein possesses ideal characteristics to be considered when designing a multivalent effective antimalarial vaccine against P. vivax.

Two 20-Residue-Long Peptides Derived from Plasmodium vivax Merozoite Surface Protein 10 EGF-Like Domains Are Involved in Binding to Human Reticulocytes

Plasmodium parasites’ invasion of their target cells is a complex, multi-step process involving many protein-protein interactions. Little is known about how complex the interaction with target cells is in Plasmodium vivax and few surface molecules related to reticulocytes’ adhesion have been described to date. Natural selection, functional and structural analysis were carried out on the previously described vaccine candidate P. vivax merozoite surface protein 10 (PvMSP10) for evaluating its role during initial contact with target cells. It has been shown here that the recombinant carboxyl terminal region (rPvMSP10-C) bound to adult human reticulocytes but not to normocytes, as validated by two different protein-cell interaction assays. Particularly interesting was the fact that two 20-residue-long regions (³⁸⁸DKEECRCRANYMPDDSVDYF⁴⁰⁷ and ⁴¹⁵KDCSKENGNCDVNAECSIDK⁴³⁴) were able to inhibit rPvMSP10-C binding to reticulocytes and rosette formation using enriched target cells. These peptides were derived from PvMSP10 epidermal growth factor (EGF)-like domains (precisely, from a well-defined electrostatic zone) and consisted of regions having the potential of being B- or T-cell epitopes. These findings provide evidence, for the first time, about the fragments governing PvMSP10 binding to its target cells, thus highlighting the importance of studying them for inclusion in a P. vivax antimalarial vaccine.

Whole genome sequencing of Plasmodium vivax isolates reveals frequent sequence and structural polymorphisms in erythrocyte binding genes

Plasmodium vivax malaria is much less common in Africa than the rest of the world because the parasite relies primarily on the Duffy antigen/chemokine receptor (DARC) to invade human erythrocytes, and the majority of Africans are Duffy negative. Recently, there has been a dramatic increase in the reporting of P. vivax cases in Africa, with a high number of them being in Duffy negative individuals, potentially indicating P. vivax has evolved an alternative invasion mechanism that can overcome Duffy negativity. Here, we analyzed single nucleotide polymorphism (SNP) and copy number variation (CNV) in Whole Genome Sequence (WGS) data from 44 P. vivax samples isolated from symptomatic malaria patients in southwestern Ethiopia, where both Duffy positive and Duffy negative individuals are found. A total of 123,711 SNPs were detected, of which 22.7% were nonsynonymous and 77.3% were synonymous mutations. The largest number of SNPs were detected on chromosomes 9 (24,007 SNPs; 19.4% of total) and 10 (16,852 SNPs, 13.6% of total). There were particularly high levels of polymorphism in erythrocyte binding gene candidates including merozoite surface protein 1 (MSP1) and merozoite surface protein 3 (MSP3.5, MSP3.85 and MSP3.9). Two genes, MAEBL and MSP3.8 related to immunogenicity and erythrocyte binding function were detected with significant signals of positive selection. Variation in gene copy number was also concentrated in genes involved in host-parasite interactions, including the expansion of the Duffy binding protein gene (PvDBP) on chromosome 6 and MSP3.11 on chromosome 10. Based on the phylogeny constructed from the whole genome sequences, the expansion of these genes was an independent process among the P. vivax lineages in Ethiopia. We further inferred transmission patterns of P. vivax infections among study sites and showed various levels of gene flow at a small geographical scale. The genomic features of P. vivax provided baseline data for future comparison with those in Duffy-negative individuals and allowed us to develop a panel of informative Single Nucleotide Polymorphic markers diagnostic at a micro-geographical scale.

Plasmodium vivax is the most geographically widespread parasite species that causes malaria in humans. Although it occurs in Africa as a member of a mix of Plasmodium species, P. vivax is dominant in other parts of the world outside of Africa (e.g., Brazil). It was previously thought that most African populations were immune to P. vivax infections due to the absence of Duffy antigen chemokine receptor (DARC) gene expression required for erythrocyte invasion. However, several recent reports have indicated the emergence and potential spread of P. vivax across human populations in Africa. Compared to Southeast Asia and South America where P. vivax is highly endemic, data on polymorphisms in erythrocyte binding gene candidates of P. vivax from Africa is limited. Filling this knowlege gap is critical for identifying functional genes in erythrocyte invasion, biomarkers for tracking the P. vivax isolates from Africa, as well as potential gene targets for vaccine development. This paper examined the level of genetic polymorphisms in a panel of 43 potential erythrocyte binding protein genes based on whole genome sequences and described transmission patterns of P. vivax infections from different study sites in Ethiopia based on the genetic variants. Our analyses showed that chromosomes 9 and 10 of the P. vivax genomes isolated in Ethiopia had the most high-quality genetic polymorphisms. Among all erythrocyte binding protein gene candidates, the merozoite surface proteins 1 and merozoite surface protein 3 showed high levels of polymorphism. MAEBL and MSP3.8 related to immunogenicity and erythrocyte binding function were detected with significant signals of positive selection. The expansion of the Duffy binding protein and merozoite surface protein 3 gene copies was an independent process among the P. vivax lineages in Ethiopia. Various levels of gene flow were observed even at a smaller geographical scale. Our study provided baseline data for future comparison with P. vivax in Duffy negative individuals and help develop a panel of genetic markers that are informative at a micro-geographical scale.

Insights into the molecular diversity of Plasmodium vivax merozoite surface protein-3γ (pvmsp3γ), a polymorphic member in the msp3 multi-gene family

Plasmodium vivax merozoite surface protein 3 (PvMSP3) is encoded by a multi-gene family. Of these, PvMSP3α, PvMSP3β and PvMSP3γ, are considered to be vaccine targets. Despite comprehensive analyses of PvMSP3α and PvMSP3β, little is known about structural and sequence diversity in PvMSP3γ. Analysis of 118 complete pvmsp3γ sequences from diverse endemic areas of Thailand and 9 reported sequences has shown 86 distinct haplotypes. Based on variation in insert domains, pvmsp3γ can be classified into 3 types, i.e. Belem, Salvador I and NR520. Imperfect nucleotide repeats were found in six regions of the gene; none encoded tandem amino acid repeats. Predicted coiled-coil heptad repeats were abundant in the protein and displayed variation in length and location. Interspersed phase shifts occurred in the heptad arrays that may have an impact on protein structure. Polymorphism in pvmsp3γ seems to be generated by intragenic recombination and driven by natural selection. Most P. vivax isolates in Thailand exhibit population structure, suggesting limited gene flow across endemic areas. Phylogenetic analysis has suggested that insert domains could have been subsequently acquired during the evolution of pvmsp3γ. Sequence and structural diversity of PvMSP3γ may complicate vaccine design due to alteration in predicted immunogenic epitopes among variants.

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.

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 characterization of Plasmodium vivax in the Kyrgyz Republic

At the end of 2016, Kyrgyz Republic was certified by the World Health Organization as a malaria-free country, while only a decade ago this disease posed a serious health threat. The progress achieved by Kyrgyz Republic provides a unique example of tertian (Plasmodium vivax) malaria elimination. This success was based on an integrated approach, including measures for the treatment of infected people and disease prevention, vector control and the development of an effective national epidemiological surveillance system. Lower P. vivax msp-1, msp-3α, csp and dbp_II genes polymorphism was revealed in Kyrgyz Republic in compare with that in Tajikistan. Molecular characterization of the causative agent found that P. vivax populations in Kyrgyz Republic was comprised by several lineages, highly divergent in the south-western and genetically homogeneous in the northern regions of Kyrgyz Republic, d. Such profile in the northern regions was compatible with several recent introductions rather than a long-term endemic circulation of the parasite. A low level of genetic variability suggested that the parasitic systems of tertian malaria, were not adapted, which, along with other factors, largely determined the possibility of malaria elimination in northern Kyrgyz Republic. Other determinants included environmental, social, and epidemiological factors that limited the spread of malaria. South-western Kyrgyz Republic, a region with a high level of interstate migration, requires considerable attention to prevent the spread of malaria.

On the Evolution and Function of Plasmodium vivax Reticulocyte Binding Surface Antigen (pvrbsa)

The RBSA protein is encoded by a gene described in Plasmodium species having tropism for reticulocytes. Since this protein is antigenic in natural infections and can bind to target cells, it has been proposed as a potential candidate for an anti-Plasmodium vivax vaccine. However, genetic diversity (a challenge which must be overcome for ensuring fully effective vaccine design) has not been described at this locus. Likewise, the minimum regions mediating specific parasite-host interaction have not been determined. This is why the rbsa gene’s evolutionary history is being here described, as well as the P. vivax rbsa (pvrbsa) genetic diversity and the specific regions mediating parasite adhesion to reticulocytes. Unlike what has previously been reported, rbsa was also present in several parasite species belonging to the monkey-malaria clade; paralogs were also found in Plasmodium parasites invading reticulocytes. The pvrbsa locus had less diversity than other merozoite surface proteins where natural selection and recombination were the main evolutionary forces involved in causing the observed polymorphism. The N-terminal end (PvRBSA-A) was conserved and under functional constraint; consequently, it was expressed as recombinant protein for binding assays. This protein fragment bound to reticulocytes whilst the C-terminus, included in recombinant PvRBSA-B (which was not under functional constraint), did not. Interestingly, two PvRBSA-A-derived peptides were able to inhibit protein binding to reticulocytes. Specific conserved and functionally important peptides within PvRBSA-A could thus be considered when designing a fully-effective vaccine against P. vivax.

Genome-wide scans for the identification of Plasmodium vivax genes under positive selection

BACKGROUND

The current trend of Plasmodium vivax cases imported from Southeast Asia into China has sharply increased recently, especially from the China-Myanmar border (CMB) area. High recombination rates of P. vivax populations associated with varied transmission intensity might cause distinct local selective pressures. The information on the genetic variability of P. vivax in this area is scant. Hence, this study assessed the genetic diversity of P. vivax genome sequence in CMB area and aimed to provide information on the positive selection of new gene loci.

RESULTS

This study reports a genome-wide survey of P. vivax in CMB area, using blood samples from local patients to identify population-specific selective processes. The result showed that considerable genetic diversity and mean pair-wise divergence among the sequenced P. vivax isolates were higher in some important gene families. Using the standardized integrated haplotype score (|iHS|) for all SNPs in chromosomal regions with SNPs above the top 1% distribution, it was observed that the top score locus involved 356 genes and most of them are associated with red blood cell invasion and immune evasion. The XP-EHH test was also applied and some important genes associated with anti-malarial drug resistance were observed in high positive scores list. This result suggests that P. vivax in CMB area is facing more pressure to survive than any other region and this has led to the strong positive selection of genes that are associated with host-parasite interactions.

CONCLUSIONS

This study suggests that greater genetic diversity in P. vivax from CMB area and positive selection signals in invasion and drug resistance genes are consistent with the history of drug use during malaria elimination programme in CMB area. Furthermore, this result also demonstrates that haplotype-based detecting selection can assist the genome-wide methods to identify the determinants of P. vivax diversity.

Plasmodium vivax msp-3α polymorphisms: analysis in the Indian subcontinent

BACKGROUND

Plasmodium vivax is the most widely distributed human malaria parasite and accounts for approximately the same number of malaria cases as Plasmodium falciparum in India. Compared with P. falciparum, P. vivax is difficult to eradicate because of its tendency to cause relapses, which impacts treatment and control strategies. The genetic diversity of these parasites, particularly of the merozoite surface protein-3 alpha (msp-3α) gene, can be used to help develop a potential vaccine. The present study aimed to investigate the genetic diversity of P. vivax using the highly polymorphic antigen gene msp-3α and to assess the suitability of using this gene for population genetic studies of P. vivax isolates and was carried out in 2004-06. No recent study has been reported for MSP 3α in the recent decade in India. Limited reports are available on the genetic diversity of the P. vivax population in India; hence, this report aimed to improve the understanding of the molecular epidemiology of the parasite by studying the P. vivax msp-3α (Pvmsp-3α) marker from P. vivax field isolates from India.

METHODS

Field isolates were collected from different sites distributed across eight states in India. A total of 182 blood samples were analysed by a nested polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) technique using the HhaI and AluI restriction enzymes to determine genetic msp-3α variation among clinical P. vivax isolates.

RESULTS

Based on the length variants of the PCR products of Pvmsp-3α gene, three allele sizes, Type A (1.8 kb), Type B (1.5 kb) and Type C (1.2 kb) were detected among the 182 samples. Type A PCR amplicon was more predominant (75.4 %) in the samples compared with the Type B (14.3 %) and Type C (10.0 %) polymorphisms. Among all of the samples analysed, 8.2 % were mixed infections detected by PCR alone. Restriction fragment length polymorphism (RFLP) analysis involving the restriction enzymes AluI and HhaI generated fragment sizes that were highly polymorphic and revealed substantial diversity at the nucleotide level.

CONCLUSIONS

The present study is the first extensive study in India using the Pvmsp-3α marker. The results indicated that Pvmps-3α, a polymorphic genetic marker of P. vivax, exhibited considerable variability in infection prevalence in field isolates from India. Additionally, the mean multiplicity of infection observed at all of the study sites indicated that P. vivax is highly diverse in nature in India, and Pvmsp-3α is likely an effective and promising epidemiological marker.

Dynamic changes of Plasmodium vivax population structure in South Korea

The vivax malaria epidemic has persisted in South Korea since its reemergence in 1993. Although there has been a significant decrease in the number of malaria cases in recent years, vivax malaria is still a major public health concern. To gain in-depth insight into the genetic makeup of Korean Plasmodium vivax, we analyzed polymorphic patterns of two major antigens, merozoite surface protein-1 (MSP-1) and MSP-3α, in 255 Korean P. vivax isolates collected over an extended period from 1998 to 2013. Combinational genetic analysis of polymorphic patterns of MSP-1 and MSP-3α in the isolates suggests that the P. vivax population in South Korea has been diversifying rapidly, with the appearance of parasites with new genotypes, despite the recent reduction of disease incidence. These results highlight the importance of molecular epidemiological investigations to supervise the genetic variation of the parasite in South Korea.

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.

Multiplicity of Infection and Disease Severity in Plasmodium vivax

Background

Multiplicity of infection (MOI) refers to the average number of distinct parasite genotypes concurrently infecting a patient. Although several studies have reported on MOI and the frequency of multiclonal infections in Plasmodium falciparum, there is limited data on Plasmodium vivax. Here, MOI and the frequency of multiclonal infections were studied in areas from South America where P. vivax and P. falciparum can be compared.

Methodology/Principal Findings

As part of a passive surveillance study, 1,328 positive malaria patients were recruited between 2011 and 2013 in low transmission areas from Colombia. Of those, there were only 38 P. vivax and 24 P. falciparum clinically complicated cases scattered throughout the time of the study. Samples from uncomplicated cases were matched in time and location with the complicated cases in order to compare the circulating genotypes for these two categories. A total of 92 P. vivax and 57 P. falciparum uncomplicated cases were randomly subsampled. All samples were genotyped by using neutral microsatellites. Plasmodium vivax showed more multiclonal infections (47.7%) than P. falciparum (14.8%). Population genetics and haplotype network analyses did not detect differences in the circulating genotypes between complicated and uncomplicated cases in each parasite. However, a Fisher exact test yielded a significant association between having multiclonal P. vivax infections and complicated malaria. No association was found for P. falciparum infections.

Conclusion

The association between multiclonal infections and disease severity in P. vivax is consistent with previous observations made in rodent malaria. The contrasting pattern between P. vivax and P. falciparum could be explained, at least in part, by the fact that P. vivax infections have lineages that were more distantly related among them than in the case of the P. falciparum multiclonal infections. Future research should address the possible role that acquired immunity and exposure may have on multiclonal infections and their association with disease severity.

Previous studies on rodent malarias and mathematical models have postulated a link between multiclonal infections and disease severity. This association has been tested in Plasmodium falciparum mostly in Africa with limited information on P. vivax. Furthermore, there is a paucity of information from areas with low transmission. Here, we used samples available from a passive surveillance carried out in Colombia, South America. We found an association between multiclonal infections and disease severity in P. vivax but not in P. falciparum. Although the number of complicated malaria cases is low, the contrasting pattern between these two species emphasizes their epidemiological differences. We discuss how this pattern could be the result of a higher divergence among the P. vivax lineages co-infecting a patient. We hypothesize that low levels of acquired immunity may play a role in the association between multiclonal infections and disease severity.

Assessment of an automated capillary system for Plasmodium vivax microsatellite genotyping

BACKGROUND

Several platforms have been used to generate the primary data for microsatellite analysis of malaria parasite genotypes. Each has relative advantages but share a limitation of being time- and cost-intensive. A commercially available automated capillary gel cartridge system was assessed in the microsatellite analysis of Plasmodium vivax diversity in the Peruvian Amazon.

METHODS

The reproducibility and accuracy of a commercially-available automated capillary system, QIAxcel, was assessed using a sequenced PCR product of 227 base pairs. This product was measured 42 times, then 27 P. vivax samples from Peruvian Amazon subjects were analyzed with this instrument using five informative microsatellites. Results from the QIAxcel system were compared with a Sanger-type sequencing machine, the ABI PRISM(®) 3100 Genetic Analyzer.

RESULTS

Significant differences were seen between the sequenced amplicons and the results from the QIAxcel instrument. Different runs, plates and cartridges yielded significantly different results. Additionally, allele size decreased with each run by 0.045, or 1 bp, every three plates. QIAxcel and ABI PRISM systems differed in giving different values than those obtained by ABI PRISM, and too many (i.e. inaccurate) alleles per locus were also seen with the automated instrument.

CONCLUSIONS

While P. vivax diversity could generally be estimated using an automated capillary gel cartridge system, the data demonstrate that this system is not sufficiently precise for reliably identifying parasite strains via microsatellite analysis. This conclusion reached after systematic analysis was due both to inadequate precision and poor reproducibility in measuring PCR product size.

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.

Malaria Molecular Epidemiology: Lessons from the International Centers of Excellence for Malaria Research Network

Molecular epidemiology leverages genetic information to study the risk factors that affect the frequency and distribution of malaria cases. This article describes molecular epidemiologic investigations currently being carried out by the International Centers of Excellence for Malaria Research (ICEMR) network in a variety of malaria-endemic settings. First, we discuss various novel approaches to understand malaria incidence and gametocytemia, focusing on Plasmodium falciparum and Plasmodium vivax. Second, we describe and compare different parasite genotyping methods commonly used in malaria epidemiology and population genetics. Finally, we discuss potential applications of molecular epidemiological tools and methods toward malaria control and elimination efforts.

Molecular genetic analysis of Plasmodium vivax isolates from Eastern and Central Sudan using pvcsp and pvmsp-3α genes as molecular markers

In Sudan, Plasmodium vivax accounts for approximately 5-10% of malaria cases. This study was carried out to determine the genetic diversity of P. vivax population from Sudan by analyzing the polymorphism of P. vivax csp (pvcsp) and pvmsp-3α genes. Blood samples (n=76) were taken from suspected malaria cases from 2012-2013 in three health centers of Eastern and Central Sudan. Parasite detection was performed by microscopy and molecular techniques, and genotyping of both genes was performed by PCR-RFLP followed by DNA sequence for only pvcsp gene (n=30). Based on microscopy analysis, 76 (%100) patients were infected with P. vivax, whereas nested-PCR results showed that 86.8% (n=66), 3.9% (n=3), and 3.9% (n=3) of tested samples had P. vivax as well as Plasmodium falciparum mono- and mixed infections, respectively. Four out of 76 samples had no results in molecular diagnosis. All sequenced samples were found to be of VK210 (100%) genotype with six distinct amino acid haplotypes, and 210A (66.7%) was the most prevalent haplotype. The Sudanese isolates displayed variations in the peptide repeat motifs (PRMs) ranging from 17 to 19 with GDRADGQPA (PRM1), GDRAAGQPA (PRM2) and DDRAAGQPA (PRM3). Also, 54 polymorphic sites with 56 mutations were found in repeat and post-repeat regions of the pvcsp and the overall nucleotide diversity (π) was 0.02149±0.00539. A negative value of dN-dS (-0.0344) was found that suggested a significant purifying selection of Sudanese pvcsp, (Z test, P<0.05). Regarding pvmsp-3α, three types were detected: types A (94.6%, 52/55), type C (3.6%, 2/55), and type B (1.8%, 1/55). No multiclonal infections were detected, and RFLP analysis identified 13 (Hha I, A1-A11, B1, and C1) and 16 (Alu I, A1-A14, B1, and C1) distinct allelic forms. In conclusion, genetic investigation among Sudanese P. vivax isolates indicated that this antigen showed limited antigenic diversity.

Heterogeneous genetic diversity pattern in Plasmodium vivax genes encoding merozoite surface proteins (MSP) -7E, -7F and -7L

Background

The msp-7 gene has become differentially expanded in the Plasmodium genus; Plasmodium vivax has the highest copy number of this gene, several of which encode antigenic proteins in merozoites.

Methods

DNA sequences from thirty-six Colombian clinical isolates from P. vivax (pv) msp-7E, −7F and -7L genes were analysed for characterizing and studying the genetic diversity of these pvmsp-7 members which are expressed during the intra-erythrocyte stage; natural selection signals producing the variation pattern so observed were evaluated.

Results

The pvmsp-7E gene was highly polymorphic compared to pvmsp-7F and pvmsp-7L which were seen to have limited genetic diversity; pvmsp-7E polymorphism was seen to have been maintained by different types of positive selection. Even though these copies seemed to be species-specific duplications, a search in the Plasmodium cynomolgi genome (P. vivax sister taxon) showed that both species shared the whole msp-7 repertoire. This led to exploring the long-term effect of natural selection by comparing the orthologous sequences which led to finding signatures for lineage-specific positive selection.

Conclusions

The results confirmed that the P. vivax msp-7 family has a heterogeneous genetic diversity pattern; some members are highly conserved whilst others are highly diverse. The results suggested that the 3′-end of these genes encode MSP-7 proteins’ functional region whilst the central region of pvmsp-7E has evolved rapidly. The lineage-specific positive selection signals found suggested that mutations occurring in msp-7s genes during host switch may have succeeded in adapting the ancestral P. vivax parasite population to humans.

Electronic supplementary material

The online version of this article (doi:10.1186/1475-2875-13-495) contains supplementary material, which is available to authorized users.

The Plasmodium vivax merozoite surface protein 3β sequence reveals contrasting parasite populations in southern and northwestern Thailand

BACKGROUND

Malaria control efforts have a significant impact on the epidemiology and parasite population dynamics. In countries aiming for malaria elimination, malaria transmission may be restricted to limited transmission hot spots, where parasite populations may be isolated from each other and experience different selection forces. Here we aim to examine the Plasmodium vivax population divergence in geographically isolated transmission zones in Thailand.

METHODOLOGY

We employed the P. vivax merozoite surface protein 3β (PvMSP3β) as a molecular marker for characterizing P. vivax populations based on the extensive diversity of this gene in Southeast Asian parasite populations. To examine two parasite populations with different transmission levels in Thailand, we obtained 45 P. vivax isolates from Tak Province, northwestern Thailand, where the annual parasite incidence (API) was more than 2%, and 28 isolates from Yala and Narathiwat Provinces, southern Thailand, where the API was less than 0.02%. We sequenced the PvMSP3β gene and examined its genetic diversity and molecular evolution between the parasite populations.

PRINCIPAL FINDINGS

Of 58 isolates containing single PvMSP3β alleles, 31 sequence types were identified. The overall haplotype diversity was 0.77 ± 0.06 and nucleotide diversity 0.0877±0.0054. The northwestern vivax malaria population exhibited extensive haplotype diversity (HD) of PvMSP3β (HD=1.0). In contrast, the southern parasite population displayed a single PvMSP3β allele (HD=0), suggesting a clonal population expansion. This result revealed that the extent of allelic diversity in P. vivax populations in Thailand varies among endemic areas.

CONCLUSION

Malaria parasite populations in a given region may vary significantly in genetic diversity, which may be the result of control and influenced by the magnitude of malaria transmission intensity. This is an issue that should be taken into account for the implementation of P. vivax control measures such as drug policy and vaccine development.

Low genetic diversity in the locus encoding the Plasmodium vivax P41 protein in Colombia's parasite population

Background

The development of malaria vaccine has been hindered by the allele-specific responses produced by some parasite antigens’ high genetic diversity. Such antigen genetic diversity must thus be evaluated when designing a completely effective vaccine. Plasmodium falciparum P12, P38 and P41 proteins have red blood cell binding regions in the s48/45 domains and are located on merozoite surface, P41 forming a heteroduplex with P12. These three genes have been identified in Plasmodium vivax and share similar characteristics with their orthologues in Plasmodium falciparum. Plasmodium vivax pv12 and pv38 have low genetic diversity but pv41 polymorphism has not been described.

Methods

The present study was aimed at evaluating the P. vivax p41 (pv41) gene’s polymorphism. DNA sequences from Colombian clinical isolates from pv41 gene were analysed for characterising and studying the genetic diversity and the evolutionary forces that produced the variation pattern so observed.

Results

Similarly to other members of the 6-Cys family, pv41 had low genetic polymorphism. pv41 3′-end displayed the highest nucleotide diversity value; several substitutions found there were under positive selection. Negatively selected codons at inter-species level were identified in the s48/45 domains; p41 would thus seem to have functional/structural constraints due to the presence of these domains.

Conclusions

In spite of the functional constraints of Pv41 s48/45 domains, immune system pressure seems to have allowed non-synonymous substitutions to become fixed within them as an adaptation mechanism; including Pv41 s48/45 domains in a vaccine should thus be carefully evaluated due to these domains containing some allele variants.

Electronic supplementary material

The online version of this article (doi:10.1186/1475-2875-13-388) contains supplementary material, which is available to authorized users.

Low genetic diversity and functional constraint in loci encoding Plasmodium vivax P12 and P38 proteins in the Colombian population

Background

Plasmodium vivax is one of the five species causing malaria in human beings, affecting around 391 million people annually. The development of an anti-malarial vaccine has been proposed as an alternative for controlling this disease. However, its development has been hampered by allele-specific responses produced by the high genetic diversity shown by some parasite antigens. Evaluating these antigens’ genetic diversity is thus essential when designing a completely effective vaccine.

Methods

The gene sequences of Plasmodium vivax p12 (pv12) and p38 (pv38), obtained from field isolates in Colombia, were used for evaluating haplotype polymorphism and distribution by population genetics analysis. The evolutionary forces generating the variation pattern so observed were also determined.

Results

Both pv12 and pv38 were shown to have low genetic diversity. The neutral model for pv12 could not be discarded, whilst polymorphism in pv38 was maintained by balanced selection restricted to the gene’s 5′ region. Both encoded proteins seemed to have functional/structural constraints due to the presence of s48/45 domains, which were seen to be highly conserved.

Conclusions

Due to the role that malaria parasite P12 and P38 proteins seem to play during invasion in Plasmodium species, added to the Pv12 and Pv38 antigenic characteristics and the low genetic diversity observed, these proteins might be good candidates to be evaluated in the design of a multistage/multi-antigen vaccine.

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.

Unrecognized fine-scale recombination can mimic the effects of adaptive radiation

Gene sequences can undergo accelerated nucleotide changes and rapid diversification. The rapid sequence changes can then potentially lead to phylogenetic incongruence. Recently, Bodilis et al. (2011) observed artificial phylogenetic incongruence using the Pseudomonas surface protein gene oprF, and hypothesized that it was the result of a long-branch attraction artifact ultimately caused by adaptive radiation. In this study, an alternative hypothesis, namely fine-scale recombination, was tested on the same dataset. The results reveal that regions in oprF are of different evolutionary origins, and the mosaic gene structure resulted in confounding phylogenetic signals. These findings demonstrate that unrecognized fine-scale recombination can confound the phylogenetic interpretation and emphasize the limitation of using whole genes as the unit of phylogenetic analysis.

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 selection in three Plasmodium vivax merozoite surface protein 7 (Pvmsp-7) genes in a Colombian population

A completely effective vaccine for malaria (one of the major infectious diseases worldwide) is not yet available; different membrane proteins involved in parasite-host interactions have been proposed as candidates for designing it. It has been found that proteins encoded by the merozoite surface protein (msp)-7 multigene family are antibody targets in natural infection; the nucleotide diversity of three Pvmsp-7 genes was thus analyzed in a Colombian parasite population. By contrast with P. falciparum msp-7 loci and ancestral P. vivax msp-7 genes, specie-specific duplicates of the latter specie display high genetic variability, generated by single nucleotide polymorphisms, repeat regions, and recombination. At least three major allele types are present in Pvmsp-7C, Pvmsp-7H and Pvmsp-7I and positive selection seems to be operating on the central region of these msp-7 genes. Although this region has high genetic polymorphism, the C-terminus (Pfam domain ID: PF12948) is conserved and could be an important candidate when designing a subunit-based antimalarial vaccine.

Molecular markers and genetic diversity of Plasmodium vivax

Enhanced understanding of the transmission dynamics and population genetics for Plasmodium vivax is crucial in predicting the emergence and spread of novel parasite phenotypes with major public health implications, such as new relapsing patterns, drug resistance and increased virulence. Suitable molecular markers are required for these population genetic studies. Here, we focus on two groups of molecular markers that are commonly used to analyse natural populations of P. vivax. We use markers under selective pressure, for instance, antigen-coding polymorphic genes, and markers that are not under strong natural selection, such as most minisatellite and microsatellite loci. First, we review data obtained using genes encoding for P. vivax antigens: circumsporozoite protein, merozoite surface proteins 1 and 3α, apical membrane antigen 1 and Duffy binding antigen. We next address neutral or nearly neutral molecular markers, especially microsatellite loci, providing a complete list of markers that have already been used in P. vivax populations studies. We also analyse the microsatellite loci identified in the P. vivax genome project. Finally, we discuss some practical uses for P. vivax genotyping, for example, detecting multiple-clone infections and tracking the geographic origin of isolates.

B cell epitope mapping and characterization of naturally acquired antibodies to the Plasmodium vivax merozoite surface protein-3α (PvMSP-3α) in malaria exposed individuals from Brazilian Amazon

The Plasmodium vivax Merozoite Surface Protein-3α (PvMSP-3α) is considered as a potential vaccine candidate. However, the detailed investigations of the type of immune responses induced in naturally exposed populations are necessary. Therefore, we aim to characterize the naturally induced antibody to PvMSP-3α in 282 individuals with different levels of exposure to malaria infections residents in Brazilian Amazon. PvMSP3 specific antibodies (IgA, IgG and IgG subclass) to five recombinant proteins and the epitope mapping by Spot-synthesis technique to full-protein sequence of amino acids (15aa sequence with overlapping sequence of 9aa) were performed. Our results indicates that PvMSP3 is highly immunogenic in naturally exposed populations, where 78% of studied individuals present IgG immune response against the full-length recombinant protein (PVMSP3-FL) and IgG subclass profile was similar to all five recombinant proteins studied with a high predominance of IgG1 and IgG3. We also observe that IgG and subclass levels against PvMSP3 are associated with malaria exposure. The PvMSP3 epitope mapping by Spot-synthesis shows a natural recognition of at least 15 antigenic determinants, located mainly in the two blocks of repeats, confirming the high immunogenicity of this region. In conclusion, PvMSP-3α is immunogenic in naturally exposed individuals to malaria infections and that antibodies to PvMSP3 are induced to several B cell epitopes. The presence of PvMSP3 cytophilic antibodies (IgG1 and IgG3), suggests that this mechanism could also occur in P. vivax.

Plasmodium falciparum and Plasmodium vivax infections in the Peruvian Amazon: propagation of complex, multiple allele-type infections without super-infection

Outcrossing potential between Plasmodium parasites is defined by the population-level diversity (PLD) and complexity of infection (COI). There have been few studies of PLD and COI in low transmission regions. Since the 1995-1998 Peruvian Amazon epidemic, there has been sustained transmission with < 0.5 P. falciparum and < 1.6 P. vivax infections/person/year. Using weekly active case detection, we described PLD by heterozygosity (H(e)) and COI using P. falciparum Pfmsp1-B2 and P. vivax Pvmsp3alpha. Not being homologous genes, we limited comparisons to within species. P. falciparum (N = 293) had low (H(e) = 0.581) and P. vivax (N = 186) had high (H(e) = 0.845) PLD. A total of 9.5% P. falciparum infections and 26.3% P. vivax infections had COI > 1. Certain allele types were in more mixed infections than expected by chance. The few appearances of new alleles could be explained by stochastic polymerase chain reaction detection or synchronization/sequestration. The results suggest propagation of mixed infections by multiple inocula, not super-infection, implying decade-long opportunity for outcrossing in these mixed infections.

Limited genetic polymorphism of the Plasmodium vivax low molecular weight rhoptry protein complex in the Colombian population

Proteins involved in parasite adhesion and invasion are considered the best candidates for the development of asexual blood-stage antimalarial vaccines. Such vaccine candidates should be accessible by the immune system and have limited diversity. Considering the promising results obtained in previous trials by immunizing monkeys with the rhoptry-associated proteins 1 and 2 (RAP-1 and RAP-2), here we assessed the genetic variability of the Plasmodium vivax rap-1 and rap-2 genes isolated from Colombian parasite populations. Limited sequence diversity was found in these genes, possibly as a result of a functional/structural restriction. The presence of several haplotypes at relatively low frequencies and the excess of singleton mutations suggests that a demographic process might be affecting the loci. Our results support the inclusion of PvRAP-1 and PvRAP-2 in the design of an antimalarial subunit-based vaccine against P. vivax, which would avoid inducing allele-specific immunity.

Molecular characterization of Plasmodium vivax clinical isolates in Pakistan and Iran using pvmsp-1, pvmsp-3alpha and pvcsp genes as molecular markers

In this study, the diversity of Plasmodium vivax populations circulating in Pakistan and Iran has been investigated by using circumsporozoite protein (csp) and merozoite surface proteins 1 and 3alpha (msp-1 and msp-3alpha) genes as genetic markers. Infected P. vivax blood samples were collected from Pakistan (n=187) and Iran (n=150) during April to October 2008, and were analyzed using nested-PCR/RFLP and sequencing methods. Genotyping pvmsp-1 (variable block 5) revealed the presence of type 1, type 2 and recombinant type 3 allelic variants, with type 1 predominant, in both study areas. The sequence analysis of 33 P. vivax isolates from Pakistan and 30 from Iran identified 16 distinct alleles each, with one allele (R-8) from Iran which was not reported previously. Genotyping pvcsp gene also showed that VK210 type is predominant in both countries. Moreover, based on the size of amplified fragment of pvmsp-3alpha, three major types: type A (1800bp), type B (1500bp) and type C (1200bp), were distinguished among the examined isolates that type A was predominant among Pakistani (72.7%) and Iranian (77.3%) parasites. PCR/RFLP products of pvmsp-3alpha with HhaI and AluI have detected 40 and 39 distinct variants among Pakistani and Iranian examined isolates, respectively. Based on these three studied genes, the rate of combined multiple genotypes were 30% and 24.6% for Pakistani and Iranian P. vivax isolates, respectively. These results indicate an extensive diversity in the P. vivax populations in both studies.

Limited global diversity of the Plasmodium vivax merozoite surface protein 4 gene

Merozoite surface proteins (MSPs) of the malaria parasites are major candidates for vaccine development targeting asexual blood stages. However, the diverse antigenic repertoire of these antigens that induce strain-specific protective immunity in human is a major challenge for vaccine design and often determines the efficacy of a vaccine. Here we further assessed the genetic diversity of Plasmodium vivax MSP4 (PvMSP4) protein using 195 parasite samples collected mostly from Thailand, Indonesia and Brazil. Overall, PvMSP4 is highly conserved with only eight amino acid substitutions. The majority of the haplotype diversity was restricted to the two short tetrapeptide repeat arrays in exon 1 and 2, respectively. Selection and neutrality tests indicated that exon 1 and the entire coding region of PvMSP4 were under purifying selection. Despite the limited nucleotide polymorphism of PvMSP4, significant genetic differentiation among the three major parasite populations was detected. Moreover, microgeographical heterogeneity was also evident in the parasite populations from different endemic areas of Thailand.

Genetic diversity of msp3alpha and msp1_b5 markers of Plasmodium vivax in French Guiana

Background

Reliable molecular typing tools are required for a better understanding of the molecular epidemiology of Plasmodium vivax. The genes msp3a and msp1_block5 are highly polymorphic and have been used as markers in many P. vivax population studies. These markers were used to assess the genetic diversity of P. vivax strains from French Guiana (South America) and to develop a molecular typing protocol.

Methods

A total of 120 blood samples from 109 patients (including 10 patients suffered from more than one malaria episode, samples were collected during each episode) with P. vivax infection were genotyped. All samples were analysed by msp3a PCR-RFLP and msp1_b5 gene sequencing was performed on 57 samples. Genotyping protocol applied to distinguish between new infection or relapse from heterologus hypnozoites and treatment failure or relapse from homologus hypnozoites was based on analysing first msp3a by PCR-RFLP and secondly, only if the genotypes of the two samples are identical, on sequencing the msp1_b5 gene.

Results

msp3a alleles of three sizes were amplified by PCR: types A, B and C. Eleven different genotypes were identified among the 109 samples analysed by msp3a PCR-RFLP. In 13.8% of cases, a mixed genotype infection was observed. The sequence of msp1_b5 gene revealed 22 unique genotypes and 12.3% of cases with mixed infection. In the 57 samples analysed by both methods, 45 genotypes were found and 21% were mixed. Among ten patients with two or three malaria episodes, the protocol allowed to identify five new infections or relapses from heterologous hypnozoites and six treatment failures of relapses from homologous hypnozoites.

Conclusion

The study showed a high diversity of msp3a and msp1_b5 genetic markers among P. vivax strains in French Guiana with a low polyclonal infection rate. These results indicated that the P. vivax genotyping protocol presented has a good discrimination power and can be used in clinical drug trials or epidemiological studies.

Clonal diversity of a lizard malaria parasite, Plasmodium mexicanum, in its vertebrate host, the western fence lizard: role of variation in transmission intensity over time and space

Within the vertebrate host, infections of a malaria parasite (Plasmodium) could include a single genotype of cells (single-clone infections) or two to several genotypes (multiclone infections). Clonal diversity of infection plays an important role in the biology of the parasite, including its life history, virulence, and transmission. We determined the clonal diversity of Plasmodium mexicanum, a lizard malaria parasite at a study region in northern California, using variable microsatellite markers, the first such study for any malaria parasite of lizards or birds (the most common hosts for Plasmodium species). Multiclonal infections are common (50-88% of infections among samples), and measures of genetic diversity for the metapopulation (expected heterozygosity, number of alleles per locus, allele length variation, and effective population size) all indicated a substantial overall genetic diversity. Comparing years with high prevalence (1996-1998 = 25-32% lizards infected), and years with low prevalence (2001-2005 = 6-12%) found fewer alleles in samples taken from the low-prevalence years, but no reduction in overall diversity (H = 0.64-0.90 among loci). In most cases, rare alleles appeared to be lost as prevalence declined. For sites chronically experiencing low transmission intensity (prevalence approximately 1%), overall diversity was also high (H = 0.79-0.91), but there were fewer multiclonal infections. Theory predicts an apparent excess in expected heterozygosity follows a genetic bottleneck. Evidence for such a distortion in genetic diversity was observed after the drop in parasite prevalence under the infinite alleles mutation model but not for the stepwise mutation model. The results are similar to those reported for the human malaria parasite, Plasmodium falciparum, worldwide, and support the conclusion that malaria parasites maintain high genetic diversity in host populations despite the potential for loss in alleles during the transmission cycle or during periods/locations when transmission intensity is low.

Genetic structures of geographically distinct Plasmodium vivax populations assessed by PCR/RFLP analysis of the merozoite surface protein 3beta gene

The recent resurgence of Plasmodium vivax malaria requires close epidemiological surveillance and monitoring of the circulating parasite populations. In this study, we developed a combination of polymerase chain reaction and restriction fragment length polymorphism (PCR/RFLP) method to investigate the genetic diversity of the P. vivax merozoite surface protein 3beta (PvMSP3beta) gene among four Asian parasite populations representing both tropical and temperate strains with dramatic divergent relapse patterns (N = 143). Using P. vivax field isolates from symptomatic patients, we have validated the feasibility of this protocol in distinguishing parasite genotypes. We have shown that PCR alone could detect three major size polymorphisms of the PvMSP3beta gene, and restriction analysis detected a total of 12 alleles within these Asian samples. Samples from different geographical areas differed dramatically in their PvMSP3beta allele composition and frequency, indicating that complex, yet different parasite genotypes were circulating in different endemic areas. This protocol allowed easy detections of multiple infections, which reached 20.5% in the samples from Thailand. It is interesting to note that samples from one temperate site in China collected during a recent outbreak of the disease also showed a high level of genetic diversity with multiple infections accounting for 5.6% of the samples. When combined with the PvMSP3alpha locus, this method provides better capability in distinguishing P. vivax genotypes and detecting mixed genotype infections.

Genetic diversity of Plasmodium vivax in Kolkata, India

Background

Plasmodium vivax malaria accounts for approximately 60% of malaria cases in Kolkata, India. There has been limited information on the genotypic polymorphism of P. vivax in this malaria endemic area. Three highly polymorphic and single copy genes were selected for a study of genetic diversity in Kolkata strains.

Methods

Blood from 151 patients with P. vivax infection diagnosed in Kolkata between April 2003 and September 2004 was genotyped at three polymorphic loci: the P. vivax circumsporozoite protein (pvcs), the merozoite surface protein 1 (pvmsp1) and the merozoite surface protein 3-alpha (pvmsp3-alpha).

Results

Analysis of these three genetic markers revealed that P. vivax populations in Kolkata are highly diverse. A large number of distinguishable alleles were found from three genetic markers: 11 for pvcs, 35 for pvmsp1 and 37 for pvmsp3-alpha. These were, in general, randomly distributed amongst the isolates. Among the 151 isolates, 142 unique genotypes were detected the commonest genotype at a frequency of less than 2% (3/151). The overall rate of mixed genotype infections was 10.6%.

Conclusion

These results indicate that the P. vivax parasite population is highly diverse in Kolkata, despite the low level of transmission. The genotyping protocols used in this study may be useful for differentiating re-infection from relapse and recrudescence in studies assessing of malarial drug efficacy in vivax malaria.

High sequence diversity and evidence of balancing selection in the Pvmsp3alpha gene of Plasmodium vivax in the Venezuelan Amazon

The genetic diversity of a defined Plasmodium vivax population from the Venezuelan Amazon was evaluated by direct sequencing of the gene encoding the P. vivax merozoite surface protein-3alpha, Pvmsp3alpha. Three allele sizes (1.9, 1.4 and 1.1kb) were amplified from 58 isolates with frequencies of 59.3%, 21.9% and 18.8%, respectively. 27 Pvmsp3alpha nucleotide sequences were determined, with nine distinct haplotypes observed. The genetic diversity (h) at this single locus was 0.774. The P. vivax population in this region exhibits significant diversity in contrast to the genetically restricted diversity of the sympatric P. falciparum population. Despite sharing vector and human hosts, different control strategies may be required for these two species in this region. Substitution patterns in the conserved C-terminus of Pvmsp3alpha showed a significant departure from neutrality, suggesting these polymorphisms are being maintained by frequency-dependent selection as the result of an effective immune response from the host. Our findings support the use of Pvmsp3alpha genotyping as a tool for monitoring interventions aimed at control of P. vivax.

Gene discovery in Plasmodium vivax through sequencing of ESTs from mixed blood stages

Despite the significance of Plasmodium vivax as the most widespread human malaria parasite and a major public health problem, gene expression in this parasite is poorly understood. To accelerate gene discovery and facilitate the annotation phase of the P. vivax genome project, we have undertaken a transcriptome approach to study gene expression in the mixed blood stages of a P. vivax field isolate. Using a cDNA library constructed from purified blood stages, we have obtained single-pass sequences for approximately 21,500 expressed sequence tags (ESTs), the largest number of transcript tags obtained so far for this species. Cluster analysis revealed that the library is highly redundant, resulting in 5407 clusters. Clustered ESTs were searched against public protein databases for functional annotation, and more than one-third showed a significant match, the majority of these to Plasmodium falciparum proteins. The most abundant clusters were to genes encoding ribosomal proteins and proteins involved in metabolism, consistent with the predominance of trophozoites in the field isolate sample. In spite of the scarcity of other parasite stages in the field isolate, we could identify genes that are expressed in rings, schizonts and gametocytes. This study should facilitate our understanding of the gene expression in P. vivax asexual stages and provide valuable data for gene prediction and annotation of the P. vivax genome sequence.