Sequenza: Allele-specific Copy Number and Mutation Profiles from Tumor Sequencing Data

Overview

Journal Ann Oncol

Publisher Elsevier

Specialty Oncology

Date 2014 Oct 17

PMID 25319062

Citations 409

Authors

F Favero

T Joshi

A M Marquard

N J Birkbak

M Krzystanek

Q Li

Z Szallasi

A C Eklund

Affiliations

Soon will be listed here.

Abstract

Background: Exome or whole-genome deep sequencing of tumor DNA along with paired normal DNA can potentially provide a detailed picture of the somatic mutations that characterize the tumor. However, analysis of such sequence data can be complicated by the presence of normal cells in the tumor specimen, by intratumor heterogeneity, and by the sheer size of the raw data. In particular, determination of copy number variations from exome sequencing data alone has proven difficult; thus, single nucleotide polymorphism (SNP) arrays have often been used for this task. Recently, algorithms to estimate absolute, but not allele-specific, copy number profiles from tumor sequencing data have been described.

Materials And Methods: We developed Sequenza, a software package that uses paired tumor-normal DNA sequencing data to estimate tumor cellularity and ploidy, and to calculate allele-specific copy number profiles and mutation profiles. We applied Sequenza, as well as two previously published algorithms, to exome sequence data from 30 tumors from The Cancer Genome Atlas. We assessed the performance of these algorithms by comparing their results with those generated using matched SNP arrays and processed by the allele-specific copy number analysis of tumors (ASCAT) algorithm.

