<p>I'll follow up to the <a href="https://bioinformatics.stackexchange.com/a/638/292">great answer</a> from <a href="https://bioinformatics.stackexchange.com/users/57/kamil-s-jaron">Kamil S Jaron</a>:</p>

<p>Regarding predicting what the variant ("mutation" is a very loaded term) will do, there are a variety of tools. Chief among these are <a href="http://annovar.openbioinformatics.org/en/latest/" rel="noreferrer">annovar</a> and <a href="http://www.ensembl.org/info/docs/tools/vep/index.html" rel="noreferrer">VEP</a>. The general idea behind these is to classify the variants according to their overlap with genes, which codons they change (if any), how big that change is (e.g., changes in charge are more likely detrimental) and so on. One could also consider conservation, since if a position is highly conserved then changes in it are more likely to be detrimental.</p>

<p>If you really want to predict how a variant will change a protein's function then that usually requires prior knowledge about the proteins in question. Eventually someone will use machine learning to cull the literature and provide good predictions, but I haven't seen that yet.</p>

CC BY-SA 3.0
2017-06-11T16:47:51.127Z
Nice summary, to add, PhastCons is the standard for conservation scores, although not the best as the Sipel lab has many newer tools alternatives and GERP++ is also better IMO. Also if you're looking at just a few individual loci, nothing beats a stare at a genome browser like UCSC
Chris_Rands
CC BY-SA 3.0
2017-06-11T10:15:56.960Z
2017-08-14T09:51:31.817Z
2017-08-14T09:51:31.817Z
Devon Ryan
<p><strong>Part 1 : how to detect mutations</strong></p>

<p>The keywords you are searching for are "variant calling". Basically you have to map sequencing reads to a reference genome (or gene) and then estimate for each position of the genome if the observed difference of mapped reads and the reference is more likely a sequencing error or a mutation (in genomic glossary - variant). </p>

<p>Popular tools for variant calling are <a href="https://software.broadinstitute.org/gatk/" rel="noreferrer">GATK</a>, <a href="https://github.com/ekg/freebayes" rel="noreferrer">FreeBayes</a> or <a href="https://samtools.github.io/bcftools/howtos/variant-calling.html" rel="noreferrer">bcftools</a> (previously part as Samtools package).</p>

<p>The question you linked asks for quick alternatives to simply looking for variants. Indeed, you can just visualise the mapped reads to the reference sequence and see if the variant is there or not. CIGAR is just a notation of read alignment used in sam files (files with mapped reads), you can find good explanation of CIGAR strings <a href="http://genome.sph.umich.edu/wiki/SAM" rel="noreferrer">here</a>.</p>

<p>The follow-up about "How to estimate an effect of mutation" is in <a href="https://bioinformatics.stackexchange.com/a/644/292">@DevonRyan 's awesome answer</a>.</p>

