ANDA: an open-source tool for automated image analysis of in vitro neuronal cells

BACKGROUND

Imaging of in vitro neuronal differentiation and measurements of cell morphologies have led to novel insights into neuronal development. Live-cell imaging techniques and large datasets of images have increased the demand for automated pipelines for quantitative analysis of neuronal morphological metrics.

RESULTS

ANDA is an analysis workflow that quantifies various aspects of neuronal morphology from high-throughput live-cell imaging screens of in vitro neuronal cell types. This tool automates the analysis of neuronal cell numbers, neurite lengths and neurite attachment points. We used chicken, rat, mouse, and human in vitro models for neuronal differentiation and have demonstrated the accuracy, versatility, and efficiency of the tool.

CONCLUSIONS

ANDA is an open-source tool that is easy to use and capable of automated processing from time-course measurements of neuronal cells. The strength of this pipeline is the capability to analyse high-throughput imaging screens.

Multi-omics approach reveals dysregulated genes during hESCs neuronal differentiation exposure to paracetamol

Prenatal paracetamol exposure has been associated with neurodevelopmental outcomes in childhood. Pharmacoepigenetic studies show differences in cord blood DNA methylation between unexposed and paracetamol-exposed neonates, however, causality and impact of long-term prenatal paracetamol exposure on brain development remain unclear. Using a multi-omics approach, we investigated the effects of paracetamol on an in vitro model of early human neurodevelopment. We exposed human embryonic stem cells undergoing neuronal differentiation with paracetamol concentrations corresponding to maternal therapeutic doses. Single-cell RNA-seq and ATAC-seq integration identified paracetamol-induced chromatin opening changes linked to gene expression. Differentially methylated and/or expressed genes were involved in neurotransmission and cell fate determination trajectories. Some genes involved in neuronal injury and development-specific pathways, such as KCNE3, overlapped with differentially methylated genes previously identified in cord blood associated with prenatal paracetamol exposure. Our data suggest that paracetamol may play a causal role in impaired neurodevelopment.

MYB regulates the SUMO protease SENP1 and its novel interaction partner UXT, modulating MYB target genes and the SUMO landscape

SUMOylation is a post-translational modification frequently found on nuclear proteins, including transcription factors (TFs) and coactivators. By controlling the activity of several TFs, SUMOylation may have far-reaching effects. MYB is an example of a developmental TF subjected to SUMO-mediated regulation, through both SUMO conjugation and SUMO binding. How SUMO affects MYB target genes is unknown. Here, we explored the global effect of reduced SUMOylation of MYB on its downstream gene programs. RNA-Seq in K562 cells after MYB knockdown and rescue with mutants having an altered SUMO status revealed a number of differentially regulated genes and distinct gene ontology term enrichments. Clearly, the SUMO status of MYB both quantitatively and qualitatively affects its regulome. The transcriptome data further revealed that MYB upregulates the SUMO protease SENP1, a key enzyme that removes SUMO conjugation from SUMOylated proteins. Given this role of SENP1 in the MYB regulome, we expanded the analysis, mapped interaction partners of SENP1, and identified UXT as a novel player affecting the SUMO system by acting as a repressor of SENP1. MYB inhibits the expression of UXT suggesting that MYB is able not only to control a specific gene program directly but also indirectly by affecting the SUMO landscape through SENP1 and UXT. These findings suggest an autoactivation loop whereby MYB, through enhancing SENP1 and reducing UXT, is itself being activated by a reduced level of repressive SUMOylation. We propose that overexpressed MYB, seen in multiple cancers, may drive this autoactivation loop and contribute to oncogenic activation of MYB.

ANDA: An open-source tool for automated image analysis of neuronal differentiation

Background

Imaging of in vitro neuronal differentiation and measurements of cell morphologies has led to novel insights into neuronal development. Live-cell imaging techniques and large datasets of images has increased the demand for automated pipelines for quantitative analysis of neuronal morphological metrics.

Results

We present ANDA, an analysis workflow for quantification of various aspects of neuronal morphology from high-throughput live-cell imaging screens. This tool automates the analysis of neuronal cell numbers, neurite lengths and neurite attachment points. We used rat, chicken and human in vitro models for neuronal differentiation and have demonstrated the accuracy, versatility, and efficiency of the tool.

Conclusions

ANDA is an open-source tool that is easy to use and capable of automated processing from time-course measurements of neuronal cells. The strength of this pipeline is the capability to analyse high-throughput imaging screens.

A multi-omics approach to visualize early neuronal differentiation from hESCs in 4D

Neuronal differentiation of pluripotent stem cells is an established method to study physiology, disease, and medication safety. However, the sequence of events in human neuronal differentiation and the ability of in vitro models to recapitulate early brain development are poorly understood. We developed a protocol optimized for the study of early human brain development and neuropharmacological applications. We comprehensively characterized gene expression and epigenetic profiles at four timepoints, because the cells differentiate from embryonic stem cells towards a heterogeneous population of progenitors, immature and mature neurons bearing telencephalic signatures. A multi-omics roadmap of neuronal differentiation, combined with searchable interactive gene analysis tools, allows for extensive exploration of early neuronal development and the effect of medications.

•

Multi-omics charting a new neuronal differentiation protocol for human ES cells

•

Single-cell analyses reveal marker genes during neuronal differentiation

•

Identified transcriptional waves similar to early human brain development

•

Searchable tools to visualize single-cell gene expression and chromatin state

Robust neuronal differentiation of human embryonic stem cells for neurotoxicology

