Abstract
Activation of the serine/threonine kinase Akt contributes to the formation, maintenance, and therapeutic resistance of cancer, which is driving development of compounds that inhibit Akt. Phosphatidylinositol ether lipid analogues (PIA) are analogues of the products of phosphoinositide-3-kinase (PI3K) that inhibit Akt activation, translocation, and the proliferation of a broad spectrum of cancer cell types. To gain insight into the mechanism of PIAs, time-dependent transcriptional profiling of five active PIAs and the PI3K inhibitor LY294002 (LY) was conducted in non–small cell lung carcinoma cells using high-density oligonucleotide arrays. Gene ontology analysis revealed that genes involved in apoptosis, wounding response, and angiogenesis were upregulated by PIAs, whereas genes involved in DNA replication, repair, and mitosis were suppressed. Genes that exhibited early differential expression were partitioned into three groups; those induced by PIAs only (DUSP1, KLF6, CENTD2, BHLHB2, and PREX1), those commonly induced by PIAs and LY (TRIB1, KLF2, RHOB, and CDKN1A), and those commonly suppressed by PIAs and LY (IGFBP3, PCNA, PRIM1, MCM3, and HSPA1B). Increased expression of the tumor suppressors RHOB (RhoB), KLF6 (COPEB), and CDKN1A (p21Cip1/Waf1) was validated as an Akt-independent effect that contributed to PIA-induced cytotoxicity. Despite some overlap with LY, active PIAs have a distinct expression signature that contributes to their enhanced cytotoxicity. Mol Cancer Ther; 10(7); 1137–48. ©2011 AACR.
Introduction
The PI3K/Akt/mTOR pathway is a promising target in cancer, as its activation promotes cellular growth and survival and contributes to tumorigenesis in vivo, whereas inhibition of the pathway promotes apoptosis in cancer cells and increases responsiveness to chemotherapy or radiation (1–5). Akt has an important role in lung cancer, as it is activated in response to tobacco components in vitro (6, 7), and the phenotypic progression of tobacco carcinogen–induced lung lesions is dependent upon activation of Akt and mTOR (8, 9). In non–small cell lung carcinoma (NSCLC), Akt activation is specific for tumor tissues versus surrounding normal lung tissues and confers a poor prognosis (10). Despite the strong rationale to target Akt, however, few Akt inhibitors exist.
To address this need, we used molecular modeling to synthesize structurally modified phosphatidylinositol ether lipid analogues (PIA) designed to interfere with the pleckstrin homology (PH) domain of Akt (11). Five PIAs were identified that rapidly inhibited Akt activation, as well as the phosphorylation of multiple downstream substrates without affecting kinases upstream of Akt (12). PIAs selectively induced apoptosis in NSCLC and breast cancer cell lines with high endogenous levels of Akt activation. Although the PIAs appeared interchangeable in their abilities to inhibit Akt and cause cell death, they induced more cell death than an established phosphoinositide-3-kinase (PI3K) inhibitor, LY294002, despite similar inhibition of the Akt pathway, which suggests that PIAs might have additional targets. Support for this hypothesis came from studies of PIAs in the NCI60 cell line panel where activity of PIAs correlated with levels of phosphorylated but not total Akt, but other targets with higher correlation coefficients were identified (13).
Microarrays have been used to query transcriptional programs that underlie processes relevant to cancer such as proliferation (14), transformation (15), senescence (16), metastasis (17), epithelial to mesenchymal transition (18), and activation of oncogenic pathways (19). Elucidation of these programs is important to the development of new therapies. For example, transcriptional profiling of normal versus tumor tissues has led to the identification of new targets and pharmacodynamic biomarkers to predict efficacy and clinical outcome (20, 21). Profiling of drug-induced gene transcription has been used to uncover the mechanism of action of novel agents, elucidate structure activity relationships, and to determine on-target versus off-target effects (22, 23).
To elucidate transcriptional changes that could mediate the cytotoxic activity of PIAs, expression profiling in NSCLC cells was conducted. We segregated changes in gene expression that were shared by both PIAs and LY294002 and therefore likely due to effects on the PI3K/Akt pathway, from those that were unique to PIAs and therefore might be considered off-target effects. Although greatly overlapping with LY294002 in suppression of cell-cycle genes, active PIAs uniquely or potently induced a number of tumor suppressor genes that might contribute to biological properties of PIAs that extend beyond inhibition of Akt. This expression profile could underlie their enhanced toxicity and could be used in pharmacodynamic studies of PIAs.
Materials and Methods
Cell lines and materials
NSCLC cell lines (H157, A549, H1703, and H1155) were obtained from National Cancer Institute (NCI)/Navy Medical Oncology (Bethesda, MD). They were maintained in RPMI 1640 medium with 10% (v/v) FBS and incubated at 37°C in a 5.0% CO2 atmosphere. All lines were recently tested and authenticated by the Core Fragment Analysis Facility (Johns Hopkins University, Baltimore, MD) using a short tandem repeat profiling in accordance with AACR best practices. The synthesis of the PIAs has previously been described (11). LY294002 was purchased from Calbiochem. Antibodies to phospho-Akt (S473), Akt1, Akt2, Akt3, HSP70, and anti-mouse or anti-rabbit secondary antibodies were purchased from Cell Signaling Technologies. The DNA primase (p49) antibody was from Lab Vision Corporation. Antibodies to KLF6, MCM3, PCNA, and IGFBP3 as well as anti-goat secondary antibody were purchased from Santa Cruz Biotechnology, Inc. RhoB antibody was purchased from Proteintech Group, Inc. Protease inhibitor cocktail tablets were obtained from Roche Diagnostics GmbH and the Micro BCA Protein Assay Kit was from Pierce. The pcDNA3-HA-RhoB was a kind gift from Dr. George Prendergast. The pCMV6-KLF6 and pCMV6-CDKN1A were from OriGene. The pcDNA3-Myr-HA-Akt1 (24) was provided by Dr. William Sellers via Addgene. RhoB, KLF6 and CDKN1A On-Target plus human short interfering RNAs (siRNA) were from Dharmacon/Thermo. Protran pure nitrocellulose membranes were purchased from Schleicher & Schuell. All cell culture reagents were purchased from Life Technologies, Inc.
Pharmacologic treatment and oligonucleotide microarray analysis
NSCLC cells were plated 2 × 105 cells per well in 6-well plates or 2 × 106 in T-75 flasks in RPMI 1640 medium containing 10% FBS and incubated for 24 hours. The medium was then changed to RPMI 1640 medium with 0.1% FBS and the cells were incubated overnight. The following morning, cells were treated with 10 μmol/L PIA6 dissolved in dimethyl sulfoxide (DMSO) for 0, 2, 6, or 12 hours, and an equal volume of DMSO was added to control samples. For the PIA comparison, 10 μmol/L PIAs (5, 6, 7, 23–25) or 10 μmol/L LY294002 (LY) were incubated with the cells for 6 hours. PIA7 (an inert analogue consisting of only the lipid side chain) was used as a control. The structures of the PIAs and LY294002 are shown in Fig. 1A. Cell viability was not affected in 0.1% FBS for the duration of these experiments. In cells cultured with 5% FBS, PIAs are highly bound to serum proteins and higher concentrations are needed to observe the same effects. Following incubation, the alterations in cellular morphology were photographed, and cells from 6-well plates were harvested for immunoblot analysis. Total RNA was extracted from cells treated in T-75 flasks using TRIzol reagent (Invitrogen) and chloroform and purified according to the RNeasy Midiprep Spin Kit Protocol (Qiagen). Oligonucleotide microarray was conducted with dye swap. Microarray chips were generated from the 34,580mer long probe set Human Genome Oligo Set Version 3.0 (Qiagen). Protocols for cDNA labeling, hybridization, and scanning are available through the National Human Genome Research Institute microarray core. The raw data were deposited in a public functional genomics data repository Gene Expression Omnibus (GSE27911; ref. 25). Immunoblotting analysis was conducted as described previously (26).
Reverse transcriptase-PCR
Semiquantitative reverse transcriptase-PCR (RT-PCR) was carried out using the SuperScript One-Step RT-PCR System with Platinum Taq DNA polymerase (Invitrogen). The following forward and reverse primers were used: (i) DUSP1, 5′-ctgccttgatcaacgtctca-3′ and 5′-acccttcctccagcattctt-3′; (ii) KLF6, 5′-ggcaacagacctgcctagag-3′ and 5′-aggattcgctgctgacatct-3′; (iii) CENTD2, 5′-gctttgaggtcaacgagagg-3′ and 5′-gaagtagcgatcgtggaagc-3′; (iv) BHLHB2, 5′-ccttgaagcatgtgaaagca-3′ and 5′-gcttggccagatactgaagc-3′; (v) PREX1, 5′-ccctggtcagtgaagagagc-3′ and 5′-tcatctccagaccccatctc-3′; (vi) TRIB1, 5′-tctggctttgaggcttgttt-3′ and 5′-cagcccagagtccttagtcg-3′; (vii) KLF2, 5′-agagggtctccctcgatgac-3′ and 5′-tctcacaaggcatcacaagc-3′; (viii) RHOB, 5′-cgacgtcattctcatgtgct-3′ and 5′-cgaggtagtcgtaggcttgg-3′; (ix) CDKN1A, 5′-atgaaattcaccccctttcc-3′ and 5′-ccctaggctgtgctcacttc-3′; (x) C21orf58, 5′-cctcttccatcacggaggta-3′ and 5′-ggcacacaggtgtccctagt-3′; (xi) IGFBP3, 5′-cagagactcgagcacagcac-3′ and 5′-gatgaccggggtttaaaggt-3′; (xii) PCNA, 5′-ggcgtgaacctcaccagtat-3′ and 5′-tctcggcatatacgtgcaaa-3′; (xiii) PRIM1, 5′-gccatacgcatcattgacag-3′ and 5′-ccaccctttacaaggctcaa-3′; (xiv) MCM3, 5′-cgcaggaaaaacgagaagag-3′ and 5′-cagaccacacagctgaggaa-3′; (xv) HSPA1B, 5′-ccgagaaggacgagtttgag-3′ and 5′-gcagcaaagtccttgagtcc-3′; and (xvi) GAPDH, 5′-gagtcaacggatttggtcgt-3′ and 5′-ttgattttggagggatctcg-3′.
Bioinformatics tools for gene clustering, visualization, and ontology
The microarray outputs were clustered and visualized by Cluster 3.0 (27) and Java TreeView (28). Gene expression dynamics was analyzed by CAGED program (Cluster Analysis of Gene Expression Dynamics; ref. 29). For gene ontology (GO) analysis, the High-Throughput GoMiner web interface (30) was used as described (31).
Cell transfection and infection
Transfection of plasmid or siRNA was carried out with a Nucleofector device using program T-16 and transfection kit V (Lonza). Cells stably expressing Myr-Akt1 were created following plasmid transfection by G418 (800 μg/mL) selection for 2 weeks. Cell lines expressing Akt isoform–specific short hairpin RNA (shRNA) were created by lentiviral infection and shRNA vectors used were from Sigma-Aldrich unless otherwise noted: Akt1, NM_005163.1-628s1c1; Akt2, NM_001626.2-1509s1c1; Akt3, NM_005465.3-671s1c1; nontargeting, pLKO-scr (Addgene). Gene overexpression or knockdown was verified by immunoblotting.
MTS assay and fluorescence-activated cell sorting analysis
The MTS assay was conducted with CellTiter 96 Aqueous One Solution Reagent (Promega) according to the manufacturer's instructions, and the cell viability was determined by measuring the absorbance at 490 nm using a BioTek ELx800 Microplate Reader. Apoptosis (sub-G1 DNA) was quantified by propidium iodide (PI) staining and analyzed using a Becton Dickinson FACSort flow cytometer and CellQuest software.
Results
Optimization of PIA treatments and microarray analysis
Preliminary experiments were carried out to optimize conditions for microarray analysis. Previously, we observed that PIAs cause profound morphologic changes in NSCLC cells, including rounding and detachment. To assess the time dependence of these changes, H157 cells were treated with PIA6 and observed over time (Fig. 1B). At 2 hours, there was little morphologic change, but by 6 hours, the cells had become highly refractile and rounded. Between 6 and 12 hours, cellular detachment occurred. Similar time-dependent changes were observed with other active PIAs, but not an inactive PIA (PIA7) or LY (data not shown). In addition, PIA exposure caused similar morphologic changes in other NSCLC cell lines but with different kinetics. For example, these changes were delayed in A549 and H1703 cells but accelerated in H1155 cells (data not shown). In H157 cells treated with PIA6, the surviving fractions measured by MTS assay at 2, 6, and 12 hours were 95%, 79%, and 48%, respectively. These experiments suggest that at treatment times up to 6 hours, cellular detachment would not confound the measurement of gene expression changes induced by PIAs.
To assess Akt inhibition, immunoblotting was carried out with parallel samples prepared from H157 cells (Fig. 1C). PIA6 inhibited Akt phosphorylation at S473 at 2, 6, and 12 hours (left). Treatment with any of the 5 active PIAs or LY also decreased S473 phosphorylation in H157 cells at 6 hours (right). PIA7, an analogue that lacks the inositol ring, did not inhibit Akt phosphorylation. To assure that RNA quality and integrity were maintained with increasing times of exposure to PIAs, analysis using a Bioanalyzer Nanochip was conducted. The 28S and 18S rRNA bands were sharp up to 12 hours and the 28S bands were more intense than 18S bands, indicating the RNA quality was adequate (Fig. 1D, left). RNA integrity was also preserved in samples treated with all PIAs or LY for 6 hours (Fig. 1D, right). On the basis of the assessment of cellular morphology, Akt inhibition, and RNA quality, 6 hours was chosen as the time point at which to compare changes in gene expression with PIAs and LY.
Following microarray analysis, 911 genes were identified that exhibited differential expression by treatment with 1 or more of the 5 active PIAs in H157 cells (using a cutoff point of a 2-fold change in expression). A hierarchical clustered heat map and cluster tree of these data is presented in Fig. 1E. The expression of genes in response to LY is shown in the far right column. Clustering of expression signatures revealed that PIA23 and PIA25 showed the highest degree of similarity to each other. Although PIA5 showed similarity to PIA24, it altered fewer genes and many less potently than PIA24. PIA6 shared the most genes in common with the other 4 according to the clustering and was relatively close to the PIA23 and PIA25 subgroup.
To quantify genes that changed in common with PIAs and LY, Venn Mapper (32) was used to calculate the number of differentially expressed genes that overlapped between each PIA treatment and LY, as well as the corresponding z-scores (Supplementary Table S1). Of genes that increased, PIA23 and PIA25 each shared 33 genes with LY, with the highest positive z-scores of 5.8 and 7.2, respectively. However, PIA5, PIA6, and PIA24 had little overlap with LY, resulting in negative z-scores. Of genes that decreased, PIA5, PIA6, PIA23, PIA24, and PIA25 had 13, 11, 43, 25, and 35 genes in common with LY, respectively, all with positive z-scores (1.0, 4.5, 9.3, 3.3, and 8.2, respectively). Given that the cutoff z-score for statistical significance (a value of P < 0.05) was 1.96 or greater (http://www.gatcplatform.nl/vennmapper/index.php), these data indicate that PIA6, PIA23, PIA24, and PIA25 shared genes that decreased in common with LY in a statistically significant manner, but only PIA23 and PIA25 shared genes that increased in common with LY in a statistically significant manner (99.9% confidence level). Interestingly, the z-scores between any 2 PIAs were positive for increased or decreased comparisons but negative at increased-decreased or decreased-increased comparisons, suggesting these compounds affected gene expression in a similar direction but to different extents. The complete list of differentially expressed genes and enlarged heat map are shown in Supplementary Table S2 and Supplementary Fig. S1.
Early changes in gene expression caused by PIAs
From the clustered heat map, 83 transcripts were identified that were similarly regulated (either positively or negatively) by all 5 active PIAs. These were partitioned into 4 gene categories (C1–4) using k-means clustering (Fig. 2A). C1 and C2 were likely due to Akt pathway inhibition, as C1 includes genes that were induced by PIAs and LY and C2 includes genes that were repressed by both PIAs and LY. C3 and C4 include PIA-specific genes that are either repressed or induced by PIAs, respectively.
To identify the earliest changes in gene expression that occurred following PIA treatment, a time course experiment was carried out. RNA was isolated after 0, 2, 6, or 12 hours of exposure to PIA6. Sixty genes were extracted from 83 transcripts similarly regulated by PIAs at the 4 time points, and their expression levels with PIA6 exposure are depicted in Fig. 2B. Because temporal patterns of gene expression can be useful to identify common regulatory mechanisms, a temporal cluster analysis was conducted using the CAGED program. This program identified 15 clusters that exhibited distinct dynamic patterns over time, with 4 genes for each cluster (Fig. 2C and Supplementary Fig. S2). Because we wanted to identify the early and constitutively up- or downregulated genes that changed with time, clusters 2, 3, 6, 11, 14, and 15 were chosen for further analysis. The 24 genes from these 6 clusters were further partitioned on the basis of data from the PIA and LY comparisons, as well as the time course experiments (10 different conditions using CAGED), which generated 3 groups that exhibited similar patterns (Fig. 2D). The genes in groups I and II were upregulated by PIAs. Group I is characterized by genes not induced by LY. Group II is characterized by genes more strongly induced by PIA treatment but commonly induced by LY. In contrast, group III genes were downregulated by PIAs and in most cases, LY. The degree of downregulation caused by PIAs was similar to that of LY. The identity of genes comprising these groups, along with attributions of their function and process, are listed in Fig. 2E.
Validation of microarray data
Changes in gene expression depicted in groups I to III were validated using semiquantitative RT-PCR (Fig. 3A). Time dependence (left) and specificity for individual PIAs versus LY (right) were assessed. From group I, DUSP1, KLF6, CENTD2, BHLHB2, and PREX1 were selected for validation. KLF6 and DUSP1 were strongly induced by PIA6 from 2 to 12 hours, CENTD2 and BHLHB2 were slightly induced at 2 hours and strongly induced from 6 to 12 hours, and PREX1 was induced beginning at 6 hours. The mRNA levels of DUSP1, CENTD2, BHLHB2, and PREX1 were increased in the DMSO-treated sample at 12 hours, indicating that these genes may be sensitive to culture conditions such as nutrient consumption and/or increases in cell density. DUSP1, KLF6, CENTD2, BHLHB2, and PREX1 were all induced by the 5 active PIAs and not by LY or PIA7. These data confirm that these genes are rapidly induced by PIAs and are selective for active PIAs.
The 4 genes chosen from group II, TRIB1, KLF2, RHOB, and CDKN1A, were induced by PIA6 from 2 to 12 hours. In particular, RHOB and KLF2 were strongly induced at early time points. The basal level of KLF2 expression was hardly detectable, and its induction appeared transient with peak expression at 2 hours. Induction of TRIB1 and KLF2 was not selective for PIAs because LY had similar effects (Fig. 3A, middle right). RHOB and CDKN1A were also induced by the 5 active PIAs but were less potently by LY. Although individual variation was observed, expression of these genes from group II was similarly regulated by PIAs and LY, suggesting that these changes were likely due to inhibition of the PI3K/Akt pathway.
Six genes were chosen from group III for validation. Expression of IGFBP3, PCNA, and PRIM1 decreased from 2 to 12 hours. Decreased expression of C21orf58, MCM3, and HSPA1B was evident from 6 to 12 hours. When assessed for specificity between PIAs and LY, group III genes were repressed by both PIAs and LY, except for C21orf58, which was not inhibited by LY. Collectively, the expression of these genes as assessed by RT-PCR is in agreement with the microarray data.
To confirm that changes in mRNA expression would lead to protein-level alterations, immunoblotting was conducted (Fig. 3B). Because of limited availability of reliable antibodies, only a subset of the genes could be assessed. From group I, KLF6 protein expression increased after PIA treatment in 4 NSCLC cell lines. From group II, PIA treatment increased CDKN1A expression and markedly induced RhoB expression. When genes from group III were analyzed, the protein levels of PIA-repressed genes such as IGFBP3, PCNA, MCM3, and HSPA1B did not appreciably decrease at 6 hours. This could be related to slow protein turnover. To test this, we conducted a longer time course and found these proteins were decreased at 12 and 24 hours (Fig. 3C). These immunoblotting experiments validate protein expression as a readout for genes that are regulated by PIAs.
Identification of biologically relevant gene categories altered by PIAs
To explore the biological significance of the microarray data, GO analysis was conducted by uploading the total gene list represented on the chip, along with the induced or repressed genes, into the High-Throughput GoMiner web interface. A summary of this analysis from the 6-hour PIA comparison experiment is presented in Table 1. The “total genes” column indicates how many genes on the chip were assigned to the respective categories. The “changed genes” column indicates how many of the genes that changed significantly (>2-fold) during treatment were assigned to the targeted GO categories. The enrichment, P value, and false discovery rate (FDR) of the GO categories were calculated and listed with an FDR cutoff point of less than 0.05. Some general categories such as regulation of cellular, physiologic, or biological processes had the lowest P values, yet these communicate little meaning. The categories with the next most significant P values were more relevant, namely apoptosis and cell death, followed by response to wounding, cell cycle, and angiogenesis. Apoptosis is a known cellular outcome of PIA treatment (12). The “under-expressed” GO category result indicates genes associated with DNA-dependent DNA replication were inhibited by PIAs. The individual genes assigned to these GO categories are listed in Supplementary Table S3.
GO category . | Total genes . | Changed genes . | Enrichment . | log10 (P) . | FDR . |
---|---|---|---|---|---|
Overexpression | |||||
GO:0050789_regulation_of_biological_process | 2,170 | 24 | 2.1105 | −4.5343 | 0.0000 |
GO:0050791_regulation_of_physiologic_process | 1,998 | 22 | 2.1012 | −4.0073 | 0.0000 |
GO:0051244_regulation_of_cellular_physiologic_process | 1,873 | 21 | 2.1395 | −3.8992 | 0.0033 |
GO:0050794_regulation_of_cellular_process | 1,963 | 21 | 2.0414 | −3.5890 | 0.0050 |
GO:0006915_apoptosis | 356 | 8 | 4.2882 | −3.3683 | 0.0080 |
GO:0012501_programmed_cell_death | 358 | 8 | 4.2642 | −3.3519 | 0.0067 |
GO:0008219_cell_death | 377 | 8 | 4.0493 | −3.2021 | 0.0086 |
GO:0016265_death | 379 | 8 | 4.0280 | −3.1868 | 0.0075 |
GO:0009611_response_to_wounding | 324 | 7 | 4.1228 | −2.8881 | 0.0300 |
GO:0048519_negative_regulation_of_biological_process | 464 | 8 | 3.2901 | −2.6194 | 0.0256 |
GO:0051243_negative_regulation_of_cellular_physiologic_process | 396 | 7 | 3.3732 | −2.3934 | 0.0310 |
GO:0007050_cell_cycle_arrest | 42 | 3 | 13.6304 | −2.8766 | 0.0270 |
GO:0000074_regulation_of_cell_cycle | 327 | 7 | 4.0849 | −2.8648 | 0.0245 |
GO:0007049_cell_cycle | 504 | 8 | 3.0290 | −2.3968 | 0.0326 |
GO:0042330_taxis | 109 | 4 | 7.0028 | −2.6160 | 0.0228 |
GO:0006935_chemotaxis | 109 | 4 | 7.0028 | −2.6160 | 0.0228 |
GO:0048514_blood_vessel_morphogenesis | 46 | 3 | 12.4451 | −2.7616 | 0.0233 |
GO:0001525_angiogenesis | 46 | 3 | 12.4451 | −2.7616 | 0.0233 |
GO:0001568_blood_vessel_development | 46 | 3 | 12.4451 | −2.7616 | 0.0233 |
GO:0001944_vasculature_development | 46 | 3 | 12.4451 | −2.7616 | 0.0233 |
Underexpression | |||||
GO:0006261_DNA-dependent_DNA_replication | 58 | 3 | 26.3207 | −3.7506 | 0.0400 |
GO category . | Total genes . | Changed genes . | Enrichment . | log10 (P) . | FDR . |
---|---|---|---|---|---|
Overexpression | |||||
GO:0050789_regulation_of_biological_process | 2,170 | 24 | 2.1105 | −4.5343 | 0.0000 |
GO:0050791_regulation_of_physiologic_process | 1,998 | 22 | 2.1012 | −4.0073 | 0.0000 |
GO:0051244_regulation_of_cellular_physiologic_process | 1,873 | 21 | 2.1395 | −3.8992 | 0.0033 |
GO:0050794_regulation_of_cellular_process | 1,963 | 21 | 2.0414 | −3.5890 | 0.0050 |
GO:0006915_apoptosis | 356 | 8 | 4.2882 | −3.3683 | 0.0080 |
GO:0012501_programmed_cell_death | 358 | 8 | 4.2642 | −3.3519 | 0.0067 |
GO:0008219_cell_death | 377 | 8 | 4.0493 | −3.2021 | 0.0086 |
GO:0016265_death | 379 | 8 | 4.0280 | −3.1868 | 0.0075 |
GO:0009611_response_to_wounding | 324 | 7 | 4.1228 | −2.8881 | 0.0300 |
GO:0048519_negative_regulation_of_biological_process | 464 | 8 | 3.2901 | −2.6194 | 0.0256 |
GO:0051243_negative_regulation_of_cellular_physiologic_process | 396 | 7 | 3.3732 | −2.3934 | 0.0310 |
GO:0007050_cell_cycle_arrest | 42 | 3 | 13.6304 | −2.8766 | 0.0270 |
GO:0000074_regulation_of_cell_cycle | 327 | 7 | 4.0849 | −2.8648 | 0.0245 |
GO:0007049_cell_cycle | 504 | 8 | 3.0290 | −2.3968 | 0.0326 |
GO:0042330_taxis | 109 | 4 | 7.0028 | −2.6160 | 0.0228 |
GO:0006935_chemotaxis | 109 | 4 | 7.0028 | −2.6160 | 0.0228 |
GO:0048514_blood_vessel_morphogenesis | 46 | 3 | 12.4451 | −2.7616 | 0.0233 |
GO:0001525_angiogenesis | 46 | 3 | 12.4451 | −2.7616 | 0.0233 |
GO:0001568_blood_vessel_development | 46 | 3 | 12.4451 | −2.7616 | 0.0233 |
GO:0001944_vasculature_development | 46 | 3 | 12.4451 | −2.7616 | 0.0233 |
Underexpression | |||||
GO:0006261_DNA-dependent_DNA_replication | 58 | 3 | 26.3207 | −3.7506 | 0.0400 |
Abbreviations: GO, gene ontology; PIA, phosphatidylinositol ether lipid analogue; FDR, false discovery rate.
Time-dependent changes in GO were also assessed (Table 2). The categories of genes induced the earliest (2h) were death related. This process continued and was further defined as apoptosis from 6 to 12 hours. By 12 hours, several categories were added including protein and macromolecule biosynthesis, wound healing, and angiogenesis. The first categories of repressed genes included DNA-dependent DNA replication and cell-cycle regulation, which were only evident beginning at 6 hours. By 12 hours, several other categories were added, similar to the observations with induced genes. Collectively, the GO analysis suggests that although PIAs induce changes in expression of a fairly small number of genes, these changes become manifest by altering many cellular processes in ways that would likely be detrimental to the growth and survival of cancer cells.
. | Changed/total genes (FDR) . | ||
---|---|---|---|
GO category . | 2 h . | 6 h . | 12 h . |
Overexpression | |||
GO:0001944_vasculature_development | 10/46 (0.0227) | ||
GO:0001568_blood_vessel_development | 10/46 (0.0227) | ||
GO:0048514_blood_vessel_morphogenesis | 10/46 (0.0227) | ||
GO:0001525_angiogenesis | 10/46 (0.0227) | ||
GO:0050878_regulation_of_body_fluids | 16/98 (0.0238) | ||
GO:0007599 hemostasis | 15/87 (0.0209) | ||
GO:0042060_wound_healing | 15/87 (0.0209) | ||
GO:0009611_response_to_wounding | 41/324 (0.0200) | ||
GO:0050817_coagulation | 15/81 (0.0150) | ||
GO:0007596_blood_coagulation | 15/81 (0.0150) | ||
GO:0006412_protein_biosynthesis | 47/371 (0.0200) | ||
GO:0009059_macromolecule_biosynthesis | 49/409 (0.0133) | ||
GO:0006915_apoptosis | 8/356 (0.030) | 43/356 (0.0129) | |
GO:0012501 _programmed_cell_death | 4/358 (0.05) | 8/358 (0.015) | 43/358 (0.0111) |
GO:0008219 cell death | 4/377 (0.04) | 8/377 (0.020) | 45/377 (0.0200) |
GO:0016265 death | 4/379 (0.03) | 8/379 (0.015) | 45/379 (0.0125) |
Underexpression | |||
GO:0006261_DNA-dependent_DNA_replication | 3/58 (0.0100) | 8/58 (0.0000) | |
GO:0006260_DNA_replication | 3/113 (0.0118) | 14/113 (0.0000) | |
GO:0006259 DNA metabolism | 5/365 (0.0075) | 23/365 (0.0000) | |
GO:0006270_DNA_replication_initiation | 5/14 (0.0000) | ||
GO:0006310 DNA recombination | 5/48 (0.0076) | ||
GO:0006139 nucleobase nucleoside nucleotide and nucleic acid metabolism | 49/1,892 (0.0000) | ||
GO:0043283_biopolymer_metabolism | 29/998 (0.0042) | ||
GO:0006281 _DNA_repair | 13/142 (0.0000) | ||
GO:0006974_response_to_DNA_damage_stimulus | 14/157 (0.0000) | ||
GO:0009719_response_to_endogenous_stimulus | 14/166 (0.0000) | ||
GO:0006950_response_to_stress | 24/755 (0.0045) | ||
GO:0007049_cell_cycle | 6/504 (0.0050) | 25/504 (0.0000) | |
GO:0000074_regulation_of_cell_cycle | 5/327 (0.0100) | 14/327 (0.0031) | |
GO:0051244_regulation_of_cellular_physiologic_process | 10/1,873 (0.0071) | 42/1,873 (0.0368) | |
GO:0000279_M_phase | 9/149 (0.0038) | ||
GO:0051301 cell division | 8/119 (0.0033) | ||
GO:0000910_cytokinesis | 8/119 (0.0033) | ||
GO:0000278_mitotic_cell_cycle | 10/177 (0.0040) | ||
GO:0000086_G2_M_transition_of_mitotic_cell_cycle | 4/37 (0.0245) | ||
GO:0000087_M_phase_of_mitotic_cell_cycle | 7/121 (0.0248) | ||
GO:0007067_mitosis | 7/119 (0.0258) |
. | Changed/total genes (FDR) . | ||
---|---|---|---|
GO category . | 2 h . | 6 h . | 12 h . |
Overexpression | |||
GO:0001944_vasculature_development | 10/46 (0.0227) | ||
GO:0001568_blood_vessel_development | 10/46 (0.0227) | ||
GO:0048514_blood_vessel_morphogenesis | 10/46 (0.0227) | ||
GO:0001525_angiogenesis | 10/46 (0.0227) | ||
GO:0050878_regulation_of_body_fluids | 16/98 (0.0238) | ||
GO:0007599 hemostasis | 15/87 (0.0209) | ||
GO:0042060_wound_healing | 15/87 (0.0209) | ||
GO:0009611_response_to_wounding | 41/324 (0.0200) | ||
GO:0050817_coagulation | 15/81 (0.0150) | ||
GO:0007596_blood_coagulation | 15/81 (0.0150) | ||
GO:0006412_protein_biosynthesis | 47/371 (0.0200) | ||
GO:0009059_macromolecule_biosynthesis | 49/409 (0.0133) | ||
GO:0006915_apoptosis | 8/356 (0.030) | 43/356 (0.0129) | |
GO:0012501 _programmed_cell_death | 4/358 (0.05) | 8/358 (0.015) | 43/358 (0.0111) |
GO:0008219 cell death | 4/377 (0.04) | 8/377 (0.020) | 45/377 (0.0200) |
GO:0016265 death | 4/379 (0.03) | 8/379 (0.015) | 45/379 (0.0125) |
Underexpression | |||
GO:0006261_DNA-dependent_DNA_replication | 3/58 (0.0100) | 8/58 (0.0000) | |
GO:0006260_DNA_replication | 3/113 (0.0118) | 14/113 (0.0000) | |
GO:0006259 DNA metabolism | 5/365 (0.0075) | 23/365 (0.0000) | |
GO:0006270_DNA_replication_initiation | 5/14 (0.0000) | ||
GO:0006310 DNA recombination | 5/48 (0.0076) | ||
GO:0006139 nucleobase nucleoside nucleotide and nucleic acid metabolism | 49/1,892 (0.0000) | ||
GO:0043283_biopolymer_metabolism | 29/998 (0.0042) | ||
GO:0006281 _DNA_repair | 13/142 (0.0000) | ||
GO:0006974_response_to_DNA_damage_stimulus | 14/157 (0.0000) | ||
GO:0009719_response_to_endogenous_stimulus | 14/166 (0.0000) | ||
GO:0006950_response_to_stress | 24/755 (0.0045) | ||
GO:0007049_cell_cycle | 6/504 (0.0050) | 25/504 (0.0000) | |
GO:0000074_regulation_of_cell_cycle | 5/327 (0.0100) | 14/327 (0.0031) | |
GO:0051244_regulation_of_cellular_physiologic_process | 10/1,873 (0.0071) | 42/1,873 (0.0368) | |
GO:0000279_M_phase | 9/149 (0.0038) | ||
GO:0051301 cell division | 8/119 (0.0033) | ||
GO:0000910_cytokinesis | 8/119 (0.0033) | ||
GO:0000278_mitotic_cell_cycle | 10/177 (0.0040) | ||
GO:0000086_G2_M_transition_of_mitotic_cell_cycle | 4/37 (0.0245) | ||
GO:0000087_M_phase_of_mitotic_cell_cycle | 7/121 (0.0248) | ||
GO:0007067_mitosis | 7/119 (0.0258) |
Abbreviations: GO, gene ontology; PIA, phosphatidylinositol ether lipid analogue; FDR, false discovery rate.
Akt dependence and responses to PIAs
Although a COMPARE analysis indicated the cytotoxicity of PIAs correlated with active Akt levels (13), a functional analysis between PIA cytotoxicity and Akt activation in an isogenic system has not been done. To address this, we created an H157 cell line stably transfected with constitutively active MyrAkt1. Although MyrAkt1-expressing cells showed lower basal levels of apoptosis as indicated by cleaved PARP and sub-G1 DNA content, apoptosis was further induced with PIA23 treatment (Fig. 4A). Similar results were observed when other apoptotic assays such as Annexin-V/PI costaining were used (data not shown). These findings were confirmed in an A549 isogenic system, in which the 3 Akt isoforms were individually stably knocked down by lentiviral infection with shRNAs. Immunoblotting confirmed Akt isoform–specific knockdown, and also showed that Akt1 was the major isoform in A549 cells, because only Akt1 knockdown decreased levels of total and phospho-Akt. Accordingly, only Akt1 knockdown resulted in significantly less apoptotic cell death with PIA treatment (Fig. 4B). These studies showed levels of active Akt, specifically Akt1, correlated with PIA cytotoxicity.
To address the Akt dependence of PIA-induced genes, we used genetic or pharmacologic approaches to modulate Akt, and measured levels of RhoB, KLF6, and p21 after PIA treatment. In H157 cells transfected with MyrAkt1 or vector, induction of RhoB, KLF6, or p21 by PIA23 was observed (Fig. 4C). Although the induction of KLF6 and p21 by PIA23 in MyrAkt1-transfected cells appeared slightly diminished compared with vector-transfected cells, this was likely an artifact related to lower levels of p42 mitogen-activated protein kinase (MAPK) under these experimental conditions, which was observed in replicate experiments and was used as a loading control. When Akt1 was knocked down in A549 cells, the induction of RhoB, KLF6, and p21 by PIA23 was not affected (Fig. 4D). To confirm these results, we pretreated H157 cells with LY for 30 minutes followed by 6 hours of treatment with PIA23. LY alone slightly induced RhoB, KLF6, and p21 protein levels, but the combination of LY with PIA23 enhanced the expression of the PIA-induced genes over either compound alone (Fig. 4E). These results indicate that induction of these tumor suppressors is only minimally dependent upon the Akt pathway.
An important question is whether any of PIA-induced genes identified contribute to the cytotoxicity of the compounds. To examine this, H157 cells were transiently transfected with RHOB, KLF6, or CDKN1A siRNAs and treated with PIA 48 hours later. Cell lysates were harvested after 6 hours to assess knockdown, and sub-G1 DNA analysis was conducted after 12 hours of PIA treatment. The results show that although the siRNAs did not entirely block the induction of their target genes, these greatly rescued H157 cells from apoptosis caused by PIA (Fig. 4F). In contrast, overexpression of these genes either individually or in combination significantly decreased the viability of H157 cells (Fig. 4G). Similar results were observed in other NSCLC cell lines such as H1155 and H2882, and in other cancer cell lines with high levels of endogenous Akt activation (data not shown). These data confirm that RhoB, KLF6, and p21 induction contribute to the cytotoxicity of PIAs.
Discussion
Using microarray analysis, we identified gene expression profiles that contribute to the biological effects of PIAs. Validation of individual gene changes using RT-PCR and immunoblotting showed that microarray results could be validated at both the mRNA and protein levels, albeit protein level decreases were delayed as compared with mRNA decreases for the PIA-suppressed genes. The fact that induction of genes by PIAs could be measured by PCR or immunoblotting suggests that these genes could serve as biomarkers for PIA administration.
To place these individual changes in gene expression in a biological framework, GO analysis was conducted and revealed that many cellular processes are altered due to PIA-induced changes in gene expression in a time-dependent manner. Early induction of apoptosis or cell death and repression of DNA replication and cell cycle were observed after PIA administration, which is consistent with development of PIAs as anticancer agents. Of the early-induced genes, KLF6, RHOB/RhoB, and CDKN1A/p21 were of particular interest because they are known tumor suppressors and their expression reduced overall cell viability and contributed to PIA-induced cytotoxicity.
RhoB is a small GTPase tumor suppressor that regulates actin organization and vesicle transport. It is required for signaling apoptosis in transformed cells that are exposed to chemotherapeutic agents, it has a negative modifier function in carcinogenesis (33), and its expression is repressed during NSCLC progression (34–36). Although RhoB has a reciprocal relationship with levels of Akt activation in cells (37, 38), our studies showed that inhibition of Akt through genetic or pharmacologic means did not significantly affect induction of RhoB by PIAs. This suggests a novel mechanism for RhoB induction by PIAs.
CDKN1A/p21 inhibits the cell cycle by binding to cyclin–cyclin-dependent kinase (CDK) complexes and PCNA in nucleus, which can be prevented by Akt via phosphorylation (39). Because induction of p21 was observed in cell lines that have either wild-type (A549) or mutant p53 (H157, H1703, and H1155), induction of p21 by PIAs is p53 independent. This observation is consistent with p53-independent induction of p21 by perifosine, an alkylphospholipid whose activity correlated highly with PIAs (13, 40). Interestingly, the induction of RhoB and p21 by PIAs might be related because PPARγ-mediated induction of p21 in anaplastic thyroid carcinoma is dependent upon upregulation of RhoB (41).
KLF6 is a member of the Krüppel-like factor family of C2H2-type zinc finger–containing transcription factors implicated in cellular differentiation and tissue development (42). The KLF6 gene encodes a family of proteins generated through alternative splicing, which results in at least 4 isoforms. Full-length KLF6 is a tumor suppressor that is frequently inactivated by LOH, somatic mutation, and/or decreased expression in human cancer. Its tumor suppressor roles are not completely known, but a number of highly relevant activities have been described such as transactivation of p21 in a p53-independent manner (43, 44), reduction of cyclin D1–CDK4 complexes via interaction with cyclin D1 (45), and induction of apoptosis through upregulation of ATF3 (46).
Genetic manipulation showed that active Akt1 is a predictor of PIA efficacy in NSCLC cells, as expression of constitutively active Akt increased cytotoxicity and knockdown of Akt1 protected against PIA-induced cytotoxicity. PIA treatment induced expression of KLF6, RHOB/RhoB, and CDKN1A/p21 in a mostly Akt-independent manner. Nonetheless, these tumor suppressors inhibited cellular viability when overexpressed and contributed to the cytotoxicity of PIAs. These studies suggest that the broad activity of PIAs against cancer cells is based upon multiple independent mechanisms that include inhibition of Akt (12) and activation of p38α (47) and AMPKα (26), as well as induction of tumor suppressor genes such as KLF6, RhoB, and p21. As a result, the development of biomarkers for PIAs will be complex and will need to account for PIA-induced changes in enzymatic activity as well as gene transcription.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Grant Support
This research was supported in part by the Intramural Research Program of the NIH, NCI, Center for Cancer Research, and in part with federal funds from the NCI, NIH, under contract NCI Contract NO1-CO-12400.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.