key: cord-0970356-0nmfk6rs
authors: Yobi, Doudou M.; Kayiba, Nadine K.; Mvumbi, Dieudonné M.; Boreux, Raphael; Kabututu, Pius Z.; Akilimali, Pierre Z.; Situakibanza, Hippolyte N. T.; De Mol, Patrick; Speybroeck, Niko; Mvumbi, Georges L.; Hayette, Marie-Pierre
title: Biennial surveillance of Plasmodium falciparum anti-malarial drug resistance markers in Democratic Republic of Congo, 2017 and 2019
date: 2022-02-10
journal: BMC Infect Dis
DOI: 10.1186/s12879-022-07112-z
sha: 40ce64e42f8547554c1ea682daf24830af311043
doc_id: 970356
cord_uid: 0nmfk6rs

BACKGROUND: Because of the loss of chloroquine (CQ) effectiveness, the Democratic Republic of Congo (DRC)’s malaria treatment policy replaced CQ by sulfadoxine–pyrimethamine (SP) as first-line treatment of uncomplicated malaria in 2003, which in turn was replaced by artemisinin-based combination therapies (ACT) in 2005. The World Health Organization (WHO) recommends monitoring of anti-malarial drug resistance every 2 years. The study aimed to provide baseline data for biennial molecular surveillance of anti-malarial drug resistance by comparing data from a study conducted in 2019 to previously published data from a similar study conducted in 2017 in the DRC. METHODS: From July to November 2019, a cross-sectional study was conducted in ten sites which were previously selected for a similar study conducted in 2017 across the DRC. P. falciparum malaria was diagnosed by a rapid diagnostic test (RDT) or by microscopy and dried blood samples (DBS) were taken from patients who had a positive test. Segments of interest in pfcrt and pfk13 genes were amplified by conventional PCR before sequencing. RESULTS: Out of 1087 enrolled patients, 906 (83.3%) were PCR-confirmed for P. falciparum. Like in the 2017-study, none of the mutations known to be associated with Artemisinine (ART) resistance in Southeast Asia was detected. However, non-synonymous (NS) mutations with unknown functions were observed among which, A578S was detected in both 2017 and 2019-studies. The overall prevalence of pfcrt-K76T mutation that confers CQ-resistance was 22.7% in 2019-study compared to 28.5% in 2017-study (p-value = 0.069), but there was high variability between sites in the two studies. Like in 2017-study, the pfcrt 72–76 SVMNT haplotype associated with resistance to amodiaquine was not detected. CONCLUSION: The study reported, within 2 years, the non-presence of molecular markers currently known to be associated with resistance to ART and to AQ in P. falciparum isolated in the DRC. However, the presence of polymorphisms with as-yet unknown functions was observed, requiring further characterization. Moreover, an overall decrease in the prevalence of CQ-resistance marker was observed in the DRC, but this prevalence remained highly variable from region to region. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s12879-022-07112-z.

Background The spread of resistance to anti-malarial drugs poses a serious public health risk and a real challenge for malaria control in endemic countries. As result of high treatment failure rate (> 10%) with chloroquine (CQ) and sulfadoxine-pyrimethamine (SP) and following the World Health Organization (WHO) recommendation, the Democratic republic of Congo (DRC) malaria national policy implemented artemisinin-based Combination Therapies (ACT) as first-line treatment of uncomplicated Plasmodium falciparum malaria in 2005 [1, 2] . WHO recommends biennial surveillance of anti-malarial drug efficacy to ensure early detection of emergence and spread of resistance. Therapeutic efficacy studies (TESs) are the gold standard for the monitoring of anti-malarial efficacy. However, today, in addition to TETs, the assessment of molecular marker in P. falciparum genes associated with anti-malarial drug resistance provides useful information about the emergence and spread of resistant parasites [3] .

