Multi-Platform Whole-Genome Microarray Analyses Refine the Epigenetic Signature of Breast Cancer Metastasis with Gene Expression and Copy Number

Background

We have previously identified genome-wide DNA methylation changes in a cell line model of breast cancer metastasis. These complex epigenetic changes that we observed, along with concurrent karyotype analyses, have led us to hypothesize that complex genomic alterations in cancer cells (deletions, translocations and ploidy) are superimposed over promoter-specific methylation events that are responsible for gene-specific expression changes observed in breast cancer metastasis.

Methodology/Principal Findings

We undertook simultaneous high-resolution, whole-genome analyses of MDA-MB-468GFP and MDA-MB-468GFP-LN human breast cancer cell lines (an isogenic, paired lymphatic metastasis cell line model) using Affymetrix gene expression (U133), promoter (1.0R), and SNP/CNV (SNP 6.0) microarray platforms to correlate data from gene expression, epigenetic (DNA methylation), and combination copy number variant/single nucleotide polymorphism microarrays. Using Partek Software and Ingenuity Pathway Analysis we integrated datasets from these three platforms and detected multiple hypomethylation and hypermethylation events. Many of these epigenetic alterations correlated with gene expression changes. In addition, gene dosage events correlated with the karyotypic differences observed between the cell lines and were reflected in specific promoter methylation patterns. Gene subsets were identified that correlated hyper (and hypo) methylation with the loss (or gain) of gene expression and in parallel, with gene dosage losses and gains, respectively. Individual gene targets from these subsets were also validated for their methylation, expression and copy number status, and susceptible gene pathways were identified that may indicate how selective advantage drives the processes of tumourigenesis and metastasis.

Conclusions/Significance

Our approach allows more precisely profiling of functionally relevant epigenetic signatures that are associated with cancer progression and metastasis.     

Non-Invasive Prenatal Detection of Trisomy 13 Using a Single Nucleotide Polymorphism- and Informatics-Based Approach

Purpose

To determine how a single nucleotide polymorphism (SNP)- and informatics-based non-invasive prenatal aneuploidy test performs in detecting trisomy 13.

Methods

Seventeen trisomy 13 and 51 age-matched euploid samples, randomly selected from a larger cohort, were analyzed. Cell-free DNA was isolated from maternal plasma, amplified in a single multiplex polymerase chain reaction assay that interrogated 19,488 SNPs covering chromosomes 13, 18, 21, X, and Y, and sequenced. Analysis and copy number identification involved a Bayesian-based maximum likelihood statistical method that generated chromosome- and sample-specific calculated accuracies.

Results

Of the samples that passed a stringent DNA quality threshold (94.1%), the algorithm correctly identified 15/15 trisomy 13 and 49/49 euploid samples, for 320/320 correct copy number calls.

Conclusions

This informatics- and SNP-based method accurately detects trisomy 13-affected fetuses non-invasively and with high calculated accuracy.     

TAFFYS: An Integrated Tool for Comprehensive Analysis of Genomic Aberrations in Tumor Samples

Background

Tumor single nucleotide polymorphism (SNP) array is a common platform for investigating the cancer genomic aberration and the functionally important altered genes. Original SNP array signals are usually corrupted by noise, and need to be de-convoluted into absolute copy number profile by analytical methods. Unfortunately, in contrast with the popularity of tumor Affymetrix SNP array, the methods that are specifically designed for this platform are still limited. The complicated characteristics of noise in signals is one of the difficulties for dissecting tumor Affymetrix SNP array data, as they inevitably blur the distinction between aberrations and create an obstacle for the copy number aberration (CNA) identification.

Results

We propose a tool named TAFFYS for comprehensive analysis of tumor Affymetrix SNP array data. TAFFYS introduce a wavelet-based de-noising approach and copy number-specific signal variance model for suppressing and modelling the noise in signals. Then a hidden Markov model is employed for copy number inference. Finally, by using the absolute copy number profile, statistical significance of each aberration region is calculated in term of different aberration types, including amplification, deletion and loss of heterozygosity (LOH). The result shows that copy number specific-variance model and wavelet de-noising algorithm fits well with the Affymetrix SNP array signals, leading to more accurate estimation for diluted tumor sample (even with only 30% of cancer cells) than other existed methods. Results of examinations also demonstrate a good compatibility and extensibility for different Affymetrix SNP array platforms. Application on the 35 breast tumor samples shows that TAFFYS can automatically dissect the tumor samples and reveal statistically significant aberration regions where cancer-related genes locate.

