Skip to main content

DNA recovery from archived RDTs for genetic characterization of Plasmodium falciparum in a routine setting in Lambaréné, Gabon



Rapid diagnostic tests (RDTs) have been described as a source of genetic material to analyse malaria parasites in proof-of-concept studies. The increasing use of RDTs (e.g., in focal or mass screening and treatment campaigns) makes this approach particularly attractive for large-scale investigations of parasite populations. In this study, the complexity of Plasmodium falciparum infections, parasite load and chloroquine resistance transporter gene mutations were investigated in DNA samples extracted from positive RDTs, obtained in a routine setting and archived at ambient temperature.


A total of 669 archived RDTs collected from malaria cases in urban, semi-urban and rural areas of central Gabon were used for P. falciparum DNA extraction. Performance of RDTs as a source of DNA for PCR was determined using: (i) amplification of a single copy merozoite surface protein 1 (msp1) gene followed by highly sensitive and automated capillary electrophoresis; (ii) genotyping of the pfcrt gene locus 72–76 using haplotype-specific-probe-based real-time PCR to characterize chloroquine resistance; and, (iii) real-time PCR targeting 18S genes to detect and quantify Plasmodium parasites.


Out of the 669 archived RDTs, amplification of P. falciparum nucleic materials had a success rate of 97% for 18S real-time PCR, and 88% for the msp1 gene. The multiplicity of infections (MOI) of the whole population was 2.6 (95% CI 2.5–2.8). The highest number of alleles detected in one infection was 11. The MOI decreased with increasing age (β = − 0.0046, p = 0.02) and residence in Lambaréné was associated with smaller MOIs (p < 0.001). The overall prevalence of mutations associated with chloroquine resistance was 78.5% and was not associated with age. In Lambaréné, prevalence of chloroquine resistance was lower compared to rural Moyen-Ogooué (β = − 0.809, p-value = 0.011).


RDT is a reliable source of DNA for P. falciparum detection and genotyping assays. Furthermore, the increasing use of RDTs allows them to be an alternative source of DNA for large-scale genetic epidemiological studies. Parasite populations in the study area are highly diverse and prevalence of chloroquine-resistant P. falciparum remains high, especially in rural areas.


Despite being treatable, malaria continues to be one of the major health problems in sub-Saharan Africa, with an estimated 219 million cases and 435,000 deaths in 2017 [1]. The clinical course can vary significantly between individuals, with many determinants remaining to be identified. Effective control and elimination require integrated multilayer strategies, including prompt diagnosis, appropriate chemotherapy, and case management [2].

Reliable, appropriate and timely diagnosis of febrile malaria patients is crucial, especially those caused by Plasmodium falciparum. This is achieved by microscopy of thick blood smears and/or using a malaria rapid diagnostic test (RDT) [2]. In addition, nucleic acid testing (NAT) is used to detect low-density Plasmodium infections in epidemiological studies aimed at surveillance of malaria control and elimination. The DNA template for NAT is usually recovered from anti-coagulant (EDTA)-treated blood or blood spotted on filter paper [3] and sometimes from thick blood smears [4]. Venous fresh blood sampling requires special storage conditions and training in venipuncture. Additionally, the use of venous blood samples leads to somewhat lower diagnostic sensitivity compared to capillary blood [5]. Capillary sampling is usually used for thick blood smears and sampling for RDTs or filter paper. Quality of DNA recovered from thick blood smears is low compared to filter paper and fresh blood [6] and the use of filter paper has to be prospectively planned, since it is not part of routine diagnostics. RDTs support rational and timely use of anti-malarial drugs in field settings, particularly if reliable microscopy is not available [7, 8]. The ease of use, low cost, and performance of RDTs has led to an increase in sales of these supplies in Africa from 240 million in 2015 to 269 million in 2016 [8]. Furthermore, successful amplification of DNA recovered from RDTs was observed [9,10,11,12] and the PCR detection rate in DNA extracted from RDTs is similar to that from filter paper [13]. Consequently, RDTs potentially represent an ideal source for large-scale retrospective analyses of parasite populations.

Parasite genotyping is essential for discrimination of new clones, from existing ones, to predict treatment outcomes (recrudescence vs re-infection). Additionally, it may provide evidence to understand the transmission dynamics in areas where there is perennial transmission of malaria. Natural infections can consist of multiple genetically distinct parasite strains co-infecting a single host. Such co-infections are reported as ‘multiplicity of infection’ (MOI). MOI is an indicator of transmission intensity [14] and may be used for molecular monitoring of any large-scale interventions. Genotyping approaches using conventional PCR targeting length polymorphic markers has been recommended to be used for assessment of efficacy of novel drugs [15]. The performance of these approaches varies between laboratories and mainly relies on the resolution of sizing systems. Despite the limitations: low resolution and time consuming, agarose gel electrophoresis is the most popular method for sizing due to its relatively low cost. A previous study showed capillary electrophoresis has better resolution and accurate estimation of genotypes than gel electrophoresis; this increased precision is critical especially in new anti-malarial drug trials to assess the efficacy of the intervention [16].

