Research Paper Volume 7, Issue 11 pp 986—999

Screening of a kinase library reveals novel pro-senescence kinases and their common NF-κB-dependent transcriptional program

Mylène Ferrand1,2,3,4, , Olivier Kirsh5, , Audrey Griveau1,2,3,4, , David Vindrieux1,2,3,4, , Nadine Martin1,2,3,4, , Pierre-Antoine Defossez5, , David Bernard1,2,3,4, ,

  • 1 Inserm U1052, Centre de Recherche en Cancérologie de Lyon, F-69373 Lyon, France
  • 2 CNRS UMR 5286, F-69373 Lyon, France
  • 3 Centre Léon Bérard, F-69373 Lyon, France
  • 4 Université de Lyon, F-69373 Lyon, France
  • 5 Epigenetics and Cell Fate, University Paris Diderot, Sorbonne Paris Cité, UMR 7216 CNRS, 75013 Paris, France

Received: September 4, 2015       Accepted: November 2, 2015       Published: November 15, 2015
How to Cite

Copyright: © 2015 Ferrand et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.


Cellular senescence results in proliferation arrest and acquisition of hallmarks such as the Senescence-Associated Secretory Phenotype (SASP). Senescence is involved in regulating numerous physio-pathological responses, including embryonic development, cancer, and several aging-related diseases. Only a few kinases, centered on the RAS signaling pathway, have been identified as inducing premature senescence. About possible other senescence-regulating kinases and signaling pathways, practically little is known. By screening a library of activated kinases, we identified 33 kinases whose constitutive expression decreases cell proliferation and induces expression of senescence markers; p16 and SASP components. Focusing on some kinases showing the strongest pro-senescence effects, we observed that they all induce expression of SASP-component genes through activation of an NF-κB-dependent transcriptional program. Furthermore, inhibition of the p53 or Rb pathway failed to prevent the SASP-inducing effect of pro-senescence kinases. Inhibition of the NF-κB, p53, or Rb pathway proved insufficient to prevent kinase-triggered cell cycle arrest. We have thus identified a repertoire of novel pro-senescence kinases and pathways. These results will open new perspectives in the understanding on the role of cellular senescence in various physio-pathological responses.


Cellular senescence is a stable form of cell cycle arrest accompanied by acquisition of a Senescence-Associated Secretory Phenotype (SASP), which notably includes production of pro-inflammatory cytokines. Cell senescence participates in various biological processes, including embryonic development, aging-related diseases, wound healing, and tumor protection [1-5].

Activation of cellular senescence in response to aberrant oncogenic activation is a well-described tumor-protective mechanism. In otherwise normal cells, oncogenic activation can trigger cell senescence, blocking the proliferation of potentially defective cells and causing them to develop a SASP that can promote their elimination by the immune system [6]. This mechanism, called Oncogene-Induced Senescence (OIS), is well characterized, and several players involved in it have been identified. Those having attracted the most attention are related to the RAS signaling pathway: oncogenic receptor tyrosine kinases upstream from RAS signaling, MEK kinases in the RAF-MEK-ERK branch downstream of RAS signaling, and the PI3K-AKT kinase branch downstream of RAS signaling have all been found to induce OIS. Despite the importance of these players, the wide range of physio-pathological responses in which cell senescence plays a role suggests that additional signaling pathways and kinases contribute to regulating cell senescence.

To identify cell-senescence-regulating kinases, we have exploited an activated kinase library containing about 200 constitutively active kinases [7]. We show that 33 kinases of this library can induce premature senescence in normal fibroblasts. This suggests that cell senescence can result from activation of diverse signaling pathways, likely to be more numerous than initially believed. We further show that the strongest pro-senescence kinases identified in this screen all induce the SASP via the NF-κB pathway, and that the observed kinase-triggered senescence cannot be alleviated by inhibiting either the Rb or the p53 pro-senescence pathway.


A kinase screen identifies new pro-senescence kinases and pathways

Little is known about the ability of activated kinases other than those acting upstream or downstream from RAS signaling to induce premature senescence. To identify novel cell-senescence-inducing kinases, we stably transduced each kinase gene of a previously described and validated library [7] into normal human fibroblasts (Supplementary Table 1). A first selection on the basis of decreased cell proliferation yielded 53 kinases (Supplementary Table 2). We next examined the effects of these 53 kinases on the levels of transcripts of four well-described SASP components: IL1A (Fig. 1A), IL1B (Fig. 1B), IL6 (Fig. 1C), and IL8 (Fig. 1D) [8,9]. Most of the tested kinases proved able to induce expression of the SASP-component genes, sometimes very strongly (in some cases, more than 100-fold induction was observed) (Fig. 1A-D). Generally, variations in the level of one transcript correlated with changes in the levels of the other three. In others words, the more strongly a kinase induced expression of a SASP-component gene, the more strongly it induced expression of the others (Fig. 1E-F). We found 44 kinases to both decrease cell proliferation and induce expression of SASP-component genes (Supplementary Table 3).

