Frequent Spread of Plasmodium vivax Malaria Maintains High Genetic Diversity at the Myanmar-China Border, Without Distance and Landscape Barriers

1-3-7 surveillance and response approach in malaria elimination: China's practice and global adaptions

There has been a significant reduction in malaria morbidity and mortality worldwide from 2000 to 2019. However, the incidence and mortality increased again in 2020 due to the disruption to services during the COVID-19 pandemic. Surveillance to reduce the burden of malaria, eliminate the disease and prevent its retransmission is, therefore, crucial. The 1-3-7 approach proposed by China has played an important role in eliminating malaria, which has been internationally popularized and adopted in some countries to help eliminate malaria. This review summarizes the experience and lessons of 1-3-7 approach in China and its application in other malaria-endemic countries, so as to provide references for its role in eliminating malaria and preventing retransmission. This approach needs to be tailored and adapted according to the region condition, considering the completion, timeliness and limitation of case-based reactive surveillance and response. It is very important to popularize malaria knowledge, train staff, improve the capacity of health centres and monitor high-risk groups to improve the performance in eliminating settings. After all, remaining vigilance in detecting malaria cases and optimizing surveillance and response systems are critical to achieving and sustaining malaria elimination.

Implementing landscape genetics in molecular epidemiology to determine drivers of vector-borne disease: A malaria case study

This study employs landscape genetics to investigate the environmental drivers of a deadly vector-borne disease, malaria caused by Plasmodium falciparum, in a more spatially comprehensive manner than any previous work. With 1804 samples from 44 sites collected in western Kenya in 2012 and 2013, we performed resistance surface analysis to show that Lake Victoria acts as a barrier to transmission between areas north and south of the Winam Gulf. In addition, Mantel correlograms clearly showed significant correlations between genetic and geographic distance over short distances (less than 70 km). In both cases, we used an identity-by-state measure of relatedness tailored to find highly related individual parasites in order to focus on recent gene flow that is more relevant to disease transmission. To supplement these results, we performed conventional population genetics analyses, including Bayesian clustering methods and spatial ordination techniques. These analyses revealed some differentiation on the basis of geography and elevation and a cluster of genetic similarity in the lowlands north of the Winam Gulf of Lake Victoria. Taken as a whole, these results indicate low overall genetic differentiation in the Lake Victoria region, but with some separation of parasite populations north and south of the Winam Gulf that is explained by the presence of the lake as a geographic barrier to gene flow. We recommend similar landscape genetics analyses in future molecular epidemiology studies of vector-borne diseases to extend and contextualize the results of traditional population genetics.

Different In Vitro Drug Susceptibility Profile of Plasmodium falciparum Isolates from Two Adjacent Areas of Northeast Myanmar and Molecular Markers for Drug Resistance

The Greater Mekong Subregion (GMS) is the epicenter of antimalarial drug resistance. We determined in vitro susceptibilities to 11 drugs of culture-adapted Plasmodium falciparum isolates from adjacent areas (Laiza and Muse) along the China−Myanmar border. Parasites from this region were highly resistant to chloroquine and pyrimethamine but relatively sensitive to other antimalarial drugs. Consistently, the Dd2-like pfcrt mutations were fixed or almost fixed in both parasite populations, and new mutations mediating piperaquine resistance were not identified. Similarly, several mutations related to pfdhfr and pfdhps were also highly prevalent. Despite their geographical proximity, malaria parasites from Laiza showed significantly higher in vitro resistance to artemisinin derivatives, naphthoquine, pyronaridine, lumefantrine, and pyrimethamine than parasites from Muse. Likewise, the pfdhfr N51I, pfdhps A581G, pfmrp1 H785N, and pfk13 F446I mutations were significantly more frequent in Laiza than in Muse (p < 0.05). For the pfmdr1 mutations, Y184F was found only in Laiza (70%), whereas F1226Y was identified only in Muse (31.8%). Parasite isolates from Laiza showed a median RSA value of 5.0%, significantly higher than the 2.4% in Muse. Altogether, P. falciparum parasite populations from neighboring regions in the GMS may diverge substantially in their resistance to several antimalarial drugs. This information about different parasite populations will guide antimalarial treatment policies to effectively manage drug resistance during malaria elimination.

