Divergent transcription: a driving force for new gene origination?

Molecular drivers of human cerebral cortical evolution

One of the most important questions in human evolutionary biology is how our ancestor has acquired an expanded volume of the cerebral cortex, which may have significantly impacted on improving our cognitive abilities. Recent comparative approaches have identified developmental features unique to the human or hominid cerebral cortex, not shared with other animals including conventional experimental models. In addition, genomic, transcriptomic, and epigenomic signatures associated with human- or hominid-specific processes of the cortical development are becoming identified by virtue of technical progress in the deep nucleotide sequencing. This review discusses ontogenic and phylogenetic processes of the human cerebral cortex, followed by the introduction of recent comprehensive approaches identifying molecular mechanisms potentially driving the evolutionary changes in the cortical development.

Transcribe this way: Rap1 confers promoter directionality by repressing divergent transcription

In eukaryotes, divergent transcription is a major source of noncoding RNAs. Recent studies have uncovered that in yeast, the transcription factor Rap1 restricts transcription in the divergent direction and thereby controls promoter directionality. Here, we summarize these findings, propose regulatory principles, and discuss the implications for eukaryotic gene regulation.

P-TEFb Regulates Transcriptional Activation in Non-coding RNA Genes

Many non-coding RNAs (ncRNAs) serve as regulatory molecules in various physiological pathways, including gene expression in mammalian cells. Distinct from protein-coding RNA expression, ncRNA expression is regulated solely by transcription and RNA processing/stability. It is thus important to understand transcriptional regulation in ncRNA genes but is yet to be known completely. Previously, we identified that a subset of mammalian ncRNA genes is transcriptionally regulated by RNA polymerase II (Pol II) promoter-proximal pausing and in a tissue-specific manner. In this study, human ncRNA genes that are expressed in the early G₁ phase, termed immediate early ncRNA genes, were monitored to assess the function of positive transcription elongation factor b (P-TEFb), a master Pol II pausing regulator for protein-coding genes, in ncRNA transcription. Our findings indicate that the expression of many ncRNA genes is induced in the G₀–G₁ transition and regulated by P-TEFb. Interestingly, a biphasic characteristic of P-TEFb-dependent transcription of serum responsive ncRNA genes was observed: Pol II carboxyl-terminal domain phosphorylated at serine 2 (S2) was largely increased in the transcription start site (TSS, -300 to +300) whereas overall, it was decreased in the gene body (GB, > +350) upon chemical inhibition of P-TEFb. In addition, the three representative, immediate early ncRNAs, whose expression is dependent on P-TEFb, metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), nuclear enriched abundant transcript 1 (NEAT1), and X-inactive specific transcript (XIST), were further analyzed for determining P-TEFb association. Taken together, our data suggest that transcriptional activation of many human ncRNAs utilizes the pausing and releasing of Pol II, and that the regulatory mechanism of transcriptional elongation in these genes requires the function of P-TEFb. Furthermore, we propose that ncRNA and mRNA transcription are regulated by similar mechanisms while P-TEFb inhibition unexpectedly increases S2 Pol II phosphorylation in the TSSs in many ncRNA genes.

One Sentence Summary: P-TEFb regulates Pol II phosphorylation for transcriptional activation in many stimulus-inducible ncRNA genes.

Shaping the nebulous enhancer in the era of high-throughput assays and genome editing

Since the 1st discovery of transcriptional enhancers in 1981, their textbook definition has remained largely unchanged in the past 37 years. With the emergence of high-throughput assays and genome editing, which are switching the paradigm from bottom-up discovery and testing of individual enhancers to top-down profiling of enhancer activities genome-wide, it has become increasingly evidenced that this classical definition has left substantial gray areas in different aspects. Here we survey a representative set of recent research articles and report the definitions of enhancers they have adopted. The results reveal that a wide spectrum of definitions is used usually without the definition stated explicitly, which could lead to difficulties in data interpretation and downstream analyses. Based on these findings, we discuss the practical implications and suggestions for future studies.

Regulation of Intronic Polyadenylation by PCF11 Impacts mRNA Expression of Long Genes

Regulation of cleavage and polyadenylation (CPA) affects gene expression and polyadenylation site (PAS) choice. Here, we report that the CPA and termination factor PCF11 modulates gene expression on the basis of gene size. Although downregulation of PCF11 leads to inhibition of short gene expression, long genes are upregulated because of suppressed intronic polyadenylation (IPA) enriched in large introns. We show that this regulatory scheme, named PCF11-mediated expression regulation through IPA (PEIPA), takes place in cell differentiation, during which downregulation of PCF11 is coupled with upregulation of long genes with functions in cell morphology, adhesion, and migration. PEIPA targets distinct gene sets in different cell contexts with similar rules. Furthermore, PCF11 is autoregulated through a conserved IPA site, the removal of which leads to global activation of PASs close to gene promotors. Therefore, PCF11 uses distinct mechanisms to regulate genes of different sizes, and its autoregulation maintains homeostasis of PAS usage in the cell.

A Molecular Portrait of De Novo Genes in Yeasts

New genes, with novel protein functions, can evolve "from scratch" out of intergenic sequences. These de novo genes can integrate the cell's genetic network and drive important phenotypic innovations. Therefore, identifying de novo genes and understanding how the transition from noncoding to coding occurs are key problems in evolutionary biology. However, identifying de novo genes is a difficult task, hampered by the presence of remote homologs, fast evolving sequences and erroneously annotated protein coding genes. To overcome these limitations, we developed a procedure that handles the usual pitfalls in de novo gene identification and predicted the emergence of 703 de novo gene candidates in 15 yeast species from 2 genera whose phylogeny spans at least 100 million years of evolution. We validated 85 candidates by proteomic data, providing new translation evidence for 25 of them through mass spectrometry experiments. We also unambiguously identified the mutations that enabled the transition from noncoding to coding for 30 Saccharomyces de novo genes. We established that de novo gene origination is a widespread phenomenon in yeasts, only a few being ultimately maintained by selection. We also found that de novo genes preferentially emerge next to divergent promoters in GC-rich intergenic regions where the probability of finding a fortuitous and transcribed ORF is the highest. Finally, we found a more than 3-fold enrichment of de novo genes at recombination hot spots, which are GC-rich and nucleosome-free regions, suggesting that meiotic recombination contributes to de novo gene emergence in yeasts.

