Arabidopsis mRNA secondary structure correlates with protein function and domains

The p53 mRNA: an integral part of the cellular stress response

A large number of signalling pathways converge on p53 to induce different cellular stress responses that aim to promote cell cycle arrest and repair or, if the damage is too severe, to induce irreversible senescence or apoptosis. The differentiation of p53 activity towards specific cellular outcomes is tightly regulated via a hierarchical order of post-translational modifications and regulated protein-protein interactions. The mechanisms governing these processes provide a model for how cells optimize the genetic information for maximal diversity. The p53 mRNA also plays a role in this process and this review aims to illustrate how protein and RNA interactions throughout the p53 mRNA in response to different signalling pathways control RNA stability, translation efficiency or alternative initiation of translation. We also describe how a p53 mRNA platform shows riboswitch-like features and controls the rate of p53 synthesis, protein stability and modifications of the nascent p53 protein. A single cancer-derived synonymous mutation disrupts the folding of this platform and prevents p53 activation following DNA damage. The role of the p53 mRNA as a target for signalling pathways illustrates how mRNA sequences have co-evolved with the function of the encoded protein and sheds new light on the information hidden within mRNAs.

Shaping the regulation of the p53 mRNA tumour suppressor: the co-evolution of genetic signatures

Structured RNA regulatory motifs exist from the prebiotic stages of the RNA world to the more complex eukaryotic systems. In cases where a functional RNA structure is within the coding sequence a selective pressure drives a parallel co-evolution of the RNA structure and the encoded peptide domain. The p53-MDM2 axis, describing the interactions between the p53 tumor suppressor and the MDM2 E3 ubiquitin ligase, serves as particularly useful model revealing how secondary RNA structures have co-evolved along with corresponding interacting protein motifs, thus having an impact on protein - RNA and protein - protein interactions; and how such structures developed signal-dependent regulation in mammalian systems. The p53(BOX-I) RNA sequence binds the C-terminus of MDM2 and controls p53 synthesis while the encoded peptide domain binds MDM2 and controls p53 degradation. The BOX-I peptide domain is also located within p53 transcription activation domain. The folding of the p53 mRNA structure has evolved from temperature-regulated in pre-vertebrates to an ATM kinase signal-dependent pathway in mammalian cells. The protein - protein interaction evolved in vertebrates and became regulated by the same signaling pathway. At the same time the protein - RNA and protein - protein interactions evolved, the p53 trans-activation domain progressed to become integrated into a range of cellular pathways. We discuss how a single synonymous mutation in the BOX-1, the p53(L22 L), observed in a chronic lymphocyte leukaemia patient, prevents the activation of p53 following DNA damage. The concepts analysed and discussed in this review may serve as a conceptual mechanistic paradigm of the co-evolution and function of molecules having roles in cellular regulation, or the aetiology of genetic diseases and how synonymous mutations can affect the encoded protein.

Conserved Secondary Structures in Viral mRNAs

