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.

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.

A multi-sensor array for detecting and analyzing nocturnal avian migration.

Avian migration has fascinated humans for centuries. Insights into the lives of migrant birds are often elusive; however, recent, standalone technological innovations have revolutionized our understanding of this complex biological phenomenon. A future challenge for following these highly mobile animals is the necessity of bringing multiple technologies together to capture a more complete understanding of their movements. Here, we designed a proof-of-concept multi-sensor array consisting of two weather surveillance radars (WSRs), one local and one regional, an autonomous moon-watching sensor capable of detecting birds flying in front of the moon, and an autonomous recording unit (ARU) capable of recording avian nocturnal flight calls. We deployed this array at a field site in central Oklahoma on select nights in March, April, and May of 2021 and integrated data from this array with wind data corresponding to this site to examine the influence of wind on the movements of spring migrants aloft across these spring nights. We found that regional avian migration intensity is statistically significantly negatively correlated with wind velocity, in line with previous research. Furthermore, we found evidence suggesting that when faced with strong, southerly winds, migrants take advantage of these conditions by adjusting their flight direction by drifting. Importantly, we found that most of the migration intensities detected by the sensors were intercorrelated, except when this correlation could not be ascertained because we lacked the sample size to do so. This study demonstrates the potential for multi-sensor arrays to reveal the detailed ways in which avian migrants move in response to changing atmospheric conditions while in flight.

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 .

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

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.

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.

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.

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.

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.

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.