In the Republic of Gabon, artemisinin-based combination therapy (ACT) replaced chloroquine (CQ) as the first-line drug for the treatment of uncomplicated P. falciparum malaria in 2005. Since then, two studies have reported on the circulating genotypes based on merozoite surface protein 1 (MSP1) genotyping from Libreville in 2011–2012 [17] and from Oyem-Owendo for years 2008–2009 [18]. For Lambaréné, where there is perennial transmission, data on parasite genotypes exist for years 1995–1996 [19, 20] and CQ resistance for years 2001–2002 [21, 22], before the introduction of ACT. A few studies recently documented a differential distribution of CQ-resistant alleles among different geographical regions within Gabon [23, 24]. Yet, another study concluded that despite the implementation of ACT and the withdrawal of CQ, no significant difference in the prevalence of CQ-susceptible wild-type haplotype CVMNK was observed in parasite isolates collected between 2011 and 2014 [25]. Studies at regular intervals are required for molecular surveillance for anti-malarial drug resistance, including CQ [26].

Against this background, this study utilized a large cohort of individuals inhabiting Lambaréné, Fougamou and the respective adjacent rural areas, and aimed at investigating the performance of archived RDTs as a source of DNA for different molecular tests. Additionally, the study aimed to characterize parasite diversity, measured by MOI, and to determine the distribution of genetic polymorphisms in the P. falciparum CQ resistance transporter (pfcrt) gene locus.


Study site and samples

The study was conducted at the Centre de Recherches Médicales de Lambaréné (CERMEL), Lambaréné, Moyen-Ogooué, Gabon in July 2018. The study region is within the equatorial rainforest and highly endemic for falciparum malaria, which is perennial with little seasonal variation [27]. Local strains of P. falciparum show high levels of resistance against CQ and sulfadoxine-pyrimethamine [22, 28].

RDTs were collected from June 2017 to July 2018 during the screening process of an ongoing clinical trial (NCT03201770) and routine activities. The inclusion and exclusion criteria can be retrieved from the clinical trial registry NCT03201770. Briefly, participants of all ages and both genders were included if they presented with uncomplicated malaria, weight > 5 kg, and provided a signed informed consent. They were excluded in case of hepatic injury, known allergy to study drugs, and pregnant or lactating women.

Participants were tested for malaria using three different RDTs (VIKIA® Malaria Ag Pf/Pan, IMAccess, Lyon, France; Paracheck Pf®, Orchid Biomedical Systems, Goa, India; SD BIOLINE Malaria Ag P.f/Pan Standard Diagnostics Inc, Hagal-Dong, Korea). All RDTs were WHO prequalified to detect P. falciparum antigens with a limit of detection of ≥ 200 parasites/µL. Once testing was done, positive RDT cassettes were stored at ambient temperature in a sealed pouch, until further use. Since the exposure to malaria varies between age groups, the recruited individuals were age-stratified in: (i) children < 5 years; (ii) children between 5 and 18 year; and, (iii) adults > 18 years old. The individuals were inhabitants of Lambaréné (semi-urban), Fougamou (semi-urban) and surrounding rural areas (radius approximately 10 km). Of the 1008 archived RDT cassettes, 669 RDT cassettes that had both demographic data and positive readable test lines (with P. falciparum) were used for further investigation.

A blood sample (with parasitaemia of P. falciparum of 6,840,000 p/mL) obtained from a malaria patient was tenfold serially diluted with malaria-free group O + blood (lowest parasitaemia 68 p/mL). Five milliliter of each dilution (6 dilutions) was spotted on new RDTs (SD BIOLINE Malaria Ag P.f/Pan) and 10 µL was spotted on each circle of filter papers (Whatman™ 903 Protein Saver Card). All 6 RDTs and filter paper (6 circles) were kept at ambient temperature overnight before DNA extraction.

DNA extraction

From RDTs Under sterile conditions, individual RDT cassettes were opened using scissors and forceps and the nitrocellulose strip was removed from the plastic case. After removal of the nitrocellulose strip, any plastic covering on the strip was stripped off. The proximal third of the test strip containing DNA was dissected and was subsequently used for DNA extraction as described elsewhere [29]. To avoid contamination, scissors and forceps were washed first with alcohol followed by DNA AWAY™ solution (Molecular BioProducts, San Diego, USA) each time and dried before opening a new RDT cassette. DNA extraction was done using QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany) following manufacturer’s instruction for filter paper with elution volume of 50 µL. A positive control from laboratory-cultured parasites (PfNF54) and a negative control devoid of parasite DNA was extracted from the RDT cassette.

