Proinflammatory phenotype of iPS cell-derived JAK2 V617F megakaryocytes induces fibrosis in 3D in vitro bone marrow niche.

The myeloproliferative disease polycythemia vera (PV) driven by the JAK2 V617F mutation can transform into myelofibrosis (post-PV-MF). It remains an open question how JAK2 V617F in hematopoietic stem cells induces MF. Megakaryocytes are major players in murine PV models but are difficult to study in the human setting. We generated induced pluripotent stem cells (iPSCs) from JAK2 V617F PV patients and differentiated them into megakaryocytes. In differentiation assays, JAK2 V617F iPSCs recapitulated the pathognomonic skewed megakaryocytic and erythroid differentiation. JAK2 V617F iPSCs had a TPO-independent and increased propensity to differentiate into megakaryocytes. RNA sequencing of JAK2 V617F iPSC-derived megakaryocytes reflected a proinflammatory, profibrotic phenotype and decreased ribosome biogenesis. In three-dimensional (3D) coculture, JAK2 V617F megakaryocytes induced a profibrotic phenotype through direct cell contact, which was reversed by the JAK2 inhibitor ruxolitinib. The 3D coculture system opens the perspective for further disease modeling and drug discovery.

Non-canonical Hedgehog signaling mediates profibrotic hematopoiesis-stroma crosstalk in myeloproliferative neoplasms.

The role of hematopoietic Hedgehog signaling in myeloproliferative neoplasms (MPNs) remains incompletely understood despite data suggesting that Hedgehog (Hh) pathway inhibitors have therapeutic activity in patients. We aim to systematically interrogate the role of canonical vs. non-canonical Hh signaling in MPNs. We show that Gli1 protein levels in patient peripheral blood mononuclear cells (PBMCs) mark fibrotic progression and that, in murine MPN models, absence of hematopoietic Gli1, but not Gli2 or Smo, significantly reduces MPN phenotype and fibrosis, indicating that GLI1 in the MPN clone can be activated in a non-canonical fashion. Additionally, we establish that hematopoietic Gli1 has a significant effect on stromal cells, mediated through a druggable MIF-CD74 axis. These data highlight the complex interplay between alterations in the MPN clone and activation of stromal cells and indicate that Gli1 represents a promising therapeutic target in MPNs, particularly that Hh signaling is dispensable for normal hematopoiesis.

EphrinA5 regulates cell motility by modulating Snhg15/DNA triplex-dependent targeting of DNMT1 to the Ncam1 promoter.

Cell-cell communication is mediated by membrane receptors and their ligands, such as the Eph/ephrin system, orchestrating cell migration during development and in diverse cancer types. Epigenetic mechanisms are key for integrating external "signals", e.g., from neighboring cells, into the transcriptome in health and disease. Previously, we reported ephrinA5 to trigger transcriptional changes of lncRNAs and protein-coding genes in cerebellar granule cells, a cell model for medulloblastoma. LncRNAs represent important adaptors for epigenetic writers through which they regulate gene expression. Here, we investigate a lncRNA-mediated targeting of DNMT1 to specific gene loci by the combined power of in silico modeling of RNA/DNA interactions and wet lab approaches, in the context of the clinically relevant use case of ephrinA5-dependent regulation of cellular motility of cerebellar granule cells. We provide evidence that Snhg15, a cancer-related lncRNA, recruits DNMT1 to the Ncam1 promoter through RNA/DNA triplex structure formation and the interaction with DNMT1. This mediates DNA methylation-dependent silencing of Ncam1, being abolished by ephrinA5 stimulation-triggered reduction of Snhg15 expression. Hence, we here propose a triple helix recognition mechanism, underlying cell motility regulation via lncRNA-targeted DNA methylation in a clinically relevant context.

Myelofibrosis at diagnosis is associated with the failure of treatment-free remission in CML patients.

