Question: Detection of viral nucleic acids in metagenomic aDNA samples
1
gravatar for stu111392
4.4 years ago by
stu11139230
Germany/Kiel
stu11139230 wrote:

Hello there,

at my Institute we are now aiming at detecting viral nucleic acids in metagenomic aDNA samples. As a guy who is totaly new to all the bioinformatics stuff at the moment I'am just looking for programs who seem to be usefull. The main Idea of the identification is to throw all the reads form the sequencer in an multiple alignment software like diamond and blast them against a db. Maybe we are going to do some de-novo assembly befor with meta velvet. Has anyone of you ever done this or is working with metagenomic samples on a familiar task? What programs are you using and is there a better approach?

With kind regards,

Julian

alignment blast sequence assembly • 1.2k views
ADD COMMENTlink modified 4.4 years ago by Mikael Huss4.6k • written 4.4 years ago by stu11139230

This paper may help you Metagenomic Detection of Viruses in Aerosol Samples from Workers in Animal Slaughterhouses

And

Metagenomics for pathogen detection in public health

ADD REPLYlink modified 4.4 years ago • written 4.4 years ago by Medhat8.2k

Thank you fior that. I'll have a look :)

ADD REPLYlink written 4.4 years ago by stu11139230
1
gravatar for Mikael Huss
4.4 years ago by
Mikael Huss4.6k
Stockholm
Mikael Huss4.6k wrote:

BLAST would definitely be too slow for mapping sequencing reads to a DB. Diamond seems much faster (I have only tried it a few times) and may be a viable alternative, although I would perhaps prefer to use it for mapping assembled contigs to a protein database. Depends on how many reads you have.

You could try something like this:

  • Try to remove host DNA by mapping to, e g, the human genome (if that is the host) using bowtie2 or bwa, and discarding the reads that match the host. Here, you can also throw in additional filters such as suspected contaminating bacteria and so on.
  • Map the remaining reads with a short-read aligner (e g bowtie2/bwa) against a (nucleotide) viral database. You could get that from GenBank, for example.
  • It is possible that you could use Diamond at this step as well to map to a (protein) viral reference database.
  • The reads that still remain at this point (which didn't map to host, viral nt database or viral aa database) can be assembled by e g IDBA-UD, MegaHit or something like that.
  • The resulting contigs can be mapped by Diamond to the NR database.
ADD COMMENTlink written 4.4 years ago by Mikael Huss4.6k

That sounds like a good approach. I'll talk about your sugesstions with the others. Thank you for your help.

ADD REPLYlink written 4.4 years ago by stu11139230
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1900 users visited in the last hour