Charting the unknown epitranscriptome.

RNA modifications can alter RNA structure-function relationships and various cellular processes. However, the genomic distribution and biological roles of most RNA modifications remain uncharacterized. Here, we propose using phage display antibody technology and direct sequencing through nanopores to facilitate systematic interrogation of the distribution, location and dynamics of RNA modifications.

A role for tRNA modifications in genome structure and codon usage.

Transfer RNA (tRNA) gene content is a differentiating feature of genomes that contributes to the efficiency of the translational apparatus, but the principles shaping tRNA gene copy number and codon composition are poorly understood. Here, we report that the emergence of two specific tRNA modifications shaped the structure and composition of all extant genomes. Through the analysis of more than 500 genomes, we identify two kingdom-specific tRNA modifications as major contributors that separated archaeal, bacterial, and eukaryal genomes in terms of their tRNA gene composition. We show that, contrary to prior observations, genomic codon usage and tRNA gene frequencies correlate in all kingdoms if these two modifications are taken into account and that presence or absence of these modifications explains patterns of gene expression observed in previous studies. Finally, we experimentally demonstrate that human gene expression levels correlate well with genomic codon composition if these identified modifications are considered.

Nano3P-seq: transcriptome-wide analysis of gene expression and tail dynamics using end-capture nanopore cDNA sequencing.

RNA polyadenylation plays a central role in RNA maturation, fate, and stability. In response to developmental cues, polyA tail lengths can vary, affecting the translation efficiency and stability of mRNAs. Here we develop Nanopore 3' end-capture sequencing (Nano3P-seq), a method that relies on nanopore cDNA sequencing to simultaneously quantify RNA abundance, tail composition, and tail length dynamics at per-read resolution. By employing a template-switching-based sequencing protocol, Nano3P-seq can sequence RNA molecule from its 3' end, regardless of its polyadenylation status, without the need for PCR amplification or ligation of RNA adapters. We demonstrate that Nano3P-seq provides quantitative estimates of RNA abundance and tail lengths, and captures a wide diversity of RNA biotypes. We find that, in addition to mRNA and long non-coding RNA, polyA tails can be identified in 16S mitochondrial ribosomal RNA in both mouse and zebrafish models. Moreover, we show that mRNA tail lengths are dynamically regulated during vertebrate embryogenesis at an isoform-specific level, correlating with mRNA decay. Finally, we demonstrate the ability of Nano3P-seq in capturing non-A bases within polyA tails of various lengths, and reveal their distribution during vertebrate embryogenesis. Overall, Nano3P-seq is a simple and robust method for accurately estimating transcript levels, tail lengths, and tail composition heterogeneity in individual reads, with minimal library preparation biases, both in the coding and non-coding transcriptome.

Comprehensive map of ribosomal 2'-O-methylation and C/D box snoRNAs in Drosophila melanogaster.

During their maturation, ribosomal RNAs (rRNAs) are decorated by hundreds of chemical modifications that participate in proper folding of rRNA secondary structures and therefore in ribosomal function. Along with pseudouridine, methylation of the 2'-hydroxyl ribose moiety (Nm) is the most abundant modification of rRNAs. The majority of Nm modifications in eukaryotes are placed by Fibrillarin, a conserved methyltransferase belonging to a ribonucleoprotein complex guided by C/D box small nucleolar RNAs (C/D box snoRNAs). These modifications impact interactions between rRNAs, tRNAs and mRNAs, and some are known to fine tune translation rates and efficiency. In this study, we built the first comprehensive map of Nm sites in Drosophila melanogaster rRNAs using two complementary approaches (RiboMethSeq and Nanopore direct RNA sequencing) and identified their corresponding C/D box snoRNAs by whole-transcriptome sequencing. We de novo identified 61 Nm sites, from which 55 are supported by both sequencing methods, we validated the expression of 106 C/D box snoRNAs and we predicted new or alternative rRNA Nm targets for 31 of them. Comparison of methylation level upon different stresses show only slight but specific variations, indicating that this modification is relatively stable in D. melanogaster. This study paves the way to investigate the impact of snoRNA-mediated 2'-O-methylation on translation and proteostasis in a whole organism.

The lncRNA Snhg11, a new candidate contributing to neurogenesis, plasticity, and memory deficits in Down syndrome.

