Discovering urinary bladder cancer risk variants: Status quo after almost ten years of genome-wide association studies

Bladder cancer course, four genetic high-risk variants, and histopathological findings

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.

N6-methyladenosine-related single-nucleotide polymorphism analyses identify oncogene RNFT2 in bladder cancer

Background

Single-nucleotide polymorphisms (SNPs) in N6-methyladenosine (m6A) related genetic locus play significant roles in tumorigenesis and development. The expression level of many oncogenes and tumour suppressor genes changed because of m6A-associated SNPs. In addition, the relationship between m6A-SNP and bladder cancer (BCa) has not been well studied.

Methods

We screened m6A-SNPs in BCa by combining m6A-SNPs data and GWAS-SNPs data. Expression quantitative trait loci (eQTL) and differential expression gene (DEGs) analyses were performed. In ring finger protein, transmembrane 2 (RNFT2), rs3088107 (C  > G) was found to have significant eQTL signals and make RNFT2 gene differentially-regulated mostly in BCa. We validated the expression level of RNFT2 in 32 pairs of BCa tissues and eight BCa cell lines by quantitative real-time PCR (qRT-PCR). Functional assays were performed to investigate the role of rs3088107 and RNFT2 in BCa in vitro.

Results

We identified 673 m6A-SNPs, which were associated with BCa. Of these m6A-SNPs, 221 showed eQTL signals, amongst which, rs3088107 in RNFT2 showed significant eQTL signals. Results of bioinformatic analyses showed that 11 genes with m6A-SNPs had a differential expression level in BCa. RNFT2 was predicted to be significantly up-regulated in BCa. The qRT-PCR results validated that RNFT2 was highly expressed in our own BCa tissues and cell lines. High expression of RNFT2 also indicated a worse overall survival. We also revealed that rs3088107 (C  > G) could inhibit the expression and m6A modification of RNFT2 by qRT-PCR, western-blot and m6A-RIP assays. Moreover, the results of functional assays indicated that RNFT2 promoted BCa cell proliferation and migration.

Conclusion

This research found that m6A-SNPs were associated with oncogene RNFT2 in BCa. Furthermore, m6A-SNPs showed great application potential as a new BCa diagnostic biomarker and prognostic indicator.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-022-02701-z.

Development of a custom next-generation sequencing panel for the determination of bladder cancer risk in a Tunisian cohort

BACKGOUND

Bladder cancer (BCa) is a heterogeneous disease caused by the interaction between environmental and genetic risk factors. The goal of this case-control study was to evaluate the implication of a selected SNP panel in the risk of BCa development in a Tunisian cohort. We were also interested in studying the interaction between this predictive panel and environmental risk factors.

METHODS

The case/control cohort was composed with 249 BCa cases and 255 controls. The designed Bladder cancer hereditary panel (BCHP) was composed of 139 selected variants. These variants were genotyped by an amplification-based targeted Next-Generation Sequencing (NGS) on the Ion Torrent Proton sequencer (Life Technologies, Ion Torrent technology).

RESULTS

We have found that rs162555, rs2228000, rs10936599, rs710521, rs3752645, rs804276, rs4639, rs4881400 and rs288980 were significantly associated with decreased risk of bladder cancer. However the homozygous genotypes for VPS37C (rs7104333, A/A), MPG (rs1013358, C/C) genes or the heterozygous genotype for ARNT gene (rs1889740, rs2228099, rs2256355, rs2864873), GSTA4 (rs17614751) and APOBR/IL27 (rs17855750) were significantly associated with increased risk of bladder cancer development compared to reference group (OR 2.53, 2.34, 1.99, 2.00, 2.00, 1.47, 1.96 and 2.27 respectively). We have also found that non-smokers patients harboring heterozygous genotypes for ARNT/rs2864873 (A > G), ARNT/ rs1889740 (C > T) or GSTA4/rs17614751 (G-A) were respectively at 2.775, 3.069 and 6.608-fold increased risk of Bca development compared to non-smokers controls with wild genotypes. Moreover the ARNT CT (rs1889740), ARNT CG (rs2228099), ARNT TC (rs2864873) and GSS GA genotypes were associated with an increased risk of BCa even in absence of professional risk factors. Finally the decision-tree analysis produced a three major BCa classes. These three classes were essentially characterized by an intensity of tobacco use more than 20 pack years (PY) and the CYP1A2 (rs762551) genotype.

