"Stripe" transcription factors provide accessibility to co-binding partners in mammalian genomes.

Regulatory elements activate promoters by recruiting transcription factors (TFs) to specific motifs. Notably, TF-DNA interactions often depend on cooperativity with colocalized partners, suggesting an underlying cis-regulatory syntax. To explore TF cooperativity in mammals, we analyze ∼500 mouse and human primary cells by combining an atlas of TF motifs, footprints, ChIP-seq, transcriptomes, and accessibility. We uncover two TF groups that colocalize with most expressed factors, forming stripes in hierarchical clustering maps. The first group includes lineage-determining factors that occupy DNA elements broadly, consistent with their key role in tissue-specific transcription. The second one, dubbed universal stripe factors (USFs), comprises ∼30 SP, KLF, EGR, and ZBTB family members that recognize overlapping GC-rich sequences in all tissues analyzed. Knockouts and single-molecule tracking reveal that USFs impart accessibility to colocalized partners and increase their residence time. Mammalian cells have thus evolved a TF superfamily with overlapping DNA binding that facilitate chromatin accessibility.

Single-Cell Analyses Reveal Megakaryocyte-Biased Hematopoiesis in Myelofibrosis and Identify Mutant Clone-Specific Targets.

Myelofibrosis is a severe myeloproliferative neoplasm characterized by increased numbers of abnormal bone marrow megakaryocytes that induce fibrosis, destroying the hematopoietic microenvironment. To determine the cellular and molecular basis for aberrant megakaryopoiesis in myelofibrosis, we performed single-cell transcriptome profiling of 135,929 CD34+ lineage- hematopoietic stem and progenitor cells (HSPCs), single-cell proteomics, genomics, and functional assays. We identified a bias toward megakaryocyte differentiation apparent from early multipotent stem cells in myelofibrosis and associated aberrant molecular signatures. A sub-fraction of myelofibrosis megakaryocyte progenitors (MkPs) are transcriptionally similar to healthy-donor MkPs, but the majority are disease specific, with distinct populations expressing fibrosis- and proliferation-associated genes. Mutant-clone HSPCs have increased expression of megakaryocyte-associated genes compared to wild-type HSPCs, and we provide early validation of G6B as a potential immunotherapy target. Our study paves the way for selective targeting of the myelofibrosis clone and illustrates the power of single-cell multi-omics to discover tumor-specific therapeutic targets and mediators of tissue fibrosis.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes.

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state Regulatory Potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbored distinctive transcription factor binding motifs that were similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we showed that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

MRP5 and MRP9 play a concerted role in male reproduction and mitochondrial function.

Multidrug Resistance Proteins (MRPs) are transporters that play critical roles in cancer even though the physiological substrates of these enigmatic transporters are poorly elucidated. In Caenorhabditis elegans, MRP5/ABCC5 is an essential heme exporter because mrp-5 mutants are unviable due to their inability to export heme from the intestine to extraintestinal tissues. Heme supplementation restores viability of these mutants but fails to restore male reproductive deficits. Correspondingly, cell biological studies show that MRP5 regulates heme levels in the mammalian secretory pathway even though MRP5 knockout (KO) mice do not show reproductive phenotypes. The closest homolog of MRP5 is MRP9/ABCC12, which is absent in C. elegans, raising the possibility that MRP9 may genetically compensate for MRP5. Here, we show that MRP5 and MRP9 double KO (DKO) mice are viable but reveal significant male reproductive deficits. Although MRP9 is highly expressed in sperm, MRP9 KO mice show reproductive phenotypes only when MRP5 is absent. Both ABCC transporters localize to mitochondrial-associated membranes, dynamic scaffolds that associate the mitochondria and endoplasmic reticulum. Consequently, DKO mice reveal abnormal sperm mitochondria with reduced mitochondrial membrane potential and fertilization rates. Metabolomics show striking differences in metabolite profiles in the DKO testes, and RNA sequencing shows significant alterations in genes related to mitochondrial function and retinoic acid metabolism. Targeted functional metabolomics reveal lower retinoic acid levels in the DKO testes and higher levels of triglycerides in the mitochondria. These findings establish a model in which MRP5 and MRP9 play a concerted role in regulating male reproductive functions and mitochondrial sufficiency.

Snapshot: a package for clustering and visualizing epigenetic history during cell differentiation.

