Simultaneous profiling of multiple chromatin proteins in the same cells.

Methods derived from CUT&RUN and CUT&Tag enable genome-wide mapping of the localization of proteins on chromatin from as few as one cell. These and other mapping approaches focus on one protein at a time, preventing direct measurements of co-localization of different chromatin proteins in the same cells and requiring prioritization of targets where samples are limiting. Here, we describe multi-CUT&Tag, an adaptation of CUT&Tag that overcomes these hurdles by using antibody-specific barcodes to simultaneously map multiple proteins in the same cells. Highly specific multi-CUT&Tag maps of histone marks and RNA Polymerase II uncovered sites of co-localization in the same cells, active and repressed genes, and candidate cis-regulatory elements. Single-cell multi-CUT&Tag profiling facilitated identification of distinct cell types from a mixed population and characterization of cell-type-specific chromatin architecture. In sum, multi-CUT&Tag increases the information content per cell of epigenomic maps, facilitating direct analysis of the interplay of different chromatin proteins.

Functional genomic screens with death rate analyses reveal mechanisms of drug action.

A common approach for understanding how drugs induce their therapeutic effects is to identify the genetic determinants of drug sensitivity. Because 'chemo-genetic profiles' are performed in a pooled format, inference of gene function is subject to several confounding influences related to variation in growth rates between clones. In this study, we developed Method for Evaluating Death Using a Simulation-assisted Approach (MEDUSA), which uses time-resolved measurements, along with model-driven constraints, to reveal the combination of growth and death rates that generated the observed drug response. MEDUSA is uniquely effective at identifying death regulatory genes. We apply MEDUSA to characterize DNA damage-induced lethality in the presence and absence of p53. Loss of p53 switches the mechanism of DNA damage-induced death from apoptosis to a non-apoptotic death that requires high respiration. These findings demonstrate the utility of MEDUSA both for determining the genetic dependencies of lethality and for revealing opportunities to potentiate chemo-efficacy in a cancer-specific manner.

Deletion of Cdkn1b in ACI rats leads to increased proliferation and pregnancy-associated changes in the mammary gland due to perturbed systemic endocrine environment.

Mammary epithelial progenitors are the normal cell-of-origin of breast cancer. We previously defined a population of p27+ quiescent hormone-responsive progenitor cells in the normal human breast whose frequency associates with breast cancer risk. Here, we describe that deletion of the Cdkn1b gene encoding the p27 cyclin-dependent kinase inhibitor in the estrogen-induced mammary tumor-susceptible ACI rat strain leads to a decrease in the relative frequencies of Cd49b+ mammary luminal epithelial progenitors and pregnancy-related differentiation. We show by comprehensive gene expression profiling of purified progenitor and differentiated mammary epithelial cell populations that p27 deletion has the most pronounced effects on luminal progenitors. Cdkn1b-/- females have decreased fertility, but rats that are able to get pregnant had normal litter size and were able to nurse their pups implying that loss of p27 in ACI rats does not completely abrogate ovarian function and lactation. Reciprocal mammary gland transplantation experiments indicate that the p27-loss-induced changes in mammary epithelial cells are not only caused by alterations in their intrinsic properties, but are likely due to altered hormonal signaling triggered by the perturbed systemic endocrine environment observed in Cdkn1b-/- females. We also observed a decrease in the frequency of mammary epithelial cells positive for progesterone receptor (Pr) and FoxA1, known direct transcriptional targets of the estrogen receptor (Erα), and an increase in phospho-Stat5 positive cells commonly induced by prolactin (Prl). Characterization of genome-wide Pr chromatin binding revealed distinct binding patterns in mammary epithelial cells of Cdkn1b+/+ and Cdkn1b-/- females and enrichment in genes with known roles in Notch, ErbB, leptin, and Erα signaling and regulation of G1-S transition. Our data support a role for p27 in regulating the pool size of hormone-responsive luminal progenitors that could impact breast cancer risk.

Calcineurin promotes adaptation to chronic stress through two distinct mechanisms.

