Gene expression profiling of long-lived dwarf mice: longevity-associated genes and relationships with diet, gender and aging

BMC Genomics20078:353

DOI: 10.1186/1471-2164-8-353

Received: 30 May 2007

Accepted: 03 October 2007

Published: 03 October 2007

Abstract

Background

Long-lived strains of dwarf mice carry mutations that suppress growth hormone (GH) and insulin-like growth factor I (IGF-I) signaling. The downstream effects of these endocrine abnormalities, however, are not well understood and it is unclear how these processes interact with aging mechanisms. This study presents a comparative analysis of microarray experiments that have measured hepatic gene expression levels in long-lived strains carrying one of four mutations (Prop1 df/df, Pit1 dw/dw, Ghrhr lit/lit, GHR-KO) and describes how the effects of these mutations relate to one another at the transcriptional level. Points of overlap with the effects of calorie restriction (CR), CR mimetic compounds, low fat diets, gender dimorphism and aging were also examined.

Results

All dwarf mutations had larger and more consistent effects on IGF-I expression than dietary treatments. In comparison to dwarf mutations, however, the transcriptional effects of CR (and some CR mimetics) overlapped more strongly with those of aging. Surprisingly, the Ghrhr lit/lit mutation had much larger effects on gene expression than the GHR-KO mutation, even though both mutations affect the same endocrine pathway. Several genes potentially regulated or co-regulated with the IGF-I transcript in liver tissue were identified, including a DNA repair gene (Snm1) that is upregulated in proportion to IGF-I inhibition. A total of 13 genes exhibiting parallel differential expression patterns among all four strains of long-lived dwarf mice were identified, in addition to 30 genes with matching differential expression patterns in multiple long-lived dwarf strains and under CR.

Conclusion

Comparative analysis of microarray datasets can identify patterns and consistencies not discernable from any one dataset individually. This study implements new analytical approaches to provide a detailed comparison among the effects of life-extending mutations, dietary treatments, gender and aging. This comparison provides insight into a broad range of issues relevant to the study of mammalian aging. In this context, 43 longevity-associated genes are identified and individual genes with the highest level of support among all microarray experiments are highlighted. These results provide promising targets for future experimental investigation as well as potential clues for understanding the functional basis of lifespan extension in mammalian systems.

Background

Long-lived dwarf strains of laboratory mice offer promising tools for advancing our understanding of aging mechanisms and the basis of extended lifespan in mammals. Dwarf models have been found to exhibit an average lifespan increase of more than 50% compared to wildtype control strains, and thus represent genetic manipulations with an impact on longevity that is comparable to calorie restricted diets [1]. In addition to lifespan extension, long-lived dwarf strains exhibit superior health at advanced ages [2], and are less susceptible to age-related declines in memory, learning ability, and locomotion [3]. Dwarf mice experience lower incidence of kidney disease, cataracts and joint disease, as well as fatal neoplastic disease, such as lymphoma and adenocarcinoma [46]. Earlier studies have also suggested that development of transplanted tumors is impeded in dwarf mice, and cancer rates are reduced following exposure to chemical carcinogens [7, 8]. The loss of immune function and progression of collagen cross-linking that normally occurs with advancing age is diminished in some dwarf strains [9], and at the cellular level, fibroblasts of dwarf mice are resistant to several forms of stress, including oxidative stress, ultraviolet light, toxic metals and heat [10, 11]. Aside from improving our knowledge of aging mechanisms, therefore, understanding the unique features of dwarf mice may provide insight into a broad range of mechanisms relevant to health and disease-prevention in mammals.

Long-lived dwarf mouse strains carry mutations that suppress the growth hormone (GH)/insulin-like growth factor I (IGF-I) endocrine pathway [1]. In Ames (Prop-1 df/df) and Snell (Pit-1 dw/dw) mutants, for instance, development of the anterior pituitary is inhibited completely. This results in reduced GH and IGF-I levels, as well as other hormonal abnormalities, such as deficiency of thyroid stimulating hormone and prolactin [12, 13]. Ames and Snell dwarfs exhibit considerably reduced body size, but also have dramatic lifespan extension of 40–69% on average. In the Little (Ghrhr lit/lit) mouse, lifespan extension is less substantial (23–25%), and only occurs when mice are provided a low fat diet [9]. Hormonal abnormality of Little mice is limited to circulating GH levels, which are reduced because the lit/lit mutation renders the pituitary unresponsive to GH releasing hormone [14]. It is somewhat surprising that genetic alterations downstream of the lit/lit mutation promote a larger lifespan increase than that observed in little mice. In particular, GH-receptor knock-out mice (GHR-KO) have elevated serum GH levels, but are growth hormone resistant and exhibit considerable lifespan increases of 38–55% on average [15]. Taken together, these findings from four long-lived dwarf models firmly establish an endocrine basis of lifespan extension in mammals. Nevertheless, it remains unclear why such endocrine abnormalities affect lifespan, since inhibition of GH/IGF-I signaling is associated with extensive downstream effects, and only a fraction of these effects may be linked to mechanisms of aging [16].

Identification of genes affected by GH and IGF-I suppression is an important step towards understanding how this axis impacts longevity. Microarray studies of long-lived dwarf mice have been especially useful in this regard [1720]. Previous studies have focused primarily on liver tissue, since this is the major manufacturing site of IGF-I. Early investigations used cDNA arrays to detect expression differences between dwarf and normal mice, and found that 13 of 265 surveyed transcripts (4.9%) were differentially expressed in Ames dwarfs [17], while 60 of 2352 (2.5%) surveyed transcripts were differentially expressed in the Snell model [18]. Subsequent studies have used Affymetrix oligonucleotide arrays to screen a larger fraction of the genome. For instance, Amador-Noguez et al. [19] found that approximately 1100 of 14-thousand Affymetrix probesets (8.1%) were differentially expressed in Ames and Little dwarfs, respectively, while 547 (4.0%) transcripts were differentially expressed in both models. Boyleston et al. [20] performed a similar study using Ames and Snell dwarf mice, but focused only on probesets differentially expressed at every age group examined (6 – 24 months). This approach highlighted 785 such probesets in the Ames dwarf (1.9% of those surveyed), along with 205 probesets (1.7%) in the Snell dwarf, with 49 probesets satisfying the criterion in both long-lived models.

Previous microarray studies of long-lived dwarf mice have used varying statistical methodologies and criteria for identifying candidate genes, which complicates comparisons among studies based upon published results. Such comparisons would be useful, since genes differentially expressed with respect to multiple long-lived dwarf models are especially likely to play a role in aging and lifespan determination [17, 18]. Endocrine abnormalities differ to some degree among dwarf models, but presumably, lifespan extension in each model results from shared alterations that affect the GH/IGF-I pathway. Genes underlying extended lifespan should therefore be identified with respect to all long-lived models, while genes involved in pathways unrelated to lifespan extension should be specific to particular models. Furthermore, from a statistical perspective, a fraction of genes identified with respect to any one long-lived model in a single study are expected to be false-positives [21]. However, since P-values combine multiplicatively when results are pooled across studies [22], genes differentially expressed in more than one study are less likely to be false positive identifications. From biological and statistical standpoints, therefore, it is desirable to utilize data from independent studies to identify genes commonly induced among multiple long-lived dwarf models.

This study presents a comparative analysis of microarray experiments that have examined hepatic gene expression differences between long-lived dwarf mice and normal controls. The results provide a detailed view of how the effects of different dwarf mutations relate to one another in terms of gene expression, and describe how these effects relate to those associated with aging, gender dimorphism, low fat diets, caloric restriction (CR) and several different CR mimetic compounds. A main goal of the analysis was to identify longevity-associated genes with the highest level of support based on currently available microarray data. Applying consistent statistical methodology, therefore, genes exhibiting parallel differential expression patterns in four long-lived dwarf models were identified (Ames, Snell, Little and GHR-KO). Since mechanisms of lifespan extension in long-lived dwarf models may overlap with those conferring increased longevity under CR [2325], genes exhibiting parallel transcriptional changes in multiple long-lived models and under CR were also identified. Further steps were taken to evaluate the role of each candidate gene with respect to aging and longevity, and to determine which genes provided the most promise for future experimental investigation.

Results

Datasets and IGF-I expression patterns

Expression data was examined from four dwarf models that have previously been associated with significant lifespan extension in laboratory studies, including the Ames (Prop1 df/df), Snell (Pit1 dw/dw), Little (Ghrhr lit/lit), and GHR-KO (GHR-KO) mutants. Genetic alterations with effects potentially related to those of life-extending dwarf mutations were also considered. For instance, the effects of two GHR knock-in models with disruptions of GHR receptor residues 391 or 569 were examined [26]. Gene expression levels of the B6.C3H-6T congenic mouse were considered because this model exhibits a 30–40% reduction in serum IGF-I levels [27]. The effects of gender were examined because female mice live longer than male mice, and also because one study has suggested that dwarf mutations induce a masculine-to-feminine shift in gene expression patterns [20]. Dietary treatments known or postulated to increase lifespan were considered for the purpose of comparison. These treatments included low fat diets, CR and several potential CR mimetic compounds (metformin, glipizide, rosiglitazone, and soy isoflavone extract). The potential CR mimetic compounds are known to either influence insulin sensitivity (metformin, glipizide, rosiglitazone) or suppress tumorigenesis (soy isoflavone extract) [28]. Metformin treatment of female mice has been found to increase mean and maximum lifespan by 8% and 13%, respectively [29].

The analysis is based on a series of contrasts comparing gene expression levels between two experimental treatments. Table 1 provides descriptions and representative symbols for all contrasts examined. Each contrast involved a pair of treatments designated A and B. The A treatment is a genotype or dietary manipulation known or postulated to be associated with increased longevity. The corresponding B treatment is an appropriate control that permitted evaluation of how treatment A affects gene expression. The analysis was based on 8525 probesets that could be reliably matched among all three platforms (using Affymetrix best match tables). For each contrast, the null hypothesis of H 0: μ Ai = μ Bi was tested for each probeset, where μ Ai and μ Bi represent the mean expression levels of gene i in treatments A and B, respectively (see Methods).
Table 1

