First Exon Length Controls Active Chromatin Signatures and Transcription

GC-AG Introns Features in Long Non-coding and Protein-Coding Genes Suggest Their Role in Gene Expression Regulation

Long non-coding RNAs (lncRNAs) are recognized as an important class of regulatory molecules involved in a variety of biological functions. However, the regulatory mechanisms of long non-coding genes expression are still poorly understood. The characterization of the genomic features of lncRNAs is crucial to get insight into their function. In this study, we exploited recent annotations by GENCODE to characterize the genomic and splicing features of long non-coding genes in comparison with protein-coding ones, both in human and mouse. Our analysis highlighted differences between the two classes of genes in terms of their gene architecture. Significant differences in the splice sites usage were observed between long non-coding and protein-coding genes (PCG). While the frequency of non-canonical GC-AG splice junctions represents about 0.8% of total splice sites in PCGs, we identified a significant enrichment of the GC-AG splice sites in long non-coding genes, both in human (3.0%) and mouse (1.9%). In addition, we found a positional bias of GC-AG splice sites being enriched in the first intron in both classes of genes. Moreover, a significant shorter length and weaker donor and acceptor sites were found comparing GC-AG introns to GT-AG introns. Genes containing at least one GC-AG intron were found conserved in many species, more prone to alternative splicing and a functional analysis pointed toward their enrichment in specific biological processes such as DNA repair. Our study shows for the first time that GC-AG introns are mainly associated with lncRNAs and are preferentially located in the first intron. Additionally, we discovered their regulatory potential indicating the existence of a new mechanism of non-coding and PCGs expression regulation.

Regulating Divergent Transcriptomes through mRNA Splicing and Its Modulation Using Various Small Compounds

Human transcriptomes are more divergent than genes and contribute to the sophistication of life. This divergence is derived from various isoforms arising from alternative splicing. In addition, alternative splicing regulated by spliceosomal factors and RNA structures, such as the RNA G-quadruplex, is important not only for isoform diversity but also for regulating gene expression. Therefore, abnormal splicing leads to serious diseases such as cancer and neurodegenerative disorders. In the first part of this review, we describe the regulation of divergent transcriptomes using alternative mRNA splicing. In the second part, we present the relationship between the disruption of splicing and diseases. Recently, various compounds with splicing inhibitor activity were established. These splicing inhibitors are recognized as a biological tool to investigate the molecular mechanism of splicing and as a potential therapeutic agent for cancer treatment. Food-derived compounds with similar functions were found and are expected to exhibit anticancer effects. In the final part, we describe the compounds that modulate the messenger RNA (mRNA) splicing process and their availability for basic research and future clinical potential.

Histone lysine demethylase KDM5B maintains chronic myeloid leukemia via multiple epigenetic actions

The histone lysine demethylase KDM5 family is implicated in normal development and stem cell maintenance by epigenetic modulation of histone methylation status. Deregulation of the KDM5 family has been reported in various types of cancers, including hematological malignancies. However, their transcriptional regulatory roles in the context of leukemia remain unclear. Here, we find that KDM5B is strongly expressed in normal CD34⁺ hematopoietic stem/progenitor cells and chronic myeloid leukemia (CML) cells. Knockdown of KDM5B in K562 CML cells reduced leukemia colony-forming potential. Transcriptome profiling of KDM5B knockdown K562 cells revealed the deregulation of genes involved in myeloid differentiation and Toll-like receptor signaling. Through the integration of transcriptome and ChIP-seq profiling data, we show that KDM5B is enriched at the binding sites of the GATA and AP-1 transcription factor families, suggesting their collaborations in the regulation of transcription. Even though the binding of KDM5B substantially overlapped with H3K4me1 or H3K4me3 mark at gene promoters, only a small subset of the KDM5B targets showed differential expression in association with the histone demethylation activity. By characterizing the interacting proteins in K562 cells, we discovered that KDM5B recruits protein complexes involved in the mRNA processing machinery, implying an alternative epigenetic action mediated by KDM5B in gene regulation. Our study highlights the oncogenic functions of KDM5B in CML cells and suggests that KDM5B is vital to the transcriptional regulation via multiple epigenetic mechanisms.

Mechanistic Dissection of RNA-Binding Proteins in Regulated Gene Expression at Chromatin Levels

Eukaryotic genomes are known to prevalently transcribe diverse classes of RNAs, virtually all of which, including nascent RNAs from protein-coding genes, are now recognized to have regulatory functions in gene expression, suggesting that RNAs are both the products and the regulators of gene expression. Their functions must enlist specific RNA-binding proteins (RBPs) to execute their regulatory activities, and recent evidence suggests that nearly all biochemically defined chromatin regions in the human genome, whether defined for gene activation or silencing, have the involvement of specific RBPs. Interestingly, the boundary between RNA- and DNA-binding proteins is also melting, as many DNA-binding proteins traditionally studied in the context of transcription are able to bind RNAs, some of which may simultaneously bind both DNA and RNA to facilitate network interactions in three-dimensional (3D) genome. In this review, we focus on RBPs that function at chromatin levels, with particular emphasis on their mechanisms of action in regulated gene expression, which is intended to facilitate future functional and mechanistic dissection of chromatin-associated RBPs.

Novel Biomarkers for Prediction of Response to Preoperative Systemic Therapies in Gastric Cancer

