Differential chromatin accessibility and transcriptional dynamics define breast cancer subtypes and their lineages
Authors: Iglesia Michael D. , Jayasinghe Reyka G. , Chen Siqi , Terekhanova Nadezhda V. , Herndon John M. , Storrs Erik , Karpova Alla , Zhou Daniel Cui , Al Deen Nataly Naser , Shinkle Andrew T. , Lu Rita Jui-Hsien , Caravan Wagma , Houston Andrew , Zhao Yanyan , Sato Kazuhito , Lal Preet , Street Cherease , Rodrigues Fernanda Martins , Southard-Smith Austin N. , Targino da Costa André Luiz N. , Zhu Houxiang , Mo Chia-Kuei , Crowson Lisa , Fulton Robert S. , Wyczalkowski Matthew A. , Fronick Catrina C. , Fulton Lucinda A. , Sun Hua , Davies Sherri R. , Appelbaum Elizabeth L. , Chasnoff Sara E. , Carmody Madelyn , Brooks Candace , Liu Ruiyang , Wendl Michael C. , Oh Clara , Bender Diane , Cruchaga Carlos , Harari Oscar , Bredemeyer Andrea , Lavine Kory , Bose Ron , Margenthaler Julie , Held Jason M. , Achilefu Samuel , Ademuyiwa Foluso , Aft Rebecca , Ma Cynthia , Colditz Graham A. , Ju Tao , Oh Stephen T. , Fitzpatrick James , Hwang E. Shelley , Shoghi Kooresh I. , Chheda Milan G. , Veis Deborah J. , Chen Feng , Fields Ryan C. , Gillanders William E. , Ding Li 
Published: November 2023

ABSTRACTBreast cancer is a heterogeneous disease, and treatment is guided by biomarker profiles representing distinct molecular subtypes. Breast cancer arises from the breast ductal epithelium, and experimental data suggests breast cancer subtypes have different cells of origin within that lineage. The precise cells of origin for each subtype and the transcriptional networks that characterize these tumor-normal lineages are not established. In this work, we applied bulk, single-cell (sc), and single-nucleus (sn) multi-omic techniques as well as spatial transcriptomics and multiplex imaging on 61 samples from 37 breast cancer patients to show characteristic links in gene expression and chromatin accessibility between breast cancer subtypes and their putative cells of origin. We applied the PAM50 subtyping algorithm in tandem with bulk RNA-seq and snRNA-seq to reliably subtype even low-purity tumor samples and confirm promoter accessibility using snATAC. Trajectory analysis of chromatin accessibility and differentially accessible motifs clearly connected progenitor populations with breast cancer subtypes supporting the cell of origin for basal-like and luminal A and B tumors. Regulatory network analysis of transcription factors underscored the importance of BHLHE40 in luminal breast cancer and luminal mature cells, and KLF5 in basal-like tumors and luminal progenitor cells. Furthermore, we identify key genes defining the basal-like (PRKCA,SOX6,RGS6,KCNQ3) and luminal A/B (FAM155A,LRP1B) lineages, with expression in both precursor and cancer cells and further upregulation in tumors. Exhausted CTLA4-expressing CD8+ T cells were enriched in basal-like breast cancer, suggesting altered means of immune dysfunction among breast cancer subtypes. We used spatial transcriptomics and multiplex imaging to provide spatial detail for key markers of benign and malignant cell types and immune cell colocation. These findings demonstrate analysis of paired transcription and chromatin accessibility at the single cell level is a powerful tool for investigating breast cancer lineage development and highlight transcriptional networks that define basal and luminal breast cancer lineages.