A bivalent chromatin structure marks key developmental genes in embryonic stem cells

Exploring the roles of noncoding RNAs in craniofacial abnormalities: A systematic review

Congenital craniofacial abnormalities are congenital anomalies of variable expressivity and severity with a recognizable set of abnormalities, which are derived from five identifiable primordial structures. They can occur unilaterally or bilaterally and include various malformations such as cleft lip with/without palate, craniosynostosis, and craniofacial microsomia. To date, the molecular etiology of craniofacial abnormalities is largely unknown. Noncoding RNAs (ncRNAs), including microRNAs, long ncRNAs, circular RNAs and PIWI-interacting RNAs, function as major regulators of cellular epigenetic hallmarks via regulation of various molecular and cellular processes. Recently, aberrant expression of ncRNAs has been implicated in many diseases, including craniofacial abnormalities. Consequently, this review focuses on the role and mechanism of ncRNAs in regulating craniofacial development in the hope of providing clues to identify potential therapeutic targets.

Sequential ChIP-Seq

ChIP-Seq has been used extensively to profile genome-wide transcription factor binding and post-translational histone modifications. A sequential ChIP assay determines the in vivo co-localization of two proteins to the same genomic locus. In this chapter, we combine the two protocols in Sequential ChIP-Seq, a method for identifying genome-wide sites of in vivo protein co-occupancy.

Selective gene expression maintains human tRNA anticodon pools during differentiation

Transfer RNAs are essential for translating genetic information into proteins. The human genome contains hundreds of predicted tRNA genes, many in multiple copies. How their expression is regulated to control tRNA repertoires is unknown. Here we combined quantitative tRNA profiling and chromatin immunoprecipitation with sequencing to measure tRNA expression following the differentiation of human induced pluripotent stem cells into neuronal and cardiac cells. We find that tRNA transcript levels vary substantially, whereas tRNA anticodon pools, which govern decoding rates, are more stable among cell types. Mechanistically, RNA polymerase III transcribes a wide range of tRNA genes in human induced pluripotent stem cells but on differentiation becomes constrained to a subset we define as housekeeping tRNAs. This shift is mediated by decreased mTORC1 signalling, which activates the RNA polymerase III repressor MAF1. Our data explain how tRNA anticodon pools are buffered to maintain decoding speed across cell types and reveal that mTORC1 drives selective tRNA expression during differentiation.

Roles and Regulation of H3K4 Methylation During Mammalian Early Embryogenesis and Embryonic Stem Cell Differentiation

From generation of germ cells, fertilization, and throughout early mammalian embryonic development, the chromatin undergoes significant alterations to enable precise regulation of gene expression and genome use. Methylation of histone 3 lysine 4 (H3K4) correlates with active regions of the genome, and it has emerged as a dynamic mark throughout this timeline. The pattern and the level of H3K4 methylation are regulated by methyltransferases and demethylases. These enzymes, as well as their protein partners, play important roles in early embryonic development and show phenotypes in embryonic stem cell self-renewal and differentiation. The various roles of H3K4 methylation are interpreted by dedicated chromatin reader proteins, linking this modification to broader molecular and cellular phenotypes. In this review, we discuss the regulation of different levels of H3K4 methylation, their distinct accumulation pattern, and downstream molecular roles with an early embryogenesis perspective.

5-Hydroxymethylcytosine Profiling of Cell-Free DNA Identifies Bivalent Genes That Are Prognostic of Survival in High-Risk Neuroblastoma

PURPOSE

Neuroblastoma is the most common extracranial solid tumor in childhood. We previously showed that circulating cell-free DNA (cfDNA) and tumor biopsy derived 5-hydroxymethylcytosime (5-hmC) profiles identified patients with neuroblastoma who experienced subsequent relapse. Here, we hypothesized that 5-hmC modifications selectively enriched in cfDNA compared with tumor biopsy samples would identify epigenetic changes associated with aggressive tumor behavior and identify novel biomarkers of outcome in patients with high-risk neuroblastoma.

METHODS

5-hmC profiles from cfDNA (n = 64) and tumor biopsies (n = 48) were compared. Two neuroblastoma cell lines underwent chromatin immunoprecipitation followed by sequencing (ChIP-Seq) for H3K27me3, H3K4me3, and H3K27ac; kethoxal-associated single-stranded DNA sequencing; hmC-Seal for 5-hmC; and RNA-sequencing (RNA-Seq). Genes enriched for both H3K27me3 and H3K4me3 in the included cell lines were defined as bivalent. Using bivalent genes defined in vitro, a bivalent signature was established in three publicly available cohorts of patients with neuroblastoma through gene set variation analysis. Differences between tumors with high or low bivalent signatures were assessed by the Kaplan-Meier method and Cox proportional hazards models.

RESULTS

In cfDNA compared with tumor biopsy derived 5-hmC profiles, we found increased 5-hmC deposition on Polycomb Repressive Complex 2 target genes, a finding previously described in the context of bivalent genes. We identified 313 genes that bore bivalent chromatin marks, were enriched for mediators of neuronal differentiation, and were transcriptionally repressed across a panel of heterogeneous neuroblastoma cell lines. In three distinct clinical cohorts, low bivalent signature was significantly and independently associated with worse clinical outcome in patients with high-risk neuroblastoma.

CONCLUSION

Low expression of bivalent genes is a biomarker of worse outcome in patients with high-risk neuroblastoma.

Transcriptional Control of Subcutaneous Adipose Tissue by the Transcription Factor CTCF Modulates Heterogeneity in Fat Distribution in Women

Determining the mechanism driving body fat distribution will provide insights into obesity-related health risks. We used functional genomics tools to profile the epigenomic landscape to help infer the differential transcriptional potential of apple- and pear-shaped women's subcutaneous adipose-derived stem cells (ADSCs). We found that CCCTC-binding factor (CTCF) expression and its chromatin binding were increased in ADSCs from pear donors compared to those from apple donors. Interestingly, the pear enriched CTCF binding sites were located predominantly at the active transcription start sites (TSSs) of genes with active histone marks and YY1 motifs and were also associated with pear enriched RNAPII binding. In contrast, apple enriched CTCF binding sites were mainly found at intergenic regions and when identified at TSS, they were enriched with the bivalent chromatin signatures. Altogether, we provide evidence that CTCF plays an important role in differential regulation of subcutaneous ADSCs gene expression and may influence the development of apple vs. pear body shape.

ChIP-DIP: A multiplexed method for mapping hundreds of proteins to DNA uncovers diverse regulatory elements controlling gene expression

Gene expression is controlled by the dynamic localization of thousands of distinct regulatory proteins to precise regions of DNA. Understanding this cell-type specific process has been a goal of molecular biology for decades yet remains challenging because most current DNA-protein mapping methods study one protein at a time. To overcome this, we developed ChIP-DIP (ChIP Done In Parallel), a split-pool based method that enables simultaneous, genome-wide mapping of hundreds of diverse regulatory proteins in a single experiment. We demonstrate that ChIP-DIP generates highly accurate maps for all classes of DNA-associated proteins, including histone modifications, chromatin regulators, transcription factors, and RNA Polymerases. Using these data, we explore quantitative combinations of protein localization on genomic DNA to define distinct classes of regulatory elements and their functional activity. Our data demonstrate that ChIP-DIP enables the generation of 'consortium level', context-specific protein localization maps within any molecular biology lab.

CRISPR-Cas12a-integrated transgenes in genomic safe harbors retain high expression in human hematopoietic iPSC-derived lineages and primary cells

