Oncotarget

Research Papers:

NFkB disrupts tissue polarity in 3D by preventing integration of microenvironmental signals

PDF |  HTML  |  Supplementary Files  |  How to cite

Oncotarget. 2013; 4:2010-2020. https://doi.org/10.18632/oncotarget.1451

Metrics: PDF 3497 views  |   HTML 4103 views  |   ?  

Sabine Becker-Weimann, Gaofeng Xiong, Saori Furuta, Ju Han, Irene Kuhn, Uri-David Akavia, Dana Pe'er, Mina J. Bissell _ and Ren Xu

Abstract

Sabine Becker-Weimann1,*, Gaofeng Xiong2, Saori Furuta1, Ju Han1, Irene Kuhn1, Uri-David Akavia4, Dana Pe’er4, Mina J Bissell1, Ren Xu2,3,*

1 Life Sciences Division, Lawrence Berkeley National Laboratory1, Berkeley, CA

2 Markey Cancer Center, University of Kentucky, Lexington

3 Department of Molecular and Biomedical Pharmacology, University of Kentucky, Lexington

4 Department of Biological Sciences, Columbia University, New York, NY

* Contributed equally to this work

Correspondence:

Mina J Bissell, email:

Ren Xu, email:

Keywords: Disorganization gene signature, p65, RelB, three-dimensional tissue structure, tissue polarity.

Received: September 26, 2013 Accepted: October 12, 2013 Published: October 14, 2013

Abstract

The microenvironment of cells controls their phenotype, and thereby the architecture of the emerging multicellular structure or tissue. We have reported more than a dozen microenvironmental factors whose signaling must be integrated in order to effect an organized, functional tissue morphology. However, the factors that prevent integration of signaling pathways that merge form and function are still largely unknown. We have identified nuclear factor kappa B (NFkB) as a transcriptional regulator that disrupts important microenvironmental cues necessary for tissue organization. We compared the gene expression of organized and disorganized epithelial cells of the HMT-3522 breast cancer progression series: the non-malignant S1 cells that form polarized spheres (‘acini’), the malignant T4-2 cells that form large tumor-like clusters, and the ‘phenotypically reverted’ T4-2 cells that polarize as a result of correction of the microenvironmental signaling. We identified 180 genes that display an increased expression in disorganized compared to polarized structures. Network, GSEA and transcription factor binding site analyses suggested that NFkB is a common activator for the 180 genes. NFkB was found to be activated in disorganized breast cancer cells, and inhibition of microenvironmental signaling via EGFR, beta1 integrin, MMPs, or their downstream signals suppressed its activation. The postulated role of NFkB was experimentally verified: Blocking the NFkB pathway with a specific chemical inhibitor or shRNA induced polarization and inhibited invasion of breast cancer cells in 3D cultures. These results may explain why NFkB holds promise as a target for therapeutic intervention: Its inhibition can reverse the oncogenic signaling involved in breast cancer progression and integrate the essential microenvironmental control of tissue architecture.


Introduction

Cancer is an organ specific disease. It is widely accepted that the loss of growth regulation is the central mechanism shared by all cancers. We have argued that the tissue architecture is at least an equally important regulator of malignant transformation, and one that may explain organ specificity of cancer [1, 2]. Disruption of polarized tissue structure [3-5] and induction of cellular invasion are accompanied by extensive changes in the cellular microenvironment; in many cases the basement membrane is degraded, cell-extracellular matrix (ECM) adhesion is altered and growth-regulatory signals are lost as a result. This microenvironmental remodeling and the consequent changes in tissue structure contribute to the development of a malignant phenotype [5, 6]. Over two decades ago we developed a 3-dimensional (3D) culture system using laminin-rich ECM (lr-ECM) to distinguish normal and malignant breast phenotypes [7]. In 3D cultures, non-malignant HMT-3522 S1 (S1) cells form polarized, growth-arrested spheres, which closely resemble the alveolar structures of healthy breast tissues referred to as acini [8]. We have extended the series to characterize premalignant phenotype of additional cell-lines referred to as S3 series [9]. We have used these assays to identify and characterize the microenvironmental factors and signals controlling acquisition of a malignant phenotype and its reversal [7, 10-12] (for further references see [2]). The isogenic malignant counterpart of S1 cells, the HMT-3522 T4-2 (T4-2) cells, display an altered profile of microenvironmental receptors and downstream signaling, and form disorganized, proliferative structures displaying enhanced cell invasion under the same conditions [10, 12] (and references cited above).

In addition to differences in tissue organization, compared to their S1 counterparts, the malignant T4-2 cells have significantly higher expression of epidermal growth factor receptor (EGFR), phosphatidylinositol 3-kinase (PI3K), mitogen-activated protein kinase (MAPK) as well as beta1-integrin and matrix metalloproteinases (MMPs) such as MMP9 among others [5, 9, 10, 12]. We have repeatedly shown that by manipulating the microenvironmental signals in T4-2 and other malignant cells [12], the cells can be induced to form polarized, non-invasive structures resembling those formed by the non-malignant S1 cells (‘T4-2R’ cells). We have referred to this phenomenon as ‘phenotypic reversion’ [10].

