SNP6 Copy number analysis (GISTIC2)
Kidney Chromophobe (Primary solid tumor)
15 July 2014  |  analyses__2014_07_15
Maintainer Information
Citation Information
Maintained by Dan DiCara (Broad Institute)
Cite as Broad Institute TCGA Genome Data Analysis Center (2014): SNP6 Copy number analysis (GISTIC2). Broad Institute of MIT and Harvard. doi:10.7908/C1BC3X8V
Overview
Introduction

GISTIC identifies genomic regions that are significantly gained or lost across a set of tumors. The pipeline first filters out normal samples from the segmented copy-number data by inspecting the TCGA barcodes and then executes GISTIC version 2.0.21 (Firehose task version: 127).

Summary

There were 66 tumor samples used in this analysis: 15 significant arm-level results, 2 significant focal amplifications, and 0 significant focal deletions were found.

Results
Focal results

Figure 1.  Genomic positions of amplified regions: the X-axis represents the normalized amplification signals (top) and significance by Q value (bottom). The green line represents the significance cutoff at Q value=0.25.

Table 1.  Get Full Table Amplifications Table - 2 significant amplifications found. Click the link in the last column to view a comprehensive list of candidate genes. If no genes were identified within the peak, the nearest gene appears in brackets.

Cytoband Q value Residual Q value Wide Peak Boundaries # Genes in Wide Peak
8q11.23 0.2211 0.2211 chr8:51697150-53706540 6
15q22.31 0.2211 0.2211 chr15:1-66428823 492
Genes in Wide Peak

This is the comprehensive list of amplified genes in the wide peak for 8q11.23.

Table S1.  Genes in bold are cancer genes as defined by The Sanger Institute: Cancer Gene Census[7].

Genes
ST18
RB1CC1
SNTG1
PCMTD1
PXDNL
FAM150A
Genes in Wide Peak

This is the comprehensive list of amplified genes in the wide peak for 15q22.31.

Table S2.  Genes in bold are cancer genes as defined by The Sanger Institute: Cancer Gene Census[7].

