Search

Single-nucleus analysis allows robust cell-type classification and helps to establish relationships between chromatin accessibility and cell-type-specific gene expression. Here, using samples from 92 women of several genetic ancestries, we developed a comprehensive chromatin accessibility and gene expression atlas of the breast tissue. Integrated analysis revealed ten distinct cell types, including three major epithelial subtypes (luminal hormone sensing, luminal adaptive secretory precursor (LASP) and basal-myoepithelial), two endothelial and adipocyte subtypes, fibroblasts, T cells, and macrophages. In addition to the known cell identity genes FOXA1(luminal hormone sensing), EHFand ELF5(LASP), TP63and KRT14(basal-myoepithelial), epithelial subtypes displayed several uncharacterized markers and inferred gene regulatory networks. By integrating breast epithelial cell gene expression signatures with spatial transcriptomics, we identified gene expression and signaling differences between lobular and ductal epithelial cells and age-associated changes in signaling networks. LASP cells and fibroblasts showed genetic ancestry-dependent variability. An estrogen receptor-positive subpopulation of LASP cells with alveolar progenitor cell state was enriched in women of Indigenous American ancestry. Fibroblasts from breast tissues of women of African and European ancestry clustered differently, with accompanying gene expression differences. Collectively, these data provide a vital resource for further exploring genetic ancestry-dependent variability in healthy breast biology.

Tumor complexity and intratumor heterogeneity contribute to subclonal diversity. Despite advances in next-generation sequencing (NGS) and bioinformatics, detecting rare mutations in primary tumors and metastases contributing to subclonal diversity is a challenge for precision genomics. Here, in order to identify rare mutations, we adapted a recently described epithelial reprograming assay for short-term propagation of epithelial cells from primary and metastatic tumors. Using this approach, we expanded minor clones and obtained epithelial cell-specific DNA/RNA for quantitative NGS analysis. Comparative Ampliseq Comprehensive Cancer Panel sequence analyses were performed on DNA from unprocessed breast tumor and tumor cells propagated from the same tumor. We identified previously uncharacterized mutations present only in the cultured tumor cells, a subset of which has been reported in brain metastatic but not primary breast tumors. In addition, whole-genome sequencing identified mutations enriched in liver metastases of various cancers, including Notch pathway mutations/chromosomal inversions in 5/5 liver metastases, irrespective of cancer types. Mutations/rearrangements in FHIT, involved in purine metabolism, were detected in 4/5 liver metastases, and the same four liver metastases shared mutations in 32 genes, including mutations of different HLA-DR family members affecting OX40 signaling pathway, which could impact the immune response to metastatic cells. Pathway analyses of all mutated genes in liver metastases showed aberrant tumor necrosis factor and transforming growth factor signaling in metastatic cells. Epigenetic regulators including KMT2C/MLL3 and ARID1B, which are mutated in >50% of hepatocellular carcinomas, were also mutated in liver metastases. Thus, irrespective of cancer types, organ-specific metastases may share common genomic aberrations. Since recent studies show independent evolution of primary tumors and metastases and in most cases mutation burden is higher in metastases than primary tumors, the method described here may allow early detection of subclonal somatic alterations associated with metastatic progression and potentially identify therapeutically actionable, metastasis-specific genomic aberrations.

ABSTRACTRNA Binding Proteins (RBPs) are a class of post-transcriptional regulatory molecules which are increasingly documented to be dysfunctional in cancer genomes. However, our current understanding of these alterations is limited. Here, we delineate the mutational landscape of ∼1300 RBPs in ∼6000 cancer genomes. Our analysis revealed that RBPs have an average of ∼3 mutations per Mb across 26 cancer types. We identified 281 RBPs to be enriched for mutations (GEMs) in at least one cancer type. GEM RBPs were found to undergo frequent frameshift and inframe deletions as well as missense, nonsense and silent mutations when compared to those that are not enriched for mutations. Functional analysis of these RBPs revealed the enrichment of pathways associated with apoptosis, splicing and translation. Using the OncodriveFM framework, we also identified more than 200 candidate driver RBPs that were found to accumulate functionally impactful mutations in at least one cancer. Expression levels of 15% of these driver RBPs exhibited significant difference, when transcriptome groups with and without deleterious mutations were compared. Functional interaction network of the driver RBPs revealed the enrichment of spliceosomal machinery, suggesting a plausible mechanism for tumorogenesis while network analysis of the protein interactions between RBPs unambiguously revealed the higher degree, betweenness and closeness centrality for driver RBPs compared to non-drivers. Analysis to reveal cancer-specific Ribonucleoprotein (RNP) mutational hotspots showed extensive rewiring even among common drivers between cancer types. Knockdown experiments on pan-cancer drivers such as SF3B1 and PRPF8 in breast cancer cell lines, revealed cancer subtype specific functions like selective stem cell features, indicating a plausible means for RBPs to mediate cancer-specific phenotypes. Hence, this study would form a foundation to uncover the contribution of the mutational spectrum of RBPs in dysregulating the post-transcriptional regulatory networks in different cancer types.

