Identifying key underlying regulatory networks and predicting targets of orphan C/D box SNORD116 snoRNAs in Prader-Willi syndrome

Prader-Willi syndrome (PWS) is a rare neurodevelopmental disorder characterized principally by initial symptoms of neonatal hypotonia and failure-to-thrive in infancy, followed by hyperphagia and obesity. It is well established that PWS is caused by loss of paternal expression of the imprinted region on chromosome 15q11-q13. While most PWS cases exhibit megabase-scale deletions of the paternal chromosome 15q11-q13 allele, several PWS patients have been identified harboring a much smaller deletion encompassing primarily SNORD116. This finding suggests SNORD116 is a direct driver of PWS phenotypes. The SNORD116 gene cluster is composed of 30 copies of individual SNORD116 C/D box small nucleolar RNAs (snoRNAs). Many C/D box snoRNAs have been shown to guide chemical modifications of other RNA molecules, often ribosomal RNA (rRNA). However, SNORD116 snoRNAs are termed 'orphans' because no verified targets have been identified and their sequences show no significant complementarity to rRNA. It is crucial to identify the targets and functions of SNORD116 snoRNAs because all reported PWS cases lack their expression. To address this, we engineered two different deletions modelling PWS in two distinct human embryonic stem cell (hESC) lines to control for effects of genetic background. Utilizing an inducible expression system enabled quick, reproducible differentiation of these lines into neurons. Systematic comparisons of neuronal gene expression across deletion types and genetic backgrounds revealed a novel list of 42 consistently dysregulated genes. Employing the recently described computational tool snoGloBe, we discovered these dysregulated genes are significantly enriched for predicted SNORD116 targeting versus multiple control analyses. Importantly, our results showed it is critical to use multiple isogenic cell line pairs, as this eliminated many spuriously differentially expressed genes. Our results indicate a novel gene regulatory network controlled by SNORD116 is likely perturbed in PWS patients.

TET3 epigenetically controls feeding and stress response behaviors via AGRP neurons

The TET family of dioxygenases promote DNA demethylation by oxidizing 5-methylcytosine to 5-hydroxymethylcytosine (5hmC). Hypothalamic agouti-related peptide-expressing (AGRP-expressing) neurons play an essential role in driving feeding, while also modulating nonfeeding behaviors. Besides AGRP, these neurons produce neuropeptide Y (NPY) and the neurotransmitter GABA, which act in concert to stimulate food intake and decrease energy expenditure. Notably, AGRP, NPY, and GABA can also elicit anxiolytic effects. Here, we report that in adult mouse AGRP neurons, CRISPR-mediated genetic ablation of Tet3, not previously known to be involved in central control of appetite and metabolism, induced hyperphagia, obesity, and diabetes, in addition to a reduction of stress-like behaviors. TET3 deficiency activated AGRP neurons, simultaneously upregulated the expression of Agrp, Npy, and the vesicular GABA transporter Slc32a1, and impeded leptin signaling. In particular, we uncovered a dynamic association of TET3 with the Agrp promoter in response to leptin signaling, which induced 5hmC modification that was associated with a chromatin-modifying complex leading to transcription inhibition, and this regulation occurred in both the mouse models and human cells. Our results unmasked TET3 as a critical central regulator of appetite and energy metabolism and revealed its unexpected dual role in the control of feeding and other complex behaviors through AGRP neurons.

Transcriptome-Wide Identification of 2'-O-Methylation Sites with RibOxi-Seq

The ability to detect 2'-O-methylation sites (Nm) in high-throughput fashion is important, as increasing evidence points to a more diverse landscape for this RNA modification as well as the possibility of yet unidentified functions. Here we describe an optimized version of RibOxi-seq, which is built upon the original published method, that not only accurately profiles ribosomal RNA (rRNA) Nm sites with minimal RNA input but is also robust enough to identify mRNA intronic and exonic sites.

Prader-Willi syndrome: reflections on seminal studies and future therapies

Prader-Willi syndrome (PWS) is caused by the loss of function of the paternally inherited 15q11-q13 locus. This region is governed by genomic imprinting, a phenomenon in which genes are expressed exclusively from one parental allele. The genomic imprinting of the 15q11-q13 locus is established in the germline and is largely controlled by a bipartite imprinting centre. One part, termed the Prader-Willi syndrome imprinting center (PWS-IC), comprises a CpG island that is unmethylated on the paternal allele and methylated on the maternal allele. The second part, termed the Angelman syndrome imprinting centre, is required to silence the PWS_IC in the maternal germline. The loss of the paternal contribution of the imprinted 15q11-q13 locus most frequently occurs owing to a large deletion of the entire imprinted region but can also occur through maternal uniparental disomy or an imprinting defect. While PWS is considered a contiguous gene syndrome based on large-deletion and uniparental disomy patients, the lack of expression of only non-coding RNA transcripts from the SNURF-SNRPN/SNHG14 may be the primary cause of PWS. Patients with small atypical deletions of the paternal SNORD116 cluster alone appear to have most of the PWS related clinical phenotypes. The loss of the maternal contribution of the 15q11-q13 locus causes a separate and distinct condition called Angelman syndrome. Importantly, while much has been learned about the regulation and expression of genes and transcripts deriving from the 15q11-q13 locus, there remains much to be learned about how these genes and transcripts contribute at the molecular level to the clinical traits and developmental aspects of PWS that have been observed.

The large repertoire of 2'-O-methylation guided by C/D snoRNAs on Trypanosoma brucei rRNA

