Polycomb contraction differentially regulates terminal human hematopoietic differentiation programs

Background

Lifelong production of the many types of mature blood cells from less differentiated progenitors is a hierarchically ordered process that spans multiple cell divisions. The nature and timing of the molecular events required to integrate the environmental signals, transcription factor activity, epigenetic modifications, and changes in gene expression involved are thus complex and still poorly understood. To address this gap, we generated comprehensive reference epigenomes of 8 phenotypically defined subsets of normal human cord blood.

Results

We describe a striking contraction of H3K27me3 density in differentiated myelo-erythroid cells that resembles a punctate pattern previously ascribed to pluripotent embryonic stem cells. Phenotypically distinct progenitor cell types display a nearly identical repressive H3K27me3 signature characterized by large organized chromatin K27-modification domains that are retained by mature lymphoid cells but lost in terminally differentiated monocytes and erythroblasts. We demonstrate that inhibition of polycomb group members predicted to control large organized chromatin K27-modification domains influences lymphoid and myeloid fate decisions of primary neonatal hematopoietic progenitors in vitro. We further show that a majority of active enhancers appear in early progenitors, a subset of which are DNA hypermethylated and become hypomethylated and induced during terminal differentiation.

Conclusion

Primitive human hematopoietic cells display a unique repressive H3K27me3 signature that is retained by mature lymphoid cells but is lost in monocytes and erythroblasts. Intervention data implicate that control of this chromatin state change is a requisite part of the process whereby normal human hematopoietic progenitor cells make lymphoid and myeloid fate decisions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01315-1.

Epigenetic Restoration of Fetal-like IGF1 Signaling Inhibits Leukemia Stem Cell Activity

Acute leukemias are aggressive malignancies of developmentally arrested hematopoietic progenitors. We sought here to explore the possibility that changes in hematopoietic stem/progenitor cells during development might alter the biology of leukemias arising from this tissue compartment. Using a mouse model of acute T cell leukemia, we found that leukemias generated from fetal liver (FL) and adult bone marrow (BM) differed dramatically in their leukemia stem cell activity with FL leukemias showing markedly reduced serial transplantability as compared to BM leukemias. We present evidence that this difference is due to NOTCH1-driven autocrine IGF1 signaling, which is active in FL cells but restrained in BM cells by EZH2-dependent H3K27 trimethylation. Further, we confirmed this mechanism is operative in human disease and show that enforced IGF1 signaling effectively limits leukemia stem cell activity. These findings demonstrate that resurrecting dormant fetal programs in adult cells may represent an alternate therapeutic approach in human cancer.

Nucleosome Density ChIP-Seq Identifies Distinct Chromatin Modification Signatures Associated with MNase Accessibility

Nucleosome position, density, and post-translational modification are widely accepted components of mechanisms regulating DNA transcription but still incompletely understood. We present a modified native ChIP-seq method combined with an analytical framework that allows MNase accessibility to be integrated with histone modification profiles. Application of this methodology to the primitive (CD34+) subset of normal human cord blood cells enabled genomic regions enriched in one versus two nucleosomes marked by histone 3 lysine 4 trimethylation (H3K4me3) and/or histone 3 lysine 27 trimethylation (H3K27me3) to be associated with their transcriptional and DNA methylation states. From this analysis, we defined four classes of promoter-specific profiles and demonstrated that a majority of bivalent marked promoters are heterogeneously marked at a single-cell level in this primitive cell type. Interestingly, extension of this approach to human embryonic stem cells revealed an altered relationship between chromatin modification state and nucleosome content at promoters, suggesting developmental stage-specific organization of histone methylation states.

Analysis of Normal Human Mammary Epigenomes Reveals Cell-Specific Active Enhancer States and Associated Transcription Factor Networks

The normal adult human mammary gland is a continuous bilayered epithelial system. Bipotent and myoepithelial progenitors are prominent and unique components of the outer (basal) layer. The inner (luminal) layer includes both luminal-restricted progenitors and a phenotypically separable fraction that lacks progenitor activity. We now report an epigenomic comparison of these three subsets with one another, with their associated stromal cells, and with three immortalized, non-tumorigenic human mammary cell lines. Each genome-wide analysis contains profiles for six histone marks, methylated DNA, and RNA transcripts. Analysis of these datasets shows that each cell type has unique features, primarily within genomic regulatory regions, and that the cell lines group together. Analyses of the promoter and enhancer profiles place the luminal progenitors in between the basal cells and the non-progenitor luminal subset. Integrative analysis reveals networks of subset-specific transcription factors.