The Komen Tissue Bank is the only biorepository in the world for normal breast tissues from women. Below we report the acquisition and processing of breast tissue from volunteer donors and describe an experimental and analysis pipeline to generate a single-cell atlas. This atlas is based on single-cell RNA-seq and is useful to derive breast epithelial cell subcluster-specific gene expression signatures, which can be applied to breast cancer gene expression data to identify putative cell-of-origin.

CXCR4, a chemokine receptor, plays an important role in breast cancer growth, invasion, and metastasis. The transcriptional targets of CXCR4 signaling are not known. Microarray analysis of CXCR4-enriched and CXCR4-low subpopulations of the MDA-MB-231 breast cancer cell line, which has a constitutively active CXCR4 signaling network, revealed differential expression of ~ 200 genes in the CXCR4-enriched subpopulation. ITF2, upregulated in CXCR4-enriched cells, was investigated further. Expression array datasets of primary breast tumors revealed higher ITF2 expression in estrogen receptor negative tumors, which correlated with reduced progression free and overall survival and suggested its relevance in breast cancer progression. CXCL12, a CXCR4 ligand, increased ITF2 expression in MDA-MB-231 cells. ITF2 is a basic helix-loop-helix transcription factor that controls the epithelial-to-mesenchymal transition and the function of the ID family (inhibitor-of-differentiation) of transcription factors, such as ID2. ID2 promotes differentiation of breast epithelial cells and its reduced expression in breast cancer is associated with an unfavorable prognosis. Both CXCR4 and ITF2 repressed ID2 expression. In xenograft studies, CXCR4-enriched cells formed large tumors and exhibited significantly elevated lung metastasis. Short interfering RNA against ITF2 reduced invasion of the CXCR4-enriched MDA-MB-231 subpopulation, whereas ITF2 overexpression restored the invasive capacity of MDA-MB-231 cells expressing CXCR4shRNA. Furthermore, overexpression of ITF2 in these cells enhanced tumor growth. We propose that ITF2 is one of the CXCR4 targets, which is involved in CXCR4-dependent tumor growth and invasion of breast cancer cells.

Activated v-AKT murine thymoma viral oncogene homolog 1 (AKT)/protein kinase B (PKB) kinase (pAKT) is localized to the plasma membrane, cytoplasm, and/or nucleus in 50% of cancers. The clinical importance of pAKT localization and the mechanism(s) controlling this compartmentalization are unknown. In this study, we examined nuclear and cytoplasmic phospho-AKT (pAKT) expression by immunohistochemistry in a breast cancer tissue microarray (n= 377) with ≈15 years follow-up and integrated these data with the expression of estrogen receptor (ER)α, progesterone receptor (PR), and FOXA1. Nuclear localization of pAKT (nuclear-pAKT) was associated with long-term survival (P= 0.004). Within the ERα+/PR+ subgroup, patients with nuclear-pAKT positivity had better survival than nuclear-pAKT–negative patients (P≤ 0.05). The association of nuclear-pAKT with the ERα+/PR+ subgroup was validated in an independent cohort (n= 145). TCL1 family proteins regulate nuclear transport and/or activation of AKT. TCL1B is overexpressed in ERα-positive compared with ERα-negative breast cancers and in lung metastasis–free breast cancers. Therefore, we examined the possible control of TCL1 family member(s) expression by the estrogen:ERα pathway. Estradiol increased TCL1B expression and increased nuclear-pAKT levels in breast cancer cells; short-interfering RNA against TCL1B reduced nuclear-pAKT. Overexpression of nuclear-targeted AKT1 in MCF-7 cells increased cell proliferation without compromising sensitivity to the anti-estrogen, tamoxifen. These results suggest that subcellular localization of activated AKT plays a significant role in determining its function in breast cancer, which in part is dependent on TCL1B expression.

