3D Culture Models with CRISPR Screens Reveal Hyperactive NRF2 as a Prerequisite for Spheroid Formation via Regulation of Proliferation and Ferroptosis.

Cancer-associated mutations that stabilize NRF2, an oxidant defense transcription factor, are predicted to promote tumor development. Here, utilizing 3D cancer spheroid models coupled with CRISPR-Cas9 screens, we investigate the molecular pathogenesis mediated by NRF2 hyperactivation. NRF2 hyperactivation was necessary for proliferation and survival in lung tumor spheroids. Antioxidant treatment rescued survival but not proliferation, suggesting the presence of distinct mechanisms. CRISPR screens revealed that spheroids are differentially dependent on the mammalian target of rapamycin (mTOR) for proliferation and the lipid peroxidase GPX4 for protection from ferroptosis of inner, matrix-deprived cells. Ferroptosis inhibitors blocked death from NRF2 downregulation, demonstrating a critical role of NRF2 in protecting matrix-deprived cells from ferroptosis. Interestingly, proteomics analyses show global enrichment of selenoproteins, including GPX4, by NRF2 downregulation, and targeting NRF2 and GPX4 killed spheroids overall. These results illustrate the value of spheroid culture in revealing environmental or spatial differential dependencies on NRF2 and reveal exploitable vulnerabilities of NRF2-hyperactivated tumors.

Simultaneous isolation of hormone receptor-positive breast cancer organoids and fibroblasts reveals stroma-mediated resistance mechanisms.

Recurrent hormone receptor-positive (HR+) breast cancer kills more than 600,000 women annually. Although HR+ breast cancers typically respond well to therapies, approximately 30% of patients relapse. At this stage, the tumors are usually metastatic and incurable. Resistance to therapy, particularly endocrine therapy is typically thought to be tumor intrinsic (e.g., estrogen receptor mutations). However, tumor-extrinsic factors also contribute to resistance. For example, stromal cells, such as cancer-associated fibroblasts (CAFs), residing in the tumor microenvironment, are known to stimulate resistance and disease recurrence. Recurrence in HR+ disease has been difficult to study due to the prolonged clinical course, complex nature of resistance, and lack of appropriate model systems. Existing HR+ models are limited to HR+ cell lines, a few HR+ organoid models, and xenograft models that all lack components of the human stroma. Therefore, there is an urgent need for more clinically relevant models to study the complex nature of recurrent HR+ breast cancer, and the factors contributing to treatment relapse. Here, we present an optimized protocol that allows a high take-rate, and simultaneous propagation of patient-derived organoids (PDOs) and matching CAFs, from primary and metastatic HR+ breast cancers. Our protocol allows for long-term culturing of HR+ PDOs that retain estrogen receptor expression and show responsiveness to hormone therapy. We further show the functional utility of this system by identifying CAF-secreted cytokines, such as growth-regulated oncogene α , as stroma-derived resistance drivers to endocrine therapy in HR+ PDOs.

Fibroblast-tumor cell signaling limits HER2 kinase therapy response via activation of MTOR and antiapoptotic pathways.

Despite the implementation of multiple HER2-targeted therapies, patients with advanced HER2+ breast cancer ultimately develop drug resistance. Stromal fibroblasts represent an abundant cell type in the tumor microenvironment and have been linked to poor outcomes and drug resistance. Here, we show that fibroblasts counteract the cytotoxic effects of HER2 kinase-targeted therapy in a subset of HER2+ breast cancer cell lines and allow cancer cells to proliferate in the presence of the HER2 kinase inhibitor lapatinib. Fibroblasts from primary breast tumors, normal breast tissue, and lung tissue have similar protective effects on tumor cells via paracrine factors. This fibroblast-mediated reduction in drug sensitivity involves increased expression of antiapoptotic proteins and sustained activation of the PI3K/AKT/MTOR pathway, despite inhibition of the HER2 and the RAS-ERK pathways in tumor cells. HER2 therapy sensitivity is restored in the fibroblast cocultures by combination treatment with inhibitors of MTOR or the antiapoptotic proteins BCL-XL and MCL-1. Expression of activated AKT in tumor cells recapitulates the effects of fibroblasts resulting in sustained MTOR signaling and poor lapatinib response. Lapatinib sensitivity was not altered by fibroblasts in tumor cells that exhibited sustained MTOR signaling due to a strong gain-of-function PI3KCA mutation. These findings indicate that in addition to tumor cell-intrinsic mechanisms that cause constitutive PI3K/AKT/MTOR pathway activation, secreted factors from fibroblasts can maintain this pathway in the context of HER2 inhibition. Our integrated proteomic-phenotypic approach presents a strategy for the discovery of protective mechanisms in fibroblast-rich tumors and the design of rational combination therapies to restore drug sensitivity.

