Bottleneck effects on vaccine-candidate antigen diversity of malaria parasites in Thailand

Genetic polymorphism of merozoite surface protein 1 and merozoite surface protein 2 in the Vietnam Plasmodium falciparum population

BACKGROUND

Plasmodium falciparum merozoite surface proteins 1 (PfMSP1) and 2 (PfMSP2) are potential candidates for malaria vaccine development. However, the genetic diversity of these genes in the global P. falciparum population presents a significant challenge in developing an effective vaccine. Hence, understanding the genetic diversity and evolutionary trends in the global P. falciparum population is crucial.

METHODS

This study analyzed the genetic variations and evolutionary changes of pfmsp1 and pfmsp2 in P. falciparum isolates from the Central Highland and South-Central regions of Vietnam. DNASTAR and MEGA7 programs were utilized for analyses. The polymorphic nature of global pfmsp1 and pfmsp2 was also investigated.

RESULTS

A total of 337 sequences of pfmsp1 and 289 sequences of pfmsp2 were obtained. The pfmsp1 and pfmsp2 from Vietnam revealed a higher degree of genetic homogeneity compared to those from other malaria-endemic countries. Remarkably, the allele diversity patterns of Vietnam pfmsp1 and pfmsp2 differed significantly from those of neighboring countries in the Greater Mekong Subregion. Declines in allele diversity and polymorphic patterns of Vietnam pfmsp1 and pfmsp2 were observed.

CONCLUSIONS

The Vietnam P. falciparum population might be genetically isolated from the parasite populations in other neighboring GMS countries, likely due to geographical barriers and distinct evolutionary pressures. Furthermore, bottleneck effects or selective sweeps may have contributed to the genetic homogeneity of Vietnam pfmsp1 and pfmsp2.

Genetic Diversity and Natural Selection of Plasmodium vivax Merozoite Surface Protein 8 in Global Populations

Plasmodium vivax Merozoite Surface Protein 8 (PvMSP8) is a promising candidate target for the development of multi-component vaccines. Therefore, determining the genetic variation pattern of Pvmsp8 is essential in providing a reference for the rational design of the P. vivax malaria vaccines. This study delves into the genetic characteristics of the Pvmsp8 gene, specifically focusing on samples from the China-Myanmar border (CMB) region, and contrasts these findings with broader global patterns. The study uncovers that Pvmsp8 exhibits a notable level of conservation across different populations, with limited polymorphisms and relatively low nucleotide diversity (0.00023-0.00120). This conservation contrasts starkly with the high polymorphisms found in other P. vivax antigens such as Pvmsp1. A total of 25 haplotypes and 14 amino acid mutation sites were identified in the global populations, and all mutation sites were confined to non-functional regions. The study also notes that most CMB Pvmsp8 haplotypes are shared among Burmese, Cambodian, Thai, and Vietnamese populations, indicating less geographical variance, but differ notably from those found in Pacific island regions or the Panama. The findings underscore the importance of considering regional genetic diversity in P. vivax when developing targeted malaria vaccines. Non departure from neutral evolution were found by Tajima's D test, however, statistically significant differences were observed between the kn/ks rates. The study's findings are crucial in understanding the evolution and population structure of the Pvmsp8 gene, particularly during regional malaria elimination efforts. The highly conserved nature of Pvmsp8, combined with the lack of mutations in its functional domain, presents it as a promising candidate for developing a broad and effective P. vivax vaccine. This research thus lays a foundation for the rational development of multivalent malaria vaccines targeting this genetically stable antigen.

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

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

Structural organization and sequence diversity of the complete nucleotide sequence encoding the Plasmodium malariae merozoite surface protein-1

The merozoite surface protein-1 (MSP1) is a prime candidate for an asexual blood stage vaccine against malaria. However, polymorphism in this antigen could compromise the vaccine’s efficacy. Although the extent of sequence variation in MSP1 has been analyzed from various Plasmodium species, little is known about structural organization and diversity of this locus in Plasmodium malariae (PmMSP1). Herein, we have shown that PmMSP1 contained five conserved and four variable blocks based on analysis of the complete coding sequences. Variable blocks were characterized by short insertion and deletion variants (block II), polymorphic nonrepeat sequences (block IV), complex repeat structure with size variation (block VI) and degenerate octapeptide repeats (block VIII). Like other malarial MSP1s, evidences of intragenic recombination have been found in PmMSP1. The rate of nonsynonymous nucleotide substitutions significantly exceeded that of synonymous nucleotide substitutions in block IV, suggesting positive selection in this region. Codon-based analysis of deviation from neutrality has identified a codon under purifying selection located in close proximity to the homologous region of the 38 kDa/42 kDa cleavage site of P. falciparum MSP1. A number of predicted linear B-cell epitopes were identified across both conserved and variable blocks of the protein. However, polymorphism in repeat-containing blocks resulted in alteration of the predicted linear B-cell epitope scores across variants. Although a number of predicted HLA-class II-binding peptides were identified in PmMSP1, all variants of block IV seemed not to be recognized by common HLA-class II alleles among Thai population, suggesting that diversity in this positive selection region could probably affect host immune recognition. The data on structural diversity in PmMSP1 could be useful for further studies such as vaccine development and strain characterization of this neglected malaria parasite.

