Conserved genes and pathways in primary human fibroblast strains undergoing replicative and radiation induced senescence
Biological Research volume 49, Article number: 34 (2016)
Cellular senescence is induced either internally, for example by replication exhaustion and cell division, or externally, for example by irradiation. In both cases, cellular damages accumulate which, if not successfully repaired, can result in senescence induction. Recently, we determined the transcriptional changes combined with the transition into replicative senescence in primary human fibroblast strains. Here, by γ-irradiation we induced premature cellular senescence in the fibroblast cell strains (HFF and MRC-5) and determined the corresponding transcriptional changes by high-throughput RNA sequencing.
Comparing the transcriptomes, we found a high degree of similarity in differential gene expression in replicative as well as in irradiation induced senescence for both cell strains suggesting, in each cell strain, a common cellular response to error accumulation. On the functional pathway level, “Cell cycle” was the only pathway commonly down-regulated in replicative and irradiation-induced senescence in both fibroblast strains, confirming the tight link between DNA repair and cell cycle regulation. However, “DNA repair” and “replication” pathways were down-regulated more strongly in fibroblasts undergoing replicative exhaustion. We also retrieved genes and pathways in each of the cell strains specific for irradiation induced senescence.
We found the pathways associated with “DNA repair” and “replication” less stringently regulated in irradiation induced compared to replicative senescence. The strong regulation of these pathways in replicative senescence highlights the importance of replication errors for its induction.
DNA is the repository of genetic information in each living cell, its integrity and stability is essential to life. DNA, however, is not inert; rather, it is subject to assault from cell-internal and environmental processes. Any resulting damage, if not repaired, will lead to mutation and possibly disease.
Cell-internally, DNA is subject to oxidative damage from metabolic byproducts such as free radicals. In addition, DNA replication is prone to error. The rate at which DNA polymerase incorporates incorrect nucleotides into newly synthesized DNA drives spontaneous mutations in an organism. While polymerase proofreading normally recognizes and corrects many of these errors, some 10−4 to 10−6 mutations per gamete for a given gene survive this process. DNA damage is also induced by the cellular environment, for example by UV-light and radiation of the cell . An individual cell can suffer up to 106 DNA changes per day.
Cells have evolved a number of mechanisms to detect and repair the various types of DNA damage, no matter whether this damage is caused by the environment or by errors in replication and cell division. If the rate of DNA damage exceeds the capacity of the cell to repair it, the accumulation of errors can overwhelm the cell [2–11] and lead to mutations and potentially to cancer. After major damage, the cell induces self-destruction by necrosis or apoptosis [12–14]. As a functional alternative to apoptosis, less damaged or replicatively exhausted but functional cells become senescent (“Hayflick limit” [15, 16]), an irreversible cell cycle arrested state experienced by all mitotically competent cells. It results from an intrinsic natural barrier to unlimited cell division exhibited by all normal somatic cells, including fibroblast [17–20]. Several mechanisms and pathways, especially the p53–p21 and p16–pRB pathways and telomere processing are involved in cellular senescence induction [15, 21–37].
The induction of apoptosis and senescence is considered to be part of a cellular cancer protection strategy . Cellular senescence arrests the growth of cells at risk for malignant transformation in culture and in vivo [39–46], in this way preventing the spread of damage to the next cell generation . Senescent cells accumulate over the life span of rodents and primates  and are found primarily in renewable tissues and in tissues that experience prolonged inflammation. Senescence-associated changes in gene expression are specific and mostly conserved within individual cell types . Most differences between the molecular signatures of pre-senescent and senescent cells entail cell-cycle and metabolism-related genes , as well as genes encoding the secretory proteins that constitute the senescence associated secretory phenotype (SASP) [50–52].
Both, the accumulation of errors internally by replication and cell division (a slow process involving changes in telomeric processing) or externally by irradiation (comparably fast, not involving telomere shortening) can induce cellular senescence of practically indistinguishable phenotypes . We therefore speculated that in both cases, the transition into senescence may correlate with the differential regulation of similar genes. Human fibroblasts are a well-established model for the investigation of cellular senescence [5, 54–56]. Recently, we determined the transcriptional changes associated with the transition into replicative senescence . Here, by γ-irradiation we induced premature (accelerated) cellular senescence  in primary human fibroblast cell strains (HFF and MRC-5), determined the corresponding transcriptional changes by high-throughput RNA sequencing and compared the results with those for replicative senescence. Indeed, for both cell strains we found a high degree of similarity in differential gene expression in replicative as well as in irradiation induced senescence. However, we also identified that the senescence induction process imprints specific differences in the two transcriptomes.
Primary human MRC-5 fibroblasts (14 weeks gestation male, fibroblasts from normal lung, normal diploid karyotype) were obtained from ATCC (LGC Standards GmbH, Wesel, Germany). HFF (primary cells, Homo sapiens, fibroblasts from foreskin, normal diploid karyotype) cells were a kind gift of T. Stamminger (University of Erlangen, ).
Cells were cultured as recommended by ATCC in Dulbeccos modified Eagles low glucose medium (DMEM) with l-glutamine (PAA Laboratories, Pasching, Austria), supplemented with 10 % fetal bovine serum (FBS) (PAA Laboratories). Cells were grown under 20 % O2 levels in a 9.5 % CO2 atmosphere at 37 °C. For sub-culturing, the remaining medium was discarded and cells were washed in 1× PBS (pH 7.4) (PAA Laboratories) and detached using trypsin/EDTA (PAA Laboratories). Primary fibroblasts were sub-cultured in a 1:4 [=2 population doublings (PDs)] or 1:2 (=1 PD) ratio. For stock purposes, cryo-conservation of the cell strains at various PDs were undertaken in cryo-conserving medium (DMEM + 10 % FBS + 5 % DMSO). Cells were immediately frozen at −80 °C and stored for 2–3 days. Afterwards, cells were transferred to liquid nitrogen for long time storage. Re-freezing and re-thawing was not performed to avoid premature senescence .
One vial of each of the two fibroblast cell strains (MRC-5 and HFF) were obtained and maintained in culture from an early PD. After obtaining enough stock on confluent growth of the fibroblasts in 75 cm2 flasks, cells were sub-cultured into three separate 75 cm2 flasks (“triplicates”) and were passaged until they were senescent in culture. We analyzed “technical” replicates in order to determine the experimental error of our technical approach. When using three samples from independent stocks (“biological” replicates), these might already differ in their transcriptome and/or proteome when starting our analysis, making it difficult to estimate the error of our experimental procedure.
Induction of cellular senescence
Cellular senescence was induced by γ-irradiation. Human fibroblast strains were irradiated by ionizing radiation in a Gamma cell GC40 (MDS Nordion, Ottawa, Canada) using the radioactive isotope 137Cs as source. Exposure time was determined by correcting the irradiation dose of 1.23 Gy/min with the decay factor time equating to roughly 62 s/Gy. Young PD fibroblast strains (MRC-5 at PD 32, HFF at PD 16) were seeded 48 h before radiation exposure. Once cells were 70 % confluent, they were subjected to different doses of γ-irradiation (0, 2, 15, 20 Gy) at room temperature (RT) and subsequently cultured at 37 °C.
Detection of SA-β galactosidase activity
The SA-β Gal assay was performed as described by  in either of the fibroblast strains at different time spans (after 0, 24, 48, 72, 96 and 120 h) after subjecting them to different doses of γ-irradiation (0, 2, 15, 20 Gy). Cells were washed in 1× PBS (pH 7.4) and fixed in 4 % paraformaldehyde (pH 7.4), 10 min at RT. After washing the cells in 1× PBS (pH 7.4), staining solution consisting of 1 mg/ml X-Gal, 8 mM citric acid/sodium phosphate pH 6,0, 5 mM K3Fe(CN)6, 5 mM K4Fe(CN)6, 150 mM NaCl, 2 mM MgCl2, was added. The enzymatic reaction occurred without CO2 for 4–16 h at 37 °C. After incubation, cells were washed in 1× PBS (pH 7.4) and, in order to visualize cell nuclei, DNA and SAHFs, mounted with 4′-6-diamidine-2-phenyl indole (DAPI) containing prolong gold antifade reagent (Invitrogen, Carlsbad, USA). Paired two-sample type 2 Student’s t-tests, assuming equal variances, were applied to determine the statistical significance of the SA-β gal assay results.
For immunoblotting, 10,000 cells/µl were used. Immunodetection was performed using 5 %-powdered milk in PBS-T (1× PBS, pH 7.4 and 1 % Tween 20) for blocking (Roth, Germany). The optimal concentration of all the primary antibodies was estimated in human fibroblasts. Primary antibodies, anti-p21 mouse antibody (OP64; Calbiochem; dilution 1:200), anti-p16 mouse antibody (550834; BD Pharmingen; 1:200), anti-IGFBP7 rabbit antibody (ab74169; Abcam; 1:500), anti-IGFBP5 rabbit antibody (ab4255; Abcam; 1:500), anti-IGFBP3 goat antibody (ab77635; Abcam; 1:500), anti-Id3 mouse antibody (ab55269; Abcam; 1:100), anti-BAX rabbit antibody (ab10813; Abcam; 1:200), anti-Caspase-3 rabbit antibody (ab2302; Abcam; 1:500) and anti-tubulin mouse antibody (T-9026; SIGMA-Aldrich; 1:5000) were diluted in 5 %-powdered milk (in PBS-T) and incubated for 1 h at RT. Washing steps were performed three times for 10 min in 1× PBS-T. The secondary horseradish peroxidase-labeled antibodies (Jackson Immuno Research Lab) were incubated for 1 h at RT. Horseradish peroxidase was detected using an ECL-detection system and radiographic film (GE Healthcare, Germany). After film development, signal intensities of immunoblot bands were quantified using Metamorph software . The signal intensity values were examined for statistical significance using paired two-sample type 2 Student’s t-tests assuming equal variances.
Total RNA was isolated using Qiazol (Qiagen, Hilden, Germany) according to the manufacturer’s protocol, with modifications. In brief, the fibroblasts were pelleted in 2 ml safe-lock tubes (Eppendorf, Hamburg, Germany). 1 ml cooled Qiazol and one 5 mm stainless steel bead (Qiagen) were added. Homogenization was performed using a TissueLyzer II (Qiagen) at 20 Hz for 1 min. After incubation for 5 min at RT, 200 ml chloroform was added. The tube was shaken for 15 s and incubated for 3 min at RT. Phase separation was achieved by centrifugation at 12,000×g for 20 min at 4 °C. The aqueous phase was transferred into a fresh cup and 10 mg of glycogen (Invitrogen, Darmstadt, Germany), 0.16 volume NaOAc (2 M, pH 4.0) and 1.1 volume isopropanol were added, mixed and incubated for 10 min at RT. The RNA was precipitated by centrifugation with 12,000×g at 4 °C for 20 min. The supernatant was removed and the pellet was washed with 80 % ethanol twice and air dried for 10 min. The RNA was re-suspended in 20 μl DEPC-treated water by pipetting up and down, followed by incubation at 65 °C for 5 min. The RNA was quantified with a NanoDrop 1000 (PeqLab, Erlangen, Germany) and stored at −80 °C until use.
To ensure appropriate RNA quality and evaluate RNA degradation, total RNA was analyzed using Agilent Bioanalyzer 2100 (Agilent Technologies, USA) and RNA 6000 Nano Kit (Agilent). An average RNA integrity number (RIN) of 8 was obtained. Total RNA was used for Illumina library preparation and RNA-seq . 2.5 µg total RNA was used for indexed library preparation using Illumina’s TruSeq™ RNA Sample Prep Kit v2 following the manufacturer’s instruction. Libraries were pooled and sequenced (five samples per lane) using a HiSeq 2000 (Illumina) in single read mode with 50 cycles using sequencing chemistry v3. Sequencing resulted in approximately 40 million reads with a length of 50 bp (base pairs) per sample. Reads were extracted in FastQ format using CASAVA v1.8.2 or v1.8.3 (Illumina).
RNA-seq data analysis
Raw sequencing data were obtained in FASTQ format. Read mapping was performed using Tophat 2.0.6  and the human genome references assembly GRCh37 (http://feb2012.archive.ensembl.org/). The resulting SAM alignment files were processed using the HTSeq Python framework and the respective GTF gene annotation, obtained from the Ensembl database . Gene counts were further processed using the R programming language  and normalized to reads per kilobase of transcript per million mapped reads (RPKM) values. In order to examine the variance and the relationship of global gene expression across the samples, different correlation coefficients were computed including Spearman’s correlation of gene counts and Pearson’s correlation of log2 RPKM values.
Subsequently, the Bioconductor packages DESeq  and edgeR  were used to identify differentially expressed genes (DEG). Both packages provide statistics for determination of differential expression in digital gene expression data using a model based on the negative binomial distribution. Here we used non-normalized gene counts since both packages include internal normalization procedures. The resulting p values were adjusted using the Benjamini and Hochberg’s approach for controlling the false discovery rate (FDR) . Genes with an adjusted p value <0.05 found by both packages were assigned as differentially expressed.
In our study, we applied DESeq [67, 68] instead of the recently presented alternative tool DESeq 2. DESeq 2 results in minor differences to DESeq, however showing a slightly lower median precision . Applying the same statistical analysis tool (DESeq) for DEG identification allows a direct comparison of results in this study with those of our recent publications [35, 49, 70, 71].
Sample clustering and analysis of variance
The variance and the relationship of global gene expression across the samples were examined by computing the Spearman correlation between all samples using genes with raw counts larger than zero. Furthermore, principal component analysis (PCA) was applied using the log2 RPKM values for genes with raw counts larger than zero.
Gene set enrichment analysis to determine the most differentially regulated pathways on aging
We used the R package gage  in order to find significantly enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. In case of our RNA-seq data, the calculation was based on the gene counts and was performed as described in the methods manual. Estimated p values were adjusted for multiple testing using the Benjamini and Hochberg’s approach for controlling false discovery rate. KEGG pathways were selected as significantly regulated if the corrected p values were smaller than 0.05.
Results and discussion
Previously, changes in global gene expression have been studied during accelerated senescence induced by oncogenes in IMR-90 fibroblast strains [73, 74] or by chemotherapeutic drugs applied to tumor cells [37, 75] and during replicative and induced senescence in skin fibroblasts derived from Li-Fraumeni syndrome patients . Here, we compared the transcriptomes of two γ-irradiation induced senescent human primary fibroblast strains with the corresponding transcriptomes of the replicatively senescent cells.
Gamma irradiation resulted in senescence induction in primary human fibroblast strains
Mild irradiation (0.5 Gy) induces low levels of DNA damage in MRC-5 fibroblasts, followed by an increase of p21 protein levels [1, 51, 56]. After 3 days, the number of p21-positive cells drops to background levels, indicating successful DNA repair and return into the cell cycle. This mild irradiation did neither result in an increase of p16 protein levels nor in the up-regulation of the cellular senescence marker SA-β Gal . After a slight time lag, the cell population continued to grow with the same rate as before, consistent with cell cycle re-entry after a transient cell cycle arrest . After high dose irradiation (20 Gy), MRC-5 fibroblasts display a high number of repair foci which during the following days hardly decrease. After this high irradiation, not only p21 but now also p16 protein levels increase, associated with a complete cell proliferation arrest and a continuous increase of SA-β Gal positive cells . Here, we subjected two different human fibroblast cell strains of different tissue origins [HFF (foreskin) and MRC-5 (embryonic lung)] to γ-irradiation, inducing premature cellular senescence. We determined the transcriptome of these irradiation-induced senescent cells in order to compare it with that of replicatively senescent cells of the same strains obtained by us before [49, 70].
MRC-5 fibroblasts were irradiated by 0, 2, 15, and 20 Gy at room temperature. Then, the percentage of SA-β Gal stained cells was determined at different time points over 5 days after irradiation treatment (Fig. 1). The highest percentage of SA-β Gal stained MRC-5 fibroblast cells (63 ± 4 %) was noted after the highest irradiation dose (20 Gy) and the longest time lapse (120 h) . Therefore, HFF strains were only irradiated by 20 Gy. After 120 h the percentage of SA-β Gal stained HFF cells (62 ± 4 %; Fig. 2) was similar to the corresponding value for the MRC-5 fibroblasts. This time lapse, with the resulting degree of SA-β Gal staining, was selected since the transcriptomes of these cells will be compared to the corresponding transcriptomes of cells in replicative senescence of the same level of SA-β Gal staining (see below). 5 days after 20 Gy irradiation, MRC-5 cells are in early; partially still reversible but not yet in irreversible deep senescence . At this time point (120 h after 20 Gy irradiation), immunoblotting revealed that increase of BAX expression [35, 78, 79] was not induced by 20 Gy irradiation. For caspase-3 , the levels of their active (cleaved, 17 kDa) form were not increased by 20 Gy irradiation. Since induction of BAX and cleavage of caspase-3 are more consistent with an induction of apoptosis, our results reveal that apoptosis was not induced (Additional file 1: Figure S1).This finding is consistent with earlier observations [1, 2, 81, 82]. Total MRC-5 and HFF sample RNA was extracted 120 h after 20 Gy irradiation and was subjected to RNA-seq.
Transcriptome analysis of fibroblast strains subjected to irradiation induced senescence
Overall, the RNA-seq data were obtained from two samples, one for each cell strain (HFF and MRC-5), with three biological replicates each. The RNA-seq results revealed transcription of 27,410 and 27,944 genes for γ-irradiated HFF and MRC-5 fibroblasts, respectively. These were compared to the corresponding RNA-seq results of non-irradiated young (PD 16) HFF and young (PD 32) MRC-5 cells (obtained by us earlier [49, 70]). First, the RNA-seq retrieved normalized transcriptome expression values were analyzed using PCA. PCA reduces (by orthogonal transformation) high-dimensional data to 2 or 3 dimensions without losing much information, thereby enabling graphical visualization of the data. PCA is done in such a way that the first component of the graph shows as much of the variation contained in the data as possible. The PCA plot (Fig. 3) indicated a clear separation of the MRC-5 and HFF strains (PC2). The triplicates in all four cases clustered closely together, indicating small experimental errors (Fig. 3). The effect of irradiation induced senescence also exhibited similarities between the two fibroblast strains, demonstrated by the location and distance of both irradiated samples relative to the non-irradiated controls (irradiated samples on the right of controls; PC1).
In order to retrieve the most significant DEG, we applied stringent selection criteria of log2 fold change >1, p < 0.001 and adherence to both of the statistical packages edgeR and DESeq. 11,000 DEG in HFF and 6000 DEG in MRC-5 satisfied these selection criteria. An additional selection criterium of minimum RPKM >10 (in each of the samples compared; as mentioned in our previous studies [49, 70, 71]) resulted in more than 500 differentially regulated genes when comparing irradiated fibroblasts with their respective non-irradiated controls. Of these DEG, 29 % of the genes were commonly up- or down-regulated between HFF and MRC-5 fibroblasts (73 commonly up- and 70 commonly down-regulated). Thus, on the gene level HFF and MRC-5 cells respond only partially similar to irradiation, to a large extent the cellular response is cell strain specific. The heatmap comparison of the common most differentially regulated genes during both, replicative and irradiation induced senescence in both MRC-5 and HFF illustrates this point (Additional 2: Figure S2). In contrast, for their transition into replicative senescence we found a strong common gene regulation between HFF and MRC-5  and among five human primary fibroblast strains (78 %) . The strain specific response to irradiation is further supported by our observation that among the fifty most differentially regulated genes in MRC-5 and HFF strains, one commonly regulated gene was found, TGFB2. TGFB2 is involved in the regulation of immune privilege, proliferation, differentiation and adhesion . Furthermore, TGFB2 is associated with senescence  and was found, as in these irradiated cells, significantly up-regulated in five replicatively senescent fibroblast strains including MRC-5 and HFF .
Identical markers involved in replicative senescence and premature senescence induced by γ-irradiation
Recently, applying the same experimental procedure, we revealed the most significantly differentially expressed common genes during replicative senescence in HFF and MRC-5 fibroblasts . As the next step, we compared these data with the irradiation results obtained here, by applying the stringency criteria of p < 0.001, and adherence to both statistical packages (edgeR and DESeq). For HFFs we found a total of 2589 commonly significantly differentially regulated genes in both, replication and irradiation induced senescent fibroblasts compared to controls. 2192 of these genes (85 %) were either up- or down-regulated in the same direction while the remaining 15 % were up-regulated in the one case but down in the other. Correspondingly, for MRC-5 we found a total of 936 commonly significantly differentially regulated genes in both, replication and irradiation induced senescent fibroblasts compared to controls. 689 of these genes (74 %) were either up- or down-regulated in the same direction. We thus found that for both fibroblast strains, the transition into replicative as well as irradiation induced senescence correlated with the common differential expression of a large number of genes and with a high degree of similarity in this common differential gene regulation. Interestingly, this common behaviour was observed for a considerably larger number of genes with a higher degree of similarity in HFF than in MRC-5. Our general conclusion is consistent with a recent study using human female lung diploid IMR-90 fibroblast strains . Using Affymetrix arrays, this study compared RNA levels of 5 Gy γ-irradiation induced with replicatively senescent IMR-90 fibroblasts and found a number of genes differentially regulated in cells either arrested by irradiation or replicative exhaustion, with a strong overlap among regulated genes or showing a general trend in the same direction . These data demonstrate the similarities in differential gene regulation between the two types of senescence induction and suggest that the majority of expression changes in replicatively senescent cells were due to proliferation arrest.
In HFF, among the most significant DEG in replication and irradiation-induced senescence were the genes EGR1, Ki67, CCNB1, CCNA2, Id3, Id1, CLDN1, LIF, FBL, CST3, GRN and TMEM47. Similarly, in MRC-5 these were EGR1, Ki67, CCNB1, CCNA2, Id3, Id1, CLDN11, LIF, FBL, CTSK, MMP3 and Wnt16 (stringency criteria: p < 0.001 and adherence to statistical packages, edgeR and DESeq). A number of these genes have cell cycle functions. GRN proteins play a role in wound healing . Ki67 is a marker for proliferation [35, 87]. CTSK is normally stimulated by inflammatory cytokines released after tissue injury . CST3 has been associated with aging-related loss of skeletal muscle (“sarcopenia”) . Down-regulation of Id1 and Id3, as observed in our primary fibroblast strains, has been detected previously in BJ foreskin, WS1 fetal skin, and LF1 lung human fibroblasts . Furthermore, Id loses function in cells transiting into senescence [91, 92]. CCNA2 is down-regulated in aging IMR-90 and WI-38 fibroblasts . Expression of CCNB1 decreases due to antibiotic treatment, resulting in senescence induction in several cell types [94–96]. Reduced expression of CCNB1 inhibits the proliferation of breast cancer cells . MMP3 up-regulation, as seen in our fibroblast strains, is reminiscent to their up-regulation during senescence in human melanocytes [98, 99]. Wnt16 is associated with senescence . Thus, these genes are associated with proliferation, cell cycle arrest or senescence. We found here that these genes commonly correlate with senescence, independent of being irradiation induced or due to replicative exhaustion. Potentially, they are functionally involved in induction of senescence. As found for IL-6 and IL-8 [50, 51], we observed a significant increase in secretion of IL-11 in the media of HFF and MRC-5 fibroblasts undergoing replicative senescence compared to young control fibroblasts (data not shown). Here, we found the mRNA expression levels of the senescence associated secretory phenotype (SASP) family members GRN, CTSK, CST3, MMP3 and IGFBP7/5/3 up-regulated in irradiation induced senescent fibroblasts. These results are consistent with the establishment of a SASP .
Several of the above genes (Ki67, CCNB1, CCNA2, LIF, FBL, CLDN1, WNT16, IGFBP3 and IGFBP7) were also among the commonly significantly differentially regulated genes during replicative and irradiation induced senescence in IMR-90 fibroblasts . However, some of the significantly differentially regulated genes retrieved in  were not identified in our study. This difference could be attributed to the difference in (1) the fibroblast strain (IMR-90 strain  compared to HFF and MRC-5 in our study), (2) the technique used to retrieve the differentially expressed genes (Affymetrix arrays compared to RNA-seq in this study), (3) differences in stringency criteria of p < 0.05 and fold change >2 in  compared to p < 0.001 and adherence to both statistical packages (edgeR and DESeq) in the present study, and finally (4) the intensity of the Gy irradiation (5 Gy  compared to 20 Gy in our study).
Impact of irradiation induced senescence on major transcription factors involved in cell survival
The transcription factors FOXM1 and E2F1 play an important role in cell survival [101–109]. As in both, MRC-5 and HFF fibroblast strains undergoing replicative senescence , FOXM1 and E2F1 were found here to be significantly (log2 fold change >1) down-regulated in irradiation induced senescent fibroblasts.
The down-regulation of FOXM1 explains the significant down-regulation of the cell cycle associated genes CENPF and CCNB2 [101, 102] in irradiation induced senescent cells. FOXM1 has been revealed to have a positive feedback loop with Polo-like-kinase 1 (Plk1) and a negative feedback loop with p53 . Furthermore, FOXM1 has been functionally associated to the expression of X-ray cross-complementing group 1 (XRCC1) involved in base excision repair and breast cancer-associated gene 2 (BRCA2) dealing with homologous recombination repair of DNA double-strand breaks . Similar to fibroblast strains undergoing replicative senescence , Plk1 mRNA expression levels were significantly down-regulated, parallel to FOXM1, in irradiated fibroblasts whereas the expression levels of p53, XRCC1 and BRCA2 were not significantly different to controls.
E2F1 is associated with senescence and cell cycle function [112, 113]. Its downstream targets p14, MMP1 and MMP3 were found here not to be significantly differentially regulated in irradiation induced senescent cells. Expression levels of other transcription factors including ATF1 [114, 115], CREB1 , NFκB1  and HSF1 [118, 119] revealed no significant differential expression on irradiation induced senescence. None of the five members of the NFκB family (NFκB1, NFκB2, RelA, RelB, c-Rel) were significantly differentially regulated in irradiation induced senescent cells. The lack of differential expression of E2F1, ATF1, CREB1, NFκB1 and HSF1 was also observed in senescence induced by replicative exhaustion. The significant differential regulation of FOXM1, E2F1, Plk1 and CENPF were also observed in the previous study  conducted in IMR-90 strains.
Interestingly, the mRNA expression levels of cyclin dependent kinase inhibitors (CDKIs) associated with senescence induction [4, 15, 27, 31, 120] were not among the significantly differentially regulated genes in irradiation induced senescent cells compared to controls in both, MRC-5 and HFF fibroblasts. However, protein expression levels of p21 and p16 were found to be significantly up-regulated in irradiation induced senescent fibroblasts compared to controls (Fig. 4). In fact, the mRNA expression level of CDKN2A (p16) was significantly down-regulated in HFF strains (Fig. 5). The selective lack of correlation of mRNA and protein expression levels has been observed before by [121–123] and by us in MRC-5 and HFF strains undergoing replicative senescence . Thus, our results reveal that protein expression of p16 and p21 is regulated by other down-stream mechanisms than at transcriptional levels.
In contrast, protein expression levels of other selected markers associated with senescence in primary human fibroblast strains (Id3, IGFBP3, IGFBP5 and IGFBP7) revealed a good correlation with mRNA expression levels (Fig. 4). The mRNA and protein expression levels of all three IGFBP family members were significantly up-regulated in both HFF and MRC-5 senescent fibroblast strains, consistent with earlier observations [124–129]. A full proteomics analysis of the transition into senescence of human primary fibroblast strains is under way in our laboratory (to be published elsewhere).
Retrieval of KEGG pathways significantly differentially regulated on irradiation induced senescence
Next, we retrieved the functional pathways significantly (p < 0.05) up- or down-regulated in irradiation induced senescent primary fibroblast strains.
After irradiation, we did not observe either an induction of BAX or a cleavage of caspase-3 (see above in “Gamma irradiation resulted in premature senescence induction in primary human fibroblast strains” section ), indicating that apoptosis was not induced. Analyzing the expression of genes involved in the “Apoptosis” KEGG pathway confirmed this finding: caspase gene family members and other genes having a role in apoptosis induction, including BAX, were not significantly up-regulated after irradiation in either of the two fibroblast strains.
We compared irradiation induced pathways with those found in replicatively senescent cells of the same two fibroblast strains [49, 70]. In HFF strains, six KEGG pathways, namely “arrhythmogenic right ventricular cardiomyopathy”, “cell adhesion molecules”, “dilated cardiomyopathy”, “ECM receptor interaction”, “PPAR signaling pathway” and “long term depression”, were found significantly up-regulated during replicative as well as irradiation-induced senescence (Additional file 3: Figure S3, Additional file 4: Figure S4, Additional file 5: Figure S5, Additional file 6: Figure S6, Additional file 7: Figure S7, Additional file 8: Figure S8). Only the “cell cycle” pathway was commonly down-regulated during senescence induced by both means (Fig. 5). In MRC-5 fibroblast strains, “NOD like receptor signaling pathway”, “cell cycle” and “TGF Beta signaling pathway” were commonly down-regulated in both, irradiation-induced and replicative senescence (Fig. 6, Additional file 11: Figures S9, Additional file 12: S10). Thus, “cell cycle” was the only pathway similarly significantly differentially down-regulated across all four cases (in replicative and irradiation-induced senescence of both fibroblast strains). This finding is consistent with earlier results obtained for IMR-90 fibroblasts  and confirms the tight link of the control of DNA repair to cell cycle regulation.
SASPs are an important hallmark and functional mediator of senescent cells . Unexpectedly, a number of cytokines and cytokine receptors (IL11, EGFR, CXCL-1,2,3,5,6,14) were significantly down-regulated on irradiation induced senescence in MRC-5 and HFF strains, resulting in a significant down-regulation of the KEGG pathway “cytokine–cytokine receptor interaction” (hsa04060) representing SASP. In contrast, TGFB2 was significantly up-regulated more than fivefold. Our results point to a heterogeneous regulation of SASP at the transcript level. Measuring protein levels directly, for example by antibody arrays , may provide a clearer picture of irradiation induced SASP regulation. Of further interest is the significant down-regulation of TGF-beta signaling pathway in both HFF and MRC-5 (Additional file 3: Figure S3, Additional file 9) strains.
The observed difference in significantly differentially regulated pathways on irradiation induced senescence between the two fibroblast strains is not due to experimental error since the strain triplicates cluster closely together (Fig. 3). The difference could be attributed to the strain differences in origin (MRC-5, embryonic lung; HFF, foreskin). Furthermore, their difference in PD numbers could also contribute to this difference: In our experiments, the MRC-5 cells (ordered from ATCC) had the starting PD 28 while we received HFF cells isolated from foreskin of primary human donors at PD 12.
While we found a strong similarity in the differential expression of genes for both senescence induction processes, we identified a difference on the level of functional pathways. Irradiation induced damage activates cellular repair processes [130, 131], often combined with p53, p21 and p16 mediated cell cycle arrest and, if repair is not successful, transition into senescence [1, 15, 132–135]. In irradiation induced senescence, only a few genes of the repair pathways were significantly down-regulated below control levels in MRC-5 fibroblasts (Additional file 10). In HFF strains we observed a down-regulation of all genes involved in the three DNA repair pathways, however, to a lesser extent than the down-regulation of other pathways. In contrast, in replicatively senescent cells all three repair pathways were significantly down-regulated in both fibroblast cell strains [49, 70]. Since we analyzed cells in early senescence, this quantitative difference could potentially indicate that replicatively senescent cells shut down repair pathways earlier than irradiation induced senescent cells during their transition into senescence.
Similarly, only a few genes of the “replication” pathway were found significantly down-regulated in irradiation induced senescent cells of both strains. Instead, “replication” was the pathway with nearly all genes significantly down-regulated in replicatively senescent fibroblast cells. Thus, during the transition into replicative senescence two essential functions, DNA repair and replication are more stringently regulated than during the transition into irradiation induced senescence. This is consistent with the view that replication errors are essential for the induction of replicative senescence while this process is not as relevant for irradiation induced senescence. Consequently, in replicative senescence the “replication” pathway is completely down-regulated.
Next, we analyzed the genes of significantly differentially regulated pathways (up- or down-regulated in either of the fibroblast strains) (Figs. 5, 6; Additional file 3: Figure S3, Additional file 4: Figure S4, Additional file 5: Figure S5, Additional file 6: Figure S6, Additional file 7: Figure S7, Additional file 8: Figure S8, Additional file 11: Figure S9, Additional file 12: Figure S10). We retrieved the expression levels of the involved genes (both up- and down-regulated) in the two fibroblast strains and identified which genes were commonly differentially regulated in both replicative and irradiation induced senescence (Figs. 5, 6; Additional file 3: Figure S3, Additional file 4: Figure S4, Additional file 5: Figure S5, Additional file 6: Figure S6, Additional file7: Figure S7, Additional file 8: Figure S8, Additional file 11: Figure S9, Additional file 12: Figure S10). The comparison with previous studies enabled us to functionally associate a number of these genes with induction of cell cycle arrest and senescence (highlighted in blue in Additional file 2). Several genes having a role in senescence induction (such as TGFB2, IGF1, Id1, Id3, Id4, IL1B, IL6 and IL8) were among those genes (highlighted in blue in Additional file 10) which were similarly differentially regulated in both replicative and irradiation induced senescence [50, 55, 84, 90–92, 136–143]. Most importantly, the list (Additional file: 10) also includes on the one hand, genes which have not previously been associated with induction of senescence and, on the other hand, genes which are differentially regulated exclusively during irradiation induced senescence. In future studies, we intend to functionally validate the role of several of these genes in senescence induction by irradiation, replicative exhaustion or both.
We compared the transcriptomes of two young and senescent human primary fibroblast strains, with the senescent state either induced by γ-irradiation or by replicative exhaustion. We found a strong similarity in the differential expression of genes for both senescence induction processes, indicating a considerably common cellular response to either internal or external damage. On the functional pathway level, “Cell cycle” was the only pathway commonly (down-) regulated in replicative and irradiation-induced senescence in both fibroblast strains, confirming the tight link between DNA repair and cell cycle regulation. In γ-irradiation induced senescence, only a few genes of the repair pathways were significantly down-regulated below control levels in MRC-5 strains. In HFF strains we observed a down-regulation of all genes involved in DNA repair pathways, however, to an extent less significant than the down-regulation of other pathways. In contrast, all three repair pathways are significantly down-regulated in replicatively senescent fibroblast cells. Furthermore, only a few genes of the “replication” pathway were found significantly down-regulated in irradiation induced senescent cells. Instead, “replication” was the pathway with nearly all genes significantly down-regulated in replicatively senescent fibroblasts. Thus, on the pathway level we identified considerable differences between both senescent states. During the transition into replicative senescence two essential functions, DNA repair and replication are more stringently regulated than during the transition into irradiation induced senescence, consistent with replication errors being essential for the induction of replicative senescence while this process is not as relevant for irradiation induced senescence.
human foreskin fibroblasts
Dulbeccos modified Eagles low glucose medium
fetal bovine serum
- CO2 :
principle component analysis
reads per kilo base per million mapped reads
false discovery rate
- SA β-Gal:
senescence associated β-Gal
high-throughput RNA sequencing
cyclin dependent kinase inhibitors
gene set enrichment for pathway analysis
Kyoto encyclopedia of genes and genomes
phosphorylated retinoblastoma protein
phosphate buffered saline
RNA integrity number
differentially expressed genes
transmembrane protein 47
insulin-like growth factor-binding protein 3
insulin-like growth factor-binding protein 5
insulin-like growth factor-binding protein 7
cyclin-dependent kinase inhibitor 2A
forkhead Box M1
activating transcription factor 1
CAMP responsive element binding protein 1
heat shock transcription factor 1
NF kappa B signaling
transforming growth factor beta 2
early growth response 1
leukemia inhibitory factor
dNA binding protein inhibitor
- K3Fe(CN)6 :
- K4Fe(CN)6 :
forkhead box protein M1
transcription factor E2F1
X-ray repair cross-complementing protein 1
breast cancer type 2 susceptibility protein
peroxisome proliferator-activated receptors
nucleotide-binding oligomerization domain receptors
interleukin 1 beta
Mirzayans R, Andrais B, Scott A, Wang YW, Murray D. Ionizing radiation-induced responses in human cells with differing TP53 status. Int J Mol Sci. 2013;14:22409–35.
Di Leonardo A, Linke SP, Clarkin K, Wahl GM. DNA damage triggers a prolonged p53-dependent G1 arrest and long-term induction of Cip1 in normal human fibroblasts. Genes Dev. 1994;8:2540–51.
Serrano M, Lin AW, McCurrach ME, Beach D, Lowe SW. Oncogenic Ras provokes premature cell senescence associated with accumulation of p53 and p16INK4a. Cell. 1997;88(5):593–602.
Robles SJ, Adami GR. Agents that cause DNA double strand breaks lead to p16INK4a enrichment and the premature senescence of normal fibroblasts. Oncogene. 1998;16(9):1113–23.
Toussaint O, Medrano EE, von Zglinicki T. Cellular and molecular mechanisms of stress-induced premature senescence (SIPS) of human diploid fibroblasts and melanocytes. Exp Gerontol. 2000;35(8):927–45.
Hwang ES. Replicative senescence and senescence-like state induced in cancer-derived cells. Mech Ageing Dev. 2002;123(12):1681–94.
Zhao J, Liu XJ, Ma JW, Zheng RL. DNA damage in healthy term neonate. Early Hum Dev. 2004;77(1–2):89–98.
Debacq-Chainiaux F, Borlon C, Pascal T, Royer V, Eliaers F, Ninane N, et al. Repeated exposure of human skin fibroblasts to UVB at sub-cytotoxic level triggers premature senescence through the TGF-β1 signaling pathway. J Cell Sci. 2005;118(4):743–58.
Havelka AM, Berndtsson M, Olofsson MH, Shoshan MC, Linder S. Mechanisms of action of DNA-damaging anticancer drugs in treatment of carcinomas: is acute apoptosis an “off-target” effect? Mini Rev Med Chem. 2007;7(10):1035–9.
Mallette FA, Gaumont-Leclerc MF, Ferbeyre G. The DNA damage signaling pathway is a critical mediator of oncogene-induced senescence. Genes Dev. 2007;21(1):43–8.
Maslov AY, Vijg J. Genome instability, cancer and aging. Biochim Biophys Acta. 2009;1790(10):963–9.
Haupt S, Berger M, Goldberg Z, Haupt Y. Apoptosis—the p53 network. J Cell Sci. 2003;116:4077–85.
Zhivotovsky B, Kroemer G. Apoptosis and genomic instability. Nat Rev Mol Cell Biol. 2004;5:752–62.
Blank M, Shiloh Y. Programs for cell death: apoptosis is only one way to go. Cell Cycle. 2007;6:686–95.
Ben-Porath I, Weinberg RA. The signals and pathways activating cellular senescence. Int J Biochem Cell Biol. 2005;37(5):961–76.
Campisi J, Sedivy J. How does proliferative homeostasis change with age? What causes it and how does it contribute to aging? J Gerontol A Biol Sci Med Sci. 2009;64(2):164–6.
Hayflick L. The limited in vitro lifetime of human diploid cell strains. Exp Cell Res. 1965;37:614–36.
Berube NG, Smith JR, Pereira-Smith OM. The genetics of cellular senescence. Am J Hum Genet. 1998;62:1015–9.
Smith JR, Pereira-Smith OM. Replicative senescence: implications for in vivo aging and tumor suppression. Science. 1996;273:63–6.
Campisi J, Dimri G, Hara E. Control of replicative senescence. In: Schneider E, Rowe J, editors. Handbook of the biology of aging. San Diego: Academic Press; 1996. p. 121–40.
Smogorzewska A, de Lange T. Different telomere damage signaling pathways in human and mouse cells. EMBO J. 2002;21(16):4338–48.
d’Adda di Fagagna F, Reaper PM, Clay-Farrace L, Fiegler H, Carr P, Von Zglinicki T, et al. A DNA damage checkpoint response in telomere-initiated senescence. Nature. 2003;426(6963):194–8.
Takai H, Smogorzewska A, de Lange T. DNA damage foci at dysfunctional telomeres. Curr Biol. 2003;13(17):1549–56.
Itahana K, Zou Y, Itahana Y, Martinez JL, Beausejour C, Jacobs JJ, et al. Control of the replicative life span of human fibroblasts by p16 and the polycomb protein Bmi-1. Mol Cell Biol. 2003;23(1):389–401.
Ohtani N, Yamakoshi K, Takahashi A, Hara E. The p16INK4a-RB pathway: molecular link between cellular senescence and tumor suppression. J Med Invest. 2004;51:146–53.
Gire V, Roux P, Wynford-Thomas D, Brondello JM, Dulic V. DNA damage checkpoint kinase Chk2 triggers replicative senescence. EMBO J. 2004;23(13):2554–63.
Herbig U, Jobling WA, Chen BPC, Chen DJ, Sedivy JM. Telomere Shortening Triggers Senescence of Human Cells through a Pathway Involving ATM, p53, and p21CIP1, but not p16INK4a. Mol Cell. 2004;14(4):501–13.
Ben-Porath I, Weinberg RA. When cells get stressed: an integrative view of cellular senescence. J Clin Invest. 2004;113(1):8–13.
Shay JW, Wright WE. Senescence and immortalization: role of telomeres and telomerase. Carcinogenesis. 2005;26(5):867–74.
Jeyapalan JC, Ferreira M, Sedivy JM, Herbig U. Accumulation of senescent cells in mitotic tissue of aging primates. Mech Ageing Dev. 2007;128(1):36–44.
Goodarzi AA, Noon AT, Deckbar D, Ziv Y, Shiloh Y, Löbrich M, et al. ATM signaling facilitates repair of DNA double-strand breaks associated with heterochromatin. Mol Cell. 2008;31(2):167–77.
Cosme-Blanco W, Chang S. Dual roles of telomere dysfunction in initiation and suppression of tumorigenesis. Exp Cell Res. 2008;314(9):1973–9.
Nakamura AJ, Redon CE, Bonner WM, Sedelnikova OA. Telomere-dependent and telomere-independent origins of endogenous DNA damage in tumor cells. Aging. 2009;1(2):212–8.
Kuilman TC, Michaloglou C, Mooi WJ, Peeper DS. The essence of senescence. Genes Dev. 2010;24(22):2463–79.
Marthandan S, Priebe S, Hemmerich P, Klement K, Diekmann S. Long-term quiescent fibroblast cells transit into senescence. PLoS One. 2014;9:e115597.
Xie Z, Jay KA, Smith DL, Zhang Y, Liu Z, Zheng J, et al. Early telomerase inactivation accelerates aging independently of telomere length. Cell. 2015;160(5):928–39.
Roninson IB. Tumor Cell Senescence in Cancer Treatment. Cancer Res. 2003;63(11):2705–15.
Campisi J. Cellular senescence: putting the paradoxes in perspective. Curr Opin Genet Dev. 2011;21:107–12.
Narita M, Lowe SW. Senescence comes of age. Nat Med. 2005;11:920–2.
Braig M, Lee S, Loddenkemper C, Rudolph C, Peters AH, Schlegelberger B, et al. Oncogene-induced senescence as an initial barrier in lymphoma development. Nature. 2005;436:660–5.
Chen Z, Trotman LC, Shaffer D, Lin H, Dotan ZA, Niki M, et al. Critical role of p53 dependent cellular senescence in suppression of Pten deficient tumourigenesis. Nature. 2005;436:725–30.
Collado M, Gil J, Efeyan A, Guerra C, Schumacher AJ, Barradas M, et al. Tumor biology: senescence in premalignant tumours. Nature. 2005;436:642.
Michaloglou C, Vredeveld LCW, Soengas MS, Denoyelle C, Kuilman T, van der Horst CM, et al. BRAFE600- associated senescence-like cell cycle arrest of human nevi. Nature. 2005;436:720–4.
Courtois-Cox S, Genther Williams SM, Reczek EE, Johnson BW, McGillicuddy LT, Johannessen CM, et al. A negative feedback signaling network underlies oncogene-induced senescence. Cancer Cell. 2006;10:459–72.
Ventura A, Kirsch DG, McLaughlin ME, Tuveson DA, Grimm J, Lintault L, et al. Restoration of p53 function leads to tumor regression in vivo. Nature. 2007;445:661–5.
Xue W, Zender L, Miething C, Dickins RA, Hernando E, Krizhanovsky V, et al. Senescence and tumor clearance is triggered by p53 restoration in murine liver carcinomas. Nature. 2007;445:656–60.
Campisi J, d’ Adda di Fagagna F. Cellular senescence: when bad things happen to good cells. Nat Rev Mol Cell Biol. 2007;8:729–40.
Dimri GP, Lee X, Basile G, Acosta M, Scott G, Roskelley C, et al. A biomarker that identifies senescent human cells in culture and in aging skin in vivo. Proc Natl Acad Sci USA. 1995;92(20):9363–7.
Marthandan S, Priebe S, Baumgart M, Groth M, Cellerino A, Guthke R, et al. Similarities in gene expression profiles during in vitro aging of primary human embryonic lung and foreskin fibroblasts. Biomed Res Int. 2015;2015:731938.
Coppé JP, Patil CK, Rodier F, Sun Y, Muñoz DP, Goldstein J, et al. Senescence-associated secretory phenotypes reveal cell-nonautonomous functions of oncogenic RAS and the p53 tumor suppressor. PLoS Biol. 2008;6(12):2853–68.
Rodier F, Coppé JP, Patil CK, Hoeijmakers WAM, Munoz D, Raza SR, et al. Persistent DNA damage signaling triggers senescence-associated inflammatory cytokine secretion. Nat Cell Biol. 2009;11:973–9.
Sun P, Yoshizuka N, New L, Moser BA, Li Y, Liao R, et al. PRAK is essential for ras-induced senescence and tumor suppression. Cell. 2007;128:295–308.
Mathon NF, Lloyd AC. Cell senescence and cancer. Nat Rev Cancer. 2001;1:203–13.
Shelton DN, Chang E, Whittier PS, Choi D, Funk WD. Microarray analysis of replicative senescence. Curr Biol. 1999;9(17):939–45.
Rodier F, Campisi J. Four faces of cellular senescence. J Cell Biol. 2011;192(4):547–56.
Schäuble S, Klement K, Marthandan S, Münch S, Heiland I, Schuster S, et al. Quantitative model of cell cycle arrest and cellular senescence in primary human fibroblasts. PLoS ONE. 2012;7:e42150.
Kronschnabl M, Stamminger T. Synergistic induction of intercellular adhesion molecule-1 by the human cytomegalovirus transactivators IE2p86 and pp71 is mediated via an Sp1-binding site. J Gen Virol. 2003;84(1):61–73.
Honda S, Hjelmeland LM, Handa JT. Oxidative stress–induced single-strand breaks in chromosomal telomeres of human retinal pigment epithelial cells in vitro. Invest Ophthalmol Vis Sci. 2001;42(9):2139–44.
Sivakumar S, Daum JR, Tipton AR, Rankin S, Gorbsky GJ. The spindle and kinetochore associated (Ska) complex enhances binding of the anaphase promoting complex/cyclosome (APC/C) to chromosomes and promote s mitotic exit. Mol Cell Biol. 2014;25(5):594–605.
Bentley DR, Balasubramanian S, Swerdlow HP, Smith GP, Milton J, Brown CG, et al. Accurate whole human genome sequencing using reversible terminator chemistry. Nature. 2008;456(7218):53–9.
Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11.
Flicek P, Amode MR, Barrell D, Beal K, Brent S, Carvalho-Silva D, et al. Ensembl 2012. Nucleic Acids Res. 2012;40:D84–90.
R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria; 2008. ISBN 3-900051-07-0, URL http://www.R-project.org.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B. 1995;57(1):289–300.
Rapaport F, Khanin R, Liang Y, Pirun M, Krek A, Zumbo P, et al. Comprehensive evaluation of differential gene expression analysis methods for RNA-seq data. Genome Biol. 2013;14(9):R95.
Seyednasrollah F, Laiho A, Elo L. Comparison of software packages for detecting differential expression in RNA-seq studies. Brief Bioinform. 2015;16(1):59–70.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Marthandan S, Baumgart M, Priebe S, Groth M, Schaer J, Kaether C, et al. Conserved senescence associated genes and pathways in primary human fibroblasts detected by RNA-seq. PLoS One. 2016;11(5):e0154531.
Marthandan S, Priebe S, Groth M, Guthke R, Platzer M, Hemmerich P, et al. Hormetic effect of rotenone in primary human fibroblasts. Immun Ageing. 2015;12:11.
Luo W, Friedman MS, Shedden K, Hankenson KD, Woolf PJ. GAGE: generally applicable gene set enrichment for pathway analysis. BMC Bioinform. 2009;10:161.
Chicas A, Wang X, Zhang C, McCurrach M, Zhao Z, Mert O, et al. Dissecting the unique role of the retinoblastoma tumor suppressor during cellular senescence. Cancer Cell. 2010;17(4):376–87.
Aksoy O, Chicas A, Zeng T, Zhao Z, McCurrach M, Wang X, et al. The atypical E2F family member E2F7 couples the p53 and RB pathways during cellular senescence. Genes Dev. 2012;26:1546–57.
Chang BD, Broude EV, Dokmanovic M, Zhu H, Ruth A, Xuan Y, et al. A senescence-like phenotype distinguishes tumor cells that undergo terminal proliferation arrest after exposure to anticancer agents. Cancer Res. 1999;59(15):3761–7.
Purcell M, Kruger A, Tainsky MA. Gene expression profiling of replication and induced senescence. Cell Cycle. 2014;13(24):3927–37.
Passos JF, Nelson G, Wang C, Richter T, Simillion C, Proctor CJ, et al. Feedback between p21 and reactive oxygen production is necessary for cell senescence. Mol Syst Biol. 2010;6:347.
Finucane DM, Bossy-Wetzel E, Waterhouse NJ, Cotter TG, Green DR. Bax-induced caspase activation and apoptosis via cytochrome c release from mitochondria in inhibitable by Bcl-xL. J Biol Chem. 1999;274(4):2225–33.
Pawlowski J, Kraft AS. Bax-induced apoptotic cell death. Proc Natl Acad Sci USA. 2000;97(2):529–31.
Parrish AB, Freel CD, Kornbluth S. Cellular mechanisms controlling caspase activation and function. Cold Spring Harb Perspect Biol. 2013;5(6):a008672.
Linke SP, Clarkin KC, Wahl GM. p53 mediates permanent arrest over multiple cell cycles in response to gamma-irradiation. Cancer Res. 1997;57(6):1171–9.
Muthna D, Soukup T, Vavrova J, Mokry J, Cmielova J, Visek B, et al. Irradiation of adult human dental pulp stem cells provokes activation of p53, cell cycle arrest and senescence but not apoptosis. Stem Cells Dev. 2010;19(12):1855–62.
Hirsch L, Nazari H, Sreekumar PG, Kannan R, Dustin L, Zhu D, et al. TGF-β2 secretion from RPE decreases with polarization and becomes apically oriented. Cytokine. 2015;71(2):394–6.
Lin S, Yang J, Elkahloun AG, Bandyopadhyay A, Wang L, Cornell JE, et al. Attenuation of TGF-β signaling suppresses premature senescence in a p21-dependent manner and promotes oncogenic Ras-mediated metastatic transformation in human mammary epithelial cells. Mol Biol Cell. 2012;23(8):1569–81.
Lackner DH, Hayashi MT, Cesare AJ, Karlseder J. A genomics approach identifies senescence-specific gene expression regulation. Aging Cell. 2014;13(5):946–50.
Smout MJ, Laha T, Mulvenna J, Sripa B, Suttiprapa S, Jones A, et al. A granulin-like growth factor secreted by the carcinogenic liver fluke, Opisthorchis viverrini, promotes proliferation of host cells. PLoS Pathog. 2009;5(10):e1000611.
Petric M, Martinez S, Acevedo F, Oddo D, Artigas R, Camus M, et al. Correlation between Ki67 and histological grade in breast cancer patients treated with preoperative chemotherapy. Asian Pac J Cancer Prev. 2014;15(23):10277–80.
Cheng XW, Kikuchi R, Ishii H, Yoshikawa D, Hu L, Takahashi R, et al. Circulating cathepsin K as a potential novel biomarker of coronary artery disease. Atherosclerosis. 2013;228(1):211–6.
de Winter CF, Echteld MA, Evenhuis HM. Chronic kidney disease in older people with intellectual disability: results of the HA-ID study. Res Dev Disabil. 2014;35(3):726–32.
Kong Y, Cui H, Zhang H. Smurf2-mediated ubiquitination and degradation of Id1 regulates p16 expression during senescence. Aging Cell. 2011;10(6):1038–46.
Hara E, Yamaguchi T, Nojima H, Ide T, Campisi J, Okayama H, et al. Id-related genes encoding helix-loop-helix proteins are required for G1 progression and are repressed in senescent human fibroblasts. J Biol Chem. 1994;269(3):2139–45.
Ling F, Kang B, Sun XH. Id proteins: small molecules, mighty regulators. Curr Top Dev Biol. 2014;110:189–216.
Chen T, Xue L, Niu J, Ma L, Li N, Cao X, et al. The retinoblastoma protein selectively represses E2F1 targets via a TAAC DNA element during cellular senescence. J Biol Chem. 2012;287(44):37540–51.
Kikuchi N, Ishii Y, Morishima Y, Yageta Y, Haraguchi N, Itoh K, et al. Nrf2 protects against pulmonary fibrosis by regulating the lung oxidant level and Th1/Th2 balance. Respir Res. 2010;11:31.
Zuryń A, Litwiniec A, Gackowska L, Pawlik A, Grzanka AA, Grzanka A. Expression of cyclin A, B1 and D1 after induction of cell cycle arrest in the Jurkat cell line exposed to doxorubicin. Cell Biol Int. 2012;36(12):1129–35.
Zuryń A, Gagat M, Grzanka AA, Gackowska L, Grzanka A. Expression of cyclin B1 after induction of senescence and cell death in non-small cell lung carcinoma A549 cells. Folia Histochem Cytobiol. 2012;50(1):58–67.
Androic I, Krämer A, Yan R, Rödel F, Gätje R, Kaufmann M, et al. Targeting cyclin B1 inhibits proliferation and sensitizes breast cancer cells to taxol. BMC Cancer. 2008;8:391.
Millis AJT, Hoyle M, McCue HM, Martini H. Differential expression of metalloproteinase and tissue inhibitor of metalloproteinase genes in aged human fibroblasts. Exp Cell Res. 1992;201(2):373–9.
Kang MK, Kameta A, Shin KH, Baluda MA, Kim HR, Park NH. Senescence-associated genes in normal human oral keratinocytes. Exp Cell Res. 2003;287(2):272–81.
Binet R, Ythier D, Robles AI, Collado M, Larrieu D, Fonti C, et al. Wnt16b is a new marker of cellular senescence that regulates p53activity and the phosphoinositide 3-Kinase/AKT pathway. Cancer Res. 2009;69(24):9183–91.
Wierstra I, Alves J. FOXM1, a typical proliferation-associated transcription factor. Biol Chem. 2007;388(12):1257–74.
Laoukili J, Kooistra MR, Brás A, Kauw J, Kerkhoven RM, Morrison A, et al. FoxM1 is required for execution of the mitotic programme and chromosome stability. Nat Cell Biol. 2005;7(2):126–36.
Myatt SS, Lam EW. The emerging roles of forkhead box (Fox) proteins in cancer. Nat Rev Cancer. 2007;7(11):847–59.
Kalinichenko VV, Major ML, Wang X, Petrovic V, Kuechle J, Yoder HM, et al. Foxm1b transcription factor is essential for development of hepatocellular carcinomas and is negatively regulated by the p19ARF tumor suppressor. Genes Dev. 2004;18(7):830–50.
Wonsey DR, Follettie MT. Loss of the forkhead transcription factor FoxM1 causes centrosome amplification and mitotic catastrophe. Cancer Res. 2005;65(12):5181–9.
Kim IM, Ackerson T, Ramakrishna S, Tretiakova M, Wang IC, Kalin TV, et al. The Forkhead Box m1 transcription factor stimulates the proliferation of tumor cells during development of lung cancer. Cancer Res. 2006;66(4):2153–61.
Kalin TV, Wang IC, Ackerson TJ, Major ML, Detrisac CJ, Kalinichenko VV, et al. Increased levels of the FoxM1 transcription factor accelerate development and progression of prostate carcinomas in both TRAMP and LADY transgenic mice. Cancer Res. 2006;66(3):1712–20.
Chan DW, Yu SY, Chiu PM, Yao KM, Liu VW, Cheung AN, et al. Over-expression of FOXM1 transcription factor is associated with cervical cancer progression and pathogenesis. J Pathol. 2008;215(3):245–52.
Wu QF, Liu C, Tai MH, Liu D, Lei L, Wang RT, et al. Knockdown of FoxM1 by siRNA interference decreases cell proliferation, induces cell cycle arrest and inhibits cell invasion in MHCC-97H cells in vitro. Acta Pharmacol Sin. 2010;31(3):361–6.
Pandit B, Halasi M, Gartel AL. p53 negatively regulates expression of FoxM1. Cell Cycle. 2009;8(20):3425–7.
Tan Y, Raychaudhuri P, Costa RH. Chk2 mediates stabilization of the FoxM1 transcription factor to stimulate expression of DNA repair genes. Mol Cell Biol. 2007;27(3):1007–16.
Dimri GP, Hara E, Campisi J. Regulation of two E2F-related genes in pre-senescent and senescent human fibroblasts. J Biol Chem. 1994;269(23):16180–6.
Dimri GP, Itahana K, Acosta M, Campisi J. Regulation of a senescence checkpoint response by the E2F1 transcription factor and p14(ARF) tumor suppressor. Mol Cell Biol. 2000;20(1):273–85.
Ors A, Grimaldi M, Kimata Y, Wilkinson CRM, Jones N, Yamano H. Protein synthesis, Post-translational modification and degradation. J Biol Chem. 2009;284:23989–94.
Gupta P, Prywes R. ATF1 phosphorylation by the ERK MAPK pathway is required for epidermal growth factor-induced c-jun expression. J Biol Chem. 2002;277:50550–6.
Chin JH, Okazaki M, Frazier JS, Hu ZW, Hoffman BB. Impaired cAMP-mediated gene expression and decreased cAMP response element binding protein in senescent cells. Am J Physiol. 1996;271(1 Pt 1):C362–71.
Vaughan S, Jat PS. Deciphering the role of nuclear factor-jB in cellular senescence. Aging. 2011;3:913–9.
Bonelli MA, Alfieri RR, Poli M, Petronini PG, Borghetti AF. Heat induced proteasomic degradation of HSF1 in sterum-starved human fibroblasts aging in vitro. Exp Cell Res. 2001;267:165–72.
Meng L, Gabai VL, Sherman MY. Heat shock transcription factor HSF1 plays a critical role in HER2-induced cellular transformation and tumorigenesis. Oncogene. 2010;29(37):5204–13.
Ressler S, Bartkova J, Niederegger H, Bartek J, Scharffetter-Kochanek K, Jansen-Dürr P, et al. p16INK4A is a robust in vivo biomarker of cellular aging in human skin. Aging Cell. 2006;5(5):379–89.
Alcorta DA, Xiong Y, Phelps D, Hannon G, Beach D, Barrett JC. Involvement of the cyclin-dependent kinase inhibitor p16 (INK4a) in replicative senescence of normal human fibroblasts. Proc Natl Acad Sci USA. 1996;93(24):13742–7.
Chen JH, Ozanne SE. Deep senescent human fibroblasts show diminished DNA damage foci but retain checkpoint capacity to oxidative stress. FEBS Lett. 2006;580:6669–73.
Kim YM, Byun HO, Jee BA, Cho H, Seo YH, Kim YS, et al. Implications of time-series gene expression profiles of replicative senescence. Aging Cell. 2013;12(4):622–34.
Yoon IK, Kim HK, Kim YK, Song IH, Kim W, Kim S, et al. Exploration of replicative senescence associated genes in human dermal fibroblasts by cDNA microarray technology. Exp Gerontol. 2004;39(9):1369–78.
Baege AC, Disbrow GL, Schlegel R. IGFBP3, a marker of cellular senescence, Is overexpressed in human papillomavirus-immortalized cervical cells and enhances IGF-1-induced mitogenesis. J Virol. 2004;78(11):5720–7.
Kim KS, Seu YB, Baek SH, Kim MJ, Kim KJ, Kim JH, et al. Induction of cellular senescence by insulin-like growth factor binding protein-5 through a p53-dependent mechanism. Mol Biol Cell. 2007;18(11):4543–52.
Johung K, Goodwin EC, DiMaio D. Human papillomavirus E7 repression in cervical carcinoma cells initiates a transcriptional cascade driven by the retinoblastoma family, resulting in senescence. J Virol. 2007;81(5):2102–16.
Dhahbi JM, Atamna H, Boffelli D, Magis W, Spindler SR, Martin DI. Deep sequencing reveals novel microRNAs and regulation of microRNA expression during cell senescence. PLoS ONE. 2011;6(5):e20509.
Wajapeyee N, Serra RW, Zhu X, Mahalingam M, Green MR. Oncogenic BRAF induces senescence and apoptosis through pathways mediated by the secreted protein IGFBP7. Cell. 2008;132(3):363–74.
Ciccia A, Elledge SJ. The DNA damage response: making it safe to play with knives. Mol Cell. 2010;40:179–204.
Larsen DH, Stucki M. Nucleolar responses to DNA double-strand breaks. Nucleic Acids Res. 2015;44(2):538–44.
Suzuki K, Mori I, Nakayama Y, Miyakoda M, Kodama S, Watanabe M. Radiation-induced senescence-like growth arrest requires TP53 function but not telomere shortening. Radiat Res. 2001;155(1 Pt 2):248–53.
Chen JH, Hales CH, Ozanne SE. DNA damage, cellular senescence and organismal ageing: casual or correlative? Nucleic Acids Res. 2007;35(22):7417–28.
Sabin RJ, Anderson RM. Cellular senescence—its role in cancer and the response to ionizing radiation. Genome Integr. 2011;2:7.
Abraham RT. Checkpoint signalling: focusing on 53BP1. Nat Cell Biol. 2002;4(12):277–9.
Kuilman T, Michaloglou C, Vredeveld LC, Douma S, van Doorn R, Desmet CJ, et al. Oncogene-induced senescence relayed by an interleukin-dependent inflammatory network. Cell. 2008;133(6):1019–31.
Bhaumik D, Scott GK, Schokrpur S, Patil CK, Orjalo AV, Rodier F, et al. MicroRNAs miR-146a/b negatively modulate the senescence-associated inflammatory mediators IL-6 and IL-8. Aging. 2009;1(4):402–11.
Bunt J, de Haas TG, Hasselt NE, Zwijnenburg DA, Koster J, Versteeg R, et al. Regulation of cell cycle genes and induction of senescence by overexpression of OTX2 in medulloblastoma cell lines. Mol Cancer Res. 2010;8(10):1344–57.
Hubackova S, Krejcikova K, Bartek J, Hodny Z. IL1- and TGFβ- Nox4 signaling, oxidative stress and DNA damage response are shared features of replicative, oncogene-induced and drug induced paracrine bystander senescence. Aging. 2012;4(12):932–51.
Acosta JC, Banito A, Wuestefeld T, Georgilis A, Janich P, Morton JP, et al. A complex secretory program orchestrated by the inflammasome controls paracrine senescence. Nat Cell Biol. 2013;15(8):978–90.
Carey JP, Knowell AE, Chinaranagari S, Chaudhary J. Id4 promotes senescence and sensitivity to doxorubicin-induced apoptosis in DU145 prostate cancer cells. Anticancer Res. 2013;33(10):4271–8.
Salama R, Sadaie M, Hoare M, Narita M. Cellular senescence and its effector programs. Genes Dev. 2014;28(2):99–114.
Hofmann JW, Zhao X, De Cecco M, Peterson AL, Pagliaroli L, Manivannan J, et al. Reduced expression of MYC increases longevity and enhances healthspan. Cell. 2015;160(3):477–88.
Takahashi Y, Moriwaki S, Sugiyama Y, Endo Y, Yamazaki K, Mori T, et al. Decreased gene expression responsible for post-ultraviolet DNA repair synthesis in aging: a possible mechanism of age-related reduction in DNA repair capacity. J Invest Dermatol. 2005;124(2):435–42.
Kovalchuk IP, Golubov A, Koturbash IV, Kutanzi K, Martin OA, Kovalchuk O. Age dependent changes in DNA repair in radiation exposed mice. Radiat Res. 2014;182(6):683–94.
Kaneko T, Tahara S, Tanno M, Taguchi T. Age-related changes in the induction of DNA polymerases in rat liver by gamma ray irradiation. Mech Ageing Dev. 2002;123(11):1521–8.
Harada H, Nakagawa H, Takaoka M, Lee J, Herlyn M, Diehl JA, et al. Cleavage of MCM2 licensing protein fosters senescence in human keratinocytes. Cell Cycle. 2008;7(22):3534–8.
Paulsen MT, Starks AM, Derheimer FA, Hanasoge S, Li L, Dixon JE, et al. The p53 targeting human phosphatase hCdc14A interacts with the Cdk1/cyclin B complex and is differentially expressed in human cancers. Mol Cancer. 2006;5:25.
Martin N, Raguz S, Dharmalingam G, Gil J. Co-regulation of senescence associated genes by oncogenic homeobox proteins and polycomb repressive complexes. Cell Cycle. 2013;12(14):2194–9.
Zhang Y, Guo L, Xing P, Chen Y, Li F, Zhu W, et al. Increased expression of oncogene induced senescence markers during cervical squamous cell cancer development. Int J Clin Exp Pathol. 2014;7(12):8911–6.
Pollok S, Bauerschmidt C, Sänger J, Nasheuer HP, Grosse F. Human cdc45 is a proliferation associated antigen. FEBS J. 2007;274(14):3669–84.
Kim HJ, Cho JH, Kim JR. Down-regulation of Polo-like kinase 1 induces cellular senescence in human primary cells through a p53 dependent pathway. J Gerontol A Biol Sci Med Sci. 2013;68(10):1145–56.
Hsu YH, Liao LJ, Yu CH, Chiang CP, Jhan JR, Chang LC, et al. Overexpression of the pituitary tumor transforming gene induces p53 dependenr senescence through activating DNA damage response pathway in normal human fibroblasts. J Biol Chem. 2010;285(29):22630–8.
Aliouat-Denis CM, Dendouga N, Van den Wyngaert I, Goehlmann H, Steller U, van de Weyer I, et al. p53 independent regulation of p21 Waf1/Cip1 expression and senescence by Chk2. Mol Cancer Res. 2005;3(11):627–34.
Kidokoro T, Tanikawa C, Furukawa Y, Katagiri T, Nakamura Y, Matsuda K. CDC20, a potential cancer therapeutic target, is negatively regulated by p53. Oncogene. 2008;27(11):1562–71.
Nakayama Y, Yamaguchi N. Role of cyclin B1 levels in DNA damage and DNA damage-induced senescence. Int Rev Cell Mol Biol. 2013;305:303–37.
Kotake Y, Naemura M, Murasaki C, Inoue Y, Okamoto H. Transcriptional regulation of the p16 tumor suppressor gene. Anticancer Res. 2015;35(8):4397–401.
Jia L, Li H, Sun Y. Induction of p21 dependent senescence by an NAE inhibitor, MLN4924, as a mechanism of growth suppression. Neoplasia. 2011;13(6):561–9.
de Jesus Bernardes. B, Blasco MA. Assessing cell and organ senescence biomarkers. Circ Res. 2012;111(1):97–109.
Wang Z, Lin H, Hua F, Hu Z. Repairing DNA damage by XRCC6/KU70 reverses TLR4 deficiency worsened HCC development via restoring senescence and autophagic flux. Autophagy. 2013;9(6):925–7.
Poehlmann A, Habold C, Walluscheck D, Reissig K, Bajbouj K, Ullrich O, et al. Cutting edge: Chk1 directs senescence and mitotic catastrophe in recovery from G2 checkpoint arrest. J Cell Mol Med. 2011;15(7):1528–41.
Sandhu C, Donovan J, Bhattacharya N, Stampfer M, Worland P, Slingerland J. Reduction of cdc25a contributes to cyclin E1-cdk2 inhibition at senescence in human mammary epithelial cells. Oncogene. 2000;19(47):5314–23.
Kaistha BP, Honstein T, Müller V, Bielak S, Sauer M, Kreider R, et al. Key role of dual specificity kinase TTK in proliferation and survival of pancreatic cancer cells. Br J Cancer. 2014;111(9):1780–7.
Tomida J, Takata K, Lange SS, Schibler AC, Yousefzadeh MJ, Bhetawal S, et al. REV7 is essential for DNA damage tolerance via two REV3L binding sites in mammalian DNA polymerase. Nucleic Acids Res. 2015;43(2):1000–11.
Baker DJ, Weaver RL, vanDeursen JM. Both attenuates and drives senescence and aging in BubR1 progeroid mice. Cell Rep. 2013;3(4):1164–74.
Mjelle R, Hegre SA, Aas PA, Slupphaug G, Drablos F, Saetrom P, et al. Cell cycle regulation of human DNA repair and chromatin remodeling genes. DNA Repair. 2015;30:53–67.
Gopinathan L, Tan SL, Padmakumar VC, Coppola V, Tessarollo L, Kaldis P. Loss of cdk2 and cyclin A2 impairs cell proliferation and tumorigenesis. Cancer Res. 2014;74(14):3870–9.
Luo P, Tresini M, Cristofalo V, Chen X, Saulewicz A, Gray MD, et al. Immortalization in a normal foreskin fibroblast culture following transduction of cyclin A2 or cdk1 genes in retroviral vectors. Exp Cell Res. 2004;294(2):406–19.
Quadri RA, Arbogast A, Phelouzat MA, Boutet S, Plastre O, Proust JJ. Age-associated decline in cdk1 activity delays cell cycle progression of human T lymphocytes. J Immunol. 1998;161(10):5203–9.
Zalzali H, Nasr B, Harajly M, Basma H, Ghamloush F, Ghayad S, et al. CDK2 transcriptional repression is an essential effector in p53 dependent cellular senescence-implications for therapeutic intervention. Mol Cancer Res. 2015;13(1):29–40.
Bartkova J, Rezaei N, Liontos M, Karakaidos P, Kletsas D, Issaeva N, et al. Oncogene-induced senescence is part of the tumorigenesis barrier imposed by DNA damage checkpoints. Nature. 2006;444(7119):633–7.
Gonzalez S, Klatt P, Delgado S, Conde E, Lopez-Rios F, Sanchez-Cespedes M, et al. Oncogenic activity of cdc6 through repression of the INK4/ARF locus. Nature. 2006;440(7084):702–6.
Gavin EJ, Song B, Wang Y, Xi Y, Ju J. Reduction of Orc6 expression sensitizes human colon cancer cells to 5-fluorouracil and cisplatin. PLoS ONE. 2008;3(12):e4054.
Choi YJ, Sicinski P. Unexpected outcomes of CDK4/6 inhibition. Oncotarget. 2013;4(2):176–7.
Zhang L, Yang Z, Ma A, Qu Y, Xia S, Xu D, et al. Growth arrest and DNA damage 45G down-regulation contributes to Janus kinase/signal transducer and activator of transcription 3 activation and cellular senescence evasion in hepatocellular carcinoma. Hepatology. 2014;59(1):178–89.
Iglesias-Ara A, Zenarruzabeitia O, Fernandez-Rueda J, Sanchez-Tillo E, Field SJ, Celada A, et al. Accelerated DNA replication in E2F1 and E2F2 deficient macrophages leads to induction of the DNA damage response and p21 (CIP1) dependent senescence. Oncogene. 2010;29(14):5579–90.
Musio A, Montagna C, Zambroni D, Indino E, Barbieri O, Citti L, et al. Inhibition of BUB1 results in genomic instability and anchorage-independent growth of normal human fibroblasts. Cancer Res. 2003;63(11):2855–63.
Dekker P, Gunn D, McBryan T, Dirks RW, van Heemst D, Lim FL, et al. Microarray based identification of age-dependent differences in gene expression of human dermal fibroblasts. Mech Ageing Dev. 2012;133:498–507.
Voutetakis K, Chatziioannou A, Gonos ES, Trougakos IP. Comparative meta-analysis of transcriptome data during cellular senescence and in vivo tissue ageing. Oxid Med Cell Longev. 2015;2015:732914.
Ferrando-Martinez S, Ruiz-Mateos E, Dudakov JA, Velardi E, Grillari J, Kreil DP, et al. Wnt signaling suppression in the senescent human thymus. J Gerontol A Biol Sci Med Sci. 2015;70(3):273–81.
Wagner W, Bork S, Horn P, Krunic D, Walenda T, Diehlmann A, et al. Aging and replicative senescence have related effects on human stem and progenitor cells. PLoS ONE. 2009;4(6):e5846.
Massoud A, Rezaei N. Immunology of Aging. Heidelberg: Springer; 2014. p. 221.
Le Morvan C, Cogne M, Drouet M. HLA-A and HLA-B transcription decrease with ageing in peripheral blood leucocytes. Clin Exp Immunol. 2001;125(2):245–50.
Gong Z, Kennedy O, Sun H, Wu Y, Williams GA, Klein L. Reductions in serum IGF-1 during aging impair health span. Aging Cell. 2014;13(3):408–18.
Kenyon C. A conserved regulatory system for Aging. Cell. 2001;105(2):165–8.
Sidler C, Woycicki R, Kovalchuk I, Kovalchuk O. WI-38 senescence is associated with global and site-specific hypomethylation. Aging. 2014;6(7):564–74.
Natarajan E, Omobono JD, Guo Z, Hopkinson S, Lazar AJ, Brenn T, et al. A keratinocyte hypermotility/growth arrest response involving laminin 5 and p16INK4A activated in wound healing and senescence. Am J Pathol. 2006;168(6):1821–37.
Pazolli E, Luo X, Brehm S, Carbery K, Chung JJ, Prior JL, et al. Senescent stromal-derived osteopontin promotes preneoplastic cell growth. Cancer Res. 2009;69(3):1230–9.
Saadane A, Mast N, Charvet CD, Omarova S, Zheng W, Huang SS, et al. Retinal and nonocular abnormalities in Cyp27a1-/-Cyp46a1-/- mice with dysfunctional metabolism of cholesterol. Am J Pathol. 2014;184(9):2403–19.
Ha MK, Chung KY, Lee JH, Bang D, Park YK, Lee KH. Expression of psoriasis-associated fatty acid binding protein in senescent human dermal microvascular endothelial cells. Exp Dermatol. 2004;13(9):543–50.
Gan Q, Huang J, Zhou R, Niu J, Zhu X, Wang J, et al. PPARγ accelerates cellular senescence by inducing p16INK4α expression in human diploid fibroblasts. J Cell Sci. 2008;121:2235–45.
Mooijaart SP, Kuningas M, Westendorp RG, Houwing-Duistermaat JJ, Slagboom PE, Rensen PC, et al. Liver X receptor alpha associates with human life span. J Gerontol A Biol Sci Med Sci. 2007;62(4):342–9.
Walter S, Atzmon G, Demerath EW, Garcia ME, Kaplan RC, Kumari M, et al. A genome wide association study of aging. Neurobiol Aging. 2011;32(11):15–28.
Choi KH, Zepp ME, Higgs BW, Weickert CS, Webster MJ. Expression profiles of schizophrenia susceptibility genes during human prefrontal cortical development. J Psychiatry Neurosci. 2009;34(6):450–8.
Yi JM, Dhir M, Van Neste L, Downing SR, Jeschke J, Glöckner SC, et al. Genomic and epigenomic integration identifies a prognostic signature in colon cancer. Clin Cancer Res. 2011;17(6):1535–45.
Cen B, Deguchi A, Weinstein IB. Activation of protein kinase G increases the expression of p21CIP1, p27KIP1 and histidine triad protein 1 through Sp1. Cancer Res. 2008;68(13):5355–62.
Lodygin D, Menssen A, Hermeking H. Induction of the Cdk inhibitor p21 by LY83583 inhibits tumor cell proliferation in a p53-independent manner. J Clin Invest. 2002;110(11):1717–27.
Kaneda A, Fujita T, Anai M, Yamamoto S, Nagae G, Morikawa M, et al. Activation of Bmp2-Smad1 signal and its regulation by coordinated alteration of H3K27 trimethylation in Ras-induced senescence. PLoS Genet. 2011;7(11):e1002359.
Coppola D, Balducci L, Chen DT, Loboda A, Nebozhyn M, Staller A, et al. Senescence associated gene signature identifies genes linked to age, prognosis and progression of human gliomas. J Geriatr Oncol. 2014;5(4):389–99.
Chien Y, Scuoppo C, Wang X, Fang X, Balgley B, Bolden JE, et al. Control of the senescence associated secretory phenotype by NF-kB promotes senescence and enhances chemosensitivity. Genes Dev. 2011;25(20):2125–36.
Rouault C, Pellegrinelli V, Schilch R, Cotillard A, Poitou C, Tordjman J, et al. Roles of chemokine ligand-2 (CXCL2) and neutrophils in influencing endothelial cell function and inflammation of human adipose tissue. Endocrinology. 2013;154(3):1069–79.
Ma Y, Adjemian S, Galluzzi L, Zitvogel L, Kroemer G. Chemokines and chemokine receptors required for optimal responses to anticancer chemotherapy. Oncoimmunology. 2014;3(1):e27663.
Kim HJ, Kim KW, Yu BP, Chung HY. The effect of age on cyclooxygenase-2 gene expression: NF-kappaB activation and IkappaBalpha degradation. Free Radic Biol Med. 2000;28(5):683–92.
Bai Y, Ding Y, Spencer S, Lasky LA, Bromberg JS. Regulation of the association between PSTPIP and CD2 in murine T cells. Exp Mol Pathol. 2001;71(2):115–24.
Manns J, Daubrawa M, Driessen S, Paasch F, Hoffmann N, Löffler A, et al. Triggering of a novel intrinsic apoptosis pathway by the kinase inhibitor staurosporine: activation of caspase-9 in the absence of Apaf-1. FASEB J. 2011;25(9):3250–61.
Thuret G, Chiquet C, Herrag S, Dumollard JM, Boudard D, Bednarz J, et al. Mechanisms of staurosporine induced apoptosis in a human corneal endothelial cell line. Br J Ophthalmol. 2003;87(3):346–52.
Johansson AC, Steen H, Ollinger K, Roberg K. Cathepsin D mediates cytochrome c release and caspase activation in human fibroblast apoptosis induced by staurosporine. Cell Death Differ. 2003;10(11):1253–9.
Kim S, Ryu S, Kang H, Choi H, Park S. Defective nuclear translocation of stress-activated signaling in senescent diploid human fibroblasts: a possible explanation for aging-associated apoptosis resistance. Apoptosis. 2011;16(8):795–807.
SM undertook the laboratory work in the laboratory of PH, CK and SD. SM wrote the manuscript together with SD. SP and UW did the bioinformatics analysis under the guidance of RG. MG undertook the high throughput RNA sequencing under the guidance of MP. All authors read the manuscript, studied it critically for its intellectual content and approved the final draft. All authors read and approved the final manuscript.
The work described here is part of the research programme of the Jena Centre for Systems Biology of Ageing-JenAge. We acknowledge JenAge funding by the German Ministry for Education and Research (Bundesministerium für Bildung und Forschung—BMBF; support code: 0315581).
The authors declare that they have no competing interests.
All reads have been deposited in the NCBI GEO under the accession number GSE77682 and will be made available at the time of publication.
No ethics approval was needed for this particular study.
Additional file 1: Figure S1. Induction or lack of apoptotic induction in primary human fibroblast strains subjected to 20 Gy irradiation. The blots show the lack of induction of two different apoptotic markers (Caspase-3 or BAX) in primary human fibroblast strains (HFF and MRC-5), 120 min after 20 Gy irradiation. Induction of either of the apoptotic markers in non-irradiated early PD fibroblasts were used as controls. As a positive control, the fibroblasts were subjected to induction of apoptosis by treatment with 2.0 μM Stauroporine for 96 h [35, 205–208].
Additional file 2: Figure S2. Heatmap showing the intersection of the most differentially expressed genes in each of the fibroblast strains (irradiated versus controls). Heatmap illustrating the log2 fold change of gene expression when comparing irradiated HFF versus controls (upper part) and irradiated MRC-5 versus controls (lower part) respectively. The horizontal axis displays the genes selected for this comparison. Genes were selected by intersecting the 200 most differentially regulated genes for each condition. This intersection contains 46 genes. The color key (top left) relates heatmap color to log2 fold change. Red color indicates a negative log2 fold change, i.e. a down-regulation under the second condition compared to the first condition, while the yellow color indicates a positive log2 fold change, i.e. an up-regulation under the second condition relative to the first condition. The dendro- gram on top of the plot clusters the genes into groups with similar expression levels for both comparisons. While most of the genes show similar log2 fold changes for both comparisons, some genes are up-regulated for one comparison, and down-regulated for the other comparison (CTGF, COLEC12, KRT19 in the middle of the plot).
Additional file 3: Figure S3. Regulation of genes of Arrhythmogenic right ventricular cardiomyopathy pathway during senescence induction in HFF strains Genes of the “Arrhythmogenic right ventricular cardiomyopathy” pathway which are significantly up- (green) and down- (red) regulated (log2 fold change >1) during irradiation induced senescence (120 h after 20 Gy irradiation) in HFF strains. Orange color signifies genes which are commonly up-regulated during both, irradiation induced and replicative senescence.
Additional file 4: Figure S4. Regulation of genes of Cell adhesion pathway during senescence induction in HFF strains. Genes of the “Cell adhesion” pathway which are significantly up- (green) and down- (red) regulated (log2 fold change >1) during irradiation induced senescence (120 h after 20 Gy irradiation) in HFF strains. Orange and blue colors signify genes which are commonly up- (orange) and down-regulated (blue) during both, irradiation induced and replicative senescence.
Additional file 5: Figure S5. Regulation of genes of Dilated cardiomyopathy pathway during senescence induction in HFF strains. Genes of the “Dilated cardiomyopathy” pathway which are significantly up- (green) and down- (red) regulated (log2 fold change >1) during irradiation induced senescence (120 h after 20 Gy irradiation) in HFF strains. Orange color signifies genes which are commonly up-regulated during both, irradiation induced and replicative senescence.
Additional file 6: Figure S6. Regulation of genes of ECM receptor interaction pathway during senescence induction in HFF strains. Genes of the “ECM receptor interaction” pathway which are significantly up- (green) and down- (red) regulated (log2 fold change >1) during irradiation induced senescence (120 h after 20 Gy irradiation) in HFF strains. Orange and blue colors signify genes which are commonly up- (orange) and down-regulated (blue) during both, irradiation induced and replicative senescence.
Additional file 7: Figure S7. Regulation of genes of PPAR signaling pathway during senescence induction in HFF strains. Genes of the “PPAR signaling” pathway which are significantly up- (green) and down- (red) regulated (log2 fold change >1) during irradiation induced senescence (120 h after 20 Gy irradiation) in HFF strains. Orange and blue colors signify genes which are commonly up- (orange) and down-regulated (blue) during both, irradiation induced and replicative senescence.
Additional file 8: Figure S8. Regulation of genes of Long term depression pathway during senescence induction in HFF strains. Genes of the “Long term depression” pathway which are significantly up- (green) and down- (red) regulated (log2 fold change >1) during irradiation induced senescence (120 h after 20 Gy irradiation) in HFF strains. Orange color signifies genes which are commonly up- regulated during both, irradiation induced and replicative senescence.
Additional file 9. Regulation and function of genes associated with repair pathways during irradiation induced senescence. Differential regulation and function of (i) genes associated with DNA repair (Nucleotide excision, Base excision, Mismatch repair) and (ii) DNA replication pathways on 20 Gy irradiation induced senescence in MRC-5 fibroblasts and HFF strains. Blue color denotes that previous studies have functionally associated these genes with induction of senescence.
Additional file 10. Function of genes associated with common differentially regulated pathways during replicative and pre-mature senescence induction. Function of each of the genes associated with the commonly differentially regulated pathways (up- and down-regulated) on replicative and 20 Gy irradiation induced senescence in the MRC-5 fibroblasts and HFF strains. Blue color denotes that previous studies have functionally associated these genes with induction of senescence.
Additional file 11: Figure S9. Regulation of genes of NOD-like receptor signaling pathway during senescence induction in MRC-5 strains. Genes of the “NOD-like receptor signaling” pathway which are significantly down- (red) regulated (log2 fold change >1) during irradiation induced senescence (120 h after 20 Gy irradiation) in MRC-5 fibroblast strains. Blue color signifies genes which are commonly down-regulated during both, irradiation induced and replicative senescence.
Additional file 12: Figure S10. Regulation of genes of TGF-beta signaling pathway during senescence induction in MRC-5 strains. Genes of the “TGF-beta signaling” pathway which are significantly up- (green) and down- (red) regulated (log2 fold change >1) during irradiation induced senescence (120 h after 20 Gy irradiation) in MRC-5 fibroblast strains. Blue color signifies genes which are commonly down-regulated during both, irradiation induced and replicative senescence.
About this article
Cite this article
Marthandan, S., Menzel, U., Priebe, S. et al. Conserved genes and pathways in primary human fibroblast strains undergoing replicative and radiation induced senescence. Biol Res 49, 34 (2016). https://doi.org/10.1186/s40659-016-0095-2