Protein restriction during pregnancy alters Cdkn1c silencing, dopamine circuitry and offspring behaviour without changing expression of key neuronal marker genes.

We tracked the consequences of in utero protein restriction in mice throughout their development and life course using a luciferase-based allelic reporter of imprinted Cdkn1c. Exposure to gestational low-protein diet (LPD) results in the inappropriate expression of paternally inherited Cdkn1c in the brains of embryonic and juvenile mice. These animals were characterised by a developmental delay in motor skills, and by behavioural alterations indicative of reduced anxiety. Exposure to LPD in utero resulted in significantly more tyrosine hydroxylase positive (dopaminergic) neurons in the midbrain of adult offspring as compared to age-matched, control-diet equivalents. Positron emission tomography (PET) imaging revealed an increase in striatal dopamine synthesis capacity in LPD-exposed offspring, where elevated levels of dopamine correlated with an enhanced sensitivity to cocaine. These data highlight a profound sensitivity of the developing epigenome to gestational protein restriction. Our data also suggest that loss of Cdkn1c imprinting and p57KIP2 upregulation alters the cellular composition of the developing midbrain, compromises dopamine circuitry, and thereby provokes behavioural abnormalities in early postnatal life. Molecular analyses showed that despite this phenotype, exposure to LPD solely during pregnancy did not significantly change the expression of key neuronal- or dopamine-associated marker genes in adult offspring.

Bioluminescence imaging of Cyp1a1-luciferase reporter mice demonstrates prolonged activation of the aryl hydrocarbon receptor in the lung.

Aryl hydrocarbon receptor (AHR) signalling integrates biological processes that sense and respond to environmental, dietary, and metabolic challenges to ensure tissue homeostasis. AHR is a transcription factor that is inactive in the cytosol but upon encounter with ligand translocates to the nucleus and drives the expression of AHR targets, including genes of the cytochrome P4501 family of enzymes such as Cyp1a1. To dynamically visualise AHR activity in vivo, we generated reporter mice in which firefly luciferase (Fluc) was non-disruptively targeted into the endogenous Cyp1a1 locus. Exposure of these animals to FICZ, 3-MC or to dietary I3C induced strong bioluminescence signal and Cyp1a1 expression in many organs including liver, lung and intestine. Longitudinal studies revealed that AHR activity was surprisingly long-lived in the lung, with sustained Cyp1a1 expression evident in discrete populations of cells including columnar epithelia around bronchioles. Our data link diet to lung physiology and also reveal the power of bespoke Cyp1a1-Fluc reporters to longitudinally monitor AHR activity in vivo.

A deep learning method for replicate-based analysis of chromosome conformation contacts using Siamese neural networks.

The organisation of the genome in nuclear space is an important frontier of biology. Chromosome conformation capture methods such as Hi-C and Micro-C produce genome-wide chromatin contact maps that provide rich data containing quantitative and qualitative information about genome architecture. Most conventional approaches to genome-wide chromosome conformation capture data are limited to the analysis of pre-defined features, and may therefore miss important biological information. One constraint is that biologically important features can be masked by high levels of technical noise in the data. Here we introduce a replicate-based method for deep learning from chromatin conformation contact maps. Using a Siamese network configuration our approach learns to distinguish technical noise from biological variation and outperforms image similarity metrics across a range of biological systems. The features extracted from Hi-C maps after perturbation of cohesin and CTCF reflect the distinct biological functions of cohesin and CTCF in the formation of domains and boundaries, respectively. The learnt distance metrics are biologically meaningful, as they mirror the density of cohesin and CTCF binding. These properties make our method a powerful tool for the exploration of chromosome conformation capture data, such as Hi-C capture Hi-C, and Micro-C.

Drug-induced loss of imprinting revealed using bioluminescent reporters of Cdkn1c.

