Skip to main content

Clearance dynamics of lactate dehydrogenase and aldolase following antimalarial treatment for Plasmodium falciparum infection

Abstract

Background

Lingering post-treatment parasite antigen in blood complicates malaria diagnosis through antigen detection. Characterization of antigen clearance dynamics is important for interpretation of positive antigen detection tests.

Results

We used a bead-based serological assay to measure lactate dehydrogenase (LDH), aldolase (Aldo), and histidine-rich protein 2 (HRP2) levels in 196 children with Plasmodium falciparum malaria treated with effective antimalarials and followed for 28 to 42 days as part of therapeutic efficacy studies in Angola. Compared to pre-treatment levels, antigen concentrations two days after treatment declined by 99.7% for LDH, 96.3% for Aldo, and 54.6% for HRP2. After Day 2, assuming a first-order kinetics clearance model, half-lives of the antigens were 1.8 days (95% CI: 1.5–2.3) for LDH, 3.2 days (95% CI: 3.0–3.4) for Aldo, and 4.8 days (95% CI: 4.7–4.9) for HRP2.

Conclusions

LDH and Aldo show substantially different clearance rates than HRP2, and their presence is largely indicative of active infection.

Background

Malaria diagnosis worldwide often relies on detection of parasite antigen using rapid diagnostic tests (RDTs) [1]. Most commonly, the target is the Plasmodium falciparum-specific histidine-rich protein 2 (HRP2), but other targets include lactate dehydrogenase (LDH) and aldolase (Aldo) [1, 2]. However, unlike a case confirmed through malaria microscopy, a positive RDT result does not necessarily imply an active infection because it may merely represent lingering antigenemia after successful clearance of parasites from the host.

As a consequence, the presence of parasite antigens in a person’s blood (leading to a possible RDT positive result) must be interpreted in the context of the clearance dynamics of these antigens. These dynamics have been well-described for HRP2, whose clearance obeys a first order kinetics model, with a half-life in the human host in the 3–5 day range [3, 4]. As a result, individuals can be HRP2-positive for weeks following adequate clearance of parasitemia [5, 6]. HRP2 clearance dynamics are highly consistent among different groups and parasite genotypes; deviations may indicate inadequate clearance of parasitemia in drug trials [3, 7]. Previous studies investigating the persistence of other P. falciparum antigens have found positive LDH-based RDT results for up to two weeks following successful parasite clearance [8, 9], though no LDH positive blood samples were identified 7 days post-treatment in an ELISA-based study [10]. To date, the longitudinal dynamics of Plasmodium Aldo in the human host has not been investigated. Information about clearance kinetics for these targets could inform interpretation of laboratory test results, particularly in persons who have recently received treatment. Here, we use a recently developed bead-based assay for laboratory antigen detection and quantification to characterize clearance dynamics of LDH and Aldo in individuals treated for Plasmodium falciparum infection.

Methods

Sample collection

We analyzed samples previously collected as part of therapeutic efficacy studies in Angola in 2015 [11] and 2017 [12]. Children with uncomplicated P. falciparum infection were treated with one of three artemisinin-based combination therapies: artemether–lumefantrine (AL), artesunate–amodiaquine (ASAQ), or dihydroartemisinin–piperaquine (DP), and followed for 28 days (AL and ASAQ) or 42 days (DP). Blood was collected on Whatman 903 filter paper (GE Healthcare, Chicago, IL, USA) on Day 0, Day 2 (2017 only), Day 3, Day 7, and weekly thereafter. Three-day treatment was given starting on Day 0. Parasite density upon enrollment in patients ranged between 2175–184,465 parasites/µl.

A total of 1500 samples from 196 subjects were included in this analysis, including 174 cases of adequate clinical and parasitological response and 22 microsatellite-confirmed recrudescences [11, 12]. Artemisinins remain efficacious in Angola and 97% of the analyzed participants had cleared microscopically-detectable parasitemia by Day 2, increasing to > 99% by Day 3.

Laboratory analysis

