Genomics reveals heterogeneous Plasmodium falciparum transmission and selection signals in Zambia

BACKGROUND

Genomic surveillance is crucial for monitoring malaria transmission and understanding parasite adaptation to interventions. Zambia lacks prior nationwide efforts in malaria genomic surveillance among African countries.

METHODS

We conducted genomic surveillance of Plasmodium falciparum parasites from the 2018 Malaria Indicator Survey in Zambia, a nationally representative household survey of children under five years of age. We whole-genome sequenced and analyzed 241 P. falciparum genomes from regions with varying levels of malaria transmission across Zambia and estimated genetic metrics that are informative about transmission intensity, genetic relatedness between parasites, and selection.

RESULTS

We provide genomic evidence of widespread within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in Zambia. Our analysis reveals country-level clustering of parasites from Zambia and neighboring regions, with distinct separation in West Africa. Within Zambia, identity by descent (IBD) relatedness analysis uncovers local spatial clustering and rare cases of long-distance sharing of closely related parasite pairs. Genomic regions with large shared IBD segments and strong positive selection signatures implicate genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Furthermore, differences in selection signatures, including drug resistance loci, are observed between eastern and western Zambian parasite populations, suggesting variable transmission intensity and ongoing drug pressure.

CONCLUSIONS

Our findings enhance our understanding of nationwide P. falciparum transmission in Zambia, establishing a baseline for analyzing parasite genetic metrics as they vary over time and space. These insights highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

Microsatellites reveal high polymorphism and high potential for use in anti-malarial efficacy studies in areas with different transmission intensities in mainland Tanzania

BACKGROUND

Tanzania is currently implementing therapeutic efficacy studies (TES) in areas of varying malaria transmission intensities as per the World Health Organization (WHO) recommendations. In TES, distinguishing reinfection from recrudescence is critical for the determination of anti-malarial efficacy. Recently, the WHO recommended genotyping polymorphic coding genes, merozoite surface proteins 1 and 2 (msp1 and msp2), and replacing the glutamate-rich protein (glurp) gene with one of the highly polymorphic microsatellites in Plasmodium falciparum to adjust the efficacy of antimalarials in TES. This study assessed the polymorphisms of six neutral microsatellite markers and their potential use in TES, which is routinely performed in Tanzania.

METHODS

Plasmodium falciparum samples were obtained from four TES sentinel sites, Kibaha (Pwani), Mkuzi (Tanga), Mlimba (Morogoro) and Ujiji (Kigoma), between April and September 2016. Parasite genomic DNA was extracted from dried blood spots on filter papers using commercial kits. Genotyping was done using six microsatellites (Poly-α, PfPK2, TA1, C3M69, C2M34 and M2490) by capillary method, and the data were analysed to determine the extent of their polymorphisms and genetic diversity at the four sites.

RESULTS

Overall, 83 (88.3%) of the 94 samples were successfully genotyped (with positive results for ≥ 50.0% of the markers), and > 50.0% of the samples (range = 47.6-59.1%) were polyclonal, with a mean multiplicity of infection (MOI) ranging from 1.68 to 1.88 among the four sites. There was high genetic diversity but limited variability among the four sites based on mean allelic richness (RS = 7.48, range = 7.27-8.03, for an adjusted minimum sample size of 18 per site) and mean expected heterozygosity (He = 0.83, range = 0.80-0.85). Cluster analysis of haplotypes using STRUCTURE, principal component analysis, and pairwise genetic differentiation (FST) did not reveal population structure or clustering of parasites according to geographic origin. Of the six markers, Poly-α was the most polymorphic, followed by C2M34, TA1 and C3M69, while M2490 was the least polymorphic.

CONCLUSION

Microsatellite genotyping revealed high polyclonality and genetic diversity but no significant population structure. Poly-α, C2M34, TA1 and C3M69 were the most polymorphic markers, and Poly-α alone or with any of the other three markers could be adopted for use in TES in Tanzania.

