This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How do you run GISTIC2 with ExomeCNV outputs?

I am using ExomeCNV to call copy numbers of some whole exome sequencing data (tumor and normal). I'd like to use GISTIC2.0 (standalone) to identify recurrent CNVs. The outputs of ExomeCNV are as follow:

  • chr
  • probe_start
  • probe_end
  • coverage
  • targeted.base
  • sequenced.base
  • copy.number
  • logR
  • ratio
  • spec
  • sens
  • average.coverage

What is the recommended way to use the ExomeCNV output to run GISTIC2.0?
More specifically, I'd like to know the answers to the following questions:
1) What should I use as Num markers (number of markers in segment)?
2) What should I use as Seg.CN (log2() -1 of copy number)?

For my question #2 it seems that I can calculate the Seg.CN value in two ways based on the GISTIC documentation (ftp://ftp.broadinstitute.org/pub/GISTIC2.0/GISTICDocumentation_standalone.htm):

Approach #1. log2(2 * ratio) - 1
Approach #2. log2(copy.number) - 1

I think approach #1 is more accurate as the raw ratio is reflected in the input data. Approach #1 is based on the post on using Sequenza outputs for GISTIC2.0 (http://crazyhottommy.blogspot.kr/2017/11/run-gistic2-with-sequenza-segmentation.html). After all, in ExomeCNV outputs there are some deletions (copy.number=1) when the ratio is 0.945 and I want to be more conservative and filter these out in GISTIC2.0.

Thanks.

exomecnv gistic cnv copy number

0 answers

No answers yet.

Log in to answer this question.