# Genome-wide identification of significant aberrations in cancer genome

- Xiguo Yuan†
^{1, 2}, - Guoqiang Yu†
^{2, 3}, - Xuchu Hou
^{2}, - Ie-Ming Shih
^{4, 8}, - Robert Clarke
^{5}, - Junying Zhang
^{1}, - Eric P Hoffman
^{6}, - Roger R Wang
^{7}, - Zhen Zhang
^{8}and - Yue Wang
^{2}Email author

**13**:342

**DOI: **10.1186/1471-2164-13-342

© Yuan et al.; licensee BioMed Central Ltd. 2012

**Received: **1 February 2012

**Accepted: **27 July 2012

**Published: **27 July 2012

## Abstract

### Background

Somatic Copy Number Alterations (CNAs) in human genomes are present in almost all human cancers. Systematic efforts to characterize such structural variants must effectively distinguish significant consensus events from random background aberrations. Here we introduce Significant Aberration in Cancer (SAIC), a new method for characterizing and assessing the statistical significance of recurrent CNA units. Three main features of SAIC include: (1) exploiting the intrinsic correlation among consecutive probes to assign a score to each CNA unit instead of single probes; (2) performing permutations on CNA units that preserve correlations inherent in the copy number data; and (3) iteratively detecting Significant Copy Number Aberrations (SCAs) and estimating an unbiased null distribution by applying an SCA-exclusive permutation scheme.

### Results

We test and compare the performance of SAIC against four peer methods (GISTIC, STAC, KC-SMART, CMDS) on a large number of simulation datasets. Experimental results show that SAIC outperforms peer methods in terms of larger area under the *Receiver Operating Characteristics* curve and increased detection power. We then apply SAIC to analyze structural genomic aberrations acquired in four real cancer genome-wide copy number data sets (ovarian cancer, metastatic prostate cancer, lung adenocarcinoma, glioblastoma). When compared with previously reported results, SAIC successfully identifies most SCAs known to be of biological significance and associated with oncogenes (e.g., KRAS, CCNE1, and MYC) or tumor suppressor genes (e.g., CDKN2A/B). Furthermore, SAIC identifies a number of novel SCAs in these copy number data that encompass tumor related genes and may warrant further studies.

### Conclusions

Supported by a well-grounded theoretical framework, SAIC has been developed and used to identify SCAs in various cancer copy number data sets, providing useful information to study the landscape of cancer genomes. Open–source and platform-independent SAIC software is implemented using C++, together with R scripts for data formatting and Perl scripts for user interfacing, and it is easy to install and efficient to use. The source code and documentation are freely available at http://www.cbil.ece.vt.edu/software.htm.

## Background

Somatic copy number alterations (CNAs) are common genetic events in the development and progression of various human cancers, and significantly contribute to tumorigenesis [1, 2]. The coverage of CNAs in tumors varies from a few hundred to several million nucleotide bases, consisting of both deletions and amplifications with highly complex patterns [3, 4]. Recent advances in oligonucleotide-based single nucleotide polymorphism (SNP) arrays have made it possible to detect regional amplifications and deletions with high resolution on a genome-wide scale [5, 6]. A critical challenge in the genome-wide analysis of CNAs is to distinguish between the “driver” mutations that allow the tumor to initiate, grow, and persist, and the “passenger” mutations that represent random somatic events accumulated during tumorigenesis [1, 3, 7]. Identification of these “driver” alterations can provide important insights into the cellular defects that cause cancer and suggest potential diagnostic, prognostic, and targeted therapeutic strategies [1, 7, 8].

By studying a sufficiently large collection of cancer samples, Significant Copy Number Aberrations (SCAs), defined as significantly recurrent CNAs that affect the same region in multiple tumors, are widely considered as informative surrogates of “driver” mutations that may help pinpoint novel cancer-causing genes [3, 9]. Past studies have detected many SCAs in a wide range of cancer types, with an impressive coverage of many known oncogenes and cancer suppressor genes [1, 2, 7]. Several methods for finding regions of SCAs using CNAs data have been described in the literature, where the task of distinguishing between sporadic CNAs and SCAs is largely a statistical significance testing. Two reviews with qualitative comparison of different methods have been published [10, 11]. Despite the use of different algorithms, a common theme in these methods is that they often adopt a four-step strategy: (1) detect CNAs and separate deletions and amplifications; (2) design and calculate ensemble test statistics associated with a genomic locus; (3) construct and/or estimate the probability distribution of test statistics under the null hypothesis; (4) perform multiple testing on a pool of genomic loci.

Significance testing for aberrant copy number (STAC) starts by converting the normalized log-ratios into a binary matrix, with zeros indicating no change and ones indicting losses and gains [12]. STAC then proposes two statistics (footprint and frequency) to define regions of SCAs while adjusting for multiple comparisons, where the null hypothesis is that the detected CNAs from single-sample analysis are the realizations of random CNA placements whose probability distribution is generated by permutations on CNA segments [13]. Genomic Identification of Significant Targets in Cancer (GISTIC) works on the real-valued step function of log-ratios that allows GISTIC to exploit both the type (amplification/deletion) and amplitude of CNAs [1, 3]. Using a semi-parametric permutation assuming independence between probes, GISTIC calculates a score that is based on both the amplitude and frequency of CNAs at each probe position and subsequently identify regions of SCAs, where amplification and deletion CNAs are handled separately, and armed-level and focal CNAs are further analyzed independently [14]. Aimed to correlate information from neighboring probes with the amplitude and frequency of CNAs at each probe position, Kernel Convolution – a Statistical Method for Aberrant Regions detection (KC-SMART) uses varying-width kernel functions to calculate the testing statistics from the original log-ratios across multiple samples, producing the kernel smoothed estimate (KSE) at each locus by locally weighted regression [15]. SCAs are selected based on a permutation-generated null distribution and Bonferroni correction. To substantially reduce computational burden in analyzing high-resolution and large-population data, correlation matrix diagonal segmentation (CMDS) identifies SCAs based on a between-chromosomal-site correlation analysis directly using the raw intensity ratios across all samples [16]. CMDS uses a correlation statistics to detect SCAs with a standard normal null distribution whose parameters are estimated directly from the data and adjusts for multiple comparisons by false discovery rate.