Genomics reveals heterogeneous Plasmodium falciparum transmission and population differentiation in Zambia and bordering countries

Genomic surveillance plays a critical role in monitoring malaria transmission and understanding how the parasite adapts in response to interventions. We conducted genomic surveillance of malaria by sequencing 241 Plasmodium falciparum genomes from regions with varying levels of malaria transmission across Zambia. We found genomic evidence of high levels of within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in the country. We identified country-level clustering of parasites from Zambia and neighboring countries, and distinct clustering of parasites from West Africa. Within Zambia, our identity by descent (IBD) relatedness analysis uncovered spatial clustering of closely related parasite pairs at the local level and rare cases of long-distance sharing. Genomic regions with large shared IBD segments and strong positive selection signatures identified genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Together, our findings enhance our understanding of P. falciparum transmission nationwide in Zambia and highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

Evidence of artemisinin partial resistance in North-western Tanzania: clinical and drug resistance markers study

Background

Artemisinin-based combination therapies (ACTs) are the recommended antimalarial drugs for the treatment of uncomplicated malaria. The recent emergence of artemisinin partial resistance (ART-R) in Rwanda, Uganda and Eritrea is of great concern. In Tanzania, a nationwide molecular malaria surveillance in 2021 showed a high prevalence of the Kelch13 (K13) 561H mutation in Plasmodium falciparum from the north-western region, close to the border with Rwanda and Uganda. This study was conducted in 2022 to evaluate the efficacy of artemether-lumefantrine (AL) and artesunate-amodiaquine (ASAQ) for the treatment of uncomplicated falciparum malaria and to confirm the presence of ART-R in Tanzania.

Methods

This single-arm study evaluated the efficacy of AL and ASAQ in eligible children aged six months to 10 years at Bukangara Dispensary in Karagwe District, Kagera Region. Clinical and parasitological responses were monitored for 28 days according to standard WHO protocol. Mutations in K13 gene and extended haplotypes with these mutations were analysed using Sanger and whole genome sequencing data, respectively.

Findings

176 children (88 in each AL and ASAQ group) were enrolled and all achieved the defined outcomes. PCR-corrected adequate clinical and parasitological response (ACPR) was 98.3% (95% CI: 90.8-100) and 100.0% (95% CI: 95.8-100) for AL and ASAQ, respectively. Parasitaemia on day 3 was observed in 11/88 (12.5%) and 17/88 (19.3%) in the AL and ASAQ groups, respectively. The half-life of parasitaemia was significantly higher (>6.5 hrs) in patients with parasitaemia on day 3 and/or mutations in K13 gene at enrolment. Most patients with parasitaemia on day 3 (8/11 = 72.7% in the AL group and 10/17 = 58.8% in the ASAQ group) had 561H mutation at enrolment. The parasites with K13 mutations were not similar to those from south-east Asia and Rwanda, but had the same core haplotype of a new 561H haplotype reported in Kagera in 2021.

Interpretation

These findings confirm the presence of ART-R in Tanzania. A context-specific strategy to respond to artemisinin partial resistance is urgently needed. Although both AL and ASAQ showed high efficacy, increased vigilance for reduced efficacy of these ACTs and detection of ART-R in other parts of the country is critical.

Emergence of Plasmodium falciparum strains with artemisinin partial resistance in East Africa and the Horn of Africa: is there a need to panic?

The emergence and spread of artemisinin partial resistance in East and Horn of Africa is alarming. However, artemisinin-based combination therapy (ACT) generally remains efficacious for the treatment of falciparum malaria. The emergence of partial artemisinin resistance does not currently meet the criteria to initiate change on treatment guidelines nor affect ACT routine procurement and distribution. It is high time for scientists and transitional researchers to be more critical and vigilant on further changes so that national programmes will be able to make informed decisions as well as remain alert and prepared for any change that may be required in the future.

Country wide surveillance reveals prevalent artemisinin partial resistance mutations with evidence for multiple origins and expansion of high level sulfadoxine-pyrimethamine resistance mutations in northwest Tanzania