Cis- and trans-factors affecting AID targeting and mutagenic outcomes in antibody diversification

Antigen receptor diversification is a hallmark of adaptive immunity which allows specificity of the receptor to particular antigen. B cell receptor (BCR) or its secreted form, antibody, is diversified through antigen-independent and antigen-dependent mechanisms. During B cell development in bone marrow, BCR is diversified via V(D)J recombination mediated by RAG endonuclease. Upon stimulation by antigen, B cell undergo somatic hypermutation (SHM) to allow affinity maturation and class switch recombination (CSR) to change the effector function of the antibody. Both SHM and CSR are initiated by activation-induced cytidine deaminase (AID). Repair of AID-initiated lesions through different DNA repair pathways results in diverse mutagenic outcomes. Here, we focus on discussing cis- and trans-factors that target AID to its substrates and factors that affect different outcomes of AID-initiated lesions. The knowledge of mechanisms that govern AID targeting and outcomes could be harnessed to elicit rare functional antibodies and develop ex vivo antibody diversification approaches with diversifying base editors.

Transcriptional activation of a chimeric retrogene PIPSL in a hominoid ancestor

Retrogenes are a class of functional genes derived from the mRNA of various intron-containing genes. PIPSL was created through a unique mechanism, whereby distinct genes were assembled at the RNA level, and the resulting chimera was then reverse transcribed and integrated into the genome by the L1 retrotransposon. Expression of PIPSL RNA via its transcription start sites (TSSs) has been confirmed in the testes of humans and chimpanzee. Here, we demonstrated that PIPSL RNA is expressed in the testis of the white-handed gibbon. The 5'-end positions of gibbon RNAs were confined to a narrow range upstream of the PIPSL start codon and overlapped with those of orangutan and human, suggesting that PIPSL TSSs are similar among hominoid species. Reporter assays using a luciferase gene and the flanking sequences of human PIPSL showed that an upstream sequence exhibits weak promoter activity in human cells. Our findings suggest that PIPSL might have acquired a promoter at an early stage of hominoid evolution before the divergence of gibbons and ultimately retained similar TSSs in all of the lineages. Moreover, the upstream sequence derived from the phosphatidylinositol-4-phosphate 5-kinase, type I, alpha 5' untranslated region and/or neighboring repetitive sequences in the genome possibly exhibits promoter activity. Furthermore, we observed that a TATA-box-like sequence has emerged by nucleotide substitution in a lineage leading to humans, with this possibly responsible for a broader distribution of the human PIPSL TSSs.

Production of Spliced Long Noncoding RNAs Specifies Regions with Increased Enhancer Activity

Active enhancers in mammals produce enhancer RNAs (eRNAs) that are bidirectionally transcribed, unspliced, and unstable. Enhancer regions are also enriched with long noncoding RNA (lncRNA) transcripts, which are typically spliced and substantially more stable. In order to explore the relationship between these two classes of RNAs, we analyzed DNase hypersensitive sites with evidence of bidirectional transcription, which we termed eRNA-producing centers (EPCs). EPCs found very close to transcription start sites of lncRNAs exhibit attributes of both enhancers and promoters, including distinctive DNA motifs and a characteristic chromatin landscape. These EPCs are associated with higher enhancer activity, driven at least in part by the presence of conserved, directional splicing signals that promote lncRNA production, pointing at a causal role of lncRNA processing in enhancer activity. Together, our results suggest that the conserved ability of some enhancers to produce lncRNAs augments their activity in a manner likely mediated through lncRNA maturation.

Divergent lncRNA GATA3-AS1 Regulates GATA3 Transcription in T-Helper 2 Cells

Long non-coding RNAs (lncRNAs) possess a diverse array of regulatory functions including activation and silencing of gene transcription, regulation of splicing, and coordinating epigenetic modifications. GATA3-AS1 is a divergent lncRNA gene neighboring GATA3. GATA3 is considered the master regulator of TH2 lineage commitment enabling TH2 effector cells to efficiently transcribe genes encoding cytokines IL-4, IL-5, and IL-13. Here, we show that the GATA3-AS1 lncRNA is selectively expressed under TH2 polarizing conditions and is necessary for efficient transcription of GATA3, IL5, and IL13 genes, while being sufficient for GATA3 transcription. GATA3-AS1 is required for formation of permissive chromatin marks, H3K27 acetylation and H3K4 di/tri-methylation, at the GATA3-AS1-GATA3 locus. Further, GATA3-AS1 binds components of the MLL methyltransferase and forms a DNA-RNA hybrid (R-loop) thus tethering the MLL methyltransferase to the gene locus. Our results indicate a novel regulatory function for a divergent lncRNA and provide new insight into the function of lncRNAs in T helper cell differentiation.

Repurposing of promoters and enhancers during mammalian evolution

Promoters and enhancers-key controllers of gene expression-have long been distinguished from each other based on their function. However, recent work suggested that common architectural and functional features might have facilitated the conversion of one type of element into the other during evolution. Here, based on cross-mammalian analyses of epigenome and transcriptome data, we provide support for this hypothesis by detecting 445 regulatory elements with signatures of activity turnover (termed P/E elements). Most events represent transformations of putative ancestral enhancers into promoters, leading to the emergence of species-specific transcribed loci or 5' exons. Distinct GC sequence compositions and stabilizing 5' splicing (U1) regulatory motif patterns may have predisposed P/E elements to regulatory repurposing, and changes in the U1 and polyadenylation signal densities and distributions likely drove the evolutionary activity switches. Our work suggests that regulatory repurposing facilitated regulatory innovation and the origination of new genes and exons during evolution.

Functional classification of long non-coding RNAs by k-mer content