Conclusions

TAFFYS provide an efficient and convenient tool for identifying the copy number alteration and allelic imbalance and assessing the recurrent aberrations for the tumor Affymetrix SNP array data.     

Whole-genome analysis of Fusarium graminearum insertional mutants identifies virulence associated genes and unmasks untagged chromosomal deletions

Background

Identifying pathogen virulence genes required to cause disease is crucial to understand the mechanisms underlying the pathogenic process. Plasmid insertion mutagenesis of fungal protoplasts is frequently used for this purpose in filamentous ascomycetes. Post transformation, the mutant population is screened for loss of virulence to a specific plant or animal host. Identifying the insertion event has previously met with varying degrees of success, from a cleanly disrupted gene with minimal deletion of nucleotides at the insertion point to multiple-copy insertion events and large deletions of chromosomal regions. Currently, extensive mutant collections exist in laboratories globally where it was hitherto impossible to identify all the affected genes.

Results

We used a whole-genome sequencing (WGS) approach using Illumina HiSeq 2000 technology to investigate DNA tag insertion points and chromosomal deletion events in mutagenised, reduced virulence F. graminearum isolates identified in disease tests on wheat (Triticum aestivum). We developed the FindInsertSeq workflow to localise the DNA tag insertions to the nucleotide level. The workflow was tested using four mutants showing evidence of single and multi-copy insertions in DNA blot analysis. FindInsertSeq was able to identify both single and multi-copy concatenation insertion sites. By comparing sequencing coverage, unexpected molecular recombination events such as large tagged and untagged chromosomal deletions, and DNA amplification were observed in three of the analysed mutants. A random data sampling approach revealed the minimum genome coverage required to survey the F. graminearum genome for alterations.

Conclusions

This study demonstrates that whole-genome re-sequencing to 22x fold genome coverage is an efficient tool to characterise single and multi-copy insertion mutants in the filamentous ascomycete Fusarium graminearum. In some cases insertion events are accompanied with large untagged chromosomal deletions while in other cases a straight-forward insertion event could be confirmed. The FindInsertSeq analysis workflow presented in this study enables researchers to efficiently characterise insertion and deletion mutants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1412-9) contains supplementary material, which is available to authorized users.     

Validation of multiple single nucleotide variation calls by additional exome analysis with a semiconductor sequencer to supplement data of whole-genome sequencing of a human population

Background

Validation of single nucleotide variations in whole-genome sequencing is critical for studying disease-related variations in large populations. A combination of different types of next-generation sequencers for analyzing individual genomes may be an efficient means of validating multiple single nucleotide variations calls simultaneously.

Results

Here, we analyzed 12 independent Japanese genomes using two next-generation sequencing platforms: the Illumina HiSeq 2500 platform for whole-genome sequencing (average depth 32.4×), and the Ion Proton semiconductor sequencer for whole exome sequencing (average depth 109×). Single nucleotide polymorphism (SNP) calls based on the Illumina Human Omni 2.5-8 SNP chip data were used as the reference. We compared the variant calls for the 12 samples, and found that the concordance between the two next-generation sequencing platforms varied between 83% and 97%.

Conclusions

Our results show the versatility and usefulness of the combination of exome sequencing with whole-genome sequencing in studies of human population genetics and demonstrate that combining data from multiple sequencing platforms is an efficient approach to validate and supplement SNP calls.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-673) contains supplementary material, which is available to authorized users.     

Distinct Genetic Alterations in Colorectal Cancer

Background

Colon cancer (CRC) development often includes chromosomal instability (CIN) leading to amplifications and deletions of large DNA segments. Epidemiological, clinical, and cytogenetic studies showed that there are considerable differences between CRC tumors from African Americans (AAs) and Caucasian patients. In this study, we determined genomic copy number aberrations in sporadic CRC tumors from AAs, in order to investigate possible explanations for the observed disparities.

Methodology/Principal Findings

