Skip to main content

JISTIC: Identification of Significant Targets in Cancer

Abstract

Background

Cancer is caused through a multistep process, in which a succession of genetic changes, each conferring a competitive advantage for growth and proliferation, leads to the progressive conversion of normal human cells into malignant cancer cells. Interrogation of cancer genomes holds the promise of understanding this process, thus revolutionizing cancer research and treatment. As datasets measuring copy number aberrations in tumors accumulate, a major challenge has become to distinguish between those mutations that drive the cancer versus those passenger mutations that have no effect.

Results

We present JISTIC, a tool for analyzing datasets of genome-wide copy number variation to identify driver aberrations in cancer. JISTIC is an improvement over the widely used GISTIC algorithm. We compared the performance of JISTIC versus GISTIC on a dataset of glioblastoma copy number variation, JISTIC finds 173 significant regions, whereas GISTIC only finds 103 significant regions. Importantly, the additional regions detected by JISTIC are enriched for oncogenes and genes involved in cell-cycle and proliferation.

Conclusions

JISTIC is an easy-to-install platform independent implementation of GISTIC that outperforms the original algorithm detecting more relevant candidate genes and regions. The software and documentation are freely available and can be found at: http://www.c2b2.columbia.edu/danapeerlab/html/software.html

Background

A comprehensive study of the genomic alterations that occur in cancer is vital for understanding this disease. Technological advances have made it possible to detect chromosomal regions of amplifications and deletions genome-wide and at high resolution. Datasets measuring such aberrations in patient tumors are accumulating at a staggering rate for multiple types of cancer [1–4]

However, tumors harbor a great number of copy number alterations and it is difficult to distinguish between driver aberrations (functional changes vital for cancer progression) and passenger aberrations (random and with no selective advantage). Thus, the distinction between driver and passenger mutations has become one of the key challenges in cancer genomics. A very successful algorithm to address this is "Genomic Identification of Significant Targets in Cancer (GISTIC)" [1], that identifies aberrant regions more likely to drive cancer pathogenesis. GISTIC calculates the background rate of random chromosomal aberrations and identifies those regions that are aberrant more often than would be expected by chance, with greater weight given to high amplitude events that are less likely to represent random aberrations. There are other algorithms that tackle this task such as GLAD [5], RAE [6] and STAC [6]. However, GISTIC is unique in its ability to combine magnitude and frequency of the alteration into a statistical score. This algorithm has been successfully applied to various datasets [2, 4, 7, 8] and the approach is becoming widespread.

GISTIC identifies those regions of the genome that are aberrant more often than would be expected by chance. While successful in most scenarios, GISTIC has trouble identifying the relevant sub-region when a very large region is amplified or deleted. Such large chromosomal aberrations frequently occur in cancer and this leaves the user with two less than optimal options - consider only a single peak within the region, or consider an entire chromosomal arm. However, we have observed that in many cases there are other small sub-regions for which the aberration is significantly stronger than in the rest of the large region. Moreover, these regions often contain known oncogenes.

To address this issue, we developed JISTIC, a tool that implements all of GISTIC's capabilities, with an additional new variant of the algorithm capable of detecting multiple significant sub-regions within large aberrant regions.

Implementation

JISTIC is based on the GISTIC algorithm [1]. JISTIC implements the previously published variants of GISTIC (standard, focal and arm-peel-off) and can also deal with LOH in the same way that the original algorithm does. More detailed information on GISTIC can be found in the Supporting Information of [1].

In brief, GISTIC calculates a statistic (G-score) which represents the strength of the aberration for each marker. The G-score for a marker m is the summation of the copy number across all samples. For this summation, the samples in which the copy number is less significant than an empirical aberration threshold (ΘAMP for amplification, ΘDEL for deletion) is set to zero. Therefore, the G-score in the case of amplifications is:

Where I(x) is the indicator function and CN(m, i) the copy number for marker m and sample i. The score is defined similarly for deletions, taking into account the change in sign. While standard GISTIC uses a fixed aberration threshold for each type of aberration, focal GISTIC uses sample-specific high-level thresholds. This threshold is set for each sample by adding the standard threshold to the maximum (minimum for deletions) of medians observed for each chromosome arm.