The parasite Trypanosoma brucei cycles between insect and mammalian hosts, and is the causative agent of sleeping sickness. Here, we performed genome-wide mapping of 2'-O-methylations (Nms) on trypanosome rRNA using three high-throughput sequencing methods; RibOxi-seq, RiboMeth-seq and 2'-OMe-seq. This is the first study using three genome-wide mapping approaches on rRNA from the same species showing the discrepancy among the methods. RibOxi-seq detects all the sites, but RiboMeth-seq is the only method to evaluate the level of a single Nm site. The sequencing revealed at least ninety-nine Nms guided by eighty-five snoRNAs among these thirty-eight Nms are trypanosome specific sites. We present the sequence and target of the C/D snoRNAs guiding on rRNA. This is the highest number of Nms detected to date on rRNA of a single cell parasite. Based on RiboMeth-seq, several Nm sites were found to be differentially regulated at the two stages of the parasite life cycle, the insect procyclic form (PCF) versus the bloodstream form (BSF) in the mammalian host.

Distinct Processing of lncRNAs Contributes to Non-conserved Functions in Stem Cells

Long noncoding RNAs (lncRNAs) evolve more rapidly than mRNAs. Whether conserved lncRNAs undergo conserved processing, localization, and function remains unexplored. We report differing subcellular localization of lncRNAs in human and mouse embryonic stem cells (ESCs). A significantly higher fraction of lncRNAs is localized in the cytoplasm of hESCs than in mESCs. This turns out to be important for hESC pluripotency. FAST is a positionally conserved lncRNA but is not conserved in its processing and localization. In hESCs, cytoplasm-localized hFAST binds to the WD40 domain of the E3 ubiquitin ligase β-TrCP and blocks its interaction with phosphorylated β-catenin to prevent degradation, leading to activated WNT signaling, required for pluripotency. In contrast, mFast is nuclear retained in mESCs, and its processing is suppressed by the splicing factor PPIE, which is highly expressed in mESCs but not hESCs. These findings reveal that lncRNA processing and localization are previously under-appreciated contributors to the rapid evolution of function.

Dynamic Imaging of RNA in Living Cells by CRISPR-Cas13 Systems

Visualizing the location and dynamics of RNAs in live cells is key to understanding their function. Here, we identify two endonuclease-deficient, single-component programmable RNA-guided and RNA-targeting Cas13 RNases (dCas13s) that allow robust real-time imaging and tracking of RNAs in live cells, even when using single 20- to 27-nt-long guide RNAs. Compared to the aptamer-based MS2-MCP strategy, an optimized dCas13 system is user friendly, does not require genetic manipulation, and achieves comparable RNA-labeling efficiency. We demonstrate that the dCas13 system is capable of labeling NEAT1, SatIII, MUC4, and GCN4 RNAs and allows the study of paraspeckle-associated NEAT1 dynamics. Applying orthogonal dCas13 proteins or combining dCas13 and MS2-MCP allows dual-color imaging of RNAs in single cells. Further combination of dCas13 and dCas9 systems allows simultaneous visualization of genomic DNA and RNA transcripts in living cells.

CIRCexplorer3: A CLEAR Pipeline for Direct Comparison of Circular and Linear RNA Expression

Sequences of circular RNAs (circRNAs) produced from back-splicing of exon(s) completely overlap with those from cognate linear RNAs transcribed from the same gene loci with the exception of their back-splicing junction (BSJ) sites. Therefore, examination of global circRNA expression from RNA-seq datasets generally relies on the detection of RNA-seq fragments spanning BSJ sites, which is different from the quantification of linear RNA expression by normalized RNA-seq fragments mapped to whole gene bodies. Thus, direct comparison of circular and linear RNA expression from the same gene loci in a genome-wide manner has remained challenging. Here, we update the previously-reported CIRCexplorer pipeline to version 3 for circular and linear RNA expression analysis from ribosomal-RNA depleted RNA-seq (CIRCexplorer3-CLEAR). A new quantitation parameter, fragments per billion mapped bases (FPB), is applied to evaluate circular and linear RNA expression individually by fragments mapped to circRNA-specific BSJ sites or to linear RNA-specific splicing junction (SJ) sites. Comparison of circular and linear RNA expression levels is directly achieved by dividing FPBcirc by FPBlinear to generate a CIRCscore, which indicates the relative circRNA expression level using linear RNA expression level as the background. Highly-expressed circRNAs with low cognate linear RNA expression background can be readily identified by CIRCexplorer3-CLEAR for further investigation. CIRCexplorer3-CLEAR is publically available at https://github.com/YangLab/CLEAR.

H19 lncRNA Promotes Skeletal Muscle Insulin Sensitivity in Part by Targeting AMPK

Skeletal muscle plays a pivotal role in regulating systemic glucose homeostasis in part through the conserved cellular energy sensor AMPK. AMPK activation increases glucose uptake, lipid oxidation, and mitochondrial biogenesis, leading to enhanced muscle insulin sensitivity and whole-body energy metabolism. Here we show that the muscle-enriched H19 long noncoding RNA (lncRNA) acts to enhance muscle insulin sensitivity, at least in part, by activating AMPK. We identify the atypical dual-specificity phosphatase DUSP27/DUPD1 as a potentially important downstream effector of H19. We show that DUSP27, which is highly expressed in muscle with previously unknown physiological function, interacts with and activates AMPK in muscle cells. Consistent with decreased H19 expression in the muscle of insulin-resistant human subjects and rodents, mice with genetic H19 ablation exhibit muscle insulin resistance. Furthermore, a high-fat diet downregulates muscle H19 via both posttranscriptional and epigenetic mechanisms. Our results uncover an evolutionarily conserved, highly expressed lncRNA as an important regulator of muscle insulin sensitivity.

