Allele-specific copy number analysis of tumor samples with aneuploidy and tumor heterogeneity

We describe a bioinformatic tool, Tumor Aberration Prediction Suite (TAPS), for the identification of allele-specific copy numbers in tumor samples using data from Affymetrix SNP arrays. It includes detailed visualization of genomic segment characteristics and iterative pattern recognition for copy number identification, and does not require patient-matched normal samples. TAPS can be used to identify chromosomal aberrations with high sensitivity even when the proportion of tumor cells is as low as 30%. Analysis of cancer samples indicates that TAPS is well suited to investigate samples with aneuploidy and tumor heterogeneity, which is commonly found in many types of solid tumors.     

Allele-specific copy number analysis of tumor samples with aneuploidy and tumor heterogeneity

We describe a bioinformatic tool, Tumor Aberration Prediction Suite (TAPS), for the identification of allele-specific copy numbers in tumor samples using data from Affymetrix SNP arrays. It includes detailed visualization of genomic segment characteristics and iterative pattern recognition for copy number identification, and does not require patient-matched normal samples. TAPS can be used to identify chromosomal aberrations with high sensitivity even when the proportion of tumor cells is as low as 30%. Analysis of cancer samples indicates that TAPS is well suited to investigate samples with aneuploidy and tumor heterogeneity, which is commonly found in many types of solid tumors.     

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.     

Quantification of Normal Cell Fraction and Copy Number Neutral LOH in Clinical Lung Cancer Samples Using SNP Array Data

Background

Technologies based on DNA microarrays have the potential to provide detailed information on genomic aberrations in tumor cells. In practice a major obstacle for quantitative detection of aberrations is the heterogeneity of clinical tumor tissue. Since tumor tissue invariably contains genetically normal stromal cells, this may lead to a failure to detect aberrations in the tumor cells.

Principal Finding

Using SNP array data from 44 non-small cell lung cancer samples we have developed a bioinformatic algorithm that accurately models the fractions of normal and tumor cells in clinical tumor samples. The proportion of normal cells in combination with SNP array data can be used to detect and quantify copy number neutral loss-of-heterozygosity (CNNLOH) in the tumor cells both in crude tumor tissue and in samples enriched for tumor cells by laser capture microdissection.

Conclusion

Genome-wide quantitative analysis of CNNLOH using the CNNLOH Quantifier method can help to identify recurrent aberrations contributing to tumor development in clinical tumor samples. In addition, SNP-array based analysis of CNNLOH may become important for detection of aberrations that can be used for diagnostic and prognostic purposes.     

Identifying Human Genome-Wide CNV,LOH and UPD by Targeted Sequencing of Selected Regions

Copy-number variations (CNV), loss of heterozygosity (LOH), and uniparental disomy (UPD) are large genomic aberrations leading to many common inherited diseases, cancers, and other complex diseases. An integrated tool to identify these aberrations is essential in understanding diseases and in designing clinical interventions. Previous discovery methods based on whole-genome sequencing (WGS) require very high depth of coverage on the whole genome scale, and are cost-wise inefficient. Another approach, whole exome genome sequencing (WEGS), is limited to discovering variations within exons. Thus, we are lacking efficient methods to detect genomic aberrations on the whole genome scale using next-generation sequencing technology. Here we present a method to identify genome-wide CNV, LOH and UPD for the human genome via selectively sequencing a small portion of genome termed Selected Target Regions (SeTRs). In our experiments, the SeTRs are covered by 99.73%~99.95% with sufficient depth. Our developed bioinformatics pipeline calls genome-wide CNVs with high confidence, revealing 8 credible events of LOH and 3 UPD events larger than 5M from 15 individual samples. We demonstrate that genome-wide CNV, LOH and UPD can be detected using a cost-effective SeTRs sequencing approach, and that LOH and UPD can be identified using just a sample grouping technique, without using a matched sample or familial information.     