Down syndrome (DS) stands as the prevalent genetic cause of intellectual disability, yet comprehensive understanding of its cellular and molecular underpinnings remains limited. In this study, we explore the cellular landscape of the hippocampus in a DS mouse model, the Ts65Dn, through single-nuclei transcriptional profiling. Our findings demonstrate that trisomy manifests as a highly specific modification of the transcriptome within distinct cell types. Remarkably, we observed a significant shift in the transcriptomic profile of granule cells in the dentate gyrus (DG) associated with trisomy. We identified the downregulation of a specific small nucleolar RNA host gene, Snhg11, as the primary driver behind this observed shift in the trisomic DG. Notably, reduced levels of Snhg11 in this region were also observed in a distinct DS mouse model, the Dp(16)1Yey, as well as in human postmortem brain tissue, indicating its relevance in Down syndrome. To elucidate the function of this long non-coding RNA (lncRNA), we knocked down Snhg11 in the DG of wild-type mice. Intriguingly, this intervention alone was sufficient to impair synaptic plasticity and adult neurogenesis, resembling the cognitive phenotypes associated with trisomy in the hippocampus. Our study uncovers the functional role of Snhg11 in the DG and underscores the significance of this lncRNA in intellectual disability. Furthermore, our findings highlight the importance of DG in the memory deficits observed in Down syndrome.

Dynamic interplay between RPL3- and RPL3L-containing ribosomes modulates mitochondrial activity in the mammalian heart.

The existence of naturally occurring ribosome heterogeneity is now a well-acknowledged phenomenon. However, whether this heterogeneity leads to functionally diverse 'specialized ribosomes' is still a controversial topic. Here, we explore the biological function of RPL3L (uL3L), a ribosomal protein (RP) paralogue of RPL3 (uL3) that is exclusively expressed in skeletal muscle and heart tissues, by generating a viable homozygous Rpl3l knockout mouse strain. We identify a rescue mechanism in which, upon RPL3L depletion, RPL3 becomes up-regulated, yielding RPL3-containing ribosomes instead of RPL3L-containing ribosomes that are typically found in cardiomyocytes. Using both ribosome profiling (Ribo-seq) and a novel orthogonal approach consisting of ribosome pulldown coupled to nanopore sequencing (Nano-TRAP), we find that RPL3L modulates neither translational efficiency nor ribosome affinity towards a specific subset of transcripts. In contrast, we show that depletion of RPL3L leads to increased ribosome-mitochondria interactions in cardiomyocytes, which is accompanied by a significant increase in ATP levels, potentially as a result of fine-tuning of mitochondrial activity. Our results demonstrate that the existence of tissue-specific RP paralogues does not necessarily lead to enhanced translation of specific transcripts or modulation of translational output. Instead, we reveal a complex cellular scenario in which RPL3L modulates the expression of RPL3, which in turn affects ribosomal subcellular localization and, ultimately, mitochondrial activity.

Ribosomes are macromolecular machines responsible for protein synthesis in all living beings. Recent studies have shown that ribosomes can be heterogeneous in their structure, possibly leading to a specialized function. Here, we focus on RPL3L, a ribosomal protein expressed exclusively in striated muscles. We find that the deletion of the Rpl3l gene in a mouse model triggers a compensation mechanism, in which the missing RPL3L protein is replaced by its paralogue, RPL3. Furthermore, we find that RPL3-containing ribosomes establish closer interactions with mitochondria, cellular organelles responsible for energy production, leading to higher energy production when compared with RPL3L-containing ribosomes. Finally, we show that the RPL3­RPL3L compensation mechanism is also triggered in heart disease conditions, such as hypertrophy and myocardial infarction.

Exploring the epitranscriptome by native RNA sequencing.

Chemical RNA modifications, collectively referred to as the "epitranscriptome," are essential players in fine-tuning gene expression. Our ability to analyze RNA modifications has improved rapidly in recent years, largely due to the advent of high-throughput sequencing methodologies, which typically consist of coupling modification-specific reagents, such as antibodies or enzymes, to next-generation sequencing. Recently, it also became possible to map RNA modifications directly by sequencing native RNAs using nanopore technologies, which has been applied for the detection of a number of RNA modifications, such as N6-methyladenosine (m6A), pseudouridine (Ψ), and inosine (I). However, the signal modulations caused by most RNA modifications are yet to be determined. A global effort is needed to determine the signatures of the full range of RNA modifications to avoid the technical biases that have so far limited our understanding of the epitranscriptome.

Spectral libraries from nucleobases and deoxyribonucleosides facilitate the identification of ribonucleosides by nano-flow liquid chromatography-tandem mass spectrometry.

