Plants respond to stresses like drought and heat through complex gene regulatory networks (GRNs). To improve resilience, understanding these is crucial, but large-scale GRNs (>100 genes) are difficult to model using ordinary differential equations (ODEs) due to the high number of parameters that have to be estimated. Here we solve this problem by introducing BADDADAN, which uses machine learning to identify gene modules—groups of co-expressed and/or co-regulated genes—and constructs an ODE model that predicts gene module dynamics under stress. By integrating time-series gene expression data with prior co-expression data it finds modules that are both coherent and interpretable. We demonstrate BADDADAN on heat and drought datasets of A. thaliana, modelling over 1,000 genes, recovering known mechanistic insights, and proposing new hypotheses. By combining machine learning with mechanistic modelling, BADDADAN deepens our understanding of stress-related GRNs in plants and potentially other organisms.

Federated learning is a technique that enables the use of distributed datasets for machine learning purposes without requiring data to be pooled, thereby better preserving privacy and ownership of the data. While supervised FL research has grown substantially over the last years, unsupervised FL methods remain scarce. This work introduces an algorithm which implements K-means clustering in a federated manner, addressing the challenges of varying number of clusters between centers, as well as convergence on less separable datasets.

Background
Many types of dementia have high heritability, which creates opportunities for DNA diagnostics. Clinicians sporadically test for causal genetic variants. However, in addition to causal genetic mutations, an increasing number of both common and rare risk factors are being identified, especially for Alzheimer’s disease (AD). Here, we describe and evaluate diagnostic performance of combining genetic risk factors for AD to assist memory clinic clinicians.

Methods
A retrospective analysis of 998 consecutive patients (mean age 62.1, 40.3% females, 63.3% dementia) was conducted over 2.5 years in a Dutch memory clinic. The patients underwent a complete genetic risk assessment, including whole-exome sequencing and array genotyping. We examined known pathogenic genetic variants for all dementia types and their correlation with clinical diagnoses. We evaluated a combined genetic score (GS) based on all genetic risk factors for AD - namely APOE genotypes, candidate risk rare variants in 11 genes, and a polygenic risk score (PRS) based on 82 common variants. Then, we analyzed the discriminatory characteristics of the GS.

Results
Causal pathogenic variants were rare, present in 3.4% of individuals, but genetic testing would have altered the diagnosis in over half of the carriers. Candidate rare risk variants were more common, identified in 31.6% of patients. Both APOE genotypes and the PRS were independently associated with AD, and gene-specific interaction was found between TREM2 and AD-PRS (β = -1.16, p = 0.015). Patients with a high GS were 7 times more likely receive an AD diagnosis compared to those with a low GS (p = 2.5E-07).

Conclusion
Overall, this study highlights the potential of integrating genetic risk factors into clinical practice to enhance AD diagnosis, though the improvement in diagnostic accuracy was moderate. The findings underscore the importance of genetic testing in diagnosis while also recognizing its limitations.

Federated learning is an upcoming machine learning paradigm which allows data from multiple sources to be used for training of classifiers without the data leaving the source it originally resides. This can be highly valuable for use cases such as medical research, where gathering data at a central location can be quite complicated due to privacy and legal concerns of the data. In such cases, federated learning has the potential to vastly speed up the research cycle. Although federated and central learning have been compared from a theoretical perspective, an extensive experimental comparison of performances and learning behavior still lacks. We have performed a comprehensive experimental comparison between federated and centralized learning. We evaluated various classifiers on various datasets exploring influences of different sample distributions as well as different class distributions across the clients. The results show similar performances under a wide variety of settings between the federated and central learning strategies. Federated learning is able to deal with various imbalances in the data distributions. It is sensitive to batch effects between different datasets when they coincide with location, similar to central learning, but this setting might go unobserved more easily. Federated learning seems to be robust to various challenges such as skewed data distributions, high data dimensionality, multiclass problems, and complex models. Taken together, the insights from our comparison gives much promise for applying federated learning as an alternative to sharing data. Code for reproducing the results in this work can be found at: https://github.com/swiergarst/FLComparison.

