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

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.

Genetic Structure of Introduced Plasmodium vivax Malaria Isolates in Greece, 2015-2019

Greece has been malaria-free since 1974, after an intense malaria control program. However, as Greece hosts migrant populations from P. vivax malaria-endemic countries, there is a risk of introducing the disease to specific vulnerable and receptive areas of the country. Knowledge of the genetic diversity of P. vivax populations is essential for understanding the dynamics of malaria disease transmission in a given region. We used nine highly polymorphic markers to genotype 124 P. vivax-infected archived DNA samples from human blood specimens referred to the NMRL from all over Greece throughout 2015-2019. The genotypic variability of the samples studied was noted, as they comprised several unique haplotypes, indicative of the importation of a large number of different P. vivax strains in the country. However, only a few events of local transmission were recorded. Genotyping revealed and confirmed the same clusters as those identified through epidemiological investigation. In only one introduction event was the index case found. No sustained/ongoing malaria transmissions in/between the studied regions or during consecutive years or additional foci of local transmission were observed. Genotyping is an important component in assisting malaria surveillance, as it provides information concerning the patterns of introduction and the effectiveness of implemented malaria control and elimination measures.

Molecular tools are crucial for malaria elimination

The eradication of Plasmodium parasites, responsible for malaria, is a daunting global public health task. It requires a comprehensive approach that addresses symptomatic, asymptomatic, and submicroscopic cases. Overcoming this challenge relies on harnessing the power of molecular diagnostic tools, as traditional methods like microscopy and rapid diagnostic tests fall short in detecting low parasitaemia, contributing to the persistence of malaria transmission. By precisely identifying patients of all types and effectively characterizing malaria parasites, molecular tools may emerge as indispensable allies in the pursuit of malaria elimination. Furthermore, molecular tools can also provide valuable insights into parasite diversity, drug resistance patterns, and transmission dynamics, aiding in the implementation of targeted interventions and surveillance strategies. In this review, we explore the significance of molecular tools in the pursuit of malaria elimination, shedding light on their key contributions and potential impact on public health.

Genome-based tools for onchocerciasis elimination: utility of the mitochondrial genome for delineating Onchocerca volvulus transmission zones

National programs in Africa have expanded their objectives from control of onchocerciasis (river blindness) as a public health problem to elimination of parasite transmission, motivated by the reduction of Onchocerca volvulus infection prevalence in many African meso- and hyperendemic areas due to mass drug administration of ivermectin (MDAi). Given the large, contiguous hypo-, meso-, and hyperendemic areas, sustainable elimination of onchocerciasis in sub-Saharan Africa requires delineation of geographic boundaries for parasite transmission zones, so that programs can consider the risk of parasite re-introduction through vector or human migration from areas with ongoing transmission when making decisions to stop MDAi. We propose that transmission zone boundaries can be delineated by characterising the parasite genetic population structure within and between potential zones. We analysed whole mitochondrial genome sequences of 189 O. volvulus adults to determine the pattern of genetic similarity across three West African countries: Ghana, Mali, and Côte d'Ivoire. Population genetic structure indicates that parasites from villages near the Pru, Daka, and Black Volta rivers in central Ghana belong to one parasite population, indicating that the assumption that river basins constitute individual transmission zones is not supported by the data. Parasites from Mali and Côte d'Ivoire are genetically distinct from those from Ghana. This research provides the basis for developing tools for elimination programs to delineate transmission zones, to estimate the risk of parasite re-introduction via vector or human movement when intervention is stopped in one area while transmission is ongoing in others, to identify the origin of infections detected post-treatment cessation, and to investigate whether persisting prevalence despite ongoing interventions in one area is due to parasites imported from others.

Advocating for PCR-RFLP as molecular tool within malaria programs in low endemic areas and low resource settings

The road to malaria elimination for low- and middle-income countries is paved with obstacles, including the complexity and high costs of advanced molecular methods for genomic analysis. The usefulness of PCR-RFLP as less complex and affordable molecular surveillance tool in low-endemic malaria regions was assessed in a cross-sectional study conducted in Suriname, currently striving for malaria elimination, but plagued by recent P. vivax outbreaks. Molecular analysis of two highly polymorphic genes Pvmsp-1 F2 and Pvmsp-3α was performed for 49 samples, collected during October 2019 through September 2021 from four different regions with varying malaria transmission risks. RFLP-profiling revealed that outbreak samples from three indigenous villages, almost exclusively, harbored a single clonal type, matching the "Palumeu" lineage previously described in 2019, despite multiple relapses and drug pressure exerted by mass drug administration events, suggesting a limited P. vivax hypnozoite reservoir in Suriname. In contrast, isolates originating from Sophie, a mining area in neighboring French Guiana displayed a highly heterogeneous parasite population consistent with its endemic malaria status, demonstrating the differentiating capacity and thus the usefulness of PCR-RFLP for P. vivax genetic diversity studies. Outbreak reconstruction emphasized the impact of undetected human movement and relapses on reintroduction and resurgence of P. vivax malaria and PCR-RFLP monitoring of circulating parasites guided the roll-out of targeted interventions. PCR-RFLP seems a suitable molecular alternative in low-endemic areas with restricted resources for outbreak analysis, for monitoring the spread or containment of circulating strains and for identification of imported cases or potential foci.

Procedure for Handling and Storage of Onchocerca volvulus Microfilariae Obtained from Skin Snips for Downstream Genetic Work

WHO and endemic countries target elimination of transmission of Onchocerca volvulus, the parasite causing onchocerciasis. Population genetic analysis of O. volvulus may provide data to improve the evidence base for decisions on when, where, and for how long to deploy which interventions and post-intervention surveillance to achieve elimination. Development of necessary methods and tools requires parasites suitable for genetic analysis. Based on our experience with microfilariae obtained from different collaborators, we developed a microfilariae transfer procedure for large-scale studies in the Democratic Republic of Congo (DRC) comparing safety and efficacy of ivermectin, the mainstay of current onchocerciasis elimination strategies, and moxidectin, a new drug. This procedure is designed to increase the percentage of microfilariae in skin snips suitable for genetic analysis, improve assignment to metadata, and minimize time and materials needed by the researchers collecting the microfilariae. Among 664 microfilariae from South Sudan, 35.7% and 39.5% failed the mitochondrial and nuclear qPCR assay. Among the 576 microfilariae from DRC, 16.0% and 16.7% failed these assays, respectively. This difference may not only be related to the microfilariae transfer procedure but also to other factors, notably the ethanol concentration in the tubes in which microfilariae were stored (64% vs. ≥75%).

Distinct Allelic Diversity of Plasmodium vivax Merozoite Surface Protein 3-Alpha (PvMSP-3α) Gene in Thailand Using PCR-RFLP

Considering the importance of merozoite surface proteins (MSPs) as vaccine candidates, this study was conducted to investigate the polymorphism and genetic diversity of Plasmodium vivax merozoite surface protein 3-alpha (PvMSP-3α) in Thailand. To analyze genetic diversity, 118 blood samples containing P. vivax were collected from four malaria-endemic areas in western and southern Thailand. The DNA was extracted and amplified for the PvMSP-3α gene using nested PCR. The PCR products were genotyped by PCR-RFLP with Hha I and Alu I restriction enzymes. The combination patterns of Hha I and Alu I RFLP were used to identify allelic variants. Genetic evaluation and phylogenic analysis were performed on 13 sequences, including 10 sequences from our study and 3 sequences from GenBank. The results revealed three major types of PvMSP-3α, 91.5% allelic type A (∼1.8 kb), 5.1% allelic type B (∼1.5 kb), and 3.4% allelic type C (∼1.2 kb), were detected based on PCR product size with different frequencies. Among all PvMSP-3α, 19 allelic subtypes with Hha I RFLP patterns were distinguished and 6 allelic subtypes with Alu I RFLP patterns were identified. Of these samples, 73 (61%) and 42 (35.6%) samples were defined as monoallelic subtype infection by Hha I and Alu I PCR-RFLP, respectively, whereas 77 (65.3%) samples were determined to be mixed-allelic subtype infection by the combination patterns of Hha I and Alu I RFLP. These results strongly indicate that PvMSP-3α gene is highly polymorphic, particularly in blood samples collected from the Thai-Myanmar border area (the western part of Thailand). The combination patterns of Hha I and Alu I RFLP of the PvMSP-3α gene could be considered for use as molecular epidemiologic markers for genotyping P. vivax isolates in Thailand.

Retrospective analysis of Plasmodium vivax genomes from a pre-elimination China inland population in the 2010s

Introduction

In malaria-free countries, imported cases are challenging because interconnections with neighboring countries with higher transmission rates increase the risk of parasite reintroduction. Establishing a genetic database for rapidly identifying malaria importation or reintroduction is crucial in addressing these challenges. This study aimed to examine genomic epidemiology during the pre-elimination stage by retrospectively reporting whole-genome sequence variation of 10 Plasmodium vivax isolates from inland China.