From filter paper Half of each circle of filter paper containing approximately 5 µL of sample was cut and used for DNA extraction. The DNA extraction was done using QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany) following manufacturer’s instruction for filter paper, elution volume was 50 µL.

Plasmodium detection and quantification

Taqman probe-based Pan-Plasmodium Real-time PCR targeting a highly conserved region of 18S rRNA was used to detect and quantify Plasmodium parasites [30] without using reverse transcriptase (RT). The utilized primer and probes for the assay are described in Additional file 1: Table S1. The assays were performed using TaqMan™ RNA-to-CT™ 1-Step Kit (Thermo Fisher Scientific, Foster City, CA, USA) following manufacturer’s instructions for non-RT PCRs in a LightCycler 480 Instrument II (Roche, Basel, Switzerland). A negative control devoid of parasite DNA and a P. falciparum NF54 positive control were integrated in each run. The results were analysed by LightCycler® 480 SW v1.5.1 software. The success of an amplification is defined by respective Cq values that are smaller or equal to 40.

pfcrt genotyping by Taqman probe-based real-time PCR

The pfcrt genotyping for codons 72–76 was performed using previously described primers and probes [31] and details were provided in Additional file 1: Table S1. To enhance the sensitivity, the pfcrt locus was pre-amplified by conventional PCR and the enriched amplicon was used as template for real-time PCR to detect pfcrt haplotypes. In brief: for conventional PCR, a total reaction volume of 20 µL was constituted by adding 2.5 µL of template to 17.5 µL of reaction mix containing: 1 × PCR buffer, 0.4 µM of each primer, 0.25 mM dNTPs, 1U of Taq polymerase (Qiagen, Hilden, Germany). The thermal conditions were: 94 °C for 5 min followed by 25 cycles of 94 °C for 30 s, 55 °C for 30 s and 72 °C for 30 s, a final annealing at 72 °C for 10 min. This PCR was done using MyCycler (BioRad, Germany). The amplicons from the conventional PCR were used as template for the multiplex real-time PCR using three haplotype specific probes namely CQ-sensitive (CVMNK) and two CQ-resistant genotypes (CVIET and SVMNT) Additional file 1: Table S1. Multiplex PCR was performed on a LightCycler 480 Instrument II using SensiMix™ II Probe Kit (Bioline GmbH., Germany) following manufacturer’s instructions. Each sample was tested in duplicates. DNA extracted from P. falciparum NF54, Dd2 and 7G8 strains were used as positive controls for genotyping CVMNK, CVVIET and SVNMT haplotypes, respectively. Non-template controls as well as positive controls were included in every PCR batch. Amplification was considered as a success or positive when Cq values were equal to or smaller than 40.

Table 1 Prevalence of pfcrt haplotypes, age distribution and msp1 genotyping result across regions

Genotyping of msp1 locus

Amplification of msp1 gene using conventional PCR The msp1 gene was chosen as a target for genotyping in this study. A nested PCR was performed using published pairs of primers [32]. Primary PCR amplified the conserved region in msp1 followed by a nested PCR which amplified the family specific block 2 of the msp1 gene: K1, MAD20 and RO33 (sequences of primers are listed in Additional file 1: Table S1). Primary PCR was set with total volume of reaction of 20 µL containing 5 µL of template mixed with 15 µL of master mix (1× AmpliTaq Gold buffer, 1.5 mM MgCl2, 0.25 mM dNTPs, 1 U of Ampli Taq Polymerase and 250 nM each primer). The thermal cycling reaction was performed using a MyCycler (BioRad, Germany) for 10 min at 94 °C followed by 35 cycles of 94 °C in 15 s, 58 °C for 30 s, 72 °C for 1 min and 72 °C for 10 min for the final extension. The PCR product of the primary PCR was used as template for nested PCR. For nested PCR, three sets of reaction mixes were made with the total volume of reaction of 20 µL for three pairs of primers. The reaction mixes contain 2.5 µL of template, other components were added with the same concentration as used in the primary PCR. Thermal conditions of the nested PCRs were also similar to that of the primary PCR except the annealing temperature was at 61 °C. DNA of P. falciparum NF54, Dd2, and 7G8 strains were used as positive controls for three msp1 families K1, MAD20, and RO33, respectively.

Amplicon sizing using capillary electrophoresis Amplicons were sized using an automated capillary gel electrophoresis QIAxcel Advance system (QIAGEN, Hilden, Germany) according to manufacturer’s instructions. In brief, the QIAxcel DNA High Resolution Kits (cat. no. 929002) were used. Samples were run on 96-well plates, with OM400 protocol using QX DNA Size Marker 50–800 bp (50 µL) v2.0 (Qiagen, cat. no. 929561) and QX Alignment Marker 15 bp/1 kb (Qiagen, cat. no. 929521). Interpretation was done with software QIAxcel ScreenGel v1.5.0. PCR product plates with positive controls that were positive with a single peak (NF54: 241 ± 3 bp, Dd2: 205 ± 3 bp and 7G8: 153 ± 3 bp) were included in the analysis. Each peak represents an allele and was excluded from the analysis if smaller than 100 bp or contributes less than 10% of total peak height. Within one sample, peaks with less than 6 bp difference were interpreted as one allele and average size was used for further analyses.

