BloodSpot 3.0: a database of gene and protein expression data in normal and malignant haematopoiesis.

BloodSpot is a specialised database integrating gene expression data from acute myeloid leukaemia (AML) patients related to blood cell development and maturation. The database and interface has helped numerous researchers and clinicians to quickly get an overview of gene expression patterns in healthy and malignant haematopoiesis. Here, we present an update to our framework that includes protein expression data of sorted single cells. With this update we also introduce datasets broadly spanning age groups, which many users have requested, with particular interest for researchers studying paediatric leukaemias. The backend of the database has been rewritten and migrated to a cloud-based environment to accommodate the growth, and provide a better user-experience for our many international users. Users can now enjoy faster transfer speeds and a more responsive interface. In conclusion, the continuing popularity of the database and emergence of new data modalities has prompted us to rewrite and futureproof the back-end, including paediatric centric views, as well as single cell protein data, allowing us to keep the database updated and relevant for the years to come. The database is freely available at www.bloodspot.eu.

PTBP1 promotes hematopoietic stem cell maintenance and red blood cell development by ensuring sufficient availability of ribosomal constituents.

Ribosomopathies constitute a range of disorders associated with defective protein synthesis mainly affecting hematopoietic stem cells (HSCs) and erythroid development. Here, we demonstrate that deletion of poly-pyrimidine-tract-binding protein 1 (PTBP1) in the hematopoietic compartment leads to the development of a ribosomopathy-like condition. Specifically, loss of PTBP1 is associated with decreases in HSC self-renewal, erythroid differentiation, and protein synthesis. Consistent with its function as a splicing regulator, PTBP1 deficiency results in splicing defects in hundreds of genes, and we demonstrate that the up-regulation of a specific isoform of CDC42 partly mimics the protein-synthesis defect associated with loss of PTBP1. Furthermore, PTBP1 deficiency is associated with a marked defect in ribosome biogenesis and a selective reduction in the translation of mRNAs encoding ribosomal proteins. Collectively, this work identifies PTBP1 as a key integrator of ribosomal functions and highlights the broad functional repertoire of RNA-binding proteins.

Quantitative single-cell proteomics as a tool to characterize cellular hierarchies.

Large-scale single-cell analyses are of fundamental importance in order to capture biological heterogeneity within complex cell systems, but have largely been limited to RNA-based technologies. Here we present a comprehensive benchmarked experimental and computational workflow, which establishes global single-cell mass spectrometry-based proteomics as a tool for large-scale single-cell analyses. By exploiting a primary leukemia model system, we demonstrate both through pre-enrichment of cell populations and through a non-enriched unbiased approach that our workflow enables the exploration of cellular heterogeneity within this aberrant developmental hierarchy. Our approach is capable of consistently quantifying ~1000 proteins per cell across thousands of individual cells using limited instrument time. Furthermore, we develop a computational workflow (SCeptre) that effectively normalizes the data, integrates available FACS data and facilitates downstream analysis. The approach presented here lays a foundation for implementing global single-cell proteomics studies across the world.

ERG Controls B Cell Development by Promoting Igh V-to-DJ Recombination.

B cell development depends on the coordinated expression and cooperation of several transcription factors. Here we show that the transcription factor ETS-related gene (ERG) is crucial for normal B cell development and that its deletion results in a substantial loss of bone marrow B cell progenitors and peripheral B cells, as well as a skewing of splenic B cell populations. We find that ERG-deficient B lineage cells exhibit an early developmental block at the pre-B cell stage and proliferate less. The cells fail to express the immunoglobulin heavy chain due to inefficient V-to-DJ recombination, and cells that undergo recombination display a strong bias against incorporation of distal V gene segments. Furthermore, antisense transcription at PAX5-activated intergenic repeat (PAIR) elements, located in the distal region of the Igh locus, depends on ERG. These findings show that ERG serves as a critical regulator of B cell development by ensuring efficient and balanced V-to-DJ recombination.

