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 V al i dated S ystematic I ntegrati on (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.

"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.

Transitions in lineage specification and gene regulatory networks in hematopoietic stem/progenitor cells over human development.

Human hematopoiesis is a dynamic process that starts in utero 18-21 days post-conception. Understanding the site- and stage-specific variation in hematopoiesis is important if we are to understand the origin of hematological disorders, many of which occur at specific points in the human lifespan. To unravel how the hematopoietic stem/progenitor cell (HSPC) compartment changes during human ontogeny and the underlying gene regulatory mechanisms, we compare 57,489 HSPCs from 5 different tissues spanning 4 developmental stages through the human lifetime. Single-cell transcriptomic analysis identifies significant site- and developmental stage-specific transitions in cellular architecture and gene regulatory networks. Hematopoietic stem cells show progression from cycling to quiescence and increased inflammatory signaling during ontogeny. We demonstrate the utility of this dataset for understanding aberrant hematopoiesis through comparison to two cancers that present at distinct time points in postnatal life-juvenile myelomonocytic leukemia, a childhood cancer, and myelofibrosis, which classically presents in older adults.

Single-cell profiling of human bone marrow progenitors reveals mechanisms of failing erythropoiesis in Diamond-Blackfan anemia.

Ribosome dysfunction underlies the pathogenesis of many cancers and heritable ribosomopathies. Here, we investigate how mutations in either ribosomal protein large (RPL) or ribosomal protein small (RPS) subunit genes selectively affect erythroid progenitor development and clinical phenotypes in Diamond-Blackfan anemia (DBA), a rare ribosomopathy with limited therapeutic options. Using single-cell assays of patient-derived bone marrow, we delineated two distinct cellular trajectories segregating with ribosomal protein genotypes. Almost complete loss of erythroid specification was observed in RPS-DBA. In contrast, we observed relative preservation of qualitatively abnormal erythroid progenitors and precursors in RPL-DBA. Although both DBA genotypes exhibited a proinflammatory bone marrow milieu, RPS-DBA was characterized by erythroid differentiation arrest, whereas RPL-DBA was characterized by preserved GATA1 expression and activity. Compensatory stress erythropoiesis in RPL-DBA exhibited disordered differentiation underpinned by an altered glucocorticoid molecular signature, including reduced ZFP36L2 expression, leading to milder anemia and improved corticosteroid response. This integrative analysis approach identified distinct pathways of erythroid failure and defined genotype-phenotype correlations in DBA. These findings may help facilitate therapeutic target discovery.

Effects of sheared chromatin length on ChIP-seq quality and sensitivity.

Chromatin immunoprecipitation followed by massively parallel, high throughput sequencing (ChIP-seq) is the method of choice for genome-wide identification of DNA segments bound by specific transcription factors or in chromatin with particular histone modifications. However, the quality of ChIP-seq datasets varies widely, with a substantial fraction being of intermediate to poor quality. Thus, it is important to discern and control the factors that contribute to variation in ChIP-seq. In this study, we focused on sonication, a user-controlled variable, to produce sheared chromatin. We systematically varied the amount of shearing of fixed chromatin from a mouse erythroid cell line, carefully measuring the distribution of resultant fragment lengths prior to ChIP-seq. This systematic study was complemented with a retrospective analysis of additional experiments. We found that the level of sonication had a pronounced impact on the quality of ChIP-seq signals. Over-sonication consistently reduced quality, while the impact of under-sonication differed among transcription factors, with no impact on sites bound by CTCF but frequently leading to the loss of sites occupied by TAL1 or bound by POL2. The bound sites not observed in low-quality datasets were inferred to be a mix of both direct and indirect binding. We leveraged these findings to produce a set of CTCF ChIP-seq datasets in rare, primary hematopoietic progenitor cells. Our observation that the amount of chromatin sonication is a key variable in success of ChIP-seq experiments indicates that monitoring the level of sonication can improve ChIP-seq quality and reproducibility and facilitate ChIP-seq in rare cell types.

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.

Normal Iron Homeostasis Requires the Transporter SLC48A1 for Efficient Heme-Iron Recycling in Mammals.

In mammals over 65% of the total body iron is located within erythrocytes in the heme moieties of hemoglobin. Iron homeostasis requires iron absorbed from the diet by the gut as well as recycling of iron after the destruction of senescent erythrocytes. Senescent erythrocytes are engulfed by reticuloendothelial system macrophages where hemoglobin is broken down in the lysosomes, releasing heme for iron recovery in the cytoplasm. We recently showed that the SLC48A1 protein is responsible for transporting heme from the lysosome to the cytoplasm. CRISPR generated SLC48A1-deficient mice accumulate heme in their reticuloendothelial system macrophages as hemozoin crystals. Here we describe additional features of SLC48A1-deficient mice. We show that visible hemozoin first appears in the reticuloendothelial system macrophages of SLC48A1-deficient mice at 8 days of age, indicating the onset of erythrocyte recycling. Evaluation of normal and SLC48A1-deficient mice on iron-controlled diets show that SLC48A1-mediated iron recycling is equivalent to at least 10 parts per million of dietary iron. We propose that mutations in human SLC48A1 could contribute to idiopathic iron disorders.

Hemozoin produced by mammals confers heme tolerance.

Free heme is cytotoxic as exemplified by hemolytic diseases and genetic deficiencies in heme recycling and detoxifying pathways. Thus, intracellular accumulation of heme has not been observed in mammalian cells to date. Here we show that mice deficient for the heme transporter SLC48A1 (also known as HRG1) accumulate over ten-fold excess heme in reticuloendothelial macrophage lysosomes that are 10 to 100 times larger than normal. Macrophages tolerate these high concentrations of heme by crystallizing them into hemozoin, which heretofore has only been found in blood-feeding organisms. SLC48A1 deficiency results in impaired erythroid maturation and an inability to systemically respond to iron deficiency. Complete heme tolerance requires a fully-operational heme degradation pathway as haplo insufficiency of HMOX1 combined with SLC48A1 inactivation causes perinatal lethality demonstrating synthetic lethal interactions between heme transport and degradation. Our studies establish the formation of hemozoin by mammals as a previously unsuspected heme tolerance pathway.

Specialized cells, known as red blood cells, are responsible for transporting oxygen to various organs in the body. Each red blood cell contains over a billion molecules of heme which make up the iron containing portion of the hemoglobin protein that binds and transports oxygen. When red blood cells reach the end of their life, they are degraded, and the heme and iron inside them is recycled to produce new red blood cells. Heme, however, is highly toxic to cells, and can cause severe tissue damage if not properly removed. Scavenger cells called macrophages perform this recycling role in the spleen, liver and bone marrow. Collectively, macrophages can process around five million red blood cells every second or about 100 trillion heme molecules. But, it is unclear how they are able to handle such enormous volumes. Macrophages isolated from human and mice have been shown to transport heme from damaged red blood cells using a protein called HRG1. To investigate the role HRG1 plays in heme-iron recycling, Pek et al. used a gene editing tool known an CRISPR/Cas9 to remove the gene for HRG1 from the macrophages of mice. If HRG1 is a major part of this process, removing the gene should result in a build-up of toxic heme and eventual death of the mouse. But, rather than dying of heme-iron overload as expected, these mutant mice managed to survive. Pek et al. found that despite being unable to recycle heme, these mice were still able to make new red blood cells as long as they had a diet that was rich in iron. However, the darkening color of the spleen, bone marrow, and liver in these HRG1 deficient mice indicated that these mice were still accumulating high levels of heme. Further experiments revealed that these mice protected themselves from toxicity by converting the excess heme into crystals called hemozoin. This method of detoxification is commonly seen in blood-feeding parasites, and this is the first time it has been observed in a mammal. These crystals invite new questions about how mammals recycle heme and what happens when this process goes wrong. The next step is to ask whether humans also start to make hemozoin if the gene for HRG1 is faulty. If so, this could open a new avenue of exploration into treatments for red blood cell diseases like anemia and iron overload.

Oxygen tension-mediated erythrocyte membrane interactions regulate cerebral capillary hyperemia.

The tight coupling between cerebral blood flow and neural activity is a key feature of normal brain function and forms the basis of functional hyperemia. The mechanisms coupling neural activity to vascular responses, however, remain elusive despite decades of research. Recent studies have shown that cerebral functional hyperemia begins in capillaries, and red blood cells (RBCs) act as autonomous regulators of brain capillary perfusion. RBCs then respond to local changes of oxygen tension (PO2) and regulate their capillary velocity. Using ex vivo microfluidics and in vivo two-photon microscopy, we examined RBC capillary velocity as a function of PO2 and showed that deoxygenated hemoglobin and band 3 interactions on RBC membrane are the molecular switch that responds to local PO2 changes and controls RBC capillary velocity. Capillary hyperemia can be controlled by manipulating RBC properties independent of the neurovascular unit, providing an effective strategy to treat or prevent impaired functional hyperemia.

All hands on DEK.

In this issue of the JCI, Capitano et al. demonstrate that the secreted form of the DNA-binding chromatin factor DEK regulates hematopoiesis. It is known that DEK can be secreted by macrophages and other cells, but no function has been attached to secreted DEK. Capitano et al. showed that extracellular DEK activates signaling through the CXCL2 receptor, which in turn enhances the proliferation of hematopoietic stem cells and decreases hematopoietic progenitor cell numbers both in vivo and in vitro. These results offer the opportunity to expand transplantable stem cells to improve outcomes in patients undergoing bone marrow transplant.

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.

Establishment of regulatory elements during erythro-megakaryopoiesis identifies hematopoietic lineage-commitment points.

BACKGROUND: Enhancers and promoters are cis-acting regulatory elements associated with lineage-specific gene expression. Previous studies showed that different categories of active regulatory elements are in regions of open chromatin, and each category is associated with a specific subset of post-translationally marked histones. These regulatory elements are systematically activated and repressed to promote commitment of hematopoietic stem cells along separate differentiation paths, including the closely related erythrocyte (ERY) and megakaryocyte (MK) lineages. However, the order in which these decisions are made remains unclear. RESULTS: To characterize the order of cell fate decisions during hematopoiesis, we collected primary cells from mouse bone marrow and isolated 10 hematopoietic populations to generate transcriptomes and genome-wide maps of chromatin accessibility and histone H3 acetylated at lysine 27 binding (H3K27ac). Principle component analysis of transcriptional and open chromatin profiles demonstrated that cells of the megakaryocyte lineage group closely with multipotent progenitor populations, whereas erythroid cells form a separate group distinct from other populations. Using H3K27ac and open chromatin profiles, we showed that 89% of immature MK (iMK)-specific active regulatory regions are present in the most primitive hematopoietic cells, 46% of which contain active enhancer marks. These candidate active enhancers are enriched for transcription factor binding site motifs for megakaryopoiesis-essential proteins, including ERG and ETS1. In comparison, only 64% of ERY-specific active regulatory regions are present in the most primitive hematopoietic cells, 20% of which containing active enhancer marks. These regions were not enriched for any transcription factor consensus sequences. Incorporation of genome-wide DNA methylation identified significant levels of de novo methylation in iMK, but not ERY. CONCLUSIONS: Our results demonstrate that megakaryopoietic profiles are established early in hematopoiesis and are present in the majority of the hematopoietic progenitor population. However, megakaryopoiesis does not constitute a "default" differentiation pathway, as extensive de novo DNA methylation accompanies megakaryopoietic commitment. In contrast, erythropoietic profiles are not established until a later stage of hematopoiesis, and require more dramatic changes to the transcriptional and epigenetic programs. These data provide important insights into lineage commitment and can contribute to ongoing studies related to diseases associated with differentiation defects.

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.