An atlas of dynamic chromatin landscapes in mouse fetal development

Novel DNA methylation changes in mouse lungs associated with chronic smoking

Smoking is a potent cause of asthma exacerbations, chronic obstructive pulmonary disease (COPD) and many other health defects, and changes in DNA methylation (DNAm) have been identified as a potential link between smoking and these health outcomes. However, most studies of smoking and DNAm have been done using blood and other easily accessible tissues in humans, while evidence from more directly affected tissues such as the lungs is lacking. Here, we identified DNAm patterns in the lungs that are altered by smoking. We used an established mouse model to measure the effects of chronic smoke exposure first on lung phenotype immediately after smoking and then after a period of smoking cessation. Next, we determined whether our mouse model recapitulates previous DNAm patterns observed in smoking humans, specifically measuring DNAm at a candidate gene responsive to cigarette smoke, Cyp1a1. Finally, we carried out epigenome-wide DNAm analyses using the newly released Illumina mouse methylation microarrays. Our results recapitulate some of the phenotypes and DNAm patterns observed in human studies but reveal 32 differentially methylated genes specific to the lungs which have not been previously associated with smoking. The affected genes are associated with nicotine dependency, tumorigenesis and metastasis, immune cell dysfunction, lung function decline, and COPD. This research emphasizes the need to study CS-mediated DNAm signatures in directly affected tissues like the lungs, to fully understand mechanisms underlying CS-mediated health outcomes.

ATRX guards against aberrant differentiation in mesenchymal progenitor cells

Alterations in the tumor suppressor ATRX are recurrently observed in mesenchymal neoplasms. ATRX has multiple epigenetic functions including heterochromatin formation and maintenance and regulation of transcription through modulation of chromatin accessibility. Here, we show in murine mesenchymal progenitor cells (MPCs) that Atrx deficiency aberrantly activated mesenchymal differentiation programs. This includes adipogenic pathways where ATRX loss induced expression of adipogenic transcription factors and enhanced adipogenic differentiation in response to differentiation stimuli. These changes are linked to loss of heterochromatin near mesenchymal lineage genes together with increased chromatin accessibility and gains of active chromatin marks. We additionally observed depletion of H3K9me3 at transposable elements, which are derepressed including near mesenchymal genes where they could serve as regulatory elements. Finally, we demonstrated that loss of ATRX in a mesenchymal malignancy, undifferentiated pleomorphic sarcoma, results in similar epigenetic disruption and de-repression of transposable elements. Together, our results reveal a role for ATRX in maintaining epigenetic states and transcriptional repression in mesenchymal progenitors and tumor cells and in preventing aberrant differentiation in the progenitor context.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state Regulatory Potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbored distinctive transcription factor binding motifs that were similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we showed that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

A conserved molecular logic for neurogenesis to gliogenesis switch in the cerebral cortex

During development, neural stem cells in the cerebral cortex, also known as radial glial cells (RGCs), generate excitatory neurons, followed by production of cortical macroglia and inhibitory neurons that migrate to the olfactory bulb (OB). Understanding the mechanisms for this lineage switch is fundamental for unraveling how proper numbers of diverse neuronal and glial cell types are controlled. We and others recently showed that Sonic Hedgehog (Shh) signaling promotes the cortical RGC lineage switch to generate cortical oligodendrocytes and OB interneurons. During this process, cortical RGCs generate intermediate progenitor cells that express critical gliogenesis genes Ascl1, Egfr, and Olig2. The increased Ascl1 expression and appearance of Egfr+ and Olig2+ cortical progenitors are concurrent with the switch from excitatory neurogenesis to gliogenesis and OB interneuron neurogenesis in the cortex. While Shh signaling promotes Olig2 expression in the developing spinal cord, the exact mechanism for this transcriptional regulation is not known. Furthermore, the transcriptional regulation of Olig2 and Egfr has not been explored. Here, we show that in cortical progenitor cells, multiple regulatory programs, including Pax6 and Gli3, prevent precocious expression of Olig2, a gene essential for production of cortical oligodendrocytes and astrocytes. We identify multiple enhancers that control Olig2 expression in cortical progenitors and show that the mechanisms for regulating Olig2 expression are conserved between the mouse and human. Our study reveals evolutionarily conserved regulatory logic controlling the lineage switch of cortical neural stem cells.

A time-resolved multi-omics atlas of transcriptional regulation in response to high-altitude hypoxia across whole-body tissues

High-altitude hypoxia acclimatization requires whole-body physiological regulation in highland immigrants, but the underlying genetic mechanism has not been clarified. Here we use sheep as an animal model for low-to-high altitude translocation. We generate multi-omics data including whole-genome sequences, time-resolved bulk RNA-Seq, ATAC-Seq and single-cell RNA-Seq from multiple tissues as well as phenotypic data from 20 bio-indicators. We characterize transcriptional changes of all genes in each tissue, and examine multi-tissue temporal dynamics and transcriptional interactions among genes. Particularly, we identify critical functional genes regulating the short response to hypoxia in each tissue (e.g., PARG in the cerebellum and HMOX1 in the colon). We further identify TAD-constrained cis-regulatory elements, which suppress the transcriptional activity of most genes under hypoxia. Phenotypic and transcriptional evidence indicate that antenatal hypoxia could improve hypoxia tolerance in offspring. Furthermore, we provide time-series expression data of candidate genes associated with human mountain sickness (e.g., BMPR2) and high-altitude adaptation (e.g., HIF1A). Our study provides valuable resources and insights for future hypoxia-related studies in mammals.

Systematic epigenome editing captures the context-dependent instructive function of chromatin modifications

Chromatin modifications are linked with regulating patterns of gene expression, but their causal role and context-dependent impact on transcription remains unresolved. Here we develop a modular epigenome editing platform that programs nine key chromatin modifications, or combinations thereof, to precise loci in living cells. We couple this with single-cell readouts to systematically quantitate the magnitude and heterogeneity of transcriptional responses elicited by each specific chromatin modification. Among these, we show that installing histone H3 lysine 4 trimethylation (H3K4me3) at promoters can causally instruct transcription by hierarchically remodeling the chromatin landscape. We further dissect how DNA sequence motifs influence the transcriptional impact of chromatin marks, identifying switch-like and attenuative effects within distinct cis contexts. Finally, we examine the interplay of combinatorial modifications, revealing that co-targeted H3K27 trimethylation (H3K27me3) and H2AK119 monoubiquitination (H2AK119ub) maximizes silencing penetrance across single cells. Our precision-perturbation strategy unveils the causal principles of how chromatin modification(s) influence transcription and dissects how quantitative responses are calibrated by contextual interactions.

Alternative splicing of a chromatin modifier alters the transcriptional regulatory programs of stem cell maintenance and neuronal differentiation

Development of embryonic stem cells (ESCs) into neurons requires intricate regulation of transcription, splicing, and translation, but how these processes interconnect is not understood. We found that polypyrimidine tract binding protein 1 (PTBP1) controls splicing of DPF2, a subunit of BRG1/BRM-associated factor (BAF) chromatin remodeling complexes. Dpf2 exon 7 splicing is inhibited by PTBP1 to produce the DPF2-S isoform early in development. During neuronal differentiation, loss of PTBP1 allows exon 7 inclusion and DPF2-L expression. Different cellular phenotypes and gene expression programs were induced by these alternative DPF2 isoforms. We identified chromatin binding sites enriched for each DPF2 isoform, as well as sites bound by both. In ESC, DPF2-S preferential sites were bound by pluripotency factors. In neuronal progenitors, DPF2-S sites were bound by nuclear factor I (NFI), while DPF2-L sites were bound by CCCTC-binding factor (CTCF). DPF2-S sites exhibited enhancer modifications, while DPF2-L sites showed promoter modifications. Thus, alternative splicing redirects BAF complex targeting to impact chromatin organization during neuronal development.