Oral vitamin C supplementation to patients with myeloid cancer on azacitidine treatment: Normalization of plasma vitamin C induces epigenetic changes.

BACKGROUND: Patients with haematological malignancies are often vitamin C deficient, and vitamin C is essential for the TET-induced conversion of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC), the first step in active DNA demethylation. Here, we investigate whether oral vitamin C supplementation can correct vitamin C deficiency and affect the 5hmC/5mC ratio in patients with myeloid cancers treated with DNA methyltransferase inhibitors (DNMTis). RESULTS: We conducted a randomized, double-blinded, placebo-controlled pilot trial (NCT02877277) in Danish patients with myeloid cancers performed during 3 cycles of DNMTi-treatment (5-azacytidine, 100 mg/m2/d for 5 days in 28-day cycles) supplemented by oral dose of 500 mg vitamin C (n = 10) or placebo (n = 10) daily during the last 2 cycles. Fourteen patients (70%) were deficient in plasma vitamin C (< 23 µM) and four of the remaining six patients were taking vitamin supplements at inclusion. Global DNA methylation was significantly higher in patients with severe vitamin C deficiency (< 11.4 µM; 4.997 vs 4.656% 5mC relative to deoxyguanosine, 95% CI [0.126, 0.556], P = 0.004). Oral supplementation restored plasma vitamin C levels to the normal range in all patients in the vitamin C arm (mean increase 34.85 ± 7.94 µM, P = 0.0004). We show for the first time that global 5hmC/5mC levels were significantly increased in mononuclear myeloid cells from patients receiving oral vitamin C compared to placebo (0.037% vs - 0.029%, 95% CI [- 0.129, - 0.003], P = 0.041). CONCLUSIONS: Normalization of plasma vitamin C by oral supplementation leads to an increase in the 5hmC/5mC ratio compared to placebo-treated patients and may enhance the biological effects of DNMTis. The clinical efficacy of oral vitamin C supplementation to DNMTis should be investigated in a large randomized, placebo-controlled clinical trial. TRIAL REGISTRATION: ClinicalTrials.gov, NCT02877277 . Registered on 9 August 2016, retrospectively registered.

Mutant CEBPA directly drives the expression of the targetable tumor-promoting factor CD73 in AML.

The key myeloid transcription factor (TF), CEBPA, is frequently mutated in acute myeloid leukemia (AML), but the direct molecular effects of this leukemic driver mutation remain elusive. To investigate CEBPA mutant AML, we performed microscale, in vivo chromatin immunoprecipitation sequencing and identified a set of aberrantly activated enhancers, exclusively occupied by the leukemia-associated CEBPA-p30 isoform. Comparing gene expression changes in human CEBPA mutant AML and the corresponding Cebpa Lp30 mouse model, we identified Nt5e, encoding CD73, as a cross-species AML gene with an upstream leukemic enhancer physically and functionally linked to the gene. Increased expression of CD73, mediated by the CEBPA-p30 isoform, sustained leukemic growth via the CD73/A2AR axis. Notably, targeting of this pathway enhanced survival of AML-transplanted mice. Our data thus indicate a first-in-class link between a cancer driver mutation in a TF and a druggable, direct transcriptional target.

The splicing factor RBM25 controls MYC activity in acute myeloid leukemia.

Cancer sequencing studies have implicated regulators of pre-mRNA splicing as important disease determinants in acute myeloid leukemia (AML), but the underlying mechanisms have remained elusive. We hypothesized that "non-mutated" splicing regulators may also play a role in AML biology and therefore conducted an in vivo shRNA screen in a mouse model of CEBPA mutant AML. This has led to the identification of the splicing regulator RBM25 as a novel tumor suppressor. In multiple human leukemic cell lines, knockdown of RBM25 promotes proliferation and decreases apoptosis. Mechanistically, we show that RBM25 controls the splicing of key genes, including those encoding the apoptotic regulator BCL-X and the MYC inhibitor BIN1. This mechanism is also operative in human AML patients where low RBM25 levels are associated with high MYC activity and poor outcome. Thus, we demonstrate that RBM25 acts as a regulator of MYC activity and sensitizes cells to increased MYC levels.

