Prognostic information is needed to balance benefits and risks of cancer treatment in older patients. Metabolomics-based scores were previously developed to predict 5- and 10-year mortality (MetaboHealth) and biological age (MetaboAge). This study aims to investigate the association of MetaboHealth and MetaboAge with 1-year mortality in older patients with solid tumors, and to study their predictive value for mortality in addition to established clinical predictors. This prospective cohort study included patients aged ≥ 70 years with a solid malignant tumor, who underwent blood sampling and a geriatric assessment before treatment initiation. The outcome was all-cause 1-year mortality. Of the 192 patients, the median age was 77 years. With each SD increase of MetaboHealth, patients had a 2.32 times increased risk of mortality (HR 2.32, 95% CI 1.59–3.39). With each year increase in MetaboAge, there was a 4% increased risk of mortality (HR 1.04, 1.01–1.07). MetaboHealth and MetaboAge showed an AUC of 0.66 (0.56–0.75) and 0.60 (0.51–0.68) for mortality prediction accuracy, respectively. The AUC of a predictive model containing age, primary tumor site, distant metastasis, comorbidity, and malnutrition was 0.76 (0.68–0.83). Addition of MetaboHealth increased AUC to 0.80 (0.74–0.87) (p = 0.09) and AUC did not change with MetaboAge (0.76 (0.69–0.83) (p = 0.89)). Higher MetaboHealth and MetaboAge scores were associated with 1-year mortality. The addition of MetaboHealth to established clinical predictors only marginally improved mortality prediction in this cohort with various types of tumors. MetaboHealth may potentially improve identification of older patients vulnerable for adverse events, but numbers were too small for definitive conclusions. The TENT study is retrospectively registered at the Netherlands Trial Register (NTR), trial number NL8107. Date of registration: 22–10-2019.@en

Children with Multisystem Inflammatory Syndrome in Children (MIS-C) can present with thrombocytopenia, which is a key feature of hemophagocytic lymphohistiocytosis (HLH). We hypothesized that thrombocytopenic MIS-C patients have more features of HLH. Clinical characteristics and routine laboratory parameters were collected from 228 MIS-C patients, of whom 85 (37%) were thrombocytopenic. Thrombocytopenic patients had increased ferritin levels; reduced leukocyte subsets; and elevated levels of ASAT and ALAT. Soluble IL-2RA was higher in thrombocytopenic children than in non-thrombocytopenic children. T-cell activation, TNF-alpha and IFN-gamma signaling markers were inversely correlated with thrombocyte levels, consistent with a more pronounced cytokine storm syndrome. Thrombocytopenia was not associated with severity of MIS-C and no pathogenic variants were identified in HLH-related genes. This suggests that thrombocytopenia in MIS-C is not a feature of a more severe disease phenotype, but the consequence of a distinct hyperinflammatory immunopathological process in a subset of children.@en

The curative benefits of autologous and allogeneic transplantation of hematopoietic stem cells (HSCs) have been proven in various diseases. However, the low number of true HSCs that can be collected from patients and the subsequent in vitro maintenance and expansion of true HSCs for genetic correction remains challenging. Addressing this issue, we here focused on optimizing culture conditions to improve ex vivo expansion of true HSCs for gene therapy purposes. In particular, we explored the use of epigenetic regulators to enhance the effectiveness of HSC-based lentiviral (LV) gene therapy. The histone deacetylase inhibitor quisinostat and bromodomain inhibitor CPI203 each promoted ex vivo expansion of functional HSCs, as validated by xenotransplantation assays and single-cell RNA sequencing analysis. We confirmed the stealth effect of LV transduction on the loss of HSC numbers in commonly used culture protocols, whereas the addition of quisinostat or CPI203 improved the expansion of HSCs in transduction protocols. Notably, we demonstrated that the addition of quisinostat improved the LV transduction efficiency of HSCs and early progenitors. Our suggested culture conditions highlight the potential therapeutic effects of epigenetic regulators in HSC biology and their clinical applications to advance HSC-based gene correction.@en

