Heterogeneity and transcriptional drivers of triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is a heterogeneous disease with limited treatment options. To characterize TNBC heterogeneity, we defined transcriptional, epigenetic, and metabolic subtypes and subtype-driving super-enhancers and transcription factors by combining functional and molecular profiling with computational analyses. Single-cell RNA sequencing revealed relative homogeneity of the major transcriptional subtypes (luminal, basal, and mesenchymal) within samples. We found that mesenchymal TNBCs share features with mesenchymal neuroblastoma and rhabdoid tumors and that the PRRX1 transcription factor is a key driver of these tumors. PRRX1 is sufficient for inducing mesenchymal features in basal but not in luminal TNBC cells via reprogramming super-enhancer landscapes, but it is not required for mesenchymal state maintenance or for cellular viability. Our comprehensive, large-scale, multiplatform, multiomics study of both experimental and clinical TNBC is an important resource for the scientific and clinical research communities and opens venues for future investigation.

JAK-STAT Signaling in Inflammatory Breast Cancer Enables Chemotherapy-Resistant Cell States.

Inflammatory breast cancer (IBC) is a difficult-to-treat disease with poor clinical outcomes due to high risk of metastasis and resistance to treatment. In breast cancer, CD44+CD24- cells possess stem cell-like features and contribute to disease progression, and we previously described a CD44+CD24-pSTAT3+ breast cancer cell subpopulation that is dependent on JAK2/STAT3 signaling. Here we report that CD44+CD24- cells are the most frequent cell type in IBC and are commonly pSTAT3+. Combination of JAK2/STAT3 inhibition with paclitaxel decreased IBC xenograft growth more than either agent alone. IBC cell lines resistant to paclitaxel and doxorubicin were developed and characterized to mimic therapeutic resistance in patients. Multi-omic profiling of parental and resistant cells revealed enrichment of genes associated with lineage identity and inflammation in chemotherapy-resistant derivatives. Integrated pSTAT3 chromatin immunoprecipitation sequencing and RNA sequencing (RNA-seq) analyses showed pSTAT3 regulates genes related to inflammation and epithelial-to-mesenchymal transition (EMT) in resistant cells, as well as PDE4A, a cAMP-specific phosphodiesterase. Metabolomic characterization identified elevated cAMP signaling and CREB as a candidate therapeutic target in IBC. Investigation of cellular dynamics and heterogeneity at the single cell level during chemotherapy and acquired resistance by CyTOF and single cell RNA-seq identified mechanisms of resistance including a shift from luminal to basal/mesenchymal cell states through selection for rare preexisting subpopulations or an acquired change. Finally, combination treatment with paclitaxel and JAK2/STAT3 inhibition prevented the emergence of the mesenchymal chemo-resistant subpopulation. These results provide mechanistic rational for combination of chemotherapy with inhibition of JAK2/STAT3 signaling as a more effective therapeutic strategy in IBC. SIGNIFICANCE: Chemotherapy resistance in inflammatory breast cancer is driven by the JAK2/STAT3 pathway, in part via cAMP/PKA signaling and a cell state switch, which can be overcome using paclitaxel combined with JAK2 inhibitors.

Smooth muscle mineralocorticoid receptor as an epigenetic regulator of vascular ageing.

AIMS: Vascular stiffness increases with age and independently predicts cardiovascular disease risk. Epigenetic changes, including histone modifications, accumulate with age but the global pattern has not been elucidated nor are the regulators known. Smooth muscle cell-mineralocorticoid receptor (SMC-MR) contributes to vascular stiffness in ageing mice. Thus, we investigated the regulatory role of SMC-MR in vascular epigenetics and stiffness. METHODS AND RESULTS: Mass spectrometry-based proteomic profiling of all histone modifications completely distinguished 3 from 12-month-old mouse aortas. Histone-H3 lysine-27 (H3K27) methylation (me) significantly decreased in ageing vessels and this was attenuated in SMC-MR-KO littermates. Immunoblotting revealed less H3K27-specific methyltransferase EZH2 with age in MR-intact but not SMC-MR-KO vessels. These ageing changes were examined in primary human aortic (HA)SMC from adult vs. aged donors. MR, H3K27 acetylation (ac), and stiffness gene (connective tissue growth factor, integrin-α5) expression significantly increased, while H3K27me and EZH2 decreased, with age. MR inhibition reversed these ageing changes in HASMC and the decline in stiffness genes was prevented by EZH2 blockade. Atomic force microscopy revealed that MR antagonism decreased intrinsic stiffness and the probability of fibronectin adhesion of aged HASMC. Conversely, ageing induction in young HASMC with H2O2; increased MR, decreased EZH2, enriched H3K27ac and MR at stiffness gene promoters by chromatin immunoprecipitation, and increased stiffness gene expression. In 12-month-old mice, MR antagonism increased aortic EZH2 and H3K27 methylation, increased EZH2 recruitment and decreased H3K27ac at stiffness genes promoters, and prevented ageing-induced vascular stiffness and fibrosis. Finally, in human aortic tissue, age positively correlated with MR and stiffness gene expression and negatively correlated with H3K27me3 while MR and EZH2 are negatively correlated. CONCLUSION: These data support a novel vascular ageing model with rising MR in human SMC suppressing EZH2 expression thereby decreasing H3K27me, promoting MR recruitment and H3K27ac at stiffness gene promoters to induce vascular stiffness and suggests new targets for ameliorating ageing-associated vascular disease.