Treatment Contrasts

Contrast Symbol

Treatment A

Treatment B

snell5a

Pit1dw/dw males, age 4–6 months, n = 4

Pit1dw/? control males, age 4–6 months, n = 4

snell25a

Pit1dw/dw males, age 24–26 months, n = 3

Pit1dw/? control males, age 24–26 months, n = 3

ames5Ab

Prop1df/df males, age 4–6 months, n = 5

Prop1+/+ control males, age 4–6 months, n = 5

ames13Ab

Prop1df/df males, age 12–14 months, n = 5

Prop1+/+ control males, age 12–14 months, n = 5

ames25Ab

Prop1df/df males, age 24–27 months, n = 5

Prop1+/+ control males, age 24–27 months, n = 5

ames3Bc

Prop1df/df males, age 3 months, n = 3

Prop1+/+ control males, age 3 months, n = 5

ames6Bc

Prop1df/df males, age 6 months, n = 3

Prop1+/+ control males, age 6 months, n = 5

ames12Bc

Prop1df/df males, age 12 months, n = 3

Prop1+/+ control males, age 12 months, n = 5

ames24Bc

Prop1df/df males, age 24 months, n = 3

Prop1+/+ control males, age 24 months, n = 5

little3c

Ghrhrlit/lit males, age 3 months, n = 3

Ghrhrlit/+ males, age 3 months, n = 3

little6c

Ghrhrlit/lit males, age 6 months, n = 3

Ghrhrlit/+ males, age 6 months, n = 3

little12c

Ghrhrlit/lit males, age 12 months, n = 3

Ghrhrlit/+ males, age 12 months, n = 3

little24c

Ghrhrlit/lit males, age 24 months, n = 3

Ghrhrlit/+ males, age 24 months, n = 3

GHR-KOd

GHR(-/-) males, age 42 days, n = 3

wild type males, age 42 days, n = 3

GHR-KI1d

GHR knock-in mutant 569, males, age 42 days, n = 3

wild type males, age 42 days, n = 3

GHR-KI2d

GHR knock-in mutant 391, males, age 42 days, n = 3

wild type males, age 42 days, n = 3

B6e

C57BL/6J (B6) females, age 2 months, n = 3 (20–30% reduced serum IGF-I)

C3H/HeJ (C3H) females, age 2 months, n = 3

Genderf

wild type females, age 3–6 months, n = 6

wild type males, age 3–6 months, n = 6

cr(2,6)g

4 months of 30% calorie restriction, initiated age 2, tissue harvest age 6, wild type females, n = 8

wild type littermates of treatment A fed on control diet, n = 7

cr(2,6)dfg

4 months of 30% calorie restriction, initiated age 2, tissue harvest age 6, Prop1 df/df females, n = 8

Prop1df/df female littermates of treatment A with control diet, n = 8

cr(20,22)h

2 months of 40% calorie restriction, initiated age 20, tissue harvest age 22, strain B6C3F1 males, n = 4

strain B6C3F1 males fed on control diet, tissue harvest age 22, n = 4

cr(5,22)h

17 months 40% calorie restriction, initiated age 5, tissue harvest age 22, strain B6C3F1 males, n = 4

strain B6C3F1 males fed on control diet, tissue harvest age 22, n = 4

meth

metformin, 2100 mg per kg diet, tissue harvest age 22, strain B6C3F1 males, n = 4

control diet, tissue harvest age 22, strain B6C3F1 males, n = 4

met(db/db)i

metformin, 400 mg/kg, C57BL/ksj – db/db males, n = 5

placebo, C57BL/ksj – db/db males, n = 5

gliph

glipizide, 1050 mg per kg diet, tissue harvest age 22, strain B6C3F1 males, n = 4

control diet, tissue harvest age 22, strain B6C3F1 males, n = 4

gmh

metformin & glipizide, met dose: 1050 mg per kg diet, glip dose: 525 mg per kg diet, tissue harvest age 22, strain B6C3F1 males, n = 4

control diet, tissue harvest age 22, strain B6C3F1 males, n = 4

rosh

rosiglitazone, 80 mg per kg diet, tissue harvest age 22, strain B6C3F1 males, n = 4

control diet, tissue harvest age 22, strain B6C3F1 males, n = 4

soyh

soy isoflavone extract, 0.25%, tissue harvest age 22, strain B6C3F1 males, n = 4

control diet, tissue harvest age 22, strain B6C3F1 males, n = 4

lowfat1j

4.5% fat diet by weight, congenic C57BL/6J males, n = 3

21% fat diet by weight, congenic C57BL/6J males, n = 3

lowfat2k

very low fat diet, strain C57BL/6J males, n = 5

normal diet, strain C57BL/6J males, n = 5

agel

age 4 months, wild type, n = 4

age 32 months, wild type, n = 4

Each contrast corresponds to a test of differential expression and set of genes identified as significantly upregulated or downregulated. The A column lists treatments known or hypothesized to be associated with increased lifespan. The B column lists comparable control treatments associated with normal lifespan. All contrasts were of the form AB, such that upregulated genes exhibit increased expression in treatment A (relative to treatment B) while downregulated genes exhibit decreased expression in treatment A (relative to treatment B). The RNA source for all treatments was liver and the value n refers to the number of independent biological replicates available for each treatment.

aGSE3129, MG-U74Av2, Boyleston et al. [20]

bGSE3150, MG 430 2.0, Boyleston et al. [20]

cEMEXP153, MOE 430A, Amador-Noguez et al. [19]

dGSE988, MG-U74Av2, Rowland et al. [26]

eGSE5959, MG-U74Av2, Adamo et al. [27]

fEMEXP347, MOE 430A, Amador-Noguez et al. [87]

gGSE1093, MG-U74Av2, Tsuchiya et al. [24]

hGSE2431, MG-U74Av2, Dhahbi et al. [28]

iEMEXP490, MG-U74Av2, Heishi et al. [88]

jGSE363, MG-U74Av2, Recinos et al. [89]

kGSE3889, MG 430 2.0, Flowers et al. [90]

lEMEXP839, MG 430 2.0, Niedernhofer et al. [91]

The Snell, Ames and Little dwarf mutants were associated with large affects on gene expression at every age examined. Of 8525 probesets, the total number of genes differentially expressed with respect to these dwarf mutants ranged from 151 (contrast snell25) to as many as 987 (contrast ames3b) (see Table 2). IGF-I expression was significantly downregulated in all contrasts involving the Snell, Ames and Little dwarfs (Fig. 1, Table 2). IGF-I suppression, however, may account for only a small fraction of genes differentially expressed in these models. This is because GHR-KO mice exhibited significant IGF-I transcript downregulation of larger magnitude than that observed in the Snell, Ames and Little models, but nevertheless, only 46 genes exhibited differential expression in GHR-KO mice (see Table 2). This suggests that a relatively small number of genes may be specifically affected by IGF-I suppression in mouse liver. It is interesting to note, moreover, that IGF-I transcript levels are significantly downregulated with respect to the GHR-KI2 contrast, but not with respect to the GHR-KI1 contrast (Fig. 1, Table 2). This difference is most likely attributable to the STAT5 signaling pathway (signal transducer and activator of transcriptions 5) downstream of GHR. GHR-KI1 mice lack GHR residue 569, but exhibit only slight loss of hepatic IGF-1 transcript because STAT5 signaling is partially functional [26]. GHR-KI2 mice, however, lack GHR residue 391 and have no STAT5 signaling capacity, resulting in substantial loss of IGF-1 transcript in liver [26].
Table 2

Overview of Differential Expression Results

Contrast Symbol

IGF1 [log2(FC)]

Upregulated

Downregulated

snell5

-3.96*

129

187

snell25

-3.91*

65

86

ames4A

-3.20*

347

273

ames10A

-3.05*

314

214

ames22A

-1.55*

229

171

ames3B

-3.79*

567

420

ames6B

-3.18*

501

470

ames12B

-2.92*

424

307

ames24B

-1.03*

305

193

little3

-2.46*

359

291

little6

-2.54*

328

314

little12

-2.23*

200

310

little24

-1.99*

206

269

GHR-KO

-4.68*

8

38

GHR-KI1

-0.40

5

5

GHR-KI2

-5.29*

24

13

B6

-0.28

1

7

Gender

0.05

243

284

cr(2, 6)

-0.08

93

67

cr(2, 6)df

-0.85*

22

24

cr(20, 22)

0.26

65

116

cr(5, 22)

0.49

279

529

met

0.24

127

172

met (db/db)

0.20

166

523

glip

0.03

80

49

gm

0.32

114

107

ros

0.11

51

65

soy

0.17

45

32

lowfat1

0.20

6

8

lowfat2

-0.31

96

54

age

0.20

494

458

For each contrast, the number of probesets significantly upregulated and downregulated at a significance level of 0.05 is listed. The log-transformed fold-change [log2(FC)] associated with the insulin-like growth factor I (IGF1) transcript is also listed.

*test for differential expression significant with P < 0.05

http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_Fig1_HTML.jpg
Figure 1

IGF-I expression. The vertical axis corresponds to log-transformed gene expression values. The 31 contrasts listed in Table 1 are ordered along the horizontal axis. Red circles correspond to expression values of replicate observations associated with the A treatment of each contrast, while black diamonds indicate replicate expression values associated with the control B treatment of each contrast (see Table 1).