A single 6 mm dried blood spot (DBS) punch (corresponding to 10 µl whole blood) was taken from each timepoint for a person’s series. Whole blood was eluted into buffer and analyzed for LDH, Aldo, and HRP2 concentration using previously-described methods [13, 14]. In brief, the 6 mm punch was placed into 200 µl blocking buffer overnight to allow elution [Buffer B: 0.5% polyvinyl alcohol (Sigma-Aldrich, St. Louis, MO, USA; P8136) 0.5% polyvinylpyrrolidine (Sigma-Aldrich; PVP360), 0.1% casein (Thermo Fisher Scientific, Waltham, MA, USA; 37528), 0.5% BSA (Sigma-Aldrich; A9418), 0.3% Tween-20, 0.05% sodium azide, and 0.01% E. coli extract to prevent non-specific binding]. For antigen detection, three unique bead regions (Bio-Plex COOH bead, BioRad, Hercules, CA, USA; 171506XXX) were individually coated by the EDC/Sulfo-NHS intermediate reaction with separate antibodies specific for each antigen to be captured: Plasmodium Aldo (12.5 µl/12.5 × 106 beads, rabbit IgG anti-Aldo, Abcam, Cambridge, UK; ab207494), Plasmodium LDH (12.5 µl/12.5 × 106 beads, mouse IgG anti-LDH, BBI Solutions, Cardiff, UK; BM355-Z8F7), and P. falciparum PfHRP2 (20 µl/12.5 × 106 beads, mouse IgG anti-HRP2, Abcam; ab9206). For the multiplex antigen assay, the three unique bead regions were mixed together in 5 ml Buffer A (PBS, 0.5% BSA, 0.05% Tween20, 0.02% NaN3) so that 1500 of each bead region would be added per well in the assay plate. The beads were washed twice with 100 µL wash buffer (PBS, 0.05% Tween20) in the assay plate, and were incubated with 50 µl sample (corresponding to 2.5 µl whole blood) run in singlets in filter bottom plates (Millipore, Burlington, MA, USA; MABVN1250). Following a 90 min incubation with samples under gentle shaking at room temperature protected from light, plates were subsequently washed three times. Beads were incubated for 45 min with a 50 µl mix of detection antibodies: anti-pAldo (1:2000×, rabbit anti-Aldo, Abcam; ab207494), anti-pLDH [1:500× of 2:1:1 mixture (BBI Solutions BM355-P4A2:BioRad Pv-pLDH HCA156:BioRad Pf-pLDH HCA158), and anti-HRP2 (1:500×, mouse IgG anti-HRP2, Abcam, ab9203)]. All detection antibodies were previously biotinylated by Thermo Fisher Scientific EZ-Link Micro Sulfo-NHS-Biotinylation Kit according to the manufacturer’s protocol. Plates were again washed three times, and after a final 30 min wash step with reagent diluent, beads were washed once and resuspended in 100 µl PBS. After one minute on the shaker, plates were read on a Bio-Plex 200 instrument (BioRad, Hercules, CA, USA) by generating the median fluorescence intensity (MFI) signal for 50 beads in each unique region and then the mean fluorescence intensity of the MFIs among duplicates. The final measure, denoted as MFI-background (bg), was reported by subtracting a background signal defined as MFI values from beads on each plate only exposed to sample diluent during the sample incubation step.

To translate between a MFI-bg value and antigen concentration, standard curves for this relationship were calculated using recombinant antigens [14]. Recombinant LDH and HRP2 antigens were provided by Microcoat Biotechnologie GmbH (Bernried, Germany), and lyophilized preparations were rehydrated according to the manufacturer’s instructions. The Plasmodium vivax-specific isoform of Aldo was produced at the Centers for Disease Control and Prevention (CDC). Concentrations above 26,621 pg/ml for LDH, 464 pg/ml for Aldo, and 84 pg/ml for HRP2 were reported as positive.

Statistical analysis

The log-transformed concentrations of each antigen were plotted over time. For each antigen, Kaplan–Meier curves were fit to estimate the time to clearance to below the level of detection (LOD) of the bead-based assay, stratifying by treatment outcome. Differences in antigen persistence by treatment outcome were compared using the log-rank test.

