Researching SNPs (PKU)

Mutants walk the earth.

There are mutants on your block. There's a mutant right next door. There's one in your shoes.

We are, face it, all mutants. Humanity amounts to a string of mutations millennia-old, tweaks that gave us an edge over the next brute down the line. Upright posture; opposable thumb; prefrontal lobes.

The all-time specifically human thing is language. Familiar yet mysterious, language is the angelic part of us angel/beasts. It is spiritual, evanescent, fugitive - even gorgeous words like those have an unearthly shimmer. - from philly.com

---

Short Task Description

In this weeks task, we will research SNPs in the PAH gene that cause or do not cause a different phenotype of the phenylalanine hydroxylase. The chosen databases are: The public release of the Human Gene Mutation Database (HGMD), dbSNP, SNPdbe, OMIM and SNPedia. You may find a detailed task description here. A very short journal, which offers a small overview over the things we did can be found here.

Databases

Overview

<figtable id="tab:modelling_scores"> Key Values of the different SNP databases

Database	Last Update	Number of Entries	Number of Entries concerning PAH	Type of information	Sources	Curation/Verification	Comment
HGMD	public after 3 years (quarterly updated)	50,129 (only mis-/nonsense)	397 (only mis-/nonsense)	all types of mutations	current literature	manual and computerised search in current literature	too much advertising
dbSNP	Oct 2011	292,031,791	2590	SNPs, short in/dels, polymorphisms, others	submitted by registered sources (labs, institutes,.. )	clustering of identical submissions by NCBI
SNPdbe	Mar 2012	1,691,464	328	nonsyn. SNPs	Swissprot, dbSNP, PMD, OMIM, 1000 genomes	cf. sources	Predictions of functional effect, experimental evidence if available in source
OMIM	June 2012	21,257 (Summary entries)	1 (64 selected SNPs)	catalog of human genes and genetic disorders and traits	current literature	manually curated
SNPedia	continuous, Wiki-style	29,058	53	SNPs	publicly edited	publicly edited	get genotyped and predicted

</figtable>

Distribution of Mutations

<figtable id="tab:mutation_distri"> Distribution of SNPs from different databases on the PAH gene.

Distribution of SNPs from HGMD on the PAH gene.

Distribution of SNPs from SNPdbe on the PAH gene.

Distribution of SNPs from dbSNP on the PAH gene.

Distribution of missense SNPs from dbSNP on the PAH gene.

Distribution of silent SNPs from dbSNP on the PAH gene.

</figtable>

Mutation Hotspots

TODO

HGMD

http://www.hgmd.cf.ac.uk/ac/gene.php?gene=PAH

We extracted 397 disease-causing mis- and non-sense SNPs from HGMD. The reference sequence is [NM_000277.1.c].

dpSNP

http://www.ncbi.nlm.nih.gov/snp/?term=PAH

As dbSNP offer differentiated search options, we extracted 29 synonymous SNPs and 358 missense SNPs separately. The reference sequence is again [NM_000277.1.c].

OMIM

http://omim.org/entry/612349

We will use OMIM to select common variants and interesting SNPs for next week's task.

SNPdbe

http://www.rostlab.org/services/snpdbe/dosearch.php?id=name&val=PAH&organism2=human&organism1=

Following the task description, we looked at the conservation score, to determine quickly, whether a mutation found here is disease causing or not. Our reasoning can be found at this subpage.

SNPedia

http://www.snpedia.com/index.php/PAH

The SNPs concerning PAH in SNPedia all appear in dbSNP and there is no additional information attached so we do not investigate them in any separate way.

Mapping

We want to present as much information in our map of SNPs on the PAH gene as possible. As it contains information from different sources and multiple annotations for the same positions, some of the annotations may appear to conflict. Consult the descriptions above and the raw material for more information in these cases.

mapping at sequence

<figure id="fig:SNPMapping">

Mapping of SNP entries found in the different databases. The bar on the left symbolizes the sequence of the PAH-gene. Every horizontal line refers to at least one entry from any of the databases. The numbers on the left of the bar refer to their amino acid position. Colorcoding in the middle section which contains the different mutations:
red: disease causing according to HGMD
orange: missense, but no entry in HGMD
green: silent mutation
Colorcoding in the last column of the picture is for better overview and is according to the meaning of the effect from SNPdbe:
red: probably disease causing
green: probably not disease causing

</figure>

mapping at structure

In order to see where the SNPs which we already mapped to the sequence are in the structure we used the PDB structure of 2PAH and marked the residues which according to the colorcoding we used above. This means: red is a residue known to be disease causing, orange is a misssense-SNP with no entry in HGMD and green is a silent mutation. Just keep in mind, that several residues have different annotations according to their substitution and that the PDB-structure is not the structure of the complete gene. The structure which can be found in the PDB starts with residue number 118 in the original sequence.

<figure id="fig:SNPMappingstructuredisease">

Mapping of the residues which are known to be disease causing.

</figure><figure id="fig:SNPMappingstructuremutation">

Mapping of the residues which are not silent, but no entry in the HGMD were found.

</figure><figure id="fig:SNPMappingstructuresilent">

Mapping of the residues with silent variants.

</figure><figure id="fig:SNPMappingstructureall">

Mapping of all residues, from <xr id="fig:SNPMappingstructuredisease"/>, <xr id="fig:SNPMappingstructuremutation"/> and <xr id="fig:SNPMappingstructuresilent"/> with the same colorcoding. Note that for each residue we only used the color which is the strongest impact. Meaning diseasecausing>mutation>silent

</figure>

For the silent and the not in HGMD found SNPs, there is nothing surprising to be found in the structure, most of them are to be found on the outside of the protein, or in the coiled regions which can be seen in <xr id="fig:SNPMappingstructuremutation"/> and <xr id="fig:SNPMappingstructuresilent"/>. The other thing we expected was the low number of silent mutations, because of the problem of finding them.
For the diseasecausing residues on the other hand, we were quite surprised, that we found so many and for them to be located in the coiled regions almost as much as in the structural regions(<xr id="fig:SNPMappingstructuredisease"/>). Most of the SNPs were in the inside of the protein, where the cosubstrate-, iron- and substratebinding site is located, so in this part we found exactly what we expected.