Use cases for genetic epidemiology in malaria elimination

Comparing newly developed SNP barcode panels with microsatellites to explore population genetics of malaria parasites in the Peruvian Amazon

Introduction

Malaria molecular surveillance (MMS) can provide insights into transmission dynamics, guiding national control programs. We previously designed AmpliSeq assays for MMS, which include different traits of interest (resistance markers and pfhrp2/3 deletions), and SNP barcodes to provide population genetics estimates of Plasmodium vivax and Plasmodium falciparum parasites in the Peruvian Amazon. The present study compares the genetic resolution of the barcodes in the AmpliSeq assays with widely used microsatellite (MS) panels to investigate population genetics of Amazonian malaria parasites.

Methods

We analyzed 51 P. vivax and 80 P. falciparum samples from three distinct areas in the Loreto region of the Peruvian Amazon: Nueva Jerusalén (NJ), Mazan (MZ), and Santa Emilia (SE). Population genetics estimates and costs were compared using the SNP barcodes (P. vivax: 40 SNPs and P. falciparum: 28 SNPs) and MS panels (P. vivax: 16 MS and P. falciparum: 7 MS).

Results

The P. vivax genetic diversity (expected heterozygosity, He) trends were similar for both markers: He MS = 0.68-0.78 (p > 0.05) and He SNP = 0.36-0.38 (p > 0.05). P. vivax pairwise genetic differentiation (fixation index, FST) was also comparable: FST-MS = 0.04-0.14 and FST-SNP = 0.03-0.12 (pairwise p > 0.05). In addition, P. falciparum genetic diversity trends (He MS = 0-0.48, p < 0.05; He SNP = 0-0.09, p < 0.05) and pairwise FST comparisons (FST-MS = 0.14-0.65, FST-SNP = 0.19-0.61, pairwise p > 0.05) were concordant between both panels. For P. vivax, no geographic clustering was observed with any panel, whereas for P. falciparum, similar population structure clustering was observed with both markers, assigning most parasites from NJ to a distinct subpopulation from MZ and SE. We found significant differences in detecting polyclonal infections: for P. vivax, MS identified a higher proportion of polyclonal infections than SNP (69% vs. 33%, p = 3.3 × 10-5), while for P. falciparum, SNP and MS detected similar rates (46% vs. 31%, p = 0.21). The AmpliSeq assay had a higher estimated per-sample cost compared to MS ($183 vs. $27-49).

Discussion

The SNP barcodes in the AmpliSeq assays offered comparable results to MS for investigating population genetics in P. vivax and P. falciparum populations, despite some discrepancies in determining polyclonality. Given both panels have their respective advantages and limitations, the choice between both should be guided by research objectives, costs, and resource availability.

India should invest in the expansion of genomic epidemiology for vector-borne diseases filariasis, malaria and visceral leishmaniasis that are targeted for elimination

Genomic epidemiology (GE) is an integration of genomics and epidemiology. The field has evolved significantly in the past decade, enhancing our understanding of genetic susceptibility, drug resistance, disease transmission patterns, outbreak surveillance, and vaccine development. It employs an arsenal of advanced tools such as whole-genome sequencing and single-nucleotide polymorphisms for analysis, tracing pathogen evolution, mapping genetic variations, and tracking drug resistance. The role of GE in infectious disease management extends beyond outbreak control to routine public health practices, precision medicine, and determining treatment policies. The expansion of GE can significantly bolster global health defenses by effectively enabling the detection and response to emerging health threats. However, challenges such as sampling bias, data quality, integration, standardization of computational pipelines, and need for trained personnel remain. To tackle these challenges, we must invest in building capacity, improving infrastructure, providing training, and fostering collaborations between scientists and public health officials. Concerted efforts must focus on overcoming existing hurdles and promoting seamless integration of basic research into public health frameworks to fully realize its potential. It is timely for India to rapidly expand its base in GE to gain valuable insights into genetic variations and disease susceptibilities. This will provide a fillip towards eliminating the three dominant vector-borne diseases in India: filariasis, malaria, and visceral leishmaniasis.

Comparing newly developed SNP barcode panels with microsatellites to explore population genetics of malaria parasites in the Peruvian Amazon

Malaria molecular surveillance (MMS) can provide insights into transmission dynamics, guiding national control/elimination programs. Considering the genetic differences among parasites from different areas in the Peruvian Amazon, we previously designed SNP barcode panels for Plasmodium vivax (Pv) and P. falciparum (Pf), integrated into AmpliSeq assays, to provide population genetics estimates of malaria parasites. These AmpliSeq assays are ideal for MMS: multiplexing different traits of interest, applicable to many use cases, and high throughput for large numbers of samples. The present study compares the genetic resolution of the SNP barcode panels in the AmpliSeq assays with widely used microsatellite (MS) panels to investigate Amazonian malaria parasites. Malaria samples collected in remote areas of the Peruvian Amazon (51 Pv & 80 Pf samples) were characterized using the Ampliseq assays and MS. Population genetics estimates (complexity of infection, genetic diversity and differentiation, and population structure) were compared using the SNP barcodes (Pv: 40 SNPs & Pf: 28 SNPs) and MS panels (Pv: 16 MS & Pf: 7 MS). The genetic diversity of Pv (expected heterozygosity, He ) was similar across the subpopulations for both makers: He MS = 0.68 - 0.78 (p = 0.23) and He SNP = 0.36 - 0.38 (p = 0.80). Pairwise genetic differentiation (fixation index, F ST ) was also comparable: F ST-MS = 0.04 - 0.14 and F ST-SNP = 0.03 - 0.12 (p = 0.34 - 0.85). No geographic clustering was observed with any panel. In addition, Pf genetic diversity trends ( He MS = 0 - 0.48 p = 0.03 - 1; He SNP = 0 - 0.09, p = 0.03 - 1) and pairwise F ST comparisons (F ST-MS = 0.14 - 0.65, F ST-SNP = 0.19 - 0.61, p = 0.24 - 0.83) were concordant between the panels. Similar population structure clustering was observed with both SNP and MS, highlighting one Pf subpopulation in an indigenous community. The SNP barcodes in the Pv AmpliSeq v2 Peru and Pf AmpliSeq v1 Peru assays offer comparable results to MS panels when investigating population genetics in Pv and Pv populations. Therefore, the AmpliSeq assays can efficiently characterize malaria transmission dynamics and population structure and support malaria elimination efforts in Peru.

Sensitive and modular amplicon sequencing of Plasmodium falciparum diversity and resistance for research and public health

Targeted amplicon sequencing is a powerful and efficient tool to interrogate the P. falciparum genome and generate actionable data from infections to complement traditional malaria epidemiology. For maximum impact, genomic tools should be multi-purpose, robust, sensitive and reproducible. We developed, characterized, and implemented MAD4HatTeR, an amplicon sequencing panel based on Multiplex Amplicons for Drug, Diagnostic, Diversity, and Differentiation Haplotypes using Targeted Resequencing, along with a bioinformatic pipeline for data analysis. MAD4HatTeR targets 165 highly diverse loci, focusing on multiallelic microhaplotypes; key markers for drug and diagnostic resistance, including duplications and deletions; and csp and potential vaccine targets. In addition, it can detect non-falciparum Plasmodium species. We used laboratory control and field sample data to demonstrate the high sensitivity and robustness of the panel. The successful implementation of this method in five laboratories, including three in malaria-endemic African countries, showcases its feasibility in generating reproducible data across laboratories. Finally, we introduce an analytical approach to detect gene duplications and deletions from amplicon sequencing data. MAD4HatTeR is thus a powerful research tool and a robust resource for malaria public health surveillance and control.

Genetic surveillance reveals low, sustained malaria transmission with clonal replacement in Sao Tome and Principe

Despite efforts to eliminate malaria in Sao Tome and Principe (STP), cases have recently increased. Understanding residual transmission structure is crucial for developing effective elimination strategies. This study collected surveillance data and generated amplicon sequencing data from 980 samples between 2010 and 2016 to examine the genetic structure of the parasite population. The mean multiplicity of infection (MOI) was 1.3, with 11% polyclonal infections, indicating low transmission intensity. Temporal trends of these genetic metrics did not align with incidence rates, suggesting that changes in genetic metrics may not straightforwardly reflect changes in transmission intensity, particularly in low transmission settings where genetic drift and importation have a substantial impact. While 88% of samples were genetically linked, continuous turnover in genetic clusters and changes in drug-resistance haplotypes were observed. Principal component analysis revealed some STP samples were genetically similar to those from Central and West Africa, indicating possible importation. These findings highlight the need to prioritize several interventions such as targeted interventions against transmission hotspots, reactive case detection, and strategies to reduce the introduction of new parasites into this island nation as it approaches elimination. This study also serves as a case study for implementing genetic surveillance in a low transmission setting.

Genomic surveillance of malaria parasites in an indigenous community in the Peruvian Amazon