Methods

The samples were collected during the last few inland outbreaks from 2011 to 2012 when China implemented a malaria control plan. After next-generation sequencing, we completed a genetic analysis of the population, explored the geographic specificity of the samples, and examined clustering of selection pressures. We also scanned genes for signals of positive selection.

Results

China's inland populations were highly structured compared to the surrounding area, with a single potential ancestor. Additionally, we identified genes under selection and evaluated the selection pressure on drug-resistance genes. In the inland population, positive selection was detected in some critical gene families, including sera, msp3, and vir. Meanwhile, we identified selection signatures in drug resistance, such as ugt, krs1, and crt, and noticed that the ratio of wild-type dhps and dhfr-ts increased after China banned sulfadoxine-pyrimethamine (SP) for decades.

Discussion

Our data provides an opportunity to investigate the molecular epidemiology of pre-elimination inland malaria populations, which exhibited lower selection pressure on invasion and immune evasion genes than neighbouring areas, but increased drug resistance in low transmission settings. Our results revealed that the inland population was severely fragmented with low relatedness among infections, despite a higher incidence of multiclonal infections, suggesting that superinfection or co-transmission events are rare in low-endemic circumstances. We identified selective signatures of resistance and found that the proportion of susceptible isolates fluctuated in response to the prohibition of specific drugs. This finding is consistent with the alterations in medication strategies during the malaria elimination campaign in inland China. Such findings could provide a genetic basis for future population studies, assessing changes in other pre-elimination countries.

Population genomics in neglected malaria parasites

Malaria elimination includes neglected human malaria parasites Plasmodium vivax, Plasmodium ovale spp., and Plasmodium malariae. Biological features such as association with low-density infection and the formation of hypnozoites responsible for relapse make their elimination challenging. Studies on these parasites rely primarily on clinical samples due to the lack of long-term culture techniques. With improved methods to enrich parasite DNA from clinical samples, whole-genome sequencing of the neglected malaria parasites has gained increasing popularity. Population genomics of more than 2200 P. vivax global isolates has improved our knowledge of parasite biology and host-parasite interactions, identified vaccine targets and potential drug resistance markers, and provided a new way to track parasite migration and introduction and monitor the evolutionary response of local populations to elimination efforts. Here, we review advances in population genomics for neglected malaria parasites, discuss how the rich genomic information is being used to understand parasite biology and epidemiology, and explore opportunities for the applications of malaria genomic data in malaria elimination practice.

THE USE OF ARCHIVED GIEMSA-STAINED BLOOD SMEARS AND RDT FOR PCR-BASED GENOTYPING OF Plasmodium vivax MEROZOITE SURFACE PROTEIN-1 IN CENTRAL KALIMANTAN PROVINCE, INDONESIA

Background:

Plasmodium vivax is transmitted most across the country of Indonesia. The country has set out a malaria elimination program by 2030. The information on genetic diversity of malarial parasites relates to malaria transmission in an endemic area may provide the information that can help the malaria control program to achieve the target. Therefore, the purpose of this study was to determine the genetic diversity of the Pvmsp-1 gene in Central Kalimantan Province.

Materials and Methods:

Samples were 140 of archived Giemsa-stained blood smear and rapid detection test. Samples were divided into the indigenous and migrant populations. After confirmation by single-step PCR, only P. vivax and mixed infection samples were amplified to nested PCR for genotyping of Pvmsp-1 allelic variation in segments F1, F2, and F3.

Results:

Genotyping of 23 PCR positive samples resulted in 13 genotypes. In segment F1, three allelic variants type A containing subtype A1 (1,050 bp), A2 (350 bp), A3 (150 bp), and type B (100 bp). In segment F2, mono genotypes were detected as variant type A (1,050 bp) and type B3 (150 bp), multiple genotypes were detected as type B containing subtype B1 (250 bp), B2 (200 bp), and B3 (150bp). In segment F3, three allelic variants generated from four mono genotypes were type A (350 bp), type B (300 bp), and two type C (250 bp).

Conclusion:

The low allelic variation of Pvmsp-1 gene may reflect the actual situation of the low malaria endemic status of the study sites.

Genetic Diversity and Natural Selection of Plasmodium vivax Duffy Binding Protein-II From China-Myanmar Border of Yunnan Province, China

Malaria incidence has declined dramatically over the past decade and China was certified malaria-free in 2021. However, the presence of malaria in border areas and the importation of cases of malaria parasites are major challenges for the consolidation of the achievements made by China. Plasmodium vivax Duffy binding protein (PvDBP) performs a significant role in erythrocyte invasion, and is considered a promising P. vivax vaccine. However, the highly polymorphic region of PvDBP (PvDBP-II) impedes the development of blood-stage vaccine against P. vivax. In this study, we investigated the genetic diversity and natural selection of PvDBP-II among 124 P. vivax isolates collected from the China-Myanmar border (CMB) in Yunnan Province, China, during 2009–2011. To compare genetic diversity, natural selection, and population structure with CMB isolates, 85 pvdbp-II sequences of eastern Myanmar isolates were obtained from GenBank. In addition, global sequences of pvdbp-II were retrieved from GenBank to establish genetic differentiation relationships and networks with the CMB isolates. In total, 22 single nucleotide polymorphisms reflected in 20 non-synonymous and two synonymous mutations were identified. The overall nucleotide diversity of PvDBP-II from the 124 CMB isolates was 0.0059 with 21 haplotypes identified (Hd = 0.91). The high ratio of non-synonymous to synonymous mutations suggests that PvDBP-II had evolved under positive selection. Population structure analysis of the CMB and eastern Myanmar isolates were optimally grouped into five sub-populations (K = 5). Polymorphisms of PvDBP-II display that CMB isolates were genetically diverse. Mutation, recombination, and positive selection promote polymorphism of PvDBP-II of P. vivax population. Although low-level genetic differentiation in eastern Myanmar was identified along with the more effective malaria control measures, the complexity of population structure in malaria parasites has maintained. In conclusion, findings from this study advance knowledge of the understanding of the dynamic of P. vivax population, which will contribute to guiding the rational design of a PvDBP-II based vaccine.

Complex infections in vivax malaria: the more you look, the more you find

The human malaria parasite Plasmodium vivax commonly causes complex multiclonal infections. Recently, Dia et al. have developed innovative methods for single-cell sequencing (SCS) of P. vivax infections by adapting an approach used previously for Plasmodium falciparum. Their studies provide fascinating new insights into P. vivax intrahost diversity and evolution.

Genetic Diversity of Plasmodium vivax Cysteine-Rich Protective Antigen (PvCyRPA) in Field Isolates from Five Different Areas of the Brazilian Amazon

The Plasmodium vivax Cysteine-Rich Protective Antigen (PvCyRPA) has an important role in erythrocyte invasion and has been considered a target for vivax malaria vaccine development. Nonetheless, its genetic diversity remains uncharted in Brazilian malaria-endemic areas. Therefore, we investigated the pvcyrpa genetic polymorphism in 98 field isolates from the Brazilian Amazon and its impact on the antigenicity of predicted B-cell epitopes. Genetic diversity parameters, population genetic analysis, neutrality test and the median-joining network were analyzed, and the potential amino acid polymorphism participation in B-cell epitopes was investigated. One synonymous and 26 non-synonymous substitutions defined fifty haplotypes. The nucleotide diversity and Tajima’s D values varied across the coding gene. The exon-1 sequence had greater diversity than those of exon-2. Concerning the prediction analysis, seven sequences were predicted as linear B cell epitopes, the majority contained in conformational epitopes. Moreover, important amino acid polymorphism was detected in regions predicted to contain residues participating in B-cell epitopes. Our data suggest that the pvcyrpa gene presents a moderate polymorphism in the studied isolates and such polymorphisms alter amino acid sequences contained in potential B cell epitopes, an important observation considering the antigen potentiality as a vaccine candidate to cover distinct P. vivax endemic areas worldwide.

Plasmodium vivax Genetic Diversity in Panama: Challenges for Malaria Elimination in Mesoamerica

