Transcribed enhancers lead waves of coordinated transcription in transitioning mammalian cells

Cap Analysis of Gene Expression (CAGE): A Quantitative and Genome-Wide Assay of Transcription Start Sites

Cap analysis of gene expression (CAGE) is an approach to identify and monitor the activity (transcription initiation frequency) of transcription start sites (TSSs) at single base-pair resolution across the genome. It has been effectively used to identify active promoter and enhancer regions in cancer cells, with potential utility to identify key factors to immunotherapy. Here, we overview a series of CAGE protocols and describe detailed experimental steps of the latest protocol based on the Illumina sequencing platform; both experimental steps (see Subheadings 3.1-3.11) and computational processing steps (see Subheadings 3.12-3.20) are described.

Lessons from eRNAs: understanding transcriptional regulation through the lens of nascent RNAs

Nascent transcription assays, such as global run-on sequencing (GRO-seq) and precision run-on sequencing (PRO-seq), have uncovered a myriad of unstable RNAs being actively produced from numerous sites genome-wide. These transcripts provide a more complete and immediate picture of the impact of regulatory events. Transcription factors recruit RNA polymerase II, effectively initiating the process of transcription; repressors inhibit polymerase recruitment. Efficiency of recruitment is dictated by sequence elements in and around the RNA polymerase loading zone. A combination of sequence elements and RNA binding proteins subsequently influence the ultimate stability of the resulting transcript. Some of these transcripts are capable of providing feedback on the process, influencing subsequent transcription. By monitoring RNA polymerase activity, nascent assays provide insights into every step of the regulated process of transcription.

The role of enhancer RNAs in epigenetic regulation of gene expression

Since the discovery that enhancers can support transcription, the roles of enhancer RNAs have remained largely elusive. We identified that enhancer RNAs interact with and augment bromodomain engagement with acetylated chromatin. Here, we discuss our recent findings and the potential mechanisms underlying the regulation and functions of enhancer RNA-bromodomain associations.

Automated minute scale RNA-seq of pluripotent stem cell differentiation reveals early divergence of human and mouse gene expression kinetics

Pluripotent stem cells retain the developmental timing of their species of origin in vitro, an observation that suggests the existence of a cell-intrinsic developmental clock, yet the nature and machinery of the clock remain a mystery. We hypothesize that one possible component may lie in species-specific differences in the kinetics of transcriptional responses to differentiation signals. Using a liquid-handling robot, mouse and human pluripotent stem cells were exposed to identical neural differentiation conditions and sampled for RNA-sequencing at high frequency, every 4 or 10 minutes, for the first 10 hours of differentiation to test for differences in transcriptomic response rates. The majority of initial transcriptional responses occurred within a rapid window in the first minutes of differentiation for both human and mouse stem cells. Despite similarly early onsets of gene expression changes, we observed shortened and condensed gene expression patterns in mouse pluripotent stem cells compared to protracted trends in human pluripotent stem cells. Moreover, the speed at which individual genes were upregulated, as measured by the slopes of gene expression changes over time, was significantly faster in mouse compared to human cells. These results suggest that downstream transcriptomic response kinetics to signaling cues are faster in mouse versus human cells, and may offer a partial account for the vast differences in developmental rates across species.

Identification and Massively Parallel Characterization of Regulatory Elements Driving Neural Induction

Epigenomic regulation and lineage-specific gene expression act in concert to drive cellular differentiation, but the temporal interplay between these processes is largely unknown. Using neural induction from human pluripotent stem cells (hPSCs) as a paradigm, we interrogated these dynamics by performing RNA sequencing (RNA-seq), chromatin immunoprecipitation sequencing (ChIP-seq), and assay for transposase accessible chromatin using sequencing (ATAC-seq) at seven time points during early neural differentiation. We found that changes in DNA accessibility precede H3K27ac, which is followed by gene expression changes. Using massively parallel reporter assays (MPRAs) to test the activity of 2,464 candidate regulatory sequences at all seven time points, we show that many of these sequences have temporal activity patterns that correlate with their respective cell-endogenous gene expression and chromatin changes. A prioritization method incorporating all genomic and MPRA data further identified key transcription factors involved in driving neural fate. These results provide a comprehensive resource of genes and regulatory elements that orchestrate neural induction and illuminate temporal frameworks during differentiation.

Transcriptional Responses to IFN-γ Require Mediator Kinase-Dependent Pause Release and Mechanistically Distinct CDK8 and CDK19 Functions

Transcriptional responses to external stimuli remain poorly understood. Using global nuclear run-on followed by sequencing (GRO-seq) and precision nuclear run-on sequencing (PRO-seq), we show that CDK8 kinase activity promotes RNA polymerase II pause release in response to interferon-γ (IFN-γ), a universal cytokine involved in immunity and tumor surveillance. The Mediator kinase module contains CDK8 or CDK19, which are presumed to be functionally redundant. We implemented cortistatin A, chemical genetics, transcriptomics, and other methods to decouple their function while assessing enzymatic versus structural roles. Unexpectedly, CDK8 and CDK19 regulated different gene sets via distinct mechanisms. CDK8-dependent regulation required its kinase activity, whereas CDK19 governed IFN-γ responses through its scaffolding function (i.e., it was kinase independent). Accordingly, CDK8, not CDK19, phosphorylates the STAT1 transcription factor (TF) during IFN-γ stimulation, and CDK8 kinase inhibition blocked activation of JAK-STAT pathway TFs. Cytokines such as IFN-γ rapidly mobilize TFs to "reprogram" cellular transcription; our results implicate CDK8 and CDK19 as essential for this transcriptional reprogramming.

Differential activity of transcribed enhancers in the prefrontal cortex of 537 cases with schizophrenia and controls

Transcription at enhancers is a widespread phenomenon which produces so-called enhancer RNA (eRNA) and occurs in an activity-dependent manner. However, the role of eRNA and its utility in exploring disease-associated changes in enhancer function, and the downstream coding transcripts that they regulate, is not well established. We used transcriptomic and epigenomic data to interrogate the relationship of eRNA transcription to disease status and how genetic variants alter enhancer transcriptional activity in the human brain. We combined RNA-seq data from 537 postmortem brain samples from the CommonMind Consortium with cap analysis of gene expression and enhancer identification, using the assay for transposase-accessible chromatin followed by sequencing (ATACseq). We find 118 differentially transcribed eRNAs in schizophrenia and identify schizophrenia-associated gene/eRNA co-expression modules. Perturbations of a key module are associated with the polygenic risk scores. Furthermore, we identify genetic variants affecting expression of 927 enhancers, which we refer to as enhancer expression quantitative loci or eeQTLs. Enhancer expression patterns are consistent across studies, including differentially expressed eRNAs and eeQTLs. Combining eeQTLs with a genome-wide association study of schizophrenia identifies a genetic variant that alters enhancer function and expression of its target gene, GOLPH3L. Our novel approach to analyzing enhancer transcription is adaptable to other large-scale, non-poly-A depleted, RNA-seq studies.

Immunity-associated long non-coding RNA and expression in response to bacterial infection in large yellow croaker (Larimichthys crocea)