Preoperative chemo- and radiotherapeutic strategies followed by surgery are currently a standard approach for treating locally advanced gastric and esophagogastric junction cancer in Western countries. However, in a large number of cases, the tumor is extremely resistant to these treatments and the patients are exposed to unnecessary toxicity and delayed surgical therapy. The current clinical trials evaluating the combination of preoperative systemic therapies with modern targeted and immunotherapeutic agents represent a unique opportunity for identifying predictive biomarkers of response to select patients that would benefit the most from these treatments. However, it is of utmost importance that these potential biomarkers are corroborated by extensive preclinical and translational research. The aim of this review article is to present the most promising biomarkers of response to classic chemotherapeutic, anti-HER2, antiangiogenic, and immunotherapeutic agents that can be potentially useful for personalized preoperative systemic therapies in gastric cancer patients.

Nascent RNA analyses: tracking transcription and its regulation

The programmes that direct an organism's development and maintenance are encoded in its genome. Decoding of this information begins with regulated transcription of genomic DNA into RNA. Although transcription and its control can be tracked indirectly by measuring stable RNAs, it is only by directly measuring nascent RNAs that the immediate regulatory changes in response to developmental, environmental, disease and metabolic signals are revealed. Multiple complementary methods have been developed to quantitatively track nascent transcription genome-wide at nucleotide resolution, all of which have contributed novel insights into the mechanisms of gene regulation and transcription-coupled RNA processing. Here we critically evaluate the array of strategies used for investigating nascent transcription and discuss the recent conceptual advances they have provided.

Exon-Mediated Activation of Transcription Starts

The processing of RNA transcripts from mammalian genes occurs in proximity to their transcription. Here, we describe a phenomenon affecting thousands of genes that we call exon-mediated activation of transcription starts (EMATS), in which the splicing of internal exons impacts promoter choice and the expression level of the gene. We observed that evolutionary gain of internal exons is associated with gain of new transcription start sites (TSSs) nearby and increased gene expression. Inhibiting exon splicing reduced transcription from nearby promoters, and creation of new spliced exons activated transcription from cryptic promoters. The strongest effects occurred for weak promoters located proximal and upstream of efficiently spliced exons. Together, our findings support a model in which splicing recruits transcription machinery locally to influence TSS choice and identify exon gain, loss, and regulatory change as major contributors to the evolution of alternative promoters and gene expression in mammals.

The SLC9A-C Mammalian Na+/H+ Exchanger Family: Molecules, Mechanisms, and Physiology

Na⁺/H⁺ exchangers play pivotal roles in the control of cell and tissue pH by mediating the electroneutral exchange of Na⁺ and H⁺ across cellular membranes. They belong to an ancient family of highly evolutionarily conserved proteins, and they play essential physiological roles in all phyla. In this review, we focus on the mammalian Na⁺/H⁺ exchangers (NHEs), the solute carrier (SLC) 9 family. This family of electroneutral transporters constitutes three branches: SLC9A, -B, and -C. Within these, each isoform exhibits distinct tissue expression profiles, regulation, and physiological roles. Some of these transporters are highly studied, with hundreds of original articles, and some are still only rudimentarily understood. In this review, we present and discuss the pioneering original work as well as the current state-of-the-art research on mammalian NHEs. We aim to provide the reader with a comprehensive view of core knowledge and recent insights into each family member, from gene organization over protein structure and regulation to physiological and pathophysiological roles. Particular attention is given to the integrated physiology of NHEs in the main organ systems. We provide several novel analyses and useful overviews, and we pinpoint main remaining enigmas, which we hope will inspire novel research on these highly versatile proteins.

Nascent RNA and the Coordination of Splicing with Transcription

At each active protein-encoding gene, nascent RNA is tethered to the DNA axis by elongating RNA polymerase II (Pol II) and is continuously altered by splicing and other processing events during its synthesis. This review discusses the development of three major methods that enable us to track the conversion of precursor messenger RNA (pre-mRNA) to messenger RNA (mRNA) products in vivo: live-cell imaging, metabolic labeling of RNA, and RNA-seq of purified nascent RNA. These approaches are complementary, addressing distinct issues of transcription rates and intron lifetimes alongside spatial information regarding the gene position of Pol II at which spliceosomes act. The findings will be placed in the context of active transcription units, each of which-because of the presence of nascent RNA, Pol II, and features of the chromatin environment-will recruit a potentially gene-specific constellation of RNA binding proteins and processing machineries.

In vivo measurements reveal a single 5'-intron is sufficient to increase protein expression level in Caenorhabditis elegans

Introns can increase gene expression levels using a variety of mechanisms collectively referred to as Intron Mediated Enhancement (IME). IME has been measured in cell culture and plant models by quantifying expression of intronless and intron-bearing reporter genes in vitro. We developed hardware and software to implement microfluidic chip-based gene expression quantification in vivo. We altered position, number and sequence of introns in reporter genes controlled by the hsp-90 promoter. Consistent with plant and mammalian studies, we determined a single, natural or synthetic, 5'-intron is sufficient for the full IME effect conferred by three synthetic introns, while a 3'-intron is not. We found coding sequence can affect IME; the same three synthetic introns that increase mcherry protein concentration by approximately 50%, increase mEGFP by 80%. We determined IME effect size is not greatly affected by the stronger vit-2 promoter. Our microfluidic imaging approach should facilitate screens for factors affecting IME and other intron-dependent processes.

Development of an MSI-positive colon tumor with aberrant DNA methylation in a PPAP patient