Panama and all nations within the Mesoamerican region have committed to eliminate malaria within this decade. With more than 90% of the malaria cases in this region caused by Plasmodium vivax, an efficient national/regional elimination plan must include a comprehensive study of this parasite's genetic diversity. Here, we retrospectively analyzed P. vivax genetic diversity in autochthonous and imported field isolates collected in different endemic regions in Panama from 2007 to 2020, using highly polymorphic markers (csp, msp-1, and msp-3α). We did the analysis using molecular techniques that are cost-effective for malaria molecular surveillance within Mesoamerica. Thus, we used molecular analyses that are feasible for malaria molecular surveillance within the region, and that can provide useful information for policy and decision making about malaria elimination. We also evaluated if haplotypes established by combining the genotypes found in these genes were associated with relevant epidemiological variables and showed structure across the transmission foci that have been observed in Panama. Ten different haplotypes were identified, some of them strongly associated with geographical origin, age, and collection year. Phylogenetic analysis of csp (central repeat domain) revealed that both major variant types (vk210 and vk247) were circulating in Panama. Variant vk247 was restricted to the eastern endemic regions, while vk210 was predominant (77.3%) and widespread, displaying higher diversity (14 alleles) and geographically biased alleles. The regional implications of these molecular findings for the control of P. vivax malaria to achieve elimination across Mesoamerica are discussed.

Temporal Changes in the Genetic Diversity of Plasmodium vivax Merozoite Surface Protein-1 in Myanmar

Despite a significant decline in the incidence of malaria in Myanmar recently, malaria is still an important public health concern in the country. Although Plasmodium falciparum is associated with the highest incidence of malaria in Myanmar, the proportion of P. vivax cases has shown a gradual increase in recent years. The genetic diversity of P. vivax merozoite surface protein-1 block 5-6 (pvmsp-1 ICB 5-6) in the P. vivax population of Myanmar was analyzed to obtain a comprehensive insight into its genetic heterogeneity and evolutionary history. High levels of genetic diversity of pvmsp-1 ICB 5-6 were identified in the P. vivax isolates collected from Myanmar between 2013 and 2015. Thirty-nine distinct haplotypes of pvmsp-1 ICB 5-6 (13 for Sal I type, 20 for recombinant type, and 6 for Belem type) were found at the amino acid level. Comparative analyses of the genetic diversity of pvmsp-1 ICB 5-6 sequences in the recent (2013–2015) and the past (2004) P. vivax populations in Myanmar revealed genetic expansion of the pvmsp-1 ICB 5-6 in recent years, albeit with a declined incidence. The recent increase in the genetic heterogeneity of Myanmar pvmsp-1 ICB 5-6 is attributed to a combination of factors, including accumulated mutations and recombination. These results suggest that the size of the P. vivax population in Myanmar is sufficient to enable the generation and maintenance of genetic diversity, warranting continuous molecular surveillance of genetic variation in Myanmar P. vivax.

Population structure and diversity of Plasmodium falciparum in children with asymptomatic malaria living in different ecological zones of Ghana

BACKGROUND

Genetic diversity in Plasmodium falciparum populations can be used to describe the resilience and spatial distribution of the parasite in the midst of intensified intervention efforts. This study used microsatellite analysis to evaluate the genetic diversity and population dynamics of P. falciparum parasites circulating in three ecological zones of Ghana.

METHODS

A total of 1168 afebrile children aged between 3 to 13 years were recruited from five (5) Primary schools in 3 different ecological zones (Sahel (Tamale and Kumbungu), Forest (Konongo) and Coastal (Ada and Dodowa)) of Ghana. Asymptomatic malaria parasite carriage was determined using microscopy and PCR, whilst fragment analysis of 6 microsatellite loci was used to determine the diversity and population structure of P. falciparum parasites.

RESULTS

Out of the 1168 samples examined, 16.1 and 39.5% tested positive for P. falciparum by microscopy and nested PCR respectively. The genetic diversity of parasites in the 3 ecological zones was generally high, with an average heterozygosity (He) of 0.804, 0.787 and 0.608 the rainy (peak) season for the Sahel, Forest and Coastal zones respectively. The mean He for the dry (off-peak) season were 0.562, 0.693 and 0.610 for the Sahel, Forest and Coastal zones respectively. Parasites from the Forest zone were more closely related to those from the Sahel than from the Coastal zone, despite the Coastal zone being closer in physical distance to the Forest zone. The fixation indexes among study sites ranged from 0.049 to 0.112 during the rainy season and 0.112 to 0.348 during the dry season.

CONCLUSION

A large asymptomatic parasite reservoir was found in the school children during both rainy and dry seasons, especially those in the Forest and Sahel savannah zones where parasites were also found to be related compared to those from the Coastal zone. Further studies are recommended to understand why despite the roll out of several malaria interventions in Ghana, high transmission still persist.

Plasmodium vivax Cysteine-Rich Protective Antigen Polymorphism at Exon-1 Shows Recombination and Signatures of Balancing Selection

Plasmodium vivax Cysteine-Rich Protective Antigen (CyRPA) is a merozoite protein participating in the parasite invasion of human reticulocytes. During natural P. vivax infection, antibody responses against PvCyRPA have been detected. In children, low anti-CyRPA antibody titers correlated with clinical protection, which suggests this protein as a potential vaccine candidate. This work analyzed the genetic and amino acid diversity of pvcyrpa in Mexican and global parasites. Consensus coding sequences of pvcyrpa were obtained from seven isolates. Other sequences were extracted from a repository. Maximum likelihood phylogenetic trees, genetic diversity parameters, linkage disequilibrium (LD), and neutrality tests were analyzed, and the potential amino acid polymorphism participation in B-cell epitopes was investigated. In 22 sequences from Southern Mexico, two synonymous and 21 nonsynonymous mutations defined nine private haplotypes. These parasites had the highest LD-R² index and the lowest nucleotide diversity compared to isolates from South America or Asia. The nucleotide diversity and Tajima's D values varied across the coding gene. The exon-1 sequence had greater diversity and Rm values than those of exon-2. Exon-1 had significant positive values for Tajima's D, β-α values, and for the Z (HA: dN > dS) and MK tests. These patterns were similar for parasites of different origin. The polymorphic amino acid residues at PvCyRPA resembled the conformational B-cell peptides reported in PfCyRPA. Diversity at pvcyrpa exon-1 is caused by mutation and recombination. This seems to be maintained by balancing selection, likely due to selective immune pressure, all of which merit further study.

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

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

Molecular epidemiology of Plasmodium falciparum by multiplexed amplicon deep sequencing in Senegal

BACKGROUND

Molecular epidemiology can provide important information regarding the genetic diversity and transmission of Plasmodium falciparum, which can assist in designing and monitoring elimination efforts. However, malaria molecular epidemiology including understanding the genetic diversity of the parasite and performing molecular surveillance of transmission has been poorly documented in Senegal. Next Generation Sequencing (NGS) offers a practical, fast and high-throughput approach to understand malaria population genetics. This study aims to unravel the population structure of P. falciparum and to estimate the allelic diversity, multiplicity of infection (MOI), and evolutionary patterns of the malaria parasite using the NGS platform.

METHODS

Multiplex amplicon deep sequencing of merozoite surface protein 1 (PfMSP1) and merozoite surface protein 2 (PfMSP2) in fifty-three P. falciparum isolates from two epidemiologically different areas in the South and North of Senegal, was carried out.

RESULTS

A total of 76 Pfmsp1 and 116 Pfmsp2 clones were identified and 135 different alleles were found, 56 and 79 belonged to the pfmsp1 and pfmsp2 genes, respectively. K1 and IC3D7 allelic families were most predominant in both sites. The local haplotype diversity (Hd) and nucleotide diversity (π) were higher in the South than in the North for both genes. For pfmsp1, a high positive Tajima's D (TD) value was observed in the South (D = 2.0453) while negative TD value was recorded in the North (D = - 1.46045) and F-Statistic (Fst) was 0.19505. For pfmsp2, non-directional selection was found with a highly positive TD test in both areas and Fst was 0.02111. The mean MOI for both genes was 3.07 and 1.76 for the South and the North, respectively, with a statistically significant difference between areas (p = 0.001).

CONCLUSION

This study revealed a high genetic diversity of pfmsp1 and pfmsp2 genes and low genetic differentiation in P. falciparum population in Senegal. The MOI means were significantly different between the Southern and Northern areas. Findings also showed that multiplexed amplicon deep sequencing is a useful technique to investigate genetic diversity and molecular epidemiology of P. falciparum infections.

Balancing selection and high genetic diversity of Plasmodium vivax circumsporozoite central region in parasites from Brazilian Amazon and Rio de Janeiro Atlantic Forest

