Benchmarking of computational methods for m6A profiling with Nanopore direct RNA sequencing.

N6-methyladenosine (m6A) is the most abundant internal eukaryotic mRNA modification, and is involved in the regulation of various biological processes. Direct Nanopore sequencing of native RNA (dRNA-seq) emerged as a leading approach for its identification. Several software were published for m6A detection and there is a strong need for independent studies benchmarking their performance on data from different species, and against various reference datasets. Moreover, a computational workflow is needed to streamline the execution of tools whose installation and execution remains complicated. We developed NanOlympicsMod, a Nextflow pipeline exploiting containerized technology for comparing 14 tools for m6A detection on dRNA-seq data. NanOlympicsMod was tested on dRNA-seq data generated from in vitro (un)modified synthetic oligos. The m6A hits returned by each tool were compared to the m6A position known by design of the oligos. In addition, NanOlympicsMod was used on dRNA-seq datasets from wild-type and m6A-depleted yeast, mouse and human, and each tool's hits were compared to reference m6A sets generated by leading orthogonal methods. The performance of the tools markedly differed across datasets, and methods adopting different approaches showed different preferences in terms of precision and recall. Changing the stringency cut-offs allowed for tuning the precision-recall trade-off towards user preferences. Finally, we determined that precision and recall of tools are markedly influenced by sequencing depth, and that additional sequencing would likely reveal additional m6A sites. Thanks to the possibility of including novel tools, NanOlympicsMod will streamline the benchmarking of m6A detection tools on dRNA-seq data, improving future RNA modification characterization.

Comprehensive transcriptome analysis reveals altered mRNA splicing and post-transcriptional changes in the aged mouse brain.

A comprehensive understanding of molecular changes during brain aging is essential to mitigate cognitive decline and delay neurodegenerative diseases. The interpretation of mRNA alterations during brain aging is influenced by the health and age of the animal cohorts studied. Here, we carefully consider these factors and provide an in-depth investigation of mRNA splicing and dynamics in the aging mouse brain, combining short- and long-read sequencing technologies with extensive bioinformatic analyses. Our findings encompass a spectrum of age-related changes, including differences in isoform usage, decreased mRNA dynamics and a module showing increased expression of neuronal genes. Notably, our results indicate a reduced abundance of mRNA isoforms leading to nonsense-mediated RNA decay and suggest a regulatory role for RNA-binding proteins, indicating that their regulation may be altered leading to the reshaping of the aged brain transcriptome. Collectively, our study highlights the importance of studying mRNA splicing events during brain aging.

Genome-wide dynamics of RNA synthesis, processing, and degradation without RNA metabolic labeling.

The quantification of the kinetic rates of RNA synthesis, processing, and degradation are largely based on the integrative analysis of total and nascent transcription, the latter being quantified through RNA metabolic labeling. We developed INSPEcT-, a computational method based on the mathematical modeling of premature and mature RNA expression that is able to quantify kinetic rates from steady-state or time course total RNA-seq data without requiring any information on nascent transcripts. Our approach outperforms available solutions, closely recapitulates the kinetic rates obtained through RNA metabolic labeling, improves the ability to detect changes in transcript half-lives, reduces the cost and complexity of the experiments, and can be adopted to study experimental conditions in which nascent transcription cannot be readily profiled. Finally, we applied INSPEcT- to the characterization of post-transcriptional regulation landscapes in dozens of physiological and disease conditions. This approach was included in the INSPEcT Bioconductor package, which can now unveil RNA dynamics from steady-state or time course data, with or without the profiling of nascent RNA.

Dynamics of transcriptional and post-transcriptional regulation.

Despite gene expression programs being notoriously complex, RNA abundance is usually assumed as a proxy for transcriptional activity. Recently developed approaches, able to disentangle transcriptional and post-transcriptional regulatory processes, have revealed a more complex scenario. It is now possible to work out how synthesis, processing and degradation kinetic rates collectively determine the abundance of each gene's RNA. It has become clear that the same transcriptional output can correspond to different combinations of the kinetic rates. This underscores the fact that markedly different modes of gene expression regulation exist, each with profound effects on a gene's ability to modulate its own expression. This review describes the development of the experimental and computational approaches, including RNA metabolic labeling and mathematical modeling, that have been disclosing the mechanisms underlying complex transcriptional programs. Current limitations and future perspectives in the field are also discussed.

An early Myc-dependent transcriptional program orchestrates cell growth during B-cell activation.

Upon activation, lymphocytes exit quiescence and undergo substantial increases in cell size, accompanied by activation of energy-producing and anabolic pathways, widespread chromatin decompaction, and elevated transcriptional activity. These changes depend upon prior induction of the Myc transcription factor, but how Myc controls them remains unclear. We addressed this issue by profiling the response to LPS stimulation in wild-type and c-myc-deleted primary mouse B-cells. Myc is rapidly induced, becomes detectable on virtually all active promoters and enhancers, but has no direct impact on global transcriptional activity. Instead, Myc contributes to the swift up- and down-regulation of several hundred genes, including many known regulators of the aforementioned cellular processes. Myc-activated promoters are enriched for E-box consensus motifs, bind Myc at the highest levels, and show enhanced RNA Polymerase II recruitment, the opposite being true at down-regulated loci. Remarkably, the Myc-dependent signature identified in activated B-cells is also enriched in Myc-driven B-cell lymphomas: hence, besides modulation of new cancer-specific programs, the oncogenic action of Myc may largely rely on sustained deregulation of its normal physiological targets.

Computational methods for RNA modification detection from nanopore direct RNA sequencing data.

The covalent modification of RNA molecules is a pervasive feature of all classes of RNAs and has fundamental roles in the regulation of several cellular processes. Mapping the location of RNA modifications transcriptome-wide is key to unveiling their role and dynamic behaviour, but technical limitations have often hampered these efforts. Nanopore direct RNA sequencing is a third-generation sequencing technology that allows the sequencing of native RNA molecules, thus providing a direct way to detect modifications at single-molecule resolution. Despite recent advances, the analysis of nanopore sequencing data for RNA modification detection is still a complex task that presents many challenges. Many works have addressed this task using different approaches, resulting in a large number of tools with different features and performances. Here we review the diverse approaches proposed so far and outline the principles underlying currently available algorithms.

Depletion of Mettl3 in cholinergic neurons causes adult-onset neuromuscular degeneration.

Motor neuron (MN) demise is a hallmark of several neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS). Post-transcriptional gene regulation can control RNA's fate, and defects in RNA processing are critical determinants of MN degeneration. N6-methyladenosine (m6A) is a post-transcriptional RNA modification that controls diverse aspects of RNA metabolism. To assess the m6A requirement in MNs, we depleted the m6A methyltransferase-like 3 (METTL3) in cells and mice. METTL3 depletion in embryonic stem cell-derived MNs has profound and selective effects on survival and neurite outgrowth. Mice with cholinergic neuron-specific METTL3 depletion display a progressive decline in motor behavior, accompanied by MN loss and muscle denervation, culminating in paralysis and death. Reader proteins convey m6A effects, and their silencing phenocopies METTL3 depletion. Among the m6A targets, we identified transactive response DNA-binding protein 43 (TDP-43) and discovered that its expression is under epitranscriptomic control. Thus, impaired m6A signaling disrupts MN homeostasis and triggers neurodegeneration conceivably through TDP-43 deregulation.

News from around the RNA world: new avenues in RNA biology, biotechnology and therapeutics from the 2022 SIBBM meeting.

Since the formalization of the Central Dogma of molecular biology, the relevance of RNA in modulating the flow of information from DNA to proteins has been clear. More recently, the discovery of a vast set of non-coding transcripts involved in crucial aspects of cellular biology has renewed the enthusiasm of the RNA community. Moreover, the remarkable impact of RNA therapies in facing the COVID19 pandemics has bolstered interest in the translational opportunities provided by this incredible molecule. For all these reasons, the Italian Society of Biophysics and Molecular Biology (SIBBM) decided to dedicate its 17th yearly meeting, held in June 2022 in Rome, to the many fascinating aspects of RNA biology. More than thirty national and international speakers covered the properties, modes of action and applications of RNA, from its role in the control of development and cell differentiation to its involvement in disease. Here, we summarize the scientific content of the conference, highlighting the take-home message of each presentation, and we stress the directions the community is currently exploring to push forward our comprehension of the RNA World 3.0.

Identification of Genes Post-Transcriptionally Regulated from RNA-seq: The Case Study of Liver Hepatocellular Carcinoma.

The field of transcriptional regulation generally assumes that changes in transcripts levels reflect changes in transcriptional status of the corresponding gene. While this assumption might hold true for a large population of transcripts, a considerable and still unrecognized fraction of the variation might involve other steps of the RNA lifecycle, that is the processing of the premature RNA, and degradation of the mature RNA. Discrimination between these layers requires complementary experimental techniques, such as RNA metabolic labeling or block of transcription experiments. Nonetheless, the analysis of the premature and mature RNA, derived from intronic and exonic read counts in RNA-seq data, allows distinguishing between transcriptionally and post-transcriptionally regulated genes, although not recognizing the specific step involved in the post-transcriptional response, that is processing, degradation, or a combination of the two. We illustrate how the INSPEcT R/Bioconductor package could be used to infer post-transcriptional regulation in TCGA RNA-seq samples for Hepatocellular Carcinoma.

Sharing the Space With the "Victim" Can Increase Help Rates. A Study With Virtual Reality.