Profiles of SASP component induction by anti-proliferative kinases

Figure 1. Profiles of SASP component induction by anti-proliferative kinases. (A-D) IMR-90 normal human fibroblasts were infected with a control vector or the indicated kinase-encoding retroviral vector. Five days after infection, RNAs were prepared and the indicated SASP-component transcript was quantified by RT-qPCR. Results were normalized with respect to the level of ACTB mRNA. (E) Relative (log2) fold change induction for each cytokine is represented on a heatmap versus the inducing kinase. (F) Correlograms and correlation matrix (non-parametric Spearman method) showing pairwise correlations between fold change induction datasets. Rho values and associated p-values are displayed.

We next focused on p16 transcripts. The p16 protein is a known cyclin-dependent kinase inhibitor (CDKI) whose expression generally increases in senescent cells [10]. The level of p16 transcripts was found to increase in response to numerous kinases, but the fold induction did not exceed about 2.5, in contrast to what was observed with SASP-component transcripts (Fig. 2A). We found 37 kinases both to decrease cell proliferation and to induce p16 gene expression (Supplementary Table 4). As shown in a Venn diagram compiling all the data, we found 33 kinases to decrease cell proliferation and to increase expression of both the p16 gene and the SASP-component genes (Fig. 2B and Table 1). Although the same 33 kinases induced expression of p16 and the SASP, p16 and SASP components appeared in different main branches after hierarchical clustering (Fig. 2C). This suggests that they might be regulated by different transcriptional programs. Still, there was a significant correlation between p16 and the different SASP components (Fig. 2D). A KEGG pathway analysis of these 33 kinases reveals over-representation of some signaling pathways (Table 2) and some of them represent new senescence regulating pathways. In conclusion, we have identified 33 kinases and some associated signaling pathways capable of inducing senescence.

Effects of anti-proliferative kinases on p16 and correlation with SASP-components

Figure 2. Effects of anti-proliferative kinases on p16 and correlation with SASP-components. (A) IMR-90 normal human fibroblasts were infected with a control vector or the indicated kinase-encoding retroviral vector. Five days after infection, RNAs were prepared and the p16 transcript was quantified by RT-qPCR. Results were normalized with respect to the level of ACTB mRNA. (B) A Venn diagram was drawn to summarize the data. (C) For p16 and for each cytokine, relative (log2) fold change induction is represented on a heatmap versus the inducing kinase. (D) Correlograms and correlation matrix (non-parametric Spearman method) showing pairwise correlations between fold change induction datasets. Rho values and associated p-values are displayed.

Table 1. List of kinases showing decreased cell proliferation and increased p16 and SASP-component induction


Table 2. Overrepresented pathways for the 33 pro-senescence kinases

CategoryTerm%PValueGenesList TotalPop HitsPop TotalFold Enrichment
KEGG_PATHWAYrno04660:T cell receptor signaling pathway1.361.32E-5MAP3K7, MAPK12, PAK4, PIK3R5, CDK4, MAP2K718109559017.09
KEGG_PATHWAYrno04930:Type II diabetes mellitus0.913.95E-4HK3, PKM2, PIK3R5, PRKCD1849559025.35
KEGG_PATHWAYrno04664:Fc epsilon RI signaling pathway0.910.0014MAPK12, PIK3R5, MAP2K7, PRKCD1876559016.34
KEGG_PATHWAYrno04620:Toll-like receptor signaling pathway0.910.0023MAP3K7, MAPK12, PIK3R5, MAP2K71890559013.80
KEGG_PATHWAYrno04722:Neurotrophin signaling pathway0.910.0060MAPK12, PIK3R5, MAP2K7, PRKCD1812655909.86
KEGG_PATHWAYrno05223:Non-small cell lung cancer0.680.0105PDPK1, PIK3R5, CDK41852559017.92
KEGG_PATHWAYrno04012:ErbB signaling pathway0.680.0268PAK4, PIK3R5, MAP2K71885559010.96
KEGG_PATHWAYrno04666:Fc gamma R-mediated phagocytosis0.680.0286LIMK1, PIK3R5, PRKCD1888559010.59
KEGG_PATHWAYrno04912:GnRH signaling pathway0.680.0323MAPK12, MAP2K7, PRKCD189455909.91
KEGG_PATHWAYrno00900:Terpenoid backbone biosynthesis0.4530.0418MVK, PMVK1814559044.36
KEGG_PATHWAYrno04010:MAPK signaling pathway0.910.0442MAP3K7, MAP3K6, MAPK12, MAP2K71826655904.67
KEGG_PATHWAYrno04910:Insulin signaling pathway0.680.0598PDPK1, HK3, PIK3R51813255907.06
KEGG_PATHWAYrno04062:Chemokine signaling pathway0.680.0937GRK6, PIK3R5, PRKCD1817155905.45