Results: Comparison between Sequenza/exome and SNP/ASCAT revealed strong correlation in cellularity (Pearson's r = 0.90) and ploidy estimates (r = 0.42, or r = 0.94 after manual inspecting alternative solutions). This performance was noticeably superior to previously published algorithms. In addition, in artificial data simulating normal-tumor admixtures, Sequenza detected the correct ploidy in samples with tumor content as low as 30%.

Conclusions: The agreement between Sequenza and SNP array-based copy number profiles suggests that exome sequencing alone is sufficient not only for identifying small scale mutations but also for estimating cellularity and inferring DNA copy number aberrations.

Citing Articles

Evolutionary fingerprints of epithelial-to-mesenchymal transition.

Perelli L, Zhang L, Mangiameli S, Giannese F, Mahadevan K, Peng F Nature. 2025; .

PMID: 40044861 DOI: 10.1038/s41586-025-08671-2.

Benchmarking copy number aberrations inference tools using single-cell multi-omics datasets.

Song M, Ma S, Wang G, Wang Y, Yang Z, Xie B Brief Bioinform. 2025; 26(2).

PMID: 40037644 PMC: 11879432. DOI: 10.1093/bib/bbaf076.

Distant origin of glioblastoma recurrence: neural stem cells in the subventricular zone serve as a source of tumor reconstruction after primary resection.

Li X, Kim H, Yoo J, Lee Y, Nam C, Park J Mol Cancer. 2025; 24(1):64.

PMID: 40033380 PMC: 11877783. DOI: 10.1186/s12943-025-02273-2.

Genomic alterations and transcriptional phenotypes in circulating free DNA and matched metastatic tumor.

Takahashi N, Pongor L, Agrawal S, Shtumpf M, Gurjar A, Rajapakse V Genome Med. 2025; 17(1):15.

PMID: 40001151 PMC: 11863907. DOI: 10.1186/s13073-025-01438-4.

Multi-omics analyses reveal biological and clinical insights in recurrent stage I non-small cell lung cancer.

Wang C, Li J, Chen J, Wang Z, Zhu G, Song L Nat Commun. 2025; 16(1):1477.

PMID: 39929832 PMC: 11811181. DOI: 10.1038/s41467-024-55068-2.

References

Ding L, Raphael B, Chen F, Wendl M . Advances for studying clonal evolution in cancer. Cancer Lett. 2013; 340(2):212-9. PMC: 3783624. DOI: 10.1016/j.canlet.2012.12.028. View

Carter S, Cibulskis K, Helman E, McKenna A, Shen H, Zack T . Absolute quantification of somatic DNA alterations in human cancer. Nat Biotechnol. 2012; 30(5):413-21. PMC: 4383288. DOI: 10.1038/nbt.2203. View

Nilsen G, Liestol K, Van Loo P, Vollan H, Eide M, Rueda O . Copynumber: Efficient algorithms for single- and multi-track copy number segmentation. BMC Genomics. 2013; 13:591. PMC: 3582591. DOI: 10.1186/1471-2164-13-591. View

. Integrated genomic analyses of ovarian carcinoma. Nature. 2011; 474(7353):609-15. PMC: 3163504. DOI: 10.1038/nature10166. View

Gerlinger M, Rowan A, Horswell S, Math M, Larkin J, Endesfelder D . Intratumor heterogeneity and branched evolution revealed by multiregion sequencing. N Engl J Med. 2012; 366(10):883-892. PMC: 4878653. DOI: 10.1056/NEJMoa1113205. View

Bengtsson H, Wirapati P, Speed T . A single-array preprocessing method for estimating full-resolution raw copy numbers from all Affymetrix genotyping arrays including GenomeWideSNP 5 & 6. Bioinformatics. 2009; 25(17):2149-56. PMC: 2734319. DOI: 10.1093/bioinformatics/btp371. View

Artez A, Bernabe R, Eerola I, Hemsley F, Jennings J, Kerr D . International network of cancer genome projects. Nature. 2010; 464(7291):993-8. PMC: 2902243. DOI: 10.1038/nature08987. View

Roschke A, Tonon G, Gehlhaus K, McTyre N, Bussey K, Lababidi S . Karyotypic complexity of the NCI-60 drug-screening panel. Cancer Res. 2003; 63(24):8634-47. View

. Comprehensive molecular characterization of clear cell renal cell carcinoma. Nature. 2013; 499(7456):43-9. PMC: 3771322. DOI: 10.1038/nature12222. View

10.

Shah S, Roth A, Goya R, Oloumi A, Ha G, Zhao Y . The clonal and mutational evolution spectrum of primary triple-negative breast cancers. Nature. 2012; 486(7403):395-9. PMC: 3863681. DOI: 10.1038/nature10933. View

11.

Popova T, Manie E, Stoppa-Lyonnet D, Rigaill G, Barillot E, Stern M . Genome Alteration Print (GAP): a tool to visualize and mine complex cancer genomic profiles obtained by SNP arrays. Genome Biol. 2009; 10(11):R128. PMC: 2810663. DOI: 10.1186/gb-2009-10-11-r128. View

12.

Ha G, Roth A, Lai D, Bashashati A, Ding J, Goya R . Integrative analysis of genome-wide loss of heterozygosity and monoallelic expression at nucleotide resolution reveals disrupted pathways in triple-negative breast cancer. Genome Res. 2012; 22(10):1995-2007. PMC: 3460194. DOI: 10.1101/gr.137570.112. View

13.

Van Loo P, Nordgard S, Lingjaerde O, Russnes H, Rye I, Sun W . Allele-specific copy number analysis of tumors. Proc Natl Acad Sci U S A. 2010; 107(39):16910-5. PMC: 2947907. DOI: 10.1073/pnas.1009843107. View

14.

Cibulskis K, Lawrence M, Carter S, Sivachenko A, Jaffe D, Sougnez C . Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat Biotechnol. 2013; 31(3):213-9. PMC: 3833702. DOI: 10.1038/nbt.2514. View

15.

Abaan O, Polley E, Davis S, Zhu Y, Bilke S, Walker R . The exomes of the NCI-60 panel: a genomic resource for cancer biology and systems pharmacology. Cancer Res. 2013; 73(14):4372-82. PMC: 4893961. DOI: 10.1158/0008-5472.CAN-12-3342. View

16.

Larson N, Fridley B . PurBayes: estimating tumor cellularity and subclonality in next-generation sequencing data. Bioinformatics. 2013; 29(15):1888-9. PMC: 3712213. DOI: 10.1093/bioinformatics/btt293. View

17.

Horswell S, Matthews N, Swanton C . Cancer heterogeneity and "the struggle for existence": diagnostic and analytical challenges. Cancer Lett. 2012; 340(2):220-6. DOI: 10.1016/j.canlet.2012.10.031. View

18.

Su X, Zhang L, Zhang J, Meric-Bernstam F, Weinstein J . PurityEst: estimating purity of human tumor samples using next-generation sequencing data. Bioinformatics. 2012; 28(17):2265-6. PMC: 3426843. DOI: 10.1093/bioinformatics/bts365. View

19.

Koboldt D, Zhang Q, Larson D, Shen D, McLellan M, Lin L . VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing. Genome Res. 2012; 22(3):568-76. PMC: 3290792. DOI: 10.1101/gr.129684.111. View

20.

Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N . The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009; 25(16):2078-9. PMC: 2723002. DOI: 10.1093/bioinformatics/btp352. View