Fitness Costs of pfhrp2 and pfhrp3 Deletions Underlying Diagnostic Evasion in Malaria Parasites

Evolution of pfhrp2 and pfhrp3 deletions in Equatorial Guinea between the pre- and post-RDT introduction

BACKGROUND

Pfhrp2 and pfhrp3 deletions are threatening Plasmodium falciparum malaria diagnosis by rapid diagnostic tests (RDT) due to false negatives. This study assesses the changes in the frequencies of pfhrp2 and pfhrp3 deletions (pfhrp2Del and pfhrp3Del, respectively) and the genes in their flaking regions, before and after RDT introduction in Equatorial Guinea.

METHODS

A total of 566 P. falciparum samples were genotyped to assess the presence of pfhrp2 and pfhrp3 deletions and their flanking genes. The specimens were collected 18 years apart from two provinces of Equatorial Guinea, North Bioko (Insular Region) and Litoral Province (Continental Region). Orthologs of pfhrp2 and pfhrp3 genes from other closely related species were used to compare sequencing data to assess pfhrp2 and pfhrp3 evolution. Additionally, population structure was studied using seven neutral microsatellites.

RESULTS

This study found that pfhrp2Del and pfhrp3Del were present before the introduction of RDT; however, they increased in frequency after their use, reaching more than 15%. Haplotype networks suggested that pfhrp2Del and pfhrp3Del emerged multiple times. Exon 2 of pfhrp2 and pfhrp3 genes had high variability, but there were no significant changes in amino acid sequences.

CONCLUSIONS

Baseline sampling before deploying interventions provides a valuable context to interpret changes in genetic markers linked to their efficacy, such as the dynamic of deletions affecting RDT efficacy.

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.

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.

Assessment of the Performance of Lactate Dehydrogenase-Based Rapid Diagnostic Test for Malaria in Djibouti in 2022-2023

Until 2020, Djiboutian health authorities relied on histidine-rich protein-2 (HRP2)-based rapid diagnostic tests (RDTs) to establish the diagnosis of Plasmodium falciparum. The rapid spread of P. falciparum histidine-rich protein-2 and -3 (pfhrp2/3) gene-deleted parasite strains in Djibouti has led the authorities to switch from HRP2-based RDTs to lactate dehydrogenase (LDH)-based RDTs targeting the plasmodial lactate dehydrogenase (pLDH) specific for P. falciparum and P. vivax (RapiGEN BIOCREDIT Malaria Ag Pf/Pv pLDH/pLDH) in 2021. This study was conducted with the primary objective of evaluating the diagnostic performance of this alternative RDT. Operational constraints related, in particular, to the implementation of this RDT during the COVID-19 pandemic were also considered. The performance of BIOCREDIT Malaria Ag Pf/Pv (pLDH/pLDH) RDT was also compared to our previously published data on the performance of two HRP2-based RDTs deployed in Djibouti in 2018-2020. The diagnosis of 350 febrile patients with suspected malaria in Djibouti city was established using two batches of RapiGEN BIOCREDIT Malaria Ag Pf/Pv (pLDH/pLDH) RDT over a two-year period (2022 and 2023) and confirmed by real-time quantitative polymerase chain reaction. The sensitivity and specificity for the detection of P. falciparum were 88.2% and 100%, respectively. For P. vivax, the sensitivity was 86.7% and the specificity was 100%. Re-training and closer supervision of the technicians between 2022 and 2023 have led to an increased sensitivity to detect P. falciparum (69.8% in 2022 versus 88.2% in 2023; p < 0.01). The receiver operating characteristic curve analysis highlighted a better performance in the diagnosis of P. falciparum with pLDH-based RDTs compared with previous HRP2-based RDTs. In Djibouti, where pfhrp2-deleted strains are rapidly gaining ground, LDH-based RDTs seem to be more suitable for diagnosing P. falciparum than HRP2-based RDTs. Awareness-raising and training for technical staff have also been beneficial.

Global risk of selection and spread of Plasmodium falciparum histidine-rich protein 2 and 3 gene deletions