Hard-to-reach communities represent Peru's main challenge for malaria elimination, but information about transmission in these areas is scarce. Here, we assessed Plasmodium vivax (Pv) and P. falciparum (Pf) transmission dynamics, resistance markers, and Pf hrp2/3 deletions in Nueva Jerusalén (NJ), a remote, indigenous community in the Peruvian Amazon with high population mobility. We collected samples from November 2019 to May 2020 by active (ACD) and passive case detection (PCD) in NJ. Parasites were identified with microscopy and PCR. Then, we analyzed a representative set of positive-PCR samples (Pv = 68, Pf = 58) using highly-multiplexed deep sequencing assays (AmpliSeq) and compared NJ parasites with ones from other remote Peruvian areas using population genetics indexes. The ACD intervention did not reduce malaria cases in the short term, and persistent malaria transmission was observed (at least one Pv infection was detected in 96% of the study days). In Nueva Jerusalen, the Pv population had modest genetic diversity (He = 0.27). Pf population had lower diversity (He = 0.08) and presented temporal clustering, one of these clusters linked to an outbreak in February 2020. Moreover, Pv and Pf parasites from NJ exhibited variable levels of differentiation (Pv Fst = 0.07-0.52 and Pf Fst = 0.11-0.58) with parasites from other remote areas. No artemisin resistance mutations but chloroquine (57%) and sulfadoxine-pyrimethamine (35-67%) were detected in NJ's Pf parasites. Moreover, pfhrp2/3 gene deletions were common (32-50% of parasites with one or both genes deleted). The persistent Pv transmission and the detection of a Pf outbreak with parasites genetically distinct from the local ones highlight the need for tailored interventions focusing on mobility patterns and imported infections in remote areas to eliminate malaria in the Peruvian Amazon.

Plasmodium vivax genomic surveillance in the Peruvian Amazon with Pv AmpliSeq assay

BACKGROUND

Plasmodium vivax is the most predominant malaria species in Latin America, constituting 71.5% of malaria cases in 2021. With several countries aiming for malaria elimination, it is crucial to prioritize effectiveness of national control programs by optimizing the utilization of available resources and strategically implementing necessary changes. To support this, there is a need for innovative approaches such as genomic surveillance tools that can investigate changes in transmission intensity, imported cases and sources of reintroduction, and can detect molecular markers associated with drug resistance.

METHODOLOGY/PRINCIPAL FINDINGS

Here, we apply a modified highly-multiplexed deep sequencing assay: Pv AmpliSeq v2 Peru. The tool targets a newly developed 41-SNP Peru barcode for parasite population analysis within Peru, the 33-SNP vivaxGEN-geo panel for country-level classification, and 11 putative drug resistance genes. It was applied to 230 samples from the Peruvian Amazon (2007-2020), generating baseline surveillance data. We observed a heterogenous P. vivax population with high diversity and gene flow in peri-urban areas of Maynas province (Loreto region) with a temporal drift using all SNPs detected by the assay (nSNP = 2909). In comparison, in an indigenous isolated area, the parasite population was genetically differentiated (FST = 0.07-0.09) with moderate diversity and high relatedness between isolates in the community. In a remote border community, a clonal P. vivax cluster was identified, with distinct haplotypes in drug resistant genes and ama1, more similar to Brazilian isolates, likely representing an introduction of P. vivax from Brazil at that time. To test its applicability for Latin America, we evaluated the SNP Peru barcode in P. vivax genomes from the region and demonstrated the capacity to capture local population clustering at within-country level.

CONCLUSIONS/SIGNIFICANCE

Together this data shows that P. vivax transmission is heterogeneous in different settings within the Peruvian Amazon. Genetic analysis is a key component for regional malaria control, offering valuable insights that should be incorporated into routine surveillance.

Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana

Plasmodium parasites, the causal agents of malaria, are eukaryotic organisms that obligately undergo sexual recombination within mosquitoes. In low transmission settings, parasites recombine with themselves, and the clonal lineage is propagated rather than broken up by outcrossing. We investigated whether stochastic/neutral factors drive the persistence and abundance of Plasmodium falciparum clonal lineages in Guyana, a country with relatively low malaria transmission, but the only setting in the Americas in which an important artemisinin resistance mutation (pfk13 C580Y) has been observed. We performed whole genome sequencing on 1,727 Plasmodium falciparum samples collected from infected patients across a five-year period (2016-2021). We characterized the relatedness between each pair of monoclonal infections (n = 1,409) through estimation of identity-by-descent (IBD) and also typed each sample for known or candidate drug resistance mutations. A total of 160 multi-isolate clones (mean IBD ≥ 0.90) were circulating in Guyana during the study period, comprising 13 highly related clusters (mean IBD ≥ 0.40). In the five-year study period, we observed a decrease in frequency of a mutation associated with artemisinin partner drug (piperaquine) resistance (pfcrt C350R) and limited co-occurence of pfcrt C350R with duplications of plasmepsin 2/3, an epistatic interaction associated with piperaquine resistance. We additionally observed 61 nonsynonymous substitutions that increased markedly in frequency over the study period as well as a novel pfk13 mutation (G718S). However, P. falciparum clonal dynamics in Guyana appear to be largely driven by stochastic factors, in contrast to other geographic regions, given that clones carrying drug resistance polymorphisms do not demonstrate enhanced persistence or higher abundance than clones carrying polymorphisms of comparable frequency that are unrelated to resistance. The use of multiple artemisinin combination therapies in Guyana may have contributed to the disappearance of the pfk13 C580Y mutation.

Meta-analysis on Plasmodium falciparum sulfadoxine-pyrimethamine resistance-conferring mutations in India identifies hot spots for genetic surveillance

BACKGROUND

India is on track to eliminate malaria by 2030 but emerging resistance to first-line antimalarials is a recognised threat. Two instances of rapid development, spread, and natural selection of drug-resistant mutant parasites in India (chloroquine across the country and artesunate + sulfadoxine-pyrimethamine [AS+SP] in the northeastern states) translated into drug policy changes for Plasmodium falciparum malaria in 2010 and 2013, respectively. Considering these rapid changes in the SP drug resistance-conferring mutation profile of P. falciparum, there is a need to systematically monitor the validated mutations in Pfdhfr and Pfdhps genes across India alongside AS+SP therapeutic efficacy studies. There has been no robust, systematic countrywide surveillance reported for these parameters in India, hence the current study was undertaken.

METHODS

Studies that reported data on WHO-validated SP resistance markers in P. falciparum across India from 2008 to January 2023 were included. Five major databases, PubMedⓇ, Web of ScienceTM, ScopusⓇ, EmbaseⓇ, and Google Scholar, were exhaustively searched. Individual and pooled prevalence estimates of mutations were obtained through random- and fixed-effect models. Data were depicted using forest plots created with a 95% confidence interval. The study is registered with PROSPERO (CRD42021236012).

RESULTS

A total of 37 publications, and 533 Pfdhfr and 134 Pfdhps National Centre of Biotechnology Information (NCBI) DNA sequences were included from >4000 samples. The study included information from 80 districts, 21 states and 3 union territories (UTs) from India. The two PfDHFR mutations, C59R (62%) and S108N (74%), were the most prevalent mutations (pooled estimates 61% and 71%, respectively) and appeared to be stabilised/fixed. Although rarest overall, the prevalence of I164L was observed to be as high as 32%. The PfDHFR double mutants were the most prevalent overall (51%; pooled 42%). The prevalence of triple and quadruple mutations was 6% and 5%, respectively, and is an immediate concern for some states. The most prevalent PfDHPS mutation was A437G (39%), followed by K540E (25%) and A581G (12%). There was a low overall prevalence of PfDHFR/PfDHPS quintuple and sextuple mutations but surveillance for these mutations is critical for some areas.

CONCLUSION

The analyses span the two critical policy changes, highlight the areas of concern, and guide policymakers in strategising and refining the anti-malaria drug policy for malaria elimination. The results of the analyses also highlight the SP-resistance hot spots, critical gaps and challenges, and indicate that focal and local malaria genetic surveillance (including drug-resistance markers) is needed until malaria is successfully eliminated.

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

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

Genomic surveillance of malaria parasites in an indigenous community in the Peruvian Amazon

Hard-to-reach communities represent Peru's main challenge for malaria elimination, but information about transmission in these areas is scarce. Here, we assessed Plasmodium vivax (Pv) and P. falciparum (Pf) transmission dynamics, resistance markers, and Pf hrp2/3 deletions in Nueva Jerusalén (NJ), a remote, indigenous community in the Peruvian Amazon with high population mobility. We collected samples from November 2019 to May 2020 by active (ACD) and passive case detection (PCD) in NJ. Parasites were identified with microscopy and PCR. Then, we analyzed a representative set of positive-PCR samples (Pv = 68, Pf = 58) using highly-multiplexed deep sequencing assays (AmpliSeq) and compared NJ parasites with ones from other remote Peruvian areas using population genetics indexes. The ACD intervention did not reduce malaria cases in the short term, and persistent malaria transmission was observed (at least one Pv infection was detected in 96% of the study days). In Nueva Jerusalen, the Pv population had modest genetic diversity (He = 0.27). Pf population had lower diversity (He = 0.08) and presented temporal clustering, one of these clusters linked to an outbreak in February 2020. Moreover, Pv and Pf parasites from NJ exhibited variable levels of differentiation (Pv Fst = -0.52 & Pf Fst = 0.11-0.58) with parasites from other remote areas. No artemisin resistance mutations but chloroquine (57%) and sulfadoxine-pyrimethamine (35-67%) were detected in NJ's Pf parasites. Moreover, pfhrp2/3 gene deletions were common (32-50% of parasites with one or both genes deleted). The persistent Pv transmission and the detection of a Pf outbreak with parasites genetically distinct from the local ones highlight the need for tailored interventions focusing on mobility patterns and imported infections in remote areas to eliminate malaria in the Peruvian Amazon.