BACKGROUND: Epigenetic modification of chromatin plays a pivotal role in regulating gene expression during cell differentiation. The scale and complexity of epigenetic data pose significant challenges for biologists to identify the regulatory events controlling cell differentiation. RESULTS: To reduce the complexity, we developed a package, called Snapshot, for clustering and visualizing candidate cis-regulatory elements (cCREs) based on their epigenetic signals during cell differentiation. This package first introduces a binarized indexing strategy for clustering the cCREs. It then provides a series of easily interpretable figures for visualizing the signal and epigenetic state patterns of the cCREs clusters during the cell differentiation. It can also use different hierarchies of cell types to highlight the epigenetic history specific to any particular cell lineage. We demonstrate the utility of Snapshot using data from a consortium project for ValIdated Systematic IntegratiON (VISION) of epigenomic data in hematopoiesis. CONCLUSION: The package Snapshot can identify all distinct clusters of genomic locations with unique epigenetic signal patterns during cell differentiation. It outperforms other methods in terms of interpreting and reproducing the identified cCREs clusters. The package of Snapshot is available at GitHub: https://github.com/guanjue/Snapshot .

An integrative view of the regulatory and transcriptional landscapes in mouse hematopoiesis.

Thousands of epigenomic data sets have been generated in the past decade, but it is difficult for researchers to effectively use all the data relevant to their projects. Systematic integrative analysis can help meet this need, and the VISION project was established for validated systematic integration of epigenomic data in hematopoiesis. Here, we systematically integrated extensive data recording epigenetic features and transcriptomes from many sources, including individual laboratories and consortia, to produce a comprehensive view of the regulatory landscape of differentiating hematopoietic cell types in mouse. By using IDEAS as our integrative and discriminative epigenome annotation system, we identified and assigned epigenetic states simultaneously along chromosomes and across cell types, precisely and comprehensively. Combining nuclease accessibility and epigenetic states produced a set of more than 200,000 candidate cis-regulatory elements (cCREs) that efficiently capture enhancers and promoters. The transitions in epigenetic states of these cCREs across cell types provided insights into mechanisms of regulation, including decreases in numbers of active cCREs during differentiation of most lineages, transitions from poised to active or inactive states, and shifts in nuclease accessibility of CTCF-bound elements. Regression modeling of epigenetic states at cCREs and gene expression produced a versatile resource to improve selection of cCREs potentially regulating target genes. These resources are available from our VISION website to aid research in genomics and hematopoiesis.

Regulation of erythrocyte Na+/K+/2Cl- cotransport by an oxygen-switched kinase cascade.

Many erythrocyte processes and pathways, including glycolysis, the pentose phosphate pathway (PPP), KCl cotransport, ATP release, Na+/K+-ATPase activity, ankyrin-band 3 interactions, and nitric oxide (NO) release, are regulated by changes in O2 pressure that occur as a red blood cell (RBC) transits between the lungs and tissues. The O2 dependence of glycolysis, PPP, and ankyrin-band 3 interactions (affecting RBC rheology) are controlled by O2-dependent competition between deoxyhemoglobin (deoxyHb), but not oxyhemoglobin (oxyHb), and other proteins for band 3. We undertook the present study to determine whether the O2 dependence of Na+/K+/2Cl- cotransport (catalyzed by Na+/K+/2Cl- cotransporter 1 [NKCC1]) might similarly originate from competition between deoxyHb and a protein involved in NKCC1 regulation for a common binding site on band 3. Using three transgenic mouse strains having mutated deoxyhemoglobin-binding sites on band 3, we found that docking of deoxyhemoglobin at the N terminus of band 3 displaces the protein with no lysine kinase 1 (WNK1) from its overlapping binding site on band 3. This displacement enabled WNK1 to phosphorylate oxidative stress-responsive kinase 1 (OSR1), which, in turn, phosphorylated and activated NKCC1. Under normal solution conditions, the NKCC1 activation increased RBC volume and thereby induced changes in RBC rheology. Because the deoxyhemoglobin-mediated WNK1 displacement from band 3 in this O2 regulation pathway may also occur in the regulation of other O2-regulated ion transporters, we hypothesize that the NKCC1-mediated regulatory mechanism may represent a general pattern of O2 modulation of ion transporters in erythrocytes.

Systematic integration of GATA transcription factors and epigenomes via IDEAS paints the regulatory landscape of hematopoietic cells.

