Disease Markers
Volume 2017 (2017), Article ID 1474560, 11 pages
https://doi.org/10.1155/2017/1474560

Research Article

Combined Genetic Biomarkers Confer Susceptibility to Risk of Urothelial Bladder Carcinoma in a Saudi Population

Nasser Attia Elhawary,1,2 Anmar Nassir,3,4 Hesham Saada,⁴ Anas Dannoun,¹ Omar Qoqandi,¹ Ammar Alsharif,⁵ and Mohammed Taher Tayeb¹

¹Department of Medical Genetics, Faculty of Medicine, Umm Al-Qura University, P.O. Box 57543, Mecca 21955, Saudi Arabia
²Department of Molecular Genetics, Medical Genetics Center, Faculty of Medicine, Ain Shams University, Cairo 11566, Egypt
³Department of Surgery, Faculty of Medicine, Umm Al-Qura University, Mecca 21955, Saudi Arabia
⁴Department of Urology, King Abdullah Medical City Specialist Hospital, Mecca 21955, Saudi Arabia
⁵Division of Internal Medicine, Al-Noor Specialist Hospital, Mecca 21955, Saudi Arabia

Correspondence should be addressed to Nasser Attia Elhawary; naelhawary@uqu.edu.sa

Received 5 December 2016; Revised 21 January 2017; Accepted 5 February 2017; Published 27 February 2017

Academic Editor: Eric A. Singer

Copyright © 2017 Nasser Attia Elhawary et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We evaluated the associations between seven single nucleotide polymorphisms and susceptibility to urothelial bladder carcinoma (UBC) in a Saudi population. Genomic DNA was taken from buccal cells of 52 patients with UBC and 104 controls for genotyping of GSTT1, GSTM1, rs4646903, rs1048943, TP53 rs1042522, rs1801133, and rs1801394 using PCR and TaqMan® assays. The rs1801133 and rs1801394 variants showed strong associations with UBC (OR = 2.3, ; OR = 2.6, , resp.). Homozygosity of Pro72 conferred a significant double risk in cases compared with controls (30.8% versus 15.4%), but the homozygote Arg/Arg had no effect on risk. Genotypic combinations of GSTM1/GSTT1, rs4646903/rs1048943, and rs1801133/rs1801394 exhibited significant linkage with the disease (, ; , ; and , , resp.). The GSTM1 and rs1042522Arg and rs1801394G variant alleles were more frequent in current smokers with UBC (52.4%, 52.5%, and 64.3%, resp.) than were the corresponding wild-types. Despite some variants having only a slight effect on UBC risk, the interaction effect of combined genetic biomarkers—or even the presence of one copy of a variant allele—is potentially much greater. Perhaps more studies regarding next-generation genetic sequencing and its utility can add to the risk of UBC.

1. Background

