RNA-Seq Uncovers Association of Endocrine-Disrupting Chemicals with Hub Genes and Transcription Factors in Aggressive Prostate Cancer

This study analyzes publicly available RNA-seq data to comprehensively include the complex heterogeneity of prostate cancer (PCa) etiology. It combines prostate and prostate cancer (PCa) cell lines, representing primary PCa cells, Gleason scores, ages, and PCa of different racial origins. Additionally, some cell lines were exposed to endocrine-disrupting chemicals (EDCs). The research aims to identify hub genes and transcription factors (TFs) of the prostate carcinogenesis pathway as molecular targets for clinical investigations to elucidate EDC-induced aggressiveness and to develop potential biomarkers for their exposure risk assessments. PCa cells rely on androgen receptor (AR)-mediated signaling to survive, develop, and function. Fifteen various RNA-seq datasets were normalized for distribution, and the significance (<i>p</i>-value < 0.05) threshold of differentially expressed genes (DEGs) was set based on |log2FC| ≥ 2 change. Through integrated bioinformatics, we applied cBioPortal, UCSC-Xena, TIMER2.0, and TRRUST platforms, among others, to associate hub genes and their TFs based on their biologically meaningful roles in aggressive prostate carcinogenesis. Among all RNA-Seq datasets, we found 75 overlapping DEGs, with BUB1B (32%) and CCNB1 (29%) genes exhibiting the highest degree of mutation, amplification, and deletion. EDC-associated CCNB1, BUB1B, and CCNA2 in PCa cells exposed to EDCs were consistently shown to be associated with high Gleason scores (≥4 + 3) and in the >60 age group of patients. Selected TFs (E2F4, MYC, and YBX1) were also significantly associated with DEGs (NCAPG, MKI67, CCNA2, CCNB1, CDK1, CCNB2, AURKA, UBE2C, BUB1B) and influenced the overall survival (<i>p</i>-value < 0.05) of PCa cases. This is one of the first comprehensive studies combining 15 publicly available RNA-seq datasets to demonstrate the association of EDC-associated hub genes and their TFs aligning with the aggressive carcinogenic pathways in the higher age group (>60 years) of patients. The findings highlight the potential of these hub genes as candidates for further studies to develop molecular biomarkers for assessing the EDC-related PCa risk, diagnosing PCa aggressiveness, and identifying therapeutic targets.