Genotypes and phenotypes of G6PD deficiency among Indonesian females across diagnostic thresholds of G6PD activity guiding safe primaquine therapy of latent malaria

Abstract

Background

Plasmodium vivax occurs as a latent infection of liver and a patent infection of red blood cells. Radical cure requires both blood schizontocidal and hypnozoitocidal chemotherapies. The hypnozoitocidal therapies available are primaquine and tafenoquine, 8-aminoquinoline drugs that can provoke threatening acute hemolytic anemia in patients having an X-linked G6PD-deficiency. Heterozygous females may screen as G6PD-normal prior to radical cure and go on to experience hemolytic crisis.

Methods & findings

This study examined G6PD phenotypes in 1928 female subjects living in malarious Sumba Island in eastern Indonesia to ascertain the prevalence of females vulnerable to diagnostic misclassification as G6PD-normal. All 367 (19%) females having <80% G6PD normal activity were genotyped. Among those, 103 (28%) were G6PD wild type, 251 (68·4%) were heterozygous, three (0·8%) were compound heterozygotes, and ten (2·7%) were homozygous deficient. The variants Vanua Lava, Viangchan, Coimbra, Chatham, and Kaiping occurred among them. Below the 70% of normal G6PD activity threshold, just 18 (8%) were G6PD-normal and 214 (92%) were G6PD-deficient. Among the 31 females with <30% G6PD normal activity were all ten homozygotes, all three compound heterozygotes, and just 18 were heterozygotes (7% of those).

Conclusions

In this population, most G6PD heterozygosity in females occurred between 30% and 70% of normal (69·3%; 183/264). The prevalence of females at risk of G6PD misclassification as normal by qualitative screening was 9·5% (183/1928). Qualitative G6PD screening prior to 8-aminoquinoline therapies against P. vivax may leave one in ten females at risk of hemolytic crisis, which may be remedied by point-of-care quantitative tests.

Introduction

Acute malaria caused by Plasmodium vivax asexual blood stage parasites provokes a debilitating febrile illness that may progress to severe and life-threatening disease syndromes associated with death [1,2]. More than 10 million of these attacks occur annually among the 2.8 billion people living at risk [3]. Each infection requires effective blood schizontocidal therapy to arrest progression of disease [4]. Unlike Plasmodium falciparum, the other major cause of human malaria, infection by P. vivax includes latency where dormant hepatic hypnozoites later awaken and provoke renewed attacks (relapses), each one threatening progression to severe disease and onward transmission. In cohorts living in endemic areas of Thailand and Papua New Guinea, hypnozoite-borne infections of blood accounted for >80% of incident P. vivax parasitemias [5,6]. Although the risk and timing of relapse varies across endemic zones [7], at least three and often five or more relapses per infectious event may be the rule [8]. In a cohort of 2,495 American soldiers repatriated from the Pacific theater of World War II, the median number of relapses over two years was 10–14 [9]. A study in eastern Indonesia followed tens of thousands of patients diagnosed with P. falciparum or P. vivax over a decade; while risk of death within 14 days was higher for P. falciparum, risk of multiple attacks, hospitalization, and premature death among P. vivax patients was two-fold higher [10]. Effective treatment of P. vivax malaria that includes therapy against latency would thus provide very substantial clinical and public health benefits [11].

In contrast to the many therapeutic options for arresting the acute attack, the 8-aminoquinoline drugs primaquine or tafenoquine remain the only hypnozoitocidal options for terminating latency of P. vivax malaria. At therapeutic hypnozoitocidal doses, primaquine invariably provokes a potentially life-threatening acute hemolytic anemia in patients having glucose-6-phosphate dehydrogenase (G6PD) deficiency, an X-linked trait affecting over 400 million people at a rate of about 8% where endemic malaria occurs [12,13]. Many dozens of specific single nucleotide polymorphisms occur all across the 13 exons of the G6PD gene [14]. Among male hemizygotes and female homozygotes, G6PD impairment almost invariably occurs at <30% of normal enzyme activity and they are considered the most vulnerable to hemolysis caused by primaquine anti-relapse therapy [15]. The World Health Organization (WHO) recognized four classes of G6PD enzyme based on its activities: Class I as the most severe, Class II as severe, Class III as intermediate and Class IV as normal. Most variants are either class II or III.

Qualitative laboratory and point-of-care G6PD screening tests detect hemizygotes and homozygotes with nearly 100% sensitivity and specificity [16–19]. Female heterozygotes, however, may express G6PD activity phenotypes ranging from fully deficient to completely normal as a consequence of mosaicism of their red blood cell populations [20]. Apparently random inactivation of one X-chromosome or the other yields that heterogeneity of phenotype (a process called Lyonization that occurs during early embryonic development). Qualitative screening cannot reliably differentiate G6PD normal patients from those having >30% to <80% of normal enzyme activity [21,22]. This problem imposes uncertainty regarding a diagnostic reading of “normal” (and eligible for primaquine therapy) by those tests among females. Patients having G6PD activities between 30–80% are under-studied, but there have been reports of a drop in hematocrit up to 30% requiring transfusion therapy following high dose of primaquine, some of those experiencing hemolytic crises as severe as homozygous females [23,24].