Smartwatches enable longitudinal and continuous data acquisition. This has the potential to remotely monitor (changes) of the health of users. However, differences among subjects (inter-subject variability) limit a model to generalize to unseen subjects. This study focused on binary classification tasks using heart rate and step counter from smartwatches, including night/day and inactive/active classification, as well as sleep and SpO2-related (oxygen saturation) tasks. To address inter-subject variability, we explored different transforming and normalization regimes for time series including per-subject and population-based strategies. We propose a modified factorized autoencoder, which separates the data into two latent spaces capturing class-specific and subject-specific information. Our proposed generalized factorized autoencoder and triplet factorized autoencoder improved classification accuracy over the baseline from 74.8 (± 10.5) to 83.1 (± 5.1) and 83.4 (± 5.3), respectively, for night/day classification, gains for inactive/active classification were modest, improving from 84.3 (± 9.4) to 86.9 (± 4.4) and 86.6 (± 4.3), respectively. Our study highlights challenges of handling inter-subject variability in smartwatch data and how factorization models can be used to enable more robust and personalized health monitoring solutions for diverse populations.

Background: Despite several randomised controlled trials (RCTs) on the use of adjuvant treatment with corticosteroids in patients with community-acquired pneumonia (CAP), the effect of this intervention on mortality remains controversial. We aimed to evaluate heterogeneity of treatment effect (HTE) of adjuvant treatment with corticosteroids on 30-day mortality in patients with CAP. Methods: In this individual patient data meta-analysis, we included RCTs published before July 1, 2024, comparing adjuvant treatment with corticosteroids versus placebo in patients hospitalised with CAP. The primary endpoint was 30-day all-cause mortality, collected across all trials, and analyses followed the intention-to-treat principle. We analysed HTE using risk and effect modelling. For risk modelling, patients were classified as having less severe or severe CAP based on the pneumonia severity index (PSI), comparing PSI class I–III versus class IV–V. For effect modelling, we trained a corticosteroid-effect model on six trials and externally validated it using data from two trials, received after model preregistration. This model classified patients into two groups: no predicted benefit and predicted benefit from adjuvant treatment with corticosteroids. The literature search was registered on PROSPERO, CRD42022380746. Findings: We included eight RCTs with 3224 patients. Across all eight trials, 246 (7·6%) patients died within 30 days (106 [6·6%] of 1618 in the corticosteroid group vs 140 [8·7%] of 1606 in the placebo group; odds ratio [OR] 0·72 [95% CI 0·56–0·94], p=0·017). The corticosteroid-effect model, which selected C-reactive protein (CRP), showed significant HTE during external validation in the two most recent trials. In these trials, 154 (11·4%) of 1355 patients died within 30 days (88 [13·1%] of 671 in the placebo group vs 66 [9·6%] of 684 in the corticosteroid group; OR 0·71 [95% CI 0·50–0·99], p=0·044). Among patients predicted to have no benefit (CRP ≤204 mg/L, n=725), no significant effect was observed (OR 0·98 [95% CI 0·63–1·50]), whereas for those with predicted benefit (CRP >204 mg/L, n=630), 39 (13·0%) of 301 patients died in the placebo group compared with 20 (6·1%) of 329 in the corticosteroid group (0·43 [0·25–0·76], p_interaction=0·026). No significant HTE was found between less severe CAP (PSI class I–III, n=229) and severe CAP (PSI class IV–V, n=1126). Corticosteroid therapy significantly increased hyperglycaemia risk (44 [12·8%] of 344 in the placebo group vs 84 [24·8%] of 339 in the corticosteroid group; OR 2·50 [95% CI 1·63–3·83], p<0·0001) and hospital re-admission risk (30 [3·7%] of 814 in the placebo group vs 57 [7·0%] of 819 in the corticosteroid group; 1·95 [1·24–3·07], p=0·0038). Interpretation: Overall, adjuvant therapy with corticosteroids significantly reduces 30-day mortality in patients hospitalised with CAP. The treatment effect varied significantly among subgroups based on CRP concentrations, with a substantial mortality reduction observed only in patients with high baseline CRP. Funding: None.