Identification of cancer genes that are independent of dominant proliferation and lineage programs.

Large, multidimensional cancer datasets provide a resource that can be mined to identify candidate therapeutic targets for specific subgroups of tumors. Here, we analyzed human breast cancer data to identify transcriptional programs associated with tumors bearing specific genetic driver alterations. Using an unbiased approach, we identified thousands of genes whose expression was enriched in tumors with specific genetic alterations. However, expression of the vast majority of these genes was not enriched if associations were analyzed within individual breast tumor molecular subtypes, across multiple tumor types, or after gene expression was normalized to account for differences in proliferation or tumor lineage. Together with linear modeling results, these findings suggest that most transcriptional programs associated with specific genetic alterations in oncogenes and tumor suppressors are highly context-dependent and are predominantly linked to differences in proliferation programs between distinct breast cancer subtypes. We demonstrate that such proliferation-dependent gene expression dominates tumor transcriptional programs relative to matched normal tissues. However, we also identified a relatively small group of cancer-associated genes that are both proliferation- and lineage-independent. A subset of these genes are attractive candidate targets for combination therapy because they are essential in breast cancer cell lines, druggable, enriched in stem-like breast cancer cells, and resistant to chemotherapy-induced down-regulation.

Identification of genes that regulate epithelial cell migration using an siRNA screening approach.

To provide a systematic analysis of genes that regulate epithelial cell migration, we performed a high throughput wound healing screen with MCF-10A breast epithelial cells, using siRNAs targeting 1,081 human genes encoding phosphatases, kinases and proteins predicted to influence cell migration and adhesion. The primary screen identified three categories of hits: those that accelerate, those that inhibit and those that impair migration with associated effects on cell proliferation or metabolism. Extensive validation of all the hits yielded 66 high confidence genes that, when downregulated, either accelerated or impaired migration; 42 of these high confidence genes have not been previously associated with motility or adhesion. Time-lapse video microscopy revealed a broad spectrum of phenotypic changes involving alterations in the extent and nature of disruption of cell-cell adhesion, directionality of motility, cell polarity and shape, and protrusion dynamics. Informatics analysis highlighted three major signalling nodes, beta-catenin, beta1-integrin and actin, and a large proportion of the genes that accelerated migration impaired cell-cell adhesion.

Glutathione synthesis in the mouse liver supports lipid abundance through NRF2 repression.

Cells rely on antioxidants to survive. The most abundant antioxidant is glutathione (GSH). The synthesis of GSH is non-redundantly controlled by the glutamate-cysteine ligase catalytic subunit (GCLC). GSH imbalance is implicated in many diseases, but the requirement for GSH in adult tissues is unclear. To interrogate this, we have developed a series of in vivo models to induce Gclc deletion in adult animals. We find that GSH is essential to lipid abundance in vivo. GSH levels are highest in liver tissue, which is also a hub for lipid production. While the loss of GSH does not cause liver failure, it decreases lipogenic enzyme expression, circulating triglyceride levels, and fat stores. Mechanistically, we find that GSH promotes lipid abundance by repressing NRF2, a transcription factor induced by oxidative stress. These studies identify GSH as a fulcrum in the liver's balance of redox buffering and triglyceride production.

Evidence for a multipotent mammary progenitor with pregnancy-specific activity.

