GST M1-T1 null Allele Frequency Patterns in Geographically Assorted Human Populations: A Phylogenetic Approach

Abstract

Genetic diversity in drug metabolism and disposition is mainly considered as the outcome of the inter-individual genetic variation in polymorphism of drug-xenobiotic metabolizing enzyme (XME). Among the XMEs, glutathione-S-transferases (GST) gene loci are an important candidate for the investigation of diversity in allele frequency, as the deletion mutations in GST M1 and T1 genotypes are associated with various cancers and genetic disorders of all major Population Affiliations (PAs). Therefore, the present population based phylogenetic study was focused to uncover the frequency distribution pattern in GST M1 and T1 null genotypes among 45 Geographically Assorted Human Populations (GAHPs). The frequency distribution pattern for GST M1 and T1 null alleles have been detected in this study using the data derived from literatures representing 44 populations affiliated to Africa, Asia, Europe, South America and the genome of PA from Gujarat, a region in western India. Allele frequency counting for Gujarat PA and scattered plot analysis for geographical distribution among the PAs were performed in SPSS-21. The GST M1 and GST T1 null allele frequencies patterns of the PAs were computed in Seqboot, Gendist program of Phylip software package (3.69 versions) and Unweighted Pair Group method with Arithmetic Mean in Mega-6 software. Allele frequencies from South African Xhosa tribe, East African Zimbabwe, East African Ethiopia, North African Egypt, Caucasian, South Asian Afghanistan and South Indian Andhra Pradesh have been identified as the probable seven patterns among the 45 GAHPs investigated in this study for GST M1-T1 null genotypes. The patternized null allele frequencies demonstrated in this study for the first time addresses the missing link in GST M1-T1 null allele frequencies among GAHPs.

About the Authors:

Senthilkumar Pitchalu Kasthurinaidu

† Deceased.

Affiliations Department of Biotechnology, Shree M. & N. Virani Science College, Rajkot, Gujarat, 360 005, India, Department of Plant Biotechnology, School of Biotechnology, Madurai Kamaraj University, Madurai, Tamilnadu, 625 021, India

Thirumurugan Ramasamy

* E-mail: [email protected] (TR)

Affiliation: Department of Animal Science, School of Life Sciences, Bharathidasan University, Tiruchirappalli, Tamilnadu, 620 024, India

Jayachitra Ayyavoo

Affiliation: Department of Plant Biotechnology, School of Biotechnology, Madurai Kamaraj University, Madurai, Tamilnadu, 625 021, India

Dhvani Kirtikumar Dave

Contributed equally to this work with: Dhvani Kirtikumar Dave, Divya Anantray Adroja

Affiliation: Department of Biotechnology, Shree M. & N. Virani Science College, Rajkot, Gujarat, 360 005, India

Divya Anantray Adroja

Contributed equally to this work with: Dhvani Kirtikumar Dave, Divya Anantray Adroja

Affiliation: Department of Biotechnology, Shree M. & N. Virani Science College, Rajkot, Gujarat, 360 005, India

Introduction

Metabolic activities play an important role in shaping the livelihood of an organism. Drug-xenobiotic compound metabolizing enzyme (XME) systems are the most investigated pathways that are involved in equilibrating the health status of an individual. Among the numerous drug related genes investigated, Glutathione-S-transferases (GST) of Phase II XMEs were found to play an important role in cellular protection and in cellular resistance to drugs by glutathione conjugation reactions. GST classes convert the active endogenous and/or exogenous carcinogenic compounds to their detoxified form. Among the GST classes, GST M1 and GST T1 were found associated to a loss of function with a structural deletion (null mutation); moreover, they were also found to modify the detoxification ability of the individual exposed to tobacco or carcinogenic pollutants in the environment [1]. Genotoxins such as aromatic hydrocarbon epoxides and products of oxidative stress such as DNA hydroperoxides, polycyclic aromatic hydrocarbon diol epoxides are catalyzed and detoxified by GST M1 while, the constituents of cigarette smoke such as alkyl halides, bezo (a) pyrene diol epoxide, acrolein are catalyzed and detoxified by GST T1 [1, 2]. Several factors such as environmental pollution, dietary habits and activity-dependent genetic differences have been reported as modulators of GST expression and susceptibility to xenobiotic compound detoxification [3]. Numerous studies in the recent past have hypothesized the difference in metabolic rate of M1 and T1 classes of GST as the risk factor associated to cancers of bladder, pancreas, upper aero digestive tract, lung, esophageal, head-neck, melanoma and also in Balken endemic nephropathy patients [4–9]. Further, the inter-individual difference in drug disposition and efficacy has been investigated by various authors [10] and the observed frequency distributions of GST M1-T1 genotypes among different populations are reported as ethnic or PAs dependent [10, 11]. Drugs are the major hope of remedy for the people around the globe with various metabolic and genetic disorders but the scenario in the past was found distressed as the effectiveness of the drugs were reported by the influence of the unidentified polymorphic patterns observed in drug metabolism genes among different ethnics or PAs [11–13]. Though researchers from different PAs are interested in analyzing the frequencies of GST M1 and T1 null genotypes and their possible risk association with various disorders, they are not able to report conclusive association in all major PAs [14, 15]. Recent advances in molecular techniques have opened a new era of pharmacogenomics and several researchers are inclined towards investigating the relationship in genetic diversity and allelic frequency of GST classes to insight genetic predisposition or susceptibility among various ethnics or PAs. In this context, probing the genetic variation in GST classes is inevitable for genomic epidemiological studies and to develop new common drugs in future to majority of PAs [9, 16, 17]. The allele frequency pattern in GST M1-T1 null genotypes of different PAs are yet to be explored to unlock several phenomenons related to a risk association with genetic diseases and drug dispositions [10]. A study including statistically valid number of subjects from various major PAs could address the issue of understanding the phenomenon for frequency distribution pattern in geographically assorted human populations (GAHPs); however, it will be tedious and might require huge population size [16]. Therefore, the present study was focused to uncover the genetic distance based ancestral origin or genetic affinity among GAHPs to address the paradigm for GST M1-T1 null allele frequency diversity. We are currently exploring how best to do this for the large number of populations in the present analysis to understand this phenomenon of frequency distribution pattern in GAHPs. GST M1 and GST T1 loci investigated in this study have been derived from literatures representing 44 different populations affiliated from Africa, Asia, Europe, South America and the genome of Gujarat PA, a region in western India. GST M1-T1 null allele frequency of 45 GAHPs were computed for phylogenesis with pair wise genetic distance based unweighted pair group method with arithmetic mean (UPGMA) and the findings of seven patterns for GST M1-T1 null allele frequency in this study have been demonstrated for the first time with highest genetic affinity. The patterns of null allele frequencies reported in this study add insights to determine a conclusive risk association of GST M1-T1 loci with several cancer or genetic disorders.

