Identifying Small Open Reading Frames in Prokaryotes with Ribosome Profiling

A survey of experimental and computational identification of small proteins

Small proteins (SPs) are typically characterized as eukaryotic proteins shorter than 100 amino acids and prokaryotic proteins shorter than 50 amino acids. Historically, they were disregarded because of the arbitrary size thresholds to define proteins. However, recent research has revealed the existence of many SPs and their crucial roles. Despite this, the identification of SPs and the elucidation of their functions are still in their infancy. To pave the way for future SP studies, we briefly introduce the limitations and advancements in experimental techniques for SP identification. We then provide an overview of available computational tools for SP identification, their constraints, and their evaluation. Additionally, we highlight existing resources for SP research. This survey aims to initiate further exploration into SPs and encourage the development of more sophisticated computational tools for SP identification in prokaryotes and microbiomes.

Ribosome profiling reveals downregulation of UMP biosynthesis as the major early response to phage infection

Bacteria have evolved diverse defense mechanisms to counter bacteriophage attacks. Genetic programs activated upon infection characterize phage-host molecular interactions and ultimately determine the outcome of the infection. In this study, we applied ribosome profiling to monitor protein synthesis during the early stages of sk1 bacteriophage infection in Lactococcus cremoris. Our analysis revealed major changes in gene expression within 5 minutes of sk1 infection. Notably, we observed a specific and severe downregulation of several pyr operons which encode enzymes required for uridine monophosphate biosynthesis. Consistent with previous findings, this is likely an attempt of the host to starve the phage of nucleotides it requires for propagation. We also observed a gene expression response that we expect to benefit the phage. This included the upregulation of 40 ribosome proteins that likely increased the host's translational capacity, concurrent with a downregulation of genes that promote translational fidelity (lepA and raiA). In addition to the characterization of host-phage gene expression responses, the obtained ribosome profiling data enabled us to identify two putative recoding events as well as dozens of loci currently annotated as pseudogenes that are actively translated. Furthermore, our study elucidated alterations in the dynamics of the translation process, as indicated by time-dependent changes in the metagene profile, suggesting global shifts in translation rates upon infection. Additionally, we observed consistent modifications in the ribosome profiles of individual genes, which were apparent as early as 2 minutes post-infection. The study emphasizes our ability to capture rapid alterations of gene expression during phage infection through ribosome profiling.

IMPORTANCE

The ribosome profiling technology has provided invaluable insights for understanding cellular translation and eukaryotic viral infections. However, its potential for investigating host-phage interactions remains largely untapped. Here, we applied ribosome profiling to Lactococcus cremoris cultures infected with sk1, a major infectious agent in dairy fermentation processes. This revealed a profound downregulation of genes involved in pyrimidine nucleotide synthesis at an early stage of phage infection, suggesting an anti-phage program aimed at restricting nucleotide availability and, consequently, phage propagation. This is consistent with recent findings and contributes to our growing appreciation for the role of nucleotide limitation as an anti-viral strategy. In addition to capturing rapid alterations in gene expression levels, we identified translation occurring outside annotated regions, as well as signatures of non-standard translation mechanisms. The gene profiles revealed specific changes in ribosomal densities upon infection, reflecting alterations in the dynamics of the translation process.

Genome-Scale Analysis of the Structure and Function of RNA Pathways and Networks in Pseudomonas aeruginosa

In recent years, several genome-wide approaches based on RNA sequencing (RNA-seq) have been developed. These methods allow a comprehensive and dynamic view of the structure and function of the multi-layered RNA pathways and networks. Many of these approaches, including the promising one of single-cell transcriptome analysis, have been successfully applied to Pseudomonas aeruginosa. However, we are only at the beginning because only a few surrounding conditions have been considered. Here, we aim to illustrate the different types of approaches based on RNA-seq that will lead us in the future to a better understanding of the dynamics of RNA biology in P. aeruginosa.

Ribosome Profiling Methods Adapted to the Study of RNA-Dependent Translation Regulation in Staphylococcus aureus

Noncoding RNAs, including regulatory RNAs (sRNAs), are instrumental in regulating gene expression in pathogenic bacteria, allowing them to adapt to various stresses encountered in their host environments. Staphylococcus aureus is a well-studied model for RNA-mediated regulation of virulence and pathogenicity, with sRNAs playing significant roles in shaping S. aureus interactions with human and animal hosts. By modulating the translation and/or stability of target mRNAs, sRNAs regulate the synthesis of virulence factors and regulatory proteins required for pathogenesis. Moreover, perturbation of the levels of RNA modifications in two other classes of noncoding RNAs, rRNAs, and tRNAs, has been proposed to contribute to stress adaptation. However, the study of how these various factors affect translation regulation has often been restricted to specific genes, using in vivo reporters and/or in vitro translation systems. Genome-wide sequencing approaches offer novel perspectives for studying RNA-dependent regulation. In particular, ribosome profiling methods provide a powerful resource for characterizing the overall landscape of translational regulation, contributing to a better understanding of S. aureus physiopathology. Here, we describe protocols that we have adapted to perform ribosome profiling in S. aureus.

