excluderanges: exclusion sets for T2T-CHM13, GRCm39, and other genome assemblies.

SUMMARY: Exclusion regions are sections of reference genomes with abnormal pileups of short sequencing reads. Removing reads overlapping them improves biological signal, and these benefits are most pronounced in differential analysis settings. Several labs created exclusion region sets, available primarily through ENCODE and Github. However, the variety of exclusion sets creates uncertainty which sets to use. Furthermore, gap regions (e.g. centromeres, telomeres, short arms) create additional considerations in generating exclusion sets. We generated exclusion sets for the latest human T2T-CHM13 and mouse GRCm39 genomes and systematically assembled and annotated these and other sets in the excluderanges R/Bioconductor data package, also accessible via the BEDbase.org API. The package provides unified access to 82 GenomicRanges objects covering six organisms, multiple genome assemblies, and types of exclusion regions. For human hg38 genome assembly, we recommend hg38.Kundaje.GRCh38_unified_blacklist as the most well-curated and annotated, and sets generated by the Blacklist tool for other organisms. AVAILABILITY AND IMPLEMENTATION: https://bioconductor.org/packages/excluderanges/. Package website: https://dozmorovlab.github.io/excluderanges/.

Erratum: Asparagine bioavailability governs metastasis in a model of breast cancer.

This corrects the article DOI: 10.1038/nature25465.

Asparagine bioavailability governs metastasis in a model of breast cancer.

Using a functional model of breast cancer heterogeneity, we previously showed that clonal sub-populations proficient at generating circulating tumour cells were not all equally capable of forming metastases at secondary sites. A combination of differential expression and focused in vitro and in vivo RNA interference screens revealed candidate drivers of metastasis that discriminated metastatic clones. Among these, asparagine synthetase expression in a patient's primary tumour was most strongly correlated with later metastatic relapse. Here we show that asparagine bioavailability strongly influences metastatic potential. Limiting asparagine by knockdown of asparagine synthetase, treatment with l-asparaginase, or dietary asparagine restriction reduces metastasis without affecting growth of the primary tumour, whereas increased dietary asparagine or enforced asparagine synthetase expression promotes metastatic progression. Altering asparagine availability in vitro strongly influences invasive potential, which is correlated with an effect on proteins that promote the epithelial-to-mesenchymal transition. This provides at least one potential mechanism for how the bioavailability of a single amino acid could regulate metastatic progression.

A model of breast cancer heterogeneity reveals vascular mimicry as a driver of metastasis.

Cancer metastasis requires that primary tumour cells evolve the capacity to intravasate into the lymphatic system or vasculature, and extravasate into and colonize secondary sites. Others have demonstrated that individual cells within complex populations show heterogeneity in their capacity to form secondary lesions. Here we develop a polyclonal mouse model of breast tumour heterogeneity, and show that distinct clones within a mixed population display specialization, for example, dominating the primary tumour, contributing to metastatic populations, or showing tropism for entering the lymphatic or vasculature systems. We correlate these stable properties to distinct gene expression profiles. Those clones that efficiently enter the vasculature express two secreted proteins, Serpine2 and Slpi, which were necessary and sufficient to program these cells for vascular mimicry. Our data indicate that these proteins not only drive the formation of extravascular networks but also ensure their perfusion by acting as anticoagulants. We propose that vascular mimicry drives the ability of some breast tumour cells to contribute to distant metastases while simultaneously satisfying a critical need of the primary tumour to be fed by the vasculature. Enforced expression of SERPINE2 and SLPI in human breast cancer cell lines also programmed them for vascular mimicry, and SERPINE2 and SLPI were overexpressed preferentially in human patients that had lung-metastatic relapse. Thus, these two secreted proteins, and the phenotype they promote, may be broadly relevant as drivers of metastatic progression in human cancer.

New generation breast cancer cell lines developed from patient-derived xenografts.