BloodSpot: a database of healthy and malignant haematopoiesis updated with purified and single cell mRNA sequencing profiles.

BloodSpot is a gene-centric database of mRNA expression of haematopoietic cells. The web-based interface to the database includes three concomitant levels of visualization for a gene query; foremost is the expression across hematopoietic cell types, second is analysis of survival of Acute Myeloid Leukaemia patients based on gene expression, and lastly, the expression visualized in an interactive developmental tree. With the introduction of single cell data we have now also included an unbiased dimensionality reduction method to show gene expression over the continuum of haematopoiesis. The webserver includes a few select analysis functionalities, like Student's t-test, identification of correlating genes and lookup of whole genetic signatures, with the aim of making generation and testing of hypotheses quick and intuitive. The visualizations have been updated to accommodate new datatypes and the database has been largely expanded with RNA-sequencing datasets, both purified in bulk and at single cell resolution, increasing the number of single samples more than 10 fold, while keeping simplicity in presentation. The database should be of interest for any researcher within leukaemia, haematopoiesis, cellular development, or stem cells. The database is freely available at www.bloodspot.eu.

Differences in Cell Cycle Status Underlie Transcriptional Heterogeneity in the HSC Compartment.

Hematopoietic stem cells (HSCs) are considered a heterogeneous cell population. To further resolve the HSC compartment, we characterized a retinoic acid (RA) reporter mouse line. Sub-fractionation of the HSC compartment in RA-CFP reporter mice demonstrated that RA-CFP-dim HSCs were largely non-proliferative and displayed superior engraftment potential in comparison with RA-CFP-bright HSCs. Gene expression analysis demonstrated higher expression of RA-target genes in RA-CFP-dim HSCs, in contrast to the RA-CFP reporter expression, but both RA-CFP-dim and RA-CFP-bright HSCs responded efficiently to RA in vitro. Single-cell RNA sequencing (RNA-seq) of >1,200 HSCs showed that differences in cell cycle activity constituted the main driver of transcriptional heterogeneity in HSCs. Moreover, further analysis of the single-cell RNA-seq data revealed that stochastic low-level expression of distinct lineage-affiliated transcriptional programs is a common feature of HSCs. Collectively, this work demonstrates the utility of the RA-CFP reporter line as a tool for the isolation of superior HSCs.

Leukemogenic nucleophosmin mutation disrupts the transcription factor hub that regulates granulomonocytic fates.

Nucleophosmin (NPM1) is among the most frequently mutated genes in acute myeloid leukemia (AML). It is not known, however, how the resulting oncoprotein mutant NPM1 is leukemogenic. To reveal the cellular machinery in which NPM1 participates in myeloid cells, we analyzed the endogenous NPM1 protein interactome by mass spectrometry and discovered abundant amounts of the master transcription factor driver of monocyte lineage differentiation PU.1 (also known as SPI1). Mutant NPM1, which aberrantly accumulates in cytoplasm, dislocated PU.1 into cytoplasm with it. CEBPA and RUNX1, the master transcription factors that collaborate with PU.1 to activate granulomonocytic lineage fates, remained nuclear; but without PU.1, their coregulator interactions were toggled from coactivators to corepressors, repressing instead of activating more than 500 granulocyte and monocyte terminal differentiation genes. An inhibitor of nuclear export, selinexor, by locking mutant NPM1/PU.1 in the nucleus, activated terminal monocytic fates. Direct depletion of the corepressor DNA methyltransferase 1 (DNMT1) from the CEBPA/RUNX1 protein interactome using the clinical drug decitabine activated terminal granulocytic fates. Together, these noncytotoxic treatments extended survival by more than 160 days versus vehicle in a patient-derived xenotransplant model of NPM1/FLT3-mutated AML. In sum, mutant NPM1 represses monocyte and granulocyte terminal differentiation by disrupting PU.1/CEBPA/RUNX1 collaboration, a transforming action that can be reversed by pharmacodynamically directed dosing of clinical small molecules.

