SYSTEMATIC REVIEW article

Front. Genet., 07 November 2022

Sec. Genetics of Common and Rare Diseases

Volume 13 - 2022 | https://doi.org/10.3389/fgene.2022.959291

Individual and combined effects of the GSTM1, GSTT1, and GSTP1 polymorphisms on type 2 diabetes mellitus risk: A systematic review and meta-analysis

  • 1. Changzhi Medical College, Changzhi, Shanxi, China

  • 2. Department of Endocrinology, Heping Hospital Affiliated to Changzhi Medical College, Changzhi, Shanxi, China

  • 3. Department of Epidemiology, School of Public Health to Southern Medical University, Guangzhou, Guangdong, China

  • 4. Institute of Evidence-Based Medicine, Heping Hospital Affiliated to Changzhi Medical College, Changzhi, China

Article metrics

View details

10

Citations

3,6k

Views

1,1k

Downloads

Abstract

Backgrounds: Compared with previously published meta-analyses, this is the first study to investigate the combined effects of glutathione-S-transferase polymorphisms (GSTM1, GSTT1 and GSTP1 IIe105Val) and type 2 diabetes mellitus (T2DM) risk; moreover, the credibility of statistically significant associations was assessed; furthermore, many new original studies were published.

Objectives: To determine the relationship between GSTM1, GSTT1, and GSTP1 polymorphisms with T2DM risk.

Methods: PubMed, Embase, Wanfang, and China National Knowledge Infrastructure Databases were searched. We quantify the relationship using crude odds ratios and their 95% confidence intervals Moreover, the Venice criteria, false-positive report probability (FPRP), and Bayesian false discovery probability (BFDP) were used to validate the significance of the results.

Results: Overall, significantly increased T2DM risk was found between individual and combined effects of GSTM1, GSTT1, and GSTP1 polymorphisms on T2DM risk, but, combined effects of the GSTT1 and GSTP1 polymorphisms was not statistically significant. GSTT1 gene polymorphism significantly increases the risk of T2DM complications, while GSTM1 and GSTP1 polymorphisms had no statistical significance. The GSTM1 null genotype was linked to a particularly increased risk of T2DM in Caucasians; the GSTT1 null genotype was connected to a significantly higher risk of T2DM in Asians and Indians; and the GSTP1 IIe105Val polymorphism was related to a substantially increased T2DM risk in Indians. Moreover, the GSTM1 and GSTT1 double null genotype was associated with substantially increased T2DM risk in Caucasians and Indians; the combined effects of GSTM1 and GSTP1 polymorphisms was associated with higher T2DM risk in Caucasians. However, all significant results were false when the Venice criteria, FPRP, and BFDP test were used (any FPRP >0.2 and BFDP value >0.8).

Conclusion: The current analysis strongly suggests that the individual and combined effects of GSTM1, GSTT1 and GSTP1 polymorphisms might not be connected with elevated T2DM risk.

Introduction

Insulin resistance, excessive or improper glucagon secretion, and insufficient insulin secretion are symptoms of type 2 diabetes mellitus (T2DM), a chronic endocrine metabolic disorder. According to the newly released International Diabetes Federation Diabetes Atlas (10th edition), there are an estimated 537 million people living with diabetes in the world now, with that number expected to rise to 784 million by 2,045 based on current trends. In recent years, the incidence of diabetic complications, such as diabetic nephropathy and diabetic cardiovascular disease, has gradually increased. With the increase of people’s lifestyle such as poor diet, sedentary and inactive, the prevalence of T2DM will continue to rise, which will undoubtedly bring a huge burden on human health and medical expenditure (https://diabetesatlas.org/en/.). T2DM pathogenesis is complicated and diverse, involving both environmental and genetic risk factors. Currently, antioxidant and detoxification gene polymorphisms are believed to play a significant role in the risk of T2DM and related complications (Wu et al., 2014; Pahwa et al., 2017; Dendup et al., 2018).

Oxidative stress is the result of imbalance of reactive oxygen species generation and clearance, which is closely related to insulin resistance and functional impairment of islet β-cells. GST is a superfamily protease encoded by multiple genes with antioxidant and detoxification functions (Tangvarasittichai, 2015; Vats et al., 2015). It catalyzes the formation of thioether bonds between glutathione and electrophiles, thereby catalyzing the reduction of hydrogen peroxide, reducing fatty acids, phospholipids and DNA bases to corresponding alcohols, and scavenging lipid free radicals and hydrogen peroxide, thus reducing body damage (Hayes et al., 2005; Tabatabaei-Malazy et al., 2017). GSTM1 and GSTT1 gene locus has been mapped on chromosome 1p13.3 and 22q11.2, respectively, resulting in the reduced activity of a functional gene product (Hayes and Strange, 2000; Gonul et al., 2012). GSTP1 gene polymorphism is a single nucleotide polymorphism, a codon 105 A-G mutation at exon 5 in GSTP1 polymorphism leads to change in isoleucine (IIe) to valine (Val), is linked to a change in decreased enzymatic activity (Zimniak et al., 1994; Mastana et al., 2013). Thus, these three gene mutations may increase the risk of T2DM based on biological effects.

Independent and combined effects of GSTM1, GSTT1 and GSTP1 IIe105Val polymorphisms with the risk of T2DM have been described in a number of original investigations, and the association of GSTs gene polymorphisms with T2DM complications has also been reported. Of note, these results were contradictory. Five published meta-analyses (Tang et al., 2013; Yi et al., 2013; Zhang et al., 2013; Saadat, 2017; Nath et al., 2019) reported the GSTM1, GSTT1 and GSTP1 IIe105Val polymorphisms with T2DM risk. Three published meta-analyses have investigated the association between GSTs gene polymorphisms and T2DM complications (Orlewski and Orlewska, 2015; Sun et al., 2015; Nath et al., 2019). Yet, their studies had certain limitations. This could be due to the limited sample size or the lack of quality assessments of the literature. More importantly, past meta-analyses did not assess the reliability of positive results. As a result, we conducted a systematic review of the literature and a latest meta-analysis to assess whether GSTM1, GSTT1, and GSTP1 gene contribute to the risk of T2DM and related complications. We are also investigated the effect of gene-gene interactions on T2DM risk.

Materials and methods

Search strategy

The present study followed the Preferred Reporting Items for Systematic Review and Meta-Analysis (PRISMA) group’s guidelines. The PubMed, Embase, Wan-fang, and China National Knowledge Infrastructure (CNKI) databases were used to conduct the literature search (deadline to 16 May 2022). The following was utilized as a search strategy (polymorphism OR variant OR variation OR mutation OR SNP OR genome-wide association study OR genetic association study OR genotype OR allele) and (diabetes OR mellitus OR diabetes mellitus) AND (glutathione S-transferase M1 OR GSTM1 OR glutathione S-transferase T1 OR GSTT1 OR glutathione S-transferase P1 OR GSTP1). References of the retrieved articles and review articles on this topic were also carefully examined in additional suitable studies. If necessary, the respective authors were contacted by e-mail. Moreover, there were no restrictions or limitations on language.

Inclusion and exclusion criteria

The inclusion criteria were: 1) case–control or cohort design; 2) studies must examine the correlation between the individual and combined effects of GST polymorphisms with T2DM or T2DM complication risk; 3) studies with odds ratio (OR) and their 95% confidence interval (CI), or enough data to calculate these numbers. Exclusion criteria were: 1) case reports, reviews, commentaries and meta-analyses; 2) studies that did not give genotype frequencies or could not determine the number of genotypes and alleles; 3) overlapping data or incomplete data.

Data extraction

Two investigators extracted and double-checked data. When their conclusions differed, through discussion were reached a consensus. The first author’s name, year, country, source of control, match, race and other information were recorded, among which race was divided into “Asian,” “African,” “European,” “Caucasian,” “Indian” and “Mixed.”

Quality score assessment

The quality assessment scale was independently assessed by two investigators (Supplementary Table S2). The total score was 18 points, studies scoring >12 were high, while scores of ≤12 were regarded as low quality in this meta-analysis. The scale was influenced by selection (two points), control and case source (four points), diagnostic criteria for T2DM (two points), Ascertainment of control (two points), matching (two points), Genotyping examination (two points), Association assessment (two points), and Hardy-Weinberg balance (HWE).

Statistical analysis

Stata 12.0 software (Stata Corporation, College Station, TX) was used to calculate all statistical analyses. The strength of the association of GSTM1, GSTT1, and GSTP1 polymorphisms with T2DM and T2DM complications risk were estimated by calculating the crude OR and their 95% CI. We used the Chi-square-based Q-test and I2 value to evaluate heterogeneity (Higgins et al., 2003). The random-effects model was used when p was less than 0.10 and/or I2 was larger than 50%. if not, the fixed-effects model was applied (DerSimonian and Laird, 2015). Furthermore, a meta-regression analysis was used to investigate sources of heterogeneity (Baker et al., 2009). Subgroups were carried out based on ethnicity, control source, control type, and matching. Sensitivity analysis was performed using the following methods: 1) excluded one study at a time; 2) removed low-quality or Hardy–Weinberg Disequilibrium (HWD) studies; 3) retained only high quality and HWE research. HWE was evaluated by Chi-square goodness-of-fit test, p < 0.05 was defined as HWE, otherwise as HWD. To determine publication bias, Begg’s funnel plot (Begg and Mazumdar, 1994) and Egger’s test (Egger et al., 1997) were conducted. If publication bias existed, the missing studies would be filled employing a nonparametric “trim and fill” method (Dual and Tweedie, 2000).

The credibility of genetic association

To assess the reliability of statistically significant correlations, the false-positive report probability (FPRP) (Wacholder et al., 2004), Bayesian False Discovery Probability (BFDP) (Wakefield, 2008), and the Venice standard (Ioannidis et al., 2008) were used. The FPRP and BFDP values were set to 0.2 and 0.8, respectively, and statistically significant correlations satisfying the following criteria were defined as “positive results” (Theodoratou et al., 2012): 1) at least two genetic models were statistically important (GSTM1 and GSTT1 polymorphisms are exempt from this requirement); 2) FPRP <0.2 and/or BFDP <0.8; 3) statistical power >80%; and 4) I2 < 50%. Otherwise, a positive result is considered less reliable.

Results

Search results and study characteristics

