Abstract 4739: Development of a prognostic ovarian cancer gene signature based on glucocorticoid receptor transcriptional activity
Authors: Ananth V. Panchamukhi, Lei Huang, Eleanore P. Gray, Yuanhang Liu, Manisha Taya, Jennifer T. Veneris, Xiaonan Hou, S John Weroha, Balázs Győrffy, Gini F. Fleming, Chen Wang, Suzanne D. Conzen
Journal: Cancer Research (2025)

Abstract Background: Epithelial ovarian cancer (EOC), commonly known as ovarian cancer (OvCa), is the deadliest gynecological malignancy. High tumor expression of either NR3C1 mRNA (encoding the glucocorticoid receptor or GR) or GR protein (IHC) are associated with poor overall survival (OS) in patients (pts) with advanced stage OvCa. We hypothesized that high GR expression and associated GR transcriptional activity favor OvCa cell survival and chemotherapy resistance, thereby reducing OS. Here we investigated whether a GR transcriptome-based gene signature (GR-sig) would be a better prognosticator than tumor NR3C1 mRNA expression alone. Methods: mRNAs up- or downregulated &amp;gt;1.5 fold after glucocorticoid treatment (Dex) in one or both cultured OvCa cell lines (OVSAHO and HEYA8) were defined as differentially expressed genes (DEGs). DEGs were then identified in n=1247 unique OvCa microarrays from 15 publicly available Gene Expression Omnibus (GEO) datasets with matched pt OS (training set). We then used the OS data to train a Cox proportional hazards model. The “best subset selection” (BeSS) algorithm and bootstrap sampling identified those GR-regulated DEGs whose tumor mRNA expression correlated most significantly with pt OS, creating a GR-sig “risk score.” Subsequently, all n=522 microarray and n=374 RNA-seq publicly available samples from the OvCa TCGA were used as two independent validation sets. Microarray and RNA-seq data were normalized using the MAS5 and FPKM-UQ algorithms respectively. The two validation sets were independently stratified relative to the median training set risk score into low- and high-risk score groups. To examine the prognostic ability of NR3C1 mRNA expression alone, the median training set NR3C1 expression value was used to stratify validation set pts into low- and high-expression groups. In all cases, OS was estimated by the Kaplan-Meier method. Results: We identified n=2288 protein-coding DEGs in one or both cell lines. Of these genes, n=1251 were expressed in each GEO microarray sample (training set). A 10-gene GR-signature (GR-Sig) was then identified using the BeSS algorithm and bootstrap sampling, yielding a median training set risk score of 1.02. In both the microarray and RNA-seq validation datasets, the 10-gene GR-sig risk score was associated with significant OS differences between low- and high-risk score groups (microarray p = 0.018; RNA-seq p = 0.00049). In contrast, the OS difference between low- and high- NR3C1 mRNA expression groups did not reach significance (microarray p = 0.12; RNA-seq p = 0.28). Conclusion: In this retrospective analysis, a 10-gene GR-sig derived from OvCa GR-regulated DEGs performed better as a prognosticator than NR3C1 mRNA expression alone. Prospective application of the GR-sig may be useful in stratifying patients for treatment response with novel agents in OvCa clinical trials. Citation Format: Ananth V. Panchamukhi, Lei Huang, Eleanore P. Gray, Yuanhang Liu, Manisha Taya, Jennifer T. Veneris, Xiaonan Hou, S John Weroha, Balázs Győrffy, Gini F. Fleming, Chen Wang, Suzanne D. Conzen. Development of a prognostic ovarian cancer gene signature based on glucocorticoid receptor transcriptional activity [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 4739.