BloodChIP: a database of comparative genome-wide transcription factor binding profiles in human blood cells

Peak Scores Significantly Depend on the Relationships between Contextual Signals in ChIP-Seq Peaks

Chromatin immunoprecipitation followed by massively parallel DNA sequencing (ChIP-seq) is a central genome-wide method for in vivo analyses of DNA-protein interactions in various cellular conditions. Numerous studies have demonstrated the complex contextual organization of ChIP-seq peak sequences and the presence of binding sites for transcription factors in them. We assessed the dependence of the ChIP-seq peak score on the presence of different contextual signals in the peak sequences by analyzing these sequences from several ChIP-seq experiments using our fully enumerative GPU-based de novo motif discovery method, Argo_CUDA. Analysis revealed sets of significant IUPAC motifs corresponding to the binding sites of the target and partner transcription factors. For these ChIP-seq experiments, multiple regression models were constructed, demonstrating a significant dependence of the peak scores on the presence in the peak sequences of not only highly significant target motifs but also less significant motifs corresponding to the binding sites of the partner transcription factors. A significant correlation was shown between the presence of the target motifs FOXA2 and the partner motifs HNF4G, which found experimental confirmation in the scientific literature, demonstrating the important contribution of the partner transcription factors to the binding of the target transcription factor to DNA and, consequently, their important contribution to the peak score.

BloodChIP Xtra: an expanded database of comparative genome-wide transcription factor binding and gene-expression profiles in healthy human stem/progenitor subsets and leukemic cells