NF-κB pathway activation is a shared characteristic of several pro-senescence kinases

To see what the identified kinases might have in common, we first interrogated the STRING database to see if at least those displaying the greatest pro-senescence effect, as judged from the level of SASP induction, might share a common protein association network. We built a network based on the seven kinases showing the strongest SASP-inducing effect (Fig. 3). In a KEGG pathway analysis of this network performed to identify overrepresented pathways, we detected strong enrichment in actors of the NF-κB signaling pathway (p=5.8-16) (green circles, Fig. 3). Interestingly, SASP induction during replicative senescence or in response to senescence triggering by the RAS signaling pathway, through RAS or MEK activation, is reported to be mediated by the NF-κB transcription factors [11,12]. These factors are thus good candidate mediators of SASP-component induction by the pro-senescence kinases identified in our screen.

Strongly SASP-inducing kinases display strong associations with the NF-κB signaling pathway

Figure 3. Strongly SASP-inducing kinases display strong associations with the NF-κB signaling pathway. The String database ( used to build networks on the basis of known and predicted protein-protein interactions was interrogated with, as entries, the indicated pro-senescence kinases (orange circles). Enrichment in KEGG pathways was used to identify common signaling pathways for the pro-senescence kinases. Proteins defined as belonging to the NF-κB signaling pathway are highlighted by green circles.

To confirm experimentally the results of this bioinformatic analysis, we first investigated the ability of the three strongest SASP and p16 inducers identified here (Fig. 2C): the MAP3K7, PRKCD, and MATK kinases to activate the NF-κB pathway using a 3κB-Luc reporter vector [13]. The luciferase activity was strongly induced by the 3 kinases indicating that these pro-senescent kinases activate the NF-κB pathway. We next investigated the role of the NF-κB pathway in mediating SASP-component induction by these kinases. We transduced each pro-senescent kinase into human fibroblasts, in combination or not with the gene encoding the IκBα super-repressor (IKBAm), a stabilized NF-κB inhibitor. As expected, the levels of IL1A, IL1B, IL6, and IL8 transcripts strongly increased upon constitutive expression of the kinases (Fig. 4B-E). Strikingly, this induction was almost completely abolished when NF-κB activity was inhibited by IKBAm (Fig. 4B-E). These data demonstrate that SASP component induction by the tested pro-senescence kinases depends on NF-κB transcription factor activity, as anticipated from the STRING bioinformatic analysis.

NF-κB activity mediates the effect of pro-senescence kinases on SASP expression

Figure 4. NF-κB activity mediates the effect of pro-senescence kinases on SASP expression. (A) 293GP cells were co-transfected by the 3κB-Luc reporter and the indicated kinases or the empty vector. Two days later, cells were lysed and Luc activity measured and normalized to protein concentration for each condition. Results are presented as fold Luc activity induced by the indicated kinase when compared to the control vector. (B-E) MRC-5 normal human fibroblasts were co-infected with the indicated retroviral vectors. Five days later, RNAs were prepared and RT-qPCR was performed to quantify (B) IL1A, (C) IL1B, (D) IL6, and (E) IL8 transcripts. Results were normalized with respect to the level of ACTB mRNA. Statistical analysis was performed with the T-test, p<0.005 ***.

We next wondered if these pro-senescence kinases might activate a broader NF-κB-dependent program affecting more than just SASP components. To answer this question, we examined the effect of our 33 kinases on transcript-level expression of the genes encoding three known direct intracellular NF-κB targets: IκBα, SOD2, and COX2. The last two genes are known to be involved in NF-κB-induced senescence [14-16]. The kinases were found to induce expression of all three genes (Fig. 5A-C). Upon hierarchical clustering of the SASP components, the intracellular NF-κB targets, and p16 on the basis of their induction profiles, the SASP components and intracellular NF-κB targets were found to cluster on the same main branch, whereas p16 was still on an independent branch (Fig. 5D). This again suggests that expression of the p16 gene is not controlled directly by an NF-κB-dependent transcriptio-nal program. Nevertheless, correlations were observed between p16, individual SASP-component, and individual intracellular NF-κB-target mRNA levels (Fig. 5E).

Intracellular NF-κB targets induction by anti-proliferative kinases correlates to SASP components

Figure 5. Intracellular NF-κB targets induction by anti-proliferative kinases correlates to SASP components. (A-C) IMR-90 normal human fibroblasts were infected with a control vector or the indicated kinase-encoding retroviral vector. Five days after infection, RNAs were prepared and RT-qPCR was performed to quantify (A) IKBA, (B) SOD2, and (C) COX2 transcripts. Results were normalized with respect to the level of ACTB mRNA. (D) For p16, each cytokine, and each NF-κB target, relative (log2) fold induction is indicated on a heatmap versus the inducing kinase. (E) Correlograms and correlation matrix (non-parametric Spearman method) showing correlations between fold change induction datasets. Rho values and associated p-values are displayed.

SASP components can propagate and amplify the senescence signal [8,9,11,17]. COX2 and SOD2 can also promote cellular senescence induced by Rel/NF-κB transcription factors [14-16]. As the responses of all these factors to our pro-senescence kinases appear to be under the transcriptional control of NF-κB, we examined whether constitutive NF-κB inhibition by IKBAm might prevent senescence induction by pro-senescence kinases. As expected, constitutive expression of PRKCD, MATK, or MAP3K7 blocked cell proliferation (Fig. 6A, upper panel) and induced p16 expression (Fig. 6B). Inhibiting NF-κB with IKBAm did not prevent the kinase-promoted proliferation arrest (Fig. 6A) and only partially inhibited p16 induction (Fig. 6B). This suggests that the tested pro-senescence kinases, although strongly inducing transcription of SASP-component and other pro-senescence regulator genes via NF-κB transcription factors, do not induce senescence solely through the NF-κB pathway.

NF-κB target induction is not the sole mediator of senescence by the pro-senescence kinases

Figure 6. NF-κB target induction is not the sole mediator of senescence by the pro-senescence kinases. (A) Fifty thousand MRC-5 normal human fibroblasts were seeded per well in 6-well plates. The next day, cells were infected with a retroviral vector encoding the indicated kinase. Ten days after infection, the cells were fixed and crystal violet stained. (B) Five days after infection, RNAs were prepared and p16 transcripts were quantified by RT-qPCR. Results were normalized with respect to the level of ACTB mRNA. Statistical analysis was performed with the T-test, p<0.01 **, p<0.05 *.

Inhibition of the p53 or Rb pathway does not prevent senescence or induction of SASP-component gene expression

As the p53 and p16/Rb pathways are important mediators of cell senescence [18,19], we wondered whether they might mediate the cell proliferation arrest and upregulation of SASP-component expression observed in response to pro-senescence kinases. We thus used protein E6 to inhibit the P53 pathway and protein E7 to inhibit the p16/Rb pathway. The ability of E6 or E7 to strongly inhibit its target pathway was demonstrated by acceleration of cell proliferation in normal human fibroblasts infected with an E6- or E7-encoding vector, as compared to cells infected with a control vector (Fig. 7A). Yet even in the presence of E6 or E7, constitutive expression of the MAP3K7, PRKCD, or MATK kinases was still found to promote cell proliferation arrest (Fig. 7B) and increased IL1A, IL1B, IL6, and IL8 transcript levels (Fig. 7C-F). Combination of p53 and NF-kB inhibition has been showed to prevent senescence induced by an oncogenic RAS [12]. We then tested whether MAP3K7-, PRKCD-, or MATK-induced senescence might be overcome by inhibition of p53 by E6 and NF-κB by IKBAm. Inhibition of p53 and NF-κB largely reverted senescence induced by MAP3K7 or PRKCD but had no effect on senescence induced by MATK (Supplemental Figure).

p53 or p16/Rb pathway inhibition did not revert kinase-induced senescence

Figure 7. p53 or p16/Rb pathway inhibition did not revert kinase-induced senescence. (A) One hundred thousand MRC-5 normal human fibroblasts were seeded per well in 6-well plates. The next day, the cells were infected with an E6- or E7-encoding retroviral vector and five days later they were fixed and crystal violet stained. (B) Fifty thousand MRC-5 normal human fibroblasts were seeded per well in 6-well plates. The next day, they were infected with the indicated retroviral vector. Thirteen days after infection, the cells were fixed and crystal violet stained. (C-F) MRC-5 cells were infected with the indicated retroviral vector and 5 days later, RNAs were prepared. RT-qPCR was performed to quantify (C) IL1A, (D) IL1B, (E) IL6 and (F) IL8 transcripts. Results were normalized with respect to the level of ACTB mRNA. Statistical analysis was performed with the T-test, p<0.005 ***.

Together, these results demonstrate that inhibition of the p16/Rb or p53 pathway is not sufficient to prevent the proliferation arrest and induction of SASP-component expression triggered by these pro-senescence kinases. Combined inhibition of NF-kB and p53 can revert senescence induced by only some of these pro-senescence kinases.


Despite the growing number of physiological responses in which cellular senescence is known to participate, little is known about the kinases and signaling pathways that induce senescence. Here we have identified 33 pro-senescence kinases, most of which were not previously known to induce senescence. This rather broad range of pro-senescence kinases might reflect the fact that senescence is a cell stress response induced in many situations where cell homeostasis is perturbed.

Our KEGG pathway analysis of these 33 kinases reveals over-representation of expected pathways, such as the “MAPK signaling pathway” itself or other pathways located downstream from receptors that activate it, such as the “T cell receptor signaling pathway”, the “Fc epsilon RI signaling pathway”, and the “Toll-like receptor signaling pathway” (Table 2). Over-representation of a “chemokine signaling pathway” also makes sense, since SASP components, particularly chemokines, are known inducers of cell senescence (Table 2) [9,11,20].

More interestingly, this analysis also reveals new putative pro-senescence pathways. Particularly appealing is the over-representation of the KEGG pathways “Type II diabetes mellitus (T2D)” and “Insulin signaling pathway”, which suggests a possible role of cell senescence in regulating insulin sensitivity and diabetes (Table 2). In line with these results, investigators have linked p16, hallmark of senescent cells, to T2D. For example, T2D susceptibility loci appear associated with the gene encoding p16, and p16 expression increases in tissues derived from diabetic patients. The role of p16 in the etiology of the disease remains unclear, however [21-25,25,26]. The step from these results to demonstrating a role of cell senescence in regulating T2D has not yet been taken, although investigators have begun to discuss this putative role and the possible advantage of targeting senescent cells in order to improve T2D [27,28]. In conclusion, our data support a potential role of cellular senescence in regulating T2D etiology.

The “terpenoid backbone biosynthesis” pathway is also strongly over-represented among our short-listed kinases (Table 2). As far as we know, this pathway has never been implicated in cellular senescence. It leads to production of numerous macromolecules, notably cholesterol, sterol, and ubiquinones. Cellular senescence is thought to be part of a general decrease in health during aging and of the aging process in general [2,27,29]. A decreased MYC level has recently been shown to increase longevity and promote a healthier lifespan in mice [30]. This improvement in the aging process was found to correlate with decreased gene expression in the cholesterol pathway, which is part of the “terpenoid backbone biosynthesis” pathway [30]. Together with our data, these results suggest a possible role for this pathway in promoting cell senescence and aging.

A functional approach on the strongest pro-senescent kinases demonstrates that they all share a common network with NF-κB transcriptional factors and this has been functionally confirmed as NF-κB inhibition is sufficient to block SASP components induction by these kinases. NF-κB is also regulating other intracellular targets showing that senescent cells display a broad chronic NF-κB activation beyond the sole SASP components induction. Focusing on two other pro-senesence pathways, the p53 and p16/Rb pathways, we further demonstrate that inactivation of either one is insufficient to prevent the senescence induced by the identified kinases in normal human fibroblasts. Nevertheless, combined inhibition of p53 and NF-κB is sufficient to reverse senescence by some but not all the pro-senescence kinases tested, suggesting common and distinct pro-senescence pathways induced by the tested kinases.

In conclusion, we have constituted a repertoire of senescence-promoting kinases, including ones not previously known to regulate senescence, and have identified some of the signaling pathways in which they are involved. This repertoire should contribute to establishing and understanding the role cellular senescence plays in the growing list of physio-pathological conditions in which it participates.

Materials and Methods

Cell culture

IMR-90 (ATCC), MRC-5 (ATCC), and GP293 cells (Clontech) were cultured in DMEM. All media were supplemented with 10% FBS (Sigma) and 1% penicillin/streptomycin (Invitrogen). Upon receipt, cells were thawed and amplified and aliquots frozen. Experiments were performed on the aliquots within a month. The cells were maintained at 37°C under a 5% CO2 atmosphere.

Vectors, transfection, and infection

The library used was the Myristoylated Kinase Library described in [7] (Kit # 1000000012, Addgene). The retroviral vectors employed were pBabe-Puro-IκBα-mut (super repressor) [7], pLXN-E6, pLXN-E7, and pLXN-E6E7 [31] (Addgene). The 3κB-Luc reporter vector and the protocol for the transactivation assays were previously described [13]. The protocols used to transfect virus-producing GP293 cells and infect target cells have been described previously [32]. The infection protocols were designed so that practically all cells were infected, as judged from the fluorescence observed after infection with a GFP-expressing retroviral vector.

Colony formation assays

Colony formation assays were carried out in 6-well plates. Depending on the experiment, 50,000 to 100,000 cells were seeded per well. One day later they were infected. Five to ten days after seeding, the cells were washed with PBS, fixed with 4% paraformaldehyde, and stained with 0.05% crystal violet (Sigma-Aldrich).

RNA extraction, reverse transcription, and PCR

To isolate RNA, we used a phenol-chloroform extraction method involving cell lysis in RNA-ISOL lysis reagent (Dutscher). PhaseLockGel tubes (Prime) were used to separate the phases. Thereafter, the Dynamo cDNA Synthesis Kit (Fisher Scientific) was used for cDNA synthesis from 1 μg total RNA. The RT reaction mixture was diluted 1/20 and the cDNA template used for qPCR analysis. TaqMan quantitative PCR analysis was carried out in the CFX96 Connect Real-Time PCR Detection System (Bio-Rad). The FastStart Essential Probes Master (Roche) was used as PCR mix. The ACTB housekeeping gene was used for normalization. Real-time intron-spanning PCR assays were designed with the ProbeFinder software (Roche Applied Science). The following primers and UPL probes were used: ACTB-Forward ATTGGCAATGAGCGGTTC and ACTB-Reverse GGATGCCACAGGACTCCAT (UPL probe: 11), IL8-F agacagcagagcacacaagc and IL8-R atggttccttccggtggt (UPL probe: 72), IL6-F caggagcccagctatgaac and IL6-R gaaggcagcaggcaacac (UPL probe: 7), IL1A-F ggttgagtttaagccaatcca and IL1A-R tgctgacctaggcttgatga (UPL probe: 6), IL1B-F tacctgtcctgcgtgttgaa and IL1B-R tctttgggtaatttttgggatct (UPL probe: 78), p16-F gtggacctggctgaggag and p16-R ctttcaatcggggatgtctg (UPL probe: 34), IKBA-F gtcaaggagctgcaggagat and IKBA-R atggccaagtgcaggaac (UPL probe: 38), COX2-F gctttatgctgaagccctatga and COX2-R tccaactctgcagacatttcc (UPL probe: 2), SOD2-F aatcaggatccactgcaagg and SOD2-R taagcgtgctcccacacat (UPL probe: 3).

Heatmaps and correlations

Heatmaps and correlation analyses were done with R (gplots, corrgram, and Hmisc package). Briefly, normalized RT-qPCR tables were loaded in R and relative fold induction was calculated for each cytokine. One in FC induction is set to Zero, maximum FC is set to One (100% induction). Details of calculation are described as supplementary Methods.


We thank the laboratory members for helpful discussions.


This work was carried out with the support of the AAP “Epigénome et Cancer” of the Plan Cancer (P030830), to PAD and DB. Work in the lab of PAD is supported by funding from Fondation ARC and by ANR-11-LABX-0071 under program ANR-11-IDEX-0005-01.

Conflicts of Interest

The authors of this manuscript have no conflict of interests to declare.


  • 1. Campisi J. Senescent cells, tumor suppression, and organismal aging: good citizens, bad neighbors. Cell. 2005; 120:513-522. [PubMed]
  • 2. Campisi J and Robert L. Cell senescence: role in aging and age-related diseases. Interdiscip. Top. Gerontol. 2014; 39:45-61. [PubMed]
  • 3. Demaria M, Ohtani N, Youssef SA, Rodier F, Toussaint W, Mitchell JR, Laberge RM, Vijg J, Van SH, Dolle ME, Hoeijmakers JH, de Bruin A, Hara E, et al. An essential role for senescent cells in optimal wound healing through secretion of PDGF-AA. Dev. Cell. 2014; 31:722-733. [PubMed]
  • 4. Munoz-Espin D, Canamero M, Maraver A, Gomez-Lopez G, Contreras J, Murillo-Cuesta S, Rodriguez-Baeza A, Varela-Nieto I, Ruberte J, Collado M, Serrano M. Programmed cell senescence during mammalian embryonic development. Cell. 2013; 155:1104-1118. [PubMed]
  • 5. Munoz-Espin D and Serrano M. Cellular senescence: from physiology to pathology. Nat. Rev. Mol. Cell Biol. 2014; 15:482-496. [PubMed]
  • 6. Kang TW, Yevsa T, Woller N, Hoenicke L, Wuestefeld T, Dauch D, Hohmeyer A, Gereke M, Rudalska R, Potapova A, Iken M, Vucur M, Weiss S, et al. Senescence surveillance of pre-malignant hepatocytes limits liver cancer development. Nature. 2011; 479:547-551. [PubMed]
  • 7. Boehm JS1, Zhao JJ, Yao J, Kim SY, Firestein R, Dunn IF, Sjostrom SK, Garraway LA, Weremowicz S, Richardson AL, Greulich H, Stewart CJ, Mulvey LA, et al. Integrative genomic approaches identify IKBKE as a breast cancer oncogene. Cell. 2007; 129:1065-1079. [PubMed]
  • 8. Coppe JP, Patil CK, Rodier F, Sun Y, Munoz DP, Goldstein J, Nelson PS, Desprez PY, Campisi J. Senescence-associated secretory phenotypes reveal cell-nonautonomous functions of oncogenic RAS and the p53 tumor suppressor. PLoS Biol. 2008; 6:2853-2868. [PubMed]
  • 9. Kuilman T, Michaloglou C, Vredeveld LCW, Douma S, van Doorn R, Desmet CJ, Aarden LA, Mooi WJ, Peeper DS. Oncogene-Induced Senescence Relayed by an Interleukin-Dependent Inflammatory Network. Cell. 2008; 133:1019-1031. [PubMed]
  • 10. Rayess H, Wang MB, Srivatsan ES. Cellular senescence and tumor suppressor gene p16. Int. J. Cancer. 2012; 130:1715-1725. [PubMed]
  • 11. Acosta JC, O'Loghlen A, Banito A, Guijarro MV, Augert A, Raguz S, Fumagalli M, Da Costa M, Brown C, Popov N, Takatsu Y, Melamed J, d'Adda di Fagagna F, et al. Chemokine Signaling via the CXCR2 Receptor Reinforces Senescence. Cell. 2008; 133:1006-1018. [PubMed]
  • 12. Chien Y, Scuoppo C, Wang X, Fang X, Balgley B, Bolden JE, Premsrirut P, Luo W, Chicas A, Lee CS, Kogan SC, Lowe SW. Control of the senescence-associated secretory phenotype by NF-kappaB promotes senescence and enhances chemosensitivity. Genes Dev. 2011; 25:2125-2136. [PubMed]
  • 13. Bernard D, Quatannens B, Vandenbunder B, Abbadie C. Rel/NF-kappaB transcription factors protect against tumor necrosis factor (TNF)-related apoptosis-inducing ligand (TRAIL)-induced apoptosis by up-regulating the TRAIL decoy receptor DcR1. J. Biol. Chem. 2001; 276:27322-27328. [PubMed]
  • 14. Bernard D, Quatannens B, Begue A, Vandenbunder B, Abbadie C. Antiproliferative and antiapoptotic effects of crel may occur within the same cells via the up-regulation of manganese superoxide dismutase. Cancer Res. 2001; 61:2656-2664. [PubMed]
  • 15. Bernard D, Gosselin K, Monte D, Vercamer C, Bouali F, Pourtier A, Vandenbunder B, Abbadie C. Involvement of Rel/nuclear factor-kappaB transcription factors in keratinocyte senescence. Cancer Res. 2004; 64:472-481. [PubMed]
  • 16. Zdanov S, Bernard D, Debacq-Chainiaux F, Martien S, Gosselin K, Vercamer C, Chelli F, Toussaint O, Abbadie C. Normal or stress-induced fibroblast senescence involves COX-2 activity. Exp Cell Res. 2007; 313:3046-3056. [PubMed]
  • 17. Acosta JC, Banito A, Wuestefeld T, Georgilis A, Janich P, Morton JP, Athineos D, Kang TW, Lasitschka F, Andrulis M, Pascual G, Morris KJ, Khan S, et al. A complex secretory program orchestrated by the inflammasome controls paracrine senescence. Nat. Cell Biol. 2013; 15:978-990. [PubMed]
  • 18. Itahana K, Campisi J, Dimri G. Mechanisms of cellular senescence in human and mouse cells. Biogerontology. 2004; 5:1-10. [PubMed]
  • 19. 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:593-602. [PubMed]
  • 20. Kuilman T and Peeper DS. Senescence-messaging secretome: SMS-ing cellular stress. Nat Rev Cancer. 2009; 9:81-94. [PubMed]
  • 21. Bantubungi K, Hannou SA, Caron-Houde S, Vallez E, Baron M, Lucas A, Bouchaert E, Paumelle R, Tailleux A, Staels B. Cdkn2a/p16Ink4a regulates fasting-induced hepatic gluconeogenesis through the PKA-CREB-PGC1alpha pathway. Diabetes. 2014; 63:3199-3209. [PubMed]
  • 22. Gonzalez-Navarro H, Vinue A, Sanz MJ, Delgado M, Pozo MA, Serrano M, Burks DJ, Andres V. Increased dosage of Ink4/Arf protects against glucose intolerance and insulin resistance associated with aging. Aging Cell. 2013; 12:102-111. [PubMed]
  • 23. Kluth O, Matzke D, Schulze G, Schwenk RW, Joost HG, Schurmann A. Differential transcriptome analysis of diabetes-resistant and -sensitive mouse islets reveals significant overlap with human diabetes susceptibility genes. Diabetes. 2014; 63:4230-4238. [PubMed]
  • 24. Saxena R, Voight BF, Lyssenko V, Burtt NP, de Bakker PI, Chen H, Roix JJ, Kathiresan S, Hirschhorn JN, Daly MJ, Hughes TE, Groop L, Altshuler D, et al. Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels. Science. 2007; 316:1331-1336. [PubMed]
  • 25. Scott LJ, Mohlke KL, Bonnycastle LL, Willer CJ, Li Y, Duren WL, Erdos MR, Stringham HM, Chines PS, Jackson AU, Prokunina-Olsson L, Ding CJ, Swift AJ, et al. A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants. Science. 2007; 316:1341-1345. [PubMed]
  • 26. Zhang YY, Guo QY, Wu MY, Zang CS, Ma FZ, Sun T, Wang WN, Miao LN, Xu ZG. p16 Expression Is Increased through 12-Lipoxygenase in High Glucose-Stimulated Glomerular Mesangial Cells and Type 2 Diabetic Glomeruli. Nephron. 2015; 130:141-150. [PubMed]
  • 27. Naylor RM, Baker DJ, van Deursen JM. Senescent cells: a novel therapeutic target for aging and age-related diseases. Clin. Pharmacol. Ther. 2013; 93:105-116. [PubMed]
  • 28. Palmer AK, Tchkonia T, LeBrasseur NK, Chini EN, Xu M, Kirkl JL. Cellular Senescence in Type 2 Diabetes: A Therapeutic Opportunity. Diabetes. 2015; 64:2289-2298. [PubMed]
  • 29. Tchkonia T, Zhu Y, van Deursen JM, Campisi J, Kirkl JL. Cellular senescence and the senescent secretory phenotype: therapeutic opportunities. J. Clin. Invest. 2013; 123:966-972. [PubMed]
  • 30. Hofmann JW, Zhao X, De Cecco M, Peterson AL, Pagliaroli L, Manivannan J, Hubbard GB, Ikeno Y, Zhang Y, Feng B, Li X, Serre T, Qi W, et al. Reduced expression of MYC increases longevity and enhances healthspan. Cell. 2015; 160:477-488. [PubMed]
  • 31. Halbert CL, Demers GW, Galloway DA. The E7 gene of human papillomavirus type 16 is sufficient for immortalization of human epithelial cells. J. Virol. 1991; 65:473-478. [PubMed]
  • 32. Vindrieux D, Augert A, Girard CA, Gitenay D, Lallet-Daher H, Wiel C, Le Calvé B, Gras B, Ferrand M, Verbeke S, de Launoit Y, Leroy X, Puisieux A, et al. PLA2R1 mediates tumor suppression by activating JAK2. Cancer Res. 2013; 73:6334-6345. [PubMed]