PvMSP-3α and PvMSP-3β genotyping reveals higher genetic diversity in Plasmodium vivax parasites from migrant workers than residents at the China-Myanmar border

BACKGROUND

The genetic diversity of malaria parasites traces the origin and spread of new variants and can be used to evaluate the effectiveness of malaria control measures. Therefore, this study aims to improve the understanding of the molecular epidemiology of Plasmodium vivax malaria at the China-Myanmar border by genotyping the PvMSP-3α and PvMSP-3β genes.

METHODS

Blood samples were collected from P. vivax malaria patients along the China-Myanmar border. The PvMSP-3α and PvMSP-3β genes were amplified by polymerase chain reaction (PCR) and the genetic polymorphism and haplotype of the two genes were analyzed.

RESULTS

A total of 422 blood samples were used for this study, of which 224 were analyzed at PvMSP-3α and 126 at PvMSP-3β. Samples mainly were from young adults aged 18-45 years, although local patients were significantly younger than migrant laborers crossing the border at Tengchong (P < 0.0001). Molecular evolutionary analysis revealed that PvMSP-3α and PvMSP-3β underwent diversifying natural selection, and intragenic recombination contributed to the diversity of the isolates. Based on the length of the genes, we identified three types of PvMSP-3α [1.9-2.0 kb (Type-A), 1.4-1.5 kb (Type-B), and 1.1-1.3 kb (Type-C)] and two types of PvMSP-3β [1.7-2.2 kb (Type-A) and 1.4-1.5 kb (Type-B)]. Migrant laborers returning to China through Tengchong bore P. vivax infections displaying significantly higher genetic diversity than local residents.

CONCLUSIONS

Both PvMSP-3 paralogs were subjected to diversifying selection in each sample population. Clustering of alleles supports ephemeral endemic differentiation of alleles, but the broader phylogeny suggests that alleles transit the globe, perhaps accelerated by movements of migrants such as those transiting Tengchong.

Multidisciplinary Investigations of Sustained Malaria Transmission in the Greater Mekong Subregion

In the course of malaria elimination in the Greater Mekong Subregion (GMS), malaria epidemiology has experienced drastic spatiotemporal changes with residual transmission concentrated along international borders and the rising predominance of Plasmodium vivax. The emergence of Plasmodium falciparum parasites resistant to artemisinin and partner drugs renders artemisinin-based combination therapies less effective while the potential spread of multidrug-resistant parasites elicits concern. Vector behavioral changes and insecticide resistance have reduced the effectiveness of core vector control measures. In recognition of these problems, the Southeast Asian International Center of Excellence for Malaria Research (ICEMR) has been conducting multidisciplinary research to determine how human migration, antimalarial drug resistance, vector behavior, and insecticide resistance sustain malaria transmission at international borders. These efforts allow us to comprehensively understand the ecology of border malaria transmission and develop population genomics tools to identify and track parasite introduction. In addition to employing in vivo, in vitro, and molecular approaches to monitor the emergence and spread of drug-resistant parasites, we also use genomic and genetic methods to reveal novel mechanisms of antimalarial drug resistance of parasites. We also use omics and population genetics approaches to study insecticide resistance in malaria vectors and identify changes in mosquito community structure, vectorial potential, and seasonal dynamics. Collectively, the scientific findings from the ICEMR research activities offer a systematic view of the factors sustaining residual malaria transmission and identify potential solutions to these problems to accelerate malaria elimination in the GMS.

Malaria Research for Tailored Control and Elimination Strategies in the Greater Mekong Subregion

The malaria landscape in the Greater Mekong Subregion has experienced drastic changes with the ramp-up of the control efforts, revealing formidable challenges that slowed down the progress toward malaria elimination. Problems such as border malaria and cross-border malaria introduction, multidrug resistance in Plasmodium falciparum, the persistence of Plasmodium vivax, the asymptomatic parasite reservoirs, and insecticide resistance in primary vectors require integrated strategies tailored for individual nations in the region. In recognition of these challenges and the need for research, the Southeast Asian International Center of Excellence for Malaria Research has established a network of researchers and stakeholders and conducted basic and translational research to identify existing and emerging problems and develop new countermeasures. The installation of a comprehensive disease and vector surveillance system at sentinel sites in border areas with the implementation of passive/active case detection and cross-sectional surveys allowed timely detection and management of malaria cases, provided updated knowledge for effective vector control measures, and facilitated the efficacy studies of antimalarials. Incorporating sensitive molecular diagnosis to expose the significance of asymptomatic parasite reservoirs for sustaining transmission helped establish the necessary evidence to guide targeted control to eliminate residual transmission. In addition, this program has developed point-of-care diagnostics to monitor the quality of artemisinin combination therapies, delivering the needed information to the drug regulatory authorities to take measures against falsified and substandard antimalarials. To accelerate malaria elimination, this program has actively engaged with stakeholders of all levels, fostered vertical and horizontal collaborations, and enabled the effective dissemination of research findings.