Discovery of genomic safe harbor sites (SHSs) is fundamental for multiple transgene integrations, such as reporter genes, chimeric antigen receptors (CARs), and safety switches, which are required for safe cell products for regenerative cell therapies and immunotherapies. Here we identified and characterized potential SHS in human cells. Using the CRISPR-MAD7 system, we integrated transgenes at these sites in induced pluripotent stem cells (iPSCs), primary T and natural killer (NK) cells, and Jurkat cell line, and demonstrated efficient and stable expression at these loci. Subsequently, we validated the differentiation potential of engineered iPSC toward CD34+ hematopoietic stem and progenitor cells (HSPCs), lymphoid progenitor cells (LPCs), and NK cells and showed that transgene expression was perpetuated in these lineages. Finally, we demonstrated that engineered iPSC-derived NK cells retained expression of a non-virally integrated anti-CD19 CAR, suggesting that several of the investigated SHSs can be used to engineer cells for adoptive immunotherapies.

Broad H3K4me3 Domain Is Associated with Spatial Coherence during Mammalian Embryonic Development

It is well known that the chromatin states play a major role in cell-fate decision and cell-identity maintenance; however, the spatial variation of chromatin states in situ remains poorly characterized. Here, by leveraging recently available spatial-CUT&Tag data, we systematically characterized the global spatial organization of the H3K4me3 profiles in a mouse embryo. Our analysis identified a subset of genes with spatially coherent H3K4me3 patterns, which together delineate the tissue boundaries. The spatially coherent genes are strongly enriched with tissue-specific transcriptional regulators. Remarkably, their corresponding genomic loci are marked by broad H3K4me3 domains, which is distinct from the typical H3K4me3 signature. Spatial transition across tissue boundaries is associated with continuous shortening of the broad H3K4me3 domains as well as expansion of H3K27me3 domains. Our analysis reveals a strong connection between the genomic and spatial variation of chromatin states, which may play an important role in embryonic development.

Minimum entropy framework identifies a novel class of genomic functional elements and reveals regulatory mechanisms at human disease loci

We introduce CoRE-BED, a framework trained using 19 epigenomic features in 33 major cell and tissue types to predict cell-type-specific regulatory function. CoRE-BED identifies nine functional classes de-novo, capturing both known and new regulatory categories. Notably, we describe a previously undercharacterized class that we term Development Associated Elements (DAEs), which are highly enriched in cell types with elevated regenerative potential and distinguished by the dual presence of either H3K4me2 and H3K9ac (an epigenetic signature associated with kinetochore assembly) or H3K79me3 and H4K20me1 (a signature associated with transcriptional pause release). Unlike bivalent promoters, which represent a transitory state between active and silenced promoters, DAEs transition directly to or from a non-functional state during stem cell differentiation and are proximal to highly expressed genes. CoRE-BED's interpretability facilitates causal inference and functional prioritization. Across 70 complex traits, distal insulators account for the largest mean proportion of SNP heritability (~49%) captured by the GWAS. Collectively, our results demonstrate the value of exploring non-conventional ways of regulatory classification that enrich for trait heritability, to complement existing approaches for cis-regulatory prediction.

Loss of cohesin regulator PDS5A reveals repressive role of Polycomb loops

Polycomb Repressive Complexes 1 and 2 (PRC1, PRC2) are conserved epigenetic regulators that promote transcriptional gene silencing. PRC1 and PRC2 converge on shared targets, catalyzing repressive histone modifications. Additionally, a subset of PRC1/PRC2 targets engage in long-range interactions whose functions in gene silencing are poorly understood. Using a CRISPR screen in mouse embryonic stem cells, we found that the cohesin regulator PDS5A links transcriptional silencing by Polycomb and 3D genome organization. PDS5A deletion impairs cohesin unloading and results in derepression of a subset of endogenous PRC1/PRC2 target genes. Importantly, derepression is not linked to loss of Polycomb chromatin domains. Instead, PDS5A removal causes aberrant cohesin activity leading to ectopic insulation sites, which disrupt the formation of ultra-long Polycomb loops. We show that these loops are important for robust silencing at a subset of PRC1/PRC2 target genes and that maintenance of cohesin-dependent genome architecture is critical for Polycomb regulation.

Electrostatic encoding of genome organization principles within single native nucleosomes

The eukaryotic genome, first packed into nucleosomes of about 150 bp around the histone core, is organized into euchromatin and heterochromatin, corresponding to the A and B compartments, respectively. Here, we asked if individual nucleosomes in vivo know where to go. That is, do mono-nucleosomes by themselves contain A/B compartment information, associated with transcription activity, in their biophysical properties? We purified native mono-nucleosomes to high monodispersity and used physiological concentrations of biological polyamines to determine their condensability. The chromosomal regions known to partition into A compartments have low condensability and vice versa. In silico chromatin polymer simulations using condensability as the only input showed that biophysical information needed to form compartments is all contained in single native nucleosomes and no other factors are needed. Condensability is also strongly anticorrelated with gene expression, and especially so near the promoter region and in a cell type dependent manner. Therefore, individual nucleosomes in the promoter know whether the gene is on or off, and that information is contained in their biophysical properties. Comparison with genetic and epigenetic features suggest that nucleosome condensability is a very meaningful axis onto which to project the high dimensional cellular chromatin state. Analysis of condensability using various condensing agents including those that are protein-based suggests that genome organization principle encoded into individual nucleosomes is electrostatic in nature. Polyamine depletion in mouse T cells, by either knocking out ornithine decarboxylase (ODC) or inhibiting ODC, results in hyperpolarized condensability, suggesting that when cells cannot rely on polyamines to translate biophysical properties of nucleosomes to control gene expression and 3D genome organization, they accentuate condensability contrast, which may explain dysfunction known to occur with polyamine deficiency.

Synaptic Mechanisms of Delay Eyeblink Classical Conditioning: AMPAR Trafficking and Gene Regulation in an In Vitro Model

An in vitro model of delay eyeblink classical conditioning was developed to investigate synaptic plasticity mechanisms underlying acquisition of associative learning. This was achieved by replacing real stimuli, such as an airpuff and tone, with patterned stimulation of the cranial nerves using an isolated brainstem preparation from turtle. Here, our primary findings regarding cellular and molecular mechanisms for learning acquisition using this unique approach are reviewed. The neural correlate of the in vitro eyeblink response is a replica of the actual behavior, and features of conditioned responses (CRs) resemble those observed in behavioral studies. Importantly, it was shown that acquisition of CRs did not require the intact cerebellum, but the appropriate timing did. Studies of synaptic mechanisms indicate that conditioning involves two stages of AMPA receptor (AMPAR) trafficking. Initially, GluA1-containing AMPARs are targeted to synapses followed later by replacement by GluA4 subunits that support CR expression. This two-stage process is regulated by specific signal transduction cascades involving PKA and PKC and is guided by distinct protein chaperones. The expression of the brain-derived neurotrophic factor (BDNF) protein is central to AMPAR trafficking and conditioning. BDNF gene expression is regulated by coordinated epigenetic mechanisms involving DNA methylation/demethylation and chromatin modifications that control access of promoters to transcription factors. Finally, a hypothesis is proposed that learning genes like BDNF are poised by dual chromatin features that allow rapid activation or repression in response to environmental stimuli. These in vitro studies have advanced our understanding of the cellular and molecular mechanisms that underlie associative learning.