Elevated hepatic expression of H19 long noncoding RNA contributes to diabetic hyperglycemia

Excessive hepatic glucose production (HGP) contributes significantly to the hyperglycemia of type 2 diabetes; however, the molecular mechanism underlying this dysregulation remains poorly understood. Here, we show that fasting temporally increases the expression of H19 long noncoding RNA (lncRNA) in nondiabetic mouse liver, whereas its level is chronically elevated in diet-induced diabetic mice, consistent with the previously reported chronic hepatic H19 increase in diabetic patients. Importantly, liver-specific H19 overexpression promotes HGP, hyperglycemia, and insulin resistance, while H19 depletion enhances insulin-dependent suppression of HGP. Using genome-wide methylation and transcriptome analyses, we demonstrate that H19 knockdown in hepatic cells alters promoter methylation and expression of Hnf4a, a master gluconeogenic transcription factor, and that this regulation is recapitulated in vivo. Our findings offer a mechanistic explanation of lncRNA H19's role in the pathogenesis of diabetic hyperglycemia and suggest that targeting hepatic H19 may hold the potential of new treatment for this disease.

An Ultraconserved Element (UCE) controls homeostatic splicing of ARGLU1 mRNA

Arginine and Glutamate-Rich protein 1 (ARGLU1) is a protein whose function is poorly understood, but may act in both transcription and pre-mRNA splicing. We demonstrate that the ARGLU1 gene expresses at least three distinct RNA splice isoforms – a fully spliced isoform coding for the protein, an isoform containing a retained intron that is detained in the nucleus, and an isoform containing an alternative exon that targets the transcript for nonsense mediated decay. Furthermore, ARGLU1 contains a long, highly evolutionarily conserved sequence known as an Ultraconserved Element (UCE) that is within the retained intron and overlaps the alternative exon. Manipulation of the UCE, in a reporter minigene or via random mutations in the genomic context using CRISPR/Cas9, changed the splicing pattern. Further, overexpression of the ARGLU1 protein shifted the splicing of endogenous ARGLU1 mRNA, resulting in an increase in the retained intron isoform and nonsense mediated decay susceptible isoform and a decrease in the fully spliced isoform. Taken together with data showing that functional protein knockout shifts splicing toward the fully spliced isoform, our data are consistent with a model in which unproductive splicing complexes assembled at the alternative exon lead to inefficient splicing and intron retention.

High-throughput and site-specific identification of 2'-O-methylation sites using ribose oxidation sequencing (RibOxi-seq)

Ribose methylation (2'-O-methylation, 2'-OMe) occurs at high frequencies in rRNAs and other small RNAs and is carried out using a shared mechanism across eukaryotes and archaea. As RNA modifications are important for ribosome maturation, and alterations in these modifications are associated with cellular defects and diseases, it is important to characterize the landscape of 2'-O-methylation. Here we report the development of a highly sensitive and accurate method for ribose methylation detection using next-generation sequencing. A key feature of this method is the generation of RNA fragments with random 3'-ends, followed by periodate oxidation of all molecules terminating in 2',3'-OH groups. This allows only RNAs harboring 2'-OMe groups at their 3'-ends to be sequenced. Although currently requiring microgram amounts of starting material, this method is robust for the analysis of rRNAs even at low sequencing depth.

Metformin alters DNA methylation genome-wide via the H19/SAHH axis

The molecular mechanisms underlying the antineoplastic properties of metformin, a first-line drug for type 2 diabetes, remain elusive. Here we report that metformin induces genome-wide alterations in DNA methylation by modulating the activity of S-adenosylhomocysteine hydrolase (SAHH). Exposing cancer cells to metformin leads to hypermethylation of tumor-promoting pathway genes and concomitant inhibition of cell proliferation. Metformin acts by upregulating microRNA let-7 through AMPK activation, leading to degradation of H19 long noncoding RNA, which normally binds to and inactivates SAHH. H19 knockdown activates SAHH, enabling DNA methyltransferase 3B to methylate a subset of genes. This metformin-induced H19 repression and alteration of gene methylation are recapitulated in endometrial cancer tissue samples obtained from patients treated with antidiabetic doses of metformin. Our findings unveil a novel mechanism of action for the drug metformin with implications for the molecular basis of epigenetic dysregulation in cancer. This novel mechanism of action also may be occurring in normal cells.

H19 lncRNA alters DNA methylation genome wide by regulating S-adenosylhomocysteine hydrolase

DNA methylation is essential for mammalian development and physiology. Here we report that the developmentally regulated H19 lncRNA binds to and inhibits S-adenosylhomocysteine hydrolase (SAHH), the only mammalian enzyme capable of hydrolysing S-adenosylhomocysteine (SAH). SAH is a potent feedback inhibitor of S-adenosylmethionine (SAM)-dependent methyltransferases that methylate diverse cellular components, including DNA, RNA, proteins, lipids and neurotransmitters. We show that H19 knockdown activates SAHH, leading to increased DNMT3B-mediated methylation of an lncRNA-encoding gene Nctc1 within the Igf2-H19-Nctc1 locus. Genome-wide methylation profiling reveals methylation changes at numerous gene loci consistent with SAHH modulation by H19. Our results uncover an unanticipated regulatory circuit involving broad epigenetic alterations by a single abundantly expressed lncRNA that may underlie gene methylation dynamics of development and diseases and suggest that this mode of regulation may extend to other cellular components.

Global Analysis of Mouse Polyomavirus Infection Reveals Dynamic Regulation of Viral and Host Gene Expression and Promiscuous Viral RNA Editing