BACKGROUND: Breast cancer is a highly heterogeneous disease characterized by multiple histologic and molecular subtypes. While a myriad of breast cancer cell lines have been developed over the past 60 years, estrogen receptor alpha (ER)+ disease and some mutations associated with this subtype remain underrepresented. Here we describe six breast cancer cell lines derived from patient-derived xenografts (PDX) and their general characteristics. METHODS: Established breast cancer PDX were processed into cell suspensions and placed into standard 2D cell culture; six emerged into long-term passageable cell lines. Cell lines were assessed for protein expression of common luminal, basal, and mesenchymal markers, growth assessed in response to estrogens and endocrine therapies, and RNA-seq and oncogenomics testing performed to compare relative transcript levels and identify putative oncogenic drivers. RESULTS: Three cell lines express ER and two are also progesterone receptor (PR) positive; PAM50 subtyping identified one line as luminal A. One of the ER+PR+ lines harbors a D538G mutation in the gene for ER (ESR1), providing a natural model that contains this endocrine-resistant genotype. The third ER+PR-/low cell line has mucinous features, a rare histologic type of breast cancer. The three other lines are ER- and represent two basal-like and a mixed ductal/lobular breast cancer. The cell lines show varied responses to tamoxifen and fulvestrant, and three were demonstrated to regrow tumors in vivo. RNA sequencing confirms all cell lines are human and epithelial. Targeted oncogenomics testing confirmed the noted ESR1 mutation in addition to other mutations (i.e., PIK3CA, BRCA2, CCND1, NF1, TP53, MYC) and amplifications (i.e., FGFR1, FGFR3) frequently found in breast cancers. CONCLUSIONS: These new generation breast cancer cell lines add to the existing repository of breast cancer models, increase the number of ER+ lines, and provide a resource that can be genetically modified for studying several important clinical breast cancer features.

The long-term prognostic and predictive capacity of cyclin D1 gene amplification in 2305 breast tumours.

BACKGROUND: Use of cyclin D1 (CCND1) gene amplification as a breast cancer biomarker has been hampered by conflicting assessments of the relationship between cyclin D1 protein levels and patient survival. Here, we aimed to clarify its prognostic and treatment predictive potential through comprehensive long-term survival analyses. METHODS: CCND1 amplification was assessed using SNP arrays from two cohorts of 1965 and 340 patients with matching gene expression array and clinical follow-up data of over 15 years. Kaplan-Meier and multivariable Cox regression analyses were used to determine survival differences between CCND1 amplified vs. non-amplified tumours in clinically relevant patient sets, within PAM50 subtypes and within treatment-specific subgroups. Boxplots and differential gene expression analyses were performed to assess differences between amplified vs. non-amplified tumours within PAM50 subtypes. RESULTS: When combining both cohorts, worse survival was found for patients with CCND1-amplified tumours in luminal A (HR = 1.68; 95% CI, 1.15-2.46), luminal B (1.37; 1.01-1.86) and ER+/LN-/HER2- (1.66; 1.14-2.41) subgroups. In gene expression analysis, CCND1-amplified luminal A tumours showed increased proliferation (P < 0.001) and decreased progesterone (P = 0.002) levels along with a large overlap in differentially expressed genes when comparing luminal A and B-amplified vs. non-amplified tumours. CONCLUSIONS: Our results indicate that CCND1 amplification is associated with worse 15-year survival in ER+/LN-/HER2-, luminal A and luminal B patients. Moreover, luminal A CCND1-amplified tumours display gene expression changes consistent with a more aggressive phenotype. These novel findings highlight the potential of CCND1 to identify patients that could benefit from long-term treatment strategies.

Separation of breast cancer and organ microenvironment transcriptomes in metastases.

