| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Regular Articles |
1 Center for Genomics and Bioinformatics and
2 Department of Cell and Molecular Biology, Karolinska Institutet, Stockholm, Sweden, and
3 Department of Pharmacology, Pharmacia, Pfizer Group, Nerviano, Italy
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Cellular senescence has been suggested to be important for cancer surveillance, because the majority of cells derived from tumors maintain their telomeres at a stable length and thereby abrogate the induction of senescence. This is achieved either by expressing telomerase, the enzyme that extends the telomeric repeat, or by using an alternative mechanism for telomere maintenance (ALT; Ref. 4 ). Recent data have shown that lymphomas enter senescence upon treatment with cyclophosphamide, indicating that induction of senescence could be an important mechanism by which chemotherapy operates, and is therefore a promising strategy for new cancer treatments (5) .
Many tumor suppressors (p21, p16, and p53) are activated during the induction of senescence, often in a cell type-specific manner, and are believed to be the downstream targets of telomeric erosion. However, the understanding of the pathway(s) leading to senescence is still far from complete. Isolation of putative genes regulating the senescence machinery downstream of p53, p16, and p21 is complicated by the fact that primary cells in culture enter senescence asynchronously over prolonged passaging. This results in difficulties in finding key genes that are early triggers or effectors of senescence, especially because they may be active only during the initial phase of the senescence transition and then turned off. Because there are cell type-specific pathways to induce senescence and a large proportion of all cancers are of epithelial origin, there is particular interest in the regulation of epithelial cell senescence.
The SV40 large T antigen cellular model applied in this study provides the possibility of studying senescence in a highly controlled manner in epithelial cells. The use of the temperature-sensitive SV40 T antigen allows inhibition of the normal induction of senescence at the permissive temperature and induction of a proliferative block that resembles senescence by raising the temperature. The appeal of this model includes its temporal synchrony and the fact that it involves activation of pathways similar to those induced by telomere erosion (p53 and RB). We studied the transcriptome of two mouse thymic epithelial cell lines immortalized with either wild-type SV40 T antigen (Epi-A1 cells) or a temperature-sensitive mutant SV40 T antigen (Epi-A1ts58 cells; Refs. 6 , 7 ). When Epi-A1ts58 cells are grown at the permissive temperature (33°C), the SV40 T antigen causes the cells to escape senescence. At 39°C, the SV40 T antigen is inactivated, and the cells enter a state that resembles senescence, as judged by morphological changes (enlarged and flattened shape), a proliferative block, and expression of senescence-associated ß-galactosidase (7) . The Epi-A1 cells (wild-type SV40 T antigen) are used to separate the effect of the inactivated SV40 T antigen from the effect of the heat shock treatment used to inactivate the temperature-sensitive SV40 T antigen.
We identified some groups of genes, e.g., genes related to transforming growth factor (TGF)-ß signaling, based on their functions and time of induction during the establishment of senescence, including several previously uncharacterized genes.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Flow Cytometry.
Epi-A1ts58 cells and Epi-A1 cells were grown to 20% confluence. Prewarmed medium at either 39°C or 33°C was added to a set of dishes from both Epi-A1ts58 cells and Epi-A1 cells, and the plates were transferred to 39°C or 33°C. Cells were harvested after 0, 6, 12, 18, 24, and 32 h for further analysis. Cells were stained with propidium iodide and analyzed using a fluorescence automated cell sorter (Becton Dickinson).
BrdUrd Assay.
Epi-A1ts58 cells were plated on glass coverslips at 50% confluence. At time 0, the medium was changed (prewarmed to 33°C or 39°C), and the cells were incubated at 33°C or 39°C. BrdUrd was added to the media at -2, 0, 2, 4, 6, or 8 h, and the cells were harvested 2 h later by fixation followed by staining for BrdUrd and Hoechst 33342, according to the instructions of the manufacturer (BrdUrd labeling kit; Roche). The BrdUrd incorporation was analyzed using a fluorescence microscope (Leica DMRXA). More than 200 cells were counted for each time point, and the experiment was repeated twice.
Heat Inactivation and mRNA Preparations.
For each heat inactivation, the Epi-A1 and Epi-A1ts58 cells were plated on day 1 at 20% confluence, and the heat inactivation was started on days 23, when cells had reached 50% confluence except for the 72-h heat inactivation of the Epi-A1 cells for which the heat inactivation was initiated at 20% confluence to avoid high confluence after 72 h. All processes of heat inactivation were initiated by changing medium to fresh medium prewarmed to 39°C and maintained by incubating cells at 39°C and 5% CO2 for the desired time (1, 2, 4, 10, or 72 h). To avoid a possible serum bias in the study, the cells for the 72-h heat inactivation were subjected to medium change 3 h before harvest. At harvest, the mRNA was extracted using the PolyAtract 1000 kit (Promega). All experiments were performed two to four times.
Expression Analysis.
Two to four independent mRNA preparations from each condition were labeled separately, resulting in independent biological replicates. After labeling, the samples were hybridized, washed, and scanned as described previously (8)
. The fragmented cRNA was hybridized to the mouse U74Av2 chip (Affymetrix).
Data Analysis.
Chips were analyzed according to the following standards. Data were extracted using the Micro Array Suite 5 (MAS5). A linear normalization approach based on the 50th percentile of each chip and a per-gene normalization to the median was used. To identify up-regulated genes, all genes at a given time of heat shock (j) of the cell types Epi-A1 (WT) and Epi-A1ts58 (TS) were filtered according to the listed criteria (each filter is combined with AND, and the number of genes that pass each added filter is shown): 1. Select genes where
2 measurements were assigned "Present" or "Marginal" by MAS5 [7430 genes/expressed sequence tags (ESTs)]; 2. TSj/TS0 > 2 (1198 genes/ESTs); 3. TSj/WTj > 2 (404 genes/ESTs); 4. TSj/WT0 > 2 (342 genes/ESTs); 5. TSj
TS0 (Students t test, P = 0.05; 303 genes/ESTs); 6. TSj
WTj (Students t test, P = 0.05; 231 genes/ESTs); and 7. TSj
WT0 (Students t test, P = 0.05; 215 genes/ESTs).
The 215 genes (with a per-gene normalization to the median applied to get a clustering, reflecting kinetic effects on single genes rather than absolute differences between different genes) were pooled and clustered using the self-organizing maps (SOM) algorithm to generate six groups (AF). Several options were tested, and the latter turned out to generate homogeneous groups regarding expression patterns. Down-regulated genes were considered positive if they passed the following filters (each filter is combined with AND and the number of genes that pass each added filter is shown): 1. Select genes where all four measurements for TS0 were assigned "Present" or "Marginal" by MAS5 (5360 genes/ESTs); 2. TS72/TS0 < 0.5 (505 genes/ESTs); 3. TS72/WT72 < 0.5 (271 genes/ESTs); 4. TS72/WT0 < 0.5 (256 genes/ESTs); 5. TS72
TS0 (Students t test, P = 0.05; 253 genes/ESTs); 6. TS72
WT72 (Students t test, P = 0.05; 226 genes/ESTs); and 7. TS72
WT0 (Students t test, P = 0.05; 198 genes/ESTs).
All down-regulated genes (with a per-gene normalization to the median to obtain a clustering, reflecting kinetic events of single genes rather than absolute differences between genes) were clustered using the SOM algorithm into two groups (G and H) as homogeneous populations of genes were generated.
We further classified all positive genes using data derived from the confluence control (con), based on the following criteria: TSj
TScon (Students t test, P = 0.05) and the maximum or minimum relative expression compared with the confluence control was calculated for each gene.
Quantitative Real-Time PCR.
Epi-A1 cells and Epi-A1ts58 cells were heat shocked for 0, 6, 18, and 48 h according to the same scheme as above. At harvest, total RNA was isolated with RNeasy Mini (Qiagen) and treated with RNase-Free DNase (Qiagen) according to the manufacturers protocol. First-strand cDNA synthesis was performed according to the manufacturers protocol (Applied Biosystems). For a reaction volume of 80 µl, 1.6 µg of total RNA were used as template. Later, 1/80 aliquots of the cDNA reactions were analyzed by quantitative real-time PCR using an ABI PRISM 7000 (Applied Biosystems) according to the manufacturers protocol. Briefly, gene-specific primers for the target genes were mixed separately with SYBR Green qPCR Mastermix Plus (MedProbe) and added to a 96-well plate containing the cDNA to be analyzed. Samples were run in triplicates, and the data obtained were analyzed with ABI Prism SDS Software (Applied Biosystems). The analysis included a melting temperature-dependent dissociation curve of the amplicon. Sequences of the primers were: Tes-1, 5'-cagcggcctcaacaatgtc-3' (forward) and 5'-cttccttgtccggctcgtt-3' (reverse); Pem, 5'-aaggagaagaattaaatggaggaaaag-3' (forward) and 5'-ttgtccccatcacccatagg-3' (reverse); p21, 5'-caaagtgtgccgttgtctcttc-3' (forward) and p21 5'-tgagcgcatcgcaatca-3' (reverse); Tgfb1, 5'-aaacggaagcgcatcgaa-3' (forward) and 5'-gggactggcgagccttagtt-3' (reverse); Tgfb3, 5'-gccaaagagatccataaattcga-3' (forward) and 5'-ggcagacggccagttcat-3' (reverse); Foxd1, 5'-acatcgcgctcatcaccat-3' (forward) and 5'-tgctgatgaactcgcagatctc-3' (reverse); Klf4, 5'-agggagaccgaggagttcaac-3' (forward) and 5'-tcctggtgggttagcgagtt-3' (reverse); Ccng1, 5'-gcgactgaagaggaaaggaatg-3' (forward) and 5'-ggtctgaaaccgtgaacctatactg-3' (reverse); Itga5, 5'-agaggagcctgtggagtacaagtc-3' (forward) and 5'-gtggagcacatgccaagatg-3' (reverse); Osmr, 5'-ttccttagtgacacagggacaaac-3' (forward) and 5'-cgagacaaagagaacagtaccaaatatt-3' (reverse); Cnn2, 5'-ccggctcctgtccaaatatg-3' (forward) and 5'-cccgtgagtccctctatcca-3 (reverse); Igfbp3, 5'-gcaggcagcctaagcaccta-3' (forward) and 5'-cctcctcggactcactgatgtt-3' (reverse); and Cyclophilin A, 5'-tgtgccagggtggtgactt-3' (forward) and 5'-caaatttctctccgtagatggacct-3' (reverse).
| RESULTS |
|---|
|
|
|---|
|
Identification of Senescence-related Genes.
Epi-A1ts58 cells and Epi-A1 cells were transferred from 33°C to 39°C, and at the time points described above, mRNA was isolated for analysis using U74Av2 chips (Affymetrix). The data obtained were analyzed to identify genes/ESTs that are induced or repressed in the Epi-A1ts58 cells but not in Epi-A1 cells upon SV40 T antigen inactivation. This assured that the observed changes were related to the growth arrest and not to the heat shock. Validity of the dataset was indicated by the expression pattern of known markers of senescence. These included Clu (Clusterin/ApoJ; Fig. 1E
), a gene found previously to be induced in the SV40 T antigen senescence model in rat (10)
, and p21 whose expression was drastically induced within 14 h (Fig. 1F)
and maintained slightly overexpressed at 72 h. The rapid induction of p21 mRNA is probably a consequence of a fast inactivation of the SV40 large T antigen, allowing p53 to activate downstream targets shortly after the heat inactivation is initiated. In addition, several genes related to cell cycle progression (e.g., cyclin B2, cyclin A2, and cdc2a) and replication (e.g., cdc6, rrm1, and mcmd2) were down-regulated at 72 h, some of which have already been identified as differentially expressed in the SV40 T antigen senescence model, such as cyclin A2 and cdc2 (12)
.
After applying our filters as described in the "Materials and Methods" section, 215 genes/ESTs were identified as induced and 198 genes/ESTs as repressed after temperature shift of the Epi-A1ts58 cells. The induced genes were clustered using the SOMs algorithm into six groups (see "Materials and Methods"; Fig. 2
): four groups were early-induced genes/ESTs (groups CF), and two groups were late-induced genes/ESTs (groups A and B). The repressed genes were clustered into two groups using the SOM algorithm (groups G and H; only genes with known names are shown in Fig. 2
).
|
To separate putative senescence-regulated genes from genes related to other types of growth arrest, we classified the induced and repressed genes as being confluence related or not by also analyzing the expression in Epi-A1ts58 cells grown to confluence. The maximum relative expression compared with the confluent culture is shown in Fig. 2
for the genes that show a difference using t-statistics. Sixty % of all identified genes/ESTs showed no difference compared with confluent cells and are therefore not senescence specific. Fewer of the up-regulated genes/ESTs (42%) compared with the down-regulated genes/ESTs (80%) showed an expression pattern similar to that of cells grown to confluence. A total of 125 up-regulated and 39 down-regulated genes/ESTs were specific for the SV40 T antigen model compared with cells grown to confluence, and some were isolated more than once. In support of this approach to classify genes as putative senescence related or not, Clusterin/ApoJ was found to be specifically up-regulated during SV40 T antigen inactivation, whereas cyclin A2 and cdc6, for example, were down-regulated in the confluent culture as well.
Genes That Are Induced Early, 14 h after SV40 T Antigen Inactivation.
To get an overview of the induced and repressed genes, we have attempted to classify the identified genes based on the expression patterns we obtained and their known functions. Among the early-induced genes (Fig. 2, CF)
at least two classes can be identified, being aware that a functional overlap between the classes could occur. The first class includes Pparg (14)
, Nedd1 (15)
, Gata3 (16)
, M-twist (17)
, Dlx2 (18)
, Pem (19)
, and Sprr1a (20)
, which have been associated with differentiation and are not induced in confluent cells. Five of these (Pparg, Gata3, M-twist, Dlx2, and Pem) are transcription factors, and the expression of Dlx2 is maintained also after 72 h. This initial induction of genes related to differentiation could indicate that senescence involves the onset of differentiation programs that also involves irreversible growth arrest. An alternative explanation is that SV40 large T antigen blocks differentiation and inactivation of SV40 large T antigen and therefore induces onset of a senescence-like differentiation program.
The second class of early-induced genes has been related to growth arrest [Ddit3 (gadd153, Chop10) (21) , Dusp6 (22) , Cnn2 (23) , Nab2 (24) , and Akap12 (SSeCKS) (25) ]. None of these but Dusp6 and Nab2 is induced in confluent cells. Nab2 and Chop10 act by inhibiting transcription, SSeCKS is a scaffolding protein that binds G1 to S signaling molecules, Dusp6 inhibits the Erk1/2 pathway, and Cnn associates to actin filaments. Nab2, Dusp6, and Cnn are overexpressed after 72 h, which may suggest that they are involved in the manifestation of the irreversibly growth-arrested phenotype, whereas the others could have a function in the initial cessation of proliferation. The lack of overlap in gene expression between senescent and confluent cultures within this class could be explained by induction of different pathways or that confluent cells have passed the initial phase of growth restriction.
A gene that does not fit into the above classes and is induced early is Ccng1 (cyclin G1). Cyclin G1 is the only cyclin that is induced upon SV40 T antigen inactivation, and the expression is maintained after 72 h. Cyclin G1 has been suggested to be part of a feedback loop because it is a putative p53-responsive gene known to affect the MDM2-mediated regulation of p53 (26) . Another gene that is induced early and whose expression is maintained after 72 h is interleukin 6. Interleukin 6 is up-regulated in SAM mice (27) , which is a model for accelerated aging. The induction of interleukin 6 (not seen in confluent cells) implies the possibility that senescence could be induced or maintained by a secreted factor.
Genes That Are Induced Late, 1072 h after SV40 T Antigen Inactivation.
Several secreted factors are found among the later-induced genes (Fig. 2, A and B)
, which suggests that secreted molecules indeed could be important for manifestation of the senescent phenotype. The induction of TGF-ß3 is the best example and attributable to the complexity of the signaling downstream of TGF-ß, which is so highly cell type- and cell stage-specific that it could act as an initiator of transcription important for senescence. A role for the TGF-ß pathway in the SV40 T antigen senescence model is further suggested by the up-regulation of known TGF-ß regulated genes. Dpt (dermatopontin; Ref. 28
) has been reported as induced by TGF-ß but is not specific for T antigen-mediated growth arrest in our experiment. Igfbp3 (insulin-like growth factor binding protein 3) is specifically induced in the T antigen model and can be induced by TGF-ß (29)
but also by p53. Igfbp3 is induced during senescence (30)
, and the mechanism of Igfbp3-induced growth arrest in a cancer cell line is dependent on the TGF-ß pathway (31)
. Other identified genes modulate the TGF-ß pathway: a non T antigen-specific induction is observed for Fmod (fibromodulin) and Dcn (decorin), which both can bind to TGF-ß and sequester it to the extracellular matrix, but decorin may have other functions related to growth arrest mediated through the Akt/protein kinase B pathway (32)
.
Among the late-induced genes (after 10 or 72 h), we could identify at least two other classes in addition to the class including TGF-ß-related genes; induced expression of extracellular matrix components is a hallmark of senescent cells. This is illustrated by the isolation of several different forms of collagen. Col3A1, Col6A1, Col6A2, and Col6A3 are all induced upon SV40 T antigen inactivation and confluence mediated growth arrest, whereas Col5A1 and Col4a1 are specific for growth arrest induced by SV40 T antigen inactivation. The second class of late-induced genes contains genes related to apoptosis and/or growth arrest, but the function is sometimes also related to terminal differentiation. Stat3 and Itga5 have been described to have growth-regulatory as well as antiapoptotic effects, Cryab (B-crystallin) to have antiapoptotic functions, and AK1 (adenylate kinase) and Klf4 (Kruppel-like factor 4) to have growth arrest-related functions. Stat3 has been associated with G0 in some cells (33)
and a constitutively active Stat3-inhibited apoptosis (34)
, but Stat3 is also induced during confluence in our experiment, indicating a function in the establishment of G0 rather than senescence-related functions. Itga5 (integrin
5) is specifically induced during SV40 T antigen-mediated growth arrest and has been shown to activate the expression of gas1 (35)
as well as to suppress apoptosis in colon cancer cells (36)
. B-crystallin can inhibit caspase-3 (37)
. AK1 has been identified as a p53-responsive gene (38)
but is only slightly induced in the SV40 T antigen model compared with a confluent culture (1.6-fold). Klf4 is a transcription factor that is specific for SV40 T antigen inactivation-mediated growth arrest and has been shown to mediate p53-mediated G1-S arrest (39)
.
A few more interesting genes that do not fit into the classes mentioned above are induced late in the SV40 T antigen model. Osmr (oncostatin M receptor) is a putative p53-induced gene that is not induced during confluence. Osmr has been shown to inhibit breast cancer cell growth (40) , sometimes accompanied by an induction of p21 (41) . Introduction of Jup (plaktoglobin) into p53-/- renal carcinoma cells suppressed tumorigenicity (42) .
Genes That Are Repressed upon SV40 T Antigen Inactivation.
There is a big overlap among the down-regulated genes/ESTs between growth arrest induced by SV40 T antigen inactivation and growth arrest induced by confluence (80%). Most of the down-regulated genes are related to the cell cycle, either as molecules important for replication (e.g., Mcmd5, Mcmd4, Mcmd2, Rrm1, Rrm2, Top2a, and Cdc6) or controlling other stages of the cell cycle [e.g., Ccnb2 (cyclin B2), Ccna2 (cyclin A2), Cdc25c, and Cdc2a]. Among these, mcmd5, Ccnb2, Cdc2a, and Top2a are only repressed in the SV40 T antigen model, whereas the other genes are repressed in confluent cells as well. Several other genes are markers for proliferation: Plk (polo-like kinase) is important for mitosis and has been described as a p21-repressed gene (together with Top2a; Refs. 43
, 44
), and Stk18 (serine/threonine kinase 18/sak, another polo-like kinase family member) is important for progression of mitosis (45
, 46)
, which is also true for Stk5 (aurora-related kinase 2; Refs. 47
, 48
). Another SV40 T antigen model-specific gene is correlated to apoptosis. Tia1 is a RNA-binding protein that regulates splicing (49)
and is an apoptosis-promoting factor (50)
.
Taken together, these results indicate that although there is a big overlap between SV40 T antigen inactivation-induced growth arrest and growth arrest induced by confluence among the down-regulated genes/ESTs, several genes with functions related to cell cycle regulation or progression are specifically repressed in the SV40 T antigen model and could therefore be important for the irreversible phenotype. However, there are several kinases that clearly associate with proliferation and are cell cycle regulated among them; hence, down-regulation is rather a consequence of growth arrest.
Validation of Differentially Expressed Genes by Quantitative Reverse Transcription-PCR.
To validate the obtained expression patterns and to obtain an estimate of the number of false positives, we selected 10 genes that we considered biologically interesting and performed quantitative real-time PCR. We included p21 as a positive control and TGF-ß1 as a negative control, because TGF-ß1 was not identified in our study, although TGF signaling was implied as a possible mediator of senescence induction. The mRNA levels in both Epi-A1ts58 and Epi-A1 cells after 0, 6, 18, and 48 h of heat inactivation of the SV40 large T antigen were measured. The time points were chosen to include an early, an intermediate, and a late time point. To obtain kinetic information as well as a relative comparison between the cell lines at different times of heat inactivation, the data were normalized to Epi-A1ts58 cells before heat inactivation. As shown in Fig. 3, AL
, 9 of 10 of our selected genes were verified (showing both a 2-fold induction as well as a 2-fold increase compared with the maximal Epi-A1 heat-related response), indicating an approximate false-positive discovery fraction of 10%. As shown in Fig. 3E
, the expression pattern of TGF-ß3 is verified. Interestingly, TGF-ß1 (Fig. 3D)
is induced in the Epi-A1ts58 but also during the heat shock in the Epi-A1 cells (this is in agreement with a microarray study, data not shown). This could indicate an overlap between the heat shock response and senescence and may suggest a risk of false-negative genes appearing in the analysis that are shared between the heat shock response and senescence. Some genes seem to be increased by heat shock (Pem, Tgfb3, and Cnn2) but never reach the levels seen in the Epi-A1ts58 cells, possibly indicating additive effects from both heat shock and SV40 large T inactivation. The expression levels of other genes in the Epi-A1ts58 and the Epi-A1 cells are not at the same level in the two cell types before the heat inactivation had been initiated (Tes-1, Pem, Foxd1, Ccng, Cnn2, or Igfbp3). This indicates some inefficiency of the temperature-sensitive form of SV40 large T antigen at 33°C. A similar conclusion can be drawn by the expression pattern for p21 (higher expression at 33°C in the Epi-A1ts58 cells compared with the Epi-A1 cells; Fig. 3C
) and appearance of spontaneous senescence of the Epi-A1ts58 cells at 33°C (<5%, data not shown), as well as the population doubling time for the Epi-A1 cells (24 h) compared with the Epi-A1ts58 cells (40 h).
|
| DISCUSSION |
|---|
|
|
|---|
In contrast to most other comparable studies, we used epithelial cells. This is important from a cancer perspective because the majority of human cancers are of epithelial origin, and induction of senescence is cell type specific. The model was established in mouse cells because mouse cells with overexpressed SV40 large T antigen do not enter crisis, which occurs when human cells that have bypassed senescence reach a critically short telomere length with genomic instability as a consequence. Using a mouse cell model to study senescence compared with a model established in human cells could affect the interpretation of the results as replicative in vitro senescence and is believed to be triggered by telomere erosion in human cells. In contrast, mouse cells have longer telomeres and normal replicative senescence in mouse cells is induced by O2 (51) . However, although the initiation signals for in vitro senescence caused by serial passaging might differ between human and mouse cells, execution of mouse senescence is likely to be similar to human senescence. This notion is supported by several findings, e.g., that other stimuli than short/eroded telomeres can induce senescence in human as well as mouse cells (e.g., overexpression of RAS or chemotherapy).
We had two primary goals in our expression study of the temperature-sensitive SV40 large T antigen model. The first goal was to study early events that lead to senescence rather than the stationary senescent phenotype, and the second goal was to separate primary genes related to the induction of the senescence phenotype from secondary genes related to the growth arrest that is a consequence of senescence. To accomplish these goals we studied early time points of senescence induction and also made a comparison with cells growth arrested by high confluence. By the addition of a later time point where the growth arrest is irreversible, we were also able to classify genes as being induced early or late during our study.
In our screen for differentially expressed genes with related functions and expression kinetics, we found that several of the identified transcripts were related to TGF-ß signaling. This is interesting because recent data indicate a role for TGF-ß in other senescence models as well. TGF-ß has been shown to be important in senescence induced by H2O2 (52) and oncogenic RAS (53) and has been correlated to the induction of Clusterin/ApoJ as well as other senescence-associated genes in H2O2-induced senescence (52) . A possible mechanism of TGF-ß in senescence is indicated by the ability of TGF-ß to activate the p38 mitogen-activated protein kinase pathway (54) as p38 has been shown to be necessary for RAS-induced senescence (55) as well as H2O2-induced senescence (56) . Interestingly, we found that although induction of TGF-ß3 seems to be specific for senescence in the SV40 large T antigen model, TGF-ß1 was induced both in senescent cells as well as in the cells used as our heat shock control. Although an overlap between heat shock and senescence has not been reported, it can be argued that both senescence and heat shock can be considered as stress responses, and an overlap might occur.
Two major efforts have been made previously to identify senescence-related genes using the SV40 T antigen model. Both these studies used rat embryo fibroblasts and applied either a subtraction-based approach (10) or a two-dimensional gel approach (9) . Our analysis differs in three major ways from the previous studies: (a) we used cells of epithelial origin (the other studies used fibroblasts); (b) we emphasized that we wanted to identify early changes in gene expression (compared with stationary senescent cells used in the other studies); and (c) we used a global gene expression approach to find senescence-associated genes. In the previous approach applying subtractive cloning, a low number of positive genes indicates a less efficient system for finding differentially expressed genes compared with monitoring of global gene expression. The two-dimensional gel approach has the advantage of studying protein levels rather than using transcription as means to approximate protein levels. The setbacks of the two-dimensional gel approach are mainly caused by the method, which is less quantifiable and enables fewer genes to be studied because of resolution problems. The genes isolated as candidate senescence genes in these two studies overlap slightly compared with the present study. The low degree of overlap could be caused by using different cell types (species), different methods, and most probably on our unique approach to identify early senescence-associated changes in gene expression.
Senescence has been studied using microarray approaches in several other models, e.g., a model where senescence is induced by adding doxorubicin to HCT116 cells (57)
or a model where a temperature-sensitive papilloma virus E2 is used to induce senescence in HeLa cells (58)
. Both of these models use human cells, but the induced senescence differs largely in terms of cell cycle profile in the senescent population that they generate. Senescence induced by introduction of doxorubicin leads to an arrest in G2, whereas induction of E2 activity leads to G1 arrest. Because our model involves mainly G1 arrest and, similar to the E2 model, does not involve any external stimuli except for the temperature change (in the E2 model the temperature is reduced), one would expect a larger overlap between the E2 model and our SV40 large T antigen model. To investigate whether this was the case and to study the similarities among the three studies in detail, we compared the genes identified as differentially expressed in these two studies to the genes identified as differentially expressed in the SV40 large T antigen model. To compare the studies we linked the differentially expressed genes from the doxorubicin and the E2 study to the corresponding genes on our chip using Netaffx. Of 703 differentially expressed genes/ESTs from the E2 study, we were able to link and find an overlap compared with our identified genes for 47 genes (represents 11% of the genes identified in the SV40 large T antigen model). Thirty of 47 genes/ESTs were down-regulated, and as expected, this group included mainly genes related to DNA replication or G2-M phase progression such as cdc20, cdc2a, Ttk, and Mad2l1. In our experiments, only Plk, Stk5, Prss1, Rad51ap1, and Cdc2a from the list of genes/ESTs shared with the E2 study were more repressed during senescence compared with confluence, and the remaining genes are therefore not senescence specific but rather a consequence of growth arrest. Seventeen genes were up-regulated in both the E2 model and the SV40 large T antigen model. Thirteen of these were senescence specific (not induced in confluent cells in our study), and among these were a few genes with possible regulatory functions [Osmr (oncostatin M receptor), Jup (plakoglobin), and Stat3]. Fifteen of the 17 shared up-regulated genes were found in clusters A and B in Fig. 2
(representing 13% of all genes in clusters A and B). Interestingly, there is an enrichment of senescence-specific genes (compared with the confluent culture) from 56% for all genes to 87% for the shared genes in clusters A and B. The overrepresentation of shared genes in clusters A and B is expected, because these clusters contain genes that are induced late in our SV40 large T antigen model and because a 72-h induction of senescence was used in the E2 study. The big overlap of senescence-specific genes (not regulated during confluence) among the induced genes could indicate that when a gene is shared between two models, it is more likely to have a senescence-specific role in either model.
When we linked the genes from the doxorubicin study with our identified genes, we found an overlap for 19 genes (
5% of our genes). The lower degree of overlap could be expected from the differences in methods to induce senescence and from the cell cycle profile generated (as discussed above). Interestingly, all of these overlapping genes were down-regulated in both studies, which indicates that there might be a difference between senescence induced by the physiological stimulus used in both the E2 and the SV40 large T antigen model compared with doxorubicin. The down-regulated pool that is shared between the SV40 large T antigen model and the doxorubicin model is characterized by genes that are associated with replication and cell cycle progression, such as Mcmd, Top2a, and Rrm1 as well as to DNA damage responses, such as Rad51. Only 8 genes/ESTs were differentially expressed (down-regulated) in all three studies [those known are Stk5, Mad2l1, Cdc2a, Bub1, Ts (thymidylate synthase), and Mcmd2], and they probably reflect that the cells are no longer dividing in any of the models because only Stk5 and Cdc2a were repressed more during senescence compared with confluence in our study (9.9-fold and 2.7-fold, respectively). The overlap between the SV40 large T antigen model and the E2 model or the doxorubicin model can be summarized in two major points:
(a) The overlap is greater compared with the E2 model than with the doxorubicin model (11% compared with 5%). This probably reflects the different methods used to induce senescence, but the percentage of shared genes is also affected by the data analysis approaches.
(b) The genes that are shared with the doxorubicin model are all down-regulated; hence, no similarities in the possible activating genes were found, which is in contrast to the E2 model where several of the shared up-regulated genes were of signaling nature and represent potentially senescence-inducing genes. The overlap of induced senescence-specific genes that are shared between the E2 model and the SV40 large T antigen model indicates coregulation of senescence in mouse and human cells.
The genes that are shared in several studies and not associated with growth arrest induced by confluence could be genes with a general role in senescence. However, further studies are needed to test these candidate genes as well as other genes identified in our study functionally for their contribution to the senescent phenotype.
| FOOTNOTES |
|---|
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.
Requests for reprints: Ola Larsson and Zicai Liang, Center for Genomics and Bioinformatics, Karolinska Institutet, Berzeliusväg 35, 171 77 Stockholm, Sweden. Phone: 46 8 728 63 81; Fax: 46 8 331 547; E-mail: ola.larsson{at}cgb.ki.se
Received 6/25/03. Revised 10/30/03. Accepted 11/ 7/03.
| REFERENCES |
|---|
|
|
|---|
-dependent and -independent pathways. J. Biol. Chem., 276: 38297-38306, 2001.
S progression by controlling the expression and cellular compartmentalization of cyclin D. Mol. Cell. Biol., 20: 7259-7272, 2000.
and IL-6 in the brain of senescence accelerated mouse (SAM) P8. Brain Res., 885: 25-31, 2000.[CrossRef][Medline]
5ß1 expression negatively regulates cell growth: reversal by attachment to fibronectin. Mol. Biol. Cell, 6: 725-740, 1995.[Abstract]
5 subunit in HT29 colon carcinoma cells suppresses apoptosis triggered by serum deprivation. Exp. Cell Res., 224: 208-213, 1996.[CrossRef][Medline]
B-crystallin negatively regulates apoptosis during myogenic differentiation by inhibiting caspase-3 activation. J. Biol. Chem., 277: 38731-38736, 2002.
genes by p21(WAF1/CIP1/SDI1). Cell Cycle, 1: 59-66, 2002.[Medline]
This article has been cited by other articles:
![]() |
O. Larsson, D. M. Perlman, D. Fan, C. S. Reilly, M. Peterson, C. Dahlgren, Z. Liang, S. Li, V. A. Polunovsky, C. Wahlestedt, et al. Apoptosis resistance downstream of eIF4E: posttranscriptional activation of an anti-apoptotic transcript carrying a consensus hairpin structure Nucleic Acids Res., September 11, 2006; 34(16): 4375 - 4386. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A. Timmons, O. Larsson, E. Jansson, H. Fischer, T. Gustafsson, P. L. Greenhaff, J. Ridden, J. Rachman, M. Peyrard-Janvid, C. Wahlestedt, et al. Human muscle gene expression responses to endurance training provide a novel perspective on Duchenne muscular dystrophy FASEB J, May 1, 2005; 19(7): 750 - 760. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||