Unraveling the battle for lysine: A review of the competition among post-translational modifications

Proteins play a critical role as key regulators in various biological systems, influencing crucial processes such as gene expression, cell cycle progression, and cellular proliferation. However, the functions of proteins can be further modified through post-translational modifications (PTMs), which expand their roles and contribute to disease progression when dysregulated. In this review, we delve into the methodologies employed for the characterization of PTMs, shedding light on the techniques and tools utilized to help unravel their complexity. Furthermore, we explore the prevalence of crosstalk and competition that occurs between different types of PTMs, specifically focusing on both histone and non-histone proteins. The intricate interplay between different modifications adds an additional layer of regulation to protein function and cellular processes. To gain insights into the competition for lysine residues among various modifications, computational systems such as MethylSight have been developed, allowing for a comprehensive analysis of the modification landscape. Additionally, we provide an overview of the exciting developments in the field of inhibitors or drugs targeting PTMs, highlighting their potential in combatting prevalent diseases. The discovery and development of drugs that modulate PTMs present promising avenues for therapeutic interventions, offering new strategies to address complex diseases. As research progresses in this rapidly evolving field, we anticipate remarkable advancements in our understanding of PTMs and their roles in health and disease, ultimately paving the way for innovative treatment approaches.

Dose rate dependent reduction in chromatin accessibility at transcriptional start sites long time after exposure to gamma radiation

Ionizing radiation (IR) impact cellular and molecular processes that require chromatin remodelling relevant for cellular integrity. However, the cellular implications of ionizing radiation (IR) delivered per time unit (dose rate) are still debated. This study investigates whether the dose rate is relevant for inflicting changes to the epigenome, represented by chromatin accessibility, or whether it is the total dose that is decisive. CBA/CaOlaHsd mice were whole-body exposed to either chronic low dose rate (2.5 mGy/h for 54 d) or the higher dose rates (10 mGy/h for 14 d and 100 mGy/h for 30 h) of gamma radiation (60Co, total dose: 3 Gy). Chromatin accessibility was analysed in liver tissue samples using Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-Seq), both one day after and over three months post-radiation (>100 d). The results show that the dose rate contributes to radiation-induced epigenomic changes in the liver at both sampling timepoints. Interestingly, chronic low dose rate exposure to a high total dose (3 Gy) did not inflict long-term changes to the epigenome. In contrast to the acute high dose rate given to the same total dose, reduced accessibility at transcriptional start sites (TSS) was identified in genes relevant for the DNA damage response and transcriptional activity. Our findings link dose rate to essential biological mechanisms that could be relevant for understanding long-term changes after ionizing radiation exposure. However, future studies are needed to comprehend the biological consequence of these findings.

An intronic LINE-1 regulates IFNAR1 expression in human immune cells

BACKGROUND

Despite their origins as selfish parasitic sequences, some transposons in the human genome have been co-opted to serve as regulatory elements, contributing to the evolution of transcriptional networks. Most well-characterized examples of transposon-derived regulatory elements derive from endogenous retroviruses (ERVs), due to the intrinsic regulatory activity of proviral long terminal repeat regions. However, one subclass of transposable elements, the Long Interspersed Nuclear Elements (LINEs), have been largely overlooked in the search for functional regulatory transposons, and considered to be broadly epigenetically repressed.

RESULTS

We examined the chromatin state of LINEs by analyzing epigenomic data from human immune cells. Many LINEs are marked by the repressive H3K9me3 modification, but a subset exhibits evidence of enhancer activity in human immune cells despite also showing evidence of epigenetic repression. We hypothesized that these competing forces of repressive and activating epigenetic marks might lead to inducible enhancer activity. We investigated a specific L1M2a element located within the first intron of Interferon Alpha/Beta Receptor 1 (IFNAR1). This element shows epigenetic signatures of B cell-specific enhancer activity, despite being repressed by the Human Silencing Hub (HUSH) complex. CRISPR deletion of the element in B lymphoblastoid cells revealed that the element acts as an enhancer that regulates both steady state and interferon-inducible expression of IFNAR1.

CONCLUSIONS

Our study experimentally demonstrates that an L1M2a element was co-opted to function as an interferon-inducible enhancer of IFNAR1, creating a feedback loop wherein IFNAR1 is transcriptionally upregulated by interferon signaling. This finding suggests that other LINEs may exhibit cryptic cell type-specific or context-dependent enhancer activity. LINEs have received less attention than ERVs in the effort to understand the contribution of transposons to the regulatory landscape of cellular genomes, but these are likely important, lineage-specific players in the rapid evolution of immune system regulatory networks and deserve further study.

H3K27me3 spreading organizes canonical PRC1 chromatin architecture to regulate developmental programs

Polycomb Repressive Complex 2 (PRC2)-mediated histone H3K27 tri-methylation (H3K27me3) recruits canonical PRC1 (cPRC1) to maintain heterochromatin. In early development, polycomb-regulated genes are connected through long-range 3D interactions which resolve upon differentiation. Here, we report that polycomb looping is controlled by H3K27me3 spreading and regulates target gene silencing and cell fate specification. Using glioma-derived H3 Lys-27-Met (H3K27M) mutations as tools to restrict H3K27me3 deposition, we show that H3K27me3 confinement concentrates the chromatin pool of cPRC1, resulting in heightened 3D interactions mirroring chromatin architecture of pluripotency, and stringent gene repression that maintains cells in progenitor states to facilitate tumor development. Conversely, H3K27me3 spread in pluripotent stem cells, following neural differentiation or loss of the H3K36 methyltransferase NSD1, dilutes cPRC1 concentration and dissolves polycomb loops. These results identify the regulatory principles and disease implications of polycomb looping and nominate histone modification-guided distribution of reader complexes as an important mechanism for nuclear compartment organization.

Highlights

The confinement of H3K27me3 at PRC2 nucleation sites without its spreading correlates with increased 3D chromatin interactions.The H3K27M oncohistone concentrates canonical PRC1 that anchors chromatin loop interactions in gliomas, silencing developmental programs.Stem and progenitor cells require factors promoting H3K27me3 confinement, including H3K36me2, to maintain cPRC1 loop architecture.The cPRC1-H3K27me3 interaction is a targetable driver of aberrant self-renewal in tumor cells.

Cellular reprogramming is driven by widespread rewiring of promoter-enhancer interactions

BACKGROUND

Long-range interactions between promoters and cis-regulatory elements, such as enhancers, play critical roles in gene regulation. However, the role of three-dimensional (3D) chromatin structure in orchestrating changes in transcriptional regulation during direct cell reprogramming is not fully understood.

RESULTS

Here, we performed integrated analyses of chromosomal architecture, epigenetics, and gene expression using Hi-C, promoter Capture Hi-C (PCHi-C), ChIP-seq, and RNA-seq during trans-differentiation of Pre-B cells into macrophages with a β-estradiol inducible C/EBPαER transgene. Within 1h of β-estradiol induction, C/EBPα translocated from the cytoplasm to the nucleus, binding to thousands of promoters and putative regulatory elements, resulting in the downregulation of Pre-B cell-specific genes and induction of macrophage-specific genes. Hi-C results were remarkably consistent throughout trans-differentiation, revealing only a small number of TAD boundary location changes, and A/B compartment switches despite significant changes in the expression of thousands of genes. PCHi-C revealed widespread changes in promoter-anchored loops with decreased interactions in parallel with decreased gene expression, and new and increased promoter-anchored interactions in parallel with increased expression of macrophage-specific genes.