Polymerase proofreading-associated polyposis (PPAP) is a disease caused by germline variations in the POLE and POLD1 genes that encode catalytic subunits of DNA polymerases. Studies of cancer genomes have identified somatic mutations in these genes, suggesting the importance of polymerase proofreading of DNA replication in suppressing tumorigenesis. Here, we identified a germline frameshift variation in the POLE gene (c.4191_4192delCT, p.Tyr1398*) in a case with multiple adenomatous polyps and three synchronous colon cancers. Interestingly, one of the colon cancers showed microsatellite instability-high (MSI-H) and another microsatellite stable. Immunohistochemical staining revealed that the MSI-H tumor cells lost the expression of MLH1 protein. Whole genome sequencing of the MSI-H tumor did not find pathogenic somatic mutations in mismatch repair genes but found frameshift mutations in the TET genes that catalyze 5-methylcytosine hydroxylation. Bisulfite sequencing of the tumor corroborated an increase in the number of hypermethylated regions including the MLH1 promoter. These data indicate that PPAP patients might develop MSI-positive tumors through epigenetic silencing of MLH1. These findings will contribute to comprehensive understanding of the molecular basis of tumors that involve deficiency of proofreading activity of DNA polymerases.

De novo emergence and potential function of human-specific tandem repeats in brain-related loci

Tandem repeats (TRs) are widespread in the genomes of all living organisms. In eukaryotes, they are found in both coding and noncoding regions and have potential roles in the regulation of cellular processes such as transcription, translation and in the modification of protein structure. Recent studies have highlighted TRs as a key regulator of gene expression and a potential contributor to human evolution. Thus, TRs are emerging as an important source of variation that can result in differential gene expression at intra- and inter-species levels. In this study, we performed a genome-wide survey to identify TRs that have emerged in the human lineage. We further examined these loci to explore their potential functional significance for human evolution. We identified 152 human-specific TR (HSTR) loci containing a repeat unit of more than ten bases, with most of them showing a repeat count of two. Gene set enrichment analysis showed that HSTR-associated genes were associated with biological functions in brain development and synapse function. In addition, we compared gene expression of human HSTR loci with orthologues from non-human primates (NHP) in seven different tissues. Strikingly, the expression level of HSTR-associated genes in brain tissues was significantly higher in human than in NHP. These results suggest the possibility that de novo emergence of TRs could have resulted in altered gene expression in humans within a short-time frame and contributed to the rapid evolution of human brain function.

Introns as Gene Regulators: A Brick on the Accelerator

A picture is beginning to emerge from a variety of organisms that for a subset of genes, the most important sequences that regulate expression are situated not in the promoter but rather are located within introns in the first kilobase of transcribed sequences. The actual sequences involved are difficult to identify either by sequence comparisons or by deletion analysis because they are dispersed, additive, and poorly conserved. However, expression-controlling introns can be identified computationally in species with relatively small introns, based on genome-wide differences in oligomer composition between promoter-proximal and distal introns. The genes regulated by introns are often expressed in most tissues and are among the most highly expressed in the genome. The ability of some introns to strongly stimulate mRNA accumulation from several hundred nucleotides downstream of the transcription start site, even when the promoter has been deleted, reveals that our understanding of gene expression remains incomplete. It is unlikely that any diseases are caused by point mutations or small deletions that reduce the expression of an intron-regulated gene unless splicing is also affected. However, introns may be particularly useful in practical applications such as gene therapy because they strongly activate expression but only affect the transcription unit in which they are located.

Pre-mRNA Splicing in the Nuclear Landscape

Eukaryotic gene expression requires the cumulative activity of multiple molecular machines to synthesize and process newly transcribed pre-messenger RNA. Introns, the noncoding regions in pre-mRNA, must be removed by the spliceosome, which assembles on the pre-mRNA as it is transcribed by RNA polymerase II (Pol II). The assembly and activity of the spliceosome can be modulated by features including the speed of transcription elongation, chromatin, post-translational modifications of Pol II and histone tails, and other RNA processing events like 5'-end capping. Here, we review recent work that has revealed cooperation and coordination among co-transcriptional processing events and speculate on new avenues of research. We anticipate new mechanistic insights capable of unraveling the relative contribution of coupled processing to gene expression.

Gender inequalities among authors who contributed equally

We analyzed 2898 scientific papers published between 1995 and 2017 in which two or more authors shared the first author position. For papers in which the first and second authors made equal contributions, mixed-gender combinations were most frequent, followed by male-male and then female-female author combinations. For mixed-gender combinations, more male authors were in the first position, although the disparity decreased over time. For papers in which three or more authors made equal contributions, there were more male authors than female authors in the first position and more all-male than all-female author combinations. The gender inequalities observed among authors who made equal contributions are not consistent with random or alphabetical ordering of authors. These results raise concerns about female authors not receiving proper credit for publications and suggest a need for journals to request clarity on the method used to decide author order among those who contributed equally.

Revealing transcription factor and histone modification co-localization and dynamics across cell lines by integrating ChIP-seq and RNA-seq data

Background

Interactions among transcription factors (TFs) and histone modifications (HMs) play an important role in the precise regulation of gene expression. The context specificity of those interactions and further its dynamics in normal and disease remains largely unknown. Recent development in genomics technology enables transcription profiling by RNA-seq and protein’s binding profiling by ChIP-seq. Integrative analysis of the two types of data allows us to investigate TFs and HMs interactions both from the genome co-localization and downstream target gene expression.

Results