Genomic imprinting is an epigenetically mediated mechanism that regulates allelic expression of genes based upon parent-of-origin and provides a paradigm for studying epigenetic silencing and release. Here, bioluminescent reporters for the maternally-expressed imprinted gene Cdkn1c are used to examine the capacity of chromatin-modifying drugs to reverse paternal Cdkn1c silencing. Exposure of reporter mouse embryonic stem cells (mESCs) to 5-Azacytidine, HDAC inhibitors, BET inhibitors or GSK-J4 (KDM6A/B inhibitor) relieved repression of paternal Cdkn1c, either selectively or by inducing biallelic effects. Treatment of reporter fibroblasts with HDAC inhibitors or GSK-J4 resulted in similar paternal Cdkn1c activation, whereas BET inhibitor-induced loss of imprinting was specific to mESCs. Changes in allelic expression were generally not sustained in dividing cultures upon drug removal, indicating that the underlying epigenetic memory of silencing was maintained. In contrast, Cdkn1c de-repression by GSK-J4 was retained in both mESCs and fibroblasts following inhibitor removal, although this impact may be linked to cellular stress and DNA damage. Taken together, these data introduce bioluminescent reporter cells as tools for studying epigenetic silencing and disruption, and demonstrate that Cdkn1c imprinting requires distinct and cell-type specific chromatin features and modifying enzymes to enact and propagate a memory of silencing.

Loss of H3K9 trimethylation alters chromosome compaction and transcription factor retention during mitosis.

Recent studies have shown that repressive chromatin machinery, including DNA methyltransferases and polycomb repressor complexes, binds to chromosomes throughout mitosis and their depletion results in increased chromosome size. In the present study, we show that enzymes that catalyze H3K9 methylation, such as Suv39h1, Suv39h2, G9a and Glp, are also retained on mitotic chromosomes. Surprisingly, however, mutants lacking histone 3 lysine 9 trimethylation (H3K9me3) have unusually small and compact mitotic chromosomes associated with increased histone H3 phospho Ser10 (H3S10ph) and H3K27me3 levels. Chromosome size and centromere compaction in these mutants were rescued by providing exogenous first protein lysine methyltransferase Suv39h1 or inhibiting Ezh2 activity. Quantitative proteomic comparisons of native mitotic chromosomes isolated from wild-type versus Suv39h1/Suv39h2 double-null mouse embryonic stem cells revealed that H3K9me3 was essential for the efficient retention of bookmarking factors such as Esrrb. These results highlight an unexpected role for repressive heterochromatin domains in preserving transcription factor binding through mitosis and underscore the importance of H3K9me3 for sustaining chromosome architecture and epigenetic memory during cell division.

Endogenous bioluminescent reporters reveal a sustained increase in utrophin gene expression upon EZH2 and ERK1/2 inhibition.

Duchenne muscular dystrophy (DMD) is an X-linked disorder caused by loss of function mutations in the dystrophin gene (Dmd), resulting in progressive muscle weakening. Here we modelled the longitudinal expression of endogenous Dmd, and its paralogue Utrn, in mice and in myoblasts by generating bespoke bioluminescent gene reporters. As utrophin can partially compensate for Dmd-deficiency, these reporters were used as tools to ask whether chromatin-modifying drugs can enhance Utrn expression in developing muscle. Myoblasts treated with different PRC2 inhibitors showed significant increases in Utrn transcripts and bioluminescent signals, and these responses were independently verified by conditional Ezh2 deletion. Inhibition of ERK1/2 signalling provoked an additional increase in Utrn expression that was also seen in Dmd-mutant cells, and maintained as myoblasts differentiate. These data reveal PRC2 and ERK1/2 to be negative regulators of Utrn expression and provide specialised molecular imaging tools to monitor utrophin expression as a therapeutic strategy for DMD.

Chromatin jets define the properties of cohesin-driven in vivo loop extrusion.