The International Human Genome Sequencing Project [1] and the International HapMap Project [2] have generated considerable data on genetic variants and candidate genes for many diseases, including cancer. Urinary bladder cancer (BC) is the second most frequently diagnosed genitourinary cancer worldwide. In Saudi Arabia, the five-year prevalence is 11.6 per 100,000 men and 2.8 per 100,000 women. There have been 14.1 million new cancer cases, 8.2 million cancer deaths, and 32.6 million individuals living with cancer worldwide [3]. Almost 83% of all diagnosed BCs are cases of urothelial bladder carcinoma (UBC) (http://www.scr.org); one-third of UBC is classified as invasive UBC with a very high risk of distant metastases [4]. Aggressive metastatic UBC is untreatable and has a five-year survival rate around 6%. After endoscopic resection, 10–15% of noninvasive UBC will be upstaged into invasive UBC progressing to recurrence, and the patients acquire additional genetic mutations [5].

A growing number of genes regulating enzymes involved in xenobiotic and folate metabolism have been investigated, giving rise to increased knowledge of allelic variants. The impact of polymorphic variants of glutathione S-transferases (GSTs), cytochrome P450s (CYPs), methylenetetrahydrofolate reductase (MTHFR), and methionine synthase reductase (MTRR) on BC risk might be reflected by an association with the frequency of somatic mutations in the TP53 gene [6].

Several polymorphisms of CYP1A1 have been described and current information can be found on the Human CYP allele nomenclature website (http://www.cypalleles.ki.se). Two functional nonsynonymous polymorphisms have been reported to be associated with risk of BC [7, 8]. The CYP1A1 gene regulates the P450-1A1 enzyme that converts polycyclic aromatic hydrocarbons, such as the procarcinogen benzo[a]pyrene, into benzo[a]pyrene diol epoxide carcinogens. These carcinogens, which might be associated with BC in aluminum workers exposed to polycyclic aromatic hydrocarbons and in coal gasification workers exposed to aromatic amines, can be converted in vivo into more hydrophilic active derivatives by the CYP450 enzyme superfamily [9]. The rs1048943 allele results in an alteration of the functional protein and thus increases enzyme activity [10].

Both GSTM1 (MIM 138350) and GSTT1 (MIM 600436) exhibit genetic polymorphisms in populations displaying a large percentage of homozygous gene deletions “null-variant.” As GSTs are involved in several chemical carcinogens, such as benzo[a]pyrene diol epoxide [11, 12], and cellular detoxification, it has been reported that all GSTs enzymes might evolve to protect cells against reactive oxygen metabolites. Inherited loss of both alleles of GSTM1 or GSTT1 is a risk factor for tumors such as bladder, breast, oral, lung, and head and neck cancers (http://www.cancerindex.org/geneweb) [13, 14].

Polymorphisms in genes encoding the folate metabolizing enzymes MTHFR and MTRR might be linked to the etiology of cancer disorders. Deficient folates can lead to decreased DNA methylation and perhaps induce genomic instability or repression of the protooncogenes MTHFR and MTRR [15–17], but few articles have focused on the association between MTHFR or MTRR genes and risk of BC [16, 18, 19]. MTHFR (MIM 607093; 1p36.22) catalyzes the irreversible conversion of 5,10-methylenetetrahydrofolate into 5-methyltetrahydrofolate, which serves as a substrate for the remethylation of homocysteine to methionine with the subsequent synthesis of S-adenosylmethionine (SAM). Cobalamin (vitamin B12) is a cofactor in this reaction, and methionine synthase reductase (MIM 602568; 5p15.31) is required for the maintenance of MTR in its active state. Methionine is then converted into SAM, and DNA methyltransferases transfer the methyl group from SAM to the DNA. The common polymorphism MTRR 66A>G (rs1801394) results in lower enzyme activity [20] and is associated with homocysteine elevation. This mechanism may promote DNA hypomethylation, which is a common feature in early carcinogenesis [21].

The tumor suppressor protein p53 (TP53, MIM 191170) is a sequence-specific transcription factor that works to maintain the integrity of the genome. On its induction in response to DNA damage, TP53 promotes cell cycle arrest in G1/S regulation point and initiates apoptosis if DNA damage proves to be irreparable [22]. An allelic nonsynonymous polymorphism of the TP53 gene (Pro72Arg) can play an important role in the carcinogenesis of many tumors, including urinary BC. This missense change (rs1042522) not only causes TP53 to no longer function as a tumor suppressor but also can affect the development of tumors [23]. It has been found that 72Arg has up to a 15-fold higher apoptotic ability than Pro72 in both inducible cell lines and cells with endogenous TP53 that is homozygous for each of the variants [24].

Historically, the Saudi population is highly heterogeneous and consists of indigenous Asians and Muslim immigrants from the Levant, Africa (who migrated in ancient Islamic times). Consequently, much intermarriage has reinforced genetic drift and gene flow to Saudi people, which might influence the development of various carcinomas, particularly urinary BC. The countries within the Arabian Peninsula region have experienced rapid economic growth over the past few decades, resulting in major changes in the population’s lifestyle. As in other Arab Gulf countries, the overall rate of consanguinity in Saudi community shows that 57.7% of the families are consanguineous; the most frequent marriages are first-degree cousin marriages (28.4%) followed by distant relative marriages (15.2%) and second-cousin marriages (14.6%) [25]. These higher rates result in increasing heterozygosity for several autosomal recessive genetic diseases [26] and even enhance the heritability of some complex multifactorial disorders [27, 28]. In an Egyptian study, parents’ consanguinity (50.8%; ) was shown to have a more prominent role than bilharziasis and smoking in the development of BC [29].

Despite a significant amount of available genetic information on susceptibility to BC, most reports are from Caucasian and Asian populations [30–34]. Data on genetic linkages and their relevance to the development of BC are still limited in Saudi patients [35].

We hypothesize that genetic variants within the GSTM1/GSTT1, CYP450, TP53, MTHFR, and MTRR genes might be each linked alone or linked in combination with the occurrence of UBC. The present study was carried out to examine the associations of GSTM1, GSTT1, rs4646903 and rs1048943, rs1042522, rs1801133, and rs1801394 genetic loci with the risk of UBC tumors in a Saudi population.

2. Subjects and Methods

2.1. Ethics Statement

This study was approved by the Institutional Biomedical Ethics Committee at Faculty of Medicine, Umm Al-Qura University (reference #HAPO-02-K-012), licensed from the National Committee of Medical & Bioethics, KACST, Riyadh (http://bioethics.kacst.edu.sa/About.aspx?lang=en-US). Details of the study were told to participants, and written informed consent was obtained before enrollment.

2.2. Study Population

The study included 52 patients (aged 49–90 years) diagnosed with UBC who were admitted to the Urology Department at King Abdullah City Hospital in Mecca between June 2014 and January 2016. Epidemiologic and clinical characteristics of gender, age at diagnosis, family history of cancer, cigarette smoking habits, alcohol consumption, pathologic tumor stage, tumor grade, and metastasis were recorded for each patient for statistical analyses. The management of the disease—through intravesical bacillus Calmette-Guérin (BCG), immunotherapy, a conservative TURBT-BCG therapy (injection), chemotherapy, chemoradiotherapy, or radical cystectomy—was also recorded. After having done a TURBT, only UBC with high-grade UBC cases was treated with BCG. Patients who previously had cancer, radiotherapy and chemotherapy, or metastasized cancer from non-Saudi or unknown origins were excluded. Patients with immune disease or any histopathologic diagnosis other than UBC were also excluded. Healthy individuals () having no evidence of any clinical phenotypes of malignancies or immune disease were selected as controls (44–89 years) in a routine follow-up at governmental hospitals in Mecca.

2.3. DNA Isolation

DNA samples were extracted from buccal mucosa using the Oragene.DNA-OGR-575 Kit (DNA Genotek Inc., Ottawa, ON, Canada) with some modifications. The full buccal cells were collected within 30 minutes, and the Oragene tube was capped immediately. The cells were incubated with the OGR-lysis buffer in a water bath at 53°C to release the DNA, which was then precipitated by ethanol and dissolved in elution buffer [36].

2.4. Duplex PCR of GSTM1 and GSTT1 Loci

Previously reported primers for the GSTM1 and GSTT1 loci were simultaneously amplified in a single assay and detect the 215-bp and 480-bp fragments, respectively [37]. The PCR amplicon was then analyzed by electrophoresis on a 2% ethidium bromide-stained NuSieve agarose gel (BMA Bioproducts, Rockland, ME). Each sample was run in duplicate, and a positive control was used for each variant. The genotypes of all DNA samples were reassessed twice to confirm the results and ensure reproducibility.

2.5. TaqMan Genotyping Analysis

We adopted TaqMan real-time PCR assays to genotype the selected five SNPs—rs4646903 (C_8879532_20), rs1048943 (C_25624888_50), rs1042522 (C_2403545_10), rs1801133 (C_1202883_20), and rs1801394 (C_3068176_10)—in cases and controls using a 7500 Fast Dx Real-Time PCR System (Applied Biosystems, Life Technologies Inc., USA). Eight negative controls and 104 DNA samples (cases and controls) were included in a 96-well plate to ensure the accuracy of genotyping results. We also repeatedly genotyped 10% of the samples, and the results were 100% concordant.

2.6. Databases of Examined SNPs

The functional consequences of the examined SNP markers were predicted with the web servers; Sorting Intolerant From Tolerant (SIFT; http://sift.jcvi.org/), and Polymorphism Phenotyping (PolyPhen2; http://genetics.bwh.harvard.edu/pph2/).

2.7. Statistical Analysis

Deviations from Hardy-Weinberg equilibrium (HWE) were examined in the study group for healthy controls via chi-squared testing using online software (http://www.oege.org/software/hwe-mr-calc.shtml).   Student’s -test and chi-square test were used to compare demographic and clinical characteristics including age, gender, cigarette smoking, alcohol consumption, tumor grade, pathologic stage, and clinical management. Associations of polymorphisms with UBC were evaluated using ORs (odds ratios), 95% CIs (confidence intervals), and a general z-test with MedCalc Statistical software version 16.4.3 (MedCalc Software bvba, Ostend, Belgium; https://www.medcalc.org; 2016). A two-sided value less than 0.05 was considered to indicate statistical significance for all analyses. We used the Power Software (Franz Faul, Universität Kiel, Germany, ver. 3.1.9.2) (http://www.psycho.uni-duesseldorf.de/abteilungen/aap/gpower3/download-and-register/) to perform a priori power analysis to calculate sufficient sample sizes to achieve adequate power for z-testing of two independent proportions. A priori sample size estimations were performed using known information on the common allele frequencies in UBC cases and healthy controls, a criterion probability of , and a power sensitivity of 80%.

3. Results

3.1. Characteristics of the Study Population

For the study, 52 eligible Saudi individuals with UBC (48 men : 4 women; 12 : 1) and 104 healthy controls (92 men : 12 women; 7.7 : 1) were enrolled. The mean age of patients was years, with no significant difference when compared with controls (). There was no significant difference between cases and controls in terms of the percentage of current cigarette smokers (80.8% versus 76.9%, resp.; ). Significant differences () were found in the proportion of cases with a high tumor grade (35 cases, 67%) and the proportion with a low tumor grade (17 cases, 32.7%). The percentages of UBC cases with specific tumor stages were 59.6% for Ta, 17.3% for T1, and 23.1% for T2 () (Table 1). The management course of UBC cases was significantly predominant with BCG immunotherapy (76.9%; , ). Conservative management with TURBT-BCG injection was significantly higher than other types of treatment (80.8% versus 9.6% for radical cystectomy, 5.8% for chemoradiotherapy, and 3.8% for chemotherapy).

Table 1: Epidemiologic and clinical characteristics in UBC cases.

Parameter
UBC cases  

(95% CI), value

Age at examination (range, years)
46–90

Mean age ± SD (years)
61.8 ± 10.63
38.5

Gender, males
48 (92.3)
28.9

Cigarette smoking
42 (80.8)
25.1

Alcohol consumption (yes)
—
—

Tumor grade:

    Low-
17 (32.7)
9.2

    High-
35 (67.3)
20.6

Pathologic stage:

31 (59.6)
18.1

9 (17.3)
4.1

12 (23.1)
6.0

Metastasis (yes)
0 (0.0)

Management:

40 (76.9)
23.8

    Conservative
42 (80.8)
25.1

    Radical
5 (9.6)
1.5 (3.2–21.0)

    Chemoradiotherap
3 (5.8)
0.3 (1.2–16.0)

2 (3.8)
0.4 (0.5–13.1)

Student’s t-test. Values are mean ± standard deviation (SD).

of patients, with percentages in parentheses.

highly significant difference ().

, bacillus Calmette-Guérin, a weakened bacterium intravesically introduced via a catheter.

therapy,” TURBT plus BCG.

significant difference ().

3.2. Hardy-Weinberg Equilibria of SNPs

All controls were in HWE at the rs4646903 T>C (, ), rs1048943 A>G (, ), rs1042522 Pro>Arg (, ), and MTHFR rs1801133 C>T (, ), and MTRR rs1801394 A>G loci (, ). HWE could not be tested for the GSTM1 and GSTT1 null genotypes because of the inability of the conventional restriction fragment length polymorphism protocol to separate heterozygous carriers of the deletion polymorphisms.

3.3. Allele Frequencies and Genotype Distributions

The frequencies of the GSTM1 null and GSTT1 null alleles were higher in UBC cases than in controls (0.46 versus 0.38; OR = 1.4, 95% CI, 0.9–2.2, , and for GSTM1; 0.12 versus 0.08; OR = 1.6, 95% CI, 0.7–3.4, , and for GSTT1) (Table 2). Genotypic distribution of the GSTM1 and GSTT1 polymorphisms showed a higher frequency of the null-variants in UBC cases than in controls (46.2% versus 38.5% for the homozygous GSTM1 null genotype and 11.5% versus 7.7% for the homozygous GSTT1 null genotype). There were no differences in the homozygous null genotypes of the GSTM1 and GSTT1 alleles between the UBC cases and the controls (, ; , ).

Table 2: Genotype distributions and allele frequencies of examined SNPs in UBC cases.

SNP ID
Cases
Controls
( value
( value OR (95% CI); ( value

(%)
(%)

GSTM1:

“+/+”
28 (53.8)
64 (61.5)
0.8
0.8

“−/−”
24 (46.2)
40 (38.5)
0.8

“+”
56 (0.54)
128 (0.62)

1.4 (0.9–2.2); 1.3

“−”
48 (0.46)
80 (0.38)

GSTT1:

“+/+”
46 (88.5)
96 (92.3)
0.6
0.6

“−/−”
6 (11.5)
8 (7.7)
0.6

“+”
92 (0.88)
192 (0.92)

1.6 (0.7–3.4); 1.1

“−”
12 (0.12)
16 (0.08)

CYP1A1 T6235C (rs4646903):

TT
34 (65.4)
60 (57.7)
0.8
1.7

TC
16 (30.8)
42 (40.4)
1.4

CC
2 (3.8)
2 (1.9)
0.5

T
84 (0.81)
162 (0.78)

0.8 (0.5–1.5); 0.6

C
20 (0.19)
46 (0.22)

CYP1A1 A4889G (rs1048943):

AA
46 (88.5)
92 (88.5)
0.0
4.5

AG
4 (30.8)
12 (11.5)
8.7

GG
2 (26.9)
0 (0.0)
30.5 (<0.0001

A
96 (0.92)
196 (0.94)

1.4 (0.5–3.4); 0.7

G
8 (0.08)
12 (0.06)

TP53 Pro72Arg (rs1042522):

Pro/Pro
16 (30.8)
16 (15.4)
5.0
5.6

Pro/Arg
22 (42.3)
60 (57.7)
3.3

Arg/Arg
14 (26.9)
28 (26.9)
0.0

Pro
54 (0.52)
92 (0.44)

0.7 (0.5–1.2); 1.3

Arg
50 (0.48)
116 (0.56)

MTHFR C677T (rs1801133):

CC
33 (63.5)
88 (84.6)
8.8
14.4

CT
14 (26.9)
16 (15.4)
2.9

TT
5 (9.6)
0 (0.0)
10.2

C
80 (0.77)
192 (0.88)

3.6 (1.8–7.1); 3.7

T
24 (0.23)
16 (0.12)

MTRR A66G (rs1801394):

AA
0 (0.0)
40 (38.5)
26.8 (<0.0001
35.6

AG
32 (61.5)
54 (51.9)
1.3

GG
20 (38.5)
10 (9.6)
18.5 (<0.0001

A
32 (0.31)
134 (0.64)

4.1 (2.5–6.7); 5.5 (<0.0001

G
72 (0.69)
74 (0.36)

“+”: present allele; “−”: null allele; “−,−”: homozygous null genotype; “+,+”: homozygous present genotype; OR: odds ratio; CI: confidence interval.

difference between two similar genotypes in cases and controls.

difference between genotypes of a specific SNP in cases compared with controls.

of allele frequencies of a specific SNP.

The differences between the variant alleles with the cases and controls were not statistically significant for three SNPs: rs4646903 (OR = 0.8, 95% CI, 0.5–1.5, , and ), rs1048943 (OR = 1.4, 95% CI, 0.5–3.4, , and ), and rs1042522 (OR = 0.7, 95% CI, 0.5–1.2, , and ). But there were significant differences in this issue for two SNPs: rs1801133 (OR = 3.6, 95% CI, 1.8–7.1, , and ) and rs1801394 (OR = 4.1, 95% CI, 2.5–6.7, , and ) (Table 2).

There were strong, statistically significant differences between the genotypes of cases and controls for the MTHFR rs1801133 SNP (, ) and the MTRR rs1801394 SNP (, ). For the rs1801394 SNP, there was an absence of AA wild-type carriers among UBC cases but a 38.5% prevalence of AA wild-type carriers among controls (, ). GG-variants were more frequent among UBC cases than among controls (38.5% versus 9.6%; , ). For the rs1801133 SNP, there were significant differences between cases and controls for the CC wild-type (, ) and TT-variants (, ) genotypes. On the other hand, there were no linkage disequilibria between the rs4646903 and rs1048943 SNPs and UBC ( and , resp.) (Table 2). For the rs1042522 SNP, the homozygote Pro/Pro conferred a double risk with a significant difference in cases compared with controls (30.8% versus 15.4%; ), but the homozygote Arg/Arg had no effect on risk (). None of the genotypes of the rs4646903 SNP or rs1048943 SNP were linked to the risk of UBC ( and , resp.).

3.4. Genotypic Combinations

GSTM1/GSTT1 genotypic combinations were significantly different in cases and controls (, ), with the “−,−/+,+” combination being the most prevalent in UBC cases (53.8% versus 38.5%; , ). There were no cases or controls with the combined null genotype “−,−/−,−". In contrast, the “+,+/+,+” combination was more prevalent in controls than in cases (53.8% versus 34.6%), with a strong statistical significance (, ). These outcomes might suppose the pathological role of the GSTM1 null genotype in the risk of the disease. There was a significant difference in the genotypic combinations of rs4646903 T>C/rs1048943 A>G in cases versus controls (, ), including a significant difference in the homozygous variant (CC/GG), which was absent in controls but seen in 3.8% of cases (, ). There was also a strong significant difference regarding rs1801133 C>T/rs1801394 A>G combinations in cases versus controls (, ), with the absence of CT/GG in controls (, ) (Table 3).

Table 3: Genotypic combinations of examined SNPs showing risk of UBC in cases and controls.

Genotypic combination
Cases  
(%)
Controls  
(%)
χ ² ( value
χ ² ( value

GSTM1/GSTT1:

+,+/+,+
36 (34.6)
112 (53.8)
10.2 (0.001)

−,−/+,+
56 (53.8)
80 (38.5)
6.6 (0.01)
10.3 (0.006)

+,+/−,−
12 (11.5)
16 (7.7)
1.2 (0.3)

−,−/−,−
0 (0.0)
0 (0.0)
NA
NA

rs4646903 T>C/ rs1048943 A>G:

TT/AA
72 (69.2)
120 (57.7)
3.9 (0.04)

TC/AG
8 (7.7)
24 (11.5)
1.1 (0.3)

TC/AA
20 (19.2)
64 (30.8)
4.7 (0.03)
13.9 (0.003)

CC/GG
4 (3.8)
0 (0.0)
4.0 (0.047)

rs1801133 C>T/ rs1801394 A>G:

CC/AA
0 (0.0)
8 (3.8)
4.2 (0.04)

CC/AG
56 (53.8)
120 (57.7)
0.4 (0.5)

CC/GG
20 (19.2)
32 (15.4)
0.7 (0.4)

CT/AA
0 (0.0)
8 (3.8)
4.2 (0.04)
20.4 (0.0004)

CT/AG
20 (19.2)
40 (19.2)
0.0 (1.0)

CT/GG
8 (7.7)
0 (0.0)
8.3 (0.0001)

“+”: present allele; “−”: null allele; “+,+”: homozygous present genotype; “−,−”: homozygote null genotype; OR: odds ratio; CI: confidence interval.

significant difference between a specific combined genotype in cases and controls.

significant difference between combined genotypes in cases compared with controls.

3.5. Effect of Smoking on Allelic Distributions

Table 4 shows the differential allelic frequency of the examined SNPs in current smokers with UBC (, 80.8%). The variant alleles of GSTM1 null, rs1042522 Arg, and MTRR rs1801394G were more prevalent in current smokers with UBC disease (52.4%, 52.5%, and 64.3%, resp.) than were the wild-type alleles. However, the difference was only statistically significant for the rs1801394G allele (, ). Significant differences were shown for other examined SNPs but linked dominantly with the wild-type alleles (Table 4).

Table 4: Effect of SNP allele percentages on phenotypes of 42 smokers with UBC.

SNP ID
Wild-type allele
(%)
Variant allele (%)
( value)

GSTM1 (+,−)
40 (47.6)
44 (52.4)
0.4 (0.5)

GSTT1 (+,−)
72 (85.7)
12 (14.3)
85.1 (<0.0001)

rs4646903 T>C
70 (83.3)
14 (16.7)
74.1 (<0.0001)

rs1048943 A>G
78 (92.9)
6 (7.1)
122.9 (<0.0001)

rs1042522 Pro>Arg
40 (47.5)
44 (52.5)
0.4 (0.5)

rs1801133 C>T
74 (88.1)
10 (11.9)
97.0 (<0.0001)

rs1801394 A>G
30 (35.7)
54 (64.3)
13.7 (0.0002)

3.6. Prediction of Examined SNP Markers

The SIFT and PolyPhen2 were utilized to predict amino acid substitution effects and to predict damaging effects. The SIFT algorithm predicted the four nonsynonymous SNPs; rs1048943 (p.I462V), rs1042522 (p.P72R), rs1801133 (p.A222V), and rs1801394 (p.I22M) were deleterious to effect protein function with scores, 0.16, 0.16, 0.03, and 0.01, respectively. The PolyPhen2 also predicted that all these SNPs were possibly damaging with a scoring range 0.265 to 0.962, except for the CYP1A1 rs1048943 polymorphism which showed benign phenotypes (Table 5).

Table 5: SIFT and PolyPhen2 functional predictions of SNP markers.

SNP ID
Allele
Amino acid position
SIFT
PolyPhen2

Score
Prediction
Confidence
Score
Prediction

rs1048943
T>C
p.I462V
0.16
Tolerated
High
0.18
Benign

rs4646903
A>G

—
—
—
—
—

rs1042522
G>C
p.P72A
0.16
Tolerated
High
0.265
Possibly damaging

rs1801133
C>T
p.A222V
0.03
Damaging
High
0.235
Possibly damaging

rs1801394
A>G
p.I22M
0.01
Damaging
High
0.962
Possibly damaging

SIFT, Sorting Intolerant From Tolerant (http://sift.jcvi.org/); PolyPhen2, Polymorphism phenotyping (http://genetics.bwh.harvard.edu/pph2/).

: not applicable as it is in the noncoding 3′ flanking untranslated region. Low confidence means that the protein alignment does not have enough sequence diversity. Because the position artificially appears to be conserved, an amino acid may incorrectly be predicted to be damaging.

4. Discussion

This hospital-based case-control study presents the first investigation, to our knowledge, of potential associations between the GSTM1 null, GSTT1 null, rs4646903 T>C, rs1048943 A>G, rs1042522 Pro>Arg, rs1801133 C>T, and rs1801394 A>G polymorphisms and UBC susceptibility in a Saudi population.

We found no deviation in our controls from HWE for all examined SNPs. The observed frequency of the MTRR rs1801394 G/G homozygote was higher than expected (38.5% versus 9.6%). The duplex PCR could not differentiate between the “+,+” homozygote or the “+,−" heterozygote carrying the null allele. However, we might consider the individual with a heterozygous genotype “+,−" as a homozygote with “+,+” genotype. This could result in a weaker impact of the null alleles on UBC risk.

Overall, our results provide strong evidence of associations between two SNPs—rs1801133 C>T and rs1801394 A>G—and risk of UBC in this population. The variants rs1801133T and rs1801394G can potentially affect the risk of UBC ( and , resp.). None of the GSTM1 null, GSTT1 null, rs4646903, rs1048943, or TP53 rs1042522 variant was linked to risk of the disease in this study. Our findings regarding combined GSTM1/GSTT1 genotypes support other studies showing that the GSTM1 null genotype is significantly associated with increased UBC risk, irrespective of the GSTT1 genotype status (−,−/+,+; ). The inherited GSTM1/GSTT1 null allele combination (−,−/−,−) was absent in our cases, which is inconsistent with studies in Turkish populations [38]. Sometimes, it is difficult to expect how the GSTM1 and GSTT1 enzymes may potentially influence BC susceptibility, as they have multifunctional roles in metabolic pathways.

Eight case-control studies have produced conflicting results regarding the association between CYP1A1 polymorphisms and BC susceptibility [7, 38–43]. However, most of these studies have observed no significant association with BC risk [44]. Among them, a few evaluated the association between the rs4646903 locus and BC risk and found slight trends toward association in German, Chinese, and North Indian populations [39, 40, 43, 44]. Despite the weak impact of the separate CYP SNPs (rs4646903 and rs1048943) in the present study, the combined genotype has a significant linkage with UBC ().

The present study found that the Pro/Arg heterozygotes of TP53 rs1042522 were prevalent among both cases and controls, with no significant difference between the two (42.3% versus 57.7%, resp.; ). However, some studies support a passive risk of BC among Arg/Arg carriers in the Saudi population [45]. Also, in our study, homozygosity of Pro72 had a double effect on risk for UBC cases when compared with controls (), but homozygosity 72Arg had no effect on our Saudi UBC risk (). It has been found that Pro72 homozygotes can improve survival after development of cancer compared with 72Arg homozygotes [46]. It is known that 72Arg increases the ability of TP53 to translocate to mitochondria and induce cell death, as Pro72 exhibits lower apoptotic potential but increases cellular arrest in G1 of the cell cycle [46]. Previous studies have demonstrated that the functional polymorphisms of the heterozygous Pro/Arg are associated with susceptibility to a variety of cancers, such as thyroid cancer, lung cancer, and cervical cancer [47–49]. Other studies have reported that Pro/Pro and Arg/Arg homozygosity may increase susceptibility to thyroid cancer and cervical cancer [45, 49, 50].

The variant homozygotes of the rs1801133 and rs1801394 SNPs can clearly influence the risk of UBC ( and , resp.). For rs1801133TT, our data outcomes are consistent with those in German (Caucasian), Chinese, Netherlander, and Turkish populations [15, 51–56]. In contrast, a meta-analysis found no association between the 677C>T polymorphism and BC in the Chinese population [57]. Besides, several studies reported no linkage disequilibrium in 677C>T with BC risk in populations from the United States, Spain, Tunis, Iran, and Taiwan [16, 18, 58–62]. Few studies have dealt with the rs1801394 SNP and those that have inconsistent results regarding association with BC [16, 18, 19].

Despite a limited trend of associations between some SNPs examined in this Saudi sample group and risk of UBC, this work clearly showed that genotypic combinations of GST (GSTM1/GSTT1) and CYP (rs4646903/rs1048943) could reinforce statistical linkage disequilibrium with the disease ( and , resp.). Moreover, genes involved in folate metabolism (rs1801133/rs1801394) clearly exhibit strong linkage with UBC (), and this outcome is consistent with previous studies regarding the same combined genotypes of MTHFR/MTRR loci in a Tunisian population [16].

The incidence of BC in men is higher in the northwestern part of West Asia (Turkey and Lebanon) than in other parts, and cigarette smoking is proposed as the most important risk factor for BC in these countries [63, 64]. Among the current smokers with UBC in our study, the variant alleles of GSTM1 null, TP53 rs1042522Arg, and MTRR rs1801394G were more frequent than the wild-type alleles. Slightly more smokers with UBC had GSTM1 null alleles than wild-type alleles. Bell et al. [65] have previously stated that the association among the GSTM1 null-variant, cigarette smoking, and BC has been controversial, indicating higher risks in smokers due to lack of GSTM1. Meta-analyses and pooled studies have found no or only weak evidence for an association between the GSTM1 null-variant and smoking habits [66–69]. Rothman et al. [69] have reported an even higher OR for nonsmokers than for ever-smokers (1.71 versus 1.47). However, only the variant G allele of the rs1801394 SNP was strongly associated with the disease (). Clearly, there was no effect of smoking habits associated with the GSTT1 null, rs4646903 T>C, rs1048943 A>G, and MTHFR 677C>T SNPs.

Occupational exposure (e.g., to diesel) has also been suggested as a risk factor for BC [64]. However, any effect of occupational history in our study was clearly passive, as these UBC patients were working in business affairs (92.3%) or as house wives (7.7%). In one study, BC patients exposed to polycyclic aromatic hydrocarbons (i.e., people with occupational histories in coal, iron, and steel industries) presented with high percentages of GSTM1 null genotypes [70, 71]. Fortunately, another study suggested that polycyclic aromatic hydrocarbons had no effect on GSTM1 once these industries had been closed, and the GSTM1 genotypes were equal in both cases and controls (GSTM1 “−” = 52%) [72].

Based on the functional prediction of SIFT and PolyPhen2 web links, three examined SNPs (TP53 rs1042522, MTHFR rs1801133, and MTRR rs1801394) are suggested to have serious effects of damaging on the functions and risks to UBC polymorphic phenotypes. But the CYP rs1048943 SNP shares a benign prediction in this issue.

Conflicting results in pinning down a genetic association for urinary BC are not uncommon, as poor replication of studies can arise from several different factors. Following are a few example of how this can happen. First, when compared with our study, some studies have had populations with admixed ethnicities, different sources of BC/control subjects, or small sample sizes, which would lessen the strength of the overall results to be consistent with our Saudi population. In addition, the literature includes other pathological types of cancer (e.g., squamous or adenocarcinoma) in addition to UBC, which would undoubtedly influence the outcomes of association of the genetic markers to BC. However, the present study did focus only on the UBC rather than other types which consequently strengthen the reliability of our outcomes of associations with the disease. Second, various technologies have been used in previous studies to examine and discriminate different genotypes within a specific SNP, resulting in a broad range of false-negative or false-positive genotypic distributions. Our study primarily used a conventional PCR protocol, as in other publications, to identify the GSTM1 and GSTT1 null-variants; hence, the slight trends of association of the null alleles with UBC risk were weaker than expected in our outcomes. Third, our post hoc statistical analyses for the MTHFR rs1801133 and MTRR rs1801394 SNPs revealed powers of 39.8% and 98%, resp., among our 156 participants. Recruiting 372 and 70 participants for these SNPs, respectively, would be sufficient to reach a power of 80%. Usually smaller sample sizes lead to lower post hoc power of detection. According to our “a priori” calculations, we would be in need of sampling sizes of 596, 882, 2842, 2554, and 611 in both cases and controls for GSTM1, GSTT1, rs4646903, rs1048943, and rs1042522, respectively, to get 80% power of detection. Overcoming limitations such as these can be difficult, but meta-analysis may offer a feasible option, as it permits surveying a wide set of subjects, thereby enhancing post hoc power and allowing for a more broadly based analysis of previously available data. Recruiting more participants within a reasonable time frame from a single hospital or clinical center would have been challenging in our study, but replication of our results through larger, multicenter genetic association studies will be of interest.

Recently, Li and Chen [43] utilized the available omics data of next-generation sequencing to evaluate the miRNA regulatory network biomarkers (e.g., miR-200a, miR-200b) to investigate BC mechanisms and design multiple drug combinations for treating BC. Moreover, the role of miRNAs in several BC scenarios has been described by Matullo et al. [73], suggesting they could become relevant clinical markers. Ward et al. [34] have utilized next-generation sequencing to amplify DNA from the urine of BC patients and have detected mutations with 70% sensitivity; 54% of the mutations were in TERT, 28% in FGFR3, 13% in PIK3CA, 13% in TP53, 5% in RXRA, and 1% in HRAS mutations. Some additional genomic alterations have included TP53 (54%) and RB1 (17%) [74].

5. Conclusions

In conclusion, this study presents evidence that the MTHFR rs1801133 and MTRR rs1801394 variants influence genes involved in folate metabolism, remarkably increasing the risk of UBC in the Saudi population. Moreover, the GSTM1, GSTT1, CYP (rs4646903 and rs1048943), and TP53 rs1042522 Pro>Arg variants showed a slight trend of linkage disequilibria but can interact additively to increase the risk of UBC. So far, these outcomes should be taken with caution, as these SNPs do not act alone to explain such complex multifactorial malignancies. More linkage outcomes based on a whole-genome approach instead of a single gene or a few candidate genes may be useful for discovering new genes associated with susceptibility to BC. The association between the GSTM1 variant and current smokers with UBC has been controversial, indicating higher risks in smokers due to the GSTM1 null-variant; however, no potential role was found for cytochrome P450, MTHFR 677C>T, or MTRR 66A>G in smoking habits in this study. These associations or linkages are predictive and intriguing, but they need to be confirmed by additional studies. Ongoing analyses of UBC patient DNA using next-generation sequencing protocols are being performed.

Competing Interests

The authors report no conflict of interests regarding the publication of this paper.

Authors’ Contributions

Nasser Attia Elhawary, Anmar Nassir, Hesham Saada, and Mohammed Taher Tayeb contributed equally to this work. Anas Dannoun, Omar Qoqandi, and Ammar Alsharif also contributed equally to this work.

Acknowledgments

The authors thank medical residents Drs. Abrar Filfilan, Amani Koshak, Areej Qadah, Ayat Aldosari, and Ayat Azhari (Department of Urology, King Abdullah Medical City, Mecca, Saudi Arabia), for collecting patients’ clinical data and sampling.

References

International Human Genome Sequencing Consortium, “Finishing the euchromatic sequence of the human genome,” Nature, vol. 431, no. 7011, pp. 931–945, 2004.

The International HapMap Consortium, “International HapMap Project,” Nature, vol. 426, no. 6968, pp. 789–796, 2003.

G. Roshandel, M. Boreiri, A. Sadjadi, and R. Malekzadeh, “A diversity of cancer incidence and mortality in west Asian populations,” Annals of Global Health, vol. 80, no. 5, pp. 346–357, 2014.

D. S. Kaufman, W. U. Shipley, and A. S. Feldman, “Bladder cancer,” The Lancet, vol. 374, no. 9685, pp. 239–249, 2009.

R. L. Siegel, K. D. Miller, and A. Jemal, “Cancer statistics, 2015,” CA Cancer Journal for Clinicians, vol. 65, no. 1, pp. 5–29, 2015.

R. Thier, T. Brüning, P. H. Roos et al., “Markers of genetic susceptibility in human environmental hygiene and toxicology: the role of selected CYP, NAT and GST genes,” International Journal of Hygiene and Environmental Health, vol. 206, no. 3, pp. 149–171, 2003.

H. Bartsch, U. Nair, A. Risch, M. Rojas, H. Wikman, and K. Alexandrov, “Genetic polymorphism of CYP genes, alone or in combination, as a risk modifier of tobacco-related cancers,” Cancer Epidemiology Biomarkers and Prevention, vol. 9, no. 1, pp. 3–28, 2000.

V. Harth, T. Brüning, J. Abel et al., “Real-time genotyping of cytochrome P4501A1 A4889G and T6235C polymorphisms,” Molecular and Cellular Probes, vol. 15, no. 2, pp. 93–97, 2001.

J. P. S. Grando, H. Kuasne, R. Losi-Guembarovski et al., “Association between polymorphisms in the biometabolism genes CYP1A1, GSTM1, GSTT1 and GSTP1 in bladder cancer,” Clinical and Experimental Medicine, vol. 9, no. 1, pp. 21–28, 2009.

K. Kawajiri, H. Eguchi, K. Nakachi, T. Sekiya, and M. Yamamoto, “Association of CYP1A1 germ line polymorphisms with mutations of the p53 gene in lung cancer,” Cancer Research, vol. 56, no. 1, pp. 72–76, 1996.

B. Ketterer, J. M. Harris, G. Talaska et al., “The human glutathione S-transferase supergene family, its polymorphism, and its effects on susceptibility to lung cancer,” Environmental Health Perspectives, vol. 98, pp. 87–94, 1992.

Y. Gao, Y. Cao, A. Tan, C. Liao, Z. Mo, and F. Gao, “Glutathione S-transferase M1 polymorphism and sporadic colorectal cancer risk: an updating meta-analysis and HuGE review of 36 case-control studies,” Annals of Epidemiology, vol. 20, no. 2, pp. 108–121, 2010.

A. Sharma, A. Pandey, S. Sardana, S. Ashok, and J. K. Sharma, “Genetic polymorphisms of GSTM1 and GSTT1 genes in Delhi and comparison with other Indian and global populations,” Asian Pacific Journal of Cancer Prevention, vol. 13, no. 11, pp. 5647–5652, 2012.

L.-P. Zhou, H. Luan, X.-H. Dong, G.-J. Jin, D.-L. Man, and H. Shang, “Genetic variants of CYP2D6 gene and cancer risk: a HuGE systematic review and meta-analysis,” Asian Pacific Journal of Cancer Prevention, vol. 13, no. 7, pp. 3165–3172, 2012.

B. T. Heijmans, J. M. A. Boer, H. E. D. Suchiman et al., “A common variant of the methylenetetrahydrofolate reductase gene (1p36) is associated with an increased risk of cancer,” Cancer Research, vol. 63, no. 6, pp. 1249–1253, 2003.

K. Rouissi, S. Ouerhani, E. Oliveira et al., “Polymorphisms in one-carbon metabolism pathway genes and risk for bladder cancer in a Tunisian population,” Cancer Genetics and Cytogenetics, vol. 195, no. 1, pp. 43–53, 2009.

F. Coppedè, F. Migheli, A. Lopomo et al., “Gene promoter methylation in colorectal cancer and healthy adjacent mucosa specimens: correlation with physiological and pathological characteristics, and with biomarkers of one-carbon metabolism,” Epigenetics, vol. 9, no. 4, pp. 621–633, 2014.

L. E. Moore, N. Malats, N. Rothman et al., “Polymorphisms in one-carbon metabolism and trans-sulfuration pathway genes and susceptibility to bladder cancer,” International Journal of Cancer, vol. 120, no. 11, pp. 2452–2458, 2007.

K. Rouissi, N. Stambouli, R. Marrakchi et al., “Smoking and polymorphisms in folate metabolizing genes and their effects on the histological stage and grade for bladder tumors,” Bulletin du Cancer, vol. 98, no. 2, pp. E1–E10, 2011.

D. L. Harmon, D. C. Shields, J. V. Woodside et al., “Methionine synthase D919G polymorphism is a significant but modest determinant of circulating homocysteine concentrations,” Genetic Epidemiology, vol. 17, no. 4, pp. 298–309, 1999.

D. Leclerc, E. Campeau, P. Goyette et al., “Human methionine synthase: cDNA cloning and identification of mutations in patients of the cblG complementation group of folate/cobalamin disorders,” Human Molecular Genetics, vol. 5, no. 12, pp. 1867–1874, 1996.

M. Hollstein, D. Sidransky, B. Vogelstein, and C. C. Harris, “p53 mutations in human cancers,” Science, vol. 253, no. 5015, pp. 49–53, 1991.

D. J. McConkey, S. Lee, W. Choi et al., “Molecular genetics of bladder cancer: emerging mechanisms of tumor initiation and progression,” Urologic Oncology: Seminars and Original Investigations, vol. 28, no. 4, pp. 429–440, 2010.

P. Dumont, J. I.-J. Leu, A. C. Della Pietra III, D. L. George, and M. Murphy, “The codon 72 polymorphic variants of p53 have markedly different apoptotic potential,” Nature Genetics, vol. 33, no. 3, pp. 357–365, 2003.

M. A. F. El-Hazmi, A. R. Al-Swailem, A. S. Warsy, A. M. Al-Swailem, R. Sulaimani, and A. A. Al-Meshari, “Consanguinity among the Saudi Arabian population,” Journal of Medical Genetics, vol. 32, no. 8, pp. 623–626, 1995.

E. H. Jiffri, N. Bogari, K. H. Zidan, S. Teama, and N. A. Elhawary, “Molecular updating of β-thalassemia mutations in the upper Egyptian population,” Hemoglobin, vol. 34, no. 6, pp. 538–547, 2010.

N. A. Elhawary, M. T. Tayeb, S. Abdel-Ghafar, M. Rashad, and A.-A. Alkhotani, “TNF-238 polymorphism may predict bronchopulmonary dysplasia among preterm infants in the Egyptian population,” Pediatric Pulmonology, vol. 48, no. 7, pp. 699–706, 2013.

N. A. Elhawary, N. Bogari, E. H. Jiffri, M. Rashad, A. Fatani, and M. Tayeb, “Transporter TAP1-637G and immunoproteasome PSMB9-60H variants influence the risk of developing vitiligo in the saudi population,” Disease Markers, vol. 2014, Article ID 260732, 8 pages, 2014.

A. H. Zarzour, M. Selim, A. A. Abd-Elsayed, D. A. Hameed, and M. A. AbdelAziz, “Muscle invasive bladder cancer in Upper Egypt: the shift in risk factors and tumor characteristics,” BMC Cancer, vol. 8, article 250, 2008.

W.-X. Sun, Y.-H. Chen, Z.-Z. Liu et al., “Association between the CYP1A2 polymorphisms and risk of cancer: a meta-analysis,” Molecular genetics and genomics, vol. 290, no. 2, pp. 709–725, 2015.

Y. Zhang, D. Yang, J.-H. Zhu, M.-B. Chen, W.-X. Shen, and J. He, “The association between NQO1 Pro187Ser polymorphism and urinary system cancer susceptibility: a meta-analysis of 22 studies,” Cancer Investigation, vol. 33, no. 2, pp. 39–46, 2015.

V. Chandra, J. J. Kim, U. Gupta, B. Mittal, and R. Rai, “Impact of DCC (rs714) and PSCA (rs2294008 and rs2976392) gene polymorphism in modulating cancer risk in asian population,” Genes, vol. 7, no. 2, article 5, 2016.

J. D. Figueroa, C. D. Middlebrooks, A. R. Banday et al., “Identification of a novel susceptibility locus at 13q34 and refinement of the 20p12.2 region as a multi-signal locus associated with bladder cancer risk in individuals of European ancestry,” Human Molecular Genetics, vol. 25, no. 6, pp. 1203–1214, 2016.

D. G. Ward, L. Baxter, N. S. Gordon et al., “Multiplex PCR and next generation sequencing for the non-invasive detection of bladder cancer,” PLoS ONE, vol. 11, no. 2, Article ID e0149756, 2016.

T. O. Al Hussain and M. Akhtar, “Molecular basis of urinary bladder cancer,” Advances in Anatomic Pathology, vol. 20, no. 1, pp. 53–60, 2013.

A. H. Arab and N. A. Elhawary, “Association between ANKK1 (rs1800497) and LTA (rs909253) genetic variants and risk of Schizophrenia,” BioMed Research International, vol. 2015, Article ID 821827, 8 pages, 2015.

R. Hezova, J. Bienertova-Vasku, M. Sachlova et al., “Common polymorphisms in GSTM1, GSTT1, GSTP1, GSTA1 and susceptibility to colorectal cancer in the Central European population,” European Journal of Medical Research, vol. 17, no. 1, article 17, 2012.

U. Berber, I. Yilmaz, O. Yilmaz et al., “Cyp1a1 (ile462val), CYP1B1 (ala119ser and val432leu), GSTM1 (null), and GSTT1 (null) polymorphisms and bladder cancer risk in a turkish population,” Asian Pacific Journal of Cancer Prevention, vol. 14, no. 6, pp. 3925–3929, 2013.

J. Brockmöller, I. Cascorbi, R. Kerb, and I. Roots, “Combined analysis of inherited polymorphisms in arylamine N-acetyltransferase 2, glutathione S-transferases M1 and T1, microsomal epoxide hydrolase, and cytochrome P450 enzymes as modulators of bladder cancer risk,” Cancer Research, vol. 56, no. 17, pp. 3915–3925, 1996.

D. S. Srivastava, A. Mandhani, and R. D. Mittal, “Genetic polymorphisms of cytochrome P450 CYP1A1 (∗2A) and microsomal epoxide hydrolase gene, interactions with tobacco-users, and susceptibility to bladder cancer: a study from North India,” Archives of Toxicology, vol. 82, no. 9, pp. 633–639, 2008.

L. Fontana, L. Delort, L. Joumard et al., “Genetic polymorphisms in CYP1A1, CYP1B1, COMT, GSTP1 and NAT2 genes and association with bladder cancer risk in a French cohort,” Anticancer Research, vol. 29, no. 5, pp. 1631–1635, 2009.

T. Öztürk, Ö. T. Kahraman, B. Toptaş et al., “The effect of CYP1A1 and GSTM1 gene polymorphisms in bladder cancer development in a Turkish population,” In Vivo, vol. 25, no. 4, pp. 663–668, 2011.

C.-W. Li and B.-S. Chen, “Network biomarkers of bladder cancer based on a genome-wide genetic and epigenetic network derived from next-generation sequencing data,” Disease Markers, vol. 2016, Article ID 4149608, 18 pages, 2016.

Y. Lu, X.-L. Zhang, L. Xie et al., “Lack of association between CYP1A1 polymorphisms and risk of bladder cancer: a meta-analysis,” Asian Pacific Journal of Cancer Prevention, vol. 15, no. 9, pp. 4071–4077, 2014.

I. Mabrouk, S. Baccouche, R. El-Abed et al., “No evidence of correlation between p53 codon 72 polymorphism and risk of bladder or breast carcinoma in Tunisian patients,” Annals of the New York Academy of Sciences, vol. 1010, pp. 764–770, 2003.

D. D. Ørsted, S. E. Bojesen, A. Tybjærg-Hansen, and B. G. Nordestgaard, “Tumor suppressor p53 Arg72Pro polymorphism and longevity, cancer survival, and risk of cancer in the general population,” Journal of Experimental Medicine, vol. 204, no. 6, pp. 1295–1301, 2007.

Y.-C. Wang, C.-Y. Chen, S.-K. Chen, Y.-Y. Chang, and P. Lin, “p53 Codon 72 polymorphism in Taiwanese lung cancer patients: association with lung cancer susceptibility and prognosis,” Clinical Cancer Research, vol. 5, no. 1, pp. 129–134, 1999.

I. Zehbe, G. Voglino, E. Wilander et al., “p53 codon 72 polymorphism and various human papillomavirus 16 E6 genotypes are risk factors for cervical cancer development,” Cancer Research, vol. 61, no. 2, pp. 608–611, 2001.

F. Granja, J. Morari, E. C. Morari, L. A. C. Correa, L. V. M. Assumpção, and L. S. Ward, “Proline homozygosity in codon 72 of p53 is a factor of susceptibility for thyroid cancer,” Cancer Letters, vol. 210, no. 2, pp. 151–157, 2004.

A. Storey, M. Thomas, A. Kalita et al., “Role of a p53 polymorphism in the development of human papilloma-virus-associated cancer,” Nature, vol. 393, no. 6682, pp. 229–234, 1998.

J. Lin, M. R. Spitz, Y. Wang et al., “Polymorphisms of folate metabolic genes and susceptibility to bladder cancer: A Case-Control Study,” Carcinogenesis, vol. 25, no. 9, pp. 1639–1647, 2004.

M. Manuguerra, G. Matullo, F. Veglia et al., “Multi-factor dimensionality reduction applied to a large prospective investigation on gene-gene and gene-environment interactions,” Carcinogenesis, vol. 28, no. 2, pp. 414–422, 2007.

D. W. Cai, X. F. Liu, R. G. Bu et al., “Genetic polymorphisms of MTHFR and aberrant promoter hypermethylation of the RASSF1A gene in bladder cancer risk in a chinese population,” Journal of International Medical Research, vol. 37, no. 6, pp. 1882–1889, 2009.

M. Wang, H. Zhu, G. Fu et al., “Polymorphisms of methylenetetrahydrofolate reductase and methionine synthase genes and bladder cancer risk: a case-control study with meta-analysis,” Clinical and Experimental Medicine, vol. 9, no. 1, pp. 9–19, 2009.

M. Izmirli, N. Inandiklioglu, D. Abat et al., “MTHFR gene polymorphisms in bladder cancer in the Turkish population,” Asian Pacific Journal of Cancer Prevention, vol. 12, no. 7, pp. 1833–1835, 2011.

W. You, Z. Li, C. Jing et al., “MTHFR C677T and A1298C polymorphisms were associated with bladder cancer risk and disease progression: a meta-analysis,” DNA and Cell Biology, vol. 32, no. 5, pp. 260–267, 2013.

W. Xu, H. Zhang, F. Wang, and H. Wang, “Quantitative assessment of the association between MHTFR C677T (rs1801133, Ala222Val) polymorphism and susceptibility to bladder cancer,” Diagnostic Pathology, vol. 8, no. 1, article 95, 2013.

S. Sanyal, F. Festa, S. Sakano et al., “Polymorphisms in DNA repair and metabolic genes in bladder cancer,” Carcinogenesis, vol. 25, no. 5, pp. 729–734, 2004.

M. R. Karagas, S. Park, H. H. Nelson et al., “Methylenetetrahydrofolate reductase (MTHFR) variants and bladder cancer: A Population-Based Case-Control Study,” International Journal of Hygiene and Environmental Health, vol. 208, no. 5, pp. 321–327, 2005.

S. Ouerhani, E. Oliveira, R. Marrakchi et al., “Methylenetetrahydrofolate reductase and methionine synthase polymorphisms and risk of bladder cancer in a Tunisian population,” Cancer Genetics and Cytogenetics, vol. 176, no. 1, pp. 48–53, 2007.

C.-J. Chung, Y.-S. Pu, C.-T. Su et al., “Polymorphisms in one-carbon metabolism pathway genes, urinary arsenic profile, and urothelial carcinoma,” Cancer Causes and Control, vol. 21, no. 10, pp. 1605–1613, 2010.

M. R. Safarinejad, N. Shafiei, and S. Safarinejad, “Genetic susceptibility of methylenetetrahydrofolate reductase (MTHFR) gene C677T, A1298C, and G1793A polymorphisms with risk for bladder transitional cell carcinoma in men,” Medical Oncology, vol. 28, no. 1, pp. S398–S412, 2011.

L. H. Kobeissi, I. A. Yassine, M. E. Jabbour, M. A. Moussa, and H. R. Dhaini, “Urinary bladder cancer risk factors: A Lebanese Case-Control Study,” Asian Pacific Journal of Cancer Prevention, vol. 14, no. 5, pp. 3205–3211, 2013.

K. Erdurak, P. E. Dundar, B. C. Ozyurt, E. Negri, C. La Vecchia, and Z. Tay, “Smoking, occupation, history of selected diseases and bladder cancer risk in Manisa, Turkey,” European Journal of Cancer Prevention, vol. 23, no. 1, pp. 58–61, 2014.

D. A. Bell, J. A. Taylor, D. F. Paulson, C. N. Robertson, J. L. Mohler, and G. W. Lucier, “Genetic risk and carcinogen exposure: a common inherited defect of the carcinogen-metabolism gene glutathione S-transferase M1 (GSTM1) that increases susceptibility to bladder cancer,” Journal of the National Cancer Institute, vol. 85, no. 14, pp. 1159–1164, 1993.

L. S. Engel, E. Taioli, R. Pfeiffer et al., “Pooled analysis and meta-analysis of glutathione S-transferase M1 and bladder cancer: a HuGE review,” American Journal of Epidemiology, vol. 156, no. 2, pp. 95–109, 2002.

M. García-Closas, N. Malats, D. Silverman et al., “NAT2 slow acetylation, GSTM1 null genotype, and risk of bladder cancer: results from the spanish bladder cancer study and meta-analyses,” Lancet, vol. 366, no. 9486, pp. 649–659, 2005.

L. E. Moore, D. R. Baris, J. D. Figueroa et al., “GSTM1 null and NAT2 slow acetylation genotypes, smoking intensity and bladder cancer risk: results from the New England bladder cancer study and NAT2 meta-analysis,” Carcinogenesis, vol. 32, no. 2, pp. 182–189, 2011.

N. Rothman, M. Garcia-Closas, N. Chatterjee et al., “A multi-stage genome-wide association study of bladder cancer identifies multiple susceptibility loci,” Nature Genetics, vol. 42, no. 11, pp. 978–984, 2010.

K. Golka, T. Reckwitz, M. Kempkes et al., “N-acetyltransferase 2 (NAT2) and glutathione S-transferase μ (GSTM1) in bladder-cancer patients in a highly industrialized area,” International Journal of Occupational and Environmental Health, vol. 3, no. 2, pp. 105–110, 1997.

K. J. Kiriluk, S. M. Prasad, A. R. Patel, G. D. Steinberg, and N. D. Smith, “Bladder cancer risk from occupational and environmental exposures,” Urologic Oncology: Seminars and Original Investigations, vol. 30, no. 2, pp. 199–211, 2012.

D. Ovsiannikov, S. Selinski, M.-L. Lehmann et al., “Polymorphic enzymes, urinary bladder cancer risk, and structural change in the local industry,” Journal of Toxicology and Environmental Health Part A: Current Issues, vol. 75, no. 8-10, pp. 557–565, 2012.

G. Matullo, A. Naccarati, and B. Pardini, “MicroRNA expression profiling in bladder cancer: the challenge of next-generation sequencing in tissues and biofluids,” International Journal of Cancer, vol. 138, no. 10, pp. 2334–2345, 2016.

J. S. Ross, K. Wang, R. N. Al-Rohil et al., “Advanced urothelial carcinoma: next-generation sequencing reveals diverse genomic alterations and targets of therapy,” Modern Pathology, vol. 27, no. 2, pp. 271–280, 2014.