CONCLUSIONS

The determined association between environmental factors, genetic variations and the risk of Bca development may provide additional information to urologists that may help them for clinical assessment and treatment decisions. Nevertheless, the underlying mechanisms through which these genes or SNPs affect the clinical behavior of BCas require further studies.

Genetic variants in N6-methyladenosine are associated with bladder cancer risk in the Chinese population

Recently N⁶-Methyladenosine (m⁶A) has been identified to guide the interaction of RNA-binding protein hnRNP C and their target RNAs, which is termed as m⁶A-switches. We systematically investigated the association between genetic variants in m⁶A-switches and bladder cancer risk. A two-stage case-control study was performed to systematically calculate the association of single nucleotide polymorphisms (SNPs) in 2798 m⁶A-switches with bladder cancer risk in 3,997 subjects. A logistic regression model was used to assess the effects of SNPs on bladder cancer risk. A series of experiments were adopted to explore the role of genetic variants of m⁶A-switches. We identified that rs5746136 (G > A) of SOD2 in m⁶A-switches was significantly associated with the reduced risk of bladder cancer (additive model in discovery stage: OR = 0.80, 95% CI 0.69-0.93, P = 3.6 × 10^-3; validation stage: adjusted OR = 0.88, 95% CI 0.79-0.99, P = 3.0 × 10^-2; combined analysis: adjusted OR = 0.85, 95% CI 0.78-0.93, P = 4.0 × 10^-4). The mRNA level of SOD2 was remarkably lower in bladder cancer tissues than the paired adjacent samples. SNP rs5746136 may affect m⁶A modification and regulate SOD2 expression by guiding the binding of hnRNP C to SOD2, which played a critical tumor suppressor role in bladder cancer cells by promoting cell apoptosis and inhibiting proliferation, migration and invasion. In conclusion, our findings suggest the important role of genetic variants in m⁶A modification. SOD2 polymorphisms may influence the expression of SOD2 via an m⁶A-hnRNP C-dependent mechanism and be promising predictors of bladder cancer risk.

Association of Glutathione S-transferase gene polymorphism with bladder Cancer susceptibility

BACKGROUND

We conducted a meta-analysis to evaluate the relationship between the glutathione S-transferase μ1 (GSTM1)- and glutathione S-transferase θ1 (GSTT1)- null genotypes and susceptibility to bladder cancer.

METHODS

We identified association reports from the databases of PubMed, Embase, the Cochrane Library and the China Biological Medicine Database (CBM disc) on July 1, 2017 and synthesized eligible investigations. Results were expressed using odds ratios (ORs) for dichotomous data, and we also calculated 95% confidence intervals (CIs).

RESULTS

In this meta-analysis, we found that the GSTM1-null genotype was associated with bladder cancer risk in the overall population, and individually in whites, Africans and Asians (overall population: OR = 1.40, 95% CI: 1.31-1.48, P<0.00001; whites: OR = 1.39, 95% CI: 1.26-1.54, P<0.00001; Africans: OR = 1.54, 95% CI: 1.16-2.05, P = 0.003; Asians: OR = 1.45, 95% CI: 1.33-1.59, P<0.00001). The GSTT1-null genotype was associated with bladder cancer risk in the overall population, but not in whites, in Africans or Asians (overall population: OR = 1.11, 95% CI: 1.01-1.22, P = 0.03; whites: OR = 1.16, 95% CI: 0.99-1.36, P = 0.07; Africans: OR = 1.07, 95% CI: 0.65-1.76, P = 0.79; Asians: OR = 1.05, 95% CI: 0.91-1.22, P = 0.51). Interestingly, a dual-null GSTM1-GSTT1 genotype was associated with bladder cancer risk in the overall population and in Asians (overall population: OR = 1.48, 95% CI: 1.15-1.92, P = 0.002; Asians: OR = 1.62, 95% CI: 1.15-2.28, P = 0.006). In conclusion, the GSTM1-null, GSTT1-null and dual-null GSTM1-GSTT1 genotypes might be associated with the onset of bladder cancer, but additional genetic-epidemiological studies should be conducted to explore this association further.