Enhancer and Transcription Factor Dynamics during Myeloid Differentiation Reveal an Early Differentiation Block in Cebpa null Progenitors.

Transcription factors PU.1 and CEBPA are required for the proper coordination of enhancer activity during granulocytic-monocytic (GM) lineage differentiation to form myeloid cells. However, precisely how these factors control the chronology of enhancer establishment during differentiation is not known. Through integrated analyses of enhancer dynamics, transcription factor binding, and proximal gene expression during successive stages of murine GM-lineage differentiation, we unravel the distinct kinetics by which PU.1 and CEBPA coordinate GM enhancer activity. We find no evidence of a pioneering function of PU.1 during late GM-lineage differentiation. Instead, we delineate a set of enhancers that gain accessibility in a CEBPA-dependent manner, suggesting a pioneering function of CEBPA. Analyses of Cebpa null bone marrow demonstrate that CEBPA controls PU.1 levels and, unexpectedly, that the loss of CEBPA results in an early differentiation block. Taken together, our data provide insights into how PU.1 and CEBPA functionally interact to drive GM-lineage differentiation.

Characterization of the enhancer and promoter landscape of inflammatory bowel disease from human colon biopsies.

Inflammatory bowel disease (IBD) is a chronic intestinal disorder, with two main types: Crohn's disease (CD) and ulcerative colitis (UC), whose molecular pathology is not well understood. The majority of IBD-associated SNPs are located in non-coding regions and are hard to characterize since regulatory regions in IBD are not known. Here we profile transcription start sites (TSSs) and enhancers in the descending colon of 94 IBD patients and controls. IBD-upregulated promoters and enhancers are highly enriched for IBD-associated SNPs and are bound by the same transcription factors. IBD-specific TSSs are associated to genes with roles in both inflammatory cascades and gut epithelia while TSSs distinguishing UC and CD are associated to gut epithelia functions. We find that as few as 35 TSSs can distinguish active CD, UC, and controls with 85% accuracy in an independent cohort. Our data constitute a foundation for understanding the molecular pathology, gene regulation, and genetics of IBD.

EZH2 is a potential therapeutic target for H3K27M-mutant pediatric gliomas.

Diffuse intrinsic pontine glioma (DIPG) is an aggressive brain tumor that is located in the pons and primarily affects children. Nearly 80% of DIPGs harbor mutations in histone H3 genes, wherein lysine 27 is substituted with methionine (H3K27M). H3K27M has been shown to inhibit polycomb repressive complex 2 (PRC2), a multiprotein complex responsible for the methylation of H3 at lysine 27 (H3K27me), by binding to its catalytic subunit EZH2. Although DIPGs with the H3K27M mutation show global loss of H3K27me3, several genes retain H3K27me3. Here we describe a mouse model of DIPG in which H3K27M potentiates tumorigenesis. Using this model and primary patient-derived DIPG cell lines, we show that H3K27M-expressing tumors require PRC2 for proliferation. Furthermore, we demonstrate that small-molecule EZH2 inhibitors abolish tumor cell growth through a mechanism that is dependent on the induction of the tumor-suppressor protein p16INK4A. Genome-wide enrichment analyses show that the genes that retain H3K27me3 in H3K27M cells are strong polycomb targets. Furthermore, we find a highly significant overlap between genes that retain H3K27me3 in the DIPG mouse model and in human primary DIPGs expressing H3K27M. Taken together, these results show that residual PRC2 activity is required for the proliferation of H3K27M-expressing DIPGs, and that inhibition of EZH2 is a potential therapeutic strategy for the treatment of these tumors.

