Question: Displaying somatic mutational signatures
1
gravatar for kezcleal
3.5 years ago by
kezcleal130
United Kingdom
kezcleal130 wrote:

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 • 2.0k views
ADD COMMENTlink modified 3.5 years ago by trausch1.3k • written 3.5 years ago by kezcleal130

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 REPLYlink written 16 months ago by Alex0
1
gravatar for trausch
3.5 years ago by
trausch1.3k
Germany
trausch1.3k wrote:

The Bioconductor package SomaticSignatures is very convenient for this:

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

 

ADD COMMENTlink written 3.5 years ago by trausch1.3k

Thanks, I will give this a try

ADD REPLYlink written 3.5 years ago by kezcleal130

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 REPLYlink written 2.4 years ago by Chirag Nepal2.2k

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 REPLYlink written 16 months ago by Alex0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 637 users visited in the last hour