Using a mobile nanopore sequencing lab for end-to-end genomic surveillance of Plasmodium falciparum: A feasibility study

Genomic epidemiology holds promise for malaria control and elimination efforts, for example by informing on Plasmodium falciparum genetic diversity and prevalence of mutations conferring anti-malarial drug resistance. Limited sequencing infrastructure in many malaria-endemic areas prevents the rapid generation of genomic data. To address these issues, we developed and validated assays for P. falciparum nanopore sequencing in endemic sites using a mobile laboratory, targeting key antimalarial drug resistance markers and microhaplotypes. Using two multiplexed PCR reactions, we amplified six highly polymorphic microhaplotypes and ten drug resistance markers. We developed a bioinformatics workflow that allows genotyping of polyclonal malaria infections, including minority clones. We validated the panels on mock dried blood spot (DBS) and rapid diagnostic test (RDT) samples and archived DBS, demonstrating even, high read coverage across amplicons (range: 580x to 3,212x median coverage), high haplotype calling accuracy, and the ability to explore within-sample diversity of polyclonal infections. We field-tested the feasibility of rapid genotyping in Zanzibar in close collaboration with the local malaria elimination program using DBS and routinely collected RDTs as sample inputs. Our assay identified haplotypes known to confer resistance to known antimalarials in the dhfr, dhps and mdr1 genes, but no evidence of artemisinin partial resistance. Most infections (60%) were polyclonal, with high microhaplotype diversity (median HE = 0.94). In conclusion, our assays generated actionable data within a few days, and we identified current challenges for implementing nanopore sequencing in endemic countries to accelerate malaria control and elimination.

Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana

Plasmodium parasites, the causal agents of malaria, are eukaryotic organisms that obligately undergo sexual recombination within mosquitoes. However, in low transmission settings where most mosquitoes become infected with only a single parasite clone, parasites recombine with themselves, and the clonal lineage is propagated rather than broken up by outcrossing. We investigated whether stochastic/neutral factors drive the persistence and abundance of Plasmodium falciparum clonal lineages in Guyana, a country with relatively low malaria transmission, but the only setting in the Americas in which an important artemisinin resistance mutation (pfk13 C580Y) has been observed. To investigate whether this clonality was potentially associated with the persistence and spatial spread of the mutation, we performed whole genome sequencing on 1,727 Plasmodium falciparum samples collected from infected patients across a five-year period (2016-2021). We characterized the relatedness between each pair of monoclonal infections (n=1,409) through estimation of identity by descent (IBD) and also typed each sample for known or candidate drug resistance mutations. A total of 160 clones (mean IBD ≥ 0.90) were circulating in Guyana during the study period, comprising 13 highly related clusters (mean IBD ≥ 0.40). In the five-year study period, we observed a decrease in frequency of a mutation associated with artemisinin partner drug (piperaquine) resistance (pfcrt C350R) and limited co-occurence of pfcrt C350R with duplications of plasmepsin 2/3, an epistatic interaction associated with piperaquine resistance. We additionally report polymorphisms exhibiting evidence of selection for drug resistance or other phenotypes and reported a novel pfk13 mutation (G718S) as well as 61 nonsynonymous substitutions that increased markedly in frequency. However, P. falciparum clonal dynamics in Guyana appear to be largely driven by stochastic factors, in contrast to other geographic regions. The use of multiple artemisinin combination therapies in Guyana may have contributed to the disappearance of the pfk13 C580Y mutation.

Molecular markers for malaria genetic epidemiology: progress and pitfalls

Over recent years, progress in molecular markers for genotyping malaria parasites has enabled informative studies of epidemiology and transmission dynamics. Results have highlighted the value of these tools for surveillance to support malaria control and elimination strategies. There are many different types and panels of markers available for malaria parasite genotyping, and for end users, the nuances of these markers with respect to 'use case', resolution, and accuracy, are not well defined. This review clarifies issues surrounding different molecular markers and their application to malaria control and elimination. We describe available marker panels, use cases, implications for different transmission settings, limitations, access, cost, and data accuracy. The information provided can be used as a guide for molecular epidemiology and surveillance of malaria.

Malaria surveillance reveals parasite relatedness, signatures of selection, and correlates of transmission across Senegal

We here analyze data from the first year of an ongoing nationwide program of genetic surveillance of Plasmodium falciparum parasites in Senegal. The analysis is based on 1097 samples collected at health facilities during passive malaria case detection in 2019; it provides a baseline for analyzing parasite genetic metrics as they vary over time and geographic space. The study's goal was to identify genetic metrics that were informative about transmission intensity and other aspects of transmission dynamics, focusing on measures of genetic relatedness between parasites. We found the best genetic proxy for local malaria incidence to be the proportion of polygenomic infections (those with multiple genetically distinct parasites), although this relationship broke down at low incidence. The proportion of related parasites was less correlated with incidence while local genetic diversity was uninformative. The type of relatedness could discriminate local transmission patterns: two nearby areas had similarly high fractions of relatives, but one was dominated by clones and the other by outcrossed relatives. Throughout Senegal, 58% of related parasites belonged to a single network of relatives, within which parasites were enriched for shared haplotypes at known and suspected drug resistance loci and at one novel locus, reflective of ongoing selection pressure.

Sampling for malaria molecular surveillance

Strategic use of Plasmodium falciparum genetic variation has great potential to inform public health actions for malaria control and elimination. Malaria molecular surveillance (MMS) begins with a strategy to identify and collect parasite samples, guided by public-health priorities. In this review we discuss sampling design practices for MMS and point out epidemiological, biological, and statistical factors that need to be considered. We present examples for different use cases, including detecting emergence and spread of rare variants, establishing transmission sources and inferring changes in malaria transmission intensity. This review will potentially guide the collection of samples and data, serve as a starting point for further methodological innovation, and enhance utilization of MMS to support malaria elimination.

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

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

The immunity modulation of transforming growth factor-β in malaria and other pathological process

The main causative agent of malaria in humans is Plasmodium falciparum, which is spread through biting Anopheles mosquitoes. Immunoregulation in the host involving the pleiotropic cytokine transforming growth factor-β (TGF-β) has a vital role in controlling the immune response to P. falciparum infection. Based on a search of the published literature, this study investigated the correlation between malaria and immune cells, specifically the role of TGF-β in the immune response. The studies analyzed showed that, when present in low amounts, TGF-β promotes inflammation, but inhibits inflammation when present in high concentrations; thus, it is an essential regulator of inflammation. It has also been shown that the quantity of TGF-β produced by the host can influence how badly the parasite affects the host. Low levels of TGF-β in the host prevent the host from being able to manage the inflammation that Plasmodium causes, which results in a pathological situation that leaves the host vulnerable to fatal infection. Additionally, the amount of TGF-β fluctuates throughout the host's Plasmodium infection. At the beginning of a Plasmodium infection, TGF-β levels are noticeably increased, and as Plasmodium multiplies quickly, they start to decline, hindering further growth. In addition, it is also involved in the growth, proliferation, and operation of various types of immune cell and correlated with levels of cytokines associated with the immune response to malaria. TGF-β levels were positively connected with the anti-inflammatory cytokine interleukin-10 (IL-10), but negatively correlated with the proinflammatory cytokines interferon-γ (IFN-γ) and IL-6 in individuals with severe malaria. Thus, TGF-β might balance immune-mediated pathological damage and the regulation and clearance of infectious pathogens. Numerous domestic and international studies have demonstrated that TGF-β maintains a dynamic balance between anti-inflammation and pro-inflammation in malaria immunity by acting as an anti-inflammatory factor when inflammation levels are too high and as a pro-inflammatory factor when inflammation levels are deficient. Such information could be of relevance to the design of urgently needed vaccines and medications to meet the emerging risks associated with the increasing spread of malaria and the development of drug resistance.

Molecular surveillance of Plasmodium falciparum drug-resistance markers in Vietnam using multiplex amplicon sequencing (2000-2016)

