A high-resolution map of active promoters in the human genome

Dynamic RNA polymerase II occupancy drives differentiation of the intestine under the direction of HNF4

Terminal differentiation requires massive restructuring of the transcriptome. During intestinal differentiation, the expression patterns of nearly 4,000 genes are altered as cells transition from progenitor cells in crypts to differentiated cells in villi. We identify dynamic occupancy of RNA polymerase II (Pol II) to gene promoters as the primary driver of transcriptomic shifts during intestinal differentiation in vivo. Changes in enhancer-promoter looping interactions accompany dynamic Pol II occupancy and are dependent upon HNF4, a pro-differentiation transcription factor. Using genetic loss-of-function, chromatin immunoprecipitation sequencing (ChIP-seq), and immunoprecipitation (IP) mass spectrometry, we demonstrate that HNF4 collaborates with chromatin remodelers and loop-stabilizing proteins and facilitates Pol II occupancy at hundreds of genes pivotal to differentiation. We also explore alternate mechanisms that drive differentiation gene expression and find that pause-release of Pol II and post-transcriptional mRNA stability regulate smaller subsets of differentially expressed genes. These studies provide insights into the mechanisms of differentiation in renewing adult tissue.

Polymorphic transposable elements provide new insights on high-altitude adaptation in the Tibetan Plateau

Several studies demonstrated that populations living in the Tibetan plateau are genetically and physiologically adapted to high-altitude conditions, showing genomic signatures ascribable to the action of natural selection. However, so far most of them relied solely on inferences drawn from the analysis of coding variants and point mutations. To fill this gap, we focused on the possible role of polymorphic transposable elements in influencing the adaptation of Tibetan and Sherpa highlanders. To do so, we compared high-altitude and middle/low-lander individuals of East Asian ancestry by performing in silico analyses and differentiation tests on 118 modern and ancient samples. We detected several transposable elements associated with high altitude, which map genes involved in cardiovascular, hematological, chem-dependent and respiratory conditions, suggesting that metabolic and signaling pathways taking part in these functions are disproportionately impacted by the effect of environmental stressors in high-altitude individuals. To our knowledge, our study is the first hinting to a possible role of transposable elements in the adaptation of Tibetan and Sherpa highlanders.

Transcriptional elongation control of hypoxic response

The release of paused RNA polymerase II (RNAPII) from promoter-proximal regions is tightly controlled to ensure proper regulation of gene expression. The elongation factor PTEF-b is known to release paused RNAPII via phosphorylation of the RNAPII C-terminal domain by its cyclin-dependent kinase component, CDK9. However, the signal and stress-specific roles of the various RNAPII-associated macromolecular complexes containing PTEF-b/CDK9 are not yet clear. Here, we identify and characterize the CDK9 complex required for transcriptional response to hypoxia. Contrary to previous reports, our data indicate that a CDK9 complex containing BRD4 but not AFF1/4 is essential for this hypoxic stress response. We demonstrate that BRD4 bromodomains (BET) are dispensable for the release of paused RNAPII at hypoxia-activated genes and that BET inhibition by JQ1 is insufficient to impair hypoxic gene response. Mechanistically, we demonstrate that the C-terminal region of BRD4 is required for Polymerase-Associated Factor-1 Complex (PAF1C) recruitment to establish an elongation-competent RNAPII complex at hypoxia-responsive genes. PAF1C disruption using a small-molecule inhibitor (iPAF1C) impairs hypoxia-induced, BRD4-mediated RNAPII release. Together, our results provide insight into potentially targetable mechanisms that control the hypoxia-responsive transcriptional elongation.

Alcohol-Associated Liver Disease Outcomes: Critical Mechanisms of Liver Injury Progression

Alcohol-associated liver disease (ALD) is a substantial cause of morbidity and mortality worldwide and represents a spectrum of liver injury beginning with hepatic steatosis (fatty liver) progressing to inflammation and culminating in cirrhosis. Multiple factors contribute to ALD progression and disease severity. Here, we overview several crucial mechanisms related to ALD end-stage outcome development, such as epigenetic changes, cell death, hemolysis, hepatic stellate cells activation, and hepatic fatty acid binding protein 4. Additionally, in this review, we also present two clinically relevant models using human precision-cut liver slices and hepatic organoids to examine ALD pathogenesis and progression.

Differences in transcription initiation directionality underlie distinctions between plants and animals in chromatin modification patterns at genes and cis-regulatory elements

Transcriptional initiation is among the first regulated steps controlling eukaryotic gene expression. High-throughput profiling of fungal and animal genomes has revealed that RNA Polymerase II often initiates transcription in both directions at the promoter transcription start site, but generally only elongates productively into the gene body. Additionally, Pol II can initiate transcription in both directions at cis-regulatory elements such as enhancers. These bidirectional RNA Polymerase II initiation events can be observed directly with methods that capture nascent transcripts, and they are also revealed indirectly by the presence of transcription-associated histone modifications on both sides of the transcription start site or cis-regulatory elements. Previous studies have shown that nascent RNAs and transcription-associated histone modifications in the model plant Arabidopsis thaliana accumulate mainly in the gene body, suggesting that transcription does not initiate widely in the upstream direction from genes in this plant. We compared transcription-associated histone modifications and nascent transcripts at both transcription start sites and cis-regulatory elements in A. thaliana, Drosophila melanogaster, and Homo sapiens. Our results provide evidence for mostly unidirectional RNA Polymerase II initiation at both promoters and gene-proximal cis-regulatory elements of A. thaliana, whereas bidirectional transcription initiation is observed widely at promoters in both D. melanogaster and H. sapiens, as well as cis-regulatory elements in Drosophila. Furthermore, the distribution of transcription-associated histone modifications around transcription start sites in the Oryza sativa (rice) and Glycine max (soybean) genomes suggests that unidirectional transcription initiation is the norm in these genomes as well. These results suggest that there are fundamental differences in transcriptional initiation directionality between flowering plant and metazoan genomes, which are manifested as distinct patterns of chromatin modifications around RNA polymerase initiation sites.

Physical Peculiarity of Two Sites in Human Promoters: Universality and Diverse Usage in Gene Function

Since the discovery of physical peculiarities around transcription start sites (TSSs) and a site corresponding to the TATA box, research has revealed only the average features of these sites. Unsettled enigmas include the individual genes with these features and whether they relate to gene function. Herein, using 10 physical properties of DNA, including duplex DNA free energy, base stacking energy, protein-induced deformability, and stabilizing energy of Z-DNA, we clarified for the first time that approximately 97% of the promoters of 21,056 human protein-coding genes have distinctive physical properties around the TSS and/or position -27; of these, nearly 65% exhibited such properties at both sites. Furthermore, about 55% of the 21,056 genes had a minimum value of regional duplex DNA free energy within TSS-centered ±300 bp regions. Notably, distinctive physical properties within the promoters and free energies of the surrounding regions separated human protein-coding genes into five groups; each contained specific gene ontology (GO) terms. The group represented by immune response genes differed distinctly from the other four regarding the parameter of the free energies of the surrounding regions. A vital suggestion from this study is that physical-feature-based analyses of genomes may reveal new aspects of the organization and regulation of genes.

Biological Properties of Cyclitols and Their Derivatives

Cyclitols are polyhydroxy cycloalkanes, each containing at least three hydroxyls attached to a different ring carbon atom. The most important cyclitol derivatives are inositols, quercitols, conduritols and pinitols, which form a group of naturally occurring polyhydric alcohols and are widely found in plants. In addition, synthetic production of cyclitols has gained importance in recent years. Cylitols are molecules synthesized in plants as a precaution against salt or water stress. They have important functions in cell functioning as they exhibit important properties such as membrane biogenesis, ion channel physiology, signal transduction, osmoregulation, phosphate storage, cell wall formation and antioxidant activity. The biological activities of these very important molecules, obtained both synthetically and from the extraction of plants, are described in this review.

