Higher microsatellite diversity in Plasmodium vivax than in sympatric Plasmodium falciparum populations in Pursat, Western Cambodia

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.

Vivax malaria: a possible stumbling block for malaria elimination in India

Plasmodium vivax is geographically the most widely dispersed human malaria parasite species. It has shown resilience and a great deal of adaptability. Genomic studies suggest that P. vivax originated from Asia or Africa and moved to the rest of the world. Although P. vivax is evolutionarily an older species than Plasmodium falciparum, its biology, transmission, pathology, and control still require better elucidation. P. vivax poses problems for malaria elimination because of the ability of a single primary infection to produce multiple relapses over months and years. P. vivax malaria elimination program needs early diagnosis, and prompt and complete radical treatment, which is challenging, to simultaneously exterminate the circulating parasites and dormant hypnozoites lodged in the hepatocytes of the host liver. As prompt surveillance and effective treatments are rolled out, preventing primaquine toxicity in the patients having glucose-6-phosphate dehydrogenase (G6PD) deficiency should be a priority for the vivax elimination program. This review sheds light on the burden of P. vivax, changing epidemiological patterns, the hurdles in elimination efforts, and the essential tools needed not just in India but globally. These tools encompass innovative treatments for eliminating dormant parasites, coping with evolving drug resistance, and the development of potential vaccines against the parasite.

Plasmodium vivax transcriptomics: What is new?

Malaria is the leading human parasitosis and is transmitted through the bite of anopheline mosquitoes infected with parasites of the genus Plasmodium spp. Among the seven species that cause malaria in humans, Plasmodium vivax is the most prevalent species in Latin America. In recent years, there have been an increasing number of reports of clinical complications caused by P. vivax infections, which were previously neglected and underestimated. P. vivax biology remains with large gaps. The emergence of next-generation sequencing technology has ensured a breakthrough in species knowledge. Coupled with this, the deposition of the P. vivax Sal-1 reference genome allowed an increase in transcriptomics projects by accessing messenger RNA. Thus, the regulation of differential gene expression according to the parasite life stage was verified, and several expressed genes were linked to different biological functions. Today, with the progress associated with RNA sequencing technologies, it is possible to detect nuances and obtain robust results. Discoveries provided by transcriptomic studies allow us to understand topics such as RNA expression and regulation and proteins and metabolic pathways involved during different stages of the parasite life cycle. The information obtained enables a better comprehension of immune system evasion mechanisms; invasion and adhesion strategies used by the parasite; as well as new vaccine targets, potential molecular markers, and others therapeutic targets. In this review, we provide new insights into P. vivax biology by summarizing recent findings in transcriptomic studies.

Plasmodium vivax From Duffy-Negative and Duffy-Positive Individuals Share Similar Gene Pools in East Africa

Plasmodium vivax malaria was thought to be rare in Africa, but an increasing number of P. vivax cases reported across Africa and in Duffy-negative individuals challenges this dogma. The genetic characteristics of P. vivax in Duffy-negative infections, the transmission of P. vivax in East Africa, and the impact of environments on transmission remain largely unknown. This study examined genetic and transmission features of P. vivax from 107 Duffy-negative and 305 Duffy-positive individuals in Ethiopia and Sudan. No clear genetic differentiation was found in P. vivax between the 2 Duffy groups, indicating between-host transmission. P. vivax from Ethiopia and Sudan showed similar genetic clusters, except samples from Khartoum, possibly due to distance and road density that inhibited parasite gene flow. This study is the first to show that P. vivax can transmit to and from Duffy-negative individuals and provides critical insights into the spread of P. vivax in sub-Saharan Africa.

Malaria in the 'Omics Era'

Genomics has revolutionised the study of the biology of parasitic diseases. The first Eukaryotic parasite to have its genome sequenced was the malaria parasite Plasmodium falciparum. Since then, Plasmodium genomics has continued to lead the way in the study of the genome biology of parasites, both in breadth-the number of Plasmodium species' genomes sequenced-and in depth-massive-scale genome re-sequencing of several key species. Here, we review some of the insights into the biology, evolution and population genetics of Plasmodium gained from genome sequencing, and look at potential new avenues in the future genome-scale study of its biology.

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.

Population genomic evidence of Plasmodium vivax Southeast Asian origin

Plasmodium vivax is the most common and widespread human malaria parasite. It was recently proposed that P. vivax originates from sub-Saharan Africa based on the circulation of its closest genetic relatives (P. vivax-like) among African great apes. However, the limited number of genetic markers and samples investigated questions the robustness of this hypothesis. Here, we extensively characterized the genomic variations of 447 human P. vivax strains and 19 ape P. vivax-like strains collected worldwide. Phylogenetic relationships between human and ape Plasmodium strains revealed that P. vivax is a sister clade of P. vivax-like, not included within the radiation of P. vivax-like By investigating various aspects of P. vivax genetic variation, we identified several notable geographical patterns in summary statistics in function of the increasing geographic distance from Southeast Asia, suggesting that P. vivax may have derived from a single area in Asia through serial founder effects.

A suitable RNA preparation methodology for whole transcriptome shotgun sequencing harvested from Plasmodium vivax-infected patients

Plasmodium vivax is a world-threatening human malaria parasite, whose biology remains elusive. The unavailability of in vitro culture, and the difficulties in getting a high number of pure parasites makes RNA isolation in quantity and quality a challenge. Here, a methodological outline for RNA-seq from P. vivax isolates with low parasitemia is presented, combining parasite maturation and enrichment with efficient RNA extraction, yielding ~ 100 pg.µL⁻¹ of RNA, suitable for SMART-Seq Ultra-Low Input RNA library and Illumina sequencing. Unbiased coding transcriptome of ~ 4 M reads was achieved for four patient isolates with ~ 51% of transcripts mapped to the P. vivax P01 reference genome, presenting heterogeneous profiles of expression among individual isolates. Amongst the most transcribed genes in all isolates, a parasite-staged mixed repertoire of conserved parasite metabolic, membrane and exported proteins was observed. Still, a quarter of transcribed genes remain functionally uncharacterized. In parallel, a P. falciparum Brazilian isolate was also analyzed and 57% of its transcripts mapped against IT genome. Comparison of transcriptomes of the two species revealed a common trophozoite-staged expression profile, with several homologous genes being expressed. Collectively, these results will positively impact vivax research improving knowledge of P. vivax biology.

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.

Genome-wide microsatellite characteristics of five human Plasmodium species, focusing on Plasmodium malariae and P. ovale curtisi

Microsatellites can be utilized to explore genotypes, population structure, and other genomic features of eukaryotes. Systematic characterization of microsatellites has not been a focus for several species of Plasmodium, including P. malariae and P. ovale, as the majority of malaria elimination programs are focused on P. falciparum and to a lesser extent P. vivax. Here, five human malaria species (P. falciparum, P. vivax, P. malariae, P. ovale curtisi, and P. knowlesi) were investigated with the aim of conducting in-depth categorization of microsatellites for P. malariae and P. ovale curtisi. Investigation of reference genomes for microsatellites with unit motifs of 1-10 base pairs indicates high diversity among the five Plasmodium species. Plasmodium malariae, with the largest genome size, displays the second highest microsatellite density (1421 No./Mbp; 5% coverage) next to P. falciparum (3634 No./Mbp; 12% coverage). The lowest microsatellite density was observed in P. vivax (773 No./Mbp; 2% coverage). A, AT, and AAT are the most commonly repeated motifs in the Plasmodium species. For P. malariae and P. ovale curtisi, microsatellite-related sequences are observed in approximately 18-29% of coding sequences (CDS). Lysine, asparagine, and glutamic acids are most frequently coded by microsatellite-related CDS. The majority of these CDS could be related to the gene ontology terms "cell parts," "binding," "developmental processes," and "metabolic processes." The present study provides a comprehensive overview of microsatellite distribution and can assist in the planning and development of potentially useful genetic tools for further investigation of P. malariae and P. ovale curtisi epidemiology.

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.

