Skip to main content
  • Research article
  • Open access
  • Published:

A barcode of multilocus nuclear DNA identifies genetic relatedness in pre- and post-Artemether/Lumefantrine treated Plasmodium falciparum in Nigeria

Abstract

Background

The decline in the efficacy of artemisinin-based combination treatment (ACT) in some endemic regions threatens the progress towards global elimination of malaria. Molecular surveillance of drug resistance in malaria-endemic regions is vital to detect the emergence and spread of mutant strains.

Methods

We observed 89 malaria patients for the efficacy of artemether-lumefantrine for the treatment of uncomplicated Plasmodium falciparum infections in Lagos, Nigeria and determined the prevalence of drug resistant strains in the population. Parasite clearance rates were determined by microscopy and the highly sensitive var gene acidic terminal sequence (varATS) polymerase chain reaction for 65 patients with samples on days 0, 1, 3, 7, 14, 21 and 28 after commencement of treatment. The genomic finger print of parasite DNA from pre- and post-treatment samples were determined using 24 nuclear single nucleotide polymorphisms (SNP) barcode for P. falciparum. Drug resistance associated alleles in chloroquine resistance transporter gene (crt-76), multidrug resistance genes (mdr1–86 and mdr1–184), dihydropteroate synthase (dhps-540), dihydrofolate reductase (dhfr-108) and kelch domain (K-13580) were genotyped by high resolution melt analysis of polymerase chain reaction (PCR) fragments.

Results

By varATS qPCR, 12 (18.5%) of the participants had detectable parasite DNA in their blood three days after treatment, while eight (12.3%) individuals presented with genotypable day 28 parasitaemia. Complexity of infection (CoI) was 1.30 on day 0 and 1.34 on day 28, the mean expected heterozygosity (HE) values across all barcodes were 0.50 ± 0.05 and 0.56 ± 0.05 on days 0 and 28 respectively. Barcode (π) pairwise comparisons showed high genetic relatedness of day 0 and day 28 parasite isolates in three (37.5%) of the eight individuals who presented with re-appearing infections. Crt-76 mutant allele was present in 38 (58.5%) isolates. The mdr1–86 mutant allele was found in 56 (86.2%) isolates. No mutation in the K-13580 was observed.

Conclusions

Persistence of DNA-detectable parasitaemia in more than 18% of cases after treatment and indications of genetic relatedness between pre- and post-treatment infections warrants further investigation of a larger population for signs of reduced ACT efficacy in Nigeria.

Peer Review reports

Background

Malaria caused by Plasmodium falciparum infection is a significant public health problem in Nigeria where over 30% of the global burden of the disease is borne despite recent mass scale-up of intervention measures [1]. Artemisinin-based combination therapy (ACT) is recommended as first line treatment against uncomplicated malaria and it remains a critical tool in malaria control programmes. The most commonly used ACT, Artemether-Lumefantrine (AL), remains largely efficacious in sub-Saharan Africa [1]. The efficacy of AL has been attributed to the fast action of artemisinin and the resultant rapid reduction in the density of asexual malaria parasites, while lumefantrine acts to prevent recrudescence [2]. However, resistance to artemisinin-based drugs has been reported in many endemic regions in South-East Asia [3]. Although evidence of resistance among African populations of Plasmodium falciparum remains scarce, there have been emerging cases of potentially adaptive strains and slow parasite clearance following treatment [4, 5]. These parasites are often below the threshold of detection and are rarely re-treated, thus contributing to disease transmission and spread of resistance [6]. Accurate detection, quantification and characterization of such parasite isolates after treatment is necessary for the proper evaluation and monitoring of drug efficacy.