Identification and partial characterization of new cell density-dependent nucleocytoplasmic shuttling proteins and open chromatin

The contact inhibition of proliferation (CIP) denotes the cell density-dependent inhibition of growth, and the loss of CIP represents a hallmark of cancer. However, the mechanism by which CIP regulates gene expression remains poorly understood. Chromatin is a highly complex structure consisting of DNA, histones, and trans-acting factors (TAFs). The binding of TAF proteins to specific chromosomal loci regulates gene expression. Therefore, profiling chromatin is crucial for gaining insight into the gene expression mechanism of CIP. In this study, using modified proteomics of TAFs bound to DNA, we identified a protein that shuttles between the nucleus and cytosol in a cell density-dependent manner. We identified TIPARP, PTGES3, CBFB, and SMAD4 as cell density-dependent nucleocytoplasmic shuttling proteins. In low-density cells, these proteins predominantly reside in the nucleus; however, upon reaching high density, they relocate to the cytosol. Given their established roles in gene regulation, our findings propose their involvement as CIP-dependent TAFs. We also identified and characterized potential open chromatin regions sensitive to changes in cell density. These findings provide insights into the modulation of chromatin structure by CIP.

Model-based characterization of the equilibrium dynamics of transcription initiation and promoter-proximal pausing in human cells

In metazoans, both transcription initiation and the escape of RNA polymerase (RNAP) from promoter-proximal pausing are key rate-limiting steps in gene expression. These processes play out at physically proximal sites on the DNA template and appear to influence one another through steric interactions. Here, we examine the dynamics of these processes using a combination of statistical modeling, simulation, and analysis of real nascent RNA sequencing data. We develop a simple probabilistic model that jointly describes the kinetics of transcription initiation, pause-escape, and elongation, and the generation of nascent RNA sequencing read counts under steady-state conditions. We then extend this initial model to allow for variability across cells in promoter-proximal pause site locations and steric hindrance of transcription initiation from paused RNAPs. In an extensive series of simulations, we show that this model enables accurate estimation of initiation and pause-escape rates. Furthermore, we show by simulation and analysis of real data that pause-escape is often strongly rate-limiting and that steric hindrance can dramatically reduce initiation rates. Our modeling framework is applicable to a variety of inference problems, and our software for estimation and simulation is freely available.

Dynamic RNA Polymerase II Recruitment Drives Differentiation of the Intestine under the direction of HNF4

Terminal differentiation requires a massive restructuring of the transcriptome. During intestinal differentiation, the expression patterns of nearly 4000 genes are altered as cells transition from progenitor cells in crypts to differentiated cells in villi. We identified dynamic recruitment of RNA Polymerase II (Pol II) to gene promoters as the primary driver of transcriptomic shifts during intestinal differentiation in vivo. Changes in enhancer-promoter looping interactions accompany dynamic Pol II recruitment and are dependent upon HNF4, a pro-differentiation transcription factor. Using genetic loss-of- function, ChIP-seq and IP mass spectrometry, we demonstrate that HNF4 collaborates with chromatin remodelers and loop-stabilizing proteins and facilitates Pol II recruitment at hundreds of genes pivotal to differentiation. We also explore alternate mechanisms which drive differentiation gene expression and find pause-release of Pol II and post- transcriptional mRNA stability regulate smaller subsets of differentially expressed genes. These studies provide insights into the mechanisms of differentiation in a renewing adult tissue.

Differences in transcription initiation directionality underlie distinctions between plants and animals in chromatin modification patterns at genes and cis-regulatory elements

Transcriptional initiation is among the first regulated steps controlling eukaryotic gene expression. High-throughput profiling of fungal and animal genomes has revealed that RNA Polymerase II (Pol II) often initiates transcription in both directions at the promoter transcription start site (TSS), but generally only elongates productively into the gene body. Additionally, Pol II can initiate transcription in both directions at cis-regulatory elements (CREs) such as enhancers. These bidirectional Pol II initiation events can be observed directly with methods that capture nascent transcripts, and they are also revealed indirectly by the presence of transcription-associated histone modifications on both sides of the TSS or CRE. Previous studies have shown that nascent RNAs and transcription-associated histone modifications in the model plant Arabidopsis thaliana accumulate mainly in the gene body, suggesting that transcription does not initiate widely in the upstream direction from genes in this plant. We compared transcription-associated histone modifications and nascent transcripts at both TSSs and CREs in Arabidopsis thaliana, Drosophila melanogaster, and Homo sapiens. Our results provide evidence for mostly unidirectional Pol II initiation at both promoters and gene-proximal CREs of Arabidopsis thaliana, whereas bidirectional transcription initiation is observed widely at promoters in both Drosophila melanogaster and Homo sapiens, as well as CREs in Drosophila. Furthermore, the distribution of transcription-associated histone modifications around TSSs in the Oryza sativa (rice) and Glycine max (soybean) genomes suggests that unidirectional transcription initiation is the norm in these genomes as well. These results suggest that there are fundamental differences in transcriptional initiation directionality between flowering plant and metazoan genomes, which are manifested as distinct patterns of chromatin modifications around RNA polymerase initiation sites.

SEM: sized-based expectation maximization for characterizing nucleosome positions and subtypes

Genome-wide nucleosome profiles are predominantly characterized using MNase-seq, which involves extensive MNase digestion and size selection to enrich for mono-nucleosome-sized fragments. Most available MNase-seq analysis packages assume that nucleosomes uniformly protect 147bp DNA fragments. However, some nucleosomes with atypical histone or chemical compositions protect shorter lengths of DNA. The rigid assumptions imposed by current nucleosome analysis packages ignore variation in nucleosome lengths, potentially blinding investigators to regulatory roles played by atypical nucleosomes. To enable the characterization of different nucleosome types from MNase-seq data, we introduce the Size-based Expectation Maximization (SEM) nucleosome calling package. SEM employs a hierarchical Gaussian mixture model to estimate the positions and subtype identity of nucleosomes from MNase-seq fragments. Nucleosome subtypes are automatically identified based on the distribution of protected DNA fragment lengths at nucleosome positions. Benchmark analysis indicates that SEM is on par with existing packages in terms of standard nucleosome-calling accuracy metrics, while uniquely providing the ability to characterize nucleosome subtype identities. Using SEM on a low-dose MNase H2B MNase-ChIP-seq dataset from mouse embryonic stem cells, we identified three nucleosome types: short-fragment nucleosomes, canonical nucleosomes, and di-nucleosomes. The short-fragment nucleosomes can be divided further into two subtypes based on their chromatin accessibility. Interestingly, the subset of short-fragment nucleosomes in accessible regions exhibit high MNase sensitivity and display distribution patterns around transcription start sites (TSSs) and CTCF peaks, similar to the previously reported "fragile nucleosomes". These SEM-defined accessible short-fragment nucleosomes are found not just in promoters, but also in enhancers and other regulatory regions. Additional investigations reveal their co-localization with the chromatin remodelers Chd6, Chd8, and Ep400. In summary, SEM provides an effective platform for distinguishing various nucleosome subtypes, paving the way for future exploration of non-standard nucleosomes.

A multilayered post-GWAS analysis pipeline defines functional variants and target genes for systemic lupus erythematosus (SLE)

Objectives

Systemic lupus erythematosus (SLE), an autoimmune disease with incompletely understood etiology, has a strong genetic component. Although genome-wide association studies (GWAS) have revealed multiple SLE susceptibility loci and associated single nucleotide polymorphisms (SNPs), the precise causal variants, target genes, cell types, tissues, and mechanisms of action remain largely unknown.

Methods

