Skip to main content

G6PD deficiency in malaria endemic areas of Nepal

Abstract

Background

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is currently a threat to malaria elimination due to risk of primaquine-induced haemolysis in G6PD deficient individuals. The World Health Organization (WHO) recommends G6PD screening before providing primaquine as a radical treatment against vivax malaria. However, evidence regarding the prevalence and causing mutations of G6PD deficiency in Nepal is scarce.

Methods

A cross-sectional, population-based, prevalence study was carried out from May to October 2016 in 12 malaria-endemic districts of Nepal. The screening survey included 4067 participants whose G6PD status was determined by G6PD Care Startâ„¢ rapid diagnostic test and genotyping.

Results

The prevalence of G6PD deficiency at the national level was 3.5% (4.1% among males and 2.1% among females). When analysed according to ethnic groups, G6PD deficiency was highest among the Janajati (6.2% overall, 17.6% in Mahatto, 7.7% in Chaudhary and 7.5% in Tharu) and low among Brahman and Chhetri (1.3%). District-wise, prevalence was highest in Banke (7.6%) and Chitwan (6.6%). Coimbra mutation (592 C>T) was found among 75.5% of the G6PD-deficient samples analysed and Mahidol (487 G>A) and Mediterranean (563 C>T) mutations were found in equal proportions in the remaining 24.5%. There was no specific geographic or ethnic distribution for the three mutations.

Conclusions

This study has identified populations with moderate to high prevalence of G6PD deficiency which provides strong evidence supporting the WHO recommendations to screen G6PD deficiency at health facility level before the use of primaquine-based radical curative regimen for Plasmodium vivax.

Background

Malaria remains a global problem with an estimated 228 million cases occurring globally and 405,000 deaths in 2018. It is still endemic in 91 countries with about half of the world’s population at risk, particularly those living in lower-income countries [1]. Ending the epidemic of malaria, AIDS, tuberculosis and other neglected tropical diseases by 2030 is related to target 3.3 of the Sustainable Development Goals. The Asia–Pacific countries, including Nepal, have shown their commitment to regional elimination of malaria by 2030 [2]. Although the global incidence of malaria has decreased by 41% between 2000 and 2015, fewer than half of the countries are on track to achieve the global strategy target of 40% case and mortality reduction by 2020 [1, 3, 4].

The presence of G6PD deficiency among the population is a huge challenge for malaria treatment due to the associated haemolytic risk with radical treatment using 8-aminoquinolines. G6PD deficiency is the most common human enzyme defect and is mostly asymptomatic. However, G6PD-deficient individuals may suffer from acute haemolytic anaemia (AHA) when exposed to several oxidative triggers, such as fava beans, infections or anti-malarial drugs such as primaquine. AHA is the most common manifestation of this deficiency associated with intra and extra-vascular haemolysis, jaundice and haemoglobinuria with acute renal failure as the most severe outcome [5,6,7,8,9,10,11]. Factors affecting individual haemolytic risk include drug dose and drug metabolization (dependent from CYP2D6 in primaquine), G6PD phenotype and genotype, and haematologic picture, particularly red blood cell age distribution [2, 10, 12]. At population level, males are historically considered more susceptible to haemolysis because the deficient phenotype is more prevalent among them, but recent evidence shows that heterozygous women with intermediate G6PD phenotypes are also at risk of severe haemolysis [13, 14].

About 400 million people are estimated to be G6PD deficient worldwide [5, 15,16,17,18], with a geographical distribution similar to that of malaria [5, 6, 19] but also with high variability within countries and regions. The affected population in malaria-endemic countries is estimated to be 220 million males and 133 million females [12, 16]. Sub-Saharan Africa has the highest average estimated prevalence (7.5%) followed by Middle East (6.0%) and Asia (4.7%). The prevalence is lowest for Pacific (2.9%), Americas (3.4%) and Europe (3.9%) [17].