Spatial enhancer activation influences inhibitory neuron identity during mouse embryonic development

The mammalian telencephalon contains distinct GABAergic projection neuron and interneuron types, originating in the germinal zone of the embryonic basal ganglia. How genetic information in the germinal zone determines cell types is unclear. Here we use a combination of in vivo CRISPR perturbation, lineage tracing and ChIP-sequencing analyses and show that the transcription factor MEIS2 favors the development of projection neurons by binding enhancer regions in projection-neuron-specific genes during mouse embryonic development. MEIS2 requires the presence of the homeodomain transcription factor DLX5 to direct its functional activity toward the appropriate binding sites. In interneuron precursors, the transcription factor LHX6 represses the MEIS2-DLX5-dependent activation of projection-neuron-specific enhancers. Mutations of Meis2 result in decreased activation of regulatory enhancers, affecting GABAergic differentiation. We propose a differential binding model where the binding of transcription factors at cis-regulatory elements determines differential gene expression programs regulating cell fate specification in the mouse ganglionic eminence.

Xenopus tropicalis osteoblast-specific open chromatin regions reveal promoters and enhancers involved in human skeletal phenotypes and shed light on early vertebrate evolution

While understanding the genetic underpinnings of osteogenesis has far-reaching implications for skeletal diseases and evolution, a comprehensive characterization of the osteoblastic regulatory landscape in non-mammalian vertebrates is still lacking. Here, we compared the ATAC-Seq profile of Xenopus tropicalis (Xt) osteoblasts to a variety of non mineralizing control tissues, and identified osteoblast-specific nucleosome free regions (NFRs) at 527 promoters and 6747 distal regions. Sequence analyses, Gene Ontology, RNA-Seq and ChIP-Seq against four key histone marks confirmed that the distal regions correspond to bona fide osteogenic transcriptional enhancers exhibiting a shared regulatory logic with mammals. We report 425 regulatory regions conserved with human and globally associated to skeletogenic genes. Of these, 35 regions have been shown to impact human skeletal phenotypes by GWAS, including one trps1 enhancer and the runx2 promoter, two genes which are respectively involved in trichorhinophalangeal syndrome type I and cleidocranial dysplasia. Intriguingly, 60 osteoblastic NFRs also align to the genome of the elephant shark, a species lacking osteoblasts and bone tissue. To tackle this paradox, we chose to focus on dlx5 because its conserved promoter, known to integrate regulatory inputs during mammalian osteogenesis, harbours an osteoblast-specific NFR in both frog and human. Hence, we show that dlx5 is expressed in Xt and elephant shark odontoblasts, supporting a common cellular and genetic origin of bone and dentine. Taken together, our work (i) unravels the Xt osteogenic regulatory landscape, (ii) illustrates how cross-species comparisons harvest data relevant to human biology and (iii) reveals that a set of genes including bnc2, dlx5, ebf3, mir199a, nfia, runx2 and zfhx4 drove the development of a primitive form of mineralized skeletal tissue deep in the vertebrate lineage.

Massively parallel reporter assays and mouse transgenic assays provide complementary information about neuronal enhancer activity

Genetic studies find hundreds of thousands of noncoding variants associated with psychiatric disorders. Massively parallel reporter assays (MPRAs) and in vivo transgenic mouse assays can be used to assay the impact of these variants. However, the relevance of MPRAs to in vivo function is unknown and transgenic assays suffer from low throughput. Here, we studied the utility of combining the two assays to study the impact of non-coding variants. We carried out an MPRA on over 50,000 sequences derived from enhancers validated in transgenic mouse assays and from multiple fetal neuronal ATAC-seq datasets. We also tested over 20,000 variants, including synthetic mutations in highly active neuronal enhancers and 177 common variants associated with psychiatric disorders. Variants with a high impact on MPRA activity were further tested in mice. We found a strong and specific correlation between MPRA and mouse neuronal enhancer activity including changes in neuronal enhancer activity in mouse embryos for variants with strong MPRA effects. Mouse assays also revealed pleiotropic variant effects that could not be observed in MPRA. Our work provides a large catalog of functional neuronal enhancers and variant effects and highlights the effectiveness of combining MPRAs and mouse transgenic assays.

Genetic architecture of trait variance in craniofacial morphology

The genetic architecture of trait variance has long been of interest in genetics and evolution. One of the earliest attempts to understand this architecture was presented in Lerner's Genetic Homeostasis (1954). Lerner proposed that heterozygotes should be better able to tolerate environmental perturbations because of functional differences between the alleles at a given locus, with each allele optimal for slightly different environments. This greater robustness to environmental variance, he argued, would result in smaller trait variance for heterozygotes. The evidence for Lerner's hypothesis has been inconclusive. To address this question using modern genomic methods, we mapped loci associated with differences in trait variance (vQTL) on 1,101 individuals from the F34 of an advanced intercross between LG/J and SM/J mice. We also mapped epistatic interactions for these vQTL in order to understand the influence of epistasis for the architecture of trait variance. We did not find evidence supporting Lerner's hypothesis, that heterozygotes tend to have smaller trait variances than homozygotes. We further show that the effects of most mapped loci on trait variance are produced by epistasis affecting trait means and that those epistatic effects account for about a half of the differences in genotypic-specific trait variances. Finally, we propose a model where the different interactions between the additive and dominance effects of the vQTL and their epistatic partners can explain Lerner's original observations but can also be extended to include other conditions where heterozygotes are not the least variable genotype.

PML-mediated nuclear loosening permits immunomodulation of mesenchymal stem/stromal cells under inflammatory conditions

Nuclear configuration plays a critical role in the compartmentalization of euchromatin and heterochromatin and the epigenetic regulation of gene expression. Under stimulation by inflammatory cytokines IFN-γ and TNF-α, human mesenchymal stromal cells (hMSCs) acquire a potent immunomodulatory function enabled by drastic induction of various effector genes, with some upregulated several magnitudes. However, whether the transcriptional upregulation of the immunomodulatory genes in hMSCs exposed to inflammatory cytokines is associated with genome-wide nuclear reconfiguration has not been explored. Here, we demonstrate that hMSCs undergo remarkable nuclear reconfiguration characterized by an enlargement of the nucleus, downregulation of LMNB1 and LMNA/C, decondensation of heterochromatin, and derepression of repetitive DNA. Interestingly, promyelocytic leukaemia-nuclear bodies (PML-NBs) were found to mediate the nuclear reconfiguration of hMSCs triggered by the inflammatory cytokines. Significantly, when PML was depleted, the immunomodulatory function of hMSCs conferred by cytokines was compromised, as reflected by the attenuated expression of effector molecules in hMSCs and their failure to block infiltration of immune cells to lipopolysaccharide (LPS)-induced acute lung injury. Our results indicate that the immunomodulatory function of hMSCs conferred by inflammatory cytokines requires PML-mediated chromatin loosening.

Increased enhancer-promoter interactions during developmental enhancer activation in mammals

Remote enhancers are thought to interact with their target promoters via physical proximity, yet the importance of this proximity for enhancer function remains unclear. Here we investigate the three-dimensional (3D) conformation of enhancers during mammalian development by generating high-resolution tissue-resolved contact maps for nearly a thousand enhancers with characterized in vivo activities in ten murine embryonic tissues. Sixty-one percent of developmental enhancers bypass their neighboring genes, which are often marked by promoter CpG methylation. The majority of enhancers display tissue-specific 3D conformations, and both enhancer-promoter and enhancer-enhancer interactions are moderately but consistently increased upon enhancer activation in vivo. Less than 14% of enhancer-promoter interactions form stably across tissues; however, these invariant interactions form in the absence of the enhancer and are likely mediated by adjacent CTCF binding. Our results highlight the general importance of enhancer-promoter physical proximity for developmental gene activation in mammals.

