- Open Access
Chromatin accessibility maps provide evidence of multilineage gene priming in hematopoietic stem cells
Epigenetics & Chromatin volume 14, Article number: 2 (2021)
Hematopoietic stem cells (HSCs) have the capacity to differentiate into vastly different types of mature blood cells. The epigenetic mechanisms regulating the multilineage ability, or multipotency, of HSCs are not well understood. To test the hypothesis that cis-regulatory elements that control fate decisions for all lineages are primed in HSCs, we used ATAC-seq to compare chromatin accessibility of HSCs with five unipotent cell types. We observed the highest similarity in accessibility profiles between megakaryocyte progenitors and HSCs, whereas B cells had the greatest number of regions with de novo gain in accessibility during differentiation. Despite these differences, we identified cis-regulatory elements from all lineages that displayed epigenetic priming in HSCs. These findings provide new insights into the regulation of stem cell multipotency, as well as a resource to identify functional drivers of lineage fate.
HSCs have higher global chromatin accessibility than any unilineage progeny
Megakaryocyte progenitors are the most closely related unipotent cell type to HSCs
B cell commitment involves de novo chromatin accessibility
Evidence of cis-element priming of lineage-specific genes in HSCs
Multipotency is a key feature of hematopoietic stem cells (HSCs) and essential for their ability to produce all types of blood and immune cells in situ and upon therapeutic stem cell transplantation. The mechanistic basis of multipotency is unclear, but previous studies have shown that the regulation of differentiation programs is achieved, in large part, through epigenetic remodeling of cis-regulatory elements (CREs) [17, 41, 46]. Thus, HSC multipotency may be enabled by accessible non-promoter CREs that keep loci competent for transcription factor binding and gene activation without active expression. Such selective “CRE priming” may underlie the developmental competence of specific cell types, which is then acted upon by inductive signals to gradually specify fate . When all CREs that drive differentiation and lineage choice are primed in stem cells, that stem cell is in a permissive state (Fig. 1a) and is competent to initiate differentiation into all mature lineages.
We sought to test two models of HSC multipotency that are based on regulation of chromatin organization: the “permissive fate model” and a “de novo activation model” (Fig. 1a). Supporting a role for the permissive model in stem cell lineage potential are observations of bivalent histone domains that maintain key developmental genes in embryonic stem cells (ESCs) poised for activation , and an overall accessible chromatin state in both ESCs and HSCs compared to lineage-restricted progenitors and mature cells [19, 20, 44]. When differentiation occurs, the genes poised for differentiation into the induced lineage are activated while CREs that would drive differentiation into alternative lineages are silenced. This has been observed in ESCs and during differentiation of ESCs into endoderm [46, 47]. Our observation of global chromatin condensation and localization of H3K9me3-marked repressed domains or heterochromatin towards the nuclear periphery during HSC differentiation also support the permissive model . Inversely, in the de novo activation model (Fig. 1a), CREs that drive lineage fate are inaccessible in HSCs. Differentiation and lineage choice occur by “unlocking” these CREs. Transcriptional and functional analyses of hematopoietic stem and progenitor cells support this de novo model, where lymphoid potential is gained in progenitor cells rather than being a consequence of CRE priming in HSCs [6, 12, 18, 32].
In order to interrogate these models and how they pertain to the regulation of competence in hematopoiesis, as well as gain a better understanding of the relationships between epigenetic, transcriptomic and functional observations, we mapped global chromatin accessibility using the Assay for Transposase Accessible Chromatin by High Throughput Sequencing (ATAC-seq) . This assay allows assessment of high resolution, genome-wide chromatin accessibility throughout differentiation programs of rare cells. The dynamics of chromatin accessibility in erythro-megakaryopoiesis  and granulocyte/macrophage development  have been highly informative. From these studies, the bulk observations gave us insight into the dynamics of lineage commitment during hematopoiesis, while single-cell analysis revealed the heterogeneity of epigenomic states and, therefore, lineage bias in progenitors throughout hematopoiesis. Based on those studies, as well as reports of global chromatin accessibility of embryonic [3, 10, 20] and hematopoietic [11, 44] stem cells, we hypothesized that HSCs are in a permissive chromatin state where CREs that control fate decisions are primed in HSCs. Here, we tested this hypothesis by performing in-depth ATAC-seq investigation of HSCs and 5 unipotent lineage cell populations representing the five main hematopoietic lineages (Fig. 1b), as defined by previously published phenotypes [4, 6].
Mapping of chromatin accessibility in HSCs and unipotent lineage cells identified a tight association of megakaryocyte progenitors to HSCs
To determine the dynamics of genome accessibility throughout hematopoiesis, we sorted six primary hematopoietic cell types (Fig. 1b) and performed ATAC-seq of libraries with expected fragment size distributions (Additional file 1: Figure S1) . We identified 70,731 peaks in HSCs, 47,363 peaks in megakaryocyte progenitors (MkPs), 38,007 in erythroid progenitors (EPs), 30,529 in granulocyte/macrophages (GMs), 70,358 in B cells, and 51,832 in T cells (Table 1). From these peak-lists we combined and filtered the peaks using the chromVAR package to only the most significant peaks, as defined by  and identified a total of 84,243 peaks, referred to as the master peak-list throughout the study (Table 1). To assess data quality, we analyzed replicate clustering and cell type relationships of all 6 cell types using principal component analysis and dimensionality reduction as a t-Distributed Stochastic Neighbor Embedding (tSNE) plot . All biological replicate samples closely associated with each other by tSNE analysis (Fig. 1c), as well as by hierarchical clustering using the chromVAR output (Fig. 1d). We observed two primary clusters in Fig. 1d: an HSC/MkP cluster and all other cell types. We also observed a distinct lymphoid cell subcluster containing only B and T cells, while GMs and EPs clustered independently. MkPs have the most similar accessibility to HSCs, with the ranking of the other cell types from most to least similar as EPs, GMs, Bs and then Ts. This is consistent with our tSNE analysis (Fig. 1c), where HSCs and MkPs closely associated with each other, and with studies that have reported a close relationship of HSCs with the megakaryocyte lineage [14, 37] and that erythropoiesis requires chromatin remodeling for differentiation to occur .
Visualization and comparison of ATAC-seq data generated in this study correlated with known expression patterns of cell type-specific genes
As another assessment of the quality and reproducibility of our ATAC-seq data, we used the Gene Expression Commons (GEXC) expression database  to generate a list of genes that were expressed only in each unipotent lineage cell type (Fig. 2a). From each list, we calculated the normalized average signal centered at the promoter of each cell-type-specific peak-list for each cell type by generating histograms using HOMER  (Fig. 2b). We observed the expected cell type-specific accessibility for each unipotent lineage with minimal signal from the other cell types. In addition, we visualized the ATAC-seq signals across promoters of some example genes with known cell type-specific expression patterns, plus a negative (expressed in none of the cell types) and a positive (expressed in all of the cell types) control: Gapdh (expressed in all cell types), Fezf2 (not expressed in any cell type), Ndn (expressed in HSCs only), Klf1 (EPs only), Gp6 (MkPs only), Ly6g (GMs only), CD19 (B cells only), and Ccr4 (T cells only) (Fig. 2c, d). Ly6g was not available in GEXC but is a well-known GM-selective gene . We observed the expected accessibility peaks in each cell type, as well as a minimal signal from cell types without expression of those genes (Fig. 2d). As an example of a well-characterized locus, we visualized our ATAC-seq data across the mouse β-globin cluster. As expected, we observed EP-selective accessibility of the hypersensitive sites in the locus control region (LCR) and of adult globin gene promoters β-major and β-minor [30, 34] (Additional file 2: Figure S2). The overall high level of reproducibility between independent sample replicates and clustering strategies (Fig. 1c, d), as well as the expected accessibility in cell type-specific genes (Fig. 2, Additional file 2: Figure S2), indicated that we had generated high-quality chromatin accessibility maps of these 6 cell types.
HSCs have greater global accessibility and undergo more extensive chromatin remodeling upon lymphoid differentiation
Using a number of quantitative, but non-sequence-specific assays, we previously reported that chromatin is progressively condensed upon HSC differentiation into unilineage and mature cells . To test whether the ATAC-seq data recapitulated these findings, we quantified the total number of distinct peaks, as well as the cumulative read-counts for all peaks, for each cell type. First, we took each cell type’s optimal peak-list from the Irreproducible Discovery Rate (IDR) analysis  and reported the number of peaks. We observed the highest number of peaks in HSCs (Fig. 3a), closely followed by B cells. In parallel, we quantified global accessibility by calculating the normalized average signal over the master peak-list for each cell type by generating histograms using HOMER . We observed similar ordering compared to the peak number, with HSCs having the highest average signal and B cells the second highest (Fig. 3b). The low signal in EPs is possibly due to widespread transcriptional silencing as the next step towards becoming highly specialized red blood cells and ejection of nuclei . Although these measurements are not completely independent, there is not a strict correlation between peak count and cumulative peak signal: for example, compared to EPs, GMs have fewer peaks (Fig. 3a) but higher cumulative readcount (Fig. 3b). Interestingly, HSCs displayed both the highest number of peaks and the greatest peak signal. These results are consistent with our previous findings of progressive chromatin condensation upon HSCs differentiation .
Comparisons of peaks gained and lost as HSCs differentiate into unilineage cells revealed an overall gain of accessibility selectively for B cell differentiation
To assess the number of peaks that changed upon HSCs differentiation, we took the IDR optimal peak-list for each cell type and performed pairwise comparisons between HSCs and the five mature/unipotent cell types (Fig. 3c). We quantified the number of peaks gained and lost by the unipotent progenitors/mature cells compared to HSCs (Fig. 3c–g). MkPs had the lowest number of peak changes (peaks gained plus lost; Fig. 3d), and therefore have the greatest proportion of peaks in common with HSCs. This was primarily driven by the low percentage of peaks gained (Fig. 3e), as opposed to peaks lost (Fig. 3f) upon HSC differentiation into MkPs. In contrast, EPs had the highest percentage of total peaks changed (Fig. 3d) due to the greatest percentage of peaks lost (Fig. 3f). This could be driven by EPs starting to shut down transcription to become highly specialized and eject their nuclei, reflected by the overall low accessibility observed (Fig. 3a, b). B cells had the highest percentage of peaks gained and the lowest percentage of peaks lost compared to the other cell types (Fig. 3e, f) and was the only cell type where the percentage of peaks gained was higher than peaks lost (Fig. 3g). This suggests that B cell fate requires chromatin remodeling to open up sites that drive B cell lineage fate.
Exclusively shared peaks between HSCs and unipotent cell types are primarily non-promoter and are enriched for known cell-type-specific transcription factors
We then turned our attention from peaks that were different between HSCs and their progeny to instead focus on elements with shared accessibility. We hypothesized that peaks that are exclusively shared between HSCs and one unipotent cell type contain elements that drive lineage commitment into that cell type. We filtered the peak-lists of all 6 cell types against each other using the HOMER mergePeaks.pl tool and annotated the peak-lists that each of unipotent lineage cell types exclusively shared with HSCs (Fig. 4a). We quantified the percentage of peaks that each unipotent cell type shared with HSCs (Fig. 4b). Consistent with the clustering profiles (Fig. 1c, d), MkPs had the highest percentage of peaks that were shared exclusively with HSCs. This similarity appeared to be primarily manifested in non-promoter elements: we annotated the exclusively shared peaks and categorized them as promoter or non-promoter peaks (Fig. 4c) and compared the distributions to the annotated peak-lists for each cell type assayed (Table 1). All of the exclusively shared peak-lists had significant enrichment (p-value < 0.001) of non-promoter peaks compared to the normal distribution of peaks in our dataset. Thus, non-promoter elements were shared between HSCs and their progeny significantly more frequently than promoter elements, especially with MkPs. Many, but likely not all, of these non-promoter accessible sites may serve as enhancers: about one-third of the non-promoter peaks overlapped with an enhancer catalog generated from chromatin immunoprecipitation (ChIP) experiments in blood cells  (Additional file 3: Figure S3A). Similar levels of overlap was observed between the ATAC-accessible peaks in our ATAC exclusively shared peak-lists with H3K4me1 modifications in HSCs, while less overlap was observed for H3K27Ac, at the aggregate and cell type-specific level (Additional file 3: Figure S3B, C).
To determine what transcription factor binding sites were present within the exclusively shared peaks, we performed motif enrichment using the HOMER package and reported the top 10 results for each cell type, sorted by p-value (Fig. 4d–h). The peaks that HSCs shared with MkPs (Fig. 4d) or EPs (Fig. 4e) were primarily enriched for Gata family transcription factors and their inhibitor TRPS1. Notably, HSC/MkP peaks also had enrichment of ERG and Runx1, which are known drivers of hematopoiesis [21, 26]. For HSC/EPs, Gata1 was the most enriched motif, with the Gata:SCL combination motif and NF-E2 and NFE2L motifs also scoring in the top ten. These factors are all known to be important in red blood cell differentiation, and NF-E2 is known to regulate SCL and Gata2 . HSC/GM peaks had enrichment of known regulators of GM cell fate, such as CEBP, PU.1, and SpiB (Fig. 4f). HSC/B cells primarily had CTCF and CTCFL motif enrichment (Fig. 4g). These motifs could be a reason for the overall high number of peaks observed in B cells (Fig. 3a, b), as 44.7% and 46.6% of the shared peaks contained CTCF or CTCFL motifs, respectively. HSC/T cell peaks were enriched for Tcf and Tbx family factors that are known to play a role in T cell development (Fig. 4h). Overall, all five HSC-shared peak-lists had enrichment of transcription factors that are known to be important for normal differentiation for each lineage.
Evidence of cis-element priming of lineage-specific genes in HSCs
Previous work on understanding multipotency and developmental competence suggests a model where competence is conferred by transcriptional priming: being competent of transcription factor binding and gene expression, without active expression . One of the suggested regulators of transcriptional priming are non-promoter cis-regulatory elements (CREs). This means that CREs that drive lineage fate for all lineages are accessible in HSCs in our permissive fate model and inaccessible in our de novo activation model. We hypothesized that CREs that are exclusively shared between HSCs and a unipotent lineage cell are potential drivers of that lineage. We utilized the GREAT tool  to annotate and predict the target genes for each exclusively shared CRE. Here we report examples of genes and a predicted CRE for each lineage that is primed in HSCs. In addition, we linked the motif enrichment with the GREAT analysis by annotating the CREs using the top 10 motifs enriched by p-value (Fig. 4d–h) for each exclusive HSC/unipotent cell type. In MkPs, a predicted CRE for Thrombin receptor like 2 (F2rl2) was found. This gene is expressed only in MkPs (Fig. 5a), while the CRE is only accessible in HSCs and MkPs (Fig. 5b). This CRE contained 9 out of the top 10 motifs, with the Runx1 motif being the only one missing (Fig. 5c). Pyruvate kinase liver and red blood cell (Pklr) was found to be expressed only in EPs (Fig. 5d), and a predicted CRE was accessible only in HSCs and EPs (Fig. 5e). Motifs for Gata2, Gata3, Gata4, and TRPS1 were found within the CRE (Fig. 5f). In GMs, Mitochondrial tumor suppressor 1 (Mtus1) was found to be primed in HSCs, with expression only in GMs (Fig. 5g), accessibility of a predicted CRE only in HSCs and GMs (Fig. 5h), and the presence of transcription factors known to play a role in GM development, such as CEBP and PU.1 (Fig. 5i). In B cells, Interferon regulatory factor 8 (Irf8), is only expressed in B cells (Fig. 5J), the predicted CRE is only accessible in both B cells and HSCs (Fig. 5k), and contained 5 out of the top 10 motifs, ZEB1/2, Slug, Ascl2, HEB, and E2A (Fig. 5l). In T cells, the gene Inducible T cell co-stimulator (Icos) is only expressed in T cells (Fig. 5m), a predicted linked CRE is accessible in both T cells and HSCs (Fig. 5n) and contains motifs for CTCF and WT1 (Fig. 5o). Taken together, these examples represent CRE priming in HSCs, along with the corresponding transcription factors that may act on each element to guide HSC fate.
MkPs and HSCs have the most similar accessibility profile
Here, we compared the genome-wide accessibility by ATAC-seq of the multipotent HSCs and unipotent lineage cell types (EPs, MkPs, GMs, B, and T cells). Through hierarchical clustering analysis, we observed erythromyeloid and lymphoid relationships that are consistent with the classical model of hematopoiesis (Fig. 1d) [4, 7, 28, 39]. By both PCA and hierarchical clustering, we observed that MkPs were the most similar to HSCs based on their accessibility profiles (Fig. 1). This relationship is reflected in a high level of overlap of peaks, as MkPs had the fewest peaks gained or lost from HSCs compared to the other cell types (Fig. 3) and had the largest percentage of peaks exclusively shared with HSCs (Fig. 4b). These findings are in agreement with recent clonal studies of hematopoiesis that reported a megakaryocyte lineage bias of HSCs [14, 37]. According to hierarchal clustering, EPs had the second closest association to HSCs (Fig. 1d) possibly supporting erythropoiesis as the default fate for hematopoiesis  under conditions where chromatin remodeling silences megakaryocyte driver elements . On the other end of the spectrum, the least similar cell types to HSCs were the lymphoid cell types (Fig. 1d). This greater difference was primarily due to a high proportion of peaks gained (Fig. 3e) rather than lost (Fig. 3f) upon differentiation from HSCs, leading to a greater ratio of peaks gained:lost for lymphoid cells than for erythromyeloid lineages (Fig. 3g).
Evidence of multilineage priming in HSCs
The priming of genes for transcription likely initiates within CREs, which can then drive the activation of promoter targets. These enhancers can act as drivers of lineage fate  and their accessibility is a putative regulator of competence in stem cells. We made the assumption that peaks that are exclusively shared between HSCs and the unipotent lineage cells contain CREs that are specific for driving differentiation into that lineage. We observed that the majority of exclusively shared peaks were non-promoter peaks (Fig. 4b) and were enriched for binding motifs of transcription factors known to be important for differentiation into each lineage (Fig. 4d–h). The enrichment of binding sites for known lineage-specific transcription factors suggests that many of the accessible sites may play functional roles. Additionally, about one-third of the exclusively shared ATAC peaks were enriched for the H3K4me1 histone modification, which is linked to a primed enhancer state , indicated that a subset are likely functional enhancers (Additional file 3: Figure S3); other ATAC-accessible elements may mark transcription start sites for non-coding genes, which are abundant and highly tissue-specific in the mouse genome . By using the GREAT tool, we made predictions for the target genes for the many ATAC-identified putative CREs that were present in the HSC/mature cell exclusive lists. The examples shown in Fig. 5 provide evidence that multilineage priming exists in HSCs.
Both permissive and de novo epigenetic mechanisms influence hematopoiesis
Analogous to other stem cell systems, multipotent HSCs with the competence to differentiate into diverse cell types reside at the top of the blood cell hierarchy. We tested two potential models of the mechanism of multipotency, the permissive fate and de novo activation (Fig. 1a). We found evidence for both. Supporting the permissive fate model are the observations that HSCs had the highest global accessibility (Fig. 3a/b), that peaks were lost in every unipotent cell type from HSCs (Fig. 3f), that every unipotent cell type shared some peaks exclusively with HSCs (Fig. 4b), and that evidence of multilineage priming of CREs were found in HSCs (Fig. 5). The de novo activation model was supported by the observation that new peaks were gained during differentiation into all five lineages (Fig. 3e), and previous studies reporting progressive upregulation of lineage-specific genes as HSCs transition into progenitors [18, 43]. Interestingly, in the β-globin locus, HS2, the strongest enhancer of globin expression [2, 16]], was highly accessible in HSCs, whereas the other HSs were not (Additional file 2: Figure S2). Thus, “priming” of this locus may occur in HSCs via HS2 (adhering to the permissive model of Fig. 1a), followed by induced accessibility (de novo model, Fig. 1a) of the other HSs and active β-globin expression upon erythroid differentiation. Thus, both permissive and de novo mechanisms likely influence hematopoietic fate decisions. Interestingly, we found evidence that the balance between the two models varies between lineages. For example, B cells, and to a lesser extent T cells, had a higher proportion of peaks gained than lost compared to erythromyeloid lineages (Fig. 3g). This may indicate that the megakaryocyte/erythroid lineage is in a more primed state in HSCs, whereas lymphopoiesis requires more extensive chromatin remodeling to both prime lymphoid CREs not accessible in HSCs and simultaneously shut down the megakaryocyte/erythrocyte trajectory. The cell output and kinetics from in vivo lineage tracing and reconstitution assays support these conclusions [4,5,6, 14, 37, 48]. Our identification of specific, putative regulatory CREs will enable functional testing of these elements.
Mice and cells
All experiments were performed using 8- to 12-week-old C57BL/6 wild-type mice in accordance with UCSC IACUC guidelines. Hematopoietic cells were isolated from BM by crushing murine femurs, tibias, hips, and sternums as previously described . Stem and progenitor cell fractions were enriched using CD117-coupled magnetic beads (Miltenyi). Cells were stained with unconjugated lineage rat antibodies (CD3, CD4, CD5, CD8, B220, Gr1, Mac1, and Ter119) followed by goat-α-rat PE-Cy5 (Invitrogen). Stem and progenitor cells were isolated using fluorescently labeled or biotinylated antibodies for the following antigens: cKit (2B8, Biolegend), Sca1 (D7, Biolegend), Slamf1(CD150) (TC15-12F12.2, Biolegend), CD41(MWReg30, Biolegend), and CD71(RI7217, Biolegend). Cells were sorted using a FACS Aria II (BD Bioscience). HSCs were defined as cKit+ Lin− Sca1+ Flk2− and Slamf1+; MkPs as cKit+Lin−Sca1−Slamf1−CD41+. Unipotent lineage cells were isolated by the following markers and as described previously [15, 29]: EPs, Lin(CD3, CD4, CD5, CD8, B220, Gr1, and Mac1)− CD71+Ter119±; GMs, Lin(CD3, CD4, CD5, CD8, B220, and Ter119)− Gr1+Mac1+ (“GM” cells were positive for both Gr1 and Mac1); T cells, Lin(CD5, B220, Gr1, Mac1, and Ter119)− CD25−CD3+CD4±CD8±; B cells, Lin(CD3, CD4, CD8, Gr1, Mac1, and Ter119)−CD43−B220+.
ATAC-seq was performed as previously described . Briefly, cells were collected after sorting into microcentrifuge tubes containing staining media (1xDPBS,1 mM EDTA with 5% serum). They were centrifuged at 500×g for 5 min at 4 ˚C to pellet the cells. The supernatant was aspirated, and the cells were washed with ice-cold 1xDPBS. Cells were centrifuged and the supernatant was discarded. Cells were immediately resuspended in ice-cold lysis buffer (10 mM Tris–HCl, pH 7.4, 10 mM NaCl, 3 mM MgCl2 and 0.1% IGEPAL CA-630) and centrifuged at 500×g for 10 min. The supernatant was aspirated, and pellets were resuspended in transposase reaction mix (25 µL 2xTD Buffer, 2.5 µL transposase (Illumina), and 22.5 µL nuclease-free water). The transposition reaction was carried out at 37 ˚C for 30 min at 600 rpm in a shaking thermomixer (Eppendorf). Immediately after completion of the transposition reaction, the samples were purified using the MinElute Reaction Clean up kit (Qiagen) and eluted into 10 µL of EB. Samples were stored at – 20 ˚C until PCR amplification step. PCR amplification was performed as previously described  using custom Nextera primers. After initial amplification, a portion of the samples were run on qPCR (ViiA7 Applied Biosystems) to determine the additional number of cycles needed for each library. The libraries were purified using the MinElute Reaction Clean up kit (Qiagen), eluted into 20 µL EB and then size selected using AmpureXP(Beckman-Coulter) beads at a ratio of 1.8:1 beads/sample, and eluted into 40 µL of nuclease-free water. Library size distribution was determined by Bioanalyzer (Agilent) capillary electrophoresis and library concentration was determined by Qubit 3 (Life Technologies). Quality of libraries was checked by shallow sequencing (1 million raw reads) on a Miseq (Illumina) at 75 × 75 paired-end sequencing. Those libraries that appeared to have size distributions similar to previous reports (Additional file 1: Figure S1) were pooled together and deep sequenced on a HiSeq2500 (Illumina) at 100 × 100 reads at the Vincent J. Coates Genomics Sequencing Laboratory at UC Berkeley.
Demultiplexed sequencing data were processed using the ENCODE ATAC-seq pipeline version 1.1.6 and 1.4.2 (https://github.com/ENCODE-DCC/atac-seq-pipeline) using the mm10 assembly and the default parameters. In version 1.4.2 changed: atac.multimapping = 0, atac.smooth_win = 150, atac.enable_idr = true, atac.idr_thresh = 0.1 to be consistent with the mapping/peak calling performed with previous versions.
Peak filtering, hierarchical clustering, and tSNE plot production were performed using the chromVAR package (https://github.com/GreenleafLab/chromVAR). First, the optimal peak-list from the IDR output for each cell type was concatenated and sorted, then used as the peak input for chromVAR. The blacklist filtered bam files for reach replicate was used as input along with the sorted peak file. The fragment counts in each peak for each replicate and GC bias was calculated, and then the peaks were filtered using filterPeaks function with the default parameters and non-overlapping = TRUE. The master peak-list was extracted at this point, which contained 84,243 peaks, and used throughout the study. The deviations were calculated using every peak, and the tSNE and correlation functions were also performed using the deviations output and the default parameters.
Annotation of peaks, generation of histogram plot, merging of peaks, and motif enrichment were performed by HOMER (http://homer.ucsd.edu/homer/). Peaks were annotated using the annotatePeaks.pl function with the mm10 assembly and default parameters. Histogram was created by first shifting the bam files using DeepTools alignmentSieve.py with the flag –ATACshift. Next, tag directories were made using the Tn5 shifted bam files using HOMER makeTagDirectory. The histogram was made using the annotatePeaks.pl function with the default settings and the flags: -size -500,500 and -hist 5. Peak lists were compared using the mergePeaks.pl function with default settings and the flags -d given, -venn, and for the unique peak-lists -prefix. Motif enrichment was performed using the findMotifsGenome.pl package with default parameters using the flag -size given and custom background peaks, which consisted of the combination of all the peak-lists for the cell types not being analyzed. Instances of motifs in non-promoter peaks were found by using the annotatePeaks.pl function with the -m flag, using custom made motif files for each cell type containing the top 10 enriched motifs found.
The GREAT tool (http://great.stanford.edu/public/html/) was used to annotate non-promoter peaks to target genes. The peak-lists were reduced to BED4 files from the HOMER annotations output and used as input. The whole mm10 genome was used as the background regions, and the association rule settings were set as Basal plus extension, proximal window 2 kb upstream, 1 kb downstream, plus distal up to 1 Mb and included curated regulatory domains. All genome track visualizations were made using the UCSC genome browser. Graphs were made in either Microsoft Excel or GraphPad Prism 8. Annotations to figures were performed using Adobe Illustrator CC and Adobe Photoshop CC.
ChIP data were handled as follows: the enhancer list from  was mapped to mm10 using the liftOver tool, then compared to the master peak-list. The raw sequencing data for H3K4me1 and H3K27Ac in LT-HSCs were downloaded from GEO and mapping to mm10 and peak calling were performed using the parameters listed in the publication .
Availability of data and materials
The datasets generated during and/or analyzed during the current study are available in the Gene Expression Omnibus (GEO), accession number GSE162949.
An X, Schulz VP, Li J, Wu K, Liu J, Xue F, Hu J, Mohandas N, Gallagher PG. Global transcriptome analyses of human and murine terminal erythroid differentiation. Blood. 2014;123:3466–77.
Bender MA, Ragoczy T, Lee J, Byron R, Telling A, Dean A, Groudine M. The hypersensitive sites of the murine β-globin locus control region act independently to affect nuclear localization and transcriptional elongation. Blood. 2012;119:3820–7.
Bernstein BE, Mikkelsen TS, Xie X, Kamal M, Huebert DJ, Cuff J, Fry B, Meissner A, Wernig M, Plath K, et al. A bivalent chromatin structure marks key developmental genes in embryonic stem cells. Cell. 2006;125:315–26.
Boyer SW, Schroeder AV, Smith-Berdan S, Forsberg EC. All hematopoietic cells develop from hematopoietic stem cells through Flk2/Flt3-positive progenitor cells. Cell Stem Cell. 2011;9:64–73.
Boyer SW, Beaudin AE, Forsberg EC. Mapping differentiation pathways from hematopoietic stem cells using Flk2/Flt3 lineage tracing. Cell Cycle. 2012;11:3180–8.
Boyer SW, Rajendiran S, Beaudin AE, Smith-berdan S, Muthuswamy PK, Perez-Cunningham J, Martin EW, Cheung C, Tsang H, Landon M, et al. Clonal and quantitative in vivo assessment of hematopoietic stem cell differentiation reveals strong erythroid potential of multipotent cells. Stem Cell Rep. 2019;12:801–15.
Bryder D, Rossi DJ, Weissman IL. Hematopoietic stem cells: the paradigmatic tissue-specific stem cell. Am J Pathol. 2006;169:338–46.
Buenrostro JD, Giresi PG, Zaba LC, Chang HY, Greenleaf WJ. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat Meth. 2013;10:1213–8.
Buenrostro JD, Corces MR, Lareau CA, Wu B, Schep AN, Aryee MJ, Majeti R, Chang HY, Greenleaf WJ. Integrated single-cell analysis maps the continuous regulatory landscape of human hematopoietic differentiation. Cell. 2018;173:1–14.
Bulut-Karslioglu A, Macrae TA, Oses-Prieto JA, Covarrubias S, Percharde M, Ku G, Diaz A, McManus MT, Burlingame AL, Ramalho-Santos M. The transcriptionally permissive chromatin state of embryonic stem cells is acutely tuned to translational output. Cell Stem Cell. 2018;22:369–83.
Cabal-Hierro L, van Galen P, Prado MA, Higby KJ, Togami K, Mowery CT, Paulo JA, Xie Y, Cejas P, Furusawa T, et al. Chromatin accessibility promotes hematopoietic and leukemia stem cell activity. Nat Commun. 2020;11:1406.
Cabezas-Wallscheid N, Klimmeck D, Hansson J, Lipka DB, Reyes A, Wang Q, Weichenhan D, Lier A, Von Paleske L, Renders S, et al. Identification of regulatory networks in HSCs and their immediate progeny via integrated proteome, transcriptome, and DNA methylome analysis. Cell Stem Cell. 2014;15:507–22.
Calo E, Wysocka J. Modification of enhancer chromatin: what, how, and why? Mol Cell. 2013;49:825–37.
Carrelha J, Meng Y, Kettyle LM, Luis TC, Norfo R, Alcolea V, Boukarabila H, Grasso F, Gambardella A, Grover A, et al. Hierarchically related lineage-restricted fates of multipotent haematopoietic stem cells. Nature. 2018;554:106–11.
Cool T, Worthington A, Poscablo D, Hussaini A, Forsberg EC. Interleukin 7 receptor is required for myeloid cell homeostasis and reconstitution by hematopoietic stem cells. Exp Hematol. 2020;90:39-45.e3.
Fiering S, Epner E, Robinson K, Zhuang Y, Telling A, Hu M, Martin DI, Enver T, Ley TJ, Groudine M. Targeted deletion of 5’HS2 of the murine beta-globin LCR reveals that it is not essential for proper regulation of the beta-globin locus. Genes Dev. 1995;15:2203–13.
Forsberg EC, Downs KM, Christensen HM, Im H, Nuzzi PA, Bresnick EH. Developmentally dynamic histone acetylation pattern of a tissue-specific chromatin domain. Proc Natl Acad Sci. 2000;97:14494–9.
Forsberg EC, Prohaska SS, Katzman S, Heffner GC, Stuart JM, Weissman IL. Differential expression of novel potential regulators in hematopoietic stem cells. PLoS Genet. 2005;1(3):e28.
Gaspar-Maia A, Alajem A, Polesso F, Sridharan R, Mason MJ, Heidersbach A, Ramalho-Santos J, McManus MT, Plath K, Meshorer E, et al. Chd1 regulates open chromatin and pluripotency of embryonic stem cells. Nature. 2009;460:863–8.
Gaspar-Maia A, Alajem A, Meshorer E, Ramalho-Santos M. Open chromatin in pluripotency and reprogramming. Nat Rev Mol Cell Biol. 2011;12:36–47.
Growney JD, Shigematsu H, Li Z, Lee BH, Adelsperger J, Rowan R, Curley DP, Kutok JL, Akashi K, Williams IR, et al. Loss of Runx1 perturbs adult hematopoiesis and is associated with a myeloproliferative phenotype. Blood. 2005;106:494–504.
Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, Cheng JX, Murre C, Singh H, Glass CK. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. 2010;38:576–89.
Hestdal K, Ruscetti FW, Ihle JN, Jacobsen SE, Dubois CM, Kopp WC, Longo DL, Keller JR. Characterization and regulation of RB6–8C5 antigen expression on murine bone marrow cells. J Immunol. 1991;147:22–8.
Heuston EF, Keller CA, Lichtenberg J, Giardine B, Anderson SM, Hardison RC, Bodine DM. Establishment of regulatory elements during erythro-megakaryopoiesis identifies hematopoietic lineage-commitment points. Epigenet Chromat. 2018;11:1–18.
Hu M, Krause D, Greaves M, Sharkis S, Dexter M, Heyworth C, Enver T. Multilineage gene expression precedes commitment in the hemopoietic system. Genes Dev. 1997;11:774–85.
Kruse EA, Loughran SJ, Baldwin TM, Josefsson EC, Ellis S, Watson DK, Nurden P, Metcalf D, Hilton DJ, Alexander WS, et al. Dual requirement for the ETS transcription factors Fli-1 and Erg in hematopoietic stem cells and the megakaryocyte lineage. Proc Natl Acad Sci. 2009;106:13814–9.
Lara-Astiaso D, Weiner A, Lorenzo-Vivas E, Zaretsky I, Jaitin DA, David E, Keren-Shaul H, Mildner A, Winter D, Jung S, et al. Chromatin state dynamics during blood formation. Science. 2014;55:1–10.
Laurenti E, Göttgens B. From haematopoietic stem cells to complex differentiation landscapes. Nature. 2018;553:418–26.
Leung GA, Cool T, Valencia CH, Worthington A, Beaudin AE, Camilla Forsberg E. The lymphoid-associated interleukin 7 receptor (IL7R) regulates tissue-resident macrophage development. Development. 2019;146:dev176180.
Li Q, Peterson KR, Fang X, Stamatoyannopoulos G. Locus control regions. Blood. 2002;100:3077–86.
Li Q, Brown JB, Huang H, Bickel PJ. Measuring reproducibility of high-throughput experiments. Ann Appl Stat. 2011;5:1752–79.
Månsson R, Hultquist A, Luc S, Yang L, Anderson K, Kharazi S, Al-Hashmi S, Liuba K, Thorén L, Adolfsson J, et al. Molecular evidence for hierarchical transcriptional lineage priming in fetal and adult stem cells and multipotent progenitors. Immunity. 2007;26:407–19.
McLean CY, Bristor D, Hiller M, Clarke SL, Schaar BT, Lowe CB, Wenger AM, Bejerano G. GREAT improves functional interpretation of cis-regulatory regions. Nat Biotechnol. 2010;28:495–501.
Palstra RJ, de Laat W, Grosveld F. Beta-globin regulation and long-range interactions. 2008;61:107–42.
Rajendiran S, Smith-Berdan S, Kunz L, Risolino M, Selleri L, Schroeder T, Forsberg EC. Ubiquitous overexpression of CXCL12 confers radiation protection and enhances mobilization of hematopoietic stem and progenitor cells. Stem Cells. 2020;38:1159–74.
Ravasi T, Suzuki H, Pang KC, Katayama S, Furuno M, Okunishi R, Fukuda S, Ru K, Frith MC, Gongora MM, Grimmond SM, Hume DA, Hayashizaki Y, Mattick JS. Experimental validation of the regulated expression of large numbers of non-coding RNAs from the mouse genome. Genome Res. 2006;16:11–9.
Rodriguez-Fraticelli AE, Wolock SL, Weinreb CS, Panero R, Patel SH, Jankovic M, Sun J, Calogero RA, Klein AM, Camargo FD. Clonal analysis of lineage fate in native haematopoiesis. Nature. 2018;553:212–6.
Schep AN, Wu B, Buenrostro JD, Greenleaf WJ. ChromVAR: Inferring transcription-factor-associated accessibility from single-cell epigenomic data. Nat Methods. 2017;14:975–8.
Seita J, Weissman IL. Hematopoietic stem cell: self-renewal versus differentiation. WIREs Syst Biol Med. 2010;2:640–53.
Seita J, Sahoo D, Rossi DJ, Bhattacharya D, Serwold T, Inlay MA, Ehrlich LIR, Fathman JW, Dill DL, Weissman IL. Gene expression commons: An open platform for absolute gene expression profiling. PLoS ONE. 2012;7:1–11.
Shivdasani RA, Fujiwara Y, McDevitt MA, Orkin SH. A lineage-selective knockout establishes the critical role of transcription factor GATA-1 in megakaryocyte growth and platelet development. EMBO J. 1997;16:3965–73.
Siegwart LC, Schwemmers S, Wehrle J, Koellerer C, Seeger T, Gründer A, Pahl HL. The transcription factor NFE2 enhances expression of the hematopoietic master regulators SCL/TAL1 and GATA2. Exp Hematol. 2020;87:42-47.e1.
Terskikh AV, Miyamoto T, Chang C, Diatchenko L, Weissman IL. Gene expression analysis of purified hematopoietic stem cells and committed progenitors. Blood. 2003;102:94–101.
Ugarte F, Sousae R, Cinquin B, Martin EW, Krietsch J, Sanchez G, Inman M, Tsang H, Warr M, Passegué E, et al. Progressive chromatin condensation and H3K9 methylation regulate the differentiation of embryonic and hematopoietic stem cells. Stem Cell Rep. 2015;5:728–40.
Waddington CH. Organisers and genes. Cambridge: University Press; 1940.
Wang A, Yue F, Li Y, Xie R, Harper T, Patel NA, Muth K, Palmer J, Qiu Y, Wang J, et al. Epigenetic priming of enhancers predicts developmental competence of hESC-derived endodermal lineage intermediates. Cell Stem Cell. 2015;16:386–99.
Xu J, Watts JA, Pope SD, Gadue P, Kamps M, Plath K, Zaret KS, Smale ST. Transcriptional competence and the active marking of tissue-specific enhancers by defined transcription factors in embryonic and induced pluripotent stem cells. Genes Dev. 2009;23:2824–38.
Yamamoto R, Morita Y, Ooehara J, Hamanaka S, Onodera M, Rudolph KL, Ema H, Nakauchi H. Clonal analysis unveils self-renewing lineage-restricted progenitors generated directly from hematopoietic stem cells. Cell. 2013;154:1112–26.
We thank Bari Nazario and the IBSC flow cytometry core for assistance and support; Sol Katzman for bioinformatic assistance; and Forsberg lab members for comments on the manuscript.
This work was supported by NIH awards (R01HL115158, DK100917, R01AG062879) to E.C.F.; by NIH/NHLBI fellowship (F31HL144115) to E.W.M.; by CIRM SCILL grant TB1-01195 to E.W.M. via San Jose State University; by CIRM Training grant TG2-01157 to J.K.; by a UCSC Genomic Sciences Graduate Training Program from NIH/NHGRI (NIH T32 HG008345) to R.E.R., by a UCSC IMSD award from NIH/NIGMS to R.S. (R25GM058903); by the Baskin School of Engineering and the Ken and Glory Levy Fund for RNA Biology to D.H.K., and by CIRM Shared Stem Cell Facilities (CL1-00506) and CIRM Major Facilities (FA1-00617-1) awards to the University of California, Santa Cruz.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Library fragment distributions for ATAC-seq samples. The library size distribution after deep-sequencing, mapping, and filtering to unique reads is shown of both replicates for A) HSCs B) EPs C) MkPs D) GMs E) B cells, and F) T cells.
Erythroid-selective accessibility of the β-globin cluster. ATAC-seq signal tracks of the six cell types in this study at the β-globin cluster (chr7: 103,792,027–103,879,340; mm10). The adult globin genes β-major (ßmaj) and β-minor (ßmin), as well as the hypersensitive sites (HS1-4,6) of the Locus Control Region (LCR) that regulates expression of the genes in this locus, displayed accessibility in EPs. HS2, but not the other HSs, and the β-major promoter were also accessible in HSCs, possibly indicating a “permissive” chromatin state. Accessibility of HS2 and HS4 in MkPs may relate to a closer relationship to HSCs and/or EPs (Fig. 1). As expected, no accessibility was observed at the fetal-specific epsilon Y globin (Ey), β-h1 (ßh1), β-h2 (ßh2) genes, or HS5. Likewise, GMs, B and T cells, that do not express β-globin genes, did not display accessibility of any of the regulatory elements in the locus.
Lineage specific, HSC-primed peaks were marked by H3K4me1 and not H3K27Ac. A) About one in three (21,085 peaks out of 71,072) of our ATAC-seq non-promoter peaks in the master peak-list overlapped with peaks designated as probable enhancers based on H3K4me1 and H3K27Ac ChIP data (Lara-Astiaso et al.). B) About one in three (32.5%) of all HSC-primed peaks for the five unipotent lineage cell types were also marked by the histone modification H3K4me1, and 2.8% were marked by H3K27Ac. C) HSC-primed peaks for each unipotent lineage were primarily marked by H3K4me1 and not H3K27Ac. Results in panel B represent the aggregate of the results shown in panel C.
About this article
Cite this article
Martin, E.W., Krietsch, J., Reggiardo, R.E. et al. Chromatin accessibility maps provide evidence of multilineage gene priming in hematopoietic stem cells. Epigenetics & Chromatin 14, 2 (2021). https://doi.org/10.1186/s13072-020-00377-1