Background

Emergence of artemisinin partial resistance (ART-R) in Plasmodium falciparum is a growing threat to the efficacy of artemisinin combination therapies (ACT) and the efforts for malaria elimination. The emergence of Plasmodium falciparum Kelch13 (K13) R561H in Rwanda raised concern about the impact in neighboring Tanzania. In addition, regional concern over resistance affecting sulfadoxine-pyrimethamine (SP), which is used for chemoprevention strategies, is high.

Methods

To enhance longitudinal monitoring, the Molecular Surveillance of Malaria in Tanzania (MSMT) project was launched in 2020 with the goal of assessing and mapping antimalarial resistance. Community and clinic samples were assessed for resistance polymorphisms using a molecular inversion probe platform.

Findings

Genotyping of 6,278 samples collected countrywide in 2021 revealed a focus of K13 561H mutants in northwestern Tanzania (Kagera) with prevalence of 7.7% (50/649). A small number of 561H mutants (about 1%) were found as far as 800 km away in Tabora, Manyara, and Njombe. Genomic analysis suggests some of these parasites are highly related to isolates collected in Rwanda in 2015, supporting regional spread of 561H. However, a novel haplotype was also observed, likely indicating a second origin in the region. Other validated resistance polymorphisms (622I and 675V) were also identified. A focus of high sulfadoxine-pyrimethamine drug resistance was also identified in Kagera with a prevalence of dihydrofolate reductase 164L of 15% (80/526).

Interpretation

These findings demonstrate the K13 561H mutation is entrenched in the region and that multiple origins of ART-R, similar as to what was seen in Southeast Asia, have occurred. Mutations associated with high levels of SP resistance are increasing. These results raise concerns about the long-term efficacy of artemisinin and chemoprevention antimalarials in the region.

Funding

This study was funded by the Bill and Melinda Gates Foundation and the National Institutes of Health.

Plasmodium falciparum resistant to artemisinin and diagnostics have emerged in Ethiopia

Diagnosis and treatment of Plasmodium falciparum infections are required for effective malaria control and are pre-requisites for malaria elimination efforts; hence we need to monitor emergence, evolution and spread of drug- and diagnostics-resistant parasites. We deep sequenced key drug-resistance mutations and 1,832 SNPs in the parasite genomes of 609 malaria cases collected during a diagnostic-resistance surveillance study in Ethiopia. We found that 8.0% (95% CI 7.0-9.0) of malaria cases were caused by P. falciparum carrying the candidate artemisinin partial-resistance kelch13 (K13) 622I mutation, which was less common in diagnostic-resistant parasites mediated by histidine-rich proteins 2 and 3 (pfhrp2/3) deletions than in wild-type parasites (P = 0.03). Identity-by-descent analyses showed that K13 622I parasites were significantly more related to each other than to wild type (P < 0.001), consistent with recent expansion and spread of this mutation. Pfhrp2/3-deleted parasites were also highly related, with evidence of clonal transmissions at the district level. Of concern, 8.2% of K13 622I parasites also carried the pfhrp2/3 deletions. Close monitoring of the spread of combined drug- and diagnostic-resistant parasites is needed.

Temporal and spatial analysis of Plasmodium falciparum genomics reveals patterns of parasite connectivity in a low-transmission district in Southern Province, Zambia

BACKGROUND

Understanding temporal and spatial dynamics of malaria transmission will help to inform effective interventions and strategies in regions approaching elimination. Parasite genomics are increasingly used to monitor epidemiologic trends, including assessing residual transmission across seasons and importation of malaria into these regions.

METHODS

In a low and seasonal transmission setting of southern Zambia, a total of 441 Plasmodium falciparum samples collected from 8 neighbouring health centres between 2012 and 2018 were genotyped using molecular inversion probes (MIPs n = 1793) targeting a total of 1832 neutral and geographically informative SNPs distributed across the parasite genome. After filtering for quality and missingness, 302 samples and 1410 SNPs were retained and used for downstream population genomic analyses.

