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

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

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

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

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

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

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

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

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

Across the lifespan, diet and physical activity profiles substantially influence immunometabolic health. DNA methylation, as a tissue-specific marker sensitive to behavioral change, may mediate these effects through modulation of transcription factor binding and subsequent gene expression. Despite this, few human studies have profiled DNA methylation and gene expression simultaneously in multiple tissues or examined how molecular levels react and interact in response to lifestyle changes. The Growing Old Together (GOTO) study is a 13-week lifestyle intervention in older adults, which imparted health benefits to participants. Here, we characterize the DNA methylation response to this intervention at over 750 thousand CpGs in muscle, adipose, and blood. Differentially methylated sites are enriched for active chromatin states, located close to relevant transcription factor binding sites, and associated with changing expression of insulin sensitivity genes and health parameters. In addition, measures of biological age are consistently reduced, with decreases in grimAge associated with observed health improvements. Taken together, our results identify responsive molecular markers and demonstrate their potential to measure progression and finetune treatment of age-related risks and 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

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

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

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

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

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

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

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

T cell development in the mouse thymus has been studied extensively, but less is known regarding T cell development in the human thymus. We used a combination of single-cell techniques and functional assays to perform deep immune profiling of human T cell development, focusing on the initial stages of prelineage commitment. We identified three thymus-seeding progenitor populations that also have counterparts in the bone marrow. In addition, we found that the human thymus physiologically supports the development of monocytes, dendritic cells, and NK cells, as well as limited development of B cells. These results are an important step toward monitoring and guiding regenerative therapies in patients after hematopoietic stem cell transplantation.

@en

Sustained night shift work is associated with various adverse health risks, including an increased risk of cardiovascular disease, type II diabetes, and susceptibility to infectious respiratory diseases. The extent of these adverse health effects, however, seems to greatly vary between night shift workers, yet the underlying reasons and the mechanisms underlying these interindividual differences remain poorly understood. Metabolomics assays in the blood have recently gained much attention as a minimally invasive biomarker platform capturing information predictive of metabolic and cardiovascular diseases. In this cross-sectional study, we explored and compared the metabolic profiles of 1010 night shift workers and 1010 age- and sex-matched day workers (non-shift workers) from the Lifelines Cohort Study. The metabolic profiles were determined using the ¹H-NMR Nightingale platform for the quantification of 250 parameters of metabolism, including routine lipids, extensive lipoprotein subclasses, fatty acid composition, and various low-molecular metabolites, including amino acids, ketone bodies, and gluconeogenesis-related metabolites. Night shift workers had an increased BMI (26.6 vs. 25.9 kg/m²) compared with day workers (non-shift workers) in both sexes, were slightly more likely to be ever smokers (only in males) (54% vs. 46%), worked on average 5.9 ± 3.7 night shifts per month, and had been working in night shifts for 18.3 ± 10.5 years on average. We observed changes in several metabolic markers in male night shift workers compared with non-shift workers, but no changes were observed in women. In men, we observed higher levels of glycoprotein acetyls (GlycA), triglycerides, and fatty acids compared with non-shift workers. The changes were seen in the ratio of triglycerides and cholesterol(esters) to total lipids in different sizes of VLDL particles. Glycoprotein acetyls (GlycAs) are of particular interest as markers since they are known as biomarkers for low-grade chronic inflammation. When the analyses were adjusted for BMI, no significant associations were observed. Further studies are needed to better understand the relationship between night shift work and metabolic profiles, particularly with respect to the role of sex and BMI in this relationship.

@en