- Research Article
- Open Access
Evolution of the unspliced transcriptome
© Engelhardt and Stadler. 2015
- Received: 16 March 2015
- Accepted: 29 July 2015
- Published: 20 August 2015
Despite their abundance, unspliced EST data have received little attention as a source of information on non-coding RNAs. Very little is know, therefore, about the genomic distribution of unspliced non-coding transcripts and their relationship with the much better studied regularly spliced products. In particular, their evolution has remained virtually unstudied.
We systematically study the evidence on unspliced transcripts available in EST annotation tracks for human and mouse, comprising 104,980 and 66,109 unspliced EST clusters, respectively. Roughly one third of these are located totally inside introns of known genes (TINs) and another third overlaps exonic regions (PINs). Eleven percent are “intergenic”, far away from any annotated gene. Direct evidence for the independent transcription of many PINs and TINs is obtained from CAGE tag and chromatin data. We predict more than 2000 3’UTR-associated RNA candidates for each human and mouse. Fifteen to twenty percent of the unspliced EST cluster are conserved between human and mouse. With the exception of TINs, the sequences of unspliced EST clusters evolve significantly slower than genomic background. Furthermore, like spliced lincRNAs, they show highly tissue-specific expression patterns.
Unspliced long non-coding RNAs are an important, rapidly evolving, component of mammalian transcriptomes. Their analysis is complicated by their preferential association with complex transcribed loci that usually also harbor a plethora of spliced transcripts. Unspliced EST data, although typically disregarded in transcriptome analysis, can be used to gain insights into this rarely investigated transcriptome component. The frequently postulated connection between lack of splicing and nuclear retention and the surprising overlap of chromatin-associated transcripts suggests that this class of transcripts might be involved in chromatin organization and possibly other mechanisms of epigenetic control.
- Transcription Start Site
- UCSC Genome Browser
- RefSeq Gene
- phastCons Score
- Unspliced Transcript
Non-coding RNA (ncRNA) constitutes a significant portion of the mammalian transcriptome [1–4]. Although a large sub-class of long ncRNAs (lncRNAs) is spliced, capped, and polyadenylated, and thus differs from their protein-coding siblings only in coding capacity , these mRNA-like lncRNAs account for only a small fraction of the striking diversity of transcripts. Nuclear retained ncRNAs are often spliced transcripts but not polyadenylated. These “dark matter RNAs”, which have remained largely un-annotated so far, can in fact be the dominating non-ribosomal RNA component in a mammalian cell [6, 7].
As a class, lncRNAs are under purifying selection [8–10] although the level of sequence conservation is typically very low [9, 11]. Comparative transcriptomics [12, 13] as well as computational studies [14, 15] showed that between one and two thirds of human lncRNAs are conserved among Eutheria, emphasizing the functional importance of these transcripts. The catalogs of RNAs used in these studies are, however, heavily biased towards spliced RNAs, and some of them use the conservation of splice sites explicitly as a means to assess conservation [12, 15–17]. Indeed, the overwhelming majority of lncRNAs for which detailed functional information is available is spliced, see e.g. , although splicing often tends to occur only after transcription and is less efficient .
By far the largest class of unspliced transcripts for which a function is known are the intron-less protein-coding genes. In human, they account for about 4.5 % of the protein-coding loci . They have received comparably little attention even though they have distinctive features e.g. related to their export pathways . On average they are expressed at lower levels, tend to be more tissue specific, evolve at faster rates, and are of relatively recent origins . The extremely well-conserved replication-dependent histone genes form a distinctive class of intron-less genes set apart by their unique 3’ end processing .
Unspliced lncRNAs fall into at least four distinct classes: (i) intronic transcripts typically associated with protein-coding loci, (ii) lncRNAs associated with long 3’-UTRs, (iii) independent unspliced RNAs found in intergenic regions, and (iv) an enigmatic class of very long macroRNAs.
Totally and partially intronic transcripts (TINs and PINs) that are usually unspliced and lack coding capacity have been reported in large numbers for both human and mouse [24–26]. This class includes many unspliced long anti-sense intronic RNAs [27, 28]. Not much is known on the biogenesis of intronic lncRNAs, although there are presumably multiple pathways. The anti-sense TIN ANRASSF1 is a pol-II transcript, capped, and polyadenylated . Sense TINs, on the other hand, might also be processing products generated from introns as described for the sisRNAs in frog oocytes . Intronic transcripts appear to have primarily regulatory functions, for which they employ several different molecular mechanisms [29, 31, 32]. LncRNAs can affect gene expression both in cis and in trans by modulating the chromatin structure . Although the first reports that RNA is an integral component of chromatin date back to the 1970s  there has not been much work on chromatin-associated RNAs (CARs). They are largely non-polyadenylated  and derive from both intronic and intergenic regions . Some of them are probably natural antisense transcripts (NATs) deriving from both coding and non-coding loci that act in cis as epigenetic regulators of gene expression and chromatin remodeling . Furthermore RNP complexes stably tethered to chromatin were characterized recently that appear to play a role in regulating higher order structures of chromatin . Many promoter-associated long ncRNAs  might also function via chromatin modification .
The 3’ untranslated regions (3’UTRs) of eukaryotic genes are well known to harbor sequence and secondary structure elements that regulate mRNA stability, localization, and translation [41–43]. Many 3’UTRs in animals can also be decoupled from the protein-coding part forming so-called UTR-associated RNAs (uaRNAs), which appear to function as non-coding RNAs in trans . Note that the abbreviation uaRNA was later also used for the unrelated class of “upstream antisense RNAs” in . In the form of independent uaRNAs, they are often detectable as unspliced ESTs. A recent study, furthermore, found parallels in sequence composition between lncRNAs and 3’ UTRs that sets both groups apart from 5’UTRs and coding regions . Since 3’UTRs typically harbor microRNA target sites they may function as sponges for microRNAs [47, 48] or RNA binding proteins .
The best-known examples of “intergenic” unspliced ncRNAs are MALAT-1 and MEN β, which organize nuclear structures known as SC35 speckles and paraspeckles, respectively [50–52]. Both transcripts are spliced only infrequently  and are rather well-conserved . This class also contains important disease-associated RNAs such as PRNCR1 .
MacroRNAs covering up to several hundred kb were recently observed as highly expressed RNAs in signaling pathways  and in cancer cells . Very similar transcripts such as Airn [57, 58] and KCNQ1OT1  are involved in the regulation of imprinted loci. These enigmatic RNAs are very poorly conserved.
Despite the wealth of unspliced EST data that is available in public databases, they have received little attention as a source of information on non-coding RNAs apart from a seminal work . In part this is probably due the difficulties involved in mining this resource. First, strand information is usually not available and the absence of splice junctions makes it impossible to guess the reading direction. Furthermore, a detailed analysis of nearly 40,000 putative ncRNAs from RIKEN’s FANTOM transcript data set showed that many of the putative intronic and intergenic transcripts might be artifacts and in fact are internally primed subsequences from even longer transcripts . The fragmentary nature of ESTs, furthermore, makes it difficult and often impossible to reliably determine transcript boundaries.
In this contribution we use unspliced EST data to obtain an overview of the unspliced transcripts in the human and the mouse genome. We focus on ESTs rather than NGS data here primarily because the much longer EST sequences are largely mapped unambiguously to the genome and are much more likely to contain splice junctions when they originate from spliced transcripts than the much shorted NGS reads. We update and extend previous investigations in the human uEST data [24, 61] and in particular provide a first overview of the similarities and differences of the situation in human and mouse.
ESTs “within range” of RefSeq genes
The majority of unspliced EST (uEST) clusters, 89 % in both human and mouse, overlaps with or lies in the vicinity of annotated RefSeq genes. A possible reason for this strong association with known pre-existing annotation could be that uEST cluster are just by-products of “normal” spliced transcripts arising from occasionally inefficient splicing or a background of not yet processed primary transcripts. We therefore compare the amount of spliced and unspliced ESTs with the range of annotated RefSeq genes (including a 5 kb flanking region), see Additional file 1: Figure S1 for human and Additional file 1: Figure S2 for mouse. The amount of ESTs is a measure that interpolates between diversity of tissues and conditions under which the locus is expressed (for low counts) and a genuine proxy for expression levels (for large counts). Despite its aggregate nature, the correlation of the amount of spliced and unspliced ESTs is indicative of the overall coupling between spliced and unspliced expression. We observe a Pearson correlation coefficient of r=0.75 and r=0.5 for spliced and unspliced RefSeq genes in human and r=0.72 and r=0.48 in mouse gene, respectively. The correlation coefficients are rather moderate, however, reaching roughly the level of correlation between mRNA and protein abundances, see e.g. [62, 63]. We take this as a strong indication that the observed unspliced EST clusters are not just a byproduct. The quantitative data suggest that a considerable fraction of unspliced transcripts is independent of spliced genes.
Human and mouse unspliced EST clusters overlapping a predicted TSS peak in their putative 5’ region and their orientation relative to the surrounding RefSeq gene. RefSeq TSS are all unspliced EST cluster that overlap any 5’UTR or that are located completely within the upstream region of a RefSeq gene. TINs and PINs are defined as in Fig. 1. All unspliced EST cluster that have an overlap with the 3’UTR but not with the 5’ UTR are interpreted as uaRNAs. “Sense” and “Antisense” are relative to the reading direction of the RefSeq gene. “Ambiguous” clusters could not be assigned to an orientation because of conflicting information on their directions. The column “All” lists the number of corresponding unspliced EST clusters without considering an overlap with CAGE-tags
The protein-coding part of RefSeq genes overlaps 39,791 unspliced EST cluster (27,897 in mouse). Less than 4 % of these cluster in each species, however, are located completely in the coding region. We analyze not only uESTs overlapping the gene body but take into account the vicinity of RefSeq genes, defined as within 5 kb of the ends of the RefSeq entry. There are 2459 (1601) unspliced EST cluster in upstream region of human (mouse) RefSeq genes and 3950 (2600) downstream, see Fig. 1. Conversely, 37,561 of 42,165 (89 %) of the annotated RefSeq genes (ignoring separately annotated haplotypes) overlap with at least one unspliced EST cluster.
Independent UTR-derived RNAs
Unspliced EST (uEST) cluster are predominately expressed in the 3’ region of RefSeq genes. There are more than 6000 human and 5000 mouse uEST clusters that overlap the 3’UTR or are located within 5000 nt downstream of 3’end of a RefSeq genes. The frequent observation of these downstream signals is consistent with the observation that the length of UTRs is often underestimated in current gene annotations . Alternative polyadenylation, and hence often dramatic changes in 3’UTR length, furthermore appear to be frequent phenomenon that may have biased the current gene annotations [65–67]. Unfortunately, we can not rule out that unspliced EST clusters of the class 3’R are only artifacts of the read-through of transcripts from the upstream gene. Therefore, we attempt to use CAGE data to gain additional evidence of the 3’R cluster to be independently transcribed or processed, for instance similar to uaRNAs.
CAGE tags  are routinely used to assay CAPed 5’ ends. The FANTOM 5 consortium predicted TSS peaks by clustering CAGE tags . These peaks can be used to predict the reading direction of an unspliced EST cluster if they are located at one of the ends of the cluster. 30 % (human) and 25 % (mouse) of the predicted TSS peaks (319,019 of 1,048,124 in human and 169,144 of 652,833 in mouse) overlap with a total of 48,435 (human) and 30,056 (mouse) uEST clusters.We are primarily interested in the TSS peaks from the forward strand that are located close to the beginning of the EST cluster, and in the TSS peaks from the reverse strand that are located close to the end of the EST cluster. The reading direction of an EST cluster can be determined unambiguously provided it contains exactly one of these two types of TSS peaks. As in previous work , we use 84 nt from either end of the EST cluster as distance cutoff. We obtained 9378 forward and 8837 reverse strand uEST cluster in human and 6732 forward and 6630 reverse in mouse. For 918 cluster in human and 719 mouse we found TSS peaks at both ends. These cases are likely bi-directionally transcribed loci. We nevertheless excluded them from further analysis, summarized in Additional file 1: Table S3.
Most of them correspond to the 5’ end or to 5’ extensions of RefSeq genes. We excluded them in the further analysis since we can not distinguish the transcription start site (TSS) of the corresponding RefSeq gene and an independent one of an unspliced EST cluster, see Table 1.
The majority (6957/9877, 70 % in human; 3841/5593, 69 % in mouse) of the remaining clusters, share the reading direction of the corresponding gene. Some of these unspliced cluster may correspond to a class of transcribed promoter-associated ncRNAs similar to the CCND1-pncRNAs , which negatively regulate CCND1 transcription by recruiting TLS to the promoter.
Enhancer RNAs are a class of non-coding RNAs that has recently gained considerable interest [70–73]. They are expressed from enhancer regions. Since these are not annotated as (RefSeq) genes, we analyzed whether intergenic unspliced EST cluster might correspond to enhancer RNAs. To this end we used a set of 38,554 predicted human enhancer RNAs compiled in  as part of the FANTOM 5 project. However, only 281 of 11,207 human IGR uEST cluster overlap with one of these enhancer RNAs.
The CAGE data collected by the FANTOM consortium can even be used to analyze the cell-specific expression of unspliced EST (uEST) cluster. Ohmiya et al.  used the FANTOM 5 dataset to predict transcriptions start sites (TSSs) in 156 different human primary cells. Overlapping these TSS with the uEST clusters provides us with an overview of cell-specific expression patterns of many uEST cluster.
Amount of unspliced EST clusters that are expressed in a certain number (1st column) of cell lines. The number in brackets is the relative amount in respect to all unspliced EST clusters of that group for which a predicted TSS was present. Percentages below 5 % are omitted. “Genic” uEST cluster are overlapping or in close vicinity of RefSeq genes
27 (5 %)
479 (6 %)
26 (5 %)
726 (6 %)
897 (11 %)
65 (13 %)
1337 (10 %)
2236 (27 %)
158 (32 %)
3319 (26 %)
We performed the same analysis for all RefSeq genes as well as long non-coding RNAs (lncRNAs) from GENCODE v14 filtered for high reliability, see Additional file 1. The distribution of RefSeq genes is similar to the unspliced EST cluster, although more uEST clusters are expressed in very few tissues and there are fewer nearly ubiquitously expressed uEST clusters compared to protein coding genes. The lncRNAs lack the peak for ubiquitous expression completely, consistent with previous observations that lncRNAs are expressed less abundantly and more specifically compared to traditional transcripts .
We tried to use the cell-specific CAGE data to find cell types with an enrichment for uaRNAs. Using all 156 cell types, normalized for the amount of predicted TSS per cell type, there is, however, no enrichment detectable for any particular cell type. We then combined related cell types into larger groups but, again, no group-wise enrichment’s were visible.
CAGE is a method developed for exactly determining the 5’-end of RNAs. A more indirect way to detect the transcription start site (TSS) of a gene is analyzing the chromatin structure.
We used a genome-segmentation track based on ENCODE data provided by ref.  based on a computational methods for unsupervised chromatin state annotation from ChIP-seq data for multiple histone modifications, general transcription factors, and chromatin accessibility assays. While the chromatin elements were computed independently for six different cell lines in  we combined all cell lines here since the definition of the EST clusters is also a composite of many different cell and tissue types.
More than 80 % of the clusters close to or within annotated genes but outside of exons (i.e., TIN, UT, DT) overlap a segment predicted to be transcriptionally active (‘TSS’ or ‘Transcribed’) based on chromatin state in at least one cell line. Although we know, of course, that unspliced EST clusters are transcribed regions, it is reassuring that the chromatin data show a consistent picture. In particular for UT and DT uEST clusters it is impossible to discern in most cases whether they are independent transcripts or processing products, or whether they belong to large, splice transcripts whose UTRs are annotated only incompletely.
Unspliced EST (uEST) cluster overlapping with chromatin elements predicted from ENCODE data . ‘TSS’ represents “Predicted promoter region including TSS”. ‘Transcribed’ refers to “Predicted transcribed region”. ‘All’ is the number of all uEST cluster of this class with the relative amount of cluster that have an overlap with ‘TSS’ or ‘Transcribed’ in parentheses
38,803 (92.6 %)
2459 (82.8 %)
3950 (87.3 %)
2821 (93.3 %)
11,207 (68.8 %)
UTR-associated RNAs (uaRNA) are identified by accumulations of CAGE tags, i.e., transcription start sites (TSS) or re-capping sites, that are independent of the underlying genes. Detectable uaRNAs overlap with the 3’UTR of their surrounding gene. Conceivably, 5’-uaRNAs might also exist. However, it is nearly impossible to distinguish the TSS (or alternative TSSs) of the surrounding gene from the 5’-end of an embedded uaRNA. In contrast, CAGE tags towards the 3’-end of a transcript exclude most or all of the coding sequence or the major part of a lncRNA, hence strongly suggesting a distinct transcript. Operationally, uaRNA candidates are defined as unspliced EST clusters that overlap a predicted TSS and the 3’UTR of a RefSeq gene but not its 5’UTR. Using TSS predicted by chromatin data as described above we identify in total of 1547 bona fide uaRNA candidates in human. Relying on TSS predicted by the FANTOM 5 consortium using CAGE tags  there are 2788 uaRNA candidates in human and 2294 in mouse. Four hundred and twelve of the uaRNA candidates in human are detected simultaneously by chromatin data and CAGE data. Of these, 160 have an ortholog unspliced EST cluster in mouse. Unfortunately, comprehensive chromatin data sets are not available for mouse.
Unannotated intergenic EST clusters
11,207 human and 6945 mouse uEST clusters are located more then 5 kb away from any RefSeq gene. In human 4986 of these overlap non-RefSeq annotated mRNAs. Removing these leaves 6222 human and 3412 mouse candidates of novel, predominantly unspliced genes.
Recent retropseudogenes  may lead to mapping artifacts since ESTs might be mapped erroneously to the intron-free locus instead of a spliced alignment to the locus of the functional gene. Although revived retrogenes exist and can be functional  they cannot be distinguished reliably from mapping artifacts. Hence we compared the sequences of the cluster against the “Retroposed Genes, Including Pseudogenes”-track from UCSC genome browser in order to identify and remove such cases from further analysis. Seven hundred and sixty loci (12 %) of the candidates are identified as deriving from protein-coding genes. In the mouse data set, the relative number is higher: 864 (25 %) of the candidates show significant similarities from coding genes. The majority of the remaining unknown transcripts are probably long extensions of 3’UTRs of spliced genes. In many of these cases we might be seeing uaRNAs.
To further investigate the 5462 intergenic unspliced EST cluster that were not recognized as likely retrogenes, we used RNAz 2.1 [80, 81] to detect signatures of purifying selection on RNA secondary structure and RNAcode  to find evidence for conserved protein-coding regions. We used multiple genome alignments taken from the UCSC genome browser as input. This data set covers 4,059,185 nt and yields 1104 RNAz hits with a classification probability P RNAz >0.5 and 371 with P RNAz >0.9. In mouse the 1,550,680 nt are covering 249 RNAz hits with a classification probability P RNAz >0.5 and 77 with P RNAz >0.9. Compared to the predicted 6880 low confidence and 2259 high confidence hits in the 30 Mb long ENCODE regions , the human result amounts to a very moderate enrichment by a factor of 1.2. It remains unclear whether this enrichment is confounded e.g. by the restriction of unspliced EST cluster to genomic regions with relatively high expression. The enrichment of structured RNA elements could be interpreted as hinting at a function of long unspliced regions in specific binding with proteins. An example for a conserved secondary structure element is shown in Additional file 1: Figure S7.
Using RNAcode, we detected 632 regions with a p-value p RNAcode <0.01 in human. Two hundred and thirty-two of these regions have a length of at least 60 nucleotides. In mouse there are 652 regions with a p-value p RNAcode <0.01. Two hundred and sixty-five have a length of at least 60 nt. Shorter regions mostly appear to be artifacts caused by short segments of coding region in retrogenes.
Overlapping all unspliced EST cluster with a recent version of GENCODE lncRNAs (v22 and vM5) we noticed that 15.1 % of human and 8.6 % of mouse uEST cluster overlap with annotated lncRNAs. Interestingly, human intergenic cluster are enriched in this set (3281, i.e., 29 % of all IGR uEST cluster). In mouse, the 642 IGR cluster (9 %) overlapping GENCODE lncRNAs closely matches the relative amount of IGR uEST cluster. The specific increase in human might be due to a better annotation compared to mouse.
Human intergenic unspliced EST cluster are the class with the most cluster that are not found to be expressed by the ENCODE RNA-seq data . In total, 96.6 % of all human unspliced EST cluster are supported by RNA-seq data, defined as being covered by at least 10 reads of ENCODE RNA-seq data. 1723, of 3519 cluster in total that are not covered, belong to the class IGR, the numbers for all classes can be found in Additional file 1: Table S4. This indicates that this class contains a higher amount of low expressed genes that are more difficult to detect.
Unsurprisingly, the class Totally EXonic (TEX) comprises the most conserved cluster (2464). One reason for that is the fact that detection of conservation between already known genes is much easier. The same goes for Totally INtronic (TIN) cluster, see Additional file 1: Figure S9 for an example. However, we also found 73 conserved intergenic cluster. One example can be seen in Additional file 1: Figure S10. 42 of the conserved uEST cluster in human and 28 in mouse overlap regions previously shown to have enhancer activity from the VISTA Enhancer Browser. The conservation between human and mouse implicitly hints at functional transcripts that have to be protected from harmful mutations.
A more global view on the sequence conservation of the unspliced EST cluster can be gained by using multiz alignments from the UCSC genome browser. In our analysis we use both the 100way-alignment with human as reference and the 60way mouse-centered alignment. For both of them phastCons scores [84, 85], which indicate the level of conservation of the reference species, are available through the UCSC genome browser. The phastCons scores, which are in the range between 0 and 1, represent probabilities of negative selection. The underlying model corrects for base composition.
Unspliced EST data, although typically disregarded in transcriptome analysis, can provide interesting insights into the structure of human transcriptomes. They outline a major component of transcriptional output that totally or at least partially escapes splicing. Nevertheless, this valuable resource has not been mined comprehensively in the past. So far, only the partially and totally intronic transcripts (PINs and TINs), which constitute more than half of the clusters of unspliced ESTs, have received attention in systematic studies [24–26]. Our analysis confirmed the conclusions drawn from the literature on these abundant class of transcripts in human and mouse. In addition, we find direct evidence of the independent transcription of PINs and TINs based on CAGE tag and chromatin data. Several promising candidates for previously unknown long RNAs could be found.
Many of the unspliced EST clusters form extensions of the UTRs of well-annotated genes. On the 5’-side they provide additional information about the transcriptional start sites (TSS) of the annotated RefSeq genes themselves. Not surprisingly, the majority of clusters with clear support for a TSS from CAGE data fall into this class. More interestingly, however, the relative majority of UTR extensions is located at the 3’-end of RefSeq genes and forms typically extensive, several kb-long extensions. In line with , many of these 3’UTR extensions contain a TSS for an independently transcribed uaRNA. We identify 2788 likely uaRNAs in human and 2294 in mouse with strong support from the CAGE data. In addition 1547 uaRNA candidates in human were detected using chromatin data. Four hundred and twelve candidates in human were predicted by both methods and 160 of them have an ortholog cluster in mouse.
CAGE data was also used to analyze the expression in 156 different primary human cell types. Most of the unspliced EST cluster are expressed in less than 10 different cell types, while only a small fraction of very abundant ones in expressed in all cell types. This distribution is similar to that of known RefSeq genes but more skewed towards specifically expressed loci.
This suggest that hundreds or even thousands of the 3’UTR extensions are EST clusters overlapping the 3’UTR of a RefSeq gene but originate rather from independent transcripts than a single 3’UTR. In fact, 227 cluster overlapping 3’UTR and downstream region of a gene and a similar number of DT uEST clusters carry histone marks characteristic for TSS, see Table 3. In addition, it is demonstrated in  that developmental-stage- and tissue-specific cleavage and subsequent secondary capping of 3’UTR is not an infrequent phenomenon. The production of such unspliced RNAs may explain the 3R uEST clusters that have CAGE tags but no promoter-like histone features.
Although TINs and PINs have been reported to be transcribed mostly independently of their surrounding gene, only a small fraction has a sufficient concentration of CAGE tags for a recognizable TSS. However there are more than 3000 human TINs overlapping a chromatin mark indicating a TSS.
The remaining set of unspliced EST clusters outside a 5 kb range around RefSeq genes comprises 11 % of the data. About a tenth (a forth in mouse) of these clusters overlaps retrogenes and retropseudogenes. For these cases, it is often impossible to distinguish between truly expressed loci and mapping artifacts of ESTs arising from the spliced original of the gene. The manual inspection of a random sample of the remaining cases shows that a large fraction of these EST clusters constitute even larger extensions of 3’UTRs. In many cases they also appear to be long 3’UTRs or uaRNAs belonging to previously unannotated, typically non-coding, spliced transcripts.
To date little is known about the evolutionary age of long ncRNAs. In our study only 15–20 % of the unspliced EST cluster are conserved between human and mouse. This fraction is lower than the numbers reported for spliced lncRNAs [12–15]. It nevertheless indicates a substantial level of conservation that is corroborated by elevated levels of sequence conservation in most subclasses of uEST clusters. Ignoring technical limitations preventing us from detecting a larger number, this pattern of conservation suggests a species or lineage-specific function of the remaining transcripts. The 5388 pairs of conserved unspliced EST clusters with the same classification in human and mouse are likely to be functionally important and we expect them be present also in other mammalian species.
In summary, the analysis of unspliced ESTs uncovers a largely unexplored realm of long transcripts. Including many previously unknown transcripts. The frequently postulated connection between lack of splicing and nuclear retention, see e.g.  suggests that this class of transcripts might be involved in chromatin organization and possibly other mechanisms of epigenetic control. Further experimental investigation of the long transcripts could provide great insight into these mechanisms.
Summary of human EST data. Analysis of the EST annotation track of hg19 and mm10 downloaded from the UCSC genome browser in April 2013. The numbers indicate successive filters. Which means just the ESTs which do not have introns larger than 30 nucleotides are included in the number of introns with more than 5 % mismatches and so on
ESTs on Haplotypes
> 30n t gaps
> 5 % mismatches
Overlap with NUMTs
Unspliced EST cluster
A possible source of contamination is nuclear encoded mitochondrial DNA (NUMT) . Since the mitochondrial transcripts remain unspliced at least in mammals , it is impossible to reliably distinguish unspliced ESTs mapping to recent NUMTs from fragments of mitochondrial transcripts. An annotation track for human and mouse NUMTs was recently provided in .
The gene locations and structure were taken from the ‘RefSeq Genes’ track  downloaded from the UCSC genome browser. Multiple genome alignments in maf format for human (46way) and mouse (60way) have also been downloaded using the UCSC table browser . Input alignments for RNAcode have to be filtered for alignments with a length of less than 21 nucleotides for computational reasons. These alignments were discarded. The genome segmentation based on ENCODE data for all six cell lines (GM12878, H1-hESC, K562, HeLa-S3, HepG2, HUVEC) was downloaded from the ENCODE public hub at the UCSC genome browser (8th of January 2014). The combined segmentation was used. Coordinates of chromatin-associated RNAs (CARs) were taken from the supplemental material of  and lifted over to hg19.
Long non-coding RNA data was downloaded from the GENCODE web page  for the version v22 in human and vM5 in mouse. Enhancer RNAs were downloaded from the supplement of Andersson et al. in the FANTOM web browser . The robust set of transcribed enhancers was used.
For most of the analyses which checked if a set of unspliced EST cluster overlaps another set of genomic intervals bedtools  (v2.20.1) were used. If the data was present the functionality of the Table Browser integrated in the UCSC genome browser was used. For the classification a custom perl script was developed. All statistical analyses were performed using the standard R function (R version 3.0.2) . The Vienna RNA package was used in version 2.1.9, RNAz in version 2.1 and RNAcode in version 0.3. Programs were combined using standard bash and awk commands. Additional information on used commands and scripts is available on request.
We define unspliced ESTs that overlap each other or that are separated by at most 30nt as members of the same cluster. In order to avoid artifacts of both the experimental and the computation procedures we only consider unspliced EST cluster comprising at least three individual ESTs, leaving 104,980 human and 66,109 mouse unspliced EST cluster for further analysis. ESTs overlapping a NUMT have been discarded before the cluster analysis. In order to computer overlaps with existing annotation we used the functionality of the Table Browser integrated in the UCSC genome browser and the “Operate on Genomic Intervals”-tools of the Galaxy Browser  as well as BEDTools .
TEX - Overlap exclusively with a RefSeq exons including UnTranslated Regions (UTR) and non-coding RNA
TIN - Overlap exclusively with a RefSeq intron
5’PIN - Overlap exclusively with an adjoining exon-intron pair, where the exon is located upstream of the intron
3’PIN - Overlap exclusively with an adjoining exon-intron pair, where the exon is located downstream of the intron
rI - Overlap exclusively with one intron and both adjacent exons
5’R - Overlap with a RefSeq exon and the 5 kb upstream region but not with downstream region or 3’UTR
3’R - Overlap with a RefSeq exon and the 5 kb downstream region but not with upstream region or 5’UTR
UT - Overlap exclusively with 5 kb upstream region
DT - Overlap exclusively with 5 kb downstream region
IGR - No overlap with any RefSeq gene or the 5 kb upstream/downstream region
NO_CLASS - None of the other classes could be applied
If the uEST cluster was determined to have the same class for all overlapping RefSeq genes, it was assigned to this class. In case of conflicting data it was assigned to the class “NO_CLASS”. An exception was made if the conflict was between one of the classes TEX, TIN, 5’ PIN, 3’ PIN, rI, 5’R and 3’R and one of the classes UT or DT. In this case the overlap with an additional upstream/downstream region was ignored. Otherwise a class assignment would be impossible in regions with a high gene density. To avoid misclassifications due to small mistakes in the borders of RefSeq gene components the overlap between the uEST cluster and the gene component had to be more than 20 nucleotides to be considered. Note that the classes are not overlapping each other. Every unspliced EST cluster is assigned to exactly one class, with all ambiguous cases collected in “NO_CLASS” an excluded from all class-specific analyses.
TSS peaks predicted by FANTOM 5 for the hg19 and mm9 assemblies were downloaded from the FANTOM 5 public hub in the UCSC genome browser. The mm9 coordinates were converted to mm10 using the UCSC liftOver tool . We used the permissive set. If an unspliced EST cluster overlaps with an predicted forward TSS peak in its first 84 nucleotides but not with an reverse TSS peak it is defined as a forward oriented transcript. If a cluster overlaps a TSS peak on the opposite strand in its last 84 nucleotides but not with a forward TSS peak its defined as reverse oriented transcript. Unspliced EST cluster which are detected as forward and reverse oriented simultaneously are disregarded.
The cell-specific CAGE data generated by FANTOM 5 was downloaded from the ’Fantom Web Resource’ [97, 102]. These comprise the supplementary data of ref. . The same analysis as above was performed individually for all 156 cell types.
Chromatin architecture analysis
Annotation for six different cell lines (GM12878, H1-hESC, K562, HeLa-S3, HepG2, HUVEC) are available. We took all cell lines together for this analysis. Every position in the genome was assigned to one of the following classes: Predicted promoter region including TSS (TSS), Predicted promoter flanking region (PF), Predicted enhancer (E), Predicted weak enhancer or open chromatin cis regulatory element (WE), CTCF enriched element (CTCF), Predicted transcribed region (T) and Predicted Repressed or Low Activity region (R).
ENCODE RNA-seq data
The full set of ENCODE RNA-seq data  was downloaded from the EBI web server (http://ftp.ebi.ac.uk/pub/databases/ensembl/encode/integration_data_jan2011/byDataType/rna_signal/jan2011/hub/). It contains 351 bigWig-files which were converted to the format bedGraph using bigWigToBedGraph, a utility from the UCSC genome browser website. BedGraph can be used as an input for bedtools. All individual files were intersected with the list of unspliced EST cluster using bedtools intersect. The output consisted of one file for each of the 351 RNA-seq libraries containing the amount of reads overlapping each unspliced EST cluster. We reported the numbers of unspliced EST cluster that overlapped with at least 10 reads in at least one library.
For detecting conserved regions we use already existing pairwise alignments between human and mouse generated by the UCSC/Penn State Bioinformatics comparative genomics alignment pipeline . The human-mouse alignment was created by aligning both genomes using blastz . The coordinates were afterwards corrected using in-house scripts. The mouse-human alignment was done by using chainSwap “to translate hg19-reference blastz alignment to mm10 into mm10-referenced chains aligned to hg19.” . Interestingly the mm10-based alignment leads to more ortholog clusters than the original one.
The chain-files were used as input to liftOver , a tool original developed to switch between genome assemblies. The human unspliced EST cluster were lifted over to mm10 by using the hg19-reference alignment, mouse unspliced EST cluster were lifted over to hg19 by using the mm10-referenced alignment, respectively. We combined all ortholog pairs detected by one of the liftOver runs.
The multiz 100way-alignment with human as reference and the 60way mouse-centered alignment from UCSC Genome Browser was used for this analysis. For every unspliced EST cluster the average phastCons score [84, 85] for the region was determined using the UCSC table browser . The phastCons scores represent probabilities of negative selection in a range from 0 to 1. They are derived from the two-state phylo-HMM and are defined as the posterior probability that the corresponding alignment column was generated by the conserved state (rather than the non-conserved state) of the phylo-HMM, given the model parameters and the multiple alignment . The phastCons score also considers the flanking columns in an alignment. The average score per cluster was used to visualize the mean phastCons score per class in Fig. 3.
In case of cluster in intergenic regions, the ones with a significant RNAcode signal have been removed. The full phastCons statistic can be found in Additional file 1: Table S5.
This work was funded by the Deutsche Forschungsgemeinschaft, Proj.No. STA 850/16-1.
- The ENCODE Project Consortium. Identification and analysis of functional elements in 1 % of the human genome by the ENCODE pilot project. Nature. 2007; 447:799–16.PubMed CentralGoogle Scholar
- Maeda N, Kasukawa T, Oyama R, Gough J, Frith M, Engström PG, et al. Transcript annotation in FANTOM3: Mouse Gene Catalog based on physical cDNAs. PLoS Genet. 2006; 2:62.Google Scholar
- Clark MB, Amaral PP, Schlesinger FJ, Dinger ME, Taft RJ, Rinn JL, et al. The reality of pervasive transcription. PLoS Biol. 2011; 9:1000625.Google Scholar
- The ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012; 489:57–74.PubMed CentralGoogle Scholar
- Guttman M, Garber M, Levin JZ, Donaghey J, Robinson J, Adiconis X, et al. Ab initio reconstruction of cell type-specific transcriptomes in mouse reveals the conserved multi-exonic structure of lincRNAs. Nat Biotechnol. 2010; 28:503–10.PubMed CentralPubMedGoogle Scholar
- Kapranov P, St Laurent G, Raz T, Ozsolak F, Reynolds CP, Sorensen PH, et al. The majority of total nuclear-encoded non-ribosomal RNA in a human cell is ‘dark matter’ un-annotated RNA. BMC Biol. 2010; 8:149.PubMed CentralPubMedGoogle Scholar
- St Laurent G, Vyatkin Y, Kapranov P. Dark matter RNA illuminates the puzzle of genome-wide association studies. BMC Med. 2014; 12:97.PubMed CentralPubMedGoogle Scholar
- Ponjavic J, Ponting CP, Lunter G. Functionality or transcriptional noise? Evidence for selection within long noncoding RNAs. Genome Res. 2007; 17:556–65.PubMed CentralPubMedGoogle Scholar
- Marques AC, Ponting CP. Catalogues of mammalian long noncoding RNAs: modest conservation and incompleteness. Genome Biol. 2009; 10:124.Google Scholar
- Guttman M, Amit I, Garber M, French C, Lin MF, Feldser D, et al. Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals. Nature. 2009; 458:223–7.PubMed CentralPubMedGoogle Scholar
- Pang KC, Frith MC, Mattick JS. Rapid evolution of noncoding RNAs: lack of conservation does not mean lack of function. Trends Genet. 2006; 22:1–5.PubMedGoogle Scholar
- Washietl S, Kellis M, Garber M. Evolutionary dynamics and tissue specificity of human long noncoding RNAs in six mammals. Genome Res. 2014; 24:616–28.PubMed CentralPubMedGoogle Scholar
- Necsulea A, Soumillon M, Warnefors M, Liechti A, Daish T, Zeller, U, et al. The evolution of lncrna repertoires and expression patterns in tetrapods. Nature. 2014; 505(7485):635–40.PubMedGoogle Scholar
- Managadze D, Lobkovsky AE, Wolf YI, Shabalina SA, Rogozin IB, Koonin EV. The vast, conserved mammalian lincRNome. PLoS Comput Biol. 2013; 9:1002917.Google Scholar
- Nitsche A, Rose D, Fasold M, Reiche K, Stadler PF. Comparison of splice sites reveals that long non-coding RNAs are evolutionarily well conserved. RNA. 2015 May; 21(5):801–12. doi:10.1261/rna.046342.114. Epub 2015 Mar 23.
- Hiller M, Findeiß S, Lein S, Marz M, Nickel C, Rose D, et al. Conserved introns reveal novel transcripts in Drosophila melanogaster. Genome Res. 2009; 19:1289–300.PubMed CentralPubMedGoogle Scholar
- Chodroff RA, Goodstadt L, Sirey TM, Oliver PL, Davies KE, Green ED, et al. Long noncoding RNA genes: conservation of sequence and brain expression among diverse amniotes. Genome Biol. 2010; 11:72.Google Scholar
- Wang KC, Chang HY. Molecular mechanisms of long noncoding RNAs. Mol Cell. 2011; 43:904–14.PubMed CentralPubMedGoogle Scholar
- Tilgner H, Knowles DG, Johnson R, Davis CA, Chakrabortty S, Djebali S, et al. Deep sequencing of subcellular RNA fractions shows splicing to be predominantly co-transcriptional in the human genome but inefficient for lncRNAs. Genome Res. 2012; 22:1616–25.PubMed CentralPubMedGoogle Scholar
- Louhichi A, Fourati A, Rebaï A. IGD: a resource for intronless genes in the human genome. Gene. 2011; 488:35–40.PubMedGoogle Scholar
- Lei H, Dias AP, Reed R. Export and stability of naturally intronless mRNAs require specific coding region sequences and the TREX mRNA export complex. Proc Natl Acad Sci USA. 2011; 108:17985–90.PubMed CentralPubMedGoogle Scholar
- Shabalina SA, Ogurtsov AY, Spiridonov AN, Novichkov PS, Spiridonov NA, Koonin EV. Distinct patterns of expression and evolution of intronless and intron-containing mammalian genes. Mol Biol Evol. 2010; 27:1745–9.PubMed CentralPubMedGoogle Scholar
- Marzluff WF, Wagner EJ, Duronio RJ. Metabolism and regulation of canonical histone mRNAs: life without a poly(A) tail. Nat Rev Genet. 2008; 9:843–54.PubMed CentralPubMedGoogle Scholar
- Nakaya HI, Amaral PP, Louro R, Lopes A, Fachel AA, Moreira YB, et al. Genome mapping and expression analyses of human intronic noncoding RNAs reveal tissue-specific patterns and enrichment in genes related to regulation of transcription. Genome Biol. 2007; 8:43.Google Scholar
- Louro R, Nakaya HI, Amaral PP, Festa F, Sogayar MC, da Silva AM, et al. Androgen responsive intronic non-coding RNAs. BMC Biol. 2007; 5:4.PubMed CentralPubMedGoogle Scholar
- Louro R, El-Jundi T, Nakaya HI, Reis EM, Verjovski-Almeida S. Conserved tissue expression signatures of intronic noncoding RNAs transcribed from human and mouse loci. Genomics. 2008; 92:18–25.PubMedGoogle Scholar
- Rinn JL, Euskirchen G, Bertone P, Martone R, Luscombe NM, Hartman S, et al. The transcriptional activity of human chromosome 22. Genes Dev. 2003; 17:529–40.PubMed CentralPubMedGoogle Scholar
- Reis EM, Nakaya HI, Louro R, Canavez FC, Flatschart AV, Almeida, GT, et al. Antisense intronic non-coding RNA levels correlate to the degree of tumor differentiation in prostate cancer. Oncogene. 2004; 23:6684–92.PubMedGoogle Scholar
- Beckedorff FC, Ayupe AC, Crocci-Souza R, Amaral MS, Nakaya HI, Soltys DT, et al. The intronic long noncoding RNA ANRASSF1 recruits PRC2 to the RASSF1A promoter, reducing the expression of RASSF1A and increasing cell proliferation. PLoS Genet. 2013; 9:1003705.Google Scholar
- Gardner EJ, Nizami ZF, Talbot Jr CC, Gall JG. Stable intronic sequence RNA (sisRNA), a new class of noncoding RNA from the oocyte nucleus of Xenopus tropicalis. Genes Dev. 2012; 26:2550–9.PubMed CentralPubMedGoogle Scholar
- Louro R, Smirnova AS, Verjovski-Almeida S. Long intronic noncoding RNA transcription: expression noise or expression choice?. Genomics. 2009; 93:291–8.PubMedGoogle Scholar
- Tahira AC, Kubrusly MS, Faria MF, Dazzani B, Fonseca RS, Maracaja-Coutinho V, et al. Long noncoding intronic RNAs are differentially expressed in primary and metastatic pancreatic cancer. Mol Cancer. 2011; 10:141.PubMed CentralPubMedGoogle Scholar
- Whitehead J, Pandey GK, Kanduri C. Regulation of the mammalian epigenome by long noncoding RNAs. Biochim Biophys Acta. 2009; 1790:936–47.PubMedGoogle Scholar
- Paul IJ, Duerksen JD. Chromatin-associated RNA content of heterochromatin and euchromatin. Mol Cell Biochem. 1975; 9:9–16.PubMedGoogle Scholar
- Rodríguez-Campos A, Azorín F. RNA is an integral component of chromatin that contributes to its structural organization. PLoS ONE. 2007; 2:1182.Google Scholar
- Mondal T, Rasmussen M, Pandey GK, Isaksson A, Kanduri C. Characterization of the RNA content of chromatin. Genome Res. 2010; 20:899–907.PubMed CentralPubMedGoogle Scholar
- Magistri M, Faghihi MA, St. Laurent 3rd G, Wahlestedt C. Regulation of chromatin structure by long noncoding RNAs: focus on natural antisense transcripts. Trends Genet. 2012; 28:389–96.Google Scholar
- Schubert T, Längst G. Changes in higher order structures of chromatin by RNP complexes. RNA Biol. 2013; 10:175–9.PubMed CentralPubMedGoogle Scholar
- Kapranov P, Cheng J, Dike S, Nix D, Duttagupta R, Willingham AT, et al. RNA maps reveal new RNA classes and a possible function for pervasive transcription. Science. 2007; 316:1484–8.PubMedGoogle Scholar
- Wang X, Arai S, Song X, Reichart D, Du K, Pascual G, et al. Induced ncRNAs allosterically modify RNA-binding proteins in cis to inhibit transcription. Nature. 2008; 454:126–30.PubMed CentralPubMedGoogle Scholar
- Hamilton RS, Davis I. Identifying and searching for conserved RNA localisation signals. Methods Mol Biol. 2011; 714:447–66.PubMed CentralPubMedGoogle Scholar
- Dassi E, Quattrone A. Tuning the engine: an introduction to resources on post-transcriptional regulation of gene expression. RNA Biol. 2012; 9:1224–32.PubMed CentralPubMedGoogle Scholar
- Jung H, Gkogkas CG, Sonenberg N, Holt CE. Remote control of gene function by local translation. Cell. 2014; 157:26–40.PubMed CentralPubMedGoogle Scholar
- Mercer TR, Wilhelm D, Dinger ME, Soldà G, Korbie DJ, Glazov EA, et al. Expression of distinct RNAs from 3’ untranslated regions. Nucleic Acids Res. 2011; 2393–2403:39.Google Scholar
- Flynn RA, Almada AE, Zamudio JR, Sharp PA. Antisense RNA polymerase II divergent transcripts are P-TEFb dependent and substrates for the RNA exosome. Proc Natl Acad Sci. 2011; 108(26):10460–5.PubMed CentralPubMedGoogle Scholar
- Niazi F, Valadkhan S. Computational analysis of functional long noncoding RNAs reveals lack of peptide-coding capacity and parallels with 3’utrs. RNA. 2012; 18:825–43.PubMed CentralPubMedGoogle Scholar
- Ebert MS, Sharp PA. Emerging roles for natural microRNA sponges. Curr Biol. 2010; 20:858–61.Google Scholar
- Bak RO, Mikkelsen JG. miRNA sponges: soaking up miRNAs for regulation of gene expression. Wiley Interdiscip Rev RNA. 2014; 5:317–33.PubMedGoogle Scholar
- Hentze MW, Preiss T. Circular RNAs: splicing’s enigma variations. EMBO J. 2013; 32:923–5.PubMed CentralPubMedGoogle Scholar
- Sasaki YTF, Ideue T, Sano M, Mituyama T, Hirose T. MEN ε/β noncoding RNAs are essential for structural integrity of nuclear paraspeckles. Proc Natl Acad Sci USA. 2009; 106:2525–30.PubMed CentralPubMedGoogle Scholar
- Sunwoo H, Dinger ME, Wilusz JE, Amaral PP, Mattick JS, Spector DL. MEN ε/β nuclear-retained non-coding RNAs are up-regulated upon muscle differentiation and are essential components of paraspeckles. Genome Res. 2009; 19:347–59.PubMed CentralPubMedGoogle Scholar
- Mao YS, Sunwoo H, Zhang B, Spector DL. Direct visualization of the co-transcriptional assembly of a nuclear body by noncoding RNAs. Nat Cell Biol. 2011; 13:95–101.PubMed CentralPubMedGoogle Scholar
- Hutchinson J, Ensminger AW, Clemson CM, Lynch CR, Lawrence JB, Chess A. A screen for nuclear transcripts identifies two linked noncoding RNAs associated with SC35 splicing domains. BMC Genomics. 2007; 8:39.PubMed CentralPubMedGoogle Scholar
- Stadler PF, editors. Evolution of the long non-coding RNAs MALAT1 and MEN β/ ε In: Ferreira CE, Miyano S, Stadler PF, editors. Advances in Bioinformatics and Computational Biology, 5th Brazilian Symposium on Bioinformatics. Lecture Notes in Computer Science, vol.6268. Heidelberg: Springer: 2010. p. 1–12.Google Scholar
- Chung S, Nakagawa H, Uemura M, Piao L, Ashikawa K, Hosono N, et al. Association of a novel long non-coding RNA in 8q24 with prostate cancer susceptibility. Cancer Sci. 2011; 102:245–52.PubMedGoogle Scholar
- Hackermüller J, Reiche K, Otto C, Hösler N, Blumert C, Brocke-Heidrich K, et al. Cell cycle, oncogenic and tumor suppressor pathways regulate numerous long and macro non-protein coding RNAs. Genome Biol. 2014; 15:48.Google Scholar
- Seidl CIM, Stricker SH, Barlow DP. The imprinted Air ncRNA is an atypical RNAPII transcript that evades splicing and escapes nuclear export. EMBO J. 2006; 25:1–11.Google Scholar
- Stricker SH, Steenpass L, Pauler FM, Santoro F, Latos PA, Huang R, et al. Silencing and transcriptional properties of the imprinted Airn ncRNA are independent of the endogenous promoter. EMBO J. 2008; 27:3116–28.PubMed CentralPubMedGoogle Scholar
- Redrup L, Branco MR, Perdeaux ER, Krueger C, Lewis A, Santos F, et al. The long noncoding RNA Kcnq1ot1 organises a lineage-specific nuclear domain for epigenetic gene silencing. Development. 2009; 136:525–30.PubMed CentralPubMedGoogle Scholar
- Nordström KJ, Mirza MA, Almén MS, Gloriam DE, Fredriksson R, Schiöth HB. Critical evaluation of the FANTOM3 non-coding RNA transcripts. Genomics. 2009; 94:169–76.PubMedGoogle Scholar
- Engelhardt J, Stadler PF. Hidden treasures in unspliced est data. Theory Biosci. 2012; 131(1):49–57.PubMedGoogle Scholar
- Greenbaum D, Colangelo C, Williams K, Gerstein M. Comparing protein abundance and mRNA expression levels on a genomic scale. Genome Biol. 2003; 4:117.PubMed CentralPubMedGoogle Scholar
- Gry M, Rimini R, Strömberg S, Asplund A, Pontén F, Uhlén M, et al. Correlations between RNA and protein expression profiles in 23 human cell lines. BMC Genomics. 2009; 10:365.PubMed CentralPubMedGoogle Scholar
- Thorrez L, Tranchevent LC, Chan HJ, Moreau Y, Schuit F. Detection of novel 3’ untranslated region extensions with 3’ expression microarrays. BMC Genomics. 2010; 11:205.PubMed CentralPubMedGoogle Scholar
- Proudfoot NJ. Ending the message: poly(A) signals then and now. Genes Dev. 2011; 25:1770–82.PubMed CentralPubMedGoogle Scholar
- Li Y, Sun Y, Fu Y, Li M, Huang G, Zhang C, et al. Dynamic landscape of tandem 3’ UTRs during zebrafish development. Genome Res. 2012; 22:1899–906.PubMed CentralPubMedGoogle Scholar
- Lianoglou S, Garg V, Yang JL, Leslie CS, Mayr C. Ubiquitously transcribed genes use alternative polyadenylation to achieve tissue-specific expression. Genes Dev. 2013; 27:2380–96.PubMed CentralPubMedGoogle Scholar
- Shiraki T, Kondo S, Katayama S, Waki K, Kasukawa T, Kawaji H, et al. Cap analysis gene expression for high-throughput analysis of transcriptional starting point and identification of promoter usage. Proc Natl Acad Sci U S A. 2003; 100(26):15776–81.PubMed CentralPubMedGoogle Scholar
- The FANTOM Consortium and others. A promoter-level mammalian expression atlas. Nature. 2014; 507(7493):462–70.Google Scholar
- Ørom UA, Shiekhattar R. Long noncoding RNAs usher in a new era in the biology of enhancers. Cell. 2013; 154:1190–3.PubMed CentralPubMedGoogle Scholar
- Lai F, Shiekhattar R. Enhancer RNAs: the new molecules of transcription. Curr Opin Genet Dev. 2014; 25:38–42.PubMed CentralPubMedGoogle Scholar
- Lam MT, Li W, Rosenfeld MG, Glass CK. Enhancer RNAs and regulated transcriptional programs. Trends Biochem Sci. 2014; 39:170–82.PubMed CentralPubMedGoogle Scholar
- Kim TK, Hemberg M, Gray JM. Enhancer RNAs: a class of long noncoding RNAs synthesized at enhancers. Cold Spring Harb Perspect Biol. 2015; 7:018622.Google Scholar
- Andersson R, Gebhard C, Miguel-Escalada I, Hoof I, Bornholdt J, Boyd M, et al. An atlas of active enhancers across human cell types and tissues. Nature. 2014; 507(7493):455–61.PubMedGoogle Scholar
- Ohmiya H, Vitezic M, Frith MC, Itoh M, Carninci P, Forrest AR, et al. RECLU: a pipeline to discover reproducible transcriptional start sites and their alternative regulation using capped analysis of gene expression (cage). BMC Genomics. 2014; 15(1):269.PubMed CentralPubMedGoogle Scholar
- Derrien T, Johnson R, Bussotti G, Tanzer A, Djebali S, Tilgner H, et al. The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome Res. 2012; 22:1775–89.PubMed CentralPubMedGoogle Scholar
- Hoffman MM, Ernst J, Wilder SP, Kundaje A, Harris RS, Libbrecht M, et al. Integrative annotation of chromatin elements from encode data. Nucleic Acids Res. 2013; 41(2):827–41.PubMed CentralPubMedGoogle Scholar
- Devor EJ, Moffat-Wilson KA. Molecular and temporal characteristics of human retropseudogenes. Hum Biol. 2003; 75:661–72.PubMedGoogle Scholar
- Yu Z, Morais MDIvanga, Harrison PM. Analysis of the role of retrotransposition in gene evolution in vertebrates. BMC Bioinformatics. 2007; 8:308.PubMed CentralPubMedGoogle Scholar
- Washietl S, Hofacker IL, Stadler PF. Fast and reliable prediction of noncoding RNAs. Proc Natl Acad Sci USA. 2005; 102:2454–9.PubMed CentralPubMedGoogle Scholar
- Gruber AR, Findeiß S, Washietl S, Hofacker IL, Stadler PF. RNAz 2.0: improved noncoding RNA detection. Pac Symp Biocomput. 2010; 15:69–79.Google Scholar
- Washietl S, Findeiß S, Müller S, Kalkhof S, von Bergen M, Hofacker IL, et al. RNAcode: robust prediction of protein coding regions in comparative genomics data. RNA. 2011; 17:578–94.PubMed CentralPubMedGoogle Scholar
- ENCODE Project Consortium. An integrated encyclopedia of dna elements in the human genome. Nature. 2012; 489(7414):57–74.Google Scholar
- Siepel A, Bejerano G, Pedersen JS, Hinrichs A, Hou M, Rosenbloom K, et al. Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes. Genome Res. 2005; 15:1034–50.PubMed CentralPubMedGoogle Scholar
- Hubisz MJ, Pollard KS, Siepel A. Phast and rphast: phylogenetic analysis with space/time models. Brief Bioinform. 2011; 12(1):41–51.PubMed CentralPubMedGoogle Scholar
- Mercer TR, Dinger ME, Bracken CP, Kolle G, Szubert JM, Korbie DJ, et al. Regulated post-transcriptional RNA cleavage diversifies the eukaryotic transcriptome. Genome Res. 2010; 20:1639–50.PubMed CentralPubMedGoogle Scholar
- Eberle AB, Hessle V, Helbig R, Dantoft W, Gimber N, Visa N. Splice-site mutations cause rrp6-mediated nuclear retention of the unspliced RNAs and transcriptional down-regulation of the splicing-defective genes. PLoS ONE. 2010; 5:11540.Google Scholar
- Karolchik D, Hinrichs AS, Furey TS, Roskin KM, Sugnet CW, Haussler D, et al. The ucsc table browser data retrieval tool. Nucleic Acids Res. 2004; 32(suppl 1):493–6.Google Scholar
- Hawkins JD. A survey on intron and exon lengths. Nucleic Acids Res. 1988; 16:9893–908.PubMed CentralPubMedGoogle Scholar
- Hong X, Scofield DG, Lynch M. Intron size, abundance, and distribution within untranslated regions of genes. Mol Biol Evol. 2006; 23:2392–404.PubMedGoogle Scholar
- Horton R, Gibson R, Coggill P, Miretti M, Allcock RJ, Almeida J, et al. Variation analysis and gene annotation of eight mhc haplotypes: the mhc haplotype project. Immunogenetics. 2008; 60(1):1–18.PubMed CentralPubMedGoogle Scholar
- Hazkani-Covo E, Zeller R, Martin W. Molecular poltergeists: Mitochondrial DNA copies (numts) in sequenced nuclear genomes. PLoS Genet. 2010; 6:1000834.Google Scholar
- Asin-Cayuela J, Gustafsson CM. Mitochondrial transcription and its regulation in mammalian cells. Trends Biochem Sci. 2007; 32:111–7.PubMedGoogle Scholar
- Tsuji J, Frith MC, Tomii K, Horton P. Mammalian NUMT insertion is non-random. Nucleic Acids Res. 2012; 40(18):9073–88.PubMed CentralPubMedGoogle Scholar
- Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007; 35:61–5.Google Scholar
- Harrow J, Frankish A, Gonzalez JM, Tapanari E, Diekhans M, Kokocinski F, et al. Gencode: the reference human genome annotation for the encode project. Genome Res. 2012; 22(9):1760–74.PubMed CentralPubMedGoogle Scholar
- Kawaji H, Severin J, Lizio M, Waterhouse A, Katayama S, Irvine KM, et al. The fantom web resource: from mammalian transcriptional landscape to its dynamic regulation. Genome Biol. 2009; 10(4):40.Google Scholar
- Quinlan AR, Hall IM. Bedtools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010; 26(6):841–2.PubMed CentralPubMedGoogle Scholar
- R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2013. http://www.r-project.org/.Google Scholar
- Goecks J, Nekrutenko A, Taylor J, Galaxy Team. Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences. Genome Biol. 2010; 11:86.Google Scholar
- Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, et al. The human genome browser at ucsc. Genome Res. 2002; 12(6):996–1006.PubMed CentralPubMedGoogle Scholar
- Kawaji H, Severin J, Lizio M, Forrest AR, van Nimwegen E, Rehli M, et al. Update of the fantom web resource: from mammalian transcriptional landscape to its dynamic regulation. Nucleic Acids Res. 2011; 39(Database issue):856–60.Google Scholar
- Chiaromonte F, Yap V, Miller W. Scoring pairwise genomic sequence alignments, 2002. In: Russ B. Altman, Keith A. Dunker, Lawrence Hunter, Kevin Lauderdale, Teri E. Klein, editors. Pac Symp Biocomput. vol. 7. USA: World Scientific Publishing Co Pte Ltd: 2002. p. 115–26.Google Scholar
- Schwartz S, Kent WJ, Smit A, Zhang Z, Baertsch R, Hardison RC, et al. Human-mouse alignments with blastz. Genome Res. 2003; 13(1):103–7.PubMed CentralPubMedGoogle Scholar
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.