Introduction
The field of forensic DNA analysis has undergone rapid advancements in recent decades. The integration of massively parallel sequencing (MPS) has notably expanded the forensic toolkit, moving beyond identity matching to predicting phenotypic traits and biogeographical ancestry. This shift is of particular significance in cases where conventional DNA profiling fails to identify a single suspect. Supplementing forensic analyses with estimated biological age may be valuable but involves a complex and time-consuming DNA methylation analysis. This study explores and validates the performance of a comprehensive forensic third-generation sequencing assay utilizing Oxford Nanopore Technologies (ONT) in an adaptive and direct sequencing approach. We incorporated the most widely used forensic markers, i.e., STRs, SNPs, InDels, mitochondrial DNA (mtDNA), and two methylation-based clock classifiers, thereby combining forensic genetic and epigenetic analysis in one single workflow.

Methods and results
In our investigation, DNA from six anonymous individuals was sequenced using the ONT standard adaptive direct sequencing approach, reaching a mean percentage of on-target reads ranging from 6.6 % to 7.7 % per sample. ONT data was compared to standard MPS data and Illumina EPIC DNA methylation profiles. Basecalling employed recommended ONT software packages. TREAT was used for ONT-based analysis of autosomal and Y-chromosome STRs, achieving 90–92 % correct calls depending on allelic read depth thresholds. InDel analyses for two lower-quality samples proved challenging due to inadequate read depth, while the remaining four samples significantly contributed to the observed percentage markers (60.9 %) and correct calls (97.8 %). SNP analysis achieved a 98 % call rate, with only two mismatches and two missed alleles. ONT-generated DNA methylation data demonstrated Pearson’s correlation coefficients with EPIC data ranging from 0.67 to 0.97 for Horvath’s clock. Additional age-associated markers exhibited Pearson’s correlation coefficients with chronological age between 0.14 (ELOVL2) and 0.96 (FHL2) at read depths of <30 and <20, respectively. Despite excluding mtDNA from our targeted sequencing approach, adaptive proof-reading fragments covered the complete mtDNA with an average read depth of 21–72, showing 100 % concordance with reference data.

Discussion
Our exploratory study using ONT adaptive sequencing for conventional forensic and age associated DNA methylation markers showed high sequencing accuracy for a significant number of markers, showcasing ONT as a promising (epi)genetic forensic method. Future studies must address three critical aspects: determining clear quantity and quality measures and detection thresholds for accuracy, optimizing input DNA quantity for forensic casework expectations, and addressing ethical considerations associated with phenotype and ancestry analysis to prevent ethnic biases.

While deep learning (DL) is used in patients’ outcome predictions, the insufficiency of patient samples limits the accuracy. In this study, we investigated how transfer learning (TL) alleviates the small sample size problem. A 2-step TL framework was constructed for a difficult task: predicting the response of the drug temozolomide (TMZ) in glioblastoma (GBM) cell cultures. The GBM is aggressive, and most patients do not benefit from the only approved chemotherapeutic agent TMZ. O6-methylguanine-DNA methyltransferase (MGMT) promoter methylation status is the only biomarker for TMZ responsiveness but has shown limited predictive power. The 2-step TL framework was built on 3 datasets: (1) the subset of the Genomics of Drug Sensitivity in Cancer (GDSC) dataset, including miscellaneous cell cultures treated by TMZ, cyclophosphamide, bortezomib, and oxaliplatin, as the source dataset; (2) the Human Glioblastoma Cell Culture (HGCC) dataset, for fine-tuning; and (3) a small target dataset GSE232173, for validation. The latter two included specifically TMZ-treated GBM cell cultures. The DL models were pretrained on the cell cultures treated by each of the 4 drugs from GDSC, respectively. Then, the DL models were refined on HGCC, where the best source drug was identified. Finally, the DL model was validated on GSE232173. Using 2-step TL with pretraining on oxaliplatin was not only superior to those without TL and with 1-step TL but also better than 3 benchmark methods, including MGMT. The oxaliplatin-based TL improved the performance probably by increasing the weights of cell cycle-related genes, which relates to the TMZ response processes. Our findings support the potential of oxaliplatin being an alternative therapy for patients with GBM and TL facilitating drug repurposing research. We recommend that following our methodology, using mixed cancers and a related drug as the source and then fine-tuning the model with the target cancer and the target drug will enhance drug response prediction.