We propose a integrative pipeline to explore the co-localization of 55 TFs and 11 HMs and its dynamics in human GM12878 and K562 by matched ChIP-seq and RNA-seq data from ENCODE. We classify TFs and HMs into three types based on their binding enrichment around transcription start site (TSS). Then a set of statistical indexes are proposed to characterize the TF-TF and TF-HM co-localizations. We found that Rad21, SMC3, and CTCF co-localized across five cell lines. High resolution Hi-C data in GM12878 shows that they associate most of the Hi-C peak loci with a specific CTCF-motif “anchor” and supports that CTCF, SMC3, and RAD2 co-localization serves important role in 3D chromatin structure. Meanwhile, 17 TF-TF pairs are highly dynamic between GM12878 and K562. We then build SVM models to correlate high and low expression level of target genes with TF binding and HM strength. We found that H3k9ac, H3k27ac, and three TFs (ELF1, TAF1, and POL2) are predictive with the accuracy about 85~92%.

Conclusion

We propose a pipeline to analyze the co-localization of TF and HM and their dynamics across cell lines from ChIP-seq, and investigate their regulatory potency by RNA-seq. The integrative analysis of two level data reveals new insight for the cooperation of TFs and HMs and is helpful in understanding cell line specificity of TF/HM interactions.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5278-5) contains supplementary material, which is available to authorized users.

Identifying sequence features that drive ribosomal association for lncRNA

BACKGROUND

With the increasing number of annotated long noncoding RNAs (lncRNAs) from the genome, researchers are continually updating their understanding of lncRNAs. Recently, thousands of lncRNAs have been reported to be associated with ribosomes in mammals. However, their biological functions or mechanisms are still unclear.

RESULTS

In this study, we tried to investigate the sequence features involved in the ribosomal association of lncRNA. We have extracted ninety-nine sequence features corresponding to different biological mechanisms (i.e., RNA splicing, putative ORF, k-mer frequency, RNA modification, RNA secondary structure, and repeat element). An [Formula: see text]-regularized logistic regression model was applied to screen these features. Finally, we obtained fifteen and nine important features for the ribosomal association of human and mouse lncRNAs, respectively.

CONCLUSION

To our knowledge, this is the first study to characterize ribosome-associated lncRNAs and ribosome-free lncRNAs from the perspective of sequence features. These sequence features that were identified in this study may shed light on the biological mechanism of the ribosomal association and provide important clues for functional analysis of lncRNAs.

Examining the Roles of H3K4 Methylation States with Systematically Characterized Antibodies

Histone post-translational modifications (PTMs) are important genomic regulators often studied by chromatin immunoprecipitation (ChIP), whereby their locations and relative abundance are inferred by antibody capture of nucleosomes and associated DNA. However, the specificity of antibodies within these experiments has not been systematically studied. Here, we use histone peptide arrays and internally calibrated ChIP (ICeChIP) to characterize 52 commercial antibodies purported to distinguish the H3K4 methylforms (me1, me2, and me3, with each ascribed distinct biological functions). We find that many widely used antibodies poorly distinguish the methylforms and that high- and low-specificity reagents can yield dramatically different biological interpretations, resulting in substantial divergence from the literature for numerous H3K4 methylform paradigms. Using ICeChIP, we also discern quantitative relationships between enhancer H3K4 methylation and promoter transcriptional output and can measure global PTM abundance changes. Our results illustrate how poor antibody specificity contributes to the "reproducibility crisis," demonstrating the need for rigorous, platform-appropriate validation.

Nucleosome Positioning of Intronless Genes in the Human Genome

Nucleosomes, the basic units of chromatin, are involved in transcription regulation and DNA replication. Intronless genes, which constitute 3 percent of the human genome, differ from intron-containing genes in evolution and function. Our analysis reveals that nucleosome positioning shows a distinct pattern in intronless and intron-containing genes. The nucleosome occupancy upstream of transcription start sites of intronless genes is lower than that of intron-containing genes. In contrast, high occupancy and well positioned nucleosomes are observed along the gene body of intronless genes, which is perfectly consistent with the barrier nucleosome model. Intronless genes have a significantly lower expression level than intron-containing genes and most of them are not expressed in CD4+ T cell lines and GM12878 cell lines, which results from their tissue specificity. However, the highly expressed genes are at the same expression level between the two types of genes. The highly expressed intronless genes require a higher density of RNA Pol II in an elongating state to compensate for the lack of introns. Additionally, 5' and 3' nucleosome depleted regions of highly expressed intronless genes are deeper than those of highly expressed intron-containing genes.

Tales around the clock: Poly(A) tails in circadian gene expression

Circadian rhythms are ubiquitous time-keeping processes in eukaryotes with a period of ~24 hr. Light is perhaps the main environmental cue (zeitgeber) that affects several aspects of physiology and behaviour, such as sleep/wake cycles, orientation of birds and bees, and leaf movements in plants. Temperature can serve as the main zeitgeber in the absence of light cycles, even though it does not lead to rhythmicity through the same mechanism as light. Additional cues include feeding patterns, humidity, and social rhythms. At the molecular level, a master oscillator orchestrates circadian rhythms and organizes molecular clocks located in most cells. The generation of the 24 hr molecular clock is based on transcriptional regulation, as it drives intrinsic rhythmic changes based on interlocked transcription/translation feedback loops that synchronize expression of genes. Thus, processes and factors that determine rhythmic gene expression are important to understand circadian rhythms. Among these, the poly(A) tails of RNAs play key roles in their stability, translational efficiency and degradation. In this article, we summarize current knowledge and discuss perspectives on the role and significance of poly(A) tails and associating factors in the context of the circadian clock. This article is categorized under: RNA Turnover and Surveillance > Regulation of RNA Stability RNA Processing > 3' End Processing.

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.

H3K4 demethylase KDM5B regulates global dynamics of transcription elongation and alternative splicing in embryonic stem cells

