Biogenesis of Circular RNAs.

Circular RNAs are generated during splicing through various mechanisms. Ashwal-Fluss et al. demonstrate that exon circularization and linear splicing compete with each other in a tissue-specific fashion, and Zhang et al. show that exon circularization depends on flanking intronic complementary sequences. Both papers show that several types of circular RNA transcripts can be produced from a single gene.

Contribution of tRNA sequence and modifications to the decoding preferences of E. coli and M. mycoides tRNAGlyUCC for synonymous glycine codons.

tRNA superwobbling, used by certain bacteria and organelles, is an intriguing decoding concept in which a single tRNA isoacceptor is used to decode all synonymous codons of a four-fold degenerate codon box. While Escherichia coli relies on three tRNAGly isoacceptors to decode the four glycine codons (GGN), Mycoplasma mycoides requires only a single tRNAGly. Both organisms express tRNAGly with the anticodon UCC, which are remarkably similar in sequence but different in their decoding ability. By systematically introducing mutations and altering the number and type of tRNA modifications using chemically synthesized tRNAs, we elucidated the contribution of individual nucleotides and chemical groups to decoding by the E. coli and M. mycoides tRNAGly. The tRNA sequence was identified as the key factor for superwobbling, revealing the T-arm sequence as a novel pivotal element. In addition, the presence of tRNA modifications, although not essential for providing superwobbling, was shown to delicately fine-tune and balance the decoding of synonymous codons. This emphasizes that the tRNA sequence and its modifications together form an intricate system of high complexity that is indispensable for accurate and efficient decoding.

SilkMeta: a comprehensive platform for sharing and exploiting pan-genomic and multi-omic silkworm data.

