Transcribed enhancers lead waves of coordinated transcription in transitioning mammalian cells

Linking FANTOM5 CAGE peaks to annotations with CAGEscan

The FANTOM5 expression atlas is a quantitative measurement of the activity of nearly 200,000 promoter regions across nearly 2,000 different human primary cells, tissue types and cell lines. Generation of this atlas was made possible by the use of CAGE, an experimental approach to localise transcription start sites at single-nucleotide resolution by sequencing the 5′ ends of capped RNAs after their conversion to cDNAs. While 50% of CAGE-defined promoter regions could be confidently associated to adjacent transcriptional units, nearly 100,000 promoter regions remained gene-orphan. To address this, we used the CAGEscan method, in which random-primed 5′-cDNAs are paired-end sequenced. Pairs starting in the same region are assembled in transcript models called CAGEscan clusters. Here, we present the production and quality control of CAGEscan libraries from 56 FANTOM5 RNA sources, which enhances the FANTOM5 expression atlas by providing experimental evidence associating core promoter regions with their cognate transcripts.

Chromatin accessibility underlies synthetic lethality of SWI/SNF subunits in ARID1A-mutant cancers

ARID1A, a subunit of the SWI/SNF chromatin remodeling complex, is frequently mutated in cancer. Deficiency in its homolog ARID1B is synthetically lethal with ARID1A mutation. However, the functional relationship between these homologs has not been explored. Here, we use ATAC-seq, genome-wide histone modification mapping, and expression analysis to examine colorectal cancer cells lacking one or both ARID proteins. We find that ARID1A has a dominant role in maintaining chromatin accessibility at enhancers, while the contribution of ARID1B is evident only in the context of ARID1A mutation. Changes in accessibility are predictive of changes in expression and correlate with loss of H3K4me and H3K27ac marks, nucleosome spacing, and transcription factor binding, particularly at growth pathway genes including MET. We find that ARID1B knockdown in ARID1A mutant ovarian cancer cells causes similar loss of enhancer architecture, suggesting that this is a conserved function underlying the synthetic lethality between ARID1A and ARID1B.

DNA Sequence Constraints Define Functionally Active Steroid Nuclear Receptor Binding Sites in Chromatin

Gene regulatory programs are encoded in the sequence of the DNA. Since the completion of the Human Genome Project, millions of gene regulatory elements have been identified in the human genome. Understanding how each of those sites functionally contributes to gene regulation, however, remains a challenge for nearly every field of biology. Transcription factors influence cell function by interpreting information contained within cis-regulatory elements in chromatin. Whereas chromatin immunoprecipitation-sequencing has been used to identify and map transcription factor-DNA interactions, it has been difficult to assign functionality to the binding sites identified. Thus, in this study, we probed the transcriptional activity, DNA-binding competence, and functional activity of select nuclear receptor mutants in cellular and animal model systems and used this information to define the sequence constraints of functional steroid nuclear receptor cis-regulatory elements. Analysis of the architecture within sNR chromatin interacting sites revealed that only a small fraction of all sNR chromatin-interacting events is associated with transcriptional output and that this functionality is restricted to elements that vary from the consensus palindromic elements by one or two nucleotides. These findings define the transcriptional grammar necessary to predict functionality from regulatory sequences, with a multitude of future implications.

Acute doses of caffeine shift nervous system cell expression profiles toward promotion of neuronal projection growth

Caffeine is a widely consumed psychoactive substance, but little is known about the effects of caffeine stimulation on global gene expression changes in neurons. Here, we conducted gene expression profiling of human neuroepithelial stem cell-derived neurons, stimulated with normal consumption levels of caffeine (3 μM and 10 μM), over a period of 9 h. We found dosage-dependent activation of immediate early genes after 1 h. Neuronal projection development processes were up-regulated and negative regulation of axon extension processes were down-regulated at 3 h. In addition, genes involved in extracellular matrix organization, response for wound healing, and regulation of immune system processes were down-regulated by caffeine at 3 h. This study identified novel genes within the neuronal projection guidance pathways that respond to acute caffeine stimulation and suggests potential mechanisms for the effects of caffeine on neuronal cells.

Transcriptome Analysis Uncovers a Growth-Promoting Activity of Orosomucoid-1 on Hepatocytes

The acute phase protein orosomucoid-1 (Orm1) is mainly expressed by hepatocytes (HPCs) under stress conditions. However, its specific function is not fully understood. Here, we report a role of Orm1 as an executer of HPC proliferation. Increases in serum levels of Orm1 were observed in patients after surgical resection for liver cancer and in mice undergone partial hepatectomy (PH). Transcriptome study showed that Orm1 became the most abundant in HPCs isolated from regenerating mouse liver tissues after PH. Both in vitro and in vivo siRNA-induced knockdown of Orm1 suppressed proliferation of mouse regenerating HPCs and human hepatic cells. Microarray analysis in regenerating mouse livers revealed that the signaling pathways controlling chromatin replication, especially the minichromosome maintenance protein complex genes were uniformly down-regulated following Orm1 knockdown. These data suggest that Orm1 is induced in response to hepatic injury and executes liver regeneration by activating cell cycle progression in HPCs.

•

Serum Orm1 levels increased approximately 1.3- to 2.5-folds in both humans and mice after partial hepatectomy.

•

Transcriptome analysis revealed that Orm1 mostly induced in hepatocytes as a regulator of mouse liver regeneration.

•

Orm1 knockdown in mice impaired liver regeneration with poor hepatocyte growth and suppressed cell cycle signaling.

Orosomucoid-1 (Orm1) is an acute phase protein mainly expressed by hepatocytes under stress conditions. Beginning from the finding that Orm1 was induced after partial hepatectomy in humans and mice, we showed enrichment of Orm1 in regenerating hepatocytes of hepatectomized mice by transcriptome analysis and following culture and animal experiments. Knockdown of Orm1 in mice resulted in decreases in hepatocyte growth accompanying suppressed signaling in controlling chromatin replication. Therefore, Orm1 would be a potential therapeutic and prognostic biomarker for liver diseases, especially after surgical resection of cancer-bearing liver, through its newly found ability to stimulate the cell cycle in regenerating hepatocytes.

Systematic analysis of transcription start sites in avian development

Cap Analysis of Gene Expression (CAGE) in combination with single-molecule sequencing technology allows precision mapping of transcription start sites (TSSs) and genome-wide capture of promoter activities in differentiated and steady state cell populations. Much less is known about whether TSS profiling can characterize diverse and non-steady state cell populations, such as the approximately 400 transitory and heterogeneous cell types that arise during ontogeny of vertebrate animals. To gain such insight, we used the chick model and performed CAGE-based TSS analysis on embryonic samples covering the full 3-week developmental period. In total, 31,863 robust TSS peaks (>1 tag per million [TPM]) were mapped to the latest chicken genome assembly, of which 34% to 46% were active in any given developmental stage. ZENBU, a web-based, open-source platform, was used for interactive data exploration. TSSs of genes critical for lineage differentiation could be precisely mapped and their activities tracked throughout development, suggesting that non-steady state and heterogeneous cell populations are amenable to CAGE-based transcriptional analysis. Our study also uncovered a large set of extremely stable housekeeping TSSs and many novel stage-specific ones. We furthermore demonstrated that TSS mapping could expedite motif-based promoter analysis for regulatory modules associated with stage-specific and housekeeping genes. Finally, using Brachyury as an example, we provide evidence that precise TSS mapping in combination with Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR)-on technology enables us, for the first time, to efficiently target endogenous avian genes for transcriptional activation. Taken together, our results represent the first report of genome-wide TSS mapping in birds and the first systematic developmental TSS analysis in any amniote species (birds and mammals). By facilitating promoter-based molecular analysis and genetic manipulation, our work also underscores the value of avian models in unravelling the complex regulatory mechanism of cell lineage specification during amniote development.

