This is a test version of Biostars. For the public version, visit https://www.biostars.org.
HPV EBV typing from cancer genome data

Hi all,

We designed several probes to enrich the virus, such as hpv, ebv, hbv, from the cancer tissue DNA samples.

Is there any tool that can directly typing these virus from the raw fq data?

Thanks.

Junfeng

hpv ebv ngs
• 1,840 views • link
updated by d-cameron d-cameron • written by J.F.Jiang J.F.Jiang

If your probes are specific you should be able to do that. There is a paper that has a pipeline to do this in RNAseq data.

• 0 views • link
written by GenoMax GenoMax

1 answer

The fastest way to type the viruses directly from the fastq data would be to use a kmer abundance tool such as kraken2.

This approach is fast as it performs neither alignment nor assembly. If you want to know the actual sequence of the virus(es) in the sample, then you'll need to do assembly.

• 0 views • link
written by d-cameron d-cameron

Thanks,

kraken, or other metagenomic tools, apply kmer strategy to search the virus status from the fq data.

Can it do further step, such as determine the HPV alleles (hpv18)?

And I will carefully review the tool, any comments and suggestions will be appreciated.

• 0 views • link
written by J.F.Jiang J.F.Jiang

Yes, kraken2 classifies according to where in the taxonomic tree it can place the reads. If that's just HPV, then it'll classify it with that taxId, if the sequence is unique enough that it can be classified as HPV 18, the it'll allocate as that taxId

• 0 views • link
written by d-cameron d-cameron

Log in to answer this question.

More posts like this