Previous Article | Next Article ![]()
Molecular and Cellular Biology, January 2008, p. 687-704, Vol. 28, No. 2
0270-7306/08/$08.00+0 doi:10.1128/MCB.01617-06
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Division of Gastroenterology, Department of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania 19104,1 Abramson Cancer Center, University of Pennsylvania, Philadelphia, Pennsylvania 19104,2 Department of Genetics, University of Pennsylvania, Philadelphia, Pennsylvania 19104,3 Division of Endocrinology, Diabetes, and Metabolism, Department of Medicine, and Institute of Diabetes, Obesity and Metabolism, University of Pennsylvania, Philadelphia, Pennsylvania 19104,4 Cancer Research Program, Hospital for Sick Children, Toronto, Ontario M5G 1X8, Canada,5 Department of Laboratory Medicine and Pathobiology, University of Toronto, Toronto, Ontario M5G 1X8, Canada,6 Penn Bioinformatics Core, Center for Bioinformatics, University of Pennsylvania, Philadelphia, Pennsylvania 19104,7 Division of Hematology/Oncology, Department of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania 191048
Received 29 August 2006/ Returned for modification 26 September 2006/ Accepted 27 October 2007
|
|
|---|
), and there was a concomitant increase in lipogenic gene expression as a downstream effector of PPAR
. Importantly, Parvin-β suppressed breast cancer growth in vivo, with associated decreased proliferation. These data suggest that Parvin-β might influence breast cancer progression. |
|
|---|
Parvin-
, -β, and -
comprise a small family of widely expressed ILK-binding proteins with tandem calponin homology domains (30, 48, 51, 57, 70). The two best-characterized members, Parvin-
and -β, interact directly with the ILK kinase domain in a mutually exclusive manner (73) and modulate both its kinase activity and connections to the actin cytoskeleton (51), although the molecular mechanisms underlying these actions are only now beginning to emerge (43, 66, 73). Less is known about Parvin-
function.
Data from several studies suggest that Parvin-
and Parvin-β may have divergent actions in the regulation of ILK signaling and cytoskeletal dynamics. For example, Parvin-
was reported to facilitate ILK-mediated phosphorylation of AKT Ser473, with subsequent protection from apoptosis (13, 73), whereas Parvin-β overexpression in HeLa cells promoted apoptosis (73). Parvin-β also inhibited ILK kinase activity and reduced AKT Ser473 and glucogen synthase kinase 3β Ser9 phosphorylation in response to epidermal growth factor stimulation, as previously reported by us (43), consistent with negative regulation of ILK signaling. In contrast to Parvin-
, Parvin-β directly bound to the actin-binding protein
-actinin and was required for proper focal adhesion formation, lamellipodium maturation, and cell spreading (69, 70). In addition to its regulation of ILK, Parvin-β was also found to activate
PIX (ARHGEF6), a GTPase exchange factor for RAC and CDC42 (41, 55). Hence, Parvin-β is also implicated in RAC- and CDC42-mediated rearrangements of the actin cytoskeleton following adhesion to the ECM.
The PARVA gene is located on human chromosome 11p15, whereas PARVB and PARVG are juxtaposed on human chromosome 22q13.31 within an approximately 1-Mb region that undergoes frequent loss of heterozygosity in sporadic breast cancers (7, 27) and mismatch repair-proficient colorectal cancers (6, 27). Mutational analysis of PARVG in sporadic breast and colorectal tumors revealed several germ line polymorphisms but no evidence of somatic mutations (8). However, we demonstrated that Parvin-β mRNA and protein levels are reduced in primary breast tumors compared with adjacent normal breast tissue and also in certain breast cancer cell lines (43). Restoration of Parvin-β expression in metastatic MDA-MB-231 breast cancer cells resulted in reduced colony-forming ability in semisolid medium, increased adhesion to type I collagen, and impaired invasion through ECM, without affecting proliferation (43).
To identify the signaling pathways regulated by Parvin-β in breast cancer cells and gene expression alterations that may mediate the perturbed in vitro behavior of Parvin-β transfectants, we propagated control and Parvin-β-expressing populations on type I collagen-coated plastic (two dimensional [2D]), within a fibrillar type I collagen gel (three-dimensional [3D] collagen), or within basement membrane-containing Matrigel (3D Matrigel) and then subjected them to expression profiling. A compelling rationale for interrogation of breast cancer cells in 3D is that it enables elucidation of the impact of extracellular matrix composition on epithelial biology (24, 28, 54, 66) and evaluation of the influence of physical stiffness on underlying biological processes mediated by Parvin-β (53). We describe herein novel findings through the demonstration that Parvin-β reexpression in MDA-MB-231 cells: (i) increased levels of transcription factors associated with epithelial differentiation (inhibitor of DNA binding 2 [ID2] and Krüppel-like factor 4 [KLF4]) and (ii) increased cyclin-dependent kinase 9 (CDK9)-mediated Ser82 phosphorylation, and transcriptional activity of the nuclear hormone receptor peroxisome proliferator-activated receptor gamma 1 (PPAR
1). There was also a selected increase in expression of PPAR
target genes involved in lipid biosynthesis, lipid droplet formation, and cholesterol efflux in Parvin-β transfectants, particularly in 3D Matrigel culture. This may be due to significantly higher levels of the PPAR
coactivator 1
(PGC-1
) in 3D compared with 2D conditions. Finally, ectopic Parvin-β expression inhibited growth of MDA-MB-231 xenografts in vivo. These novel findings may start to explain how Parvin-β loss may influence breast cancer progression.
|
|
|---|
agonist rosiglitazone (Cayman Chemical Co., Ann Arbor, MI). U0126 (Sigma) and SP600125 (A.G. Scientific Inc., San Diego, CA) were reconstituted as recommended, and each was used at 10 µM. Gene expression profiling and bioinformatics analysis. MDA-MB-231 vector control (231_VC) and MDA-MB-231 Parvin-β (231_PARVβ) cells were propagated on type I collagen in monolayer (2D), in type I collagen gels (3D collagen), or in growth factor-reduced Matrigel gels (3D Matrigel) for 7 days. Cells were directly lysed in Trizol reagent (Invitrogen, Carlsbad, CA), and total RNA was isolated according to the manufacturer's instructions for expression profiling. Briefly, 50 ng of total RNA was converted to first-strand cDNA using Superscript II reverse transcriptase (Invitrogen) primed by a poly(T) oligomer that incorporated the T7 promoter. Second-strand cDNA synthesis was followed by in vitro transcription for linear amplification of each transcript. The resulting cRNA (200 ng) was used as template for randomly primed cDNA synthesis and a second round of in vitro transcription, which incorporated biotinylated CTP and UTP. The cRNA products were fragmented to 200 bp or less, heated at 99°C for 5 min, and hybridized for 16 h at 45°C to U133Plus 2.0 oligonucleotide microarrays (Affymetrix Inc., Santa Clara, CA). Microarrays were subsequently washed at low (6x SSPE [1x SSPE is 0.18 M NaCl, 10 mM NaH2PO4, and 1 mM EDTA {pH 7.7}]) and high (100 mM morpholineethanesulfonic acid [MES], 0.1 M NaCl) stringency and stained with streptavidin-phycoerythrin. The fluorescence signal was amplified by addition of biotinylated antistreptavidin and an additional aliquot of streptavidin-phycoerythrin stain. A confocal scanner was used to acquire the fluorescent signal after excitation (570 nm).
Affymetrix microarray suite 5.0 (MAS5) was used to quantitate mRNA expression levels; default values provided by Affymetrix were applied to all analysis parameters. The number of probe pairs meeting the default discrimination threshold (tau = 0.015) was used to assign a call of absent, present, or marginal for each assayed gene, and a P value was calculated to reflect confidence in the detection call. A weighted mean of probe fluorescence (corrected for nonspecific signal by subtracting the mismatch probe value) was calculated using the one-step Tukey biweight estimate. Global scaling was applied to allow comparison of gene signals across multiple microarrays. All signal values from one microarray were then multiplied by the appropriate scaling factor. The MAS5 algorithm was used to determine gene expression ratios (flagged as I [increased], D [decreased], MI [marginally increased], MD [marginally decreased], or NC [no change]) and P values for the change. Ratios of mRNA expression levels between 231_VC and 231_PARVβ cells were calculated by determining the inverse log2 of the signal log ratio for each differentially expressed gene. Array results from 2D and 3D experiments were compared, and Venn diagrams were generated using Genespring GX software (Agilent Technologies, Palo Alto, CA). Comparisons of the 3D microarray results to the myoepithelial and luminal epithelial gene expression signatures derived by Grigoriadis and colleagues (16) were performed with Genespring GX. Pathway analysis and functional grouping of genes were carried out using Ingenuity Systems (Redwood City, CA) software.
RT-PCR. Analysis of mRNA expression by reverse transcription-PCR (RT-PCR) was performed as previously described (47). For oligonucleotide sequences of the PCR primer pairs used, see the supplemental materials and methods. Normal human colon was obtained from the Cooperative Human Tissue Network at the hospital of the University of Pennsylvania, and total RNA was isolated from the purified epithelium as described previously (27).
qPCR.
mRNA expression levels were determined by SYBR green (Applied Biosystems, Warrington, United Kingdom) quantitative real-time RT-PCRs (qPCRs) with 231_VC and 231_PARVβ cells as described previously (43). Optimal PCR conditions were determined by performing primer matrix reactions and generating standard curves for each primer pair used. Assays for both the gene of interest and the β-actin internal control were performed in triplicate in an ABI PRISM 7000 sequence detection system, and results were analyzed as previously described (43). Two separate RNA preparations from each growth condition were evaluated and data combined. Thus, means and standard deviations were calculated from a total of six independent PCRs on cDNA synthesized from two separate total RNA preparations. For the primer pairs used see the supplementary materials and methods. TaqMan gene expression assays were used to determine expression levels of total PPAR
, PPAR
transcript variant 2 (PPAR
2), PGC-1
, and ABCA1 and were in accordance with the manufacturer's recommendations (Applied Biosystems, Foster City, CA). Student's t test for means was applied for evaluation of expression level differences using SPSS (version 12) statistical software (SPSS Inc., Chicago, IL). Frozen primary human breast tumors were obtained from the Cooperative Human Tissue Network at the hospital of the University of Pennsylvania, and cDNA was synthesized as previously described (27).
Western blot analysis.
Cells were isolated from Matrigel using cell recovery solution (BD Biosciences) according to the manufacturer's instructions. Whole-cell, nuclear, and cytoplasmic extracts were prepared as described previously (43). For analysis of PPAR
expression, 30 to 60 µg of whole-cell lysate (15 to 20 µg of nuclear extract) was used, whereas 15 to 25 µg of whole-cell lysate was used for analysis of other proteins. Polyvinylidene difluoride membranes were incubated with primary antibody overnight (4°C), followed by the appropriate horseradish peroxidase-conjugated secondary antibody (Amersham, GE Healthcare UK, Little Chalfont, United Kingdom) for 1 h at room temperature. Either β-actin (Sigma), GAPDH (glyceraldehyde-3-phosphate dehydrogenase; Ambion, Austin, TX), or lamin A/C (Cell Signaling Technologies) was used as a loading control. The primary antibodies used were as follows: PPAR
, mouse monoclonal antibody, clone E-8 (Santa Cruz Biotechnology, Santa Cruz, CA); phospho-PPAR
(pPPAR
; Ser82), rabbit polyclonal antibody (Upstate, Charlottesville, VA); green fluorescent protein, rabbit polyclonal antibody (Clontech); CDK9, rabbit polyclonal antibody (Santa Cruz Biotechnology); CDK7, rabbit polyclonal antibody (Santa Cruz Biotechnology); total mitogen-activated protein kinase (MAPK), rabbit polyclonal antibody (Upstate); and phospho-MAPK (pMAPK; T202, Y204), rabbit polyclonal antibody; phospho-Jun N-terminal protein kinase (pJNK; T183, Y185), rabbit polyclonal antibody; total JNK, rabbit polyclonal antibody; phospho-c-JUN (T183, Y185), rabbit polyclonal antibody, total c-JUN, rabbit polyclonal antibody; lamin A/C, rabbit polyclonal antibody; and hemagglutinin, mouse monoclonal antibody, clone 6E2 (all from Cell Signaling Technology Inc., Danvers, MA). Antibodies recognizing Parvin-
, Parvin-β, Parvin-
, KLF4, and PGC-1
are described elsewhere (35, 43, 71).
Dual luciferase assays.
Dual luciferase assays (Promega, Madison, WI) were conducted as previously described (26). Briefly, cells were seeded into 24-well plates, and each well was transiently transfected the following morning with the PPAR
response element (PPRE) firefly luciferase reporter vector (1 µg) and either 2 ng of the pRL-CMV Renilla control vector or 20 ng of the pRL-SV40 Renilla control vector using 2 µl/well of FuGene6 reagent (Roche). The respective firefly and Renilla luciferase activities of cell lysates (20 to 40 µl) were acquired sequentially using an Orion microplate luminometer (Berthold Detection Systems, Pforzheim, Germany).
Transient transfection of CDK9 and CDK7 siRNA. Cells were seeded into Vitrogen-coated six-well plates and transiently transfected at 50% to 70% confluence with Lipofectamine 2000 reagent (Invitrogen) mixed with prevalidated small interfering RNAs (siRNAs) to either human CDK9 or CDK7 or scrambled negative control siRNA no. 5 (all from Ambion) at a final concentration of 40 nM. Cells were propagated for 2 to 3 days before harvesting without a medium change.
Generation of firefly luciferase- and tdTomato-expressing MDA-MB-231 derivatives.
The pBabe-puro_luciferase expression vector (a generous gift from Wafik El-Deiry) was used to produce retroviral supernatants as described previously (14). Supernatant (400 µl) was used to transduce MDA-MB-231 derivatives by the spin infection method as previously described. Six separate populations of luciferase-transduced 231_VC and 231_PARVβ cells were then assayed for firefly luciferase activity in duplicate, and activity was normalized to total protein levels. One population from each cell line was selected for expansion in puromycin (5 µg/ml). tdTomato cDNA in pRSET-B (a generous gift from Roger Tsien) was PCR amplified and subcloned (BamHI/EcoRI) into the pBabe-BlaS expression vector to produce pBabe-BlaS_tdTomato (45, 58). Luciferase-transduced 231_VC and 231_PARVβ cells were transduced further with pBabe-BlaS_tdTomato as described above and selected in Blasticidin S (20 µg/ml; Sigma Chemical Co.). Percentages of tdTomato expression and fluorescence intensity levels were measured by fluorescence microscopy and fluorometry, respectively. Fluorescence intensity (excitation filter
, 530/25 nm; emission filter
, 580/35 nm; Synergy HT; Biotek Instruments Inc., Winooski, VT) was normalized to protein levels and expressed as mean intensity/µg protein ± standard deviation (n = 5).
Subcutaneous injections and in vivo optical imaging. MDA-MB-231 cells (1 x 106) were resuspended in serum-free DME medium containing 30% Matrigel and injected subcutaneously into the flanks (four injection sites per mouse) of anesthetized nonlethally irradiated immunodeficient mice (NCRNU-M; Taconic, Hudson, NY) using a 27.5-gauge needle. Mice were maintained with the approval of and according to the guidelines produced by the University of Pennsylvania Institutional Animal Care and Use and Committee and Small Animal Imaging Facility. In vivo bioluminescent imaging (BLI) was performed using an IVIS 100 series optical imaging system (Xenogen Corporation, Hopkinton, MA) as described previously (14). Briefly, following intraperitoneal injection of D-luciferin, three or four anesthetized mice were placed 25 cm from the camera, and luminescence was measured every 1 min (30-s exposure, medium bin) until the maximum luminescence intensity was reached (approximately 15 min postinjection). For quantitation, counts were converted to flux (photons of light/s) keeping the same-size field for each time point. A mixed-effects generalized linear regression model using a random effect for tumor (that is, accounting for the repeated measurements of each tumor) was used to compare 231_VC and 231_PARVβ growth rates over time (SAS version 9.1; SAS Institute, Cary, NC). A P value of <0.05 was taken to represent a significant difference in growth rates between the two cell lines. In vivo fluorescence images of tdTomato were acquired using Maestro instrumentation (Cambridge Research and Instrumentation, Woburn, MA) with DsRed filter settings. Anesthetized mice were scanned (550 to 700 nm) using a 5-nm step, and background fluorescence was eliminated by employing the spectral unmixing capability.
Immunohistochemistry. Tumors were fixed in 4% paraformaldehyde for 8 h. Sections (5 µm) were subjected to antigen retrieval in a citrate buffer. Immunohistochemistry was performed as previously described (25). The primary antibodies used were antipancytokeratin (Sigma Chemical Co.), antivimentin (Novus Biologicals, Littleton, CO), anti-Ki-67 (Vector Laboratories, Burlingame, CA), anti-cyclin D1 (Lab Vision, Fremont, CA), and anti-CD34 (BD Biosciences). Slides were imaged using a Nikon TE2000 Eclipse microscope with a QiCam (Q Imaging, Surrey, Canada) camera and IPLab imaging software (BD Bioscience Bioimaging, Rockville, MD). The Ki-67 proliferation index was obtained by choosing five random fields from each of the two 231_VC and 231_PARVb tumors stained (total n = 10) and calculating (following blinded counting) the average number of Ki-67-positive nuclei.
Microarray data accession number. Microarray data were deposited in the Gene Expression Omnibus database with accession no. GSE9747.
|
|
|---|
As expected, Parvin-β mRNA itself was higher in 231_PARVβ cells than control cells in each of the three growth conditions (Fig. 1A), and this was corroborated by Western blotting in 2D (Fig. 1B). Interestingly, endogenous Parvin-
mRNA levels were decreased in 231_PARVβ cells cultured in both 3D conditions (Fig. 1B), whereas Parvin-
protein levels were reduced by approximately 50% in both 2D and 3D Matrigel growth conditions (Fig. 1B). Parvin-
mRNA and protein levels were not changed. The number of genes with differential expression between the two cell lines cultured under each of the three conditions is summarized using Venn diagrams (Fig. 1C). From a global perspective, more genes were repressed in 231_PARVβ cells (5,535) than were activated (2,785) (see Tables S1 and S2 in the supplemental material). This provided the basis for interrogation of specific families of genes that might relate to Parvin-β mediated function(s).
![]() View larger version (25K): [in a new window] |
FIG. 1. Expression profiling of MDA-MB-231 stable transfectants cultured in 2D and 3D conditions (A) Expression ratios of Parvin- and Parvin-β mRNAs in 231_PARVβ cells versus 231_VC cells grown on type I collagen-coated plastic (2D), in a 3D fibrillar type I collagen gel (3D Col), or embedded in 40% Matrigel (3D Mat). (B) Western blot analysis of Parvin- , Parvin- , and Parvin-β expression in 231_VC and 231_PARVβ cells grown in 2D or 3D Mat conditions. Numbers indicate relative Parvin- expression following normalization to GAPDH (bottom panel). The lower Parvin-β band may result from translation initiation from a downstream in-frame ATG. Molecular weights (M.W.) are in thousands. (C) Venn diagrams depicting numbers of genes significantly increased (Inc.) or decreased (Dec.) in 231_PARVβ cells compared with 231_VC cells following gene expression profiling using Affymetrix U133_Plus2 gene chips (see Materials and Methods). The numbers of genes increased or decreased in each growth condition are shown. For a complete description of the genes altered, see Table S1 in the supplemental material. The total numbers of genes significantly upregulated and downregulated in 231_PARVβ cells were 2,785 and 5,535, respectively.
|
|
View this table: [in a new window] |
TABLE 1. Altered gene expression in 231_PARVβ cells following 2D or 3D culturea
|
![]() View larger version (37K): [in a new window] |
FIG. 2. Parvin-β induces transcription factors linked to epithelial differentiation. Levels of ID2 mRNA (A) and KLF4 mRNA (B) expression were evaluated by SYBR green qPCR analysis with 231_VC and 231_PARVβ cells cultured in 2D, 3D collagen (3D col), or 3D Matrigel (3D Mat). Graphs show the mean relative gene expression levels and standard deviations, and expression in 231_VC cells was set to 1 in each condition. N/S, not significant; *, P < 0.05. (C) Western blot analysis of KLF4 expression in 231_VC and 231_PARVβ cells. Both nuclear extracts (n.e.; 15 µg, 2D and 3D Mat conditions) and whole-cell lysates (w.c.l.; 80 µg, 3D Mat only) were examined (top). The membrane was then reprobed for lamin A/C (middle), and β-actin (bottom), as loading controls. Molecular weights (M.W.) are in thousands. (D and E) Comparison of gene expression alterations in 231_PARVβ cells cultured in 3D with gene expression signatures of purified human luminal and myoepithelial mammary epithelial cells (16) (D) Significantly more luminal epithelial transcripts (n = 44) were upregulated in 231_PARVβ 3D cultures than were myoepithelial transcripts (n = 30; P = 0.0392). (E) Significantly more myoepithelial transcripts were downregulated (n = 105) in 231_PARVβ 3D cultures than were luminal epithelial transcripts (n = 24; P < 0.0001).
|
mRNA, encoding a pleiotrophic nuclear hormone receptor implicated in breast cancer progression (9, 31, 34, 65), were increased in both 3D cultures (2.0-fold for 3D collagen and 3.5-fold for 3D Matrigel) of 231_PARVβ cells, and the increases were greater than that for 2D cultures (1.6-fold) (Fig. 3A; Table 1). The human PPARG gene contains three major promoters that together generate up to four discrete transcript variants, which may ultimately produce two protein isoforms (see Fig. S2 and S3 in the supplemental material), namely, epithelium-enriched PPAR
1 (synthesized from mRNAs transcribed from promoters 1 and 3) and adipocyte-enriched PPAR
2 (generated from promoter 2) (9, 34, 65). Thus, expression of each PPAR
transcript variant was investigated individually. Variant 2 was expressed at higher levels in 231_PARVβ cells cultured in both 2D and 3D Matrigel (Fig. 3B and C); however, its relative expression in MDA-MB-231 cells was significantly lower than that of variants 1 and 3 (data not shown), as expected (9). Whereas expression of PPAR
variant 1 appeared unchanged between control and 231_PARVβ cells, importantly, variant 3 mRNA was augmented in 231_PARVβ cells grown in 3D Matrigel (Fig. 3D). Of note, differential expression of variant 3 was not evident in 2D cultures (Fig. 3D).
![]() View larger version (50K): [in a new window] |
FIG. 3. Parvin-β induces PPAR mRNA expression. Levels of total PPAR mRNA (A) and PPAR 2 mRNA (B) expression in 231_VC and 231_PARVβ cells cultured in 2D, 3D collagen (3D col), or 3D Matrigel (3D Mat) were evaluated by TaqMan qPCR analysis. Graphs show the mean relative gene expression levels (normalized to β-actin mRNA) and standard deviations, and expression in 231_VC cells was set to 1 in each condition. N/S, not significant; *, P < 0.05. (C) Semiquantitative RT-PCR analysis of PPAR 2 mRNA expression (28 cycles) in 231_VC and 231_PARVβ cells grown in 2D or 3D Matrigel. β2-Microglobulin expression (23 cycles) was measured as an internal control. (D) Semiquantitative RT-PCR analysis (28 cycles) of PPAR variant 1 (top) and variant 3 (middle) levels in cells cultured in either 2D or 3D Matrigel conditions. RNA extracted from purified normal human colonic epithelial cells (N.C.E.) was used as a positive control for both transcripts. β2-Microglobulin expression (23 cycles) was measured as an internal control (bottom). (E) Relationship between Parvin-β and PPAR mRNA expression levels in primary breast tumors. Parvin-β and total PPAR mRNA levels in 19 primary breast ductal adenocarcinomas were determined by SYBR green and TaqMan qPCR, respectively, and normalized to β-actin mRNA levels. The linearized CT values for Parvin-β and PPAR are represented on the y axis. Dotted lines are used to link the Parvin-β and PPAR data points of the six specimens with the highest Parvin-β mRNA expression and the six specimens with the lowest. Solid lines are used for the seven remaining tumors. Significance of the correlation (P = 0.013) was calculated by simple linear regression (r = 0.56).
|
mRNA, we compared Parvin-β and PPAR
mRNA expression levels in sporadic primary breast tumors by qPCR. Analysis of 19 primary ductal adenocarcinomas revealed a statistically significant correlation (P = 0.013; correlation coefficient [r] = 0.56) between Parvin-β mRNA and total PPAR
mRNA levels (Fig. 3E).
Increased PPAR
Ser82 phosphorylation, transcriptional activity, and lipogenic gene expression in 231_PARVβ cells.
Since PPAR
mRNA was upregulated in 231_PARVβ cells, PPAR
protein expression was also assessed. Since posttranslational modification of PPAR
had been shown to regulate its function (1, 4, 5, 21), we also examined PPAR
phosphorylation and transcriptional activity. As expected, only the PPAR
1 isoform was detected in MDA-MB-231 cells (Fig. 4A), and levels of unphosphorylated PPAR
1 in 231_VC and 231_PARVβ cells did not appear appreciably different. However, the abundance of a slower-migrating pPPAR
1 species was increased in 231_PARVβ cells (22), regardless of the culture conditions (Fig. 4A and data not shown). Serine 82 of PPAR
1 (Ser112 in PPAR
2) was identified as the target residue for enhanced phosphorylation in 231_PARVβ cells through use of a pPPAR
(Ser82/112)-specific antibody (Fig. 4A, bottom). The transcriptional activity of PPAR
was evaluated using a PPRE-driven reporter gene assay. Basal PPAR
transcriptional activity was increased in 231_PARVβ cells compared to control cells (Fig. 4B). Incubation with the synthetic PPAR
agonist rosiglitazone increased reporter gene activity in both cell lines in a dose-dependent manner (33); however, rosiglitazone-stimulated reporter activity in 231_PARVβ cells was approximately 1.7-fold higher than in control cells at each concentration used (Fig. 4B). Increased PPAR
1 activity in rosiglitazone-stimulated 231_PARVβ cells was not due to modulation of pPPAR
1 abundance or phosphorylation by rosiglitazone (Fig. 4C).
![]() View larger version (52K): [in a new window] |
FIG. 4. Expression and activity of PPAR protein in 231_VC and 231_PARVβ cells. (A, top two panels) Western blot analysis of PPAR 1 and pPPAR 1 expression in 231_VC and 231_PARVβ cells cultured in 2D or 3D Matrigel (3D Mat; embedded in 40% growth factor-reduced Matrigel [40%] or layered on 80% Matrigel in medium containing 2% Matrigel [2%]) using anti-total PPAR (E-8) monoclonal antibody. MDA-MB-231 cells transduced with pMSCV encoding murine Ppar 2 or pMSCV alone served as a control to distinguish pPPAR 1 from PPAR 2. β-Actin was used as a loading control. Bottom two panels) A pPPAR (Ser82/112)-specific antibody was used in Western blot analysis of whole-cell lysates from 2D or 3D Matrigel cultures and shows higher pPPAR in 231_PARVβ cells than in 231_VC cells. β-Actin was used as a loading control. Molecular weights (M.W.) are in thousands. (B) Transcriptional activity of PPAR 1 in 231_VC (VC) and 231_PARVβ (PARVβ) cells cultured in 2D. Cells were cotransfected with the PPREx3 vector and the Renilla control vector as described in Materials and Methods. Cells were then treated with dimethyl sulfoxide alone (vehicle) or increasing concentrations of the selective PPAR agonist rosiglitazone (Rosi). The empty firefly luciferase vector (pGL3basic) was used as a negative control. The graph shows the mean normalized luciferase activities and standard errors from four independent transfections. The experiment was conducted three times. *, P < 0.05; **, P < 0.005; ***, P < 0.001. (C) Rosiglitazone treatment did not augment PPAR 1 phosphorylation in 231_VC or 231_PARVβ cells grown in either 2D or 3D Matrigel (data not shown) conditions. Cells were treated with rosiglitazone (1.0 µM) for 48 h, followed by Western blot analysis of whole-cell lysates. β-Actin was used as a loading control (bottom).
|
activates a transcriptional cascade promoting both lipogenesis and adipogenesis (9, 37, 40, 63). Since PPAR
-dependent transcriptional activity was higher in 231_PARVβ cells compared to control cells, as measured by a reporter gene assay (Fig. 4B), the array results were interrogated for evidence of differential expression of potential PPAR
targets. Importantly, five genes encoding enzymes with direct roles in the lipid biosynthetic pathway (citrate transporter gene [SLC25A1], ATP-citrate lyase gene [ACLY], fatty acid synthase gene [FASN], stearoyl-coenzyme A [CoA] desaturase gene [SCD], and enoyl-CoA hydratase 1 gene [ECH1]) were found to be coordinately upregulated in 231_PARVβ cells, but only when cultured in 3D Matrigel and to a lesser extent in 3D collagen (Table 1). Increased FASN and SCD mRNA levels in 231_PARVβ cells were confirmed by qPCR (Fig. 5A and B; Table 1). In addition, mRNA levels of adipocyte differentiation-related protein, vimentin, and caveolin 2, each of which plays a structural role in intracellular lipid droplet formation (20, 38, 42), were also upregulated after 3D culture in Matrigel (Fig. 5C; Table 1). Finally, PPAR
was shown to orchestrate activation of a cholesterol efflux pathway in macrophage foam cells through PPAR
-mediated upregulation of liver X receptor-
(LXR
), which in turn induced expression of the cholesterol transporter ABCA1 (10, 44). Both LXR
and ABCA1 mRNAs were increased in 231_PARVβ cells, and again only when cultured in the 3D conditions (Fig. 5D; Table 1). Taken together, these data indicate that only the 3D Matrigel culture system was permissive for consistent and robust activation of putative or known key PPAR
-regulated genes and cellular processes.
![]() View larger version (50K): [in a new window] |
FIG. 5. PPAR -regulated pathways involved in lipogenesis, lipid droplet formation, and cholesterol efflux are induced in 231_PARVβ cells cultured in 3D conditions only. Levels of fatty acid synthase (FASN) mRNA (A), stearoyl-CoA desaturase (SCD) mRNA (B), adipocyte differentiation-related protein (ADFP) mRNA (C), and cholesterol transporter (ABCA1) mRNA (D) in 231_VC and 231_PARVβ cells cultured in 2D, 3D collagen (3D col), and 3D Matrigel (3D Mat) were evaluated by TaqMan qPCR analysis. (E) Evaluation of PGC-1 mRNA expression in 231_VC and 231_PARVβ cells by TaqMan qPCR. Graphs show the mean relative gene expression levels (normalized to β-actin mRNA) and standard deviations, and expression in 231_VC cells was set to 1 in each condition. N/S, not significant; *, P < 0.05. "A" (absent) and "P" (present) refer to the call tags given (PGC-1 probe set 219195_at) by the Affymetrix MAS5.0 algorithm (see Materials and Methods). (F) Western blot analysis of PGC-1 protein expression in nuclear extracts (25 µg) of 231_VC and 231_PARVβ cells grown in 2D and 3D Mat conditions. Lysate from 293T cells engineered to overexpress PGC-1 was used as a positive control (35). Two different exposures are shown, and β-actin was used as a loading control (bottom). Molecular weights (M.W.) are in thousands.
|
transcriptional coactivators or corepressors between the 2D and 3D systems. Interrogation of the microarray data suggested that the PPAR
coactivator PGC-1
was induced when MDA-MB-231 cells were cultured in either 3D growth condition (data not shown). Indeed, qPCR analysis confirmed that PGC-1
mRNA levels were approximately 10- to 20-fold upregulated in both 231_VC and 231_PARVβ cells grown in either 3D system (type I collagen or Matrigel) compared to the 2D system (Fig. 5E). PGC-1
protein levels were also shown to be induced by 3D Matrigel culture (Fig. 5F).
Assessment of the role of CDK9 in Parvin-β-induced phosphorylation of PPAR
.
Considering that hyperphosphorylation of PPAR
1 was observed in 231_PARVβ cells and that MAPK and JNK had been demonstrated to phosphorylate PPAR
1/2 on Ser82/112 (1, 4, 5, 21), MAPK and JNK pathway activation was evaluated as natural candidates initially for mediating the phosphorylation. Both pMAPK and pJNK levels were suppressed in cells cultured in 3D compared with 2D culture, and no differences in pMAPK and pJNK levels between control and 231_PARVβ cells were found, suggesting that MAPK and/or JNK may not be responsible for the increased pPPAR
1 in 231_PARVβ cells (see Fig. S5 in the supplemental material). Furthermore, pharmacologically mediated inhibition of MAPK or JNK was unrevealing (see Fig. S4 in the supplemental material). These data indicated that a kinase(s) other than MAPK or JNK was responsible for phosphorylating PPAR
in the 231_PARVβ cells.
When bound to cyclin T partners, CDK9 comprises part of the basal machinery for RNA polymerase II-directed gene transcription (36). Mouse Cdk9 was shown recently to phosphorylate Ppar
2 on Ser112 during adipogenic differentiation of 3T3-L1 cells (22), a modification that was associated with enhanced Ppar
2 transcriptional activity. We evaluated initially whether a small-molecule CDK9 inhibitor, DRB, could abrogate PPAR
1 phosphorylation in breast cancer cells (39). A DRB dose-response experiment revealed that treatment with 25 µM DRB for 24 h attenuated PPAR
1 phosphorylation by approximately 50%, whereas higher concentrations (35 to 45 µM) appeared to completely inhibit phosphorylation in 231_PARVβ cells (Fig. 6A). The inhibitory effect of DRB (50 µM) on pPPAR
levels was observed by 16 h in control cells and by 9 h in Parvin-β-transfected cells, and pPPAR
levels continued to decline over time (Fig. 6B). Importantly, DRB treatment did not alter pMAPK or phospho-c-JUN levels in either cell line (Fig. 6B). We next perturbed CDK9 function using a genetic approach. Transfection with CDK9-specific siRNA was used to abrogate endogenous CDK9 function. CDK9 protein was reduced to negligible levels in both cell lines following transfection with CDK9-specific siRNA but not with a scrambled (control) siRNA (Fig. 6D). Diminished CDK9 expression abrogated PPAR
1 phosphorylation in 231_PARVβ cells, which was associated with a simultaneous increase in levels of the unphosphorylated protein (Fig. 6C), thereby establishing that Parvin-β mediates phosphorylation of PPAR
via CDK9. Since the data suggested that Parvin-β could influence CDK9 function, we hypothesized that Parvin-β protein may be present in the nucleus. Indeed, subcellular fractionation revealed that both CDK9 (Fig. 6C), and Parvin-β (Fig. 6E) were contained in the nuclei of MDA-MB-231 cells, previously not described for Parvin-β.
![]() View larger version (44K): [in a new window] |
FIG. 6. CDK9 mediates increased PPAR 1 Ser82 phosphorylation in 231_PARVβ cells. (A, top) Cells cultured in 2D were treated with the indicated concentrations of DRB for 24 h, and nuclear extracts were evaluated for total PPAR expression by Western blotting. pPPAR 1 levels were reduced by approximately 50% following treatment with 25 µM DRB. Lamin A/C was used as a loading control (bottom). (B) Cells were incubated in DRB (50 µM) for the indicated times followed by Western blot analysis for total PPAR expression. Reduced pPPAR 1 levels were observed at 16 h of treatment in 231_VC cells and at 9 h of treatment in 231_PARVβ cells (top). DRB treatment did not affect pMAPK or pJNK levels (middle). β-Actin was used as a loading control (bottom). (C) 231_VC and 231_PARVβ cells were transfected with siRNA specific to CDK9 or with scrambled (control) siRNA, followed by Western blotting of nuclear extracts for PPAR levels (top). Lamin A/C was used as an internal control (bottom). (D) Knockdown of CDK9 in 231_VC and 231_PARVβ cells transfected with CDK9-specific siRNA, scrambled siRNA, or vehicle alone (control) was confirmed by Western blotting for CDK9 expression (top). β-Actin was used as a loading control. (E) Nuclear extracts were prepared from 231_VC and 231_PARVβ cells cultured in 2D or 3D Matrigel (3D Mat) and analyzed for Parvin-β expression (top). Lamin A/C was used as an internal control (bottom). For panels A and C to E, molecular weights (M.W.) are in thousands.
|
![]() View larger version (44K): [in a new window] |
FIG. 7. 231_VC and 231_PARVβ cells were engineered to express both firefly luciferase and tdTomato. (A, i) Luciferase activity was measured in firefly luciferase-transduced cells in duplicate and normalized to protein levels. (ii) tdTomato fluorescence intensity was measured using a fluorometer and normalized to protein levels. Results are expressed as mean fluorescence intensities/µg protein ± standard deviations (n = 5). (B) Luciferase activities of subcutaneously implanted tumors (left, two 231_VC tumors; right, two 231_PARVβ tumors) from a representative mouse at four time points following injection. Images were acquired using the IVIS system (see Materials and Methods). The color bar to the right represents luminescence intensity in flux (light photons/s). Greyscale and color images were superimposed for visualization of tumor location. 231_PARVβ tumors demonstrate reduced growth rate compared with 231_VC tumors. (C) Qualitative correlation of bioluminescent and fluorescent signals from a representative mouse 12 days and 18 days following subcutaneous implantation, respectively (left, two 231_VC tumors; right, two 231_PARVβ tumors). One 231_PARVβ tumor (front right flank) shows marked growth retardation compared to its 231_VC control (front left flank).
|
![]() View larger version (22K): [in a new window] |
FIG. 8. Growth rates of 231_VC- and 231_PARVβ-derived tumors in vivo. The growth of tumors from three mice (231_VC, n = 6; 231_PARVβ, n = 6) was followed longitudinally for a period of 22 days using in vivo BLI. (A) Tumors derived from 231_VC cells increased in size over time (P < 0.0001). (B) Tumors derived from 231_PARVβ cells did not increase in size over time (P = 0.4372). (C) Statistical comparison of the curves from 231_VC and 231_PARVβ cells (see Materials and Methods) demonstrated that the rate of change in tumor size was significantly greater in 231_VC tumors than in 231_PARVβ tumors (P = 0.0003).
|
![]() View larger version (109K): [in a new window] |
FIG. 9. Analysis of tumors for histology and immunohistochemistry. Four tumors (two of 231_VC and two of 231_PARVβ) were excised from the mouse shown in Fig. 7B and analyzed by immunohistochemistry for a marker of epithelial intermediate filaments (pancytokeratin [CK]; A and B), a marker of mesenchymal intermediate filaments (vimentin; C and D), a marker of mouse vascular endothelial cells (CD34; E and F), and markers of proliferation (cyclin D1 [G and H] and Ki-67 [I and J]) and counterstained with hematoxylin. One representative image is shown for each marker (original magnification, x200). Scale bars, 50 µm. The number of Ki-67-positive nuclei was significantly lower for 231_PARVβ-derived tumors than for 231_VC-derived tumors (P = 0.03, Student's t test; data not shown).
|
|
|
|---|
Global gene expression profiling revealed major differences between Parvin-β-regulated transcriptional profiles of cells grown on type I collagen (2D) and those grown embedded in either type I collagen (3D collagen) or basement membrane proteins (3D Matrigel). It is well documented that 3D growth more closely mimics both the normal tissue microenvironment and tumor microenvironment than 2D culture (29, 53, 66), in terms of both chemical properties (allowing proper cell-ECM interactions) and physical properties (appropriate tissue stiffness). We wished to interrogate how differences between 2D and 3D systems, as well as ECM composition, affected Parvin-β-mediated functions. For example, entire classes of genes were deregulated in 231_PARVβ cells compared to control cells and were restricted to 3D growth conditions only (see Tables S1 and S2 in the supplemental material). Examples include metabolic genes and those encoding chromatin-modifying factors, DNA repair proteins, and replication factors (see Table S1 in the supplemental material). Additionally, coordinated regulation of multiple genes involved in a specific cellular process or signaling pathway was also observed, again particularly when cells were propagated in 3D. First, evidence for increased expression of genes involved in epithelial differentiation in the Parvin-β-transfected cells was obtained. ID2, encoding a basic helix-loop-helix transcriptional regulator (3), was upregulated in 231_PARVβ cells in embedded 3D culture (Fig. 2A). ID2 plays a critical role in differentiation of mammary epithelial cells in vitro (26, 30) and is also required for mammary epithelial cell proliferation and differentiation during pregnancy in vivo (46, 42). Likewise, KLF4, a zinc finger transcription factor expressed in many differentiated epithelial cell types, regulates cell cycle exit and induction of differentiation programs in both columnar and squamous epithelia (12, 15) and was specifically upregulated in 3D culture (Fig. 2B and C). Comparison of transcriptional profiles in a global fashion showed that Parvin-β was more likely to upregulate luminal epithelium-expressed genes than those found in myoepithelial cells (Fig. 2D) and was more likely to downregulate myoepithelium-enriched transcripts than luminal epithelial ones. Thus, it is conceivable that Parvin-β is involved in "early" differentiation of breast cancer cells in 3D but is not sufficient by itself to trigger "terminal" differentiation, as another gene(s) is needed to orchestrate expression of the full repertoire of differentiation-linked genes and a differentiated morphology.
Increased mRNA expression, Ser82 phosphorylation, and activity of PPAR
were observed in 231_PARVβ cells and serve as a key novel component of our findings. How exactly might Parvin-β regulate PPAR
? The first effect may be transcriptional in nature. PPAR
3 was specifically increased in 231_PARVβ cells cultured in 3D Matrigel. Alignment of the proximal sequence of PPARG promoter 3 from several species revealed conserved consensus
EF1/TCF8 (56), AP-1, C/EBP β/
(2), Sp1 (62), and RREB-1 (61) DNA binding sites (see Fig. S3 in the supplemental material). Indeed, the proximal and overlapping putative C/EBP β/
and Sp1 binding sites were almost totally conserved among humans, chimpanzees, rhesus monkeys, mice, and rats, as were the adjacent putative
EF1/TCF8 and AP-1 binding sites. Accordingly, qPCR analysis of C/EBP β mRNA levels revealed a twofold increase in 231_PARVβ cells cultured in 3D Matrigel only (data not shown). Similarly, mRNA levels of ras-responsive element binding protein 1 (RREB-1), a dual transcriptional activator/repressor (61), were found to be decreased in 231_PARVβ cells cultured in 3D Matrigel by array analysis (see Table S1 in the supplemental material). Therefore, C/EBP β and RREB-1 are candidate transcriptional regulators of PPARG promoter 3 in 3D cultures and require further evaluation. Surprisingly, despite an increase in mRNA levels, overall PPAR
protein levels were not higher in 231_PARVβ cells, suggesting that MDA-MB-231 cells exhibit tight regulation of PPAR
protein levels. Our data demonstrated that Parvin-β's effect(s) on PPAR
is posttranslational in nature. CDK9 was a particularly attractive candidate kinase, and pharmacologic and siRNA strategies were employed to perturb CDK9 function. DRB (25 µM) reduced pPARG1 levels by about 50% (Fig. 6A). Since the concentration of DRB required to inhibit CDK9 activity by 50% in vitro (IC50) was reported to be 3 µM (37) and the in vivo (using intact cells) IC50 was predicted to be in the vicinity of 30 µM, these results supported the notion of CDK9 involvement in regulation of PPAR
1 phosphorylation. Indeed, siRNA-mediated knockdown of CDK9 reduced pPPAR
1 levels in 231_PARVβ cells (Fig. 6C), thus demonstrating clearly a role for CDK9 in Parvin-β-mediated PPAR
Ser82 phosphorylation. Ectopically expressed Parvin-β was also resident in the nuclei of MDA-MB-231 cells, therefore suggesting potential regulation of CDK9 activity.
It was intriguing that potential PPAR
target genes such as those involved in lipid biosynthesis, lipid droplet formation, and cholesterol efflux were essentially activated only in 231_PARVβ cells maintained in 3D Matrigel culture, and to a lesser extent in 3D collagen (Table 1). PPAR
transcriptional activity is controlled not only by level of expression, availability of ligand, and phosphorylation but also by the availability of transcriptional cofactors (17, 31). Levels of PGC-1
, a coactivator for PPAR
, were substantially elevated in both 3D growth conditions (Fig. 5E and F), whereas expression of other PPAR
coactivators, such as NCOA1/SRC-1 (steroid receptor coactivator 1), and PGC-1β was unchanged (data not shown). In this model, induction of PGC-1
expression by 3D culture may then permit ligand-independent transactivation of PPAR
target genes. This induction may contribute to the apparently enhanced ability of Ser82-phosphorylated PPAR
in the Parvin-β transfectants to activate endogenous target genes when the cell-basement membrane interaction is maintained. We extended our 3D studies to in vivo xenotransplantation experiments with nude athymic mice. Indeed, Parvin-β transfectants exhibited suppressed tumor growth, as determined by two different optical imaging modalities, and a reduced proliferation index. Perhaps the suppression of proliferation in vivo is related to the induction of a specific early differentiation program observed in 3D cultures.
Could the induction of pPPAR
by Parvin-β play a role in mammary epithelial cell differentiation as a consequence? PPAR
1 is expressed highly in differentiated colonic epithelial cells and is implicated in the colonic epithelial differentiation program (15). More recently, the generation of conditional knockout mice revealed a requirement for Pparg in the normal differentiation of airway epithelial cells in the lung (59), where Ppar
1-dependent regulation of both differentiation-linked (for example, moesin) and developmentally regulated (for example, cathepsin B) gene expression was invoked. Prior work also implicated PPAR
1 in the differentiation of breast cancer cells since treatment of several breast cancer cell lines with PPAR
agonists stimulated lipid synthesis and also reduced proliferation and clonogenic growth in soft agar (47). However, no apparent phenotype was observed in the mammary tissues of mice lacking Ppar
expression in mammary epithelium (11), indicating that Ppar
-dependent effects may be subtle or that another factor can compensate in its absence. One unanswered question is whether loss of Pparg expression enhances or suppresses oncogene-induced breast carcinogenesis in vivo (52). A potential clue was obtained recently from the observation that in a large prospective study of patients with diabetes mellitus receiving the glucose-lowering PPAR
agonist pioglitazone (Actos) for at least 2.5 years, a significant reduction (P = 0.034) in the incidence of breast cancer was observed only in the group receiving pioglitazone (n = 3/2,605) compared with the blinded placebo control group (n = 11/2,633) (http://www.proactive-results.com/).
antibody and protein, and Mitch Lazar for the pMSCV_Ppar
2 construct and helpful reading of the manuscript. We thank members of the Rustgi and Lazar labs for helpful discussions. This work was supported in part by NIH/NIDDK grant R01-DK056645 (to A.K.R. and C.N.J.), a grant from the National Colorectal Cancer Research Alliance/EIF (to A.K.R.), a grant from the Irving A. Hansen Foundation (to A.K.R.), a Pennsylvania Department of Health Fellowship in Basic Cancer Research from the American Association for Cancer Research (to C.N.J.), a Concept Award (W81XWH-04-1-0658) from the U.S. Department of Defense Breast Cancer Research Program (to C.N.J.), and in part by grants from the National Cancer Institute of Canada (with funds from the Terry Fox Run) and the Canadian Institutes of Health Research (both to G.E.H.). G.E.H. was a CIHR scholar. This study was also supported by the Penn Center for Molecular Studies in Digestive and Liver Diseases and its core facilities (Morphology, Molecular Biology, and Cell Culture) through NIH/NIDDK grant P30 DK50306, the Penn Diabetes and Endocrinology Research Center (NIH/NIDDK grant DK19525), and the Optical/Bioluminescence Core of the Penn Small Animal Imaging Facility (supported in part by NIH grant CA105008).
Published ahead of print on 12 November 2007. ![]()
Supplemental material for this article may be found at http://mcb.asm.org/. ![]()
|
|
|---|
-crystallin enhancer-binding protein
EF1 is a repressor of E2-box-mediated gene activation. Mol. Cell. Biol. 14:5692-5700.This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»