Spatiotemporal Changes in Plasmodium vivax msp142 Haplotypes in Southern Mexico: From the Control to the Pre-Elimination Phase

For 20 years, Plasmodium vivax has been the only prevalent malaria species in Mexico, and cases have declined significantly and continuously. Spatiotemporal genetic studies can be helpful for understanding parasite dynamics and developing strategies to weaken malaria transmission, thus facilitating the elimination of the parasite. The aim of the current contribution was to analyze P. vivax-infected blood samples from patients in southern Mexico during the control (1993–2007) and pre-elimination phases (2008–2011). Nucleotide and haplotype changes in the pvmsp142 fragment were evaluated over time. The majority of multiple genotype infections occurred in the 1990s, when the 198 single nucleotide sequences exhibited 57 segregating sites, 64 mutations, and 17 haplotypes. Nucleotide and genetic diversity parameters showed subtle fluctuations from across time, in contrast to the reduced haplotype diversity and the increase in the R2 index and Tajima’s D value from 2008 to 2011. The haplotype network consisted of four haplogroups, the geographical distribution of which varied slightly over time. Haplogroup-specific B-cell epitopes were predicted. Since only high-frequency and divergent haplotypes persisted, there was a contraction of the parasite population. Given that 84% of haplotypes were exclusive to Mesoamerica, P. vivax flow is likely circumscribed to this region, representing important information for parasite surveillance.

Diversity and natural selection of Merozoite surface Protein-1 in three species of human malaria parasites: Contribution from South-East Asian isolates

The present study aimed to examine the genetic diversity of human malaria parasites (i.e., P. falciparum, P. vivax and P. knowlesi) in Malaysia and southern Thailand targeting the 19-kDa C-terminal region of Merozoite Surface Protein-1 (MSP-1₁₉). This region is essential for the recognition and invasion of erythrocytes and it is considered one of the leading candidates for asexual blood stage vaccines. However, the genetic data of MSP-1₁₉ among human malaria parasites in Malaysia is limited and there is also a need to update the current sequence diversity of this gene region among the Thailand isolates. In this study, genomic DNA was extracted from 384 microscopy-positive blood samples collected from patients who attended the hospitals or clinics in Malaysia and malaria clinics in Thailand from the year 2008 to 2016. The MSP-1₁₉ was amplified using PCR followed by bidirectional sequencing. DNA sequences identified in the present study were subjected to Median-joining network analysis with sequences of MSP-1₁₉ obtained from GenBank. DNA sequence analysis revealed that PfMSP-1₁₉ of Malaysian and Thailand isolates was not genetically conserved as high number of haplotypes were detected and positive selection was prevalent in PfMSP-1₁₉, hence questioning its suitability to be used as a vaccine candidate. A novel haplotype (Q/TNG/L) was also detected in Thailand P. falciparum isolate. In contrast, PvMSP-1₁₉ was highly conserved, however for the first time, a non-synonymous substitution (A1657S) was reported among Malaysian isolates. As for PkMSP-1₁₉, the presence of purifying selection and low nucleotide diversity indicated that it might be a potential vaccine target for P. knowlesi.

Structural diversity, natural selection and intragenic recombination in the Plasmodium vivax merozoite surface protein 9 locus in Thailand

