Depletion of Ribosomal RNA for Mosquito Gut Metagenomic RNA-seq

Summary

A ribosomal RNA (rRNA) depletion protocol was developed to enrich messenger RNA (mRNA) for RNA-seq of the mosquito gut metatranscriptome. Sample specific rRNA probes, which were used to remove rRNA via subtraction, were created from the mosquito and its gut microbes. Performance of the protocol can result in the removal of approximately 90-99% of rRNA.

Abstract

The mosquito gut accommodates dynamic microbial communities across different stages of the insect's life cycle. Characterization of the genetic capacity and functionality of the gut community will provide insight into the effects of gut microbiota on mosquito life traits. Metagenomic RNA-Seq has become an important tool to analyze transcriptomes from various microbes present in a microbial community. Messenger RNA usually comprises only 1-3% of total RNA, while rRNA constitutes approximately 90%. It is challenging to enrich messenger RNA from a metagenomic microbial RNA sample because most prokaryotic mRNA species lack stable poly(A) tails. This prevents oligo d(T) mediated mRNA isolation. Here, we describe a protocol that employs sample derived rRNA capture probes to remove rRNA from a metagenomic total RNA sample. To begin, both mosquito and microbial small and large subunit rRNA fragments are amplified from a metagenomic community DNA sample. Then, the community specific biotinylated antisense ribosomal RNA probes are synthesized in vitro using T7 RNA polymerase. The biotinylated rRNA probes are hybridized to the total RNA. The hybrids are captured by streptavidin-coated beads and removed from the total RNA. This subtraction-based protocol efficiently removes both mosquito and microbial rRNA from the total RNA sample. The mRNA enriched sample is further processed for RNA amplification and RNA-Seq.

Introduction