The functions of most long non-coding RNAs (lncRNAs) are unknown. In contrast to proteins, lncRNAs with similar functions often lack linear sequence homology; thus, the identification of function in one lncRNA rarely informs the identification of function in others. We developed a sequence comparison method to deconstruct linear sequence relationships in lncRNAs and evaluate similarity based on the abundance of short motifs called kmers. We found that lncRNAs of related function often had similar kmer profiles despite lacking linear homology, and that kmer profiles correlated with protein binding to lncRNAs and with their subcellular localization. Using a novel assay to quantify Xist-like regulatory potential, we directly demonstrated that evolutionarily unrelated lncRNAs can encode similar function through different spatial arrangements of related sequence motifs. Kmer-based classification is a powerful approach to detect recurrent relationships between sequence and function in lncRNAs.

Long noncoding RNA repertoire and targeting by nuclear exosome, cytoplasmic exonuclease, and RNAi in fission yeast

Long noncoding RNAs (lncRNAs), which are longer than 200 nucleotides but often unstable, contribute a substantial and diverse portion to pervasive noncoding transcriptomes. Most lncRNAs are poorly annotated and understood, although several play important roles in gene regulation and diseases. Here we systematically uncover and analyze lncRNAs in Schizosaccharomyces pombe. Based on RNA-seq data from twelve RNA-processing mutants and nine physiological conditions, we identify 5775 novel lncRNAs, nearly 4× the previously annotated lncRNAs. The expression of most lncRNAs becomes strongly induced under the genetic and physiological perturbations, most notably during late meiosis. Most lncRNAs are cryptic and suppressed by three RNA-processing pathways: the nuclear exosome, cytoplasmic exonuclease, and RNAi. Double-mutant analyses reveal substantial coordination and redundancy among these pathways. We classify lncRNAs by their dominant pathway into cryptic unstable transcripts (CUTs), Xrn1-sensitive unstable transcripts (XUTs), and Dicer-sensitive unstable transcripts (DUTs). XUTs and DUTs are enriched for antisense lncRNAs, while CUTs are often bidirectional and actively translated. The cytoplasmic exonuclease, along with RNAi, dampens the expression of thousands of lncRNAs and mRNAs that become induced during meiosis. Antisense lncRNA expression mostly negatively correlates with sense mRNA expression in the physiological, but not the genetic conditions. Intergenic and bidirectional lncRNAs emerge from nucleosome-depleted regions, upstream of positioned nucleosomes. Our results highlight both similarities and differences to lncRNA regulation in budding yeast. This broad survey of the lncRNA repertoire and characteristics in S. pombe, and the interwoven regulatory pathways that target lncRNAs, provides a rich framework for their further functional analyses.

Control of human gene expression: High abundance of divergent transcription in genes containing both INR and BRE elements in the core promoter

Background

DNA sequence elements in the core promoter can play a central role in regulation of gene expression. Core elements (e.g. INR and TATA box) are located within ~50bp of the transcription start site and both upstream and downstream elements are known. Although all can affect the level of gene expression, their mechanism of action has yet to be fully defined. The studies described here are focused on two core promoter elements, INR and BRE, in the human genome. The locations of the two elements were determined in a large number of human promoters and the results were interpreted in terms of overall promoter function.

Results

A total of 13,406 promoters were collected from the reference version of the human genome and found to contain 62,891 INR sequences and 32,290 BRE. An INR sequence was found in the core region of 1231 (9.2%) promoters and a BRE in 2592 (19.3%); 121 promoters (0.9%) have both INR and BRE elements. Counts support the view that most human promoters lack an INR or BRE element in the core promoter. Further analysis was carried out with the aligned aggregate of promoters from each chromosome. The results showed distinct INR distributions in separate chromosome groups indicating a degree of chromosome specificity to the way core promoter elements are deployed in the genome. The rare promoters with both INR and BRE elements were found to be enriched among the genes with divergent transcription. Enrichment raises the possibility that core promoter elements can have a function in chromosome organization as well as in initiation of transcription.

Transcriptional regulation by promoters with enhancer function

Promoters with enhancer activity have been described recently. In this point of view, we will discuss current findings highlighting the commonality of this type of regulatory elements, their genetic and epigenetic characteristics, their potential biological roles in the regulation of gene expression and the underlining molecular mechanisms.

ABBREVIATIONS

TSS: transcription start site; IFN: interferon; STARR-seq: Self-Transcribing Active Regulatory Region sequencing; MPRA: Massively Parallel Reporter Assay; ChIP: chromatin immunoprecipitation; CRISPR: Clustered Regularly Interspaced Short Palindromic Repeats; lncRNA: long non-coding RNA.

Expression and replication of virus-like circular DNA in human cells

The consumption of bovine milk and meat is considered a risk factor for colon- and breast cancer formation, and milk consumption has also been implicated in an increased risk for developing Multiple Sclerosis (MS). A number of highly related virus-like DNAs have been recently isolated from bovine milk and sera and from a brain sample of a MS patient. As a genetic activity of these Acinetobacter-related bovine milk and meat factors (BMMFs) is unknown in eukaryotes, we analyzed their expression and replication potential in human HEK293TT cells. While all analyzed BMMFs show transcriptional activity, the MS brain isolate MSBI1.176, sharing homology with a transmissible spongiform encephalopathy-associated DNA molecule, is transcribed at highest levels. We show expression of a replication-associated protein (Rep), which is highly conserved among all BMMFs, and serological tests indicate a human anti-Rep immune response. While the cow milk isolate CMI1.252 is replication-competent in HEK293TT cells, replication of MSBI1.176 is complemented by CMI1.252, pointing at an interplay during the establishment of persistence in human cells. Transcriptome profiling upon BMMF expression identified host cellular gene expression changes related to cell cycle progression and cell viability control, indicating potential pathways for a pathogenic involvement of BMMFs.

Widespread transcriptional pausing and elongation control at enhancers

In this study, Henriques et al. demonstrate that transcription is a nearly universal feature of enhancers in Drosophila and mammalian cells and that nascent RNA sequencing strategies are optimal for identification of both enhancers and superenhancers. Their findings provide insights into the unique characteristics of superenhancers, which stimulate high-level gene expression through rapid pause release; interestingly, this property renders associated genes resistant to loss of factors that stabilize paused RNAPII.