Mouse polyomavirus (MPyV) lytically infects mouse cells, transforms rat cells in culture, and is highly oncogenic in rodents. We have used deep sequencing to follow MPyV infection of mouse NIH3T6 cells at various times after infection and analyzed both the viral and cellular transcriptomes. Alignment of sequencing reads to the viral genome illustrated the transcriptional profile of the early-to-late switch with both early-strand and late-strand RNAs being transcribed at all time points. A number of novel insights into viral gene expression emerged from these studies, including the demonstration of widespread RNA editing of viral transcripts at late times in infection. By late times in infection, 359 host genes were seen to be significantly upregulated and 857 were downregulated. Gene ontology analysis indicated transcripts involved in translation, metabolism, RNA processing, DNA methylation, and protein turnover were upregulated while transcripts involved in extracellular adhesion, cytoskeleton, zinc finger binding, SH3 domain, and GTPase activation were downregulated. The levels of a number of long noncoding RNAs were also altered. The long noncoding RNA MALAT1, which is involved in splicing speckles and used as a marker in many late-stage cancers, was noticeably downregulated, while several other abundant noncoding RNAs were strongly upregulated. We discuss these results in light of what is currently known about the MPyV life cycle and its effects on host cell growth and metabolism.

Mouse polyomavirus (MPyV) is a small 5.3kb circular double-stranded DNA virus capable of causing tumors in a variety of tissues in immunocompromised mice. It has been a subject of study for over 60 years, yielding insights into a number of processes including tumorigenesis, cell cycle signaling, and transformation. This study serves to provide a global view of the MPyV infection by utilizing Illumina sequencing to observe changes in total RNA from both the virus and the host cell as well as applying new methods to more directly confirm the extent of A-to-I editing of viral RNA by host ADAR enzymes. This allows for a simultaneous observation of both host and viral transcriptional changes that occur as a result of early gene expression and the viral switch from early to late genes that occurs coincident with the initiation of DNA replication.

Attenuated Innate Immunity in Embryonic Stem Cells and Its Implications in Developmental Biology and Regenerative Medicine

Embryonic stem cells (ESCs) represent a promising cell source for regenerative medicine. Intensive research over the past 2 decades has led to the feasibility of using ESC-differentiated cells (ESC-DCs) in regenerative medicine. However, increasing evidence indicates that ESC-DCs generated by current differentiation methods may not have equivalent cellular functions to their in vivo counterparts. Recent studies have revealed that both human and mouse ESCs as well as some types of ESC-DCs lack or have attenuated innate immune responses to a wide range of infectious agents. These findings raise important concerns for their therapeutic applications since ESC-DCs, when implanted to a wound site of a patient, where they would likely be exposed to pathogens and inflammatory cytokines. Understanding whether an attenuated immune response is beneficial or harmful to the interaction between host and grafted cells becomes an important issue for ESC-based therapy. A substantial amount of recent evidence has demonstrated that the lack of innate antiviral responses is a common feature to ESCs and other types of pluripotent cells. This has led to the hypothesis that mammals may have adapted different antiviral mechanisms at different stages of organismal development. The underdeveloped innate immunity represents a unique and uncharacterized property of ESCs that may have important implications in developmental biology, immunology, and in regenerative medicine.

The H19/let-7 double-negative feedback loop contributes to glucose metabolism in muscle cells

The H19 lncRNA has been implicated in development and growth control and is associated with human genetic disorders and cancer. Acting as a molecular sponge, H19 inhibits microRNA (miRNA) let-7. Here we report that H19 is significantly decreased in muscle of human subjects with type-2 diabetes and insulin resistant rodents. This decrease leads to increased bioavailability of let-7, causing diminished expression of let-7 targets, which is recapitulated in vitro where H19 depletion results in impaired insulin signaling and decreased glucose uptake. Furthermore, acute hyperinsulinemia downregulates H19, a phenomenon that occurs through PI3K/AKT-dependent phosphorylation of the miRNA processing factor KSRP, which promotes biogenesis of let-7 and its mediated H19 destabilization. Our results reveal a previously undescribed double-negative feedback loop between sponge lncRNA and target miRNA that contributes to glucose regulation in muscle cells.

SnoVectors for nuclear expression of RNA

Many long noncoding RNAs (lncRNAs) are constrained to the nucleus to exert their functions. However, commonly used vectors that were designed to express mRNAs have not been optimized for the study of nuclear RNAs. We reported recently that sno-lncRNAs are not capped or polyadenylated but rather are terminated on each end by snoRNAs and their associated proteins. These RNAs are processed from introns and are strictly confined to the nucleus. Here we have used these features to design expression vectors that can stably express virtually any sequence of interest and constrain its accumulation to the nucleus. Further, these RNAs appear to retain normal nuclear associations and function. SnoVectors should be useful in conditions where nuclear RNA function is studied or where export to the cytoplasm needs to be avoided.

Long noncoding RNAs in imprinting and X chromosome inactivation

The field of long noncoding RNA (lncRNA) research has been rapidly advancing in recent years. Technological advancements and deep-sequencing of the transcriptome have facilitated the identification of numerous new lncRNAs, many with unusual properties, however, the function of most of these molecules is still largely unknown. Some evidence suggests that several of these lncRNAs may regulate their own transcription in cis, and that of nearby genes, by recruiting remodeling factors to local chromatin. Notably, lncRNAs are known to exist at many imprinted gene clusters. Genomic imprinting is a complex and highly regulated process resulting in the monoallelic silencing of certain genes, based on the parent-of-origin of the allele. It is thought that lncRNAs may regulate many imprinted loci, however, the mechanism by which they exert such influence is poorly understood. This review will discuss what is known about the lncRNAs of major imprinted loci, and the roles they play in the regulation of imprinting.