Existing methods have several limitations. When working with unprocessed raw intensity ratios [13, 15, 16], most methods are oblivious to noise clutter that can significantly confound estimation of the null distribution about true yet sporadic CNAs [9, 17]. Furthermore, these methods cannot distinguish between contributions of amplifications and deletions to the calculated overall test statistics that may affect the power to detect SCAs. While some effort has been made to incorporate correlation among neighboring probes into the test statistics, most methods assign a score to, and test the significance at, each individual probe locus [14, 15]. In addition, while it is widely accepted that CNAs signals at adjacent probes are highly correlated [9, 13–15], the assumption of probe independence is often adopted in constructing and learning the null distribution, probably for mathematical convenience [3, 16]. Moreover, existing permutation experiments using multiple samples cannot distinguish between the contributions of sporadic CNAs (obeying null distribution) and actual SCAs (deviating from null distribution) to the estimation of null distributions, resulting in theoretically conservative estimations especially when the number of true SCAs participating in the permutation is large.

We now report Significant Aberration in Cancer (SAIC), a carefully motivated method for accurately identifying SCAs using CNAs data from multiple samples. To distinguish between different biological roles of CNAs types and between noise and sporadic CNAs, we use discretized CNAs data and separately analyze copy number amplifications and deletions. By exploiting the intrinsic correlation among consecutive probes, we calculate and assign a score (test statistics) to each CNA unit instead of each single probe, based on both the amplitude and frequency of CNAs within the unit. To accurately estimate the null distribution governing sporadic CNAs, we perform random positional permutations on CNA units that preserve correlations inherent to the copy number data. More importantly, to minimize the unwanted participation of true SCAs in determining the null distribution [3, 14], we iteratively detect SCAs and estimate an unbiased null distribution by an SCA-exclusive permutation scheme.

We tested SAIC on extensive simulation data sets, observing significantly improved performance with larger areas under the *Receiver Operating Characteristics* (ROC) curves and higher sensitivities at acceptable low false discovery rates, as compared to four popular peer methods (GISTIC, STAC, KC-SMART, and CMDS). We then applied SAIC to four real benchmark data sets, successfully identified the majority (84%) of previously reported SCAs harboring regions associated with well-known tumor-causing genes, and more importantly, detected some novel SCAs partially validated by the presence of known cancer-related genes.

## Methods

### Data format and definitions

*N*×

*M*matrix

*X*. Each entry

*x*

_{ nm }represents DNA copy number (in log2-ratio) for sample

*n*at probe

*m*, where each row

*X*

_{ n }corresponds to copy number for

*n*th sample at

*M*probes. Copy number amplifications and deletions are analyzed separately. We use the indicator function to divide matrix

*X*into two matrices

*X*=

*X*

_{amplification}+

*X*

_{deletion}, where

with *θ*_{amplification} and *θ*_{deletion} being the pre-specified thresholds. For brevity, we focus all subsequent discussion on *X*_{amplification} and make comments on *X*_{deletion} when necessary.

#### Definition 1

Any copy number probe *m* whose associated copy number is amplified or deleted in at least one of *N* samples is called a CNA probe.

*ρ*

_{ij}between CNA probes

*i*and

*j*is then calculated for $\left\{i\ne j\right\}\in M$:

*s*

_{ i }and

*s*

_{ j }are the estimated means and standard deviations of copy numbers at probes

*i*and

*j*across

*N*samples, respectively. If

*ρ*

_{ij}is less than a pre-specified threshold

*θρ*, a breakpoint occurs between probes

*i*and

*j*.

#### Definition 2

A sequence of consecutive CNA probes with no breakpoints is defined as a CNA unit, denoted by *u* (*k, L*) with *k* being the starting probe index and *L* being the length of the CNA unit.

Intuitively, a CNA unit consists of a sequence of highly correlated consecutive CNA probes. Figure 1 illustrates the concepts of CNA region and CNA unit, where two CNA regions contain 10 and 3 CNA probes, respectively, and the first CNA region is further split into three CNA units due to two breakpoints within the CNA region.

### Summary statistics and significance assessment

*U*score) that incorporates both the frequencies of occurrence and the amplitudes of the CNA probes within the unit, leading to the unit summary statistics given by

Second, the method assesses the statistical significance of each CNA unit by comparing the observed statistic to the *U* scores that would be expected by chance.

Sporadic CNA units often occur throughout the genome, so a null distribution for *U*_{
k, L
} under the hypothesis that no SCAs are present, can be estimated by randomly permuting the overall pattern of presumed all-sporadic CNA units across the genome [3, 9, 12, 15]. Though various permutation schemes can be adopted, due to different rates of CNA and different percentages of normal tissue contamination in tumor samples [18], permutation of CNA units across rows/samples should be avoided. As aforementioned, permutation should be performed on CNA units (instead of single CNA probes) that preserve correlations inherent to the copy number data, even if the CNA units are sporadic [3, 9, 15]. Another subtle but conveniently ignored issue is the different background rates of CNA units with varying lengths [1]. Short CNA units occur at a frequency inversely related to their lengths and long CNA units occur approximately 30 times more frequently than would be expected by the inverse-length distribution. This observation is seen across all cancer types, is applicable to both copy gains and losses, and is supported by the calculated genome-average background rates for CNAs as a function of length [1]. These considerations motivate our carefully designed SAIC permutation scheme.