A Bayesian Markov Chain Monte Carlo (MCMC) algorithm was used to impute antigen concentrations for time points that were below the LOD of the bead-based assay for LDH (1111/1185 data points below the LOD), Aldo (723/1208), and HRP2 (7/945). In brief, the concentration at each time was modelled to follow a log-normal distribution, and a first-order kinetics model was fit to each individual’s antigen concentration time series data. After visual inspection showed an evident difference in clearance dynamics pre- and post-Day 2, separate pre-Day 2 (\(\lambda_{1}\)) and post-Day 2 (\(\lambda_{2}\)) clearance constants were estimated. The clearance constants, the standard deviation (\(\sigma\)), and the missing observations were sampled using an MCMC algorithm using Gibbs sampling. Uniform priors were used for \(\lambda_{1}\), \(\lambda_{2}\), and σ. The median and 95% credible intervals (CIs) were estimated from 9000 samples from the MCMC output, after discarding the first 1000 ‛burn-in’ samples. Data from patients with undetectable antigen at day 0 (n = 2/196 for LDH; n = 1/196 for Aldo; and n = 0/196 for HRP2) or time courses that were not monotonically declining (n = 24/196 for LDH; n = 43/196 for Aldo; and n = 73/196 for HRP2), defined as a > 25% relative or > 1000 pg/ml absolute increase in adjacent time points, were excluded from the estimation of clearance rates.

Statistical analysis was carried out in R version 3.3.2 (R Foundation for Statistical Computing, Vienna, Austria).

Results

Pre-treatment Day 0 concentrations of all antigens were strongly correlated with initial microscopy-determined parasite densities, with the relationship strongest for LDH (R2 = 0.316), followed by Aldo (R2 = 0.245) and HRP2 (R2 = 0.168) (Additional file 1: Figure S1). Day 0 concentrations of Aldo (median: 35,758 pg/ml; geometric mean: 30,374 pg/ml) were orders of magnitude lower than concentrations of LDH (median: 555,255 pg/ml; geometric mean: 681,374 pg/ml) or HRP2 (median: 699,256 pg/ml; geometric mean: 787,603 pg/ml) (Table 1). On Day 0, the median ratio of HRP2 to LDH concentration was 1.2 (IQR: 0.6–2.0), of HRP2 to Aldo was 26 (8–84), and of LDH to Aldo was 21 (12–34).

Table 1 Clearance parameters for post-treatment dynamics of LDH, Aldo, and HRP2 proteins from P. falciparum-infected children in Angola

While the average decline between Day 0 and Day 2 HRP2 concentration was only 54.6%, it was 99.7% for LDH and 96.3% for Aldo (Fig. 1, Table 1). Nearly all patients had detectable HRP2 over the entire duration of follow-up, but levels of LDH and Aldo positivity declined over time (Fig. 2, Additional file 1: Figure S2). For both LDH and Aldo, rates of antigen persistence were higher in individuals that eventually recrudesced compared to those who did not, although this did not reach statistical significance (log-rank test P-value 0.2 for both LDH and Aldo).

Fig. 1
figure 1

Decline in absolute malaria antigen following antimalarial treatment in Angolan children with uncomplicated P. falciparum infection. Dashed red line represents the limit of detection for each antigen

Fig. 2
figure 2

Persistence of malaria antigens following antimalarial treatment in Angolan children with uncomplicated P. falciparum infection. Dashed lines represent confidence intervals of the survival function

Based on first-order kinetics models fit to post-Day 2 antigen concentrations, long-term clearance of LDH (half-life: 1.8 days, 95% CI: 1.5–2.3) was quickest, followed by Aldo (3.2 days, 95% CI: 3.0–3.4) and then HRP2 (4.8 days, 95% CI: 4.7–4.9) (Table 1). When compared to clearance after Day 2, clearance before Day 2 was: 7.4 times quicker for LDH, 8.1 times quicker for Aldo, and 2.7 times quicker for HRP2.

Discussion