Circumsporozoite protein (CSP) is the primary pre-erythrocytic vaccine target in Plasmodium species. Knowledge about their genetic diversity can help predict vaccine efficacy and the spread of novel parasite variants. Thus, we investigated pvcsp gene polymorphisms in 219 isolates (136 from Brazilian Amazon [BA], 71 from Rio de Janeiro Atlantic Forest [AF], and 12 from non-Brazilian countries [NB]). Forty-eight polymorphic sites were detected, 46 in the central repeat region (CR), and two in the C-terminal region. Also, the CR presents InDels and a variable number of repeats. All samples correspond to the VK210 variant, and 24 VK210 subtypes based on CR. Nucleotide diversity (π = 0.0135) generated a significant number of haplotypes (168) with low genetic differentiation between the Brazilian regions (Fst = 0.208). The haplotype network revealed similar distances among the BA and AF regions. The linkage disequilibrium indicates that recombination does not seem to be acting in diversity, reinforcing natural selection's role in accelerating adaptive evolution. The high diversity (low Fst) and polymorphism frequencies could be indicators of balancing selection. Although malaria in BA and AF have distinct vector species and different host immune pressures, consistent genetic signature was found in two regions. The immunodominant B-cell epitope mapped in the CR varies from seven to 19 repeats. The CR T-cell epitope is conserved only in 39 samples. Concerning to C-terminal region, the Th2R epitope presented nonsynonymous SNP only in 6% of Brazilian samples, and the Th3R epitope remained conserved in all studied regions. We conclude that, although the uneven distribution of alleles may jeopardize the deployment of vaccines directed to a specific variable locus, a unique vaccine formulation could protect populations in all Brazilian regions.

SNP barcodes provide higher resolution than microsatellite markers to measure Plasmodium vivax population genetics

BACKGROUND

Genomic surveillance of malaria parasite populations has the potential to inform control strategies and to monitor the impact of interventions. Barcodes comprising large numbers of single nucleotide polymorphism (SNP) markers are accurate and efficient genotyping tools, however may need to be tailored to specific malaria transmission settings, since 'universal' barcodes can lack resolution at the local scale. A SNP barcode was developed that captures the diversity and structure of Plasmodium vivax populations of Papua New Guinea (PNG) for research and surveillance.

METHODS

Using 20 high-quality P. vivax genome sequences from PNG, a total of 178 evenly spaced neutral SNPs were selected for development of an amplicon sequencing assay combining a series of multiplex PCRs and sequencing on the Illumina MiSeq platform. For initial testing, 20 SNPs were amplified in a small number of mono- and polyclonal P. vivax infections. The full barcode was then validated by genotyping and population genetic analyses of 94 P. vivax isolates collected between 2012 and 2014 from four distinct catchment areas on the highly endemic north coast of PNG. Diversity and population structure determined from the SNP barcode data was then benchmarked against that of ten microsatellite markers used in previous population genetics studies.

RESULTS

From a total of 28,934,460 reads generated from the MiSeq Illumina run, 87% mapped to the PvSalI reference genome with deep coverage (median = 563, range 56-7586) per locus across genotyped samples. Of 178 SNPs assayed, 146 produced high-quality genotypes (minimum coverage = 56X) in more than 85% of P. vivax isolates. No amplification bias was introduced due to either polyclonal infection or whole genome amplification (WGA) of samples before genotyping. Compared to the microsatellite panels, the SNP barcode revealed greater variability in genetic diversity between populations and geographical population structure. The SNP barcode also enabled assignment of genotypes according to their geographic origins with a significant association between genetic distance and geographic distance at the sub-provincial level.

CONCLUSIONS

High-throughput SNP barcoding can be used to map variation of malaria transmission dynamics at sub-national resolution. The low cost per sample and genotyping strategy makes the transfer of this technology to field settings highly feasible.

Population genomics reveals the expansion of highly inbred Plasmodium vivax lineages in the main malaria hotspot of Brazil

BACKGROUND

Plasmodium vivax is a neglected human malaria parasite that causes significant morbidity in the Americas, the Middle East, Asia, and the Western Pacific. Population genomic approaches remain little explored to map local and regional transmission pathways of P. vivax across the main endemic sites in the Americas, where great progress has been made towards malaria elimination over the past decades.

METHODOLOGY/PRINCIPAL FINDINGS

We analyze 38 patient-derived P. vivax genome sequences from Mâncio Lima (ML)-the Amazonian malaria hotspot next to the Brazil-Peru border-and 24 sequences from two other sites in Acre State, Brazil, a country that contributes 23% of malaria cases in the Americas. We show that the P. vivax population of ML is genetically diverse (π = 4.7 × 10-4), with a high polymorphism particularly in genes encoding proteins putatively involved in red blood cell invasion. Paradoxically, however, parasites display strong genome-wide linkage disequilibrium, being fragmented into discrete lineages that are remarkably stable across time and space, with only occasional recombination between them. Using identity-by-descent approaches, we identified a large cluster of closely related sequences that comprises 16 of 38 genomes sampled in ML over 26 months. Importantly, we found significant ancestry sharing between parasites at a large geographic distance, consistent with substantial gene flow between regional P. vivax populations.

CONCLUSIONS/SIGNIFICANCE

We have characterized the sustained expansion of highly inbred P. vivax lineages in a malaria hotspot that can seed regional transmission. Potential source populations in hotspots represent a priority target for malaria elimination in the Amazon.

On the survival of 48 h Plasmodium vivax Aotus monkey-derived ex vivo cultures: the role of leucocytes filtration and chemically defined lipid concentrate media supplementation

Background

Filtration of leukocytes (WBCs) is a standard practice of malaria ex vivo cultures. To date, few studies have considered the effect of filtration or the lack thereof on the survival of Plasmodium vivax ex vivo cultures through one cycle of maturation. This study investigates the effect of WBC filtration and culture media supplementation on the survival of 48–72 h ex vivo cultures.

Methods

Using parasitaemia density, the study compares the survival of Plasmodipur^® filtered, filter-retained or washed ex vivo cultures, maintained with McCoy’s5A medium supplemented with 25% serum alone or 20% in combination with 5% chemically defined lipid concentrate (CDLC), and in washed ex vivo cultures plus GlutaMAX™, benchmarked against IMDM™ or AIM-V™ media; also, assessed the survival of ex vivo cultures co-cultivated with human red blood cells (hRBCs).

Results

After 48 h of incubation a statistically significant difference was detected in the survival proportions of filtered and the filter-retained ex vivo cultures supplemented with serum plus CDLC (p = 0.0255), but not with serum alone (p = 0.1646). To corroborate these finding, parasitaemias of washed ex vivo cultures maintained with McCoy’s5A complete medium were benchmarked against IMDM™ or AIM-V™ media; again, a statistically significant difference was detected in the cultures supplemented with CDLC and GlutaMAX™ (p = 0.03), but not when supplemented with either alone; revealing a pattern of McCoy’s5A medium supplementation for Aotus-derived P. vivax cultures as follows: serum < serum + GlutaMAX™ < serum + CDLC < serum + CDLC + GlutaMAX™; confirming a key role of CDLC in combination with GlutaMAX™ in the enhanced survival observed. Lastly, results showed that co-cultivation with malaria-naïve hRBCs improved the survival of ex vivo cultures.

Conclusions

This study demonstrates that WBC filtration is not essential for the survival of P. vivax ex vivo cultures. It also demonstrates that McCoy’s5A complete medium improves the survival of Aotus-derived P. vivax ex vivo cultures, with no significant difference in survival compared to IMDM and AIM-V media. Finally, the study demonstrates that co-cultivation with hRBCs enhances the survival of ex vivo cultures. These findings are expected to help optimize seeding material for long-term P. vivax in vitro culture.

Genetic diversity and genotype multiplicity of Plasmodium falciparum infection in patients with uncomplicated malaria in Chewaka district, Ethiopia

Background

Genetic diversity in Plasmodium falciparum poses a major threat to malaria control and elimination interventions. Characterization of the genetic diversity of P. falciparum strains can be used to assess intensity of parasite transmission and identify potential deficiencies in malaria control programmes, which provides vital information to evaluating malaria elimination efforts. This study investigated the P. falciparum genetic diversity and genotype multiplicity of infection in parasite isolates from cases with uncomplicated P. falciparum malaria in Southwest Ethiopia.

Methods

A total of 80 P. falciparum microscopy and qPCR positive blood samples were collected from study participants aged 6 months to 60 years, who visited the health facilities during study evaluating the efficacy of artemether-lumefantrine from September–December, 2017. Polymorphic regions of the msp-1 and msp-2 were genotyped by nested polymerase chain reactions (nPCR) followed by gel electrophoresis for fragment analysis.

Results

Of 80 qPCR-positive samples analysed for polymorphisms on msp-1 and msp-2 genes, the efficiency of msp-1 and msp-2 gene amplification reactions with family-specific primers were 95% and 98.8%, respectively. Allelic variation of 90% (72/80) for msp-1 and 86.2% (69/80) for msp-2 were observed. K1 was the predominant msp-1 allelic family detected in 20.8% (15/72) of the samples followed by MAD20 and RO33. Within msp-2, allelic family FC27 showed a higher frequency (26.1%) compared to IC/3D7 (15.9%). Ten different alleles were observed in msp-1 with 6 alleles for K1, 3 alleles for MAD20 and 1 allele for RO33. In msp-2, 19 individual alleles were detected with 10 alleles for FC27 and 9 alleles for 3D7. Eighty percent (80%) of isolates had multiple genotypes and the overall mean multiplicity of infection was 3.2 (95% CI 2.87–3.46). The heterozygosity indices were 0.43 and 0.85 for msp-1 and msp-2, respectively. There was no significant association between multiplicity of infection and age or parasite density.