We applied genome-wide array comparative genome hybridization (aCGH) using a 105k chip to identify copy number aberrations in samples from 15 AAs. In addition, we did a population comparative analysis with aCGH data in Caucasians as well as with a widely publicized list of colon cancer genes (CAN genes). There was an average of 20 aberrations per patient with more amplifications than deletions. Analysis of DNA copy number of frequently altered chromosomes revealed that deletions occurred primarily in chromosomes 4, 8 and 18. Chromosomal duplications occurred in more than 50% of cases on chromosomes 7, 8, 13, 20 and X. The CIN profile showed some differences when compared to Caucasian alterations.

Conclusions/Significance

Chromosome X amplification in male patients and chromosomes 4, 8 and 18 deletions were prominent aberrations in AAs. Some CAN genes were altered at high frequencies in AAs with EXOC4, EPHB6, GNAS, MLL3 and TBX22 as the most frequently deleted genes and HAPLN1, ADAM29, SMAD2 and SMAD4 as the most frequently amplified genes. The observed CIN may play a distinctive role in CRC in AAs.     

seqCNA: an R package for DNA copy number analysis in cancer using high-throughput sequencing

Background

Deviations in the amount of genomic content that arise during tumorigenesis, called copy number alterations, are structural rearrangements that can critically affect gene expression patterns. Additionally, copy number alteration profiles allow insight into cancer discrimination, progression and complexity. On data obtained from high-throughput sequencing, improving quality through GC bias correction and keeping false positives to a minimum help build reliable copy number alteration profiles.

Results

We introduce seqCNA, a parallelized R package for an integral copy number analysis of high-throughput sequencing cancer data. The package includes novel methodology on (i) filtering, reducing false positives, and (ii) GC content correction, improving copy number profile quality, especially under great read coverage and high correlation between GC content and copy number. Adequate analysis steps are automatically chosen based on availability of paired-end mapping, matched normal samples and genome annotation.

Conclusions

seqCNA, available through Bioconductor, provides accurate copy number predictions in tumoural data, thanks to the extensive filtering and better GC bias correction, while providing an integrated and parallelized workflow.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-178) contains supplementary material, which is available to authorized users.     

SNP rs10248565 in HDAC9 as a novel genomic aberration biomarker of lung adenocarcinoma in non-smoking women

Background

Numerous efforts have been made to elucidate the etiology and improve the treatment of lung cancer, but the overall five-year survival rate is still only 15%. Although cigarette smoking is the primary risk factor for lung cancer, only 7% of female lung cancer patients in Taiwan have a history of smoking. Since cancer results from progressive accumulation of genetic aberrations, genomic rearrangements may be early events in carcinogenesis.

Results

In order to identify biomarkers of early-stage adenocarcinoma, the genome-wide DNA aberrations of 60 pairs of lung adenocarcinoma and adjacent normal lung tissue in non-smoking women were examined using Affymetrix Genome-Wide Human SNP 6.0 arrays. Common copy number variation (CNV) regions were identified by ≥30% of patients with copy number beyond 2 ± 0.5 of copy numbers for each single nucleotide polymorphism (SNP) and at least 100 continuous SNP variant loci. SNPs associated with lung adenocarcinoma were identified by McNemar’s test. Loss of heterozygosity (LOH) SNPs were identified in ≥18% of patients with LOH in the locus. Aberration of SNP rs10248565 at HDAC9 in chromosome 7p21.1 was identified from concurrent analyses of CNVs, SNPs, and LOH.

Conclusion

The results elucidate the genetic etiology of lung adenocarcinoma by demonstrating that SNP rs10248565 may be a potential biomarker of cancer susceptibility.     

Reliable Single Cell Array CGH for Clinical Samples

Background

Disseminated cancer cells (DCCs) and circulating tumor cells (CTCs) are extremely rare, but comprise the precursors cells of distant metastases or therapy resistant cells. The detailed molecular analysis of these cells may help to identify key events of cancer cell dissemination, metastatic colony formation and systemic therapy escape.

Methodology/Principal Findings

Using the Ampli1™ whole genome amplification (WGA) technology and high-resolution oligonucleotide aCGH microarrays we optimized conditions for the analysis of structural copy number changes. The protocol presented here enables reliable detection of numerical genomic alterations as small as 0.1 Mb in a single cell. Analysis of single cells from well-characterized cell lines and single normal cells confirmed the stringent quantitative nature of the amplification and hybridization protocol. Importantly, fixation and staining procedures used to detect DCCs showed no significant impact on the outcome of the analysis, proving the clinical usability of our method. In a proof-of-principle study we tracked the chromosomal changes of single DCCs over a full course of high-dose chemotherapy treatment by isolating and analyzing DCCs of an individual breast cancer patient at four different time points.