We appreciate the opportunity to further clarify our findings in response to the insightful comments from Shota Yamamoto and colleagues and Luis Felipe Reyes and Ignacio Martin-Loeches regarding our recent community-acquired pneumonia (CAP) study. [...]

Rheumatoid arthritis (RA) is a heterogeneous disease with variable symptoms, prognosis, and treatment response, necessitating refined patient classification. We applied multimodal deep learning and clustering to identify distinct RA phenotypes using baseline clinical data from 1,387 patients in the Leiden Rheumatology clinic. Four Joint Involvement Patterns (JIP) emerged: foot-predominant arthritis, seropositive oligoarticular disease, seronegative hand arthritis, and polyarthritis. Findings were validated in clinical trial data (n = 307) and an independent secondary care cohort (n = 515). Clusters showed high stability and significant differences in remission rates (P = 0.007) and methotrexate failure (P < 0.001). JIP-hand patients had superior outcomes (particularly in ACPA-positive patients) versus JIP-foot (HR:0.37, P < 0.001) and JIP-poly (HR:0.33, P = 0.005), independent of baseline disease activity and clinical markers. Synovial histology analysis (n = 194) revealed distinct inflammatory patterns across clusters, hinting at different underlying biological mechanisms. These validated RA phenotypes based on joint involvement patterns may enable targeted research into disease mechanisms and personalized treatment strategies.

A polygenic score (PGS) for Alzheimer’s disease (AD) was derived recently from data on genome-wide significant loci in European ancestry populations. We applied this PGS to populations in 17 European countries and observed a consistent association with the AD risk, age at onset and cerebrospinal fluid levels of AD biomarkers, independently of apolipoprotein E locus (APOE). This PGS was also associated with the AD risk in many other populations of diverse ancestries. A cross-ancestry polygenic risk score improved the association with the AD risk in most of the multiancestry populations tested when the APOE region was included. Finally, we found that the PGS/polygenic risk score captured AD-specific information because the association weakened as the diagnosis was broadened. In conclusion, a simple PGS captures the AD-specific genetic information that is common to populations of different ancestries, although studies of more diverse populations are still needed to better characterize the genetics of AD.

