Genome-wide RNA polymerase II: not genes only!

Early in an SV40 infection, histone modifications correlate with the presence or absence of RNAPII and direction of transcription

Since epigenetic regulation seemed likely to be involved in SV40 early transcription following infection, we have analyzed the organization of nucleosomes carrying histone modifications (acetyl-H3, acetyl-H4, H3K9me1, H3K9me3, H3K4me1, H3K4me3, H3K27me3, H4K20me1) at 30 min and 2 h post infection in SV40 minichromosomes prepared in the absence or presence of the transcription inhibitor dichloro-1-beta-d-ribofuranosyl benzimidazole. The former condition was used to determine how SV40 chromatin structure changed during early transcription, and the latter was used to determine the role of active transcription. The location of RNAPII was used as a marker to identify where histone modifications were most likely to be involved in regulation. Acetyl-H3 acted like epigenetic memory by being present at sites subsequently bound by RNAPII, while H3K9me1 and H3K27me3 were reorganized to the late side of the SV40 regulatory region apparently to repress late transcription. The organization of acetyl-H3 and H3K9me1 but not H3K27me3 required active transcription.

Androgen Receptor-Mediated Transcription in Prostate Cancer

Androgen receptor (AR)-mediated transcription is critical in almost all stages of prostate cancer (PCa) growth and differentiation. This process involves a complex interplay of coregulatory proteins, chromatin remodeling complexes, and other transcription factors that work with AR at cis-regulatory enhancer regions to induce the spatiotemporal transcription of target genes. This enhancer-driven mechanism is remarkably dynamic and undergoes significant alterations during PCa progression. In this review, we discuss the AR mechanism of action in PCa with a focus on how cis-regulatory elements modulate gene expression. We explore emerging evidence of genetic variants that can impact AR regulatory regions and alter gene transcription in PCa. Finally, we highlight several outstanding questions and discuss potential mechanisms of this critical transcription factor.

EnhancerAtlas 2.0: an updated resource with enhancer annotation in 586 tissue/cell types across nine species