Synthetic reversed sequences reveal default genomic states

Pervasive transcriptional activity is observed across diverse species. The genomes of extant organisms have undergone billions of years of evolution, making it unclear whether these genomic activities represent effects of selection or 'noise'1-4. Characterizing default genome states could help understand whether pervasive transcriptional activity has biological meaning. Here we addressed this question by introducing a synthetic 101-kb locus into the genomes of Saccharomyces cerevisiae and Mus musculus and characterizing genomic activity. The locus was designed by reversing but not complementing human HPRT1, including its flanking regions, thus retaining basic features of the natural sequence but ablating evolved coding or regulatory information. We observed widespread activity of both reversed and native HPRT1 loci in yeast, despite the lack of evolved yeast promoters. By contrast, the reversed locus displayed no activity at all in mouse embryonic stem cells, and instead exhibited repressive chromatin signatures. The repressive signature was alleviated in a locus variant lacking CpG dinucleotides; nevertheless, this variant was also transcriptionally inactive. These results show that synthetic genomic sequences that lack coding information are active in yeast, but inactive in mouse embryonic stem cells, consistent with a major difference in 'default genomic states' between these two divergent eukaryotic cell types, with implications for understanding pervasive transcription, horizontal transfer of genetic information and the birth of new genes.

The growth factor EPIREGULIN promotes basal progenitor cell proliferation in the developing neocortex

Neocortex expansion during evolution is linked to higher numbers of neurons, which are thought to result from increased proliferative capacity and neurogenic potential of basal progenitor cells during development. Here, we show that EREG, encoding the growth factor EPIREGULIN, is expressed in the human developing neocortex and in gorilla cerebral organoids, but not in the mouse neocortex. Addition of EPIREGULIN to the mouse neocortex increases proliferation of basal progenitor cells, whereas EREG ablation in human cortical organoids reduces proliferation in the subventricular zone. Treatment of cortical organoids with EPIREGULIN promotes a further increase in proliferation of gorilla but not of human basal progenitor cells. EPIREGULIN competes with the epidermal growth factor (EGF) to promote proliferation, and inhibition of the EGF receptor abrogates the EPIREGULIN-mediated increase in basal progenitor cells. Finally, we identify putative cis-regulatory elements that may contribute to the observed inter-species differences in EREG expression. Our findings suggest that species-specific regulation of EPIREGULIN expression may contribute to the increased neocortex size of primates by providing a tunable pro-proliferative signal to basal progenitor cells in the subventricular zone.

Omics-based construction of regulatory variants can be applied to help decipher pig liver-related traits

Genetic variants can influence complex traits by altering gene expression through changes to regulatory elements. However, the genetic variants that affect the activity of regulatory elements in pigs are largely unknown, and the extent to which these variants influence gene expression and contribute to the understanding of complex phenotypes remains unclear. Here, we annotate 90,991 high-quality regulatory elements using acetylation of histone H3 on lysine 27 (H3K27ac) ChIP-seq of 292 pig livers. Combined with genome resequencing and RNA-seq data, we identify 28,425 H3K27ac quantitative trait loci (acQTLs) and 12,250 expression quantitative trait loci (eQTLs). Through the allelic imbalance analysis, we validate two causative acQTL variants in independent datasets. We observe substantial sharing of genetic controls between gene expression and H3K27ac, particularly within promoters. We infer that 46% of H3K27ac exhibit a concomitant rather than causative relationship with gene expression. By integrating GWAS, eQTLs, acQTLs, and transcription factor binding prediction, we further demonstrate their application, through metabolites dulcitol, phosphatidylcholine (PC) (16:0/16:0) and published phenotypes, in identifying likely causal variants and genes, and discovering sub-threshold GWAS loci. We provide insight into the relationship between regulatory elements and gene expression, and the genetic foundation for dissecting the molecular mechanism of phenotypes.

Epigenomic tomography for probing spatially defined chromatin state in the brain

Spatially resolved epigenomic profiling is critical for understanding biology in the mammalian brain. Single-cell spatial epigenomic assays were developed recently for this purpose, but they remain costly and labor intensive for examining brain tissues across substantial dimensions and surveying a collection of brain samples. Here, we demonstrate an approach, epigenomic tomography, that maps spatial epigenomes of mouse brain at the scale of centimeters. We individually profiled neuronal and glial fractions of mouse neocortex slices with 0.5 mm thickness. Tri-methylation of histone 3 at lysine 27 (H3K27me3) or acetylation of histone 3 at lysine 27 (H3K27ac) features across these slices were grouped into clusters based on their spatial variation patterns to form epigenomic brain maps. As a proof of principle, our approach reveals striking dynamics in the frontal cortex due to kainic-acid-induced seizure, linked with transmembrane ion transporters, exocytosis of synaptic vesicles, and secretion of neurotransmitters. Epigenomic tomography provides a powerful and cost-effective tool for characterizing brain disorders based on the spatial epigenome.

The characterization of protein lactylation in relation to cardiac metabolic reprogramming in neonatal mouse hearts

In mammals, the neonatal heart can regenerate upon injury within a short time after birth, while adults lose this ability. Metabolic reprogramming has been demonstrated to be critical for cardiomyocyte proliferation in the neonatal heart. Here, we reveal that cardiac metabolic reprogramming could be regulated by altering global protein lactylation. By performing 4D label-free proteomics and lysine lactylation (Kla) omics analyses in mouse hearts at postnatal days 1, 5, and 7, 2297 Kla sites from 980 proteins are identified, among which 1262 Kla sites from 409 proteins are quantified. Functional clustering analysis reveals that the proteins with altered Kla sites are mainly involved in metabolic processes. The expression and Kla levels of proteins in glycolysis show a positive correlation while a negative correlation in fatty acid oxidation. Furthermore, we verify the Kla levels of several differentially modified proteins, including ACAT1, ACADL, ACADVL, PFKM, PKM, and NPM1. Overall, our study reports a comprehensive Kla map in the neonatal mouse heart, which will help to understand the regulatory network of metabolic reprogramming and cardiac regeneration.

Atlas of cardiac endothelial cell enhancer elements linking the mineralocorticoid receptor to pathological gene expression

Endothelial cells play crucial roles in physiology and are increasingly recognized as therapeutic targets in cardiovascular disease. Here, we analyzed the regulatory landscape of cardiac endothelial cells by assessing chromatin accessibility, histone modifications, and 3D chromatin organization and confirmed the functional relevance of enhancer-promoter interactions by CRISPRi-mediated enhancer silencing. We used this dataset to explore mechanisms of transcriptional regulation in cardiovascular disease and compared six different experimental models of heart failure, hypertension, or diabetes. Enhancers that regulate gene expression in diseased endothelial cells were enriched with binding sites for a distinct set of transcription factors, including the mineralocorticoid receptor (MR), a known drug target in heart failure and hypertension. For proof of concept, we applied endothelial cell-specific MR deletion in mice to confirm MR-dependent gene expression and predicted direct MR target genes. Overall, we have compiled here a comprehensive atlas of cardiac endothelial cell enhancer elements that provides insight into the role of transcription factors in cardiovascular disease.

Dynamic enhancer landscapes in human craniofacial development

The genetic basis of human facial variation and craniofacial birth defects remains poorly understood. Distant-acting transcriptional enhancers control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development. However, a lack of accurate maps of the genomic locations and cell type-resolved activities of craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combine histone modification, chromatin accessibility, and gene expression profiling of human craniofacial development with single-cell analyses of the developing mouse face to define the regulatory landscape of facial development at tissue- and single cell-resolution. We provide temporal activity profiles for 14,000 human developmental craniofacial enhancers. We find that 56% of human craniofacial enhancers share chromatin accessibility in the mouse and we provide cell population- and embryonic stage-resolved predictions of their in vivo activity. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