Epigenetic regulation of chromatin plays a critical role in controlling embryonic stem (ES) cell self-renewal and pluripotency. However, the roles of histone demethylases and activating histone modifications such as trimethylated histone 3 lysine 4 (H3K4me3) in transcriptional events such as RNA polymerase II (RNAPII) elongation and alternative splicing are largely unknown. In this study, we show that KDM5B, which demethylates H3K4me3, plays an integral role in regulating RNAPII occupancy, transcriptional initiation and elongation, and alternative splicing events in ES cells. Depletion of KDM5B leads to altered RNAPII promoter occupancy, and decreased RNAPII initiation and elongation rates at active genes and at genes marked with broad H3K4me3 domains. Moreover, our results demonstrate that spreading of H3K4me3 from promoter to gene body regions, which is mediated by depletion of KDM5B, modulates RNAPII elongation rates and RNA splicing in ES cells. We further show that KDM5B is enriched nearby alternatively spliced exons, and depletion of KDM5B leads to altered levels of H3K4 methylation in alternatively spliced exon regions, which is accompanied by differential expression of these alternatively splice exons. Altogether, our data indicate an epigenetic role for KDM5B in regulating RNAPII elongation and alternative splicing, which may support the diverse mRNA repertoire in ES cells.

Nuclear retention of mRNAs - quality control, gene regulation and human disease

Nuclear retention of incompletely spliced or mature mRNAs emerges as a novel, previously underappreciated layer of gene regulation, which enables the cell to rapidly respond to stress, viral infection, differentiation cues or changing environmental conditions. Focusing on mammalian cells, we discuss recent insights into the mechanisms and functions of nuclear retention, describe retention-promoting features in protein-coding transcripts and propose mechanisms for their regulated release into the cytoplasm. Moreover, we discuss examples of how aberrant nuclear retention of mRNAs may lead to human diseases.

Defining the location of promoter-associated R-loops at near-nucleotide resolution using bisDRIP-seq

R-loops are features of chromatin consisting of a strand of DNA hybridized to RNA, as well as the expelled complementary DNA strand. R-loops are enriched at promoters where they have recently been shown to have important roles in modifying gene expression. However, the location of promoter-associated R-loops and the genomic domains they perturb to modify gene expression remain unclear. To resolve this issue, we developed a bisulfite-based approach, bisDRIP-seq, to map R-loops across the genome at near-nucleotide resolution in MCF-7 cells. We found the location of promoter-associated R-loops is dependent on the presence of introns. In intron-containing genes, R-loops are bounded between the transcription start site and the first exon-intron junction. In intronless genes, the 3' boundary displays gene-specific heterogeneity. Moreover, intronless genes are often associated with promoter-associated R-loop formation. Together, these studies provide a high-resolution map of R-loops and identify gene structure as a critical determinant of R-loop formation.

RNA sequencing identifies novel non-coding RNA and exon-specific effects associated with cigarette smoking

BACKGROUND

Cigarette smoking is the leading modifiable risk factor for disease and death worldwide. Previous studies quantifying gene-level expression have documented the effect of smoking on mRNA levels. Using RNA sequencing, it is possible to analyze the impact of smoking on complex regulatory phenomena (e.g. alternative splicing, differential isoform usage) leading to a more detailed understanding of the biology underlying smoking-related disease.

METHODS

We used whole-blood RNA sequencing to describe gene and exon-level expression differences between 229 current and 286 former smokers in the COPDGene study. We performed differential gene expression and differential exon usage analyses using the voom/limma and DEXseq R packages. Samples from current and former smokers were compared while controlling for age, gender, race, lifetime smoke exposure, cell counts, and technical covariates.

RESULTS

At an adjusted p-value <0.05, 171 genes were differentially expressed between current and former smokers. Differentially expressed genes included 7 long non-coding RNAs that have not been previously associated with smoking: LINC00599, LINC01362, LINC00824, LINC01624, RP11-563D10.1, RP11-98G13.1, AC004791.2. Secondary analysis of acute smoking (having smoked within 2-h) revealed 5 of the 171 smoking genes demonstrated an acute response above the baseline effect of chronic smoking. Exon-level analyses identified 9 exons from 8 genes with significant differential usage by smoking status, suggesting smoking-induced changes in isoform expression.

CONCLUSIONS

Transcriptomic changes at the gene and exon levels from whole blood can refine our understanding of the molecular mechanisms underlying the response to smoking.

The histone variant H2A.Z promotes splicing of weak introns

In this study, Nissen et al. investigated the function of the highly conserved histone variant H2A.Z in pre-mRNA splicing using the intron-rich model yeast S. pombe. The findings suggest that H2A.Z occupancy promotes cotranscriptional splicing of suboptimal introns that may otherwise be discarded via proofreading ATPases.

Multiple lines of evidence implicate chromatin in the regulation of premessenger RNA (pre-mRNA) splicing. However, the influence of chromatin factors on cotranscriptional splice site usage remains unclear. Here we investigated the function of the highly conserved histone variant H2A.Z in pre-mRNA splicing using the intron-rich model yeast Schizosaccharomyces pombe. Using epistatic miniarray profiles (EMAPs) to survey the genetic interaction landscape of the Swr1 nucleosome remodeling complex, which deposits H2A.Z, we uncovered evidence for functional interactions with components of the spliceosome. In support of these genetic connections, splicing-specific microarrays show that H2A.Z and the Swr1 ATPase are required during temperature stress for the efficient splicing of a subset of introns. Notably, affected introns are enriched for H2A.Z occupancy and more likely to contain nonconsensus splice sites. To test the significance of the latter correlation, we mutated the splice sites in an affected intron to consensus and found that this suppressed the requirement for H2A.Z in splicing of that intron. These data suggest that H2A.Z occupancy promotes cotranscriptional splicing of suboptimal introns that may otherwise be discarded via proofreading ATPases. Consistent with this model, we show that overexpression of splicing ATPase Prp16 suppresses both the growth and splicing defects seen in the absence of H2A.Z.

