SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Coverage and Read Depth Recommendations by Sequencing Application Genohub Bioinformatics 1 05-12-2015 07:34 PM
How does read length affect sequencing depth for paired end RNA-Seq? bob-loblaw RNA Sequencing 2 02-15-2013 03:02 PM
picard error: Mismatch between read length and quals length writing read shawpa Bioinformatics 0 08-20-2012 05:52 AM
Does the read depth distribution(not read depth) matters yuhao Bioinformatics 0 08-09-2012 05:22 AM
Read depth recommendations dpryan RNA Sequencing 2 09-30-2011 11:15 AM

Reply
 
Thread Tools
Old 10-14-2016, 05:14 PM   #1
TimK
Member
 
Location: Australia

Join Date: Jul 2014
Posts: 13
Default Metatranscriptome recommendations: read length and depth

Hi all,

I have a couple of questions regarding metatranscriptomics of environmental samples.

1. What is the recommended read length?
I'm aware that for transcriptome data sets where a reference genome exists read length is often less important. However, given that metatranscriptomics include assembly of transcripts and subsequent mapping of reads back to the assembly I would assume that, as for metagenomics, metatranscriptome reads should be as long as possible.

2. Overlapping reads or not?
Am I correct that read overlapping, especially long/complete overlapping regions, are only useful for application cases where accuracy should be maximized, e.g. rRNA aplicon sequencing?

3. Sequencing depth:
I guess sequencing depth is always a tricky one and depends not only on the question but also on the environment (high/low diversity/complexity). But is there a better rule of thumb than "as much as possible"?

Thanks,
Tim

Last edited by TimK; 10-14-2016 at 05:16 PM.
TimK is offline   Reply With Quote
Old 10-15-2016, 04:22 AM   #2
Ali May
Member
 
Location: Netherlands

Join Date: Aug 2016
Posts: 12
Default

I have some experience in comparative metatranscriptomics (e.g. disease vs healthy). From my point of view you have two options: HiSeq 2500 50 bp or 150x2 bp. The rule of thumb is that more depth never hurts in environmental sequencing, so if you have access to an instrument with even more yield, go for it. I try to answer your questions below:

1- If you want to do metatranscriptome assembly, needless to say, the read length becomes very important, as it is likely with very short reads to end up having many chimeric assemblies. But there is a trade-off between read length and sequencing depth, so an educated choice here is necessary. Also, depth as you mentioned is very closely linked to the environment of interest. You don't need very deep sequencing for a vaginal sample for instance, but the same doesn't go for oral samples, or soil etc.

2- Although I'm not very sure, I do not think that this is important.

3- In my opinion read length and sequencing depth should be considered together. If you have 40 different bacteria (2.5 -3 MB genomes) in your sample, I would suggest at least 30x. But, to be really honest, I am not sure if anyone can show you what the difference would be between 30x and 10x.

Well I didn't really say much than you already did, but I hope these help. On a different note: one thing to consider is to not do assembly... I've read and heard people saying that it's super important to assemble, because it lets you map many more reads to the transcriptomes. Fair enough... But I've never assembled and it worked just fine for me. especially if you have access to something like KEGG.
Ali May is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 06:01 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO