Differential chromatin accessibility and transcriptional dynamics define breast cancer subtypes and their lineages.

Breast 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.

Epigenetic regulation during cancer transitions across 11 tumour types.

Chromatin accessibility is essential in regulating gene expression and cellular identity, and alterations in accessibility have been implicated in driving cancer initiation, progression and metastasis1-4. Although the genetic contributions to oncogenic transitions have been investigated, epigenetic drivers remain less understood. Here we constructed a pan-cancer epigenetic and transcriptomic atlas using single-nucleus chromatin accessibility data (using single-nucleus assay for transposase-accessible chromatin) from 225 samples and matched single-cell or single-nucleus RNA-sequencing expression data from 206 samples. With over 1 million cells from each platform analysed through the enrichment of accessible chromatin regions, transcription factor motifs and regulons, we identified epigenetic drivers associated with cancer transitions. Some epigenetic drivers appeared in multiple cancers (for example, regulatory regions of ABCC1 and VEGFA; GATA6 and FOX-family motifs), whereas others were cancer specific (for example, regulatory regions of FGF19, ASAP2 and EN1, and the PBX3 motif). Among epigenetically altered pathways, TP53, hypoxia and TNF signalling were linked to cancer initiation, whereas oestrogen response, epithelial-mesenchymal transition and apical junction were tied to metastatic transition. Furthermore, we revealed a marked correlation between enhancer accessibility and gene expression and uncovered cooperation between epigenetic and genetic drivers. This atlas provides a foundation for further investigation of epigenetic dynamics in cancer transitions.

A phase II study of palbociclib plus letrozole plus trastuzumab as neoadjuvant treatment for clinical stages II and III ER+ HER2+ breast cancer (PALTAN).

Patients with ER+/HER2+ breast cancer (BC) are less likely to achieve pathological complete response (pCR) after chemotherapy with dual HER2 blockade than ER-/HER2+ BC. Endocrine therapy plus trastuzumab is effective in advanced ER+/HER2+ BC. Inhibition of CDK4/6 and HER2 results in synergistic cell proliferation reduction. We combined palbociclib, letrozole, and trastuzumab (PLT) as a chemotherapy-sparing regimen. We evaluated neoadjuvant PLT in early ER+/HER2+ BC. Primary endpoint was pCR after 16 weeks. Research biopsies were performed for whole exome and RNA sequencing, PAM50 subtyping, and Ki67 assessment for complete cell cycle arrest (CCCA: Ki67 ≤ 2.7%). After 26 patients, accrual stopped due to futility. pCR (residual cancer burden-RCB 0) was 7.7%, RCB 0/I was 38.5%. Grade (G) 3/4 treatment-emergent adverse events occurred in 19. Among these, G3/4 neutropenia was 50%, hypertension 26.9%, and leucopenia 7.7%. Analysis indicated CCCA in 85% at C1 day 15 (C1D15), compared to 27% at surgery after palbociclib was discontinued. Baseline PAM50 subtyping identified 31.2% HER2-E, 43.8% Luminal B, and 25% Luminal A. 161 genes were differentially expressed comparing C1D15 to baseline. MKI67, TK1, CCNB1, AURKB, and PLK1 were among the genes downregulated, consistent with CCCA at C1D15. Molecular Signatures Database gene-sets analyses demonstrated downregulated processes involved in proliferation, ER and mTORC1 signaling, and DNA damage repair at C1D15, consistent with the study drug's mechanisms of action. Neoadjuvant PLT showed a pCR of 7.7% and an RCB 0/I rate of 38.5%. RNA sequencing and Ki67 data indicated potent anti-proliferative effects of study treatments. ClinicalTrials.gov- NCT02907918.

Spatially restricted drivers and transitional cell populations cooperate with the microenvironment in untreated and chemo-resistant pancreatic cancer.

Pancreatic ductal adenocarcinoma is a lethal disease with limited treatment options and poor survival. We studied 83 spatial samples from 31 patients (11 treatment-naïve and 20 treated) using single-cell/nucleus RNA sequencing, bulk-proteogenomics, spatial transcriptomics and cellular imaging. Subpopulations of tumor cells exhibited signatures of proliferation, KRAS signaling, cell stress and epithelial-to-mesenchymal transition. Mapping mutations and copy number events distinguished tumor populations from normal and transitional cells, including acinar-to-ductal metaplasia and pancreatic intraepithelial neoplasia. Pathology-assisted deconvolution of spatial transcriptomic data identified tumor and transitional subpopulations with distinct histological features. We showed coordinated expression of TIGIT in exhausted and regulatory T cells and Nectin in tumor cells. Chemo-resistant samples contain a threefold enrichment of inflammatory cancer-associated fibroblasts that upregulate metallothioneins. Our study reveals a deeper understanding of the intricate substructure of pancreatic ductal adenocarcinoma tumors that could help improve therapy for patients with this disease.