Decreasing activation of any of the above pathways down-modulates their downstream signals, normalizes the signaling through all the other pathways and reverses the malignant phenotype of T4-2 cells in 3D culture. We know also that changes in glucose uptake play an equally important role in both disorganization and reversion [13]. Thus many seemingly disparate agents and mechanisms induce a common phenotypic change (reversion of the chaotic malignant phenotype to an organized acinar-like structure) if the cells are in a physiologically relevant 3D microenvironment.

This observation led us to hypothesize that depending on the microenvironmental signaling, the cells must experience changes in their transcriptional programs that accompany the observed changes in colony architecture. We also hypothesized that there would be one or more common transcriptional regulators that integrate the signals to effect this switch. To find such regulators, we used microarray data of non-malignant S1 and malignant T4-2 cells, as well as data that had been obtained for T4-2 cells ‘phenotypically reverted’ by the widest variety of inhibitors, and compared the common denominators. These analyses identified the transcription factor nuclear factor-kappa B (NFkB) as a prominent switch for regulating the change between disorganization and polarity in response to microenvironmental changes.

NFkB regulates a number of cancer-related processes including immune-response, inflammation, cell survival and cancer progression [14]. The NFkB family consists of five structurally and functionally conserved members in mammalian cells, including RELA (p65), RELB (p60), c-REL, NFkB1 (p105 and p50), and NFkB2 (p100 and p52) [15]. Elevated NFkB binding activity has been observed in both primary human breast cancer tissues and breast cancer cell lines, and contributes to the activation of CYCLIN D1, c-MYC, and MUC1 [16-18]. Despite considerable literature suggesting the importance of NFkB in cancer in general and in breast cancer in particular, the reasons behind its dominant activation in tumors are still undefined. We find that NFkB acts as a transcriptional rheostat the level of which directs a phenotype of order and its aberrant expression spells architectural disorder in response to microenvironmental changes. This finding could provide a fine-tuning of a framework for designing more compelling therapeutics.

Results

Comparison of malignant disorganized colonies with organized acini identified a set of genes common in reversion.

We showed previously that the modification of aberrant microenvironmental signaling in malignant human mammary epithelial cells (HMECs) reprograms them to form polarized, growth-arrested structures in 3D cultures [10, 12] (Figure 1A). To identify the alterations in gene expression profiles associated with such reprogramming events, we utilized microarrays of S1, T4-2, and T4-2R cells reverted by inhibitors of EGFR, beta1-integrin, MAPK, PI3K and by broadband inhibitors of MMPs as well as by reduced glucose metabolism (Figure 1B and Supplemental Table 1, GEO database, GSE50444). We obtained a gene expression matrix with data from 6560 genes and 24 microarrays – 5 of S1, 5 of T4-2 and 14 of T4-2, reverted with different agents (Supplemental Table 2). Unsupervised hierarchical clustering analysis showed that S1, T4-2 and T4-2R cells each cluster together (Supplemental Figure S1).

We identified a set of 345 genes that were differentially expressed in cells with different 3D phenotypes – polarized or disorganized - (LIMMA/R Bioconductor,p=0.005). The genes upregulated in T4-2 cells (180 genes) were called ‘disorganization gene signature’ (Supplemental Table 3). Genes that were upregulated in S1 cells and reverted T4-2 cells (165 genes) were referred to as ‘organization gene signature’. We used the disorganization gene signature for further analysis in order to find transcriptional regulators involved in the disturbance of acinar structure in T4-2 cells. The expression and hierarchical clustering of these genes is depicted in Fig. 1C.

Using the statistics gene ontology tool DAVID [19], we identified genes encoding proteins located in adherens junctions and the actin cytoskeleton as well as genes involved in processes such as cell adhesion, focal adhesion and actin cytoskeleton organization, that were significantly enriched in the disorganization gene signature (Figure 1D). We also found that genes involved in cell migration were overrepresented, consistent with a previous report correlating disruption of tissue polarity with endothelial cell migration [20]. Furthermore, gene set enrichment analysis (GSEA, [21]) of the gene expression matrix confirmed an increase in expression of actin organization in T4-2 cells (p=0.002).

A cluster of genes revealed by microarray analysis is associated with tissue polarity in HMT3522 cell lines.

Figure 1: A cluster of genes revealed by microarray analysis is associated with tissue polarity in HMT3522 cell lines. A) Immunofluorescence analysis of basal marker alpha6-integrin in S1, T4-2, and T4-2 cells phenotypically reverted in 3D culture. S1 and reverted T4-2 cells formed polarized spheroids, but T4-2 cell formed disorganized structures. Blue:DAPI, green: alpha6 integrin. Scale: 25µm B) Scheme showing the pathways that were blocked to reverse the malignant phenotypes of T4-2 cells in 3D culture for the microarray analysis. C) Unsupervised clustering of the 180 genes upregulated in T4-2 cells as compared to S1 cells and reverted T4-2 cells. D) Gene Ontology terms (biological processes, cell components, molecular functions) and KEGG pathways overrepresented in the disorganization gene signature.

NFkB associated with the expression of disorganization genes.