Materials and Methods

Subjects

The present investigation includes GST M1 and GST T1 null allele frequency of 45 GAHPs from 39 studies. Null allele frequency of Gujarat population was investigated in this study from 504 healthy unrelated volunteers of Gujarati origin with a mean age of 60 years. After signing the informed consent to participate in the study, blood samples of 2 mL were collected from each subject. Data of the remaining 44 populations were collected from different populations investigated by various authors (Table 1). Several studies of same ethnicity were also gathered in the study to fulfill the statistical significance of the study and to minimize the varying frequency of polymorphism among the ethnic groups while, the data gathered from Naveen et al. [18] had allele frequencies of combined Tamilnadu and Pondicherry PAs. The study was approved by the institutional ethics committee of Shrimathi Vasantben Ratilal Desai Cancer Research Centre, Rajkot Cancer Society - India.

[Figure omitted. See PDF.]

Table 1. GST M1-T1 null allele frequency of geographically assorted human populations.

https://doi.org/10.1371/journal.pone.0118660.t001

DNA isolation and Genotyping

Lahiri and Nurnberger method was used to isolate genomic DNA from whole blood [19]; the Huang et al., method of multiplex polymerase chain reaction was performed to identify GST M1 and T1 polymorphism with albumin gene as internal control [20]. Amplified products of PCR were visualized in 2% agarose gel and the band patterns were analyzed for polymorphism.

Geographically assorted human populations (GAHPs)

