Search

Background: 1 H-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 1 H-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., Competing Interests: Declaration of interests Authors declare no competing interests., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

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., (© 2024. The Author(s).)

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., Competing Interests: Declaration of Competing Interest The authors declare that they have no competing interests., (Copyright © 2024. Published by Elsevier B.V.)

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., Competing Interests: Declaration of competing interest The authors declare no conflict of interests., (Copyright © 2024. Published by Elsevier Inc.)

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., (© 2024. The Author(s).)

Abstract: Allogeneic stem cell transplantation (alloSCT) is a curative treatment for hematological malignancies. After HLA-matched alloSCT, antitumor immunity is caused by donor T cells recognizing polymorphic peptides, designated minor histocompatibility antigens (MiHAs), that are presented by HLA on malignant patient cells. However, T cells often target MiHAs on healthy nonhematopoietic tissues of patients, thereby inducing side effects known as graft-versus-host disease. Here, we aimed to identify the dominant repertoire of HLA-I-restricted MiHAs to enable strategies to predict, monitor or modulate immune responses after alloSCT. To systematically identify novel MiHAs by genome-wide association screening, T-cell clones were isolated from 39 transplanted patients and tested for reactivity against 191 Epstein-Barr virus transformed B cell lines of the 1000 Genomes Project. By discovering 81 new MiHAs, we more than doubled the antigen repertoire to 159 MiHAs and demonstrated that, despite many genetic differences between patients and donors, often the same MiHAs are targeted in multiple patients. Furthermore, we showed that one quarter of the antigens are cryptic, that is translated from unconventional open reading frames, for example long noncoding RNAs, showing that these antigen types are relevant targets in natural immune responses. Finally, using single cell RNA-seq data, we analyzed tissue expression of MiHA-encoding genes to explore their potential role in clinical outcome, and characterized 11 new hematopoietic-restricted MiHAs as potential targets for immunotherapy. In conclusion, we expanded the repertoire of HLA-I-restricted MiHAs and identified recurrent, cryptic and hematopoietic-restricted antigens, which are fundamental to predict, follow or manipulate immune responses to improve clinical outcome after alloSCT., (© 2024 American Society of Hematology. Published by Elsevier Inc. Licensed under Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0), permitting only noncommercial, nonderivative use with attribution. All other rights reserved.)

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., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism 1-7 . This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases 8-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., (© 2024. The Author(s).)

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., (© 2024. The Author(s).)

1 H-NMR metabolomics data is increasingly used to track health and disease. Nightingale Health, a major supplier of 1 H-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.

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., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Gerontological Society of America.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Cordes, Pike-Overzet, Van Den Akker, Staal and Canté-Barrett.)

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., (© 2023. Springer Nature Limited.)

Objective: To investigate whether established genetic predictors for rheumatoid arthritis (RA) differentiate healthy controls, patients with clinically suspect arthralgia (CSA), and RA patients., Methods: Using analyses of variance, chi-square tests, and mean risk difference analyses, we investigated the association of an RA polygenic risk score (PRS) and HLA shared epitope (HLA-SE) with all participant groups, both unstratified and stratified for anti-citrullinated protein antibody (ACPA) status. We used 3 separate data sets sampled from the same Dutch population (1,015 healthy controls, 479 CSA patients, and 1,146 early classified RA patients). CSA patients were assessed for conversion to inflammatory arthritis over a period of 2 years, after which they were classified as either CSA converters (n = 84) or CSA nonconverters (n = 395)., Results: The PRS was increased in RA patients (mean ± SD PRS 1.31 ± 0.96) compared to the complete CSA group (1.07 ± 0.94) and compared to CSA converters (1.12 ± 0.94). In ACPA- strata, PRS distributions differed strongly when comparing the complete CSA group (mean ± SD PRS 1.05 ± 0.94) and CSA converters (0.97 ± 0.87) to RA patients (1.20 ± 0.94), while in the ACPA+ strata, the complete CSA group (1.25 ± 0.99) differed clearly from healthy controls (1.05 ± 0.94) and RA patients (1.41 ± 0.96). HLA-SE was more prevalent in the RA group (prevalence 0.64) than the complete CSA group (0.45), with small differences between RA patients and CSA converters (0.64 versus 0.60) and larger differences between CSA converters and CSA nonconverters (0.60 versus 0.42). HLA-SE prevalence differed more strongly within the ACPA+ strata as follows: healthy controls (prevalence 0.43), CSA nonconverters (0.48), complete CSA group (0.59), CSA converters (0.66), and RA patients (0.79)., Conclusion: We observed that genetic predisposition increased across pre-RA participant groups. The RA PRS differed in early classified RA and inflammatory pre-disease stages, regardless of ACPA stratification. HLA-SE prevalence differed between arthritis patients, particularly ACPA+ patients, and healthy controls. Genetics seem to fulfill different etiologic roles., (© 2022 The Authors. Arthritis & Rheumatology published by Wiley Periodicals LLC on behalf of American College of Rheumatology.)

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 1 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 2 ) 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.

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.

