Eukaryotic snoRNAs: a paradigm for gene expression flexibility

Protein kinase R and its cellular regulators in cancer: An active player or a surveillant?

Protein kinase R (PKR), originally known as an antiviral protein, senses various stresses as well as pathogen-driven double-stranded RNAs. Thereby activated PKR provokes diverse downstream events, including eIF2α phosphorylation and nuclear factor kappa-light-chain-enhancer of activated B cells activation. Consequently, PKR induces apoptosis and inflammation, both of which are highly important in cancer as much as its original antiviral role. Therefore, cellular proteins and RNAs should tightly control PKR activity. PKR and its regulators are often dysregulated in cancer and it is undoubted that such dysregulation contributes to tumorigenesis. However, PKR's precise role in cancer is still in debate, due to incomprehensible and even contradictory data. In this review, we introduce important cellular PKR regulators and discuss about their roles in cancer. Among them, we pay particular attention to nc886, a PKR repressor noncoding RNA that has been identified relatively recently, because its expression pattern in cancer can explain interesting yet obscure oncologic aspects of PKR. Based on nc886 and its regulation of PKR, we have proposed a tumor surveillance model, which reconciles contradictory data about PKR in cancer. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

H/ACA box small nucleolar RNA 7B acts as an oncogene and a potential prognostic biomarker in breast cancer

BACKGROUND

Breast cancer (BC) is the most frequent malignancy occurring in women worldwide. Emerging evidence indicates that small nucleolar RNAs (snoRNAs) play a role in tumor development. In the current study, we evaluated expression profiles and functions of snoRNAs associated with BC.

METHODS

We analyzed the expression levels of snoRNAs between breast cancer and normal tissues in TCGA database and found that SNORA7B is upregulated in BC. We confirmed this result in clinical cancer tissues and BC cell lines via qRT-PCR. Then, we investigated clinical significance in public datasets and biological function of SNORA7B using a series of in vitro gain- and loss-of-function experiments.

RESULTS

SNORA7B expression was significantly upregulated in samples from patients with BC in both public database and our clinical tissues compared to its expression in normal tissues. Meanwhile, patients with high SNORA7B expression have worse prognosis. Inhibition of SNORA7B expression impaired cell growth, proliferation, migration, and invasion via inducing apoptosis.

CONCLUSIONS

SNORA7B functions as an important oncogenic snoRNA in BC and may serve as a potential prognosis biomarker for BC.

The cellular landscape of mid-size noncoding RNA

Noncoding RNA plays an important role in all aspects of the cellular life cycle, from the very basic process of protein synthesis to specialized roles in cell development and differentiation. However, many noncoding RNAs remain uncharacterized and the function of most of them remains unknown. Mid-size noncoding RNAs (mncRNAs), which range in length from 50 to 400 nucleotides, have diverse regulatory functions but share many fundamental characteristics. Most mncRNAs are produced from independent promoters although others are produced from the introns of other genes. Many are found in multiple copies in genomes. mncRNAs are highly structured and carry many posttranscriptional modifications. Both of these facets dictate their RNA-binding protein partners and ultimately their function. mncRNAs have already been implicated in translation, catalysis, as guides for RNA modification, as spliceosome components and regulatory RNA. However, recent studies are adding new mncRNA functions including regulation of gene expression and alternative splicing. In this review, we describe the different classes, characteristics and emerging functions of mncRNAs and their relative expression patterns. Finally, we provide a portrait of the challenges facing their detection and annotation in databases. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

Two Proximally Close Priority Candidate Genes for diplopodia-1, an Autosomal Inherited Craniofacial-Limb Syndrome in the Chicken: MRE11 and GPR83

Next-generation sequencing (NGS) and expression technologies were utilized to investigate the genes and sequence elements in a 586 kb region of chicken chromosome 1 associated with the autosomal recessive diplopodia-1 (dp-1) mutation. This mutation shows a syndromic phenotype similar to known human developmental abnormalities (e.g., cleft palate, polydactyly, omphalocele [exposed viscera]). Toward our goal to ascertain the variant responsible, the entire 586 kb region was sequenced following utilization of a specifically designed capture array and to confirm/validate fine-mapping results. Bioinformatic analyses identified a total of 6142 sequence variants, which included SNPs, indels, and gaps. Of these, 778 SNPs, 146 micro-indels, and 581 gaps were unique to the UCD-Dp-1.003 inbred congenic line; those found within exons and splice sites were studied for contribution to the mutant phenotype. Upon further validation with additional mutant samples, a smaller subset (of variants [51]) remains linked to the mutation. Additionally, utilization of specific samples in the NGS technology was advantageous in that fine-mapping methodologies eliminated an additional 326 kb of sequence information on chromosome 1. Predicted and confirmed protein-coding genes within the smaller 260 kb region were assessed for their developmental expression patterns over several stages of early embryogenesis in regions/tissues of interest (e.g., digits, craniofacial region). Based on these results and known function in other vertebrates, 2 genes within 5 kb of each other, MRE11 and GPR83, are proposed as high-priority candidates for the dp-1 mutation.

Expression profiling of snoRNAs in normal hematopoiesis and AML

Small nucleolar RNAs (snoRNAs) are noncoding RNAs that contribute to ribosome biogenesis and RNA splicing by modifying ribosomal RNA and spliceosome RNAs, respectively. We optimized a next-generation sequencing approach and a custom analysis pipeline to identify and quantify expression of snoRNAs in acute myeloid leukemia (AML) and normal hematopoietic cell populations. We show that snoRNAs are expressed in a lineage- and development-specific fashion during hematopoiesis. The most striking examples involve snoRNAs located in 2 imprinted loci, which are highly expressed in hematopoietic progenitors and downregulated during myeloid differentiation. Although most snoRNAs are expressed at similar levels in AML cells compared with CD34⁺, a subset of snoRNAs showed consistent differential expression, with the great majority of these being decreased in the AML samples. Analysis of host gene expression, splicing patterns, and whole-genome sequence data for mutational events did not identify transcriptional patterns or genetic alterations that account for these expression differences. These data provide a comprehensive analysis of the snoRNA transcriptome in normal and leukemic cells and should be helpful in the design of studies to define the contribution of snoRNAs to normal and malignant hematopoiesis.

The complex enzymology of mRNA decapping: Enzymes of four classes cleave pyrophosphate bonds