Malaria Molecular Epidemiology: An Evolutionary Genetics Perspective

Malaria is a vector-borne disease that involves multiple parasite species in a variety of ecological settings. However, the parasite species causing the disease, the prevalence of subclinical infections, the emergence of drug resistance, the scale-up of interventions, and the ecological factors affecting malaria transmission, among others, are aspects that vary across areas where malaria is endemic. Such complexities have propelled the study of parasite genetic diversity patterns in the context of epidemiologic investigations. Importantly, molecular studies indicate that the time and spatial distribution of malaria cases reflect epidemiologic processes that cannot be fully understood without characterizing the evolutionary forces shaping parasite population genetic patterns. Although broad in scope, this review in the Microbiology Spectrum Curated Collection: Advances in Molecular Epidemiology highlights the need for understanding population genetic concepts when interpreting parasite molecular data. First, we discuss malaria complexity in terms of the parasite species involved. Second, we describe how molecular data are changing our understanding of malaria incidence and infectiousness. Third, we compare different approaches to generate parasite genetic information in the context of epidemiologically relevant questions related to malaria control. Finally, we describe a few Plasmodium genomic studies as evidence of how these approaches will provide new insights into the malaria disease dynamics. *This article is part of a curated collection.

Polymorphic markers for identification of parasite population in Plasmodium malariae

Background

Molecular genotyping in Plasmodium serves many aims including providing tools for studying parasite population genetics and distinguishing recrudescence from reinfection. Microsatellite typing, insertion-deletion (INDEL) and single nucleotide polymorphisms is used for genotyping, but only limited information is available for Plasmodium malariae, an important human malaria species. This study aimed to provide a set of genetic markers to facilitate the study of P. malariae population genetics.

Methods

Markers for microsatellite genotyping and pmmsp1 gene polymorphisms were developed and validated in symptomatic P. malariae field isolates from Myanmar (N = 37). Fragment analysis was used to determine allele sizes at each locus to calculate multiplicity of infections (MOI), linkage disequilibrium, heterozygosity and construct dendrograms. Nucleotide diversity (π), number of haplotypes, and genetic diversity (H_d) were assessed and a phylogenetic tree was constructed. Genome-wide microsatellite maps with annotated regions of newly identified markers were constructed.

Results

Six microsatellite markers were developed and tested in 37 P. malariae isolates which showed sufficient heterozygosity (0.530–0.922), and absence of linkage disequilibrium (I_A^S=0.03, p value > 0.05) (N = 37). In addition, a tandem repeat (VNTR)-based pmmsp1 INDEL polymorphisms marker was developed and assessed in 27 P. malariae isolates showing a nucleotide diversity of 0.0976, haplotype gene diversity of 0.698 and identified 14 unique variants. The size of VNTR consensus repeat unit adopted as allele was 27 base pairs. The markers Pm12_426 and pmmsp1 showed greatest diversity with heterozygosity scores of 0.920 and 0.835, respectively. Using six microsatellites markers, the likelihood that any two parasite strains would have the same microsatellite genotypes was 8.46 × 10⁻⁴ and was further reduced to 1.66 × 10⁻⁴ when pmmsp1 polymorphisms were included.

Conclusions

Six novel microsatellites genotyping markers and a set of pmmsp1 VNTR-based INDEL polymorphisms markers for P. malariae were developed and validated. Each marker could be independently or in combination employed to access genotyping of the parasite. The newly developed markers may serve as a useful tool for investigating parasite diversity, population genetics, molecular epidemiology and for distinguishing recrudescence from reinfection in drug efficacy studies.

Genetic Diversity of Plasmodium vivax in Clinical Isolates from Southern Thailand using PvMSP1, PvMSP3 (PvMSP3α, PvMSP3β) Genes and Eight Microsatellite Markers

Plasmodium vivax is usually considered morbidity in endemic areas of Asia, Central and South America, and some part of Africa. In Thailand, previous studies indicated the genetic diversity of P. vivax in malaria-endemic regions such as the western part of Thailand bordering with Myanmar. The objective of the study is to investigate the genetic diversity of P. vivax circulating in Southern Thailand by using 3 antigenic markers and 8 microsatellite markers. Dried blood spots were collected from Chumphon, Phang Nga, Ranong and, Surat Thani provinces of Thailand. By PCR, 3 distinct sizes of PvMSP3α, 2 sizes of PvMSP3β and 2 sizes of PvMSP1 F2 were detected based on the length of PCR products, respectively. PCR/RFLP analyses of these antigen genes revealed high levels of genetic diversity. The genotyping of 8 microsatellite loci showed high genetic diversity as indicated by high alleles per locus and high expected heterozygosity (HE). The genotyping markers also showed multiple-clones of infection. Mixed genotypes were detected in 4.8% of PvMSP3α, 29.1% in PvMSP3β and 55.3% of microsatellite markers. These results showed that there was high genetic diversity of P. vivax isolated from Southern Thailand, indicating that the genetic diversity of P. vivax in this region was comparable to those observed other areas of Thailand.

Structural basis for neutralization of Plasmodium vivax by naturally acquired human antibodies that target DBP

The Plasmodium vivax Duffy-binding protein (DBP) is a prime target of the protective immune response and a promising vaccine candidate for P. vivax malaria. Naturally acquired immunity (NAI) protects against malaria in adults residing in infection-endemic regions, and the passive transfer of malarial immunity confers protection. A vaccine that replicates NAI will effectively prevent disease. Here, we report the structures of DBP region II in complex with human-derived, neutralizing monoclonal antibodies obtained from an individual in a malaria-endemic area with NAI. We identified protective epitopes using X-ray crystallography, hydrogen-deuterium exchange mass spectrometry, mutational mapping and P. vivax invasion studies. These approaches reveal that naturally acquired human antibodies neutralize P. vivax by targeting the binding site for Duffy antigen receptor for chemokines (DARC) and the dimer interface of P. vivax DBP. Antibody binding is unaffected by polymorphisms in the vicinity of epitopes, suggesting that the antibodies have evolved to engage multiple polymorphic variants of DBP. The human antibody epitopes are broadly conserved and are distinct from previously defined epitopes for broadly conserved murine monoclonal antibodies. A library of globally conserved epitopes of neutralizing human antibodies offers possibilities for rational design of strain-transcending DBP-based vaccines and therapeutics against P. vivax.

Identification and Characterization of Functional Human Monoclonal Antibodies to Plasmodium vivax Duffy-Binding Protein