RESULTS

The analyses revealed most (67%, n = 202) infections harboured one clone (monogenomic) with some variation at local level suggesting low, but heterogenous malaria transmission. Relatedness identity-by-descent (IBD) analysis revealed variable distribution of IBD segments across the genome and 6% of pairs were highly-related (IBD ≥ 0.25). Some of the highly-related parasite populations persisted across multiple seasons, suggesting that persistence of malaria in this low-transmission region is fueled by parasites "seeding" across the dry season. For recent years, clusters of clonal parasites were identified that were dissimilar to the general parasite population, suggesting parasite populations were increasingly fragmented at small spatial scales due to intensified control efforts. Clustering analysis using PCA and t-SNE showed a lack of substantial parasite population structure.

CONCLUSION

Leveraging both genomic and epidemiological data provided comprehensive picture of fluctuations in parasite populations in this pre-elimination setting of southern Zambia over 7 years.

Genomic surveillance of SARS-CoV-2 in a university community: insights into tracking variants, transmission, and spread of Gamma (P.1) variant

Background

Using a combination of data from routine surveillance, genomic sequencing, and phylogeographic analysis we tracked the spread and introduction events of SARS-CoV-2 variants focusing on a large university community.

Methods

Here, we sequenced and analyzed 677 high-quality SARS-CoV-2 genomes from positive RNA samples collected from Purdue University students, faculty, and staff who tested positive for the virus between January 2021 and May 2021, comprising an average of 32% of weekly cases across the time frame.

Results

Our analysis of circulating SARS-CoV-2 variants over time revealed periods when Variant of Concern (VOC) Alpha (B.1.1.7) and Iota (B.1.526) reached rapid dominance and documented that VOC Gamma (P.1) was increasing in frequency as campus surveillance was ending. Phylodynamic analysis of Gamma genomes from campus alongside a subsampling of &gt;20,000 previously published P.1 genomes revealed ten independent introductions of this variant into the Purdue community, predominantly from elsewhere in the United States, with introductions from within the state of Indiana and from Illinois, and possibly Washington and New York, suggesting a degree of domestic spread.

Conclusions

We conclude that a robust and sustained active and passive surveillance program coupled with genomic sequencing during a pandemic offers important insights into the dynamics of pathogen arrival and spread in a campus community and can help guide mitigation measures.

Performance evaluation of RDT, light microscopy, and PET-PCR for detecting Plasmodium falciparum malaria infections in the 2018 Zambia National Malaria Indicator Survey

Background

Zambia continues to advance on the path to elimination with significant reductions in malaria morbidity and mortality. Crucial components that have contributed to progress thus far and are necessary for achieving the national malaria elimination goals include properly identifying and treating all malaria cases through accurate diagnosis. This study sought to compare and assess the diagnostic performance of Rapid Diagnostic Tests (RDT) and Light Microscopy (LM) with photo-induced electron transfer polymerase chain reaction (PET-PCR) as the gold standard using 2018 Malaria Indicator Survey (MIS) data across Zambia to better understand diagnostic accuracy metrics and how these vary across a transmission gradient.

Methods

Cross-sectional samples collected in a nationally representative survey from 7 provinces in Zambia were tested for the presence of malaria parasites by light microscopy (LM), rapid diagnostic test (RDT) and the gold standard PET-PCR. Diagnostic performance was assessed including sensitivity, specificity, negative- and positive-predictive values across a wide malaria transmission spectrum. Diagnostic accuracy metrics were measured, and statistically significant differences were calculated between test methods for different outcome variables.

Results