Ribosome profiling reveals the fine-tuned response of Escherichia coli to mild and severe acid stress

IMPORTANCE

Bacteria react very differently to survive in acidic environments, such as the human gastrointestinal tract. Escherichia coli is one of the extremely acid-resistant bacteria and has a variety of acid-defense mechanisms. Here, we provide the first genome-wide overview of the adaptations of E. coli K-12 to mild and severe acid stress at both the transcriptional and translational levels. Using ribosome profiling and RNA sequencing, we uncover novel adaptations to different degrees of acidity, including previously hidden stress-induced small proteins and novel key transcription factors for acid defense, and report mRNAs with pH-dependent differential translation efficiency. In addition, we distinguish between acid-specific adaptations and general stress response mechanisms using denoising autoencoders. This workflow represents a powerful approach that takes advantage of next-generation sequencing techniques and machine learning to systematically analyze bacterial stress responses.

Small protein mediates inhibition of ammonium transport in Methanosarcina mazei-an ancient mechanism?

IMPORTANCE

Small proteins containing fewer than 70 amino acids, which were previously disregarded due to computational prediction and biochemical detection challenges, have gained increased attention in the scientific community in recent years. However, the number of functionally characterized small proteins, especially in archaea, is still limited. Here, by using biochemical and genetic approaches, we demonstrate a crucial role of the small protein sP36 in the nitrogen metabolism of M. mazei, which modulates the ammonium transporter AmtB1 according to nitrogen availability. This modulation might represent an ancient archaeal mechanism of AmtB1 inhibition, in contrast to the well-studied uridylylation-dependent regulation in bacteria.

Recent insights into the world of dual-function bacterial sRNAs

Dual-function sRNAs refer to a small subgroup of small regulatory RNAs that merges base-pairing properties of antisense RNAs with peptide-encoding properties of mRNA. Both functions can be part of either same or in another metabolic pathway. Here, we want to update the knowledge of to the already known dual-function sRNAs and review the six new sRNAs found since 2017 regarding their structure, functional mechanisms, evolutionary conservation, and role in the regulation of distinct biological/physiological processes. The increasing identification of dual-function sRNAs through bioinformatics approaches, RNomics and RNA-sequencing and the associated increase in regulatory understanding will likely continue to increase at the same rate in the future. This may improve our understanding of the physiology, virulence and resistance of bacteria, as well as enable their use in technical applications. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

The power of the small: the underestimated role of small proteins in bacterial and archaeal physiology

Small proteins encoded by small open-reading frames (sORFs) (≤70 aa) were overlooked for decades due to methodological reasons and are thus often missing in genome annotations. Novel detection methods such as ribosome profiling (Ribo-Seq) and mass spectrometry optimized for small proteins (peptidomics) have opened up a new field of interest and several catalogs of small proteins in bacteria and archaea have been recently reported. Many translated sORFs have been discovered in genomic locations previously thought to be noncoding, such as 5' or 3' untranslated regions or well-studied regulatory small RNAs (sRNAs). Even within longer ORFs, additional functional sORFs have been detected. Today, only a small proportion is characterized, but those small proteins indicate important and diverse functions in cellular physiology. Here, we summarize recently characterized small proteins involved in microbial metabolism.

Two for one: regulatory RNAs that encode small proteins

RNAs are commonly categorized as being either protein-coding mRNAs or noncoding RNAs. However, an increasing number of transcripts, in organisms ranging from bacteria to humans, are being found to have both coding and noncoding functions. In some cases, the sequences encoding the protein and the regulatory RNA functions are separated, while in other cases the sequences overlap. The protein and RNA can regulate similar or distinct pathways. Here we describe examples illustrating how these dual-function (also denoted bifunctional or dual-component) RNAs are identified and their mechanisms of action and cellular roles. We also discuss the synergy or competition between coding and RNA activity and how these regulators evolved, as well as how more dual-function RNAs might be discovered and exploited.

StORF-Reporter: finding genes between genes