Enhancing Malaria Research, Surveillance, and Control in Endemic Areas of Kenya and Ethiopia

Malaria control programs in Africa encounter daunting challenges that hinder progressive steps toward elimination of the disease. These challenges include widespread insecticide resistance in mosquito vectors, increasing outdoor malaria transmission, lack of vector surveillance and control tools suitable for outdoor biting vectors, weakness in malaria surveillance, and an inadequate number of skilled healthcare personnel. Ecological and epidemiological changes induced by environmental modifications resulting from water resource development projects pose additional barriers to malaria control. Cognizant of these challenges, our International Center of Excellence for Malaria Research (ICEMR) works in close collaboration with relevant government ministries and agencies to align its research efforts with the objectives and strategies of the national malaria control and elimination programs for the benefit of local communities. Our overall goal is to assess the impact of water resource development projects, shifting agricultural practices, and vector interventions on Plasmodium falciparum and P. vivax malaria in Kenya and Ethiopia. From 2017 to date, the ICEMR has advanced knowledge of malaria epidemiology, transmission, immunology, and pathogenesis, and developed tools to enhance vector surveillance and control, improved clinical malaria surveillance and diagnostic methods, and strengthened the capacity of local healthcare providers. Research findings from the ICEMR will inform health policy and strategic planning by ministries of health in their quest to sustain malaria control and achieve elimination goals.

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.

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.

A retrospective analysis of malaria epidemiological characteristics in Yingjiang County on the China-Myanmar border

Yingjiang County, which is on the China–Myanmar border, is the main focus for malaria elimination in China. The epidemiological characteristics of malaria in Yingjiang County were analysed in a retrospective analysis. A total of 895 malaria cases were reported in Yingjiang County between 2013 and 2019. The majority of cases occurred in males (70.7%) and individuals aged 19–59 years (77.3%). Plasmodium vivax was the predominant species (96.6%). The number of indigenous cases decreased gradually and since 2017, no indigenous cases have been reported. Malaria cases were mainly distributed in the southern and southwestern areas of the county; 55.6% of the indigenous cases were reported in Nabang Township, which also had the highest risk of imported malaria. The “1–3–7” approach has been implemented effectively, with 100% of cases reported within 24 h, 88.9% cases investigated and confirmed within 3 days and 98.5% of foci responded to within 7 days. Although malaria elimination has been achieved in Yingjiang County, sustaining elimination and preventing the re-establishment of malaria require the continued strengthening of case detection, surveillance and response systems targeting the migrant population in border areas.

The changing epidemiology of Plasmodium vivax: Insights from conventional and novel surveillance tools

Sarah Auburn and co-authors discuss the unique biology and epidemiology of P. vivax and current evidence on conventional and new approaches to surveillance.

Ecological drivers of genetic connectivity for African malaria vectors Anopheles gambiae and An. arabiensis

Anopheles gambiae and An. arabiensis are major malaria vectors in sub-Saharan Africa. Knowledge of how geographical factors drive the dispersal and gene flow of malaria vectors can help in combatting insecticide resistance spread and planning new vector control interventions. Here, we used a landscape genetics approach to investigate population relatedness and genetic connectivity of An. gambiae and An. arabiensis across Kenya and determined the changes in mosquito population genetic diversity after 20 years of intensive malaria control efforts. We found a significant reduction in genetic diversity in An. gambiae, but not in An. arabiensis as compared to prior to the 20-year period in western Kenya. Significant population structure among populations was found for both species. The most important ecological driver for dispersal and gene flow of An. gambiae and An. arabiensis was tree cover and cropland, respectively. These findings highlight that human induced environmental modifications may enhance genetic connectivity of malaria vectors.

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

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

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

