High HSF4 expression is an independent indicator of poor overall survival and recurrence free survival in patients with primary colorectal cancer.

Heat shock factor 4 (HSF4) is a member of the HSF family. In this study, by using data from the Cancer Genome Atlas-Colorectal Cancer (TCGA-CRC), we investigated the expression profile and the prognostic value of the HSF4 in terms of overall survival (OS) and recurrence free survival (RFS) in CRC patients. RNA-Seq data showed that HSF4 RNA expression was significantly higher in CRC tissues (N = 380) than in the corresponding normal tissues (N = 51) (mean ± SD: 3.56 ± 1.28 vs. 1.85 ± 0.87, P < 0.0001). High HSF4 expression group had significantly higher ratio of stages III/IV patients (52/86, 60.5%) than low HSF4 expression group (110/264, 41.7%; P = 0.0024). Besides, the high HSF4 expression group also had significantly increased expression of CEA (CEA ≥ 5, 26/51, 51.0% vs. 64/186, 34.4%), higher proportion of recurrence (32/86, 37.2% vs. 48/254, 18.9%, P = 0.0005) and death (36/90, 40.0% vs. 49/277, 17.7%, P < 0.0001) compared with the low HSF4 expression group. Multivariate analysis confirmed that high HSF4 expression was an independent prognostic factor of poor OS (HR = 2.111, 95%CI: 1.350-3.302, P = 0.001) and RFS (HR = 1.958, 95%CI: 1.224-3.131, P = 0.005). Bioinformatic analysis showed that HSF4 can directly interact with DUSP26, ZBED8, and MAPK14. It is also coexpressed with PTGER1, COL11A2, CLPS, and ARSA and colocalized with PTGER1, ADRB1, PEX12, CLPS, PSEN2, KCNJ5, CPA1, ARSA, PNLIP, IRX4, CPA2, IDUA, BCKDHA, and CTRL. We hypothesized that HSF4 might exert its oncogenic effects in CRC via some of these genes.