Splicing and transcription touch base: co-transcriptional spliceosome assembly and function

Several macromolecular machines collaborate to produce eukaryotic messenger RNA. RNA polymerase II (Pol II) translocates along genes that are up to millions of base pairs in length and generates a flexible RNA copy of the DNA template. This nascent RNA harbours introns that are removed by the spliceosome, which is a megadalton ribonucleoprotein complex that positions the distant ends of the intron into its catalytic centre. Emerging evidence that the catalytic spliceosome is physically close to Pol II in vivo implies that transcription and splicing occur on similar timescales and that the transcription and splicing machineries may be spatially constrained. In this Review, we discuss aspects of spliceosome assembly, transcription elongation and other co-transcriptional events that allow the temporal coordination of co-transcriptional splicing.

The determinants of alternative RNA splicing in human cells

Alternative splicing represents an important level of the regulation of gene function in eukaryotic organisms. It plays a critical role in virtually every biological process within an organism, including regulation of cell division and cell death, differentiation of tissues in the embryo and the adult organism, as well as in cellular response to diverse environmental factors. In turn, studies of the last decade have shown that alternative splicing itself is controlled by different mechanisms. Unfortunately, there is no clear understanding of how these diverse mechanisms, or determinants, regulate and constrain the set of alternative RNA species produced from any particular gene in every cell of the human body. Here, we provide a consolidated overview of alternative splicing determinants including RNA-protein interactions, epigenetic regulation via chromatin remodeling, coupling of transcription-to-alternative splicing, effect of secondary structures in pre-RNA, and function of the RNA quality control systems. We also extensively and critically discuss some mechanistic insights on coordinated inclusion/exclusion of exons during the formation of mature RNA molecules. We conclude that the final structure of RNA is pre-determined by a complex interplay between cis- and trans-acting factors. Altogether, currently available empirical data significantly expand our understanding of the functioning of the alternative splicing machinery of cells in normal and pathological conditions. On the other hand, there are still many blind spots that require further deep investigations.

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.

How introns enhance gene expression

In many eukaryotes, including mammals, plants, yeast, and insects, introns can increase gene expression without functioning as a binding site for transcription factors. This phenomenon was termed 'intron-mediated enhancement'. Introns can increase transcript levels by affecting the rate of transcription, nuclear export, and transcript stability. Moreover, introns can also increase the efficiency of mRNA translation. This review discusses the current knowledge about these processes. The role of splicing in IME and the significance of intron position relative to the sites of transcription and translation initiation are elaborated. Particular emphasis is placed on the question why different introns, present at the same location of the same genes and spliced at a similar high efficiency, can have very different impacts on expression - from almost no effect to considerable stimulation. This situation can be at least partly accounted for by the identification of splicing-unrelated intronic elements with a special ability to enhance mRNA accumulation or translational efficiency. The many factors that could lead to the large variation observed between the impact of introns in different genes and experimental systems are highlighted. It is suggested that there is no sole, definite answer to the question "how do introns enhance gene expression". Rather, each intron-gene combination might undergo its own unique mixture of processes that lead to the perceptible outcome.

Single-cell, single-mRNA analysis of Ccnb1 promoter regulation

Promoter activation drives gene transcriptional output. Here we report generating site-specifically integrated single-copy promoter transgenes and measuring their expression to indicate promoter activities at single-mRNA level. mRNA counts, Pol II density and Pol II firing rates of the Ccnb1 promoter transgene resembled those of the native Ccnb1 gene both among asynchronous cells and during the cell cycle. We observed distinct activation states of the Ccnb1 promoter among G1 and G2/M cells, suggesting cell cycle-independent origin of cell-to-cell variation in Ccnb1 promoter activation. Expressing a dominant-negative mutant of NF-YA, a key transcriptional activator of the Ccnb1 promoter, increased its “OFF”/“ON” time ratios but did not alter Pol II firing rates during the “ON” period. Furthermore, comparing H3K4me2 and H3K79me2 levels at the Ccnb1 promoter transgene and the native Ccnb1 gene indicated that the enrichment of these two active histone marks did not predispose higher transcriptional activities. In summary, this experimental system enables bridging transcription imaging with molecular analysis to provide novel insights into eukaryotic transcriptional regulation.

A saga of cancer epigenetics: linking epigenetics to alternative splicing

The discovery of an increasing number of alternative splicing events in the human genome highlighted that ∼94% of genes generate alternatively spliced transcripts that may produce different protein isoforms with diverse functions. It is now well known that several diseases are a direct and indirect consequence of aberrant splicing events in humans. In addition to the conventional mode of alternative splicing regulation by 'cis' RNA-binding sites and 'trans' RNA-binding proteins, recent literature provides enormous evidence for epigenetic regulation of alternative splicing. The epigenetic modifications may regulate alternative splicing by either influencing the transcription elongation rate of RNA polymerase II or by recruiting a specific splicing regulator via different chromatin adaptors. The epigenetic alterations and aberrant alternative splicing are known to be associated with various diseases individually, but this review discusses/highlights the latest literature on the role of epigenetic alterations in the regulation of alternative splicing and thereby cancer progression. This review also points out the need for further studies to understand the interplay between epigenetic modifications and aberrant alternative splicing in cancer progression.

PHF13: A new player involved in RNA polymerase II transcriptional regulation and co-transcriptional splicing