Mapping catalytically engaged TOP2B in neurons reveals the principles of topoisomerase action within the genome

We trapped catalytically engaged topoisomerase IIβ (TOP2B) in covalent DNA cleavage complexes (TOP2Bccs) and mapped their positions genome-wide in cultured mouse cortical neurons. We report that TOP2Bcc distribution varies with both nucleosome and compartmental chromosome organization. While TOP2Bccs in gene bodies correlate with their level of transcription, highly expressed genes that lack the usually associated chromatin marks, such as H3K36me3, show reduced TOP2Bccs, suggesting that histone posttranslational modifications regulate TOP2B activity. Promoters with high RNA polymerase II occupancy show elevated TOP2B chromatin immunoprecipitation sequencing signals but low TOP2Bccs, indicating that TOP2B catalytic engagement is curtailed at active promoters. Surprisingly, either poisoning or inhibiting TOP2B increases nascent transcription at most genes and enhancers but reduces transcription within long genes. These effects are independent of transcript length and instead correlate with the presence of intragenic enhancers. Together, these results clarify how cells modulate the catalytic engagement of topoisomerases to affect transcription.

Massively parallel disruption of enhancers active in human neural stem cells

Changes in gene regulation have been linked to the expansion of the human cerebral cortex and to neurodevelopmental disorders, potentially by altering neural progenitor proliferation. However, the effects of genetic variation within regulatory elements on neural progenitors remain obscure. We use sgRNA-Cas9 screens in human neural stem cells (hNSCs) to disrupt 10,674 genes and 26,385 conserved regions in 2,227 enhancers active in the developing human cortex and determine effects on proliferation. Genes with proliferation phenotypes are associated with neurodevelopmental disorders and show biased expression in specific fetal human brain neural progenitor populations. Although enhancer disruptions overall have weaker effects than gene disruptions, we identify enhancer disruptions that severely alter hNSC self-renewal. Disruptions in human accelerated regions, implicated in human brain evolution, also alter proliferation. Integrating proliferation phenotypes with chromatin interactions reveals regulatory relationships between enhancers and their target genes contributing to neurogenesis and potentially to human cortical evolution.

Data integration and inference of gene regulation using single-cell temporal multimodal data with scTIE

Single-cell technologies offer unprecedented opportunities to dissect gene regulatory mechanisms in context-specific ways. Although there are computational methods for extracting gene regulatory relationships from scRNA-seq and scATAC-seq data, the data integration problem, essential for accurate cell type identification, has been mostly treated as a standalone challenge. Here we present scTIE, a unified method that integrates temporal multimodal data and infers regulatory relationships predictive of cellular state changes. scTIE uses an autoencoder to embed cells from all time points into a common space by using iterative optimal transport, followed by extracting interpretable information to predict cell trajectories. Using a variety of synthetic and real temporal multimodal data sets, we show scTIE achieves effective data integration while preserving more biological signals than existing methods, particularly in the presence of batch effects and noise. Furthermore, on the exemplar multiome data set we generated from differentiating mouse embryonic stem cells over time, we show scTIE captures regulatory elements highly predictive of cell transition probabilities, providing new potentials to understand the regulatory landscape driving developmental processes.

Mapping Cell Atlases at the Single-Cell Level

Recent advancements in single-cell technologies have led to rapid developments in the construction of cell atlases. These atlases have the potential to provide detailed information about every cell type in different organisms, enabling the characterization of cellular diversity at the single-cell level. Global efforts in developing comprehensive cell atlases have profound implications for both basic research and clinical applications. This review provides a broad overview of the cellular diversity and dynamics across various biological systems. In addition, the incorporation of machine learning techniques into cell atlas analyses opens up exciting prospects for the field of integrative biology.

Evolutionary innovation in conserved regulatory elements across the mammalian tree of life

Transcriptional enhancers orchestrate cell type- and time point-specific gene expression programs. Evolution of enhancer sequences can alter target gene expression without causing detrimental misexpression in other contexts. It has long been thought that this modularity allows evolutionary changes in enhancers to escape pleiotropic constraints, which is especially important for evolutionary constrained developmental patterning genes. However, there is still little data supporting this hypothesis. Here we identified signatures of accelerated evolution in conserved enhancer elements across the mammalian phylogeny. We found that pleiotropic genes involved in gene regulatory and developmental processes were enriched for accelerated sequence evolution within their enhancer elements. These genes were associated with an excess number of enhancers compared to other genes, and due to this they exhibit a substantial degree of sequence acceleration over all their enhancers combined. We provide evidence that sequence acceleration is associated with turnover of regulatory function. We studied one acceleration event in depth and found that its sequence evolution led to the emergence of a new enhancer activity domain that may be involved in the evolution of digit reduction in hoofed mammals. Our results provide tangible evidence that enhancer evolution has been a frequent contributor to modifications involving constrained developmental signaling genes in mammals.

The Dynamics of Histone Modifications during Mammalian Zygotic Genome Activation

Mammalian fertilization initiates the reprogramming of oocytes and sperm, forming a totipotent zygote. During this intricate process, the zygotic genome undergoes a maternal-to-zygotic transition (MZT) and subsequent zygotic genome activation (ZGA), marking the initiation of transcriptional control and gene expression post-fertilization. Histone modifications are pivotal in shaping cellular identity and gene expression in many mammals. Recent advances in chromatin analysis have enabled detailed explorations of histone modifications during ZGA. This review delves into conserved and unique regulatory strategies, providing essential insights into the dynamic changes in histone modifications and their variants during ZGA in mammals. The objective is to explore recent advancements in leading mechanisms related to histone modifications governing this embryonic development phase in depth. These considerations will be useful for informing future therapeutic approaches that target epigenetic regulation in diverse biological contexts. It will also contribute to the extensive areas of evolutionary and developmental biology and possibly lay the foundation for future research and discussion on this seminal topic.

HOPX-associated molecular programs control cardiomyocyte cell states underpinning cardiac structure and function

Genomic regulation of cardiomyocyte differentiation is central to heart development and function. This study uses genetic loss-of-function human-induced pluripotent stem cell-derived cardiomyocytes to evaluate the genomic regulatory basis of the non-DNA-binding homeodomain protein HOPX. We show that HOPX interacts with and controls cardiac genes and enhancer networks associated with diverse aspects of heart development. Using perturbation studies in vitro, we define how upstream cell growth and proliferation control HOPX transcription to regulate cardiac gene programs. We then use cell, organoid, and zebrafish regeneration models to demonstrate that HOPX-regulated gene programs control cardiomyocyte function in development and disease. Collectively, this study mechanistically links cell signaling pathways as upstream regulators of HOPX transcription to control gene programs underpinning cardiomyocyte identity and function.

Atlas of fetal metabolism during mid-to-late gestation and diabetic pregnancy

Mounting evidence suggests metabolism instructs stem cell fate decisions. However, how fetal metabolism changes during development and how altered maternal metabolism shapes fetal metabolism remain unexplored. We present a descriptive atlas of in vivo fetal murine metabolism during mid-to-late gestation in normal and diabetic pregnancy. Using 13C-glucose and liquid chromatography-mass spectrometry (LC-MS), we profiled the metabolism of fetal brains, hearts, livers, and placentas harvested from pregnant dams between embryonic days (E)10.5 and 18.5. Our analysis revealed metabolic features specific to a hyperglycemic environment and signatures that may denote developmental transitions during euglycemic development. We observed sorbitol accumulation in fetal tissues and altered neurotransmitter levels in fetal brains isolated from hyperglycemic dams. Tracing 13C-glucose revealed disparate fetal nutrient sourcing depending on maternal glycemic states. Regardless of glycemic state, histidine-derived metabolites accumulated in late-stage fetal tissues. Our rich dataset presents a comprehensive overview of in vivo fetal tissue metabolism and alterations due to maternal hyperglycemia.