BACKGROUND: The seed and soil hypothesis was proposed over a century ago to describe why cancer cells (seeds) grow in certain organs (soil). Since then, the genetic properties that define the cancer cells have been heavily investigated; however, genomic mediators within the organ microenvironment that mediate successful metastatic growth are less understood. These studies sought to identify cancer- and organ-specific genomic programs that mediate metastasis. METHODS: In these studies, a set of 14 human breast cancer patient-derived xenograft (PDX) metastasis models was developed and then tested for metastatic tropism with two approaches: spontaneous metastases from mammary tumors and intravenous injection of PDX cells. The transcriptomes of the cancer cells when growing as tumors or metastases were separated from the transcriptomes of the microenvironment via species-specific separation of the genomes. Drug treatment of PDX spheroids was performed to determine if genes activated in metastases may identify targetable mediators of viability. RESULTS: The experimental approaches that generated metastases in PDX models were identified. RNA sequencing of 134 tumors, metastases, and normal non-metastatic organs identified cancer- and organ-specific genomic properties that mediated metastasis. A common genomic response of the liver microenvironment was found to occur in reaction to the invading PDX cells. Genes within the cancer cells were found to be either transiently regulated by the microenvironment or permanently altered due to clonal selection of metastatic sublines. Gene Set Enrichment Analyses identified more than 400 gene signatures that were commonly activated in metastases across basal-like PDXs. A Src signaling signature was found to be extensively upregulated in metastases, and Src inhibitors were found to be cytotoxic to PDX spheroids. CONCLUSIONS: These studies identified that during the growth of breast cancer metastases, there were genomic changes that occurred within both the cancer cells and the organ microenvironment. We hypothesize that pathways upregulated in metastases are mediators of viability and that simultaneously targeting changes within different cancer cell pathways and/or different tissue compartments may be needed for inhibition of disease progression.

Sphingosine kinase 1 activation by estrogen receptor α36 contributes to tamoxifen resistance in breast cancer.

In breast cancer, 17ß-estradiol (E2) plays critical roles mainly by binding to its canonical receptor, estrogen receptor (ER) α66, and eliciting genomic effects. E2 also triggers rapid, nongenomic responses. E2 activates sphingosine kinase 1 (SphK1), increasing sphingosine-1-phosphate (S1P) that binds to its receptors, leading to important breast cancer signaling. However, the E2 receptor responsible for SphK1 activation has not yet been identified. Here, we demonstrate in triple-negative breast cancer cells, which lack the canonical ERα66 but express the novel splice variant ERα36, that ERα36 is the receptor responsible for E2-induced activation of SphK1 and formation and secretion of S1P and dihydro-S1P, the ligands for S1PRs. Tamoxifen, the first-line endocrine therapy for breast cancer, is an antagonist of ERα66, but an agonist of ERα36, and, like E2, activates SphK1 and markedly increases secretion of S1P. A major problem with tamoxifen therapy is development of acquired resistance. We found that tamoxifen resistance correlated with increased SphK1 and ERα36 expression in tamoxifen-resistant breast cancer cells, in patient-derived xenografts, and in endocrine-resistant breast cancer patients. Our data also indicate that targeting this ERα36 and SphK1 axis may be a therapeutic option to circumvent endocrine resistance and improve patient outcome.

Characterizing the efficacy of cancer therapeutics in patient-derived xenograft models of metastatic breast cancer.

PURPOSE: Basal-like breast cancers are aggressive and often metastasize to vital organs. Treatment is largely limited to chemotherapy. This study aims to characterize the efficacy of cancer therapeutics in vitro and in vivo within the primary tumor and metastatic setting, using patient-derived xenograft (PDX) models. METHODS: We employed two basal-like, triple-negative PDX models, WHIM2 and WHIM30. PDX cells, obtained from mammary tumors grown in mice, were treated with twelve cancer therapeutics to evaluate their cytotoxicity in vitro. Four of the effective drugs-carboplatin, cyclophosphamide, bortezomib, and dacarbazine-were tested in vivo for their efficacy in treating mammary tumors, and metastases generated by intracardiac injection of tumor cells. RESULTS: RNA sequencing showed that global gene expression of PDX cells grown in the mammary gland was similar to those tested in culture. In vitro, carboplatin was cytotoxic to WHIM30 but not WHIM2, whereas bortezomib, dacarbazine, and cyclophosphamide were cytotoxic to both lines. Yet, these drugs were ineffective in treating both primary and metastatic WHIM2 tumors in vivo. Carboplatin and cyclophosphamide were effective in treating WHIM30 mammary tumors and reducing metastatic burden in the brain, liver, and lungs. WHIM2 and WHIM30 metastases showed distinct patterns of cytokeratin and vimentin expression, regardless of treatment, suggesting that different tumor cell subpopulations may preferentially seed in different organs. CONCLUSIONS: This study highlights the utility of PDX models for studying the efficacy of therapeutics in reducing metastatic burden in specific organs. The differential treatment responses between two PDX models of the same intrinsic subtype, in both the primary and metastatic setting, recapitulates the challenges faced in treating cancer patients and highlights the need for combination therapies and predictive biomarkers.