In order to identify common signaling pathways and/or transcription factors involved in the disorganization of the 3D structures, we employed several bioinformatics tools. Comparison of HMT-3522 expression profiles with identified gene sets (Gene Set Enrichment Analysis, [21]) revealed that genes related to the ‘positive regulation of IKK and NF-κB cascade’ are activated in disorganized T4-2 cells (p<0.09, Figure 2A). Using the network analysis tool Cytoscape, we identified a regulatory network associated with the differentially expressed genes. p65, i.e. one of the subunits of NFkB, was one of the major hubs of the network with more than 100 neighbors (Figure 2B). Most notably, transcription factor binding site analysis of the disorganization cluster using TransFind (transfind.sys-bio.net/, [22]) identified NFkB as a potential transcription factor for the disorganization genes with the third highest significance (p<0.005); 24 of the 180 genes were identified as potential NFkappaB target genes (Supplemental Table 3, marked in yellow). This corresponds to a percentage 1.8 times as high as in the total set of 6560 genes.

These results suggest that NFkB is a potential regulator of the disorganization genes, and that activation of NFkB may be associated with disruption of polarized mammary tissue structure.

GSEA and network analyses reveal NFkB as a potential regulator of disorganization genes.

Figure 2: GSEA and network analyses reveal NFkB as a potential regulator of disorganization genes. A) Gene set enrichment analysis [21] finds the gene set ‘positive regulation of IKK and NFkB cascade’ to be enriched among genes that show a high expression in unorganized cells (1:’negatively correlated’) and low expression in organized cells in the microarrays (0:’positively correlated’) (p=0.09). B) A regulatory network consisting of the 180 genes of the disorganization signature as well as connecting nodes that interact with at least two disorganization genes was identified using Cytoscape. Red nodes: disorganization genes, light blue nodes: connecting genes, green edges: interactions found in NCBI, dark blue node: p65.

Expression and activation of NFkB is associated with the disorganized tissue phenotype.

To confirm that NFkB is indeed activated in disorganized MECs compared to polarized cellular structures, we assessed the expression of the NFkB subunit p65 in S1, T4-2, and reverted T4-2 cells. We found that mRNA levels of p65 were significantly upregulated in disorganized T4-2 cells compared to polarized S1 cells or reverted T4-2 cells. Interestingly, p65 is one of the 180 genes of the disorganization signature in T4-2 cells (Supplemental Figure S2). Among the 25 genes with the highest correlation of expression with p65, there were three genes involved in ECM-receptor interaction, alpha6-integrin, beta4-integrin and thrombospondin (with correlation coefficients of 0.91, 0.87 and 0.86 respectively), indicating a tight connection between their expression and p65 (data not shown). p65 was also increased at the protein level in untreated T4-2 cells compared to S1 cells. Reversion of T4-2 cells by blocking of EGFR, beta1-integrin, MMP, and their downstream signals drastically reduced the expression of p65 in those cells (Figure 3A), suggesting that microenvironmental cues contribute to NFkB activation in the disorganized HMECs.

It has been reported that p65 is phosphorylated at Ser536 by IkappaB kinase (IKK), NFkB activating kinase, and Ribosomal S6 kinase 1 (RSK1), which in turn induces nuclear translocation of p65 [23]. Thus, we analyzed the phosphorylation at Ser536 and nuclear levels of p65. Quantification of immunoblotting results showed that Ser536 phosphorylation of p65 was increased in T4-2 cells and repressed in S1 and reverted T4-2 cells (Figure 3B). In addition, nuclear translocation of p65 was consistent with phosphorylation status (Figure 3C). In order to assess the NFkB activity in T4-2 and reverted T4-2 cells directly, we measured the luciferase activity in cells stably transfected with a luciferase reporter construct that is controlled by NFkB response elements. We found that the transactivation activities of NFkB were significantly reduced in T4-2 cells reverted by the EGFR inhibitor (Figure 3D). In summary, we find the malignant phenotype associated with higher levels of p65 mRNA, p65 protein levels as well as increased NFkappaB transcriptional activity.

The NFkB pathway is activated in disorganized T4-2 cells compared to polarized S1 and reverted T4 cells.

Figure 3: The NFkB pathway is activated in disorganized T4-2 cells compared to polarized S1 and reverted T4 cells. A) Immunoblotting analysis of p65 expression in S1, T4-2, and reverted T4-2 cells in 3D culture. Laminin A/C was used as a loading control. T4-2 cells were reverted by tyrphostin (Tyr), an EGFR inhibitor; LY294002 (LY), a PI3K inhibitor; GM6001 (GM), a MMP inhibitor; 2DG, a glucose metabolism inhibitor; VPA, a HDAC inhibitor; AIIB2, a β1 integrin blocking antibody, MAB225 (MAB), an EGFR blocking antibody. B) Immunoblotting analysis of phosphorylated p65 in S1, T4-2, and T4-2 cells reverted by tyrphostin. The blot results were quantified by AlphaEaseFC software, and expressed as relative levels of phosphorylated p65 to total p65. C) Immunoblotting of nuclear protein extracted from S1, T4-2, and T4-2 cells reverted with tyrphostin. D) Luciferase analysis of NFkB transcription activities in T4-2 and reverted T4-2 cells. Graph displays average ± SEM; n=4.

Downmodulation of NFkB leads to phenotypic reversion.