Isolated loss of the AUTS2 long isoform, brain-wide or targeted to Calbindin-lineage cells, generates a specific suite of brain, behavioral, and molecular pathologies

Rearrangements within the AUTS2 region are associated with a rare syndromic disorder with intellectual disability, developmental delay, and behavioral abnormalities as core features. In addition, smaller regional variants are linked to wide range of neuropsychiatric disorders, underscoring the gene's essential role in brain development. Like many essential neurodevelopmental genes, AUTS2 is large and complex, generating distinct long (AUTS2-l) and short (AUTS2-s) protein isoforms from alternative promoters. Although evidence suggests unique isoform functions, the contributions of each isoform to specific AUTS2-linked phenotypes have not been clearly resolved. Furthermore, Auts2 is widely expressed across the developing brain, but cell populations most central to disease presentation have not been determined. In this study, we focused on the specific roles of AUTS2-l in brain development, behavior, and postnatal brain gene expression, showing that brain-wide AUTS2-l ablation leads to specific subsets of the recessive pathologies associated with mutations in 3' exons (exons 8-19) that disrupt both major isoforms. We identify downstream genes that could explain expressed phenotypes including hundreds of putative direct AUTS2-l target genes. Furthermore, in contrast to 3' Auts2 mutations which lead to dominant hypoactivity, AUTS2-l loss-of-function is associated with dominant hyperactivity and repetitive behaviors, phenotypes exhibited by many human patients. Finally, we show that AUTS2-l ablation in Calbindin 1-expressing cell lineages is sufficient to yield learning/memory deficits and hyperactivity with abnormal dentate gyrus granule cell maturation, but not other phenotypic effects. These data provide new clues to in vivo AUTS2-l functions and novel information relevant to genotype-phenotype correlations in the human AUTS2 region.

Cell-type differential targeting of SETDB1 prevents aberrant CTCF binding, chromatin looping, and cis-regulatory interactions

SETDB1 is an essential histone methyltransferase that deposits histone H3 lysine 9 trimethylation (H3K9me3) to transcriptionally repress genes and repetitive elements. The function of differential H3K9me3 enrichment between cell-types remains unclear. Here, we demonstrate mutual exclusivity of H3K9me3 and CTCF across mouse tissues from different developmental timepoints. We analyze SETDB1 depleted cells and discover that H3K9me3 prevents aberrant CTCF binding independently of DNA methylation and H3K9me2. Such sites are enriched with SINE B2 retrotransposons. Moreover, analysis of higher-order genome architecture reveals that large chromatin structures including topologically associated domains and subnuclear compartments, remain intact in SETDB1 depleted cells. However, chromatin loops and local 3D interactions are disrupted, leading to transcriptional changes by modifying pre-existing chromatin landscapes. Specific genes with altered expression show differential interactions with dysregulated cis-regulatory elements. Collectively, we find that cell-type specific targets of SETDB1 maintain cellular identities by modulating CTCF binding, which shape nuclear architecture and transcriptomic networks.

The Cross-Regulation Between Set1, Clr4, and Lsd1/2 in Schizosaccharomyces pombe

Eukaryotic chromatin is organized into either silenced heterochromatin or relaxed euchromatin regions, which controls the accessibility of transcriptional machinery and thus regulates gene expression. In fission yeast, Schizosaccharomyces pombe, Set1 is the sole H3K4 methyltransferase and is mainly enriched at the promoters of actively transcribed genes. In contrast, Clr4 methyltransferase initiates H3K9 methylation, which has long been regarded as a hallmark of heterochromatic silencing. Lsd1 and Lsd2 are two highly conserved H3K4 and H3K9 demethylases. As these histone-modifying enzymes perform critical roles in maintaining histone methylation patterns and, consequently, gene expression profiles, cross-regulations among these enzymes are part of the complex regulatory networks. Thus, elucidating the mechanisms that govern their signaling and mutual regulations remains crucial. Here, we demonstrated that C-terminal truncation mutants, lsd1-ΔHMG and lsd2-ΔC, do not compromise the integrity of the Lsd1/2 complex but impair their chromatin-binding capacity at the promoter region of target genomic loci. We identified protein-protein interactions between Lsd1/2 and Raf2 or Swd2, which are the subunits of the Clr4 complex (CLRC) and Set1-associated complex (COMPASS), respectively. We showed that Clr4 and Set1 modulate the protein levels of Lsd1 and Lsd2 in opposite ways through the ubiquitin-proteasome-dependent pathway. During heat stress, the protein levels of Lsd1 and Lsd2 are upregulated in a Set1-dependent manner. The increase in protein levels is crucial for differential gene expression under stress conditions. Together, our results support a cross-regulatory model by which Set1 and Clr4 methyltransferases control the protein levels of Lsd1/2 demethylases to shape the dynamic chromatin landscape.

Ablation of placental REST deregulates fetal brain metabolism and impacts gene expression of the offspring brain at the postnatal and adult stages

In this study, the transcriptional repressor REST (Repressor Element 1 Silencing Transcription factor) was ablated in the mouse placenta to investigate molecular and cellular impacts on the offspring brain at different life stages. Ablation of placental REST deregulated several brain metabolites, including glucose and lactate that fuel brain energy, vitamin C (ascorbic acid) that functions in the epigenetic programming of the brain during postnatal development, and glutamate and creatine that help the brain to respond to stress conditions during adult life. Bulk RNA-seq analysis showed that a lack of placental REST persistently altered multiple transport genes, including those related to oxygen transportation in the offspring brain. While metabolic genes were impacted in the postnatal brain, different stress response genes were activated in the adult brain. DNA methylation was also impacted in the adult brain due to the loss of placental REST, but in a sex-biased manner. Single-nuclei RNA-seq analysis showed that specific cell types of the brain, particularly those of the choroid plexus and ependyma, which play critical roles in producing cerebrospinal fluid and maintaining metabolic homeostasis, were significantly impacted due to the loss of placental REST. These cells showed significant differential expression of genes associated with the metabotropic (G coupled protein) and ionotropic (ligand-gated ion channel) glutamate receptors, suggesting an impact of ablation of placental REST on the glutamatergic signaling of the offspring brain. The study expands our understanding of placental influences on the offspring brain.

Chromatin state transitions in the Drosophila intestinal lineage identify principles of cell-type specification

Tissue homeostasis relies on rewiring of stem cell transcriptional programs into those of differentiated cells. Here, we investigate changes in chromatin occurring in a bipotent adult stem cells. Combining mapping of chromatin-associated factors with statistical modeling, we identify genome-wide transitions during differentiation in the adult Drosophila intestinal stem cell (ISC) lineage. Active, stem-cell-enriched genes transition to a repressive heterochromatin protein-1-enriched state more prominently in enteroendocrine cells (EEs) than in enterocytes (ECs), in which the histone H1-enriched Black state is preeminent. In contrast, terminal differentiation genes associated with metabolic functions follow a common path from a repressive, primed, histone H1-enriched Black state in ISCs to active chromatin states in EE and EC cells. Furthermore, we find that lineage priming has an important function in adult ISCs, and we identify histone H1 as a mediator of this process. These data define underlying principles of chromatin changes during adult multipotent stem cell differentiation.

Sequence characteristics and an accurate model of abundant hyperactive loci in the human genome