Treg depletion potentiates checkpoint inhibition in claudin-low breast cancer.

Claudin-low breast cancer is an aggressive subtype that confers poor prognosis and is found largely within the clinical triple-negative group of breast cancer patients. Here, we have shown that intrinsic and immune cell gene signatures distinguish the claudin-low subtype clinically as well as in mouse models of other breast cancer subtypes. Despite adaptive immune cell infiltration in claudin-low tumors, treatment with immune checkpoint inhibitory antibodies against cytotoxic T lymphocyte-associated protein 4 (CTLA-4) and programmed death receptor 1 (PD-1) were ineffective in controlling tumor growth. CD4+FoxP3+ Tregs represented a large proportion of the tumor-infiltrating lymphocytes (TILs) in claudin-low tumors, and Tregs isolated from tumor-bearing mice were able to suppress effector T cell responses. Tregs in the tumor microenvironment highly expressed PD-1 and were recruited partly through tumor generation of the chemokine CXCL12. Antitumor efficacy required stringent Treg depletion combined with checkpoint inhibition; delays in tumor growth were not observed using therapies that modestly diminished the number of Tregs in the tumor microenvironment. This study provides evidence that the recruitment of Tregs to the tumor microenvironment inhibits an effective antitumor immune response and highlights early Treg recruitment as a possible mechanism for the lack of response to immune checkpoint blockade antibodies in specific subtypes of cancer that are heavily infiltrated with adaptive immune cells.

Claudin-low bladder tumors are immune infiltrated and actively immune suppressed.

We report the discovery of a claudin-low molecular subtype of high-grade bladder cancer that shares characteristics with the homonymous subtype of breast cancer. Claudin-low bladder tumors were enriched for multiple genetic features including increased rates of RB1, EP300, and NCOR1 mutations; increased frequency of EGFR amplification; decreased rates of FGFR3, ELF3, and KDM6A mutations; and decreased frequency of PPARG amplification. While claudin-low tumors showed the highest expression of immune gene signatures, they also demonstrated gene expression patterns consistent with those observed in active immunosuppression. This did not appear to be due to differences in predicted neoantigen burden, but rather was associated with broad upregulation of cytokine and chemokine levels from low PPARG activity, allowing unopposed NFKB activity. Taken together, these results define a molecular subtype of bladder cancer with distinct molecular features and an immunologic profile that would, in theory, be primed for immunotherapeutic response.

Assembly-based inference of B-cell receptor repertoires from short read RNA sequencing data with V'DJer.

MOTIVATION: B-cell receptor (BCR) repertoire profiling is an important tool for understanding the biology of diverse immunologic processes. Current methods for analyzing adaptive immune receptor repertoires depend upon PCR amplification of VDJ rearrangements followed by long read amplicon sequencing spanning the VDJ junctions. While this approach has proven to be effective, it is frequently not feasible due to cost or limited sample material. Additionally, there are many existing datasets where short-read RNA sequencing data are available but PCR amplified BCR data are not. RESULTS: We present here V'DJer, an assembly-based method that reconstructs adaptive immune receptor repertoires from short-read RNA sequencing data. This method captures expressed BCR loci from a standard RNA-seq assay. We applied this method to 473 Melanoma samples from The Cancer Genome Atlas and demonstrate V'DJer's ability to accurately reconstruct BCR repertoires from short read mRNA-seq data. AVAILABILITY AND IMPLEMENTATION: V'DJer is implemented in C/C ++, freely available for academic use and can be downloaded from Github: https://github.com/mozack/vdjer CONTACT: benjamin_vincent@med.unc.edu or parkerjs@email.unc.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Genomic Analysis of Immune Cell Infiltrates Across 11 Tumor Types.