Conclusions/Significance

The protocol enables detailed genome analysis of DCCs and thereby assessment of the clonal evolution during the natural course of the disease and under selection pressures. The results from an exemplary patient provide evidence that DCCs surviving selective therapeutic conditions may be recruited from a pool of genomically less advanced cells, which display a stable subset of specific genomic alterations.     

Systematic biases in DNA copy number originate from isolation procedures

Background

The ability to accurately detect DNA copy number variation in both a sensitive and quantitative manner is important in many research areas. However, genome-wide DNA copy number analyses are complicated by variations in detection signal.

Results

While GC content has been used to correct for this, here we show that coverage biases are tissue-specific and independent of the detection method as demonstrated by next-generation sequencing and array CGH. Moreover, we show that DNA isolation stringency affects the degree of equimolar coverage and that the observed biases coincide with chromatin characteristics like gene expression, genomic isochores, and replication timing.

Conclusion

These results indicate that chromatin organization is a main determinant for differential DNA retrieval. These findings are highly relevant for germline and somatic DNA copy number variation analyses.     

ExCyto PCR Amplification

Background

ExCyto PCR cells provide a novel and cost effective means to amplify DNA transformed into competent bacterial cells. ExCyto PCR uses host E. coli with a chromosomally integrated gene encoding a thermostable DNA polymerase to accomplish robust, hot-start PCR amplification of cloned sequences without addition of exogenous enzyme.

Results

Because the thermostable DNA polymerase is stably integrated into the bacterial chromosome, ExCyto cells can be transformed with a single plasmid or complex library, and then the expressed thermostable DNA polymerase can be used for PCR amplification. We demonstrate that ExCyto cells can be used to amplify DNA from different templates, plasmids with different copy numbers, and master mixes left on ice for up to two hours. Further, PCR amplification with ExCyto cells is comparable to amplification using commercial DNA polymerases. The ability to transform a bacterial strain and use the endogenously expressed protein for PCR has not previously been demonstrated.

Conclusions

ExCyto PCR reduces pipetting and greatly increases throughput for screening EST, genomic, BAC, cDNA, or SNP libraries. This technique is also more economical than traditional PCR and thus broadly useful to scientists who utilize analysis of cloned DNAs in their research.     

EGFR protein overexpression correlates with chromosome 7 polysomy and poor prognostic parameters in clear cell renal cell carcinoma

Background

The role of epidermal growth factor (EGF) and its receptor (EGFR) in the pathogenesis and progression of various malignant tumors has long been known, but there is still disagreement concerning prognostic significance of EGFR expression in clear cell renal cell carcinoma (CCRCC). The present study was designed to analyze more objectively the protein EGFR expression in CCRCC and to compare its value with EGFR gene copy number changes and clinicopathologic characteristics including patient survival.

Methods

The protein EGFR expression was analyzed immunohistochemically on 94 CCRCC, and gene copy number alterations of EGFR by FISH analysis on 41 CCRCC selected according to distinct membrane EGFR staining.

Results

Membrane EGFR expression in tumor cells was heterogeneous with respect to the proportion of positive cells and staining intensity. FISH analysis did not reveal EGFR gene amplification, while polysomy of chromosome 7 found in 41% was associated with higher EGFR membrane expression. Moreover, EGFR overexpression was associated with a higher nuclear grade, larger tumor size and shorter patient''s survival, while there was no connection with pathological stage.

Conclusion

In conclusion, the protein expression of EGFR had an impact on prognosis in patients with CCRCC, while an increased copy number of chromosome 7 could be the possible reason for EGFR protein overexpression in the absence of gene amplification.     

A meta-analysis of array-CGH studies implicates antiviral immunity pathways in the development of hepatocellular carcinoma

Background

The development and progression of hepatocellular carcinoma (HCC) is significantly correlated to the accumulation of genomic alterations. Array-based comparative genomic hybridization (array CGH) has been applied to a wide range of tumors including HCCs for the genome-wide high resolution screening of DNA copy number changes. However, the relevant chromosomal variations that play a central role in the development of HCC still are not fully elucidated.