The 5' ends of most RNAs are chemically modified to enable protection from nucleases. In bacteria, this is often achieved by keeping the triphosphate terminus originating from transcriptional initiation, while most eukaryotic mRNAs and small nuclear RNAs have a 5'→5' linked N7 -methyl guanosine (m7 G) cap added. Several other chemical modifications have been described at RNA 5' ends. Common to all modifications is the presence of at least one pyrophosphate bond. To enable RNA turnover, these chemical modifications at the RNA 5' end need to be reversible. Dependent on the direction of the RNA decay pathway (5'→3' or 3'→5'), some enzymes cleave the 5'→5' cap linkage of intact RNAs to initiate decay, while others act as scavengers and hydrolyse the cap element of the remnants of the 3'→5' decay pathway. In eukaryotes, there is also a cap quality control pathway. Most enzymes involved in the cleavage of the RNA 5' ends are pyrophosphohydrolases, with only a few having (additional) 5' triphosphonucleotide hydrolase activities. Despite the identity of their enzyme activities, the enzymes belong to four different enzyme classes. Nudix hydrolases decap intact RNAs as part of the 5'→3' decay pathway, DXO family members mainly degrade faulty RNAs, members of the histidine triad (HIT) family are scavenger proteins, while an ApaH-like phosphatase is the major mRNA decay enzyme of trypanosomes, whose RNAs have a unique cap structure. Many novel cap structures and decapping enzymes have only recently been discovered, indicating that we are only beginning to understand the mechanisms of RNA decapping. This article is categorized under: RNA Turnover and Surveillance > Turnover/Surveillance Mechanisms RNA Turnover and Surveillance > Regulation of RNA Stability RNA Processing > Capping and 5' End Modifications.

RNA-Mediated Disease Mechanisms in Neurodegenerative Disorders

RNA is accurately entangled in virtually all pathways that maintain cellular homeostasis. To name but a few, RNA is the "messenger" between DNA encoded information and the resulting proteins. Furthermore, RNAs regulate diverse processes by forming DNA::RNA or RNA::RNA interactions. Finally, RNA itself can be the scaffold for ribonucleoprotein complexes, for example, ribosomes or cellular bodies. Consequently, disruption of any of these processes can lead to disease. This review describes known and emerging RNA-based disease mechanisms like interference with regular splicing, the anomalous appearance of RNA-protein complexes and uncommon RNA species, as well as non-canonical translation. Due to the complexity and entanglement of the above-mentioned pathways, only few drugs are available that target RNA-based disease mechanisms. However, advances in our understanding how RNA is involved in and modulates cellular homeostasis might pave the way to novel treatments.

Guide snoRNAs: Drivers or Passengers in Human Disease?

In every domain of life, RNA-protein interactions play a significant role in co- and post-transcriptional modifications and mRNA translation. RNA performs diverse roles inside the cell, and therefore any aberrancy in their function can cause various diseases. During maturation from its primary transcript, RNA undergoes several functionally important post-transcriptional modifications including pseudouridylation and ribose 2′-O-methylation. These modifications play a critical role in the stability of the RNA. In the last few decades, small nucleolar RNAs (snoRNAs) were revealed to be one of the main components to guide these modifications. Due to their active links to the nucleoside modification, deregulation in the snoRNA expressions can cause multiple disorders in humans. Additionally, host genes carrying snoRNA-encoding sequences in their introns also show differential expression in disease. Although few reports support a causal link between snoRNA expression and disease manifestation, this emerging field will have an impact on the way we think about biomarkers or identify novel targets for therapy. This review focuses on the intriguing aspect of snoRNAs that function as a guide in post-transcriptional RNA modification, and regulation of their host genes in human disease.

Small Nucleolar RNAs Tell a Different Tale

Transcribing RNA Polymerase II interacts with multiple factors that orchestrate maturation and stabilisation of messenger RNA. For the majority of noncoding RNAs, the polymerase complex employs entirely different strategies, which usually direct the nascent transcript to ribonucleolytic degradation. However, some noncoding RNA classes use endo- and exonucleases to achieve functionality. Here we review processing of small nucleolar RNAs that are transcribed by RNA Polymerase II as precursors, and whose 5' and 3' ends undergo processing to release mature, functional molecules. The maturation strategies of these noncoding RNAs in various organisms follow a similar pattern but employ different factors and are strictly correlated with genomic organisation of their genes.

Differential expression of snoRNAs in oral squamous cell carcinomas: new potential diagnostic markers

BACKGROUND

Small nucleolar RNAs (snoRNAs) are small non-coding RNA sequences whose most studied function is ribosome biogenesis. The altered expression of snoRNA is observed in tumoral processes such as breast cancer and multiple myeloma. However, we have not found any references to snoRNAs in oral squamous cell carcinomas (OSCC) in the literature at the time this article was written.

MATERIAL AND METHODS

We have analyzed snoRNA expression in frozen OSCC tissue samples and have compared them to healthy controls. RNA was extracted from a total of eight OSCC samples and eight control samples, measuring the differential expression of small RNAs with the Affymetrix® miRNA 4.1 Array Plate microarray platform.

RESULTS

Results were analyzed using the Transcriptome Analysis Console 3.0 (TAC) software. We obtained a total of 16 deregulated snoRNAs of which one was over expressed and 15 were under expressed. SnoRNAs expression was altered in OSCC and could serve as a diagnostic marker.

Knock-Down of a Novel snoRNA in Tetrahymena Reveals a Dual Role in 5.8S rRNA Processing and Generation of a 26S rRNA Fragment

In eukaryotes, 18S, 5.8S, and 28S rRNAs are transcribed as precursor molecules that undergo extensive modification and nucleolytic processing to form the mature rRNA species. Central in the process are the small nucleolar RNAs (snoRNAs). The majority of snoRNAs guide site specific chemical modifications but a few are involved in defining pre-rRNA cleavages. Here, we describe an unusual snoRNA (TtnuCD32) belonging to the box C/D subgroup from the ciliate Tetrahymena thermophila. We show that TtnuCD32 is unlikely to function as a modification guide snoRNA and that it is critical for cell viability. Cell lines with genetic knock-down of TtnuCD32 were impaired in growth and displayed two novel and apparently unrelated phenotypes. The most prominent phenotype is the accumulation of processing intermediates of 5.8S rRNA. The second phenotype is the decrease in abundance of a ~100 nt 26S rRNA fragment of unknown function. Sequence analysis demonstrated that TtnuCD32 share features with the essential snoRNA U14 but an alternative candidate (TtnuCD25) was more closely related to other U14 sequences. This, together with the fact that the observed rRNA processing phenotypes were not similar to what has been observed in U14 depleted cells, suggests that TtnuCD32 is a U14 homolog that has gained novel functions.

