Summary

RIBO-seq in Bacteria: a Sample Collection and Library Preparation Protocol for NGS Sequencing

Published: August 07, 2021
doi:

Summary

Here we describe the stages of sample collection and preparation for RIBO-seq in bacteria. Sequencing of the libraries prepared according to these guidelines results in sufficient data for comprehensive bioinformatic analysis. The protocol we present is simple, uses standard laboratory equipment and takes seven days from lysis to obtaining libraries.

Abstract

The ribosome profiling technique (RIBO-seq) is currently the most effective tool for studying the process of protein synthesis in vivo. The advantage of this method, in comparison to other approaches, is its ability to monitor translation by precisely mapping the position and number of ribosomes on a mRNA transcript.

In this article, we describe the consecutive stages of sample collection and preparation for RIBO-seq method in bacteria, highlighting the details relevant to the planning and execution of the experiment.

Since the RIBO-seq relies on intact ribosomes and related mRNAs, the key step is rapid inhibition of translation and adequate disintegration of cells. Thus, we suggest filtration and flash-freezing in liquid nitrogen for cell harvesting with an optional pretreatment with chloramphenicol to arrest translation in bacteria. For the disintegration, we propose grinding frozen cells with mortar and pestle in the presence of aluminum oxide to mechanically disrupt the cell wall. In this protocol, sucrose cushion or a sucrose gradient ultracentrifugation for monosome purification is not required. Instead, mRNA separation using polyacrylamide gel electrophoresis (PAGE) followed by the ribosomal footprint excision (28-30 nt band) is applied and provides satisfactory results. This largely simplifies the method as well as reduces the time and equipment requirements for the procedure. For library preparation, we recommend using the commercially available small RNA kit for Illumina sequencing from New England Biolabs, following manufacturer's guidelines with some degree of optimization.

The resulting cDNA libraries present appropriate quantity and quality required for next generation sequencing (NGS). Sequencing of the libraries prepared according to the described protocol results in 2 to 10 mln uniquely mapped reads per sample providing sufficient data for comprehensive bioinformatic analysis. The protocol we present is quick and relatively easy and can be performed with standard laboratory equipment.

Introduction

The ribosome profiling technique (RIBO-seq) was developed in the laboratory of Jonathan Weissman at the University of California, San Francisco1. In comparison to other methods used to study gene expression at the translational level, RIBO-seq focuses on each ribosome binding to mRNA and provides information about its location and the relative number of ribosomes on a transcript. It enables monitoring the process of protein synthesis in vivo and can provide single codon resolution and accuracy allowing the measurement of the ribosome density on both, the individual mRNA and along the entire transcriptome in the cell. At the foundation of the RIBO-seq technique lies the fact that during translation the ribosome binds the mRNA molecule and thus protects the buried fragment of the transcript from a ribonuclease digestion. Upon addition of the ribonuclease, the unprotected mRNA is digested and the fragments enclosed by ribosomes – typically of ~28-30 nt long – remain intact. These fragments, called ribosomal footprints (RF), can then be isolated, sequenced and mapped onto the transcript they originated from resulting in the detection of the exact position of the ribosomes. In fact, the ribosome ability to protect mRNA fragments has been used since the 1960s to study ribosomal binding and translation initiation sites (TIS)2,3,4. However, with the advancement in deep sequencing technology, RIBO-seq has become a gold standard for translation monitoring5 which, through the ribosome engagement, can provide a genome-wide information on protein synthesis6. Ribosome profiling filled the technological gap that existed between quantifying the transcriptome and the proteome6.

To conduct ribosome profiling we need to obtain cell lysate of the organism that had grown under the investigated conditions. Disrupting these conditions during cell collection and lysis may provide unreliable data. To prevent this, translation inhibitors, rapid harvesting and flash freezing in liquid nitrogen are commonly used. Cells can be lysed by cryogenic grinding in a mechanical homogenizer like a mixer mill7,8 or a bead beater9, and by trituration through a pipette10 or with a needle11. The lysis buffer can be added just before or shortly after pulverization of the cells. In our protocol we use liquid nitrogen to precool mortar and pestle, as well as aluminum oxide as a gentler approach to disruption of the bacterial cell wall, which prevents RNA shearing often encountered when methods such as sonification are applied. After pulverization, we add an ice-cold lysis buffer into the cooled contents of the mortar. Selection of an appropriate lysis buffer is important for obtaining the best resolution of ribosomal footprints. Since ionic strength affects both the RF size and the reading frame precision, it is currently recommended to use lysis buffers with low ionic strength and buffer capacity, even if it appears that buffer composition does not affect ribosomal occupancy on mRNAs11,12. Important components of the lysis buffer are magnesium ions, the presence of which prevents dissociation of the ribosomal subunits and inhibits spontaneous conformational changes in the bacterial ribosomes11,13. Calcium ions also play a significant role and are essential for the activity of micrococcal nuclease (MNase) used in the bacterial ribosome profiling method14. Addition of guanosine 5′-[β,γ-imido]triphosphate (GMP-PNP), a non-hydrolyzable analog of GTP, together with chloramphenicol inhibits translation during lysis15.

When the lysate is obtained, it is clarified by centrifugation and divided into two portions, each for a RIBO-seq and a high-throughput total mRNA sequencing (RNA-seq) since they are performed simultaneously (Figure 1). RNA-seq provides a point of reference which enables the comparison of data from both RIBO-seq and RNA-seq during data analysis. The investigated translatome is defined by normalization of ribosomal footprints to mRNA abundance16. Data from RNA-seq can also help identify cloning or sequencing artifacts17.

Figure 1
Figure 1. Schematics of mRNA sample preparation for RIBO-seq and RNA-seq. For RIBO-seq library preparation, RNA is digested with MNase (A), followed by the size selection of RF of ~28-30 nt length (B); for RNA-seq RNA is isolated (C), depleted of rRNA (D), and the resulting mRNA is randomly fragmented into fragments of varying lengths (E). Please click here to view a larger version of this figure.

Initial steps of the procedure of sample preparation for RIBO-seq and RNA-seq differ slightly (Figure 1). For the ribosomal profiling, the lysate needs to be digested by a specific endonuclease to degrade the mRNA molecules not protected by the ribosomes. In standard protocols, the obtained monosomes are recovered by a sucrose cushion ultracentrifugation or a sucrose gradient ultracentrifugation8,14. In this article, we show that this step is not necessary to isolate RF required for the RIBO-seq in bacteria, likewise for eukaryotic cells18, and that size selection of the appropriate length mRNA fragments from the polyacrylamide gel is sufficient.

For RNA-seq, mRNA is obtained by the depletion of rRNA from the total RNA – rRNA molecules hybridize to the biotinylated oligonucleotide probes which bind to the streptavidin-coated magnetic beads. The rRNA-oligonucleotide-bead complexes are then removed from the sample with a magnet resulting in a rRNA depleted sample19,20. The purified mRNA molecules are then randomly fragmented by alkaline hydrolysis. The obtained fragments of mRNA as well as the ribosomal footprints are converted into cDNA libraries and prepared for deep sequencing (Figure 2). This involves ends repair needed after alkaline hydrolysis (for mRNA) and enzymatic digestion (for RF): dephosphorylation of 3' ends followed by phosphorylation of 5' ends. The next steps are adaptors ligation and the reverse transcription to create cDNA inserts framed by sequences required for the next generation sequencing (NGS) using Illumina platform. The last phase of library preparation is a PCR reaction in which the constructs are amplified and labelled with sample specific barcodes to allow multiplexing and sequencing various samples on one channel. Before sequencing, the quality and quantity of the libraries are assessed by the high-sensitivity DNA on-chip electrophoresis. cDNA libraries with appropriate parameters can then be pooled and sequenced. Sequencing can be performed on different Illumina platforms, such as MiSeq, NextSeq or HighSeq, depending on the number of libraries, required read length and sequencing depth. After sequencing, the bioinformatic analysis is performed.