Regulation by gene-distal enhancers is critical for cell type-specific and condition-specific patterns of gene expression. Thus, to understand the basis of gene activity in a given cell type or tissue, we must identify the precise locations of enhancers and functionally characterize their behaviors. Here, we demonstrate that transcription is a nearly universal feature of enhancers in Drosophila and mammalian cells and that nascent RNA sequencing strategies are optimal for identification of both enhancers and superenhancers. We dissect the mechanisms governing enhancer transcription and discover remarkable similarities to transcription at protein-coding genes. We show that RNA polymerase II (RNAPII) undergoes regulated pausing and release at enhancers. However, as compared with mRNA genes, RNAPII at enhancers is less stable and more prone to early termination. Furthermore, we found that the level of histone H3 Lys4 (H3K4) methylation at enhancers corresponds to transcriptional activity such that highly active enhancers display H3K4 trimethylation rather than the H3K4 monomethylation considered a hallmark of enhancers. Finally, our work provides insights into the unique characteristics of superenhancers, which stimulate high-level gene expression through rapid pause release; interestingly, this property renders associated genes resistant to the loss of factors that stabilize paused RNAPII.

High-resolution transcription maps reveal the widespread impact of roadblock termination in yeast

Transcription termination delimits transcription units but also plays important roles in limiting pervasive transcription. We have previously shown that transcription termination occurs when elongating RNA polymerase II (RNAPII) collides with the DNA-bound general transcription factor Reb1. We demonstrate here that many different DNA-binding proteins can induce termination by a similar roadblock (RB) mechanism. We generated high-resolution transcription maps by the direct detection of RNAPII upon nuclear depletion of two essential RB factors or when the canonical termination pathways for coding and non-coding RNAs are defective. We show that RB termination occurs genomewide and functions independently of (and redundantly with) the main transcription termination pathways. We provide evidence that transcriptional readthrough at canonical terminators is a significant source of pervasive transcription, which is controlled to a large extent by RB termination. Finally, we demonstrate the occurrence of RB termination around centromeres and tRNA genes, which we suggest shields these regions from RNAPII to preserve their functional integrity.

The Importance of Controlling Transcription Elongation at Coding and Noncoding RNA Loci

Here we discuss current paradigms for how transcription initiation and elongation control are achieved in mammalian cells, and how they differ at protein-coding mRNA genes versus noncoding RNA (ncRNA) loci. We present a model for the function of ncRNAs wherein the act of transcription is regulatory, rather than the ncRNA products themselves. We further describe how the establishment of transcriptionally engaged, but paused, RNA polymerase II impacts chromatin structure around divergent transcription start sites, and how this can influence transcription factor binding and mRNA gene activity in the region.

The functions and unique features of long intergenic non-coding RNA

Long intergenic non-coding RNA (lincRNA) genes have diverse features that distinguish them from mRNA-encoding genes and exercise functions such as remodelling chromatin and genome architecture, RNA stabilization and transcription regulation, including enhancer-associated activity. Some genes currently annotated as encoding lincRNAs include small open reading frames (smORFs) and encode functional peptides and thus may be more properly classified as coding RNAs. lincRNAs may broadly serve to fine-tune the expression of neighbouring genes with remarkable tissue specificity through a diversity of mechanisms, highlighting our rapidly evolving understanding of the non-coding genome.

Spurious transcription and its impact on cell function

Most RNA polymerases can initiate transcription from diverse DNA template sequences with relatively few outright sequence restraints. Recent reports have demonstrated that failure to subdue the promiscuity of RNA polymerase in vivo can severely impede cell function. This phenomenon appears common to all cell types with undesirable effects ranging from growth inhibition in prokaryotes to cancer in higher organisms. Here we discuss similarities and differences in strategies employed by cells to minimise spurious transcription across life's domains.

Not Just Noise: Genomics and Genetics Bring Long Noncoding RNAs into Focus

Thousands of long noncoding RNAs (lncRNAs) have been annotated, yet their functions remain unclear. Recent studies-including Schlackow et al. (2017) in this issue of Molecular Cell-using orthogonal methods investigated the expression and functions of lncRNAs, resulting in deeper appreciation for the salient differences between lncRNAs and mRNAs and the roles lncRNAs serve.

RPA Interacts with HIRA and Regulates H3.3 Deposition at Gene Regulatory Elements in Mammalian Cells

The histone chaperone HIRA is involved in depositing histone variant H3.3 into distinct genic regions, including promoters, enhancers, and gene bodies. However, how HIRA deposits H3.3 to these regions remains elusive. Through a short hairpin RNA (shRNA) screening, we identified single-stranded DNA binding protein replication protein A (RPA) as a regulator of the deposition of newly synthesized H3.3 into chromatin. We show that RPA physically interacts with HIRA to form RPA-HIRA-H3.3 complexes, and it co-localizes with HIRA and H3.3 at gene promoters and enhancers. Depletion of RPA1, the largest subunit of the RPA complex, dramatically reduces both HIRA association with chromatin and the deposition of newly synthesized H3.3 at promoters and enhancers and leads to altered transcription at gene promoters. These results support a model whereby RPA, best known for its role in DNA replication and repair, recruits HIRA to promoters and enhancers and regulates deposition of newly synthesized H3.3 to these regulatory elements for gene regulation.

Gene regulation of mammalian long non-coding RNA

RNA polymerase II (Pol II) transcribes two classes of RNAs, protein-coding and non-protein-coding (ncRNA) genes. ncRNAs are also synthesized by RNA polymerases I and III (Pol I and III). In humans, the number of ncRNA genes exceeds more than twice that of protein-coding genes. However, the history of studying Pol II-synthesized ncRNA is relatively short. Since early 2000s, important biological and pathological functions of these ncRNA genes have begun to be discovered and intensively studied. And transcription mechanisms of long non-coding RNA (lncRNA) have been recently reported. Transcription of lncRNAs utilizes some transcription factors and mechanisms shared in that of protein-coding genes. In addition, tissue specificity in lncRNA gene expression has been shown. LncRNAs play essential roles in regulating the expression of neighboring or distal genes through different mechanisms. This leads to the implication of lncRNAs in a wide variety of biological pathways and pathological development. In this review, the newly discovered transcription mechanisms, characteristics, and functions of lncRNA are discussed.