Non-Coding RNAs and Hepatitis C Virus-Induced Hepatocellular Carcinoma

Hepatitis C virus (HCV) infection is a worldwide health problem and is one of the main causes of chronic hepatitis, liver cirrhosis, and hepatocellular carcinoma (HCC). Despite recent improvements, effective treatments for HCC are still missing and new tools for early detection are needed. Non-coding RNAs (ncRNAs) have emerged as important regulators of gene expression and key players in human carcinogenesis, including HCC. Aberrant expression of ncRNAs is associated with HCC metastasis, invasion, dissemination, and recurrence. This review will focus on the recent advances in ncRNA expression profiles, their dysregulation in HCV-related HCC, and the clinical perspective of ncRNA signatures for the early detection of HCC.

NanoPARE: parallel analysis of RNA 5' ends from low-input RNA

Diverse RNA 5′ ends are generated through both transcriptional and post-transcriptional processes. These important modes of gene regulation often vary across cell types and can contribute to the diversification of transcriptomes and thus cellular differentiation. Therefore, the identification of primary and processed 5′ ends of RNAs is important for their functional characterization. Methods have been developed to profile either RNA 5′ ends from primary transcripts or the products of RNA degradation genome-wide. However, these approaches either require high amounts of starting RNA or are performed in the absence of paired gene-body mRNA-seq data. This limits current efforts in RNA 5′ end annotation to whole tissues and can prevent accurate RNA 5′ end classification due to biases in the data sets. To enable the accurate identification and precise classification of RNA 5′ ends from standard and low-input RNA, we developed a next-generation sequencing-based method called nanoPARE and associated software. By integrating RNA 5′ end information from nanoPARE with gene-body mRNA-seq data from the same RNA sample, our method enables the identification of transcription start sites at single-nucleotide resolution from single-cell levels of total RNA, as well as small RNA-mediated cleavage events from at least 10,000-fold less total RNA compared to conventional approaches. NanoPARE can therefore be used to accurately profile transcription start sites, noncapped RNA 5′ ends, and small RNA targeting events from individual tissue types. As a proof-of-principle, we utilized nanoPARE to improve Arabidopsis thaliana RNA 5′ end annotations and quantify microRNA-mediated cleavage events across five different flower tissues.

Noncanonical RNA-capping: Discovery, mechanism, and physiological role debate

Recently a new type of 5'-RNA cap was discovered. In contrast to the specialized eukaryotic m⁷ G cap, the novel caps are abundant cellular cofactors like NAD⁺ . RNAs capped with cofactors are found in prokaryotes and eukaryotes. Unlike m⁷ G cap, installed by specialized enzymes, cofactors are attached by main enzyme of transcription, RNA polymerase (RNAP). Cofactors act as noncanonical initiating substrates, provided cofactor's nucleoside base-pairs with template DNA at the transcription start site. Adenosine-containing NAD(H), flavin adenine dinucleotide (FAD), and CoA modify transcripts on promoters starting with +1A. Similarly, uridine-containing cell wall precursors, for example, uridine diphosphate-N-acetylglucosamine were shown to cap RNA in vitro on +1U promoters. Noncanonical capping is a universal feature of evolutionary unrelated RNAPs-multisubunit bacterial and eukaryotic RNAPs, and single-subunit mitochondrial RNAP. Cellular concentrations of cofactors, for example, NAD(H) are significantly higher than their K_m in transcription. Yet, only a small proportion of a given cellular RNA is noncanonically capped (if at all). This proportion is a net balance between capping, seemingly stochastic, and decapping, possibly determined by RNA folding, protein binding and transcription rate. NUDIX hydrolases in bacteria and eukaryotes, and DXO family proteins eukaryotes act as decapping enzymes for noncanonical caps. The physiological role of noncanonical RNA capping is only starting to emerge. It was demonstrated to affect RNA stability in vivo in bacteria and eukaryotes and to stimulate RNAP promoter escape in vitro in Escherichia coli. NAD⁺ /NADH capping ratio may connect transcription to cellular redox state. Potentially, noncanonical capping affects mRNA translation, RNA-protein binding and RNA localization. This article is categorized under: RNA Processing > Capping and 5' End Modifications RNA Export and Localization > RNA Localization RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry.

Simultaneous sequencing of coding and noncoding RNA reveals a human transcriptome dominated by a small number of highly expressed noncoding genes

Comparing the abundance of one RNA molecule to another is crucial for understanding cellular functions but most sequencing techniques can target only specific subsets of RNA. In this study, we used a new fragmented ribodepleted TGIRT sequencing method that uses a thermostable group II intron reverse transcriptase (TGIRT) to generate a portrait of the human transcriptome depicting the quantitative relationship of all classes of nonribosomal RNA longer than 60 nt. Comparison between different sequencing methods indicated that FRT is more accurate in ranking both mRNA and noncoding RNA than viral reverse transcriptase-based sequencing methods, even those that specifically target these species. Measurements of RNA abundance in different cell lines using this method correlate with biochemical estimates, confirming tRNA as the most abundant nonribosomal RNA biotype. However, the single most abundant transcript is 7SL RNA, a component of the signal recognition particle. Structured noncoding RNAs (sncRNAs) associated with the same biological process are expressed at similar levels, with the exception of RNAs with multiple functions like U1 snRNA. In general, sncRNAs forming RNPs are hundreds to thousands of times more abundant than their mRNA counterparts. Surprisingly, only 50 sncRNA genes produce half of the non-rRNA transcripts detected in two different cell lines. Together the results indicate that the human transcriptome is dominated by a small number of highly expressed sncRNAs specializing in functions related to translation and splicing.

Eukaryotic RNA 5'-End NAD+ Capping and DeNADding

A hallmark of eukaryotic mRNAs has long been the 5'-end m7G cap. This paradigm was recently amended by recent reports that Saccharomyces cerevisiae and mammalian cells also contain mRNAs carrying a novel nicotinamide adenine dinucleotide (NAD+) cap at their 5'-end. The presence of an NAD+ cap on mRNA uncovers a previously unknown mechanism for controlling gene expression through nucleotide metabolite-directed mRNA turnover. In contrast to the m7G cap that stabilizes mRNA, the NAD+ cap targets RNA for rapid decay in mammalian cells through the DXO non-canonical decapping enzyme which removes intact NAD+ from RNA in a process termed 'deNADding'. This review highlights the identification of NAD+ caps, their mode of addition, and their functional significance in cells.