RATIONALE: The study addresses the challenge of identifying RNA post-transcriptional modifications when commercial standards are not available to generate reference spectral libraries. It proposes employing homologous nucleobases and deoxyribonucleosides as alternative reference spectral libraries to aid in identifying modified ribonucleosides and distinguishing them from their positional isomers when the standards are unavailable. METHODS: Complete sets of ribonucleoside, deoxyribonucleoside and nucleobase standards were analyzed using high-performance nano-flow liquid chromatography coupled to an Orbitrap Eclipse Tribrid mass spectrometer. Spectral libraries were constructed from homologous nucleobases and deoxyribonucleosides using targeted MS2 and neutral-loss-triggered MS3 methods, and collision energies were optimized. The feasibility of using these libraries for identifying modified ribonucleosides and their positional isomers was assessed through comparison of spectral fragmentation patterns. RESULTS: Our analysis reveals that both MS2 and neutral-loss-triggered MS3 methods yielded rich spectra with similar fragmentation patterns across ribonucleosides, deoxyribonucleosides and nucleobases. Moreover, we demonstrate that spectra from nucleobases and deoxyribonucleosides, generated at optimized collision energies, exhibited sufficient similarity to those of modified ribonucleosides to enable their use as reference spectra for accurate identification of positional isomers within ribonucleoside families. CONCLUSIONS: The study demonstrates the efficacy of utilizing homologous nucleobases and deoxyribonucleosides as interchangeable reference spectral libraries for identifying modified ribonucleosides and their positional isomers. This approach offers a valuable solution for overcoming limitations posed by the unavailability of commercial standards, enhancing the analysis of RNA post-transcriptional modifications via mass spectrometry.

Molecular barcoding of native RNAs using nanopore sequencing and deep learning.

Nanopore sequencing enables direct measurement of RNA molecules without conversion to cDNA, thus opening the gates to a new era for RNA biology. However, the lack of molecular barcoding of direct RNA nanopore sequencing data sets severely affects the applicability of this technology to biological samples, where RNA availability is often limited. Here, we provide the first experimental protocol and associated algorithm to barcode and demultiplex direct RNA nanopore sequencing data sets. Specifically, we present a novel and robust approach to accurately classify raw nanopore signal data by transforming current intensities into images or arrays of pixels, followed by classification using a deep learning algorithm. We demonstrate the power of this strategy by developing the first experimental protocol for barcoding and demultiplexing direct RNA sequencing libraries. Our method, DeePlexiCon, can classify 93% of reads with 95.1% accuracy or 60% of reads with 99.9% accuracy. The availability of an efficient and simple multiplexing strategy for native RNA sequencing will improve the cost-effectiveness of this technology, as well as facilitate the analysis of lower-input biological samples. Overall, our work exemplifies the power, simplicity, and robustness of signal-to-image conversion for nanopore data analysis using deep learning.

Human tRNAs with inosine 34 are essential to efficiently translate eukarya-specific low-complexity proteins.

The modification of adenosine to inosine at the wobble position (I34) of tRNA anticodons is an abundant and essential feature of eukaryotic tRNAs. The expansion of inosine-containing tRNAs in eukaryotes followed the transformation of the homodimeric bacterial enzyme TadA, which generates I34 in tRNAArg and tRNALeu, into the heterodimeric eukaryotic enzyme ADAT, which modifies up to eight different tRNAs. The emergence of ADAT and its larger set of substrates, strongly influenced the tRNA composition and codon usage of eukaryotic genomes. However, the selective advantages that drove the expansion of I34-tRNAs remain unknown. Here we investigate the functional relevance of I34-tRNAs in human cells and show that a full complement of these tRNAs is necessary for the translation of low-complexity protein domains enriched in amino acids cognate for I34-tRNAs. The coding sequences for these domains require codons translated by I34-tRNAs, in detriment of synonymous codons that use other tRNAs. I34-tRNA-dependent low-complexity proteins are enriched in functional categories related to cell adhesion, and depletion in I34-tRNAs leads to cellular phenotypes consistent with these roles. We show that the distribution of these low-complexity proteins mirrors the distribution of I34-tRNAs in the phylogenetic tree.

ModPhred: an integrative toolkit for the analysis and storage of nanopore sequencing DNA and RNA modification data.

MOTIVATION: DNA and RNA modifications can now be identified using nanopore sequencing. However, we currently lack a flexible software to efficiently encode, store, analyze and visualize DNA and RNA modification data. RESULTS: Here, we present ModPhred, a versatile toolkit that facilitates DNA and RNA modification analysis from nanopore sequencing reads in a user-friendly manner. ModPhred integrates probabilistic DNA and RNA modification information within the FASTQ and BAM file formats, can be used to encode multiple types of modifications simultaneously, and its output can be easily coupled to genomic track viewers, facilitating the visualization and analysis of DNA and RNA modification information in individual reads in a simple and computationally efficient manner. AVAILABILITY AND IMPLEMENTATION: ModPhred is available at https://github.com/novoalab/modPhred, is implemented in Python3, and is released under an MIT license. Docker images with all dependencies preinstalled are also provided. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

