Advancing protein design is crucial for breakthroughs in medicine and biotechnology. Traditional approaches for protein sequence representation often rely solely on the 20 canonical amino acids, limiting the representation of non-canonical amino acids and residues that undergo post-translational modifications. This work explores discrete diffusion models for generating novel protein sequences using the all-atom chemical representation SELFIES. By encoding the atomic composition of each amino acid in the protein, this approach expands the design possibilities beyond standard sequence representations. Using a modified ByteNet architecture within the discrete diffusion D3PM framework, we evaluate the impact of this all-atom representation on protein quality, diversity, and novelty, compared to conventional amino acid-based models. To this end, we develop a comprehensive assessment pipeline to determine whether generated SELFIES sequences translate into valid proteins containing both canonical and non-canonical amino acids. Additionally, we examine the influence of two noise schedules within the diffusion process—uniform (random replacement of tokens) and absorbing (progressive masking)—on generation performance. While models trained on the all-atom representation struggle to consistently generate fully valid proteins, the successfully generated proteins show improved novelty and diversity compared to their amino acid-based model counterparts. Furthermore, the all-atom representation achieves structural foldability results comparable to those of amino acid-based models. Lastly, our results highlight the absorbing noise schedule as the most effective for both representations. Data and code are available at https://github.com/Intelligent-molecular-systems/All-Atom-Protein-Sequence-Generation.

NPM1-mutated AML is one of the largest entities in international classification systems of myeloid neoplasms, which are based on integrating morphologic and clinical data with genomic data. Previous research, however, indicates that bulk transcriptomics-based subtyping may improve prognostication and therapy guidance. Here, we characterized the heterogeneity in NPM1-mutated AML by performing single-cell RNA-sequencing and spectral flow cytometry on 16 AML belonging to three distinct subtypes previously identified by bulk transcriptomics. Using single-cell expression profiling we generated a comprehensive atlas of NPM1-mutated AML, collectively reconstituting complete myelopoiesis. The three NPM1-mutated transcriptional subtypes showed consistent differences in the proportions of myeloid cell clusters with distinct patterns in lineage commitment and maturational arrest. In all samples, leukemic cells were detected across different myeloid cell clusters, indicating that NPM1-mutated AML are heavily skewed but not fully arrested in myelopoiesis. Same-sample multi-color spectral flow cytometry recapitulated these skewing patterns, indicating that the three NPM1-mutated subtypes can be consistently identified across platforms. Moreover, our analyses highlighted differences in the abundance of rare hematopoietic stem cells suggesting that skewing occurs early in myelopoiesis. To conclude, by harnessing single-cell RNA-sequencing and spectral flow cytometry, we provide a detailed description of three distinct and reproducible patterns in lineage skewing in NPM1-mutated AML that may have potential relevance for prognosis and treatment of patients with NPM1-mutated AML. (Figure presented.)

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.

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.

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.

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
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.

Background: The accumulation of amyloid-β1−42 (Aβ1−42) peptides and phosphorylated-Tau181 (p-Tau181) tangles from the preclinical stages of Alzheimer’s disease (AD) has led to a biological definition of the disease. However, among Aβ1−42-positive individuals, cognitive decline onset varies, and some never develop symptoms. Genetic influences on molecular pathways and their interactions with proteinopathy may underlie this heterogeneity. Leveraging data from a large sample of cognitively intact older adults in the European Prevention of Alzheimer Dementia (EPAD) cohort, we examined how AD-related pathophysiological changes (i.e., Aβ1−42 and p-Tau181), polygenic pathways and their interaction are associated with WM micro- and macrostructural properties. Methods: We selected 803 individuals (mean age = 64.7 ± 7.3 years, 458 [57.0%] females, 275 [34.2%] APOE-ε4 carriers) with CSF-Aβ1−42 and p-Tau181 measurements available, full genotyping, and structural and diffusion MRI. Polygenic risk scores (PRSs) were computed using 85 AD-related genetic variants. These were mapped to their corresponding genes and, after excluding those belonging to the APOE locus, clustered by function into six pathway-specific PRSs (i.e., immune activation, signal transduction, inflammation, lipid, amyloid, and clearance pathways). Diffusion MRIs were processed through the fixel-based analysis framework to derive fiber density (FD) and fiber cross-section (FC) metrics, which were averaged within WM tracts. Linear models assessed the effects of AD-related pathophysiological changes, global and pathway-specific PRSs, and their interactions on FD and FC at both the tract and fixel levels. Models were corrected for multiple comparisons. Results: P-Tau181 was primarily associated with greater FD. The lipid pathway was associated with greater FD and FC, with these effects predominantly occurring in the left hemisphere, consistent with evidence of hemispheric dominance. The clearance pathway moderated the effect of Aβ1−42 on FD, with a positive slope in A + compared to A- individuals. The immune activation pathway moderated the effect of p-Tau181 on FD, with a negative slope in T + compared to T- individuals. Conclusions: Pathway-specific genetic vulnerability to AD is associated with alterations in WM tracts both directly and by moderating the effects of AD-related pathophysiological changes. AD-associated genetic risk should be integrated into the AD diagnostic framework to enable targeted screening and intervention for future preclinical trials aimed at specific biological pathways.