Enhancers are distal cis-regulatory elements that activate the transcription of their target genes. They regulate a wide range of important biological functions and processes, including embryogenesis, development, and homeostasis. As more and more large-scale technologies were developed for enhancer identification, a comprehensive database is highly desirable for enhancer annotation based on various genome-wide profiling datasets across different species. Here, we present an updated database EnhancerAtlas 2.0 (http://www.enhanceratlas.org/indexv2.php), covering 586 tissue/cell types that include a large number of normal tissues, cancer cell lines, and cells at different development stages across nine species. Overall, the database contains 13 494 603 enhancers, which were obtained from 16 055 datasets using 12 high-throughput experiment methods (e.g. H3K4me1/H3K27ac, DNase-seq/ATAC-seq, P300, POLR2A, CAGE, ChIA-PET, GRO-seq, STARR-seq and MPRA). The updated version is a huge expansion of the first version, which only contains the enhancers in human cells. In addition, we predicted enhancer–target gene relationships in human, mouse and fly. Finally, the users can search enhancers and enhancer–target gene relationships through five user-friendly, interactive modules. We believe the new annotation of enhancers in EnhancerAtlas 2.0 will facilitate users to perform useful functional analysis of enhancers in various genomes.

Differential expression networks and inheritance patterns of long non-coding RNAs in castor bean seeds

Long non-coding RNAs (lncRNAs) serve as versatile regulators of plant growth and development. The potential functions and inheritance patterns of lncRNAs, as well as the epigenetic regulation of lncRNA itself, remain largely uncharacterized in plant seeds, especially in the persistent endosperm of the dicotyledons. In this study, we investigated diverse RNA-seq data and catalogued 5356 lncRNAs in castor bean seeds. A small fraction of lncRNAs were transcribed from the same direction as the promoters of protein-coding genes (PCgenes) and exhibited strongly coordinated expression with the nearby PCgene. Co-expression analysis with weighted gene co-expression network analysis (WGCNA) showed these lncRNAs to be involved in differential transcription networks between the embryo and endosperm in the early developing seed. Genomic DNA methylation analyses revealed that the expression level of lncRNAs was tightly linked to DNA methylation and that endosperm hypomethylation could promote the expression of linked lncRNAs. Intriguingly, upon hybridization, most lncRNAs with divergent genome sequences between two parents could be reconciled and were expressed according to their parental genome contribution; however, some deviation in the expression of allelic lncRNAs was observed and found to be partially dependent on parental effects. In triploid endosperm, the expression of most lncRNAs was not dosage sensitive, as only 20 lncRNAs had balanced dosage. Our findings not only demonstrate that lncRNAs play potential roles in regulating the development of castor bean endosperm and embryo, but also provide novel insights into the parental effects, allelic expression and epigenetic regulation of lncRNAs in dicotyledonous seeds.

Genome-wide methylome and chromatin interactome identify abnormal enhancer to be risk factor of breast cancer

Enhancer is critical cis regulatory elements in gene expression. To understand whether and how the aberrant enhancer activation may contribute to cancer risk, the differentially methylated enhancers (eDMRs) in normal and malignant breast tissues were identified and analyzed. By incorporating genome-wide chromatin interaction, integrated analysis of eDMRs and target gene expression identified 1,272 enhancer-promoter pairs. Surprisingly, two functionally distinct groups of genes were identified in these pairs, one showing better correlation to enhancer methylation (eRGs) and the other showing better correlation to promoter methylation (pRGs), and the former group is functionally enriched with cancer related genes. Moreover, enhancer methylation based clustering of breast cancer samples is capable of discriminating basal breast cancer from other subtypes. By correlating enhancer methylation status to patient survival, 345 enhancers show the impact on the disease outcome and the majority of their target genes are important regulators of cell survival pathways including known cancer related genes. Together, these results suggest reactivation of enhancers in cancer cells has the add-on effect and contributes to cancer risk in combination.

Mapping the Phosphorylation Pattern of Drosophila melanogaster RNA Polymerase II Carboxyl-Terminal Domain Using Ultraviolet Photodissociation Mass Spectrometry

Phosphorylation of the C-terminal domain of RNA polymerase II (CTD) plays an essential role in eukaryotic transcription by recruiting transcriptional regulatory factors to the active polymerase. However, the scarcity of basic residues and repetitive nature of the CTD sequence impose a huge challenge for site-specific characterization of phosphorylation, hindering our understanding of this crucial biological process. Herein, we apply LC-UVPD-MS methods to analyze post-translational modification along native sequence CTDs. Application of our method to the Drosophila melanogaster CTD reveals the phosphorylation pattern of this model organism for the first time. The divergent nature of fly CTD allows us to derive rules defining how flanking residues affect phosphorylation choice by CTD kinases. Our data support the use of LC-UVPD-MS to decipher the CTD code and determine rules that program its function.

Long non-coding RNAs: Mechanism of action and functional utility

Recent RNA sequencing studies have revealed that most of the human genome is transcribed, but very little of the total transcriptomes has the ability to encode proteins. Long non-coding RNAs (lncRNAs) are non-coding transcripts longer than 200 nucleotides. Members of the non-coding genome include microRNA (miRNA), small regulatory RNAs and other short RNAs. Most of long non-coding RNA (lncRNAs) are poorly annotated. Recent recognition about lncRNAs highlights their effects in many biological and pathological processes. LncRNAs are dysfunctional in a variety of human diseases varying from cancerous to non-cancerous diseases. Characterization of these lncRNA genes and their modes of action may allow their use for diagnosis, monitoring of progression and targeted therapies in various diseases. In this review, we summarize the functional perspectives as well as the mechanism of action of lncRNAs.

The Hierarchy of Transcriptional Activation: From Enhancer to Promoter

Regulatory elements (enhancers) that are remote from promoters play a critical role in the spatial, temporal, and physiological control of gene expression. Studies on specific loci, together with genome-wide approaches, suggest that there may be many common mechanisms involved in enhancer-promoter communication. Here, we discuss the multiprotein complexes that are recruited to enhancers and the hierarchy of events taking place between regulatory elements and promoters.

Dephosphorylating eukaryotic RNA polymerase II

The phosphorylation state of the C-terminal domain of RNA polymerase II is required for the temporal and spatial recruitment of various factors that mediate transcription and RNA processing throughout the transcriptional cycle. Therefore, changes in CTD phosphorylation by site-specific kinases/phosphatases are critical for the accurate transmission of information during transcription. Unlike kinases, CTD phosphatases have been traditionally neglected as they are thought to act as passive negative regulators that remove all phosphate marks at the conclusion of transcription. This over-simplified view has been disputed in recent years and new data assert the active and regulatory role phosphatases play in transcription. We now know that CTD phosphatases ensure the proper transition between different stages of transcription, balance the distribution of phosphorylation for accurate termination and re-initiation, and prevent inappropriate expression of certain genes. In this review, we focus on the specific roles of CTD phosphatases in regulating transcription. In particular, we emphasize how specificity and timing of dephosphorylation are achieved for these phosphatases and consider the various regulatory factors that affect these dynamics.

A class of genes in the HER2 regulon that is poised for transcription in breast cancer cell lines and expressed in human breast tumors

HER2-positive breast cancer accounts for 25% of all cases and has a poor prognosis. Although progress has been made in understanding signal transduction, little is known of how HER2 achieves gene regulation. We performed whole genome expression analysis on a HER2⁺ and HER2⁻ breast cancer cell lines and compared these results to expression in 812 primary tumors stratified by their HER2 expression level. Chip-on-chip with anti-RNA polymerase II was compared among breast cancer cell lines to identify genes that are potentially activated by HER2. The expression levels of these HER2-dependent POL II binding genes were determined for the 812 HER2+/− breast cancer tissues. Genes differentially expressed between HER2+/− cell lines were generally regulated in the same direction as in breast cancer tissues. We identified genes that had POLII binding in HER2⁺ cell lines, but without significant gene expression. Of 737 such genes “poised” for expression in cell lines, 113 genes were significantly differentially expressed in breast tumors in a HER2-dependent manner. Pathway analysis of these 113 genes revealed that a large group of genes were associated with stem cell and progenitor cell control as indicated by networks centered on NANOG, SOX2, OCT3/4. HER2 directs POL II binding to a large number of genes in breast cancer cells. A “poised” class of genes in HER2⁺ cell lines with POLII binding and low RNA expression but is differentially expressed in primary tumors, strongly suggests a role of the microenvironment and further suggests a role for stem cells proliferation in HER2-regulated breast cancer tissue.

Unravelling the RNA-Binding Properties of SAFB Proteins in Breast Cancer Cells

Scaffold attachment factor B1 (SAFB1) and SAFB2 proteins are oestrogen (ER) corepressors that bind to and modulate ER activity through chromatin remodelling or interaction with the basal transcription machinery. SAFB proteins also have an internal RNA-recognition motif but little is known about the RNA-binding properties of SAFB1 or SAFB2. We utilised crosslinking and immunoprecipitation (iCLIP) coupled with high-throughput sequencing to enable a transcriptome-wide mapping of SAFB1 protein-RNA interactions in breast cancer MCF-7 cells. Analysis of crosslinking frequency mapped to transcript regions revealed that SAFB1 binds to coding and noncoding RNAs (ncRNAs). The highest proportion of SAFB1 crosslink sites mapped to ncRNAs, followed by intergenic regions, open reading frames (ORFs), introns, and 3′ or 5′ untranslated regions (UTR). Furthermore, we reveal that SAFB1 binds directly to RNA and its binding is particularly enriched at purine-rich sequences not dissimilar to the RNA-binding motifs for SR proteins. Using RNAi, we also show, for the first time, that single depletion of either SAFB1 or SAFB2 leads to an increase in expression of the other SAFB protein in both MCF-7 and MDA-MD231 breast cancer cells.

ATF4 licenses C/EBPβ activity in human mesenchymal stem cells primed for adipogenesis

A well-established cascade of transcription factor (TF) activity orchestrates adipogenesis in response to chemical cues, yet how cell-intrinsic determinants of differentiation such as cell shape and/or seeding density inform this transcriptional program remain enigmatic. Here, we uncover a novel mechanism licensing transcription in human mesenchymal stem cells (hMSCs) adipogenically primed by confluence. Prior to adipogenesis, confluency promotes heterodimer recruitment of the bZip TFs C/EBPβ and ATF4 to a non-canonical C/EBP DNA sequence. ATF4 depletion decreases both cell-density-dependent transcription and adipocyte differentiation. Global profiling in hMSCs and a novel cell-free assay reveals that ATF4 requires C/EBPβ for genomic binding at a motif distinct from that bound by the C/EBPβ homodimer. Our observations demonstrate that C/EBPβ bridges the transcriptional programs in naïve, confluent cells and early differentiating pre-adipocytes. Moreover, they suggest that homo- and heterodimer formation poise C/EBPβ to execute diverse and stage-specific transcriptional programs by exploiting an expanded motif repertoire.

Kisspeptin Induces Dynamic Chromatin Modifications to Control GnRH Gene Expression

In vitro studies have demonstrated an increase in GnRH gene expression associated with an elevated secretory response to kisspeptin administration, suggesting that kisspeptin mediates GnRH expression at both the secretory and pretranslational levels. However, the kisspeptin-mediated intracellular mechanisms associated with the dynamic chromatin modifications modulating GnRH gene expression are unclear. The studies in this manuscript describe specific histone modifications on the enhancer and promoter of the mouse GnRH (mGnRH) gene induced by kisspeptin in GnRH neuronal cell lines (GT1-7 cells). ChIP assays followed by quantitative real-time PCR (qPCR) demonstrate that 15 and 45 min of 10(-9) M kisspeptin significantly increased histone 3 acetylation (H3Ac) at the kisspeptin response element (KsRE) contained between -3446 and -2806 bp of the mGnRH enhancer (GnRHen) in GT1-7 cells, while no changes were observed in the downstream neuron-specific element (NSE). Moreover, kisspeptin specifically induced acetylation of H3AcK14 and K27 and trimethylation of H3 lysine 4 at the KsRE (markers of active chromatin) and no changes in dimethylation of H3K9 (a marker associated with gene repression). Occupancy of RNA Pol II (RNAPII) and a differential carboxyl-terminal domain (CTD) phosphorylation pattern was observed. An interaction between the NSE and the KsRE via a chromatin loop in the mGnRH gene by kisspeptin was detected by the chromosome conformation capture assay (3C). In conclusion, these results demonstrate that kisspeptin induces histone acetylation/methylation and consequently enhances the formation of a chromatin loop in the mGnRH gene which results in known increase in kisspeptin-dependent mGnRH expression.

Genomic redistribution of GR monomers and dimers mediates transcriptional response to exogenous glucocorticoid in vivo

Glucocorticoids (GCs) are commonly prescribed drugs, but their anti-inflammatory benefits are mitigated by metabolic side effects. Their transcriptional effects, including tissue-specific gene activation and repression, are mediated by the glucocorticoid receptor (GR), which is known to bind as a homodimer to a palindromic DNA sequence. Using ChIP-exo in mouse liver under endogenous corticosterone exposure, we report here that monomeric GR interaction with a half-site motif is more prevalent than homodimer binding. Monomers colocalize with lineage-determining transcription factors in both liver and primary macrophages, and the GR half-site motif drives transcription, suggesting that monomeric binding is fundamental to GR's tissue-specific functions. In response to exogenous GC in vivo, GR dimers assemble on chromatin near ligand-activated genes, concomitant with monomer evacuation of sites near repressed genes. Thus, pharmacological GCs mediate gene expression by favoring GR homodimer occupancy at classic palindromic sites at the expense of monomeric binding. The findings have important implications for improving therapies that target GR.

The distal upstream promoter in Ly49 genes, Pro1, is active in mature NK cells and T cells, does not require TATA boxes, and displays enhancer activity

Missing self recognition of MHC class I molecules is mediated in murine species primarily through the stochastic expression of CD94/NKG2 and Ly49 receptors on NK cells. Previous studies have suggested that the stochastic expression of Ly49 receptors is achieved through the use of an alternate upstream promoter, designated Pro1, that is active only in immature NK cells and operates via the mutually exclusive binding of transcription initiation complexes to closely opposed forward and reverse TATA boxes, with forward transcription being transiently required to activate the downstream promoters, Pro2/Pro3, that are subsequently responsible for transcription in mature NK cells. In this study, we report that Pro1 transcripts are not restricted to immature NK cells but are also found in mature NK cells and T cells, and that Pro1 fragments display strong promoter activity in mature NK cell and T cell lines as well as in immature NK cells. However, the strength of promoter activity in vitro does not correlate well with Ly49 expression in vivo and forward promoter activity is generally weak or undetectable, suggesting that components outside of Pro1 are required for efficient forward transcription. Indeed, conserved sequences immediately upstream and downstream of the core Pro1 region were found to inhibit or enhance promoter activity. Most surprisingly, promoter activity does not require either the forward or reverse TATA boxes, but is instead dependent on residues in the largely invariant central region of Pro1. Importantly, Pro1 displays strong enhancer activity, suggesting that this may be its principal function in vivo.

Histone cross-talk connects protein phosphatase 1α (PP1α) and histone deacetylase (HDAC) pathways to regulate the functional transition of bromodomain-containing 4 (BRD4) for inducible gene expression

Transcription elongation has been recognized as a rate-limiting step for the expression of signal-inducible genes. Through recruitment of positive transcription elongation factor P-TEFb, the bromodomain-containing protein BRD4 plays critical roles in regulating the transcription elongation of a vast array of inducible genes that are important for multiple cellular processes. The diverse biological roles of BRD4 have been proposed to rely on its functional transition between chromatin targeting and transcription regulation. The signaling pathways and the molecular mechanism for regulating this transition process, however, are largely unknown. Here, we report a novel role of phosphorylated Ser(10) of histone H3 (H3S10ph) in governing the functional transition of BRD4. We identified that the acetylated lysines 5 and 8 of nucleosomal histone H4 (H4K5ac/K8ac) is the BRD4 binding site, and the protein phosphatase PP1α and class I histone deacetylase (HDAC1/2/3) signaling pathways are essential for the stress-induced BRD4 release from chromatin. In the unstressed state, phosphorylated H3S10 prevents the deacetylation of nucleosomal H4K5ac/K8ac by HDAC1/2/3, thereby locking up the majority of BRD4 onto chromatin. Upon stress, PP1α-mediated dephosphorylation of H3S10ph allows the deacetylation of nucleosomal H4K5ac/K8ac by HDAC1/2/3, thereby leading to the release of chromatin-bound BRD4 for subsequent recruitment of P-TEFb to enhance the expression of inducible genes. Therefore, our study revealed a novel mechanism that the histone cross-talk between H3S10ph and H4K5ac/K8ac connects PP1α and HDACs to govern the functional transition of BRD4. Combined with previous studies on the regulation of P-TEFb activation, the intricate signaling network for the tight control of transcription elongation is established.

miRNAs and lncRNAs in vascular injury and remodeling

Vascular injury, remodeling, as well as angiogenesis, are the leading causes of coronary or cerebrovascular disease. The blood vessel functional imbalance trends to induce atherosclerosis, hypertension, and pulmonary arterial hypertension. As several genes have been identified to be dynamically regulated during vascular injury and remodeling, it is becoming widely accepted that several types of non-coding RNA, such as microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), are involved in regulating the endothelial cell and vascular smooth muscle cell (VSMC) behaviors. Here, we review the progress of the extant studies on mechanistic, clinical and diagnostic implications of miRNAs and lncRNAs in vascular injury and remodeling, as well as angiogenesis, emphasizing the important roles of miRNAs and lncRNAs in vascular diseases. Furthermore, we introduce the interaction between miRNAs and lncRNAs, and highlight the mechanism through which lncRNAs are regulating the miRNA function. We envisage that continuous in-depth research of non-coding RNAs in vascular disease will have significant implications for the treatment of coronary or cerebrovascular diseases.