CONCLUSIONS

Overall, our data demonstrate that C/EBPα-induced trans-differentiation involves few changes in genome architecture at the level of TADs and A/B compartments, in contrast with widespread reorganization of thousands of promoter-anchored loops in association with changes in gene expression and cell identity.

The role of noncoding RNAs in pancreatic birth defects

Congenital defects in the pancreas can cause severe health issues such as pancreatic cancer and diabetes which require lifelong treatment. Regenerating healthy pancreatic cells to replace malfunctioning cells has been considered a promising cure for pancreatic diseases including birth defects. However, such therapies are currently unavailable in the clinic. The developmental gene regulatory network underlying pancreatic development must be reactivated for in vivo regeneration and recapitulated in vitro for cell replacement therapy. Thus, understanding the mechanisms driving pancreatic development will pave the way for regenerative therapies. Pancreatic progenitor cells are the precursors of all pancreatic cells which use epigenetic changes to control gene expression during differentiation to generate all of the distinct pancreatic cell types. Epigenetic changes involving DNA methylation and histone modifications can be controlled by noncoding RNAs (ncRNAs). Indeed, increasing evidence suggests that ncRNAs are indispensable for proper organogenesis. Here, we summarize recent insight into the role of ncRNAs in the epigenetic regulation of pancreatic development. We further discuss how disruptions in ncRNA biogenesis and expression lead to developmental defects and diseases. This review summarizes in vivo data from animal models and in vitro studies using stem cell differentiation as a model for pancreatic development.

Histone H3 lysine 27 trimethylation suppresses jasmonate biosynthesis and signaling to affect male fertility under high temperature in cotton

High-temperature (HT) stress causes male sterility in crops, thus decreasing yields. To explore the possible contribution of histone modifications to male fertility under HT conditions, we defined the histone methylation landscape for the marks histone H3 lysine 27 trimethylation (H3K27me3) and histone H3 lysine 4 trimethylation (H3K4me3) by chromatin immunoprecipitation sequencing (ChIP-seq) in two differing upland cotton (Gossypium hirsutum) varieties. We observed a global disruption in H3K4me3 and H3K27me3 modifications, especially H3K27me3, in cotton anthers subjected to HT. HT affected the bivalent H3K4me3-H3K27me3 modification more than either monovalent modification. We determined that removal of H3K27me3 at the promoters of jasmonate-related genes increased their expression, maintaining male fertility under HT in the HT-tolerant variety at the anther dehiscence stage. Modulating jasmonate homeostasis or signaling resulted in an anther indehiscence phenotype under HT. Chemical suppression of H3K27me3 deposition increased jasmonic acid contents and maintained male fertility under HT. In summary, our study provides new insights into the regulation of male fertility by histone modifications under HT and suggests a potential strategy for improving cotton HT tolerance.

A mouse model with high clonal barcode diversity for joint lineage, transcriptomic, and epigenomic profiling in single cells

Cellular lineage histories and their molecular states encode fundamental principles of tissue development and homeostasis. Current lineage-recording mouse models have insufficient barcode diversity and single-cell lineage coverage for profiling tissues composed of millions of cells. Here, we developed DARLIN, an inducible Cas9 barcoding mouse line that utilizes terminal deoxynucleotidyl transferase (TdT) and 30 CRISPR target sites. DARLIN is inducible, generates massive lineage barcodes across tissues, and enables the detection of edited barcodes in ∼70% of profiled single cells. Using DARLIN, we examined fate bias within developing hematopoietic stem cells (HSCs) and revealed unique features of HSC migration. Additionally, we established a protocol for joint transcriptomic and epigenomic single-cell measurements with DARLIN and found that cellular clonal memory is associated with genome-wide DNA methylation rather than gene expression or chromatin accessibility. DARLIN will enable the high-resolution study of lineage relationships and their molecular signatures in diverse tissues and physiological contexts.

Modeling of covalent modifications of histones to estimate the binding affinity

Covalent histone modifications such as methylation, acetylation, phosphorylation, and other epigenetic modifications of the chromatin play an essential role in regulating eukaryotic cells of which most of these reactions are catalyzed by the enzymes. The binding energy of enzymes is often determined by experimental data via mathematical and statistical models due to specific modifications. Many theoretical models have been introduced to study histone modifications and reprogramming experiments in mammalian cells, in which all efforts in determining the affinity binding are essential part of the work. Here, we introduce a one-dimensional statistical Potts model to accurately determine the enzyme's binding free energy using the experimental data for various types of cells. We study the methylation of lysine 4 and 27 on histone H3 and suppose that each histone has one modification site with one of the seven states: H3K27me3, H3K27me2, H3K27me1, unmodified, H3K4me1, H3K4me2, and H3K4me3. Based on this model, the histone covalent modification is described. Moreover, by using simulation data, the histone's binding free energy and the energy of chromatin states are determined, when they are subject to changes from unmodified to active or repressive states, by finding the probability of the transition.

Corticosteroid-induced chromatin loop dynamics at the FKBP5 gene

FKBP5 is a 115-kb-long glucocorticoid-inducible gene implicated in psychiatric disorders. To investigate the complexities of chromatin interaction frequencies at the FKBP5 topologically associated domain (TAD), we deployed 15 one-to-all chromatin capture viewpoints near gene promoters, enhancers, introns, and CTCF-loop anchors. This revealed a "one-TAD-one-gene" structure encompassing the FKBP5 promoter and its enhancers. The FKBP5 promoter and its two glucocorticoid-stimulated enhancers roam the entire TAD while displaying subtle cell type-specific interactomes. The FKBP5 TAD consists of two nested CTCF loops that are coordinated by one CTCF site in the eighth intron of FKBP5 and another beyond its polyadenylation site, 61 kb further. Loop extension correlates with transcription increases through the intronic CTCF site. This is efficiently compensated for, since the short loop is restored even under high transcription regimes. The boundaries of the FKBP5 TAD consist of divergent CTCF site patterns, harbor multiple smaller genes, and are resilient to glucocorticoid stimulation. Interestingly, both FKBP5 TAD boundaries harbor H3K27me3-marked heterochromatin blocks that may reinforce them. We propose that cis-acting genetic and epigenetic polymorphisms underlying FKBP5 expression variation are likely to reside within a 240-kb region that consists of the FKBP5 TAD, its left sub-TAD, and both its boundaries.

Leveraging dominant-negative histone H3 K-to-M mutations to study chromatin during differentiation and development

Histone modifications are associated with regulation of gene expression that controls a vast array of biological processes. Often, these associations are drawn by correlating the genomic location of a particular histone modification with gene expression or phenotype; however, establishing a causal relationship between histone marks and biological processes remains challenging. Consequently, there is a strong need for experimental approaches to directly manipulate histone modifications. A class of mutations on the N-terminal tail of histone H3, lysine-to-methionine (K-to-M) mutations, was identified as dominant-negative inhibitors of histone methylation at their respective and specific residues. The dominant-negative nature of K-to-M mutants makes them a valuable tool for studying the function of specific methylation marks on histone H3. Here, we review recent applications of K-to-M mutations to understand the role of histone methylation during development and homeostasis. We highlight important advantages and limitations that require consideration when using K-to-M mutants, particularly in a developmental context.

Transposable Elements Are Co-opted as Oncogenic Regulatory Elements by Lineage-Specific Transcription Factors in Prostate Cancer

