Introduction (original) (raw)

EXCLI J EXCLI Journal 1611-2156 Leibniz Research Centre for Working Environment and Human Factors 2023-5862 10.17179/excli2023-5862 Doc867 Original article Bladder cancer course, four genetic high-risk variants, and histopathological findings Kadhum Thura 1 2 Selinski Silvia 1 Blaszkewicz Meinolf 1 Reinders Jörg 1 Roth Emanuel 3 Volkert Frank 3 Ovsiannikov Daniel 4 Moormann Oliver 4 Gerullis Holger 5 Barski Dimitri 5 Otto Thomas * 5 Höhne Svetlana 1 Hengstler Jan G. 1 Golka Klaus 1 1Leibniz Research Centre for Working Environment and Human Factors at TU Dortmund (IfADo), Dortmund, Germany 2Specialist Clinic for Psychosomatic Rehabilitation, Mittelrhein-Klinik, Boppard - Bad Salzig, Germany 3Department of Urology, Evangelic Hospital, Paul-Gerhardt Foundation, Lutherstadt Wittenberg, Germany 4St.-Josefs-Hospital Dortmund-Hoerde, Dortmund, Germany 5Rheinland Klinikum Lukaskrankenhaus Neuss, Neuss, Germany *To whom correspondence should be addressed: Thomas Otto, Rheinland Klinikum Lukaskrankenhaus Neuss, Neuss, Germany; Tel: +49 (0)2131 888-2401, E-mail: thomas.otto@rheinlandklinikum.de 18 08 2023 2023 22 867 879 23 01 2023 07 08 2023 Copyright © 2023 Kadhum et al. 2023