Here, we describe a protocol for rapid neuronal differentiation from human embryonic stem cells (hESCs) toward a heterogenous population of telencephalic progenitors, immature and mature neurons, for drug-screening and early-brain differentiation studies. hESC neuronal differentiation depends on adhesion and minimal cell-passaging to avert monolayer cross-connectivity rupture. In this protocol, we detail optimized cell-seeding densities and coating conditions with high cell viability suitable for neurotoxicology and high-resolution single-cell omics studies. Daily media changes reduce compound instability and degradation for optimal screening.

For complete details on the use and execution of this protocol, please refer to Samara et al. (2022).

•

Day-by-day visual aids to monitor morphological changes during differentiation

•

Optimized cell numbers and passages for neurotoxicology studies

•

Web tools available for single-cell omics data

•

Ideal for early brain development exploration from human embryonic stem cells

Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Paracetamol perturbs neuronal arborization and disrupts the cytoskeletal proteins SPTBN1 and TUBB3 in both human and chicken in vitro models

Epidemiological studies have linked long-term/high-dose usage of paracetamol (N-acetyl-para-aminophenol, APAP) during pregnancy to adverse neuropsychiatric outcomes, primarily attention-deficit hyperactive disorder (ADHD), in the offspring. Though variable, ADHD has been associated with phenotypic alterations characterized by reductions in grey matter densities and aberrations in structural connectivity, effects which are thought to originate in neurodevelopment. We used embryonic chicken cerebellar granule neurons (CGNs) and neuronally differentiating human NTERA2 cells (NT2Ns) to investigate the in vitro effects of APAP on cell viability, migration, neuritogenesis, and the intracellular levels of various proteins involved in neurodevelopment as well as in the maintenance of the structure and function of neurites. Exposure to APAP ranging from 100 to 1600 μM yielded concentration- and time-dependent reductions in cell viability and levels of neurite arborization, as well as reductions in the levels of the cytoskeletal protein β2-spectrin, with the highest APAP concentration resulting in between 50 and 75% reductions in the aforementioned metrics over the course of 72 h. Exposure to APAP also reduced migration in the NT2Ns but not CGNs. Moreover, we found concentration- and time-dependent increases in punctate aggregation of the cytoskeletal protein β3-tubulin following exposure to APAP in both cell model systems, with the highest APAP concentration approximately doubling the number of aggregates over 72-120 h. Our findings demonstrate that APAP negatively perturbs neurite arborization degree, with concurrent reductions in the protein levels of β2-spectrin and disruption of the integrity of β3-tubulin, both proteins of which play important roles in neuronal structure and function.

Pioneer transcription factors are associated with the modulation of DNA methylation patterns across cancers

Methylation of cytosines on DNA is a prominent modification associated with gene expression regulation. Aberrant DNA methylation patterns have recurrently been linked to dysregulation of the regulatory program in cancer cells. To shed light on the underlying molecular mechanism driving this process, we hypothesised that aberrant methylation patterns could be controlled by the binding of specific transcription factors (TFs) across cancer types. By combining DNA methylation arrays and gene expression data with TF binding sites (TFBSs), we explored the interplay between TF binding and DNA methylation in 19 cancer types. We performed emQTL (expression–methylation quantitative trait loci) analyses independently in each cancer type and identified 13 TFs whose expression levels are correlated with local DNA methylation patterns around their binding sites in at least 2 cancer types. The 13 TFs are mainly associated with local demethylation and are enriched for pioneer function, suggesting a specific role for these TFs in modulating chromatin structure and transcription in cancer patients. Furthermore, we confirmed that de novo methylation is precluded across cancers at CpGs lying in genomic regions enriched for TF binding signatures associated with SP1, CTCF, NRF1, GABPA, KLF9, and/or YY1. The modulation of DNA methylation associated with TF binding was observed at cis-regulatory regions controlling immune- and cancer-associated pathways, corroborating that the emQTL signals were derived from both cancer and tumor-infiltrating cells. As a case example, we experimentally confirmed that FOXA1 knock-down is associated with higher methylation in regions bound by FOXA1 in breast cancer MCF-7 cells. Finally, we reported physical interactions between FOXA1 with TET1 and TET2 both in an in vitro setup and in vivo at physiological levels in MCF-7 cells, adding further support for FOXA1 attracting TET1 and TET2 to induce local demethylation in cancer cells.

ShinyArchR.UiO: user-friendly,integrative and open-source tool for visualization of single-cell ATAC-seq data using ArchR

MOTIVATION

Mapping of chromatin accessibility landscapes in single-cells and the integration with gene expression enables a better understanding of gene regulatory mechanisms defining cell identities and cell-fate determination in development and disease. Generally, raw data generated from single-cell Assay for Transposase-Accessible Chromatin sequencing (scATAC-seq) are deposited in repositories that are generally inaccessible due to lack of in-depth knowledge of computational programming.

RESULTS

We have developed ShinyArchR.UiO, an R-based shiny app, that facilitates scATAC-seq data accessibility and visualization in a user-friendly, interactive and open-source web interface. ShinyArchR.UiO is an application that can streamline collaborative efforts for interpretation of massive chromatin accessibility datasets and allow for open access data sharing for wider audiences.

AVAILABILITY AND IMPLEMENTATION

https://Github.com/EskelandLab/ShinyArchRUiO and a demo server with a hematopoietic tutorial dataset https://cancell.medisin.uio.no/ShinyArchR.UiO.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Chromatin occupancy and target genes of the haematopoietic master transcription factor MYB