In the current study we aimed to evaluate prevalent G6PD genotypes in relation to G6PD activity phenotype in a large sample of females living in a malaria-endemic area of Sumba Island in eastern Indonesia where we have conducted population surveys of G6PD heterogeneity [25,26]. This evidence directly informs the distribution of risk with G6PD screening in the context of primaquine therapy against latent malaria, and may serve to mitigate risk to G6PD heterozygotes. Reversing the historic neglect of this serious clinical and public health problem of poor access to hypnozoitocidal therapy should benefit males and females alike [27], and achieving that requires specific attention to the complexity of the heterozygous state of G6PD deficiency.

Methods

Ethics statement

A protocol detailing this survey was reviewed and approved by Eijkman Institute Research Ethics Committee (17th March 2015, EIREC Approval No. 81). Written informed consent was obtained from each participant and written informed consent was signed by guardian or parent of child participants.

Study population

A cross-sectional prevalence survey of G6PD deficiency, anemia, and malaria status was undertaken in April 2015 at two sub-districts of Southwest Sumba: villages of Wainyapu and Waiha of Kodi Balaghar sub-district and Umbu Ngedo of Kodi Bangedo sub-district. In that year, these locales had respective annual parasite incidences of 58·16, 85·24 and 26·85/cases per 1000 residents respectively according to local health district records (P.falciparum and P. vivax). Only healthy (without fever > 37.5°C or history of fever within 48 hours) females older or equal to 6 years old who provided informed individual or parental/guardian consents were included in the survey. Based on our previous studies in Sumba Island, the prevalence of homozygous and heterozygous G6PD deficient females was 0.19 and thus to get the sample size we applied this equation:

Z is the statistic corresponding to the level of confidence of 1.96 (95% confidence level) and d for precision (2%), our sample size was thus calculated:

This value was 1345, and we increased that by 30% as a means of ensuring and adequate number heterozygous females for our research purposes, yielding a target sample size of 1921.

Venipuncture and malaria screening

Willing subjects who signed written informed consent were asked to allow taking 3 ml of venous blood collected in EDTA tubes. Residual blood in the syringe was used for microscopic examinations for malaria parasites (thick and thin blood smears) in the field laboratory according to standard protocol by technicians certified as competent according to WHO specifications in doing so. Participants found to be positive for malaria infection were offered treatment with therapies as stipulated by the Indonesian Ministry of Health guidance the next day. Blood in sealed EDTA tubes was put on ice immediately after venipuncture and stored at 4°C within a few hours. Within 3 days these samples were transported on cold packs by air to Jakarta for the laboratory analyses detailed here.

G6PD qualitative and quantitative tests

G6PD qualitative testing was done at point-of-care using the Carestart G6PD rapid diagnostic test (CSG, Accessbio, USA) according to the manufacturer’s protocol. In brief, blood placed into test cassette with liquid reagents migrates across a white cellulose wick within a maximum of minutes: G6PD-normal blood causes a distinct purple color to develop, whereas with G6PD-deficient blood either remains white or develops only a very slight purple hue. Technicians performing the test were instructed and trained to classify no color or only lightly colored tests as G6PD-deficient. Refresher training for this took place for 2 days prior to going to the field. Most CSG were read by single technician. A second technician was conferred when the result was found difficult to interpret.

Blood samples in Jakarta were held at 4°C and analyzed less than 24hr arrival from Sumba. These were first examined for hemoglobin (HemoCue Hb 201, Sweden) according to the manufacturer’s instruction. Samples having hemoglobin measurements of less than 10 g/dL were excluded from further analyses in order to avoid the impact of anemia causing falsely elevated G6PD readings [22].