GISTIC assesses the significance of the G-score for each marker by comparing this score to results expected by chance using genome-wide permutation testing. This significance is then corrected using False Discovery Rate (FDR) using Benjamini and Hochberg method [9], and a q-value is obtained. All regions with a q-value below a threshold (0.25 in previously published articles) are deemed significant. For large aberrations, the sub-region with a minimal q-value is identified as a peak driver region.

To identify independent peaks within a region and discard spurious peaks, GISTIC uses a peel-off algorithm. The algorithm iteratively finds the most significant peak and then, for each sample, if the peak is included in the aberration, it sets to zero all markers consecutive with the peak, thus removing all aberrations overlapping the peak. Once a peak has been detected, the peak region is extended by leaving each sample out in turn, and recalculating the peak boundaries.

Typical to cancer genomes is the presence of broad copy number aberrations, defined as aberrant regions at least as large as one half of a chromosomal arm [1]. Peel-off on broad regions using standard GISTIC usually results in identifying only a single peak. GISTIC also has a focal variant of the algorithm which can potentially capture more aberrations, but misses many peaks in practice.

The crux of the matter is to distinguish between genuine multiple peaks and a single peak within noisy fluctuations of microarray intensities. Focal peel-off is designed to deal with this issue, but since the focal threshold is defined according to the highest broad aberration genome-wide, this variant not only requires the aberration to have a strong focal behavior, but also depends on the strength of other broad aberrations across the genome. In Figure 1 we demonstrate how different thresholds (determined by broad aberrations in other chromosomes not shown in the figure) can lead to either detection or failure to detect a second peak. In this example, standard GISTIC only captures a single peak and misses a second peak that seems equally significant, with just a slight difference in G-score. In this type of situations standard GISTIC is doomed to fail. The detection of the second peak with focal GISTIC highly depends on the focal threshold. This threshold is set for each sample by the highest broad aberration genome-wide and there are cases, such as C in Figure 1, in which the focal events that define the aberration are masked by broad aberrations in other chromosome (not shown in Figure). In the next section we will see that this problem is prevalent in real datasets.

Figure 1
figure 1

GISTIC applied to example. Represented is a toy chromosome for 5 samples. The X axis corresponds to consecutive markers across the chromosome and each bar represents the G-score contribution of that marker. Two possible cut-offs determined by focal GISTIC are represented by dotted blue and red lines, these cut-offs are determined based on broad regions in other chromosomes (which are not shown in the Figure for the sake of clarity). Three results from GISTIC are contemplated: standard (A), focal based on the blue threshold (B) and focal based on the red threshold (C). Reported peaks are displayed as green bars at the top, Focal GISTIC's ability to capture additional peaks is dependent on the threshold determined by aberrations in other chromosomes.

For JISTIC, we developed a new variant called limited peel-off, designed to tackle the problem posed in Figure 1. While GISTIC sets to zero any aberration containing the peak, our algorithm "peels off" only part of the aberration. Limited peel-off uses the global G-score to determine the extent of "peel-off. Given a single aberration, we expect G(m) would decrease as we get further from it until it reaches the noise level. The idea behind limited peel-off is to decompose the G-score for a marker in two parts, one that represents what remains from the peak that we are peeling off (Gr(m)) and another that depicts contribution independent of the peak (Gn(m)). We use a threshold value (Gthres) to estimate whether Gn(m) contains only noise or the aberration is likely due to an independent peak.

Denote G(m, i) as the G-score contribution of sample i in marker m, which can be defined as:

For each sample i, the algorithm iteratively calculates on both sides of a newly discovered peak the magnitude G r (m, i). For the right side of the peak:

the left side of the peak is the symmetric equation:

Gr(m, i) represents the amount of aberration that remains from the primary peak in marker m. JISTIC considers any reduction in the aberration that is consistent for at least a minimal number of markers, s. This parameter was introduced after observing that transient fluctuations in copy number result in a greater number of false positives. For s = 0 a total of 147 peaks are detected with limited peel-off. However, many of those peaks are located next to each other in the genome and the biological assessment suggested that those were spurious peaks. The number of reported peaks decreases when the parameter s increases, eliminating most of those spurious peaks.

In order to determine whether a new peak exists, we need the complementary component of Gr(m, i), Gn(m, i), measuring the component of the aberration that is independent from the primary peak, for each sample:

The condition that JISTIC checks in order to abort the peel-off at a marker m is