IGF-I expression was far more sensitive to the dwarf mutations than any of the dietary treatments examined, including CR (see Fig. 1). IGF-I was significantly downregulated in all contrasts involving dwarf mutants, but only one of four CR treatments resulted in significant IGF-I downregulation [contrast cr(2,6)df]. For all other dietary treatments, such as low fat diets and CR mimetics, IGF-I transcript levels were not significantly affected, and in fact, were slightly upregulated in most cases (Fig. 1, Table 2). Several genes exhibited expression patterns among all contrasts that were positively or negatively related to IGF-I expression (Fig. 2 and Additional files 1 and 2). These genes are potentially regulated by or coregulated with the IGF-I transcript in mouse liver tissue.
http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_Fig2_HTML.jpg
Figure 2

Potential IGF-I regulated or co-regulated genes. The vertical axis corresponds to log-transformed fold-change. The 31 contrasts listed in Table 1 are ordered along the horizontal axis. The thick black line represents fold-changes associated with IGF-I across the 31 contrasts. Red lines represent five genes for which fold-changes across contrasts are most positively associated with those of IGF-I (Mup3, Es31, Igfals, Keg1, Socs2). Green lines represent five genes for which fold-changes across contrasts are most negatively associated with those of IGF-I (Scd2, Slc16a7, Pcp4l1, Snm1, Igfbp1). Additional files 1 and 2 provide plots for the top 40 genes most positively and negatively associated with the IGF-I expression pattern among all contrasts.

Differential expression signatures

The differential expression signature associated with a given contrast was determined by testing H 0: μ Ai = μ Bi for all i = 1,..., 8525 probesets included in the analysis. Genes for which H 0: μ Ai = μ Bi was not rejected were assigned a score of 0. Genes for which H 0: μ Ai = μ Bi was rejected were assigned a score of ± 1, depending on whether the gene was upregulated (mean expression higher in treatment A relative to B) or downregulated (mean expression lower in treatment A relative to B). The pattern of 0, 1, and -1 scores among all i = 1, .., 8525 genes defined the differential expression signature of each contrast.

Three of the four dwarf mutants (Snell, Ames, Little) were associated with similar differential expression signatures (Fig. 3). While the GHR-KO mutation was associated with only a small number of transcriptional changes, most of these changes were also found in the other three long-lived models (Fig. 3). This similarity among long-lived models is reflected by a hierarchical cluster analysis of differential expression signatures (see Methods). All signatures associated with long-lived dwarf models clustered together in a single branch, and sub-branches within this single branch are joined at high levels of similarity (see Fig. 4). Dietary treatments, on the other hand, clustered together at lower levels of similarity and in some cases were placed in separate branches. It is noteworthy, for instance, that signatures associated with contrasts cr(2,6)df and cr(2,6) clustered in a branch apart from the signatures of cr(5,22) and cr(20,22). This reflects the degree to which the affect of CR on gene expression depends on how CR was carried out and, more particularly, on the age at which necropsies are performed.
http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_Fig3_HTML.jpg
Figure 3

Differential Expression Signatures. Each row corresponds to one of the contrasts listed in Table 1. Each column corresponds to one of 2192 genes differentially expressed with respect to more than one contrast. Rows have been ordered to correspond to the dendrogram shown in Figure 4. Red coloring indicates that a gene is upregulated (P < 0.05), while green coloring indicates that a gene is downregulated (P < 0.05).

http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_Fig4_HTML.jpg
Figure 4

Hierarchical cluster analysis. Each branch corresponds to a differential expression signature shown in Figure 3. The horizontal axis indicates the average level of similarity at which two clusters were joined (see Equation 1 in Methods).

There were many cases in which the similarity between signatures associated with two different contrasts was significantly greater than expected on the basis of chance (Fig. 5) (see Methods). This was generally true with regard to comparisons between any two dwarf mutation contrasts (e.g., snell5 vs. little3), but was less often the case with respect to comparisons between any two dietary treatments (e.g., ros vs. lowfat1). For every contrast, however, there was some evidence that its signature overlapped significantly with that of at least one other contrast (Fig. 5). The effects of gender overlapped significantly with those of the dwarf mutations as well as the four CR treatments. Additionally, differential expression signatures associated with lowfat diet treatments were significantly similar to those of dwarf mutations (especially the Ghrhr lit/lit mutation) and multiple CR treatments (see Fig. 5).
http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_Fig5_HTML.jpg
Figure 5

Differential expression signature similarity matrix. In the upper-right triangle region, dark red colors indicate high similarity between signatures associated with two contrasts (indicated by row and column labels). This similarity is defined by Equation (1) in the Methods section. The binary coding in the lower-left triangle region indicates whether signatures associated with two contrasts exhibit a significant level of similarity. The statistical procedure used to evaluate similarity is described in the Methods section. Contrast pairs with significant similarity (P < 0.05) are coded dark red, while pairs with non-significant similarity have no coloring.

The "age" contrast provides an indication of how young livers (4 months) differ from aged livers (32 months) at the gene expression level. Genetic alterations or dietary manipulations with signatures similar to that of the age contrast can be viewed as inducing a "reversal" of the aging process, at least insofar as gene expression levels in mouse liver are concerned. From Figure 3, it is clear that the age signature has fairly weak overlap with those of other contrasts. Overall, however, the age signature clusters together with dietary manipulations in a branch that includes CR and CR mimetic compounds (Fig. 3). This suggests that, in comparison to the dwarf mutations, the effects of CR and some CR mimetics more closely resemble a shift of the liver transcriptome towards a youthful state. Additionally, whereas the differential expression signature of the age contrast was significantly similar to those of some CR treatments, this wasn't the case for any contrasts involving dwarf mutants (see Fig. 5).

Longevity-associated genes I: long-lived dwarf mutants

A total of 13 genes were differentially expressed with respect to all four long-lived mouse models (Fig. 6). Most of these genes (10/13) were downregulated with respect to all four models, while only three genes were upregulated (Hao3, Sult2a2, Spink3). Six of the identified genes are localized to extracellular space, which is a significantly unlikely result based on gene ontology cellular component analysis (P < 0.01). Within the extracellular space, these six genes influence in a diversity of processes, including protein binding and transport (Mup4, Spink3), signaling pathways that regulate cell proliferation (Igf1, Lifr, Igfals) and detoxification (Es31). Significantly overrepresented biological process gene ontology terms included steroid metabolism, positive regulation of cell proliferation, positive regulation of cellular process, enzyme linked receptor protein signaling pathway, electron transport, and cellular morphogenesis (P < 0.05). The molecular function gene ontologies pheromone binding and oxiodreductase activity (acting on CH-OH group of donors) were also significantly overrepresented (P < 0.01).
http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_Fig6_HTML.jpg
Figure 6

Longevity-associated genes I. Listed genes are those that are differentially expressed with respect to each of four-long lived dwarf models (Snell, Ames, Little, GHR-KO). Each row corresponds to an individual candidate gene, while each column corresponds to one of the contrasts listed in Table 1. Red squares indicate significant upregulation, while green squares indicate significant downregulation.

Williams et al. [30] have recently reported on baseline expression levels of liver organs from 31 BxD mouse strains (GEO Series GSE6621). The lifespan of 21 of these strains had previously been measured by Gelman et al. [31]. Using these two data sources, the relationship between expression levels of candidate genes and mean lifespan across BxD mouse strains was examined (Fig. 7). This analysis provided some slight, though non-significant, additional support for a role of Hao3 in lifespan determination (Fig. 7a). Hao3 is a peroxisome-targeted hydroxyacid oxidase [32], and was significantly upregulated with respect to all four long-lived dwarf models. Correspondingly, long-lived BxD strains exhibited higher Hao3 expression (r s = 0.366). Taken alone, this correlation was marginally significant (P = 0.051), but was non-signficant following Benjamini-Hochberg adjustments for multiple testing (P = 0.322).
http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_Fig7_HTML.jpg
Figure 7

Candidate gene expression versus mean lifespan I. The expression level of nine candidate genes was examined among 21 BxD mouse strains. Expression data was generated by Williams et al. [30] (GEO series GSE6621). Lifespans of BxD strains were assayed by Gelman et al. [31]. The dashed horizontal line indicates the average gene expression level for each gene, while the solid line represents the least-squares regression estimate. Individual plots are shown for (A) Hao3, (B) Sult2a2, (C) Spink3, (D) Socs2, (E) Mup4, (F) Igfals, (G) Lifr, (H) Igf1 and (I) Efgr. The spearman rank correlation between expression and mean lifespan is shown in the upper right corner of each plot.

Longevity-associated genes II: dwarf mutants and CR

The same mechanisms underlying extended lifespan in dwarf mutants may contribute to the longevity extension that results from CR [2325]. Of genes differentially expressed in all four long-lived models, only three were also differentially expressed with respect to at least one CR contrast (Igf1, Igfals and Lifr) (Fig. 6). To the extent that longevity extension in dwarf mutants has a common basis with longevity extension by CR, these three genes are the most well supported longevity-associated genes identified by this analysis.

Previous criteria was altered to obtain an expanded list of gene candidates that are differentially expressed with respect to long-lived dwarf models and under CR. In particular, genes differentially expressed with respect to at least three of four long-lived models, as well as with respect to at least one of the four CR contrasts were identified, which yielded 30 additional candidate genes (Fig. 8). Significantly overrepresented biological process gene ontology terms among these genes included electron transport, secondary metabolism, amine metabolism, steroid biosynthesis, and establishment of localization (P < 0.05). The cellular components endoplasmic reticulum, microsome, membrane fraction and extracellular space were also significantly overrepresented (P < 0.01). A wide range of molecular function gene ontology terms were overrepresented, the most significant of which were monooxygenase activity, heme binding, iron ion binding and insulin-like growth factor binding (P < 0.01).
http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_Fig8_HTML.jpg
Figure 8

Longevity-associated genes II. Listed genes are those that are differentially expressed with respect to at least three of four long-lived dwarf models (Snell, Ames, Little, GHR-KO), and with respect to at least one of the four caloric restriction contrasts. Each row corresponds to an individual candidate gene, while each column corresponds to one of the contrasts listed in Table 1. Red squares indicate significant upregulation and green squares indicate significant downregulation.