The transcription factor MYB is a master regulator in haematopoietic progenitor cells and a pioneer factor affecting differentiation and proliferation of these cells. Leukaemic transformation may be promoted by high MYB levels. Despite much accumulated molecular knowledge of MYB, we still lack a comprehensive understanding of its target genes and its chromatin action. In the present work, we performed a ChIP-seq analysis of MYB in K562 cells accompanied by detailed bioinformatics analyses. We found that MYB occupies both promoters and enhancers. Five clusters (C1-C5) were found when we classified MYB peaks according to epigenetic profiles. C1 was enriched for promoters and C2 dominated by enhancers. C2-linked genes were connected to hematopoietic specific functions and had GATA factor motifs as second in frequency. C1 had in addition to MYB-motifs a significant frequency of ETS-related motifs. Combining ChIP-seq data with RNA-seq data allowed us to identify direct MYB target genes. We also compared ChIP-seq data with digital genomic footprinting. MYB is occupying nearly a third of the super-enhancers in K562. Finally, we concluded that MYB cooperates with a subset of the other highly expressed TFs in this cell line, as expected for a master regulator.

The SUMO protease SENP1 and the chromatin remodeler CHD3 interact and jointly affect chromatin accessibility and gene expression

The small ubiquitin-like modifier (SUMO) post-translationally modifies lysine residues of transcription factors and co-regulators and thereby contributes to an important layer of control of the activities of these transcriptional regulators. Likewise, deSUMOylation of these factors by the sentrin-specific proteases (SENPs) also plays a role in gene regulation, but whether SENPs functionally interact with other regulatory factors that control gene expression is unclear. In the present work, we focused on SENP1, specifically, on its role in activation of gene expression investigated through analysis of the SENP1 interactome, which revealed that SENP1 physically interacts with the chromatin remodeler chromodomain helicase DNA-binding protein 3 (CHD3). Using several additional methods, including GST pulldown and co-immunoprecipitation assays, we validated and mapped this interaction, and using CRISPR-Cas9-generated CHD3- and SENP1-KO cells (in the haploid HAP1 cell line), we investigated whether these two proteins are functionally linked in regulating chromatin remodeling and gene expression. Genome-wide ATAC-Seq analysis of the CHD3- and SENP1-KO cells revealed a large degree of overlap in differential chromatin openness between these two mutant cell lines. Moreover, motif analysis and comparison with ChIP-Seq profiles in K562 cells pointed to an association of CHD3 and SENP1 with CCCTC-binding factor (CTCF) and SUMOylated chromatin-associated factors. Lastly, genome-wide RNA-Seq also indicated that these two proteins co-regulate the expression of several genes. We propose that the functional link between chromatin remodeling by CHD3 and deSUMOylation by SENP1 uncovered here provides another level of control of gene expression.

The pioneer factor activity of c-Myb involves recruitment of p300 and induction of histone acetylation followed by acetylation-induced chromatin dissociation

Background

The concept of pioneer transcription factors is emerging as an essential part of the epigenetic regulation, taking place during cell development and differentiation. However, the precise molecular mechanism underlying pioneer factor activity remains poorly understood. We recently reported that the transcription factor c-Myb acts as a pioneer factor in haematopoiesis, and a point mutation in its DNA binding domain, D152V, is able to abrogate this function.

Results

Here, we show that specific histone modifications, including H3K27ac, prevent binding of c-Myb to histone tails, representing a novel effect of histone modifications, namely restricting binding of a pioneer factor to chromatin. Furthermore, we have taken advantage of the pioneer-defect D152V mutant to investigate mechanisms of c-Myb’s pioneer factor activity. We show that c-Myb-dependent transcriptional activation of a gene in inaccessible chromatin relies on c-Myb binding to histones, as well as on c-Myb interacting with the histone acetyltransferase p300. ChIP assays show that both wild type and the D152V mutant of c-Myb bind to a selected target gene at its promoter and enhancer, but only wild-type c-Myb causes opening and activation of the locus. Enhancement of histone acetylation restores activation of the same gene in the absence of c-Myb, suggesting that facilitating histone acetylation is a crucial part of the pioneer factor function of c-Myb.

Conclusions

We suggest a pioneer factor model in which c-Myb binds to regions of closed chromatin and then recruits histone acetyltransferases. By binding to histones, c-Myb facilitates histone acetylation, acting as a cofactor for p300 at c-Myb bound sites. The resulting H3K27ac leads to chromatin opening and detachment of c-Myb from the acetylated chromatin. We propose that the latter phenomenon, acetylation-induced chromatin dissociation, represents a mechanism for controlling the dynamics of pioneer factor binding to chromatin.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0208-y) contains supplementary material, which is available to authorized users.

A c-Myb mutant causes deregulated differentiation due to impaired histone binding and abrogated pioneer factor function

The transcription factor c-Myb is involved in early differentiation and proliferation of haematopoietic cells, where it operates as a regulator of self-renewal and multi-lineage differentiation. Deregulated c-Myb plays critical roles in leukaemias and other human cancers. Due to its role as a master regulator, we hypothesized it might function as a pioneer transcription factor. Our approach to test this was to analyse a mutant of c-Myb, D152V, previously reported to cause haematopoietic defects in mice by an unknown mechanism. Our transcriptome data from K562 cells indicates that this mutation specifically affects c-Myb's ability to regulate genes involved in differentiation, causing failure in c-Myb's ability to block differentiation. Furthermore, we see a major effect of this mutation in assays where chromatin opening is involved. We show that each repeat in the minimal DNA-binding domain of c-Myb binds to histones and that D152V disrupts histone binding of the third repeat. ATAC-seq data indicates this mutation impairs the ability of c-Myb to cause chromatin opening at specific sites. Taken together, our findings support that c-Myb acts as a pioneer factor and show that D152V impairs this function. The D152V mutant is the first mutant of a transcription factor specifically destroying pioneer factor function.