Down-syndrome-induced senescence disrupts the nuclear architecture of neural progenitors.

Down syndrome (DS) is a genetic disorder driven by the triplication of chromosome 21 (T21) and characterized by a wide range of neurodevelopmental and physical disabilities. Transcriptomic analysis of tissue samples from individuals with DS has revealed that T21 induces a genome-wide transcriptional disruption. However, the consequences of T21 on the nuclear architecture and its interplay with the transcriptome remain unknown. In this study, we find that unlike human induced pluripotent stem cells (iPSCs), iPSC-derived neural progenitor cells (NPCs) exhibit genome-wide "chromosomal introversion," disruption of lamina-associated domains, and global chromatin accessibility changes in response to T21, consistent with the transcriptional and nuclear architecture changes characteristic of senescent cells. Treatment of T21-harboring NPCs with senolytic drugs alleviates the transcriptional, molecular, and cellular dysfunctions associated with DS. Our findings provide a mechanistic link between T21 and global transcriptional disruption and indicate that senescence-associated phenotypes may play a key role in the neurodevelopmental pathogenesis of DS.

Proteomic profiling dataset of chemical perturbations in multiple biological backgrounds.

While gene expression profiling has traditionally been the method of choice for large-scale perturbational profiling studies, proteomics has emerged as an effective tool in this context for directly monitoring cellular responses to perturbations. We previously reported a pilot library containing 3400 profiles of multiple perturbations across diverse cellular backgrounds in the reduced-representation phosphoproteome (P100) and chromatin space (Global Chromatin Profiling, GCP). Here, we expand our original dataset to include profiles from a new set of cardiotoxic compounds and from astrocytes, an additional neural cell model, totaling 5300 proteomic signatures. We describe filtering criteria and quality control metrics used to assess and validate the technical quality and reproducibility of our data. To demonstrate the power of the library, we present two case studies where data is queried using the concept of "connectivity" to obtain biological insight. All data presented in this study have been deposited to the ProteomeXchange Consortium with identifiers PXD017458 (P100) and PXD017459 (GCP) and can be queried at https://clue.io/proteomics .

Causal interactions from proteomic profiles: Molecular data meet pathway knowledge.

We present a computational method to infer causal mechanisms in cell biology by analyzing changes in high-throughput proteomic profiles on the background of prior knowledge captured in biochemical reaction knowledge bases. The method mimics a biologist's traditional approach of explaining changes in data using prior knowledge but does this at the scale of hundreds of thousands of reactions. This is a specific example of how to automate scientific reasoning processes and illustrates the power of mapping from experimental data to prior knowledge via logic programming. The identified mechanisms can explain how experimental and physiological perturbations, propagating in a network of reactions, affect cellular responses and their phenotypic consequences. Causal pathway analysis is a powerful and flexible discovery tool for a wide range of cellular profiling data types and biological questions. The automated causation inference tool, as well as the source code, are freely available at http://causalpath.org.

Epigenetic silencing by SETDB1 suppresses tumour intrinsic immunogenicity.