The FANTOM5 collection, a data series underpinning mammalian transcriptome atlases in diverse cell types

The latest project from the FANTOM consortium, an international collaborative effort initiated by RIKEN, generated atlases of transcriptomes, in particular promoters, transcribed enhancers, and long-noncoding RNAs, across a diverse set of mammalian cell types. Here, we introduce the FANTOM5 collection, bringing together data descriptors, articles and analyses of FANTOM5 data published across the Nature Research journals. Associated data are openly available for reuse by all.

FANTOM5 CAGE profiles of human and mouse samples

In the FANTOM5 project, transcription initiation events across the human and mouse genomes were mapped at a single base-pair resolution and their frequencies were monitored by CAGE (Cap Analysis of Gene Expression) coupled with single-molecule sequencing. Approximately three thousands of samples, consisting of a variety of primary cells, tissues, cell lines, and time series samples during cell activation and development, were subjected to a uniform pipeline of CAGE data production. The analysis pipeline started by measuring RNA extracts to assess their quality, and continued to CAGE library production by using a robotic or a manual workflow, single molecule sequencing, and computational processing to generate frequencies of transcription initiation. Resulting data represents the consequence of transcriptional regulation in each analyzed state of mammalian cells. Non-overlapping peaks over the CAGE profiles, approximately 200,000 and 150,000 peaks for the human and mouse genomes, were identified and annotated to provide precise location of known promoters as well as novel ones, and to quantify their activities.

Distinct transcriptional responses of lymphatic endothelial cells to VEGFR-3 and VEGFR-2 stimulation

Vascular endothelial growth factors (VEGFs) and their receptors play crucial roles in the formation of blood and lymphatic vessels during embryogenesis, and also under pathologic conditions in the adult. Despite intensive efforts over the last decades to elucidate the precise functions of VEGFs, transcriptional responses to VEGF receptor stimulation are still not fully characterized. To investigate the specific transcriptional effects of VEGFR-2 and VEGFR-3 activation, we performed a correlation analysis of previously published CAGE sequencing and microarray data of human lymphatic endothelial cells (LECs) stimulated with distinct VEGFs acting through either VEGFR-2 or VEGFR-3. We identified that specific activation of VEGFR-3 by VEGF-C156S results in the downregulation of many genes involved in immune regulation and inflammation, suggesting that VEGFR-3 stimulation has direct anti-inflammatory effects. Comparing CAGE and microarray data sets, we furthermore identified a small number of genes that showed a receptor-dependent response in LECs, demonstrating that these receptors, despite activating very similar signaling pathways, fulfill overlapping but not identical functions within the same cell type (LECs).

FANTOM5 CAGE profiles of human and mouse reprocessed for GRCh38 and GRCm38 genome assemblies

The FANTOM5 consortium described the promoter-level expression atlas of human and mouse by using CAGE (Cap Analysis of Gene Expression) with single molecule sequencing. In the original publications, GRCh37/hg19 and NCBI37/mm9 assemblies were used as the reference genomes of human and mouse respectively; later, the Genome Reference Consortium released newer genome assemblies GRCh38/hg38 and GRCm38/mm10. To increase the utility of the atlas in forthcoming researches, we reprocessed the data to make them available on the recent genome assemblies. The data include observed frequencies of transcription starting sites (TSSs) based on the realignment of CAGE reads, and TSS peaks that are converted from those based on the previous reference. Annotations of the peak names were also updated based on the latest public databases. The reprocessed results enable us to examine frequencies of transcription initiations on the recent genome assemblies and to refer promoters with updated information across the genome assemblies consistently.

An integrated expression atlas of miRNAs and their promoters in human and mouse

MicroRNAs (miRNAs) are short non-coding RNAs with key roles in cellular regulation. As part of the fifth edition of the Functional Annotation of Mammalian Genome (FANTOM5) project, we created an integrated expression atlas of miRNAs and their promoters by deep-sequencing 492 short RNA (sRNA) libraries, with matching Cap Analysis Gene Expression (CAGE) data, from 396 human and 47 mouse RNA samples. Promoters were identified for 1,357 human and 804 mouse miRNAs and showed strong sequence conservation between species. We also found that primary and mature miRNA expression levels were correlated, allowing us to use the primary miRNA measurements as a proxy for mature miRNA levels in a total of 1,829 human and 1,029 mouse CAGE libraries. We thus provide a broad atlas of miRNA expression and promoters in primary mammalian cells, establishing a foundation for detailed analysis of miRNA expression patterns and transcriptional control regions.

A Novel and Multivalent Role of Pax6 in Cerebellar Development

Pax6 is a prominent gene in brain development. The deletion of Pax6 results in devastated development of eye, olfactory bulb, and cortex. However, it has been reported that the Pax6-null Sey cerebellum only has minor defects involving granule cells despite Pax6 being expressed throughout cerebellar development. The present work has uncovered a requirement of Pax6 in the development of all rhombic lip (RL) lineages. A significant downregulation of Tbr1 and Tbr2 expression is found in the Sey cerebellum, these are cell-specific markers of cerebellar nuclear (CN) neurons and unipolar brush cells (UBCs), respectively. The examination of Tbr1 and Lmx1a immunolabeling and Nissl staining confirmed the loss of CN neurons from the Sey cerebellum. CN neuron progenitors are produced in the mutant but there is an enhanced death of these neurons as shown by increased presence of caspase-3-positive cells. These data indicate that Pax6 regulates the survival of CN neuron progenitors. Furthermore, the analysis of experimental mouse chimeras suggests a cell-extrinsic role of Pax6 in CN neuron survival. For UBCs, using Tbr2 immunolabeling, these cells are significantly reduced in the Sey cerebellum. The loss of UBCs in the mutant is due partly to cell death in the RL and also to the reduced production of progenitors from the RL. These results demonstrate a critical role for Pax6 in regulating the generation and survival of UBCs. This and previous work from our laboratory demonstrate a seminal role of Pax6 in the development of all cerebellar glutamatergic neurons.

SIGNIFICANCE STATEMENT

Pax6 is a key molecule in development. Pax6 is best known as the master control gene in eye development with mutations causing aniridia in humans. Pax6 also plays important developmental roles in the cortex and olfactory bulb. During cerebellar development, Pax6 is robustly expressed in the germinal zone of all glutamatergic neurons [cerebellar nuclear (CN) neurons, granule cells, and unipolar brush cells (UBCs)]. Past work has not found abnormalities in the CN and UBC populations. Our study reveals that the Pax6-null mutation dramatically affects these cells and identifies Pax6 as a key regulator of cell survival in CN neurons and of cell production in UBCs. The present study shows how Pax6 is key to the development of glutamatergic cells in the cerebellum.

Transcriptional response to stress is pre-wired by promoter and enhancer architecture

Programs of gene expression are executed by a battery of transcription factors that coordinate divergent transcription from a pair of tightly linked core initiation regions of promoters and enhancers. Here, to investigate how divergent transcription is reprogrammed upon stress, we measured nascent RNA synthesis at nucleotide-resolution, and profiled histone H4 acetylation in human cells. Our results globally show that the release of promoter-proximal paused RNA polymerase into elongation functions as a critical switch at which a gene's response to stress is determined. Highly transcribed and highly inducible genes display strong transcriptional directionality and selective assembly of general transcription factors on the core sense promoter. Heat-induced transcription at enhancers, instead, correlates with prior binding of cell-type, sequence-specific transcription factors. Activated Heat Shock Factor 1 (HSF1) binds to transcription-primed promoters and enhancers, and CTCF-occupied, non-transcribed chromatin. These results reveal chromatin architectural features that orient transcription at divergent regulatory elements and prime transcriptional responses genome-wide.Heat Shock Factor 1 (HSF1) is a regulator of stress-induced transcription. Here, the authors investigate changes to transcription and chromatin organization upon stress and find that activated HSF1 binds to transcription-primed promoters and enhancers, and to CTCF occupied, untranscribed chromatin.

