Serum proteomics reveals hemophagocytic lymphohistiocytosis-like phenotype in a subset of patients with multisystem inflammatory syndrome in children.

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.

A transcriptomic based deconvolution framework for assessing differentiation stages and drug responses of AML.

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.

Expanding the repertoire reveals recurrent, cryptic, and hematopoietic HLA class I minor histocompatibility antigens.

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.

Mapping AML heterogeneity - multi-cohort transcriptomic analysis identifies novel clusters and divergent ex-vivo drug responses.

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.

Quantitative multiple fragment monitoring with enhanced in-source fragmentation/annotation mass spectrometry.

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.

Single-cell immune profiling reveals thymus-seeding populations, T cell commitment, and multilineage development in the human thymus.

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.

Metabolomic profiles predict individual multidisease outcomes.

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.

Clonal Hematopoiesis Analyses in Clinical, Epidemiologic, and Genetic Aging Studies to Unravel Underlying Mechanisms of Age-Related Dysfunction in Humans.

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.

T and NK Cells in IL2RG-Deficient Patient 50 Years After Hematopoietic Stem Cell Transplantation.

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

A framework for employing longitudinally collected multicenter electronic health records to stratify heterogeneous patient populations on disease history.

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.

Comprehensive diagnostics of acute myeloid leukemia by whole transcriptome RNA sequencing.

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.

Broad phenotype of cysteine-altering NOTCH3 variants in UK Biobank: CADASIL to nonpenetrance.

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.

Human Plasma N-glycosylation as Analyzed by Matrix-Assisted Laser Desorption/Ionization-Fourier Transform Ion Cyclotron Resonance-MS Associates with Markers of Inflammation and Metabolic Health.

Glycosylation is an abundant co- and post-translational protein modification of importance to protein processing and activity. Although not template-defined, glycosylation does reflect the biological state of an organism and is a high-potential biomarker for disease and patient stratification. However, to interpret a complex but informative sample like the total plasma N-glycome, it is important to establish its baseline association with plasma protein levels and systemic processes. Thus far, large-scale studies (n >200) of the total plasma N-glycome have been performed with methods of chromatographic and electrophoretic separation, which, although being informative, are limited in resolving the structural complexity of plasma N-glycans. MS has the opportunity to contribute additional information on, among others, antennarity, sialylation, and the identity of high-mannose type species.Here, we have used matrix-assisted laser desorption/ionization (MALDI)-Fourier transform ion cyclotron resonance (FTICR)-MS to study the total plasma N-glycome of 2144 healthy middle-aged individuals from the Leiden Longevity Study, to allow association analysis with markers of metabolic health and inflammation. To achieve this, N-glycans were enzymatically released from their protein backbones, labeled at the reducing end with 2-aminobenzoic acid, and following purification analyzed by negative ion mode intermediate pressure MALDI-FTICR-MS. In doing so, we achieved the relative quantification of 61 glycan compositions, ranging from Hex4HexNAc2 to Hex7HexNAc6dHex1Neu5Ac4, as well as that of 39 glycosylation traits derived thereof. Next to confirming known associations of glycosylation with age and sex by MALDI-FTICR-MS, we report novel associations with C-reactive protein (CRP), interleukin 6 (IL-6), body mass index (BMI), leptin, adiponectin, HDL cholesterol, triglycerides (TG), insulin, gamma-glutamyl transferase (GGT), alanine aminotransferase (ALT), and smoking. Overall, the bisection, galactosylation, and sialylation of diantennary species, the sialylation of tetraantennary species, and the size of high-mannose species proved to be important plasma characteristics associated with inflammation and metabolic health.

Fibroblast growth factor 21 reflects liver fat accumulation and dysregulation of signalling pathways in the liver of C57BL/6J mice.