Next-generation sequencing technology has greatly advanced metagenomics study by allowing to assess the taxonomic composition and genetic functionality of a microbial assemblage. RNA-Sequencing (RNA-Seq)¹ can bypass culture-based methods to investigate microbial metatranscriptomes in different contexts ^2-5. A major obstacle to microbial RNA-seq is the difficulty in enriching mRNA, as the prokaryotic mRNA species are not stably polyadenylated. Therefore, oligo d(T) mediated messenger enrichment is not applicable. Removal of abundant rRNA is an alternative approach to enrich mRNA. Commercial rRNA depletion kits, such as Microexpress Bacterial mRNA Enrichment kit (Ambion), RiboMinus Transcriptome Isolation Kit (Bacteria) (Life Technologies), and mRNA-ONLY Prokaryotic mRNA Isolation kit (Epicentre) that preferentially degrades rRNA with an exonuclease, have been used for removing rRNA ^6-8. However, the capture probes in Microexpress or RiboMinus are good for removing known rRNA from typical Gram-positive and Gram-negative bacteria (see manufacturers' specifications), but less compatible with rRNA from unknown microbes. Consequently, the removal may be less efficient for metagenomic samples^8-10. Besides, the mRNA abundance fidelity was questionable when the exonuclease treatment was applied ¹¹. Overall, subtraction-based rRNA depletion was less biased and more effective in mRNA enrichment in metagenomic settings ^10-13.

The mosquito gut accommodates a dynamic microbial community¹⁴. We are interested in characterizing the function of the mosquito gut microbiome by using RNA-Seq. In an RNA sample isolated from the mosquito guts, both mosquito and microbial RNA are present. Here, we describe a modified protocol to use community specific rRNA probes to efficiently deplete mosquito and microbial rRNA by subtractive hybridization. The resultant mRNA enriched samples are appropriate for RNA-Seq. The overall workflow is depicted in Figure 1.

Protocol

Procedure

1. Mosquito Rearing

1. Rear mosquito Anopheles gambiae G3 strain in an insectary at 27.5 °C with 80% humidity and 12:12 hr cycle of light/dark.
2. Feed larvae with ground cat food with brewer's yeast at ratio 1:1.
3. Feed adult mosquitoes on mice blood at day 3 post-emergence for egg production.

2. Mosquito Gut Dissection

1. Autoclave dissection tools (slides and forceps).
2. Collect 50 mosquitoes using an aspirator and place them on CO₂ flow bed (Ultimate Flypad). Rinse a mosquito specimen sequentially in three Petri dishes containing 70% ethanol to clean the mosquito body surface.
3. Place the specimen on a glass slide under a stereomicroscopy. Remove the gut.
4. Collect 50 guts per condition for metagenomic DNA isolation. Collect 50 guts per condition for metagenomic RNA isolation.

3. Metagenomic DNA Isolation

Note: Metagenomic DNA is isolated using the Meta-G-Nome DNA Isolation Kit (Epicentre Biotechnologies) with the modification described below.

1. Place 50 mosquito guts in 300 μl TE. Homogenize them for 1 min using Bio-Gen PRO200 homogenizer (PRO Scientific Inc, USA) at speed 2,000 rpm on ice.
2. Add 2 μl of Ready-Lyse Lysozyme Solution and 1 μl of RNase A to the cell suspension. Mix by vortexing and incubate at 37 °C for 30 min.
3. Add 300 μl of Meta-Lysis Solution (2X) and 1 μl of Proteinase K to the tube. Mix by vortexing. Briefly pulse-centrifuge the tube to ensure that all of the solution is in the bottom of the tube, and incubate at 65 °C for 15 min, cool to room temperature, then place on ice for 3-5 min.
4. Add 350 μl of MPC Protein Precipitation Reagent to the tube and mix by vortexing vigorously for 10 sec.
5. Pellet the debris by centrifugation for 10 min at 12,000 x g at 4 °C.
6. Transfer the supernatant to a clean 2-ml tube, and discard the pellet.
7. Add 570 μl of isopropanol to the supernatant. Mix by inverting the tube several times.
8. Pellet the DNA by centrifugation for 10 min at 12,000 x g at 4 °C. Remove the isopropanol without dislodging the DNA pellet.
9. Add 500 μl of 75% ethanol to wash the pellet. Centrifuge for 5 min at 12,000 x g at 4 °C. Remove the ethanol without disturbing the DNA pellet. Air-dry the pellet at room temperature for 2 min. Note: Do not over-dry the DNA pellet.
10. Resuspend the DNA pellet in 50 μl of TE Buffer.
11. Validate the size of the isolated DNA by comparison to the Fosmid Control DNA (40 kb; 100 ng/μl) provided in the kit, via gel electrophoresis on a 1% agarose gel.

4. Total RNA Isolation

Note: Clean the working area with RNaseZap to minimize RNase contamination. Metagenomic RNA was isolated from mosquito gut specimens using the TriPure Isolation Reagent (Roche) with a minor modification.

1. Put 50 mosquito guts in a 2-ml tube with 500 μl TriPure Isolation Reagent. Homogenize for 1 min at speed 2,000 rpm on ice using a homogenizer. Let sit at RT for 5 min to allow dissociation of nucleoprotein complexes.
2. Add 100 μl chloroform and vortex for 12 sec, and let sit at RT for 2 min.
3. Centrifuge at 10,000 x g for 10 min. Carefully take out the tube without disturbing separated phases.
4. Transfer aqueous phase to a new tube, add 250 μl isopropanol, mix well, and put at -20 °C for 20 min.
5. Centrifuge at 12,000 x g for 10 min at 4 °C to precipitate RNA.
6. Discard the supernatant without dislodging pellet. Wash the pellet with 750 μl 75% ethanol.
7. Pipette off the supernatant without disturbing the pellet. Air dry the tube for 2 min.
8. Resuspend the RNA with 30 μl nuclease free water.
9. Treat the total RNA with DNase at 37 °C for 20 min, and then inactivate the DNase by incubating at 75 °C for 15 min. Precipitate RNA with ethanol and resuspend RNA in 30 μl nuclease free water.
10. Determine the quantity of total RNA using a NanoDrop.

5. Ribosomal RNA Depletion

Note: The protocol was developed based on previously described methods^12,13,15.

1. PCR amplification of ribosomal RNA gene fragments
   This step creates sample-specific rRNA pools, which will be used as templates for in vitro transcription to produce anti-sense ribosomal RNA (arRNA) probes that are complimentary to rRNA in the total RNA sample.
   1. Design and synthesize primer sets for rRNA fragment PCR (Table 1).
   2. Run PCR in 50 μl reaction using Taq DNA polymerase (Qiagen) with 50 ng DNA template, 1 x PCR buffer, 1.5mM Mg₂Cl, 0.2 μM primers with 35 cycles of denaturing at 94 °C for 10 sec, annealing 15 sec at 40 °C for bacterial 23S amplicon, and 50 °C for the other amplicons, and extending at 72 °C for 1 min (final extension at 72 °C for 5 min).
   3. View PCR products on 1% agarose gel.
   4. Purify PCR products using the QIAquick PCR purification kit with elution in 50 μl elution buffer. Quantify the concentration using NanoDrop.
2. In vitro transcription of biotin-labeled anti-sense rRNA probes (arRNA). Note: In this step, sample specific probes from various rRNA amplicons are synthesized using the MEGAscript T7 kit.
   1. Set up a 20 μl reaction for each sample specific rRNA amplicon by mixing the below ingredients (Table 2). Incubate overnight at 37 °C. Note: usually the reaction yields >50 μg RNA.
   2. Add 1 μl DNase I (included in the kit) to the reaction and incubate at 37 °C for 20 min to remove the DNA template.
   3. Add 100 μl 100% ethanol to the reaction, centrifuge at 12,000 x g for 10 min at 4 °C to precipitate synthesized arRNA. Discard supernatant and wash the pellet with 750 μl 75% ethanol. Resuspend RNA probes in 50 μl nuclease free water.
   4. Quantify RNA concentration using a NanoDrop. Store probes at -80 °C.
3. rRNA subtraction with biotinylated arRNA. Note: This step employs the biotinylated arRNA to hybridize to rRNA in the total RNA sample. The rRNA-arRNA hybrids are captured by streptavidin-coated magnetic beads. The procedure is based on the protocol described in reference 13 with minor modification.
   1. Reagents
      1. Prepare 0.1 N NaOH, 20X sodium chloride-citrate (SSC) and 1 X SSC with 20% formamide.
   2. Hybridization
      1. Mix 1 μg total RNA with bacterial 16S and 23S probes (0.75 μg each), mosquito 18S and 28S probes (0.75 μg each), 1 μl RNase inhibitor, 2.5 μl 20X SSC buffer, 10 μl 100% formamide in a 200 μl PCR tube, add nuclease-free water to 50 μl.
      2. Place the reaction tube in a thermal cycler and incubate at 70 °C for 5 minutes to open secondary structure and ramp down to 25 °C using 5 °C increments for 1 min each °C to allow the hybridization of rRNA and arRNA probes.
   3. Removal of rRNA
      1. Transfer 300 μl streptavidin coated beads into a 1.7 ml tube. Place the tube on a magnetic separation rack to immobilize the beads. Remove supernatant. Resuspend the beads in 300 μl of 0.1N NaOH and mix well. Place the tube back to the magnetic rack and pipette off the supernatant. Resuspend the beads in 300 μl of 1X SSC, mix well, immobilize the beads and pipette off the supernatant. Repeat the wash with 300 μl of 1X SSC one more time.
      2. Make 3 aliquots of beads, 100 μl each, in 1.7 ml tubes. Place the tube on the magnetic rack, and remove the supernatant by pipetting. Note: The aliquots of beads are used for 3 rounds of capturing biotinylated rRNA-arRNA hybrids.
      3. Add 50 μl 1X SSC with 20% formamide into the 50 μl hybridization reaction. Transfer the reaction into the first bead aliquot tube, and mix well. Incubate at RT for 10 min to allow the biotinylated probe-rRNA hybrids to be captured by streptavidin beads. Flick the tube every 2 minutes to mix during the incubation.
      4. Place the tube on the magnetic rack to immobilize the beads, transfer the supernatant to the second aliquot of beads, mix well and incubate as above. Conduct the third capture by transferring supernatant into the third tube of beads. Mix and incubate.
      5. Transfer non-rRNA supernatant into a new tube and conduct purification using RNeasy MinElute kit to remove formamide following the manufacturer's manual.
      6. Check the rRNA depletion efficiency by using an Agilent RNA 6000 Pico chip kit on a Bioanalyzer (Figure 4).

Note: The rRNA depleted RNA samples are subject to mRNA amplification if necessary⁸.

Results

The protocol includes three sections: (1) preparation of metagenomic DNA templates for rRNA PCR; (2) creation of sample specific capture rRNA probes; (3) depletion of rRNA from total RNA by subtractive hybridization. Isolation of high quality metagenomic DNA and RNA is essential to the entire process. The modified Meta-G-Nome DNA isolation protocol yields high-quality metagenomic DNA from mosquito guts, as shown in Figure 2. It can be challenging to isolate high-quality total RNA from mosquito guts. This...

Discussion

A complex microbial community resides in the mosquito gut ecosystem ^14,16,17. Metatranscriptomic sequencing (RNA-seq) can reveal context dependent functional information by interrogating the entire microbial transcriptome^4,18. Technically, oligo-d(T) mediated enrichment of prokaryotic mRNA is not applicable due to the absence of stable poly (A) tails of the messengers. Alternatively, rRNA depletion has been used for mRNA enrichment. Here, we developed a subtraction-based protocol to deplete rRNA usi...

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
Reagent/Material
Meta-G-Nome DNA isolation kit	Epicentre	MGN0910	Metagenomic DNA isolation
TriPure	Roche	11667165001	Mdetagenomic RNA isolation
MEGAscript T7 kit	Ambion	AM1334	In vitro synthesis of RNA probes
RNaseZap	Ambion	AM9780	RNase free working area
Biotin-16-UTP	Roche	11388908910	In vitro synthesis of RNA
Biotin-11-CTP	Roche	4739205001	In vitro synthesis of RNA
Streptavidin magnetic beads	NEB	S1420S	Capture of rRNA hybrids
Magnetic separation rack	NEB	S1506S	Capture of rRNA hybrids
RNeasy mini kit	QIAGEN	74104	Purification of subtracted RNA
RNase-Free DNase Set	QIAGEN	79254	Removal DNA contamination
Agilent RNA 6000 Pico Kit	Agilent Technologies Inc.	5067-1513	Electropherogram of RNA
Equipment
Bio-Gen PRO200 Homogenizer	PRO Scientific	01-01200	Mosquito gut tissue homogenization
NanoDrop 1000 Spectrophotometer	Thermo Scientific		DNA & RNA quantitation
2100 Bioanalyzer	Agilent Technologies Inc.	G2940CA	Electropherogram of RNA