Nuclear fate of yeast snoRNA is determined by co-transcriptional Rnt1 cleavage

Small nucleolar RNA (snoRNA) are conserved and essential non-coding RNA that are transcribed by RNA Polymerase II (Pol II). Two snoRNA classes, formerly distinguished by their structure and ribonucleoprotein composition, act as guide RNA to target RNA such as ribosomal RNA, and thereby introduce specific modifications. We have studied the 5'end processing of individually transcribed snoRNA in S. cerevisiae to define their role in snoRNA biogenesis and functionality. Here we show that pre-snoRNA processing by the endonuclease Rnt1 occurs co-transcriptionally with removal of the m7G cap facilitating the formation of box C/D snoRNA. Failure of this process causes aberrant 3'end processing and mislocalization of snoRNA to the cytoplasm. Consequently, Rnt1-dependent 5'end processing of box C/D snoRNA is critical for snoRNA-dependent methylation of ribosomal RNA. Our results reveal that the 5'end processing of box C/D snoRNA defines their distinct pathway of maturation.

Small nucleolar RNA U2_19 promotes hepatocellular carcinoma progression by regulating Wnt/β-catenin signaling

Emerging evidence suggests that small nucleolar RNAs (snoRNAs) have malfunctioning roles in oncogenesis. In the present study, we investigated the role of box C/D small nucleolar RNA U2_19 (snoU2_19) in the tumorigenesis of hepatocellular carcinoma (HCC). Recently, we screened snoRNAs differential signatures by performing high-throughput small RNA sequence in HCC tissues and validated that upregulated snoU2_19 was associated with aggressive phenotypes in HCC patients. Aberrant snoU2_19 facilitated HCC cell proliferation, inhibited apoptosis and induced cell cycle progression in vitro analyses. We globally investigated the molecular mechanisms of snoU2_19 in HCC and found that snoU2_19 knockdown inhibited Wnt/β-catenin signaling pathway through inducing the translocation of β-catenin in cytoplasm. We concluded that snoU2_19 plays a pathological role in the development and progression of HCC, and is a potential therapeutic target for HCC.

RNA Cloaking by Reversible Acylation

We describe a selective and mild chemical approach for controlling RNA hybridization, folding, and enzyme interactions. Reaction of RNAs in aqueous buffer with an azide-substituted acylating agent (100-200 mm) yields several 2'-OH acylations per RNA strand in as little as 10 min. This poly-acylated ("cloaked") RNA is strongly blocked from hybridization with complementary nucleic acids, from cleavage by RNA-processing enzymes, and from folding into active aptamer structures. Importantly, treatment with a water-soluble phosphine triggers a Staudinger reduction of the azide groups, resulting in spontaneous loss of acyl groups ("uncloaking"). This fully restores RNA folding and biochemical activity.

Identification of small non-coding RNA classes expressed in swine whole blood during HP-PRRSV infection

It has been established that reduced susceptibility to porcine reproductive and respiratory syndrome virus (PRRSV) has a genetic component. This genetic component may take the form of small non-coding RNAs (sncRNA), which are molecules that function as regulators of gene expression. Various sncRNAs have emerged as having an important role in the immune system in humans. The study uses transcriptomic read counts to profile the type and quantity of both well and lesser characterized sncRNAs, such as microRNAs and small nucleolar RNAs to identify and quantify the classes of sncRNA expressed in whole blood between healthy and highly pathogenic PRRSV-infected pigs. Our results returned evidence on nine classes of sncRNA, four of which were consistently statistically significantly different based on Fisher's Exact Test, that can be detected and possibly interrogated for their effect on host dysregulation during PRRSV infections.

The circadian dynamics of small nucleolar RNA in the mouse liver

The circadian regulation of gene expression allows plants and animals to anticipate predictable environmental changes. While the influence of the circadian clock has recently been shown to extend to ribosome biogenesis, the dynamics and regulation of the many small nucleolar RNA that are required in pre-ribosomal RNA folding and modification are unknown. Using a novel computational method, we show that 18S and 28S pre-rRNA are subject to circadian regulation in a nuclear RNA sequencing time course. A population of snoRNA with circadian expression is identified that is functionally associated with rRNA modification. More generally, we find the abundance of snoRNA known to modify 18S and 28S to be inversely correlated with the abundance of their target. Cyclic patterns in the expression of a number of snoRNA indicate a coordination with rRNA maturation, potentially through an upregulation in their biogenesis, or their release from mature rRNA at the end of the previous cycle of rRNA maturation, in antiphase with the diurnal peak in pre-rRNA. Few cyclic snoRNA have cyclic host genes, indicating the action of regulatory mechanisms in addition to transcriptional activation of the host gene. For highly expressed independently transcribed snoRNA, we find a characteristic RNA polymerase II and H3K4me3 signature that correlates with mean snoRNA expression over the day.

Pancreatic gene expression during recovery after pancreatitis reveals unique transcriptome profiles

It is well known that pancreatic recovery after a single episode of injury such as an isolated bout of pancreatitis occurs rapidly. It is unclear, however, what changes are inflicted in such conditions to the molecular landscape of the pancreas. In the caerulein hyperstimulation model of pancreatitis, the murine pancreas has the ability to recover within one week based on histological appearance. In this study, we sought to characterize by RNA-sequencing (RNA-seq) the transcriptional profile of the recovering pancreas up to two weeks post-injury. We found that one week after injury there were 319 differentially expressed genes (DEGs) compared with baseline and that after two weeks there were 53 DEGs. Forty (12.5%) of the DEGs persisted from week one to week two, and another 13 DEGs newly emerged in the second week. Amongst the top up-regulated DEGs were several trypsinogen genes (trypsinogen 4, 5, 12, 15, and 16). To our knowledge, this is the first characterization of the transcriptome during pancreatic recovery by deep sequencing, and it reveals on a molecular basis that there is an ongoing recovery of the pancreas even after apparent histological resolution. The findings also raise the possibility of an emerging novel transcriptome upon pancreatic recovery.

DE-kupl: exhaustive capture of biological variation in RNA-seq data through k-mer decomposition

We introduce a k-mer-based computational protocol, DE-kupl, for capturing local RNA variation in a set of RNA-seq libraries, independently of a reference genome or transcriptome. DE-kupl extracts all k-mers with differential abundance directly from the raw data files. This enables the retrieval of virtually all variation present in an RNA-seq data set. This variation is subsequently assigned to biological events or entities such as differential long non-coding RNAs, splice and polyadenylation variants, introns, repeats, editing or mutation events, and exogenous RNA. Applying DE-kupl to human RNA-seq data sets identified multiple types of novel events, reproducibly across independent RNA-seq experiments.

