- Open Access
Prevalence and genetic variants of G6PD deficiency among two Malagasy populations living in Plasmodium vivax-endemic areas
Malaria Journal volume 16, Article number: 139 (2017)
The prevalence and variants of G6PD deficiency in the Plasmodium vivax-endemic zones of Madagascar remain unknown. The admixed African-Austronesian origins of the Malagasy population make it probable that a heterogeneous mix of genetic variants with a spectrum of clinical severity will be circulating. This would have implications for the widespread use of P. vivax radical cure therapy. Two study populations in the P. vivax-endemic western foothills region of Madagascar were selected for G6PD screening. Both the qualitative fluorescent spot test and G6PD genotyping were used to screen all participants.
A total of 365 unrelated male volunteers from the Tsiroanomandidy, Mandoto, and Miandrivazo districts of Madagascar were screened and 12.9% were found to be phenotypically G6PD deficient. Full gene sequencing of 95 samples identified 16 single nucleotide polymorphisms, which were integrated into a genotyping assay. Genotyping (n = 291) found one individual diagnosed with the severe G6PD Mediterranean C563T mutation, while the remaining G6PD deficient samples had mutations of African origin, G6PD A- and G6PD A.
Deployment of P. vivax radical cure in Madagascar must be considerate of the risks presented by the observed prevalence of G6PDd prevalence. The potential morbidity associated with cumulative episodes of P. vivax clinical relapses requires a strategy for increasing access to safe radical cure. The observed dominance of African G6PDd haplotypes is surprising given the known mixed African-Austronesian origins of the Malagasy population; more widespread surveying of G6PDd epidemiology across the island would be required to characterize the distribution of G6PD haplotypes across Madagascar.
Malaria remains a major public health burden in Madagascar, with the entire Malagasy population living at risk of infection and suffering an estimated 2.4 million (1.5–4.0 million) cases of malaria in 2015 . Temporal trends in reported health facility data suggest significant increases in malaria burden from 2010 to 2015 . The last major national review of the health status of the Malagasy population reported malaria to be the second most common cause of death in under-5 year olds in district hospitals, and the fourth leading cause of all outpatient consultations across all age groups . The island is biogeographically diverse, with varied malaria transmission patterns across different parts of the island ranging from epidemic, to highly seasonal, to sustained year-round transmission . Malaria in Madagascar is predominantly caused by Plasmodium falciparum, though the other human malaria species are also present [4, 5] and likely under-reported. Plasmodium vivax transmission is reported from sites across the island, but most commonly in the fringe zone between the west coast and the central highlands [6, 7]. Parasites in these regions were found to be infecting and causing clinical P. vivax malaria in Duffy negative individuals, the blood group previously considered refractory to infection due to the absence of Duffy antigen expression on the surface of red blood cells [6, 8, 9]. A cross-sectional study using PCR-based molecular diagnosis in this western highland fringe region found an 8% prevalence of P. vivax, the same as that for P. falciparum .
Madagascar’s current national strategic plan (2013–2017) calls for dual applications of the anti-malarial primaquine: (i) in low dose regimens for blocking P. falciparum transmission in epidemic-prone zones and in low transmission areas targeting pre-elimination status, and (ii) in its 14-day formulation as P. vivax radical cure . However, neither application of primaquine is deployed, in large part due to uncertainties relating to safety in glucose-6-phosphate dehydrogenase deficient (G6PDd) patients (National Malaria Control Programme, Pers. Comm.) .
G6PD deficiency is highly diverse, with at least 217 genetic variants expressing a spectrum of enzyme phenotypes of differing susceptibilities to toxic substances including primaquine . At the global scale, distinct regional patterns of G6PDd epidemiology emerge . Given the contemporary Malagasy population’s predominantly African and Austronesian origins , it could be anticipated that a unique combination of alleles originating from these regions may be present. G6PD deficiency tends to be more common but less clinically severe and heterogeneous in African populations than in Asian . For instance, 27 African countries are predicted to have mean frequencies of >10% in males, while prevalence is estimated at 7.1% in Indonesian males (50% CI 5.3–9.4%) . A single G6PDd allelic variant predominates in African communities (G6PD A-202A/376G) which is generally less severe than most of the diverse variants found in Asian populations (including, among many, G6PD Viangchan 871A, G6PD Vanua Lava 383C and G6PD Chatham 1003A) [17, 18].
Recognizing the uncertainty around the potentially complex G6PD allele inheritance in Madagascar, this study therefore aims to begin strengthening the evidence base of epidemiological data available to the Malagasy National Malaria Control Programme (NMCP) regarding G6PD deficiency. Providing this insight will contribute to improved case management by increasing access to safe primaquine treatment, notably in its higher dosing for P. vivax radical cure. This study directly addresses a stated operational research priority in the country’s National Strategic Plan (2013–2017). To date, the only published survey of the Malagasy population was from the east coast city of Manakara in 1966 which found 21% of males with G6PDd (n tested = 70) .
This study protocol was approved by the ethical review panels of University Hospitals Case Medical Center, Cleveland, Ohio, USA and the National Institutes of Health, USA (DMID Protocol #13-0067), and the Ministry of Public Health, Madagascar (Nº099-MSANP/CE). Written informed consent was obtained from all subjects, or subject guardians, enrolled into the study.
Two populations were surveyed which were considered genetically representative of populations targeted for the two regimens of primaquine therapy (Fig. 1). First, rural communities resident in the P. vivax-endemic highland fringe ecozone of Mandoto and Miandrivazo districts were surveyed. Second, recruitment was conducted among a migrant population resident in the rural Ampasimpotsy community of Tsiroanomandidy district where P. vivax is also endemic . This migrant community was recently (within at most one generation) relocated from urban Antananarivo in the central highlands. The central highlands are an area where single-dose primaquine interventions are scheduled to prevent P. falciparum transmission and reinforce elimination efforts . This recently relocated population was therefore considered equivalent to a population sample from the city of Antananarivo itself.
To facilitate diagnosis of G6PD enzyme activity and increase the reliability of this study’s prevalence estimates, only males were surveyed. The G6PD gene’s position on the X-chromosome means that allele frequency corresponds directly to the prevalence of deficiency in males, though not in females [14, 20]. Sampling was random and included males of all age groups. To limit double sampling of genetically related individuals, a single individual per family was included. A census was available for ensuring that individuals surveyed were consistent with this collection strategy in the Ampasimpotsy region, while in Mandoto/Miandrivazo districts this relied on the enrolment interviews. Demographic information was obtained for all participants by questionnaire. Sampling took place in Mandoto and Miandrivazo districts in July 2015, and in Ampasimpotsy in June and August 2015.
G6PDd allele frequency in Madagascar has been predicted to be around 19.4% (50% CI 11.5–30.3%) . A sample size of 235 was determined from Mandoto/Miandrivazo districts to ensure 95% confidence around a 5% precision interval of the G6PDd prevalence estimate . The smaller population size and high degree of relatedness within the closed migrant community in Ampasimpotsy constrained the sample size to 135, which corresponded to approximately 90% confidence level in the estimate. The study design is summarised in Fig. 2.
Field-based blood processing
All assays used finger-prick capillary blood . First, rapid diagnosis of Plasmodium blood-stage parasitaemia performed using the CareStartTM Malaria HRP2/pLDH(Pf/PAN) Combo assay (AccessBio, Seoul, Korea), according to the manufacturer’s instructions. Second, the fluorescent spot test (catalogue number 203-A; Trinity Biotech, Wicklow, Ireland) was used for qualitative assessment of G6PD enzyme activity. Finally, 200 µL of blood was collected into K+EDTA microtainers for molecular diagnosis.
Qualitative G6PD phenotype field-based assay
The fluorescent spot test (FST) was used according to the manufacturer’s recommended protocol with minor modifications as described here. Prior to field studies, diagnostic performance was compared between the manufacturer’s protocol and the conditions anticipated during field trials using the normal, intermediate, and deficient controls provided by the supplier (catalogue numbers G6888, G5029, and G5888, respectively). First, to address the lack of electricity to power incubators in the field, controlled laboratory tests varied the incubation temperature (37 °C vs ambient temperature range in the field [19–31 °C]) and assessed the impact on the test reaction outcomes. Incubation of samples from the supplier’s controls at ambient field site temperatures instead of 37 °C resulted in all samples (n = 9) being correctly diagnosed (“normal”, “intermediate” or “deficient”) by visual qualitative assessment under UV light following incubation at the lower temperatures. Furthermore, in the field, controls were run at several time points during the day to reflect the fluctuating ambient temperatures, and diagnoses were based on comparison with the fluorescence intensities of controls performed during the same time/temperature window. Second, diagnostic test volumes of blood and reactant were experimentally reduced to 50 and 25% of the protocol-recommended volumes to allow a greater number of tests to be performed. When reaction volumes were reduced to 25% of recommended volumes, it was not possible to distinguish G6PD normal from deficient blood as no fluorescence could be seen, even from normal samples. However, reducing the volumes by 50% resulted in 100% concordance between the visual qualitative diagnoses and the supplier control classifications; visual inspections of results stemming from these modifications were performed by three independent observers and their evaluations were in complete agreement. All field studies were therefore conducted at ambient air incubation temperatures that ranged from 19 to 31 °C using 50% of recommended reaction volumes, with the comparator controls performed at least 3 times daily at different time points during the day.
Briefly, 10 µL of whole blood was added to 200 µL of G6PDH Assay Solution and mixed in a K+EDTA microtainer to suspend RBCs. From this reaction volume, 60 µL was spotted onto Whatman 3 MM filter paper following 0, 5 and 10 min of incubation; a further 15 min spot was tested to see whether the absence of stable 37 °C incubation would require longer assay reaction time. After the last time point was spotted, filter papers were air dried for a minimum of 15 min. Fluorescence of NADPH, which is proportional to G6PD activity, was evaluated under long-wave UV light (365 nm). Samples were classified as normal (moderate/strong fluorescence at 5 min, strong fluorescence at 10 min), intermediate (weak fluorescence at 5 min, weak/moderate fluorescence at 10 min), and deficient (faint/no fluorescence at 10 min) by two blinded readers (REH, TF, AR, PAZ), and any discrepancies resolved by discussion with a third adjudicator.
Whole blood collected in K+EDTA microtainers was stored at 4 °C until DNA extraction could be performed. DNA was extracted from the whole blood samples using the Qiagen QIAmp 96 DNA Blood Kit following recommended protocols with a starting blood volume of 200 µL and an elution volume of 100 µL of Buffer AE (Valencia, CA).
G6PD PCR amplification
Three primer pairs (IDT, Coralville, IA) were designed to amplify three G6PD-specific products. The 812 base pair (bp) exon 2 fragment was amplified using primer set A (Forward A: 5′-AGCATTCCCTGTTTTCCCCA-3′ and Reverse A: 5′-GTGAGACCCCAGAGGAACTC-3′); hg19 nucleotide coordinates 154545441–154546253, GenBank X55448.1 coordinates 3253–4065, GRCh38.p2 coordinates 155316129–155316941. A 3200 bp exon 3–8 fragment was amplified using primer set B (Forward B: 5′-GGTTCTACTGAATTGAATTTTCAGG-3′ and Reverse B: 5′-AGTCTATTCTGATGAACAAGCTGAGG-3′); hg19 coordinates 154533315–154536430, GenBank X55448.1 coordinates 13,075–16,188, GRCh38.p2 coordinates 155304002–155307117. Finally, a 3115 bp exon 9–13 fragment was amplified using primer set C (Forward C: 5′-GTTCATCAGAATAGACTCGAGATGG-3′ and Reverse C: 5′-AACCAGGACAAGAGGACTAAAACC-3′); hg19 coordinates 154530131–154533331, GenBank X55448.1 coordinates 16,172–19,372, GRCh38.p2 coordinates 155300818–155304018. The total length of the human G6PD gene analysed in our studies included 7111 bp (812 bp inclusive of exon 2 and flanking intronic sequence [hg19 nucleotide coordinates 154545441–154546253]; 6299 bp inclusive of exons 3–13 and associated intronic sequence [154530131–154536430]).
For exon 2-specific amplification, PCR was performed using the Promega GoTaq® Flexi DNA Polymerase Kit (Madison, WI) in a final reaction mixture of 50 μL; containing 3 μL of genomic DNA, 0.1 μM of each primer, 200 μM of each dNTP, 1.0 mM of MgCl2, and 1.25 units of GoTaq DNA polymerase. Amplification conditions were 95 °C for 2 min, followed by 30 cycles of 95 °C for 45 s, 59 °C for 45 s, and 72 °C for 45 s and a final extension at 72 °C for 5 min. For fragments B (exons 3–8) and C (exons 9–13), PCR amplification was performed using the Roche Expand Long Range dNTPack Kit (Indianapolis, IN) in a final reaction mixture of 50 μL; containing 3 μL of genomic DNA, 0.3 μM of each primer (primer set B) or 0.45 μM of each primer (primer set C), 500 μM of each dNTP, 2.5 mM of MgCl2, 1.25% DMSO and 3.5 units of Expand Long Range enzyme mix. Amplification conditions were 92 °C for 2 min, followed by 10 cycles of 92 °C for 10 s, 55 °C for 15 s, and 68 °C for 3 min; followed by 25 cycles of 92 °C for 10 s, 55 °C for 15 s, and 68 °C for 3 min increasing 20 s every cycle and a final extension at 68 °C for 7 min.
G6PD gene sequencing
PCR products were sheared on a Covaris machine for an estimated size of 300 bp. Next, 95 samples were prepared for sequencing according to the Nextera XT DNA Library Preparation Kit (Illumina, Inc., San Diego, CA). Samples were then sequenced on an Illumina MiSeq multiplex run for 2 × 300 bp paired-end reads.
G6PD sequence data analysis
Illumina MiSeq-generated sequences were first evaluated through FastQC for initial quality control. Sequence reads were trimmed for quality (Q > 20) and adapter sequences using TrimGalore! Those reads that passed the quality filter were aligned to the G6PD locus on human reference genome GRCh38 (X:154,526,391–154,552,572) using the Bowtie 2 alignment algorithm . Genotype calls at each base position were determined using mpileup (SAMtools)  and custom Perl scripts. Only positions with a coverage greater than 50 reads were considered. Genotypes were then determined based on reference allele frequency, where reference allele frequencies greater than 0.7 were determined to be hemizygous reference, frequencies less than 0.3 were determined to be hemizygous alternate.
G6PD genotyping focused on four single nucleotide polymorphism positions identified through Illumina sequence analysis (202G → A, 376A → G, 563C → T and 1311C → T; nucleotide reference information based on wild-type cDNA sequence; GenBank X03674).
Following PCR amplification, products were further processed by a ligation detection reaction (LDR) that has been described in detail for a variety of additional studies [6, 25–27]. This LDR was performed in a reaction mixture (15 μL) containing 20 mM Tris–HCl buffer (pH 7.6), 25 mM potassium acetate, 10 mM magnesium acetate, 1 mM NAD+ , 10 mM DTT, 0.1% Triton X-100, 13 nM each LDR probe, 1 μL of PCR product, and 2 units of Taq DNA ligase (New England BioLabs, Ipswich, MA). LDR probes consisted of eight allele-specific oligonucleotides and three fluorescently labelled conserved-sequence oligonucleotides. Each allele-specific ‘classification’ probe contained a unique 5′TAG sequence for further hybridization with complementary sequence oligonucleotides bound to Luminex FlexMAP fluorescent microspheres (Luminex Corp., Austin, TX). The conserved-sequence ‘reporter’ probes were phosphorylated at the 5′ end and biotinylated at the 3′ end. Sequences of the classification and reporter probes are provided in Table 1.
All post-PCR LDR-FMA methods include the same basic three-step procedure: (i) ligation of specific oligonucleotides to target single- or multiple-nucleotide polymorphisms, (ii) FlexMAP microsphere and streptavidin-R-phycoerythrin (SA:PE) labelling of allele-specific ligation products, and (iii) detection and analysis of the specific fluorescent signals using the BioPlex suspension array system and Bio-Plex Manager analytical software (Bio-Rad Laboratories, Hercules, CA).
Reaction mixtures were initially heated for 2 min at 95 °C, followed by 32 cycles of 95 °C for 15 s and 58 °C for 2 min (annealing and ligation). The LDR product (5 μL) was then added to 60 μL of hybridization solution [3 M tetramethylammonium chloride (TMAC), 50 mMTris–HCl (pH 8.0), 3 mM EDTA (pH 8.0), 0.10% SDS containing Luminex FlexMAP microspheres from each single nucleotide polymorphism (SNP)-specific set (total number of SNP-specific microspheres, n = 8)]. Mixtures were heated to 95 °C for 90 s and incubated at 37 °C for 40 min to allow hybridization between SNP-specific LDR products and microsphere-specific anti-TAG oligonucleotides.
Following hybridization, 6 μL of streptavidin-R-phycoerythrin (Molecular Probes/ThermoFisher) in TMAC hybridization solution (20 ng/μL) was added to the post-LDR mixture and incubated at 37 °C for 40 min in Costar 6511 M polycarbonate 96-well V-bottom plates (Corning). Hybrid complexes consisting of SNP-specific LDR products and microsphere-labelled anti-TAG probes were detected using a Bio-Plex array reader (Bio-Rad Laboratories); the plate temperature was set to 37 °C throughout detection. All fluorescence data were collected using Bio-Rad software, BioPlex Manager 5.0.
G6PDd phenotype prevalence
Recruitment of unrelated males in the Ampasimpotsy community (Tsiroanomandidy district) and Mandoto/Miandrivazo districts of Madagascar identified a total of 365 study participants. Mean age (SD) was 23 years (18.7) and 15 years (15.5) in each site, respectively. All were screened for their G6PD phenotype using the FST qualitative test. Overall, 47 study participants (12.9%) were diagnosed as being G6PD deficient. In Mandoto and Miandrivazo districts, 28 of the 230 individuals tested (12.1%) were deficient, with 61% (n = 17) of the deficient having a severe deficiency, and 39% (n = 11) intermediate deficiency. In Ampasimpotsy, 19 of the 135 study participants were G6PD deficient, including 63% (n = 12) severely and 37% (n = 7) intermediate deficient. Overall G6PDd prevalence was not significantly different between the two study populations surveyed (Chi squared test, p = 0.7178).
G6PD gene sequencing and genotyping
A subset of 95 samples was selected for full G6PD gene sequencing following PCR amplification of exons 2–13 and Illumina short-read sequencing. These samples included 15 of the 29 severely deficient, 14 of the 18 intermediate and 66 of the 318 G6PD normal individuals, as diagnosed by the FST. Overall, 16 SNPs were identified across the G6PD genomic DNA sequence examined (Table 2). This included six exon-based (three non-synonymous and three synonymous) and ten non-coding SNPs (seven intron and three 3′ non-coding); with the exception of the SNP identified at nucleotide X:154,546,029 (intron 2), all of the SNPs identified in this study had RefSNP accession ID numbers. The six coding region SNPs included nucleotide 202G > A (amino acid 68V > M), 376A > G (126N > D), 563C > T (188S > F), 1116G > A (372Q > Q), 1311C > T (437Y > Y), and 1431C > T (477P > P). The ten non-coding SNPs were identified in introns 2, 5, 7, 8, 11 and in the 3′ untranslated region (Table 2). These SNPs were observed in 20 different haplotypes resulting from Illumina sequencing. The haplotypes included the wild-type G6PD B allele and 13 additional variants predicted to share the G6PD normal enzyme activity phenotype. Haplotypes also included phenotypically deficient alleles G6PD A-202A/376G (two haplotypes), G6PD A 376G (three haplotypes) and G6PD Mediterranean 563T (one haplotype) alleles (haplotype details presented in Table 3). Based on these findings an LDR-FMA genotyping assay was developed to evaluate the presence of the identified coding region SNPs for all study samples for which there was sufficient blood available for DNA extraction.
Ligation detection reaction-FMA genotyping was applied to 208 samples, including 121 from Mandoto/Miandrivazo and 87 from Ampasimpotsy; 11 samples were analysed by both Illumina sequencing and LDR-FMA (all 11 samples demonstrated 100% concordance between methods). When combined with Illumina sequencing results, 291 individuals were surveyed for G6PD DNA sequence polymorphism (Fig. 2). As all individuals studied were males, all results correspond to single haplotypes given male X chromosome hemizygosity. Overall results revealed seven haplotypes characterized by exon-based nucleotide sequence changes across the study sites (Table 3). The majority of the haplotypes were distributed similarly across the study locations, although there was unique distribution of some lower frequency alleles. Of note, the 202A SNP was associated with the 376G among individuals across the Mandoto/Miandrivazo study sites while in Ampasimpotsy, 202A was independent of 376G in 5 of 11 individuals (45.5%). Of further interest, we observed no G6PD alleles characteristic of Asian origins or alternative G6PD A- mutations observed previously in other African populations (e.g. G6PD A-680T/376G, G6PD Santamaria 542T/376G, G6PD Betica-Selma 968C/376 ).
G6PD phenotype-genotype diagnostic concordance
Four of the observed haplotypes (G6PD B, B 1311T, A 376G and A 376G/1311T) were phenotypically normal (95%), while three (G6PD A-202A/376G, A-202A and Mediterranean 563T) were exclusively associated with a deficient phenotype (Table 3). When G6PD haplotypes were compared with FST phenotypes, results were 95.5% concordant. Overall sensitivity of the FST results was 0.951 with one G6PD B individual diagnosed as phenotypically deficient and 12 individuals who were genetically normal but characterized with intermediate enzyme activity. The genetically normal but phenotypically deficient individual was diagnosed positive for malaria by RDT at the time of G6PD screening (June 2015). Blinded retesting two months later (August 2015) resulted in a normal G6PD phenotype diagnosis. Specificity of the FST was 1.00 as all 24 individuals with G6PDd alleles were diagnosed with intermediate (n = 4) or deficient (n = 20) enzyme activity.
Across both sites, malaria RDT results were available from 360 participants (Fig. 2). The overall rate of infection was 2.5%, but significantly different across both sites (Chi squared test, p = 0.003), at 0.4% (n = 1/230) in the Mandoto/Miandrivazo population in July 2015 and 6.2% (n = 8/130) from the Ampasimpotsy sample in June and August 2015. There was no significant association between G6PD phenotype status and malaria infection in the sample (Chi squared test, p > 0.5).
The Malagasy population’s origins are unique, though the details relating to the island’s human colonization continue to intrigue and elude archaeologists today [29, 30]. Diverse sources of evidence clearly point to Austronesian origins tracing back to transoceanic trading routes from modern-day Indonesia (specifically, Borneo ) and admixture with Bantu communities from continental Africa . For example, the Malagasy language is rooted in the Austronesian languages  but also shows the major influence of Bantu dialects [34, 35]. Radiocarbon dating of ancient crop traces suggests that Malagasy crops from the 7th to 12th centuries were of Asian origin ; there is also, however, evidence of close ties with the Swahili trading system in western coastal areas [30, 37]. Genome studies are beginning to explore the resulting spatio-genetic structure of the contemporary population . At a coarse level, observations from mitochondrial and Y chromosome studies among Malagasy sub-populations have indicated genetic discontinuity between highland and coastal populations, with highland populations having mainly Indonesian origins and coastal communities being more African [15, 32]. The communities surveyed in the present study were in the highland fringe region (Fig. 1) and self-identified as being predominantly of the Merina highland ethnic group.
The areas selected for this G6PDd screening study were endemic with P. vivax transmission. Although the 14 day primaquine regimen is recommended for the radical cure of P. vivax by Madagascar’s National Strategic Plan (2013–2017) , as in most malaria-endemic regions , primaquine is not prescribed. Its numerous limitations as a drug  contribute to this, but the potential risk of inducing haemolytic anaemia in patients with a deficiency in glucose-6-phosphate dehydrogenase (G6PDd) enzyme activity levels is the primary limitation. Fears of causing harm by triggering haemolysis, together with logical difficulties of procuring primaquine tables of appropriate posology, have also prevented the application of single-dose treatment to block gametocyte transmission of P. falciparum, even in the context of mortal epidemics (Roll Back Malaria committee in Madagascar, November 2016, Pers Comm.). A greater understanding of the relative risks associated with the two applications of primaquine could impact on the burden of malaria in Madagascar. Madagascar is a rare country in 2016 for having substantially increased (>20%) malaria-attributable morbidity and mortality between 2010 and 2015 . All efforts must be made to ensure the toolbox of available interventions is not constrained by poor information. This study, therefore, addresses a key barrier to P. vivax radical cure through a study of the epidemiology of G6PDd in a P. vivax-endemic zone.
An overall G6PDd phenotype prevalence among males of 12.9% was diagnosed using the fluorescence spot test (12.2% in Mandoto/Miandrivazo and 14.1% in Ampasimpotsy). Global compilations of community-based G6PDd prevalence surveys and of genetic variants have been mapped and published [16, 18]. The methodology followed in the present study was consistent with that of these previous efforts, enabling direct comparisons. The observed G6PDd prevalence in males of 12.9% falls within the 50% confidence interval (11.5–30.3%) of previous predictions based on geostatistical analysis of available population surveys in 2012 . Madagascar was notable in that analysis for being one of the countries with the highest uncertainty surrounding its predicted G6PDd prevalence, emphasizing the need for additional surveys to inform the sparse evidence from the Malagasy population. Setting this result in its geographic context, with the aim of identifying neighbouring areas with potentially similar G6PDd epidemiology, could help to guide programmatic roll-out of single-dose primaquine in Madagascar, a country where no evaluation of this intervention have yet been carried out . The observed G6PDd prevalence in Madagascar was generally consistent with modelled estimates across continental Africa where 15 countries had estimated population-weighted G6PDd allele frequencies ≥15%, and prevalence in 19 countries was <10% . Observed prevalence in Madagascar was considerably lower than that estimated from Mozambique, Madagascar’s closest geographic neighbour in continental Africa, which had a modelled predicted prevalence of 21.1% (50% CI 14.7–29.8%) . To the east across the Indian Ocean, Indonesia had a predicted national prevalence of 7.1% (50% CI 5.3–9.4%).
Genetic diversity across the G6PDd Malagasy study population was low, dominated by variants commonly identified among African populations ; this was unexpectedly homogenous given the population’s origins. The limited sample size, which was powered to assess phenotypic prevalence, restricted the scope for comprehensive characterization of variant diversity in these populations. However, the lack of an Indonesian footprint on the observed sample was notable. G6PDd gene diversity in Indonesian populations is among the greatest globally, including many endemic variants ; these did not appear in the Malagasy population screened by the present study. An Indonesian genetic footprint would only be expected to be identifiable if the genetic trait was present in the founding population. Current evidence from genetic studies suggests that the main Indonesian arrival took place in the eighth century CE. This potentially falls after the emergence of some of the most prevalent mutations, estimated to be 1000-6357 years ago for G6PD A- [40, 41], 3330 years ago for the G6PD Mediterranean 563T variant , and 1575 years ago for the G6PD Mahidol 487A variant common across Asia . Though only a small subset of variants, and not including the main Indonesian variants, these dates would suggest that an ancestral Indonesian G6PDd footprint could be expected in the contemporary Malagasy population. A larger number of surveys are required to draw any firm conclusions about an Indonesian contribution to the G6PD gene pool in this region of Madagascar.
The present findings carry implications for public health policy. A prevalence of 12.9% G6PDd primarily caused by G6PD A-202A/376G is similar to that of many sub-Saharan African countries where P. falciparum transmission blocking is the main application of primaquine. The WHO has endorsed the use of the single low-dose regimen as a gametocytocide (0.25 mg/kg), a dose considered safe for patients with G6PDd (with the exception of pregnant women, infants aged <6 months and women breastfeeding infants aged <6 months) [43–46]. Certain areas of Madagascar, including the central highlands (from where the Ampasimpotsy population sample originate) and the arid south of the island are prone to local epidemics . In the event of rapidly increasing transmission, the single-dose primaquine regimen could be used in these areas to limit transmission and contribute to curtailing epidemic transmission.
The principal risks with primaquine come from the substantially higher dosing (0.5 mg/kg daily for 14 days) required for the radical cure of P. vivax. Given the observed presence of G6PDd in this region, appropriate risk assessments and diagnostic capacity need to be implemented to safely avert morbidity from cumulative relapses of P. vivax clinical episodes . While most G6PDd individuals tested for this study carried variants of relatively moderate severity in whom drug-induced haemolysis may be self-limiting , there was also a severe mutation observed, G6PD Mediterranean 563T, for whom primaquine-induced haemolysis would be much more dangerous  and likely require transfusion to ensure recovery. The WHO guidelines for primaquine and G6PDd screening are deliberately flexible, allowing for differing G6PDd epidemiology and public health infrastructure scenarios. Countries must determine their own risk assessments and reviews of existing experience with primaquine therapy. In the present context of Madagascar, the relatively poorly developed infrastructure for emergency medicine would advise caution regarding the blind use of primaquine. Similarly, given the observation of relatively common (>10%) G6PDd frequency, including occurrence of the severe variant G6PD Mediterranean 563T, it would be strongly advisable to screen all P. vivax patients for G6PD enzyme activity prior to administering radical cure doses of primaquine.
The present study has limitations. Anaemia is a known confounder of G6PD screening , but was not adjusted for in this study. The high concordance between phenotype and genotype results, however, suggests that this was not a major problem. The limited geographic range of the study, which was restricted to a single malaria transmission ecozone, restricts extrapolation of results across the country; more extensive sampling across different regions would help to gain perspectives on the national-level G6PDd epidemiology as well as the population’s genetic relatedness across the island. The area studied was nevertheless the country’s principal focus of P. vivax transmission and hence concerned the population most frequently implicated in potential primaquine drug administration. The present study was limited to male participants. This ensured clearly binary G6PD diagnoses and avoided the spectrum of phenotypes associated with heterozygous G6PDd allele carriage. While predictions based on theoretical genetic inheritance principles can estimate frequencies of G6PDd phenotypes and genotypes among females at the population level, large-scale public health programmes increasing access to P. vivax radical cure will require the introduction of a quantitative enzyme activity assay that can detect pre-determined acceptable thresholds for safely administering the drug to females . To date, no well-adapted point of care tests suited to this application exist. The 14 days of treatment with primaquine affords opportunities for discontinuing treatment in response to side-effects. In contrast, the drug tafenoquine, which is currently in Phase 3 clinical trials as a single-dose therapy, carries exacerbated risks as a haemolytic agent due to its longer half-life, and will therefore require a G6PDd screening diagnostic with appropriate sensitivity.
The study presented here explores the epidemiology of G6PDd among previously unstudied populations in P. vivax-endemic regions of Madagascar. Prevalence of deficiency being greater than 10% in males carries implications for the administration of primaquine therapy, emphasizing the need for appropriate G6PDd screening ahead of radical cure therapy. The hypothesized G6PD genetic diversity anticipated from the diverse origins of the Malagasy population was not reflected in the G6PD gene variants despite evidence from other genes, such as the Duffy blood group, suggesting mixed inheritance [6, 8]. More widespread investigations into the epidemiology of G6PDd would provide insight into malaria treatment risk assessments, as well as contribute evidence to ongoing investigations of human archaeology in Madagascar.
WHO. World Malaria Report. Geneva. Switzerland: World Health Organization; 2016. p. 2016.
Howes RE, Mioramalala SA, Ramiranirina B, Franchard T, Rakotorahalahy AJ, Bisanzio D, et al. Contemporary epidemiological overview of malaria in Madagascar: operational utility of reported routine case data for malaria control planning. Malar J. 2016;15:502.
Service des Statistiques Sanitaires-Ministère de la Santé Publique. Annuaire des Statistiques du Secteur Santé de Madagascar. Antananarivo: Ministry of Health; 2012.
Kesteman T, Randrianarivelojosia M, Mattern C, Raboanary E, Pourette D, Girond F, et al. Nationwide evaluation of malaria infections, morbidity, mortality, and coverage of malaria control interventions in Madagascar. Malar J. 2014;13:465.
Institut National de la Statistique (INSTAT), Programme National de Lutte contre le Paludisme (PNLP), Institut Pasteur de Madagascar (IPM), and ICF International. Madagascar Malaria Indicator Survey 2013 [Enquête sur les Indicateurs du Paludisme (EIPM)]. Calverton: INSTAT, PNLP, IPM and ICF International; 2013.
Menard D, Barnadas C, Bouchier C, Henry-Halldin C, Gray LR, Ratsimbasoa A, et al. Plasmodium vivax clinical malaria is commonly observed in Duffy-negative Malagasy people. Proc Natl Acad Sci USA. 2010;107:5967–71.
Barnadas C, Ratsimbasoa A, Tichit M, Bouchier C, Jahevitra M, Picot S, et al. Plasmodium vivax resistance to chloroquine in Madagascar: clinical efficacy and polymorphisms in pvmdr1 and pvcrt-o genes. Antimicrob Agents Chemother. 2008;52:4233–40.
Howes RE, Patil AP, Piel FB, Nyangiri OA, Kabaria CW, Gething PW, et al. The global distribution of the Duffy blood group. Nat Commun. 2011;2:266.
Zimmerman PA, Ferreira MU, Howes RE, Mercereau-Puijalon O. Red blood cell polymorphism and susceptibility to Plasmodium vivax. Adv Parasitol. 2013;81:27–76.
Howes RE, Franchard T, Rakotomanga TA, Rosenjack J, Zikursh M, Ramanatiaray T, et al. Prevalence of sub-microscopic infections: the role for molecular diagnostics in Madagascar. In: Annual meeting of the American society of tropical medicine and hygiene. Philadelphia: Poster P1511; 2015.
National malaria control programme of Madagascar. National strategic plan for malaria control in Madagascar 2013–2017: Consolidating the gains with a view to elimination of malaria from Madagascar, 2015–2017 revision. 2015.
Chen I, Poirot E, Newman M, Kandula D, Shah R, Hwang J, et al. An assessment of the supply, programmatic use, and regulatory issues of single low-dose primaquine as a Plasmodium falciparum gametocytocide for sub-Saharan Africa. Malar J. 2015;14:204.
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.
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.
Tofanelli S, Bertoncini S, Castri L, Luiselli D, Calafell F, Donati G, et al. On the origins and admixture of Malagasy: new evidence from high-resolution analyses of paternal and maternal lineages. Mol Biol Evol. 2009;26:2109–24.
Howes RE, Piel FB, Patil AP, Nyangiri OA, Gething PW, Hogg MM, et al. G6PD deficiency prevalence and estimates of affected populations in malaria endemic countries: a geostatistical model-based map. PLoS Med. 2012;9:e1001339.
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;9:e0003602.
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.
Mayer G, Mayoux A. Recherche d’anomalies sanguines génétiques dans une population de la côte est de Madagascar. Arch Inst Pasteur Madag. 1966;35:209–11.
Cappellini MD, Fiorelli G. Glucose-6-phosphate dehydrogenase deficiency. Lancet. 2008;371:64–74.
OpenEpi: Open source epidemiologic statistics for public health, Version 3 [www.OpenEpi.com].
Bancone G, Chu CS, Chowwiwat N, Somsakchaicharoen R, Wilaisrak P, Charunwatthana P, et al. Suitability of capillary blood for quantitative assessment of G6PD activity and performances of G6PD point-of-care tests. Am J Trop Med Hyg. 2015;92:818–24.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.
Benish RL, Rodriguez B, Zimmerman PA, Mehlotra RK. Comparative description of haplotype structure and genetic diversity of MDR1 (ABCB1) in HIV-positive and HIV-negative populations. Infect Genet Evol. 2010;10:60–7.
Mehlotra RK, Cheruvu VK, Blood Zikursh MJ, Benish RL, Lederman MM, Salata RA, et al. Chemokine (C-C motif) receptor 5–2459 genotype in patients receiving highly active antiretroviral therapy: race-specific influence on virologic success. J Infect Dis. 2011;204:291–8.
Mehlotra RK, Bockarie MJ, Zimmerman PA. CYP2B6 983T > C polymorphism is prevalent in West Africa but absent in Papua New Guinea: implications for HIV/AIDS treatment. Br J Clin Pharmacol. 2007;64:391–5.
Clark TG, Fry AE, Auburn S, Campino S, Diakite M, Green A, et al. Allelic heterogeneity of G6PD deficiency in West Africa and severe malaria susceptibility. Eur J Hum Genet. 2009;17:1080–5.
Radimilahy CM, Crossland Z. Situating Madagascar: Indian Ocean dynamics and archaeological histories. Azania Archaeol Res Afr. 2015;50:495–518.
Campbell G. Early exchange between Africa and the Wider Indian Ocean World. London: Palgrave Macmillan; 2016.
Brucato N, Kusuma P, Cox MP, Pierron D, Purnomo GA, Adelaar A, et al. Malagasy genetic ancestry comes from an historical malay trading post in southeast Borneo. Mol Biol Evol. 2016;33:2396–400.
Hodgson JA. A genomic investigation of the malagasy confirms the highland-coastal divide, and the lack of Middle Eastern gene flow. In: Campbell GH, editor. Early exchange between Africa and the Wider Indian Ocean World. London: Palgrave Macmillan; 2016.
Kusuma P, Brucato N, Cox MP, Pierron D, Razafindrazaka H, Adelaar A, et al. Contrasting linguistic and genetic origins of the asian source populations of malagasy. Sci Rep. 2016;6:26066.
Serva M, Petroni F, Volchenkov D, Wichmann S. Malagasy dialects and the peopling of Madagascar. J R Soc Interface. 2012;9:54–67.
Blench R, Walsh M. Faunal names in Malagasy: their etymologies and implications for the prehistory of the East African coast. In: Proc 11th Int Conf on Autronesian Linguistics (11 ICAL). France: Aussois; 2009.
Crowther A, Lucas L, Helm R, Horton M, Shipton C, Wright HT, et al. Ancient crops provide first archaeological signature of the westward Austronesian expansion. Proc Natl Acad Sci USA. 2016;113:6635–40.
Radimilahy C. Mahilaka: an archaeological investigation of an early town in northwestern Madagascar. Uppsala: Societas Archaeologica Uppsaliensis; 1998.
National Malaria Control Programme of Madagascar. National strategic plan for malaria control in Madagascar 2013–2017. Reinforcing the gains of malaria control with a perspective towards malaria elimination in Madagascar. 2012.
Baird JK. Primaquine toxicity forestalls effective therapeutic management of the endemic malarias. Int J Parasitol. 2012;42:1049–54.
Tishkoff SA, Varkonyi R, Cahinhinan N, Abbes S, Argyropoulos G, Destro-Bisol G, et al. Haplotype diversity and linkage disequilibrium at human G6PD: recent origin of alleles that confer malarial resistance. Science. 2001;293:455–62.
Sabeti PC, Reich DE, Higgins JM, Levine HZ, Richter DJ, Schaffner SF, et al. Detecting recent positive selection in the human genome from haplotype structure. Nature. 2002;419:832–7.
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.
WHO. Guidelines for the treatment of malaria. 3rd ed. Geneva: World Health Organization; 2015.
Malaria policy advisory committee meeting (MPAC). WHO. Evidence Review Group Report: the safety and effectiveness of single dose primaquine as a P. falciparum gametocytocide. Geneva: World Health Organization; 2012.
Bancone G, Chowwiwat N, Somsakchaicharoen R, Poodpanya L, Moo PK, Gornsawun G, et al. Single low dose primaquine (0.25 mg/kg) does not cause clinically significant haemolysis in G6PD deficient subjects. PLoS ONE. 2016;11:e0151898.
Eziefula AC, Bousema T, Yeung S, Kamya M, Owaraganise A, Gabagaya G, et al. Single dose primaquine for clearance of Plasmodium falciparum gametocytes in children with uncomplicated malaria in Uganda: a randomised, controlled, double-blind, dose-ranging trial. Lancet Infect Dis. 2014;14:130–9.
Douglas NM, Anstey NM, Buffet PA, Poespoprodjo JR, Yeo TW, White NJ, et al. The anaemia of Plasmodium vivax malaria. Malar J. 2012;11:135.
Dern RJ, Beutler E, Alving AS. The hemolytic effect of primaquine. II. The natural course of the hemolytic anemia and the mechanism of its self-limited character. J Lab Clin Med. 1954;44:171–6.
Piomelli S, Corash LM, Davenport DD, Miraglia J, Amorosi EL. In vivo lability of glucose-6-phosphate dehydrogenase in GdA− and GdMediterranean deficiency. J Clin Invest. 1968;47:940–8.
Ley B, Luter N, Espino FE, Devine A, Kalnoky M, Lubell Y, et al. The challenges of introducing routine G6PD testing into radical cure: a workshop report. Malar J. 2015;14:377.
REH, AR and PAZ conceived the study; REH wrote the first draft of the paper with PAZ. REH and AR oversaw field activities. TAR, TF and BR ran field collections. SS and MZ conducted all laboratory assays. REH, JG, ERC and PAZ analysed the results. All authors read and approved the final manuscript.
The authors thank all study participants, local health officials and field doctors for their participation and support. Study permission from Association ASA (Ankohonana Sahirana Arenina, www.asa-madagascar.org) and logistical support Pact Madagascar (www.pact-madagascar.org) is also acknowledged.
The authors declare that they have no competing interests.
Availability of data and materials
Data are available on request from the authors.
Consent for publication
All authors read and approved the final manuscript’s submission for publication.
Ethics approval and consent to participate
This study’s protocol was approved by the ethical review panels of University Hospitals Case Medical Centre, Cleveland, Ohio, USA (DMID Protocol #13-0067); the National Institutes of Health, USA (DMID Protocol 13-0067); and the Ministry of Public Health, Madagascar (Nº099-MSANP/CE). Written informed consent was obtained from all subjects, or subject guardians, enrolled into the study.
This study was funded by an NIH Grant (R01 AI097366) to PAZ and a Core Utilization Pilot Grant from the Case Western Reserve University Clinical and Translational Science Collaborative Fund (UL1TR000439) to ERC. Additional support to ERC comes from the Computational Genomic Epidemiology of Cancer Fellowship (CoGEC, 14R25CA094186).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Howes, R.E., Chan, E.R., Rakotomanga, T.A. et al. Prevalence and genetic variants of G6PD deficiency among two Malagasy populations living in Plasmodium vivax-endemic areas. Malar J 16, 139 (2017). https://doi.org/10.1186/s12936-017-1771-6
- Glucose-6-phosphate dehydrogenase deficiency
- G6PDd genotypes
- Plasmodium vivax