Complex genomes show intricate organization in three-dimensional (3D) nuclear space. Current models posit that cohesin extrudes loops to form self-interacting domains delimited by the DNA binding protein CTCF. Here, we describe and quantitatively characterize cohesin-propelled, jet-like chromatin contacts as landmarks of loop extrusion in quiescent mammalian lymphocytes. Experimental observations and polymer simulations indicate that narrow origins of loop extrusion favor jet formation. Unless constrained by CTCF, jets propagate symmetrically for 1-2 Mb, providing an estimate for the range of in vivo loop extrusion. Asymmetric CTCF binding deflects the angle of jet propagation as experimental evidence that cohesin-mediated loop extrusion can switch from bi- to unidirectional and is controlled independently in both directions. These data offer new insights into the physiological behavior of in vivo cohesin-mediated loop extrusion and further our understanding of the principles that underlie genome organization.

Cohesin couples transcriptional bursting probabilities of inducible enhancers and promoters.

Innate immune responses rely on inducible gene expression programmes which, in contrast to steady-state transcription, are highly dependent on cohesin. Here we address transcriptional parameters underlying this cohesin-dependence by single-molecule RNA-FISH and single-cell RNA-sequencing. We show that inducible innate immune genes are regulated predominantly by an increase in the probability of active transcription, and that probabilities of enhancer and promoter transcription are coordinated. Cohesin has no major impact on the fraction of transcribed inducible enhancers, or the number of mature mRNAs produced per transcribing cell. Cohesin is, however, required for coupling the probabilities of enhancer and promoter transcription. Enhancer-promoter coupling may not be explained by spatial proximity alone, and at the model locus Il12b can be disrupted by selective inhibition of the cohesinopathy-associated BET bromodomain BD2. Our data identify discrete steps in enhancer-mediated inducible gene expression that differ in cohesin-dependence, and suggest that cohesin and BD2 may act on shared pathways.

Learning representations of chromatin contacts using a recurrent neural network identifies genomic drivers of conformation.

Despite the availability of chromatin conformation capture experiments, discerning the relationship between the 1D genome and 3D conformation remains a challenge, which limits our understanding of their affect on gene expression and disease. We propose Hi-C-LSTM, a method that produces low-dimensional latent representations that summarize intra-chromosomal Hi-C contacts via a recurrent long short-term memory neural network model. We find that these representations contain all the information needed to recreate the observed Hi-C matrix with high accuracy, outperforming existing methods. These representations enable the identification of a variety of conformation-defining genomic elements, including nuclear compartments and conformation-related transcription factors. They furthermore enable in-silico perturbation experiments that measure the influence of cis-regulatory elements on conformation.

Epigenetic changes induced by in utero dietary challenge result in phenotypic variability in successive generations of mice.

Transmission of epigenetic information between generations occurs in nematodes, flies and plants, mediated by specialised small RNA pathways, modified histones and DNA methylation. Similar processes in mammals can also affect phenotype through intergenerational or trans-generational mechanisms. Here we generate a luciferase knock-in reporter mouse for the imprinted Dlk1 locus to visualise and track epigenetic fidelity across generations. Exposure to high-fat diet in pregnancy provokes sustained re-expression of the normally silent maternal Dlk1 in offspring (loss of imprinting) and increased DNA methylation at the somatic differentially methylated region (sDMR). In the next generation heterogeneous Dlk1 mis-expression is seen exclusively among animals born to F1-exposed females. Oocytes from these females show altered gene and microRNA expression without changes in DNA methylation, and correct imprinting is restored in subsequent generations. Our results illustrate how diet impacts the foetal epigenome, disturbing canonical and non-canonical imprinting mechanisms to modulate the properties of successive generations of offspring.

Cohesin-dependence of neuronal gene expression relates to chromatin loop length.

