2018 Volume 246 Issue 1 Pages 35-44
Glutathione S-transferases (GSTs), a superfamily of multifunctional enzymes, play an important role in the onset and progression of renal cell carcinoma (RCC). However, novel GST omega class (GSTO), consisting of GSTO1-1 and GSTO2-2 isoenzymes, has not been studied in RCC yet. Two coding single nucleotide polymorphisms (SNPs) supposedly affect their functions: GSTO1*C419A (rs4925) causing alanine to aspartate substitution (*A140D) and GSTO2*A424G (rs156697) causing asparagine to aspartate substitution (*N142D), and have been associated with several neurodegenerative diseases and cancers. Functional relevance of yet another GSTO2 polymorphism, identified at the 5′ untranslated (5′UTR) gene region (GSTO2*A183G, rs2297235), has not been clearly discerned so far. Therefore, we aimed to assess the effect of specific GSTO1 and GSTO2 gene variants, independently and in interaction with established risk factors (smoking, obesity and hypertension) on the risk for the most aggressive RCC subtype, the clear cell RCC (ccRCC). Genotyping was performed in 239 ccRCC patients and 350 matched controls, while plasma levels of 8-hydroxy-2′-deoxyguanosine (8-OHdG), a biomarker of oxidative DNA damage, were determined by ELISA. As a result, combined effect of all three variant genotypes exhibited almost 3-fold risk of RCC development. Additionally, this association was confirmed at the haplotype level [variant GSTO1*A/GSTO2*G (rs156697)/GSTO2*G (rs2297235) haplotype], suggesting a potential role of those variants in propensity to RCC. Regarding the gene-environment interactions, variant GSTO2*G (rs156697) homozygous smokers are at higher ccRCC risk. Association in terms of oxidative DNA damage was found for GSTO2 polymorphism in 5′UTR and 8-OHdG. In conclusion, the concomitance of GSTO polymorphisms may influence ccRCC risk.
Glutathione S-transferases (GSTs), a large family of multifunctional enzymes, play a significant role in the onset and progression of renal cell carcinoma (RCC) (Pljesa-Ercegovac et al. 2008; Searchfield et al. 2011). Cytosolic GSTs are involved in detoxification of electrophilic metabolites of endogenous or exogenous origin by conjugation with glutathione (Strange et al. 2001). A vast number of carcinogens implicated in the occurrence of RCC, such as polycyclic aromatic hydrocarbon diol-epoxides and halogenated solvents (Sweeney et al. 2000), as well as chemotherapeutic agents are metabolized by GSTs (Wu and Dong 2012). Recent studies suggest that the possible mechanisms underlying the association of GSTs with cancer development and progression might be contributed to their role in cellular redox regulation (Tew and Townsend 2012). It is noteworthy to mention that concerning redox status, RCC belongs to a group of solid tumors with a certain level of oxidative distress. Specifically, prooxidative environment might be important in the early phase of RCC development, while in the course of cancer progression, a shift towards more reduced state occurs (Lusini et al. 2001). Lately, other biologically important roles of these enzymes have been specified, including protein-protein interactions. Namely, GSTs act as modulators of the mitogen-activated protein kinase (MAPK) signaling pathway involved in cellular survival and apoptosis (Tew and Townsend 2012). Furthermore, polymorphisms in GST genes could affect both individual response to carcinogen exposure and pharmacogenomic-based cancer treatment (Lo and Ali-Osman 2007). So far the association of common GST gene polymorphisms (GSTM1, GSTT1, GSTP1 and GSTA1), independently or in interaction with environmental factors, with increased risk of developing RCC has been shown (Coric et al. 2016; Zhong et al. 2018). Novel omega class GST consists of two isoenzymes, GSTO1-1 and GSTO2-2. In contrast to all other GSTs, they exhibit the whole range of specific activities, including a thioltransferase, dehydroascorbate reductase and monomethylarsonat reductase activity (Board and Menon 2016). In this line, both enzymes are considered to be involved in the regulation of cellular redox balance: GSTO1-1 by its deglutathionylation activity and GSTO2-2 as the enzyme with the highest dehydroascorbate-reductase (DHAR) activity that is responsible for maintaining high cellular vitamin C levels. GSTO1-1 also possesses numerous regulatory roles, such as modulation of ryanodine receptors, posttranslational processing of interleukin-1β (IL1-β) and proposed anti-apoptotic role (Dulhunty et al. 2001; Laliberte et al. 2003; Piaggi et al. 2010). Furthermore, it has been suggested that GSTO1-1 could also contribute to tumor chemoresistance. As shown for cisplatin-resistant sublines of ovarian cancer and HeLa cells, increased GSTO1-1 expression was associated with the activation of survival signaling pathways (Akt and ERK1/2) and inhibition of apoptotic MAPK pathway (JNK1) (Yan et al. 2007; Piaggi et al. 2010). Numerous studies which investigated apoptosis and survival pathways in terms of RCC (Banumathy and Cairns 2010) have confirmed the therapeutic potential of phosphatidylinositol-3-Kinase (PI3K)/Akt inhibitors (Lin et al. 2006; Zhong et al. 2016). In that context, it seems that this novel class of GST enzymes could contribute to RCC progression, affecting both the proliferation capacity of tumor cells and response to therapy.
Two GSTO actively transcribed genes (GSTO1 and GSTO2) were identified in the human population (Fig. 1), located 1.5 kb apart on the long arm of chromosome 10 (10q25.1) (Board and Menon 2016). Mukherjee et al. (2006) decribed a total of 31 polymorphisms in GSTO1 and 66 polymorphisms of GSTO2 gene, two of which are commonly studied: GSTO1*C419A (rs4925) causing alanine to aspartate substitution in amino acid 140 (*A140D) and GSTO2*A424G (rs156697) which causes an asparagine to aspartate substitution in amino acid 142 (*N142D). Moreover, in recent years several studies investigated another polymorphism of GSTO2, identified within the 5′ untranslated (5′ UTR) gene region (GSTO2*A183G, rs2297235) (Wang et al. 2009). Strong linkage disequilibrium has been demonstrated for all three mentioned polymorphisms (Wang et al. 2009). Moreover, it has been shown that the most investigated GSTO1 polymorphism leads to change in aforementioned deglutathionylase and thioltransferase activity, however, without affecting their monomethylarsonat reductase activity (Tanaka-Kagawa et al. 2003; Whitbread et al. 2003; Menon and Board 2013). Regarding GSTO2 rs156697 polymorphism, a strong association was indicated between variant GSTO2*G allele and its protein levels (Mukherjee et al. 2006; Allen et al. 2012). Therefore, it came as no surprise that polymorphisms of GST omega class members have been investigated in regard to a number of clinical disorders, including Alzheimer’s disease, Parkinson’s disease, vascular dementia and stroke, amyotrophic lateral sclerosis, chronic obstructive pulmonary disease (Board and Menon 2016) and cancer (Xu et al. 2014). Although the particular association of these polymorphisms with several neurodegenerative diseases might be explained by modulating role of GSTO1-1 and GSTO2-2 in redox balance regulation (Allen et al. 2012; Board and Menon 2016), their functional relevance in cancer has not been clearly discerned (Xu et al. 2014).
Due to observed expression of GST omega class enzymes in kidney parenchyma (Whitbread et al. 2005), we hypothesized that GSTO1 and GSTO2 polymorphisms could either alone or in association with known risk factors (smoking, obesity and hypertension) influence the risk of clear cell RCC (ccRCC), the most aggressive RCC subtype. Therefore, in this study, we for the first time evaluated the modifying effect of three specific GSTO1 and GSTO2 gene variants (GSTO1*C419A rs4925, GSTO2*A424G rs156697 and GSTO2*A183G rs2297235), independently and in interaction with known risk factors, on propensity to develop ccRCC. To obtain a clear insight into the functional significance of these findings, we also studied the association of distinct GSTO1 and GSTO2 polymorphisms with the systemic level of oxidative DNA damage, determined as the level of 8-OHdG.

