Identification of bicalutamide resistance-related genes and prognosis prediction in patients with prostate cancer

Background: Prostate cancer (PCa) is the second most common type of cancer and the fifth leading cause of cancer-related death in men. Androgen deprivation therapy (ADT) has become the first-line therapy for inhibiting PCa progression; however, nearly all patients receiving ADT eventually progress to castrate-resistant prostate cancer. Therefore, this study aimed to identify hub genes related to bicalutamide resistance in PCa and provide new insights into endocrine therapy resistance.

Methods: The data were obtained from public databases. Weighted correlation network analysis was used to identify the gene modules related to bicalutamide resistance, and the relationship between the samples and disease-free survival was analyzed. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses were performed, and hub genes were identified. The LASSO algorithm was used to develop a bicalutamide resistance prognostic model in patients with PCa, which was then verified. Finally, we analyzed the tumor mutational heterogeneity and immune microenvironment in both groups.

Results: Two drug resistance gene modules were identified. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses revealed that both modules are involved in RNA splicing. The protein–protein interaction network identified 10 hub genes in the brown module LUC7L3, SNRNP70, PRPF3, LUC7L, CLASRP, CLK1, CLK2, U2AF1L4, NXF1, and THOC1) and 13 in the yellow module (PNN, PPWD1, SRRM2, DHX35, DMTF1, SALL4, MTA1, HDAC7, PHC1, ACIN1, HNRNPH1, DDX17, and HDAC6). The prognostic model composed of RNF207, REC8, DFNB59, HOXA2, EPOR, PILRB, LSMEM1, TCIRG1, ABTB1, ZNF276, ZNF540, and DPY19L2 could effectively predict patient prognosis. Genomic analysis revealed that the high- and low-risk groups had different mutation maps. Immune infiltration analysis showed a statistically significant difference in immune infiltration between the high- and low-risk groups, and that the high-risk group may benefit from immunotherapy.

Conclusion: In this study, bicalutamide resistance genes and hub genes were identified in PCa, a risk model for predicting the prognosis of patients with PCa was constructed, and the tumor mutation heterogeneity and immune infiltration in high- and low-risk groups were analyzed. These findings offer new insights into ADT resistance targets and prognostic prediction in patients with PCa.