Population Genetics of Plasmodium vivax in Four High Malaria Endemic Areas in Thailand

Extensive genetic diversity in Plasmodium vivax from Sudan and its genetic relationships with other geographical isolates

Plasmodium vivax, traditionally overlooked has experienced a notable increase in cases in East Africa. This study investigated the geographical origin and genetic diversity of P. vivax in Sudan using 14 microsatellite markers. A total of 113 clinical P. vivax samples were collected from two different ecogeographical zones, New Halfa and Khartoum, in Sudan. Additionally, 841 geographical samples from the database were incorporated for a global genetic analysis to discern genetic relationships among P. vivax isolates on regional and worldwide scales. On the regional scale, our findings revealed 91 unique and 8 shared haplotypes among the Sudan samples, showcasing a remarkable genetic diversity compared to other geographical isolates and supporting the hypothesis that P. vivax originated from Africa. On a global scale, distinct genetic clustering of P. vivax isolates from Africa, South America, and Asia (including Papua New Guinea and Solomon Island) was observed, with limited admixture among the three clusters. Principal component analysis emphasized the substantial contribution of African isolates to the observed global genetic variation. The Sudanese populations displayed extensive genetic diversity, marked by significant multi-locus linkage disequilibrium, suggesting an ancestral source of P. vivax variation globally and frequent recombination among the isolates. Notably, the East African P. vivax exhibited similarity with some Asian isolates, indicating potential recent introductions. Overall, our results underscore the effectiveness of utilizing microsatellite markers for implementing robust control measures, given their ability to capture extensive genetic diversity and linkage disequilibrium patterns.

Plasmodium vivax populations in the western Greater Mekong Subregion evaluated using a genetic barcode

An improved understanding of the Plasmodium vivax populations in the Great Mekong Subregion (GMS) is needed to monitor the progress of malaria elimination. This study aimed to use a P. vivax single nucleotide polymorphism (SNP) barcode to evaluate the population dynamics and explore the gene flow among P. vivax parasite populations in the western GMS (China, Myanmar and Thailand). A total of 315 P. vivax patient samples collected in 2011 and 2018 from four regions of the western GMS were genotyped for 42 SNPs using the high-throughput MassARRAY SNP genotyping technology. Population genetic analysis was conducted to estimate the genetic diversity, effective population size, and population structure among the P. vivax populations. Overall, 291 samples were successfully genotyped at 39 SNPs. A significant difference was observed in the proportion of polyclonal infections among the five P. vivax populations (P = 0.0012, Pearson Chi-square test, χ2 = 18.1), with western Myanmar having the highest proportion (96.2%, 50/52) in 2018. Likewise, the average complexity of infection was also highest in western Myanmar (1.31) and lowest in northeast Myanmar (1.01) in 2018. The older samples from western China in 2011 had the highest pairwise nucleotide diversity (π, 0.388 ± 0.046), expected heterozygosity (He, 0.363 ± 0.02), and the largest effective population size. In comparison, in the neighboring northeast Myanmar, the more recent samples in 2018 showed the lowest values (π, 0.224 ± 0.036; He, 0.220 ± 0.026). Furthermore, the 2018 northeast Myanmar parasites showed high and moderate genetic differentiation from other populations with FST values of 0.162-0.252, whereas genetic differentiation among other populations was relatively low (FST ≤ 0.059). Principal component analysis, phylogeny, and STRUCTURE analysis showed that the P. vivax population in northeast Myanmar in 2018 substantially diverged from other populations. Although the 42 SNP barcode is a valuable tool for tracking parasite origins of worldwide parasite populations, a more extended barcode with additional SNPs is needed to distinguish the more related parasite populations in the western GMS.

Low Genetic Diversity of Plasmodium vivax Circumsporozoite Surface Protein in Clinical Isolates from Southern Thailand

The genetic diversity within the circumsporozoite surface protein (PvCSP) of Plasmodium vivax, the predominant malaria species in Thailand, is primarily observed in the northwestern region along the Thailand-Myanmar border. However, as P. vivax cases shift to southern provinces, particularly Yala Province near the Thailand-Malaysia border, PvCSP diversity remains understudied. Between 2018 and 2020, 89 P. vivax isolates were collected in Yala Province, a significant malaria hotspot. Employing polymerase chain reaction amplification, restriction fragment length polymorphism (PCR-RFLP), and DNA sequencing, the gene encoding PvCSP (Pvcsp) was analyzed. All Yala P. vivax isolates belonged to the VK210 type, distinct from strains in the western region near the Myanmar border. The central repeat region of Pvcsp revealed two common peptide repeat motifs-GDRADGQPA and GDRAAGQPA-across all southern isolates. Sequence analysis identified two subtypes, with S1 more prevalent (92%) than S2 (8%). This study underscores the limited diversity of VK210 variants of P. vivax populations in southern Thailand. These baseline findings facilitate monitoring for potential new parasite variants, aiding in the future control and management of P. vivax in the region.

Comparison of Vector Trapping Methods for Outdoor Biting Malaria Vector Surveillance in Thailand and Vietnam

The performances of the human-baited double net trap (HDNT) and the human-baited host decoy trap (HDT) methods were compared against the outdoor human landing catch (OHLC) method in Thailand and Vietnam. Two study sites were selected in each country: a rural village and a nearby forest setting. The three outdoor trap methods were rotated nightly between three set trapping positions, in a pre-assigned Latin square design. Volunteers were rotated following the trap rotation to avoid bias. The greatest number of adult mosquitoes was collected from the forest sites in both countries, showing Anopheles minimus (s.s.) Theobald (96.54%) and Anopheles dirus (s.s.) Peyton & Harrison (25.71%) as the primary malaria vectors in Thailand and Vietnam, respectively. At the Thai forest site, OHLC collected significantly more anopheline mosquitoes per trap night than HDNT and HDT, with mean ± standard error values of 14.17 ± 4.42, 4.83 ± 1.56, and 4.44 ± 1.45, respectively, whilst HDNT and HDT were significantly less productive at 0.34 times and 0.31 times, respectively, than OHLC in capturing anopheline mosquitoes. However, there were no significant differences among the three methods of trapping malaria vectors for the village site. At the Vietnamese forest site, HDNT achieved the highest performance in collecting Anopheline mosquitoes at 1.54 times compared to OHLC, but there was no significant difference between the two traps. The results suggested HDNT could be a possible alternative trap to OHLC in this area. Although HDT was less efficient at attracting Anopheline mosquitoes, it was highly efficient at trapping culicine mosquitoes.

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.