The merozoite surface protein 9 (MSP9) of malarial parasite forms co-ligand complex with the 19 kDa fragment of merozoite surface protein 1 (MSP1) prior to erythrocyte invasion. Interruption of this process could hamper subsequent asexual erythrocytic development of malaria parasites; therefore, these proteins are considered potential vaccine candidates. In Plasmodium vivax, MSP9 (PvMSP9) contains both conserved and polymorphic repetitive domains that were immunogenic upon natural malaria exposure and conferred protection in vaccination studies in animal models. To investigate the extent of sequence diversity at this locus, 104 P. vivax isolates from 4 major malaria endemic areas of Thailand were analyzed. Results revealed that pvmsp9 contained 3 repeat domains (R1-R3) flanked by conserved domains. Repeat domains exhibit extensive sequence and length variation, in which 14, 39 and 16 haplotypes for domains R1-R3, respectively, circulated in this country. Sequence diversity in pvmsp9 among P. vivax isolates from each endemic area displayed population structure. The extent of sequence diversity in pvmsp9 isolates from the provinces of Tak, Chanthaburi, Ubon Ratchathani and Prachuap Khiri Khan in northwestern, eastern, northeastern and southwestern areas, respectively, was almost comparable and was remarkably higher than that from Yala/Narathiwat population in southern Thailand. Evidence for intragenic recombination in this locus was observed within each P. vivax population except among isolates from Yala and Narathiwat. Synonymous nucleotide diversity significantly exceeded nonsynonymous nucleotide diversity in domains R2 and R3, indicating purifying selection. However, micro-scale signatures of positive and negative selections occurred in both conserved and repeat domains, implying two opposing forces, probably from functional or structural constraint and host immune pressure, could have influenced diversity at this locus. The immunodominant T and B cell epitopes so far identified were invariant or highly conserved across isolates. Further analysis of global isolates is warranted for vaccine design based on this protein.

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.

Multiple Novel Mutations in Plasmodium falciparum Chloroquine Resistance Transporter Gene during Implementation of Artemisinin Combination Therapy in Thailand

Mutations in the chloroquine resistance transporter gene of Plasmodium falciparum (Pfcrt) are associated with drug susceptibility status of chloroquine and other antimalarials that interfere with heme detoxification process including artemisinin. We aim to investigate whether an increase in duration of artemisinin combination therapy (ACT) in Thailand could affect mutations in Pfcrt. The complete coding sequences of Pfcrt and dihydrofolate reductase (Pfdhfr), and size polymorphisms of the merozoite surface proteins-1 and 2 (Pfmsp-1 and Pfmsp-2) of 189 P. falciparum isolates collected during 1991 and 2016 were analyzed. In total, 12 novel amino acid substitutions and 13 novel PfCRT haplotypes were identified. The most prevalent haplotype belonged to the Dd2 sequence and no wild type was found. A significant positive correlation between the frequency of Pfcrt mutants and the year of sample collection was observed during nationwide ACT implementation (r = 0.780; P = 0.038). The number of haplotypes and nucleotide diversity of isolates collected during 3-day ACT (2009-2016) significantly outnumbered those collected before this treatment regimen. Positive Darwinian selection occurred in the transmembrane domains only among isolates collected during 3-day ACT but not among those collected before this period. No remarkable change was observed in the molecular indices for other loci analyzed when similar comparisons were performed. An increase in the duration of artesunate in combination therapy in Thailand could exert selective pressure on the Pfcrt locus, resulting in emergence of novel variants. The impact of these novel haplotypes on antimalarial susceptibilities requires further study.

Polymorphism in merozoite surface protein-7E of Plasmodium vivax in Thailand: Natural selection related to protein secondary structure

Merozoite surface protein 7 (MSP-7) is a multigene family expressed during malaria blood-stage infection. MSP-7 forms complex with MSP-1 prior to merozoite egress from erythrocytes, and could affect merozoite invasion of erythrocytes. To characterize sequence variation in the orthologue in P. vivax (PvMSP-7), a gene member encoding PvMSP-7E was analyzed among 92 Thai isolates collected from 3 major endemic areas of Thailand (Northwest: Tak, Northeast: Ubon Ratchathani, and South: Yala and Narathiwat provinces). In total, 52 distinct haplotypes were found to circulate in these areas. Although population structure based on this locus was observed between each endemic area, no genetic differentiation occurred between populations collected from different periods in the same endemic area, suggesting spatial but not temporal genetic variation. Sequence microheterogeneity in both N- and C- terminal regions was predicted to display 4 and 6 α-helical domains, respectively. Signals of purifying selection were observed in α-helices II-X, suggesting structural or functional constraint in these domains. By contrast, α-helix-I spanning the putative signal peptide was under positive selection, in which amino acid substitutions could alter predicted CD4+ T helper cell epitopes. The central region of PvMSP-7E comprised the 5’-trimorphic and the 3’-dimorphic subregions. Positive selection was identified in the 3’ dimorphic subregion of the central domain. A consensus of intrinsically unstructured or disordered protein was predicted to encompass the entire central domain that contained a number of putative B cell epitopes and putative protein binding regions. Evidences of intragenic recombination were more common in the central region than the remainders of the gene. These results suggest that the extent of sequence variation, recombination events and selective pressures in the PvMSP-7E locus seem to be differentially affected by protein secondary structure.

Genetic Diversity of the Plasmodium falciparum Glutamate-Rich Protein R2 Region Before and Twelve Years after Introduction of Artemisinin Combination Therapies among Febrile Children in Nigeria