Schematic representation of the investigated polymorphisms in GSTO1 and GSTO2 genes.
GSTO1 and GSTO2 genes are located 1.5 kb apart on the long arm of chromosome 10 (10q25.1). Exons that encode the open read frame are represented as black rectangles, whereas white rectangles represent untranslated regions (UTRs). The SNPs are indicated as position with the allele variants followed by the rs number. These polymorphisms are in strong linkage disequilibrium.
We enrolled 239 patients (162 men, 77 women; average age 58.94 ± 11.64 years) from the Urology Clinic, Clinical Center of Serbia, Belgrade, diagnosed with ccRCC that was histologically confirmed by pathologists specialized in uropathology. The control group comprised 350 sex- and age-matched individuals (217 men, 133 women; average age 60.16 ± 11.11 years) recruited at the same clinical center, with no previous personal history of cancer. Upon obtaining informed consent, each subject was interviewed using a standard questionnaire, composed at the Institute of Epidemiology, University of Belgrade Faculty of Medicine (UBFM). The information regarding demographic characteristics, as well as known RCC risk factors (smoking history, hypertension and obesity) were collected.
Ethical approvalAll procedures performed in the study were in accordance with the ethical standards of the 1964 Declaration of Helsinki revised in 2013. The study was approved by the Institutional Ethical board (October13, 2011, approval number 29/X-3, Faculty of Medicine, University of Belgrade, Serbia).
Informed consentInformed consent was obtained from all individual participants included in the study.
GST genotypingDNA from the whole blood was isolated using QIAamp DNA mini kit (Qiagen,USA). GSTO1*C419A (rs4925), GSTO2*A424G (rs156697) and GSTO2*A183G (rs2297235) polymorphisms were determined by quantitative polymerase chain reaction (qPCR), performed on Mastercycler ep realplex (Eppendorf, Germany) using Applied Biosystems TaqMan SNP Genotyping assay (Thermo Fisher Scientific, USA, assay ID: C_11309430_30, C_3223136_1 and C_3223142_1, respectively).
Determination of 8-hydroxy-2′-deoxyguanosine (8-OHdG)The quantification of 8-OHdG (ng/ml) in plasma samples of ccRCC patients and controls was performed using the OxiSelect Oxidative DNA Damage ELISA Kit (Cell Biolabs, Inc., USA).
Statistical analysisStatistical analysis was performed using the SPSS software version 17.0 (SPSS Inc., Chicago, IL, USA). Selected characteristics among ccRCC patients and controls were compared using the χ2 test for categorical variables, whereas depending on data distribution, the Student’s t test or Mann-Whitney ranksum test was used for continuous variables. Distribution was tested by using graphical methods, as well as Kolmogorov-Smirnov and Shapiro-Wilk tests. χ2 test was used to test deviation of the genotype distribution from Hardy-Weinberg equilibrium. Estimation of linkage disequilibrium (LD) between pairs of SNPs was examined using the SNPStats (Solé et al. 2006). The strength of LD was expressed in terms of D′ = D/Dmax. Effects of genetic variants on ccRCC risk were computed by odds ratios (OR) and 95% confidence intervals (CI) using logistic regression analysis adjusted for age, sex, smoking, hypertension and obesity.
Baseline characteristics of 239 patients with ccRCC and 350 controls are shown in Table 1. There was no significant difference between patients and controls with regards to age, sex, smoking and obesity (P > 0.05). We found that hypertensive subjects were at 3.54-fold higher risk of developing ccRCC compared to normotensive subjects (95% CI: 2.35-5.32, P < 0.001).
Distribution of the GSTO1 (rs4925) and GSTO2 (rs156697, rs2297235) genotypes in patients and controls is presented in Table 2. As shown, carriers of GSTO1*A/A and GSTO2*G/G (rs156697) variant genotypes were at higher risk of ccRCC development when compared to referent genotypes, but this association did not reach statistical significance (OR = 1.35, 95% CI: 0.70-2.61, P = 0.364 and OR = 1.78, 95% CI: 0.91-3.50, P = 0.092, respectively). However, when GSTO1*C419A and GSTO2*A424G polymorphisms were analyzed in combination, the significant association was obtained. Namely, individuals with combined variant GSTO1*A/A and GSTO2*G/G (rs156697) genotypes had 2.6-fold higher risk of developing cancer compared to referent genotype combination (95% CI: 1.09-6.19, P = 0.031). Combined effect of all three polymorphisms showed no further increase in the risk of ccRCC development (Table 2). It is important to note that our results on combined effects of GSTO polymorphisms were also confirmed in haplotype analysis. Namely, since both GSTO1 and GSTO2 genes are located on the same chromosome, just 1.5 kb apart, we performed a linkage disequilibrium analysis. In this LD analysis, we found D′ of 0.64 for GSTO1*C419A and GSTO2*A424G (P < 0.001), 0.83 for GSTO1*C419A and GSTO2*A183G (P < 0.001), and 0.80 for GSTO2*A424G and GSTO2*A183G (P < 0.001), indicating high LD between these SNPs. The most prevalent haplotype among controls (52%) and patients (56%) is H1, consisting of GSTO1*C, GSTO2*A (rs156697) and GSTO2*A (rs2297235) wild-type alleles. Carriers of H2 haplotype, represented by one copy of variant GSTO1*A, GSTO2*G (rs156697) and GSTO2*G (rs2297235), had significantly increased risk of 1.46 for ccRCC (95% CI: 1.02-2.09, P = 0.041) (Table 3).
We also investigated possible modifying effect of GSTO genotypes in conjunction with hypertension, obesity and smoking, as recognized risk factors, on ccRCC risk. We found no modifying effect of GSTO genotypes with hypertension (Table 4) and obesity (Table 5), confirming the significance of hypertension as ccRCC independent risk factor, regardless of GSTO genotypes. However, significant modifying effect on ccRCC risk conferred by smoking has been found only in GSTO2*G/G (rs156697) carriers (OR = 2.44, 95% CI: 1.04-5.71, P = 0.040) (Table 6), while another two studied polymorphisms (GSTO1 rs4925 and GSTO2 rs2297235) did not show any association with smoking.
Regarding the degree of oxidative DNA damage in plasma samples, 8-OHdG levels were higher in ccRCC patients compared to controls (1.00 ng/ml vs. 0.70 ng/ml, respectively). In attempt to discern functional role of GSTO1 and GSTO2 polymorphisms, we assessed the degree of oxidative DNA damage in ccRCC patients stratified according to GSTO genotype. The results obtained showed that levels of 8-OHdG, a biomarker of oxidative DNA damage, were significantly higher in patients with GSTO2*G/G variant genotype (rs2297235) than in carriers of at least one GSTO2*A referent allele (1.41 ng/ml vs. 0.99 ng/ml, respectively P = 0.042). No association was found with plasma 8-OHdG levels and GSTO1 or another GSTO2 (rs156697) genotypes (Table 7).