RNA polymerase III limits longevity downstream of TORC1

Three distinct RNA polymerases (Pols) transcribe different classes of genes in the eukaryotic nucleus1. Pol III is the essential, evolutionarily conserved enzyme that generates short, non-coding RNAs, including transfer RNAs (tRNAs) and 5S ribosomal RNA (rRNA)2. Historical focus on transcription of protein-coding genes has left the roles of Pol III in organismal physiology relatively unexplored. The prominent regulator of Pol III activity, Target of Rapamycin kinase Complex 1 (TORC1), is an important longevity determinant3, raising the question of Pol III’s involvement in ageing. Here we show that Pol III limits lifespan downstream of TORC1. We find that a reduction in Pol III extends chronological lifespan in yeast and organismal lifespan in worms and flies. Inhibiting Pol III activity in the adult worm or fly gut is sufficient to extend lifespan, and in flies, longevity can be achieved by Pol III inhibition specifically in the intestinal stem cells (ISCs). The longevity phenotype is associated with amelioration of age-related gut pathology and functional decline, dampened protein synthesis and increased tolerance of proteostatic stress. Importantly, Pol III acts downstream of TORC1 for lifespan and limiting Pol III activity in the adult gut achieves the full longevity benefit of systemic TORC1 inhibition. Hence, Pol III is a pivotal output of this key nutrient signalling network for longevity; Pol III’s growth-promoting, anabolic activity mediates the acceleration of ageing by TORC1. The evolutionary conservation of Pol III affirms its potential as a therapeutic target.

The Long Noncoding RNA Landscape of the Ischemic Human Left Ventricle

BACKGROUND

The discovery of functional classes of long noncoding RNAs (lncRNAs) has expanded our understanding of the variety of RNA species that exist in cells. In the heart, lncRNAs have been implicated in the regulation of development, ischemic and dilated cardiomyopathy, and myocardial infarction. Nevertheless, there is a limited description of expression profiles for these transcripts in human subjects.

METHODS AND RESULTS

We obtained left ventricular tissue from human patients undergoing cardiac surgery and used RNA sequencing to describe an lncRNA profile. We then identified a list of lncRNAs that were differentially expressed between pairs of samples before and after the ischemic insult of cardiopulmonary bypass. The expression of some of these lncRNAs correlates with ischemic time. Coding genes in close proximity to differentially expressed lncRNAs and coding genes that have coordinated expression with these lncRNAs are enriched in functional categories related to myocardial infarction, including heart function, metabolism, the stress response, and the immune system.

CONCLUSIONS

We describe a list of lncRNAs that are differentially expressed after ischemia in the human heart. These genes are predicted to function in pathways consistent with myocardial injury. As a result, lncRNAs may serve as novel diagnostic and therapeutic targets for ischemic heart disease.

CLINICAL TRIAL REGISTRATION

URL: http://www.clinicaltrials.gov. Unique identifier: NCT00985049.

Androgen-responsive non-coding small RNAs extend the potential of HCG stimulation to act as a bioassay of androgen sufficiency

BACKGROUND

It is unclear whether a short-term change in circulating androgens is associated with changes in the transcriptome of the peripheral blood mononuclear cells (PBMC).

AIMS AND METHODS

To explore the effect of hCG stimulation on the PBMC transcriptome, 12 boys with a median age (range) of 0.7 years (0.3, 11.2) who received intramuscular hCG 1500u on 3 consecutive days as part of their investigations underwent transcriptomic array analysis on RNA extracted from peripheral blood mononuclear cells before and after hCG stimulation.

RESULTS

Median pre- and post-hCG testosterone for the overall group was 0.7 nmol/L (<0.5, 6) and 7.9 nmol/L (<0.5, 31.5), respectively. Of the 12 boys, 3 (25%) did not respond to hCG stimulation with a pre and post median serum testosterone of <0.5 nmol/L and <0.5 nmol/L, respectively. When corrected for gene expression changes in the non-responders to exclude hCG effects, all 9 of the hCG responders consistently demonstrated a 20% or greater increase in the expression of piR-37153 and piR-39248, non-coding PIWI-interacting RNAs (piRNAs). In addition, of the 9 responders, 8, 6 and 4 demonstrated a 30, 40 and 50% rise, respectively, in a total of 2 further piRNAs. In addition, 3 of the responders showed a 50% or greater rise in the expression of another small RNA, SNORD5. On comparing fold-change in serum testosterone with fold-change in the above transcripts, a positive correlation was detected for SNORD5 (P = 0.01).

CONCLUSIONS

The identification of a dynamic and androgen-responsive PBMC transcriptome extends the potential value of the hCG test for the assessment of androgen sufficiency.

Repertoire of noncoding RNAs in corpus luteum of early pregnancy in buffalo (Bubalus bubalis)

Aim:

The present study was designed to identify other noncoding RNAs (ncRNAs) in the corpus luteum (CL) during early pregnancy in buffalo.

Materials and Methods:

For this study, CL (n=2) from two buffalo gravid uteri, obtained from the slaughter house, was transported to laboratory after snap freezing in liquid nitrogen (−196°C). The stage of pregnancy was determined by measuring the crown-rump region of the fetus. This was followed by isolation of RNA and deep sequencing. Post-deep sequencing, the obtained reads were checked and aligned against various ncRNA databases (GtRNA, RFAM, and deep guide). Various parameters, namely, frequency of specific ncRNAs, length, mismatch, and genomic location target in several model species were deciphered.

Results:

Frequency of piwi-interacting RNAs (piwi-RNAs), having target location in rodents and human genomes, were significantly higher compared to other piwi-RNAs and ncRNAs. Ribosomal RNAs (rRNAs) deduced had nucleotides (nts) ranging from 17 to 50 nts, but the occurrence of small length rRNAs was more than lengthier fragments. The target on 16S rRNA species confirms the conservation of 16S rRNA across species. With respect to transfer RNA (tRNA), the abundantly occurring tRNAs were unique with no duplication. Small nucleolar RNAs (snoRNAs), identified in this study, showed a strong tendency for coding box C/D snoRNAs in comparison to H/ACA snoRNAs. Regulatory and evolutionary implications of these identified ncRNAs are yet to be delineated in many species, including buffaloes.