Microgeographic Epidemiology of Malaria Parasites in an Irrigated Area of Western Kenya by Deep Amplicon Sequencing

To improve food security, investments in irrigated agriculture are anticipated to increase throughout Africa. However, the extent to which environmental changes from water resource development will impact malaria epidemiology remains unclear. This study was designed to compare the sensitivity of molecular markers used in deep amplicon sequencing for evaluating malaria transmission intensities and to assess malaria transmission intensity at various proximities to an irrigation scheme. Compared to ama1, csp, and msp1 amplicons, cpmp required the smallest sample size to detect differences in infection complexity between transmission risk zones. Transmission intensity was highest within 5 km of the irrigation scheme by polymerase chain reaction positivity rate, infection complexity, and linkage disequilibrium. The irrigated area provided a source of parasite infections for the surrounding 2- to 10-km area. This study highlights the suitability of the cpmp amplicon as a measure for transmission intensities and the impact of irrigation on microgeographic epidemiology of malaria parasites.

New Plasmodium vivax Genomes From the China-Myanmar Border

Plasmodium vivax is increasingly the dominant species of malaria in the Greater Mekong Subregion (GMS), which is pursuing regional malaria elimination. P. vivax lineages in the GMS are poorly characterized. Currently, P. vivax reference genomes are scarce due to difficulties in culturing the parasite and lack of high-quality samples. In addition, P. vivax is incredibly diverse, necessitating the procurement of reference genomes from different geographical regions. Here we present four new P. vivax draft genomes assembled de novo from clinical samples collected in the China-Myanmar border area. We demonstrate comparable length and content to existing genomes, with the majority of structural variation occurring around subtelomeric regions and exported proteins, which we corroborated with detection of copy number variations in these regions. We predicted peptides from all PIR gene subfamilies, except for PIR D. We confirmed that proteins classically labeled as PIR D family members are not identifiable by PIR motifs, and actually bear stronger resemblance to DUF (domain of unknown function) family DUF3671, potentially pointing to a new, closely related gene family. Further, phylogenetic analyses of MSP7 genes showed high variability within the MSP7-B family compared to MSP7-A and -C families, and the result was comparable to that from whole genome analyses. The new genome assemblies serve as a resource for studying P. vivax within the GMS.

Population genomics identifies a distinct Plasmodium vivax population on the China-Myanmar border of Southeast Asia

Plasmodium vivax has become the predominant malaria parasite and a major challenge for malaria elimination in the Greater Mekong Subregion (GMS). Yet, our knowledge about the evolution of P. vivax populations in the GMS is fragmental. We performed whole genome sequencing on 23 P. vivax samples from the China-Myanmar border (CMB) and used 21 high-coverage samples to compare to over 200 samples from the rest of the GMS. Using genome-wide single nucleotide polymorphisms (SNPs), we analyzed population differentiation, genetic structure, migration and potential selection using an array of methods. The CMB parasites displayed a higher proportion of monoclonal infections, and 52% shared over 90% of their genomes in identity-by-descent segments with at least one other sample from the CMB, suggesting preferential expansion of certain parasite strains in this region, likely resulting from the P. vivax outbreaks occurring during this study period. Principal component, admixture, fixation index and phylogenetic analyses all identified that parasites from the CMB were genetically distinct from parasites from eastern parts of the GMS (Cambodia, Laos, Vietnam, and Thailand), whereas the eastern GMS parasite populations were largely undifferentiated. Such a genetic differentiation pattern of the P. vivax populations from the GMS parasite was largely explainable through geographic distance. Using the genome-wide SNPs, we narrowed down to a set of 36 SNPs for differentiating parasites from different areas of the GMS. Genome-wide scans to determine selection in the genome with two statistical methods identified genes potentially under drug selection, including genes associated with antifolate resistance and genes linked to chloroquine resistance in Plasmodium falciparum.

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.

Molecular evidence of sustained urban malaria transmission in Amazonian Brazil, 2014-2015