Members of the GATA family of transcription factors play key roles in the differentiation of specific cell lineages by regulating the expression of target genes. Three GATA factors play distinct roles in hematopoietic differentiation. In order to better understand how these GATA factors function to regulate genes throughout the genome, we are studying the epigenomic and transcriptional landscapes of hematopoietic cells in a model-driven, integrative fashion. We have formed the collaborative multi-lab VISION project to conduct ValIdated Systematic IntegratiON of epigenomic data in mouse and human hematopoiesis. The epigenomic data included nuclease accessibility in chromatin, CTCF occupancy, and histone H3 modifications for 20 cell types covering hematopoietic stem cells, multilineage progenitor cells, and mature cells across the blood cell lineages of mouse. The analysis used the Integrative and Discriminative Epigenome Annotation System (IDEAS), which learns all common combinations of features (epigenetic states) simultaneously in two dimensions-along chromosomes and across cell types. The result is a segmentation that effectively paints the regulatory landscape in readily interpretable views, revealing constitutively active or silent loci as well as the loci specifically induced or repressed in each stage and lineage. Nuclease accessible DNA segments in active chromatin states were designated candidate cis-regulatory elements in each cell type, providing one of the most comprehensive registries of candidate hematopoietic regulatory elements to date. Applications of VISION resources are illustrated for the regulation of genes encoding GATA1, GATA2, GATA3, and Ikaros. VISION resources are freely available from our website http://usevision.org.

Defective erythropoiesis caused by mutations of the thyroid hormone receptor α gene.

Patients with mutations of the THRA gene exhibit classical features of hypothyroidism, including erythroid disorders. We previously created a mutant mouse expressing a mutated TRα1 (denoted as PV; Thra1PV/+ mouse) that faithfully reproduces the classical hypothyroidism seen in patients. Using Thra1PV/+ mice, we explored how the TRα1PV mutant acted to cause abnormalities in erythropoiesis. Thra1PV/+ mice exhibited abnormal red blood cell indices similarly as reported for patients. The total bone marrow cells and erythrocytic progenitors were markedly reduced in the bone marrow of Thra1PV/+ mice. In vitro terminal differentiation assays showed a significant reduction of mature erythrocytes in Thra1PV/+ mice. In wild-type mice, the clonogenic potential of progenitors in the erythrocytic lineage was stimulated by thyroid hormone (T3), suggesting that T3 could directly accelerate the differentiation of progenitors to mature erythrocytes. Analysis of gene expression profiles showed that the key regulator of erythropoiesis, the Gata-1 gene, and its regulated genes, such as the Klf1, ß-globin, dematin genes, CAII, band3 and eALAS genes, involved in the maturation of erythrocytes, was decreased in the bone marrow cells of Thra1PV/+ mice. We further elucidated that the Gata-1 gene was a T3-directly regulated gene and that TRα1PV could impair erythropoiesis via repression of the Gata-1 gene and its regulated genes. These results provide new insights into how TRα1 mutants acted to cause erythroid abnormalities in patients with mutations of the THRA gene. Importantly, the Thra1PV/+ mouse could serve as a preclinical mouse model to identify novel molecular targets for treatment of erythroid disorders.

Molecular convergence in ex vivo models of Diamond-Blackfan anemia.

Diamond-Blackfan anemia (DBA) is a congenital bone marrow failure syndrome characterized by erythroid hypoplasia, usually without perturbation of other hematopoietic lineages. Approximately 65% of DBA patients with autosomal dominant inheritance have heterozygous mutations or deletions in ribosomal protein (RP) genes while <1% of patients with X-linked inheritance have been identified with mutations in the transcription factor GATA1 Erythroid cells from patients with DBA have not been well characterized, and the mechanisms underlying the erythroid specific effects of either RP or GATA1 associated DBA remain unclear. We have developed an ex vivo culture system to expand peripheral blood CD34+ progenitor cells from patients with DBA and differentiate them into erythroid cells. Cells from patients with RP or GATA1 mutations showed decreased proliferation and delayed erythroid differentiation in comparison with controls. RNA transcript analyses of erythroid cells from controls and patients with RP or GATA1 mutations showed distinctive differences, with upregulation of heme biosynthesis genes prominently in RP-mediated DBA and failure to upregulate components of the translational apparatus in GATA1-mediated DBA. Our data show that dysregulation of translation is a common feature of DBA caused by both RP and GATA1 mutations. This trial was registered at www.clinicaltrials.gov as #NCT00106015.

SigSeeker: a peak-calling ensemble approach for constructing epigenetic signatures.