Conclusions

The study revealed high levels of genetic diversity and mixed-strain infections of P. falciparum populations in Chewaka district, Ethiopia, suggesting that both endemicity level and malaria transmission remain high and that strengthened control efforts are needed in Ethiopia.

Dynamics of Plasmodium vivax populations in border areas of the Greater Mekong sub-region during malaria elimination

Background

Countries within the Greater Mekong Sub-region (GMS) of Southeast Asia have committed to eliminating malaria by 2030. Although the malaria situation has greatly improved, malaria transmission remains at international border regions. In some areas, Plasmodium vivax has become the predominant parasite. To gain a better understanding of transmission dynamics, knowledge on the changes of P. vivax populations after the scale-up of control interventions will guide more effective targeted control efforts.

Methods

This study investigated genetic diversity and population structures in 206 P. vivax clinical samples collected at two time points in two international border areas: the China-Myanmar border (CMB) (n = 50 in 2004 and n = 52 in 2016) and Thailand-Myanmar border (TMB) (n = 50 in 2012 and n = 54 in 2015). Parasites were genotyped using 10 microsatellite markers.

Results

Despite intensified control efforts, genetic diversity remained high (H_E = 0.66–0.86) and was not significantly different among the four populations (P > 0.05). Specifically, H_E slightly decreased from 0.76 in 2004 to 0.66 in 2016 at the CMB and increased from 0.80 in 2012 to 0.86 in 2015 at the TMB. The proportions of polyclonal infections varied significantly among the four populations (P < 0.05), and showed substantial decreases from 48.0% in 2004 to 23.7 at the CMB and from 40.0% in 2012 to 30.7% in 2015 at the TMB, with corresponding decreases in the multiplicity of infection. Consistent with the continuous decline of malaria incidence in the GMS over time, there were also increases in multilocus linkage disequilibrium, suggesting more fragmented and increasingly inbred parasite populations. There were considerable genetic differentiation and sub-division among the four tested populations. Temporal genetic differentiation was observed at each site (F_ST = 0.081 at the CMB and F_ST = 0.133 at the TMB). Various degrees of clustering were evident between the older parasite samples collected in 2004 at the CMB and the 2016 CMB and 2012 TMB populations, suggesting some of these parasites had shared ancestry. In contrast, the 2015 TMB population was genetically distinctive, which may reflect a process of population replacement. Whereas the effective population size (N_e) at the CMB showed a decrease from 4979 in 2004 to 3052 in 2016 with the infinite allele model, the N_e at the TMB experienced an increase from 6289 to 10,259.

Conclusions

With enhanced control efforts on malaria, P. vivax at the TMB and CMB showed considerable spatial and temporal differentiation, but the presence of large P. vivax reservoirs still sustained genetic diversity and transmission. These findings provide new insights into P. vivax transmission dynamics and population structure in these border areas of the GMS. Coordinated and integrated control efforts on both sides of international borders are essential to reach the goal of regional malaria elimination.

Evolution of the Plasmodium vivax multidrug resistance 1 gene in the Greater Mekong Subregion during malaria elimination

Background

The malaria elimination plan of the Greater Mekong Subregion (GMS) is jeopardized by the increasing number of Plasmodium vivax infections and emergence of parasite strains with reduced susceptibility to the frontline drug treatment chloroquine/primaquine. This study aimed to determine the evolution of the P. vivax multidrug resistance 1 (Pvmdr1) gene in P. vivax parasites isolated from the China–Myanmar border area during the major phase of elimination.

Methods

Clinical isolates were collected from 275 P. vivax patients in 2008, 2012–2013 and 2015 in the China–Myanmar border area and from 55 patients in central China. Comparison was made with parasites from three border regions of Thailand.

Results

Overall, genetic diversity of the Pvmdr1 was relatively high in all border regions, and over the seven years in the China–Myanmar border, though slight temporal fluctuation was observed. Single nucleotide polymorphisms previously implicated in reduced chloroquine sensitivity were detected. In particular, M908L approached fixation in the China–Myanmar border area. The Y976F mutation sharply decreased from 18.5% in 2008 to 1.5% in 2012–2013 and disappeared in 2015, whereas F1076L steadily increased from 33.3% in 2008 to 77.8% in 2015. While neutrality tests suggested the action of purifying selection on the pvmdr1 gene, several likelihood-based algorithms detected positive as well as purifying selections operating on specific amino acids including M908L, T958M and F1076L. Fixation and selection of the nonsynonymous mutations are differently distributed across the three border regions and central China. Comparison with the global P. vivax populations clearly indicated clustering of haplotypes according to geographic locations. It is noteworthy that the temperate-zone parasites from central China were completely separated from the parasites from other parts of the GMS.

Conclusions

This study showed that P. vivax populations in the China–Myanmar border has experienced major changes in the Pvmdr1 residues proposed to be associated with chloroquine resistance, suggesting that drug selection may play an important role in the evolution of this gene in the parasite populations.

Microsatellite analysis reveals connectivity among geographically distant transmission zones of Plasmodium vivax in the Peruvian Amazon: A critical barrier to regional malaria elimination

Despite efforts made over decades by the Peruvian government to eliminate malaria, Plasmodium vivax remains a challenge for public health decision-makers in the country. The uneven distribution of its incidence, plus its complex pattern of dispersion, has made ineffective control measures based on global information that lack the necessary detail to understand transmission fully. In this sense, population genetic tools can complement current surveillance. This study describes the genetic diversity and population structure from September 2012 to March 2015 in three geographically distant settlements, Cahuide (CAH), Lupuna (LUP) and Santa Emilia (STE), located in the Peruvian Amazon. A total 777 P. vivax mono-infections, out of 3264, were genotyped. Among study areas, LUP showed 19.7% of polyclonal infections, and its genetic diversity (H_exp) was 0.544. Temporal analysis showed a significant increment of polyclonal infections and H_exp, and the introduction and persistence of a new parasite population since March 2013. In STE, 40.1% of infections were polyclonal, with H_exp = 0.596. The presence of four genetic clusters without signals of clonal expansion and infections with lower parasite densities compared against the other two areas were also found. At least four parasite populations were present in CAH in 2012, where, after June 2014, malaria cases decreased from 213 to 61, concomitant with a decrease in polyclonal infections (from 0.286 to 0.18), and expectedly variable H_exp. Strong signals of gene flow were present in the study areas and wide geographic distribution of highly diverse parasite populations were found. This study suggests that movement of malaria parasites by human reservoirs connects geographically distant malaria transmission areas in the Peruvian Amazon. The maintenance of high levels of parasite genetic diversity through human mobility is a critical barrier to malaria elimination in this region.

Plasmodium vivax transmission is heterogeneous and discontinuous in the Peruvian Amazon. Such heterogeneity is the result of factors that include, but are not restricted to, the environment, public policies, and characteristics of the parasite, the vector, and human activities. All these factors make P. vivax transmission resilient to interventions. In order to achieve the goals of control and local elimination, P. vivax surveillance must inform how those factors sustain disease transmission in order to focalize and synchronize control strategies. In this study, we implemented molecular surveillance complemented with population genetic tools in the areas of Cahuide, Lupuna, and Santa Emilia located in the Peruvian Amazon. In particular, we characterize the transmission and the parasite genetic variation in these sites from September 2012 to March 2015. The changes in parasite diversity, the wide geographic dispersion of parasite subpopulation and the introduction of a new parasite clone or subpopulation in Lupuna documented in this study suggest that connectivity among the different endemic areas, likely due to human mobility, sustains disease transmission in the region hindering the success of control measures. This information must be considered in the design of current control strategies.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

Genetic polymorphism of Merozoite Surface Protein 1 (msp1) and 2 (msp2) genes and multiplicity of Plasmodium falciparum infection across various endemic areas in Senegal

INTRODUCTION

Despite a significant decline in Senegal, malaria remains a burden in various parts of the country. Assessment of multiplicity of Plasmodium falciparum infection and genetic diversity of parasites population could help in monitoring of malaria control.

OBJECTIVE

To assess genetic diversity and multiplicity of infection in P. falciparum isolates from three areas in Senegal with different malaria transmissions.

METHODS

136 blood samples were collected from patients with uncomplicated P. falciparum malaria in Pikine, Kedougou and Thies. Polymorphic loci of msp1 and 2 (Merozoite surface protein-1 and 2) genes were amplified by nested PCR.