Long non-coding RNA refers to an RNA transcript of a non-coding protein with a sequence length greater than 200 bp. More and more reports indicated that lncRNA was involved in the regulation of gene expression as a signalling molecule, an inducing molecule, a leader molecule and a scaffold molecule. Previous studies have sequenced the draft genome and several transcriptome data sets for protein-coding genes of the large yellow croaker (Larimichthys crocea), but little is known about the expression and function of lncRNAs in this species. In order to obtain a catalogue of lncRNAs for this croaker, Vibrio parahaemolyticus infection challenge experiment was conducted and long non-coding RNA sequences were obtained. Using high-throughput sequencing of lncRNA, a total of 73,233 high-confidence transcripts were reconstructed in 32,726 loci, recovering most of the expressed reference transcripts, and 6473 novel expressed loci were identified. The tissue expression profile revealed that most lacunas were specifically enriched in distinct tissues. A set of 163 lncRNAs were identified as being specifically expressed in the spleen and may be involved in the immune response. It is the first time to identify specific lncRNAs in the L. crocea systematically in this croaker, aiming to benefit the future genomic study of this species.

Detailed Functional Characterization of a Waist-Hip Ratio Locus in 7p15.2 Defines an Enhancer Controlling Adipocyte Differentiation

We combined CAGE sequencing in human adipocytes during differentiation with data from genome-wide association studies to identify an enhancer in the SNX10 locus on chromosome 7, presumably involved in body fat distribution. Using reporter assays and CRISPR-Cas9 gene editing in human cell lines, we characterized the role of the enhancer in adipogenesis. The enhancer was active during adipogenesis and responded strongly to insulin and isoprenaline. The allele associated with increased waist-hip ratio in human genetic studies was associated with higher enhancer activity. Mutations of the enhancer resulted in less adipocyte differentiation. RNA sequencing of cells with disrupted enhancer showed reduced expression of established adipocyte markers, such as ADIPOQ and LPL, and identified CHI3L1 on chromosome 1 as a potential gene involved in adipocyte differentiation. In conclusion, we identified and characterized an enhancer in the SNX10 locus and outlined its plausible mechanisms of action and downstream targets.

Interrogation of nonconserved human adipose lincRNAs identifies a regulatory role of linc-ADAL in adipocyte metabolism

Long intergenic noncoding RNAs (lincRNAs) have emerged as important modulators of cellular functions. Most lincRNAs are not conserved among mammals, raising the fundamental question of whether nonconserved adipose-expressed lincRNAs are functional. To address this, we performed deep RNA sequencing of gluteal subcutaneous adipose tissue from 25 healthy humans. We identified 1001 putative lincRNAs expressed in all samples through de novo reconstruction of noncoding transcriptomes and integration with existing lincRNA annotations. One hundred twenty lincRNAs had adipose-enriched expression, and 54 of these exhibited peroxisome proliferator-activated receptor γ (PPARγ) or CCAAT/enhancer binding protein α (C/EBPα) binding at their loci. Most of these adipose-enriched lincRNAs (~85%) were not conserved in mice, yet on average, they showed degrees of expression and binding of PPARγ and C/EBPα similar to those displayed by conserved lincRNAs. Most adipose lincRNAs differentially expressed (n = 53) in patients after bariatric surgery were nonconserved. The most abundant adipose-enriched lincRNA in our subcutaneous adipose data set, linc-ADAL, was nonconserved, up-regulated in adipose depots of obese individuals, and markedly induced during in vitro human adipocyte differentiation. We demonstrated that linc-ADAL interacts with heterogeneous nuclear ribonucleoprotein U (hnRNPU) and insulin-like growth factor 2 mRNA binding protein 2 (IGF2BP2) at distinct subcellular locations to regulate adipocyte differentiation and lipogenesis.

MicroRNA-27a/b-3p and PPARG regulate SCAMP3 through a feed-forward loop during adipogenesis

MicroRNAs (miRNA) modulate gene expression through feed-back and forward loops. Previous studies identified miRNAs that regulate transcription factors, including Peroxisome Proliferator Activated Receptor Gamma (PPARG), in adipocytes, but whether they influence adipogenesis via such regulatory loops remain elusive. Here we predicted and validated a novel feed-forward loop regulating adipogenesis and involved miR-27a/b-3p, PPARG and Secretory Carrier Membrane Protein 3 (SCAMP3). In this loop, expression of both PPARG and SCAMP3 was independently suppressed by miR-27a/b-3p overexpression. Knockdown of PPARG downregulated SCAMP3 expression at the late phase of adipogenesis, whereas reduction of SCAMP3 mRNA levels increased PPARG expression at early phase in differentiation. The latter was accompanied with upregulation of adipocyte-enriched genes, including ADIPOQ and FABP4, suggesting an anti-adipogenic role for SCAMP3. PPARG and SCAMP3 exhibited opposite behaviors regarding correlations with clinical phenotypes, including body mass index, body fat mass, adipocyte size, lipolytic and lipogenic capacity, and secretion of pro-inflammatory cytokines. While adipose PPARG expression was associated with more favorable metabolic phenotypes, SCAMP3 expression was linked to increased fat mass and insulin resistance. Together, we identified a feed-forward loop through which miR-27a/b-3p, PPARG and SCAMP3 cooperatively fine tune the regulation of adipogenesis, which potentially may impact whole body metabolism.

Identification of novel cerebellar developmental transcriptional regulators with motif activity analysis

Background

The work of the FANTOM5 Consortium has brought forth a new level of understanding of the regulation of gene transcription and the cellular processes involved in creating diversity of cell types. In this study, we extended the analysis of the FANTOM5 Cap Analysis of Gene Expression (CAGE) transcriptome data to focus on understanding the genetic regulators involved in mouse cerebellar development.

Results

We used the HeliScopeCAGE library sequencing on cerebellar samples over 8 embryonic and 4 early postnatal times. This study showcases temporal expression pattern changes during cerebellar development. Through a bioinformatics analysis that focused on transcription factors, their promoters and binding sites, we identified genes that appear as strong candidates for involvement in cerebellar development. We selected several candidate transcriptional regulators for validation experiments including qRT-PCR and shRNA transcript knockdown. We observed marked and reproducible developmental defects in Atf4, Rfx3, and Scrt2 knockdown embryos, which support the role of these genes in cerebellar development.

Conclusions

The successful identification of these novel gene regulators in cerebellar development demonstrates that the FANTOM5 cerebellum time series is a high-quality transcriptome database for functional investigation of gene regulatory networks in cerebellar development.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6063-9) contains supplementary material, which is available to authorized users.

NET-CAGE characterizes the dynamics and topology of human transcribed cis-regulatory elements

Promoters and enhancers are key cis-regulatory elements, but how they operate to generate cell type-specific transcriptomes is not fully understood. We developed a simple and robust method, native elongating transcript-cap analysis of gene expression (NET-CAGE), to sensitively detect 5' ends of nascent RNAs in diverse cells and tissues, including unstable transcripts such as enhancer-derived RNAs. We studied RNA synthesis and degradation at the transcription start site level, characterizing the impact of differential promoter usage on transcript stability. We quantified transcription from cis-regulatory elements without the influence of RNA turnover, and show that enhancer-promoter pairs are generally activated simultaneously on stimulation. By integrating NET-CAGE data with chromatin interaction maps, we show that cis-regulatory elements are topologically connected according to their cell type specificity. We identified new enhancers with high sensitivity, and delineated primary locations of transcription within super-enhancers. Our NET-CAGE dataset derived from human and mouse cells expands the FANTOM5 atlas of transcribed enhancers, with broad applicability to biomedical research.

Are TADs supercoiled?