The response to lifestyle intervention studies is often heterogeneous, especially in older adults. Subtle responses that may represent a health gain for individuals are not always detected by classical health variables, stressing the need for novel biomarkers that detect intermediate changes in metabolic, inflammatory, and immunity-related health. Here, our aim was to develop and validate a molecular multivariate biomarker maximally sensitive to the individual effect of a lifestyle intervention; the Personalized Lifestyle Intervention Status (PLIS). We used 1 H-NMR fasting blood metabolite measurements from before and after the 13-week combined physical and nutritional Growing Old TOgether (GOTO) lifestyle intervention study in combination with a fivefold cross-validation and a bootstrapping method to train a separate PLIS score for men and women. The PLIS scores consisted of 14 and four metabolites for females and males, respectively. Performance of the PLIS score in tracking health gain was illustrated by association of the sex-specific PLIS scores with several classical metabolic health markers, such as BMI, trunk fat%, fasting HDL cholesterol, and fasting insulin, the primary outcome of the GOTO study. We also showed that the baseline PLIS score indicated which participants respond positively to the intervention. Finally, we explored PLIS in an independent physical activity lifestyle intervention study, showing similar, albeit remarkably weaker, associations of PLIS with classical metabolic health markers. To conclude, we found that the sex-specific PLIS score was able to track the individual short-term metabolic health gain of the GOTO lifestyle intervention study. The methodology used to train the PLIS score potentially provides a useful instrument to track personal responses and predict the participant's health benefit in lifestyle interventions similar to the GOTO study., (© 2022 The Authors. The FASEB Journal published by Wiley Periodicals LLC on behalf of Federation of American Societies for Experimental Biology.)

Risk stratification is critical for the early identification of high-risk individuals and disease prevention. Here we explored the potential of nuclear magnetic resonance (NMR) spectroscopy-derived metabolomic profiles to inform on multidisease risk beyond conventional clinical predictors for the onset of 24 common conditions, including metabolic, vascular, respiratory, musculoskeletal and neurological diseases and cancers. Specifically, we trained a neural network to learn disease-specific metabolomic states from 168 circulating metabolic markers measured in 117,981 participants with ~1.4 million person-years of follow-up from the UK Biobank and validated the model in four independent cohorts. We found metabolomic states to be associated with incident event rates in all the investigated conditions, except breast cancer. For 10-year outcome prediction for 15 endpoints, with and without established metabolic contribution, a combination of age and sex and the metabolomic state equaled or outperformed established predictors. Moreover, metabolomic state added predictive information over comprehensive clinical variables for eight common diseases, including type 2 diabetes, dementia and heart failure. Decision curve analyses showed that predictive improvements translated into clinical utility for a wide range of potential decision thresholds. Taken together, our study demonstrates both the potential and limitations of NMR-derived metabolomic profiles as a multidisease assay to inform on the risk of many common diseases simultaneously., (© 2022. The Author(s).)

Motivation: 1H-NMR metabolomics is rapidly becoming a standard resource in large epidemiological studies to acquire metabolic profiles in large numbers of samples in a relatively low-priced and standardized manner. Concomitantly, metabolomics-based models are increasingly developed that capture disease risk or clinical risk factors. These developments raise the need for user-friendly toolbox to inspect new 1H-NMR metabolomics data and project a wide array of previously established risk models., Results: We present MiMIR (Metabolomics-based Models for Imputing Risk), a graphical user interface that provides an intuitive framework for ad hoc statistical analysis of Nightingale Health's 1H-NMR metabolomics data and allows for the projection and calibration of 24 pre-trained metabolomics-based models, without any pre-required programming knowledge., Availability and Implementation: The R-shiny package is available in CRAN or downloadable at https://github.com/DanieleBizzarri/MiMIR, together with an extensive user manual (also available as Supplementary Documents to the article)., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2022. Published by Oxford University Press.)

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 dim CD27 + 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., (© 2022. The Author(s).)

