A population survey of the glucose-6-phosphate dehydrogenase (G6PD) 563C>T (Mediterranean) mutation in Afghanistan.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is a common inherited enzyme defect and an important problem in areas with Plasmodium vivax infection because of the risk of haemolysis following administration of primaquine to treat the liver forms of the parasite. We undertook a genotypic survey of 713 male individuals across nine provinces of Afghanistan in which malaria is found, four in the north and five in the east. RFLP typing at nucleotide position 563 detected 40 individuals with the Mediterranean mutation 563C>T, an overall prevalence of 5.6%. This varied according to self-reported ethnicity, with prevalence in the Pashtun/Pashai group of 33/369 (8.9%) compared to 7/344 individuals in the rest of the population (2.0%; p<0.001, Chi-squared test). Multivariate analysis of ethnicity and geographical location indicated an adjusted odds ratio of 3.50 (95% CI 1.36-9.02) for the Pashtun/Pashai group, while location showed only a trend towards higher prevalence in eastern provinces (adjusted odds ratio = 1.73, 0.73-4.13). Testing of known polymorphic markers (1311C>T in exon 11, and C93T in intron XI) in a subset of 82 individuals wild-type at C563 revealed a mixture of 3 haplotypes in the background population and was consistent with data from the 1000 Genomes Project and published studies. By comparison individuals with G6PD deficiency showed a highly skewed haplotype distribution, with 95% showing the CT haplotype, a finding consistent with relatively recent appearance and positive selection of the Mediterranean variant in Afghanistan. Overall, the data confirm that the Mediterranean variant of G6PD is common in many ethnic groups in Afghanistan, indicating that screening for G6PD deficiency is required in all individuals before radical treatment of P. vivax with primaquine.

Introduction

G6PD deficiency is the most common enzyme defect of humans [1]. Many G6PD-deficient individuals remain asymptomatic and unaware of their enzyme deficiency, with no reduction in life expectancy [2]. Clinical manifestations include acute haemolysis following oxidative stress triggered by drugs (such as primaquine, dapsone, chloramphenicol and ciprofloxacin), infection, or the ingestion of fava beans. G6PD deficiency is also associated with an increased risk of neonatal jaundice, and rarely with chronic non-spherocytic haemolytic anaemia and gallstones. The degree of enzyme deficiency and severity of clinical complications depend on the exact G6PD mutation involved. The Mediterranean and certain southeast Asian variants are associated with less than 10% residual enzyme activity (Class II deficiency) [3] and generally more severe clinical manifestations than the African A- form which provides 10–60% of residual activity (Class III). Since the G6PD gene is found on the X-chromosome, there is a higher risk of haemolytic crisis in males (hemizygous) and homozygous females with mutations than in heterozygous females, although heterozygotes are also at some risk due to X-chromosome inactivation [4].

The term ‘Mediterranean’ is used to describe the 563C>T mutation in exon 6 of the human G6PD gene (changing a serine to a phenylalanine residue at position 188 of the protein product), reflecting the first description of this variant in countries such as Italy and Cyprus. However the mutation is also the dominant molecular determinant of G6PD deficiency in many Middle Eastern countries and the Indian subcontinent and has been documented as far east as China [5], Malaysia [6] and Singapore [7]. In south Asia this variant of G6PD deficiency is of practical relevance as there is substantial geographical overlap with malaria. The administration of primaquine to clear the hypnozoite forms of P. vivax [8] (radical cure) is compromised by the lack of readily available tests to assess G6PD activity in most locations. G6PD deficiency hence represents a major obstacle to malaria control efforts in Afghanistan and the wider region.