Statistical analysis

All statistical analyses were computed with R software (version 3.5.1). Where applicable, 95% confidence intervals are given. Parasitaemia was calculated by extrapolating cycle threshold (Ct) values and standard curves based on linear regression analysis (assuming the amount of blood spotted on each RDT and calibrators was identical: 5 µL). All positive PCRs with Ct values > 35 and < 40 were considered negative (0 parasite/mL) using a regression model and excluded from analysis. Paired sample t-test was used to compare Ct-values of PCRs run on RDT and filter paper serial dilutions.

The MOI was defined as the mean number of P. falciparum genotypes per infected individual. The MOI was calculated as the proportion of the total number of P. falciparum msp1 genotypes and the total number of PCR positive isolates. Isolates with more than one genotype were considered as a polyclonal infection while the presence of a single allele was considered as a monoclonal infection.

Mean number of alleles detected was displayed, but comparisons of MOI between age or location groups were done using the non-parametric Kruskal–Wallis H and Wilcoxon rank-sum tests. These tests were also used to compare log10 transformed parasitaemia density between age groups. Logistic regression models were used in this study: to check whether there was an effect of age and location on MOI (Poisson regression model), on pfcrt genotype (Binomial regression model) and to analyse the relationship between time of sample storage and positivity of qPCR (Binomial regression model). A linear regression model was used to assess the relationship between age and parasite density. A two-sided p-value of < 0.05 was considered statistically significant.


Demographic details of studied population

Of the 669 individuals from whom the positive RDTs were collected, 332 (50%) were female and the median age was 8 years (interquartile range (IQR): 4–15). There were 185 (28%) children younger than 5 years, 346 (52%) aged between 5 and 18 years and 138 (21%) aged > 18 years. A total of 125 (18.7%) samples were from Lambaréné, 122 (18%) from Fougamou and 422 (63%) from rural areas as shown in Table 1. The mean time from sample collection to extraction was 113.2 days (range 9–231 days).

Performance of RDTs as a source of DNA for quantitative PCR

Two series of diluted samples including 6 RDTs and 6 circles of filter paper were run on 18S real-time PCR in triplicates in order to compare the performance of these DNA sources. PCRs were positive with all samples including lowest parasitaemia samples (70 parasite/mL). The PCR results showed no differences in Ct-values between two groups of samples: RDT and filter paper (df = 5, p-value = 0.1) (Fig. 1).

Fig. 1

Comparison of Ct-values for different parasitaemia between two groups of samples RDT (red) and filter paper (blue). The PCRs were done in triplicate. The dots represent the average Ct-values of triplications. The smooth lines visualized the linear models between Ct-values and log10 transformed parasitaemia of RDT group (red, adjusted R2: 0.9944) and filter paper group (blue, adjusted R2: 0.9892)

Median log10 parasitaemia in the age group < 5 years was 3.50 (IQR: 2.5–4.5), in the other two age groups 5 to 18 years and > 18 years, were 2.91 (IQR: 2.2–3.9) and 2.75 (IQR: 1.8–3.6), respectively (Fig. 2) and parasitaemia in the first group was significantly higher than the two latter ones. PCR amplification success rates were 97, 88.5 and 87% for 18S, msp1 and pfcrt, respectively (Table 2). Longer time of sample storage did not lower the performance of RDTs as a source of DNA for 18 s PCR (p-value = 0.78) (Table 2).

Fig. 2

Parasitaemia across age groups. Y axis represents the log 10 transformed parasitaemia. The significance levels were shown as: significant (***: < 0.001, ****: < 0.0001) and non-significant (ns)

Table 2 Proportion of positivity of molecular assays using DNA extracted from RDTs

msp1 genotyping reveal high population diversity of Plasmodium falciparum

Out of 592 isolates positive for msp1, the K1, MAD20 and RO33 alleles were observed at the following proportions K1: 75% (n = 446), MAD20: 47% (n = 275) and RO33 51% (n = 305). The highest number of alleles detected in a single infection was 11. The MOI of the studied population was: 2.6 (95% CI 2.5–2.8).

Based on msp1 genotyping, 27.5% (n = 163) were monoclonal infections Table 1. The major proportion of infections were with msp1/K1 family, whereas polyclonal infections were most frequently with msp1/K1 + RO33 and K1 + Ro33 + MAD20, Table 1. Allele sizes of msp1 were 133–374 bp (K1), 133–311 bp (MAD20) and 111–258 bp (RO33) (Fig. 3). There was no distinct allelic pattern for a specific age group. There were 8 samples from Lambaréné carrying K1 alleles (125–155 bp), which were not observed in other regions. The mean log parasitaemia of the 8 samples from Lambaréné was 4.5 (95% CI 3.8–5.2), which is higher than the mean log parasitaemia of the population.