In the thirteen years since the first report of pfhrp2-deleted parasites in 2010, the World Health Organization (WHO) has found that 40 of 47 countries surveyed worldwide have reported pfhrp2/3 gene deletions. Due to a high prevalence of pfhrp2/3 deletions causing false-negative HRP2 RDTs, in the last five years, Eritrea, Djibouti and Ethiopia have switched or started switching to using alternative RDTs, that target pan-specific-pLDH or P. falciparum specific-pLDH alone of in combination with HRP2. However, manufacturing of alternative RDTs has not been brought to scale and there are no WHO prequalified combination tests that use Pf-pLDH instead of HRP2 for P. falciparum detection. For these reasons, the continued spread of pfhrp2/3 deletions represents a growing public health crisis that threatens efforts to control and eliminate P. falciparum malaria. National malaria control programmes, their implementing partners and test developers desperately seek pfhrp2/3 deletion data that can inform their immediate and future resource allocation. In response, we use a mathematical modelling approach to evaluate the global risk posed by pfhrp2/3 deletions and explore scenarios for how deletions will continue to spread in Africa. We incorporate current best estimates of the prevalence of pfhrp2/3 deletions and conduct a literature review to estimate model parameters known to impact the selection of pfhrp2/3 deletions for each malaria endemic country. We identify 20 countries worldwide to prioritise for surveillance and future deployment of alternative RDT, based on quickly selecting for pfhrp2/3 deletions once established. In scenarios designed to explore the continued spread of deletions in Africa, we identify 10 high threat countries that are most at risk of deletions both spreading to and subsequently being rapidly selected for. If HRP2-based RDTs continue to be relied on for malaria case management, we predict that the major route for pfhrp2 deletions to spread is south out from the current hotspot in the Horn of Africa, moving through East Africa over the next 20 years. We explore the variation in modelled timelines through an extensive parameter sensitivity analysis and despite wide uncertainties, we identify three countries that have not yet switched RDTs (Senegal, Zambia and Kenya) that are robustly identified as high risk for pfhrp2/3 deletions. These results provide a refined and updated prediction model for the emergence of pfhrp2/3 deletions in an effort to help guide pfhrp2/3 policy and prioritise future surveillance efforts and innovation.

Low Prevalence of Plasmodium falciparum Histidine-Rich Protein 2 and 3 Gene Deletions—A Multiregional Study in Central and West Africa

Plasmodium falciparum parasites carrying deletions of histidine-rich protein 2 and 3 genes, pfhrp2 and pfhrp3, respectively, are likely to escape detection via HRP2-based rapid diagnostic tests (RDTs) and, consequently, treatment, posing a major risk to both the health of the infected individual and malaria control efforts. This study assessed the frequency of pfhrp2- and pfhrp3-deleted strains at four different study sites in Central Africa (number of samples analyzed: Gabon N = 534 and the Republic of Congo N = 917) and West Africa (number of samples analyzed: Nigeria N = 466 and Benin N = 120) using a highly sensitive multiplex qPCR. We found low prevalences for pfhrp2 (1%, 0%, 0.03% and 0) and pfhrp3 single deletions (0%, 0%, 0.03% and 0%) at all study sites (Gabon, the Republic of Congo, Nigeria and Benin, respectively). Double-deleted P. falciparum were only found in Nigeria in 1.6% of all internally controlled samples. The results of this pilot investigation do not point towards a high risk for false-negative RDT results due to pfhrp2/pfhrp3 deletions in Central and West African regions. However, as this scenario can change rapidly, continuous monitoring is essential to ensure that RDTs remain a suitable tool for the malaria diagnostic strategy.

Portable and cost-effective genetic detection and characterization of Plasmodium falciparum hrp2 using the MinION sequencer

The prevalence of Plasmodium falciparum hrp2 (pfhrp2)-deleted parasites threatens the efficacy of the most used and sensitive malaria rapid diagnostic tests and highlights the need for continued surveillance for this gene deletion. While PCR methods are adequate for determining pfhrp2 presence or absence, they offer a limited view of its genetic diversity. Here, we present a portable sequencing method using the MinION. Pfhrp2 amplicons were generated from individual samples, barcoded, and pooled for sequencing. To overcome potential crosstalk between barcodes, we implemented a coverage-based threshold for pfhrp2 deletion confirmation. Amino acid repeat types were then counted and visualized with custom Python scripts following de novo assembly. We evaluated this assay using well-characterized reference strains and 152 field isolates with and without pfhrp2 deletions, of which 38 were also sequenced on the PacBio platform to provide a standard for comparison. Of 152 field samples, 93 surpassed the positivity threshold, and of those samples, 62/93 had a dominant pfhrp2 repeat type. PacBio-sequenced samples with a dominant repeat-type profile from the MinION sequencing data matched the PacBio profile. This field-deployable assay can be used alone for surveilling pfhrp2 diversity or as a sequencing-based addition to the World Health Organization's existing deletion surveillance protocol.