The expression signature of in vitro senescence resembles mouse but not human aging
© Wennmalm et al.; licensee BioMed Central Ltd. 2005
Received: 16 May 2005
Accepted: 17 November 2005
Published: 16 December 2005
The biological mechanisms that underlie aging have not yet been fully identified. Senescence, a phenomenon occurring in vitro, limits the number of cell divisions in mammalian cell cultures and has been suggested to contribute to aging.
We investigated whether the changes in gene expression that occur during mammalian aging and induction of cellular senescence are similar. We compared changes of gene expression in seven microarray datasets from aging human, mouse and rat, as well as four microarray datasets from senescent cells of man and mouse. The datasets were publicly available or obtained from other laboratories. Correlation measures were used to establish similarities of the expression profiles and gene ontology analyses to identify functional groups of genes that are co-regulated. Robust similarities were established between aging in different species and tissues, indicating that there is an aging transcriptome. Although some cross-species comparisons displayed high correlation, intra-species similarities were more reliable. Similarly, a senescence transcriptome was demonstrated that is conserved across cell types. A similarity between the expression signatures of cellular senescence and aging could be established in mouse, but not in human.
Our study is the first to use microarray data from several studies and laboratories for dissection of a complex biological phenotype. We demonstrate the presence of a mammalian aging transcriptome, and discuss why similarity between cellular senescence and aging is apparent in aging mice only.
In vitro senescence was discovered as a phenomenon that limits the replicative lifespan of cells grown in culture, and was immediately hypothesised to be a possible cause of human aging [1, 2]. During several decades, attempts have been made to establish a link between the age and the replicative potential of cells derived from humans . Although several studies have shown a correlation between in vivo and in vitro life span, these results have been questioned and re-evaluated . A more convincing approach to establish a link between cellular senescence and aging would be to detect senescent cells in vivo. Unfortunately, no good senescence marker exists, except for senescence-associated β-galactosidase (SAβ-gal) staining . This technique has been used extensively but its specificity has been questioned . A correlation between human age and an accumulation of SA-β GAL staining cells has been reported  but, for unknown reasons, could not be repeated . In mice, increased SAβ-GAL staining has been reported in animals harbouring the Werner mutation and short telomeres, which gives rise to an aging phenotype , following chemotherapeutic treatment , in liver regeneration after partial hepatectomy in animals with short telomeres  and upon aging . These studies indicate that senescent cells could exist in vivo but do not provide any measure of their relative contribution to the aging process.
In the current study, publicly available microarray data describing both aging and in vitro senescence were used. By comparing the expression changes that occur during aging in human, mouse and rat across several tissues, we conclude that there is a mammalian aging expression signature. This expression signature is mainly conserved within species, independent of organ, but also highly conserved between some organs between species. Using a similar approach, we also establish that the expression signature of in vitro senescence has similarity to that of in vivo aging in mouse but not human. This indicates differences between human and mouse aging and suggests that senescence might contribute to aging in mice but not humans.
A shared expression profile of aging in mammals
Included aging studies
Welle et al. 
HG U133 A+B
Old (n = 16) 66-77 years
Young (n = 15) 21-24 years
Lu et al. 
HG U95 Av2
Frontal cortex of brain
Old (n = 11) 73-106 years
Young (n = 10) 26-42 years
Rodwell et al. 
HG U133 A+B
Kidney cortex and medulla
Old (n = 28) 72-92 years
Young (n = 11) 27-45 years
Lee et al. 
Cerebellum and neocortex of brain
Old (n = 2/3*) 30 months
Young (n = 3) 5 months
Lee et al. 
MG U74 A
Old (n = 5) 30 months
Young (n = 5) 5 months
Blalock et al. 
RG U34 A
Old (n = 10) 24 months
Young (n = 9) 4 months
3c). The strongest determinant of similarity was species, as this was the only group demonstrating an average Z-score that was larger then zero (P = 0.02). Z-scores were preferred in the comparisons as they contain a significance estimate, although similar results were obtained when using correlation measures. The results regarding average correlations within species and tissues indicate a species-specific transcription profile in aging, although some organs show similarities across species (Figure 3c).
We also assessed whether the relationships differed between species and across tissues. We compared the correlations, with regard to their displacement (in standard deviations) from the distribution of simulated correlations, across both species and tissues (Figure
Comparing the transcriptomes of in vitrosenescence and aging
Included senescence studies
Induction of senescence
Zhang et al. 
Schwarze et al. 
Prostate epithelial cells
Zhang et al. 
Mammary epithelial cells
Larsson et al. 
MG U74 Av2
Thymus epithelial cells
Heat inactivation of SV40 large T antigen
Although these data would suggest that there are no similarities between senescence and aging, we suspected that this might reflect the models used in the senescence studies. In the case of the aging datasets, we compared young and old tissues, where extensive cell division does not take place. In the case of the senescence studies, however, the four senescent transcription profiles were generated by comparisons of logarithmic phase growing cultures to non-growing senescent cultures. We speculated, therefore, that the difference in terms of growth-associated transcriptional changes might obscure the senescence transcriptome signature. To clarify this we studied Gene Ontology categories and demonstrated an overrepresentation of categories related to proliferation in all senescence studies but not in the aging studies (data not shown). To test the hypothesis that cellular senescence is indeed similar to aging (if proliferation is not a component of the comparison), we used two of the datasets where non-growing quiescent cultures had been used to identify senescence specific - not growth-related - gene expression changes. The similarities between the senescent datasets that included a quiescent control were maintained after exchanging the proliferating control for the quiescent control (the correlation was 0.22 when compared to proliferating cells and 0.12 when compared to growth arrested cells).
Gene Ontology categories correlating in six comparisons of senescence and mouse aging data sets
Gene Ontology category
Aging mouse cerebellum
Aging mouse neocortex
Aging mouse heart
Aging mouse cerebellum
Aging mouse neocortex
Aging mouse heart
P = 0.00648
P = 0.019
P = 0.00262
P = 0.0135
P = 0.0222
P = 0.0284
P = 0.00701
P = 0.00293
P = 0.00394
Pre-mRNA splicing factor activity
P = 0.0229
P = 0.00199
P = 0.00533
P = 0.000848
Hydrogen ion transporter activity
P = 0.0134
P = 0.00532
P = 0.000814
P = 0.023
Mitotic sister chromatid segregation
P = 0.000947
P = 0.0462
P = 0.0103
Nuclear mRNA splicing via spliceosome
P = 0.022
P < 0.000001
P = 0.000918
P = 0.000124
P = 0.00589
P = 0.0199
P = 0.00545
P = 0.00641
P = 0.0449
P = 0.0469
P = 3.47e-06
P = 1.25e-05
Regulation of cell cycle
P = 1.74e-05
P = 0.00372
P = 0.0083
P = 0.0012
P = 0.0115
P = 0.00277
P = 0.000188
P = 0.00143
P < 0.000001
P = 0.0447
P < 0.000001
P < 0.000001
Negative regulation of cell proliferation
P = 0.00615
P = 0.00636
P = 0.00246
P = 0.0126
P = 0.00713
P = 0.0042
P = 0.0402
P = 0.0372
P = 0.0272
P = 0.0347
P = 0.0308
P = 0.0251
P = 0.00305
P = 0.0112
P = 4.89e-05
P = 0.0137
P = 0.00965
P = 5.49e-05
Respiratory gaseous exchange
P = 6.25e-05
P = 8.77e-06
P < 0.000001
This is the first microarray meta-analysis where the components of a complex biological phenotype are dissected into its in vitro parts. An expression signature was identified that is common to the aging transcriptomes of mouse, man and rat, and conserved mainly within species, but also within organs across species. Evidence of a similarity between the gene expression signatures of aging and cellular senescence was apparent in mouse but, interestingly, not found in human.
The fundamental tool used in this report was estimates of correlations between gene expression changes in aging and senescence. The validity and the accuracy of the significance estimations are, therefore, important. Here, we calculated a ratio of change in gene expression that concurs with the transition from one biological state to another. By comparing changes in gene expression rather than absolute measures several advantages were achieved. Effects on absolute measures of expression level - such as array type differences in ability to quantify transcript abundance, or species and tissue differences in background gene expression - were eliminated. As a consequence, biased correlations within species or between identical arrays are probably of limited concern. Also, all the expression data were used in the comparison, which would not have been the case if genes differentially expressed above a certain threshold had been compared. Thus, the changes in gene expression could be compared to changes in other experiments, regardless of species and the microarray platform used. The ratio approach that we used cannot be applied without considering the biological conditions studied. This is demonstrated by the fact that we were unable to establish similarity between senescence and aging when a gene expression ratio of senescent to proliferating cells was used to describe a senescent transcriptome signature. Similarity became evident when the proliferating cell data were replaced with quiescent cell data, that is, when cell-growth differences no longer obscured the comparison.
Distributions of random correlations were produced for each study-to-study comparison and used to assess the significance of the correlation in that comparison. It might be argued that such simulations underestimate the variance of the true distribution as they assume that all probes (or genes) will display similar changes in signal intensity. In a study  where genes were randomly paired within quantiles of overall hybridization intensity, however, the resulting distributions of correlations did not show a significantly increased variance. We consider the statistical methods applied in this study solid, and conclude that many of the correlations obtained in this study are highly unlikely to occur by chance. Yet, it is possible that circumstantial factors may have contributed to unexpected correlations. Any unintentional source of effects on gene expression, transcript abundance or hybridization efficiency that is unevenly distributed between the case and control groups, in two or more datasets, could promote 'false' correlations. The occurrence of such false positive correlations between processes with no known similarity seems to be a lesser issue, however, as a study in which a similar approach to identify similarities of changes in global expression was applied identified no correlations >0.06 when comparing 320 unrelated studies . The correlations of all our summaries that we consider significant were above 0.06. To assess a biological similarity we used several studies and averaged the correlation in the aging-aging, aging-senescence, and senescence-senescence comparison groups. This should further limit the risk of a false positive leading to a false conclusion. A non-significant correlation between processes that indeed are similar would be caused by variation in the individual studies. The main source of variation on the Affymetrix platform, when the normalization approach we used is applied, appears to be biological (compared to technical) . In our approach we carefully assessed the technical quality of the individual studies, thereby limiting the risk of a non-significant result because of technical variation. In addition, each study we used was well replicated and able to show similarities to at least one other study in our comparisons, indicating that the biological variation did not limit the appearance of positive correlations. A full understanding of the risk of false negative correlations would need a large, controlled study where the expression changes of a priori known similar processes were compared across, for example, laboratories, platforms and labelling protocols. In our study, the calculated P values for the group of expression patterns supported average similarity, indicating that the data should be considered robust. We claim that a valid result can be obtained if the unmodified and robust Monte Carlo simulation is used to approximate significance and the similarities obtained are present in repeated comparisons of similar biological states.
Using these methods, we observed a shared expression signature in mammalian aging. Furthermore, we found that aging is more different across species than across organs of the same species. This finding is intuitive, given the differences between the organisms in terms of life span. Humans have a lifespan that is about 40 times that of mice. It seems possible that the aging phenotype that is allowed to develop during an 80 year life span would be different from that produced during a two to three year life span. In this perspective, it might even seem surprising that some of the cross species but within tissue comparisons give rise to significant correlations. Our findings argue for complexity in the aging program: one component is species specific and is manifested in all organs; another component is organ specific and can give rise to correlation within organs, across species.
An interesting finding in the present study is the detection of a senescence-like component within the murine aging gene expression pattern. This correlation could either occur as a consequence of an accumulation of senescent cells or by a phenotypic change of normal cells towards a senescent phenotype. In either case, the tissue seems to become more similar to senescent cells as it ages. It may appear contradictory that mice, with a much shorter life span and significantly longer telomeres than humans, would display a senescent like phenotype, as senescence is believed to occur as a result of extended cell division and telomere erosion. Cellular stress is an alternative mechanism of senescence induction, however, and the overlap between aging and senescence that we observe could be caused by stress-induced senescence. There are indications in the literature that stress is a general pathway towards senescence; for example, human primary cells can senesce because of replicative exhaustion or cellular stress . A contribution of stress is also observed during RAS induced senescence. Over-expression of RAS seems to induce senescence in cells with high levels of p16 by further increasing the p16 expression. In contrast, cells with low levels of p16 only attain moderate levels of p16 after RAS over-expression and, therefore, do not senesce . Thus, senescence can be induced by both telomeres, via p53/p21, and cellular stress independent of telomere length, characterised by p16 over-expression. The signalling pathways of telomere induced senescence and stress induced senescence are believed to overlap. This is supported by our study where the transcriptional changes with senescence induced by telomere shortening and stress correlate highly. Furthermore, if stress-induced senescence occurs in vivo in mice, one would expect p16 to accumulate with age. This is supported by a recent study where p16 was identified as a biomarker of aging in mice in a range of organs . Taken together, these results indicate that aging mouse cells could be stressed in vivo and that senescence could occur as a consequence of stress.
Intriguingly, we were unable to establish any similarities between cellular senescence and human tissue aging. One might argue that the low availability of human senescence studies where senescent cells could be compared to quiescent cells limits our analysis. This is unlikely, however, as not only the mouse senescence study, but also the human senescence study, showed similarities to mouse aging. This indicates that the similarities between aging and senescence are robust and less dependent on species. When the quiescence controlled senescence studies and human aging studies were compared, only one comparison showed a significant similarity. Given the robustness of the senescence to mouse aging comparisons, and the overlap between mouse aging and human senescence, it seems likely that if senescence were a pronounced phenotype of human aging, it should have been observed with the technique used presently.
A possible explanation for the difference between the correlation between senescence and aging in mouse and humans, respectively, is the difference in sensitivity to reactive oxygen species. Mice produce more reactive oxygen species . Hence, reactive oxygen species may elicit more pronounced damage in mice than in humans. This idea is supported by the fact that human cells can proliferate indefinitely if telomerase is expressed  whereas mouse cells with telomerase eventually enter senescence.
Accumulation of reactive oxygen species could, therefore, induce expression of p16, which in turn might be part of the process of aging in mice. In humans the oxidative stress effect could be expected to be less pronounced and may not induce senescence. The contribution of cellular senescence in human aging may, therefore, be restricted to the telomere driven replicative senescence. One could speculate that senescence indeed occurs in vivo in humans, but to a much lesser extent compared to mice. The effect of the senescent cells on the overall function of the organ, and thereby their contribution to aging, would then be expected to be less important. Other possibilities are that apoptosis is the end point of telomere erosion in humans, and that senescence down-stream of telomere erosion is an in vitro phenomenon, or that senescence does not generally occur in tissues but only affects stem cells when these fail to escape stress induced senescence or are unable to maintain their telomeres. Regardless of the explanation, our data provide indications that the senescent phenotype does not dominate the aging tissues of humans.
By comparing changes in gene expression, it is possible to identify conditions that share expression profiles. We found a shared expression profile in aging mammals that seems to be more robust within the species tested but also pronounced within some organs across species. By comparing changes of gene expression that occur during in vitro senescence to the aging transcriptomes, senescence like changes in gene expression characterise mouse aging but we could not detect such a similarity in human tissue. Our data provide indications that the senescent phenotype does not dominate the aging tissues of humans.
Materials and methods
Data collection and selection
PubMed was screened for studies of senescence or aging where microarrays had been used to monitor global changes in gene expression. Only studies with data derived from human, mouse and rat were included. From the 34 publications that were considered, we were able to use data from 10 (excluded studies are listed in Additional data file 1). Studies using microarrays containing less then 6,000 features were excluded, as subsequent cross-array and cross-species comparisons would become less reliable due to relatively few data points. Other studies were excluded beause we were unable to obtain CEL files, robust multi-array average (RMA) normalized data, or dChip normalized data from Affymetrix studies, or raw data from cDNA-chip studies. Three studies were excluded because of insufficient data quality (details on reasons for exclusion of individual studies are listed in Additional data file 1). The data from the selected studies were either downloaded [18, 19] or obtained directly from the authors.
In several cases, we excluded samples from individual datasets if the correlations to replicate samples were dramatically different from other within-replicate correlations (details on data selection can be found in Additional data file 2).
A summary of the normalization methods used for each dataset can be found in Additional data file 2. We normalized all datasets but one, for which we obtained dChip normalized data. All Affymetrix data were normalized with RMA, all cDNA data were normalized using loess normalization with default settings in the Bioconductor package LIMMA.
Allocation to case and control groups
For subsequent analysis, data needed to be divided into aged and young groups and senescent and control groups. For the aging data, this was straightforward, and the groups are summarized in Table 1. Data on senescence were divided into senescent and control groups. The control group consisted of proliferating or quiescent cells. In the datasets by Zhang et al. [20, 21], samples from senescent and control populations had been hybridized to the same microarray so that the differential expression between senescent and control populations is represented by a single value (the data belong to a single group). In the dataset by Larsson et al. , a temperature sensitive inhibitor of senescence (SV40 large T antigen) was used. To subtract differential expression caused by the temperature shift, wild-type (with a temperature insensitive SV40 large T antigen) cells were used as a reference, yielding four groups instead of two: cells with or without the temperature sensitive senescence inhibitor, each divided into permissive temperature groups and restrictive temperature groups. We used data for the 0 h and 72 h time points because the other datasets contained samples from senescent cells harvested days after senescence induction. Additional information on specific datasets is available in Additional file 3.
We used the RESOURCERER tool  for annotation and cross-referencing of platforms. Individual features in all datasets were assigned unique and unambiguous identifiers, specifying the gene and its orthologs, thus making cross array and species comparisons possible. For the datasets using Affymetrix MU6500 and U74A chips (annotation information not available within the RESOURCERER tool), additional annotation information was downloaded  and used for cross-reference to the identifiers given by RESOURCERER (for details see Additional data file 3).
Calculation of differential gene expression measures
Normalized log2 expression measures were used to produce an expression ratio:
Ri = (1/nj × Σ xj) - (1/nk × Σ xk)
where i is a specific feature on a microarray and x represents the log2-transformed expression measures for that feature on all chips in the experiment. Values for x belong to the aging (or senescent; j) group, or to the young (or senescence control; k) group. nj and nk indicate the number of chips in the j and k groups. Thus, the geometric mean of the expression measures from the aging (or senescent) samples was divided by the geometric mean of the expression measures from the young (or control) samples, giving rise to a measure of differential expression with aging or senescence. This removes all expression differences produced by irrelevant variables (for example, chip type and species) because these variables will be identical in the two groups. Ratios from different features, probing the same gene, were averaged. This way of calculating a differential expression ratio was applied to all Affymetrix studies and to all cDNA microarray data where a common reference sample was used for competitive hybridization. In the Zhang et al. [20, 21] studies, senescent and control samples were competitively hybridized to the same chip, thus no individual expression measures for senescent or control groups can be derived. In this case, the average normalized expression measure for an individual feature across arrays in the experiment was considered an equivalent of Ri and only subjected to averaging within the same gene.
In the study by Larsson et al. , four groups were used for calculating differential expression, and the expression of the ratio was adapted:
Ri = ((1/nj × Σ xj) - (1/nk × Σ xk)) - ((1/nl × Σ xl) - (1/nm × Σ xm))
Where j, k, l and m correspond to the four groups: cells with a temperature sensitive (j, k) or insensitive (l, m) SV40 large T antigen divided into restrictive (j, l) and permissive (k, m) temperature groups.
Having annotated all datasets using one set of identifiers, and transformed all data for a single gene in an individual dataset to a ratio of differential expression, two or more datasets could be assessed with regards to their correlation. We used Pearson, Spearman and Kendall models of correlation for all comparisons, and they produced consistent results. We have chosen to present Pearson correlations only.
Monte Carlo simulations
To assess the significance of individual correlations, we iteratively (10,000 times) and randomly paired all genes between the two datasets compared, and calculated the Pearson, Spearman and Kendall correlations. Thus, a distribution of random correlations for each comparison was produced, serving as a reference for each correlation measure.
Calculation of Pvalues for correlation measures
The correlation values produced by the Monte Carlo simulations were normally distributed. We used values for means and standard deviations of these distributions to calculate P values for the actual correlation in each comparison. We used Student's t tests to assess the means of groups of comparisons, with regards to their P values and 95% confidence intervals.
Confluent cells as a reference for senescence
In the Zhang et al.  and Larsson et al.  studies, quiescent/confluent - rather than proliferating - cells had also been used to provide alternative controls for the senescent cell models. To remove effects of cell proliferation differences (between control data in the aging and senescence studies) on the gene expression signature of senescence, we used the data derived from quiescent cells. For the Zhang et al. dataset, where the senescence and control samples have been hybridized to the same array, we simply used the senescence versus quiescence data. For the Larsson et al. dataset, the log2 expression measures from the temperature permissive group was replaced by the expression measures from the confluent group, in the k entry of the ratio expression above.
Analysis of Gene Ontology categories giving rise to correlation between mouse aging and mouse and human senescence
We analysed the six comparisons of senescence (human and mouse) and mouse aging that showed average correlation. For all Gene Ontology categories, with more than two genes present in a specific comparison, we calculated the partial summation of the Pearson Correlation:
rj = Σ (xj - μx)(yj - μy)/nσxσy
where j is the set of ortholog pairs associated with that Gene Ontology category, μx, μy and σx, σy are the means and standard deviations of expression ratios for all genes or orthologs appearing in the comparison (the entire study to study comparison), and n is the total number of genes. Statistical significance was assessed by performing Monte Carlo simulations for individual Gene Ontology categories. Random pairing of genes was followed by random assignment to gene ontology categories, and repeated 10,000 times. Calculations of P values were performed as described above for the calculation of P values of correlations between entire datasets.
Additional data files
The following additional data are available with the online version of this paper. Additional data file 1 lists the excluded studies and the reasons for their exclusion. Additional data file 2 provides details on data selection and the normalization methods used for each data set. Additional data file 3 provides information on how data was allocated to case (aging or senescent) and control (young, proliferating or quiescent) groups as well as on the cross-referencing of microarrays.
- Hayflick L, Moorhead PS: The serial cultivation of human diploid cell strains. Exp Cell Res. 1961, 25: 585-621. 10.1016/0014-4827(61)90192-6.PubMedView ArticleGoogle Scholar
- Hayflick L: The limited in vitro lifetime of human diploid cell strains. Exp Cell Res. 1965, 37: 614-636. 10.1016/0014-4827(65)90211-9.PubMedView ArticleGoogle Scholar
- Baird DM, Kipling D: The extent and significance of telomere loss with age. Ann NY Acad Sci. 2004, 1019: 265-268. 10.1196/annals.1297.044.PubMedView ArticleGoogle Scholar
- Cristofalo VJ, Allen RG, Pignolo RJ, Martin BG, Beck JC: Relationship between donor age and the replicative lifespan of human cells in culture: a reevaluation. Proc Natl Acad Sci USA. 1998, 95: 10614-10619. 10.1073/pnas.95.18.10614.PubMedPubMed CentralView ArticleGoogle Scholar
- Dimri GP, Lee X, Basile G, Acosta M, Scott G, Roskelley C, Medrano EE, Linskens M, Rubelj I, Pereira-Smith O, et al: A biomarker that identifies senescent human cells in culture and in aging skin in vivo. Proc Natl Acad Sci USA. 1995, 92: 9363-9367.PubMedPubMed CentralView ArticleGoogle Scholar
- Severino J, Allen RG, Balin S, Balin A, Cristofalo VJ: Is beta-galactosidase staining a marker of senescence in vitro and in vivo?. Exp Cell Res. 2000, 257: 162-171. 10.1006/excr.2000.4875.PubMedView ArticleGoogle Scholar
- Chang S, Multani AS, Cabrera NG, Naylor ML, Laud P, Lombard D, Pathak S, Guarente L, DePinho RA: Essential role of limiting telomeres in the pathogenesis of Werner syndrome. Nat Genet. 2004, 36: 877-882. 10.1038/ng1389.PubMedView ArticleGoogle Scholar
- Schmitt CA, Fridman JS, Yang M, Lee S, Baranov E, Hoffman RM, Lowe SW: A senescence program controlled by p53 and p16INK4a contributes to the outcome of cancer therapy. Cell. 2002, 109: 335-346. 10.1016/S0092-8674(02)00734-1.PubMedView ArticleGoogle Scholar
- Satyanarayana A, Wiemann SU, Buer J, Lauber J, Dittmar KE, Wustefeld T, Blasco MA, Manns MP, Rudolph KL: Telomere shortening impairs organ regeneration by inhibiting cell cycle re-entry of a subpopulation of cells. EMBO J. 2003, 22: 4003-4013. 10.1093/emboj/cdg367.PubMedPubMed CentralView ArticleGoogle Scholar
- Krishnamurthy J, Torrice C, Ramsey MR, Kovalev GI, Al-Regaiey K, Su L, Sharpless NE: Ink4a/Arf expression is a biomarker of aging. J Clin Invest. 2004, 114: 1299-1307. 10.1172/JCI200422475.PubMedPubMed CentralView ArticleGoogle Scholar
- McCarroll SA, Murphy CT, Zou S, Pletcher SD, Chin CS, Jan YN, Kenyon C, Bargmann CI, Li H: Comparing genomic expression patterns across species identifies shared transcriptional profile in aging. Nat Genet. 2004, 36: 197-204. 10.1038/ng1291.PubMedView ArticleGoogle Scholar
- Timmons JA, Larsson O, Jansson E, Fischer H, Gustafsson T, Greenhaff PL, Ridden J, Rachman J, Peyrard-Janvid M, Sundberg CJ: Human skeletal muscle exercise transcriptome provides a novel perspective on Duchenne muscular dystrophy. FASEB J.Google Scholar
- Zakharkin SO, Kim K, Mehta T, Chen L, Barnes S, Scheirer KE, Parrish RS, Allison DB, Page GP: Sources of variation in Affymetrix microarray experiments. BMC Bioinformatics. 2005, 6: 214-10.1186/1471-2105-6-214.PubMedPubMed CentralView ArticleGoogle Scholar
- Herbig U, Jobling WA, Chen BP, Chen DJ, Sedivy JM: Telomere shortening triggers senescence of human cells through a pathway involving ATM, p53, and p21(CIP1), but not p16(INK4a). Mol Cell. 2004, 14: 501-513. 10.1016/S1097-2765(04)00256-4.PubMedView ArticleGoogle Scholar
- Benanti JA, Galloway DA: Normal human fibroblasts are resistant to RAS-induced senescence. Mol Cell Biol. 2004, 24: 2842-2852. 10.1128/MCB.24.7.2842-2852.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Ku HH, Brunk UT, Sohal RS: Relationship between mitochondrial superoxide and hydrogen peroxide production and longevity of mammalian species. Free Radic Biol Med. 1993, 15: 621-627. 10.1016/0891-5849(93)90165-Q.PubMedView ArticleGoogle Scholar
- Bodnar AG, Ouellette M, Frolkis M, Holt SE, Chiu CP, Morin GB, Harley CB, Shay JW, Lichtsteiner S, Wright WE: Extension of life-span by introduction of telomerase into normal human cells. Science. 1998, 279: 349-352. 10.1126/science.279.5349.349.PubMedView ArticleGoogle Scholar
- Gene Expression Omnibus. [http://www.ncbi.nih.gov/geo]
- Stanford MicroArray Database. [http://genome-www5.stanford.edu/]
- Zhang H, Herbert BS, Pan KH, Shay JW, Cohen SN: Disparate effects of telomere attrition on gene expression during replicative senescence of human mammary epithelial cells cultured under different conditions. Oncogene. 2004, 23: 6193-6198. 10.1038/sj.onc.1207834.PubMedView ArticleGoogle Scholar
- Zhang H, Pan KH, Cohen SN: Senescence-specific gene expression fingerprints reveal cell-type-dependent physical clustering of up-regulated chromosomal loci. Proc Natl Acad Sci USA. 2003, 100: 3251-3256. 10.1073/pnas.2627983100.PubMedPubMed CentralView ArticleGoogle Scholar
- Larsson O, Scheele C, Liang Z, Moll J, Karlsson C, Wahlestedt C: Kinetics of senescence-associated changes of gene expression in an epithelial, temperature-sensitive SV40 large T antigen model. Cancer Res. 2004, 64: 482-489. 10.1158/0008-5472.CAN-03-1872.PubMedView ArticleGoogle Scholar
- Animal Resourcerer. [http://www.tigr.org/tigr-scripts/magic/r1.pl]
- Affymetrix. [http://www.affymetrix.com]
- Lee CK, Weindruch R, Prolla TA: Gene-expression profile of the ageing brain in mice. Nat Genet. 2000, 25: 294-297. 10.1038/77046.PubMedView ArticleGoogle Scholar
- Welle S, Brooks A, Thornton CA: Senescence-related changes in gene expression in muscle: similarities and differences between mice and men. Physiol Genomics. 2001, 5: 67-73.PubMedGoogle Scholar
- Lu T, Pan Y, Kao SY, Li C, Kohane I, Chan J, Yankner BA: Gene regulation and DNA damage in the ageing human brain. Nature. 2004, 429: 883-891. 10.1038/nature02661.PubMedView ArticleGoogle Scholar
- Rodwell GE, Sonu R, Zahn JM, Lund J, Wilhelmy J, Wang L, Xiao W, Mindrinos M, Crane E, Segal E, et al: A transcriptional profile of aging in the human kidney. PLoS Biol. 2004, 2: e427-10.1371/journal.pbio.0020427.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee CK, Allison DB, Brand J, Weindruch R, Prolla TA: Transcriptional profiles associated with aging and middle age-onset caloric restriction in mouse hearts. Proc Natl Acad Sci USA. 2002, 99: 14988-14993. 10.1073/pnas.232308999.PubMedPubMed CentralView ArticleGoogle Scholar
- Blalock EM, Chen KC, Sharrow K, Herman JP, Porter NM, Foster TC, Landfield PW: Gene microarrays in hippocampal aging: statistical profiling identifies novel processes correlated with cognitive impairment. J Neurosci. 2003, 23: 3807-3819.PubMedGoogle Scholar
- Schwarze SR, DePrimo SE, Grabert LM, Fu VX, Brooks JD, Jarrard DF: Novel pathways associated with bypassing cellular senescence in human prostate epithelial cells. J Biol Chem. 2002, 277: 14877-14883. 10.1074/jbc.M200373200.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.