Conclusion:

This is the first report of identification of other ncRNAs in CL of early pregnancy in buffalo.

Abf1 and other general regulatory factors control ribosome biogenesis gene expression in budding yeast

Ribosome biogenesis in Saccharomyces cerevisiae involves a regulon of >200 genes (Ribi genes) coordinately regulated in response to nutrient availability and cellular growth rate. Two cis-acting elements called PAC and RRPE are known to mediate Ribi gene repression in response to nutritional downshift. Here, we show that most Ribi gene promoters also contain binding sites for one or more General Regulatory Factors (GRFs), most frequently Abf1 and Reb1, and that these factors are enriched in vivo at Ribi promoters. Abf1/Reb1/Tbf1 promoter association was required for full Ribi gene expression in rich medium and for its modulation in response to glucose starvation, characterized by a rapid drop followed by slow recovery. Such a response did not entail changes in Abf1 occupancy, but it was paralleled by a quick increase, followed by slow decrease, in Rpd3L histone deacetylase occupancy. Remarkably, Abf1 site disruption also abolished Rpd3L complex recruitment in response to starvation. Extensive mutational analysis of the DBP7 promoter revealed a complex interplay of Tbf1 sites, PAC and RRPE in the transcriptional regulation of this Ribi gene. Our observations point to GRFs as new multifaceted players in Ribi gene regulation both during exponential growth and under repressive conditions.

Small non-coding RNA and cancer

The ENCODE project has reported that at least 80% of the human genome is biologically active, yet only a small part of human DNA encodes for protein. The massive amount of RNA transcribed but not translated into protein can be classified as housekeeping RNA (such as rRNA, tRNA) and regulatory RNA (such as miRNA, piRNA, lncRNA). Small non-coding RNAs, in particular, have been the focus of many studies in the last 20 years and their fundamental role in many human diseases is currently well established. Inter alia, their role in cancer development and progression, as well as in drug resistance, is being increasingly investigated. In this review, focusing our attention on recent research results, we provide an overview of the four large classes of small non-coding RNAs, namely, miRNAs, piRNAs, snoRNA and the new class of tRNA-derived fragments, highlighting their fundamental role in cancer and their potential as diagnostic and prognostic biomarkers.

Protein coding genes as hosts for noncoding RNA expression

With the emergence of high-throughput sequence characterization methods and the subsequent improvements in gene annotations, it is becoming increasingly clear that a large proportion of eukaryotic protein-coding genes (as many as 50% in human) serve as host genes for non-coding RNA genes. Amongst the most extensively characterized embedded non-coding RNA genes, small nucleolar RNAs and microRNAs represent abundant families. Encoded individually or clustered, in sense or antisense orientation with respect to their host and independently expressed or dependent on host expression, the genomic characteristics of embedded genes determine their biogenesis and the extent of their relationship with their host gene. Not only can host genes and the embedded genes they harbour be co-regulated and mutually modulate each other, many are functionally coupled playing a role in the same cellular pathways. And while host-non-coding RNA relationships can be highly conserved, mechanisms have been identified, and in particular an association with transposable elements, allowing the appearance of copies of non-coding genes nested in host genes, or the migration of embedded genes from one host gene to another. The study of embedded non-coding genes and their relationship with their host genes increases the complexity of cellular networks and provides important new regulatory links that are essential to properly understand cell function.

Meta-analysis of the prognostic value of lncRNA ZFAS1 in patients with solid tumors

Expression of ZFAS1, a newly identified long noncoding RNA (lncRNA), is dysregulated in several types of cancer. Here we assessed the prognostic value of ZFAS1 in solid tumors. A comprehensive literature search was performed by screening the PubMed, EMBASE, MEDLINE, Cochrane Library, CNKI, and Wanfang databases. A total of 874 patients from 10 studies were included. The pooled analysis demonstrated that patients with high ZFAS1 expression had a significantly shorter overall survival (OS) (HR, 1.58; 95% CI, 1.28–1.97; P < 0.001) and recurrence-free survival (RFS) (HR, 1.90; 95% CI, 1.29–2.79; P = 0.001). Moreover, elevated ZFAS1 expression correlated with tumor size, tumor-node-metastasis (TNM) stage, and lymph node metastasis (LNM). These results demonstrate that increased ZFAS1 expression correlates with a poor prognosis in cancer patients, which suggests ZFAS1 might be useful as a potential prognostic biomarker in patients with solid tumors.

5' End Nicotinamide Adenine Dinucleotide Cap in Human Cells Promotes RNA Decay through DXO-Mediated deNADding

Eukaryotic mRNAs generally possess a 5' end N7 methyl guanosine (m⁷G) cap that promotes their translation and stability. However, mammalian mRNAs can also carry a 5' end nicotinamide adenine dinucleotide (NAD⁺) cap that, in contrast to the m⁷G cap, does not support translation but instead promotes mRNA decay. The mammalian and fungal noncanonical DXO/Rai1 decapping enzymes efficiently remove NAD⁺ caps, and cocrystal structures of DXO/Rai1 with 3'-NADP⁺ illuminate the molecular mechanism for how the "deNADding" reaction produces NAD⁺ and 5' phosphate RNA. Removal of DXO from cells increases NAD⁺-capped mRNA levels and enables detection of NAD⁺-capped intronic small nucleolar RNAs (snoRNAs), suggesting NAD⁺ caps can be added to 5'-processed termini. Our findings establish NAD⁺ as an alternative mammalian RNA cap and DXO as a deNADding enzyme modulating cellular levels of NAD⁺-capped RNAs. Collectively, these data reveal that mammalian RNAs can harbor a 5' end modification distinct from the classical m⁷G cap that promotes rather than inhibits RNA decay.

Phylogenetic distribution of plant snoRNA families

Background

Small nucleolar RNAs (snoRNAs) are one of the most ancient families amongst non-protein-coding RNAs. They are ubiquitous in Archaea and Eukarya but absent in bacteria. Their main function is to target chemical modifications of ribosomal RNAs. They fall into two classes, box C/D snoRNAs and box H/ACA snoRNAs, which are clearly distinguished by conserved sequence motifs and the type of chemical modification that they govern. Similarly to microRNAs, snoRNAs appear in distinct families of homologs that affect homologous targets. In animals, snoRNAs and their evolution have been studied in much detail. In plants, however, their evolution has attracted comparably little attention.

Results

