In silico analysis of gene expression profiles in the olfactory mucosae of aging senescence-accelerated mice

We utilized high-density Affymetrix oligonucleotide arrays to investigate gene expression in the olfactory mucosae of near age-matched aging senescence-accelerated mice (SAM). The senescence-prone (SAMP) strain has a significantly shorter lifespan than does the senescence-resistant (SAMR) strain. To analyze our data, we applied biostatistical methods that included a correlation analysis to evaluate sources of methodologic and biological variability; a two-sided t-test to identify a subpopulation of Present genes with a biologically relevant P-value <0.05; and a false discovery rate (FDR) analysis adjusted to a stringent 5% level that yielded 127 genes with a P-value of <0.001 that were differentially regulated in near age-matched SAMPs (SAMP-Os; 13.75 months) compared to SAMRs (SAMR-Os, 12.5 months). Volcano plots related the variability in the mean hybridization signals as determined by the two-sided t-test to fold changes in gene expression. The genes were categorized into the six functional groups used previously in gene profiling experiments to identify candidate genes that may be relevant for senescence at the genomic and cellular levels in the aging mouse brain (Lee et al. [2000] Nat Genet 25:294-297) and in the olfactory mucosa (Getchell et al. [2003] Ageing Res Rev 2:211-243), which serves several functions that include chemosensory detection, immune barrier function, xenobiotic metabolism, and neurogenesis. Because SAMR-Os and SAMP-Os have substantially different median lifespans, we related the rate constant a in the Gompertz equation on aging to intrinsic as opposed to environmental mechanisms of senescence based on our analysis of genes modulated during aging in the olfactory mucosa.