Characterization of the human thyroid epigenome

The thyroid gland, necessary for normal human growth and development, functions as an essential regulator of metabolism by the production and secretion of appropriate levels of thyroid hormone. However, assessment of abnormal thyroid function may be challenging suggesting a more fundamental understanding of normal function is needed. One way to characterize normal gland function is to study the epigenome and resulting transcriptome within its constituent cells. This study generates the first published reference epigenomes for human thyroid from four individuals using ChIP-seq and RNA-seq. We profiled six histone modifications (H3K4me1, H3K4me3, H3K27ac, H3K36me3, H3K9me3, H3K27me3), identified chromatin states using a hidden Markov model, produced a novel quantitative metric for model selection and established epigenomic maps of 19 chromatin states. We found that epigenetic features characterizing promoters and transcription elongation tend to be more consistent than regions characterizing enhancers or Polycomb-repressed regions and that epigenetically active genes consistent across all epigenomes tend to have higher expression than those not marked as epigenetically active in all epigenomes. We also identified a set of 18 genes epigenetically active and consistently expressed in the thyroid that are likely highly relevant to thyroid function. Altogether, these epigenomes represent a powerful resource to develop a deeper understanding of the underlying molecular biology of thyroid function and provide contextual information of thyroid and human epigenomic data for comparison and integration into future studies.

Genome-Wide Profiles of Extra-cranial Malignant Rhabdoid Tumors Reveal Heterogeneity and Dysregulated Developmental Pathways

Malignant rhabdoid tumors (MRTs) are rare lethal tumors of childhood that most commonly occur in the kidney and brain. MRTs are driven by SMARCB1 loss, but the molecular consequences of SMARCB1 loss in extra-cranial tumors have not been comprehensively described and genomic resources for analyses of extra-cranial MRT are limited. To provide such data, we used whole-genome sequencing, whole-genome bisulfite sequencing, whole transcriptome (RNA-seq) and microRNA sequencing (miRNA-seq), and histone modification profiling to characterize extra-cranial MRTs. Our analyses revealed gene expression and methylation subgroups and focused on dysregulated pathways, including those involved in neural crest development.

Integrative analysis of 111 reference human epigenomes

The reference human genome sequence set the stage for studies of genetic variation and its association with human disease, but epigenomic studies lack a similar reference. To address this need, the NIH Roadmap Epigenomics Consortium generated the largest collection so far of human epigenomes for primary cells and tissues. Here we describe the integrative analysis of 111 reference human epigenomes generated as part of the programme, profiled for histone modification patterns, DNA accessibility, DNA methylation and RNA expression. We establish global maps of regulatory elements, define regulatory modules of coordinated activity, and their likely activators and repressors. We show that disease- and trait-associated genetic variants are enriched in tissue-specific epigenomic marks, revealing biologically relevant cell types for diverse human traits, and providing a resource for interpreting the molecular basis of human disease. Our results demonstrate the central role of epigenomic information for understanding gene regulation, cellular differentiation and human disease.

DNA barcoding reveals diverse growth kinetics of human breast tumour subclones in serially passaged xenografts

Genomic and phenotypic analyses indicate extensive intra- as well as intertumoral heterogeneity in primary human malignant cell populations despite their clonal origin. Cellular DNA barcoding offers a powerful and unbiased alternative to track the number and size of multiple subclones within a single human tumour xenograft and their response to continued in vivo passaging. Using this approach we find clone-initiating cell frequencies that vary from ~1/10 to ~1/10,000 cells transplanted for two human breast cancer cell lines and breast cancer xenografts derived from three different patients. For the cell lines, these frequencies are negatively affected in transplants of more than 20,000 cells. Serial transplants reveal five clonal growth patterns (unchanging, expanding, diminishing, fluctuating or of delayed onset), whose predominance is highly variable both between and within original samples. This study thus demonstrates the high growth potential and diverse growth properties of xenografted human breast cancer cells.

