Alternative processing as evolutionary mechanism for the origin of novel nonprotein coding RNAs

Reference Genes across Nine Brain Areas of Wild Type and Prader-Willi Syndrome Mice: Assessing Differences in Igfbp7, Pcsk1, Nhlh2 and Nlgn3 Expression

Prader−Willi syndrome (PWS) is a complex neurodevelopmental disorder caused by the deletion or inactivation of paternally expressed imprinted genes at the chromosomal region 15q11−q13. The PWS-critical region (PWScr) harbors tandemly repeated non-protein coding IPW-A exons hosting the intronic SNORD116 snoRNA gene array that is predominantly expressed in brain. Paternal deletion of PWScr is associated with key PWS symptoms in humans and growth retardation in mice (PWScr model). Dysregulation of the hypothalamic−pituitary axis (HPA) is thought to be causally involved in the PWS phenotype. Here we performed a comprehensive reverse transcription quantitative PCR (RT-qPCR) analysis across nine different brain regions of wild-type (WT) and PWScr mice to identify stably expressed reference genes. Four methods (Delta Ct, BestKeeper, Normfinder and Genorm) were applied to rank 11 selected reference gene candidates according to their expression stability. The resulting panel consists of the top three most stably expressed genes suitable for gene-expression profiling and comparative transcriptome analysis of WT and/or PWScr mouse brain regions. Using these reference genes, we revealed significant differences in the expression patterns of Igfbp7, Nlgn3 and three HPA associated genes: Pcsk1, Pcsk2 and Nhlh2 across investigated brain regions of wild-type and PWScr mice. Our results raise a reasonable doubt on the involvement of the Snord116 in posttranscriptional regulation of Nlgn3 and Nhlh2 genes. We provide a valuable tool for expression analysis of specific genes across different areas of the mouse brain and for comparative investigation of PWScr mouse models to discover and verify different regulatory pathways affecting this complex disorder.

A Comprehensive Review of Genetically Engineered Mouse Models for Prader-Willi Syndrome Research

Prader-Willi syndrome (PWS) is a neurogenetic multifactorial disorder caused by the deletion or inactivation of paternally imprinted genes on human chromosome 15q11-q13. The affected homologous locus is on mouse chromosome 7C. The positional conservation and organization of genes including the imprinting pattern between mice and men implies similar physiological functions of this locus. Therefore, considerable efforts to recreate the pathogenesis of PWS have been accomplished in mouse models. We provide a summary of different mouse models that were generated for the analysis of PWS and discuss their impact on our current understanding of corresponding genes, their putative functions and the pathogenesis of PWS. Murine models of PWS unveiled the contribution of each affected gene to this multi-facetted disease, and also enabled the establishment of the minimal critical genomic region (PWScr) responsible for core symptoms, highlighting the importance of non-protein coding genes in the PWS locus. Although the underlying disease-causing mechanisms of PWS remain widely unresolved and existing mouse models do not fully capture the entire spectrum of the human PWS disorder, continuous improvements of genetically engineered mouse models have proven to be very powerful and valuable tools in PWS research.

Circular RNA Encoded Amyloid Beta peptides-A Novel Putative Player in Alzheimer's Disease

Alzheimer's disease (AD) is an age-related detrimental dementia. Amyloid beta peptides (Aβ) play a crucial role in the pathology of AD. In familial AD, Aβ are generated from the full-length amyloid beta precursor protein (APP) via dysregulated proteolytic processing; however, in the case of sporadic AD, the mechanism of Aβ biogenesis remains elusive. circRNAs are a class of transcripts preferentially expressed in brain. We identified a circRNA harboring the Aβ-coding region of the APP gene termed circAβ-a. This circular RNA was detected in the brains of AD patients and non-dementia controls. With the aid of our recently established approach for analysis of circRNA functions, we demonstrated that circAβ-a is efficiently translated into a novel Aβ-containing Aβ175 polypeptide (19.2 KDa) in both cultured cells and human brain. Furthermore, Aβ175 was shown to be processed into Aβ peptides-a hallmark of AD. In summary, our analysis revealed an alternative pathway of Aβ biogenesis. Consequently, circAβ-a and its corresponding translation product could potentially represent novel therapeutic targets for AD treatment. Importantly, our data point to yet another evolutionary route for potentially increasing proteome complexity by generating additional polypeptide variants using back-splicing of primary transcripts that yield circular RNA templates.

Ectopic expression of Snord115 in choroid plexus interferes with editing but not splicing of 5-Ht2c receptor pre-mRNA in mice

