GRHL2 motif is associated with intratumor heterogeneity of cis-regulatory elements in luminal breast cancer

In breast cancer patients, tumor heterogeneity is associated with prognosis and therapeutic response; however, the epigenetic diversity that exists in primary tumors remains unknown. Using a single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq), we obtained the chromatin accessibility profiles of 12,452 cells from 16 breast cancer patients including 11 luminal, 1 luminal-HER2, 1 HER2(+), and 3 triple-negative subtypes. Via this profiling process, tumors were classified into cancer cells and the tumor microenvironment, highlighting the heterogeneity of disease-related pathways including estrogen receptor (ER) signaling. Furthermore, the coexistence of cancer cell clusters with different ER binding motif enrichments was identified in a single ER(+) tumor. In a cluster with reduced ER motif enrichment, we identified GRHL2, a transcription factor, as the most enriched motif, and it cooperated with FOXA1 to initiate endocrine resistance. Coaccessibility analysis revealed that GRHL2 binding elements potentially regulate genes associated with endocrine resistance, metastasis, and poor prognosis in patients that received hormonal therapy. Overall, our study suggests that epigenetic heterogeneity could lead to endocrine resistance and poor prognosis in breast cancer patients and it offers a large-scale resource for further cancer research.