Plasmodium vivax invasion of reticulocytes relies on distinct receptor-ligand interactions between the parasite and host erythrocytes. Engagement of the highly polymorphic domain II of the P. vivax Duffy-binding protein (DBPII) with the erythrocyte's Duffy Ag receptor for chemokines (DARC) is essential. Some P. vivax-exposed individuals acquired Abs to DBPII that block DBPII-DARC interaction and inhibit P. vivax reticulocyte invasion, and Ab levels correlate with protection against P. vivax malaria. To better understand the functional characteristics and fine specificity of protective human Abs to DBPII, we sorted single DBPII-specific IgG⁺ memory B cells from three individuals with high blocking activity to DBPII. We identified 12 DBPII-specific human mAbs from distinct lineages that blocked DBPII-DARC binding. All mAbs were P. vivax strain transcending and targeted known binding motifs of DBPII with DARC. Eleven mAbs competed with each other for binding, indicating recognition of the same or overlapping epitopes. Naturally acquired blocking Abs to DBPII from individuals with high levels residing in different P. vivax-endemic areas worldwide competed with mAbs, suggesting broadly shared recognition sites. We also found that mAbs inhibited P. vivax entry into reticulocytes in vitro. These findings suggest that IgG⁺ memory B cell activity in individuals with P. vivax strain-transcending Abs to DBPII display a limited clonal response with inhibitory blocking directed against a distinct region of the molecule.

Paleopathological Considerations on Malaria Infection in Korea before the 20th Century

Malaria, one of the deadliest diseases in human history, still infects many people worldwide. Among the species of the genus Plasmodium, P. vivax is commonly found in temperate-zone countries including South Korea. In this article, we first review the history of malarial infection in Korea by means of studies on Joseon documents and the related scientific data on the evolutionary history of P. vivax in Asia. According to the historical records, malarial infection was not unusual in pre-20th-century Korean society. We also found that certain behaviors of the Joseon people might have affected the host-vector-pathogen relationship, which could explain why malarial infection prevalence was so high in Korea at that time. In our review of genetic studies on P. vivax, we identified substantial geographic differentiation among continents and even between neighboring countries. Based on these, we were able to formulate a strategy for future analysis of ancient Plasmodium strains in Korea.

Plasmodium vivax Biology: Insights Provided by Genomics, Transcriptomics and Proteomics

During the last decade, the vast omics field has revolutionized biological research, especially the genomics, transcriptomics and proteomics branches, as technological tools become available to the field researcher and allow difficult question-driven studies to be addressed. Parasitology has greatly benefited from next generation sequencing (NGS) projects, which have resulted in a broadened comprehension of basic parasite molecular biology, ecology and epidemiology. Malariology is one example where application of this technology has greatly contributed to a better understanding of Plasmodium spp. biology and host-parasite interactions. Among the several parasite species that cause human malaria, the neglected Plasmodium vivax presents great research challenges, as in vitro culturing is not yet feasible and functional assays are heavily limited. Therefore, there are gaps in our P. vivax biology knowledge that affect decisions for control policies aiming to eradicate vivax malaria in the near future. In this review, we provide a snapshot of key discoveries already achieved in P. vivax sequencing projects, focusing on developments, hurdles, and limitations currently faced by the research community, as well as perspectives on future vivax malaria research.

Increasingly inbred and fragmented populations of Plasmodium vivax associated with the eastward decline in malaria transmission across the Southwest Pacific

The human malaria parasite Plasmodium vivax is more resistant to malaria control strategies than Plasmodium falciparum, and maintains high genetic diversity even when transmission is low. To investigate whether declining P. vivax transmission leads to increasing population structure that would facilitate elimination, we genotyped samples from across the Southwest Pacific region, which experiences an eastward decline in malaria transmission, as well as samples from two time points at one site (Tetere, Solomon Islands) during intensified malaria control. Analysis of 887 P. vivax microsatellite haplotypes from hyperendemic Papua New Guinea (PNG, n = 443), meso-hyperendemic Solomon Islands (n = 420), and hypoendemic Vanuatu (n = 24) revealed increasing population structure and multilocus linkage disequilibrium yet a modest decline in diversity as transmission decreases over space and time. In Solomon Islands, which has had sustained control efforts for 20 years, and Vanuatu, which has experienced sustained low transmission for many years, significant population structure was observed at different spatial scales. We conclude that control efforts will eventually impact P. vivax population structure and with sustained pressure, populations may eventually fragment into a limited number of clustered foci that could be targeted for elimination.

Plasmodium vivax is a major human malaria parasite, common in endemic areas outside sub-Saharan Africa, and more difficult to control than other malaria parasite species. The distinct lifecycle biology of P. vivax is thought to contribute to its more stable and efficient transmission allowing the maintenance of high diversity and potentially, gene flow. Independent studies are therefore needed to understand how P. vivax populations respond to changing transmission levels, in order to inform malaria control and elimination efforts. Here we have determined parasite population genetic structure in three countries of the Southwest Pacific, an island chain with a natural west to east decline in transmission intensity (Papua New Guinea > Solomon Islands > Vanuatu). With declining transmission, P. vivax populations experience only a modest decline in diversity but a significant increase in multilocus linkage disequilibrium and population structure, indicating that parasite populations become more inbred and begin to fragment into clustered foci. Analysis of two time points in one study area (Tetere, Solomon Islands) also show similar changes in association with intensifying malaria control. The results indicate that with long term sustained malaria control P. vivax populations will eventually fracture into population clusters that could be targeted for elimination.

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

Background

In Myanmar, civil unrest and the establishment of internally displaced person (IDP) settlements along the Myanmar-China border have impacted malaria transmission.

Methods

Microsatellite markers were used to examine source-sink dynamics for Plasmodium vivax between IDP settlements and surrounding villages in the border region. Genotypic structure and diversity were compared across the 3 years following the establishment of IDP settlements, to infer demographic history. We investigated whether human migration and landscape heterogeneity contributed to P. vivax transmission.

Results

P. vivax from IDP settlements and local communities consistently exhibited high genetic diversity within populations but low polyclonality within individuals. No apparent genetic structure was observed among populations and years. P. vivax genotypes in China were similar to those in Myanmar, and parasite introduction was unidirectional. Landscape factors, including distance, elevation, and land cover, do not appear to impede parasite gene flow.

Conclusions

The admixture of P. vivax genotypes suggested that parasite gene flow via human movement contributes to the spread of malaria both locally in Myanmar and across the international border. Our genetic findings highlight the presence of large P. vivax gene reservoirs that can sustain transmission. Thus, it is important to reinforce and improve existing control efforts along border areas.

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 Genetics of Plasmodium vivax in Four High Malaria Endemic Areas in Thailand

Recent trends of malaria in Thailand illustrate an increasing proportion of Plasmodium vivax, indicating the importance of P. vivax as a major causative agent of malaria. P. vivax malaria is usually considered a benign disease so the knowledge of this parasite has been limited, especially the genetic diversity and genetic structure of isolates from different endemic areas. The aim of this study was to examine the population genetics and structure of P. vivax isolates from 4 provinces with different malaria endemic settings in Thailand using 6 microsatellite markers. Total 234 blood samples from P. vivax mono-infected patients were collected. Strong genetic diversity was observed across all study sites; the expected heterozygosity values ranged from 0.5871 to 0.9033. Genetic variability in this study divided P. vivax population into 3 clusters; first was P. vivax isolates from Mae Hong Son and Kanchanaburi Provinces located on the western part of Thailand; second, Yala isolates from the south; and third, Chanthaburi isolates from the east. P. vivax isolates from patients having parasite clearance time (PCT) longer than 24 hr after the first dose of chloroquine treatment had higher diversity when compared with those having PCT within 24 hr. This study revealed a clear evidence of different population structure of P. vivax from different malaria endemic areas of Thailand. The findings provide beneficial information to malaria control programme as it is a useful tool to track the source of infections and current malaria control efforts.