INTRODUCTION: The mouse mammary gland provides a powerful model system for studying processes involved in epithelial tissue development. Although markers that enrich for mammary stem cells and progenitors have been identified, our understanding of the mammary developmental hierarchy remains incomplete. METHODS: We used the MMTV promoter linked to the reverse tetracycline transactivator to induce H2BGFP expression in the mouse mammary gland. Mammary epithelial cells (MECs) from virgin mice were sorted by flow cytometry for expression of the mammary stem cell/progenitor markers CD24 and CD29, and H2BGFP. Sorted populations were analyzed for in vivo repopulation ability, expression of mammary lineage markers, and differential gene expression. RESULTS: The reconstituting activity of CD24⁺/CD29⁺ cells in cleared fat pad transplantation assays was not distinguished in GFP⁺ compared to GFP⁻ subpopulations. However, within the CD24⁺/CD29(lo) luminal progenitor-enriched population, H2BGFP⁺, but not H2BGFP⁻, MECs formed mammary structures in transplantation assays; moreover, this activity was dramatically enhanced in pregnant recipients. These outgrowths contained luminal and myoepithelial mammary lineages and produced milk, but lacked the capacity for serial transplantation. Transcriptional microarray analysis revealed that H2BGFP⁺/CD24⁺/CD29(lo) MECs are distinct from H2BGFP⁻/CD24⁺/CD29(lo) MECs and enriched for gene expression signatures with both the stem cell (CD24⁺/CD29⁺) and luminal progenitor (CD24⁺/CD29(lo)/CD61⁺) compartments. CONCLUSIONS: We have identified a population of MECs containing pregnancy-activated multipotent progenitors that are present in the virgin mammary gland and contribute to the expansion of the mammary gland during pregnancy.

Dose-dependent induction of distinct phenotypic responses to Notch pathway activation in mammary epithelial cells.

Aberrant activation of Notch receptors has been implicated in breast cancer; however, the mechanisms contributing to Notch-dependent transformation remain elusive because Notch displays dichotomous functional activities, promoting both proliferation and growth arrest. We investigated the cellular basis for the heterogeneous responses to Notch pathway activation in 3D cultures of MCF-10A mammary epithelial cells. Expression of a constitutively active Notch-1 intracellular domain (NICD) was found to induce two distinct types of 3D structures: large, hyperproliferative structures and small, growth-arrested structures with reduced cell-to-matrix adhesion. Interestingly, we found that these heterogeneous phenotypes reflect differences in Notch pathway activation levels; high Notch activity caused down-regulation of multiple matrix-adhesion genes and inhibition of proliferation, whereas low Notch activity maintained matrix adhesion and provoked a strong hyperproliferative response. Moreover, microarray analyses implicated NICD-induced p63 down-regulation in loss of matrix adhesion. In addition, a reverse-phase protein array-based analysis and subsequent loss-of-function studies identified STAT3 as a dominant downstream mediator of the NICD-induced outgrowth. These results indicate that the phenotypic responses to Notch are determined by the dose of pathway activation; and this dose affects the balance between growth-stimulative and growth-suppressive effects. This unique feature of Notch signaling provides insights into mechanisms that contribute to the dichotomous effects of Notch during development and tumorigenesis.

Glutathione supports lipid abundance in vivo.

Cells rely on antioxidants to survive. The most abundant antioxidant is glutathione (GSH). The synthesis of GSH is non-redundantly controlled by the glutamate-cysteine ligase catalytic subunit (GCLC). GSH imbalance is implicated in many diseases, but the requirement for GSH in adult tissues is unclear. To interrogate this, we developed a series of in vivo models to induce Gclc deletion in adult animals. We find that GSH is essential to lipid abundance in vivo. GSH levels are reported to be highest in liver tissue, which is also a hub for lipid production. While the loss of GSH did not cause liver failure, it decreased lipogenic enzyme expression, circulating triglyceride levels, and fat stores. Mechanistically, we found that GSH promotes lipid abundance by repressing NRF2, a transcription factor induced by oxidative stress. These studies identify GSH as a fulcrum in the liver's balance of redox buffering and triglyceride production.

Statistical methods for analysis of high-throughput RNA interference screens.