Mobilization of LINE-1 retrotransposons is restricted by Tex19.1 in mouse embryonic stem cells

Mobilization of retrotransposons to new genomic locations is a significant driver of mammalian genome evolution, but these mutagenic events can also cause genetic disorders. In humans, retrotransposon mobilization is mediated primarily by proteins encoded by LINE-1 (L1) retrotransposons, which mobilize in pluripotent cells early in development. Here we show that TEX19.1, which is induced by developmentally programmed DNA hypomethylation, can directly interact with the L1-encoded protein L1-ORF1p, stimulate its polyubiquitylation and degradation, and restrict L1 mobilization. We also show that TEX19.1 likely acts, at least in part, through promoting the activity of the E3 ubiquitin ligase UBR2 towards L1-ORF1p. Moreover, loss of Tex19.1 increases L1-ORF1p levels and L1 mobilization in pluripotent mouse embryonic stem cells, implying that Tex19.1 prevents de novo retrotransposition in the pluripotent phase of the germline cycle. These data show that post-translational regulation of L1 retrotransposons plays a key role in maintaining trans-generational genome stability in mammals.

Abstract LB-097: DNA methylation at enhancers distinguishes distinct breast cancer lineages

Breast cancers exhibit genome-wide aberrant DNA methylation patterns. To investigate how these affect the transcriptome and which changes are linked to transformation or progression we applied genome-wide expression-methylation quantitative trait loci (emQTL) analysis between 5meCpG and gene expression. We find that on a whole genome scale, DNA methylation and gene expression have remarkably and reproducibly conserved patterns of association, both in cis and in trans, which lead to identify distinct transcriptional networks and specific disease cell lineages deviating from the normal tissue. In three breast cancer cohorts (n=104, n=253 and n=330), we invariably identify enhancers whose DNA methylation tethers the binding of three key transcription factors in breast cancer in a lineage specific manner. Our emQTL analysis also identifies tumor infiltrating immune cell signatures. Using ChromHMM segmentation and ChIP-seq associated TF binding regions, we identify enhancers and TF binding sites around emQTLs-CpGs. We functionally validate the genes associated with these CpGs using knock-down by siRNA and applying GRO-seq analysis after transcriptional stimulation with estrogen. DNA methylation at TF binding regions is considered an early event during normal breast cell transformation into estrogen dependent breast cancer. Two distinct key gene regulatory networks reported here are prominently altered by DNA methylation.