The genetic diversity of glutamate-rich protein (GLURP) R2 region in Plasmodium falciparum isolates collected before and 12 years after the introduction of artemisinin combination treatment of malaria in Osogbo, Osun State, Nigeria, was compared in this study. Blood samples were collected on filter paper in 2004 and 2015 from febrile children from ages 1-12 years. The R2 region of the GLURP gene was genotyped using nested polymerase chain reaction and by nucleotide sequencing. In all, 12 GLURP alleles were observed in a total of 199 samples collected in the two study years. The multiplicity of infection (MOI) marginally increased over the two study years; however, the differences were statistically insignificant (2004 samples MOI = 1.23 versus 2015 samples MOI = 1.47). Some alleles were stable in their prevalence, whereas two GLURP alleles, VIII and XI, showed considerable variability between both years. This variability was replicated when GLURP sequences from other regions were compared with ours. The expected heterozygosity (He) values (He = 0.87) were identical for the two groups. High variability in the rearrangement of the amino acid repeat units in the R2 region were observed, with the amino acid repeat sequence DKNEKGQHEIVEVEEILPE more prevalent in both years, compared with the two other repeat sequences observed in the study. The parasite population characterized in this study displayed extensive genetic diversity. The detailed genetic profile of the GLURP R2 region has the potential to help guide further epidemiological studies aimed toward the rational design of novel chemotherapies that are antagonistic toward malaria.

Analysis of Pvama1 genes from China-Myanmar border reveals little regional genetic differentiation of Plasmodium vivax populations

Background

With the premise of diminishing parasite genetic diversity following the reduction of malaria incidence, the analysis of polymorphic antigenic markers may provide important information about the impact of malaria control on local parasite populations. Here we evaluated the genetic diversity of Plasmodium vivax apical membrane antigen 1 (Pvama1) gene in a parasite population from the China-Myanmar border and compared it with global P. vivax populations.

Methods

We performed evolutionary analysis to examine the genetic diversity, natural selection, and population differentiation of 73 Pvama1 sequences acquired from the China-Myanmar border as well as 615 publically available Pvama1 sequences from seven global P. vivax populations.

Results

A total of 308 Pvama1 haplotypes were identified among the global P. vivax isolates. The overall nucleotide diversity of Pvama1 gene among the 73 China-Myanmar border parasite isolates was 0.008 with 41 haplotypes being identified (Hd = 0.958). Domain I (DI) harbored the majority (26/33) of the polymorphic sites. The McDonald Kreitman test showed a significant positive selection across the ectodomain and the DI of Pvama1. The fixation index (F_ST) estimation between the China-Myanmar border, Thailand (0.01) and Myanmar (0.10) showed only slight geographical genetic differentiation. Notably, the Sal-I haplotype was not detected in any of the analyzed global isolates, whereas the Belem strain was restricted to the Thai population. The detected mutations are mapped outside the overlapped region of the predicted B-cell epitopes and intrinsically unstructured/disordered regions.

Conclusions

This study revealed high levels of genetic diversity of Pvama1 in the P. vivax parasite population from the China-Myanmar border with DI displaying stronger diversifying selection than other domains. There were low levels of population subdivision among parasite populations from the Greater Mekong Subregion.

Electronic supplementary material

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

Selective sweep suggests transcriptional regulation may underlie Plasmodium vivax resilience to malaria control measures in Cambodia

Cambodia, in which both Plasmodium vivax and Plasmodium falciparum are endemic, has been the focus of numerous malaria-control interventions, resulting in a marked decline in overall malaria incidence. Despite this decline, the number of P vivax cases has actually increased. To understand better the factors underlying this resilience, we compared the genetic responses of the two species to recent selective pressures. We sequenced and studied the genomes of 70 P vivax and 80 P falciparum isolates collected between 2009 and 2013. We found that although P falciparum has undergone population fracturing, the coendemic P vivax population has grown undisrupted, resulting in a larger effective population size, no discernable population structure, and frequent multiclonal infections. Signatures of selection suggest recent, species-specific evolutionary differences. Particularly, in contrast to P falciparum, P vivax transcription factors, chromatin modifiers, and histone deacetylases have undergone strong directional selection, including a particularly strong selective sweep at an AP2 transcription factor. Together, our findings point to different population-level adaptive mechanisms used by P vivax and P falciparum parasites. Although population substructuring in P falciparum has resulted in clonal outgrowths of resistant parasites, P vivax may use a nuanced transcriptional regulatory approach to population maintenance, enabling it to preserve a larger, more diverse population better suited to facing selective threats. We conclude that transcriptional control may underlie P vivax's resilience to malaria control measures. Novel strategies to target such processes are likely required to eradicate P vivax and achieve malaria elimination.

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.

