Skip to content

Annotation

Learning outcomes

After having completed this chapter you will be able to:

  • Describe the aims of variant annotation
  • Explain how variants are ranked in order of importance
  • Explain how splice variation affects variant annotation
  • Perform a variant annotation with snpEff
  • Interpret the report generated by snpEff
  • Explain how variant annotation can be added to a vcf file

Material

Download the presentation

Exercises

To use the human genome as a reference, we have downloaded the database with:

No need to download, it’s already downloaded for you

# don't run this. It's already downloaded for you
snpEff download -v GRCh38.99

You can run snpEff like so:

mkdir annotation

snpEff -Xmx4g \
-v \
-dataDir /data/ \
GRCh38.99 \
<filtered_variants.vcf> \
> <annotated_variants.vcf>

Exercise: Run the command on the filtered vcf (trio.filtered.vcf) using a script called C18_annotate_snpEff.sh. Check out the html file (snpEff_summary.html). Try to answer these questions:

A. How many effects were calculated?

B. How many variants are in the vcf?

C. Why is this different?

D. How many effects result in a missense mutation?

Answer

Your script:

C18_annotate_snpEff.sh
#!/usr/bin/env bash

cd ~/project/results/variants

snpEff -Xmx4g \
-v \
-dataDir /data/ \
GRCh38.99 \
trio.filtered.vcf \
> trio.filtered.snpeff.vcf

A. There were 10,357 effects calculated.

B. There are only 556 variants in the vcf.

C. This means that there are multiple effects per variant. snpEff calculates effects for each splice variant, and therefore the number of effects are a multitude of the number of variants.

D. Two effects result in a missense mutation.

You can (quick and dirty) query the annotated vcf for the missense mutation with grep.

Exercise: Find the variant causing the missense mutation (the line contains the string missense). And answer the following questions:

Hint
grep missense trio.filtered.snpeff.vcf 

Run the command and have a look at the SnpEff ANN field documentation. Answer the following questions:

A. How are the SNP annotations stored in the vcf?

B. What are the genotypes of the individuals?

C. Which amino acid change does it cause?

Answer

Find the line with the missense mutation like this:

grep missense annotation/trio.filtered.snpeff.vcf

This results in (truncated long line, scroll to the right to see more):

chr20   10049540        .       T       A       220.29  PASS    AC=1;AF=0.167;AN=6;BaseQRankSum=-6.040e-01;DP=85;ExcessHet=3.0103;FS=0.000;MLEAC=1;MLEAF=0.167;MQ=60.00;MQRankSum=0.00;QD=8.16;ReadPosRankSum=0.226;SOR=0.951;ANN=A|missense_variant|MODERATE|ANKEF1|ENSG00000132623|transcript|ENST00000378392.6|protein_coding|7/11|c.971T>A|p.Leu324Gln|1426/5429|971/2331|324/776||,A|missense_variant|MODERATE|ANKEF1|ENSG00000132623|transcript|ENST00000378380.4|protein_coding|6/10|c.971T>A|p.Leu324Gln|1300/5303|971/2331|324/776||       GT:AD:DP:GQ:PL  0/0:34,0:34:99:0,102,1163       0/1:17,10:27:99:229,0,492       0/0:24,0:24:72:0,72,811

A. SNP annotations are stored in the INFO field, starting with ANN=

B. The genotypes are homozygous reference for the father and son, and heterozygous for the mother. (find the order of the samples with grep ^#CHROM)

C. The triplet changes from cTg to cAg, resulting in a change from Leu (Leucine) to Gln (Glutamine).

Questions & Answers