Normalization of Illumina Infinium whole-genome SNP data improves copy number estimates and allelic intensity ratios

Background  

Illumina Infinium whole genome genotyping (WGG) arrays are increasingly being applied in cancer genomics to study gene copy number alterations and allele-specific aberrations such as loss-of-heterozygosity (LOH). Methods developed for normalization of WGG arrays have mostly focused on diploid, normal samples. However, for cancer samples genomic aberrations may confound normalization and data interpretation. Therefore, we examined the effects of the conventionally used normalization method for Illumina Infinium arrays when applied to cancer samples.     

Assessing the significance of conserved genomic aberrations using high resolution genomic microarrays

Genomic aberrations recurrent in a particular cancer type can be important prognostic markers for tumor progression. Typically in early tumorigenesis, cells incur a breakdown of the DNA replication machinery that results in an accumulation of genomic aberrations in the form of duplications, deletions, translocations, and other genomic alterations. Microarray methods allow for finer mapping of these aberrations than has previously been possible; however, data processing and analysis methods have not taken full advantage of this higher resolution. Attention has primarily been given to analysis on the single sample level, where multiple adjacent probes are necessarily used as replicates for the local region containing their target sequences. However, regions of concordant aberration can be short enough to be detected by only one, or very few, array elements. We describe a method called Multiple Sample Analysis for assessing the significance of concordant genomic aberrations across multiple experiments that does not require a-priori definition of aberration calls for each sample. If there are multiple samples, representing a class, then by exploiting the replication across samples our method can detect concordant aberrations at much higher resolution than can be derived from current single sample approaches. Additionally, this method provides a meaningful approach to addressing population-based questions such as determining important regions for a cancer subtype of interest or determining regions of copy number variation in a population. Multiple Sample Analysis also provides single sample aberration calls in the locations of significant concordance, producing high resolution calls per sample, in concordant regions. The approach is demonstrated on a dataset representing a challenging but important resource: breast tumors that have been formalin-fixed, paraffin-embedded, archived, and subsequently UV-laser capture microdissected and hybridized to two-channel BAC arrays using an amplification protocol. We demonstrate the accurate detection on simulated data, and on real datasets involving known regions of aberration within subtypes of breast cancer at a resolution consistent with that of the array. Similarly, we apply our method to previously published datasets, including a 250K SNP array, and verify known results as well as detect novel regions of concordant aberration. The algorithm has been fully implemented and tested and is freely available as a Java application at http://www.cbil.upenn.edu/MSA.     

Genome-Wide Loss of Heterozygosity and DNA Copy Number Aberration in HPV-Negative Oral Squamous Cell Carcinoma and Their Associations with Disease-Specific Survival

Oral squamous cell cancer of the oral cavity and oropharynx (OSCC) is associated with high case-fatality. For reasons that are largely unknown, patients with the same clinical and pathologic staging have heterogeneous response to treatment and different probability of recurrence and survival, with patients with Human Papillomavirus (HPV)-positive oropharyngeal tumors having the most favorable survival. To gain insight into the complexity of OSCC and to identify potential chromosomal changes that may be associated with OSCC mortality, we used Affymtrix 6.0 SNP arrays to examine paired DNA from peripheral blood and tumor cell populations isolated by laser capture microdissection to assess genome-wide loss of heterozygosity (LOH) and DNA copy number aberration (CNA) and their associations with risk factors, tumor characteristics, and oral cancer-specific mortality among 75 patients with HPV-negative OSCC. We found a highly heterogeneous and complex genomic landscape of HPV-negative tumors, and identified regions in 4q, 8p, 9p and 11q that seem to play an important role in oral cancer biology and survival from this disease. If confirmed, these findings could assist in designing personalized treatment or in the creation of models to predict survival in patients with HPV-negative OSCC.     

Inferring loss-of-heterozygosity from unpaired tumors using high-density oligonucleotide SNP arrays