Natural selection of K13 mutants of Plasmodium falciparum in response to artemisinin combination therapies in Thailand

Resistance of Plasmodium falciparum to artemisinin combination therapy (ACT) in Southeast Asia can have a devastating impact on chemotherapy and control measures. In this study, the evolution of artemisinin-resistant P. falciparum in Thailand was assessed by exploring mutations in the K13 locus believed to confer drug resistance phenotype. P. falciparum-infected blood samples were obtained from patients in eight provinces of Thailand over two decades (1991-2014; n = 904). Analysis of the K13 gene was performed by either sequencing the complete coding region (n = 259) or mutation-specific PCR-restriction fragment length polymorphism method (n = 645). K13 mutations related to artesunate resistance were detected in isolates from Trat province bordering Cambodia in 1991, about 4 years preceding widespread deployment of ACT in Thailand and increased in frequency over time. Nonsynonymous nucleotide diversity exceeded synonymous nucleotide diversity in the propeller region of the K13 gene, supporting the hypothesis that this diversity was driven by natural selection. No single mutant appeared to be favoured in every population, and propeller-region mutants were rarely observed in linkage with each other in the same haplotype. On the other hand, there was a highly significant association between the occurrence of a propeller mutant and the insertion of two or three asparagines after residue 139 of K13. Whether this insertion plays a compensatory role for deleterious effects of propeller mutants on the function of the K13 protein requires further investigation. However, modification of duration of ACT from 2-day to 3-day regimens in 2008 throughout the country does not halt the increase in frequency of mutants conferring artemisinin resistance phenotype.

Plasmodium falciparum population structure in Sudan post artemisinin-based combination therapy

Over the past decade, Sudan has stepped up malaria control backed by WHO, and this has resulted in significant reduction in parasite rate, malaria morbidity and mortality. The present study analyzed Plasmodium falciparum parasites in four geographical separated areas, to examine whether the success in malaria control following the use of artemisinin-based combination therapy (ACT) has disrupted the population structure and evolution of the parasite. We examined 319 P. falciparum isolates collected between October 2009 and October 2012 in four different areas in Sudan (Jazira [central Sudan], Southern Darfur [western Sudan], Upper Nile [southern Sudan] and Kasala [eastern Sudan]). Twelve microsatellites were analyzed for allelic diversity, multi-locus haplotype and inter-population differentiation. Level of diversity was compared to that detected for three of the above microsatellites among P. falciparum parasites in central and eastern Sudan in 1999, prior to introduction of ACT. Diversity at each locus (unbiased heterozygosity [H]) was high in all areas (Jazira, H=0.67), (Southern Darfur, H=0.71), (Upper Nile, H=0.71), and (Kasala, H=0.63). Microsatellites were distributed widely and private alleles, detected in a single population, were rare. The extent of diversity in the above sites was similar to that seen, in 1999, in central (Khartoum, H=0.73) and eastern Sudan (Gedaref, H=0.75). Significant Linkage disequilibrium (LD) was observed between the microsatellites in all populations. Pairwise FST analysis revealed that parasites in the four areas could be considered as one population. However, the parasites in Sudan clustered away from parasites in West Africa and the Arabian Peninsula. Despite marked reduction in malaria risk in Sudan, the extent of diversity and parasite genetic structure are indicative of a large population size. Further considerable reduction in transmission would be needed before fragmented sub-population can be seen. In addition, the large divergence of P. falciparum in Sudan from West Africa and Arabian Peninsula populations may result from differential evolutionary pressures acting at the population level, which shall be considered in eradication plans.

Genetic polymorphism and amino acid sequence variation in Plasmodium falciparum GLURP R2 repeat region in Assam, India, at an interval of five years

Background

The R2 repeat region of GLURP has been reported as a good genetic marker for Plasmodium falciparum genotyping. Proper knowledge of the extent and nature of P. falciparum genetic diversity using highly immunogenic R2 repeat region in malaria-endemic areas is a crucial element to understand various aspects related to immunity acquisition and disease pathogenesis.

Methods

Population diversity of P. falciparum GLURP and amino acid sequence repeats in GLURP R2 region was studied in malaria-endemic Assam state, northeast India and compared at an interval of five years during 2005 (Group-A) and 2011 (Group-B).

Results

