Abstract Current data using gene expression and DNA modification profiles yield limited information for evaluating the potential resistance to therapies. Recent studies show that the chromatin landscape is dynamic and defines cell identity and disease status. Thus, changes in open chromatin sites could provide important information on disease progression and identify potential therapeutic targets. The majority of breast cancers express estrogen receptor (ER) and initially respond to endocrine therapies blocking ER activity. However, endocrine therapy resistance (ETR) occurs de novo or follows an initial response. Recent studies have identified new drivers of ETR, such as mTOR and CDK4/6. Therapy with these inhibitors in combination with endocrine therapy have improved patient survival, but the molecular mechanism of ETR is not fully understood and additional therapeutic options for ETR are needed. To investigate the molecular mechanisms of ETR, we established in vitro Long-Term Estrogen Deprivation (LTED) model using multiple human ER-positive breast cancer cell lines. These cells showed differential time-dependent patterns of ESR1 expression upon LTED with ESR1 upregulation in some cell lines and downregulation in others. LTEDs with ESR1 upregulation were sensitive to fulvestrant, while the others were resistant. We analyzed time-course of chromatin landscape transitions during acquisition of ETR in ESR1 up- and downregulated LTEDs cells by Assay for Transposase Accessible Chromatin (ATAC-seq) combined with RNA-seq. Because transcription factors play a primary key role in gene expression and in cell-fate changes including response to therapies, we performed bioinformatic analysis to uncover differential open chromatin sites, select enhancer signatures and identify digital footprints of transcription factors. Unlike ChIP-seq, which identifies binding sites for a known transcription factor, ATAC-seq permits analysis of all sites potentially accessible to the transcriptional machinery and predict occupancy by ~700 different transcription factors at once. With these algorithms, we found that Grainyhead-like 2 (GRHL2) is associated with development of both ESR1 overexpressed and downregulated types of LTEDs. Furthermore, RNA-seq analysis showed that GRHL2 target genes are upregulated in LTEDs. Furthermore, in the METABRIC database, high GRHL2 expression is significantly associated with poor outcome in ER+ breast cancer patients treated with endocrine therapy. Taken together, these data suggest that GRHL2 may be a novel important transcription factor in endocrine therapy resistant breast cancer regardless of fulvestrant sensitivity. Our studies provide evidence that chromatin landscape analysis, coupled with the transcription factor network algorithm, is capable of identifying novel biomarkers or therapeutic targets, This approach will expand the range of translational research as it is applicable to many cancers and diseases. Citation Format: Saori Fujiwara, Songjoon Baek, Kaustubh Wagh, Diana Stavreva, Lyuba Varticovski, Gordon Hager. Chromatin landscape analysis based on ATAC-seq and RNA- seq reveals that GRHL2 is a novel key transcription factor for endocrine therapy resistance [abstract]. In: Proceedings of the 2022 San Antonio Breast Cancer Symposium; 2022 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2023;83(5 Suppl):Abstract nr P3-11-06.
Read full abstract