Before commencement of this study, there had been no published studies of G6PD deficiency from within the borders of Afghanistan. The 563C>T mutation has been shown to be the specific molecular cause of most cases of phenotypic G6PD deficiency in two locations in Pakistan [9], [10], and has also been shown to be responsible for the majority of cases of phenotypic G6PD deficiency in Pashtun refugees in Pakistan originating from Afghanistan [11]. G6PD deficiency appears to be at relatively low prevalence in countries north of Afghanistan [12], [13]. It was therefore of interest to investigate the prevalence of the 563C>T allele in Afghanistan itself in both the northern and eastern regions of the country. We also examined two silent polymorphic markers in order to explore the evolutionary history of the Mediterranean mutation in Afghanistan.

Materials and Methods

Ethical Approval

Ethical approval for the study was obtained from the Ethics Committee of the Faculty of Tropical Medicine, Mahidol University, Thailand and the Institutional Review Board of the Afghanistan National Public Health Institute. Participants or their parents (in the case of children) provided written informed consent to participate in the study.

Sites

The study was undertaken in one urban centre in each of nine provinces: Jalalabad (Nangarhar province), Mehtherlam (Laghman), Asadabad (Kunar), Maimana (Faryab), Taloqan (Takhar), Imamsahib (Kunduz), Pulikhumri (Baghlan), Pulialam (Lowgar) and the capital, Kabul (Kabul province). Provinces were chosen on the basis of practical accessibility combined with having a high risk of vivax malaria (>1 case/1000 population/year), the exceptions being Kabul and Lowgar which have low transmission [14].

Samples

Only males were studied given that the problem of haemolysis after primaquine administration more commonly affects males, and determination of haplotypes is relatively straightforward with single X chromosomes. A convenience approach to sampling was used for practical reasons related to security and resources. EDTA blood was obtained from male adults and children attending outpatient medical laboratories or neighbouring vaccine administration clinics in 2009. No individuals with malaria or infants less than two months of age were included so that malaria and neonatal jaundice were excluded from consideration. There is extensive evidence showing that G6PD deficiency is not associated with other commonly encountered medical conditions that might lead to bias in the selection of cases [1], [2], [15], [16]. Participants or their parents provided written informed consent to participate in the study. Volunteers were asked to state their ethnic group. Samples were transferred to filter-paper blood spots and stored in plastic zip-lock bags with silica gel at room temperature before transport to Bangkok for genotyping.

Genomic DNA Extraction

Genomic DNA was extracted from dried blood spots using the QIAamp® DNA Mini Kit (QIAGEN, Germany), following the manufacturer’s instructions. Eluted genomic DNA samples were frozen at −20°C until amplification via PCR.

PCR-RFLP for Detection of G6PD Mediterranean Variant

All samples were subjected to PCR-RFLP to assess the Mediterranean variant in exon 6, based on the protocol of Samilchuk et al. [17] with modifications including use of primers F14948 and R15158 at 250 nM, a volume of 2 µl of each genomic DNA template, a final reaction volume of 30 µl containing 20 mM Tris-HCl (pH 8.4), 50 mM KCl, 1 mM MgCl2, 125 µM 4-deoxynucleotide triphosphate (dNTPs), and 0.05 units Platinum®Taq DNA polymerase (Invitrogen, Brazil). After pre-denaturation at 95°C for 5 min, 45 PCR cycles were undertaken involving denaturation at 94°C for 1 min, annealing at 55°C for 1 min and extension at 72°C for 1 min, with post-extension at 72°C for 7 min, using a MyCycler™ thermal cycler (Bio-Rad Laboratories, U.S.A.). To examine the 563C>T mutation (Mediterranean variant, Table 1), 10 µl of each PCR product was digested with 10 units of the restriction enzyme MboII (New England Biolabs Inc.) at 37°C for 3 hours and visualized by 2% agarose gel electrophoresis. A product pattern of 104+98+28 bp indicates G6PD Mediterranean mutation, whereas the digestion pattern 202+28 bp indicates wild-type.

Table 1

G6PD molecular markers examined in the study.

