Displaying somatic mutational signatures
1
1
Entering edit mode
8.2 years ago
kezcleal ▴ 160

Sorry if this is a very basic question. Ive been following the methodology of Alexandrov et al.,2013 (see http://cancer.sanger.ac.uk/cosmic/signatures) to find the mutational signature of a cancer sample and have become confused about how mutational signatures are supposed to be displayed on a graph, relating to the normalising of mutation counts. From the cosmic website:

Mutational signatures are displayed and reported based on the observed trinucleotide frequency of the human genome, i.e., representing the relative proportions of mutations generated by each signature based on the actual trinucleotide frequencies of the reference human genome version GRCh37.

Do I need to count the frequency of trinucleotide repeats in the human genome and then somehow normalise to this? I am a bit confused about what this sentence means.

Hypothetical example for a C>A substitution type and ACA context;

If the mutation frequency in my sample is 1% and the trinucleotide frequency of ACA in the human genome is 3.5%, what do I do next? Many thanks

snp genome • 3.7k views
ADD COMMENT
0
Entering edit mode

I am currently also confused by this statement. The way I currently interpret this is the following:

Each published signatures seems to sum up to 1 over the 96 dimensions. This, I assume, refers to the "relative proportion of the mutations generated by each signature". This leaves the normalization by the trinucleotide frequency: to ensure the normalization to 1 per signature still holds, I think this was done before the above.

I.e. given an un-normalized mutational signature S and the trinucleotide frequency T, the published signature would be: (S/T) / sum(S/T).

ADD REPLY
1
Entering edit mode
8.2 years ago
trausch ★ 1.9k

The Bioconductor package SomaticSignatures is very convenient for this:

https://www.bioconductor.org/packages/3.3/bioc/html/SomaticSignatures.html

ADD COMMENT
0
Entering edit mode

Thanks, I will give this a try

ADD REPLY
0
Entering edit mode

Does SomaticSignature package tells whether the identified signature is similar to signature-1 or signature-2 ........signature-30 from http://cancer.sanger.ac.uk/cosmic/signatures

I did not find in the manual where it said how similar it is to cosmic signatures.

ADD REPLY
0
Entering edit mode

In particular, https://www.bioconductor.org/packages/release/bioc/vignettes/SomaticSignatures/inst/doc/SomaticSignatures-vignette.html#nik-zainal_mutational_2012 should be helpful with regards to normalizing between e.g. exome and whole genome.

ADD REPLY

Login before adding your answer.

Traffic: 2671 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6