Topologically associating domains (TADs) are megabase-sized building blocks of interphase chromosomes in higher eukaryotes. TADs are chromosomal regions with increased frequency of internal interactions. On average a pair of loci separated by a given genomic distance contact each other 2–3 times more frequently when they are in the same TAD as compared to a pair of loci located in two neighbouring TADs. TADs are also functional blocks of chromosomes as enhancers and their cognate promoters are normally located in the same TAD, even if their genomic distance from each other can be as large as a megabase. The internal structure of TADs, causing their increased frequency of internal interactions, is not established yet. We survey here experimental studies investigating presence of supercoiling in interphase chromosomes. We also review numerical simulation studies testing whether transcription-induced supercoiling of chromatin fibres can explain how TADs are formed and how they can assure very efficient interactions between enhancers and their cognate promoters located in the same TAD.

JASPAR 2018: update of the open-access database of transcription factor binding profiles and its web framework

JASPAR (http://jaspar.genereg.net) is an open-access database of curated, non-redundant transcription factor (TF)-binding profiles stored as position frequency matrices (PFMs) and TF flexible models (TFFMs) for TFs across multiple species in six taxonomic groups. In the 2018 release of JASPAR, the CORE collection has been expanded with 322 new PFMs (60 for vertebrates and 262 for plants) and 33 PFMs were updated (24 for vertebrates, 8 for plants and 1 for insects). These new profiles represent a 30% expansion compared to the 2016 release. In addition, we have introduced 316 TFFMs (95 for vertebrates, 218 for plants and 3 for insects). This release incorporates clusters of similar PFMs in each taxon and each TF class per taxon. The JASPAR 2018 CORE vertebrate collection of PFMs was used to predict TF-binding sites in the human genome. The predictions are made available to the scientific community through a UCSC Genome Browser track data hub. Finally, this update comes with a new web framework with an interactive and responsive user-interface, along with new features. All the underlying data can be retrieved programmatically using a RESTful API and through the JASPAR 2018 R/Bioconductor package.

Transcriptional control by enhancers and enhancer RNAs

The regulation of gene expression is a fundamental cellular process and its misregulation is a key component of disease. Enhancers are one of the most salient regulatory elements in the genome and help orchestrate proper spatiotemporal gene expression during development, in homeostasis, and in response to signaling. Notably, molecular aberrations at enhancers, such as translocations and single nucleotide polymorphisms, are emerging as an important source of human variation and susceptibility to disease. Herein we discuss emerging paradigms addressing how genes are regulated by enhancers, common features of active enhancers, and how non-coding enhancer RNAs (eRNAs) can direct gene expression programs that underlie cellular phenotypes. We survey the current evidence, which suggests that eRNAs can bind to transcription factors, mediate enhancer-promoter interactions, influence RNA Pol II elongation, and act as decoys for repressive cofactors. Furthermore, we discuss current methodologies for the identification of eRNAs and novel approaches to elucidate their functions.

Phenotypic impacts of CSF1R deficiencies in humans and model organisms

Mϕ proliferation, differentiation, and survival are controlled by signals from the Mϕ CSF receptor (CSF1R). Mono-allelic gain-of-function mutations in CSF1R in humans are associated with an autosomal-dominant leukodystrophy and bi-allelic loss-of-function mutations with recessive skeletal dysplasia, brain disorders, and developmental anomalies. Most of the phenotypes observed in these human disease states are also observed in mice and rats with loss-of-function mutations in Csf1r or in Csf1 encoding one of its two ligands. Studies in rodent models also highlight the importance of genetic background and likely epistatic interactions between Csf1r and other loci. The impacts of Csf1r mutations on the brain are usually attributed solely to direct impacts on microglial number and function. However, analysis of hypomorphic Csf1r mutants in mice and several other lines of evidence suggest that primary hydrocephalus and loss of the physiological functions of Mϕs in the periphery contribute to the development of brain pathology. In this review, we outline the evidence that CSF1R is expressed exclusively in mononuclear phagocytes and explore the mechanisms linking CSF1R mutations to pleiotropic impacts on postnatal growth and development.

The many lives of KATs - detectors, integrators and modulators of the cellular environment

Research over the past three decades has firmly established lysine acetyltransferases (KATs) as central players in regulating transcription. Recent advances in genomic sequencing, metabolomics, animal models and mass spectrometry technologies have uncovered unexpected new roles for KATs at the nexus between the environment and transcriptional regulation. Thousands of reversible acetylation sites have been mapped in the proteome that respond dynamically to the cellular milieu and maintain major processes such as metabolism, autophagy and stress response. Concurrently, researchers are continuously uncovering how deregulation of KAT activity drives disease, including cancer and developmental syndromes characterized by severe intellectual disability. These novel findings are reshaping our view of KATs away from mere modulators of chromatin to detectors of the cellular environment and integrators of diverse signalling pathways with the ability to modify cellular phenotype.

The RNA Polymerase II Factor RPAP1 Is Critical for Mediator-Driven Transcription and Cell Identity

The RNA polymerase II-associated protein 1 (RPAP1) is conserved across metazoa and required for stem cell differentiation in plants; however, very little is known about its mechanism of action or its role in mammalian cells. Here, we report that RPAP1 is essential for the expression of cell identity genes and for cell viability. Depletion of RPAP1 triggers cell de-differentiation, facilitates reprogramming toward pluripotency, and impairs differentiation. Mechanistically, we show that RPAP1 is essential for the interaction between RNA polymerase II (RNA Pol II) and Mediator, as well as for the recruitment of important regulators, such as the Mediator-specific RNA Pol II factor Gdown1 and the C-terminal domain (CTD) phosphatase RPAP2. In agreement, depletion of RPAP1 diminishes the loading of total and Ser5-phosphorylated RNA Pol II on many genes, with super-enhancer-driven genes among the most significantly downregulated. We conclude that Mediator/RPAP1/RNA Pol II is an ancient module, conserved from plants to mammals, critical for establishing and maintaining cell identity.

•

RPAP1 is an RNA Pol II regulator, conserved from plants to mammals

•

RPAP1 depletion erases cell identity gene expression, triggering de-differentiation

•

Mechanistically, RPAP1 is critical for the Mediator-RNA Pol II interaction

•

RPAP1 preferentially contributes to enhancer-driven gene transcription

Genome-wide enhancer annotations differ significantly in genomic distribution, evolution, and function

Background

Non-coding gene regulatory enhancers are essential to transcription in mammalian cells. As a result, a large variety of experimental and computational strategies have been developed to identify cis-regulatory enhancer sequences. Given the differences in the biological signals assayed, some variation in the enhancers identified by different methods is expected; however, the concordance of enhancers identified by different methods has not been comprehensively evaluated. This is critically needed, since in practice, most studies consider enhancers identified by only a single method. Here, we compare enhancer sets from eleven representative strategies in four biological contexts.

Results

All sets we evaluated overlap significantly more than expected by chance; however, there is significant dissimilarity in their genomic, evolutionary, and functional characteristics, both at the element and base-pair level, within each context. The disagreement is sufficient to influence interpretation of candidate SNPs from GWAS studies, and to lead to disparate conclusions about enhancer and disease mechanisms. Most regions identified as enhancers are supported by only one method, and we find limited evidence that regions identified by multiple methods are better candidates than those identified by a single method. As a result, we cannot recommend the use of any single enhancer identification strategy in all settings.

Conclusions

Our results highlight the inherent complexity of enhancer biology and identify an important challenge to mapping the genetic architecture of complex disease. Greater appreciation of how the diverse enhancer identification strategies in use today relate to the dynamic activity of gene regulatory regions is needed to enable robust and reproducible results.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5779-x) contains supplementary material, which is available to authorized users.