Large regions of prokaryotic genomes are currently without any annotation, in part due to well-established limitations of annotation tools. For example, it is routine for genes using alternative start codons to be misreported or completely omitted. Therefore, we present StORF-Reporter, a tool that takes an annotated genome and returns regions that may contain missing CDS genes from unannotated regions. StORF-Reporter consists of two parts. The first begins with the extraction of unannotated regions from an annotated genome. Next, Stop-ORFs (StORFs) are identified in these unannotated regions. StORFs are open reading frames that are delimited by stop codons and thus can capture those genes most often missing in genome annotations. We show this methodology recovers genes missing from canonical genome annotations. We inspect the results of the genomes of model organisms, the pangenome of Escherichia coli, and a set of 5109 prokaryotic genomes of 247 genera from the Ensembl Bacteria database. StORF-Reporter extended the core, soft-core and accessory gene collections, identified novel gene families and extended families into additional genera. The high levels of sequence conservation observed between genera suggest that many of these StORFs are likely to be functional genes that should now be considered for inclusion in canonical annotations.

Exposing the small protein load of bacterial life

The ever-growing repertoire of genomic techniques continues to expand our understanding of the true diversity and richness of prokaryotic genomes. Riboproteogenomics laid the foundation for dynamic studies of previously overlooked genomic elements. Most strikingly, bacterial genomes were revealed to harbor robust repertoires of small open reading frames (sORFs) encoding a diverse and broadly expressed range of small proteins, or sORF-encoded polypeptides (SEPs). In recent years, continuous efforts led to great improvements in the annotation and characterization of such proteins, yet many challenges remain to fully comprehend the pervasive nature of small proteins and their impact on bacterial biology. In this work, we review the recent developments in the dynamic field of bacterial genome reannotation, catalog the important biological roles carried out by small proteins and identify challenges obstructing the way to full understanding of these elusive proteins.

Small proteins in Gram-positive bacteria

Small proteins comprising less than 100 amino acids have been often ignored in bacterial genome annotations. About 10 years ago, focused efforts started to investigate whole peptidomes, which resulted in the discovery of a multitude of small proteins, but only a number of them have been characterized in detail. Generally, small proteins can be either membrane or cytosolic proteins. The latter interact with larger proteins, RNA or even metal ions. Here, we summarize our current knowledge on small proteins from Gram-positive bacteria with a special emphasis on the model organism Bacillus subtilis. Our examples include membrane-bound toxins of type I toxin-antitoxin systems, proteins that block the assembly of higher order structures, regulate sporulation or modulate the RNA degradosome. We do not consider antimicrobial peptides. Furthermore, we present methods for the identification and investigation of small proteins.

Bacterial acid stress response: from cellular changes to antibiotic tolerance and phenotypic heterogeneity

Most bacteria are neutralophiles but can survive fluctuations in pH in their environment. Herein, we provide an overview of the adaptation of several human, soil, and food bacteria to acid stress, mainly based on next-generation sequencing studies, highlighting common and specific strategies. We also discuss the interplay between acid stress response and antibiotic tolerance, as well as the response of individual cells.

Little reason to call them small noncoding RNAs

Hundreds of different species of small RNAs can populate a bacterial cell. This small transcriptome contains important information for the adaptation of cellular physiology to environmental changes. Underlying cellular networks involving small RNAs are RNA-RNA and RNA-protein interactions, which are often intertwined. In addition, small RNAs can function as mRNAs. In general, small RNAs are referred to as noncoding because very few are known to contain translated open reading frames. In this article, we intend to highlight that the number of small RNAs that fall within the set of translated RNAs is bound to increase. In addition, we aim to emphasize that the dynamics of the small transcriptome involve different functional codes, not just the genetic code. Therefore, since the role of small RNAs is always code-driven, we believe that there is little reason to continue calling them small noncoding RNAs.

Initiation at AUGUG and GUGUG sequences can lead to translation of overlapping reading frames in E. coli

During initiation, the ribosome is tasked to efficiently recognize open reading frames (ORFs) for accurate and fast translation of mRNAs. A critical step is start codon recognition, which is modulated by initiation factors, mRNA structure, a Shine Dalgarno (SD) sequence and the start codon itself. Within the Escherichia coli genome, we identified more than 50 annotated initiation sites harboring AUGUG or GUGUG sequence motifs that provide two canonical start codons, AUG and GUG, in immediate proximity. As these sites may challenge start codon recognition, we studied if and how the ribosome is accurately guided to the designated ORF, with a special focus on the SD sequence as well as adenine at the fourth coding sequence position (A4). By in vitro and in vivo experiments, we characterized key requirements for unambiguous start codon recognition, but also discovered initiation sites that lead to the translation of both overlapping reading frames. Our findings corroborate the existence of an ambiguous translation initiation mechanism, implicating a multitude of so far unrecognized ORFs and translation products in bacteria.

