Stool metatranscriptomics: A technical guideline for mRNA stabilisation and isolation
Open Access
- 4 July 2015
- journal article
- research article
- Published by Springer Science and Business Media LLC in BMC Genomics
- Vol. 16 (1), 1-18
- https://doi.org/10.1186/s12864-015-1694-y
Abstract
The complex microbiome of the gut has an enormous impact on human health. Analysis of the transcriptional activity of microorganisms through mRNA sequencing (metatranscriptomics) opens a completely new window into their activity in vivo, but it is highly challenging due to numerous technical and bioinformatical obstacles. Here we present an optimized pipeline for extraction of high quality mRNA from stool samples. Comparison of three commercially available RNA extraction kits with the method of Zoetendal revealed that the Powermicrobiome Kit (MoBio) performed best with respect to RNA yield and purity. Next, the influence of the stabilization reagent during sample storage for up to 15 days was studied. RIN analysis and qRT-PCR of spiked-in and indigenous genes revealed that RNA Later preserved mRNA integrity most efficiently, while samples conserved in RNA Protect showed substantial mRNA decay. Using the optimized pipeline developed here, recovery rates for spiked-in E.coli cells expressing fluorescing proteins were 8.7-9.7 % for SuperfolderGFP and 14.7-17.8 % for mCherry. The mRNA of stabilized stool samples as well as of snap-frozen controls was sequenced with Illumina Hiseq, yielding on average 74 million reads per sample. PCoA analysis, taxonomic classification using Kraken and functional classification using bwa showed that the transcriptomes of samples conserved in RNA Later were unchanged for up to 6 days even at room temperature, while RNA Protect was inefficient for storage durations exceeding 24 h. However, our data indicate that RNA Later introduces a bias which is then maintained throughout storage, while RNA Protect conserved samples are initially more similar to the snap frozen controls. RNA Later conserved samples had a reduced abundance of e.g. Prevotellaceae transcripts and were depleted for e.g. COG category “Carbohydrate transport and metabolism”. Since the overall similarity between all stool transcriptional profiles studied here was >0.92, these differences are unlikely to affect global comparisons, but should be taken into account when rare but critically important members of the stool microbiome are being studied.This publication has 34 references indexed in Scilit:
- General properties of transcriptional time series in Escherichia coliNature Genetics, 2011
- Validation of two ribosomal RNA removal methods for microbial metatranscriptomicsNature Methods, 2010
- Quantitative analysis of a deeply sequenced marine microbial metatranscriptomeThe ISME Journal, 2010
- Extraction of mRNA from SoilApplied and Environmental Microbiology, 2010
- Organismal, genetic, and transcriptional variation in the deeply sequenced gut microbiomes of identical twinsProceedings of the National Academy of Sciences of the United States of America, 2010
- The Sequence Alignment/Map format and SAMtoolsBioinformatics, 2009
- Fast and accurate short read alignment with Burrows–Wheeler transformBioinformatics, 2009
- Characterizing a model human gut microbiota composed of members of its two dominant bacterial phylaProceedings of the National Academy of Sciences of the United States of America, 2009
- The Human Microbiome ProjectNature, 2007
- A new mathematical model for relative quantification in real-time RT-PCRNucleic Acids Research, 2001