The expression level of flavin-containing monoxygenase 3 (Fmo3) was positively correlated with mean lifespan among 21 BxD strains (r s = 0.410, see Fig. 9A). This correlation was significant when taken alone (P = 0.033), but not following p-value adjustment for multiple testing among all 33 genes considered (P = 0.366). In addition to Fmo3, five other genes listed in Figure 7 exhibited marginally significant correlations prior to multiple testing adjustment (Ero11b, Serpina12, Hes6, Cyp2f2 and Cyp4a14) (0.063 < P < 0.108; see Fig. 9B–9F). In each case, the relationship between baseline expression level and mean lifespan was consistent with that expected based on differential expression patterns in dwarf models and under CR (compare Fig. 8 to Fig. 9).
http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_Fig9_HTML.jpg
Figure 9

Candidate gene expression versus mean lifespan II. The expression level of six candidate genes was examined among 21 BxD mouse strains (see Fig. 6 caption). The dashed horizontal line indicates the average gene expression level for each gene, while the solid line indicates the estimates least-squares regression estimate. Individual plots are shown for (A) Fmo3, (B) Ero11b, (C) Serpina12, (D) Hes6, (E) Cyp2f2 and (F) Cyp4a14. The spearman rank correlation between expression and mean lifespan is shown in the upper right corner of each plot.

Discussion

Comparative analysis of multiple microarray datasets may provide insight not obtainable through analysis of any one dataset individually. This study presented a side-by-side comparison of microarray datasets generated from more than ten different studies of mouse liver organs, along with an analysis of each dataset using a consistent statistical methodology. The results provide a comprehensive view of how the effects of life-extending dwarf mutations, caloric restriction (CR), CR mimetic compounds, low fat diets, gender and aging are related at the level of gene expression. Comparisons among these various effects, based on expression patterns of more than 8500 genes, are of interest from multiple perspectives and shed light on a range of issues related to mammalian aging [e.g., see discussions in [16, 19, 23, 24, 3335]]. A total of 43 candidate longevity-associated genes were identified based upon common differential expression patterns among four long-lived dwarf models, or between dwarf models and CR treatments. Each gene was identified across multiple independent experiments, and is therefore very unlikely to be a false positive identification. The set of longevity-associated genes identified in this analysis therefore provides well-supported targets for future experimental investigation as well as potential clues for understanding the functional basis of lifespan extension in mammalian systems.

The GHR-KO mutation was associated with much smaller transcriptional effects than mutations carried by Snell, Ames and Little mice. This result is surprising, particularly since GHR-KO and Little mice both carry mutations that specifically affect GH signaling. Little mice are GH deficient because the pituitary is unresponsive to GHRH, while GHR-KO mice have elevated GH levels and lack GH receptor. Transcriptional changes associated with both mutants should therefore reflect inhibition of GH signaling and, on this basis, it might be expected that the lit and GHR-KO mutations have similar effects on gene expression patterns. This was not the case, however, given that 500–600 genes were differentially expressed with respect to Little mice, while fewer than 50 genes were differentially expressed with respect to GHR-KO mice. This difference does not reflect sample size and statistical power disparities between experiments (see Table 2), and moreover, agrees with an earlier study that found very few genes (10 of 2352) differentially expressed between GHR-KO mice and normal controls [23]. This observation is also consistent with phenotypic data, which indicate that Little and GHR-KO mice exhibit dissimilar magnitudes of lifespan extension [9, 15], as well as differing developmental weight-gain patterns [15]. It is possible that some disparity between differential expression signatures of GHR-KO and Little mice can be attributed to technical differences in sample preparation and array hybridizations. Another possibility is that the difference is due to levels of circulating growth hormone, which is nearly absent in Little mice and elevated in GHR-KO mice [9, 15]. This difference would influence gene expression patterns in liver if GH has systemic or local effects independent of the GH receptor. These considerations suggest that, in comparison to other dwarf mutants, GHR-KO mice may provide a more useful model for the purpose of elucidating the mechanisms of longevity extension in mammals. GHR-KO mice exhibit lifespan extension that is comparable to that associated with Snell, Ames and Little mice. However, since the downstream transcriptional effects of the GHR-KO mutation are much less extensive than those of other long-lived models, understanding how this mutation contributes to longevity extension may prove less difficult.

Dwarf mutations and CR both extend lifespan in mice, but the degree to which this effect is mediated by common mechanisms remains unclear [2325]. With regard to mouse liver tissue, this study reveals both similarities and differences between the effects of dwarf mutations and CR. The overall similarity between differential expression signatures associated with dwarf mutations and CR is significantly larger than expected on the basis of chance (Fig. 5). This is reflected by the 33 genes that exhibit parallel differential expression patterns in most long-lived mutants and under at least one CR treatment (Fig. 8). On the other hand, results of this analysis suggest that dwarf mutations have a much larger and more consistent effect on IGF-I expression than CR treatments (as well as potential CR mimetic compounds). This result was surprising, since several studies have shown that IGF-I expression and protein levels decline under CR [3638]. It was further surprising that of the four CR treatments, IGF-I expression was only decreased significantly when CR was applied to dwarf mice (which already exhibit IGF-I deficiency). Taken together, these results suggest that life-extending dwarf mutations have mechanisms in common with CR, but that IGF-I suppression may not necessarily be included among these shared mechanisms. A limitation related to this inference is that the effects of CR treatments on expression patterns vary considerably among different studies, depending on the duration of caloric restriction, age at which necropsies are performed and the laboratory in which CR is carried out (see Figs. 3 and 4). A robust generalization regarding how the transcriptional effects of dwarf mutations relate to those of CR would therefore require even more data from CR mice than was analyzed in the present study.

Candidate longevity-associated genes were identified based on one of either two premises. The first premise is that lifespan extension in Snell, Ames, Little and GHR-KO mice is due to shared GH/IGF-I alterations, such that genes exhibiting parallel transcriptional changes in all four models are likely to play an important role in longevity. The second premise is that longevity extension in dwarf models and under CR is achieved through common mechanisms, suggesting that genes with similar expression changes in multiple dwarf models and under CR are potentially important for longevity. There is evidence in favor of both of these premises, and both have been endorsed in previous analyses [17, 18, 20, 23]. It should be noted, however, that since the functional basis of mammalian lifespan is not well understood, the validity of each premise has not been established with certainty. Furthermore, all data examined in this study were generated from liver tissue. Liver tissue is the primary manufacturing site of IGF-I, but it is nonetheless possible that expression patterns in other tissues are of equal or greater importance in determining longevity. For instance, fibroblasts of long-lived mice are more stress-resistant than those of normal controls, suggesting that dwarf mutations affect a wide range of tissue and cell types [10, 11]. In fact, reduced IGF-I levels in circulation may be the main factor behind extended lifespan in dwarf mice, such that with respect to liver, the only expression change consequential for longevity determination is that of the IGF-I transcript. These considerations should be weighed when evaluating the potential role of candidate genes in determining lifespan and rates of aging in mice and other mammalian species.

Only three genes were differentially expressed with respect to all four long-lived dwarf models with corresponding effects under at least one CR treatment (IGF-I, Igfals and Lifr). Both Igf1 and Igfals have a well-documented role in lifespan determination [39, 40], but a potential role of leukemia inhibitory factor receptor (Lifr) in longevity has not been widely explored. Lifr expression was downregulated at all ages in four long-lived dwarf models, by short-term CR and by low-fat diet. Determining whether Lifr downregulation in liver contributes to longevity extension represents a difficult task, since Lifr is highly pleiotropic and has been associated with a wide range of biological effects. While there are several molecules that Lifr may interact with [4143], Lifr expression is a primary determinant of cellular responsiveness to leukemia inhibitory factor (Lif) [44]. The Lif glycoprotein is a member of the IL-6 type cytokine family and exhibits a wide range of effects among different cell types. At a systemic level, excess Lif has been associated with low body weight, hypermotility, overgrowth of bone, calcification of several organs, loss of spermatocytes in males and severe loss of adipose tissue [44], and in the liver specifically, Lif has been found to stimulate triglyceride secretion [45]. Interestingly, however, Lif(-/-) mutants have been found to develop normally and possess generally good health, although pregnancy does not occur in females due to problems with blastocyst implantation [46]. To some degree, inhibition of Lif/Lifr signaling may be countered in long-lived mice by down-regulation of Socs2 (suppressor of cytokine signaling) [47, 48], which was also down-regulated in all four long-lived dwarf models. Declines in Lifr expression may increase sensitivity to drug-induced liver disease [49], which in general is contrary to the stress resistance characteristics of long-lived mice [11], but is consistent with one study that found decreased resistance to acetaminophen toxicity in multiple dwarf strains [50].

Fmo3 and Cyp2f2 were differentially expressed with respect to most (3/4) long-lived dwarf models and also under the majority (3/4) of CR treatments that were examined. In addition, the expression level of both genes varied with mean lifespan among 21 BxD recombinant inbred mouse strains in a direction consistent with differential expression analyses. Since both genes are monooxygenases with established roles in drug metabolism, their expression in liver could mediate the life-extending effects of dwarf-mutations and CR through similar mechanisms. Fmo3 is a flavin containing monoxygenase for which expression levels are elevated in long-lived dwarf mutants, female mice, CR mice and mice treated with each of three different CR mimetic compounds (metaformin, glipizide and rosiglitazone). The Fmo3 gene has received considerable attention in studies of human populations, since it exhibits extensive polymorphism among individuals and may be a factor promoting differential drug response [51]. With regard to aging processes, a plausible hypothesis is that elevated Fmo3 levels in liver contributes to extended longevity by increasing stores of glutathione (GSH) and glutathione-S-transferase (GST), which then leads to enhanced resistance to oxidative stress. This hypothesis is appealing, since Fmo3 has an important role in sulfoxidation of methionine [52, 53], and it has been found that long-lived Ames mice exhibit elevated methionine metabolism [54], which may underlie increased GSH, GST and oxidative stress resistance in long-lived models [54].