Single-cell profiling reveals that eRNA accumulation at enhancer-promoter loops is not required to sustain transcription

Enhancers are intergenic DNA elements that regulate the transcription of target genes in response to signaling pathways by interacting with promoters over large genomic distances. Recent studies have revealed that enhancers are bi-directionally transcribed into enhancer RNAs (eRNAs). Using single-molecule fluorescence in situ hybridization (smFISH), we investigated the eRNA-mediated regulation of transcription during estrogen induction in MCF-7 cells. We demonstrate that eRNAs are localized exclusively in the nucleus and are induced with similar kinetics as target mRNAs. However, eRNAs are mostly nascent at enhancers and their steady-state levels remain lower than those of their cognate mRNAs. Surprisingly, at the single-allele level, eRNAs are rarely co-expressed with their target loci, demonstrating that active gene transcription does not require the continuous transcription of eRNAs or their accumulation at enhancers. When co-expressed, sub-diffraction distance measurements between nascent mRNA and eRNA signals reveal that co-transcription of eRNAs and mRNAs rarely occurs within closed enhancer–promoter loops. Lastly, basal eRNA transcription at enhancers, but not E2-induced transcription, is maintained upon depletion of MLL1 and ERα, suggesting some degree of chromatin accessibility prior to signal-dependent activation of transcription. Together, our findings suggest that eRNA accumulation at enhancer–promoter loops is not required to sustain target gene transcription.

Initial high-resolution microscopic mapping of active and inactive regulatory sequences proves non-random 3D arrangements in chromatin domain clusters

BACKGROUND

The association of active transcription regulatory elements (TREs) with DNAse I hypersensitivity (DHS[+]) and an 'open' local chromatin configuration has long been known. However, the 3D topography of TREs within the nuclear landscape of individual cells in relation to their active or inactive status has remained elusive. Here, we explored the 3D nuclear topography of active and inactive TREs in the context of a recently proposed model for a functionally defined nuclear architecture, where an active and an inactive nuclear compartment (ANC-INC) form two spatially co-aligned and functionally interacting networks.

RESULTS

Using 3D structured illumination microscopy, we performed 3D FISH with differently labeled DNA probe sets targeting either sites with DHS[+], apparently active TREs, or DHS[-] sites harboring inactive TREs. Using an in-house image analysis tool, DNA targets were quantitatively mapped on chromatin compaction shaped 3D nuclear landscapes. Our analyses present evidence for a radial 3D organization of chromatin domain clusters (CDCs) with layers of increasing chromatin compaction from the periphery to the CDC core. Segments harboring active TREs are significantly enriched at the decondensed periphery of CDCs with loops penetrating into interchromatin compartment channels, constituting the ANC. In contrast, segments lacking active TREs (DHS[-]) are enriched toward the compacted interior of CDCs (INC).

CONCLUSIONS

Our results add further evidence in support of the ANC-INC network model. The different 3D topographies of DHS[+] and DHS[-] sites suggest positional changes of TREs between the ANC and INC depending on their functional state, which might provide additional protection against an inappropriate activation. Our finding of a structural organization of CDCs based on radially arranged layers of different chromatin compaction levels indicates a complex higher-order chromatin organization beyond a dichotomic classification of chromatin into an 'open,' active and 'closed,' inactive state.

Circular RNAs promote TRPM3 expression by inhibiting hsa-miR-130a-3p in coronary artery disease patients

We investigated the differential expression of circular RNAs (circRNAs) in plasma samples from three coronary artery disease (CAD) patients to identify putative therapeutic targets. We identified 24 differentially expressed circRNAs (18 up-regulated and 6 down-regulated) and 7 differentially expressed mRNAs (6 up-regulated and 1 down-regulated) in CAD patients based on competing endogenous RNA (ceRNA) microarray analysis. MiR-221(p = 0.001), miR-155(p = 0.049), and miR-130a (p = 0.001) were downregulated in CAD patients based on qRT-PCR analysis of another independent population of 932 study subjects (648 CAD subjects and 284 controls). We constructed a hsa-miR-130a-3p-mediated circRNA-mRNA ceRNA network using the miRanda database. This included 9 circRNAs (hsa_circ_0089378, hsa_circ_0083357, hsa_circ_0082824, hsa_circ_0068942, hsa_circ_0057576, hsa_circ_0054537, hsa_circ_0051172, hsa_circ_0032970, and hsa_circ_0006323) and 1 mRNA (transient receptor potential cation channel subfamily M member 3 [TRPM3]). We have shown that 9 circRNAs promote TRPM3 expression by inhibiting hsa-miR-130a-3p in CAD patients.

Direct targets of pSTAT5 signalling in erythropoiesis

Erythropoietin (EPO) acts through the dimeric erythropoietin receptor to stimulate proliferation, survival, differentiation and enucleation of erythroid progenitor cells. We undertook two complimentary approaches to find EPO-dependent pSTAT5 target genes in murine erythroid cells: RNA-seq of newly transcribed (4sU-labelled) RNA, and ChIP-seq for pSTAT5 30 minutes after EPO stimulation. We found 302 pSTAT5-occupied sites: ~15% of these reside in promoters while the rest reside within intronic enhancers or intergenic regions, some >100kb from the nearest TSS. The majority of pSTAT5 peaks contain a central palindromic GAS element, TTCYXRGAA. There was significant enrichment for GATA motifs and CACCC-box motifs within the neighbourhood of pSTAT5-bound peaks, and GATA1 and/or KLF1 co-occupancy at many sites. Using 4sU-RNA-seq we determined the EPO-induced transcriptome and validated differentially expressed genes using dynamic CAGE data and qRT-PCR. We identified known direct pSTAT5 target genes such as Bcl2l1, Pim1 and Cish, and many new targets likely to be involved in driving erythroid cell differentiation including those involved in mRNA splicing (Rbm25), epigenetic regulation (Suv420h2), and EpoR turnover (Clint1/EpsinR). Some of these new EpoR-JAK2-pSTAT5 target genes could be used as biomarkers for monitoring disease activity in polycythaemia vera, and for monitoring responses to JAK inhibitors.

Models of global gene expression define major domains of cell type and tissue identity

The current classification of cells in an organism is largely based on their anatomic and developmental origin. Cells types and tissues are traditionally classified into those that arise from the three embryonic germ layers, the ectoderm, mesoderm and endoderm, but this model does not take into account the organization of cell type-specific patterns of gene expression. Here, we present computational models for cell type and tissue specification derived from a collection of 921 RNA-sequencing samples from 272 distinct mouse cell types or tissues. In an unbiased fashion, this analysis accurately predicts the three known germ layers. Unexpectedly, this analysis also suggests that in total there are eight major domains of cell type-specification, corresponding to the neurectoderm, neural crest, surface ectoderm, endoderm, mesoderm, blood mesoderm, germ cells and the embryonic domain. Further, we identify putative genes responsible for specifying the domain and the cell type. This model has implications for understanding trans-lineage differentiation for stem cells, developmental cell biology and regenerative medicine.

A Pentanucleotide ATTTC Repeat Insertion in the Non-coding Region of DAB1, Mapping to SCA37, Causes Spinocerebellar Ataxia