Here, we report a comprehensive post-GWAS analysis using extensive bioinformatics, molecular modeling, and integrative functional genomic and epigenomic analyses to optimize fine-mapping. We compile and cross-reference immune cell-specific expression quantitative trait loci ( cis - and trans -eQTLs) with promoter-capture Hi-C, allele-specific chromatin accessibility, and massively parallel reporter assay data to define predisposing variants and target genes. We experimentally validate a predicted locus using CRISPR/Cas9 genome editing, qPCR, and Western blot.

Results

Anchoring on 452 index SNPs, we selected 9,931 high-linkage disequilibrium (r 2 >0.8) SNPs and defined 182 independent non-HLA SLE loci. 3,746 SNPs from 143 loci were identified as regulating 564 unique genes. Target genes are enriched in lupus-related tissues and associated with other autoimmune diseases. Of these, 329 SNPs (106 loci) showed significant allele-specific chromatin accessibility and/or enhancer activity, indicating regulatory potential. Using CRISPR/Cas9, we validated rs57668933 as a functional variant regulating multiple targets, including SLE risk gene ELF1 , in B-cells.

Conclusion

We demonstrate and validate post-GWAS strategies for utilizing multi-dimensional data to prioritize likely causal variants with cognate gene targets underlying SLE pathogenesis. Our results provide a catalog of significantly SLE-associated SNPs and loci, target genes, and likely biochemical mechanisms, to guide experimental characterization.

FOXG1 targets BMP repressors and cell cycle inhibitors in human neural progenitor cells

FOXG1 is a critical transcription factor in human brain where loss-of-function mutations cause a severe neurodevelopmental disorder, while increased FOXG1 expression is frequently observed in glioblastoma. FOXG1 is an inhibitor of cell patterning and an activator of cell proliferation in chordate model organisms but different mechanisms have been proposed as to how this occurs. To identify genomic targets of FOXG1 in human neural progenitor cells (NPCs), we engineered a cleavable reporter construct in endogenous FOXG1 and performed chromatin immunoprecipitation (ChIP) sequencing. We also performed deep RNA sequencing of NPCs from two females with loss-of-function mutations in FOXG1 and their healthy biological mothers. Integrative analyses of RNA and ChIP sequencing data showed that cell cycle regulation and Bone Morphogenic Protein (BMP) repression gene ontology categories were over-represented as FOXG1 targets. Using engineered brain cell lines, we show that FOXG1 specifically activates SMAD7 and represses CDKN1B. Activation of SMAD7 which inhibits BMP signaling may be one way that FOXG1 patterns the forebrain, while repression of cell cycle regulators such as CDKN1B may be one way that FOXG1 expands the NPC pool to ensure proper brain size. Our data reveal novel mechanisms on how FOXG1 may control forebrain patterning and cell proliferation in human brain development.

Characterization of Gene Regulatory Elements in Human Fetal Cortical Development: Enhancing Our Understanding of Neurodevelopmental Disorders and Evolution

The neocortex is the region that most distinguishes human brain from other mammals and primates [Annu Rev Genet. 2021 Nov;55(1):555-81]. Studying the development of human cortex is important in understanding the evolutionary changes occurring in humans relative to other primates, as well as in elucidating mechanisms underlying neurodevelopmental disorders. Cortical development is a highly regulated process, spatially and temporally coordinated by expression of essential transcriptional factors in response to signaling pathways [Neuron. 2019 Sep;103(6):980-1004]. Enhancers are the most well-understood cis-acting, non-protein-coding regulatory elements that regulate gene expression [Nat Rev Genet. 2014 Apr;15(4):272-86]. Importantly, given the conservation of both DNA sequence and molecular function of the majority of proteins across mammals [Genome Res. 2003 Dec;13(12):2507-18], enhancers [Science. 2015 Mar;347(6226):1155-9], which are far more divergent at the sequence level, likely account for the phenotypes that distinguish the human brain by changing the regulation of gene expression. In this review, we will revisit the conceptual framework of gene regulation during human brain development, as well as the evolution of technologies to study transcriptional regulation, with recent advances in genome biology that open a window allowing us to systematically characterize cis-regulatory elements in developing human brain [Hum Mol Genet. 2022 Oct;31(R1):R84-96]. We provide an update on work to characterize the suite of all enhancers in the developing human brain and the implications for understanding neuropsychiatric disorders. Finally, we discuss emerging therapeutic ideas that utilize our emerging knowledge of enhancer function.

Genome-Wide Identification of Open Chromatin in Plants Using MH-Seq

Functional cis-regulatory elements (CREs) act as precise transcriptional switches for fine-tuning gene transcription. Identification of CREs is critical for understanding regulatory mechanisms of gene expression associated with various biological processes in eukaryotes. It is well known that CREs reside in open chromatin that exhibits hypersensitivity to enzyme cleavage and physical shearing. Currently, high-throughput methodologies, such as DNase-seq, ATAC-seq, and FAIRE-seq, have been widely applied in mapping open chromatin in various eukaryotic genomes. More recently, differential MNase (micrococcal nuclease) treatment has been successfully employed to map open chromatin in addition to profiling nucleosome landscape in both mammalian and plant species. We have developed a MNase hypersensitivity sequencing (MH-seq) technique in plants. The MH-seq procedure includes plant nuclei fixation and purification, differential treatments of purified nuclei with MNase, specific recovery of MNase-trimmed small DNA fragments within 20~100 bp in length, and MH-seq library construction followed by Illumina sequencing and data analysis. MH-seq has been successfully applied for global identification of open chromatin in both Arabidopsis thaliana and maize. It has been proven to be an attractive alternative for profiling open chromatin. Thus, MH-seq is expected to be valuable in probing chromatin accessibility on a genome-wide scale for other plants with sequenced genomes. Moreover, MHS data allow to implement footprinting assays to unveil binding sites of transcription factors.

Shaking up the silence: consequences of HMGN1 antagonizing PRC2 in the Down syndrome brain

Intellectual disability is a well-known hallmark of Down Syndrome (DS) that results from the triplication of the critical region of human chromosome 21 (HSA21). Major studies were conducted in recent years to gain an understanding about the contribution of individual triplicated genes to DS-related brain pathology. Global transcriptomic alterations and widespread changes in the establishment of neural lineages, as well as their differentiation and functional maturity, suggest genome-wide chromatin organization alterations in trisomy. High Mobility Group Nucleosome Binding Domain 1 (HMGN1), expressed from HSA21, is a chromatin remodeling protein that facilitates chromatin decompaction and is associated with acetylated lysine 27 on histone H3 (H3K27ac), a mark correlated with active transcription. Recent studies causatively linked overexpression of HMGN1 in trisomy and the development of DS-associated B cell acute lymphoblastic leukemia (B-ALL). HMGN1 has been shown to antagonize the activity of the Polycomb Repressive Complex 2 (PRC2) and prevent the deposition of histone H3 lysine 27 trimethylation mark (H3K27me3), which is associated with transcriptional repression and gene silencing. However, the possible ramifications of the increased levels of HMGN1 through the derepression of PRC2 target genes on brain cell pathology have not gained attention. In this review, we discuss the functional significance of HMGN1 in brain development and summarize accumulating reports about the essential role of PRC2 in the development of the neural system. Mechanistic understanding of how overexpression of HMGN1 may contribute to aberrant brain cell phenotypes in DS, such as altered proliferation of neural progenitors, abnormal cortical architecture, diminished myelination, neurodegeneration, and Alzheimer's disease-related pathology in trisomy 21, will facilitate the development of DS therapeutic approaches targeting chromatin.

Multi-omics peripheral and core regions of cancer