Subtyping of acute myeloid leukaemia (AML) is predominantly based on recurrent genetic abnormalities, but recent literature indicates that transcriptomic phenotyping holds immense potential to further refine AML classification. Here we integrated five AML transcriptomic datasets with corresponding genetic information to provide an overview (n = 1224) of the transcriptomic AML landscape. Consensus clustering identified 17 robust patient clusters which improved identification of CEBPA-mutated patients with favourable outcomes, and uncovered transcriptomic subtypes for KMT2A rearrangements (2), NPM1 mutations (5), and AML with myelodysplasia-related changes (AML-MRC) (5). Transcriptomic subtypes of KMT2A, NPM1 and AML-MRC showed distinct mutational profiles, cell type differentiation arrests and immune properties, suggesting differences in underlying disease biology. Moreover, our transcriptomic clusters show differences in ex-vivo drug responses, even when corrected for differentiation arrest and superiorly capture differences in drug response compared to genetic classification. In conclusion, our findings underscore the importance of transcriptomics in AML subtyping and offer a basis for future research and personalised treatment strategies. Our transcriptomic compendium is publicly available and we supply an R package to project clusters to new transcriptomic studies.

@en

The diagnostic spectrum for AML patients is increasingly based on genetic abnormalities due to their prognostic and predictive value. However, information on the AML blast phenotype regarding their maturational arrest has started to regain importance due to its predictive power for drug responses. Here, we deconvolute 1350 bulk RNA-seq samples from five independent AML cohorts on a single-cell healthy BM reference and demonstrate that the morphological differentiation stages (FAB) could be faithfully reconstituted using estimated cell compositions (ECCs). Moreover, we show that the ECCs reliably predict ex-vivo drug resistances as demonstrated for Venetoclax, a BCL-2 inhibitor, resistance specifically in AML with CD14+ monocyte phenotype. We validate these predictions using LUMC proteomics data by showing that BCL-2 protein abundance is split into two distinct clusters for NPM1-mutated AML at the extremes of CD14+ monocyte percentages, which could be crucial for the Venetoclax dosing patients. Our results suggest that Venetoclax resistance predictions can also be extended to AML without recurrent genetic abnormalities and possibly to MDS-related and secondary AML. Lastly, we show that CD14+ monocytic dominated Ven/Aza treated patients have significantly lower overall survival. Collectively, we propose a framework for allowing a joint mutation and maturation stage modeling that could be used as a blueprint for testing sensitivity for new agents across the various subtypes of AML.

@en

Background
1H-NMR metabolomics and DNA methylation in blood are widely known biomarkers predicting age-related physiological decline and mortality yet exert mutually independent mortality and frailty signals.

Methods
Leveraging multi-omics data in four Dutch population studies (N = 5238, ∼40% of which male) we investigated whether the mortality signal captured by 1H-NMR metabolomics could guide the construction of DNA methylation-based mortality predictors.

Findings
We trained DNA methylation-based surrogates for 64 metabolomic analytes and found that analytes marking inflammation, fluid balance, or HDL/VLDL metabolism could be accurately reconstructed using DNA-methylation assays. Interestingly, a previously reported multi-analyte score indicating mortality risk (MetaboHealth) could also be accurately reconstructed. Sixteen of our derived surrogates, including the MetaboHealth surrogate, showed significant associations with mortality, independent of relevant covariates.

Interpretation
The addition of our metabolic analyte-derived surrogates to the well-established epigenetic clock GrimAge demonstrates that our surrogates potentially represent valuable mortality signal.

Funding
BBMRI-NL, X-omics, VOILA, Medical Delta, NWO, ERC.@en

Existing metabolomic clocks exhibit deficiencies in capturing the heterogeneous aging rates among individuals with the same chronological age. Yet, the modifiable and non-modifiable factors in metabolomic aging have not been systematically studied. Here, a new aging measure—MetaboAgeMort—is developed using metabolomic profiles from 239,291 UK Biobank participants for 10-year all-cause mortality prediction. The MetaboAgeMort showed significant associations with all-cause mortality, cause-specific mortality, and diverse incident diseases. Adding MetaboAgeMort to a conventional risk factors model improved the predictive ability of 10-year mortality. A total of 99 modifiable factors across seven categories are identified for MetaboAgeMort. Among these, 16 factors representing pulmonary function, body composition, socioeconomic status, dietary quality, smoking status, alcohol intake, and disease status showed quantitatively stronger associations. The genetic analyses revealed 99 genomic risk loci and 271 genes associated with MetaboAgeMort. The tissue-enrichment analysis showed significant enrichment in liver. While the external validation of the MetaboAgeMort is required, this study illuminates heterogeneous metabolomic aging across the same age, providing avenues for identifying high-risk individuals, developing anti-aging therapies, and personalizing interventions, thus promoting healthy aging and longevity.@en