Malaria Epidemiology at the Clone Level

Genotyping to distinguish between parasite clones is nowadays a standard in many molecular epidemiological studies of malaria. It has become crucial in drug trials and to follow individual clones in epidemiological studies, and to understand how drug resistance emerges and spreads. Here, we review the applications of the increasingly available genotyping tools and whole-genome sequencing data, and argue for a better integration of population genetics findings into malaria-control strategies.

Transmission dynamics of co-endemic Plasmodium vivax and P. falciparum in Ethiopia and prevalence of antimalarial resistant genotypes

Ethiopia is one of the few African countries where Plasmodium vivax is co-endemic with P. falciparum. Malaria transmission is seasonal and transmission intensity varies mainly by landscape and climate. Although the recent emergence of drug resistant parasites presents a major issue to malaria control in Ethiopia, little is known about the transmission pathways of parasite species and prevalence of resistant markers. This study used microsatellites to determine population diversity and gene flow patterns of P. falciparum (N = 226) and P. vivax (N = 205), as well as prevalence of drug resistant markers to infer the impact of gene flow and existing malaria treatment regimes. Plasmodium falciparum indicated a higher rate of polyclonal infections than P. vivax. Both species revealed moderate genetic diversity and similar population structure. Populations in the northern highlands were closely related to the eastern Rift Valley, but slightly distinct from the southern basin area. Gene flow via human migrations between the northern and eastern populations were frequent and mostly bidirectional. Landscape genetic analyses indicated that environmental heterogeneity and geographical distance did not constrain parasite gene flow. This may partly explain similar patterns of resistant marker prevalence. In P. falciparum, a high prevalence of mutant alleles was detected in codons related to chloroquine (pfcrt and pfmdr1) and sulfadoxine-pyrimethamine (pfdhps and pfdhfr) resistance. Over 60% of the samples showed pfmdr1 duplications. Nevertheless, no mutation was detected in pfK13 that relates to artemisinin resistance. In P. vivax, while sequences of pvcrt-o were highly conserved and less than 5% of the samples showed pvmdr duplications, over 50% of the samples had pvmdr1 976F mutation. It remains to be tested if this mutation relates to chloroquine resistance. Monitoring the extent of malaria spread and markers of drug resistance is imperative to inform policy for evidence-based antimalarial choice and interventions. To effectively reduce malaria burden in Ethiopia, control efforts should focus on seasonal migrant populations.

Genetic micro-epidemiology of malaria in Papua Indonesia: Extensive P. vivax diversity and a distinct subpopulation of asymptomatic P. falciparum infections

BACKGROUND

Genetic analyses of Plasmodium have potential to inform on transmission dynamics, but few studies have evaluated this on a local spatial scale. We used microsatellite genotyping to characterise the micro-epidemiology of P. vivax and P. falciparum diversity to inform malaria control strategies in Timika, Papua Indonesia.

METHODS

Genotyping was undertaken on 713 sympatric P. falciparum and P. vivax isolates from a cross-sectional household survey and clinical studies conducted in Timika. Standard population genetic measures were applied, and the data was compared to published data from Kalimantan, Bangka, Sumba and West Timor.

RESULTS

Higher diversity (HE = 0.847 vs 0.625; p = 0.017) and polyclonality (46.2% vs 16.5%, p<0.001) were observed in P. vivax versus P. falciparum. Distinct P. falciparum substructure was observed, with two subpopulations, K1 and K2. K1 was comprised solely of asymptomatic infections and displayed greater relatedness to isolates from Sumba than to K2, possibly reflecting imported infections.

CONCLUSIONS

The results demonstrate the greater refractoriness of P. vivax versus P. falciparum to control measures, and risk of distinct parasite subpopulations persisting in the community undetected by passive surveillance. These findings highlight the need for complimentary new surveillance strategies to identify transmission patterns that cannot be detected with traditional malariometric methods.

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.

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

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

Dissecting malaria biology and epidemiology using population genetics and genomics

Molecular approaches have an increasingly recognized utility in surveillance of malaria parasite populations, not only in defining prevalence and incidence with higher sensitivity than traditional methods, but also in monitoring local and regional parasite transmission patterns. In this review, we provide an overview of population genetic and genomic studies of human-infecting Plasmodium species, highlighting recent advances in the field. In accordance with the renewed impetus for malaria eradication, many studies are now using genetic and genomic epidemiology to support local evidence-based intervention strategies. Microsatellite genotyping remains a popular approach for both Plasmodium falciparum and Plasmodium vivax. However, with the increasing availability of whole genome sequencing data enabling effective single nucleotide polymorphism-based panels tailored to a given study question and setting, this approach is gaining popularity. The availability of new reference genomes for Plasmodium malariae and Plasmodium ovale should see a surge in similar molecular studies on these currently neglected species. Genomic studies are revealing new insights into important adaptive mechanisms of the parasite including antimalarial drug resistance. The advent of new methodologies such as selective whole genome amplification for dealing with extensive human DNA in low density field isolates should see genome-wide approaches becoming routine for parasite surveillance once the economic costs outweigh the current cost benefits of targeted approaches.

Implications of Plasmodium vivax Biology for Control, Elimination, and Research

This paper summarizes our current understanding of the biology of Plasmodium vivax, how it differs from Plasmodium falciparum, and how these differences explain the need for P. vivax-tailored interventions. The article further pinpoints knowledge gaps where investments in research are needed to help identify and develop such specific interventions. The principal obstacles to reduce and eventually eliminate P. vivax reside in 1) its higher vectorial capacity compared with P. falciparum due to its ability to develop at lower temperature and over a shorter sporogonic cycle in the vector, allowing transmission in temperate zones and making it less sensitive to vector control measures that are otherwise effective on P. falciparum; 2) the presence of dormant liver forms (hypnozoites), sustaining multiple relapsing episodes from a single infectious bite that cannot be diagnosed and are not susceptible to any available antimalarial except primaquine, with routine deployment restricted by toxicity; 3) low parasite densities, which are difficult to detect with current diagnostics leading to missed diagnoses and delayed treatments (and protracted transmission), coupled with 4) transmission stages (gametocytes) occurring early in acute infections, before infection is diagnosed.

Where chloroquine still works: the genetic make-up and susceptibility of Plasmodium vivax to chloroquine plus primaquine in Bhutan

BACKGROUND

Bhutan has made substantial progress in reducing malaria incidence. The national guidelines recommend chloroquine (CQ) and primaquine (PQ) for radical cure of uncomplicated Plasmodium vivax, but the local efficacy has not been assessed. The impact of cases imported from India on the genetic make-up of the local vivax populations is currently unknown.

METHODS

Patients over 4 years of age with uncomplicated P. vivax mono-infection were enrolled into a clinical efficacy study and molecular survey. Study participants received a standard dose of CQ (25 mg/kg over 3 days) followed by weekly review until day 28. On day 28 a 14-day regimen of PQ (0.25 mg/kg/day) was commenced under direct observation. After day 42, patients were followed up monthly for a year. The primary and secondary endpoints were risk of treatment failure at day 28 and at 1 year. Parasite genotyping was undertaken at nine tandem repeat markers, and standard population genetic metrics were applied to examine population diversity and structure in infections thought to be acquired inside or outside of Bhutan.