CIB1 depletion impairs cell survival and tumor growth in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is an aggressive breast cancer subtype with generally poor prognosis and no available targeted therapies, highlighting a critical unmet need to identify and characterize novel therapeutic targets. We previously demonstrated that CIB1 is necessary for cancer cell survival and proliferation via regulation of two oncogenic signaling pathways, RAF-MEK-ERK and PI3K-AKT. Because these pathways are often upregulated in TNBC, we hypothesized that CIB1 may play a broader role in TNBC cell survival and tumor growth. Methods utilized include inducible RNAi depletion of CIB1 in vitro and in vivo, immunoblotting, clonogenic assay, flow cytometry, RNA-sequencing, bioinformatics analysis, and Kaplan-Meier survival analysis. CIB1 depletion resulted in significant cell death in 8 of 11 TNBC cell lines tested. Analysis of components related to PI3K-AKT and RAF-MEK-ERK signaling revealed that elevated AKT activation status and low PTEN expression were key predictors of sensitivity to CIB1 depletion. Furthermore, CIB1 knockdown caused dramatic shrinkage of MDA-MB-468 xenograft tumors in vivo. RNA sequence analysis also showed that CIB1 depletion in TNBC cells activates gene programs associated with decreased proliferation and increased cell death. CIB1 expression levels per se did not predict TNBC susceptibility to CIB1 depletion, and CIB1 mRNA expression levels did not associate with TNBC patient survival. Our data are consistent with the emerging concept of non-oncogene addiction, where a large subset of TNBCs depend on CIB1 for cell survival and tumor growth, independent of CIB1 expression levels. Our data establish CIB1 as a novel therapeutic target for TNBC.

Maintenance of hormone responsiveness in luminal breast cancers by suppression of Notch.

Luminal breast cancers express estrogen (ER) and/or progesterone (PR) receptors and respond to hormone therapies. Basal-like "triple negative" cancers lack steroid receptors but are cytokeratin (CK) 5-positive and require chemotherapy. Here we show that more than half of primary ER(+)PR(+) breast cancers contain an ER(-)PR(-)CK5(+) "luminobasal" subpopulation exceeding 1% of cells. Starting from ER(+)PR(+) luminal cell lines, we generated lines with varying luminal to luminobasal cell ratios and studied their molecular and biological properties. In luminal disease, luminobasal cells expand in response to antiestrogen or estrogen withdrawal therapies. The phenotype and gene signature of the hormone-resistant cells matches that of clinical triple negative basal-like and claudin-low disease. Luminobasal cell expansion in response to hormone therapies is regulated by Notch1 signaling and can be blocked by γ-secretase inhibitors. Our data establish a previously unrecognized plasticity of ER(+)PR(+) luminal breast cancers that, without genetic manipulation, mobilizes outgrowth of hormone-resistant basal-like disease in response to treatment. This undesirable outcome can be prevented by combining endocrine therapies with Notch inhibition.

Vascular endothelial growth factor C promotes breast cancer progression via a novel antioxidant mechanism that involves regulation of superoxide dismutase 3.