Emergence and spread of Plasmodium falciparum resistance to artemisinin-based combination therapies (ACT) is a major challenge for Greater Mekong Subregion countries in their goal to eliminate malaria by 2030. Tools to efficiently monitor drug resistance beyond resource-demanding therapeutic efficacy studies are necessary. A custom multiplex amplicon sequencing assay based on Illumina technology was designed to target the marker of partial resistance to artemisinin (K13), five candidate modulators of artemisinin resistance, the marker of resistance to chloroquine (crt), and four neutral microsatellite loci. The assay was used to genotype 635 P. falciparum-positive blood samples collected across seven provinces of Vietnam and one of Cambodia between 2000 and 2016. Markers of resistance to artemisinin partner-drugs piperaquine (copy number of plasmepsin-2) and mefloquine (copy number of multidrug-resistance 1) were determined by qPCR. Parasite population structure was further assessed using a 101-SNP barcode. Validated mutations of artemisinin partial resistance in K13 were found in 48.1% of samples, first detection was in 2000, and by 2015 prevalence overcame > 50% in Central Highlands and Binh Phuoc province. K13-C580Y variant became predominant country-wide, quickly replacing an outbreak of K13-I543T in Central Highlands. Mutations in candidate artemisinin resistance modulator genes paralleled the trends of K13 mutants, whereas resistance to piperaquine and mefloquine remained low (≈ 10%) by 2015-2016. Genomic tools applied to malaria surveillance generate comprehensive information on dynamics of drug resistance and population structure and reflect drug efficacy profiles from in vivo studies.

Malaria surveillance reveals parasite relatedness, signatures of selection, and correlates of transmission across Senegal

Parasite genetic surveillance has the potential to play an important role in malaria control. We describe here an analysis of data from the first year of an ongoing, nationwide program of genetic surveillance of Plasmodium falciparum parasites in Senegal, intended to provide actionable information for malaria control efforts. Looking for a good proxy for local malaria incidence, we found that the best predictor was the proportion of polygenomic infections (those with multiple genetically distinct parasites), although that relationship broke down in very low incidence settings (r = 0.77 overall). The proportion of closely related parasites in a site was more weakly correlated ( r = -0.44) with incidence while the local genetic diversity was uninformative. Study of related parasites indicated their potential for discriminating local transmission patterns: two nearby study areas had similarly high fractions of relatives, but one area was dominated by clones and the other by outcrossed relatives. Throughout the country, 58% of related parasites proved to belong to a single network of relatives, within which parasites were enriched for shared haplotypes at known and suspected drug resistance loci as well as at one novel locus, reflective of ongoing selection pressure.

Molecular epidemiology of continued Plasmodium falciparum disease transmission after an outbreak in Ecuador

To better understand the factors underlying the continued incidence of clinical episodes of falciparum malaria in E-2025 countries targeting elimination, we characterized the molecular epidemiology of Plasmodium falciparum disease transmission after a clonal outbreak in Ecuador. Here we study disease transmission by documenting the diversity and population structure of the major variant surface antigen of the blood stages of P. falciparum encoded by the var multigene family. We used a high-resolution genotyping method, "varcoding", involving targeted amplicon sequencing to fingerprint the DBLα encoding region of var genes to describe both antigenic var diversity and var repertoire similarity or relatedness in parasite isolates from clinical cases. We identified nine genetic varcodes in 58 P. falciparum isolates causing clinical disease in 2013-2015. Network analyses revealed that four of the varcodes were highly related to the outbreak varcode, with identification of possible diversification of the outbreak parasites by recombination as seen in three of those varcodes. The majority of clinical cases in Ecuador were associated with parasites with highly related or recombinant varcodes to the outbreak clone and due to local transmission rather than recent importation of parasites from other endemic countries. Sharing of types in Ecuadorian varcodes to those sampled in South American varcodes reflects historical parasite importation of some varcodes, especially from Colombia and Peru. Our findings highlight the translational application of varcoding for outbreak surveillance in epidemic/unstable malaria transmission, such as in E-2025 countries, and point to the need for surveillance of local reservoirs of infection in Ecuador to achieve the malaria elimination goal by 2025.

Malaria Molecular Surveillance in the Peruvian Amazon with a Novel Highly Multiplexed Plasmodium falciparum AmpliSeq Assay

Molecular surveillance for malaria has great potential to support national malaria control programs (NMCPs). To bridge the gap between research and implementation, several applications (use cases) have been identified to align research, technology development, and public health efforts. For implementation at NMCPs, there is an urgent need for feasible and cost-effective tools. We designed a new highly multiplexed deep sequencing assay (Pf AmpliSeq), which is compatible with benchtop sequencers, that allows high-accuracy sequencing with higher coverage and lower cost than whole-genome sequencing (WGS), targeting genomic regions of interest. The novelty of the assay is its high number of targets multiplexed into one easy workflow, combining population genetic markers with 13 nearly full-length resistance genes, which is applicable for many different use cases. We provide the first proof of principle for hrp2 and hrp3 deletion detection using amplicon sequencing. Initial sequence data processing can be performed automatically, and subsequent variant analysis requires minimal bioinformatic skills using any tabulated data analysis program. The assay was validated using a retrospective sample collection (n = 254) from the Peruvian Amazon between 2003 and 2018. By combining phenotypic markers and a within-country 28-single-nucleotide-polymorphism (SNP) barcode, we were able to distinguish different lineages with multiple resistance haplotypes (in dhfr, dhps, crt and mdr1) and hrp2 and hrp3 deletions, which have been increasing in recent years. We found no evidence to suggest the emergence of artemisinin (ART) resistance in Peru. These findings indicate a parasite population that is under drug pressure but is susceptible to current antimalarials and demonstrate the added value of a highly multiplexed molecular tool to inform malaria strategies and surveillance systems. IMPORTANCE While the power of next-generation sequencing technologies to inform and guide malaria control programs has become broadly recognized, the integration of genomic data for operational incorporation into malaria surveillance remains a challenge in most countries where malaria is endemic. The main obstacles include limited infrastructure, limited access to high-throughput sequencing facilities, and the need for local capacity to run an in-country analysis of genomes at a large-enough scale to be informative for surveillance. In addition, there is a lack of standardized laboratory protocols and automated analysis pipelines to generate reproducible and timely results useful for relevant stakeholders. With our standardized laboratory and bioinformatic workflow, malaria genetic surveillance data can be readily generated by surveillance researchers and malaria control programs in countries of endemicity, increasing ownership and ensuring timely results for informed decision- and policy-making.

Measurably recombining malaria parasites

Genomic epidemiology has guided research and policy for various viral pathogens and there has been a parallel effort towards using genomic epidemiology to combat diseases that are caused by eukaryotic pathogens, such as the malaria parasite. However, the central concept of viral genomic epidemiology, namely that of measurably mutating pathogens, does not apply easily to sexually recombining parasites. Here we introduce the related but different concept of measurably recombining malaria parasites to promote convergence around a unifying theoretical framework for malaria genomic epidemiology. Akin to viral phylodynamics, we anticipate that an inferential framework developed around recombination will help guide practical research and thus realize the full public health potential of genomic epidemiology for malaria parasites and other sexually recombining pathogens.

Resolving drug selection and migration in an inbred South American Plasmodium falciparum population with identity-by-descent analysis

The human malaria parasite Plasmodium falciparum is globally widespread, but its prevalence varies significantly between and even within countries. Most population genetic studies in P. falciparum focus on regions of high transmission where parasite populations are large and genetically diverse, such as sub-Saharan Africa. Understanding population dynamics in low transmission settings, however, is of particular importance as these are often where drug resistance first evolves. Here, we use the Pacific Coast of Colombia and Ecuador as a model for understanding the population structure and evolution of Plasmodium parasites in small populations harboring less genetic diversity. The combination of low transmission and a high proportion of monoclonal infections means there are few outcrossing events and clonal lineages persist for long periods of time. Yet despite this, the population is evolutionarily labile and has successfully adapted to changes in drug regime. Using newly sequenced whole genomes, we measure relatedness between 166 parasites, calculated as identity by descent (IBD), and find 17 distinct but highly related clonal lineages, six of which have persisted in the region for at least a decade. This inbred population structure is captured in more detail with IBD than with other common population structure analyses like PCA, ADMIXTURE, and distance-based trees. We additionally use patterns of intra-chromosomal IBD and an analysis of haplotypic variation to explore past selection events in the region. Two genes associated with chloroquine resistance, crt and aat1, show evidence of hard selective sweeps, while selection appears soft and/or incomplete at three other key resistance loci (dhps, mdr1, and dhfr). Overall, this work highlights the strength of IBD analyses for studying parasite population structure and resistance evolution in regions of low transmission, and emphasizes that drug resistance can evolve and spread in small populations, as will occur in any region nearing malaria elimination.

Drug resistance of Plasmodium falciparum and Plasmodium vivax isolates in Indonesia