To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

@en

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

Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1–7. This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases8–11. Here we present a genome-wide association study for 233 circulating metabolic traits quantified by nuclear magnetic resonance spectroscopy in up to 136,016 participants from 33 cohorts. We identify more than 400 independent loci and assign probable causal genes at two-thirds of these using manual curation of plausible biological candidates. We highlight the importance of sample and participant characteristics that can have significant effects on genetic associations. We use detailed metabolic profiling of lipoprotein- and lipid-associated variants to better characterize how known lipid loci and novel loci affect lipoprotein metabolism at a granular level. We demonstrate the translational utility of comprehensively phenotyped molecular data, characterizing the metabolic associations of intrahepatic cholestasis of pregnancy. Finally, we observe substantial genetic pleiotropy for multiple metabolic pathways and illustrate the importance of careful instrument selection in Mendelian randomization analysis, revealing a putative causal relationship between acetone and hypertension. Our publicly available results provide a foundational resource for the community to examine the role of metabolism across diverse diseases.@en

Aging is a multifaceted and intricate physiological process characterized by a gradual decline in functional capacity, leading to increased susceptibility to diseases and mortality. While chronological age serves as a strong risk factor for age-related health conditions, considerable heterogeneity exists in the aging trajectories of individuals, suggesting that biological age may provide a more nuanced understanding of the aging process. However, the concept of biological age lacks a clear operationalization, leading to the development of various biological age predictors without a solid statistical foundation. This paper addresses these limitations by proposing a comprehensive operationalization of biological age, introducing the “AccelerAge” framework for predicting biological age, and introducing previously underutilized evaluation measures for assessing the performance of biological age predictors. The AccelerAge framework, based on Accelerated Failure Time (AFT) models, directly models the effect of candidate predictors of aging on an individual’s survival time, aligning with the prevalent metaphor of aging as a clock. We compare predictors based on the AccelerAge framework to a predictor based on the GrimAge predictor, which is considered one of the best-performing biological age predictors, using simulated data as well as data from the UK Biobank and the Leiden Longevity Study. Our approach seeks to establish a robust statistical foundation for biological age clocks, enabling a more accurate and interpretable assessment of an individual’s aging status.

@en

Traditionally, flow cytometry has been the preferred method to characterize immune cells at the single-cell level. Flow cytometry is used in immunology mostly to measure the expression of identifying markers on the cell surface, but—with good antibodies—can also be used to assess the expression of intracellular proteins. The advent of single-cell RNA-sequencing has paved the road to study immune development at an unprecedented resolution. Single-cell RNA-sequencing studies have not only allowed us to efficiently chart the make-up of heterogeneous tissues, including their most rare cell populations, it also increasingly contributes to our understanding how different omics modalities interplay at a single cell resolution. Particularly for investigating the immune system, this means that these single-cell techniques can be integrated to combine and correlate RNA and protein data at the single-cell level. While RNA data usually reveals a large heterogeneity of a given population identified solely by a combination of surface protein markers, the integration of different omics modalities at a single cell resolution is expected to greatly contribute to our understanding of the immune system.

@en

1H-NMR metabolomics data is increasingly used to track health and disease. Nightingale Health, a major supplier of 1H-NMR metabolomics, has recently updated the quantification strategy to further align with clinical standards. Such updates, however, might influence backward replicability, particularly affecting studies with repeated measures. Using data from BBMRI-NL consortium (~28,000 samples from 28 cohorts), we compared Nightingale data, originally released in 2014 and 2016, with a re-quantified version released in 2020, of which both versions were based on the same NMR spectra. Apart from two discontinued and twenty-three new analytes, we generally observe a high concordance between quantification versions with 73 out of 222 (33%) analytes showing a mean ρ > 0.9 across all cohorts. Conversely, five analytes consistently showed lower Spearman’s correlations (ρ < 0.7) between versions, namely acetoacetate, LDL-L, saturated fatty acids, S-HDL-C, and sphingomyelins. Furthermore, previously trained multi-analyte scores, such as MetaboAge or MetaboHealth, might be particularly sensitive to platform changes. Whereas MetaboHealth replicated well, the MetaboAge score had to be retrained due to use of discontinued analytes. Notably, both scores in the re-quantified data recapitulated mortality associations observed previously. Concluding, we urge caution in utilizing different platform versions to avoid mixing analytes, having different units, or simply being discontinued.@en

