Population study of 1311 C/T polymorphism of Glucose 6 Phosphate Dehydrogenase gene in Pakistan – an analysis of 715 X-chromosomes
© Moiz et al; licensee BioMed Central Ltd. 2009
Received: 2 April 2009
Accepted: 30 July 2009
Published: 30 July 2009
Nucleotide 1311 polymorphism at exon 11 of G6PD gene is widely prevalent in various populations of the world. The aim of the study was to evaluate 1311 polymorphism in subjects carrying G6PD Mediterranean gene and in general population living in Pakistan.
Patients already known to be G6PD deficient were tested for 563C-T (G6PD Mediterranean) and 1311 C-T mutation through RFLP based PCR and gene sequencing. A control group not known to be G6PD deficient was tested for 1311C/T only.
C-T transition at nt 1311 was detected in 60/234 X-chromosomes with 563 C-T mutation (gene frequency of 0.26) while in 130 of normal 402 X-chromosomes (gene frequency of 0.32).
We conclude that 1311 T is a frequent polymorphism both in general populations and in subjects with G6PD Mediterranean gene in Pakistan. The prevalence is higher compared to most of the populations of the world. The present study will help in understanding genetic basis of G6PD deficiency in Pakistani population and in developing ancestral links of its various ethnic groups.
Glucose 6-phosphate dehydrogenase (G6PD) deficiency is one of the most frequent red cell enzymopathies affecting some 400 million people globally. G6PD gene located on human chromosome Xq28 , is associated with at least 140 distinct disease producing mutations . These mutations are significant in causing hemolysis, favism, neonatal hyerpbilirubinemia . Besides the disease producing mutations, numerous silent mutations unassociated with any conformational change in enzyme structure have also been demonstrated. Such polymorphisms have greatly helped the geneticists for evaluating clonality of various tumors including uterine myomas , lymphomas  and other tumors . Further to it, silent mutations of G6PD gene have proved to be invaluable tools in the study of human genome by providing markers for gene mapping and linkage analysis. There is an extensive list of G6PD genetic polymorphism that has been analyzed so far including IVS5 at NT 611C-G, IVS2 NT 9722 A-G, IVS11 NT 93 T-C, IVS7 NT 175 C-T and IVS8 NT 163 C-T . Unfortunately, many of these polymorphisms had limited utilization as were exclusively observed in African populations.
Of special interest is the polymorphism seen at 1311 nucleotide of exon 11 of G6PD gene. This is widely prevalent in multiethnic non African population of the world and is an important marker in formal genetic analysis of X-linked loci. The nucleotide produces two alleles (1311 C and 1311T) that are observed in both G6PD replete and deficient populations. Polymorphism of 1311 T has been studied extensively in several populations [12, 13]. It does not abolish or create any specific site for enzyme restriction and therefore can be detected by creating mismatched antisense primer which produces a site for Bcl I digestion when 1311T is present .
Pakistan with a population of over 165 million people is a home to 18 ethnic groups speaking 60 different languages . The ethnic diversity of this region is attributed to repeated invasions by Aryans , Macedonians , Arabs and Mongols in its long history . Majority of its people belong to five major ethnic groups: Punjabis, Sindhis, Pathans, Baluch (Balochis) and Mohajirs . Several distinct ethnic groups exist in each major group which has generated complex social, cultural and lingual blend.
The primary objectives were to evaluate NT 1311 polymorphism of G6PD gene in various ethnic groups living in Pakistan and to evaluate its association with 563C-T mutation in G6PD deficient subjects. The secondary objective was to compare the results of 1311 polymorphism with that reported from other populations of the world.
The Aga Khan University Hospital Clinical Laboratory uses brilliant cresyl blue dye test (Trinity biotech plc, Wicklow, Ireland) for screening individuals with G6PD deficiency. Subjects reported as G6PD deficient from January 2006 to December 2008 were identified using a computerized institutional database search utilizing International Classification of Disease (v9, American Health Information Management Association, AHIMA, USA). They were contacted via telephone and were requested for participation in the study.
In addition, a control group consisting of healthy subjects who visited the hospital for blood donation was included. These subjects were classified in to various ethnic groups based on subject's self classification. They were verbally scanned through an in-house questionnaire to rule out the possibility of any disorder but were not tested for G6PD deficiency. As the frequency of consanguineous marriages is reported as 31.1% to 58.7% from different places of the country [19, 20], some ethnic mix is expected in our cultural groups.
From each participant, 5 ml of blood was collected in an EDTA containing tube afterward white blood cells were separated and stored for DNA extraction according to manufacturer's instructions.
DNA extraction and PCR/RFLP for screening G6PD 563-T and 1311/T mutations
Genomic DNA was isolated from peripheral blood leucocytes by Wizard® Genomic DNA Purification Kit (Wisconsin, USA). PCR/RFLP technique was used for the detection of G6PD 563C-T and 1311C/T with primers published elsewhere . Each PCR reaction in 25 μl volume contained 34 mM Tris-HCl pH 8.8, 8.3 mM ammonium sulfate, 1.5 mM MgCl2, 0.2 mM of each dNTP, 75 ng of each primer, 2 U of Taq polymerase (Promega, Wisconsin, USA) and 200 ng of genomic DNA. After 30 cycles consisting of 1 min at 94°C, 30 sec at 58/62°C (respectively for1311C/T and 563 C-T) and 40 sec at 72°C, PCR products were digested overnight at 37°C with either Mbo II or Bcl I restriction enzymes. The digested products were subsequently analyzed electrophorically; Mbo II produced 277, 119 and 110 bp fragments in G6PD Mediterranean samples whereas 377 and 119 bp fragments were observed in normal subjects. DNA samples containing C at position 1311 showed 23 and 203 bp fragments with Bcl I digestion whereas replacement of C by T at this nucleotide position resulted in a 180 bp product.
Undigested PCR products were sent to Macrogen® (Seoul, Korea) for gene sequencing.
The frequency of 1311 polymorphism was evaluated and compared in two groups; (a) patients carrying G6PD Mediterranean gene (b) healthy individuals.
All calculations were computed using statistical package SPSS version 16.0 (SPSS inc.; Chicago, IL, USA). Pearson's chi-square test was used for the comparison of statistical difference between the groups and the threshold of significance (p-value) was set at 0.05.
The study was approved by institutional ethical review committee and informed consent was obtained from all subjects prior to their enrollment.
Analysis of X-chromosomes from G6PD deficient subjects with 563 C-T mutation (n = 214).
mean ± SD
16.8 ± 16.6
8.2 ± 12.7
9.3 ± 13.4
Demographic details and gene frequency for 1311 C/T polymorphism of healthy subjects with respect to their ethnicity (n = 400, M: 398, F: 2)
No. of X Chromosomes
n (gene frequency)
N (gene frequency)
No significant statistical differences were assessed in 1311 T polymorphism in compared groups: Pathans and Jews (p-value 0.29), Mohajirs and Indians (p-value 0.32) and Baluchis and Iranians (p-value 0.71).
We performed an extensive analysis of 1311 C/T polymorphism of G6PD gene within normal Pakistani population and in subjects with G6PD deficiency. It was found that this polymorphism is quite frequent in normal and individuals with G6PD Mediterranean gene. The work enabled us to understand the genetic basis of G6PD deficiency in Pakistani population and to compare T allele frequency with various world populations including those that are considered traditionally to be origin or source of various ethnic groups residing in Pakistan.
Pakistan is a multilingual and multiethnic nation. Its various ethnic groups are unequally distributed through out the country. Punjabis and Pathans which are mainly located in the North Pakistan constitute 60% and 13% of the population respectively  while Sindhis and Baluchis occupy Southern Pakistan. There is considerable ethnic admixture among Punjabis because most of the incoming invaders in the past 2000 years passed through Punjab and mixed with previous settlers. In contrast, Pathans are claimed to be descendants of European soldiers who accompanied Alexander the Great  and are also considered to be the lost tribe of Jews . Sindhis and Baluchis that are main occupants of Southern Pakistan comprise respectively of 13 and 8% of Pakistani population[34, 35]. While former have mixed ethnic background, the Baluchis are believed to have originated in Aleppo (Syria) who migrated via Iran to this land . Mohajirs constitute some 8% of population and are a group of people who migrated to Pakistan from near east – the present day India after its independence as a separate state in 1947 . They lack the tribal based ethnic identity in contrast to the rest of Pakistanis.
Polymorphism of 1311C/T with and without 563C-T mutation in Pakistani G6PD deficient population
C-T transition at NT 1311 was detected in 60/234 X-chromosomes carrying G6PD Mediterranean gene (equaling a gene frequency of 0.26). The frequency of this polymorphism was not statistically different for patients whose DNA did not manifest G6PD 563C-T as was detected in 11/72 X-chromosomes (0.15 gene frequency; p-value 0.189). A lower incidence of 0.13 was observed by Saha et al in 1994 for this silent mutation among subjects with G6PD Mediterranean in Pakistan.
Association of 1311 C/T and 93 T/C polymorphisms
Recently, 1311 C/T polymorphism was observed in association with another silent mutation 93 T/C at intron 11 in South East Asians [36, 37]and Chinese populations [38, 39]. The association of two is significant in reducing G6PD enzyme activity and hence has clinical implications . Unfortunately, gene sequencing results were available for limited number of G6PD deficient patients in the present study which showed a gene frequency of 0.66 for combined mutations of 1311C-T and 93 T-C.
G6PD NT1311 Polymorphism in various ethnicities in Pakistan
G6PD NT 1311 C-T was observed uniformly in various ethnic groups in Pakistan with an overall frequency of 0.32. Different results seen in local Hindu and Sariki populations may be secondary to small sample size demanding further evaluation. A lower though not statistically significant frequency for T allele (0.20) was observed by Saha N et al in 1994 for Pashtoons and Punjabi population  in contrast to 0.3 in the present work for the same population (p-value 0.133).
Comparison of NT 1311 polymorphism with global data
Global Distribution of T allele gene in various populations
No of X-chromosomes studied
1311 T gene frequency
Assessment of 1311 polymorphism with Source Population
Qamar R et al in 2002 through their extensive work on Y-chromosomal DNA supported origin of Parsis in Iran, Hazaras from Genghis Khan's army and Negroid Makrani in Africa . However, they denied the Syrian, Tibetan, Greek or Jewish origin for Baluch, Balti and Pathan population respectively. Subsequent work revealed limited Greek contribution to Pakistani Pathan Population . Our work was limited in investigating a single allele which by no means allows us to compare the populations. Also, statistical comparison between two unrelated studies can not be conclusive. However, it was interesting to note that Pathans and Jews, Mohajir and Indians and Baluch and Irani revealed insignificant statistical differences in 1311 T frequency. Other groups like Punjabis and Sindhis were not compared with any population because of their 'mix' origin.
In this study, we report the presence of silent mutation in normal G6PD alleles in a frequency of 0.32 which apart from Asian Indians  is highest reported from any population.
We conclude that G6PD nucleotide 1311 T polymorphism is relatively common in Pakistani population both with and without G6PD Mediterranean. It is expected that the work will help in understanding genetic basis of G6PD deficiency in Pakistani population and establishing ancestral links between various ethnicities residing in this region.
This work was supported by grants from The Aga Khan University, Grant No. SM# 050905. We appreciate the technical and moral support provided by Dr Rumina Hasan, our departmental chair and Dr Anwar Siddiqui-Associate Dean of Research office. The authors are grateful to the blood donors for providing their blood samples for DNA analysis.
- Beutler E: The genetics of glucose-6-phosphate dehydrogenase deficiency. Semin Hematol. 1990, 27 (2): 137-164.PubMedGoogle Scholar
- Pai GS, Sprenkle JA, Do TT, Mareni CE, Migeon BR: Localization of loci for hypoxanthine phosphoribosyltransferase and glucose-6-phosphate dehydrogenase and biochemical evidence of nonrandom X chromosome expression from studies of a human X-autosome translocation. Proc Natl Acad Sci USA. 1980, 77 (5): 2810-2813. 10.1073/pnas.77.5.2810.PubMed CentralView ArticlePubMedGoogle Scholar
- Beutler E, Vulliamy TJ: Hematologically important mutations: glucose-6-phosphate dehydrogenase. Blood Cells Mol Dis. 2002, 28 (2): 93-103. 10.1006/bcmd.2002.0490.View ArticlePubMedGoogle Scholar
- Beutler E: G6PD deficiency. Blood. 1994, 84 (11): 3613-3636.PubMedGoogle Scholar
- Linder D, Gartler SM: Glucose-6-phosphate dehydrogenase mosaicism: utilization as a cell marker in the study of leiomyomas. Science. 1965, 150 (692): 67-69. 10.1126/science.150.3692.67.View ArticlePubMedGoogle Scholar
- Beutler E, Collins Z, Irwin LE: Value of genetic variants of glucose-6-phosphate dehydrogenase in tracing the origin of malignant tumors. N Engl J Med. 1967, 276 (7): 389-391.View ArticlePubMedGoogle Scholar
- Fialkow PJ: Clonal origin of human tumors. Biochim Biophys Acta. 1976, 458 (3): 283-321.PubMedGoogle Scholar
- Pack SD, Kolonin MG, Borodin PM, Searle JB, Serov OL: Gene mapping in the common shrew (Sorex araneus; Insectivora) by shrew-rodent cell hybrids: chromosome localization of the loci for HPRT, TK, LDHA, MDH1, G6PD, PGD, and ADA. Mamm Genome. 1995, 6 (11): 784-787. 10.1007/BF00539004.View ArticlePubMedGoogle Scholar
- De bruyn A, Raeymaekers P, Mendelbaum K, Sandkuijl LA, Raes G, Delvenne V, Hirsch D, Staner L, Mendlewicz J, Van Broeckhoven C: Linkage analysis of bipolar illness with X-chromosome DNA markers: a susceptibility gene in Xq27-q28 cannot be excluded. Am J Med Genet. 1994, 54 (4): 411-419. 10.1002/ajmg.1320540423.View ArticlePubMedGoogle Scholar
- Fey MF, Wainscoat JS, Mukwala EC, Falusi AG, Vulliamy TJ, Luzzatto L: A PvuII restriction fragment length polymorphism of the glucose-6-phosphate dehydrogenase gene is an African-specific marker. Hum Genet. 1990, 84 (5): 471-472. 10.1007/BF00195822.View ArticlePubMedGoogle 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 USA. 1991, 88 (19): 8568-8571. 10.1073/pnas.88.19.8568.PubMed CentralView ArticlePubMedGoogle Scholar
- Beutler E, Kuhl W: The NT 1311 polymorphism of G6PD: G6PD Mediterranean mutation may have originated independently in Europe and Asia. Am J Hum Genet. 1990, 47 (6): 1008-1012.PubMed CentralPubMedGoogle Scholar
- Kurdi-Haidar B, Mason PJ, Berrebi A, Ankra-Badu G, al-Ali A, Oppenheim A, Luzzatto L: Origin and spread of the glucose-6-phosphate dehydrogenase variant (G6PD-Mediterranean) in the Middle East. Am J Hum Genet. 1990, 47 (6): 1013-1019.PubMed CentralPubMedGoogle Scholar
- Newcomb L: The Islamic Republic of Pakistan: country profile. Int Demogr. 1986, 5 (7): 1-8.Google Scholar
- Bernhard W: Ethnogenesis of South Asia with special reference to India. Anthropol Anz. 1983, 41 (2): 93-110.PubMedGoogle Scholar
- Birdwood: Review: A History of the Pathans: Review. The Geographical Journal. 1959, 125 (3/4): 414-10.2307/1791129.View ArticleGoogle Scholar
- Lapidus IM: A history of Islamic societies. 2002, Cambridge; New York: Cambridge University PressGoogle Scholar
- Black C: Pakistan. The people. 2003, New York: Crabtree Pubb. CoGoogle Scholar
- Hussain R, Bittles AH: The prevalence and demographic characteristics of consanguineous marriages in Pakistan. Journal of biosocial science. 1998, 30 (2): 261-275. 10.1017/S0021932098002612.View ArticlePubMedGoogle Scholar
- Wahab A, Ahmad M: Biosocial perspective of consanguineous marriages in rural and urban Swat, Pakistan. Journal of biosocial science. 1996, 28 (3): 305-313. 10.1017/S0021932000022379.View ArticlePubMedGoogle Scholar
- Saha N, Ramzan M, Tay JS, Low PS, Basair JB, Khan FM: Molecular characterisation of red cell glucose-6-phosphate dehydrogenase deficiency in north-west Pakistan. Hum Hered. 1994, 44 (2): 85-89. 10.1159/000154196.View ArticlePubMedGoogle Scholar
- Ronald AR, Underwood BA, Woodward TE: Glucose-6-phosphate dehydrogenase deficiency in Pakistani males. Trans R Soc Trop Med Hyg. 1968, 62 (4): 531-533. 10.1016/0035-9203(68)90139-9.View ArticlePubMedGoogle Scholar
- McCurdy PR, Mahmood L: Red cell glucose-6-phosphate dehydrogenase deficiency in Pakistan. J Lab Clin Med. 1970, 76 (6): 943-948.PubMedGoogle Scholar
- Hashmi JA, Farzana F, Ahmed M: Abnormal hemoglobins, thalasemia trait & G6PD deficiency in young Pakistani males. J Pak Med Assoc. 1976, 26 (1): 2-4.PubMedGoogle Scholar
- Bollinger RC: Glucose 6 phosphate dehyrogenase deficiency in Lahore area. Pak J Med Res. 1985, 24: 85-87.Google Scholar
- Khattak MFDM, Saleem M: The prevalence of glucose 6 phosphate dehydrogenase deficiency in Nortehrn Pakisatn. Pak Armed Forces MEd J. 1992, 42 (1): 35-38.Google Scholar
- Bouma MJ, Goris M, Akhtar T, Khan N, Khan N, Kita E: Prevalence and clinical presentation of glucose-6-phosphate dehydrogenase deficiency in Pakistani Pathan and Afghan refugee communities in Pakistan; implications for the use of primaquine in regional malaria control programmes. Trans R Soc Trop Med Hyg. 1995, 89 (1): 62-64. 10.1016/0035-9203(95)90661-4.View ArticlePubMedGoogle Scholar
- Khan M: Glucose 6 phosphate dehydrogenase deficiency in adults. J Coll Physicians Surg Pak. 2004, 14 (7): 400-403.PubMedGoogle Scholar
- Khan TAAS, Anwar M, Ayub M: The frequancy of glucose 6 phosphate deficeint in Punjabis and Pathans. J Post Grad Med Inst. 2004, 18 (4): 592-597.Google Scholar
- Ali N, Anwar M, Ayyub M, Bhatti FA, Nadeem M, Nadeem A: Frequency of glucose-6-phosphate dehydrogenase deficiency in some ethnic groups of Pakistan. J Coll Physicians Surg Pak. 2005, 15 (3): 137-141.PubMedGoogle Scholar
- Wirsing RG: The Baluchis and Pathans. London. 1987Google Scholar
- Firasat S, Khaliq S, Mohyuddin A, Papaioannou M, Tyler-Smith C, Underhill PA, Ayub Q: Y-chromosomal evidence for a limited Greek contribution to the Pathan population of Pakistan. Eur J Hum Genet. 2007, 15 (1): 121-126. 10.1038/sj.ejhg.5201726.PubMed CentralView ArticlePubMedGoogle Scholar
- Papiha SS, Deka R, Chakraborty R: Genomic diversity: applications in human population genetics. 1999, New York: Kluwer Academic/Plenum PublishersView ArticleGoogle Scholar
- Wilber DN: The nations of Asia. 1966, New York: Hart Pub. CoGoogle Scholar
- Rengel M: Pakistan: a primary source cultural guide. 2003, New York: Rosen Pub. GroupGoogle Scholar
- Nuchprayoon I, Louicharoen C, Charoenvej W: Glucose-6-phosphate dehydrogenase mutations in Mon and Burmese of southern Myanmar. Journal of human genetics. 2008, 53 (1): 48-54. 10.1007/s10038-007-0217-3.View ArticlePubMedGoogle Scholar
- Matsuoka H, Thuan DT, van Thien H, Kanbe T, Jalloh A, Hirai M, Arai M, Dung NT, Kawamoto F: Seven different glucose-6-phosphate dehydrogenase variants including a new variant distributed in Lam Dong Province in southern Vietnam. Acta medica Okayama. 2007, 61 (4): 213-219.PubMedGoogle Scholar
- Yu GL, Jiang WY, Du CS, Chen LM, Lin QD, Tian QH, Zeng JB, Li SG: [Identification of G6PD gene variants from Hakka population in Guangdong province]. Zhonghua yi xue yi chuan xue za zhi = Zhonghua yixue yichuanxue zazhi = Chinese journal of medical genetics. 2004, 21 (5): 448-451.PubMedGoogle Scholar
- Yu GL, Jiang WY, Du CS, Lin QD, Chen LM, Tian QH, Li SG, Zeng JB: [Complex mutations of 1311 C-->T in exon 11 and 93 T-->C in intron 11 in G6PD gene]. Zhonghua xue ye xue za zhi = Zhonghua xueyexue zazhi. 2004, 25 (10): 610-612.PubMedGoogle Scholar
- Yang Z, Chu J, Ban G, Huang X, Xu S, Li M: [The genotype analysis of glucose-6-phosphate dehydrogenase deficiency in Yunnan province]. Zhonghua yi xue yi chuan xue za zhi = Zhonghua yixue yichuanxue zazhi = Chinese journal of medical genetics. 2001, 18 (4): 259-263.PubMedGoogle Scholar
- Nabeel Al Momen SSAA, Ahmed Al Alawi A: Molecular Homogenity of G6PD deficiency. Bahrain Medical Bulletin. 2004, 26: 1-7.Google Scholar
- Mortazavi Y, Chopra R, Gordon-Smith EC, Rutherford TR: Frequency of the G6PD nt 1311 C/T polymorphism in English and Iranian populations: relevance to studies of X chromosome inactivation. J Med Genet. 1997, 34 (12): 1028-1029. 10.1136/jmg.34.12.1028.PubMed CentralView ArticlePubMedGoogle Scholar
- Qamar R, Ayub Q, Mohyuddin A, Helgason A, Mazhar K, Mansoor A, Zerjal T, Tyler-Smith C, Mehdi SQ: Y-chromosomal DNA variation in Pakistan. Am J Hum Genet. 2002, 70 (5): 1107-1124. 10.1086/339929.PubMed CentralView 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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.