Advances in human genetics in recent years have largely been driven by next-generation sequencing (NGS); however, the discovery of disease-related gene mutations has been biased toward the exome because the large and very repetitive regions that characterize the non-coding genome remain difficult to reach by that technology. For autosomal-dominant spinocerebellar ataxias (SCAs), 28 genes have been identified, but only five SCAs originate from non-coding mutations. Over half of SCA-affected families, however, remain without a genetic diagnosis. We used genome-wide linkage analysis, NGS, and repeat analysis to identify an (ATTTC)_n insertion in a polymorphic ATTTT repeat in DAB1 in chromosomal region 1p32.2 as the cause of autosomal-dominant SCA; this region has been previously linked to SCA37. The non-pathogenic and pathogenic alleles have the configurations [(ATTTT)_7-400] and [(ATTTT)_60-79(ATTTC)_31-75(ATTTT)_58-90], respectively. (ATTTC)_n insertions are present on a distinct haplotype and show an inverse correlation between size and age of onset. In the DAB1-oriented strand, (ATTTC)_n is located in 5' UTR introns of cerebellar-specific transcripts arising mostly during human fetal brain development from the usage of alternative promoters, but it is maintained in the adult cerebellum. Overexpression of the transfected (ATTTC)₅₈ insertion, but not (ATTTT)_n, leads to abnormal nuclear RNA accumulation. Zebrafish embryos injected with RNA of the (AUUUC)₅₈ insertion, but not (AUUUU)_n, showed lethal developmental malformations. Together, these results establish an unstable repeat insertion in DAB1 as a cause of cerebellar degeneration; on the basis of the genetic and phenotypic evidence, we propose this mutation as the molecular basis for SCA37.

Constrained release of lamina-associated enhancers and genes from the nuclear envelope during T-cell activation facilitates their association in chromosome compartments

The 3D organization of the genome changes concomitantly with expression changes during hematopoiesis and immune activation. Studies have focused either on lamina-associated domains (LADs) or on topologically associated domains (TADs), defined by preferential local chromatin interactions, and chromosome compartments, defined as higher-order interactions between TADs sharing functionally similar states. However, few studies have investigated how these affect one another. To address this, we mapped LADs using Lamin B1-DamID during Jurkat T-cell activation, finding significant genome reorganization at the nuclear periphery dominated by release of loci frequently important for T-cell function. To assess how these changes at the nuclear periphery influence wider genome organization, our DamID data sets were contrasted with TADs and compartments. Features of specific repositioning events were then tested by fluorescence in situ hybridization during T-cell activation. First, considerable overlap between TADs and LADs was observed with the TAD repositioning as a unit. Second, A1 and A2 subcompartments are segregated in 3D space through differences in proximity to LADs along chromosomes. Third, genes and a putative enhancer in LADs that were released from the periphery during T-cell activation became preferentially associated with A2 subcompartments and were constrained to the relative proximity of the lamina. Thus, lamina associations influence internal nuclear organization, and changes in LADs during T-cell activation may provide an important additional mode of gene regulation.

Krüppel-like factors compete for promoters and enhancers to fine-tune transcription

Krüppel-like factors (KLFs) are a family of 17 transcription factors characterized by a conserved DNA-binding domain of three zinc fingers and a variable N-terminal domain responsible for recruiting cofactors. KLFs have diverse functions in stem cell biology, embryo patterning, and tissue homoeostasis. KLF1 and related family members function as transcriptional activators via recruitment of co-activators such as EP300, whereas KLF3 and related members act as transcriptional repressors via recruitment of C-terminal Binding Proteins. KLF1 directly activates the Klf3 gene via an erythroid-specific promoter. Herein, we show KLF1 and KLF3 bind common as well as unique sites within the erythroid cell genome by ChIP-seq. We show KLF3 can displace KLF1 from key erythroid gene promoters and enhancers in vivo. Using 4sU RNA labelling and RNA-seq, we show this competition results in reciprocal transcriptional outputs for >50 important genes. Furthermore, Klf3-/- mice displayed exaggerated recovery from anemic stress and persistent cell cycling consistent with a role for KLF3 in dampening KLF1-driven proliferation. We suggest this study provides a paradigm for how KLFs work in incoherent feed-forward loops or networks to fine-tune transcription and thereby control diverse biological processes such as cell proliferation.

Chromatin-enriched lncRNAs can act as cell-type specific activators of proximal gene transcription

We recently described a new class of long noncoding RNA defined by especially tight chromatin association, whose presence is strongly correlated with expression of nearby genes in HEK293 cells. Here we critically examine the generality and cis-enhancer mechanism of this class of chromatin enriched RNA (cheRNA). CheRNA are largely cell-type specific, and remain the most effective chromatin signature for predicting cis-gene transcription in all cell types examined. Targeted depletion of three cheRNAs decreases gene expression of their neighbors, indicating potential co-activator function. Single-molecule FISH of one cheRNA-distal target gene pair suggests spatial overlap consistent with a role in chromosome looping. In another example, the cheRNA HIDALGO stimulates the fetal hemoglobin HBG1 gene during erythroid differentiation by promoting contacts to a downstream enhancer. Our results suggest that many cheRNAs activate proximal, lineage-specific gene transcription.

Relatively frequent switching of transcription start sites during cerebellar development

BACKGROUND

Alternative transcription start site (TSS) usage plays important roles in transcriptional control of mammalian gene expression. The growing interest in alternative TSSs and their role in genome diversification spawned many single-gene studies on differential usages of tissue-specific or temporal-specific alternative TSSs. However, exploration of the switching usage of alternative TSS usage on a genomic level, especially in the central nervous system, is largely lacking.

RESULTS