Figure 2
Figure 2. Library preparation. Library preparation includes the ends repair, adapters ligation, reverse transcription and amplification with barcoding. Please click here to view a larger version of this figure.

The ribosome profiling is a universal method which can be easily modified and adjusted according to the scientific question. Originally it was used in yeast1, but shortly after it was applied to bacterial cells21 as well as eukaryotic model organisms including mouse10, zebrafish22, fruit fly23 and Arabidopsis thaliana24. It was also used for studying different ribosome types: cytoplasmic, mitochondrial25,26 and chloroplast27,28. In eukaryotes RIBO-seq is commonly adapted and refined to investigate specific aspects of translation, including initiation10,11,29,30,31,32, elongation1,10,11,31,33, ribosome stalling33 and conformation change33. Most of the modifications involve the use of different translation inhibitors. In bacteria however, analogous studies have been difficult to conduct because of the paucity of inhibitors with the required mechanism of action34. The most commonly used translation inhibitor in bacteria is chloramphenicol (CAM) which binds to the peptidyl transferase center (PTC) and prevents correct positioning of the aminoacyl-tRNA in the A-site. As a result, CAM prevents the formation of a peptide bond which leads to arresting the elongating ribosomes35. Other examples of translation inhibitors in bacteria are tetracycline (TET)36, retapamulin (RET)34 and Onc11237 which have been used to investigate translation initiation sites. TET, which prevents tRNA delivery to the ribosome by directly overlapping with the anticodon stem-loop of tRNA at the A-site, was originally applied to verify the results obtained from CAM treatment since they are both antibiotics inhibiting translation elongation38. TET was found to detect primary TIS, however was unable to reveal internal TIS36. RET binds in the PTC of the bacterial ribosome, and prevents formation of the first peptide bond by interfering with an elongator aminoacyl-tRNA in the A site. Applying RET results in ribosomes arrest at both primary as well as internal TISs34. Onc112, a proline-rich antimicrobial peptide, binds in the exit tunnel and blocks aminoacyl-tRNA binding in the ribosomal A site. As a result, Onc112 prevents initiation complexes from entering the elongation phase37

The main information ribosome profiling provides is ribosomes density and their position on the mRNA. It was successfully applied to investigate differential gene expression at the level of translation in various growth conditions1,6, measure translational efficiency1,38,39 and detect translation regulation events such as ribosomal pausing10. RIBO-seq also allows for uncovering the translation of annotated ncRNA, pseudogenes and unannotated small open reading frames (ORF) leading to the identification of novel and/or very short protein coding genes10,12,22,30,37. In such cases, RIBO-seq can fine-tune and improve genome annotation. With its high sensitivity for the identification of translated ORFs and its quantitative nature, ribosome profiling can also serve as a proxy for the proteome determination or in aiding proteomics studies31,34,39. By mapping TIS, ribosome profiling reveals N-terminally extended and truncated isoforms of known proteins10,32. RIBO-seq was also adapted to study co-translational folding of proteins14,21,24. This method enables measuring of elongation rates1,10,39 or decoding speeds of individual codons6 and helps in developing quantitative models of translation17. The ribosome profiling method is also capable of providing mechanistic insights into the ribosome pausing in bacteria7,15,17, frameshifting40, stop-codon readthrough21, termination/recycling defects41,42 and ribosomal conformation changes33 in eukaryotes. RIBO-seq was also adapted to examine the impact of specific trans-acting factors on translation such as miRNAs6 and RNA-binding proteins in eukaryotes16,43. However, it is important to acknowledge that the experimental design and the obtained resolution of RIBO-seq determine the amount of information that can be extracted from the resulting sequencing data12.

Protocol

1. Sample collection

  1. Prepare a bacterial culture. We recommend a culture volume of 100 mL per sample.
  2. Prepare equipment and reagents for sample collection: two scoopulas per sample, sterile 0.45 µm mixed cellulose esters membrane (MCE) filters, 50 mL sterile tubes, liquid nitrogen, 50 mg/mL chloramphenicol in 70% (vol/vol) ethanol (optional). Puncture the lid of 50 mL tubes to allow liquid nitrogen evaporation and prevent the explosion of closed tubes. Decontaminate scoopulas with laboratory detergent and 70% ethanol.
  3. Prewarm the filtering equipment and one of the scoopulas to growth temperature of the bacterial culture.
    NOTE: We recommend an all-glass filtration apparatus with a funnel, a fritted base, a ground joint flask and a 47 mm Ø spring clamp.
  4. Before sample harvesting, add chloramphenicol into the bacterial culture to the final concentration of 100 µg/mL and incubate 1 minute (optional).
  5. Pour liquid nitrogen into 50 mL sterile tube and place the second scoopula inside the tube to cool.
    CAUTION! Liquid nitrogen can cause closed containers to explode due to the pressure change upon evaporation. To prevent this, it is necessary to create a few holes in the lid of 50 mL tubes to allow liquid nitrogen evaporation.
  6. Collect cells by filtration. Stop filtering when medium has passed through the membrane, but do not allow the filter to dry completely.
  7. Collect bacterial pellets by rapidly scraping the cells off of the filter disc using a prewarmed scoopula. Immediately place the entire scoopula with the harvested cells in the 50 mL tube filled with liquid nitrogen. The harvested pellet should be completely covered in liquid nitrogen.
  8. Let the pellet freeze thoroughly and dislodge frozen cells using a previously cooled second scoopula. Close the punctured lid and transfer the tube to -80°C. STOP point
    ​NOTE: Disinfect scoopulas after each sample harvesting.

2. Cell lysis

  1. Prepare 0.1 M GMP-PNP in 10 mM Tris pH 8, DNase I and 1.5 mL reaction tubes for lysates and keep them on ice. Cool the centrifuge to 4 °C.
  2. Prepare lysis buffer (20 mM TRIS pH 7.6, 10 mM MgAcet, 150 mM KAcet, 0.4% TRITON X-100, 6 mM β-mercaptoethanol, 5 mM CaCl2 and optional 1 mM chloramphenicol). Aliquot 500 µL of the lysis buffer per sample into 1.5 mL reaction tubes and place them on ice.
  3. Decontaminate mortar and pestle with a laboratory detergent and 70% ethanol, and dry them. Cool mortar and pestle by pouring liquid nitrogen into the mortar.
  4. Transfer frozen pellet into the pre-chilled mortar and grind it to a powder. With a spatula, add approximately 1 volume of aluminum oxide and continue grinding. Keep mortar, pestle and cells cool by pouring liquid nitrogen when needed, do not let the contents of the mortar to thaw.
  5. Just before using lysis buffer, add GMP-PNP and DNase I into the lysis buffer aliquot to the final concentration of 2 mM and 100 U/mL respectively. Transfer the solution into the mortar with cells and aluminum oxide and continue grinding. Let the lysate thaw slowly while grinding and transfer the mixture into the pre-cooled 1.5 mL reaction tube and immediately place on ice.
  6. Centrifuge lysates at 20 000 x g for 5 minutes at 4 °C. Transfer supernatants into new, pre-cooled 1,5 mL reaction tubes and keep them on ice.
  7. Measure the RNA concentration in each sample with a NanoDrop. Use 1:10 dilutions of samples and 1:10 dilution of lysis buffer in nuclease-free water as a blank.
    NOTE: Be aware that chloramphenicol shows significant absorption at 260 nm.
  8. Divide each lysate into two portions: one for RIBO-seq (0.5 – 1 mg of RNA) and the second for RNA-seq (the rest).
  9. Clean the samples for RNA-seq using a commercially available RNA clean-up kit according to the manufacturer's protocol.
    ​NOTE: We recommend using Zymo RNA Clean & Concentrate -25 kit (see Table of Materials). We also recommend adding 4.5 volumes of ethanol (compared to sample volume) into the mix of sample and Binding Buffer for ribosomal footprints and fragmented mRNA, in order to increase the purification efficiency of short RNA fragments.
  10. Store the samples intended for RNA-seq at -80 °C. The procedure for RNA-seq will continue from step 5.