Population-scale expression profiling studies can provide valuable insights into biological and disease-underlying mechanisms. The availability of phenotypic traits is essential for studying clinical effects. Therefore, missing, incomplete, or inaccurate phenotypic information can make analyses challenging and prevent RNA-seq or other omics data to be reused. A possible solution are predictors that infer clinical or behavioral phenotypic traits from molecular data. While such predictors have been developed based on different omics data types and are being applied in various studies, metabolomics-based surrogates are less commonly used than predictors based on DNA methylation profiles.In this study, we inferred 17 traits, including diabetes status and exposure to lipid medication, using previously trained metabolomic predictors. We evaluated whether these metabolomic surrogates can be used as an alternative to reported information for studying the respective phenotypes using expression profiling data of four population cohorts. For the majority of the 17 traits, the metabolomic surrogates performed similarly to the reported phenotypes in terms of effect sizes, number of significant associations, replication rates, and significantly enriched pathways.The application of metabolomics-derived surrogate outcomes opens new possibilities for reuse of multi-omics data sets. In studies where availability of clinical metadata is limited, missing or incomplete information can be complemented by these surrogates, thereby increasing the size of available data sets. Additionally, the availability of such surrogates could be used to correct for potential biological confounding. In the future, it would be interesting to further investigate the use of molecular predictors across different omics types and cohorts., (© 2022. The Author(s).)

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., (© The Author(s) 2022. Published by Oxford University Press on behalf of the American Medical Informatics Association.)

Aging is characterized by increased mortality, functional decline, and exponential increases in the incidence of diseases such as cancer, stroke, cardiovascular disease, neurological disease, respiratory disease, etc. Though the role of aging in these diseases is widely accepted and considered to be a common denominator, the underlying mechanisms are largely unknown. A significant age-related feature observed in many population cohorts is somatic mosaicism, the detectable accumulation of somatic mutations in multiple cell types and tissues, particularly those with high rates of cell turnover (e.g., skin, liver, and hematopoietic cells). Somatic mosaicism can lead to the development of cellular clones that expand with age in otherwise normal tissues. In the hematopoietic system, this phenomenon has generally been referred to as "clonal hematopoiesis of indeterminate potential" (CHIP) when it applies to a subset of clones in which mutations in driver genes of hematologic malignancies are found. Other mechanisms of clonal hematopoiesis, including large chromosomal alterations, can also give rise to clonal expansion in the absence of conventional CHIP driver gene mutations. Both types of clonal hematopoiesis (CH) have been observed in studies of animal models and humans in association with altered immune responses, increased mortality, and disease risk. Studies in murine models have found that some of these clonal events are involved in abnormal inflammatory and metabolic changes, altered DNA damage repair and epigenetic changes. Studies in long-lived individuals also show the accumulation of somatic mutations, yet at this advanced age, carriership of somatic mutations is no longer associated with an increased risk of mortality. While it remains to be elucidated what factors modify this genotype-phenotype association, i.e., compensatory germline genetics, cellular context of the mutations, protective effects to diseases at exceptional age, it points out that the exceptionally long-lived are key to understand the phenotypic consequences of CHIP mutations. Assessment of the clinical significance of somatic mutations occurring in blood cell types for age-related outcomes in human populations of varied life and health span, environmental exposures, and germline genetic risk factors will be valuable in the development of personalized strategies tailored to specific somatic mutations for healthy aging., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Walsh, Raghavachari, Kerr, Bick, Cummings, Druley, Dunbar, Genovese, Goodell, Jaiswal, Maciejewski, Natarajan, Shindyapina, Shuldiner, Van Den Akker and Vijg.)

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 1 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 1 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)., Competing Interests: Declaration of interests The authors declare that there are no competing interests., (Copyright © 2021 The Author(s). Published by Elsevier B.V. All rights reserved.)