Glycans are an essential structural component of immunoglobulin G (IgG) that modulate its structure and function. However, regulatory mechanisms behind this complex posttranslational modification are not well known. Previous genome-wide association studies (GWAS) identified 29 genomic regions involved in regulation of IgG glycosylation, but only a few were functionally validated. One of the key functional features of IgG glycosylation is the addition of galactose (galactosylation), a trait which was shown to be associated with ageing. We performed GWAS of IgG galactosylation (N=13,705) and identified 16 significantly associated loci, indicating that IgG galactosylation is regulated by a complex network of genes that extends beyond the galactosyltransferase enzyme that adds galactose to IgG glycans. Gene prioritization identified 37 candidate genes. Using a recently developed CRISPR/dCas9 system we manipulated gene expression of candidate genes in the in vitro IgG expression system. Upregulation of three genes, EEF1A1, MANBA and TNFRSF13B, changed the IgG glycome composition, which confirmed that these three genes are involved in IgG galactosylation in this in vitro expression system.

@en

Biological age captures a person's age-related risk of unfavorable outcomes using biophysiological information. Multivariate biological age measures include frailty scores and molecular biomarkers. These measures are often studied in isolation, but here we present a large-scale study comparing them. In 2 prospective cohorts (n = 3 222), we compared epigenetic (DNAm Horvath, DNAm Hannum, DNAm Lin, DNAm epiTOC, DNAm PhenoAge, DNAm DunedinPoAm, DNAm GrimAge, and DNAm Zhang) and metabolomic-based (MetaboAge and MetaboHealth) biomarkers in reflection of biological age, as represented by 5 frailty measures and overall mortality. Biomarkers trained on outcomes with biophysiological and/or mortality information outperformed age-trained biomarkers in frailty reflection and mortality prediction. DNAm GrimAge and MetaboHealth, trained on mortality, showed the strongest association with these outcomes. The associations of DNAm GrimAge and MetaboHealth with frailty and mortality were independent of each other and of the frailty score mimicking clinical geriatric assessment. Epigenetic, metabolomic, and clinical biological age markers seem to capture different aspects of aging. These findings suggest that mortality-trained molecular markers may provide novel phenotype reflecting biological age and strengthen current clinical geriatric health and well-being assessment.

@en

Analytical techniques with high sensitivity and selectivity are essential to the quantitative analysis of clinical samples. Liquid chromatography coupled to tandem mass spectrometry is the gold standard in clinical chemistry. However, tandem mass spectrometers come at high capital expenditure and maintenance costs. We recently showed that it is possible to generate very similar results using a much simpler single mass spectrometry detector by performing enhanced in-source fragmentation/annotation (EISA) combined with correlated ion monitoring. Here we provide a step-by-step protocol for optimizing the analytical conditions for EISA, so anyone properly trained in liquid chromatography–mass spectrometry can follow and apply this technique for any given analyte. We exemplify the approach by using 2-hydroxyglutarate (2-HG) which is a clinically relevant metabolite whose d-enantiomer is considered an ‘oncometabolite’, characteristic of cancers associated with mutated isocitrate dehydrogenases 1 or 2 (IDH1/2). We include procedures for determining quantitative robustness, and show results of these relating to the analysis of dl-2-hydroxyglutarate in cells, as well as in serum samples from patients with acute myeloid leukemia that contain the IDH1/2 mutation. This EISA–mass spectrometry protocol is a broadly applicable and low-cost approach for the quantification of small molecules that has been developed to work well for both single-quadrupole and time-of-flight mass analyzers.

@en

Skeletal muscles are composed of different myofiber types characterized by the expression of myosin heavy chain isoforms, which can be affected by physical activity, aging, and pathological conditions. Here, we present a step-by-step high-throughput semi-automated approach for performing myofiber type quantification of entire human or mouse muscle tissue sections, including immunofluorescence staining, image acquisition, processing, and quantification. For complete details on the use and execution of this protocol, please refer to Abbassi-Daloii et al. (2022).¹

@en

