This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Removal of host reads from shotgun metagenomics data.

The presence of host sequence in shotgun metagenome data is considered as contamination and it is removed by mapping with a reference genome in different tools like bowtie2. Can anyone kindly through some light on how to create the reference genome when the whole genome sequence data is not available for the host organism? I am currently working on the gut microbiota of a lepidopteran species and its whole genome is not sequenced yet.

Any response is much appreciated.

metagenomics shotgun
• 3,635 views • link
updated by GenoMax GenoMax • written by dharitrisaikia162 • dharitrisaikia162

Thank you for your reply GenoMax , @colindaven and @Mensur Dlakic. I will follow your suggestions and let you know.

• 0 views • link
written by dharitrisaikia162 • dharitrisaikia162

3 answers

You should use a tool like bbsplit.sh to bin the reads. (LINK).

gut microbiota of a lepidopteran organism and its whole genome is not sequenced yet.

That is going to limit what you can identify and separate for host genome.

• 0 views • link
written by GenoMax GenoMax

You should just be able to use cat *.fa > draft_genome.fa to combine the contigs to a draft genome. Make sure you check it's usefulness using samtools faidx x.fa . Also check whole genome size is as expected.

Then as GenoMax says you can use bbsplit.sh to proceed.

• 0 views • link
written by colindaven colindaven

You can remove the host contigs after the assembly, and that way you will get to assemble some of the insect genome. It should be plenty different from whatever is in its gut, and I expect it will easily separate when you bin the contigs based on tetranucleotide frequencies. The only downside is that it may make the assembly difficult depending on the proportion of host / metagenomic reads.

• 0 views • link
written by Mensur Dlakic Mensur Dlakic

Log in to answer this question.

More posts like this