3. MNase digestion of samples for RIBO-seq

  1. To 1 mg of RNA add 3.8 µL of 187.5 U/µL MNase in 10 mM Tris pH 8 and top it up with lysis buffer to the total volume of 500 µL.
  2. Incubate at 25 °C, 300 RPM, for 45 minutes in a thermomixer.
  3. Clean the samples with a commercially available RNA clean up kit (as in step 2.9).
  4. Store the samples for RIBO-seq at -80 °C (STOP point) or proceed to size selection.

4. Polyacrylamide gel electrophoresis (PAGE) and size selection of samples for RIBO-seq

  1. Prepare a 15% polyacrylamide-TBE gel with 8 M urea and place it in a tank with TBE buffer. Pre-run for minimum 10 minutes at a constant voltage of 200 V.
  2. Mix samples with TBE-Urea Sample Buffer, denature at 95 °C for 1 minute and place them immediately on ice.
  3. Wash out the urea by injecting TBE buffer into gel wells using a syringe. Load the samples leaving one well space between them to separate each sample and prevent cross-contamination. Use 29 nt oligonucleotide and a mix of 26 nt and 32 nt oligonucleotides as markers. Run electrophoresis at a constant voltage of 180 V.
  4. Prepare sterile buffer for overnight incubation (5 mM EDTA, 10 mM NaOAc pH 5).
  5. After electrophoresis, stain the gel in a SYBR Gold-bath for about 2-3 minutes and rinse the gel with nuclease-free water.
  6. Excise fragments of gel between 26 and 32 nt with the sterile needles or the razor blades and place the gel fragments in separate 1.5 mL tubes for each sample. Change the needle or the razor blade between the samples.
  7. Add 200 µL of the overnight incubation buffer to each reaction tube.
  8. Incubate the samples at 10 °C, 1000 RPM overnight in a thermomixer.
  9. The next day, clean the samples as in step 2.9. Elute the footprints in 80 µL of nuclease-free water.
  10. Store the samples at -80 °C. STOP point. The procedure for RIBO-seq will continue from step 7.

5. Purification of bacterial mRNA by rRNA depletion from samples for RNA-seq

  1. Remove rRNA from samples with a bacterial mRNA purification kit (e.g., Invitrogen MICROBExpress). Follow the manufacturer's protocol.
  2. Clean the samples as in step 2.9. Elute the mRNA in 50 µL of nuclease-free water.
  3. Store the samples for RNA-seq at -80 °C (STOP point) or continue to alkaline fragmentation.

6. Alkaline fragmentation of samples for RNA-seq

  1. Prepare alkaline hydrolysis buffer (mix 220 µL of 0.1 M NaHCO3, 30 µL of 0.1 M Na2CO3 and 1 µL of 0.5 M EDTA). Mix 1 volume (50 µL) of alkaline buffer with 1 volume (50 µL) of the sample and incubate at 95 °C for 25 minutes.
  2. Add 5 µL of 3 M NaOAc pH 5.5 to stop the reaction.
  3. Clean the samples as in step 2.9 and elute the samples in 80 µL of nuclease-free water.
  4. Possible STOP point: store RNA-seq samples at -80 °C. However, we recommend to go directly to step 7 in order to avoid unnecessary freezing and thawing.

7. Dephosphorylation and phosphorylation of samples for both RNA- and RIBO-seq

  1. Add 10 µL of 10x reaction buffer PNK and 5 µL T4 PNK to each sample. Incubate at 37 °C for 1.5 hour in order to dephosphorylate the 3' ends.
  2. Add 3 µL of 1 mM ATP and incubate at 37 °C for 1 hour in order to phosphorylate the 5' ends.
  3. Clean the samples as in step 2.9, elute in 30 µL nuclease-free water.
  4. Store the samples at -80 °C. STOP point.

8. Library preparation using NEBNext Multiplex Small RNA Library Prep Set for Illumina

  1. Prepare 100-1000 ng of input RNA (obtained mRNA fragments and ribosomal footprints) in 6 µL of nuclease-free water and add 1 µL of 3' SR Adaptor. Incubate according to the manufacturer's protocol.
  2. Add 10 µL of 3' Ligation Reaction Buffer (2X) and 3 µL of 3' Ligation Enzyme Mix and incubate at 37 °C for 2.5 hour, instead of 1 hour at 25 °C as the manufacturer's protocol specifies.
  3. Hybridize the Reverse Transcription Primer according to the manufacturer's protocol with a modified program: 5 minutes at 75 °C, 30 minutes at 37 °C and hold at 4 °C.
  4. Ligate the 5´ SR Adaptor. Follow the manufacturer's protocol with one change: perform the incubation at 37 °C for 2.5 hour, instead of 1 hour at 25 °C.
  5. Perform Reverse Transcription according to the manufacturer's protocol.
  6. Possible STOP point: incubate the samples containing synthesized cDNAs at 70 °C for 15 minutes to inactivate the reverse transcriptase and store them at -80 °C. However, we recommend to proceed directly to the next steps to avoid unnecessary freezing and thawing of the samples.
  7. Store half of each cDNA library at -80 °C as a back-up.
  8. Perform PCR amplification according to the manufacturer's protocol. Use half of the reaction mix. Store the obtained libraries at -80 °C. STOP point.

9. Size selection of cDNA libraries using PAGE

  1. Prepare 6% polyacrylamide-TBE gel and place it in a tank with TBE buffer. Pre-run for minimum 10 minutes at a constant voltage of 200 V.
  2. Mix the samples with Gel Loading Dye, Blue (6X) and load them on the gel. Use Quick-Load pBR322 DNA-MspI Digest as a ladder. At the beginning, run electrophoresis at a constant voltage of 120 V to allow DNA to enter the gel and then change the voltage to 180 V.
  3. When electrophoresis is finished, stain the gel in a SYBR Gold-bath for about 2-3 minutes and rinse the gel with nuclease-free water.
  4. Excise gel fragments containing the libraries. For RNA-seq samples excise between 135-180 nt and for RIBO-seq between 135-170 nt. Use sterile needle or razor blade and place the excised gel fragments in separate 1.5 mL reaction tubes. Remember to change the needle or razor blade between the samples.
    NOTE: The adapters and barcode from NEBNext Multiplex Small RNA Library Prep Set for Illumina have 119 nt in total which determines the lower excision cut-off.
  5. Add 100 µL of nuclease-free water to each excised gel fragment.
  6. Incubate the gel fragments at 10 °C, 450 RPM overnight in a thermomixer.
  7. Clean and concentrate the cDNA libraries with a commercial DNA clean-up kit according to the manufacturer's protocol.
    NOTE: We recommend using Zymo DNA Clean & Concentrate -5 kit (see Table of Materials).
  8. Store purified cDNA libraries at -80 °C. STOP point.