Baseline characteristics of ccRCC patients and controls.
amean ± SD.
bOR, odds ratio adjusted to age, smoking status, hypertension, obesity.
cEvery day smoking during a minimum of 60-day period prior to the study onset.
dOR, odds ratio adjusted to age, sex, hypertension, obesity.
eMedian (min-max).
fOR, odds ratio adjusted to age, sex, smoking status, obesity.
gBMI, body mass index; Obese participants were defined as individuals with BMI above 30.
hOR, odds ratio adjusted to age, sex, smoking status, hypertension; CI, confidence interval.
iAvailable data on patients’ tumor grade and stage, depending on the type of surgery and pathohistology diagnostics.
P < 0.05 was considered to be statistically significant.

Distribution of individual GSTO1 (rs4925) and GSTO2 (rs156697 and rs2297235), as well as combined GSTO1/GSTO2 genotypes in ccRCC patients and controls.
aOR, odds ratio adjusted to age, sex, smoking status, hypertension, obesity; CI, confidence interval.
bFor GSTO1 rs4925, genotyping was successful in 236 of 239 patients and 338 of 350 controls.
cFor GSTO2 rs156697, genotyping was successful in all recruited patients and 333 of 350 controls.
dFor GSTO2 rs2297235, genotyping was successful in 231 of 239 patients and 338 of 350 controls.
P < 0.05 was considered to be statistically significant.