RESULTS

A total of 24 patients were enrolled in the clinical study between April 2013 and October 2015. Eight patients (33.3 %) were lost to follow-up in the first 6 months and another eight patients lost between 6 and 12 months. No (0/24) treatment failures occurred by day 28 and no (0/8) parasitaemia was detected following PQ treatment. Some 95.8 % (23/24) of patients were aparasitaemic by day 2. There were no haemolytic or serious events. Genotyping was undertaken on parasites from 12 autochthonous cases and 16 suspected imported cases. Diversity was high (H E 0.87 and 0.90) in both populations. There was no notable differentiation between the autochthonous and imported populations.

CONCLUSIONS

CQ and PQ remains effective for radical cure of P. vivax in Bhutan. The genetic analyses indicate that imported infections are sustaining the local vivax population, with concomitant risk of introducing drug-resistant strains.

Complexity of Infection and Genetic Diversity in Cambodian Plasmodium vivax

BACKGROUND

Plasmodium vivax is the most widely distributed human malaria parasite with 2.9 billion people living in endemic areas. Despite intensive malaria control efforts, the proportion of cases attributed to P. vivax is increasing in many countries. Genetic analyses of the parasite population and its dynamics could provide an assessment of the efficacy of control efforts, but, unfortunately, these studies are limited in P. vivax by the lack of informative markers and high-throughput genotyping methods.

METHODOLOGY/PRINCIPAL FINDINGS

We developed a sequencing-based assay to simultaneously genotype more than 100 SNPs and applied this approach to ~500 P. vivax-infected individuals recruited across nine locations in Cambodia between 2004 and 2013. Our analyses showed that the vast majority of infections are polyclonal (92%) and that P. vivax displays high genetic diversity in Cambodia without apparent geographic stratification. Interestingly, our analyses also revealed that the proportion of monoclonal infections significantly increased between 2004 and 2013, possibly suggesting that malaria control strategies in Cambodia may be successfully affecting the parasite population.

CONCLUSIONS/SIGNIFICANCE

Our findings demonstrate that this high-throughput genotyping assay is efficient in characterizing P. vivax diversity and can provide valuable insights to assess the efficacy of malaria elimination programs or to monitor the spread of specific parasites.

Microsatellite Genotyping of Plasmodium vivax Isolates from Pregnant Women in Four Malaria Endemic Countries

Plasmodium vivax is the most widely distributed human parasite and the main cause of human malaria outside the African continent. However, the knowledge about the genetic variability of P. vivax is limited when compared to the information available for P. falciparum. We present the results of a study aimed at characterizing the genetic structure of P. vivax populations obtained from pregnant women from different malaria endemic settings. Between June 2008 and October 2011 nearly 2000 pregnant women were recruited during routine antenatal care at each site and followed up until delivery. A capillary blood sample from the study participants was collected for genotyping at different time points. Seven P. vivax microsatellite markers were used for genotypic characterization on a total of 229 P. vivax isolates obtained from Brazil, Colombia, India and Papua New Guinea. In each population, the number of alleles per locus, the expected heterozygosity and the levels of multilocus linkage disequilibrium were assessed. The extent of genetic differentiation among populations was also estimated. Six microsatellite loci on 137 P. falciparum isolates from three countries were screened for comparison. The mean value of expected heterozygosity per country ranged from 0.839 to 0.874 for P. vivax and from 0.578 to 0.758 for P. falciparum. P. vivax populations were more diverse than those of P. falciparum. In some of the studied countries, the diversity of P. vivax population was very high compared to the respective level of endemicity. The level of inter-population differentiation was moderate to high in all P. vivax and P. falciparum populations studied.

Further Evidence of Increasing Diversity of Plasmodium vivax in the Republic of Korea in Recent Years

BACKGROUND

Vivax malaria was successfully eliminated from the Republic of Korea (ROK) in the late 1970s but re-emerged in 1993. Two decades later as the ROK enters the final stages of malaria elimination, dedicated surveillance of the local P. vivax population is critical. We apply a population genetic approach to gauge P. vivax transmission dynamics in the ROK between 2010 and 2012.

METHODOLOGY/PRINCIPAL FINDINGS

P. vivax positive blood samples from 98 autochthonous cases were collected from patients attending health centers in the ROK in 2010 (n = 27), 2011 (n = 48) and 2012 (n = 23). Parasite genotyping was undertaken at 9 tandem repeat markers. Although not reaching significance, a trend of increasing population diversity was observed from 2010 (HE = 0.50 ± 0.11) to 2011 (HE = 0.56 ± 0.08) and 2012 (HE = 0.60 ± 0.06). Conversely, linkage disequilibrium declined during the same period: IAS = 0.15 in 2010 (P = 0.010), 0.09 in 2011 (P = 0.010) and 0.05 in 2012 (P = 0.010). In combination with data from other ROK studies undertaken between 1994 and 2007, our results are consistent with increasing parasite divergence since re-emergence. Polyclonal infections were rare (3% infections) suggesting that local out-crossing alone was unlikely to explain the increased divergence. Cases introduced from an external reservoir may therefore have contributed to the increased diversity. Aside from one isolate, all infections carried a short MS20 allele (142 or 149 bp), not observed in other studies in tropical endemic countries despite high diversity, inferring that these regions are unlikely reservoirs.

CONCLUSIONS

Whilst a number of factors may explain the observed population genetic trends, the available evidence suggests that an external geographic reservoir with moderate diversity sustains the majority of P. vivax infection in the ROK, with important implications for malaria elimination.

Population Genetics of Plasmodium vivax in Four Rural Communities in Central Vietnam

Background

The burden of malaria in Vietnam has drastically reduced, prompting the National Malaria Control Program to officially engage in elimination efforts. Plasmodium vivax is becoming increasingly prevalent, remaining a major problem in the country's central and southern provinces. A better understanding of P. vivax genetic diversity and structure of local parasite populations will provide baseline data for the evaluation and improvement of current efforts for control and elimination. The aim of this study was to examine the population genetics and structure of P. vivax isolates from four communities in Tra Leng commune, Nam Tra My district in Quang Nam, Central Vietnam.

Methodology/Principal Findings

P. vivax mono infections collected from 234 individuals between April 2009 and December 2010 were successfully analyzed using a panel of 14 microsatellite markers. Isolates displayed moderate genetic diversity (He = 0.68), with no significant differences between study communities. Polyclonal infections were frequent (71.4%) with a mean multiplicity of infection of 1.91 isolates/person. Low but significant genetic differentiation (F_ST value from -0.05 to 0.18) was observed between the community across the river and the other communities. Strong linkage disequilibrium (IAS = 0.113, p < 0.001) was detected across all communities, suggesting gene flow within and among them. Using multiple approaches, 101 haplotypes were grouped into two genetic clusters, while 60.4% of haplotypes were admixed.

Conclusions/Significance

In this area of Central Vietnam, where malaria transmission has decreased significantly over the past decade, there was moderate genetic diversity and high occurrence of polyclonal infections. Local human populations have frequent social and economic interactions that facilitate gene flow and inbreeding among parasite populations, while decreasing population structure. Findings provide important information on parasites populations circulating in the study area and are relevant to current malaria elimination efforts.