Epigenetic dysregulation is a defining feature of tumorigenesis that is implicated in immune escape1,2. Here, to identify factors that modulate the immune sensitivity of cancer cells, we performed in vivo CRISPR-Cas9 screens targeting 936 chromatin regulators in mouse tumour models treated with immune checkpoint blockade. We identified the H3K9 methyltransferase SETDB1 and other members of the HUSH and KAP1 complexes as mediators of immune escape3-5. We also found that amplification of SETDB1 (1q21.3) in human tumours is associated with immune exclusion and resistance to immune checkpoint blockade. SETDB1 represses broad domains, primarily within the open genome compartment. These domains are enriched for transposable elements (TEs) and immune clusters associated with segmental duplication events, a central mechanism of genome evolution6. SETDB1 loss derepresses latent TE-derived regulatory elements, immunostimulatory genes, and TE-encoded retroviral antigens in these regions, and triggers TE-specific cytotoxic T cell responses in vivo. Our study establishes SETDB1 as an epigenetic checkpoint that suppresses tumour-intrinsic immunogenicity, and thus represents a candidate target for immunotherapy.

Noncanonical open reading frames encode functional proteins essential for cancer cell survival.

Although genomic analyses predict many noncanonical open reading frames (ORFs) in the human genome, it is unclear whether they encode biologically active proteins. Here we experimentally interrogated 553 candidates selected from noncanonical ORF datasets. Of these, 57 induced viability defects when knocked out in human cancer cell lines. Following ectopic expression, 257 showed evidence of protein expression and 401 induced gene expression changes. Clustered regularly interspaced short palindromic repeat (CRISPR) tiling and start codon mutagenesis indicated that their biological effects required translation as opposed to RNA-mediated effects. We found that one of these ORFs, G029442-renamed glycine-rich extracellular protein-1 (GREP1)-encodes a secreted protein highly expressed in breast cancer, and its knockout in 263 cancer cell lines showed preferential essentiality in breast cancer-derived lines. The secretome of GREP1-expressing cells has an increased abundance of the oncogenic cytokine GDF15, and GDF15 supplementation mitigated the growth-inhibitory effect of GREP1 knockout. Our experiments suggest that noncanonical ORFs can express biologically active proteins that are potential therapeutic targets.

Avant-garde: an automated data-driven DIA data curation tool.

Several challenges remain in data-independent acquisition (DIA) data analysis, such as to confidently identify peptides, define integration boundaries, remove interferences, and control false discovery rates. In practice, a visual inspection of the signals is still required, which is impractical with large datasets. We present Avant-garde as a tool to refine DIA (and parallel reaction monitoring) data. Avant-garde uses a novel data-driven scoring strategy: signals are refined by learning from the dataset itself, using all measurements in all samples to achieve the best optimization. We evaluate the performance of Avant-garde using benchmark DIA datasets and show that it can determine the quantitative suitability of a peptide peak, and reach the same levels of selectivity, accuracy, and reproducibility as manual validation. Avant-garde is complementary to existing DIA analysis engines and aims to establish a strong foundation for subsequent analysis of quantitative mass spectrometry data.

The nucleosome acidic patch and H2A ubiquitination underlie mSWI/SNF recruitment in synovial sarcoma.

Interactions between chromatin-associated proteins and the histone landscape play major roles in dictating genome topology and gene expression. Cancer-specific fusion oncoproteins, which display unique chromatin localization patterns, often lack classical DNA-binding domains, presenting challenges in identifying mechanisms governing their site-specific chromatin targeting and function. Here we identify a minimal region of the human SS18-SSX fusion oncoprotein (the hallmark driver of synovial sarcoma) that mediates a direct interaction between the mSWI/SNF complex and the nucleosome acidic patch. This binding results in altered mSWI/SNF composition and nucleosome engagement, driving cancer-specific mSWI/SNF complex targeting and gene expression. Furthermore, the C-terminal region of SSX confers preferential affinity to repressed, H2AK119Ub-marked nucleosomes, underlying the selective targeting to polycomb-marked genomic regions and synovial sarcoma-specific dependency on PRC1 function. Together, our results describe a functional interplay between a key nucleosome binding hub and a histone modification that underlies the disease-specific recruitment of a major chromatin remodeling complex.

Synthetic Lethal and Resistance Interactions with BET Bromodomain Inhibitors in Triple-Negative Breast Cancer.