In this study, We have prepared a unique set of time-course data for the developing cerebellum, as part of the FANTOM5 consortium ( http://fantom.gsc.riken.jp/5/ ) that uses their innovative capturing of 5' ends of all transcripts followed by Helicos next generation sequencing. We analyzed the usage of all transcription start sites (TSSs) at each time point during cerebellar development that provided information on multiple RNA isoforms that emerged from the same gene. We developed a mathematical method that systematically compares the expression of different TSSs of a gene to identify temporal crossover and non-crossover switching events. We identified 48,489 novel TSS switching events in 5433 genes during cerebellar development. This includes 9767 crossover TSS switching events in 1511 genes, where the dominant TSS shifts over time.

CONCLUSIONS

We observed a relatively high prevalence of TSS switching in cerebellar development where the resulting temporally-specific gene transcripts and protein products can play important regulatory and functional roles.

Chromatin loops and causality loops: the influence of RNA upon spatial nuclear architecture

An intrinsic and essential trait exhibited by cells is the properly coordinated and integrated regulation of an astoundingly large number of simultaneous molecular decisions and reactions to maintain biochemical homeostasis. This is especially true inside the cell nucleus, where the recognition of DNA and RNA by a vast range of nucleic acid-interacting proteins organizes gene expression patterns. However, this dynamic system is not regulated by simple "on" or "off" signals. Instead, transcription factor and RNA polymerase recruitment to DNA are influenced by the local chromatin and epigenetic environment, a gene's relative position within the nucleus and the action of noncoding RNAs. In addition, major phase-separated structural features of the nucleus, such as nucleoli and paraspeckles, assemble in direct response to specific transcriptional activities and, in turn, influence global genomic function. Currently, the interpretation of these data is trapped in a causality dilemma reminiscent of the "chicken and the egg" paradox as it is unclear whether changes in nuclear architecture promote RNA function or vice versa. Here, we review recent advances that suggest a complex and interdependent interaction network between gene expression, chromatin topology, and noncoding RNA function. We also discuss the functional links between these essential nuclear processes from the nanoscale (gene looping) to the macroscale (sub-nuclear gene positioning and nuclear body function) and briefly highlight some of the challenges that researchers may encounter when studying these phenomena.

Identification of genetic variants affecting vitamin D receptor binding and associations with autoimmune disease

Large numbers of statistically significant associations between sentinel SNPs and case-control status have been replicated by genome-wide association studies. Nevertheless, few underlying molecular mechanisms of complex disease are currently known. We investigated whether variation in binding of a transcription factor, the vitamin D receptor (VDR), whose activating ligand vitamin D has been proposed as a modifiable factor in multiple disorders, could explain any of these associations. VDR modifies gene expression by binding DNA as a heterodimer with the Retinoid X receptor (RXR). We identified 43,332 genetic variants significantly associated with altered VDR binding affinity (VDR-BVs) using a high-resolution (ChIP-exo) genome-wide analysis of 27 HapMap lymphoblastoid cell lines. VDR-BVs are enriched in consensus RXR::VDR binding motifs, yet most fell outside of these motifs, implying that genetic variation often affects the binding affinity only indirectly. Finally, we compared 341 VDR-BVs replicating by position in multiple individuals against background sets of variants lying within VDR-binding regions that had been matched in allele frequency and were independent with respect to linkage disequilibrium. In this stringent test, these replicated VDR-BVs were significantly (q < 0.1) and substantially (>2-fold) enriched in genomic intervals associated with autoimmune and other diseases, including inflammatory bowel disease, Crohn's disease and rheumatoid arthritis. The approach's validity is underscored by RXR::VDR motif sequence being predictive of binding strength and being evolutionarily constrained. Our findings are consistent with altered RXR::VDR binding contributing to immunity-related diseases. Replicated VDR-BVs associated with these disorders could represent causal disease risk alleles whose effect may be modifiable by vitamin D levels.

The identification and functional annotation of RNA structures conserved in vertebrates

Structured elements of RNA molecules are essential in, e.g., RNA stabilization, localization, and protein interaction, and their conservation across species suggests a common functional role. We computationally screened vertebrate genomes for conserved RNA structures (CRSs), leveraging structure-based, rather than sequence-based, alignments. After careful correction for sequence identity and GC content, we predict ∼516,000 human genomic regions containing CRSs. We find that a substantial fraction of human–mouse CRS regions (1) colocalize consistently with binding sites of the same RNA binding proteins (RBPs) or (2) are transcribed in corresponding tissues. Additionally, a CaptureSeq experiment revealed expression of many of our CRS regions in human fetal brain, including 662 novel ones. For selected human and mouse candidate pairs, qRT-PCR and in vitro RNA structure probing supported both shared expression and shared structure despite low abundance and low sequence identity. About 30,000 CRS regions are located near coding or long noncoding RNA genes or within enhancers. Structured (CRS overlapping) enhancer RNAs and extended 3′ ends have significantly increased expression levels over their nonstructured counterparts. Our findings of transcribed uncharacterized regulatory regions that contain CRSs support their RNA-mediated functionality.

Biotagging of Specific Cell Populations in Zebrafish Reveals Gene Regulatory Logic Encoded in the Nuclear Transcriptome

Interrogation of gene regulatory circuits in complex organisms requires precise tools for the selection of individual cell types and robust methods for biochemical profiling of target proteins. We have developed a versatile, tissue-specific binary in vivo biotinylation system in zebrafish termed biotagging that uses genetically encoded components to biotinylate target proteins, enabling in-depth genome-wide analyses of their molecular interactions. Using tissue-specific drivers and cell-compartment-specific effector lines, we demonstrate the specificity of the biotagging toolkit at the biochemical, cellular, and transcriptional levels. We use biotagging to characterize the in vivo transcriptional landscape of migratory neural crest and myocardial cells in different cellular compartments (ribosomes and nucleus). These analyses reveal a comprehensive network of coding and non-coding RNAs and cis-regulatory modules, demonstrating that tissue-specific identity is embedded in the nuclear transcriptomes. By eliminating background inherent to complex embryonic environments, biotagging allows analyses of molecular interactions at high resolution.

•

Biotagging enables cell- and compartment-specific in vivo biotinylation in zebrafish

•

Technique yields comprehensive nuclear transcriptional analysis of cardiomyocytes

•

Biotagging finds bidirectionally transcribed neural crest cis-regulatory modules

•

System reveals tissue-specific regulation of noncoding RNA species

Enhancers and chromatin structures: regulatory hubs in gene expression and diseases

Gene expression requires successful communication between enhancer and promoter regions, whose activities are regulated by a variety of factors and associated with distinct chromatin structures; in addition, functionally related genes and their regulatory repertoire tend to be arranged in the same subchromosomal regulatory domains. In this review, we discuss the importance of enhancers, especially clusters of enhancers (such as super-enhancers), as key regulatory hubs to integrate environmental cues and encode spatiotemporal instructions for genome expression, which are critical for a variety of biological processes governing mammalian development. Furthermore, we emphasize that the enhancer–promoter interaction landscape provides a critical context to understand the aetiologies and mechanisms behind numerous complex human diseases and provides new avenues for effective transcription-based interventions.

Identification of Gene Transcription Start Sites and Enhancers Responding to Pulmonary Carbon Nanotube Exposure in Vivo

Increased use of nanomaterials in industry, medicine, and consumer products has raised concerns over their toxicity. To ensure safe use of nanomaterials, understanding their biological effects at the molecular level is crucial. In particular, the regulatory mechanisms responsible for the cascade of genes activated by nanomaterial exposure are not well-characterized. To this end, we profiled the genome-wide usage of gene transcription start sites and linked active enhancer regions in lungs of C57BL/6 mice 24 h after intratracheal instillation of a single dose of the multiwalled carbon nanotube (MWCNT) Mitsui-7. Our results revealed a massive gene regulatory response, where expression of key inflammatory genes (e.g., Csf3, Il24, and Fgf23) was increased >100-fold 24 h after Mitsui-7 exposure. Many of the Mitsui-7-responsive transcription start sites were alternative transcription start sites for known genes, and the number of alternative transcription start sites used in a given gene was correlated with overall Mitsui-7 response. Strikingly, genes that were up-regulated after Mitsui-7 exposure only through their main annotated transcription start site were linked to inflammatory and defense responses, while genes up-regulated only through alternative transcription start sites were functionally heterogeneous and not inflammation-associated. Furthermore, we identified almost 12 000 active enhancers, many of which were Mitsui-7-responsive, and we identified similarly responding putative target genes. Overall, our study provides the location and activity of Mitsui-7-induced enhancers and transcription start sites, providing a useful resource for targeted experiments elucidating the biological effects of nanomaterials and the identification of biomarkers for early detection of MWCNT-induced inflammation.

Correlation of EGFR or KRAS mutation status with 18F-FDG uptake on PET-CT scan in lung adenocarcinoma

Background

¹⁸F-fluoro-2-deoxy-glucose (¹⁸F-FDG) positron emission tomography (PET) is a functional imaging modality based on glucose metabolism. The correlation between EGFR or KRAS mutation status and the standardized uptake value (SUV) of ¹⁸F-FDG PET scanning has not been fully elucidated.

Methods

Correlations between EGFR or KRAS mutation status and clinicopathological factors including SUV_max were statistically analyzed in 734 surgically resected lung adenocarcinoma patients. Molecular causal relationships between EGFR or KRAS mutation status and glucose metabolism were then elucidated in 62 lung adenocarcinomas using cap analysis of gene expression (CAGE), a method to determine and quantify the transcription initiation activities of mRNA across the genome.

Results

EGFR and KRAS mutations were detected in 334 (46%) and 83 (11%) of the 734 lung adenocarcinomas, respectively. The remaining 317 (43%) patients had wild-type tumors for both genes. EGFR mutations were more frequent in tumors with lower SUV_max. In contrast, no relationship was noted between KRAS mutation status and SUV_max. CAGE revealed that 4 genes associated with glucose metabolism (GPI, G6PD, PKM2, and GAPDH) and 5 associated with the cell cycle (ANLN, PTTG1, CIT, KPNA2, and CDC25A) were positively correlated with SUV_max, although expression levels were lower in EGFR-mutated than in wild-type tumors. No similar relationships were noted with KRAS mutations.

Conclusions

EGFR-mutated adenocarcinomas are biologically indolent with potentially lower levels of glucose metabolism than wild-type tumors. Several genes associated with glucose metabolism and the cell cycle were specifically down-regulated in EGFR-mutated adenocarcinomas.

Super-Enhancer-Driven Transcriptional Dependencies in Cancer

Transcriptional deregulation is one of the core tenets of cancer biology and is underpinned by alterations in both protein-coding genes and noncoding regulatory elements. Large regulatory elements, so-called super-enhancers (SEs), are central to the maintenance of cancer cell identity and promote oncogenic transcription to which cancer cells become highly addicted. Such dependence on SE-driven transcription for proliferation and survival offers an Achilles heel for the therapeutic targeting of cancer cells. Indeed, inhibition of the cellular machinery required for the assembly and maintenance of SEs dampens oncogenic transcription and inhibits tumor growth. In this article, we review the organization, function, and regulation of oncogenic SEs and their contribution to the cancer cell state.

Psip1/p52 regulates posterior Hoxa genes through activation of lncRNA Hottip

Long noncoding RNAs (lncRNAs) have been implicated in various biological functions including the regulation of gene expression, however, the functionality of lncRNAs is not clearly understood and conflicting conclusions have often been reached when comparing different methods to investigate them. Moreover, little is known about the upstream regulation of lncRNAs. Here we show that the short isoform (p52) of a transcriptional co-activator-PC4 and SF2 interacting protein (Psip1), which is known to be involved in linking transcription to RNA processing, specifically regulates the expression of the lncRNA Hottip-located at the 5' end of the Hoxa locus. Using both knockdown and knockout approaches we show that Hottip expression is required for activation of the 5' Hoxa genes (Hoxa13 and Hoxa10/11) and for retaining Mll1 at the 5' end of Hoxa. Moreover, we demonstrate that artificially inducing Hottip expression is sufficient to activate the 5' Hoxa genes and that Hottip RNA binds to the 5' end of Hoxa. By engineering premature transcription termination, we show that it is the Hottip lncRNA molecule itself, not just Hottip transcription that is required to maintains active expression of posterior Hox genes. Our data show a direct role for a lncRNA molecule in regulating the expression of developmentally-regulated mRNA genes in cis.

Interplay between genetic and epigenetic mechanisms in rheumatoid arthritis

Genetic and environmental factors contribute to the risk for rheumatoid arthritis (RA), with epigenetics serving as a possible interface through which risk factors contribute to RA. High-throughput technologies for interrogating genome and epigenome, and the availability of genetic and epigenetic datasets across a diversity of cell types, enable the identification of candidate causal genetic variants for RA to study their function in core RA processes. To date, RA risk variants were studied in the immune cells but not joint resident cells, for example, synovial fibroblasts. Synovial fibroblasts from different joints are distinct, anatomically specialized cells, defined by joint-specific transcriptomes, epigenomes and phenotypes. Cell type-specific analysis of epigenetic changes, together with genetic fine mapping and interrogation of chromatin 3D interactions may identify new disease relevant pathways, potential therapeutic targets and biomarkers for RA progression or therapy response.

TT-seq captures enhancer landscapes immediately after T-cell stimulation

To monitor transcriptional regulation in human cells, rapid changes in enhancer and promoter activity must be captured with high sensitivity and temporal resolution. Here, we show that the recently established protocol TT-seq ("transient transcriptome sequencing") can monitor rapid changes in transcription from enhancers and promoters during the immediate response of T cells to ionomycin and phorbol 12-myristate 13-acetate (PMA). TT-seq maps eRNAs and mRNAs every 5 min after T-cell stimulation with high sensitivity and identifies many new primary response genes. TT-seq reveals that the synthesis of 1,601 eRNAs and 650 mRNAs changes significantly within only 15 min after stimulation, when standard RNA-seq does not detect differentially expressed genes. Transcription of enhancers that are primed for activation by nucleosome depletion can occur immediately and simultaneously with transcription of target gene promoters. Our results indicate that enhancer transcription is a good proxy for enhancer regulatory activity in target gene activation, and establish TT-seq as a tool for monitoring the dynamics of enhancer landscapes and transcription programs during cellular responses and differentiation.

Analysis of the human monocyte-derived macrophage transcriptome and response to lipopolysaccharide provides new insights into genetic aetiology of inflammatory bowel disease

The FANTOM5 consortium utilised cap analysis of gene expression (CAGE) to provide an unprecedented insight into transcriptional regulation in human cells and tissues. In the current study, we have used CAGE-based transcriptional profiling on an extended dense time course of the response of human monocyte-derived macrophages grown in macrophage colony-stimulating factor (CSF1) to bacterial lipopolysaccharide (LPS). We propose that this system provides a model for the differentiation and adaptation of monocytes entering the intestinal lamina propria. The response to LPS is shown to be a cascade of successive waves of transient gene expression extending over at least 48 hours, with hundreds of positive and negative regulatory loops. Promoter analysis using motif activity response analysis (MARA) identified some of the transcription factors likely to be responsible for the temporal profile of transcriptional activation. Each LPS-inducible locus was associated with multiple inducible enhancers, and in each case, transient eRNA transcription at multiple sites detected by CAGE preceded the appearance of promoter-associated transcripts. LPS-inducible long non-coding RNAs were commonly associated with clusters of inducible enhancers. We used these data to re-examine the hundreds of loci associated with susceptibility to inflammatory bowel disease (IBD) in genome-wide association studies. Loci associated with IBD were strongly and specifically (relative to rheumatoid arthritis and unrelated traits) enriched for promoters that were regulated in monocyte differentiation or activation. Amongst previously-identified IBD susceptibility loci, the vast majority contained at least one promoter that was regulated in CSF1-dependent monocyte-macrophage transitions and/or in response to LPS. On this basis, we concluded that IBD loci are strongly-enriched for monocyte-specific genes, and identified at least 134 additional candidate genes associated with IBD susceptibility from reanalysis of published GWA studies. We propose that dysregulation of monocyte adaptation to the environment of the gastrointestinal mucosa is the key process leading to inflammatory bowel disease.

Macrophages are immune cells that form the first line of defense against pathogens, but also mediate tissue damage in inflammatory disease. Macrophages initiate inflammation by recognising and responding to components of bacterial cells. Macrophages of the wall of the gut are constantly replenished from the blood. Upon entering the intestine, newly-arrived cells modulate their response to stimuli derived from the bacteria in the wall of the gut. This process fails in chronic inflammatory bowel diseases (IBD). Both the major forms of IBD, Crohn’s disease and ulcerative colitis, run in families. The inheritance is complex, involving more than 200 different regions of the genome. We hypothesised that the genetic risk of IBD is associated specifically with altered regulation of genes that control the development of macrophages. In this study, we used the comprehensive transcriptome dataset produced by the FANTOM5 consortium to identify the sets of promoters and enhancers that are involved in adaptation of macrophages to the gut wall, their response to bacterial stimuli, and how their functions are integrated. A reanalysis of published genome-wide association data based upon regulated genes in monocytes as candidates strongly supports the view that susceptibility to IBD arises from a primary defect in macrophage differentiation.

An atlas of human long non-coding RNAs with accurate 5' ends

Long non-coding RNAs (lncRNAs) are largely heterogeneous and functionally uncharacterized. Here, using FANTOM5 cap analysis of gene expression (CAGE) data, we integrate multiple transcript collections to generate a comprehensive atlas of 27,919 human lncRNA genes with high-confidence 5' ends and expression profiles across 1,829 samples from the major human primary cell types and tissues. Genomic and epigenomic classification of these lncRNAs reveals that most intergenic lncRNAs originate from enhancers rather than from promoters. Incorporating genetic and expression data, we show that lncRNAs overlapping trait-associated single nucleotide polymorphisms are specifically expressed in cell types relevant to the traits, implicating these lncRNAs in multiple diseases. We further demonstrate that lncRNAs overlapping expression quantitative trait loci (eQTL)-associated single nucleotide polymorphisms of messenger RNAs are co-expressed with the corresponding messenger RNAs, suggesting their potential roles in transcriptional regulation. Combining these findings with conservation data, we identify 19,175 potentially functional lncRNAs in the human genome.

Discriminative identification of transcriptional responses of promoters and enhancers after stimulus

Promoters and enhancers regulate the initiation of gene expression and maintenance of expression levels in spatial and temporal manner. Recent findings stemming from the Cap Analysis of Gene Expression (CAGE) demonstrate that promoters and enhancers, based on their expression profiles after stimulus, belong to different transcription response subclasses. One of the most promising biological features that might explain the difference in transcriptional response between subclasses is the local chromatin environment. We introduce a novel computational framework, PEDAL, for distinguishing effectively transcriptional profiles of promoters and enhancers using solely histone modification marks, chromatin accessibility and binding sites of transcription factors and co-activators. A case study on data from MCF-7 cell-line reveals that PEDAL can identify successfully the transcription response subclasses of promoters and enhancers from two different stimulations. Moreover, we report subsets of input markers that discriminate with minimized classification error MCF-7 promoter and enhancer transcription response subclasses. Our work provides a general computational approach for identifying effectively cell-specific and stimulation-specific promoter and enhancer transcriptional profiles, and thus, contributes to improve our understanding of transcriptional activation in human.

Retinoic Acid Negatively Impacts Proliferation and MCTC Specific Attributes of Human Skin Derived Mast Cells, but Reinforces Allergic Stimulability

The Vitamin-A-metabolite retinoic acid (RA) acts as a master regulator of cellular programs. Mast cells (MCs) are primary effector cells of type-I-allergic reactions. We recently uncovered that human cutaneous MCs are enriched with RA network components over other skin cells. Yet, direct experimental evidence on the significance of the RA-MC axis is limited. Here, skin-derived cultured MCs were exposed to RA for seven days and investigated by flow-cytometry (BrdU incorporation, Annexin/PI, FcεRI), microscopy, RT-qPCR, histamine quantitation, protease activity, and degranulation assays. We found that while MC size and granularity remained unchanged, RA potently interfered with MC proliferation. Conversely, a modest survival-promoting effect from RA was noted. The granule constituents, histamine and tryptase, remained unaffected, while RA had a striking impact on MC chymase, whose expression dropped by gene and by peptidase activity. The newly uncovered MRGPRX2 performed similarly to chymase. Intriguingly, RA fostered allergic MC degranulation, in a way completely uncoupled from FcεRI expression, but it simultaneously restricted MRGPRX2-triggered histamine release in agreement with the reduced receptor expression. Vitamin-A-derived hormones thus re-shape skin-derived MCs numerically, phenotypically, and functionally. A general theme emerges, implying RA to skew MCs towards processes associated with (allergic) inflammation, while driving them away from the skin-imprinted MC_TC (“MCs containing tryptase and chymase”) signature (chymase, MRGPRX2). Collectively, MCs are substantial targets of the skin retinoid network.

Altered enhancer transcription underlies Huntington's disease striatal transcriptional signature

Epigenetic and transcriptional alterations are both implicated in Huntington’s disease (HD), a progressive neurodegenerative disease resulting in degeneration of striatal neurons in the brain. However, how impaired epigenetic regulation leads to transcriptional dysregulation in HD is unclear. Here, we investigated enhancer RNAs (eRNAs), a class of long non-coding RNAs transcribed from active enhancers. We found that eRNAs are expressed from many enhancers of mouse striatum and showed that a subset of those eRNAs are deregulated in HD vs control mouse striatum. Enhancer regions producing eRNAs decreased in HD mouse striatum were associated with genes involved in striatal neuron identity. Consistently, they were enriched in striatal super-enhancers. Moreover, decreased eRNA expression in HD mouse striatum correlated with down-regulation of associated genes. Additionally, a significant number of RNA Polymerase II (RNAPII) binding sites were lost within enhancers associated with decreased eRNAs in HD vs control mouse striatum. Together, this indicates that loss of RNAPII at HD mouse enhancers contributes to reduced transcription of eRNAs, resulting in down-regulation of target genes. Thus, our data support the view that eRNA dysregulation in HD striatum is a key mechanism leading to altered transcription of striatal neuron identity genes, through reduced recruitment of RNAPII at super-enhancers.

Guidance of regulatory T cell development by Satb1-dependent super-enhancer establishment

Most Foxp3+ regulatory T (Treg) cells develop in the thymus as a functionally mature T cell subpopulation specialized for immune suppression. Their cell fate appears to be determined before Foxp3 expression; yet molecular events that prime Foxp3- Treg precursor cells are largely obscure. We found that Treg cell-specific super-enhancers (Treg-SEs), which were associated with Foxp3 and other Treg cell signature genes, began to be activated in Treg precursor cells. T cell-specific deficiency of the genome organizer Satb1 impaired Treg-SE activation and the subsequent expression of Treg signature genes, causing severe autoimmunity due to Treg cell deficiency. These results suggest that Satb1-dependent Treg-SE activation is crucial for Treg cell lineage specification in the thymus and that its perturbation is causative of autoimmune and other immunological diseases.

Sex-specific chromatin landscapes in an ultra-compact chordate genome

Background

In multicellular organisms, epigenome dynamics are associated with transitions in the cell cycle, development, germline specification, gametogenesis and inheritance. Evolutionarily, regulatory space has increased in complex metazoans to accommodate these functions. In tunicates, the sister lineage to vertebrates, we examine epigenome adaptations to strong secondary genome compaction, sex chromosome evolution and cell cycle modes.

Results

Across the 70 MB Oikopleura dioica genome, we profiled 19 histone modifications, and RNA polymerase II, CTCF and p300 occupancies, to define chromatin states within two homogeneous tissues with distinct cell cycle modes: ovarian endocycling nurse nuclei and mitotically proliferating germ nuclei in testes. Nurse nuclei had active chromatin states similar to other metazoan epigenomes, with large domains of operon-associated transcription, a general lack of heterochromatin, and a possible role of Polycomb PRC2 in dosage compensation. Testis chromatin states reflected transcriptional activity linked to spermatogenesis and epigenetic marks that have been associated with establishment of transgenerational inheritance in other organisms. We also uncovered an unusual chromatin state specific to the Y-chromosome, which combined active and heterochromatic histone modifications on specific transposable elements classes, perhaps involved in regulating their activity.

Conclusions

Compacted regulatory space in this tunicate genome is accompanied by reduced heterochromatin and chromatin state domain widths. Enhancers, promoters and protein-coding genes have conserved epigenomic features, with adaptations to the organization of a proportion of genes in operon units. We further identified features specific to sex chromosomes, cell cycle modes, germline identity and dosage compensation, and unusual combinations of histone PTMs with opposing consensus functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0110-4) contains supplementary material, which is available to authorized users.

