Factors affecting haemoglobin dynamics in African children with acute uncomplicated Plasmodium falciparum malaria treated with single low-dose primaquine or placebo

Background Single low-dose primaquine (SLDPQ) effectively blocks the transmission of Plasmodium falciparum malaria, but anxiety remains regarding its haemolytic potential in patients with glucose-6-phopshate dehydrogenase (G6PD) deficiency. We, therefore, examined the independent effects of several factors on haemoglobin (Hb) dynamics in falciparum-infected children with a particular interest in SLDPQ and G6PD status. Methods This randomised, double-blind, placebo-controlled, safety trial was conducted in Congolese and Ugandan children aged 6 months–11 years with acute uncomplicated P. falciparum and day (D) 0 Hbs ≥ 6 g/dL who were treated with age-dosed SLDPQ/placebo and weight-dosed artemether lumefantrine (AL) or dihydroartemisinin piperaquine (DHAPP). Genotyping defined G6PD (G6PD c.202T allele), haemoglobin S (HbS), and α-thalassaemia status. Multivariable linear and logistic regression assessed factor independence for continuous Hb parameters and Hb recovery (D42 Hb > D0 Hb), respectively. Results One thousand one hundred thirty-seven children, whose median age was 5 years, were randomised to receive: AL + SLDPQ (n = 286), AL + placebo (286), DHAPP + SLDPQ (283), and DHAPP + placebo (282). By G6PD status, 284 were G6PD deficient (239 hemizygous males, 45 homozygous females), 119 were heterozygous females, 418 and 299 were normal males and females, respectively, and 17 were of unknown status. The mean D0 Hb was 10.6 (SD 1.6) g/dL and was lower in younger children with longer illnesses, lower mid-upper arm circumferences, splenomegaly, and α-thalassaemia trait, who were either G6PDd or heterozygous females. The initial fractional fall in Hb was greater in younger children with higher D0 Hbs and D0 parasitaemias and longer illnesses but less in sickle cell trait. Older G6PDd children with lower starting Hbs and greater factional falls were more likely to achieve Hb recovery, whilst lower D42 Hb concentrations were associated with younger G6PD normal children with lower fractional falls, sickle cell disease, α-thalassaemia silent carrier and trait, and late treatment failures. Ten blood transfusions were given in the first week (5 SLDPQ, 5 placebo). Conclusions In these falciparum-infected African children, posttreatment Hb changes were unaffected by SLDPQ, and G6PDd patients had favourable posttreatment Hb changes and a higher probability of Hb recovery. These reassuring findings support SLDPQ deployment without G6PD screening in Africa. Trial registration The trial is registered at ISRCTN 11594437. Supplementary Information The online version contains supplementary material available at 10.1186/s12916-023-03105-0.