Representative Results

The exemplary results presented here were obtained in a study examining translation regulation in sporulating WT Bacillus subtilis cells. Overnight cultures were diluted to OD600 equal to 0.1 in 100 mL of rich medium and incubated at 37 °C with vigorous shaking until OD600 reached 0.5-0.6. The rich medium was then replaced with minimal medium to induce sporulation process and the incubation was continued for up to four hours. Cells were harvested every hour beginning with T0 – sporulation induction – resulting in five time points. One minute before harvesting, the cultures were treated with chloramphenicol as described in the protocol above. Cells were collected by filtration in a prewarmed glass filtration system using 0.45 µm mixed cellulose esters membrane (MCE) filters and flash-frozen in liquid nitrogen.

The amount of ribonucleic acid obtained in the lysate depends on growth conditions, culture volume and density. Following our protocol, 100 mL of the bacterial culture (Bacillus subtilis) yields 1-4 mg of RNA on average. Examples of lysates obtained from the experiment described above are shown in Table 1.

RNA concentration [ng/µL] RNA amount [mg] A260/A280 A260/A230
WT0 5845 3.21 1.57 0.83
WT1 6275 3.45 1.58 0.79
WT2 3525 1.94 1.38 0.56
WT3 3717 2.04 1.37 0.61
WT4 1653 0.91 1.11 0.43

Table 1. RNA concentration and amount in representative samples of B. subtilis lysates. The table includes samples from the wild type (WT) B. subtilis which were harvested before the sporulation induction (0) and one hour (1), two (2), three (3) or four (4) hours post sporulation induction. Lysis was performed with 550 µL of lysis buffer containing chloramphenicol. It can be noticed that as the sporulation proceeds the amount of obtained RNA decreases.

When the amount of RNA is lower than 1 mg, 0.25 mg of RNA is used for RIBO-seq instead of 1-0.5 mg and the amount of MNase is then adjusted respectively. The representative RNA amount and concentration after nuclease digestion is shown in Table 2. The average amount of RNA after digestion and cleaning is 50 µg from 0.5 mg of RNA input.

RNA concentration [ng/µL] RNA amount [µg] A260/A280 A260/A230
WT0-R 1015 50.75 2.12 1.84
WT1-R 1157.4 57.87 2.10 2.19
WT2-R 847.2 42.36 2.10 2.00
WT3-R 983.6 49.18 2.09 2.16
WT4-R 206 10.30 2.13 1.00

Table 2. The amount and concentration of RNA after nuclease digestion of the representative samples for RIBO-seq. The table includes samples from the wild type (WT) B. subtilis which were harvested before sporulation induction (0) and one hour (1), two (2), three (3) or four (4) hours post sporulation induction; the letter "R" refers to samples for RIBO-seq. Because RNA concentration of WT4-R sample was low, 0.25 mg of RNA was digested, while 0.5 mg of the RNA for the remaining samples was used for ribosomal footprint generation. The amount of MNase for digestion was 178.125 U for 0.25 mg of RNA in WT4-R, and 356.25 U for 0.5 mg of RNA in other samples. After digestion the samples were cleaned with the commercial RNA clean-up kit and ribonucleic acid concentration was measured with the NanoDrop.

After MNase digestion, size selection is performed with the aid of PAGE. An example of a 15% TBE-urea polyacrylamide gel with the digested samples is shown in Figure 3. Fragments of gel between 26 and 32 nt were excised with a sterile razor blade and incubated overnight in the overnight incubation buffer. The next day, ribosomal footprints were cleaned with the commercial RNA clean-up kit.

Figure 3
Figure 3. 15% TBE-urea polyacrylamide gel with representative samples after MNase digestion. The figure includes samples from the wild type (WT) B. subtilis which were harvested before sporulation induction (0) and one hour (1), two (2), three (3) or four (4) hours post sporulation induction; the letter "R" refers to samples for RIBO-seq. 15 µL of denatured samples were loaded into the gel wells in order: WT2-R, WT3-R, WT4-R, WT0-R, WT1-R and the remaining volume was stored at -80 °C as a back-up. 29 nt oligonucleotide and mix of 26 nt and 32 nt oligonucleotides were used as markers. Electrophoresis was performed as described above and the gel was stained in a SYBR Gold-bath. Please click here to view a larger version of this figure.

Parallel to nuclease digestion, the lysates for RNA-seq were cleaned with the commercial RNA clean-up kit and the obtained concentrations of RNA were measured with NanoDrop. The representative results are presented in Table 3. After cleaning, the amount of ribonucleic acid decreased to 40-500 µg.

RNA concentration [ng/µL] RNA amount [µg] A260/A280 A260/A230
WT0-m 10274.1 513.70 2.14 2.35
WT1-m 10118.3 505.92 2.12 2.33
WT2-m 5599.7 279.98 2.15 2.29
WT3-m 6562.2 328.11 2.08 2.31
WT4-m 757 37.85 2.01 1.95

Table 3. RNA amount and concentration of samples for RNA-seq after RNA purification. The table includes samples from the wild type (WT) B. subtilis which were harvested before sporulation induction (0) and one hour (1), two (2), three (3) or four (4) hours post sporulation induction; the letter "m" refers to samples for RNA-seq (mRNA).

The obtained RNA concentrations of the samples for RNA-seq were used to define RNA input for rRNA depletion. 8 µg of RNA from each RNA-seq sample were used with MICROBExpress bacterial mRNA purification kit (see Table of Materials) according to the manufacturer's protocol and were cleaned with the commercial RNA clean-up kit afterwards. The representative results of rRNA depletion are shown in Table 4.

RNA concentration [ng/µL] A260/A280 A260/A230 RNA amount [µg]
WT0-m 62.7 1.99 1.75 3.1
WT1-m 59.3 2.00 2.16 2.9
WT2-m 93.7 1.97 1.69 4.7
WT3-m 68.3 1.99 2.33 3.4
WT4-m 85.0 1.97 2.00 4.2

Table 4. The representative RNA amount and concentration after rRNA depletion. The table includes samples from the wild type (WT) B. subtilis which were harvested before sporulation induction (0) and one hour (1), two (2), three (3) or four (4) hours post sporulation induction; the letter "m" refers to samples for RNA-seq. 8 µg of RNA from samples for RNA-seq were used for rRNA depletion and cleaned with a commercial RNA clean-up kit afterwards.

Next, the obtained mRNA was fragmented by alkaline hydrolysis and cleaned with a commercial kit as described in the protocol. The fragmented mRNA and ribosomal footprints were dephosphorylated at 3' ends and phosphorylated at 5' ends according to the above protocol. The samples were then used to construct cDNA libraries for Illumina sequencing, as described in the protocol. Figure 4 shows an example of a 6% polyacrylamide gel with the obtained libraries. Gel fragments in the range of 135-180 nt for RNA-seq and 135-170 nt for RIBO-seq were excised with sterile razor blades and incubated overnight in nuclease-free water. Libraries were cleaned with a commercial DNA clean-up kit and quality and quantity was assessed by a high-sensitivity DNA on-chip electrophoresis using the Agilent 2100 Bioanalyzer.