BET bromodomain inhibitors (BBDIs) are candidate therapeutic agents for triple-negative breast cancer (TNBC) and other cancer types, but inherent and acquired resistance to BBDIs limits their potential clinical use. Using CRISPR and small-molecule inhibitor screens combined with comprehensive molecular profiling of BBDI response and resistance, we identified synthetic lethal interactions with BBDIs and genes that, when deleted, confer resistance. We observed synergy with regulators of cell cycle progression, YAP, AXL, and SRC signaling, and chemotherapeutic agents. We also uncovered functional similarities and differences among BRD2, BRD4, and BRD7. Although deletion of BRD2 enhances sensitivity to BBDIs, BRD7 loss leads to gain of TEAD-YAP chromatin binding and luminal features associated with BBDI resistance. Single-cell RNA-seq, ATAC-seq, and cellular barcoding analysis of BBDI responses in sensitive and resistant cell lines highlight significant heterogeneity among samples and demonstrate that BBDI resistance can be pre-existing or acquired.

Dual protease type XIII/pepsin digestion offers superior resolution and overlap for the analysis of histone tails by HX-MS.

The N-terminal regions of histone proteins (tails) are dynamic elements that protrude from the nucleosome and are involved in many aspects of chromatin organization. Their epigenetic role is well-established, and post-translational modifications (PTMs) present on these regions contribute to transcriptional regulation. While hydrogen/deuterium exchange mass spectrometry (HX-MS) is well-suited for the analysis of dynamic structures, it has seldom been employed to analyze histones due to the poor N-terminal coverage obtained using pepsin. Here, we test the applicability of a dual protease type XIII/pepsin digestion column to this class of proteins. We optimize online digestion conditions using the H4 monomer, and extend the method to the analysis of histones in monomeric states and nucleosome core particles (NCPs). We show that the dual protease column generates many short and overlapping N-terminal peptides. We evaluate our method by performing hydrogen exchange experiments of NCPs for different time points and present full coverage of the tails at excellent resolution. We further employ electron transfer dissociation and showcase an unprecedented degree of overlap across multiple peptides that is several fold higher than previously reported methods. The method we report here may be readily applied to the HX-MS investigation of histone dynamics and to the footprints of histone binding proteins on nucleosomes.

Phosphorylation of RAB7 by TBK1/IKKε Regulates Innate Immune Signaling in Triple-Negative Breast Cancer.

Triple-negative breast cancer (TNBC) is a heterogeneous disease enriched for mutations in PTEN and dysregulation of innate immune signaling. Here, we demonstrate that Rab7, a recently identified substrate of PTEN phosphatase activity, is also a substrate of the innate immune signaling kinases TANK-binding kinase 1 (TBK1)/IκB kinase ε (IKKε) on the same serine-72 (S72) site. An unbiased search for novel TBK1/IKKε substrates using stable isotope labeling with amino acids in cell culture phosphoproteomic analysis identified Rab7-S72 as a top hit. PTEN-null TNBC cells expressing a phosphomimetic version of Rab7-S72 exhibited diffuse cytosolic Rab7 localization and enhanced innate immune signaling, in contrast to a kinase-resistant version, which localized to active puncta that promote lysosomal-mediated stimulator of interferon genes (STING) degradation. Thus, convergence of PTEN loss and TBK1/IKKε activation on Rab7-S72 phosphorylation limited STING turnover and increased downstream production of IRF3 targets including CXCL10, CCL5, and IFNß. Consistent with this data, PTEN-null TNBC tumors expressed higher levels of STING, and PTEN-null TNBC cell lines were hyperresponsive to STING agonists. Together, these findings begin to uncover how innate immune signaling is dysregulated downstream of TBK1/IKKε in a subset of TNBCs and reveals previously unrecognized cross-talk with STING recycling that may have implications for STING agonism in the clinic. SIGNIFICANCE: These findings identify Rab7 as a substrate for TBK1 for regulation of innate immune signaling, thereby providing important insight for strategies aimed at manipulating the immune response to enhance therapeutic efficacy in TNBC.

High-fat diet fuels prostate cancer progression by rewiring the metabolome and amplifying the MYC program.

Systemic metabolic alterations associated with increased consumption of saturated fat and obesity are linked with increased risk of prostate cancer progression and mortality, but the molecular underpinnings of this association are poorly understood. Here, we demonstrate in a murine prostate cancer model, that high-fat diet (HFD) enhances the MYC transcriptional program through metabolic alterations that favour histone H4K20 hypomethylation at the promoter regions of MYC regulated genes, leading to increased cellular proliferation and tumour burden. Saturated fat intake (SFI) is also associated with an enhanced MYC transcriptional signature in prostate cancer patients. The SFI-induced MYC signature independently predicts prostate cancer progression and death. Finally, switching from a high-fat to a low-fat diet, attenuates the MYC transcriptional program in mice. Our findings suggest that in primary prostate cancer, dietary SFI contributes to tumour progression by mimicking MYC over expression, setting the stage for therapeutic approaches involving changes to the diet.

