Reproducible Analysis of Sequencing-Based RNA Structure Probing Data with User-Friendly Tools

Light-Dependent Translation Change of Arabidopsis psbA Correlates with RNA Structure Alterations at the Translation Initiation Region

mRNA secondary structure influences translation. Proteins that modulate the mRNA secondary structure around the translation initiation region may regulate translation in plastids. To test this hypothesis, we exposed Arabidopsis thaliana to high light, which induces translation of psbA mRNA encoding the D1 subunit of photosystem II. We assayed translation by ribosome profiling and applied two complementary methods to analyze in vivo RNA secondary structure: DMS-MaPseq and SHAPE-seq. We detected increased accessibility of the translation initiation region of psbA after high light treatment, likely contributing to the observed increase in translation by facilitating translation initiation. Furthermore, we identified the footprint of a putative regulatory protein in the 5' UTR of psbA at a position where occlusion of the nucleotide sequence would cause the structure of the translation initiation region to open up, thereby facilitating ribosome access. Moreover, we show that other plastid genes with weak Shine-Dalgarno sequences (SD) are likely to exhibit psbA-like regulation, while those with strong SDs do not. This supports the idea that changes in mRNA secondary structure might represent a general mechanism for translational regulation of psbA and other plastid genes.

Detection of internal N7-methylguanosine (m7G) RNA modifications by mutational profiling sequencing

Methylation of guanosine on position N7 (m⁷G) on internal RNA positions has been found in all domains of life and have been implicated in human disease. Here, we present m⁷G Mutational Profiling sequencing (m⁷G-MaP-seq), which allows high throughput detection of m⁷G modifications at nucleotide resolution. In our method, m⁷G modified positions are converted to abasic sites by reduction with sodium borohydride, directly recorded as cDNA mutations through reverse transcription and sequenced. We detect positions with increased mutation rates in the reduced and control samples taking the possibility of sequencing/alignment error into account and use replicates to calculate statistical significance based on log likelihood ratio tests. We show that m⁷G-MaP-seq efficiently detects known m⁷G modifications in rRNA with mutational rates up to 25% and we map a previously uncharacterised evolutionarily conserved rRNA modification at position 1581 in Arabidopsis thaliana SSU rRNA. Furthermore, we identify m⁷G modifications in budding yeast, human and arabidopsis tRNAs and demonstrate that m⁷G modification occurs before tRNA splicing. We do not find any evidence for internal m⁷G modifications being present in other small RNA, such as miRNA, snoRNA and sRNA, including human Let-7e. Likewise, high sequencing depth m⁷G-MaP-seq analysis of mRNA from E. coli or yeast cells did not identify any internal m⁷G modifications.

Time-Resolved Hydroxyl Radical Footprinting of RNA with X-Rays

RNA footprinting by hydroxyl radical cleavage provides 'snapshots' of RNA tertiary structure or protein interactions that bury the RNA backbone. Generation of hydroxyl radicals with a high-flux synchrotron X-ray beam provides analysis on a short timescale (5-100 msec), which enables the structures of folding intermediates or other transient conformational states to be determined in biochemical solutions or cells. This article provides protocols for using synchrotron beamlines for hydroxyl radical footprinting. © 2018 by John Wiley & Sons, Inc.

PROBer Provides a General Toolkit for Analyzing Sequencing-Based Toeprinting Assays

A number of sequencing-based transcriptase drop-off assays have recently been developed to probe post-transcriptional dynamics of RNA-protein interaction, RNA structure, and RNA modification. Although these assays survey a diverse set of epitranscriptomic marks, we use the term toeprinting assays since they share methodological similarities. Their interpretation is predicated on addressing a similar computational challenge: how to learn isoform-specific chemical modification profiles in the face of complex read multi-mapping. We introduce PROBer, a statistical model and associated software, that addresses this challenge for the analysis of toeprinting assays. PROBer takes sequencing data as input and outputs estimated transcript abundances and isoform-specific modification profiles. Results on both simulated and biological data demonstrate that PROBer significantly outperforms individual methods tailored for specific toeprinting assays. Since the space of toeprinting assays is ever expanding and these assays are likely to be performed and analyzed together, we believe PROBer's unified data analysis solution will be valuable to the RNA community.

High-Throughput Nuclease Probing of RNA Structures Using FragSeq

High-throughput sequencing of cDNA (RNA-Seq) can be used to generate nuclease accessibility data for many distinct transcripts in the same mixture simultaneously. Such assays accelerate RNA structure analysis and provide researchers with new technologies to tackle biological questions on a transcriptome-wide scale. FragSeq is an experimental assay for transcriptome-wide RNA structure probing using RNA-Seq, coupled with data analysis tools that allow quantitative determination of nuclease accessibility at single-base resolution. We provide a practical guide to designing and carrying out FragSeq experiments and data analysis.