From the individuals included in the MIS, the overall prevalence of Plasmodium falciparum malaria was 32.9% by RDT, 19.4% by LM, and 23.2% by PET-PCR. Herein, RDT and LM diagnostic performance was compared against gold standard PET-PCR with LM displaying a higher diagnostic accuracy than RDTs (91.3% vs. 84.6% respectively) across the transmission spectrum in Zambia. However, the performance of both diagnostics was significantly reduced in low parasitaemia samples. Consistent with previous studies, RDT diagnostic accuracy was predominantly affected by a high rate of false positives.

Conclusions

RDTs and LM both perform well across a range of transmission intensities within their respective target applications, i.e., in the community, for the former, where ease of use and speed of result is critical, and at the health facility, for the latter, where accuracy is prioritized. However, the performance of both diagnostic methods is adversely affected by low parasitaemia infections. As Zambia moves towards elimination more sensitive tools may be required to identify the last cases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03917-6.

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.

SNP barcodes provide higher resolution than microsatellite markers to measure Plasmodium vivax population genetics

BACKGROUND

Genomic surveillance of malaria parasite populations has the potential to inform control strategies and to monitor the impact of interventions. Barcodes comprising large numbers of single nucleotide polymorphism (SNP) markers are accurate and efficient genotyping tools, however may need to be tailored to specific malaria transmission settings, since 'universal' barcodes can lack resolution at the local scale. A SNP barcode was developed that captures the diversity and structure of Plasmodium vivax populations of Papua New Guinea (PNG) for research and surveillance.

METHODS

Using 20 high-quality P. vivax genome sequences from PNG, a total of 178 evenly spaced neutral SNPs were selected for development of an amplicon sequencing assay combining a series of multiplex PCRs and sequencing on the Illumina MiSeq platform. For initial testing, 20 SNPs were amplified in a small number of mono- and polyclonal P. vivax infections. The full barcode was then validated by genotyping and population genetic analyses of 94 P. vivax isolates collected between 2012 and 2014 from four distinct catchment areas on the highly endemic north coast of PNG. Diversity and population structure determined from the SNP barcode data was then benchmarked against that of ten microsatellite markers used in previous population genetics studies.

RESULTS

From a total of 28,934,460 reads generated from the MiSeq Illumina run, 87% mapped to the PvSalI reference genome with deep coverage (median = 563, range 56-7586) per locus across genotyped samples. Of 178 SNPs assayed, 146 produced high-quality genotypes (minimum coverage = 56X) in more than 85% of P. vivax isolates. No amplification bias was introduced due to either polyclonal infection or whole genome amplification (WGA) of samples before genotyping. Compared to the microsatellite panels, the SNP barcode revealed greater variability in genetic diversity between populations and geographical population structure. The SNP barcode also enabled assignment of genotypes according to their geographic origins with a significant association between genetic distance and geographic distance at the sub-provincial level.

CONCLUSIONS

High-throughput SNP barcoding can be used to map variation of malaria transmission dynamics at sub-national resolution. The low cost per sample and genotyping strategy makes the transfer of this technology to field settings highly feasible.

Higher Complexity of Infection and Genetic Diversity of Plasmodium vivax Than Plasmodium falciparum Across All Malaria Transmission Zones of Papua New Guinea

Plasmodium falciparum and Plasmodium vivax have varying transmission dynamics that are informed by molecular epidemiology. This study aimed to determine the complexity of infection and genetic diversity of P. vivax and P. falciparum throughout Papua New Guinea (PNG) to evaluate transmission dynamics across the country. In 2008-2009, a nationwide malaria indicator survey collected 8,936 samples from all 16 endemic provinces of PNG. Of these, 892 positive P. vivax samples were genotyped at PvMS16 and PvmspF3, and 758 positive P. falciparum samples were genotyped at Pfmsp2. The data were analyzed for multiplicity of infection (MOI) and genetic diversity. Overall, P. vivax had higher polyclonality (71%) and mean MOI (2.32) than P. falciparum (20%, 1.39). These measures were significantly associated with prevalence for P. falciparum but not for P. vivax. The genetic diversity of P. vivax (PvMS16: expected heterozygosity = 0.95, 0.85-0.98; PvMsp1F3: 0.78, 0.66-0.89) was higher and less variable than that of P. falciparum (Pfmsp2: 0.89, 0.65-0.97). Significant associations of MOI with allelic richness (rho = 0.69, P = 0.009) and expected heterozygosity (rho = 0.87, P < 0.001) were observed for P. falciparum. Conversely, genetic diversity was not correlated with polyclonality nor mean MOI for P. vivax. The results demonstrate higher complexity of infection and genetic diversity of P. vivax across the country. Although P. falciparum shows a strong association of these parameters with prevalence, a lack of association was observed for P. vivax and is consistent with higher potential for outcrossing of this species.