The results critically depend on the threshold Gthres. If Gthres is too high, the algorithm will obtain exactly the same result than the standard variant of GISTIC. If the value is too low, the algorithm will obtain many spurious peaks. For consistency, we use the G-score value that corresponds to the threshold q-value when running the focal variant, thus using the same threshold both for a peak on baseline copy number and a sub-peak within a broader region.

Figure 2 shows how limited peel-off successfully detects the second peak in the example from Figure 1. Contrary to the focal variant, this success is independent of broad regions in other chromosomes. It is important to note that limited peel-off fully replaces the standard variant, as any peak captured by the latter, will also be captured by the former. On the other hand, focal and limited peel-off can be considered complementary. While limited peel-off looks at the average of the aberration across samples, focal concentrates on the most acute aberrations for each sample.

Figure 2
figure 2

Limited peel-off applied to example. Limited peel-off for the same example shown in Figure 1. The x axis corresponds to consecutive markers across the chromosome while each bar represents the G-score contribution of that marker. Reported peaks are displayed as green bars at the top. Limited peel-off successfully detects two peaks in the example, independently of aberrations in other chromosomes. The figure also illustrates how thresholds in limited peel-off are applied. In this example s = 1. (A) The G-score contribution for each sample is decomposed into the G-score for the primary peak (green), Gr (white) and Gn (red). (B) Each bar illustrates the total Gn of the marker. The threshold Gthres represents the cut-off that determines whether the peel-off should be aborted. In this example the peel-off is aborted in marker 14, allowing the detection of an additional peak.

Results

We tested JISTIC on 178 glioblastoma samples that had been previously used to test the latest version of GISTIC [10]. This dataset was generated by Harvard Medical School using Agilent 244 K Arrays.

The parameter s was tuned using the histogram of segments sizes to obtain candidate values and the optimal value was chosen based on the distribution of minimum distance between peaks. In the test dataset, the histogram of segment size showed that the number of segments decreases with the segment size, as it was expected. We observed a stronger decrease in the number of segments for 3 different segment sizes (3, 8 and 12); those sizes were considered as candidate values for s. The distribution of minimum distance between peaks was used to estimate the number of spurious peaks for each of those values. This distribution shows a considerable change when s was increased to 8, as the number of adjacent peaks decreased. The distribution did not show any significant change when s was increased to 12 and we adopted a conservative approach by setting s = 10, which corresponds to 126.87 kb on average. Limited peel-off detected 124 peaks using s = 10.

Table 1 shows the peaks and genes obtained with JISTIC in the three variants: standard, focal and limited peel-off. Note that standard and focal are equivalent to implementations of the GISTIC algorithm [10]. When looking at peaks by chromosome (Table 2), we observe that most of the novel peaks appear in chromosomes with a broad region for which standard detects only a single peak and focal detects one or no peaks.

Table 1 Summary of peaks and genes for the three variants.
Table 2 Peaks by chromosome.

Figure 3 demonstrates an example of such a broad peak in chromosome 19. Standard (A) captures only a single peak and focal (B) fails to capture any aberration. On the other hand, limited peel-off (C) successfully captures multiple relevant regions. The ultimate assessment is whether the additional peaks identified are indeed drivers. To evaluate this we assessed whether the new peaks in chromosome 19 contained solid candidates for oncogenes. Standard GISTIC identified a single peak containing the gene (ZNF626), whose role is unknown. In comparison, limited peel-off detected 7 additional regions, where each includes at least one oncogene and multiple genes related to cell replication. Examples of oncogenes include AP1 M2 [11], Cyclin E1 [12], LYPD3 [13], see Table 3 for list of potential oncogenes in each peak.

Table 3 Novel amplification peaks in chrom. 19 reported by limited peel-off.
Figure 3
figure 3

Amplification aberrations for chromosome 19. Represented are standard (A), focal (B) and limited peel-off (C) results for amplifications in the whole chromosome 19. The x axis corresponds to consecutive markers across the chromosome and the y axis corresponds to the q-value in logarithmic scale. The significance threshold for aberrations (q-value = 0.25) is represented by a blue line. Reported peak regions are illustrated in green. Limited peel-off detects 6 novel peaks.

Figure 4 shows another example, a deletion of the entire chromosome 22. Standard finds a single peak and focal finds no peaks, whereas limited peel-off finds 3 additional peaks. The peak identified by standard includes FBLN1, a gene that has been previously selected as candidate tumor suppressor [14]. Each additional peak identified by limited peek-off has at least one tumor suppressor, such as SMARCB1 [15], HIRA [16] and LARGE [17], see Table 4 for list of potential tumor suppressors in each peak.