The relative contribution of imported vs. locally acquired infections to urban malaria burden remains largely unexplored in Latin America, the most urbanised region in the developing world. Here we use a simple molecular epidemiology framework to examine the transmission dynamics of Plasmodium vivax in Mâncio Lima, the Amazonian municipality with the highest malaria incidence rate in Brazil. We prospectively genotyped 177 P. vivax infections diagnosed in urban residents between June 2014 and July 2015 and showed that local parasites are structured into several lineages of closely related microsatellite haplotypes, with the largest genetic cluster comprising 32% of all infections. These findings are very unlikely under the hypothesis of multiple independent imports of parasite strains from the rural surroundings. Instead, the presence of an endemic near-clonal parasite lineage circulating over 13 consecutive months is consistent with a local P. vivax transmission chain in the town, with major implications for malaria elimination efforts in this and similar urban environments across the Amazon.

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.

Genomic Epidemiology in Filarial Nematodes: Transforming the Basis for Elimination Program Decisions

Onchocerciasis and lymphatic filariasis are targeted for elimination, primarily using mass drug administration at the country and community levels. Elimination of transmission is the onchocerciasis target and global elimination as a public health problem is the end point for lymphatic filariasis. Where program duration, treatment coverage, and compliance are sufficiently high, elimination is achievable for both parasites within defined geographic areas. However, transmission has re-emerged after apparent elimination in some areas, and in others has continued despite years of mass drug treatment. A critical question is whether this re-emergence and/or persistence of transmission is due to persistence of local parasites-i.e., the result of insufficient duration or drug coverage, poor parasite response to the drugs, or inadequate methods of assessment and/or criteria for determining when to stop treatment-or due to re-introduction of parasites via human or vector movement from another endemic area. We review recent genetics-based research exploring these questions in Onchocerca volvulus, the filarial nematode that causes onchocerciasis, and Wuchereria bancrofti, the major pathogen for lymphatic filariasis. We focus in particular on the combination of genomic epidemiology and genome-wide associations to delineate transmission zones and distinguish between local and introduced parasites as the source of resurgence or continuing transmission, and to identify genetic markers associated with parasite response to chemotherapy. Our ultimate goal is to assist elimination efforts by developing easy-to-use tools that incorporate genetic information about transmission and drug response for more effective mass drug distribution, surveillance strategies, and decisions on when to stop interventions to improve sustainability of elimination.

High Plasmodium falciparum genetic diversity and temporal stability despite control efforts in high transmission settings along the international border between Zambia and the Democratic Republic of the Congo

Background

While the utility of parasite genotyping for malaria elimination has been extensively documented in low to moderate transmission settings, it has been less well-characterized in holoendemic regions. High malaria burden settings have received renewed attention acknowledging their critical role in malaria elimination. Defining the role for parasite genomics in driving these high burden settings towards elimination will enhance future control programme planning.

Methods

Amplicon deep sequencing was used to characterize parasite population genetic diversity at polymorphic Plasmodium falciparum loci, Pfama1 and Pfcsp, at two timepoints in June–July 2016 and January–March 2017 in a high transmission region along the international border between Luapula Province, Zambia and Haut-Katanga Province, the Democratic Republic of the Congo (DRC).

Results

High genetic diversity was observed across both seasons and in both countries. No evidence of population structure was observed between parasite populations on either side of the border, suggesting that this region may be one contiguous transmission zone. Despite a decline in parasite prevalence at the sampling locations in Haut-Katanga Province, no genetic signatures of a population bottleneck were detected, suggesting that larger declines in transmission may be required to reduce parasite genetic diversity. Analysing rare variants may be a suitable alternative approach for detecting epidemiologically important genetic signatures in highly diverse populations; however, the challenge is distinguishing true signals from potential artifacts introduced by small sample sizes.

Conclusions

Continuing to explore and document the utility of various parasite genotyping approaches for understanding malaria transmission in holoendemic settings will be valuable to future control and elimination programmes, empowering evidence-based selection of tools and methods to address pertinent questions, thus enabling more efficient resource allocation.

Increasing trends of malaria in a border area of the Greater Mekong Subregion

Background

Intensive malaria transmission along international borders is a significant impediment to malaria elimination in the Greater Mekong Subregion (GMS) of Southeast Asia. Passive case detection (PCD) was used to study the dynamics and trends of malaria transmission at the China–Myanmar border to provide epidemiologic information for improved malaria control.

Methods

