The presence of rRNA sequences in polyadenylated RNA and its potential functions

RNA marker modifications reveal the necessity for rigorous preparation protocols to avoid artifacts in epitranscriptomic analysis

The accurate definition of an epitranscriptome is endangered by artefacts resulting from RNA degradation after cell death, a ubiquitous yet little investigated process. By tracing RNA marker modifications through tissue preparation protocols, we identified a major blind spot from daily lab routine, that has massive impact on modification analysis in small RNAs. In particular, m6,6A and Am as co-varying rRNA marker modifications, appeared in small RNA fractions following rRNA degradation in vitro and in cellulo. Analysing mouse tissue at different time points post mortem, we tracked the progress of intracellular RNA degradation after cell death, and found it reflected in RNA modification patterns. Differences were dramatic between liver, where RNA degradation commenced immediately after death, and brain, yielding essentially undamaged RNA. RNA integrity correlated with low amounts of co-varying rRNA markers. Thus validated RNA preparations featured differentially modified tRNA populations whose information content allowed a distinction even among the related brain tissues cortex, cerebellum and hippocampus. Inversely, advanced cell death correlated with high rRNA marker content, and correspondingly little with the naïve state of living tissue. Therefore, unless RNA and tissue preparations are executed with utmost care, interpretation of modification patterns in tRNA and small RNA are prone to artefacts.

New Insights into the Genome Organization of Yeast Double-Stranded RNA LBC Viruses

The yeasts Torulaspora delbrueckii (Td) and Saccharomyces cerevisiae (Sc) may show a killer phenotype that is encoded in dsRNA M viruses (V-M), which require the helper activity of another dsRNA virus (V-LA or V-LBC) for replication. Recently, two TdV-LBCbarr genomes, which share sequence identity with ScV-LBC counterparts, were characterized by high-throughput sequencing (HTS). They also share some similar characteristics with Sc-LA viruses. This may explain why TdV-LBCbarr has helper capability to maintain M viruses, whereas ScV-LBC does not. We here analyze two stretches with low sequence identity (LIS I and LIS II) that were found in TdV-LBCbarr Gag-Pol proteins when comparing with the homologous regions of ScV-LBC. These stretches may result from successive nucleotide insertions or deletions (indels) that allow compensatory frameshift events required to maintain specific functions of the RNA-polymerase, while modifying other functions such as the ability to bind V-M (+)RNA for packaging. The presence of an additional frameshifting site in LIS I may ensure the synthesis of a certain amount of RNA-polymerase until the new compensatory indel appears. Additional 5'- and 3'-extra sequences were found beyond V-LBC canonical genomes. Most extra sequences showed high identity to some stretches of the canonical genomes and can form stem-loop structures. Further, the 3'-extra sequence of two ScV-LBC genomes contains rRNA stretches. The origin and possible functions of these extra sequences are here discussed.

Genome Organization of a New Double-Stranded RNA LA Helper Virus From Wine Torulaspora delbrueckii Killer Yeast as Compared With Its Saccharomyces Counterparts

Wine killer yeasts such as killer strains of Torulaspora delbrueckii and Saccharomyces cerevisiae contain helper large-size (4.6 kb) dsRNA viruses (V-LA) required for the stable maintenance and replication of killer medium-size dsRNA viruses (V-M) which bear the genes that encode for the killer toxin. The genome of the new V-LA dsRNA from the T. delbrueckii Kbarr1 killer yeast (TdV-LAbarr1) was characterized by high-throughput sequencing (HTS). The canonical genome of TdV-LAbarr1 shares a high sequence identity and similar genome organization with its Saccharomyces counterparts. It contains all the known conserved motifs predicted to be necessary for virus translation, packaging, and replication. Similarly, the Gag-Pol amino-acid sequence of this virus contains all the features required for cap-snatching and RNA polymerase activity, as well as the expected regional variables previously found in other LA viruses. Sequence comparison showed that two main clusters (99.2-100% and 96.3-98.8% identity) include most LA viruses from Saccharomyces, with TdV-LAbarr1 being the most distant from all these viruses (61.5-62.5% identity). Viral co-evolution and cross transmission between different yeast species are discussed based on this sequence comparison. Additional 5' and 3' sequences were found in the TdV-LAbarr1 genome as well as in some newly sequenced V-LA genomes from S. cerevisiae. A stretch involving the 5' extra sequence of TdV-LAbarr1 is identical to a homologous stretch close to the 5' end of the canonical sequence of the same virus (self-identity). Our modeling suggests that these stretches can form single-strand stem loops, whose unpaired nucleotides could anneal to create an intramolecular kissing complex. Similar stem loops are also found in the 3' extra sequence of the same virus as well as in the extra sequences of some LA viruses from S. cerevisiae. A possible origin of these extra sequences as well as their function in obviating ssRNA degradation and allowing RNA transcription and replication are discussed.