Serotonin 5-HT2C receptor is a G-protein coupled excitatory receptor that regulates several biochemical pathways and has been implicated in obesity, mental state, sleep cycles, autism, neuropsychiatric disorders and neurodegenerative diseases. The activity of 5-HT2CR is regulated via alternative splicing and A to I editing of exon Vb of its pre-mRNA. Snord115 is a small nucleolar RNA that is expressed in mouse neurons and displays an 18-nucleotide base complementary to exon Vb of 5-HT2CR pre-mRNA. For almost two decades this putative guide element of Snord115 has wandered like a ghost through the literature in attempts to elucidate the biological significance of this complementarity. In mice, Snord115 is expressed in neurons and absent in the choroid plexus where, in contrast, 5-Ht2cr mRNA is highly abundant. Here we report the analysis of 5-Ht2cr pre-mRNA posttranscriptional processing via RNA deep sequencing in a mouse model that ectopically expresses Snord115 in the choroid plexus. In contrast to previous reports, our analysis demonstrated that Snord115 does not control alternative splicing of 5-Ht2cr pre-mRNA in vivo. We identified a modest, yet statistically significant reduction of 5-Ht2cr pre-mRNA A to I editing at the major A, B, C and D sites. We suggest that Snord115 and exon Vb of 5Ht2cr pre-mRNA form a double-stranded structure that is subject to ADAR-mediated A to I editing. To the best of our knowledge, this is the first comprehensive Snord115 gain-of-function analysis based on in vivo mouse models.

Exaptation at the molecular genetic level

The realization that body parts of animals and plants can be recruited or coopted for novel functions dates back to, or even predates the observations of Darwin. S.J. Gould and E.S. Vrba recognized a mode of evolution of characters that differs from adaptation. The umbrella term aptation was supplemented with the concept of exaptation. Unlike adaptations, which are restricted to features built by selection for their current role, exaptations are features that currently enhance fitness, even though their present role was not a result of natural selection. Exaptations can also arise from nonaptations; these are characters which had previously been evolving neutrally. All nonaptations are potential exaptations. The concept of exaptation was expanded to the molecular genetic level which aided greatly in understanding the enormous potential of neutrally evolving repetitive DNA-including transposed elements, formerly considered junk DNA-for the evolution of genes and genomes. The distinction between adaptations and exaptations is outlined in this review and examples are given. Also elaborated on is the fact that such distinctions are sometimes more difficult to determine; this is a widespread phenomenon in biology, where continua abound and clear borders between states and definitions are rare.

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.

BC1 RNA motifs required for dendritic transport in vivo

BC1 RNA is a small brain specific non-protein coding RNA. It is transported from the cell body into dendrites where it is involved in the fine-tuning translational control. Due to its compactness and established secondary structure, BC1 RNA is an ideal model for investigating the motifs necessary for dendritic localization. Previously, microinjection of in vitro transcribed BC1 RNA mutants into the soma of cultured primary neurons suggested the importance of RNA motifs for dendritic targeting. These ex vivo experiments identified a single bulged nucleotide (U22) and a putative K-turn (GA motif) structure required for dendritic localization or distal transport, respectively. We generated six transgenic mouse lines (three founders each) containing neuronally expressing BC1 RNA variants on a BC1 RNA knockout mouse background. In contrast to ex vivo data, we did not find indications of reduction or abolition of dendritic BC1 RNA localization in the mutants devoid of the GA motif or the bulged nucleotide. We confirmed the ex vivo data, which showed that the triloop terminal sequence had no consequence on dendritic transport. Interestingly, changing the triloop supporting structure completely abolished dendritic localization of BC1 RNA. We propose a novel RNA motif important for dendritic transport in vivo.

Unusual Novel SnoRNA-Like RNAs in Drosophila melanogaster

A computational screen for novel small nucleolar RNAs in Drosophila melanogaster uncovered 15 novel snoRNAs and snoRNA-like long non-coding RNAs. In contrast to earlier surverys, the novel sequences are mostly poorly conserved and originate from unusual genomic locations. The majority derive from precurors antisense to well-known protein-coding genes, and four of the candidates are produced from exon-coding regions. Only a minority of the new sequences appears to have canonical target sites in ribosomal or small nuclear RNAs. Taken together, these evolutionary young, poorly conserved, and genomically atypical sequences point at a class of snoRNA-like transcripts with predominantly regulatory functions in the fruit fly genome.

Evidence for a Novel Mechanism of Influenza Virus-Induced Type I Interferon Expression by a Defective RNA-Encoded Protein