Enhancers and promoters are classically considered to be bound by a small set of TFs in a sequence-specific manner. This assumption has come under increasing skepticism as the datasets of ChIP-seq assays of TFs have expanded. In particular, high-occupancy target (HOT) loci attract hundreds of TFs with seemingly no detectable correlation between ChIP-seq peaks and DNA-binding motif presence. Here, we used a set of 1,003 TF ChIP-seq datasets (HepG2, K562, H1) to analyze the patterns of ChIP-seq peak co-occurrence in combination with functional genomics datasets. We identified 43,891 HOT loci forming at the promoter (53%) and enhancer (47%) regions. HOT promoters regulate housekeeping genes, whereas HOT enhancers are involved in tissue-specific process regulation. HOT loci form the foundation of human super-enhancers and evolve under strong negative selection, with some of these loci being located in ultraconserved regions. Sequence-based classification analysis of HOT loci suggested that their formation is driven by the sequence features, and the density of mapped ChIP-seq peaks across TF-bound loci correlates with sequence features and the expression level of flanking genes. Based on the affinities to bind to promoters and enhancers we detected 5 distinct clusters of TFs that form the core of the HOT loci. We report an abundance of HOT loci in the human genome and a commitment of 51% of all TF ChIP-seq binding events to HOT locus formation thus challenging the classical model of enhancer activity and propose a model of HOT locus formation based on the existence of large transcriptional condensates.

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.

Mapping Chromatin Occupancy of Ppp1r1b-lncRNA Genome-Wide Using Chromatin Isolation by RNA Purification (ChIRP)-seq

Long non-coding RNA (lncRNA) mediated transcriptional regulation is increasingly recognized as an important gene regulatory mechanism during development and disease. LncRNAs are emerging as critical regulators of chromatin state; yet the nature and the extent of their interactions with chromatin remain to be fully revealed. We have previously identified Ppp1r1b-lncRNA as an essential epigenetic regulator of myogenic differentiation in cardiac and skeletal myocytes in mice and humans. We further demonstrated that Ppp1r1b-lncRNA function is mediated by the interaction with the chromatin-modifying complex polycomb repressive complex 2 (PRC2) at the promoter of myogenic differentiation transcription factors, TBX5 and MyoD1. Herein, we employed unbiased chromatin isolation by RNA purification (ChIRP) and high throughput sequencing to map the repertoire of Ppp1r1b-lncRNA chromatin occupancy genome-wide in the mouse muscle myoblast cell line. We uncovered a total of 99732 true peaks corresponding to Ppp1r1b-lncRNA binding sites at high confidence (p-value < 1E-5) and enrichment score ≥ 10). The Ppp1r1b-lncRNA-binding sites averaged 558 bp in length and were distributed widely within the coding and non-coding regions of the genome. Approximately 46% of these true peaks were mapped to gene elements, of which 1180 were mapped to experimentally validated promoter sequences. Importantly, the promoter-mapped binding sites were enriched in myogenic transcription factors and heart development while exhibiting focal interactions with known motifs of proximal promoters and transcription initiation by RNA Pol-II, including TATA-box, transcription initiator motif, CCAAT-box, and GC-box, supporting Ppp1r1b-lncRNA role in transcription initiation of myogenic regulators. Remarkably, nearly 40% of Ppp1r1b-lncRNA-binding sites mapped to gene introns were enriched with the Homeobox family of transcription factors and exhibited TA-rich motif sequences, suggesting potential motif-specific Ppp1r1b-lncRNA-bound introns. Lastly, more than 136521 enhancer sequences were detected in Ppp1r1b-lncRNA-occupancy sites at high confidence. Among these enhancers, 3390 (12%) exhibited cell type/tissue-specific enrichment in fetal heart and muscles. Together, our findings provide further insights into the genome-wide Ppp1r1b-lncRNA: Chromatin interactome that may dictate its function in myogenic differentiation and potentially other cellular and biological processes.

Single-cell analysis of chromatin accessibility in the adult mouse brain

Recent advances in single-cell technologies have led to the discovery of thousands of brain cell types; however, our understanding of the gene regulatory programs in these cell types is far from complete1-4. Here we report a comprehensive atlas of candidate cis-regulatory DNA elements (cCREs) in the adult mouse brain, generated by analysing chromatin accessibility in 2.3 million individual brain cells from 117 anatomical dissections. The atlas includes approximately 1 million cCREs and their chromatin accessibility across 1,482 distinct brain cell populations, adding over 446,000 cCREs to the most recent such annotation in the mouse genome. The mouse brain cCREs are moderately conserved in the human brain. The mouse-specific cCREs-specifically, those identified from a subset of cortical excitatory neurons-are strongly enriched for transposable elements, suggesting a potential role for transposable elements in the emergence of new regulatory programs and neuronal diversity. Finally, we infer the gene regulatory networks in over 260 subclasses of mouse brain cells and develop deep-learning models to predict the activities of gene regulatory elements in different brain cell types from the DNA sequence alone. Our results provide a resource for the analysis of cell-type-specific gene regulation programs in both mouse and human brains.

Pleiotropic influence of DNA methylation QTLs on physiological and ageing traits

DNA methylation is influenced by genetic and non-genetic factors. Here, we chart quantitative trait loci (QTLs) that modulate levels of methylation at highly conserved CpGs using liver methylome data from mouse strains belonging to the BXD family. A regulatory hotspot on chromosome 5 had the highest density of trans-acting methylation QTLs (trans-meQTLs) associated with multiple distant CpGs. We refer to this locus as meQTL.5a. Trans-modulated CpGs showed age-dependent changes and were enriched in developmental genes, including several members of the MODY pathway (maturity onset diabetes of the young). The joint modulation by genotype and ageing resulted in a more 'aged methylome' for BXD strains that inherited the DBA/2J parental allele at meQTL.5a. Further, several gene expression traits, body weight, and lipid levels mapped to meQTL.5a, and there was a modest linkage with lifespan. DNA binding motif and protein-protein interaction enrichment analyses identified the hepatic nuclear factor, Hnf1a (MODY3 gene in humans), as a strong candidate. The pleiotropic effects of meQTL.5a could contribute to variations in body size and metabolic traits, and influence CpG methylation and epigenetic ageing that could have an impact on lifespan.

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.

Histone deacetylation and cytosine methylation compartmentalize heterochromatic regions in the genome organization of Neurospora crassa

Chromosomes must correctly fold in eukaryotic nuclei for proper genome function. Eukaryotic organisms hierarchically organize their genomes, including in the fungus Neurospora crassa, where chromatin fiber loops compact into Topologically Associated Domain-like structures formed by heterochromatic region aggregation. However, insufficient data exist on how histone posttranslational modifications (PTMs), including acetylation, affect genome organization. In Neurospora, the HCHC complex [composed of the proteins HDA-1, CDP-2 (Chromodomain Protein-2), Heterochromatin Protein-1, and CHAP (CDP-2 and HDA-1 Associated Protein)] deacetylates heterochromatic nucleosomes, as loss of individual HCHC members increases centromeric acetylation, and alters the methylation of cytosines in DNA. Here, we assess whether the HCHC complex affects genome organization by performing Hi-C in strains deleted of the cdp-2 or chap genes. CDP-2 loss increases intra- and interchromosomal heterochromatic region interactions, while loss of CHAP decreases heterochromatic region compaction. Individual HCHC mutants exhibit different patterns of histone PTMs genome-wide, as CDP-2 deletion increases heterochromatic H4K16 acetylation, yet smaller heterochromatic regions lose H3K9 trimethylation and gain interheterochromatic region interactions; CHAP loss produces minimal acetylation changes but increases heterochromatic H3K9me3 enrichment. Loss of both CDP-2 and the DIM-2 DNA methyltransferase causes extensive genome disorder as heterochromatic-euchromatic contacts increase despite additional H3K9me3 enrichment. Our results highlight how the increased cytosine methylation in HCHC mutants ensures genome compartmentalization when heterochromatic regions become hyperacetylated without HDAC activity.