On the basis of interest in allele frequency patterns and availability of GST M1-T1 null allele frequency data shared by all populations, we choose 45 representative geographically assorted human populations around the world from 38 investigations reported by various authors from Asia [10, 14, 15, 18, 21–38], Europe [4–7, 9, 10, 39–42], Africa [21, 37, 43–50] and South America [49, 50] as summarized in Table 1. Of these 45 GAHPs, 4 were chosen from Eastern Asia (Japan, Korea, China and Mongolia [10, 21]), 4 from South Eastern Asia (Vietnam, Philippines, Indonesia and Singapore-Malay [10, 22–25]), 8 from Southern Asia India (Tamilnadu, Kerala, Karnataka, Andhra Pradesh, Maharashtra, West Bengal, Uttar Pradesh and Gujarat [14, 15, 18, 26–35, present study]), 3 from Southern Asia (Afghanistan, Iran, Pakistan [36–38]), 4 from Northern Europe (Sweden, Finland, Denmark and UK [5, 10, 39]), 4 from Southern Europe (Italy, Spain, Slovenia and Greece [5–7, 10]), 5 from Eastern Europe (Czech Republic, Bulgaria, Poland, Slovakia and Russia [9, 10, 39–42]), 3 from Western Europe (Netherlands, Germany and France [4, 5, 10]), 8 from Africa (Egypt, Nigeria, Xhosa tribe, Namibia, Cameroon, Ethiopia, Somalia and Zimbabwe [21, 37, 43–50]), 1 from South American Brazil [49, 50] and Caucasian (Americans and Canadians [10]). The “Caucasian” population used in this study was arbitrarily termed as “West Asian Caucasians” (wAs_Cau) to precise the geographical region and allele frequency patterns. All the PAs were grouped into continental regions as per the guidelines of Statistics Division of the United Nations (http://unstats.un.org/unsd/methods/m49/m49regin.htm [51]). Initially, 20 different continental region populations, as summarized in Table 2, were used in phylogenetic analysis for GST M1-T1 null allele frequency to minimize the effect of inbreeding [52]. Of these 20 continental region populations, the null allele frequency of thirteen populations such as Afghanistan, Iran, Pakistan (3 from South Asia), Brazil (1 from South America), Egypt (1 from North Africa), Nigeria (1 from West Africa), Xhosa tribe, Namibia (2 from South Africa), Cameroon (1 from Middle Africa), Ethiopia, Somalia and Zimbabwe (3 from East Africa), Caucasian (West Asia, as described earlier) and the average null allele frequency of seven continental regions such as East Asia (Japan, Korea, China, Mongolia), South East Asia (Vietnam, Philippines, Indonesia, Singapore-Malay), South Asia India (Tamilnadu, Kerala, Karnataka, Andhra Pradesh, Maharashtra, West Bengal, Uttar Pradesh, Gujarat), North Europe (Sweden, Finland, Denmark, UK), West Europe (Netherlands, Germany, France), South Europe (Italy, Spain, Slovenia, Greece) and East Europe (Czech Republic, Bulgaria, Poland, Slovakia, Russia) were used in the phylogenetic analyses, as summarized in Table 1 for GST M1-T1 null allele frequency. Since, some of the populations have considerable gene admixture from other PAs [52] another phylogenetic analysis was performed using all the 45 GAHPs for GST M1-T1 null allele frequency (Table 3).

[Figure omitted. See PDF.]

Table 2. GST M1-T1 null allele frequency based genetic distance between 20 different continental region populations.

https://doi.org/10.1371/journal.pone.0118660.t002

[Figure omitted. See PDF.]

Fig 1. Phylogenetic tree of 20 different continental regions population affiliations for GST M1-T1 null allele frequency.

The tree was produced by the UPGMA method from DST values in Table 2 and cluster with more than 50% of 1000 bootstrap replicates were included in the consensus tree obtained by Felsenstein (1989) phylogeny interference package. Major group of GST M1-T1 null allele frequencies were from population of Xhosa tribe, Zimbabwe, Ethiopia, Egypt, Afghanistan and Caucasian. Abbreviations used were same as those in Table 1.

https://doi.org/10.1371/journal.pone.0118660.g001

[Figure omitted. See PDF.]

Table 3. GST M1-T1 null allele frequency based genetic distance between 45 geographically assorted human populations (GAHPs).

https://doi.org/10.1371/journal.pone.0118660.t003

[Figure omitted. See PDF.]

Fig 2. Phylogenetic tree of 45 geographically assorted human populations for GST M1-T1 null allele frequency.

This tree was based on DST values in Table 3. Other aspects were same as those in Fig. 1. Major group of GST M1-T1 null allele frequencies were from population of Xhosa tribe, Zimbabwe, Ethiopia, Egypt, Afghanistan, Caucasian and Andhra Pradesh. Abbreviations used were same as those in Table 1.

https://doi.org/10.1371/journal.pone.0118660.g002

Statistical analysis

Distributions of GST M1 and GST T1 null alleles in Gujarati population were calculated by frequency counting method in SPSS-21 (4-27AEA) for windows. The standard genetic distance (DST) between different PAs for GST M1-T1 null allele frequencies were calculated by Nei’s (1972) method in Phylip 3.69 version [52, 53]. Least DST values between the PAs were used to compute clades with more than 50% of 1000 bootstrap replicates by Felsenstein (1989) method and then the phylogenetic trees were constructed in Mega-6 software by UPGMA method [54–56]. Finally, the clusters of PAs split found among the geographically assorted human populations in phylogenetic tree were used in the scattered plot to analyze their geographical distribution. The longitude (X-axis) and latitude (Ys-axis) of different continental regions were used to construct the scattered plot in SPSS-21 as summarized in Table 4. The online web source world atlas was used to compute the latitude and longitude of the respective geographic locations (http://www.worldatlas.com/aatlas/latitude_and_longitude_finder.htm [57]).

[Figure omitted. See PDF.]

Table 4. Geographical location of population affiliations from 20 different continental regions used in scattered plot analyses.

https://doi.org/10.1371/journal.pone.0118660.t004

Results

Phylogenetic tree for GST M1-T1 null allele frequency in GAHPs

The frequency of GST M1 and T1 null genotypes in Gujarat populations of India was observed as 0.200 and 0.355 respectively in this present study. The pair wise genetic distance matrix computed for GST M1-T1 null allele frequency between the respective PAs using Nei’s (1972) standard genetic distance method (DST) in gendist program of Phylip-3.69 software for the 20 different continental regions PAs and the 45 GAHPs described elsewhere were summarized in Tables 2 and 3 respectively. The phylogenetic analyses of 20 different continental regions (Fig. 1) and 45 GAHPs (Fig. 2) for GST M1-T1 null allele frequency were performed using the pair wise genetic distance matrix by UPGMA method in MEGA-6 software [52, 56]. Consense program that clustered more than 50% of 1000 bootstrap replicates was used to assess the reliability of the constructed phylogenetic trees [53, 54]. The Nei’s DST value varies from 0.0001 to 0.007 (Table 2) and 0.0006 to 0.008 (Table 3) for GST M1-T1 null allele frequency among the 20 different continental regions PAs and the 45 GAHPs respectively. The phylogenetic trees in Figs. 1 and 2 showed consistent clusters or patterns for GST M1-T1 null allele frequency from populations such as South African Xhosa tribe, South African Zimbabwe, East African Ethiopia, North African Egypt, Caucasian and South Asian Afghanistan. The least DST value of PAs from South Asian Iran (0.000882) and East African Somalia (0.00552) with East African Ethiopian allele frequency grouped them as one of the clusters in the phylogenetic trees for GST M1-T1 null allele frequency. Few of the African PAs such as West African Nigeria, Middle African Cameroon and South African Namibia were clustered with South African Xhosa tribe allele frequency with least DST value that ranged from 0.006642 to 0.007702. In addition, the South Asian Indian Gujarat PAs (0.00214) was observed with highest affinity to South African Xhosa cluster than any other clusters (Fig. 2). The average GST M1-T1 null allele frequency in South Asian Indian PA was observed as a cluster to East African Zimbabwe allele frequency (Fig. 1) with least DST value (0.007203) out of 20 continental regions. Further, the least DST value to East African Zimbabwe allele frequency was observed with South Asian Indian PAs (Table 3) from Tamilnadu (0.0035) and Karnataka (0.003729). Furthermore, PA from Karnataka was observed with least DST value to Maharashtra allele frequency (0.003307), which was clustered with allele frequency from West Bengal (0.000574), Uttar Pradesh (0.003826) and Kerala PAs of South Asian India (0.006243) as shown in Fig. 2. Nevertheless, the South Asian Indian Andhra Pradesh allele frequency was observed with least DST value only to Uttar Pradesh PA (0.005068) than any other Indian PAs for GST M1-T1 null allele frequency and stood as a separate cluster in the phylogenetic tree with highest genetic affinity to South America (Brazil, 0.001318) and South East Asia allele frequency (Vietnam, 0.004738) from geographically distant continental regions.

PAs from China (0.000757), Japan (0.001457) and Korea (0.001522) of East Asia and Pakistan (0.001511) of South Asia were observed as another cluster with least DST value to Caucasians [Americans and Canadians (10)] for GST M1-T1 null allele frequency out of 45 combinations (Fig. 2). However, Mongolia of East Asia was observed with least DST value to Pakistan allele frequency (0.000686) than the Caucasians (0.046671). The least DST value between European continental regions and South Asian Afghanistan allele frequency that ranged from 0.002064 to 0.004708 was clustered together for GST M1-T1 null allele frequency in the phylogenetic tree of 20 different continental regions (Fig. 1). Nevertheless, the phylogenetic tree of 45 GAHPs (Fig. 2) clustered only 13 European PAs (Sweden, Finland, Denmark, Netherlands, Germany, France, Italy, Spain, Greece, Bulgaria, Poland, Slovakia and Russia) out of 16 investigated in this study with South Asian Afghanistan allele frequency (least DST value that ranged from 0.001176 to 0.00723) while, the other 3 PAs [Slovenia (0.002462), Czech Republic (0.004116) and UK (0.006718)] were clustered with North African Egypt allele frequency. Singapore-Malay and Indonesia PAs from South East Asia were observed with least DST value to East Asia (China, 0.001246) and Caucasian (0.003788) respectively for GST M1-T1 null allele frequency. Nevertheless, the other counter parts from same continental region were observed as the most diverse PAs with PA admixture from North Africa (Egypt, 0.002571) and East Africa (Ethiopia, 0.007943) for Philippines; South Asia India (Andhra Pradesh, 0.004738) and East Africa (Ethiopia, 0.004922) for Vietnam among the 45 GAHPs investigated in this study for GST M1-T1 null allele frequency as shown in Table 3 and Fig. 2 respectively.

GST M1-T1 null allele frequency patterns among the GAHPs

The effect of isolation by geographical distance in population differentiation [51] was validated in a scattered plot with respect to the phylogenetic clusters of 45 GAHPs for GST M1-T1 null allele frequency that corresponds to the latitudes and longitudes of 20 different continental regions representing PAs from Africa, Asia, Europe and America (Table 4). The scattered plot illustrated in Fig. 3 suggest three major geographical split for the seven GST M1-T1 null allele frequency clusters or patterns observed in the phylogenetic tree of 45 GAHPs (Fig. 2). South African Xhosa allele frequency pattern (I) observed mostly in continental regions of Africa suggest an “Africa” split in the scattered plot with least population differentiation to Nigeria (West Africa), Cameroon (Middle Africa) and Namibia (South Africa). However, the GST M1-T1 null allele frequency patterns of PAs from other African continental region such as East African Zimbabwe (II), East African Ethiopia (III) and North African Egypt (IV) were observed with least population differentiation to PAs from non African continental regions such as South Asia (India and Iran), South East Asia (Philippines and Vietnam), Southern Europe (Slovenia), Eastern Europe (Czech Republic) and Northern Europe (UK) irrespective of the geographical isolation suggest an “out of Africa” split. Finally, the remaining three non-African GST M1-T1 null allele frequency patterns observed from Caucasian (V), South Asian Afghanistan (VI) and South Asian Indian Andhra Pradesh (VII) were geographically distributed in different continental regions such as Asia, Europe and America with the exception of Africa suggest an “other than Africa” split among the 45 GAHPs in the scattered plot.

[Figure omitted. See PDF.]

Fig 3. Geographic Location of GST M1-T1 null allele frequency patterns in Scattered Plot.

This scenario is largely a speculation based on Nei’s DST based phylogenetic trees in Figs. 1 and 2, which suggest seven patterns for GST M1-T1 null allele frequency among the 45 geographically assorted human populations. The scattered plot shows the genetic affinity and geographical distribution of the seven probable patterns with three major split among the continental regions such as an “African” split with South African Xhosa pattern (I), an “out of Africa” split with East African Zimbabwe (II), East African Ethiopia (III), North African Egypt (IV) patterns and an “other than Africa” split with South Asian Afghanistan (V), West Asian Caucasian (VI), South Asian Indian Andhra Pradesh patterns (VII). Abbreviations used were same as those in Table 1.

https://doi.org/10.1371/journal.pone.0118660.g003

Discussion

Understanding the genetic variation and diversity among the geographically assorted human populations (GAHPs) is an interesting topic in population genetics with a wide range of neutral genetic markers and adaptive markers being employed to uncover the patterns of genetic diversity [58]. To address the extent of diversity in allele frequency distribution among populations from different ethnicity, region, country or continent is difficult, nevertheless understanding this phenomenon is inevitable [8, 11, 16, 17] and recent advances in the molecular techniques excel the perspective of inter-individual genetic variations in GAHPs. Allele frequencies of large number of neutral markers or of even few candidate markers that duplicate or decay to favor new environments and lead to rapid adaptations are often used for investigating the patterns [58]. The paradigm of allele frequency among the populations holds the key to unlock the existing problem of inter-individual genetic variation in xenobiotic metabolizing enzymes (XMEs) and in particular the decay or null allele frequency of Glutathione-S-transferase’s classes such as Mu 1 (GST M1) and Theta 1 (GST T1), which are considered as the major risk factor for various diseases including several types of cancers [8–10]. Therefore, the present investigation analyzed the pattern for GST M1-T1 null allele frequency among GAHPs using a phylogenetic approach. A set of 20 different continental region PAs (Table 2) and 45 GAHPs (Table 3) were recruited for GST M1-T1 null allele frequency data from 38 previously reported works and genomic data of Indian Gujarat PA in this study (Table 1) and the respective phylogenetic trees (Figs. 1 and 2) have been constructed by UPGMA method based on Nei’s (1972) standard genetic distance (DST) with clusters more than 50% of 1000 bootstrap replicates obtained by Felsenstein (1989) program [52, 54, 56]. In addition to the ancestral origin or genetic affinity based clusters for GST M1-T1 null allele frequency demonstrated in the phylogenetic trees, a positive correlation between the genetic distance and geographical distance were analyzed for the effect of isolation in population differentiation by distance in a scattered plot (Fig. 3). Indeed, the observations from phylogenetic trees and scattered plot of different PAs constructively reveals the findings of seven probable patterns for GST M1-T1 null allele frequency among the GAHPs in concordance to the reports of archeological signatures, ancient gene flows and sex-specific components [59–62].

The genetic affinity and geographical distribution of 20 different continental regions that included 45 GAHPs investigated in this study (Figs. 1–3) revealed the findings of an allele frequency pattern for GST M1-T1 null genotypes among Namibia (South Africa), Nigeria (West Africa), Cameroon (Middle Africa), Gujarat (South Asian Indian) and Xhosa tribe (South Africa) for the first time. We report here, the findings of Xhosa allele frequency (I) with major genetic affinity towards populations from Africa (Namibia, Nigeria, Cameroon) as an “Africa” split pattern for GST M1-T1 null allele frequency in agreement to the reports of linkage disequilibrium computed for loss of variants in GST classes by Polimanti et al. (2013). The observations in phylogenetic trees (Figs. 1 and 2) and scattered plot analysis (Fig. 3) demonstrated the findings of another three patterns such as East Africa Zimbabwean allele frequency - II in population from India (South Asia), East Africa Ethiopian allele frequency - III in populations from Iran (South Asia) and Somalia (East Africa) and North Africa Egyptian allele frequency - IV in populations from Slovenia (Southern Europe), Czech Republic (Eastern Europe) and UK (Northern Europe) for GST M1-T1 null genotypes. The findings of Ethiopian - III and Egyptian - IV allele frequency pattern for GST M1-T1 null genotypes in this study are in concordance to the earlier reports of genome wide diversity study in the Levant by Haber et al., (2013), who found two major groups with one close to Africans and Middle Easterners and the other closer to modern day Europeans [61]. Further, the findings of these Zimbabwean, Ethiopian and Egyptian patterns from African populations with high genetic affinity towards non-African populations for GST M1-T1 null allele frequency have been reported as an “Out of Africa” split in this study in corroborate to the findings of Templeton (2002), who reported the out-of-Africa theory of migration and the ancestral root of allele frequency admixture [63].

The allele frequency of population from South Asian Afghanistan with high genetic affinity to majority of European PAs investigated in this study (Tables 2 and 3 and Figs. 2 and 3) has been reported as pattern - V for GST M1-T1 null genotypes in accordance to the reports of various authors [60, 62, 64]. Further, population in Pakistan (South Asia) has been reported with Afghanistan (South Asia) pattern for GST M1-T1 null allele frequency (Table 3) in this study, though it was found with genetic affinity to PAs from Mongolia (East Asia), Europe (South, East and West) and Andhra Pradesh (South India) in corroborate to the earlier reports of Templeton (2002), who stated the findings of considerable overlap among East Asians, Europeans and South Indian populations [64]. Moreover, the pattern of allele frequency from Caucasians (Americans and Canadians) found among East Asians (Fig. 1) in this study has been identified as pattern - VI for GST M1-T1 null allele frequency and reported for the first time. Finally, the allele frequency from South Indian Andhra Pradesh PA was found with least genetic distance (Table 3) to populations from Pakistan (South Asia), Vietnam (South East Asia) and Brazil (South America) irrespective of the phenomenon of population differentiation by geographical isolation [52] and has been reported as pattern - VII for GST M1-T1 null genotypes (Figs. 2 and 3) among the 45 GAHPs in this study. These observations of South East Asian and South American PAs with the null allele frequency pattern from South Indian Andhra Pradesh PA are in agreement to the reports of agro-pastoral system in South India that acted as agricultural center and source of dispersion to lineages from different preexisting populations [60]. Furthermore, the reported East Africa patterns from Zimbabwe (II), Ethiopia (III) among India, Iran (South Asia) and South Asian pattern from South Indian Andhra Pradesh (VII) among Vietnam (South East Asia) populations in this study for GST M1-T1 null allele frequency are in concordance to the reports of migration pattern of Homo sapiens from East Africa with the demographic expansions by Field and Lahr (2006), who investigated the geographic information systems during oxygen isotope stage 4 [62]. GST M1-T1 null allele frequency from South East Asian PAs has been reported as the complex admixture of Zimbabwe (II), Ethiopia (III) and Andhra Pradesh (VII) patterns in this study. Finally, the scattered plot analysis (Fig. 3), clearly demonstrates the findings of allele frequency patterns from South Asian Afghanistan - V, Caucasian - VI and South Indian Andhra Pradesh - VII as an “Other than Africa” split among 45 GAHPs for GST M1-T1 null genotypes with respect to their geographical distribution. This observation of other than Africa split in this study has been reported here in agreement to the concepts of later migration of the populations in regions other than Africa [60, 64]. In conclusion, the data of seven patterns for GST M1-T1 null allele frequency from Xhosa tribe (I), Zimbabwe (II), Ethiopia (III), Egypt (IV), Afghanistan (V), Caucasian (VI) and South Indian Andhra Pradesh (VII) reported in this study compare constructively with the earlier studies that suggested the PAs of relatively recent origin show comparatively small genetic differences and high genetic affinity among them [11, 46, 52, 58]. Findings of these seven patterns (I-VII) for GST M1-T1 null allele frequency reported here, would shed some light to address the missing link in most of the genomic epidemiological studies that lacks conclusive risk association [9, 16, 17]. The “Africa” (I), “Out of Africa” (II, III and IV) and “Other than Africa” (V, VI and VII) split among the 45 GAHPs reported in this study have to be explored further to rationalize the GST M1-T1 null allele’s frequency patterns in world populations.

Acknowledgments

The authors wish to thank the participants of this study for their cooperation; Dr. V. K. Gupta, Director Shree Gulabkuvar Talakchand Sheth Cancer hospital, Rajkot for providing technical guidance; Shree Manibhai Virani and Shremathi Navalben Virani Science College, Rajkot, Gujarat, India and Department of Animal Science, Bharathidasan University, Tiruchirappalli, Tamilnadu, India for providing the support in sample collection and facilities respectively. The Authors are grateful to Dr. Akbar Ali Khan Pathan, Assistant Professor, Genome Research Chair, Department of Biochemistry, College of Science, King Saud University, Riyadh-11451, Kingdom of Saudi Arabia for facilitating data analysis and for Dr. K Balakrishnan, Associate Professor, Department of Immunology, Madurai Kamaraj University, Madurai, Tamilnadu - 625 021, India for his critical suggestions.

Author Contributions

Conceived and designed the experiments: SPK DKD DAA JA TR. Performed the experiments: SPK DKD DAA. Analyzed the data: SPK JA TR. Contributed reagents/materials/analysis tools: SPK DKD DAA. Wrote the paper: SPK JA TR.

Citation: Kasthurinaidu SP, Ramasamy T, Ayyavoo J, Dave DK, Adroja DA (2015) GST M1-T1 null Allele Frequency Patterns in Geographically Assorted Human Populations: A Phylogenetic Approach. PLoS ONE 10(4): e0118660. https://doi.org/10.1371/journal.pone.0118660

References

1. Strange RC, Spiter MA,Ramachandran S,Fryer AA. Glutathione-S-transferase family of enzymes. Mutat Res. 2001;482:21–26. pmid:11535245

2. Kim WJ, Kim H,Kim CH,Lee MS,Oh BR,Lee HM, et al. GSTT1-null genotype is a protective factor against bladder cancer. Urology 2002;60:913–918. pmid:12429337

3. Thier R, Bruning T, Roos PH, Rihs H, Golka K, Ko Y, et al. Markers of genetic susceptibility in human environmental hygiene and toxicology: the role of selected CYP, NAT and GST genes. Int J Hyg Environ Health 2003;206:149–171. pmid:12872524

4. Abbas A, Delvinquiere K, Lechevrel M, Lebailly P, Gauduchon P, Launoy G, et al. GSTM1, GSTT1, GSTP1 and Cyp1A1 genetic polymorphisms and susceptibility to esophageal cancer in a French population: different pattern of squamous cell carcinoma and adenocarcinoma. World J Gastroenterol 2004;10:3389–3393. pmid:15526353

5. Bu H, Rosdahl I, Holmdahl-Kallen K, Sun XF, Zhang. Significance of glutathione S-transferase M1, T1 and P1 polymorphisms in Swedish melanoma patients. Oncol Rep. 2007;17:859–864. pmid:17342328

6. Dialyna IA, Miyakis S, Georgatou N, Spandidos DA. Genetic polymorphisms of CYP1A1, GSTM1 and GSTT1 genes and lung cancer risk. Oncol Rep. 2003;10:1829–1835. pmid:14534704

7. Garcia-Closas M, Malats N, Silverman D, Dosemeci M,Kogevinas M,Hein DW, et al. NAT2 slow Acetylation, GSTM1 null genotype and risk of bladder cancer: results from the Spanish Bladder Cancer Study and meta-analyses. Lancet 2005;366:649–659. pmid:16112301

8. Senthilkumar KP, Thirumurugan R. Risk Modulation of GSTM1 - GSTT1 interactions to Head and Neck Cancer in tobacco users. Mol Bio Rep. 2014;41:5635–5644. pmid:25015263

9. Toncheva DL, Von AN, Atanasova SY, Dimitrow TG, Armstrong VW, Oellerich M. Identification of NQ01 and GSTs genotype frequencies in Bulgarian patients with Balken endemic nephropathy. J Nephrol. 2004;17:384–389. pmid:15365958

10. Kurose K, Sugiyama E, Saito Y. Population Differences in Major Functional Polymorphisms of Pharmacokinetics / pharmacodynamics-related Genes in Eastern Asians and Europeans: Implications in the Clinical Trials for Novel Drug Development. Drug Metab Pharmacokinetics 2012;27:9–54. pmid:22123129

11. Polimanti R, Carboni C, Baesso I, Piacentini S, Iorio A, Stefano DGF, et al. Genetic variability of glutathione S-transferase enzymes in human populations: Functional inter-ethnic differences in detoxification systems. Gene 2013;512:102–107. pmid:23043933

12. Evans WE, Johnson JA. Pharmacogenomics: The Inherited Basis for Interindividual Differences in Drug Response. Annu Rev Genomics Hum Genet. 2001;2:9–39. pmid:11701642

13. O'Donnell PH, Dolan ME. Cancer Pharmacoethnicity: Ethnic Differences in Susceptibility to the Effects of Chemotherapy. Clin Cancer Res. 2009;15:4806–4814. pmid:19622575

14. Ghosh T, Gupta S, Bajpai P, Agarwal D, Agarwal M, Gupta OP, et al. Association of CYP1A1, GSTM1, and GSTT1 gene polymorphism with risk of oral submucous fibrosis in a section of North Indian population. Mol Biol Rep. 2012;39:9383–9389. pmid:22752755

15. Shukla D, Kale AD, Hallikerimath S, Vivekanandhan S, Venkatakanthaiah Y. Genetic polymorphism of drug metabolizing enzymes (GSTM1 and CYP1A1) as risk factors for oral premalignant lesions and oral cancer. Biomed Pap Med Fac Univ Palacky Olomouc Czech Repub. 2012;156:253–259. pmid:22660220

16. Saijo N. The Role of Pharmacoethnicity in the Development of Cytotoxic and Molecular Targeted Drugs in Oncology. Yonsei Med J. 2013;54:1–14. pmid:23225792

17. Salah GB, Kallabi F, Maatoug S, Mkaouar-Rebai E, Fourati A, Fakhfakh F, et al. Polymorphisms of glutathione S-transferases M1, T1, P1 and A1 genes in the Tunisian population: An intra and interethnic comparative approach. Gene 2012;498:317–322. pmid:22313526

18. Naveen AT, Adithan C, Padmaja N, Shashindran CH, Abraham BK, Sathyanarayanamoorthy K, et al. Glutathione S-transferase M1 and T1 null genotype distribution in South Indians. Eur J Clin Pharmacol. 2004;60:403–406. pmid:15197518

19. Lahiri DK, Nurnberger JI. A rapid non-enzymatic method for the preparation of HMW DNA from blood for RFLP studies. Nucleic Acids Res. 1991;19:5444. pmid:1681511

20. Huang K, Sandler RS, Millikan RC, Schroeder JC, North KE, Hu J. GSTM1 and GSTT1 polymorphisms, cigarette smoking, and risk of colon cancer: a population-based case-control study in North Carolina (United States). Cancer causes & control 2006;17:385–394.

21. Hamdy SI, Hiratsuka M, Narahara K, Endo N, El-Enany M, Moursi N, et al. Genotype and allele frequencies of TPMT, NAT2, GST, SULT1A1 and MDR-1 in the Egyptian population. Br J Clin Pharmacol. 2003;55:560–569. pmid:12814450

22. Agusa T, Iwata H, Fujihara J, Kunito T, Takeshita H, Minh TB, et al. Genetic polymorphisms in glutathione S-transferase (GST) superfamily and arsenic metabolism in residents of the Red River Delta, Vietnam. Toxicol Appl Pharmacol. 2010;242:353–362.

23. Rimando MG, Chu MN, Yuson E, de Castro-Berna G, Okamoto T. Prevalence of GSTT1, GSTM1 and NQ01(609C>T) in Filipino children with ALL (acute lymphoblastic leukemia). Biosci Rep. 2008;28:117–124. pmid:18444911

24. Amtha R, Ching CS, Zang R, Razak IA, Basuki B, Roeslan BO, et al. GSTM1, GSTT1 and CYP1A1 polymorphisms and risk of oral cancer: a case-contro study in Jakarta, Indonesia. Asian Pac J Cancer Prev. 2009;10:21–26. pmid:19469619

25. Lee EJ, Wong JY, Yoeh PN, Gong NH. Glutathione S transferase-theta (GSTT1) genetic polymorphism among Chinese, Malaya and Indians in Singapore. Pharmacogenetics 1995;5:332–334. pmid:8563775

26. Anantharaman D, Chaubal PM, Kannan S, Bhisey RA, Mahimkar MB. Susceptibility to oral cancer by genetic polymorphisms at CYP1A1, GSTM1 and GSTT1 loci among Indians: tobacco exposure as a risk modulator. Carcinogenesis 2007;28:1455–1462. pmid:17307803

27. Buch S, Kotekar A, Kawle D, Bhisey R (2001). Polymorphisms at CYP and GST gene loci. Prevalence in the Indian population. Eur J Clin Pharmacol 57 (6–7), 553–555. pmid:11699625

28. Buch SC, Notani PN, Bhisey RA. Polymorphism at GSTM1, GSTM3 and GSTT1 gene loci and susceptibility to oral cancer in an Indian population. Carcinogenesis 2002;23:803–807. pmid:12016153

29. Jain M, Kumar S, Rastogi N, Lal P, Ghoshal UC, Tiwari A, et al. GSTT1, GSTM1 and GSTP1 genetic polymorphisms and interaction with tobacco, alcohol and occupational exposure in esophageal cancer patients from North India. Cancer Lett. 2006;242:60–67. pmid:16338071

30. Konwar R, Manchanda PK, Chaudhary P, Nayak VL, Singh V, Bid HK. Glutathione S-transferase (GST) gene variants and risk of benign prostatic hyperplasia: a report in a North Indian population. Asian Pac J Cancer Prev. 2010;11:1067–1072. pmid:21133626

31. Mishra DK, Kumar A, Srivastava DS, Mittal RD. Allelic variation of GSTT1, GSTM1 and GSTP1 genes in North Indian population. Asian Pac J Cancer Prev. 2004;5:362–365. pmid:15546237

32. Singh M, Shah PP, Singh AP, Ruwali M, Mathur N, Pant MC, et al. Association of genetic polymorphisms in glutathione S-transferases and susceptibility to head and neck cancer. Mutat Res. 2008;638:184–194. pmid:18035380

33. Sabitha K, Reddy MVV, Jamil K. GST genotypes in head and neck cancer patients and its clinical implications. Afr J Biotechnol. 2008;7:3853–3859.

34. Sikdar N, Datta S, Dey B, Paul RR, Panda CK, Roy B. Homozygous Null Genotype at Glutathione S-transferase M1 Locus as a Risk Factor for Oral Squamous Cell Carcinoma in Indian Tobacco Users. Int J Hum Genet. 2005;5:37–44.

35. Sreelekha TT, Ramadas K, Pandey M, Thomas G, Nalinakumari KR, Pillai MR. Genetic polymorphism of CYP1A1, GSTM1 and GSTT1 genes in Indian oral cancer. Oral Oncol. 2001;37:593–598. pmid:11564581

36. Saify K, Saadat I, Saadat M. Genetic polymorphisms of glutathione S-transferase T1 (GSTT1) and M1 (GSTM1) in selected populations of Afghanistan. Mol Biol Rep. 2012;39:7855–7859. pmid:22539183

37. Saadat M. GSTM1 null genotype associated with age standardized cancer mortality rate in 45 countries from five continents: an ecologic study. Int J Cancer Res. 2007;3:74–91.

38. Shaikh RS, Amir M, Masood AI, Sohail A, Athar HU, Siraj S, et al. Frequency distribution of GSTM1 and GSTT1 null allele in Pakistani population and risk of disease incidence. Environ Toxicol Pharmacol. 2010;30:76–9 pmid:21787632

39. Garte S, Gaspari L, Alexandrie A, Ambrosone C, Autrup H, Autrup JL, et al. Metabolic Gene Polymorphism Frequencies in Control populations. Cancer Epidemiol Biomarkers Prev. 2001;10:1239–1248. pmid:11751440

40. Binkova B, Smerhovsky Z, Strejc P, Boubelik O Stavkova Z, Chvatalova I, et al. DNA-adducts and atherosclerosis: a study of accidental and sudden death males in the Czech Republic. Mutat Res. 2002;501:115–128. pmid:11934443

41. Gajecka M, Rydzaniacz M, Jaskula-Sztul R, Kujawski M, Szyfter W, Szyfter K. CYP1A1, CYP2D6, CYP2E1, NAT2, GSTM1 and GSTT1 polymorphisms or their combinations are associated with increased risk of the laryngeal squamous cell carcinoma. Mutat Res. 2005;574:112–123. pmid:15914211

42. Gra O, Mityaeva O, Berdichevets I, Kozhekbaeva Z, Fesenko D, Goldenkova-Pavlova I, et al. Microarray-Based Detection of CYP1A1, CYP2C9, CYP2C19, CYP2D6, GSTT1, GSTM1, MTHFR, MTRR, NQ01, NAT2, HLA-DQA1, and ABO allele frequencies in native Russians. Genet Test Mol Biomarkers 2010;14:329–342. pmid:20373852

43. Ebeshi BU, Bolaji OO, Masimirembwa CM. Glutothione-S-transferase (M1 and T1) polymorphisms in Nigerian populations. J Med Genet Genomics 2011;3:56–60.

44. Adams CH, Werely CJ, Victor TC, Hoal EG, Rossouw G, Heldin VPD. Allele frequencies for glutathione S-transferase and N-acetyltransferase 2 differ in African population groups and may be associated with oesophageal cancer or tuberculosis incidence. Clin Chem Lab Med. 2003;41:600–605. pmid:12747608

45. Fujihara J, Yasuda T, Iida R, Takatsuka H, Fujii Y, Takeshita H. Cytochrome P450 1A1, glutathione S-transferases M1 and T1 polymorphisms in Ovambos and Mongolians. Leg Med (Tokyo) 2009;Suppl 1:S408–S410. pmid:19264525

46. Piacentini S, Polimanti R, Porreca F, Martı´nez-Labarga C, Stefano GFD, Fuciarelli M. GSTT1 and GSTM1 gene polymorphisms in European and African populations. Mol Biol Rep. 2011;38:1225–1230. pmid:20563854

47. Buchard A, Sanchez JJ, Dalhoff K, Morling N. Multiplex PCR detection of GSTM1, GSTT1, and GSTP1 gene variants: simultaneously detecting GSTM1 and GSTT1 gene copy number and the allelic status of the GSTP1 Ile105Val genetic variant. J Mol Diagn. 2007;9:612–617. pmid:17916600

48. Dandara C, Sayi J, Masimirembwa CM, Magimba A, Kaaya S, Sommers DK, et al. Genetic polymorphism of cytochrome P450 1A1 (Cyp1A1) and glutathione transferases (M1, T1 and P1) among Africans. Clin Chem Lab Med. 2002;40:953–957.

49. Rossini A, Rapozo CM, Amorim MF, Macedo JM, Medina R, Neto JF, et al. Frequencies of GSTM1, GSTT1, and GSTP1 polymorphisms in a Brazilian population. Genet Mol Res. 2002;1:233–240. pmid:14963830

50. Magno LAV, Talbot J, Talbot T, Santos BAM, Souza RP, Marin LJ, et al. Glutathione-S-Transferase variants in a Brazilian Population. Pharmacology 2009;83:231–236. pmid:19258736

51. Composition of macro geographical (continental) regions, geographical sub-regions, and selected economic and other groupings. Available: http://unstats.un.org/unsd/methods/m49/m49regin.htm

52. Nei M. Genetic distance between populations. Am Nat. 1972;106:283–292.

53. Felsenstein J. Phylip - Phylogeny Inference Package (Version 3.6). Distributed by the author, Department of Genome Sciences, University of Washington, Seattle; 2004.

54. Felsenstein J. Phylip - Phylogeny Inference Package (Version 3.2). Cladistics 1989;5: 164–166.

55. Tamura K, Stecher G, Peterson D, Filipski A, and Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30:2725–2729. pmid:24132122

56. Sneath PHA and Sokal RR. Numerical Taxonomy. Freeman, San Francisco; 1973.

57. Latitude and Longitude Finder. Available: http://www.worldatlas.com/aatlas/latitude_and_longitude_finder.htm

58. Hoffmann AA, Willi Y. Detecting genetic responses to environmental change. Nat Rev Genet. 2008;9:421–432. pmid:18463665

59. Jorde L, Wooding S. Genetic variation, classification and’race'. Nat Genet. 2004;36:28–33.

60. ArunKumar G, Soria-Hernanz DF, Kavitha VJ, Arun VS, Syama A, Ashokan KS, et al. Population Differentiation of Southern Indian Male Lineages Correlates with Agricultural Expansions Predating the Caste System. PLoS ONE 2012;7: e50269. pmid:23209694

61. Haber M, Gauguier D, Youhanna S, Patterson N, Moorjani P, Botigue LR, et al. Genome-Wide Diversity in the Levant Reveals Recent Structuring by Culture. PLoS Genet. 2013;9: e1003316. pmid:23468648

62. Field J, Lahr MM. Assessment of the Southern Dispersal: GIS-Based Analyses of Potential Routes at Oxygen Isotopic Stage 4. J World Prehist. 2006;19:1–45.

63. Templeton A. Out of Africa again and again. Nature 2002;416:45–51. pmid:11882887

64. Bamshad M, Kivisild T, Watkins WS, Dixon ME, Ricker CE, Naidu JM, et al. Genetic Evidence on the Origins of Indian Caste Populations. Genome Research 2001;11:994–1104. pmid:11381027