Thousands of genes are perturbed by cancer, and these disturbances can be seen in transcriptome, methylation, somatic mutation, and copy number variation omics studies. Understanding their connectivity patterns as an omnigenic neighbourhood in a molecular interaction network (interactome) is a key step towards advancing knowledge of the molecular mechanisms underlying cancers. Here, we introduce a unified connectivity line (CLine) to pinpoint omics-specific omnigenic patterns across 15 curated cancers. Taking advantage of the universality of CLine, we distinguish the peripheral and core genes for each omics aspect. We propose a network-based framework, multi-omics periphery and core (MOPC), to combine peripheral and core genes from different omics into a button-like structure. On the basis of network proximity, we provide evidence that core genes tend to be specifically perturbed in one omics, but the peripheral genes are diversely perturbed in multiple omics. And the core of one omics is regulated by multiple omics peripheries. Finally, we take the MOPC as an omnigenic neighbourhood, describe its characteristics, and explore its relative contribution to network-based mechanisms of cancer. We were able to present how multi-omics perturbations percolate through the human interactome and contribute to an integrated periphery and core.

Palmitate Potentiates Lipopolysaccharide-Induced IL-6 Production via Coordinated Acetylation of H3K9/H3K18, p300, and RNA Polymerase II

IL-6 is elevated in obese individuals and participates in the metabolic dysfunction associated with that condition. However, the mechanisms that promote IL-6 expression in obesity are incompletely understood. Because elevated levels of palmitate and LPS have been reported in obesity, we investigated whether these agents interact to potentiate IL-6 production. In this study, we report that LPS induces higher levels of IL-6 in human monocytes in the presence of palmitate. Notably, the priming effect of palmitate is associated with enhanced p300 binding and transcription factor recruitment to Il6 promoter regions. Gene silencing of p300 blocks this action of palmitate. RNA polymerase II recruitment was also enhanced at the Il6 promoter in palmitate/LPS-exposed cells. Acetylation levels of H3K9 and H3K18 were increased in monocytes treated with palmitate. Moreover, LPS stimulation of palmitate-treated cells led to increased levels of the transcriptionally permissive acetylation marks H3K9/H3K18 in the Il6 promoter compared with LPS alone. The effect of palmitate on LPS-induced IL-6 production was suppressed by the inhibition of histone acetyltransferases. Conversely, histone deacetylase inhibitors trichostatin A or sodium butyrate can substitute for palmitate in IL-6 production. Esterification of palmitate with CoA was involved, whereas β-oxidation and ceramide biosynthesis were not required, for the induction of IL-6 and H3K9/H3K18 acetylation. Monocytes of obese individuals showed significantly higher H3K9/H3K18 acetylation and Il6 expression. Overall, our findings support a model in which increased levels of palmitate in obesity create a setting for LPS to potentiate IL-6 production via chromatin remodeling, enabling palmitate to contribute to metabolic inflammation.

Conserved Control Path in Multilayer Networks

The determination of directed control paths in complex networks is important because control paths indicate the structure of the propagation of control signals through edges. A challenging problem is to identify them in complex networked systems characterized by different types of interactions that form multilayer networks. In this study, we describe a graph pattern called the conserved control path, which allows us to model a common control structure among different types of relations. We present a practical conserved control path detection method (CoPath), which is based on a maximum-weighted matching, to determine the paths that play the most consistent roles in controlling signal transmission in multilayer networks. As a pragmatic application, we demonstrate that the control paths detected in a multilayered pan-cancer network are statistically more consistent. Additionally, they lead to the effective identification of drug targets, thereby demonstrating their power in predicting key pathways that influence multiple cancers.

GoPeaks: histone modification peak calling for CUT&Tag

Genome-wide mapping of histone modifications is critical to understanding transcriptional regulation. CUT&Tag is a new method for profiling histone modifications, offering improved sensitivity and decreased cost compared with ChIP-seq. Here, we present GoPeaks, a peak calling method specifically designed for histone modification CUT&Tag data. We compare the performance of GoPeaks against commonly used peak calling algorithms to detect histone modifications that display a range of peak profiles and are frequently used in epigenetic studies. We find that GoPeaks robustly detects genome-wide histone modifications and, notably, identifies a substantial number of H3K27ac peaks with improved sensitivity compared to other standard algorithms.

The chromatin remodeler Ino80 mediates RNAPII pausing site determination

BACKGROUND

Promoter-proximal pausing of RNA polymerase II (RNAPII) is a critical step for the precise regulation of gene expression. Despite the apparent close relationship between promoter-proximal pausing and nucleosome, the role of chromatin remodeler governing this step has mainly remained elusive.

RESULTS

Here, we report highly confined RNAPII enrichments downstream of the transcriptional start site in Saccharomyces cerevisiae using PRO-seq experiments. This non-uniform distribution of RNAPII exhibits both similar and different characteristics with promoter-proximal pausing in Schizosaccharomyces pombe and metazoans. Interestingly, we find that Ino80p knockdown causes a significant upstream transition of promoter-proximal RNAPII for a subset of genes, relocating RNAPII from the main pausing site to the alternative pausing site. The proper positioning of RNAPII is largely dependent on nucleosome context. We reveal that the alternative pausing site is closely associated with the + 1 nucleosome, and nucleosome architecture around the main pausing site of these genes is highly phased. In addition, Ino80p knockdown results in an increase in fuzziness and a decrease in stability of the + 1 nucleosome. Furthermore, the loss of INO80 also leads to the shift of promoter-proximal RNAPII toward the alternative pausing site in mouse embryonic stem cells.

CONCLUSIONS

Based on our collective results, we hypothesize that the highly conserved chromatin remodeler Ino80p is essential in establishing intact RNAPII pausing during early transcription elongation in various organisms, from budding yeast to mouse.

Distinguishing between models of mammalian gene expression: telegraph-like models versus mechanistic models

Two-state models (telegraph-like models) have a successful history of predicting distributions of cellular and nascent mRNA numbers that can well fit experimental data. These models exclude key rate limiting steps, and hence it is unclear why they are able to accurately predict the number distributions. To answer this question, here we compare these models to a novel stochastic mechanistic model of transcription in mammalian cells that presents a unified description of transcriptional factor, polymerase and mature mRNA dynamics. We show that there is a large region of parameter space where the first, second and third moments of the distributions of the waiting times between two consecutively produced transcripts (nascent or mature) of two-state and mechanistic models exactly match. In this region: (i) one can uniquely express the two-state model parameters in terms of those of the mechanistic model, (ii) the models are practically indistinguishable by comparison of their transcript numbers distributions, and (iii) they are distinguishable from the shape of their waiting time distributions. Our results clarify the relationship between different gene expression models and identify a means to select between them from experimental data.

Targeting the Transcriptome Through Globally Acting Components

Transcription is a step in gene expression that defines the identity of cells and its dysregulation is associated with diseases. With advancing technologies revealing molecular underpinnings of the cell with ever-higher precision, our ability to view the transcriptomes may have surpassed our knowledge of the principles behind their organization. The human RNA polymerase II (Pol II) machinery comprises thousands of components that, in conjunction with epigenetic and other mechanisms, drive specialized programs of development, differentiation, and responses to the environment. Parts of these programs are repurposed in oncogenic transformation. Targeting of cancers is commonly done by inhibiting general or broadly acting components of the cellular machinery. The critical unanswered question is how globally acting or general factors exert cell type specific effects on transcription. One solution, which is discussed here, may be among the events that take place at genes during early Pol II transcription elongation. This essay turns the spotlight on the well-known phenomenon of promoter-proximal Pol II pausing as a step that separates signals that establish pausing genome-wide from those that release the paused Pol II into the gene. Concepts generated in this rapidly developing field will enhance our understanding of basic principles behind transcriptome organization and hopefully translate into better therapies at the bedside.

The Core Promoter Is a Regulatory Hub for Developmental Gene Expression