Cyp2f2 is a cytochrome P450 enzyme, and in contrast to Fmo3, its expression was downregulated in long-lived and CR mice. Several other cytochrome P450 enzymes exhibited differential expression, with slightly less consistency, among multiple long-lived models and under CR treatments (Cyp4a14, Cyp2c38, Cyp2c39 and Cyp4a12b). Although Cyp2f2 displayed the most consistent differential expression pattern, it is unlikely that Cyp2f2 would have a unique role in mammalian aging apart from other cytochrome P450 enzymes. Cytochrome P450 enzymes are involved in xenobiotic detoxification and are catalysts for a large number of metabolic reactions. It is unclear how this system may influence aging in mammals, but multiple P450 enzymes are also differentially expressed in long-lived C. elegans daf-2 mutants, which has led to the suggestion that P450 enzymes slow aging by reducing damage generated from toxic compounds [55, 56]. From the standpoint of mammalian aging, it is of interest that some P450 enzymes, including Cyp2f2, are regulated by peroxisome proliferator-activated receptors (PPARs) [57, 58]. Cytochrome P450 enzymes may therefore be part of a broader cellular response that has previously been associated with longevity in long-lived and CR mice [59]. An interesting pattern is that, in many cases, P450 enzyme expression changes in long-lived mice and under CR were mirrored by expression differences in females relative to males. In a comprehensive analysis of 41 cytochrome P450 enzyme genes, it was found that 14 genes exhibited this pattern to some degree (data not shown). This result may be attributable to the role of GH as a factor regulating cytochrome P450 expression [60].

The correlation structure of expression patterns among genes may be as informative as differential expression analyses but is often an underexploited aspect of microarray datasets. Similarity among expression patterns across many conditions can be used, for example, to make inferences and generate testable hypotheses regarding interactions between genes [61, 62]. Along these lines, genes with expression patterns among contrasts that closely corresponded with those of the IGF-I transcript were identified in this study. This analysis identified Snm1 (DNA cross-link repair 1a; also called Dclre1a) as a potential IGF-I regulated or co-regulated gene (see Fig. 2 and additional file 2). The induction pattern of Snm1 is opposite that of IGF-I across contrasts, such that its expression is increased in proportion to IGF-I inhibition. Snm1(-/-) mice have reduced lifespan, which is primarily due to elevated mortality from bacterial infection and cancer, suggesting that Snm1 is a tumor suppressor with an immunological role [63]. At the cellular level, overexpression of Snm1 in yeast increases resistance to genotoxic stress agents that induce DNA cross-links [64], and stem cells derived from Snm1(-/-) mice are sensitive to the cross-linking agent mitomycin C [65]. These results suggest that elevated Snm1 expression may be a causal factor underlying resistance to the DNA-alkylating agent methyl methanesulfonate that has previously been found in cell lines derived from long-lived mice [10, 11]. It would thus be worthwhile to evaluate whether Snm1 expression is elevated in such cell lines (as in hepatic tissue), and if so, whether these cells are resistant to DNA cross-link inducing agents (e.g., nitrogen mustard, cisplatin). Interestingly, recent studies report that DNA repair mutations leading to accelerated senescence phenotypes are, like long-lived dwarf mutations, associated with IGF-I inhibition [66]. It is therefore surprising that in long-lived and CR mice, IGF-I inhibition is associated with elevated Snm1 expression, which may enhance DNA repair and promote genomic stability.

The present study has demonstrated analytical approaches for the comparative analysis of microarray datasets that may have application in other contexts. Microarray analyses based upon individual datasets often identify an exceptionally large number of genes, which limits the utility of microarray data as a tool for selecting candidates in follow-up studies. Recently, however, public repositories of high-quality microarray data have been established [67, 68], and statistical methods aimed at comparative analysis of these resources continue to be developed [6972]. Increasingly, therefore, comparative analysis provides a feasible approach for filtering out false-positive identifications and identifying transcripts most consistently supported across multiple experiments. This yields a set of candidate genes that is necessarily smaller and more tractable for subsequent experimental investigation, and moreover, each identified gene is more likely to represent a statistically significant finding. Furthermore, apart from the identification of individual gene candidates, comparison of expression datasets provides unique insight into genome-wide patterns among studies.

The maximum human lifespan is approximately twice as large as the maximum chimpanzee lifespan [73], which is more than ten times larger than the maximum mouse lifespan [2]. The evolutionary lineage connecting humans to chimpanzees to mouse may therefore represent a remarkable instance of lifespan extension and delayed rates of aging. Inhibition of IGF-I signaling is (at present) the only known genetic manipulation that extends lifespan in multiple species, so it is tempting to ask whether this pathway has contributed to evolutionary extension of lifespan within the mammalian lineage. It is interesting to note that, contrary to longevity extension via nearly all IGF-I signaling mutations, lifespan extension in the mouse-chimp-human lineage has been accompanied by an increase in body size. This observation is encouraging with regard to the possibility of developing interventions that delay the onset of age-related disease in humans without undesirable consequences (e.g., dwarfism). Clarke et al. [74] used sequence data and d N /d S ratios to identify mouse-chimp-human orthologs that exhibit more rapid evolutionary change than expected based upon neutral substitution models (i.e., positive selection). Supplemental data from their study provides no indication of positive selection with respect to IGF-I, providing little indication that IGF-I sequence changes have been consequential during human evolution. It is interesting to note, however, that Dio1 (iodothyronine deiodinase), which is downregulated in most dwarf mouse strains and under CR (see Fig. 8), has been associated with positive selection and accelerated evolution within the mouse-chimp-human lineage (see supplemental data from [74]). This is also the case for Papp-A (pregnancy associated plasma protein A) [74], which was not identified in the present study, but is of importance since Papp-A-KO mice exhibit diminished IGF-I bioactivity and extended lifespan [75]. These results are suggestive, but given the many phenotypic differences between mice, chimpanzees and humans (besides lifespan), this evolutionary criterion does not establish the functional significance of Dio1, Papp-A or other positively selected genes.

These findings provide a useful reference point for future experimental studies of long-lived dwarf mice and mammalian aging. It would be of interest, for example, to determine whether genes identified in the present study are differentially expressed with respect to other long-lived mouse strains for which data is not currently available. Similar gene expression changes may be found in liver tissue of p66shc knockouts [76], IGF-IR knockouts [77], Klotho transgenic mice [78] and Papp-A-KO mice [75], since each of these long-lived models carry mutations that also inhibit the IGF-I signaling pathway. Because these mutations inhibit IGF-I signals at points further downstream than those considered in this study, expression data generated from these models could be combined with results of the present study to assemble a more comprehensive picture of hepatic IGF-I signaling. Since the currently known life-extending mutations are maintained on differing genetic backgrounds [1], evaluating potential background effects on lifespan and gene expression patterns will also be critical for elucidating hepatic IGF-1 signaling pathways [79]. Ultimately, however, it is important to develop a systemic model of GH/IGF-I signaling, with further evaluation of the respective roles of IGF-I and GH in extended longevity. IGF-I is manufactured in a wide variety of tissue types, and in some cases IGF-I production may be independent of GH signaling [80]. It is therefore important to consider not only local effects of GH reduction on gene expression in liver, but also effects of GH/IGF-I in other tissue types. Studies of long-lived C. elegans mutants, for instance, have suggested that IGF-I signaling in nervous tissue may be consequential in determining lifespan [81]. It would therefore be useful to evaluate whether genes identified in the present study are also differentially expressed in non-hepatic tissues of Ames, Snell, Little or GHR-KO mice.

Conclusion

Several mutations are known to increase longevity in mouse, but most of these mutations are associated with pleiotropic effects, which include traits that are undesirable from a therapeutic perspective (e.g., dwarfism). A key challenge in future work towards potential longevity-promoting therapeutic compounds is to decouple the positive life-extending effects of IGF-I inhibition from diminutions in growth and body size. Determining whether this is a possibility requires a more complete functional understanding of longevity extension in long-lived mutant mouse models. The present study has taken steps in this direction by presenting comparisons between the transcriptional effects of dwarf mutations and those of dietary treatments, gender and aging. Among other findings, these comparisons reveal that lit/lit and GHR-KO mutations exhibit transcriptional effects of surprisingly different magnitude, transcriptional effects of CR (and some CR mimetics) resemble aging more closely then those of dwarf mutations, and that the transcriptional effects of gender and lowfat diets overlap significantly with those of dwarf mutations and CR. A total of 43 genes with the highest levels of support as longevity-associated transcripts were identified in this context. The main value of identified genes is the potential that, as putative downstream elements in the IGF-I signaling pathway, some genes may be more directly (and functionally) involved in longevity determination than upstream IGF-I signaling components. Experimental studies focusing on identified genes may therefore enhance our functional understanding of how mutations affecting IGF-I signaling lead to extended lifespan and deceleration of aging in mammals.

Methods

Gene expression datasets analyzed in this study were selected from those available in the Gene Expression Omnibus (GEO) and ArrayExpress depositories for MIAME-compliant microarray data [67, 68]. Expression datasets were generated using one of three Affymetrix microarray platforms (MG-U74A, MOE430A, 430 2.0). Expression values were calculated by either the MAS 5.0 algorithm or Robust Multichip Average [82]. When raw CEL files were available, MAS 5.0 generated datasets were re-normalized using RMA.

Differential expression testing was performed using the Limma linear modeling package available in the R Bioconductor software suite [83]. This approach fits a linear model to expression values associated with each individual gene, and the distribution of sample residual variances among all genes is used to stabilize the residual variance estimates of individual genes by shrinkage towards a prior value. This limits false-positive gene identifications arising from underestimated residual variances for datasets with low levels of replication in each experimental treatment. Separate linear model analyses were conducted for contrast groups sharing the same superscript in Table 1. Expression values were not combined, therefore, when expression was assayed using different platforms or when experiments were performed in different laboratories. For each contrast, P-values were adjusted across genes using the Benjamini-Hochberg method [84]. When multiple contrasts were specified for a single dataset, the nested-F test approach was used to evaluate the significance of moderated t-statistics associated with each individual contrast [83].