This review article aims to investigate the genotypic profiles of Plasmodium falciparum and Plasmodium vivax isolates collected across a wide geographic region and their association with resistance to anti-malarial drugs used in Indonesia. A systematic review was conducted between 1991 and date. Search engines, such as PubMed, Science Direct, and Google Scholar, were used for articles published in English and Indonesian to search the literature. Of the 471 initially identified studies, 61 were selected for 4316 P. falciparum and 1950 P. vivax individual infections. The studies included 23 molecular studies and 38 therapeutic efficacy studies. K76T was the most common pfcrt mutation. K76N (2.1%) was associated with the haplotype CVMNN. By following dihydroartemisinin-piperaquine (DHA-PPQ) therapy, the mutant pfmdr1 alleles 86Y and 1034C were selected. Low prevalence of haplotype N86Y/Y184/D1246Y pfmdr1 reduces susceptibility to AS-AQ. SNP mutation pvmdr1 Y976F reached 96.1% in Papua and East Nusa Tenggara. Polymorphism analysis in the pfdhfr gene revealed 94/111 (84.7%) double mutants S108N/C59R or S108T/A16V in Central Java. The predominant pfdhfr haplotypes (based on alleles 16, 51, 59,108, 164) found in Indonesia were ANCNI, ANCSI, ANRNI, and ANRNL. Some isolates carried A437G (35.3%) or A437G/K540E SNPs (26.5%) in pfdhps. Two novel pfdhps mutant alleles, I588F/G and K540T, were associated with six pfdhps haplotypes. The highest prevalence of pvdhfr quadruple mutation (F57L/S58R/T61M/S117T) (61.8%) was detected in Papua. In pvdhps, the only polymorphism before and after 2008 was 383G mutation with 19% prevalence. There were no mutations in the pfk13 gene reported with validated and candidate or associated k13 mutation. An increased copy number of pfpm2, associated with piperaquine resistance, was found only in cases of reinfection. Meanwhile, mutation of pvk12 and pvpm4 I165V is unlikely associated with ART and PPQ drug resistance. DHA-PPQ is still effective in treating uncomplicated falciparum and vivax malaria. Serious consideration should be given to interrupt local malaria transmission and dynamic patterns of resistance to anti-malarial drugs to modify chemotherapeutic policy treatment strategies. The presence of several changes in pfk13 in the parasite population is of concern and highlights the importance of further evaluation of parasitic ART susceptibility in Indonesia.

Genomic approaches for monitoring transmission dynamics of malaria: A case for malaria molecular surveillance in Sub-Saharan Africa

Transmission dynamics is an important indicator for malaria control and elimination. As we move closer to eliminating malaria in Sub-Saharan Africa (sSA), transmission indices with higher resolution (genomic approaches) will complement our current measurements of transmission. Most of the present programmatic knowledge of malaria transmission patterns are derived from assessments of epidemiologic and clinical data, such as case counts, parasitological estimates of parasite prevalence, and Entomological Inoculation Rates (EIR). However, to eliminate malaria from endemic areas, we need to track changes in the parasite population and how they will impact transmission. This is made possible through the evolving field of genomics and genetics, as well as the development of tools for more in-depth studies on the diversity of parasites and the complexity of infections, among other topics. If malaria elimination is to be achieved globally, country-specific elimination activities should be supported by parasite genomic data from regularly collected blood samples for diagnosis, surveillance and possibly from other programmatic interventions. This presents a unique opportunity to track the spread of malaria parasites and shed additional light on intervention efficacy. In this review, various genetic techniques are highlighted along with their significance for an enhanced understanding of transmission patterns in distinct topological settings throughout Sub-Saharan Africa. The importance of these methods and their limitations in malaria surveillance to guide control and elimination strategies, are explored.

Temporal dynamics of Plasmodium falciparum population in Metehara, east-central Ethiopia

BACKGROUND

Plasmodium falciparum is the most serious, genetically most complex and fastest-evolving malaria parasite. Information on genetic diversity of this parasite would guide policy decision and malaria elimination endeavors. This study explored the temporal dynamics of P. falciparum population in two time points in Metehara, east-central Ethiopia.

METHODS

The participants were quantitative real-time polymerase chain reaction-confirmed patients who were recruited for uncomplicated falciparum malaria therapeutic efficacy test in 2015 and 2019. Dry blood spot samples were analysed by the nested PCR to genotype P. falciparum merozoite surface protein (msp1, msp2) and glutamate-rich protein (glurp) genes.

RESULTS

While msp1, msp2 and glurp genotypes were successfully detected in 26(89.7%), 24(82.8%) and 14(48.3%) of 2015 samples (n = 29); the respective figures for 2019 (n = 41) were 31(68.3%), 39(95.1%), 25(61.0%). In 2015, the frequencies of K1, MAD20 and RO33 allelic families of msp1, and FC27 and IC/3D7 of msp2 were 19(73.1%), 8(30.6%), 14(53.8%), 21(87.5%), 12(50.5%); and in 2019 it was 15(48.4%), 19(61.3%), 15(48.4%), 30(76.9%), 27(69.2%) respectively. MAD20 has shown dominance over both K1 and RO33 in 2019 compared to the proportion in 2015. Similarly, although FC27 remained dominant, there was shifting trend in the frequency of IC/3D7 from 50.5% in 2015 to 69.2% in 2019. The multiplicity of infection (MOI) and expected heterozygosity index (He) in 2015 and 2019 were respectively [1.43 ± 0.84] and [1.15 ± 0.91], 0.3 and 0.03 for msp1. However, there was no significant association between MOI and age or parasitaemia in both time points.

CONCLUSION

The lower genetic diversity in P. falciparum population in the two time points and overall declining trend as demonstrated by the lower MOI and He may suggest better progress in malaria control in Metehara. But, the driving force and selective advantage of switching to MAD20 dominance over the other two msp1 allelic families, and the dynamics within msp2 alleles needs further investigation.

Novel highly-multiplexed AmpliSeq targeted assay for Plasmodium vivax genetic surveillance use cases at multiple geographical scales