The management of patients with chronic myeloid leukemia (CML) has been revolutionized by the introduction of tyrosine kinase inhibitors (TKIs), which induce deep molecular responses so that treatment can eventually be discontinued, leading to treatment-free remission (TFR) in a subset of patients. Unfortunately, leukemic stem cells (LSCs) often persist and a fraction of these can again expand in about half of patients that attempt TKI discontinuation. In this study, we show that presence of myelofibrosis (MF) at the time of diagnosis is a factor associating with TFR failure. Fibrotic transformation is governed by the action of several cytokines, and interestingly, some of them have also been described to support LSC persistence. At the cellular level, these could be produced by both malignant cells and by components of the bone marrow (BM) niche, including megakaryocytes (MKs) and mesenchymal stromal cells (MSCs). In our cohort of 57 patients, around 40% presented with MF at diagnosis and the number of blasts in the peripheral blood and BM was significantly elevated in patients with higher grade of MF. Employing a CML transgenic mouse model, we could observe higher levels of alpha-smooth muscle actin (α-SMA) in the BM when compared to control mice. Short-term treatment with the TKI nilotinib, efficiently reduced spleen weight and BCR::ABL1 mRNA levels, while α-SMA expression was only partially reduced. Interestingly, the number of MKs was increased in the spleen of CML mice and elevated in both BM and spleen upon nilotinib treatment. Analysis of human CML-vs healthy donor (HD)-derived MSCs showed an altered expression of gene signatures reflecting fibrosis as well as hematopoietic support, thus suggesting MSCs as a potential player in these two processes. Finally, in our cohort, 12 patients qualified for TKI discontinuation, and here we observed that all patients who failed TFR had BM fibrosis at diagnosis, whereas this was only the case in 25% of patients with achieved TFR, further supporting the link between fibrosis and LSC persistence.

KIT D816V Mast Cells Derived from Induced Pluripotent Stem Cells Recapitulate Systemic Mastocytosis Transcriptional Profile.

Mast cells (MCs) represent a population of hematopoietic cells with a key role in innate and adaptive immunity and are well known for their detrimental role in allergic responses. Yet, MCs occur in low abundance, which hampers their detailed molecular analysis. Here, we capitalized on the potential of induced pluripotent stem (iPS) cells to give rise to all cells in the body and established a novel and robust protocol for human iPS cell differentiation toward MCs. Relying on a panel of systemic mastocytosis (SM) patient-specific iPS cell lines carrying the KIT D816V mutation, we generated functional MCs that recapitulate SM disease features: increased number of MCs, abnormal maturation kinetics and activated phenotype, CD25 and CD30 surface expression and a transcriptional signature characterized by upregulated expression of innate and inflammatory response genes. Therefore, human iPS cell-derived MCs are a reliable, inexhaustible, and close-to-human tool for disease modeling and pharmacological screening to explore novel MC therapeutics.

Next-Generation Morphometry for pathomics-data mining in histopathology.

Pathology diagnostics relies on the assessment of morphology by trained experts, which remains subjective and qualitative. Here we developed a framework for large-scale histomorphometry (FLASH) performing deep learning-based semantic segmentation and subsequent large-scale extraction of interpretable, quantitative, morphometric features in non-tumour kidney histology. We use two internal and three external, multi-centre cohorts to analyse over 1000 kidney biopsies and nephrectomies. By associating morphometric features with clinical parameters, we confirm previous concepts and reveal unexpected relations. We show that the extracted features are independent predictors of long-term clinical outcomes in IgA-nephropathy. We introduce single-structure morphometric analysis by applying techniques from single-cell transcriptomics, identifying distinct glomerular populations and morphometric phenotypes along a trajectory of disease progression. Our study provides a concept for Next-generation Morphometry (NGM), enabling comprehensive quantitative pathology data mining, i.e., pathomics.

HIF1α-AS1 is a DNA:DNA:RNA triplex-forming lncRNA interacting with the HUSH complex.