Fundamental differences in the clearance rates of LDH, Aldo, and HRP2 following microscopic clearance of P. falciparum infection were observed. On one extreme, nearly all LDH antigen (> 99%) is cleared within the first two days following antimalarial treatment, with the remaining amount being cleared at a much slower rate, with a half-life of 1.8 days. On the other extreme, more than half of the Day 0 HRP2 antigen is still present at Day 2 and HRP2 is cleared at an even slower rate after that (half-life of 4.8 days). These post-Day 2 clearance estimates for the HRP2 antigen are quite similar to previous estimates [3], confirming the long-term persistence of this antigen in the human host. Aldo takes an intermediate position in the spectrum between the other two antigens, and these differences hint at inherent biological differences in the clearance of these exogenous proteins from humans. LDH and Aldo are both expressed inside the parasite and are crucial for the functioning of the parasite cell [15, 16]. In contrast, HRP2 is largely secreted outside the parasite [17], and its biological function is not known. Its persistence is largely due to accumulated HRP2 in once-infected red blood cells [18]. In this context, it is not surprising that the vast majority of LDH and Aldo is cleared concurrent with clearance of the parasites themselves, and that the unbound extracellular HRP2 does not show such a dramatic decline in the first two days. Moreover, the ratios of pre- and post-day 2 clearance rates for LDH and Aldo are very similar, suggesting the two antigens might share a similar mechanism of clearance. The link between LDH and Aldo concentration and the presence of the parasite itself is consistent with the previous observation that LDH and Aldo levels are more closely associated with current parasite density, compared to HRP2, which can be interpreted more accurately as a measure of cumulative parasite load over the course of the infection [4, 19, 20].

A further nuance is introduced by apparent differences in antigen expression. The pre-treatment concentration of Aldo was 20–30 times lower than LDH and HRP2, which confirms previous results showing lower expression of Aldo compared to LDH and HRP2 [1]. These lower absolute levels of the protein during blood stage infection should be considered in the design of tests using Aldo as a diagnostic target.

Ultimately, antigen positivity in individuals living in endemic areas should be interpreted in the context of the differences in post-treatment clearance. Previous analysis of outpatients presenting to Angola health facilities showed that all patients with detectable levels of all three antigens were also positive for parasite nucleic acids, whereas PCR positivity rates were lower in individuals just positive for Aldo and HRP2 and much lower for HRP2 positivity alone [14]. Moreover, previous analytical work has shown that while presence of HRP2 is not necessarily synonymous with acute febrile disease, presence of LDH and Aldo at levels detectable through the bead-based assay used here are predictive of parasite densities high enough to provoke clinical disease [20].

These results are only from one country and are not necessarily generalizable. Rapid clearance of LDH and Aldo, combined with a higher assay LOD for these antigens, limit the ability of a definitive statement regarding the clearance dynamics of these antigens to be drawn from this study. In particular, it was not possible to test whether a second-order kinetics model provided a better fit than a first-order kinetics model. However, in relation to the slow clearance dynamics of the HRP2 antigen, these two pan-Plasmodium targets would clearly be more indicative of an active infection since they were absent in almost all patients within 7 days of starting an ACT regimen. Finally, the observation of a substantial number of patients showing post-day 7 spikes in antigen concentration could be due to false positive readings on the bead-based assay, or evidence of undetected new infections or recrudescences during the period of follow-up, and warrants further investigation.

Conclusions

Differential clearance of malaria antigens used as targets for diagnostic testing raises the possibility of inferring past history of the disease course in an individual based on the presence and quantity of the different antigens. HIV programs have for years capitalized on the differences in the dynamics of how certain antibodies appear in the host to distinguish recent infections from older infections, allowing estimation of incidence from cross-sectional survey data [21]. Analogously, specimens from household community surveys in malaria-endemic countries could be potentially analyzed by the multi-antigen bead-based assay to distinguish past infections from active ones, and possibly make inferences regarding treatment rates preceding the survey. Findings from this study should be confirmed in other settings and countries.

Availability of data and materials

All data are available from the authors upon request.

Abbreviations

AL:

artemether–lumefantrine

Aldo:

aldolase

ASAQ:

artesunate–amodiaquine

CDC:

Centers for Disease Control and Prevention

DP:

dihydroartemisinin–piperaquine

HRP2:

histidine-rich protein 2

LDH:

lactate dehydrogenase

LOD:

level of detection

MCMC:

Markov chain Monte Carlo

MFI-bg:

median fluorescence intensity-background