Genes
BUB1B
TCF12
FLJ27352
C15orf55
hsa-mir-4311
hsa-mir-1272
hsa-mir-422a
hsa-mir-190
hsa-mir-2116
hsa-mir-628
hsa-mir-1266
hsa-mir-147b
hsa-mir-1282
hsa-mir-627
hsa-mir-4310
hsa-mir-626
hsa-mir-1233-2
hsa-mir-1233-1
hsa-mir-211
hsa-mir-1268
hsa-mir-3118-6
hsa-mir-3118-4
ACTC1
ADAM10
ANXA2
APBA2
AQP9
B2M
NBEAP1
BNIP2
CA12
CAPN3
CHRM5
CHRNA7
CKMT1B
CYP19A1
DUT
EPB42
FBN1
FGF7
GABPB1
GABRA5
GABRB3
GABRG3
GALK2
GANC
GATM
GCHFR
PDIA3
GTF2A2
HDC
ONECUT1
IPW
ITPKA
IVD
LIPC
LTK
MAP1A
MEIS2
MFAP1
TRPM1
MYO1E
MYO5A
NDN
NEDD4
OAZ2
OCA2
PLCB2
PPIB
MAPK6
RAB27A
RAD51
RORA
RYR3
SCG5
SLC12A1
SNRPN
SNX1
SORD
SPINT1
SRP14
THBS1
TJP1
TP53BP1
TPM1
TYRO3
UBE3A
MKRN3
SLC30A4
PAR5
CILP
EIF3J
JMJD7-PLA2G4B
RAB11A
SNAP23
ALDH1A2
HERC2
HERC1
SNURF
USP8
CCNB2
SLC28A2
SLC24A1
CCPG1
GCNT3
COPS2
TRIP4
TGM5
PIGB
IGDCC3
PPIP5K1
AQR
SECISBP2L
KIAA0101
ARHGAP11A
LCMT2
USP3
SLC12A6
BCL2L10
PDCD7
RASGRP1
SERF2
DENND4A
GNB5
ARPP19
CLPX
SLC27A2
GPR176
CHP
OIP5
BAHD1
FAN1
CEP152
MAPKBP1
GOLGA8A
ZNF609
RTF1
CYFIP1
MGA
DMXL2
VPS39
FAM189A1
AP4E1
CCNDBP1
DAPK2
EID1
C15orf2
C15orf63
TMEM87A
RPAP1
DKFZP434L187
PYGO1
PLDN
GREM1
FOXB1
RPUSD2
TUBGCP4
SCG3
TMOD3
TMOD2
EHD4
DUOX2
MYEF2
RPS27L
NDUFAF1
RSL24D1
NUSAP1
TMEM85
RASL12
SPG21
SPTBN5
PTPLAD1
CTDSPL2
KLF13
RAB8B
DUOX1
CSNK1G1
MAGEL2
DLL4
INO80
FAM63B
RNF111
ZNF280D
TRPM7
VPS13C
PPP1R14D
DPP8
MTMR10
PARP16
ZNF770
FLJ10038
HAUS2
FAM82A2
DNAJC17
MNS1
NOP10
MYO5C
NDNL2
FAM214A
C15orf24
PAK6
DTWD1
CASC5
AVEN
ATP10A
GJD2
STARD9
VPS18
IGDCC4
SQRDL
ZFP106
RFX7
SPATA5L1
CHAC1
NARG2
C15orf29
SLTM
SNX22
ATP8B4
WDR76
TMEM62
SEMA6D
PIF1
SPG11
ELL3
PLEKHO2
POLR2M
C15orf44
NIPA2
APH1B
TLN2
FAM96A
C15orf48
MEGF11
C15orf41
SPPL2A
ZFYVE19
CGNL1
FRMD5
DISP2
CHRFAM7A
ULK4P3
ULK4P1
ARHGAP11B
ATPBD4
C15orf57
C15orf23
BMF
SHF
DUOXA1
SNORD107
LDHAL6B
CHST14
CASC4
LACTB
TUBGCP5
TGM7
CATSPER2
LEO1
MTFMT
OSTBETA
NIPA1
PLA2G4E
TRIM69
PAR1
C15orf43
LOC145663
C2CD4A
FAM81A
GCOM1
LOC145783
LOC145845
TMCO5A
ZSCAN29
TTBK2
CDAN1
STRC
DYX1C1
CSNK1A1P1
OTUD7A
SPRED1
PGBD4
ADAL
EXD1
FSIP1
RHOV
C15orf33
MGC15885
UBR1
PATL2
LPCAT4
PLA2G4F
LRRC57
LYSMD2
WDR72
HMGN2P46
SLC24A5
PRTG
LOC283663
LOC283683
OR4N4
LOC283710
FAM98B
PLA2G4D
HERC2P3
GOLGA6L1
GOLGA8G
GOLGA8IP
FBXL22
SNORD108
SNORD109A
SNORD109B
SNORD115-1
WHAMMP3
POTEB
GLDN
FMN1
SNORD64
PAR4
PAR-SN
RBPMS2
ANKDD1A
LOC348120
MRPL42P5
USP50
TEX9
C15orf52
TNFAIP8L3
C2CD4B
GOLGA8E
OR4M2
OR4N3P
KBTBD13
UBAP1L
SHC4
CTXN2
HERC2P2
C15orf53
C15orf54
DUOXA2
MIR190A
MIR211
NF1P2
HERC2P9
WHAMMP2
GOLGA8B
EIF2AK4
CATSPER2P1
UNC13C
LOC503519
CKMT1A
SERINC4
C15orf62
C15orf56
PHGR1
LOC645212
CHEK2P2
LOC646214
CXADRP2
LOC646278
REREP3
LOC653061
LOC653075
HSP90AB4P
MIR626
MIR627
MIR628
ANP32AP1
SNORD116-19
GOLGA6L6
LOC727924
LOC728758
OIP5-AS1
GOLGA8C
PWRN1
PWRN2
SNORD116-1
SNORD116-2
SNORD116-3
SNORD116-4
SNORD116-5
SNORD116-6
SNORD116-7
SNORD116-8
SNORD116-9
SNORD116-10
SNORD116-11
SNORD116-12
SNORD116-13
SNORD116-14
SNORD116-15
SNORD116-16
SNORD116-17
SNORD116-18
SNORD116-20
SNORD116-21
SNORD116-22
SNORD116-23
SNORD116-24
SNORD116-25
SNORD115-2
SNORD116-26
SNORD116-27
SNORD115-3
SNORD115-4
SNORD115-5
SNORD115-6
SNORD115-7
SNORD115-8
SNORD115-9
SNORD115-10
SNORD115-11
SNORD115-12
SNORD115-13
SNORD115-14
SNORD115-15
SNORD115-16
SNORD115-17
SNORD115-18
SNORD115-19
SNORD115-20
SNORD115-21
SNORD115-22
SNORD115-23
SNORD115-25
SNORD115-26
SNORD115-29
SNORD115-30
SNORD115-31
SNORD115-32
SNORD115-33
SNORD115-34
SNORD115-35
SNORD115-36
SNORD115-37
SNORD115-38
SNORD115-39
SNORD115-40
SNORD115-41
SNORD115-42
SNORD115-43
SNORD115-44
SNORD116-28
SNORD116-29
SNORD115-48
SNORD115-24
SNORD115-27
SNORD115-28
SNORD115-45
SNORD115-47
MIR147B
LOC100128714
LOC100129387
LOC100130855
LOC100131089
ANKRD63
HERC2P7
GOLGA8F
LOC100132724
GOLGA8DP
JMJD7
PLA2G4B
ULK4P2
LOC100288615
LOC100288637
LOC100289656
MIR1233-1
MIR1272
MIR1266
MIR1282
LOC100306975
MIR2116
MIR1233-2
MIR4311
MIR4310
MIR3942
LOC100505648
LOC100507466
SERF2-C15ORF63
DYX1C1-CCPG1
MIR4509-1
MIR4509-2
MIR4508
MIR4510
MIR4716
MIR4713
MIR4511
MIR4509-3
MIR4712
MIR4715
TMCO5B
MYZAP