RNA secondary structure in untranslated and protein coding regions has been shown to play an important role in regulatory processes and the viral replication cycle. While structures in non-coding regions have been investigated extensively, a thorough overview of the structural repertoire of protein coding mRNAs, especially for viruses, is lacking. Secondary structure prediction of large molecules, such as long mRNAs remains a challenging task, as the contingent of structures a sequence can theoretically fold into grows exponentially with sequence length. We applied a structure prediction pipeline to Viral Orthologous Groups that first identifies the local boundaries of potentially structured regions and subsequently predicts their functional importance. Using this procedure, the orthologous groups were split into structurally homogenous subgroups, which we call subVOGs. This is the first compilation of potentially functional conserved RNA structures in viral coding regions, covering the complete RefSeq viral database. We were able to recover structural elements from previous studies and discovered a variety of novel structured regions. The subVOGs are available through our web resource RNASIV (RNA structure in viruses; http://rnasiv.bio.wzw.tum.de).

Adenine Enrichment at the Fourth CDS Residue in Bacterial Genes Is Consistent with Error Proofing for +1 Frameshifts

Beyond selection for optimal protein functioning, coding sequences (CDSs) are under selection at the RNA and DNA levels. Here, we identify a possible signature of “dual-coding,” namely extensive adenine (A) enrichment at bacterial CDS fourth sites. In 99.07% of studied bacterial genomes, fourth site A use is greater than expected given genomic A-starting codon use. Arguing for nucleotide level selection, A-starting serine and arginine second codons are heavily utilized when compared with their non-A starting synonyms. Several models have the ability to explain some of this trend. In part, A-enrichment likely reduces 5′ mRNA stability, promoting translation initiation. However T/U, which may also reduce stability, is avoided. Further, +1 frameshifts on the initiating ATG encode a stop codon (TGA) provided A is the fourth residue, acting either as a frameshift “catch and destroy” or a frameshift stop and adjust mechanism and hence implicated in translation initiation. Consistent with both, genomes lacking TGA stop codons exhibit weaker fourth site A-enrichment. Sequences lacking a Shine–Dalgarno sequence and those without upstream leader genes, that may be more error prone during initiation, have greater utilization of A, again suggesting a role in initiation. The frameshift correction model is consistent with the notion that many genomic features are error-mitigation factors and provides the first evidence for site-specific out of frame stop codon selection. We conjecture that the NTG universal start codon may have evolved as a consequence of TGA being a stop codon and the ability of NTGA to rapidly terminate or adjust a ribosome.

RNA thermosensors: how might animals exploit their regulatory potential?

The secondary and tertiary orders of RNA structure are crucial for a suite of RNA-related functions, including regulation of translation, gene expression and RNA turnover. The temperature sensitivity of RNA secondary and tertiary structures is exploited by bacteria to fabricate RNA thermosensing systems that allow a rapid adaptive response to temperature change. RNA thermometers (RNATs) present in non-coding regions of certain mRNAs of pathogenic bacteria enable rapid upregulation of translation of virulence proteins when the temperature of the bacterium rises after entering a mammalian host. Rapid upregulation of translation of bacterial heat-shock proteins likewise is governed in part by RNATs. Turnover of mRNA may be regulated by temperature-sensitive RNA structures. Whereas the roles of temperature-sensitive RNA structures similar to RNATs in Eukarya and Archaea are largely unknown, there would appear to be a potential for all taxa to adaptively regulate their thermal physiology through exploitation of RNA-based thermosensory responses akin to those of bacteria. In animals, these responses might include regulation of translation of stress-induced proteins, alternative splicing of messenger RNA precursors, differential expression of allelic proteins, modulation of activities of small non-coding RNAs, regulation of mRNA turnover and control of RNA editing. New methods for predicting, detecting and experimentally modifying RNA secondary structure offer promising windows into these fascinating aspects of RNA biochemistry. Elucidating whether animals too have exploited the types of RNA thermosensing tools that are used so effectively by bacteria seems likely to provide exciting new insights into the mechanisms of evolutionary adaptation and acclimatization to temperature.

The Evolutionary Basis of Translational Accuracy in Plants

Mistranslation errors compromise fitness by wasting resources on nonfunctional proteins. In order to reduce the cost of mistranslations, natural selection chooses the most accurately translated codons at sites that are particularly important for protein structure and function. We investigated the determinants underlying selection for translational accuracy in several species of plants belonging to three clades: Brassicaceae, Fabidae, and Poaceae. Although signatures of translational selection were found in genes from a wide range of species, the underlying factors varied in nature and intensity. Indeed, the degree of synonymous codon bias at evolutionarily conserved sites varied among plant clades while remaining uniform within each clade. This is unlikely to solely reflect the diversity of tRNA pools because there is little correlation between synonymous codon bias and tRNA abundance, so other factors must affect codon choice and translational accuracy in plant genes. Accordingly, synonymous codon choice at a given site was affected not only by the selection pressure at that site, but also its participation in protein domains or mRNA secondary structures. Although these effects were detected in all the species we analyzed, their impact on translation accuracy was distinct in evolutionarily distant plant clades. The domain effect was found to enhance translational accuracy in dicot and monocot genes with a high GC content, but to oppose the selection of more accurate codons in monocot genes with a low GC content.

JNSViewer-A JavaScript-based Nucleotide Sequence Viewer for DNA/RNA secondary structures

Many tools are available for visualizing RNA or DNA secondary structures, but there is scarce implementation in JavaScript that provides seamless integration with the increasingly popular web computational platforms. We have developed JNSViewer, a highly interactive web service, which is bundled with several popular tools for DNA/RNA secondary structure prediction and can provide precise and interactive correspondence among nucleotides, dot-bracket data, secondary structure graphs, and genic annotations. In JNSViewer, users can perform RNA secondary structure predictions with different programs and settings, add customized genic annotations in GFF format to structure graphs, search for specific linear motifs, and extract relevant structure graphs of sub-sequences. JNSViewer also allows users to choose a transcript or specific segment of Arabidopsis thaliana genome sequences and predict the corresponding secondary structure. Popular genome browsers (i.e., JBrowse and BrowserGenome) were integrated into JNSViewer to provide powerful visualizations of chromosomal locations, genic annotations, and secondary structures. In addition, we used StructureFold with default settings to predict some RNA structures for Arabidopsis by incorporating in vivo high-throughput RNA structure profiling data and stored the results in our web server, which might be a useful resource for RNA secondary structure studies in plants. JNSViewer is available at http://bioinfolab.miamioh.edu/jnsviewer/index.html.

Insights into the adaptive response of Arabidopsis thaliana to prolonged thermal stress by ribosomal profiling and RNA-Seq

BACKGROUND

Environmental stress puts organisms at risk and requires specific stress-tailored responses to maximize survival. Long-term exposure to stress necessitates a global reprogramming of the cellular activities at different levels of gene expression.

RESULTS

Here, we use ribosome profiling and RNA sequencing to globally profile the adaptive response of Arabidopsis thaliana to prolonged heat stress. To adapt to long heat exposure, the expression of many genes is modulated in a coordinated manner at a transcriptional and translational level. However, a significant group of genes opposes this trend and shows mainly translational regulation. Different secondary structure elements are likely candidates to play a role in regulating translation of those genes.

CONCLUSIONS

Our data also uncover on how the subunit stoichiometry of multimeric protein complexes in plastids is maintained upon heat exposure.

Genome-Wide Analysis of RNA Secondary Structure

Single-stranded RNA molecules fold into extraordinarily complicated secondary and tertiary structures as a result of intramolecular base pairing. In vivo, these RNA structures are not static. Instead, they are remodeled in response to changes in the prevailing physicochemical environment of the cell and as a result of intermolecular base pairing and interactions with RNA-binding proteins. Remarkable technical advances now allow us to probe RNA secondary structure at single-nucleotide resolution and genome-wide, both in vitro and in vivo. These data sets provide new glimpses into the RNA universe. Analyses of RNA structuromes in HIV, yeast, Arabidopsis, and mammalian cells and tissues have revealed regulatory effects of RNA structure on messenger RNA (mRNA) polyadenylation, splicing, translation, and turnover. Application of new methods for genome-wide identification of mRNA modifications, particularly methylation and pseudouridylation, has shown that the RNA "epitranscriptome" both influences and is influenced by RNA structure. In this review, we describe newly developed genome-wide RNA structure-probing methods and synthesize the information emerging from their application.

The Conservation and Function of RNA Secondary Structure in Plants

RNA transcripts fold into secondary structures via intricate patterns of base pairing. These secondary structures impart catalytic, ligand binding, and scaffolding functions to a wide array of RNAs, forming a critical node of biological regulation. Among their many functions, RNA structural elements modulate epigenetic marks, alter mRNA stability and translation, regulate alternative splicing, transduce signals, and scaffold large macromolecular complexes. Thus, the study of RNA secondary structure is critical to understanding the function and regulation of RNA transcripts. Here, we review the origins, form, and function of RNA secondary structure, focusing on plants. We then provide an overview of methods for probing secondary structure, from physical methods such as X-ray crystallography and nuclear magnetic resonance (NMR) imaging to chemical and nuclease probing methods. Combining these latter methods with high-throughput sequencing has enabled them to scale across whole transcriptomes, yielding tremendous new insights into the form and function of RNA secondary structure.

p53 mRNA and p53 Protein Structures Have Evolved Independently to Interact with MDM2

The p53 tumor suppressor and its key regulator MDM2 play essential roles in development, ageing, cancer, and cellular stress responses in mammals. Following DNA damage, MDM2 interacts with p53 mRNA in an ATM kinase-dependent fashion and stimulates p53 synthesis, whereas under normal conditions, MDM2 targets the p53 protein for degradation. The peptide- and RNA motifs that interact with MDM2 are encoded by the same conserved BOX-I sequence, but how these interactions have evolved is unknown. Here, we show that a temperature-sensitive structure in the invertebrate Ciona intestinalis (Ci) p53 mRNA controls its interaction with MDM2. We also show that a nonconserved flanking region of Ci-BOX-I domain prevents the p53-MDM2 protein-protein interaction. These results indicate that the temperature-regulated p53 mRNA-MDM2 interaction evolved to become kinase regulated in the mammalian DNA damage response. The data also suggest that the negative regulation of p53 by MDM2 via protein-protein interaction evolved in vertebrates following changes in the BOX-I flanking sequence.

Chemical Modifications Mark Alternatively Spliced and Uncapped Messenger RNAs in Arabidopsis

Posttranscriptional chemical modification of RNA bases is a widespread and physiologically relevant regulator of RNA maturation, stability, and function. While modifications are best characterized in short, noncoding RNAs such as tRNAs, growing evidence indicates that mRNAs and long noncoding RNAs (lncRNAs) are likewise modified. Here, we apply our high-throughput annotation of modified ribonucleotides (HAMR) pipeline to identify and classify modifications that affect Watson-Crick base pairing at three different levels of the Arabidopsis thaliana transcriptome (polyadenylated, small, and degrading RNAs). We find this type of modifications primarily within uncapped, degrading mRNAs and lncRNAs, suggesting they are the cause or consequence of RNA turnover. Additionally, modifications within stable mRNAs tend to occur in alternatively spliced introns, suggesting they regulate splicing. Furthermore, these modifications target mRNAs with coherent functions, including stress responses. Thus, our comprehensive analysis across multiple RNA classes yields insights into the functions of covalent RNA modifications in plant transcriptomes.

Transcriptome-wide measurement of plant RNA secondary structure

RNAs fold into intricate and precise secondary structures. These structural patterns regulate multiple steps of the RNA lifecycle, while also conferring catalytic and scaffolding functions to certain transcripts. Therefore, a full understanding of RNA posttranscriptional regulation requires a comprehensive picture of secondary structure. Here, we review several high throughput sequencing-based methods to globally survey plant RNA secondary structure. These methods are more accurate than computational prediction, and more scalable than physical techniques such as crystallography. We note hurdles to reliably measuring secondary structure, including RNA-binding proteins, RNA base modifications, and intramolecular duplexes. Finally, we survey the functional knowledge that has been gleaned from each of these methods, and identify some unanswered questions that remain.

Toward a next-generation atlas of RNA secondary structure

RNA structure plays a crucial role in gene maturation, regulation and function. Determining the form and frequency of RNA folds is essential for a better understanding of how RNA exerts its functions. Low-throughput studies have focused on RNA primary sequences and expression levels, but with an emphasis on relatively small numbers of transcripts. However, with the recent advent of high-throughput technologies, it is realistic to begin analyzing RNA secondary structures on a genome-wide scale. Here, we review genome-wide RNA secondary structure profiles as well as advances in computational structure predictions. We further discuss the novel characteristics of RNA secondary structure across messenger RNAs. Probing RNA secondary structure by high-throughput sequencing will enable us to build atlases of RNA secondary structures, an important step in helping us to understand the versatility of RNA functions in diverse cellular processes.

Structural imprints in vivo decode RNA regulatory mechanisms

Visualizing the physical basis for molecular behavior inside living cells is a grand challenge in biology. RNAs are central to biological regulation, and RNA’s ability to adopt specific structures intimately controls every step of the gene expression program¹. However, our understanding of physiological RNA structures is limited; current in vivo RNA structure profiles view only two of four nucleotides that make up RNA^2,3. Here we present a novel biochemical approach, In Vivo Click SHAPE (icSHAPE), that enables the first global view of RNA secondary structures of all four bases in living cells. icSHAPE of mouse embryonic stem cell transcriptome versus purified RNA folded in vitro shows that the structural dynamics of RNA in the cellular environment distinguishes different classes of RNAs and regulatory elements. Structural signatures at translational start sites and ribosome pause sites are conserved from in vitro, suggesting that these RNA elements are programmed by sequence. In contrast, focal structural rearrangements in vivo reveal precise interfaces of RNA with RNA binding proteins or RNA modification sites that are consistent with atomic-resolution structural data. Such dynamic structural footprints enable accurate prediction of RNA-protein interactions and N6-methyladenosine (m⁶A) modification genome-wide. These results open the door for structural genomics of RNA in living cells and reveal key physiological structures controlling gene expression.

Insights into RNA structure and function from genome-wide studies

A comprehensive understanding of RNA structure will provide fundamental insights into the cellular function of both coding and non-coding RNAs. Although many RNA structures have been analysed by traditional biophysical and biochemical methods, the low-throughput nature of these approaches has prevented investigation of the vast majority of cellular transcripts. Triggered by advances in sequencing technology, genome-wide approaches for probing the transcriptome are beginning to reveal how RNA structure affects each step of protein expression and RNA stability. In this Review, we discuss the emerging relationships between RNA structure and the regulation of gene expression.

Landes Highlights

CRM1 and NMD3 localize to the nucleolus and affect rRNA synthesis

mRNA secondary structure correlates with protein function and domains

Transcriptional changes induced by tumor dormancy-associated microRNA-190

Ribonuclease binase: A potential anticancer therapeutic