https://orcid.org/0000-0002-5734-0845
Figures
Abstract
Background
Glucose-6-phosphate dehydrogenase (G6PD) activity is dependent upon G6PD genotype and age of the red blood cell (RBC) population, with younger RBCs having higher activity. Peripheral parasitemia with Plasmodium spp. induces hemolysis, replacing older RBCs with younger cells with higher G6PD activity. This study aimed to assess whether G6PD activity varies between individuals with and without malaria or a history of malaria.
Methods and findings
Individuals living in the Chittagong Hill Tracts of Bangladesh were enrolled into 3 complementary studies: (i) a prospective, single-arm clinical efficacy trial of patients (n = 175) with uncomplicated malaria done between 2014 and 2015, (ii) a cross-sectional survey done between 2015 and 2016 (n = 999), and (iii) a matched case–control study of aparasitemic individuals with and without a history of malaria done in 2020 (n = 506). G6PD activity was compared between individuals with and without malaria diagnosed by microscopy, rapid diagnostic test (RDT), or polymerase chain reaction (PCR), and in aparasitemic participants with and without a history of malaria. In the cross-sectional survey and clinical trial, 15.5% (182/1,174) of participants had peripheral parasitemia detected by microscopy or RDT, 3.1% (36/1,174) were positive by PCR only, and 81.4% (956/1,174) were aparasitemic. Aparasitemic individuals had significantly lower G6PD activity (median 6.9 U/g Hb, IQR 5.2–8.6) than those with peripheral parasitemia detected by microscopy or RDT (7.9 U/g Hb, IQR 6.6–9.8, p < 0.001), but G6PD activity similar to those with parasitemia detected by PCR alone (submicroscopic parasitemia) (6.1 U/g Hb, IQR 4.8–8.6, p = 0.312). In total, 7.7% (14/182) of patients with malaria had G6PD activity < 70% compared to 25.0% (248/992) of participants with submicroscopic or no parasitemia (odds ratio [OR] 0.25, 95% CI 0.14–0.44, p < 0.001). In the case–control study, the median G6PD activity was 10.3 U/g Hb (IQR 8.8–12.2) in 253 patients with a history of malaria and 10.2 U/g Hb (IQR 8.7–11.8) in 253 individuals without a history of malaria (p = 0.323). The proportion of individuals with G6PD activity < 70% was 11.5% (29/253) in the cases and 15.4% (39/253) in the controls (OR 0.7, 95% CI 0.41–1.23, p = 0.192). Limitations of the study included the non-contemporaneous nature of the clinical trial and cross-sectional survey.
Conclusions
Patients with acute malaria had significantly higher G6PD activity than individuals without malaria, and this could not be accounted for by a protective effect of G6PD deficiency. G6PD-deficient patients with malaria may have higher than expected G6PD enzyme activity and an attenuated risk of primaquine-induced hemolysis compared to the risk when not infected.
Author summary
Why was this study done?
- Primaquine and tafenoquine are the only drugs that prevent relapses of the vivax malaria parasite.
- Primaquine and tafenoquine cause significant side effects in patients with low activity of the G6PD enzyme.
- Approximately 400 million people worldwide are affected by low G6PD enzyme activity.
What did the researchers do and find?
- We measured G6PD activity among individuals with and without malaria.
- We found G6PD activity to be significantly higher among malaria patients.
- We then enrolled healthy individuals with and without a history of malaria to assess whether people with low G6PD activity were less likely to contract malaria.
- G6PD activity did not differ significantly between people with and without a history of malaria who were free of malaria at the time point of G6PD measurement.
What do these findings mean?
- We conclude that G6PD activity increases during a malaria infection and that the observed clinically relevant difference in G6PD activity between those with and without malaria is not due to people with lower G6PD activity being less likely to contract malaria.
- These findings suggest that standard treatment for vivax malaria may be safer than assumed.
- These findings will need to be confirmed.
Citation: Ley B, Alam MS, Kibria MG, Marfurt J, Phru CS, Ami JQ, et al. (2021) Glucose-6-phosphate dehydrogenase activity in individuals with and without malaria: Analysis of clinical trial, cross-sectional and case–control data from Bangladesh. PLoS Med 18(4): e1003576. https://doi.org/10.1371/journal.pmed.1003576
Academic Editor: Kenji Hirayama, Institute of Tropical Medicine (NEKKEN), Nagasaki University, JAPAN
Received: September 14, 2020; Accepted: March 1, 2021; Published: April 23, 2021
Copyright: © 2021 Ley 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 included as supporting information.
Funding: This study was funded by the Bill and Melinda Gates Foundation (OPP1054404) and the Wellcome Trust (Senior Fellowship in Clinical Science awarded to RNP, 200909). SA is supported by a Georgina Sweet Award for Women in Quantitative Biomedical Science and KT is CSL Centenary fellow. No funding bodies had any role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Abbreviations: AMM, adjusted male median; DSS, demographic surveillance system; G6PD, glucose-6-phosphate dehydrogenase; G6PDd, glucose-6-phosphate dehydrogenase deficiency; Hb, hemoglobin; icddr,b, International Centre for Diarrhoeal Disease Research, Bangladesh; OR, odds ratio; PQ, primaquine; RBC, red blood cell; RDT, rapid diagnostic test; TQ, tafenoquine
Introduction
Over the last 2 decades, significant progress has been made in reducing the burden of malaria. However, the impact of intense malaria control activities has been far greater for Plasmodium falciparum than P. vivax, which is now the predominant malaria species in many co-endemic countries [1,2]. Unlike P. falciparum, P. vivax forms dormant liver stages (hypnozoites) that can reactivate weeks to months after the first infection, causing recurrent episodes of febrile illness associated with significant morbidity and mortality [3]. To date, the only drugs that can kill hypnozoites are the 8-aminoquinolines compounds, primaquine (PQ), and tafenoquine (TQ), which, although generally well tolerated, can cause severe hemolysis in individuals with glucose-6-phosphate dehydrogenase (G6PD) deficiency [4]. In view of the risk of severe drug-induced hemolysis, the World Health Organization (WHO) recommends that all patients should be tested for G6PD deficiency (G6PDd) prior to prescribing PQ or TQ [5].
G6PDd is an inherited X-linked enzymopathy with highest prevalence in malaria endemic settings [6,7]. The strong geographic overlap between current and historic areas of malaria transmission and G6PDd suggest an interaction between these exposures; however, investigations of this interaction are conflicting [8]. The correlation between G6PD genotype and phenotype is poor [9–12], since G6PD activity is also determined by a range of factors including the age of the red blood cell (RBC) population, with younger RBCs exhibiting higher levels of G6PD activity than older ones, and, in heterozygous females, the degree of Lyonization [13].
In the presence of oxidative stress, triggered by certain foods, drugs, and infection, G6PD-deficient RBCs hemolyze to a far greater degree than G6PD-normal RBCs [14]. Plasmodia are intraerythrocytic parasites that induce intravascular hemolysis, and this occurs in both parasitized and unparasitized RBCs [15,16]. Whilst P. vivax preferentially invades younger RBCs, hemolysis of non-parasitized cells may be particularly marked in older RBCs with lower G6PD activity [17]. We hypothesize that G6PD activity is elevated in patients with malaria compared to individuals without malaria, that this enzyme elevation is attributable to acute parasitemia, and that this enzyme elevation occurs to a greater degree than that attributable to inherent protection from infection in these individuals. The aim of this study was to determine whether G6PD activity varied in individuals with and without malaria and whether this was attributable to a reduced risk of malaria.
Methods
Malaria in Bangladesh is highly seasonal, with more than 80% of all cases reported during the rainy season from May to October, and 80% to 90% of all reported cases are attributed to P. falciparum, with the majority of other cases due to P. vivax [18].
Individuals were enrolled into 3 studies at the same location in the Chittagong Hill Tracts, Bangladesh: a clinical trial of patients presenting with uncomplicated P. falciparum or P. vivax malaria [19], a cross-sectional survey [20], and a case–control study. First, we determined whether G6PD activity differed between patients with symptomatic malaria, recruited from the efficacy trial, and afebrile and malaria-negative participants from the cross-sectional survey. We then compared G6PD activities between participants with and without a history of malaria, all of whom were afebrile and free of malaria at the time point of G6PD measurement. This study is reported as per the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guideline (S1 File).
Clinical trial
The prospective, single-arm clinical trial [19] was conducted between September 2014 and February 2015 in Alikadam in the district of Bandarban (Chittagong Hill Tracts) and was registered with ClinicalTrials.org as NCT02389374. Patients with uncomplicated malaria diagnosed by microscopy were enrolled and treated according to national treatment guidelines. Venous blood samples were collected for phenotypic and genotypic G6PD testing on day of enrollment, prior to administration of PQ, which was commenced 48 hours later.
Cross-sectional survey
The cross-sectional survey [20,21] was conducted in the area around Alikadam, district of Bandarban, between August 2015 and January 2016 to assess the prevalence of G6PDd in the local population. Venous blood, collected from 1 household member of randomly selected households, was analyzed for peripheral parasitemia and quantified by microscopy and malaria rapid diagnostic test (RDT) (FalciVax, Zephyr Biomedicals, Goa, India). Parasite species was confirmed by PCR [20]. G6PD status was determined phenotypically and genotypically as described below. Participants with microscopic parasitemia provided a sample for laboratory analysis prior to being referred to the closest healthcare facility, where they were treated according to national treatment guidelines.
Case–control study
In 2009, a prospective demographic surveillance system (DSS) was established in the Bandarban district amongst a population of 20,558 and is currently part of the INDEPTH network [22–24]. Coordinates of every household within the surveillance area were recorded, and all households were visited between 2009 and 2010 to enroll participants [25]. All participants were assigned a unique identifier (ID) at the initial visit and then reviewed every 3 months, until December 2014, then every 4 months in 2015, and biannually since 2016. At the end of 2018, a total of 19,117 people were under surveillance. At each follow-up visit, a history of malaria was recorded, along with in- and outmigration, births, deaths, and changes of address. Febrile participants were asked to notify the study clinical staff, who then provided a malaria RDT (FalciVax, Zephyr Biomedicals, Goa, India). RDT-positive participants were prescribed free antimalarial treatment according to national treatment guidelines, the episode of malaria was recorded, and the participant was reviewed at the end of treatment to ensure clinical recovery.
This cohort was used to generate a case–control study comprising individuals who were all afebrile and aparasitemic at the time of enrollment. Cases were selected randomly from DSS participants who had had at least 1 episode of malaria confirmed by microscopy during their monitoring within the DSS. Controls were selected from participants with no recorded episodes of malaria during their monitoring within the DSS. Cases and controls were selected in a ratio of 1:1, and matched according to the duration under surveillance, sex, age, and village of residence.
The recorded history of malaria was confirmed by all participants before enrollment, and the absence of parasitemia was determined by RDT (CareStart Malaria HRP2/pLDH [Pf/Pv] Combo RDT, Access Bio, Somerset, NJ, US) and confirmed by microscopic blood film examination. Thick and thin microscopy films were stained with Giemsa, parasites were counted per 200 white blood cells or 2,000 RBCs by 2 independent readers, and the mean of both readings was recorded. In the case of discrepant results, a third reading was done by an expert microscopist, and this finding was considered final.
After informed written consent, individuals were included in the study if no peripheral parasitemia was detected by either RDT or microscopy, and they had an axillary temperature of less than 37.5°C and no history of fever in the preceding 48 hours. Individuals not fulfilling these criteria were excluded and replaced by participants fitting all criteria until the desired sample size was reached. A venous blood sample was collected from all participants to measure hemoglobin (Hb) concentration, measure G6PD activity by spectrophotometry, and genotyping of the G6PD variant.
Hb measurement and G6PD phenotypic testing
Hb concentration was measured immediately after venipuncture, using a HemoCue 301 (Angelholm, Sweden). Venous samples were stored at 4°C, and a cold chain was maintained until samples reached the reference laboratory at the International Centre for Diarrhoeal Disease Research, Bangladesh (icddr,b) in the capital Dhaka.
G6PD activity was measured using a UV-1800 Shimadzu spectrophotometer (Shimadzu, Kyoto, Japan). Samples from the efficacy trial and the cross-sectional survey were tested using kits and controls (normal and deficient) from Randox (Crumlin, UK) at 37°C [19,26]. Samples from the case–control study were analyzed at 30°C using kits from Pointe Scientific (Canton, MI, US); quality control was done using normal, intermediate, and deficient controls from Analytical Control Systems (Fishers, IN, US) [19,21]. Controls were run daily before testing, which only commenced if the control readings were within the recommended G6PD activity range. All measurements were run in duplicate and repeated if the difference between the first and second measurement was greater than 10%. Since absolute values between different spectrophotometry assays can vary, results from the Randox and Pointe Scientific assays were not compared directly but rather used to categorize participants according to the percentage of normal enzyme activity [27].
G6PD genotypic testing
All samples from patients with G6PD activity less than 70% of the adjusted male median (AMM) [28] and a randomly selected subset of samples with higher G6PD activity were genotyped by PCR at a reference laboratory in Dhaka for the variants Mahidol, Viangchan, Mediterranean, and Orissa, as described previously [29]. All PCR products were also sequenced to identify the Kalyan–Kerala variant [20,30].
Sample size
G6PD data were available from 175 patients with uncomplicated malaria enrolled into the efficacy trial [19] and 999 individuals from the cross-sectional survey [20,26]. The case–control study was powered to determine a minimal difference of 0.75 U/g Hb between cases and controls with a 95% 2-sided confidence interval and 80% power, which required 253 cases and 253 controls. The total sample size of the 3 studies was 1,680 participants.
Statistical analysis
Clinical and laboratory data were entered on paper-based case report forms and digitalized using EpiData version 3.1 (EpiData, Odense, Denmark) [31]. All analyses were done using Stata version 14 (StataCorp, College Station, TX, US).
G6PD activity values, measured by spectrophotometry (in U/dl), were divided by the corresponding Hb measurement (in g/dl) collected in the field to provide a result in U/g Hb. The results from Randox control samples generated from the efficacy trial and cross-sectional survey were compared to assess whether the measurements from the trial and survey were comparable.
The AMM was calculated for each spectrophotometry assay as described previously and defined as 100% activity [28]. Current treatment guidelines recommend that individuals with G6PD activity < 30% be categorized as G6PD deficient, those with G6PD activity from 30% to <70% as G6PD intermediate, and those with activity ≥70% as G6PD normal [5,32,33].
Peripheral parasitemia was defined according to 3 categories: (i) individuals with peripheral parasitemia diagnosed by blood film examination or RDT and confirmed by PCR (microscopic parasitemia), (ii) individuals with submicroscopic parasitemia diagnosed by PCR but not microscopy or RDT (submicroscopic parasitemia), and (iii) individuals without detectable peripheral parasitemia assessed by microscopy, RDT, or PCR (aparasitemia).
Normally distributed continuous data were compared using Student t test or paired Student t test. Non-normally distributed continuous data were compared using the Wilcoxon signed-rank test, the Hodges–Lehmann estimator, or the Kruskal–Wallis equality-of-populations rank test. Unpaired categorical data were compared using the chi squared test or Fisher exact test as appropriate. For each spectrophotometry assay, a multivariable regression model was developed to predict non-normalized G6PD activity (in U/dl), including Hb, sex, and G6PD genotype as biologically plausible independent variables. In addition, malaria status (defined as negative, positive only by PCR, or positive by microscopy or RDT) was included in the model comparing participants with and without parasitemia, while history of malaria was included as a predictor in the model corresponding to the case–control study.
Ethics
All protocols were approved by the Ethics Review Committee of the icddr,b, Bangladesh, and the Human Research Ethics Committee of the Northern Territory Department of Health and Menzies School of Health Research, Australia. Protocol-specific identification numbers together with the exact title of the protocol can be found in S1 Table. Written informed consent was collected from all participants before enrollment, and, in addition, written assent was collected from all minors above the age of 11 years.
Results
Cross-sectional survey and clinical trial
Quality control samples were run on 75 normal and deficient control isolates during the efficacy trial and 255 control isolates during the cross-sectional survey. The median reading from normal controls in the efficacy trial was 8.7 U/g Hb (IQR 8.2–9.2) compared to 8.3 U/g Hb (IQR 7.8–8.7) in the cross-sectional survey (difference 0.4 U/g Hb, 95% CI 0.2–0.6). In deficient controls the median was 1.3 U/g Hb (IQR 1.0–1.7) in the efficacy trial and 1.6 U/g Hb (IQR 1.4–1.8) in the cross-sectional survey (difference −0.3 U/g Hb, 95% CI −0.4 to −0.2). After stratifying all results for type of control, a Hodges–Lehmann estimation showed no significant difference between the controls of the 2 studies (p = 0.962) (Fig 1).
G6PD, glucose-6-phosphate dehydrogenase.
A total of 1,174 participants were enrolled in the clinical trial and cross-sectional survey and assessed by spectrophotometry, of whom 182 (15.5%) had microscopic parasitemia, 36 (3.0%) had submicroscopic parasitemia, and 956 (81.4%) were aparasitemic (Fig 2). Of the 182 participants with microscopic parasitemia, 117 (64.3%) were infected with P. falciparum, 54 (29.7%) with P. vivax, and 11 (6.0%) with a mixed P. falciparum and P. vivax infection. Among the 36 participants with submicroscopic parasitemia, 17 (47.2%) had a P. falciparum infection, 16 (44.4%) had a P. vivax infection, and 3 (8.3%) were positive for both species.
*Four participants had a P. falciparum parasitemia detected by RDT, but not by microscopy. PCR, polymerase chain reaction; RDT, rapid diagnostic test.
Overall, 25.2% (134/531) of males had microscopic parasitemia, compared to 7.5% (48/643) of females (p < 0.001). All subsequent analyses were therefore stratified by sex. Participants in the treatment efficacy study differed significantly from those in the survey: They were younger, had lower body weight, were more likely to be febrile, and had lower Hb concentration (Table 1).
Individuals with microscopic parasitemia had a median G6PD activity of 7.86 U/g Hb (IQR 6.62–9.84), compared to 6.88 U/g Hb (IQR 5.16–8.59) in aparasitemic participants (p < 0.001). Conversely, there was no significant difference in G6PD activity between submicroscopic and aparasitemic individuals (p > 0.05) (Table 2; Figs 3 and S1).
Red vertical lines indicate 30%, 70%, and 100% G6PD activity of the adjusted male median (from left to right). 100% G6PD activity = 7.4 U/g Hb. Participants with submicroscopic parasitemia included in malaria-negative cohort. G6PD, glucose-6-phosphate dehydrogenase.
The AMM was defined as 7.4 U/g Hb (IQR 5.68–9.07), based on 496 normal males. A total of 77.7% (912/1,174) of participants were categorized as G6PD normal (≥70%), 14.3% (168/1,174) had intermediate G6PD activity (30% to <70%), and 8.0% (94/1,174) were G6PD deficient (<30%).
There was no significant difference in G6PD status between aparasitemic and subpatent parasitemic individuals (Table 2), and these groups were therefore pooled (Table 3). Overall, 9.2% (92/992) of those without microscopic parasitemia had G6PD activities below 30%, compared to 1.1% (2/182) of individuals with microscopic parasitemia (odds ratio [OR] 9.2, 95% CI 2.43–77.71, p < 0.001). At the 70% G6PD threshold, 33.3% (244/992) of participants without microscopic parasitemia were G6PD deficient compared to 7.6% (14/182) of patients with microscopic parasitemia (OR = 4.0, 95%CI: 2.26 to 7.61, p<0.001) (Table 2). This difference was particularly apparent in females, of whom 23.0% (137/595) with submicroscopic parasitemia or aparasitemia were intermediate or severely deficient compared to 2.1% (1/48) of those with microscopic parasitemia (OR 14.1, 95% CI 2.4–570.9, p < 0.001). The corresponding proportions for males were 9.7% (13/134) and 28.0% (111/397) (OR 3.6, 95% CI 1.9–7.3, p < 0.001).
In the 178 individuals with microscopic parasitemia, the geometric mean parasitemia was 5,107 (95% CI 3,632–7,180) parasites μl−1. There was no significant correlation between peripheral parasitemia and G6PD activity (rs = 0.001, p = 0.988), and this did not change after stratifying by species (all p > 0.05).
Genotyping was undertaken in 29.6% (348/1,174) of patients, including all 94 patients with G6PD activity < 30%, 79.8% (134/168) of randomly selected patients with intermediate G6PD activity, and 13.2% (120/912) of randomly selected G6PD-normal patients. A known G6PD variant was identified in 31.6% (110/348) of samples; variants showed a considerable range of G6PD activity (Tables 4 and S2).
The median non-normalized G6PD activity was 90.8 U/dl (IQR 72.4–111.5). In a multivariable regression model (F[9, 338] = 24.78, p < 0.001, adjusted R2 = 0.3815) including all 348 participants who were genotyped, the following independent factors were correlated negatively with G6PD activity: G6PD variant, female sex, and Hb concentration. Being positive for either species by microscopy, but not submicroscopic parasite carriage, was correlated positively with G6PD enzyme activity (Table 5).
Case–control study
A total of 506 participants were enrolled in the case–control study, 253 in each arm, which were well matched (Table 6). Before enrollment, cases had been under observation for a median time of 8.7 years (IQR 8.4–8.9, range 1.3–9.5) and controls for a median time of 8.6 years (IQR 8.5–8.9, range 1.9–9.5; p = 0.669).
The median G6PD activity was 10.2 U/g Hb (IQR 8.7–11.8) for cases and 10.3 U/g Hb (IQR 8.8–12.2) for controls (p = 0.323) (Fig 4). The combined AMM for cases and controls was 10.7 U/g Hb (IQR 9.0–12.4). There was no significant difference between cases and controls in the proportion of patients categorized as G6PD normal, intermediate, or deficient (p = 0.328); although cases were at lower risk than controls of having a G6PD deficiency defined at the 30% threshold (OR = 0.60, 95%CI: 0.27 to 1.29) and at the 70% threshold (OR = 0.71, 95%CI: 0.41 to 1.23) this did not reach statistical significance (p = 0.155 and p = 0.241 respectively) (Fig 4; Table 7).
Red vertical lines indicate 30%, 70%, and 100% G6PD activity of the adjusted male median. G6PD, glucose-6-phosphate dehydrogenase.
All but 1 of the G6PD-deficient (<30% activity) participants were genotyped (n = 33), and a known variant was identified in 97.0% (32/33); also, all 34 G6PD-intermediate participants were genotyped, and a variant was identified in 23.5% (8/34) of participants (S3 Table).
The median non-normalized G6PD activity was 147.4 U/dl (IQR 123.0–171.5). In a multivariable regression model (F[7, 60] = 27.91, p < 0.001, adjusted R2 = 0.7377, n = 68), being hemi- or homozygous for the Mahidol or Orissa variant and being heterozygous for the Mediterranean variant were found to be significant predictors of non-normalized G6PD activity, but Hb (p = 0.169), sex (p = 0.693), and a history of malaria (p = 0.725) were not significant predictors (Table 8).
Discussion
Our study highlights a significant difference in G6PD activity between patients with clinical malaria and asymptomatic individuals who were either aparasitemic or had very low-level parasitemia that was only detected by PCR. This difference could not be accounted for by a protective effect against malaria from low G6PD activity [8]. Although individuals with a history of malaria were less likely to be G6PD deficient or G6PD intermediate, the ORs were only 0.6 and 0.7, respectively, and did not reach statistical significance.
Within an individual, G6PD activity varies with the age of the RBC population, with younger RBCs, reticulocytes, and nucleated normoblasts having higher G6PD activity [13,34–36]. Malaria induces hemolysis, to which the host responds through reactive erythropoiesis, increasing the number of reticulocytes and, in severe cases, normoblasts within the bloodstream [36–39]. As a consequence, the mean age of the host RBC population decreases, while the mean G6PD activity increases. It is therefore plausible that G6PD activity rises during or shortly after an episode of acute malaria [40].
The WHO malaria treatment guidelines recommend that patients with malaria who have G6PD enzyme activity less than 30% should not be treated with a standard 14-day regimen of PQ, but rather PQ should be administered weekly for 8 weeks [33]. These guidelines were formulated following studies suggesting that, except for the Mediterranean variant, PQ-induced hemolysis is self-limiting and that an initial dose of PQ would lead to a replacement of older RBCs with younger cells, conferring transitional low-grade resistance against PQ-induced hemolysis [32,33,40–43]. Acute malaria may have a similar effect.
The hemolytic risk associated with PQ rises as G6PD activity falls [43]. G6PD-deficient individuals with activity below 30% can be identified with qualitative diagnostics that can be used in remote settings [44]. Heterozygous females with intermediate G6PD activity (30% to <70%) are also at a significant risk of PQ-induced hemolysis, but their diagnosis requires a quantitative enzyme assay [45]. Our analysis highlights that, compared to females with subpatent parasitemia or aparasitemia, females with microscopic parasitemia were 14 times less likely to be G6PD intermediate or deficient; the corresponding OR for males was 3.6. We hypothesize that this is due to preferential parasite-induced hemolysis of G6PD-deficient RBCs in mosaic heterozygous females. If confirmed, our findings suggest that the risk of PQ-induced hemolysis in heterozygous females with malaria would be considerably less than that predicted by cross-sectional surveys of healthy populations.
Our study has some important limitations. We relied on non-contemporaneous collection of specimens from a prospective, single-arm clinical trial and cross-sectional survey, rather than clinical observation of the same patients over multiple time points with and without malaria. However, when comparing results from laboratory quality control testing, results did not differ significantly, suggesting that the results of the studies are comparable.
In the case–control study, participants were tested using a combination RDT for both P. falciparum and P. vivax, and this is known to have lower sensitivity for identifying patients with P. vivax malaria [46]. It is also possible that some participants had malaria that was not detected due to treatment-seeking outside of the surveillance system. However, patients must pay for services and drugs provided at private clinics, pharmacies, and, to some extent, public hospitals; it is unlikely that many participants would prefer this option over free treatment delivered at home as part of the DSS. A subset of 100 participants from each case and control arm were included in a separate study of the perceptions of malaria [47]. A total of 6 participants from the control arm recalled having had malaria previously, whereas only 1 case could not recall any previous episodes of malaria. While these findings are vulnerable to recall and response bias, they suggest that the majority of participants were categorized correctly.
There was a tendency for individuals without a history of malaria to be G6PD intermediate or deficient, as has been observed by others [8,48–50]. A recent study among more than 17,000 infants and neonates in Kenya found the odds in heterozygous females of being admitted to hospital with a malaria-positive microscopy slide to be lower (OR 0.79, 95% CI 0.68–0.93, p = 0.004) compared to G6PD-normal individuals; a recent article from Afghanistan reported an even greater protective effect against P. vivax malaria for the Mediterranean variant (76% in hemizygous males and 55% in females) [51,52]. In our analysis the protective effect was modest and did not reach statistical significance. If a protective effect of G6PDd against malaria in Bangladesh was present, it is likely to be a minor contributor to G6PD status and would not have accounted for asymptomatic and parasite-free individuals being 4- to 9-fold more likely to be diagnosed as G6PD deficient in the clinical trial and cross-sectional survey.
More than 30% of participants were genotyped against a panel of G6PD variants known to be prevalent in the study area [4,53,54]. However, a genetic variant could not be identified in a large proportion of phenotypically G6PD-deficient individuals, perhaps indicative of the genetic diversity in the area [20]. Only patients with known variants were included in the regression models, and, accordingly, these were restricted to 25% of all study participants, which could have affected their validity. Reassuringly, in the multivariable models, the phenotypic observations highlighted that malaria was a significant predictor of phenotypic G6PD activity, whereas previous history of malaria was not.
In conclusion, our study highlights the potential impact of acute malaria on G6PD status, which may influence treatment allocation of 8-aminoquinoline drugs. If our findings are confirmed prospectively, 8-aminoquinoline-based radical cure may be safer than assumed in patients with clinical malaria, compared to individuals who have submicroscopic parasitemia or who are aparasitemic.
Supporting information
S1 Fig. Distribution of G6PD activity (in U/g Hb) in patients with and without peripheral parasitemia.
Red vertical lines indicate 30%, 70%, and 100% G6PD activity of the adjusted male median. 100% G6PD activity = 7.4 U/g Hb.
https://doi.org/10.1371/journal.pmed.1003576.s001
(TIF)
S1 Table. Protocol names and ethics identifiers.
https://doi.org/10.1371/journal.pmed.1003576.s002
(DOCX)
S2 Table. Distribution of G6PD activity among genotyped participants.
https://doi.org/10.1371/journal.pmed.1003576.s003
(DOCX)
S3 Table. Distribution of G6PD variants among cases and controls and corresponding G6PD activity.
https://doi.org/10.1371/journal.pmed.1003576.s004
(DOCX)
S2 File. Database. Control results of clinical trial and cross-sectional survey.
https://doi.org/10.1371/journal.pmed.1003576.s006
(XLS)
S3 File. Database. Clinical trial and cross-sectional survey.
https://doi.org/10.1371/journal.pmed.1003576.s007
(XLS)
Acknowledgments
We would like to thank all participants who agreed to participate, and all staff involved in making this work happen. icddr,b is grateful to the governments of Bangladesh, Canada, Sweden, and the UK for providing core/unrestricted support.
References
- 1. Battle KE, Lucas TCD, Nguyen M, Howes RE, Nandi AK, Twohig KA, et al. Mapping the global endemicity and clinical burden of Plasmodium vivax, 2000–17: a spatial and temporal modelling study. Lancet. 2019;394(10195):332–43. pmid:31229233
- 2. Price RN, Commons RJ, Battle KE, Thriemer K, Mendis K. Plasmodium vivax in the era of the shrinking P. falciparum map. Trends Parasitol. 2020;36(6):560–70. pmid:32407682
- 3. Dini S, Douglas NM, Poespoprodjo JR, Kenangalem E, Sugiarto P, Plumb ID, et al. The risk of morbidity and mortality following recurrent malaria in Papua, Indonesia: a retrospective cohort study. BMC Med. 2020;18(1):28. pmid:32075649
- 4. Phru C, Kibria M, Thriemer K, Chowdhury M, Jahan N, Aktaruzzaman MM, et al. Case report: a case of primaquine-induced hemoglobinuria in glucose-6-phosphate dehydrogenase deficient malaria patient in southeastern Bangladesh. Am J Trop Med Hyg. 2020;102(1):156–8. pmid:31701865
- 5.
World Health Organization. Guide to G6PD deficiency rapid diagnostic testing to support P. vivax radical cure. Geneva: World Health Organization; 2018.
- 6. Cappellini MD, Fiorelli G. Glucose-6-phosphate dehydrogenase deficiency. Lancet. 2008;371(9606):64–74. pmid:18177777
- 7. Howes RE, Battle KE, Satyagraha AW, Baird JK, Hay SI. G6PD deficiency: global distribution, genetic variants and primaquine therapy. Adv Parasitol. 2013;81:133–201. pmid:23384623
- 8. Mbanefo EC, Ahmed AM, Titouna A, Elmaraezy A, Trang NT, Phuoc Long N, et al. Association of glucose-6-phosphate dehydrogenase deficiency and malaria: a systematic review and meta-analysis. Sci Rep. 2017;7:45963. pmid:28382932
- 9. Robinson KM, Yang W, Haidar CE, Hankins JS, Jay DW, Kornegay N, et al. Concordance between glucose-6-phosphate dehydrogenase (G6PD) genotype and phenotype and rasburicase use in patients with hematologic malignancies. Pharmacogenomics J. 2019;19(3):305–14. pmid:30206300
- 10. Roh ME, Oyet C, Orikiriza P, Wade M, Mwanga-Amumpaire J, Boum Y 2nd, et al. Screening for glucose-6-phosphate dehydrogenase deficiency using three detection methods: a cross-sectional survey in southwestern Uganda. Am J Trop Med Hyg. 2016;95(5):1094–9. pmid:27672207
- 11. Shah SS, Macharia A, Makale J, Uyoga S, Kivinen K, Craik R, et al. Genetic determinants of glucose-6-phosphate dehydrogenase activity in Kenya. BMC Med Genet. 2014;15:93. pmid:25201310
- 12. Johnson MK, Clark TD, Njama-Meya D, Rosenthal PJ, Parikh S. Impact of the method of G6PD deficiency assessment on genetic association studies of malaria susceptibility. PLoS ONE. 2009;4(9):e7246. pmid:19789650
- 13. Chu CS, Bancone G, Nosten F, White NJ, Luzzatto L. Primaquine-induced haemolysis in females heterozygous for G6PD deficiency. Malar J. 2018;17(1):101. pmid:29499733
- 14. Luzzatto L, Arese P. Favism and glucose-6-phosphate dehydrogenase deficiency. N Engl J Med. 2018;378(1):60–71. pmid:29298156
- 15. White NJ. Malaria parasite clearance. Malar J. 2017;16(1):88. pmid:28231817
- 16. Anstey NM, Russell B, Yeo TW, Price RN. The pathophysiology of vivax malaria. Trends Parasitol. 2009;25(5):220–7. pmid:19349210
- 17. Douglas NM, Anstey NM, Buffet PA, Poespoprodjo JR, Yeo TW, White NJ, et al. The anaemia of Plasmodium vivax malaria. Malar J. 2012;11:135. pmid:22540175
- 18. Noe A, Zaman SI, Rahman M, Saha AK, Aktaruzzaman MM, Maude RJ. Mapping the stability of malaria hotspots in Bangladesh from 2013 to 2016. Malar J. 2018;17(1):259. pmid:29996835
- 19. Ley B, Alam MS, Thriemer K, Hossain MS, Kibria MG, Auburn S, et al. G6PD deficiency and antimalarial efficacy for uncomplicated malaria in Bangladesh: a prospective observational study. PLoS ONE. 2016;11(4):e0154015. pmid:27128675
- 20. Ley B, Kibria M, Khan W, Auburn S, Phru C, Jahan N, et al. Wide range of G6PD activities found among ethnic groups of the Chittagong Hill Tracts, Bangladesh. PLoS Negl Trop Dis. 2020;14(9):e0008697. pmid:32925910
- 21. Ley B, Alam MS, O’Donnell JJ, Hossain MS, Kibria MG, Jahan N, et al. A comparison of three quantitative methods to estimate G6PD activity in the Chittagong Hill Tracts, Bangladesh. PLoS ONE. 2017;12(1):e0169930. pmid:28121993
- 22.
Scobie H, Ahmed S, Ram M, Hossain MS, Prue CS, Khyang J, et al. Initial demographic surveillance system report for the icddr,b Bandarban field site. Scientific Report no. 119. Dhaka: International Centre for Diarrhoeal Disease Research, Bangladesh; 2012.
- 23. Khan WA, Sack DA, Ahmed S, Prue CS, Alam MS, Haque R, et al. Mapping hypoendemic, seasonal malaria in rural Bandarban, Bangladesh: a prospective surveillance. Malar J. 2011;10:124. pmid:21569599
- 24. Streatfield PK, Khan WA, Bhuiya A, Alam N, Sie A, Soura AB, et al. Cause-specific mortality in Africa and Asia: evidence from INDEPTH health and demographic surveillance system sites. Glob Health Action. 2014;7:25362. pmid:25377324
- 25. Galagan SR, Prue CS, Khyang J, Khan WA, Ahmed S, Ram M, et al. The practice of jhum cultivation and its relationship to Plasmodium falciparum infection in the Chittagong Hill Districts of Bangladesh. Am J Trop Med Hyg. 2014;91(2):374–83. pmid:24821843
- 26. Alam MS, Kibria MG, Jahan N, Thriemer K, Hossain MS, Douglas NM, et al. Field evaluation of quantitative point of care diagnostics to measure glucose-6-phosphate dehydrogenase activity. PLoS ONE. 2018;13(11):e0206331. pmid:30388146
- 27. Pfeffer DA, Ley B, Howes RE, Adu P, Alam MS, Bansil P, et al. Quantification of glucose-6-phosphate dehydrogenase activity by spectrophotometry: a systematic review and meta-analysis. PLoS Med. 2020;17(5):e1003084. pmid:32407380
- 28. Domingo GJ, Satyagraha AW, Anvikar A, Baird K, Bancone G, Bansil P, et al. G6PD testing in support of treatment and elimination of malaria: recommendations for evaluation of G6PD tests. Malar J. 2013;12:391. pmid:24188096
- 29. Kashmoola MA, Eissa AA, Al-Takay DT, Al-Allawi NA. Molecular characterization of G6PD deficient variants in Nineveh Province, Northwestern Iraq. Indian J Hematol Blood Transfus. 2015;31(1):133–6. pmid:25548459
- 30. Sarker SK, Islam MT, Eckhoff G, Hossain MA, Qadri SK, Muraduzzaman AK, et al. Molecular analysis of glucose-6-phosphate dehydrogenase gene mutations in Bangladeshi individuals. PLoS ONE. 2016;11(11):e0166977. pmid:27880809
- 31. Ley B, Rijal KR, Marfurt J, Adhikari NR, Banjara MR, Shrestha UT, et al. Analysis of erroneous data entries in paper based and electronic data collection. BMC Res Notes. 2019;12(1):537. pmid:31439025
- 32. von Seidlein L, Auburn S, Espino F, Shanks D, Cheng Q, McCarthy J, et al. Review of key knowledge gaps in glucose-6-phosphate dehydrogenase deficiency detection with regard to the safe clinical deployment of 8-aminoquinoline treatment regimens: a workshop report. Malar J. 2013;12:112. pmid:23537118
- 33.
World Health Organization. Guidelines for the treatment of malaria. 3rd edition. Geneva: World Health Organization; 2015.
- 34. Luzzatto L, Nannelli C, Notaro R. Glucose-6-phosphate dehydrogenase deficiency. Hematol Oncol Clin North Am. 2016;30(2):373–93. pmid:27040960
- 35. Minucci A, Giardina B, Zuppi C, Capoluongo E. Glucose-6-phosphate dehydrogenase laboratory assay: how, when, and why? IUBMB Life. 2009;61(1):27–34. pmid:18942156
- 36. Watson J, Taylor WR, Menard D, Kheng S, White NJ. Modelling primaquine-induced haemolysis in G6PD deficiency. Elife. 2017;6:e23061. pmid:28155819
- 37. Haldar K, Mohandas N. Malaria, erythrocytic infection, and anemia. Hematology Am Soc Hematol Educ Program. 2009;2009:87–93.
- 38. Lamikanra AA, Brown D, Potocnik A, Casals-Pascual C, Langhorne J, Roberts DJ. Malarial anemia: of mice and men. Blood. 2007;110(1):18–28. pmid:17341664
- 39. Hillman RS. Characteristics of marrow production and reticulocyte maturation in normal man in response to anemia. J Clin Invest. 1969;48(3):443–53. pmid:5773082
- 40. Beutler E, Dern RJ, Alving AS. The hemolytic effect of primaquine. IV. The relationship of cell age to hemolysis. J Lab Clin Med. 1954;44(3):439–42. pmid:13201853
- 41. Alving AS, Johnson CF, Tarlov AR, Brewer GJ, Kellermeyer RW, Carson PE. Mitigation of the haemolytic effect of primaquine and enhancement of its action against exoerythrocytic forms of the Chesson strain of Plasmodium vivax by intermittent regimens of drug administration: a preliminary report. Bull World Health Organ. 1960;22:621–31. pmid:13793053
- 42. Dern RJ, Beutler E, Alving AS. The hemolytic effect of primaquine. II. The natural course of the hemolytic anemia and the mechanism of its self-limited character. J Lab Clin Med. 1954;44(2):171–6. pmid:13184224
- 43. Ashley EA, Recht J, White NJ. Primaquine: the risks and the benefits. Malar J. 2014;13:418. pmid:25363455
- 44. Ley B, Winasti Satyagraha A, Rahmat H, von Fricken ME, Douglas NM, Pfeffer DA, et al. Performance of the Access Bio/CareStart rapid diagnostic test for the detection of glucose-6-phosphate dehydrogenase deficiency: a systematic review and meta-analysis. PLoS Med. 2019;16(12):e1002992. pmid:31834890
- 45. Chu CS, Bancone G, Moore KA, Win HH, Thitipanawan N, Po C, et al. Haemolysis in G6PD heterozygous females treated with primaquine for Plasmodium vivax malaria: a nested cohort in a trial of radical curative regimens. PLoS Med. 2017;14(2):e1002224. pmid:28170391
- 46. Gunasekera W, Premaratne RG, Weerasena O, Premawansa WS, Handunnetti SM, Fernando SD. Utility of pf/pan RDT for diagnosis in the prevention of re-establishment of malaria in Sri Lanka. Pathog Glob Health. 2018;112(7):360–7. pmid:30355068
- 47. Matin MA, Sarkar NDP, Phru CS, Ley B, Thriemer K, Price RN, et al. Precarity at the margins of malaria control in the Chittagong Hill Tracts in Bangladesh: a mixed-methods study. Pathogens. 2020;9(10):840. pmid:33066621
- 48. Ruwende C, Khoo SC, Snow RW, Yates SN, Kwiatkowski D, Gupta S, et al. Natural selection of hemi- and heterozygotes for G6PD deficiency in Africa by resistance to severe malaria. Nature. 1995;376(6537):246–9. pmid:7617034
- 49. Guindo A, Fairhurst RM, Doumbo OK, Wellems TE, Diallo DA. X-linked G6PD deficiency protects hemizygous males but not heterozygous females against severe malaria. PLoS Med. 2007;4(3):e66. pmid:17355169
- 50. Clarke GM, Rockett K, Kivinen K, Hubbart C, Jeffreys AE, Rowlands K, et al. Characterisation of the opposing effects of G6PD deficiency on cerebral malaria and severe malarial anaemia. Elife. 2017;6:e15085. pmid:28067620
- 51. Awab GR, Aaram F, Jamornthanyawat N, Suwannasin K, Pagornrat W, Watson JA, et al. Protective effect of Mediterranean-type glucose-6-phosphate dehydrogenase deficiency against Plasmodium vivax malaria. Elife. 2021;10:e62448. pmid:33543710
- 52. Uyoga S, Macharia AW, Ndila CM, Nyutu G, Shebe M, Awuondo KO, et al. Glucose-6-phosphate dehydrogenase deficiency and susceptibility to childhood diseases in Kilifi, Kenya. Blood Adv. 2020;4(23):5942–50. pmid:33275767
- 53. Howes RE, Dewi M, Piel FB, Monteiro WM, Battle KE, Messina JP, et al. Spatial distribution of G6PD deficiency variants across malaria-endemic regions. Malar J. 2013;12:418. pmid:24228846
- 54. Plewes K, Soontarawirat I, Ghose A, Bancone G, Kingston HW, Herdman MT, et al. Genotypic and phenotypic characterization of G6PD deficiency in Bengali adults with severe and uncomplicated malaria. Malar J. 2017;16(1):134. pmid:28356147