STATR: A simple analysis pipeline of Ribo-Seq in bacteria

Gene expression changes in response to diverse environmental stimuli to regulate numerous cellular functions. Genes are expressed into their functional products with the help of messenger RNA (mRNA). Thus, measuring levels of mRNA in cells is important to understand cellular functions. With advances in next-generation sequencing (NGS), the abundance of cellular mRNA has been elucidated via transcriptome sequencing. However, several studies have found a discrepancy between mRNA abundance and protein levels induced by translational regulation, including different rates of ribosome entry and translational pausing. As such, the levels of mRNA are not necessarily a direct representation of the protein levels found in a cell. To determine a more precise way to measure protein expression in cells, the analysis of the levels of mRNA associated with ribosomes is being adopted. With an aid of NGS techniques, a single nucleotide resolution footprint of the ribosome was determined using a method known as Ribo-Seq or ribosome profiling. This method allows for the high-throughput measurement of translation in vivo , which was further analyzed to determine the protein synthesis rate, translational pausing, and cellular responses toward a variety of environmental changes. Here, we describe a simple analysis pipeline for Ribo-Seq in bacteria, so-called simple translatome analysis tool for Ribo-Seq (STATR). STATR can be used to carry out the primary processing of Ribo-Seq data, subsequently allowing for multiple levels of translatome study, from experimental validation to in-depth analyses. A command-by-command explanation is provided here to allow a broad spectrum of biologists to easily reproduce the analysis.

[1]  D. Sturdevant,et al.  RNAIII-independent target gene control by the agr quorum-sensing system: insight into the evolution of virulence regulation in Staphylococcus aureus. , 2008, Molecular cell.

[2]  Gene-Wei Li,et al.  The anti-Shine-Dalgarno sequence drives translational pausing and codon choice in bacteria , 2012, Nature.

[3]  Kristin Reiche,et al.  The primary transcriptome of the major human pathogen Helicobacter pylori , 2010, Nature.

[4]  Byung-Kwan Cho,et al.  Adaptive laboratory evolution of a genome-reduced Escherichia coli , 2019, Nature Communications.

[5]  Broňa Brejová,et al.  Massive programmed translational jumping in mitochondria , 2014, Proceedings of the National Academy of Sciences.

[6]  Rachel Green,et al.  High-precision analysis of translational pausing by ribosome profiling in bacteria lacking EFP. , 2015, Cell reports.

[7]  Christopher A. Vakulskas,et al.  Global role of the bacterial post-transcriptional regulator CsrA revealed by integrated transcriptomics , 2017, Nature Communications.

[8]  R. Green,et al.  A systematically-revised ribosome profiling method for bacteria reveals pauses at single-codon resolution , 2019, eLife.

[9]  Eric Rivals,et al.  Ribo-seq enlightens codon usage bias , 2017, DNA research : an international journal for rapid publication of reports on genes and genomes.

[10]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[11]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[12]  Jos Vanderleyden,et al.  RNA-binding proteins involved in post-transcriptional regulation in bacteria , 2015, Front. Microbiol..

[13]  G. Weinstock,et al.  Subtle genetic changes enhance virulence of methicillin resistant and sensitive Staphylococcus aureus , 2007, BMC Microbiology.

[14]  M. Gerstein,et al.  The Transcriptional Landscape of the Yeast Genome Defined by RNA Sequencing , 2008, Science.

[15]  Nicholas T. Ingolia,et al.  Genome-Wide Analysis in Vivo of Translation with Nucleotide Resolution Using Ribosome Profiling , 2009, Science.

[16]  V. Gladyshev,et al.  Genome-wide ribosome profiling reveals complex translational regulation in response to oxidative stress , 2012, Proceedings of the National Academy of Sciences.

[17]  Min Woo Kim,et al.  The dynamic transcriptional and translational landscape of the model antibiotic producer Streptomyces coelicolor A3(2) , 2016, Nature Communications.

[18]  B. Palsson,et al.  Genome-scale analysis of Methicillin-resistant Staphylococcus aureus USA300 reveals a tradeoff between pathogenesis and drug resistance , 2018, Scientific Reports.

[19]  C. Burge,et al.  Widespread regulation of translation by elongation pausing in heat shock , 2013, Molecular cell.

[20]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[21]  B. Palsson,et al.  A streamlined ribosome profiling protocol for the characterization of microorganisms. , 2015, BioTechniques.

[22]  L. Steinmetz,et al.  Extensive transcriptional heterogeneity revealed by isoform profiling , 2013, Nature.

[23]  Aaron R. Quinlan,et al.  BIOINFORMATICS APPLICATIONS NOTE , 2022 .

[24]  Pascale Cossart,et al.  Term-seq reveals abundant ribo-regulation of antibiotics resistance in bacteria , 2016, Science.

[25]  J. Puglisi,et al.  Dynamic pathways of -1 translational frameshifting , 2014, Nature.

[26]  M. Gilmour,et al.  Tedizolid: A Novel Oxazolidinone with Potent Activity Against Multidrug-Resistant Gram-Positive Pathogens , 2015, Drugs.

[27]  A. Renslo Antibacterial oxazolidinones: emerging structure–toxicity relationships , 2010, Expert review of anti-infective therapy.

[28]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[29]  Colin N. Dewey,et al.  RNA-Seq gene expression estimation with read mapping uncertainty , 2009, Bioinform..

[30]  Karl W Barber,et al.  Heterogeneity of Stop Codon Readthrough in Single Bacterial Cells and Implications for Population Fitness. , 2017, Molecular cell.

[31]  B. Larkins,et al.  Polyribosomes from Peas: II. Polyribosome Metabolism during Normal and Hormone-induced Growth. , 1973, Plant physiology.