Fig. 3

Frequency of alleles across age groups and region. The upper lane: frequency of msp1 alleles across age groups. The lower lane: msp1 alleles frequency across 3 studied areas. RR rural areas, FGM Fougamou, LA Lambaréné

The MOI for ages < 5 years was 2.8 (95% CI 2.6–3.1); for ages 5 to 18 years 2.7 (95% CI 2.5–2.9) and for ages > 18 years 2.1 (95% CI 1.9–2.4). The MOI for Lambaréné was 2.0 (95% CI 1.8–2.3), Fougamou was 2.8 (95% CI 2.5–3.1) and rural areas was 2.8 (95% CI 2.6–2.9) (Fig. 4). Both age and location were observed to associate with MOI as computed in the Poisson regression model (age: β = − 0.0046, p-values = 0.019 and location: β = − 0.314, p = 0.0004).

Fig. 4

MOI across studied regions (a) and age groups (b). MOI was shown as number of allele per individual. The box plots showed the medians and inter quartile ranges (IQR). Number of allele (countable number) is visualized by jitter point. The significance levels were shown as: significant (***: < 0.001, ****: < 0.0001) and non-significant (ns). RR rural areas, FGM Fougamou, LA Lambaréné

Prevalence of CQ-resistant haplotypes

In the group aged < 5 years, 26% of the samples were CQ-sensitive, 48% were resistant (CVIET genotype) and 27% were mixed infections carrying both sensitive and resistant strains (Table 1). The prevalence of sensitive, CVIET and mixed strains in other groups was 18.8, 44.2 and 37%, respectively for the age group from 5 to 18 years; 22.6, 47.0 and 30.4%, respectively for adults (Fig. 5). The overall prevalence of CQ-resistant strains was 78.5%. There was no statistically significant difference in the prevalence of CQ-resistant strains between age groups (p-value = 0.2). A statistically significant difference was observed in the prevalence of CQ-resistant genotypes between Lambaréné and the other areas (Fougamou and rural areas), 65.7 vs 81.3% (p-value = 0.0017) and there was no difference observed between Fougamou and rural areas Table 1. This result was further confirmed by logistic regression, in which living in Lambaréné was significantly associated with resistance of P. falciparum (β = − 0.809, p-value = 0.011 after correction for age).

Fig. 5

Prevalence of Pfcrt haplotypes across age groups. Mix infections contain both sensitive and resistant genotypes. Haplotype SVMNT was not detected in this population


This study retrospectively investigated the utility of archived RDTs with a large sample size for molecular characterization of P. falciparum population diversity, MOI and the distribution of CQ-resistant haplotypes. The main findings of this study are a high amplification success rate of PCRs using DNA extracted from archived RDTs which were not initially stored for this purpose, diversity of P. falciparum parasites and the persistence of CQ-resistant strains.

Results from pan-Plasmodium PCR showed that archived RDTs are a reliable source of DNA for real-time PCR, comparable to filter paper [13]. PCR using DNA extracted from both sources could detect Plasmodium parasites at a low density. Even approximate quantification of malaria parasites from RDT samples is possible, which adds important information beyond the RDT result. However, quantification can be biased by the unequal blood volume spotted on RDT strips during RDT preparation [33]. The log transformed parasitaemia result in this work, and therefore was used only in comparison analyses with the assumption that blood volumes were similar (on field samples and controls). The storage duration of archived samples (up to 6 months) did not compromise either with DNA quality or the amplification success rate. Age of the host and the intensity of exposure to the parasites are the major factors which influence the risk of developing parasitaemia and disease by inducing and regulating the host immunity to malaria, therefore higher parasitaemia in children than in adults has been observed [34].

There have been many approaches used to estimate the diversity of malaria parasite populations; those based on PCR and targeting the length polymorphic genes have been commonly used. The performance of these approaches relies mainly on the resolution of the sizing system and could be enhanced by using capillary electrophoresis. This study is the first one performed in Lambaréné using the high resolution, automated capillary gel electrophoresis by QIAxcel Advanced System to investigate the genetic diversity of P. falciparum msp1 gene in field isolates. This approach provides a highly discriminatory power since the two fragments can be distinguished by three base-pairs difference, and is able to detect up to 11 alleles in an infection.