Subcellular relocalization and nuclear redistribution of the RNA methyltransferases TRMT1 and TRMT1L upon neuronal activation.

RNA modifications are dynamic chemical entities that expand the RNA lexicon and regulate RNA fate. The most abundant modification present in mRNAs, N6-methyladenosine (m6A), has been implicated in neurogenesis and memory formation. However, whether additional RNA modifications may be playing a role in neuronal functions and in response to environmental queues is largely unknown. Here we characterize the biochemical function and cellular dynamics of two human RNA methyltransferases previously associated with neurological dysfunction, TRMT1 and its homolog, TRMT1-like (TRMT1L). Using a combination of next-generation sequencing, LC-MS/MS, patient-derived cell lines and knockout mouse models, we confirm the previously reported dimethylguanosine (m2,2G) activity of TRMT1 in tRNAs, as well as reveal that TRMT1L, whose activity was unknown, is responsible for methylating a subset of cytosolic tRNAAla(AGC) isodecoders at position 26. Using a cellular in vitro model that mimics neuronal activation and long term potentiation, we find that both TRMT1 and TRMT1L change their subcellular localization upon neuronal activation. Specifically, we observe a major subcellular relocalization from mitochondria and other cytoplasmic domains (TRMT1) and nucleoli (TRMT1L) to different small punctate compartments in the nucleus, which are as yet uncharacterized. This phenomenon does not occur upon heat shock, suggesting that the relocalization of TRMT1 and TRMT1L is not a general reaction to stress, but rather a specific response to neuronal activation. Our results suggest that subcellular relocalization of RNA modification enzymes may play a role in neuronal plasticity and transmission of information, presumably by addressing new targets.

Elucidation of Codon Usage Signatures across the Domains of Life.

Because of the degeneracy of the genetic code, multiple codons are translated into the same amino acid. Despite being "synonymous," these codons are not equally used. Selective pressures are thought to drive the choice among synonymous codons within a genome, while GC content, which is typically attributed to mutational drift, is the major determinant of variation across species. Here, we find that in addition to GC content, interspecies codon usage signatures can also be detected. More specifically, we show that a single amino acid, arginine, is the major contributor to codon usage bias differences across domains of life. We then exploit this finding and show that domain-specific codon bias signatures can be used to classify a given sequence into its corresponding domain of life with high accuracy. We then wondered whether the inclusion of codon usage codon autocorrelation patterns, which reflects the nonrandom distribution of codon occurrences throughout a transcript, might improve the classification performance of our algorithm. However, we find that autocorrelation patterns are not domain-specific, and surprisingly, are unrelated to tRNA reusage, in contrast to previous reports. Instead, our results suggest that codon autocorrelation patterns are a by-product of codon optimality throughout a sequence, where highly expressed genes display autocorrelated "optimal" codons, whereas lowly expressed genes display autocorrelated "nonoptimal" codons.

The RNA modification landscape in human disease.

RNA modifications have been historically considered as fine-tuning chemo-structural features of infrastructural RNAs, such as rRNAs, tRNAs, and snoRNAs. This view has changed dramatically in recent years, to a large extent as a result of systematic efforts to map and quantify various RNA modifications in a transcriptome-wide manner, revealing that RNA modifications are reversible, dynamically regulated, far more widespread than originally thought, and involved in major biological processes, including cell differentiation, sex determination, and stress responses. Here we summarize the state of knowledge and provide a catalog of RNA modifications and their links to neurological disorders, cancers, and other diseases. With the advent of direct RNA-sequencing technologies, we expect that this catalog will help prioritize those RNA modifications for transcriptome-wide maps.

Homologous trans-editing factors with broad tRNA specificity prevent mistranslation caused by serine/threonine misactivation.