A subpopulation (CD44+/CD24-) of breast cancer cells has been reported to have stem/progenitor cell properties. The aim of this study was to investigate whether this subpopulation of cancer cells has the unique ability to invade, home, and proliferate at sites of metastasis. CD44 and CD24 expression was determined by flow cytometry. Northern blotting was used to determine the expression of proinvasive and 'bone and lung metastasis signature' genes. A matrigel invasion assay and intracardiac inoculation into nude mice were used to evaluate invasion, and homing and proliferation at sites of metastasis, respectively. Five among 13 breast cancer cell lines examined (MDA-MB-231, MDA-MB-436, Hs578T, SUM1315, and HBL-100) contained a higher percentage (>30%) of CD44+/CD24-cells. Cell lines with high CD44+/CD24-cell numbers express basal/mesenchymal or myoepithelial but not luminal markers. Expression levels of proinvasive genes (IL-1α, IL-6, IL-8, and urokinase plasminogen activator [UPA]) were higher in cell lines with a significant CD44+/CD24-population than in other cell lines. Among the CD44+/CD24--positive cell lines, MDA-MB-231 has the unique property of expressing a broad range of genes that favor bone and lung metastasis. Consistent with previous studies in nude mice, cell lines with CD44+/CD24-subpopulation were more invasive than other cell lines. However, only a subset of CD44+/CD24--positive cell lines was able to home and proliferate in lungs. Breast cancer cells with CD44+/CD24-subpopulation express higher levels of proinvasive genes and have highly invasive properties. However, this phenotype is not sufficient to predict capacity for pulmonary metastasis.

Macrophage inhibitory cytokine 1 (MIC-1), a divergent member of the transforming growth factor beta superfamily, plays a role in the progression of a number of cancers, including breast, gastric, prostate and colorectal carcinomas. Serum MIC-1 levels are elevated in patients with metastatic prostate, breast and colorectal carcinomas. In vitro studies have revealed a cell type-specific role for MIC-1 in senescence and apoptosis. MIC-1 activates the survival kinase AKT/PKB in neuronal cells. Depending on the cell type, it activates or represses the MAP kinases ERK1/2. Mechanisms responsible for an increased MIC-1 expression in cancers and the consequences of MIC-1 overexpression, however, are not known. In this study, we show that AKT/PKB directly regulates the expression of MIC-1 in breast cancer cells. Sequences within −88 to +30 of the MIC-1 promoter are required for the AKT-mediated induction of MIC-1. This region of the promoter contains two SP-1 binding sites (SP-1B and SP-1C), which bind to the SP-1 and SP-3 proteins. Mutation of SP-1C but not SP-1B reduced the AKT-mediated activation of MIC-1. MIC-1 increased the basal ERK1 phosphorylation and prolonged the estrogen-stimulated ERK1 phosphorylation in MCF-7 breast cancer cells without altering the phosphorylation status of AKT/PKB. Immunohistochemistry with MIC-1 antibody revealed an MIC-1 expression within the cancer cells of primary breast cancer and in the MCF-7 xenografts. Furthermore, a limited analysis of RNA from primary breast cancers revealed an overexpression of MIC-1 in tumors, compared with normal tissues. These results suggest that AKT/PKB through MIC-1 could regulate the ERK1 activity and the MIC-1 expression levels may serve as a surrogate marker for the AKT activation in tumors.

Estrogen receptors (ERs) mediate most of the biological effects of estrogen in mammary and uterine epithelial cells by binding to estrogen response elements in the promoter region of target genes or through protein-protein interactions. Anti-estrogens such as tamoxifen inhibit the growth of ER-positive breast cancers by reducing the expression of estrogen-regulated genes. However, anti-estrogen-resistant growth of ER-positive tumors remains a significant clinical problem. Here we show that phosphatidylinositol (PI) 3-kinase and AKT activate ERα in the absence of estrogen. Although PI 3-kinase increased the activity of both estrogen-independent activation function 1 (AF-1) and estrogen-dependent activation function 2 (AF-2) of ERα, AKT increased the activity of only AF-1. PTEN and a catalytically inactive AKT decreased PI 3-kinase-induced AF-1 activity, suggesting that PI 3-kinase utilizes AKT-dependent and AKT-independent pathways in activating ERα. The consensus AKT phosphorylation site Ser-167 of ERα is required for phosphorylation and activation by AKT. In addition, LY294002, a specific inhibitor of the PI 3-kinase/AKT pathway, reduced phosphorylation of ERα in vivo. Moreover, AKT overexpression led to up-regulation of estrogen-regulated pS2 gene, Bcl-2, and macrophage inhibitory cytokine 1. We demonstrate that AKT protects breast cancer cells from tamoxifen-induced apoptosis. Taken together, these results define a molecular link between activation of the PI 3-kinase/AKT survival pathways, hormone-independent activation of ERα, and inhibition of tamoxifen-induced apoptotic regression.

