PT  - JOURNAL ARTICLE
AU  - Alan Le Goallec
AU  - Sasha Collin
AU  - Samuel Diai
AU  - Jean-Baptiste Prost
AU  - M’Hamed Jabri
AU  - Théo Vincent
AU  - Chirag J. Patel
TI  - Analyzing the multidimensionality of biological aging with the tools of deep learning across diverse image-based and physiological indicators yields robust age predictors
AID  - 10.1101/2021.04.25.21255767
DP  - 2021 Jan 01
TA  - medRxiv
PG  - 2021.04.25.21255767
4099  - http://medrxiv.org/content/early/2021/04/26/2021.04.25.21255767.short
4100  - http://medrxiv.org/content/early/2021/04/26/2021.04.25.21255767.full
AB  - 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 StatementThe authors have declared no competing interest.Funding StatementNIEHS 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 DeclarationsI confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.YesThe details of the IRB/oversight body that provided approval or exemption for the research described are given below:IRB#: Harvard University IRB16-2145All necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.YesI 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.YesWe 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 UKBUK BiobankNHANESNational Health and Nutrition Examination SurveyECGelectrocardiogramEEGBrain electroencephalographyMRImagnetic resonance imageOCToptical coherence tomographyDXAdual-energy X-ray absorptiometryPWApulse wave analysisCAchronological ageBAbiological ageAAaccelerated agingCPUcentral processing unitGPUgraphics processing unitSLURMSimple Linux Utility for Resource ManagementCVcross-validationOCVouter cross-validationICVinner cross-validationNCVnested cross-validationENelastic netGBMgradient-boosted machineNNneural networkDNNdeep neural networkCNNconvolutional neural network1D CNNone-dimensional convolutional neural network2D CNNconvolutional neural network3D CNNconvolutional neural networkRNNrecurrent neural networkLSTMlong short-term memory neural networkGRUgated recurrent units neural networkReLUrectified linear unitSELUscaled exponential linear unitCAMclass activation mappingRAMregression activation mappingGrad-CAMgradient-weighted class activation mappingGrad-RAMgradient-weighted regression activation mappingTPETree-structured Parzen Estimator ApproachR2R-SquaredMSEmean squared errorRMSEroot mean squared errorMAEmean absolute errorGWASgenome wide association studySNPsingle-nucleotide polymorphismLDlinkage disequilibriumQ-Q plotquantile-quantile plotBWASbiomarkers wide association studyEWASenvironmental wide association studyPEWASphenotypic and environmental wide association studyXWASX-wide association studyCVDcardiovascular disease