Cohesin and CTCF are major drivers of 3D genome organization, but their role in neurons is still emerging. Here, we show a prominent role for cohesin in the expression of genes that facilitate neuronal maturation and homeostasis. Unexpectedly, we observed two major classes of activity-regulated genes with distinct reliance on cohesin in mouse primary cortical neurons. Immediate early genes (IEGs) remained fully inducible by KCl and BDNF, and short-range enhancer-promoter contacts at the IEGs Fos formed robustly in the absence of cohesin. In contrast, cohesin was required for full expression of a subset of secondary response genes characterized by long-range chromatin contacts. Cohesin-dependence of constitutive neuronal genes with key functions in synaptic transmission and neurotransmitter signaling also scaled with chromatin loop length. Our data demonstrate that key genes required for the maturation and activation of primary cortical neurons depend on cohesin for their full expression, and that the degree to which these genes rely on cohesin scales with the genomic distance traversed by their chromatin contacts.

Neuronal genes deregulated in Cornelia de Lange Syndrome respond to removal and re-expression of cohesin.

Cornelia de Lange Syndrome (CdLS) is a human developmental disorder caused by mutations that compromise the function of cohesin, a major regulator of 3D genome organization. Cognitive impairment is a universal and as yet unexplained feature of CdLS. We characterize the transcriptional profile of cortical neurons from CdLS patients and find deregulation of hundreds of genes enriched for neuronal functions related to synaptic transmission, signalling processes, learning and behaviour. Inducible proteolytic cleavage of cohesin disrupts 3D genome organization and transcriptional control in post-mitotic cortical mouse neurons, demonstrating that cohesin is continuously required for neuronal gene expression. The genes affected by acute depletion of cohesin belong to similar gene ontology classes and show significant numerical overlap with genes deregulated in CdLS. Interestingly, reconstitution of cohesin function largely rescues altered gene expression, including the expression of genes deregulated in CdLS.

The order and logic of CD4 versus CD8 lineage choice and differentiation in mouse thymus.

CD4 and CD8 mark helper and cytotoxic T cell lineages, respectively, and serve as coreceptors for MHC-restricted TCR recognition. How coreceptor expression is matched with TCR specificity is central to understanding CD4/CD8 lineage choice, but visualising coreceptor gene activity in individual selection intermediates has been technically challenging. It therefore remains unclear whether the sequence of coreceptor gene expression in selection intermediates follows a stereotypic pattern, or is responsive to signaling. Here we use single cell RNA sequencing (scRNA-seq) to classify mouse thymocyte selection intermediates by coreceptor gene expression. In the unperturbed thymus, Cd4+Cd8a- selection intermediates appear before Cd4-Cd8a+ selection intermediates, but the timing of these subsets is flexible according to the strength of TCR signals. Our data show that selection intermediates discriminate MHC class prior to the loss of coreceptor expression and suggest a model where signal strength informs the timing of coreceptor gene activity and ultimately CD4/CD8 lineage choice.

Identifying proteins bound to native mitotic ESC chromosomes reveals chromatin repressors are important for compaction.

Epigenetic information is transmitted from mother to daughter cells through mitosis. Here, to identify factors that might play a role in conveying epigenetic memory through cell division, we report on the isolation of unfixed, native chromosomes from metaphase-arrested cells using flow cytometry and perform LC-MS/MS to identify chromosome-bound proteins. A quantitative proteomic comparison between metaphase-arrested cell lysates and chromosome-sorted samples reveals a cohort of proteins that were significantly enriched on mitotic ESC chromosomes. These include pluripotency-associated transcription factors, repressive chromatin-modifiers such as PRC2 and DNA methyl-transferases, and proteins governing chromosome architecture. Deletion of PRC2, Dnmt1/3a/3b or Mecp2 in ESCs leads to an increase in the size of individual mitotic chromosomes, consistent with de-condensation. Similar results were obtained by the experimental cleavage of cohesin. Thus, we identify chromosome-bound factors in pluripotent stem cells during mitosis and reveal that PRC2, DNA methylation and Mecp2 are required to maintain chromosome compaction.

STATegra, a comprehensive multi-omics dataset of B-cell differentiation in mouse.