Loss of heterozygosity (LOH) of chromosomal regions bearing tumor suppressors is a key event in the evolution of epithelial and mesenchymal tumors. Identification of these regions usually relies on genotyping tumor and counterpart normal DNA and noting regions where heterozygous alleles in the normal DNA become homozygous in the tumor. However, paired normal samples for tumors and cell lines are often not available. With the advent of oligonucleotide arrays that simultaneously assay thousands of single-nucleotide polymorphism (SNP) markers, genotyping can now be done at high enough resolution to allow identification of LOH events by the absence of heterozygous loci, without comparison to normal controls. Here we describe a hidden Markov model-based method to identify LOH from unpaired tumor samples, taking into account SNP intermarker distances, SNP-specific heterozygosity rates, and the haplotype structure of the human genome. When we applied the method to data genotyped on 100 K arrays, we correctly identified 99% of SNP markers as either retention or loss. We also correctly identified 81% of the regions of LOH, including 98% of regions greater than 3 megabases. By integrating copy number analysis into the method, we were able to distinguish LOH from allelic imbalance. Application of this method to data from a set of prostate samples without paired normals identified known regions of prevalent LOH. We have developed a method for analyzing high-density oligonucleotide SNP array data to accurately identify of regions of LOH and retention in tumors without the need for paired normal samples.     

GPHMM: an integrated hidden Markov model for identification of copy number alteration and loss of heterozygosity in complex tumor samples using whole genome SNP arrays

There is an increasing interest in using single nucleotide polymorphism (SNP) genotyping arrays for profiling chromosomal rearrangements in tumors, as they allow simultaneous detection of copy number and loss of heterozygosity with high resolution. Critical issues such as signal baseline shift due to aneuploidy, normal cell contamination, and the presence of GC content bias have been reported to dramatically alter SNP array signals and complicate accurate identification of aberrations in cancer genomes. To address these issues, we propose a novel Global Parameter Hidden Markov Model (GPHMM) to unravel tangled genotyping data generated from tumor samples. In contrast to other HMM methods, a distinct feature of GPHMM is that the issues mentioned above are quantitatively modeled by global parameters and integrated within the statistical framework. We developed an efficient EM algorithm for parameter estimation. We evaluated performance on three data sets and show that GPHMM can correctly identify chromosomal aberrations in tumor samples containing as few as 10% cancer cells. Furthermore, we demonstrated that the estimation of global parameters in GPHMM provides information about the biological characteristics of tumor samples and the quality of genotyping signal from SNP array experiments, which is helpful for data quality control and outlier detection in cohort studies.     

Integrated study of copy number states and genotype calls using high-density SNP arrays

We propose a statistical framework, named genoCN, to simultaneously dissect copy number states and genotypes using high-density SNP (single nucleotide polymorphism) arrays. There are at least two types of genomic DNA copy number differences: copy number variations (CNVs) and copy number aberrations (CNAs). While CNVs are naturally occurring and inheritable, CNAs are acquired somatic alterations most often observed in tumor tissues only. CNVs tend to be short and more sparsely located in the genome compared with CNAs. GenoCN consists of two components, genoCNV and genoCNA, designed for CNV and CNA studies, respectively. In contrast to most existing methods, genoCN is more flexible in that the model parameters are estimated from the data instead of being decided a priori. GenoCNA also incorporates two important strategies for CNA studies. First, the effects of tissue contamination are explicitly modeled. Second, if SNP arrays are performed for both tumor and normal tissues of one individual, the genotype calls from normal tissue are used to study CNAs in tumor tissue. We evaluated genoCN by applications to 162 HapMap individuals and a brain tumor (glioblastoma) dataset and showed that our method can successfully identify both types of copy number differences and produce high-quality genotype calls.     

DNA Copy Number Aberrations in Breast Cancer by Array Comparative Genomic Hybridization