A conserved noncoding sequence can function as a spermatocyte-specific enhancer and a bidirectional promoter for a ubiquitously expressed gene and a testis-specific long noncoding RNA

Tissue-specific gene expression is tightly regulated by various elements such as promoters, enhancers, and long noncoding RNAs (lncRNAs). In the present study, we identified a conserved noncoding sequence (CNS1) as a novel enhancer for the spermatocyte-specific mouse testicular cell adhesion molecule 1 (Tcam1) gene. CNS1 was located 3.4kb upstream of the Tcam1 gene and associated with histone H3K4 mono-methylation in testicular germ cells. By the in vitro reporter gene assay, CNS1 could enhance Tcam1 promoter activity only in GC-2spd(ts) cells, which were derived from mouse spermatocytes. When we integrated the 6.9-kb 5'-flanking sequence of Tcam1 with or without a deletion of CNS1 linked to the enhanced green fluorescent protein gene into the chromatin of GC-2spd(ts) cells, CNS1 significantly enhanced Tcam1 promoter activity. These results indicate that CNS1 could function as a spermatocyte-specific enhancer. Interestingly, CNS1 also showed high bidirectional promoter activity in the reporter assay, and consistent with this, the Smarcd2 gene and lncRNA, designated lncRNA-Tcam1, were transcribed from adjacent regions of CNS1. While Smarcd2 was ubiquitously expressed, lncRNA-Tcam1 expression was restricted to testicular germ cells, although this lncRNA did not participate in Tcam1 activation. Ubiquitous Smarcd2 expression was correlated to CpG hypo-methylation of CNS1 and partially controlled by Sp1. However, for lncRNA-Tcam1 transcription, the strong association with histone acetylation and histone H3K4 tri-methylation also appeared to be required. The present data suggest that CNS1 is a spermatocyte-specific enhancer for the Tcam1 gene and a bidirectional promoter of Smarcd2 and lncRNA-Tcam1.