Let $\mathbb{L}$ denote the integer set containing the lengths of all the observed CNA units in *X*, $\mathbb{K}$ denote the integer set containing the starting probe indices of all the observed CNA units in *X*, and *X*^{(t)} be the random positional permutation of *X* for *t* = 1,2,…,*T*, with *T* being the total number of permutations. We now describe our method for analyzing CNA units for evidence of significant alteration in cancer, where we account for the difference in background rates between CNA units of different lengths by considering them adaptively.

#### Algorithm 1

*U*

_{ k, L }

- (1)
Perform

*T*random within-row positional permutations*X*^{(1)},*X*^{(2)}, …,*X*^{(T)}of the data matrix*X*on CNA units; - (2)
Compute the value of summary statistic ${U}_{k\text{,}\phantom{\rule{0.25em}{0ex}}L}\left({X}^{\left(t\right)}\right)$ for each permuted data set

*t*= 1,2,…,*T*, and for each starting probe $k=1,2,\dots ,M-L+1$ and each length $L\in \mathbb{L}$; - (3)Calculate and assign a P-value to each observed CNA unit
*u*(*k, L*) for $k\in \mathbb{K}$ based on the extreme right-hand tail probability given by [9, 19]$P\left({U}_{k\text{,}\phantom{\rule{0.25em}{0ex}}L}\left(X\right)\right)=\frac{1+\sum _{t=1}^{T}I\left(\underset{k\text{'}}{max}{U}_{k\text{'}\text{,}\phantom{\rule{0.25em}{0ex}}L}\left({X}^{\left(t\right)}\right)\ge {U}_{k\text{,}\phantom{\rule{0.25em}{0ex}}L}\left(X\right)\right)}{T+1}\text{,}$(4)

where $I\left(\xb7\right)$is the indicator function.

The empirical P-values on *X*_{deletion} are calculated by the extreme left-hand tail probabilities and reversing the inequality in Eq. (4). Both definitions produce P-values that are easy to interpret, and the “max” operation automatically adjusted P-values for multiple comparisons across CNA units thus controls the family-wise error rate [9].

In algorithm 1, it is important to note that when we generate a randomly permuted dataset based on the observed data, we do not re-define the CNA units but re-use the already-defined CNA units. Specifically, in each permutation, we randomly place the already-defined CNA units over the whole genome or each chromosome within each sample, and calculate the summary *U* score for each length of CNA units. Thus, independent of the unit length, the observed CNA units will always be retained (implicitly) in the permuted dataset. Moreover, when the number of permutations is sufficiently large, the p-values of observed CNA units can be accurately estimated. More precisely, to assess the p-value associated with an observed CNA unit of length *L*, we calculate the *U* scores for any consecutive *L* probes (probes do not need to reside within the same unit) across the genome, and compare the maximum score with the score of the observed CNA unit.

### Iterative estimation of unbiased null distribution

One important issue concerning Algorithm 1 is the presence of true SCAs (departing from null distribution) in cancer genomes that presumably contribute high copy number deviations to the estimation of overall null distribution (governing only sporadic CNAs), potentially reducing power to detect less-extreme SCAs due to theoretical conservativeness [9, 14]. Loss of power is particularly critical in real-world applications where the number of true SCAs in cancer genomes may be large. Thus, to minimize the unwanted participation of true SCAs in determining the null distribution, we iteratively detect SCAs and estimate an unbiased null distribution by applying an SCA-exclusive permutation scheme. SAIC assesses the ‘new’ SCAs conditional on having found the ‘existing’ SCAs, successively correcting for true SCAs in order to better dissect and detect SCAs. Specifically, the CNA units associated with the ‘existing’ SCAs are masked as zeros after each iteration, resulting in a new data set *X*_{-SCAs} in which already-detected SCAs becomes null.

#### Algorithm 2

*U*

_{ k, L }

- (1)
Perform Algorithm 1;

- (2)
Check whether ‘new’ SCAs are detected. If ‘yes’, continue; if “no”, stop and re-calculate the P-values for all SCAs using truth converging null distribution;

- (3)
Mask the CNA units associated with newly detected SCAs as zeros and let $X={X}_{-\text{SCAs}}$, then go to step (1).