To show that inhibition of NFkB indeed leads to phenotypic reversion, we blocked NFkB activation in T4-2 cells by two means. T4-2 cells cultured on top of 3D laminin-rich ECM (lrECM) gels for four days were treated with Wedelolactone, an IKK inhibitor [24].The Wedelolactone-treated cells formed small spheroid structures with basal polarity (alpha6-integrin), whereas the vehicle-treated control cells remained as disorganized, larger structures (Figure 4A). We also utilized shRNA to reduce p65 expression in T4-2 cells (Figure 4B). Reducing p65 expression reprogrammed T4-2 cells to form polarized spheroids structures in 3D cultures (Figure 4A), and cell invasion was significantly suppressed (Figure 4C). Using quantitative RT-PCR, we measured several NFkB -targeted genes to compare the NFkB activity in control vs. p65 knockdown cells. Our results showed significant reduction in mRNA levels of SERPINA1, TNFSF12A, PLAT, DUSP, MMP19, and TNFSF9 after downregulation of p65 (Figure 4D). As expected from our previous data, we found that the protein levels of beta1-integrin and EGFR were also decreased in p65 knockdown cells (Figure 4B). NFkB provides a ‘disintegration node’ in a 3D tissue-like setting.

Inhibition of the NFkB pathway reverses the malignant phenotypes of T4-2 cells in 3D culture.

Figure 4: Inhibition of the NFkB pathway reverses the malignant phenotypes of T4-2 cells in 3D culture. A) Phase and immunofluorescence images of T4-2 cells. Cells were either treated with NFkB inhibitor (Wedelolactone, 10 µM) and cultured in 3D for four days orinfected with p65 shRNA lentivirus before plating in 3D culture. Blocking the NFkB pathway reprograms the cells to form polarized spheroids structure. Green: alpha6-integrin. Scale: 25µm. B) Immunoblotting analysis of beta1-integrin, EGFR and p65 expression in control and p65 shRNA-expressing T4-2 cells. C) Transwell invasion analysis of control and p65 shRNA-expressing T4-2 cells. Graph displays average ± SEM; n=3. D) Quantification of mRNA levels of NFkB target genes by real-time PCR. Control and p65 shRNA-expressing T4-2 cells were cultured in 3D for 4 days before RNA extraction.

The expression of the disorganization genes correlates also with the phenotype of other breast cancer cell lines in 3D.

The 3D morphologies of a substantial panel of human breast cancer cell lines correlate with their gene expression profiles in 3D cultures and tumor cell invasiveness [25]. Since mRNA levels of the 180 disorganization genes are associated with tissue structure in the HMT-3522 cell lines, we sought to determine whether these genes are similarly associated with 3D morphologies of other breast cancer cell lines. Comparing the average expression of the disorganization genes in 24 breast cancer cells lines, we found a higher expression in basal cell lines than in luminal cell lines (Supplemental figure S3). Unsupervised hierarchical clustering analysis of the disorganization gene set revealed a correlation of expression with subtype as well as 3D morphology. Basal B cell lines clustered together with higher expression of disorganization genes compared to luminal cell lines. Within the basal B cluster (I), the majority of the cancer cells with stellate morphology grouped together (Figure 5A). Most of the cell lines in the second cluster (II) are luminal and basal A subtypes with spheroid 3D, mass or grape-like phenotypes. The majority of basal B cells form stellate structure in 3D culture and lack E-cadherin expression, reflecting decreased cell–cell interactions and enhanced invasion [25]. Therefore, the up-regulation of disorganization genes seems associated with aggressive morphologies and increased cell invasion in 3D cultures.

Expression of disorganization genes and activation of the NF-&#x3ba;B pathway are associated with aggressive 3D phenotypes of multiple breast cancer cell lines.

Figure 5: Expression of disorganization genes and activation of the NF-κB pathway are associated with aggressive 3D phenotypes of multiple breast cancer cell lines. A) Unsupervised hierarchical clustering of breast cancer cell lines in 3D culture using the 180 disorganization genes. Each row represents a gene; each column represents a cell line. B) Bar graph of RelB mRNA levels in basal and luminal types of breast cancer cell lines. RelB mRNA expression is assessed by Affymetrix microarray and measured as log2 (probe intensities). C) Phase and immunofluorescence images of BT549 cells. Cells were either treated with Wedelolactone in 3D culture for 4 days or infected with RelB shRNA retrovirus before being seeded in 3D culture. D) Transwell invasion analysis of BT549 cells. The cells were treated with NF-κB inhibitor (Wedelolactone, 10 µM) or transfected with RelB shRNA vector.

To examine whether NFkB activity is functionally involved in the 3D phenotypes of other breast cancer cell lines, we chose BT549, a basal cancer cell line with a stellate morphology. It was treated with Wedelolactone in 3D culture. Most of the Wedelolactone-treated cells formed round or mass structures in 3D lrECM gels (Figure 5C). The colony structures were more organized in drug-treated cells as shown by beta-catenin and nuclear staining, although the metastatic cells were not fully polarized (Figure 5C), as was the case for MDA-MB-231 cells which require more than one inhibitory factor to be fully reverted [12]. The dramatic changes in tissue structures were accompanied by a significant reduction in cell invasion (Figure 5D). Since mRNA levels of the NFkB subunit RELB are elevated in basal-type breast cancer cells (Figure 5B), NFkB signals in those cells can also be blocked with shRNA targeting the RELB gene. We found that in 3D culture, RELB-silenced BT549 cells form structures with reduced invasive branches which we termed ‘mass’ (Figure 5C and supplemental Figure S4) [25] . In addition, knockdown of RELB significantly inhibited invasion of BT549 cells (Figure 5D). These results indicate the involvement of NFkB pathway in invasion of basal-type breast cancer cells in 3D gels.