Divergent transcription is associated with promoters of transcriptional regulators

Background

Divergent transcription is a wide-spread phenomenon in mammals. For instance, short bidirectional transcripts are a hallmark of active promoters, while longer transcripts can be detected antisense from active genes in conditions where the RNA degradation machinery is inhibited. Moreover, many described long non-coding RNAs (lncRNAs) are transcribed antisense from coding gene promoters. However, the general significance of divergent lncRNA/mRNA gene pair transcription is still poorly understood. Here, we used strand-specific RNA-seq with high sequencing depth to thoroughly identify antisense transcripts from coding gene promoters in primary mouse tissues.

Results

We found that a substantial fraction of coding-gene promoters sustain divergent transcription of long non-coding RNA (lncRNA)/mRNA gene pairs. Strikingly, upstream antisense transcription is significantly associated with genes related to transcriptional regulation and development. Their promoters share several characteristics with those of transcriptional developmental genes, including very large CpG islands, high degree of conservation and epigenetic regulation in ES cells. In-depth analysis revealed a unique GC skew profile at these promoter regions, while the associated coding genes were found to have large first exons, two genomic features that might enforce bidirectional transcription. Finally, genes associated with antisense transcription harbor specific H3K79me2 epigenetic marking and RNA polymerase II enrichment profiles linked to an intensified rate of early transcriptional elongation.

Conclusions

We concluded that promoters of a class of transcription regulators are characterized by a specialized transcriptional control mechanism, which is directly coupled to relaxed bidirectional transcription.

RNA in unexpected places: long non-coding RNA functions in diverse cellular contexts

The increased application of transcriptome-wide profiling approaches has led to an explosion in the number of documented long non-coding RNAs (lncRNAs). While these new and enigmatic players in the complex transcriptional milieu are encoded by a significant proportion of the genome, their functions are mostly unknown. Early discoveries support a paradigm in which lncRNAs regulate transcription via chromatin modulation, but new functions are steadily emerging. Given the biochemical versatility of RNA, lncRNAs may be used for various tasks, including post-transcriptional regulation, organization of protein complexes, cell-cell signalling and allosteric regulation of proteins.

Long non-coding RNA identification over mouse brain development by integrative modeling of chromatin and genomic features

In silico prediction of genomic long non-coding RNAs (lncRNAs) is prerequisite to the construction and elucidation of non-coding regulatory network. Chromatin modifications marked by chromatin regulators are important epigenetic features, which can be captured by prevailing high-throughput approaches such as ChIP sequencing. We demonstrate that the accuracy of lncRNA predictions can be greatly improved when incorporating high-throughput chromatin modifications over mouse embryonic stem differentiation toward adult Cerebellum by logistic regression with LASSO regularization. The discriminating features include H3K9me3, H3K27ac, H3K4me1, open reading frames and several repeat elements. Importantly, chromatin information is suggested to be complementary to genomic sequence information, highlighting the importance of an integrated model. Applying integrated model, we obtain a list of putative lncRNAs based on uncharacterized fragments from transcriptome assembly. We demonstrate that the putative lncRNAs have regulatory roles in vicinity of known gene loci by expression and Gene Ontology enrichment analysis. We also show that the lncRNA expression specificity can be efficiently modeled by the chromatin data with same developmental stage. The study not only supports the biological hypothesis that chromatin can regulate expression of tissue-specific or developmental stage-specific lncRNAs but also reveals the discriminating features between lncRNA and coding genes, which would guide further lncRNA identifications and characterizations.

A novel method to predict regulatory regions based on histone mark landscapes in macrophages

