ABSTRACT Background The expression of CYP19A1 has implications for the prognosis of female bladder cancer. However, this study aimed to explore the association between single nucleotide polymorphisms (SNPs) in CYP19A1 and bladder cancer risk, as no prior research has addressed this association. Research design and methods We selected and genotyped five CYP19A1 SNPs (rs4646, rs6493487, rs1062033, rs17601876, and rs3751599) in 217 patients and 550 controls using the Agena MassARRAY system. Logistic regression analysis was employed to calculate the odds ratio (OR) and 95% confidence intervals (CIs). Bioinformatics predicted SNP functions and CYP19A1 involving pathways. Results Our study revealed a significant association between bladder cancer risk and four SNPs (rs4646 (AC vs. CC: OR = 1.71, FDR-p = 0.005), rs6493487 (G vs. A: OR = 0.68, FDR-p = 0.011), rs1062033 (G vs. C: OR = 0.36, FDR-p < 0.001), and rs17601876 (GA vs. GG: OR = 1.66, FDR-p = 0.008)) in CYP19A1. The three SNPs (rs4646, rs1062033, and rs17601876) were significantly correlated with CYP19A1 expression levels in normal whole blood (p < 0.05). Moreover, CYP19A1 was found to primarily participate in the steroid hormone biosynthesis and metabolic pathways. Conclusions Consequently, CYP19A1 gene polymorphisms may play a crucial role in the genetic susceptibility to bladder cancer.