Influenza A virus (IAV) defective RNAs are generated as byproducts of error-prone viral RNA replication. They are commonly derived from the larger segments of the viral genome and harbor deletions of various sizes resulting in the generation of replication incompatible viral particles. Furthermore, small subgenomic RNAs are known to be strong inducers of pattern recognition receptor RIG-I-dependent type I interferon (IFN) responses. The present study identifies a novel IAV-induced defective RNA derived from the PB2 segment of A/Thailand/1(KAN-1)/2004 (H5N1). It encodes a 10 kDa protein (PB2_∆) sharing the N-terminal amino acid sequence of the parental PB2 protein followed by frame shift after internal deletion. PB2_∆ induces the expression of IFNβ and IFN-stimulated genes by direct interaction with the cellular adapter protein MAVS, thereby reducing viral replication of IFN-sensitive viruses such as IAV or vesicular stomatitis virus. This induction of IFN is completely independent of the defective RNA itself that usually serves as pathogen-associated pattern and thus does not require the cytoplasmic sensor RIG-I. These data suggest that not only defective RNAs, but also some defective RNA-encoded proteins can act immunostimulatory. In this particular case, the KAN-1-induced defective RNA-encoded protein PB2_∆ enhances the overwhelming immune response characteristic for highly pathogenic H5N1 viruses, leading to a more severe phenotype in vivo.

Error-prone polymerase function of RNA viruses can result in expression of defective RNAs harboring internal deletions of various sizes. Small subgenomic RNAs are strong inducers of the antiviral response by serving as pathogen-associated patterns that are predominantly detected by cellular sensors. Recently, it has been shown that influenza A virus defective RNAs are not only generated upon passages in cell culture, but also in infected humans, indicating that these subgenomic RNAs may also be relevant in infections in vivo. Here, we characterize a novel defective RNA derived from the PB2 segment of a highly pathogenic H5N1 influenza A virus. This RNA encodes a 10 kDa peptide (PB2_Δ) which activates type I interferon (IFN) responses through direct interaction with the adapter protein MAVS, a key component of the RIG-I-dependent IFN induction. This is the first time that such a function was described for a defective RNA-encoded protein, a finding that has several important implications with regard to deciphering viral protein functions and options for immunostimulatory approaches. Furthermore, this is an example of how influenza viruses may acquire novel polypeptides with altered functions from its limited genome.

Differential regulation of non-protein coding RNAs from Prader-Willi Syndrome locus

Prader-Willi Syndrome (PWS) is a neurogenetic disorder caused by the deletion of imprinted genes on the paternally inherited human chromosome 15q11-q13. This locus harbours a long non-protein-coding RNA (U-UBE3A-ATS) that contains six intron-encoded snoRNAs, including the SNORD116 and SNORD115 repetitive clusters. The 3′-region of U-UBE3A-ATS is transcribed in the cis-antisense direction to the ubiquitin-protein ligase E3A (UBE3A) gene. Deletion of the SNORD116 region causes key characteristics of PWS. There are few indications that SNORD115 might regulate serotonin receptor (5HT2C) pre-mRNA processing. Here we performed quantitative real-time expression analyses of RNAs from the PWS locus across 20 human tissues and combined it with deep-sequencing data derived from Cap Analysis of Gene Expression (CAGE-seq) libraries. We found that the expression profiles of SNORD64, SNORD107, SNORD108 and SNORD116 are similar across analyzed tissues and correlate well with SNORD116 embedded U-UBE3A-ATS exons (IPW116). Notable differences in expressions between the aforementioned RNAs and SNORD115 together with the host IPW115 and UBE3A cis-antisense exons were observed. CAGE-seq analysis revealed the presence of potential transcriptional start sites originated from the U-UBE3A-ATS spanning region. Our findings indicate novel aspects for the expression regulation in the PWS locus.

The persistent contributions of RNA to eukaryotic gen(om)e architecture and cellular function

Currently, the best scenario for earliest forms of life is based on RNA molecules as they have the proven ability to catalyze enzymatic reactions and harbor genetic information. Evolutionary principles valid today become apparent in such models already. Furthermore, many features of eukaryotic genome architecture might have their origins in an RNA or RNA/protein (RNP) world, including the onset of a further transition, when DNA replaced RNA as the genetic bookkeeper of the cell. Chromosome maintenance, splicing, and regulatory function via RNA may be deeply rooted in the RNA/RNP worlds. Mostly in eukaryotes, conversion from RNA to DNA is still ongoing, which greatly impacts the plasticity of extant genomes. Raw material for novel genes encoding protein or RNA, or parts of genes including regulatory elements that selection can act on, continues to enter the evolutionary lottery.