Amplicon sequencing as a potential surveillance tool for complexity of infection and drug resistance markers in Plasmodium falciparum asymptomatic infections

Detection of low frequency artemisinin resistance mutations, C469Y, P553L and A675V, and fixed antifolate resistance mutations in asymptomatic primary school children in Kenya

BACKGROUND

To understand the emergence and spread of drug-resistant parasites in malaria-endemic areas, accurate assessment and monitoring of antimalarial drug resistance markers is critical. Recent advances in next-generation sequencing (NGS) technologies have enabled the tracking of drug-resistant malaria parasites.

METHODS

In this study, we used Targeted Amplicon Deep Sequencing (TADS) to characterise the genetic diversity of the Pfk13, Pfdhfr, Pfdhps, and Pfmdr1 genes among primary school-going children in 15 counties in Kenya (Bungoma, Busia, Homa Bay, Migori, Kakamega, Kilifi, Kirinyaga, Kisii, Kisumu, Kwale, Siaya, Tana River, Turkana, Vihiga and West Pokot). A total of 920 dried blood spot (DBS) samples collected from 121 selected primary schools within the country were used to extract genomic DNA. A nested polymerase chain reaction (PCR) was used to generate amplicons that were sequenced to determine the prevalence of known and novel polymorphisms.

RESULTS

Pfk13 mutations associated with artemisinin resistance were present as mixed genotype infections for the C469Y mutation in 23 samples (4%), the A675V mutation in 2 samples (1.7%), and the P553L mutation in 7 samples (1.2%). The A578S mutation, was also identified in mixed infections, appearing in 15.2% of the 87 samples analysed. The Pfdhfr 51I and 108 N pyrimethamine-resistance mutations were at fixation (100% frequency), and the Pfmdr1 Y184F mutation, which is linked to reduced susceptibility to several antimalarial drugs, especially those used in combination therapies for malaria treatment, was detected in 97.5% of the samples as mixed-genotype infections.

CONCLUSION

The genomic surveillance of asymptomatic school children in Kenya provides an early warning signal of at least 1 of the 3 validated artemisinin resistance mutations circulating in all regions in Western Kenya sampled except Homa Bay and Kisii Counties. These signals in asymptomatic and mixed infections would have been missed without deep sequencing.

Automated total nucleic acid extraction with magnetic beads for the detection of Plasmodium falciparum in large study cohorts

BACKGROUND

Molecular methods play an important role in clinical trials assessing anti-malarial drugs and vaccines, as well as in epidemiological studies aimed at detecting Plasmodium species, especially when dealing with large sample sizes. Molecular techniques are more sensitive and generally have a higher throughput compared to the gold standard microscopy. Further optimization can be achieved with automation of nucleic acid isolation, allowing for rapid and precise extraction. This study evaluated the isolation of total nucleic acids from Plasmodium falciparum mocked samples using an automated extraction method with a magnetic bead-based kit compared to a manual silica column-based kit. Additionally, two different RNA preservation solutions were compared.

METHODS

Plasmodium falciparum Dd2 parasites were serially diluted and spiked into whole blood. The dilutions were stored in two different RNA preservation solutions and total nucleic acids extracted with an automated magnetic bead-based kit and a manual silica column-based kit. Subsequently, a reverse transcription (RT) qPCR for Plasmodium detection targeting Plasmodium 18S rRNA and DNA in a single reaction was performed and the quantification cycle (Cq) values across the different sample groups were compared.

RESULTS

Comparable Cq values across the various sample preparations were obtained, suggesting minimal influence from RNA preservation solutions (p = 0.686) or extraction methods (p = 0.119) on RT-qPCR outcomes. Automated nucleic acids extraction allowed processing numerous samples in a shorter timeframe and showed similar efficiency in detecting Plasmodium in blood samples by RT-qPCR as manual extraction.

CONCLUSIONS

The automated method for nucleic acid isolation is a valuable tool for the detection of Plasmodium infections in large-scale studies. It is efficient, reliable, and cost-effective. Its potential applications extend to other molecular surveillance studies to support malaria control measures.

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.

ONT sequencing identifies a high prevalence of crt sensitive, triple mutant dhfr and single mutant dhps parasites within an ANC population in Nigeria

Background

Malaria in pregnancy is a major public health issue, particularly among vulnerable populations in malaria-endemic sub-Saharan African countries. To mitigate its risks, WHO recommends sulphadoxine-pyrimethamine (SP) for chemoprevention and artemisinin-based combination therapy (ACT) to treat uncomplicated Plasmodium falciparum malaria. These interventions have helped to alleviate the risk associated with malaria in pregnancy; however, in the context of the emergence of SP- and ACT-resistant P. falciparum, maintained efficacy is under threat. Molecular surveillance is a reliable tool to monitor the emergence of resistance where molecular markers are known. Thus, the objective of the study was to use a multiplexed amplicon Oxford Nanopore sequencing approach to assess the molecular markers for antimalarial resistance among pregnant women in Nigeria.