Characterising the genetic basis of immune response variation to identify causal mechanisms underlying disease susceptibility

Over the last 10 years, genome-wide association studies (GWAS) have identified hundreds of susceptibility loci for autoimmune diseases. However, despite increasing power for the detection of both common and rare coding variants affecting disease susceptibility, a large fraction of disease heritability has remained unexplained. In addition, a majority of the identified loci are located in noncoding regions, and translation of disease-associated loci into new biological insights on the etiology of immune disorders has been lagging. This highlights the need for a better understanding of noncoding variation and new strategies to identify causal genes at disease loci. In this review, I will first detail the molecular basis of gene expression and review the various mechanisms that contribute to alter gene activity at the transcriptional and post-transcriptional level. I will then review the findings from 10 years of functional genomics studies regarding the genetics on gene expression, in particular in the context of infection. Finally, I will discuss the extent to which genetic variants that modulate gene expression at transcriptional and post-transcriptional level contribute to disease susceptibility and present strategies to leverage this information for the identification of causal mechanisms at disease loci in the era of whole genome sequencing.

Novel Bioinformatics Approach Identifies Transcriptional Profiles of Lineage-Specific Transposable Elements at Distinct Loci in the Human Dorsolateral Prefrontal Cortex

Expression of transposable elements (TE) is transiently activated during human preimplantation embryogenesis in a developmental stage- and cell type-specific manner and TE-mediated epigenetic regulation is intrinsically wired in developmental genetic networks in human embryos and embryonic stem cells. However, there are no systematic studies devoted to a comprehensive analysis of the TE transcriptome in human adult organs and tissues, including human neural tissues. To investigate TE expression in the human Dorsolateral Prefrontal Cortex (DLPFC), we developed and validated a straightforward analytical approach to chart quantitative genome-wide expression profiles of all annotated TE loci based on unambiguous mapping of discrete TE-encoded transcripts using a de novo assembly strategy. To initially evaluate the potential regulatory impact of DLPFC-expressed TE, we adopted a comparative evolutionary genomics approach across humans, primates, and rodents to document conservation patterns, lineage-specificity, and colocalizations with transcription factor binding sites mapped within primate- and human-specific TE. We identified 654,665 transcripts expressed from 477,507 distinct loci of different TE classes and families, the majority of which appear to have originated from primate-specific sequences. We discovered 4,687 human-specific and transcriptionally active TEs in DLPFC, of which the prominent majority (80.2%) appears spliced. Our analyses revealed significant associations of DLPFC-expressed TE with primate- and human-specific transcription factor binding sites, suggesting potential cross-talks of concordant regulatory functions. We identified 1,689 TEs differentially expressed in the DLPFC of Schizophrenia patients, a majority of which is located within introns of 1,137 protein-coding genes. Our findings imply that identified DLPFC-expressed TEs may affect human brain structures and functions following different evolutionary trajectories. On one side, hundreds of thousands of TEs maintained a remarkably high conservation for ∼8 My of primates’ evolution, suggesting that they are likely conveying evolutionary-constrained primate-specific regulatory functions. In parallel, thousands of transcriptionally active human-specific TE loci emerged more recently, suggesting that they could be relevant for human-specific behavioral or cognitive functions.

Computational Biology Solutions to Identify Enhancers-target Gene Pairs

Enhancers are non-coding regulatory elements that are distant from their target gene. Their characterization still remains elusive especially due to challenges in achieving a comprehensive pairing of enhancers and target genes. A number of computational biology solutions have been proposed to address this problem leveraging the increasing availability of functional genomics data and the improved mechanistic understanding of enhancer action. In this review we focus on computational methods for genome-wide definition of enhancer-target gene pairs. We outline the different classes of methods, as well as their main advantages and limitations. The types of information integrated by each method, along with details on their applicability are presented and discussed. We especially highlight the technical challenges that are still unresolved and hamper the effective achievement of a satisfactory and comprehensive solution. We expect this field will keep evolving in the coming years due to the ever-growing availability of data and increasing insights into enhancers crucial role in regulating genome functionality.

Gene regulation and the architecture of complex human traits in the genomics era

Virtually all human psychological and behavioral traits are at least partially heritable. For nearly a century, classical genetic studies have sought to understand how genetic variation contributes to human variation in these traits. More recently, genome-wide association studies have identified large numbers of specific genetic variants linked with complex traits. Many of these variants fall outside of protein-coding genes, in putative gene regulatory elements. This suggests that some fraction of causal human genetic variation acts through gene regulation. New developments in the field of regulatory genomics offer resources and methods to understand how genetic variants that alter gene expression contribute to human psychology and risk for psychiatric disease.

Functional characterization of the C7ORF76 genomic region, a prominent GWAS signal for osteoporosis in 7q21.3

Genome-wide association studies (GWAS) have repeatedly identified genetic variants associated with bone mineral density (BMD) and osteoporotic fracture in non-coding regions of C7ORF76, a poorly studied gene of unknown function. The aim of the present study was to elucidate the causality and molecular mechanisms underlying the association. We re-sequenced the genomic region in two extreme BMD groups from the BARCOS cohort of postmenopausal women to search for functionally relevant variants. Eight selected variants were tested for association in the complete cohort and 2 of them (rs4342521 and rs10085588) were found significantly associated with lumbar spine BMD and nominally associated with osteoporotic fracture. cis-eQTL analyses of these 2 SNPs, together with SNP rs4727338 (GWAS lead SNP in Estrada et al., Nat Genet. 44:491-501, 2012), performed in human primary osteoblasts, disclosed a statistically significant influence on the expression of the proximal neighbouring gene SLC25A13 and a tendency on the distal SHFM1. We then studied the functionality of a putative upstream regulatory element (UPE), containing rs10085588. Luciferase reporter assays showed transactivation capability with a strong allele-dependent effect. Finally, 4C-seq experiments in osteoblastic cell lines showed that the UPE interacted with different tissue-specific enhancers and a lncRNA (LOC100506136) in the region. In summary, this study is the first one to analyse in depth the functionality of C7ORF76 genomic region. We provide functional regulatory evidence for the rs10085588, which may be a causal SNP within the 7q21.3 GWAS signal for osteoporosis.

Detection of condition-specific marker genes from RNA-seq data with MGFR