MSP1 of P. falciparum is a major protein, which plays an important role in the invasion of the parasite into the host’s red blood cell. Genotyping malaria parasites based on merozoite surface genes may contribute to the understanding of selection pressure. This work shows the differences in multiplicity of malaria infections between age groups and geographic regions and a high level of polymorphisms of msp1 in the studied population. The overall MOI in the studied population was 2.6 (2.0 for Lambaréné and 2.8 for Fougamou and rural areas), which is similar to previous studies conducted in Gabon [35] and lower compared to 2000 [36], before the application of ACT as the first-line therapy. Innate immunity against malaria developed after exposure to different P. falciparum infections, which directly related to age and transmission intensity, therefore, allelic diversity can decrease with increasing age of the infected individual. After being corrected for the effect of age, difference in geographic location has an effect on MOI. One reason could be that less diversity of parasites within the mosquito or different drug pressures might cause disappearance and re-appearance of parasite clones in patients [36].

Product of the msp1/K1 and msp1/MAD20 gene families have been identified as candidates for P. falciparum vaccines [37]. In this studied population, K1 was predominantly detected followed by RO33 and MAD20 families. RO33 appeared to be not as polymorphic as K1 and MAD20 families and there were no regional or age specific alleles. This finding is consistent with results from other studies performed in different settings and geographic regions [36, 38]. The predominant allele size of K1 P. falciparum was around 240 bp (26.7%), similar to the length of the msp1 gene of strain NF54 (241 bp), which is being used as a vaccine candidate in some clinical trials in Lambaréné.

Gabon is a country where prevalence of biomarkers of CQ resistance remained high although a partial reduction was noticed in the last decade [39]. In Moyen-Ogooué, CQ-resistant haplotypes were significantly less common in Lambaréné and surrounding areas (radius approximately 10 km) than in Fougamou and rural areas. Prevalence of CQ-resistant haplotypes was age independent. All of the CQ-resistant haplotypes were CVIET (triple mutation at codons 74, 75 and 76), no sample carried haplotype SVMNT (double mutation at codons 72 and 76). This result highlights the need for further analyses to better understand the mechanism behind the maintained high prevalence of CQ resistance. Amodiaquine has been linked to the persistence of CQ-resistant strains since its use in combination with artesunate may have resulted in continued selection of mutant pfcrt haplotypes [39].

One limitation of this work is that it did not include negative field RDTs, therefore infections with low parasite density, which were negative by RDT could not be analysed. That limitation could hinder documenting the existence of circulating parasite genotypes which might also contribute towards transmission intensity and drug resistance.


This study showed that RDTs are a reliable source of DNA for P. falciparum detection and genotyping analyses using PCR and capillary electrophoresis. The parasite population in the studied regions is highly diverse and prevalence of pfcrt in Gabon remains high. Prevalence of molecular markers of CQ resistance is age independent, but regionally dependent, which may indicate the effect of transmission intensity on the prevalence of resistance rather than acquired immunity. The work also suggests the possibility of quantifying Plasmodium parasites in RDTs, especially when the volume of blood used is uniform.

Availability of data and materials

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



artemisinin-based combination therapy


Centre de Recherches Médicales de Lambaréné




cycle threshold


deoxyribonucleic acid


deoxynucleotide Triphosphates


ethylenediaminetetraacetic acid






multiplicity of infection

msp1 :

merozoite surface protein 1


nucleic acid test


polymerase chain reaction


Plasmodium falciparum chloroquine resistance transporter


rapid diagnostic test


rural region


reverse transcriptase