RESULTS

For msp1gene, K1 allelic family was predominant with frequency of 71%. Concerning msp2 gene, IC3D7 allelic family was the most represented with frequency of 83%. Multiclonal isolates found were 36% and 31% for msp1et msp2 genes respectively. The MOI found in all areas was 2.56 and was statistically different between areas (P=0.024). Low to intermediate genetic diversity were found with heterozygosity range (He=0,394-0,637) and low genetic differentiation (Fst msp1= 0.011; Fst msp2=0.017) were observed between P. falciparum population within the country.

CONCLUSION

Low to moderate genetic diversity of P.falciparum strains and MOI disparities were found in Senegal.

Malaria Situation in Latin America and the Caribbean: Residual and Resurgent Transmission and Challenges for Control and Elimination

Despite recent progress toward malaria elimination in Latin America and the Caribbean, with an overall 62% decrease in incidence between 2000 and 2015, malaria remains endemic to 21 countries and territories in the region, where 120 million people are exposed to some risk of infection. Here we review recent epidemiologic trends, highlight current challenges, and briefly discuss the relative role of traditional and novel strategies for better malaria control and elimination across the continent.

Genetic Diversity of Plasmodium vivax Causing Epidemic Malaria in the Republic of Korea

Plasmodium vivax is more challenging to control and eliminate than P. falciparum due to its more asymptomatic infections with low parasite densities making diagnosis more difficult, in addition to its unique biological characteristics. The potential re-introduction of incidence cases, either through borders or via human migrations, is another major hurdle to sustained control and elimination. The Republic of Korea has experienced re-emergence of vivax malaria in 1993 but is one of the 32 malaria-eliminating countries to-date. Despite achieving successful nationwide control and elimination of vivax malaria, the evolutionary characteristics of vivax malaria isolates in the Republic of Korea have not been fully understood. In this review, we present an overview of the genetic variability of such isolates to increase understanding of the epidemiology, diversity, and dynamics of vivax populations in the Republic of Korea.

Plasmodium genomics: an approach for learning about and ending human malaria

Malaria causes high levels of morbidity and mortality in human beings worldwide. According to the World Health Organization (WHO), about half a million people die of this disease each year. Malaria is caused by six species of parasites belonging to the Plasmodium genus: P. falciparum, P. knowlesi, P. vivax, P. malariae, P. ovale curtisi, and P. ovale wallikeri. Currently, malaria is being kept under control with varying levels of elimination success in different countries. The development of new molecular tools as well as the use of next-generation sequencing (NGS) technologies and novel bioinformatic approaches has improved our knowledge of malarial epidemiology, diagnosis, treatment, vaccine development, and surveillance strategies. In this work, the genetics and genomics of human malarias have been analyzed. Since the first P. falciparum genome was sequenced in 2002, various population-level genetic and genomic surveys, together with transcriptomic and proteomic studies, have shown the importance of molecular approaches in supporting malaria elimination.

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.

Seasonal variations in Plasmodium falciparum genetic diversity and multiplicity of infection in asymptomatic children living in southern Ghana

Background

Genetic diversity in Plasmodium falciparum (P. falciparum) parasites is a major hurdle to the control of malaria. This study monitored changes in the genetic diversity and the multiplicity of P. falciparum parasite infection in asymptomatic children living in southern Ghana at 3 month intervals between April 2015 and January 2016.

Methods

Filter paper blood spots (DBS) were collected quarterly from children living in Obom, a community with perennial malaria transmission and Abura, a community with seasonal malaria transmission. Genomic DNA was extracted from the DBS and used in polymerase chain reaction (PCR)-based genotyping of the merozoite surface protein 1 (msp 1) and merozoite surface protein 2 (msp 2) genes.

Results

Out of a total of 787 samples that were collected from the two study sites, 59.2% (466/787) tested positive for P. falciparum. The msp 1 and msp 2 genes were successfully amplified from 73.8% (344/466) and 82.5% (385/466) of the P. falciparum positive samples respectively. The geometric mean MOI in Abura ranged between 1.17 (95% CI: 1.08–1.28) and 1.48 (95% CI: 1.36–1.60) and was significantly lower (p < 0.01, Dunn’s multiple comparison test) than that determined in Obom, where the geometric mean MOI ranged between 1.82 (95% CI: 1.58–2.08) and 2.50 (95% CI: 2.33–2.678) over the study period. Whilst the msp 1 R033:MAD20:KI allelic family ratio was dynamic, the msp 2 3D7:FC27 allelic family ratio remained relatively stable across the changing seasons in both sites.

Conclusions

This study shows that seasonal variations in parasite diversity in these communities can be better estimated by msp 1 rather than msp 2 due to the constantly changing relative intra allelic frequencies observed in msp 1 and the fact that the dominance of any msp 2 allele was dependent on the transmission setting but not on the season as opposed to the dominance of any msp 1 allele, which was dependent on both the season and the transmission setting.

Electronic supplementary material

The online version of this article (10.1186/s12879-018-3350-z) contains supplementary material, which is available to authorized users.

Seasonal variations in Plasmodium falciparum parasite prevalence assessed by varying diagnostic tests in asymptomatic children in southern Ghana

Plasmodium falciparum infections presenting either as symptomatic or asymptomatic may contain sexual stage parasites (gametocytes) that are crucial to malaria transmission. In this study, the prevalence of microscopic and submicroscopic asexual and gametocyte parasite stages were assessed in asymptomatic children from two communities in southern Ghana. Eighty children aged twelve years and below, none of whom exhibited signs of clinical malaria living in Obom and Cape Coast were sampled twice, one during the rainy (July 2015) and subsequently during the dry (January 2016) season. Venous blood was used to prepare thick and thin blood smears, spot a rapid malaria diagnostic test (PfHRP2 RDT) as well as prepare filter paper blood spots. Blood cell pellets were preserved in Trizol for RNA extraction. Polymerase chain reaction (PCR) and semi-quantitative real time reverse transcriptase PCR (qRT-PCR) were used to determine submicroscopic parasite prevalence. In both sites 87% (95% CI: 78-96) of the asymptomatic individuals surveyed were parasites positive during the 6 month study period. The prevalence of asexual and gametocyte stage parasites in the rainy season were both significantly higher in Obom than in Cape Coast (P < 0.001). Submicroscopic gametocyte prevalence was highest in the rainy season in Obom but in the dry season in Cape Coast. Parasite prevalence determined by PCR was similar to that determined by qRT-PCR in Obom but significantly lower than that determined by qRT-PCR in Cape Coast. Communities with varying parasite prevalence exhibit seasonal variations in the prevalence of gametocyte carriers. Submicroscopic asymptomatic parasite and gametocyte carriage is very high in southern Ghana, even during the dry season in communities with low microscopic parasite prevalence and likely to be missed during national surveillance exercises.

Genetic diversity of the Plasmodium vivax phosphatidylinositol 3-kinase gene in two regions of the China-Myanmar border

Artemisinin resistance in Plasmodium falciparum was associated with mutations in the propeller domain of the PfK13 gene and increased phosphatidylinositol-3'-kinase (PfPI3K) activity. Assessment of the genetic diversity of the PfK13 ortholog PvK12 in Plasmodium vivax field samples from the same hotspots of P. falciparum artemisinin resistance revealed a limited genetic diversity of PvK12. Following the same logic, we analyzed genetic variations of the PvPI3K gene in 188 P. vivax field isolates from two geographic locations along the China-Myanmar border. Overall, high genetic diversity of PvPI3K was observed; parasites from Yunnan's Tengchong County had higher genetic diversity than those from Laiza Township, Kachin State, Myanmar. Almost all the neutrality tests applied detected statistically significant deviation from zero. The negative Tajima's D values in both populations implicated that PvPI3K gene might have experienced either a directional selection or an expansion in population size. There was low linkage disequilibrium between the PvPI3K mutations in both populations, suggesting the existence of large, almost panmictic, parasite populations that enabled effective recombination. This later result was confirmed by the detection of a minimum of five recombination events in each population with two major breakpoints. Multiple tests for selection confirmed a signature of purifying selection on PvPI3K. All the amino acid mutations were predicted to be neutral for the PI3K protein's function. These findings provide insights on the genetic diversity of P. vivax populations along the China-Myanmar border.

Microbial Pre-exposure and Vectorial Competence of Anopheles Mosquitoes