DNA:DNA:RNA triplexes that are formed through Hoogsteen base-pairing of the RNA in the major groove of the DNA duplex have been observed in vitro, but the extent to which these interactions occur in cells and how they impact cellular functions remains elusive. Using a combination of bioinformatic techniques, RNA/DNA pulldown and biophysical studies, we set out to identify functionally important DNA:DNA:RNA triplex-forming long non-coding RNAs (lncRNA) in human endothelial cells. The lncRNA HIF1α-AS1 was retrieved as a top hit. Endogenous HIF1α-AS1 reduces the expression of numerous genes, including EPH Receptor A2 and Adrenomedullin through DNA:DNA:RNA triplex formation by acting as an adapter for the repressive human silencing hub complex (HUSH). Moreover, the oxygen-sensitive HIF1α-AS1 is down-regulated in pulmonary hypertension and loss-of-function approaches not only result in gene de-repression but also enhance angiogenic capacity. As exemplified here with HIF1α-AS1, DNA:DNA:RNA triplex formation is a functionally important mechanism of trans-acting gene expression control.

Induction of senescence upon loss of the Ash2l core subunit of H3K4 methyltransferase complexes.

Gene expression is controlled in part by post-translational modifications of core histones. Methylation of lysine 4 of histone H3 (H3K4), associated with open chromatin and gene transcription, is catalyzed by type 2 lysine methyltransferase complexes that require WDR5, RBBP5, ASH2L and DPY30 as core subunits. Ash2l is essential during embryogenesis and for maintaining adult tissues. To expand on the mechanistic understanding of Ash2l, we generated mouse embryo fibroblasts (MEFs) with conditional Ash2l alleles. Upon loss of Ash2l, methylation of H3K4 and gene expression were downregulated, which correlated with inhibition of proliferation and cell cycle progression. Moreover, we observed induction of senescence concomitant with a set of downregulated signature genes but independent of SASP. Many of the signature genes are FoxM1 responsive. Indeed, exogenous FOXM1 was sufficient to delay senescence. Thus, although the loss of Ash2l in MEFs has broad and complex consequences, a distinct set of downregulated genes promotes senescence.

LINC00152 Drives a Competing Endogenous RNA Network in Human Hepatocellular Carcinoma.

Genomic and epigenomic studies revealed dysregulation of long non-coding RNAs in many cancer entities, including liver cancer. We identified an epigenetic mechanism leading to upregulation of the long intergenic non-coding RNA 152 (LINC00152) expression in human hepatocellular carcinoma (HCC). Here, we aimed to characterize a potential competing endogenous RNA (ceRNA) network, in which LINC00152 exerts oncogenic functions by sponging miRNAs, thereby affecting their target gene expression. Database and gene expression data of human HCC were integrated to develop a potential LINC00152-driven ceRNA in silico. RNA immunoprecipitation and luciferase assay were used to identify miRNA binding to LINC00152 in human HCC cells. Functionally active players in the ceRNA network were analyzed using gene editing, siRNA or miRNA mimic transfection, and expression vectors in vitro. RNA expression in human HCC in vivo was validated by RNA in situ hybridization. Let-7c-5p, miR-23a-3p, miR-125a-5p, miR-125b-5p, miR-143a-3p, miR-193-3p, and miR-195-5p were detected as new components of the potential LINC00152 ceRNA network in human HCC. LINC00152 was confirmed to sponge miR143a-3p in human HCC cell lines, thereby limiting its binding to their respective target genes, like KLC2. KLC2 was identified as a central mediator promoting pro-tumorigenic effects of LINC00152 overexpression in HCC cells. Furthermore, co-expression of LINC00152 and KLC2 was observed in human HCC cohorts and high KLC2 expression was associated with shorter patient survival. Functional assays demonstrated that KLC2 promoted cell proliferation, clonogenicity and migration in vitro. The LINC00152-miR-143a-3p-KLC2 axis may represent a therapeutic target in human HCC.

Single-cell analysis of cultured bone marrow stromal cells reveals high similarity to fibroblasts in situ.