Figure 4
Figure 4. 6% polyacrylamide gel of cDNA libraries before and after size selection. Quick-Load pBR322 DNA-MspI Digest was used as a ladder. 25 µL of samples were loaded into the gel in the following order: three RIBO-seq libraries and six RNA-seq libraries. The remaining volumes of samples were stored at -80 °C as a back-up. The figure includes samples from the wild type (WT) B. subtilis which were harvested one hour (1), two (2), three (3) or four (4) hours post sporulation induction; the letter "R" refers to samples for RIBO-seq and the letter "m" refers to samples for RNA-seq. Please click here to view a larger version of this figure.

The examples of virtual gel images and electropherograms obtained from the Agilent 2100 Bioanalyzer are shown in Figure 5. Bands and peaks representing the RIBO-seq libraries are more narrow and better defined compared to these representing RNA-seq libraries. According to the on-chip electrophoresis, libraries are approximately 150 bp long, which is an expected library size. The adapters and the barcode used in library preparation are 119 nt long and the inserts are either 29-30 nt long (for RIBO-seq) or in a range between 25 and 50 nt (for RNA-seq). The additional peaks closer to the 200 bp obtained from RIBO-seq libraries may indicate PCR artefacts which can be discarded during bioinformatic analysis when the data obtained from sequencing is trimmed and rRNA/tRNA is filtered. The average quantity of DNA is 32 ng providing sufficient quantity of material required for the next generation sequencing.

Figure 5
Figure 5. The examples of electropherograms and virtual gel images from a high-sensitivity DNA on-chip electrophoresis. The figure includes samples from the wild type (WT) B. subtilis which were harvested before sporulation induction (0) and one hour (1), two (2), three (3) or four (4) hours post sporulation induction; the letter "R" refers to RIBO-seq libraries and the letter "m" refers to RNA-seq libraries. Please click here to view a larger version of this figure.

After quality control by the on-chip electrophoresis, libraries were pooled for single-end 50 bp sequencing on the Illumina's NextSeq500 platform (10 libraries per channel, with a minimum of 350 mln reads per channel). Sequencing yielded 30-57 million reads per sample for RNA-seq libraries and 30-50 million reads per sample for RIBO-seq libraries. Trimming the adaptors and poor quality sequences resulted in 24-47 mln reads per sample for RNA-seq samples and 25-50 mln reads per sample for RIBO-seq samples. Mapping, performed with STAR44, yielded 2.4-9.6 million of uniquely mapped reads per sample for RNA-seq samples (which constitute 8-16.8% of the total number of reads) and 2.3-10.4 million uniquely mapped reads (7.7-22.1%) for RIBO-seq samples. The obtained data was quality control checked with FastQC45. Examples of the quality plots produced by FastQC for the trimmed sequences are shown in Figure 6. The sequences of the length of interest, which is <32 nt for RIBO-seq and <50 nt for RNA-seq, presented very good quality. The examples of plots of GC content and distributions of sequence lengths over all trimmed sequences are shown in Figure 7. The additional peak at 72% in GC distribution plot of RIBO-seq library may potentially indicate rRNA contamination which can be removed during bioinformatic analysis by rRNA filtration46,47,48. The sequence lengths distribution is very narrow for RIBO-seq libraries, with a peak at 26-27 nt whereas the length distribution for RNA-seq libraries is broader and ranges between 15 and 50 nt.

Figure 6
Figure 6. The examples of box-and-whisker plots of quality scores across all bases for RIBO-seq and RNA-seq libraries after trimming. The figure includes samples from the wild type (WT) B. subtilis which were harvested before sporulation induction (0) and four (4) hours post sporulation induction; the letter "R" refers to RIBO-seq samples and the letter "m" refers to RNA-seq samples. Please click here to view a larger version of this figure.

Figure 7
Figure 7. The examples of plots of GC content and sequence length distribution of all sequences for RIBO-seq and RNA-seq libraries after trimming. The figure shows samples from the wild type (WT) B. subtilis which were harvested four (4) hours post sporulation induction; the letter "R" refers to sample for RIBO-seq and the letter "m" refers to sample for RNA-seq. Please click here to view a larger version of this figure.

To verify whether the RIBO-seq data provides information on true ribosomal footprints rather than random mRNA fragments of the selected length, we compared sequencing data from RIBO-seq and RNA-seq using RiboToolkit, a web server for RIBO-seq data analysis49. Ribo-seq data exhibits triplet periodicity and a tall, narrow peak corresponding to the initiating ribosomes and characteristic of the RF as shown in Figure 8A, which is not observed in the RNA-seq data (Figure 8B). Moreover, the metagene plot showing the profiles of average coverage of RF and mRNA fragments on the coding sequences (CDSs) revealed a higher proportion of reads mapped to the CDSs in the RIBO-seq dataset (Figure 8C), as expected.

Figure 8
Figure 8. Metagene periodicity plots and average coverage of RF and mRNA fragments on CDSs. Metagene profiles were obtained with RiboToolkit, an integrated web server49. CDS length was normalized to a value of 1 and each CDS was divided into 100 equal bins. The figure shows samples from the wild type (WT) B. subtilis which were harvested one (1) hour post sporulation induction; the letter "R" refers to sample for RIBO-seq and the letter "m" refers to sample for RNA-seq. Please click here to view a larger version of this figure.

To compare whether the RIBO-seq protocol presented here yields results similar to those obtained by standard method of RF recovery (sucrose gradient ultracentrifugation), we paralleled the sequencing results from libraries prepared following the two methods. The experiments investigating translation regulation during sporulation in the WT B. subtilis were conducted according to the presented RIBO-seq protocol, as well as following the standard protocol which includes monosomes recovery by a sucrose gradient ultracentrifugation. The results of ribosome coverage profiles obtained from the two experiments are shown in Figure 9. The visualisation of mapped RFs is very similar between the two experiments and yielded similar results upon bioinformatic analysis.

Figure 9
Figure 9. The visualisation of obtained RF mapped onto the genome. The figure includes RF obtained from wild type B. subtilis harvested four hours post sporulation induction and prepared according to the standard protocol (WT4-R-s, upper panel) or the presented protocol (WT4-R, lower panel). Please click here to view a larger version of this figure.

Discussion

The key technical challenge of the ribosome profiling is the need to rapidly inhibit translation in order to capture a snapshot of ribosomes on mRNAs at a particular physiological state of interest. To accomplish this, translation inhibitors, rapid harvesting and flash freezing in liquid nitrogen are commonly used. Applying antibiotics is optional since they can cause artifacts. Chloramphenicol is a commonly used drug to arrest elongating ribosomes in bacterial RIBO-seq. However, it does not prevent initiation, resulting in an accumulation of ribosomes about six codons downstream of the translation initiation site14. Moreover, its ability to inhibit the peptidyltransferase reaction can be dependent on the nature of the ribosomal P- and A-site substrates. Presumably, CAM arrests ribosomes more effectively with Gly, Ala, Ser, Cys, or Thr as penultimate amino acid of the nascent peptide17,38,50. This is worth taking under consideration during experimental design and bioinformatic analysis, especially when studying ribosomal pausing17. CAM has also been shown to influence translational accuracy by promoting a stop codon readthrough and frameshifting51. Still, CAM inhibition is relatively fast and sufficient for typical RIBO-seq research14 and does not induce misincorporation51. Some of the artifacts caused by antibiotics can be avoided by increasing the drug amount9. CAM concentration applied in our protocol is the typical amount used in RIBO-seq experiments7,14,21, however, the concentration may vary for different microorganisms and should be optimized before sample harvesting. Due to the possibility of an occurrence of artifacts caused by applying translation inhibitors, it is recommended to avoid the inhibitor treatment if not necessary for the research question and if the rapid harvesting is possible14. Fast cells collection is crucial, especially if it is performed without the inhibitor treatment, since prolonged harvesting can result in a loss of polysomes and/or change in the translatome as a response to the changing environmental conditions14. Thus, we recommend using filtration as a harvesting method as it can be performed under the same or similar conditions as the culture growth, including temperature and/or shaking. During filtration cells collected at the filter surface are constantly flushed with growth medium, impeding them from sensing shifts in cell densities, oxygenation and level of amino acids and other nutrients14. Filtration is gentler and may be faster than centrifugation, if the density of the culture is moderate, otherwise a large number of cells can clog the filter pores and slow down the process. If there is a probability for a prolonged filtration, we recommend a chloramphenicol pre-treatment of the culture, which results in a translation arrest and allows to extend harvesting time. If filtration becomes difficult to complete due to the filter clogging, we suggest to stop the filtration, remove the excess culture and flash freeze fraction of the culture collected on the filter. We observed that 50 mL of dense bacterial culture is sufficient to perform RIBO-seq. The subsequent flash freezing of cells in liquid nitrogen is the most universal and effective approach to immediately stop translation8,21,52.