G6PD quantitative test from Trinity Biotech (Cat # 345-B) was used as the reference test for G6PD activity screening. The test relies on the principle of NADPH formation that is proportional to the G6PD activity directly measured using a spectrophotometer (Shimadzu UV-1800 series) at 340 nm absorbance, precisely as prescribed by the test manufacturer. Briefly, 10 μl of blood was added to a snap cap tube containing 1.0 ml of G-6-PDH Assay Solution, mixed thoroughly and incubated at 30°C for 5–10 min before adding 2.0 ml of G-6-PDH Substrate Solution. The tube and was mixed gently by several inversions. The contents were transferred to a cuvette and incubated in the temperature-controlled spectrophotometer at 30°C for five min. The cuvette was then read at 340 nm using water as blank and marked as INITIAL A. The reaction was incubated for another five min at 30°C where the contents were read again and marked as FINAL A. G6PD activity was then calculated as U/g Hb by subtracting FINAL A with INITIAL A and divided by five as the manufacturer’s manual instructs. The reading was done once per sample. The absolute value for normal G6PD activity (as 100%) was calculated as the mean of G6PD activity among subjects having ≥5U/g Hb.

G6PD genotyping

G6PD genotyping was restricted to samples showing <80% of normal G6PD activity. DNA was extracted from EDTA blood samples using a modified salting out method as previously described [28]. Extracted DNA was amplified in PCR and digested with specific restriction enzymes for common variants (Vanua Lava, Viangchan, Chatham, Coimbra and Kaiping) in these regions [25,26]. Samples showing no digestion by those enzymes were whole-gene or next-generation sequenced using primers described by Saunders et al [29].

Statistical analysis

We described the female sample characteristics including age median and range, hemoglobin measurements, G6PD point-of-care screening results (G6PD deficient, G6PD normal), and malaria microscopic findings (negative/positive for malaria, Plasmodium species). The primary outcome was the proportion of G6PD deficient females in the female sample population. We analyzed G6PD enzyme activities using quantitative G6PD Trinity Biotech for the samples with hemoglobin values > 10 g/dl. Mean, median, standard deviations and range of G6PD enzymatic activities were calculated to determine the reference values in normal and deficient subjects. We assessed 100% G6PD activity as the mean of the G6PD activity ≥5U/g Hb as our cut off, and set our 30%, 70% and 80% diagnostic thresholds based on this value. The 5 U/g Hb cut-off was considered the lowest limit of truly normal G6PD activity and protected the estimate of 100% enzyme activity from the diminishing bias of G6PD-deficiency. The diagnostic thresholds reflected those considered effective for CSG (<30%) or safe for administration of 8-aminoquinolines (>70% or >80%, depending on authoritative recommendations) [25,26] We classified genotype characteristics (wild type G6PD genotype, heterozygotes, compound heterozygotes, and homozygotes deficient) for those females with G6PD deficiency. We measured the performance of CareStart G6PD rapid test specific to genotypes and G6PD activity ranges among female subjects below 80% of normal. We also measured sensitivity, specificity, negative and positive predicted value of CSG against spectrophotometric reference test at 10%, 30%, 70% and 80%. A 95% confidence of each indicator was measured. Statistical significance of G6PD prevalence was evaluated by Fisher’s exact test. Scattered boxplot was used to show how the G6PD activities are within G6PD variants identified in the study. Data were analyzed using R software.

Results

Field survey

A total of 2056 females were screened. Table 1 lists their demographic characteristics along with G6PD point-of-care screening, hemoglobin measurements, and malaria blood film examination findings. The prevalence of microscopy positivity for malaria was 6·7% (129/1928), with P. falciparum being the dominant species by a margin of 2:1. P. vivax was found in 53 females, including 11 of 12 mixed species infections (prevalence of 2·6%). A total of 1928 females (93·8%) had hemoglobin levels ≥10g/dL and were included in the study (Fig 1). The point-of-care G6PD rapid test identified 1814 women as normal and 113 women as deficient (1 sample did not have CSG data), indicating a prevalence of 6·2% for G6PD deficiency using that device.

[Figure omitted. See PDF.]

Fig 1. Survey process and sample analysis.

https://doi.org/10.1371/journal.pntd.0009610.g001

[Figure omitted. See PDF.]

Table 1. Survey sample demographics and screening findings.

https://doi.org/10.1371/journal.pntd.0009610.t001

Laboratory findings

Quantitative G6PD testing was accomplished for 1928 samples, with all 128 exclusions being due to hemoglobin values below 10g/dl (Fig 1). Fig 2 illustrates the frequency distribution of subjects across the range of observed G6PD activity values. The mean value among those ≥5 U/g Hb was 11·04 U/g Hb and set as 100% of normal enzyme activity. The thresholds of 10%, 30%, 70%, and 80% are illustrated in Fig 2. The distribution of malaria-positive subjects by species of diagnosis, also in Fig 2, appeared independent of G6PD activity level.

[Figure omitted. See PDF.]

Fig 2. G6PD activity distribution and malaria microscopy positivity in the survey.

https://doi.org/10.1371/journal.pntd.0009610.g002

Fig 3 illustrates G6PD measurements among the 367 subjects having <80% of G6PD normal and according to genotype. Vanua Lava and Viangchan variants co-dominated, with Chatham, Coimbra, and Kaiping variants in minority representation. The very few homozygotes (10) were represented by Vanua Lava, Viangchan, and Coimbra, and all occurred at the lowest spectrum of G6PD activity. The 3 compound heterozygotes represented by Viangchan/Vanua Lava and Chatham/Vanua Lava also occurred at the lower end of G6PD activity (<2 U/g Hb). The majority of G6PD deficiency among females occurred as heterozygotes having between 30% and 70% of normal activity, with few exceptions below the lower threshold and relatively more above the higher threshold; whereas 103 G6PD-normal subjects (Normal and ND) occurred below the 80% threshold, just 18 did so below the 70% threshold. The ND (not determined) lane of Fig 3 represents subjects with <80% of normal G6PD activity negative for RFLP analysis for the SNPs analyzed but who were not whole gene sequenced. Excepting the two subjects below 5 U/g Hb, ND subjects were presumed G6PD-normal.

[Figure omitted. See PDF.]

Fig 3. Genotypes identified among the subjects having <80% of normal G6PD activity and that phenotype illustrated in relation to 30% and 70% of normal activity thresholds.

CO, VC, and VL Hom are Coimbra, Viangchan, and Vanua Lava homozygous respectively. CO, VC, CT and VL Het are Coimbra, Viangchan, Chatham and Vanua Lava heterozygous respectively. VL & CT and VL & VC represent the 3 compound heterozygotes of those variants. The ND lane represent subjects negative for the G6PD variant SNPs analyzed, but not whole gene sequenced, whereas the Normal lane represents those whole gene sequenced as wild type G6PD. The black dots represent outliers calculated from the boxplot whereas the pink outliers are calculated from the scatterplot of same samples.

https://doi.org/10.1371/journal.pntd.0009610.g003

The activity values among heterozygotes appeared normally distributed. Among the 251 heterozygotes detected in this survey, 73% (183) occurred between 30% and 70% of G6PD normal activity, within ±1 standard deviation of the 50% mean predicted for a random distribution between 0 and 100%. Fig 4 illustrates the distribution of G6PD genotypes between 0 and 70% of normal G6PD activity. The relative predominance of G6PD-normal subjects (Normal and ND) above the 70% threshold and paucity of the same below that threshold may be seen. Heterozygotes overwhelmingly dominated the 30% to 70% range. Likewise, below the 30% threshold, homozygotes and compound heterozygotes dominated the lowest values therein, with most of the heterozygotes just below 30%.

[Figure omitted. See PDF.]

Fig 4. Distribution of G6PD genotypes among phenotypic activity thresholds.

ND-not determined, CO–Coimbra, VC–Viangchan, VL–Vanua Lava, CT–Chatham, and KP–Kaiping. Hom–homozygous mutant and Het–heterozygous. VL & CT, and VL & VC represent the compound heterozygotes.

https://doi.org/10.1371/journal.pntd.0009610.g004

Diagnostic result of CareStart G6PD rapid test

Table 2 summarizes the essential features of diagnostic result of the CareStart G6PD rapid test (CSG) relative to the defining Trinity quantitative spectrophotometric assay and genotype result among the 1928 non-anemic subjects (>10 gHb/dL) at varied thresholds of G6PD activity (% of normal). The CSG detected only 33 of 183 (18·0%) heterozygotes between 30% and 70% of normal G6PD activity as G6PD-deficient. The overall sensitivity, specificity, positive predictive value, and negative predictive value of the CSG relative to quantitative testing at the 70% of normal activity threshold was 25.1% (95%CI, 19.7–31.2), 96.8% (95%CI, 95.8–97.6), 21.4% (95%CI, 16.2–27.7), and 97.4% (95%CI, 97.2–97.5), respectively.

[Figure omitted. See PDF.]

Table 2. Summary of diagnostic performance of the CSG specific to genotypes and G6PD activity ranges among individual subjects below 80% of normal.

https://doi.org/10.1371/journal.pntd.0009610.t002

Discussion

The findings of this survey in malarious Sumba Island in eastern Indonesia affirms a G6PD activity threshold of 70% of normal in identifying most G6PD-deficient females for exclusion from 8-aminoquinoline therapies with minimal exclusion of G6PD-normal females. Among the 232 subjects below that threshold, 214 (92%) had genotype-confirmed G6PD deficiency. At a threshold of 80% of normal activity 264 out of 367 (72%) subjects were G6PD-deficient. The 50 G6PD-deficient subjects having <80% but >70% represented 19% of G6PD-deficient females present in the sample. Given the relatively high degree of enzyme activity among those females, which reflect a high proportion of G6PD-normal red blood cells, they presumably represented the least vulnerable to 8-aminoquinoline hemolytic crisis among G6PD-deficient females. While an 80% threshold would exclude those 50 potentially susceptible females from therapy, it would also have excluded 103 G6PD-normal women representing 5% of the population as a whole. The unnecessary exclusion of G6PD-normal women from 8-aminoquinoline therapy is reduced to 0.9% of the population at a threshold of 70%. The actual safety threshold of G6PD activity for anti-relapse therapy with primaquine or tafenoquine has not been determined in practice or clinical research on heterozygotes. Any threshold given remains a supposition and our findings do not offer the assurance of safety at a 70% activity threshold. Instead, our findings only inform the distribution of G6PD heterozygotes at these thresholds of diagnostic performance (<30%) or those of safety (70% or 80%).

As elsewhere in much of the malaria endemic world, most acute malaria in eastern Indonesia occurs at the under-resourced periphery of healthcare delivery. That setting may not accommodate sustainable quantitative G6PD assessments allowing the application of a 70% threshold defining 8-aminoquinoline therapy proceed vs. do not proceed rules. We assessed the likely impacts of qualitative vs. quantitative G6PD assessments in the Sumba population of females surveyed. As reported by others and described from our laboratory [30,31], the CSG qualitative test becomes increasingly insensitive to deficiency at G6PD activities above 30% as it is not designed to distinguish 30–80% from >80% individuals. The sensitivity, specificity, negative and positive predicted values at cut off 70% were 25.1%, 96.8%, 21.4% and 97.4% respectively. The sensitivity and positive predicted values further dropped at 80% cut off (S1 Table and S1 Fig). The current survey showed that 82% of the 183 G6PD-deficient females having G6PD activities between 30% and 70% screened as G6PD-normal by the CSG. In other words, 9·5% (183/1928) of all females at the study sites on Sumba would be at high risk of being cleared for potentially dangerous 8-aminoquinoline therapies with qualitative G6PD screening.

Varied national drug regulators have registered tafenoquine as a single-dose alternative to the standard 14 days daily dosing with primaquine for presumptive anti-relapse therapy (PART) following a diagnosis of P. vivax malaria [32,33]. The administration of tafenoquine is not recommended in patients having <70% of normal G6PD activity; thus, imposing quantitative testing as a necessity corroborated by our findings from females residing in Sumba. Indeed, a Kenyan female A- (Class III) heterozygote accidently dosed with 400mg tafenoquine required transfusion therapy [34]. Among healthy ethnic Karin women in Thailand heterozygous for the moderate Class III Mahidol variant of G6PD deficiency and having between 40% and 60% of normal G6PD activity, all three subjects receiving the 300mg single dose of tafenoquine hemolyzed approximately 8% of their red blood cells [35]. All of the five variants identified at Sumba (Vanua Lava, Viangchan, Chatham, Coimbra, and Kaiping) were Class II and perhaps more vulnerable to 8-aminoquinoline toxicity than Mahidol Class III variant.

Qualitative G6PD screening prior to PART applying low-dose primaquine (0·25mg/kg/d x 14d) may be viewed as adequately safe for females because the prolonged dosing allows mitigation of potential harm by cessation of therapy after the onset of symptoms of acute hemolytic anemia. However, that onset occurs late and abruptly; typically, a day after the 3rd or 4th dose with a sudden drop of hemoglobin accompanied by haemoglobinuria and jaundice, i.e., with hemolytic crisis in progress [36,37]. There is no validated means of monitoring for those events in averting serious harm, i.e., specification of any specific clinical or laboratory parameter at a designated time point during treatment reliably indicating cessation of treatment before serious harm is done. In eastern Indonesia severe Class II variants are both dominant among G6PD-deficient people and highly prevalent in the general population [38]. Dosing G6PD unknowns, especially when using the more efficacious high-dose regimens of primaquine (0·5mg/kg/d x 14d or 1·0mg/kg/d x 7d), cannot be recommended at sites like Sumba without quantitative test screening.

Primaquine efficacy is further complicated by both impacts of partner blood schizontocides and natural polymorphisms of cytochrome P450 2D6 impairing the necessary metabolic processing of primaquine [39–42]. Inadequately treated P. vivax latency in Sumba—where most people at risk of malaria live in impoverished rural settings with chronic malnourishment and co-endemic neglected tropical infections–likely exerts the same insidiously harmful effects observed elsewhere in eastern Indonesia [43]. Furthermore, studies [44,45] have shown the increased risk of anemia due to relapses in no or low dose (0·25 mg/kg) of primaquine. Access to efficacious primaquine therapy against latent malaria at sites like Sumba will require either clinically validated cease-dosing criteria or quantitative G6PD screening. Currently, there are several point-of-care quantitative G6PD screening tests available, but none has yet been validated as reliable and practical for lay users in village settings. In Indonesia during 2021, we found one such commercially available kit costs USD 824 for the instrument, USD12 per test strip and USD 104 for controls, with promising performance that may be used in these instances in lieu of reference test [46,47]. However, the high cost per strip and 1 year expiration date may hinder its wide use in Indonesia.

Primaquine has been in clinical use for nearly 70 years despite of its hemolytic toxicity in G6PD-deficient patients. Global health authorities and national malaria control programmes alike advise using primaquine to prevent relapse of latent malaria applying varying doses, dosing strategies, and G6PD precautions [48]. As the findings reported here illustrate, the therapeutic dilemma of 8-aminoquinoline anti-relapse therapy is complex and involves significant proportions of local populations being vulnerable to potential harm. Protecting the vulnerable minority by withholding therapy exposes all infected to the harm of repeated acute attacks. The health authorities in eastern Indonesia, as in many endemic areas, cannot confidently attack the latent reservoir of that harm because rational and practical guidance for safely and effectively doing so is lacking. The findings reported here aimed to aid in formulating that guidance with specific respect to G6PD screening parameters that help to mitigate risk of hemolysis by 8-aminoquinolines. Practical guidance for efficacious radical cure without quantitative G6PD testing may require validated parameters of clinical monitoring for acute hemolytic anemia.

Study limitation

This study is limited in the biochemical aspect since it did not compare the Kcat of G6PD from genotypically G6PD normal and heterozygous females with G6PD activities below 80%. This knowledge on top of G6PD genotype and activity can help fine tune the hemolytic risk of G6PD heterozygous females in oxidative stress.

Supporting information

S1 Table. Performance of CSG G6PD Test against reference test Pointe Scientific quantitative assay.

https://doi.org/10.1371/journal.pntd.0009610.s001

(DOCX)

S1 Fig. CSG G6PD results relative to G6PD activities derived from the reference test. Red boxes are deficient results by CSG and blue boxes are normal results by CSG.

https://doi.org/10.1371/journal.pntd.0009610.s002

(TIF)

Acknowledgments

We thank the participating subjects at villages Wainyapu, Waiha and Umbu Ngedo. The authors thank Dedi Sudiana as field coordinator, Sunardi as field microscopist, Yustus Selan as the Southwest Sumba Regency Malaria Program Officer, Lisa Mone and Nathan as Malaria Program Officers at Panenggo Ede Health Center who helped during field sampling. The authors thank Rosalie Elvira, Jeny, Lia Waslia and Denny Feriandika for their technical expertise in the lab/field and Fitria Wulandari for her excellent administrative help. The authors also thank Dr. Bob Taylor and Prof. Lucio Luzzatto for their thorough critical review of this manuscript.

Citation: Satyagraha AW, Sadhewa A, Panggalo LV, Subekti D, Elyazar I, Soebianto S, et al. (2021) Genotypes and phenotypes of G6PD deficiency among Indonesian females across diagnostic thresholds of G6PD activity guiding safe primaquine therapy of latent malaria. PLoS Negl Trop Dis 15(7): e0009610. https://doi.org/10.1371/journal.pntd.0009610

References

1. Anstey NM, Douglas NM, Poespoprodjo JR, Price RN. Plasmodium vivax: clinical spectrum, risk factors and pathogenesis. Adv Parasitol. 2012;80:151–201. pmid:23199488

2. Baird JK. Evidence and Implications of Mortality Associated with Acute Plasmodium vivax Malaria. Clin Microbiol Rev. 2013 Jan;26(1):36–57. pmid:23297258

3. 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. The Lancet. 2019 Jul 27;394(10195):332–43. pmid:31229233

4. Baird JK, Valecha N, Duparc S, White NJ, Price RN. Diagnosis and Treatment of Plasmodium vivax Malaria. Am J Trop Med Hyg. 2016 Dec 28;95(6 Suppl):35–51. pmid:27708191

5. Adekunle AI, Pinkevych M, McGready R, Luxemburger C, White LJ, Nosten F, et al. Modeling the Dynamics of Plasmodium vivax Infection and Hypnozoite Reactivation In Vivo. PLoS Negl Trop Dis. 2015 Mar 17;9(3):e0003595. pmid:25780913

6. Robinson LJ, Wampfler R, Betuela I, Karl S, White MT, Suen CSNLW, et al. Strategies for Understanding and Reducing the Plasmodium vivax and Plasmodium ovale Hypnozoite Reservoir in Papua New Guinean Children: A Randomised Placebo-Controlled Trial and Mathematical Model. PLOS Med. 2015 Oct 27;12(10):e1001891. pmid:26505753

7. Battle KE, Karhunen MS, Bhatt S, Gething PW, Howes RE, Golding N, et al. Geographical variation in Plasmodium vivax relapse. Malar J. 2014 Apr 15;13:144. pmid:24731298

8. White NJ, Imwong M. Relapse. Adv Parasitol. 2012;80:113–50. pmid:23199487

9. Kitchen SF. Chapters 40–43. In: Boyd MF (ed), Malariology: a comprehensive survey of all aspects of this group of diseases from a global standpoint. Philadelphia, PA: W.B. Saunders Company; 1949. p. 966–1045.

10. 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 Feb 20;18(1):28. pmid:32075649

11. World Health Organization. Control and elimination of Plasmodium vivax malaria: a technical brief. Geneva: World Health Organization; 2015.

12. Ashley EA, Recht J, White NJ. Primaquine: the risks and the benefits. Malar J. 2014;13:418. pmid:25363455

13. Howes RE, Piel FB, Patil AP, Nyangiri OA, Gething PW, Dewi M, et al. G6PD Deficiency Prevalence and Estimates of Affected Populations in Malaria Endemic Countries: A Geostatistical Model-Based Map. PLOS Med. 2012 Nov 13;9(11):e1001339. pmid:23152723

14. Luzzatto L, Nannelli C, Notaro R. Glucose-6-Phosphate Dehydrogenase Deficiency. Hematol Oncol Clin North Am. 2016 Apr;30(2):373–93. pmid:27040960

15. Luzzatto L, Seneca E. G6PD deficiency: a classic example of pharmacogenetics with on-going clinical implications. Br J Haematol. 2014 Feb 1;164(4):469–80. pmid:24372186

16. Wojnarski B, Lon C, Sea D, Sok S, Sriwichai S, Chann S, et al. Evaluation of the CareStartTM glucose-6-phosphate dehydrogenase (G6PD) rapid diagnostic test in the field settings and assessment of perceived risk from primaquine at the community level in Cambodia. PLOS ONE. 2020 Jan 31;15(1):e0228207. pmid:32004348

17. Oo NN, Bancone G, Maw LZ, Chowwiwat N, Bansil P, Domingo GJ, et al. Validation of G6PD Point-of-Care Tests among Healthy Volunteers in Yangon, Myanmar. PLOS ONE. 2016 Apr 1;11(4):e0152304. pmid:27035821

18. Espino FE, Bibit J-A, Sornillo JB, Tan A, Seidlein L von, Ley B. Comparison of Three Screening Test Kits for G6PD Enzyme Deficiency: Implications for Its Use in the Radical Cure of Vivax Malaria in Remote and Resource-Poor Areas in the Philippines. PLOS ONE. 2016 Feb 5;11(2):e0148172. pmid:26849445

19. Ley B, Satyagraha AW, Rahmat H, Fricken ME von, 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 Dec 13;16(12):e1002992. pmid:31834890

20. Shah SS, Diakite SAS, Traore K, Diakite M, Kwiatkowski DP, Rockett KA, et al. A novel cytofluorometric assay for the detection and quantification of glucose-6-phosphate dehydrogenase deficiency. Sci Rep. 2012 Mar 5;2:299. pmid:22393475

21. 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 [Internet]. 2019 Dec 13 [cited 2020 Dec 15];16(12). Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6910667/ pmid:31834890

22. Bancone G, Gornsawun G, Chu CS, Porn P, Pal S, Bansil P, et al. Validation of the quantitative point-of-care CareStart biosensor for assessment of G6PD activity in venous blood. Schallig HDFH, editor. PLOS ONE. 2018 May 8;13(5):e0196716. pmid:29738562

23. 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 Feb 7;14(2):e1002224. pmid:28170391

24. Pamba A, Richardson ND, Carter N, Duparc S, Premji Z, Tiono AB, et al. Clinical spectrum and severity of hemolytic anemia in glucose 6-phosphate dehydrogenase-deficient children receiving dapsone. Blood. 2012 Nov 15;120(20):4123–33. pmid:22993389

25. Satyagraha AW, Sadhewa A, Baramuli V, Elvira R, Ridenour C, Elyazar I, et al. G6PD Deficiency at Sumba in Eastern Indonesia Is Prevalent, Diverse and Severe: Implications for Primaquine Therapy against Relapsing Vivax Malaria. PLoS Negl Trop Dis. 2015 Mar 6;9(3):e0003602. pmid:25746733

26. Satyagraha AW, Sadhewa A, Elvira R, Elyazar I, Feriandika D, Antonjaya U, et al. Assessment of Point-of-Care Diagnostics for G6PD Deficiency in Malaria Endemic Rural Eastern Indonesia. PLoS Negl Trop Dis. 2016 Feb 19;10(2):e0004457. pmid:26894297

27. Brummaier T, Gilder ME, Gornsawun G, Chu CS, Bancone G, Pimanpanarak M, et al. Vivax malaria in pregnancy and lactation: a long way to health equity. Malar J. 2020 Jan 22;19(1):40. pmid:31969155

28. Ciulla TA, Sklar RM, Hauser SL. A simple method for DNA purification from peripheral blood. Anal Biochem. 1988 Nov 1;174(2):485–8. pmid:3239751

29. Saunders MA, Hammer MF, Nachman MW. Nucleotide Variability at G6pd and the Signature of Malarial Selection in Humans. Genetics. 2002 Dec 1;162(4):1849–61. pmid:12524354

30. Baird JK, Dewi M, Subekti D, Elyazar I, Satyagraha AW. Noninferiority of glucose-6-phosphate dehydrogenase deficiency diagnosis by a point-of-care rapid test vs the laboratory fluorescent spot test demonstrated by copper inhibition in normal human red blood cells. Transl Res. 2015 Jun;165(6):677–88. pmid:25312015

31. LaRue N, Kahn M, Murray M, Leader BT, Bansil P, McGray S, et al. Comparison of Quantitative and Qualitative Tests for Glucose-6-Phosphate Dehydrogenase Deficiency. Am J Trop Med Hyg. 2014 Oct 1;91(4):854–61. pmid:25071003

32. Chu CS, Bancone G, Kelley M, Advani N, Domingo GJ, Cutiongo-de la Paz EM, et al. Optimizing G6PD testing for Plasmodium vivax case management: why sex, counseling, and community engagement matter. Wellcome Open Res. 2020 Feb 7;5:21. pmid:32766454

33. Haston JC, Hwang J, Tan KR. Guidance for Using Tafenoquine for Prevention and Antirelapse Therapy for Malaria—United States, 2019. Morb Mortal Wkly Rep. 2019 Nov 22;68(46):1062–8. pmid:31751320

34. Shanks GD, Oloo AJ, Aleman GM, Ohrt C, Klotz FW, Braitman D, et al. A New Primaquine Analogue, Tafenoquine (WR 238605), for Prophylaxis against Plasmodium falciparum Malaria. Clin Infect Dis. 2001 Dec 15;33(12):1968–74. pmid:11700577

35. Rueangweerayut R, Bancone G, Harrell EJ, Beelen AP, Kongpatanakul S, Möhrle JJ, et al. Hemolytic Potential of Tafenoquine in Female Volunteers Heterozygous for Glucose-6-Phosphate Dehydrogenase (G6PD) Deficiency (G6PD Mahidol Variant) versus G6PD-Normal Volunteers. Am J Trop Med Hyg. 2017 Jul 24;tpmd160779. pmid:28749773

36. Chu CS, Bancone G, Soe NL, Carrara VI, Gornsawun G, Nosten F. The impact of using primaquine without prior G6PD testing: a case series describing the obstacles to the medical management of haemolysis. Wellcome Open Res. 2019 Apr 9;4:25. pmid:31069260

37. Brito-Sousa JD, Santos TC, Avalos S, Fontecha G, Melo GC, Val F, et al. Clinical Spectrum of Primaquine-induced Hemolysis in Glucose-6-Phosphate Dehydrogenase Deficiency: A 9-Year Hospitalization-based Study From the Brazilian Amazon. Clin Infect Dis Off Publ Infect Dis Soc Am. 2019 Sep 27;69(8):1440–2. pmid:30753364

38. 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

39. Baird JK. 8-Aminoquinoline Therapy for Latent Malaria. Clin Microbiol Rev. 2019 Sep 18;32(4). pmid:31366609

40. Bennett JW, Pybus BS, Yadava A, Tosh D, Sousa JC, McCarthy WF, et al. Primaquine Failure and Cytochrome P-450 2D6 in Plasmodium vivax Malaria. N Engl J Med. 2013 Oct 3;369(14):1381–2. pmid:24088113

41. Marcsisin SR, Reichard G, Pybus BS. Primaquine pharmacology in the context of CYP 2D6 pharmacogenomics: Current state of the art. Pharmacol Ther. 2016 May 1;161:1–10. pmid:27016470

42. Baird JK, Louisa M, Noviyanti R, Ekawati L, Elyazar I, Subekti D, et al. Association of Impaired Cytochrome P450 2D6 Activity Genotype and Phenotype With Therapeutic Efficacy of Primaquine Treatment for Latent Plasmodium vivax Malaria. JAMA Netw Open. 2018 Aug 31;1(4):e181449–e181449. pmid:30646129

43. 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 [Internet]. 2020 Feb 20 [cited 2020 Dec 15];18. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7031957/ pmid:32075649

44. Setyadi A, Arguni E, Kenangalem E, Hasanuddin A, Lampah DA, Thriemer K, et al. Safety of primaquine in infants with Plasmodium vivax malaria in Papua, Indonesia. Malar J. 2019 Apr 2;18(1):111. pmid:30940140

45. 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 Aug 29;14(8):e1002379. pmid:28850568

46. Pal S, Bansil P, Bancone G, Hrutkay S, Kahn M, Gornsawun G, et al. Evaluation of a Novel Quantitative Test for Glucose-6-Phosphate Dehydrogenase Deficiency: Bringing Quantitative Testing for Glucose-6-Phosphate Dehydrogenase Deficiency Closer to the Patient. Am J Trop Med Hyg. 2019 Jan;100(1):213–21. pmid:30350771

47. 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 Nov 2;13(11):e0206331. pmid:30388146

48. Recht J, Ashley EA, White NJ. Use of primaquine and glucose-6-phosphate dehydrogenase deficiency testing: Divergent policies and practices in malaria endemic countries. PLoS Negl Trop Dis. 2018 Apr 19;12(4):e0006230. pmid:29672516