Innate immunity in pluripotent human cells: attenuated response to interferon-β

Type I interferon (IFN-α/β) binds to cell surface receptors IFNAR1 and IFNAR2 and triggers a signaling cascade that leads to the transcription of hundreds of IFN-stimulated genes. This response is a crucial component in innate immunity in that it establishes an "antiviral state" in cells and protects them against further damage. Previous work demonstrated that, compared with their differentiated counterparts, pluripotent human cells have a much weaker response to cytoplasmic double-stranded RNA (dsRNA) and are only able to produce a minimal amount of IFN-β. We show here that human embryonic stem cells (hESCs) and human induced pluripotent stem cells (hiPSCs) also exhibit an attenuated response to IFN-β. Even though all known type I IFN signaling components are expressed in these cells, STAT1 phosphorylation is greatly diminished upon IFN-β treatment. This attenuated response correlates with a high expression of suppressor of cytokine signaling 1 (SOCS1). Upon differentiation of hESCs into trophoblasts, cells acquire the ability to respond to IFN-β, and this is accompanied by a significant induction of STAT1 phosphorylation as well as a decrease in SOCS1 expression. Furthermore, SOCS1 knockdown in hiPSCs enhances their ability to respond to IFN-β. Taken together, our results suggest that an attenuated cellular response to type I IFNs may be a general feature of pluripotent human cells and that this is associated with high expression of SOCS1.

Long noncoding RNAs with snoRNA ends

We describe the discovery of sno-lncRNAs, a class of nuclear-enriched intron-derived long noncoding RNAs (lncRNAs) that are processed on both ends by the snoRNA machinery. During exonucleolytic trimming, the sequences between the snoRNAs are not degraded, leading to the accumulation of lncRNAs flanked by snoRNA sequences but lacking 5' caps and 3' poly(A) tails. Such RNAs are widely expressed in cells and tissues and can be produced by either box C/D or box H/ACA snoRNAs. Importantly, the genomic region encoding one abundant class of sno-lncRNAs (15q11-q13) is specifically deleted in Prader-Willi Syndrome (PWS). The PWS region sno-lncRNAs do not colocalize with nucleoli or Cajal bodies, but rather accumulate near their sites of synthesis. These sno-lncRNAs associate strongly with Fox family splicing regulators and alter patterns of splicing. These results thus implicate a previously unannotated class of lncRNAs in the molecular pathogenesis of PWS.

Genome-wide studies reveal that Lin28 enhances the translation of genes important for growth and survival of human embryonic stem cells

Lin28 inhibits the expression of let-7 microRNAs but also exhibits let-7-independent functions. Using immunoprecipitation and deep sequencing, we show here that Lin28 preferentially associates with a small subset of cellular mRNAs. Of particular interest are those for ribosomal proteins and metabolic enzymes, the expression levels of which are known to be coupled to cell growth and survival. Polysome profiling and reporter analyses suggest that Lin28 stimulates the translation of many or most of these targets. Moreover, Lin28-responsive elements were found within the coding regions of all target genes tested. Finally, a mutant Lin28 that still binds RNA but fails to interact with RNA helicase A (RHA), acts as a dominant-negative inhibitor of Lin28-dependent stimulation of translation. We suggest that Lin28, working in concert with RHA, enhances the translation of genes important for the growth and survival of human embryonic stem cells.

Genomewide characterization of non-polyadenylated RNAs

BACKGROUND

RNAs can be physically classified into poly(A)+ or poly(A)- transcripts according to the presence or absence of a poly(A) tail at their 3' ends. Current deep sequencing approaches largely depend on the enrichment of transcripts with a poly(A) tail, and therefore offer little insight into the nature and expression of transcripts that lack poly(A) tails.

RESULTS

We have used deep sequencing to explore the repertoire of both poly(A)+ and poly(A)- RNAs from HeLa cells and H9 human embryonic stem cells (hESCs). Using stringent criteria, we found that while the majority of transcripts are poly(A)+, a significant portion of transcripts are either poly(A)- or bimorphic, being found in both the poly(A)+ and poly(A)- populations. Further analyses revealed that many mRNAs may not contain classical long poly(A) tails and such messages are overrepresented in specific functional categories. In addition, we surprisingly found that a few excised introns accumulate in cells and thus constitute a new class of non-polyadenylated long non-coding RNAs. Finally, we have identified a specific subset of poly(A)- histone mRNAs, including two histone H1 variants, that are expressed in undifferentiated hESCs and are rapidly diminished upon differentiation; further, these same histone genes are induced upon reprogramming of fibroblasts to induced pluripotent stem cells.

CONCLUSIONS

We offer a rich source of data that allows a deeper exploration of the poly(A)- landscape of the eukaryotic transcriptome. The approach we present here also applies to the analysis of the poly(A)- transcriptomes of other organisms.

Nuclear Editing of mRNA 3'-UTRs

Hundreds of human genes express mRNAs that contain inverted repeat sequences within their 3'-UTRs. When expressed, these sequences can be promiscuously edited by ADAR enzymes, leading to the retention of mRNAs in nuclear paraspeckles. Here we discuss how this retention system can be used to regulate gene expression.

Molecular basis for an attenuated cytoplasmic dsRNA response in human embryonic stem cells

