- Open Access
A novel multi-tissue RNA diagnostic of healthy ageing relates to cognitive health status
Genome Biology volume 16, Article number: 185 (2015)
The Correspondence to this article has been published in Genome Biology 2019 20:152
The Correspondence to this article has been published in Genome Biology 2018 19:97
Diagnostics of the human ageing process may help predict future healthcare needs or guide preventative measures for tackling diseases of older age. We take a transcriptomics approach to build the first reproducible multi-tissue RNA expression signature by gene-chip profiling tissue from sedentary normal subjects who reached 65 years of age in good health.
One hundred and fifty probe-sets form an accurate classifier of young versus older muscle tissue and this healthy ageing RNA classifier performed consistently in independent cohorts of human muscle, skin and brain tissue (n = 594, AUC = 0.83–0.96) and thus represents a biomarker for biological age. Using the Uppsala Longitudinal Study of Adult Men birth-cohort (n = 108) we demonstrate that the RNA classifier is insensitive to confounding lifestyle biomarkers, while greater gene score at age 70 years is independently associated with better renal function at age 82 years and longevity. The gene score is ‘up-regulated’ in healthy human hippocampus with age, and when applied to blood RNA profiles from two large independent age-matched dementia case–control data sets (n = 717) the healthy controls have significantly greater gene scores than those with cognitive impairment. Alone, or when combined with our previously described prototype Alzheimer disease (AD) RNA ‘disease signature’, the healthy ageing RNA classifier is diagnostic for AD.
We identify a novel and statistically robust multi-tissue RNA signature of human healthy ageing that can act as a diagnostic of future health, using only a peripheral blood sample. This RNA signature has great potential to assist research aimed at finding treatments for and/or management of AD and other ageing-related conditions.
It is anticipated that novel genomic diagnostics that predict future health risks will help guide targeted preventative measures and enable the evaluation of individualized treatment strategies for many prevalent diseases of older age. So far, use of individual molecular biomarkers in healthy populations has offered modest performance [1, 2] compared with traditional, more integrated disease markers (e.g., blood pressure) or chronological age . For example, in people with cardiovascular disease, circulating cystatin C concentration, a parameter that estimates renal function, was related to 10-year mortality but was insufficient to predict cardiovascular deaths in healthy older subjects . Global RNA [5–9] and DNA methylation profiling [10–12] have been recently utilized to study the biology of chronological age. These existing signatures will incorporate influences of age-related disease and drug treatment. For example, Hannum et al. and Horvath et al. built distinct multi-tissue linear models, fitting age-related changes in DNA methylation with chronological age [13, 14]. These models have a statistical association with long-term health in the elderly  but the associations are not substantive enough to make it a practical diagnostic. In fact, as there are no molecular diagnostics of ‘healthy’ ageing status in humans, we hypothesized that a molecular profile may be useful at distinguishing people at risk for a variety of age-related diseases.
The shift in population demographics in the coming decades will mean that more than 1.2 billion people will be aged 65 years or older worldwide . Approximately 7 % of this population will have dementia, with at least 60 % of these having Alzheimer’s disease (AD). AD is the single largest healthcare cost  and there are currently no drug treatments that halt or cure it . Consensus is that only the earliest possible intervention is likely to significantly impact on AD and thus we need to identify those at greatest risk. The available validated diagnostics for AD are neither scalable for mass population screening nor sufficiently cost-effective to be practical . For example, brain imaging can provide clear evidence of neurodegeneration but is restricted to specialist centers  and an imaging-based public health screening program would not be affordable [19, 21]. There is a pressing need to stratify the older healthy population, using simple and cost-effective methods, to, for example, identify those appropriate to enrich clinical trials of novel AD treatments. Prototype blood diagnostics can be 75–85 % accurate at distinguishing AD patients from controls; however, these have not been validated using independently processed samples or have failed to replicate in independent studies . For example, blood-based protein signatures can diagnose mild cognitive impairment (MCI) and/or AD from controls in single studies [22–26], yet a common set of proteins has not been found across multiple studies. Further, the candidate AD marker proteins included cytokines and other markers of metabolic or cardiovascular disease  and thus these will not be clinically specific for AD when applied to older populations .
The expression of RNA is under genetic [29, 30], epigenetic  and environmental control [31, 32] and so the abundance of individual RNA molecules in blood cells reflects the integration of a variety of influences, whether or not blood directly interacts with a diseased organ. Thus, blood RNA [33–41] has been used to distinguish controls from MCI and/or AD, where variations in blood RNA expression should reflect the shared genetic, epigenetic  and environmental influences with the brain. For some prototype RNA diagnostics, the performances reported have been remarkably high (~95 %), but the same samples have been used during model building and validation [37, 38] and thus these represent examples of extreme over-fitting. In general there is always a danger that a classification model, when built using a specific set of cases and control samples from a single study, reflects unknown specific features of that particular cohort and thus is not generalizable.
In the present study we developed a RNA classifier of ‘healthy’ ageing starting with human muscle, with the hypothesis that this gene expression pattern may provide reliable genomic predictors for risk of age-related disease. We built the RNA classifier using human muscle global gene expression profiles because it has proven a useful tissue for predicting systemic physiological traits in humans  and because we can define healthy physiological status with ease . When the RNA classifier was related to cognitive health, this ‘healthy ageing gene score’ had the advantage of being hypothesis driven, and built using a paradigm and samples entirely distinct from clinical case–control samples. When applied to blood RNA, we established good validation for AD diagnosis and selectivity over common age-related pathologies. The results of the present study further support the idea that analysis of peripheral blood RNA would be a fruitful strategy for developing biomarkers of cognitive health and prove that a common healthy ageing gene-expression program is detectable across multiple tissues.
Results and discussion
Identification of a reproducible RNA signature for age of human muscle, brain and skin
Our objective was to discover a pattern of RNA expression that could be reliably used as a biomarker for ‘health status’ in older subjects — one that differed substantially in terms of ability to stratify health, and one that was more informative than chronological age. We applied machine-learning methods to RNA expression data to distinguish between healthy 25-year-old and healthy 65-year-old individuals. We took a simple classifier approach  without ad hoc a priori filtering to identify a consistent set of RNA markers of ageing across tissue types because standard differential expression is unable to provide a common multi-tissue set of discriminatory RNA molecules . We selected muscle tissue gene-chip profiles from 15 sedentary young and 15 sedentary older subjects with good aerobic fitness (Gene Expression Omnibus (GEO) accession [GSE59880]) [31, 44] and who were free of diabetes [42, 44]. Specifically, we utilized a k-nearest neighbor (kNN) classification approach because this captures data features that share non-linear interactions with robust performance  and is a method consistent with strategies recommended by the Microarray Quality Control consortium . This first data set — called the ‘training data-set’ — was used only once to select genes (Affymetrix probe-sets) and direction of gene expression change, and was then discarded from the project (Fig. 1). Expression differences of ~54,000 probe-sets were ranked using an empirical Bayesian statistic and a leave-one-out cross-validation (LOOCV) process (see “Materials and methods”). Probe-sets that targeted multiple genomic loci were removed and a 150 probe-set list, each gene having a nominal performance of 90 % or better, was selected for further study (Additional file 1). The extended list of probe-sets with a 70 % or better performance is also included in Additional file 1.
We checked that the 150 RNAs were not differentially expressed to any measurable extent in human muscle by exercise or a number of other common diseases that impact on skeletal muscle, using our previously published gene-chip data [8, 31, 44, 46]. We later confirmed this lack of association with lifestyle disease using a sensitive gene-set approach. Use of fully independent training and validation data sets allows for genuine external validation to be demonstrated (see “Materials and methods”). Using the ‘Campbell’ muscle data set [GEO:GSE9419]  as the samples of known identity, we demonstrated that additional young and old muscle samples selected from four additional muscle data sets (‘Trappe’ [GEO:GSE28422] , ‘Hoffman’ [GEO:GSE38718] , and ‘Kraus’ [GEO:GSE47969] and ‘Derby’ [GEO:GSE47881] ) could be classified with an average ~93 % accuracy (70–100 %) using only the 150 probe-sets selected at the start of the project. Substitution of the Campbell data set with the other muscle data sets worked equally as well. These data shared a common microarray platform (Affymetrix HGU133plus2) but, as we demonstrate below, the classifier remains robust in the face of alternative platforms. Receiver operating characteristic (ROC) curves for kNN = 5 demonstrating classifier performance for a number of tissue types are presented in Fig. 2.
Remarkably, the muscle-derived 150 RNA profile performed very well in classifying brain tissue by age. Using data from the HGU133Plus2 microarray platform for old and young samples of ectodermal origin (I, e, brain, n = 120)  we confirmed that the 150 RNA ‘healthy ageing’ genes selected in muscle could also distinguish the age of human brain one sample at a time, with a classification success rate up to 91 % (Fig. 2). Four brain regions were evaluated (postcentral gyrus, entorhinal cortex, hippocampus and superior frontal gyrus; [GEO:GSE11882]) and while they were confirmed disease-free by histopathology in the original study , unlike our muscle cohorts, their true functional status remains unknown. The postcentral gyrus samples were classified with 86 % sensitivity and 89 % specificity. In this cohort, older hippocampal regions were often misclassified using the 150 genes (33 % sensitivity) as ‘young’. This higher misclassification rate may relate to the substantial neurogenesis known to take place in the adult hippocampus or delays in tissue processing. We evaluated whether the 150 genes could accurately classify the age of tissue of mesodermal origin (skin) using gene expression data in a total of 279 human skin samples, of which there were up to three technical replicates per clinical sample . Notably, these data originated from a different technology platform (Illumina Human HT-12 V3, Array-express: E-TABM-1140), adding variability above that derived from a distinct tissue and potentially limiting the classification process. The two gene-chip technologies had 129 genes in common, and we observed excellent classification of human skin age [n = 131, area under the curve (AUC) = 0.85; Fig. 2]. The classification success was similar for all three replicates (71–78 raw classification success). Thus, the technical performance of the 150-gene healthy ageing classifier was excellent, providing accurate tissue classification despite inter-laboratory technical variation, different gene-chip platforms and antemortem issues. We were able, therefore, to conclude that we have identified a reliable multi-tissue RNA signature of healthy tissue ageing in humans, something that has not been previously demonstrated [8, 9].
A healthy ageing gene score that is distinct from chronological age and unrelated to lifestyle regulated phenotypes in the ULSAM study
In order to examine specificity for ‘healthy ageing’, we examined the relationship between the classifier genes, chronological age and markers of lifestyle-associated genes. We collapsed the expression pattern of all genes into a single score for each sample (see “Materials and methods”). The distribution of scores was examined for ~70-year-old males (subjects born in Uppsala within a 1-year period) and the gene ranking score was also correlated with markers of lifestyle-associated disease (Fig. 3). The gene expression profiles from 108 muscle samples from ~70-year-old male subjects from the Uppsala Longitudinal Study of Adult Men (ULSAM) cohort  were produced using Affymetrix arrays (Human Exon 1.0 ST Array). We ranked each subject for each of the 150 genes, taking the direction of gene expression change from the original classifier model into account (85 % down-regulated; see “Materials and methods”). We then converted the individual gene rankings into a summed median gene score for each subject. We demonstrated that despite all subjects being ~70 years of age at the time of the RNA sample, there was a very wide distribution in gene score (Fig. 3a). Thus, the healthy ageing gene score in muscle was very distinct from chronological age. The healthy ageing gene score was regressed against a variety of continuous clinical variables (variables listed in Additional file 2). The gene score at chronological age ~70 years was unrelated to conventional lifestyle regulated biomarkers (e.g., blood pressure, glucose, cholesterol, or renal function; Fig. 3b). This confirmed that the 150 gene expression markers were not reflecting a variety of lifestyle regulated biomarkers and diseases (e.g., exercise, diabetes) and tissue ‘healthy ageing status’ could not be derived from a simpler clinical biomarker.
Despite the limited sample size of the ULSAM cohort (n = 108), we were also able to demonstrate that subjects with the highest muscle healthy ageing gene score at age 70 years had significantly better renal function 12 years later (at age 82 years, p = 0.009). Remarkably, the healthy ageing gene score in muscle at ~70 years was also independently related to 20-year survival (p = 0.0295; Fig. S1a in Additional file 3) in a logistic regression model that included factors listed in Additional file 2). While this observation should be interpreted cautiously, to illustrate the temporal relationship between the healthy ageing gene score and death, we divided the gene score into quartiles and applied a Cox-regression model (Fig. S1b in Additional file 3) and found a significant difference between the first versus the fourth quartile (p = 0.04). In contrast to the healthy ageing gene score, a median gene rank score based on inflammatory gene (GO:0006954) or mitochondrial gene (GO:0005739) expression in muscle demonstrated no relationship with health or mortality (data not shown). The significant relationship between the healthy ageing gene score and organ function demonstrates that the gene expression pattern most similar to the healthy 65-year profile in the classifier model (i.e., the largest gene score in the ranking system) was associated with better health in the ULSAM cohort.
A greater healthy ageing gene expression score is associated with better cognitive health
Neurocognitive pathology (e.g., AD) becomes more pronounced with age and is often apparent in individuals who are otherwise healthy. Our analysis of the relationship between lifestyle factors and the healthy ageing gene score in the ULSAM cohort suggested that the gene score was robust to confounding effects of lifestyle disease. We next examined whether the healthy ageing gene score (median rank sum of the 150 RNA markers) was selectively useful in relation to identifying neurocognitive disease over lifestyle disease. To support this analysis, we utilized a large publically available gene-chip data set derived from healthy human brain samples of various ages . The BrainEac.org gene-chip resource  [GEO:GSE60862] comprises ten post-mortem brain regions from 134 subjects representing 1231 samples (Additional file 1). Using the same ranking approach as applied to the ULSAM cohort, the median sum of the rank score was calculated for each anatomical brain region (Fig. 4a). As before, in healthy older individuals the ‘age’ signature was ‘switched on’ (yielding a greater ranking score) compared with younger subjects. Regulation of the healthy ageing gene score increased in a distinct manner across individual healthy brain regions with chronological age, especially in the hippocampus (p = 0.00000002), as well as other regions (putamen, thalamus, substantia nigra, and the occipital, frontal, and temporal cortex regions (all at least p < 0.002 by Holm adjusted Mann–Whitney test).
Our primary hypothesis was that, compared with control subjects of similar chronological age and gender, patients with AD would have a lower median healthy ageing gene score, but the score would not distinguish diabetes or vascular (i.e. lifestyle influenced) disease patients from matched controls. We used two independent case–control studies of AD and two case–control studies of lifestyle disease with RNA profiles derived from blood. The first AD cohort has been previously used to study disease pathway changes in blood [41, 53] and we have deposited this data set (cohort 1 [GEO:GSE63060]) and a second analysis (cohort 2 [GEO:GSE63061]) at the GEO. We first used a maximum possible subset of subjects from each entire cohort, so that gender and chronological age could be precisely balanced (~70 years) remove these as potentially confounding factors. From cohort 1, 113 subjects were ranked for gene score, while 111 subjects were ranked in cohort 2 (Table 1). We checked for overlap between the 150 healthy ageing gene markers and previous genomic and genetic disease markers of AD (Additional file 1). Only three genes were in common and none were from previously validated AD diagnostics. Their inclusion or exclusion did not impact our analyses.
Blood RNA from AD case–control cohort 1 was profiled on Illumina HT-12 V3 bead-chips. We first mapped the appropriate probes from Affymetrix to Illumina, yielding 128 genes from the original 150-gene list. The relative median rank score for AD patients was significantly lower than for the age- and gender-matched controls (p = 0.004; Fig. 4b) based on Wilcoxon rank sum test. Blood RNA from the second AD case–control cohort was profiled on the Illumina HT-12 V4 platform and in this case 122 genes were in common with the 150-gene healthy ageing gene signature. As before, the median rank healthy ageing gene score for AD patients in cohort 2 was significantly lower than in the control group (p = 0.009; Fig. 4b). Furthermore, for both cohort 1 and cohort 2, the age-matched controls had a higher median gene score than subjects diagnosed with MCI (Fig. 4b; p = 0.00005 and p = 0.003 for cohorts 1 and 2, respectively). It is important to note that the control samples used for comparison with MCI overlapped with those used for comparison with AD and that the MCI analysis cannot, therefore, be considered a fully independent observation. As expected from the ULSAM analysis, the healthy ageing gene score was not related to diabetes or vascular disease status using blood profiles from 366 individuals (Additional file 4).
We formally evaluated whether the healthy ageing signature could act as a diagnostic for AD using ROC analysis and found that it had robust independent performance (AUC = 0.66–0.73; Fig. 5). We have previously published a whole blood RNA-based prototype AD diagnostic  consisting of 48 genes identified using machine learning methods applied to cohort 1 samples. We demonstrated that this prototype ‘RNA disease signature’ was independently validated in cohort 2 using LOOCV. Further, when we combined the two independently produced and validated gene expression classifiers we yielded an improved AD diagnostic (AUC = 0.73–0.86; Fig. 5) that matches best in class  for blood-based AD diagnostics validated using independent data, while our RNA-based analysis uses a technology platform more suited to reproducible high-throughput diagnostics.
Biological features of the healthy ageing diagnostic
We were interested in whether the healthy ageing diagnostic identified any particular biological processes that might be open to therapeutic targeting. The 150-gene list (Additional file 1) was evaluated using both Ingenuity pathway analysis and R-based gene ontology (GO) analysis. Ingenuity analysis (where a total of 127 genes were annotated in the database) revealed a few marginal functional associations (e.g., nervous system development genes) but these did not remain significant following Benjamini and Hochberg correction. The top ranked database network (genes with published interactions) was defined as ‘cell death and survival’ and contained 31 molecules. In Fig. 6a the density curves of p values for each one of 10,000 hypergeometric tests using a randomly sampled gene set (n = 150 in size) are plotted (black), along with the density curve of the p values from the healthy ageing 150-gene set (red). The profile of ontological enrichment in the healthy ageing diagnostic was not different from a random sample of 150 genes from the gene-chip, of which more than 99 % of the 54,000 probe-sets had no ability to discriminate tissue age in our training model. Manual searching of PubMed and the Online Mendelian Inheritance in Man (OMIM) database yielded some plausible connections with age-related and disease processes (Additional file 1) but such analysis is subjective and it cannot be concluded that these biological functions appear in the ‘healthy ageing’ diagnostic more than by simple proportionality. We did note that the 150 genes included some previously identified ‘ageing’ genes, for example, LMNA (linked with Hutchinson-Gilford progeria syndrome), Unc-13 homolog (UNC13C; linked with beta-amyloid biology), as well as COL1A1 (thought to change in skin ageing).
We also examined whether the 150 age-related genes were over represented at genomic loci using positional enrichment analysis  but found no significant associations. Using the top 670 genes from the first stage of the project (>70 % success in the training model) there were a number of significant findings (Additional file 1) with three genes originating from the top 150. In this analysis, 11q made a significantly greater contribution (adjusted p value = 0.005–0.007) to the enlarged prototype classifier than would be expected by chance (Fig. 6b), and there were a total of 15 genes from the 11q13 (ALDH3B1, CAPN1, CDC42EP2, CORO1B, LTBP3, NRXN2, PPP1R14B, RCE1, RCOR2, SART1, SYT12, and ZDHHC24; p = 0.0005) and 11q23 (FXYD2, SCN2B, and TMPRSS13; p = 0.0009) over-represented genomic locations. Interestingly, 11q23 is the location for age-related genetic interactions, namely the apolipoprotein A family [56, 57] as well as a region containing genetic association single nucleotide polymorphisms (SNPs) which modify the age of onset of colorectal cancer [58, 59]. Furthermore, 11q13 harbors SNPs associated with age of onset of renal cell carcinoma and prostate cancer and modulating age-related disease emergence by 5 years [60–62]. While the lack of an apparent specific biological dialogue may be considered disappointing, the extensive independent clinical results strongly support that the novel 150-gene healthy ageing signature is an important marker of healthy ageing in humans. Therefore, regulation of this gene expression program may in time reveal itself to be an important mechanism for maintaining human health and thereby a new opportunity for target development.
The molecular mechanisms that define healthy ageing remain elusive in both human and animal models . Many of the molecular mechanisms which extend the lifespan of laboratory animals also appear to extend health-span or disease-free ageing in these models . However, it remains unclear whether any of these mechanisms are central to human ageing [8, 9, 65] or define healthy ageing in humans. Our approach was novel because we first sought to define a set of genes associated with healthy ageing in ‘normal’ 65-year-old subjects rather than gene expression associated with disease or extreme longevity. This is an important distinction — ageing is thought to be a continuous physiological process that could be expected to have a gene expression signature distinct from lifestyle related (e.g., type II diabetes) or mutation driven (e.g., cancer) pathologies, thus explaining its independent prognostic rather than specific diagnostic capacity. It is however of potential importance that regulation of the healthy ageing signature in human brain is most evident in those regions associated with neurodegeneration. In contrast, it is thought that the cerebellar cortex is not subject to substantial age-related anatomical changes  and this was consistent with our new model of human healthy ageing (Fig. 4a).
The 65-year-old subjects used to build the RNA model were in good health despite leading a normal sedentary lifestyle. Rather than using individual differential expression values to define discriminatory genes, we selected a group of genes that would act together to make a ‘majority vote’. Indeed, we were able to demonstrate that the 150 healthy ageing genes are consistently modulated in several tissue types, but to very differing degrees in people of the same chronological age (e.g., Fig. 3a). Thus, the healthy ageing gene score fulfilled the first main criteria for being a novel diagnostic of healthy (or biological) ageing. Including the ULSAM analysis (males only), we have demonstrated in three independent clinical cohorts that greater healthy ageing gene score associates with better health in men and women, suggesting that promotion of this gene expression profile may be beneficial and could reflect an adaptive compensatory response. The present RNA diagnostic could be used to facilitate the evaluation of anti-ageing-related treatments in middle-aged humans, screen for long-term safety during drug development, or augment clinical decision-making that currently inputs chronological age rather than ‘biological’ age into treatment algorithms. Future efforts should focus on discovering strategies to modulate the healthy ageing gene signature to establish if it is causally determining health or just acting as a robust biomarker of a more complex set of molecular interactions.
The multi-tissue healthy ageing gene score is predictive of health in older subjects
Exceptional longevity is driven by a measurable genetic contribution [67, 68], while being active and healthy at age 65 years is a more common occurrence, likely to reflect complex molecular factors [64, 69], and is less obviously linked to only variations in DNA sequence. We profiled RNA from healthy members of the ULSAM cohort at age 70 years and analyzed follow-up data over two decades. In 1992, these 70-year-old Swedish men had normal levels of physical activity ‘for their age’ and most demonstrated longevity to 90 years, which is not exceptional in the Swedish population . The healthy ageing score demonstrated a four-fold range (Fig. 3a) while chronological age varied by no more than one year across the group. A greater gene score was associated with better cognitive function, and better renal function across a 12 year span and both cognitive decline and renal function are important determinants of all-cause mortality [71, 72]. A concurrent reduction in cognitive and renal function is clinically observed, suggesting both are subject to a general age-related decline in organ function . It is perceived that lifestyle-regulated diseases, such as type II diabetes, causally increase AD . This relationship did not appear causal when the rates of emergence of diabetes (or not) and number of emerging cases of AD were compared in the Framingham or Baltimore Age studies [75, 76]. This suggests that risk of developing type II diabetes and AD may share some concurrent risk factors, e.g., aerobic capacity ; a physiological capacity defined by a large genetic and gene–environmental interaction . Additionally, type II diabetes and AD may share epigenetic or genetic risk factors. Interpretation of such associations is further complicated by the interaction between type II diabetes, vascular disease, and other types of dementia (which complicate the diagnosis of AD).
Neurological decline is predicted to contribute substantially to the economic burden of healthcare in the coming decades. AD is a multi-factorial disease  with around 22 genetic loci potentially associated with disease risk or progression of symptoms. The strongest and most reproducible genomic association, APOE-ε4, is a modifier of risk, contributing to the variance in age of onset of the disease by 3.7 % . The remaining approximately nine reproducible risk loci for late-onset AD (the most common form) contribute a further 2.2 % of the variance in age of onset . In short, these DNA sequence variants will not be clinically useful for diagnosing or managing AD or even assessing risk in the majority of people. Differential gene expression analysis and molecular classification have found disease-related RNA markers of AD, using patient materials to build the model . However, such diagnostics can be biased by unknown features of the training data (the data used to select the RNA markers). In contrast our healthy ageing genes were selected via a hypothesis-driven strategy that then relied on a validation process that included seven independent tissue cohorts including multiple RNA detection technologies (so ruling out some unknown technology platform bias). Thus, our healthy ageing gene expression signature has the key advantage of being a signature built using a paradigm and samples entirely distinct from AD case–control samples.
The healthy ageing gene score allowed us to demonstrate that patients diagnosed with AD have an altered healthy ageing RNA expression signature in blood that demonstrates significant association with disease. Furthermore, the muscle or blood gene score was unrelated to lifestyle diseases such as type II diabetes and thus may be more clinically specific than earlier AD biomarkers [20, 22–26, 41, 81], most of which did not replicate in independent clinical studies. We were able to provide independent validation for our earlier AD-related ‘disease’ diagnostic ; like many AD disease biomarkers , however, it includes pro-inflammatory markers and oxidative stress, features that can be common to several diseases and thus it may not be specific in clinical practice. Nevertheless, when we combined the Lunnon et al.  AD biomarker (even after removing the eight genes we found to be regulated in blood by diabetes or vascular disease) with the healthy ageing genes we yielded an improved diagnostic for AD over and above either diagnostic alone (Fig. 5). Ultimately, formal diagnosis of AD will continue to rely on a combination of diagnostics, including invasive cerebrospinal fluid sampling, positron emission tomography (PET) imaging and magnetic resonance imaging (MRI). However, given the scale of screening required (e.g., more than 1 million people in 2015/2016 to deliver sufficient numbers of at-risk subjects for AD clinical trials ) a blood-based diagnostic will be extremely useful for pre-screening ahead of invasive and costly follow-up analysis. Enrichment of prevention trials with asymptomatic people most at risk for AD is required to ensure that event rates are sufficiently high to evaluate the multitude of drug trials being considered for AD .
Like many genomic diagnostics, the full clinical utility of ours will only emerge when combined with additional data and clinical insight. While we could also demonstrate that patients with MCI had a significantly lower healthy ageing gene score it remains to be shown that this can be converted into a diagnostic for future cognitive health (i.e., a blood sample from older healthy subjects or those with recently diagnosed MCI combined with 5–10 years of follow-up data to prove they did or did not develop AD). Epidemiological efforts to build long-range (~36 year) forecasting of dementia risk (AD or vascular) using clinical demographics  (CAIDE score) provide assessment of risk at middle-age (~45 years and over) and can assign patients into a low (9 %) or high (29 %) risk with un-validated ROC AUC = 0.74. In diabetes patients, age is by far the more powerful predictor of future dementia rather than severity of the diabetes measured using glycosylated hemoglobin A1 (HbA1)  and it will be informative to replace age with our healthy ageing gene diagnostic for many conditions. These examples highlight that, clinically, various decision trees exist and our healthy ageing score could be integrated to help decide which middle-aged subjects could be offered entry into a preventative clinical trial many years before the clinical expression of AD.
As can be observed in Fig. 5, we obtained an independently validated ROC AUC of 0.73 using default classification settings. This is not an optimized or over-fitted model, in that it is most likely possible to tune ROC parameters to yield an improved performance using the same list of genes or a subset thereof. Some previous authors have reported ROC AUC scores of >0.8, but as mentioned these do not represent valid scores, being derived from data over-fitted to a single data set [33–40]. High ‘scores’ encompass shared technical variance in the ‘test’ and ‘validation’ data. For example, a microRNA dementia diagnostic (relying on 2–17 microRNA real-time PCR assays) was validated in samples used to build the initial model, yielding inflated specificity and sensitivity values . Other practical factors must be considered, such as the complexity of the laboratory test and costs. Cheng et al.  used a complex process to isolate serum exosome microRNAs and a split-cohort partial validation approach. Their best AD model (16 microRNAs) was 87 % sensitive and 77 % specific but could not diagnose MCI. DiaMir Inc. claimed 95 % specificity and sensitivity for MCI (their model did not work in AD) but this failed to replicate in a second study .
During the past 5 years several projects have worked with whole-blood mRNA and produced 20–225-gene assays to classify samples during the training phase, comparing controls with patients (MCI or AD). The RNA classifier from DiaGenic A/S used TaqMan assays and this remains the only replicated blood-based AD diagnostic performing to a similar level as the present study , while there were no genes in common with our own ‘AD disease’ or ‘healthy ageing’ RNA signatures. To date, few blood-based protein and metabolite diagnostics have been replicated using a fully independent process. Doecke et al.  combined a panel of eight protein markers with age, gender, and APOE genotype in the Australian Imaging, Biomarkers and Lifestyle (AIBL) cohort and found an AUC of 0.84 using a subset of the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort. However, their protein markers were a priori known to be regulated in the ADNI cohort, and it remains unclear how many of the eight protein markers contribute to the model, including age, gender, and genotype, in the validation process.
Comparison with known markers of human ageing or longevity
Other approaches have been utilized in humans to understand the molecular determinants of human ageing, but not ‘healthy’ ageing. Genome-wide association analysis has shown 281 DNA variants linked with exceptional longevity, and collectively explaining 17 % variance in humans  with an AUC value of 0.6. This remains to be replicated and this list of genes did not overlap with our healthy ageing gene list. In addition, long-lived humans appear to have a similar genetic burden for common DNA disease variants, suggesting the human exceptional longevity model may not be reflective of the processes that determine average longevity . There have been several linear molecular models of chronological age [5, 8, 11, 14] but the variance captured by these across chronological age is limited and the disease status of samples used to build or validate the model unclear; thus, it is uncertain if such models reflect ageing or age-related disease and drug treatment. There was no overlap between the genes in our healthy ageing RNA classifier and the quasi-linear methylation model derived by Horvath et al. .
We found four genes in common between our healthy ageing RNA classifier and the two gene lists identified by Hannum et al. in separate DNA methylation models of ageing (n = 94 and n = 326): one gene from their primary model (PKM2) and three genes from their RNA methylation association analysis (ANKRD13B, RUNX3, and TCF3) . Neither the Horvath nor the Hannum models generate sufficient distinction from chronological age to provide a useable ‘size effect’ when considering longevity . Passtoors et al.  reported that a set of 21 RNA molecules ‘marked out’ familial longevity in blood RNA, but this was a weak correlation with no discriminatory capacity as a diagnostic, possibly because it reflects a mixture of ageing, disease, and drug therapy. Furthermore, none of those age-related blood RNA changes were consistently correlated with age in human brain or muscle [8, 85], indicating that these 21 RNAs do not represent universal markers of human ageing (they were also not part of our 150 healthy ageing gene list). We did not note any significant ontology pathway enrichment within our healthy ageing diagnostic gene lists (Fig. 6a). Thus, we cannot neatly place the genes that contribute to the healthy physiological age diagnostic into a convenient canonical signaling pathway.
Materials and methods
Informed consent was obtained from all volunteers and ethical approval received from Institutional Research Ethics Committee as reported in primary clinical publications [8, 9, 31, 44, 46–50, 52] and all studies included in this work were conducted under the auspices of the declaration of Helsinki. For new gene-chip tissue profiles and hence new GEO deposits, the Institutional Research Ethics Committee approvals were as follows: ULSAM (Regional Ethical Review Boards Uppsala Ethical applications 09-154M/2010/400), STRRIDE (Duke Medical School IRB, Pro00012628) and AddNeuroMed/DCR (SLaM/IOP 30/07/2006/SLaM/IOP 30/04/2008).
A summary of the analysis strategy was as follows. The first aim was to generate a reliable RNA classifier of healthy older muscle tissue (healthy ageing gene score). We utilized k-nearest neighbor (kNN) classification methods because they capture data features that share non-linear interactions and have robust performance using methods consistent with the Microarray Quality Control Consortium . The probe-set level intensities of each set of independent microarrays were normalized using the Robust Multi-array Analysis (RMA) method implemented within the R statistical software environment using the ‘affy’ package, and then scaled and centered (Bioconductor project [86, 87]). When Affymetrix gene-chips originated from independent laboratories, we used Frozen Robust Multi-array Analysis (fRMA) [88, 89]. Having identified a healthy ageing gene score comprising 150 RNA markers (probe-sets), we established that these 150 RNAs could reliably classify multiple independent sets of human muscle and brain tissue using external validation. External validation uses independent training and validation data sets. Finally, we examined if the healthy ageing gene score in blood was related to cognitive health, alone or in combination with our prototype blood marker of early AD. Figure 1 presents the project analysis scheme.
Production and independent external validation of the healthy ageing gene score
We identified 150 RNA markers of muscle ageing using samples [31, 44] and gene-chip profiles [GEO:GSE59880] from 15 young (aged 19–28 years) and 15 older subjects (aged 59–77 years) free from metabolic and signs of cardiovascular disease and validated this observation in more than 500 independent samples. The older subjects were sedentary (did not do any regular sport/exercise) but nevertheless were free of diabetes and had good levels of aerobic capacity, a marker of general health into older age . The RNA markers were selected using a nested-loop, holding out two arrays at any one time to estimate two parameters from the data. The first of these was the conventional classification result; i.e., was the ‘unknown’ sample correctly classified, yes or no? The second parameter was used to calculate the performance of the probe-sets contributing to the decision. We selected 200 probe-sets during each of the inner-most loops by ranking gene expression differences using an empirical Bayesian statistic (implemented as eBayes in the ‘limma’ package) . Following iterative assessment of all probe-sets and all samples, involving ~180,000 permutations, a list of ~800 probe-sets was identified as having good performance (>70 % correct). We removed probe-sets that targeted multiple genomic loci and selected the top ranked 150 probe-sets (involved in >90 % correct decisions) for further study. Classifier performance was assessed using ROC analysis and the R package ROCR .
We implemented fully independent external validation of the 150-probe-set healthy ageing classifier, a process that requires both independent ‘known samples’ to define the expression space and independent test gene-chips . When combined with LOOCV methods, this represents a gold standard approach to validation of a classification model. A new set of young and old muscle profiles (selected from the Campbell data set; n = 66 chips ; [GEO:GSE9419]) was used to represent the new expression space of known samples. We then carried out evaluation of sets of independent gene-chip profiles from young and old human muscle (all Affymetrix U133+2) normalized using fRMA. The various fully independent samples were obtained from GEO or produced from our own clinical samples . For each dataset a subset of samples were selected to belong to either the young (~25 years) or older group (~65 years) from a larger collection of samples. The sets of young and older samples were selected from the Trappe  [GEO:GSE28422]; n = 48), Hoffman  [GEO:GSE38718]; n = 22), Derby  [GEO:GSE47881]; n = 26) and Kraus [GEO:GSE47969]; n = 33) data sets. For the Kraus data set total RNA was extracted from frozen muscle biopsy samples (vastas lateralis) using TRIzol reagent and in vitro transcription was performed using the Bioarray high yield RNA transcript labeling kit (P/N 900182, Affymetrix, Inc.) as previously described . For all data sets, arrays were examined using hierarchical clustering and normalized unscaled standard error (NUSE). In cases where we identified a small number of gene-chips (two to three) that had evidence of technical defects, these were removed prior to any analysis.
To assess if human brain and skin also demonstrated the same 150 age-related gene expression signature as healthy older muscle, we used young and old samples from the brain-bank array source (n = 120; [GEO:GSE11882]) and the MuTHER cohort skin data set (n = 279, which includes a subset of three replicates, n = 131, n = 124, and n = 24). The skin data were produced using the Illumina Human HT-12 V3 Bead chip (Array-express: E-TABM-1140) and log-2 transformed signals were normalized using quantile normalization. The 150 Affymetrix probe-sets were mapped to the Illumina platform (giving 129 probes). Due to differences in gene-chip technology, a LOOCV approach was used to classify the age of each skin sample, using only the 150 probes selected at the start of the project. For skin, individuals aged 45 years or less were defined as young, and those aged 70 years or older as old to ensure balanced numbers of young and old samples existed to fairly assess the classifier performance. The three sets of technical replicates were analyzed separately and confirmed the intra-study technical reproducibility of the classifier using repeated RNA profiles of a single clinical sample (data not shown).
The healthy muscle ageing gene score differs substantially from chronological age
We used a set of tissue samples from a birth cohort of men, such that the same chronological age (~70 years) could be contrasted with the variation in healthy ageing gene score. The ULSAM cohort comprises men born in 1920–1924 and living in Uppsala, Sweden and was used to compare a constant chronological age (and similar environment) with the healthy muscle age gene score across individuals . Dual-energy X-ray absorptiometry (DXA) scan measurements were performed during the last decade of the study and muscle mass status varied between −15 % to +10 % between age 70 years and 88 years and was unrelated to physical activity scores (recorded at 82 years and 88 years of age, with 80 % being recorded as being moderately active). Renal function was estimated using cystatin C, which is a marker of glomerular filtration rate . We had access to 129 skeletal muscle biopsies taken at age 70 years (in 1992) and we processed these in 2012 with the majority having excellent NUSE plot profiles. Total RNA was extracted from frozen muscle biopsy samples (vastas lateralis) using TRIzol reagent as previously described . A total of 113 samples provided sufficient RNA and 50 ng total RNA was amplified using Ambion’s WT expression kit to produce cDNA. The cDNA was fragmented and labeled with GeneChip WT Terminal labeling kit (Affymetrix, Inc.). Unincorporated nucleotides from the in vitro transcription reaction were removed using an RNeasy column (QIAGEN Inc.). Hybridization, washing, staining, and scanning of the arrays were performed according to the manufacturer’s instructions (Affymetrix, Inc. Santa Clara, USA).
One hundred and eight samples passed gene-chip quality control procedures (see above). A cumulative gene ranking-based score was calculated using each of the 150 gene expression values for each of the 108 male subjects and the final score was compared in a linear fashion with a number of clinical parameters. For an RNA down-regulated in the original training classification data set (i.e., down-regulated between 25 years to 65 years) the ULSAM subject with the highest expression was assigned a score of 1 and the subject with the lowest expression 108. For genes up-regulated in the original age classification model, the opposite strategy was used. Thus, both feature selection (genes) and direction of regulation were taken from the original model. The median sum of these rank scores (reflecting the 150 probe-sets) was calculated and that represented the healthy ageing gene score for each individual in the ULSAM cohort. Median rank ensured each gene provided equal weighting and regression analysis was used to study the variation in gene score in these men, all of who had approximately the same chronological age.
The relationship between the gene score at age 70 years and a number of clinical features was carried out using multi-factor models. Model selection was executed using a forward selection approach, with p > 0.1 as stop criterion (backwards elimination yielded identical results). Clinical variables, previously reported , were added to the baseline model one at a time, and selected based on p value  (Additional file 2). Over the observation period mortality rate was 18 % (19 events) and the relationship between mortality and gene score was analyzed as a continuous variable. Both the Cox-regression and the logistic regression model were implemented in R. For the Cox model we used the latest ‘survival package’ whereas the logistic regression model was estimated using the glm (generalized linear model) function and ‘logit’ model, which models the log odds of the outcome as a linear combination of the predictor variables. For the Kaplan–Meier plots, gene score was divided into quartiles and the plot was produced using the plot-survfit function in the survival package. All three approaches yielded consistent results.
Relationship between the healthy ageing gene score in blood and disease status using age- and gender-matched case–control analysis
Demonstration that the healthy ageing gene score was clearly demonstrable in neuro-muscular tissue suggested that it might also relate to cognitive health. Indeed, to provide additional support for the observations in human brain, we used the BrainEac.org gene-chip resource , which comprises ten post-mortem brain regions from 134 subjects representing 1,231 samples [GEO:GSE60862] (Additional file 1). For each brain region, and for a down-regulated gene in the original model, the subject with highest expression was assigned a score of 1 and the subject with the lowest expression was assigned a score of 134 (upper and lower score depends on total number of samples; Additional file 1). The median sum of the rank score was calculated for each anatomical brain region in the same manner as described above, with the ULSAM cohort. The healthy ageing gene score differed across the brain regions with chronological age, as determined by a Kruskal–Wallis test. A Kruskal–Wallis test was used as we were comparing unequal observations per age group (Additional file 1). Post hoc Mann–Whitney test with correction for multiple comparisons (Holm) was used to confirm regulation of the ageing signature genes in each region.
We used blood RNA profiles from subjects from the AddNeuroMed consortium, a large cross-European AD biomarker study and a follow-on Dementia Case Register (DCR) cohort in London. Patient selection, design, and clinical data have been reported previously [53, 97]. AD data sets have been deposited under accessions [GEO:GSE63060] and [GEO:GSE63061]. A summary of the cohort characteristics can be found in Table 1. Briefly, subjects were excluded from the study if they had neurological or psychiatric illness other than AD, unstable systematic illness or organ failure, or a geriatric depression rating scale score ≥ 4/5. AD was diagnosed using the National Institute of Neurological and Communicative Disease and Stroke and Alzheimer’s disease (NINCDS-ADRDA) and Diagnostic and Statistical Manual of Mental Disorders (DSM-IV) criteria for possible or probable AD. All MCI subjects reported problems with memory, corroborated by an informant, but had normal activities of daily living as specified in the Petersen’s criteria for amnestic MCI [41, 97]. All subjects underwent a structured interview and a battery of neuropsychological assessments, including the mini-mental state examination (MMSE). Control and MCI subjects were further assessed using the CERAD battery and detailed information on subject recruitment and assessments can be found in other published studies describing the AddNeuroMed consortium [54, 97]. RNA was obtained from whole venous blood and it was collected from the subjects who had fasted 2 hours prior to collection into a PAXgene™ Blood RNA tube (Becton & Dickenson, QIAGEN Inc., Valencia, CA, USA). The tubes were frozen at −20 °C overnight prior to long-term storage at −80 °C. RNA was extracted using PAXgene™ Blood RNA Kit (QIAGEN) according to the manufacturer’s instructions.
We used two independently produced gene-chip data sets from the AddNeuroMed/DCR consortia, one data-set produced in a UK gene-chip facility and another produced in the USA. Gene expression data was produced using Illumina Human HT-12 v.3 Expression BeadChips for the first case–control study (USA; cohort 1) and Illumina Human HT-12 v4 Expression BeadChips for the second case–control study (UK; cohort 2). cDNA was synthesized from 200 ng total RNA using the TotalPrep™ RNA Amplification Kit (Ambion), which was followed by amplification and biotinylation of cRNA and hybridization. The expression data were first transformed using variance stabilization and then quantile normalized using the LUMI package in R. The 150 probe-sets were mapped from the Affymetrix platform to the Illumina platform. For our primary analysis, control subjects were matched in a manner that created the largest possible group with the same chronological age and gender balance as the AD or MCI groups. Thus, our analysis was carried out on a subset of subjects deposited at the GEO, with each case–control group having a similar median chronological age as the ULSAM cohort. In total 297 samples were utilized (batch 1 Control = 67, MCI = 39, AD = 49; batch 2 Control = 72, MCI = 30, AD = 40). Retrospective inclusion of the entire cohort (n = 717) did not alter the outcome of our analysis.
We also utilized two additional large gene-chip clinical studies; one comparing blood RNA in type II diabetes with control  and the other from our laboratory comparing blood RNA in people with and without coronary artery disease . For each case–control comparison the ranking metric was computed in the exact same manner as for the ULSAM subjects and AD patients (see above). A Wilcoxon rank sum test from the R stats package was used to test if the median gene score ranks between groups were significantly different or not. For data presentation, ranking scores were scaled to the total number of samples being ranked to ensure each data plot was on the same scale.
The bioinformatics tool Ingenuity Pathway Analysis (IPA)  was used to explore the biology of the age classifier genes. HUGO gene name identifiers were uploaded into IPA and queried against the verified IPA knowledge database. To establish the GO profile of the 150 genes, we generated a null distribution of GO enrichment p values by randomly sampling 10,000 lists of 150 probe-sets from the hgu133plus2 chip and testing each list for the GO term molecular function using the GOstats package in R. The entire population of probe-sets on the hgu133plus2 microarray was used as the background population for these tests. The resulting p values for each tested probe-set list were corrected using the method of Benjamini and Hochberg. The 150 healthy ageing genes were then tested for GO term molecular function and the p values Benjamini and Hochberg corrected. Positional gene enrichment analysis was used to identify whether the classification genes (or the classifier network genes) were significantly enriched within given chromosomal regions  as previously implemented .
area under the curve
Dementia Case Register
Frozen Robust Multi-array Analysis
Gene Expression Omnibus
Ingenuity Pathway Analysis
mild cognitive impairment
normalized unscaled standard error
Robust Multi-array Analysis
receiver operating characteristic
single nucleotide polymorphism
Uppsala Longitudinal Study of Adult Men
Wallentin L, Zethelius B, Berglund L, Eggers KM, Lind L, Lindahl B, et al. GDF-15 for prognostication of cardiovascular and cancer morbidity and mortality in men. PLoS One. 2013;8:e78797.
Daniels LB, Clopton P, Laughlin GA, Maisel AS, Barrett-Connor E. Growth-differentiation factor-15 is a robust, independent predictor of 11-year mortality risk in community-dwelling older adults: the Rancho Bernardo Study. Circulation. 2011;123:2101–10.
Rohatgi A, Patel P, Das SR, Ayers CR, Khera A. Association of growth differentiation factor-15 with coronary atherosclerosis and mortality in a young, multiethnic population: observations from the dallas heart study. Clin Chem. 2014;58:172–82.
Inker LA, Schmid CH, Tighiouart H, Eckfeldt JH, Feldman HI, Greene T, et al. Estimating glomerular filtration rate from serum creatinine and cystatin C. N Engl J Med. 2012;367:20–9.
Passtoors WM, Boer JM, Goeman JJ, Akker EB, Deelen J, Zwaan BJ, et al. Transcriptional profiling of human familial longevity indicates a role for ASF1A and IL7R. PLoS One. 2012;7:e27759.
Passtoors WM, Beekman M, Deelen J, van der Breggen R, Maier AB, Guigas B, et al. Gene expression analysis of mTOR pathway: association with human longevity. Aging Cell. 2013;12:24–31.
Gheorghe M, Snoeck M, Emmerich M, Bäck T, Goeman JJ, Raz V. Major aging-associated RNA expressions change at two distinct age-positions. BMC Genomics. 2014;15:132.
Phillips BE, Williams JP, Gustafsson T, Bouchard C, Rankinen T, Knudsen S, et al. Molecular networks of human muscle adaptation to exercise and age. PLoS Genet. 2013;9:e1003389.
Glass D, Viñuela A, Davies MN, Ramasamy A, Parts L, Knowles D, et al. Gene expression changes with age in skin, adipose tissue, blood and brain. Genome Biol. 2013;14:R75.
Christensen BC, Houseman EA, Marsit CJ, Zheng S, Wrensch MR, Wiemels JL, et al. Aging and environmental exposures alter tissue-specific DNA methylation dependent upon CpG island context. PLoS Genet. 2009;5:e1000602.
Horvath S. DNA methylation age of human tissues and cell types DNA methylation age of human tissues and cell types. Genome Biol. 2013;14:R115.
Bell JT, Tsai P-C, Yang T-P, Pidsley R, Nisbet J, Glass D, et al. Epigenome-wide scans identify differentially methylated regions for age and age-related phenotypes in a healthy ageing population. PLoS Genet. 2012;8:e1002629.
Horvath S, Zhang Y, Langfelder P, Kahn RS, Boks MP, van Eijk K, et al. Aging effects on DNA methylation modules in human brain and blood tissue. Genome Biol. 2012;13:R97.
Hannum G, Guinney J, Zhao L, Zhang L, Hughes G, Sadda S, et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell. 2013;49:359–67.
Marioni RE, Shah S, McRae AF, Chen BH, Colicino E, Harris SE, et al. DNA methylation age of blood predicts all-cause mortality in later life. Genome Biol. 2015;16:25.
Harper S. Economic and social implications of aging societies. Science. 2014;346:587–91.
Bynum JPW. The long reach of alzheimer’s disease: Patients, practice, and policy. Health Aff. 2014;33:534–40.
Salloway S, Sperling R, Fox NC, Blennow K, Klunk W, Raskind M, et al. Two phase 3 trials of bapineuzumab in mild-to-moderate Alzheimer’s disease. N Engl J Med. 2014;370:322–33.
Biasutti M, Dufour N, Ferroud C, Dab W, Temime L. Cost-effectiveness of magnetic resonance imaging with a new contrast agent for the early diagnosis of Alzheimer’s disease. PLoS One. 2012;7:e35559.
Laske C, Sohrabi HR, Frost SM, López-de-Ipiña K, Garrard P, Buscema M, et al. Innovative diagnostic tools for early detection of Alzheimer’s disease. Alzheimer’s Dement. 2014;11:1–18.
Furiak NM, Klein RW, Kahle-Wrobleski K, Siemers ER, Sarpong E, Klein TM. Modeling screening, prevention, and delaying of Alzheimer’s disease: an early-stage decision analytic model. BMC Med Inform Decis Mak. 2010;10:24.
Ray S, Britschgi M, Herbert C, Takeda-Uchimura Y, Boxer A, Blennow K, et al. Classification and prediction of clinical Alzheimer’s diagnosis based on plasma signaling proteins. Nat Med. 2007;13:1359–62.
Hye A, Riddoch-Contreras J, Baird AL, Ashton NJ, Bazenet C, Leung R, et al. Plasma proteins predict conversion to dementia from prodromal disease. Alzheimers Dement. 2014;44:1–9.
O’Bryant SE, Xiao G, Barber R, Huebinger R, Wilhelmsen K, Edwards M, et al. A blood-based screening tool for Alzheimer’s disease that spans serum and plasma: findings from TARC and ADNI. PLoS One. 2011;6:e28092.
Hu WT, Holtzman DM, Fagan AM, Shaw LM, Perrin R, Arnold SE, et al. Plasma multianalyte profiling in mild cognitive impairment and Alzheimer disease. Neurology. 2012;79:897–905.
Sattlecker M, Kiddle SJ, Newhouse S, Proitsi P, Nelson S, Williams S, et al. Alzheimer’s disease biomarker discovery using SOMAscan multiplexed protein technology. Alzheimers Dement. 2014;10:724–34.
Doecke JD, Laws SM, Faux NG, Wilson W, Burnham SC, Lam C-P, et al. Blood-based protein biomarkers for diagnosis of Alzheimer disease. Arch Neurol. 2012;69:1318–25.
Dudley JT, Butte AJ. Identification of discriminating biomarkers for human disease using integrative network biology. Pac Symp Biocomput. 2009;27–38.
Schwanhäusser B, Busse D, Li N, Dittmar G, Schuchhardt J, Wolf J, et al. Global quantification of mammalian gene expression control. Nature. 2011;473:337–42.
Westra H-J, Franke L. From genome to function by studying eQTLs. Biochim Biophys Acta. 2014;1842:1896–902.
Keller P, Vollaard NBJ, Gustafsson T, Gallagher IJ, Sundberg CJ, Rankinen T, et al. A transcriptional map of the impact of endurance exercise training on skeletal muscle phenotype. J Appl Physiol. 2011;110:46–59.
Larrouy D, Barbe P, Valle C, Déjean S, Pelloux V, Thalamas C, et al. Gene expression profiling of human skeletal muscle in response to stabilized weight loss. Am J Clin Nutr. 2008;88:125–32.
Rye PD, Booij BB, Grave G, Lindahl T, Kristiansen L, Andersen H-M, et al. A novel blood test for the early detection of Alzheimer’s disease. J Alzheimers Dis. 2011;23:121–9.
Booij BB, Lindahl T, Wetterberg P, Skaane NV, Sæbø S, Feten G, et al. A gene expression pattern in blood for the early detection of Alzheimer’s disease. J Alzheimers Dis. 2011;23:109–19.
Fehlbaum-Beurdeley P, Sol O, Désiré L, Touchon J, Dantoine T, Vercelletto M, et al. Validation of AclarusDxTM, a blood-based transcriptomic signature for the diagnosis of Alzheimer’s disease. J Alzheimers Dis. 2012;32:169–81.
Roed L, Grave G, Lindahl T, Rian E, Horndalsveen PO, Lannfelt L, et al. Prediction of mild cognitive impairment that evolves into Alzheimer’s disease dementia within two years using a gene expression signature in blood: a pilot study. J Alzheimers Dis. 2013;35:611–21.
Kumar P, Dezso Z, MacKenzie C, Oestreicher J, Agoulnik S, Byrne M, et al. Circulating miRNA biomarkers for Alzheimer’s disease. PLoS One. 2013;8:e69807.
Leidinger P, Backes C, Deutscher S, Schmitt K, Mueller SC, Frese K, et al. A blood based 12-miRNA signature of Alzheimer disease patients. Genome Biol. 2013;14:R78.
Cheng L, Doecke JD, Sharples RA, Villemagne VL, Fowler CJ, Rembach A, et al. Prognostic serum miRNA biomarkers associated with Alzheimer’s disease shows concordance with neuropsychological and neuroimaging assessment. Mol Psychiatry. 2014. doi:10.1038/mp.2014.127.
Sheinerman KS, Tsivinsky VG, Abdullah L, Crawford F, Samuil R. Plasma microRNA biomarkers for detection of mild cognitive impairment: biomarker validation study. Aging. 2013;5:925–38.
Lunnon K, Sattlecker M, Furney SJ, Coppola G, Simmons A, Proitsi P, et al. A blood gene expression marker of early Alzheimer’s disease. J Alzheimers Dis. 2013;33:737–53.
Timmons JA, Knudsen S, Rankinen T, Koch LG, Sarzynski M, Jensen T, et al. Using molecular classification to predict gains in maximal aerobic capacity following endurance exercise training in humans. J Appl Physiol. 2010;108:1487–96.
Parry RM, Jones W, Stokes TH, Phan JH, Moffitt RA, Fang H, et al. k-Nearest neighbor models for microarray gene expression analysis and clinical outcome prediction. Pharmacogenomics J. 2010;10:292–309.
Gallagher IJ, Scheele C, Keller P, Nielsen AR, Remenyi J, Fischer CP, et al. Integration of microRNA changes in vivo identifies novel molecular features of muscle insulin resistance in type 2 diabetes. Genome Med. 2010;2:9.
Haury A-C, Gestraud P, Vert J-P. The influence of feature selection methods on accuracy, stability and interpretability of molecular signatures. PLoS One. 2011;6:e28210.
Gallagher IJ, Stephens NA, MacDonald AJ, Skipworth RJ, Husi H, Greig C, et al. Suppression of skeletal muscle turnover in cancer cachexia: evidence from the transcriptome in sequential human muscle biopsies. Clin Cancer Res. 2012;18:2817–27.
Thalacker-Mercer AE, Fleet JC, Craig BA, Campbell WW. The skeletal muscle transcript profile reflects accommodative responses to inadequate protein intake in younger and older males. J Nutr Biochem. 2010;21:1076–82.
Raue U, Trappe TA, Estrem ST, Qian H-R, Helvering LM, Smith RC, et al. Transcriptome signature of resistance exercise adaptations: mixed muscle and fiber type specific profiles in young and old adults. J Appl Physiol. 2012;112:1625–36.
Liu D, Sartor MA, Nader GA, Pistilli EE, Tanton L, Lilly C, et al. Microarray analysis reveals novel features of the muscle aging process in men and women. J Gerontol A Biol Sci Med Sci. 2013;1–10.
Berchtold NC, Cribbs DH, Coleman PD, Rogers J, Head E, Kim R, et al. Gene expression changes in the course of normal brain aging are sexually dimorphic. Proc Natl Acad Sci U S A. 2008;105:15605–10.
Dunder K, Lind L, Zethelius B, Berglund L, Lithell H. Evaluation of a scoring scheme, including proinsulin and the apolipoprotein B/apolipoprotein A1 ratio, for the risk of acute coronary events in middle-aged men: Uppsala Longitudinal Study of Adult Men (ULSAM). Am Heart J. 2004;148:596–601.
Ramasamy A, Trabzuni D, Guelfi S, Varghese V, Smith C, Walker R, et al. Genetic variability in the regulation of gene expression in ten regions of the human brain. Nat Neurosci. 2014;17:1418–28.
Lunnon K, Ibrahim Z, Proitsi P, Lourdusamy A, Newhouse S, Sattlecker M, et al. Mitochondrial dysfunction and immune activation are detectable in early Alzheimer’s disease blood. J Alzheimers Dis. 2012;30:685–710.
Snyder HM, Carrillo MC, Grodstein F, Henriksen K, Jeromin A, Lovestone S, et al. Developing novel blood-based biomarkers for Alzheimer’s disease. Alzheimers Dement. 2014;10:109–14.
De Preter K, Barriot R, Speleman F, Vandesompele J, Moreau Y. Positional gene enrichment analysis of gene sets for high-resolution identification of overrepresented chromosomal regions. Nucleic Acids Res. 2008;36:e43.
Garasto S, Rose G, Derango F, Berardelli M, Corsonello A, Feraco E, et al. The study of APOA1, APOC3 and APOA4 variability in healthy ageing people reveals another paradox in the oldest old subjects. Ann Hum Genet. 2003;67:54–62.
Feitosa MF, Wojczynski MK, Straka R, Kammerer CM, Lee JH, Kraja AT, et al. Genetic analysis of long-lived families reveals novel variants influencing high density-lipoprotein cholesterol. Front Genet. 2014;5:159.
Talseth-Palmer BA, Wijnen JT, Brenne IS, Jagmohan-Changur S, Barker D, Ashton KA, et al. Combined analysis of three Lynch syndrome cohorts confirms the modifying effects of 8q23.3 and 11q23.1 in MLH1 mutation carriers. Int J Cancer. 2013;132:1556–64.
Lubbe SJ, Di Bernardo MC, Broderick P, Chandler I, Houlston RS. Comprehensive evaluation of the impact of 14 genetic variants on colorectal cancer phenotype and risk. Am J Epidemiol. 2012;175:1–10.
Audenet F, Cancel-Tassin G, Bigot P, Audouin M, Gaffory C, Ondet V, et al. Germline genetic variations at 11q13 and 12p11 locus modulate age at onset for renal cell carcinoma. J Urol. 2014;191:487–92.
Lange EM, Salinas CA, Zuhlke KA, Ray AM, Wang Y, Lu Y, et al. Early onset prostate cancer has a significant genetic component. Prostate. 2012;72:147–56.
Jin G, Lu L, Cooney KA, Ray AM, Zuhlke KA, Lange EM, et al. Validation of prostate cancer risk-related loci identified from genome-wide association studies using family-based association analysis: evidence from the International Consortium for Prostate Cancer Genetics (ICPCG). Hum Genet. 2012;131:1095–103.
Gierman HJ, Fortney K, Roach JC, Coles NS, Li H, Glusman G, et al. Whole-genome sequencing of the world’s oldest people. PLoS One. 2014;9:e112430.
Kenyon CJ. The genetics of ageing. Nature. 2010;464:504–12.
Beltran Valls MR, Wilkinson DJ, Narici MV, Smith K, Phillips BE, Caporossi D, et al. Protein carbonylation and heat shock proteins in human skeletal muscle: relationships to age and sarcopenia. J Gerontol A Biol Sci Med Sci. 2014;70:174–81.
Raz N, Lindenberger U, Rodrigue KM, Kennedy KM, Head D, Williamson A, et al. Regional brain changes in aging healthy adults: General trends, individual differences and modifiers. Cereb Cortex. 2005;15:1676–89.
Puca A, Daly MJ, Brewster SJ, Matise TC, Barrett J, Shea-Drinkwater M, et al. A genome-wide scan for linkage to human exceptional longevity identifies a locus on chromosome 4. Proc Natl Acad Sci U S A. 2001;98:10505–8.
Sebastiani P, Solovieff N, Dewan AT, Walsh KM, Puca A, Hartley SW, et al. Genetic signatures of exceptional longevity in humans. PLoS One. 2012;7:e29848.
Sabia S, Singh-Manoux A, Hagger-Johnson G, Cambois E, Brunner EJ, Kivimaki M. Influence of individual and combined healthy behaviours on successful aging. CMAJ. 2012;184:1985–92.
Danielsson M, Talbäck M. Public health: an overview: Health in Sweden: The National Public Health Report 2012. Chapter 1. Scand J Public Health. 2012;40:6.
Zethelius B, Berglund L, Sundström J, Ingelsson E, Basu S, Larsson A, et al. Use of multiple biomarkers to improve the prediction of death from cardiovascular causes. N Engl J Med. 2008;358:2107–16.
Raichlen DA, Alexander GE. Exercise, APOE genotype, and the evolution of the human lifespan. Trends Neurosci. 2014;37:247–55.
Romijn MDM, van Marum RJ, Emmelot-Vonk MH, Verhaar HJJ, Koek HL. Mild chronic kidney disease is associated with cognitive function in patients presenting at a memory clinic. Int J Geriatr Psychiatry. 2014;30:758–65.
Norton S, Matthews FE, Barnes DE, Yaffe K, Brayne C. Potential for primary prevention of Alzheimer’s disease: an analysis of population-based data. Lancet Neurol. 2014;13:788–94.
Thambisetty M, Jeffrey Metter E, Yang A, Dolan H, Marano C, Zonderman AB, et al. Glucose intolerance, insulin resistance, and pathological features of Alzheimer disease in the Baltimore Longitudinal Study of Aging. JAMA Neurol. 2013;70:1167–72.
Akomolafe A, Beiser A, Meigs JB, Au R, Green RC, Farrer LA, et al. Diabetes mellitus and risk of developing Alzheimer disease: results from the Framingham Study. Arch Neurol. 2006;63:1551–5.
Liu R, Sui X, Laditka JN, Church TS, Colabianchi N, Hussey J, et al. Cardiorespiratory fitness as a predictor of dementia mortality in men and women. Med Sci Sports Exerc. 2012;44:253–9.
Bouchard C, Rankinen T, Timmons JA. Genomics and genetics in the biology of adaptation to exercise. Compr Physiol. 2011;1:1603–48.
Hampel H, Lista S, Teipel SJ, Garaci F, Nisticò R, Blennow K, et al. Perspective on future role of biological markers in clinical therapy trials of Alzheimer’s disease: a long-range point of view beyond 2020. Biochem Pharmacol. 2014;88:426–49.
Naj AC, Jun G, Reitz C, Kunkle BW, Perry W, Park YS, et al. Effects of multiple genetic loci on age at onset in late-onset Alzheimer disease a genome-wide association study. JAMA Neurol. 2014;71:1394–404.
Lotz M, Martel-Pelletier J, Christiansen C, Brandi M-L, Bruyère O, Chapurlat R, et al. Value of biomarkers in osteoarthritis: current status and perspectives. Ann Rheum Dis. 2013;72:1756–63.
The international Alzheimer’s disease research portfolio. http://iadrp.nia.nih.gov
Exalto LG, Quesenberry CP, Barnes D, Kivipelto M, Biessels GJ, Whitmer RA. Midlife risk score for the prediction of dementia four decades later. Alzheimers Dement. 2014;10:562–70.
Exalto LG, Biessels GJ, Karter AJ, Huang ES, Katon WJ, Minkoff JR, et al. Risk score for prediction of 10 year dementia risk in individuals with type 2 diabetes: a cohort study. Lancet Diabetes Endocrinol. 2013;1:183–90.
Glorioso C, Oh S, Douillard GG, Sibille E. Brain molecular aging, promotion of neurological disease and modulation by sirtuin 5 longevity gene polymorphism. Neurobiol Dis. 2011;41:279–90.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80.
Gautier L, Cope L, Bolstad BM, Irizarry RA. affy--analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004;20:307–15.
Leek JT, Scharpf RB, Bravo HC, Simcha D, Langmead B, Johnson WE, et al. Tackling the widespread and critical impact of batch effects in high-throughput data. Nat Rev Genet. 2010;11:733–9.
Leek JT, Storey JD. Capturing heterogeneity in gene expression studies by surrogate variable analysis. PLoS Genet. 2007;3:1724–35.
Myers J, Prakash M, Froelicher V, Do D, Partington S, Atwood JE. Exercise capacity and mortality among men referred for exercise testing. N Engl J Med. 2002;346:793–801.
Smyth GK. Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:3.
Sing T, Sander O, Beerenwinkel N, Lengauer T. ROCR: visualizing classifier performance in R. Bioinformatics. 2005;21:3940–1.
Shao L, Fan X, Cheng N, Wu L, Cheng Y. Determination of minimum training sample size for microarray-based cancer outcome prediction-an empirical assessment. PLoS One. 2013;8:e68579.
Slentz CA, Bateman LA, Willis LH, Shields AT, Tanner CJ, Piner LW, et al. The Effects of Aerobic versus Resistance Training on Visceral and Liver Fat Stores, Liver Enzymes and HOMA from STRRIDE AT/RT: A Randomized Trial. Am J Physiol Endocrinol Metab. 2011;301:E1033–9.
Timmons JA, Larsson O, Jansson E, Fischer H, Gustafsson T, Greenhaff PL, et al. Human muscle gene expression responses to endurance training provide a novel perspective on Duchenne muscular dystrophy. FASEB J. 2005;19:750–60.
Hagström E, Ingelsson E, Sundström J, Hellman P, Larsson TE, Berglund L, et al. Plasma parathyroid hormone and risk of congestive heart failure in the community. Eur J Heart Fail. 2010;12:1186–92.
Lovestone S, Francis P, Kloszewska I, Mecocci P, Simmons A, Soininen H, et al. AddNeuroMed--the European collaboration for the discovery of novel biomarkers for Alzheimer’s disease. Ann N Y Acad Sci. 2009;1180:36–46.
Tabassum R, Cunningham L, Stephens EH, Sturdivant K, Martin GS, Brigham KL, et al. A longitudinal study of health improvement in the Atlanta CHDWB Wellness Cohort. J Pers Med. 2014;4:489–507.
Sinnaeve PR, Donahue MP, Grass P, Seo D, Vonderscher J, Chibout SD, et al. Gene expression patterns in peripheral blood correlate with the extent of coronary artery disease. PLoS One. 2009;4:e7037.
Ingenuity pathway analysis. http://www.ingenuity.com
Huang X, Sjögren P, Arnlöv J, Cederholm T, Lind L, Stenvinkel P, et al. Serum fatty acid patterns, insulin sensitivity and the metabolic syndrome in individuals with chronic kidney disease. J Intern Med. 2013;275:71–83.
We would like to thank Iwona Kloszewska, Patrizia Mecocci, Hilkka Soininen, Magda Tsolaki, Bruno Vellas, and Simon Lovestone for playing a major role in developing the AD cohort resource. We would like to thank Rufina Leung, Megan Pritchard, Belinda Martin, Kathryn Lord, Nicola Dunlop, and Catherine Tunnard for sample collection and neuropsychiatric assessments at the London site. We would like to thank Nabila Kazmi for providing provisional coding input into the project under the supervision of IJG and JAT. We thank Dr Steen Knudsen (MPI) and Dr Helge and Dr Dela (Panum Institute) for discussions around the project validation studies in 2012. Parts of this work were supported by InnoMed, (Innovative Medicines in Europe) an Integrated Project funded by the European Union of the Sixth Framework program priority FP6-2004-LIFESCIHEALTH-5, Alzheimer’s Research UK, The John and Lucille van Geest Foundation, the National Institute of Health Research Biomedical Research Centre for Mental Health at the South London and Maudsley NHS Foundation Trust and [Institute of Psychiatry] Kings College London, Psychiatry Research Trust, Institute of Psychiatry and the Rosetrees Trust. AH receives support from the Innovative Medicines Initiative Joint Undertaking under EMIF grant agreement number 115372, resources of which are composed of financial contribution from the European Union’s Seventh Framework Programme (FP7/2007-2013) and EFPIA companies in kind contribution. The clinical validation studies were supported by a grant on ageing to JT from the Medical Research Council, UK (G1100015/1). Aspects of this study were also supported by the Wallenberg Foundation (TG), Karolinska Institutet (JT), Swedish Medical Research Council Grant D0328602 (TG) and an Affymetrix Translational Medicine Grant (JT). ER is partly supported by a grant from Swedish Society for Medical Research (SSMF). BEK (STRRIDE Study) was supported by the NIH.
A patent application has been filed on the age diagnostic and JT, TG, BP, ER, TJ, and IJG are associated with and/or are shareholders in XRGenomics Ltd. No other competing interests exist.
The original idea and data analysis approach for the healthy ageing diagnostic was developed by JAT. The broader project objectives were discussed and refined with PJA, TG, WEK, LJCvL, AH, and RH. Laboratory studies were carried out by BEP, TJ, LL, TC, WEK, TG, AH, KL, AK, and ER. Informatics analysis strategy was developed by JAT, IJG, and SS and informatics analysis carried out by SS, IJG, ER, and JAT. HC, JAT and PJA carried out PubMed analysis. The manuscript was drafted by JAT, SS, AH, and WEK and edited by AH, RH, JAT, SS, IJG, ER, HC, BEP, LJCvL, LL, WEK, TC, PJA, and TG. All authors read and approved the final manuscript.
Thomas Gustafsson Angela Hodges and James A Timmons are Joint Senior Authors.
An Excel spreadsheet containing data related to our study with six tabs: 1) analysis of the top 150 genes from age prototype classifier in PUBMED; 2) the top 670 genes from the first stage of the project; 3) phenodata for the training data set and validation data sets in our study; 4) a list of prior markers in the literature for Alzheimer’s disease; 5) positional gene enrichment analysis; 6) sample information for the BrainEac study. (XLSX 221 kb)
Table S1. Clinical regression data from the ULSAM cohort. (DOCX 127 kb)
Figure S1. A cumulative ranking metric of the healthy ageing metric was prognostic for mortality over a 20-year follow-up period. One-hundred and eight subjects provided a healthy tissue biopsy in 1992 that was suitable for RNA profiling and the fully annotated mortality data, covering 2009–2011, was retrieved from the Swedish national health registry. a The rank score for healthy ageing gene expression was calculated from the top 150 genes of the healthy ageing prototype classifier (n = 108, male subjects all ~70 years of age). Logistic regression analysis performed using the cumulative ranking metric of the top 150 genes from original prototype was prognostic for mortality. It showed that those subjects with the lowest median healthy ageing gene score had a much higher probability of death during the 20-year follow-up period (p = 0.0295). In contrast, members of the inflammatory response (GO:0006954) and mitochondrion (GO:0005739) gene ontology families - selected from ENSEMBL (BioMart) - showed no significant relationship with health during the 20-year follow-up period (p = 0.34 and p = 0.17). b The rank score for healthy ageing gene expression was calculated from the top 150 genes of the healthy ageing prototype classifier (n = 108, male subjects all ~70 years of age) and Kaplan–Meier plots were used to illustrate the temporal pattern of survival. Gene score was divided into quartiles and the plot was produced using the plot-survfit function in the R survival package. The plot allows us to compare overall survival rates between the four quartiles for gene score. The third and fourth quartiles differed from the first quartile (p < 0.04). (PDF 46 kb)
Figure S2. Diabetes and vascular disease plots. The healthy aging signature activation was studied in blood samples from two independent large case–control studies of diabetes and vascular disease. Applying a Wilcoxon rank sum test, neither diabetes nor vascular disease was related to the healthy ageing gene score. This is consistent with our original hypothesis, and methods, that the healthy ageing gene score is not related to lifestyle factors and it is also consistent with the results observed in the ULSAM cohort (Fig. 3). a The diabetes data (94 controls versus 50 cases, group mean age = 66 years) originates from Tabassum et al.  (using Illumina Human HT.12.V4 arrays). b The vascular disease data (112 controls and 110 cases, group age = 53.3 years) originates from Sinnaeve et al.  (using Affymetrix HG-U133A arrays). (PDF 44 kb)
About this article
Cite this article
Sood, S., Gallagher, I.J., Lunnon, K. et al. A novel multi-tissue RNA diagnostic of healthy ageing relates to cognitive health status. Genome Biol 16, 185 (2015). https://doi.org/10.1186/s13059-015-0750-x