World Health Organization


  1. 1.

    WHO. World Malaria Report. Geneva: World Health Organization; 2018. 238 p.

  2. 2.

    WHO. Guidelines for the treatment of malaria. 3rd Edn. Geneva: World Health Organization; 2015.

  3. 3.

    Kain KC, Lanar DE. Determination of genetic variation within Plasmodium falciparum by using enzymatically amplified DNA from filter paper disks impregnated with whole blood. J Clin Microbiol. 1991;29:1171–4.

    CAS  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Edoh D, Steiger S, Genton B, Beck HP. PCR amplification of DNA from malaria parasites on fixed and stained thick and thin blood films. Trans R Soc Trop Med Hyg. 1997;91:361–3.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Mischlinger J, Pitzinger P, Veletzky L, Groger M, Zoleko-Manego R, Adegnika AA, et al. Use of capillary blood samples leads to higher parasitemia estimates and higher diagnostic sensitivity of microscopic and molecular diagnostics of malaria than venous blood samples. J Infect Dis. 2018;218:1296–305.

    Article  PubMed  Google Scholar 

  6. 6.

    Scopel KKG, Fontes CJF, Nunes AC, Horta MF, Braga EM. Low sensitivity of nested PCR using Plasmodium DNA extracted from stained thick blood smears: an epidemiological retrospective study among subjects with low parasitaemia in an endemic area of the Brazilian Amazon region. Malar J. 2004;3:8.

    Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Wongsrichanalai C, Barcus MJ, Muth S, Sutamihardja A, Wernsdorfer WH. A review of malaria diagnostic tools: microscopy and rapid diagnostic test (RDT). Am J Trop Med Hyg. 2007;77:119–27.

    Article  PubMed  Google Scholar 

  8. 8.

    WHO. Rapid diagnostic tests. Geneva: World Health Organization; 2018.

  9. 9.

    Veron V, Carme B. Recovery and use of Plasmodium DNA from malaria rapid diagnostic tests. Am J Trop Med Hyg. 2006;74:941–3.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Papa Mze N, Ndiaye YD, Diedhiou CK, Rahamatou S, Dieye B, Daniels RF, et al. RDTs as a source of DNA to study Plasmodium falciparum drug resistance in isolates from Senegal and the Comoros Islands. Malar J. 2015;14:373.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Ishengoma DS, Lwitiho S, Madebe RA, Nyagonde N, Persson O, Vestergaard LS, et al. Using rapid diagnostic tests as source of malaria parasite DNA for molecular analyses in the era of declining malaria prevalence. Malar J. 2011;10:6.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Cnops L, Boderie M, Gillet P, Van Esbroeck M, Jacobs J. Rapid diagnostic tests as a source of DNA for Plasmodium species-specific real-time PCR. Malar J. 2011;10:67.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Morris U, Aydin-Schmidt B, Shakely D, Mårtensson A, Jörnhagen L, Ali AS, et al. Rapid diagnostic tests for molecular surveillance of Plasmodium falciparum malaria -assessment of DNA extraction methods and field applicability. Malar J. 2013;12:106.

    Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Kolakovich KA, Ssengoba A, Wojcik K, Tsuboi T, Al-Yaman F, Alpers M, et al. Plasmodium vivax: favored gene frequencies of the merozoite surface protein-1 and the multiplicity of infection in a malaria endemic region. Exp Parasitol. 1996;83:11–8.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    WHO. Methods and techniques for clinical trials on antimalarial drug efficacy: genotyping to identify parasite populations. Geneva, World Health Organization, 2008;.

  16. 16.

    Gupta V, Dorsey G, Hubbard AE, Rosenthal PJ, Greenhouse B. Gel versus capillary electrophoresis genotyping for categorizing treatment outcomes in two anti-malarial trials in Uganda. Malar J. 2010;9:19.

    Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Bouyou-Akotet MK, M’Bondoukwé NP, Mawili-Mboumba DP. Genetic polymorphism of merozoite surface protein-1 in Plasmodium falciparum isolates from patients with mild to severe malaria in Libreville, Gabon. Parasite. 2015;22:12.

    Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Mawili-Mboumba DP, Mbondoukwe N, Adande E, Bouyou-Akotet MK. Allelic diversity of msp1 gene in Plasmodium falciparum from rural and urban areas of Gabon. Korean J Parasitol. 2015;53:413–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Kun JF, Schmidt-Ott RJ, Lehman LG, Lell B, Luckner D, Greve B, et al. Merozoite surface antigen 1 and 2 genotypes and rosetting of Plasmodium falciparum in severe and mild malaria in Lambarene, Gabon. Trans R Soc Trop Med Hyg. 1998;92:110–4.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Mayengue PI, Luty AJF, Rogier C, Baragatti M, Kremsner PG, Ntoumi F. The multiplicity of Plasmodium falciparum infections is associated with acquired immunity to asexual blood stage antigens. Microbes Infect. 2009;11:108–14.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Borrmann S, Binder RK, Adegnika AA, Missinou MA, Issifou S, Ramharter M, et al. Reassessment of the resistance of Plasmodium falciparum to chloroquine in Gabon: implications for the validity of tests in vitro vs. in vivo. Trans R Soc Trop Med Hyg. 2002;96:660–3.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Ramharter M, Wernsdorfer WH, Kremsner PG. In vitro activity of quinolines against Plasmodium falciparum in Gabon. Acta Trop. 2004;90:55–60.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Lekana-Douki JB, Dinzouna Boutamba SD, Zatra R, Zang Edou SE, Ekomy H, Bisvigou U, et al. Increased prevalence of the Plasmodium falciparum Pfmdr1 86 N genotype among field isolates from Franceville, Gabon after replacement of chloroquine by artemether–lumefantrine and artesunate–mefloquine. Infect Genet Evol. 2011;11:512–7.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Zatra R, Lekana-douki JB, Lekoulou F, Bisvigou U, Ngoungou EB, Ndouo FST. In vitro antimalarial susceptibility and molecular markers of drug resistance in Franceville, Gabon. BMC Infect Dis. 2012;12:307.

    Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Voumbo-Matoumona DF, Kouna LC, Madamet M, Maghendji-Nzondo S, Pradines B, Lekana-Douki JB. Prevalence of Plasmodium falciparum antimalarial drug resistance genes in Southeastern Gabon from 2011 to 2014. Infect Drug Resist. 2018;11:1329–38.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Takala-harrison S, Laufer MK. Antimalarial drug resistance in Africa: key lessons for the future. Ann N Y Acad Sci. 2015;1342:62–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Manego RZ, Mombo-Ngoma G, Witte M, Held J, Gmeiner M, Gebru T, et al. Demography, maternal health and the epidemiology of malaria and other major infectious diseases in the rural department Tsamba-Magotsi, Ngounie Province, in central African Gabon. BMC Public Health. 2017;17:130.

    Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Mombo-Ngoma G, Oyakhirome S, Ord R, Gabor JJ, Greutélaers KC, Profanter K, et al. High prevalence of dhfr triple mutant and correlation with high rates of sulphadoxine-pyrimethamine treatment failures in vivo in Gabonese children. Malar J. 2011;10:123.

    Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    WWARN. DNA Extraction: Preparing Rapid Diagnostic Tests. 2014.

  30. 30.

    Mordmüller B, Surat G, Lagler H, Chakravarty S, Ishizuka AS, Lalremruata A, et al. Sterile protection against human malaria by chemoattenuated PfSPZ vaccine. Nature. 2017;542:445–9.

    Article  PubMed  Google Scholar 

  31. 31.

    Bushman M, Morton L, Duah N, Quashie N, Abuaku B, Koram KA, et al. Within-host competition and drug resistance in the human malaria parasite Plasmodium falciparum. Proc R Soc B Biol Sci. 2016;283:20153038.

    Article  Google Scholar 

  32. 32.

    Messerli C, Hofmann NE, Beck HP, Felger I. Critical evaluation of molecular monitoring in malaria drug efficacy trials and pitfalls of length-polymorphic markers. Antimicrob Agents Chemother. 2017;61:1–11.

    Article  Google Scholar 

  33. 33.

    Luchavez J, Lintag ME, Coll-Black M, Baik F, Bell D. An assessment of various blood collection and transfer methods used for malaria rapid diagnostic tests. Malar J. 2007;6:149.

    Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Doolan DL, Dobaño C, Baird JK. Acquired immunity to malaria. Clin Microbiol Rev. 2009;22:13–36.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Ndong Ngomo JM, M’Bondoukwe NP, Yavo W, Bongho Mavoungou LC, Bouyou-Akotet MK, Mawili-Mboumba DP. Spatial and temporal distribution of Pfmsp1 and Pfmsp2 alleles and genetic profile change of Plasmodium falciparum populations in Gabon. Acta Trop. 2018;178:27–33.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Aubouy A, Migot-Nabias F, Deloron P. Polymorphism in two merozoite surface proteins of Plasmodium falciparum isolates from Gabon. Malar J. 2003;2:12.

    Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Conway DJ, Cavanagh DR, Tanabe K, Roper C, Mikes ZS, Sakihama N, et al. A principal target of human immunity to malaria identified by molecular population genetic and immunological analyses. Nat Med. 2000;6:689–92.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Mayengue P, Ndounga M, Malonga F, Bitemo M, Ntoumi F. Genetic polymorphism of merozoite surface protein-1 and merozoite surface protein-2 in Plasmodium falciparum isolates from Brazzaville, Republic of Congo. Malar J. 2011;10:276.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Frank M, Lehners N, Mayengue PI, Gabor J, Dal-Bianco M, Kombila DU, et al. A thirteen-year analysis of Plasmodium falciparum populations reveals high conservation of the mutant pfcrt haplotype despite the withdrawal of chloroquine from national treatment guidelines in Gabon. Malar J. 2011;10:304.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