RDT:

rapid diagnostic test

References

  1. WHO. Malaria rapid diagnostic test performance: results of WHO product testing of malaria RDTs: round 7 (2015–2016). Geneva: World Health Organization; 2017.

    Google Scholar 

  2. Chiodini PL. Malaria diagnostics: now and the future. Parasitology. 2014;141:1873–9.

    Article  Google Scholar 

  3. Plucinski MM, Dimbu PR, Fortes F, Abdulla S, Ahmed S, Gutman J, et al. Post-treatment HRP2 clearance in patients with uncomplicated Plasmodium falciparum malaria. J Infect Dis. 2018;217:685–92.

    Article  CAS  Google Scholar 

  4. Dondorp AM, Desakorn V, Pongtavornpinyo W, Sahassananda D, Silamut K, Chotivanich K, et al. Estimation of the total parasite biomass in acute falciparum malaria from plasma PfHRP2. PLoS Med. 2005;2:e204.

    Article  Google Scholar 

  5. Bell DR, Wilson DW, Martin LB. False-positive results of a Plasmodium falciparum histidine-rich protein 2-detecting malaria rapid diagnostic test due to high sensitivity in a community with fluctuating low parasite density. Am J Trop Med Hyg. 2005;73:199–203.

    Article  CAS  Google Scholar 

  6. Tjitra E, Suprianto S, McBroom J, Currie BJ, Anstey NM. Persistent ICT malaria Pf/Pv panmalarial and HRP2 antigen reactivity after treatment of Plasmodium falciparum malaria is associated with gametocytemia and results in false-positive diagnoses of Plasmodium vivax in convalescence. J Clin Microbiol. 2001;39:1025–31.

    Article  CAS  Google Scholar 

  7. Koita OA, Krogstad DJ. Converting a liability to an asset: using the clearance of a malaria parasite protein from the blood of infected subjects to predict the outcome of treatment. J Infect Dis. 2018;217:683–4.

    Article  Google Scholar 

  8. Iqbal J, Siddique A, Jameel M, Hira PR. Persistent histidine-rich protein 2, parasite lactate dehydrogenase, and panmalarial antigen reactivity after clearance of Plasmodium falciparum monoinfection. J Clin Microbiol. 2004;42:4237–41.

    Article  CAS  Google Scholar 

  9. Aydin-Schmidt B, Mubi M, Morris U, Petzold M, Ngasala BE, Premji Z, et al. Usefulness of Plasmodium falciparum-specific rapid diagnostic tests for assessment of parasite clearance and detection of recurrent infections after artemisinin-based combination therapy. Malar J. 2013;12:349.

    Article  Google Scholar 

  10. Bashir IM, Otsyula N, Awinda G, Spring M, Schneider P, Waitumbi JN. Comparison of PfHRP-2/pLDH ELISA, qPCR and microscopy for the detection of Plasmodium events and prediction of sick visits during a malaria vaccine study. PLoS ONE. 2013;8:e56828.

    Article  CAS  Google Scholar 

  11. Plucinski MM, Dimbu PR, Macaia AP, Ferreira CM, Samutondo C, Quivinja J, et al. Efficacy of artemether–lumefantrine, artesunate–amodiaquine, and dihydroartemisinin–piperaquine for treatment of uncomplicated Plasmodium falciparum malaria in Angola, 2015. Malar J. 2017;16:62.

    Article  Google Scholar 

  12. Davlantes E, Dimbu PR, Ferreira CM, Joao MF, Pode D, Félix J, et al. Efficacy and safety of artemether–lumefantrine, artesunate–amodiaquine, and dihydroartemisinin–piperaquine for the treatment of uncomplicated Plasmodium falciparum malaria in three provinces in Angola, 2017. Malar J. 2018;17:144.

    Article  Google Scholar 

  13. Rogier E, Plucinski M, Lucchi N, Mace K, Chang M, Lemoine JF, et al. Bead-based immunoassay allows sub-picogram detection of histidine-rich protein 2 from Plasmodium falciparum and estimates reliability of malaria rapid diagnostic tests. PLoS ONE. 2017;12:e0172139.

    Article  Google Scholar 

  14. Plucinski MM, Herman C, Jones S, Dimbu R, Fortes F, Ljolje D, et al. Screening for Pfhrp2/3-deleted Plasmodium falciparum, non-falciparum, and low-density malaria infections by a multiplex antigen assay. J Infect Dis. 2018;219:437–47.

    Article  Google Scholar 

  15. Wanidworanun C, Nagel RL, Shear HL. Antisense oligonucleotides targeting malarial aldolase inhibit the asexual erythrocytic stages of Plasmodium falciparum. Mol Biochem Parasitol. 1999;102:91–101.

    Article  CAS  Google Scholar 

  16. Deck LM, Royer RE, Chamblee BB, Hernandez VM, Malone RR, Torres JE, et al. Selective inhibitors of human lactate dehydrogenases and lactate dehydrogenase from the malarial parasite Plasmodium falciparum. J Med Chem. 1998;41:3879–87.

    Article  CAS  Google Scholar 

  17. Howard RJ, Uni S, Aikawa M, Aley SB, Leech JH, Lew AM, et al. Secretion of a malarial histidine-rich protein (Pf HRP II) from Plasmodium falciparum-infected erythrocytes. J Cell Biol. 1986;103:1269–77.

    Article  CAS  Google Scholar 

  18. Poti KE, Balaban AE, Pal P, Kobayashi T, Goldberg DE, Sinnis P, et al. In vivo compartmental kinetics of Plasmodium falciparum histidine-rich protein II in the blood of humans and in BALB/c mice infected with a transgenic Plasmodium berghei parasite expressing histidine-rich protein II. Malar J. 2019;18:78.

    Article  Google Scholar 

  19. Ochola L, Marsh K, Lowe B, Gal S, Pluschke G, Smith T. Estimation of the sequestered parasite load in severe malaria patients using both host and parasite markers. Parasitology. 2005;131:449–58.

    Article  CAS  Google Scholar 

  20. Plucinski MM, Rogier E, Dimbu R, Fortes F, Halsey ES, Aidoo M, et al. Performance of antigen concentration thresholds for attributing fever to malaria among outpatients in Angola. J Clin Microbiol. 2019;57:01901–18.

    Article  Google Scholar 

  21. Janssen RS, Satten GA, Stramer SL, Rawal BD, Oʼbrien TR, Weiblen BJ, et al. New testing strategy to detect early HIV-1 infection for use in incidence estimates and for clinical and prevention purposes. JAMA. 1998;280:42–8.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention.

