Bilirubin is not associated with urinary bladder cancer risk and prognosis: A Mendelian Randomization Study in the UK Biobank ============================================================================================================================= * Nadezda Lipunova * Richard T Bryan * Maurice Zeegers ## ABSTRACT **Background** Mutations in *UGT1A* gene have been associated with the development and prognosis of urinary bladder cancer (UBC). UGT1A proteins are involved in a spectrum of detoxification processes, hence the biological mechanism between *UGT1A* and UBC is difficult to elucidate. Concurrently, mild hyperbilirubinemia, caused by alterations in *UGT1A*, has been associated with multiple health outcomes. We have investigated the potential effect of mild hyperbilirubinemia on UBC risk and prognosis, using a Mendelian Randomization (MR) approach in the UK Biobank. **Methods** Data on 1,281 UBC patients and 4,071 controls was available for a two-stage least squares MR estimation with rs6742078 as an instrumental variable. First, linear regression was fitted to establish the relationship between the rs6742078 and bilirubin levels (total and unconjugated). Secondly, bilirubin values were used to predict tested outcomes under a logistic model. Both stages were adjusted for participant sex, smoking status, and age. **Results** MR analysis showed no significant effects of bilirubin levels on UBC risk (total bilirubin: OR=1.02, 95% CI: 0.99-1.04; unconjugated bilirubin: OR=1.02, 95% CI: 0.99-1.05). No effects were observed for events of UBC recurrence, progression, or survival. **Conclusion** Our study suggests mild hyperbilirubinemia is not associated with urinary bladder cancer risk and prognosis. Keywords * bilirubin * bladder cancer * risk * prognosis * Mendelian Randomization ## INTRODUCTION Urinary bladder cancer (UBC) poses a great burden for patients and healthcare systems, but the clinical management has seen little change over the last few decades [1, 2], Genetic discoveries are a useful tool for investigating UBC development and enhance patient care. However, elucidating the biological mechanisms behind the effect of genetics that would help identify clinical targets offers a challenge. Previous studies have identified the relevance of single nucleotide polymorphisms (SNPs) in the *UGT1A* gene in UBC development, namely rs11892031 [3], rs17863783 [4], and rs7571337 [5]. Moreover, reports for *UGT1A* effects being modified by smoking [6] is especially relevant for the smoking population, as impaired metabolite clearance may result in a longer exposure to tobacco-related carcinogens. Importantly, the evidence is lacking whether *UGT1A* mutations are also relevant for UBC prognosis. One study assessing SNPs in *UGT1A* (rs11892031 and rs17863783) for UBC prognosis did not find any significant results [7]. Interestingly, a mutation in the *UGT1A*28* locus is a principal cause for Gilbert’s syndrome (GS). GS is a benign condition with global prevalence of about 5-10%, clinically present as mildly elevated bilirubinaemia [8, 9]. Among affected, the expression of the bilirubin UDP-glucuronosyltransferase (UGT) is reduced to about 30% of normal hepatic glucuronidating activity [10, 11]. Rs6742078, located close to the *UGT1A*28* TATA box polymorphism was found to explain roughly 18% of the total serum bilirubin and can be used as a proxy in genetic analyses [12], Mildly elevated bilirubinaemia has been associated with various health outcomes, specifically protection against heart disease, diabetes, and overall mortality [9, 13]. It is difficult to establish whether hyperbilirubinemia is the actual causal agent for these benefits or merely a correlated marker, but the use of Mendelian Randomization (MR) makes it possible to decipher its’ role in human health. MR studies on bilirubinaemia and cardiovascular disease [14], stroke [15], and non-fatty liver disease [16] showed bilirubin is unlikely to play a role in these traits. On the other hand, a report on diabetes risk indicates increased bilirubin levels exhibit a protective effect [17], It has to be noted most MR studies (with the exception of Luo *et al*. [16]) have only considered total serum bilirubin (TB) levels. Instead, high TB levels may reflect impaired liver function, and hence only unconjugated bilirubin (UB) is hypothesized to offer health benefits, if any [9]. Despite *UGT1A* being associated with bladder cancer, the potential effect of mild hyperbilirubinemia on UBC risk and outcomes have not yet been assessed. In the current MR study, we aimed to investigate the association between TB and UB serum levels and the risk of developing UBC, as well as UBC prognosis (recurrence, progression, and survival) in the UK Biobank cohort [18, 19]. ## METHODS ### Study population UK Biobank is one of the largest population-based cohorts in the UK. Genetic and clinical data is collected for over 500,000 participants, aged 40-69 at the time of recruitment (2006-2010). Detailed descriptions on UK Biobank design, data collection, and processing are described in detail elsewhere [18, 19]. For the analysis on UBC prognosis, the sample was limited to UBC patients only (International Classification of Diseases (ICD) codes of C67.0, C67.1, C67.2, C67.3, C67.4, C67.5, C67.6, C67.7, C67.8, C67.9, D09.0 (ICD10) and 1880, 1882, 1884, 1886, 1888, 1889, 2337 (ICD9). For analyses on UBC risk, a random sample of 5000 participants were drawn as controls. After data cleaning procedures, there were 1,281 UBC patients and 4,071 controls with clinical and genetic data available for analysis. ### SNP selection Rs6742078 was chosen as the instrumental variable (IV) for bilirubin levels. Rs6742078 shows high correlation (r2>0.8) with the causative locus for GS and has been successfully used for other MR studies on various phenotypes [15, 16, 17], ### Outcomes and exposures The overall risk of developing UBC was treated as a binary outcome. Death was modelled as an overall (death vs no death) or a UBC-specific event (death vs no death, when primary cause of death was assigned C67- (ICD10) or 188-related (ICD9) codes). Prognostic endpoints (recurrence, progression) were modelled using data from Hospital Episode Statistics (HES) and Death Registry data. For recurrence, three conditions to be representative of an event. First, a transurethral resection of a bladder tumour (TURBT) (OPCS4 code M42) is regarded to be enough to signify a UBC event. Secondly, a time period longer than 4 months between chemotherapeutic treatments into urinary bladder (OPCS4 codes M494/M495) was considered to be substantial to correspond to two independent events. Thirdly, a UBC diagnosis if an examination of the urinary bladder (OPCS4 code M45) was led by a relevant intervention within 6 months. Relevant interventions were chemotherapeutic treatments into urinary bladder, cystectomy, radiotherapy, and chemotherapy (OPCS4 codes of M494/M495; M34; X65; X72; X292, X298, X308, X352, respectively). UBC progression was considered to have taken place if either a TURBT (OPCS4 code M42) or examination of the urinary bladder (OPCS4 code M45) was followed by interventions of cystectomy (OPCS4 code M34) and/or radiotherapy (OPCS4 code X65) within 6 months. Two recurrence and/or progression events were considered independent of one another if time in between the records was greater than 3 months. Bilirubin levels (micromoles per litre, umol/L) have been measured in the UK Biobank by photometric colour analysis on a Beckman Coulter AU5800. ### Ethics and consent All UK Biobank participants have provided informed consent. Current research has been conducted using the UK Biobank Resource under Application Number 42772. ### Genotype data quality control (QC) and imputation Detailed QC procedures on imputation in the UK Biobank are described elsewhere [19]. Our analyses were restricted White British participants to avoid population stratification bias [19]. ### Statistical analysis F statistic was calculated to establish the instrument strength for our analyses [20]. The rule of thumb is to consider an instrument acceptable if the F statistic is >10 [20]. Firstly, have tested the association between TB and UB levels with the following outcomes: UBC risk, death, recurrence, and progression. For all analyses, logistic regression was fitted while adjusting for participant sex, smoking status (ever vs never), and age at the time of diagnosis (or recruitment for controls). Secondly, we have used a two-stage least squares (2SLS) MR estimation on all outcomes with rs6742078 as an IV. In the first stage, linear regression was fitted to establish the relationship between the IV and bilirubin levels (separately for TB and UB) assuming an additive model of inheritance. The value of bilirubin from the first step was used to predict UBC outcomes in the second stage under a logistic model. Both stages were adjusted for participant sex, smoking status (ever vs never), and age. The difference between baseline characteristics between cases and controls was tested using Chi-square for categorical variables and t-test for continuous outcomes. Statistical significance was assumed at two-tailed p<0.05 throughout the analysis. ## RESULTS The rs6742078 proved to be a valid instrument for the analyses, as the F statistic has been calculated as 1,701 for TB and 1,751 for UB (both p<0.001). Baseline characteristics of the cases and controls are presented in **Table 1**. UBC cases were significantly older (mean age=61.3 years) than the controls (mean age=57.4 years), and were more often male (cases=81.2%, controls=49.8%) as well as having history of smoking (cases=68.4%, controls=45.7%). However, bilirubin levels and rs6742028 status were comparable between bladder cancer patients and controls (TB among cases=9.92 umol/L, TB among controls=9.66 umol/L, p=0.07; Chi-square p value for independence of rs6742028 genotype between cases and controls was 0.52). View this table: [Table 1.](http://medrxiv.org/content/early/2020/08/14/2020.08.13.20174102/T1) Table 1. Descriptive characteristic of bladder cancer patients and controls in the UK Biobank. For the subgroup analyses within the UBC cases, there were 512 recurrence, 50 progression, and 176 death events. Non-muscle-invasive bladder cancer (NMIBC) comprised a total of 1,181 out of 1,281 UBC cases. A multivariable-adjusted model for testing the effect of TB on overall UBC risk resulted in an odds ratio (OR) of 0.99 (95% confidence interval (CI): 0.98-1.01, p=0.43) **(Table 2)** The estimate was equivalent for UB (OR=0.99, 95% CI:0.97-1.01, p=0.36). For the events of UBC recurrence, progression, and death, both TB and UB have indicated neutral effects with all ORs being close to 1 **(Table 2)**. View this table: [Table 2.](http://medrxiv.org/content/early/2020/08/14/2020.08.13.20174102/T2) Table 2. Mendelian randomization analysis for bilirubin and urinary bladder cancer risk and prognosis in the UK Biobank. In the MR analysis on UBC risk, causal OR estimates were identical for TB and UB levels (TB: 1.02, 95% CI: 0.99-1.04, p=0.18; UB: 1.02, 95% CI: 0.99-1.05, p=0.18) **(Table 2)**. In the additional analyses on clinically-relevant phenotypes of NMIBC recurrence, progression, and muscle-invasive bladder cancer (MIBC) death, no association was observed in either observational model or MR analysis (data not shown). ## DISCUSSION In the current study, we have investigated the causal role of TB and UB in the development of bladder cancer, as well as UBC prognosis. Rs6742028 showed to be a valid instrument for estimating causality, however no association was observed with any of the tested outcomes. Mutations in the *UGT1A* cluster, located in 2q37.1, have been associated with UBC in previous genetic studies [3, 4, 21, 22], SNPs in *UGT1A* have shown both harmful [21] and protective effects for the risk of developing *de novo* UBC [4, 22], but there is uncertainty if *UGT1A* influences the disease course, especially in NMIBC. Lacombe *et al* evaluated the effect of *UGT1A1*28* SNP *rs8175347* in a group of 189 NMIBC cases [23] and concluded that it significantly increased the risk of recurrence, with a gene-smoking interaction present (higher recurrence risk for ever smokers). Nevertheless, the many strata in the study resulted in small subgroup sample sizes [23]; therefore, the findings must be replicated in other cohorts to establish validity. The UGTs, coded by the region, are involved in various processes of detoxification. Until now, the most common hypothesis suggests that variations in the activity of UGTs may confer to an altered speed and efficiency of carcinogen clearance, and thus have a role in cancer development[24]. In fact, a mapping study of the T allele of rs17863783 has been shown it results in an increased mRNA expression of the UGT1A6 in human liver tissue cells[4]. As such, it is hypothesized increased UGT1A6 function helps to facilitate faster removal of carcinogens from the bladder epithelium[4, 6, 24], However, the UGTs may have more than one mode of action, due to a variety of function. Given that mild hyperbilirubinemia has been implicated to play a role in many conditions, it was useful to investigate whether the alterations in *UGT1A*28*, resulting high blood concentrations of bilirubin may also have an effect on UBC development and prognosis. However, our study suggests such effects are unlikely and tobacco-related carcinogen detoxification remains the most likely hypothesis for explaining the effect of *UGT1A* mutations on UBC. We acknowledge our study carries important limitations. First, prognostic events have been modelled using an algorithm that makes use of HES, but none of the events have been confirmed on a case-by-case basis. Moreover, some events had low sample sizes (e.g. progression), and likely underpowered to detect an existing association. ## CONCLUSIONS To summarise, we have investigated the potential role of mild hyperbilirubinemia in UBC risk and prognosis. Using an MR methods for estimating causal relationships, our study suggests *UGT1A* locus contributes to bladder cancer via mechanisms other than elevated bilirubin. ## Data Availability UK Biobank data are available through a procedure described at http://www.ukbiobank.ac.uk/using-the-resource/. [http://www.ukbiobank.ac.uk/using-the-resource/](http://www.ukbiobank.ac.uk/using-the-resource/) ## Additional information ### Ethics Approval and consent to participate All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments, or comparable ethical standards. Informed Consent: Informed consent was obtained from all individual participants included in the study. ## Data availability UK Biobank data are available through a procedure described at [http://www.ukbiobank.ac.uk/using-the-resource/](http://www.ukbiobank.ac.uk/using-the-resource/). ### Conflict of Interest The authors declare they have no conflict of interest. ### Funding No funding to declare. ### Author Contribution Statement NL designed and carried out the study and wrote the main manuscript; RBT, MZ, and NL have reviewed the manuscript. * Received August 13, 2020. * Revision received August 13, 2020. * Accepted August 14, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## REFERENCES 1. 1.Park JC, Hahn NM. Bladder cancer: a disease ripe for major advances. Clinical advances in hematology & oncology : H&O. 2014;12(12):838–45. 2. 2.Berdik C. Unlocking bladder cancer. Nature. 2017;551(7679):S34-s5. 3. 3.Rothman N, Garcia-Closas M, Chatterjee N, Malats N, Wu X, Figueroa JD, et al. A multi-stage genome-wide association study of bladder cancer identifies multiple susceptibility loci. Nature genetics. 2010;42(11):978–84. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.687&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20972438&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F14%2F2020.08.13.20174102.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000283540500015&link_type=ISI) 4. 4.Tang W, Fu YP, Figueroa JD, Malats N, Garcia-Closas M, Chatterjee N, et al. Mapping of the UGT1A locus identifies an uncommon coding variant that affects mRNA expression and protects from bladder cancer. Hum Mol Genet. 2012;21(8):1918–30. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddr619&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22228101&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F14%2F2020.08.13.20174102.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000302302400020&link_type=ISI) 5. 5.Wang J, Wu X, Kamat A, Barton Grossman H, Dinney CP, Lin J. Fluid intake, genetic variants of UDP-glucuronosyltransferases, and bladder cancer risk. British journal of cancer. 2013;108(11):2372–80. 6. 6.Garcia-Closas M, Rothman N, Figueroa JD, Prokunina-Olsson L, Han SS, Baris D, et al. Common genetic polymorphisms modify the effect of smoking on absolute risk of bladder cancer. Cancer Res. 2013;73(7):2211–20. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FucmVzIjtzOjU6InJlc2lkIjtzOjk6IjczLzcvMjIxMSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA4LzE0LzIwMjAuMDguMTMuMjAxNzQxMDIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 7. 7.Grotenhuis AJ, Dudek AM, Verhaegh GW, Witjes JA, Aben KK, van der Marel SL, et al. Prognostic relevance of urinary bladder cancer susceptibility loci. PloS one. 2014;9(2):e89164. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0089164&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24586564&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F14%2F2020.08.13.20174102.atom) 8. 8.Hirschfield GM, Alexander GJ. Gilbert’s syndrome: an overview for clinical biochemists. Annals of clinical biochemistry. 2006;43(Pt 5):340-3. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1258/000456306778520034&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17022875&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F14%2F2020.08.13.20174102.atom) 9. 9.Wagner KH, Shiels RG, Lang CA, Seyed Khoei N, Bulmer AC. Diagnostic criteria and contributors to Gilbert’s syndrome. Critical reviews in clinical laboratory sciences. 2018;55(2):129–39. 10. 10.Bosma PJ, Chowdhury JR, Bakker C, Gantla S, de Boer A, Oostra BA, et al. The genetic basis of the reduced expression of bilirubin UDP-glucuronosyltransferase 1 in Gilbert’s syndrome. The New England journal of medicine. 1995;333(18):1171–5. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJM199511023331802&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7565971&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F14%2F2020.08.13.20174102.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995TB56000002&link_type=ISI) 11. 11.Sugatani J. Function, genetic polymorphism, and transcriptional regulation of human UDP-glucuronosyltransferase (UGT) 1A1. Drug metabolism and pharmacokinetics. 2013;28(2):83–92. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2133/dmpk.DMPK-12-RV-096&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23089802&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F14%2F2020.08.13.20174102.atom) 12. 12.Johnson AD, Kavousi M, Smith AV, Chen MH, Dehghan A, Aspelund T, et al. Genome- wide association meta-analysis for total serum bilirubin levels. Human Molecular Genetics. 2009;18(14):2700–10. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddp202&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19414484&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F14%2F2020.08.13.20174102.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000267441200017&link_type=ISI) 13. 13.Wagner KH, Wallner M, Molzer C, Gazzin S, Bulmer AC, Tiribelli C, et al. Looking to the horizon: the role of bilirubin in the development and prevention of age-related chronic diseases. Clinical science (London, England : 1979). 2015;129(1):1–25. 14. 14.McArdle PF, Whitcomb BW, Tanner K, Mitchell BD, Shuldiner AR, Parsa A. Association between bilirubin and cardiovascular disease risk factors: using Mendelian randomization to assess causal inference. BMC cardiovascular disorders. 2012;12:16. 15. 15.Lee SJ, Jee YH, Jung KJ, Hong S, Shin ES, Jee SH. Bilirubin and Stroke Risk Using a Mendelian Randomization Design. Stroke. 2017;48(5):1154–60. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToic3Ryb2tlYWhhIjtzOjU6InJlc2lkIjtzOjk6IjQ4LzUvMTE1NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA4LzE0LzIwMjAuMDguMTMuMjAxNzQxMDIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 16. 16.Luo L, An P, Jia X, Yue X, Zheng S, Liu S, et al. Genetically Regulated Bilirubin and Risk of Non-alcoholic Fatty Liver Disease: A Mendelian Randomization Study. Frontiers in Genetics. 2018;9:662. 17. 17.Abbasi A, Deetman PE, Corpeleijn E, Gansevoort RT, Gans RO, Hillege HL, et al. Bilirubin as a potential causal factor in type 2 diabetes risk: a Mendelian randomization study. Diabetes. 2015;64(4): 1459-69. 18. 18.Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS medicine. 2015;12(3):e1001779. 19. 19.Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018;562(7726):203–9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-018-0579-z&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30305743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F14%2F2020.08.13.20174102.atom) 20. 20.Davies NM, Holmes MV, Davey Smith G. Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians. BMJ (Clinical research ed). 2018;362:k601. [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE2OiIzNjIvanVsMTJfMi9rNjAxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDgvMTQvMjAyMC4wOC4xMy4yMDE3NDEwMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 21. 21.Selinski S, Lehmann ML, Blaszkewicz M, Ovsiannikov D, Moormann O, Guballa C, et al. Rs11892031[A] on chromosome 2q37 in an intronic region of the UGT1A locus is associated with urinary bladder cancer risk. Arch Toxicol. 2012;86(9):1369–78. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00204-012-0854-y&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22532026&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F14%2F2020.08.13.20174102.atom) 22. 22.Zhang Y, Sun Y, Chen T, Hu H, Xie W, Qiao Z, et al. Genetic variations rs11892031 and rs401681 are associated with bladder cancer risk in a Chinese population. Int J Mol Sci. 2014;15(11): 19330-41. 23. 23.Lacombe L, Fradet V, Levesque E, Pouliot F, Larue H, Bergeron A, et al. Phase II Drug- Metabolizing Polymorphisms and Smoking Predict Recurrence of Non-Muscle-Invasive Bladder Cancer: A Gene-Smoking Interaction. Cancer prevention research (Philadelphia, Pa). 2016;9(2):189–95. 24. 24.Sundararaghavan VL, Sindhwani P, Hinds TD, , Jr. Glucuronidation and UGT isozymes in bladder: new targets for the treatment of uroepithelial carcinomas? Oncotarget. 2017;8(2):3640–8.