This article has Open Peer Review reports available.
Genetic determinants of glucose-6-phosphate dehydrogenase activity in Kenya
© Shah et al.; licensee BioMed Central Ltd. 2014
Received: 31 March 2014
Accepted: 31 July 2014
Published: 9 September 2014
The relationship between glucose-6-phosphate dehydrogenase (G6PD) deficiency and clinical phenomena such as primaquine-sensitivity and protection from severe malaria remains poorly defined, with past association studies yielding inconsistent and conflicting results. One possibility is that examination of a single genetic variant might underestimate the presence of true effects in the presence of unrecognized functional allelic diversity.
We systematically examined this possibility in Kenya, conducting a fine-mapping association study of erythrocyte G6PD activity in 1828 Kenyan children across 30 polymorphisms at or around the G6PD locus.
We demonstrate a strong functional role for c.202G>A (rs1050828), which accounts for the majority of variance in enzyme activity observed (P=1.5 × 10-200, additive model). Additionally, we identify other common variants that exert smaller, intercorrelated effects independent of c.202G>A, and haplotype analyses suggest that each variant tags one of two haplotype motifs that are opposite in sequence identity and effect direction. We posit that these effects are of biological and possible clinical significance, specifically noting that c.376A>G (rs1050829) augments 202AG heterozygote risk for deficiency trait by two-fold (OR = 2.11 [1.12 - 3.84], P=0.014).
Our results suggest that c.202G>A is responsible for the majority of the observed prevalence of G6PD deficiency trait in Kenya, but also identify a novel role for c.376A>G as a genetic modifier which marks a common haplotype that augments the risk conferred to 202AG heterozygotes, suggesting that variation at both loci merits consideration in genetic association studies probing G6PD deficiency-associated clinical phenotypes.
Glucose-6-phosphate dehydrogenase (G6PD) is a housekeeping enzyme encoded on chromosome X. It is ubiquitously expressed in all body tissues and found in eukaryotes and prokaryotes alike. G6PD catalyses the rate-limiting first step in the pentose phosphate pathway, generating NADPH, an electron donor important in a variety of metabolic processes, including fatty and nucleic acid synthesis, nitric oxide production, and the regeneration of antioxidant glutathione. Deficiency of G6PD is the most common enzymopathy in man, with prevalence especially high in sub-Saharan Africa and southeast Asia. It is a trait that is thought to be both harmful and beneficial. While clinical penetrance is generally mild, G6PD deficiency can predispose to hemolytic anemia and neonatal jaundice, generally in the context of oxidant stressors generated by certain foods, medications, and infections. Counterbalancing these deleterious effects is evidence that suggests G6PD deficiency confers partial protection against malaria, and that such protection has conferred a selective advantage to G6PD-deficient individuals in malaria-endemic regions, as reflected in the close geographical correlation between the global distribution of G6PD deficiency and that of historical malaria endemicity.
Despite decades of research, however, our understanding of both beneficial and deleterious clinical phenotypes associated with G6PD deficiency remains incomplete. Association studies of severe malaria have disagreed over the strength and specificity of protection, especially with respect to gender, with past studies demonstrating protection for only males , only females , both ,, or neither . Additionally, there is geographical heterogeneity in association testing results as well, with the only compelling evidence thus far coming from west Africa. Similarly, the precise determinants of the hemolytic response to oxidative stressors remain unknown. One especially important drug reaction is that involving the antimalarial primaquine ,. It is currently the only medication capable of radical cure of Plasmodium vivax infection and one of few drugs that acts against the Plasmodium falciparum gametocyte; it therefore has high value in eliminating reservoirs of infection and reducing transmission rates, tasks that are crucial to the success of any malaria elimination program . While G6PD deficiency alleles are strongly associated with increased risk to primaquine-induced hemolysis, there exists significant unexplained variance in this association, especially for mosaic female heterozygotes .
One challenging aspect to designing and interpreting clinical association studies lies in how G6PD deficiency is defined. While many genetic association studies at G6PD implicitly assume that one nonsynonymous coding variant has both full positive and negative predictive power for G6PD deficiency trait, this approach belies the considerable allelic heterogeneity present at the locus ,,. Indeed, G6PD deficiency is perhaps best understood as a family of heritable traits exhibiting widely varying prevalence, as well as differential biochemical and clinical penetrance, all factors which influence the power of an association test . Additionally, it is possible that unrecognized genetic modifiers of G6PD activity have yet to be identified, and that these may play a role in determining an individual’s G6PD deficiency status.
To inform our understanding of how genetic variation at G6PD impacts clinical outcomes, we sought to examine its more immediate relationship with molecular phenotypes. We took a comprehensive two-part approach, first resequencing the G6PD gene to validate known and discover novel variants, and then testing those polymorphisms for association with G6PD enzyme activity in Kenya.
Ethics approval was granted by the Kenya Medical Research Institute (KEMRI) review board, and written informed consent was obtained from a parent or guardian prior to enrollment.
Capillary blood was drawn by heel prick from 2100 healthy infants (less than 6 months old) of predominantly Mijikenda ethnicity, recruited into an ongoing prospective birth cohort study run by KEMRI and the Wellcome Trust in Kilifi, Kenya (SCC1192: The genomic epidemiology of childhood diseases in Kilifi, Kenya). Samples were stored in EDTA at 4°C and processed for activity assays and DNA extraction within four days. Of note, serial assay measurements taken on the same samples remained consistent even up to 10-14 days from collection.
Enzyme activity assay
Sequencing and genotyping
Resequencing was performed as part of an exon sequencing project called ExoSeq at Wellcome Trust Sanger Institute (WTSI). A total of 288 individuals comprised our sample set, with 48 individuals each from Burkina Faso, Cameroon, Nigeria, Kenya, Vietnam, and Papua New Guinea. A total of 28 amplicons approximately 600-800 bp in length were designed via automated algorithm (i.e. in silico PCR), covering about two-thirds of the 15 kb G6PD gene, with full coverage of all exons, but a few gaps in noncoding sequence, specifically in introns 1 and 2. Capillary sequencing of the PCR products generated yielded chromatograms that were input into an automated SNP calling pipeline, which was supplemented with manual calls from examination of individual sequence traces. A total of 72 SNPs, including forty-nine novel variants were identified.
Genotyping was performed on the Sequenom®; MASSarray®; iPLEX®; platform. A total of three multiplex genotyping reactions were designed, assaying a total of 68 SNPs. The SNPs input to the design software included those discovered during resequencing (N = 72), as well as nonsynonymous variants from the literature (N = 25), and several hundred variants within 1 Mb of G6PD extracted from dbSNP and Ensembl databases. Inclusion in the multiplex reactions was priority-weighted according to predicted functional consequence and estimated minor allele frequency criteria. Template DNA used for genotyping was genomic DNA amplified by primer extension pre-amplification . Genotype calls were manually curated by examining cluster plots to identify ambiguous genotypes. Ancestral and derived alleles were inferred according to the Ensembl EPO pipeline.
Data curation and analysis
Enzyme activity data was filtered to keep only individuals with positive absorbance readings at both 340 and 540 nm, and checked for batch effects and extreme outliers. Genotype data was filtered using custom Perl and R scripts, along with PLINK  and snpMatrix . Individuals exhibiting greater than 10% missing data or ambiguous gender (based on chromosome X genotypes) were excluded. SNPs were removed if they showed greater than 5% missing data, minor allele frequency < 1%, or extreme deviation from Hardy-Weinberg equilibrium (P<10-4, Additional file 1). In total, 1828 of 2100 individuals and 30 out of 70 SNPs met these inclusion criteria (Additional file 2).
Statistical analyses were conducted in R . Single marker tests of association for each derived allele were conducted via multiple linear regression. Covariates included sex, and in the controlled tests, c.202G>A genotype. Three different modes of inheritance are considered- in the additive model, male genotypes were coded 0 or 2, and females were 0, 1, or 2; in the dominant and recessive models, males were coded 0/1, and females were coded 0/1/1, or 0/0/1, respectively. Haplotypes were inferred using fastPHASE  under default settings with known male haplotypes used as a reference input to reduce search space and minimize phasing error rate in females (Additional file 3). All haplotypes with frequency > 1% (26 full-length, 13 partial) were tested for association; individuals were coded 0 or 1 based on absence/presence of the test haplotype, with females represented by one haplotype selected at random. Odds ratios (OR) for qualitative tests were calculated via logistic regression. Effect sizes (ES) are presented here with 95% confidence interval in brackets. Allele sharing distance was calculated using a Euclidean distance metric for each individual from both H+ and H- haplotypes. An individual is categorized as either H+ or H- if their distance metric falls into the lowest quintile of the respective empirical distribution. In linear modeling of the additive effect of allele sharing distance on enzyme activity, this categorization is coded as -1 (H-), 0 (neither), or +1 (H+).
The variants studied here were ascertained using a combination of sources, including our own resequencing efforts, published accounts in the literature, and online variation databases. Resequencing was carried out as part of an ongoing laboratory effort to catalogue global genetic variation at the G6PD locus (manuscript in preparation). A total of 15 kb at the G6PD locus was sequenced in a worldwide panel comprised of 288 African and Asian individuals, including 48 from Kenya. This work identified a total of 72 SNPs, including 49 novel variants. Published accounts in the literature yielded an additional 25 nonsporadic, nonsynonymous variants associated with G6PD deficiency in Africa and Asia, including 5 that have been identified in sub-Saharan Africa : the c.202G>A, c.376A>G, c.542A>T, c.680G>T, and c.968T>C polymorphisms (number refers to cDNA coordinate). Finally, several hundred polymorphisms (mostly noncoding SNPs) within 1 Mb of G6PD were also identified in dbSNP and Ensembl to supplement our list.
This list of several hundred SNPs was then priority-weighted according to predicted functional consequence (e.g. nonsynonymous, splice site, etc.), estimated allele frequency, and proximity to G6PD in order to decide which would be prioritized for inclusion in genotyping assays. A total of 68, including all five nonsynonymous African variants, were typed in the current study using the Sequenom®; MASSarray®; iPLEX®; platform. Of the SNPs typed, 30 met final study inclusion criteria, which included having derived allele frequency (DAF) greater than 0.01 and less than 5% missing data. It should be noted that most excluded SNPs failed based on DAF criteria, as many were private to west Africa or Asia. This included the c.542A>T, c.680G>T, and c.968T>C variants, which were monomorphic (DAF = 0%) in Kenya. The only two nonsynonymous variants at the locus that were found to be polymorphic in Kenya were c.202G>A and c.376A>G, exhibiting DAF of 17% and 40%, respectively. All SNPs included in the association study are shown in Additional file 2.
The c.202G>A polymorphism is the major determinant of G6PD activity in Kenya
Modifiers of G6PD activity that are independent of c.202G>A
Minimal adequate genetic model of G6PD activity
Allelic ES (95% CI)
-0.198 (-0.213 - -0.183)
+0.022 (+0.011 - +0.034)
-0.017 (-0.030 - -0.005)
While the association of c.376A>G with a modest reduction in enzyme activity (approximately 10-20% in male hemizygotes) has been reported previously ,, and this has been generally thought to be the result of changes to protein structure leading to enzyme instability, it should be emphasized that neither the controlled association scan nor subsequent haplotype testing was suggestive of a uniquely causal role for c.376A>G. Indeed, what our work suggests instead is the possibility that this nonsynonymous variant is merely a marker, as it is only one of nine SNPs tagging the H- haplotype and one of twelve overall that exhibits modest effects independent of c.202G>A.
Both c.202G>A and c.376A>G contribute to risk of G6PD deficiency state
While it is known that the distribution of enzyme activity in 202AG female heterozygotes is broad, with some individuals exhibiting an effectively deficient biochemical phenotype, little is known about the genetic factors that might contribute to this risk. We tested the possibility that our H- tag SNP, c.376A>G, might reduce enzyme activity sufficiently to augment risk of phenotypic G6PD deficiency state in such individuals. The 202A allele is known to obligately lie on a 376G background in African populations to form the `A-’ haplotype ,,, and thus 202AG individuals can be either 376GA or 376GG. We found that among 202AG heterozygotes, those that were 376GG exhibited not only reduced enzyme activity as compared to 376GA (ES = -0.06 [-0.10 - -0.03], P=0.0009, Figure 4b), a trend consistent with that observed earlier in 202 wild-type individuals, but were also notably at increased risk for deficiency state (OR = 2.11 [1.12 - 3.84], P=0.014). This latter finding has not been previously reported, and represents the first demonstration of a genetic modifier that influences risk to deficiency state in individuals heterozygous for the A- form of G6PD deficiency.
The work presented here represents the first fine-mapping association study of enzyme activity at the G6PD locus, and provides new insights into the genetic basis for quantitative and qualitative variation in this clinically-important molecular phenotype. We found that the largest proportion of variance in G6PD activity could be explained by a single polymorphism, c.202G>A, a coding variant thought to be a target of recent positive selection ,, but one that has yielded conflicting results in previous association studies of severe malaria. In addition, we found that certain common variants at the locus modulate G6PD activity independent of c.202G>A, and identify among these a role for c.376A>G as a genetic modifier which augments by two-fold the risk of G6PD deficiency state in 202AG heterozygote females. Taken together, these results suggest that future clinical association studies in Africa should consider the effects of both of these variants in order to arrive at a more accurate genetic definition of phenotypic G6PD deficiency state.
Our study is the first to systematically consider the possibility that other SNPs, even those found several hundred kilobases away from G6PD, might be associated with G6PD deficiency, and that c.202G>A might only be indirectly associated with deficiency, a finding which could explain discrepancies in previous association studies. Although the number of SNPs surveyed here is relatively small and does not include rare variants, we did not observe any SNPs more strongly associated with G6PD deficiency in Kenya than the c.202G>A variant, strongly suggestive of a causal role. Additionally, while all known functional variants in sub-Saharan Africa were genotyped as part of our study, only c.202G>A was found to be polymorphic in Kenya (DAF = 17%). This allelic homogeneity contrasts with the situation in some areas of west Africa , and suggests that given a similar allele frequency and biochemical penetrance, clinical association studies of malaria in each region should be equivalently powered if indeed the G6PD deficiency trait is directly responsible for protection. The strongest evidence implicating G6PD deficiency in malaria protection would ideally be consistent with an additive genetic model, with hemizygous males and homozygous females exhibiting similar clinical phenotypes, and heterozygous females displaying intermediate phenotypes.
Our work also demonstrates the power of a haplotype-based approach in resolving multiple correlated association signals in a region. Though we found two haplotype motifs that were common and of significant effect independent of c.202G>A, we were unable to further localize the signal, even when examining all possible two- and three-locus haplotypes (data not shown), suggesting that these haplotypes are merely markers for the true causal variants. Interestingly, among the twelve variants of significant effect independent of c.202G>A were c.376A>G, a nonsynonymous coding variant, as well as rs915942 and rs111827785, predicted 5’ UTR splice-site variants within RPL10 (a large-subunit ribosomal gene) and G6PD, respectively. While we can only speculate as to the precise nature of the structure-function relationship mediating these effects, the biology of the locus makes for interesting possibilities for interaction between c.202G>A and the H+/H- variants, including the possibility of cell-cell interactions synergistically impacting enzyme activity distribution in 202AG mosaic females.
At one level, the work presented here reveals how common, standing genetic variation superimposed on a Mendelian trait can significantly impact expression of that trait. More important, however, is how this finding might translate into a better understanding of how G6PD deficiency is linked to key clinical phenomena. Specifically, it will be worth considering whether two site genetic testing is more predictive of clinical outcome than single SNP testing in women heterozygous for the common A- form of G6PD deficiency in Africa. For example, in considering a public health intervention such as mass primaquine administration as part of a malaria elimination campaign, two-site testing might help to identify the subset of A- heterozygotes more likely to be susceptible to acute hemolytic anemia. Lastly, at a more general level, our work highlights the importance of measuring intermediate molecular phenotypes that can bridge the gulf between variance at the genetic level and variance in clinical phenotypes. This is especially important for association studies of G6PD deficiency given its considerable allelic heterogeneity, varying biochemical penetrance, and long history of inconsistent findings.
We have surveyed genetic variation at the G6PD locus in Kenya via resequencing, and conducted a fine-mapping association study of erythrocyte G6PD activity. We find that while the c.202G>A polymorphism accounts for most of the variance in enzyme activity seen in this population, several other common polymorphisms are also relevant, specifically noting that c.376A>G acts as a genetic modifier which increases by two-fold the risk of deficiency state in 202AG heterozygotes. Our results suggest that multiple alleles at the G6PD locus play a role in determining G6PD deficiency state and therefore merit consideration in genetic association studies of clinical phenotypes such as severe malaria or drug sensitivity.
SS, KR, TEW, DPK and TNW designed the study. SS and KR designed and validated the enzyme activity and genotyping assays. AM, JM and SU performed enzyme activity assays. SS, RC and CH performed genotyping. KK performed sequencing. SS conducted statistical analyses. SS, KR, TEW, DPK and TNW wrote the paper. All authors read and approved the final manuscript.
The authors thank all the children and their families who participated in this study. We also thank Metrine Tendwa, Adan Mohammed, and Emily Nyatichi for their efforts in support of this work. This paper is published with the permission of the Director of the Kenya Medical Research Institute. Work was supported in part by the Wellcome Trust (WT077383/Z/05/Z) and the Bill and Melinda Gates Foundation through the Foundations of the National Institutes of Health (566) as part of the Grand Challenges in Global Health Initiative, and the Intramural Research Program of the National Institute of Allergy and Infectious Diseases at the National Institutes of Health (NIH). The Wellcome Trust also provides core awards to Wellcome Trust Centre for Human Genetics (075491/Z/04; 090532/Z/09/Z) and the Wellcome Trust Sanger Institute (077012/Z/05/Z). TNW is funded by the Wellcome Trust (WT Senior Research Fellowship 091758) and the European Union EVIMALAR FP7 Network of Excellence. Travel funds for SS were furnished by Wolfson College and Oxford University Exploration Club. SS is also supported by the Medical Scientist Training Program at the NIH. DPK receives support from the Medical Research Council (G19/9).
- Guindo A, Fairhurst RM, Doumbo OK, Wellems TE, Diallo DA: X-linked G6PD deficiency protects hemizygous males but not heterozygous females against severe malaria. PLoS Med. 2007, 4 (3): 66-10.1371/journal.pmed.0040066.View ArticleGoogle Scholar
- Bienzle U, Ayeni O, Lucas A, 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 (7742): 107-10.1016/S0140-6736(72)90676-9.View ArticlePubMedGoogle Scholar
- Ruwende C, Khoo SC, Snow RW, Yates SN, Kwiatkowski D, Gupta S, Warn P, Allsopp CE, Gilbert SC, Peschu N: Natural selection of hemi- and heterozygotes for G6PD deficiency in Africa by resistance to severe malaria. Nature. 1995, 376 (6537): 246-249. 10.1038/376246a0.View ArticlePubMedGoogle Scholar
- Clark T, Fry A, Auburn S, Campino S, Diakite M, Green A, Richardson A, Teo Y, Small K, Wilson J, Jallow M, Sisay-Joof F, Pinder M, Sabeti P, Kwiatkowski D, Rockett K: Allelic heterogeneity of G6PD deficiency in West Africa and severe malaria susceptibility. Eur J Hum Genet. 2009, 17 (8): 1080-1085. 10.1038/ejhg.2009.8.View ArticlePubMedPubMed CentralGoogle Scholar
- Jallow M, Teo Y, Small K, Rockett K, Deloukas P, Clark T, Kivinen K, Bojang K, Conway D, Pinder M, Sirugo G, Sisay-Joof F, Usen S, Auburn S, Bumpstead S, Campino S, Coffey A, Dunham A, Fry A, Green A, Gwilliam R, Hunt S, Inouye M, Jeffreys A, Mendy A, Palotie A, Potter S, Ragoussis J, Rogers J, Rowlands K, et al: Genome-wide and fine-resolution association analysis of malaria in West Africa. Nat Genet. 2009, 41 (6): 657-665. 10.1038/ng.388.View ArticlePubMedPubMed CentralGoogle Scholar
- Dern RJ, Beutler E, Alving A: The hemolytic effect of primaquine. II. The natural course of the hemolytic anemia and the mechanism of its self-limited character. J Lab Clin Med. 1954, 44 (2): 171-176.PubMedGoogle Scholar
- Carson PE, Flanagan C, Ickes C, Alving A: Enzymatic deficiency in primaquine-sensitive erythrocytes. Science. 1956, 124 (3220): 484-485. 10.1126/science.124.3220.484-a.View ArticlePubMedGoogle Scholar
- Howes RE, Piel FB, Patil AP, Nyangiri OA, Gething PW, Dewi M, Hogg MM, Battle KE, Padilla CD, Baird JK, Hay SI: G6PD deficiency prevalence and estimates of affected populations in malaria endemic countries: a geostatistical model-based map. PLoS Med. 2012, 9 (11): 1001339-10.1371/journal.pmed.1001339.View ArticleGoogle Scholar
- Shekalaghe SA, ter Braak R, Daou M, Kavishe R, van den Bijllaardt W, van den Bosch S, Koenderink JB, Luty AJF, Whitty CJM, Drakeley C, Sauerwein RW, Bousema T: In Tanzania, hemolysis after a single dose of primaquine coadministered with an artemisinin is not restricted to glucose-6-phosphate dehydrogenase-deficient (G6PD A-) individuals. Antimicrob Agents Chemother. 2010, 54 (5): 1762-1768. 10.1128/AAC.01135-09.View ArticlePubMedPubMed CentralGoogle Scholar
- Beutler E: G6PD: population genetics and clinical manifestations. Blood Rev. 1996, 10 (1): 45-52. 10.1016/S0268-960X(96)90019-3.View ArticlePubMedGoogle Scholar
- Ganczakowski M, Town M, Bowden DK, Vulliamy TJ, Kaneko A, Clegg JB, Weatherall DJ, Luzzatto L: Multiple glucose 6-phosphate dehydrogenase-deficient variants correlate with malaria endemicity in the Vanuatu archipelago (southwestern Pacific). Am J Hum Genet. 1995, 56 (1): 294-301.PubMedPubMed CentralGoogle Scholar
- Zhang L, Cui X, Schmitt K, Hubert R, Navidi W, Arnheim N: Whole genome amplification from a single cell: implications for genetic analysis. Proc Natl Acad Sci U S A. 1992, 89 (13): 5847-5851. 10.1073/pnas.89.13.5847.View ArticlePubMedPubMed CentralGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, de Bakker PIW, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81 (3): 559-575. 10.1086/519795.View ArticlePubMedPubMed CentralGoogle Scholar
- Clayton D, Leung H-T: An R package for analysis of whole-genome association studies. Hum Hered. 2007, 64 (1): 45-51. 10.1159/000101422.View ArticlePubMedGoogle Scholar
- R Development Core Team (2008): R: A language and environment for statistical computing. R Foundation for Statistical Computing: Vienna, Austria;2008. ISBN 3-900051-07-0, ., [http://www.R-project.org]
- Scheet P, Stephens M: A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. Am J Hum Genet. 2006, 78 (4): 629-644. 10.1086/502802.View ArticlePubMedPubMed CentralGoogle Scholar
- Beutler E, Kuhl W, Vives-Corrons JL, Prchal JT: Molecular heterogeneity of glucose-6-phosphate dehydrogenase A-. Blood. 1989, 74 (7): 2550-2555.PubMedGoogle Scholar
- May J, Meyer CG, Grossterlinden L, Ademowo OG, Mockenhaupt FP, Olumese PE, Falusi AG, Luzzatto L, Bienzle U: Red cell glucose-6-phosphate dehydrogenase status and pyruvate kinase activity in a Nigerian population. Trop Med Int Health. 2000, 5 (2): 119-123. 10.1046/j.1365-3156.2000.00529.x.View ArticlePubMedGoogle Scholar
- Johnson MK, Clark TD, Njama-Meya D, Rosenthal PJ, Parikh S: Impact of the method of G6PD deficiency assessment on genetic association studies of malaria susceptibility. PLoS ONE. 2009, 4 (9): 7246-10.1371/journal.pone.0007246.View ArticleGoogle Scholar
- Nance WE: Quantitative studies of glucose-6-phosphate dehydrogenase. Am J Hum Genet. 1977, 29 (5): 537-544.PubMedPubMed CentralGoogle Scholar
- Battistuzzi G, Esan GJ, Fasuan FA, Modiano G, Luzzatto L: Comparison of GdA and GdB activities in Nigerians. A study of the variation of the G6PD activity. Am J Hum Genet. 1977, 29 (1): 31-36.PubMedPubMed CentralGoogle Scholar
- Hirono A, Beutler E: Molecular cloning and nucleotide sequence of cDNA for human glucose-6-phosphate dehydrogenase variant A(-). Proc Natl Acad Sci U S A. 1988, 85 (11): 3951-3954. 10.1073/pnas.85.11.3951.View ArticlePubMedPubMed CentralGoogle Scholar
- Vulliamy TJ, Othman A, Town M, Nathwani A, Falusi AG, Mason PJ, Luzzatto L: Polymorphic sites in the African population detected by sequence analysis of the glucose-6-phosphate dehydrogenase gene outline the evolution of the variants A and A-. Proc Natl Acad Sci U S A. 1991, 88 (19): 8568-8571. 10.1073/pnas.88.19.8568.View ArticlePubMedPubMed CentralGoogle Scholar
- Tishkoff SA, Varkonyi R, Cahinhinan N, Abbes S, Argyropoulos G, Destro-Bisol G, Drousiotou A, Dangerfield B, Lefranc G, Loiselet J, Piro A, Stoneking M, Tagarelli A, Tagarelli G, Touma EH, Williams SM, Clark AG: Haplotype diversity and linkage disequilibrium at human G6PD: recent origin of alleles that confer malarial resistance. Science. 2001, 293 (5529): 455-462. 10.1126/science.1061573.View ArticlePubMedGoogle Scholar
- Sabeti PC, Reich DE, Higgins JM, Levine HZP, Richter DJ, Schaffner SF, Gabriel SB, Platko JV, Patterson NJ, McDonald GJ, Ackerman HC, Campbell SJ, Altshuler D, Cooper R, Kwiatkowski D, Ward R, Lander ES: Detecting recent positive selection in the human genome from haplotype structure. Nature. 2002, 419 (6909): 832-837. 10.1038/nature01140.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.