CC BY-SA 3.0
2017-06-10T23:58:21.667Z
2017-06-12T11:36:07.287Z
2017-06-12T11:36:07.287Z
Kamil S Jaron
<h2>Part 3</h2>
<blockquote>
<p>I thought I'd add a <strong>part 3</strong> . Having found via VEP a truncation (missense or frameshift) or a missense, how do you find what are the consequences at the protein level?</p>
</blockquote>
<p>The effect a mutation has on a protein can be diverse and how much imbalance in function can the cell's equilibrium absorb depends on the component. As a result you get recessive, dominant and embryonically lethal mutations.
In this figure are shown some of the effects, divided into core and surface mutations:
<a href="https://i.sstatic.net/W8Faw.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/W8Faw.png" alt="effects" /></a></p>
<h2>Core variants</h2>
<p>Core mutations are the most black-and-white to analyse, but require a protein structure available. A severe truncation and a highly deleterious missense will result in a protein that will likely be degraded. The difference in folding Gibbs free energy (∆∆G, where <span class="math-container">$\Delta \Delta G = \Delta G_{\rm{mutant}} – \Delta G_{\rm{WT}}$</span>)) can be calculated for a missense, but how large this must be to result in a certain percentage being degraded depends on factors that cannot be calculated analytically. Alternatively, the Gibbs free energy of ligand binding is affected resulting in loss of catalysis.</p>
<p>To calculate the ∆∆G effect (stability), there are many tools:</p>
<ul>
<li>Online <a href="http://marid.bioc.cam.ac.uk/sdm2" rel="nofollow noreferrer">SDM</a> uses ML to estimate the ∆∆G quickly (do note that positive kcal/mol normally is a loss of stability, but here the sign is inverted so positive is a gain in stability</li>
<li>Rosetta toolkit. A highly flexible panel of movers  (samplers) can calculate ∆∆G</li>
<li>FoldX</li>
<li>MD runs</li>
<li><a href="http://www.sbg.bio.ic.ac.uk/%7Emissense3d/" rel="nofollow noreferrer">Missense3D</a> assesses a variant based on several structural criteria</li>
</ul>
<p>In the case of ligand binding modelling as above but using the ligand bound structure and calculating the binding energy as the difference between bound and unbound (far far away). For altered catalysis, modelling as above but using the transition state and assessing the difference in the transition hill can be done. However, an active site variant is very likely to affect either ligand binding, transition state binding (catalysis) or product binding that it's surely deleterious.</p>
<p>This is all (mostly) straightforward. If the domain is a regulatory domain however, there is loss of inhibition, which results in a gain of function (dominant variant).</p>
<h2>Surface variants</h2>
<p>Then stuff gets messier when you get to surface variants and these can cause dominant effects. However, if you are close to a phosphorylation, ubiquitination or other PTM site, as identified in <a href="https://www.phosphosite.org/homeAction.action" rel="nofollow noreferrer">phosphosite plus</a>, you likely have a case of deregulation. For linear motifs, <a href="http://elm.eu.org/" rel="nofollow noreferrer">ELM</a> is useful, but requires a lot of literature scanning to confirm any find. Alternative a loss of an interface may break the protein function, however, what binds where can only be ascertained by available structures, cross-linking mass spec or in some cases well defined Western experiments, so this data is rare.
For this type of variant there are two pipelines:</p>
<ul>
<li><a href="http://miscast.broadinstitute.org/" rel="nofollow noreferrer">Miscast</a> collates information on the protein</li>
<li><a href="https://venus.sgc.ox.ac.uk/" rel="nofollow noreferrer">VENUS (beta)</a> calculates ∆∆G and finds any relevant close-by PTM or gnomAD variants [disclaimer: this is mine].</li>
</ul>
<h2>Footnote</h2>
<p>As the variants can have a variety of effects, none of these <em>score</em> (unlike the genetics parts), but instead try to <em>inform</em> as to what may be going on.</p>

CC BY-SA 4.0
2021-01-05T15:49:51.597Z
2021-01-05T15:49:51.597Z
Matteo Ferla
<p>I read lecture notes about <a href="http://mcb.berkeley.edu/courses/mcb142/lecture%20topics/Dernburg/Lecture6_Chapter8_screenviewing.pdf" rel="nofollow noreferrer">mutations</a>, and wondered what kind of algorithms are there to detect mutations?</p>
<p>How do we know if the gene is mutated or whether it is a sequencing error?</p>
<p>I saw <a href="https://bioinformatics.stackexchange.com/questions/98/how-to-quickly-determine-mutations-in-a-read-of-a-sam-file">this</a> thread which is related, but I am not sure how to use the CIGAR information, e.g., is it 100% accurate? What is the underline mechanism to detect mutation? Is there a way to predict what the mutation will cause?</p>

<variant-calling><gene>
How to detect a mutation and predict its consequence?
# How to detect a mutation and predict its consequence?

(Asked by: username_0 on 2017-06-10)

I read lecture notes about [mutations](http://mcb.berkeley.edu/courses/mcb142/lecture%20topics/Dernburg/Lecture6_Chapter8_screenviewing.pdf), and wondered what kind of algorithms are there to detect mutations?

How do we know if the gene is mutated or whether it is a sequencing error?

I saw [this](https://bioinformatics.stackexchange.com/questions/98/how-to-quickly-determine-mutations-in-a-read-of-a-sam-file) thread which is related, but I am not sure how to use the CIGAR information, e.g., is it 100% accurate? What is the underline mechanism to detect mutation? Is there a way to predict what the mutation will cause?

**Tags:** variant-calling, gene

---

## 3 answers

### Answer by username_1 on 2017-06-11. Score: 10

I'll follow up to the [great answer](https://bioinformatics.stackexchange.com/a/638/292) from [Kamil S Jaron](https://bioinformatics.stackexchange.com/users/57/kamil-s-jaron):

Regarding predicting what the variant ("mutation" is a very loaded term) will do, there are a variety of tools. Chief among these are [annovar](http://annovar.openbioinformatics.org/en/latest/) and [VEP](http://www.ensembl.org/info/docs/tools/vep/index.html). The general idea behind these is to classify the variants according to their overlap with genes, which codons they change (if any), how big that change is (e.g., changes in charge are more likely detrimental) and so on. One could also consider conservation, since if a position is highly conserved then changes in it are more likely to be detrimental.

If you really want to predict how a variant will change a protein's function then that usually requires prior knowledge about the proteins in question. Eventually someone will use machine learning to cull the literature and provide good predictions, but I haven't seen that yet.

> **Comments:**
> **username_2**: Nice summary, to add, PhastCons is the standard for conservation scores, although not the best as the Sipel lab has many newer tools alternatives and GERP++ is also better IMO. Also if you're looking at just a few individual loci, nothing beats a stare at a genome browser like UCSC  

---

### Answer by username_3 on 2017-06-10. Score: 9

**Part 1 : how to detect mutations**

The keywords you are searching for are "variant calling". Basically you have to map sequencing reads to a reference genome (or gene) and then estimate for each position of the genome if the observed difference of mapped reads and the reference is more likely a sequencing error or a mutation (in genomic glossary - variant).

Popular tools for variant calling are [GATK](https://software.broadinstitute.org/gatk/), [FreeBayes](https://github.com/ekg/freebayes) or [bcftools](https://samtools.github.io/bcftools/howtos/variant-calling.html) (previously part as Samtools package).

The question you linked asks for quick alternatives to simply looking for variants. Indeed, you can just visualise the mapped reads to the reference sequence and see if the variant is there or not. CIGAR is just a notation of read alignment used in sam files (files with mapped reads), you can find good explanation of CIGAR strings [here](http://genome.sph.umich.edu/wiki/SAM).

The follow-up about "How to estimate an effect of mutation" is in [@DevonRyan 's awesome answer](https://bioinformatics.stackexchange.com/a/644/292).


---

### Answer by username_4 on 2021-01-05. Score: 3

## Part 3

> I thought I'd add a **part 3** . Having found via VEP a truncation (missense or frameshift) or a missense, how do you find what are the consequences at the protein level?

The effect a mutation has on a protein can be diverse and how much imbalance in function can the cell's equilibrium absorb depends on the component. As a result you get recessive, dominant and embryonically lethal mutations.
In this figure are shown some of the effects, divided into core and surface mutations:
[![effects](https://i.sstatic.net/W8Faw.png)](https://i.sstatic.net/W8Faw.png)

## Core variants

Core mutations are the most black-and-white to analyse, but require a protein structure available. A severe truncation and a highly deleterious missense will result in a protein that will likely be degraded. The difference in folding Gibbs free energy (∆∆G, where $\Delta \Delta G = \Delta G\_{\rm{mutant}} – \Delta G\_{\rm{WT}}$)) can be calculated for a missense, but how large this must be to result in a certain percentage being degraded depends on factors that cannot be calculated analytically. Alternatively, the Gibbs free energy of ligand binding is affected resulting in loss of catalysis.

To calculate the ∆∆G effect (stability), there are many tools:

* Online [SDM](http://marid.bioc.cam.ac.uk/sdm2) uses ML to estimate the ∆∆G quickly (do note that positive kcal/mol normally is a loss of stability, but here the sign is inverted so positive is a gain in stability
* Rosetta toolkit. A highly flexible panel of movers (samplers) can calculate ∆∆G
* FoldX
* MD runs
* [Missense3D](http://www.sbg.bio.ic.ac.uk/%7Emissense3d/) assesses a variant based on several structural criteria

In the case of ligand binding modelling as above but using the ligand bound structure and calculating the binding energy as the difference between bound and unbound (far far away). For altered catalysis, modelling as above but using the transition state and assessing the difference in the transition hill can be done. However, an active site variant is very likely to affect either ligand binding, transition state binding (catalysis) or product binding that it's surely deleterious.

This is all (mostly) straightforward. If the domain is a regulatory domain however, there is loss of inhibition, which results in a gain of function (dominant variant).

## Surface variants

Then stuff gets messier when you get to surface variants and these can cause dominant effects. However, if you are close to a phosphorylation, ubiquitination or other PTM site, as identified in [phosphosite plus](https://www.phosphosite.org/homeAction.action), you likely have a case of deregulation. For linear motifs, [ELM](http://elm.eu.org/) is useful, but requires a lot of literature scanning to confirm any find. Alternative a loss of an interface may break the protein function, however, what binds where can only be ascertained by available structures, cross-linking mass spec or in some cases well defined Western experiments, so this data is rare.
For this type of variant there are two pipelines:

* [Miscast](http://miscast.broadinstitute.org/) collates information on the protein
* [VENUS (beta)](https://venus.sgc.ox.ac.uk/) calculates ∆∆G and finds any relevant close-by PTM or gnomAD variants [disclaimer: this is mine].

## Footnote

As the variants can have a variety of effects, none of these *score* (unlike the genetics parts), but instead try to *inform* as to what may be going on.


---