Adaptation to environmental stress requires coordination between stress-defense programs and cell cycle progression. The immediate response to many stressors has been well characterized, but how cells survive in challenging environments long-term is unknown. Here, we investigate the role of the stress-activated phosphatase calcineurin (CN) in adaptation to chronic CaCl2 stress in Saccharomyces cerevisiae. We find that prolonged exposure to CaCl2 impairs mitochondrial function and demonstrate that cells respond to this stressor using two CN-dependent mechanisms - one that requires the downstream transcription factor Crz1 and another that is Crz1-independent. Our data indicate that CN maintains cellular fitness by promoting cell cycle progression and preventing CaCl2-induced cell death. When Crz1 is present, transient CN activation suppresses cell death and promotes adaptation despite high levels of mitochondrial loss. However, in the absence of Crz1, prolonged activation of CN prevents mitochondrial loss and further cell death by upregulating glutathione (GSH) biosynthesis genes thereby mitigating damage from reactive oxygen species. These findings illustrate how cells maintain long-term fitness during chronic stress and suggest that CN promotes adaptation in challenging environments by multiple mechanisms.

Calcineurin promotes adaptation to chronic stress through two distinct mechanisms.

Adaptation to environmental stress requires coordination between stress-defense programs and cell cycle progression. The immediate response to many stressors has been well characterized, but how cells survive in challenging environments long term is unknown. Here, we investigate the role of the stress-activated phosphatase calcineurin (CN) in adaptation to chronic CaCl2 stress in Saccharomyces cerevisiae. We find that prolonged exposure to CaCl2 impairs mitochondrial function and demonstrate that cells respond to this stressor using two CN-dependent mechanisms-one that requires the downstream transcription factor Crz1 and another that is Crz1 independent. Our data indicate that CN maintains cellular fitness by promoting cell cycle progression and preventing CaCl2-induced cell death. When Crz1 is present, transient CN activation suppresses cell death and promotes adaptation despite high levels of mitochondrial loss. However, in the absence of Crz1, prolonged activation of CN prevents mitochondrial loss and further cell death by upregulating glutathione biosynthesis genes thereby mitigating damage from reactive oxygen species. These findings illustrate how cells maintain long-term fitness during chronic stress and suggest that CN promotes adaptation in challenging environments by multiple mechanisms.

p53 controls choice between apoptotic and non-apoptotic death following DNA damage.

DNA damage can activate apoptotic and non-apoptotic forms of cell death; however, it remains unclear what features dictate which type of cell death is activated. We report that p53 controls the choice between apoptotic and non-apoptotic death following exposure to DNA damage. In contrast to the conventional model, which suggests that p53-deficient cells should be resistant to DNA damage-induced cell death, we find that p53-deficient cells die at high rates following DNA damage, but exclusively using non-apoptotic mechanisms. Our experimental data and computational modeling reveal that non-apoptotic death in p53-deficient cells has not been observed due to use of assays that are either insensitive to cell death, or that specifically score apoptotic cells. Using functional genetic screening - with an analysis that enables computational inference of the drug-induced death rate - we find in p53-deficient cells that DNA damage activates a mitochondrial respiration-dependent form of cell death, called MPT-driven necrosis. Cells deficient for p53 have high basal respiration, which primes MPT-driven necrosis. Finally, using metabolite profiling, we identified mitochondrial activity-dependent metabolic vulnerabilities that can be targeted to potentiate the lethality of DNA damage specifically in p53-deficient cells. Our findings reveal how the dual functions of p53 in regulating mitochondrial activity and the DNA damage response combine to facilitate the choice between apoptotic and non-apoptotic death.

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.

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.

Breast cancer prevention by short-term inhibition of TGFß signaling.

Cancer prevention has a profound impact on cancer-associated mortality and morbidity. We previously identified TGFß signaling as a candidate regulator of mammary epithelial cells associated with breast cancer risk. Here, we show that short-term TGFBR inhibitor (TGFBRi) treatment of peripubertal ACI inbred and Sprague Dawley outbred rats induces lasting changes and prevents estrogen- and carcinogen-induced mammary tumors, respectively. We identify TGFBRi-responsive cell populations by single cell RNA-sequencing, including a unique epithelial subpopulation designated secretory basal cells (SBCs) with progenitor features. We detect SBCs in normal human breast tissues and find them to be associated with breast cancer risk. Interactome analysis identifies SBCs as the most interactive cell population and the main source of insulin-IGF signaling. Accordingly, inhibition of TGFBR and IGF1R decrease proliferation of organoid cultures. Our results reveal a critical role for TGFß in regulating mammary epithelial cells relevant to breast cancer and serve as a proof-of-principle cancer prevention strategy.