Methods

In present study, in order to further characterize the copy number alterations (CNAs) important to HCC development, we conducted a meta-analysis of four published independent array-CGH datasets including total 159 samples.

Results

Eighty five significant gains (frequency ≥25%) were mostly mapped to five broad chromosomal regions including 1q, 6p, 8q, 17q and 20p, as well as two narrow regions 5p15.33 and 9q34.2-34.3. Eighty eight significant losses (frequency ≥25%) were most frequently present in 4q, 6q, 8p, 9p, 13q, 14q, 16q, and 17p. Significant correlations existed between chromosomal aberrations either located on the same chromosome or the different chromosomes. HCCs with different etiologies largely exhibited surprisingly similar profiles of chromosomal aberrations with only a few exceptions. Furthermore, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated that the genes affected by these chromosomal aberrations were significantly enriched in 31 canonical pathways with the highest enrichment observed for antiviral immunity pathways.

Conclusions

Taken together, our findings provide novel and important clues for the implications of antiviral immunity-related gene pathways in the pathogenesis and progression of HCC.     

Locked Nucleic Acid Pentamers as Universal PCR Primers for Genomic DNA Amplification

Background

Multiplexing technologies, which allow for simultaneous detection of multiple nucleic acid sequences in a single reaction, can save a lot of time, cost and labor compared to traditional single reaction detection methods. However, the multiplexing method currently used requires precise handiwork and many complicated steps, making a new, simpler technique desirable. Oligonucleotides containing locked nucleic acid residues are an attractive tool because they have strong affinities for their complementary targets, they have been used to avoid dimer formation and mismatch hybridization and to enhance efficient priming. In this study, we aimed to investigate the use of locked nucleic acid pentamers for genomic DNA amplification and multiplex genotyping.

Results

We designed locked nucleic acid pentamers as universal PCR primers for genomic DNA amplification. The locked nucleic acid pentamers were able to prime amplification of the selected sequences within the investigated genomes, and the resulting products were similar in length to those obtained by restriction digest. In Real Time PCR of genomic DNA from three bacterial species, locked nucleic acid pentamers showed high priming efficiencies. Data from bias tests demonstrated that locked nucleic acid pentamers have equal affinities for each of the six genes tested from the Klebsiella pneumoniae genome. Combined with suspension array genotyping, locked nucleic acid pentamer-based PCR amplification was able to identify a total of 15 strains, including 3 species of bacteria, by gene- and species-specific probes. Among the 32 species used in the assay, 28 species and 50 different genes were clearly identified using this method.

Conclusion

As a novel genomic DNA amplification, the use of locked nucleic acid pentamers as universal primer pairs in conjunction with suspension array genotyping, allows for the identification of multiple distinct genes or species with a single amplification procedure. This demonstrates that locked nucleic acid pentamer-based PCR can be utilized extensively in pathogen identification.     

Relationships linking amplification level to gene over-expression in gliomas

Background

Gene amplification is thought to promote over-expression of genes favouring tumour development. Because amplified regions are usually megabase-long, amplification often concerns numerous syntenic or non-syntenic genes, among which only a subset is over-expressed. The rationale for these differences remains poorly understood.

Methodology/Principal Finding

To address this question, we used quantitative RT-PCR to determine the expression level of a series of co-amplified genes in five xenografted and one fresh human gliomas. These gliomas were chosen because we have previously characterised in detail the genetic content of their amplicons. In all the cases, the amplified sequences lie on extra-chromosomal DNA molecules, as commonly observed in gliomas. We show here that genes transcribed in non-amplified gliomas are over-expressed when amplified, roughly in proportion to their copy number, while non-expressed genes remain inactive. When specific antibodies were available, we also compared protein expression in amplified and non-amplified tumours. We found that protein accumulation barely correlates with the level of mRNA expression in some of these tumours.

Conclusions/Significance

Here we show that the tissue-specific pattern of gene expression is maintained upon amplification in gliomas. Our study relies on a single type of tumour and a limited number of cases. However, it strongly suggests that, even when amplified, genes that are normally silent in a given cell type play no role in tumour progression. The loose relationships between mRNA level and protein accumulation and/or activity indicate that translational or post-translational events play a key role in fine-tuning the final outcome of amplification in gliomas.     

Whole-genome reconstruction and mutational signatures in gastric cancer

Background

Gastric cancer is the second highest cause of global cancer mortality. To explore the complete repertoire of somatic alterations in gastric cancer, we combined massively parallel short read and DNA paired-end tag sequencing to present the first whole-genome analysis of two gastric adenocarcinomas, one with chromosomal instability and the other with microsatellite instability.

Results

Integrative analysis and de novo assemblies revealed the architecture of a wild-type KRAS amplification, a common driver event in gastric cancer. We discovered three distinct mutational signatures in gastric cancer - against a genome-wide backdrop of oxidative and microsatellite instability-related mutational signatures, we identified the first exome-specific mutational signature. Further characterization of the impact of these signatures by combining sequencing data from 40 complete gastric cancer exomes and targeted screening of an additional 94 independent gastric tumors uncovered ACVR2A, RPL22 and LMAN1 as recurrently mutated genes in microsatellite instability-positive gastric cancer and PAPPA as a recurrently mutated gene in TP53 wild-type gastric cancer.

Conclusions

These results highlight how whole-genome cancer sequencing can uncover information relevant to tissue-specific carcinogenesis that would otherwise be missed from exome-sequencing data.     

Global copy number analyses by next generation sequencing provide insight into pig genome variation

Background

Copy number variations (CNVs) confer significant effects on genetic innovation and phenotypic variation. Previous CNV studies in swine seldom focused on in-depth characterization of global CNVs.

Results

Using whole-genome assembly comparison (WGAC) and whole-genome shotgun sequence detection (WSSD) approaches by next generation sequencing (NGS), we probed formation signatures of both segmental duplications (SDs) and individualized CNVs in an integrated fashion, building the finest resolution CNV and SD maps of pigs so far. We obtained copy number estimates of all protein-coding genes with copy number variation carried by individuals, and further confirmed two genes with high copy numbers in Meishan pigs through an enlarged population. We determined genome-wide CNV hotspots, which were significantly enriched in SD regions, suggesting evolution of CNV hotspots may be affected by ancestral SDs. Through systematically enrichment analyses based on simulations and bioinformatics analyses, we revealed CNV-related genes undergo a different selective constraint from those CNV-unrelated regions, and CNVs may be associated with or affect pig health and production performance under recent selection.

Conclusions

Our studies lay out one way for characterization of CNVs in the pig genome, provide insight into the pig genome variation and prompt CNV mechanisms studies when using pigs as biomedical models for human diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-593) contains supplementary material, which is available to authorized users.     

Using the realized relationship matrix to disentangle confounding factors for the estimation of genetic variance components of complex traits

Background

In the analysis of complex traits, genetic effects can be confounded with non-genetic effects, especially when using full-sib families. Dominance and epistatic effects are typically confounded with additive genetic and non-genetic effects. This confounding may cause the estimated genetic variance components to be inaccurate and biased.

Methods

In this study, we constructed genetic covariance structures from whole-genome marker data, and thus used realized relationship matrices to estimate variance components in a heterogenous population of ~ 2200 mice for which four complex traits were investigated. These mice were genotyped for more than 10,000 single nucleotide polymorphisms (SNP) and the variances due to family, cage and genetic effects were estimated by models based on pedigree information only, aggregate SNP information, and model selection for specific SNP effects.

Results and conclusions

We show that the use of genome-wide SNP information can disentangle confounding factors to estimate genetic variances by separating genetic and non-genetic effects. The estimated variance components using realized relationship were more accurate and less biased, compared to those based on pedigree information only. Models that allow the selection of individual SNP in addition to fitting a relationship matrix are more efficient for traits with a significant dominance variance.     

Deciphering clonality in aneuploid breast tumors using SNP array and sequencing data

Intra-tumor heterogeneity concerns the existence of genetically different subclones within the same tumor. Single sample quantification of heterogeneity relies on precise determination of chromosomal copy numbers throughout the genome, and an assessment of whether identified mutation variant allele fractions match clonal or subclonal copy numbers. We discuss these issues using data from SNP arrays, whole exome sequencing and pathologist purity estimates on several breast cancers characterized by ERBB2 amplification. We show that chromosomal copy numbers can only be estimated from SNP array signals or sequencing depths for subclonal tumor samples with simple subclonal architectures under certain assumptions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0470-7) contains supplementary material, which is available to authorized users.