To mount an adequate immune response against pathogens, stepwise mutation and selection processes are crucial functions of the adaptive immune system. To better characterize a successful vaccination response, we performed longitudinal (days 0, 5, 7, 10, and 14 after Boostrix vaccination) analysis of the single-cell transcriptome as well as the B-cell receptor (BCR) repertoire (scBCR-rep) in plasma cells of an immunized donor and compared it with baseline B-cell characteristics as well as flow cytometry findings. Based on the flow cytometry knowledge and literature findings, we discriminated individual B-cell subsets in the transcriptomics data and traced over-time maturation of plasmablasts/plasma cells (PB/PCs) and identified the pathways associated with the plasma cell maturation. We observed that the repertoire in PB/PCs differed from the baseline B-cell repertoire e.g., regarding expansion of unique clones in post-vaccination visits, high usage of IGHG1 in expanded clones, increased class-switching events post-vaccination represented by clonotypes spanning multiple IGHC classes and positive selection of CDR3 sequences over time. Importantly, the Variable gene family-based clustering of BCRs represented a similar measure as the gene-based clustering, but certainly improved the clustering of BCRs, as BCRs from duplicated Variable gene families could be clustered together. Finally, we developed a query tool to dissect the immune response to the components of the Boostrix vaccine. Using this tool, we could identify the BCRs related to anti-tetanus and anti-pertussis toxoid BCRs. Collectively, we developed a bioinformatic workflow which allows description of the key features of an ongoing (longitudinal) immune response, such as activation of PB/PCs, Ig class switching, somatic hypermutation, and clonal expansion, all of which are hallmarks of antigen exposure, followed by mutation & selection processes.

Immunoglobulin (IG) loci harbor inter-individual allelic variants in many different germline IG variable, diversity and joining genes of the IG heavy (IGH), kappa (IGK) and lambda (IGL) loci, which together form the genetic basis of the highly diverse antigen-specific B-cell receptors. These allelic variants can be shared between or be specific to human populations. The current immunogenetics resources gather the germline alleles, however, lack the population specificity of the alleles which poses limitations for disease-association studies related to immune responses in different human populations. Therefore, we systematically identified germline alleles from 26 different human populations around the world, profiled by "1000 Genomes" data. We identified 409 IGHV, 179 IGKV, and 199 IGLV germline alleles supported by at least seven haplotypes. The diversity of germline alleles is the highest in Africans. Remarkably, the variants in the identified novel alleles show strikingly conserved patterns, the same as found in other IG databases, suggesting over-time evolutionary selection processes. We could relate the genetic variants to population-specific immune responses, e.g. IGHV1-69 for flu in Africans. The population matched IG (pmIG) resource will enhance our understanding of the SHM-related B-cell receptor selection processes in (infectious) diseases and vaccination within and between different human populations., (© 2021. The Author(s).)

Acute myeloid leukemia (AML) is caused by genetic aberrations that also govern the prognosis of patients and guide risk-adapted and targeted therapy. Genetic aberrations in AML are structurally diverse and currently detected by different diagnostic assays. This study sought to establish whole transcriptome RNA sequencing as single, comprehensive, and flexible platform for AML diagnostics. We developed HAMLET (Human AML Expedited Transcriptomics) as bioinformatics pipeline for simultaneous detection of fusion genes, small variants, tandem duplications, and gene expression with all information assembled in an annotated, user-friendly output file. Whole transcriptome RNA sequencing was performed on 100 AML cases and HAMLET results were validated by reference assays and targeted resequencing. The data showed that HAMLET accurately detected all fusion genes and overexpression of EVI1 irrespective of 3q26 aberrations. In addition, small variants in 13 genes that are often mutated in AML were called with 99.2% sensitivity and 100% specificity, and tandem duplications in FLT3 and KMT2A were detected by a novel algorithm based on soft-clipped reads with 100% sensitivity and 97.1% specificity. In conclusion, HAMLET has the potential to provide accurate comprehensive diagnostic information relevant for AML classification, risk assessment and targeted therapy on a single technology platform.