Mapping Chromatin Occupancy of Ppp1r1b-lncRNA Genome-Wide Using Chromatin Isolation by RNA Purification (ChIRP)-seq

Long non-coding RNA (lncRNA) mediated transcriptional regulation is increasingly recognized as an important gene regulatory mechanism during development and disease. LncRNAs are emerging as critical regulators of chromatin state; yet the nature and the extent of their interactions with chromatin remain to be fully revealed. We have previously identified Ppp1r1b-lncRNA as an essential epigenetic regulator of myogenic differentiation in cardiac and skeletal myocytes in mice and humans. We further demonstrated that Ppp1r1b-lncRNA function is mediated by the interaction with the chromatin-modifying complex polycomb repressive complex 2 (PRC2) at the promoter of myogenic differentiation transcription factors, TBX5 and MyoD1. Herein, we employed an unbiased chromatin isolation by RNA purification (ChIRP) and high throughput sequencing to map the repertoire of Ppp1r1b-lncRNA chromatin occupancy genome-wide in the mouse muscle myoblast cell line. We uncovered a total of 99732 true peaks corresponding to Ppp1r1b-lncRNA binding sites at high confidence (P-value < 1e-5 and enrichment score ≥ 10). The Ppp1r1b-lncRNA-binding sites averaged 558 bp in length and were distributed widely within the coding and non-coding regions of the genome. Approximately 46% of these true peaks were mapped to gene elements, of which 1180 were mapped to experimentally validated promoter sequences. Importantly, the promoter-mapped binding sites were enriched in myogenic transcription factors and heart development while exhibiting focal interactions with known motifs of proximal promoters and transcription initiation by RNA polII, including TATA, transcription initiator, CCAAT-box, and GC-box, supporting Ppp1r1b-lncRNA role in transcription initiation of myogenic regulators. Remarkably, nearly 40% of Ppp1r1b-lncRNA-binding sites mapped to gene introns, were enriched with the Homeobox family of transcription factors, and exhibited TA-rich motif sequences, suggesting potential motif specific Ppp1r1b-lncRNA-bound introns. Lastly, more than 136521enhancer sequences were detected in Ppp1r1b-lncRNA-occupancy sites at high confidence. Among these enhancers,12% exhibited cell type/tissue-specific enrichment in fetal heart and muscles. Together, our findings provide further insights into the genome-wide Ppp1r1b-lncRNA: Chromatin interactome that may potentially dictate its function in myogenic differentiation and potentially other cellular and biological processes.

eQTL mapping in fetal-like pancreatic progenitor cells reveals early developmental insights into diabetes risk

The impact of genetic regulatory variation active in early pancreatic development on adult pancreatic disease and traits is not well understood. Here, we generate a panel of 107 fetal-like iPSC-derived pancreatic progenitor cells (iPSC-PPCs) from whole genome-sequenced individuals and identify 4065 genes and 4016 isoforms whose expression and/or alternative splicing are affected by regulatory variation. We integrate eQTLs identified in adult islets and whole pancreas samples, which reveal 1805 eQTL associations that are unique to the fetal-like iPSC-PPCs and 1043 eQTLs that exhibit regulatory plasticity across the fetal-like and adult pancreas tissues. Colocalization with GWAS risk loci for pancreatic diseases and traits show that some putative causal regulatory variants are active only in the fetal-like iPSC-PPCs and likely influence disease by modulating expression of disease-associated genes in early development, while others with regulatory plasticity likely exert their effects in both the fetal and adult pancreas by modulating expression of different disease genes in the two developmental stages.

Synaptic Activity Causes Minute-scale Changes in BAF Complex Composition and Function

Genes encoding subunits of the SWI/SNF or BAF ATP-dependent chromatin remodeling complex are among the most enriched for deleterious de novo mutations in intellectual disabilities and autism spectrum disorder, but the causative molecular pathways are not fully known 1,2 . Synaptic activity in neurons is critical for learning and memory and proper neural development 3 . Neural activity prompts calcium influx and transcription within minutes, facilitated in the nucleus by various transcription factors (TFs) and chromatin modifiers 4 . While BAF is required for activity-dependent developmental processes such as dendritic outgrowth 5-7 , the immediate molecular consequences of neural activity on BAF complexes and their functions are unknown. Here we mapped minute-scale biochemical consequences of neural activity, modeled by membrane depolarization of embryonic mouse primary cortical neurons, on BAF complexes. We used acute chemical perturbations of BAF ATPase activity and kinase signaling to define the activity-dependent effects on BAF complexes and activity-dependent BAF functions. Our studies found that BAF complexes change in subunit composition and are selectively phosphorylated within 10 minutes of depolarization. Increased levels of the core PBAF subunit Baf200/ Arid2 , uniquely containing an RFX-like DNA-binding domain, are concurrent with ATPase-dependent opening of chromatin at RFX/X-box motifs. Changes in BAF composition and phosphorylation lead to the regulation of chromatin accessibility for critical neurogenesis TFs. These biochemical effects are a convergent phenomenon downstream of multiple growth factor signaling pathways in mouse neurons and fibroblasts suggesting that BAF integrates signaling information from the membrane. In support of such a membrane-to-nucleus signaling cascade, we also identified a BAF-interacting kinase, Dclk2, whose inhibition attenuates BAF phosphorylation selectively. Our findings support a direct role of BAF complexes in responding to synaptic activity to regulate TF binding and transcription.

Iron Oxide Nanoparticles with and without Cobalt Functionalization Provoke Changes in the Transcription Profile via Epigenetic Modulation of Enhancer Activity

Despite the progress in the field of nanotoxicology, much about the cellular mechanisms that mediate the adverse effects of nanoparticles (NPs) and, in particular, the possible role of epigenetics in nanotoxicity, remains to be clarified. Therefore, we studied the changes occurring in the genome-wide distribution of H3K27ac, H3K4me1, H3K9me2, and H3K27me3 histone modifications and compared them with the transcriptome after exposing NIH3T3 cells to iron-based magnetic NPs (i.e., Fe2O3 and Fe2O3@Co NPs). We found that the transcription response is mainly due to changes in the genomic distribution of H3K27ac that can modulate the activity of enhancers. We propose that alteration of the epigenetic landscape is a key mechanism in defining the gene expression program changes resulting in nanotoxicity. With this approach, it is possible to construct a data set of genomic regions that could be useful for defining toxicity in a manner that is more comprehensive than what is possible with the present toxicology assays.

Mouse methylation profiles for leukocyte cell types, and estimation of leukocyte fractions in inflamed gastrointestinal DNA samples