This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (http://creativecommons.org/licenses/by/4.0/) You are free to copy, distribute and transmit the work, provided the original author and source are credited.

This article is available from https://www.excli.de/vol22/excli2023-5862.pdf

Urinary bladder cancer, a smoking and occupation related disease, was subject of several genome-wide association studies (GWAS). However, studies on the course of the disease based on GWAS findings differentiating between muscle invasive bladder cancer (MIBC) and non-muscle invasive bladder cancer (NMIBC) are rare. Thus we investigated 4 single nucleotide polymorphisms (SNPs) detected in GWAS, related to the genes coding for TACC3 (transforming, acidic coiled-coil containing protein 3), for FGFR3 (fibroblast growth factor receptor 3), for PSCA (prostate stem cell antigen) and the genes coding for CBX6 (chromobox homolog 6) and APOBEC3A (apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like 3A). This study is based on 712 bladder cancer patients and 875 controls from 3 different case control studies in Germany. The 4 SNPs of interest (PSCA rs2294008 and rs2978974, FGFR3-TACC3 rs798766, and CBX6-APOBEC3A rs1014971) were determined by real-time polymerase chain reaction. The distribution of the 4 SNPs does not vary significantly between cases and controls in the entire study group and in the 3 local subgroups, including two former highly industrialized areas and a region without such history. Also, no significant differences in the bladder cancer subgroups of MIBC and NMIBC were observed. The 4 investigated SNPs do not noticeably contribute differently to the bladder cancer risk for the bladder cancer subgroups of MIBC and NMIBC.

muscle invasive bladder cancer (MIBC) non-muscle invasive bladder cancer (NMIBC) PSCA gene rs2294008 and rs2978974 FGFR3-TACC3 gene region rs798766 CBX6-APOBEC3A gene region rs1014971 Introduction

Bladder cancer (BC) is the 4th common cancer in males and the 14th common cancer in females in Germany with an estimated 22.430 men and 7.100 women newly diagnosed in 2014 (Robert Koch Institute, 2018[38]). World-wide, it is the 6th common cancer in males and the 18th common cancer in females with an estimated 440.864 men and 132.414 women newly diagnosed in 2020 (Sung et al., 2021[45]). This smoking and occupation related disease was a subject of several genome-wide association studies (GWAS; overview in Selinski, 2017[42]; de Maturana et al., 2018[5]). Furthermore, a considerable number of studies investigating the impact of genetically based risk factors on bladder cancer was performed (overview: de Maturana et al., 2018[5]; Golka et al., 2011[11]). Regarding the impact on the course of the disease, multiple loci have been identified, aiming at different endpoints, e.g., T category, grade, recurrence, and mortality yet demonstrating little consensus (Lipunova et al., 2019[28]).

Thus, we investigated 4 single nucleotide polymorphisms (SNPs) detected in GWAS (Kiemeney et al., 2010[18]; Rothman et al., 2010[40]; Wu et al., 2009[52]; Fu et al., 2012[7]) related to the region of the genes coding for TACC3 and for FGFR3 (transforming, acidic coiled-coil containing protein - fibroblast growth factor receptor 3) located on 4p16.3, the gene coding for PSCA (prostate stem cell antigen) located on 8q24.3, and the region of the genes coding for CBX6 (chromobox homolog 6) and for APOBEC3A (apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like 3A) located on 22q13.1. Furthermore, T category and grading in muscle invasive bladder cancer (MIBC) and non-muscle invasive bladder cancer (NMIBC) and the course of the disease, i.e., the number of the relapses and their histopathological findings, were observed.

The region of the genes coding for TACC3 (transforming, acidic coiled-coil containing protein 3) and for FGFR3 (fibroblast growth factor receptor 3) is located on 4p16.3, a known Huntington disease region (Lee et al., 2012[23]). The reference SNP rs798766(T) is associated with bladder cancer risk (Kiemeney et al., 2010[18]; Meng et al., 2017[32]). rs798766 is located in an intron of TACC3, 70 kb from FGFR3 (Kiemeney et al., 2010[18]). FGFR3 belongs to the family of tyrosine kinases. FGFR3 signaling is involved in development, differentiation, cell survival, migration, angiogenesis, and carcinogenesis (Turner and Grose, 2010[46]). FGFR3 mutations in bladder cancer are mostly activating, followed by gene rearrangements and amplification (Cancer Genome Atlas Research Network, 2014[1]; Helsten et al., 2016[13]). They are mainly detected in genetically stable bladder cancer (van Rhijn et al., 2003[47]) and have been associated with progression in bladder cancer (Parker et al., 2014[36]). FGFR3 gene rearrangements generate constitutively activated and oncogenic FGFR3 kinase protein products, and cellular dependence on these drivers confers sensitivity to selective FGFR inhibition (Williams et al., 2013[51]; Wu et al., 2013[53]; Kim et al., 2018[20]). FGFR alterations are targeted in clinical routine already by the kinase inhibitor erdafitinib, which was approved by FDA in April 2019 (Loriot et al., 2019[30]). Most recently, it was shown in a phase 3 study that erdafitinib is superior to chemotherapy in the therapy of advanced or metastatic bladder cancer in patients with unresectable, advanced or metastatic urothelial cancer and select FGFR3/2 receptor alterations (mutations/fusions) (Loriot et al., 2023[29]) According to the TCGA data base, 15.43 % (29 out of 188) bladder cancer cases provided somatic mutations in their tumors (NIH, 2023[34]).

PSCA (prostate stem cell antigen) was initially identified as a prostate-specific cell-surface antigen (Reiter et al., 1998[37]), but later found to be expressed in many human tissues (Fu et al., 2012[7]). The gene coding for PSCA is located on 8q24.3. This gene encodes a glycosylphosphatidylinositol-anchored cell membrane glycoprotein. It is highly expressed in the prostate and, to lesser extent, in the bladder, placenta, colon, kidney, and stomach. It is up-regulated in many prostate cancers and expressed in cancers of the bladder and pancreas. Its polymorphisms result in an upstream start codon in some individuals (Weizmann Institute of Science, 2022[50]). rs2978974 is located 10 kb upstream of rs2294008, within an alternative untranslated first exon of PSCA (Fu et al., 2012[7]). The variant rs2294008(T) (Wu et al., 2009[52]; Fu et al., 2012[7]; Deng et al., 2019[6]; Cui et al., 2019[4]) is associated with bladder cancer risk. In a recent meta-analysis, rs2294008 was associated with bladder cancer (OR = 1.15, 95 % CI 1.11-1.18) (Cui et al., 2019[4]). Furthermore, this SNP was strongly associated with gastric cancer (OR = 1.32, 95 % CI 1.27-1.39). rs2978974 was significant, but to a clearly lower extent, associated with bladder cancer (OR 1.09; 95 % CI 1.03-1.15) (Cui et al., 2019[4]). It is noteworthy, that Fu et al. (2012[7]) reported a significant multiplicative interaction between these two SNPs (P = 0.035). The non-risk allele G of rs2978974 showed strong interaction with nuclear proteins from five cell lines tested, implying a regulatory function. Thus, Fu et al. (2012[7]) concluded that a joint effect of the two PSCA SNPs may be important for bladder cancer susceptibility. According to the TCGA data base, 0.46 % (1 out of 216) bladder cancer cases provided somatic mutations in their tumors (NIH, 2023[34]).

The gene region coding for CBX6 (chromobox homolog 6) and for APOBEC3A (apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like 3A) is located on 22q13.1. This locus is located in a non-genic region approximately 25 kb centromeric of APOBEC3A and 64 kb telomeric of CBX6 (Rothman et al., 2010[40]). APOBEC is a protein family encoded by eleven genes. Subtypes of APOBEC3 can cause specific mutations in RNA and DNA at distinct preferred nucleotide contexts in human cancer (Cao and Wu, 2018[2]). APOBEC enzymes are a major source of mutation in bladder cancer. Tumors enriched for APOBEC mutagenesis have better survival (Glaser et al., 2017[8]). CBX6 is a protein coding gene. It is a component of a Polycomb Group (PcG) multiprotein PRC1 (polycomb repressive complex 1)-like complex, a complex class required to maintain the transcriptionally repressive state of many genes, including Hox genes, throughout development (Vandamme et al., 2011[48]). PcG PRC1 complex acts via chromatin remodeling and modification of histones; it mediates monoubiquitination of histone H2A 'Lys-119', rendering chromatin heritably changed in its expressibility (Weizmann Institute of Science, 2022[49]). The variant rs1014971(T) is associated with bladder cancer risk (de Maturana et al., 2018[5]; Rothman et al., 2010[40]). CBX gene variants are also associated with risk of other cancers (Li et al., 2020[25]; Lin et al., 2020[26]; Xie et al., 2020[54]). According to the TCGA data base, 0.74 % (3 out of 408) bladder cancer cases provided somatic mutations in APOBEC3A. Likewise, 0.74 % (3 out of 408) bladder cancer cases provided somatic mutations in CBX6 (NIH, 2023[34]).

Materials and Methods

This study is based on 712 bladder cancer patients and 875 controls from 3 different case control studies: Evangelic Hospital, Paul-Gerhardt Foundation Lutherstadt Wittenberg, Lutherstadt Wittenberg, Germany (203 cases and 210 controls (Roth et al., 2012[39])), St.-Josefs-Hospital Dortmund-Hoerde, Dortmund, Germany (156 cases and 237 controls (Ovsiannikov et al., 2012[35])), and Rheinland Klinikum Lukaskrankenhaus Neuss, Neuss, Germany (353 cases and 428 controls (Höhne et al., 2017[14])). Controls were from the same respective hospitals without a known history of cancer. Enrolment of patients was as follows: Lutherstadt Wittenberg from December 1995 to January 1999, Dortmund from July 2009 to December 2010, and Neuss from June 2009 to November 2011. Follow-up derived from medical records was as follows: Lutherstadt Wittenberg from September 2008 to June 2009, Dortmund from May 2012 to August 2012, Neuss from July 2013 to February 2014 (Selinski et al. 2016[44]; Roth et al., 2012[39]).

Analytical methods

DNA of cases and controls was isolated from venous blood samples according to standard methods in the Dortmund Institute (Selinski et al., 2017[43]). The 4 SNPs of interest (PSCA rs2294008 and rs2978974, gene region FGFR3-TACC3 rs798766, and gene region CBX6-APOBEC3A rs1014971) were detected by real-time polymerase chain reaction (rt-PCR) using a TaqMan® assay (Saravana Devi et al., 2008[41]; Selinski et al., 2017[43]). The missings for the detected SNPs for cases and controls respectively were as follows: 0 and 0 for rs2294008, 2 and 1 for rs2978974, 0 and 0 for rs798766, and 0 and 4 for rs1014971.

Statistical methods

Hardy-Weinberg Equilibrium was checked using a chi-squared test with one degree of freedom. Chi-squared tests were used testing differences between categorial variables, logistic regression was used to compute Likelihood Ratio and Wald tests as well as odds ratios (OR) and 95 % confidence intervals (95 % CI). The analysis was conducted searching for differences between cases vs controls, recurrent vs non-recurrent and MIBC vs NMIBC and stratified for gender, invasiveness, and smoking habits. False discovery rate (FDR) was used for multiple testing correction.

The studies in the 3 departments of urology were approved by the ethics committee of the Dortmund Institute (24/2009, 31/2009, 50/2011, 74/2014) and by the Institutional Review Board.

Results

The total study group contains 3 subgroups from different regions in Germany with comparable characteristics regarding gender and age among cases as well as among controls. In bladder cancer patients, the portions of the described characteristics were as follows: Total study group: 79 % male, median age 70.5 years. Wittenberg: 86 % male, median age 65.4 years. Dortmund: 75 % male, median age 71.0 years. Neuss: 77 % male, median age 73.2 years. Detailed information, including that of the controls is presented in Table 1(Tab. 1) and 2(Tab. 2).

A description of the clinical parameters of the study group is presented in Table 3a, b(Tab. 3). As expected, in the combined group of muscle invasive bladder cancer (MIBC) and non-muscle invasive bladder cancer (NMIBC), significant differences were observed regarding the smoking habits. A total of 16 % of the male patients, in contrast to 49 % of the female patients were non-smokers (p <0.0001). The portions of the different T categories did not differ significantly between male and female patients.

Male patients showed a higher percentage of NMIBC (82 % vs 79 %), whereas female patients showed a higher percentage of MIBC (21 % vs 18 %). Grading differed significantly between the genders in 2 subgroups (Table 3b(Tab. 3)). There was no significant gender-related difference in the distribution of NMIBC and MIBC (Table 3a(Tab. 3)). A stratification of the clinical parameters for NMIBC and MIBC is presented in Table 3b(Tab. 3).

The frequency of the 4 investigated SNPs in the combined group of NMIBC and MIBC cases and in controls, additionally stratified for the 3 local subgroups, is presented in Table 4(Tab. 4).

The distribution of the 4 investigated SNPs, stratified for NMIBC and MIBC cases does not vary significantly between the 3 local subgroups, including two formerly highly industrialized areas with industries associated with elevated bladder cancer risk (Dortmund (Ovsiannikov et al., 2012[35]) and Lutherstadt Wittenberg (Roth et al., 2012[39])) and a region with no industries with known bladder cancer risk (Neuss (Höhne et al., 2017[14])). No significant differences between cases and controls and between the subgroups of NMIBC and MIBC bladder cancer were observed (Table 5(Tab. 5)).

The subgroup of cases with progression stratified for the histopathological findings for T category and grading are presented in Table 6(Tab. 6). Due to the small number of cases with progression, a statistical evaluation of the impact of the 4 investigated SNPs was not performed.

Discussion

Since the first papers on the impact of polymorphism on the metabolism of the antitubercular agent isoniazid (INH) in the 1950s (Hughes et al., 1954[15]; Mitchell et al., 1957[33]), a considerable number of studies have been published on the impact of genetic polymorphisms on bladder cancer risk (overview: Golka et al., 2011[11]). However, with the first paper on genome-wide association studies on bladder cancer risk (Kiemeney et al., 2008[19]), a new era started, showing the impact of genetics on bladder cancer of loci, mostly not associated with bladder cancer before (Golka et al., 2011[11]).

To date, about 40 single nucleotide polymorphisms are shown to be associated with bladder cancer risk, mostly within the range between 1.1 and 1.5, mostly identified and/or confirmed in GWAS (Selinski, 2017[43]; de Maturana et al., 2018[5]; Koutros et al., 2023[22]). This is much lower than the risk conferred by the portion of N-acetyltransferase 2 (NAT2) slow acetylators reported in the first studies on bladder cancer risk due to high exposure to highly carcinogenic aromatic amines in the chemical industry at that time (Lewalter and Miksche, 1992[24]; Cartwright et al., 1982[3]; Golka et al., 2002[10]) or conferred by the portion of GSTM1 negatives reported formerly in highly industrialized areas (Golka et al., 1998[9]; Hung et al., 2004[16]).

The first study on variants detected in GWAS on bladder cancer course was published by Roth et al. (2012[39]), showing an association of longer relapse-free survival in GSTT1 positives (HR = 0.60, 95 % CI = 0.42-0.87) in the formerly highly industrialized area of Lutherstadt Wittenberg. GSTM1 negatives tended to be associated with a better prognosis, whereas the NAT2 risk alleles rs710521(A) and rs9642880(T), detected in the first GWAS on bladder cancer (Kiemeney et al., 2008[19]) tended to be associated with a poor prognosis. Since then, only few studies on the impact of genetic polymorphisms detected in GWAS on the course of the bladder cancer disease have been conducted (Grotenhuis et al., 2014[12]; Lipunova et al., 2019[27]), most possibly due to the enormous logistic effort. Suitable candidates for an impact on the course of the disease are those genetic polymorphisms and/or loci, which have a proven influence on bladder cancer risk (Roth et al., 2012[39]). Another promising approach is the investigation of the effect of SNP combinations or, even better, a large panel of SNPs as recently applied for the determination of the risk to contract bladder cancer (Koutros et al., 2023[22]). Although our group has experience with the impact of SNP combinations on bladder cancer risk regarding smoking habits (Selinski et al., 2017[43]), we decided not to apply this approach for this study because the available case numbers are by far too small when two or more polymorphisms are combined.

It is now established for a long time that bladder cancer is no more one tumor entity, but two: Non-muscle invasive bladder cancer (NMIBC) and muscle-invasive bladder cancer (MIBC) (Knowles and Hurst, 2015[21]). It may be noted, that genetic analysis has resulted in a new biological classification for MIBC discriminating 6 different subtypes (Kamoun et al., 2020[17]). Similar activities for NMIBC are undergoing (Marzouka et al., 2022[31]). Due to the very different course of these two tumor entities, we decided to investigate the impact on the subgroups of NMIBC and the clinically much more important MIBC separately in larger study group now including in addition to Roth et al. (2012[39]), patients from two other studies (Ovsiannikov et al., 2012[35]; Höhne et al., 2017[14]). The 4 investigated SNPs showed no noticeable differences in the frequency between the NMIBC and MIBC subgroups. This is in line with the findings of Grotenhuis et al. (2014[12]) for PSCA gene related rs2294008, TACC3 - FGFR gene region related rs798766, and CBX6 - APOBEC3A gene region related rs1014971. This may be due to the minimal impact of a single SNP. Thus, in future studies an impact of SNP combinations should be investigated on large study groups, as this procedure had successfully shown that a 4 SNP combination confers increased bladder cancer risk particularly in never smokers (Selinski et al., 2017[43]). Recently, this approach was applied by a GWAS using a 24 SNP marker panel (Koutros et al., 2023[22]). Another promising approach would be to investigate the impact of selected SNPs and/or their combinations in the clinically relevant group of symptomatic patients showing micro and/or gross hematuria.

DeclarationAuthor contributions

Conceptualization T.K., S.S., T.O., J.G.H., and K.G., methodology S.S., M.B., and J.R., software S.S., laboratory analysis M.B., J.R., data analysis T.K., S.S., investigation T.K., S.S., T.O., and K.G., resources F.V., O.M., T.O., and J.G.H., data curation E.R., D.O., H.G., D.B., S.H., and S.S., writing-original draft preparation T.K., S.S., M.B., and K.G., writing-review and editing T.K., T.O., and K.G,. visualization T.K., K.G., supervision F.V., O.M., T.O., and J.G.H, project administration J.G.H., K.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional review board statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board and by the Ethics Committee of Leibniz Research Centre for Working Environment and Human Factors at TU Dortmund (IfADo).

Informed consent statement

Informed consent was obtained from all subjects involved in the study.

Data availability statement

The data are available on reasonable request from the corresponding author.

Conflicts of interest

The authors declare no conflict of interest.

1 Cancer Genome Atlas Research Network Comprehensive molecular characterization of urothelial bladder carcinoma Nature2014 507 7492 315 322 10.1038/nature12965 2 Cao W Wu W Apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like gene expression, RNA editing, and microRNAs regulation Methods Mol Biol2018 1699 75 81 10.1007/978-1-4939-7435-1_5 3 Cartwright RA Glashan RW Rogers HJ Ahmad RA Barham-Hall D Higgins E Role of N-acetyltransferase phenotypes in bladder carcinogenesis: a pharmacogenetic epidemiological approach to bladder cancer Lancet1982 2 8303 842 845 10.1016/s0140-6736(82)90810-8 4 Cui H Tang M Zhang M Liu S Chen S Zeng Z Variants in the PSCA gene associated with risk of cancer and nonneoplastic diseases: systematic research synopsis, meta-analysis and epidemiological evidence Carcinogenesis2019 40 70 83 10.1093/carcin/bgy151 5 de Maturana EL Rava M Anumudu C Sáez O Alonso D Malats N Bladder cancer genetic susceptibility. A systematic review Bladder Cancer2018 4 215 226 10.3233/BLC-170159 6 Deng S Ren ZJ Jin T Yang B Dong Q Contribution of prostate stem cell antigen variation rs2294008 to the risk of bladder cancer Medicine (Baltimore)2019 98 e15179 10.1097/MD.0000000000015179 7 Fu YP Kohaar I Rothman N Earl J Figueroa JD Ye Y Common genetic variants in the PSCA gene influence gene expression and bladder cancer risk Proc Natl Acad Sci U S A2012 109 4974 4979 10.1073/pnas.1202189109 8 Glaser AP Fantini D Wang Y Yu Y Rimar KJ Podojil JR APOBEC-mediated mutagenesis in urothelial carcinoma is associated with improved survival, mutations in DNA damage response genes, and immune response Oncotarget2017 9 4537 4538 10.18632/oncotarget.23344 9 Golka K Bandel T Schlaefke S Reich SE Reckwitz T Urfer W Urothelial cancer of the bladder in an area of former coal, iron, and steel industries in Germany: a case-control study Int J Occup Environ Health1998 4 79 84 10.1179/oeh.1998.4.2.79 10 Golka K Prior V Blaszkewicz M Bolt HM The enhanced bladder cancer susceptibility of NAT2 slow acetylators towards aromatic amines: a review considering ethnic differences Toxicol Lett2002 128 229 241 10.1016/s0378-4274(01)00544-6 11 Golka K Selinski S Lehmann ML Blaszkewicz M Marchan R Ickstadt K Genetic variants in urinary bladder cancer: collective power of the "wimp SNPs" Arch Toxicol2011 85 539 554 10.1007/s00204-011-0676-3 12 Grotenhuis AJ Dudek AM Verhaegh GW Witjes JA Aben KK van der Marel SL Prognostic relevance of urinary bladder cancer susceptibility loci PLoS One2014 One PLoS 9 e89164 10.1371/journal.pone.0089164 13 Helsten T Elkin S Arthur E Tomson BN Carter J Kurzrock R The FGFR landscape in cancer: analysis of 4.853 tumors by next-generation sequencing Clin Cancer Res2016 22 259 267 10.1158/1078-0432.CCR-14-3212 14 Höhne S Gerullis H Blaszkewicz M Selinski S Hengstler JG Otto T N-acetyltransferase 1*10 genotype in bladder cancer patients J Toxicol Environ Health A2017 80 417 422 10.1080/10937404.2017.1304727 15 Hughes HB Biehl JP Jones AP Schmidt LH Metabolism of isoniazid in man as related to the occurrence of peripheral neuritis Am Rev Tuberc1954 70 266 273 10.1164/art.1954.70.2.266 16 Hung RJ Boffetta P Brennan P Malaveille C Hautefeuille A Donato F GST, NAT, SULT1A1, CYP1B1 genetic polymorphisms, interactions with environmental exposures and bladder cancer risk in a high-risk population Int J Cancer2004 110 598 604 10.1002/ijc.20157 17 Kamoun A de Reyniès A Allory Y Sjödahl G Robertson AG Seiler R A consensus molecular classification of muscle-invasive bladder cancer Eur Urol2020 77 420 433 10.1016/j.eururo.2019.09.006 18 Kiemeney LA Sulem P Besenbacher S Vermeulen SH Sigurdsson A Thorleifsson G A sequence variant at 4p16.3 confers susceptibility to urinary bladder cancer Nat Genet2010 42 415 419 10.1038/ng.558 19 Kiemeney LA Thorlacius S Sulem P Geller F Aben KK Stacey SN Sequence variant on 8q24 confers susceptibility to urinary bladder cancer Nat Genet2008 40 1307 1312 10.1038/ng.229 20 Kim YS Kim K Kwon GY Lee SJ Park SH Fibroblast growth factor receptor 3 (FGFR3) aberrations in muscle-invasive urothelial carcinoma BMC Urol2018 18 68 10.1186/s12894-018-0380-1 21 Knowles MA Hurst CD Molecular biology of bladder cancer: new insights into pathogenesis and clinical diversity Nat Rev Cancer2015 15 25 41 10.1038/nrc3817 22 Koutros S Kiemeney LA Pal Choudhury P Milne RL Lopez de Maturana E Ye Y Genome-wide association study of bladder cancer reveals new biological and translational insights Eur Urol2023 84 127 137 10.1016/j.eururo.2023.04.020 23 Lee JM Gillis T Mysore JS Ramos EM Myers RH Hayden MR Common SNP-based haplotype analysis of the 4p16.3 Huntington disease gene region Am J Hum Genet2012 90 434 444 10.1016/j.ajhg.2012.01.005 24 Lewalter J Miksche L Empfehlungen zur arbeitsmedizinischen Prävention expositions- und dispositionsbedingter Arbeitsstoff-Beanspruchungen Verh Dt Ges Arbeitsmed1992 31 135–9 25 Li X Gou J Li H Yang X Bioinformatic analysis of the expression and prognostic value of chromobox family proteins in human breast cancer Sci Rep2020 10 17739 10.1038/s41598-020-74792-5 26 Lin K Zhu J Hu C Bu F Luo C Zhu X Comprehensive analysis of the prognosis for chromobox family in gastric cancer J Gastrointest Oncol2020 11 932 951 10.21037/jgo-20-208 27 Lipunova N Wesselius A Cheng KK van Schooten FJ Bryan RT Cazier JB Genome-wide association study for tumour stage, grade, size, and age at diagnosis of non-muscle-invasive bladder cancer Eur Urol Oncol2019 2 381 389 10.1016/j.euo.2018.08.020 28 Lipunova N Wesselius A Cheng KK van Schooten FJ Cazier JB Bryan RT Systematic review: genetic associations for prognostic factors of urinary bladder cancer Biomark Cancer2019 11 1179299X19897255 10.1177/1179299X19897255 29 Loriot Y Matsubara N Park SH Huddart RA Burgess EF Houede N Phase 3 THOR study: Results of erdafitinib (erda) versus chemotherapy (chemo) in patients (pts) with advanced or metastatic urothelial cancer (mUC) with select fibroblast growth factor receptor alterations (FGFRalt) J Clin Oncol2023 41 17, Suppl LBA4619 10.1200/JCO.2023.41.17_suppl.LBA4619 30 Loriot Y Necchi A Park SH Garcia-Donas J Huddart R Burgess E Erdafitinib in locally advanced or metastatic urothelial carcinoma N Engl J Med2019 381 338 348 10.1056/NEJMoa1817323 31 Marzouka NA Eriksson P Bernardo C Hurst CD Knowles MA Sjödahl G The Lund Molecular Taxonomy applied to non-muscle-invasive urothelial carcinoma J Mol Diagn2022 24 992 1008 10.1016/j.jmoldx.2022.05.006 32 Meng XY Shi MJ Chen JF Liao Y Hu BW Hireche A Association between the TACC3 rs798766 polymorphism and risk of urinary bladder cancer: a synthesis based on current evidence Dis Markers2017 2017 7850708 10.1155/2017/7850708 33 Mitchell R Bell JC Clinical implications of isoniazid blood levels in pulmonary tuberculosis N Engl J Med1957 257 1066 1070 10.1056/NEJM195711282572202 34 NIH, National Cancer Institute GDC data portal 2023 12 June 2023 Available from: https://portal.gdc.cancer.gov/ 35 Ovsiannikov D Selinski S Lehmann ML Blaszkewicz M Moormann O Haenel MW Polymorphic enzymes, urinary bladder cancer risk, and structural change in the local industry J Toxicol Environ Health A2012 75 557 565 10.1080/15287394.2012.675308 36 Parker BC Engels M Annala M Zhang W Emergence of FGFR family gene fusions as therapeutic targets in a wide spectrum of solid tumours J Pathol2014 232 4–15 10.1002/path.4297 37 Reiter RE Gu Z Watabe T Thomas G Szigeti K Davis E Prostate stem cell antigen: a cell surface marker overexpressed in prostate cancer Proc Natl Acad Sci U S A1998 95 1735 1740 10.1073/pnas.95.4.1735 38 Robert Koch Institute and the Association of Population-based Cancer Registries in Germany Cancer in Germany 2013/20142018 11th ed Berlin Robert Koch Institute 104 107 10.25646/5889 39 Roth E Selinski S Schikowsky C Seidel T Volkert F Blaszkewicz M Bladder cancer survival in a former industrial area in Saxony-Anhalt, Germany J Toxicol Environ Health A2012 75 1216 1225 10.1080/15287394.2012.709168 40 Rothman N Garcia-Closas M Chatterjee N Malats N Wu X Figueroa JD A multi-stage genome-wide association study of bladder cancer identifies multiple susceptibility loci Nat Genet2010 42 978 984 10.1038/ng.687 41 Saravana Devi S Vinayagamoorthy N Agrawal M Biswas A Biswas R Naoghare P Distribution of detoxifying genes polymorphism in Maharastrian population of central India Chemosphere2008 70 1835–9 10.1016/j.chemosphere.2007.08.008 42 Selinski S Discovering urinary bladder cancer risk variants: Status quo after almost ten years of genome-wide association studies EXCLI J2017 16 1288 1296 10.17179/excli2017-1000 43 Selinski S Blaszkewicz M Ickstadt K Gerullis H Otto T Roth E Identification and replication of the interplay of four genetic high-risk variants for urinary bladder cancer Carcinogenesis2017 38 1167 1179 10.1093/carcin/bgx102 44 Selinski S Bürger H Blaszkewicz M Otto T Volkert F Moormann O Occupational risk factors for relapse-free survival in bladder cancer patients J Toxicol Environ Health A2016 79 1136 1143 10.1080/15287394.2016.1219606 45 Sung H Ferlay J Siegel RL Laversanne M Soerjomataram I Jemal A Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries CA Cancer J Clin2021 71 209 249 10.3322/caac.21660 46 Turner N Grose R Fibroblast growth factor signalling: from development to cancer Nat Rev Cancer2010 10 116 129 10.1038/nrc2780 47 van Rhijn BW Vis AN van der Kwast TH Kirkels WJ Radvanyi F Ooms EC Molecular grading of urothelial cell carcinoma with fibroblast growth factor receptor 3 and MIB-1 is superior to pathologic grade for the prediction of clinical outcome J Clin Oncol2003 21 1912 1921 10.1200/JCO.2003.05.073 48 Vandamme J Völkel P Rosnoblet C Le Faou P Angrand PO Interaction proteomics analysis of polycomb proteins defines distinct PRC1 complexes in mammalian cells Mol Cell Proteomics2011 10 M110 10.1074/mcp.M110.002642 49 Weizmann Institute of Science Gene Cards. The Human Gene Database. CBX6 Gene (Protein Coding) Chromobox 6 2022 27 December 2022 Available from: https://www.genecards.org/cgi-bin/carddisp.pl?gene=CBX6 50 Weizmann Institute of Science GeneCards. The Human Gene Database. PSCA Gene. Prostate Stem Cell Antigen 2022 27 December 2022 Available from: https://www.genecards.org/cgi-bin/carddisp.pl?gene=PSCA 51 Williams SV Hurst CD Knowles MA Oncogenic FGFR3 gene fusions in bladder cancer Hum Mol Genet2013 22 795–803 10.1093/hmg/dds486 52 Wu X Ye Y Kiemeney LA Sulem P Rafnar T Matullo G Genetic variation in the prostate stem cell antigen gene PSCA confers susceptibility to urinary bladder cancer Nat Genet2009 41 991 995 10.1038/ng.421 53 Wu YM Su F Kalyana-Sundaram S Khazanov N Ateeq B Cao X Identification of targetable FGFR gene fusions in diverse cancers Cancer Discov2013 3 636–47 10.1158/2159-8290.CD-13-0050 54 Xie X Ning Y Long J Wang H Chen X Diverse CBX family members as potential prognostic biomarkers in non-small-cell lung cancer FEBS Open Bio2020 10 2206 2215 10.1002/2211-5463.12971 Table 1 Age distribution in the total study group of bladder cancer cases and the subgroups of non-muscle invasive bladder cancer (NMIBC) and muscle-invasive bladder cancer (MIBC), and controls Table 2 Gender distribution in the total study group of bladder cancer cases and the subgroups of non-muscle invasive bladder cancer (NMIBC) and muscle-invasive bladder cancer (MIBC), and controls Table 3 Description of the total study group of bladder cancer (BC) cases (a) and the subgroups of non-muscle invasive bladder cancer (NMIBC) and muscle-invasive bladder cancer (MIBC), stratified for gender, smoking habits, T category, and grading (b). Table 4 Frequency of the double mutated (MM), heterozygous (MN) and double native (NN) status of the 4 investigated single nucleotide polymorphisms (SNPs) in the total group of bladder cancer cases and in controls stratified for the 3 local subgroups Table 5 Frequency of double mutated (MM), heterozygous (MN) and double native (NN) status of the 4 investigated single nucleotide polymorphisms (SNPs) in bladder cancer (BC) cases, controls and in the BC cases stratified for non-muscle invasive (NMIBC) and muscle invasive (MIBC) bladder cancer Table 6 T category, staging and grading at progression in the investigated bladder cancer cases