INTRODUCTION: Triple-negative breast cancers, particularly the claudin-low subtype, are highly aggressive and exhibit increased tumor-initiating cell (TIC) characteristics. In this study, we demonstrate that vascular endothelial growth factor C (VEGF-C) is highly expressed in the claudin-low breast cancer subtype and also that it mediates tumor progression, not only through its role in lymphangiogenesis but also through regulating TIC characteristics and the response to reactive oxygen species (ROS). METHODS: VEGF C expression was examined in breast cancer subtypes, and a VEGF C expression signature was derived. VEGF C expression and/or its associated signature was correlated with TIC and chemoresistance signatures. In vitro and in vivo assays were performed to determine whether VEGF-C expression alters TIC characteristics and the response of breast cancer cells to chemotherapy and oxidative stress. Array analysis was used to identify a downstream effector of VEGF-C, superoxide dismutase 3 (Sod3), which was tested for its involvement in VEGF-C-mediated resistance to oxidative stress and enhancement of in vivo metastasis. The VEGF-C-associated receptor neuropilin 2 (Nrp2) was knocked down to determine whether it is required for the observed effects of VEGF-C. Expression of VEGF C and Sod3 was assessed in human breast cancers. RESULTS: VEGF C is highly expressed in claudin-low breast cancers, and VEGF C and the VEGF C signature are associated with TIC-related gene signatures. VEGF-C-knockdown in mammary carcinoma cells decreases TIC properties in vitro and in vivo, sensitizing cells to oxidative stress and chemotherapy. We identified Sod3 as a target of VEGF-C in breast cancer cells by demonstrating that it is required for VEGF-C-mediated cell survival in response to oxidative stress and for VEGF-C-mediated metastasis. We demonstrate that Nrp2 is the VEGF-C-associated receptor that mediates alterations in Sod3 expression and the response of tumor cells to oxidative stress. We show that VEGF C and Sod3 are positively associated in human breast cancer. CONCLUSIONS: We describe a novel mechanism by which VEGF-C contributes to metastasis via its ability to enhance TIC-associated characteristics, particularly the response to ROS. We identified Sod3 as a critical mediator of VEGF-C-induced metastasis, and we provide evidence that the VEGF-C-Sod3 axis plays a role in human breast cancers.

Identification of a stable molecular signature in mammary tumor endothelial cells that persists in vitro.

Long-term, in vitro propagation of tumor-specific endothelial cells (TEC) allows for functional studies and genome-wide expression profiling of clonally derived, well-characterized subpopulations. Using a genetically engineered mouse model of mammary adenocarcinoma, we have optimized an isolation procedure and defined growth conditions for long-term propagation of mammary TEC. The isolated TEC maintain their endothelial specification and phenotype in culture. Furthermore, gene expression profiling of multiple TEC subpopulations revealed striking, persistent overexpression of several candidate genes including Irx2 and Zfp503 (transcription factors), Alcam and Cd133 (cell surface markers), Ccl4 and neurotensin (Nts) (angiocrine factors), and Gpr182 and Cnr2 (G protein-coupled receptors). Taken together, we have developed an effective method for isolating and culture-expanding mammary TEC, and uncovered several new TEC-selective genes whose overexpression persists even after long-term in vitro culture. These results suggest that the tumor microenvironment may induce changes in vascular endothelium in vivo that are stably transmittable in vitro.

Primary breast tumor-derived cellular models: characterization of tumorigenic, metastatic, and cancer-associated fibroblasts in dissociated tumor (DT) cultures.

Our goal was to establish primary cultures from dissociation of breast tumors in order to provide cellular models that may better recapitulate breast cancer pathogenesis and the metastatic process. Here, we report the characterization of six cellular models derived from the dissociation of primary breast tumor specimens, referred to as "dissociated tumor (DT) cells." In vitro, DT cells were characterized by proliferation assays, colony formation assays, protein, and gene expression profiling, including PAM50 predictor analysis. In vivo, tumorigenic and metastatic potential of DT cultures was assessed in NOD/SCID and NSG mice. These cellular models differ from recently developed patient-derived xenograft models in that they can be used for both in vitro and in vivo studies. PAM50 predictor analysis showed DT cultures similar to their paired primary tumor and as belonging to the basal and Her2-enriched subtypes. In vivo, three DT cultures are tumorigenic in NOD/SCID and NSG mice, and one of these is metastatic to lymph nodes and lung after orthotopic inoculation into the mammary fat pad, without excision of the primary tumor. Three DT cultures comprised of cancer-associated fibroblasts (CAFs) were isolated from luminal A, Her2-enriched, and basal primary tumors. Among the DT cells are those that are tumorigenic and metastatic in immunosuppressed mice, offering novel cellular models of ER-negative breast cancer subtypes. A group of CAFs provide tumor subtype-specific components of the tumor microenvironment (TME). Altogether, these DT cultures provide closer-to-primary cellular models for the study of breast cancer pathogenesis, metastasis, and TME.

Machine and deep learning methods for predicting 3D genome organization.