Insights into Immune Escape During Tumor Evolution and Response to Immunotherapy Using a Rat Model of Breast Cancer.

Animal models are critical for the preclinical validation of cancer immunotherapies. Unfortunately, mouse breast cancer models do not faithfully reproduce the molecular subtypes and immune environment of the human disease. In particular, there are no good murine models of estrogen receptor-positive (ER+) breast cancer, the predominant subtype in patients. Here, we show that Nitroso-N-methylurea-induced mammary tumors in outbred Sprague-Dawley rats recapitulate the heterogeneity for mutational profiles, ER expression, and immune evasive mechanisms observed in human breast cancer. We demonstrate the utility of this model for preclinical studies by dissecting mechanisms of response to immunotherapy using combination TGFBR inhibition and PD-L1 blockade. Short-term treatment of early-stage tumors induced durable responses. Gene expression profiling and spatial mapping classified tumors as inflammatory and noninflammatory, and identified IFNγ, T-cell receptor (TCR), and B-cell receptor (BCR) signaling, CD74/MHC II, and epithelium-interacting CD8+ T cells as markers of response, whereas the complement system, M2 macrophage phenotype, and translation in mitochondria were associated with resistance. We found that the expression of CD74 correlated with leukocyte fraction and TCR diversity in human breast cancer. We identified a subset of rat ER+ tumors marked by expression of antigen-processing genes that had an active immune environment and responded to treatment. A gene signature characteristic of these tumors predicted disease-free survival in patients with ER+ Luminal A breast cancer and overall survival in patients with metastatic breast cancer receiving anti-PD-L1 therapy. We demonstrate the usefulness of this preclinical model for immunotherapy and suggest examination to expand immunotherapy to a subset of patients with ER+ disease. See related Spotlight by Roussos Torres, p. 672.

The impact of tumor epithelial and microenvironmental heterogeneity on treatment responses in HER2+ breast cancer.

Despite the availability of multiple human epidermal growth factor receptor 2-targeted (HER2-targeted) treatments, therapeutic resistance in HER2+ breast cancer remains a clinical challenge. Intratumor heterogeneity for HER2 and resistance-conferring mutations in the PIK3CA gene (encoding PI3K catalytic subunit α) have been investigated in response and resistance to HER2-targeting agents, while the role of divergent cellular phenotypes and tumor epithelial-stromal cell interactions is less well understood. Here, we assessed the effect of intratumor cellular genetic heterogeneity for ERBB2 (encoding HER2) copy number and PIK3CA mutation on different types of neoadjuvant HER2-targeting therapies and clinical outcome in HER2+ breast cancer. We found that the frequency of cells lacking HER2 was a better predictor of response to HER2-targeted treatment than intratumor heterogeneity. We also compared the efficacy of different therapies in the same tumor using patient-derived xenograft models of heterogeneous HER2+ breast cancer and single-cell approaches. Stromal determinants were better predictors of response than tumor epithelial cells, and we identified alveolar epithelial and fibroblastic reticular cells as well as lymphatic vessel endothelial hyaluronan receptor 1-positive (Lyve1+) macrophages as putative drivers of therapeutic resistance. Our results demonstrate that both preexisting and acquired resistance to HER2-targeting agents involve multiple mechanisms including the tumor microenvironment. Furthermore, our data suggest that intratumor heterogeneity for HER2 should be incorporated into treatment design.

ELP-dependent expression of MCL1 promotes resistance to EGFR inhibition in triple-negative breast cancer cells.