The orphan receptors COUP-TFI and COUP-TFII play an important role in development and differentiation by activating specific genes and by modulating the activity of nuclear receptors including estrogen receptor alpha (ERα) and retinoic acid receptors (RARs). Previously, it was demonstrated that the expression and activity of ERα and RARs are lost or impaired in anti-estrogen-resistant breast cancers. Here we show that, similar to ERα and RARs, the expression of COUP-TFII but not COUP-TFI is reduced in ∼30% of breast cancer cell lines. Introduction of COUP-TFII to MDA-MB-435 cells resulted in reduced growth and plating efficiency. Interestingly, COUP-TFII increased the expression of cyclin D1 and p21WAF1/CIP1 in MDA-MB-435 cells. Although parental and COUP-TFII-transduced cells progressed through the G1-S phase at a similar rate, progression of COUP-TFII cells through the G2/M transition phase was delayed. The activity of cdk2 required for G2/M progression was reduced in COUP-TFII cells compared to parental cells. This property of COUP-TFII is distinct from that of ERα and RARs, which usually modulate the G1 phase of breast cancer cells. Furthermore, these results reveal an important physiological function of COUP-TFII, which correlates with its ability to induce gene expression rather than modulation of nuclear receptor activity.

Cancer-induced skeletal muscle defects show sex-specific differences in severity with men performing poorly compared to women. Hormones and sex chromosomal differences are suggested to mediate these differences, but the functional skeletal muscle markers to document these differences are unknown. We show that the myogenic microRNA miR-486 is a marker of sex-specific differences in cancer-induced skeletal muscle defects. Cancer-induced loss of circulating miR-486 was more severe in men with bladder, lung, and pancreatic cancers compared to women with the same cancer types. In a syngeneic model of pancreatic cancer, circulating and skeletal muscle loss of miR-486 was more severe in male mice compared to female mice. Estradiol (E2) and the clinically used selective estrogen receptor modulator toremifene increased miR-486 in undifferentiated and differentiated myoblast cell line C2C12 and E2-inducible expression correlated with direct binding of estrogen receptor alpha (ERα) to the regulatory region of the miR-486gene. E2and toremifene reduced the actions of cytokines such as myostatin, transforming growth factor β, and tumor necrosis factor α, which mediate cancer-induced skeletal muscle wasting. E2- and toremifene-treated C2C12 myoblast/myotube cells contained elevated levels of active protein kinase B (AKT) with a corresponding decrease in the levels of its negative regulator PTEN, which is a target of miR-486. We propose an ERα:E2-miR-486-AKT signaling axis, which reduces the deleterious effects of cancer-induced cytokines/chemokines on skeletal muscle mass and/or function.

The transcription factor NF-kappaB regulates the expression of genes involved in cancer cell invasion, metastasis, angiogenesis, and resistance to chemotherapy. In normal cells NF-kappaB is maintained in the cytoplasm by protein-protein interaction with inhibitor IkappaBs. In contrast, in cancer cells a substantial amount of NF-kappaB is in the nucleus and constitutively activates target genes. To understand the mechanisms of constitutive NF-kappaB activation, we have analyzed the function of IkappaBalpha and IkappaBbeta in breast cancer cells. In most cases, constitutive NF-kappaB DNA binding correlated with reduced levels of either IkappaBalpha or IkappaBbeta isoforms. Overexpression of IkappaBalpha but not IkappaBbeta1 resulted in reduced constitutive DNA binding of NF-kappaB in MDA-MB-231 cells. Unexpectedly, IkappaBbeta1 overexpression moderately increased 12-O-tetradecanoylphorbol-13-acetate- and interleukin-1-inducible NF-kappaB DNA binding. 12-O-Tetradecanoylphorbol-13-acetate- and interleukin-1-induced transactivation by NF-kappaB, however, was lower in IkappaBbeta1-overexpressing cells. Mutants of IkappaBbeta1 lacking the C-terminal casein kinase II phosphorylation sites, which form a stable complex with DNA bound NF-kappaB without inhibiting its transactivation in other cell types, repressed the transactivation by NF-kappaB in MDA-MB-231 cells. Consistent with the results of transient transfections, the expression of urokinase plasminogen activator, an NF-kappaB target gene, was reduced in IkappaBbeta1-overexpressing cells. These results suggest that depending on the cell type, IkappaBbeta1 represses the expression of NF-kappaB-regulated genes by inhibiting either DNA binding or transactivation function of NF-kappaB.

