SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
How to extract assembled transcript sequence from RNA-seq data instead of ref genome? lzu Bioinformatics 10 02-24-2014 01:56 AM
De novo assembly from DNA sequence and RNA seq jjjscuedu De novo discovery 3 11-11-2013 08:28 PM
Gene name assignment/ Merge genome-based assembled transcripts with de novo transcrip htk De novo discovery 0 04-01-2013 03:32 AM
Finding alternative splicing using assembled transcripts sadiexiaoyu RNA Sequencing 0 04-01-2013 02:37 AM
Finding *new* regions of DNA in genome assemblies green tree De novo discovery 5 02-20-2012 03:19 PM

Reply
 
Thread Tools
Old 12-16-2013, 01:54 PM   #1
cwisch88
Member
 
Location: St. Louis

Join Date: Jan 2012
Posts: 15
Default Finding 16s DNA Sequence in de novo assembled genome

From time to time I find I'm working with a bunch of de novo assemblies. I've been interested in doing some comparative genomics and perhaps making some group-specific 16s rRNA primers.

I'm just wondering if there is an accepted way of extracting 16s DNA from de novo assemblies that I'm just not finding.

I've used BLAST and other Local alignment programs using universal primers, but this doesn't strike me as very efficient and I'm sure I'm not the only one trying do something like this. Can anyone point me toward any papers or personal experience doing something like this.

Thanks!
cwisch88 is offline   Reply With Quote
Old 12-16-2013, 02:32 PM   #2
boetsie
Senior Member
 
Location: NL, Leiden

Join Date: Feb 2010
Posts: 245
Default

You can use a tool like barrnap

HTML Code:
http://www.vicbioinformatics.com/software.barrnap.shtml
or RNAmmer
HTML Code:
http://www.cbs.dtu.dk/services/RNAmmer/
to predict rRNA locations on your contig sequences.

Boetsie
boetsie is offline   Reply With Quote
Old 12-17-2013, 09:59 AM   #3
cwisch88
Member
 
Location: St. Louis

Join Date: Jan 2012
Posts: 15
Default Thanks

I will definitely check these out!
cwisch88 is offline   Reply With Quote
Old 12-17-2013, 10:02 AM   #4
rhinoceros
Senior Member
 
Location: sub-surface moon base

Join Date: Apr 2013
Posts: 372
Default

Blast the contigs against a 16S db?
__________________
savetherhino.org
rhinoceros is offline   Reply With Quote
Old 12-17-2013, 10:09 AM   #5
cwisch88
Member
 
Location: St. Louis

Join Date: Jan 2012
Posts: 15
Default

Quote:
Originally Posted by rhinoceros View Post
Blast the contigs against a 16S db?
That is a very straightforward approach and I have done that in past. I was just wondering if there was something that could quickly extract the sequences without reinventing the wheel.

It looks like the barrnap program does a quick annotation and then I can use something like bedtools to extract the sequences using the produced GFF.
cwisch88 is offline   Reply With Quote
Reply

Tags
16s, alignment, de novo assembly, dna, illumina

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 07:55 AM.


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