Ribonuclease digestion is a critical step required to obtain good quality RF. Selecting the right enzyme is important since the ribosomal subunits are composed of RNA molecules that can be digested, disrupting the integrity of the ribosome, causing rRNA contamination and destroying RF44. Originally, RNase I was used for digestion of mRNA and RF generation in eukaryotes1. Since this enzyme seemed inactive in bacteria, it was replaced with a micrococcal nuclease (MNase) from Staphylococcus aureus21. It was later showed that RNAse I can in fact perform digestion of bacterial samples, however MNase remains commonly used for this purpose53. The disadvantage of MNase lies in the sequence bias of its catalytic activity, which is 30-fold increased proximal to A or T, resulting in 80% of mRNA fragments starting with A or T at the 5′ end14. However, this has been overcome by the observation that most of the variation in length occurs at the 5′-end of RF, and assigning ribosome density to the 3′-end yields precise and accurate information about the ribosome position7. The MNase activity is affected by the enzyme concentration, pH and digestion time and should be determined for each new MNase batch. An increased nuclease activity results in increased rRNA contamination, while decreased activity causes less accurate cleavage of ribosomal footprints. This may reduce the overall yield of RF in the size-selective gel purification procedure and provide less accurate information on ribosome positions on mRNA14. Thus, the required amount of MNase should be determined carefully. For our experiments we apply 712.5 U of micrococcal nuclease per 1 mg of RNA while the concentration reported in the literature varies between 50 U38, 450 U34, 1000 U14 and 1500 U15,21,36 per 1 mg of RNA. As mentioned above, MNase requires presence of calcium ions for its activity. In standard protocols EGTA, a chelating agent for divalent ions which has a higher affinity for calcium compared to magnesium54, is used for stopping the digestion by binding calcium and inactivating MNase. However, in our protocol the digestion reaction is stopped by using RNA cleaning kit directly after incubation with the nuclease.

In standard protocols, the next step is ribosomes recovery by a sucrose cushion ultracentrifugation or a sucrose gradient ultracentrifugation8,14. A sucrose gradient centrifugation enables the selective separation of the monosomes from the subunits and polysomes, however it requires specific equipment such as a gradient mixer and an ultracentrifuge. Also, the throughput of this method is limited to 500-700 µL of the lysate. On the other hand, sucrose cushion ultracentrifugation pellets most ribosomal particles, enables processing of over dozen milliliters of lysate and is less demanding in terms of the instrumentation needed. However, both approaches are lengthy as they take minimum 4 to 6.5 hours14 and require an ultracentrifuge. In our protocol this step is omitted as we proceed directly to size selection in 15% TBE-urea polyacrylamide gel. The electrophoresis is performed under denaturing conditions in order to unfold RNA's secondary structures and prevent intramolecular interactions. We use 26 nt, 29 nt and 32 nt oligonucleotides as markers for excising a narrow range of ribosomal footprints. We assume that the concentration of the actual RF significantly exceeds potential contamination with false RF considered as random, unprotected by the ribosome mRNA or rRNA fragments of length in that range. The results obtained using our protocol indeed showed enriched fraction of the RF in the selected size range. The quality check with the Agilent 2100 Bioanalyzer shows RIBO-seq libraries as a single, well define bands and peaks, in contrast to RNA-seq libraries which are visualized as significantly wider smears and peaks (Figure 5). Also, obtained Ribo-seq data exhibits triplet periodicity (Figure 8A), which is not observed in the RNA-seq data (Figure 8B). Moreover, the visualizations of the obtained RF mapped onto the genome show very similar expression patterns regardless of the footprint isolation method which was further confirmed in our bioinformatical analyses (Figure 9). Still, we are aware that false RF could map to the genome which may disturb the obtained signal. However, the fact that false RF are randomly generated implies that they should not accumulate to one place on the genome and hence, should not result in artifacts and false result. Also, it has to be noted that there is no consensus in bacterial studies regarding the size of RFs that should be isolated53,55 and that selecting only a subfraction of length range may lead to misinterpretations of the results17. It is suggested that a wide range of lengths characteristic for bacterial footprints is a consequence of the prokaryotic ribosome behaviour rather than MNase digestion55. For example, complementary interaction between the 3′ end of 16S rRNA in small ribosomal subunit and a Shine-Dalgarno motif in mRNA can result in longer RF56. Considering a wide range of lengths of bacterial footprints, Mohammad et al. recommend to perform size selection between 15 and 40 nt17,55. This should lead to the isolation of potentially an entire range of RFs and ensure obtaining a comprehensive representation of ribosomes in various stages of the translational cycle. Because of the ongoing discussion, the excision range should be considered carefully and it is important to remember that the incorrect choice may exclude some subfractions of RF and lead to bias in the subsequent analysis.

Due to the size selection performed directly after nuclease digestion and with omitting the monosome isolation by ultracentrifugation, our protocol is quick, relatively easy and can be conducted with a standard laboratory equipment. Using dedicated kits for cleaning of the RNA and DNA, purification of mRNA and library preparation ensures standardization and repeatability. Appling more volume of ethanol during the RNA cleaning increases the purification efficiency of <200 nt RNA fragments57. Modifications to the library preparation – extending the incubation time and increasing the temperature during adaptors ligation – are applied to decrease ligation efficiency for methylated RNA, such as rRNA, and to reduce rRNA contamination58. Our protocol has been used to investigate translation regulation in various growth conditions (publications under preparation) but can be applied to study other aspects of translation like TIS detecting or to improve genome annotation and assist proteomic studies. An additional modification to the protocol, dependent on the research question, can be easily included, e.g. crosslinking the ribosome with the protein of interest or affinity purification of the ribosomes of interest resulting in the enriched fractions of ribosomes. With our protocol the RIBO-seq procedure is more accessible for scientists and thus, may stimulate new discoveries and accelerate developments in the field.

Disclosures

The authors have nothing to disclose.

Acknowledgements

ALS would like to acknowledge the financial support of EMBO Installation Grants IG 3914, and POIR. 04.04.00-00-3E9C/17-00 carried out within the First TEAM programme of the Foundation for Polish Science co-financed by the European Union under the European Regional Development Fund.

Materials

