SAMSA: a comprehensive metatranscriptome analysis pipeline.
BACKGROUND
Although metatranscriptomics-the study of diverse microbial population activity based on RNA-seq data-is rapidly growing in popularity, there are limited options for biologists to analyze this type of data. Current approaches for processing metatranscriptomes rely on restricted databases and a dedicated computing cluster, or metagenome-based approaches that have not been fully evaluated for processing metatranscriptomic datasets. We created a new bioinformatics pipeline, designed specifically for metatranscriptome dataset analysis, which runs in conjunction with Metagenome-RAST (MG-RAST) servers. Designed for use by researchers with relatively little bioinformatics experience, SAMSA offers a breakdown of metatranscriptome transcription activity levels by organism or transcript function, and is fully open source. We used this new tool to evaluate best practices for sequencing stool metatranscriptomes.


RESULTS
Working with the MG-RAST annotation server, we constructed the Simple Annotation of Metatranscriptomes by Sequence Analysis (SAMSA) software package, a complete pipeline for the analysis of gut microbiome data. SAMSA can summarize and evaluate raw annotation results, identifying abundant species and significant functional differences between metatranscriptomes. Using pilot data and simulated subsets, we determined experimental requirements for fecal gut metatranscriptomes. Sequences need to be either long reads (longer than 100 bp) or joined paired-end reads. Each sample needs 40-50 million raw sequences, which can be expected to yield the 5-10 million annotated reads necessary for accurate abundance measures. We also demonstrated that ribosomal RNA depletion does not equally deplete ribosomes from all species within a sample, and remaining rRNA sequences should be discarded. Using publicly available metatranscriptome data in which rRNA was not depleted, we were able to demonstrate that overall organism transcriptional activity can be measured using mRNA counts. We were also able to detect significant differences between control and experimental groups in both organism transcriptional activity and specific cellular functions.


CONCLUSIONS
By making this new pipeline publicly available, we have created a powerful new tool for metatranscriptomics research, offering a new method for greater insight into the activity of diverse microbial communities. We further recommend that stool metatranscriptomes be ribodepleted and sequenced in a 100 bp paired end format with a minimum of 40 million reads per sample.
Department of Molecular and Cellular Biology, University of California - Davis, Davis, CA, USA.|Genome Center, University of California - Davis, Davis, CA, USA.
Samuel T
ST
Westreich
Department of Molecular and Cellular Biology, University of California - Davis, Davis, CA, USA.|Genome Center, University of California - Davis, Davis, CA, USA.
Ian
I
Korf
Department of Food Science and Technology, University of California - Davis, Davis, CA, USA.
David A
DA
Mills
Genome Center, University of California - Davis, Davis, CA, USA. dglemay@ucdavis.edu.
Danielle G
DG
Lemay
Best practices; Big data; Metagenome; Metatranscriptome; Microbiome; Pipeline; RNA-seq; Software package
J Immunol. 2015 Nov 15;195(10):5011-24
26432894
Nucleic Acids Res. 2014 Jan;42(Database issue):D553-9
24316578
Proc Natl Acad Sci U S A. 2014 Jun 3;111(22):E2329-38
24843156
PLoS Comput Biol. 2015 Jan 08;11(1):e1004008
25569221
Nucleic Acids Res. 2014 Jan;42(Database issue):D206-14
24293654
BMC Genomics. 2015 Jul 04;16:494
26140923
PLoS One. 2011 Apr 20;6(4):e18493
21533138
Genome Biol. 2012;13(3):R23
22455878
BMC Genomics. 2013 Aug 02;14:530
23915218
Bioinformatics. 2014 Aug 1;30(15):2114-20
24695404
Nat Rev Immunol. 2009 May;9(5):313-23
19343057
Nat Rev Genet. 2012 Mar 13;13(4):260-70
22411464
PLoS One. 2011 Mar 08;6(3):e17447
21408168
J Microbiol Methods. 2008 Oct;75(2):172-6
18582973
Clin Vaccine Immunol. 2006 Nov;13(11):1231-6
16988006
Dig Dis Sci. 1997 Dec;42(12):2577-80
9440640
BMC Bioinformatics. 2008 Sep 19;9:386
18803844
OMICS. 2011 Jul-Aug;15(7-8):419-30
21121704
Bioinformatics. 2011 Nov 1;27(21):2957-63
21903629
ISME J. 2014 Apr;8(4):757-67
24152715
Genome Biol. 2014;15(12):550
25516281