Slide-seq: A scalable technology for measuring genome-wide expression at high spatial resolution

Spatial positions of cells in tissues strongly influence function, yet a high-throughput, genome-wide readout of gene expression with cellular resolution is lacking. We developed Slide-seq, a method for transferring RNA from tissue sections onto a surface covered in DNA-barcoded beads with known positions, allowing the locations of the RNA to be inferred by sequencing. Using Slide-seq, we localized cell types identified by single-cell RNA sequencing datasets within the cerebellum and hippocampus, characterized spatial gene expression patterns in the Purkinje layer of mouse cerebellum, and defined the temporal evolution of cell type-specific responses in a mouse model of traumatic brain injury. These studies highlight how Slide-seq provides a scalable method for obtaining spatially resolved gene expression data at resolutions comparable to the sizes of individual cells.

The Ribosome as a Missing Link in Prebiotic Evolution III: Over-Representation of tRNA- and rRNA-Like Sequences and Plieofunctionality of Ribosome-Related Molecules Argues for the Evolution of Primitive Genomes from Ribosomal RNA Modules

We propose that ribosomal RNA (rRNA) formed the basis of the first cellular genomes, and provide evidence from a review of relevant literature and proteonomic tests. We have proposed previously that the ribosome may represent the vestige of the first self-replicating entity in which rRNAs also functioned as genes that were transcribed into functional messenger RNAs (mRNAs) encoding ribosomal proteins. rRNAs also encoded polymerases to replicate itself and a full complement of the transfer RNAs (tRNAs) required to translate its genes. We explore here a further prediction of our “ribosome-first” theory: the ribosomal genome provided the basis for the first cellular genomes. Modern genomes should therefore contain an unexpectedly large percentage of tRNA- and rRNA-like modules derived from both sense and antisense reading frames, and these should encode non-ribosomal proteins, as well as ribosomal ones with key cell functions. Ribosomal proteins should also have been co-opted by cellular evolution to play extra-ribosomal functions. We review existing literature supporting these predictions. We provide additional, new data demonstrating that rRNA-like sequences occur at significantly higher frequencies than predicted on the basis of mRNA duplications or randomized RNA sequences. These data support our “ribosome-first” theory of cellular evolution.

Changes in transcript abundance relating to colony collapse disorder in honey bees (Apis mellifera)

Colony collapse disorder (CCD) is a mysterious disappearance of honey bees that has beset beekeepers in the United States since late 2006. Pathogens and other environmental stresses, including pesticides, have been linked to CCD, but a causal relationship has not yet been demonstrated. Because the gut acts as a primary interface between the honey bee and its environment as a site of entry for pathogens and toxins, we used whole-genome microarrays to compare gene expression between guts of bees from CCD colonies originating on both the east and west coasts of the United States and guts of bees from healthy colonies sampled before the emergence of CCD. Considerable variation in gene expression was associated with the geographical origin of bees, but a consensus list of 65 transcripts was identified as potential markers for CCD status. Overall, elevated expression of pesticide response genes was not observed. Genes involved in immune response showed no clear trend in expression pattern despite the increased prevalence of viruses and other pathogens in CCD colonies. Microarray analysis revealed unusual ribosomal RNA fragments that were conspicuously more abundant in the guts of CCD bees. The presence of these fragments may be a possible consequence of picorna-like viral infection, including deformed wing virus and Israeli acute paralysis virus, and may be related to arrested translation. Ribosomal fragment abundance and presence of multiple viruses may prove to be useful diagnostic markers for colonies afflicted with CCD.