Point mutations in the chloroquine resistance transporter gene (pfcrt), notably mutations in amino acids position 72 to 76 have been found to confer resistance to both CQ and amodiaquine (AQ) [4] . The amino acid substitution from lysine (K) to threonine (T) at position 76 of the PFCRT protein (K76T) has been shown to confer CQ resistance [5] , while the 72-76 PFCRT haplotype, namely SVMNT (Serine-Valine-Methionine-Asparagine-Threonine) has been linked to AQ resistance [6] . K76T mutation, which is considered as the most reliable molecular marker of CQ resistance, is found in all CQ-resistant isolates of P. falciparum [5] . The decline in this CQ-resistance marker following the official CQ removal from national treatment guidelines varies significantly between countries. In several countries, it was reported that after the discontinuance of CQ use, K76T-mutant parasites have been replaced by wild ones, resulting in the decrease in prevalence of K76T marker [7] [8] [9] [10] . In other countries, studies reported the persistence of resistant strains several years after stopping the use of CQ as first-line treatment of uncomplicated malaria [11, 12] . In the DRC, the prevalence of K76T marker was 55.4% in a study conducted with samples obtained from 2007 Demographic and Health Survey (2007-DHS) [13] and 63.9% in a study conducted in 2014 [14] . Moreover, a large variability in this prevalence of K76T was observed across the DRC in 2017, with presence of very high prevalence (89.5%) in one province and very low (1.5%) in another [15] .

Single-nucleotide polymorphisms (SNPs) occurring within P. falciparum Kelch 13 gene (pfk13), precisely in the propeller domain region, have been described as molecular markers of ART resistance [16, 17] . The list of these SNPs categorized as validated mutations and candidate mutations of ART resistance is updated by the WHO [18] . Originally, these mutations have emerged and spread in Southeast Asia, but nowadays, some of them are increasingly detected in sub-Saharan African countries providing evidence for de novo emergence of resistance to ART in Sub-Saharan African countries [19] [20] [21] [22] .

To help improve long-term monitoring of anti-malarial drug resistance, the present study aimed to provide baseline data for biennial molecular surveillance of P. falciparum anti-malarial drug resistance by comparing data from the study conducted in 2019 to previously published data from a similar study conducted in 2017 in the DRC [15, 23] .

From July to November 2019, a cross-sectional study was conducted in ten sites (sentinel sites of National Malaria Control Programme) located in ten provinces: Kingasani in Kinshasa, Kabondo in Kisangani/Tshopo, Lubumbashi in Haut-Katanga, Bolenge in Equateur, Karawa in Nord-Ubangi, Vanga in Kwilu, Kalima in Maniema, Kamina in Lomami, Fungurume in Lwalaba and Katana in Sud-Kivu ( Fig. 1) . P. falciparum malaria was diagnosed using a rapid diagnostic test (RDT) or the microscopy in all patients with fever or history of fever in the last 24 h, who presented at one of the selected health structures of each study site. The RDT used was the one available on the site during sample collection (SD Bioline Malaria Ag P.f/ Standard Diagnostics/South Korea or CareStart Malaria Pf/Access Bio/South Korea) while the light microscopy of Giemsa-stained thick blood test was used in some referral health structures where it was feasible. To participate in the study, informed consent was obtained from each adult patient (≥ 18 years old) or from a parent or legal guardian for children and adolescents. Blood samples were collected on filter paper as dried blood spots (DBS) and then transported at the University of Kinshasa for storage at -20 °C until use.

Extraction of genomic DNA from DBS was carried out using the QIAamp DNA Mini Kit (QIAGEN GmbH, Keywords: Biennial, Surveillance, Plasmodium falciparum, Anti-malarial, Resistance, Markers, Democratic Republic of Congo Germany) according to the manufacturer's instruction (https:// www. qiagen. com). Briefly, the DNA purification was performed following the manufacturer's manual procedure which included four steps: cell lysis, DNA binding to the QIAamp membrane, washing contaminants and DNA elution. The extracted DNA was kept at − 20 °C until time of use. To detect P. falciparum, a real-time Pf-PCR assay was performed according to a previously described procedure [24] , which was previously modified [15] . Samples with PCR-confirmed presence of P. falciparum were used to amplify segments of interest on pfk13 and pfcrt genes. The procedures previously described were used to amplify the interest segment on pfcrt gene [25] and that on the propeller region of pfk13 gene [23] .