PCD was conducted in one hospital and 12 clinics near the Laiza town in northeast Myanmar from 2011 to 2016. Clinical malaria was diagnosed by microscopy and demographic information was captured using a structured questionnaire at the time of the patient’s presentation for care.

Results

Over the study period, 6175 (19.7%) malaria cases were confirmed by microscopy from 31,326 suspected cases. The four human malaria parasite species were all identified, with Plasmodium vivax and Plasmodium falciparum accounting for 5607 (90.8%) and 481 (7.8%) of the confirmed cases, respectively. In contrast to the steady decline of malaria in the general GMS, the study site had an upward trend of malaria incidence with vivax malaria outbreaks in 2013 and 2016. Adult males, children under the age of 15, and those with occupations such as farming, being a soldier or student, had significantly higher risks of clinical malaria compared to having fevers from other aetiologies. A self-reported history of clinical malaria was also associated with a higher risk of confirmed malaria.

Conclusions

The China–Myanmar border area has experienced an overall upward trend of malaria incidence in recent years with P. vivax becoming the predominant species. Evidence-based control strategies need to focus on high-risk populations.

Limited differentiation among Plasmodium vivax populations from the northwest and to the south Pacific Coast of Colombia: A malaria corridor?

Background

Malaria remains endemic in several countries of South America with low to moderate transmission intensity. Regional human migration through underserved endemic areas may be responsible for significant parasite dispersion making the disease resilient to interventions. Thus, the genetic characterization of malarial parasites is an important tool to assess how endemic areas may connect via the movement of infected individuals. Here, four sites in geographically separated areas reporting 80% of the malaria morbidity in Colombia were studied. The sites are located on an imaginary transect line of 1,500 km from the northwest to the south Pacific Coast of Colombia with a minimal distance of 500 km between populations that display noticeable ethnic, economic, epidemiological, and ecological differences.

Methodology/Principal findings

A total of 624 Plasmodium vivax samples from the four populations were genotyped by using eight microsatellite loci. Although a strong geographic structure was expected between these populations, only moderate evidence of genetic differentiation was observed using a suite of population genetic analyses. High genetic diversity, shared alleles, and low linkage disequilibrium were also found in these P. vivax populations providing no evidence for a bottleneck or clonal expansions as expected from recent reductions in the transmission that could have been the result of scaling up interventions or environmental changes. These patterns are consistent with a disease that is not only endemic in each site but also imply that there is gene flow among these populations across 1,500 km.

Conclusion /Significance

The observed patterns in P. vivax are consistent with a “corridor” where connected endemic areas can sustain a high level of genetic diversity locally and can restore parasite-subdivided populations via migration of infected individuals even after local interventions achieved a substantial reduction of clinical cases. The consequences of these findings in terms of control and elimination are discussed.

The regional movements of infected individuals that connect suitable transmission areas make malaria resilient to control efforts. Those movements are expected to leave genetic signatures in the parasite populations that can be detected using analytical tools. In this study, the genetic makeups of Plasmodium vivax populations were characterized to assess whether the most endemic areas in Colombia were connected. Samples were collected from passive surveillance studies in four locations across an imaginary transect line of 1,500 km from the northwest to the south Pacific Coast of Colombia (South America). Considering the distance, and contrary to expectations, we found weak levels of genetic differentiation between these parasite populations with no evidence indicating that their genetic diversity has been eroded as expected whenever the prevalence of the disease is successfully reduced, e.g., through control programs or environmental changes. Although the sampling lacks the geographic and temporal detail to describe how the dispersion of parasite lineages occurred, the observed patterns are consistent with a series of infected populations that are connected in space by human movements allowing the parasite to diffuse across this 1,500 km transect. This malaria corridor needs to be characterized to achieve elimination.

The Plasmodium knowlesi Pk41 surface protein diversity, natural selection, sub population and geographical clustering: a 6-cysteine protein family member

Introduction

The zoonotic malaria parasite Plasmodium knowlesi has currently become the most dominant form of infection in humans in Malaysia and is an emerging infectious disease in most Southeast Asian countries. The P41 is a merozoite surface protein belonging to the 6-cysteine family and is a well-characterized vaccine candidate in P. vivax and P. falciparum; however, no study has been done in the orthologous gene of P. knowlesi. This study investigates the level of polymorphism, haplotypes and natural selection of pk41 genes in clinical isolates from Malaysia.