We recently identified PHF13 as an H3K4me2/3 chromatin reader and transcriptional co-regulator. We found that PHF13 interacts with RNAPIIS5P and PRC2 stabilizing their association with active and bivalent promoters. Furthermore, mass spectrometry analysis identified ∼50 spliceosomal proteins in PHF13s interactome. Here, we will discuss the potential role of PHF13 in RNAPII pausing and co-transcriptional splicing.

Intron-Mediated Enhancement: A Tool for Heterologous Gene Expression in Plants?

Many plant promoters were characterized and used for transgene expression in plants. Even though these promoters drive high levels of transgene expression in plants, the expression patterns are rarely constitutive but restricted to some tissues and developmental stages. In terms of crop improvement not only the enhancement of expression per se but, in particular, tissue-specific and spatial expression of genes plays an important role. Introns were used to boost expression in transgenic plants in the field of crop improvement for a long time. However, the mechanism behind this so called intron-mediated enhancement (IME) is still largely unknown. This review highlights the complexity of IME on the levels of its regulation and modes of action and gives an overview on IME methodology, examples in fundamental research and models of proposed mechanisms. In addition, the application of IME in heterologous gene expression is discussed.

DGW: an exploratory data analysis tool for clustering and visualisation of epigenomic marks

Background

Functional genomic and epigenomic research relies fundamentally on sequencing based methods like ChIP-seq for the detection of DNA-protein interactions. These techniques return large, high dimensional data sets with visually complex structures, such as multi-modal peaks extended over large genomic regions. Current tools for visualisation and data exploration represent and leverage these complex features only to a limited extent.

Results

We present DGW, an open source software package for simultaneous alignment and clustering of multiple epigenomic marks. DGW uses Dynamic Time Warping to adaptively rescale and align genomic distances which allows to group regions of interest with similar shapes, thereby capturing the structure of epigenomic marks. We demonstrate the effectiveness of the approach in a simulation study and on a real epigenomic data set from the ENCODE project.

Conclusions

Our results show that DGW automatically recognises and aligns important genomic features such as transcription start sites and splicing sites from histone marks. DGW is available as an open source Python package.

The 5'UTR Intron of Arabidopsis GGT1 Aminotransferase Enhances Promoter Activity by Recruiting RNA Polymerase II

Photorespiration is essential for the detoxification of glycolate and recycling of carbon to the Calvin Benson Bassham cycle. Enzymes participating in the pathway have been identified, and investigations now focus on the regulation of photorespiration by transporters and metabolites. However, regulation of photorespiration on the gene level has not been intensively studied. Here, we show that maximum transcript abundance of Glu:glyoxylate aminotransferase 1 (GGT1) is regulated by intron-mediated enhancement (IME) of the 5' leader intron rather than by regulatory elements in the 5' upstream region. The intron is rich in CT-stretches and contains the motif TGTGATTTG that is highly similar to the IME-related motif TTNGATYTG. The GGT1 intron also confers leaf-specific expression of foreign promoters. Quantitative PCR analysis and GUS activity measurements revealed that IME of the GGT1 5'UTR intron is controlled on the transcriptional level. IME by the GGT1 5'UTR intron was at least 2-fold. Chromatin immunoprecipitation experiments showed that the abundance of RNA polymerase II binding to the intron-less construct is reduced.

TALE-directed local modulation of H3K9 methylation shapes exon recognition

In search for the function of local chromatin environment on pre-mRNA processing we established a new tool, which allows for the modification of chromatin using a targeted approach. Using Transcription Activator-Like Effector domains fused to histone modifying enzymes (TALE-HME), we show locally restricted alteration of histone methylation modulates the splicing of target exons. We provide evidence that a local increase in H3K9 di- and trimethylation promotes inclusion of the target alternative exon, while demethylation by JMJD2D leads to exon skipping. We further demonstrate that H3K9me3 is localized on internal exons genome-wide suggesting a general role in splicing. Consistently, targeting of the H3K9 demethylase to a weak constitutive exon reduced co-transcriptional splicing. Together our data show H3K9 methylation within the gene body is a factor influencing recognition of both constitutive and alternative exons.

The chicken erythrocyte epigenome

BACKGROUND

Transcriptional regulation is impacted by multiple layers of genome organization. A general feature of transcriptionally active chromatin is sensitivity to DNase I and association with acetylated histones. However, very few of these active DNase I-sensitive domains, such as the chicken erythrocyte β-globin domain, have been identified and characterized. In chicken polychromatic erythrocytes, dynamically acetylated histones associated with DNase I-sensitive, transcriptionally active chromatin prevent histone H1/H5-induced insolubility at physiological ionic strength.

RESULTS

Here, we identified and mapped out all the transcriptionally active chromosomal domains in the chicken polychromatic erythrocyte genome by combining a powerful chromatin fractionation method with next-generation DNA and RNA sequencing. Two classes of transcribed chromatin organizations were identified on the basis of the extent of solubility at physiological ionic strength. Highly transcribed genes were present in multigenic salt-soluble chromatin domains ranging in length from 30 to over 150 kb. We identified over 100 highly expressed genes that were organized in broad dynamically highly acetylated, salt-soluble chromatin domains. Highly expressed genes were associated with H3K4me3 and H3K27ac and produced discernible antisense transcripts. The moderately- and low-expressing genes had highly acetylated, salt-soluble chromatin regions confined to the 5' end of the gene.

CONCLUSIONS

Our data provide a genome-wide profile of chromatin signatures in relation to expression levels in chicken polychromatic erythrocytes.

Splicing of Nascent RNA Coincides with Intron Exit from RNA Polymerase II