Background and ObjectivesIdentifying genetic causes of dementia in patients visiting memory clinics is important for patient care and family planning. Traditional clinical selection criteria for genetic testing may miss carriers of pathogenic variants in dementia-related genes. This study aimed identify how many carriers we are missing and to optimize criteria for selecting patients for genetic counseling in memory clinics.MethodsIn this clinical cohort study, we retrospectively genetically tested patients during 2.5 years (2010-2012) visiting the Alzheimer Center Amsterdam, a specialized memory clinic. Genetic tests consisted of a 54-gene dementia panel, focusing on Class IV/V variants per American College of Medical Genetics and Genomics guidelines, including APP duplications and the C9ORF72 repeat expansion. We determined the prevalence of pathogenic variants and propose new eligibility criteria for genetic testing in memory clinics. The eligibility criteria were prospectively applied for 1 year (2021-2022), and results were compared with the retrospective cohort.ResultsGenetic tests were retrospectively performed in in 1,022 of 1,138 patients (90%) who consecutively visited the memory clinic. Among these, 1,022 patients analyzed (mean age 62.1 ± 8.9 years; 40.4% were female), 34 pathogenic variant carriers were identified (3.3%), with 24 being symptomatic. Previous clinical criteria would have identified only 15 carriers (44% of all carriers, 65% of symptomatic carriers). The proposed criteria increased identification to 22 carriers (62.5% of all carriers, 91% of symptomatic carriers). In the prospective cohort, 148 (28.7%) of 515 patients were eligible for testing under the new criteria. Of the 90 eligible patients who consented to testing, 13 pathogenic carriers were identified, representing a 73% increase compared with the previous criteria.DiscussionWe found that patients who visit a memory clinic and carry a pathogenic genetic variant are often not eligible for genetic testing. The proposed new criteria improve the identification of patients with a genetic cause for their cognitive complaints. In systems without practical or financial barriers to genetic testing, the new criteria can enhance personalized care. In other countries where the health care systems differs and in other genetic ancestry groups, the performance of the criteria may be different.

Background
Switching from controlled to assisted ventilation is crucial in the trajectory of intensive care unit (ICU) stay, but no guidelines exist. We described current practices, analyzed patient characteristics associated with switch success or failure, and explored the feasibility to predict switch failure.

Methods
In this retrospective study, we obtained highly granular longitudinal ICU data sets from three medical centers, covering demographics, severity scores, vital signs, ventilation, and laboratory parameters. The primary endpoint was switch success, considering a switch attempt to be successful if a patient did not return to controlled ventilation for the next 72 h while alive, and to be failed otherwise. We compared the characteristics of patients with successful vs. failed first switch attempts at ICU admission, immediately before, and 3 h after the attempt. We trained LASSO logistic regression models to predict switch failure.

Results
In 4524/6715 (67%) patients attempting a switch, the first attempt failed. The first switch attempt, regardless of success or failure, was generally made at normalized PaCO2 and pH levels, with PEEP < 10 cmH2O and PaO2/FiO2 indicating mild injury. Despite very similar baseline disease severity, switch failure was associated with significantly worse outcomes, including a 28-day mortality of 27% vs. 16% and median ventilator-free days of 16 vs. 22 (p < 0.001). Failed attempts were initiated significantly earlier than successful ones (median 1.8 vs. 1.3 days, p < 0.001). Before the switch, PaO2/FiO2, if measured at PEEP > 10 cmH2O, and respiratory system compliance was lower in patients with switch failure (median 185 vs. 205 mmHg, p < 0.001; 39 vs. 41 mL/cmH2O, P = 0.001), and post-switch, patients with switch failure experienced greater deterioration in gas exchange and minimal improvement in ventilatory parameters post-switch. Contrary to our hypotheses, patient characteristics for failed vs. successful switches were surprisingly similar, resulting in prediction models with limited discriminative performance.

Conclusions
Approximately two-thirds of attempts to switch patients to assisted ventilation fail, which are associated with significantly worse clinical outcomes, despite similar baseline disease severity. Contrary to our hypotheses, patients with successful and failed attempts showed similar characteristics, making switch failure difficult to predict. These findings underscore the importance of preventing switch failures and, given the retrospective nature of this study, highlight the need for prospective studies to better understand the reasons for switch failure and when spontaneous breathing can be safely initiated.