Transposable elements hold regulatory functions that impact cell fate determination by controlling gene expression. However, little is known about the transcriptional machinery engaged at transposable elements in pluripotent and mature versus oncogenic cell states. Through positional analysis over repetitive DNA sequences of H3K27ac chromatin immunoprecipitation sequencing data from 32 normal cell states, we report pluripotent/stem and mature cell state-specific "regulatory transposable elements." Pluripotent/stem elements are binding sites for pluripotency factors (e.g., NANOG, SOX2, OCT4). Mature cell elements are docking sites for lineage-specific transcription factors, including AR and FOXA1 in prostate epithelium. Expanding the analysis to prostate tumors, we identify a subset of regulatory transposable elements shared with pluripotent/stem cells, including Tigger3a. Using chromatin editing technology, we show how such elements promote prostate cancer growth by regulating AR transcriptional activity. Collectively, our results suggest that oncogenesis arises from lineage-specific transcription factors hijacking pluripotent/stem cell regulatory transposable elements.

SIGNIFICANCE

We show that oncogenesis relies on co-opting transposable elements from pluripotent stem cells as regulatory elements altering the recruitment of lineage-specific transcription factors. We further discover how co-option is dependent on active chromatin states with important implications for developing treatment options against drivers of oncogenesis across the repetitive DNA. This article is featured in Selected Articles from This Issue, p. 2293.

Noncoding rules of survival: epigenetic regulation of normal and malignant hematopoiesis

Hematopoiesis is an essential process for organismal development and homeostasis. Epigenetic regulation of gene expression is critical for stem cell self-renewal and differentiation in normal hematopoiesis. Increasing evidence shows that disrupting the balance between self-renewal and cell fate decisions can give rise to hematological diseases such as bone marrow failure and leukemia. Consequently, next-generation sequencing studies have identified various aberrations in histone modifications, DNA methylation, RNA splicing, and RNA modifications in hematologic diseases. Favorable outcomes after targeting epigenetic regulators during disease states have further emphasized their importance in hematological malignancy. However, these targeted therapies are only effective in some patients, suggesting that further research is needed to decipher the complexity of epigenetic regulation during hematopoiesis. In this review, an update on the impact of the epigenome on normal hematopoiesis, disease initiation and progression, and current therapeutic advancements will be discussed.

The mesodermal and myogenic specification of hESCs depend on ZEB1 and are inhibited by ZEB2

Human embryonic stem cells (hESCs) can differentiate into any cell lineage. Here, we report that ZEB1 and ZEB2 promote and inhibit mesodermal-to-myogenic specification of hESCs, respectively. Knockdown and/or overexpression experiments of ZEB1, ZEB2, or PAX7 in hESCs indicate that ZEB1 is required for hESC Nodal/Activin-mediated mesodermal specification and PAX7+ human myogenic progenitor (hMuP) generation, while ZEB2 inhibits these processes. ZEB1 downregulation induces neural markers, while ZEB2 downregulation induces mesodermal/myogenic markers. Mechanistically, ZEB1 binds to and transcriptionally activates the PAX7 promoter, while ZEB2 binds to and activates the promoter of the neural OTX2 marker. Transplanting ZEB1 or ZEB2 knocked down hMuPs into the muscles of a muscular dystrophy mouse model, showing that hMuP engraftment and generation of dystrophin-positive myofibers depend on ZEB1 and are inhibited by ZEB2. The mouse model results suggest that ZEB1 expression and/or downregulating ZEB2 in hESCs may also enhance hESC regenerative capacity for human muscular dystrophy therapy.

Unbiased transcription factor CRISPR screen identifies ZNF800 as master repressor of enteroendocrine differentiation

Enteroendocrine cells (EECs) are hormone-producing cells residing in the epithelium of stomach, small intestine (SI), and colon. EECs regulate aspects of metabolic activity, including insulin levels, satiety, gastrointestinal secretion, and motility. The generation of different EEC lineages is not completely understood. In this work, we report a CRISPR knockout screen of the entire repertoire of transcription factors (TFs) in adult human SI organoids to identify dominant TFs controlling EEC differentiation. We discovered ZNF800 as a master repressor for endocrine lineage commitment, which particularly restricts enterochromaffin cell differentiation by directly controlling an endocrine TF network centered on PAX4. Thus, organoid models allow unbiased functional CRISPR screens for genes that program cell fate.

A genetic screen identifies BEND3 as a regulator of bivalent gene expression and global DNA methylation

Epigenetic mechanisms are essential to establish and safeguard cellular identities in mammals. They dynamically regulate the expression of genes, transposable elements and higher-order chromatin structures. Consequently, these chromatin marks are indispensable for mammalian development and alterations often lead to disease, such as cancer. Bivalent promoters are especially important during differentiation and development. Here we used a genetic screen to identify new regulators of a bivalent repressed gene. We identify BEND3 as a regulator of hundreds of bivalent promoters, some of which it represses, and some of which it activates. We show that BEND3 is recruited to a CpG-containg consensus site that is present in multiple copies in many bivalent promoters. Besides having direct effect on the promoters it binds, the loss of BEND3 leads to genome-wide gains of DNA methylation, which are especially marked at regions normally protected by the TET enzymes. DNA hydroxymethylation is reduced in Bend3 mutant cells, possibly as consequence of altered gene expression leading to diminished alpha-ketoglutarate production, thus lowering TET activity. Our results clarify the direct and indirect roles of an important chromatin regulator, BEND3, and, more broadly, they shed light on the regulation of bivalent promoters.

Dicer-like 3a mediates intergenerational resistance against root-knot nematodes in rice via hormone responses

In a continuously changing and challenging environment, passing down the memory of encountered stress factors to offspring could provide an evolutionary advantage. In this study, we demonstrate the existence of "intergenerational acquired resistance" in the progeny of rice (Oryza sativa) plants attacked by the belowground parasitic nematode Meloidogyne graminicola. Transcriptome analyses revealed that genes involved in defense pathways are generally downregulated in progeny of nematode-infected plants under uninfected conditions but show a stronger induction upon nematode infection. This phenomenon was termed "spring loading" and depends on initial downregulation by the 24-nucleotide (nt) siRNA biogenesis gene dicer-like 3a (dcl3a) involved in the RNA-directed DNA methylation pathway. Knockdown of dcl3a led to increased nematode susceptibility and abolished intergenerational acquired resistance, as well as jasmonic acid/ethylene spring loading in the offspring of infected plants. The importance of ethylene signaling in intergenerational resistance was confirmed by experiments on a knockdown line of ethylene insensitive 2 (ein2b), which lacks intergenerational acquired resistance. Taken together, these data indicate a role for DCL3a in regulating plant defense pathways during both within-generation and intergenerational resistance against nematodes in rice.

Chromatin-modifying enzymes as modulators of nuclear size during lineage differentiation

The mechanism of nuclear size determination and alteration during normal lineage development and cancer pathologies which is not fully understood. As recently reported, chromatin modification can change nuclear morphology. Therefore, we screened a range of pharmacological chemical compounds that impact the activity of chromatin-modifying enzymes, in order to get a clue of the specific types of chromatin-modifying enzymes that remarkably effect nuclear size and shape. We found that interrupted activity of chromatin-modifying enzymes is associated with nuclear shape abnormalities. Furthermore, the activity of chromatin-modifying enzymes perturbs cell fate determination in cellular maintenance and lineage commitment. Our results indicated that chromatin-modifying enzyme regulates cell fate decision during lineage differentiation and is associate with nuclear size alteration.