In Vietnam, Plasmodium vivax (P. vivax) is the second most frequent human malaria parasite and a major obstacle to countrywide malaria elimination. Knowing the local parasite structure is useful for elimination efforts. Therefore, we analyzed, with a panel of 14 microsatellite markers, 234 P. vivax mono infections in blood samples collected from 4 communities in central Vietnam. Genetic diversity in the population was moderate; a high occurrence of polyclonal infections and significant linkage disequilibrium were detected, suggesting inbreeding or recombination between highly related haplotypes. In addition, both genetic differentiation and population structure was low and only detected between communities at each side of the river. Those results suggest gene flow between study communities with the river defining a moderate geographical barrier. Future studies should determine how this genetic variation is maintained in an area of extremely low transmission.

Population Genetics of Plasmodium vivax in the Peruvian Amazon

Background

Characterizing the parasite dynamics and population structure provides useful information to understand the dynamic of transmission and to better target control interventions. Despite considerable efforts for its control, vivax malaria remains a major health problem in Peru. In this study, we have explored the population genetics of Plasmodium vivax isolates from Iquitos, the main city in the Peruvian Amazon, and 25 neighbouring peri-urban as well as rural villages along the Iquitos-Nauta Road.

Methodology/ Results

From April to December 2008, 292 P. vivax isolates were collected and successfully genotyped using 14 neutral microsatellites. Analysis of the molecular data revealed a similar proportion of monoclonal and polyclonal infections in urban areas, while in rural areas monoclonal infections were predominant (p = 0.002). Multiplicity of infection was higher in urban (MOI = 1.5–2) compared to rural areas (MOI = 1) (p = 0.003). The level of genetic diversity was similar in all areas (He = 0.66–0.76, p = 0.32) though genetic differentiation between areas was substantial (PHI_PT = 0.17, p<0.0001). Principal coordinate analysis showed a marked differentiation between parasites from urban and rural areas. Linkage disequilibrium was detected in all the areas (IAs = 0.08–0.49, for all p<0.0001). Gene flow among the areas was stablished through Bayesian analysis of migration models. Recent bottleneck events were detected in 4 areas and a recent parasite expansion in one of the isolated areas. In total, 87 unique haplotypes grouped in 2 or 3 genetic clusters described a sub-structured parasite population.

Conclusion/Significance

Our study shows a sub-structured parasite population with clonal propagation, with most of its components recently affected by bottleneck events. Iquitos city is the main source of parasite spreading for all the peripheral study areas. The routes of transmission and gene flow and the reduction of the parasite population described are important from the public health perspective as well for the formulation of future control policies.

We present the population genetics of malaria vivax parasites in a large area of the Peruvian Amazon. Our results showed that the parasite population had a predominant clonal propagation, reproducing themselves with identically or closely related parasites; therefore, the same genetic characteristics are maintained in the offspring. The clonal propagation may favour the higher levels of genetic differentiation among the parasites from isolated areas compared to areas where human migration is common. The patterns of gene flow have been established, finding Iquitos city as a reservoir of parasite genetic variability. Moreover, a recent reduction of the parasite population was observed in areas where recent control activities were performed. This research provides a picture of the nature and dynamics of the parasite population which have a significant impact in the malaria epidemiology; therefore, this knowledge is crucial for the development of efficient control policies.

Targeting vivax malaria in the Asia Pacific: The Asia Pacific Malaria Elimination Network Vivax Working Group

The Asia Pacific Malaria Elimination Network (APMEN) is a collaboration of 18 country partners committed to eliminating malaria from within their borders. Over the past 5 years, APMEN has helped to build the knowledge, tools and in-country technical expertise required to attain this goal. At its inaugural meeting in Brisbane in 2009, Plasmodium vivax infections were identified across the region as a common threat to this ambitious programme; the APMEN Vivax Working Group was established to tackle specifically this issue. The Working Group developed a four-stage strategy to identify knowledge gaps, build regional consensus on shared priorities, generate evidence and change practice to optimize malaria elimination activities. This case study describes the issues faced and the solutions found in developing this robust strategic partnership between national programmes and research partners within the Working Group. The success of the approach adopted by the group may facilitate similar applications in other regions seeking to deploy evidence-based policy and practice.

Whole Genome Sequencing of Field Isolates Reveals Extensive Genetic Diversity in Plasmodium vivax from Colombia

Plasmodium vivax is the most prevalent malarial species in South America and exerts a substantial burden on the populations it affects. The control and eventual elimination of P. vivax are global health priorities. Genomic research contributes to this objective by improving our understanding of the biology of P. vivax and through the development of new genetic markers that can be used to monitor efforts to reduce malaria transmission. Here we analyze whole-genome data from eight field samples from a region in Cordóba, Colombia where malaria is endemic. We find considerable genetic diversity within this population, a result that contrasts with earlier studies suggesting that P. vivax had limited diversity in the Americas. We also identify a selective sweep around a substitution known to confer resistance to sulphadoxine-pyrimethamine (SP). This is the first observation of a selective sweep for SP resistance in this species. These results indicate that P. vivax has been exposed to SP pressure even when the drug is not in use as a first line treatment for patients afflicted by this parasite. We identify multiple non-synonymous substitutions in three other genes known to be involved with drug resistance in Plasmodium species. Finally, we found extensive microsatellite polymorphisms. Using this information we developed 18 polymorphic and easy to score microsatellite loci that can be used in epidemiological investigations in South America.

Variation in Complexity of Infection and Transmission Stability between Neighbouring Populations of Plasmodium vivax in Southern Ethiopia

Background

P. vivax is an important public health burden in Ethiopia, accounting for almost half of all malaria cases. Owing to heterogeneous transmission across the country, a stronger evidence base on local transmission dynamics is needed to optimise allocation of resources and improve malaria interventions.

Methodology and Principal Findings

In a pilot evaluation of local level P. vivax molecular surveillance in southern Ethiopia, the diversity and population structure of isolates collected between May and November 2013 were investigated. Blood samples were collected from microscopy positive P. vivax patients recruited to clinical and cross-sectional surveys from four sites: Arbaminch, Halaba, Badawacho and Hawassa. Parasite genotyping was undertaken at nine tandem repeat markers. Eight loci were successfully genotyped in 197 samples (between 36 and 59 per site). Heterogeneity was observed in parasite diversity and structure amongst the sites. Badawacho displayed evidence of unstable transmission, with clusters of identical clonal infections. Linkage disequilibrium in Badawacho was higher (I_AS = 0.32, P = 0.010) than in the other populations (I_AS range = 0.01–0.02) and declined markedly after adjusting for identical infections (I_AS = 0.06, P = 0.010). Other than Badawacho (H_E = 0.70), population diversity was equivalently high across the sites (H_E = 0.83). Polyclonal infections were more frequent in Hawassa (67%) than the other populations (range: 8–44%). Despite the variable diversity, differentiation between the sites was low (F_ST range: 5 x 10⁻³–0.03).

Conclusions

Marked variation in parasite population structure likely reflects differing local transmission dynamics. Parasite genotyping in these heterogeneous settings has potential to provide important complementary information with which to optimise malaria control interventions.

Multiple-clone infections of Plasmodium vivax: definition of a panel of markers for molecular epidemiology

Background

Plasmodium vivax infections commonly contain multiple genetically distinct parasite clones. The detection of multiple-clone infections depends on several factors, such as the accuracy of the genotyping method, and the type and number of the molecular markers analysed. Characterizing the multiplicity of infection has broad implications that range from population genetic studies of the parasite to malaria treatment and control. This study compared and evaluated the efficiency of neutral and non-neutral markers that are widely used in studies of molecular epidemiology to detect the multiplicity of P. vivax infection.