Genome-wide analysis of differential transcriptional and epigenetic variability across human immune cell types.

BACKGROUND: A healthy immune system requires immune cells that adapt rapidly to environmental challenges. This phenotypic plasticity can be mediated by transcriptional and epigenetic variability. RESULTS: We apply a novel analytical approach to measure and compare transcriptional and epigenetic variability genome-wide across CD14+CD16- monocytes, CD66b+CD16+ neutrophils, and CD4+CD45RA+ naïve T cells from the same 125 healthy individuals. We discover substantially increased variability in neutrophils compared to monocytes and T cells. In neutrophils, genes with hypervariable expression are found to be implicated in key immune pathways and are associated with cellular properties and environmental exposure. We also observe increased sex-specific gene expression differences in neutrophils. Neutrophil-specific DNA methylation hypervariable sites are enriched at dynamic chromatin regions and active enhancers. CONCLUSIONS: Our data highlight the importance of transcriptional and epigenetic variability for the key role of neutrophils as the first responders to inflammatory stimuli. We provide a resource to enable further functional studies into the plasticity of immune cells, which can be accessed from: http://blueprint-dev.bioinfo.cnio.es/WP10/hypervariability .

SNHG5 promotes colorectal cancer cell survival by counteracting STAU1-mediated mRNA destabilization.

We currently have limited knowledge of the involvement of long non-coding RNAs (lncRNAs) in normal cellular processes and pathologies. Here, we identify and characterize SNHG5 as a stable cytoplasmic lncRNA with up-regulated expression in colorectal cancer. Depletion of SNHG5 induces cell cycle arrest and apoptosis in vitro and limits tumour outgrowth in vivo, whereas SNHG5 overexpression counteracts oxaliplatin-induced apoptosis. Using an unbiased approach, we identify 121 transcript sites interacting with SNHG5 in the cytoplasm. Importantly, knockdown of key SNHG5 target transcripts, including SPATS2, induces apoptosis and thus mimics the effect seen following SNHG5 depletion. Mechanistically, we suggest that SNHG5 stabilizes the target transcripts by blocking their degradation by STAU1. Accordingly, depletion of STAU1 rescues the apoptosis induced after SNHG5 knockdown. Hence, we characterize SNHG5 as a lncRNA promoting tumour cell survival in colorectal cancer and delineate a novel mechanism in which a cytoplasmic lncRNA functions through blocking the action of STAU1.

MicroRNA-941 Expression in Polymorphonuclear Granulocytes Is Not Related to Granulomatosis with Polyangiitis.

Jumonji Domain-Containing Protein 3 (JMJD3)/lysine demethylase 6B (KDM6B) is an epigenetic modulator that removes repressive histone marks on genes. Expression of KDM6B mRNA is elevated in leukocytes from patients with ANCA-associated vasculitis (AAV) and has been suggested to be the reason for higher proteinase 3 (PR3) mRNA expression in these cells due to derepression of PRTN3 gene transcription. MicroRNA-941 (miR-941) has been shown to target KDM6B mRNA and inhibit JMJD3 production. We therefore investigated whether polymorphonuclear granulocytes (PMNs) from patients suffering from granulomatosis with polyangiitis (GPA) have lower expression of miR-941 than healthy control donors as a biological cause for higher JMJD3 levels. We found no significant difference in the degree of maturation of PMNs from GPA patients (n = 8) and healthy controls (n = 11) as determined from cell surface expression of the neutrophil maturation marker CD16 and gene expression profile of FCGR3B. The expression of PRTN3 and KDM6B mRNAs and miR-941 was not significantly different in GPA patients and healthy controls. Transfection of pre-miR-941 into the neutrophil promyelocyte cell line PLB-985 cells did not result in reduction of the KDM6B mRNA level as shown previously in a hepatocellular carcinoma cell line. The amount of PR3 in PMNs from GPA patients and healthy controls was comparable. In conclusion, we found that PRTN3 mRNA, KDM6B mRNA, and miR-941 expression levels in PMNs do not differ between GPA patients and healthy controls, and that miR-941 does not uniformly regulate KDM6B mRNA levels by inducing degradation of the transcript. Thus, decreased miR-941 expression in PMNs cannot be part of the pathogenesis of GPA.

