This is a test version of Biostars. For the public version, visit https://www.biostars.org.
visualize CNV data based on WGS

Dear all,

we have been calling CNV (copy-number variations) in cancer genomes based on whole-genome sequencing data (WGS). Please would you advise -- what would be the best way to visualize WGS in order inspect visually the CNV predictions ?
The BAM files of germline and tumor samples are big (> 60-80 GB).

thanks a lot !

-- bogdan

cnv
• 7,338 views • link
updated by clealk • clealk • written by Bogdan Bogdan

Thank you all for all your comments and suggestions ;)

• 0 views • link
written by Bogdan Bogdan

I would like to suggest to keep this gender neutral.

Please use ADD COMMENT or ADD REPLY to answer to previous reactions, as such this thread remains logically structured and easy to follow. I have now moved your post but as you can see it's not optimal. Adding an answer should only be used for providing a solution to the question asked.

If an answer was helpful you should upvote it, if the answer resolved your question you should mark it as accepted. Upvote|Bookmark|Accept

• 0 views • link
written by WouterDeCoster WouterDeCoster

thank you Wouter for your recommendations.

• 0 views • link
written by Bogdan Bogdan

6 answers

I've written : http://lindenb.github.io/jvarkit/LowResBam2Raster.html , you're just limited by the memory...

it works for my needs, e.g:

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Hi, Bogdan.

Have you tried IGV? You can load the BAMs and the BED files to check.

Leandro

• 0 views • link
written by Leandro Lima Leandro Lima

yes, however we encountered a bit of problem with IGV : when zooming-out, the BAM tracks were not visible anymore in IGV.

• 0 views • link
written by Bogdan Bogdan

You can generate a multi-page PDF of scatter plots focused on individual genomic regions (genes, chromosomes, etc.) with cnvkit.py scatter --range-list.

At UCSF you might be able to get access to the group license for Biodiscovery Nexus Copy Number, which is commonly used by pathologists for reviewing array CGH data. You can export CNVkit .cnr files to that program's format with the cnvkit.py export nexus-basic command, or include SNP allele frequencies from a VCF file with export nexus-ogt.

• 0 views • link
written by Eric T. Eric T.

Dear Eric, thank you very much for your suggestions !

• 0 views • link
written by Bogdan Bogdan

gw was designed for this purpose https://github.com/kcleal/gw

You can dynamically visualise whole chromosomes without memory problems

• 0 views • link
written by clealk • clealk

You can try Alamut We use it and it works quite well for large BAM files.

• 0 views • link
written by NB NB

Thanks to @Medhat for "promoting" my tool. However, also ASCIIGenome and any interactive tool working directly with BAM files will be slow when handling very large regions like CNVs (I have to check out Pierre's solution though!). In my opinion, best option is to convert BAM to bigWig or tdf and work with these instead. Have a look at deepTools (for BAM -> bigWig) and igvtools for BAM to tdf.

• 0 views • link
written by dariober dariober

It would be great If you could implement method to go directly to specific region in genome.

• 0 views • link
written by Medhat Medhat

Hi- I think this is already available with the -r/--region option. To go directly to chr:start-end:

ASCIIGenome -r chr7:1000-2000 aln.bam

Is this what you mean? (By the way, ASCIIGenome, like IGV, will intentionally not show BAM tracks for very large regions).

• 0 views • link
written by dariober dariober

This is what I mean, My bad :)

• 0 views • link
written by Medhat Medhat

No Worries- feel free to post questions, suggestions, bugs, whatever!

• 0 views • link
written by dariober dariober

Log in to answer this question.

More posts like this