Molecular markers for malaria genetic epidemiology: progress and pitfalls

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

BACKGROUND

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.

METHODS

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.

RESULTS

Here we show that the mean multiplicity of infection (MOI) is 1.3, with 11% polyclonal infections, indicating low transmission intensity. Temporal trends of these genetic metrics do 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 are genetically linked, continuous turnover in genetic clusters and changes in drug-resistance haplotypes are observed. Principal component analysis reveals some STP samples are genetically similar to those from Central and West Africa, indicating possible importation.

CONCLUSIONS

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.

Var genes, strain hyperdiversity, and malaria transmission dynamics

The microbiological paradigm for surveillance of diverse pathogens requires knowledge of the variation of the major surface antigen under the most intense immune selection as immune responses to these antigens drive transmission dynamics. This creates a pathway for population genetics/genomics to be combined with mathematical modelling to describe transmission dynamics to inform public health policy. Here we consider how we can bring population genetics and population dynamics together for a highly recombining pathogen like Plasmodium falciparum. We do this through the lens of what has been recently learnt about the population genetics of the var multigene family encoding the major surface antigen of the blood stages of Plasmodium falciparum, known as PfEMP1.

Detection of novel Plasmodium falciparum haplotypes under treatment pressure in paediatric severe malaria

Background. In Africa, the clearance time for Plasmodium falciparum severe malaria varies significantly, likely due to the complexity of P. falciparum infections and the sequestration phenomenon exhibited by this parasite. This study aims to evaluate different methods to study the intra-host dynamics of polyclonal infections during parasite clearance under antimalarial treatment. Additionally, it seeks to determine the association between parasite clearance rate following artesunate or quinine treatment and the genetic complexity of P. falciparum in Beninese children with severe malaria.Methods. Sixty-five P. falciparum severe malaria individuals diagnosed by microscopy and treated with artesunate or quinine were sampled every 8 h for 24 h. Using whole-genome sequencing (WGS) data, we estimated the multiplicity of infection (MOI) with three algorithms (Fws, THE REAL McCOIL and RoH). We then characterized the P. falciparum genetic complexity in WGS-identified polyclonal infections using amplicon sequencing (AmpSeq) on DNA extracted from plasma and the red blood cell pellet.Results. AmpSeq demonstrated greater sensitivity in detecting multiple genomes within isolates compared to WGS methods. The MOI from AmpSeq was significantly higher in red blood cell pellets compared to plasma (2.4 vs. 1.8 distinct microhaplotypes per isolate). However, at parasitaemia over 1,000 parasites per microlitre, the same MOI was detected in both plasma and pellet samples in 85.4% of the isolates. We observed a high variability in parasite clearance rate among participants, but it was not associated with parasite MOI at diagnosis. Interestingly, in 60.9% of participants, previously undetected microhaplotypes appeared in circulation 16 h after treatment initiation.Conclusion. These findings demonstrate that combining different haplotyping techniques effectively determines parasite genetic complexity. Additionally, plasma can be effectively used for parasite genotyping at sufficient parasitaemia levels. The parasite clearance rate of severe malaria is independent of parasite MOI. However, genotyping a single blood sample upon hospital admission does not capture the full spectrum of parasite genotypes present in the infection.

Detecting imported malaria infections in endemic settings using molecular surveillance: current state and challenges

The Global Technical Strategy for Malaria 2016-2030 targets eliminating malaria from at least 35 countries and reducing case incidence by 90% globally. The importation of parasites due to human mobilization poses a significant obstacle to achieve malaria elimination as it can undermine the effectiveness of local interventions. Gaining a comprehensive understanding of parasite importation is essential to support control efforts and advance progress toward elimination. Parasite genetic data is widely used to investigate the spatial and temporal dynamics of imported infections. In this context, this systematic review aimed to aggregate evidence on the application of parasite genetic data for mapping imported malaria and the analytical methods used to analyze it. We discuss the advantages and limitations of the genetic approaches employed and propose a suitable type of genetic data along with an analytical framework to discriminate imported malaria infections from local infections. The findings offer potential actionable insights for national control programs, enabling them select the most effective methods for detecting imported cases. This also may aid in the evaluation and refinement of elimination programs by identifying high-risk areas and enabling the targeted allocation of resources to these regions.

Performance of Molecular Inversion Probe DR23K and Paragon MAD4HatTeR Amplicon Sequencing Panels for Detection of Plasmodium falciparum Mutations Associated with Antimalarial Drug Resistance

Background

Molecular surveillance of drug-resistant Plasmodium falciparum is crucial for malaria control in endemic regions. Two targeted-resequencing tools, the Molecular Inversion Probe (MIP) drug resistance panel DR23K and the Multiplexed Amplicons for Drugs, Diagnostics, Diversity, and Differentiation using High-Throughput Targeted Resequencing (MAD4HatTeR) panel, are widely used to detect resistance genotypes. However, comparisons of their performance for genotyping drug resistance polymorphisms in malaria parasites and their comparative utility for other use cases is lacking.

Methods

To compare the performance of DR23K and MAD4HatTeR in terms of sequencing depth, sensitivity to minor alleles, and precision, each platform was used to evaluate SNP alleles and microhaplotypes in double- and triple-strain mixtures of well-characterized laboratory parasites at densities of 10, 100, 1,000, and 10,000 parasites/μL. In addition, 67 Ugandan field samples collected in 2022 were genotyped using each platform to assess performance and concordance.

