- Research article
Population study of 1311 C/T polymorphism of Glucose 6 Phosphate Dehydrogenase gene in Pakistan – an analysis of 715 X-chromosomes
BMC Geneticsvolume 10, Article number: 41 (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.
Two hundred and seventy four G6PD deficient subjects including 235 males and 39 females agreed to participate in the study. Of them 563 C-T mutation was observed in 187 males (79.6%) and in 27 females (69.2%). Thus, cumulative gene frequency for G6PD Mediterranean was 0.74 in both sexes detected in 234 of 313 X-chromosomes. Associated 1311 T allele gene frequency was 60/234 or 0.26. Gene sequencing was performed on 21 G6PD deficient DNA samples (three with 1311T and 18 with 1311C) which showed additional 93T/C polymorphism in intron 11: seven had 93 C and 14 had 93 T. Only two cases of combined mutations of 1311C-T and 93 T-C were observed. Both were males and their ages were 4 and 8 days at the time of diagnosis. Their associated G6PD mutation, if any remained undetected. RFLP-PCR results for 563 C-T and 1311 C/T are shown in figures 1 &2 and demographics and gender based results are detailed in Table 1.
The X-chromosomes of 400 unrelated subjects belonging to five major and 11 minor ethnicities were analyzed for 1311C/T polymorphism of G6PD gene. All except two were males with age ranging from 17 to 52 years (median ± SD 27.0 ± 7.0). G6PD levels of the participants were not performed. It was assumed based on clinical history that the enrolled subjects were not suffering from G6PD deficiency. Moreover, frequency of G6PD deficiency noted in general population of the country is 1.8–8.6% [22–30]. There fore, it was expected that the results would not be skewed even if deficient samples were included in the study. The distribution of 1311 T with respect to various ethnicities in Pakistan is summarized in Table 2 which demonstrated striking similarities among various ethnic groups living in Pakistan. It was observed that the polymorphism was frequent with an incidence of 0.28–0.38 in five main ethnic groups. Highest prevalence was observed in Saraiki subjects with a frequency of 0.5. In contrast, no polymorphism was observed in local Hindu population.
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 the silent mutation within normal population is summarized in Table 3. The reported frequency of 1311T is surprisingly high (0.45) among Asian Indians compared to other populations of the world . Central and South America, Jews, Sicilians and Africans had considerably low gene frequency for this silent mutation ranging from 0.1 to 0.25 . Beutler and Kuhl in 1990 reported frequency of 0.05 for the same in Orientals . Recently the haplotype is recognized with increasing frequency in Chinese population [38, 40] but the same is not true for Burmese and Mons . Variable frequencies have been reported from Bahrain  Saudi Arabia  and Iran  as 0.30, 0.14 and 0.22 respectively. Similar to Americans, English (Caucasians) also showed a lower frequency of 0.137 for this polymorphism . We found that haplotype1311T is prevalent extensively in our various ethnic groups with the overall frequency of 0.32 in general population.
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.
Beutler E: The genetics of glucose-6-phosphate dehydrogenase deficiency. Semin Hematol. 1990, 27 (2): 137-164.
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.
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.
Beutler E: G6PD deficiency. Blood. 1994, 84 (11): 3613-3636.
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.
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.
Fialkow PJ: Clonal origin of human tumors. Biochim Biophys Acta. 1976, 458 (3): 283-321.
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.
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.
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.
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.
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.
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.
Newcomb L: The Islamic Republic of Pakistan: country profile. Int Demogr. 1986, 5 (7): 1-8.
Bernhard W: Ethnogenesis of South Asia with special reference to India. Anthropol Anz. 1983, 41 (2): 93-110.
Birdwood: Review: A History of the Pathans: Review. The Geographical Journal. 1959, 125 (3/4): 414-10.2307/1791129.
Lapidus IM: A history of Islamic societies. 2002, Cambridge; New York: Cambridge University Press
Black C: Pakistan. The people. 2003, New York: Crabtree Pubb. Co
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.
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.
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.
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.
McCurdy PR, Mahmood L: Red cell glucose-6-phosphate dehydrogenase deficiency in Pakistan. J Lab Clin Med. 1970, 76 (6): 943-948.
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.
Bollinger RC: Glucose 6 phosphate dehyrogenase deficiency in Lahore area. Pak J Med Res. 1985, 24: 85-87.
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.
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.
Khan M: Glucose 6 phosphate dehydrogenase deficiency in adults. J Coll Physicians Surg Pak. 2004, 14 (7): 400-403.
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.
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.
Wirsing RG: The Baluchis and Pathans. London. 1987
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.
Papiha SS, Deka R, Chakraborty R: Genomic diversity: applications in human population genetics. 1999, New York: Kluwer Academic/Plenum Publishers
Wilber DN: The nations of Asia. 1966, New York: Hart Pub. Co
Rengel M: Pakistan: a primary source cultural guide. 2003, New York: Rosen Pub. Group
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.
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.
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.
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.
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.
Nabeel Al Momen SSAA, Ahmed Al Alawi A: Molecular Homogenity of G6PD deficiency. Bahrain Medical Bulletin. 2004, 26: 1-7.
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.
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.
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.
BM designed the study, analyzed the data and wrote the paper. AN performed bench work, analyze the results and wrote part of paper. TM wrote part of paper and critically evaluated the manuscript. ZAN did bench work and analyzed the results. MK critically evaluated the paper. All authors approved final form of manuscript before submission.