Funding

MMP and ESH were supported by the U.S. President’s Malaria Initiative.

Author information

Authors and Affiliations

Authors

Contributions

MMP and ER designed the study. MMP, PRD, and FF collected the samples. DN and ER performed laboratory analyses. MMP, PDM, PRD, FF, ESH, and ER analyzed the data. MMP and ER wrote the first draft. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Mateusz M. Plucinski.

Ethics declarations

Ethics approval and consent to participate

The original therapeutic efficacy studies were approved by human subjects review boards at the Centers for Disease Control and Prevention and the Angolan Ministry of Health. Additional antigen detection in the previously-collected samples was approved by the office of the Associate Director of Science in the Center for Global Health at the Centers for Disease Control and Prevention.

Consent for publication

All authors read the final manuscript and provided consent for publication.

Competing interests

The authors declare they do not have any commercial or other associations that might pose a competing interests.

Additional information

Publisher's Note

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

Additional file

Additional file 1: Figure S1.

Relationship between pre-treatment (Day 0) microscopically-determined parasite density and LDH, Aldo, and HRP2 concentrations in Angolan children with uncomplicated P. falciparum infection. Figure S2. Decline in absolute malaria antigen following antimalarial treatment in Angolan children with uncomplicated P. falciparum infection. Dashed red line represents the limit of detection for each antigen. Each line represents a time course for a single participant.

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

Plucinski, M.M., McElroy, P.D., Dimbu, P.R. et al. Clearance dynamics of lactate dehydrogenase and aldolase following antimalarial treatment for Plasmodium falciparum infection. Parasites Vectors 12, 293 (2019). https://doi.org/10.1186/s13071-019-3549-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-019-3549-x

Keywords