The silkworm Bombyx mori is a domesticated insect that serves as an animal model for research and agriculture. The silkworm super-pan-genome dataset, which we published last year, is a unique resource for the study of global genomic diversity and phenotype-genotype association. Here we present SilkMeta (http://silkmeta.org.cn), a comprehensive database covering the available silkworm pan-genome and multi-omics data. The database contains 1082 short-read genomes, 546 long-read assembled genomes, 1168 transcriptomes, 294 phenotype characterizations (phenome), tens of millions of variations (variome), 7253 long non-coding RNAs (lncRNAs), 18 717 full length transcripts and a set of population statistics. We have compiled publications on functional genomics research and genetic stock deciphering (mutant map). A range of bioinformatics tools is also provided for data visualization and retrieval. The large batch of omics data and tools were integrated in twelve functional modules that provide useful strategies and data for comparative and functional genomics research. The interactive bioinformatics platform SilkMeta will benefit not only the silkworm but also the insect biology communities.

Identification of Up47 in three thermophilic archaea, one mesophilic archaeon, and one hyperthermophilic bacterium.

Analysis of the profile of the tRNA modifications in several Archaea allowed us to observe a novel modified uridine in the V-loop of several tRNAs from two species: Pyrococcus furiosus and Sulfolobus acidocaldarius Recently, Ohira and colleagues characterized 2'-phosphouridine (Up) at position 47 in tRNAs of thermophilic Sulfurisphaera tokodaii, as well as in several other archaea and thermophilic bacteria. From the presence of the gene arkI corresponding to the RNA kinase responsible for Up47 formation, they also concluded that Up47 should be present in tRNAs of other thermophilic Archaea Reanalysis of our earlier data confirms that the unidentified residue in tRNAs of both P. furiosus and S. acidocaldarius is indeed 2'-phosphouridine followed by m5C48. Moreover, we find this modification in several tRNAs of other Archaea and of the hyperthermophilic bacterium Aquifex aeolicus.

Anionic G•U pairs in bacterial ribosomal rRNAs.

Wobble GU pairs (or G•U) occur frequently within double-stranded RNA helices interspersed between standard G=C and A-U Watson-Crick pairs. Another type of G•U pair interacting via their Watson-Crick edges has been observed in the A site of ribosome structures between a modified U34 in the tRNA anticodon triplet and G + 3 in the mRNA. In such pairs, the electronic structure of the U is changed with a negative charge on N3(U), resulting in two H-bonds between N1(G)…O4(U) and N2(G)…N3(U). Here, we report that such pairs occur in other highly conserved positions in ribosomal RNAs of bacteria in the absence of U modification. An anionic cis Watson-Crick G•G pair is also observed and well conserved in the small subunit. These pairs are observed in tightly folded regions.

SARS-CoV-2 NSP1 induces mRNA cleavages on the ribosome.

In severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the non-structural protein NSP1 inhibits translation of host mRNAs by binding to the mRNA entry channel of the ribosome and, together with the 5'-untranslated region (UTR) of the viral mRNAs, allows the evasion of that inhibition. Here, we show that NSP1 mediates endonucleolytic cleavages of both host and viral mRNAs in the 5'UTR, but with different cleavage patterns. The first pattern is observed in host mRNAs with cleavages interspersed regularly and close to the 5' cap (6-11 nt downstream of the cap). Those cleavage positions depend more on the position relative to the 5' cap than on the sequence itself. The second cleavage pattern occurs at high NSP1 concentrations and only in SARS-CoV-2 RNAs, with the cleavages clustered at positions 45, 46 and 49. Both patterns of cleavage occur with the mRNA and NSP1 bound to the ribosome, with the SL1 hairpin at the 5' end sufficient to protect from NSP1-mediated degradation at low NSP1 concentrations. We show further that the N-terminal domain of NSP1 is necessary and sufficient for efficient cleavage. We suggest that in the ribosome-bound NSP1 protein the catalytic residues of the N-terminal domain are unmasked by the remodelling of the α1- and α2-helices of the C-terminal domain.

Mistranslation of the genetic code by a new family of bacterial transfer RNAs.

The correct coupling of amino acids with transfer RNAs (tRNAs) is vital for translating genetic information into functional proteins. Errors during this process lead to mistranslation, where a codon is translated using the wrong amino acid. While unregulated and prolonged mistranslation is often toxic, growing evidence suggests that organisms, from bacteria to humans, can induce and use mistranslation as a mechanism to overcome unfavorable environmental conditions. Most known cases of mistranslation are caused by translation factors with poor substrate specificity or when substrate discrimination is sensitive to molecular changes such as mutations or posttranslational modifications. Here we report two novel families of tRNAs, encoded by bacteria from the Streptomyces and Kitasatospora genera, that adopted dual identities by integrating the anticodons AUU (for Asn) or AGU (for Thr) into the structure of a distinct proline tRNA. These tRNAs are typically encoded next to a full-length or truncated version of a distinct isoform of bacterial-type prolyl-tRNA synthetase. Using two protein reporters, we showed that these tRNAs translate asparagine and threonine codons with proline. Moreover, when expressed in Escherichia coli, the tRNAs cause varying growth defects due to global Asn-to-Pro and Thr-to-Pro mutations. Yet, proteome-wide substitutions of Asn with Pro induced by tRNA expression increased cell tolerance to the antibiotic carbenicillin, indicating that Pro mistranslation can be beneficial under certain conditions. Collectively, our results significantly expand the catalog of organisms known to possess dedicated mistranslation machinery and support the concept that mistranslation is a mechanism for cellular resiliency against environmental stress.

Evaluation of the stereochemical quality of predicted RNA 3D models in the RNA-Puzzles submissions.

In silico prediction is a well-established approach to derive a general shape of an RNA molecule based on its sequence or secondary structure. This paper reports an analysis of the stereochemical quality of the RNA three-dimensional models predicted using dedicated computer programs. The stereochemistry of 1052 RNA 3D structures, including 1030 models predicted by fully automated and human-guided approaches within 22 RNA-Puzzles challenges and reference structures, is analyzed. The evaluation is based on standards of RNA stereochemistry that the Protein Data Bank requires from deposited experimental structures. Deviations from standard bond lengths and angles, planarity, or chirality are quantified. A reduction in the number of such deviations should help in the improvement of RNA 3D structure modeling approaches.

Correlated sequence signatures are present within the genomic 5'UTR RNA and NSP1 protein in coronaviruses.

The 5'UTR part of coronavirus genomes plays key roles in the viral replication cycle and translation of viral mRNAs. The first 75-80 nt, also called the leader sequence, are identical for genomic mRNA and subgenomic mRNAs. Recently, it was shown that cooperative actions of a 5'UTR segment and the nonstructural protein NSP1 are essential for both the inhibition of host mRNAs and for specific translation of viral mRNAs. Here, sequence analyses of both the 5'UTR RNA segment and the NSP1 protein have been done for several coronaviruses, with special attention to the betacoronaviruses. The conclusions are: (i) precise specific molecular signatures can be found in both the RNA and the NSP1 protein; (ii) both types of signatures correlate between each other. Indeed, definite sequence motifs in the RNA correlate with sequence motifs in the protein, indicating a coevolution between the 5'UTR and NSP1 in betacoronaviruses. Experimental mutational data on 5'UTR and NSP1 from SARS-CoV-2 using cell-free translation extracts support these conclusions and show that some conserved key residues in the amino-terminal half of the NSP1 protein are essential for evasion to the inhibitory effect of NSP1 on translation.

The dynamic landscapes of RNA architecture.

A wealth of information on RNA folding and ribonucleoprotein assembly has emerged from analyses of structures and from the use of innovative biophysical tools. Although integrating data obtained from static structures with dynamic measurements presents major challenges, such efforts are opening new vistas on the RNA folding landscape.

Eukaryotic tRNA sequences present conserved and amino acid-specific structural signatures.

Metazoan organisms have many tRNA genes responsible for decoding amino acids. The set of all tRNA genes can be grouped in sets of common amino acids and isoacceptor tRNAs that are aminoacylated by corresponding aminoacyl-tRNA synthetases. Analysis of tRNA alignments shows that, despite the high number of tRNA genes, specific tRNA sequence motifs are highly conserved across multicellular eukaryotes. The conservation often extends throughout the isoacceptors and isodecoders with, in some cases, two sets of conserved isodecoders. This study is focused on non-Watson-Crick base pairs in the helical stems, especially GoU pairs. Each of the four helical stems may contain one or more conserved GoU pairs. Some are amino acid specific and could represent identity elements for the cognate aminoacyl tRNA synthetases. Other GoU pairs are found in more than a single amino acid and could be critical for native folding of the tRNAs. Interestingly, some GoU pairs are anticodon-specific, and others are found in phylogenetically-specific clades. Although the distribution of conservation likely reflects a balance between accommodating isotype-specific functions as well as those shared by all tRNAs essential for ribosomal translation, such conservations may indicate the existence of specialized tRNAs for specific translation targets, cellular conditions, or alternative functions.

Bacterial translation machinery for deliberate mistranslation of the genetic code.

Inaccurate expression of the genetic code, also known as mistranslation, is an emerging paradigm in microbial studies. Growing evidence suggests that many microbial pathogens can deliberately mistranslate their genetic code to help invade a host or evade host immune responses. However, discovering different capacities for deliberate mistranslation remains a challenge because each group of pathogens typically employs a unique mistranslation mechanism. In this study, we address this problem by studying duplicated genes of aminoacyl-transfer RNA (tRNA) synthetases. Using bacterial prolyl-tRNA synthetase (ProRS) genes as an example, we identify an anomalous ProRS isoform, ProRSx, and a corresponding tRNA, tRNAProA, that are predominately found in plant pathogens from Streptomyces species. We then show that tRNAProA has an unusual hybrid structure that allows this tRNA to mistranslate alanine codons as proline. Finally, we provide biochemical, genetic, and mass spectrometric evidence that cells which express ProRSx and tRNAProA can translate GCU alanine codons as both alanine and proline. This dual use of alanine codons creates a hidden proteome diversity due to stochastic Ala→Pro mutations in protein sequences. Thus, we show that important plant pathogens are equipped with a tool to alter the identity of their sense codons. This finding reveals the initial example of a natural tRNA synthetase/tRNA pair for dedicated mistranslation of sense codons.

Assessment of three-dimensional RNA structure prediction in CASP15.

The prediction of RNA three-dimensional structures remains an unsolved problem. Here, we report assessments of RNA structure predictions in CASP15, the first CASP exercise that involved RNA structure modeling. Forty-two predictor groups submitted models for at least one of twelve RNA-containing targets. These models were evaluated by the RNA-Puzzles organizers and, separately, by a CASP-recruited team using metrics (GDT, lDDT) and approaches (Z-score rankings) initially developed for assessment of proteins and generalized here for RNA assessment. The two assessments independently ranked the same predictor groups as first (AIchemy_RNA2), second (Chen), and third (RNAPolis and GeneSilico, tied); predictions from deep learning approaches were significantly worse than these top ranked groups, which did not use deep learning. Further analyses based on direct comparison of predicted models to cryogenic electron microscopy (cryo-EM) maps and x-ray diffraction data support these rankings. With the exception of two RNA-protein complexes, models submitted by CASP15 groups correctly predicted the global fold of the RNA targets. Comparisons of CASP15 submissions to designed RNA nanostructures as well as molecular replacement trials highlight the potential utility of current RNA modeling approaches for RNA nanotechnology and structural biology, respectively. Nevertheless, challenges remain in modeling fine details such as noncanonical pairs, in ranking among submitted models, and in prediction of multiple structures resolved by cryo-EM or crystallography.

Neocles B. Leontis (1955 - 2020).

None.

Comparative study on tertiary contacts and folding of RNase P RNAs from a psychrophilic, a mesophilic/radiation-resistant, and a thermophilic bacterium.

In most bacterial type A RNase P RNAs (P RNAs), two major loop-helix tertiary contacts (L8-P4 and L18-P8) help to orient the two independently folding S- and C-domains for concerted recognition of precursor tRNA substrates. Here, we analyze the effects of mutations in these tertiary contacts in P RNAs from three different species: (i) the psychrophilic bacterium Pseudoalteromonas translucida (Ptr), (ii) the mesophilic radiation-resistant bacterium Deinococcus radiodurans (Dra), and (iii) the thermophilic bacterium Thermus thermophilus (Tth). We show by UV melting experiments that simultaneous disruption of these two interdomain contacts has a stabilizing effect on all three P RNAs. This can be inferred from reduced RNA unfolding at lower temperatures and a more concerted unfolding at higher temperatures. Thus, when the two domains tightly interact via the tertiary contacts, one domain facilitates structural transitions in the other. P RNA mutants with disrupted interdomain contacts showed severe kinetic defects that were most pronounced upon simultaneous disruption of the L8-P4 and L18-P8 contacts. At 37°C, the mildest effects were observed for the thermostable Tth RNA. A third interdomain contact, L9-P1, makes only a minor contribution to P RNA tertiary folding. Furthermore, D. radiodurans RNase P RNA forms an additional pseudoknot structure between the P9 and P12 of its S-domain. This interaction was found to be particularly crucial for RNase P holoenzyme activity at near-physiological Mg2+ concentrations (2 mM). We further analyzed an exceptionally stable folding trap of the G,C-rich Tth P RNA.

An RNA-centric historical narrative around the Protein Data Bank.

Some of the amazing contributions brought to the scientific community by the Protein Data Bank (PDB) are described. The focus is on nucleic acid structures with a bias toward RNA. The evolution and key roles in science of the PDB and other structural databases for nucleic acids illustrate how small initial ideas can become huge and indispensable resources with the unflinching willingness of scientists to cooperate globally. The progress in the understanding of the molecular interactions driving RNA architectures followed the rapid increase in RNA structures in the PDB. That increase was consecutive to improvements in chemical synthesis and purification of RNA molecules, as well as in biophysical methods for structure determination and computer technology. The RNA modeling efforts from the early beginnings are also described together with their links to the state of structural knowledge and technological development. Structures of RNA and of its assemblies are physical objects, which, together with genomic data, allow us to integrate present-day biological functions and the historical evolution in all living species on earth.

Data, data, burning deep, in the forests of the net.

Continuous and imaginative technological developments are leading to a massive accumulation of various types of data in all areas of biological research. As a result, the central importance of databases is increasing. Databases related to biology must not only be structured using controlled vocabularies, but also be fully integrated into the whole biological domain. To achieve this goal, they must be systematically grounded in biological evolution and exploit the available tools of evolutionary systematics to contribute to our understanding of life processes.

Unusual tertiary pairs in eukaryotic tRNAAla.

tRNA molecules have well-defined sequence conservations that reflect the conserved tertiary pairs maintaining their architecture and functions during the translation processes. An analysis of aligned tRNA sequences present in the GtRNAdb database (the Lowe Laboratory, University of California, Santa Cruz) led to surprising conservations on some cytosolic tRNAs specific for alanine compared to other tRNA species, including tRNAs specific for glycine. First, besides the well-known G3oU70 base pair in the amino acid stem, there is the frequent occurrence of a second wobble pair at G30oU40, a pair generally observed as a Watson-Crick pair throughout phylogeny. Second, the tertiary pair R15/Y48 occurs as a purine-purine R15/A48 pair. Finally, the conserved T54/A58 pair maintaining the fold of the T-loop is observed as a purine-purine A54/A58 pair. The R15/A48 and A54/A58 pairs always occur together. The G30oU40 pair occurs alone or together with these other two pairs. The pairing variations are observed to a variable extent depending on phylogeny. Among eukaryotes, insects display all variations simultaneously, whereas mammals present either the G30oU40 pair or both R15/A48 and A54/A58. tRNAs with the anticodon 34A(I)GC36 are the most prone to display all those pair variations in mammals and insects. tRNAs with anticodon Y34GC36 have preferentially G30oU40 only. These unusual pairs are not observed in bacterial, nor archaeal, tRNAs, probably because of the avoidance of A34-containing anticodons in four-codon boxes. Among eukaryotes, these unusual pairing features were not observed in fungi and nematodes. These unusual structural features may affect, besides aminoacylation, transcription rates (e.g., 54/58) or ribosomal translocation (30/40).

The viral protein NSP1 acts as a ribosome gatekeeper for shutting down host translation and fostering SARS-CoV-2 translation.

SARS-CoV-2 coronavirus is responsible for Covid-19 pandemic. In the early phase of infection, the single-strand positive RNA genome is translated into non-structural proteins (NSP). One of the first proteins produced during viral infection, NSP1, binds to the host ribosome and blocks the mRNA entry channel. This triggers translation inhibition of cellular translation. In spite of the presence of NSP1 on the ribosome, viral translation proceeds however. The molecular mechanism of the so-called viral evasion to NSP1 inhibition remains elusive. Here, we confirm that viral translation is maintained in the presence of NSP1. The evasion to NSP1-inhibition is mediated by the cis-acting RNA hairpin SL1 in the 5'UTR of SARS-CoV-2. NSP1-evasion can be transferred on a reporter transcript by SL1 transplantation. The apical part of SL1 is only required for viral translation. We show that NSP1 remains bound on the ribosome during viral translation. We suggest that the interaction between NSP1 and SL1 frees the mRNA accommodation channel while maintaining NSP1 bound to the ribosome. Thus, NSP1 acts as a ribosome gatekeeper, shutting down host translation or fostering SARS-CoV-2 translation depending on the presence of the SL1 5'UTR hairpin. SL1 is also present and necessary for translation of sub-genomic RNAs in the late phase of the infectious program. Consequently, therapeutic strategies targeting SL1 should affect viral translation at early and late stages of infection. Therefore, SL1 might be seen as a genuine 'Achille heel' of the virus.

Comparative patterns of modified nucleotides in individual tRNA species from a mesophilic and two thermophilic archaea.

To improve and complete our knowledge of archaeal tRNA modification patterns, we have identified and compared the modification pattern (type and location) in tRNAs of three very different archaeal species, Methanococcus maripaludis (a mesophilic methanogen), Pyrococcus furiosus (a hyperthermophile thermococcale), and Sulfolobus acidocaldarius (an acidophilic thermophilic sulfolobale). Most abundant isoacceptor tRNAs (79 in total) for each of the 20 amino acids were isolated by two-dimensional gel electrophoresis followed by in-gel RNase digestions. The resulting oligonucleotide fragments were separated by nanoLC and their nucleotide content analyzed by mass spectrometry (MS/MS). Analysis of total modified nucleosides obtained from complete digestion of bulk tRNAs was also performed. Distinct base- and/or ribose-methylations, cytidine acetylations, and thiolated pyrimidines were identified, some at new positions in tRNAs. Novel, some tentatively identified, modifications were also found. The least diversified modification landscape is observed in the mesophilic Methanococcus maripaludis and the most complex one in Sulfolobus acidocaldarius Notable observations are the frequent occurrence of ac4C nucleotides in thermophilic archaeal tRNAs, the presence of m7G at positions 1 and 10 in Pyrococcus furiosus tRNAs, and the use of wyosine derivatives at position 37 of tRNAs, especially those decoding U1- and C1-starting codons. These results complete those already obtained by others with sets of archaeal tRNAs from Methanocaldococcus jannaschii and Haloferax volcanii.