Haplotypes of GSTO1 (rs4925) and GSTO2 (rs156697 and rs2297235) in relation to the risk of ccRCC.
Global haplotype association P-value: < 0.0001.
aOR, odds ratio adjusted to age, sex, smoking status, hypertension, obesity; CI, confidence interval.
P < 0.05 was considered to be statistically significant.

Distribution of GSTO1 (rs4925) and GSTO2 (rs156697 and rs2297235) genotypes in relation to hypertension in ccRCC patients and controls.
aOR, odds ratio adjusted to age, sex, smoking status, obesity; CI, confidence interval.
P < 0.05 was considered to be statistically significant.

Distribution of GSTO1 (rs4925) and GSTO2 (rs156697 and rs2297235) genotypes in relation to obesity in ccRCC patients and controls.
aOR, odds ratio adjusted to age, sex, smoking status, hypertension; CI, confidence interval.
bObese participants were defined as individuals with BMI (body mass index) above 30.
P < 0.05 was considered to be statistically significant.

Distribution of GSTO1 (rs4925) and GSTO2 (rs156697 and rs2297235) genotypes in relation to smoking status in ccRCC patients and controls.
aOR, odds ratio adjusted to age, sex, hypertension, obesity; CI, confidence interval.
bSmoking status was categorized into non-smokers and smokers with respect to the limit of a minimum of 60-day period of every day smoking prior to their enrollment in the study.
P < 0.05 was considered to be statistically significant.