Methods

Dried blood spots (DBS) were collected from pregnant women who received IPTp-SP at the enrollment and follow-up visits. P. falciparum genomic DNA was extracted by the Chelex® method and Pf18S qPCR was used to detect parasite DNA in each sample. With nested PCR assays, fragments of Pfdhps, Pfdhfr, Pfmdr1, Pfcrt, Pfk13 and Pfama1 genes were amplified and multiplexed amplicon-based sequencing was conducted on the minION Oxford Nanopore Technology.

Result

In total, 251 pregnant women were enrolled in the study and 457 DBS samples were collected. P. falciparum genomic DNA was detected in 12% (56/457) of the samples, 31 at baseline and the remaining during the follow-up visits. Pfama1, pfk13, Pfdhps, Pfdhfr, Pfmdr1 and Pfcrt were successfully sequenced in a single run. Notably, k13 artemisinin resistance mutations were absent, the frequencies of Pfdhfr and Pfdhps SP resistance haplotypes, IRN for pyrimethamine resistance and ISGKA/IAGKA associated with sulphadoxine resistance were 82% (36/44) and 64% (27/42), respectively, and the Pfcrt CVIET resistant haplotype was at approximately 22% (7/32).

Conclusion and recommendations

Here a multiplexed amplicon-based ONT assay established that triple mutant Pfdfhr-IRN, double mutant Pfdhps-SG haplotypes and the chloroquine sensitive strain were prevalent among pregnant women in Nigeria.

Integrating targeted genetic markers to genotyping-by-sequencing for an ultimate genotyping tool

New selection methods, using trait-specific markers (marker-assisted selection (MAS)) and/or genome-wide markers (genomic selection (GS)), are becoming increasingly widespread in breeding programs. This new era requires innovative and cost-efficient solutions for genotyping. Reduction in sequencing cost has enhanced the use of high-throughput low-cost genotyping methods such as genotyping-by-sequencing (GBS) for genome-wide single-nucleotide polymorphism (SNP) profiling in large breeding populations. However, the major weakness of GBS methodologies is their inability to genotype targeted markers. Conversely, targeted methods, such as amplicon sequencing (AmpSeq), often face cost constraints, hindering genome-wide genotyping across a large cohort. Although GBS and AmpSeq data can be generated from the same sample, an efficient method to achieve this is lacking. In this study, we present the Genome-wide & Targeted Amplicon (GTA) genotyping platform, an innovative way to integrate multiplex targeted amplicons into the GBS library preparation to provide an all-in-one cost-effective genotyping solution to breeders and research communities. Custom primers were designed to target 23 and 36 high-value markers associated with key agronomical traits in soybean and barley, respectively. The resulting multiplex amplicons were compatible with the GBS library preparation enabling both GBS and targeted genotyping data to be produced efficiently and cost-effectively. To facilitate data analysis, we have introduced Fast-GBS.v3, a user-friendly bioinformatic pipeline that generates comprehensive outputs from data obtained following sequencing of GTA libraries. This high-throughput low-cost approach will greatly facilitate the application of DNA markers as it provides required markers for both MAS and GS in a single assay.

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.

A retrospective analysis of P. falciparum drug resistance markers detects an early (2016/17) high prevalence of the k13 C469Y mutation in asymptomatic infections in Northern Uganda

The emergence of drug-resistant Plasmodium falciparum parasites in sub-Saharan Africa will substantially challenge malaria control. Here, we evaluated the frequency of common drug resistance markers among adolescents from Northern Uganda with asymptomatic infections. We used an established amplicon deep sequencing strategy to screen dried blood spot samples collected from 2016 to 2017 during a reported malaria epidemic within the districts of Kitgum and Pader in Northern Uganda. We screened single-nucleotide polymorphisms within: kelch13 (Pfk13), dihydropteroate synthase (Pfdhps), multidrug resistance-1 (Pfmdr1), dihydrofolate reductase (Pfdhfr), and apical membrane antigen (Pfama1) genes. Within the study population, the median age was 15 years (14.3-15.0, 95% CI), and 54.9% (78/142) were Plasmodium positive by 18S rRNA qPCR, which were subsequently targeted for sequencing analysis. We observed a high frequency of resistance markers particularly for sulfadoxine-pyrimethamine (SP), with no wild-type-only parasites observed for Pfdhfr (N51I, C59R, and S108N) and Pfdhps (A437G and K540E) mutations. Within Pfmdr1, mixed infections were common for NF/NY (98.5%). While for artemisinin resistance, in kelch13, there was a high frequency of C469Y (34%). Using the pattern for Pfama1, we found a high level of polygenomic infections with all individuals presenting with complexity of infection greater than 2 with a median of 6.9. The high frequency of the quintuple SP drug-resistant parasites and the C469Y artemisinin resistance-associated mutation in asymptomatic individuals suggests an earlier high prevalence than previously reported from symptomatic malaria surveillance studies (in 2016/2017). Our data demonstrate the urgency for routine genomic surveillance programs throughout Africa and the value of deep sequencing.

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.