Cancer cells within the same tumour are heterogeneous in their tumorigenic potential, differentiation status and sensitivity to treatments. Here Nguyen et al. use a sensitive DNA barcoding method to characterize the diversity of clonal growth behaviour within human breast tumours.

HISTONE DEACETYLASE 9 represses seedling traits in Arabidopsis thaliana dry seeds

Plant life is characterized by major phase changes. We studied the role of histone deacetylase (HDAC) activity in the transition from seed to seedling in Arabidopsis. Pharmacological inhibition of HDAC stimulated germination of freshly harvested seeds. Subsequent analysis revealed that histone deacetylase 9 (hda9) mutant alleles displayed reduced seed dormancy and faster germination than wild-type plants. Transcriptome meta-analysis comparisons between the hda9 dry seed transcriptome and published datasets demonstrated that transcripts of genes that are induced during imbibition in wild-type prematurely accumulated in hda9-1 dry seeds. This included several genes associated with photosynthesis and photoautotrophic growth such as RuBisCO and RuBisCO activase (RCA). Chromatin immunoprecipitation experiments demonstrated enhanced histone acetylation levels at their loci in young hda9-1 seedlings. Our observations suggest that HDA9 negatively influences germination and is involved in the suppression of seedling traits in dry seeds, probably by transcriptional repression via histone deacetylation. Accordingly, HDA9 transcript is abundant in dry seeds and becomes reduced during imbibition in wild-type seeds. The proposed function of HDA9 is opposite to that of its homologous genes HDA6 and HDA19, which have been reported to repress embryonic properties in germinated seedlings.

E-cadherin is required for the proper activation of the Lifr/Gp130 signaling pathway in mouse embryonic stem cells

The leukemia inhibitory factor (Lif) signaling pathway is a crucial determinant for mouse embryonic stem (mES) cell self-renewal and pluripotency. One of the hallmarks of mES cells, their compact growth morphology, results from tight cell adhesion mediated through E-cadherin, β-catenin (Ctnnb1) and α-catenin with the actin cytoskeleton. β-catenin is also involved in canonical Wnt signaling, which has also been suggested to control mES cell stemness. Here, we analyze Ctnnb1(-/-) mES cells in which cell adhesion is preserved by an E-cadherin-α-catenin (Eα) fusion protein (Ctnnb1(-/-)Eα mES cells), and show that mimicking only the adhesive function of β-catenin is necessary and sufficient to maintain the mES cell state, making β-catenin/Wnt signaling obsolete in this process. Furthermore, we propose a role for E-cadherin in promoting the Lif signaling cascade, showing an association of E-cadherin with the Lifr-Gp130 receptor complex, which is most likely facilitated by the extracellular domain of E-cadherin. Without Eα, and thus without maintained cell adhesion, Ctnnb1(-/-) mES cells downregulate components of the Lif signaling pathway, such as Lifr, Gp130 and activated Stat3, as well as pluripotency-associated markers. From these observations, we hypothesize that the changes in gene expression accompanying the loss of pluripotency are a direct consequence of dysfunctional cell adhesion. Supporting this view, we find that the requirement for intact adhesion can be circumvented by the forced expression of constitutively active Stat3. In summary, we put forward a model in which mES cells can be propagated in culture in the absence of Ctnnb1, as long as E-cadherin-mediated cell adhesion is preserved.

Arabidopsis paired amphipathic helix proteins SNL1 and SNL2 redundantly regulate primary seed dormancy via abscisic acid-ethylene antagonism mediated by histone deacetylation

Histone (de)acetylation is a highly conserved chromatin modification that is vital for development and growth. In this study, we identified a role in seed dormancy for two members of the histone deacetylation complex in Arabidopsis thaliana, SIN3-LIKE1 (SNL1) and SNL2. The double mutant snl1 snl2 shows reduced dormancy and hypersensitivity to the histone deacetylase inhibitors trichostatin A and diallyl disulfide compared with the wild type. SNL1 interacts with HISTONE DEACETYLASE19 in vitro and in planta, and loss-of-function mutants of SNL1 and SNL2 show increased acetylation levels of histone 3 lysine 9/18 (H3K9/18) and H3K14. Moreover, SNL1 and SNL2 regulate key genes involved in the ethylene and abscisic acid (ABA) pathways by decreasing their histone acetylation levels. Taken together, we showed that SNL1 and SNL2 regulate seed dormancy by mediating the ABA-ethylene antagonism in Arabidopsis. SNL1 and SNL2 could represent a cross-link point of the ABA and ethylene pathways in the regulation of seed dormancy.