Anopheles female mosquitoes can transmit Plasmodium, the malaria parasite. During their aquatic life, wild Anopheles mosquito larvae are exposed to a huge diversity of microbes present in their breeding sites. Later, adult females often take successive blood meals that might also carry different micro-organisms, including parasites, bacteria, and viruses. Therefore, prior to Plasmodium ingestion, the mosquito biology could be modulated at different life stages by a suite of microbes present in larval breeding sites, as well as in the adult environment. In this article, we highlight several naturally relevant scenarios of Anopheles microbial pre-exposure that we assume might impact mosquito vectorial competence for the malaria parasite: (i) larval microbial exposures; (ii) protist co-infections; (iii) virus co-infections; and (iv) pathogenic bacteria co-infections. In addition, significant behavioral changes in African Anopheles vectors have been associated with increasing insecticide resistance. We discuss how these ethological modifications may also increase the repertoire of microbes to which mosquitoes could be exposed, and that might also influence their vectorial competence. Studying Plasmodium–Anopheles interactions in natural microbial environments would efficiently contribute to refining the transmission risks.

Nationwide genetic surveillance of Plasmodium vivax in Papua New Guinea reveals heterogeneous transmission dynamics and routes of migration amongst subdivided populations

The Asia Pacific Leaders in Malaria Alliance (APLMA) have committed to eliminate malaria from the region by 2030. Papua New Guinea (PNG) has the highest malaria burden in the Asia-Pacific region but with the intensification of control efforts since 2005, transmission has been dramatically reduced and Plasmodium vivax is now the dominant malaria infection in some parts of the country. To gain a better understanding of the transmission dynamics and migration patterns of P. vivax in PNG, here we investigate population structure in eight geographically and ecologically distinct regions of the country. A total of 219 P. vivax isolates (16-30 per population) were successfully haplotyped using 10 microsatellite markers. A wide range of genetic diversity (H_e=0.37-0.87, R_s=3.60-7.58) and significant multilocus linkage disequilibrium (LD) was observed in six of the eight populations (I_A^S=0.08-0.15 p-value<0.05) reflecting a spectrum of transmission intensities across the country. Genetic differentiation between regions was evident (Jost's D=0.07-0.72), with increasing divergence of populations with geographic distance. Overall, P. vivax isolates clustered into three major genetic populations subdividing the Mainland lowland and coastal regions, the Islands and the Highlands. P. vivax gene flow follows major human migration routes, and there was higher gene flow amongst Mainland parasite populations than among Island populations. The Central Province (samples collected in villages close to the capital city, Port Moresby), acts as a sink for imported infections from the three major endemic areas. These insights into P. vivax transmission dynamics and population networks will inform targeted strategies to contain malaria infections and to prevent the spread of drug resistance in PNG.

Population-level estimates of the proportion of Plasmodium vivax blood-stage infections attributable to relapses among febrile patients attending Adama Malaria Diagnostic Centre, East Shoa Zone, Oromia, Ethiopia

Background

Malaria is ranked as the leading communicable disease in Ethiopia, where Plasmodium falciparum and Plasmodium vivax are co-endemic. The incidence of P. vivax is usually considered to be less seasonal than P. falciparum. Clinical cases of symptomatic P. falciparum exhibit notable seasonal variation, driven by rainfall-dependent variation in the abundance of Anopheles mosquitoes. A similar peak of clinical cases of P. vivax is usually observed during the rainy season. However, the ability of P. vivax to relapse causing new blood-stage infections weeks to months after an infectious mosquito bite can lead to substantial differences in seasonal patterns of clinical cases. These cannot be detected with currently available diagnostic tools and are not cleared upon treatment with routinely administered anti-malarial drugs.

Methods

A health- facility based cross-sectional study was conducted in Adama malaria diagnostic centre from May 2015 to April 2016. Finger-prick blood samples were collected for thin and thick blood film preparation from participants seeking treatment for suspected cases of febrile malaria. Informed consent was obtained from each study participant or their guardians. Seasonal patterns in malaria cases were analysed using statistical models, identifying the peaks in cases, and the seasonally varying proportion of P. vivax cases attributable to relapses.

Results

The proportion of patients with malaria detectable by light microscopy was 36.1% (1141/3161) of which P. vivax, P. falciparum, and mixed infections accounted for 71.4, 25.8 and 2.8%, respectively. Of the febrile patients diagnosed, 2134 (67.5%) were males and 1919 (60.7%) were urban residents. The model identified a primary peak in P. falciparum and P. vivax cases from August to October, as well as a secondary peak of P. vivax cases from February to April attributable to cases arising from relapses. During the secondary peak of P. vivax cases approximately 77% (95% CrI 68, 84%) of cases are estimated to be attributable to relapses. During the primary peak from August to October, approximately 40% (95% CrI 29, 57%) of cases are estimated to be attributable to relapses.

Discussion

It is not possible to diagnose whether a P. vivax case has been caused by blood-stage infection from a mosquito bite or a relapse. However, differences in seasonal patterns of P. falciparum and P. vivax cases can be used to estimate the population-level proportion of P. vivax cases attributable to relapses. These observations have important implications for the epidemiological assessment of vivax malaria, and initiating therapy that is effective against both blood stages and relapses.

Electronic supplementary material

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

Identification of the Plasmodium species in clinical samples from children residing in five epidemiological strata of malaria in Cameroon

Background

Malaria in Cameroon was previously known to be caused solely by Plasmodium falciparum but today, evidence points to other Plasmodium species including P. vivax, P. ovale and P. malariae. The purpose of this study was to identify the Plasmodium species in clinical samples from children residing in five epidemiological strata of malaria in Cameroon, so as to advise control policies.

Methods

One thousand six hundred nine febrile children (≤15 years) were recruited from five epidemiological strata of malaria including the Sudano-sahelian (SS) strata, the High inland plateau (HIP) strata, the South Cameroonian Equatorial forest (SCEF) strata, the High western plateau (HWP) strata and the Coastal (C) strata. Malaria parasites were detected by Giemsa microscopy (GM) while a multiplex polymerase chain reaction (PCR) was used to identify the Plasmodium species. Statistical analysis performed included the Pearson chi-square test, and statistical significance was set at p < 0.05.

Results

The PCR-adjusted prevalence of malaria was 17.6%. The detection rate of PCR was higher than GM (p = 0.05). However, GM demonstrated a high sensitivity (85.5%) and specificity (100%) and, overall, a perfectly correlated agreement with PCR (97.5%). The prevalence of malaria was significantly higher in children between 60 and 119 months (p < 0.001) and in Limbe (in the Coastal strata) (p < 0.001). Contrariwise, the prevalence of malaria was not associated with gender (p = 0.239). P. falciparum was identified in all (100%) the cases of malaria; P. ovale, P. vivax, P. malariae and P. knowlesi were all absent. No case of mixed infection was identified.

Conclusions

P. falciparum was the only species causing clinical malaria in the target population, which is contrary to studies that have reported P. vivax, P. malariae and P. ovale as causing clinical malaria in Cameroon.

Genomic variation in Plasmodium vivax malaria reveals regions under selective pressure

BACKGROUND

Although Plasmodium vivax contributes to almost half of all malaria cases outside Africa, it has been relatively neglected compared to the more deadly P. falciparum. It is known that P. vivax populations possess high genetic diversity, differing geographically potentially due to different vector species, host genetics and environmental factors.

RESULTS

We analysed the high-quality genomic data for 46 P. vivax isolates spanning 10 countries across 4 continents. Using population genetic methods we identified hotspots of selection pressure, including the previously reported MRP1 and DHPS genes, both putative drug resistance loci. Extra copies and deletions in the promoter region of another drug resistance candidate, MDR1 gene, and duplications in the Duffy binding protein gene (PvDBP) potentially involved in erythrocyte invasion, were also identified. For surveillance applications, continental-informative markers were found in putative drug resistance loci, and we show that organellar polymorphisms could classify P. vivax populations across continents and differentiate between Plasmodia spp.

CONCLUSIONS

This study has shown that genomic diversity that lies within and between P. vivax populations can be used to elucidate potential drug resistance and invasion mechanisms, as well as facilitate the molecular barcoding of the parasite for surveillance applications.

Higher Complexity of Infection and Genetic Diversity of Plasmodium vivax Than Plasmodium falciparum Across All Malaria Transmission Zones of Papua New Guinea

Plasmodium falciparum and Plasmodium vivax have varying transmission dynamics that are informed by molecular epidemiology. This study aimed to determine the complexity of infection and genetic diversity of P. vivax and P. falciparum throughout Papua New Guinea (PNG) to evaluate transmission dynamics across the country. In 2008-2009, a nationwide malaria indicator survey collected 8,936 samples from all 16 endemic provinces of PNG. Of these, 892 positive P. vivax samples were genotyped at PvMS16 and PvmspF3, and 758 positive P. falciparum samples were genotyped at Pfmsp2. The data were analyzed for multiplicity of infection (MOI) and genetic diversity. Overall, P. vivax had higher polyclonality (71%) and mean MOI (2.32) than P. falciparum (20%, 1.39). These measures were significantly associated with prevalence for P. falciparum but not for P. vivax. The genetic diversity of P. vivax (PvMS16: expected heterozygosity = 0.95, 0.85-0.98; PvMsp1F3: 0.78, 0.66-0.89) was higher and less variable than that of P. falciparum (Pfmsp2: 0.89, 0.65-0.97). Significant associations of MOI with allelic richness (rho = 0.69, P = 0.009) and expected heterozygosity (rho = 0.87, P < 0.001) were observed for P. falciparum. Conversely, genetic diversity was not correlated with polyclonality nor mean MOI for P. vivax. The results demonstrate higher complexity of infection and genetic diversity of P. vivax across the country. Although P. falciparum shows a strong association of these parameters with prevalence, a lack of association was observed for P. vivax and is consistent with higher potential for outcrossing of this species.

VivaxGEN: An open access platform for comparative analysis of short tandem repeat genotyping data in Plasmodium vivax populations

BACKGROUND

The control and elimination of Plasmodium vivax will require a better understanding of its transmission dynamics, through the application of genotyping and population genetics analyses. This paper describes VivaxGEN (http://vivaxgen.menzies.edu.au), a web-based platform that has been developed to support P. vivax short tandem repeat data sharing and comparative analyses.

RESULTS

The VivaxGEN platform provides a repository for raw data generated by capillary electrophoresis (FSA files), with fragment analysis and standardized allele calling tools. The query system of the platform enables users to filter, select and differentiate samples and alleles based on their specified criteria. Key population genetic analyses are supported including measures of population differentiation (FST), expected heterozygosity (HE), linkage disequilibrium (IAS), neighbor-joining analysis and Principal Coordinate Analysis. Datasets can also be formatted and exported for application in commonly used population genetic software including GENEPOP, Arlequin and STRUCTURE. To date, data from 10 countries, including 5 publicly available data sets have been shared with VivaxGEN.

CONCLUSIONS

VivaxGEN is well placed to facilitate regional overviews of P. vivax transmission dynamics in different endemic settings and capable to be adapted for similar genetic studies of P. falciparum and other organisms.

Size polymorphism and low sequence diversity in the locus encoding the Plasmodium vivax rhoptry neck protein 4 (PvRON4) in Colombian isolates

Background

Designing a vaccine against Plasmodium vivax has focused on selecting antigens involved in invasion mechanisms that must have domains with low polymorphism for avoiding allele-specific immune responses. The rhoptry neck protein 4 (RON4) forms part of the tight junction, which is essential in the invasion of hepatocytes and/or erythrocytes; however, little is known about this locus’ genetic diversity.

Methods

DNA sequences from 73 Colombian clinical isolates from pvron4 gene were analysed for characterizing their genetic diversity; pvron4 haplotype number and distribution, as well as the evolutionary forces determining diversity pattern, were assessed by population genetics and molecular evolutionary approaches.

Results

ron4 has low genetic diversity in P. vivax at sequence level; however, a variable amount of tandem repeats at the N-terminal region leads to extensive size polymorphism. This region seems to be exposed to the immune system. The central region has a putative esterase/lipase domain which, like the protein’s C-terminal fragment, is highly conserved at intra- and inter-species level. Both regions are under purifying selection.

Conclusions

pvron4 is the locus having the lowest genetic diversity described to date for P. vivax. The repeat regions in the N-terminal region could be associated with immune evasion mechanisms while the central region and the C-terminal region seem to be under functional or structural constraint. Bearing such results in mind, the PvRON4 central and/or C-terminal portions represent promising candidates when designing a subunit-based vaccine as they are aimed at avoiding an allele-specific immune response, which might limit vaccine efficacy.

Electronic supplementary material

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

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 Structure and Inhibition of a Malaria Drug Target: Geranylgeranyl Diphosphate Synthase

We report a molecular dynamics investigation of the structure, function, and inhibition of geranylgeranyl diphosphate synthase (GGPPS), a potential drug target, from the malaria parasite Plasmodium vivax. We discovered several GGPPS inhibitors, benzoic acids, and determined their structures crystallographically. We then used molecular dynamics simulations to investigate the dynamics of three such inhibitors and two bisphosphonate inhibitors, zoledronate and a lipophilic analogue of zoledronate, as well as the enzyme's product, GGPP. We were able to identify the main motions that govern substrate binding and product release as well as the molecular features required for GGPPS inhibition by both classes of inhibitor. The results are of broad general interest because they represent the first detailed investigation of the mechanism of action, and inhibition, of an important antimalarial drug target, geranylgeranyl diphosphate synthase, and may help guide the development of other, novel inhibitors as new drug leads.

Population genetics structure of Plasmodium vivax circumsporozoite protein during the elimination process in low and unstable malaria transmission areas, southeast of Iran

In Iran, the prevalence of Plasmodium falciparum and Plasmodium vivax has dropped after a national malaria elimination program was launched. To estimate the likelihood of success and to measure the outcome of malaria intervention tools during elimination programs (2008-2012), the population genetic surveys of Iranian P. vivax isolates (n=60) were carried out using the CSP genetic marker. The results were compared with a similar work that was carried out during a control phase (2000-2003) in the same study areas. Based on PCR-RFLP analysis, 49 (81.67%) of 60 studied samples were VK210 and 11 (18.33%) were VK247 with no mixed genotypes. However, 10.97% of P. vivax isolates of control phase harbored the mixed genotypes. Sequencing analysis of 50 pvcsp gene showed 14 distinct haplotypes, of which 11 and 3 were VK210 and VK247 types, respectively. However, during the control phase, 19 distinct subtypes (11 VK210 and 8 VK247) were reported. Also, 7 of 11 VK210 and the VK247F subtypes were new, and 3 out of 7 new VK210 and VK247F were isolated from the patients with Pakistani nationality. The lower nucleotide diversity per site (π=0.02017±0.00436 and π=0.04525±0.00255) and haplotype diversity (Hd=0.513±0.093 and Hd=0.691±0.128) as well as lower In/Del haplotype [Hd(i)=0.243 and 0] and nucleotide diversity [π(i)=0.00078 and 0] were recorded for VK210 and VK247of the elimination samples, respectively. In conclusion, the comparison of PRMs and RATs in CRR along with the polymorphism analysis of the sequence lengths, SNPs, and In/Del polymorphisms in all analyzed samples showed lower genetic diversity for PvCSP in the elimination samples. Also, although there is a turnover of P. vivax parasite genotypes in the study areas, reduction in genetic diversity and transmission was detected due to scaling-up of the intervention tools during an elimination program in Iran. This notable challenge of the elimination program must be taken into account and controlled by active surveillance for limiting both reintroductions of new allelic forms as well as the spread of drug-resistant parasite to prevent any disease outbreaks.

Worldwide population genetic analysis and natural selection in the Plasmodium vivax Generative Cell Specific 1 (PvGCS1) as a transmission-blocking vaccine candidate

GENERATIVE CELL SPECIFIC 1 (GCS1) is one of the Transmission Blocking Vaccine (TBV) candidate antigens, which is expressed on the surface of male gametocytes and gametes of Plasmodium species. Since antigenic diversity could inhibit the successful development of a malaria vaccine, it is crucial to determine the diversity of gcs1 gene in global malaria-endemic areas. Therefore, gene diversity and selection of gcs1 gene were analyzed in Iranian Plasmodium vivax isolates (n=52) and compared with the corresponding sequences from worldwide clinical P. vivax isolates available in PlasmoDB database. Totally 12 SNPs were detected in the pvgcs1 sequences as compared to Sal-1 sequence. Five out of 12 SNPs including three synonymous (T797C, G1559A, and G1667T) and two amino acid replacements (Y133S and Q634P) were detected in Iranian pvgcs1 sequences. According to four amino acid replacements (Y133S, N575S, Q634P and D637N) observed in all world PvGCS1 sequences, totally 5 PvGCS1 haplotypes were detected in the world, that three of them observed in Iranian isolates including the PvGCS-A (133S/634Q, 92.3%), PvGCS-B (133Y/634Q, 5.8%), and PvGCS-C (133S/634P, 1.9%). The overall nucleotide diversity (π) for all 52 sequences of Iranian pvgcs1 gene was 0.00018±0.00006, and the value of dN-dS (-0.00031) were negative, however, it was not statistically significant. In comparison with global isolates, Iranian and PNG pvgcs1 sequences had the lowest nucleotide and haplotype diversity, while the highest nucleotide and haplotype diversity was observed in China population. Moreover, epitope prediction in this antigen showed that all B-cell epitopes were located in conserved regions. However, Q634P (in one Iranian isolate) and D637N (observed in Thailand, China, Vietnam and North Korea) mutations are involved in predicted IURs. The obtained results in this study could be used in development of PvGCS1 based malaria vaccine.