Fibroblast growth factor 21 (Fgf21) has emerged as a potential plasma marker to diagnose non-alcoholic fatty liver disease (NAFLD). To study the molecular processes underlying the association of plasma Fgf21 with NAFLD, we explored the liver transcriptome data of a mild NAFLD model of aging C57BL/6J mice at 12, 24, and 28 months of age. The plasma Fgf21 level significantly correlated with intrahepatic triglyceride content. At the molecular level, elevated plasma Fgf21 levels were associated with dysregulated metabolic and cancer-related pathways. The up-regulated Fgf21 levels in NAFLD were implied to be a protective response against the NAFLD-induced adverse effects, e.g. lipotoxicity, oxidative stress and endoplasmic reticulum stress. An in vivo PPARα challenge demonstrated the dysregulation of PPARα signalling in the presence of NAFLD, which resulted in a stochastically increasing hepatic expression of Fgf21. Notably, elevated plasma Fgf21 was associated with declining expression of Klb, Fgf21's crucial co-receptor, which suggests a resistance to Fgf21. Therefore, although liver fat accumulation is a benign stage of NAFLD, the elevated plasma Fgf21 likely indicated vulnerability to metabolic stressors that may contribute towards progression to end-stage NAFLD. In conclusion, plasma levels of Fgf21 reflect liver fat accumulation and dysregulation of metabolic pathways in the liver.

Metabolic effects of a 13-weeks lifestyle intervention in older adults: The Growing Old Together Study.

For people in their 40s and 50s, lifestyle programs have been shown to improve metabolic health. For older adults, however, it is not clear whether these programs are equally healthy. In the Growing Old Together study, we applied a 13-weeks lifestyle program, with a target of 12.5% caloric restriction and 12.5% increase in energy expenditure through an increase in physical activity, in 164 older adults (mean age=63.2 years; BMI=23-35 kg/m2). Mean weight loss was 4.2% (SE=2.8%) of baseline weight, which is comparable to a previous study in younger adults. Fasting insulin levels, however, showed a much smaller decrease (0.30 mU/L (SE=3.21)) and a more heterogeneous response (range=2.0-29.6 mU/L). Many other parameters of metabolic health, such as blood pressure, and thyroid, glucose and lipid metabolism improved significantly. Many 1H-NMR metabolites changed in a direction previously associated with a low risk of type 2 diabetes and cardiovascular disease and partially independently of weight loss. In conclusion, 25% reduction in energy balance for 13 weeks induced a metabolic health benefit in older adults, monitored by traditional and novel metabolic markers.

Exome and whole genome sequencing in aging and longevity.

Calendar age is the major risk factor for common disease. It is therefore expected that understanding the aging process will eventually lead to promotion of better health conditions in elderly populations. Such insight may be obtained by identifying the genetic determinants of familial and exceptional longevity and age-related disease. Research of these determinants has been performed in candidate gene, genome-wide association and linkage studies. Because exploration of the common variation in the genome did not explain much of the variation in the rate of aging and longevity, researchers in the field have only recently started to investigate the contribution of rare genetic variants to these traits. The increased throughput and decreased costs of next generation sequencing (NGS) have resulted in a great deal of novel applications for sequencing sets of candidate genes, whole exomes, and whole genomes of individuals. Most of the successful NGS applications are as yet those focused on genetic syndromes and cancers for which causal mutations are readily being identified. In this chapter, we discuss the genetic and phenomic aspects of human aging research and the use of NGS data to identify genes relevant for age-related diseases and lifespan regulation, and to investigate the accumulation of somatic genetic variation during the course of life.

Integrating protein-protein interaction networks with gene-gene co-expression networks improves gene signatures for classifying breast cancer metastasis.

Multiple studies have illustrated that gene expression profiling of primary breast cancers throughout the final stages of tumor development can provide valuable markers for risk prediction of metastasis and disease sub typing. However, the identification of a biologically interpretable and universally shared set of markers proved to be difficult. Here, we propose a method for de novo grouping of genes by dissecting the protein-protein interaction network into disjoint sub networks using pair wise gene expression correlation measures. We show that the obtained sub networks are functionally coherent and are consistently identified when applied on a compendium composed of six different breast cancer studies. Application of the proposed method using different integration approaches underlines the robustness of the identified sub network related to cell cycle and identifies putative new sub network markers for metastasis related to cell-cell adhesion, the proteasome complex and JUN-FOS signalling. Although gene selection with the proposed method does not directly improve upon previously reported cross study classification performances, it shows great promises for applications in data integration and result interpretation.