Macrophages as phagocytes and professional antigen presenting cells play critical roles in both innate and adaptive immunity. Main transcription factors acting during their differentiation and function are known, but the behavior and co-operation of these factors still remained unexplored. We introduce a new approach to map nucleosome-free regions using exclusively active enhancer and core promoter marking histone modification ChIP-seq data. We could detect approximately 56,000 potential active enhancers/promoters showing different lengths and histone modification shapes. Beside the highly enriched PU.1 and C/EBP sites, we could also detect binding sites for RUNX and AP-1, as well as for the MiT (MITF-TFE) family and MEF2 proteins. The PU.1 and C/EBP transcription factors are known for transforming cells into macrophages. The other transcription factors found in this study can play a role in macrophages as well, since it is known that the MiT family proteins are responsible for phagocytic activity and the MEF2 proteins specify monocytic differentiation over the granulocyte direction. Our results imply that this method can provide novel information about transcription factor organization at enhancers and core promoters as well as about the histone modifications surrounding regulatory regions in any immune or other cell types.

GBSA: a comprehensive software for analysing whole genome bisulfite sequencing data

High-throughput sequencing is increasingly being used in combination with bisulfite (BS) assays to study DNA methylation at nucleotide resolution. Although several programmes provide genome-wide alignment of BS-treated reads, the resulting information is not readily interpretable and often requires further bioinformatic steps for meaningful analysis. Current post-alignment BS-sequencing programmes are generally focused on the gene-specific level, a restrictive feature when analysis in the non-coding regions, such as enhancers and intergenic microRNAs, is required. Here, we present Genome Bisulfite Sequencing Analyser (GBSA—http://ctrad-csi.nus.edu.sg/gbsa), a free open-source software capable of analysing whole-genome bisulfite sequencing data with either a gene-centric or gene-independent focus. Through analysis of the largest published data sets to date, we demonstrate GBSA’s features in providing sequencing quality assessment, methylation scoring, functional data management and visualization of genomic methylation at nucleotide resolution. Additionally, we show that GBSA’s output can be easily integrated with other high-throughput sequencing data, such as RNA-Seq or ChIP-seq, to elucidate the role of methylated intergenic regions in gene regulation. In essence, GBSA allows an investigator to explore not only known loci but also all the genomic regions, for which methylation studies could lead to the discovery of new regulatory mechanisms.

A multi-parametric flow cytometric assay to analyze DNA-protein interactions

Interactions between DNA and transcription factors (TFs) guide cellular function and development, yet the complexities of gene regulation are still far from being understood. Such understanding is limited by a paucity of techniques with which to probe DNA–protein interactions. We have devised magnetic protein immobilization on enhancer DNA (MagPIE), a simple, rapid, multi-parametric assay using flow cytometric immunofluorescence to reveal interactions among TFs, chromatin structure and DNA. In MagPIE, synthesized DNA is bound to magnetic beads, which are then incubated with nuclear lysate, permitting sequence-specific binding by TFs, histones and methylation by native lysate factors that can be optionally inhibited with small molecules. Lysate protein–DNA binding is monitored by flow cytometric immunofluorescence, which allows for accurate comparative measurement of TF-DNA affinity. Combinatorial fluorescent staining allows simultaneous analysis of sequence-specific TF-DNA interaction and chromatin modification. MagPIE provides a simple and robust method to analyze complex epigenetic interactions in vitro.

Noncoding transcription at enhancers: general principles and functional models

Mammalian genomes are extensively transcribed outside the borders of protein-coding genes. Genome-wide studies recently demonstrated that cis-regulatory genomic elements implicated in transcriptional control, such as enhancers and locus-control regions, represent major sites of extragenic noncoding transcription. Enhancer-templated transcripts provide a quantitatively small contribution to the total amount of cellular nonribosomal RNA; nevertheless, the possibility that enhancer transcription and the resulting enhancer RNAs may, in some cases, have functional roles, rather than represent mere transcriptional noise at accessible genomic regions, is supported by an increasing amount of experimental data. In this article we review the current knowledge on enhancer transcription and its functional implications.

Research resource: whole-genome estrogen receptor α binding in mouse uterine tissue revealed by ChIP-seq

To advance understanding of mechanisms leading to biological and transcriptional endpoints related to estrogen action in the mouse uterus, we have mapped ERα and RNA polymerase II (PolII) binding sites using chromatin immunoprecipitation followed by sequencing of enriched chromatin fragments. In the absence of hormone, 5184 ERα-binding sites were apparent in the vehicle-treated ovariectomized uterine chromatin, whereas 17,240 were seen 1 h after estradiol (E₂) treatment, indicating that some sites are occupied by unliganded ERα, and that ERα binding is increased by E₂. Approximately 15% of the uterine ERα-binding sites were adjacent to (<10 kb) annotated transcription start sites, and many sites are found within genes or are found more than 100 kb distal from mapped genes; however, the density (sites per base pair) of ERα-binding sites is significantly greater adjacent to promoters. An increase in quantity of sites but no significant positional differences were seen between vehicle and E₂-treated samples in the overall locations of ERα-binding sites either distal from, adjacent to, or within genes. Analysis of the PolII data revealed the presence of poised promoter-proximal PolII on some highly up-regulated genes. Additionally, corecruitment of PolII and ERα to some distal enhancer regions was observed. A de novo motif analysis of sequences in the ERα-bound chromatin confirmed that estrogen response elements were significantly enriched. Interestingly, in areas of ERα binding without predicted estrogen response element motifs, homeodomain transcription factor-binding motifs were significantly enriched. The integration of the ERα- and PolII-binding sites from our uterine sequencing of enriched chromatin fragments data with transcriptional responses revealed in our uterine microarrays has the potential to greatly enhance our understanding of mechanisms governing estrogen response in uterine and other estrogen target tissues.

An RNA-independent linkage of noncoding transcription to long-range enhancer function

The detection of noncoding transcription at multiple enhancers within the mammalian genome raises critical questions regarding whether and how this activity contributes to enhancer function. Here, using in vivo analysis of a human growth hormone (hGH) transgene locus, we report that activation of a domain of noncoding transcription adjacent to the long-range hGH-N enhancer, HSI, is established by the enhancer independent of any interactions with its target promoter. We further demonstrate that the appearance of this enhancer-linked noncoding transcription is temporally and spatially concordant with induction of hGH-N in the embryonic pituitary. Finally, we show that the level of transcriptional enhancement of hGH-N by HSI is directly related to the intensity of HSI-dependent noncoding transcription and is fully independent of the structure of the locally transcribed RNA. These data extend our understanding of the relationship of long-range enhancer activity to enhancer-dependent noncoding transcription and establish a model that may be of general relevance to additional mammalian loci.

Chromatin signatures of active enhancers

Gene-distal cis-regulatory sequences, such as enhancers, are key contributors of tissue-specific gene expression. In particular, enhancers can be located up to hundreds of kilobases from the promoters that they control, making their identification challenging. Thanks to the recent technological advances to map histone modifications and chromatin-associated factors genome-wide, several studies have begun to characterize chromatin signatures of active enhancers. Here, we discuss some of these results and how they provide new insights into the tissue-specific organization of enhancer repertoires.

Cell-type specific and combinatorial usage of diverse transcription factors revealed by genome-wide binding studies in multiple human cells

Cell-type diversity is governed in part by differential gene expression programs mediated by transcription factor (TF) binding. However, there are few systematic studies of the genomic binding of different types of TFs across a wide range of human cell types, especially in relation to gene expression. In the ENCODE Project, we have identified the genomic binding locations across 11 different human cell types of CTCF, RNA Pol II (RNAPII), and MYC, three TFs with diverse roles. Our data and analysis revealed how these factors bind in relation to genomic features and shape gene expression and cell-type specificity. CTCF bound predominantly in intergenic regions while RNAPII and MYC preferentially bound to core promoter regions. CTCF sites were relatively invariant across diverse cell types, while MYC showed the greatest cell-type specificity. MYC and RNAPII co-localized at many of their binding sites and putative target genes. Cell-type specific binding sites, in particular for MYC and RNAPII, were associated with cell-type specific functions. Patterns of binding in relation to gene features were generally conserved across different cell types. RNAPII occupancy was higher over exons than adjacent introns, likely reflecting a link between transcriptional elongation and splicing. TF binding was positively correlated with the expression levels of their putative target genes, but combinatorial binding, in particular of MYC and RNAPII, was even more strongly associated with higher gene expression. These data illuminate how combinatorial binding of transcription factors in diverse cell types is associated with gene expression and cell-type specific biology.

Initiating RNA polymerase II and TIPs as hallmarks of enhancer activity and tissue-specificity

In past years, many efforts were invested to define epigenetic features associated with enhancers of transcription. We propose that both transcription initiation and the H3K4me3 histone modification are among the best hallmarks of active enhancers in several primary tissues and extend the concept of large transcription initiation platforms (TIPs).

H3K4 tri-methylation provides an epigenetic signature of active enhancers

Combinations of post-translational histone modifications shape the chromatin landscape during cell development in eukaryotes. However, little is known about the modifications exactly delineating functionally engaged regulatory elements. For example, although histone H3 lysine 4 mono-methylation (H3K4me1) indicates the presence of transcriptional gene enhancers, it does not provide clearcut information about their actual position and stage-specific activity. Histone marks were, therefore, studied here at genomic loci differentially expressed in early stages of T-lymphocyte development. The concomitant presence of the three H3K4 methylation states (H3K4me1/2/3) was found to clearly reflect the activity of bona fide T-cell gene enhancers. Globally, gain or loss of H3K4me2/3 at distal genomic regions correlated with, respectively, the induction or the repression of associated genes during T-cell development. In the Tcrb gene enhancer, the H3K4me3-to-H3K4me1 ratio decreases with the enhancer's strength. Lastly, enhancer association of RNA-polymerase II (Pol II) correlated with the presence of H3K4me3 and Pol II accumulation resulted in local increase of H3K4me3. Our results suggest the existence of functional links between Pol II occupancy, H3K4me3 enrichment and enhancer activity.

Allele-specific distribution of RNA polymerase II on female X chromosomes

While the distribution of RNA polymerase II (PolII) in a variety of complex genomes is correlated with gene expression, the presence of PolII at a gene does not necessarily indicate active expression. Various patterns of PolII binding have been described genome wide; however, whether or not PolII binds at transcriptionally inactive sites remains uncertain. The two X chromosomes in female cells in mammals present an opportunity to examine each of the two alleles of a given locus in both active and inactive states, depending on which X chromosome is silenced by X chromosome inactivation. Here, we investigated PolII occupancy and expression of the associated genes across the active (Xa) and inactive (Xi) X chromosomes in human female cells to elucidate the relationship of gene expression and PolII binding. We find that, while PolII in the pseudoautosomal region occupies both chromosomes at similar levels, it is significantly biased toward the Xa throughout the rest of the chromosome. The general paucity of PolII on the Xi notwithstanding, detectable (albeit significantly reduced) binding can be observed, especially on the evolutionarily younger short arm of the X. PolII levels at genes that escape inactivation correlate with the levels of their expression; however, additional PolII sites can be found at apparently silenced regions, suggesting the possibility of a subset of genes on the Xi that are poised for expression. Consistent with this hypothesis, we show that a high proportion of genes associated with PolII-accessible sites, while silenced in GM12878, are expressed in other female cell lines.

Transcription initiation platforms and GTF recruitment at tissue-specific enhancers and promoters

Recent work has shown that RNA polymerase (Pol) II can be recruited to and transcribe distal regulatory regions. Here we analyzed transcription initiation and elongation through genome-wide localization of Pol II, general transcription factors (GTFs) and active chromatin in developing T cells. We show that Pol II and GTFs are recruited to known T cell-specific enhancers. We extend this observation to many new putative enhancers, a majority of which can be transcribed with or without polyadenylation. Importantly, we also identify genomic features called transcriptional initiation platforms (TIPs) that are characterized by large areas of Pol II and GTF recruitment at promoters, intergenic and intragenic regions. TIPs show variable widths (0.4-10 kb) and correlate with high CpG content and increased tissue specificity at promoters. Finally, we also report differential recruitment of TFIID and other GTFs at promoters and enhancers. Overall, we propose that TIPs represent important new regulatory hallmarks of the genome.

Epigenetic Control of MicroRNA Expression and Aging

MicroRNAs are a major category among the noncoding RNA fraction that negatively regulate gene expression at the post-transcriptional level, by either degrading the target messages or inhibiting their translation. MicroRNAs may be referred to as ‘dimmer switches’ of gene expression, because of their ability to repress gene expression without completely silencing it. Whether through up-regulating specific groups of microRNAs to suppress unwanted gene expressions, or by down-regulating other microRNAs whose target genes’ expression is necessary for cellular function, such as cell proliferation, apoptosis, or differentiation, these regulatory RNAs play pivotal roles in a wide variety of cellular processes. The equilibrium between these two groups of microRNA expressions largely determines the function of particular cell types. Our recent results with several model systems show that upon aging, there is a trend of up-regulation of microRNA expression, with concomitant inverse down-regulation of target genes. This review addresses molecular mechanisms that may provide the underlying control for this up-regulating trend, focusing on activation by various microRNAs’ own promoters, through binding with pivotal transcription factors, stress response, methylation of clustered DNA domains, etc. Thus, epigenomic control of aging may be due in part to heightened promoter activation of unwanted microRNA expressions, which in turn down-regulate their target gene products. Overriding and dampening the activation of these noncoding RNAs may prove to be a new frontier for future research, to delay aging and extend healthy life-span.

Pervasive transcription - Lessons from yeast

Pervasive transcription is now accepted to be a general feature of eukaryotic genomes, generating short and long non-coding RNAs (ncRNAs). Growing number of examples have shown that regulatory ncRNAs can control gene expression and chromatin domain formation. In this review, we discuss recent reports that show that Saccharomyces cerevisiae's genome also supports pervasive transcription, which is strongly controlled by RNA decay pathways and nucleosome positioning. We therefore propose that S. cerevisiae is an excellent model for studying large ncRNAs, which has already provided important examples of antisense-mediated transcriptional silencing.

Processing ChIP-chip data: from the scanner to the browser

High-density tiling microarrays are increasingly used in combination with chromatin immunoprecipitation (ChIP) assays to delineate the regulation of gene expression. Besides the technical challenges inherent to such complex biological assays, a critical, often daunting issue is the correct interpretation of the sheer amount of raw data generated by utilizing computational methods. Here, we go through the main steps of this intricate process, including optimized chromatin immunoprecipitation on chip (ChIP-chip) data normalization, peak detection, as well as quality control reports. We also describe convenient standalone software suites, including our own, CoCAS, which works on the latest generation of Agilent high-density arrays, allows dye-swap, replicate correlation, and easy connection with genome browsers for results interpretation, or with, e.g., other peak detection algorithms. Overall, the guidelines described herein provide an effective introduction to ChIP-chip technology and analysis.

MED29, a component of the mediator complex, possesses both oncogenic and tumor suppressive characteristics in pancreatic cancer

Mediator complex subunit 29 (MED29) is part of a large multiprotein coactivator complex that mediates regulatory signals from gene-specific activators to general transcription machinery in RNA polymerase II mediated transcription. We previously found that MED29 is amplified and overexpressed in pancreatic cancer and that MED29 silencing leads to decreased cell survival in PANC-1 pancreatic cancer cells with high MED29 expression. Here we further demonstrate decreased migration, invasion and colony formation in PANC-1 cells after MED29 silencing. Unexpectedly, lentiviral-based overexpression of MED29 led to decreased proliferation of NIH/3T3 cells as well as MIAPaCa-2 pancreatic cancer cells with low endogenous expression. More importantly, subcutaneous inoculation of the MED29-transduced pancreatic cancer cells into immuno-compromised mice resulted in dramatic tumor suppression. The mock-control mice developed large tumors, whereas the animals with MED29-xenografts showed both decreased tumor incidence and a major reduction in tumor size. Gene expression analysis in the MED29-transduced pancreatic cancer cells revealed differential expression of genes involved in control of cell cycle and cell division. The observed gene expression changes are expected to modulate the cell cycle in a way that leads to reduced cell growth, explaining the in vivo tumor suppressive phenotype. Taken together, these data implicate MED29 as an important regulator of key cellular functions in pancreatic cancer with both oncogenic and tumor suppressive characteristics. Such a dualistic role appears to be more common than previously thought and is likely to depend on the genetic background of the cancer cells and their surrounding environment.

Functional and mechanistic diversity of distal transcription enhancers

Biological differences among metazoans and between cell types in a given organism arise in large part due to differences in gene expression patterns. Gene-distal enhancers are key contributors to these expression patterns, exhibiting both sequence diversity and cell type specificity. Studies of long-range interactions indicate that enhancers are often important determinants of nuclear organization, contributing to a general model for enhancer function that involves direct enhancer-promoter contact. However, mechanisms for enhancer function are emerging that do not fit solely within such a model, suggesting that enhancers as a class of DNA regulatory element may be functionally and mechanistically diverse.

Dynamic chromatin modifications control GnRH gene expression during neuronal differentiation and protein kinase C signal transduction

GnRH, a neuropeptide produced by rare, specialized hypothalamic secretory neurons, is critical for reproduction. During development, GnRH gene expression increases as neurons migrate from the olfactory placode to the hypothalamus, with highest levels in the mature, postmitotic state. While neuronal differentiation is known to be controlled by chromatin modulations, the role of chromatin dynamics in GnRH gene regulation has not been studied. Here, we use mature and immature GnRH neuronal cell models to show that both neuron-specific and protein kinase C regulation of GnRH expression are mediated by chromatin structure and histone modifications. Only in GT1-7 mature GnRH neuronal cells did GnRH regulatory elements display high sensitivity to DNase and enrichment of active histone markers histone-H3 acetylation and H3 lysine 4 trimethylation (H3K4-Me3), as well as RNA polymerase II (RNAPII) binding and enhancer RNA transcription. In contrast, H3K9-Me2, a marker of inactive chromatin, was highest in nonneuronal cells, low in GT1-7 cells, and intermediate in immature GnRH neuronal cells. The chromatin of the GnRH gene was therefore active in mature GnRH neuronal cells, inactive in nonneuronal cells, but not fully inactive in immature GnRH neuronal cells. Activation of protein kinase C (PKC) potently represses GnRH expression. PKC activation caused closing of the chromatin and decreased RNAPII occupancy at the GnRH minimal promoter (-278/-97). At GnRH-Enhancer-1 (-2404/-2100), PKC activation decreased phosphorylated-RNAPII binding, enhancer RNA transcription, and H3 acetylation, and reciprocally increased H3K9-Me2. Chromatin modifications therefore participate in the dynamic regulation and specification of GnRH expression to differentiated hypothalamic neurons.

The long arm of long noncoding RNAs: roles as sensors regulating gene transcriptional programs

A major surprise arising from genome-wide analyses has been the observation that the majority of the genome is transcribed, generating noncoding RNAs (ncRNAs). It is still an open question whether some or all of these ncRNAs constitute functional networks regulating gene transcriptional programs. However, in light of recent discoveries and given the diversity and flexibility of long ncRNAs and their abilities to nucleate molecular complexes and to form spatially compact arrays of complexes, it becomes likely that many or most ncRNAs act as sensors and integrators of a wide variety of regulated transcriptional responses and probably epigenetic events. Because many RNA-binding proteins, on binding RNAs, show distinct allosteric conformational alterations, we suggest that a ncRNA/RNA-binding protein-based strategy, perhaps in concert with several other mechanistic strategies, serves to integrate transcriptional, as well as RNA processing, regulatory programs.

Nonconventional initiation complex assembly by STAT and NF-kappaB transcription factors regulates nitric oxide synthase expression

Transcriptional regulation of the Nos2 gene encoding inducible nitric oxide synthase (iNOS) requires type I interferon (IFN-I) signaling and additional signals emanating from pattern recognition receptors. Here we showed sequential and cooperative contributions of the transcription factors ISGF3 (a complex containing STAT1, STAT2, and IRF9 subunits) and NF-κB to the transcriptional induction of the Nos2 gene in macrophages infected with the intracellular bacterial pathogen Listeria monocytogenes. NF-κB preceded ISGF3 at the Nos2 promoter and generated a transcriptional memory effect by depositing basal transcription factor TFIIH with the associated CDK7 kinase for serine 5 phosphorylation of the RNA polymerase II (pol II) carboxyterminal domain (CTD). Subsequent to TFIIH deposition by NF-κB, ISGF3 attracted the pol II enzyme and phosphorylation at CTD S5 occurred. Thus, STATs and NF-κB cooperate through pol II promoter recruitment and the phosphorylation of its CTD, respectively, as a prerequisite for productive elongation of iNOS mRNA.

Functional analysis and identification of cis-regulatory elements of human chromosome 21 gene promoters

Given the inherent limitations of in silico studies relying solely on DNA sequence analysis, the functional characterization of mammalian promoters and associated cis-regulatory elements requires experimental support, which demands cloning and analysis of putative promoter regions. Focusing on human chromosome 21, we cloned 182 gene promoters of 2500 bp in length and conducted reporter gene assays on transfected-cell arrays. We found 56 promoters that were active in HEK293 cells, while another 49 promoters could be activated by treatment of cells with Trichostatin A or depletion of serum. We observed high correlations between promoter activities and endogenous transcript levels, RNA polymerase II occupancy, CpG islands and core promoter elements. Truncation of a subset of 62 promoters to ∼500 bp revealed that truncation rarely resulted in loss of activity, but rather in loss of responses to external stimuli, suggesting the presence of cis-regulatory response elements within distal promoter regions. In these regions, we found a strong enrichment of transcription factor binding sites that could potentially activate gene expression in the presence of stimuli. This study illustrates the modular functional architecture of chromosome 21 promoters and helps to reveal the complex mechanisms governing transcriptional regulation.

Widespread transcription at neuronal activity-regulated enhancers

We used genome-wide sequencing methods to study stimulus-dependent enhancer function in neurons. We identified ∼12,000 neuronal activity-regulated enhancers that are bound by the general transcriptional co-activator CBP in an activity-dependent manner. A function of CBP at enhancers may be to recruit RNA polymerase II (RNAPII), as we also observed activity-regulated RNAPII binding to thousands of enhancers. Remarkably, RNAPII at enhancers transcribes bi-directionally a novel class of enhancer RNAs (eRNAs) within enhancer domains defined by the presence of histone H3 that is mono-methylated at lysine 4 (H3K4me1). The level of eRNA expression at neuronal enhancers positively correlates with the level of mRNA synthesis at nearby genes, suggesting that eRNA synthesis occurs specifically at enhancers that are actively engaged in promoting mRNA synthesis. These findings reveal that a widespread mechanism of enhancer activation involves RNAPII binding and eRNA synthesis.

c-Myc is required for the CHREBP-dependent activation of glucose-responsive genes

Glucose regulates programs of gene expression that orchestrate changes in cellular phenotype in several metabolically active tissues. Carbohydrate response element-binding protein (ChREBP) and its binding partner, Mlx, mediate glucose-regulated gene expression by binding to carbohydrate response elements on target genes, such as the prototypical glucose-responsive gene, liver-type pyruvate kinase (Pklr). c-Myc is also required for the glucose response of the Pklr gene, although the relationship between c-Myc and ChREBP has not been defined. Here we describe the molecular events of the glucose-mediated activation of Pklr and determine the effects of decreasing the activity or abundance of c-Myc on this process. Time-course chromatin immunoprecipitation revealed a set of transcription factors [hepatocyte nuclear factor (HNF)1alpha, HNF4alpha, and RNA polymerase II (Pol II)] constitutively resident on the Pklr promoter, with a relative enrichment of acetylated histones 3 and 4 in the same region of the gene. Glucose did not affect HNF1alpha binding or the acetylation of histones H3 or H4. By contrast, glucose promoted the recruitment of ChREBP and c-Myc and increased the occupancy of HNF4alpha and RNA Pol II, which were coincident with the glucose-mediated increase in transcription as determined by a nuclear run-on assay. Depletion of c-Myc activity using a small molecule inhibitor (10058-F4/1RH) abolished the glucose-mediated recruitment of HNF4alpha, ChREBP, and RNA Pol II, without affecting basal gene expression, histone acetylation, and HNF1alpha or basal HNF4alpha occupancy. The activation and recruitment of ChREBP to several glucose-responsive genes were blocked by 1RH, indicating a general necessity for c-Myc in this process.

Life without RNAi: noncoding RNAs and their functions in Saccharomyces cerevisiae

There are a number of well-characterized and fundamental roles for noncoding RNAs (ncRNAs) in gene regulation in all kingdoms of life. ncRNAs, such as ribosomal RNAs, transfer RNAs, small nuclear RNAs, small nucleolar RNAs, and small interfering RNAs, can serve catalytic and scaffolding functions in transcription, messenger RNA processing, translation, and RNA degradation. Recently, our understanding of gene expression has been dramatically challenged by the identification of large and diverse populations of novel ncRNAs in the eukaryotic genomes surveyed thus far. Studies carried out using the budding yeast Saccharomyces cerevisiae indicate that at least some coding genes are regulated by these novel ncRNAs. S. cerevisiae lacks RNA interference (RNAi) and, thus, provides an ideal system for studying the RNAi-independent mechanisms of ncRNA-based gene regulation. The current picture of gene regulation is one of great unknowns, in which the transcriptional environment surrounding a given locus may have as much to do with its regulation as its DNA sequence or local chromatin structure. Drawing on the recent research in S. cerevisiae and other organisms, this review will discuss the identification of ncRNAs, their origins and processing, and several models that incorporate ncRNAs into the regulation of gene expression and chromatin structure.