Figure 2.  Genomic positions of deleted regions: the X-axis represents the normalized deletion signals (top) and significance by Q value (bottom). The green line represents the significance cutoff at Q value=0.25.

Arm-level results

Table 2.  Get Full Table Arm-level significance table - 15 significant results found. The significance cutoff is at Q value=0.25.

Arm # Genes Amp Frequency Amp Z score Amp Q value Del Frequency Del Z score Del Q value
1p 2121 0.08 -1.89 0.995 0.82 8.52 0
1q 1955 0.07 -1.95 0.995 0.80 8.47 0
2p 924 0.10 -1.62 0.995 0.72 8.41 0
2q 1556 0.10 -1.89 0.995 0.72 7.51 1.09e-13
3p 1062 0.12 -2.39 0.995 0.17 -1.68 1
3q 1139 0.14 -2.27 0.995 0.14 -2.27 1
4p 489 0.37 2.54 0.184 0.02 -3.23 1
4q 1049 0.37 1.9 0.312 0.02 -3.54 1
5p 270 0.14 -1.46 0.995 0.17 -0.944 1
5q 1427 0.14 -2.36 0.995 0.17 -1.9 1
6p 1173 0.07 -1.79 0.995 0.78 9.29 0
6q 839 0.07 -1.67 0.995 0.78 9.82 0
7p 641 0.37 2.36 0.184 0.02 -3.32 1
7q 1277 0.37 1.66 0.326 0.02 -3.67 1
8p 580 0.29 0.793 0.57 0.20 -0.669 1
8q 859 0.30 0.746 0.57 0.18 -1.15 1
9p 422 0.16 -1.2 0.995 0.19 -0.7 1
9q 1113 0.16 -1.76 0.995 0.19 -1.29 1
10p 409 0.18 -0.537 0.995 0.78 10.3 0
10q 1268 0.12 -1.47 0.995 0.77 8.74 0
11p 862 0.25 -0.0195 0.995 0.14 -1.93 1
11q 1515 0.25 -0.603 0.995 0.14 -2.4 1
12p 575 0.30 1.06 0.524 0.04 -3.17 1
12q 1447 0.31 0.383 0.78 0.02 -3.97 1
13q 654 0.09 -1.76 0.995 0.67 7.96 3.55e-15
14q 1341 0.33 0.845 0.57 0.04 -3.53 1
15q 1355 0.33 0.832 0.57 0.04 -3.53 1
16p 872 0.34 1.49 0.337 0.09 -2.57 1
16q 702 0.34 1.76 0.312 0.11 -2.11 1
17p 683 0.00 -2.28 0.995 0.76 9.65 0
17q 1592 0.00 -2.58 0.995 0.76 8.27 4.93e-16
18p 143 0.29 1.4 0.36 0.16 -0.918 1
18q 446 0.29 0.926 0.57 0.20 -0.558 1
19p 995 0.29 0.549 0.686 0.02 -3.76 1
19q 1709 0.27 -0.524 0.995 0.06 -3.65 1
20p 355 0.32 1.76 0.312 0.09 -2.29 1
20q 753 0.33 1.54 0.337 0.07 -2.84 1
21q 509 0.13 -1.41 0.995 0.56 6.07 1.94e-09
22q 921 0.33 1.2 0.458 0.17 -1.38 1
Xq 1312 0.22 -0.659 0.995 0.65 6.4 2.51e-10
Methods & Data
Input
Description
  • Segmentation File: The segmentation file contains the segmented data for all the samples identified by GLAD, CBS, or some other segmentation algorithm. (See GLAD file format in the Genepattern file formats documentation.) It is a six column, tab-delimited file with an optional first line identifying the columns. Positions are in base pair units.The column headers are: (1) Sample (sample name), (2) Chromosome (chromosome number), (3) Start Position (segment start position, in bases), (4) End Position (segment end position, in bases), (5) Num markers (number of markers in segment), (6) Seg.CN (log2() -1 of copy number).

  • Markers File: The markers file identifies the marker names and positions of the markers in the original dataset (before segmentation). It is a three column, tab-delimited file with an optional header. The column headers are: (1) Marker Name, (2) Chromosome, (3) Marker Position (in bases).

  • Reference Genome: The reference genome file contains information about the location of genes and cytobands on a given build of the genome. Reference genome files are created in Matlab and are not viewable with a text editor.

  • CNV Files: There are two options for the cnv file. The first option allows CNVs to be identified by marker name. The second option allows the CNVs to be identified by genomic location. Option #1: A two column, tab-delimited file with an optional header row. The marker names given in this file must match the marker names given in the markers file. The CNV identifiers are for user use and can be arbitrary. The column headers are: (1) Marker Name, (2) CNV Identifier. Option #2: A 6 column, tab-delimited file with an optional header row. The 'CNV Identifier' is for user use and can be arbitrary. 'Narrow Region Start' and 'Narrow Region End' are also not used. The column headers are: (1) CNV Identifier, (2) Chromosome, (3) Narrow Region Start, (4) Narrow Region End, (5) Wide Region Start, (6) Wide Region End

  • Amplification Threshold: Threshold for copy number amplifications. Regions with a log2 ratio above this value are considered amplified.

  • Deletion Threshold: Threshold for copy number deletions. Regions with a log2 ratio below the negative of this value are considered deletions.

  • Cap Values: Minimum and maximum cap values on analyzed data. Regions with a log2 ratio greater than the cap are set to the cap value; regions with a log2 ratio less than -cap value are set to -cap. Values must be positive.

  • Broad Length Cutoff: Threshold used to distinguish broad from focal events, given in units of fraction of chromosome arm.

  • Remove X-Chromosome: Flag indicating whether to remove data from the X-chromosome before analysis. Allowed values= {1,0} (1: Remove X-Chromosome, 0: Do not remove X-Chromosome.

  • Confidence Level: Confidence level used to calculate the region containing a driver.

  • Join Segment Size: Smallest number of markers to allow in segments from the segmented data. Segments that contain fewer than this number of markers are joined to the neighboring segment that is closest in copy number.

  • Arm Level Peel Off: Flag set to enable arm-level peel-off of events during peak definition. The arm-level peel-off enhancement to the arbitrated peel-off method assigns all events in the same chromosome arm of the same sample to a single peak. It is useful when peaks are split by noise or chromothripsis. Allowed values= {1,0} (1: Use arm level peel off, 0: Use normal arbitrated peel-off).

  • Maximum Sample Segments: Maximum number of segments allowed for a sample in the input data. Samples with more segments than this threshold are excluded from the analysis.

  • Gene GISTIC: When enabled (value = 1), this option causes GISTIC to analyze deletions using genes instead of array markers to locate the lesion. In this mode, the copy number assigned to a gene is the lowest copy number among the markers that represent the gene.

Values

List of inputs used for this run of GISTIC2. All files listed should be included in the archived results.

  • Segmentation File = /xchip/cga/gdac-prod/tcga-gdac/jobResults/GDAC_MergeDataFilesPipeline/KICH-TP/9825221/GDAC_MergeDataFiles_3125456/KICH-TP.snp__genome_wide_snp_6__broad_mit_edu__Level_3__segmented_scna_minus_germline_cnv_hg19__seg.seg.txt

  • Markers File = /xchip/cga/reference/gistic2/genome.info.6.0_hg19.na31_minus_frequent_nan_probes_sorted_2.1.txt

  • Reference Genome = /xchip/cga/reference/gistic2/hg19_with_miR_20120227.mat

  • CNV Files = /xchip/cga/reference/gistic2/CNV.hg19.bypos.111213.txt

  • Amplification Threshold = 0.1

  • Deletion Threshold = 0.1

  • Cap Values = 1.5

  • Broad Length Cutoff = 0.7

  • Remove X-Chromosome = 0

  • Confidence Level = 0.99

  • Join Segment Size = 4

  • Arm Level Peel Off = 1

  • Maximum Sample Segments = 2000

  • Gene GISTIC = 1

Table 3.  Get Full Table First 10 out of 66 Input Tumor Samples.

Tumor Sample Names
TCGA-KL-8323-01A-21D-2308-01
TCGA-KL-8324-01A-11D-2308-01
TCGA-KL-8325-01A-11D-2308-01
TCGA-KL-8326-01A-11D-2308-01
TCGA-KL-8327-01A-11D-2308-01
TCGA-KL-8328-01A-11D-2308-01
TCGA-KL-8329-01A-11D-2308-01
TCGA-KL-8330-01A-11D-2308-01
TCGA-KL-8331-01A-11D-2308-01
TCGA-KL-8332-01A-11D-2308-01

Figure 3.  Segmented copy number profiles in the input data

Output
All Lesions File (all_lesions.conf_##.txt, where ## is the confidence level)

The all lesions file summarizes the results from the GISTIC run. It contains data about the significant regions of amplification and deletion as well as which samples are amplified or deleted in each of these regions. The identified regions are listed down the first column, and the samples are listed across the first row, starting in column 10.

Region Data

Columns 1-9 present the data about the significant regions as follows:

  1. Unique Name: A name assigned to identify the region.

  2. Descriptor: The genomic descriptor of that region.

  3. Wide Peak Limits: The 'wide peak' boundaries most likely to contain the targeted genes. These are listed in genomic coordinates and marker (or probe) indices.

  4. Peak Limits: The boundaries of the region of maximal amplification or deletion.

  5. Region Limits: The boundaries of the entire significant region of amplification or deletion.

  6. Q values: The Q value of the peak region.

  7. Residual Q values: The Q value of the peak region after removing ('peeling off') amplifications or deletions that overlap other, more significant peak regions in the same chromosome.

  8. Broad or Focal: Identifies whether the region reaches significance due primarily to broad events (called 'broad'), focal events (called 'focal'), or independently significant broad and focal events (called 'both').

  9. Amplitude Threshold: Key giving the meaning of values in the subsequent columns associated with each sample.

Sample Data

Each of the analyzed samples is represented in one of the columns following the lesion data (columns 10 through end). The data contained in these columns varies slightly by section of the file. The first section can be identified by the key given in column 9 - it starts in row 2 and continues until the row that reads 'Actual Copy Change Given.' This section contains summarized data for each sample. A '0' indicates that the copy number of the sample was not amplified or deleted beyond the threshold amount in that peak region. A '1' indicates that the sample had low-level copy number aberrations (exceeding the low threshold indicated in column 9), and a '2' indicates that the sample had high-level copy number aberrations (exceeding the high threshold indicated in column 9).The second section can be identified the rows in which column 9 reads 'Actual Copy Change Given.' The second section exactly reproduces the first section, except that here the actual changes in copy number are provided rather than zeroes, ones, and twos.The final section is similar to the first section, except that here only broad events are included. A 1 in the samples columns (columns 10+) indicates that the median copy number of the sample across the entire significant region exceeded the threshold given in column 9. That is, it indicates whether the sample had a geographically extended event, rather than a focal amplification or deletion covering little more than the peak region.

Amplification Genes File (amp_genes.conf_##.txt, where ## is the confidence level)

The amp genes file contains one column for each amplification peak identified in the GISTIC analysis. The first four rows are:

  1. Cytoband

  2. Q value

  3. Residual Q value

  4. Wide Peak Boundaries

These rows identify the lesion in the same way as the all lesions file.The remaining rows list the genes contained in each wide peak. For peaks that contain no genes, the nearest gene is listed in brackets.

Deletion Genes File (del_genes.conf_##.txt, where ## is the confidence level)

The del genes file contains one column for each deletion peak identified in the GISTIC analysis. The file format for the del genes file is identical to the format for the amp genes file.

Gistic Scores File (scores.gistic)

The scores file lists the Q values [presented as -log10(q)], G scores, average amplitudes among aberrant samples, and frequency of aberration, across the genome for both amplifications and deletions. The scores file is viewable with the Genepattern SNPViewer module and may be imported into the Integrated Genomics Viewer (IGV).

Segmented Copy Number (raw_copy_number.{fig|pdf|png} )

The segmented copy number is a pdf file containing a colormap image of the segmented copy number profiles in the input data.

Amplification Score GISTIC plot (amp_qplot.{fig|pdf|png|v2.pdf})

The amplification pdf is a plot of the G scores (top) and Q values (bottom) with respect to amplifications for all markers over the entire region analyzed.

Deletion Score GISTIC plot (del_qplot.{fig|pdf|png|v2.pdf})

The deletion pdf is a plot of the G scores (top) and Q values (bottom) with respect to deletions for all markers over the entire region analyzed.

Tables (table_{amp|del}.conf_##.txt, where ## is the confidence level)

Tables of basic information about the genomic regions (peaks) that GISTIC determined to be significantly amplified or deleted. These describe three kinds of peak boundaries, and list the genes contained in two of them. The region start and region end columns (along with the chromosome column) delimit the entire area containing the peak that is above the significance level. The region may be the same for multiple peaks. The peak start and end delimit the maximum value of the peak. The extended peak is the peak determined by robust, and is contained within the wide peak reported in {amp|del}_genes.txt by one marker.

Broad Significance Results (broad_significance_results.txt)

A table of per-arm statistical results for the data set. Each arm is a row in the table. The first column specifies the arm and the second column counts the number of genes known to be on the arm. For both amplification and deletion, the table has columns for the frequency of amplification or deletion of the arm, and a Z score and Q value.

Broad Values By Arm (broad_values_by_arm.txt)

A table of chromosome arm amplification levels for each sample. Each row is a chromosome arm, and each column a sample. The data are in units of absolute copy number -2.

All Data By Genes (all_data_by_genes.txt)

A gene-level table of copy number values for all samples. Each row is the data for a gene. The first three columns name the gene, its NIH locus ID, and its cytoband - the remaining columns are the samples. The copy number values in the table are in units of (copy number -2), so that no amplification or deletion is 0, genes with amplifications have positive values, and genes with deletions are negative values. The data are converted from marker level to gene level using the extreme method: a gene is assigned the greatest amplification or the least deletion value among the markers it covers.

Broad Data By Genes (broad_data_by_genes.txt)

A gene-level table of copy number data similar to the all_data_by_genes.txt output, but using only broad events with lengths greater than the broad length cutoff. The structure of the file and the methods and units used for the data analysis are otherwise identical to all_data_by_genes.txt.

Focal Data By Genes (focal_data_by_genes.txt)

A gene-level table of copy number data similar to the all_data_by_genes.txt output, but using only focal events with lengths greater than the focal length cutoff. The structure of the file and the methods and units used for the data analysis are otherwise identical to all_data_by_genes.txt.

All Thresholded By Genes (all_thresholded.by_genes.txt)

A gene-level table of discrete amplification and deletion indicators at for all samples. There is a row for each gene. The first three columns name the gene, its NIH locus ID, and its cytoband - the remaining columns are the samples. A table value of 0 means no amplification or deletion above the threshold. Amplifications are positive numbers: 1 means amplification above the amplification threshold; 2 means amplifications larger to the arm level amplifications observed for the sample. Deletions are represented by negative table values: -1 represents deletion beyond the threshold; -2 means deletions greater than the minimum arm-level deletion observed for the sample.

Sample Cutoffs (sample_cutoffs.txt)

A table of the per-sample threshold cutoffs (in units of absolute copy number -2) used to distinguish the high level amplifications (+/-2) from ordinary amplifications (+/-1) in the all_thresholded.by_genes.txt output file. The table contains three columns: the sample identifier followed by the low (deletion) and high (amplification) cutoff values. The cutoffs are calculated as the minimum arm-level amplification level less the deletion threshold for deletions and the maximum arm-level amplification plus the amplification threshold for amplifications.

Focal Input To Gistic (focal_input.seg.txt)

A list of copy number segments describing just the focal events present in the data. The segment amplification/deletion levels are in units of (copy number -2), with amplifications positive and deletions negative numbers. This file may be viewed with IGV.

Gene Counts vs. Copy Number Alteration Frequency (freqarms_vs_ngenes.{fig|pdf})

An image showing the correlation between gene counts and frequency of copy number alterations.

Confidence Intervals (regions_track.conf_##.bed, where ## is the confidence level)

A file indicating the position of the confidence intervals around GISTIC peaks that can be loaded as a track in a compatible viewer browser such as IGV or the UCSC genome browser.

GISTIC

GISTIC identifies genomic regions that are significantly gained or lost across a set of tumors. It takes segmented copy number ratios as input, separates arm-level events from focal events, and then performs two tests: (i) identifies significantly amplified/deleted chromosome arms; and (ii) identifies regions that are significantly focally amplified or deleted. For the focal analysis, the significance levels (Q values) are calculated by comparing the observed gains/losses at each locus to those obtained by randomly permuting the events along the genome to reflect the null hypothesis that they are all 'passengers' and could have occurred anywhere. The locus-specific significance levels are then corrected for multiple hypothesis testing. The arm-level significance is calculated by comparing the frequency of gains/losses of each arm to the expected rate given its size. The method outputs genomic views of significantly amplified and deleted regions, as well as a table of genes with gain or loss scores. A more in depth discussion of the GISTIC algorithm and its utility is given in [1], [3], and [5].

CNV Description

Regions of the genome that are prone to germ line variations in copy number are excluded from the GISTIC analysis using a list of germ line copy number variations (CNVs). A CNV is a DNA sequence that may be found at different copy numbers in the germ line of two different individuals. Such germ line variations can confound a GISTIC analysis, which finds significant somatic copy number variations in cancer. A more in depth discussion is provided in [6]. GISTIC currently uses two CNV exclusion lists. One is based on the literature describing copy number variation, and a second one comes from an analysis of significant variations among the blood normals in the TCGA data set.

Download Results

In addition to the links below, the full results of the analysis summarized in this report can also be downloaded programmatically using firehose_get, or interactively from either the Broad GDAC website or TCGA Data Coordination Center Portal.

References
[1] Beroukhim et al, Assessing the significance of chromosomal aberrations in cancer: Methodology and application to glioma, Proc Natl Acad Sci U S A. Vol. 104:50 (2007)
[3] Mermel et al, GISTIC2.0 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers, Genome Biology Vol. 12:4 (2011)
[5] Beroukhim et al., The landscape of somatic copy-number alteration across human cancers, Nature Vol. 463:7283 (2010)
[6] McCarroll, S. A. et al., Integrated detection and population-genetic analysis of SNPs and copy number variation, Nat Genet Vol. 40(10):1166-1174 (2008)