Recommendations for bacterial ribosome profiling experiments based on bioinformatic evaluation of published data

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 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.

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.

Limited evidence for protein products of non-coding transcripts in the HEK293T cellular cytosol

Ribosome profiling has revealed translation outside of canonical coding sequences (CDSs) including translation of short upstream ORFs, long non-coding RNAs, overlapping ORFs, ORFs in UTRs or ORFs in alternative reading frames. Studies combining mass spectrometry, ribosome profiling and CRISPR-based screens showed that hundreds of ORFs derived from non-coding transcripts produce (micro)proteins, while other studies failed to find evidence for such types of non-canonical translation products. Here, we attempted to discover translation products from non-coding regions by strongly reducing the complexity of the sample prior to mass spectrometric analysis. We used an extended database as the search space and applied stringent filtering of the identified peptides to find evidence for novel translation events. We show that, theoretically our strategy facilitates the detection of translation events of transcripts from non-coding regions, but experimentally only find 19 peptides that might originate from such translation events. Finally, Virotrap based interactome analysis of two N-terminal proteoforms originating from non-coding regions finally showed the functional potential of these novel proteins.

Comparison of rRNA depletion methods for efficient bacterial mRNA sequencing

Current methods of high-throughput RNA sequencing of prokaryotes, including transcriptome analysis or ribosomal profiling, need deep sequencing to achieve sufficient numbers of effective reads (e.g., mapping to mRNA) in order to also find weakly expressed genetic elements. The fraction of high-quality reads mapping to coding RNAs (i.e., mRNA) is mainly influenced by the large content of rRNA and, to a lesser extent, tRNA in total RNA. Thus, depletion of rRNA increases coverage and thus sequencing costs. RiboZero, a depletion kit based on probe hybridisation and rRNA-removal was found to be most efficient in the past, but it was discontinued in 2018. To facilitate comparability with previous experiments and to help choose adequate replacements, we compare three commercially available rRNA depletion kits also based on hybridization and magnetic beads, i.e., riboPOOLs, RiboMinus and MICROBExpress, with the former RiboZero. Additionally, we constructed biotinylated probes for magnetic bead capture and rRNA depletion in this study. Based on E. coli, we found similar efficiencies in rRNA depletion for riboPOOLs and the self-made depletion method; both comparable to the former RiboZero, followed by RiboMinus, succeeded by MICROBExpress. Further, our in-house protocol allows customized species-specific rRNA or even tRNA depletion or depletion of other RNA targets. Both, the self-made biotinylated probes and riboPOOLs, were most successful in reducing the rRNA content and thereby increasing sequencing depth concerning mRNA reads. Additionally, the number of reads matching to weakly expressed genes are increased. In conclusion, the self-made specific biotinylated probes and riboPOOLs are an adequate replacement for the former RiboZero. Both are very efficient in depleting rRNAs, increasing mRNA reads and thus sequencing efficiency.

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.

Translation complex stabilization on messenger RNA and footprint profiling to study the RNA responses and dynamics of protein biosynthesis in the cells

During protein biosynthesis, ribosomes bind to messenger (m)RNA, locate its protein-coding information, and translate the nucleotide triplets sequentially as codons into the corresponding sequence of amino acids, forming proteins. Non-coding mRNA features, such as 5' and 3' untranslated regions (UTRs), start sites or stop codons of different efficiency, stretches of slower or faster code and nascent polypeptide interactions can alter the translation rates transcript-wise. Most of the homeostatic and signal response pathways of the cells converge on individual mRNA control, as well as alter the global translation output. Among the multitude of approaches to study translational control, one of the most powerful is to infer the locations of translational complexes on mRNA based on the mRNA fragments protected by these complexes from endonucleolytic hydrolysis, or footprints. Translation complex profiling by high-throughput sequencing of the footprints allows to quantify the transcript-wise, as well as global, alterations of translation, and uncover the underlying control mechanisms by attributing footprint locations and sizes to different configurations of the translational complexes. The accuracy of all footprint profiling approaches critically depends on the fidelity of footprint generation and many methods have emerged to preserve certain or multiple configurations of the translational complexes, often in challenging biological material. In this review, a systematic summary of approaches to stabilize translational complexes on mRNA for footprinting is presented and major findings are discussed. Future directions of translation footprint profiling are outlined, focusing on the fidelity and accuracy of inference of the native in vivo translation complex distribution on mRNA.

HflX is a GTPase that controls hypoxia-induced replication arrest in slow-growing mycobacteria

GTPase high frequency of lysogenization X (HflX) is highly conserved in prokaryotes and acts as a ribosome-splitting factor as part of the heat shock response in Escherichia coli. Here we report that HflX produced by slow-growing Mycobacterium bovis bacillus Calmette-Guérin (BCG) is a GTPase that plays a critical role in the pathogen's transition to a nonreplicating, drug-tolerant state in response to hypoxia. Indeed, HflX-deficient M. bovis BCG (KO) replicated markedly faster in the microaerophilic phase of a hypoxia model that resulted in premature entry into dormancy. The KO mutant displayed hallmarks of nonreplicating mycobacteria, including phenotypic drug resistance, altered morphology, low intracellular ATP levels, and overexpression of Dormancy (Dos) regulon proteins. Mice nasally infected with HflX KO mutant displayed increased bacterial burden in the lungs, spleen, and lymph nodes during the chronic phase of infection, consistent with the higher replication rate observed in vitro in microaerophilic conditions. Unlike fast growing mycobacteria, M. bovis BCG HlfX was not involved in antibiotic resistance under aerobic growth. Proteomics, pull-down, and ribo-sequencing approaches supported that mycobacterial HflX is a ribosome-binding protein that controls translational activity of the cell. With HflX fully conserved between M. bovis BCG and M. tuberculosis, our work provides further insights into the molecular mechanisms deployed by pathogenic mycobacteria to adapt to their hypoxic microenvironment.

Are Antisense Proteins in Prokaryotes Functional?

Many prokaryotic RNAs are transcribed from loci outside of annotated protein coding genes. Across bacterial species hundreds of short open reading frames antisense to annotated genes show evidence of both transcription and translation, for instance in ribosome profiling data. Determining the functional fraction of these protein products awaits further research, including insights from studies of molecular interactions and detailed evolutionary analysis. There are multiple lines of evidence, however, that many of these newly discovered proteins are of use to the organism. Condition-specific phenotypes have been characterized for a few. These proteins should be added to genome annotations, and the methods for predicting them standardized. Evolutionary analysis of these typically young sequences also may provide important insights into gene evolution. This research should be prioritized for its exciting potential to uncover large numbers of novel proteins with extremely diverse potential practical uses, including applications in synthetic biology and responding to pathogens.