Prominent use of distal 5' transcription start sites and discovery of a large number of additional exons in ENCODE regions

The Evolution of Ultraconserved Elements in Vertebrates

Ultraconserved elements were discovered two decades ago, arbitrarily defined as sequences that are identical over a length ≥ 200 bp in the human, mouse, and rat genomes. The definition was subsequently extended to sequences ≥ 100 bp identical in at least three of five mammalian genomes (including dog and cow), and shown to have undergone rapid expansion from ancestors in fish and strong negative selection in birds and mammals. Since then, many more genomes have become available, allowing better definition and more thorough examination of ultraconserved element distribution and evolutionary history. We developed a fast and flexible analytical pipeline for identifying ultraconserved elements in multiple genomes, dedUCE, which allows manipulation of minimum length, sequence identity, and number of species with a detectable ultraconserved element according to specified parameters. We suggest an updated definition of ultraconserved elements as sequences ≥ 100 bp and ≥97% sequence identity in ≥50% of placental mammal orders (12,813 ultraconserved elements). By mapping ultraconserved elements to ∼200 species, we find that placental ultraconserved elements appeared early in vertebrate evolution, well before land colonization, suggesting that the evolutionary pressures driving ultraconserved element selection were present in aquatic environments in the Cambrian-Devonian periods. Most (>90%) ultraconserved elements likely appeared after the divergence of gnathostomes from jawless predecessors, were largely established in sequence identity by early Sarcopterygii evolution-before the divergence of lobe-finned fishes from tetrapods-and became near fixed in the amniotes. Ultraconserved elements are mainly located in the introns of protein-coding and noncoding genes involved in neurological and skeletomuscular development, enriched in regulatory elements, and dynamically expressed throughout embryonic development.

Alternative isoform expression of key thermogenic genes in human beige adipocytes

Background

The beneficial effect of thermogenic adipocytes in maintaining body weight and protecting against metabolic disorders has raised interest in understanding the regulatory mechanisms defining white and beige adipocyte identity. Although alternative splicing has been shown to propagate adipose browning signals in mice, this has yet to be thoroughly investigated in human adipocytes.

Methods

We performed parallel white and beige adipogenic differentiation using primary adipose stem cells from 6 unrelated healthy subjects and assessed differential gene and isoform expression in mature adipocytes by RNA sequencing.

Results

We find 777 exon junctions with robust differential usage between white and beige adipocytes in all 6 subjects, mapping to 562 genes. Importantly, only 10% of these differentially spliced genes are also differentially expressed, indicating that alternative splicing constitutes an additional layer of gene expression regulation during beige adipocyte differentiation. Functional classification of alternative isoforms points to a gain of function for key thermogenic transcription factors such as PPARG and CITED1, and enzymes such as PEMT, or LPIN1. We find that a large majority of the splice variants arise from differential TSS usage, with beige-specific TSSs being enriched for PPARγ and MED1 binding compared to white-specific TSSs. Finally, we validate beige specific isoform expression at the protein level for two thermogenic regulators, PPARγ and PEMT.

Discussion

These results suggest that differential isoform expression through alternative TSS usage is an important regulatory mechanism for human adipocyte thermogenic specification.

RTCpredictor: identification of read-through chimeric RNAs from RNA sequencing data

Read-through chimeric RNAs are being recognized as a means to expand the functional transcriptome and contribute to cancer tumorigenesis when mis-regulated. However, current software tools often fail to predict them. We have developed RTCpredictor, utilizing a fast ripgrep tool to search for all possible exon-exon combinations of parental gene pairs. We also added exonic variants allowing searches containing common SNPs. To our knowledge, it is the first read-through chimeric RNA specific prediction method that also provides breakpoint coordinates. Compared with 10 other popular tools, RTCpredictor achieved high sensitivity on a simulated and three real datasets. In addition, RTCpredictor has less memory requirements and faster execution time, making it ideal for applying on large datasets.

Evidence for widespread existence of functional novel and non-canonical human transcripts

BACKGROUND

Fraction of functional sequence in the human genome remains a key unresolved question in Biology and the subject of vigorous debate. While a plethora of studies have connected a significant fraction of human DNA to various biochemical processes, the classical definition of function requires evidence of effects on cellular or organismal fitness that such studies do not provide. Although multiple high-throughput reverse genetics screens have been developed to address this issue, they are limited to annotated genomic elements and suffer from non-specific effects, arguing for a strong need to develop additional functional genomics approaches.

RESULTS

In this work, we established a high-throughput lentivirus-based insertional mutagenesis strategy as a forward genetics screen tool in aneuploid cells. Application of this approach to human cell lines in multiple phenotypic screens suggested the presence of many yet uncharacterized functional elements in the human genome, represented at least in part by novel exons of known and novel genes. The novel transcripts containing these exons can be massively, up to thousands-fold, induced by specific stresses, and at least some can represent bi-cistronic protein-coding mRNAs.

CONCLUSIONS

Altogether, these results argue that many unannotated and non-canonical human transcripts, including those that appear as aberrant splice products, have biological relevance under specific biological conditions.

Complexities of mammalian transcriptome revealed by targeted RNA enrichment techniques

Studies using highly sensitive targeted RNA enrichment methods have shown that a large portion of the human transcriptome remains to be discovered and that most of the genome is transcribed in a complex, interleaved fashion characterized by a complex web of transcripts emanating from protein coding and noncoding loci. These results resonate with those from single-cell transcriptome profiling endeavors that reveal the existence of multiple novel, cell type-specific transcripts and clearly demonstrate that our understanding of the complexities of the human transcriptome is far from being complete. Here, we review the current status of the targeted RNA enrichment techniques, their application to the discovery of novel cell type-specific transcripts, and their impact on our understanding of the human genome and transcriptome.

RTCpredictor: Identification of Read-Through Chimeric RNAs from RNA Sequencing Data