Results

Across the four parasite densities, MAD4HatTeR exhibited superior sequencing depth (mean reads per locus: 144, 992, 1,153, and 1,300) compared to DR23K (mean unique molecular identifiers [UMIs] per locus: 1, 4, 49, and 364). For SNP detection, MAD4HatTeR achieved 100% sensitivity at 2% within-sample allele frequency (WSAF) at 1,000 and 10,000 parasites/μL, whereas DR23K achieved 100% sensitivity at 40% and 5% WSAF at these densities, respectively. Microhaplotype sensitivity was lower for both assays; MAD4HatTeR reached 69% sensitivity at 10 parasites/μL when WSAF was ≥ 10%, increasing to 100% sensitivity at 2% WSAF and 100 parasites/μL. DR23K had < 50% sensitivity at 10 and 100 parasites/μL. In field samples, which commonly contain polyclonal infections, high concordance was observed between the two methods for all SNPs (94%, 1,848/1,969) and polymorphic SNPs (88%, 898/1,019). All discrepancies were attributed to varied detection of minority alleles in mixed genotype infections.

Conclusions

MAD4HatTeR demonstrated higher sensitivity than DR23K, particularly at low parasite densities. Both assays showed strong concordance for genotyping key resistance mutations in field samples, supporting their reliability. These findings suggest MAD4HatTeR as the preferred assay for low-density parasite studies and microhaplotype analysis, while DR23K may be appropriate for specific applications with high-parasite density samples, where detection of minority alleles is not prioritized, or when more comprehensive genome coverage is required.

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.

Imported malaria cases by Plasmodium falciparum and Plasmodium vivax in Mexican territory: Potential impact of the migration crisis

BACKGROUND

As the migratory flow to the USA has intensified in recent months, health problems associated have been identified. The aim of this work was the identification of malaria cases imported into Mexican territory.

METHODS

Operational definitions of suspected and confirmed cases were used for investigation of malaria cases. Detection of parasitic entities by thick blood smear and molecular biology served as a confirmatory test. With the characteristics of the cases, a heat map was made to determine common clinical pictures. Finally, epidemiological analysis of cases was performed for the construction of timelines of imported malaria and the tracing of migratory routes.

RESULTS

Twelve migrants from four countries were treated for presenting clinical symptoms with suspected dengue or malaria. Malaria was confirmed and two Plasmodium species were identified. From the epidemiological dates of arrival in Mexico, onset of symptoms and migratory routes, we speculate that ten cases acquired P. vivax during their crossing through Honduras, El Salvador or Guatemala. For the Guinea cases, we conclude that there was African importation of P. falciparum.

CONCLUSION

The epidemiological panorama of malaria cases imported into Mexico show the need to join efforts to ensure universal access to health services, with the objective of timely detection of imported cases.

Review of MrsFreqPhase methods: methods designed to estimate statistically malaria parasite multiplicity of infection, relatedness, frequency and phase

Malaria parasites are haploid within humans, but infections often contain genetically distinct groups of clonal parasites. When the per-infection number of genetically distinct clones (i.e., the multiplicity of infection, MOI) exceeds one, and per-infection genetic data are generated in bulk, important information are obfuscated. For example, the MOI, the phases of the haploid genotypes of genetically distinct clones (i.e., how the alleles concatenate into sequences), and their frequencies. This complicates many downstream analyses, including relatedness estimation. MOIs, parasite sequences, their frequencies, and degrees of relatedness are used ubiquitously in malaria studies: for example, to monitor anti-malarial drug resistance and to track changes in transmission. In this article, MrsFreqPhase methods designed to estimate statistically malaria parasite MOI, relatedness, frequency and phase are reviewed. An overview, a historical account of the literature, and a statistical description of contemporary software is provided for each method class. The article ends with a look towards future method development, needed to make best use of new data types generated by cutting-edge malaria studies reliant on MrsFreqPhase methods.

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.

Lineage-informative microhaplotypes for recurrence classification and spatio-temporal surveillance of Plasmodium vivax malaria parasites

Challenges in classifying recurrent Plasmodium vivax infections constrain surveillance of antimalarial efficacy and transmission. Recurrent infections may arise from activation of dormant liver stages (relapse), blood-stage treatment failure (recrudescence) or reinfection. Molecular inference of familial relatedness (identity-by-descent or IBD) can help resolve the probable origin of recurrences. As whole genome sequencing of P. vivax remains challenging, targeted genotyping methods are needed for scalability. We describe a P. vivax marker discovery framework to identify and select panels of microhaplotypes (multi-allelic markers within small, amplifiable segments of the genome) that can accurately capture IBD. We evaluate panels of 50-250 microhaplotypes discovered in a global set of 615 P. vivax genomes. A candidate global 100-microhaplotype panel exhibits high marker diversity in the Asia-Pacific, Latin America and horn of Africa (median HE = 0.70-0.81) and identifies 89% of the polyclonal infections detected with genome-wide datasets. Data simulations reveal lower error in estimating pairwise IBD using microhaplotypes relative to traditional biallelic SNP barcodes. The candidate global panel also exhibits high accuracy in predicting geographic origin and captures local infection outbreak and bottlenecking events. Our framework is open-source enabling customised microhaplotype discovery and selection, with potential for porting to other species or data resources.

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.