The introduction of double stranded RNA (dsRNA) into the cytoplasm of mammalian cells usually leads to a potent antiviral response resulting in the rapid induction of interferon beta (IFNβ). This response can be mediated by a number of dsRNA sensors, including TLR3, MDA5, RIG-I and PKR. We show here that pluripotent human cells (human embryonic stem (hES) cells and induced pluripotent (iPS) cells) do not induce interferon in response to cytoplasmic dsRNA, and we have used a variety of approaches to learn the underlying basis for this phenomenon. Two major cytoplasmic dsRNA sensors, TLR3 and MDA5, are not expressed in hES cells and iPS cells. PKR is expressed in hES cells, but is not activated by transfected dsRNA. In addition, RIG-I is expressed, but fails to respond to dsRNA because its signaling adapter, MITA/STING, is not expressed. Finally, the interferon-inducible RNAse L and oligoadenylate synthetase enzymes are also expressed at very low levels. Upon differentiation of hES cells into trophoblasts, cells acquire the ability to respond to dsRNA and this correlates with a significant induction of expression of TLR3 and its adaptor protein TICAM-1/TRIF. Taken together, our results reveal that the lack of an interferon response may be a general characteristic of pluripotency and that this results from the systematic downregulation of a number of genes involved in cytoplasmic dsRNA signaling.

Long noncoding RNAs in mammalian cells: what, where, and why?

Not all long, polyadenylated cellular RNAs encode polypeptides. In recent years, it has become apparent that a number of organisms express abundant amounts of transcripts that lack open reading frames or that are retained in the nucleus. Rather than accumulating silently in the cell, we now know that many of these long noncoding RNAs (lncRNAs) play important roles in nuclear architecture or in the regulation of gene expression. Here, we discuss some recent progress in our understanding of the functions of a number of important lncRNAs in mammalian cells.

Decoding the function of nuclear long non-coding RNAs

Long non-coding RNAs (lncRNAs) are mRNA-like, non-protein-coding RNAs that are pervasively transcribed throughout eukaryotic genomes. Rather than silently accumulating in the nucleus, many of these are now known or suspected to play important roles in nuclear architecture or in the regulation of gene expression. In this review, we highlight some recent progress in how lncRNAs regulate these important nuclear processes at the molecular level.

Altered nuclear retention of mRNAs containing inverted repeats in human embryonic stem cells: functional role of a nuclear noncoding RNA

In many cells, mRNAs containing inverted repeats (Alu repeats in humans) in their 3' untranslated regions (3'UTRs) are inefficiently exported to the cytoplasm. Nuclear retention correlates with adenosine-to-inosine editing and is in paraspeckle-associated complexes containing the proteins p54(nrb), PSF, and PSP1 alpha. We report that robust editing activity in human embryonic stem cells (hESCs) does not lead to nuclear retention. p54(nrb), PSF, and PSP1 alpha are all expressed in hESCs, but paraspeckles are absent and only appear upon differentiation. Paraspeckle assembly and function depend on expression of a long nuclear-retained noncoding RNA, NEAT1. This RNA is not detectable in hESCs but is induced upon differentiation. Knockdown of NEAT1 in HeLa cells results both in loss of paraspeckles and in enhanced nucleocytoplasmic export of mRNAs containing inverted Alu repeats. Taken together, these results assign a biological function to a large noncoding nuclear RNA in the regulation of mRNA export.

RNA processing in the polyoma virus life cycle

Not only is gene regulation in polyoma interesting, but it has also proven to be highly informative and illustrative of a number of novel concepts in gene regulation. Of special interest and importance are the mechanisms by which this virus switches from the expression of early gene products to late gene products after the onset of viral DNA replication. This switch is mediated at least in part by changes in transcription elongation and polyadenylation in the late region, and by the formation and editing of dsRNA in the nucleus. In this review we will summarize the regulation of RNA synthesis and processing during polyoma infection, and will point out in particular those aspects that have been most novel.

Gene regulation by sense-antisense overlap of polyadenylation signals

We show here that expression of genes from convergent transcription units can be regulated by the formation of double-stranded RNA (dsRNA) in the region of overlapping polyadenylation signals. The model system employed is the mouse polyomavirus. The early and late genes of polyomavirus are transcribed from opposite strands of the circular viral genome. At early times after infection, the early genes are expressed predominantly. Late gene expression increases dramatically upon the onset of DNA replication, when a major defect in polyadenylation of the late primary transcripts generates multigenomic RNAs that are precursors to the mature late mRNAs. Embedded in these late pre-mRNAs are sequences complementary to the early RNAs that act to down-regulate early gene expression via A-to-I editing of dsRNAs. In this system, the defective polyadenylation, and consequently the production of multigenomic late RNAs, depends on the context, and perhaps also, on the A-to-I editing of the poly(A) signal that overlaps the 3'-end of early transcripts.

Gene regulation by SINES and inosines: biological consequences of A-to-I editing of Alu element inverted repeats

The Alu elements are conserved approximately 300 nucleotide long repeat sequences that belong to the SINE family of retrotransposons found abundantly in primate genomes. Although the vast majority of Alu elements appear to be genetically inert, it has been tempting to consider the great majority of them as "junk DNA." However, a growing line of evidence suggests that transcribed Alu RNAs are in fact functionally involved in a number of diverse biological processes. Pairs of inverted Alu repeats in RNA can form duplex structures that lead to A-to-I editing by the ADAR enzymes. In this review we discuss the possible biological effects of Alu editing, with particular focus on the regulation of gene expression by inverted Alu repeats in the 3'-UTR regions of mRNAs.

On the mechanism of induction of heterochromatin by the RNA-binding protein vigilin