Distinct transcriptomic signatures define febrile malaria depending on initial infective states, asymptomatic or uninfected

BACKGROUND

Cumulative malaria parasite exposure in endemic regions often results in the acquisition of partial immunity and asymptomatic infections. There is limited information on how host-parasite interactions mediate the maintenance of chronic symptomless infections that sustain malaria transmission.

METHODS

Here, we determined the gene expression profiles of the parasite population and the corresponding host peripheral blood mononuclear cells (PBMCs) from 21 children (< 15 years). We compared children who were defined as uninfected, asymptomatic and those with febrile malaria.

RESULTS

Children with asymptomatic infections had a parasite transcriptional profile characterized by a bias toward trophozoite stage (~ 12 h-post invasion) parasites and low parasite levels, while early ring stage parasites were characteristic of febrile malaria. The host response of asymptomatic children was characterized by downregulated transcription of genes associated with inflammatory responses, compared with children with febrile malaria,. Interestingly, the host responses during febrile infections that followed an asymptomatic infection featured stronger inflammatory responses, whereas the febrile host responses from previously uninfected children featured increased humoral immune responses.

CONCLUSIONS

The priming effect of prior asymptomatic infection may explain the blunted acquisition of antibody responses seen to malaria antigens following natural exposure or vaccination in malaria endemic areas.

Drug resistance and vaccine target surveillance of Plasmodium falciparum using nanopore sequencing in Ghana

Malaria results in over 600,000 deaths annually, with the highest burden of deaths in young children living in sub-Saharan Africa. Molecular surveillance can provide important information for malaria control policies, including detection of antimalarial drug resistance. However, genome sequencing capacity in malaria-endemic countries is limited. We designed and implemented an end-to-end workflow to detect Plasmodium falciparum antimalarial resistance markers and diversity in the vaccine target circumsporozoite protein (csp) using nanopore sequencing in Ghana. We analysed 196 clinical samples and showed that our method is rapid, robust, accurate and straightforward to implement. Importantly, our method could be applied to dried blood spot samples, which are readily collected in endemic settings. We report that P. falciparum parasites in Ghana are mostly susceptible to chloroquine, with persistent sulfadoxine-pyrimethamine resistance and no evidence of artemisinin resistance. Multiple single nucleotide polymorphisms were identified in csp, but their significance is uncertain. Our study demonstrates the feasibility of nanopore sequencing for malaria genomic surveillance in endemic countries.

Population structure and genetic connectivity of Plasmodium falciparum in pre-elimination settings of Southern Africa

To accelerate malaria elimination in the Southern African region by 2030, it is essential to prevent cross-border malaria transmission. However, countries within the region are highly interconnected due to human migration that aids in the movement of the parasite across geographical borders. It is therefore important to better understand Plasmodium falciparum transmission dynamics in the region, and identify major parasite source and sink populations, as well as cross-border blocks of high parasite connectivity. We performed a meta-analysis using collated parasite allelic data generated by microsatellite genotyping of malaria parasites from Namibia, Eswatini, South Africa, and Mozambique (N = 5,314). The overall number of unique alleles was significantly higher (P ≤ 0.01) in Namibia (mean A = 17.3 ± 1.46) compared to South Africa (mean A = 12.2 ± 1.22) and Eswatini (mean A = 13.3 ± 1.27, P ≤ 0.05), whilst the level of heterozygosity was not significantly different between countries. The proportion of polyclonal infections was highest for Namibia (77%), and lowest for Mozambique (64%). A was significant population structure was detected between parasites from the four countries, and patterns of gene flow showed that Mozambique was the major source area and Eswatini the major sink area of parasites between the countries. This study showed strong signals of parasite population structure and genetic connectivity between malaria parasite populations across national borders. This calls for strengthening the harmonization of malaria control and elimination efforts between countries in the southern African region. This data also proves its potential utility as an additional surveillance tool for malaria surveillance on both a national and regional level for the identification of imported cases and/or outbreaks, as well as monitoring for the potential spread of anti-malarial drug resistance as countries work towards malaria elimination.