Longevity-associated genes were identified based on common differential expression patterns among contrasts using criteria described in the Results section. The GOstats R Bioconductor package was used evaluate whether particular gene ontology terms were overrepresented with respect to identified genes [85]. The gene universe in each test was defined as the 8525 probesets included in the analysis, i.e., those that could be matched among the MG-U74A, MOE430A and 430 2.0 Affymetrix platforms [85].

Differential expression signature similarity metric

A similarity measure was developed to compare gene expression signatures associated with different contrasts. This similarity is intuitive in terms of Venn diagrams, where the similarity between two signatures is proportional to the overlap between sets of genes differentially expressed with respect to each contrast individually. This notion was generalized in order to define a measure of similarity that can be used for clustering signatures associated with different contrasts. Consider two signatures that have been generated by two contrasts α and β. Using the notations defined in Table 3, the signatures contain n +,+ + n -,- genes with identical differential expression patterns, n +,- + n -,+ genes with opposite differential expression patterns, and n +,0 + n 0,+ + n -,0 + n 0,- genes differentially expressed with respect to just one of the two contrasts. Given these values, Equation (1) defines a measure of similarity (s) between two different signatures, where 0 ≤ s ≤ 1.
Table 3

Notations associated with Equation (1)

  

Contrast α

  

Upregulated

Downregulated

H0 not rejected

Contrast β

Upregulated

n +,+

n -,+

n 0,+

 

Downregulated

n +,-

n -,-

n 0,-

 

H0 not rejected

n +,0

n -,0

n 0,0

The null hypothesis H 0: μ Ai = μ Bi was evaluated for N genes with respect to contrasts α and β (see text). With respect to each contrast individually, genes are either upregulated, downregulated, or H 0 is not rejected. With respect to both contrasts α and β, therefore, all N genes are classified into one of nine categories. The number of genes assigned to each category is indicated by the values of n given in the table.

http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_Equ1_HTML.gif

There are several plausible similarity measures between two different signatures. A main advantage of the similarity measure defined by Equation (1) is that similarity not increased by the n 0,0 genes that are not differentially expressed with respect to either contrast. Thus, emphasis is placed only on the minority of genes that are differentially expressed. Additionally, the denominator of Equation (1) is configured such that two signatures may be similar even if they differ greatly in the total number of genes that are differentially expressed, provided that the contrast associated with less differential expression yields approximately n +,+ + n -,- differentially expressed genes. This is sensible in consideration of the fact that the number of differentially expressed genes generated by a given contrast depends on the sample sizes used in the experimental treatments being compared. For signature comparisons in the present study, the value n +,- + n -,+ was negligible and therefore not included in Equation (1). If n +,- + n -,+ had been large, however, it would have been appropriate to add this term to the denominator of Equation (1) when evaluating similarity between signatures.

Differential expression signatures for each contrast listed in Table 1 were compared using the similarity measure defined by Equation (1). The distance between two signatures was defined as 1 – s and a hierarchical cluster analysis of signatures associated with the 31 contrasts specified in Table 1 was performed (Fig. 4). Groups of signatures were joined using the average distance method. This yielded a dendrogram providing an indication of which contrasts were associated with similar differential expression patterns.

Statistical evaluation of overlap between differential expression signatures

A statistical procedure was developed to evaluate whether differential expression signatures associated with two different contrasts exhibited a significant level of similarity. Among the 31 contrasts, the correspondence between all 31(31-1)/2 = 465 pairwise combinations of signatures was evaluated. The procedure used is similar to that proposed by Smid et al. [86]. The null hypothesis for this statistical test is that, for any two contrasts α and β, the probability that a gene is differentially expressed with respect to contrast α is independent of the probability that a gene is differentially expressed with respect to contrast β. It should be noted that the three contrasts with a d superscript in Table 1 were defined using a common control treatment as a reference. It is thus expected a priori that differential expression signatures of these three contrasts will exhibit association beyond that stated by the null hypothesis defined above. This is also true of the seven contrasts carrying an h superscript in Table 1.

The value T defined below is proportional to s and serves as the test statistic.T = n +,+ + n -,-

A total of N genes are considered in the analysis. It is given that http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_IEq1_HTML.gif genes upregulated with respect to contrast α, http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_IEq2_HTML.gif genes downregulated with respect to contrast α, http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_IEq3_HTML.gif genes upregulated with respect to contrast β, and http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_IEq3_HTML.gif genes downregulated with respect to contrast β. Given these quantities, http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_IEq4_HTML.gif is the probability of upregulation with respect to contrast α, http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_IEq5_HTML.gif is the probability of downregulation with respect to contrast α, http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_IEq6_HTML.gif is the probability of upregulation with respect to contrast β, and http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_IEq7_HTML.gif is the probability of downregulation with respect to contrast β. Under the null hypothesis stated above, http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_IEq4_HTML.gif is independent of http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_IEq6_HTML.gif and http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_IEq5_HTML.gif is independent of http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_IEq7_HTML.gif . Consequently, the probability that a gene is jointly upregulated with respect to contrasts α and β is given by p 1, while the probability that a gene is jointly downregulated with respect to contrasts α and β is given by p 2.
http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_Equ2_HTML.gif
http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_Equ3_HTML.gif
Random variable X denotes the number of jointly upregulated genes among all N genes, and random variable Y equals the number of jointly downregulated genes among all N genes, such that T = X + Y. The distributions of X and Y are binomial under the null hypothesis.
http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_Equ4_HTML.gif
http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-8-353/MediaObjects/12864_2007_Article_1066_Equ5_HTML.gif

When Np 1 and Np 2 are large, P(X = x) and P(Y = y) and their convolution T = X + Y are approximately Normal under H 0. When Np 1 and Np 2 are small, P(X = x) and P(Y = y) and their convolution T = X + Y are approximately Poisson under H 0. The significance of observed T values for each pair of contrasts by comparison to either Normal or Poisson cumulative distribution functions, depending on the magnitude of Np 1 and Np 2. This test was carried out for all pairwise combinations of the 31 contrasts evaluated, yielding a total of 31(31-1)/2 = 465 p-values. All p-values were confirmed by simulation analyses. To control for multiple testing, p-values were adjusted using the Benjamini-Hochberg method.

Declarations

Acknowledgements

This work was supported by NIA training grant AG000114 and the University of Michigan Department of Pathology. Dr. Richard A. Miller and three anonymous reviewers provided helpful comments and suggestions on this manuscript. The contributors of primary microarray data sources are thanked for providing their data to public databases.

Authors’ Affiliations

(1)
Department of Pathology, University of Michigan