The analysis of P. falciparum genes polymorphisms was done by DNA bi-directional sequencing using Sanger method as previously described [15, 23] . Nucleotides sequences were analyzed using Vector NTI (Thermo Fisher Scientific, US) and aligned to the reference annotated PF3D7 (PF3D7_0709000 for PFCRT sequences and PF3D7_1343700 for PFK13 sequences) to assess SNPs.

Data were encoded and verified in a Microsoft Excel 2010 database and imported in SPSS Statistics 20.0 for analysis. Only isolates successfully sequenced have been considered to calculate the prevalence of each SNP. Frequencies were statistically compared using Chisquared tests. Secondary analyses were conducted using logistic regression model to determine if the marginal associations from the univariate model matched those of the multivariate model. The p-values less than 0.05 were considered significant.

In total, 1087 patients were enrolled in 2019-study in comparison to 1070 in 2017-study. Their age ranged from 0 to 77 years old. Table 1 presents the distribution of patients according to different characteristics in comparison to 2017-study patients. The distribution per age range showed a low difference (p-value = 0.01) between 2017-study and 2019-study while no difference was observed when considered per sex. Among enrolled patients, 906 (83.3%) were confirmed positive for P. falciparum by real-time PCR in 2019-study versus 805 (75.3%) in 2017-study (p < 0.001).

As shown in Table 2 , the prevalence of CQ-resistance marker was not different between age groups and sex categories (p > 0.05). Globally, this prevalence decreased 

In 2019-study, among 822 isolates successfully sequenced on pfk13-propeller gene, following nonsynonymous mutations have been detected: M472L, D516A, V568M, S576T, A578S, D584E and Y588C (Additional file 1: Table S1 ), A578S and D584E were each found in two isolates. A578S mutation was the only one detected in both 2017 and 2019-studies in Fungurume and in Kinshasa while all other mutations were detected once (p < 0.001). On pfcrt gene, when considered in the same site, the prevalence of K76T mutation did not shown significant difference between 2019-study and 2017-study (p > 0.05), except in Kinshasa (Table 3) , while when comparing site by site, this prevalence remained highly variable in the two studies ranging from 1.5 to 89.5% in 2017-study and from 1.4 to 84.2% in 2019-study respectively in Fungurume and in Katana (p < 0.001) (Fig. 2) . The 72-76 PfCRT CVIET haplotype was the most prevalent among the K76Tmutant parasites (93.6% in 2019-study versus 95.9% in 2017-study), while the SVMNT haplotype was detected neither by the 2019-study nor the 2017-study.

Unlike univariate model, the multivariate logistic regression model showed no statistical difference between the prevalence of CQ-resistance from 2017study and that from 2019-study (p = 0.069; OR = 0.76; IC 95%: 0.56-1.02). However, the between sites variability of this prevalence has been confirmed by the logistic regression model, as shown in Table 4 .

Like in 2017-study [23] , the 2019-study did not detect any polymorphisms that are currently known to be associated with resistance to ART. However, other coding substitutions observed, such as A578S, should be considered. A578S mutation, commonly reported in African countries including the DRC [23, 26, 27] , was the only one which was detected in both 2019-study and 2017-study. Although, the amino acid substitution in A578S mutation is likely to alter the function of PFK13 protein [28] , previous studies have not shown the link with ART resistance [17, 27] . The function of the pfk13 gene is largely unknown and the phenotypic expression of this gene could involve other determinants which could vary from one geographical region to another due to the genetic diversity of P falciparum. Other unusual polymorphisms have been observed by both the 2019-study and 2017study; some of them (M472I, V534A and A569T) were previously reported in African countries [26] . Since the polymorphisms associated with ART resistance are numerous and occur on a background of genetic diversity, the mutations that are increasingly observed in Africa merit further characterization [29, 30] . Globally, the prevalence of the K76T mutation known to be associated with CQ-resistance continues to decline in the DRC. Indeed, compared to previously reported prevalence which was more than 50% [13, 14, 25] , the present study has shown that the prevalence of CQresistance marker has decreased to less than 30% (28.5% in 2017-study and 22.7% in 2019-study) in the DRC. However, this prevalence remained very variable from one region to another, with the highest prevalence in Katana (89.5% in 2017-study and 84.2% in 2019-study) and the lowest in Fungurume (1.5% in 2017-study and 1.4% in 2019-study). A possible illicit use of CQ that contributes in maintaining CQ selective pressure on the local P. falciparum population could explain the persistence of high CQ resistance rate in some regions. In addition, the selection for CQ-resistance is largely dependent on the genetic structure of parasite populations, which have been shown to vary greatly between regions of the DRC [30] ; this may partly explain the observed regional variability in prevalence of CQ-resistance rate. Other potential explanations for the heterogeneity in prevalence of CQ-resistance marker in different regions of the DRC have been addressed in the previous 2017-study [15] . The return of CQ-susceptible malaria may present a novel opportunity for reintroducing this molecule to prevent malaria, especially in vulnerable populations. Indeed, if the proportion of CQ-resistant parasites declines in African endemic countries to an undetectable level, the reintroduction of CQ to be used alone or in combination with other anti-malarial drugs for prophylaxis or treatment may possibly be considered. The absence of drug pressure is supposed to be the key driver in the reemergence of CQ-sensitivity parasites in the field setting [31] . Unfortunately, with the advent of the COVID-19 pandemic, CQ has been largely used for its supposed antiviral properties in affected countries including the DRC [32] . In this context, CQ has been made widely available in private drugstores with the presence of falsified version containing low amount of CQ, which was discovered in the DRC [33] . This is likely to maintain CQ pressure on P. falciparum populations, thus compromising the hope towards the return of CQ-susceptible malaria in endemic countries particularly in the DRC.

