- Open Access
Glucose-6-phosphate dehydrogenase mutations in malaria endemic area of Thailand by multiplexed high‐resolution melting curve analysis
Malaria Journal volume 20, Article number: 194 (2021)
Glucose-6-phosphate dehydrogenase (G6PD) deficiency, the most common enzymopathy in humans, is prevalent in tropical and subtropical areas where malaria is endemic. Anti-malarial drugs, such as primaquine and tafenoquine, can cause haemolysis in G6PD-deficient individuals. Hence, G6PD testing is recommended before radical treatment against vivax malaria. Phenotypic assays have been widely used for screening G6PD deficiency, but in heterozygous females, the random lyonization causes difficulty in interpreting the results. Over 200 G6PD variants have been identified, which form genotypes associated with differences in the degree of G6PD deficiency and vulnerability to haemolysis. This study aimed to assess the frequency of G6PD mutations using a newly developed molecular genotyping test.
A multiplexed high-resolution melting (HRM) assay was developed to detect eight G6PD mutations, in which four mutations can be tested simultaneously. Validation of the method was performed using 70 G6PD-deficient samples. The test was then applied to screen 725 blood samples from people living along the Thai–Myanmar border. The enzyme activity of these samples was also determined using water-soluble tetrazolium salts (WST-8) assay. Then, the correlation between genotype and enzyme activity was analysed.
The sensitivity of the multiplexed HRM assay for detecting G6PD mutations was 100 % [95 % confidence interval (CI): 94.87–100 %] with specificity of 100 % (95 % CI: 87.66–100 %). The overall prevalence of G6PD deficiency in the studied population as revealed by phenotypic WST-8 assay was 20.55 % (149/725). In contrast, by the multiplexed HRM assay, 27.17 % (197/725) of subjects were shown to have G6PD mutations. The mutations detected in this study included four single variants, G6PD Mahidol (187/197), G6PD Canton (4/197), G6PD Viangchan (3/197) and G6PD Chinese-5 (1/197), and two double mutations, G6PD Mahidol + Canton (1/197) and G6PD Chinese-4 + Viangchan (1/197). A broad range of G6PD enzyme activities were observed in individuals carrying G6PD Mahidol, especially in females.
The multiplexed HRM-based assay is sensitive and reliable for detecting G6PD mutations. This genotyping assay can facilitate the detection of heterozygotes, which could be useful as a supplementary approach for high-throughput screening of G6PD deficiency in malaria endemic areas before the administration of primaquine and tafenoquine.
Glucose-6-phosphate dehydrogenase (G6PD) deficiency is an inherited genetic defect and the most common enzymopathy, affecting approximately 500 million people worldwide with more than 200 variants have been identified . G6PD deficiency is prevalent in tropical and subtropical areas where malaria is endemic, including Africa and Southeast Asia . Evidence has suggested that G6PD deficiency confers protection against malaria infection [3,4,5]. However, this is still controversial because several studies have yielded contradictory results with some claiming that the protective effects of G6PD deficiency were observed in male hemizygotes only, in female heterozygotes only, or in both [6,7,8,9]. The major clinical concern associated with G6PD deficiency is haemolysis upon exposure to oxidant drugs, including anti-malarials such as 8-aminoquinolines (primaquine and tafenoquine) [10,11,12,13]. Primaquine and tafenoquine are the only medications capable of killing Plasmodium vivax and Plasmodium ovale at the dormant liver stage (hypnozoite). The World Health Organization (WHO) recommends that G6PD activity be measured before efforts to perform radical treatment of malaria .
G6PD deficiency can be diagnosed by either phenotypic or genotypic assay. Phenotypic tests are based on the assessment of G6PD activity, measuring the production of reduced nicotinamide adenine dinucleotide phosphate (NADPH), which can be done quantitatively. The standard quantitative method is spectrophotometry, in which NADPH production is monitored at 340 nm . This method is accurate and reliable, but is laborious, time-consuming, and requires complicated sample preparation and technical skills; as such, it is not commonly used for field-based screening. A colorimetric G6PD assay, based on water-soluble tetrazolium salts (WST-8), was developed as an alternative to the gold standard of spectrophotometry . In this approach, no sample preparation is required and whole blood or dried blood spots can be used to test G6PD activity . The WST-8-based assay can be used as a quantitative method or a qualitative one by the naked eye, offering the possibility of performing mass screening of G6PD deficiency in the context of malaria elimination using primaquine and tafenoquine [16, 18]. Although not the standard method for measuring G6PD activity, the sensitivity of WST-8 for detecting NAD(P)H was found to be five-fold greater than that of the spectrophotometric assay. Moreover, results obtained by measuring dehydrogenase activities in biological samples using WST-8 assay were in parallel with the standard method . For G6PD testing, WST-8 was applied, in 96-well format, to the screening of G6PD deficiency in different populations [20,21,22]. The sensitivity and specificity of WST-8 for detecting G6PD activity < 30 % were 55 % and 98 %, respectively, compared with the spectrophotometric method . In addition, sensitivity of 72 % and specificity of 98 % were reported for WST-8, in comparison with the standard quantitative G6PD assay (R&D Diagnostics) . This suggests that WST-8 could be a key tool for G6PD testing, but it requires further development before deployment in the field.
G6PD diagnostic tests are currently available, including qualitative tests such as fluorescent spot test (FST) and CareStart™ G6PD rapid diagnostic test, as well as quantitative point-of-care tests such as CareStart™ G6PD biosensor and STANDARD™ G6PD test. Unfortunately, these tests are not widely used for G6PD testing because they are too expensive and can be difficult to interpret [23,24,25]. Qualitative tests are reliable for identifying G6PD deficiency in hemizygous males and homozygous females, but are unable to identify heterozygous females [26,27,28]. This is because, in heterozygous females, a wide range of G6PD activities are observed as a result of the random X-chromosome inactivation or lyonization . To date, over 200 G6PD variants have been identified, which form genotypes associated with differences in the degree of deficiency and vulnerability to haemolysis . Moreover, G6PD activities vary among G6PD-deficient individuals carrying the same genotype [31, 32].
G6PD genotyping can be performed using restriction fragment length polymorphism [33, 34], amplification refractory mutation system [35, 36], gold nanoparticles-based assay , high resolution melting (HRM) curve analysis [38, 39] and DNA sequencing [40, 41]. Additionally, multiplex genotyping systems are currently available. DiaPlexC™ G6PD Genotyping Kit (Asian type) can detect eight mutations, namely, G6PD Vanua Lava, G6PD Mahidol, G6PD Mediterranean, G6PD Coimbra, G6PD Viangchan, G6PD Union, G6PD Canton, and G6PD Kaiping. Thus, this assay offers high-throughput screening of G6PD mutations by one-step PCR . However, after PCR amplification, an additional gel electrophoresis step is required to check the size of the amplified PCR products, which is impractical for large population screening. The HRM assay is a powerful and reliable tool that has been widely used in the detection of gene mutations [43,44,45]. Previously, HRM assays were applied to detect G6PD mutations in different population groups [38, 46,47,48]. However, previous HRM assays could detect only one or two mutations at a time. Although a multiplexed system to detect six mutations in four reactions was later described, the assay system and interpretation of results were complex .
The prevalence of G6PD deficiency in Thailand ranges between 5 and 18 %, depending on the geographical area [50,51,52,53,54]. More than 20 G6PD variants have been identified in the country, among which the most common is G6PD Viangchan, followed by G6PD Mahidol, G6PD Canton, G6PD Union, G6PD Kaiping, G6PD Gaohe, G6PD Chinese-4, G6PD Chinese-5, G6PD Valladolid, G6PD Coimbra and G6PD Aures. Along the Thai–Myanmar border, a malaria endemic area, prevalence of G6PD deficiency of 9–18 % was reported in males . Moreover, a rate of G6PD deficiency of 7.4 % was reported from the screening of 1,340 newborns . G6PD Mahidol was shown to be the most common variant in this population, accounting for 88 % of all variants, followed by G6PD Chinese-4, G6PD Viangchan, and G6PD Mediterranean. Generally, to avoid the risk of haemolysis upon malaria treatment, G6PD testing is recommended before the administration of primaquine and tafenoquine. The aim of this study was to develop a molecular diagnostic test to enable an accurate, reliable and high-throughput platform for detecting G6PD mutations, which can be used as a supplement to the screening of G6PD deficiency, especially in heterozygous females. To validate the method, 70 G6PD-deficient and 28 non-deficient samples were tested and the results were compared with the findings obtained by direct DNA sequencing. The potential utility of the developed HRM test for the detection of G6PD variants in a study area in Thailand was then examined. The correlation between genotype and the phenotype of enzyme activity (as determined using WST-8) was also analysed.
Blood samples were collected in ethylenediaminetetraacetic acid (EDTA) tubes and transported to the laboratory under storage at 4 °C. Thereafter, samples were stored at −20 °C until use, for approximately 1−3 months. Under these conditions, the integrity of samples for phenotypic analysis was maintained as it was recently reported that blood samples were stable for up to 7–12 months when stored in EDTA tubes at − 20 °C .
For the validation of HRM assays, 70 G6PD-deficient and 28 non-deficient blood samples were collected from healthy Thai volunteers at the Faculty of Medicine Ramathibodi Hospital. All samples were spectrophotometrically tested for G6PD activity and genotyped by DNA sequencing. Ethical approval for the study was provided by the Committee on Human Rights Related to Research Involving Human Subjects, Faculty of Medicine Ramathibodi Hospital, Mahidol University, Bangkok, Thailand (approval number MURA 2018/252).
For the screening of G6PD deficiency, 725 blood samples (from 368 males and 357 females) were collected in EDTA tubes from residents living along the Thai–Myanmar border, a malaria endemic area, namely, in Tha Song Yang District, Tak Province, Thailand. Ethical approval for the study was provided by the Human Ethics Committee, Faculty of Tropical Medicine, Mahidol University (approval number MUTM 2019-016-01).
Phenotypic screening of G6PD deficiency using WST-8 assay
WST-8 is not a standard method for measuring G6PD activity. However, this assay was used for phenotypic screening in this study because its performance was found to be indistinguishable from that of the spectrophotometric method involving measurement of the absorbance of NAD(P)H at 340 nm . The method showed high accuracy with % relative error of 0.7–0.25. For precision, % coefficient of variation for within-run and between-run of the WST-8 method ranged between 0.6 and 4.5. WST-8 also exhibited excellent reproducibility with Z′ values of 0.90–0.99. Although WST-8 provides advantages regarding the diagnosis of G6PD deficiency, this method will require further development before being deployed in a clinical context .
Reaction mixtures of 100 µl, consisting of 20 mM Tris-HCl pH 8.0, 10 mM MgCl2, 500 µM glucose-6-phosphate (G6P), 100 µM NADP+, and 100 µM WST-8 (Sigma-Aldrich, Darmstadt, Germany), were mixed with a blood sample of 2 µl in a 96-well plate. The absorbance was measured at 450 nm with a reference at 650 nm using a microplate reader (Sunrise; Tecan, Männedorf, Switzerland). The absorbance at 450 nm of a reaction mixture set up in the absence of G6P substrate was used for background subtraction. The G6PD activity was calculated using an NADPH calibration curve. Haemoglobin concentration was measured using Drabkin’s reagent (Sigma-Aldrich). G6PD activity was reported as units (U) per gram of haemoglobin (gHb). Experiments were performed in triplicate.
DNA extraction was performed using the QIAsymphony DNA Mini Kit (QIAGEN, Hilden, Germany), in accordance with the manufacturer’s instructions. Blood samples of 100 µl were extracted and eluted into a final volume of 50 µl. DNA concentration was measured using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA).
Primers were designed to detect eight common G6PD mutations in the Thai population: G6PD Gaohe (A95G), G6PD Chinese-4 (G392T), G6PD Mahidol (G487A), G6PD Viangchan (G871A), G6PD Chinese-5 (C1024T), G6PD Union (C1360T), G6PD Canton (G1376T) and G6PD Kaiping (G1388A; Table 1). The primers were designed to detect the mutations by generating PCR products with distinctive melting temperatures (Tm; Fig. 1).
PCR amplification and melting curve analysis
Assay conditions, including primer concentrations, assay protocol, and detection conditions, were optimized to maximize the sensitivity and specificity of the assay and to minimize the cross-reactivity. Multiplexed HRM assay was performed in a total volume of 12.5 µl, containing 6.25 µl of 2× HRM Type-It mix (QIAGEN), various concentrations of each primer (Table 1), molecular-grade water and 2.5 µl of the gDNA template (3–10 ng/µl). PCR amplification and melting curve analysis were performed using the Rotor-Gene Q (QIAGEN) with the following conditions: 1 cycle of 95 °C for 5 min, and then 30 cycles of 95 °C for 10 s, 63 °C for 30 s, and 72 °C for 10 s. Subsequently, HRM analysis was performed by melting from 75 to 90 °C, reading at every 0.1 °C step with 2 s of stabilization. Positive (gDNA with known mutations, confirmed by DNA sequencing) and negative controls (gDNA of G6PD wild-type (WT), confirmed by DNA sequencing) were included in every run. Data analysis was carried out using the Rotor-Gene Q software. Experiments were performed in triplicate.
PCR amplification and DNA sequencing
To validate the HRM results, PCR and sequencing primers were designed, as shown in Table 2. For DNA amplification, extracted gDNA was used as a template. The g6pd gene was amplified using four primer sets (Exon2F−Exon2R, Exon3F−Exon5R, Exon6F−Exon8R, and Exon9F−Exon13R), which cover all 13 exons. The PCR reaction was set up in a final volume of 50 µl, containing 1⋅ Taq Buffer with (NH4)2SO4, 2.5 mM MgCl2, 200 µM of each dNTP, 0.25 µM of each primer, 50 ng gDNA and 1.25 U of Taq DNA polymerase (Thermo Fisher Scientific). The thermal cycling profile was as follows: initial denaturation at 95 °C for 3 min; 35 cycles of denaturation at 95 °C for 30 s, annealing for 30 s, and extension at 72 °C for 1 min; followed by final extension at 72 °C for 10 min. The annealing temperature was 60 °C for the primers Exon2F−Exon2R, Exon3F−Exon5R, and Exon6F−Exon8R, and 65 °C for Exon9F−Exon13R. PCR products were subjected to gel purification and sequenced (Bio Basic, Ontario, Canada).
Data are presented as mean ± SD. Statistical analyses and plotting of graphs were performed using GraphPad Prism (GraphPad Software, La Jolla, CA, USA). To assess the performance of multiplexed HRM in the detection of G6PD mutations, the numbers of true positives, true negatives, false positives, and false negatives were determined. The following parameters were calculated: sensitivity = [true positives/(true positives + false negatives)] ⋅100; specificity = [true negatives/(true negatives + false positives)] ⋅100; positive predictive value = [true positives/(true positives + false positives)] ⋅100; and negative predictive value = [true negatives/(true positives + false negatives)] ⋅100.
Development and validation of 4-plex HRM assay
A multiplexed HRM assay was developed to detect eight G6PD variants that are common in Thailand in two reactions (Fig. 1). By using a specific primer pair for each mutation, reaction 1 simultaneously detects four mutations [G6PD Gaohe (A95G), G6PD Mahidol (G487A), G6PD Viangchan (G871A), and G6PD Canton (G1376T)]. Reaction 2 concurrently detects another four mutations [G6PD Chinese-4 (G392T), G6PD Chinese-5 (C1024T), G6PD Union (C1360T), and G6PD Kaiping (G1388A)]. The assay is based on a single fluorescent dye, EvaGreen, without the need for a quenching probe. The primers were designed to detect the mutations by generating PCR products with distinctive melting temperatures, Tm. In contrast, no amplification occurred in WT samples. A peak at the corresponding Tm reveals the genotype of each sample. The gDNA of known G6PD mutations was used as positive controls. Overall, 70 G6PD-deficient samples and 28 non-deficient samples were used to evaluate the performance of the developed 4-plex HRM assay, while direct DNA sequencing was used as a reference test (Table 3).
In comparison to direct DNA sequencing, the 4-plex HRM assay was 100 % sensitive [95 % confidence interval (CI): 94.87–100 %] and 100 % specific (95 % CI: 87.66–100 %), with no cross-reactivity for the detection of G6PD mutations (Table 4). Additionally, the multiplexed HRM assay could correctly identify the double mutations (G6PD Mahidol + Canton and G6PD Gaohe + Kaiping). This indicates that the developed method is reliable for detecting G6PD mutations.
Phenotypic screening of G6PD deficiency by WST-8 assay
The prevalence of G6PD deficiency in people living in a malaria endemic area in Thailand, namely, Tha Song Yang District, Tak Province, was determined by G6PD activity assay (WST-8). Figure 2 indicates the G6PD enzyme activity of 725 samples measured by WST-8. The average G6PD activity in males and females was 9.99 ± 4.14 and 10.35 ± 3.81 U/gHb, respectively. The adjusted male median (AMM) value was determined (10.31 ± 3.81 U/gHb) and defined as 100 % G6PD activity . The WHO defined G6PD activity of less than 30 % as deficient and G6PD activity ranging between 30 and 80 % as intermediate .
Nonetheless, G6PD activity of 70 % was used as a threshold for tafenoquine prescription [58, 59]. In this study, G6PD activity levels of less than 30 % and 30–70 % of the AMM were thus considered as deficient and intermediate, respectively. Subjects with G6PD activity over 70 % of the AMM were defined as normal. Based on the WST-8 assay, the prevalence of G6PD deficiency in the studied population was 20.55 % (149/725; Table 5). Prevalence rates of G6PD deficiency of 20.11 % (74/368) and 21.01 % (75/357) were observed in males and females, respectively. In addition, average G6PD activity of deficient males and females was 1.59 ± 0.89 and 1.69 ± 0.77 U/gHb, respectively. Intermediate G6PD activity (30–70 %) was found in 7.34 % (27/368) of males and 16.25 % (58/357) of females. Average G6PD activity of non-deficient (> 70 %) cases was 11.78 ± 2.11 U/gHb in males and 11.89 ± 2.49 U/gHb in females. The frequency distribution of G6PD activity of the 725 individuals measured by WST-8 is shown in Fig. 3a. The majority of the enzyme activities were distributed between 7 and 16 U/gHb. The frequency distribution of G6PD activity by sex is illustrated in Fig. 3b, c. A broader distribution of G6PD activities was seen in females than in males.
Genotypic screening of G6PD deficiency using the multiplexed HRM assay
The developed 4-plex HRM assay was applied to screen for G6PD mutations in 725 blood samples. This assay identified 197 of the 725 (27.17 %) individuals as possessing at least one mutation with an adverse effect on function (Table 6). The prevalence of subjects carrying at least one G6PD mutation was 20.11 % (74/368) in males and 34.45 % (123/357) in females. The most common G6PD mutation detected in the studied population was G6PD Mahidol, accounting for 94.92 % of the total (n = 187; 72 in males and 115 in females). Other single mutations observed in the study included G6PD Canton (2.03 %; 4 in females), G6PD Viangchan (1.52 %; 1 in a male and 2 in females), and G6PD Chinese-5 (0.51 %; 1 in a male). The HRM assay could also detect the double mutant variants, which were G6PD Mahidol + Canton (0.51 %; 1 in a female) and G6PD Chinese-4 + Viangchan (0.51 %; 1 in a female). Figure 4 shows the G6PD activity of deficient and normal samples identified by HRM for males and females. G6PD enzyme activity of deficient subjects, especially in females, spanned from the deficient region (< 30 %) to the normal region (> 70 %). A large distribution of G6PD enzyme activities in females is caused by genetic mosaicism as a result of X-inactivation. The distribution of G6PD activity by mutation type is illustrated in Fig. 5. Non-variant individuals are also included in this plot. Variation of G6PD activities among the different mutations was observed. Moreover, compared with the G6PD enzyme activity in males with the same mutation, that in females was greater. Enzyme activity of 0.89 and 6.16 U/gHb was observed for G6PD Viangchan in males and females, respectively. Interestingly, G6PD Mahidol, a Class III variant with mild deficiency, which was the most prevalent variant in the studied population, exhibited a wide range of G6PD activities, in both males (range: 0.10–10.73 U/gHb, mean: 3.20 ± 2.46 U/gHb) and females (range: 0.10–17.72 U/gHb, mean: 7.72 ± 4.24 U/gHb). Notably, G6PD enzyme activity in the double mutant variants (G6PD Mahidol + Canton and G6PD Chinese-4 + Viangchan) was significantly decreased compared with that of the single mutants.
HRM assays have been widely used to detect gene mutations [43,44,45]. However, for G6PD genotyping, most of the developed HRM assays are singleplex or duplex, which can only detect one or two mutations simultaneously [38, 46, 47]. To enable the detection of multiple mutations, more than one fluorescent dye must be included in the reaction mixture, which usually makes the reaction more expensive [39, 60, 61].
As reported in this paper, a 4-plex HRM assay for the detection of G6PD mutations common in Thailand and elsewhere in Asia was developed, using a single dye with a run time of 80 min. Evaluation of the 4-plex HRM assay using 70 G6PD-deficient samples indicated that it was accurate and reliable for detecting G6PD mutations (with specificity and sensitivity of 100 %) compared with DNA sequencing. Among the 70 deficient assay validation samples, G6PD Viangchan was the most prevalent variant, followed by G6PD Canton and G6PD Kaiping. This is in accordance with previous reports showing that G6PD Viangchan was the most common variant in Thais [28, 51]. Two double mutations, G6PD Mahidol + Canton and G6PD Gaohe + Kaiping, were also identified in the studied population. G6PD Mahidol + Canton was first identified in people living along the Thai–Myanmar border where Karen and Burman are the major population groups .
After validation, the multiplexed HRM assay was applied to screen G6PD mutations of 725 people living in a malaria endemic area along the Thai–Myanmar border. The prevalence of G6PD deficiency in this population was also determined by phenotypic enzyme activity assay using WST-8. Considering a 30 % activity cut-off, the overall prevalence of G6PD deficiency was 8.83 % by WST-8 assay. If the upper limit of a 70 % cut-off was considered, the overall prevalence increased to 20.55 %. By sex, at the 70 % cut-off, the prevalence of G6PD deficiency was 20.11 % in males and 21.01 % in females. In contrast, by the multiplexed HRM assay, the frequency of G6PD mutations was 20.11 and 34.45 % in males and females, respectively. Thus, the prevalence of G6PD deficiency in males is equivalent between these two assays, but in females, genetic analysis using HRM indicated a high frequency of g6pd gene mutations (34.45 %), that is notably greater than the prevalence of G6PD deficiency measured by WST-8 (21.01 %). The multiplexed HRM assay identified four single mutants (95.11 % G6PD Mahidol, 2.03 % G6PD Canton, 1.52 % G6PD Viangchan, and 0.51 % G6PD Chinese-5) and two double mutants (0.51 % G6PD Mahidol + Canton and 0.51 % G6PD Chinese-4 + Viangchan) in the studied population. In good agreement with previous reports, G6PD Mahidol was the most common variant in the Karen population [26, 62]. The double mutant G6PD Chinese-4 + Viangchan was identified here for the first time.
A broad range of G6PD activities were observed among the different genotypes. G6PD Mahidol showed ranges of enzyme activity of 0.10–10.73 and 0.10–17.72 U/gHb in males and females, respectively. For G6PD Canton, the range of enzyme activity was 6.50–10.48 U/gHb in females. A wider distribution of enzyme activities was observed in females carrying G6PD Mahidol than in males. Additionally, 5 males and 58 females carrying G6PD Mahidol and 3 females carrying G6PD Canton showed G6PD activity even greater than the 70 % activity cut-off (7 U/gHb). Similar findings were previously described in genotype–phenotype association studies [31, 32]. This is mainly attributable to the fact that the degree of deficiency and vulnerability to haemolysis vary substantially among the different genotypes. Furthermore, because G6PD deficiency is an X-linked genetic defect, males exhibit hemizygous deficiency while females can exhibit either homozygous or heterozygous deficiency. In heterozygous females, a wide range of G6PD activities can be observed, in which the activity of the normal erythrocyte population may compensate for the lost activity of the deficient erythrocyte population. In addition, abnormally high G6PD activity observed in persons carrying G6PD mutations could also be attributable to the following factors. The first is the limited performance quality of the WST-8 used in this study. The second is elevated reticulocyte count in blood samples. Young red blood cells usually exhibit greater G6PD enzyme activity than aged red blood cells. The last is the presence of leukocytes in tested samples (in this study, G6PD activity was measured from whole blood samples). This is because leukocytes retain their nucleus and therefore can continuously synthesize the G6PD enzyme.
G6PD-deficient individuals have increased susceptibility to haemolysis upon exposure to oxidative agents, including primaquine and tafenoquine which are the only medications effective for the radical treatment of infection by P. vivax and P. ovale. To ensure successful and safe treatment, a single dose of 300 mg of tafenoquine or 0.5 mg/kg/day primaquine for 14 days should be prescribed only in patients with at least 70 % G6PD activity . The effect of the primaquine dose on haemolysis was reported to differ between G6PD normal and G6PD Mahidol heterozygous individuals . These heterozygous females were identified as normal by FST assay while the quantitative assay revealed G6PD activity of 62 and 99 %. Tafenoquine was also reported to cause dose-dependent haemolysis in G6PD Mahidol heterozygous individuals with enzyme activity of 61–80 % . As such, drug-induced haemolysis associated with G6PD deficiency depends on two major factors: the first is the level of G6PD activity, which is determined by G6PD genotype, and the second is the exposure to oxidative stress, namely, metabolites of antimalarial drugs (8-aminoquinolines). Therefore, the drug dose and the ability to metabolize the parent compounds also contribute to the severity of haemolysis in malaria treatment.
Currently, the cut-off of 70 % AMM is widely accepted as an appropriate threshold to determine whether or not to administer tafenoquine [58, 59, 63]. However, the cut-off value should be carefully defined and tested in each population group. Based on the obtained results, to enable accurate determination of the prevalence of G6PD deficiency using WST-8, different cut-off values are required in males and females. The upper limit of 70 % of AMM is recommended for males. However, in heterozygous females, neither the lower (30 %) nor the upper (70 %) limit is reliable for screening G6PD deficiency. It should be noted that the results reported here are based on the WST-8 assay, which is an alternative to the standard method for measuring G6PD activity. Upon using other G6PD tests, the results might be different. For male populations, phenotype tests are useful for G6PD deficiency screening and to enable safe treatment. In contrast, in heterozygous females in whom a wide range of G6PD activities are observed, phenotypic enzyme assay alone might be insufficient to identify G6PD deficiency. Hence, alternative approaches such as genetic analysis could be useful for determining whether drugs should be administered in populations suspected of having G6PD deficiency. The multiplexed HRM assay developed here could be useful for identifying G6PD variants in Thai populations. Although other multiplex systems for genetic analysis are currently available, they might be unsuitable for large population screening. The G6PD gene chip kit can detect 13–14 mutations common in Chinese populations, but must be combined with a hybridization kit, which is time-consuming [32, 65]. The DiaPlexC™ (Asian type) can simultaneously detect eight mutations, but requires an additional gel electrophoresis step to check the amplified products . Additionally, the kit might not be applicable for deployment in regions where populations carry other G6PD mutations (e.g., G6PD Gaohe, G6PD Chinese-4, and G6PD Chinese-5), for which the kit cannot test.
It should be mentioned that, in this study, no mutation was detected in 12 females who were considered likely to be G6PD-deficient because the observed G6PD activity was lower than 7 U/gHb. This might have been because of the limited performance quality (sensitivity of 55–72 %) of the WST-8 assay used in this study [20, 21]. Alternatively, these subjects might carry G6PD mutations for which the multiplexed HRM assay cannot test. DNA sequencing of the whole g6pd gene might be required to detect mutations in such individuals. Additionally, to enable mutational screening in more diverse population groups, the assay should be expanded to include other mutations, such as G6PD Mediterranean, G6PD Valladolid, G6PD Coimbra, and G6PD Aures . It should also be noted that the multiplexed HRM assays developed here are not able to identify zygosity of the samples and, thus, require further development before being deployed. Nevertheless, the HRM assays could be of great use for analysing G6PD mutations in supplement to phenotypic G6PD screening in heterozygous females as well as in populations suspected of having G6PD deficiency. Primarily, G6PD genotyping is being done in G6PD-deficient individuals. However, more data on the genotype–phenotype association of G6PD deficiency in diverse population groups should be obtained, which requires a high-throughput screening platform.
A multiplexed HRM assay for the detection of eight common G6PD mutations in Thailand was developed. The performance of the assay was excellent, with 100 % specificity and 100 % sensitivity. The prevalence of G6PD mutations in 725 people living in a malaria endemic area along the Thai–Myanmar border was determined to be 27.17 % by HRM, which is greater than the prevalence of G6PD deficiency determined by the WST-8 phenotypic assay (20.55 %). Performing a phenotypic assay alone might thus be inadequate and the result might not be an accurate predictor of G6PD deficiency, especially in heterozygous females. As an option to overcome this problem, the multiplexed HRM assay is rapid, accurate and reliable for detecting G6PD mutations, enabling high-throughput screening. This assay could be useful as a supplementary approach for high-throughput screening of G6PD deficiency before the administration of 8-aminoquinolones in malaria endemic areas.
Availability of data and materials
All data analysed during the study are included in this published article.
Adjusted male median
Fluorescent spot test
Reduced nicotinamide adenine dinucleotide phosphate
World Health Organization
Water-soluble tetrazolium salts
Luzzatto L, Ally M, Notaro R. Glucose-6-phosphate dehydrogenase deficiency. Blood. 2020;136:1225–40.
Nkhoma ET, Poole C, Vannappagari V, Hall SA, Beutler E. The global prevalence of glucose-6-phosphate dehydrogenase deficiency: a systematic review and meta-analysis. Blood Cells Mol Dis. 2009;42:267–78.
Leslie T, Briceno M, Mayan I, Mohammed N, Klinkenberg E, Sibley CH, et al. The impact of phenotypic and genotypic G6PD deficiency on risk of Plasmodium vivax infection: a case-control study amongst Afghan refugees in Pakistan. PLoS Med. 2010;7:e1000283.
Louicharoen C, Patin E, Paul R, Nuchprayoon I, Witoonpanich B, Peerapittayamongkol C, et al. Positively selected G6PD-Mahidol mutation reduces Plasmodium vivax density in Southeast Asians. Science. 2009;326:1546–9.
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.
Mbanefo EC, Ahmed AM, Titouna A, Elmaraezy A, Trang NT, Phuoc Long A, et al. Association of glucose-6-phosphate dehydrogenase deficiency and malaria: a systematic review and meta-analysis. Sci Rep. 2017;7:45963.
Bienzle U, Ayeni O, Lucas AO, Luzzatto L. Glucose-6-phosphate dehydrogenase and malaria. Greater resistance of females heterozygous for enzyme deficiency and of males with non-deficient variant. Lancet. 1972;1:107–10.
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:246–9.
Uyoga S, Ndila CM, Macharia AW, Nyutu G, Shah S, Peshu N, et al. Glucose-6-phosphate dehydrogenase deficiency and the risk of malaria and other diseases in children in Kenya: a case-control and a cohort study. Lancet Haematol. 2015;2:e437-44.
Carmona-Fonseca J, Alvarez G, Maestre A. Methemoglobinemia and adverse events in Plasmodium vivax malaria patients associated with high doses of primaquine treatment. Am J Trop Med Hyg. 2009;80:188–93.
Beutler E. Glucose-6-phosphate dehydrogenase deficiency: a historical perspective. Blood. 2008;111:16–24.
Cappellini MD, Fiorelli G. Glucose-6-phosphate dehydrogenase deficiency. Lancet. 2008;371:64–74.
Rueangweerayut R, Bancone G, Harrell EJ, Beelen AP, Kongpatanakul S, Mohrle 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;97:702–11.
WHO. Testing for G6PD deficiency for safe use of primaquine in radical cure of P. vivax and P. ovale malaria. Geneva: World Health Organization; 2016.
WHO. Standardization of procedures for the study of glucose-6-phosphate dehydrogenase. Report of a WHO Scientific Group. Geneva, World Health Organization. 1967:1–53.
Tantular IS, Kawamoto F. An improved, simple screening method for detection of glucose-6-phosphate dehydrogenase deficiency. Trop Med Int Health. 2003;8:569–74.
Arai M, Kosuge K, Kawamoto F, Matsuoka H. Reactivity of blood samples spotted onto filter papers in the WST-8 method for screening of G6PD deficiency. Acta Med Okayama. 2006;60:127–34.
Jalloh A, Tantular IS, Pusarawati S, Kawilarang AP, Kerong H, Lin K, et al. Rapid epidemiologic assessment of glucose-6-phosphate dehydrogenase deficiency in malaria-endemic areas in Southeast Asia using a novel diagnostic kit. Trop Med Int Health. 2004;9:615–23.
Chamchoy K, Pakotiprapha D, Pumirat P, Leartsakulpanich U, Boonyuen U. Application of WST-8 based colorimetric NAD(P)H detection for quantitative dehydrogenase assays. BMC Biochem. 2019;20:4.
Ley B, Alam MS, O’Donnell JJ, Hossain MS, Kibria MG, Jahan N, et al. Comparison of three quantitative methods to estimate G6PD activity in the Chittagong Hill Tracts, Bangladesh. PLoS ONE. 2017;12:e0169930.
De Niz M, Eziefula AC, Othieno L, Mbabazi E, Nabukeera D, Ssemmondo E, et al. Tools for mass screening of G6PD deficiency: validation of the WST8/1-methoxy-PMS enzymatic assay in Uganda. Malar J. 2013;12:210.
Kuwahata M, Wijesinghe R, Ho MF, Pelecanos A, Bobogare A, Landry L, et al. Population screening for glucose-6-phosphate dehydrogenase deficiencies in Isabel Province, Solomon Islands, using a modified enzyme assay on filter paper dried bloodspots. Malar J. 2010;9:223.
Brito MA, Peixoto HM, Almeida AC, Oliveira MR, Romero GA, Moura-Neto JP, et al. Validation of the rapid test Carestart™ G6PD among malaria vivax-infected subjects in the Brazilian Amazon. Rev Soc Bras Med Trop. 2016;49:446–55.
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;10:e0004457.
Oo NN, Bancone G, Maw LZ, Chowwiwat N, Bansil P, Domingo DJ, et al. Validation of G6PD point-of-care tests among healthy volunteers in Yangon, Myanmar. PLoS ONE. 2016;11:e0152304.
Bancone G, Chu CS, Somsakchaicharoen R, Chowwiwat N, Parker DM, Charunwatthana P, et al. Characterization of G6PD genotypes and phenotypes on the northwestern Thailand-Myanmar border. PLoS ONE. 2014;9:e116063.
Thielemans L, Gornsawun G, Hanboonkunupakarn B, Paw MK, Porn P, Moo PK, et al. Diagnostic performances of the fluorescent spot test for G6PD deficiency in newborns along the Thailand-Myanmar border: a cohort study. Wellcome Open Res. 2018;3:1.
Nantakomol D, Paul R, Palasuwan A, Day NP, White NJ, Imwong M. Evaluation of the phenotypic test and genetic analysis in the detection of glucose-6-phosphate dehydrogenase deficiency. Malar J. 2013;12:289.
Harper PS. Mary Lyon and the hypothesis of random X chromosome inactivation. Hum Genet. 2011;130:169–74.
Gomez-Manzo S, Marcial-Quino J, Vanoye-Carlo A, Serrano-Posada H, Ortega-Cuellar D, Gonzalez-Valdez A, et al. Glucose-6-phosphate dehydrogenase: update and analysis of new mutations around the world. Int J Mol Sci. 2016;17:2069.
Powers JL, Best DH, Grenache DG. Genotype–phenotype correlations of glucose-6-phosphate–deficient variants throughout an activity distribution. J Appl Lab Med. 2018;2:841–50.
He Y, Zhang Y, Chen X, Wang Q, Ling L, Xu Y. Glucose-6-phosphate dehydrogenase deficiency in the Han Chinese population: molecular characterization and genotype-phenotype association throughout an activity distribution. Sci Rep. 2020;10:17106.
D’Urso M, Luzzatto L, Perroni L, Ciccodicola A, Gentile G, Peluso I, et al. An extensive search for RFLP in the human glucose-6-phosphate dehydrogenase locus has revealed a silent mutation in the coding sequence. Am J Hum Genet. 1988;42:735–41.
Yoshida A, Takizawa T, Prchal JT. RFLP of the X chromosome-linked glucose-6-phosphate dehydrogenase locus in blacks. Am J Hum Genet. 1988;42:872–6.
Maffi D, Pasquino MT, Caprari P, Caforio MP, Cianciulli P, Sorrentino F, et al. Identification of G6PD Mediterranean mutation by amplification refractory mutation system. Clin Chim Acta. 2002;321:43–7.
Du CS, Ren X, Chen L, Jiang W, He Y, Yang M. Detection of the most common G6PD gene mutations in Chinese using amplification refractory mutation system. Hum Hered. 1999;49:133–8.
Seow N, Lai PS, Yung LY. Gold nanostructures for the multiplex detection of glucose-6-phosphate dehydrogenase gene mutations. Anal Biochem. 2014;451:56–62.
Yan JB, Xu HP, Xiong C, Ren ZR, Tian GL, Zeng F, et al. Rapid and reliable detection of glucose-6-phosphate dehydrogenase (G6PD) gene mutations in Han Chinese using high-resolution melting analysis. J Mol Diagn. 2010;12:305–11.
Liu Z, Yu C, Li Q, Cai R, Qu Y, Wang W, et al. Chinese newborn screening for the incidence of G6PD deficiency and variant of G6PD gene from 2013 to 2017. Hum Mutat. 2020;41:212–21.
Vulliamy TJ, D’Urso M, Battistuzzi G, Estrada M, Foulkes NS, Martini G, et al. Diverse point mutations in the human glucose-6-phosphate dehydrogenase gene cause enzyme deficiency and mild or severe hemolytic anemia. Proc Natl Acad Sci USA. 1988;85:5171–5.
Poggi V, Town M, Foulkes NS, Luzzatto L. Identification of a single base change in a new human mutant glucose-6-phosphate dehydrogenase gene by polymerase-chain-reaction amplification of the entire coding region from genomic DNA. Biochem J. 1990;271:157–60.
Lee J, Kim TI, Kang JM, Jun H, Le HG, Thai TL, et al. Prevalence of glucose-6-phosphate dehydrogenase (G6PD) deficiency among malaria patients in Upper Myanmar. BMC Infect Dis. 2018;18:131.
Twist GP, Gaedigk R, Leeder JS, Gaedigk A. High-resolution melt analysis to detect sequence variations in highly homologous gene regions: application to CYP2B6. Pharmacogenomics. 2013;14:913–22.
Thomas V, Mazard B, Garcia C, Lacan P, Gagnieu MC, Joly P. UGT1A1 (TA)n genotyping in sickle-cell disease: high resolution melting (HRM) curve analysis or direct sequencing, what is the best way? Clin Chim Acta. 2013;424:258–60.
Edwards T, Sasaki S, Williams C, Hobbs G, Feasey NA, Evans K, et al. Speciation of common Gram-negative pathogens using a highly multiplexed high resolution melt curve assay. Sci Rep. 2018;8:1114.
Joly P, Lacan P, Garcia C, Martin C, Francina A. Rapid genotyping of two common G6PD variants, African (A-) and Mediterranean, by high-resolution melting analysis. Clin Biochem. 2010;43:193–7.
Pan M, Lin M, Yang L, Wu J, Zhan X, Zhao Y, et al. Glucose-6-phosphate dehydrogenase (G6PD) gene mutations detection by improved high-resolution DNA melting assay. Mol Biol Rep. 2013;40:3073–82.
Yang H, Wang Q, Zheng L, Zhan XF, Lin M, Lin F, et al. Incidence and molecular characterization of glucose-6-phosphate dehydrogenase deficiency among neonates for newborn screening in Chaozhou, China. Int J Lab Hematol. 2015;37:410–9.
Fan Z, Weng X, Huang G, Pan Z, Long Z, Fan Q, et al. STARD-rapid screening for the 6 most common G6PD gene mutations in the Chinese population using the amplification refractory mutation system combined with melting curve analysis. Medicine. 2018;97:e0426.
Tanphaichitr VS, Pung-amritt P, Yodthong S, Soongswang J, Mahasandana C, Suvatte V. Glucose-6-phosphate dehydrogenase deficiency in the newborn: its prevalence and relation to neonatal jaundice. Southeast Asian J Trop Med Public Health. 1995;26:137–41.
Nuchprayoon I, Sanpavat S, Nuchprayoon S. Glucose-6-phosphate dehydrogenase (G6PD) mutations in Thailand: G6PD Viangchan (871G > A) is the most common deficiency variant in the Thai population. Hum Mutat. 2002;19:185.
Laosombat V, Sattayasevana B, Janejindamai W, Viprakasit V, Shirakawa T, Nishiyama K, et al. Molecular heterogeneity of glucose-6-phosphate dehydrogenase (G6PD) variants in the south of Thailand and identification of a novel variant (G6PD Songklanagarind). Blood Cells Mol Dis. 2005;34:191–6.
Phompradit P, Kuesap J, Chaijaroenkul W, Rueangweerayut R, Hongkaew Y, Yamnuan R, et al. Prevalence and distribution of glucose-6-phosphate dehydrogenase (G6PD) variants in Thai and Burmese populations in malaria endemic areas of Thailand. Malar J. 2011;10:368.
Tuchinda S, Rucknagel DL, Na-Nakorn S, Wasi P. The Thai variant and the distribution of alleles of 6-phosphogluconate dehydrogenase and the distribution of glucose 6-phosphate dehydrogenase deficiency in Thailand. Biochem Genet. 1968;2:253–64.
Chamchoy K, Praoparotai A, Pakparnich P, Sudsumrit S, Swangsri T, Chamnanchanunt S, et al. The integrity and stability of specimens under different storage conditions for glucose-6-phosphate dehydrogenase deficiency screening using WST-8. Acta Trop. 2021;217:105864.
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.
World Health Organization. Technical specifications series for submission to WHO prequalification–diagnostic assessment: in vitro diagnostic medical devices to identify glucose-6-phosphate dehydrogenase (G6PD) activity. Geneva, World Health Organization, 2016.
Llanos-Cuentas A, Lacerda MVG, Hien TT, Velez ID, Namaik-Larp C, Chu CS, et al. Tafenoquine versus primaquine to prevent relapse of Plasmodium vivax malaria. N Engl J Med. 2019;380:229–41.
Lacerda MVG, Llanos-Cuentas A, Krudsood S, Lon C, Saunders DL, Mohammed R, et al. Single-dose tafenoquine to prevent relapse of Plasmodium vivax malaria. N Engl J Med. 2019;380:215–28.
Chen Y, Xiu W, Dong Y, Wang J, Zhao H, Su Y, et al. Mutation of glucose-6-phosphate dehydrogenase deficiency in Chinese Han children in eastern Fujian. Medicine. 2018;97:e11553.
Xia Z, Chen P, Tang N, Yan T, Zhou Y, Xiao Q, et al. Rapid detection of G6PD mutations by multicolor melting curve analysis. Mol Genet Metab. 2016;119:168–73.
Bancone G, Gilder ME, Chowwiwat N, Gornsawun G, Win E, Cho WW, et al. Prevalences of inherited red blood cell disorders in pregnant women of different ethnicities living along the Thailand-Myanmar border. Wellcome Open Res. 2017;2:72.
Commons RJ, McCarthy JS, Price RN. Tafenoquine for the radical cure and prevention of malaria: the importance of testing for G6PD deficiency. Med J Aust. 2020;212:152–3.
Chu CS, Bancone G, Nosten F, White NJ, Luzzatto L. Primaquine-induced haemolysis in females heterozygous for G6PD deficiency. Malar J. 2018;17:101.
Hu R, Lin M, Ye J, Zheng BP, Jiang LX, Zhu JJ, et al. Molecular epidemiological investigation of G6PD deficiency by a gene chip among Chinese Hakka of southern Jiangxi province. Int J Clin Exp Pathol. 2015;8:15013–8.
The authors are greatly indebted to the communities for their participation in this research.
This research was supported by Thailand Science Research and Innovation (RTA6280006), Mahidol University (Goal-oriented Research Project Grant) and Newton Fund Researcher Links Travel Grant. PP acknowledges Thailand Graduate Institute of Science and Technology (TGIST).
Ethics approval and consent to participate
Ethical approval for the study was provided by the Committee on Human Rights Related to Research Involving Human Subjects, Faculty of Medicine Ramathibodi Hospital, Mahidol University, Bangkok, Thailand (approval number MURA2018/252) and the Human Ethics Committee, Faculty of Tropical Medicine, Mahidol University (approval number MUTM 2019-016-01).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Boonyuen, U., Songdej, D., Tanyaratsrisakul, S. et al. Glucose-6-phosphate dehydrogenase mutations in malaria endemic area of Thailand by multiplexed high‐resolution melting curve analysis. Malar J 20, 194 (2021). https://doi.org/10.1186/s12936-021-03731-0
- G6PD deficiency
- G6PD mutations
- High-resolution melting
- G6PD enzyme activity
- G6PD genotype