Peptides are versatile tools in immunotherapy, serving as vaccines and targets for specific immunotherapeutic strategies. Peptides engage immune cells like macrophages and T cells, enabling precise modulation of immune responses. In this context, we highlight the utility of macrophages, innate immune cells involved in constant surveillance, for detecting their phagolysosomal content as a minimally-invasive biomarker strategy. Analyzing proteolytic patterns in phagolysosomes offers a high-sensitivity approach to assess tissue homeostasis and tissue disruption, such as in cancer. Despite their potential, a major challenge lies in the lack of comprehensive tools for predicting cutting sites across phagolysosomal proteases. Therefore, we developed the computational tool PredLyP (abbreviation for “prediction of lysosomal proteases”) to identify cutting sites of phagolysosomal proteases, which are essential enzymes involved in protein degradation within (phago)lysosomes, to predict the potential peptides generated from the input proteins. Unlike existing tools, PredLyP utilizes Position Specific Scoring Matrices derived from amino acid sequences, physical (charge and hydropathy) and structural (secondary structure and solvent accessibility) features. Moreover, it incorporates a sequential cutter functionality that mimics the ordered action of proteases, providing predictive insights into substrate fragment generation. Comparisons with other tools demonstrate the superior sensitivity of PredLyP, enabling accurate prediction of complete and partial digestion fragments, a critical requirement for real-world applications in proteomics, antibody development, and immune system research. Overall, PredLyP represents a robust tool for advancing our understanding of proteolytic processes in phagolysosomes and their implications in health and disease.

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.

Background: Fractional exhaled nitric oxide (FeNO) is a noninvasive method to determine the degree of airway inflammation. Handheld devices such as the Vivatmo Me are used for home monitoring. Differences were found between the Vivatmo Me and standard measurements with the NIOX VERO. Therefore, we aimed to determine the accuracy of the Vivatmo Me for FeNO measurements. Methods: Adult patients with an appointment for FeNO-measurement according to regular care, were invited to perform the FeNO measurement with both devices. From these measurements the FeNO values were compared, and the device user-friendliness was determined. Results: One hundred and sixty-four patients were included. The number of attempts needed for a successful measurement and the failure rate were higher with the Vivatmo Me. Although the measurements were highly correlated, a significant difference (p < 0.001) was found between FeNO values measured with both devices. From the Vivatmo measurements, 32% did not fall within the claimed accuracy ranges. A linear correction on the FeNO values reduced this number. Conclusion: Our findings indicate that the Vivatmo Me does not comply with the claimed accuracy of clinical FeNO measurements and the measurement is challenging to perform. By applying the proposed correction, the comparative validity of the FeNO measurement improves and therefore its clinical usefulness.

Background
Protein truncating variants (PTVs) in SORL1 are observed almost exclusively in Alzheimer’s Disease (AD) cases, but the effect of rare SORL1 missense variants is unclear.

Methods
To identify high-priority missense variants (HPVs), we applied ‘domain mapping of disease mutations’ for the 637 unique coding SORL1 variants detected in 18,959 AD-cases and 21,893 non-demented controls.

Results
In this sample, PTVs and HPVs associated with respectively a 35- and 10-fold increased risk of early onset AD and 17- and 6-fold increased risk of overall AD. The median age at onset (AAO) of PTV- and HPV-carriers was 62 and 64 years, and APOE-genotype contributed to AAO-variability. The median AAO of PTV- and HPV-carriers is ~8–10 years earlier than wild-type SORL1 carriers, matched for APOE-genotype. Specific HPVs are highly penetrant and lead to earlier AAOs than PTVs, suggesting possible dominant negative effects.