NF-Y is a heterotrimeric transcription factor that specifically recognizes a CCAAT box motif found in a variety of eukaryotic promoter and enhancer elements. The subunit association and DNA binding properties of the NF-Y complex were examined as a function of redox state using recombinant NF-YA, NF-YB, and NF-YC subunits. Reduction of NF-YB by dithiothreitol (DTT) was essential for reconstitution of specific NF-Y CCAAT box DNA binding activity in vitro. Approximately 30% of the Escherichia coli-derived NF-YB subunit existed as intermolecular disulfide-linked dimers. NF-YB mutants in which the highly conserved cysteine residues at positions 85 and 89 had been converted to serines existed only as monomers and did not require DTT for functional NF-Y DNA binding activity. DTT was required, however, for the functional association of NF-YC with wild-type NF-YB but not with the NF-YB cysteine mutants. The cellular redox factors Ref-1 and adult T-cell leukemia-derived factor stimulated the DNA binding activity of recombinant NF-Y in the absence of DTT. Cells treated with 1-chloro-2,4-dinitrobenzene, an irreversible inhibitor of thioredoxin reductase, exhibited reduced endogenous NF-Y DNA binding activity. Together these results suggest that the cellular redox environment of mammalian cells is an important posttranscriptional regulator of NF-Y subunit association and DNA binding activities.

Breast cancers often progress from a hormone-dependent, nonmetastatic, antiestrogen-sensitive phenotype to a hormone-independent, antiestrogen- and chemotherapy-resistant phenotype with highly invasive and metastatic growth properties. This progression is usually accompanied by altered function of the estrogen receptor (ER) or outgrowth of ER-negative cancer cells. To understand the molecular mechanisms responsible for metastatic growth of ER-negative breast cancers, the activities of the transcription factor NF-κB (which modulates the expression of genes involved in cell proliferation, differentiation, apoptosis, and metastasis) were compared in ER-positive (MCF-7 and T47-D) and ER-negative (MDA-MB-231 and MDA-MB-435) human breast cancer cell lines. NF-κB, which is usually maintained in an inactive state by protein-protein interaction with inhibitor IκBs, was found to be constitutively active in ER-negative breast cancer cell lines. Constitutive DNA binding of NF-κB was also observed with extracts from ER-negative, poorly differentiated primary breast tumors. Progression of the rat mammary carcinoma cell line RM22-F5 from an ER-positive, nonmalignant phenotype (E phenotype) to an ER-negative, malignant phenotype (F phenotype) was also accompanied by constitutive activation of NF-κB. Analysis of individual subunits of NF-κB revealed that all ER-negative cell lines, including RM22-F5 cells of F phenotype, contain a unique 37-kDa protein which is antigenically related to the RelA subunit. Cell-type-specific differences in IκBα, -β, and -γ were also observed. In transient-transfection experiments, constitutive activity of an NF-κB-dependent promoter was observed in MDA-MB-231 and RM22-F5 cells of F phenotype, and this activity was efficiently repressed by cotransfected ER. Since ER inhibits the constitutive as well as inducible activation function of NF-κB in a dose-dependent manner, we propose that breast cancers that lack functional ER overexpress NF-κB-regulated genes. Furthermore, since recent data indicate that NF-κB protects cells from tumor necrosis factor alpha-, ionizing radiation-, and chemotherapeutic agent daunorubicin-mediated apoptosis, our results provide an explanation for chemotherapeutic resistance in ER-negative breast cancers.

Single-cell RNA sequencing (scRNA-seq) is an evolving technology used to elucidate the cellular architecture of adult organs. Previous scRNA-seq on breast tissue utilized reduction mammoplasty samples, which are often histologically abnormal. We report a rapid tissue collection/processing protocol to perform scRNA-seq of breast biopsies of healthy women and identify 23 breast epithelial cell clusters. Putative cell-of-origin signatures derived from these clusters are applied to analyze transcriptomes of ~3,000 breast cancers. Gene signatures derived from mature luminal cell clusters are enriched in ~68% of breast cancers, whereas a signature from a luminal progenitor cluster is enriched in ~20% of breast cancers. Overexpression of luminal progenitor cluster-derived signatures in HER2+, but not in other subtypes, is associated with unfavorable outcome. We identify TBX3 and PDK4 as genes co-expressed with estrogen receptor (ER) in the normal breasts, and their expression analyses in >550 breast cancers enable prognostically relevant subclassification of ER+ breast cancers.