Objective: To identify common genetic variants associated with the presence of brain microbleeds (BMBs)., Methods: We performed genome-wide association studies in 11 population-based cohort studies and 3 case-control or case-only stroke cohorts. Genotypes were imputed to the Haplotype Reference Consortium or 1000 Genomes reference panel. BMBs were rated on susceptibility-weighted or T2*-weighted gradient echo MRI sequences, and further classified as lobar or mixed (including strictly deep and infratentorial, possibly with lobar BMB). In a subset, we assessed the effects of APOE ε2 and ε4 alleles on BMB counts. We also related previously identified cerebral small vessel disease variants to BMBs., Results: BMBs were detected in 3,556 of the 25,862 participants, of which 2,179 were strictly lobar and 1,293 mixed. One locus in the APOE region reached genome-wide significance for its association with BMB (lead single nucleotide polymorphism rs769449; odds ratio [OR] any BMB [95% confidence interval (CI)] 1.33 [1.21-1.45]; p = 2.5 × 10 -10 ). APOE ε4 alleles were associated with strictly lobar (OR [95% CI] 1.34 [1.19-1.50]; p = 1.0 × 10 -6 ) but not with mixed BMB counts (OR [95% CI] 1.04 [0.86-1.25]; p = 0.68). APOE ε2 alleles did not show associations with BMB counts. Variants previously related to deep intracerebral hemorrhage and lacunar stroke, and a risk score of cerebral white matter hyperintensity variants, were associated with BMB., Conclusions: Genetic variants in the APOE region are associated with the presence of BMB, most likely due to the APOE ε4 allele count related to a higher number of strictly lobar BMBs. Genetic predisposition to small vessel disease confers risk of BMB, indicating genetic overlap with other cerebral small vessel disease markers., (Copyright © 2020 The Author(s). Published by Wolters Kluwer Health, Inc. on behalf of the American Academy of Neurology.)