The PFCRT 72-76 haplotype, SVMNT has been correlated with high-level resistance to the AQ metabolite desethylamodiaquine (DEAQ) in in vitro tests [34] . Like in 2017-study, the SVMNT haplotype associated with AQ-resistance was not detected in 2019-study. AQ is one of the partner drugs of ACTs used in many endemic countries like the DRC for the treatment of uncomplicated P. falciparum malaria. Generally, AQ was found to be effective in Africa, even in the face of increasing CQ-resistance [35] that would constitute a risk factor of emergence of AQ-resistance in Africa. However, high prevalence of SVMNT has been reported in P. falciparum from South America and Asia [36] and the presence of this haplotype was also reported in Tanzania [37] and in Angola [38] . AQ is structurally related to CQ with both belonging to the group of 4-aminoquinolines. There is a possibility that AQ could continue to contribute to The present study contributed to ongoing malaria surveillance providing baseline data for the biennial surveillance of P. falciparum molecular markers associated with resistance to anti-malarial drugs in the DRC. Regular studies with size of samples taking into account the regional variability and a broader exploration of the genome of P. falciparum could help to provide more information about the landscape of P. falciparum drug resistance in the DRC.

The study did not report, within 2 years, the presence of molecular markers currently known to be associated with resistance to ART and to AQ in P. falciparum isolated in the DRC. Moreover, the overall prevalence of CQ-resistance marker continues to decrease in the DRC, but there was observed a high variability from region to region with persistence of high prevalence of parasites carrying CQ-resistant haplotype. Further characterization of unusual detected polymorphisms is required and regular surveillance should continue in order to ensure early detection and containment of the emerging anti-malarial drug resistance in the DRC. 

Assessing the efficacy of chloroquine and sulfadoxine-pyrimethamine for treatment of uncomplicated Plasmodium falciparum malaria in the Democratic Republic of Congo

Artesunate+amodiaquine and artesunate+sulphadoxine-pyrimethamine for treatment of uncomplicated malaria in Democratic Republic of Congo: a clinical trial with determination of sulphadoxine and pyrimethamine-resistant haplotypes

Surveillance, suivi et évaluation. Un manuel de référence. Organisation Mondiale de la Santé

PfCRT and its role in antimalarial drug resistance

A molecular marker for chloroquine-resistant falciparum malaria

Amodiaquine resistance in Plasmodium falciparum malaria in Afghanistan is associated with the pfcrt SVMNT Allele at Codons 72 to 76

Return of chloroquine antimalarial efficacy in Malawi

Trends in chloroquine resistance marker, Pfcrt-K76T mutation ten years after chloroquine withdrawal in Tanzania

Impacts of antimalarial drugs on Plasmodium falciparum drug resistance markers