A Worldwide Map of Plasmodium falciparum K13-Propeller Polymorphisms

BACKGROUND

Recent gains in reducing the global burden of malaria are threatened by the emergence of Plasmodium falciparum resistance to artemisinins. The discovery that mutations in portions of a P. falciparum gene encoding kelch (K13)-propeller domains are the major determinant of resistance has provided opportunities for monitoring such resistance on a global scale.

METHODS

We analyzed the K13-propeller sequence polymorphism in 14,037 samples collected in 59 countries in which malaria is endemic. Most of the samples (84.5%) were obtained from patients who were treated at sentinel sites used for nationwide surveillance of antimalarial resistance. We evaluated the emergence and dissemination of mutations by haplotyping neighboring loci.

RESULTS

We identified 108 nonsynonymous K13 mutations, which showed marked geographic disparity in their frequency and distribution. In Asia, 36.5% of the K13 mutations were distributed within two areas--one in Cambodia, Vietnam, and Laos and the other in western Thailand, Myanmar, and China--with no overlap. In Africa, we observed a broad array of rare nonsynonymous mutations that were not associated with delayed parasite clearance. The gene-edited Dd2 transgenic line with the A578S mutation, which expresses the most frequently observed African allele, was found to be susceptible to artemisinin in vitro on a ring-stage survival assay.

CONCLUSIONS

No evidence of artemisinin resistance was found outside Southeast Asia and China, where resistance-associated K13 mutations were confined. The common African A578S allele was not associated with clinical or in vitro resistance to artemisinin, and many African mutations appear to be neutral. (Funded by Institut Pasteur Paris and others.).

Diagnostic performance of rapid diagnostic tests for the diagnosis of malaria at public health facilities in north-west Ethiopia

OBJECTIVE

To assess the performance of RDTs against nested polymerase chain reaction (nPCR) for the diagnosis of malaria in public health facilities in north-western Ethiopia.

METHODS

Cross-sectional study at public health facilities in North Gondar, Ethiopia, of 359 febrile patients with signs and symptoms consistent with malaria. Finger prick blood samples were collected for testing in a P. falciparum/pan-malaria RDTs and for molecular analysis. Sensitivity, specificity and predictive values were determined for the RDTs using nPCR as reference diagnostic method. Kappa value was determined to demonstrate the consistency of the results between the diagnostic tools.

RESULTS

By RDTs, 22.28% (80/359) of patients tested positive for malaria, and by nPCR, 27.02% (97/359) did. In nPCR, 1.67% (6/359) and 0.28% (1/359) samples were positive for P. ovale and P. malariae, which had almost all tested negative in the RDTs. The sensitivity, specificity, positive and negative predictive values of RDTs for the diagnosis of malaria were 62.9%, 92.7%, 76.3% and 87.1%, respectively, with 0.589 measurement agreement between RDTs and nPCR. The sensitivity and specificity of RDTs for P. falciparum identification only were 70.8% and 95.2%, and 65.2% and 93.1% for P. vivax.

CONCLUSION

Although RDTs are commonly used at health posts in resource-limited environments, their sensitivity and specificity for the detection and species identification of Plasmodium parasites were poor compared to nPCR, suggesting caution in interpreting RDTs results. Particularly, in the light of expanded efforts to eliminate malaria in the country, more sensitive diagnostic procedures will be needed.