Analyzing histone ChIP-seq data with a bin-based probability of being signal

Histone ChIP-seq is one of the primary methods for charting the cellular epigenomic landscape, the components of which play a critical regulatory role in gene expression. Analyzing the activity of regulatory elements across datasets and cell types can be challenging due to shifting peak positions and normalization artifacts resulting from, for example, differing read depths, ChIP efficiencies, and target sizes. Moreover, broad regions of enrichment seen in repressive histone marks often evade detection by commonly used peak callers. Here, we present a simple and versatile method for identifying enriched regions in ChIP-seq data that relies on estimating a gamma distribution fit to non-overlapping 5kB genomic bins to establish a global background. We use this distribution to assign a probability of being signal (PBS) between zero and one to each 5 kB bin. This approach, while lower in resolution than typical peak-calling methods, provides a straightforward way to identify enriched regions and compare enrichments among multiple datasets, by transforming the data to values that are universally normalized and can be readily visualized and integrated with downstream analysis methods. We demonstrate applications of PBS for both broad and narrow histone marks, and provide several illustrations of biological insights which can be gleaned by integrating PBS scores with downstream data types.

Epigenetic dynamics during capacitation of naïve human pluripotent stem cells

Human pluripotent stem cells (hPSCs) are of fundamental relevance in regenerative medicine. Naïve hPSCs hold promise to overcome some of the limitations of conventional (primed) hPSCs, including recurrent epigenetic anomalies. Naïve-to-primed transition (capacitation) follows transcriptional dynamics of human embryonic epiblast and is necessary for somatic differentiation from naïve hPSCs. We found that capacitated hPSCs are transcriptionally closer to postimplantation epiblast than conventional hPSCs. This prompted us to comprehensively study epigenetic and related transcriptional changes during capacitation. Our results show that CpG islands, gene regulatory elements, and retrotransposons are hotspots of epigenetic dynamics during capacitation and indicate possible distinct roles of specific epigenetic modifications in gene expression control between naïve and primed hPSCs. Unexpectedly, PRC2 activity appeared to be dispensable for the capacitation. We find that capacitated hPSCs acquire an epigenetic state similar to conventional hPSCs. Significantly, however, the X chromosome erosion frequently observed in conventional female hPSCs is reversed by resetting and subsequent capacitation.

Tissue-Specific Tumour Suppressor and Oncogenic Activities of the Polycomb-like Protein MTF2

The Polycomb repressive complex 2 (PRC2) is a conserved chromatin-remodelling complex that catalyses the trimethylation of histone H3 lysine 27 (H3K27me3), a mark associated with gene silencing. PRC2 regulates chromatin structure and gene expression during organismal and tissue development and tissue homeostasis in the adult. PRC2 core subunits are associated with various accessory proteins that modulate its function and recruitment to target genes. The multimeric composition of accessory proteins results in two distinct variant complexes of PRC2, PRC2.1 and PRC2.2. Metal response element-binding transcription factor 2 (MTF2) is one of the Polycomb-like proteins (PCLs) that forms the PRC2.1 complex. MTF2 is highly conserved, and as an accessory subunit of PRC2, it has important roles in embryonic stem cell self-renewal and differentiation, development, and cancer progression. Here, we review the impact of MTF2 in PRC2 complex assembly, catalytic activity, and spatiotemporal function. The emerging paradoxical evidence suggesting that MTF2 has divergent roles as either a tumour suppressor or an oncogene in different tissues merits further investigations. Altogether, our review illuminates the context-dependent roles of MTF2 in Polycomb group (PcG) protein-mediated epigenetic regulation. Its impact on disease paves the way for a deeper understanding of epigenetic regulation and novel therapeutic strategies.

SMCHD1 has separable roles in chromatin architecture and gene silencing that could be targeted in disease

The interplay between 3D chromatin architecture and gene silencing is incompletely understood. Here, we report a novel point mutation in the non-canonical SMC protein SMCHD1 that enhances its silencing capacity at endogenous developmental targets. Moreover, it also results in enhanced silencing at the facioscapulohumeral muscular dystrophy associated macrosatellite-array, D4Z4, resulting in enhanced repression of DUX4 encoded by this repeat. Heightened SMCHD1 silencing perturbs developmental Hox gene activation, causing a homeotic transformation in mice. Paradoxically, the mutant SMCHD1 appears to enhance insulation against other epigenetic regulators, including PRC2 and CTCF, while depleting long range chromatin interactions akin to what is observed in the absence of SMCHD1. These data suggest that SMCHD1's role in long range chromatin interactions is not directly linked to gene silencing or insulating the chromatin, refining the model for how the different levels of SMCHD1-mediated chromatin regulation interact to bring about gene silencing in normal development and disease.

The default and directed pathways of hepatoblast differentiation involve distinct epigenomic mechanisms

The effectiveness of multiomics analyses in defining cell differentiation pathways during development is ambiguous. During liver development, hepatoblasts follow a default or directed pathway to differentiate into hepatocytes or cholangiocytes, respectively, and this provides a practical model to address this issue. Our study discovered that promoter-associated histone modifications and chromatin accessibility dynamics, rather than enhancer-associated histone modifications, effectively delineated the "default vs. directed" process of hepatoblast differentiation. Histone H3K27me3 on bivalent promoters is associated with this asymmetric differentiation strategy in mice and humans. We demonstrated that Ezh2 and Jmjd3 exert opposing regulatory roles in hepatoblast-cholangiocyte differentiation. Additionally, active enhancers, regulated by P300, correlate with the development of both hepatocytes and cholangiocytes. This research proposes a model highlighting the division of labor between promoters and enhancers, with promoter-associated chromatin modifications governing the "default vs. directed" differentiation mode of hepatoblasts, whereas enhancer-associated modifications primarily dictate the progressive development processes of hepatobiliary lineages.

Analysis of regulatory network modules in hundreds of human stem cell lines reveals complex epigenetic and genetic factors contribute to pluripotency state differences between subpopulations

Stem cells exist in vitro in a spectrum of interconvertible pluripotent states. Analyzing hundreds of hiPSCs derived from different individuals, we show the proportions of these pluripotent states vary considerably across lines. We discovered 13 gene network modules (GNMs) and 13 regulatory network modules (RNMs), which were highly correlated with each other suggesting that the coordinated co-accessibility of regulatory elements in the RNMs likely underlied the coordinated expression of genes in the GNMs. Epigenetic analyses revealed that regulatory networks underlying self-renewal and pluripotency have a surprising level of complexity. Genetic analyses identified thousands of regulatory variants that overlapped predicted transcription factor binding sites and were associated with chromatin accessibility in the hiPSCs. We show that the master regulator of pluripotency, the NANOG-OCT4 Complex, and its associated network were significantly enriched for regulatory variants with large effects, suggesting that they may play a role in the varying cellular proportions of pluripotency states between hiPSCs. Our work captures the coordinated activity of tens of thousands of regulatory elements in hiPSCs and bins these elements into discrete functionally characterized regulatory networks, shows that regulatory elements in pluripotency networks harbor variants with large effects, and provides a rich resource for future pluripotent stem cell research.

NEUROD1 reinforces endocrine cell fate acquisition in pancreatic development