The association between GSTO genotypes and the levels of 8-OHdG in ccRCC patients.
aMedian (min-max).
P < 0.05 was considered to be statistically significant.
The results of this study have shown that individuals with combined variant GSTO1 and GSTO2 genotypes exhibit higher risk of ccRCC development compared to the carriers of both referent genotypes. Regarding the gene-environment interactions and RCC risk, variant GSTO2 homozygous smokers are at higher risk when compared to non-smokers with at least one referent allele. Association in terms of oxidative phenotype was found only for GSTO2 polymorphism in 5′UTR and 8-OHdG.
Until now, significant association between GSTO1*C419A polymorphism (rs4925) and risk of acute lymphoblastic leukemia, hepatocellular, breast, bile duct and non-small cell lung cancer has been reported (Xu et al. 2014), whereas GSTO2*A424G polymorphism (rs156697) increased the risk of ovarian (Pongstaporn et al. 2006), breast (Xu et al. 2014) and bladder cancer (Djukic et al. 2015). In line with these findings, we also showed that GSTO2*A424G polymorphism (rs156697) increases the risk of ccRCC, but the effect was statistically significant when it was analyzed in combination with GSTO1 polymorphism or smoking. Considering the observed significant LD of these genes, we also focused on evaluation of possible effect of haplotypes on ccRCC risk. We found that the carriers of H2 haplotype, comprising GSTO1*A (rs4925), GSTO2*G (rs156697) and GSTO2*G (rs2297235) allelic variants, are at increased risk for ccRCC.
In the course of RCC development, complex changes of cellular redox regulation contribute to survival of tumor cells and disease progression (Pljesa-Ercegovac et al. 2008). Based on specific set of antioxidant and regulatory activities of the GSTO enzymes, the polymorphisms of GSTO1 and GSTO2 might influence inter-individual susceptibility to oxidative stress. It has been shown that GSTO1*C allelic variant has higher deglutathionylation and thioltransferase activity in contrast to *A variant (Tanaka-Kagawa et al. 2003; Menon and Board 2013). Since glutathionylation, as post-translational modification, can influence the activity of many proteins involved in tumor growth, the presence of GSTO1 allelic variants with altered activity could provide a plausible mechanism to explain the associations between this genetic polymorphism and different cancers (Xu et al. 2014). In addition to GSTO1, it seems that polymorphisms in GSTO2 could also affect primarily its antioxidant dehydroascorbate reductase (DHAR) activity (Whitbread et al. 2005; Piacentini et al. 2013), which is responsible for regeneration of ascorbic acid (vitamin C). Moreover, increased uptake of dehydroascorbate via the GLUT1 transporters has been shown for some colorectal tumor cell lines (Yun et al. 2015). Since increased expression of GLUT1 transporters is common feature of different solid tumors, including RCC, this phenomenon of higher dehydroascorbate uptake needs further elucidation. Nevertheless, low DHAR activity in subjects with both variant GSTO2 alleles might result in deficient vitamin C regeneration and accumulation of oxidized form, dehydroascorbate, contributing to disruption of redox homeostasis.
An interesting approach to study relevance of antioxidant gene polymorphisms is to compare the level of oxidative stress byproducts in individuals stratified according to their antioxidant gene variants. In the study of oxidative DNA modifications among malignant diseases, 8-OHdG seems to be the most appropriate biomarker. Indeed, bladder cancer patients, carriers of GSTM1-null and GSTA1-variant genotypes, exhibited higher urinary 8-OHdG levels (Savic-Radojevic et al. 2013). In this line, our results on the higher levels of 8-OHdG, a biomarker of free radical-induced DNA-oxidative damage, among ccRCC patients with GSTO2*G/G variant genotype (rs2297235) seem plausible. Taken together, the presence of H2 haplotype, comprising low deglutathionylation and low DHAR activity, might underlie the impaired redox regulation and increased RCC risk. Regarding joint effect of GST polymorphisms and smoking on ccRCC risk, it seems that smoking, as an important source of ROS (Valavanidis et al. 2009), contributes to genotype-associated ccRCC risk in carriers of GSTO2-variant genotype. Indeed, we found that smokers carriers of GSTO2*G/G (rs156697) were at almost 2.5-fold increased RCC risk.
Potential interplay between disrupted redox homeostasis, hypoxia/pseudohypoxia, diminished carcinogenic detoxification and the hypoxia inducible factor (HIF) signaling pathway emerged as important aspect of RCC pathogenesis. Interestingly, oxygen-dependent protein hydroxylases, which represent main regulators of HIF activity by marking HIFα for ubiquitinylation and consequent proteasomal degradation, use Fe(II), 2-oxoglutarate and vitamin C as cofactors. In that context, it might be speculated that vitamin C-dependent inhibition of the HIF pathway may provide additional approach for controlling tumor progression and inflammation (Li and Schellhorn 2007). Inefficient regeneration of vitamin C, as a potential result of GSTO2 polymorphism, might affect HIFα hydroxylation, resulting in HIFα stabilization and downstream overexpression of genes involved in angiogenesis, proliferation, cell migration and invasion, metabolic shift towards glycolysis, survival, erythropoiesis, which eventually contribute to ccRCC tumorigenesis (Mehdi and Riazalhosseini 2017). Recently, HIF-α inhibitors are proposed as potential therapeutic strategy in ccRCC after significant results in inhibiting tumor growth and angiogenesis (Wallace et al. 2016).
Assuming the specific roles of GSTO enzymes in these processes, our results support the hypothesis that GSTO polymorphisms might be associated to the risk of ccRCC, with special emphasis on GSTO2-variant genotype. In the future, investigations of these polymorphisms in conjunction with exposure to specific environmental factors, known to be associated to RCC development, conducted on larger cohort of patients would be warranted.
This work was supported by the Grant 175052 from the Serbian Ministry of Education, Science and Technological Development.
The authors declare no conflict of interest.