Background: Financial codes are often used to extract diagnoses from electronic health records. This approach is prone to false positives. Alternatively, queries are constructed, but these are highly center and language specific. A tantalizing alternative is the automatic identification of patients by employing machine learning on format-free text entries., Objective: The aim of this study was to develop an easily implementable workflow that builds a machine learning algorithm capable of accurately identifying patients with rheumatoid arthritis from format-free text fields in electronic health records., Methods: Two electronic health record data sets were employed: Leiden (n=3000) and Erlangen (n=4771). Using a portion of the Leiden data (n=2000), we compared 6 different machine learning methods and a naïve word-matching algorithm using 10-fold cross-validation. Performances were compared using the area under the receiver operating characteristic curve (AUROC) and the area under the precision recall curve (AUPRC), and F1 score was used as the primary criterion for selecting the best method to build a classifying algorithm. We selected the optimal threshold of positive predictive value for case identification based on the output of the best method in the training data. This validation workflow was subsequently applied to a portion of the Erlangen data (n=4293). For testing, the best performing methods were applied to remaining data (Leiden n=1000; Erlangen n=478) for an unbiased evaluation., Results: For the Leiden data set, the word-matching algorithm demonstrated mixed performance (AUROC 0.90; AUPRC 0.33; F1 score 0.55), and 4 methods significantly outperformed word-matching, with support vector machines performing best (AUROC 0.98; AUPRC 0.88; F1 score 0.83). Applying this support vector machine classifier to the test data resulted in a similarly high performance (F1 score 0.81; positive predictive value [PPV] 0.94), and with this method, we could identify 2873 patients with rheumatoid arthritis in less than 7 seconds out of the complete collection of 23,300 patients in the Leiden electronic health record system. For the Erlangen data set, gradient boosting performed best (AUROC 0.94; AUPRC 0.85; F1 score 0.82) in the training set, and applied to the test data, resulted once again in good results (F1 score 0.67; PPV 0.97)., Conclusions: We demonstrate that machine learning methods can extract the records of patients with rheumatoid arthritis from electronic health record data with high precision, allowing research on very large populations for limited costs. Our approach is language and center independent and could be applied to any type of diagnosis. We have developed our pipeline into a universally applicable and easy-to-implement workflow to equip centers with their own high-performing algorithm. This allows the creation of observational studies of unprecedented size covering different countries for low cost from already available data in electronic health record systems., (©Tjardo D Maarseveen, Timo Meinderink, Marcel J T Reinders, Johannes Knitza, Tom W J Huizinga, Arnd Kleyer, David Simon, Erik B van den Akker, Rachel Knevel. Originally published in JMIR Medical Informatics (http://medinform.jmir.org), 30.11.2020.)

Background: Established rheumatoid arthritis (RA) patients display differentially expressed genes coding for cytokine/chemokine-mediated immunity compared to healthy controls. It is unclear, however, if in the pre-arthritis phase of clinically suspect arthralgia (CSA) expression of immune genes differ between patients who do or do not develop clinically evident inflammatory arthritis (IA)., Methods: Two hundred thirty-six consecutive patients presenting with arthralgia clinically suspected for progression to RA were followed until IA development or else for median 24 months (IQR 12-26). Baseline whole blood RNA expression was determined for a previously identified set of 133 genes associated with the innate and adaptive immune system by dual-color reverse-transcription multiplex ligation-dependent probe amplification (dcRT-MLPA) profiling. Cox proportional hazard models were used., Results: Twenty percent of CSA patients developed IA. After correction for multiple testing, expression levels of six genes (IFNG, PHEX, IGF-1, IL-7R, CD19, CCR7) at the time of presentation were associated with progression to IA. PHEX and IGF-1 were highly correlated with each other (ρ = 0.97). In multivariable analysis correcting for the different genes, expressions of IL-7R and IGF-1 were independently associated with IA development (p = 0.025, p = 0.046, respectively). Moreover, IL-7R and IGF-1 remained significantly associated even after correction for known predictors (ACPA, CRP, imaging-detected subclinical joint inflammation; p = 0.039, p = 0.005, respectively). These genes are also associated with RA development., Conclusions: IL-7R and IGF-1 were differentially expressed between CSA patients who did or did not progress to IA, independent from regularly used predictors. These biomarkers may become helpful in prognostication of CSA patients. Furthermore, because both genes are associated with T cell functioning, T cell dysregulation may mediate progression from arthralgia to arthritis.

Background: The blood metabolome incorporates cues from the environment and the host's genetic background, potentially offering a holistic view of an individual's health status., Methods: We have compiled a vast resource of proton nuclear magnetic resonance metabolomics and phenotypic data encompassing over 25 000 samples derived from 26 community and hospital-based cohorts., Results: Using this resource, we constructed a metabolomics-based age predictor (metaboAge) to calculate an individual's biological age. Exploration in independent cohorts demonstrates that being judged older by one's metabolome, as compared with one's chronological age, confers an increased risk on future cardiovascular disease, mortality, and functionality in older individuals. A web-based tool for calculating metaboAge (metaboage.researchlumc.nl) allows easy incorporation in other epidemiological studies. Access to data can be requested at bbmri.nl/samples-images-data., Conclusions: In summary, we present a vast resource of metabolomics data and illustrate its merit by constructing a metabolomics-based score for biological age that captures aspects of current and future cardiometabolic health.

Objective: To determine the small vessel disease spectrum associated with cysteine-altering NOTCH3 variants in community-dwelling individuals by analyzing the clinical and neuroimaging features of UK Biobank participants harboring such variants., Methods: The exome and genome sequencing datasets of the UK Biobank (n = 50,000) and cohorts of cognitively healthy elderly (n = 751) were queried for cysteine-altering NOTCH3 variants. Brain MRIs of individuals harboring such variants were scored according to Standards for Reporting Vascular Changes on Neuroimaging criteria, and clinical information was extracted with ICD-10 codes. Clinical and neuroimaging data were compared to age- and sex-matched UK Biobank controls and clinically diagnosed patients from the Dutch cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) registry., Results: We identified 108 individuals harboring a cysteine-altering NOTCH3 variant (2.2 of 1,000), of whom 75% have a variant that has previously been reported in CADASIL pedigrees. Almost all variants were located in 1 of the NOTCH3 protein epidermal growth factor-like repeat domains 7 to 34. White matter hyperintensity lesion load was higher in individuals with NOTCH3 variants than in controls ( p = 0.006) but lower than in patients with CADASIL with the same variants ( p < 0.001). Almost half of the 24 individuals with brain MRI had a Fazekas score of 0 or 1 up to age 70 years. There was no increased risk of stroke., Conclusions: Although community-dwelling individuals harboring a cysteine-altering NOTCH3 variant have a higher small vessel disease MRI burden than controls, almost half have no MRI abnormalities up to age 70 years. This shows that NOTCH3 cysteine altering variants are associated with an extremely broad phenotypic spectrum, ranging from CADASIL to nonpenetrance., (Copyright © 2020 The Author(s). Published by Wolters Kluwer Health, Inc. on behalf of the American Academy of Neurology.)