SOX9 chromatin folding domains correlate with its real and putative distant cis-regulatory elements

Evolutionary conserved transcription factor SOX9, encoded by the dosage sensitive SOX9 gene on chromosome 17q24.3, plays an important role in development of multiple organs, including bones and testes. Heterozygous point mutations and genomic copy-number variant (CNV) deletions involving SOX9 have been reported in patients with campomelic dysplasia (CD), a skeletal malformation syndrome often associated with male-to-female sex reversal. Balanced and unbalanced structural genomic variants with breakpoints mapping up to 1.3 Mb up- and downstream to SOX9 have been described in patients with milder phenotypes, including acampomelic campomelic dysplasia, sex reversal, and Pierre Robin sequence. Based on the localization of breakpoints of genomic rearrangements causing different phenotypes, 5 genomic intervals mapping upstream to SOX9 have been defined. We have analyzed the publically available database of high-throughput chromosome conformation capture (Hi-C) in multiple cell lines in the genomic regions flanking SOX9. Consistent with the literature data, chromatin domain boundaries in the SOX9 locus exhibit conservation across species and remain largely constant across multiple cell types. Interestingly, we have found that chromatin folding domains in the SOX9 locus associate with the genomic intervals harboring real and putative regulatory elements of SOX9, implicating that variation in intra-domain interactions may be critical for dynamic regulation of SOX9 expression in a cell type-specific fashion. We propose that tissue-specific enhancers for other transcription factor genes may similarly utilize chromatin folding sub-domains in gene regulation.