Peak-valley-peak pattern of histone modifications delineates active regulatory elements and their directionality.

Formation of nucleosome free region (NFR) accompanied by specific histone modifications at flanking nucleosomes is an important prerequisite for enhancer and promoter activity. Due to this process, active regulatory elements often exhibit a distinct shape of histone signal in the form of a peak-valley-peak (PVP) pattern. However, different features of PVP patterns and their robustness in predicting active regulatory elements have never been systematically analyzed. Here, we present PARE, a novel computational method that systematically analyzes the H3K4me1 or H3K4me3 PVP patterns to predict NFRs. We show that NFRs predicted by H3K4me1 and me3 patterns are associated with active enhancers and promoters, respectively. Furthermore, asymmetry in the height of peaks flanking the central valley can predict the directionality of stable transcription at promoters. Using PARE on ChIP-seq histone modifications from four ENCODE cell lines and four hematopoietic differentiation stages, we identified several enhancers whose regulatory activity is stage specific and correlates positively with the expression of proximal genes in a particular stage. In conclusion, our results demonstrate that PVP patterns delineate both the histone modification landscape and the transcriptional activities governed by active enhancers and promoters, and therefore can be used for their prediction. PARE is freely available at http://servers.binf.ku.dk/pare.

BloodSpot: a database of gene expression profiles and transcriptional programs for healthy and malignant haematopoiesis.

Research on human and murine haematopoiesis has resulted in a vast number of gene-expression data sets that can potentially answer questions regarding normal and aberrant blood formation. To researchers and clinicians with limited bioinformatics experience, these data have remained available, yet largely inaccessible. Current databases provide information about gene-expression but fail to answer key questions regarding co-regulation, genetic programs or effect on patient survival. To address these shortcomings, we present BloodSpot (www.bloodspot.eu), which includes and greatly extends our previously released database HemaExplorer, a database of gene expression profiles from FACS sorted healthy and malignant haematopoietic cells. A revised interactive interface simultaneously provides a plot of gene expression along with a Kaplan-Meier analysis and a hierarchical tree depicting the relationship between different cell types in the database. The database now includes 23 high-quality curated data sets relevant to normal and malignant blood formation and, in addition, we have assembled and built a unique integrated data set, BloodPool. Bloodpool contains more than 2000 samples assembled from six independent studies on acute myeloid leukemia. Furthermore, we have devised a robust sample integration procedure that allows for sensitive comparison of user-supplied patient samples in a well-defined haematopoietic cellular space.

ERG promotes the maintenance of hematopoietic stem cells by restricting their differentiation.

The balance between self-renewal and differentiation is crucial for the maintenance of hematopoietic stem cells (HSCs). Whereas numerous gene regulatory factors have been shown to control HSC self-renewal or drive their differentiation, we have relatively few insights into transcription factors that serve to restrict HSC differentiation. In the present work, we identify ETS (E-twenty-six)-related gene (ERG) as a critical factor protecting HSCs from differentiation. Specifically, loss of Erg accelerates HSC differentiation by >20-fold, thus leading to rapid depletion of immunophenotypic and functional HSCs. Molecularly, we could demonstrate that ERG, in addition to promoting the expression of HSC self-renewal genes, also represses a group of MYC targets, thereby explaining why Erg loss closely mimics Myc overexpression. Consistently, the BET domain inhibitor CPI-203, known to repress Myc expression, confers a partial phenotypic rescue. In summary, ERG plays a critical role in coordinating the balance between self-renewal and differentiation of HSCs.