Precise analysis of tissue DNA and RNA samples is often hampered by contaminating non-target cells whose amounts are highly variable. DNA methylation profiles are specific to cell types, and can be utilized for assessment of the fraction of such contaminating non-target cells. Here, we aimed 1) to identify methylation profiles specific to multiple types of mouse leukocytes, and 2) to estimate the fraction of leukocytes infiltrating inflamed tissues using DNA samples. First, genome-wide DNA methylation analysis was conducted for three myeloid-lineage cells and four lymphoid-lineage cells isolated by fluorescence-activated cell sorting after magnetic-activated cell sorting from leukocytes in the spleen. Clustering analysis using CpG sites within enhancers separated the three myeloid-lineage cells and four lymphoid-lineage cells while that using promoter CpG islands (TSS200CGIs) did not. Among the 266,108 CpG sites analyzed, one CpG site was specifically hypermethylated (β value ≥ 0.7) in B cells, and four, seven, 183, and 34 CpG sites were specifically hypomethylated (β value < 0.2) in CD4+ T cells, CD8+ T cells, B cells, and NK cells, respectively. Importantly, cell type-specific hypomethylated CpG sites were located at genes involved in cell type-specific biological functions. Then, marker CpG sites to estimate the leukocyte fraction in a tissue with leukocyte infiltration were selected, and an estimation algorithm was established. The fractions of infiltrating leukocytes were estimated to be 1.6-12.4% in the stomach (n = 10) with Helicobacter pylori-induced inflammation and 1.5-4.3% in the colon with dextran sulfate sodium-induced colitis (n = 4), and the fractions were highly correlated with those estimated histologically using Cd45-stained tissue sections [R = 0.811 (p = 0.004)]. These results showed that mouse methylation profiles at CpG sites within enhancers reflected leukocyte cell lineages, and the use of marker CpG sites successfully estimated the leukocyte fraction in inflamed gastric and colon tissues.

Histone deacetylation and cytosine methylation compartmentalize heterochromatic regions in the genome organization of Neurospora crassa

Chromosomes must correctly fold in eukaryotic nuclei for proper genome function. Eukaryotic organisms hierarchically organize their genomes, including in the fungus Neurospora crassa, where chromatin fiber loops compact into Topologically Associated Domain (TAD)-like structures formed by heterochromatic region aggregation. However, insufficient data exists on how histone post-translational modifications, including acetylation, affect genome organization. In Neurospora, the HCHC complex (comprised of the proteins HDA-1, CDP-2, HP1, and CHAP) deacetylates heterochromatic nucleosomes, as loss of individual HCHC members increases centromeric acetylation and alters the methylation of cytosines in DNA. Here, we assess if the HCHC complex affects genome organization by performing Hi-C in strains deleted of the cdp-2 or chap genes. CDP-2 loss increases intra- and inter-chromosomal heterochromatic region interactions, while loss of CHAP decreases heterochromatic region compaction. Individual HCHC mutants exhibit different patterns of histone post-translational modifications genome-wide: without CDP-2, heterochromatic H4K16 acetylation is increased, yet smaller heterochromatic regions lose H3K9 trimethylation and gain inter-heterochromatic region interactions; CHAP loss produces minimal acetylation changes but increases heterochromatic H3K9me3 enrichment. Loss of both CDP-2 and the DIM-2 DNA methyltransferase causes extensive genome disorder, as heterochromatic-euchromatic contacts increase despite additional H3K9me3 enrichment. Our results highlight how the increased cytosine methylation in HCHC mutants ensures genome compartmentalization when heterochromatic regions become hyperacetylated without HDAC activity.

Integrative analysis of single-cell RNA-seq and ATAC-seq reveals heterogeneity of induced pluripotent stem cell-derived hepatic organoids

To gain deeper insights into transcriptomes and epigenomes of organoids, liver organoids from two states (expandable and more differentiated) were subjected to single-cell RNA-seq (scRNA-seq) and single-cell ATAC-seq (scATAC-seq) analyses. Mitochondrial gene expression was higher in differentiated than in non-differentiated hepatocytes, with ATAC-seq peaks increasing near the mitochondrial control region. Differentiation of liver organoids resulted in the expression of transcription factors that act as enhancers and repressors. In addition, epigenetic mechanisms regulating the expression of alpha-fetoprotein (AFP) and albumin (ALB) differed in liver organoids and adult liver. Knockdown of PDX1, an essential transcription factor for pancreas development, led to the hepatic maturation of liver organoids through regulation of AFP and ALB expression. This integrative analysis of the transcriptomes and epigenomes of liver organoids at the single-cell level may contribute to a better understanding of the regulatory networks during liver development and the further development of mature in vitro human liver models.

CHARGE syndrome-associated CHD7 acts at ISL1-regulated enhancers to modulate second heart field gene expression

AIMS

Haploinsufficiency of the chromo-domain protein CHD7 underlies most cases of CHARGE syndrome, a multisystem birth defect including congenital heart malformation. Context specific roles for CHD7 in various stem, progenitor, and differentiated cell lineages have been reported. Previously, we showed severe defects when Chd7 is absent from cardiopharyngeal mesoderm (CPM). Here, we investigate altered gene expression in the CPM and identify specific CHD7-bound target genes with known roles in the morphogenesis of affected structures.

METHODS AND RESULTS

We generated conditional KO of Chd7 in CPM and analysed cardiac progenitor cells using transcriptomic and epigenomic analyses, in vivo expression analysis, and bioinformatic comparisons with existing datasets. We show CHD7 is required for correct expression of several genes established as major players in cardiac development, especially within the second heart field (SHF). We identified CHD7 binding sites in cardiac progenitor cells and found strong association with histone marks suggestive of dynamically regulated enhancers during the mesodermal to cardiac progenitor transition of mESC differentiation. Moreover, CHD7 shares a subset of its target sites with ISL1, a pioneer transcription factor in the cardiogenic gene regulatory network, including one enhancer modulating Fgf10 expression in SHF progenitor cells vs. differentiating cardiomyocytes.

CONCLUSION

We show that CHD7 interacts with ISL1, binds ISL1-regulated cardiac enhancers, and modulates gene expression across the mesodermal heart fields during cardiac morphogenesis.

ATRX guards against aberrant differentiation in mesenchymal progenitor cells

Alterations in the tumor suppressor ATRX are recurrently observed in several cancer types including sarcomas, which are mesenchymal neoplasms. ATRX has multiple epigenetic functions including heterochromatin formation and maintenance and regulation of transcription through modulation of chromatin accessibility. Here, we show in murine mesenchymal progenitor cells (MPCs) that Atrx deficiency aberrantly activated mesenchymal differentiation programs. This includes adipogenic pathways where ATRX loss induced expression of adipogenic transcription factors (Pparγ and Cebpα) and enhanced adipogenic differentiation in response to differentiation stimuli. These changes are linked to loss of heterochromatin near mesenchymal lineage genes together with increased chromatin accessibility and gains of active chromatin marks at putative enhancer elements and promoters. Finally, we observed depletion of H3K9me3 at transposable elements, which are derepressed including near mesenchymal genes where they could serve as regulatory elements. Our results demonstrate that ATRX functions to buffer against differentiation in mesenchymal progenitor cells, which has implications for understanding ATRX loss of function in sarcomas.

Neuron-specific chromatin disruption at CpG islands and aging-related regions in Kabuki syndrome mice

Many Mendelian developmental disorders caused by coding variants in epigenetic regulators have now been discovered. Epigenetic regulators are broadly expressed, and each of these disorders typically exhibits phenotypic manifestations from many different organ systems. An open question is whether the chromatin disruption - the root of the pathogenesis - is similar in the different disease-relevant cell types. This is possible in principle, since all these cell-types are subject to effects from the same causative gene, that has the same kind of function (e.g. methylates histones) and is disrupted by the same germline variant. We focus on mouse models for Kabuki syndrome types 1 and 2, and find that the chromatin accessibility abnormalities in neurons are mostly distinct from those in B or T cells. This is not because the neuronal abnormalities occur at regulatory elements that are only active in neurons. Neurons, but not B or T cells, show preferential chromatin disruption at CpG islands and at regulatory elements linked to aging. A sensitive analysis reveals that the regions disrupted in B/T cells do exhibit chromatin accessibility changes in neurons, but these are very subtle and of uncertain functional significance. Finally, we are able to identify a small set of regulatory elements disrupted in all three cell types. Our findings reveal the cellular-context-specific effect of variants in epigenetic regulators, and suggest that blood-derived "episignatures" may not be well-suited for understanding the mechanistic basis of neurodevelopment in Mendelian disorders of the epigenetic machinery.