Of the 66 samples, a total of 55 samples showed positive PCR bands for GLURP R2 region and altogether ten types of alleles with size ranging from 501 bp to 1,050 bp (50 bp bin) were observed and coded as genotypes I-X. In Group-A (n = 29), 24 samples were found infected with single, four with double and one with triple P. falciparum genotype, while in Group-B (n = 26), single genotype was found in 23 samples, double in two samples and triple in one sample. Genotype IV showed significant increase (p = 0.002) during 2011 (Group-B). Genotypes I to V were more common in Group-B (62%), however genotypes VI to X were more frequently distributed in Group-A. The expected heterozygosity was found slightly higher in Group-A (H_E = 0.87) than Group-B (H_E = 0.85), whereas multiplicity of infection (MOI) in Group-A (MOI = 1.21 ± 0.49) and Group-B (MOI = 1.12 ± 0.43) did not display significant variation. The amino acid repeat sequence unit (AAU) DKNEKGQHEIVEVEEILPE (called ‘a’) was more frequent in the well-conserved part of R2 repeat region.

Conclusion

The present study is the first extensive study in India which has generated substantial data for understanding the type and distribution of naturally evolved genetic polymorphism at amino acid sequence level in GLURP R2 repeat region in P. falciparum. There was decrease in the PCR amplicon size as well as the number of AAU [amino acid repeat unit] in Group-B displaying the bottleneck effect. The present study described a new type of AAU ‘d’ which varied from the other previous known AAUs.

Electronic supplementary material

The online version of this article (doi:10.1186/1475-2875-13-450) 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.

Spatial variation in genetic diversity and natural selection on the thrombospondin-related adhesive protein locus of Plasmodium vivax (PvTRAP)

Thrombospondin-related adhesive protein (TRAP) of malaria parasites is essential for sporozoite motility and invasions into mosquito’s salivary gland and vertebrate’s hepatocyte; thereby, it is a promising target for pre-erythrocytic vaccine. TRAP of Plasmodium vivax (PvTRAP) exhibits sequence heterogeneity among isolates, an issue relevant to vaccine development. To gain insights into variation in the complete PvTRAP sequences of parasites in Thailand, 114 vivax malaria patients were recruited in 2006–2007 from 4 major endemic provinces bordering Myanmar (Tak in the northwest, n = 30 and Prachuap Khirikhan in the southwest, n = 25), Cambodia (Chanthaburi in the east, n = 29) and Malaysia (Yala and Narathiwat in the south, n = 30). In total, 26 amino acid substitutions were detected and 9 of which were novel, resulting in 44 distinct haplotypes. Haplotype and nucleotide diversities were lowest in southern P. vivax population while higher levels of diversities were observed in other populations. Evidences of positive selection on PvTRAP were demonstrated in domains II and IV and purifying selection in domains I, II and VI. Genetic differentiation was significant between each population except that between populations bordering Myanmar where transmigration was common. Regression analysis of pairwise linearized Fst and geographic distance suggests that P. vivax populations in Thailand have been isolated by distance. Sequence diversity of PvTRAP seems to be temporally stable over one decade in Tak province based on comparison of isolates collected in 1996 (n = 36) and 2006–2007. Besides natural selection, evidences of intragenic recombination have been supported in this study that could maintain and further generate diversity in this locus. It remains to be investigated whether amino acid substitutions in PvTRAP could influence host immune responses although several predicted variant T cell epitopes drastically altered the epitope scores. Knowledge on geographic diversity in PvTRAP constitutes an important basis for vaccine design provided that vaccination largely confers variant-specific immunity.

Genetic diversity of MSP1 Block 2 of Plasmodium vivax isolates from Manaus (central Brazilian Amazon)

The diversity of MSP1 in both Plasmodium falciparum and P. vivax is presumed be associated to parasite immune evasion. In this study, we assessed genetic diversity of the most variable domain of vaccine candidate N-terminal PvMSP1 (Block 2) in field isolates of Manaus. Forty-seven blood samples the polymorphism of PvMSP1 Block 2 generates four fragment sizes. In twenty-eight of them, sequencing indicated seven haplotypes of PvMSP1 Block 2 circulating among field isolates. Evidence of striking exchanges was observed with two stretches flanking the repeat region and two predicted recombination sites were described. Single nucleotide polymorphisms determined with concurrent infections per patient indicated that nonsynonymous substitutions occurred preferentially in the repeat-rich regions which also were predicted as B-cell epitopes. The comprehensive understanding of the genetic diversity of the promising Block 2 associated with clinical immunity and a reduced risk of infection by Plasmodium vivax would be important for the rationale of malaria vaccine designs.

Low level of sequence diversity at merozoite surface protein-1 locus of Plasmodium ovale curtisi and P. ovale wallikeri from Thai isolates

Background

The merozoite surface protein-1 (MSP-1) is a candidate target for the development of blood stage vaccines against malaria. Polymorphism in MSP-1 can be useful as a genetic marker for strain differentiation in malarial parasites. Although sequence diversity in the MSP-1 locus has been extensively analyzed in field isolates of Plasmodium falciparum and P. vivax, the extent of variation in its homologues in P. ovale curtisi and P. ovale wallikeri, remains unknown.