Array comparative genomic hybridization (CGH) has been popularly used for an-alyzing DNA copy number variations in diseases like cancer. In this study, we investigated 82 sporadic samples from 49 breast cancer patients using 1-Mb reso-lution bacterial artificial chromosome CGH arrays. A number of highly frequent genomic aberrations were discovered, which may act as "drivers" of tumor pro-gression. Meanwhile, the genomic profiles of four "normal" breast tissue samples taken at least 2 cm away from the primary tumor sites were also found to have some genomic aberrations that recurred with high frequency in the primary tu-mors, which may have important implications for clinical therapy. Additionally, we performed class comparison and class prediction for various clinicopathological pa-rameters, and a list of characteristic genomic aberrations associated with different clinicopathological phenotypes was compiled. Our study provides clues for further investigations of the underlying mechanisms of breast carcinogenesis.     

Clonal and non-clonal chromosome aberrations and genome variation and aberration.

The theoretical view that genome aberrations rather than gene mutations cause a majority of cancers has gained increasing support from recent experimental data. Genetic aberration at the chromosome level is a key aspect of genome aberration and the systematic definition of chromosomal aberrations with their impact on genome variation and cancer genome evolution is of great importance. However, traditionally, efforts have focused on recurrent clonal chromosome aberrations (CCAs). The significance of stochastic non-clonal chromosome aberrations (NCCAs) is discussed in this paper with emphasis on the simple types of NCCAs that have until recently been considered "non-significant background". Comparison of various subtypes of transitional and late-stage CCAs with simple and complex types of NCCAs has uncovered a dynamic relationship among NCCAs, CCAs, overall genomic instability, and karyotypic evolution, as well as the stochastic nature of cancer evolution. Here, we review concepts and methodologies to measure NCCAs and discuss the possible causative mechanism and consequences of NCCAs. This study raises challenging questions regarding the concept of cancer evolution driven by stochastic chromosomal aberration mediated genome irregularities that could have repercussions reaching far beyond cancer and organismal genomes.     

Micro-Scale Genomic DNA Copy Number Aberrations as Another Means of Mutagenesis in Breast Cancer

Introduction

In breast cancer, the basal-like subtype has high levels of genomic instability relative to other breast cancer subtypes with many basal-like-specific regions of aberration. There is evidence that this genomic instability extends to smaller scale genomic aberrations, as shown by a previously described micro-deletion event in the PTEN gene in the Basal-like SUM149 breast cancer cell line.

Methods

We sought to identify if small regions of genomic DNA copy number changes exist by using a high density, gene-centric Comparative Genomic Hybridizations (CGH) array on cell lines and primary tumors. A custom tiling array for CGH (244,000 probes, 200 bp tiling resolution) was created to identify small regions of genomic change, which was focused on previously identified basal-like-specific, and general cancer genes. Tumor genomic DNA from 94 patients and 2 breast cancer cell lines was labeled and hybridized to these arrays. Aberrations were called using SWITCHdna and the smallest 25% of SWITCHdna-defined genomic segments were called micro-aberrations (<64 contiguous probes, ∼ 15 kb).

Results

Our data showed that primary tumor breast cancer genomes frequently contained many small-scale copy number gains and losses, termed micro-aberrations, most of which are undetectable using typical-density genome-wide aCGH arrays. The basal-like subtype exhibited the highest incidence of these events. These micro-aberrations sometimes altered expression of the involved gene. We confirmed the presence of the PTEN micro-amplification in SUM149 and by mRNA-seq showed that this resulted in loss of expression of all exons downstream of this event. Micro-aberrations disproportionately affected the 5′ regions of the affected genes, including the promoter region, and high frequency of micro-aberrations was associated with poor survival.

Conclusion

Using a high-probe-density, gene-centric aCGH microarray, we present evidence of small-scale genomic aberrations that can contribute to gene inactivation. These events may contribute to tumor formation through mechanisms not detected using conventional DNA copy number analyses.     

Identifying patterns of DNA for tumor diagnosis using capillary electrophoresis-amplified fragment length polymorphism (CE-AFLP) screening