NEUROD1 is a transcription factor that helps maintain a mature phenotype of pancreatic β cells. Disruption of Neurod1 during pancreatic development causes severe neonatal diabetes; however, the exact role of NEUROD1 in the differentiation programs of endocrine cells is unknown. Here, we report a crucial role of the NEUROD1 regulatory network in endocrine lineage commitment and differentiation. Mechanistically, transcriptome and chromatin landscape analyses demonstrate that Neurod1 inactivation triggers a downregulation of endocrine differentiation transcription factors and upregulation of non-endocrine genes within the Neurod1-deficient endocrine cell population, disturbing endocrine identity acquisition. Neurod1 deficiency altered the H3K27me3 histone modification pattern in promoter regions of differentially expressed genes, which resulted in gene regulatory network changes in the differentiation pathway of endocrine cells, compromising endocrine cell potential, differentiation, and functional properties.

ChromGene: gene-based modeling of epigenomic data

Various computational approaches have been developed to annotate epigenomes on a per-position basis by modeling combinatorial and spatial patterns within epigenomic data. However, such annotations are less suitable for gene-based analyses. We present ChromGene, a method based on a mixture of learned hidden Markov models, to annotate genes based on multiple epigenomic maps across the gene body and flanks. We provide ChromGene assignments for over 100 cell and tissue types. We characterize the mixture components in terms of gene expression, constraint, and other gene annotations. The ChromGene method and annotations will provide a useful resource for gene-based epigenomic analyses.

G-quadruplexes associated with R-loops promote CTCF binding

CTCF is a critical regulator of genome architecture and gene expression that binds thousands of sites on chromatin. CTCF genomic localization is controlled by the recognition of a DNA sequence motif and regulated by DNA modifications. However, CTCF does not bind to all its potential sites in all cell types, raising the question of whether the underlying chromatin structure can regulate CTCF occupancy. Here, we report that R-loops facilitate CTCF binding through the formation of associated G-quadruplex (G4) structures. R-loops and G4s co-localize with CTCF at many genomic regions in mouse embryonic stem cells and promote CTCF binding to its cognate DNA motif in vitro. R-loop attenuation reduces CTCF binding in vivo. Deletion of a specific G4-forming motif in a gene reduces CTCF binding and alters gene expression. Conversely, chemical stabilization of G4s results in CTCF gains and accompanying alterations in chromatin organization, suggesting a pivotal role for G4 structures in reinforcing long-range genome interactions through CTCF.

Altered H3K4me3 profile at the TFAM promoter causes mitochondrial alterations in preadipocytes from first-degree relatives of type 2 diabetics

BACKGROUND

First-degree relatives of type 2 diabetics (FDR) exhibit a high risk of developing type 2 diabetes (T2D) and feature subcutaneous adipocyte hypertrophy, independent of obesity. In FDR, adipose cell abnormalities contribute to early insulin-resistance and are determined by adipocyte precursor cells (APCs) early senescence and impaired recruitment into the adipogenic pathway. Epigenetic mechanisms signal adipocyte differentiation, leading us to hypothesize that abnormal epigenetic modifications cause adipocyte dysfunction and enhance T2D risk. To test this hypothesis, we examined the genome-wide histone profile in APCs from the subcutaneous adipose tissue of healthy FDR.

RESULTS

Sequencing-data analysis revealed 2644 regions differentially enriched in lysine 4 tri-methylated H3-histone (H3K4me3) in FDR compared to controls (CTRL) with significant enrichment in mitochondrial-related genes. These included TFAM, which regulates mitochondrial DNA (mtDNA) content and stability. In FDR APCs, a significant reduction in H3K4me3 abundance at the TFAM promoter was accompanied by a reduction in TFAM mRNA and protein levels. FDR APCs also exhibited reduced mtDNA content and mitochondrial-genome transcription. In parallel, FDR APCs exhibited impaired differentiation and TFAM induction during adipogenesis. In CTRL APCs, TFAM-siRNA reduced mtDNA content, mitochondrial transcription and adipocyte differentiation in parallel with upregulation of the CDKN1A and ZMAT3 senescence genes. Furthermore, TFAM-siRNA significantly expanded hydrogen peroxide (H2O2)-induced senescence, while H2O2 did not affect TFAM expression.

CONCLUSIONS

Histone modifications regulate APCs ability to differentiate in mature cells, at least in part by modulating TFAM expression and affecting mitochondrial function. Reduced H3K4me3 enrichment at the TFAM promoter renders human APCs senescent and dysfunctional, increasing T2D risk.

Transcriptional pausing factor M1BP regulates cellular homeostasis by suppressing autophagy and apoptosis in Drosophila eye

During organogenesis cellular homeostasis plays a crucial role in patterning and growth. The role of promoter proximal pausing of RNA polymerase II, which regulates transcription of several developmental genes by GAGA factor or Motif 1 Binding Protein (M1BP), has not been fully understood in cellular homeostasis. Earlier, we reported that M1BP, a functional homolog of ZKSCAN3, regulates wingless and caspase-dependent cell death (apoptosis) in the Drosophila eye. Further, blocking apoptosis does not fully rescue the M1BPRNAi phenotype of reduced eye. Therefore, we looked for other possible mechanism(s). In a forward genetic screen, members of the Jun-amino-terminal-(NH2)-Kinase (JNK) pathway were identified. Downregulation of M1BP ectopically induces JNK, a pro-death pathway known to activate both apoptosis and caspase-independent (autophagy) cell death. Activation of JNK pathway components can enhance M1BPRNAi phenotype and vice-versa. Downregulation of M1BP ectopically induced JNK signaling, which leads to apoptosis and autophagy. Apoptosis and autophagy are regulated independently by their genetic circuitry. Here, we found that blocking either apoptosis or autophagy alone rescues the reduced eye phenotype of M1BP downregulation; whereas, blocking both apoptosis and autophagy together significantly rescues the M1BP reduced eye phenotype to near wild-type in nearly 85% progeny. This data suggests that the cellular homeostasis response demonstrated by two independent cell death mechanisms, apoptosis and autophagy, can be regulated by a common transcriptional pausing mechanism orchestrated by M1BP. Since these fundamental processes are conserved in higher organisms, this novel functional link between M1BP and regulation of both apoptosis and autophagy can be extrapolated to humans.

Shaping faces: genetic and epigenetic control of craniofacial morphogenesis

Major differences in facial morphology distinguish vertebrate species. Variation of facial traits underlies the uniqueness of human individuals, and abnormal craniofacial morphogenesis during development leads to birth defects that significantly affect quality of life. Studies during the past 40 years have advanced our understanding of the molecular mechanisms that establish facial form during development, highlighting the crucial roles in this process of a multipotent cell type known as the cranial neural crest cell. In this Review, we discuss recent advances in multi-omics and single-cell technologies that enable genes, transcriptional regulatory networks and epigenetic landscapes to be closely linked to the establishment of facial patterning and its variation, with an emphasis on normal and abnormal craniofacial morphogenesis. Advancing our knowledge of these processes will support important developments in tissue engineering, as well as the repair and reconstruction of the abnormal craniofacial complex.

Inhibition of Ezh2 redistributes bivalent domains within transcriptional regulators associated with WNT and Hedgehog pathways in osteoblasts