Read-through chimeric RNAs are gaining attention in cancer and other research fields, yet current tools often fail in predicting them. We have thus developed the first read-through chimeric RNA specific prediction method, RTCpredictor, utilizing a fast ripgrep algorithm to search for all possible exon-exon combinations of parental gene pairs. Compared with other ten popular tools, RTCpredictor achieved top performance on both simulated and real datasets. We randomly selected up to 30 candidate read-through chimeras predicted from each software method and experimentally validated a total of 109 read-throughs and on this set, RTCpredictor outperformed all the other methods. In addition, RTCpredictor ( https://github.com/sandybioteck/RTCpredictor ) has less memory requirements and faster execution time.

Hotspots of single-strand DNA “breakome” are enriched at transcriptional start sites of genes

Single-strand breaks (SSBs) represent one of the most common types of DNA damage, yet not much is known about the genome landscapes of this type of DNA lesions in mammalian cells. Here, we found that SSBs are more likely to occur in certain positions of the human genome—SSB hotspots—in different cells of the same cell type and in different cell types. We hypothesize that the hotspots are likely to represent biologically relevant breaks. Furthermore, we found that the hotspots had a prominent tendency to be enriched in the immediate vicinity of transcriptional start sites (TSSs). We show that these hotspots are not likely to represent technical artifacts or be caused by common mechanisms previously found to cause DNA cleavage at promoters, such as apoptotic DNA fragmentation or topoisomerase type II (TOP2) activity. Therefore, such TSS-associated hotspots could potentially be generated using a novel mechanism that could involve preferential cleavage at cytosines, and their existence is consistent with recent studies suggesting a complex relationship between DNA damage and regulation of gene expression.

AZFa Y gene, DDX3Y, evolved novel testis transcript variants in primates with proximal 3´UTR polyadenylation for germ cell specific translation

Translational control is a major level of gene expression regulation in the male germ line. DDX3Y located in the AZFa region of the human Y chromosome encodes a conserved RNA helicase important for translational control at the G1-S phase of the cell cycle. In human, DDX3Y protein is expressed only in premeiotic male germ cells. In primates, DDX3Y evolved a second promoter producing novel testis-specific transcripts. Here, we show primate species-specific use of alternative polyadenylation (APA) sites for these testis-specific DDX3Y transcript variants. They have evolved subsequently in the 3´UTRs of the primates´ DDX3Y transcripts. Whereas a distal APA site (PAS4) is still used for polyadenylation of most DDX3Y testis transcripts in Callithrix jacchus; two proximal APAs (PAS1; PAS2) are used predominantly in Macaca mulatta, in Pan trogloydates and in human. This shift corresponds with a significant increase of DDX3Y protein expression in the macaque testis tissue. In chimpanzee and human, shift to predominant use of the most proximal APA site (PAS1) is associated with translation of these DDX3Y transcripts in only premeiotic male germ cells. We therefore assume evolution of a positive selection process for functional DDX3Y testis transcripts in these primates which increase their stability and translation efficiency to promote its cell cycle balancing function in the human male germ line.

Very long intergenic non-coding (vlinc) RNAs directly regulate multiple genes in cis and trans

BACKGROUND

The majority of the human genome is transcribed in the form of long non-coding (lnc) RNAs. While these transcripts have attracted considerable interest, their molecular mechanisms of function and biological significance remain controversial. One of the main reasons behind this lies in the significant challenges posed by lncRNAs requiring the development of novel methods and concepts to unravel their functionality. Existing methods often lack cross-validation and independent confirmation by different methodologies and therefore leave significant ambiguity as to the authenticity of the outcomes. Nonetheless, despite all the caveats, it appears that lncRNAs may function, at least in part, by regulating other genes via chromatin interactions. Therefore, the function of a lncRNA could be inferred from the function of genes it regulates. In this work, we present a genome-wide functional annotation strategy for lncRNAs based on identification of their regulatory networks via the integration of three distinct types of approaches: co-expression analysis, mapping of lncRNA-chromatin interactions, and assaying molecular effects of lncRNA knockdowns obtained using an inducible and highly specific CRISPR/Cas13 system.

RESULTS

We applied the strategy to annotate 407 very long intergenic non-coding (vlinc) RNAs belonging to a novel widespread subclass of lncRNAs. We show that vlincRNAs indeed appear to regulate multiple genes encoding proteins predominantly involved in RNA- and development-related functions, cell cycle, and cellular adhesion via a mechanism involving proximity between vlincRNAs and their targets in the nucleus. A typical vlincRNAs can be both a positive and negative regulator and regulate multiple genes both in trans and cis. Finally, we show vlincRNAs and their regulatory networks potentially represent novel components of DNA damage response and are functionally important for the ability of cancer cells to survive genotoxic stress.

CONCLUSIONS

This study provides strong evidence for the regulatory role of the vlincRNA class of lncRNAs and a potentially important role played by these transcripts in the hidden layer of RNA-based regulation in complex biological systems.

Gene Fusions Derived by Transcriptional Readthrough are Driven by Segmental Duplication in Human

Gene fusion occurs when two or more individual genes with independent open reading frames becoming juxtaposed under the same open reading frame creating a new fused gene. A small number of gene fusions described in detail have been associated with novel functions, for example, the hominid-specific PIPSL gene, TNFSF12, and the TWE-PRIL gene family. We use Sequence Similarity Networks and species level comparisons of great ape genomes to identify 45 new genes that have emerged by transcriptional readthrough, that is, transcription-derived gene fusion. For 35 of these putative gene fusions, we have been able to assess available RNAseq data to determine whether there are reads that map to each breakpoint. A total of 29 of the putative gene fusions had annotated transcripts (9/29 of which are human-specific). We carried out RT-qPCR in a range of human tissues (placenta, lung, liver, brain, and testes) and found that 23 of the putative gene fusion events were expressed in at least one tissue. Examining the available ribosome foot-printing data, we find evidence for translation of three of the fused genes in human. Finally, we find enrichment for transcription-derived gene fusions in regions of known segmental duplication in human. Together, our results implicate chromosomal structural variation brought about by segmental duplication with the emergence of novel transcripts and translated protein products.

ChiTaRS 5.0: the comprehensive database of chimeric transcripts matched with druggable fusions and 3D chromatin maps

Chimeric RNA transcripts are formed when exons from two genes fuse together, often due to chromosomal translocations, transcriptional errors or trans-splicing effect. While these chimeric RNAs produce functional proteins only in certain cases, they play a significant role in disease phenotyping and progression. ChiTaRS 5.0 (http://chitars.md.biu.ac.il/) is the latest and most comprehensive chimeric transcript repository, with 111 582 annotated entries from eight species, including 23 167 known human cancer breakpoints. The database includes unique information correlating chimeric breakpoints with 3D chromatin contact maps, generated from public datasets of chromosome conformation capture techniques (Hi-C). In this update, we have added curated information on druggable fusion targets matched with chimeric breakpoints, which are applicable to precision medicine in cancers. The introduction of a new section that lists chimeric RNAs in various cell-lines is another salient feature. Finally, using text-mining techniques, novel chimeras in Alzheimer's disease, schizophrenia, dyslexia and other diseases were collected in ChiTaRS. Thus, this improved version is an extensive catalogue of chimeras from multiple species. It extends our understanding of the evolution of chimeric transcripts in eukaryotes and contributes to the analysis of 3D genome conformational changes and the functional role of chimeras in the etiopathogenesis of cancers and other complex diseases.

Regulation of mTOR signaling by long non-coding RNA

The mechanistic target of rapamycin (mTOR) is a major signaling hub that coordinates cellular and organismal responses, such as cell growth, proliferation, apoptosis, and metabolism. Dysregulation of mTOR signaling occurs in many human diseases, and there are significant ongoing efforts to pharmacologically target this pathway. Long noncoding RNAs (lncRNA), defined by a length > 200 nucleotides and absence of a long open-reading-frame, are a class of non-protein-coding RNAs. Mutations and dysregulations of lncRNAs are directly linked to the development and progression of many diseases, including cancer, diabetes, and neurologic disorders. Recent findings reveal diverse functions for lncRNA that include transcriptional regulation, organization of nuclear domains, and regulation of proteins or RNA molecules. Despite considerable development in our understanding of lncRNA over the past decade, only a fraction of annotated lncRNAs has been examined for biological function. In addition, lncRNAs have emerged as therapeutic targets due to their ability to modulate multiple pathways, including mTOR signaling. This review will provide an up-to-date summary of lncRNAs that are involved in regulating mTOR pathway.

A methyl-sensitive element induces bidirectional transcription in TATA-less CpG island-associated promoters

How TATA-less promoters such as those within CpG islands (CGI) control gene expression is still a subject of active research. Here, we have identified the "CGCG element", a ten-base pair motif with a consensus sequence of TCTCGCGAGA present in a group of promoter-associated CGI-enriched in ribosomal protein and housekeeping genes. This element is evolutionarily conserved in vertebrates, found in DNase-accessible regions and employs RNA Pol II to activate gene expression. Through analysis of capped-nascent transcripts and supporting evidence from reporter assays, we demonstrate that this element activates bidirectional transcription through divergent start sites. Methylation of this element abrogates the associated promoter activity. When coincident with a TATA-box, directional transcription remains CGCG-dependent. Because the CGCG element is sufficient to drive transcription, we propose that its unmethylated form functions as a heretofore undescribed promoter element of a group of TATA-less CGI-associated promoters.

Decoding the transcriptional programs activated by psychotropic drugs in the brain

Analysis of drug-induced gene expression in the brain has long held the promise of revealing the molecular mechanisms of drug actions as well as predicting their long-term clinical efficacy. However, despite some successes, this promise has yet to be fulfilled. Here, we present an overview of the current state of understanding of drug-induced gene expression in the brain and consider the obstacles to achieving a robust prediction of the properties of psychoactive compounds based on gene expression profiles. We begin with a comprehensive overview of the mechanisms controlling drug-inducible transcription and the complexity resulting from expression of noncoding RNAs and alternative gene isoforms. Particular interest is placed on studies that examine the associations within drug classes with regard to the effects on gene transcription, alterations in cell signaling and neuropharmacological drug properties. While the ability of gene expression signatures to distinguish specific clinical classes of psychotropic and addictive drugs remains unclear, some reports show that under specific constraints, drug properties can be predicted based on gene expression. Such signatures offer a simple and effective way to classify psychotropic drugs and screen novel psychoactive compounds. Finally, we note that the amount of data regarding molecular programs activated in the brain by drug treatment has grown exponentially in recent years and that future advances may therefore come in large part from integrating the currently available high-throughput data sets.

Gene Regulatory Network Perturbation by Genetic and Epigenetic Variation

Gene regulatory networks underlie biological function and cellular physiology. Alternative splicing (AS) is a fundamental step in gene regulatory networks and plays a key role in development and disease. In addition to the identification of aberrant AS events, an increasing number of studies are focusing on molecular determinants of AS, including genetic and epigenetic regulators. We review here recent efforts to identify various deregulated AS events as well as their molecular determinants that alter biological functions, and discuss clinical features of AS and their druggable potential.

Modeling Transcriptional Rewiring in Neutrophils Through the Course of Treated Juvenile Idiopathic Arthritis

Neutrophils in children with the polyarticular form of juvenile idiopathic arthritis (JIA) display abnormal transcriptional patterns linked to fundamental metabolic derangements. In this study, we sought to determine the effects of therapy on mRNA and miRNA expression networks in polyarticular JIA. Using exon and miRNA microarrays, we studied children with untreated active JIA (ADU, n = 35), children with active disease on therapy with methotrexate ± etanercept (ADT, n = 26), and children with inactive disease also on therapy (ID, n = 14). We compared the results to findings from healthy control children (HC, n = 35). We found substantial re-ordering of mRNA and miRNA expression networks after the initiation of therapy. Each disease state was associated with a distinct transcriptional profile, with the ADT state differing the most from HC, and ID more strongly resembling HC. Changes at the mRNA level were mirrored in changes in miRNA expression patterns. The analysis of the expression dynamics from differentially expressed genes across three disease states indicated that therapeutic response is a complex process. This process does not simply involve genes slowly correcting in a linear fashion over time. Computational modeling of miRNA and transcription factor (TF) co-regulatory networks demonstrated that combinational regulation of miRNA and TF might play an important role in dynamic transcriptome changes.

Transcriptome sequencing of the human pathogen Corynebacterium diphtheriae NCTC 13129 provides detailed insights into its transcriptional landscape and into DtxR-mediated transcriptional regulation

Background

The human pathogen Corynebacterium diphtheriae is the causative agent of diphtheria. In the 1990s a large diphtheria outbreak in Eastern Europe was caused by the strain C. diphtheriae NCTC 13129. Although the genome was sequenced more than a decade ago, not much is known about its transcriptome. Our aim was to use transcriptome sequencing (RNA-Seq) to close this knowledge gap and gain insights into the transcriptional landscape of a C. diphtheriae tox⁺ strain.

Results

We applied two different RNA-Seq techniques, one to retrieve 5′-ends of primary transcripts and the other to characterize the whole transcriptional landscape in order to gain insights into various features of the C. diphtheriae NCTC 13129 transcriptome. By examining the data we identified 1656 transcription start sites (TSS), of which 1202 were assigned to genes and 454 to putative novel transcripts. By using the TSS data promoter regions recognized by the housekeeping sigma factor σ^A and its motifs were analyzed in detail, revealing a well conserved −10 but an only weakly conserved −35 motif, respectively. Furthermore, with the TSS data 5’-UTR lengths were explored. The observed 5’-UTRs range from zero length (leaderless transcripts), which make up 20% of all genes, up to over 450 nt long leaders, which may harbor regulatory functions. The C. diphtheriae transcriptome consists of 471 operons which are further divided into 167 sub-operon structures. In a differential expression analysis approach, we discovered that genetic disruption of the iron-sensing transcription regulator DtxR, which controls expression of diphtheria toxin (DT), causes a strong influence on general gene expression. Nearly 15% of the genome is differentially transcribed, indicating that DtxR might have other regulatory functions in addition to regulation of iron metabolism and DT. Furthermore, our findings shed light on the transcriptional landscape of the DT encoding gene tox and present evidence for two tox antisense RNAs, which point to a new way of transcriptional regulation of toxin production.

Conclusions

This study presents extensive insights into the transcriptome of C. diphtheriae and provides a basis for future studies regarding gene characterization, transcriptional regulatory networks, and regulation of the tox gene in particular.

Electronic supplementary material

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

RNAseq analysis of α-proteobacterium Gluconobacter oxydans 621H

Background

The acetic acid bacterium Gluconobacter oxydans 621H is characterized by its exceptional ability to incompletely oxidize a great variety of carbohydrates in the periplasm. The metabolism of this α-proteobacterium has been characterized to some extent, yet little is known about its transcriptomes and related data. In this study, we applied two different RNAseq approaches. Primary transcriptomes enriched for 5′-ends of transcripts were sequenced to detect transcription start sites, which allow subsequent analysis of promoter motifs, ribosome binding sites, and 5´-UTRs. Whole transcriptomes were sequenced to identify expressed genes and operon structures.

Results

Sequencing of primary transcriptomes of G. oxydans revealed 2449 TSSs, which were classified according to their genomic context followed by identification of promoter and ribosome binding site motifs, analysis of 5´-UTRs including validation of predicted cis-regulatory elements and correction of start codons. 1144 (41%) of all genes were found to be expressed monocistronically, whereas 1634 genes were organized in 571 operons. Together, TSSs and whole transcriptome data were also used to identify novel intergenic (18), intragenic (328), and antisense transcripts (313).

Conclusions

This study provides deep insights into the transcriptional landscapes of G. oxydans. The comprehensive transcriptome data, which we made publicly available, facilitate further analysis of promoters and other regulatory elements. This will support future approaches for rational strain development and targeted gene expression in G. oxydans. The corrections of start codons further improve the high quality genome reference and support future proteome analysis.

Electronic supplementary material

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

Difficulty in obtaining the complete mRNA coding sequence at 5' region (5' end mRNA artifact): Causes, consequences in biology and medicine and possible solutions for obtaining the actual amino acid sequence of proteins (Review)

The known difficulty in obtaining the actual full length, complete sequence of a messenger RNA (mRNA) may lead to the erroneous determination of its coding sequence at the 5' region (5' end mRNA artifact), and consequently to the wrong assignment of the translation start codon, leading to the inaccurate prediction of the encoded polypeptide at its amino terminus. Among the known human genes whose study was affected by this artifact, we can include disco interacting protein 2 homolog A (DIP2A; KIAA0184), Down syndrome critical region 1 (DSCR1), SON DNA binding protein (SON), trefoil factor 3 (TFF3) and URB1 ribosome biogenesis 1 homolog (URB1; KIAA0539) on chromosome 21, as well as receptor for activated C kinase 1 (RACK1, also known as GNB2L1), glutaminyl‑tRNA synthetase (QARS) and tyrosyl-DNA phosphodiesterase 2 (TDP2) along with another 474 loci, including interleukin 16 (IL16). In this review, we discuss the causes of this issue, its quantitative incidence in biomedical research, the consequences in biology and medicine, and the possible solutions for obtaining the actual amino acid sequence of proteins in the post-genomics era.