Conclusion
Our results justify a debate on whether HPV carriers should be considered for clinical counseling.

While biological age scores have been shown to characterize aging by estimating chronological age based on physiological biomarkers, interactions between different age scores are largely unknown. To study this, large-scale multi-modal data are crucial. However, such data are scarce as population-based cohorts are generally restricted in sharing their data. Here, we employ federated learning to study the relationship between the two types of biological age scores: BrainAge based on brain MRI and MetaboAge based on metabolites. Using three large population-based cohorts, we trained a federated deep learning model to estimate BrainAge and compared its performance to models trained in a single cohort. The federated BrainAge model yielded significantly lower error for age prediction across the cohorts than locally trained models. Harmonizing the age interval between cohorts further improved BrainAge accuracy. Subsequently, we compared BrainAge and MetaboAge by performing association analysis and survival analysis for dementia and mortality prediction to further characterize both scores. The association analysis showed a weak association between BrainAge and MetaboAge, while the survival analysis indicated complementary predictive values for the mortality risk of the two scores. Federated learning has been shown to be a valuable technique for enabling the use of research cohorts that are restricted in data sharing. We conclude that BrainAge and MetaboAge act synergetically for the prediction of time to all-cause mortality, and both aging scores capture different aspects of the aging process.

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
Genome-wide association studies (GWAS) linked TMEM106B variants to susceptibility for neurodegenerative diseases, but the causal genetic elements remain unclear.

Method
We used genotyping data from 5,792 Alzheimer disease cases and controls, and applied COJO to identify haplotypes in the TMEM106B locus that independently associated with AD. Then, we used long-read sequencing data from 513 individuals to annotate these haplotypes with structural variations that map into them.

Results
Analysis of the genotyping data revealed that the TMEM106B locus consists of four major haplotypes: HA/Ha (covering the coding region), and HB/Hb (covering the upstream regulatory region). These combine into four combinations with varying population-frequencies: HAB (57%), HaB (34%), Hab (9%), and HAb (<1%). Long-read sequencing of 513 individuals showed that HA haplotypes (marked by 185-Threonine) carry unique methylated CpG sites and an AluYb8-retrotransposon in the 3' UTR, while the Ha haplotypes are marked by the 185-Serine allele. Hb haplotypes carry several structural variants (SVs) in nearby distal enhancers, including a 19 Kbp rearrangement, absent in all other haplotypes. Joint association models revealed that the HAB combination (AluYb8+185-Threonine) is risk-increasing, while Hab (SVs+185-Serine) confers the protective effect. HaB (185-Serine only) is neutral, while HAb was too rare to assess. Relative to middle-aged non-demented controls, cognitively healthy centenarians were more enriched with Hab (OR=1.49, padj=2.18×10-2) than with HaB (OR=1.23, padj=5.06×10-2). Proteomic analysis of temporal cortex tissues (n = 182) indicated that relative to the neutral HaB combination, the protective Hab is associated with 1.1-fold lower TMEM106B C-terminal peptide abundance, while the risk-increasing HAB is associated with 1.16-fold higher abundance.

Conclusion
Our data indicates that the genetic structure underlying the association of the TMEM106B locus with neurodegenerative diseases is driven by the effect of multiple haplotypes.

We constructed a polygenic protective score specific to Alzheimer’s disease (AD PPS) based on the current literature among the participants enrolled in five studies of healthy aging and extreme longevity in the USA, Europe, and Asia. This AD PPS did not include variants on apolipoprotein E (APOE) gene. Comparisons of AD PPS in different data sets of healthy agers and centenarians showed that centenarians have stronger genetic protection against AD compared to individuals without familial longevity. The current study also shows evidence that this genetic protection increases with increasingly older ages in centenarians (centenarians who died before reaching age 105 years, semi-supercentenarians who reached age 105 to 109 years, and supercentenarians who reached age 110 years and older). However, the genetic protection was of modest size: the average increase in AD PPS was approximately one additional protective allele per 5 years of gained lifetime. Additionally, we show that the higher AD PPS was associated with better cognitive function and decreased mortality. Taken together, this analysis suggests that individuals who achieve the most extreme ages, on average, have the greatest protection against AD. This finding is robust to different genetic backgrounds with important implications for universal applicability of therapeutics that target this AD PPS.

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.

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.

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.

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. [...]