Amplified Fragment Length Polymorphism (AFLP) screening is a genome-wide genotyping strategy that has been widely used in plants and bacteria, but little has been reported concerning its use in humans. We investigated if the AFLP procedure could be coupled with high-throughput capillary electrophoresis (CE) for use in tumor diagnosis and classification. Using CE-AFLP, a series of molecular 'fingerprints' were generated for a set of gastric tumor and normal genomic DNA samples. The CE-AFLP procedure was qualitatively and quantitatively robust, and a variety of clustering tools were used to identify a specific DNA marker 'pattern' of 20 features that classified the tumor and normal samples to reasonable degrees of accuracy (Sensitivity 95%, Specificity 80%). The CE-AFLP-based approach also correctly classified 16 tumor samples, which in a previous study had exhibited no detectable genomic aberrations by comparative genome hybridization (CGH). This is the first reported application of CE-AFLP screening in tumor diagnosis. As the procedure is relatively inexpensive and requires minimal prior sequence knowledge and biological material, we suggest that CE-AFLP-based protocols may represent a promising new approach for DNA-based cancer screening and diagnosis.     

Molekulare Karyotypisierung in der genetischen Diagnostik

Molecular karyotyping by array comparative genomic hybridization (CGH) and single nucleotide polymorphism (SNP) arrays allows for a high resolution scan of the entire genome. It detects gains or losses (copy number variants) that might be the underlying cause of a genetic disorder. This technique is mainly applied to cases with syndromal or non-syndromal impaired (intellectual) development and is used to characterize genetic aberrations of tumor samples. Furthermore, molecular karyotyping might be useful for resolving prenatal cases with abnormal ultrasound findings. The purpose of this article is to explain the basic techniques, their limitations and strengths and to provide an outlook on future prospects.     

Visualization of genomic aberrations using Affymetrix SNP arrays

MOTIVATION: DNA copy number aberrations are frequently found in different types of cancer. Recent developments of microarray-based approaches have broadened the knowledge on number and structure of such aberrations. High-density single nucleotide polymorphism (SNP) microarrays provide an extremely high resolution with up to 500,000 SNPs per genome. Owing to the enormous amount of data the detection of common aberrations in large datasets is a great challenge. We describe a novel open source software tool--IdeogramBrowser--which was specifically designed for use with the Affymetrix SNP arrays. It provides an interactive karyotypic visualization of multiple aberration profiles and direct links to GeneCards. Visualization of consensus regions together with gene representation allows the explorative assessment of the data. AVAILABILITY: IdeogramBrowser and its source code are freely available under a creative commons license and can be obtained from http://www.informatik.uni-ulm.de/ni/staff/HKestler/ideo/. IdeogramBrowser is a platform independent Java application.     

Oncogenomic Approaches in Exploring Gain of Function of Mutant p53

Cancer is caused by the spatial and temporal accumulation of alterations in the genome of a given cell. This leads to the deregulation of key signalling pathways that play a pivotal role in the control of cell proliferation and cell fate. The p53 tumor suppressor gene is the most frequent target in genetic alterations in human cancers. The primary selective advantage of such mutations is the elimination of cellular wild type p53 activity. In addition, many evidences in vitro and in vivo have demonstrated that at least certain mutant forms of p53 may possess a gain of function, whereby they contribute positively to cancer progression. The fine mapping and deciphering of specific cancer phenotypes is taking advantage of molecular-profiling studies based on genome-wide approaches. Currently, high-throughput methods such as array-based comparative genomic hybridization (CGH array), single nucleotide polymorphism array (SNP array), expression arrays and ChIP-on-chip arrays are available to study mutant p53-associated alterations in human cancers. Here we will mainly focus on the integration of the results raised through oncogenomic platforms that aim to shed light on the molecular mechanisms underlying mutant p53 gain of function activities and to provide useful information on the molecular stratification of tumor patients.