Understanding DNA Copy Number Data

Transcription

1 Understanding DNA Copy Number Data Adam B. Olshen Department of Epidemiology and Biostatistics Helen Diller Family Comprehensive Cancer Center University of California, San Francisco May 20, 2010

2 Background The DNA sequence copy number at any locus in a genome is the number of copies of genomic DNA. The normal copy number is two for human autosomes. Copy number alterations are gains or losses of DNA They modify the function and/or expression of genes They are common in cancer: copy number is - Increased at sites of oncogenes; - Decreased at sites of tumor suppressor genes.

3 Array CGH: Hybridization and Analysis Forward hybridization Reverse hybridization Reference DNA Tumor DNA Reference DNA Tumor DNA Cy5 Cy3 Cy3 Cy5 Co-hybridization Co-hybridization Scan Scan Genomic markers (BAC, cdna, or oligo) Normal Loss Genomic markers (BAC, cdna, or oligo) Normal Gain Gain Loss Array analysis - resulting data are normalized log test over reference intensities for genomic markers

4 DNA Copy Number Arrays-Array CGH Platform # Colors Type of Markers # of Markers Commercial BAC 2 BACs ( kb) Yes 31,000 cdna 2 cdnas (1kb) ,000 No ROMA 2 long oligo 85,000 No 400,000 Agilent 2 long oligo 244,000 Yes 1,000,000 Nimblegen 2 long oligo 400,000 Yes Affymetrix 1 short oligo 500,000 Yes 1,800,000 Mips 1 short oligo 50,000 Yes Illumina beads 1,200,000 Yes Review: Pinkel & Albertson (2005, Nature Genetics)

5 Example Breast ROMA array with 9820 probes log 2 (Test Reference) Copies 7 Copies 6 Copies 5 Copies 4 Copies 3 Copies 2 Copies 1 Copy Genomic Position

6 Additional Complications 1. Tumor samples are a mixture of tumor and normal cells, and the amount of normal contamination is often unknown. 2. Tumor cells are not homogeneous. Gains and losses may occur in differing proportions of cells. 3. There is (usually) no gold standard of similar resolution.

7 Analysis Goals 1. Reconstruct the copy number state for the entire genome 2. Identify the regions of gain and loss 3. Divide the genome into regions of equal copy number

8 Analysis Approaches Smoothing (Hsu et al.; Tibshirani and Wang) Hidden Markov models (Fridlyand et al.; Guha et al.) Segmentation (Picard et al.; Venkatraman and Olshen) Methods compared by Lai et al. (2005, Bioinformatics)

9 Segmentation Approach log 2 (T R) Gains Losses Position in Chromosome

10 Circular Binary Segmentation View the data as if on a circle and segment into two arcs. Hence named circular binary segmentation (CBS). This results in two or three segments of the original data. Test statistic: T = max T ij, where 0 <i<j m Ȳ ij Z ij T ij = s ij (j i) 1 +(m j + i) 1, =(X i X j )/(j i), Ȳ ij Z ij =(X X i + X j X m )/(m j + i), and s 2 ij is the corresponding mean square error. Split data if P (T >T obs ) α (probability under null of no change-point) and recurse until no further splits. Estimate probability by permutation.

11 Permutations Real, T=30.9 Permuted 1, T=3.1 log 2 (T R) log 2 (T R) Position in Chromosome Position in Chromosome Permuted 2, T=3.1 Permuted 3, T=2.9 log 2 (T R) log 2 (T R) Position in Chromosome Position in Chromosome

12 Circular?

13 We Made it Faster 9820 probes, max probe count 824, time: 347s vs 13s. log2(t/r) Genomic Position Cyan permutation Black hybrid with early stopping

14 Gains and Losses via Plateau Plot

15 Gains and Losses vis Plateau Plot

16 Looking Across Samples

17 Clustering Samples

18 Defining Regions Patients _ Markers

19 Defining Regions cont. Patients _ Markers

20 Defining Regions cont. Identifying regions of high frequency gain or lossminimal common regions (MCRs) Genomic Identification of Significant Targets in Cancer (GISTIC) integration of magnitude as well as frequencies (Beroukhim, PNAS, 2007) Significance Testing for Aberrant Copy number (STAC) (Diskin et al., Genome Res., 2006) Multiple Sample Analysis (MSA) (Guttman et al., PLoS Genetics, 2007)

21 Advanced Topics 1. SNP Arrays, Genoytping and Allele-Specific Copy Number 2. Copy Number Variation 3. Clonality

22 SNP Arrays and Genotyping SNP arrays can be used for tradional copy number analysis. They can also used by used for genotyping: genome-wide association studies integration of copy number and loss of heterozygosity (LOH)

23 Two Regions of Normal Copy Number? Copy Number Copy Number Position Position