Vigilin is an RNA-binding protein localized to both the cytoplasm and the nucleus and has been previously implicated in heterochromatin formation and chromosome segregation. We demonstrate here that the C-terminal domain of human vigilin binds to the histone methyltransferase SUV39H1 in vivo. This association is independent of RNA and maps to a site on vigilin that is not involved in its interaction with several other known protein partners. Cells that express high levels of the C-terminal fragment display chromosome segregation defects, and ChIP analyses show changes in the status of pericentric beta-satellite and rDNA chromatin from heterochromatic to more euchromatic form. Finally, a cell line with inducible expression of the vigilin C-terminal fragment displays inducible alterations in beta-satellite chromatin. These and other results lead us to present a new model for vigilin-mediated, RNA-induced gene silencing.

How a small DNA virus uses dsRNA but not RNAi to regulate its life cycle

Mouse polyomavirus contains a circular DNA genome, with early and late genes transcribed from opposite strands. At early times after infection, genes encoded from the early transcription unit are predominantly expressed. After the onset of viral DNA replication, expression of genes encoded from the late transcription unit increases dramatically. At late times, late primary transcripts are inefficiently polyadenylated, leading to the generation of multigenomic RNAs that are precursors to mature mRNAs. These transcripts contain sequences complementary to the early RNAs and downregulate early-strand gene expression by inducing RNA editing. Our recent work leads to a model where the production of the multigenomic late RNAs is also controlled by the editing of poly(A) signals, directed by overlapping primary transcripts.

ADAR editing wobbles the microRNA world

Recent work reveals that adenosine-to-inosine editing occurs in a number of cellular microRNAs (miRNAs). Such editing is shown to diminish the expression of one miRNA and alter the target specificity of another. This changes our current views significantly by not only increasing the repertoire of miRNAs and their potential targets, but also providing mechanisms for how to regulate them and direct them to alternative targets.

Evidence for variation in abundance of antisense transcripts between multicellular animals but no relationship between antisense transcriptionand organismic complexity

Given that humans have about the same number of genes as mice and not so many more than worm, what makes us more complex? Antisense transcripts are implicated in many aspects of gene regulation. Is there a functional connection between antisense transcription and organismic complexity, that is, is antisense regulation especially prevalent in humans? We used the same robust protocol to identify antisense transcripts in humans and five other metazoan genomes (mouse, rat, chicken, fruit fly, and nematode), and found that the estimated proportions of genes involved in antisense transcription are highly sensitive to the number of transcripts included in the analysis. By controlling for transcript abundance, we find that the probability that any given transcript is putatively involved in sense-antisense regulation is no higher in humans than in other vertebrates but appears unusually high in flies and especially low in nematodes. Similarly, there is no evidence that the proportion of sense-antisense transcripts is especially higher in humans than other vertebrates in a given subset of transcript sequences such as mRNAs, coding sequences, conserved, or nonconserved transcripts. Although antisense transcription might be enriched in mammalian brains compared with nonbrain tissues, it is no more enriched in human brain than in mouse brain. Overall, therefore, while we see striking variation between multicellular animals in the abundance of antisense transcripts, there is no evidence for a link between antisense transcription and organismic complexity. More particularly, we see no evidence that humans are in any way unusual among the vertebrates in this regard. Instead, our results suggest that antisense transcription might be prevalent in almost all metazoan genomes, nematodes being an unexplained exception.

Evidence for a preferential targeting of 3'-UTRs by cis-encoded natural antisense transcripts

Although both the 5'- and 3'-untranslated regions (5'- and 3'-UTRs) of eukaryotic mRNAs may play a crucial role in posttranscriptional gene regulation, we observe that cis-encoded natural antisense RNAs have a striking preferential complementarity to the 3'-UTRs of their target genes in mammalian (human and mouse) genomes. A null neutral model, evoking differences in the rate of 3'-UTR and 5'-UTR extension, could potentially explain high rates of 3'-to-3' overlap compared with 5'-to-5' overlap. However, employing a simulation model we show that this null model probably cannot explain the finding that 3'-to-3' overlapping pairs have a much higher probability (>5 times) of conservation in both mouse and human genomes with the same overlapping pattern than do 5'-to-5' overlaps. Furthermore, it certainly cannot explain the finding that overlapping pairs seen in both genomes have a significantly higher probability of having co-expression and inverse expression (i.e. characteristic of sense-antisense regulation) than do overlapping pairs seen in only one of the two species. We infer that the function of many 3'-to-3' overlaps is indeed antisense regulation. These findings underscore the preference for, and conservation of, 3'-UTR-targeted antisense regulation, and the importance of 3'-UTRs in gene regulation.

Retention and repression: fates of hyperedited RNAs in the nucleus

Double-stranded RNA (dsRNA) is often formed in the nuclei of mammalian cells, but in this compartment it does not induce the effects characteristic of cytoplasmic dsRNA. Rather, recent work has suggested that nuclear dsRNA is a target for the ADAR class of enzymes, which deaminate adenosines to inosines. Further, there are a number of distinct fates of such edited RNA, including nuclear retention and perhaps also gene silencing.

Genome-wide analysis of coordinate expression and evolution of human cis-encoded sense-antisense transcripts

Is sense-antisense (SA) pairing of transcripts a common mode of gene regulation in the human genome? Although >20% of human genes might form SA pairs, the extent to which they are involved in antisense regulation is unknown. Simultaneous expression of paired sense and antisense genes is an essential step and an important indicator of antisense regulation. In this article, we demonstrate that human SA pairs tend to be co-expressed and/or inversely expressed more frequently than expected by chance. Moreover, co-expressed and inversely expressed SA pairs exhibit a striking pattern of evolutionary conservation. These findings suggest that antisense regulation is a common and important mechanism of gene regulation in the human genome.