Methods

The performance of six markers was evaluated using 11 mixtures of DNA with well-defined proportions of two different parasite genotypes for each marker. These mixtures were generated by mixing cloned PCR products or patient-derived genomic DNA. In addition, 51 samples of natural infections from the Brazil were genotyped for all markers. The PCR-capillary electrophoresis-based method was used to permit direct comparisons among the markers. The criteria for differentiating minor peaks from artifacts were also evaluated.

Results

The analysis of DNA mixtures showed that the tandem repeat MN21 and the polymorphic blocks 2 (msp1B2) and 10 (msp1B10) of merozoite surface protein-1 allowed for the estimation of the expected ratio of both alleles in the majority of preparations. Nevertheless, msp1B2 was not able to detect the majority of multiple-clone infections in field samples; it identified only 6 % of these infections. The merozoite surface protein-3 alpha and microsatellites (PvMS6 and PvMS7) did not accurately estimate the relative clonal proportions in artificial mixtures, but the microsatellites performed well in detecting natural multiple-clone infections. Notably, the use of a less stringent criterion to score rare alleles significantly increased the sensitivity of the detection of multi-clonal infections.

Conclusions

Depending on the type of marker used, a considerable amplification bias was observed, which may have serious implications for the characterization of the complexity of a P. vivax infection. Based on the performance of markers in artificial mixtures of DNA and natural infections, a minimum panel of four genetic markers (PvMS6, PvMS7, MN21, and msp1B10) was defined, and these markers are highly informative regarding the genetic variability of P. vivax populations.

Electronic supplementary material

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

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

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

Plasmodium vivax Diversity and Population Structure across Four Continents

Plasmodium vivax is the geographically most widespread human malaria parasite. To analyze patterns of microsatellite diversity and population structure across countries of different transmission intensity, genotyping data from 11 microsatellite markers was either generated or compiled from 841 isolates from four continents collected in 1999–2008. Diversity was highest in South-East Asia (mean allelic richness 10.0–12.8), intermediate in the South Pacific (8.1–9.9) Madagascar and Sudan (7.9–8.4), and lowest in South America and Central Asia (5.5–7.2). A reduced panel of only 3 markers was sufficient to identify approx. 90% of all haplotypes in South Pacific, African and SE-Asian populations, but only 60–80% in Latin American populations, suggesting that typing of 2–6 markers, depending on the level of endemicity, is sufficient for epidemiological studies. Clustering analysis showed distinct clusters in Peru and Brazil, but little sub-structuring was observed within Africa, SE-Asia or the South Pacific. Isolates from Uzbekistan were exceptional, as a near-clonal parasite population was observed that was clearly separated from all other populations (F_ST>0.2). Outside Central Asia F_ST values were highest (0.11–0.16) between South American and all other populations, and lowest (0.04–0.07) between populations from South-East Asia and the South Pacific. These comparisons between P. vivax populations from four continents indicated that not only transmission intensity, but also geographical isolation affect diversity and population structure. However, the high effective population size results in slow changes of these parameters. This persistency must be taken into account when assessing the impact of control programs on the genetic structure of parasite populations.

Plasmodium vivax is the predominant malaria parasite in Latin America, Asia and the South Pacific. Different factors are expected to shape diversity and population structure across continents, e.g. transmission intensity which is much lower in South America as compared to Southeast-Asia and the South Pacific, or geographical isolation of P. vivax populations in the South Pacific. We have compiled data from 841 isolates from South and Central America, Africa, Central Asia, Southeast-Asia and the South Pacific typed with a panel of 11 microsatellite markers. Diversity was highest in Southeast-Asia, where transmission is intermediate-high and migration of infected hosts is high, and lowest in South America and Central Asia where malaria transmission is low and focal. Reducing the panel of microsatellites showed that 2–6 markers are sufficient for genotyping for most drug trials and epidemiological studies, as these markers can identify >90% of all haplotypes. Parasites clustered according to continental origin, with high population differentiation between South American and Central Asian populations and the other populations, and lowest differences between Southeast-Asia and the South Pacific. Current attempts to reduce malaria transmission might change this pattern, but only after transmission is reduced for an extended period of time.

Contrasting Transmission Dynamics of Co-endemic Plasmodium vivax and P. falciparum: Implications for Malaria Control and Elimination

BACKGROUND

Outside of Africa, P. falciparum and P. vivax usually coexist. In such co-endemic regions, successful malaria control programs have a greater impact on reducing falciparum malaria, resulting in P. vivax becoming the predominant species of infection. Adding to the challenges of elimination, the dormant liver stage complicates efforts to monitor the impact of ongoing interventions against P. vivax. We investigated molecular approaches to inform the respective transmission dynamics of P. falciparum and P. vivax and how these could help to prioritize public health interventions.

METHODOLOGY/PRINCIPAL FINDINGS

Genotype data generated at 8 and 9 microsatellite loci were analysed in 168 P. falciparum and 166 P. vivax isolates, respectively, from four co-endemic sites in Indonesia (Bangka, Kalimantan, Sumba and West Timor). Measures of diversity, linkage disequilibrium (LD) and population structure were used to gauge the transmission dynamics of each species in each setting. Marked differences were observed in the diversity and population structure of P. vivax versus P. falciparum. In Bangka, Kalimantan and Timor, P. falciparum diversity was low, and LD patterns were consistent with unstable, epidemic transmission, amenable to targeted intervention. In contrast, P. vivax diversity was higher and transmission appeared more stable. Population differentiation was lower in P. vivax versus P. falciparum, suggesting that the hypnozoite reservoir might play an important role in sustaining local transmission and facilitating the spread of P. vivax infections in different endemic settings. P. vivax polyclonality varied with local endemicity, demonstrating potential utility in informing on transmission intensity in this species.

CONCLUSIONS/SIGNIFICANCE

Molecular approaches can provide important information on malaria transmission that is not readily available from traditional epidemiological measures. Elucidation of the transmission dynamics circulating in a given setting will have a major role in prioritising malaria control strategies, particularly against the relatively neglected non-falciparum species.

Uncovering the transmission dynamics of Plasmodium vivax using population genetics

Population genetic analysis of malaria parasites has the power to reveal key insights into malaria epidemiology and transmission dynamics with the potential to deliver tools to support control and elimination efforts. Analyses of parasite genetic diversity have suggested that Plasmodium vivax populations are more genetically diverse and less structured than those of Plasmodium falciparum indicating that P. vivax may be a more ancient parasite of humans and/or less susceptible to population bottlenecks, as well as more efficient at disseminating its genes. These population genetic insights into P. vivax transmission dynamics provide an explanation for its relative resilience to control efforts. Here, we describe current knowledge on P. vivax population genetic structure, its relevance to understanding transmission patterns and relapse and how this information can inform malaria control and elimination programmes.

Plasmodium vivax populations are more genetically diverse and less structured than sympatric Plasmodium falciparum populations

Introduction

The human malaria parasite, Plasmodium vivax, is proving more difficult to control and eliminate than Plasmodium falciparum in areas of co-transmission. Comparisons of the genetic structure of sympatric parasite populations may provide insight into the mechanisms underlying the resilience of P. vivax and can help guide malaria control programs.

Methodology/Principle findings