Protein-coding genes in eukaryotes are transcribed by RNA polymerase II (Pol II) and introns are removed from pre-mRNA by the spliceosome. Understanding the time lag between Pol II progression and splicing could provide mechanistic insights into the regulation of gene expression. Here, we present two single-molecule nascent RNA sequencing methods that directly determine the progress of splicing catalysis as a function of Pol II position. Endogenous genes were analyzed on a global scale in budding yeast. We show that splicing is 50% complete when Pol II is only 45 nt downstream of introns, with the first spliced products observed as introns emerge from Pol II. Perturbations that slow the rate of spliceosome assembly or speed up the rate of transcription caused splicing delays, showing that regulation of both processes determines in vivo splicing profiles. We propose that matched rates streamline the gene expression pathway, while allowing regulation through kinetic competition.

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.

Dynamic Histone Acetylation of H3K4me3 Nucleosome Regulates MCL1 Pre-mRNA Splicing

Pre-mRNA splicing is a cotranscriptional process affected by the chromatin architecture along the body of coding genes. Recruited to the pre-mRNA by splicing factors, histone deacetylases (HDACs) and K-acetyltransferases (KATs) catalyze dynamic histone acetylation along the gene. In colon carcinoma HCT 116 cells, HDAC inhibition specifically increased KAT2B occupancy as well as H3 and H4 acetylation of the H3K4 trimethylated (H3K4me3) nucleosome positioned over alternative exon 2 of the MCL1 gene, an event paralleled with the exclusion of exon 2. These results were reproduced in MDA-MB-231, but not in MCF7 breast adenocarcinoma cells. These later cells have much higher levels of demethylase KDM5B than either HCT 116 or MDA-MB-231 cells. We show that H3K4me3 steady-state levels and H3K4me3 occupancy at the end of exon 1 and over exon 2 of the MCL1 gene were lower in MCF7 than in MDA-MB-231 cells. Furthermore, in MCF7 cells, there was minimal effect of HDAC inhibition on H3/H4 acetylation and H3K4me3 levels along the MCL1 gene and no change in pre-mRNA splicing choice. These results show that, upon HDAC inhibition, the H3K4me3 mark plays a critical role in the exclusion of exon 2 from the MCL1 pre-mRNA. J. Cell. Physiol. 231: 2196-2204, 2016. © 2016 Wiley Periodicals, Inc.

Histone H3K4 trimethylation: dynamic interplay with pre-mRNA splicing

Histone H3 lysine 4 trimethylation (H3K4me3) is often stated as a mark of transcriptionally active promoters. However, closer study of the positioning of H3K4me3 shows the mark locating primarily after the first exon at the 5′ splice site and overlapping with a CpG island in mammalian cells. There are several enzyme complexes that are involved in the placement of the H3K4me3 mark, including multiple protein complexes containing SETD1A, SETD1B, and MLL1 enzymes (writers). CXXC1, which is associated with SETD1A and SETD1B, target these enzymes to unmethylated CpG islands. Lysine demethylases (KDM5 family members, erasers) demethylate H3K4me3. The H3K4me3 mark is recognized by several proteins (readers), including lysine acetyltransferase complexes, chromatin remodelers, and RNA bound proteins involved in pre-mRNA splicing. Interestingly, attenuation of H3K4me3 impacts pre-mRNA splicing, and inhibition of pre-mRNA splicing attenuates H3K4me3.

Discover hidden splicing variations by mapping personal transcriptomes to personal genomes

RNA-seq has become a popular technology for studying genetic variation of pre-mRNA alternative splicing. Commonly used RNA-seq aligners rely on the consensus splice site dinucleotide motifs to map reads across splice junctions. Consequently, genomic variants that create novel splice site dinucleotides may produce splice junction RNA-seq reads that cannot be mapped to the reference genome. We developed and evaluated an approach to identify ‘hidden’ splicing variations in personal transcriptomes, by mapping personal RNA-seq data to personal genomes. Computational analysis and experimental validation indicate that this approach identifies personal specific splice junctions at a low false positive rate. Applying this approach to an RNA-seq data set of 75 individuals, we identified 506 personal specific splice junctions, among which 437 were novel splice junctions not documented in current human transcript annotations. 94 splice junctions had splice site SNPs associated with GWAS signals of human traits and diseases. These involve genes whose splicing variations have been implicated in diseases (such as OAS1), as well as novel associations between alternative splicing and diseases (such as ICA1). Collectively, our work demonstrates that the personal genome approach to RNA-seq read alignment enables the discovery of a large but previously unknown catalog of splicing variations in human populations.

The enduring mystery of intron-mediated enhancement

Within two years of their discovery in 1977, introns were found to have a positive effect on gene expression. Numerous examples of stimulatory introns have been described since then in very diverse organisms, including plants. In some cases, the mechanism through which the intron affects expression is readily understood. However, many introns that affect expression increase mRNA accumulation through an unknown mechanism, referred to as intron-mediated enhancement (IME). Despite several decades of research into IME, and the clear benefits of using introns to increase transgene expression, little progress has been made in understanding the mechanism of IME. Several fundamental questions regarding the role of transcription and splicing, the sequences responsible for IME, the involvement of other factors, and the relationship between introns and promoters remain unanswered. The more we learn about the properties of stimulating introns, the clearer it becomes that the effects of introns are unfamiliar and difficult to reconcile with conventional views of how transcription is controlled. We hypothesize that introns increase transcript initiation upstream of themselves by creating a localized region of accessible chromatin. Introns might represent a novel kind of downstream regulatory element for genes transcribed by RNA polymerase II.