SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
separate the unique and redundant reads from the raw data vishwesh Bioinformatics 2 04-11-2016 06:53 AM
Bfast alignement with paired end reads in separate files david.tamborero Bioinformatics 2 11-29-2011 08:49 AM
How to separate coverage of forward and reverse reads on same axis? Kennels Bioinformatics 6 05-04-2011 01:12 AM
PubMed: De novo assembly of a 40 Mb eukaryotic genome from short sequence reads: Sord Newsbot! Literature Watch 0 04-14-2010 03:01 AM
De novo Assembly of a 40 Mb Eukaryotic Genome from Short Sequence Reads Michael.James.Clark Literature Watch 1 04-09-2010 01:16 PM

Reply
 
Thread Tools
Old 04-24-2017, 09:50 AM   #1
SaiPrabhas
Junior Member
 
Location: Lubbock,TX.

Join Date: Apr 2017
Posts: 1
Default How to separate eukaryotic sequence reads from my metagenome data?

Hello!

I have Illumina Hiseq data from leafhopper Empoasca fabae. It has two symbionts, Sulcia muelleri, and Nasuia deltocephalinicola. In my project, we sequenced head and bacteriome regions of 4 different leafhopper specimens. We want to know if they are any differences in bacterial populations between these specimens and also want to see if they are any new bacteria present in the head samples of leafhopper. As the genome of the leafhopper is not sequenced yet, I cannot map the reads to a reference genome and filter out the eukaryotic reads so that I can analyze the bacterial population. How can I solve this problem? Thank you very much.

Regards,
Sai
SaiPrabhas is offline   Reply With Quote
Old 04-24-2017, 10:04 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,083
Default

Want to make you aware of BBSplit from BBMap.

Following is not completely thought through so take with a grain of salt.

There seem to be too many unknowns in this design to be able to use BBsplit directly. You could bbsplit all samples against each other to see if you can get as much of the leafhopper out as you can but it is going to take out some of the common bacteria as well.

Is there a related genome (to leafhopper) available? That would be the second best bet.

Last edited by GenoMax; 04-24-2017 at 10:08 AM.
GenoMax is offline   Reply With Quote
Old 04-25-2017, 06:26 AM   #3
WhatsOEver
Senior Member
 
Location: Germany

Join Date: Apr 2012
Posts: 215
Default

1) bacteria often have a higher GC-content than eukaryots, so this might be a starting point
2) You could filter your dataset for the known symbionts (which genomes are available, at least acc to a 5s google research) and try an assembly with the rest (depending on what this Hiseq data is: WGS/RNA-Seq?)
3) If you could get a collection of genomes of known insect symbionts, try merging all these and map against this set
4) Contact the people belonging to the abstract "Hybrid De Novo Assembly of the Leafhopper Genome with SMRT Sequencing Method" I also found via google
WhatsOEver 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 07:21 AM.


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