Gene Mapping & Cancer Driver Annotation

Gene Mapping & SIDG's

All genes have an internal ID, allowing mapping to current and previous HGNC gene symbols, Ensembl Gene IDs (v91) and other external gene identifiers. All genes with an HGNC-approved symbol as of April 2018 are currently included, including those without a protein product.

Any dataset values that were mapped to genes without an official gene symbol have been discarded from processing, but continue to be available in raw data downloads.

Cancer Driver Annotation: Genes

The list of driver genes is the union of two complementary gene sets: intOGen & COSMIC Tier 1 .

intOGen Genes

The intOGen analysis pipeline uses seven different methods to identify cancer driver genes:

• Mutation count bias: dNdScv, CBaSE & MutPanning.
• Mutations in protein:
  • sequence: OncodriveCLUSTL
  • structure: HotMAPS
  • functional domains: smRegions
• Functional impact of mutations: OncodriveFML

COSMIC Tier 1 Genes

COSMIC Tier 1 gene requirements:

• Mutation patterns which strongly support involvement in cancer aetiology.
• Evidence of how gene’s dysfunction impacts hallmarks of cancer.
• Publications from two independent groups describing mutations in at least one type of cancer.

Method of Action

The method of action of a gene can be one of four values:

Given its roles assigned by intOGen and COSMIC, the method of action for each gene was chosen using to the following rules:

1. Ignore fusion unless it is the only method seen.
2. If there’s only one method of action across both intOGen and COSMIC, use that provided.
3. Otherwise choose ambiguous.

The original methods of action assigned by intOGen (intogen_moa) and COSMIC (cosmic_moa) are also available.

Cancer Driver Annotation: Mutations

The driver mutations table is the union of four datasets:

• CGI catalog of validated oncogenic mutations
• MSKCC Hotspots version 2
  • Includes a set of in-frame indels unique to this dataset
• BoostDM dataset
  • A new machine learning based classification of mutation in the intOGen cohorts
• CPV dataset
  • A catalog of germline cancer predisposition variants from Huang et.al., Cell, 2018

Processing of Mutation Datasets

The CGI, MSKCC and CPV datasets have chromosome coordinates recorded on GRCh37, and were remapped to GRCh38 using UCSC’s hg19 to hg38 chain file and liftOver utility. BoostDM is provided on GRCh38.

The files for each dataset were processed to extract the information required for a row in the table:

CGI catalog

Each row in the catalog_of_validated_oncogenic_mutations.tsv file contains a separate row for each nucleotide mutation. Mutations with the same effect (e.g. "missense", “ess_splice”) at the same location are merged together.

MSKCC Hotspots

Each entry in the mskcc_hotspots.json file (downloaded from the MSKCC API) groups mutations at an amino acid position in a gene. Missense and nonsense mutations in the same amino acid are split into separate rows in the table. Splice site mutations are grouped by amino acid position, i.e. by intron, which were split into separate donor and acceptor entries by grouping contiguous ranges of coordinates into new rows. The SNV-hotspots and INDEL-hotspots sheets of hotspots_v2.xls were parsed in parallel to add information required, such as genomic coordinates, which is not present in the mskcc_hotspots.json file. (Other information, such as EnsEMBL transcript IDs, is only found in mskcc_hotspots.json.)

BoostDM dataset

The BoostDM data consists of a separate .tsv file for each cancer type. Any mutation which has a “boostDM class” of “True” (which is a score > 0.5) in any tissue is loaded.

CPV dataset

Germ line predisposition variants taken from the S2A.Pathogenic_variants sheet in supplementary data file Germline predisposition variants from Huang et.al.

Post Processing and QC

Overlapping mutations with the same effect were merged together. The EnsEMBL transcript IDs were checked and updated where they were no longer current, and the genomic strand was added. Amino acid coordinates were mapped back to genomic coordinates to check that they match the recorded coordinates, and mutations had their footprint on the genomic sequence expanded to span whole codons. Mutations in codons which straddle an intron were split into two separate rows with separate genomic coordinates.

Driver Mutation Data Field Definitions