P. vivax isolates representing the parasite populations of four areas on the north coast of Papua New Guinea (PNG) were genotyped using microsatellite markers and compared with previously published microsatellite data from sympatric P. falciparum isolates. The genetic diversity of P. vivax (H_e = 0.83–0.85) was higher than that of P. falciparum (H_e = 0.64–0.77) in all four populations. Moderate levels of genetic differentiation were found between P. falciparum populations, even over relatively short distances (less than 50 km), with 21–28% private alleles and clear geospatial genetic clustering. Conversely, very low population differentiation was found between P. vivax catchments, with less than 5% private alleles and no genetic clustering observed. In addition, the effective population size of P. vivax (30353; 13043–69142) was larger than that of P. falciparum (18871; 8109–42986).

Conclusions/Significance

Despite comparably high prevalence, P. vivax had higher diversity and a panmictic population structure compared to sympatric P. falciparum populations, which were fragmented into subpopulations. The results suggest that in comparison to P. falciparum, P. vivax has had a long-term large effective population size, consistent with more intense and stable transmission, and limited impact of past control and elimination efforts. This underlines suggestions that more intensive and sustained interventions will be needed to control and eventually eliminate P. vivax. This research clearly demonstrates how population genetic analyses can reveal deeper insight into transmission patterns than traditional surveillance methods.

The neglected human malaria parasite Plasmodium vivax is responsible for a large proportion of the global malaria burden. Efforts to control malaria have revealed that P. vivax is more resilient than the other major human malaria parasite, Plasmodium falciparum. This study utilised population genetics to compare patterns of P. vivax and P. falciparum transmission in Papua New Guinea, a region where infection rates of the two species are similar. The results demonstrated that P. vivax populations are more genetically diverse than those of P. falciparum suggestive of a parasite population that is more resilient to environmental challenges, undergoing higher levels of interbreeding locally and between distant parasite populations. Unique characteristics of P. vivax such as relapse, which allows different strains from past infections to produce subsequent infections, may provide more opportunities for the exchange and dissemination of genetic material. The contrasting patterns observed for the two species may be the result of a differential impact of past elimination attempts and indicate that more rigorous interventions will be needed in efforts to control and eventually eliminate P. vivax.

Plasmodium vivax and Plasmodium falciparum at the crossroads of exchange among islands in Vanuatu: implications for malaria elimination strategies

Understanding the transmission and movement of Plasmodium parasites is crucial for malaria elimination and prevention of resurgence. Located at the limit of malaria transmission in the Pacific, Vanuatu is an ideal candidate for elimination programs due to low endemicity and the isolated nature of its island setting. We analyzed the variation in the merozoite surface protein 1 (msp1) and the circumsporozoite protein (csp) of P. falciparum and P. vivax populations to examine the patterns of gene flow and population structures among seven sites on five islands in Vanuatu. Genetic diversity was in general higher in P. vivax than P. falciparum from the same site. In P. vivax, high genetic diversity was likely maintained by greater extent of gene flow among sites and among islands. Consistent with the different patterns of gene flow, the proportion of genetic variance found among islands was substantially higher in P. falciparum (28.81–31.23%) than in P. vivax (-0.53–3.99%). Our data suggest that the current island-by-island malaria elimination strategy in Vanuatu, while adequate for P. falciparum elimination, might need to be complemented with more centrally integrated measures to control P. vivax movement across islands.

Using Amplicon Deep Sequencing to Detect Genetic Signatures of Plasmodium vivax Relapse

Plasmodium vivax infections often recur due to relapse of hypnozoites from the liver. In malaria-endemic areas, tools to distinguish relapse from reinfection are needed. We applied amplicon deep sequencing to P. vivax isolates from 78 Cambodian volunteers, nearly one-third of whom suffered recurrence at a median of 68 days. Deep sequencing at a highly variable region of the P. vivax merozoite surface protein 1 gene revealed impressive diversity-generating 67 unique haplotypes and detecting on average 3.6 cocirculating parasite clones within individuals, compared to 2.1 clones detected by a combination of 3 microsatellite markers. This diversity enabled a scheme to classify over half of recurrences as probable relapses based on the low probability of reinfection by multiple recurring variants. In areas of high P. vivax diversity, targeted deep sequencing can help detect genetic signatures of relapse, key to evaluating antivivax interventions and achieving a better understanding of relapse-reinfection epidemiology.

Genetic diversity of Plasmodium vivax over time and space: a community-based study in rural Amazonia

To examine how community-level genetic diversity of the malaria parasite Plasmodium vivax varies across time and space, we investigated the dynamics of parasite polymorphisms during the early phases of occupation of a frontier settlement in the Amazon Basin of Brazil. Microsatellite characterization of 84 isolates of P. vivax sampled over 3 years revealed a moderate-to-high genetic diversity (mean expected heterozygosity, 0.699), with a large proportion (78.5%) of multiple-clone infections (MCI), but also a strong multilocus linkage disequilibrium (LD) consistent with rare outcrossing. Little temporal and no spatial clustering was observed in the distribution of parasite haplotypes. A single microsatellite haplotype was shared by 3 parasites collected during an outbreak; all other 81 haplotypes were recovered only once. The lowest parasite diversity, with the smallest proportion of MCI and the strongest LD, was observed at the time of the outbreak, providing a clear example of epidemic population structure in a human pathogen. Population genetic parameters returned to pre-outbreak values during last 2 years of study, despite the concomitant decline in malaria incidence. We suggest that parasite genotyping can be useful for tracking the spread of new parasite strains associated with outbreaks in areas approaching malaria elimination.

Distinct patterns of diversity, population structure and evolution in the AMA1 genes of sympatric Plasmodium falciparum and Plasmodium vivax populations of Papua New Guinea from an area of similarly high transmission

Background

As Plasmodium falciparum and Plasmodium vivax co-exist in most malaria-endemic regions outside sub-Saharan Africa, malaria control strategies in these areas must target both species in order to succeed. Population genetic analyses can predict the effectiveness of interventions including vaccines, by providing insight into patterns of diversity and evolution. The aim of this study was to investigate the population genetics of leading malaria vaccine candidate AMA1 in sympatric P. falciparum and P. vivax populations of Papua New Guinea (PNG), an area of similarly high prevalence (Pf = 22.3 to 38.8%, Pv = 15.3 to 31.8%).

Methods

A total of 72 Pfama1 and 102 Pvama1 sequences were collected from two distinct areas, Madang and Wosera, on the highly endemic PNG north coast.

Results

Despite a greater number of polymorphic sites in the AMA1 genes of P. falciparum (Madang = 52; Wosera = 56) compared to P. vivax (Madang = 36, Wosera = 34), the number of AMA1 haplotypes, haplotype diversity (Hd) and recombination (R) was far lower for P. falciparum (Madang = 12, Wosera = 20; Hd ≤0.92, R ≤45.8) than for P. vivax (Madang = 50, Wosera = 38; Hd = 0.99, R = ≤70.9). Balancing selection was detected only within domain I of AMA1 for P. vivax, and in both domains I and III for P. falciparum.

Conclusions

Higher diversity in the genes encoding P. vivax AMA1 than in P. falciparum AMA1 in this highly endemic area has important implications for development of AMA1-based vaccines in PNG and beyond. These results also suggest a smaller effective population size of P. falciparum compared to P. vivax, a finding that warrants further investigation. Differing patterns of selection on the AMA1 genes indicate that critical antigenic sites may differ between the species, highlighting the need for independent investigations of these two leading vaccine candidates.