A scheme showing the role of NF-&#x3ba;B as potential integrator of various microenvironmental cues and regulator of multicellular organization.

Figure 6: A scheme showing the role of NF-κB as potential integrator of various microenvironmental cues and regulator of multicellular organization.

Discussion

Loss of tissue organization is a trait of every malignant tumor, but growth rather than architecture has traditionally been accepted to be the most important hallmark of cancer. Yet for more than a decade we have shown repeatedly that the architecture of the breast acini in 3D cultures can be dominant over the malignant genome, i.e. the cancer cell phenotype and growth status can be reversibly altered if the architecture is restored via microenvironmental signaling [1, 2]. This reversibility in the absence of genetic changes [26] implies that the phenotype is linked to the regulatory status and relatively unaffected by the genetic composition of the cells. Interestingly, within the same malignant cells, the inhibition of many different signaling pathways leads to a similar organized phenotype [3] . Because of reversibility and commonality of phenotype, we postulated that the dramatic change seen with multiple reverting agents could indicate the existence of specific molecular ‘nodes’ that control architectural reorganization in response to this great variety of microenvironmental factors.

Our analyses revealed that the extensive common differences in the gene expression profile of cells with an organized and disorganized phenotype in mammary epithelial cells are associated with at least one transcriptional regulator, NFkB. This makes NFkB one of the critical nodes through which the microenvironmental signaling pathways pass/integrate. When NFkB is suppressed, the architecture is maintained and tissues remain functional. Should the microenvironmental signaling be disrupted, NFkB becomes activated, and both tissue polarity and normal functions become imbalanced. In our studies blocking NFkB with a chemical inhibitor or shRNA was sufficient for the reestablishment of basal polarity in malignant cells; conversely in all conditions where malignant breast cancer cells were ‘reverted’, NFkB levels were decreased.

NFkB has been shown to be involved in the development and progression of breast cancer: It is frequently highly activated in mammary tumors, in particular within disorganized ER-negative tumors where it has been reported to promote cancer progression by enhancing cancer cell survival and proliferation [27, 28]. To this list we now add an additional important function of NFkB as an important player in loss of tissue polarity and organization. It was shown previously that modulation of many different oncogenic pathways including EGFR, PI3K and integrin signaling leads to changes in NFkB activity [29-35]. Also, activation of the NFkB pathway is required for induction and maintenance of EMT in mammary tumors [36]. In our 3D model, down-modulation of all these pathways leads to a reorganized phenotype, and our analysis suggests that NFkB is a major molecular linker between these pathways and tissue architecture: Of the 24 genes that are upregulated in disorganized cell clusters and that have potential NFkB binding sites in their promotor region (Supplemental Table 3, marked in yellow), four genes were associated with actin cytsoskeletal organization, cell junction and cell adhesion in the GO database (ACTN1,EPB41L2, GNE, PANX1).

This has strong implications for NFkB in cancer prevention and as target for therapeutic efforts: If keeping NFkB balanced is the key to healthy mammary tissues, signaling gone awry in malignant cells can be compensated through reducing other stimuli activating NFkB, including reductions in stress and inflammation. As long as NFkB is high, integration of the aberrant signaling cannot be restored to normal levels. As a consequence, compensatory signaling pathways can be easily activated in cells that have become imbalanced. We suggest that inhibition of NFkB in breast cells could be one crucial change that would allow these cells to reorganize and thereby to re-integrate their signaling. The additional details of how exactly the microenvironmental signals reach NFkB and what are the other arms by which the switch occurs require further investigation.

Materials and Methods:

Cell culture.

HMT-3522 S1 and T4-2 cells were maintained on tissue culture plastic as previously described [7]. BT549 cells were propagated in DMEM/F12 (Sigma) with 10% fetal bovine serum (Invitrogen). Three dimensional laminin-rich extracellular matrix (3D lrECM) on-top cultures [37] were prepared by trypsinization of cells from tissue culture plastic, seeding of single cells on top of a thin gel of Engelbreth-Holm-Swarm (EHS) tumor extract (Matrigel: BD Biosciences; Cultrex BME: Trevigen), and addition of medium containing 5% EHS. S1 cells were seeded at a density of 3.1×104 cells per cm2; T4 and BT549 cell lines were seeded at 2.1 ×104 cells per cm2. S1 and T4-2 were maintained in their propagation medium with media change every 2 days. BT549 cells were maintained in H14 medium [7] with 1% fetal bovine serum.

Microarray hybridization and analysis.