One recent technique for the ultra-sensitive detection of low parasitaemia involves the amplification of the var gene family present in the sub-telomere of the parasite [7]. Each parasite isolate comprises about 50–150 var genes [8], which possess acidic terminal sequence (ATS) with well-conserved domains that are targeted. Once detected, parasite isolates before and after treatment can be genetically fingerprinted with a single nucleotide polymorphism (SNP) barcode, which is another sensitive assay that requires only a small amount of input DNA and can be used to distinguish recrudescence from re-infection [9]. Prompt detection of mutations in loci implicated in drug resistance such as multidrug resistance gene (mdr1), chloroquine resistance transporter (crt), dihydrofolate reductase (dhfr), dihydropteroate synthase (dhps) and kelch-13 (K13) genes provides molecular evidence of antimalarial resistant parasites [10, 11]. For accurate surveillance of ACT resistance in a malaria-endemic setting as Nigeria, there is a need for continuous therapeutic efficacy studies using sensitive techniques that enable early detection of tolerant parasite strains. In this study, we assessed the efficacy of artemether-lumefantrine, determined genetic relatedness of parasites isolated before and after treatment and the prevalence of markers of drug resistance in the parasite population.

Methods

Study design and sample collection

This observational study of a 28-day follow-up was conducted in Ijede and Agbowa General Hospitals (GH) in Ikorodu, a peri-urban settlement in Lagos, South-West Nigeria, from August – November, 2016 following the revised WHO drug efficacy protocol [12]. Ijede GH is one of the sentinel sites used by the country’s Federal Ministry of Health (FMoH) for antimalarial drug efficacy trials. Individuals (1–70 years) presenting at the hospital with symptoms of uncomplicated malaria were screened for Plasmodium falciparum, first by rapid diagnostic test (RDT), followed by microscopy (Olympus BX53M, UK). The recruitment criteria included fever in the preceding two days, no history of antimalarial intake in the previous four weeks and P. falciparum parasitaemia > 2000/μl of blood on presentation [12]. Individuals with signs of complicated malaria or mixed infections were censored and those that developed concurrent ailments and complications during follow-up were promptly attended to by clinicians and hospital staff and therefore withdrawn from the study. Confirmed cases of uncomplicated malaria were treated with appropriate doses of AL (Coartem, Novartis, Switzerland) according to national policy on malaria treatment [13]. Based on an expected maximum treatment failure rate of 5% in a related population [14] and the desired confidence level (95%) and precision (5%), a minimum of 73 patients were required [12]. An additional 20% was added as provision for patients lost to follow up. Participants and/or their parents/guardians were instructed on how and when to take the evening dose. They were also advised to take fatty meals before or after drug intake to ensure proper absorption of the drug. Following consent, participants were followed up on days 1, 3, 7, 14, 21 and 28 post-treatment (Fig. 1). Confirmation of drug administration at home was done during the follow-up visits by tablet count. Patients were advised to return on any day during the follow-up period if symptoms returned and not to wait for our next scheduled visit day.

Fig. 1
figure 1

Study flow chart

The primary endpoint was clinical and parasitological response within 28 days of commencement of treatment. The secondary endpoint was prevalence of drug resistance markers and genetic relatedness of pre- and post-treatment infections. Participants failing treatment, defined as the reappearance of parasitaemia after an initial complete clearance or failure of complete clearance on day 28, were retreated with AL whether or not they developed symptoms and after samples had been obtained for microscopy and molecular analyses. Before the first treatment was administered, dried blood spots from finger pricks were collected from individuals who gave informed consent and/or assent and also met the recruitment criteria.

varATS qPCR assessment of antimalarial treatment

Plasmodium falciparum DNA was extracted from pre- and post-treatment dried blood spots using the QiaAmp DNA minikit (Qiagen, Germany). The var gene acidic terminal sequence (varATS) quantitative PCR was used to detect multi-copy genomic sequences of low-density infections in follow-up samples [7]. The primer/probe sequences and the cycling conditions are described in Additional file 1: Table S1. Briefly, 1 μl of PCR water, 10 μl of 2× Taqman Universal PCR Mastermix (Applied Biosystems, New Jersey, USA), 1.6 μl of 10uM forward and reverse primers, 0.8 μl of 10 μM probe and 5 μl of parasite DNA were vortexed and run on CFX 96 Touch™ Real-Time System (Bio-Rad Laboratories, CA, USA).