Table 4 Novel deletion peaks in chrom. 22 reported by limited peel-off.
Figure 4
figure 4

Deletion aberrations for chromosome 22. Represented are standard (A), focal (B) and limited peel-off (C) results for deletions in the whole chromosome 22. The x axis corresponds to consecutive markers across the chromosome and the y axis corresponds to the q-value in logarithmic scale. Note that only markers for which data is available are plotted. The significance threshold for aberrations (q-value = 0.25) is represented by a blue line. Reported peak regions are illustrated in green. Limited peel-off detects 3 novel peaks.

For a systematic evaluation we designed an automated statistical test to study the limited peel-off's specificity and establish that the increased number of peaks detected by limited peel-off does not increase the false positive rate. Our permutation based test estimates the number of candidate cancer driver genes expected in random regions matching in size with the regions detected by JISTIC. We compiled a list of 1880 likely driver genes based on two different sources:

  • Genes listed in the Sanger Cancer Gene Census [18].

  • Genes annotated in GO [19] for at least one of the following processes: cell cycle, cell proliferation, cell death and neurogenesis.

We generated 10,000 random outputs with similar characteristics to the real output. A random output is generated by iterating over all detected peaks and shifting the peak's location to a random position in the genome. For each output, we count the number of regions containing at least one gene in the list of driver candidates above. This provides us with a distribution for the number of potential driver genes expected in the null model.

To compare focal and limited peel-off, we tested three sets of peaks: the peaks that are exclusive to each of the two methods and the peaks that are common to both. Note that all the peaks obtained by standard GISTIC are also obtained by limited peel-off. The results of the 10,000 iterations for each of the three sets are shown in Table 5.

Table 5 Results for permutation test for focal and limited peel-off peaks.

As expected, the p-value obtained for the overlapping peaks is more significant than the p-value for the exclusive peaks for each method.

For the 78 peaks obtained exclusively by limited peel-off, 25 contained likely driver genes, only 23/10,000 permutations reached the value obtained by the true output, giving a p-value of 0.0023. In comparison, for exclusive peaks for the focal variant, 9 of the 49 peaks contained likely driver genes. In 73/10,000 permutations the test obtained as many peaks with likely driver genes, giving a p-value of 0.0073.

Conclusions

The analysis performed on limited peel-off against standard and focal GISTIC demonstrates the superiority of limited peel-off to achieve both better specificity and dramatically increase recall by obtaining a large number of novel peaks.

In conclusion, JISTIC is a significantly improved algorithm to distinguish between driver and passenger copy number aberrations in cancer genomes. Importantly, it detects a significant number of additional driver regions while maintaining a similar false positive rate. We conclude that both limited peel-off and focal GISTIC should be used together as they provide complementary and significant results.

The tool is implemented in Java, has been tested on Linux and Windows. It can be downloaded from http://www.c2b2.columbia.edu/danapeerlab/html/software.html. It does not have dependencies to external libraries and can be downloaded as a single Java JAR file. The execution time for the glioblastoma dataset on a standard desktop computer (Intel Xeon W3505 @ 2.53 GHz, 3GB of RAM) was 8 minutes for all variants.

Matlab scripts are provided in order to visualize the output and obtain different statistics. The tool can also convert to the format used by the open-source visualization tool IGV [20], used to display cancer genomic data using a user-friendly interface, demonstrated in Figure 5.

Figure 5
figure 5

IGV display of JISTIC results. Segmented copy number and JISTIC results can be simultaneously displayed using the visualization tool IGV [20]. The main window displays the segmented copy number data used as JISTIC input. Two tracks at the bottom show JISTIC's output (G-score and peak regions respectively).

Availability and Requirements

Conflict of interests

The authors declare that they have no competing interests.

