This is a test version of Biostars. For the public version, visit https://www.biostars.org.
TMB Tumor Mutation Burden

HI,

Can someone help me with the formulae or a tool to identify the Tumor Mutation Burden from a Whole Exome Sequence?

Thanks,

Abilesh

tumor mutation burden cancer exome
• 14,390 views • link
updated by lincaijin1994 • lincaijin1994 • written by abibioinformatics • abibioinformatics

I am also interested in calculating the TMB. I have the list of "PASS" somatic mutations from GATK. Does anyone has used above software or other way to calculate the TMB?

I am looking for recommendations based on success to count the TMB. If you have calculated TMB, please describe the method that worked best for you.

• 0 views • link
written by sutturka sutturka

The TMB is defined as the total number of nonsynonymous mutations per coding area of a tumor genome. Initially, it was determined using whole exome sequencing, but due to the high costs and long turnaround time of this method, targeted panel sequencing is currently being explored to measure TMB.

samtools flagstat can be used for total base calculation.

• 0 views • link
written by Shicheng Guo Shicheng Guo

3 answers

You can start by this paper : https://genomemedicine.biomedcentral.com/articles/10.1186/s13073-017-0424-2

In the method part :

TMB was defined as the number of somatic, coding, base substitution, and indel mutations per megabase of genome examined.

• 0 views • link
written by Nicolas Rosewick Nicolas Rosewick

I think the "number of mutations per magabase" could be a misleading estimate of TMB. The deeper you sequence the more mutations per magabase you find since you detect more and more variants with low allele frequency. Maybe one should weight a variant by its frequency in order to compute TMB. (NB, I just glanced through the paper.)

• 0 views • link
written by dariober dariober

I think they sequence to ridiculously high median depth (500X) so they might have started to saturate by then?

That having been said it's a crude metric for 2018 to say the least. I'm surprised it even has a name.

• 1 views • link
written by Jeremy Leipzig Jeremy Leipzig

Note that the human exome size is ~30Mb. So you can take the number of somatic mutations in a given tumor sample an divide that by 30 to obtain the Mut/Mb value (normally > 4-6 is considered 'hyper mutation')

Or maybe you can try this tool

• 0 views • link
written by roy.granit roy.granit

Instead of using 30 as denominator, I would use tools like GATK CallableLoci to get the exact number of bases.

• 0 views • link
written by markus.riester markus.riester

I believe standard WXS capture sizes are around 50 mb. (eg; Agilent sure select)

• 0 views • link
written by poisonAlien poisonAlien

I would go with either the GATK CallableLoci or the capture size , that should be the one in use rather than something arbitrary.

• 0 views • link
written by ivivek_ngs ivivek_ngs

Hi. How to get the capture size of WXS. I have some bed files but I dont know which of the below was used to get the capture size.

[design ID]_Regions.bed - This BED file contains a single track of the target regions of interest that SureDesign used to select the probes. You can use this track to see the exact regions that the program was attempting to cover when selecting the probes.

[design ID]_Covered.bed - This BED file contains a single track of the genomic regions that are covered by one or more probes in the design. The fourth column of the file contains annotation information. You can use this file for assessing coverage metrics.
• 0 views • link
written by lincaijin1994 • lincaijin1994

Just to add a little more confusion to this topic, there is another method implemented in Varlociraptor:

Varlociraptor enables an uncertainty aware computation of the tumor mutational burden (TMB). TMB is usually defined as the number of somatic, non-synonymous coding mutations per megabase of the measured coding genome. ... the TMB is calculated as expected value over the posterior probabilities for each variant to be somatic. Hence, the TMB estimate properly considers the uncertainty in the data. Moreover, as we show a TMB estimate for each minimum allele frequency, it becomes possible to reason over the clonal structure of the tumor, instead of considering only a single overall number. We expect this to increase the predictive power of the TMB.

• 0 views • link
written by igor igor

Log in to answer this question.

More posts like this