Studies on G6PD deficiency among a small number of individuals of Nepalese origin were conducted as early as 1960s in Calcutta [20]. Studies conducted until the early 2000s, among Nepalese healthy individuals, found extreme variability in G6PD prevalence among ethnic groups and geographic regions; Tharus living in Banke and Bardiya district showed an exceptionally high prevalence (14.0%) compared to most of other ethnicities and regions. The mutation responsible for the deficiency was identified in one study to be the Mediterranean variant [21,22,23,24]. In a recent and first population-based study conducted in six high-risk malaria-endemic districts in Nepal in 2013 [25], G6PD deficiency screening was carried out among 1341 individuals using two screening tests: BinaxNOW and CareStart™. The prevalence of G6PD deficiency was found to be 7.2% according to BinaxNOW and 6.0% by CareStart™. A higher prevalence was observed among males of Rajbanshi and Tharu ethnic groups, and in Jhapa and Morang districts.

Nepal has a long history of malaria control programmes starting from 1954 [26] and its commitment to malaria elimination has been reflected in Nepal Malaria Strategic Plan (2014–2025) which aims for malaria-free Nepal by 2026 [26, 27]. The high risk areas of malaria in the country are the foothills with river belts, forest fringe areas in the Tarai, hill river valleys and inner Tarai areas [26]; for recent maps of malaria endemicity see [28]. In Nepal, malaria is decreasing sharply (from 12,750 cases in 2002 to 991 in 2016) and it is caused predominantly by Plasmodium vivax [26, 29]. However, P. vivax cases have been recorded in non-endemic areas representing a challenge for malaria elimination efforts. The Epidemiology and Disease Control Division (EDCD) has adopted the World Health Organization (WHO) policy to screen malaria patients for G6PD deficiency in its National Malaria treatment protocol published in 2015 [30].

Understanding the burden of G6PD deficiency and identifying vulnerable populations is crucial for the success of the malaria control programme in Nepal. In this regard, this study aims to estimate the prevalence of G6PD deficiency and to describe the underlying mutations in areas where malaria transmission is endemic. It is anticipated that the study findings will inform strategic plans for the malaria elimination programme in Nepal.

Methods

A cross-sectional prevalence survey was conducted among 4067 study participants of 12 malaria-endemic districts of Nepal. Data were collected from May 2016 to October 2016. The study population comprised both male and female participants above 6 years of age and residing in the study districts.

Study sites