Targeted therapeutics for cancer generally exploit "oncogene addiction," a phenomenon in which the growth and survival of tumor cells depend on the activity of a particular protein. However, the efficacy of oncogene-targeted therapies varies substantially. For instance, targeting epidermal growth factor receptor (EGFR) signaling is effective in some non-small cell lung cancer (NSCLC) but not in triple-negative breast cancer (TNBC), although these cancers show a similar degree of increase in EGFR activity. Using a genome-wide CRISPR-Cas9 genetic knockout screen, we found that the Elongator (ELP) complex mediates insensitivity to the EGFR inhibitor erlotinib in TNBC cells by promoting the synthesis of the antiapoptotic protein Mcl-1. Depleting ELP proteins promoted apoptotic cell death in an EGFR inhibition-dependent manner. Pharmacological inhibition of Mcl-1 synergized with EGFR inhibition in a panel of genetically diverse TNBC cells. The findings indicate that TNBC "addiction" to EGFR signaling is masked by the ELP complex and that resistance to EGFR inhibitors in TNBC might be overcome by cotargeting Mcl-1.

EN1 Is a Transcriptional Dependency in Triple-Negative Breast Cancer Associated with Brain Metastasis.

To define transcriptional dependencies of triple-negative breast cancer (TNBC), we identified transcription factors highly and specifically expressed in primary TNBCs and tested their requirement for cell growth in a panel of breast cancer cell lines. We found that EN1 (engrailed 1) is overexpressed in TNBCs and its downregulation preferentially and significantly reduced viability and tumorigenicity in TNBC cell lines. By integrating gene expression changes after EN1 downregulation with EN1 chromatin binding patterns, we identified genes involved in WNT and Hedgehog signaling, neurogenesis, and axonal guidance as direct EN1 transcriptional targets. Quantitative proteomic analyses of EN1-bound chromatin complexes revealed association with transcriptional repressors and coactivators including TLE3, TRIM24, TRIM28, and TRIM33. High expression of EN1 correlated with short overall survival and increased risk of developing brain metastases in patients with TNBC. Thus, EN1 is a prognostic marker and a potential therapeutic target in TNBC. SIGNIFICANCE: These findings show that the EN1 transcription factor regulates neurogenesis-related genes and is associated with brain metastasis in triple-negative breast cancer.

Subclonal cooperation drives metastasis by modulating local and systemic immune microenvironments.

Most human tumours are heterogeneous, composed of cellular clones with different properties present at variable frequencies. Highly heterogeneous tumours have poor clinical outcomes, yet the underlying mechanism remains poorly understood. Here, we show that minor subclones of breast cancer cells expressing IL11 and FIGF (VEGFD) cooperate to promote metastatic progression and generate polyclonal metastases composed of driver and neutral subclones. Expression profiling of the epithelial and stromal compartments of monoclonal and polyclonal primary and metastatic lesions revealed that this cooperation is indirect, mediated through the local and systemic microenvironments. We identified neutrophils as a leukocyte population stimulated by the IL11-expressing minor subclone and showed that the depletion of neutrophils prevents metastatic outgrowth. Single-cell RNA-seq of CD45+ cell populations from primary tumours, blood and lungs demonstrated that IL11 acts on bone-marrow-derived mesenchymal stromal cells, which induce pro-tumorigenic and pro-metastatic neutrophils. Our results indicate key roles for non-cell-autonomous drivers and minor subclones in metastasis.

Perturbed myoepithelial cell differentiation in BRCA mutation carriers and in ductal carcinoma in situ.

Myoepithelial cells play key roles in normal mammary gland development and in limiting pre-invasive to invasive breast tumor progression, yet their differentiation and perturbation in ductal carcinoma in situ (DCIS) are poorly understood. Here, we investigated myoepithelial cells in normal breast tissues of BRCA1 and BRCA2 germline mutation carriers and in non-carrier controls, and in sporadic DCIS. We found that in the normal breast of non-carriers, myoepithelial cells frequently co-express the p63 and TCF7 transcription factors and that p63 and TCF7 show overlapping chromatin peaks associated with differentiated myoepithelium-specific genes. In contrast, in normal breast tissues of BRCA1 mutation carriers the frequency of p63+TCF7+ myoepithelial cells is significantly decreased and p63 and TCF7 chromatin peaks do not overlap. These myoepithelial perturbations in normal breast tissues of BRCA1 germline mutation carriers may play a role in their higher risk of breast cancer. The fraction of p63+TCF7+ myoepithelial cells is also significantly decreased in DCIS, which may be associated with invasive progression.