The development of multicellular organisms and the uniqueness of each cell are achieved by distinct transcriptional programs. Multiple processes that regulate gene expression converge at the core promoter region, an 80 bp region that directs accurate transcription initiation by RNA polymerase II (Pol II). In recent years, it has become apparent that the core promoter region is not a passive DNA component, but rather an active regulatory module of transcriptional programs. Distinct core promoter compositions were demonstrated to result in different transcriptional outputs. In this mini-review, we focus on the role of the core promoter, particularly its downstream region, as the regulatory hub for developmental genes. The downstream core promoter element (DPE) was implicated in the control of evolutionarily conserved developmental gene regulatory networks (GRNs) governing body plan in both the anterior-posterior and dorsal-ventral axes. Notably, the composition of the basal transcription machinery is not universal, but rather promoter-dependent, highlighting the importance of specialized transcription complexes and their core promoter target sequences as key hubs that drive embryonic development, differentiation and morphogenesis across metazoan species. The extent of transcriptional activation by a specific enhancer is dependent on its compatibility with the relevant core promoter. The core promoter content also regulates transcription burst size. Overall, while for many years it was thought that the specificity of gene expression is primarily determined by enhancers, it is now clear that the core promoter region comprises an important regulatory module in the intricate networks of developmental gene expression.

Identification and Characterisation of cis-Regulatory Elements Upstream of the Human Receptor Tyrosine Kinase Gene MERTK

BACKGROUND:

MERTK encodes a receptor tyrosine kinase that regulates immune homeostasis via phagocytosis of apoptotic cells and cytokine-mediated immunosuppression. MERTK is highly expressed in the central nervous system (CNS), specifically in myeloid derived innate immune cells and its dysregulation is implicated in CNS pathologies including the autoimmune disease multiple sclerosis (MS).

OBJECTIVE:

While the cell types and tissues that express MERTK have been well described, the genetic elements that define the gene’s promoter and regulate specific transcription domains remain unknown. The primary objective of this study was to define and characterise the human MERTK promoter region.

METHODS:

We cloned and characterized the 5’ upstream region of MERTK to identify cis-acting DNA elements that promote gene transcription in luciferase reporter assays. In addition, promoter regions were tested for sensitivity to the anti-inflammatory glucocorticoid dexamethasone.

RESULTS:

This study identified identified both proximal and distal-acting DNA elements that promote transcription. The strongest promoter activity was identified in an ∼850 bp region situated 3 kb upstream of the MERTK transcription start site. Serial deletions of this putative enhancer revealed that the entire region is essential for expression activity. Using in silico analysis, we identified several candidate transcription factor binding sites. Despite a well-established upregulation of MERTK in response to anti-inflammatory glucocorticoids, no DNA region within the 5 kb putative promoter was found to directly respond to dexamethasone treatment.

CONCLUSIONS:

Elucidating the genetic mechanisms that regulate MERTK expression gives insights into gene regulation during homeostasis and disease, providing potential targets for therapeutic modulation of MERTK transcription.

Computational identification and experimental characterization of preferred downstream positions in human core promoters

Metazoan core promoters, which direct the initiation of transcription by RNA polymerase II (Pol II), may contain short sequence motifs termed core promoter elements/motifs (e.g. the TATA box, initiator (Inr) and downstream core promoter element (DPE)), which recruit Pol II via the general transcription machinery. The DPE was discovered and extensively characterized in Drosophila, where it is strictly dependent on both the presence of an Inr and the precise spacing from it. Since the Drosophila DPE is recognized by the human transcription machinery, it is most likely that some human promoters contain a downstream element that is similar, though not necessarily identical, to the Drosophila DPE. However, only a couple of human promoters were shown to contain a functional DPE, and attempts to computationally detect human DPE-containing promoters have mostly been unsuccessful. Using a newly-designed motif discovery strategy based on Expectation-Maximization probabilistic partitioning algorithms, we discovered preferred downstream positions (PDP) in human promoters that resemble the Drosophila DPE. Available chromatin accessibility footprints revealed that Drosophila and human Inr+DPE promoter classes are not only highly structured, but also similar to each other, particularly in the proximal downstream region. Clustering of the corresponding sequence motifs using a neighbor-joining algorithm strongly suggests that canonical Inr+DPE promoters could be common to metazoan species. Using reporter assays we demonstrate the contribution of the identified downstream positions to the function of multiple human promoters. Furthermore, we show that alteration of the spacing between the Inr and PDP by two nucleotides results in reduced promoter activity, suggesting a spacing dependency of the newly discovered human PDP on the Inr. Taken together, our strategy identified novel functional downstream positions within human core promoters, supporting the existence of DPE-like motifs in human promoters.

Transcription of genes by the RNA polymerase II enzyme initiates at a genomic region termed the core promoter. The core promoter is a regulatory region that may contain diverse short DNA sequence motifs/elements that confer specific properties to it. Interestingly, core promoter motifs can be located both upstream and downstream of the transcription start site. Variable compositions of core promoter elements were identified. The initiator (Inr) motif and the downstream core promoter element (DPE) is a combination of elements that has been identified and extensively characterized in fruit flies. Although a few Inr+DPE -containing human promoters were identified, the presence of transcriptionally important downstream core promoter positions within human promoters has been a matter of controversy in the literature. Here, using a newly-designed motif discovery strategy, we discovered preferred downstream positions in human promoters that resemble fruit fly DPE. Clustering of the corresponding sequence motifs in eight additional species indicated that such promoters could be common to multicellular non-plant organisms. Importantly, functional characterization of the newly discovered preferred downstream positions supports the existence of Inr+DPE-containing promoters in human genes.

Comparison of Histone H3K4me3 between IVF and ICSI Technologies and between Boy and Girl Offspring

Epigenetics play a vital role in early embryo development. Offspring conceived via assisted reproductive technologies (ARTs) have a three times higher risk of epigenetic diseases than naturally conceived children. However, investigations into ART-associated placental histone modifications or sex-stratified analyses of ART-associated histone modifications remain limited. In the current study, we carried out immunohistochemistry, chip-sequence analysis, and a series of in vitro experiments. Our results demonstrated that placentas from intra-cytoplasmic sperm injection (ICSI), but not in vitro fertilization (IVF), showed global tri-methylated-histone-H3-lysine-4 (H3K4me3) alteration compared to those from natural conception. However, for acetylated-histone-H3-lysine-9 (H3K9ac) and acetylated-histone-H3-lysine-27 (H3K27ac), no significant differences between groups could be found. Further, sex -stratified analysis found that, compared with the same-gender newborn cord blood mononuclear cell (CBMC) from natural conceptions, CBMC from ICSI-boys presented more genes with differentially enriched H3K4me3 (n = 198) than those from ICSI-girls (n = 79), IVF-girls (n = 5), and IVF-boys (n = 2). We also found that varying oxygen conditions, RNA polymerase II subunit A (Polr2A), and lysine demethylase 5A (KDM5A) regulated H3K4me3. These findings revealed a difference between IVF and ICSI and a difference between boys and girls in H3K4me3 modification, providing greater insight into ART-associated epigenetic alteration.

A gain-of-function single nucleotide variant creates a new promoter which acts as an orientation-dependent enhancer-blocker

Many single nucleotide variants (SNVs) associated with human traits and genetic diseases are thought to alter the activity of existing regulatory elements. Some SNVs may also create entirely new regulatory elements which change gene expression, but the mechanism by which they do so is largely unknown. Here we show that a single base change in an otherwise unremarkable region of the human α-globin cluster creates an entirely new promoter and an associated unidirectional transcript. This SNV downregulates α-globin expression causing α-thalassaemia. Of note, the new promoter lying between the α-globin genes and their associated super-enhancer disrupts their interaction in an orientation-dependent manner. Together these observations show how both the order and orientation of the fundamental elements of the genome determine patterns of gene expression and support the concept that active genes may act to disrupt enhancer-promoter interactions in mammals as in Drosophila. Finally, these findings should prompt others to fully evaluate SNVs lying outside of known regulatory elements as causing changes in gene expression by creating new regulatory elements.