In order to chart the phylogenetic distribution of individual snoRNA families in plants, we applied a sophisticated approach for identifying homologs of known plant snoRNAs across the plant kingdom. In response to the relatively fast evolution of snoRNAs, information on conserved sequence boxes, target sequences, and secondary structure is combined to identify additional snoRNAs. We identified 296 families of snoRNAs in 24 species and traced their evolution throughout the plant kingdom. Many of the plant snoRNA families comprise paralogs. We also found that targets are well-conserved for most snoRNA families.

Conclusions

The sequence conservation of snoRNAs is sufficient to establish homologies between phyla. The degree of this conservation tapers off, however, between land plants and algae. Plant snoRNAs are frequently organized in highly conserved spatial clusters. As a resource for further investigations we provide carefully curated and annotated alignments for each snoRNA family under investigation.

Electronic supplementary material

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

Plasticity of archaeal C/D box sRNA biogenesis

Archaeal and eukaryotic organisms contain sets of C/D box s(no)RNAs with guide sequences that determine ribose 2'-O-methylation sites of target RNAs. The composition of these C/D box sRNA sets is highly variable between organisms and results in varying RNA modification patterns which are important for ribosomal RNA folding and stability. Little is known about the genomic organization of C/D box sRNA genes in archaea. Here, we aimed to obtain first insights into the biogenesis of these archaeal C/D box sRNAs and analyzed the genetic context of more than 300 archaeal sRNA genes. We found that the majority of these genes do not possess independent promoters but are rather located at positions that allow for co-transcription with neighboring genes and their start or stop codons were frequently incorporated into the conserved boxC and D motifs. The biogenesis of plasmid-encoded C/D box sRNA variants was analyzed in vivo in Sulfolobus acidocaldarius. It was found that C/D box sRNA maturation occurs independent of their genetic context and relies solely on the presence of intact RNA kink-turn structures. The observed plasticity of C/D box sRNA biogenesis is suggested to enable their accelerated evolution and, consequently, allow for adjustments of the RNA modification landscape.

Identification of Non-Coding RNAs in the Candida parapsilosis Species Group

The Candida CTG clade is a monophyletic group of fungal species that translates CTG as serine, and includes the pathogens Candida albicans and Candida parapsilosis. Research has typically focused on identifying protein-coding genes in these species. Here, we use bioinformatic and experimental approaches to annotate known classes of non-coding RNAs in three CTG-clade species, Candida parapsilosis, Candida orthopsilosis and Lodderomyces elongisporus. We also update the annotation of ncRNAs in the C. albicans genome. The majority of ncRNAs identified were snoRNAs. Approximately 50% of snoRNAs (including most of the C/D box class) are encoded in introns. Most are within mono- and polycistronic transcripts with no protein coding potential. Five polycistronic clusters of snoRNAs are highly conserved in fungi. In polycistronic regions, splicing occurs via the classical pathway, as well as by nested and recursive splicing. We identified spliceosomal small nuclear RNAs, the telomerase RNA component, signal recognition particle, RNase P RNA component and the related RNase MRP RNA component in all three genomes. Stem loop IV of the U2 spliceosomal RNA and the associated binding proteins were lost from the ancestor of C. parapsilosis and C. orthopsilosis, following the divergence from L. elongisporus. The RNA component of the MRP is longer in C. parapsilosis, C. orthopsilosis and L. elongisporus than in S. cerevisiae, but is substantially shorter than in C. albicans.

Profiling of Small Nucleolar RNAs by Next Generation Sequencing: Potential New Players for Breast Cancer Prognosis

One of the most abundant, yet least explored, classes of RNA is the small nucleolar RNAs (snoRNAs), which are well known for their involvement in post-transcriptional modifications of other RNAs. Although snoRNAs were only considered to perform housekeeping functions for a long time, recent studies have highlighted their importance as regulators of gene expression and as diagnostic/prognostic markers. However, the prognostic potential of these RNAs has not been interrogated for breast cancer (BC). The objective of the current study was to identify snoRNAs as prognostic markers for BC. Small RNA sequencing (Illumina Genome Analyzer IIx) was performed for 104 BC cases and 11 normal breast tissues. Partek Genomics Suite was used for analyzing the sequencing files. Two independent and proven approaches were used to identify prognostic markers: case-control (CC) and case-only (CO). For both approaches, snoRNAs significant in the permutation test, following univariate Cox proportional hazards regression model were used for constructing risk scores. Risk scores were subsequently adjusted for potential confounders in a multivariate Cox model. For both approaches, thirteen snoRNAs were associated with overall survival and/or recurrence free survival. Patients belonging to the high-risk group were associated with poor outcomes, and the risk score was significant after adjusting for confounders. Validation of representative snoRNAs (SNORD46 and SNORD89) using qRT-PCR confirmed the observations from sequencing experiments. We also observed 64 snoRNAs harboring piwi-interacting RNAs and/or microRNAs that were predicted to target genes (mRNAs) involved in tumorigenesis. Our results demonstrate the potential of snoRNAs to serve (i) as novel prognostic markers for BC and (ii) as indirect regulators of gene expression.

Telomerase Regulation from Beginning to the End

The vast body of literature regarding human telomere maintenance is a true testament to the importance of understanding telomere regulation in both normal and diseased states. In this review, our goal was simple: tell the telomerase story from the biogenesis of its parts to its maturity as a complex and function at its site of action, emphasizing new developments and how they contribute to the foundational knowledge of telomerase and telomere biology.

Chemical Tools for Dissecting the Role of lncRNAs in Epigenetic Regulation

Proper control and maintenance of gene expression is critical for cellular identity and maintenance. Transcription of RNA from the genome is intimately controlled by post-translational chemical modification of histone tails and DNA. Recent studies have demonstrated that chromatin-remodeling complexes seek out their target genomic loci through the help of noncoding RNA molecules. Within this Review, we will outline how the use of biochemical techniques has shed light on the mechanisms employed by RNA to guide these complexes and therefore control gene expression.

Inhibition of 13-cis retinoic acid-induced gene expression of reactive-resistance genes by thalidomide in glioblastoma tumours in vivo