Although the power of genetic surveillance tools has been acknowledged widely, there is an urgent need in malaria endemic countries for feasible and cost-effective tools to implement in national malaria control programs (NMCPs) that can generate evidence to guide malaria control and elimination strategies, especially in the case of Plasmodium vivax. Several genetic surveillance applications ('use cases') have been identified to align research, technology development, and public health efforts, requiring different types of molecular markers. Here we present a new highly-multiplexed deep sequencing assay (Pv AmpliSeq). The assay targets the 33-SNP vivaxGEN-geo panel for country-level classification, and a newly designed 42-SNP within-country barcode for analysis of parasite dynamics in Vietnam and 11 putative drug resistance genes in a highly multiplexed NGS protocol with easy workflow, applicable for many different genetic surveillance use cases. The Pv AmpliSeq assay was validated using: 1) isolates from travelers and migrants in Belgium, and 2) routine collections of the national malaria control program at sentinel sites in Vietnam. The assay targets 229 amplicons and achieved a high depth of coverage (mean 595.7 ± 481) and high accuracy (mean error-rate of 0.013 ± 0.007). P. vivax parasites could be characterized from dried blood spots with a minimum of 5 parasites/µL and 10% of minority-clones. The assay achieved good spatial specificity for between-country prediction of origin using the 33-SNP vivaxGEN-geo panel that targets rare alleles specific for certain countries and regions. A high resolution for within-country diversity in Vietnam was achieved using the designed 42-SNP within-country barcode that targets common alleles (median MAF 0.34, range 0.01-0.49. Many variants were detected in (putative) drug resistance genes, with different predominant haplotypes in the pvmdr1 and pvcrt genes in different provinces in Vietnam. The capacity of the assay for high resolution identity-by-descent (IBD) analysis was demonstrated and identified a high rate of shared ancestry within Gia Lai Province in the Central Highlands of Vietnam, as well as between the coastal province of Binh Thuan and Lam Dong. Our approach performed well in geographically differentiating isolates at multiple spatial scales, detecting variants in putative resistance genes, and can be easily adjusted to suit the needs in other settings in a country or region. We prioritize making this tool available to researchers and NMCPs in endemic countries to increase ownership and ensure data usage for decision-making and malaria policy.

Receptivity to malaria: meaning and measurement

"Receptivity" to malaria is a construct developed during the Global Malaria Eradication Programme (GMEP) era. It has been defined in varied ways and no consistent, quantitative definition has emerged over the intervening decades. Despite the lack of consistency in defining this construct, the idea that some areas are more likely to sustain malaria transmission than others has remained important in decision-making in malaria control, planning for malaria elimination and guiding activities during the prevention of re-establishment (POR) period. This manuscript examines current advances in methods of measurement. In the context of a decades long decline in global malaria transmission and an increasing number of countries seeking to eliminate malaria, understanding and measuring malaria receptivity has acquired new relevance.

Integrating genomic and epidemiologic data to accelerate progress toward schistosomiasis elimination

The global community has adopted ambitious goals to eliminate schistosomiasis as a public health problem, and new tools are needed to achieve them. Mass drug administration programs, for example, have reduced the burden of schistosomiasis, but the identification of hotspots of persistent and reemergent transmission threaten progress toward elimination and underscore the need to couple treatment with interventions that reduce transmission. Recent advances in DNA sequencing technologies make whole-genome sequencing a valuable and increasingly feasible option for population-based studies of complex parasites such as schistosomes. Here, we focus on leveraging genomic data to tailor interventions to distinct social and ecological circumstances. We consider two priority questions that can be addressed by integrating epidemiological, ecological, and genomic information: (1) how often do non-human host species contribute to human schistosome infection? and (2) what is the importance of locally acquired versus imported infections in driving transmission at different stages of elimination? These questions address processes that can undermine control programs, especially those that rely heavily on treatment with praziquantel. Until recently, these questions were difficult to answer with sufficient precision to inform public health decision-making. We review the literature related to these questions and discuss how whole-genome approaches can identify the geographic and taxonomic sources of infection, and how such information can inform context-specific efforts that advance schistosomiasis control efforts and minimize the risk of reemergence.

Genetic surveillance for monitoring the impact of drug use on Plasmodium falciparum populations

The use of antimalarial drugs is an effective strategy in the fight against malaria. However, selection of drug resistant parasites is a constant threat to the continued use of this approach. Antimalarial drugs are used not only to treat infections but also as part of population-level strategies to reduce malaria transmission toward elimination. While there is strong evidence that the ongoing use of antimalarial drugs increases the risk of the emergence and spread of drug-resistant parasites, it is less clear how population-level use of drug-based interventions like seasonal malaria chemoprevention (SMC) or mass drug administration (MDA) may contribute to drug resistance or loss of drug efficacy. Critical to sustained use of drug-based strategies for reducing the burden of malaria is the surveillance of population-level signals related to transmission reduction and resistance selection. Here we focus on Plasmodium falciparum and discuss the genetic signatures of a parasite population that are correlated with changes in transmission and related to drug pressure and resistance as a result of drug use. We review the evidence for MDA and SMC contributing to malaria burden reduction and drug resistance selection and examine the use and impact of these interventions in Senegal. Throughout we consider best strategies for ongoing surveillance of both population and resistance signals in the context of different parasite population parameters. Finally, we propose a roadmap for ongoing surveillance during population-level drug-based interventions to reduce the global malaria burden.

Time to scale up molecular surveillance for anti-malarial drug resistance in sub-saharan Africa

Artemisinin resistance has emerged and spread in the Greater Mekong Sub-region (GMS), followed by artemisinin-based combination therapy failure, due to both artemisinin and partner drug resistance. More worrying, artemisinin resistance has been recently reported and confirmed in Rwanda. Therefore, there is an urgent need to strengthen surveillance systems beyond the GMS to track the emergence or spread of artemisinin and partner drug resistance in other endemic settings. Currently, anti-malarial drug efficacy is monitored primarily through therapeutic efficacy studies (TES). Even though essential for anti-malarial drug policy change, these studies are difficult to conduct, expensive, and may not detect the early emergence of resistance. Additionally, results from TES may take years to be available to the stakeholders, jeopardizing their usefulness. Molecular markers are additional and useful tools to monitor anti-malarial drug resistance, as samples collected on dried blood spots are sufficient to monitor known and validated molecular markers of resistance, and could help detecting and monitoring the early emergence of resistance. However, molecular markers are not monitored systematically by national malaria control programmes, and are often assessed in research studies, but not in routine surveillance. The implementation of molecular markers as a routine tool for anti-malarial drug resistance surveillance could greatly improve surveillance of anti-malarial drug efficacy, making it possible to detect resistance before it translates to treatment failures. When possible, ex vivo assays should be included as their data could be useful complementary, especially when no molecular markers are validated.

Molecular surveillance of malaria scales up

Parasite and vector genetic data can guide malaria control, and technological advances are enabling more informative genetic data generation at unprecedented scales. Jacob et al. employ multiplexed amplicon sequencing to profile parasite genetic diversity from thousands of malaria samples, illuminating spatiotemporal patterns of drug resistance to inform regional drug policy change.

Advances and opportunities in malaria population genomics

Almost 20 years have passed since the first reference genome assemblies were published for Plasmodium falciparum, the deadliest malaria parasite, and Anopheles gambiae, the most important mosquito vector of malaria in sub-Saharan Africa. Reference genomes now exist for all human malaria parasites and nearly half of the ~40 important vectors around the world. As a foundation for genetic diversity studies, these reference genomes have helped advance our understanding of basic disease biology and drug and insecticide resistance, and have informed vaccine development efforts. Population genomic data are increasingly being used to guide our understanding of malaria epidemiology, for example by assessing connectivity between populations and the efficacy of parasite and vector interventions. The potential value of these applications to malaria control strategies, together with the increasing diversity of genomic data types and contexts in which data are being generated, raise both opportunities and challenges in the field. This Review discusses advances in malaria genomics and explores how population genomic data could be harnessed to further support global disease control efforts.

Linked surveillance and genetic data uncovers programmatically relevant geographic scale of Guinea worm transmission in Chad

Background

Guinea worm (Dracunculus medinensis) was detected in Chad in 2010 after a supposed ten-year absence, posing a challenge to the global eradication effort. Initiation of a village-based surveillance system in 2012 revealed a substantial number of dogs infected with Guinea worm, raising questions about paratenic hosts and cross-species transmission.

Methodology/principal findings

We coupled genomic and surveillance case data from 2012-2018 to investigate the modes of transmission between dog and human hosts and the geographic connectivity of worms. Eighty-six variants across four genes in the mitochondrial genome identified 41 genetically distinct worm genotypes. Spatiotemporal modeling revealed worms with the same genotype (‘genetically identical’) were within a median range of 18.6 kilometers of each other, but largely within approximately 50 kilometers. Genetically identical worms varied in their degree of spatial clustering, suggesting there may be different factors that favor or constrain transmission. Each worm was surrounded by five to ten genetically distinct worms within a 50 kilometer radius. As expected, we observed a change in the genetic similarity distribution between pairs of worms using variants across the complete mitochondrial genome in an independent population.

Conclusions/significance

In the largest study linking genetic and surveillance data to date of Guinea worm cases in Chad, we show genetic identity and modeling can facilitate the understanding of local transmission. The co-occurrence of genetically non-identical worms in quantitatively identified transmission ranges highlights the necessity for genomic tools to link cases. The improved discrimination between pairs of worms from variants identified across the complete mitochondrial genome suggests that expanding the number of genomic markers could link cases at a finer scale. These results suggest that scaling up genomic surveillance for Guinea worm may provide additional value for programmatic decision-making critical for monitoring cases and intervention efficacy to achieve elimination.

The global eradication effort for Guinea worm disease has dramatically decreased the global burden of the disease and enabled 187 countries to be certified by the World Health Organization to be free of endemic transmission. Despite this progress, several countries continue to have endemic transmission. In Chad, a long absence of reported cases was interrupted with the identification of new Guinea worm cases, prompting a substantial scale up of surveillance and intervention efforts. Here, we study the value of increasing genomic surveillance as a tool for programmatic evaluation of surveillance and intervention efforts in Chad. Linking surveillance and genomic samples, parsimonious spatial models help reveal a consistent geographic clustering of similar genetic sequences across Chad. We also demonstrate that expanding the sequencing can offer better resolution for distinguishing Guinea worm samples. In this retrospective study, we found evidence that scaling up genomic surveillance can be an important monitoring and evaluation tool for the eradication program in Chad.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

Antimalarial Drug Resistance and Implications for the WHO Global Technical Strategy

PURPOSE OF REVIEW

Five years have passed since the World Health Organization released its Global Technical Strategy for Malaria (GTS). In that time, progress against malaria has plateaued. This review focuses on the implications of antimalarial drug resistance for the GTS and how interim progress in parasite genomics and antimalarial pharmacology offer a bulwark against it.

RECENT FINDINGS

For the first time, drug resistance-conferring genes have been identified and validated before their global expansion in malaria parasite populations. More efficient methods for their detection and elaboration have been developed, although low-density infections and polyclonality remain a nuisance to be solved. Clinical trials of alternative regimens for multidrug-resistant malaria have delivered promising results. New agents continue down the development pipeline, while a nascent infrastructure in sub-Saharan Africa for conducting phase I trials and trials of transmission-blocking agents has come to fruition after years of preparation.

SUMMARY

These and other developments can help inform the GTS as the world looks ahead to the next two decades of its implementation. To remain ahead of the threat that drug resistance poses, wider application of genomic-based surveillance and optimization of existing and forthcoming antimalarial drugs are essential.

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.

Multi-locus genotyping reveals established endemicity of a geographically distinct Plasmodium vivax population in Mauritania, West Africa

BACKGROUND

Plasmodium vivax has been recently discovered as a significant cause of malaria in Mauritania, although very rare elsewhere in West Africa. It has not been known if this is a recently introduced or locally remnant parasite population, nor whether the genetic structure reflects epidemic or endemic transmission.

METHODOLOGY/PRINCIPAL FINDINGS

To investigate the P. vivax population genetic structure in Mauritania and compare with populations previously analysed elsewhere, multi-locus genotyping was undertaken on 100 clinical isolates, using a genome-wide panel of 38 single nucleotide polymorphisms (SNPs), plus seven SNPs in drug resistance genes. The Mauritanian P. vivax population is shown to be genetically diverse and divergent from populations elsewhere, indicated consistently by genetic distance matrix analysis, principal components analyses, and fixation indices. Only one isolate had a genotype clearly indicating recent importation, from a southeast Asian source. There was no linkage disequilibrium in the local parasite population, and only a small number of infections appeared to be closely genetically related, indicating that there is ongoing genetic recombination consistent with endemic transmission. The P. vivax diversity in a remote mining town was similar to that in the capital Nouakchott, with no indication of local substructure or of epidemic population structure. Drug resistance alleles were virtually absent in Mauritania, in contrast with P. vivax in other areas of the world.

CONCLUSIONS/SIGNIFICANCE

The molecular epidemiology indicates that there is long-standing endemic transmission that will be very challenging to eliminate. The virtual absence of drug resistance alleles suggests that most infections have been untreated, and that this endemic infection has been more neglected in comparison to P. vivax elsewhere.

Describing the current status of Plasmodium falciparum population structure and drug resistance within mainland Tanzania using molecular inversion probes

High-throughput Plasmodium genomic data is increasingly useful in assessing prevalence of clinically important mutations and malaria transmission patterns. Understanding parasite diversity is important for identification of specific human or parasite populations that can be targeted by control programmes, and to monitor the spread of mutations associated with drug resistance. An up-to-date understanding of regional parasite population dynamics is also critical to monitor the impact of control efforts. However, this data is largely absent from high-burden nations in Africa, and to date, no such analysis has been conducted for malaria parasites in Tanzania countrywide. To this end, over 1,000 P. falciparum clinical isolates were collected in 2017 from 13 sites in seven administrative regions across Tanzania, and parasites were genotyped at 1,800 variable positions genome-wide using molecular inversion probes. Population structure was detectable among Tanzanian P. falciparum parasites, approximately separating parasites from the northern and southern districts and identifying genetically admixed populations in the north. Isolates from nearby districts were more likely to be genetically related compared to parasites sampled from more distant districts. Known drug resistance mutations were seen at increased frequency in northern districts (including two infections carrying pfk13-R561H), and additional variants with undetermined significance for antimalarial resistance also varied by geography. Malaria Indicator Survey (2017) data corresponded with genetic findings, including average region-level complexity-of-infection and malaria prevalence estimates. The parasite populations identified here provide important information on extant spatial patterns of genetic diversity of Tanzanian parasites, to which future surveys of genetic relatedness can be compared.

Identity-by-descent with uncertainty characterises connectivity of Plasmodium falciparum populations on the Colombian-Pacific coast

Characterising connectivity between geographically separated biological populations is a common goal in many fields. Recent approaches to understanding connectivity between malaria parasite populations, with implications for disease control efforts, have used estimates of relatedness based on identity-by-descent (IBD). However, uncertainty around estimated relatedness has not been accounted for. IBD-based relatedness estimates with uncertainty were computed for pairs of monoclonal Plasmodium falciparum samples collected from five cities on the Colombian-Pacific coast where long-term clonal propagation of P. falciparum is frequent. The cities include two official ports, Buenaventura and Tumaco, that are separated geographically but connected by frequent marine traffic. Fractions of highly-related sample pairs (whose classification using a threshold accounts for uncertainty) were greater within cities versus between. However, based on both highly-related fractions and on a threshold-free approach (Wasserstein distances between parasite populations) connectivity between Buenaventura and Tumaco was disproportionally high. Buenaventura-Tumaco connectivity was consistent with transmission events involving parasites from five clonal components (groups of statistically indistinguishable parasites identified under a graph theoretic framework). To conclude, P. falciparum population connectivity on the Colombian-Pacific coast abides by accessibility not isolation-by-distance, potentially implicating marine traffic in malaria transmission with opportunities for targeted intervention. Further investigations are required to test this hypothesis. For the first time in malaria epidemiology (and to our knowledge in ecological and epidemiological studies more generally), we account for uncertainty around estimated relatedness (an important consideration for studies that plan to use genotype versus whole genome sequence data to estimate IBD-based relatedness); we also use threshold-free methods to compare parasite populations and identify clonal components. Threshold-free methods are especially important in analyses of malaria parasites and other recombining organisms with mixed mating systems where thresholds do not have clear interpretation (e.g. due to clonal propagation) and thus undermine the cross-comparison of studies.

Changes in the frequencies of Plasmodium falciparum dhps and dhfr drug-resistant mutations in children from Western Kenya from 2005 to 2018: the rise of Pfdhps S436H

BACKGROUND

Sulfadoxine-pyrimethamine (SP) is the only anti-malarial drug formulation approved for intermittent preventive treatment in pregnancy (IPTp). However, mutations in the Plasmodium falciparum dhfr (Pfdhfr) and dhps (Pfdhps) genes confer resistance to pyrimethamine and sulfadoxine, respectively. Here, the frequencies of SP resistance-associated mutations from 2005 to 2018 were compared in samples from Kenyan children with malaria residing in a holoendemic transmission region.

METHODS

Partial sequences of the Pfdhfr and Pfdhps genes were amplified and sequenced from samples collected in 2005 (n = 81), 2010 (n = 95), 2017 (n = 43), and 2018 (n = 55). The frequency of known mutations conferring resistance to pyrimethamine and sulfadoxine were estimated and compared. Since artemisinin-based combination therapy (ACT) is the current first-line treatment for malaria, the presence of mutations in the propeller domain of P. falciparum kelch13 gene (Pfk13) linked to ACT-delayed parasite clearance was studied in the 2017/18 samples.

RESULTS

Among other changes, the point mutation of Pfdhps S436H increased in frequency from undetectable in 2005 to 28% in 2017/18. Triple Pfdhfr mutant allele (CIRNI) increased in frequency from 84% in 2005 to 95% in 2017/18, while the frequency of Pfdhfr double mutant alleles declined (allele CICNI from 29% in 2005 to 6% in 2017/18, and CNRNI from 9% in 2005 to undetectable in 2010 and 2017/18). Thus, a multilocus Pfdhfr/Pfdhps genotype with six mutations (HGEAA/CIRNI), including Pfdhps S436H, increased in frequency from 2010 to 2017/18. Although none of the mutations associated with ACT-delayed parasite clearance was observed, the Pfk13 mutation A578S, the most widespread Pfk13 SNP found in Africa, was detected in low frequency (2.04%).

CONCLUSIONS

There were changes in SP resistance mutant allele frequencies, including an increase in the Pfdhps S436H. Although these patterns seem consistent with directional selection due to drug pressure, there is a lack of information to determine the actual cause of such changes. These results suggest incorporating molecular surveillance of Pfdhfr/Pfdhps mutations in the context of SP efficacy studies for intermittent preventive treatment in pregnancy (IPTp).

Ultra-sensitive RDT performance and antigen dynamics in a high-transmission Plasmodium falciparum setting in Mali

BACKGROUND

The recent expansion of tools designed to accurately quantify malaria parasite-produced antigens has enabled us to evaluate the performance of rapid diagnostic tests (RDTs) as a function of the antigens they detect-typically histidine rich protein 2 (HRP2) or lactate dehydrogenase (LDH).

METHODS

For this analysis, whole blood specimens from a longitudinal study in Bancoumana, Mali were used to evaluate the performance of the ultra-sensitive HRP2-based Alere™ Malaria Ag P.f RDT (uRDT). The samples were collected as part of a transmission-blocking vaccine trial in a high transmission region for Plasmodium falciparum malaria. Furthermore, antigen dynamics after successful anti-malarial drug treatment were evaluated in these samples using the Q-Plex Human Malaria Array (4-Plex) to quantify antigen concentrations.

RESULTS

The uRDT had a 50% probability of a positive result at 207 pg/mL HRP2 [95% credible interval (CrI) 160-268]. Individuals with symptomatic infection remained positive by uRDT for a median of 33 days [95% confidence interval (CI) 28-47] post anti-malarial drug treatment. Biphasic exponential decay models accurately captured the population level post-treatment dynamics of both HRP2 and Plasmodium LDH (pLDH), with the latter decaying more rapidly. Motivated by these differences in rates of decay, a novel algorithm that used HRP2:pLDH ratios to predict if an individual had active versus recently cleared P. falciparum infection was developed. The algorithm had 77.5% accuracy in correctly classifying antigen-positive individuals as those with and without active infection.

CONCLUSIONS

These results characterize the performance of the ultra-sensitive RDT and demonstrate the potential for emerging antigen-quantifying technologies in the field of malaria diagnostics to be helpful tools in distinguishing between active versus recently cleared malaria infections.

Sensitive, highly multiplexed sequencing of microhaplotypes from the Plasmodium falciparum heterozygome

Background

Targeted next-generation sequencing offers the potential for consistent, deep coverage of information-rich genomic regions to characterize polyclonal Plasmodium falciparum infections. However, methods to identify and sequence these genomic regions are currently limited.

Methods

A bioinformatic pipeline and multiplex methods were developed to identify and simultaneously sequence 100 targets and applied to dried blood spot (DBS) controls and field isolates from Mozambique. For comparison, whole-genome sequencing data were generated for the same controls.

Results

Using publicly available genomes, 4465 high-diversity genomic regions suited for targeted sequencing were identified, representing the P. falciparum heterozygome. For this study, 93 microhaplotypes with high diversity (median expected heterozygosity = 0.7) were selected along with 7 drug resistance loci. The sequencing method achieved very high coverage (median 99%), specificity (99.8%), and sensitivity (90% for haplotypes with 5% within sample frequency in dried blood spots with 100 parasites/µL). In silico analyses revealed that microhaplotypes provided much higher resolution to discriminate related from unrelated polyclonal infections than biallelic single-nucleotide polymorphism barcodes.

Conclusions

The bioinformatic and laboratory methods outlined here provide a flexible tool for efficient, low-cost, high-throughput interrogation of the P. falciparum genome, and can be tailored to simultaneously address multiple questions of interest in various epidemiological settings.

Genetic evidence for imported malaria and local transmission in Richard Toll, Senegal

BACKGROUND

Malaria elimination efforts can be undermined by imported malaria infections. Imported infections are classified based on travel history.

METHODS

A genetic strategy was applied to better understand the contribution of imported infections and to test for local transmission in the very low prevalence region of Richard Toll, Senegal.

RESULTS

Genetic relatedness analysis, based upon molecular barcode genotyping data derived from diagnostic material, provided evidence for both imported infections and ongoing local transmission in Richard Toll. Evidence for imported malaria included finding that a large proportion of Richard Toll parasites were genetically related to parasites from Thiès, Senegal, a region of moderate transmission with extensive available genotyping data. Evidence for ongoing local transmission included finding parasites of identical genotype that persisted across multiple transmission seasons as well as enrichment of highly related infections within the households of non-travellers compared to travellers.

CONCLUSIONS

These data indicate that, while a large number of infections may have been imported, there remains ongoing local malaria transmission in Richard Toll. These proof-of-concept findings underscore the value of genetic data to identify parasite relatedness and patterns of transmission to inform optimal intervention selection and placement.

Implementing parasite genotyping into national surveillance frameworks: feedback from control programmes and researchers in the Asia-Pacific region

The Asia-Pacific region faces formidable challenges in achieving malaria elimination by the proposed target in 2030. Molecular surveillance of Plasmodium parasites can provide important information on malaria transmission and adaptation, which can inform national malaria control programmes (NMCPs) in decision-making processes. In November 2019 a parasite genotyping workshop was held in Jakarta, Indonesia, to review molecular approaches for parasite surveillance and explore ways in which these tools can be integrated into public health systems and inform policy. The meeting was attended by 70 participants from 8 malaria-endemic countries and partners of the Asia Pacific Malaria Elimination Network. The participants acknowledged the utility of multiple use cases for parasite genotyping including: quantifying the prevalence of drug resistant parasites, predicting risks of treatment failure, identifying major routes and reservoirs of infection, monitoring imported malaria and its contribution to local transmission, characterizing the origins and dynamics of malaria outbreaks, and estimating the frequency of Plasmodium vivax relapses. However, the priority of each use case varies with different endemic settings. Although a one-size-fits-all approach to molecular surveillance is unlikely to be applicable across the Asia-Pacific region, consensus on the spectrum of added-value activities will help support data sharing across national boundaries. Knowledge exchange is needed to establish local expertise in different laboratory-based methodologies and bioinformatics processes. Collaborative research involving local and international teams will help maximize the impact of analytical outputs on the operational needs of NMCPs. Research is also needed to explore the cost-effectiveness of genetic epidemiology for different use cases to help to leverage funding for wide-scale implementation. Engagement between NMCPs and local researchers will be critical throughout this process.

Malaria in Venezuela: changes in the complexity of infection reflects the increment in transmission intensity

Background

Malaria incidence has reached staggering numbers in Venezuela. Commonly, Bolívar State accounted for approximately 70% of the country cases every year. Most cases cluster in the Sifontes municipality, a region characterized by an extractive economy, including gold mining. An increase in migration to Sifontes, driven by gold mining, fueled a malaria spillover to the rest of the country and the region. Here samples collected in 2018 were compared with a previous study of 2003/2004 to describe changes in the parasites population structures and the frequency of point mutations linked to anti-malarial drugs.

Methods

A total of 88 Plasmodium falciparum and 94 Plasmodium vivax isolates were collected in 2018 and compared with samples from 2003/2004 (106 P. falciparum and 104 P. vivax). For P. falciparum, mutations linked to drug resistance (Pfdhfr, Pfdhps, and Pfcrt) and the Pfk13 gene associated with artemisinin delayed parasite clearance, were analysed. To estimate the multiplicity of infection (MOI), and perform P. falciparum and P. vivax population genetic analyses, the parasites were genotyped by using eight standardized microsatellite loci.

Results

The P. falciparum parasites are still harbouring drug-resistant mutations in Pfdhfr, Pfdhps, and Pfcrt. However, there was a decrease in the frequency of highly resistant Pfdhps alleles. Mutations associated with artemisinin delayed parasite clearance in the Pfk13 gene were not found. Consistent with the increase in transmission, polyclonal infections raised from 1.9% in 2003/2004 to 39% in 2018 in P. falciparum and from 16.3 to 68% in P. vivax. There is also a decrease in linkage disequilibrium. Bayesian clustering yields two populations linked to the time of sampling, showing that the parasite populations temporarily changed. However, the samples from 2003/2004 and 2018 have several alleles per locus in common without sharing multi-locus genotypes.

Conclusions

The frequency of mutations linked with drug resistance in P. falciparum shows only changes in Pfdhps. Observations presented here are consistent with an increase in transmission from the previously circulating parasites. Following populations longitudinally, using molecular surveillance, provides valuable information in cases such as Venezuela with a fluid malaria situation that is affecting the regional goals toward elimination.

Molecular surveillance over 14 years confirms reduction of Plasmodium vivax and falciparum transmission after implementation of Artemisinin-based combination therapy in Papua, Indonesia

Genetic epidemiology can provide important insights into parasite transmission that can inform public health interventions. The current study compared long-term changes in the genetic diversity and structure of co-endemic Plasmodium falciparum and P. vivax populations. The study was conducted in Papua Indonesia, where high-grade chloroquine resistance in P. falciparum and P. vivax led to a universal policy of Artemisinin-based Combination Therapy (ACT) in 2006. Microsatellite typing and population genetic analyses were undertaken on available isolates collected between 2004 and 2017 from patients with uncomplicated malaria (n = 666 P. falciparum and n = 615 P. vivax). The proportion of polyclonal P. falciparum infections fell from 28% (38/135) before policy change (2004-2006) to 18% (22/125) at the end of the study (2015-2017); p<0.001. Over the same period, polyclonal P. vivax infections fell from 67% (80/119) to 35% (33/93); p<0.001. P. falciparum strains persisted for up to 9 years compared to 3 months for P. vivax, reflecting higher rates of outbreeding in the latter. Sub-structure was observed in the P. falciparum population, but not in P. vivax, confirming different patterns of outbreeding. The P. falciparum population exhibited 4 subpopulations that changed in frequency over time. Notably, a sharp rise was observed in the frequency of a minor subpopulation (K2) in the late post-ACT period, accounting for 100% of infections in late 2016-2017. The results confirm epidemiological evidence of reduced P. falciparum and P. vivax transmission over time. The smaller change in P. vivax population structure is consistent with greater outbreeding associated with relapsing infections and highlights the need for radical cure to reduce recurrent infections. The study emphasizes the challenge in disrupting P. vivax transmission and demonstrates the potential of molecular data to inform on the impact of public health interventions.

The resistome and genomic reconnaissance in the age of malaria elimination

Malaria is an infectious disease caused by parasitic protozoa in the Plasmodium genus. A complete understanding of the biology of these parasites is challenging in view of their need to switch between the vertebrate and insect hosts. The parasites are also capable of becoming highly motile and of remaining dormant for decades, depending on the stage of their life cycle. Malaria elimination efforts have been implemented in several endemic countries, but the parasites have proven to be resilient. One of the major obstacles for malaria elimination is the development of antimalarial drug resistance. Ineffective treatment regimens will fail to remove the circulating parasites and to prevent the local transmission of the disease. Genomic epidemiology of malaria parasites has become a powerful tool to track emerging drug-resistant parasite populations almost in real time. Population-scale genomic data are instrumental in tracking the hidden pockets of Plasmodium in nationwide elimination efforts. However, genomic surveillance data can be useful in determining the threat only when combined with a thorough understanding of the malarial resistome – the genetic repertoires responsible for causing and potentiating drug resistance evolution. Even though long-term selection has been a standard method for drug target identification in laboratories, its implementation in large-scale exploration of the druggable space in Plasmodium falciparum, along with genome-editing technologies, have enabled mapping of the genetic repertoires that drive drug resistance. This Review presents examples of practical use and describes the latest technology to show the power of real-time genomic epidemiology in achieving malaria elimination.

Summary: This Review discusses the challenges in malaria elimination and how implementation of national-scale genomic surveillance programmes in combination with resistome analyses could provide a powerful solution.