RiboPlotR: a visualization tool for periodic Ribo-seq reads

Background Ribo-seq has revolutionized the study of genome-wide mRNA translation. High-quality Ribo-seq data display strong 3-nucleotide (nt) periodicity, which corresponds to translating ribosomes deciphering three nts at a time. While 3-nt periodicity has been widely used to study novel translation events such as upstream ORFs in 5′ untranslated regions and small ORFs in presumed non-coding RNAs, tools that allow the visualization of these events remain underdeveloped. Results RiboPlotR is a visualization package written in R that presents both RNA-seq coverage and Ribo-seq reads in genomic coordinates for all annotated transcript isoforms of a gene. Specifically, for individual isoform models, RiboPlotR plots Ribo-seq data in the context of gene structures, including 5′ and 3′ untranslated regions and introns, and it presents the reads for all three reading frames in three different colors. The inclusion of gene structures and color-coding the reading frames facilitate observing new translation events and identifying potential regulatory mechanisms. Conclusions RiboPlotR is freely available (https://github.com/hsinyenwu/RiboPlotR and https://sourceforge.net/projects/riboplotr/) and allows the visualization of translated features identified in Ribo-seq data.

RiboPlotR: a Visualization Tool for Periodic Ribo-seq Reads

Background: Ribo-seq has revolutionized the study of genome-wide mRNA translation. High-quality Ribo-seq data display strong 3-nucleotide (nt) periodicity, which corresponds to translating ribosomes deciphering three nts at a time. While 3-nt periodicity has been widely used to study novel translation events such as upstream ORFs in 5’ untranslated regions and small ORFs in presumed non-coding RNAs, tools that allow the visualization of these events remain underdeveloped.Results: RiboPlotR is a visualization package written in R that presents both RNA-seq coverage and Ribo-seq reads in genomic coordinates for all annotated transcript isoforms of a gene. Specifically, for individual isoform models, RiboPlotR plots Ribo-seq data related to splice junctions and presents the reads for all three reading frames in three different colors. Moreover, RiboPlotR shows Ribo-seq reads in upstream ORFs, 5' and 3' untranslated regions and introns, which is critical for observing new translation events and identifying potential regulatory mechanisms.Conclusions: RiboPlotR is freely available (https://github.com/hsinyenwu/RiboPlotR and https://sourceforge.net/projects/riboplotr/) and allows the visualization of translated features identified in Ribo-seq data.

Mycobacterium tuberculosis VapC4 toxin engages small ORFs to initiate an integrated oxidative and copper stress response

The Mycobacterium tuberculosis (Mtb) VapBC4 toxin–antitoxin system is essential for the establishment of Mtb infection. Using a multitier, systems-level approach, we uncovered the sequential molecular events triggered by the VapC4 toxin that activate a circumscribed set of critical stress survival pathways which undoubtedly underlie Mtb virulence. VapC4 exclusively inactivated the sole transfer RNACys (tRNACys) through cleavage at a single site within the anticodon sequence. Depletion of the pool of tRNACys led to ribosome stalling at Cys codons within actively translating messenger RNAs. Genome mapping of these Cys-stalled ribosomes unexpectedly uncovered several unannotated Cys-containing open reading frames (ORFs). Four of these are small ORFs (sORFs) encoding Cys-rich proteins of fewer than 50 amino acids that function as Cys-responsive attenuators that engage ribosome stalling at tracts of Cys codons to control translation of downstream genes. Thus, VapC4 mimics a state of Cys starvation, which then activates Cys attenuation at sORFs to globally redirect metabolism toward the synthesis of free Cys. The resulting newly enriched pool of Cys feeds into the synthesis of mycothiol, the glutathione counterpart in this pathogen that is responsible for maintaining cellular redox homeostasis during oxidative stress, as well as into a circumscribed subset of cellular pathways that enable cells to defend against oxidative and copper stresses characteristically endured by Mtb within macrophages. Our ability to pinpoint activation or down-regulation of pathways that collectively align with Mtb virulence–associated stress responses and the nonreplicating persistent state brings to light a direct and vital role for the VapC4 toxin in mediating these critical pathways.

Identifying small ORFs in prokaryotes with ribosome profiling

Small proteins encoded by ORFs shorter than 50 codons (sORFs) are often overlooked by annotation engines and are difficult to characterize using traditional biochemical techniques. Ribosome profiling has tremendous potential to empirically improve the annotations of prokaryotic genomes. Recent improvements in ribosome profiling methods for bacterial model organisms have revealed many new sORFs in well-characterized genomes. Antibiotics that trap ribosomes just after initiation have played a key role in these developments by allowing unambiguous identification of the start codons (and hence the reading frame) for novel ORFs. Here we describe these new methods and highlight critical controls and considerations for adapting ribosome profiling to different prokaryotic species.

RiboNT: A Noise-Tolerant Predictor of Open Reading Frames from Ribosome-Protected Footprints

Ribo-seq, also known as ribosome profiling, refers to the sequencing of ribosome-protected mRNA fragments (RPFs). This technique has greatly advanced our understanding of translation and facilitated the identification of novel open reading frames (ORFs) within untranslated regions or non-coding sequences as well as the identification of non-canonical start codons. However, the widespread application of Ribo-seq has been hindered because obtaining periodic RPFs requires a highly optimized protocol, which may be difficult to achieve, particularly in non-model organisms. Furthermore, the periodic RPFs are too short (28 nt) for accurate mapping to polyploid genomes, but longer RPFs are usually produced with a compromise in periodicity. Here we present RiboNT, a noise-tolerant ORF predictor that can utilize RPFs with poor periodicity. It evaluates RPF periodicity and automatically weighs the support from RPFs and codon usage before combining their contributions to identify translated ORFs. The results demonstrate the utility of RiboNT for identifying both long and small ORFs using RPFs with either good or poor periodicity. We implemented the pipeline on a dataset of RPFs with poor periodicity derived from membrane-bound polysomes of Arabidopsis thaliana seedlings and identified several small ORFs (sORFs) evolutionarily conserved in diverse plant species. RiboNT should greatly broaden the application of Ribo-seq by minimizing the requirement of RPF quality and allowing the use of longer RPFs, which is critical for organisms with complex genomes because these RPFs can be more accurately mapped to the position from which they were derived.

Geminiviruses encode additional small proteins with specific subcellular localizations and virulence function

Geminiviruses are plant viruses with limited coding capacity. Geminivirus-encoded proteins are traditionally identified by applying a 10-kDa arbitrary threshold; however, it is increasingly clear that small proteins play relevant roles in biological systems, which calls for the reconsideration of this criterion. Here, we show that geminiviral genomes contain additional ORFs. Using tomato yellow leaf curl virus, we demonstrate that some of these small ORFs are expressed during the infection, and that the encoded proteins display specific subcellular localizations. We prove that the largest of these additional ORFs, which we name V3, is required for full viral infection, and that the V3 protein localizes in the Golgi apparatus and functions as an RNA silencing suppressor. These results imply that the repertoire of geminiviral proteins can be expanded, and that getting a comprehensive overview of the molecular plant-geminivirus interactions will require the detailed study of small ORFs so far neglected.

Define protein variant functions with high-complexity mutagenesis libraries and enhanced mutation detection software

Open reading frame (ORF) variant libraries have advanced our ability to query the functions of a large number of variants of a protein simultaneously. A variant library targeting a full-length ORF typically consists of all possible single-amino-acid substitutions and a stop codon at each amino-acid position. Because a single codon variation separates a variant from the template ORF, variant quantification presents the most profound challenge. Efforts such as dividing a library into sub-libraries for direct sequencing or using a tag-directed subassembly approach are practical only with small ORFs. Our approach, on the other hand, features single-pool libraries for all genes up to 3600bp (EGFR), and an enhanced variant-detecting toolkit. Having succeeded in processing screens of ~20 ORF variant libraries, this tool calls variants reliably, and also presents variant annotations in datafiles to enable analyses that have in turn reshaped our strategies governing library design, screen deconvolution, sequencing and its analysis.

Identification of a micropeptide and multiple secondary cell genes that modulate Drosophila male reproductive success

Even in well-characterized genomes, many transcripts are considered noncoding RNAs (ncRNAs) simply due to the absence of large open reading frames (ORFs). However, it is now becoming clear that many small ORFs (smORFs) produce peptides with important biological functions. In the process of characterizing the ribosome-bound transcriptome of an important cell type of the seminal fluid-producing accessory gland of Drosophila melanogaster, we detected an RNA, previously thought to be noncoding, called male-specific abdominal (msa). Notably, msa is nested in the HOX gene cluster of the Bithorax complex and is known to contain a micro-RNA within one of its introns. We find that this RNA encodes a “micropeptide” (9 or 20 amino acids, MSAmiP) that is expressed exclusively in the secondary cells of the male accessory gland, where it seems to accumulate in nuclei. Importantly, loss of function of this micropeptide causes defects in sperm competition. In addition to bringing insights into the biology of a rare cell type, this work underlines the importance of small peptides, a class of molecules that is now emerging as important actors in complex biological processes.

Absence of 4-Formylaminooxyvinylglycine Production by Pseudomonas fluorescens WH6 Results in Resource Reallocation from Secondary Metabolite Production to Rhizocompetence

Pseudomonas fluorescens WH6 produces the non-proteinogenic amino acid 4-formylaminooxyvinylglycine (FVG), a secondary metabolite with antibacterial and pre-emergent herbicidal activities. The gvg operon necessary for FVG production encodes eight required genes: one regulatory (gvgR), two of unknown functional potential (gvgA and C), three with putative biosynthetic function (gvgF, H, and I), and two small ORFs (gvgB and G). To gain insight into the role of GvgA and C in FVG production, we compared the transcriptome of knockout (KO) mutants of gvgR, A, and C to wild type (WT) to test two hypotheses: (1) GvgA and GvgC play a regulatory role in FVG production and (2) non-gvg cluster genes are regulated by GvgA and GvgC. Our analyses show that, collectively, 687 genes, including the gvg operon, are differentially expressed in all KO strains versus WT, representing >10% of the genome. Fifty-one percent of these genes were similarly regulated in all KO strains with GvgC having the greatest number of uniquely regulated genes. Additional transcriptome data suggest cluster regulation through feedback of a cluster product. We also discovered that FVG biosynthesis is regulated by L-glu, L-asp, L-gln, and L-asn and that resources are reallocated in KO strains to increase phenotypes involved in rhizocompetence including motility, biofilm formation, and denitrification. Altogether, differential transcriptome analyses of mutants suggest that regulation of the cluster is multifaceted and the absence of FVG production or its downregulation can dramatically shift the lifestyle of WH6.

Geminiviral genomes encode additional proteins with specific subcellular localizations and virulence function

ABSTRACTGeminiviruses are plant viruses with limited coding capacity. Geminivirus-encoded proteins were identified applying a 10-kDa arbitrary threshold; however, it is increasingly clear that small proteins play relevant roles in biological systems, which calls for the reconsideration of this criterion. Here, we show that geminiviral genomes contain additional ORFs. Using tomato yellow leaf curl virus, we demonstrate that some of these novel ORFs are expressed during the infection, and that the encoded proteins display specific subcellular localizations. We prove that the largest of these new ORFs, which we name V3, is required for full viral infection, and that the V3 protein localizes in the Golgi apparatus and functions as an RNA silencing suppressor. These results imply that the repertoire of geminiviral proteins can be expanded, and that getting a comprehensive overview of the molecular plant-geminivirus interactions will require the detailed study of small ORFs so far neglected.