Novel RNA variants in colorectal cancers

With an annual estimated incidence of 1.4 million, and a five-year survival rate of 60%, colorectal cancer (CRC) is a major clinical burden. To identify novel RNA variants in CRC, we analyzed exon-level microarray expression data from a cohort of 202 CRCs. We nominated 25 genes with increased expression of their 3′ parts in at least one cancer sample each. To efficiently investigate underlying transcript structures, we developed an approach using rapid amplification of cDNA ends followed by high throughput sequencing (RACE-seq). RACE products from the targeted genes in 23 CRC samples were pooled together and sequenced. We identified VWA2-TCF7L2, DHX35-BPIFA2 and CASZ1-MASP2 as private fusion events, and novel transcript structures for 17 of the 23 other candidate genes. The high-throughput approach facilitated identification of CRC specific RNA variants. These include a recurrent read-through fusion transcript between KLK8 and KLK7, and a splice variant of S100A2. Both of these were overrepresented in CRC tissue and cell lines from external RNA-seq datasets.

ChimPipe: accurate detection of fusion genes and transcription-induced chimeras from RNA-seq data

Background

Chimeric transcripts are commonly defined as transcripts linking two or more different genes in the genome, and can be explained by various biological mechanisms such as genomic rearrangement, read-through or trans-splicing, but also by technical or biological artefacts. Several studies have shown their importance in cancer, cell pluripotency and motility. Many programs have recently been developed to identify chimeras from Illumina RNA-seq data (mostly fusion genes in cancer). However outputs of different programs on the same dataset can be widely inconsistent, and tend to include many false positives. Other issues relate to simulated datasets restricted to fusion genes, real datasets with limited numbers of validated cases, result inconsistencies between simulated and real datasets, and gene rather than junction level assessment.