S1, T4-2 and reverted T4 cells were isolated from 3D cultures with PBS/EDTA as previously described [37]. Purified total cellular RNA was extracted using RNeasy Mini Kit with on-column DNase digestion (Qiagen). RNA was quantified by measuring optical density at A260 and quality was verified by agarose gel electrophoresis. Affymetrix microarray analysis was performed using the Affymetrix HG-U133A High Throughput Array (HTA) GeneChip system, generating two sets of data from two separate array platforms (GPL3921 and GPL4685, GEO database GSE50444). Preprocessing, normalization and filtering was performed using R Bioconductor [38].MAS5 normalization was performed separately on the two sets of expression data. ProbeIDs were filtered out if the average log2 expression value was lower than 6.5. For genes with more than one probeID, we only considered the one with the highest variance. The two sets of expression values were row-centered (mean and variance) separately and then merged to one gene expression matrix. Unsupervised hierarchical clustering analysis was performed with Cluster (uncentered correlation, average linkage), and results were visualized with TreeView (http://www.eisenlab.org/eisen/?page_id=42). Gene set enrichment analysis was performed with GSEA v2.07 (http://www.broadinstitute.org/gsea/index.jsp, [21]). Network analysis was performed with Cytoscape 2.7.0. Transcription factor binding site analysis was performed using TransFind (-300 to 100 bp, top 500 genes, [22]). Potential targets of NFkB transcriptional regulation were retrieved by repeated transcription factor binding site analysis with subsets of the disorganization signature using Transfind. DAVID Gene Ontology was used to determine overrepresented GO and KEGG terms [19].

Virus production and transduction.

For p65 shRNA production, a double-stranded DNA oligonucleotide was generated from the following sequences: sense,5’-GATCCGGACATATGAGACCTTCAACTTCCTGTCAGATATATCTCTCCTTCCACACTTTTTG-3’; antisense,5’-AATTCAAAAAGGACATATGAGACCTTCAATCTGACAGGAAGTATATCTCTCCTTCCACACG-3’ (target sequence underlined; BamH1/EcoR1 cohesive ends italicized). Both oligonucleotides were annealed and ligated into BamH1/EcoR1 site of pGreen puro lentival vector (System Biosciences). Lentivirus production and transduction of target cells were conducted following guideline by System Biosciences. Briefly, lentivirus vector and packaging plasmid mix (System Biosciences) were transfected into 293FT cells (Invitrogen) using Lipofectamine® 2000. After 48 hours, medium was harvested, filtered and used to infect target T4-2 cells with the addition of polybrene (10 µg/ml). After 24 hours medium was replaced. At 72 hours post-infection puromycin (0.5 µg/ml) was added for selection and maintained throughout the culturing period. shRNA vector against RelB were from Open Biosystems. Retrovirus production and infection protocol from the company were followed.

Immunofluorescence and image acquisition.

Cells in lrECM gel were smeared on slides, dried briefly, and fixed with 4% paraformaldehyde and permeabilized with 0.5 % Triton X-100. Stained samples were imaged using a Spot RT camera attached to a Zeiss upright epifluorescence microscope or a Stanford Photonics XR/Mega-10 ICCD camera attached to a Solamere Technology Group (Salt Lake City, UT) spinning disk confocal system comprised of a Zeiss Axiovert 200M inverted microscope. Pictures were taken using a 63×oil immersion objective with QED InVivo imaging software at room temperature. The digital images were pseudocolored, overlayed and merged using ImageJ 1.38 or Adobe Photoshop 7.0.

RT-PCR and real-time PCR.

Total RNA was extracted from cells using RNeasy Mini Kit (Qiagen). cDNA was synthesized using Superscript first strand synthesis kit (Invitrogen) from 0.5-1.0 μg RNA samples. Quantitative real-time PCR analysis was performed with the Lightcycler System using the Lightcycler FastStart DNA Master SYBR Green I kit (Roche). The following Lightcycler PCR amplification protocol was used: 95°C for 10 min (initial denaturation), and 45 amplification cycles (95°C for 5 s, 60°C for 10 s, 72°C for 5 s). Amplification was followed by melting curve analysis to verify the presence of a single PCR product [39].

Luciferase assays and western blot analysis.

Three copies of NF-κB response element were linked to TK minimal promoter and cloned into a reporter vector pGL3 (Promega). T4-2 cells were co-transfected with pGL-NF-κB and pNeo plasmid (1:10). Stably transfected cells were isolated by G418 selection and cultured in 3D in the presence of Tyrphostin. Following treatment, equal amounts of cell lysates were assayed for luciferase activity.

S1 and T4-2 cells were cultured in 3D for 10 days. Total and nuclear protein was extracted from S1, T4-2, and reverted T4 cells or nuclei. Western blot experiments were performed as previously described [40]. After sonication, insoluble material was removed by centrifugation at 15,000 g for 10 min. Proteins (20 µg) from each sample were subjected to SDS gel electrophoresis and then transferred to nitrocellulose membrane (Schleicher & Schuell). The membrane was subsequently incubated in blocking buffer containing primary antibody and horseradish peroxidase-conjugated secondary antibodies. The membrane was then subjected to enhanced chemiluminescence (ECL) using the SuperSignal chemiluminescent substrate (Pierce, Rockford, IL).

Acknowledgements

This work was supported by grants from the U.S. Department of Energy, the Office of Biological and Environmental Research (contract no. DE-AC02-05CH1123) and a Distinguished Fellow Award from the above (contract no. 03-76SF00098), from the National Cancer Institute (awards R01CA064786, U54CA143836 and U54CA112970) to MJB and SB, the U.S. Department of Defense Medical and Materiel Command innovator awards (contract no.W81XWH0810736) to MJB, in part by a grant from the Breast Cancer Research Foundation to MJB, and in part by AHA (12SDG8600000) and ACS grants (IRG 85-001-22) to R.X. R.X. thanks Dr. Catherine Anthony for editorial assistance. We (MJB and SB) thank the Physical Sciences Oncology Center for providing support for completion of this work.

References

1. Hagios C, Lochter A and Bissell MJ. Tissue architecture: the ultimate regulator of epithelial function? Philosophical transactions of the Royal Society of London Series B, Biological sciences. 1998; 353(1370):857-870.

2. Bissell MJ and Hines WC. Why don’t we get more cancer? A proposed role of the microenvironment in restraining cancer progression. Nat Med. 2011; 17(3):320-329.

3. Bissell MJ, Kenny PA and Radisky DC. Microenvironmental regulators of tissue structure and function also regulate tumor induction and progression: the role of extracellular matrix and its degrading enzymes. Cold Spring Harbor symposia on quantitative biology. 2005; 70:343-356.

4. Aranda V, Haire T, Nolan ME, Calarco JP, Rosenberg AZ, Fawcett JP, Pawson T and Muthuswamy SK. Par6-aPKC uncouples ErbB2 induced disruption of polarized epithelial organization from proliferation control. Nat Cell Biol. 2006; 8(11):1235-1245.

5. Beliveau A, Mott JD, Lo A, Chen EI, Koller AA, Yaswen P, Muschler J and Bissell MJ. Raf-induced MMP9 disrupts tissue architecture of human breast cells in three-dimensional culture and is necessary for tumor growth in vivo. Genes Dev. 2010; 24(24):2800-2811.

6. Bissell MJ and Radisky D. Putting tumours in context. Nature reviews Cancer. 2001; 1(1):46-54.

7. Petersen OW, Ronnov-Jessen L, Howlett AR and Bissell MJ. Interaction with basement membrane serves to rapidly distinguish growth and differentiation pattern of normal and malignant human breast epithelial cells. Proc Natl Acad Sci U S A. 1992; 89(19):9064-9068.

8. Briand P, Nielsen KV, Madsen MW and Petersen OW. Trisomy 7p and malignant transformation of human breast epithelial cells following epidermal growth factor withdrawal. Cancer Res. 1996; 56(9):2039-2044.

9. Rizki A, Weaver VM, Lee SY, Rozenberg GI, Chin K, Myers CA, Bascom JL, Mott JD, Semeiks JR, Grate LR, Mian IS, Borowsky AD, Jensen RA, Idowu MO, Chen F, Chen DJ, et al. A human breast cell model of preinvasive to invasive transition. Cancer Res. 2008; 68(5):1378-1387.

10. Weaver VM, Petersen OW, Wang F, Larabell CA, Briand P, Damsky C and Bissell MJ. Reversion of the malignant phenotype of human breast cells in three-dimensional culture and in vivo by integrin blocking antibodies. J Cell Biol. 1997; 137(1):231-245.

11. Lelievre SA, Weaver VM, Nickerson JA, Larabell CA, Bhaumik A, Petersen OW and Bissell MJ. Tissue phenotype depends on reciprocal interactions between the extracellular matrix and the structural organization of the nucleus. Proc Natl Acad Sci U S A. 1998; 95(25):14711-14716.

12. Wang F, Weaver VM, Petersen OW, Larabell CA, Dedhar S, Briand P, Lupu R and Bissell MJ. Reciprocal interactions between beta1-integrin and epidermal growth factor receptor in three-dimensional basement membrane breast cultures: a different perspective in epithelial biology. Proc Natl Acad Sci U S A. 1998; 95(25):14821-14826.

13. Onodera Y, Nam J and Bissell MJ. Increased sugar uptake and metabolism lead to oncogenic activation via EPAC-Rap1 and O-GlcNAc pathways. J Clin Invest. 2013; in press.

14. Aggarwal BB. Nuclear factor-kappaB: the enemy within. Cancer Cell. 2004; 6(3):203-208.

15. Perkins ND. Integrating cell-signalling pathways with NF-kappaB and IKK function. Nat Rev Mol Cell Biol. 2007; 8(1):49-62.

16. Westerheide SD, Mayo MW, Anest V, Hanson JL and Baldwin AS, Jr. The putative oncoprotein Bcl-3 induces cyclin D1 to stimulate G(1) transition. Mol Cell Biol. 2001; 21(24):8428-8436.

17. Kim DW, Gazourian L, Quadri SA, Romieu-Mourez R, Sherr DH and Sonenshein GE. The RelA NF-kappaB subunit and the aryl hydrocarbon receptor (AhR) cooperate to transactivate the c-myc promoter in mammary cells. Oncogene. 2000; 19(48):5498-5506.

18. Lagow EL and Carson DD. Synergistic stimulation of MUC1 expression in normal breast epithelia and breast cancer cells by interferon-gamma and tumor necrosis factor-alpha. J Cell Biochem. 2002; 86(4):759-772.

19. Huang DW, Sherman BT and Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009; 4(1):44-57.

20. Chen A, Cuevas I, Kenny PA, Miyake H, Mace K, Ghajar C, Boudreau A, Bissell M and Boudreau N. Endothelial Cell Migration and Vascular Endothelial Growth Factor Expression Are the Result of Loss of Breast Tissue Polarity. Cancer Res. 2009; 69(16):6721-6729.

21. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES and Mesirov JP. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005; 102(43):15545-15550.

22. Kielbasa SM, Klein H, Roider HG, Vingron M and Bluthgen N. TransFind-predicting transcriptional regulators for gene sets. Nucleic Acids Res. 2010; 38:W275-W280.

23. Perkins ND. Post-translational modifications regulating the activity and function of the nuclear factor kappa B pathway. Oncogene. 2006; 25(51):6717-6730.

24. Kobori M, Yang Z, Gong D, Heissmeyer V, Zhu H, Jung YK, Gakidis MA, Rao A, Sekine T, Ikegami F, Yuan C and Yuan J. Wedelolactone suppresses LPS-induced caspase-11 expression by directly inhibiting the IKK complex. Cell Death Differ. 2004; 11(1):123-130.

25. Kenny PA, Lee GY, Myers CA, Neve RM, Semeiks JR, Spellman PT, Lorenz K, Lee EH, Barcellos-Hoff MH, Petersen OW, Gray JW and Bissell MJ. The morphologies of breast cancer cell lines in three-dimensional assays correlate with their profiles of gene expression. Mol Oncol. 2007; 1(1):84-96.

26. Rizki A, Weaver VM, Lee SY, Rozenberg GI, Chin K, Myers CA, Bascom JL, Semeiks JR, Grate LR, Man IS, Borowsky AD, Jensen RA, Idowu MO, Chen F, Chen DJ, Petersen OW, et al. A human breast cell model of preinvasive to invasive transition. Cancer Res. 2008; 68(5):1378-1387.

27. Biswas DK, Cruz AP, Gansberger E and Pardee AB. Epidermal growth factor-induced nuclear factor kappa B activation: A major pathway of cell-cycle progression in estrogen-receptor negative breast cancer cells. Proc Natl Acad Sci U S A. 2000; 97(15):8542-8547.

28. Connelly L, Barham W, Onishko HM, Sherrill T, Chodosh LA, Blackwell TS and Yull FE. Inhibition of NF-kappa B activity in mammary epithelium increases tumor latency and decreases tumor burden. Oncogene. 2010.

29. Bhat-Nakshatri P, Sweeney CJ and Nakshatri H. Identification of signal transduction pathways involved in constitutive NF-kappaB activation in breast cancer cells. Oncogene. 2002; 21(13):2066-2078.

30. Faraldo MM, Deugnier MA, Tlouzeau S, Thiery JP and Glukhova MA. Perturbation of beta1-integrin function in involuting mammary gland results in premature dedifferentiation of secretory epithelial cells. Mol Biol Cell. 2002; 13(10):3521-3531.

31. Fiorini E, Marchisio PC, Scupoli MT, Poffe O, Tagliabue E, Brentegani M, Colombatti M, Santini F, Tridente G and Ramarli D. Adhesion of immature and mature T cells induces in human thymic epithelial cells (TEC) activation of IL-6 gene trascription factors (NF-kappaB and NF-IL6) and IL-6 gene expression: role of alpha3beta1 and alpha6beta4 integrins. Dev Immunol. 2000; 7(2-4):195-208.

32. Friedland JC, Lakins JN, Kazanietz MG, Chernoff J, Boettiger D and Weaver VM. alpha 6 beta 4 integrin activates Rac-dependent p21-activated kinase 1 to drive NF-kappa B-dependent resistance to apoptosis in 3D mammary acini. J Cell Sci. 2007; 120(20):3700-3712.

33. Lamb LE, Zarif JC and Miranti CK. The androgen receptor induces integrin alpha6beta1 to promote prostate tumor cell survival via NF-kappaB and Bcl-xL Independently of PI3K signaling. Cancer Res. 2011; 71(7):2739-2749.

34. Weaver VM, Lelievre S, Lakins JN, Chrenek MA, Jones JC, Giancotti F, Werb Z and Bissell MJ. beta4 integrin-dependent formation of polarized three-dimensional architecture confers resistance to apoptosis in normal and malignant mammary epithelium. Cancer Cell. 2002; 2(3):205-216.

35. Zhou BP, Hu MC, Miller SA, Yu Z, Xia W, Lin SY and Hung MC. HER-2/neu blocks tumor necrosis factor-induced apoptosis via the Akt/NF-kappaB pathway. J Biol Chem. 2000; 275(11):8027-8031.

36. Huber MA, Azoitei N, Baumann B, Grunert S, Sommer A, Pehamberger H, Kraut N, Beug H and Wirth T. NF-kappaB is essential for epithelial-mesenchymal transition and metastasis in a model of breast cancer progression. J Clin Invest. 2004; 114(4):569-581.

37. Lee GY, Kenny PA, Lee EH and Bissell MJ. Three-dimensional culture models of normal and malignant breast epithelial cells. Nat Methods. 2007; 4(4):359-365.

38. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge YC, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004; 5(10).

39. Xu R, Spencer VA and Bissell MJ. Extracellular Matrix-regulated Gene Expression Requires Cooperation of SWI/SNF and Transcription Factors. J Biol Chem. 2007; 282(20):14992-14999.

40. Muschler J, Lochter A, Roskelley CD, Yurchenco P and Bissell MJ. Division of labor among the alpha6beta4 integrin, beta1 integrins, and an E3 laminin receptor to signal morphogenesis and beta-casein expression in mammary epithelial cells. Mol Biol Cell. 1999; 10(9):2817-2828.


Creative Commons License All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 4.0 License.
PII: 1451