Chasing Tails: Cathepsin-L Improves Structural Analysis of Histones by HX-MS.

The N-terminal regions (tails) of histone proteins are dynamic elements that protrude from the nucleosome and are involved in many aspects of chromatin organization. Their epigenetic role is well-established, and post-translational modifications present on these regions contribute to transcriptional regulation. Considering their biological significance, relatively few structural details have been established for histone tails, mainly because of their inherently disordered nature. Although hydrogen/deuterium exchange mass spectrometry (HX-MS) is well-suited for the analysis of dynamic structures, it has seldom been employed in this context, presumably because of the poor N-terminal coverage provided by pepsin. Inspired from histone-clipping events, we profiled the activity of cathepsin-L under HX-MS quench conditions and characterized its specificity employing the four core histones (H2A, H2B, H3 and H4). Cathepsin-L demonstrated cleavage patterns that were substrate- and pH-dependent. Cathepsin-L generated overlapping N-terminal peptides about 20 amino acids long for H2A, H3, and H4 proving its suitability for the analysis of histone tails dynamics. We developed a comprehensive HX-MS method in combination with pepsin and obtained full sequence coverage for all histones. We employed our method to analyze histones H3 and H4. We observe rapid deuterium exchange of the N-terminal tails and cooperative unfolding (EX1 kinetics) in the histone-fold domains of histone monomers in-solution. Overall, this novel strategy opens new avenues for investigating the dynamic properties of histones that are not apparent from the crystal structures, providing insights into the structural basis of the histone code.

Histone demethylase KDM6A directly senses oxygen to control chromatin and cell fate.

Oxygen sensing is central to metazoan biology and has implications for human disease. Mammalian cells express multiple oxygen-dependent enzymes called 2-oxoglutarate (OG)-dependent dioxygenases (2-OGDDs), but they vary in their oxygen affinities and hence their ability to sense oxygen. The 2-OGDD histone demethylases control histone methylation. Hypoxia increases histone methylation, but whether this reflects direct effects on histone demethylases or indirect effects caused by the hypoxic induction of the HIF (hypoxia-inducible factor) transcription factor or the 2-OG antagonist 2-hydroxyglutarate (2-HG) is unclear. Here, we report that hypoxia promotes histone methylation in a HIF- and 2-HG-independent manner. We found that the H3K27 histone demethylase KDM6A/UTX, but not its paralog KDM6B, is oxygen sensitive. KDM6A loss, like hypoxia, prevented H3K27 demethylation and blocked cellular differentiation. Restoring H3K27 methylation homeostasis in hypoxic cells reversed these effects. Thus, oxygen directly affects chromatin regulators to control cell fate.

A Curated Resource for Phosphosite-specific Signature Analysis.

Signaling pathways are orchestrated by post-translational modifications (PTMs) such as phosphorylation. However, pathway analysis of PTM data sets generated by mass spectrometry (MS)-based proteomics is typically performed at a gene-centric level because of the lack of appropriately curated PTM signature databases and bioinformatic tools that leverage PTM site-specific information. Here we present the first version of PTMsigDB, a database of modification site-specific signatures of perturbations, kinase activities and signaling pathways curated from more than 2,500 publications. We adapted the widely used single sample Gene Set Enrichment Analysis approach to utilize PTMsigDB, enabling PTMSignature Enrichment Analysis (PTM-SEA) of quantitative MS data. We used a well-characterized data set of epidermal growth factor (EGF)-perturbed cancer cells to evaluate our approach and demonstrated better representation of signaling events compared with gene-centric methods. We then applied PTM-SEA to analyze the phosphoproteomes of cancer cells treated with cell-cycle inhibitors and detected mechanism-of-action specific signatures of cell cycle kinases. We also applied our methods to analyze the phosphoproteomes of PI3K-inhibited human breast cancer cells and detected signatures of compounds inhibiting PI3K as well as targets downstream of PI3K (AKT, MAPK/ERK) covering a substantial fraction of the PI3K pathway. PTMsigDB and PTM-SEA can be freely accessed at https://github.com/broadinstitute/ssGSEA2.0.