Genotyping of antimalarial drug resistance SNPs

The prevalence of antimalarial drug resistance markers was determined using High Resolution Melt analysis of PCR fragments and probes targeting SNPs [15]. Six SNPs marking variant codons for antimalarial drug resistance were analysed: Chloroquine resistance transporter (crt-76), multidrug resistance (mdr1–86 and mdr1–184), dihydropteroate synthase (dhps-540), dihydrofolate reductase (dhfr-108) and Kelch-13 (K-13580), (see Additional file 1: Table S2 for primer sequences and reaction conditions). Individual samples were traced through time points and the phenotype-genotype associations were profiled for day 0 and day 28 parasites. SNPs were genotyped for all samples positive for parasite DNA by PCR before treatment and on follow-up days.

Molecular SNP barcoding

Molecular fingerprinting and changes in infection complexity before and after drug application was determined by SNP barcode assay. Day 0 and 28 parasitaemia barcodes were compared to determine if there was failure to clear the original parasites by the drug (recrudescence) or re-infection by other parasite strains. A panel of twenty-four unlinked neutral SNP barcodes with minor allele frequencies (MAF) ranging between 0.2–0.5 and representatively spread across the parasite genome (Additional file 1: Table S3) were used. Assays were performed as described by Daniels et al. [9], using the same primer and probe sequences. For pre- and post-treatment parasites, pairwise analysis, barcodes with a mixed allele call for more than one locus were considered as polygenomic and omitted from the analysis. Samples missing allele calls for more than five SNPs were also eliminated from the analysis.

Data and genetic analyses