BACKGROUND: Immune infiltration of the tumor microenvironment has been associated with improved survival for some patients with solid tumors. The precise makeup and prognostic relevance of immune infiltrates across a broad spectrum of tumors remain unclear. METHODS: Using mRNA sequencing data from The Cancer Genome Atlas (TCGA) from 11 tumor types representing 3485 tumors, we evaluated lymphocyte and macrophage gene expression by tissue type and by genomic subtypes defined within and across tumor tissue of origin (Cox proportional hazards, Pearson correlation). We investigated clonal diversity of B-cell infiltrates through calculating B-cell receptor (BCR) repertoire sequence diversity. All statistical tests were two-sided. RESULTS: High expression of T-cell and B-cell signatures predicted improved overall survival across many tumor types including breast, lung, and melanoma (breast CD8_T_Cells hazard ratio [HR] = 0.36, 95% confidence interval [CI] = 0.16 to 0.81, P = .01; lung adenocarcinoma B_Cell_60gene HR = 0.71, 95% CI = 0.58 to 0.87, P = 7.80E-04; melanoma LCK HR = 0.86, 95% CI = 0.79 to 0.94, P = 6.75E-04). Macrophage signatures predicted worse survival in GBM, as did B-cell signatures in renal tumors (Glioblastoma Multiforme [GBM]: macrophages HR = 1.62, 95% CI = 1.17 to 2.26, P = .004; renal: B_Cell_60gene HR = 1.17, 95% CI = 1.04 to 1.32, P = .009). BCR diversity was associated with survival beyond gene segment expression in melanoma (HR = 2.67, 95% CI = 1.32 to 5.40, P = .02) and renal cell carcinoma (HR = 0.36, 95% CI = 0.15 to 0.87, P = .006). CONCLUSIONS: These data support existing studies suggesting that in diverse tissue types, heterogeneous immune infiltrates are present and typically portend an improved prognosis. In some tumor types, BCR diversity was also associated with survival. Quantitative genomic signatures of immune cells warrant further testing as prognostic markers and potential biomarkers of response to cancer immunotherapy.

Prognostic B-cell signatures using mRNA-seq in patients with subtype-specific breast and ovarian cancer.

PURPOSE: Lymphocytic infiltration of tumors predicts improved survival in patients with breast cancer. Previous studies have suggested that this survival benefit is confined predominantly to the basal-like subtype. Immune infiltration in ovarian tumors is also associated with improved prognosis. Currently, it is unclear what aspects of the immune response mediate this improved outcome. EXPERIMENTAL DESIGN: Using The Cancer Genome Atlas mRNA-seq data and a large microarray dataset, we evaluated adaptive immune gene expression by genomic subtype in breast and ovarian cancer. To investigate B-cells observed to be prognostic within specific subtypes, we developed methods to analyze B-cell population diversity and degree of somatic hypermutation (SHM) from B-cell receptor (BCR) sequences in mRNA-seq data. RESULTS: Improved metastasis-free/progression-free survival was correlated with B-cell gene expression signatures, which were restricted mainly to the basal-like and HER2-enriched breast cancer subtypes and the immunoreactive ovarian cancer subtype. Consistent with a restricted epitope-driven response, a subset of basal-like and HER2-enriched breast tumors and immunoreactive ovarian tumors showed high expression of a low-diversity population of BCR gene segments. More BCR segments showed improved prognosis with increased expression in basal-like breast tumors and immunoreactive ovarian tumors compared with other subtypes. Basal-like and HER2-enriched tumors exhibited more BCR sequence variants in regions consistent with SHM. CONCLUSION: Taken together, these data suggest the presence of a productive and potentially restricted antitumor B-cell response in basal-like breast and immunoreactive ovarian cancers. Immunomodulatory therapies that support B-cell responses may be a promising therapeutic approach to targeting these B-cell infiltrated tumors.

Tumor-specific retargeting of an oncogenic transcription factor chimera results in dysregulation of chromatin and transcription.

Chromosomal translocations involving transcription factor genes have been identified in an increasingly wide range of cancers. Some translocations can create a protein "chimera" that is composed of parts from different proteins. How such chimeras cause cancer, and why they cause cancer in some cell types but not others, is not understood. One such chimera is EWS-FLI, the most frequently occurring translocation in Ewing Sarcoma, a malignant bone and soft tissue tumor of children and young adults. Using EWS-FLI and its parental transcription factor, FLI1, we created a unique experimental system to address questions regarding the genomic mechanisms by which chimeric transcription factors cause cancer. We found that in tumor cells, EWS-FLI targets regions of the genome distinct from FLI1, despite identical DNA-binding domains. In primary endothelial cells, however, EWS-FLI and FLI1 demonstrate similar targeting. To understand this mistargeting, we examined chromatin organization. Regions targeted by EWS-FLI are normally repressed and nucleosomal in primary endothelial cells. In tumor cells, however, bound regions are nucleosome depleted and harbor the chromatin signature of enhancers. We next demonstrated that through chimerism, EWS-FLI acquired the ability to alter chromatin. Expression of EWS-FLI results in nucleosome depletion at targeted sites, whereas silencing of EWS-FLI in tumor cells restored nucleosome occupancy. Thus, the EWS-FLI chimera acquired chromatin-altering activity, leading to mistargeting, chromatin disruption, and ultimately, transcriptional dysregulation.