We would like to thank Lidwine Badjina (project coordinator—CERMEL) and PYRAMAX study team at CERMEL for their help with sample collection and management of sample’s information.

All P. falciparum strains (3D7, NF54, Dd2, 7G8) used as controls for this study were obtained from Malaria Research and Reference Reagent Resource Center (MRR4). TPV is a member of CANTAM (EDCTP-RegNet2015-1045) and PANDORA-ID-Net (RIA2016E-1609) networks. We acknowledge support by Deutsche Forschungsgemeinschaft and Open Access Publishing Fund of University of Tübingen.


Not applicable.

Author information




GMN, TTN and BN initiated the study and generated the study protocol. Samples were collected by RZM, LBDM, EK. TTN, BN and AL conducted the study with contribution of all authors. Data was analyzed by TTN. GMN, RZM, LBDM, AAA, STA, BL participated to the field study. TTN, TPV interpreted the data and drafted the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to The Trong Nguyen or Benjamin Mordmüller.

Ethics declarations

Ethics approval and consent to participate

The study was approved by Institutional Ethics Committee of CERMEL. Signed informed written consent and assent where applicable were obtained from the screened individuals.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary information

Additional file 1: Table S1.

Primer and probe sequences.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Nguyen, T.T., Nzigou Mombo, B., Lalremruata, A. et al. DNA recovery from archived RDTs for genetic characterization of Plasmodium falciparum in a routine setting in Lambaréné, Gabon. Malar J 18, 336 (2019).

Download citation


  • Malaria
  • Plasmodium falciparum
  • RDT
  • msp1
  • pfcrt
  • Gabon
  • Capillary electrophoresis