Bivalent epigenomic regulatory domains containing both activating histone 3 lysine 4 (H3K4me3) and repressive lysine 27 (H3K27me3) trimethylation are associated with key developmental genes. These bivalent domains repress transcription in the absence of differentiation signals but maintain regulatory genes in a poised state to allow for timely activation. Previous studies demonstrated that enhancer of zeste homolog 2 (Ezh2), a histone 3 lysine 27 (H3K27) methyltransferase, suppresses osteogenic differentiation and that inhibition of Ezh2 enhances commitment of osteoblast progenitors in vitro and bone formation in vivo. Here, we examined the mechanistic effects of Tazemetostat (EPZ6438), an Food and Drug Administration approved Ezh2 inhibitor for epithelioid sarcoma treatment, because this drug could potentially be repurposed to stimulate osteogenesis for clinical indications. We find that Tazemetostat reduces H3K27me3 marks in bivalent domains in enhancers required for bone formation and stimulates maturation of MC3T3 preosteoblasts. Furthermore, Tazemetostat activates bivalent genes associated with the Wingless/integrated (WNT), adenylyl cyclase (cAMP), and Hedgehog (Hh) signaling pathways based on transcriptomic (RNA-seq) and epigenomic (chromatin immunoprecipitation [ChIP]-seq) data. Functional analyses using selective pathway inhibitors and silencing RNAs demonstrate that the WNT and Hh pathways modulate osteogenic differentiation after Ezh2 inhibition. Strikingly, we show that loss of the Hh-responsive transcriptional regulator Gli1, but not Gli2, synergizes with Tazemetostat to accelerate osteoblast differentiation. These studies establish epigenetic cooperativity of Ezh2, Hh-Gli1 signaling, and bivalent regulatory genes in suppressing osteogenesis. Our findings may have important translational ramifications for anabolic applications requiring bone mass accrual and/or reversal of bone loss.

Refined read-out: The hUHRF1 Tandem-Tudor domain prefers binding to histone H3 tails containing K4me1 in the context of H3K9me2/3

UHRF1 is an essential chromatin protein required for DNA methylation maintenance, mammalian development, and gene regulation. We investigated the Tandem-Tudor domain (TTD) of human UHRF1 that is known to bind H3K9me2/3 histones and is a major driver of UHRF1 localization in cells. We verified binding to H3K9me2/3 but unexpectedly discovered stronger binding to H3 peptides and mononucleosomes containing K9me2/3 with additional K4me1. We investigated the combined binding of TTD to H3K4me1-K9me2/3 versus H3K9me2/3 alone, engineered mutants with specific and differential changes of binding, and discovered a novel read-out mechanism for H3K4me1 in an H3K9me2/3 context that is based on the interaction of R207 with the H3K4me1 methyl group and on counting the H-bond capacity of H3K4. Individual TTD mutants showed up to a 10,000-fold preference for the double-modified peptides, suggesting that after a conformational change, WT TTD could exhibit similar effects. The frequent appearance of H3K4me1-K9me2 regions in human chromatin demonstrated in our TTD chromatin pull-down and ChIP-western blot data suggests that it has specific biological roles. Chromatin pull-down of TTD from HepG2 cells and full-length murine UHRF1 ChIP-seq data correlate with H3K4me1 profiles indicating that the H3K4me1-K9me2/3 interaction of TTD influences chromatin binding of full-length UHRF1. We demonstrate the H3K4me1-K9me2/3 specific binding of UHRF1-TTD to enhancers and promoters of cell-type-specific genes at the flanks of cell-type-specific transcription factor binding sites, and provided evidence supporting an H3K4me1-K9me2/3 dependent and TTD mediated downregulation of these genes by UHRF1. All these findings illustrate the important physiological function of UHRF1-TTD binding to H3K4me1-K9me2/3 double marks in a cellular context.

Genome-wide high-throughput chromosome conformation capture analysis reveals hierarchical chromatin interactions during early somatic embryogenesis

Somatic embryogenesis (SE), like zygotic embryo development, is a progressive process. Early SE is the beginning of a switch from a somatic to an embryogenic state and is an important stage for initiating chromatin reprogramming of SE. Previous studies suggest that changes in chromatin accessibility occur during early SE, although information on the 3D structure of chromatin is not yet available. Here, we present a chromosome-level genome assembly of longan (Dimocarpus longan) using PacBio combined with high-through chromosome conformation capture scaffolding, which resulted in a 446 Mb genome assembly anchored onto 15 scaffolds. During early SE, chromatin was concentrated and then decondensed, and a large number of long terminal repeat retrotransposons (LTR-RTs) were enriched in the local chromatin interaction region, suggesting LTR-RTs were involved in chromatin reorganization. Early SE was accompanied by the transformation from A to B compartments, and the interactions between B compartments were enhanced. Results from chromatin accessibility, monomethylation of histone H3 at lysine 4 (H3K4me1) modification, and transcription analyses further revealed a gene regulatory network for cell wall thickening during SE. Particularly, we found that the H3K4me1 differential peak binding motif showed abnormal activation of ethylene response factor transcription factors and participation in SE. The chromosome-level genomic and multiomics analyses revealed the 3D conformation of chromatin during early SE, providing insight into the molecular mechanisms underlying cell wall thickening and the potential regulatory networks of TFs during early SE in D. longan. These results provide additional clues for revealing the molecular mechanisms of plant SE.

Circulating Biomarkers Associated with the Diagnosis and Prognosis of B-Cell Progenitor Acute Lymphoblastic Leukemia

Acute lymphoblastic leukemia (ALL) is a hematological disease characterized by the dysfunction of the hematopoietic system that leads to arrest at a specific stage of stem cells development, suppressing the average production of cellular hematologic components. BCP-ALL is a neoplasm of the B-cell lineage progenitor. BCP-ALL is caused and perpetuated by several mechanisms that provide the disease with its tumor potential and genetic and cytological characteristics. These pathological features are used for diagnosis and the prognostication of BCP-ALL. However, most of these paraclinical tools can only be obtained by bone marrow aspiration, which, as it is an invasive study, can delay the diagnosis and follow-up of the disease, in addition to the anesthetic risk it entails for pediatric patients. For this reason, it is crucial to find noninvasive and accessible ways to supply information concerning diagnosis, prognosis, and the monitoring of the disease, such as circulating biomarkers. In oncology, a biomarker is any measurable indicator that demonstrates the presence of malignancy, tumoral behavior, prognosis, or responses to treatments. This review summarizes circulating molecules associated with BCP-ALL with potential diagnostic value, classificatory capacity during monitoring specific clinic features of the disease, and/or capacity to identify each BCP-ALL stage regarding its evolution and outcome of the patients with BCP-ALL. In the same way, we provide and classify biomarkers that may be used in further studies focused on clinical approaches or therapeutic target identification for BCP-ALL.

Extravillous trophoblast cell lineage development is associated with active remodeling of the chromatin landscape

The extravillous trophoblast cell lineage is a key feature of placentation and successful pregnancy. Knowledge of transcriptional regulation driving extravillous trophoblast cell development is limited. Here, we map the transcriptome and epigenome landscape as well as chromatin interactions of human trophoblast stem cells and their transition into extravillous trophoblast cells. We show that integrating chromatin accessibility, long-range chromatin interactions, transcriptomic, and transcription factor binding motif enrichment enables identification of transcription factors and regulatory mechanisms critical for extravillous trophoblast cell development. We elucidate functional roles for TFAP2C, SNAI1, and EPAS1 in the regulation of extravillous trophoblast cell development. EPAS1 is identified as an upstream regulator of key extravillous trophoblast cell transcription factors, including ASCL2 and SNAI1 and together with its target genes, is linked to pregnancy loss and birth weight. Collectively, we reveal activation of a dynamic regulatory network and provide a framework for understanding extravillous trophoblast cell specification in trophoblast cell lineage development and human placentation.