Multi-omics approaches use a diversity of high-throughput technologies to profile the different molecular layers of living cells. Ideally, the integration of this information should result in comprehensive systems models of cellular physiology and regulation. However, most multi-omics projects still include a limited number of molecular assays and there have been very few multi-omic studies that evaluate dynamic processes such as cellular growth, development and adaptation. Hence, we lack formal analysis methods and comprehensive multi-omics datasets that can be leveraged to develop true multi-layered models for dynamic cellular systems. Here we present the STATegra multi-omics dataset that combines measurements from up to 10 different omics technologies applied to the same biological system, namely the well-studied mouse pre-B-cell differentiation. STATegra includes high-throughput measurements of chromatin structure, gene expression, proteomics and metabolomics, and it is complemented with single-cell data. To our knowledge, the STATegra collection is the most diverse multi-omics dataset describing a dynamic biological system.

Chromatinization of Escherichia coli with archaeal histones.

Nucleosomes restrict DNA accessibility throughout eukaryotic genomes, with repercussions for replication, transcription, and other DNA-templated processes. How this globally restrictive organization emerged during evolution remains poorly understood. Here, to better understand the challenges associated with establishing globally restrictive chromatin, we express histones in a naive system that has not evolved to deal with nucleosomal structures: Escherichia coli. We find that histone proteins from the archaeon Methanothermus fervidus assemble on the E. coli chromosome in vivo and protect DNA from micrococcal nuclease digestion, allowing us to map binding footprints genome-wide. We show that higher nucleosome occupancy at promoters is associated with lower transcript levels, consistent with local repressive effects. Surprisingly, however, this sudden enforced chromatinization has only mild repercussions for growth unless cells experience topological stress. Our results suggest that histones can become established as ubiquitous chromatin proteins without interfering critically with key DNA-templated processes.

Building gene regulatory networks from scATAC-seq and scRNA-seq using Linked Self Organizing Maps.

Rapid advances in single-cell assays have outpaced methods for analysis of those data types. Different single-cell assays show extensive variation in sensitivity and signal to noise levels. In particular, scATAC-seq generates extremely sparse and noisy datasets. Existing methods developed to analyze this data require cells amenable to pseudo-time analysis or require datasets with drastically different cell-types. We describe a novel approach using self-organizing maps (SOM) to link scATAC-seq regions with scRNA-seq genes that overcomes these challenges and can generate draft regulatory networks. Our SOMatic package generates chromatin and gene expression SOMs separately and combines them using a linking function. We applied SOMatic on a mouse pre-B cell differentiation time-course using controlled Ikaros over-expression to recover gene ontology enrichments, identify motifs in genomic regions showing similar single-cell profiles, and generate a gene regulatory network that both recovers known interactions and predicts new Ikaros targets during the differentiation process. The ability of linked SOMs to detect emergent properties from multiple types of highly-dimensional genomic data with very different signal properties opens new avenues for integrative analysis of heterogeneous data.

Epigenomic signatures underpin the axonal regenerative ability of dorsal root ganglia sensory neurons.

Axonal injury results in regenerative success or failure, depending on whether the axon lies in the peripheral or the CNS, respectively. The present study addresses whether epigenetic signatures in dorsal root ganglia discriminate between regenerative and non-regenerative axonal injury. Chromatin immunoprecipitation for the histone 3 (H3) post-translational modifications H3K9ac, H3K27ac and H3K27me3; an assay for transposase-accessible chromatin; and RNA sequencing were performed in dorsal root ganglia after sciatic nerve or dorsal column axotomy. Distinct histone acetylation and chromatin accessibility signatures correlated with gene expression after peripheral, but not central, axonal injury. DNA-footprinting analyses revealed new transcriptional regulators associated with regenerative ability. Machine-learning algorithms inferred the direction of most of the gene expression changes. Neuronal conditional deletion of the chromatin remodeler CCCTC-binding factor impaired nerve regeneration, implicating chromatin organization in the regenerative competence. Altogether, the present study offers the first epigenomic map providing insight into the transcriptional response to injury and the differential regenerative ability of sensory neurons.