Identification of low-frequency variants of UGT1A3 associated with bladder cancer risk by next-generation sequencing

Although genome-wide association studies (GWASs) have successfully revealed many common risk variants for bladder cancer, the heritability is still largely unexplained. We hypothesized that low-frequency variants involved in bladder cancer risk could reveal the unexplained heritability. Next-generation sequencing of 113 patients and 118 controls was conducted on 81 genes/regions of known bladder cancer GWAS loci. A two-stage validation comprising 3,350 cases and 4,005 controls was performed to evaluate the effects of low-frequency variants on bladder cancer risk. Biological experiments and techniques, including electrophoretic mobility shift assays, CRISPR/Cas9, RNA-Seq, and bioinformatics approaches, were performed to assess the potential functions of low-frequency variants. The low-frequency variant rs28898617 was located in the first exon of UGT1A3 and was significantly associated with increased bladder cancer risk (odds ratio = 1.50, P = 3.10 × 10^-6). Intriguingly, rs28898617 was only observed in the Asian population, but monomorphism was observed in the European population. The risk-associated G allele of rs28898617 increased UGT1A3 expression, facilitated UGT1A3 transcriptional activity, and enhanced the binding activity. In addition, UGT1A3 deletion significantly inhibited the proliferation, invasion, and migration of bladder cancer cells and xenograft tumor growth. Mechanistically, UGT1A3 induced LAMC2 expression by binding CBP and promoting histone acetylation, which remarkably promoted the progression of bladder cancer. This is the first targeted sequencing study to reveal that the novel low-frequency variant rs28898617 and its associated gene UGT1A3 are involved in bladder cancer development, providing new insights into the genetic architecture of bladder cancer.