Results

Here we present ChimPipe, a modular and easy-to-use method to reliably identify fusion genes and transcription-induced chimeras from paired-end Illumina RNA-seq data. We have also produced realistic simulated datasets for three different read lengths, and enhanced two gold-standard cancer datasets by associating exact junction points to validated gene fusions. Benchmarking ChimPipe together with four other state-of-the-art tools on this data showed ChimPipe to be the top program at identifying exact junction coordinates for both kinds of datasets, and the one showing the best trade-off between sensitivity and precision. Applied to 106 ENCODE human RNA-seq datasets, ChimPipe identified 137 high confidence chimeras connecting the protein coding sequence of their parent genes. In subsequent experiments, three out of four predicted chimeras, two of which recurrently expressed in a large majority of the samples, could be validated. Cloning and sequencing of the three cases revealed several new chimeric transcript structures, 3 of which with the potential to encode a chimeric protein for which we hypothesized a new role. Applying ChimPipe to human and mouse ENCODE RNA-seq data led to the identification of 131 recurrent chimeras common to both species, and therefore potentially conserved.

Conclusions

ChimPipe combines discordant paired-end reads and split-reads to detect any kind of chimeras, including those originating from polymerase read-through, and shows an excellent trade-off between sensitivity and precision. The chimeras found by ChimPipe can be validated in-vitro with high accuracy.

Electronic supplementary material

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

Markers of disease and steroid responsiveness in paediatric idiopathic nephrotic syndrome: Whole-transcriptome sequencing of peripheral blood mononuclear cells

Objective

To identify markers of disease and steroid responsiveness in paediatric idiopathic nephrotic syndrome.

Methods

Whole-transcriptome sequencing was performed of peripheral blood mononuclear cells (PBMCs) from patients with NS. Differentially expressed genes (DEGs) were identified in patients with active NS vs those in remission, and those with steroid-sensitive NS (SSNS) vs steroid-resistant NS (SRNS).

Results

A total of 1065 DEGs were identified in patients with NS (n = 10) vs those in remission (n = 9). These DEGs correlated with cytokine and/or immune system signalling and the extracellular matrix. Comparisons between SSNS (n = 6) and SRNS (n = 4) identified 1890 DEGs. These markers of steroid responsiveness were enriched with genes related to the cell cycle, targets of microRNAs, and genes related to cytokines.

Conclusions

Meaningful DEGs were identified. Additional studies with larger numbers of patients will provide more comprehensive data.

ChiTaRS-3.1-the enhanced chimeric transcripts and RNA-seq database matched with protein-protein interactions