Tubulin tyrosine ligase like 12, a TTLL family member with SET- and TTL-like domains and roles in histone and tubulin modifications and mitosis

hTTLL12 is a member of the tubulin tyrosine ligase (TTL) family that is highly conserved in phylogeny. It has both SET-like and TTL-like domains, suggesting that it could have histone methylation and tubulin tyrosine ligase activities. Altered expression of hTTLL12 in human cells leads to specific changes in H4K20 trimethylation, and tubulin detyrosination, hTTLL12 does not catalyse histone methylation or tubulin tyrosination in vitro, as might be expected from the lack of critical amino acids in its SET-like and TTLL-like domains. hTTLL12 misexpression increases mitotic duration and chromosome numbers. These results suggest that hTTLL12 has non-catalytic functions related to tubulin and histone modification, which could be linked to its effects on mitosis and chromosome number stability.

Control and consequences of chromatin compaction during seed maturation in Arabidopsis thaliana

Nuclei in dry Arabidopsis thaliana seeds are very small and have highly condensed chromatin. Nuclear shrinkage and chromatin compaction occur during seed maturation and have been shown to be independent, developmentally controlled processes. To confirm this genetically, we studied chromatin compaction in a mutant of the seed developmental regulator ABA INSENSITIVE 3, and in a double mutant of the nuclear matrix proteins LITTLE NUCLEI 1 and 2. Our results indicated that the nuclear shrinking and chromatin condensation during seed maturation can be genetically uncoupled, confirming that these are independent processes. In addition, we demonstrated that transcript levels of siliques toward the end of seed maturation are comparable to those in vegetative tissues, despite the highly compacted chromatin, small nuclear volume and low hydration status of seeds.

Head and neck squamous cell carcinoma transcriptome analysis by comprehensive validated differential display

Head and neck squamous cell carcinoma (HNSCC) is common worldwide and is associated with a poor rate of survival. Identification of new markers and therapeutic targets, and understanding the complex transformation process, will require a comprehensive description of genome expression, that can only be achieved by combining different methodologies. We report here the HNSCC transcriptome that was determined by exhaustive differential display (DD) analysis coupled with validation by different methods on the same patient samples. The resulting 820 nonredundant sequences were analysed by high throughput bioinformatics analysis. Human proteins were identified for 73% (596) of the DD sequences. A large proportion (>50%) of the remaining unassigned sequences match ESTs (expressed sequence tags) from human tumours. For the functionally annotated proteins, there is significant enrichment for relevant biological processes, including cell motility, protein biosynthesis, stress and immune responses, cell death, cell cycle, cell proliferation and/or maintenance and transport. Three of the novel proteins (TMEM16A, PHLDB2 and ARHGAP21) were analysed further to show that they have the potential to be developed as therapeutic targets.

PipeAlign: A new toolkit for protein family analysis

PipeAlign is a protein family analysis tool integrating a five step process ranging from the search for sequence homologues in protein and 3D structure databases to the definition of the hierarchical relationships within and between subfamilies. The complete, automatic pipeline takes a single sequence or a set of sequences as input and constructs a high-quality, validated MACS (multiple alignment of complete sequences) in which sequences are clustered into potential functional subgroups. For the more experienced user, the PipeAlign server also provides numerous options to run only a part of the analysis, with the possibility to modify the default parameters of each software module. For example, the user can choose to enter an existing multiple sequence alignment for refinement, validation and subsequent clustering of the sequences. The aim is to provide an interactive workbench for the validation, integration and presentation of a protein family, not only at the sequence level, but also at the structural and functional levels. PipeAlign is available at http://igbmc.u-strasbg.fr/PipeAlign/.