The role of promoters as potential insulator elements has been largely unexplored in mammals. Here the authors show that a single nucleotide variant in the α-globin locus forms a new promoter and acts as an orientation-dependent enhancer-blocking insulator element.

CDK9 keeps RNA polymerase II on track

Cyclin-dependent kinase 9 (CDK9), the kinase component of positive transcription elongation factor b (P-TEFb), is essential for transcription of most protein-coding genes by RNA polymerase II (RNAPII). By releasing promoter-proximally paused RNAPII into gene bodies, CDK9 controls the entry of RNAPII into productive elongation and is, therefore, critical for efficient synthesis of full-length messenger (m)RNAs. In recent years, new players involved in P-TEFb-dependent processes have been identified and an important function of CDK9 in coordinating elongation with transcription initiation and termination has been unveiled. As the regulatory functions of CDK9 in gene expression continue to expand, a number of human pathologies, including cancers, have been associated with aberrant CDK9 activity, underscoring the need to properly regulate CDK9. Here, I provide an overview of CDK9 function and regulation, with an emphasis on CDK9 dysregulation in human diseases.

Promoter switching in response to changing environment and elevated expression of protein-coding genes overlapping at their 5' ends

Despite the number of studies focused on sense-antisense transcription, the key question of whether such organization evolved as a regulator of gene expression or if this is only a byproduct of other regulatory processes has not been elucidated to date. In this study, protein-coding sense-antisense gene pairs were analyzed with a particular focus on pairs overlapping at their 5’ ends. Analyses were performed in 73 human transcription start site libraries. The results of our studies showed that the overlap between genes is not a stable feature and depends on which TSSs are utilized in a given cell type. An analysis of gene expression did not confirm that overlap between genes causes downregulation of their expression. This observation contradicts earlier findings. In addition, we showed that the switch from one promoter to another, leading to genes overlap, may occur in response to changing environment of a cell or tissue. We also demonstrated that in transfected and cancerous cells genes overlap is observed more often in comparison with normal tissues. Moreover, utilization of overlapping promoters depends on particular state of a cell and, at least in some groups of genes, is not merely coincidental.

Deletion of a non-canonical regulatory sequence causes loss of Scn1a expression and epileptic phenotypes in mice

BACKGROUND

Genes with multiple co-active promoters appear common in brain, yet little is known about functional requirements for these potentially redundant genomic regulatory elements. SCN1A, which encodes the NaV1.1 sodium channel alpha subunit, is one such gene with two co-active promoters. Mutations in SCN1A are associated with epilepsy, including Dravet syndrome (DS). The majority of DS patients harbor coding mutations causing SCN1A haploinsufficiency; however, putative causal non-coding promoter mutations have been identified.

METHODS

To determine the functional role of one of these potentially redundant Scn1a promoters, we focused on the non-coding Scn1a 1b regulatory region, previously described as a non-canonical alternative transcriptional start site. We generated a transgenic mouse line with deletion of the extended evolutionarily conserved 1b non-coding interval and characterized changes in gene and protein expression, and assessed seizure activity and alterations in behavior.

RESULTS

Mice harboring a deletion of the 1b non-coding interval exhibited surprisingly severe reductions of Scn1a and NaV1.1 expression throughout the brain. This was accompanied by electroencephalographic and thermal-evoked seizures, and behavioral deficits.

CONCLUSIONS

This work contributes to functional dissection of the regulatory wiring of a major epilepsy risk gene, SCN1A. We identified the 1b region as a critical disease-relevant regulatory element and provide evidence that non-canonical and seemingly redundant promoters can have essential function.

The Transposable Element Environment of Human Genes Differs According to Their Duplication Status and Essentiality

Transposable elements (TEs) are major components of eukaryotic genomes and represent approximately 45% of the human genome. TEs can be important sources of novelty in genomes and there is increasing evidence that TEs contribute to the evolution of gene regulation in mammals. Gene duplication is an evolutionary mechanism that also provides new genetic material and opportunities to acquire new functions. To investigate how duplicated genes are maintained in genomes, here, we explored the TE environment of duplicated and singleton genes. We found that singleton genes have more short-interspersed nuclear elements and DNA transposons in their vicinity than duplicated genes, whereas long-interspersed nuclear elements and long-terminal repeat retrotransposons have accumulated more near duplicated genes. We also discovered that this result is highly associated with the degree of essentiality of the genes with an unexpected accumulation of short-interspersed nuclear elements and DNA transposons around the more-essential genes. Our results underline the importance of taking into account the TE environment of genes to better understand how duplicated genes are maintained in genomes.

Profiling single-cell histone modifications using indexing chromatin immunocleavage sequencing

Recently, multiple single-cell assays were developed for detecting histone marks at the single-cell level. These techniques are either limited by the low cell throughput or sparse reads which limit their applications. To address these problems, we introduce indexing single-cell immunocleavage sequencing (iscChIC-seq), a multiplex indexing method based on TdT terminal transferase and T4 DNA ligase-mediated barcoding strategy and single-cell ChIC-seq, which is capable of readily analyzing histone modifications across tens of thousands of single cells in one experiment. Application of iscChIC-seq to profiling H3K4me3 and H3K27me3 in human white blood cells (WBCs) enabled successful detection of more than 10,000 single cells for each histone modification with 11 K and 45 K nonredundant reads per cell, respectively. Cluster analysis of these data allowed identification of monocytes, T cells, B cells, and NK cells from WBCs. The cell types annotated from H3K4me3 single-cell data are specifically correlated with the cell types annotated from H3K27me3 single-cell data. Our data indicate that iscChIC-seq is a reliable technique for profiling histone modifications in a large number of single cells, which may find broad applications in studying cellular heterogeneity and differentiation status in complex developmental and disease systems.

Topoisomerase IIα represses transcription by enforcing promoter-proximal pausing

Accumulation of topological stress in the form of DNA supercoiling is inherent to the advance of RNA polymerase II (Pol II) and needs to be resolved by DNA topoisomerases to sustain productive transcriptional elongation. Topoisomerases are therefore considered positive facilitators of transcription. Here, we show that, in contrast to this general assumption, human topoisomerase IIα (TOP2A) activity at promoters represses transcription of immediate early genes such as c-FOS, maintaining them under basal repressed conditions. Thus, TOP2A inhibition creates a particular topological context that results in rapid release from promoter-proximal pausing and transcriptional upregulation, which mimics the typical bursting behavior of these genes in response to physiological stimulus. We therefore describe the control of promoter-proximal pausing by TOP2A as a layer for the regulation of gene expression, which can act as a molecular switch to rapidly activate transcription, possibly by regulating the accumulation of DNA supercoiling at promoter regions.

Understanding transcription across scales: From base pairs to chromosomes

The influence of genome organization on transcription is central to our understanding of cell type specification. Higher-order genome organization is established through short- and long-range DNA interactions. Coordination of these interactions, from single atoms to entire chromosomes, plays a fundamental role in transcriptional control of gene expression. Loss of this coupling can result in disease. Analysis of transcriptional regulation typically involves disparate experimental approaches, from structural studies that define angstrom-level interactions to cell-biological and genomic approaches that assess mesoscale relationships. Thus, to fully understand the mechanisms that regulate gene expression, it is critical to integrate the findings gained across these distinct size scales. In this review, I illustrate fundamental ways in which cells regulate transcription in the context of genome organization.

Enhancers predominantly regulate gene expression during differentiation via transcription initiation