207 records were returned by PubMed, EMBASE, Wan-fang and CNKI databases. By carefully reading titles and abstracts, 147 publications were eliminated. Two authors independently read the remaining 60 publications in their entirety. 12 is excluded from the full text. As a result, 48 studies (Fujita et al., 2000; Yang et al., 2004; Hayek et al., 2006; Wang et al., 2006; Hori et al., 2007; Yalin et al., 2007; Oniki et al., 2008; Nowier et al., 2009; Tiwari et al., 2009; Bid et al., 2010; Datta et al., 2010; Amer et al., 2011; Ramprasath et al., 2011; Tsai et al., 2011; Amer et al., 2012; Cilenšek et al., 2012; Gonul et al., 2012; Jana and Petrovic., 2012; Moasser et al., 2012; Dadbinpour et al., 2013; Grubisa et al., 2013; Mastana et al., 2013; Pinheiro et al., 2013; Vats et al., 2013; Abbasi et al., 2014; Al-Badran and Al-Mayah, 2014; Moasser et al., 2014; Purkait et al., 2014; Rao et al., 2014; Raza et al., 2014; Afrand et al., 2015; Rasheed et al., 2015; Stoian et al., 2015; Zaki et al., 2015; Etemad et al., 2016; Mergani et al., 2016; Mir et al., 2016; Rasheed et al., 2016; Ahmed and Al-Bachary, 2017; Azarova et al., 2018; de Lima et al., 2018; Abbas et al., 2019; Osman et al., 2019; Klusek et al., 2020; Pourkeramati et al., 2020; Gusti et al., 2021; Albeladi et al., 2022; Jamil et al., 2022) were included (Figure 1). There were 28 articles (involving 4,878 cases and 4,621 controls, Table 1) on the GSTM1 present/null polymorphism, 28 articles (involving 4,710 cases and 4,471 controls, Table 2) on the GSTT1 present/null polymorphism, 24 studies on the GSTP1 IIe105Val polymorphism (including 4,297 cases and 4,244 controls, Table 3),17 studies including 3,035 cases and 3,241 controls concerned the combined GSTM1 and GSTT1 gene (Table 4), nine studies with 1,953 cases and 1,905 controls concerned the combined GSTM1 and GSTP1 gene (Table 5), nine papers for the combined of GSTT1 and GSTP1 gene studies, 1,953 cases and 1,905 controls were involved (Table 6), there were seven studies for the combined of GSTM1, GSTT1, and GSTP1 gene, included 1,299 cases and 1,334 controls (Table 7). Finally, there were five high-quality investigations on GSTM1, GSTT1 and GSTP1, the combined of GSTM1 and GSTT1 had four high-quality investigations, three high-quality researches combined of GSTM1 and GSTP1 were found, the combined of GSTT1 and GSTP1 had two high-quality studies, as well as two high-quality research on the combined effects of GSTM1, GSTT1, and GSTP1 (Supplementary Table S3). 17 studies on the correlation between GSTs gene and complications of T2DM (Supplementary Table S8).

FIGURE 1

FIGURE 1

Flow diagram for searching studies in the current meta-analysis.

TABLE 1

VariablenCases/ControlsTest of associationTest of heterogeneityModel
Or (95%CI)PhI2 (%)
Overall284,878/4,6211.36 (1.08–1.72)<0.00183.8Random-effect
Ethnicity
 Indian81,689/1,4071.37 (0.85–2.20)<0.00188.7Random-effect
 Asian5647/1,0401.46 (0.81–2.64)<0.00186.3Random-effect
 Caucasian132,194/19401.44 (1.04–2.01)<0.00179.3Random-effect
Source of control
 HB254,330/4,1621.37 (1.08–1.74)<0.00183.1Random-effect
 PB3521/4591.29 (0.45–3.73)<0.00191.1Random-effect
Type of control
 HC254,683/3,8901.24 (0.99–1.56)<0.00180.9Random-effect
 NDC3275/7311.97 (1.01–3.86)0.00978.9Random-effect
Matching
 Yes183,769/3,2221.36 (1.02–1.82)<0.00187.3Random-effect
 No101,109/1,3991.36 (0.93–1.99)<0.00171.5Random-effect
geographic region
 Asia193,584/3,2541.48 (1.13–1.99)<0.00174.0Random-effect
 Africa3183/1671.10 (0.48–2.51)0.02268.8Random-effect
 Europe41,011/9661.29 (0.93–1.78)<0.00183.4Random-effect
Sensitivity analysis
 Quality score>12
  Overall51,633/1,3991.23 (0.65–2.32)<0.00193.5Random-effect
 Ethnicity
  Indian3933/7791.58 (0.60–4.19)<0.00194.9Random-effect
  Caucasian2700/6200.84 (0.39–1.81)0.01284.3Random-effect
 geographic region
  Asia3933/7791.58 (0.60–4.19)<0.00194.9Random-effect
 Source of control
  HB31,212/9901.21 (0.52–2.85)<0.00194.6Random-effect
  PB2421/4091.22 (0.26–5.70)<0.00195.5Random-effect
Matching and Quality score>12
  Overall51,633/1,3991.23 (0.65–2.32)<0.00193.5Random-effect

Meta-analysis of the association of GSTM1 polymorphism with risk of T2DM.

TABLE 2

VariablenCases/ControlsTest of associationTest of heterogeneityModel
OR (95% CI)PhI2 (%)
Overall284710/44711.45 (1.18–1.79)<0.00171.5Random-effect
Ethnicity
Indian71582/13551.71 (1.20–2.43)<0.00166.9Random-effect
Asian4562/8421.46 (1.23–1.88)0.9490.0Random-effect
Caucasian152345/13551.21 (0.86–1.70)<0.00175.8Random-effect
Source of control
HB244138/39621.60 (1.33–1.92)<0.00157.4Random-effect
PB4572/5090.69 (0.23–2.06)<0.00191.1Random-effect
Type of control
HC254420/38881.43 (1.13–1.80)<0.00174.3Random-effect
NDC3290/5831.63 (1.15–2.30)0.6200.0Random-effect
Matching
Yes193793/32181.28 (0.98–1.67)<0.00177.5Random-effect
No9917/11991.98 (1.36–2.89)0.0360.0Random-effect
geographic region
Asia183363/32041.55 (1.16-2.08)<0.00184.8Random-effect
Africa4283/2170.63 (0.43-1.00)0.6220.0Random-effect
Europe41011/9661.50 (0.91-2.46)<0.00183.4Random-effect
Sensitivity analysis
Quality score >12
Overall51633/13991.34 (0.79–2.28)<0.00187.0Random-effect
Ethnicity
Indian3933/7791.69 (0.86–3.31)0.00186.1Random-effect
Caucasian2700/6200.92 (0.29–2.92)0.00190.6Random-effect
geographic region
Asia3933/7791.69 (0.86-3.31)0.00186.1Random-effect
Source of control
HB31212/9901.66 (0.86–3.23)0.00186.8Random-effect
PB2421/4090.94 (0.28–3.09)0.00190.7Random-effect
Matching and quality score >12
Overall51633/13991.34 (0.79–2.28)<0.00187.0Random-effect

Meta-analysis of the association of GSTT1 polymorphism with risk of T2DM.

HB = hospital-based studies, PB = population-based studies, HC = Healthy control, NDC = Non-diabetic controls

TABLE 3

Variablen (Cases/Controls)Val/Val vs. IIe/IIeIIe/Val vs. IIe/IIeVal/Val vs. IIe/IIe + IIe/ValVal/Val + IIe/Val vs. IIe/IIeVal vs. IIe
Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)
Overall24 (4,297/4,244)1.37 (0.99–1.88)0.002/53.21.24 (1.02–1.50)0.000/70.01.27 (0.98–1.61)0.104/29.21.27 (1.05–1.53)0.000/72.11.22 (1.05–1.42)0.000/71.1
Ethnicity
 Caucasian14 (2049/2026)1.11 (0.76–1.62)0.057/42.71.13 (0.88–1.46)0.000/67.61.00 (0.79–1.26)0.534/0.01.16 (0.90–1.48)0.000/70.11.11 (0.92–1.33)0.000/64.1
 Indian6 (1751/1,328)2.17 (1.59–2.96)0.441/0.01.54 (1.22–1.94)0.081/48.91.84 (1.37–2.47)0.431/0.01.61 (1.32–1.97)0.167/36.01.47 (1.27–1.70)0.238/26.3
Source of control
 HB18 (3,605/3,182)1.24 (0.89–1.73)0.007/52.91.16 (0.91–1.47)0.000/75.31.14 (0.90–1.45)0.224/19.71.18 (0.93–1.49)0.000/77.61.33 (0.96–1.34)0.000/74.0
 PB3 (496/902)2.56 (1.47–4.44)1.42 (0.98–2.07)0.220/34.12.29 (1.36–3.85)1.45 (1.13–1.85)0.540/0.01.48 (1.20–1.83)0.566/0.0
Type of control
 HC20 (3,652/3,394)1.38 (0.96–1.98)0.001/59.31.34 (1.09–1.64)0.000/69.51.23 (0.94–1.63)0.046/39.81.37 (1.11–1.69)0.000/73.51.28 (1.09–1.51)0.000/72.5
 NDC4 (645/850)1.28 (0.65–2.51)0.635/0.00.77 (0.52–1.15)0.168/44.01.41 (0.73–2.72)0.790/0.00.89 (0.64–1.24)0.171/40.10.91 (0.66–1.25)0.174/42.9
Matching
 Yes14 (3,080/2,847)1.49 (1.01–2.20)0.001/63.21.35 (1.00–1.72)0.000/74.71.32 (0.99–1.77)0.051/42.71.40 (1.09–1.79)0.000/77.81.31 (1.09–1.57)0.000/76.3
 No10 (1,217/1,397)1.09 (0.61–1.96)0.228/25.21.04 (0.76–1.43)0.021/55.51.12 (0.72–1.76)0.420/1.31.07 (0.81–1.42)0.020/54.41.07 (0.83–1.36)0.025/54.3
Geographic region
 Asia15 (3,967/3,513)1.76 (1.24–2.50)0.112/32.91.38 (1.10–1.73)0.000/64.91.52 (1.45–2.00)0.299/14.11.42 (1.34–1.78)0.000/66.51.32 (1.12–1.56)0.001/60.4
 Europe5 (1,151/1,066)0.99 (0.62–1.54)0.124/44.60.88 (0.68–1.14)0.140/42.21.00 (0.71–1.43)0.269/22.80.92 (0.68–1.24)0.036/61.10.98 (0.77–1.26)0.016/67.1