Mutation	SNP code	Location	Ancestral allele	Derived allele	Amino acid effect
c.563C>T	rs5030868	Exon 6	C	T	S188F (G6PD Med)
c.1311C>T	rs2230037	Exon 11	C	T	Y437Y
IVSXI C93T*	rs2071429	Intron 11	C	T	Non-coding

Numbers refer to the short isoform version of G6PD. *Commonly referred to as T93C, but correctly annotated as C93T in NCBI dbSNP, consistent with 1000 Genomes and Pan spp. data. This mutation could also be formally referred to as c.1365-13C>T.

Polymorphic Markers

In order to study polymorphic markers relevant to the Asian setting (1311C>T [18] and IVSXI C93T [19], Table 1), in all samples with 563C>T and a subset of samples wild-type at C563, PCR and sequencing of exons 11 and 12 including introns 11 and 12 was undertaken using the protocol of Tang et al. [20]; the forward primer was GAAGCCGGGCATGTTCTTCAAC and the reverse CCAGGGCTCAGAGCTTGTG. All amplification fragments spanning exons 11 and 12 were subjected to gel electrophoresis and purified using FavorPrep™ GEL/PCR Purification kit (Favorgen, Austria). Purified PCR products were assessed by gel electrophoresis and PCR products more concentrated than 50 ng/µl were sequenced at Macrogen, Republic of Korea and analysed using BioEdit version 7.1.3.0. These molecular data were supplemented with relevant data on the same markers either published or downloaded from the 1000 Genomes Project [21], using chromosomes as the denominator. A literature review was also undertaken examining the 1311C>T polymorphism in the context of the 563C>T mutation (minimum of 10 563C>T chromosomes per study); in virtually all these datasets the IVSXI C93T polymorphism was not examined.

Power and Statistical Analysis

Before the study the available evidence indicated that males from the Pashtun ethnic group had a prevalence of G6PD deficiency of around 10%. Assuming an even mixture of Pashtuns and other groups, a sample size of 512 had 90% power to detect an absolute difference of 7% between Pashtuns and other groups (pooled). In order to ensure sufficient numbers in each of these groups, the sampling was extended to more than 700 patients across nine provinces (covering both northern and eastern regions).

Data were stored and statistical tests undertaken in Stata v12.0.

Results

Population Sample

Samples were obtained from 713 individuals across nine provinces of Afghanistan (Table 2, Figure 1), consisting of four in the north (Baghlan, Faryab, Kunduz and Takhar) and five in the east (Kabul, Kunar, Laghman, Lowgar, Nangarhar), these two regions being separated by the Hindu Kush mountain range. There were eight self-reported groups of ethnicity: Arab, Baluch, Hazara, Pashai, Pashtun, Tajik, Turkmen and Uzbek (Table 2) while in 6 individuals ethnicity was not described. As expected, there was a clear link between province and ethnicity with Pashtuns and Pashais generally concentrated in the eastern provinces and Tajiks, Uzbeks and Turkmen generally found in the north.

Table 2

Number of individuals with the 563C>T mutation among 713 male individuals studied.

	Ethnicity
Province	Arab	Baluch	Hazara	Pashai	Pashtun	Tajik	Turkmen	Uzbek	Unknown	Total
Northern
Baghlan	0	0	0/5	0	0/4	0/16	0	01	0/6	0/32 (0%)
Faryab	0	0	0	0	0/4	0/6	0	4//94	0	4/104 (3.8%)
Kunduz	0/3	0	0/12	0	3/57	0/19	0/9	0/7	0	3/107 (2.8%)
Takhar	0	0	0	0	1/10	1/70	0	0/16	0	2/96 (2.1%)
Eastern
Kabul	0	0/1	0/10	0/1	0/10	0/9	0/1	0/3	0	0/35 (0%)
Kunar	0/1	0	0	3/29	2/32	0/5	0	0	0	5/67 (7.5%)
Laghman	0	0	0	1/10	14/76	0/10	0	0	0	15/96 (15.6%)
Lowgar	0/1	0	0	0	4/74	0/25	0	0	0	4/100 (4%)
Nangarhar	2/11	0	0	2/9	3/53	0/3	0	0	0	7/76 (9.2%)
Total	2/16(12.5%)	0/1(0%)	0/27(0%)	6/49(12.2%)	27/320(8.4%)	1/163 (0.6%)	0/10(0%)	4/121(3.3%)	0/6(0%)	40/713 (5.6%)