The return of chloroquine-susceptible Plasmodium falciparum malaria in Zambia

High prevalence of PfCRT K76T mutation in Plasmodium falciparum isolates in Ghana

High prevalence of pfcrt-CVIET haplotype in isolates from asymptomatic and symptomatic patients in south-central Oromia Ethiopia

A crosssectional survey of Plasmodium falciparum pfcrt mutant haplotypes in the Democratic Republic of Congo

Molecular surveillance of Plasmodium falciparum resistance to artemisinin-based combination therapies in the Democratic Republic of Congo

Molecular surveillance of anti-malarial drug resistance in Democratic Republic of Congo: high variability of chloroquinoresistance and lack of amodiaquinoresistance

Spread of artemisinin resistance in Plasmodium falciparum malaria

A molecular marker of artemisinin-resistant Plasmodium falciparum malaria

World Health Organization/Global Malaria programme

Emergence and clonal expansion of in vitro artemisinin-resistant Plasmodium falciparum kelch13 R561H mutant parasites in Rwanda

Artemisinin resistance-associated K13 polymorphisms of Plasmodium falciparum in Southern Rwanda

Artemisinin-resistant Plasmodium falciparum with high survival rates

Spatial and molecular mapping of Pfkelch13 gene polymorphism in Africa in the era of emerging Plasmodium falciparum resistance to artemisinin: a systematic review

The lack of K13-propeller mutations associated with artemisinin resistance in Plasmodium falciparum in Democratic Republic of Congo (DRC)

Validation of a four-primer real-time PCR as a diagnostic tool for single and mixed Plasmodium infections

Assessment of pfcrt 72-76 haplotypes eight years after chloroquine withdrawal in Kinshasa, Democratic Republic of Congo

A worldwide map of Plasmodium falciparum K13-propeller polymorphisms

Association of mutations in the Plasmodium falciparum (Pf3D7_1343700) with parasite clearance rates after artemisinin-based treatments-a WWARN individual patient data meta-analysis

Mutations in Plasmodium falciparum K13 propeller gene from Bangladesh (2009-2013)

K13-propeller mutations confer artemisinin resistance in Plasmodium falciparum clinical isolates

The impact of antimalarial resistance on the genetic structure of Plasmodium falciparum in the DRC

Re-emergence of chloroquine-sensitive Plasmodium falciparum malaria after cessation of chloroquine use in Malawi

Off-label use of chloroquine and hydroxychloroquine for COVID-19 treatment in Africa against WHO recommendation

Identification of falsified chloroquine tablets in Africa at the time of the COVID-19 pandemic

Geographic patterns of Plasmodium falciparum drug resistance distinguished by differential responses to amodiaquine and chloroquine

Amodiaquine for treating malaria. Cochrane Database Syst Rev

Discordant patterns of genetic variation at two chloroquine resistance loci in worldwide populations of the malaria parasite Plasmodium

The study was conducted in collaboration with the National Malaria Control Programme (NMCP). We would like to thank all patients and their guardians, as well as the staff of all collecting sites for their participation in this study.

The online version contains supplementary material available at https:// doi. org/ 10. 1186/ s12879-022-07112-z.Additional file 1: Table S1 . Plasmodium falciparum K13 nucleotide sequences from 2019-study.Authors' contributions DMY, DMM, PDM, GLM and MPH initiated the study. DMY and NKK carried out the data and sample collection. DMY and RB performed the molecular analysis. PZA carried out the statistical analysis. DMY designed the study and drafted the manuscript. NKK, DMM, RB, PZK, PZA, HNTS, PDM, NS, GLM and MPH contributed to the interpretation of data and they substantively revised the manuscript. All authors read and approved the final manuscript.

This study was supported by the Académie de Recherche et d'Enseignement Supérieur (ARES), Belgium. The funding body had no role in the study design, interpretation of data and in writing and publication of the manuscript.

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

The protocol and the informed consent form were approved by the Ethics Committee of the Faculty of Medicine, University of Kinshasa. All participants involved in the study signed an informed consent form. For participants aged less than 18 years old, the consent form was approved and signed by their parent or guardian.

Not applicable.

The authors declare that they have no competing interests. 

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.