HWE
 Yes4 (609/556)1.34 (0.51–3.51)0.038/64.51.52 (0.81–2.84)0.003/78.41.27 (0.62–2.61)0.162/41.61.50 (0.82–2.75)0.003/79.01.38 (0.88–2.17)0.005/77.0
 No19 (3,625/3,193)1.35 (0.96–1.90)0.007/51.71.19 (0.97–1.43)0.000/66.91.24 (0.95–1.60)0.121/30.31.23 (1.00–1.51)0.000/71.31.19 (1.01–1.40)0.000/69.6
Sensitivity analysis
 HWD
  Overall19 (3,625/3,193)1.35 (0.96–1.90)0.007/51.71.19 (0.97–1.43)0.000/66.91.24 (0.95–1.60)0.121/30.31.23 (1.00–1.51)0.000/71.31.19 (1.01–1.40)0.000/69.6
 Ethnicity
  Indian4 (1,303/899)2.17 (1.15–3.12)0.639/0.01.41 (1.06–1.86)0.111/50.11.88 (1.33–2.65)0.560/0.01.52 (1.16–1.99)0.116/49.31.44 (1.17–1.78)0.117/49.0
  Caucasian13 (1989/1986)1.06 (0.72–1.55)0.065/42.61.14 (0.88–1.47)0.000/65.90.97 (0.76–1.23)0.456/0.01.15 (0.89–1.49)0.000/68.91.10 (0.92–1.32)0.002/62.2
 geographic region
  Asia11 (2060/1778)1.63 (1.07–2.50)0.091/40.01.22 (0.94–1.59)0.001/66.31.44 (1.03–2.01)0.253/20.71.27 (0.97–1.66)0.000/69.31.21 (0.99–1.46)0.003/62.9
  Europe5 (1,151/1,066)0.99 (0.62–1.54)0.124/44.60.88 (0.68–1.14)0.140/42.21.00 (0.71–1.43)0.269/22.80.92 (0.68–1.24)0.036/61.10.98 (0.77–1.26)0.016/67.1
 Source of control
 HB13 (2,971/2,468)1.20 (0.84–1.71)0.021/49.71.17 (0.91–1.51)0.000/72.31.10 (0.86–1.40)0.311/13.31.20 (0.93–1.54)0.000/75.91.14 (0.95–1.37)0.000/72.4
 PB2 (433/497)2.56 (1.47–4.44)1.42 (0.98–2.07)0.218/34.12.29 (1.36–3.85)1.51 (1.15–1.99)0.412/0.01.48 (1.20–1.83)0.566/0.0
 Type of control
 HC17 (3,963/3,999)1.43 (0.90–2.20)0.001/62.61.32 (1.05–1.67)0.000/68.61.21 (0.89–1.65)0.053/42.41.37 (1.07–1.75)0.000/73.81.28 (1.06–1.55)0.000/73.4
 NDC3 (582/445)1.28 (0.65–2.51)0.635/0.00.77 (0.52–1.15)0.168/44.01.41 (0.73–2.72)0.790/0.00.81 (0.54–1.22)0.142/48.90.91 (0.66–1.25)0.174/42.9
 Matching
 Yes12 (2,632/2,418)1.40 (0.92–2.14)0.002/65.51.24 (0.94–1.64)0.000/75.61.26 (0.93–1.72)0.067/43.71.32 (1.01–1.74)0.000/78.31.26 (1.03–1.55)0.000/77.2
 No7 (1,165/1826)1.14 (0.61–2.12)0.321/14.51.03 (0.77–1.38)0.115/41.51.15 (0.65–2.06)0.350/10.31.05 (0.79–1.39)0.112/41.81.04 (0.83–1.31)0.154/36.0
 Quality score>12
 Overall5 (1,645/1,487)1.64 (0.76–3.57)0.001/81.91.39 (0.99–1.96)0.002/76.81.46 (0.77–2.76)0.006/75.81.47 (1.00–2.15)0.000/82.51.39 (1.01–1.93)0.000/85.0
 Ethnicity
 Caucasian2 (712/708)0.72 (0.48–1.10)1.21 (0.57–2.55)0.012/84.20.77 (0.52–1.16)1.19 (0.55–2.59)0.008/85.71.15 (0.59–2.25)0.011/84.5
 Indian3 (933/779)2.21 (1.51–3.24)0.444/0.01.53 (1.09–2.15)0.103/56.11.88 (1.30–2.73)0.357/2.801.65 (1.21–2.26)0.123/52.31.54 (1.18–2.00)0.087/59.0
 Geographic region
 Asia3 (933/779)2.21 (1.51–3.24)0.444/0.01.53 (1.09–2.15)0.103/56.11.88 (1.30–2.73)0.357/2.801.65 (1.21–2.26)0.123/52.31.54 (1.18–2.00)0.087/59.0
 Source of control
 HB3 (1,212/990)1.35 (0.54–3.39)0.010/78.41.38 (0.80–2.39)0.001/86.21.16 (0.59–2.27))0.060/64.41.43 (0.78–2.62)0.000/89.51.33 (0.82–2.18)0.000/90.1
 PB2 (443/497)2.56 (1.47–4.44)1.42 (0.98–2.07)0.218/34.12.29 (1.36–3.85)1.51 (1.15–1.99)0.412/0.01.48 (1.20–1.83)0.566/0.0

Meta-analysis of the association of GSTP1 polymorphism with risk of T2DM.

HB, hospital-based studies; PB, population-based studies; HC, healthy control; NDC, Non-diabetic controls.

TABLE 4

VariableN (case/Control)Model 1Model 2Model 3Model 4Model 5Model 6
Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)
Overall17 (3,035/3,241)1.25 (0.82–1.92)0.016/57.62.28 (1.72–3.04)0.004/54.72.19 (1.64–2.93)0.005/54.01.41 (1.15–1.72)0.004/56.21.89 (1.06–1.33)0.536/0.01.75 (1.41–2.18)0.124/30.6
Ethnicity
 Caucasian8 (1,638/1,364)0.88 (0.50–1.55)0.184/38.01.55 (0.99–2.44)0.034/65.42.12 (1.31–3.43)0.005/65.71.47 (1.02–2.13)0.037/57.81.25 (1.02–1.53)0.399/2.701.39 (0.98–1.99)0.203/30.9
 Indian5 (969/988)1.30 (0.60–2.83)0.158/49.81.06 (0.44–2.54)0.081/67.12.39 (1.40–4.09)0.034/61.61.37 (0.93–2.01)0.007/71.71.14 (0.97–1.35)0.377/5.202.08 (1.40–3.08)0.175/37.0
Source of control
 HB14 (2,614/2,782)1.18 (0.75–1.86)0.014/60.31.28 (0.91–1.79)0.007/64.01.90 (1.41–2.56)0.033/46.61.27 (1.03–1.58)0.220/50.61.13 (1.00–1.29)0.570/0.01.59 (1.29–1.96)0.339/10.8
 PB3 (521/459)3.92 (2.56–6.00)0.592/0.02.04 (1.54–2.70)0.652/0.01.40 (1.11–1.77)0.612/0.02.68 (1.68–4.29)0.296/17.9
Type of control
 HC15 (2,938/2,681)1.16 (0.74–1.82)0.019/58.41.25 (0.90–1.74)0.010/62.22.15 (1.58–2.92)0.004/56.41.36 (1.10–1.69)0.008/54.11.19 (1.05–1.33)0.460/0.01.69 (1.37–2.08)0.204/23.1
Matching
 Yes14 (2,847/2,615)1.22 (0.74–2.01)0.014/62.21.22 (0.83–1.79)0.007/66.42.34 (1.73–3.17)0.028/47.81.40 (1.11–1.76)0.002/63.21.18 (1.03–1.33)0.413/3.31.90 (1.53–2.36)0.337/11.0
 No3 (288/626)1.37 (0.45–4.17)0.087/65.91.65 (1.07–2.53)0.452/0.01.60 (0.63–4.05)0.018/75.11.45 (0.98–2.14)0.363/1.21.25 (0.89–1.77)0.495/0.01.28 (0.70–2.35)0.068/62.8
Geographic region
 Asia11 (2,761/2,943)1.18 (0.87–1.60)0.368/7.92.25 (1.70–3.0)0.136/34.01.42 (0.96–2.10)0.006/66.91.42 (1.10–1.81)0.007/60.41.18 (1.04–1.34)0.457/0.01.89 (1.53–2.27)0.552/0.0
 Europe3 (811/745)2.52 (0.85–7.47)0.013/77.0
Sensitivity analysis
 Quality score
 >124 (1,121/1,129)3.40 (2.42–4.79)0.669/0.02.00 (1.54–2.59)0.616/0.01.39 (1.12–1.74)0.61/0.02.21 (1.40–3.50)0.299/17.1
 Sample size
 >20011 (2,451/2,644)1.40 (0.87–2.25)0.030/59.51.17 (0.86–1.59)0.096/46.52.30 (1.78–2.97)0.162/29.91.34 (1.05–1.70)0.016/57.61.16 (1.02–1.32)0.518/0.01.80 (1.47–2.19)0.768/0.0
 HWE
 Yes15 (2,589/2,812)1.31 (0.80–2.14)0.014/60.32.34 (1.68–3.26)0.006/55.61.42 (1.04–1.94)0.040/52.31.56 (1.25–1.94)0.040/46.11.27 (1.11-1-46)0.667/0.01.69 (1.31–2.19)0.118/34.0

Meta-analysis of the combined effects of GSTM1 present/null and GSTT1 present/null on T2DM risk.

Model 1 = M1 present/T1 null vs. M1 present/T1 present, Model 2 = M1 null/T1 present vs. M1 present/T1 present, Model 3 = M1 null/T1 null vs. M1 present/T1 present, Model 4 = All one risk genotypes vs. M1 present/T1 present, Model 5 = All risk genotypes vs. M1 present/T1 present, Model 6 = M1 null/T1 null vs. M1 present/T1 present + M1 present/T1 null + M1 null/T1 present, HB, hospital-based studies; PB, population-based studies.

The bold values indicated that these results are statistically significant.

TABLE 5

VariableSample sizeModel 1Model 2Model 3Model 4Model 5Model 6
Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)
Overall9 (1953/1905)1.40 (0.94–2.08)0.002/69.80.98 (0.58–1.65)0.001/73.31.15 (0.82–1.50)0.114/41.51.84 (1.13–3.00)0.001/70.11.26 (0.89–1.76)0.007/66.21.37 (1.03–1.82)0.290/18.3
Ethnicity
 Caucasian4 (836/767)1.11 (0.60–2.07)0.192/39.51.05 (0.49–2.22)0.139/49.31.02 (0.67–1.56)0.992/0.02.72 (1.78–4.16)0.970/0.01.31 (0.87–1.96)0.944/0.02.11 (1.29–3.44)0.944/0.0
 Indian4 (869/938)1.45 (0.74–2.87)0.001/82.30.93 (0.44–1.96)0.000/83.71.18 (0.78–1.77)0.023/68.41.28 (0.53–3.06)0.00/85.21.21 (0.71–2.04)0.001/82.91.17 (0.86–1.59)0.348/9.1