Gene transcription occurs via a cycle of linked events, including initiation, promoter-proximal pausing, and elongation of RNA polymerase II (Pol II). A key question is how transcriptional enhancers influence these events to control gene expression. Here, we present an approach that evaluates the level and change in promoter-proximal transcription (initiation and pausing) in the context of differential gene expression, genome-wide. This combinatorial approach shows that in primary cells, control of gene expression during differentiation is achieved predominantly via changes in transcription initiation rather than via release of Pol II pausing. Using genetically engineered mouse models, deleted for functionally validated enhancers of the α- and β-globin loci, we confirm that these elements regulate Pol II recruitment and/or initiation to modulate gene expression. Together, our data show that gene expression during differentiation is regulated predominantly at the level of initiation and that enhancers are key effectors of this process.

LINE retrotransposons characterize mammalian tissue-specific and evolutionarily dynamic regulatory regions

BACKGROUND

To investigate the mechanisms driving regulatory evolution across tissues, we experimentally mapped promoters, enhancers, and gene expression in the liver, brain, muscle, and testis from ten diverse mammals.

RESULTS

The regulatory landscape around genes included both tissue-shared and tissue-specific regulatory regions, where tissue-specific promoters and enhancers evolved most rapidly. Genomic regions switching between promoters and enhancers were more common across species, and less common across tissues within a single species. Long Interspersed Nuclear Elements (LINEs) played recurrent evolutionary roles: LINE L1s were associated with tissue-specific regulatory regions, whereas more ancient LINE L2s were associated with tissue-shared regulatory regions and with those switching between promoter and enhancer signatures across species.

CONCLUSIONS

Our analyses of the tissue-specificity and evolutionary stability among promoters and enhancers reveal how specific LINE families have helped shape the dynamic mammalian regulome.

Histone Modifying Enzymes in Gynaecological Cancers

Simple Summary

Epigenetics is a process that allows genetic control, without the involvement of sequence changes to DNA or genes. In cancer, epigenetics is a key event in tumour development that can alter the expression of cancer driver genes and result in genomic instability. Due to the critical role of epigenetics in malignant transformation, therapies that target these processes have been developed to treat cancer. Here, we provide a summary of the epigenetic changes that have been described in a variety of gynaecological cancers. We then highlight how these changes are being targeted in preclinical models and clinical trials for gynaecological cancers.

Abstract

Genetic and epigenetic factors contribute to the development of cancer. Epigenetic dysregulation is common in gynaecological cancers and includes altered methylation at CpG islands in gene promoter regions, global demethylation that leads to genome instability and histone modifications. Histones are a major determinant of chromosomal conformation and stability, and unlike DNA methylation, which is generally associated with gene silencing, are amenable to post-translational modifications that induce facultative chromatin regions, or condensed transcriptionally silent regions that decondense resulting in global alteration of gene expression. In comparison, other components, crucial to the manipulation of chromatin dynamics, such as histone modifying enzymes, are not as well-studied. Inhibitors targeting DNA modifying enzymes, particularly histone modifying enzymes represent a potential cancer treatment. Due to the ability of epigenetic therapies to target multiple pathways simultaneously, tumours with complex mutational landscapes affected by multiple driver mutations may be most amenable to this type of inhibitor. Interrogation of the actionable landscape of different gynaecological cancer types has revealed that some patients have biomarkers which indicate potential sensitivity to epigenetic inhibitors. In this review we describe the role of epigenetics in gynaecological cancers and highlight how it may exploited for treatment.

A machine learning-based framework for modeling transcription elongation

RNA polymerase II (Pol II) generally pauses at certain positions along gene bodies, thereby interrupting the transcription elongation process, which is often coupled with various important biological functions, such as precursor mRNA splicing and gene expression regulation. Characterizing the transcriptional elongation dynamics can thus help us understand many essential biological processes in eukaryotic cells. However, experimentally measuring Pol II elongation rates is generally time and resource consuming. We developed PEPMAN (polymerase II elongation pausing modeling through attention-based deep neural network), a deep learning-based model that accurately predicts Pol II pausing sites based on the native elongating transcript sequencing (NET-seq) data. Through fully taking advantage of the attention mechanism, PEPMAN is able to decipher important sequence features underlying Pol II pausing. More importantly, we demonstrated that the analyses of the PEPMAN-predicted results around various types of alternative splicing sites can provide useful clues into understanding the cotranscriptional splicing events. In addition, associating the PEPMAN prediction results with different epigenetic features can help reveal important factors related to the transcription elongation process. All these results demonstrated that PEPMAN can provide a useful and effective tool for modeling transcription elongation and understanding the related biological factors from available high-throughput sequencing data.

The Role of H3K4 Trimethylation in CpG Islands Hypermethylation in Cancer

CpG methylation in transposons, exons, introns and intergenic regions is important for long-term silencing, silencing of parasitic sequences and alternative promoters, regulating imprinted gene expression and determining X chromosome inactivation. Promoter CpG islands, although rich in CpG dinucleotides, are unmethylated and remain so during all phases of mammalian embryogenesis and development, except in specific cases. The biological mechanisms that contribute to the maintenance of the unmethylated state of CpG islands remain elusive, but the modification of established DNA methylation patterns is a common feature in all types of tumors and is considered as an event that intrinsically, or in association with genetic lesions, feeds carcinogenesis. In this review, we focus on the latest results describing the role that the levels of H3K4 trimethylation may have in determining the aberrant hypermethylation of CpG islands in tumors.

A holistic view of mouse enhancer architectures reveals analogous pleiotropic effects and correlation with human disease

BACKGROUND

Efforts to elucidate the function of enhancers in vivo are underway but their vast numbers alongside differing enhancer architectures make it difficult to determine their impact on gene activity. By systematically annotating multiple mouse tissues with super- and typical-enhancers, we have explored their relationship with gene function and phenotype.

RESULTS

Though super-enhancers drive high total- and tissue-specific expression of their associated genes, we find that typical-enhancers also contribute heavily to the tissue-specific expression landscape on account of their large numbers in the genome. Unexpectedly, we demonstrate that both enhancer types are preferentially associated with relevant 'tissue-type' phenotypes and exhibit no difference in phenotype effect size or pleiotropy. Modelling regulatory data alongside molecular data, we built a predictive model to infer gene-phenotype associations and use this model to predict potentially novel disease-associated genes.

CONCLUSION

Overall our findings reveal that differing enhancer architectures have a similar impact on mammalian phenotypes whilst harbouring differing cellular and expression effects. Together, our results systematically characterise enhancers with predicted phenotypic traits endorsing the role for both types of enhancers in human disease and disorders.

Comprehensive multi-omics analysis uncovers a group of TGF-β-regulated genes among lncRNA EPR direct transcriptional targets

Long non-coding RNAs (lncRNAs) can affect multiple layers of gene expression to control crucial cellular functions. We have previously demonstrated that the lncRNA EPR, by controlling gene expression at different levels, affects cell proliferation and migration in cultured mammary gland cells and impairs breast tumor formation in an orthotopic transplant model in mice. Here, we used ChIRP-Seq to identify EPR binding sites on chromatin of NMuMG mammary gland cells overexpressing EPR and identified its trans binding sites in the genome. Then, with the purpose of relating EPR/chromatin interactions to the reshaping of the epitranscriptome landscape, we profiled histone activation marks at promoter/enhancer regions by ChIP-Seq. Finally, we integrated data derived from ChIRP-Seq, ChIP-Seq as well as RNA-Seq in a comprehensive analysis and we selected a group of bona fide direct transcriptional targets of EPR. Among them, we identified a subset of EPR targets whose expression is controlled by TGF-β with one of them—Arrdc3—being able to modulate Epithelial to Mesenchymal Transition. This experimental framework allowed us to correlate lncRNA/chromatin interactions with the real outcome of gene expression and to start defining the gene network regulated by EPR as a component of the TGF-β pathway.