The R-loop is a common chromatin feature of the Arabidopsis genome

R-loops are functional structures in chromatin comprising one single-stranded DNA and a DNA:RNA hybrid. Here, we report ssDRIP-seq, a single-strand DNA ligation-based library preparation technique for genome-wide identification of R-loops. When applied in Arabidopsis, ssDRIP-seq exhibits high efficiency, low bias and strand specificity. We found that Arabidopsis R-loops are enriched by both AT and GC skews, and are formed in the sense and antisense orientations. R-loops are strongly enriched in gene promoters and gene bodies, and are highly associated with noncoding RNA and repetitive genomic regions. Furthermore, R-loops are negatively correlated with CG DNA hypermethylation, and are prevalent in regions with multiple chromatin modifications, showing strong correlations with the activated and repressed gene loci. Our analyses indicate that R-loops are common features in the Arabidopsis genome and suggest that the R-loops play diverse roles in genome organization and gene regulation, thereby providing insights into plant nuclear genome formation and function.

A subset of conserved mammalian long non-coding RNAs are fossils of ancestral protein-coding genes

Background

Only a small portion of human long non-coding RNAs (lncRNAs) appear to be conserved outside of mammals, but the events underlying the birth of new lncRNAs in mammals remain largely unknown. One potential source is remnants of protein-coding genes that transitioned into lncRNAs.

Results

We systematically compare lncRNA and protein-coding loci across vertebrates, and estimate that up to 5% of conserved mammalian lncRNAs are derived from lost protein-coding genes. These lncRNAs have specific characteristics, such as broader expression domains, that set them apart from other lncRNAs. Fourteen lncRNAs have sequence similarity with the loci of the contemporary homologs of the lost protein-coding genes. We propose that selection acting on enhancer sequences is mostly responsible for retention of these regions. As an example of an RNA element from a protein-coding ancestor that was retained in the lncRNA, we describe in detail a short translated ORF in the JPX lncRNA that was derived from an upstream ORF in a protein-coding gene and retains some of its functionality.

Conclusions

We estimate that ~ 55 annotated conserved human lncRNAs are derived from parts of ancestral protein-coding genes, and loss of coding potential is thus a non-negligible source of new lncRNAs. Some lncRNAs inherited regulatory elements influencing transcription and translation from their protein-coding ancestors and those elements can influence the expression breadth and functionality of these lncRNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1293-0) contains supplementary material, which is available to authorized users.

The punctilious RNA polymerase II core promoter

This review by Vo ngoc et al. expands the view of the RNA polymerase II core promoter, which is comprised of classical DNA sequence motifs, sequence-specific DNA-binding transcription factors, chromatin signals, and DNA structure.

The signals that direct the initiation of transcription ultimately converge at the core promoter, which is the gateway to transcription. Here we provide an overview of the RNA polymerase II core promoter in bilateria (bilaterally symmetric animals). The core promoter is diverse in terms of its composition and function yet is also punctilious, as it acts with strict rules and precision. We additionally describe an expanded view of the core promoter that comprises the classical DNA sequence motifs, sequence-specific DNA-binding transcription factors, chromatin signals, and DNA structure. This model may eventually lead to a more unified conceptual understanding of the core promoter.

Pervasive coexpression of spatially proximal genes is buffered at the protein level

Genes are not randomly distributed in the genome. In humans, 10% of protein-coding genes are transcribed from bidirectional promoters and many more are organised in larger clusters. Intriguingly, neighbouring genes are frequently coexpressed but rarely functionally related. Here we show that coexpression of bidirectional gene pairs, and closeby genes in general, is buffered at the protein level. Taking into account the 3D architecture of the genome, we find that co-regulation of spatially close, functionally unrelated genes is pervasive at the transcriptome level, but does not extend to the proteome. We present evidence that non-functional mRNA coexpression in human cells arises from stochastic chromatin fluctuations and direct regulatory interference between spatially close genes. Protein-level buffering likely reflects a lack of coordination of post-transcriptional regulation of functionally unrelated genes. Grouping human genes together along the genome sequence, or through long-range chromosome folding, is associated with reduced expression noise. Our results support the hypothesis that the selection for noise reduction is a major driver of the evolution of genome organisation.

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

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

Heat Shock Factor 1 (HSF1) is a regulator of stress-induced transcription. Here, the authors investigate changes to transcription and chromatin organization upon stress and find that activated HSF1 binds to transcription-primed promoters and enhancers, and to CTCF occupied, untranscribed chromatin.

The Ground State and Evolution of Promoter Region Directionality

Eukaryotic promoter regions are frequently divergently transcribed in vivo, but it is unknown whether the resultant antisense RNAs are a mechanistic by-product of RNA polymerase II (Pol II) transcription or biologically meaningful. Here, we use a functional evolutionary approach that involves nascent transcript mapping in S. cerevisiae strains containing foreign yeast DNA. Promoter regions in foreign environments lose the directionality they have in their native species. Strikingly, fortuitous promoter regions arising in foreign DNA produce equal transcription in both directions, indicating that divergent transcription is a mechanistic feature that does not imply a function for these transcripts. Fortuitous promoter regions arising during evolution promote bidirectional transcription and over time are purged through mutation or retained to enable new functionality. Similarly, human transcription is more bidirectional at newly evolved enhancers and promoter regions. Thus, promoter regions are intrinsically bidirectional and are shaped by evolution to bias transcription toward coding versus non-coding RNAs.

GT-rich promoters can drive RNA pol II transcription and deposition of H2A.Z in African trypanosomes