The total number of individuals in each province/ethnicity grouping is presented as denominator (zero indicates that there were no individuals present in a given category).

Figure 1

Location, ethnicity and number of Mediterranean variants for the individuals studied.

Black circles placed within the relevant segment represent individuals of a particular ethnic group with the 563C>T Mediterranean mutation. The number of individuals sampled in each province is shown next to the province name.

Mediterranean Variant

The Mediterranean variant (563C>T mutation) of G6PD was found in 40 of the 713 (5.6%) sampled male individuals (Table 2). Prevalence was clearly associated with both province (Pearson chi2 = 29.4070, p<0.001) and ethnic group (21.8520, p = 0.005). Grouping the data according to region, the prevalence in northern and eastern provinces was 9/339 (2.7%) and 31/374 (8.3%) respectively. Grouping according to self-reported ethnicity, Pashtuns/Pashais had a G6PD Mediterranean variant prevalence of 33/369 (8.9%) compared to the rest of the population (7/344 = 2.0%). A multivariate analysis undertaken using these methods of grouping indicated a dominant role for ethnicity over geographical location (adjusted odds ratio = 3.50 (95% CI 1.36–9.02) for the Pashtun/Pashai group); location in eastern provinces showed only a trend towards higher prevalence (adjusted odds ratio = 1.73, 0.73–4.13). There was a non-significant trend only towards a lower prevalence in Tajiks compared to Uzbeks (Fisher’s exact p = 0.167).

Polymorphic Markers and Haplotypes

Two polymorphic markers were examined in all 40 individuals with the Mediterranean mutation and 82 individuals wild-type at C563. These were the synonymous coding sequence polymorphism 1311C>T [18], [22] and a non-coding polymorphism present in the eleventh intron, IVSXI C93T [19]. Both markers are thought to be of no functional significance. Examination of corresponding human 1000 Genomes Project data [23] and available published data from India [24] illustrated that, contrary to many published studies, the ancestral haplotype for the two markers is likely to be C1311/IVSXI C93, (CC haplotype); this is based on the predominance of the CC haplotype in African populations as well as its presence in Pan spp.

The 2-marker haplotypes derived from examination of these positions in Afghanistan are shown in Table 3 and, along with the 1000 Genomes data, in Figure 2. In the 82 C563 individuals (wild-type at the Mediterranean locus) there was a mixture of CT (46) and TC (33) haplotypes, with only 3 ancestral CC haplotypes, consistent with previous data from India [24]. There was no association with ethnic group; among the 38 Pashtun/Pashai individuals studied 1, 21 and 16 had the CC, CT and TC haplotypes respectively, with a very similar pattern in the rest of the population (2, 25 and 17 respectively; p = 0.928 in Fisher’s exact test). Two CT individuals were incidentally found to have a novel synonymous C to T mutation at coding position 1398 in exon 12 (preserving a threonine residue) that is not reported in the 1000 Genomes data or in any publication; this SNP therefore appears so far to be private to Afghanistan.

Table 3

Frequency of 1311C>T/IVSXI C93T haplotypes according to allele at 563C>T (Mediterranean).

Allele at 563 Exon 6	Allele at 1311 Exon 11	Allele at IVSXI 93 Intron 11	Number of individuals
C (wild-type)	C	C	3
	C	T	46*
	T	C	33
	T	T	0
T (Med variant)	C	C	0
	C	T	36
	T	C	2
	T	T	0

*Includes 2 individuals with a novel synonymous 1398C>T mutation.

Figure 2