RNA interference (RNAi) has become a powerful technique for reverse genetics and drug discovery, and in both of these areas large-scale high-throughput RNAi screens are commonly performed. The statistical techniques used to analyze these screens are frequently borrowed directly from small-molecule screening; however, small-molecule and RNAi data characteristics differ in meaningful ways. We examine the similarities and differences between RNAi and small-molecule screens, highlighting particular characteristics of RNAi screen data that must be addressed during analysis. Additionally, we provide guidance on selection of analysis techniques in the context of a sample workflow.

Lineage-specific silencing of PSAT1 induces serine auxotrophy and sensitivity to dietary serine starvation in luminal breast tumors.

A major challenge of targeting metabolism for cancer therapy is pathway redundancy, in which multiple sources of critical nutrients can limit the effectiveness of some metabolism-targeted therapies. Here, we analyze lineage-dependent gene expression in human breast tumors to identify differences in metabolic gene expression that may limit pathway redundancy and create therapeutic vulnerabilities. We find that the serine synthesis pathway gene PSAT1 is the most depleted metabolic gene in luminal breast tumors relative to basal tumors. Low PSAT1 prevents de novo serine biosynthesis and sensitizes luminal breast cancer cells to serine and glycine starvation in vitro and in vivo. This PSAT1 expression disparity preexists in the putative cells of origin of basal and luminal tumors and is due to luminal-specific hypermethylation of the PSAT1 gene. Our data demonstrate that luminal breast tumors are auxotrophic for serine and may be uniquely sensitive to therapies targeting serine availability.

Clonal populations of a human TNBC model display significant functional heterogeneity and divergent growth dynamics in distinct contexts.

Intratumoral heterogeneity has been described for various tumor types and models of human cancer, and can have profound effects on tumor progression and drug resistance. This study describes an in-depth analysis of molecular and functional heterogeneity among subclonal populations (SCPs) derived from a single triple-negative breast cancer cell line, including copy number analysis, whole-exome and RNA sequencing, proteome analysis, and barcode analysis of clonal dynamics, as well as functional assays. The SCPs were found to have multiple unique genetic alterations and displayed significant variation in anchorage independent growth and tumor forming ability. Analyses of clonal dynamics in SCP mixtures using DNA barcode technology revealed selection for distinct clonal populations in different in vitro and in vivo environmental contexts, demonstrating that in vitro propagation of cancer cell lines using different culture conditions can contribute to the establishment of unique strains. These analyses also revealed strong enrichment of a single SCP during the development of xenograft tumors in immune-compromised mice. This SCP displayed attenuated interferon signaling in vivo and reduced sensitivity to the antiproliferative effects of type I interferons. Reduction in interferon signaling was found to provide a selective advantage within the xenograft microenvironment specifically. In concordance with the previously described role of interferon signaling as tumor suppressor, these findings suggest that similar selective pressures may be operative in human cancer and patient-derived xenograft models.

A human breast atlas integrating single-cell proteomics and transcriptomics.

The breast is a dynamic organ whose response to physiological and pathophysiological conditions alters its disease susceptibility, yet the specific effects of these clinical variables on cell state remain poorly annotated. We present a unified, high-resolution breast atlas by integrating single-cell RNA-seq, mass cytometry, and cyclic immunofluorescence, encompassing a myriad of states. We define cell subtypes within the alveolar, hormone-sensing, and basal epithelial lineages, delineating associations of several subtypes with cancer risk factors, including age, parity, and BRCA2 germline mutation. Of particular interest is a subset of alveolar cells termed basal-luminal (BL) cells, which exhibit poor transcriptional lineage fidelity, accumulate with age, and carry a gene signature associated with basal-like breast cancer. We further utilize a medium-depletion approach to identify molecular factors regulating cell-subtype proportion in organoids. Together, these data are a rich resource to elucidate diverse mammary cell states.

Clinical evaluation of BCL-2/XL levels pre- and post- HER2-targeted therapy.