Method

Thirty-five full-length pk41 sequences from clinical isolates of Malaysia along with four laboratory lines (along with H-strain) were downloaded from public databases. For comparative analysis between species, orthologous P41 genes from P. falciparum, P. vivax, P. coatneyi and P. cynomolgi were also downloaded. Genetic diversity, polymorphism, haplotype and natural selection were determined using DnaSP 5.10 software. Phylogenetic relationships between Pk41 genes were determined using MEGA 5.0 software.

Results

Analysis of 39 full-length pk41 sequences along with the H-strain identified 36 SNPs (20 non-synonymous and 16 synonymous substitutions) resulting in 31 haplotypes. Nucleotide diversity across the full-length gene was low and was similar to its ortholog in P. vivax; pv41. Domain-wise amino acid analysis of the two s48/45 domains indicated low level of polymorphisms for both the domains, and the glutamic acid rich region had extensive size variations. In the central domain, upstream to the glutamate rich region, a unique two to six (K-E)_n repeat region was identified within the clinical isolates. Overall, the pk41 genes were indicative of negative/purifying selection due to functional constraints. Domain-wise analysis of the s48/45 domains also indicated purifying selection. However, analysis of Tajima’s D across the genes identified non-synonymous SNPs in the s48/45 domain II with high positive values indicating possible epitope binding regions. All the 6-cysteine residues within the s48/45 domains were conserved within the clinical isolates indicating functional conservation of these regions. Phylogenetic analysis of full-length pk41 genes indicated geographical clustering and identified three subpopulations of P. knowlesi; one originating in the laboratory lines and two originating from Sarawak, Malaysian Borneo.

Conclusion

This is the first study to report on the polymorphism and natural selection of pk41 genes from clinical isolates of Malaysia. The results reveal that there is low level of polymorphism in both s48/45 domains, indicating that this antigen could be a potential vaccine target. However, genetic and molecular immunology studies involving higher number of samples from various parts of Malaysia would be necessary to validate this antigen’s candidacy as a vaccine target for P. knowlesi.

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.

Chloroquine efficacy for Plasmodium vivax in Myanmar in populations with high genetic diversity and moderate parasite gene flow

BACKGROUND

Plasmodium vivax malaria remains a major public health burden in Myanmar. Resistance to chloroquine (CQ), the first-line treatment for P. vivax, has been reported in the country and has potential to undermine local control efforts.

METHODS

Patients over 6 years of age with uncomplicated P. vivax mono-infection were enrolled into clinical efficacy studies in Myawaddy in 2014 and Kawthoung in 2012. Study participants received a standard dose of CQ (25 mg/kg over 3 days) followed by weekly review until day 28. Pvmdr1 copy number (CN) and microsatellite diversity were assessed on samples from the patients enrolled in the clinical study and additional cross-sectional surveys undertaken in Myawaddy and Shwegyin in 2012.

RESULTS

A total of 85 patients were enrolled in the CQ clinical studies, 25 in Myawaddy and 60 in Kawthoung. One patient in Myawaddy (1.2%) had an early treatment failure and two patients (2.3%) in Kawthoung presented with late treatment failures on day 28. The day 28 efficacy was 92.0% (95% CI 71.6-97.9) in Myawaddy and 98.3% (95% CI 88.7-99.8) in Kawthoung. By day 2, 92.2% (23/25) in Myawaddy and 85.0% (51/60) in Kawthoung were aparasitaemic. Genotyping and pvmdr1 CN assessment was undertaken on 43, 52 and 46 clinical isolates from Myawaddy, Kawthoung and Shwegyin respectively. Pvmdr1 amplification was observed in 3.2% (1/31) of isolates in Myawaddy, 0% (0/49) in Kawthoung and 2.5% (1/40) in Shwegyin. Diversity was high in all sites (H _E 0.855-0.876), with low inter-population differentiation (F _ST 0.016-0.026, P < 0.05).

CONCLUSIONS

Treatment failures after chloroquine were observed following chloroquine monotherapy, with pvmdr1 amplification present in both Myawaddy and Shwegyin. The results emphasize the importance of ongoing P. vivax drug resistance surveillance in Myanmar, particularly given the potential connectivity between parasite population at different sites.