Genome‐wide transcription studies are revealing an increasing number of “dispersed promoters” that, unlike “focused promoters”, lack well‐conserved sequence motifs and tight regulation. Dispersed promoters are nevertheless marked by well‐defined chromatin structures, suggesting that specific sequence elements must exist in these unregulated promoters. Here, we have analyzed regions of transcription initiation in the eukaryotic parasite Trypanosoma brucei, in which RNA polymerase II transcription initiation occurs over broad regions without distinct promoter motifs and lacks regulation. Using a combination of site‐specific and genome‐wide assays, we identified GT‐rich promoters that can drive transcription and promote the targeted deposition of the histone variant H2A.Z in a genomic context‐dependent manner. In addition, upon mapping nucleosome occupancy at high resolution, we find nucleosome positioning to correlate with RNA pol II enrichment and gene expression, pointing to a role in RNA maturation. Nucleosome positioning may thus represent a previously unrecognized layer of gene regulation in trypanosomes. Our findings show that even highly dispersed, unregulated promoters contain specific DNA elements that are able to induce transcription and changes in chromatin structure.

From intronization to intron loss: How the interplay between mRNA-associated processes can shape the architecture and the expression of eukaryotic genes

Transcription-coupled processes such as capping, splicing, and cleavage/polyadenylation participate in the journey from genes to proteins. Although they are traditionally thought to serve only as steps in the generation of mature mRNAs, a synthesis of available data indicates that these processes could also act as a driving force for the evolution of eukaryotic genes. A theoretical framework for how mRNA-associated processes may shape gene structure and expression has recently been proposed. Factors that promote splicing and cleavage/polyadenylation in this framework compete for access to overlapping or neighboring signals throughout the transcription cycle. These antagonistic interactions allow mechanisms for intron gain and splice site recognition as well as common trends in eukaryotic gene structure and expression to be coherently integrated. Here, I extend this framework further. Observations that largely (but not exclusively) revolve around the formation of DNA-RNA hybrid structures, called R loops, and promoter directionality are integrated. Additionally, the interplay between splicing factors and cleavage/polyadenylation factors is theorized to also affect the formation of intragenic DNA double-stranded breaks thereby contributing to intron loss. The most notable prediction in this proposition is that RNA molecules can mediate intron loss by serving as a template to repair DNA double-stranded breaks. The framework presented here leverages a vast body of empirical observations, logically extending previous suggestions, and generating verifiable predictions to further substantiate the view that the intracellular environment plays an active role in shaping the structure and the expression of eukaryotic genes.

Primate-specific Long Non-coding RNAs and MicroRNAs

Non-coding RNAs (ncRNAs) are critical regulators of gene expression in essentially all life forms. Long ncRNAs (lncRNAs) and microRNAs (miRNAs) are two important RNA classes possessing regulatory functions. Up to date, many primate-specific ncRNAs have been identified and investigated. Their expression specificity to primate lineage suggests primate-specific roles. It is thus critical to elucidate the biological significance of primate or even human-specific ncRNAs, and to develop potential ncRNA-based therapeutics. Here, we have summarized the studies regarding regulatory roles of some key primate-specific lncRNAs and miRNAs.

Tandem duplications lead to novel expression patterns through exon shuffling in Drosophila yakuba

One common hypothesis to explain the impacts of tandem duplications is that whole gene duplications commonly produce additive changes in gene expression due to copy number changes. Here, we use genome wide RNA-seq data from a population sample of Drosophila yakuba to test this ‘gene dosage’ hypothesis. We observe little evidence of expression changes in response to whole transcript duplication capturing 5′ and 3′ UTRs. Among whole gene duplications, we observe evidence that dosage sharing across copies is likely to be common. The lack of expression changes after whole gene duplication suggests that the majority of genes are subject to tight regulatory control and therefore not sensitive to changes in gene copy number. Rather, we observe changes in expression level due to both shuffling of regulatory elements and the creation of chimeric structures via tandem duplication. Additionally, we observe 30 de novo gene structures arising from tandem duplications, 23 of which form with expression in the testes. Thus, the value of tandem duplications is likely to be more intricate than simple changes in gene dosage. The common regulatory effects from chimeric gene formation after tandem duplication may explain their contribution to genome evolution.

The enclosed work shows that whole gene duplications rarely affect gene expression, in contrast to widely held views that the adaptive value of duplicate genes is related to additive changes in gene expression due to gene copy number. We further explain how tandem duplications that create shuffled gene structures can force upregulation of gene sequences, de novo gene creation, and multifold changes in transcript levels. These results show that tandem duplications can produce new genes that are a source of immediate novelty associated with more extreme expression changes than previously suggested by theory. Further, these gene expression changes are a potential source of both beneficial and pathogenic mutations, immediately relevant to clinical and medical genetics in humans and other metazoans.

Underrepresentation of active histone modification marks in evolutionarily young genes

It is known that evolutionarily new genes can rapidly evolve essential roles in fundamental biological processes. Nevertheless, the underlying molecular mechanism of how they acquire their novel transcriptional pattern is less characterized except for the role of cis-regulatory evolution. Epigenetic modification offers an alternative possibility. Here, we examined how histone modifications have changed among different gene age groups in Drosophila melanogaster by integrative analyses of an updated new gene dataset and published epigenomic data. We found a robust pattern across various datasets where both the coverage and intensity of active histone modifications, histone 3 lysine 4 trimethylation and lysine 36 trimethylation, increased with evolutionary age. Such a temporal correlation is negative and much weaker for the repressive histone mark, lysine 9 trimethylation, which is expected given its major association with heterochromatin. By further comparison with neighboring old genes, the depletion of active marks of new genes could be only partially explained by the local epigenetic context. All these data are consistent with the observation that older genes bear relatively higher expression levels and suggest that the evolution of histone modifications could be implicated in transcriptional evolution after gene birth.

Integrative classification of human coding and noncoding genes through RNA metabolism profiles

Pervasive transcription of the human genome results in a heterogeneous mix of coding RNAs and long noncoding RNAs (lncRNAs). Only a small fraction of lncRNAs have demonstrated regulatory functions, thus making functional lncRNAs difficult to distinguish from nonfunctional transcriptional byproducts. This difficulty has resulted in numerous competing human lncRNA classifications that are complicated by a steady increase in the number of annotated lncRNAs. To address these challenges, we quantitatively examined transcription, splicing, degradation, localization and translation for coding and noncoding human genes. We observed that annotated lncRNAs had lower synthesis and higher degradation rates than mRNAs and discovered mechanistic differences explaining slower lncRNA splicing. We grouped genes into classes with similar RNA metabolism profiles, containing both mRNAs and lncRNAs to varying extents. These classes exhibited distinct RNA metabolism, different evolutionary patterns and differential sensitivity to cellular RNA-regulatory pathways. Our classification provides an alternative to genomic context-driven annotations of lncRNAs.