The identification of condition-specific genes is key to advancing our understanding of cell fate decisions and disease development. Differential gene expression analysis (DGEA) has been the standard tool for this task. However, the amount of samples that modern transcriptomic technologies allow us to study, makes DGEA a daunting task. On the other hand, experiments with low numbers of replicates lack the statistical power to detect differentially expressed genes. We have previously developed MGFM, a tool for marker gene detection from microarrays, that is particularly useful in the latter case. Here, we have adapted the algorithm behind MGFM to detect markers in RNA-seq data. MGFR groups samples with similar gene expression levels and flags potential markers of a sample type if their highest expression values represent all replicates of this type. We have benchmarked MGFR against other methods and found that its proposed markers accurately characterize the functional identity of different tissues and cell types in standard and single cell RNA-seq datasets. Then, we performed a more detailed analysis for three of these datasets, which profile the transcriptomes of different human tissues, immune and human blastocyst cell types, respectively. MGFR’s predicted markers were compared to gold-standard lists for these datasets and outperformed the other marker detectors. Finally, we suggest novel candidate marker genes for the examined tissues and cell types. MGFR is implemented as a freely available Bioconductor package (https://doi.org/doi:10.18129/B9.bioc.MGFR), which facilitates its use and integration with bioinformatics pipelines.

Enhancer RNAs: Insights Into Their Biological Role

Enhancers play a central role in the transcriptional regulation of metazoans. Almost a decade ago, the discovery of their pervasive transcription into noncoding RNAs, termed enhancer RNAs (eRNAs), opened a whole new field of study. The presence of eRNAs correlates with enhancer activity; however, whether they act as functional molecules remains controversial. Here we review direct experimental evidence supporting a functional role of eRNAs in transcription and provide a general pipeline that could help in the design of experimental approaches to investigate the function of eRNAs. We propose that induction of transcriptional activity at enhancers promotes an increase in its activity by an RNA-mediated titration of regulatory proteins that can impact different processes like chromatin accessibility or chromatin looping. In a few cases, transcripts originating from enhancers have acquired specific molecular functions to regulate gene expression. We speculate that these transcripts are either nonannotated long noncoding RNAs (lncRNAs) or are evolving toward functional lncRNAs. Further work will be needed to comprehend better the biological activity of these transcripts.

Sex-Biased lncRNAs Inversely Correlate With Sex-Opposite Gene Coexpression Networks in Diversity Outbred Mouse Liver

Sex differences in liver gene expression are determined by pituitary growth hormone secretion patterns, which regulate sex-dependent liver transcription factors and establish sex-specific chromatin states. Hypophysectomy (hypox) identifies two major classes of liver sex-biased genes, defined by their sex-dependent positive or negative responses to pituitary hormone ablation. However, the mechanisms that underlie each hypox-response class are unknown. We sought to discover candidate, regulatory, long noncoding RNAs (lncRNAs) controlling responsiveness to hypox. We characterized gene structures and expression patterns for 15,558 mouse liver-expressed lncRNAs, including many sex-specific lncRNAs regulated during postnatal development or subject to circadian regulation. Using the high natural allelic variance of Diversity Outbred (DO) mice, we discovered tightly coexpressed clusters of sex-specific protein-coding genes (gene modules) in male and female DO liver. Remarkably, many gene modules were strongly enriched for sex-specific genes within a single hypox-response class, indicating that the genetic heterogeneity of DO mice encompasses responsiveness to hypox. Moreover, several distant gene modules were enriched for gene subsets of the same hypox-response class, highlighting the complex regulation of hypox-responsiveness. Finally, we identified eight sex-specific lncRNAs with strong negative regulatory potential, as indicated by their strong negative correlation of expression across DO mouse livers with that of protein-coding gene modules enriched for genes of the opposite sex bias and inverse hypox-response class. These findings reveal an important role for genetic factors in regulating responsiveness to hypox, and present testable hypotheses for the roles of sex-biased liver lncRNAs in controlling the sex-bias of liver gene expression.

Nestin in immature embryonic neurons affects axon growth cone morphology and Semaphorin3a sensitivity

Correct wiring in the neocortex requires that responses to an individual guidance cue vary among neurons in the same location, and within the same neuron over time. Nestin is an atypical intermediate filament expressed strongly in neural progenitors and is thus used widely as a progenitor marker. Here we show a subpopulation of embryonic cortical neurons that transiently express nestin in their axons. Nestin expression is thus not restricted to neural progenitors, but persists for 2-3 d at lower levels in newborn neurons. We found that nestin-expressing neurons have smaller growth cones, suggesting that nestin affects cytoskeletal dynamics. Nestin, unlike other intermediate filament subtypes, regulates cdk5 kinase by binding the cdk5 activator p35. Cdk5 activity is induced by the repulsive guidance cue Semaphorin3a (Sema3a), leading to axonal growth cone collapse in vitro. Therefore, we tested whether nestin-expressing neurons showed altered responses to Sema3a. We find that nestin-expressing newborn neurons are more sensitive to Sema3a in a roscovitine-sensitive manner, whereas nestin knockdown results in lowered sensitivity to Sema3a. We propose that nestin functions in immature neurons to modulate cdk5 downstream of the Sema3a response. Thus, the transient expression of nestin could allow temporal and/or spatial modulation of a neuron's response to Sema3a, particularly during early axon guidance.

Competitive endogenous RNA is an intrinsic component of EMT regulatory circuits and modulates EMT

The competitive endogenous RNA (ceRNA) hypothesis suggests an intrinsic mechanism to regulate biological processes. However, whether the dynamic changes of ceRNAs can modulate miRNA activities remains controversial. Here, we examine the dynamics of ceRNAs during TGF-β-induced epithelial-to-mesenchymal transition (EMT). We observe that TGFBI, a transcript highly induced during EMT in A549 cells, acts as the ceRNA for miR-21 to modulate EMT. We further identify FN1 as the ceRNA for miR-200c in the canonical SNAIL-ZEB-miR200 circuit in MCF10A cells. Experimental assays and computational simulations demonstrate that the dynamically induced ceRNAs are directly coupled with the canonical double negative feedback loops and are critical to the induction of EMT. These results help to establish the relevance of ceRNA in cancer EMT and suggest that ceRNA is an intrinsic component of the EMT regulatory circuit and may represent a potential target to disrupt EMT during tumorigenesis.

Competitive endogenous RNAs help to regulate biological processes by regulating miRNA activity levels. Here the author show TGFBI acts as a ceRNA for miR-21 in the epithelial-to-mesenchymal transition.

Treg cells in autoimmunity: from identification to Treg-based therapies

Regulatory (Treg) cells are key regulators of inflammation and important for immune tolerance and homeostasis. A major progress has been made in the identification and classification of Treg cells. Due to technological advances, we have gained deep insights in the epigenetic regulation of Treg cells. The use of fate reporter mice allowed addressing the functional consequences of loss of Foxp3 expression. Depending on the environment Treg cells gain effector functions upon loss of Foxp3 expression. However, the traditional view that Treg cells become necessarily pathogenic by gaining effector functions was challenged by recent findings and supports the notion of Treg cell lineage plasticity. Treg cell stability is also a major issue for Treg cell therapies. Clinical trials are designed to use polyclonal Treg cells as therapeutic tools. Here, we summarize the role of Treg cells in selected autoimmune diseases and recent advances in the field of Treg targeted therapies.

Differential Targeting of c-Maf, Bach-1, and Elmo-1 by microRNA-143 and microRNA-365 Promotes the Intracellular Growth of Mycobacterium tuberculosis in Alternatively IL-4/IL-13 Activated Macrophages

Mycobacterium tuberculosis (Mtb) can subvert the host defense by skewing macrophage activation toward a less microbicidal alternative activated state to avoid classical effector killing functions. Investigating the molecular basis of this evasion mechanism could uncover potential candidates for host directed therapy against tuberculosis (TB). A limited number of miRNAs have recently been shown to regulate host-mycobacterial interactions. Here, we performed time course kinetics experiments on bone marrow-derived macrophages (BMDMs) and human monocyte-derived macrophages (MDMs) alternatively activated with IL-4, IL-13, or a combination of IL-4/IL-13, followed by infection with Mtb clinical Beijing strain HN878. MiR-143 and miR-365 were highly induced in Mtb-infected M(IL-4/IL-13) BMDMs and MDMs. Knockdown of miR-143 and miR-365 using antagomiRs decreased the intracellular growth of Mtb HN878, reduced the production of IL-6 and CCL5 and promoted the apoptotic death of Mtb HN878-infected M(IL-4/IL-13) BMDMs. Computational target prediction identified c-Maf, Bach-1 and Elmo-1 as potential targets for both miR-143 and miR-365. Functional validation using luciferase assay, RNA-pulldown assay and Western blotting revealed that c-Maf and Bach-1 are directly targeted by miR-143 while c-Maf, Bach-1, and Elmo-1 are direct targets of miR-365. Knockdown of c-Maf using GapmeRs promoted intracellular Mtb growth when compared to control treated M(IL-4/IL-13) macrophages. Meanwhile, the blocking of Bach-1 had no effect and blocking Elmo-1 resulted in decreased Mtb growth. Combination treatment of M(IL-4/IL-13) macrophages with miR-143 mimics or miR-365 mimics and c-Maf, Bach-1, or Elmo-1 gene-specific GapmeRs restored Mtb growth in miR-143 mimic-treated groups and enhanced Mtb growth in miR-365 mimics-treated groups, thus suggesting the Mtb growth-promoting activities of miR-143 and miR-365 are mediated at least partially through interaction with c-Maf, Bach-1, and Elmo-1. We further show that knockdown of miR-143 and miR-365 in M(IL-4/IL-13) BMDMs decreased the expression of HO-1 and IL-10 which are known targets of Bach-1 and c-Maf, respectively, with Mtb growth-promoting activities in macrophages. Altogether, our work reports a host detrimental role of miR-143 and miR-365 during Mtb infection and highlights for the first time the role and miRNA-mediated regulation of c-Maf, Bach-1, and Elmo-1 in Mtb-infected M(IL-4/IL-13) macrophages.

Dismissal of RNA Polymerase II Underlies a Large Ligand-Induced Enhancer Decommissioning Program

Nuclear receptors induce both transcriptional activation and repression programs responsible for development, homeostasis, and disease. Here, we report a previously overlooked enhancer decommissioning strategy underlying a large estrogen receptor alpha (ERα)-dependent transcriptional repression program. The unexpected signature for this E₂-induced program resides in indirect recruitment of ERα to a large cohort of pioneer factor basally active FOXA1-bound enhancers that lack cognate ERα DNA-binding elements. Surprisingly, these basally active estrogen-repressed (BAER) enhancers are decommissioned by ERα-dependent recruitment of the histone demethylase KDM2A, functioning independently of its demethylase activity. Rather, KDM2A tethers the E3 ubiquitin-protein ligase NEDD4 to ubiquitylate/dismiss Pol II to abrogate eRNA transcription, with consequent target gene downregulation. Thus, our data reveal that Pol II ubiquitylation/dismissal may serve as a potentially broad strategy utilized by indirectly bound nuclear receptors to abrogate large programs of pioneer factor-mediated, eRNA-producing enhancers.

Arginine to Glutamine Variant in Olfactomedin Like 3 (OLFML3) Is a Candidate for Severe Goniodysgenesis and Glaucoma in the Border Collie Dog Breed

Goniodysgenesis is a developmental abnormality of the anterior chamber of the eye. It is generally considered to be congenital in dogs (Canis lupus familiaris), and has been associated with glaucoma and blindness. Goniodysgenesis and early-onset glaucoma initially emerged in Border Collies in Australia in the late 1990s and have subsequently been found in this breed in Europe and the USA. The objective of the present study was to determine the genetic basis of goniodysgenesis in Border Collies. Clinical diagnosis was based on results of examinations by veterinary ophthalmologists of affected and unaffected dogs from eleven different countries. Genotyping using the Illumina high density canine single nucleotide variant genotyping chip was used to identify a candidate genetic region. There was a highly significant peak of association over chromosome 17, with a p-value of 2 × 10-13 Expression profiles and evolutionary conservation of candidate genes were assessed using public databases. Whole genome sequences of three dogs with glaucoma, three severely affected by goniodysgenesis and three unaffected dogs identified a missense variant in the olfactomedin like 3 (OLFML3) gene in all six affected animals. This was homozygous for the risk allele in all nine cases with glaucoma and 12 of 14 other severely affected animals. Of 67 reportedly unaffected animals, only one was homozygous for this variant (offspring of parents both with goniodysgenesis who were also homozygous for the variant). Analysis of pedigree information was consistent with an autosomal recessive mode of inheritance for severe goniodysgenesis (potentially leading to glaucoma) in this breed. The identification of a candidate genetic region and putative causative variant will aid breeders to reduce the frequency of goniodysgenesis and the risk of glaucoma in the Border Collie population.

Prediction of RNA-protein interactions using conjoint triad feature and chaos game representation

RNA-protein interactions (RPIs) play a very important role in a wide range of post-transcriptional regulations, and identifying whether a given RNA-protein pair can form interactions or not is a vital prerequisite for dissecting the regulatory mechanisms of functional RNAs. Currently, expensive and time-consuming biological assays can only determine a very small portion of all RPIs, which calls for computational approaches to help biologists efficiently and correctly find candidate RPIs. Here, we integrated a successful computing algorithm, conjoint triad feature (CTF), and another method, chaos game representation (CGR), for representing RNA-protein pairs and by doing so developed a prediction model based on these representations and random forest (RF) classifiers. When testing two benchmark datasets, RPI369 and RPI2241, the combined method (CTF+CGR) showed some superiority compared with four existing tools. Especially on RPI2241, the CTF+CGR method improved prediction accuracy (ACC) from 0.91 (the best record of all published works) to 0.95. When independently testing a newly constructed dataset, RPI1449, which only contained experimentally validated RPIs released between 2014 and 2016, our method still showed some generalization capability with an ACC of 0.75. Accordingly, we believe that our hybrid CTF+CGR method will be an important tool for predicting RPIs in the future.

Comprehensive profiling of the fission yeast transcription start site activity during stress and media response

Fission yeast, Schizosaccharomyces pombe, is an attractive model organism for transcriptional and chromatin biology research. Such research is contingent on accurate annotation of transcription start sites (TSSs). However, comprehensive genome-wide maps of TSSs and their usage across commonly applied laboratory conditions and treatments for S. pombe are lacking. To this end, we profiled TSS activity genome-wide in S. pombe cultures exposed to heat shock, nitrogen starvation, hydrogen peroxide and two commonly applied media, YES and EMM2, using Cap Analysis of Gene Expression (CAGE). CAGE-based annotation of TSSs is substantially more accurate than existing PomBase annotation; on average, CAGE TSSs fall 50-75 bp downstream of PomBase TSSs and co-localize with nucleosome boundaries. In contrast to higher eukaryotes, dispersed TSS distributions are not common in S. pombe. Our data recapitulate known S. pombe stress expression response patterns and identify stress- and media-responsive alternative TSSs. Notably, alteration of growth medium induces changes of similar magnitude as some stressors. We show a link between nucleosome occupancy and genetic variation, and that the proximal promoter region is genetically diverse between S. pombe strains. Our detailed TSS map constitutes a central resource for S. pombe gene regulation research.

Epigenetic control of early dendritic cell lineage specification by the transcription factor IRF8 in mice

Dendritic cells (DCs), which are vital for immune responses, are derived from bone marrow hematopoietic stem cells via common DC progenitors (CDPs). DC lineage fate decisions occurring at stages much earlier than CDPs have recently been recognized, yet the mechanism remains elusive. By single-cell RNA-sequencing, in vivo cell transfer experiments, and an assay for transposase-accessible chromatin sequencing using wild-type, IRF8-GFP chimera knock-in or IRF8-knockout mice, we demonstrate that IRF8 regulates chromatin at the lymphoid-primed multipotent progenitor (LMPP) stage to induce early commitment toward DCs. A low but significant expression of IRF8, a transcription factor essential for DC and monocyte development, was initiated in a subpopulation within LMPPs. These IRF8⁺ LMPPs were derived from IRF8^- LMPPs and predominantly produced DCs, especially classical DC1s, potentially via known progenitors, such as monocyte-DC progenitors, CDPs, and preclassical DCs. IRF8⁺ LMPPs did not generate significant numbers of monocytes, neutrophils, or lymphocytes. Although IRF8^- and IRF8⁺ LMPPs displayed very similar global gene expression patterns, the chromatin of enhancers near DC lineage genes was more accessible in IRF8⁺ LMPPs than in IRF8^- LMPPs, an epigenetic change dependent on IRF8. The majority of the genes epigenetically primed by IRF8 were still transcriptionally inactive at the LMPP stage, but were highly expressed in the downstream DC lineage populations such as CDPs. Therefore, early expression of the key transcription factor IRF8 changes chromatin states in otherwise multipotent progenitors, biasing their fate decision toward DCs.

Twenty-Seven Tamoxifen-Inducible iCre-Driver Mouse Strains for Eye and Brain, Including Seventeen Carrying a New Inducible-First Constitutive-Ready Allele

To understand gene function, the cre/loxP conditional system is the most powerful available for temporal and spatial control of expression in mouse. However, the research community requires more cre recombinase expressing transgenic mouse strains (cre-drivers) that restrict expression to specific cell types. To address these problems, a high-throughput method for large-scale production that produces high-quality results is necessary. Further, endogenous promoters need to be chosen that drive cell type specific expression, or we need to further focus the expression by manipulating the promoter. Here we test the suitability of using knock-ins at the docking site 5′ of Hprt for rapid development of numerous cre-driver strains focused on expression in adulthood, using an improved cre tamoxifen inducible allele (icre/ERT2), and testing a novel inducible-first, constitutive-ready allele (icre/f3/ERT2/f3). In addition, we test two types of promoters either to capture an endogenous expression pattern (MaxiPromoters), or to restrict expression further using minimal promoter element(s) designed for expression in restricted cell types (MiniPromoters). We provide new cre-driver mouse strains with applicability for brain and eye research. In addition, we demonstrate the feasibility and applicability of using the locus 5′ of Hprt for the rapid generation of substantial numbers of cre-driver strains. We also provide a new inducible-first constitutive-ready allele to further speed cre-driver generation. Finally, all these strains are available to the research community through The Jackson Laboratory.

Identifying cis Elements for Spatiotemporal Control of Mammalian DNA Replication

The temporal order of DNA replication (replication timing [RT]) is highly coupled with genome architecture, but cis-elements regulating either remain elusive. We created a series of CRISPR-mediated deletions and inversions of a pluripotency-associated topologically associating domain (TAD) in mouse ESCs. CTCF-associated domain boundaries were dispensable for RT. CTCF protein depletion weakened most TAD boundaries but had no effect on RT or A/B compartmentalization genome-wide. By contrast, deletion of three intra-TAD CTCF-independent 3D contact sites caused a domain-wide early-to-late RT shift, an A-to-B compartment switch, weakening of TAD architecture, and loss of transcription. The dispensability of TAD boundaries and the necessity of these "early replication control elements" (ERCEs) was validated by deletions and inversions at additional domains. Our results demonstrate that discrete cis-regulatory elements orchestrate domain-wide RT, A/B compartmentalization, TAD architecture, and transcription, revealing fundamental principles linking genome structure and function.

Transcription-driven chromatin repression of Intragenic transcription start sites

Progression of RNA polymerase II (RNAPII) transcription relies on the appropriately positioned activities of elongation factors. The resulting profile of factors and chromatin signatures along transcription units provides a “positional information system” for transcribing RNAPII. Here, we investigate a chromatin-based mechanism that suppresses intragenic initiation of RNAPII transcription. We demonstrate that RNAPII transcription across gene promoters represses their function in plants. This repression is characterized by reduced promoter-specific molecular signatures and increased molecular signatures associated with RNAPII elongation. The conserved FACT histone chaperone complex is required for this repression mechanism. Genome-wide Transcription Start Site (TSS) mapping reveals thousands of discrete intragenic TSS positions in fact mutants, including downstream promoters that initiate alternative transcript isoforms. We find that histone H3 lysine 4 mono-methylation (H3K4me1), an Arabidopsis RNAPII elongation signature, is enriched at FACT-repressed intragenic TSSs. Our analyses suggest that FACT is required to repress intragenic TSSs at positions that are in part characterized by elevated H3K4me1 levels. In sum, conserved and plant-specific chromatin features correlate with the co-transcriptional repression of intragenic TSSs. Our insights into TSS repression by RNAPII transcription promise to inform the regulation of alternative transcript isoforms and the characterization of gene regulation through the act of pervasive transcription across eukaryotic genomes.

Genes represent DNA elements that are transcribed into mRNA. However, the position where transcription actually starts can be dynamically regulated to expand the diversity of RNA isoforms produced from a single gene. Functionally, alternative Transcription Start Sites (TSSs) may generate protein isoforms with differing N-terminal regions and distinct cellular functions. In plants, light signaling regulates protein isoforms largely through regulated TSS selection, emphasizing the biological significance of this mechanism. Despite the importance of alternative TSS selection, little is known about the underlying molecular mechanisms. Here, we characterize for the first time how transcription initiation from an upstream promoter represses alternative downstream promoter activity in plants. This repression mechanism is associated with chromatin changes that are required to maintain precise gene expression control. Specific chromatin signatures are established during transcription via dynamic interactions between the transcription machinery and associated factors. The conserved histone chaperone complex FACT is one such factor involved in regulating the chromatin environment along genes during transcription. We find that mutant plants with reduced FACT activity specifically initiate transcription from thousands of intragenic positions, thus expanding RNA isoform diversity. Overall, our study reveals conserved and plant-specific chromatin features associated with the co-transcriptional repression of downstream intragenic TSSs. These findings promise to help inform the molecular mechanism underlying environmentally-triggered TSS regulation in plants.

Extracellular signal-regulated kinase-1 phosphorylates early growth response-1 at serine 26

Egr-1, an immediate-early gene product and master regulator was originally described as a phosphoprotein following its discovery in the 1980s. However specific residue(s) phosphorylated in Egr-1 remain elusive. Here we phosphorylated recombinant Egr-1 in vitro with ERK1 prior to mass spectrometry, which identified phosphorylation of Ser¹² and Ser²⁶ with the latter ∼12 times more abundant than Ser¹². Phosphorylation of wild-type recombinant Egr-1 (as compared with Ser²⁶>Ala²⁶ mutant Egr-1) revealed that Ser²⁶ accounts for the majority of phosphorylation of Egr-1 by ERK1. N-FGSFPH(pS)PTMDNYC-C was used as an antigen to generate mouse monoclonal antibodies (pS26 MAb). pS26 MAb recognised ERK1-phosphorylated Egr-1 but not Egr-1 bearing a point mutation at Ser²⁶. pS26 MAb recognised inducible ∼75 kDa and 100 kDa species in nuclear extracts of cells exposed to FGF-2. Peptide blocking revealed both inducible species were phosphosite-specific. Immunoprecipitation of nuclear extracts of cells exposed to FGF-2 with pS26 MAb followed by SDS-PAGE and mass spectrometry identified Egr-1 sequences corresponding to the ∼75 kDa species but not ∼100 kDa species. This study identifies a specific amino acid phosphorylated in endogenous Egr-1.

C1 CAGE detects transcription start sites and enhancer activity at single-cell resolution

Single-cell transcriptomic profiling is a powerful tool to explore cellular heterogeneity. However, most of these methods focus on the 3′-end of polyadenylated transcripts and provide only a partial view of the transcriptome. We introduce C1 CAGE, a method for the detection of transcript 5′-ends with an original sample multiplexing strategy in the C1^TM microfluidic system. We first quantifiy the performance of C1 CAGE and find it as accurate and sensitive as other methods in the C1 system. We then use it to profile promoter and enhancer activities in the cellular response to TGF-β of lung cancer cells and discover subpopulations of cells differing in their response. We also describe enhancer RNA dynamics revealing transcriptional bursts in subsets of cells with transcripts arising from either strand in a mutually exclusive manner, validated using single molecule fluorescence in situ hybridization.

Single-cell transcriptomic profiling allows the exploration of cellular heterogeneity but commonly focuses on the 3′-end of the transcript. Here the authors introduce C1 CAGE, which detects the 5′ transcript end in a multiplexed microfluidic system.

Antisense Transcription in Loci Associated to Hereditary Neurodegenerative Diseases

Natural antisense transcripts are common features of mammalian genes providing additional regulatory layers of gene expression. A comprehensive description of antisense transcription in loci associated to familial neurodegenerative diseases may identify key players in gene regulation and provide tools for manipulating gene expression. We take advantage of the FANTOM5 sequencing datasets that represent the largest collection to date of genome-wide promoter usage in almost 2000 human samples. Transcription start sites (TSSs) are mapped at high resolution by the use of a modified protocol of cap analysis of gene expression (CAGE) for high-throughput single molecule next-generation sequencing with Helicos (hCAGE). Here we present the analysis of antisense transcription at 17 loci associated to hereditary Alzheimer’s disease, Frontotemporal Dementia, Parkinson’s disease, Amyotrophic Lateral Sclerosis, and Huntington’s disease. We focused our analysis on libraries derived from brain tissues and primary cells. We also screened libraries from total blood and blood cell populations in the quest for peripheral biomarkers of neurodegenerative diseases. We identified 63 robust promoters in antisense orientation to genes associated to familial neurodegeneration. When applying a less stringent cutoff, this number increases to over 400. A subset of these promoters represents alternative TSSs for 24 FANTOM5 annotated long noncoding RNA (lncRNA) genes, in antisense orientation to 13 of the loci analyzed here, while the remaining contribute to the expression of additional transcript variants. Intersection with GWAS studies, sample ontology, and dynamic expression reveals association to specific genetic traits as well as cell and tissue types, not limited to neurodegenerative diseases. Antisense transcription was validated for a subset of genes, including those encoding for Microtubule-Associated Protein Tau, α-synuclein, Parkinsonism-associated deglycase DJ-1, and Leucin-Rich Repeat Kinase 2. This work provides evidence for the existence of additional regulatory mechanisms of the expression of neurodegenerative disease-causing genes by previously not-annotated and/or not-validated antisense long noncoding RNAs.

Regulatory functions of the Mediator kinases CDK8 and CDK19

The Mediator-associated kinases CDK8 and CDK19 function in the context of three additional proteins: CCNC and MED12, which activate CDK8/CDK19 kinase function, and MED13, which enables their association with the Mediator complex. The Mediator kinases affect RNA polymerase II (pol II) transcription indirectly, through phosphorylation of transcription factors and by controlling Mediator structure and function. In this review, we discuss cellular roles of the Mediator kinases and mechanisms that enable their biological functions. We focus on sequence-specific, DNA-binding transcription factors and other Mediator kinase substrates, and how CDK8 or CDK19 may enable metabolic and transcriptional reprogramming through enhancers and chromatin looping. We also summarize Mediator kinase inhibitors and their therapeutic potential. Throughout, we note conserved and divergent functions between yeast and mammalian CDK8, and highlight many aspects of kinase module function that remain enigmatic, ranging from potential roles in pol II promoter-proximal pausing to liquid-liquid phase separation.

Enhancer RNAs: a missing regulatory layer in gene transcription

Enhancers and super-enhancers exert indispensable roles in maintaining cell identity through spatiotemporally regulating gene transcription. Meanwhile, active enhancers and super-enhancers also produce transcripts termed enhancer RNAs (eRNAs) from their DNA elements. Although enhancers have been identified for more than 30 years, widespread transcription from enhancers are just discovered by genome-wide sequencing and considered as the key to understand longstanding questions in gene transcription. RNA-transcribed enhancers are marked by histone modifications such as H3K4m1/2 and H3K27Ac, and enriched with transcription regulatory factors such as LDTFs, P300, CBP, BRD4 and MED1. Those regulatory factors might constitute a Mega-Trans-like complex to potently activate enhancers. Compared to mRNAs, eRNAs are quite unstable and play roles at local. Functionally, it has been shown that eRNAs promote formation of enhancer-promoter loops. Several studies also demonstrated that eRNAs help the binding of RNA polymerase II (RNAPII) or transition of paused RNAPII by de-association of the negative elongation factor (NELF) complex. Nevertheless, these proposed mechanisms are not universally accepted and still under controversy. Here, we comprehensively summarize the reported findings and make perspectives for future exploration. We also believe that super-enhancer derived RNAs (seRNAs) might be informative to understand the nature of super-enhancers.

TELS: A Novel Computational Framework for Identifying Motif Signatures of Transcribed Enhancers

In mammalian cells, transcribed enhancers (TrEns) play important roles in the initiation of gene expression and maintenance of gene expression levels in a spatiotemporal manner. One of the most challenging questions is how the genomic characteristics of enhancers relate to enhancer activities. To date, only a limited number of enhancer sequence characteristics have been investigated, leaving space for exploring the enhancers’ DNA code in a more systematic way. To address this problem, we developed a novel computational framework, Transcribed Enhancer Landscape Search (TELS), aimed at identifying predictive cell type/tissue-specific motif signatures of TrEns. As a case study, we used TELS to compile a comprehensive catalog of motif signatures for all known TrEns identified by the FANTOM5 consortium across 112 human primary cells and tissues. Our results confirm that combinations of different short motifs characterize in an optimized manner cell type/tissue-specific TrEns. Our study is the first to report combinations of motifs that maximize classification performance of TrEns exclusively transcribed in one cell type/tissue from TrEns exclusively transcribed in different cell types/tissues. Moreover, we also report 31 motif signatures predictive of enhancers’ broad activity. TELS codes and material are publicly available at http://www.cbrc.kaust.edu.sa/TELS.

An evolutionary framework for measuring epigenomic information and estimating cell-type-specific fitness consequences

Here we ask the question "How much information do epigenomic datasets provide about human genomic function?" We consider nine epigenomic features across 115 cell types and measure information about function as a reduction in entropy under a probabilistic evolutionary model fitted to human and nonhuman primate genomes. Several epigenomic features yield more information in combination than they do individually. We find that the entropy in human genetic variation predominantly reflects a balance between mutation and neutral drift. Our cell-type-specific FitCons scores reveal relationships among cell types and suggest that around 8% of nucleotide sites are constrained by natural selection.