The first successful European hematopoietic stem cell transplantation (HSCT) was performed in 1968 as treatment in a newborn with IL2RG deficiency using an HLA-identical sibling donor. Because of declining naive T and natural killer (NK) cells, and persistent human papilloma virus (HPV)-induced warts, the patient received a peripheral stem cell boost at the age of 37 years. NK and T cells were assessed before and up to 14 years after the boost by flow cytometry. The boost induced renewed reconstitution of functional NK cells that were 14 years later enriched for CD56^dimCD27⁺ NK cells. T-cell phenotype and T-cell receptor (TCR) repertoire were simultaneously analyzed by including TCR Vβ antibodies in the cytometry panel. Naive T-cell numbers with a diverse TCR Vβ repertoire were increased by the boost. Before and after the boost, clonal expansions with a homogeneous TIGIT and PD-1 phenotype were identified in the CD27⁻ and/or CD28⁻ memory population in the patient, but not in the donor. TRB sequencing was applied on sorted T-cell subsets from blood and on T cells from skin biopsies. Abundant circulating CD8 memory clonotypes with a chronic virus-associated CD57⁺KLRG1⁺CX3CR1⁺ phenotype were also present in warts, but not in healthy skin of the patient, suggesting a link with HPV. In conclusion, we demonstrate in this IL2RG-deficient patient functional NK cells, a diverse and lasting naive T-cell compartment, supported by a stem cell boost, and an oligoclonal memory compartment half a century after HSCT.

@en

Objective: To facilitate patient disease subset and risk factor identification by constructing a pipeline which is generalizable, provides easily interpretable results, and allows replication by overcoming electronic health records (EHRs) batch effects. Material and Methods: We used 1872 billing codes in EHRs of 102 880 patients from 12 healthcare systems. Using tools borrowed from single-cell omics, we mitigated center-specific batch effects and performed clustering to identify patients with highly similar medical history patterns across the various centers. Our visualization method (PheSpec) depicts the phenotypic profile of clusters, applies a novel filtering of noninformative codes (Ranked Scope Pervasion), and indicates the most distinguishing features. Results: We observed 114 clinically meaningful profiles, for example, linking prostate hyperplasia with cancer and diabetes with cardiovascular problems and grouping pediatric developmental disorders. Our framework identified disease subsets, exemplified by 6 "other headache"clusters, where phenotypic profiles suggested different underlying mechanisms: migraine, convulsion, injury, eye problems, joint pain, and pituitary gland disorders. Phenotypic patterns replicated well, with high correlations of ≥0.75 to an average of 6 (2-8) of the 12 different cohorts, demonstrating the consistency with which our method discovers disease history profiles. Discussion: Costly clinical research ventures should be based on solid hypotheses. We repurpose methods from single-cell omics to build these hypotheses from observational EHR data, distilling useful information from complex data. Conclusion: We establish a generalizable pipeline for the identification and replication of clinically meaningful (sub)phenotypes from widely available high-dimensional billing codes. This approach overcomes datatype problems and produces comprehensive visualizations of validation-ready phenotypes.

@en

Background: Missing or incomplete phenotypic information can severely deteriorate the statistical power in epidemiological studies. High-throughput quantification of small-molecules in bio-samples, i.e. ‘metabolomics’, is steadily gaining popularity, as it is highly informative for various phenotypical characteristics. Here we aim to leverage metabolomics to impute missing data in clinical variables routinely assessed in large epidemiological and clinical studies. Methods: To this end, we have employed ∼26,000 ¹H-NMR metabolomics samples from 28 Dutch cohorts collected within the BBMRI-NL consortium, to create 19 metabolomics-based predictors for clinical variables, including diabetes status (AUC_{5-Fold CV} = 0·94) and lipid medication usage (AUC_{5-Fold CV} = 0·90). Findings: Subsequent application in independent cohorts confirmed that our metabolomics-based predictors can indeed be used to impute a wide array of missing clinical variables from a single metabolomics data resource. In addition, application highlighted the potential use of our predictors to explore the effects of totally unobserved confounders in omics association studies. Finally, we show that our predictors can be used to explore risk factor profiles contributing to mortality in older participants. Interpretation: To conclude, we provide ¹H-NMR metabolomics-based models to impute clinical variables routinely assessed in epidemiological studies and illustrate their merit in scenarios when phenotypic variables are partially incomplete or totally unobserved. Funding: BBMRI-NL, X-omics, VOILA, Medical Delta and the Dutch Research Council (NWO-VENI).

@en