Abstract
Background
Sarcomatoid differentiation of RCC (sRCC) is associated with poor survival. Recent studies showed marked response of sRCC to immune checkpoint blockade (ICB). While distinctive patterns of gene expression in sRCC have been identified, the gene regulatory programs and TFs that drive SD remain unknown. The aim of this study is to nominate TFs responsible for SD and to investigate their association with the clinical outcomes of patients with RCC.
Methods
Chromatin immunoprecipitation and sequencing (ChIP-seq) for H3K27ac – a histone modification associated with active regulatory elements – was performed on pathologically reviewed sRCC and non-sRCC samples collected at the Dana-Farber Cancer Institute. Regulatory elements that were differentially active between the two groups were identified based on levels of H3K27ac (Benjamini-Hochberg q<0.01, log-fold change (LFC) threshold=1). Enrichment of specific TF binding motifs at activated regulatory elements in sRCC was assessed using HOMER. Differential gene expression analysis of TFs was performed using DESeq2 on RNA-seq data from TCGA. A Mann-Whitney U test was performed on RNA-seq data from the IMmotion151 and Javelin Renal 101 clinical trials to compare mean expression level of TFs in transcript per million (TPM) in these trials. Patients with non-sRCC enrolled in the IMmotion151 trials were divided into quartiles based on gene expression levels of candidate TFs. Progression-free survival (PFS) was compared between non-sRCC patients stratified by expression quartiles as well as patients with sRCC using a multivariable Cox proportional-hazards model accounting for age and IMDC risk score. To validate these findings, a similar analysis was performed in the Javelin Renal 101 trial.
Results
We obtained high-quality H3K27ac ChIP-seq profiles for 9 sRCC and 17 non-sRCC samples. We identified 278 candidate regulatory elements with increased H3K27ac levels in sRCC vs. non-sRCC. These regulatory elements were enriched for nucleotide motifs bound by the TFs FOSL1 and E2F7. Differential expression analysis between 48 sRCC vs. 493 non-sRCC samples showed that FOSL1 and E2F7 were significantly upregulated in sRCC vs. non-sRCC (LFC=1.7, q=5e-11; LFC=1.8, q=1.3e-20; resp.). Mean TPMs of FOSL1 and E2F7 were significantly increased in sRCC vs. non-sRCC in IMmotion151 cohort and Javelin Renal 101 (all p<0.001). Among patients who received sunitinib, those with the highest quartile of FOSL1 and E2F7 expression showed significantly shorter PFS in IMmotion151 patients compared to patients with the lowest quartile of expression (HR=1.6, 95%CI=1.3-2.2, p=0.008 & HR=2.6, 95%CI=1.8-3.7, p<0.001; resp.). Furthermore, patients with highest quartile of expression showed similar PFS compared to patients with sRCC (p=0.56 and p=0.64; resp.). These results were validated in the sunitinib arm of the Javelin Renal 101 cohort (Figure).
Figure: Kaplan-Meier curves of progression-free survival (PFS) in patients with in the sunitinib arm of Javelin Renal 101 by sarcomatoid differentiation and FOSL1 expression levels.
Conclusions
This is the first study to characterize the epigenomic landscape of sRCC by integrating ChIP-seq and RNA-seq data. Our findings implicated FOSL1 and E2F7 as transcriptional regulators of SD with prognostic relevance. These TFs seems to be associated with aggressive behavior in non-sRCC as well. Further studies are underway to functionally validate these results.