It has been shown experimentally that additional power to detect SCAs can be gained by removing the effect of newly detected SCAs after each iteration [9]. However, an iterative SCA-exclusive permutation scheme raises another subtle yet critical issue concerning the convergence of null distribution learning and potential bias due to the expected false positive SCAs under the truth-converging null distribution. Fortunately, based on the careful design of Algorithm 2, the following theorem shows that, if we apply a significance level $\alpha \text{'}=\alpha /\left(1+\alpha \right)$ where *α* is the targeted false positive rate (FPR), an unbiased estimation and detection results can be readily obtained using Algorithm 2 (see formal proof in Appendix A).

#### Theorem 1

*Suppose that Algorithm 2 is used to iteratively detect SCAs and estimate truth converging null distribution*. *Let α* *be the targeted FPR and*$\alpha \text{'}=\alpha /\left(1+\alpha \right)$*be the significance level used to detect SCAs. Then an unbiased truth converging null distribution can be obtained together with a theoretical FPR α.*

### SAIC algorithm and data preprocessing

## Results

In the absence of definitive ground truth about the recurrent CNAs in the cancer genomes, the validation of a new method for detecting SCAs is always problematic [9, 13, 16, 18, 24]. We first validate SAIC on multiple realistic simulation data sets and then proceed to evaluate the method using real CNA data sets. All data sets were analyzed according to the algorithm described in Figure 2. We tested SAIC and the four peer methods (GISTIC, STAC, KC-SMART, CMDS) on realistic simulation data sets. Comparative performance was based on the ground truth in terms of detection power [18] and the *Receiver Operating Characteristics* (ROC) curves [16]. When applied to real CNA data, we compared and discussed biological plausibility of the implicated SCAs, and examined relative SCAs coverage between SAIC and GISTIC on benchmark data sets using Venn diagrams. To assure a meaningful and differential comparison, we emphasized experiment suitability when choosing algorithm parameter settings. For example, the algorithm parameter settings cannot be too “simple” (if there are only a few arm-level SCAs, all methods may perform equally well) or too “complex” (if there are many weak focal SCAs, no method will perform consistently well) [14].

### Simulation studies

*U*

_{ k, L }(

*X*) in

*X*that was significant at

*α*= 0.05 level [9]. Values of the observed FWER in Table 1 (0.0497) suggest that SAIC is almost perfect when compared with slightly conservative values (0.0452) by similar method [9].

**Empirical type 1 error rate for simulated data sets under the null hypothesis**

Null simulation model | Empirical FWER atα= 0.05level |
---|---|

Copy number data | 0.0488 |

Clumped copy number data (25%) | 0.0500 |

Clumped copy number data (50%) | 0.0493 |

Clumped copy number data (75%) | 0.0505 |

*N*= 40 ~ 80 samples and each sample contains

*M*= 5,000 probes. To replicate the effect of inevitable normal cell contamination [18], the copy numbers at every probes are simulated by a mixture of normal and tumor genomes, where the normal cell fraction

*λ*is randomly drawn from a normal distribution $\mathcal{N}\left({\mu}_{\lambda},{\sigma}_{\lambda}\right)$ with

*μ*

_{ λ }and

*σ*

_{ λ }being the mean and standard deviation of normal cell fraction in the sample. Each sample contains two sporadic CNA regions, one deletion and one amplification randomly drawn from integer sets {0, 1} and {3, 4,…,8}, respectively. Each data set contains two recurrent CNA regions that are contributed from a fraction of samples according to a specified frequency

*ω*, one deletion and one amplification similarly designed as aforementioned. The length of both sporadic and recurrent CNA regions is randomly assigned from 150 to 250 probes, realistically reflecting the estimated background rate of focal CNAs in a typical cancer sample genome [1]. To equally assess the power in detecting deletion or amplification SCAs, we calculate the detection power of SAIC or GISTIC as the rate of successfully detecting inserted, deleted or amplified SCAs across 100 data sets. Table 2 summarizes the comparative detection power of SAIC and GISTIC for a total of 19 parameter settings across 1,900 data sets. These comparative experimental results consistently show that SAIC outperforms GISTIC with significantly increased detection power in 18 out of 19 simulations.

**Power to detect SCAs by SAIC and GISTIC in simulation studies**

N = 60,ω = 0.2,μ | 0.15 | 0.2 | 0.25 | 0.3 | 0.35 |
---|---|---|---|---|---|

GISTIC | 89% | 86% | 79% | 74% | 72% |

SAIC | 96% | 94% | 86% | 86% | 82% |

| 0.4 | 0.5 | 0.6 | 0.7 | 0.8 |

GISTIC | 83% | 81% | 82% | 72% | 79% |

SAIC | 93% | 91% | 87% | 79% | 74% |

| 40 | 50 | 60 | 70 | 80 |

GISTIC | 58% | 73% | 79% | 86% | 89% |

SAIC | 65% | 83% | 87% | 93% | 94% |

| 0.1 | 0.15 | 0.2 | 0.25 | |

GISTIC | 30% | 58% | 80% | 92% | |

SAIC | 37% | 72% | 87% | 97% |

*N*= 50 samples and each sample contains

*M*= 5,000 probes. The log-ratios at every probe are simulated by a mixture of normal and tumor genomes, with the normal cell fraction

*λ*being randomly drawn from a uniform distribution $\mathcal{U}\left(0.2,0.8\right)$. Zero-mean Gaussian noise is then added to each sample with three levels of standard deviation

*σ*randomly drawn from uniform distributions $\mathcal{U}\left(0.2,0.4\right)$$\mathcal{U}\left(0.4,0.6\right)$, and $\mathcal{U}\left(0.6,0.8\right)$. To make the simulations more realistic, for each simulated sample genome, we insert 2 to 10 randomly located background CNA regions with the lengths varying from 10 to 50 probes. There are three ‘amplification’ (

*L*= 30, 20, 10) and one ‘deletion’ (

*L*= 20) ground truth SCAs embedded in each of the simulation data sets with a baseline frequency

*ω*= 0.1. The copy numbers associated with amplification SCAs are 3, 4 and 5, and deletion SCAs are 0 and 1. In our simulation software, we use two parameters

*β*

_{ L }and

*β*

_{ ω }to modify the length and frequency of these SCAs. Other parameter settings include

*θ*

_{ ρ }= 0.75,

*θ*

_{amplification}= 0.1 and

*θ*

_{deletion}= −0.1 (default setting by GISTIC and CBS) for defining CNAs probes and units. Based on the estimated true positive rate (TPR) and corresponding FPR at different significance levels, Figure 3 presents ROC curves of SAIC and peer methods derived from the simulation studies. These comparative experimental results consistently show that SAIC outperforms the peer methods in terms of larger areas

*A*

_{ z }under the ROC curves or increased sensitivity at low FPR. More simulation studies are given in Additional file 1, where we report the power in detecting the boundaries of SCAs by these methods, and once again, showing outperformance of SAIC as compared to the peer methods [3, 14].

### Application to four real cancer copy number data sets

We applied SAIC to four real cancer copy number data sets and identified many SCAs that encompass established or potentially novel cancer ‘driver’ genes. The data sets are from ovarian cancer [26, 27], prostate cancer [8, 18], lung adenocarcinoma [1, 7], and glioblastoma [1, 3]. Due to their distinct biological functions in cancer development, SAIC analyzes separately chromosomes 1–22 and chromosome X/Y. To account for the different background CNA rates across chromosomes, we identify SCAs by performing SAIC on individual chromosomes. Other parameter settings include *T* = 1000 and *α* = 0.05 (theoretical significance level or FPR/FWER). To provide a somewhat independent verification, we compared the SCAs detected by SAIC with what reproduced by GISTIC on lung adenocarcinoma and glioblastoma data sets that have been previously reported [3, 7].

### Results on the ovarian cancer data set

*N*= 63 tumor samples [26–28]. Copy number signals were acquired using the Affymetrix Human Mapping 250 K Sty SNP Array platform [1]. Each sample contains a total of 238,230 probes across the whole genome. Other algorithm parameter settings include

*θ*

_{ ρ }= 0.95,

*θ*

_{amplification}= 0.263 (2.4 copies) and

*θ*

_{deletion}= −0.322 (1.6 copies) [14]. The genome-wide landscapes (via -log

_{10}

*P*) of recurrent or sporadic CNAs observed in the data sets are given in Figure 4, where amplifications and deletions are separately shown (left and right sides). SAIC detected several SCAs (both amplification and deletion), many of which are biologically plausible and include known oncogenes (e.g., KRAS, CCNE1 and CCND2) and tumor suppressor genes (e.g., CDKN2A and CDKN2B) [26, 27, 29, 30]. Full lists of the genes covered by these SCAs are given in Additional file 2 (ST 2). SAIC also identified many other cancer driver genes within individual chromosomes (ST 3), such as SKIL, CDK4, PIK3CA, PTEN, FGD4, FGFR1.

### Results on the metastatic prostate cancer dataset

*N*= 55 clustered metastatic tumor samples, obtained from 13 prostate cancer patients. Copy number signals were acquired using Affymetrix Genome-Wide Human SNP Array 6.0 [8, 18]. Each sample contains a total of 1,868,857 probes across the whole genome. To discount the potential bias due to imbalanced subject-cluster sampling [8], we chose to analyze the

*N*= 13 representative samples and to detect global recurrent CNAs by SAIC. Other algorithm parameter settings include

*θ*

_{ ρ }= 0.95,

*θ*

_{amplification}= 0.263 and

*θ*

_{deletion}= −0.322, the same as used in analyzing ovarian cancer data. The genome-wide landscape of recurrent or sporadic CNAs observed in metastatic prostate cancer data is given in Figure 5, where SAIC detected 15 amplification SCAs (318 genes) and 21 deletion SCAs (756genes). Full list of the genes covered by these SCAs are given in Additional file 3 (ST 4). Many of these genes are cancer related (e.g., EGFR, BRCA2, TP53, ATBF1, MYC and RB1). In individual chromosome analysis of the data set, SAIC identified many other SCAs involved with cancer driver genes, such as PTEN (ST 5).

### Results on the lung adenocarcinoma and glioblastoma datasets

*N*= 371 tumor samples, publicly available at http://www.broad.mit.edu/cancer/pub/tsp[7]. Copy number signals were acquired using Affymetrix 250K Sty SNP Array, where each sample contains a total of 216,327 probes across the whole genome [7]. To assure the general comparability of the results produced by SAIC and GISTIC, we adopted similar algorithm parameter settings used by GISTIC for detecting focal SCAs:

*θ*

_{amplification}= 0.848 and

*θ*

_{deletion}= −1.15, in addition to

*θ*

_{ ρ }= 0.9. The genome-wide landscape of recurrent or sporadic CNAs observed in lung adenocarcinoma data is given in Figure 6, where SAIC detected 23 amplification SCAs and 26 deletion SCAs (after combining some of 98 recurrent CNAs within the same cytobands). Full list of the genes covered by these SCAs is given in Additional file 4 (ST 6). The Venn diagram in Figure 7 reveals the numbers of common and distinctive SCAs detected by SAIC and GISTIC. It can be seen that SAIC successfully detected most (87% amplification and 75% deletion regions) of the SCAs that have been detected by GISTIC, while also revealing many additional SCAs (10 amplification and 23 deletion regions) [7]. In addition, the result from within-chromosome analysis of the data set is listed in Additional file 4 (ST 7).

*N*= 141 tumor samples, publicly available at http://www.broad.mit.edu/cancer/pub/GISTIC, where each sample contains a total of 115,593 probes across the whole genome [3]. Once again, we adopted the similar algorithm parameter settings used by GISTIC for detecting focal SCAs. The genome-wide landscape of recurrent or sporadic CNAs observed in glioblastoma data is given in Figure 8, where SAIC detected 15 amplification SCAs and 30 deletion SCAs (after combining some of 67 recurrent CNAs within the same cytobands). Full list of the genes covered by these SCAs are given in Additional file 5 (ST 8). The Venn diagram in Figure 9 reveals the numbers of common and distinctive SCAs detected by SAIC and GISTIC. It can be seen that SAIC successfully detected most (88% amplification and 75% deletion regions) of the SCAs that have been detected by GISTIC, while it also revealed many additional SCAs (8 amplification and 27 deletion regions) [3]. In addition, the result from within-chromosome analysis of the data set is listed in Additional file 5 (ST 9).

The common SCAs regions (e.g., 7p11.2, 12p12.1, 9p21.3, etc.) are highly consistent with previous reports, and largely encompass well-known oncogenes or tumor suppressor genes. For example, EGFR (epidermal growth factor receptor) is an oncogene within 7p11.2 whose mutations or amplifications have been shown to contribute to uncontrolled cell division (a predisposition for cancer) [31]. Many additional SCAs regions (e.g., 8p23.2, 21q22.2) contain or adjacent to disease-related genes (e.g., CSMD1 and TMPRSS3) and may warrant further study.

## Discussion

SAIC is similar to many peer methods in that it assesses statistical significance of SCAs using a permutation-based null distribution [9, 12, 14–16]. However, in contrast to the existing procedures, the CNA units used by SAIC preserve the essential correlation structures of serial probes whose estimated average correlation coefficient can be as high as 0.985 [32]. Moreover, by automatically adjusting P-values for multiple comparisons [33, 34] and iteratively re-estimating the null distribution exclusive of detected SCAs [9], SAIC can preserve the intrinsic false positive rate, without compromising detection power to resort to sometimes overly conservative schemes [3, 14–16]. Theoretic analysis and extensive experimental results show that SAIC preserves both type 1 error and detection power, see Tables 12. Furthermore, the novel concept of CNA unit and associated scoring and permutation scheme neatly parallels many considerations in the revised GISTIC2.0 [14], for example, serial probes covering driver events should be more highly correlated than probes covering only passengers and thus more likely to identify the target genes. The flexible length-adaptive significance assessment of CNA units via Eq. (4) automatically accounts for distinct background rates according to their lengths and thus more likely to detect independent SCAs.

As for the *θ*_{amplification} and *θ*_{deletion} parameters in the SAIC algorithm, there is no general guideline about how to select their values [14], since different types of cancers usually have different rates and magnitudes of background CNAs [14, 26, 35]. In addition, various degrees of normal cell contamination [18] and intratumor heterogeneity [35, 36] occur in many samples and these further complicate the selection of parameter values. Practically, lower thresholds were used to define broad (arm-level) CNAs while higher thresholds were used to define focal CNAs [3, 14]. A newly proposed strategy is to apply joint magnitude-length thresholds [14] and to correct normal cell contamination using BACOM [18]. Since our main objective here is to identify focal CNAs, we have largely adopted the same strategy used in [3, 14], i.e., we used relatively higher thresholds to define focal CNAs for subsequent analyses. Specifically, based on the observation that the magnitude of CNAs in ovarian and prostate cancers is relatively low, we used relatively lower and commonly used thresholds (2.0 ± 0.4), i.e., 2.4 copies for amplification and 1.6 copies for deletion. In contrast, on the datasets of lung adenocarcinoma and glioblastoma, we applied relatively higher thresholds (2.0 + 1.6, 2.0–1.1), i.e., 3.6 copies for amplification and 0.9 copies for deletion, that are similar to the thresholds used by GISTIC algorithms [3, 14].

Similar situation occurs to the selection of *θ*_{
ρ
} in defining CNA units [9]. Lower values of *θ*_{
ρ
} often produce longer CNA units while higher values of *θ*_{
ρ
} often produce shorter CNA units. It has been reported that the average successive probe correlation of the segmented data can be as high as 0.985 [9, 32]. In our experience in analyzing real cancer datasets, a value of ${\theta}_{\rho}$ taking between 0.7 and 0.95 would be a suitable choice.

It is important to note that the general conclusion on the relative performance of our SAIC and peer methods, at least based on the extensive simulation studies, remains largely true. We have used the same parameter values in all methods so that a fair comparison on their relative performances can be assured. Based on our analysis of real datasets using current parameter settings, it appears that SAIC performs well when compared to peer methods. In addition, the results of extensive simulation studies, performed under a variety of probe correlation schemes, show that SAIC preserves well the expected type 1 error, even when the probes follow non-stationary correlation structures similar to those found in real data [9].

SAIC currently can perform either genome-wide (except X/Y chromosome due to its distinct biological role) or chromosome-based CNA unit permutations. In the application of SAIC to real cancer data sets, we performed genome-wide, autosome-based, and X/Y-chromosome-based permutations. The combined results from using different permutation schemes contain more SCAs that may involve novel cancer driver genes. By exploiting the novel concepts of CNA probe, CNA unit, and multiscale permutation, experimental results show that SAIC can accurately detect the boundaries of SCAs with different lengths, see Additional file 1.

We have also performed simulation studies (data not shown) that indicate that detection power of SAIC can be further improved by correcting for normal tissue contamination using a recently developed BACOM method [18]. However, the current version of BACOM requires paired tumor-normal sampling, availability of two-channel signals, and existence of deletion CNAs. Thus, we leave the combination of SAIC and BACOM as an extension for future research.

## Conclusions

We have presented a novel approach to accurately detect significant recurrent CNAs in cancer genomes which is both statistically-principled and which, as illustrated by real examples, can be very effective at revealing SCAs within data. The concepts of CNA unit and iterative permutation are relatively simple to interpret, yet still convey considerable novel mathematical insights into data structure and bias correction.

It is worth noting that there are three novel features associated with SAIC. First, we define CNA unit to capture the intrinsic correlation structure in copy number data. Second, we perform iterative SCA-exclusive permutation to produce an unbiased null distribution. Third, we apply SAIC to real cancer copy number datasets and detect most previously reported SCAs covering well-known cancer genes.

Two important pending issues with the present algorithm are the expected significant impact of intratumor heterogeneity and normal cell contamination [18, 35, 36]. We are currently investigating applications of BACOM based normal cell correction [18] and hierarchical bi-clustering that optimize critical steps such as the selection of various thresholds and identification of subtype-specific copy number alterations.

## Appendix A

*Proof of theorem 1.*Let α’ be the significance level used in each iteration to detect SCAs in Algorithm 2. Under the truth converging null distribution, we have

for iterations $r=1,2,\dots ,\infty $ since SAIC assesses the ‘new’ SCAs at the *r* th iteration conditional on having found the ‘existing’ SCAs at the (*r*-1)th iteration.Considering

*r*th iteration,

and

*α*be the targeted FPR, we have

## Additional files

## Notes

## Declarations

### Acknowledgements

This work was supported in part by the US National Institutes of Health under Grants CA160036, CA149147, NS029525, and GM085665, and the Project Supported by Natural Science Basic Research Plan in Shaanxi Province of China (Program No. 2012JQ8027), and the Fundamental Research Funds for the Central Universities (No.K50511030002), and the Natural Science Foundation of China under Grants 61070137, 91130006, and 60933009.

## Authors’ Affiliations

## References

- Beroukhim R, Mermel CH, Porter D, Wei G, Raychaudhuri S, Donovan J, Barretina J, Boehm JS, Dobson J, Urashima M, et al: The landscape of somatic copy-number alteration across human cancers. Nature. 2010, 463 (7283): 899-905. 10.1038/nature08822.PubMed CentralView ArticlePubMed
- Leary RJ, Lin JC, Cummins J, Boca S, Wood LD, Parsons DW, Jones S, Sjoblom T, Park BH, Parsons R, et al: Integrated analysis of homozygous deletions, focal amplifications, and sequence alterations in breast and colorectal cancers. Proc Natl Acad Sci U S A. 2008, 105 (42): 16224-16229. 10.1073/pnas.0808041105.PubMed CentralView ArticlePubMed
- Beroukhim R, Getz G, Nghiemphu L, Barretina J, Hsueh T, Linhart D, Vivanco I, Lee JC, Huang JH, Alexander S, et al: Assessing the significance of chromosomal aberrations in cancer: methodology and application to glioma. Proc Natl Acad Sci U S A. 2007, 104 (50): 20007-20012. 10.1073/pnas.0710052104.PubMed CentralView ArticlePubMed
- Wood LD, Parsons DW, Jones S, Lin J, Sjoblom T, Leary RJ, Shen D, Boca SM, Barber T, Ptak J, et al: The genomic landscapes of human breast and colorectal cancers. Science. 2007, 318 (5853): 1108-1113. 10.1126/science.1145720.View ArticlePubMed
- Zhao X, Li C, Paez JG, Chin K, Janne PA, Chen TH, Girard L, Minna J, Christiani D, Leo C, et al: An integrated view of copy number and allelic alterations in the cancer genome using single nucleotide polymorphism arrays. Cancer Res. 2004, 64 (9): 3060-3071. 10.1158/0008-5472.CAN-03-3308.View ArticlePubMed
- Pinkel D, Albertson DG: Array comparative genomic hybridization and its applications in cancer. Nat Genet. 2005, 37 (Suppl): S11-S17.View ArticlePubMed
- Weir BA, Woo MS, Getz G, Perner S, Ding L, Beroukhim R, Lin WM, Province MA, Kraja A, Johnson LA, et al: Characterizing the cancer genome in lung adenocarcinoma. Nature. 2007, 450 (7171): 893-898. 10.1038/nature06358.PubMed CentralView ArticlePubMed
- Liu W, Laitinen S, Khan S, Vihinen M, Kowalski J, Yu G, Chen L, Ewing CM, Eisenberger MA, Carducci MA, et al: Copy number analysis indicates monoclonal origin of lethal metastatic prostate cancer. Nat Med. 2009, 15 (5): 559-565. 10.1038/nm.1944.PubMed CentralView ArticlePubMed
- Walter V, Nobel AB, Wright FA: DiNAMIC: a method to identify recurrent DNA copy number aberrations in tumors. Bioinformatics. 2011, 27 (5): 678-685. 10.1093/bioinformatics/btq717.PubMed CentralView ArticlePubMed
- Shah SP, Lam WL, Ng RT, Murphy KP: Modeling recurrent DNA copy number alterations in array CGH data. Bioinformatics. 2007, 23 (13): i450-i458. 10.1093/bioinformatics/btm221.View ArticlePubMed
- Rueda OM, Diaz-Uriarte R: Finding recurrent copy number alteration regions: a review of methods. Curr Bioinforma. 2010, 5: 17-View Article
- Diskin SJ, Eck T, Greshock J, Mosse YP, Naylor T, Stoeckert CJ, Weber BL, Maris JM, Grant GR: STAC: A method for testing the significance of DNA copy number aberrations across multiple array-CGH experiments. Genome Res. 2006, 16 (9): 1149-1158. 10.1101/gr.5076506.PubMed CentralView ArticlePubMed
- Ivakhno S, Tavare S: CNAnova: a new approach for finding recurrent copy number abnormalities in cancer SNP microarray data. Bioinformatics. 2010, 26 (11): 1395-1402. 10.1093/bioinformatics/btq145.View ArticlePubMed
- Mermel CH, Schumacher SE, Hill B, Meyerson ML, Beroukhim R, Getz G: GISTIC2.0 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers. Genome Biol. 2011, 12 (4): R41-10.1186/gb-2011-12-4-r41.PubMed CentralView ArticlePubMed
- Klijn C, Holstege H, de Ridder J, Liu X, Reinders M, Jonkers J, Wessels L: Identification of cancer genes using a statistical framework for multiexperiment analysis of nondiscretized array CGH data. Nucleic Acids Res. 2008, 36 (2): e13-PubMed CentralView ArticlePubMed
- Zhang Q, Ding L, Larson DE, Koboldt DC, McLellan MD, Chen K, Shi X, Kraja A, Mardis ER, Wilson RK, et al: CMDS: a population-based method for identifying recurrent DNA copy number aberrations in cancer from high-resolution data. Bioinformatics. 2010, 26 (4): 464-469. 10.1093/bioinformatics/btp708.PubMed CentralView ArticlePubMed
- Choi H, Qin ZS, Ghosh D: A double-layered mixture model for the joint analysis of DNA copy number and gene expression data. J Comput Biol. 2010, 17 (2): 121-137. 10.1089/cmb.2009.0019.PubMed CentralView ArticlePubMed
- Yu G, Zhang B, Bova GS, Xu J, Shih IM, Wang Y: BACOM: in silico detection of genomic deletion types and correction of normal cell contamination in copy number data. Bioinformatics. 2011, 27 (11): 1473-1480. 10.1093/bioinformatics/btr183.PubMed CentralView ArticlePubMed
- Fu JC, Johnson BC, Chang Y-M: Approximating the extreme right-hand tail probability for the distribution of the number of patterns in a sequence of multi-state trials. Journal of Statistical Planning and Inference. 2011, 142 (2): 473-480.View Article
- Li C, Wong WH: Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. Proc Natl Acad Sci U S A. 2001, 98 (1): 31-36.PubMed CentralView ArticlePubMed
- Lin M, Wei LJ, Sellers WR, Lieberfarb M, Wong WH, Li C: dChipSNP: significance curve and clustering of SNP-array-based loss-of-heterozygosity data. Bioinformatics. 2004, 20 (8): 1233-1240. 10.1093/bioinformatics/bth069.View ArticlePubMed
- Venkatraman ES, Olshen AB: A faster circular binary segmentation algorithm for the analysis of array CGH data. Bioinformatics. 2007, 23 (6): 657-663. 10.1093/bioinformatics/btl646.View ArticlePubMed
- Olshen AB, Venkatraman ES, Lucito R, Wigler M: Circular binary segmentation for the analysis of array-based DNA copy number data. Biostatistics. 2004, 5 (4): 557-572. 10.1093/biostatistics/kxh008.View ArticlePubMed
- Willenbrock H, Fridlyand J: A comparison study: applying segmentation to array CGH data for downstream analyses. Bioinformatics. 2005, 21 (22): 4084-4091. 10.1093/bioinformatics/bti677.View ArticlePubMed
- Chen L, Yu G, Langefeld CD, Miller DJ, Guy RT, Raghuram J, Yuan X, Herrington DM, Wang Y: Comparative analysis of methods for detecting interacting loci. BMC Genomics. 2011, 12: 344-10.1186/1471-2164-12-344.PubMed CentralView ArticlePubMed
- Kuo KT, Guan B, Feng Y, Mao TL, Chen X, Jinawath N, Wang Y, Kurman RJ, Shih Ie M, Wang TL: Analysis of DNA copy number alterations in ovarian serous tumors identifies new molecular genetic changes in low-grade and high-grade carcinomas. Cancer Res. 2009, 69 (9): 4036-4042. 10.1158/0008-5472.CAN-08-3913.PubMed CentralView ArticlePubMed
- Kuo KT, Mao TL, Chen X, Feng Y, Nakayama K, Wang Y, Glas R, Ma MJ, Kurman RJ, Shih Ie M, et al: DNA copy numbers profiles in affinity-purified ovarian clear cell carcinoma. Clin Cancer Res. 2010, 16 (7): 1997-2008. 10.1158/1078-0432.CCR-09-2105.PubMed CentralView ArticlePubMed
- Feng Y, Yu G, Wang T-L, Shih I-M, Wang Y: Analyzing DNA copy number changes using fused margin regression. Intl J of Functional Informatics and Personalized Medicine. 2010, 3 (1): 3-15. 10.1504/IJFIPM.2010.033242.View Article
- Bast RC, Hennessy B, Mills GB: The biology of ovarian cancer: new opportunities for translation. Nat Rev Cancer. 2009, 9 (6): 415-428. 10.1038/nrc2644.PubMed CentralView ArticlePubMed
- Integrated genomic analyses of ovarian carcinoma. Nature. 2011, 474 (7353): 609-615. 10.1038/nature10166.
- Lynch TJ, Bell DW, Sordella R, Gurubhagavatula S, Okimoto RA, Brannigan BW, Harris PL, Haserlat SM, Supko JG, Haluska FG, et al: Activating mutations in the epidermal growth factor receptor underlying responsiveness of non-small-cell lung cancer to gefitinib. N Engl J Med. 2004, 350 (21): 2129-2139. 10.1056/NEJMoa040938.View ArticlePubMed
- Kotliarov Y, Steed ME, Christopher N, Walling J, Su Q, Center A, Heiss J, Rosenblum M, Mikkelsen T, Zenklusen JC, et al: High-resolution global genomic survey of 178 gliomas reveals novel regions of copy number alteration and allelic imbalances. Cancer Res. 2006, 66 (19): 9428-9436. 10.1158/0008-5472.CAN-06-1691.PubMed CentralView ArticlePubMed
- Westfall PH, Young SS: Resampling-based multiple testing : examples and methods for P-value adjustment. 1993, New York, Wiley
- Shaffer JP: Multiple hypothesis testing. Annu Rev Psychol. 1995, 46: 24-View Article
- Shibata D: Cancer. Heterogeneity and tumor history. Science. 2012, 336 (6079): 304-305. 10.1126/science.1222361.View ArticlePubMed
- Marusyk A, Almendro V, Polyak K: Intra-tumour heterogeneity: a looking glass for cancer?. Nat Rev Cancer. 2012, 12 (5): 323-334. 10.1038/nrc3261.View ArticlePubMed

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.