24 Two Regions of Normal Copy Number? Copy Number Copy Number Position Position TRUE TRUE Homozygotes Homozygotes FALSE FALSE Position Position

25 Allele-Specific Copy Number Traditional methods measure the sum of the copy numbers from the two parental chromosomes. This is total copy number. SNP arrays can be used for allele-specific copy number (PSCN). A copy number of 2 for a SNP could mean: 1+1(normal) or 0+2(both altered).

26 Raw Data Total Copy Number A Copy Number SQRT Copy Number SQRT Copy Number Genomic Position Genomic Position B Copy Number SQRT Copy Number Genomic Position

27 PSCBS Segmentation

28 An Example Sample data for PSCN of 1 (paternal chrom.) and 2 (maternal chrom.) Observe A B A B A B A B A B A B A B A B Maternal A A A B B B A A Paternal A B A B A A A A SNP SNP Genotype A Copy Number B Copy Number Total Copy Number 1 AA AB AA BB AB AB AA AA

29 PSCBS 1. Genotyping 2. Two Rounds of Segmentation 3. Calling Copy Number

30 Finding Heterozygotes Minimum of A Allele and B Allele Heterozygotes Homozygotes Minimize p l, p is points in the window, l in length of window

31 PSCBS 1. Genotyping 2. Two Rounds of Segmentation 3. Calling Copy Number

32 First Round of Segmentation If the total copy number has not changed, generally neither PSCN has changed. Start by running CBS on total copy number.

33 Second Round of Segmentation

34 Second Round of Segmentation The B-allele frequency (BAF) is B/(A+B)=B/Total CN. If PSCN changes, so should the BAF of the heterozygotes (based on paired normal). We use BAF from TumorBoost (Bengtsson et al., in press), which adjusts the tumor BAF based on the paired normal BAF. Segment on the decrease in heterozygosity ρ =2 BAF 0.5.

35 Second Round of Segmentation

36 PSCBS 1. Genotyping 2. Two Rounds of Segmentation 3. Calling Copy Number

37 Copy Number within Segments Four distinct clusters if unequal parent-specific copy number, three if equal copy number, two if LOH B Allele Copy Number Chromosome 8p A Allele Copy Number B Allele Copy Number Chromosome 12q A Allele Copy Number B Allele Copy Number Chromosome 1p A Allele Copy Number

38 Densities of BAFs

39 Parent-Specific CBS Algorithm 1. Assign genotypes based on minimum allele. 2. Run CBS on total copy number to get an estimate of the sum of the parent-specific copy numbers. 3. Run CBS within segments on 2 BAF 0.5 of heterozygotes. 4. Test for regions of LOH. 5. Test whether the parent-specific copy numbers are equal, whichisif the mode of the BAF = If not LOH and unequal, estimate the difference of the parent-specific copy numbers within every CBS segment using only the heterozygotes. 7. From 2) and 6) estimate the parent-specific copy numbers.

40 Reconstruction

41 Advanced Topics 1. SNP Arrays, Genotyping and Allele-Specific Copy Number 2. Copy Number Variation 3. Clonality

42 Copy Number Variation Copy number variations are gains or losses in the germ line >1Kb (Redon et al., 2006). They have been associated with familial cancer (Lucito et al., 2007) and other complex disease (Sebat et al., 2007). When analyzing cancer samples it is important to distinguish between variations and cancer aberrations. Large regions of gain or loss are aberrations; small regions could be either.

43 The Cancer Genome Atlas (TCGA) TCGA is an NIH-funded project whose goal is better understanding of cancer though large-scale sequencing. To decide which genes to sequence, Cancer Genome Characterization Centers (CGCCs) were set up. Each CGCC applies a different array technology to the same hundreds of glioblastoma, lung, and ovarian samples. TCGA samples are supposed to have matching normals.

44 ACGH Matched Sample log2(test Reference) Tumor 8 Copies 7 Copies 6 Copies 5 Copies 4 Copies 3 Copies 2 Copies 1 Copy log2(test Reference) Normal 8 Copies 7 Copies 6 Copies 5 Copies 4 Copies 3 Copies 2 Copies 1 Copy

