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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Deubiquitinases Maintain Protein Homeostasis and Survival of Cancer Cells upon Glutathione Depletion.

Cells are subjected to oxidative stress during the initiation and progression of tumors, and this imposes selective pressure for cancer cells to adapt mechanisms to tolerate these conditions. Here, we examined the dependency of cancer cells on glutathione (GSH), the most abundant cellular antioxidant. While cancer cell lines displayed a broad range of sensitivities to inhibition of GSH synthesis, the majority were resistant to GSH depletion. To identify cellular pathways required for this resistance, we carried out genetic and pharmacologic screens. Both approaches revealed that inhibition of deubiquitinating enzymes (DUBs) sensitizes cancer cells to GSH depletion. Inhibition of GSH synthesis, in combination with DUB inhibition, led to an accumulation of polyubiquitinated proteins, induction of proteotoxic stress, and cell death. These results indicate that depletion of GSH renders cancer cells dependent on DUB activity to maintain protein homeostasis and cell viability and reveal a potentially exploitable vulnerability for cancer therapy.

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.

CRB3 and the FERM protein EPB41L4B regulate proliferation of mammary epithelial cells through the release of amphiregulin.

Numerous observations have suggested a connection between the maintenance of cell polarity and control of cell proliferation; however, the mechanisms underlying these connections remain poorly understood. Here we found that ectopic expression of CRB3, which was previously shown to restore tight junctions and membrane polarity in MCF-10A cells, induced a hyperproliferative phenotype, with significantly enlarged acini in basement membrane culture, similar to structures induced by expression of proliferative oncogenes such as cyclinD1. We found that CRB3-induced proliferation is epidermal growth factor (EGF)-independent and occurs through a mechanism that involves secretion of the EGF-family ligand, amphiregulin (AREG). The increase in AREG secretion is associated with an increase in the number and size of both early and late endosomes. Both the proliferative and endocytic phenotypes associated with CRB3 expression require the FERM-binding domain (FBD) but not the PDZ-binding domain of CRB3, arguing that this proliferative phenotype is independent of the PDZ-dependent polarity signaling by CRB3. We identified the FBD-containing protein, EPB41L4B, as an essential mediator of CRB3-driven proliferation and observed that the CRB3-dependent changes in endocytic trafficking were also dependent on EPB41L4B. Taken together, these data reveal a previously uncharacterized role for CRB3 in regulating proliferation in mammalian cells that is associated with changes in the endocytic trafficking machinery.

Cancer Cells Co-opt the Neuronal Redox-Sensing Channel TRPA1 to Promote Oxidative-Stress Tolerance.

Cancer cell survival is dependent on oxidative-stress defenses against reactive oxygen species (ROS) that accumulate during tumorigenesis. Here, we show a non-canonical oxidative-stress defense mechanism through TRPA1, a neuronal redox-sensing Ca2+-influx channel. In TRPA1-enriched breast and lung cancer spheroids, TRPA1 is critical for survival of inner cells that exhibit ROS accumulation. Moreover, TRPA1 promotes resistance to ROS-producing chemotherapies, and TRPA1 inhibition suppresses xenograft tumor growth and enhances chemosensitivity. TRPA1 does not affect redox status but upregulates Ca2+-dependent anti-apoptotic pathways. NRF2, an oxidant-defense transcription factor, directly controls TRPA1 expression, thus providing an orthogonal mechanism for protection against oxidative stress together with canonical ROS-neutralizing mechanisms. These findings reveal an oxidative-stress defense program involving TRPA1 that could be exploited for targeted cancer therapies.

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.

Erratum: Niche-localized tumor cells are protected from HER2-targeted therapy via upregulation of an anti-apoptotic program in vivo.

[This corrects the article DOI: 10.1038/s41523-017-0020-z.].

Systems analysis of apoptotic priming in ovarian cancer identifies vulnerabilities and predictors of drug response.

The lack of effective chemotherapies for high-grade serous ovarian cancers (HGS-OvCa) has motivated a search for alternative treatment strategies. Here, we present an unbiased systems-approach to interrogate a panel of 14 well-annotated HGS-OvCa patient-derived xenografts for sensitivity to PI3K and PI3K/mTOR inhibitors and uncover cell death vulnerabilities. Proteomic analysis reveals that PI3K/mTOR inhibition in HGS-OvCa patient-derived xenografts induces both pro-apoptotic and anti-apoptotic signaling responses that limit cell killing, but also primes cells for inhibitors of anti-apoptotic proteins. In-depth quantitative analysis of BCL-2 family proteins and other apoptotic regulators, together with computational modeling and selective anti-apoptotic protein inhibitors, uncovers new mechanistic details about apoptotic regulators that are predictive of drug sensitivity (BIM, caspase-3, BCL-XL) and resistance (MCL-1, XIAP). Our systems-approach presents a strategy for systematic analysis of the mechanisms that limit effective tumor cell killing and the identification of apoptotic vulnerabilities to overcome drug resistance in ovarian and other cancers.High-grade serous ovarian cancers (HGS-OvCa) frequently develop chemotherapy resistance. Here, the authors through a systematic analysis of proteomic and drug response data of 14 HGS-OvCa PDXs demonstrate that targeting apoptosis regulators can improve response of these tumors to inhibitors of the PI3K/mTOR pathway.

Niche-localized tumor cells are protected from HER2-targeted therapy via upregulation of an anti-apoptotic program in vivo.

Several lines of evidence suggest that components of the tumor microenvironment, specifically basement membrane and extracellular matrix proteins, influence drug sensitivities. We previously reported differential drug sensitivity of tumor cells localized adjacent to laminin-rich extracellular matrix in three-dimensional tumor spheroid cultures. To evaluate whether differential intra-tumor responses to targeted therapy occur in vivo, we examined the sensitivity of human epidermal growth factor receptor 2-positive tumors to lapatinib using a previously described ductal carcinoma in situ-like model characterized by tumor cell confinement within ductal structures surrounded by an organized basement membrane. Here we show that tumor cells localized to a 'niche' in the outer layer of the intraductal tumors adjacent to myoepithelial cells and basement membrane are resistant to lapatinib. We found that the pro-survival protein BCL2 is selectively induced in the niche-protected tumor cells following lapatinib treatment, and combined inhibition of HER2 and BCL-2/XL enhanced targeting of these residual tumor cells. Elimination of the niche-protected tumor cells was achieved with the HER2 antibody-drug conjugate T-DM1, which delivers a chemotherapeutic payload. Thus, these studies provide evidence that subpopulations of tumor cells within specific microenvironmental niches can adapt to inhibition of critical oncogenic pathways, and furthermore reveal effective strategies to eliminate these resistant subpopulations.