Principles for RNA metabolism and alternative transcription initiation within closely spaced promoters

Mammalian transcriptomes are complex and formed by extensive promoter activity. In addition, gene promoters are largely divergent and initiate transcription of reverse-oriented promoter upstream transcripts (PROMPTs). Although PROMPTs are commonly terminated early, influenced by polyadenylation sites, promoters often cluster so that the divergent activity of one might impact another. Here we found that the distance between promoters strongly correlates with the expression, stability and length of their associated PROMPTs. Adjacent promoters driving divergent mRNA transcription support PROMPT formation, but owing to polyadenylation site constraints, these transcripts tend to spread into the neighboring mRNA on the same strand. This mechanism to derive new alternative mRNA transcription start sites (TSSs) is also evident at closely spaced promoters supporting convergent mRNA transcription. We suggest that basic building blocks of divergently transcribed core promoter pairs, in combination with the wealth of TSSs in mammalian genomes, provide a framework with which evolution shapes transcriptomes.

Open questions in the study of de novo genes: what, how and why

The study of de novo protein-coding genes is maturing from the ad hoc reporting of individual cases to the systematic analysis of extensive genomic data from several species. We identify three key challenges for this emerging field: understanding how best to identify de novo genes, how they arise and why they spread. We highlight the intellectual challenges of understanding how a de novo gene becomes integrated into pre-existing functions and becomes essential. We suggest that, as with protein sequence evolution, antagonistic co-evolution may be key to de novo gene evolution, particularly for new essential genes and new cancer-associated genes.

R-loop: an emerging regulator of chromatin dynamics

The dynamic structure of chromatin, which exists in two conformational states: heterochromatin and euchromatin, alters the accessibility of the DNA to regulatory factors during transcription, replication, recombination, and DNA damage repair. Chemical modifications of histones and DNA, as well as adenosine triphospahate-dependent nucleosome remodeling, have been the major focus of research on chromatin dynamics over the past two decades. However, recent studies using a DNA-RNA hybrid-specific antibody and next-generation sequencing approaches have revealed that the formation of R-loops, one of the most common non-canonical DNA structures, is an emerging regulator of chromatin states. This review focuses on recent insights into the interplay between R-loop formation and the epigenetic modifications of chromatin in normal and disease states.

The Role of Transcription Factors at Antisense-Expressing Gene Pairs in Yeast

Genes encoded close to one another on the chromosome are often coexpressed, by a mechanism and regulatory logic that remain poorly understood. We surveyed the yeast genome for tandem gene pairs oriented tail-to-head at which expression antisense to the upstream gene was conserved across species. The intergenic region at most such tandem pairs is a bidirectional promoter, shared by the downstream gene mRNA and the upstream antisense transcript. Genomic analyses of these intergenic loci revealed distinctive patterns of transcription factor regulation. Mutation of a given transcription factor verified its role as a regulator in trans of tandem gene pair loci, including the proximally initiating upstream antisense transcript and downstream mRNA and the distally initiating upstream mRNA. To investigate cis-regulatory activity at such a locus, we focused on the stress-induced NAD(P)H dehydratase YKL151C and its downstream neighbor, the metabolic enzyme GPM1. Previous work has implicated the region between these genes in regulation of GPM1 expression; our mutation experiments established its function in rich medium as a repressor in cis of the distally initiating YKL151C sense RNA, and an activator of the proximally initiating YKL151C antisense RNA. Wild-type expression of all three transcripts required the transcription factor Gcr2. Thus, at this locus, the intergenic region serves as a focal point of regulatory input, driving antisense expression and mediating the coordinated regulation of YKL151C and GPM1. Together, our findings implicate transcription factors in the joint control of neighboring genes specialized to opposing conditions and the antisense transcripts expressed between them.

Long Terminal Repeats: From Parasitic Elements to Building Blocks of the Transcriptional Regulatory Repertoire

The life cycle of endogenous retroviruses (ERVs), also called long terminal repeat (LTR) retrotransposons, begins with transcription by RNA polymerase II followed by reverse transcription and re-integration into the host genome. While most ERVs are relics of ancient integration events, "young" proviruses competent for retrotransposition-found in many mammals, but not humans-represent an ongoing threat to host fitness. As a consequence, several restriction pathways have evolved to suppress their activity at both transcriptional and post-transcriptional stages of the viral life cycle. Nevertheless, accumulating evidence has revealed that LTR sequences derived from distantly related ERVs have been exapted as regulatory sequences for many host genes in a wide range of cell types throughout mammalian evolution. Here, we focus on emerging themes from recent studies cataloging the diversity of ERV LTRs acting as important transcriptional regulatory elements in mammals and explore the molecular features that likely account for LTR exaptation in developmental and tissue-specific gene regulation.

Enhancement of Transcription by a Splicing-Competent Intron Is Dependent on Promoter Directionality

Enhancement of transcription by a splicing-competent intron is an evolutionarily conserved feature among eukaryotes. The molecular mechanism underlying the phenomenon, however, is not entirely clear. Here we show that the intron is an important regulator of promoter directionality. Employing strand-specific transcription run-on (TRO) analysis, we show that the transcription of mRNA is favored over the upstream anti-sense transcripts (uaRNA) initiating from the promoter in the presence of an intron. Mutation of either the 5′ or 3′ splice site resulted in the reversal of promoter directionality, thereby suggesting that it is not merely the 5′ splice site but the entire splicing-competent intron that regulates transcription directionality. ChIP analysis revealed the recruitment of termination factors near the promoter region in the presence of an intron. Removal of intron or the mutation of splice sites adversely affected the promoter localization of termination factors. We have earlier demonstrated that the intron-mediated enhancement of transcription is dependent on gene looping. Here we show that gene looping is crucial for the recruitment of termination factors in the promoter-proximal region of an intron-containing gene. In a looping-defective mutant, despite normal splicing, the promoter occupancy of factors required for poly(A)-dependent termination of transcription was compromised. This was accompanied by a concomitant loss of transcription directionality. On the basis of these results, we propose that the intron-dependent gene looping places the terminator-bound factors in the vicinity of the promoter region for termination of the promoter-initiated upstream antisense transcription, thereby conferring promoter directionality.