Within the heterogenous pool of bone marrow stromal cells, mesenchymal stromal cells (MSCs) are of particular interest because of their hematopoiesis-supporting capacities, contribution to disease progression, therapy resistance, and leukemic initiation. Cultured bone marrow-derived stromal cells (cBMSCs) are used for in vitro modeling of hematopoiesis-stroma interactions, validation of disease mechanisms, and screening for therapeutic targets. Here, we place cBMSCs (mouse and human) in a bone marrow tissue context by systematically comparing the transcriptome of plastic-adherent cells on a single-cell level with in vivo counterparts. Cultured BMSCs encompass a rather homogenous cell population, independent of the isolation method used and, although still possessing hematopoiesis-supporting capacity, are distinct from freshly isolated MSCs and more akin to in vivo fibroblast populations. Informed by combined cell trajectories and pathway analyses, we illustrate that TGFb inhibition in vitro can preserve a more "MSC"-like phenotype.

The spatial self-organization within pluripotent stem cell colonies is continued in detaching aggregates.

Colonies of induced pluripotent stem cells (iPSCs) reveal aspects of self-organization even under culture conditions that maintain pluripotency. To investigate the dynamics of this process under spatial confinement, we used either polydimethylsiloxane (PDMS) pillars or micro-contact printing of vitronectin. There was a progressive upregulation of OCT4, E-cadherin, and NANOG within 70 µm from the outer rim of iPSC colonies. Single-cell RNA-sequencing and spatial reconstruction of gene expression demonstrated that OCT4high subsets, residing at the edge of the colony, have pronounced up-regulation of the TGF-ß pathway, particularly of NODAL and its inhibitor LEFTY. Interestingly, after 5-7 days, iPSC colonies detached spontaneously from micro-contact printed substrates to form 3D aggregates. This new method allowed generation of embryoid bodies (EBs) of controlled size without enzymatic or mechanical treatment. Within the early 3D aggregates, radial organization and differential gene expression continued in analogy to the changes observed during self-organization of iPSC colonies. Early self-detached aggregates revealed up-regulated germline-specific gene expression patterns as compared to conventional EBs. However, there were no marked differences after further directed differentiation toward hematopoietic, mesenchymal, and neuronal lineages. Our results provide further insight into the gradual self-organization within iPSC colonies and at their transition into EBs.

Genetic barcoding systematically compares genes in del(5q) MDS and reveals a central role for CSNK1A1 in clonal expansion.

How genetic haploinsufficiency contributes to the clonal dominance of hematopoietic stem cells (HSCs) in del(5q) myelodysplastic syndrome (MDS) remains unresolved. Using a genetic barcoding strategy, we performed a systematic comparison on genes implicated in the pathogenesis of del(5q) MDS in direct competition with each other and wild-type (WT) cells with single-clone resolution. Csnk1a1 haploinsufficient HSCs expanded (oligo)clonally and outcompeted all other tested genes and combinations. Csnk1a1-/+ multipotent progenitors showed a proproliferative gene signature and HSCs showed a downregulation of inflammatory signaling/immune response. In validation experiments, Csnk1a1-/+ HSCs outperformed their WT counterparts under a chronic inflammation stimulus, also known to be caused by neighboring genes on chromosome 5. We therefore propose a crucial role for Csnk1a1 haploinsufficiency in the selective advantage of 5q-HSCs, implemented by creation of a unique competitive advantage through increased HSC self-renewal and proliferation capacity, as well as increased fitness under inflammatory stress.

Cyclin E1 in Murine and Human Liver Cancer: A Promising Target for Therapeutic Intervention during Tumour Progression.