Haplotypes of 1311C>T and IVSXI C93T in G6PD-normal individuals.

Data are expressed as total chromosomes and are derived from the 1000 genomes project [23], one published dataset from India [24](INDIA) and this study (AFG). Population abbreviations: ASW, people with African ancestry in Southwest United States; CEU, Utah residents with ancestry from Northern and Western Europe; CHB, Han Chinese in Beijing, China; CHS, Han Chinese South, China; CLM, Colombians in Medellin, Colombia; FIN, Finnish in Finland; GBR, British from England and Scotland, UK; IBS, Iberian populations in Spain; LWK, Luhya in Webuye, Kenya; JPT, Japanese in Tokyo, Japan; MXL, people with Mexican ancestry in Los Angeles, California; PUR, Puerto Ricans in Puerto Rico; TSI, Toscani in Italia; YRI, Yoruba in Ibadan, Nigeria.

In our samples from Afghanistan, the Mediterranean 563T mutation was almost exclusively associated with the CT haplotype (36/38 cases) with only two TC individuals detected; compared to the C563 (wild-type) population this was a highly significant difference (Fisher’s exact p<0.001). We examined these data in the context of relevant published studies (also summarized previously [25]) describing allele frequencies at 1311C>T in individuals with Mediterranean G6PD deficiency in Asia, with work covering the Middle East [17], [25]–[28], the Indian subcontinent [10], [29]–[31] and Malaysia [6](Figure 3). The data obtained here for Afghanistan match those from north-west Pakistan [10] where the wild-type allele (C1311) is the predominant background in Mediterranean variant individuals, in marked contrast to studies from the Middle East where the 563C>T mutation is strongly associated with a 1311C>T background [22].

Figure 3

Allele frequencies at 1311C>T marker in individuals carrying 563C>T mutation in southwest Asia.

Data are derived from published studies of individuals (expressed as chromosomes) and this study. Countries of study and references are as follows: 1 = Palestine [25], 2 = Saudi Arabia [18], 3 = Iraq (Kurds) [26], 4 = Iraq (Baghdad) [27], 5 = Kuwait [17], 6 = Iran (Kurds) [28], 7 = Afghanistan (this study), 8 = Pakistan (northwest) [10], 9 = Pakistan (Karachi) [30], 10 and 11 = India (Bombay) [29], [31], 12 = Malaysia [6].

Discussion

Overall Findings

This survey was undertaken primarily in order to assess the prevalence of the Mediterranean G6PD variant in a range of provinces within Afghanistan, spanning regions in the north and east of the country where primaquine is an important part of efforts to treat and control P. vivax infection. During this work, an independent study by Leslie et al., taking an analogous approach, but also including phenotypic assessment of G6PD status, was published [32]. The two studies have produced similar findings. Neither study was designed to provide a comprehensive measure of the overall prevalence of G6PD deficiency across Afghanistan, but in both studies the prevalence of the Mediterranean variant 563C>T in males was approximately 5%. In our study the Mediterranean mutation was more common among individuals identifying themselves as Pashtun or Pashai (8.9%); the phenotype-based study reported a prevalence of deficiency of 10.0% in Pashtuns [32]. Individuals reporting themselves to be of Tajik, Uzbek, Hazara or Turkmen ethnicity have lower prevalence of 563C>T in the region of 1–3%. We saw a non-significant trend towards higher prevalence in Uzbeks compared to Tajiks that is consistent with work from within Pakistan [33] and Tajikistan [12], but there was no evidence for this in the other large survey from within Afghanistan [32].

The work reported here involved only genotypic assessment of the Mediterranean allele, and so did not assess whether there are other molecular causes of G6PD deficiency in Afghanistan. In the phenotype-based survey from Afghanistan [32], only 2 of the 44 individuals found to be G6PD deficient by the fluorescent spot test (and successfully genotyped) proved to be wild-type at C563, so it is unlikely that our C563-focused approach missed a significant number of other molecular causes for G6PD deficiency. Nevertheless, a more systematic genotyping strategy might have revealed rare cases of G6PD deficiency with an alternative molecular cause in Afghanistan.

