Abstract
It is hypothesized that there are inter-individual differences in biological aging; however, differences in aging among (heart images vs. electrophysiology) and across (e.g., brain vs heart) physiological dimensions have not been systematically evaluated and compared. We analyzed 676,787 samples from 502,211 UK Biobank participants aged 37-82 years with deep learning approaches to build a total of 331 chronological age predictors on different data modalities such as videos (e.g. heart magnetic resonance imaging [MRI]), images (e.g. brain, liver and pancreas MRIs), time-series (e.g. electrocardiograms [ECGs], wrist accelerometer data) and scalar data (e.g. blood biomarkers) to characterize the multiple dimensions of aging. We combined these age predictors into 11 main aging dimensions, 31 subdimensions and 84 sub-subdimensions ensemble models based on specific organ systems. Heart dimension features predict chronological age with a testing root mean squared error (RMSE) and standard error of 2.83±0.04 years and musculoskeletal dimension features predict age with a RMSE of 2.65±0.04 years. We defined “accelerated” agers as participants whose predicted age was greater than their chronological age and computed the correlation between these different definitions of accelerated aging. We found that most aging dimensions are modestly correlated (average correlation=.139±.090) but that dimensions that are biologically related tend to be more positively correlated. For example, we found that heart anatomical (from MRI) accelerated aging and heart electrical (from ECG) accelerated aging are correlated (average Pearson of .249±.005). Overall, most dimensions of aging are complex traits with both genetic and non-genetic correlates. We identified 9,697 SNPs in 3,318 genes associated with accelerated aging and found an average GWAS-based heritability for accelerated aging of 26.1±7.42% (e.g. heart aging: 35.2±1.6%). We used GWAS summary statistics to estimate genetic correlation between aging dimensions and we found that most aging dimensions are genetically not correlated (average correlation=.104±.149). However, on the other hand, specific dimensions were genetically correlated, such as heart anatomical and electrical accelerated aging (Pearson rho .508±.089 correlated [r_g]). Finally, we identified biomarkers, clinical phenotypes, diseases, family history, environmental variables and socioeconomic variables associated with accelerated aging in each aging dimension and computed the correlation between the different aging dimensions in terms of these associations. We found that environmental and socioeconomic variables are similarly associated with accelerated aging across aging dimensions (average correlations of respectively .639±.180 and .607±.309). Dimensions are weakly correlated with each other, highlighting the multidimensionality of the aging process. Our results can be interactively explored on the following website: https://www.multidimensionality-of-aging.net/
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
NIEHS R00 ES023504 NIEHS R21 ES25052. NIAID R01 AI127250 NSF 163870 MassCATS, Massachusetts Life Science Center Sanofi The funders had no role in the study design or drafting of the manuscript(s).
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
IRB#: Harvard University IRB16-2145
All necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.
Yes
Footnotes
↵+ Co-second authors
Contact Information Chirag J Patel 10 Shattuck St. Boston, MA. 02215 (617) 432 1195 chirag_patel{at}hms.harvard.edu
Data Availability
We used the UK Biobank. We will make the phenotypes we generated available upon publication in a peer-reviewed journal.
https://www.multidimensionality-of-aging.net/
https://github.com/alanlegoallec/Multidimensionality_of_Aging
Abbreviations
- UKB
- UK Biobank
- NHANES
- National Health and Nutrition Examination Survey
- ECG
- electrocardiogram
- EEG
- Brain electroencephalography
- MRI
- magnetic resonance image
- OCT
- optical coherence tomography
- DXA
- dual-energy X-ray absorptiometry
- PWA
- pulse wave analysis
- CA
- chronological age
- BA
- biological age
- AA
- accelerated aging
- CPU
- central processing unit
- GPU
- graphics processing unit
- SLURM
- Simple Linux Utility for Resource Management
- CV
- cross-validation
- OCV
- outer cross-validation
- ICV
- inner cross-validation
- NCV
- nested cross-validation
- EN
- elastic net
- GBM
- gradient-boosted machine
- NN
- neural network
- DNN
- deep neural network
- CNN
- convolutional neural network
- 1D CNN
- one-dimensional convolutional neural network
- 2D CNN
- convolutional neural network
- 3D CNN
- convolutional neural network
- RNN
- recurrent neural network
- LSTM
- long short-term memory neural network
- GRU
- gated recurrent units neural network
- ReLU
- rectified linear unit
- SELU
- scaled exponential linear unit
- CAM
- class activation mapping
- RAM
- regression activation mapping
- Grad-CAM
- gradient-weighted class activation mapping
- Grad-RAM
- gradient-weighted regression activation mapping
- TPE
- Tree-structured Parzen Estimator Approach
- R2
- R-Squared
- MSE
- mean squared error
- RMSE
- root mean squared error
- MAE
- mean absolute error
- GWAS
- genome wide association study
- SNP
- single-nucleotide polymorphism
- LD
- linkage disequilibrium
- Q-Q plot
- quantile-quantile plot
- BWAS
- biomarkers wide association study
- EWAS
- environmental wide association study
- PEWAS
- phenotypic and environmental wide association study
- XWAS
- X-wide association study
- CVD
- cardiovascular disease