Cyclin E1 (CCNE1) is a regulatory subunit of Cyclin-dependent kinase 2 (CDK2) and is thought to control the transition of quiescent cells into cell cycle progression. Recently, we identified CCNE1 and CDK2 as key factors for the initiation of hepatocellular carcinoma (HCC). In the present study, we dissected the contributions of CCNE1 and CDK2 for HCC progression in mice and patients. Therefore, we generated genetically modified mice allowing inducible deletion of Ccne1 or Cdk2. After initiation of HCC, using the hepatocarcinogen diethylnitrosamine (DEN), we deleted Ccne1 or Cdk2 and subsequently analysed HCC progression. The relevance of CCNE1 or CDK2 for human HCC progression was investigated by in silico database analysis. Interventional deletion of Ccne1, but not of Cdk2, substantially reduced the HCC burden in mice. Ccne1-deficient HCCs were characterised by attenuated proliferation, impaired DNA damage response and downregulation of markers for stemness and microinvasion. Additionally, the tumour microenvironment of Ccne1-deficient mice showed a reduction in immune mediators, myeloid cells and cancer-associated fibroblasts. In sharp contrast, Cdk2 was dispensable for HCC progression in mice. In agreement with our mouse data, CCNE1 was overexpressed in HCC patients independent of risk factors, and associated with reduced disease-free survival, a common signature for enhanced chromosomal instability, proliferation, dedifferentiation and invasion. However, CDK2 lacked diagnostic or prognostic value in HCC patients. In summary, CCNE1 drives HCC progression in a CDK2-independent manner in mice and man. Therefore, interventional inactivation of CCNE1 represents a promising strategy the treatment of liver cancer.

CALR frameshift mutations in MPN patient-derived iPSCs accelerate maturation of megakaryocytes.

Calreticulin (CALR) mutations are driver mutations in myeloproliferative neoplasms (MPNs), leading to activation of the thrombopoietin receptor and causing abnormal megakaryopoiesis. Here, we generated patient-derived CALRins5- or CALRdel52-positive induced pluripotent stem cells (iPSCs) to establish an MPN disease model for molecular and mechanistic studies. We demonstrated myeloperoxidase deficiency in granulocytic cells derived from homozygous CALR mutant iPSCs, rescued by repairing the mutation using CRISPR/Cas9. iPSC-derived megakaryocytes showed characteristics of primary megakaryocytes such as formation of demarcation membrane system and cytoplasmic pro-platelet protrusions. Importantly, CALR mutations led to enhanced megakaryopoiesis and accelerated megakaryocytic development in a thrombopoietin-independent manner. Mechanistically, our study identified differentially regulated pathways in mutated versus unmutated megakaryocytes, such as hypoxia signaling, which represents a potential target for therapeutic intervention. Altogether, we demonstrate key aspects of mutated CALR-driven pathogenesis dependent on its zygosity, and found novel therapeutic targets, making our model a valuable tool for clinical drug screening in MPNs.

From cell to cell: Identification of actionable targets in bone marrow fibrosis using single-cell technologies.

Single-cell technologies have rapidly developed in recent years and have already had a significant impact on the research of myeloproliferative neoplasms. The increasing number of publicly available data sets allows characterization of the bone marrow niche in patients and mouse models at unprecedented resolution. Single-cell RNA sequencing has successfully been used to identify and characterize disease-driving cell populations and to identify the alarmin S100A8/A9 as an important mediator of myelofibrosis and potent therapeutic target. It is now possible to execute a streamlined set of experiments to specifically identify and validate actionable target genes functionally with the advance of reliable in vivo models and the possibility of conducting single-cell analyses with a minimal amount of patient material. The advent of large-scale analyses of both hematopoietic and non-hematopoietic bone marrow cells will allow comprehensive network analyses guiding an increasingly detailed mapping of the MPN interactome.

The Expression of the Cancer-Associated lncRNA Snhg15 Is Modulated by EphrinA5-Induced Signaling.