MOTIVATION: Epigenetic data are invaluable when determining the regulatory programs governing a cell. Based on use of next-generation sequencing data for characterizing epigenetic marks and transcription factor binding, numerous peak-calling approaches have been developed to determine sites of genomic significance in these data. Such analyses can produce a large number of false positive predictions, suggesting that sites supported by multiple algorithms provide a stronger foundation for inferring and characterizing regulatory programs associated with the epigenetic data. Few methodologies integrate epigenetic based predictions of multiple approaches when combining profiles generated by different tools. RESULTS: The SigSeeker peak-calling ensemble uses multiple tools to identify peaks, and with user-defined thresholds for peak overlap and signal strength it retains only those peaks that are concordant across multiple tools. Peaks predicted to be co-localized by only a very small number of tools, discovered to be only marginally overlapping, or found to represent significant outliers to the approximation model are removed from the results, providing concise and high quality epigenetic datasets. SigSeeker has been validated using established benchmarks for transcription factor binding and histone modification ChIP-Seq data. These comparisons indicate that the quality of our ensemble technique exceeds that of single tool approaches, enhances existing peak-calling ensembles, and results in epigenetic profiles of higher confidence. AVAILABILITY AND IMPLEMENTATION: http://sigseeker.org. CONTACT: lichtenbergj@mail.nih.gov. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Reversible binding of hemoglobin to band 3 constitutes the molecular switch that mediates O2 regulation of erythrocyte properties.

Functional studies have shown that the oxygenation state of the erythrocyte regulates many important pathways, including glucose metabolism, membrane mechanical stability, and cellular adenosine triphosphate (ATP) release. Deoxyhemoglobin (deoxyHb), but not oxyhemoglobin, binds avidly and reversibly to band 3, the major erythrocyte membrane protein. Because band 3 associates with multiple metabolic, solute transport, signal transduction, and structural proteins, the hypothesis naturally arises that the O2-dependent regulation of erythrocyte properties might be mediated by the reversible association of deoxyHb with band 3. To explore whether the band 3-deoxyHb interaction constitutes a "molecular switch" for regulating erythrocyte biology, we have generated transgenic mice with mutations in the deoxyHb-binding domain of band 3. One strain of mouse contains a "humanized" band 3 in which the N-terminal 45 residues of mouse band 3 are replaced by the homologous sequence from human band 3, including the normal human band 3 deoxyHb-binding site. The second mouse contains the same substitution as the first, except the deoxyHb site on band 3 (residues 12-23) has been deleted. Comparison of these animals with wild-type mice demonstrates that the following erythrocyte properties are controlled by the O2-dependent association of hemoglobin with band 3: (1) assembly of a glycolytic enzyme complex on the erythrocyte membrane which is associated with a shift in glucose metabolism between the pentose phosphate pathway and glycolysis, (2) interaction of ankyrin with band 3 and the concomitant regulation of erythrocyte membrane stability, and (3) release of ATP from the red cell which has been linked to vasodilation.

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.

The COP9 signalosome is vital for timely repair of DNA double-strand breaks.

The DNA damage response is vigorously activated by DNA double-strand breaks (DSBs). The chief mobilizer of the DSB response is the ATM protein kinase. We discovered that the COP9 signalosome (CSN) is a crucial player in the DSB response and an ATM target. CSN is a protein complex that regulates the activity of cullin ring ubiquitin ligase (CRL) complexes by removing the ubiquitin-like protein, NEDD8, from their cullin scaffold. We find that the CSN is physically recruited to DSB sites in a neddylation-dependent manner, and is required for timely repair of DSBs, affecting the balance between the two major DSB repair pathways-nonhomologous end-joining and homologous recombination repair (HRR). The CSN is essential for the processivity of deep end-resection-the initial step in HRR. Cullin 4a (CUL4A) is recruited to DSB sites in a CSN- and neddylation-dependent manner, suggesting that CSN partners with CRL4 in this pathway. Furthermore, we found that ATM-mediated phosphorylation of CSN subunit 3 on S410 is critical for proper DSB repair, and that loss of this phosphorylation site alone is sufficient to cause a DDR deficiency phenotype in the mouse. This novel branch of the DSB response thus significantly affects genome stability.

Lineage and species-specific long noncoding RNAs during erythro-megakaryocytic development.