Source of control
 HB8 (1,632/1,596)1.23 (0.85–1.78)0.025/58.50.91 (0.47–1.76)0.001/75.7.31.05 (0.82–1.33)0.378/6.11.67 (0.96–2.91)0.005/67.51.11 (0.82–1.50)0.133/40.81.36 (0.91–2.03)0.198/31.7
Matching
 Yes8 (1899/1854)1.48 (0.97–2.25)0.002/71.90.91 (0.52–1.61)0.001/77.01.14 (0.84–1.55)0.070/50.91.75 (1.04–2.97)0.001/73.91.25 (0.85–1.82)0.003/71.81.33 (0.98–1.79)0.260/23.2
Sensitivity analysis
Quality score
 >123 (1,021/1,029)2.57 (1.53–4.31)0.204/37.90.70 (0.16–2.98)0.002/89.81.41 (0.88–2.56)0.139/54.31.88 (0.84–4.18)0.005/80.91.48 (0.69–3.19)0.011/84.61.10 (0.58–2.06)0.115/59.6
HWE
  Yes7 (1,505/1,476)1.72 (1.15–2.58)0.037/57.80.90 (0.48–1.70)0.008/70.81.30 (1.01–1.67)0.400/1.202.24 (1.41–3.54)0.050/54.81.42 (1.01–2.01)0.112/46.61.48 (1.03–2.15)0.217/30.7

Meta-analysis of the combined effects of GSTM1 present/null and GSTP1 IIe105Val on T2DM risk.

Model 1 = M1 null/P1 IIe/IIe, vs. M1 present/P1 IIe/IIe, Model 2 = M1 present/P1 Val* vs. M1 present/P1 IIe/IIe, Model 3 = (M1 null/P1 IIe/IIe + M1 present/P1 Val*) vs. M1 present/P1 IIe/IIe; Model 4 = M1 null/P1 Val* vs. M1 present/P1 IIe/IIe, Model 5 = All risk genotypes vs. M1 present/P1 IIe/IIe, Model 6 = M1 null/P1 Val* vs. (M1 present/P1 IIe/IIe + M1 null/P1 IIe/IIe + M1 Present/P1 Val*), HB, hospital-based studies; PB, population-based studies.

TABLE 6

VariableSample sizeModel 1Model 2Model 3Model 4Model 5Model 6
Or (95%CI)Ph/I2Or (95%CI)Ph/I2Or (95%CI)Ph/I2Or (95%CI)Ph/I2Or (95%CI)Ph/I2Or (95%CI)Ph/I2
Overall9 (1953/1905)1.20 (0.72–1.20)0.002/69.21.02 (0.68–1.53)0.009/65.00.98 (0.73–1.31)0.054/51.61.68 (0.90–3.15)0.002/71.21.12 (0.79–1.59)0.003/70.11.47 (0.94–2.28)0.098/43.9
Ethnicity
 Caucasian3 (236/247)0.79 (0.34–1.86)0.118/53.21.41 (0.92–2.17)0.534/0.001.14 (0.78–1.66)0.832/0.001.22 (0.48–3.10)0.129/51.21.15 (0.80–1.66)0.911/0.01.04 (0.45–2.42)0.172/43.1
 Indian4 (869/938)1.45 (0.59–3.55)0.002/76.90.82 (0.46–1.48)0.005/76.90.89 (0.56–1.42)0.009/74.02.04 (0.92–4.52)0.006/76.21.09 (0.63–1.90)0.000/84.91.75 (1.10–2.78)0.203/34.9
Source of control
 HB7 (1,032/1,076)0.85 (0.61–1.18)0.434/0.001.01 (0.60–1.72)0.005/70.50.90 (0.65–1.25)0.099/46.01.46 (0.68–3.12)0.004/70.61.03 (0.70–1.52)0.010/67.01.36 (0.77–2.41)0.083/48.7
Matching
 Yes7 (1,299/1,334)1.29 (0.67–2.49)0.002/72.90.97 (0.63–1.50)0.006/69.20.98 (0.71–1.36)0.031/59.41.58 (0.78–3.20)0.001/75.91.11 (0.76–1.64)0.001/75.11.40 (0.85–2.30)0.065/52.0
Sensitivity analysis
Quality score
 ≤125 (680/676)0.85 (0.56–1.28)0.303/17.51.30 (0.95–1.78)0.598/0.001.07 (0.82–1.38)0.1962/0.01.87 (0.92–3.81)0.029/62.91.27 (1.00–1.61)0.947/0.001.47 (0.77–2.81)0.060/55.8
HWE
 Yes6 (905/956)1.30 (0.65–2.60)0.001/74.61.00 (0.58–1.74)0.004/74.30.96 (0.62–1.48)0.015/67.61.24 (0.50–3.03)0.002/76.21.02 (0.61–1.71)0.001/78.91.23 (0.71–2.13)0.147/41.2

Meta-analysis of the combined effects of GSTT1 present/null and GSTP1 IIe105Val on T2DM risk.

Model 1 = T1 null/P1 IIe/IIe, vs. T1 present/P1 IIe/IIe, Model 2 = T1 present/P1 Val* vs. T1 present/P1 IIe/IIe, Model 3 = (T1 null/P1 IIe/IIe + T1 present/P1 Val*) vs. T1 present/P1 IIe/IIe, Model 4 = T1 null/P1 Val* vs. T1 present/P1 IIe/IIe, Model 5 = All risk genotypes vs. T1 present/P1 IIe/IIe, Model 6 = T1 null/P1 Val* vs. (T1 present/P1 IIe/IIe + T1 null/P1 IIe/IIe + T1 Present/P1 Val*), HB, hospital-based studies; PB, population-based studies.

The bold values indicated that these results are statistically significant.

TABLE 7