Methodology/Principal Findings

Analysis of the mitochondrial cytochrome b sequences of 10 P. ovale isolates from symptomatic malaria patients from diverse endemic areas of Thailand revealed co-existence of P. ovale curtisi (n = 5) and P. ovale wallikeri (n = 5). Direct sequencing of the PCR-amplified products encompassing the entire coding region of MSP-1 of P. ovale curtisi (PocMSP-1) and P. ovale wallikeri (PowMSP-1) has identified 3 imperfect repeated segments in the former and one in the latter. Most amino acid differences between these proteins were located in the interspecies variable domains of malarial MSP-1. Synonymous nucleotide diversity (π_S) exceeded nonsynonymous nucleotide diversity (π_N) for both PocMSP-1 and PowMSP-1, albeit at a non-significant level. However, when MSP-1 of both these species was considered together, π_S was significantly greater than π_N (p<0.0001), suggesting that purifying selection has shaped diversity at this locus prior to speciation. Phylogenetic analysis based on conserved domains has placed PocMSP-1 and PowMSP-1 in a distinct bifurcating branch that probably diverged from each other around 4.5 million years ago.

Conclusion/Significance

The MSP-1 sequences support that P. ovale curtisi and P. ovale wallikeri are distinct species. Both species are sympatric in Thailand. The low level of sequence diversity in PocMSP-1 and PowMSP-1 among Thai isolates could stem from persistent low prevalence of these species, limiting the chance of outcrossing at this locus.

Genetic variation in the Plasmodium falciparum circumsporozoite protein in India and its relevance to RTS,S malaria vaccine

RTS,S is the most advanced malaria vaccine candidate, currently under phase-III clinical trials in Africa. This Plasmodium falciparum vaccine contains part of the central repeat region and the complete C-terminal T cell epitope region (Th2R and Th3R) of the circumsporozoite protein (CSP). Since naturally occurring polymorphisms at the vaccine candidate loci are critical determinants of the protective efficacy of the vaccines, it is imperative to investigate these polymorphisms in field isolates. In this study we have investigated the genetic diversity at the central repeat, C-terminal T cell epitope (Th2R and Th3R) and N-terminal T cell epitope regions of the CSP, in P. falciparum isolates from Madhya Pradesh state of India. These isolates were collected through a 5-year prospective study aimed to develop a well-characterized field-site for the future evaluation of malaria vaccine in India. Our results revealed that the central repeat (63 haplotypes, n = 161) and C-terminal Th2R/Th3R epitope (24 haplotypes, n = 179) regions were highly polymorphic, whereas N-terminal non-repeat region was less polymorphic (5 haplotypes, n = 161) in this population. We did not find any evidence of the role of positive natural selection in maintaining the genetic diversity at the Th2R/Th3R regions of CSP. Comparative analysis of the Th2R/Th3R sequences from this study to the global isolates (n = 1160) retrieved from the GenBank database revealed two important points. First, the majority of the sequences (∼61%, n = 179) from this study were identical to the Dd2/Indochina type, which is also the predominant Th2R/Th3R haplotype in Asia (∼59%, n = 974). Second, the Th2R/Th3R sequences in Asia, South America and Africa are geographically distinct with little allele sharing between continents. In conclusion, this study provides an insight on the existing polymorphisms in the CSP in a parasite population from India that could potentially influence the efficacy of RTS,S vaccine in this region.

Sahel, savana, riverine and urban malaria in West Africa: Similar control policies with different outcomes

The study sites for the West African ICEMR are in three countries (The Gambia, Senegal, Mali) and are located within 750 km of each other. In addition, the National Malaria Control Programmes of these countries have virtually identical policies: (1) Artemisinin Combination Therapies (ACTs) for the treatment of symptomatic Plasmodium falciparum infection, (2) Long-Lasting Insecticide-treated bed Nets (LLINs) to reduce the Entomololgic Inoculation Rate (EIR), and (3) sulfadoxine-pyrimethamine for the Intermittent Preventive Treatment of malaria during pregnancy (IPTp). However, the prevalence of P. falciparum malaria and the status of malaria control vary markedly across the four sites with differences in the duration of the transmission season (from 4-5 to 10-11 months), the intensity of transmission (with EIRs from unmeasurably low to 4-5 per person per month), multiplicity of infection (from a mean of 1.0 to means of 2-5) and the status of malaria control (from areas which have virtually no control to areas that are at the threshold of malaria elimination). The most important priority is the need to obtain comparable data on the population-based prevalence, incidence and transmission of malaria before new candidate interventions or combinations of interventions are introduced for malaria control.