10X TBE (powder) Invitrogen AM9864
2-Mercaptoethanol, 99%, pure Acros Organics 125472500
Adenosine 5'-Triphosphate (ATP) New England Biolabs P0756S
Aluminium oxide calcinated pure p.a. Chempur 114560600
Calcium chloride dihydrate Sigma-Aldrich C3881-500G
Chloramphenicol MP Biomedicals 190321
DNA Clean & Concentrator -5 Zymo Research D4004
Dnase I recombinant, Rnase-free Roche 4716728001
EDTA disodium salt Fisher Scientific E/P140/48
Ethyl Alcohol Absolut 99,8%  Pure-P.A.-Basic POCH Avantor Performance Materials Poland S.A BA6480111
Filtration apparatus VWR Collection 511-0265 all-glass filtration apparatus, with funnel, fritted base, cap, 47 mm Ø spring clamp and ground joint flask
Gel 40 (19:1) Rotiphorese 3030.1
Gel Loading Dye, Blue, 6X New England Biolabs E6138G
Guanosine 5′-[β,γ-imido]triphosphate trisodium salt hydrate Sigma-Aldrich G0635-25MG
labZAP A&A Biotechnology 040-500
Magnesium acetate tetrahydrate Sigma-Aldrich M5661-250G
MCE membrane fiter Alfatec Technology M47MCE45GWS pore size: 0.45um
MICROBExpress Bacterial mRNA Purification Invitrogen AM1905
Multiplex Small RNA Library Prep Set for Illumina New England Biolabs E7300S
Nuclease-Free Water Ambion AM9937
Potassium Acetate Anhydrous Pure P.A. POCH Avantor Performance Materials Poland S.A 744330113
Quick-Load pBR322 DNA-MspI Digest New England Biolabs E7323A
RNA Clean & Concentrator -25 Zymo Research R1018
Sodium acetate Sigma-Aldrich S2889-250G
Sodium carbonate Sigma-Aldrich 223530-500G
Sodium hydrogen carbonate pure p.a. POCH Avantor Performance Materials Poland S.A 810530115
SYBR Gold nucleic acid gel stain Life Technologies S11494
T4 Polynucleotide Kinase New England Biolabs M0201L
T4 Polynucleotide Kinase Reaction Buffer New England Biolabs B0201S
TBE-Urea Sample Buffer (2x) Invitrogen LC6876
Tris(hydroxymethyl)amino-methane, ultrapure, 99,9% AlfaAesar J65594
Triton X-100, 98% Acros Organics 327371000
Urea G.R. lach:ner 40096-AP0

