Cancer Genome Annotation
Exclusion of hypermutated samples:
Cancer genomes can be a very noisy source of data. It is estimated that an individual's tumour is caused by 5-10 driver mutations, but genome resequencing regularly reveals over 10,000 somatic mutations per tumour, with much larger numbers not unusual in hypermutated samples (we've seen samples with over 100,000 mutations each, the greatest being 178,763).
Across studies from different groups using different techniques, it is unclear whether these huge numbers reflect true hypermutation, substantial germline variation or technical artefacts. To try to improve the value of these data, we exclude any sample with over 15,000 mutations, as this immediately introduces huge noise into COSMIC.
We have also been developing and refining strategies to highlight the most significant mutations. Initially we flagged variants previosuly identified as SNPs, and added FATHMM scores to predict the relative significance of mutations. These methods have now been superseded by the The Cancer Mutation Census (CMC), which was released in v92. This identifies coding variants which drive different types of cancer. The CMC integrates all coding somatic mutations collected by COSMIC with biological and biochemical information from multiple sources, combining data obtained from manual curation and computational analyses. Metrics like ClinVar significance, dN/dS ratios, and variant frequencies in normal populations (gnomAD) have been integrated into this resource.
Ultimately, we aim to identify the most significant high-value data within cancer genomes, making it much easier to identify actionable biomarkers.
Copy Number Variants (CNV)
For Cancer Genome Project data (including the Cell Lines Project) copy number analysis was carried out using the Affymetrix SNP6.0 array in conjunction with a bespoke algorithm (PICNIC: Predicting Integral Copy Numbers In Cancer).
Where available, copy number data from TCGA and ICGC have been included in COSMIC (for samples already present in the database ie samples with mutations). All TCGA data included in COSMIC has been reanalysed using ASCAT 2.4. Please refer to Peter Van Loo et al. PNAS 107:16910-16915, 2010. for more information.
Definition of Minor Allele and Copy Number in tables:
- Minor Allele: the number of copies of the least frequent allele eg if ABB, minor allele = A ( 1 copy) and major allele = B ( 2 copies)
- Copy Number: the sum of the major and minor allele counts eg if ABB, copy number = 3
Definition of Gain and Loss:
We have introduced filtering thresholds to only display CNVs which are high level amplifications, homozygous deletions, or where there has been 'substantial loss' within an otherwise duplicated genome. We also use a higher threshold for amplification if genome duplication has occurred. We use average ploidy > 2.7 to define genome duplication.
- ICGC: Gain and Loss as defined in the original data.
- TCGA: (reanalysed with ASCAT 2.4,Peter Van Loo et al. PNAS 107:16910-16915, 2010. ) and
Cell Lines Project ( Affymetrix SNP6.0 array data analysed with PICNIC)
- average genome ploidy <= 2.7 AND total copy number >= 5
- OR average genome ploidy > 2.7 AND total copy number >= 9
- average genome ploidy <= 2.7 AND total copy number = 0
- OR average genome ploidy > 2.7 AND total copy number < ( average genome ploidy - 2.7 )
Where available, copy number data from TCGA and ICGC have been included in COSMIC (for samples already present in the database ie samples with mutations). All TCGA data included in COSMIC has been reanalysed using ASCAT 2.4. Please refer to Peter Van Loo et al. PNAS 107:16910-16915, 2010 for more information.
Gene expression level 3 data has been downloaded from the publicly accessible TCGA portal. The platform codes currently used to produce the COSMIC gene expression values are: IlluminaHiSeq_RNASeqV2, IlluminaGA_RNASeqV2, IlluminaHiSeq_RNASeq, and IlluminaGA_RNASeq.
Please note that as from COSMIC v71 we no longer show results from the array platforms AgilentG4502A_07_2 and AgilentG4502A_07_3. By using only RNAseq data we can show more results. This is because disagreement between the array and RNAseq data was quite common and resulted in the exclusion of data (see 'Qualitative merging of results' below).
For the RNASeq platforms we used the .trimmed.annotated.gene.quantification.txt, files which contain Level 3 expression data and used RPKM as a method of quantifying gene expression from RNA sequencing data by normalizing for total read length and the number of sequencing reads.[https://wiki.nci.nih.gov/display/TCGA/RNASeq]
For the RNASeqV2 platforms, the files used were rsem.genes.normalized_results, which contain Level 3 expression data produced using MapSplice to do the alignment and RSEM to perform the quantitation. [https://wiki.nci.nih.gov/display/TCGA/RNASeq+Version+2]
Qualitative merging of results
We downloaded methylation data for TCGA studies from the ICGC portal that were produced using the Infinium HumanMethylation450 beadchip. Only TCGA studies were downloaded as they include normal samples which are used to predict differential methylation. For the statistical test to be valid only studies with > 19 normal samples were analysed.
GRCh37/Hg19 genomic coordinates were derived from the probe description file from illumina. We have used hgLiftOver to map these loci on to the new GRCh28/Hg38 genome assembly.
Background (TCGA literature)
Then we corrected for multiple testing using the Bonferroni correction as follows:
the p-value of each locus (CpG) is multiplied by the total number of CpGs in the list.
If the corrected p-value is still below the error rate, the locus will be considered significant:
Corrected P-value= p-value * n (number of CpGs in the test) <0.05
In practice this means that a p-value < 0.0000001655 is significant.
Qualitative Representation of Results
Details of the anaysis performed can be found in Alexandrov L.B et al., Nature. 22;500(7463):415-21 (2013)
Reanalysis of TCGA data by the Cancer Genome Project (CGP), Sanger Institute
- Data was downloaded from CGhub as compressed summary TSV files.
The following settings were used:
- By Library Type = WXS
- By Platform = Illumina
- State = Live
- Remapping was performed using BWA-MEM.
- Pindel and CaVEMan were run on the tumour normal pairs.
- MAF, VCF and MAGE-TAB files were generated from the Pindel and CaVEMan output.
- Mutations from the VCF files were imported into COSMIC.