Three-Dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, Topologically Associating Domains (TADs), and A/B compartments play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers, Transcription Factor Binding Site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, TAD boundaries) and analyze their pros and cons. We also point out obstacles of computational prediction of 3D interactions and suggest future research directions.

The Human Intermediate Prolactin Receptor I-tail Contributes Breast Oncogenesis by Targeting Ras/MAPK Pathway.

Prolactin and its receptor (PRLr) in humans are significantly involved in breast cancer pathogenesis. The intermediate form of human PRLr (hPRLrI) is produced by alternative splicing and has a novel 13 amino acid tail ("I-tail") gain. hPRLrI induces significant proliferation and anchorage-independent growth of normal mammary epithelia in vitro when coexpressed with the long form hPRLr (hPRLrL). hPRLrL and hPRLrI coexpression is necessary to induce the transformation of mammary epithelia in vivo. The I-tail is associated with the ubiquitin-like protein neural precursor cell expressed developmentally downregulated protein 8. Treatment with the neural precursor cell expressed developmentally downregulated protein 8-activating enzyme inhibitor pevonedistat resulted in increased hPRLrL and the death of breast cancer cells. The goal of this study was to determine the function of the hPRLrI I-tail in hPRLrL/hPRLrI-mediated mammary transformation. hPRLrL/hPRLrI and hPRLrL/hPRLrIΔ13 (I-tail removal mutant) were delivered to MCF10AT cells. Cell proliferation was decreased when hPRLrI I-tail was removed. I-tail deletion decreased anchorage-independent growth and attenuated cell migration. The I-tail was involved in Ras/MAPK signaling but not PI3K/Akt signaling pathway as shown by western blot. I-tail removal resulted in decreased hPRLrI stability. RNA-sequencing data revealed that I-tail removal resulted in differential gene expression induced by prolactin. Ingenuity Pathway Analysis revealed that the activity of ERK was attenuated. Treatment of breast cancer cells with ERK1/2 inhibitor ulixertinib resulted in decreased colony-forming ability and less proliferation. These studies suggest that the hPRLrI I-tail contributed to breast oncogenesis and may be a promising target for the development of new breast cancer therapies.

The Transformative Role of 3D Culture Models in Triple-Negative Breast Cancer Research.

Advancements in cell culturing techniques have allowed the development of three-dimensional (3D) cell culture models sourced directly from patients' tissues and tumors, faithfully replicating the native tissue environment. These models provide a more clinically relevant platform for studying disease progression and treatment responses compared to traditional two-dimensional (2D) models. Patient-derived organoids (PDOs) and patient-derived xenograft organoids (PDXOs) emerge as innovative 3D cancer models capable of accurately mimicking the tumor's unique features, enhancing our understanding of tumor complexities, and predicting clinical outcomes. Triple-negative breast cancer (TNBC) poses significant clinical challenges due to its aggressive nature, propensity for early metastasis, and limited treatment options. TNBC PDOs and PDXOs have significantly contributed to the comprehension of TNBC, providing novel insights into its underlying mechanism and identifying potential therapeutic targets. This review explores the transformative role of various 3D cancer models in elucidating TNBC pathogenesis and guiding novel therapeutic strategies. It also provides an overview of diverse 3D cell culture models, derived from cell lines and tumors, highlighting their advantages and culturing challenges. Finally, it delves into live-cell imaging techniques, endpoint assays, and alternative cell culture media and methodologies, such as scaffold-free and scaffold-based systems, essential for advancing 3D cancer model research and development.

Estrogen switches pure mucinous breast cancer to invasive lobular carcinoma with mucinous features.

