SEQanswers

Go Back   SEQanswers > Applications Forums > Metagenomics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Metagenomics/Metatrancriptomics Sequencing newBioinfo Core Facilities 11 09-20-2013 11:41 AM
Metagenomics/Metatrancriptomics Sequencing newBioinfo Metagenomics 0 09-13-2013 10:51 AM
16s rRNA gene sequencing(metagenomics) newBioinfo Sample Prep / Library Generation 0 03-26-2012 05:08 PM
16s rRNA gene sequencing(metagenomics) newBioinfo Illumina/Solexa 0 03-26-2012 02:24 PM

Reply
 
Thread Tools
Old 04-22-2016, 05:50 AM   #1
JQL
Member
 
Location: MO, USA

Join Date: Apr 2011
Posts: 76
Default sequencing QC and metagenomics

Hi,

I just started my first metagenomics shot-gun sequencing analysis. I hope I can get some advice here.

1. QC process: I am looking at removing duplicate sequences and trimming low quality bases, according to this HMP's SOP. Is this adequate? What is considered a good QC practice?

2. Data analysis: I am first trying MetaPhlAn by Huttenhower's group. Is there any 'standard' procedure/workflow that this community recommends?

Thanks for your help,
John
JQL is offline   Reply With Quote
Old 04-22-2016, 11:31 AM   #2
fanli
Senior Member
 
Location: California

Join Date: Jul 2014
Posts: 192
Default

Quote:
Originally Posted by JQL View Post
1. QC process: I am looking at removing duplicate sequences and trimming low quality bases, according to this HMP's SOP. Is this adequate? What is considered a good QC practice?
You may want to remove human/host DNA as well. kneaddata from Curtis Huttonhower's group does this, but uses bowtie2's --un-conc output to filter reads, which misses repetitive DNA that doesn't align concordantly. We've taken to manually filtering any read pair that aligns in any fashion to the human reference.

Quote:
Originally Posted by JQL View Post
2. Data analysis: I am first trying MetaPhlAn by Huttenhower's group. Is there any 'standard' procedure/workflow that this community recommends?
There are also k-mer based approaches (e.g. kraken and CLARK) that were recently shown to outperform MetaPhlAn: see here. Generally speaking, I don't think shotgun is as well figured out as 16S, so there is still some exploring to do. You may also be interested in PanPhlAn, ConStrains, etc. that do strain-level profiling.

We've just gotten into shotgun too, so happy to share any advice.
fanli is offline   Reply With Quote
Old 04-22-2016, 12:33 PM   #3
JQL
Member
 
Location: MO, USA

Join Date: Apr 2011
Posts: 76
Default

Thank you very much for your reply. At this point, any thing will help me get started.

have a nice weekend,
John

Quote:
Originally Posted by fanli View Post
You may want to remove human/host DNA as well. kneaddata from Curtis Huttonhower's group does this, but uses bowtie2's --un-conc output to filter reads, which misses repetitive DNA that doesn't align concordantly. We've taken to manually filtering any read pair that aligns in any fashion to the human reference.



There are also k-mer based approaches (e.g. kraken and CLARK) that were recently shown to outperform MetaPhlAn: see here. Generally speaking, I don't think shotgun is as well figured out as 16S, so there is still some exploring to do. You may also be interested in PanPhlAn, ConStrains, etc. that do strain-level profiling.

We've just gotten into shotgun too, so happy to share any advice.
JQL is offline   Reply With Quote
Old 07-20-2016, 07:42 AM   #4
bastianwur
Member
 
Location: Germany/Netherlands

Join Date: Feb 2014
Posts: 98
Default

Other things you'll need to do:
- assembly, best cross-assembly if you have multiple samples
- genomic binning (e.g. maxBin)
- gene prediction (prodigal has a setting for meta, but only bacteria; ncRNAs can be predicted with rnammer, trnascan-SE and rfam)
- function prediction (via KO, InterproScan, PRIAM, dbCAN, etc)
bastianwur 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:36 PM.


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