Aminoacyl-tRNA synthetases (ARSs) establish the rules of the genetic code, whereby each amino acid is attached to a cognate tRNA. Errors in this process lead to mistranslation, which can be toxic to cells. The selective forces exerted by species-specific requirements and environmental conditions potentially shape quality-control mechanisms that serve to prevent mistranslation. A family of editing factors that are homologous to the editing domain of bacterial prolyl-tRNA synthetase includes the previously characterized trans-editing factors ProXp-ala and YbaK, which clear Ala-tRNA(Pro) and Cys-tRNA(Pro), respectively, and three additional homologs of unknown function, ProXp-x, ProXp-y, and ProXp-z. We performed an in vivo screen of 230 conditions in which an Escherichia coli proXp-y deletion strain was grown in the presence of elevated levels of amino acids and specific ARSs. This screen, together with the results of in vitro deacylation assays, revealed Ser- and Thr-tRNA deacylase function for this homolog. A similar activity was demonstrated for Bordetella parapertussis ProXp-z in vitro. These proteins, now renamed "ProXp-ST1" and "ProXp-ST2," respectively, recognize multiple tRNAs as substrates. Taken together, our data suggest that these free-standing editing domains have the ability to prevent mistranslation errors caused by a number of ARSs, including lysyl-tRNA synthetase, threonyl-tRNA synthetase, seryl-tRNA synthetase, and alanyl-tRNA synthetase. The expression of these multifunctional enzymes is likely to provide a selective growth advantage to organisms subjected to environmental stresses and other conditions that alter the amino acid pool.

SwiSpot: modeling riboswitches by spotting out switching sequences.

MOTIVATION: Riboswitches are cis-regulatory elements in mRNA, mostly found in Bacteria, which exhibit two main secondary structure conformations. Although one of them prevents the gene from being expressed, the other conformation allows its expression, and this switching process is typically driven by the presence of a specific ligand. Although there are a handful of known riboswitches, our knowledge in this field has been greatly limited due to our inability to identify their alternate structures from their sequences. Indeed, current methods are not able to predict the presence of the two functionally distinct conformations just from the knowledge of the plain RNA nucleotide sequence. Whether this would be possible, for which cases, and what prediction accuracy can be achieved, are currently open questions. RESULTS: Here we show that the two alternate secondary structures of riboswitches can be accurately predicted once the 'switching sequence' of the riboswitch has been properly identified. The proposed SwiSpot approach is capable of identifying the switching sequence inside a putative, complete riboswitch sequence, on the basis of pairing behaviors, which are evaluated on proper sets of configurations. Moreover, it is able to model the switching behavior of riboswitches whose generated ensemble covers both alternate configurations. Beyond structural predictions, the approach can also be paired to homology-based riboswitch searches. AVAILABILITY AND IMPLEMENTATION: SwiSpot software, along with the reference dataset files, is available at: http://www.iet.unipi.it/a.bechini/swispot/Supplementary information: Supplementary data are available at Bioinformatics online. CONTACT: a.bechini@ing.unipi.it.

Analogs of natural aminoacyl-tRNA synthetase inhibitors clear malaria in vivo.

Malaria remains a major global health problem. Emerging resistance to existing antimalarial drugs drives the search for new antimalarials, and protein translation is a promising pathway to target. Here we explore the potential of the aminoacyl-tRNA synthetase (ARS) family as a source of antimalarial drug targets. First, a battery of known and novel ARS inhibitors was tested against Plasmodium falciparum cultures, and their activities were compared. Borrelidin, a natural inhibitor of threonyl-tRNA synthetase (ThrRS), stands out for its potent antimalarial effect. However, it also inhibits human ThrRS and is highly toxic to human cells. To circumvent this problem, we tested a library of bioengineered and semisynthetic borrelidin analogs for their antimalarial activity and toxicity. We found that some analogs effectively lose their toxicity against human cells while retaining a potent antiparasitic activity both in vitro and in vivo and cleared malaria from Plasmodium yoelii-infected mice, resulting in 100% mice survival rates. Our work identifies borrelidin analogs as potent, selective, and unexplored scaffolds that efficiently clear malaria both in vitro and in vivo.

Ancestral AlaX editing enzymes for control of genetic code fidelity are not tRNA-specific.

Accurate protein synthesis requires the hydrolytic editing of tRNAs incorrectly aminoacylated by aminoacyl-tRNA synthetases (ARSs). Recognition of cognate tRNAs by ARS is less error-prone than amino acid recognition, and, consequently, editing domains are generally believed to act only on the tRNAs cognate to their related ARSs. For example, the AlaX family of editing domains, including the editing domain of alanyl-tRNA synthetase and the related free-standing trans-editing AlaX enzymes, are thought to specifically act on tRNA(Ala), whereas the editing domains of threonyl-tRNA synthetases are specific for tRNA(Thr). Here we show that, contrary to this belief, AlaX-S, the smallest of the extant AlaX enzymes, deacylates Ser-tRNA(Thr) in addition to Ser-tRNA(Ala) and that a single residue is important to determine this behavior. Our data indicate that promiscuous forms of AlaX are ancestral to tRNA-specific AlaXs. We propose that former AlaX domains were used to maintain translational fidelity in earlier stages of genetic code evolution when mis-serylation of several tRNAs was possible.