Citation Format: Thomas Fleischer, Xavier Tekpli, Anthony Mathelier, Shixiong Wang, Daniel Nebdal, Hari P. Dhakal, Kristine Kleivi Sandberg, Ellen Schliting, Oslo Breast Cancer Research Consortium (OSBREAC), Anne-Lise Børresen-Dale, Elin Borgen, Bjørn Naume, Ragnhild Eskeland, Arnoldo Frigessi, Jörg Tost, Antoni Hurtado, Vessela N. Kristensen. DNA methylation at enhancers distinguishes distinct breast cancer lineages [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr LB-097. doi:10.1158/1538-7445.AM2017-LB-097

GSuite HyperBrowser: integrative analysis of dataset collections across the genome and epigenome

Background

Recent large-scale undertakings such as ENCODE and Roadmap Epigenomics have generated experimental data mapped to the human reference genome (as genomic tracks) representing a variety of functional elements across a large number of cell types. Despite the high potential value of these publicly available data for a broad variety of investigations, little attention has been given to the analytical methodology necessary for their widespread utilisation.

Findings

We here present a first principled treatment of the analysis of collections of genomic tracks. We have developed novel computational and statistical methodology to permit comparative and confirmatory analyses across multiple and disparate data sources. We delineate a set of generic questions that are useful across a broad range of investigations and discuss the implications of choosing different statistical measures and null models. Examples include contrasting analyses across different tissues or diseases. The methodology has been implemented in a comprehensive open-source software system, the GSuite HyperBrowser. To make the functionality accessible to biologists, and to facilitate reproducible analysis, we have also developed a web-based interface providing an expertly guided and customizable way of utilizing the methodology. With this system, many novel biological questions can flexibly be posed and rapidly answered.

Conclusions

Through a combination of streamlined data acquisition, interoperable representation of dataset collections, and customizable statistical analysis with guided setup and interpretation, the GSuite HyperBrowser represents a first comprehensive solution for integrative analysis of track collections across the genome and epigenome. The software is available at: https://hyperbrowser.uio.no.

In the loop: promoter-enhancer interactions and bioinformatics

Enhancer-promoter regulation is a fundamental mechanism underlying differential transcriptional regulation. Spatial chromatin organization brings remote enhancers in contact with target promoters in cis to regulate gene expression. There is considerable evidence for promoter-enhancer interactions (PEIs). In the recent years, genome-wide analyses have identified signatures and mapped novel enhancers; however, being able to precisely identify their target gene(s) requires massive biological and bioinformatics efforts. In this review, we give a short overview of the chromatin landscape and transcriptional regulation. We discuss some key concepts and problems related to chromatin interaction detection technologies, and emerging knowledge from genome-wide chromatin interaction data sets. Then, we critically review different types of bioinformatics analysis methods and tools related to representation and visualization of PEI data, raw data processing and PEI prediction. Lastly, we provide specific examples of how PEIs have been used to elucidate a functional role of non-coding single-nucleotide polymorphisms. The topic is at the forefront of epigenetic research, and by highlighting some future bioinformatics challenges in the field, this review provides a comprehensive background for future PEI studies.

Spatial genome organization: contrasting views from chromosome conformation capture and fluorescence in situ hybridization

Although important for gene regulation, most studies of genome organization use either fluorescence in situ hybridization (FISH) or chromosome conformation capture (3C) methods. FISH directly visualizes the spatial relationship of sequences but is usually applied to a few loci at a time. The frequency at which sequences are ligated together by formaldehyde cross-linking can be measured genome-wide by 3C methods, with higher frequencies thought to reflect shorter distances. FISH and 3C should therefore give the same views of genome organization, but this has not been tested extensively. We investigated the murine HoxD locus with 3C carbon copy (5C) and FISH in different developmental and activity states and in the presence or absence of epigenetic regulators. We identified situations in which the two data sets are concordant but found other conditions under which chromatin topographies extrapolated from 5C or FISH data are not compatible. We suggest that products captured by 3C do not always reflect spatial proximity, with ligation occurring between sequences located hundreds of nanometers apart, influenced by nuclear environment and chromatin composition. We conclude that results obtained at high resolution with either 3C methods or FISH alone must be interpreted with caution and that views about genome organization should be validated by independent methods.

H4K16 acetylation marks active genes and enhancers of embryonic stem cells, but does not alter chromatin compaction

Compared with histone H3, acetylation of H4 tails has not been well studied, especially in mammalian cells. Yet, H4K16 acetylation is of particular interest because of its ability to decompact nucleosomes in vitro and its involvement in dosage compensation in flies. Here we show that, surprisingly, loss of H4K16 acetylation does not alter higher-order chromatin compaction in vivo in mouse embryonic stem cells (ESCs). As well as peaks of acetylated H4K16 and KAT8 histone acetyltransferase at the transcription start sites of expressed genes, we report that acetylation of H4K16 is a new marker of active enhancers in ESCs and that some enhancers are marked by H3K4me1, KAT8, and H4K16ac, but not by acetylated H3K27 or EP300, suggesting that they are novel EP300 independent regulatory elements. Our data suggest a broad role for different histone acetylation marks and for different histone acetyltransferases in long-range gene regulation.

Lamin A/C-promoter interactions specify chromatin state-dependent transcription outcomes

The nuclear lamina is implicated in the organization of the eukaryotic nucleus. Association of nuclear lamins with the genome occurs through large chromatin domains including mostly, but not exclusively, repressed genes. How lamin interactions with regulatory elements modulate gene expression in different cellular contexts is unknown. We show here that in human adipose tissue stem cells, lamin A/C interacts with distinct spatially restricted subpromoter regions, both within and outside peripheral and intra-nuclear lamin-rich domains. These localized interactions are associated with distinct transcriptional outcomes in a manner dependent on local chromatin modifications. Down-regulation of lamin A/C leads to dissociation of lamin A/C from promoters and remodels repressive and permissive histone modifications by enhancing transcriptional permissiveness, but is not sufficient to elicit gene activation. Adipogenic differentiation resets a large number of lamin-genome associations globally and at subpromoter levels and redefines associated transcription outputs. We propose that lamin A/C acts as a modulator of local gene expression outcome through interaction with adjustable sites on promoters, and that these position-dependent transcriptional readouts may be reset upon differentiation.

Oslo Epigenetics Symposium 2012. Oslo, Norway, 8-9 November 2012

The Oslo Epigenetics Symposium 2012 held in Oslo, Norway, brought together ten speakers from several European countries and the USA for an evening public lecture and a full day of presentations on emerging topics in the field of epigenetics, gene regulation and organization of the cell nucleus.

Anterior-posterior differences in HoxD chromatin topology in limb development

A late phase of HoxD activation is crucial for the patterning and growth of distal structures across the anterior-posterior (A-P) limb axis of mammals. Polycomb complexes and chromatin compaction have been shown to regulate Hox loci along the main body axis in embryonic development, but the extent to which they have a role in limb-specific HoxD expression, an evolutionary adaptation defined by the activity of distal enhancer elements that drive expression of 5' Hoxd genes, has yet to be fully elucidated. We reveal two levels of chromatin topology that differentiate distal limb A-P HoxD activity. Using both immortalised cell lines derived from posterior and anterior regions of distal E10.5 mouse limb buds, and analysis in E10.5 dissected limb buds themselves, we show that there is a loss of polycomb-catalysed H3K27me3 histone modification and a chromatin decompaction over HoxD in the distal posterior limb compared with anterior. Moreover, we show that the global control region (GCR) long-range enhancer spatially colocalises with the 5' HoxD genomic region specifically in the distal posterior limb. This is consistent with the formation of a chromatin loop between 5' HoxD and the GCR regulatory module at the time and place of distal limb bud development when the GCR participates in initiating Hoxd gene quantitative collinearity and Hoxd13 expression. This is the first example of A-P differences in chromatin compaction and chromatin looping in the development of the mammalian secondary body axis (limb).

Identification of c-Myb Target Genes in K562 Cells Reveals a Role for c-Myb as a Master Regulator

The c-Myb transcription factor is an important regulator of hematopoietic cell development. c-Myb is expressed in immature hematopoietic cells and plays a direct role in lineage fate selection, cell cycle progression, and differentiation of myeloid as well as B- and T-lymphoid progenitor cells. As a DNA-binding transcription factor, c-Myb regulates specific gene programs through activation of target genes. Still, our understanding of these programs is incomplete. Here, we report a set of novel c-Myb target genes, identified using a combined approach: specific c-Myb knockdown by 2 different siRNAs and subsequent global expression profiling, combined with the confirmation of direct binding of c-Myb to the target promoters by ChIP assays. The combination of these 2 approaches, as well as additional validation such as cloning and testing the promoters in reporter assays, confirmed that MYADM, LMO2, GATA2, STAT5A, and IKZF1 are target genes of c-Myb. Additional studies, using chromosome conformation capture, demonstrated that c-Myb target genes may directly interact with each other, indicating that these genes may be coordinately regulated. Of the 5 novel target genes identified, 3 are transcription factors, and one is a transcriptional co-regulator, supporting a role of c-Myb as a master regulator controlling the expression of other transcriptional regulators in the hematopoietic system.

Histone H2A mono-ubiquitination is a crucial step to mediate PRC1-dependent repression of developmental genes to maintain ES cell identity

Two distinct Polycomb complexes, PRC1 and PRC2, collaborate to maintain epigenetic repression of key developmental loci in embryonic stem cells (ESCs). PRC1 and PRC2 have histone modifying activities, catalyzing mono-ubiquitination of histone H2A (H2AK119u1) and trimethylation of H3 lysine 27 (H3K27me3), respectively. Compared to H3K27me3, localization and the role of H2AK119u1 are not fully understood in ESCs. Here we present genome-wide H2AK119u1 maps in ESCs and identify a group of genes at which H2AK119u1 is deposited in a Ring1-dependent manner. These genes are a distinctive subset of genes with H3K27me3 enrichment and are the central targets of Polycomb silencing that are required to maintain ESC identity. We further show that the H2A ubiquitination activity of PRC1 is dispensable for its target binding and its activity to compact chromatin at Hox loci, but is indispensable for efficient repression of target genes and thereby ESC maintenance. These data demonstrate that multiple effector mechanisms including H2A ubiquitination and chromatin compaction combine to mediate PRC1-dependent repression of genes that are crucial for the maintenance of ESC identity. Utilization of these diverse effector mechanisms might provide a means to maintain a repressive state that is robust yet highly responsive to developmental cues during ES cell self-renewal and differentiation.

Polycomb-group (PcG) proteins play essential roles in the epigenetic regulation of gene expression during development. PcG proteins form two distinct multimeric complexes, PRC1 and PRC2. In the widely accepted hierarchical model, PRC2 is recruited to specific genomic locations and catalyzes trimethylation of H3 lysine 27 (H3K27me3), thereby creating binding sites for PRC1, which then catalyzes mono-ubiquitination of histone H2A (H2AK119u1). Recently, PRC1 has been shown to be able to compact chromatin structure at target loci independently of its histone ubiquitination activity. Therefore, the role of H2AK119u1 still remains unclear. To gain insight into this issue, we used ChIP-on-chip analysis to map H2AK119u1 genome-wide in mouse ES cells (ESCs). The data demonstrate that H2AK119u1 occupies a distinctive subset of genes with H3K27me3 enrichment. These genes are the central targets of Polycomb silencing to maintain ESC identity. We further show that the H2A ubiquitination activity of PRC1 is dispensable for its target binding and its activity to compact chromatin at Hox loci, but is indispensable for efficient repression of target genes and therefore ESC maintenance. We propose that multiple effector mechanisms including H2A ubiquitination and chromatin compaction combine to mediate PRC1-dependent repression of developmental genes to maintain the identity of ESCs.

PRC1 and PRC2 are not required for targeting of H2A.Z to developmental genes in embryonic stem cells

The essential histone variant H2A.Z localises to both active and silent chromatin sites. In embryonic stem cells (ESCs), H2A.Z is also reported to co-localise with polycomb repressive complex 2 (PRC2) at developmentally silenced genes. The mechanism of H2A.Z targeting is not clear, but a role for the PRC2 component Suz12 has been suggested. Given this association, we wished to determine if polycomb functionally directs H2A.Z incorporation in ESCs. We demonstrate that the PRC1 component Ring1B interacts with multiple complexes in ESCs. Moreover, we show that although the genomic distribution of H2A.Z co-localises with PRC2, Ring1B and with the presence of CpG islands, H2A.Z still blankets polycomb target loci in the absence of Suz12, Eed (PRC2) or Ring1B (PRC1). Therefore we conclude that H2A.Z accumulates at developmentally silenced genes in ESCs in a polycomb independent manner.

Histone acetylation and the maintenance of chromatin compaction by Polycomb repressive complexes

Mechanisms controlling higher-order chromatin structure or chromatin compaction and linking this to gene regulation are poorly understood. Previously, we had shown that the PRC1 Polycomb repressive complex is required to maintain a compact chromatin state at Polycomb target loci in embryonic stem cells (ESCs) of the mouse and that this activity, together with the ability to repress target gene expression, is surprisingly independent of the histone ubiquitination activity of the Ring1B component of PRC1. Here we investigate and discuss the role of another histone modification--histone acetylation--in Polycomb function. We show that inhibition of histone deacetylases leads to some decompaction of Hox loci and suggest that histone deacetylation has a role in the pathway of PRC1-mediated chromatin compaction. We discuss whether PRC1 and histone hypoacetylation function together to establish a chromatin template at which stable nucleosomes act to antagonize transcriptional elongation.

Ring1B compacts chromatin structure and represses gene expression independent of histone ubiquitination

How polycomb group proteins repress gene expression in vivo is not known. While histone-modifying activities of the polycomb repressive complexes (PRCs) have been studied extensively, in vitro data have suggested a direct activity of the PRC1 complex in compacting chromatin. Here, we investigate higher-order chromatin compaction of polycomb targets in vivo. We show that PRCs are required to maintain a compact chromatin state at Hox loci in embryonic stem cells (ESCs). There is specific decompaction in the absence of PRC2 or PRC1. This is due to a PRC1-like complex, since decompaction occurs in Ring1B null cells that still have PRC2-mediated H3K27 methylation. Moreover, we show that the ability of Ring1B to restore a compact chromatin state and to repress Hox gene expression is not dependent on its histone ubiquitination activity. We suggest that Ring1B-mediated chromatin compaction acts to directly limit transcription in vivo.

Interaction of HP1 and Brg1/Brm with the globular domain of histone H3 is required for HP1-mediated repression

The heterochromatin-enriched HP1 proteins play a critical role in regulation of transcription. These proteins contain two related domains known as the chromo- and the chromoshadow-domain. The chromo-domain binds histone H3 tails methylated on lysine 9. However, in vivo and in vitro experiments have shown that the affinity of HP1 proteins to native methylated chromatin is relatively poor and that the opening of chromatin occurring during DNA replication facilitates their binding to nucleosomes. These observations prompted us to investigate whether HP1 proteins have additional histone binding activities, envisioning also affinity for regions potentially occluded by the nucleosome structure. We find that the chromoshadow-domain interacts with histone H3 in a region located partially inside the nucleosomal barrel at the entry/exit point of the nucleosome. Interestingly, this region is also contacted by the catalytic subunits of the human SWI/SNF complex. In vitro, efficient SWI/SNF remodeling requires this contact and is inhibited in the presence of HP1 proteins. The antagonism between SWI/SNF and HP1 proteins is also observed in vivo on a series of interferon-regulated genes. Finally, we show that SWI/SNF activity favors loading of HP1 proteins to chromatin both in vivo and in vitro. Altogether, our data suggest that HP1 chromoshadow-domains can benefit from the opening of nucleosomal structures to bind chromatin and that HP1 proteins use this property to detect and arrest unwanted chromatin remodeling.

HP1 proteins are transcriptional regulators frequently associated with gene silencing, a phenomenon involving masking of promoter DNA by dense chromatin. Owing to their chromo-domain, these proteins can read and bind an epigenetic mark that on many non-expressed genes is present on histone H3 at the surface of the nucleosome (the fundamental packing unit of chromatin). However, the binding to this mark does not explain the repressing activity of HP1 proteins. Here, we show that these proteins can establish a second contact with histone H3, independently of the epigenetic mark. This second contact site is located inside the nucleosome, in a position likely to be inaccessible. Interestingly, this site is also contacted by a subunit of the SWI/SNF complex and this contact is required for the ATP-dependent chromatin remodeling catalyzed by SWI/SNF. We provide evidence suggesting that HP1 proteins use the SWI/SNF chromatin remodeling to gain access to the contact site inside the nucleosome and to prevent further remodeling by competing with SWI/SNF for binding at this position. These observations lead us to suggest that HP1 proteins function as gatekeepers on promoters, detecting and stopping unwanted exposure of internal nucleosomal sites.

HP1 binding to chromatin methylated at H3K9 is enhanced by auxiliary factors

A large portion of the eukaryotic genome is packaged into transcriptionally silent heterochromatin. Several factors that play important roles during the establishment and maintenance of this condensed form have been identified. Methylation of lysine 9 within histone H3 and the subsequent binding of the chromodomain protein heterochromatin protein 1 (HP1) are thought to initiate heterochromatin formation in vivo and to propagate a heterochromatic state lasting through several cell divisions. For the present study we analyzed the binding of HP1 to methylated chromatin in a fully reconstituted system. In contrast to its strong binding to methylated peptides, HP1 binds only weakly to methylated chromatin. However, the addition of recombinant SU(VAR) protein, such as ACF1 or SU(VAR)3-9, facilitates HP1 binding to chromatin methylated at lysine 9 within the H3 N terminus (H3K9). We propose that HP1 has multiple target sites that contribute to its recognition of chromatin, only one of them being methylated at H3K9. These findings have implications for the mechanisms of recognition of specific chromatin modifications in vivo.

The Drosophila G9a gene encodes a multi-catalytic histone methyltransferase required for normal development

Mammalian G9a is a histone H3 Lys-9 (H3–K9) methyltransferase localized in euchromatin and acts as a co-regulator for specific transcription factors. G9a is required for proper development in mammals as g9a⁻/g9a⁻ mice show growth retardation and early lethality. Here we describe the cloning, the biochemical and genetical analyses of the Drosophila homolog dG9a. We show that dG9a shares the structural organization of mammalian G9a, and that it is a multi-catalytic histone methyltransferase with specificity not only for lysines 9 and 27 on H3 but also for H4. Surprisingly, it is not the H4–K20 residue that is the target for this methylation. Spatiotemporal expression analyses reveal that dG9a is abundantly expressed in the gonads of both sexes, with no detectable expression in gonadectomized adults. In addition we find a low but clearly observable level of dG9a transcript in developing embryos, larvae and pupae. Genetic and RNAi experiments reveal that dG9a is involved in ecdysone regulatory pathways.

Identification of a specific inhibitor of the histone methyltransferase SU(VAR)3-9

Histone methylation plays a key role in establishing and maintaining stable gene expression patterns during cellular differentiation and embryonic development. Here, we report the characterization of the fungal metabolite chaetocin as the first inhibitor of a lysine-specific histone methyltransferase. Chaetocin is specific for the methyltransferase SU(VAR)3-9 both in vitro and in vivo and may therefore be used to study heterochromatin-mediated gene repression.

Histone mRNAs do not accumulate during S phase of either mitotic or endoreduplicative cycles in the chordate Oikopleura dioica

Metazoan histones are generally classified as replication-dependent or replacement variants. Replication-dependent histone genes contain cell cycle-responsive promoter elements, their transcripts terminate in an unpolyadenylated conserved stem-loop, and their mRNAs accumulate sharply during S phase. Replacement variant genes lack cell cycle-responsive promoter elements, their polyadenylated transcripts lack the stem-loop, and they are expressed at low levels throughout the cell cycle. During early development of some organisms with rapid cleavage cycles, replication-dependent mRNAs are not fully S phase restricted until complete cell cycle regulation is achieved. The accumulation of polyadenylated transcripts during this period has been considered incompatible with metazoan development. We show here that histone metabolism in the urochordate Oikopleura dioica does not accord with some key tenets of the replication-dependent/replacement variant paradigm. During the premetamorphic mitotic phase of development, expressed variants shared characteristics of replication-dependent histones, including the 3' stem-loop, but, in contrast, were extensively polyadenylated. After metamorphosis, when cells in many tissues enter endocycles, there was a global downregulation of histone transcript levels, with most variant transcripts processed at the stem-loop. Contrary to the 30-fold S-phase upregulation of histone transcripts described in common metazoan model organisms, we observed essentially constant histone transcript levels throughout both mitotic and endoreduplicative cell cycles.

The N-Terminus of Drosophila SU(VAR)3−9 Mediates Dimerization and Regulates Its Methyltransferase Activity†

In most eukaryotes, the histone methyltransferase SU(VAR)3-9 and its orthologues play a major role in the function of centromeric heterochromatin. Although the methyltransferase domain is required for the formation of a fully functional centromere, mutations within other regions of the gene such as the N-terminus also have a strong impact on its in vivo function. To analyze the contribution of the N-terminus on the methyltransferase activity, we have expressed the full-length Drosophila SU(VAR)3-9 (dSU(VAR)3-9) together with various N-terminal deletions in Escherichia coli and analyzed the structural and enzymatic properties of the purified recombinant enzymes. Full-length dSU(VAR)3-9 specifically methylates lysine 9 within histone H3 on peptides, on intact histones, and, to a lesser extent, on nucleosomes. A detailed analysis of the reaction products shows that dSU(VAR)3-9 adds two methyl groups to an unmethylated H3 tail peptide in a nonprocessive manner. The full-length enzyme elutes with an apparent molecular weight of 160 kDa from a gel filtration column, which indicates the formation of a dimer. This property is dependent on an intact N-terminus. In contrast to the full-length enzymes, proteins lacking the N-terminus fail to dimerize, and show a 10-fold lower specific activity and a linear dependence of methyltransferase activity on enzyme concentration. A N-terminal peptide containing amino acids 1-152 of dSU(VAR)3-9 is sufficient to mediate this interaction in vitro. The dimerization of dSU(VAR)3-9 and the subsequent increase of its methyltransferase activity provide a starting point to understand the molecular details of the formation of heterochromatic structures in vivo.

Histone gene complement, variant expression, and mRNA processing in a urochordate Oikopleura dioica that undergoes extensive polyploidization

Considerable data exist on coding sequences of histones in a wide variety of organisms. Much more restricted information is available on total histone gene complement, gene organization, transcriptional regulation, and histone mRNA processing. In particular, there is a significant phylogenetic gap in information for the urochordates, a subphylum near the invertebrate-vertebrate transition. In this study, we show that the appendicularian Oikopleura dioica has a histone gene complement that is similar to that of humans, though its genome size is 40- to 50-fold smaller. At a total length of 3.5 kb, the H3, H4, H1, H2A, and H2B quintet cluster is the most compact described thus far, but despite very rapid early developmental cleavage cycles, no extensive tandem repeats of the cluster were present. The high degree of variation within each of the complements of O. dioica H2A and H2B subtypes resembled that found in plants as opposed to more closely related vertebrate and invertebrate species, and developmental stage-specific expression of different subtypes was observed. The linker histone H1 was present in relatively few copies per haploid genome and contained short N- and C-terminal tails, a feature similar to that of copepods but different from many standard model organisms. The 3'UTRs of the histone genes contained both the consensus stem-loop sequence and the polyadenylation signals but lacked the consensus histone downstream element that is involved in the processing of histone mRNAs in echinoderms and vertebrates. Two types of transcripts were found, i.e., those containing both the stem-loop and a polyA tail as well as those cleaved at the normal site just 3' of the stem-loop. The O. dioica data are an important addition to the limited number of eukaryotes for which sufficiently extensive information on histone gene complements is available. Increasingly, it appears that understanding the evolution of histone gene organization, transcriptional regulation, and mRNA processing will depend at least as much on comparative analysis of constraints imposed by certain life history features and cell biological characteristics as on projections based on simple phylogenetic relationships.