Our previous pre-clinical work defined BCL-2 induction as a critical component of the adaptive response to lapatinib-mediated inhibition of HER2. To determine whether a similar BCL-2 upregulation occurs in lapatinib-treated patients, we evaluated gene expression within tumor biopsies, collected before and after lapatinib or trastuzumab treatment, from the TRIO-B-07 clinical trial (NCT#00769470). We detected BCL2 mRNA upregulation in both HER2+/ER- as well as HER2+/ER+ patient tumors treated with lapatinib or trastuzumab. To address whether mRNA expression correlated with protein expression, we evaluated pre- and post-treatment tumors for BCL-2 via immunohistochemistry. Despite BCL2 mRNA upregulation within HER2+/ER- tumors, BCL-2 protein levels were undetectable in most of the lapatinib- or trastuzumab-treated HER2+/ER- tumors. BCL-2 upregulation was evident within the majority of lapatinib-treated HER2+/ER+ tumors and was often coupled with increased ER expression and decreased proliferation. Comparable BCL-2 upregulation was not observed within the trastuzumab-treated HER2+/ER+ tumors. Together, these results provide clinical validation of the BCL-2 induction associated with the adaptive response to lapatinib and support evaluation of BCL-2 inhibitors within the context of lapatinib and other HER2-targeted receptor tyrosine kinase inhibitors.

Screensaver: an open source lab information management system (LIMS) for high throughput screening facilities.

BACKGROUND: Shared-usage high throughput screening (HTS) facilities are becoming more common in academe as large-scale small molecule and genome-scale RNAi screening strategies are adopted for basic research purposes. These shared facilities require a unique informatics infrastructure that must not only provide access to and analysis of screening data, but must also manage the administrative and technical challenges associated with conducting numerous, interleaved screening efforts run by multiple independent research groups. RESULTS: We have developed Screensaver, a free, open source, web-based lab information management system (LIMS), to address the informatics needs of our small molecule and RNAi screening facility. Screensaver supports the storage and comparison of screening data sets, as well as the management of information about screens, screeners, libraries, and laboratory work requests. To our knowledge, Screensaver is one of the first applications to support the storage and analysis of data from both genome-scale RNAi screening projects and small molecule screening projects. CONCLUSIONS: The informatics and administrative needs of an HTS facility may be best managed by a single, integrated, web-accessible application such as Screensaver. Screensaver has proven useful in meeting the requirements of the ICCB-Longwood/NSRB Screening Facility at Harvard Medical School, and has provided similar benefits to other HTS facilities.

Aging-Associated Alterations in Mammary Epithelia and Stroma Revealed by Single-Cell RNA Sequencing.

Aging is closely associated with increased susceptibility to breast cancer, yet there have been limited systematic studies of aging-induced alterations in the mammary gland. Here, we leverage high-throughput single-cell RNA sequencing to generate a detailed transcriptomic atlas of young and aged murine mammary tissues. By analyzing epithelial, stromal, and immune cells, we identify age-dependent alterations in cell proportions and gene expression, providing evidence that suggests alveolar maturation and physiological decline. The analysis also uncovers potential pro-tumorigenic mechanisms coupled to the age-associated loss of tumor suppressor function and change in microenvironment. In addition, we identify a rare, age-dependent luminal population co-expressing hormone-sensing and secretory-alveolar lineage markers, as well as two macrophage populations expressing distinct gene signatures, underscoring the complex heterogeneity of the mammary epithelia and stroma. Collectively, this rich single-cell atlas reveals the effects of aging on mammary physiology and can serve as a useful resource for understanding aging-associated cancer risk.

Organoid cultures from normal and cancer-prone human breast tissues preserve complex epithelial lineages.

Recently, organoid technology has been used to generate a large repository of breast cancer organoids. Here we present an extensive evaluation of the ability of organoid culture technology to preserve complex stem/progenitor and differentiated cell types via long-term propagation of normal human mammary tissues. Basal/stem and luminal progenitor cells can differentiate in culture to generate mature basal and luminal cell types, including ER+ cells that have been challenging to maintain in culture. Cells associated with increased cancer risk can also be propagated. Single-cell analyses of matched organoid cultures and native tissues by mass cytometry for 38 markers provide a higher resolution representation of the multiple mammary epithelial cell types in the organoids, and demonstrate that protein expression patterns of the tissue of origin can be preserved in culture. These studies indicate that organoid cultures provide a valuable platform for studies of mammary differentiation, transformation, and breast cancer risk.

Transient commensal clonal interactions can drive tumor metastasis.

The extent and importance of functional heterogeneity and crosstalk between tumor cells is poorly understood. Here, we describe the generation of clonal populations from a patient-derived ovarian clear cell carcinoma model which forms malignant ascites and solid peritoneal tumors upon intraperitoneal transplantation in mice. The clonal populations are engineered with secreted Gaussia luciferase to monitor tumor growth dynamics and tagged with a unique DNA barcode to track their fate in multiclonal mixtures during tumor progression. Only one clone, CL31, grows robustly, generating exclusively malignant ascites. However, multiclonal mixtures form large solid peritoneal metastases, populated almost entirely by CL31, suggesting that transient cooperative interclonal interactions are sufficient to promote metastasis of CL31. CL31 uniquely harbors ERBB2 amplification, and its acquired metastatic activity in clonal mixtures is dependent on transient exposure to amphiregulin, which is exclusively secreted by non-tumorigenic clones. Amphiregulin enhances CL31 mesothelial clearance, a prerequisite for metastasis. These findings demonstrate that transient, ostensibly innocuous tumor subpopulations can promote metastases via "hit-and-run" commensal interactions.

Direct stimulation of NADP+ synthesis through Akt-mediated phosphorylation of NAD kinase.

Nicotinamide adenine dinucleotide phosphate (NADP+) is essential for producing NADPH, the primary cofactor for reductive metabolism. We find that growth factor signaling through the phosphoinositide 3-kinase (PI3K)-Akt pathway induces acute synthesis of NADP+ and NADPH. Akt phosphorylates NAD kinase (NADK), the sole cytosolic enzyme that catalyzes the synthesis of NADP+ from NAD+ (the oxidized form of NADH), on three serine residues (Ser44, Ser46, and Ser48) within an amino-terminal domain. This phosphorylation stimulates NADK activity both in cells and directly in vitro, thereby increasing NADP+ production. A rare isoform of NADK (isoform 3) lacking this regulatory region exhibits constitutively increased activity. These data indicate that Akt-mediated phosphorylation of NADK stimulates its activity to increase NADP+ production through relief of an autoinhibitory function inherent to its amino terminus.

Combined MEK and BCL-2/XL Inhibition Is Effective in High-Grade Serous Ovarian Cancer Patient-Derived Xenograft Models and BIM Levels Are Predictive of Responsiveness.

Most patients with late-stage high-grade serous ovarian cancer (HGSOC) initially respond to chemotherapy but inevitably relapse and develop resistance, highlighting the need for novel therapies to improve patient outcomes. The MEK/ERK pathway is activated in a large subset of HGSOC, making it an attractive therapeutic target. Here, we systematically evaluated the extent of MEK/ERK pathway activation and efficacy of pathway inhibition in a large panel of well-annotated HGSOC patient-derived xenograft models. The vast majority of models were nonresponsive to the MEK inhibitor cobimetinib (GDC-0973) despite effective pathway inhibition. Proteomic analyses of adaptive responses to GDC-0973 revealed that GDC-0973 upregulated the proapoptotic protein BIM, thus priming the cells for apoptosis regulated by BCL2-family proteins. Indeed, combination of both MEK inhibitor and dual BCL-2/XL inhibitor (ABT-263) significantly reduced cell number, increased cell death, and displayed synergy in vitro in most models. In vivo, GDC-0973 and ABT-263 combination was well tolerated and resulted in greater tumor growth inhibition than single agents. Detailed proteomic and correlation analyses identified two subsets of responsive models-those with high BIM at baseline that was increased with MEK inhibition and those with low basal BIM and high pERK levels. Models with low BIM and low pERK were nonresponsive. Our findings demonstrate that combined MEK and BCL-2/XL inhibition has therapeutic activity in HGSOC models and provide a mechanistic rationale for the clinical evaluation of this drug combination as well as the assessment of the extent to which BIM and/or pERK levels predict drug combination effectiveness in chemoresistant HGSOC.