Discovery of chimeric RNAs, which are produced by chromosomal translocations as well as the joining of exons from different genes by trans-splicing, has added a new level of complexity to our study and understanding of the transcriptome. The enhanced ChiTaRS-3.1 database (http://chitars.md.biu.ac.il) is designed to make widely accessible a wealth of mined data on chimeric RNAs, with easy-to-use analytical tools built-in. The database comprises 34 922 chimeric transcripts along with 11 714 cancer breakpoints. In this latest version, we have included multiple cross-references to GeneCards, iHop, PubMed, NCBI, Ensembl, OMIM, RefSeq and the Mitelman collection for every entry in the ‘Full Collection’. In addition, for every chimera, we have added a predicted chimeric protein–protein interaction (ChiPPI) network, which allows for easy visualization of protein partners of both parental and fusion proteins for all human chimeras. The database contains a comprehensive annotation for 34 922 chimeric transcripts from eight organisms, and includes the manual annotation of 200 sense-antiSense (SaS) chimeras. The current improvements in the content and functionality to the ChiTaRS database make it a central resource for the study of chimeric transcripts and fusion proteins.

Long non-coding RNAs in human early embryonic development and their potential in ART

BACKGROUND

Human long non-coding RNAs (lncRNAs) are an emerging category of transcripts with increasingly documented functional roles during development. LncRNAs and roles during human early embryo development have recently begun to be unravelled.

OBJECTIVE AND RATIONALE

This review summarizes the most recent knowledge on lncRNAs and focuses on their expression patterns and role during early human embryo development and in pluripotent stem cells (PSCs). Public mRNA sequencing (mRNA-seq) data were used to illustrate these expression signatures.

SEARCH METHODS

The PubMed and EMBASE databases were first interrogated using specific terms, such as 'lncRNAs', to get an extensive overview on lncRNAs up to February 2016, and then using 'human lncRNAs' and 'embryo', 'development', or 'PSCs' to focus on lncRNAs involved in human embryo development or in PSC.Recently published RNA-seq data from human oocytes and pre-implantation embryos (including single-cell data), PSC and a panel of normal and malignant adult tissues were used to describe the specific expression patterns of some lncRNAs in early human embryos.

OUTCOMES

The existence and the crucial role of lncRNAs in many important biological phenomena in each branch of the life tree are now well documented. The number of identified lncRNAs is rapidly increasing and has already outnumbered that of protein-coding genes. Unlike small non-coding RNAs, a variety of mechanisms of action have been proposed for lncRNAs. The functional role of lncRNAs has been demonstrated in many biological and developmental processes, including cell pluripotency induction, X-inactivation or gene imprinting. Analysis of RNA-seq data highlights that lncRNA abundance changes significantly during human early embryonic development. This suggests that lncRNAs could represent candidate biomarkers for developing non-invasive tests for oocyte or embryo quality. Finally, some of these lncRNAs are also expressed in human cancer tissues, suggesting that reactivation of an embryonic lncRNA program may contribute to human malignancies.

WIDER IMPLICATIONS

LncRNAs are emerging potential key players in gene expression regulation. Analysis of RNA-seq data from human pre-implantation embryos identified lncRNA signatures that are specific to this critical step. We anticipate that further studies will show that these new transcripts are major regulators of embryo development. These findings might also be used to develop new tests/treatments for improving the pregnancy success rate in IVF procedures or for regenerative medicine applications involving PSC.

Intergenically Spliced Chimeric RNAs in Cancer

Gene fusions and their encoded products (fusion RNAs and proteins) are viewed as one of the hallmarks of cancer. Traditionally, they were thought to be generated solely by chromosomal rearrangements. However, recent discoveries of trans-splicing and cis-splicing events between neighboring genes, suggest that there are other mechanisms to generate chimeric fusion RNAs without corresponding changes in DNA. In addition, chimeric RNAs have been detected in normal physiology, complicating the use of fusions in cancer detection and therapy. On the other hand, "intergenically spliced" fusion RNAs represent a new repertoire of biomarkers and therapeutic targets. Here, we review current knowledge on chimeric RNAs and implications for cancer detection and treatment, and discuss outstanding questions for the advancement of the field.

Complexity and Specificity of the Neutrophil Transcriptomes in Juvenile Idiopathic Arthritis

NIH projects such as ENCODE and Roadmap Epigenomics have revealed surprising complexity in the transcriptomes of mammalian cells. In this study, we explored transcriptional complexity in human neutrophils, cells generally regarded as nonspecific in their functions and responses. We studied distinct human disease phenotypes and found that, at the gene, gene isoform, and miRNA level, neutrophils exhibit considerable specificity in their transcriptomes. Thus, even cells whose responses are considered non-specific show tailoring of their transcriptional repertoire toward specific physiologic or pathologic contexts. We also found that miRNAs had a global impact on neutrophil transcriptome and are associated with innate immunity in juvenile idiopathic arthritis (JIA). These findings have important implications for our understanding of the link between genes, non-coding transcripts and disease phenotypes.

Discovering and understanding oncogenic gene fusions through data intensive computational approaches

Although gene fusions have been recognized as important drivers of cancer for decades, our understanding of the prevalence and function of gene fusions has been revolutionized by the rise of next-generation sequencing, advances in bioinformatics theory and an increasing capacity for large-scale computational biology. The computational work on gene fusions has been vastly diverse, and the present state of the literature is fragmented. It will be fruitful to merge three camps of gene fusion bioinformatics that appear to rarely cross over: (i) data-intensive computational work characterizing the molecular biology of gene fusions; (ii) development research on fusion detection tools, candidate fusion prioritization algorithms and dedicated fusion databases and (iii) clinical research that seeks to either therapeutically target fusion transcripts and proteins or leverages advances in detection tools to perform large-scale surveys of gene fusion landscapes in specific cancer types. In this review, we unify these different-yet highly complementary and symbiotic-approaches with the view that increased synergy will catalyze advancements in gene fusion identification, characterization and significance evaluation.

Functional annotation of the vlinc class of non-coding RNAs using systems biology approach

Functionality of the non-coding transcripts encoded by the human genome is the coveted goal of the modern genomics research. While commonly relied on the classical methods of forward genetics, integration of different genomics datasets in a global Systems Biology fashion presents a more productive avenue of achieving this very complex aim. Here we report application of a Systems Biology-based approach to dissect functionality of a newly identified vast class of very long intergenic non-coding (vlinc) RNAs. Using highly quantitative FANTOM5 CAGE dataset, we show that these RNAs could be grouped into 1542 novel human genes based on analysis of insulators that we show here indeed function as genomic barrier elements. We show that vlinc RNAs genes likely function in cisto activate nearby genes. This effect while most pronounced in closely spaced vlinc RNA-gene pairs can be detected over relatively large genomic distances. Furthermore, we identified 101 vlinc RNA genes likely involved in early embryogenesis based on patterns of their expression and regulation. We also found another 109 such genes potentially involved in cellular functions also happening at early stages of development such as proliferation, migration and apoptosis. Overall, we show that Systems Biology-based methods have great promise for functional annotation of non-coding RNAs.

Non-coding RNAs: An Introduction

For many years the main role of RNA, it addition to the housekeeping functions of for example tRNAs and rRNAs, was believed to be a messenger between the genes encoded on the DNA and the functional units of the cell, the proteins. This changed drastically with the identification of the first small non-coding RNA, termed microRNA, some 20 years ago. This discovery opened the field of regulatory RNAs with no or little protein-coding potential. Since then many new classes of regulatory non-coding RNAs, including endogenous small interfering RNAs (endo-siRNAs), PIWI-associated RNAs (piRNAs), and long non-coding RNAs, have been identified and we have made amazing progress in elucidating their expression, biogenesis, mechanisms and mode of action, and function in many, if not all, biological processes. In this chapter we provide an introduction about the current knowledge of the main classes of non-coding RNAs, what is know about their biogenesis and mechanism of function.

Abundant and broad expression of transcription-induced chimeras and protein products in mammalian genomes

The expression of transcription-induced chimeras (TICs) was underestimated due to strategic and logical reasons. In order to thoroughly examine TICs, systematic survey of TIC events was conducted in mammalian genomes using ESTs, followed by experimental validation using RT-PCR and real-time quantitative PCR (qPCR). The expression of ∼98% TIC events in at least one tissue or cell line from both mouse and human was verified. Besides, ∼40% TICs were broadly expressed, and ∼33% of TICs showed expression levels comparable to or higher than their upstream parental genes. We further identified putative chimeric proteins in public databases and validated two using Western blotting. GO analysis revealed that proteins resided in one multi-protein complex or functioning in metabolic or signaling pathway tended to produce fused products. Taken together, we have shown substantial evidence for the underestimated TIC events; and TICs could be a novel regulated way to further increases the proteome complexity in mammalian genomes. Possible regulation mechanisms and evolution of TICs were also discussed.

Identification and analysis of the promoter region of the STGC3 gene

INTRODUCTION

Nasopharyngeal carcinoma (NPC) is a common malignant tumor of the head and neck. The STGC3 gene is related to development of nasopharyngeal cancer. The aim of this study is to explore the promoter region of the STGC3 gene.

MATERIAL AND METHODS

The bioinformatic technique was applied to predict its promoter region and construct the gene promoter region luciferase for the gene vector and transfection of the human embryonic kidney epithelial 293T cell line, human nasopharyngeal carcinoma CNE2 cell line and immortalized nasopharyngeal epithelial NP69 cell line. The recombinant plasmid pGL3-en283, pGL3-en281, pGL3-en571, empty plasmid pGL3-control, negative control pGL3-enhance and internal control of marine intestine luciferase expression vector pRL-SV40 were transfected into NP69 cells, 293T cells and CNE2 cells. Dual luciferase activity detection showed luciferase luminescence values and marine intestine luciferase luminescence values. Relative luciferase activity (RLA) in each cell was calculated.

RESULTS

We observed strong promoter activity of plasmid pGL3-en283, pGL3-en281 and pGL3-en571 in NP69, 293T and CNE2 cells compared with the negative control pGL3-enhance plasmid. Among them, pGL3-en281 showed the strongest promoter activity, and these three kinds of recombinant plasmids showed stronger promoter activity in 293T cells than in CNE2 cells.

CONCLUSIONS

The pGL3-en281 plasmid showed stronger promoter activity than pGL3-en571 in the three cells, indicating that -11048 bp to -653 bp might be the core promoter region.

Long non-coding RNA SOX2OT: expression signature, splicing patterns, and emerging roles in pluripotency and tumorigenesis

SOX2 overlapping transcript (SOX2OT) is a long non-coding RNA which harbors one of the major regulators of pluripotency, SOX2 gene, in its intronic region. SOX2OT gene is mapped to human chromosome 3q26.3 (Chr3q26.3) locus and is extended in a high conserved region of over 700 kb. Little is known about the exact role of SOX2OT; however, recent studies have demonstrated a positive role for it in transcription regulation of SOX2 gene. Similar to SOX2, SOX2OT is highly expressed in embryonic stem cells and down-regulated upon the induction of differentiation. SOX2OT is dynamically regulated during the embryogenesis of vertebrates, and delimited to the brain in adult mice and human. Recently, the disregulation of SOX2OT expression and its concomitant expression with SOX2 have become highlighted in some somatic cancers including esophageal squamous cell carcinoma, lung squamous cell carcinoma, and breast cancer. Interestingly, SOX2OT is differentially spliced into multiple mRNA-like transcripts in stem and cancer cells. In this review, we are describing the structural and functional features of SOX2OT, with an emphasis on its expression signature, its splicing patterns and its critical function in the regulation of SOX2 expression during development and tumorigenesis.

JAFFA: High sensitivity transcriptome-focused fusion gene detection

Genomic instability is a hallmark of cancer and, as such, structural alterations and fusion genes are common events in the cancer landscape. RNA sequencing (RNA-Seq) is a powerful method for profiling cancers, but current methods for identifying fusion genes are optimised for short reads. JAFFA (https://github.com/Oshlack/JAFFA/wiki) is a sensitive fusion detection method that outperforms other methods with reads of 100 bp or greater. JAFFA compares a cancer transcriptome to the reference transcriptome, rather than the genome, where the cancer transcriptome is inferred using long reads directly or by de novo assembling short reads.

The Landscape of long noncoding RNA classification

Advances in the depth and quality of transcriptome sequencing have revealed many new classes of long noncoding RNAs (lncRNAs). lncRNA classification has mushroomed to accommodate these new findings, even though the real dimensions and complexity of the noncoding transcriptome remain unknown. Although evidence of functionality of specific lncRNAs continues to accumulate, conflicting, confusing, and overlapping terminology has fostered ambiguity and lack of clarity in the field in general. The lack of fundamental conceptual unambiguous classification framework results in a number of challenges in the annotation and interpretation of noncoding transcriptome data. It also might undermine integration of the new genomic methods and datasets in an effort to unravel the function of lncRNA. Here, we review existing lncRNA classifications, nomenclature, and terminology. Then, we describe the conceptual guidelines that have emerged for their classification and functional annotation based on expanding and more comprehensive use of large systems biology-based datasets.

Computational characterisation of cancer molecular profiles derived using next generation sequencing

Our current understanding of cancer genetics is grounded on the principle that cancer arises from a clone that has accumulated the requisite somatically acquired genetic aberrations, leading to the malignant transformation. It also results in aberrent of gene and protein expression. Next generation sequencing (NGS) or deep sequencing platforms are being used to create large catalogues of changes in copy numbers, mutations, structural variations, gene fusions, gene expression, and other types of information for cancer patients. However, inferring different types of biological changes from raw reads generated using the sequencing experiments is algorithmically and computationally challenging. In this article, we outline common steps for the quality control and processing of NGS data. We highlight the importance of accurate and application-specific alignment of these reads and the methodological steps and challenges in obtaining different types of information. We comment on the importance of integrating these data and building infrastructure to analyse it. We also provide exhaustive lists of available software to obtain information and point the readers to articles comparing software for deeper insight in specialised areas. We hope that the article will guide readers in choosing the right tools for analysing oncogenomic datasets.

ChiTaRS 2.1--an improved database of the chimeric transcripts and RNA-seq data with novel sense-antisense chimeric RNA transcripts

Chimeric RNAs that comprise two or more different transcripts have been identified in many cancers and among the Expressed Sequence Tags (ESTs) isolated from different organisms; they might represent functional proteins and produce different disease phenotypes. The ChiTaRS 2.1 database of chimeric transcripts and RNA-Seq data (http://chitars.bioinfo.cnio.es/) is the second version of the ChiTaRS database and includes improvements in content and functionality. Chimeras from eight organisms have been collated including novel sense–antisense (SAS) chimeras resulting from the slippage of the sense and anti-sense intragenic regions. The new database version collects more than 29 000 chimeric transcripts and indicates the expression and tissue specificity for 333 entries confirmed by RNA-seq reads mapping the chimeric junction sites. User interface allows for rapid and easy analysis of evolutionary conservation of fusions, literature references and experimental data supporting fusions in different organisms. More than 1428 cancer breakpoints have been automatically collected from public databases and manually verified to identify their correct cross-references, genomic sequences and junction sites. As a result, the ChiTaRS 2.1 collection of chimeras from eight organisms and human cancer breakpoints extends our understanding of the evolution of chimeric transcripts in eukaryotes as well as their functional role in carcinogenic processes.

De-repressing LncRNA-Targeted Genes to Upregulate Gene Expression: Focus on Small Molecule Therapeutics

Non-protein coding RNAs (ncRNAs) make up the overwhelming majority of transcripts in the genome and have recently gained attention for their complex regulatory role in cells, including the regulation of protein-coding genes. Furthermore, ncRNAs play an important role in normal development and their expression levels are dysregulated in several diseases. Recently, several long noncoding RNAs (lncRNAs) have been shown to alter the epigenetic status of genomic loci and suppress the expression of target genes. This review will present examples of such a mechanism and focus on the potential to target lncRNAs for achieving therapeutic gene upregulation by de-repressing genes that are epigenetically silenced in various diseases. Finally, the potential to target lncRNAs, through their interactions with epigenetic enzymes, using various tools, such as small molecules, viral vectors and antisense oligonucleotides, will be discussed. We suggest that small molecule modulators of a novel class of drug targets, lncRNA-protein interactions, have great potential to treat some cancers, cardiovascular disease, and neurological disorders.

DDX3X, the X homologue of AZFa gene DDX3Y, expresses a complex pattern of transcript variants only in the male germ line

DDX3X, the functional X homologue of the major AZFa gene, DDX3Y, belongs to the highly conserved PL10-subfamily of DEAD-box RNA helicase genes which are functionally conserved from yeast to man. They are mainly involved in cell cycle control and translation initiation control of gene transcripts with long 5'UTR extensions containing complex secondary structures. Interestingly, in humans both gene copies were found to be expressed at different phases of human spermatogenesis. Whereas DDX3Y transcripts are translated only in premeiotic male germ cells, the DDX3X protein is expressed only in postmeiotic spermatids. In this study, we found that the major class of DDX3X transcripts in human testis become activated first after meiosis and at a specific core promoter not active in somatic tissues and not present upstream of the DDX3Y homologue. Two alternative 5'UTR transcript lengths are subsequently produced by an additional testis-specific 5'UTR splicing event. Both transcripts are mainly processed for polyadenylation in their proximal 3'UTR. A minor transcript class starting at the same male germ line-specific core promoter produces primary transcripts with an extremely long 3'UTR (∼ 17 kb), which is subsequently spliced at distinct sites resulting in six short 3'UTR splice variants (I-VI). Comparative analyses of the DDX3X transcripts in mouse and primates revealed that this complex pattern of male germ line-specific transcript variants first evolved in primates. Our data thus suggest complex translational control mechanism(s) for the human DDX3X gene locus functioning only in the male germ line and resulting in expression of its protein only in the postmeiotic spermatids.

Tandem RNA chimeras contribute to transcriptome diversity in human population and are associated with intronic genetic variants

Chimeric RNAs originating from two or more different genes are known to exist not only in cancer, but also in normal tissues, where they can play a role in human evolution. However, the exact mechanism of their formation is unknown. Here, we use RNA sequencing data from 462 healthy individuals representing 5 human populations to systematically identify and in depth characterize 81 RNA tandem chimeric transcripts, 13 of which are novel. We observe that 6 out of these 81 chimeras have been regarded as cancer-specific. Moreover, we show that a prevalence of long introns at the fusion breakpoint is associated with the chimeric transcripts formation. We also find that tandem RNA chimeras have lower abundances as compared to their partner genes. Finally, by combining our results with genomic data from the same individuals we uncover intronic genetic variants associated with the chimeric RNA formation. Taken together our findings provide an important insight into the chimeric transcripts formation and open new avenues of research into the role of intronic genetic variants in post-transcriptional processing events.

Chimeric RNAs generated by intergenic splicing in normal and cancer cells

A hallmark of many neoplasias is chromosomal rearrangement, an event that commonly results in the fusion of two separate genes. The RNA and protein resulting from these gene fusions often play critical roles in cancer development, maintenance, and progression. Traditionally, these fusion products are thought to be produced solely due to DNA level changes and are therefore considered unique to cancer. Recent advances in microarray and deep-sequencing have revealed many more fusion transcripts. Surprisingly, some are without detectable rearrangement at the DNA level. Reports have demonstrated that at least some of these chimeric RNAs are generated via intergenic splicing. In this review, we highlight three examples of these noncanonical chimeric transcripts that are formed by trans-splicing or cis-splicing of adjacent genes and summarize the knowledge we have regarding these noncanonical fusions. We discuss the implications of the chimeric RNAs in both cancer and normal physiology, as some of these fusion transcripts are found in normal, noncancerous cells with sequences identical to those generated by canonical chromosomal translocation found in cancer cells. Finally, we present methods that are currently being used to discover additional chimeric RNAs.

Retrotransposons in pluripotent cells: Impact and new roles in cellular plasticity

Transposable Elements are pieces of DNA able to mobilize from one location to another within genomes. Although they constitute more than 50% of the human genome, they have been classified as selfish DNA, with the only mission to spread within genomes and generate more copies of themselves that will ensure their presence over generations. Despite their remarkable prevalence, only a minor group of transposable elements remain active in the human genome and can sporadically be associated with the generation of a genetic disorder due to their ongoing mobility. Most of the transposable elements identified in the human genome corresponded to fixed insertions that no longer move in genomes. As selfish DNA, transposable element insertions accumulate in cell types where genetic information can be passed to the next generation. Indeed, work from different laboratories has demonstrated that the main heritable load of TE accumulation in humans occurs during early embryogenesis. Thus, active transposable elements have a clear impact on our pluripotent genome. However, recent findings suggest that the main proportion of fixed non-mobile transposable elements might also have emerging roles in cellular plasticity. In this concise review, we provide an overview of the impact of currently active transposable elements in our pluripotent genome and further discuss new roles of transposable elements (active or not) in regulating pluripotency. This article is part of a Special Issue entitled: Stress as a fundamental theme in cell plasticity.

Dark matter RNA illuminates the puzzle of genome-wide association studies

In the past decade, numerous studies have made connections between sequence variants in human genomes and predisposition to complex diseases. However, most of these variants lie outside of the charted regions of the human genome whose function we understand; that is, the sequences that encode proteins. Consequently, the general concept of a mechanism that translates these variants into predisposition to diseases has been lacking, potentially calling into question the validity of these studies. Here we make a connection between the growing class of apparently functional RNAs that do not encode proteins and whose function we do not yet understand (the so-called 'dark matter' RNAs) and the disease-associated variants. We review advances made in a different genomic mapping effort - unbiased profiling of all RNA transcribed from the human genome - and provide arguments that the disease-associated variants exert their effects via perturbation of regulatory properties of non-coding RNAs existing in mammalian cells.

Epigenetic reprogramming in Mist1(-/-) mice predicts the molecular response to cerulein-induced pancreatitis

Gene expression is affected by modifications to histone core proteins within chromatin. Changes in these modifications, or epigenetic reprogramming, can dictate cell fate and promote susceptibility to disease. The goal of this study was to determine the extent of epigenetic reprogramming in response to chronic stress that occurs following ablation of MIST1 (Mist1^−/−), which is repressed in pancreatic disease. Chromatin immunoprecipitation for trimethylation of lysine residue 4 on histone 3 (H3K4Me3) in purified acinar cells from wild type and Mist1^−/− mice was followed by Next Generation sequencing (ChIP-seq) or ChIP-qPCR. H3K4Me3-enriched genes were assessed for expression by qRT-PCR in pancreatic tissue before and after induction of cerulein-induced pancreatitis. While most of H3K4Me3-enrichment is restricted to transcriptional start sites, >25% of enrichment sites are found within, downstream or between annotated genes. Less than 10% of these sites were altered in Mist1^−/− acini, with most changes in H3K4Me3 enrichment not reflecting altered gene expression. Ingenuity Pathway Analysis of genes differentially-enriched for H3K4Me3 revealed an association with pancreatitis and pancreatic ductal adenocarcinoma in Mist1^−/− tissue. Most of these genes were not differentially expressed but several were readily induced by acute experimental pancreatitis, with significantly increased expression in Mist1^−/− tissue relative to wild type mice. We suggest that the chronic cell stress observed in the absence of MIST1 results in epigenetic reprogramming of genes involved in promoting pancreatitis to a poised state, thereby increasing the sensitivity to events that promote disease.

Comprehensive analysis of the Corynebacterium glutamicum transcriptome using an improved RNAseq technique

Background

The use of RNAseq to resolve the transcriptional organization of an organism was established in recent years and also showed the complexity and dynamics of bacterial transcriptomes. The aim of this study was to comprehensively investigate the transcriptome of the industrially relevant amino acid producer and model organism Corynebacterium glutamicum by RNAseq in order to improve its genome annotation and to describe important features for transcription and translation.

Results

RNAseq data sets were obtained by two methods, one that focuses on 5′-ends of primary transcripts and another that provides the overall transcriptome with an improved resolution of 3′-ends of transcripts. Subsequent data analysis led to the identification of more than 2,000 transcription start sites (TSSs), the definition of 5′-UTRs (untranslated regions) for annotated protein-coding genes, operon structures and many novel transcripts located between or in antisense orientation to protein-coding regions. Interestingly, a high number of mRNAs (33%) is transcribed as leaderless transcripts. From the data, consensus promoter and ribosome binding site (RBS) motifs were identified and it was shown that the majority of genes in C. glutamicum are transcribed monocistronically, but operons containing up to 16 genes are also present.

Conclusions

The comprehensive transcriptome map of C. glutamicum established in this study represents a major step forward towards a complete definition of genetic elements (e.g. promoter regions, gene starts and stops, 5′-UTRs, RBSs, transcript starts and ends) and provides the ideal basis for further analyses on transcriptional regulatory networks in this organism. The methods developed are easily applicable for other bacteria and have the potential to be used also for quantification of transcriptomes, replacing microarrays in the near future.

Functional insights into the role of nuclear-retained long noncoding RNAs in gene expression control in mammalian cells

The mammalian genome harbors thousands of long noncoding RNA (lncRNA) genes. Recent studies have indicated the involvement of several of these lncRNAs in the regulation of gene expression. lncRNAs play crucial roles in various biological processes ranging from epigenetic gene regulation, transcriptional control,to post-transcriptional regulation. lncRNAs are localized in various subcellular compartments, and major proportion of these are retained in the cell nucleus and could be broadly classified as nuclear-retained lncRNAs (nrRNAs). Based on the identified functions,members of the nrRNAs execute diverse roles, including providing architectural support to the hierarchical subnuclear organization and influencing the recruitment of chromatin modifier factors to specific chromatin sites. In this review, we will summarize the recently described roles of mammalian nrRNAs in controlling gene expression by influencing chromatin organization, transcription,pre-mRNA processing, nuclear organization, and their involvement in disease.

Systematic classification of non-coding RNAs by epigenomic similarity

Background

Even though only 1.5% of the human genome is translated into proteins, recent reports indicate that most of it is transcribed into non-coding RNAs (ncRNAs), which are becoming the subject of increased scientific interest. We hypothesized that examining how different classes of ncRNAs co-localized with annotated epigenomic elements could help understand the functions, regulatory mechanisms, and relationships among ncRNA families.

Results

We examined 15 different ncRNA classes for statistically significant genomic co-localizations with cell type-specific chromatin segmentation states, transcription factor binding sites (TFBSs), and histone modification marks using GenomeRunner (http://www.genomerunner.org). P-values were obtained using a Chi-square test and corrected for multiple testing using the Benjamini-Hochberg procedure. We clustered and visualized the ncRNA classes by the strength of their statistical enrichments and depletions.

We found piwi-interacting RNAs (piRNAs) to be depleted in regions containing activating histone modification marks, such as H3K4 mono-, di- and trimethylation, H3K27 acetylation, as well as certain TFBSs. piRNAs were further depleted in active promoters, weak transcription, and transcription elongation regions, and enriched in repressed and heterochromatic regions. Conversely, transfer RNAs (tRNAs) were depleted in heterochromatin regions and strongly enriched in regions containing activating H3K4 di- and trimethylation marks, H2az histone variant, and a variety of TFBSs. Interestingly, regions containing CTCF insulator protein binding sites were associated with tRNAs. tRNAs were also enriched in the active, weak and poised promoters and, surprisingly, in regions with repetitive/copy number variations.

Conclusions

Searching for statistically significant associations between ncRNA classes and epigenomic elements permits detection of potential functional and/or regulatory relationships among ncRNA classes, and suggests cell type-specific biological roles of ncRNAs.

Intron retention and transcript chimerism conserved across mammals: Ly6g5b and Csnk2b-Ly6g5b as examples

Background

Alternative splicing (AS) is a major mechanism for modulating gene expression of an organism, allowing the synthesis of several structurally and functionally distinct mRNAs and protein isoforms from a unique gene. Related to AS is the Transcription Induced Chimerism (TIC) or Tandem Chimerism, by which chimeric RNAs between adjacent genes can be found, increasing combinatorial complexity of the proteome. The Ly6g5b gene presents particular behaviours in its expression, involving an intron retention event and being capable to form RNA chimera transcripts with the upstream gene Csnk2b. We wanted to characterise these events more deeply in four tissues in six different mammals and analyse their protein products.

Results

While canonical Csnk2b isoform was widely expressed, Ly6g5b canonical isoform was less ubiquitous, although the Ly6g5b first intron retained transcript was present in all the tissues and species analysed. Csnk2b-Ly6g5b chimeras were present in all the samples analysed, but with restricted expression patterns. Some of these chimeric transcripts maintained correct structural domains from Csnk2b and Ly6g5b. Moreover, we found Csnk2b, Ly6g5b, and Csnk2b-Ly6g5b transcripts that present exon skipping, alternative 5' and 3' splice site and intron retention events. These would generate truncated or aberrant proteins whose role remains unknown. Some chimeric transcripts would encode CSNK2B proteins with an altered C-terminus, which could affect its biological function broadening its substrate specificity. Over-expression of human CSNK2B, LY6G5B, and CSNK2B-LY6G5B proteins, show different patterns of post-translational modifications and cell distribution.

Conclusions

Ly6g5b intron retention and Csnk2b-Ly6g5b transcript chimerism are broadly distributed in tissues of different mammals.