Mammals express thousands of long noncoding (lnc) RNAs, a few of which are known to function in tissue development. However, the entire repertoire of lncRNAs in most tissues and species is not defined. Indeed, most lncRNAs are not conserved, raising questions about function. We used RNA sequencing to identify 1109 polyadenylated lncRNAs expressed in erythroblasts, megakaryocytes, and megakaryocyte-erythroid precursors of mice, and 594 in erythroblasts of humans. More than half of these lncRNAs were unannotated, emphasizing the opportunity for new discovery through studies of specialized cell types. Analysis of the mouse erythro-megakaryocytic polyadenylated lncRNA transcriptome indicates that ~75% arise from promoters and 25% from enhancers, many of which are regulated by key transcription factors including GATA1 and TAL1. Erythroid lncRNA expression is largely conserved among 8 different mouse strains, yet only 15% of mouse lncRNAs are expressed in humans and vice versa, reflecting dramatic species-specificity. RNA interference assays of 21 abundant erythroid-specific murine lncRNAs in primary mouse erythroid precursors identified 7 whose knockdown inhibited terminal erythroid maturation. At least 6 of these 7 functional lncRNAs have no detectable expression in human erythroblasts, suggesting that lack of conservation between mammalian species does not predict lack of function.

Functional Niche Competition Between Normal Hematopoietic Stem and Progenitor Cells and Myeloid Leukemia Cells.

Hematopoietic stem and progenitor cells (HSPCs) reside in a specialized niche that regulates their proliferative capacity and their fate. There is increasing evidence for similar roles of marrow niches on controlling the behavior of leukemic cells; however, whether normal hematopoietic stem cell (HSC) and leukemic cells reside in or functionally compete for the same marrow niche is unclear. We used the mixed lineage leukemia-AF9 (MLL-AF9) murine acute myeloid leukemia (AML) in a competitive repopulation model to investigate whether normal HSPC and leukemic cells functionally compete for the same marrow niches. Irradiated recipient mice were transplanted with fixed numbers of MLL-AF9 cells mixed with increasing doses of normal syngeneic whole bone marrow (WBM) or with purified HSPC (LSK). Survival was significantly increased and leukemic progression was delayed proportional to increasing doses of normal WBM or normal LSK cells in multiple independent experiments, with all doses of WBM or LSK cells studied above the threshold for rapid and complete hematopoietic reconstitution in the absence of leukemia. Confocal microscopy demonstrated nests of either leukemic cells or normal hematopoietic cells but not both in the marrow adjacent to endosteum. Early following transplantation, leukemic cells from animals receiving lower LSK doses were cycling more actively than in those receiving higher doses. These results suggest that normal HSPC and AML cells compete for the same functional niche. Manipulation of the niche could impact on response to antileukemic therapies, and the numbers of normal HSPC could impact on leukemia outcome, informing approaches to cell dose in the context of stem cell transplantation.

Genome-wide DNA methylation profiles in hematopoietic stem and progenitor cells reveal overrepresentation of ETS transcription factor binding sites.

DNA methylation is an essential epigenetic mark that is required for normal development. Knockout of the DNA methyltransferase enzymes in the mouse hematopoietic compartment reveals that methylation is critical for hematopoietic differentiation. To better understand the role of DNA methylation in hematopoiesis, we characterized genome-wide DNA methylation in primary mouse hematopoietic stem cells (HSCs), common myeloid progenitors (CMPs), and erythroblasts (ERYs). Methyl binding domain protein 2 (MBD) enrichment of DNA followed by massively parallel sequencing (MBD-seq) was used to map genome-wide DNA methylation. Globally, DNA methylation was most abundant in HSCs, with a 40% reduction in CMPs, and a 67% reduction in ERYs. Only 3% of peaks arise during differentiation, demonstrating a genome-wide decline in DNA methylation during erythroid development. Analysis of genomic features revealed that 98% of promoter CpG islands are hypomethylated, while 20%-25% of non-promoter CpG islands are methylated. Proximal promoter sequences of expressed genes are hypomethylated in all cell types, while gene body methylation positively correlates with gene expression in HSCs and CMPs. Elevated genome-wide DNA methylation in HSCs and the positive association between methylation and gene expression demonstrates that DNA methylation is a mark of cellular plasticity in HSCs. Using de novo motif discovery, we identified overrepresented transcription factor consensus binding motifs in methylated sequences. Motifs for several ETS transcription factors, including GABPA and ELF1, are overrepresented in methylated regions. Our genome-wide survey demonstrates that DNA methylation is markedly altered during myeloid differentiation and identifies critical regions of the genome and transcription factor programs that contribute to hematopoiesis.