Inflammatory bowel disease (IBD) chronicity results from memory T helper cell (Tmem) reactivation. Identifying patient-specific immunotypes is crucial for tailored treatment. We conducted a comprehensive study integrating circulating immune proteins and circulating Tmem, with intestinal tissue histology and mRNA analysis, in therapy-naïve pediatric IBD (Crohn's disease, CD: n = 62; ulcerative colitis, UC: n = 20; age-matched controls n = 43), and after 10–12 weeks’ induction therapy. At diagnosis, plasma protein profiles unveiled two UC and three CD clusters with distinct disease courses. UC patients displayed unchanged circulating Tmem, while CD exhibited increased frequencies of gut-homing ex-Th17, known for high IFN-γ production. UC#2 had elevated Th17/neutrophil-pathway-related proteins and severe disease, with higher endoscopic and histological damage and Th17/neutrophil infiltration. Although both UC#1 and UC#2 responded to therapy, UC#2 required earlier immunomodulation. CD#3 had lower plasma protein concentrations, especially IFN-γ pathway proteins, fewer gut-homing ex-Th17 and clinically milder disease, confirmed by intestinal gene expression. CD#1 and CD#2 had comparably high Th1-related immune profiles, but CD#1 exhibited higher concentrations of proteins previously associated with poorer prognosis. Both CD clusters responded to induction therapy, with similar one-year outcomes. This study highlights feasibility of discriminating patient-specific immunotypes in IBD, advancing our understanding of immune pathogenesis, needed for tailored treatment strategies.

RNA Velocity allows the inference of cellular differentiation trajectories from single-cell RNA sequencing (scRNA-seq) data. It would be highly interesting to study these differentiation dynamics in the spatial context of tissues. Estimating spatial RNA velocities is, however, limited by the inability to spatially capture spliced and unspliced mRNA molecules in high-resolution spatial transcriptomics. We present SIRV, a method to spatially infer RNA velocities at the single-cell resolution by enriching spatial transcriptomics data with the expression of spliced and unspliced mRNA from reference scRNA-seq data. We used SIRV to infer spatial differentiation trajectories in the developing mouse brain, including the differentiation of midbrain-hindbrain boundary cells and marking the forebrain origin of the cortical hem and diencephalon cells. Our results show that SIRV reveals spatial differentiation patterns not identifiable with scRNA-seq data alone. Additionally, we applied SIRV to mouse organogenesis data and obtained robust spatial differentiation trajectories. Finally, we verified the spatial RNA velocities obtained by SIRV using 10x Visium data of the developing chicken heart and MERFISH data from human osteosarcoma cells. Altogether, SIRV allows the inference of spatial RNA velocities at the single-cell resolution to facilitate studying tissue development.

Introduction
Unraveling how Alzheimer's disease (AD) genetic risk is related to neuropathological heterogeneity, and whether this occurs through specific biological pathways, is a key step toward precision medicine.

Methods
We computed pathway-specific genetic risk scores (GRSs) in non-demented individuals and investigated how AD risk variants predict cerebrospinal fluid (CSF) and imaging biomarkers reflecting AD pathology, cardiovascular, white matter integrity, and brain connectivity.

Results
CSF amyloidbeta and phosphorylated tau were related to most GRSs. Inflammatory pathways were associated with cerebrovascular disease, whereas quantitative measures of white matter lesion and microstructure integrity were predicted by clearance and migration pathways. Functional connectivity alterations were related to genetic variants involved in signal transduction and synaptic communication.

Discussion
This study reveals distinct genetic risk profiles in association with specific pathophysiological aspects in predementia stages of AD, unraveling the biological substrates of the heterogeneity of AD-associated endophenotypes and promoting a step forward in disease understanding and development of personalized therapies.

Highlights
- Polygenic risk for Alzheimer's disease encompasses six biological pathways that can be quantified with pathway-specific genetic risk scores, and differentially relate to cerebrospinal fluid and imaging biomarkers.
- Inflammatory pathways are mostly related to cerebrovascular burden.
- White matter health is associated with pathways of clearance and membrane integrity, whereas functional connectivity measures are related to signal transduction and synaptic communication pathways.

Molecular effects of lifestyle interventions are typically studied in a single tissue. Here, we perform a secondary analysis on the sex-specific effects of the Growing Old TOgether trial (GOTO, trial registration number GOT NL3301 (https://onderzoekmetmensen.nl/nl/trial/27183), NL-OMON27183, primary outcomes have been previously reported in ref. 1), a moderate 13-week combined lifestyle intervention on the transcriptomes of postprandial blood, subcutaneous adipose tissue (SAT) and muscle tissue in healthy older adults, the overlap in effect between tissues and their relation to whole-body parameters of metabolic health. The GOTO intervention has virtually no effect on the postprandial blood transcriptome, while the SAT and muscle transcriptomes respond significantly. In SAT, pathways involved in HDL remodeling, O2/CO2 exchange and signaling are overrepresented, while in muscle, collagen and extracellular matrix pathways are significantly overexpressed. Additionally, we find that the effects of the SAT transcriptome closest associates with gains in metabolic health. Lastly, in males, we identify a shared variation between the transcriptomes of the three tissues. We conclude that the GOTO intervention has a significant effect on metabolic and muscle fibre pathways in the SAT and muscle transcriptome, respectively. Aligning the response in the three tissues revealed a blood transcriptome component which may act as an integrated health marker for metabolic intervention effects across tissues.

Cardiovascular disease (CVD) is the most important cause of morbidity and mortality worldwide. Early detection, prevention or even prediction is of pivotal importance to reduce the burden of cardiovascular disease and its associated costs. Low cost, consumer-grade smartwatches have the potential to revolutionize cardiovascular medicine by enabling continuous monitoring of heart rate and activity. When combined with machine learning(ML), the resulting large amounts of time series data hold the potential of detection, or exclusion of CVD. However, analyzing such large datasets is challenging due to the sparse presence of informative segments. Efficient selection of these segments is essential for developing predictive models for clinical deployment. The objective of this paper was to investigate the potential of an acceleration-deceleration curvebased ML model as a novel clinical indicator for the detection of cardiovascular diseases. We used data from the ME-TIME study; 42 participants from which 21 have a cardiovascular disease and 21 are health controls. Data from each subject was normalized to decrease inter-subject variability. A neural network model aggregated predictions per week. We showed that per-subject normalization by the peak value of curves during inactivity, aggregation of model predictions over a week, and using a contrastive loss, resulted in a predictive model with 99 % ± 3 % specificity and 40 % ± 49 % sensitivity on the development set, and 100 % specificity with 67 % ± 47 % sensitivity on the test set. Acceleration-deceleration curves are effective patterns for ruling out the presence of cardiovascular disease, but caution must be taken to properly pre-process the curves and carefully choosing a model that reduces the variability in the extracted curves.

Biological age uses biophysiological information to capture a person’s age-related risk of adverse outcomes. MetaboAge and MetaboHealth are metabolomics-based biomarkers of biological age trained on chronological age and mortality risk, respectively. Lifestyle factors contribute to the extent chronological and biological age differ. The association of lifestyle factors with MetaboAge and MetaboHealth, potential sex differences in these associations, and MetaboAge’s and MetaboHealth’s sensitivity to lifestyle changes have not been studied yet.

Linear regression analyses and mixed-effect models were used to examine the cross-sectional and longitudinal associations of scaled lifestyle factors with scaled MetaboAge and MetaboHealth in 24,332 middle-aged participants from the Doetinchem Cohort Study, Rotterdam Study, and UK Biobank. Random-effect meta-analyses were performed across cohorts. Repeated metabolomics measurements had a ten-year interval in the Doetinchem Cohort Study and a five-year interval in the UK Biobank.

In the first study incorporating longitudinal information on MetaboAge and MetaboHealth, we demonstrate associations between current smoking, sleeping ≥8 hours/day, higher BMI, and larger waist circumference were associated with higher MetaboHealth, the latter two also with higher MetaboAge. Furthermore, adhering to the dietary and physical activity guidelines were inversely associated with MetaboHealth. Lastly, we observed sex differences in the associations between alcohol use and MetaboHealth.