ARPf has been present in SE Asia for the past 15 years and has seen significant declines in efficacy of the partner drugs used in artemisinin-based combination treatments (ACTs) [7][8][9][10][11].ARPf is now gaining a foothold in eastern Africa [12,13] with concomitant evidence of a decline in the sensitivity of lumefantrine, a commonly used ACT [14,15].SLDPQ offers an opportunity to counter this threat, but concern remains that SLDPQ could increase the risk of life-threatening haemolysis and blood transfusions (primaquine's most feared toxicity) in patients with glucose-6-phosphate dehydrogenase deficiency (G6PDd) with a consequential loss of credibility of this strategy by malaria control programmes (MCPs).
There is growing evidence that SLDPQ is well tolerated across the age spectrum in patients, healthy individuals, and asymptomatic P. falciparum carriers with the African A-and the more severe Southeast Asian G6PD variants [27][28][29][30][31][32][33].The mean difference in nadir Hb concentrations between SLDPQ treated G6PDd and G6PD normal malaria patients is ~ 1 g/dL [27,28,31], but at the higher primaquine doses used by Shekalaghe et al., 0.75 − 1.0 mg/kg, the mean difference was 2 g/dL between the G6PDd and G6PD normal children [34].
Our group recently published the largest study to date on the safety of a bespoke, age-dosed regimen of SLDPQ [35] in young G6PDd children with uncomplicated P. falciparum.SLDPQ was associated with the same blood transfusion rate as placebo (0.9%) and a very similar tolerability profile [36].Extending our analysis of this study, we sought to determine the independent effects of several factors, including G6PD and sickle cell status, thalassaemia, and SLDPQ on the course of haemoglobin over time.
The study sites, the Kinshasa Mahidol Oxford Research Unit (outside Kinshasa) and the Mbale Regional Referral Hospital (Mbale, SE Uganda), are in areas of hyperendemic, perennial malaria transmission of P. falciparum and both sites serve urban and rural populations; children < 5 years are the main risk group for malaria.

Statistical analysis
The analysis population included all patients who received at least one dose of study drug and had at least one posttreatment Hb measurement.As per the WHO, anaemia of any degree was defined as a Hb < 11 g/dL in children < 5 years and < 11.5 mg/dL in those aged 5-11 years [43]; adapting the WHO definition, we defined moderate anaemia as a Hb < 8 g/dL for all ages.
Multivariable linear regression models assessed a range of clinically relevant variables that, collectively, have been assessed in previous studies [21-29, 44, 45]; these were age, sex, length of illness, mid upper arm circumference (MUAC), D0 Hb, D0 parasitaemia, splenomegaly, hepatomegaly, G6PD, HbS and α-thalassaemia genotype, ACT, SDLPQ/placebo, Hb fractional fall, and late treatment failure.We assessed their independent effects on the (i) D0 Hb concentration, (ii) initial decline in Hb concentration from D0 to the day of the nadir Hb in the first 14 days, (iii) nadir and D42 Hb concentrations, and (iv) increase in Hb from the nadir to D42 concentration, the total malaria attributable fall in Hb following treatment (MAFt, Additional file 1: Fig S1) [45], excluding variables where needed, e.g.ACT, SLDPQ/placebo, and late treatment failure were not used to assess the baseline Hb.In linear regression models, the slopes and the corresponding 95% confidence intervals for slopes are reported.Collinearity was assessed using the variance inflation factor (VIF) and variables with VIF > 5 would be excluded from the models.However, all variables had low VIFs (all < 3).
The factors associated with Hb recovery, defined as a D42 Hb > D0 Hb, were determined by logistic regression (age, sex, length of illness, MUAC, D0 Hb, D0 parasitaemia, splenomegaly, hepatomegaly, G6PD, HbS and α-thalassaemia genotype, ACT, SDLPQ/placebo, Hb fractional fall and late treatment failure).In these models, the odds ratios and the corresponding 95% confidence intervals for odds ratios are reported.
The times to Hb recovery by G6PD status were determined by the Kaplan-Meier survival approach, and potential explanatory variables for times to Hb recovery were explored by the Cox proportional hazards regression model.In these models, we used the scaled Schoenfeld residuals to test the proportional hazards assumption.Variables not meeting the proportional hazards assumptions were excluded from the analysis.The hazards ratios and the corresponding 95% confidence intervals for hazard ratios are reported.P-values have also been reported for all models considered.
A subgroup analysis (n = 258) is presented of patients with PK data to assess the effects of (i) the maximum primaquine concentration (C max ), (ii) primaquine exposure, i.e. area under the drug concentration time curve from 0 to the last PK sample (AUC 0-t ), and (iii) mg/kg dose and CYP2D6 activity score, inferred from the CYP2D6 genotype [46]; C max and AUC 0-t were determined by non-compartmental analysis.All statistical analyses were carried out using Stata v17 (Stata Corporation, Texas, USA) and were two sided.A P-value < 0.05 denoted statistical significance.
The mean Hb fell initially to reach its mean nadir concentration on D2 and recovered by D14 in all children and when stratified by G6PD status (Fig. 1); the nadir Hb occurred as early as 6 h in 4.7% of children (54/1137, Fig. 2).A total of 3/1137 (0.3%) of children dropped their Hb to < 5 g/dL; all were transfused.

Factors associated with baseline, nadir, and fractional fall in Hb
The baseline Hb was lower in younger children with a longer illness history, a lower MUAC, splenomegaly, who were homozygous α-thalassaemia, and either G6PDd or heterozygous females (Table 2).Although HbSS was associated with a significantly lower baseline Hb concentration, this result should be interpreted with caution because only 3 children were HbSS.The nadir Hb was higher with a higher baseline Hb, increasing age, and HbAS and lower in those with longer illness duration and higher baseline parasitaemias (Table 2).
The median fractional fall in Hb was 11.7% (IQR 6.8 − 17.5) with maximum falls of 37% (heterozygous females), 40.4% (G6PDd), and 57.9% (G6PD normal); these falls were greater in younger children with higher D0 Hbs and D0 parasitaemias and longer illnesses but less in HbAS (Table 2).The D0 Hb had the greatest influence on the fractional fall, but HbAS had a greater protective effect than increasing age.

Factors associated with D42 Hb concentration and Hb recovery
A more robust MAFt response was seen in older G6PDd children with a longer illness duration and higher D0 Hb and D0 parasitaemia but was lower with a greater fractional fall in Hb and sickle cell disease (n = 3) vs. Hb AA.A higher D42 Hb concentration was seen in older G6PDd children with a greater fractional fall but was lower in patients with sickle cell disease (vs.Hb AA), heterozygous and homozygous α-thalassaemia, and those with late treatment failure.By D42, 807 of 1066 (75.7%) children recovered their Hb and were more likely to be older with G6PDd, lower D0 Hbs, and greater Hb fractional falls.The time to Hb recovery ranged from 6 h − 42 days for a median of 14 (IQR 2-21) days and 7 (2 − 21) days in the G6PDd group.However, there was an interaction between baseline Hb and G6PD status on the time to haemoglobin recovery, suggesting a differential time to recovery in the G6PD-deficient group that depended on Hb concentration (Additional file 2: Tab S1).Therefore, we stratified these analyses by the presence/absence of anaemia of any degree to understand better the effect of the observed significant interaction.The time to Hb recovery was significantly faster in the G6PDd children without anaemia, but there was no difference by G6PD status in the anaemic group (Fig. 3).This was confirmed in the Cox proportional hazards regression model, which also showed that there were no other explanatory variables (Additional file 2: Tab S1).MUAC, fractional fall, splenomegaly, and log baseline parasitaemia violated the proportional hazards assumption in the univariate analyses and were excluded from the multivariable analysis (all were significant in the univariate analysis).
Being on SLDPQ was not an independent explanatory factor for the time to and probability of Hb recovery by D42 and the D42 Hb concentration.

Effect of primaquine pharmacokinetics and activity score
For all the postbaseline parameters of haemoglobin dynamics, there was no effect of primaquine when modelled either as a categorical factor, i.e.ACT + SLDPQ vs. ACT alone in all patients (Table 2), or in the PK subgroup (Additional file 3: Tab S2) when expressed as the mg/kg dose with CYP2D6 activity score (Additional file 3: Tab S3-8), AUC 0-t (Additional file 3: Tab S9-14), or as the C max (Additional file 3: Tab S15-20).

Hb dynamics in children with normal Hbs or moderate anaemia at baseline
There were 407 patients without baseline anaemia and 76 with moderate anaemia (Hb < 8 g/dL).G6PD status is detailed in Fig. 4.There was a very small, initial decline in the mean Hb and more rapid recovery and robust mean MAFt in the moderately anaemic children, with a greater response in the G6PDd group.

Blood transfusions
There were 10 blood transfusions (0.9%) in the first week, equally distributed between SLDPQ and placebo, and one transfusion on D31 in a patient with recurrent parasitaemia (Table 3).The 3 children with Hbs < 5 g/dL had fractional falls of 23 to 58%.None of the HbAS or HbSS children were transfused.

Discussion
MCPs with residual anxiety regarding the safety of SLDPQ should be reassured by the substantial evidence from our earlier randomised, placebo-controlled trial, which showed that the toxicity of age-dosed SLDPQ was similar to that of placebo and significantly reduced gametocyte carriage [36].Moreover, despite the resulting higher PQ exposure [37] compared to SLDPQ in asymptomatic P. falciparum carriers [47], SLDPQ was well tolerated and, as we show here, had no deleterious effects on any of the posttreatment Hb parameters.
In this analysis, we detail the factors influencing Hb dynamics to guide clinicians and MCPs on identifying particular patients who may be at increased risk of developing clinically significant and, possibly, life-threatening anaemia and those who may not recover their Hb by 6 weeks.
At baseline, children with a lower Hb concentration were younger, had longer illness durations, lower MUACs, an increased prevalence of splenomegaly, and were more frequently affected by either homozygous α-thalassaemia, HbSS, or G6PDd (either as heterozygous females or hemizygous males/homozygous females).Consistent with previous clinical descriptions [48], baseline Hb appeared lowest among children with HbSS, although the small numbers affected precluded meaningful further analysis.Given the high frequency of HbSS of ~ 1-2% at birth [49,50], the low frequency of HbSS among children recruited to this trial (3/1137; 0.26%) is striking, particularly given our recent observation that ~ 30% of children recruited to a recent trial of severe anaemia (one site was Mbale) had HbSS [51].This low HbSS rate can be explained by several factors.Children with Hbs < 6 g/dL were excluded; HbSS-affected children [51,52], like those with HbAS [53], are innately resistant to P. falciparum infections but suffer high rates (> 50%) of early-life mortality [54].Some of our findings, notably age, length of illness, and splenomegaly, were also significant factors for baseline Hb or anaemia in African patients, especially in children < 5 years [21], in children and adults in Thailand [26], and in adults in Indonesia [45].We found no association between baseline Hb and baseline parasite count in contrast to WWARN who reported an inverse relationship until a parasitaemia of 10,000/μL, followed by a positive relationship [22], whilst Zwang et al. reported a weak positive relationship between baseline Hb and baseline parasitaemia in children < 5 years (mean 0.08 g/dL/year) but an inverse relationship in patients aged ≥ 5 years [21].
The majority of patients with uncomplicated falciparum infections experience an initial fall in Hb and reach their nadir concentration on day 2 [21,22] and the first posttreatment week is the time of highest risk for a blood transfusion [23-25, 36, 55].Therefore, the fractional decline and the nadir Hbs are important parameters to consider in conjunction with clinical signs of severity before deciding the need for a blood transfusion.Previous work has also identified hyperparasitaemia (≥ 150,000/μL) [23] and giving CDA or CD to G6PDd males and heterozygous females [24,25] as risk factors for a blood transfusion.
We found that a lower nadir Hb was associated with a lower baseline Hb in younger children with longer illnesses and higher baseline parasite counts; these factors were also associated with greater fractional falls in Hb except that a higher baseline Hb was associated with a greater fractional fall.Importantly, SLDPQ, G6PDd, and α-thalassaemia genotypes were not explanatory factors for both parameters.HbAS was associated with a slightly higher nadir Hb and lower fractional fall compared to those with normal Hbs that may be related to greater pretreatment haemolysis and, therefore, fewer red cells are available to haemolyse after treatment is given.Our findings concur with previous work showing that younger age, higher baseline Hbs and parasitaemias, and length of illness are significantly associated with a greater mean fractional fall in Hb [21,22,26,27,45], although Zwang et al. found age was not a factor [21].
Few studies have examined the effects of inherited blood disorders on malaria dynamics.Genotypically defined G6PDd and CYP2D6 genotypes were not factors in the initial absolute fall in falciparum-infected Tanzanian patients (median age 6.4 years), who were treated with SLDPQ [56], and α-thalassaemia genotype was not a factor in the D7 fractional fall in Hb in Tanzanian children (median age 4 years) [34].In a small study from Cambodia, there was no association between the fractional fall in Hb and G6PD Viangchan and HbE status [27].
Although we observed a rapid median time to Hb recovery (14 days), ~ 25% of children did not recover their Hb.A higher baseline Hb resulted in a longer time to Hb recovery (Fig. 4), consistent with Zwang et al. [21], a lower probability of achieving Hb recovery, and a lower MAFt (Fig. 4); the latter was also reported in a small study of falciparum-infected Papuan adults [45].By contrast, our multivariable analysis found a higher D0 Hb was associated with a higher, albeit very small, mean increase in MAFt (Table 2), and this discrepancy might be a chance finding.Increasing age, a higher initial fractional fall and being G6PDd were associated with achieving a higher D42 Hb concentration and, in addition to a lower D0 Hb, a greater probability of Hb recovery by D42.G6PDd children had a more rapid time to Hb recovery, only if they had normal baseline Hbs, ~ 70% increased likelihood of D42 Hb recovery compared to G6PD normal children, and achieved a modest increase in mean D42 Hb concentration of 0.2 g/dL.
Despite the apparent disadvantage of a lower baseline Hb, the D42 Hb concentration was not affected, in contrast to a small study of adults from Papua [45].However, a lower D42 Hb concentration was associated with either α-thalassaemia genotype (potentially reflecting normal physiology), sickle cell disease, or failing treatment; whilst the latter has been previously reported in two studies [26,45], it was not a factor in Zwang et al. 's meta-analysis of posttreatment anaemia in patients of all ages [21].
Our study had several limitations.The number of children with moderate anaemia was quite low (~ 7%); there were only three children with HbSS (cautioning interpretation of significant findings), and the majority of children had good nutritional status.As a result, we cannot generalise our findings to children in the community with Hb concentrations of 4-6 g/dL and uncomplicated disease, those with HbSS, or those with moderate or severe acute malnutrition.We dosed SLDPQ by age, which is typically associated with broader mg/kg doses compared to weight-based dosing [35], so our findings may differ from studies based on weight-based regimens.Finally, we conducted many analyses, meaning that some significant findings (e.g. the association between the baseline Hb and the MAFt) might have reflected chance.

Conclusions
To conclude, our large study has shown clearly that SLDPQ did not affect any of the posttreatment parameters of Hb dynamics and has provided additional reassuring evidence of its tolerability.Moreover, G6PDd was only associated with one negative finding-a lower mean baseline Hb concentration.Thereafter, it had a positive posttreatment effect on achieving Hb recovery and a higher D42 Hb but, crucially, did not affect the nadir Hb or the initial fall in Hb.Moreover, children with a Hb < 8 g/dL had a robust posttreatment Hb recovery, allaying anxieties that a lower baseline Hb is necessarily harmful.
Our findings strongly support the notion that using SDLPQ in G6PDd patients is safe and does not increase the risk of subsequent blood transfusion.Indeed, equal numbers of G6PDd and G6PD normal patients were transfused and our overall transfusion rate was just under 1% [36].SLDPQ should be deployed more widely in Africa as part of a global strategy to eliminate ARPf.

Fig. 1 Fig. 2
Fig. 1 Changes in mean haemoglobin concentrations over time by G6PD status

Fig. 3
Fig. 3 Times to haemoglobin recovery by G6PD status as a function of any degree of baseline anaemia (A = not anaemic, B = anaemic)

Table 2
Factors associated with the haemoglobin concentrations at baseline, day of nadir, and on day 42 in 1137 patients and changes over time