Population Genetics

A number of polymorphic markers for G6PD have been described and used to explore the evolution of the gene in various contexts, and the abnormal functional variants thought to offer protection against malaria in particular. Several markers are only polymorphic in Africa (with mutation producing the A and A- forms) and only the markers 1311C>T and IVSXI C93T are of relevance in Asia. We noted that in the published literature, this second marker is generally interpreted as a T93C mutation, the reason being a historical one in that the first studies of this polymorphism were undertaken in Asia, where the T allele predominates in this position. However a global analysis of haplotypes indicates that both polymorphisms are ancient mutations that almost certainly predate human migration out of Africa, and the ancestral allele is clearly a C at both positions; it is also labelled as such at the NCBI SNP database. In Eurasian populations the ancestral CC genotype has become rare, and the single mutant IVSXI C93T haplotype (C1311, IVSXI 93T, referred to as CT) predominates, presumably as a result of genetic drift associated with population bottlenecks during human migrations into Europe and Asia. The TC haplotype (mutation at 1311C>T only) is rarely found in the general population in the Mediterranean and the Far East but is found at intermediate levels in Africa, India and (as shown here) Afghanistan (Figure 2). It is unsurprising that the double mutant TT haplotype has never been recorded, since this would have required recombination between two markers separated by approximately 150 DNA base pairs.

Understanding this haplotype background is a useful basis for exploring the history of the Mediterranean variant in Europe and Asia. The 563C>T mutation is likely to have been relatively recently introduced into southern Europe/Middle East [34] where it is strongly associated with the 1311C>T mutation (as far east as Iran and Iraq). In contrast, our study found a strong predominance of the wild-type allele at position C1311 in association (‘hitchhiking’) with the Mediterranean mutation in Afghanistan, findings that are consistent with previous data from Pakistan [9], [10] and India [22], [31]. The data therefore provide further support for the idea of distinct origins for the 563C>T mutation in southern Europe and the Middle East compared to Afghanistan and Pakistan [22], with little admixture between these populations in the subsequent period (despite the relatively small intervening distance).

It is not possible with the data obtained to comment in robust way on the driving forces and timescale for the emergence of the 563C>T Mediterranean allele in Afghanistan; this would require examination of markers over a wider length of the X-chromosome to determine the extent to which homozygosity extends beyond the polymorphic markers already examined [35], and/or a study of nearby microsatellite repeats [34]. The data do show that in individuals in Afghanistan without the Mediterranean mutation, background polymorphic markers are evenly mixed among various ethnic groups, a finding that is consistent with the late development of population structure between Afghanistan’s various ethnic groups after the Bronze age [36]. The relatively high prevalence of 563C>T reached in Pashtuns (8.4%) and the clear differences in its prevalence between different ethnic groups, also described in another large survey [32], are consistent with positive selection occurring after that period, in an analogous manner to the emergence of both the 563C>T variant in Europe and the Middle East [34] and the G6PD Mahidol variant in southeast Asia [37]. The detailed mechanisms by which G6PD deficiency might provide a fitness advantage remains unclear, but whereas P. falciparum appears to be a driving force for selection of the A- variant in Africa [38], recent studies in Asian populations indicate that protection against P. vivax infection may be more relevant [11], [37].

Conclusions

This genotypic study of the Mediterranean variant of G6PD deficiency in males in nine provinces in Afghanistan confirms that the 563C>T mutation reaches its highest frequencies in Pashtuns and Pashais (8.9%), with lower prevalence in groups historically associated with northern provinces and countries to the north of Afghanistan. The findings illustrate the need to assess G6PD status in all individuals (irrespective of stated ethnicity) prior to administration of primaquine for radical treatment of P. vivax infection, and also shed further light on the evolutionary history of G6PD variants in humans.