Transcriptional Dynamics During Human Adipogenesis and Its Link to Adipose Morphology and Distribution

White adipose tissue (WAT) can develop into several phenotypes with different pathophysiological impact on type 2 diabetes. To better understand the adipogenic process, the transcriptional events that occur during in vitro differentiation of human adipocytes were investigated and the findings linked to WAT phenotypes. Single-molecule transcriptional profiling provided a detailed map of the expressional changes of genes, enhancers, and long noncoding RNAs, where different types of transcripts share common dynamics during differentiation. Common signatures include early downregulated, transient, and late induced transcripts, all of which are linked to distinct developmental processes during adipogenesis. Enhancers expressed during adipogenesis overlap significantly with genetic variants associated with WAT distribution. Transiently expressed and late induced genes are associated with hypertrophic WAT (few but large fat cells), a phenotype closely linked to insulin resistance and type 2 diabetes. Transcription factors that are expressed early or transiently affect differentiation and adipocyte function and are controlled by several well-known upstream regulators such as glucocorticosteroids, insulin, cAMP, and thyroid hormones. Taken together, our results suggest a complex but highly coordinated regulation of adipogenesis.

Application of a CAGE Method to an Avian Development Study

Cap analysis of gene expression (CAGE) is a convenient approach for genome-wide identification of promoter regions at single base-pair resolution level and accurate expression estimation of the corresponding transcripts. Depending on the initial biomaterial amount and sequencing technology, different computational pipelines for data processing are available, as well as variations of the CAGE protocol that improve sensitivity and accuracy. Therefore, this chapter elucidates the key steps of sample preparation, sequencing, and data analysis via an example of a promoter expression estimation study in chicken development. We also describe the applicability of this approach for studying other avian and reptilian species.

