I have used biomart quite a bit in the past to convert gene ID's or get dN/dS values, but I can't find any attribute in any obvious mart or dataset to retrieve background mutation rate.
This may just be a terminology problem, but I have dug around to no avail. Is there a way to access background mutation rate from biomart given a gene ID (say, TP53)?
Thanks very much for your help,
Marcus G.
•
written
by
marcusmgoldberg •
1 answer
I think background mutation rates are typically estimated from the data, not from population estimates. The reason for this is that different cancers will have different mutational profiles--there isn't a "normal human mutation rate" in that sense.
There are a number of tools for performing such analysis looking for evidence of mutational burden. MutSig and genome music are a couple of them. You might start by looking at these tools before rolling your own.
•
written
by
Sean Davis
Log in to answer this question.
-
retrieving dN and dS in BiomaRt
written by va.pedro •
-
BioMart : Retrieve set_gene_id (ENSG) and set_transcript_id (ENST) from SNPs rsIDs
written by amandinelecerfdefer •
-
Mysterious genes in my Biomart results (genes that were not part of original query)
written by adam.faranda
-
why does my biomaRt query return inconsistent dataset lists?
written by adam.faranda
-
TCGA calling pipeline query
written by adampennycuick
-
getBM error in R while retrieving IDs and missing GC values
written by manu.tiwari •
-
Ensembl biomart attribute for number of synonymous vs. non-synonymou sites
written by Max
-
dn/ds ratio for somatic mutation
written by deepaksrna •
-
Retrieving Dn Or Ds From Ensembl Using Biomart
written by Max
-
Obtaining Cds Reference Genome Sequences And Coordinance Using Biomart
written by Max
What do you mean by background mutation rate? In what context are you working (e.g., cancer)?
Sorry - I should I have more specific. I want to compare mutation rate in cancer cells vs normal human mutation rates (specifically synonymous, silent substitutions).
Background mutation rate over the entire genome would be sufficient, but it would obviously be better if it is possible to get gene specific rates, as well.