Flow-Seq Method: Features and Application in Bacterial Translation Studies

The Flow-seq method is based on using reporter construct libraries, where a certain element regulating the gene expression of fluorescent reporter proteins is represented in many thousands of variants. Reporter construct libraries are introduced into cells, sorted according to their fluorescence level, and then subjected to next-generation sequencing. Therefore, it turns out to be possible to identify patterns that determine the expression efficiency, based on tens and hundreds of thousands of reporter constructs in one experiment. This method has become common in evaluating the efficiency of protein synthesis simultaneously by multiple mRNA variants. However, its potential is not confined to this area. In the presented review, a comparative analysis of the Flow-seq method and other alternative approaches used for translation efficiency evaluation of mRNA was carried out; the features of its application and the results obtained by Flow-seq were also considered.

A Practical Guide to Small Protein Discovery and Characterization Using Mass Spectrometry

Small proteins of up to ∼50 amino acids are an abundant class of biomolecules across all domains of life. Yet due to the challenges inherent in their size, they are often missed in genome annotations, and are difficult to identify and characterize using standard experimental approaches. Consequently, we still know few small proteins even in well-studied prokaryotic model organisms. Mass spectrometry (MS) has great potential for the discovery, validation, and functional characterization of small proteins. However, standard MS approaches are poorly suited to the identification of both known and novel small proteins due to limitations at each step of a typical proteomics workflow, i.e., sample preparation, protease digestion, liquid chromatography, MS data acquisition, and data analysis. Here, we outline the major MS-based workflows and bioinformatic pipelines used for small protein discovery and validation. Special emphasis is placed on highlighting the adjustments required to improve detection and data quality for small proteins. We discuss both the unbiased detection of small proteins and the targeted analysis of small proteins of interest. Finally, we provide guidelines to prioritize novel small proteins, and an outlook on methods with particular potential to further improve comprehensive discovery and characterization of small proteins.

Small Proteins; Big Questions

In recent years, there has been increased appreciation that a whole category of proteins, small proteins of around 50 amino acids or fewer in length, has been missed by annotation as well as by genetic and biochemical assays. With the increased recognition that small proteins are stable within cells and have regulatory functions, there has been intensified study of these proteins. As a result, important questions about small proteins in bacteria and archaea are coming to the fore. Here, we give an overview of these questions, the initial answers, and the approaches needed to address these questions more fully. More detailed discussions of how small proteins can be identified by ribosome profiling and mass spectrometry approaches are provided by two accompanying reviews (N. Vazquez-Laslop, C. M. Sharma, A. S. Mankin, and A. R. Buskirk, J Bacteriol 204:e00294-21, 2022, https://doi.org/10.1128/JB.00294-21; C. H. Ahrens, J. T. Wade, M. M. Champion, and J. D. Langer, J Bacteriol 204:e00353-21, 2022, https://doi.org/10.1128/JB.00353-21). We are excited by the prospects of new insights and possible therapeutic approaches coming from this emerging field.

RiboReport - benchmarking tools for ribosome profiling-based identification of open reading frames in bacteria

Small proteins encoded by short open reading frames (ORFs) with 50 codons or fewer are emerging as an important class of cellular macromolecules in diverse organisms. However, they often evade detection by proteomics or in silico methods. Ribosome profiling (Ribo-seq) has revealed widespread translation in genomic regions previously thought to be non-coding, driving the development of ORF detection tools using Ribo-seq data. However, only a handful of tools have been designed for bacteria, and these have not yet been systematically compared. Here, we aimed to identify tools that use Ribo-seq data to correctly determine the translational status of annotated bacterial ORFs and also discover novel translated regions with high sensitivity. To this end, we generated a large set of annotated ORFs from four diverse bacterial organisms, manually labeled for their translation status based on Ribo-seq data, which are available for future benchmarking studies. This set was used to investigate the predictive performance of seven Ribo-seq-based ORF detection tools (REPARATION_blast, DeepRibo, Ribo-TISH, PRICE, smORFer, ribotricer and SPECtre), as well as IRSOM, which uses coding potential and RNA-seq coverage only. DeepRibo and REPARATION_blast robustly predicted translated ORFs, including sORFs, with no significant difference for ORFs in close proximity to other genes versus stand-alone genes. However, no tool predicted a set of novel, experimentally verified sORFs with high sensitivity. Start codon predictions with smORFer show the value of initiation site profiling data to further improve the sensitivity of ORF prediction tools in bacteria. Overall, we find that bacterial tools perform well for sORF detection, although there is potential for improving their performance, applicability, usability and reproducibility.