Vigilins bind to promiscuously A-to-I-edited RNAs and are involved in the formation of heterochromatin

The fate of double-stranded RNA (dsRNA) in the cell depends on both its length and location . The expression of dsRNA in the nucleus leads to several distinct consequences. First, the promiscuous deamination of adenosines to inosines by dsRNA-specific adenosine deaminase (ADAR) can lead to the nuclear retention of edited transcripts . Second, dsRNAs might induce heterochromatic gene silencing through an RNAi-related mechanism . Is RNA editing also connected to heterochromatin? We report that members of the conserved Vigilin class of proteins have a high affinity for inosine-containing RNAs. In agreement with other work , we find that these proteins localize to heterochromatin and that mutation or depletion of the Drosophila Vigilin, DDP1, leads to altered nuclear morphology and defects in heterochromatin and chromosome segregation. Furthermore, nuclear Vigilin is found in complexes containing not only the editing enzyme ADAR1 but also RNA helicase A and Ku86/70. In the presence of RNA, the Vigilin complex recruits the DNA-PKcs enzyme, which appears to phosphorylate a discrete set of targets, some or all of which are known to participate in chromatin silencing. These results are consistent with a mechanistic link between components of the DNA-repair machinery and RNA-mediated gene silencing.

Human antisense genes have unusually short introns: evidence for selection for rapid transcription

We suggest that small introns are favored in genes requiring a minimal response time ('nimble' genes), and that antisense genes might be prime candidates for such nimble genes. In this article, we show that antisense genes have significantly shorter introns than all other gene categories, which supports our proposal that the short introns of antisense genes might be functionally important and integrally related to their potential role in efficient gene regulation.

SINEs point to abundant editing in the human genome

Recent bioinformatic analyses suggest that almost all human transcripts are edited by adenosine deaminases; most of this editing is in Alu elements.

Recent bioinformatic analyses suggest that almost all human transcripts are edited by adenosine deaminases (ADARs), converting adenosines to inosines. Most of this editing is in Alu element transcripts, which are unique to primates. This editing might have no function or might be involved in functions such as the regulation of splicing, chromatin or nuclear localization of transcripts.

Effects of length and location on the cellular response to double-stranded RNA

Since double-stranded RNA (dsRNA) has not until recently generally been thought to be deliberately expressed in cells, it has commonly been assumed that the major source of cellular dsRNA is viral infections. In this view, the cellular responses to dsRNA would be natural and perhaps ancient antiviral responses. While the cell may certainly react to some dsRNAs as an antiviral response, this does not represent the only response or even, perhaps, the major one. A number of recent observations have pointed to the possibility that dsRNA molecules are not seen only as evidence of viral infection or recognized for degradation because they cannot be translated. In some instances they may also play important roles in normal cell growth and function. The purpose of this review is to outline our current understanding of the fate of dsRNA in cells, with a focus on the apparent fact that their fates and functions appear to depend critically not only on where in the cell dsRNA molecules are found, but also on how long they are and perhaps on how abundant they are.

Methods for the analysis of adenosine-to-inosine editing in RNA

In this work we describe methods for the analysis of RNAs that have been edited by the double-strand RNA-specific adenosine deaminase, ADAR. These RNAs contain inosine residues that can be detected and quantified by a variety of approaches, including base hydrolysis and thin-layer chromatography, reverse transcription polymerase chain reaction, primer extension, and inosine-specific base cleavage. The most common method for the analysis of editing will be described here. This method involves complete hydrolysis of edited RNAs to nucleoside monophosphates, followed by separation of the products using thin-layer chromatography.

The fate of dsRNA in the nucleus: a p54(nrb)-containing complex mediates the nuclear retention of promiscuously A-to-I edited RNAs

How do cells discriminate between selectively edited mRNAs that encode new protein isoforms, and dsRNA-induced, promiscuously edited RNAs that encode nonfunctional, mutant proteins? We have developed a Xenopus oocyte model system which shows that a variety of hyperedited, inosine-containing RNAs are specifically retained in the nucleus. To uncover the mechanism of inosine-induced retention, HeLa cell nuclear extracts were used to isolate a multiprotein complex that binds specifically and cooperatively to inosine-containing RNAs. This complex contains the inosine-specific RNA binding protein p54(nrb), the splicing factor PSF, and the inner nuclear matrix structural protein matrin 3. We provide evidence that one function of the complex identified here is to anchor hyperedited RNAs to the nuclear matrix, while allowing selectively edited mRNAs to be exported.

Intronless mRNA transport elements may affect multiple steps of pre-mRNA processing

We have reported recently that a small element within the mouse histone H2a-coding region permits efficient cytoplasmic accumulation of intronless beta-globin cDNA transcripts. This sequence lowers the levels of spliced products from intron-containing constructs and can functionally replace Rev and the Rev-responsive element (RRE) in the nuclear export of unspliced HIV-1-related mRNAs. In work reported here, we further investigate the molecular mechanisms by which this element might work. We demonstrate here through both in vivo and in vitro assays that, in addition to promoting mRNA nuclear export, this element acts as a polyadenylation enhancer and as a potent inhibitor of splicing. Surprisingly, two other described intronless mRNA transport elements (from the herpes simplex virus thymidine kinase gene and hepatitis B virus) appear to function in a similar manner. These findings prompt us to suggest that a general feature of intronless mRNA transport elements might be a collection of phenotypes, including the inhibition of splicing and the enhancement of both polyadenylation and mRNA export.