References

  1. Liang H, Masoro EJ, Nelson JF, Strong R, McMahan CA, Richardson A: Genetic mouse models of extended lifespan. Experimental Gerontology 2003, 38:1353–1364.View ArticlePubMed
  2. Miller RA: Genetic approaches to the study of aging. J Amer Ger Soc 2005, 53:S284-S286.View Article
  3. Kinney BA, Meliska CJ, Steger RW, Bartke A: Evidence that Ames dwarf mice age differently from their normal siblings in behavioral and learning and memory parameters. Horm Behav 2001, 39:277–284.View ArticlePubMed
  4. Silberberg R: Articular aging and osteroarthritis in dwarf mice. Pathol Microbiol (Basel) 1972, 38:417–430.
  5. Ikeno Y, Bronson RT, Hubbard GB, Lee S, Bartke A: Delayed occurrence of fatal neoplastic diseases in Ames dwarf mice: correlation to extended longevity. J Gerontology 2003, 58A:291–296.
  6. Vergara M, Smith-Wheelock M, Harper JM, Sigler R, Miller RA: Hormone-treated snell dwarf mice regain fertility but remain long-lived and disease resistance. J Gerontology 2004, 59A:1244–1250.
  7. Bielschowsky F, Bielschowsky M: Carcinogenesis in the pituitary of dwarf mouse. The response to dimethylbenzanthracene applied to the skin. Br J Cancer 1961, 15:257–262.View Article
  8. Rennels EG, Anigstein DM, Anigstein L: A cumulative study of the growth of sarcoma 180 in anterior pituitary dwarf mice. Tex Rep Biol Med 1965, 23:776–781.PubMed
  9. Flurkey K, Papaconstantinou J, Miller RA, Harrison DE: Lifespan extension and delayed immune and collagen aging in mutant mice with defects in growth hormone production. Proc Natl Acad Sci USA 2001, 98:6736–6741.View ArticlePubMed
  10. Murakami S, Salmon A, Miller RA: Multiplex stress resistance in cells from long-lived dwarf mice. FASEB J 2003, 17:1565–1566.PubMed
  11. Salmon AB, Murakami S, Bartke A, Kopchick J, Yasumura K, Miller RA: Fibroblast cell lines from young adult mice of long-lived mutant strains are resistant to multiple forms of stress. Am J Physiol Endocrinol Metab 2005, 289:E23-E29.View ArticlePubMed
  12. Brown-Borg HM, Borg KE, Meliska CJ, Bartke A: Dwarf mice and the aging process. Nature 1996, 384:33.View ArticlePubMed
  13. Flurkey K, Papaconstantinou J, Harrison DE: The snell dwarf mutation Pit1(dw) can increase lifespan in mice. Mech Ageing Dev 2002, 123:121–130.View ArticlePubMed
  14. Godfrey P, Rahal JO, Beamer WG, Copeland NG, Jenkins NA, Mayo KE: Ghrh receptor of little mice contains a missense mutation in the extracellular domain that disrupts receptor function. Nat Genet 1993, 4:227–232.View ArticlePubMed
  15. Coschigano KT, Clemmons D, Bellush LL, Kopchick JJ: Assessment of growth parameters and lifespan of GHR/BP gene-disrupted mice. Endocrinology 2000, 141:2608–2613.View ArticlePubMed
  16. Carter CS, Ramsey MM, Sonntag WE: A critical analysis of the role of growth hormone and IGF-I in aging and lifespan. Trends in Genetics 2002, 18:295–301.View ArticlePubMed
  17. Dozmorov IA, Bartke A, Miller RA: Array-based expression analysis of mouse liver genes: effect of age and the longevity mutant Prop1 df . Journal of Gerontology 2001, 56A:B72-B80.
  18. Dozmorov IA, Galecki A, Chang Y, Krzesicki R, Vergara M, Miller RA: Gene expression profile of long-lived Snell dwarf mice. Journal of Gerontology 2002, 57A:B99-B108.
  19. Amador-Noguez D, Yagi K, Venable S, Darlington G: Gene expression profile of long-lived Ames dwarf mice and Little mice. Aging Cell 2004, 3:423–441.View ArticlePubMed
  20. Boylston WH, DeFord JH, Popaconstantinou J: Identification of longevity-associated genes in long-lived Snell and Ames dwarf mice. Age 2006, 28:125–144.View ArticlePubMed
  21. Miller RA, Galecki A, Shmookler-Reis RJ: Interpretation, design, and analysis of gene array expression experiments. J Gerontology 2001, 56A:B52-B57.
  22. Fisher RA: Combining independent tests of significance. American Statistician 1948, 2:30.View Article
  23. Miller RA, Chang Y, Galecki AT, Al-Regaiey K, Kopchick JJ, Bartke A: Gene expression patterns in calorically restricted mice: partial overlap with long-lived mutant mice. Molecular Endocrinology 2002, 16:2657–2666.View ArticlePubMed
  24. Tsuchiya T, Dhahbi JM, Cui X, Mote PL, Bartke A, Spindler SR: Additive regulation of hepatic gene expression by dwarfism and caloric restriction. Physiol Genomics 2004, 17:307–315.View ArticlePubMed
  25. Bonkowski MS, Rocha JS, Masternak MM, Al Regaiey KA, Bartke A: Targeted disruption of growth hormone receptor interferes with the beneficial actions of calorie restriction. Proc Natl Acad Sci 2006, 103:7901–7905.View ArticlePubMed
  26. Rowland JE, Lichanska AM, Kerr LM, White M, d'Aniello EM, Maher SL, Brown R, Teasdale RD, Noakes PG, Waters MJ: In vivo analysis of growth hormone receptor signaling domains and their associated transcripts. Molecular and Cellular Biology 2005, 25:66–77.View ArticlePubMed
  27. Adamo ML, Ma X, Ackert-Bicknell CL, Donahue LR, Beamer WG, Rosen CJ: Genetic increase in serum insulin-like growth factor-I (IGF-I) in C3H/HeJ compared with C57BL/6J mice is associated with increased transcription from the IGF-I exon 2 promoter. Endocrinology 2006, 147:2944–2955.View ArticlePubMed
  28. Dhahbi JM, Mote PL, Fahy GM, Spindler SR: Identification of potential caloric restriction mimetics by microarray profiling. Physiol Genomics 2005, 23:343–350.View ArticlePubMed
  29. Anisimov VN, Berstein LM, Egormin PA, Piskunova TS, Popovich IG, Zabezhinski MA, Kovalenko IG, Poroshina TE, Semenchenko AV, Provinciali M: Effect of metformin on life span and on the development of spontaneous mammary tumors in HER-2/neu transgenic mice. Exp Gerontol 2005, 40:685–693.View ArticlePubMed
  30. Williams RBH, Cotsapas CJ, Cowley MJ, Chan E, Nott DJ, Little PFR: Normalization procedures and detection of linkage signal in genetical-genomics experiments. Nat Genetics 2006, 38:855–856.View Article
  31. Gelman R, Watson A, Bronson R, Yunis E: Murine chromosomal regions correlated with longevity. Genetics 1998, 118:693–704.
  32. Jones JM, Morrell JC, Gould SJ: Identification and characterization of HAOX1, HAOX2 and HAOX3, three human peroxisomal 2-hydroxy acid oxidases. J Biol Chem 2000, 275:12590–12597.View ArticlePubMed
  33. Park S-K, Prolla TA: Lessons learned from gene expression profile studies of aging and caloric restriction. Ageing Research Reviews 2005, 4:55–65.View ArticlePubMed
  34. Ingram DK, Zhu M, Mamczarz J, Zou S, Lane MA, Roth GS, deCabo R: Calorie restriction mimetics: an emerging research field. Aging Cell 2006, 5:97–108.View ArticlePubMed
  35. Flurkey C, Currer JM, Harrison DE: The mouse in Aging Research. The Mouse in Biomedical Research 2 Edition (Edited by: Fox JG, Barthold S, Davisson M, Newcomer CE, Quimby FW, Smith A). Massachusetts: Elsevier 2007, 3:637–672.View Article
  36. Maiter D, Fliesen T, Underwood LE, Maes M, Gerard G, Davenport ML, Ketelslegers JM: Dietary protein restriction decreases insulin-like growth factor I independent of insulin and liver growth hormone binding. Endocrinology 1989, 124:2604–2611.View ArticlePubMed
  37. Oster MH, Fielder PJ, Levin N, Cronin MJ: Adaptation of the growth hormone and insulin-like growth factor-1 axis to chronic and severe calorie or protein malnutrition. J Clin Invest 1995, 95:2258–2265.View ArticlePubMed
  38. Goya L, García-Segura LM, Ramos S, Pascual-Leone AM, Argente J, Martin MA, Chowen JA: Interaction between malnutrition and ovarian hormones on the systemic IGF-I axis. Eur J Endocrin 2002, 147:417–424.View Article
  39. Hsieh CC, DeFord JH, Flurkey K, Harrison DE, Papaconstantinou J: Implications of the insulin signaling pathway in Snell dwarf mouse longevity: a similarity with the C. elegans longevity paradigm. Mech Ageing Dev 2002, 123:1229–1243.PubMed
  40. Bartke A: Minireview: role of growth hormone/insulin-like growth factor system in mammalian aging. Endocrinology 2005, 146:3718–3723.View ArticlePubMed
  41. Di Marco A, Gloaguen I, Demartis A, Saggio I, Graziani R, Paonessa G, Laufer R: Agonistic and antagonistic variants of ciliary neurotrophic factor (CNTF) reveal functional differences between membrane-bound and soluble CNTF alpha-receptor. J Biol Chem 1997, 272:23069–23075.View ArticlePubMed
  42. Robledo O, Fourcin M, Chevalier S, Guillet C, Auguste P, Pouplard-Barthelaix A, Pennica D, Gascan H: Signaling of the cardiotrophin-1 receptor. Evidence for a third receptor component. J Biol Chem 1997, 272:4855–4863.View ArticlePubMed
  43. Nakashima K, Taga T: gp130 and the IL-6 family of cytokines: signaling mechanisms and thrombopoietic activities. Semin Hematol 1998, 35:210–221.PubMed
  44. Metcalf D: The unsolved enigmas of leukemia inhibitory factor. Stem Cells 2003, 21:5–14.View ArticlePubMed
  45. Nonogaki K, Pan XM, Moser AH, Shigenaga J, Staprans I, Sakamoto N, Grunfeld C, Feingold KR: LIF and CNTF, which share the gp130 transduction system, stimulate hepatic lipid metabolism in rats. Am J Physiol 1996, 271:E521-E528.PubMed
  46. Stewart CL, Kaspar P, Brunet LJ, Bhatt H, Gadi I, Kontgen F, Abbondanzo SJ: Blastocyst implantation depends on maternal expression of leukemia inhibitory factor. Nature 1992, 359:76–79.View ArticlePubMed
  47. Starr R, Willson TA, Viney EM, Murray LJ, Rayner JR, Jenkins BJ, Gonda TJ, Alexander WS, Metcalf D, Nicola NA, Hilton DJ: A family of cytokine-inducible inhibitors of signalling. Nature 1997, 387:917–921.View ArticlePubMed
  48. Nicholson SE, Willson TA, Farley A, Starr R, Zhang JG, Baca M, Alexander WS, Metcalf D, Hilton DJ, Nicola NA: Mutational analyses of the SOCS proteins suggest a dual domain requirement but distinct mechanisms for inhibition of LIF and IL-6 signal transduction. EMBO J 1999, 18:375–385.View ArticlePubMed
  49. Masubuchi Y, Bourdi M, Reilly TP, Graf ML, George JW, Pohl LR: Role of interleukin-6 in hepatic heat shock protein expression and protection against acetaminophen-induced liver disease. Biochem Biophys Res Commun 2003, 304:207–212.View ArticlePubMed
  50. Harper JM, Salmon AB, Chang Y, Bonkowski M, Bartke A, Miller RA: Stress resistance and aging: influence of genes and nutrition. Mech Age Dev 2006, 127:687–694.View Article
  51. Zhou J, Shephard EA: Mutation, polymorphism and perspectives for the future of human flavin-containing monooxygenase 3. Mutat Res 2006, 612:165–171.View ArticlePubMed
  52. Duescher RJ, Lawton MP, Philpot RM, Elfarra AA: Flavin-containing monooxygenase (FMO)-dependent metabolism of methionine and evidence for FMO3 being the major FMO involved in methionine sulfoxidation in rabbit liver and kidney microsomes. J Biol Chem 1994, 269:17525–17530.PubMed
  53. Krause RJ, Ripp SL, Sausen PJ, Overby LH, Philpot RM, Elfarra AA: Characterization of the methionine S-oxidase activity of rat liver and kidney microsomes: immunochemical and kinetic evidence for FMO3 being the major catalyst. Arch Biochem Biophys 1996, 333:109–116.View ArticlePubMed
  54. Uthus EO, Brown-Borg HM: Methionine flux to transsulfuration is enhanced in the long living Ames dwarf mouse. Mech Age Dev 2006, 127:444–450.
  55. McElwee JJ, Schuster E, Blanc E, Thomas JH, Gems D: Shared transcriptional signature in Caenorhabditis elegans dauer larvae and long-lived daf-2 mutants implicates detoxification system in longevity assurance. J Biol Chem 2004, 279:44533–44543.View ArticlePubMed
  56. Gems D, McElwee JJ: Broad spectrum detoxification: the major longevity assurance process regulated by insulin/IGF-I signaling. Mech Age Dev 2005, 126:381–387.View Article
  57. Ye X, Lu L, Gill SS: Suppression of cytochrome P450 Cyp2f2 mRNA levels in mice by peroxisome proliferator diethylhexylphthalate. Biochem Biophys Res Commun 1997, 239:660–665.View ArticlePubMed
  58. Ito O, Nakamura Y, Tan L, Ishizuka T, Sasaki Y, Minami N, Kanazawa M, Ito S, Sasano H, Kohzuki M: Expression of cytochrome P-450 enzymes in the kidney and liver: regulation by PPAR and species-difference between rat and human. Mol Cell Biochem 2006, 284:141–148.View ArticlePubMed
  59. Corton JC, Apte U, Anderson SP, Limaye P, Yoon L, Latendresse J, Dunn C, Everitt JI, Voss KA, Swanson C, Kimbrough C, Wong JS, Gill SS, Chandraratna RAS, Kwak M-K, Kensler TW, Stulnig TM, Steffensen KR, Gustafsson JÅ, Mehendale HM: Mimetics of caloric restriction include agonists of lipid-activated nuclear receptors. J Biol Chem 2004, 279:46204–46212.View ArticlePubMed
  60. Waxman DJ, O'Connor C: Growth hormone regulation of sex-dependent liver gene expression. Mol Endocrinol 2006, 20:2613–2629.View ArticlePubMed
  61. Oldham MC, Horvath S, Geschwind DH: Conservation and evolution of gene coexpression networks in human and chimpanzee brains. Proc Natl Acad Sci 2006, 103:17973–17978.View ArticlePubMed
  62. Yip AM, Horvath S: Gene network interconnectedness and the generalized topological overlap measure. BMC Bioinform 2007, 8:22.View Article
  63. Ahkter S, Richie CT, Zhang N, Behringer RR, Zhu C, Legerski RJ: Snm1-deficient mice exhibit accelerated tumorigenesis and susceptibility to infection. Mol Cell Biol 2005, 25:10071–10078.View ArticlePubMed
  64. Henriques JA, Brendel M: The role of PSO and SNM genes in DNA repair of the yeast Saccharomyces cerevisiae . Curr Genet 1990, 18:387–393.View ArticlePubMed
  65. Dronkert MLG, De Wit J, Boeve M, Vasconcelos ML, Van Steeg H, Tan TLR, Hoeijmakers JHJ, Kanaar R: Disruption of mouse Snm1 causes increased sensitivity to the DNA interstrand cross-linking agent mitomycin C. Mol Cell Biol 2000, 20:4553–4561.View ArticlePubMed
  66. Van de Ven M, Andressoo J-O, Holcomb VB, von Lindern M, Jong WMC, De Zeeuw CI, Suh Y, Hasty P, Hoeijmakers JHJ, van der Horst GTJ, Mitchell JR: Adaptive stress response in segmental progeria resembles long-lived dwarfism and calorie restriction in mice. PLoS Genetics 2006, 2:e192.View ArticlePubMed
  67. Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Edgar R: NCBI GEO: mining tens of millions of expression profiles – database and tools update. Nucleic Acids Research 2007, 35:D760-D765.View ArticlePubMed
  68. Parkinson H, Kapushesky M, Shojatalab M, Abeygunawardena N, Coulson R, Farne A, Holloway E, Kolesnykov N, Lilja P, Lukk M, Mani R, Rayner T, Sharma A, William E, Sarkans U, Brazma A: ArrayExpress – a public database of microarray experiments and gene expression profiles. Nucleic Acids Research 2007, 35:D747-D750.View ArticlePubMed
  69. Stevens JR, Doerge RW: Combining affymetrix microarray results. BMC Bioinformatics 2005, 6:57.View ArticlePubMed
  70. Hong F, Breitling R, McEntee CW, Wittner BS, Nemhauser JL, Chory J: RankProd: a bioconductor package for detecting differentially expressed genes in meta-analysis. Bioinformatics 2006, 22:2825–2827.View ArticlePubMed
  71. Larsson O, Wennmalm K, Sandberg R: Comparative microarray analysis. Omics 2006, 10:381–397.View ArticlePubMed
  72. Park T, Sung-Gon Y, Shin YK, Lee SY: Combining multiple microarrays in the presence of controlling variables. Bioinformatics 2006, 22:1682–1689.View ArticlePubMed
  73. Hill K, Boesch C, Goodall J, Pusey A, Williams J, Wrangham R: Mortality rates among wild chimpanzees. J Hum Evol 2001, 40:437–450.View ArticlePubMed
  74. Clark AG, Glanowski S, Nielsen R, Thomas PD, Kejariwal A, Todd MA, Tanenbaum DM, Civello D, Lu F, Murphy B, Ferriera S, Wang G, Zheng X, White TJ, Sninsky JJ, Adams MD, Cargill M: Inferring nonneutral evolution from human-chimp-mouse orthologous gene trios. Science 2003, 302:1960–1963.View ArticlePubMed
  75. Conover CA, Bale LK: Loss of pregnancy-associated plasma protein A extends lifespan in mice. Aging Cell, in press.
  76. Migliaccio E, Giorgio M, Mele S, Pelicci G, Reboldi P, Pandolfi PP, Lanfrancone L, Pelicci PG: The p66 shc adaptor protein controls oxidative stress resistance response and lifespan in mammals. Nature 1999, 402:309–313.View ArticlePubMed
  77. Holzenberger M, Dupont J, Ducos B, Leneuve P, Geloen A, Even PC, Cervera P, Le Bouc Y: IGF-I receptor regulates lifespan and resistance to oxidative stress in mice. Nature 2003, 421:182–187.View ArticlePubMed
  78. Kurosu H, Yamamoto M, Clark JD, Pastor JV, Nandi A, Gurnani P, McGuinness OP, Chikuda H, Yamaguchi M, Kawaguchi H, Shimomur I, Takayama Y, Herz J, Kahn CR, Rosenblatt KP, Kuro-o M: Suppression of aging in mice by the hormone klotho. Science 2005, 309:1829–1833.View ArticlePubMed
  79. Spencer CC, Howell CE, Wright AR, Promislow DE: Testing an 'aging gene' in long-lived Drosophila strains: increased longevity depends on sex and genetic background. Aging Cell 2003, 2:123–130.View ArticlePubMed
  80. Butler AA, LeRoith D: Minireview: tissue-specific versus generalized gene targeting of the igf1 and igf1r genes and their roles in insulin-like growth factor physiology. Endocrinology 2001, 142:1685–1688.View ArticlePubMed
  81. Wolkow CA, Kimura KD, Lee M-S, Ruvkun G: Regulation of C. elegans life-span by insulin-like signaling in the nervous system. Science 2000, 290:147–150.View ArticlePubMed
  82. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Atonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 2003, 4:249–264.View ArticlePubMed
  83. Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol 2004, 3:Article 3.
  84. Benjamini Y, Hochberg Y: Controlling the false discovery rate: a powerful and practical approach to multiple testing. J Roy Stat Soc B 1995, 57:289–300.
  85. Falcon S, Gentleman R: Using GOstats to test gene lists for GO term association. Bioinformatics 2007, 23:257–258.View ArticlePubMed
  86. Smid M, Dorssers LCJ, Jenster G: Venn Mapping: clustering of heterologous microarray data based on the number of co-occurring differentially expressed genes. Bioinformatics 2003, 19:2065–2071.View ArticlePubMed
  87. Amador-Noguez D, Zimmerman J, Venable S, Darlington G: Gender-specific alterations in gene expression and loss of liver sexual dimorphism in the long-lived Ames dwarf mice. Biochemical and Biophysical Research Communications 2005, 332:1086–1100.View ArticlePubMed
  88. Heishi M, Ichihara J, Teramoto R, Itakura Y, Hayashi K, Ishikawa H, Gomi H, Sakai J, Kanaoka M, Taiji M, Kimura T: Global gene expression analysis in liver of obese diabetic db/db mice treated with metformin. Diabetologia 2006, 49:1647–1655.View ArticlePubMed
  89. Recinos A III, Carr BK, Bartos DB, Boldogh I, Carmical JR, Belalcazar LM, Brasier AR: Liver gene expression associated with diet and lesion development in atherosclerosis-prone mice: induction of components of alternative complement pathway. Physiol Genomics 2004, 19:131–142.View ArticlePubMed
  90. Flowers MT, Groen AK, Oler AT, Keller MP, Choi Y, Schueler KL, Richards OC, Lan H, Miyazaki M, Kuipers F, Kendziorski CM, Ntambi JM, Attie AD: Cholestasis and hypercholesterolemia in SCD1-deficient mice fed a low-fat, high-carbohydrate diet. Journal of Lipid Research 2006, 47:2668–2680.View ArticlePubMed
  91. Niedernhofer LJ, Garinis GA, Raams A, Lalai AS, Robinson AR, Appeldoorn E, Odijk H, Oostendorp R, Ahmad A, van Leeuwen W, Theil AF, Vermeulen W, van der Horst GTJ, Meinecke P, Kleijer WJ, Vijg J, Jaspers NGJ, Hoeijmakers JHJ: A new progeroid syndrome reveals that genotoxic stress suppresses the somatotroph axis. Nature 2006, 444:1038–1043.View ArticlePubMed

Copyright

© Swindell. 2007

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.

Advertisement