Plasmodium vivax populations revisited: mitochondrial genomes of temperate strains in Asia suggest ancient population expansion

BACKGROUND

Plasmodium vivax is the most widely distributed human malaria parasite outside of Africa, and its range extends well into the temperate zones. Previous studies provided evidence for vivax population differentiation, but temperate vivax parasites were not well represented in these analyses. Here we address this deficit by using complete mitochondrial (mt) genome sequences to elucidate the broad genetic diversity and population structure of P. vivax from temperate regions in East and Southeast Asia.

RESULTS

From the complete mtDNA sequences of 99 clinical samples collected in China, Myanmar and Korea, a total of 30 different haplotypes were identified from 26 polymorphic sites. Significant differentiation between different East and Southeast Asian parasite populations was observed except for the comparison between populations from Korea and southern China. Haplotype patterns and structure diversity analysis showed coexistence of two different groups in East Asia, which were genetically related to the Southeast Asian population and Myanmar population, respectively. The demographic history of P. vivax, examined using neutrality tests and mismatch distribution analyses, revealed population expansion events across the entire P. vivax range and the Myanmar population. Bayesian skyline analysis further supported the occurrence of ancient P. vivax population expansion.

CONCLUSIONS

This study provided further resolution of the population structure and evolution of P. vivax, especially in temperate/warm-temperate endemic areas of Asia. The results revealed divergence of the P. vivax populations in temperate regions of China and Korea from other populations. Multiple analyses confirmed ancient population expansion of this parasite. The extensive genetic diversity of the P. vivax populations is consistent with phenotypic plasticity of the parasites, which has implications for malaria control.

Understanding the population genetics of Plasmodium vivax is essential for malaria control and elimination

Traditionally, infection with Plasmodium vivax was thought to be benign and self-limiting, however, recent evidence has demonstrated that infection with P. vivax can also result in severe illness and death. Research into P. vivax has been relatively neglected and much remains unknown regarding the biology, pathogenesis and epidemiology of this parasite. One of the fundamental factors governing transmission and immunity is parasite diversity. An understanding of parasite population genetic structure is necessary to understand the epidemiology, diversity, distribution and dynamics of natural P. vivax populations. In addition, studying the population structure of genes under immune selection also enables investigation of the dynamic interplay between transmission and immunity, which is crucial for vaccine development. A lack of knowledge regarding the transmission and spread of P. vivax has been particularly highlighted in areas where malaria control and elimination programmes have made progress in reducing the burden of Plasmodium falciparum, yet P. vivax remains as a substantial obstacle. With malaria elimination back on the global agenda, mapping of global and local P. vivax population structure is essential prior to establishing goals for elimination and the roll-out of interventions. A detailed knowledge of the spatial distribution, transmission and clinical burden of P. vivax is required to act as a benchmark against which control targets can be set and measured. This paper presents an overview of what is known and what is yet to be fully understood regarding P. vivax population genetics, as well as the importance and application of P. vivax population genetics studies.

Tracking origins and spread of sulfadoxine-resistant Plasmodium falciparum dhps alleles in Thailand

The emergence and spread of drug-resistant Plasmodium falciparum have been a major impediment for the control of malaria worldwide. Earlier studies have shown that similar to chloroquine (CQ) resistance, high levels of pyrimethamine resistance in P. falciparum originated independently 4 to 5 times globally, including one origin at the Thailand-Cambodia border. In this study we describe the origins and spread of sulfadoxine-resistance-conferring dihydropteroate synthase (dhps) alleles in Thailand. The dhps mutations and flanking microsatellite loci were genotyped for P. falciparum isolates collected from 11 Thai provinces along the Burma, Cambodia, and Malaysia borders. Results indicated that resistant dhps alleles were fixed in Thailand, predominantly being the SGEGA, AGEAA, and SGNGA triple mutants and the AGKAA double mutant (mutated codons are underlined). These alleles had different geographical distributions. The SGEGA alleles were found mostly at the Burma border, while the SGNGA alleles occurred mainly at the Cambodia border and nearby provinces. Microsatellite data suggested that there were two major genetic lineages of the triple mutants in Thailand, one common for SGEGA/SGNGA alleles and another one independent for AGEAA. Importantly, the newly reported SGNGA alleles possibly originated at the Thailand-Cambodia border. All parasites in the Yala province (Malaysia border) had AGKAA alleles with almost identical flanking microsatellites haplotypes. They were also identical at putatively neutral loci on chromosomes 2 and 3, suggesting a clonal nature of the parasite population in Yala. In summary, this study suggests multiple and independent origins of resistant dhps alleles in Thailand.