References

  1. 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 USA 2007, 104(50):20007–20012. 10.1073/pnas.0710052104

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. 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: 893–898. 10.1038/nature06358

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Pollack JR, Sørlie T, Perou CM, Rees CA, Jeffrey SS, Lonning PE, Tibshirani R, Botstein D, Børresen-Dale A-L, Brown PO: Microarray analysis reveals a major direct role of DNA copy number alteration in the transcriptional program of human breast tumors. Proceedings of the National Academy of Sciences of the United States of America 2002, 99(20):12963–12968. 10.1073/pnas.162471999

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Northcott PA, Nakahara Y, Wu X, Feuk L, Ellison DW, Croul S, Mack S, Kongkham PN, Peacock J, Dubuc A, et al.: Multiple recurrent genetic events converge on control of histone lysine methylation in medulloblastoma. Nat Genet 2009, 41(4):465–472. 10.1038/ng.336

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Hupe P, Stransky N, Thiery J-P, Radvanyi F, Barillot E: Analysis of array CGH data: from signal ratio to gain and loss of DNA regions. Bioinformatics 2004, 20(18):3413–3422. 10.1093/bioinformatics/bth418

    Article  CAS  PubMed  Google Scholar 

  6. Taylor BS, Barretina J, Socci ND, DeCarolis P, Ladanyi M, Meyerson M, Singer S, Sander C: Functional Copy-Number Alterations in Cancer. PLoS ONE 2008, 3(9):e3179. 10.1371/journal.pone.0003179

    Article  PubMed  PubMed Central  Google Scholar 

  7. Firestein R, Bass AJ, Kim SY, Dunn IF, Silver SJ, Guney I, Freed E, Ligon AH, Vena N, Ogino S, et al.: CDK8 is a colorectal cancer oncogene that regulates [bgr]-catenin activity. Nature 2008, 455(7212):547–551. 10.1038/nature07179

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Salvesen HB, Carter SL, Mannelqvist M, Dutt A, Getz G, Stefansson IM, Raeder MB, Sos ML, Engelsen IB, Trovik J, et al.: Integrated genomic profiling of endometrial carcinoma associates aggressive tumors with indicators of PI3 kinase activation. Proceedings of the National Academy of Sciences 2009, 106(12):4834–4839. 10.1073/pnas.0806514106

    Article  CAS  Google Scholar 

  9. Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B (Methodological) 1995, 57(1):289–300.

    Google Scholar 

  10. Cancer Genome Atlas Research Network: Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature 2008, 455(7216):1061–1068. 10.1038/nature07385

    Article  Google Scholar 

  11. Rohan S, Tu JJ, Kao J, Mukherjee P, Campagne F, Zhou XK, Hyjek E, Alonso MA, Chen YT: Gene expression profiling separates chromophobe renal cell carcinoma from oncocytoma and identifies vesicular transport and cell junction proteins as differentially expressed genes. Clin Cancer Res 2006, 12(23):6937–6945. 10.1158/1078-0432.CCR-06-1268

    Article  CAS  PubMed  Google Scholar 

  12. Jeon HM, Jin X, Lee JS, Oh SY, Sohn YW, Park HJ, Joo KM, Park WY, Nam DH, DePinho RA, et al.: Inhibitor of differentiation 4 drives brain tumor-initiating cell genesis through cyclin E and notch signaling. Genes Dev 2008, 22(15):2028–2033. 10.1101/gad.1668708

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Jacobsen B, Ploug M: The urokinase receptor and its structural homologue C4.4A in human cancer: expression, prognosis and pharmacological inhibition. Curr Med Chem 2008, 15(25):2559–2573. 10.2174/092986708785909012

    Article  CAS  PubMed  Google Scholar 

  14. Cheng YY, Jin H, Liu X, Siu JMT, Wong YP, Ng EKO, Yu J, Leung Wk, Sung JJY, Chan FKL: Fibulin 1 is downregulated through promoter hypermethylation in gastric cancer. Br J Cancer 2008, 99(12):2083–2087. 10.1038/sj.bjc.6604760

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Bug S, Durig J, Oyen F, Klein-Hitpass L, Martin-Subero JI, Harder L, Baudis M, Arnold N, Kordes U, Duhrsen U, et al.: Recurrent loss, but lack of mutations, of the SMARCB1 tumor suppressor gene in T-cell prolymphocytic leukemia with TCL1A-TCRAD juxtaposition. Cancer Genet Cytogenet 2009, 192(1):44–47. 10.1016/j.cancergencyto.2009.03.001

    Article  CAS  PubMed  Google Scholar 

  16. Hall C, Nelson DM, Ye X, Baker K, DeCaprio JA, Seeholzer S, Lipinski M, Adams PD: HIRA, the Human Homologue of Yeast Hir1p and Hir2p, Is a Novel Cyclin-cdk2 Substrate Whose Expression Blocks S-Phase Progression. Mol Cell Biol 2001, 21(5):1854–1865. 10.1128/MCB.21.5.1854-1865.2001

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. de Bernabe DB-V, Inamori K-i, Yoshida-Moriguchi T, Weydert CJ, Harper HA, Willer T, Henry MD, Campbell KP: Loss of {alpha}-Dystroglycan Laminin Binding in Epithelium-derived Cancers Is Caused by Silencing of LARGE. J Biol Chem 2009, 284(17):11279–11284. 10.1074/jbc.C900007200

    Article  PubMed  PubMed Central  Google Scholar 

  18. Futreal PA, Coin L, Marshall M, Down T, Hubbard T, Wooster R, Rahman N, Stratton MR: A census of human cancer genes. Nat Rev Cancer 2004, 4(3):177–183. 10.1038/nrc1299

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al.: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 2000, 25(1):25–29. 10.1038/75556

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Gey S, Nedelec E: Model selection for CART regression trees. Information Theory, IEEE Transactions on 2005, 51(2):658–670. 10.1109/TIT.2004.840903

    Article  Google Scholar 

  21. Bushell M, Stoneley M, Kong YW, Hamilton TL, Spriggs KA, Dobbyn HC, Qin X, Sarnow P, Willis AE: Polypyrimidine tract binding protein regulates IRES-mediated gene expression during apoptosis. Mol Cell 2006, 23(3):401–412. 10.1016/j.molcel.2006.06.012

    Article  CAS  PubMed  Google Scholar 

  22. Su J, Chen X, Kanekura T: A CD147-targeting siRNA inhibits the proliferation, invasiveness, and VEGF production of human malignant melanoma cells by down-regulating glycolysis. Cancer Lett 2009, 273(1):140–147. 10.1016/j.canlet.2008.07.034

    Article  CAS  PubMed  Google Scholar 

  23. Glinsky GV, Glinskii AB, Stephenson AJ, Hoffman RM, Gerald WL: Gene expression profiling predicts clinical outcome of prostate cancer. J Clin Invest 2004, 113(6):913–923.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Yam JW, Jin DY, So CW, Chan LC: Identification and characterization of EBP, a novel EEN binding protein that inhibits Ras signaling and is recruited into the nucleus by the MLL-EEN fusion protein. Blood 2004, 103(4):1445–1453. 10.1182/blood-2003-07-2452

    Article  CAS  PubMed  Google Scholar 

  25. Guan D, Altan-Bonnet N, Parrott AM, Arrigo CJ, Li Q, Khaleduzzaman M, Li H, Lee C-G, Pe'ery T, Mathews MB: Nuclear Factor 45 (NF45) Is a Regulatory Subunit of Complexes with NF90/110 Involved in Mitotic Control. Mol Cell Biol 2008, 28(14):4629–4641. 10.1128/MCB.00120-08

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Suzuki C, Murakumo Y, Kawase Y, Sato T, Morinaga T, Fukuda N, Enomoto A, Ichihara M, Takahashi M: A novel GDNF-inducible gene, BMZF3, encodes a transcriptional repressor associated with KAP-1. Biochemical and Biophysical Research Communications 2008, 366(1):226–232. 10.1016/j.bbrc.2007.11.118

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank Gad Getz for help with implementing the GISTIC algorithm and useful discussions.

Funding: This research was supported by the National Institutes of Health Roadmap Initiative, NIH Director's New Innovator Award Program through grant number 1-DP2-OD002414-01 and National Centers for Biomedical Computing Grant 1U54CA121852-01A1. DP holds a Career Award at the Scientific Interface from the Burroughs Welcome Fund.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Felix Sanchez-Garcia.

Additional information

Authors' contributions

DP, FSG and UDA designed research. FSG and UDA developed the limited-peel off algorithm EM and FSG implemented the JAVA code. FSG and UDA analyzed and evaluated the algorithm. FSG and DP wrote the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

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.

Reprints and permissions

About this article

Cite this article

Sanchez-Garcia, F., Akavia, U.D., Mozes, E. et al. JISTIC: Identification of Significant Targets in Cancer. BMC Bioinformatics 11, 189 (2010). https://doi.org/10.1186/1471-2105-11-189

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-11-189

Keywords