The cell differentiation potential of 13-cis retinoic acid (RA) has not succeeded in the clinical treatment of glioblastoma (GBM) so far. However, RA may also induce the expression of resistance genes such as HOXB7 which can be suppressed by Thalidomide (THAL). Therefore, we tested if combined treatment with RA+THAL may inhibit growth of glioblastoma in vivo. Treatment with RA+THAL but not RA or THAL alone significantly inhibited tumour growth. The synergistic effect of RA and THAL was corroborated by the effect on proliferation of glioblastoma cell lines in vitro. HOXB7 was not upregulated but microarray analysis validated by real-time PCR identified four potential resistance genes (IL-8, HILDPA, IGFBPA, and ANGPTL4) whose upregulation by RA was suppressed by THAL. Furthermore, genes coding for small nucleolar RNAs (snoRNA) were identified as a target for RA for the first time, and their upregulation was maintained after combined treatment. Pathway analysis showed upregulation of the Ribosome pathway and downregulation of pathways associated with proliferation and inflammation. In conclusion, combined treatment with RA + THAL delayed growth of GBM xenografts and suppressed putative resistance genes associated with hypoxia and angiogenesis. This encourages further pre-clinical and clinical studies of this drug combination in GBM.

tRid, an enabling method to isolate previously inaccessible small RNA fractions

Detection of rare small RNA species whose sizes are overlapping with tRNAs often suffers from insufficient sensitivity due to the overwhelming abundance of tRNAs. We here report a method, named tRid (tRNA rid), for removing abundant tRNAs from small RNA fractions regardless of tRNA sequence species. By means of tRid, we are able to selectively enrich small RNAs which have been previously difficult to access due to mass existence of tRNAs in such fractions. A flexible tRNA-acylation ribozyme, known as flexizyme, is a key tool where the total tRNAs are aminoacylated with N-biotinylated phenylalanine regardless of tRNA sequences, and therefore the biotin-tagged tRNAs could be readily removed from the small RNA fractions by the use of streptavidin-immobilized magnetic beads. Next generation sequencing of the isolated small RNA fraction revealed that small RNAs with less than 200nt were effectively enriched, allowing us to identify previously unknown small RNAs in HeLa and E. coli.

Promoter architecture and transcriptional regulation of Abf1-dependent ribosomal protein genes in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, ribosomal protein gene (RPG) promoters display binding sites for either Rap1 or Abf1 transcription factors. Unlike Rap1-associated promoters, the small cohort of Abf1-dependent RPGs (Abf1-RPGs) has not been extensively investigated. We show that RPL3, RPL4B, RPP1A, RPS22B and RPS28A/B share a common promoter architecture, with an Abf1 site upstream of a conserved element matching the sequence recognized by Fhl1, a transcription factor which together with Ifh1 orchestrates Rap1-associated RPG regulation. Abf1 and Fhl1 promoter association was confirmed by ChIP and/or gel retardation assays. Mutational analysis revealed a more severe requirement of Abf1 than Fhl1 binding sites for RPG transcription. In the case of RPS22B an unusual Tbf1 binding site promoted both RPS22B and intron-hosted SNR44 expression. Abf1-RPG down-regulation upon TOR pathway inhibition was much attenuated at defective mutant promoters unable to bind Abf1. TORC1 inactivation caused the expected reduction of Ifh1 occupancy at RPS22B and RPL3 promoters, but unexpectedly it entailed largely increased Abf1 association with Abf1-RPG promoters. We present evidence that Abf1 recruitment upon nutritional stress, also observed for representative ribosome biogenesis genes, favours RPG transcriptional rescue upon nutrient replenishment, thus pointing to nutrient-regulated Abf1 dynamics at promoters as a novel mechanism in ribosome biogenesis control.

Cajal bodies are linked to genome conformation

The mechanisms underlying nuclear body (NB) formation and their contribution to genome function are unknown. Here we examined the non-random positioning of Cajal bodies (CBs), major NBs involved in spliceosomal snRNP assembly and their role in genome organization. CBs are predominantly located at the periphery of chromosome territories at a multi-chromosome interface. Genome-wide chromosome conformation capture analysis (4C-seq) using CB-interacting loci revealed that CB-associated regions are enriched with highly expressed histone genes and U small nuclear or nucleolar RNA (sn/snoRNA) loci that form intra- and inter-chromosomal clusters. In particular, we observed a number of CB-dependent gene-positioning events on chromosome 1. RNAi-mediated disassembly of CBs disrupts the CB-targeting gene clusters and suppresses the expression of U sn/snoRNA and histone genes. This loss of spliceosomal snRNP production results in increased splicing noise, even in CB-distal regions. Therefore, we conclude that CBs contribute to genome organization with global effects on gene expression and RNA splicing fidelity.

Nuclear bodies can nucleate at sites of active transcription and are beneficial for efficient gene expression. Here, the authors show that Cajal bodies, a prominent type of nuclear body, contribute to genome organization with global effects on gene expression and RNA splicing fidelity.

Introns: The Functional Benefits of Introns in Genomes

The intron has been a big biological mystery since it was first discovered in several aspects. First, all of the completely sequenced eukaryotes harbor introns in the genomic structure, whereas no prokaryotes identified so far carry introns. Second, the amount of total introns varies in different species. Third, the length and number of introns vary in different genes, even within the same species genome. Fourth, all introns are copied into RNAs by transcription and DNAs by replication processes, but intron sequences do not participate in protein-coding sequences. The existence of introns in the genome should be a burden to some cells, because cells have to consume a great deal of energy to copy and excise them exactly at the correct positions with the help of complicated spliceosomal machineries. The existence throughout the long evolutionary history is explained, only if selective advantages of carrying introns are assumed to be given to cells to overcome the negative effect of introns. In that regard, we summarize previous research about the functional roles or benefits of introns. Additionally, several other studies strongly suggesting that introns should not be junk will be introduced.

Identifying transcription start sites and active enhancer elements using BruUV-seq

BruUV-seq utilizes UV light to introduce transcription-blocking DNA lesions randomly in the genome prior to bromouridine-labeling and deep sequencing of nascent RNA. By inhibiting transcription elongation, but not initiation, pre-treatment with UV light leads to a redistribution of transcription reads resulting in the enhancement of nascent RNA signal towards the 5′-end of genes promoting the identification of transcription start sites (TSSs). Furthermore, transcripts associated with arrested RNA polymerases are protected from 3′–5′ degradation and thus, unstable transcripts such as putative enhancer RNA (eRNA) are dramatically increased. Validation of BruUV-seq against GRO-cap that identifies capped run-on transcripts showed that most BruUV-seq peaks overlapped with GRO-cap signal over both TSSs and enhancer elements. Finally, BruUV-seq identified putative enhancer elements induced by tumor necrosis factor (TNF) treatment concomitant with expression of nearby TNF-induced genes. Taken together, BruUV-seq is a powerful new approach for identifying TSSs and active enhancer elements genome-wide in intact cells.