Mucinous breast cancer (MBC) is mainly a disease of postmenopausal women. Pure MBC is rare and augurs a good prognosis. In contrast, MBC mixed with other histological subtypes of invasive disease loses the more favorable prognosis. Because of the relative rarity of pure MBC, little is known about its cell and tumor biology and relationship to invasive disease of other subtypes. We have now developed a human breast cancer cell line called BCK4, in which we can control the behavior of MBC. BCK4 cells were derived from a patient whose poorly differentiated primary tumor was treated with chemotherapy, radiation and tamoxifen. Malignant cells from a recurrent pleural effusion were xenografted in mammary glands of a nude mouse. Cells from the solid tumor xenograft were propagated in culture to generate the BCK4 cell line. Multiple marker and chromosome analyses demonstrate that BCK4 cells are human, near diploid and luminal, expressing functional estrogen, androgen, and progesterone receptors. When xenografted back into immunocompromised cycling mice, BCK4 cells grow into small pure MBC. However, if mice are supplemented with continuous estradiol, tumors switch to invasive lobular carcinoma (ILC) with mucinous features (mixed MBC), and growth is markedly accelerated. Tamoxifen prevents the expansion of this more invasive component. The unexpected ability of estrogens to convert pure MBC into mixed MBC with ILC may explain the rarity of the pure disease in premenopausal women. These studies show that MBC can be derived from lobular precursors and that BCK4 cells are new, unique models to study the phenotypic plasticity, hormonal regulation, optimal therapeutic interventions, and metastatic patterns of MBC.

Characterization of cell lines derived from breast cancers and normal mammary tissues for the study of the intrinsic molecular subtypes.

Five molecular subtypes (luminal A, luminal B, HER2-enriched, basal-like, and claudin-low) with clinical implications exist in breast cancer. Here, we evaluated the molecular and phenotypic relationships of (1) a large in vitro panel of human breast cancer cell lines (BCCLs), human mammary fibroblasts (HMFs), and human mammary epithelial cells (HMECs); (2) in vivo breast tumors; (3) normal breast cell subpopulations; (4) human embryonic stem cells (hESCs); and (5) bone marrow-derived mesenchymal stem cells (hMSC). First, by integrating genomic data of 337 breast tumor samples with 93 cell lines we were able to identify all the intrinsic tumor subtypes in the cell lines, except for luminal A. Secondly, we observed that the cell lines recapitulate the differentiation hierarchy detected in the normal mammary gland, with claudin-low BCCLs and HMFs cells showing a stromal phenotype, HMECs showing a mammary stem cell/bipotent progenitor phenotype, basal-like cells showing a luminal progenitor phenotype, and luminal B cell lines showing a mature luminal phenotype. Thirdly, we identified basal-like and highly migratory claudin-low subpopulations of cells within a subset of triple-negative BCCLs (SUM149PT, HCC1143, and HCC38). Interestingly, both subpopulations within SUM149PT were enriched for tumor-initiating cells, but the basal-like subpopulation grew tumors faster than the claudin-low subpopulation. Finally, claudin-low BCCLs resembled the phenotype of hMSCs, whereas hESCs cells showed an epithelial phenotype without basal or luminal differentiation. The results presented here help to improve our understanding of the wide range of breast cancer cell line models through the appropriate pairing of cell lines with relevant in vivo tumor and normal cell counterparts.

Evaluation of the NRF1-proteasome axis as a therapeutic target in breast cancer.

Proteasomes are multi-subunit complexes that specialize in protein degradation. Cancer cells exhibit a heightened dependence on proteasome activity, presumably to support their enhanced proliferation and other cancer-related characteristics. Here, a systematic analysis of TCGA breast cancer datasets revealed that proteasome subunit transcript levels are elevated in all intrinsic subtypes (luminal, HER2-enriched, and basal-like/triple-negative) when compared to normal breast tissue. Although these observations suggest a pan-breast cancer utility for proteasome inhibitors, our further experiments with breast cancer cell lines and patient-derived xenografts (PDX) pointed to triple-negative breast cancer (TNBC) as the most sensitive subtype to proteasome inhibition. Finally, using TNBC cells, we extended our studies to in vivo xenograft experiments. Our previous work has firmly established a cytoprotective role for the transcription factor NRF1 via its ability to upregulate proteasome genes in response to proteasome inhibition. In further support of this notion, we show here that NRF1 depletion significantly reduced tumor burden in an MDA-MB-231 TNBC xenograft mouse model treated with carfilzomib. Taken together, our results point to TNBC as a particularly vulnerable breast cancer subtype to proteasome inhibition and provide a proof-of-principle for targeting NRF1 as a viable means to increase the efficacy of proteasome inhibitors in TNBC tumors.