Histone modifications, DNA methylation, and the epigenetic code of alcohol use disorder

Alcohol use disorder (AUD) is a leading cause of morbidity and mortality. Despite AUD's substantial contributions to lost economic productivity and quality of life, there are only a limited number of approved drugs for treatment of AUD in the United States. This chapter will update progress made on the epigenetic basis of AUD, with particular focus on histone post-translational modifications and DNA methylation and how these two epigenetic mechanisms interact to contribute to neuroadaptive processes leading to initiation, maintenance and progression of AUD pathophysiology. We will also evaluate epigenetic therapeutic strategies that have arisen from preclinical models of AUD and epigenetic biomarkers that have been discovered in human populations with AUD.

Altered Enhancer and Promoter Usage Leads to Differential Gene Expression in the Normal and Failed Human Heart

BACKGROUND

The failing heart is characterized by changes in gene expression. However, the regulatory regions of the genome that drive these gene expression changes have not been well defined in human hearts.

METHODS

To define genome-wide enhancer and promoter use in heart failure, cap analysis of gene expression sequencing was applied to 3 healthy and 4 failed human hearts to identify promoter and enhancer regions used in left ventricles. Healthy hearts were derived from donors unused for transplantation and failed hearts were obtained as discarded tissue after transplantation.

RESULTS

Cap analysis of gene expression sequencing identified a combined potential for ≈23 000 promoters and ≈5000 enhancers active in human left ventricles. Of these, 17 000 promoters and 1800 enhancers had additional support for their regulatory function. Comparing promoter usage between healthy and failed hearts highlighted promoter shifts which altered aminoterminal protein sequences. Enhancer usage between healthy and failed hearts identified a majority of differentially used heart failure enhancers were intronic and primarily localized within the first intron, revealing this position as a common feature associated with tissue-specific gene expression changes in the heart.

CONCLUSIONS

This data set defines the dynamic genomic regulatory landscape underlying heart failure and serves as an important resource for understanding genetic contributions to cardiac dysfunction. Additionally, regulatory changes contributing to heart failure are attractive therapeutic targets for controlling ventricular remodeling and clinical progression.

Spirits in the Material World: Enhancer RNAs in Transcriptional Regulation

Responses to developmental and environmental cues depend on precise spatiotemporal control of gene transcription. Enhancers, which comprise DNA elements bound by regulatory proteins, can activate target genes in response to these external signals. Recent studies have shown that enhancers are transcribed to produce enhancer RNAs (eRNAs). Do eRNAs play a functional role in activating gene expression or are they non-functional byproducts of nearby transcription machinery? The unstable nature of eRNAs and over-reliance on knockdown approaches have made elucidating the possible functions of eRNAs challenging. We focus here on studies using cloned eRNAs to study their function as transcripts, revealing roles for eRNAs in enhancer-promoter looping, recruiting transcriptional machinery, and facilitating RNA polymerase pause-release to regulate gene expression.

Enhancer dependence of cell-type-specific gene expression increases with developmental age

Gene regulatory logic reflects the occupancy of cis elements by transcription factors and the configuration of promoters and enhancers. As the majority of genome-wide analyses have focused on adult cells, scant attention has been paid to embryonic cells, other than embryonic stem cells. Focusing on genome-wide comparative analyses of two stages of erythroblasts, we discovered that regulation of embryonic-specific genes is promoter-centric through Gata1, whereas adult-specific control is combinatorial enhancer-driven and requires Myb, which is confirmed by increased enhancer–promoter interactions of adult specific genes. Extending genome-wide comparative analyses more broadly to available datasets of diverse mouse and human cells and tissues, we conclude that the progressively increased enhancer dependence of cell-type–specific genes with developmental age is conserved during development.

How overall principles of cell-type–specific gene regulation (the “logic”) may change during ontogeny is largely unexplored. We compared transcriptomic, epigenomic, and three-dimensional (3D) genomic profiles in embryonic (EryP) and adult (EryD) erythroblasts. Despite reduced chromatin accessibility compared to EryP, distal chromatin of EryD is enriched in H3K27ac, Gata1, and Myb occupancy. EryP-/EryD-shared enhancers are highly correlated with red blood cell identity genes, whereas cell-type–specific regulation employs different cis elements in EryP and EryD cells. In contrast to EryP-specific genes, which exhibit promoter-centric regulation through Gata1, EryD-specific genes rely more on distal enhancers for regulation involving Myb-mediated enhancer activation. Gata1 HiChIP demonstrated an overall increased enhancer–promoter interactions at EryD-specific genes, whereas genome editing in selected loci confirmed distal enhancers are required for gene expression in EryD but not in EryP. Applying a metric for enhancer dependence of transcription, we observed a progressive reliance on cell-specific enhancers with increasing ontogenetic age among diverse tissues of mouse and human origin. Our findings highlight fundamental and conserved differences at distinct developmental stages, characterized by simpler promoter-centric regulation of cell-type–specific genes in embryonic cells and increased combinatorial enhancer-driven control in adult cells.

Interactions With Histone H3 & Tools to Study Them

Histones are an integral part of chromatin and thereby influence its structure, dynamics, and functions. The effects of histone variants, posttranslational modifications, and binding proteins is therefore of great interest. From the moment that they are deposited on chromatin, nucleosomal histones undergo dynamic changes in function of the cell cycle, and as DNA is transcribed and replicated. In the process, histones are not only modified and bound by various proteins, but also shuffled, evicted, or replaced. Technologies and tools to study such dynamic events continue to evolve and better our understanding of chromatin and of histone proteins proper. Here, we provide an overview of H3.1 and H3.3 histone dynamics throughout the cell cycle, while highlighting some of the tools used to study their protein–protein interactions. We specifically discuss how histones are chaperoned, modified, and bound by various proteins at different stages of the cell cycle. Established and select emerging technologies that furthered (or have a high potential of furthering) our understanding of the dynamic histone–protein interactions are emphasized. This includes experimental tools to investigate spatiotemporal changes on chromatin, the role of histone chaperones, histone posttranslational modifications, and histone-binding effector proteins.

MTTFsite: cross-cell type TF binding site prediction by using multi-task learning

Motivation

The prediction of transcription factor binding sites (TFBSs) is crucial for gene expression analysis. Supervised learning approaches for TFBS predictions require large amounts of labeled data. However, many TFs of certain cell types either do not have sufficient labeled data or do not have any labeled data.

Results

In this paper, a multi-task learning framework (called MTTFsite) is proposed to address the lack of labeled data problem by leveraging on labeled data available in cross-cell types. The proposed MTTFsite contains a shared CNN to learn common features for all cell types and a private CNN for each cell type to learn private features. The common features are aimed to help predicting TFBSs for all cell types especially those cell types that lack labeled data. MTTFsite is evaluated on 241 cell type TF pairs and compared with a baseline method without using any multi-task learning model and a fully shared multi-task model that uses only a shared CNN and do not use private CNNs. For cell types with insufficient labeled data, results show that MTTFsite performs better than the baseline method and the fully shared model on more than 89% pairs. For cell types without any labeled data, MTTFsite outperforms the baseline method and the fully shared model by more than 80 and 93% pairs, respectively. A novel gene expression prediction method (called TFChrome) using both MTTFsite and histone modification features is also presented. Results show that TFBSs predicted by MTTFsite alone can achieve good performance. When MTTFsite is combined with histone modification features, a significant 5.7% performance improvement is obtained.

Availability and implementation

The resource and executable code are freely available at http://hlt.hitsz.edu.cn/MTTFsite/ and http://www.hitsz-hlt.com:8080/MTTFsite/.

Supplementary information

Supplementary data are available at Bioinformatics online.