A typical protocol for the psychological study of helping behavior features two core roles: a help seeker suffering from some personal or situational emergency (often called "victim") and a potential helper. The setting of these studies is such that the victim and the helper often share the same space. We wondered whether this spatial arrangement might affect the help rate. Thus, we designed a simple study with virtual reality in which space sharing could be manipulated. The participant plays the role of a potential helper; the victim is a humanoid located inside the virtual building. When the request for help is issued, the participant can be either in the same spatial region as the victim (the virtual building) or outside it. The effect of space was tested in two kinds of emergencies: a mere request for help and a request for help during a fire. The analysis shows that, in both kinds of emergencies, the participants were more likely to help the victim when sharing the space with it. This study suggests controlling the spatial arrangement when investigating helping behavior. It also illustrates the expediency of virtual reality to further investigate the role of space on pro-social behavior during emergencies.

INSPEcT-GUI Reveals the Impact of the Kinetic Rates of RNA Synthesis, Processing, and Degradation, on Premature and Mature RNA Species.

The abundance of RNA species and their response to perturbations are set by the kinetics rates of RNA synthesis, processing, and degradation. However, the visualization, interpretation, and manipulation of these data require familiarity with mathematical modeling and command line tools. INSPEcT-GUI is an R-Shiny interface that allows researchers without specific training to effortlessly explore how the fine kinetic regulation of the RNA life cycle can shape gene expression programs. In particular, it allows to: (i) interactively visualize gene-level RNA dynamics; (ii) refine the model fit of experimental data; (iii) test alternative regulatory models; (iv) explore, independently from the availability of data, how the combined action of the RNA kinetic rates impacts on premature and mature RNA. INSPEcT-GUI is freely available within the R/Bioconductor package INSPEcT at http://bioconductor.org/packages/INSPEcT/. An HTML vignette including documentation on the tool startup and usage, executable examples, and a video demonstration, are available at: http://bioconductor.org/packages/release/bioc/vignettes/INSPEcT/inst/doc/INSPEcT_GUI.html.

Direct RNA Sequencing for the Study of Synthesis, Processing, and Degradation of Modified Transcripts.

It has been known for a few decades that transcripts can be marked by dozens of different modifications. Yet, we are just at the beginning of charting these marks and understanding their functional impact. High-quality methods were developed for the profiling of some of these marks, and approaches to finely study their impact on specific phases of the RNA life-cycle are available, including RNA metabolic labeling. Thanks to these improvements, the most abundant marks, including N6-methyladenosine, are emerging as important determinants of the fate of marked RNAs. However, we still lack approaches to directly study how the set of marks for a given RNA molecule shape its fate. In this perspective, we first review current leading approaches in the field. Then, we propose an experimental and computational setup, based on direct RNA sequencing and mathematical modeling, to decipher the functional consequences of RNA modifications on the fate of individual RNA molecules and isoforms.

Anti-Inflammatory Activity of Exopolysaccharides from Phormidium sp. ETS05, the Most Abundant Cyanobacterium of the Therapeutic Euganean Thermal Muds, Using the Zebrafish Model.

The Euganean Thermal District (Italy) represents the oldest and largest thermal center in Europe, and its therapeutic mud is considered a unique product whose beneficial effects have been documented since Ancient Roman times. Mud properties depend on the heat and electrolytes of the thermal water, as well as on the bioactive molecules produced by its biotic component, mainly represented by cyanobacteria. The investigation of the healing effects of compounds produced by the Euganean cyanobacteria represents an important goal for scientific validation of Euganean mud therapies and for the discovering of new health beneficial biomolecules. In this work, we evaluated the therapeutic potential of exopolysaccharides (EPS) produced by Phormidium sp. ETS05, the most abundant cyanobacterium of the Euganean mud. Specifically, Phormidium EPS resulted in exerting anti-inflammatory and pro-resolution activities in chemical and injury-induced zebrafish inflammation models as demonstrated using specific transgenic zebrafish lines and morphometric and expression analyses. Moreover, in vivo and in vitro tests showed no toxicity at all for the EPS concentrations tested. The results suggest that these EPS, with their combined anti-inflammatory and pro-resolution activities, could be one of the most important therapeutic molecules present in the Euganean mud and confirm the potential of these treatments for chronic inflammatory disease recovery.

m6A-Dependent RNA Dynamics in T Cell Differentiation.

N6-methyladenosine (m6A) is the most abundant RNA modification. It has been involved in the regulation of RNA metabolism, including degradation and translation, in both physiological and disease conditions. A recent study showed that m6A-mediated degradation of key transcripts also plays a role in the control of T cells homeostasis and IL-7 induced differentiation. We re-analyzed the omics data from that study and, through the integrative analysis of total and nascent RNA-seq data, we were able to comprehensively quantify T cells RNA dynamics and how these are affected by m6A depletion. In addition to the expected impact on RNA degradation, we revealed a broader effect of m6A on RNA dynamics, which included the alteration of RNA synthesis and processing. Altogether, the combined action of m6A on all major steps of the RNA life-cycle closely re-capitulated the observed changes in the abundance of premature and mature RNA species. Ultimately, our re-analysis extended the findings of the initial study, focused on RNA stability, and proposed a yet unappreciated role for m6A in RNA synthesis and processing dynamics.