45 Current Methods for Handling CNVs Compile data from Toronto database ( Then Screen out CNV regions from analysis Screen out regions within samples that may be due to CNVs Screen out regions defined by multiple samples that may be due to CNVs Ignore the problem completely

46 Predicting Copy Number Variation Segment both the tumor and matching normal sample. An observation is a small region of gain or loss in the tumor. The class label is whether the region is a CNV in the normal sample. The predictors are whether the region was found in the literature, the length of the region, etc. We used 43 matched pairs for training and 36 as a test set. Once the classification model is fit, it can be used to predict in cases where there is no matching sample.

47 Univariate Results Significant predictors Segment mean Gains and losses Other patients Literature (count) Literature (subjects) Length Segmental duplication region Not significant predictors Near centromere Near telomere

48 Multivariate Results in.lit< in.lit>=2.606 absmeansd< absmeansd>= in.lit< in.lit>= segdup< 0.5 segdup>=

49 Multivariate Results Best in literature models 80% accurate CART models 82% accurate Random Forests models 84% accurate Matched normals 94% accurate (no CNVs) Normals 70% accurate (all CNVs)

50 Advanced Topics 1. SNP Arrays, Genotyping and Allele-Specific Copy Number 2. Copy Number Variation 3. Clonality

51 Second Cancers When a second cancer appears, it could be a 1. Second primary - independent origin 2. Metastasis - clonal origin Distinguishing between the two has clinical importance: a metastasis is more serious Pathological decision criteria include histological type, stage, and anatomic location Molecular markers may improve decision making Begg with others have used LOH data for this purpose Can array DNA copy number data be used?

52 LOH and Clonality Clonal Case Ambiguous Case Independent Case Locus L R L R L R 1p q p q q p NA NA 11p NA NA NA NA NA NA 11q q q p q q q NA NA

53 DNA Copy Number Example Chromosome Arm (Genomic Order)

54 Lung Data 20 patients with paired non-small cell lung cancers (Girard et al., 2009) Samples hybridized to Agilent 244K array. Clinical and mutation data included. Clinical decision based on second primaries or metastases. The former may be spared adjuvant chemotherapy.

55 A Clonal Example (LR2=7.9E+23)

56 An Indep. Example (LR2=7.3E-06)

57 An Equivocal Example (LR2=0.3)

58 Lung Summary Copy number classification contradicted clinical classification in 4 of 20 cases. Additional 4 cases called equivocal by copy number. Three clinical indep. called clonal and supported by matching somatic mutations; one clinical clonal called indep. and no mutation data. Copy number data may be useful in clinical decision making.

59 Software CBS can be found in an R library called DNAcopy that is part of Bioconductor ( PSCBS can be found in the library PSCBS on R-Forge (r-forge-r-project.org). CNV and clonality methods can be found on Irina Ostrovnaya s web page.

60 References 1. Begg, C. et al. (2006). Statistical Tests for Clonality. Biometrics. 2. Beroukhim, R. et al. (2007). Assessing the significance of chromosomal aberrations in cancer: Methodology and application to glioma. PNAS Diksin, S. et al. (2004). STAC: A method for testing the significance DNA copy number aberrations across multiple array-cgh experiments. Genome Res Fridlyand, J. et al. (2004). Understanding Array CGH data. JMVA Girard, N. et al. (2009). Genomic and mutational profiling to assess clonal relationships between multiple non-small cell lung cancers. Clin Cancer Res Guha, S. et al. (2006). Bayesian Hidden Markov Modeling of Array CGH Data. Harvard University Biostatistics Working Paper Series. Working paper Guttman, M. et al. (2007). Assessing the significance of conserved genomic aberrations using high resolution genomic microarrays. PLoS Genetics 3 e Hardenbol, P. et al. (2003). Multiplexed genotyping with sequence-tagged molecular inversion probes. Nat. Biotechnol Hsu, L. et al. (2005). Denoising array-based comparative genomic hybridization data using wavelets. Bioinformatics Lai, WR. et al. (2005). Comparative analysis of algorithms for identifying amplifications and deletions in array CGH data. Bioinformatics Lucito, R. et al. (2007). Copy-Number variants in patients with a strong family history of pancreatic cancer. Cancer Biol Ther. Epub ahead of print.

61 References 12. Olshen, A. et al. (2010). Extension of circular binary segmentation to parent-specic copy number. Submitted. 13. Olshen, A., Venkatraman, E., Lucito, R. and Wigler, M. (2004). Circular Binary Segmentation for the analysis of array-based DNA copy number data. Biostatistics Ostrovnaya, A. et al. (2010). A metastasis or a second independent cancer? Evaluating the clonal origin of tumors using array copy number data. Stat Med. Epub ahead of print. 15. Ostrovnaya, A. et al. (2010). A classification model for distinguishing copy number variants from cancer-related alterations BMC Bioinformatics In press. 16. Picard, F. (2005). A statistical approach for array CGH data analysis. BMC Bioinformatics Pinkel, D. and Albertson, D.G. (2005). Array comparative genomic hybridization and its applications in cancer. Nat Genet. 37 S11-S Redon, R. et al. (2006). Global variation in copy number in the human genome. Nature Sebat, J. et al. (2007). Strong association of de novo copy number mutations with autism. Science Tibshirani, R. and Wang, P. (2008). Spatial smoothing and hotspot detection for CGH data using the fused lasso. Biostatistics Venkatraman, E.S. and Olshen, A.B. (2007). A faster circular binary segmentation algorithm for the analysis of array CGH data. Bioinformatics