The study districts were Banke, Bardiya, Chitwan, Dadeldhura, Kailali, Kanchanpur, Kapilvastu, Makwanpur, Nawalparasi, Rautahat, Sindhuli, and Surkhet. These districts were selected according to malaria endemicity data based on the malaria micro-stratification report 2013 (http://edcd.gov.np/resource-detail/report-on-malaria-microstratification-2013).

Sampling method

Two-stage, cluster-sampling method was used for drawing participants in the sample. Among the 1254 Village Development Committees (VDCs), 54 were identified to be at high risk of malaria based on the malaria micro-stratification report 2013. These VDCs were further divided into clusters and 30 of them were selected using systematic random sampling with probability inclusion proportional to size (PPS) method from the VDC.

The field teams visited each selected cluster to prepare a list of households for the study and 132 households were selected by systematic random sampling method from each cluster. From this list, a male participant was selected from the first and second household and a female participant was selected from the third household; this was repeated for every block of three successive households resulting in the selection of a total of (30 × 132) = 3960 participants with a male to female ratio of 2:1. The total final sample size of the study was 4067 due to additional samples taken in some sites.

Data and samples collection methods

A structured questionnaire was used for data collection, including socio-demographic characteristics. Interviews were conducted in Nepali language by researchers of the same gender of interviewed participants. Screening of G6PD status was done on capillary blood using CareStart™ (Access Bio. Inc., New Jersey, USA). CareStart™ retraining was supported by investigators of the Shoklo Malaria Research Unit (SMRU) prior to survey activities. The test was run following manufacturer’s instructions. Test interpretation was based on colour change in the reading window, with a pink purple colour interpreted as ‘G6PD normal’ result and a faint or no colour development interpreted as ‘G6PD deficient’ result. Blood for human genotyping was spotted on filter papers, dried, and stored at − 20 °C.

Study variables

The study variables were socio-demographic variables (age, gender, ethnicity, permanent residence), previous history of malaria and G6PD screening results. Since the ethnicity and caste classification of Nepal is extremely complex [31], both the ethnicity and the Ministry of Health larger group classification were used in the analysis.

Data were collected in the field by district level laboratory staff, peripheral level laboratory and health workers supervised by District Public Health Office (DPHO) team and EDCD team. Data were quality controlled to assure uniformity across all sites.

Laboratory analysis

Dried blood spots of subjects who tested deficient were shipped on ice packs to SMRU in Mae Sot (Thailand). DNA was extracted by mini-column (Favorgen Biotech, Taiwan) or saponin-chelex method for PCR–RFLP and gene sequencing from exon 2 to exon 13 of G6PD gene. Primer list and PCR conditions are reported in Additional file 1: Table S1 and S2. Amplified fragments were sent to external sequencing platforms and results were analysed by CLC Main Workbench 6 (Qiagen Bioinformatics). Based on published literature at the time of the study, Mediterranean (563 C>T) and Coimbra (592 C>T) mutations were expected to be found in the population, so all samples were analysed first for exon 6 which contains both variants and Mahidol (487 G>A). When samples were found to harbour no variants in exon 6, the other exons were analysed.

Data management and statistical analysis

Survey data were entered in Microsoft Excel and then transported into IBM SPSS version 20. Double entry was done to ensure accuracy. Prevalence of G6PD deficiency was calculated from the pooled number of subjects with deficient or suspected deficient phenotype over the total of subjects included. Allelic frequencies were calculated as the number of mutated alleles (one in hemizygous males, one in heterozygous females and two in homozygous mutated females) over the total number of alleles (where females contribute two alleles each and males one allele each). Allelic frequencies among groups were compared by Chi square test. A P < 0.05 was considered significant.

Results

The screening of G6PD deficiency was done among 4067 participants. Ten participants whose data were missing were removed from further analysis, thus 4057 participants were included in the final analysis.

Demographic characteristics of the participants

The demographic characteristics of the participants are summarized in Table 1. As per study design, almost two-thirds of the participants (2564, 63.2%) were male. The age of the participants ranged from 6 years to 90 years with a mean (standard deviation, SD) of 36.2 (18.0) years. Most of the participants were from the age group 20–39 years (36.0%) followed by 40–59 years (29.0%). The age of participants might not have been representative of the general population as participants were those present in the house at the time of study visit. The Tharu ethnic group was the most prevalent representing alone 24.7% of the whole population. Following the Ministry of Health classification, almost half of the population analysed belonged to the Janajati ethnic group (48%) followed by Brahmin/Chhetri (30.3%). One-third of participants were from Kailali district (30.2%) followed by Nawalparasi (13.6%) and Bardiya (10.3%). A total of 7.7% of the study participants mentioned that they had previous history of malaria.

Table 1 Demographic characteristics of participants, n (%)

G6PD phenotypes

The G6PD phenotypic status among the participants is presented in Tables 2, 3 and 4. The overall prevalence of G6PD deficiency in the study population was 3.5%, with a higher prevalence among males (4.4%) compared to females (2.1%, P < 0.01) as expected by an X-linked condition. This estimate includes 16 tests (0.4%) that were interpreted as suspected or borderline deficient. The prevalence of G6PD deficiency when classified by ethnic groups (Table 3) was highest among Janajati MOH ethnic group (overall 6.2%); in particular, in the Mahatto group (17.6%) it was significantly higher compared to any other group and it was similar in the Chaudhary (7.7%) and Tharu (7.5%), lower than the Mahatto but significantly higher compared to all the remaining groups.

Table 2 Distribution of G6PD phenotypes by gender
Table 3 G6PD phenotypes by ethnic group
Table 4 Prevalence of G6PD deficiency by district

Prevalence of G6PD deficiency was highest in Banke district (7.6%) followed by Chitwan (6.6%), Nawalparasi and Bardiya (4.3%), and Kailali (4.0%) (Table 4 and Fig. 1). The other districts showed prevalence lower than 4% and nil in Makwanpur and Dadeldhura. When analysed by larger developmental region (Additional file 1: Table S3), G6PD deficiency prevalence ranged from 2.9% to 4.2% and was not different among the regions.

Fig. 1
figure 1

Prevalence of G6PD deficiency in malaria-endemic districts of Nepal

Prevalence of G6PD deficiency across districts within the same ethnic group was highly variable; the most numerous ethnic group Tharu had a prevalence ranging from 1.6% in Kanchanpur to 23.4% in Banke district. Chhetri and Brahman had a prevalence ranging from nil to 7.1% and 9.1%, respectively, in different districts (Additional file 1: Table S4 and Fig. S1).

G6PD genotypes

Out of 143 deficient samples screened in the field, 118 had an available sample to perform genotyping. Overall the DNA quality of the samples analysed was acceptable for both standard PCR–RFLP protocol and sequencing of short-length fragments. Long-length fragments such as exon 9 and 13 required instead more work and numerous attempts of amplification and sequencing. A total of three samples could not be amplified in these two exons and one sample could not be amplified at any exon.

Among the 117 analysed samples the genotypes were distributed as shown in Table 5. The 19 samples that did not have Coimbra, Mahidol or Mediterranean mutations in exon 6 were analysed in the remaining exons. No mutations were found in the other exons and samples were therefore considered wild type. Coimbra accounted for the vast majority of mutations (74.8% of all mutated alleles), Mahidol and Mediterranean were found at same frequency (12.6% each).

Table 5 G6PD genotypes

As already seen in the prevalence of phenotypic deficiency, the cumulative allelic frequency of the three mutations was higher in the Mahatto, Chaudhary and Tharu groups who all also harboured the three mutations (Table 6). No specific geographic or ethnic distribution pattern was observed for the three mutations, with Coimbra being present in all districts where G6PD deficiency was observed (Additional file 1: Table S5) and at the higher allelic frequency in all ethnic groups with the exception of the Brahman where it was absent.

Table 6 Estimated allelic frequencies (%) by ethnic group

Performance of CareStartâ„¢ rapid diagnostic test

The performances of the CareStart™ rapid diagnostic test (RDT) are shown in Table 7. Although it was not possible to calculate sensitivity and specificity of the G6PD CareStart™ because the genotyping was carried out only on subjects with deficient phenotype, a relatively large number of RDT-deficient samples (12/101 and 5/16 suspected deficient) had a wild type genotype confirmed by full gene sequence. This could be due to failure in developing a strong colour in the reading window or problems in the quality of the tests. The CareStart™ G6PD test in its current formulation does not have a control band.

Table 7 Carestart™ RDT results by G6PD genotype

Discussion and conclusion

In the current study, the overall prevalence of G6PD deficiency was 3.5% at the population level, with a wide geographical variation (from nil to 7.6%) and across ethnicities (from nil to 17.6%). This prevalence is consistent with the worldwide prevalence (4.9%), and the prevalence in Asia (4.7%) and across much of the Indian Sub-continent (2–5%) [5, 12]. In a similar survey in Nepal in 2013, the prevalence of G6PD deficiency was higher, especially in the most eastern districts of the Terai regions [25] among Rajbanshi (11.7%, classified as Janajati by MOH) and Tharus (5.6%); prevalence was only 1.4% among Brahmins and Chhetris. A similar prevalence among Brahmins and Chhetris was found in the current study and higher prevalence in Tharus, Mahatto and Chaudhary of the Janajati group.

The sensitivity and specificity of the G6PD CareStartâ„¢ RDT could not be formally assessed in this study; the inconsistencies observed in this context hint to a more general consideration on the necessity to use fully validated and quality-controlled G6PD tests for malaria elimination strategies to be effective.

The three G6PD variants found had a broad distribution geographically and among ethnicities in the studied districts. Interestingly, Mahidol is a mutation found mainly in the most western countries of Southeast Asia (Myanmar and Thailand [32, 33]) and Coimbra is a mutation found at low frequencies predominantly in the Indian region of Uttar Pradesh [34], in Tharus living in India [35] and rarely in Myanmar [36] and Malaysia [37]. The relatively high allelic frequency of Coimbra variant found in the Mahatto (8.5%) and in the Tharu (3.9%) in the current study seems to exceed all previous published reports. The Mediterranean variant is found in Europe and Asia, being the most frequent G6PD mutation in India among the caste groups [38, 39]. Overall the genetic data of the most numerous ethnic group in the current study (the Tharus) align to previous evidence of their dual Southeast Asian and Indian ancestry [40, 41].

Primaquine and tafenoquine are the only currently available drugs effective against hypnozoites of P. vivax [12]. The adverse effects of 8-aminoquinolines among G6PD-deficient individuals are widely documented [2, 11, 42] although more evidence is still needed on haemolytic risk associated with different regimens and populations [12, 43]. The WHO has recommended G6PD screening among populations in regions where G6PD deficiency prevalence is relatively high, and P. vivax radical cure treatment regimens be based on the level of G6PD activity [44, 45].

The findings of the current study have therefore very clear programmatic and policy-making implications. They show that G6PD deficiency is present broadly in the Terai regions of Nepal and is caused by mutations associated to described primaquine-induced haemolysis. Some ethnic groups are more affected but the prevalence of deficiency within the same ethnicity varies widely in different districts. Since the risk of haemolysis due to G6PD deficiency cannot be excluded on the basis of the geographic region or ethnicity alone, G6PD testing will be always necessary before deployment of primaquine. The EDCD in Nepal has provided directives for G6PD screening among the population at health facility level [46]. However, its implementation in the routine health system is still not in practice. Further evidence might also be required in assessing the cost-effectiveness of introducing screening kits. In the context of the available evidence and following the WHO recommendation, stakeholders, including government and academia are invited to: (i) conduct further studies in malaria-endemic as well as non-endemic districts to determine the prevalence, genetic variant of G6PD and the risk population; (ii) ensure the availability of G6PD screening kits; (iii) demand training of health workers to ensure that the revised malaria treatment protocol is followed across primary health facilities; (iv) develop a pharmacovigilance system to report the adverse events with the use of primaquine; and, (v) manage the adverse events associated with the drug through proper clinical training of health workers and availability of referral system.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

References

  1. WHO. World malaria report 2016. Geneva: World Health Organization; 2017.

    Google Scholar 

  2. Thriemer K, Ley B, Bobogare A, Dysoley L, Alam MS, Pasaribu AP, et al. Challenges for achieving safe and effective radical cure of Plasmodium vivax: a round table discussion of the APMEN Vivax Working Group. Malar J. 2017;16:141.

    PubMed  PubMed Central  Google Scholar 

  3. United Nations. The sustainable development goals report 2017. New York: United Nations; 2017.

    Google Scholar 

  4. WHO. Global technical strategy for malaria 2016–2030. Geneva: World Health Organization; 2015.

    Google Scholar 

  5. Cappellini MD, Fiorelli G. Glucose-6-phosphate dehydrogenase deficiency. Lancet. 2008;371:64–74.

    CAS  PubMed  Google Scholar 

  6. Luzzatto L, Nannelli C, Notaro R. Glucose-6-phosphate dehydrogenase deficiency. Hematol Oncol Clin North Am. 2016;30:373–93.

    PubMed  Google Scholar 

  7. Manganelli G, Masullo U, Passarelli S, Filosa S. Glucose-6-phosphate dehydrogenase deficiency: disadvantages and possible benefits. Cardiovasc Hematol Disord: Drug Targets. 2013;13:73–82.

    CAS  Google Scholar 

  8. Monteiro WM, Franca GP, Melo GC, Queiroz AL, Brito M, Peixoto HM, et al. Clinical complications of G6PD deficiency in Latin American and Caribbean populations: systematic review and implications for malaria elimination programmes. Malar J. 2014;13:70.

    PubMed  PubMed Central  Google Scholar 

  9. Luzzatto L, Seneca E. G6PD deficiency: a classic example of pharmacogenetics with on-going clinical implications. Br J Haematol. 2014;164:469–80.

    CAS  PubMed  PubMed Central  Google Scholar 

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

    CAS  PubMed  Google Scholar 

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

    PubMed  PubMed Central  Google Scholar 

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

    PubMed  Google Scholar 

  13. 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:e1002224.

    PubMed  PubMed Central  Google Scholar 

  14. Chu CS, Bancone G, Nosten F, White NJ, Luzzatto L. Primaquine-induced haemolysis in females heterozygous for G6PD deficiency. Malar J. 2018;17:101.

    PubMed  PubMed Central  Google Scholar 

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

    Google Scholar 

  16. 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;9:e1001339.

    CAS  PubMed  PubMed Central  Google Scholar 

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

    CAS  PubMed  Google Scholar 

  18. Peters AL, Van Noorden CJ. Glucose-6-phosphate dehydrogenase deficiency and malaria: cytochemical detection of heterozygous G6PD deficiency in women. J Histochem Cytochem. 2009;57:1003–11.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Galappaththy GN, Tharyan P, Kirubakaran R. Primaquine for preventing relapse in people with Plasmodium vivax malaria treated with chloroquine. Cochrane Database Syst Rev. 2013;10:CD004389.

    Google Scholar 

  20. Chatterjea J. Haemoglobinopathies, glucose-6-phosphate dehydrogenase deficiency and allied problems in the Indian subcontinent. Bull World Health Organ. 1966;35:837.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. Henderson A, Baronne A. Glucose-6-phosphate dehydrogenase deficiency in healthy Nepalese patients. Trans R Soc Trop Med Hyg. 1987;81:543.

    CAS  PubMed  Google Scholar 

  22. Modiano G, Morpurgo G, Terrenato L, Novelletto A, Di Rienzo A, Colombo B, et al. Protection against malaria morbidity: near-fixation of the alpha-thalassemia gene in a Nepalese population. Am J Hum Genet. 1991;48:390–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Matsuoka H, Jichun W, Hirai M, Yoshida S, Arai M, Ishii A, et al. Two cases of glucose-6-phophate dehydrogenase-deficient Nepalese belonging to the G6PD Mediterranean-type, not India-Pakistan sub-type but Mediterranean-Middle East sub-type. J Hum Genet. 2003;48:275–7.

    CAS  PubMed  Google Scholar 

  24. Suzuki A, Hamano S, Shirakawa T, Watanabe K, Endo T, Sharma S, et al. The distribution of hereditary erythrocytic disorders associated with malaria, in a lowland area of Nepal: a micro-epidemiological study. Ann Trop Med Parasitol. 2007;101:113–22.

    CAS  PubMed  Google Scholar 

  25. Ghimire P, Singh N, Ortega L, Rijal KR, Adhikari B, Thakur GD, et al. Glucose-6-phosphate dehydrogenase deficiency in people living in malaria endemic districts of Nepal. Malar J. 2017;16:214.

    PubMed  PubMed Central  Google Scholar 

  26. MoH. Annual report 2015/2016. Kathmandu: Department of Health Services, Ministry of Health; 2016.

    Google Scholar 

  27. EDCD. Nepal malaria strategic plan 2014–2025. Kathmandu: EDCD; 2014.

    Google Scholar 

  28. Rijal KR, Adhikari B, Adhikari N, Dumre SP, Banjara MS, Shrestha UT, et al. Micro-stratification of malaria risk in Nepal: implications for malaria control and elimination. Trop Med Health. 2019;47:21.

    PubMed  PubMed Central  Google Scholar 

  29. Dhimal M, Ahrens B, Kuch U. Malaria control in Nepal 1963–2012: challenges on the path towards elimination. Malar J. 2014;13:241.

    PubMed  PubMed Central  Google Scholar 

  30. EDCD. National malaria treatment protocol. Kathmandu: EDCD; 2015.

    Google Scholar 

  31. Nepal Health Sector Support Programme. Nepal Population Report 2011. http://www.nhssp.org.np/NHSSP_Archives/monitoring/Nepal_Population_Report_2011.pdf.

  32. Bancone G, Menard D, Khim N, Kim S, Canier L, Nguong C, et al. Molecular characterization and mapping of glucose-6-phosphate dehydrogenase (G6PD) mutations in the Greater Mekong Subregion. Malar J. 2019;18:20.

    PubMed  PubMed Central  Google Scholar 

  33. Matsuoka H, Wang J, Hirai M, Arai M, Yoshida S, Kobayashi T, et al. Glucose-6-phosphate dehydrogenase (G6PD) mutations in Myanmar: G6PD Mahidol (487 G>A) is the most common variant in the Myanmar population. J Hum Genet. 2004;49:544–7.

    CAS  PubMed  Google Scholar 

  34. Tripathi P, Agarwal S, Muthuswamy S. Prevalence and genetic characterization of glucose-6-phosphate dehydrogenase deficiency in anemic subjects from Uttar Pradesh, India. J Pediatr Genet. 2019;8:47–53.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Sarkar S, Biswas NK, Dey B, Mukhopadhyay D, Majumder PP. A large, systematic molecular-genetic study of G6PD in Indian populations identifies a new non-synonymous variant and supports recent positive selection. Infect Genet Evol. 2010;10:1228–36.

    CAS  PubMed  Google Scholar 

  36. Nuchprayoon I, Louicharoen C, Charoenvej W. Glucose-6-phosphate dehydrogenase mutations in Mon and Burmese of southern Myanmar. J Hum Genet. 2008;53:48–54.

    CAS  PubMed  Google Scholar 

  37. Iwai K, Hirono A, Matsuoka H, Kawamoto F, Horie T, Lin K, et al. Distribution of glucose-6-phosphate dehydrogenase mutations in Southeast Asia. Hum Genet. 2001;108:445–9.

    CAS  PubMed  Google Scholar 

  38. Sukumar S, Mukherjee MB, Colah RB, Mohanty D. Molecular basis of G6PD deficiency in India. Blood Cells Mol Dis. 2004;33:141–5.

    CAS  PubMed  Google Scholar 

  39. Mohanty D, Mukherjee MB, Colah RB. Glucose-6-phosphate dehydrogenase deficiency in India. Indian J Pediatr. 2004;71:525–9.

    PubMed  Google Scholar 

  40. Chaubey G, Singh M, Crivellaro F, Tamang R, Nandan A, Singh K, et al. Unravelling the distinct strains of Tharu ancestry. Eur J Hum Genet. 2014;22:1404–12.

    PubMed  PubMed Central  Google Scholar 

  41. Passarino G, Semino O, Pepe G, Shrestha SL, Modiano G, Santachiara Benerecetti AS. MtDNA polymorphisms among Tharus of eastern Terai (Nepal). Gene Geogr. 1992;6:139–47.

    CAS  PubMed  Google Scholar 

  42. Rajapakse S, Rodrigo C, Fernando SD. Tafenoquine for preventing relapse in people with Plasmodium vivax malaria. Cochrane Database Syst Rev. 2015;4:CD010458.

    Google Scholar 

  43. White NJ, Qiao LG, Qi G, Luzzatto L. Rationale for recommending a lower dose of primaquine as a Plasmodium falciparum gametocytocide in populations where G6PD deficiency is common. Malar J. 2012;11:418.

    PubMed  PubMed Central  Google Scholar 

  44. Bosman A, Cunningham J. Proposal for an Evidence Review Group (ERG) on G6PD testing to support increased access to primaquine for radical cure of Plasmodium vivax and for malaria chemoprophylaxis. Geneva: World Health Organization; 2013.

    Google Scholar 

  45. WHO. Guidelines for the treatment of malaria. Geneva: World Health Organization; 2015.

    Google Scholar 

  46. EDCD. Guidelines for case based investigation. Kathmandu: EDCD, Government of Nepal; 2017.

    Google Scholar 

Download references

Acknowledgements

The authors wish to acknowledge all participants of the study and the staff who supported field and laboratory activities. We would like to thank Mr. Saw Myo Chit Minh (SMRU) for producing the map.

Funding

The study was funded by the Government of Nepal, Ministry of Health and population, Department of Health Services, Epidemiology and Disease Control Division.

Author information

Authors and Affiliations

Authors

Contributions

BM, BKL, BPA: supervised project writing, implementation and data analysis. ST, KRA, BB: contributed to project writing, implemented field activities and data collection. PK: conceived the project, performed preliminary analysis and wrote the first draft of manuscript. SN: contributed to project writing and data analysis. SNJ, SA, SI, MK, DK, SB, RKM, AC, PG: collected data and conducted field activities. RGM, RB: performed data analysis and manuscript writing. GG, PP: contributed to field activities and samples analysis. GB: contributed to samples analysis; performed data analysis and interpretation; manuscript writing. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Bijay Bajracharya or Germana Bancone.

Ethics declarations

Ethics approval and consent to participate

The study was approved by Nepal Health Research Council (NHRC) with the (Reference number: 23/2016). All study participants gave consent to participate in the study; they were provided with pre-test counseling and results of testing. Blood testing was conducted by qualified lab personnel.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1.

 Tables S1–S5 and Figure S1.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Marasini, B., Lal, B.K., Thapa, S. et al. G6PD deficiency in malaria endemic areas of Nepal. Malar J 19, 287 (2020). https://doi.org/10.1186/s12936-020-03359-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-020-03359-6

Keywords