Eukaryotic genes differ from their prokaryotic counterparts in having intervening non-coding sequences called introns. The precise biological role of introns in eukaryotic systems remains unclear even more than forty years after their initial discovery. One function of intron that has been remarkably conserved during evolution is their ability to enhance the transcription of genes that harbor them. How does the intron regulate transcription, however, is not known. Here we show that the intron enhances gene expression by affecting direction of the promoter-initiated transcription. In the presence of an intron, polymerase tends to transcribe the downstream coding region producing mRNA, while in the absence of a splicing-competent intron polymerase starts transcribing promoter upstream region producing upstream antisense RNA (uaRNA). Intron-mediated promoter directionality was dependent on gene looping, which is the interaction off the promoter and terminator region of a gene in a transcription-dependent manner. We show that the intron-dependent gene looping facilitates the recruitment of termination factors in the promoter-proximal region. The recruited termination factors stop uaRNA synthesis thereby conferring directionality to the promoter-bound polymerase.

Peak-valley-peak pattern of histone modifications delineates active regulatory elements and their directionality

Formation of nucleosome free region (NFR) accompanied by specific histone modifications at flanking nucleosomes is an important prerequisite for enhancer and promoter activity. Due to this process, active regulatory elements often exhibit a distinct shape of histone signal in the form of a peak-valley-peak (PVP) pattern. However, different features of PVP patterns and their robustness in predicting active regulatory elements have never been systematically analyzed. Here, we present PARE, a novel computational method that systematically analyzes the H3K4me1 or H3K4me3 PVP patterns to predict NFRs. We show that NFRs predicted by H3K4me1 and me3 patterns are associated with active enhancers and promoters, respectively. Furthermore, asymmetry in the height of peaks flanking the central valley can predict the directionality of stable transcription at promoters. Using PARE on ChIP-seq histone modifications from four ENCODE cell lines and four hematopoietic differentiation stages, we identified several enhancers whose regulatory activity is stage specific and correlates positively with the expression of proximal genes in a particular stage. In conclusion, our results demonstrate that PVP patterns delineate both the histone modification landscape and the transcriptional activities governed by active enhancers and promoters, and therefore can be used for their prediction. PARE is freely available at http://servers.binf.ku.dk/pare.

Unique features of long non-coding RNA biogenesis and function

Long non-coding RNAs (lncRNAs) are a diverse class of RNAs that engage in numerous biological processes across every branch of life. Although initially discovered as mRNA-like transcripts that do not encode proteins, recent studies have revealed features of lncRNAs that further distinguish them from mRNAs. In this Review, we describe special events in the lifetimes of lncRNAs - before, during and after transcription - and discuss how these events ultimately shape the unique characteristics and functional roles of lncRNAs.

Divergent lncRNAs Regulate Gene Expression and Lineage Differentiation in Pluripotent Cells

Divergent lncRNAs that are transcribed in the opposite direction to nearby protein-coding genes comprise a significant proportion (∼20%) of total lncRNAs in mammalian genomes. Through genome-wide analysis, we found that the distribution of this lncRNA class strongly correlates with essential developmental regulatory genes. In pluripotent cells, divergent lncRNAs regulate the transcription of nearby genes. As an example, the divergent lncRNA Evx1as promotes transcription of its neighbor gene, EVX1, and regulates mesendodermal differentiation. At a single-cell level, early broad expression of Evx1as is followed by a rapid, high-level transcription of EVX1, supporting the idea that Evx1as plays an upstream role to facilitate EVX1 transcription. Mechanistically, Evx1as RNA binds to regulatory sites on chromatin, promotes an active chromatin state, and interacts with Mediator. Based on our analyses, we propose that the biological function of thousands of uncharacterized lncRNAs of this class may be inferred from the role of their neighboring adjacent genes.

Enhancers as non-coding RNA transcription units: recent insights and future perspectives

Networks of regulatory enhancers dictate distinct cell identities and cellular responses to diverse signals by instructing precise spatiotemporal patterns of gene expression. However, 35 years after their discovery, enhancer functions and mechanisms remain incompletely understood. Intriguingly, recent evidence suggests that many, if not all, functional enhancers are themselves transcription units, generating non-coding enhancer RNAs. This observation provides a fundamental insight into the inter-regulation between enhancers and promoters, which can both act as transcription units; it also raises crucial questions regarding the potential biological roles of the enhancer transcription process and non-coding enhancer RNAs. Here, we review research progress in this field and discuss several important, unresolved questions regarding the roles and mechanisms of enhancers in gene regulation.

Bilaterian-like promoters in the highly compact Amphimedon queenslandica genome

The regulatory systems underlying animal development must have evolved prior to the emergence of eumetazoans (cnidarians and bilaterians). Although representatives of earlier-branching animals - sponges ctenophores and placozoans - possess most of the developmental transcription factor families present in eumetazoans, the DNA regulatory elements that these transcription factors target remain uncharted. Here we characterise the core promoter sequences, U1 snRNP-binding sites (5' splice sites; 5'SSs) and polyadenylation sites (PASs) in the sponge Amphimedon queenslandica. Similar to unicellular opisthokonts, Amphimedon's genes are tightly packed in the genome and have small introns. In contrast, its genes possess metazoan-like core promoters populated with binding motifs previously deemed to be specific to vertebrates, including Nrf-1 and Krüppel-like elements. Also as in vertebrates, Amphimedon's PASs and 5'SSs are depleted downstream and upstream of transcription start sites, respectively, consistent with non-elongating transcripts being short-lived; PASs and 5'SSs are more evenly distributed in bidirectional promoters in Amphimedon. The presence of bilaterian-like regulatory DNAs in sponges is consistent with these being early and essential innovations of the metazoan gene regulatory repertoire.