http://orcid.org/0000-0002-7822-7059
Figures
Abstract
Background
In Vietnam, the importance of vivax malaria relative to falciparum during the past decade has steadily increased to 50%. This, together with the spread of multidrug-resistant Plasmodium falciparum, is a major challenge for malaria elimination. A 2-year prospective cohort study to assess P. vivax morbidity after radical cure treatment and related risk factors was conducted in Central Vietnam.
Methods and findings
The study was implemented between April 2009 and December 2011 in four neighboring villages in a remote forested area of Quang Nam province. P. vivax-infected patients were treated radically with chloroquine (CQ; 25 mg/kg over 3 days) and primaquine (PQ; 0.5 mg/kg/day for 10 days) and visited monthly (malaria symptoms and blood sampling) for up to 2 years. Time to first vivax recurrence was estimated by Kaplan–Meier survival analysis, and risk factors for first and recurrent infections were identified by Cox regression models. Among the 260 P. vivax patients (61% males [159/260]; age range 3–60) recruited, 240 completed the 10-day treatment, 223 entered the second month of follow-up, and 219 were followed for at least 12 months. Most individuals (76.78%, 171/223) had recurrent vivax infections identified by molecular methods (polymerase chain reaction [PCR]); in about half of them (55.61%, 124/223), infection was detected by microscopy, and 84 individuals (37.67%) had symptomatic recurrences. Median time to first recurrence by PCR was 118 days (IQR 59–208). The estimated probability of remaining free of recurrence by month 24 was 20.40% (95% CI [14.42; 27.13]) by PCR, 42.52% (95% CI [35.41; 49.44]) by microscopy, and 60.69% (95% CI [53.51; 67.11]) for symptomatic recurrences. The main risk factor for recurrence (first or recurrent) was prior P. falciparum infection. The main limitations of this study are the age of the results and the absence of a comparator arm, which does not allow estimating the proportion of vivax relapses among recurrent infections.
Conclusion
A substantial number of P. vivax recurrences, mainly submicroscopic (SM) and asymptomatic, were observed after high-dose PQ treatment (5.0 mg/kg). Prior P. falciparum infection was an important risk factor for all types of vivax recurrences. Malaria elimination efforts need to address this largely undetected P. vivax transmission by simultaneously tackling the reservoir of P. falciparum and P. vivax infections.
Author summary
Why was this study done?
- WHO currently recommends a 14-day primaquine (PQ) regimen at 0.5 mg/kg/day together with chloroquine (CQ; 25 mg/kg over 3 days) for the radical cure of P. vivax infections in Southeast Asia.
- Between 2007 and 2009, Vietnam introduced a high-dose PQ regimen shortened to 10 days in an attempt to improve patient compliance.
- A 2-year one-arm cohort study was conducted between May 2009 and December 2011 in four remote villages of Central Vietnam to determine the incidence and risk factors for vivax recurrence after the new radical cure regimen.
What did the researchers do and find?
- 260 P. vivax-infected patients were enrolled between May 2009 and December 2010, received directly observed 10-day treatment, and were followed up monthly for up to 24 months.
- During scheduled and unscheduled visits, participants were examined, and blood samples were systematically taken for light microscopy (LM) and polymerase chain reaction (PCR) analysis. All LM-detected vivax infections (patent) were treated with CQ alone, while PQ was given at the end of individual follow-up.
- Time-to-event analysis showed that by month 24, only 20% of patients remained free of vivax recurrences detected by PCR, while this was 40% for patent (by LM) recurrences. The majority of recurrences were submicroscopic and asymptomatic because only 38% of patients had symptomatic recurrences.
- The main risk factors for P. vivax recurrences were prior P. falciparum infections, especially submicroscopic, and prior treatment of vivax infections with CQ.
What do these findings mean?
- Despite the age (10 years) of the data collected, the results of this cohort study remain relevant to the current malaria situation in Vietnam, where the elimination of vivax malaria and of multidrug-resistant P. falciparum are major challenges.
- A substantial proportion of the vivax malaria reservoir is not detected by current surveillance strategies in Vietnam, i.e., passive and re-active case detection based on LM/rapid diagnostic tests (RDTs).
- Prior P. falciparum infections (including submicroscopic) and vivax infections treated with CQ alone were significant risk factors for vivax recurrences.
- To further improve vivax elimination efforts in Vietnam, new treatment strategies (i.e., single-dose radical treatment) and more sensitive diagnostic tools (i.e., field PCR/ultrasensitive RDTs) simultaneously tackling falciparum and vivax malaria infections are needed.
Citation: Pham TV, Nguyen HV, Aguirre AR, Nguyen VV, A. Cleves M, Nguyen XX, et al. (2019) Plasmodium vivax morbidity after radical cure: A cohort study in Central Vietnam. PLoS Med 16(5): e1002784. https://doi.org/10.1371/journal.pmed.1002784
Academic Editor: Elizabeth A. Ashley, Mahidol University, THAILAND
Received: May 31, 2018; Accepted: March 19, 2019; Published: May 17, 2019
Copyright: © 2019 Pham et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: The authors received no specific funding for this work.
Competing interests: The authors have declared that no competing interests exist.
Abbreviations: ACT, artemisinin-based combination therapy; CGT, Conditional Gap Time; CHC, Commune Health Center; CQ, chloroquine; DHA-PQ, dihydroartemisinin-piperaquine; FPBS, filter paper blood sample; G6PDd, Glucose-6-Phosphate-dehydronage deficient; Hb, Hemoglobin; HHW, Hamlet Health Worker; HR, Hazard Ratio; ITMA, Institute of Tropical Medicine, Antwerp; LCF, late clinical failure; LLIN, long-lasting insecticidal net; LM, light microscopy; LPF, late parasitological failure; NIMPE, National Institute of Malariology, Parasitology, and Entomology; NMCEP, National Malaria Control and Elimination Program; PCD, passive case detection; PCR, polymerase chain reaction; PH, Proportional Hazard; PQ, primaquine; RDT, rapid diagnostic test; Re-ACD, re-active case detection; SM, submicroscopic; SnM, Seminested Multiplex; TB, tuberculosis; WBC, white blood cell; WHO, World Health Organization
Introduction
Radical cure of P. vivax malaria, the most prevalent species outside sub-Saharan Africa [1], remains a challenge because preventing subsequent relapses and onward transmission requires treatment of both circulating blood (schizonticidal drugs) and the dormant liver (hypnozoiticidal drug) stages [2]. Chloroquine (CQ) and primaquine (PQ) have been combined for this purpose since the late 1950s, and PQ—an 8-aminoquinoline licensed in 1952—remained the only antirelapse therapy available until July 2018 [3]. Its use implies several challenges, including the risk of acute hemolysis in Glucose-6-Phosphate-dehydrogrenase–deficient (G6PDd) patients and poor adherence to the 14-day treatment schedule [2]. Moreover, the rapid spread of P. vivax CQ resistance will require alternative blood schizonticides and thus new combinations for radical cure to be tested [4]. Over the past 60 years, various PQ regimens have been used [5], and currently the World Health Organization (WHO) guidelines recommend PQ at 0.25 mg/kg/day for 14 days, together with CQ or an artemisinin-based combination therapy (ACT) [6]. In Southeast Asia and Oceania, where the vivax Chesson strain is prevalent, a higher PQ daily dose (0.5 mg/kg/day; 7.0 mg/kg total dose) is recommended.
In Vietnam, the importance of P. vivax malaria has recently increased following the successful control efforts against P. falciparum [7] and the recent country’s commitment to eliminate malaria by 2030 [8]. Currently, malaria transmission remains confined in forested areas of Central Vietnam, with the highest incidence reported along the international borders with Laos and Cambodia [9]. Between 2011 and 2015, the annual incidence of confirmed malaria cases dropped by 44% (i.e., from 16,612 to 9,331) and malaria deaths by 79% (3 deaths recorded in 2015). Concomitantly, the annual incidence of confirmed P. falciparum cases decreased by 57% (from 10,101 to 4,327 cases) while P. vivax cases decreased only by 15% (from 5,602 to 4,756 cases). Consequently, the frequency of P. vivax relative to P. falciparum increased from 34% to 51%, a significant change because P. vivax consistently accounted for 20%–30% of all malaria cases since 1991 [10].
In the framework of the current National Malaria Control and Elimination Program (NMCEP) launched in 2011, surveillance strategies include both passive case detection (PCD) using light microscopy (LM) or rapid diagnostic tests (RDTs) (at the health centers, hospitals, and private clinics) and re-active case detection (Re-ACD) following the detection of any malaria case (“index case”) by PCD at community level [11]. Re-ACD consists of screening all residents within 20–30 houses around the index case for malaria (using RDTs/LM) and treating all positive individuals. Re-ACD is done by a mobile team (including health staff and Hamlet Health Workers [HHWs]) as part of the “2-4-7 surveillance and response strategy” adapted from China [12]. All antimalarial treatments are provided free of charge by the NMCEP, and the current guidelines for the radical treatment of P. vivax recommend the use of PQ at 0.25 mg/kg/day for 14 days together with CQ for 3 days [13].
However, over the past 20 years, treatment guidelines for vivax malaria in Vietnam have changed several times, illustrating the difficult balance between maximizing the effectiveness of and compliance to PQ while minimizing the risk of acute hemolysis in G6PDd patients. Therefore, after the 5-day PQ regimen at 0.5 mg/kg/day introduced in 1997, guidelines were changed to 0.25 mg/kg/day for 10 days in 2003 and to 0.50 mg/kg/day for 10 days between 2007 and 2009, the latter of which was used in the present study [14]. However, in 2009, this high-dose regimen was again changed to PQ at 0.25 mg/kg/day for 14 days [15] following increased fear of hemolysis in the absence of readily available G6PD testing. In remote areas of Central Vietnam, where treatment monitoring is challenging, adherence to PQ treatment is reportedly low because of both fear by the health staff of acute hemolysis in G6PDd patients and poor patient compliance to the 14-day course [16]. These challenges related to PQ treatment are commonly reported in all vivax-endemic areas [17,18].
We report here the results of a cohort study in which P. vivax patients were treated radically with the previously recommended high-dose PQ regimen (0.50 mg/kg/day for 10 days) and followed up monthly for up to 2 years, with the aim of evaluating the post-treatment transmission dynamics of recurrent vivax infections.
Materials and methods
Data were collected following a prospective study protocol that did not include a detailed prospective statistical analysis plan (S1 Text).
Study design
This is a 2-year prospective cohort study carried out in Central Vietnam in which P. vivax-infected patients were followed up monthly after radical treatment. At each visit, a clinical examination and blood sampling for malaria parasite identification by LM and polymerase chain reaction (PCR) were performed. During the follow-up, vivax recurrent infections identified by LM and/or symptomatic infections were treated with CQ alone. This aimed at better determining the vivax parasite reservoir in real-life conditions because PQ was not usually administered to these patients at the time this project was implemented.
Study site and population
The study was implemented between April 2009 and December 2011 in four neighboring villages (Villages 1–3 in the Tra Leng commune and Village 4 in the Tra Don commune) in the mountainous and forested district of Nam Tra My in Quang Nam province. Villages were organized in several scattered clusters of about 4 to 45 households each (total population of 1,810 individuals according to the March 2009 census). A detailed description of the sociodemographic characteristics and malariometric indices of this study area has been published elsewhere [19]. Briefly, the study population belonged to the M’nong (Villages 1–3) and Cadong (Village 4) ethnic groups, living mainly off slash-and-burn agriculture and cinnamon plantations in forest fields. Malaria transmission is perennial with two peaks, May–June and October–November, with two main vector species: Anopheles dirus sensu stricto and A. minimus sensu stricto [20,21]. In April 2009, at baseline, malaria prevalence in the four study villages was 7.8% by LM and 23.6% by PCR, with 60% P. falciparum and the remaining P. vivax infections [19]. The G6PD genetic polymorphism (Vianchiang mutation) was estimated at less than 1.5% in both males and females, with no difference between ethnic groups.
Data collection
P. vivax-infected patients were first identified during the baseline survey (April–May 2009) [19] and then by PCD at the Tra Leng Community Health Center (CHC) (May 2009–Dec 2010) until the target sample size was reached (see below). Inclusion criteria were age ≥3 years (following contemporary guidelines [14]; changed to ≥6 months in the current guidelines [15]) and ≤60 years old, axillary temperature ≥37.5 °C and/or history of fever in the previous 48 hours, microscopically confirmed P. vivax monoinfection (asexual stage), permanent residency in the study area, and ability and willingness to participate in the study confirmed by a written informed consent (parents/guardians for patients <18 years). Patients were excluded if they had any danger signs or severe malaria, any concurrent infection or underlying chronic condition (e.g., tuberculosis [TB], HIV, epilepsy, etc.) requiring specific treatment, severe malnutrition, known allergy or intolerance to study drugs, prior PQ treatment within the past month, or if they were pregnant (as confirmed by rapid test) or breastfeeding. Patients were treated daily with CQ (25 mg/kg over 3 days) and PQ (0.5 mg/kg/day for 10 days) according to the contemporary national guidelines [15]. During the first 10 days of follow-up, patients were examined daily, and treatment was directly observed. Patients were then asked to attend the CHC at day 14, 21, and 28 or if ill between scheduled visits. Any patent vivax recurrence was considered as treatment failure; patients received rescue treatment (dihydroartemisinin-piperaquine [DHA-PQ]) and PQ (10 days), and follow-up was stopped. From day 28 onwards, monthly home visits were done either until completing 2-year follow-up or the end of study in December 2011. Between scheduled visits, patients were advised to consult their HHW or the study team at the CHC if unwell. At each visit (scheduled or unscheduled), they were systematically interviewed, body temperature was collected, and a finger prick blood sample was taken for LM (thick and thin film) and later molecular analysis (filter paper blood sample [FPBS]). Any LM-confirmed P. vivax infection identified during the monthly follow-up was treated with a 3-day course of CQ; PQ was given to all patients at the end of the follow-up. In addition, LM-confirmed P. falciparum infections (either mono- or mixed infections) were treated with a 3-day course DHA-PQ as per national guidelines. All patients were closely followed up by the study team for the detection of adverse events, which were managed as per study protocol.
Definitions
The main study outcomes were the incidence and risk factors for vivax recurrence as measured by time-to-event analysis. A recurrence was defined as a vivax infection detected after the initial post-treatment parasite clearance. Because the parasite reservoir in low-transmission settings is mostly represented by submicroscopic (SM) and asymptomatic infections [19,22], vivax recurrences were systematically detected by both PCR (Seminested Multiplex [SnM] PCR detecting densities below 1 parasite/μl [23,24]) and LM (detection limit of 10–50 parasites/μl for experienced technicians). Therefore, three outcome variables were defined as follows: i) all recurrences identified by PCR (= “all PCR-detected recurrences”); ii) all recurrences identified by PCR and LM (= “patent recurrences”); and iii) all patent recurrences with malaria symptoms (i.e., body temperature ≥37.5 °C and/or other symptoms; = “symptomatic recurrences”).
The following covariates were considered for the risk factor analysis: village, ethnic groups (M’nong, Ca Dong), gender, age, occupation (farmers, other), bed net in house (none, at least one), economic level (low, medium, high), year of recruitment (2009, 2010), season of recruitment (dry, rainy), parasite density at day 0 (parasite/μl), and “prior P. falciparum infection within 2 months before P. vivax recurrence” (i.e., all PCR-detected P. falciparum infections). In addition, we also explored whether being treated with an antimalarial (ACT or CQ) during the previous visit was associated with a higher risk of P. vivax recurrence. Among risk factors, the economic level was computed using three different variables for livestock ownership, i.e., number of i) buffaloes, ii) cows, and iii) pigs, and by using principal component analysis as described previously [19].
Laboratory procedures
Parasite species and density were determined by LM. Thin films were fixed with methanol for 15–30 seconds; both thin and thick smears were stained with 3% Giemsa solution for 45 minutes. P. vivax and P. falciparum asexual stages were counted against 200 white blood cells (WBCs), assuming a mean WBC count of 8,000/μl. Quality control was done at the National Institute of Malariology, Parasitology, and Entomology (NIMPE), Hanoi, Vietnam by a senior technician blinded to patients’ details. All positive blood slides and 10% of randomly chosen negative slides were double-checked; in case of discrepancy, slides were read by a third senior technician.
FPBSs were first stored at 4 °C until transferred to NIMPE where they were stored at −20 °C until processing. DNA extraction was done using the QIAamp DNA Micro Kit (Qiagen, Hilden Germany); a species-specific SnM-PCR to detect the four human Plasmodium species was performed [25]. Quality control was done at the Institute of Tropical Medicine, Antwerp (ITMA) by a senior technician who blindly reread 10% of randomly chosen blood samples.
Sample size
The sample size was calculated on the basis of P. vivax CQ failure rates in different sentinel sites in Central Vietnam, ranging from 0% to 5.7% [26]. Assuming 5% treatment failure at day 28 post-treatment and 10% loss to follow-up, 205 patients would provide 3% precision at a 5% significance level (“CSample” command/Epi Info6). Beyond 28 days post-treatment, it was assumed that 20% of the remaining patients would be lost to follow-up. Therefore, an additional security margin was added to start with a cohort of 250 patients, resulting in 205 patients with complete follow-up data points (i.e., D0 + 24 monthly records).
Data analysis
Data were double entered and checked in Epidata version 3.1 (The EpiData Association, Odense M, Denmark), and analyzed using STATA version 11 (Stata Corp, College Station, TX, USA). The study design was taken into account at all steps of the analysis by setting the data with household as primary sampling unit and village as strata (alpha = 0.05) using the survey (“svy”) command in STATA.
All patients entered the follow-up at D0 until censoring, which occurred when patients were lost follow-up, had a treatment failure during the first 28 days, or at the end of the monthly follow-up. Any patent recurrence after day 28 was systematically treated with CQ (3 days) and considered as an independent event from the next one. However, the SM recurrences that could be detected for several consecutive months in the same individual (but were not treated because they were identified months later) were considered as “positive PCR person-month.”
Summary statistics (accounting for study design) were used to describe the baseline characteristics at inclusion and at subsequent visits. Incidence rate of vivax recurrences were calculated by dividing the number of recurrent episodes by the person-years at risk and expressed as the number of recurrences per person-years at risk. The Wilcoxon sign rank test was used to compare the median time intervals between consecutive recurrences identified by PCR.
Time to first P. vivax recurrence post-treatment was analyzed by Kaplan–Meier survival analysis, which estimated the probability for patients to remain free of any vivax recurrence (PCR-detected, patent, or symptomatic) as well as of any P. falciparum infection (detected by PCR). To minimize the risk of vivax infections acquired outside the study area, individuals absent from the study area for two or more consecutive visits were censored at their last visit. In case of a single monthly visit missed, this was not considered as a missing month if the individual had remained within the study site, while if s/he had left the study area, it was considered as a person-time at risk censored at the time s/he was absent.
The time at risk to first recurrence was counted as days from D0 to first vivax recurrence, while for any subsequent recurrence, the time at risk was counted as days since the previous recurrence. Patients treated with either CQ (patent P. vivax recurrence) or DHA-PQ (patent P. falciparum infection) were censored for 14 days post-treatment.
The risk factor analysis was conducted both for the first and recurrent vivax recurrences. A uni- and multivariable-adjusted Cox Proportional Hazard (PH) regression model (Hazard Ratio [HR]) was used to identify significant risk factors for the first PCR-detected P. vivax recurrence (first recurrence per subject). A Cox Conditional Gap Time (CGT) model was used to assess the risk for recurrent infections per subject over time [27,28]. In this model, each recurrence was analyzed separately and stratified by recurrence order, the underlying assumption being that a subject was not at risk of a second recurrence until the first recurrence had occurred. Thus, the conditional risk set at time for recurrence k (tk) was made up for all subjects under observation at time they had had recurrence k − 1. This method considers the total number of recurrences analyzed and the time at risk between each recurrence. The extra correlation due to repeated observations within patients was accounted for (using the vce (cluster ID) option within the “stcox” command in the CGT model). Interactions were systematically checked up to order two. The PH assumption was tested for each Cox model.
Ethical clearance
Ethical clearance was obtained from the ethical committees of both the NIMPE, Hanoi and the University of Antwerp. The fundamental principles of ethics in research on human participants were upheld throughout the project. The study objectives and methods were first explained to the community leaders for their approval. Each study participant provided a written informed consent following explanation of the study objectives and follow-up procedures as well as their right to withdraw without prejudice for themselves or their families. For minors (<18 years old), parents/guardians signed the individual informed consent, while individuals aged 12–18 were asked to provide individual assent.
Results
Cohort characteristics
A total of 260 P. vivax-infected patients were identified and enrolled (21 during initial survey; 239 during PCD). Of these, 20 (7.69%) patients withdrew consent and did not complete the 10-day radical cure (Fig 1). By day 28 (first monthly visit), 9 additional patients had withdrawn consent, and 8 had treatment failure, were treated with DHA-PQ (3 days) and PQ (10 days), and were excluded as per protocol. Therefore, after day 28, the cohort comprised 223 patients, of whom 4 withdrew consent (moved out of study area), 219 completed at least 12 months, and 107 completed 24 months of follow-up. Overall, 33 patients (33/260 = 12.7%) withdrew consent, mainly (20/33 = 60.6%) during the 10-day PQ treatment. Total individual follow-up times ranged from 9 to 719 days, with a median of 628 days (IQR 508–718), and none of the participants were censored because of prolonged absence outside the study area.
CQ, chloroquine; d, day; LCF, Late Clinical Failure; LM, light microscopy; LPF, Late Parasitological Failure; m, month; PCR, polymerase chain reaction; PQ, primaquine.
Most patients were males (sex ratio M/F = 1.6), and children 3–9 years old represented almost half (43.08%) of the cohort (Table 1). The M’nong ethnic group (78.46%) in Villages 1–3 was the most represented. All adults were farmers. Most patients were very poor (lowest income category = 56.54%) and without a bed net at home (70.77%). Beside fever and/or history of fever, the most frequent symptoms at enrollment were headache and fatigue. The spleen rate was 6.15% (16/260). The mean parasite density at enrollment was 2,754.1/μl, and gametocytes were found in most infections (86.15%). By day 3, all patients had cleared infection and were symptom-free. About half of the patients were recruited in 2009 (48.46%), and the rest in 2010. The sociodemographic characteristics of patients recruited in 2009 were similar to those recruited in 2010 except for village, ethnicity, and bed net availability (S1 Table).
Characteristics of vivax recurrences
About two-thirds (76.68%, 171/223) of individuals had at least one PCR-detected recurrence (median 4 recurrences per patient [IQR 2–6]), 124 (55.61%) at least one patent recurrence (median = 2/patient; IQR 1–4), and 84 (37.67%) at least one symptomatic recurrence (median = 1; IQR 1–2) (Table 2). Most patients with PCR-detected (78.95%) and patent recurrences (65.32%) had multiple events, while those with symptomatic recurrences (63.10%) experienced mostly one event. Among the 4,604 visits during the study period, 71 (1.54%) were unscheduled because of symptomatic vivax recurrences.
Among the 748 PCR-detected P. vivax recurrences, 45% (337/748) were patent, and among these, 41% (139/337) were symptomatic. The incidence rate was 1.98 per person-years (95% CI [1.84; 2.12]) for PCR-detected, 0.89 per person-years (95% CI [0.80; 0.99]) for patent, and 0.37 per person-years (95% CI [0.31; 0.43]) for symptomatic recurrences (Table 3). Additionally, a total of 285 P. falciparum infections were detected by PCR (median: 2 episodes/person [IQR 1–3]; maximum 8/person) of which 240 (84.21%) were patent (median: 1 episode/person [IQR 1–3]; maximum 8), including 149 (/240 = 62.08%) symptomatic infections (median = 1/person [IQR 1–2]; maximum 6).
Probability of remaining free of vivax recurrence
The Kaplan–Meier probability of remaining free of recurrences by month 24 was 20.40% (95% CI [14.42; 27.13]) by PCR, 42.52% (95% CI [35.41; 49.44]) by LM, and 60.69% (95% CI [53.51; 67.11]) for symptomatic recurrences (Fig 2A). For P. falciparum recurrences, such probability was 55.61% (95% CI [48.79; 61.90]) by PCR. The median time of remaining free of PCR-detected recurrence was 174 days—i.e., half of the cohort had experienced a first recurrence by the sixth month of follow-up, while this was more than double (356 days) for patent recurrences. Most of the first events for the three types of recurrences occurred during the first year, with little changes in the cumulative risk afterwards.
Wilcoxon sign rank P-values were computed. *Second time interval was compared to first among patients with at least two recurrences; **third time interval was compared to first among patients with at least three recurrences; ***fourth interval compared to first among patients with at least four recurrences. LM, light microscopy; PCR, polymerase chain reaction; S, survival probability.
The median time to first recurrence was 118 days (IQR 59–208) by PCR, 141 days (IQR 86–237) by LM, and 181 days (IQR 97–316) for symptomatic recurrence. Moreover, the median time to first recurrence tended to be shorter (88 days [IQR 58–1,784]) among patients with multiple PCR recurrences compared to those with a single PCR recurrence (163 days [IQR 58–291]; Wilcoxon rank sum P < 0.001 [S1 Data]). Time between PCR-detected recurrences tended to decrease with increasing number of events (Fig 2B); median times between first and second recurrence (median = 78 days; IQR 43–120) and between second and third recurrence (median = 60 days; IQR 40–120) were significantly shorter than that to first recurrence (patients with at least 2 recurrences: median = 88 days, IQR 58–174, P = 0.024; patients with at least 3 recurrences: median = 88 days, IQR 58–148, P = 0.002). From the fourth to the ninth recurrence, the median time intervals remained stable at around 60 days, and from the 10th to the 13th recurrence, at 30 days. A similar pattern was observed with patent recurrences, the median time intervals stabilizing around 76 days from the third recurrence onwards (S1 Data).
Fig 3 shows a heatmap of all patients’ recurrences by month of follow-up ordered by the timing of the first recurrence. SM recurrences were frequent in almost all patients and persisted for 2–3 months (up to 7 months) before either clearing spontaneously or becoming patent or symptomatic. The frequency of symptomatic and patent recurrences decreased with increasing follow-up time, with a large majority of SM recurrences detected between month 18 and 24 of follow-up.
ID, identifier; PCR, polymerase chain reaction.
Risk factors for P. vivax recurrences
Potential risk factors for a first vivax recurrence identified by the univariable analysis, regardless of the detection method, were village, ethnicity, bed nets, year of recruitment, and having a prior (1–2 months) PCR-detected P. falciparum infection (Table 3). This applied also to all P. vivax recurrences, with the additional strong effect of a prior (1–2 months) antimalarial treatment (ACT or CQ). Since village and ethnicity were highly correlated and only Village 4 (i.e., Cadong ethnicity) had a significant effect, only ethnicity was kept in the final model.
In the final Cox PH regression model, only ethnicity, prior P. falciparum infection, and year of recruitment remained significantly associated with the risk of first PCR-detected recurrence (Table 4). In addition, for the Cox PH assumption to be verified, the analysis had to be stratified by year of recruitment. The main risk factor identified for both first and all vivax recurrences was a prior P. falciparum infection, which increased by 5- to 7-fold the hazard of vivax recurrence. However, this effect was not homogenously distributed among all patients because there were significant interactions with ethnicity/Village 4 (final model for first recurrence) and prior antimalaria treatment (final model for all recurrences). Indeed, among patients recruited in 2009, the hazard of having first PCR-detected recurrence after a prior falciparum infection was 3-fold higher among Cadong (HR = 2.98; 95% CI [1.48; 5.99]; P = 0.002) and at least 6-fold higher among M’nong (HR = 99.19; 95% CI [6.16; 1,597.55]; P = 0.001) patients. Conversely, the hazard of first PCR recurrence was 2-fold higher in Cadong than in M’nong patients (HR = 2.21; 95% CI [1.33; 3.67]; P = 0.002), but only in those without a prior falciparum infection. In 2010, few Cadong patients were recruited, and thus the effect of ethnicity/Village 4 was no longer significant; thus, only prior falciparum infection remained significantly associated with the hazard of both first PCR-detected and patent recurrence. None of the patients recruited in 2010 had a falciparum infection before their first symptomatic recurrence.
Similarly, the final Cox CGT models for all recurrences among patients recruited in 2009 and 2010 showed that the hazard of PCR-detected vivax recurrences after a prior falciparum infection (versus no P. falciparum) was 5- to 6-fold higher, but only in patients without previous antimalaria treatment (i.e., SM falciparum infections). Conversely, the effect of prior antimalarial treatment (versus no treatment) was very strong in both recruitment years (respectively, HR = 6.35 and HR = 8.94; P < 0.001), but only in individuals without a prior falciparum infection, i.e., patent recurrences treated with CQ. For patent and symptomatic recurrences, a similar pattern was observed for patients recruited in 2009.
The hazard of patent recurrences was significantly higher (HR = 6.34; P < 0.001) in patients with a prior untreated falciparum infection, while the opposite effect was observed in patients who were treated (HR = 0.23; P < 0.001). For patients recruited in 2010, only prior antimalarial treatment remained significantly associated with the risk of patent (HR = 13.83; P < 0.001) and symptomatic recurrences (HR = 6.05; P < 0.001). Fig 4A shows the study timeline with corresponding climatic data; Fig 4B shows the evolution of the cohort size by calendar month together with the proportion of study patients with vivax recurrences (PCR-detected, patent, and symptomatic) as well as the contemporary monthly incidence risk of P. falciparum infections (Fig 4B) as detected by PCD and ACD (LM) in the four study villages. The first mass distribution of long-lasting insecticidal nets (LLINs) in the study area occurred in September 2010 and was followed by a steady decrease in both the risk of P. falciparum infections in the study area and of P. vivax recurrences among cohort patients (Fig 4B). The monthly incidence risk of P. falciparum decreased to almost zero from February 2011 until the end of follow-up; similarly, the risk of all PCR-detected and patent vivax recurrences steadily decreased to 7% and 1%, respectively, by December 2011.
ACD, active case detection; LM, light microscopy; PCD, passive case detection.
Discussion
Post-treatment vivax recurrences were extremely frequent, especially SM infections, which tended to occur earlier than patent or symptomatic infections. Prior falciparum infection was the main risk factor for all types of vivax recurrence (first and all recurrences; PCR-detected, patent, and symptomatic recurrences) together with prior CQ treatment for the risk of repeated recurrences. The drop-out rate was low given the length of follow-up, and most withdrawals occurred during the 10-day PQ treatment.
Although collected about ten years ago, our results remain relevant for Vietnam where the current challenges for malaria elimination include i) the increasing relative frequency of P. vivax malaria, from 25% of all cases in 2009 to 51% in 2015 [9]; and ii) the spread of artemisinin resistance [29]. In addition, the downward trend of annual malaria incidence has changed as the total number of confirmed cases started to increase in 2017, with a 16% increase observed in 2018 (4,813 cases) compared to 2016 (4,161 cases) [30,31], mostly due to P. falciparum; cases were observed mainly in a few provinces in Central Vietnam where artemisinin resistance has been reported [32].
Despite PQ directly observed treatment, 8% of patients (most withdrawals) withdrew their consent during the 10-day treatment, emphasizing the difficulty of achieving good compliance for PQ radical treatment [18,33]. In our cohort, PQ was only used at the time of recruitment and at the end of each patient’s follow-up. This is probably similar to the true situation in remote areas in Vietnam, including our study area, where the health staff is reluctant to administer PQ in the absence of diagnostic tests for G6PD deficiency [16]. Tafenoquine, a new 8-aminoquinoline with prolonged half-life recently approved by the US Food and Drug Administration [3], will not change this situation because it will be administered only after performing a quantitative G6PD rapid test [34], although the single-dose scheme will improve compliance.
Incidence of PCR-detected vivax recurrence was >2-fold higher than that of patent recurrence and 5-fold higher than symptomatic recurrence, suggesting a much larger infection reservoir than estimated by standard diagnostic tools (LM/RDT) and current surveillance strategies (PCD and Re-ACD). Although some of the PCR-detected infections at different time points may be the same infection, possibly overestimating the incidence of recurrent infections, their persistence for long periods of time (up to 7 months) contributes to maintaining transmission. Our results corroborate the findings of a recent cohort study in Central Vietnam in which the median duration of SM P. vivax infections was 6 months (IQR 3–9), with parasite densities oscillating from ultralow to high density between monthly samples [35]. Since SM vivax infections with gametocytes can successfully infect mosquitoes [36,37], all PCR-detected vivax infections should be considered as potentially infectious, regardless of parasite density or presence of gametocytes at time of sampling, and effectively treated in order to interrupt transmission.
Although the annual malaria incidence reported by PCD has further decreased since 2011, with no falciparum malaria cases since 2016 (S2 Table), the current proportion of SM vivax infections may be similar to that observed in our study. Indeed, during the last 6 months of our cohort study, in 2011, when malaria transmission had decreased substantially, PCR-detected vivax recurrences were significantly more frequent than patent or symptomatic infections. In a recent systematic review [38], SM vivax infections probably represent about 67% of all PCR-detected vivax infections and are negatively correlated, as for P. falciparum [22], with the prevalence of patent vivax infections. This suggests that, despite the observed reduction of patent vivax and falciparum cases since 2011, the proportion of SM vivax infections may not have changed; this should be further confirmed by large-scale prevalence surveys using both PCR and LM.
P. falciparum infections were identified as the main risk factor for first and repeated vivax recurrences. The triggering effect of symptomatic P. falciparum infections on P. vivax relapses has been known since the 1920s [39] and has been reported in other cohort studies in Asia [40–42]. Interestingly reported and for the first time, to our knowledge, in our study, falciparum infections were a strong risk factor for repeated vivax recurrences, but only in patients untreated (with either CQ or DHA-PQ) at the previous visit, while if treated, falciparum infection tended to be protective. This suggests that chronic SM falciparum infections can trigger vivax recurrences, while antimalarial treatment, by clearing falciparum infections, prevents vivax recurrence. Though it remains unclear whether P. falciparum triggered vivax relapses or suppressed vivax parasites during coinfections, our results suggest that chronic SM falciparum infections could trigger vivax recurrences, hence contributing to the human reservoir of P. vivax infections. This, together with previous evidence [40–42], calls for new and integrated treatment strategies tackling both species simultaneously, e.g., treating all malaria infections, regardless of the species, with one ACT and antirelapse drug (PQ or tafenoquine) [43]. The increased risk of vivax recurrence related to previous antimalarial treatment in patients without prior falciparum infection suggests CQ vivax resistance as confirmed earlier in the same cohort [44].
It is not possible to distinguish between new infections and relapses within our cohort. Nevertheless, assuming a similar incidence of P. falciparum and P. vivax new infections, the difference between the PCR-detected vivax and falciparum infections’ Kaplan–Meier curves should provide the relative contribution of vivax relapses, i.e., around 50%. This compares well with the results of a trial comparing artesunate alone or combined with PQ in Papua New Guinea, where P. falciparum and P. vivax sporozoite rates in the local vector population were comparable [45].
The total PQ dose administered (5.0 mg/kg) was lower than that recommended by WHO, i.e., 7.0 mg/kg [46]. A recent systematic review reported that a total PQ dose ≥5.0 mg/kg had the lowest risk of recurrence (median = 0%; range 0%–15%) than the very low (≤2.5 mg/kg) and low total dose (>2.5 and <5.0 mg/kg) [5]. This confirmed a previous review concluding that 0.5 mg/kg/day for 14 days was the recommended dose in nonpregnant subjects without G6PD deficiency [47]. Although it is impossible to know whether a 2.0 mg/kg difference in PQ dosage would have affected the risk of first recurrence, the synergistic effect of combining CQ with PQ could have been hampered by P. vivax CQ resistance, which was observed in our cohort [44]. This is also suggested by the strong deleterious effect of prior CQ treatment on all types of vivax recurrence in 2009 and for the PCR-detected recurrences in 2010. CQ resistance in P. vivax should be carefully monitored in Vietnam, and treatment for vivax patients changed to an ACT in case of substantial resistance [48,49].
A major limitation of this study is the absence of a comparator arm in which all vivax recurrences would have been systematically treated with PQ; this would have enabled estimating the relative contribution of vivax relapses to recurrent infections [45]. However, such a design was not chosen given the low incidence rates of symptomatic vivax infections (main inclusion criteria), which would have required an even longer recruitment period that could not be supported by the available budget.
As mentioned earlier, the age of the data represents another limitation because malaria continued to decrease steadily after the end of the study, with no falciparum cases reported in 2016 and 2017 (S2 Table). However, as mentioned earlier, P. falciparum is on the rise in Central Vietnam, and the province of Quang Nam also reported an increase in P. falciparum cases in 2018 (+4% compared to 2017). Therefore, even though the data presented in this paper are not directly comparable to the current situation in the study area, the challenges they underline, i.e., the elimination of P. vivax and of artemisinin-resistant P. falciparum in Central Vietnam, remain more than ever relevant.
Because hemoglobin (Hb) levels were not monitored during the follow-up, it was not possible to determine the risk and severity of anemia related to vivax recurrences, especially for the prolonged SM vivax infections, which was another major limitation. Finally, in our statistical analysis, the different sampling schemes underlying scheduled and unscheduled visits were not considered. Nevertheless, the related bias is probably minor because of the intense monthly sampling and the moderate number of unscheduled visits (4.533 scheduled visits versus 71 unscheduled visits). Further research should focus on developing appropriate methodology taking into account the outcome-dependent sampling underlying unscheduled events, especially in case of larger time gaps between scheduled visits [50].
In conclusion, this large cohort study identified a high burden of post-treatment vivax infections, mainly asymptomatic and SM, despite a high-dose PQ regimen (total 5 mg/kg) administered under direct observation. The main risk factors for vivax recurrence were prior P. falciparum infections—particularly the SM ones, which were untreated because they were not detected immediately after sampling and persisted for several months—and prior CQ monotherapy, indicating vivax drug resistance. When considering Vietnam reverted to the standard PQ dose of 0.25 mg/kg/day for 14 days, which is often not prescribed or supervised, the current vivax parasite reservoir, despite the decreasing number of reported cases, may be larger than estimated by standard diagnostic tools. This calls for improved treatment (e.g., single-dose tafenoquine with G6PD testing) and surveillance strategies (including the detection of SM infections) tackling simultaneously P. falciparum and P. vivax infections if malaria elimination is to be achieved by 2030 as endorsed by the Vietnamese government.
Supporting information
S1 Table. Baseline characteristics of the 260 patients enrolled in the cohort split by the year of recruitment.
https://doi.org/10.1371/journal.pmed.1002784.s001
(XLSX)
S2 Table. Incidence of malaria in four study villages from 2012 to 2018.
https://doi.org/10.1371/journal.pmed.1002784.s002
(XLSX)
S1 Data. Data of Tables 2, 3 and 4 and Figs 2 and 3.
https://doi.org/10.1371/journal.pmed.1002784.s003
(XLSX)
S1 Text. Protocol: P. vivax epidemiology and control in Vietnam.
https://doi.org/10.1371/journal.pmed.1002784.s006
(DOC)
S2 Text. STROBE statement.
STROBE Statement is a checklist of 22 items that we consider essential for good reporting of cohort studies. These items relate to the article′s title and abstract (item 1), the introduction (items 2 and 3), methods (items 4–12), results (items 13–17) and discussion sections (items 18–21), and other information (item 22 on funding). STROBE, Strengthening The Reporting of OBservational Studies in Epidemiology.
https://doi.org/10.1371/journal.pmed.1002784.s007
(DOC)
Acknowledgments
We would like to thank the M’nong and Cadong communities in the Tra Leng and Tra Don communes for their contribution to the present study, as well as the local and provincial authorities for their strong and continuous support.
References
- 1. Guerra CA, Howes RE, Patil AP, Gething PW, Van Boeckel TP, Temperley WH, et al. The international limits and population at risk of Plasmodium vivax transmission in 2009. PLoS Negl Trop Dis. 2010;4(8):e774. pmid:20689816
- 2. Baird KJ, Maguire JD, Price RN. Diagnosis and treatment of Plasmodium vivax malaria. Adv Parasitol. 2012;80:203–70. pmid:23199489
- 3.
Fleming M. FDA Approves Tafenoquine, First New P Vivax Malaria Treatment in 60 Years. Contagion Live. JUL 23, 2018 [cited 2018 Aug 16]. https://www.contagionlive.com/news/fda-approves-tafenoquine-first-new-treatment-of-p-vivax-malaria-in-60-years.
- 4. Baird JK. Resistance to chloroquine unhinges vivax malaria therapeutics. Antimicrob Agents Chemother. 2011;55(5):1827–30. pmid:21383088
- 5. John GK, Douglas NM, von Seidlein L, Nosten F, Baird JK, White NJ, et al. Primaquine radical cure of Plasmodium vivax: a critical review of the literature. Malar J. 2012;11:280. pmid:22900786
- 6.
WHO. Guidelines for the treatment of malaria. Third edition. April 2015 [cited 2018 Mar 10]. http://www.who.int/malaria/publications/atoz/9789241549127/en/.
- 7.
WHO. World Malaria Report 2017 [cited 2018 Aug 16]. http://apps.who.int/iris/bitstream/handle/10665/259492/9789241565523-eng.pdf;jsessionid=0F04E7242F1B3AEFFCF70870D09674A1?sequence=1.
- 8.
The Government of Vietnam. National Strategy for malaria control and elimination 2011–2020 and orientation for 2030. Decision No. 1920/QĐ-TTg dated of 27 October 2011 of the Prime Minister [cited 2017 Dec 21]. http://www.wpro.who.int/vietnam/topics/malaria/vnm_national_strategy_for_malaria_controlnelimination.pdf.
- 9.
National Malaria Control and Elimination Program (NMECP) in Vietnam. Malaria report from 2011 to 2015, Intenal Report. NIMPE, Hanoi 2016.
- 10.
National Malaria Control and Elimination Program (NMECP) in Vietnam. Malaria report from 1999 to 2016, Internal Report. NIMPE, Hanoi. 2016.
- 11.
Vietnam Ministry of Health: Guideline Surveillance and prevention of malaria disease. Hanoi, 2016 [cited 2017 Dec 21]. http://vncdc.gov.vn/files/Notice/2016/3/quyet-dinh-741_qdbyt-ngay-0232016-huong-dan-giam-sat-va-phong-chong-sot-ret(2).pdf.
- 12. Cao J, Sturrock HJ, Cotter C, Zhou S, Zhou H, Liu Y, et al. Communicating and monitoring surveillance and response activities for malaria elimination: China′s "1-3-7" strategy. PLoS Med. 2014;11(5):e1001642. pmid:24824170
- 13.
Ministry of Health. Guideline of diagnosis, treatment and prevention malaria disease 2016 [cited 2017 Dec 25]. https://thuvienphapluat.vn/van-ban/The-thao-Y-te/Quyet-dinh-4845-QD-BYT-Huong-dan-chan-doan-dieu-tri-benh-Sot-ret-2016-321533.aspx.
- 14.
Ministry of Health. Guideline of diagnosis, treatment and prevention malaria disease 2007 [cited 2017 Dec 25]. https://thuvienphapluat.vn/van-ban/The-thao-Y-te/Quyet-dinh-339-QD-BYT-Huong-dan-chuan-doan-dieu-tri-benh-sot-ret-93442.aspx.
- 15.
Ministry of Health. Guideline of diagnosis, treatment and prevention malaria disease 2009 [cited 2017 Dec 25]. https://vanbanphapluat.co/quyet-dinh-3232-qd-byt-2013-huong-dan-chan-doan-dieu-tri-benh-sot-ret.
- 16.
Thanh PV. Epidemiology of Plasmodium vivax malaria in Central Vietnam, PhD thesis 2018 [cited 2018 Oct 10]. https://dial.uclouvain.be/pr/boreal/en/object/boreal%3A202970.
- 17. Grietens KP, Soto V, Erhart A, Ribera JM, Toomer E, Tenorio A, et al. Adherence to 7-day primaquine treatment for the radical cure of P. vivax in the Peruvian Amazon. Am J Trop Med Hyg. 2010;82(6):1017–23. pmid:20519594
- 18. Douglas NM, Poespoprodjo JR, Patriani D, Malloy MJ, Kenangalem E, Sugiarto P, et al. Unsupervised primaquine for the treatment of Plasmodium vivax malaria relapses in southern Papua: A hospital-based cohort study. PLoS Med. 2017;14(8):e1002379. pmid:28850568
- 19. Thanh PV, Van Hong N, Van Van N, Van Malderen C, Obsomer V, Rosanas-Urgell A, et al. Epidemiology of forest malaria in Central Vietnam: the hidden parasite reservoir. Malar J. 2015;14:86. pmid:25880664
- 20. Garros C, Van Nguyen C, Trung HD, Van Bortel W, Coosemans M, Manguin S. Distribution of Anopheles in Vietnam, with particular attention to malaria vectors of the Anopheles minimus complex. Malar J. 2008;7:11. pmid:18190697
- 21. Van Bortel W, Trung HD, Hoi le X, Van Ham N, Van Chut N, Luu ND, et al. Malaria transmission and vector behaviour in a forested malaria focus in central Vietnam and the implications for vector control. Malar J. 2010;9:373. pmid:21182774
- 22. Okell LC, Bousema T, Griffin JT, Ouedraogo AL, Ghani AC, Drakeley CJ. Factors determining the occurrence of submicroscopic malaria infections and their relevance for control. Nat Commun. 2012;3:1237. pmid:23212366
- 23. Rubio JM, Post RJ, van Leeuwen WM, Henry MC, Lindergard G, Hommel M. Alternative polymerase chain reaction method to identify Plasmodium species in human blood samples: the semi-nested multiplex malaria PCR (SnM-PCR). Trans R Soc Trop Med Hyg. 2002;96 Suppl 1:S199–204.
- 24. Van den Eede P, Van HN, Van Overmeir C, Vythilingam I, Duc TN, Hung le X, et al. Human Plasmodium knowlesi infections in young children in central Vietnam. Malar J. 2009;8:249. pmid:19878553
- 25. Nguyen HV, van den Eede P, van Overmeir C, Thang ND, Hung le X, D′Alessandro U, et al. Marked age-dependent prevalence of symptomatic and patent infections and complexity of distribution of human Plasmodium species in central Vietnam. Am J Trop Med Hyg. 2012;87(6):989–95. pmid:23128294
- 26.
Tien NT. Reviewing the monitoring of antimalarial drug efficacy and resistance in sentinel sites in Vietnam. International Colloquium ITM/NIMPE 2007 [cited 2017 Dec 21]. Abstract book p. 71 [71]. https://studylib.net/doc/8823884/abstract-book.
- 27. Amorim LD, Cai J. Modelling recurrent events: a tutorial for analysis in epidemiology. Int J Epidemiol. 2015;44(1):324–33. pmid:25501468
- 28.
Cleves M. Analysis of multiple failure-time survival data [cited 2017 Dec 21]. http://www.stata.com/support/faqs/statistics/multiple-failure-time-data/.
- 29.
WHO. Status report on artemisinin resistance and ACT efficacy (August 2018) 2018 [cited 2018 Aug 16]. https://www.who.int/malaria/publications/atoz/artemisinin-resistance-august2018/en/.
- 30.
National Malaria Control and Elimination Program (NMECP) in Vietnam. Annual Malaria report, Internal report. NIMPE, Hanoi. 2018.
- 31.
National Malaria Control and Elimination Program (NMECP) in Vietnam. Annual Malaria report, Internal report. NIMPE, Hanoi. 2017.
- 32. Phuc BQ, Rasmussen C, Duong TT, Dong LT, Loi MA, Menard D, et al. Treatment Failure of Dihydroartemisinin/Piperaquine for Plasmodium falciparum Malaria, Vietnam. Emerg Infect Dis. 2017;23(4):715–7. pmid:28322709
- 33. Thriemer K, Bobogare A, Ley B, Gudo CS, Alam MS, Anstey NM, et al. Quantifying primaquine effectiveness and improving adherence: a round table discussion of the APMEN Vivax Working Group. Malar J. 2018;17(1):241. pmid:29925430
- 34.
Medicines for Malaria Venture (MMV). Tafenoquine 2018 [cited 2018 Oct 20]. https://www.mmv.org/related-story-type/tafenoquine.
- 35. Nguyen TN, von Seidlein L, Nguyen TV, Truong PN, Hung SD, Pham HT, et al. The persistence and oscillations of submicroscopic Plasmodium falciparum and Plasmodium vivax infections over time in Vietnam: an open cohort study. Lancet Infect Dis. 2018;18(5):565–72. pmid:29398388
- 36. Sattabongkot J, Maneechai N, Phunkitchar V, Eikarat N, Khuntirat B, Sirichaisinthop J, et al. Comparison of artificial membrane feeding with direct skin feeding to estimate the infectiousness of Plasmodium vivax gametocyte carriers to mosquitoes. Am J Trop Med Hyg. 2003;69(5):529–35. pmid:14695091
- 37. Collins WE, Jeffery GM, Roberts JM. A retrospective examination of the effect of fever and microgametocyte count on mosquito infection on humans infected with Plasmodium vivax. Am J Trop Med Hyg. 2004;70(6):638–41. pmid:15211005
- 38. Cheng Q, Cunningham J, Gatton ML. Systematic review of sub-microscopic P. vivax infections: prevalence and determining factors. PLoS Negl Trop Dis. 2015;9(1):e3413. pmid:25569135
- 39. Sinton JA. Some lacunae in our knowledge of the malaria parasite. Indian J Malariol. 1955;9(4):229–45. pmid:13306281
- 40. White NJ. Determinants of relapse periodicity in Plasmodium vivax malaria. Malar J. 2011;10:297. pmid:21989376
- 41. Douglas NM, Nosten F, Ashley EA, Phaiphun L, van Vugt M, Singhasivanon P, et al. Plasmodium vivax recurrence following falciparum and mixed species malaria: risk factors and effect of antimalarial kinetics. Clin Infect Dis. 2011;52(5):612–20. pmid:21292666
- 42. Looareesuwan S, White NJ, Chittamas S, Bunnag D, Harinasuta T. High rate of Plasmodium vivax relapse following treatment of falciparum malaria in Thailand. Lancet. 1987;2(8567):1052–5. pmid:2889965
- 43. Price RN, Douglas NM, Anstey NM, von Seidlein L. Plasmodium vivax treatments: what are we looking for? Curr Opin Infect Dis. 2011;24(6):578–85. pmid:21986614
- 44. Thanh PV, Hong NV, Van NV, Louisa M, Baird K, Xa NX, et al. Confirmed Plasmodium vivax Resistance to Chloroquine in Central Vietnam. Antimicrobial agents and chemotherapy. 2015;59(12):7411–9. pmid:26392501
- 45. Betuela I, Rosanas-Urgell A, Kiniboro B, Stanisic DI, Samol L, de Lazzari E, et al. Relapses contribute significantly to the risk of Plasmodium vivax infection and disease in Papua New Guinean children 1–5 years of age. J Infect Dis. 2012;206(11):1771–80. pmid:22966124
- 46.
WHO. Guidelines for the treatment of malaria. Second Edition. 2010 [cited 2017 Dec 21]. http://apps.who.int/medicinedocs/documents/s19105en/s19105en.pdf.]
- 47. Baird JK, Hoffman SL. Primaquine therapy for malaria. Clin Infect Dis. 2004;39(9):1336–45. pmid:15494911
- 48. Gogtay N, Kannan S, Thatte UM, Olliaro PL, Sinclair D. Artemisinin-based combination therapy for treating uncomplicated Plasmodium vivax malaria. Cochrane Database Syst Rev. 2013(10):CD008492. pmid:24163021
- 49. Douglas NM, Anstey NM, Angus BJ, Nosten F, Price RN. Artemisinin combination therapy for vivax malaria. Lancet Infect Dis. 2010;10(6):405–16. pmid:20510281
- 50.
Mugenyi L, Herzog SA, Hens N, Abrams S. Modelling longitudinal binary outcomes with outcome- dependent observation times: an application to a malaria cohort study. BioRxiv 252411 [Preprint]. 2018 [cited 2018 Oct 30]. https://www.biorxiv.org/content/10.1101/252411v1.