VariableSample sizeModel 1Model 2Model 3Model 4Model 5Model 6Model 7Model 8Model 9Model 10
Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)Or (95%CI)Ph/I2 (%)
Overall7 (1,299/1,334)1.47 (0.10–1.96)0.218/27.60.80 (0.55–1.14)0.242/24.50.73 (0.44–1.23)0.009/67.20.91 (0.70–1.19)0.143/39.41.61 (0.82–3.17)0.010/64.51.17 (0.75–1.82)0.098/46.11.06 (0.42–2.66)0.002/74.01.18 (0.83–1.08)0.065/51.92.71 (1.56–4.72)0.276/20.92.69 (1.74–4.17)0.382/5.4
Sensitivity analysis
 Quality score
  >122 (421/509)1.81 (1.02–3.22)0.185/43.20.64 (0.17–2.40)0.063/71.10.44 (0.30–0.64)0.474/0.00.84 (0.38–1.86)0.015/83.31.08 (0.45–2.60)0.555/0.00.81 (0.54–1.22)0.993/0.00.66 (0.37–1.18)0.682/0.00.80 (0.56–1.1200.870/0.01.75 (1.00–3.07)0.697/0.02.24 (1.32–3.80)0.502/0.0
 HWE
  Yes5 (851/905)1.61 (1.15–2.25)0.256/24.80.76 (0.42–1.39)0.140/42.20.59 (0.33–1.04)0.096/52.70.89 (0.57–1.38)0.071/57.32.02 (0.65–6.31)0.004/74.21.32 (0.66–2.64)0.028/67.20.59 (0.35–1.00)0.796/0.01.01 (0.64–1.60)0.128/47.21.87 (1.12–3.11)0.757/0.02.19 (1.35–3.54)0.577/0.0

Meta-analysis of the combined effects of GSTM1 present/null, GSTT1 present/null and GSTP1 present/null on T2DM risk.

Model 1 = M1 null/T1 present/P1 IIe/IIe, vs. M1 present/T1 present/P1 IIe/IIe, Model 2 = M1 present/T1 null/P1 IIe/IIe, vs. M1 present/T1 present/P1 IIe/IIe, Model 3 = M1 present/T1 present/P1 Val 1 vs. M1 present/T1 present/P1 IIe/IIe, Model 4 = all one high-risk genotype vs. vs. M1 present/T1 present/P1 IIe/IIe, Model 5 = M1 null/T1 null/P1 IIe/IIe, vs. M1 present/T1 present/P1 IIe/IIe, Model 6 = M1 null/T1 present/P1 Val 1 vs. M1 present/T1 present/P1 IIe/IIe, Model 7 = M1 present/T1 null/P1 Val1 vs. M1 present/T1 present/P1 IIe/IIe, Model 8 = all two high-risk genotype vs. M1 present/T1 present/P1 IIe/IIe, Model 9 = M1 null/T1 null/P1 Val 1 vs. M1 present/T1 present/P1 IIe/IIe, Model 10 = M1 null/T1 null/P1 Val 1 vs. M1 present/T1 present/P1 IIe/IIe + all one high-risk genotype + all two high-risk genotypes.

The bold values indicated that these results are statistically significant.

Quantitative synthesis

Overall, the individuals who carried the GSTM1 null genotype had a significantly increased T2DM risk (OR = 1.36, 95% CI: 1.04–1.72, Table 1). When subgroup analyses were performed, the same substantial association were found in Caucasians (OR = 1.44, 95% CI: 1.04–2.01), hospital-based studies (OR = 1.37, 95% CI: 1.08–1.74), matching (OR = 1.36, 95% CI: 1.02–1.82), non-diabetic controls (OR = 1.97, 95% CI: 1.01–3.86), and Asia (OR = 1.48, 95% CI: 1.13–1.99).

Overall, a notable connection was discovered between the GSTT1 polymorphism and an increased risk of T2DM (OR = 1.45, 95% CI: 1.18–1.79, Table 2). Similarly, high correlation was found in subgroup analyses in Indians (OR = 1.71, 95% CI: 1.20–2.43), Asians (OR = 1.46, 95% CI: 1.23–1.88), hospital-based studies (OR = 1.60, 95% CI: 1.33–1.92), healthy controls (OR = 1.43, 95% CI: 1.13–1.80), non-diabetic controls (OR = 1.63, 95% CI: 1.15–2.30), non-matching (OR = 1.98, 95% CI: 1.36–2.89), and Asia (OR = 1.55, 95% CI: 1.16–2.08).

The GSTP1 polymorphism were also found to be associated with increased risk of T2DM (IIe/Val vs. IIe/IIe: OR = 1.24, 95% CI = 1.02–1.50, Val/Val + IIe/Val vs. IIe/IIe: OR = 1.27, 95% CI = 1.05–1.53 and Val vs. IIe: OR = 1.22, 95% CI = 1.05–1.42, Table 3). As for the subgroup analysis of populations in different Ethnicity, high associations were observed in Caucasians (Val/Val vs. IIe/IIe: OR = 2.17, 95% CI = 1.59–2.96, IIe/Val vs. IIe/IIe: OR = 1.54, 95% CI = 1.22–1.94, Val/Val vs. IIe/IIe + IIe/Val: OR = 1.84, 95% CI = 1.37–2.47, Val/Val + IIe/Val vs. IIe/IIe: OR = 1.61, 95% CI = 1.32–1.97, Val vs. IIe: OR = 1.47, 95% CI = 1.27–1.70). In addition, similar result was found in the population-based studies, healthy controls and Asia.

This study revealed a strong association combination of GSTM1 and GSTT1 polymorphisms and T2DM susceptibility in overall analysis (model 2: OR = 2.28, 95% CI = 1.72–3.04; model 3: OR = 2.19, 95% CI = 1.64–2.93; model 4: OR = 1.41, 95% CI = 1.15–1.72; model 5: OR = 1.89, 95% CI = 1.06–1.33; model 6: OR = 1.75, 95% CI = 1.41–2.18, Table 4). Furthermore, substantially association were showed in Caucasians (model 3: OR = 2.12, 95% CI = 1.31–3.43; model 4: OR = 1.47, 95% CI = 1.02–2.13; model 5: OR = 1.25, 95% CI = 1.02–1.53), Indians, hospital-based studies, population-based studies, healthy controls, non-matching, matching and Asia.

An increased risk of T2DM was yielded on the combined effects of both GSTM1 and GSTP1 (model 4: OR = 1.84, 95% CI = 1.13–3.00, model 6: OR = 1.37, 95% CI = 1.03–1.82, Table 5) in overall analysis. subgroups also showed an increased T2DM risk in Caucasians (model 4: OR = 2.72, 95% CI = 1.78–4.16, model 6: OR = 2.11, 95% CI = 1.29–3.44) and matching controls (model 4: OR = 1.75, 95% CI = 1.04–2.97).

Combinations of GSTT1 and GSTP1 polymorphism were not associated with T2DM risk in overall analysis. However, when subgroup analysis for ethnicity was stratified, important association were found in Indians (model 6: OR = 1.75, 95% CI = 1.10–2.78, Table 6).

In overall populations, combinations of GSTM1, GSTT1 and GSTP1 polymorphisms was considerably associated with T2DM risk (model 1: OR = 1.47, 95% CI = 1.10–1.96, model 9: OR = 2.71, 95% CI = 1.56–4.72, model 10: OR = 2.69, 95% CI = 1.74–4.17, Table 7).

Overall, a notable association was discovered between the GSTT1 null genotype and an increased risk of T2DM complications (OR = 1.52, 95% CI: 1.16–1.99, Figure 2). No association was found between GSTM1 and GSTP1 polymorphisms and the risk of T2DM complications (Figures 3, 4).

FIGURE 2

FIGURE 2

Forest plot of the association between GSTT1 polymorphism and T2DM complications risk.

FIGURE 3

FIGURE 3

Forest plot of the association between GSTM1 polymorphism and T2DM complications risk.

FIGURE 4

FIGURE 4

Forest plot of the association between GSTPI polymorphism and T2DM complications risk [(A): VV vs. II; (B): IV vs. II; (C): VV vs. IV + II; (D): VV + IV vs. II; (E): V vs. I].

Heterogeneity and sensitivity analyses

We assessed heterogeneity source by applying a meta-regression analysis. The source of heterogeneity for the GSTT1 null genotype and the combined effects of GSTT1 and GSTP1 was revealed to be the source of controls (p = 0.028 and model 1: p = 0.037). HWE (model 7: p = 0.045) was the source of heterogeneity between the combined effects of GSTM1, GSTT1, and GSTP1 gene and T2DM risk.

Sensitivity analysis was estimated by applying three methods. First, results did not change when removing a single study each time to appraise the robustness. Second, there was no significant connection individual GSTM1 and GSTT1 polymorphisms when studies of low quality and non-matching were eliminated (Tables 1, 2). Moreover, concerning the GSTP1 polymorphism, a noteworthy association with T2DM risk was found in the overall analysis (Val vs. IIe: OR = 1.39, 95% CI = 1.01–1.93, Table 3), Indians (Val/Val vs. IIe/IIe: OR = 2.21, 95% CI = 1.51–3.24, IIe/Val vs. IIe/IIe: OR = 1.53, 95% CI = 1.09–2.15, Val/Val vs. IIe/IIe + IIe/Val: OR = 1.88, 95% CI = 1.30–2.73, Val/Val + IIe/Val vs. IIe/IIe: OR = 1.65, 95% CI = 1.21–2.26, Val vs. IIe: OR = 1.54, 95% CI = 1.18–2.00), population-based studies (Val/Val vs. IIe/IIe: OR = 2.56, 95% CI = 1.47–4.44, Val/Val vs. IIe/IIe + IIe/Val: OR = 2.29, 95% CI = 1.36–3.85, Val/Val + IIe/Val vs. IIe/IIe: OR = 1.51, 95% CI = 1.15–1.99, Val vs. IIe: OR = 1.48, 95% CI = 1.20–1.83) and Asia when we retained research of low quality and HWD in control.

When we restricted only high-quality studies, notably increased T2DM risk was found the combinations of GSTM1 and GSTT1 polymorphisms (model 3: OR = 3.40, 95% CI = 2.42–4.79; model 4: OR = 2.00, 95% CI = 1.54–2.59; model 5: OR = 1.39, 95% CI = 1.12–1.74; model 6: OR = 2.21, 95% CI = 1.40–3.50, Table 4). Significantly increased T2DM risk was observed the GSTM1-GSTT1 polymorphisms (model 2: OR = 2.34, 95% CI = 1.68–3.26; model 3: OR = 1.42, 95% CI = 1.04–1.94; model 4: OR = 1.56, 95% CI = 1.25–1.94; model 5: OR = 1.27, 95% CI = 1.11–1.46; model 6: OR = 1.69, 95% CI = 1.31–2.19, Table 4) When we excluded HWD.

When we removed low-quality studies, significantly increased T2DM risk was found in overall analysis (model 1: OR = 2.57, 95% CI = 1.53–4.31, Table 5), and when we removed HWD, significantly association was revealed in the overall analysis (model 1: OR = 1.72, 95% CI = 1.15–2.58; model 3: OR = 1.30, 95% CI = 1.01–1.67; model 4: OR = 2.24, 95% CI = 1.41–3.54; model 5: OR = 1.42, 95% CI = 1.01–2.01; model 6: OR = 1.48, 95% CI = 1.03–2.15) combined effects on GSTM1 and GSTP1 polymorphisms.

Concerning the combination effects of GSTM1, GSTT1, and GSTP1 IIe105Val polymorphisms, considerably increased T2DM risk was observed (model 1: OR = 2.81, 95% CI = 1.02–3.22, model 9: OR = 1.75, 95% CI = 1.00–3.07, model 10: OR = 2.24, 95% CI = 1.32–3.80, Table 7) when we excluded low-quality studies. Significantly increased T2DM risk was revealed (model 1: OR = 1.61, 95% CI = 1.15–2.25, model 9: OR = 1.87, 95% CI = 1.12–3.11, model 10: OR = 2.19, 95% CI = 1.35–3.54) when HWD was removed in overall analysis. On the contrary, no notable was detected on the combined effects of GSTT1 and GSTP1 polymorphisms (Table 6).

Publication bias

No publication bias was observed between GSTs polymorphisms with T2DM risk through the Begg’s funnel plot and Egger’s test (Figures not shown).

Discussion

Oxidative stress is considered to play an important role in the pathogenesis of diabetes mellitus and related complications. GST is a member of phase II metabolic isoenzyme group, which has the effect of regulating various cytotoxicity, genotoxicity and antioxidant. The risk genotypes expressed by each gene (GSTM1 null, GSTT1 null, and GSTP1 Val/Val) may decrease the antioxidant activity of the enzyme, thus increasing the susceptibility of T2DM and other diseases. There was much important evidence showing that GSTM1, GSTT1, and GSTP1 gene polymorphisms were potential genetic factors for T2DM. Yet, previous research remains contradictory. In addition, seven previously published meta-analyses also had certain limitations, therefore the association of the GSTs polymorphisms with T2DM risk needs to be further evaluated. Furthermore, to our knowledge, this research is the first large-scale meta-analysis to investigate the relationship between individual and combined GSTM1, GSTT1, and GSTP1 genes polymorphism of T2DM susceptibility, and utilizing FPRP and BFDP to explore the validity of positive results.

Overall, on the individual and combined effects of GSTM1, GSTT1, and GSTP1 polymorphisms, a statistically significantly increased T2DM risk was found. We found that GSTT1 null genotype increased the risk of T2DM complications. However, GSTM1 and GSTP1 polymorphism were not associated with T2DM complications risk. The GSTM1 null genotype was correlated to a significantly increased risk of T2DM in Caucasians, Asia, hospital-based studies, and non-diabetic controls. The GSTT1 null genotype was connected with an increased risk of T2DM in Asians, Indians, Asia, and so on. GSTP1 IIe105Val polymorphism was correlated with an increased risk of T2DM in Indians, Asia, population-based studies, and healthy controls. The GSTM1-GSTT1 double null genotype significantly increased the risk of T2DM in both subgroup and sensitivity analysis. The effects of GSTM1-GSTP1 gene polymorphisms were related to the risk of T2DM in Caucasians and matching. However, no relationship was identified between GSTT1-GSTP1 polymorphisms and T2DM susceptibility, and, interestingly, Indians had a significantly increased risk. For the triple genotype combinations, we also discovered an association with T2DM. These findings showed that the same gene has variable susceptibility to T2DM in different race groups. It is worth considering that the obvious inconsistencies could be due to sample size, GSTs prevalence and environmental factors (lifestyle). Additionally, because T2DM is a chronic metabolic disease with multiple factors and genes, and different genetic backgrounds may lead to such differences. We found that GSTs gene polymorphism increased the risk of T2DM in Asia, but not in Europe. Most of the original research were conducted in Asia, which may lead to this difference. Therefore, a larger sample size is needed to study the relationship between GSTs gene polymorphism and T2DM in Europe, America, Africa and other geographic regions. The GSTM1- GSTT1, GSTM1- GSTP1 and GSTM1-GSTT1-GSTP1 gene polymorphisms combination are all associated with the increased susceptibility of T2DM. These results indicate that polygenes may have cumulative effects, but the potential gene-gene interactions still need to be studied with a larger sample size.

However, the present study used several subgroups and diverse genetic models, resulting in multiple comparisons, which necessitated adjusting the pooled p value (Attia et al., 2003). FPRP is regarded an acceptable tool for calculating the likelihood of positive findings in multiple hypothesis testing of molecular epidemiology investigations (Wacholder et al., 2004). Wakefield provided a much more reliable Bayesian metric of false positives in genetic epidemiology research (Wakefield, 2008). Potential errors and bias may be caused by genotyping errors or high inter-study heterogeneity (I2 > 50%) (Page et al., 2003; Clayton et al., 2005; Pompanon et al., 2005; Balding, 2006; Zhang et al., 2013), while large volumes evidence (statistical power >80%) have stricter statistical significance levels or lower false-discovery rate (Benjamini and Hochberg, 1995). Therefore, we consider using FPRP tests, BFDP tests and the Venice criteria to evaluate the false findings in present studies. The current meta-analysis found that all substantial relationship were deemed as less credible during did the reliability analysis. According to the results of the meta-regression analysis, the HWD studies and sources of controls were the main causes of heterogeneity. Bias and errors are prevalent in some low-quality, small-sample HWD studies, making the conclusions of these original studies undependable, molecular biology and disease susceptibility study, in particular. Furthermore, research with a small sample size and substantial results may be likely to swallow than studies with negative outcomes. When they produce positive results, however, their research might not have been rigorous and typically of poor quality.

Therefore, we assessed the sensitivity analysis using matching, high-quality, and HWE.

These three GSTs gene polymorphisms and their interactions were discovered to be strongly related to diabetes susceptibility in this investigation. From 2012 to 2019, a total of five related meta-studies were published. Recent meta-studies by Nath et al. (2019) showed that the GSTM1, GSTT1 null genotype and GSTM1-GSTT1 double null genotype increased the risk of T2DM in both Asian and Caucasian populations, the results are consistent with this study. An analysis of eleven studies by Yi et al. (2013) revealed that the GSTM1 gene polymorphisms increases the susceptibility to T2DM in Asians, Caucasians, and Africans. However, Saadat (2017) selected 18 studies (2,595 patients and 2,888 controls) and suggested that GSTP1 IIe105Val polymorphism was not connected with the risk of T2DM, which was contrary to the results of this study. Tang et al. (2013) also analyzed the three genes and suggested that GSTM1 and GSTT1 null genotypes were correlated to an increased risk of T2DM in Asians, while GSTP1 gene polymorphism was not. For the GSTP1 gene, the opposite result was obtained. Three related meta-analysis were published from 2015 to 2019. Nath et al. showed that GSTT1 null genotype increased the risk of T2DM related complications, but GSTM1 null genotype did not, which was consistent with this study. On the contrary, Orlewski and Orlewska (2015) found that GSTM1 gene polymorphisms increased the risk of diabetic nephropathy, while GSTT1 and GSTP1 gene polymorphisms were not statistically significant. Sun et al. (2015) suggested that a strong association GSTM1 and GSTP1 polymorphisms and diabetic retinopathy susceptibility. Further investigate the effect of gene-gene interaction on T2DM, two previous studies reported the correlation between GSTM1-GSTT1 gene polymorphism and T2DM susceptibility, but did not study the combined effect between the remaining genes. In this study, all except GSTT1-GSTP1 combination were found to be related. Moreover, these published meta-analyses are obviously inconsistent with current meta-analyses in terms of racial classification. The variability in results may be due to a variety of factors, such as small sample sizes, ethnic differences, failure to assess quality in most eligible studies and failure to develop a more complete genetic model.

There are several advantages to the newer meta-analysis. First, The FPRP, BFDP test, and Venice criterion were used to assess the credibility with positive outcomes. Secondly, quality assessment of qualified research is carried out. Third, the sample size was much larger and more detailed data and information were collected. Fourth, quality assessment of qualified research. Fifthly, a relatively complete genetic model is established. Sixth, more subgroup analyses were performed. Seventh, to our knowledge, this is the first meta-analysis to explore the combined effects of GSTM1-GSTP1, GSTT1-GSTP1, and GSTM1-GSTT1-GSTP1 gene polymorphisms on the risk of T2DM. There are some limitations of the present meta-analysis. To begin with, latest meta-analyses only include published articles, and positive results are frequently more likely to be published than negative consequences. The zero effect of GSTM1 may be overestimated if negative results are included. Second, because we could not calculate the HWE of these two genes, we did not evaluate whether the genotype distributions of GSTM1 and GSTT1 polymorphisms in the control group was in HWE. Third, heterogeneity between studies was much large in overall and among several subgroups. According to the results of the meta-regression analysis, the HWD studies and sources of controls were the main causes of heterogeneity, and adequate statistical correction and analysis did not significantly improve the heterogeneity. The source of heterogeneity may also be due to some other factors, such as study area, lifestyle and so on. Fourth, we conducted a subgroup analysis of race and region, and population heterogeneity has been resolved. However, T2DM is closely associated with the environment, the original study did not provide data on lifestyle, so we failed to solve the heterogeneity caused by lifestyle, and did not explore the effect of the gene-environment combination. Fifth, we failed to further investigate the effect of GSTs gene polymorphisms on the risk of T2DM related complications.

In conclusion, the current study revealed the relationship between GSTM1, GSTT1, and GSTP1 gene polymorphisms with the risk of T2DM, which may be attributed to false positive results rather than actual correlations or biological variables. Larger epidemiological studies should be carried out in the future to confirm or deny our findings.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Author contributions

L-sL: designed research, performed research, collected data, analyzed data, wrote paper. DW and RT: collected data. QW, LZ and JW: check the data. YL and X-fH: designed research and revised article.

Funding

Nature Science Foundation of China (81500648); Shanxi Basic Research Program (20210302123076; 20210302123012).

Acknowledgments

The authors of all of the original studies included in this meta-analysis are gratefully acknowledged.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Correction note

A correction has been made to this article. Details can be found at: 10.3389/fgene.2025.1749934.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.959291/full#supplementary-material

References

  • 1

    AbbasS.RazaS. T.MirS.SiddiqiZ.MahdiF. (2019). No association of SNP 313A→G in GSTP1 with nephropathy, hypertension and dyslipidemia in type 2 diabetes mellitus. Br. J. Biomed. Sci.76 (3), 153155. 10.1080/09674845.2019.1595870

  • 2

    AbbasiN.SalehiZ.AlizadehY. (2014). Genetic variation of GSTP1 in diabetic retinopathy [inPersian]. Arak Med. Uni J.17, 7076.

  • 3

    AfrandM.KhalilzadehS.BashardoostN.SheikhhaM. H. (2015). Evaluation of glutathione S-transferase T1 deletion polymorphism on type 2 diabetes mellitus risk in Zoroastrian females in Yazd, Iran. Indian J. Endocrinol. Metab.19 (1), 124128. 10.4103/2230-8210.146867

  • 4

    AhmedF. A. W.Al-BacharyH. S. (2017). Genetic polymorphism of glutathione S-transferase gene (GSTP1) in type 2 di-abetes mellitus patients in Basra province/Iraq. ijSciences.6 (5), 17. 10.18483/ijsci.1265

  • 5

    Al-BadranA. I.Al-MayahM. K. (2014). Association between GSTT1 and GSTM1 genes polymorphisms Type II diabetes miletus patients in Basra Iraq. Int. J. Curr. Microbiol. App Sci.3 (11), 288299.

  • 6

    AlbeladiF. I.MostafaM. M.ZayedM. A.AttaH. (2022). Association of polymorphisms in antioxidant enzyme-encoding genes with diabetic nephropathy in a group of Saudi arabian patients with type II diabetes mellitus. Int. J. Gen. Med.15, 59195928. 10.2147/IJGM.S367673

  • 7

    AmerM. A.GhattasM. H.Abo-ElmattyD. M.Abou-El-ElaS. H. (2011). Influence of glutathione S-transferase polymorphisms on type-2 diabetes mellitus risk. Genet. Mol. Res.10 (4), 37223730. 10.4238/2011.October.31.14

  • 8

    AmerM. A.GhattasM. H.Abo-ElmattyD. M.Abou-El-ElaS. H. (2012). Evaluation of glutathione S-transferase P1 genetic variants affecting type-2 diabetes susceptibility and glycemic control. Arch. Med. Sci.8, 631636. 10.5114/aoms.2012.30286

  • 9

    AttiaJ.ThakkinstianA.D'EsteC. (2003). Meta-analyses of molecular association studies: Methodologic lessons for genetic epidemiology. J. Clin. Epidemiol.56 (4), 297303. 10.1016/s0895-4356(03)00011-8

  • 10

    AzarovaI.BushuevaO.KonoplyaA.PolonikovA. (2018). Glutathione S-transferase genes and the risk of type 2 diabetes mellitus: Role of sexual dimorphism, gene-gene and gene-smoking interactions in disease susceptibility. J. Diabetes10 (5), 398407. 10.1111/1753-0407.12623

  • 11

    BakerW. L.WhiteC. M.CappelleriJ. C.KlugerJ.ColemanC. I.Health Outcomes, Policy, and Economics (HOPE) Collaborative Group (2009). Understanding heterogeneity in meta-analysis: The role of meta-regression. Int. J. Clin. Pract.63 (10), 14261434. 10.1111/j.1742-1241.2009.02168.x

  • 12

    BaldingD. J. (2006). A tutorial on statistical methods for population association studies. Nat. Rev. Genet.7 (10), 781791. 10.1038/nrg1916

  • 13

    BeggC. B.MazumdarM. (1994). Operating characteristics of a rank correlation test for publication bias. Biometrics50 (4), 10881101. 10.2307/2533446

  • 14

    BenjaminiY.HochbergY. (1995). Controlling the false discovery rate-a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B57, 289300. 10.1111/j.2517-6161.1995.tb02031.x

  • 15

    BidH. K.KonwarR.SaxenaM.ChaudhariP.AgrawalC. G.BanerjeeM. (2010). Association of glutathione S-transferase (GSTM1, T1 and P1) gene polymorphisms with type 2 diabetes mellitus in north Indian population. J. Postgrad. Med.56 (3), 176181. 10.4103/0022-3859.68633

  • 16

    CilenšekI.MankočS.PetrovičM. G.PetrovičD. (2012). GSTT1 null genotype is a risk factor for diabetic retinopathy in Caucasians with type 2 diabetes, whereas GSTM1 null genotype might confer protection against retinopathy. Dis. Markers32 (2), 9399. 10.3233/DMA-2011-0863

  • 17

    ClaytonD. G.WalkerN. M.SmythD. J.PaskR.CooperJ. D.MaierL. M.et al (2005). Population structure, differential bias and genomic control in a large-scale, case-control association study. Nat. Genet.37 (11), 12431246. 10.1038/ng1653

  • 18

    DadbinpourA.SheikhhaM. H.DarbouyM.Afkhami-ArdekaniM. (2013). Investigating GSTT1 and GSTM1 null genotype as the risk factor of diabetes type 2 retinopathy. J. Diabetes Metab. Disord.12 (1), 48. 10.1186/2251-6581-12-48

  • 19

    DattaS. K.KumarV.PathakR.TripathiA. K.AhmedR. S.KalraO. P.et al (2010). Association of glutathione S-transferase M1 and T1 gene polymorphism with oxidative stress in diabetic and nondiabetic chronic kidney disease. Ren. Fail.32 (10), 11891195. 10.3109/0886022X.2010.517348

  • 20

    de LimaR. M.Dos AnjosL. R. B.AlvesT. B.CoelhoA. S. G.PedrinoG. R.da Silva SantosR.et al (2018). Do GST polymorphisms influence in the pathogenesis of diabetic nephropathy?Mol. Cell. Endocrinol.478, 1016. 10.1016/j.mce.2018.07.001

  • 21

    DendupT.FengX.ClinganS.Astell-BurtT. (2018). Environmental risk factors for developing type 2 diabetes mellitus: A systematic review. Int. J. Environ. Res. Public Health15 (1), 78. 10.3390/ijerph15010078

  • 22

    DerSimonianR.LairdN. (2015). Meta-analysis in clinical trials revisited. Contemp. Clin. Trials45, 139145. 10.1016/j.cct.2015.09.002

  • 23

    DualS.TweedieR. (2000). A nonparametric “trim and fill” method of accounting for publication bias in meta-analysis. J. Am. Stat. Assoc.95, 8998. 10.2307/2669529

  • 24

    EggerM.Davey SmithG.SchneiderM.MinderC. (1997). Bias in meta-analysis detected by a simple, graphical test. BMJ315 (7109), 629634. 10.1136/bmj.315.7109.629

  • 25

    EtemadA.VasudevanR.AzizA. F.YusofA. K.KhazaeiS.FawziN.et al (2016). Analysis of selected glutathione S-transferase gene polymorphisms in Malaysian type 2 diabetes mellitus patients with and without cardiovascular disease. Genet. Mol. Res.15 (2), gmr5845. 10.4238/gmr.15025845

  • 26

    FujitaH.NaritaT.MeguroH.ShimoTomaiT.KitazatoH.KagayaE.et al (2000). No association of glutathione S-transferase M1 gene polymorphism with diabetic nephropathy in Japanese type 2 diabetic patients. Ren. Fail.22 (4), 479486. 10.1081/jdi-100100889

  • 27

    GonulN.KadiogluE.KocabasN. A.OzkayaM.KarakayaA. E.KarahalilB. (2012). The role of GSTM1, GSTT1, GSTP1, and OGG1 polymorphisms in type 2 diabetes mellitus risk: A case-control study in a Turkish population. Gene505 (1), 121127. 10.1016/j.gene.2012.05.025

  • 28

    GrubisaI.OtasevicP.DespotovicN.DedicV.MilašinJ.VucinicN. (2013). Genetic polymorphism of glutathion S-transferase P1 (GSTP1) Ile105Val and susceptibility to atherogenesis in patients with type 2 diabetes mellitus. Genetika45, 227236. 10.2298/gensr1301227g

  • 29

    GustiA. M. T.QustiS. Y.BahijriS. M.ToraihE. A.BokhariS.AttallahS. M.et al (2021). Glutathione S-transferase (GSTT1 rs17856199) and nitric oxide synthase (NOS2 rs2297518) genotype combination as potential oxidative stress-related molecular markers for type 2 diabetes mellitus. Diabetes Metab. Syndr. Obes.14, 13851403. 10.2147/DMSO.S300525

  • 30

    HayekT.StephensJ. W.HubbartC. S.AcharyaJ.CaslakeM. J.HaweE.et al (2006). A common variant in the glutathione S transferase gene is associated with elevated markers of inflammation and lipid peroxidation in subjects with diabetes mellitus. Atherosclerosis184 (2), 404412. 10.1016/j.atherosclerosis.2005.05.017

  • 31

    HayesJ. D.JuF.JowseyI. R. (2005). Glutathione transferases. Annu. Rev. Pharmacol. Toxicol.45, 5188. 10.1146/annurev.pharmtox.45.120403.095857

  • 32

    HayesJ. D.StrangeR. C. (2000). Glutathione S-transferase polymorphisms and their biological consequences. Pharmacology61 (3), 154166. 10.1159/000028396

  • 33

    HigginsJ. P.ThompsonS. G.DeeksJ. J.AltmanD. G. (2003). Measuring inconsistency in meta-analyses. BMJ327 (7414), 557560. 10.1136/bmj.327.7414.557

  • 34

    HoriM.OnikiK.UedaK.GotoS.MiharaS.MarubayashiT.et al (2007). Combined glutathione S-transferase T1 and M1 positive genotypes afford protection against type 2 diabetes in Japanese. Pharmacogenomics8 (10), 13071314. 10.2217/14622416.8.10.1307

  • 35

    IoannidisJ. P.BoffettaP.LittleJ.O'BrienT. R.UitterlindenA. G.VineisP.et al (2008). Assessment of cumulative evidence on genetic associations: Interim guidelines. Int. J. Epidemiol.37 (1), 120132. 10.1093/ije/dym159

  • 36

    JamilH.AwanA.AkbarA.BabarM.AkhtarS.IqbalR. K.et al (2022). A study of association between presence or absence of GSTT1 and GSTM1 and/or single nucleotide polymorphism in FABP2 and GSTP1 with incidence of diabetes type 2: A case-control study. J. Pak. Med. Assoc.72 (4), 714720. 10.47391/JPMA.1337

  • 37

    JanaM.PetrovicD. (2012). Diabetic nephropathy in type 2 diabetes: MPO T-764C genotype is associated with oxidative stress. Open Life Sci.7 (6), 964972. 10.2478/s11535-012-0105-5

  • 38

    KlusekJ.Błońska-SikoraE.WitczakB.OrlewskaK.GluszekS.OrlewskaE. (2020). Glutathione S-transferases gene polymorphism influence on the age of diabetes type 2 onset. BMJ Open Diabetes Res. Care8 (2), e001773. 10.1136/bmjdrc-2020-001773

  • 39

    MastanaS. S.KaurA.HaleR.LindleyM. R. (2013). Influence of glutathione S-transferase polymorphisms (GSTT1, GSTM1, GSTP1) on type-2 diabetes mellitus (T2D) risk in an endogamous population from north India. Mol. Biol. Rep.40 (12), 71037110. 10.1007/s11033-013-2833-7

  • 40

    MerganiA.MansourA. A.AskarT.ZahranR. N.Mus-tafaA. M.MohammedM. A.et al (2016). Glutathione S-trans-ferase pi-ile 105 val polymorphism and susceptibility to T2DM in population from turabah region of Saudi arabia. Biochem. Genet.54, 544551. 10.1007/s10528-016-9740-2

  • 41

    MirR.BhatM. A.JavaidJ.Shah,N.,KumarP.SharmaE.et al (2016). Glutathione S-transferase M1 and T1 (rs4025935 and rs71748309) null genotypes are associated with increased susceptibility to coronary artery disease in Indian populations. Acta Cardiol.71 (6), 678684. 10.2143/AC.71.6.3178186

  • 42

    MoasserE.AzarpiraN.ShiraziB.SaadatM.GeramizadehB. (2014). Genetic polymorphisms of glutathione-s-transferase M1 and T1 genes with risk of diabetic retinopathy in Iranian population. Iran. J. Basic Med. Sci.17 (5), 351356.

  • 43

    MoasserE.Kazemi-NezhadS. R.SaadatM.AzarpiraN. (2012). Study of the association between glutathione S-transferase (GSTM1, GSTT1, GSTP1) polymorphisms with type II diabetes mellitus in southern of Iran. Mol. Biol. Rep.39 (12), 1018710192. 10.1007/s11033-012-1893-4

  • 44

    NathS.DasS.BhowmikA.GhoshS. K.ChoudhuryY. (2019). The GSTM1 and GSTT1 null genotypes increase the risk for type 2 diabetes mellitus and the subsequent development of diabetic complications: A meta-analysis. Curr. Diabetes Rev.15 (1), 3143. 10.2174/1573399814666171215120228

  • 45

    NowierS. R.KashmiryN. K.RasoolH. A.MoradH.IsmailS. (2009). Association of Type 2 diabetes mellitus and glutathione S-transferase (GSTM1 and GSTT1) genetic polymorphism. Res. J. Med. Med. Sci.4 (2), 181188.

  • 46

    OnikiK.UmemotoY.NagataR.HoriM.MiharaS.MarubayashiT.et al (2008). Glutathione S-transferase A1 polymorphism as a risk factor for smoking-related type 2 diabetes among Japanese. Toxicol. Lett.178, 143145. 10.1016/j.toxlet.2008.03.004

  • 47

    OrlewskiJ.OrlewskaE. (2015). Effects of genetic polymorphisms of glutathione S-transferase genes (GSTM1, GSTT1, GSTP1) on the risk of diabetic nephropathy: A meta-analysis. Pol. Arch. Med. Wewn.125 (9), 649658. 10.20452/pamw.3045

  • 48

    OsmanH.OsmanM. Y.AlsaidyA. A. (2019). GSTT1, Calpain 10 SNP 19 and indices of glycaemia in type 2 diabetes. Br. J. Biomed. Sci.76 (4), 205207. 10.1080/09674845.2019.1660501

  • 49

    PageG. P.GeorgeV.GoR. C.PageP. Z.AllisonD. B. (2003). Are we there yet?”: Deciding when one has demonstrated specific genetic causation in complex diseases and quantitative traits. Am. J. Hum. Genet.73 (4), 711719. 10.1086/378900

  • 50

    PahwaS.SharmaR.SinghB. (2017). Role of glutathione S-transferase in coronary artery disease patients with and without type 2 diabetes mellitus. J. Clin. Diagn. Res.11 (1), BC05BC08. 10.7860/JCDR/2017/23846.9281

  • 51

    PinheiroD. S.RochaF. C. R.MundimC. A.JuniorP. d. M.UlhoaC. J.ReisA. A. S.et al (2013). Evaluation of glutathione S-transferase GSTM1 and GSTT1 deletion polymorphisms on type-2 diabetes mellitus risk. PLoS One8 (10), e76262. 10.1371/journal.pone.0076262

  • 52

    PompanonF.BoninA.BellemainE.TaberletP. (2005). Genotyping errors: Causes, consequences and solutions. Nat. Rev. Genet.6 (11), 847859. 10.1038/nrg1707

  • 53

    PourkeramatiA.Zare MehrjardiE.Dehghan TezerjaniM.SeifatiS. M. (2020). Association of GSTP1, GSTT1 and GSTM1 gene variants with coronary artery disease in Iranian population: A case-control study. Int. J. Gen. Med.13, 249259. 10.2147/IJGM.S252552

  • 54

    PurkaitP.HalderK.RoyA. (2014). GSTM1 null genotype associated with type 2 diabetic nephropathy patients among Indian population. World J. Pharm. Res.3, 44524463.

  • 55

    RamprasathT.Senthil MuruganP.PrabakaranA. D.GomathiP.RathinavelA.SelvamG. S. (2011). Potential risk modifications of GSTT1, GSTM1 and GSTP1 (glutathione-S-transferases) variants and their association to CAD in patients with type-2 diabetes. Biochem. Biophys. Res. Commun.407 (1), 4953. 10.1016/j.bbrc.2011.02.097

  • 56

    RaoD. K.ShaikN. A.ImranA.MurthyD. K.GantiE.ChintaC.et al (2014). Variations in the GST activity are associated with single and combinations of GST genotypes in both male and female diabetic patients. Mol. Biol. Rep.41 (2), 841848. 10.1007/s11033-013-2924-5

  • 57

    RasheedM. N.HasanO. M.MahmoodA. S. (2015). Association of glutathione S-transferase (GSTM1, T1) gene polymorphisms with type 2 diabetes mellitus (T2DM) in the Iraqi patients. Iraqi J. Biotechnol.14 (1), 7076.

  • 58

    RasheedR. H.Al-EssaN. E.ShihabB. A. (2016). Association of glutathione S-transferase (GSTP1) genetic polymor-phism in Iraqi patients with diabetes mellitus type 2. Baghdad Sci. J.13, 3643. 10.21123/bsj.13.1.36-43

  • 59

    RazaS.AbbasS.AhmadA.AhmedF.ZaidiZ.MahdiF. (2014). Association of Glutathione-S-Transferase (GSTM1 and GSTT1) and FTO gene polymorphisms with type 2 diabetes mellitus cases in northern India. Balk. J. Med. Genet.17 (1), 4754. 10.2478/bjmg-2014-0027

  • 60

    SaadatM. (2017). Evaluation of glutathione S-transferase P1 (GSTP1) Ile105Val polymorphism and susceptibility to type 2 diabetes mellitus, a meta-analysis. EXCLI J.16, 11881197. 10.17179/excli2017-828

  • 61

    StoianA.BanescuC.BalasaR. I.MotataianuA.StoianM.MoldovanV. G.et al (2015). Influence of GSTM1, GSTT1, and GSTP1 polymorphisms on type 2 diabetes mellitus and diabetic sensorimotor peripheral neuropathy risk. Dis. Markers2015, 638693. 10.1155/2015/638693

  • 62

    SunL.ZhangY.XiongY. (2015). GSTM1 and GSTT1 null genotype and diabetic retinopathy: A meta-analysis. Int. J. Clin. Exp. Med.8 (2), 16771683.

  • 63

    Tabatabaei-MalazyO.KhodaeianM.BitarafanF.LarijaniB.AmoliM. (2017). Polymorphisms of antioxidant genes as a target for diabetes management. Int. J. Mol. Cell. Med.6 (3), 135147. 10.22088/acadpub.BUMS.6.3.135

  • 64

    TangS. T.WangC. J.TangH. Q.ZhangQ.WangY. (2013). Evaluation of glutathione S-transferase genetic variants affecting type 2 diabetes susceptibility: A meta- analysis. Gene530 (2), 301308. 10.1016/j.gene.2013.08.043

  • 65

    TangvarasittichaiS. (2015). Oxidative stress, insulin resistance, dyslipidemia and type 2 diabetes mellitus. World J. Diabetes6 (3), 456480. 10.4239/wjd.v6.i3.456

  • 66

    TheodoratouE.MontazeriZ.HawkenS.AllumG. C.GongJ.TaitV.et al (2012). Systematic meta-analyses and field synopsis of genetic association studies in colorectal cancer. J. Natl. Cancer Inst.104 (19), 14331457. 10.1093/jnci/djs369

  • 67

    TiwariA. K.PrasadP.KumarK. M. P.AmminiA. C.GuptaA.et al (2009). Oxidative stress pathway genes and chronic renal insufficiency in Asian Indians with Type 2 diabetes. J. Diabetes Complicat.23 (2), 102111. 10.1016/j.jdiacomp.2007.10.003

  • 68

    TsaiJ. P.YangS. F.WuS. W.HungT. W.TsaiH. C.LianJ. D.et al (2011). Glutathione S-transferase gene polymorphisms are not major risks for susceptibility to posttransplantation diabetes mellitus in Taiwan renal transplant recipients. J. Clin. Lab. Anal.25 (6), 432435. 10.1002/jcla.20498

  • 69

    VatsP.ChandraH.BanerjeeM. (2013). Glutathione S-transferase and Catalase gene polymorphisms with Type 2 diabetes mellitus. Dis. Mol. Med.1 (3), 4653. 10.5455/dmm.20131027101207

  • 70

    VatsP.SagarN.SinghT. P.BanerjeeM. (2015). Association of Superoxide dismutases (SOD1 and SOD2) and Glutathione peroxidase 1 (GPx1) gene polymorphisms with type 2 diabetes mellitus. Free Radic. Res.49 (1), 1724. 10.3109/10715762.2014.971782

  • 71

    WacholderS.ChanockS.Garcia-ClosasM.El GhormliL.RothmanN. (2004). Assessing the probability that a positive report is false: An approach for molecular epidemiology studies. J. Natl. Cancer Inst.96 (6), 434442. 10.1093/jnci/djh075

  • 72

    WakefieldJ. (2008). A Bayesian measure of the probability of false discovery in genetic epidemiology studies. Am. J. Hum. Genet.81 (2), 208227. 10.1086/519024.1

  • 73

    WangG.ZhangL.LiQ. (2006). Genetic polymorphisms of GSTT1, GSTM1, and NQO1 genes and diabetes mellitus risk in Chinese population. Biochem. Biophys. Res. Commun.341 (2), 310313. 10.1016/j.bbrc.2005.12.195

  • 74

    WuY.DingY.TanakaY.ZhangW. (2014). Risk factors contributing to type 2 diabetes and recent advances in the treatment and prevention. Int. J. Med. Sci.11 (11), 11851200. 10.7150/ijms.10001

  • 75

    YalinS.HatungilR.TamerL.AtesN. A.DogruerN.YildirimH.et al (2007). Glutathione S-transferase gene polymorphisms in Turkish patients with diabetes mellitus. Cell biochem. Funct.25 (5), 509513. 10.1002/cbf.1339

  • 76

    YangY.KaoM. T.ChangC. C.ChungS. Y.ChenC. M.TsaiJ. J. P.et al (2004). Glutathione S-transferase T1 deletion is a risk factor for developing end-stage renal disease in diabetic patients. Int. J. Mol. Med.14 (5), 855859. 10.3892/ijmm.14.5.855

  • 77

    YiR.LiuB.DongQ. (2013). Assessment of the association between GSTM1 null genotype and risk of type 2 diabetes. Mol. Biol. Rep.40 (6), 39613967. 10.1007/s11033-012-2473-3

  • 78

    ZakiM. A.MoghazyT. F.El-DeebM. M. K.MohamedA. H.MohamedN. A. A. (2015). Glutathione S-transferase M1, T1 and P1 gene polymorphisms and the risk of developing type 2 diabetes mellitus in Egyptian diabetic patients with and without diabetic vascular complications. Alexandria J. Med.51 (1), 7382. 10.1016/j.ajme.2014.03.003

  • 79

    ZhangJ.LiuH.YanH.HuangG.WangB. (2013). Null genotypes of GSTM1 and GSTT1 contribute to increased risk of diabetes mellitus: A meta-analysis. Gene518 (2), 405411. 10.1016/j.gene.2012.12.086

  • 80

    ZimniakP.NanduriB.PikułaS.Bandorowicz-PikulaJ.SinghalS. S.SrivastavaS. K.et al (1994). Naturally occurring human glutathione S-transferase GSTP1-1 isoforms with isoleucine and valine in position 104 differ in enzymic properties. Eur. J. Biochem.224 (3), 893899. 10.1111/j.1432-1033.1994.00893.x

Summary

Keywords

meta-analysis, genetic polymorphism, GSTM1, GSTP1, GSTT1, T2DM

Citation

Liu L-s, Wang D, Tang R, Wang Q, Zheng L, Wei J, Li Y and He X-f (2022) Individual and combined effects of the GSTM1, GSTT1, and GSTP1 polymorphisms on type 2 diabetes mellitus risk: A systematic review and meta-analysis. Front. Genet. 13:959291. doi: 10.3389/fgene.2022.959291

Received

01 June 2022

Accepted

14 October 2022

Published

07 November 2022

Corrected

01 December 2025

Volume

13 - 2022

Edited by

Claudia Banescu, University of Medicine, Pharmacy, Sciences and Technology of TârguMureş, Romania

Reviewed by

Stoian Adina, George Emil Palade University of Medicine, Pharmacy, Sciences and Technology of Târgu Mureş, Romania

Florin Tripon, George Emil Palade University of Medicine, Pharmacy, Sciences and Technology of Târgu Mureş, Romania

Updates

Copyright

*Correspondence: Yan Li, ; Xiao-feng He,

This article was submitted to Genetics of Common and Rare Diseases, a section of the journal Frontiers in Genetics

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics