Research Paper Volume 5, Issue 8 pp 599—606

Serum from humans on long-term calorie restriction enhances stress resistance in cell culture

Daniela Omodei1,2, , Danilo Licastro3, , Francesco Salvatore2,4, , Seth D. Crosby5, , Luigi Fontana1,2,6, ,

  • 1 Division of Geriatrics and Nutritional Science, Washington University School of Medicine, St. Louis, MO 63130, USA
  • 2 CEINGE-Biotecnologie Avanzate scarl, Napoli, Italy
  • 3 CBM Scrl - Genomics, Area Science Park, Basovizza, Trieste, Italy
  • 4 Fondazione SDN-IRCCS, Napoli, Italy
  • 5 Department of Genetics, Washington University School of Medicine, St. Louis, MO 63130, USA
  • 6 Department of Medicine, Salerno University School of Medicine, Salerno; Italy
* Equal contribution

Received: July 10, 2013       Accepted: July 24, 2013       Published: July 27, 2013      

https://doi.org/10.18632/aging.100584
How to Cite

Copyright: © 2013 Omodei 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.

Abstract

Calorie restriction (CR) without malnutrition is the most robust intervention to slow aging and extend healthy lifespan in experimental model organisms. Several metabolic and molecular adaptations have been hypothesized to play a role in mediating the anti-aging effects of CR, including enhanced stress resistance, reduced oxidative stress and several neuroendocrine modifications. However, little is known about the independent effect of circulating factors in modulating key molecular pathways. In this study, we used sera collected from individuals practicing long-term CR and from age- and sex-matched individuals on a typical US diet to culture human primary fibroblasts and assess the effects on gene expression and stress resistance. We show that treatment of cultured cells with CR sera caused increased expression of stress-response genes and enhanced tolerance to oxidants. Cells cultured in serum from CR individuals showed a 30% increase in resistance to H2O2 damage. Consistently, SOD2 and GPX1 mRNA, two key endogenous antioxidant enzymes, were increased by 2 and 2.5 folds respectively in cells cultured with CR sera. These cellular and molecular adaptations mirror some of the key effects of CR in animals, and further suggest that circulating factors contribute to the CR-mediated protection against oxidative stress and stress-response in humans as well.

Introduction

Calorie restriction (CR) without malnutrition is the most potent and reproducible intervention for slowing aging and protecting against cancer in mammals [1,2]. Some of the beneficial effects of CR are mediated by down-regulation of inflammatory pathways and up-regulation of genes promoting stress resistance, DNA and cellular repair, and antioxidant enzymatic activity [3]. It has been hypothesized that changes in circulating factors may contribute, at least in part, to some of these beneficial effects of CR in rodents [4]. However, nothing is known on the molecular and cellular adaptations induced by circulating factors in men and women practicing long-term severe CR.

The importance of circulating factors in modulating health and longevity is supported by data obtained by heterochronic parabiosis experiments in mice [5]. Exposure of old mice to circulating factors present in young serum increased neurogenesis and regeneration after muscle injury [6,7]. It is well known that senescent stromal cells secrete a number of factors that can promote cancer and aging by increasing proliferation and tumorigenesis of epithelial cells, and by inhibiting stem cell function and antioxidant capacity [8]. It is our working hypothesis that CR by changing factors present in the systemic milieu (i.e. serum) can influence cell function by modulating key signaling pathways.

To test this hypothesis we have evaluated the effects of serum collected from volunteers practicing long-term CR on global gene expression and on stress resistance in cell culture experiments. In these in vitro experiments, we compared the molecular and cellular effects on human primary fibroblasts cultured of media supplemented with human sera obtained from: 1) healthy lean and weight-stable volunteers who were consuming CR diets, with more than adequate intake of protein and vitamins, for an average of 7 years, and 2) age-, sex-matched, non-obese individuals eating Western diets (WD).

Results

Gene expression analysis shows up-regulation of genes involved in NRF2-mediated oxidative stress response

Mammalian fibroblasts have been extensively used to assess global transcriptional changes induced by circulating factors (i.e. serum) on cell function [9]. To test the hypothesis that the CR-induced changes in circulating factors can influence cell function, we compared the transcriptional profiles of human primary fibroblasts treated with CR or WD serum. Differentially expressed genes were selected using a moderate t-test approach [10] and multiple test correction [11], and these genes were used to identify transcriptionally altered pathways. Ingenuity Pathways Analysis software (IPA) showed a significant association of a subgroup of 83 genes (Supplementary Table S1) with many oxidative response-related pathways (Fig. 1). In particular, NF-E2-related factor 2 (Nrf2)-mediated oxidative stress response contained 30 genes that were dysregulated in the treated fibroblasts (BH p-value 3.33E-02), and many of them were direct expression effectors of Nrf2 activity (Fig. 2). Meta-analysis performed using a competitive gene set test [12] pinpoints 11 molecular signatures significantly down-regulated in a ‘self-contained’ gene test (Table 2 red Highlights). To further investigate the result of microarray data and confirm the association with oxidative response-related pathways, we performed quantitative real-time PCR for the following genes which were shown by microarray to be detectable and significantly dysregulated: GSTK1, GPX1, SOD2, GSTT1, and IDH2. As shown in Figure 3, relative abundance of message for these genes was significantly higher in the CR group than the control group (p<0.0005) with a fold change difference ≥ 2.2-fold.

Pathway-based analysis. Circos diagram shows the relationship between CR differentially expressed genes and oxidative-related pathways. Ribbon ends represent links between genes and pathways while the width of the ribbon correlates with the number of genes involved. Segments in the outer ring indicate the total number of genes involved in the corresponding pathways while thin internal segments indicate the number of up-regulated (light green) and down-regulated (orange) genes.

Figure 1. Pathway-based analysis. Circos diagram shows the relationship between CR differentially expressed genes and oxidative-related pathways. Ribbon ends represent links between genes and pathways while the width of the ribbon correlates with the number of genes involved. Segments in the outer ring indicate the total number of genes involved in the corresponding pathways while thin internal segments indicate the number of up-regulated (light green) and down-regulated (orange) genes.

Connectivity map of NRF2-mediated oxidative stress response adapted from IPA software. Pathways analysis identified the NRF2-mediated oxidative stress response pathway with a statistically significant value (BH p-value 3.33E-02). Green colored shapes represent up-regulated genes while red ones represent down-regulated genes of fibroblast cell lines treated 48 hours with CR serum.

Figure 2. Connectivity map of NRF2-mediated oxidative stress response adapted from IPA software. Pathways analysis identified the NRF2-mediated oxidative stress response pathway with a statistically significant value (BH p-value 3.33E-02). Green colored shapes represent up-regulated genes while red ones represent down-regulated genes of fibroblast cell lines treated 48 hours with CR serum.

Table 1. Characteristics of the study subjects

CR group (n=12)WD group (n=13)P value
Age (yrs)57.8±8.455.7±8.7ns
Sex (M/F)10/211/2
Height (m)1.73±0.071.76±0.12ns
Weight (Kg)57.1±6.679.8±14.00.0002
BMI (kg/m2)19.1±1.225.7±2.3P<0.0001
Total body fat (%)9.8±4.825.9±5.6P<0.0001
Values are means ± SD

Table 2. Gene Set Enrichment CAMERA and ROAST results

Gene Set Enrichment results.Correlation Adjusted Mean Rank GeneRotation gene set testing for linear models
Molecular SignatureNGenesCorDownUpTwoSidedAdj.P.Value.DownAdj.P.Value.UpAdj.P.Value.MixedActive.Proportion.DownActive.Proportion.UpActive.Proportion.Mixed
PROTEIN_MATURATION170.0270.99760.00240.001.000.100.000.180.590.76
DNA_FRAGMENTATION_DURING_APOPTOSIS31−0.0020.99510.00490.011.000.070.000.130.420.55
AMYLOID_PRECURSOR_PROTEIN_METABOLIC_PROCESS240.0050.99480.00520.011.000.070.000.080.380.46
HYDROLASE_ACTIVITY_ACTING_ON_CARBON_NITROGE N_BUT_N0T_PEPT1DE_B0NDS_IN_LINEAR_AMIDES240.0050.99250.00750.021.000.100.000.130.290.42
MITOT1C_ CELL_ CYCLE_ CHECKPOINT450.0120.00760.99240.020.061.000.010.310.090.40
RAS_ GTPASE_ BINDING35−0.0100.00740.99260.010.041.000.010.340.140.49
HOMEOS TASIS_ OF_NUMBER_ OF_ CELLS36−0.0090.00690.99310.010.031.000.000.310.080.39
U12 DEPENDENT SPLICEOSOME10−0.0690.00630.99370.010.031.000.010.400.200.60
APICAL_PLASMA_MEMBRANE25−0.0040.00560.99440.010.041.000.040.240.040.28
Stress_response_CR110.0020.00520.99480.010.031.000.000.550.090.64
INSOLUBLE_FRACTION28−0.0030.00460.99540.010.031.000.010.290.040.32
DNA_DAMAGE_RESPONSE_SIGNAL_TRANSDUCTION_BY _P53_CLASS_MEDIAT0R31−0.0020.00430.99570.010.031.000.010.320.100.42
APICAL_PART_OF_CELL29−0.0040.00420.99580.010.031.000.060.240.030.28
SMALL_CONJUGA TING_PROTElN_BINDING21−0.0150.00400.99600.010.031.000.020.330.050.38
AMINOPEPTIDASE_ACTIVITY190.0050.00150.9985>0.010.031.000.000.470.000.47
UBIQUITIN_BINDING14−0.0160.00080.9992>0.010.031.000.010.500.070.57
Real Time qPCR analysis of GSTK1, GSTT1, IDH2, GPX1 and SOD2 expression in BJ cells cultured 48hs with CR and WD sera. (A-E) Data are expressed as the gene of interest abundance relative to the GAPDH (Mean ± SEM); n = 12 CR, n = 13 WD. Student's t test was used to determine P-values. (F) Data are expressed as fold changes relative to CR cultures.

Figure 3. Real Time qPCR analysis of GSTK1, GSTT1, IDH2, GPX1 and SOD2 expression in BJ cells cultured 48hs with CR and WD sera. (A-E) Data are expressed as the gene of interest abundance relative to the GAPDH (Mean ± SEM); n = 12 CR, n = 13 WD. Student's t test was used to determine P-values. (F) Data are expressed as fold changes relative to CR cultures.

Incubation with human CR serum increases stress resistance to hydrogen peroxide in vitro

To confirm the gene expression findings with a functional test, human primary fibroblasts cultured for 48 hours with CR or WD sera were exposed to increasing concentrations of the powerful oxidizing agent, hydrogen peroxide (H2O2) for 24 hours. As shown in Figure 4, treatment with H2O2 produced a dose-dependent decline in survival, but fibroblasts cultured with CR serum were significantly protected from the cytotoxic effect of the oxidant. At 600 μM about 67% of cells pretreated with CR serum were viable, while only 37% of cells pretreated with WD serum survived (p< 0.001), and even at higher H2O2 concentrations (800 μM) 32% of the cells pretreated with CR serum survived, while only 10% of the cells exposed to WD serum were still alive (p< 0.01).

Differential Stress Resistance to hydrogen peroxide. A human fibroblast cell line (BJ) was cultured 48-hrs with sera from CR and WD individuals. Cell viability (WST1 assay) was determined after a 24 h treatment with H2O2 (0 – 1000 M). Data are represented as means ± SEM. P- values were determined with 2way ANOVA with Bonferroni post test. ***P

Figure 4. Differential Stress Resistance to hydrogen peroxide. A human fibroblast cell line (BJ) was cultured 48-hrs with sera from CR and WD individuals. Cell viability (WST1 assay) was determined after a 24 h treatment with H2O2 (0 – 1000 M). Data are represented as means ± SEM. P- values were determined with 2way ANOVA with Bonferroni post test. ***P < 0.001; **P < 0.01

Discussion

We report here that some molecular and functional adaptations of CR, such enhanced stress resistance, can be mimicked in vitro by culturing cells with serum from individuals practicing long-term CR without malnutrition. Our findings show that the incubation of human fibroblasts with CR sera induces a transcriptional activation of several anti-oxidant genes along the Nrf2 pathway. In addition, we found that the incubation of fibroblasts with human CR serum is sufficient to increase stress resistance to hydrogen peroxide in vitro. These molecular and cellular adaptations are consistent with findings in CR rodents and monkeys, and may contribute at least in part to the protective anti-aging effects of CR [1,3].

In particular, accumulating data indicate that aging is associated with a down-regulation of stress-response genes and pathways [3]. To test the hypothesis that CR can counteract some of these age-associated detrimental effects in vitro, we first sought to determine if circulating factors of CR sera modify gene expression profiles of fibroblasts, which are the most abundant cell type in connective tissues and form the structural framework of tissues through their secretion of extracellular matrix components. To this end, we cultured human primary fibroblasts with sera from 12 middle-aged (57.8±8.4 yrs), weight-stable very lean (BMI=19.1±1.2 kg/m2) individuals who have been practicing ~30% CR with adequate nutrition (at least 100% of RDI for each nutrient), for 3-20 years, and a control group of 13 non-obese (BMI=25.7±2.3 kg/m2) age-matched sedentary controls eating a typical Western diet. Remarkably, we found that culturing cells with human CR sera for 48-hrs was sufficient to trigger a rapid transcriptional up-regulation of important genes and pathways that control stress resistance. In particular, we found that a highly significant number of transcripts along the Nrf2 pathway were altered by CR, and genes encoding important antioxidant enzymes, such as GSTK1, GSTT1, IDH2, GPX1 and SOD2, were up-regulated at least two fold in the CR group. Nrf2 is a key transcription factor that binds to the antioxidant response element (ARE) in response to environmental carcinogens or low levels of insulin, and increases the transcription of a variety of antioxidative and carcinogen-detoxification enzymes that boost protection against cellular and molecular damage and cancer [14-16]. It might be also possible that CR causes oxidative stress and therefore induce a NRF2-dependent response, supporting the hypothesis that oxidative damage is not a cause of aging [17-19].

In order to further investigate the role of CR sera in modifying the resistance of cells to oxidative stress, we next performed a series of experiments in which we exposed human primary fibroblasts to increasing concentrations of H2O2, a powerful oxidizing agent. Our data indicate that incubating fibroblasts with CR sera exerts a rapid and profound protection against hydrogen peroxide-induced cytotoxicity. These findings are consistent with data from studies in mice, suggesting that these selected cellular adaptations induced by circulating factors are conserved among species [4]. Indeed, most of the metabolic and hormonal parameters examined in CR humans are similar to those reported in CR rodents [20-22].

In conclusion, our findings reported here show that enhanced stress responsiveness known to occur in animals on CR following oxidative stress can be reproduced in vitro by culturing cells using serum from CR patients. Our findings provide support for a possible neuroendocrine based mechanism of CR and offer a model for the screening and investigation of CR mimetic agents. We anticipate that this novel in vitro model will be instrumental in the elucidation of biochemical, molecular and cellular mechanisms underlying the anti-aging effects of CR in humans.

Materials and Methods

Human subjects

Human serum samples were collected in the morning after an overnight fast from twelve individuals of the Calorie Restriction Society who had been practicing severe CR with adequate nutrition for an average of 7 years (range 3-20 years). Thirteen healthy, non-obese and sedentary (regular exercise < 1 h per week) individuals eating typical Western Diets (WD) were used as a sedentary comparison group; they were matched with the CR group in terms of age, sex and height. Whole blood was collected into sterile vacutainer tubes, allowed to clot for 30 minutes and then centrifuged. The serum was then stored at −80°C. The characteristics of the participants are shown in Table 1. All of the study participants were weight stable, i.e. < 2 kg weight change in the preceding 6 months. None of the participants had evidence of chronic disease (including cardiovascular, lung, gastrointestinal, type 2 diabetes and cancer). None were smokers. This study was approved by the Human Studies Committee of Washington University School of Medicine, and all participants gave informed consent before their participation.

Cell culture

The BJ human fibroblast cell line was obtained from ATCC and used in these experiments. BJ cells were cultured in 75 cm2 culture flasks containing basic culture medium: 86% of E-MEM medium supplemented with 10% fetal bovine serum (Sigma-Aldrich, St. Louis, MO, USA), 1% Sodium Pyruvate, 1% non-essential amino acids, 1% GlutaMAX and 1% antibiotics (GIBCO, Carlsbad, CA, USA) until the time of treatment. Cells were cultured in a humidified incubator at 37 °C and 5% CO2. For the experiments, FBS was replaced with human serum from CR and WD participants. Both FBS and human serum were heat inactivated for 30 minutes at 55 °C.

Microarray analysis

BJ cells were seeded at a concentration of 4.000 cells/cm2 in 25 T-25 flasks; when they reached 60% of confluence, the cells were starved for 24h without FBS. After starvation, 10% of human serum from 12 CR and 13 WD individuals was added to the cultures for 48h. Whole cells were detached from the flask with trypsin/EDTA (Sigma-Aldrich, St. Louis, MO, USA) and Total RNA was extracted using the mirVANA (Ambion, Inc, Austin, TX, USA) according to manufacturer's protocol. Integrity of total RNA was evaluated using capillary electrophoresis (Bioanalyzer 2100, Agilent technologies, Santa Clara, CA, USA) and quantified using a Nanodrop 1000 (Nanodrop, Wilmington, DE). Aliquots of RNA (400 ng) samples were amplified according to the specifications of the Illumina® TotalPrep™ RNA Amplification Kit (Ambion, Austin, TX, USA) to produce a pool of biotin-labeled RNA corresponding to the polyadenylated (mRNA) fraction. The cRNA samples were applied to the arrays of Illumina whole-genome HumanHT-12 v 3.0 (Illumina, San Diego, CA, USA) and hybridized according to the manufacturer's specification. Each array on the BeadChip targets over 25.000 annotated genes using 3-micron beads bearing covalently attached 50-base oligonucleotide probes. Each probe interrogates a single gene, and each bead type is represented with an average 15-fold redundancy on every array. The BeadChips were scanned with the Illumina's Beadarray system scanner (Illumina, San Diego, CA, USA). The hybridization images signal intensity was extracted and background subtracted using Illumina Inc. BeadStudio software version 3.3.7. The produced data were checked for the Illumina internal quality control and loaded into Bioconductor software [23,24] for statistical analysis using lumi package [25] (see below). The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus [26] and are accessible through GEO Series accession number GSE41790.(http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=nnonbsicqouiilw&acc=GSE41790).

Real Time PCR

Real-time PCR was used to verify microarray data. Total RNA was DNase I treated and reverse transcribed using the High Capacity RNA-to-cDNA kit (Applied Biosystems, Carlsbad, CA, USA). The cDNA was then used to evaluate the relative expression of select genes including SOD2, GPX1, IDH2, GSTK1 and GSTT1 using specific TaqMan Gene Expression Assays (Applied Biosystems, Carlsbad, CA, USA), which are pre-designed and pre-optimized gene-specific probe sets. DNA amplification was performed using the Applied Biosystems 7500 Fast Real-Time PCR system according to the manufacturer's protocol. The relative amount was calculated from the threshold cycles with the instrument's software (SDS 2.0, Applied Biosystems, Carlsbad, CA, USA), according to the manufacturer's instructions. We used the comparative cycle threshold (Ct) method (2 - CT) to calculate relative gene expression under experimental and control conditions normalized to GAPDH. The results are expressed either as the abundance of the gene of interest relative to the house keeping gene (2 - CT), or as fold change over control values (2 - CT).

Cell Proliferation assays in response to Hydrogen Peroxide

To determine the Differential Stress Resistance to hydrogen peroxide of serum from CR and WD individuals, human fibroblasts were seeded in 96-well plates at a concentration of 3000 cells/well in 86% D-MEM low glucose w/o phenol red supplemented with 1% Sodium Pyruvate, 1% non-essential amino acids, 1% GlutaMAX and 1% antibiotics (GIBCO, Carlsbad, CA, USA). After cells had seeded, media was removed and replaced with media containing 10% of human serum from 11 CR and 11 WD individuals was added in triplicate to the plates for 48hs. After incubation, freshly prepared hydrogen peroxide was added to the plates at a concentration of 200, 400, 600, 800, 1000 μM and cells were challenged for 24hs.

Cell viability was determined by the addition of 20 μl of the ready-to-use WST1 reagent (Roche diagnostics, Mannheim, Germany) to each well; the cells were incubated for approximately 3 hours at 37 °C and 5% CO2 in a humidified incubator. After this incubation period, the absorbance of the formazan dye formed was measured at a wavelength of 450 nm using a microplate spectrometer. The measured absorbance directly correlates to the number of viable cells.

Statistical analysis

Microarray data produced were quality checked with arrayQualityMetrics [24] package under R version 2.15.1 Bioconductor software [23] and undetectable and low detection value probes were filtered out from further analysis accordingly to the lumi package standards [25]. To identify differentially expressed genes, based on a moderate t-test, the Limma package [10] was used and genes were selected based on a p-value cut-off <0.05 (after FDR adjustment) and absolute fold change >1 following the Macroarrays Quality Control Program criteria [11]. IPA software was then used (www.ingenuity.com) to detect the enrichment of biofunctions and networks in the resulting list. Core analysis was performed using the following settings: Reference set: HumanHT-12 v 4.0 [v3 is mentioned in the methods]; Relationship to include: Direct and Indirect; Includes Endogenous Chemicals; Filter Summary: Consider only molecules and/or relationships where (species = Uncategorized (e.g. chemicals) OR Human) AND (confidence = Experimentally Observed). Canonical pathways analysis has been performed using as scoring method B-H Multiple Testing Correction p-values set to a threshold value of 0.05.

Meta-analysis was performed using the CAMERA function implemented in the Bioconductor Limma package for a ‘competitive gene set test’ [12]. The resulting top Molecular signatures have been further analyzed using a ‘focused’ gene set test ROAST [27]. All the list of genes organized as Signature and present in The Molecular Signature Database C2 collection (MSigDB) (www.broadinstitute.org/gsea/msigdb/) were updated to official symbols using the illumina annotation file HumanHT-12_V3_0_R2_11283641_A, while the CR signatures were created according to the Table 2 of 3 and M.musculus to H.sapiens conversion obtained using the NCBI HomoloGene tool (http://www.ncbi.nlm.nih.gov/homologene).

Representation of differentially expressed genes in their Genomic localization and visual representation of their relationships with the different pathways were produced using CIRCOS software [28]. Statistical analysis of real time PCR and cell proliferation data was performed using Student's t test to compare the transcripts relative abundance in the CR group and WD group. P values <0.05 were considered statistically significant.

Acknowledgments

This research is supported by grants from the Bakewell Foundation, AFAR, Glenn Foundation, the Longer Life Foundation (an RGA/Washington University Partnership), the Scott and Annie Appleby Charitable Trust, the National Center for Research Resources (UL1 RR024992), the National Institute of Diabetes And Digestive And Kidney Diseases (P30DK056341), the Italian Ministry of University and Research (PON02_00677, BIOGENE, potenziamento lab.8, 2012), and the Regione Campania (POR Campania FSE 2007/2013, CAMPUS-Biofreme). We thank the Genome Technology Access Center in the Department of Genetics for help with genomic analysis. The Center is partially supported by NCI Cancer Center Support Grant #P30 CA91842 to the Siteman Cancer Center and by ICTS/CTSA Grant# UL1RR024992 from the National Center for Research Resources (NCRR), a component of the National Institutes of Health (NIH), and NIH Roadmap for Medical Research. This publication is solely the responsibility of the authors and does not necessarily represent the official view of NCRR or NIH.

Conflicts of Interest

We declare that we have no conflicts of interest.

References

  • 1. Anderson RM, Shanmuganayagam D, Weindruch R. Caloric restriction and aging: studies in mice and monkeys. Toxicol Pathol. 2009; 37:47-51. [PubMed]
  • 2. Fontana L, Partridge L, Longo VD. Extending healthy life span--from yeast to humans. Science. 2010; 328:321-326. [PubMed]
  • 3. Lee CK, Klopp RG, Weindruch R, Prolla TA. Gene expression profile of aging and its retardation by caloric restriction. Science. 1999; 285:1390-1393. [PubMed]
  • 4. de Cabo R, Furer-Galban S, Anson RM, Gilman C, Gorospe M, Lane MA. An in vitro model of caloric restriction. Exp Gerontol. 2003; 38:631-639. [PubMed]
  • 5. Conboy IM and Rando TA. Heterochronic parabiosis for the study of the effects of aging on stem cells and their niches. Cell Cycle. 2012; 11:2260-2267. [PubMed]
  • 6. Villeda SA, Luo J, Mosher KI, Zou B, Britschgi M, Bieri G, Stan TM, Fainberg N, Ding Z, Eggel A, Lucin KM, Czirr E, Park JS, Couillard-Després S, Aigner L, Li G, et al. The ageing systemic milieu negatively regulates neurogenesis and cognitive function. Nature. 2011; 477:90-94. [PubMed]
  • 7. Conboy IM, Conboy MJ, Wagers AJ, Girma ER, Weissman IL, Rando TA. Rejuvenation of aged progenitor cells by exposure to a young systemic environment. Nature. 2005; 433:760-764. [PubMed]
  • 8. Campisi J, Andersen JK, Kapahi P, Melov S. Cellular senescence: a link between cancer and age-related degenerative disease? Semin Cancer Biol. 2011; 21:354-359. [PubMed]
  • 9. Iyer VR, Eisen MB, Ross DT, Schuler G, Moore T, Lee JC, Trent JM, Staudt LM, Hudson J Jr, Boguski MS, Lashkari D, Shalon D, Botstein D, Brown PO. The transcriptional program in the response of human fibroblasts to serum. Science. 1999; 283:83-87. [PubMed]
  • 10. Wettenhall JM, Simpson KM, Satterley K, Smyth GK. affylmGUY: a graphical user interface for linear modeling of single channel microarray data. Bioinformatics. 2006; 22:897-899. [PubMed]
  • 11. Shi L, Reid LH, Jones WD, Shippy R, Warrington JA, Baker SC, Collins PJ, de Longueville F, Kawasaki ES, Lee KY, Luo Y, Sun YA, Willey JC, Setterquist RA, Fischer GM, Tong W, et al. The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements. Nat Biotechnol. 2006; 24:1151-1161. [PubMed]
  • 12. Wu D and Smyth GK. Camera: a competitive gene set test accounting for inter-gene correlation. Nucleic Acids Res. 2012; 40:e133 [PubMed]
  • 13. Finkel T and Holbrook NJ. Oxidants, oxidative stress and the biology of ageing. Nature. 2000; 408:239-247. [PubMed]
  • 14. Kannan S and Jaiswal AK. Low and high dose UVB regulation of transcription factor NF-E2-related factor 2. Cancer Res. 2006; 66:8421-8429. [PubMed]
  • 15. Kim SK and Novak RF. The role of intracellular signaling in insulin-mediated regulation of drug metabolizing enzyme gene and protein expression. Pharmacol Ther. 2007; 113:88-120. [PubMed]
  • 16. Kwak MK, Itoh K, Yamamoto M, Kensler TW. Enhanced expression of the transcription factor Nrf2 by cancer chemopreventive agents: role of antioxidant response element-like sequences in the nrf2 promoter. Mol Cell Biol. 2002; 22:2883-2892. [PubMed]
  • 17. Gems D and Partridge L. Genetics of longevity in model organisms: debates and paradigm shifts. Annu Rev Physiol. 2013; 75:621-644. [PubMed]
  • 18. Gems D and de la Guardia Y. Alternative Perspectives on Aging in Caenorhabditis elegans: Reactive Oxygen Species or Hyperfunction? Antioxid Redox Signal. 2013; 19:321-329. [PubMed]
  • 19. Blagosklonny MV. Answering the ultimate question “what is the proximal cause of aging?”. Aging (Albany NY). 2012; 4:861-277. [PubMed]
  • 20. Cangemi R, Friedmann AJ, Holloszy JO, Fontana L. Effects of long-term calorie restriction on serum sex hormones concentration in men. Aging Cell. 2010; 9:236-242. [PubMed]
  • 21. Fontana L and Klein S. Aging, adiposity and calorie restriction. JAMA. 2007; 297:986-994. [PubMed]
  • 22. Fontana L, Klein S, Holloszy JO. Effects of long-term calorie restriction and endurance exercise on glucose tolerance, insulin action and adipokine production. Age. 2010; 32:97-108. [PubMed]
  • 23. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004; 5:R80 [PubMed]
  • 24. Ritchie ME, Diyagama D, Neilson J, van Laar R, Dobrovic A, Holloway A, Smyth GK. Empirical array quality weights in the analysis of microarray data. BMC Bioinformatics. 2006; 7:261 [PubMed]
  • 25. Du P, Kibbe WA, Lin SM. lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008; 24:1547-1548. [PubMed]
  • 26. Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002; 30:207-210. [PubMed]
  • 27. Wu D, Lim E, Vaillant F, Asselin-Labat ML, Visvader JE, Smyth GK. ROAST: rotation gene set tests for complex microarray experiments. Bioinformatics. 2010; 26:2176-2182. [PubMed]
  • 28. Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA. Circos: an information aesthetic for comparative genomics. Genome Res. 2009; 19:1639-1645. [PubMed]