References

  1. Ingolia, N. T., Ghaemmaghami, S., Newman, J. R., Weissman, J. S. Genome-wide analysis in vivo of translation with nucleotide resolution using ribosome profiling. Science. 324 (5924), 218-223 (2009).
  2. Takanami, M., Yan, Y., Jukes, T. H. Studies on the site of ribosomal binding of f2 bacteriophage RNA. Journal of Molecular Biology. 12 (3), 761-773 (1965).
  3. Steitz, J. A. Polypeptide Chain Initiation: Nucleotide Sequences of the Three Ribosomal Binding Sites in Bacteriophage R17 RNA. Nature. 224 (5223), 957-964 (1969).
  4. Wolin, S. L., Walter, P. Ribosome pausing and stacking during translation of a eukaryotic mRNA. The EMBO journal. 7 (11), 3559-3569 (1988).
  5. Argüello, R. J., et al. SunRiSE – measuring translation elongation at single-cell resolution by means of flow cytometry. Journal of Cell Science. 131 (10), 214346 (2018).
  6. Michel, A. M., Baranov, P. V. Ribosome profiling: a Hi-Def monitor for protein synthesis at the genome-wide scale. Wiley Interdisciplinary Reviews: RNA. 4 (5), 473-490 (2013).
  7. Woolstenhulme, C. J., Guydosh, N. R., Green, R., Buskirk, A. R. High-precision analysis of translational pausing by ribosome profiling in bacteria lacking EFP. Cell reports. 11 (1), 13-21 (2015).
  8. McGlincy, N. J., Ingolia, N. T. Transcriptome-wide measurement of translation by ribosome profiling. Methods (San Diego, Calif). 126, 112-129 (2017).
  9. Gerashchenko, M. V., Gladyshev, V. N. Translation inhibitors cause abnormalities in ribosome profiling experiments. Nucleic acids research. 42 (17), 134 (2014).
  10. Ingolia, N. T., Lareau, L. F., Weissman, J. S. Ribosome Profiling of Mouse Embryonic Stem Cells Reveals the Complexity of Mammalian Proteomes. Cell. 147 (4), 789-802 (2011).
  11. Ingolia, N. T., Brar, G. A., Rouskin, S., McGeachy, A. M., Weissman, J. S. The ribosome profiling strategy for monitoring translation in vivo by deep sequencing of ribosome-protected mRNA fragments. Nature protocols. 7 (8), 1534-1550 (2012).
  12. Hsu, P. Y., et al. Super-resolution ribosome profiling reveals unannotated translation events in Arabidopsis. Proceedings of the National Academy of Sciences of the United States of America. 113 (45), 7126-7135 (2016).
  13. Blanchard, S. C., Kim, H. D., Gonzalez, R. L., Puglisi, J. D., Chu, S. tRNA dynamics on the ribosome during translation. Proceedings of the National Academy of Sciences of the United States of America. 101 (35), 12893-12898 (2004).
  14. Becker, A. H., Oh, E., Weissman, J. S., Kramer, G., Bukau, B. Selective ribosome profiling as a tool for studying the interaction of chaperones and targeting factors with nascent polypeptide chains and ribosomes. Nature protocols. 8 (11), 2212-2239 (2013).
  15. Li, G. W., Oh, E., Weissman, J. S. The anti-Shine-Dalgarno sequence drives translational pausing and codon choice in bacteria. Nature. 484 (7395), 538-541 (2012).
  16. King, H. A., Gerber, A. P. Translatome profiling: methods for genome-scale analysis of mRNA translation. Briefings in functional genomics. 15 (1), 22-31 (2014).
  17. Mohammad, F., Woolstenhulme, C. J., Green, R., Buskirk, A. R. Clarifying the Translational Pausing Landscape in Bacteria by Ribosome Profiling. Cell reports. 14 (4), 686-694 (2016).
  18. Reid, D. W., Shenolikar, S., Nicchitta, C. V. Simple and inexpensive ribosome profiling analysis of mRNA translation. Methods (San Diego, Calif). 91, 69-74 (2015).
  19. Sanz, E., Yang, L., Su, T., Morris, D. R., McKnight, G. S., Amieux, P. S. Cell-type-specific isolation of ribosome-associated mRNA from complex tissues. Proceedings of the National Academy of Sciences of the United States of America. 106 (33), 13939-13944 (2009).
  20. Petrova, O. E., Garcia-Alcalde, F., Zampaloni, C., Sauer, K. Comparative evaluation of rRNA depletion procedures for the improved analysis of bacterialbiofilm and mixed pathogen culture transcriptomes. Scientific Reports. 7, 41114 (2017).
  21. Oh, E. Selective ribosome profiling reveals the cotranslational chaperone action of trigger factor in vivo. Cell. 147 (6), 1295-1308 (2011).
  22. Chew, G. L. Ribosome profiling reveals resemblance between long non-coding RNAs and 5′ leaders of coding RNAs. Development. 140 (13), 2828-2834 (2013).
  23. Dunn, J. G., Foo, C. K., Belletier, N. G., Gavis, E. R., Weissman, J. S. Ribosome profiling reveals pervasive and regulated stop codon readthrough in Drosophila melanogaster. eLife. 2, 01179 (2013).
  24. Liu, M. J., et al. Translational landscape of photomorphogenic Arabidopsis. The Plant cell. 25 (10), 3699-3710 (2013).
  25. Rooijers, K., Loayza-Puch, F., Nijtmans, L. G., Agami, R. Ribosome profiling reveals features of normal and disease-associated mitochondrial translation. Nature Communications. 4, 2886 (2013).
  26. Danielle, L., et al. Fidelity of translation initiation is required for coordinated respiratory complex assembly. Science Advances. 5 (12), 2118 (2019).
  27. Zoschke, R., Watkins, K. P., Barkan, A. A rapid ribosome profiling method elucidates chloroplast ribosome behavior in vivo. The Plant cell. 25 (6), 2265-2275 (2013).
  28. Chotewutmontri, P., Barkan, A. Dynamics of chloroplast translation during chloroplast differentiation in maize. PLoS genetics. 12 (7), 1006106 (2016).
  29. Lee, S., Liu, B., Huang, S. X., Shen, B., Qian, S. B. Global mapping of translation initiation sites in mammalian cells at single-nucleotide resolution. Proceedings of the National Academy of Sciences of the United States of America. 109 (37), 2424-2432 (2012).
  30. Martinez, T. F., Chu, Q., Donaldson, C., Tan, D., Shokhirev, M. N., Saghatelian, A. Accurate annotation of human protein-coding small open reading frames. Nature chemical biology. 16 (4), 458-468 (2020).
  31. Simsek, D., et al. The Mammalian Ribo-interactome Reveals Ribosome Functional Diversity and Heterogeneity. Cell. 169 (6), 1051-1065 (2017).
  32. Fritsch, C., et al. Genome-wide search for novel human uORFs and N-terminal protein extensions using ribosomal footprinting. Genome research. 22 (11), 2208-2218 (2012).
  33. Lareau, L. F., Hite, D. H., Hogan, G. J., Brown, P. O. Distinct stages of the translation elongation cycle revealed by sequencing ribosome-protected mRNA fragments. eLife. 3, 01257 (2014).
  34. Meydan, S., et al. Retapamulin-Assisted Ribosome Profiling Reveals the Alternative Bacterial Proteome. Molecular Cell. 74 (3), 481-493 (2019).
  35. Wilson, D. N. The A-Z of bacterial translation inhibitors. Critical Reviews in Biochemistry and Molecular Biology. 44, 393-433 (2009).
  36. Nakahigashi, K. Comprehensive identification of translation start sites by tetracycline-inhibited ribosome profiling. DNA research : an international journal for rapid publication of reports on genes and genomes. 23 (3), 193-201 (2016).
  37. Weaver, J., Mohammad, F., Buskirk, A. R., Storz, G. Identifying Small Proteins by Ribosome Profiling with Stalled Initiation Complexes. mBio. 10 (2), 02819 (2019).
  38. Nakahigashi, K. Effect of codon adaptation on codon-level and gene-level translation efficiency in vivo. BMC genomics. 15 (1), 1115 (2014).
  39. Li, G. W., Burkhardt, D., Gross, C., Weissman, J. S. Quantifying absolute protein synthesis rates reveals principles underlying allocation of cellular resources. Cell. 157 (3), 624-635 (2014).
  40. Michel, A. M., et al. Observation of dually decoded regions of the human genome using ribosome profiling data. Genome research. 22 (11), 2219-2229 (2012).
  41. Guydosh, N. R., Green, R. Dom34 rescues ribosomes in 3′ untranslated regions. Cell. 156 (5), 950-962 (2014).
  42. Young, D. J., Guydosh, N. R., Zhang, F., Hinnebusch, A. G., Green, R. Rli1/ABCE1 Recycles Terminating Ribosomes and Controls Translation Reinitiation in 3’UTRs In Vivo. Cell. 162 (4), 872-884 (2015).
  43. Ingolia, N. T. Ribosome profiling: new views of translation, from single codons to genome scale. Nature Reviews Genetics. 15 (3), 205-213 (2014).
  44. Dobin, A., et al. ultrafast universal RNA seq aligner. Bioinformatics. 29 (1), 15-21 (2013).
  45. . . Babraham Bioinformatics. Fast QC. , (2021).
  46. Gerashchenko, M. V., Gladyshev, V. N. Ribonuclease selection for ribosome profiling. Nucleic acids research. 45 (2), 6 (2017).
  47. Verbruggen, S., Menschaert, G. mQC: A post-mapping data exploration tool for ribosome profiling. Computer Methods and Programs in Biomedicine. 181, 104806 (2019).
  48. Cui, H., Hu, H., Zeng, J., Chen, T. DeepShape: estimating isoform-level ribosome abundance and distribution with Ribo-seq data. BMC bioinformatics. 20, 678 (2019).
  49. Choi, J. . RiboToolkit: an integrated platform for analysis and annotation of ribosome profiling data to decode RNA translation at codon resolution. , (2021).
  50. Choi, J. Dynamics of the context-specific translation arrest by chloramphenicol and linezolid. Nature chemical biology. 16 (3), 310-317 (2020).
  51. Tompson, J., O’Connor, M., Mills, J. A., Dahlberg, A. E. The protein synthesis inhibitors, oxazolidinones and chloramphenicol, cause extensive translational inaccuracy in vivo. Journal of Molecular Biology. 322 (2), 273-279 (2002).
  52. Brar, G. A., Weissman, J. S. Ribosome profiling reveals the what, when, where and how of protein synthesis. Nature reviews. Molecular cell biology. 16 (11), 651-664 (2015).
  53. Glaub, A., Huptas, C., Neuhaus, K., Ardern, Z. Recommendations for bacterial ribosome profiling experiments based on bioinformatic evaluation of published data. Journal of Biological Chemistry. 295, 8999-9011 (2020).
  54. Sorour, M. H., Hani, H. A., Shaalan, H. F., El-Sayed, M. M. H. Experimental screening of some chelating agents for calcium and magnesium removal from saline solutions. Desalination and Water Treatment. 57 (48-49), 22799-22808 (2015).
  55. Mohammad, F., Green, R., Buskirk, A. R. A systematically-revised ribosome profiling method for bacteria reveals pauses at single-codon resolution. Elife. 8, 42591 (2019).
  56. O’Connor, P. B., Li, G. W., Weissman, J. S., Atkins, J. F., Baranov, P. V. rRNA:mRNA pairing alters the length and the symmetry of mRNA-protected fragments in ribosome profiling experiments. Bioinformatics. 29 (12), 1488-1491 (2013).
  57. Protocol for RNA Clean & Concentrator -25. Zymo Research Available from: https://files.zymoresearch.com/protocols/_r1017_r1018_rna_clean_concentrator-25.pdf (2021)
  58. Protocol for use with NEBNext Small RNA Library Prep Set for Illumina (E7300, E7580, E7560, E7330). New England Biolabs Available from: https://www.international.neb.com/protocols/2018/03/27/protocol-for-use-with-nebnext-small-rna-library-prep-set-for-illumina-e7300-e7580-e7560-e7330 (2021)

Play Video

Cite This Article
Kopik, N., Chrobak, O., Latoch, P., Kovalenko, M., Starosta, A. L. RIBO-seq in Bacteria: a Sample Collection and Library Preparation Protocol for NGS Sequencing. J. Vis. Exp. (174), e62544, doi:10.3791/62544 (2021).

View Video