The Eph receptor tyrosine kinases and their respective ephrin-ligands are an important family of membrane receptors, being involved in developmental processes such as proliferation, migration, and in the formation of brain cancer such as glioma. Intracellular signaling pathways, which are activated by Eph receptor signaling, are well characterized. In contrast, it is unknown so far whether ephrins modulate the expression of lncRNAs, which would enable the transduction of environmental stimuli into our genome through a great gene regulatory spectrum. Applying a combination of functional in vitro assays, RNA sequencing, and qPCR analysis, we found that the proliferation and migration promoting stimulation of mouse cerebellar granule cells (CB) with ephrinA5 diminishes the expression of the cancer-related lncRNA Snhg15. In a human medulloblastoma cell line (DAOY) ephrinA5 stimulation similarly reduced SNHG15 expression. Computational analysis identified triple-helix-mediated DNA-binding sites of Snhg15 in promoters of genes found up-regulated upon ephrinA5 stimulation and known to be involved in tumorigenic processes. Our findings propose a crucial role of Snhg15 downstream of ephrinA5-induced signaling in regulating gene transcription in the nucleus. These findings could be potentially relevant for the regulation of tumorigenic processes in the context of glioma.

Endothelium-derived stromal cells contribute to hematopoietic bone marrow niche formation.

Bone marrow stromal cells (BMSCs) play pivotal roles in tissue maintenance and regeneration. Their origins, however, remain incompletely understood. Here we identify rare LNGFR+ cells in human fetal and regenerative bone marrow that co-express endothelial and stromal markers. This endothelial subpopulation displays transcriptional reprogramming consistent with endothelial-to-mesenchymal transition (EndoMT) and can generate multipotent stromal cells that reconstitute the bone marrow (BM) niche upon transplantation. Single-cell transcriptomics and lineage tracing in mice confirm robust and sustained contributions of EndoMT to bone precursor and hematopoietic niche pools. Interleukin-33 (IL-33) is overexpressed in subsets of EndoMT cells and drives this conversion process through ST2 receptor signaling. These data reveal generation of tissue-forming BMSCs from mouse and human endothelial cells and may be instructive for approaches to human tissue regeneration.

Deep learning-based clustering approaches for bioinformatics.

Clustering is central to many data-driven bioinformatics research and serves a powerful computational method. In particular, clustering helps at analyzing unstructured and high-dimensional data in the form of sequences, expressions, texts and images. Further, clustering is used to gain insights into biological processes in the genomics level, e.g. clustering of gene expressions provides insights on the natural structure inherent in the data, understanding gene functions, cellular processes, subtypes of cells and understanding gene regulations. Subsequently, clustering approaches, including hierarchical, centroid-based, distribution-based, density-based and self-organizing maps, have long been studied and used in classical machine learning settings. In contrast, deep learning (DL)-based representation and feature learning for clustering have not been reviewed and employed extensively. Since the quality of clustering is not only dependent on the distribution of data points but also on the learned representation, deep neural networks can be effective means to transform mappings from a high-dimensional data space into a lower-dimensional feature space, leading to improved clustering results. In this paper, we review state-of-the-art DL-based approaches for cluster analysis that are based on representation learning, which we hope to be useful, particularly for bioinformatics research. Further, we explore in detail the training procedures of DL-based clustering algorithms, point out different clustering quality metrics and evaluate several DL-based approaches on three bioinformatics use cases, including bioimaging, cancer genomics and biomedical text mining. We believe this review and the evaluation results will provide valuable insights and serve a starting point for researchers wanting to apply DL-based unsupervised methods to solve emerging bioinformatics research problems.

CrossTalkeR: analysis and visualization of ligand-receptorne tworks.

MOTIVATION: Ligand-receptor (LR) network analysis allows the characterization of cellular crosstalk based on single cell RNA-seq data. However, current methods typically provide a list of inferred LR interactions and do not allow the researcher to focus on specific cell types, ligands or receptors. In addition, most of these methods cannot quantify changes in crosstalk between two biological phenotypes. RESULTS: CrossTalkeR is a framework for network analysis and visualization of LR interactions. CrossTalkeR identifies relevant ligands, receptors and cell types contributing to changes in cell communication when contrasting two biological phenotypes, i.e. disease versus homeostasis. A case study on scRNA-seq of human myeloproliferative neoplasms reinforces the strengths of CrossTalkeR for characterization of changes in cellular crosstalk in disease. AVAILABILITY AND IMPLEMENTATION: CrosstalkeR is an R package available at: Github: https://github.com/CostaLab/CrossTalkeR. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.