The FANTOM5 Computation Ecosystem: Genomic Information Hub for Promoters and Active Enhancers

The Functional Annotation of the Mammalian Genome 5 (FANTOM5) project conducted transcriptome analysis of various mammalian cell types and provided a comprehensive resource to understand transcriptome and transcriptional regulation in individual cellular states encoded in the genome.FANTOM5 used cap analysis of gene expression (CAGE) with single-molecule sequencing to map transcription start sites (TSS) and measured their expression in a diverse range of samples. The main results from FANTOM5 were published as a promoter-level mammalian expression atlas and an atlas of active enhancers across human cell types. The FANTOM5 dataset is composed of raw experimental data and the results of bioinformatics analyses. In this chapter, we give a detailed description of the content of the FANTOM5 dataset and elaborate on different computing applications developed to publish the data and enable reproducibility and discovery of new findings. We present use cases in which the FANTOM5 dataset has been reused, leading to new findings.

On-the-fly selection of cell-specific enhancers, genes, miRNAs and proteins across the human body using SlideBase

Genomics consortia have produced large datasets profiling the expression of genes, micro-RNAs, enhancers and more across human tissues or cells. There is a need for intuitive tools to select subsets of such data that is the most relevant for specific studies. To this end, we present SlideBase, a web tool which offers a new way of selecting genes, promoters, enhancers and microRNAs that are preferentially expressed/used in a specified set of cells/tissues, based on the use of interactive sliders. With the help of sliders, SlideBase enables users to define custom expression thresholds for individual cell types/tissues, producing sets of genes, enhancers etc. which satisfy these constraints. Changes in slider settings result in simultaneous changes in the selected sets, updated in real time. SlideBase is linked to major databases from genomics consortia, including FANTOM, GTEx, The Human Protein Atlas and BioGPS.Database URL: http://slidebase.binf.ku.dk.

Genome-wide mapping of autonomous promoter activity in human cells

Previous methods to systematically characterize sequence-intrinsic activity of promoters have been limited by relatively low throughput and the length of the sequences that could be tested. Here we present 'survey of regulatory elements' (SuRE), a method that assays more than 10⁸ DNA fragments, each 0.2-2 kb in size, for their ability to drive transcription autonomously. In SuRE, a plasmid library of random genomic fragments upstream of a 20-bp barcode is constructed, and decoded by paired-end sequencing. This library is used to transfect cells, and barcodes in transcribed RNA are quantified by high-throughput sequencing. When applied to the human genome, we achieve 55-fold genome coverage, allowing us to map autonomous promoter activity genome-wide in K562 cells. By computational modeling we delineate subregions within promoters that are relevant for their activity. We show that antisense promoter transcription is generally dependent on the sense core promoter sequences, and that most enhancers and several families of repetitive elements act as autonomous transcription initiation sites.

DNA Hypomethylation in Intragenic and Intergenic Enhancer Chromatin of Muscle-Specific Genes Usually Correlates with their Expression

Tissue-specific enhancers are critical for gene regulation. In this study, we help elucidate the contribution of muscle-associated differential DNA methylation to the enhancer activity of highly muscle-specific genes. By bioinformatic analysis of 44 muscle-associated genes, we show that preferential gene expression in skeletal muscle (SkM) correlates with SkM-specific intragenic and intergenic enhancer chromatin and overlapping foci of DNA hypomethylation. Some genes, e.g., CASQ1 and FBXO32, displayed broad regions of both SkM- and heart-specific enhancer chromatin but exhibited focal SkM-specific DNA hypomethylation. Half of the genes had SkM-specific super-enhancers. In contrast to simple enhancer/gene-expression correlations, a super-enhancer was associated with the myogenic MYOD1 gene in both SkM and myoblasts even though SkM has < 1 percent as much MYOD1 expression. Local chromatin differences in this super-enhancer probably contribute to the SkM/myoblast differential expression. Transfection assays confirmed the tissue-specificity of the 0.3-kb core enhancer within MYOD1's super-enhancer and demonstrated its repression by methylation of its three CG dinucleotides. Our study suggests that DNA hypomethylation increases enhancer tissue-specificity and that SkM super-enhancers sometimes are poised for physiologically important, rapid up-regulation.