Biostar Beta. Not for public use.
visualize CNV data based on WGS
1
Entering edit mode
13 months ago
Bogdan • 780
Palo Alto, CA, USA

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 • 1.4k views
ADD COMMENTlink
1
Entering edit mode

Did you tried using terminal genome viewer like https://github.com/dariober/ASCIIGenome maybe it will help.

ADD REPLYlink
0
Entering edit mode

Thank you all for all your comments and suggestions ;)

ADD REPLYlink
2
Entering edit mode
12 months ago
France/Nantes/Institut du Thorax - INSE…

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

it works for my needs, e.g:

ADD COMMENTlink
1
Entering edit mode
12 months ago
Leandro Lima • 920
San Francisco, CA

Hi, Bogdan.

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

Leandro

ADD COMMENTlink
0
Entering edit mode

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

ADD REPLYlink
1
Entering edit mode
13 months ago
Eric T. ♦ 2.4k
San Francisco, CA

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.

ADD COMMENTlink
0
Entering edit mode

Dear Eric, thank you very much for your suggestions !

ADD REPLYlink
0
Entering edit mode
11 months ago
Nandini • 810
London

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

ADD COMMENTlink
0
Entering edit mode
12 months ago
WCIP | Glasgow | UK

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.

ADD COMMENTlink
0
Entering edit mode

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

ADD REPLYlink
0
Entering edit mode

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).

ADD REPLYlink
0
Entering edit mode

This is what I mean, My bad :)

ADD REPLYlink
0
Entering edit mode

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

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1