GraphPad Prism 7 (CA, USA) software was used for statistical analyses. Parasite clearance rates were determined by the Parasite Clearance Estimator (PCE) tool (http://www.wwarn.org). The varATS qPCR survival estimates of parasites at different time-points were determined using Kaplan-Meier survival analysis. Association between age and parasite counts on days 0, 3 and 28 was assessed using Spearman’s correlation coefficient (rs). The relationships between enrollment parasite density on day 0 and parasitaemia on days 3 and 28 were also determined using Spearman’s correlation coefficient. Bonferroni correction was used to adjust multiple testing of age and parasitaemia on days 0, 3 and 28. To determine the genetic identity of each pair of samples on day 0 and day 28, we calculated pairwise barcode similarity index (π), counting ambiguous or missing calls as mismatches [16]. To avoid biasing the percent similarity because of matching major alleles, we limited the percent similarity calculations to include only sites where the minor alleles were present in pairwise barcodes of day 0 and day 28 isolates. The similarity index, π, was calculated as the number of sites where the minor allele matched divided by the total number of sites where the minor allele was present. π has a potential range from 0 (no genetic relatedness) to 1(high genetic relatedness). The allele frequencies per SNP barcode, allelic diversity in pre-treatment and day 28 parasite populations and allele frequencies per barcode per population were calculated using GENALEX 6 [17]. Complexity of infection per SNP barcode was estimated using COIL [18]. Allelic diversity was calculated for each of SNP barcodes based on the allele frequencies, using the formula Nei’s unbiased genetic diversity index for expected heterozygosity, HE = [\( \left(\frac{n}{n-1}\right)\left(1-\sum {P}^2\right)\Big] \), where n is the number of isolates analyzed and P represents the frequency of each allele of the SNP. HE provides an indication of the probability that two samples will be different. It has a potential range from 0 (no allele diversity) to 1 (high allelic diversity). All statistical assumptions were made at 95% confidence limit.

Results

Baseline information on the study participants

In this study, 561 individuals were screened, 150 of them tested positive for falciparum malaria. Eighty-nine patients met the inclusion criteria, were enrolled and followed up. Sixty-five participants had complete sample information for days 1, 3, 7, 14, 21 and 28 after commencement of treatment and were used for genetic analyses. The demographic profiles of the participants are presented in Table 1. Participants were stratified into three age groups: 1-5 yrs.; 6–14 years and ≥ 15 years old (Table 1).

Table 1 Summary of participant characteristics at baseline

Parasite clearance rates

Microscopic analysis showed that none of the participants presented with parasitaemia on day 3 while one of the participants presented with parasitaemia on day 28 (Fig. 2). Based on the detection of parasite DNA in peripheral blood by varATS PCR assay, four main parasite clearance profiles could be identified (Fig. 3): Group A- individuals who cleared infection by day 1, (ii) Group B-those who remained positive until day 3, Group C-individuals who remained positive until day 7, and Group D-those who tested positive beyond day 7. varATS results showed that 12 (18.5%) of participants presented with parasitaemia 3 days post-treatment while eight (12.3%) patients presented with genotypable day 28 parasitaemia (Fig. 4).

Fig. 2
figure 2

Parasite count/μl by microscopy before and after artemether/lumefantrine treatment. (After day 1, the mean concentration was the same up to day 14 with re-appearing infections on days 21 and 28)

Fig. 3
figure 3

Kaplan-Meier estimates of parasite reduction following artemether-lumefantrine treatment (Group A represents % reduction of parasitaemia on day 1, (ii) Group B represents % reduction of parasitaemia on day 3, Group C represents % reduction of parasitaemia on day 7, and Group D represents % parasitaemia beyond day 7

Fig. 4
figure 4

Log10 Parasite count/μl by varATS qPCR before and after artemether/lumefantrine treatment (Low-grade residual infections persist on days 1, 3, 7, 4 and 28)

Relationship between enrolment parasite density and post-treatment parasitaemia

The relationship between parasite persistence after treatment and the levels of parasitaemia at enrolment was examined by varATS qPCR using Spearman’s correlation coefficient (rs). There was no linear association (P = 0.79) between the starting parasite density and day 3 parasitaemia. Similarly, there was no correlation (P = 0.69) between the enrolment parasitaemia and day 28 parasite count. There was a weak linear association between parasite densities on day 3 and reappearance on day 28 (rs = 0.12) even though the relationship was not statistically significant (P = 0.33). Analysis of parasitaemia and patient characteristics (Table 2) did not indicate any relationship between age of the participant and parasite density on day 0 (P = 0.25). There was a weak inverse association between the age of patients and (i) persistence of parasitaemia on day 3 (rs = − 0.24; P = 0.056) and (ii) re-appearance of parasitaemia on day 28 (rs = − 0.41; P = 0.0008). After adjusting for multiple testing using Bonferroni method and with a new cut-off of 0.0167 for the tests between age and days 0, 3 and 28 parasitaemia, we observed an inverse relationship between parasitaemia on day 28 and age of the patient. None of the patients with positive parasitaemia post-treatment reported with clinical symptoms of malaria.

Table 2 Mean (and standard deviation) of parasitaemia (by varATS) per age group

Genetic description of day 0 and day 28 parasite populations

Allele frequencies at each of the barcodes in the day 0 and day 28 populations are presented in Additional file 2: Datasheet S1. The complexity of infection in the day 0 population was 1.30 while the CoI on day 28 marginally increased to 1.34 on day 28 (Additional file 2: Datasheet S2), however, this difference was not statistically significant (P = 0.42). The mean HE values across all barcodes on day 0 was 0.50 ± 0.05 and 0.56 ± 0.05 on day 28. There was no evidence of significant linkage disequilibrium (Additional file 1: Table S4), although the \( {I}_A^S \) value was higher in day 0 (0.08) than in day 28 populations (0.03).

Signatures of relatedness between individual pre- and post-treatment infections

We used the analysis of barcode sequence identity to estimate relatedness between pairwise parasite isolates in day 0 and day 28 infections (Additional file 2: Datasheet S3). The barcode (π) pairwise comparisons showed high genetic relatedness of day 0 and re-appearing day 28 parasite isolates in three (37.5%) of the eight participants who presented with day 28 parasitaemia (Table 3).

Table 3 Age of participants and barcode similarity (π) indices of parasites before and after drug treatment

Prevalence of drug resistance polymorphisms

When the pre-treatment parasites were taken together, the crt-76 T mutant allele was present in 38 (58.5%) isolates. The mdr1–86 mutant allele was found in 56 (86.2%) of isolates. We observed 94.4% mutant and 5.6% wildtype alleles of Dhps-540 and 100% mutant alleles of Dhfr-108 in the population. No mutant allele of the K-13 C580Y was observed (Fig. 5). Pairwise phenotype-genotype comparisons showed no association of mdr1 (N86Y and Y184F) and crt K76 T mutations with persistence of parasitaemia on day 3 (Table 4) and re-appearance on day 28 (Table 5).

Fig. 5
figure 5

Prevalence of polymorphisms in Plasmodium falciparum drug resistance genes in pre-treatment parasite population (K13: Kelch 13 gene; Dhps: Dihydropteroate synthase; Pfcrt: Plasmodium falciparum chloroquine resistance transporter gene; Pfmdr1: Plasmodium falciparum chloroquine resistance transporter gene; Dhfr: Dihydrofolate reductase)

Table 4 Phenotype-genotype association of infections on day 3
Table 5 Pairwise comparison of drug resistance genotypes in pre-treatment and re-appearing P. falciparum infections

Discussion

One of the key approaches to malaria control is prompt treatment of patients for an improvement in health and prevention of onward transmission. With rising reports of drug-resistant Plasmodium falciparum, surveillance of ACT efficacy becomes important so that recrudescent infections can be promptly detected, treated and documented. In this study, we reported in vivo re-appearance of parasites following appropriate treatment of uncomplicated P. falciparum with artemether-lumefantrine in Nigeria. The genomic barcode analysis of pre-treatment and re-appearing infections suggested recrudescence despite the absence of mutation in K-13 C580Y. This suggests that the potency of the most promising antimalarial drug, artemisinin-based combination, in treating P. falciparum infection in Nigeria may be progressively waning.

High genetic diversity of the parasites was observed in the population. This observation is in line with expectations from high transmission settings where multiple parasite genomes in each sample and highly variable combinations of barcodes among individuals in the population have been observed [19]. We also observed a marginal increase in CoI 28 days after treatment suggesting a possible lack of clonal expansion following treatment. This agrees with the observation in Kenya where, contrary to expectations of reduced CoIs caused by the persistence of only drug-resistant strains in post-treatment infections, CoI on day 28 was higher than day 0 [20].

Age and immunity have previously been demonstrated as indicators of varying therapeutic responses [21]. Thus, older individuals living in a high transmission area who have a more developed acquired immunity than children in the same setting would be expected to have fewer persistent infections following treatment. While our finding could not directly correlate age with the persistence of parasitaemia on the third day of treatment, it revealed re-appearance of infections on day 28 only among children less than fifteen years old. This observation possibly explains the role of acquired immunity in complementing antimalarial drug action in parasite clearance. However, additional studies are required to understand the influence of immunity on therapeutic response and drug resistance by correlating parasite clearance rates with malaria antibody titres in the individuals.

Although we did not observe any clinical outcome of persistent infections and there was no direct association between re-appearance of infection and day 0 parasite densities, it is important to continue monitoring the prevalence of residual infections following treatment, even though they were observed in low quantities in this study. Since parasitaemia as low as 0.5/μl of blood corresponds to as high as approximately one million parasites in the body [21], low-grade infections observed after treatment in this study area may sustain a reservoir for silent transmission of the disease. Studies that involve longer follow-up durations and laboratory adaptation of the sub-microscopic parasites may provide more information on the transmission potential of these post-treatment parasites.

A high frequency of mutant alleles of crt and mdr1 (N86Y and Y134F) was reported in this study, but these mutations have not been directly associated with post- treatment re-appearance of infections. In addition, contrary to a previous report showing that polymorphisms in mdr1 genes were associated with decreased sensitivity to lumefantrine [22], our phenotype-genotype association study provided no evidence that the presence of 86Y allele sensitized parasites to lumefantrine as there was no correlation between the distribution of N86 and 134F alleles and presence of parasites on day 28. An increase in the frequency of Pfcrt 76 T is associated with reduced sensitivity to chloroquine [23] and the removal of drug pressure is expected to cause a decline in resistance conferring mutations [24]. The prevalence of mutant 76 T allele in this present investigation is lower than previously reported before the proscription of CQ for mild malaria treatment in Nigeria [25]. This may be a pointer to ongoing selection for chloroquine-susceptible parasites following widespread use of artemisinin. Kelch-13 sequencing information from recent therapeutic efficacy studies conducted in Cambodia, Lao and Vietnam has revealed cysteine to tyrosine mutation in codon 580 of K13 gene (C580Y) as the dominant polymorphism associated with artemisinin resistance [3]. However, our investigation revealed the absence of this mutation in re-appearing parasites.

We propose further investigations involving complete genome-wide associations between pre-treatment and persistent parasites to circumvent the limitation of our study which did not preclude the influence of re-infection with the same parasite type or residual parasite DNA and non-replicating mature gametocytes that can persist at low concentrations after treatment. Computational analyses that involve the combination of genomic information with epidemiological simulations will present complementary tools for interpreting population-level impact of drug interventions. In addition, since P. falciparum is thought to exhibit drug-induced developmental arrest until the activity of artemisinin in the host has diminished [26], details about quiescence in post-treatment infections using in vitro models may further explain the underlying mechanisms of re-appearance of parasites after treatment in the study area. A broader line of investigation on the efficacy of ACTs in larger populations in the country and across the borders is necessary to understand spatial and temporal heterogeneity of infectious reservoirs post-treatment.

Conclusions

We have reported in vivo re-appearance of parasites following appropriate treatment of uncomplicated Plasmodium falciparum with Artemether/Lumefantrine (AL) in Nigeria. The genomic barcode analysis of pre-treatment and re-appearing infections suggested recrudescence despite the absence of mutation in K-13 C580Y. This finding implies that the Nigerian parasites could be deploying a more complex drug evasion mechanism beyond the C580Y polymorphism. Further investigation of a larger population of parasites in the country is required for signs of reduced efficacy of artemisinin-based treatment of malaria.

Abbreviations

ACT:

Artemisinin-based Combination Treatment

AL:

Artemether/Lumefantrine

CoI:

Complexity of Infection

crt :

Chloroquine Resistance Transporter

dhfr :

Dihydrofolate Reductase

dhps :

Dihydropteroate Synthase

HE :

Expected Heterozygosity

HRM:

High Resolution Melting

K-13:

Kelch-13 Domain

MAF:

Minor Allele Frequency

mdr :

Multidrug Resistance

PCE:

Parasite Clearance Estimator

PCR:

Polymerase Chain Reaction

RDT:

Rapid Diagnostic Test

SNP:

Single Nucleotide Polymorphism

varATS:

var gene Acidic Terminal Sequence

References

  1. WHO, 2016: Summary. Geneva: World Health Organization; 2017 (WHO/HTM/GMP/2017.4).

  2. Adjuik M, Babiker A, Garner P, Olliaro P, Taylor W, White N et al. Artesunate combinations for treatment of malaria: meta-analysis. Lancet. 2004; 3:363(9402):9–17.

  3. Ashley EA, Dhorda M, Fairhurst RM, Amaratunga C, Lim P, Suon S, et al. Tracking resistance to artemisinin collaboration (TRAC). Spread of artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2014;371:411–23.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Lu F, Culleton R, Zhang M, Ramaprasad A, von Seidlein L, Zhou, et al. Emergence of indigenous artemisinin-resistant Plasmodium falciparum in Africa. N. Engl J Med. 2017;376(10):991–3.

    Article  PubMed  Google Scholar 

  5. Teklemariam M, Assefa A, Kassa M, Mohammed H, Mamo H. Therapeutic efficacy of artemether-lumefantrine against uncomplicated Plasmodium falciparum malaria in a high-transmission area in Northwest Ethiopia. PLoS One. 2017;12(4):e0176004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. White NJ. Malaria parasite clearance. Malar J. 2017;16:88.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Hofmann N, Mwingira F, Shekalaghe S, Robinson LJ, Mueller I, Felger I. Ultra-Sensitive Detection of Plasmodium falciparum by Amplification of Multi-Copy Subtelomeric Targets. PLoS Med. 2015; https://doi.org/10.1371/journal.pmed.1001788.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Thompson J, Fernandez-Reyes D, Sharling L, Moore SG, Eling WM, Kyes SA, et al. Plasmodium cysteine repeat modular proteins 1-4: complex proteins with roles throughout the malaria parasite life cycle. Cell Microbiol. 2017;9:1466–80.

    Article  CAS  Google Scholar 

  9. Daniels R, Volkman SK, Milner DA, Mahesh N, Neafsey DE, Park DJ, et al. A general SNP-based molecular barcode for Plasmodium falciparum identification and tracking. Malar J. 2008;7:223.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. MalariaGEN Plasmodium falciparum Community Project (2016) eLife, 5: e08714.

  11. Daniels R, Chang HH, Séne PD, Park DC, Neafsey DE, Schaffner SF, et al. Genetic surveillance detects both clonal and epidemic transmission of malaria following enhanced intervention in Senegal. PLoS One. 2013;8(4):e60780.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  12. WHO. Methods for surveillance of anti-malarial drug efficacy. Geneva, Switzerland: World Health Organization; 2009.

    Google Scholar 

  13. Federal Ministry of Health. National Antimalarial treatment guidelines; 2005. p. 5–22.

    Google Scholar 

  14. Ojurongbe O, Adegbosin OO, Taiwo SS, Alli OA, Olowe OA, Ojurongbe TA, et al. Assessment of clinical diagnosis, microscopy, rapid diagnostic tests, and polymerase chain reaction in the diagnosis of Plasmodium falciparum in Nigeria. Malaria Research and Treatment. 2013;2013:1–5.

    Article  CAS  Google Scholar 

  15. Daniels RF, Schaffner SF, Wenger EA, Proctor JL, Chang HH, Wong W, et al. Modeling malaria genomics reveals transmission decline and rebound in Senegal. Proc Natl Acad Sci U S A. 2015;112:7067–72.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  16. Baniecki ML, Faust AL, Schaffner SF, Park DJ, Galinsky K, Daniels RF, et al. Development of a single nucleotide polymorphism barcode to genotype Plasmodium vivax infections. PLoS Negl Trop Dis. 2015;9(3):e0003539.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Peakall R, Smouse PE. GENALEX 6: genetic analysis in excel. Population genetic software for teaching and research. Mol Ecol. 2006;6:288–95.

    Article  Google Scholar 

  18. Galinsky K, Valim C, Salmier A, de Thoisy B, Musset L, Legrand E, et al. COIL: a methodology for evaluating malarial complexity of infection using likelihood from single nucleotide polymorphism data. Malar J. 2015;14:4.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Volkman SK, Neafsey DE, Schaffner SF, Park DJ, Wirth DF. Harnessing genomics and genome biology to understand malaria biology. Nat Rev Gen. 2012;13:315–28.

    Article  CAS  Google Scholar 

  20. Chang HH, Meibalan E, Zelin J, Daniels R, Eziefula AC, Meyer EC, et al. Persistence of Plasmodium falciparum parasitemia after artemisinin combination therapy: evidence from a randomized trial in Uganda. Sci Rep. 2016;6:26330.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  21. White NJ, Pukrittayakamee S, Hien TT, Faiz MA, Mokuolu OA, Dondorp AM. Malaria. Lancet. 2013;383:723–35.

    Article  PubMed  Google Scholar 

  22. Sisowath C, Strömberg J, Mårtensson A, Msellem M, Obondo C, Björkman A, et al. In vivo selection of Plasmodium falciparum pfmdr1 86N coding alleles by Artemether-Lumefantrine (Coartem). J Infect Dis. 2005;191(6):1014–7.

    Article  CAS  PubMed  Google Scholar 

  23. Djimdé A, Doumbo OK, Cortese JF, Kayentao K, Doumbo S, Diourté Y, et al. A molecular marker for chloroquine-resistant falciparum malaria. N Engl J Med. 2001;344:257–63.

    Article  PubMed  Google Scholar 

  24. Kublin JG, Cortese JF, Njunju EM, Mukadam RA, Wirima JJ, Kazembe PN, et al. Reemergence of chloroquine-sensitive Plasmodium falciparum malaria after cessation of chloroquine use in Malawi. J Infect Dis. 2003;187:1870–5.

    Article  PubMed  Google Scholar 

  25. Olukosi YA, Oyebola KM, Ajibaye O, Orok BA, Aina OO, Agomo CO, et al. Persistence of markers of chloroquine resistance among P falciparum isolates recovered from two Nigerian communities. Mal World J. 2014;5:3.

    Google Scholar 

  26. Witkowski B, Lelièvre J, Barragán MJ, Laurent V, Su XZ, Berry A, et al. Increased tolerance to artemisinin in Plasmodium falciparum is mediated by a quiescence mechanism. Antimicr Ag Chem. 2010;54(5):1872–7.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The authors appreciate the co-operation of participants and the efforts of the staff and management of the health facilities visited. Special recognition goes to Dr. S. Babatimilehin for coordinating the patients on follow-up. We thank the Department of Immunology and Infectious Diseases, Harvard TH Chan School of Public Health, Boston, MA, USA for the training on Molecular SNP Barcoding and HRM analysis. The authors acknowledge Dr. Rachel Daniels for her comments on the manuscript.

Funding

Dr. Kolapo Muyiwa Oyebola was supported by a DELTAS Africa grant (DEL-15-007: Awandare). The DELTAS Africa Initiative is an independent funding scheme of the African Academy of Sciences (AAS)‘s Alliance for Accelerating Excellence in Science in Africa (AESA) and supported by the New Partnership for Africa’s Development Planning and Coordinating Agency (NEPAD Agency) with funding from the Wellcome Trust (107,755/Z/15/ZZ: Awandare) and the UK government. The views expressed in this publication are those of the authors and not necessarily those of AAS, NEPAD Agency, Wellcome Trust or the UK government.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article and its additional files.

Author information

Authors and Affiliations

Authors

Contributions

Conceived and designed the study: KMO, GAA and AA-N. Implemented the field study: KMO, OOA, OSA and YAO. Carried out the laboratory studies: KMO and AA-N. Collated and analysed the data: KMO and AA-N. Drafted the manuscript: KMO, ETI, OAO, AA-N and GAA. Edited the manuscript: OOA, TSA, YAO, ETI, GAA, and AA-N. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Kolapo Muyiwa Oyebola.

Ethics declarations

Ethical approval and consent to participate

Individuals were requested to provide a written informed consent or assent (by parents/guardians of children below 16 years) before progressing on any study-related procedure. Ethical approval for the study was obtained from the Nigerian Institute of Medical Research Ethics Committee (IRB/17/019) and the Management Board of the health facilities.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Table S1. Primer sequences and qPCR conditions for varATS assay. Table S2. High Resolution Melting Drug Resistance Assay. Table 3 Relationship between day 0 parasite density and percentage of individuals with day 3 parasitaemia. Table S4. Linkage disequilibrium analysis for pre- and post-treatment parasite populations. (DOCX 26 kb)

Additional file 2:

Datasheet S1. Allele frequencies at each of the barcodes in the day 0 and day 28 populations. Datasheet S2. Complexity of infection in day 0 and day 28 populations. Datasheet S3. Pairwise sequence comparison of day 0 and day 28 parasite isolates. (ZIP 94 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Oyebola, K.M., Aina, O.O., Idowu, E.T. et al. A barcode of multilocus nuclear DNA identifies genetic relatedness in pre- and post-Artemether/Lumefantrine treated Plasmodium falciparum in Nigeria. BMC Infect Dis 18, 392 (2018). https://doi.org/10.1186/s12879-018-3314-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-018-3314-3

Keywords