The BloodChIP Xtra database (http://bloodchipXtra.vafaeelab.com/) facilitates genome-wide exploration and visualization of transcription factor (TF) occupancy and chromatin configuration in rare primary human hematopoietic stem (HSC-MPP) and progenitor (CMP, GMP, MEP) cells and acute myeloid leukemia (AML) cell lines (KG-1, ME-1, Kasumi1, TSU-1621-MT), along with chromatin accessibility and gene expression data from these and primary patient AMLs. BloodChIP Xtra features significantly more datasets than our earlier database BloodChIP (two primary cell types and two cell lines). Improved methodologies for determining TF occupancy and chromatin accessibility have led to increased availability of data for rare primary cell types across the spectrum of healthy and AML hematopoiesis. However, there is a continuing need for these data to be integrated in an easily accessible manner for gene-based queries and use in downstream applications. Here, we provide a user-friendly database based around genome-wide binding profiles of key hematopoietic TFs and histone marks in healthy stem/progenitor cell types. These are compared with binding profiles and chromatin accessibility derived from primary and cell line AML and integrated with expression data from corresponding cell types. All queries can be exported to construct TF-gene and protein-protein networks and evaluate the association of genes with specific cellular processes.

TAL1 hijacks MYCN enhancer that induces MYCN expression and dependence on mevalonate pathway in T-cell acute lymphoblastic leukemia

A hallmark of T-cell acute lymphoblastic leukemia (T-ALL) is the dysregulated expression of oncogenic transcription factors (TFs), including TAL1, NOTCH1 and MYC. Rewiring of the transcriptional program disrupts the tightly controlled spatiotemporal expression of downstream target genes, thereby contributing to leukemogenesis. In this study, we first identify an evolutionarily conserved enhancer element controlling the MYCN oncogene (named enhMYCN) that is aberrantly activated by the TAL1 complex in T-ALL cells. TAL1-positive T-ALL cells are highly dependent on MYCN expression for their maintenance in vitro and in xenograft models. Interestingly, MYCN drives the expression of multiple genes involved in the mevalonate pathway, and T-ALL cells are sensitive to inhibition of HMG-CoA reductase (HMGCR), a rate-limiting enzyme of this pathway. Importantly, MYC and MYCN regulate the same targets and compensate for each other. Thus, MYCN-positive T-ALL cells display a dual dependence on the TAL1-MYCN and NOTCH1-MYC pathways. Together, our results demonstrate that enhMYCN-mediated MYCN expression is required for human T-ALL cells and implicate the TAL1-MYCN-HMGCR axis as a potential therapeutic target in T-ALL.

EVI1-mediated Programming of Normal and Malignant Hematopoiesis

Ecotropic viral integration site 1 (EVI1), encoded at the MECOM locus, is an oncogenic zinc finger transcription factor with diverse roles in normal and malignant cells, most extensively studied in the context of hematopoiesis. EVI1 interacts with other transcription factors in a context-dependent manner and regulates transcription and chromatin remodeling, thereby influencing the proliferation, differentiation, and survival of cells. Interestingly, it can act both as a transcriptional activator as well as a transcriptional repressor. EVI1 is expressed, and fulfills important functions, during the development of different tissues, including the nervous system and hematopoiesis, demonstrating a rigid spatial and temporal expression pattern. However, EVI1 is regularly overexpressed in a variety of cancer entities, including epithelial cancers such as ovarian and pancreatic cancer, as well as in hematologic malignancies like myeloid leukemias. Importantly, EVI1 overexpression is generally associated with a very poor clinical outcome and therapy-resistance. Thus, EVI1 is an interesting candidate to study to improve the prognosis and treatment of high-risk patients with "EVI1high" hematopoietic malignancies.

Transcription factor genetics and biology in predisposition to bone marrow failure and hematological malignancy

Transcription factors (TFs) play a critical role as key mediators of a multitude of developmental pathways, with highly regulated and tightly organized networks crucial for determining both the timing and pattern of tissue development. TFs can act as master regulators of both primitive and definitive hematopoiesis, tightly controlling the behavior of hematopoietic stem and progenitor cells (HSPCs). These networks control the functional regulation of HSPCs including self-renewal, proliferation, and differentiation dynamics, which are essential to normal hematopoiesis. Defining the key players and dynamics of these hematopoietic transcriptional networks is essential to understanding both normal hematopoiesis and how genetic aberrations in TFs and their networks can predispose to hematopoietic disease including bone marrow failure (BMF) and hematological malignancy (HM). Despite their multifaceted and complex involvement in hematological development, advances in genetic screening along with elegant multi-omics and model system studies are shedding light on how hematopoietic TFs interact and network to achieve normal cell fates and their role in disease etiology. This review focuses on TFs which predispose to BMF and HM, identifies potential novel candidate predisposing TF genes, and examines putative biological mechanisms leading to these phenotypes. A better understanding of the genetics and molecular biology of hematopoietic TFs, as well as identifying novel genes and genetic variants predisposing to BMF and HM, will accelerate the development of preventative strategies, improve clinical management and counseling, and help define targeted treatments for these diseases.

Germline RUNX1 variation and predisposition to childhood acute lymphoblastic leukemia

Genetic alterations in the RUNX1 gene are associated with benign and malignant blood disorders, particularly of megakaryocyte and myeloid lineages. The role of RUNX1 in acute lymphoblastic leukemia (ALL) is less clear, particularly how germline genetic variation influences the predisposition to this type of leukemia. Sequencing 4,836 children with B-ALL and 1,354 cases of T-ALL, we identified 31 and 18 germline RUNX1 variants, respectively. RUNX1 variants in B-ALL consistently showed minimal damaging effects. By contrast, 6 T-ALL-related variants result in drastic loss of RUNX1 activity as a transcription activator in vitro. Ectopic expression of dominant-negative RUNX1 variants in human CD34+ cells repressed differentiation into erythroid, megakaryocytes, and T cells, while promoting myeloid cell development. Chromatin immunoprecipitation sequencing of T-ALL models showed distinctive patterns of RUNX1 binding by variant proteins. Further whole genome sequencing identified JAK3 mutation as the most frequent somatic genomic abnormality in T-ALL with germline RUNX1 variants. Co-introduction of RUNX1 variant and JAK3 mutation in hematopoietic stem and progenitor cells in mice gave rise to T-ALL with early T-cell precursor phenotype. Taken together, these results indicated that RUNX1 is an important predisposition gene for T-ALL and pointed to novel biology of RUNX1-mediated leukemogenesis in the lymphoid lineages.

Atypical 3q26/MECOM rearrangements genocopy inv(3)/t(3;3) in acute myeloid leukemia

Acute myeloid leukemia (AML) with inv(3)/t(3;3)(q21q26) is a distinct World Health Organization recognized entity, characterized by its aggressive course and poor prognosis. In this subtype of AML, the translocation of a GATA2 enhancer (3q21) to MECOM (3q26) results in overexpression of the MECOM isoform EVI1 and monoallelic expression of GATA2 from the unaffected allele. The full-length MECOM transcript, MDS1-EVI1, is not expressed as the result of the 3q26 rearrangement. Besides the classical inv(3)/t(3;3), a number of other 3q26/MECOM rearrangements with poor treatment response have been reported in AML. Here, we demonstrate, in a group of 33 AML patients with atypical 3q26 rearrangements, MECOM involvement with EVI1 overexpression but no or low MDS1-EVI1 levels. Moreover, the 3q26 translocations in these AML patients often involve superenhancers of genes active in myeloid development (eg, CD164, PROM1, CDK6, or MYC). In >50% of these cases, allele-specific GATA2 expression was observed, either by copy-number loss or by an unexplained allelic imbalance. Altogether, atypical 3q26 recapitulate the main leukemic mechanism of inv(3)/t(3;3) AML, namely EVI1 overexpression driven by enhancer hijacking, absent MDS1-EVI1 expression and potential GATA2 involvement. Therefore, we conclude that both atypical 3q26/MECOM and inv(3)/t(3;3) can be classified as a single entity of 3q26-rearranged AMLs. Routine analyses determining MECOM rearrangements and EVI1 and MDS1-EVI1 expression are required to recognize 3q-rearranged AML cases.

ChIPSummitDB: a ChIP-seq-based database of human transcription factor binding sites and the topological arrangements of the proteins bound to them

ChIP-seq reveals genomic regions where proteins, e.g. transcription factors (TFs) interact with DNA. A substantial fraction of these regions, however, do not contain the cognate binding site for the TF of interest. This phenomenon might be explained by protein–protein interactions and co-precipitation of interacting gene regulatory elements. We uniformly processed 3727 human ChIP-seq data sets and determined the cistrome of 292 TFs, as well as the distances between the TF binding motif centers and the ChIP-seq peak summits. ChIPSummitDB enables the analysis of ChIP-seq data using multiple approaches. The 292 cistromes and corresponding ChIP-seq peak sets can be browsed in GenomeView. Overlapping SNPs can be inspected in dbSNPView. Most importantly, the MotifView and PairShiftView pages show the average distance between motif centers and overlapping ChIP-seq peak summits and distance distributions thereof, respectively. In addition to providing a comprehensive human TF binding site collection, the ChIPSummitDB database and web interface allows for the examination of the topological arrangement of TF complexes genome-wide. ChIPSummitDB is freely accessible at http://summit.med.unideb.hu/summitdb/. The database will be regularly updated and extended with the newly available human and mouse ChIP-seq data sets.

Somatic Mutations Drive Specific, but Reversible, Epigenetic Heterogeneity States in AML

Epigenetic allele diversity is linked to inferior prognosis in acute myeloid leukemia (AML). However, the source of epiallele heterogeneity in AML is unknown. Herein we analyzed epiallele diversity in a genetically and clinically annotated AML cohort. Notably, AML driver mutations linked to transcription factors and favorable outcome are associated with epigenetic destabilization in a defined set of susceptible loci. In contrast, AML subtypes linked to inferior prognosis manifest greater abundance and highly stochastic epiallele patterning. We report an epiallele outcome classifier supporting the link between epigenetic diversity and treatment failure. Mouse models with TET2 or IDH2 mutations show that epiallele diversity is especially strongly induced by IDH mutations, precedes transformation to AML, and is enhanced by cooperation between somatic mutations. Furthermore, epiallele complexity was partially reversed by epigenetic therapies in AML driven by TET2/IDH2, suggesting that epigenetic therapy might function in part by reducing population complexity and fitness of AMLs. SIGNIFICANCE: We show for the first time that epigenetic clonality is directly linked to specific mutations and that epigenetic allele diversity precedes and potentially contributes to malignant transformation. Furthermore, epigenetic clonality is reversible with epigenetic therapy agents.This article is highlighted in the In This Issue feature, p. 1775.

Endogenous retroviruses are a source of enhancers with oncogenic potential in acute myeloid leukaemia

Acute myeloid leukemia (AML) is characterised by a series of genetic and epigenetic alterations that result in deregulation of transcriptional networks. One understudied source of transcriptional regulators are transposable elements (TEs), whose aberrant usage could contribute to oncogenic transcriptional circuits. However, the regulatory influence of TEs and their links to AML pathogenesis remain unexplored. Here we identify six endogenous retrovirus (ERV) families with AML-associated enhancer chromatin signatures that are enriched in binding of key regulators of hematopoiesis and AML pathogenesis. Using both locus-specific genetic editing and simultaneous epigenetic silencing of multiple ERVs, we demonstrate that ERV deregulation directly alters the expression of adjacent genes in AML. Strikingly, deletion or epigenetic silencing of an ERV-derived enhancer suppresses cell growth by inducing apoptosis in leukemia cell lines. This work reveals that ERVs are a previously unappreciated source of AML enhancers that may be exploited by cancer cells to help drive tumour heterogeneity and evolution.

Transposable elements are a potential source of transcriptional regulators, but how these sequences contribute to oncogenesis remains poorly understood. Here, the authors identify endogenous retroviruses (ERVs) with acute myeloid leukemia (AML)-associated enhancer chromatin signatures, and provide evidence that ERV activation provides an additional layer of gene regulation in AML.

Bayesian modelling of high-throughput sequencing assays with malacoda

NGS studies have uncovered an ever-growing catalog of human variation while leaving an enormous gap between observed variation and experimental characterization of variant function. High-throughput screens powered by NGS have greatly increased the rate of variant functionalization, but the development of comprehensive statistical methods to analyze screen data has lagged. In the massively parallel reporter assay (MPRA), short barcodes are counted by sequencing DNA libraries transfected into cells and the cell's output RNA in order to simultaneously measure the shifts in transcription induced by thousands of genetic variants. These counts present many statistical challenges, including overdispersion, depth dependence, and uncertain DNA concentrations. So far, the statistical methods used have been rudimentary, employing transformations on count level data and disregarding experimental and technical structure while failing to quantify uncertainty in the statistical model. We have developed an extensive framework for the analysis of NGS functionalization screens available as an R package called malacoda (available from github.com/andrewGhazi/malacoda). Our software implements a probabilistic, fully Bayesian model of screen data. The model uses the negative binomial distribution with gamma priors to model sequencing counts while accounting for effects from input library preparation and sequencing depth. The method leverages the high-throughput nature of the assay to estimate the priors empirically. External annotations such as ENCODE data or DeepSea predictions can also be incorporated to obtain more informative priors-a transformative capability for data integration. The package also includes quality control and utility functions, including automated barcode counting and visualization methods. To validate our method, we analyzed several datasets using malacoda and alternative MPRA analysis methods. These data include experiments from the literature, simulated assays, and primary MPRA data. We also used luciferase assays to experimentally validate several hits from our primary data, as well as variants for which the various methods disagree and variants detectable only with the aid of external annotations.

Loss of lenalidomide-induced megakaryocytic differentiation leads to therapy resistance in del(5q) myelodysplastic syndrome

Interstitial deletion of the long arm of chromosome 5 (del(5q)) is the most common structural genomic variant in myelodysplastic syndromes (MDS)¹. Lenalidomide (LEN) is the treatment of choice for patients with del(5q) MDS, but half of the responding patients become resistant² within 2 years. TP53 mutations are detected in ~20% of LEN-resistant patients³. Here we show that patients who become resistant to LEN harbour recurrent variants of TP53 or RUNX1. LEN upregulated RUNX1 protein and function in a CRBN- and TP53-dependent manner in del(5q) cells, and mutation or downregulation of RUNX1 rendered cells resistant to LEN. LEN induced megakaryocytic differentiation of del(5q) cells followed by cell death that was dependent on calpain activation and CSNK1A1 degradation^4,5. We also identified GATA2 as a LEN-responsive gene that is required for LEN-induced megakaryocyte differentiation. Megakaryocytic gene-promoter analyses suggested that LEN-induced degradation of IKZF1 enables a RUNX1-GATA2 complex to drive megakaryocytic differentiation. Overexpression of GATA2 restored LEN sensitivity in the context of RUNX1 or TP53 mutations by enhancing LEN-induced megakaryocytic differentiation. Screening for mutations that block LEN-induced megakaryocytic differentiation should identify patients who are resistant to LEN.

A multi-centre study on the performance of the molecular genotyping platform ID RHD XT for resolving serological weak RhD phenotype in routine clinical practice

BACKGROUND AND OBJECTIVES

There is concern regarding the lack of prevention of unnecessary transfusion of RhD negative red cells and unnecessary administration of Rh immunoglobulin (RhIG) to pregnant women. In this study, performance of ID RHD XT, a genotyping assay for identification of six RHD allelic variants and human platelet antigens HPA-1a/1b was assessed.

MATERIALS AND METHODS

Whole blood samples presenting weak, discrepant or inconclusive D phenotype results were genotyped with ID RHD XT and compared to reference molecular tests. Candidacy for RhIG prophylaxis was determined by analysing samples from pregnant women. Hands-on time to complete the procedures was measured.

RESULTS

Overall, 167 samples were tested (55 donors, 56 patients, 52 pregnant women and four newborns). Agreement between ID RHD XT and the reference method was 100% (51% weak D type 1, 2 or 3; 35·5% weak D Types 1, 2 or 3 not detected; 4% RHD deletion; 1% RHD*Pseudogene; 1% RHD*DIIIa-CE(3-7)-D; and 4% no amplification variant detected for RHD genotype; and 64% HPA-1a/a; 30% HPA-1a/b; and 3% HPA-1b/b for HPA-1 genotype). Call rate was 98·2%. ID RHD XT identified 40% of the pregnant women that would not have required RhIG prophylaxis. Overall hands-on time was 25-45 min to process a batch of 24 samples, and four hours for total assay time.

CONCLUSION

ID RHD XT yielded reproducible results for RHD typing in serologically weak D phenotype individuals. ID RHD XT was proven useful for the correct management of patients with RhD serological discrepancies and the rational use of RhIG in pregnancy.

Transcriptional networks in acute myeloid leukemia

Acute myeloid leukemia (AML) is a complex disease characterized by a diverse range of recurrent molecular aberrations that occur in many different combinations. Components of transcriptional networks are a common target of these aberrations, leading to network-wide changes and deployment of novel or developmentally inappropriate transcriptional programs. Genome-wide techniques are beginning to reveal the full complexity of normal hematopoietic stem cell transcriptional networks and the extent to which they are deregulated in AML, and new understandings of the mechanisms by which AML cells maintain self-renewal and block differentiation are starting to emerge. The hope is that increased understanding of the network architecture in AML will lead to identification of key oncogenic dependencies that are downstream of multiple network aberrations, and that this knowledge will be translated into new therapies that target these dependencies. Here, we review the current state of knowledge of network perturbation in AML with a focus on major mechanisms of transcription factor dysregulation, including mutation, translocation, and transcriptional dysregulation, and discuss how these perturbations propagate across transcriptional networks. We will also review emerging mechanisms of network disruption, and briefly discuss how increased knowledge of network disruption is already being used to develop new therapies.

A Forward Genetic Screen Targeting the Endothelium Reveals a Regulatory Role for the Lipid Kinase Pi4ka in Myelo- and Erythropoiesis

Given its role as the source of definitive hematopoietic cells, we sought to determine whether mutations initiated in the hemogenic endothelium would yield hematopoietic abnormalities or malignancies. Here, we find that endothelium-specific transposon mutagenesis in mice promotes hematopoietic pathologies that are both myeloid and lymphoid in nature. Frequently mutated genes included previously recognized cancer drivers and additional candidates, such as Pi4ka, a lipid kinase whose mutation was found to promote myeloid and erythroid dysfunction. Subsequent validation experiments showed that targeted inactivation of the Pi4ka catalytic domain or reduction in mRNA expression inhibited myeloid and erythroid cell differentiation in vitro and promoted anemia in vivo through a mechanism involving deregulation of AKT, MAPK, SRC, and JAK-STAT signaling. Finally, we provide evidence linking PI4KAP2, previously considered a pseudogene, to human myeloid and erythroid leukemia.

FUSE binding protein 1 (FUBP1) expression is upregulated by T-cell acute lymphocytic leukemia protein 1 (TAL1) and required for efficient erythroid differentiation

During erythropoiesis, haematopoietic stem cells (HSCs) differentiate in successive steps of commitment and specification to mature erythrocytes. This differentiation process is controlled by transcription factors that establish stage- and cell type-specific gene expression. In this study, we demonstrate that FUSE binding protein 1 (FUBP1), a transcriptional regulator important for HSC self-renewal and survival, is regulated by T-cell acute lymphocytic leukaemia 1 (TAL1) in erythroid progenitor cells. TAL1 directly activates the FUBP1 promoter, leading to increased FUBP1 expression during erythroid differentiation. The binding of TAL1 to the FUBP1 promoter is highly dependent on an intact GATA sequence in a combined E-box/GATA motif. We found that FUBP1 expression is required for efficient erythropoiesis, as FUBP1-deficient progenitor cells were limited in their potential of erythroid differentiation. Thus, the finding of an interconnection between GATA1/TAL1 and FUBP1 reveals a molecular mechanism that is part of the switch from progenitor- to erythrocyte-specific gene expression. In summary, we identified a TAL1/FUBP1 transcriptional relationship, whose physiological function in haematopoiesis is connected to proper erythropoiesis.

Lymphocyte-Specific Chromatin Accessibility Pre-determines Glucocorticoid Resistance in Acute Lymphoblastic Leukemia

Glucocorticoids play a critical role in the treatment of lymphoid malignancies. While glucocorticoid efficacy can be largely attributed to lymphocyte-specific apoptosis, its molecular basis remains elusive. Here, we studied genome-wide lymphocyte-specific open chromatin domains (LSOs), and integrated LSOs with glucocorticoid-induced RNA transcription and chromatin modulation using an in vivo patient-derived xenograft model of acute lymphoblastic leukemia (ALL). This led to the identification of LSOs critical for glucocorticoid-induced apoptosis. Glucocorticoid receptor cooperated with CTCF at these LSOs to mediate DNA looping, which was inhibited by increased DNA methylation in glucocorticoid-resistant ALL and non-lymphoid cell types. Our study demonstrates that lymphocyte-specific epigenetic modifications pre-determine glucocorticoid resistance in ALL and may account for the lack of glucocorticoid sensitivity in other cell types.

A novel role for Lyl1 in primitive erythropoiesis

Stem cell leukemia (Scl or Tal1) and lymphoblastic leukemia 1 (Lyl1) encode highly related members of the basic helix-loop-helix family of transcription factors that are co-expressed in the erythroid lineage. Previous studies have suggested that Scl is essential for primitive erythropoiesis. However, analysis of single-cell RNA-seq data of early embryos showed that primitive erythroid cells express both Scl and Lyl1 Therefore, to determine whether Lyl1 can function in primitive erythropoiesis, we crossed conditional Scl knockout mice with mice expressing a Cre recombinase under the control of the Epo receptor, active in erythroid progenitors. Embryos with 20% expression of Scl from E9.5 survived to adulthood. However, mice with reduced expression of Scl and absence of Lyl1 (double knockout; DKO) died at E10.5 because of progressive loss of erythropoiesis. Gene expression profiling of DKO yolk sacs revealed loss of Gata1 and many of the known target genes of the SCL-GATA1 complex. ChIP-seq analyses in a human erythroleukemia cell line showed that LYL1 exclusively bound a small subset of SCL targets including GATA1. Together, these data show for the first time that Lyl1 can maintain primitive erythropoiesis.

RUNX1c Regulates Hematopoietic Differentiation of Human Pluripotent Stem Cells Possibly in Cooperation with Proinflammatory Signaling

Runt-related transcription factor 1 (Runx1) is a master hematopoietic transcription factor essential for hematopoietic stem cell (HSC) emergence. Runx1-deficient mice die during early embryogenesis due to the inability to establish definitive hematopoiesis. Here, we have used human pluripotent stem cells (hPSCs) as model to study the role of RUNX1 in human embryonic hematopoiesis. Although the three RUNX1 isoforms a, b, and c were induced in CD45+ hematopoietic cells, RUNX1c was the only isoform induced in hematoendothelial progenitors (HEPs)/hemogenic endothelium. Constitutive expression of RUNX1c in human embryonic stem cells enhanced the appearance of HEPs, including hemogenic (CD43+) HEPs and promoted subsequent differentiation into blood cells. Conversely, specific deletion of RUNX1c dramatically reduced the generation of hematopoietic cells from HEPs, indicating that RUNX1c is a master regulator of human hematopoietic development. Gene expression profiling of HEPs revealed a RUNX1c-induced proinflammatory molecular signature, supporting previous studies demonstrating proinflammatory signaling as a regulator of HSC emergence. Collectively, RUNX1c orchestrates hematopoietic specification of hPSCs, possibly in cooperation with proinflammatory signaling. Stem Cells 2017;35:2253-2266.

Identification of a functional genetic variant driving racially dimorphic platelet gene expression of the thrombin receptor regulator, PCTP

Platelet activation in response to stimulation of the Protease Activated Receptor 4 (PAR4) receptor differs by race. One factor that contributes to this difference is the expression level of Phosphatidylcholine Transfer Protein (PCTP), a regulator of platelet PAR4 function. We have conducted an expression Quantitative Trait Locus (eQTL) analysis that identifies single nucleotide polymorphisms (SNPs) linked to the expression level of platelet genes. This analysis revealed 26 SNPs associated with the expression level of PCTP at genome-wide significance (p < 5×10^-8). Using annotation from ENCODE and other public data we prioritised one of these SNPs, rs2912553, for functional testing. The allelic frequency of rs2912553 is racially-dimorphic, in concordance with the racially differential expression of PCTP. Reporter gene assays confirmed that the single nucleotide change caused by rs2912553 altered the transcriptional potency of the surrounding genomic locus. Electromobility shift assays, luciferase assays, and overexpression studies indicated a role for the megakaryocytic transcription factor GATA1. In summary, we have integrated multi-omic data to identify and functionalise an eQTL. This, along with the previously described relationship between PCTP and PAR4 function, allows us to characterise a genotype-phenotype relationship through the mechanism of gene expression.

A quantitative proteomics approach identifies ETV6 and IKZF1 as new regulators of an ERG-driven transcriptional network

Aberrant stem cell-like gene regulatory networks are a feature of leukaemogenesis. The ETS-related gene (ERG), an important regulator of normal haematopoiesis, is also highly expressed in T-ALL and acute myeloid leukaemia (AML). However, the transcriptional regulation of ERG in leukaemic cells remains poorly understood. In order to discover transcriptional regulators of ERG, we employed a quantitative mass spectrometry-based method to identify factors binding the 321 bp ERG +85 stem cell enhancer region in MOLT-4 T-ALL and KG-1 AML cells. Using this approach, we identified a number of known binders of the +85 enhancer in leukaemic cells along with previously unknown binders, including ETV6 and IKZF1. We confirmed that ETV6 and IKZF1 were also bound at the +85 enhancer in both leukaemic cells and in healthy human CD34+ haematopoietic stem and progenitor cells. Knockdown experiments confirmed that ETV6 and IKZF1 are transcriptional regulators not just of ERG, but also of a number of genes regulated by a densely interconnected network of seven transcription factors. At last, we show that ETV6 and IKZF1 expression levels are positively correlated with expression of a number of heptad genes in AML and high expression of all nine genes confers poorer overall prognosis.

Biological Databases for Hematology Research

With the advances of genome-wide sequencing technologies and bioinformatics approaches, a large number of datasets of normal and malignant erythropoiesis have been generated and made public to researchers around the world. Collection and integration of these datasets greatly facilitate basic research and clinical diagnosis and treatment of blood disorders. Here we provide a brief introduction of the most popular omics data resources of normal and malignant hematopoiesis, including some integrated web tools, to help users get better equipped to perform common analyses. We hope this review will promote the awareness and facilitate the usage of public database resources in the hematology research.

Cistrome Data Browser: a data portal for ChIP-Seq and chromatin accessibility data in human and mouse

Chromatin immunoprecipitation, DNase I hypersensitivity and transposase-accessibility assays combined with high-throughput sequencing enable the genome-wide study of chromatin dynamics, transcription factor binding and gene regulation. Although rapidly accumulating publicly available ChIP-seq, DNase-seq and ATAC-seq data are a valuable resource for the systematic investigation of gene regulation processes, a lack of standardized curation, quality control and analysis procedures have hindered extensive reuse of these data. To overcome this challenge, we built the Cistrome database, a collection of ChIP-seq and chromatin accessibility data (DNase-seq and ATAC-seq) published before January 1, 2016, including 13 366 human and 9953 mouse samples. All the data have been carefully curated and processed with a streamlined analysis pipeline and evaluated with comprehensive quality control metrics. We have also created a user-friendly web server for data query, exploration and visualization. The resulting Cistrome DB (Cistrome Data Browser), available online at http://cistrome.org/db, is expected to become a valuable resource for transcriptional and epigenetic regulation studies.

Hematopoietic (stem) cell development - how divergent are the roads taken?

The development of the hematopoietic system during early embryonic stages occurs in spatially and temporally distinct waves. Hematopoietic stem cells (HSC), the most potent and self‐renewing cells of this system, are produced in the final ‘definitive’ wave of hematopoietic cell generation. In contrast to HSCs in the adult, which differentiate via intermediate progenitor populations to produce functional blood cells, the generation of hematopoietic cells in the embryo prior to HSC generation occurs in the early waves by producing blood cells without intermediate progenitors (such as the ‘primitive’ hematopoietic cells). The lineage relationship between the early hematopoietic cells and the cells giving rise to HSCs, the genetic networks controlling their emergence, and the precise temporal determination of HSC fate remain topics of intense research and debate. This Review article discusses the current knowledge on the step‐wise embryonic establishment of the adult hematopoietic system, examines the roles of pivotal intrinsic regulators in this process, and raises questions concerning the temporal onset of HSC fate determination.

CardioTF, a database of deconstructing transcriptional circuits in the heart system

Background: Information on cardiovascular gene transcription is fragmented and far behind the present requirements of the systems biology field. To create a comprehensive source of data for cardiovascular gene regulation and to facilitate a deeper understanding of genomic data, the CardioTF database was constructed. The purpose of this database is to collate information on cardiovascular transcription factors (TFs), position weight matrices (PWMs), and enhancer sequences discovered using the ChIP-seq method.

Methods: The Naïve-Bayes algorithm was used to classify literature and identify all PubMed abstracts on cardiovascular development. The natural language learning tool GNAT was then used to identify corresponding gene names embedded within these abstracts. Local Perl scripts were used to integrate and dump data from public databases into the MariaDB management system (MySQL). In-house R scripts were written to analyze and visualize the results.

Results: Known cardiovascular TFs from humans and human homologs from fly, Ciona, zebrafish, frog, chicken, and mouse were identified and deposited in the database. PWMs from Jaspar, hPDI, and UniPROBE databases were deposited in the database and can be retrieved using their corresponding TF names. Gene enhancer regions from various sources of ChIP-seq data were deposited into the database and were able to be visualized by graphical output. Besides biocuration, mouse homologs of the 81 core cardiac TFs were selected using a Naïve-Bayes approach and then by intersecting four independent data sources: RNA profiling, expert annotation, PubMed abstracts and phenotype.

Discussion: The CardioTF database can be used as a portal to construct transcriptional network of cardiac development.

Availability and Implementation: Database URL: http://www.cardiosignal.org/database/cardiotf.html.

SBR-Blood: systems biology repository for hematopoietic cells

Extensive research into hematopoiesis (the development of blood cells) over several decades has generated large sets of expression and epigenetic profiles in multiple human and mouse blood cell types. However, there is no single location to analyze how gene regulatory processes lead to different mature blood cells. We have developed a new database framework called hematopoietic Systems Biology Repository (SBR-Blood), available online at http://sbrblood.nhgri.nih.gov, which allows user-initiated analyses for cell type correlations or gene-specific behavior during differentiation using publicly available datasets for array- and sequencing-based platforms from mouse hematopoietic cells. SBR-Blood organizes information by both cell identity and by hematopoietic lineage. The validity and usability of SBR-Blood has been established through the reproduction of workflows relevant to expression data, DNA methylation, histone modifications and transcription factor occupancy profiles.

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.

Role of ChIP-seq in the discovery of transcription factor binding sites, differential gene regulation mechanism, epigenetic marks and beyond

Many biologically significant processes, such as cell differentiation and cell cycle progression, gene transcription and DNA replication, chromosome stability and epigenetic silencing etc. depend on the crucial interactions between cellular proteins and DNA. Chromatin immunoprecipitation (ChIP) is an important experimental technique for studying interactions between specific proteins and DNA in the cell and determining their localization on a specific genomic locus. In recent years, the combination of ChIP with second generation DNA-sequencing technology (ChIP-seq) allows precise genomic functional assay. This review addresses the important applications of ChIP-seq with an emphasis on its role in genome-wide mapping of transcription factor binding sites, the revelation of underlying molecular mechanisms of differential gene regulation that are governed by specific transcription factors, and the identification of epigenetic marks. Furthermore, we also describe the ChIP-seq data analysis workflow and a perspective for the exciting potential advancement of ChIP-seq technology in the future.

ING3 protein expression profiling in normal human tissues suggest its role in cellular growth and self-renewal

Members of the INhibitor of Growth (ING) family of proteins act as readers of the epigenetic code through specific recognition of the trimethylated form of lysine 4 of histone H3 (H3K4Me3) by their plant homeodomains. The founding member of the family, ING1, was initially identified as a tumor suppressor with altered regulation in a variety of cancer types. While alterations in ING1 and ING4 levels have been reported in a variety of cancer types, little is known regarding ING3 protein levels in normal or transformed cells due to a lack of reliable immunological tools. In this study we present the characterization of a new monoclonal antibody we have developed against ING3 that specifically recognizes human and mouse ING3. The antibody works in western blots, immunofluorescence, immunoprecipitation and immunohistochemistry. Using this antibody we show that ING3 is most highly expressed in small intestine, bone marrow and epidermis, tissues in which cells undergo rapid proliferation and renewal. Consistent with this observation, we show that ING3 is expressed at significantly higher levels in proliferating versus quiescent epithelial cells. These data suggest that ING3 levels may serve as a surrogate for growth rate, and suggest possible roles for ING3 in growth and self renewal and related diseases such as cancer.

Whole-transcriptome analysis of endothelial to hematopoietic stem cell transition reveals a requirement for Gpr56 in HSC generation

Using highly sensitive RNAseq to examine the whole transcriptome of enriched aortic hematopoietic stem cells and endothelial cells, the authors find G-protein–coupled receptor, Gpr56, is required to generate the first HSCs during endothelial to hematopoietic cell transition.

Hematopoietic stem cells (HSCs) are generated via a natural transdifferentiation process known as endothelial to hematopoietic cell transition (EHT). Because of small numbers of embryonal arterial cells undergoing EHT and the paucity of markers to enrich for hemogenic endothelial cells (ECs [HECs]), the genetic program driving HSC emergence is largely unknown. Here, we use a highly sensitive RNAseq method to examine the whole transcriptome of small numbers of enriched aortic HSCs, HECs, and ECs. Gpr56, a G-coupled protein receptor, is one of the most highly up-regulated of the 530 differentially expressed genes. Also, highly up-regulated are hematopoietic transcription factors, including the “heptad” complex of factors. We show that Gpr56 (mouse and human) is a target of the heptad complex and is required for hematopoietic cluster formation during EHT. Our results identify the processes and regulators involved in EHT and reveal the surprising requirement for Gpr56 in generating the first HSCs.

Opposing regulation of BIM and BCL2 controls glucocorticoid-induced apoptosis of pediatric acute lymphoblastic leukemia cells

Glucocorticoids are critical components of combination chemotherapy regimens in pediatric acute lymphoblastic leukemia (ALL). The proapoptotic BIM protein is an important mediator of glucocorticoid-induced apoptosis in normal and malignant lymphocytes, whereas the antiapoptotic BCL2 confers resistance. The signaling pathways regulating BIM and BCL2 expression in glucocorticoid-treated lymphoid cells remain unclear. In this study, pediatric ALL patient-derived xenografts (PDXs) inherently sensitive or resistant to glucocorticoids were exposed to dexamethasone in vivo. Microarray analysis showed that KLF13 and MYB gene expression changes were significantly greater in dexamethasone-sensitive than -resistant PDXs. Chromatin immunoprecipitation (ChIP) analysis detected glucocorticoid receptor (GR) binding at the KLF13 promoter to trigger KLF13 expression only in sensitive PDXs. Next, KLF13 bound to the MYB promoter, deactivating MYB expression only in sensitive PDXs. Sustained MYB expression in resistant PDXs resulted in maintenance of BCL2 expression and inhibition of apoptosis. ChIP sequencing analysis revealed a novel GR binding site in a BIM intronic region (IGR) that was engaged only in dexamethasone-sensitive PDXs. The absence of GR binding at the BIM IGR was associated with BIM silencing and dexamethasone resistance. This study has identified novel mechanisms of opposing BCL2 and BIM gene regulation that control glucocorticoid-induced apoptosis in pediatric ALL cells in vivo.

Overexpression of ERG in cord blood progenitors promotes expansion and recapitulates molecular signatures of high ERG leukemias

High expression of the ETS family transcription factor ERG is associated with poor clinical outcome in acute myeloid leukemia (AML) and acute T-cell lymphoblastic leukemia (T-ALL). In murine models, high ERG expression induces both T-ALL and AML. However, no study to date has defined the effect of high ERG expression on primary human hematopoietic cells. In the present study, human CD34+ cells were transduced with retroviral vectors to elevate ERG gene expression to levels detected in high ERG AML. RNA sequencing was performed on purified populations of transduced cells to define the effects of high ERG on gene expression in human CD34+ cells. Integration of the genome-wide expression data with other data sets revealed that high ERG drives an expression signature that shares features of normal hematopoietic stem cells, high ERG AMLs, early T-cell precursor-ALLs and leukemic stem cell signatures associated with poor clinical outcome. Functional assays linked this gene expression profile to enhanced progenitor cell expansion. These results support a model whereby a stem cell gene expression network driven by high ERG in human cells enhances the expansion of the progenitor pool, providing opportunity for the acquisition and propagation of mutations and the development of leukemia.

CODEX: a next-generation sequencing experiment database for the haematopoietic and embryonic stem cell communities

CODEX (http://codex.stemcells.cam.ac.uk/) is a user-friendly database for the direct access and interrogation of publicly available next-generation sequencing (NGS) data, specifically aimed at experimental biologists. In an era of multi-centre genomic dataset generation, CODEX provides a single database where these samples are collected, uniformly processed and vetted. The main drive of CODEX is to provide the wider scientific community with instant access to high-quality NGS data, which, irrespective of the publishing laboratory, is directly comparable. CODEX allows users to immediately visualize or download processed datasets, or compare user-generated data against the database's cumulative knowledge-base. CODEX contains four types of NGS experiments: transcription factor chromatin immunoprecipitation coupled to high-throughput sequencing (ChIP-Seq), histone modification ChIP-Seq, DNase-Seq and RNA-Seq. These are largely encompassed within two specialized repositories, HAEMCODE and ESCODE, which are focused on haematopoiesis and embryonic stem cell samples, respectively. To date, CODEX contains over 1000 samples, including 221 unique TFs and 93 unique cell types. CODEX therefore provides one of the most complete resources of publicly available NGS data for the direct interrogation of transcriptional programmes that regulate cellular identity and fate in the context of mammalian development, homeostasis and disease.

Transcriptional network control of normal and leukaemic haematopoiesis

Transcription factors (TFs) play a key role in determining the gene expression profiles of stem/progenitor cells, and defining their potential to differentiate into mature cell lineages. TF interactions within gene-regulatory networks are vital to these processes, and dysregulation of these networks by TF overexpression, deletion or abnormal gene fusions have been shown to cause malignancy. While investigation of these processes remains a challenge, advances in genome-wide technologies and growing interactions between laboratory and computational science are starting to produce increasingly accurate network models. The haematopoietic system provides an attractive experimental system to elucidate gene regulatory mechanisms, and allows experimental investigation of both normal and dysregulated networks. In this review we examine the principles of TF-controlled gene regulatory networks and the key experimental techniques used to investigate them. We look in detail at examples of how these approaches can be used to dissect out the regulatory mechanisms controlling normal haematopoiesis, as well as the dysregulated networks associated with haematological malignancies.

Hematopoietic transcriptional mechanisms: from locus-specific to genome-wide vantage points

Hematopoiesis is an exquisitely regulated process in which stem cells in the developing embryo and the adult generate progenitor cells that give rise to all blood lineages. Master regulatory transcription factors control hematopoiesis by integrating signals from the microenvironment and dynamically establishing and maintaining genetic networks. One of the most rudimentary aspects of cell type-specific transcription factor function, how they occupy a highly restricted cohort of cis-elements in chromatin, remains poorly understood. Transformative technologic advances involving the coupling of next-generation DNA sequencing technology with the chromatin immunoprecipitation assay (ChIP-seq) have enabled genome-wide mapping of factor occupancy patterns. However, formidable problems remain; notably, ChIP-seq analysis yields hundreds to thousands of chromatin sites occupied by a given transcription factor, and only a fraction of the sites appear to be endowed with critical, non-redundant function. It has become en vogue to map transcription factor occupancy patterns genome-wide, while using powerful statistical tools to establish correlations to inform biology and mechanisms. With the advent of revolutionary genome editing technologies, one can now reach beyond correlations to conduct definitive hypothesis testing. This review focuses on key discoveries that have emerged during the path from single loci to genome-wide analyses, specifically in the context of hematopoietic transcriptional mechanisms.

The 2014 Nucleic Acids Research Database Issue and an updated NAR online Molecular Biology Database Collection

The 2014 Nucleic Acids Research Database Issue includes descriptions of 58 new molecular biology databases and recent updates to 123 databases previously featured in NAR or other journals. For convenience, the issue is now divided into eight sections that reflect major subject categories. Among the highlights of this issue are six databases of the transcription factor binding sites in various organisms and updates on such popular databases as CAZy, Database of Genomic Variants (DGV), dbGaP, DrugBank, KEGG, miRBase, Pfam, Reactome, SEED, TCDB and UniProt. There is a strong block of structural databases, which includes, among others, the new RNA Bricks database, updates on PDBe, PDBsum, ArchDB, Gene3D, ModBase, Nucleic Acid Database and the recently revived iPfam database. An update on the NCBI’s MMDB describes VAST+, an improved tool for protein structure comparison. Two articles highlight the development of the Structural Classification of Proteins (SCOP) database: one describes SCOPe, which automates assignment of new structures to the existing SCOP hierarchy; the other one describes the first version of SCOP2, with its more flexible approach to classifying protein structures. This issue also includes a collection of articles on bacterial taxonomy and metagenomics, which includes updates on the List of Prokaryotic Names with Standing in Nomenclature (LPSN), Ribosomal Database Project (RDP), the Silva/LTP project and several new metagenomics resources. The NAR online Molecular Biology Database Collection, http://www.oxfordjournals.org/nar/database/c/, has been expanded to 1552 databases. The entire Database Issue is freely available online on the Nucleic Acids Research website (http://nar.oxfordjournals.org/).