Review Article | Open Access
A. Mesut Erzurumluoglu, Santiago Rodriguez, Hashem A. Shihab, Denis Baird, Tom G. Richardson, Ian N. M. Day, Tom R. Gaunt, "Identifying Highly Penetrant Disease Causal Mutations Using Next Generation Sequencing: Guide to Whole Process", BioMed Research International, vol. 2015, Article ID 923491, 16 pages, 2015. https://doi.org/10.1155/2015/923491
Identifying Highly Penetrant Disease Causal Mutations Using Next Generation Sequencing: Guide to Whole Process
Recent technological advances have created challenges for geneticists and a need to adapt to a wide range of new bioinformatics tools and an expanding wealth of publicly available data (e.g., mutation databases, and software). This wide range of methods and a diversity of file formats used in sequence analysis is a significant issue, with a considerable amount of time spent before anyone can even attempt to analyse the genetic basis of human disorders. Another point to consider that is although many possess “just enough” knowledge to analyse their data, they do not make full use of the tools and databases that are available and also do not fully understand how their data was created. The primary aim of this review is to document some of the key approaches and provide an analysis schema to make the analysis process more efficient and reliable in the context of discovering highly penetrant causal mutations/genes. This review will also compare the methods used to identify highly penetrant variants when data is obtained from consanguineous individuals as opposed to nonconsanguineous; and when Mendelian disorders are analysed as opposed to common-complex disorders.
Next generation sequencing (NGS) and other high throughput technologies have brought new challenges concomitantly. The colossal amount of information that is produced has led researchers to look for ways of reducing the time and effort it takes to analyse the resulting data whilst also keeping up with the storage needs of the resulting files, which are in the magnitude of gigabytes each. The recently emerged variant call format (VCF) has somewhat provided a way out of this complex issue . Using a reference sequence and comparing it with the query sequence, only the differences (i.e., variants) between the two are encoded into a VCF file. Not only are VCF files substantially smaller in size (e.g., for whole-exome data, <300x in relation to BAM files which store all raw read alignments), they also make the data relatively easy to analyse since there are many bioinformatics tools (e.g., annotation and mutation effect prediction) which accept the VCF format as standard input. The genome analysis toolkit (GATK) made available by the Broad Institute also provides useful suggestions to bring a universal standard for the annotation and filtering of variants in VCF files . The abovementioned reasons have made the VCF the established format for the sharing of genetic variation produced from large sequencing projects (e.g., 1000 Genomes Project, NHLBI Exome Project, also known as EVS). However the VCF does have some disadvantages. The files can be information dense and initially difficult to understand and parse. Comprehensive information about the VCF and its companion software VFCtools  is available online (http://vcftools.sourceforge.net/).
Because of the substantial decrease in the price of DNA sequencing and SNP chip arrays , there has been a sharp increase in the number of genetic association studies being carried out, especially in the form of genome-wide association studies (GWAS, statistics available at http://www.genome.gov/gwastudies/). As whole genome sequencing (WGS) is prohibitively expensive for large genetic association studies [4–6], whole exome sequencing (WES) has emerged as the attractive alternative, where only the protein coding region of the genome (i.e., exome) is targeted and sequenced . The decision to carry out WES over WGS is not solely influenced by the cost which currently stands at one-third in comparison , but also by the fact that most of the known Mendelian disorders (~85%) are caused by mutations in the exome ; and reliably interpreting variation outside of the exome is still challenging as there is little consensus on interpreting their functional effects (even with ENCODE data  and noncoding variant effect prediction tools such as CADD , FATHMM-MKL , and GWAVA ). For complex diseases, WES can provide more evidence for causality compared to GWAS, assuming that the causal variants are exonic. This is because the latter uses linkage disequilibrium (LD) patterns between common markers  whereas WES directly associates the variant itself with the trait/disorder. Therefore using GWAS, especially in gene-dense regions, one cannot usually make conclusive judgements about which gene(s) is causal without further sequencing or functional analysis. WES has been successfully used in identifying and/or verifying over 300 causal variants for Mendelian disorders (statistics from http://omim.org/) (also see references [15, 16] for discussion of the use and benefits of WES in clinical genetics). WES currently stands at approximately $1000 for 50x read depth (variable prices, less for larger studies). However since there is a great deal of variation in the human genome , finding the causal variant(s), especially ones with low penetrance, is not going to be trivial. This problem can be exacerbated by the nature of the disorder(s) analysed. It is relatively easier to map variants causing rare monogenic diseases (when several affected individuals/families are available for analysis), as there is most likely to be a single variant present in the cases that is not in the controls; but in contrast, common complex (polygenic) disorders are much harder to dissect when searching for causal variants.
In this paper, our aims are to (i) provide a guide for genetic association studies dealing with sequencing data to identify highly penetrant variants (ii) compare the different approaches taken when data is obtained from unrelated or consanguineous individuals, and (iii) make suggestions about how to rank single nucleotide variation (SNV) and/or insertion/deletions (indels) following the standard filtering/ranking steps if there are several candidate variants, using annotated variants within VCF files as examples. To aid the process of analysing sequencing data obtained from consanguineous individuals, we have also made available an autozygosity mapping algorithm (AutoZplotter) which takes VCF files as input and enables manual identification of regions that have longer stretches of homozygosity than would be expected by chance.
2. Stage 1: Quality Control and Variant Calling
Before any genetic analysis, it is important to understand how the raw data were produced and processed to make better judgements about the reliability of the data received. Thorough quality control steps are required to ensure the reliability of the dataset. Lack of adequate prior quality control will inevitably lead to loss of statistical power and increase false positive and false negative findings. Fully comprehending each step during the creation of the dataset will have implications on the interpretation stage, where genotyping errors (also known as “phantom” mutations ) may turn out to be statistically associated (e.g., batch effects between case and control batches) or the causal variant may not be identified due to poorly applied quality control (QC) and/or filtering methods. The most fitting example for this comes from a recent primary ciliary dyskinesia (PCD) study , where the causal variant was only detected after the authors manually noticed an absence of reads in the relevant region of the genome (personal communication with authors). The subsequent variant was not only missing in the VCF files, but also in the initial BAM files, requiring remapping of reads. Another point of consideration from this finding would be that the authors knew where to look because the RSPH9 gene (the p.Lys268del mutation) was one of their a priori candidates . This is also an example demonstrating the importance of deep prior knowledge and screening for known variants as it is impossible for one to manually check the whole exome (or the genome) for sequencing and/or mapping errors.
2.1. Targeted Sequencing
As far as WES projects are concerned, questions about coverage arise right from the start (Figure 1). Since knowledge concerning exons in our own genome is far from complete, there are differing definitions about the human exome coordinates. Therefore, the targeted regions by the commonly used and commercially available Agilent SureSelect  and the Nimblegen SeqCap EZ  exome capture kits are not entirely overlapping . Thus it is possible that the missing regions of the exome due to the chosen probe kit may turn out to have the functional region in relation to the disorder analysed. One must also bear in mind that the kits available for targeting the exome are not fully efficient due to a certain quantity of poorly synthesized and/or designed probes not being able to hybridize to the target DNA. Next step is target enrichment where high coverage is vital as NGS machines produce more erroneous base calls compared to other techniques ; therefore, especially for rare variant analyses, it is important to have data with high average read depth (i.e., ≥50x).
2.2. Mapping Sequence Reads
The raw reads produced should then be aligned to a reference genome (e.g., GRCh38, see NCBI Genome Reference Consortium) and there are many open sources and widely applied tools (Table 1). However, solely depending on automated methods and software can leave many reads spanning insertions and deletions (indels) misaligned; therefore postreviewing the data for mismapping is always a good practice, especially in the candidate regions. Attempting to remap misaligned reads with a lower stringency using software such as Pindel would be an ideal way to go about solving such a problem . GATK also provides a base recalibration and indel realignment algorithm for this purpose.
|These are some of the many tools built for aligning reads produced from high throughput sequencing. Some have made speed their main purpose whereas others have paid more attention to annotating the files produced (such as mapping quality).|
Effective variant calling depends on accurate mapping to a dependable reference sequence. If available, using a population specific reference genome would be most ideal to filter out known neutral SNPs existing within the region of origin of the analysed subjects (e.g., East-Asian reference genome for subjects of Japanese origin). Inclusion of ambiguity codes (e.g., IUPAC codes) for known polyallelic variants to create a composite reference genome can also be useful (although not essential).
2.3. Variant Calling
There are many tools available for the identification of SNVs, indels, splice-site variants, and CNVs present in the query sequence(s). Each variant calling tool has advantages and disadvantages and has made compromises relating to issues such as speed of analysis, annotation, and reliability of the output file (Table 2). Separating true variation from sequencing artefacts still represents a considerable challenge. When dealing with very rare disorders, the candidate regions in the output VCF (or BAM) files should be reviewed either by reviewing the QC scores in the VCF or by visualising the alignments in IGV . Performing this step could highlight sequencing errors such as overcoverage (due to greater abundance of capture probes for the region or double capturing due to poorly discriminated probes hybridising to the same region) or undercoverage (due to probes not hybridising because of high variability in the region). For rare Mendelian disorders, since there is going to be a single causal variant it is more important to make sure that the variants in the dataset are reliable. Therefore setting strict parameters for read depth (e.g., ≥10x), base quality score (e.g., ≥100), and genotype quality scores (e.g., ≥100) initially can eliminate wrong base and genotype calls. This can then be adjusted subsequently (i.e., made less stringent) if no variants with a strong candidacy for causality are found after filtering (also see Best Practices section of GATK documentation for variant analysis).
|GATK, SOAPsnp, and SAMTools have constantly been cited in large genetic association projects indicating their ease of use, reliability, and functionality. However, this is also helped by the fact that they have additional features. There are other tools such as Beagle , IMPUTE2 , and MaCH  which have modules for SNP and genotype calling but are mostly used for their main purpose such as imputation and haplotype phasing.|
As mentioned above, there are many tools available for the identification of variants present in the query sequence (see Table 2). GATK  is one of the most established SNP discovery and genome analysis toolkits, with extensive documentation and helpful forums. It is a structured programming framework which makes use of the programming philosophy of MapReduce to solve the data management challenge of NGS by separating data access patterns from analysis algorithms. GATK is constantly updated and cited and also has a vibrant forum which is maintained continually.
SAMtools  is a variant caller which uses a Bayesian approach and has been used in many WGS and WES projects including the 1000 Genomes Project . SAMtools also offers many additional features such as alignment viewing and conversion to a BAM file. A recent study has compared GATK, SAMtools, and Atlas2 and found GATK to perform best in many settings (see reference  for details). However all three were highly consistent with an overlapping rate of ~90%. SOAPsnp is another highly used SNP and genotype caller and is part of the reliable SOAP family of bioinformatics tools (http://soap.genomics.org.cn/).
2.4. Additional Checks of Autozygosity
For data obtained from consanguineous families, confirming expected autozygosity (i.e., homozygous for alleles inherited from a common ancestor) would be an additional check worth carrying out. If the individual is the offspring of first cousins then the level of autozygosity would be approximately 6.25% () and 12.5% () for offspring of double first cousins (or uncle-niece unions, see Figure S1 in Supplementary Material available online at http://dx.doi.org/10.1155/2015/923491 for a depiction of these). These values will be higher in endogamous populations (e.g., for offspring of first cousins: 6.25% + autozygosity brought about due to endogamy, see Figure S3 for an example). Autozygosity could be checked by inspecting long runs of homozygosity (LRoH) for each individual by using tools such as Plink (for SNP chip data) , EXCLUDEAR (for SNP chip data) , AgilentVariantMapper (for WES data) , and AutoSNPa (for SNP chip data)  and dividing total autozygous regions by total length of autosomes in the human genome (can be obtained from http://www.ensembl.org/Homo_sapiens/Location/Genome). AutoZplotter (available to download in Supplementary Materials) that we developed takes VCF files as input, enabling easy and reliable visualisation and analysis of LRoH for any type of data (WGS, WES or SNP chip). The code (written in the Python programming language) can also be adapted relatively easily for use in analyses of other species.
3. Stage 2: Filtering/Ranking of Variants
Once the quality control process is complete and VCF files are deemed “analysis ready,” the approach taken will depend on the type of disorder analysed. For rare Mendelian disorders, many filtering and/or ranking steps can be taken to reduce the thousands of variants to a few strong candidates. Screening previously identified genes for causal variants is a good starting point. Carrying out this simple check will allow the identification of the causal variant even from a single proband thus saving time, effort, and funding. If no previously identified variant is found in the proband analysed, there are several steps which can be taken to identify novel mutations.
3.1. Using Prior Information to Rank/Filter Variants
Locus specific databases (see http://www.hgvs.org/dblist/glsdb.html for a comprehensive list) and “whole-genome” mutation databases such as HGMD , ClinVar , LOVD , and OMIM  are very informative resources for this task. Finding no previously identified variants indicates a novel variant in the proband analysed. For rare Mendelian disorders, the look for the variant can begin by removal of known neutral and/or common variants (≥0.1%) as this would provide a smaller subset of potentially causal variants. This is a pragmatic choice as Mendelian disease causal variants are likely to be very rare in the population or unique to the proband. If the latter is true, the variant should be absent from public databases. For this process to be thorough, an automated annotation tool such as Ensembl VEP or ANNOVAR can be used (see reference  for a review on the caveats of using these consequence predictors). Ensembl VEP enables incorporation of allele frequency (labelled as GMAF, global minor allele frequency) information from the EVS and the 1000 Genomes Project (see Supplementary Material and Methods for details).
3.2. Using Effect Prediction Algorithms to Rank/Filter Variants
Ranking this subset of variants based on consequence (e.g., stop gains would rank higher than missense) and scores derived from mutation prediction tools (e.g., “probably damaging” variants would rank higher than “possibly damaging” according to Polyphen-2 prediction) would enable assessment of the predicted impact of all rare mutations. It is important to understand what is assumed at each filtering/ranking stage; and comments are included about each assumption and their caveats in Figure 2.
For individuals of European ancestry, a VCF file will have between eighty and ninety thousand variants for WES (more for individuals with African ancestry ); and approximately a tenth will be variants with “predicted high impact” (also known as variants, i.e., rare nonsense, missense, splice-site acceptor or donor variants, exonic indels, and start losses ). There are many algorithms which predict the functional effect of these variants (Table 3). A large proportion of these algorithms utilize sequence conservation within a multiple sequence alignment (MSA) of homologous sequences to identify intolerant substitutions, for example, a substitution falling within a conserved region of the alignment is less likely to be tolerated than a substitution falling within a diverse region of the alignment (see reference  for a review). A handful of these algorithms also utilize structural properties, such as the protein secondary structure and solvent accessible surface area, in order to boost performance. Well known examples of a sequence-based and structure-based algorithm are SIFT  and PolyPhen , respectively. Newer software such as FATHMM  and MutPred , which use state-of-the-art hidden Markov models and machine learning paradigms, are worth using for their performance. There are also several tools such as CONDEL-2  which combine the output of several prediction tools to produce a consensus deleteriousness score. Although SIFT and Polyphen are highly cited tools, comparative analyses carried out by Thusberg et al. and Shihab et al. found FATHMM, MutPred, and SNPs&GO to perform better using the VariBench benchmarking dataset containing missense mutations [43, 46]. For predicting the effects of noncoding variants, FATHMM-MKL , GWAVA , and/or CADD  should be used. Also Human Splice Finder (latest: v3.0) can be used for intronic variants which predicts whether splicing is affected by the variant or not . Many of these tools can be incorporated into the analyses through the Ensembl website (http://www.ensembl.org/info/docs/tools/vep/index.html) where VCF files are annotated .
|Many methods have been developed to predict the functional effect of variants in the genome. Many of the tools listed above use different features and datasets to predict these effects. This is not an exhaustive list of all prediction tools but a collection of the most used/cited ones.|
*Comprehensive information about the prediction tool including accuracy, specificity, and sensitivity available in [43, 46]. N/A: not applicable. MCC: Matthew’s correlation coefficient. MCCs are obtained from .
These prediction algorithms are, as their name suggests, only there to make predictions about whether a variant is expected to be functionally disruptive or not. Thus their main purpose is to enable researchers to rank certain variants higher than others in order for them to be studied in a systematic way. Thus they do not “prove” anything about the causality of the variant. The variants predicted “deleterious” still require following up through replication and/or functional studies. Also disagreements amongst different tools can be observed which can lead to different interpretations about the evolutionary history of the variant (e.g., same function conserved throughout different species or a recently acquired function). Users of prediction algorithms should be aware of how these algorithms derive their predictions and then decide whether the tool can be generalized to their datasets. For example, those interested in somatic mutations should choose cancer-specific algorithms for example, FATHMM-Cancer  and SPF-Cancer , given that germline variant prediction algorithms are incapable of discriminating between cancer driver mutations and other germline mutations.
3.3. Further Filtering/Ranking
With current knowledge, there are approximately fifty synonymous mutations with proven causality, complex traits and Mendelian disorders combined . This is a very small proportion when compared to the thousands of published clinically relevant nonsynonymous (i.e., missense and nonsense) mutations. Therefore, when filtering variants for rare monogenic disorders, not taking noncoding variants and synonymous variants into account in the initial stages is a pragmatic choice. If ranking is preferred, then tools such as SilVA , which ranks all synonymous variants, and CADD  which ranks all variants (including synonymous variants) in the VCF files should be used.
Highly penetrant (Mendelian or common-complex) disease causal variants are expected to be very rare; therefore most of them should not appear in publicly available datasets. However filtering all variants present in dbSNP which is common practice should not be carried out as amplification and/or sequencing errors as well as potentially causal variants are known to make their way into this database (see references [53, 54] for details). Thus use of a MAF threshold (e.g., ≤0.1% in 1000 genomes and/or EVS) is a wiser choice in contrast to using absence in dbSNP as a filter. Upon completion of these steps, a smaller subset of variants with strong candidacy will remain for further follow up to determine causality.
Another initially pragmatic choice is to filter out all the annotations except for the “canonical” transcripts (i.e., longest transcript of a gene, if several exist) as this can reduce the amount of variants present in the Ensembl VEP (or ANNOVAR) annotated files considerably (~5x fold). However, this can be a problem for genes where the canonical transcript does not contain all the exons present within the gene, as a mutation which falls in an exon which is not present in the canonical transcript will not be observed in the filtered file (coded “CANONICAL” in Ensembl VEP annotated variants).
As many online tools are expected to keep logs of the processes undergoing in their servers, to protect the confidentiality of genetic information, downloading a local version of the chosen tools (or the VEP cache from the Ensembl website) is recommended. VEP also enables the incorporation of many other annotations (e.g., conservation scores, is variant position present in HGMD public version, whether variant is cited in PubMed), which will make the screening and filtering steps more manageable.
4. Stage 3: Building Evidence for Causality
Figure 3 suggests an example route to take to help differentiate causal variant(s) from noncausal ones for Mendelian disorders. At this stage one must gather all information that is available about the disorder and use it to determine which inheritance pattern fits the data and what complications there might be (e.g., the possibility of compound heterozygotes in disorders which show allelic heterogeneity). Figure S2 can be used to observe the contrast between the routes taken when analysing Mendelian (Figure 3) and complex disorders.
4.1. Public Data as a Source of Evidence
Having a candidate gene list based on previously published literature (e.g., by using OMIM or a disease/pathway specific database such as the Ciliome database ) and knowledge about the biology of the disorder (e.g., biological pathways) is useful. Software such as STRING and KEGG predicts protein-protein interactions using a variety of sources [56, 57]. SNPs3D has a user friendly interface which is designed to suggest candidates for different disorders . UCSC Gene Sorter (accessible from https://genome.ucsc.edu/) is another useful tool for collating a candidate gene list as it groups genes according to several features such as protein homology, co-expression and gene ontology (GO) similarity. Uniprot’s (http://www.uniprot.org/) Blast and Align functions can provide essential information about the crucial role a certain residue plays within a protein if it is highly conserved throughout many species. This is especially important for SNVs where the SNV loci itself should be causal (e.g., missense mutations, excludeing nonsense mutations as they truncate the gene product, thus the deleted segment of the protein requires further follow-up to prove causality, not just the loci where the mutation occurred as in other SNVs).
An example of the filtering process for an autosomal recessive disorder such as PCD is depicted in Figure 5. If several variants pass the filtering steps, information about the relevant genes should be gathered using databases such as GeneCards (http://www.genecards.org/) and NCBI Gene (http://www.ncbi.nlm.nih.gov/gene) for functional information, GEO Profiles (http://www.ncbi.nlm.nih.gov/geoprofiles) and Unigene (http://www.ncbi.nlm.nih.gov/unigene) for translational data about the gene’s product; and if available, one can check if a homologue is present in different species using databases such as HomoloGene (http://www.ncbi.nlm.nih.gov/homologene) and whether a similar phenotype is observed in model organisms. For example, if the disorder affects the cerebral cortex but the gene product is only active in the tissues located in the foot, then one cannot make a good argument about the identified variant in the respective gene as being “causal.”
There are many complications that may arise depending on the disorder such as genetic (locus) heterogeneity , allelic heterogeneity  and incomplete penetrance . Therefore gathering as many cases from the same family is helpful. However for very rare Mendelian disorders this may not be possible, thus it is important to seek other lines of evidence for causality (e.g., animal models, molecular analyses).
4.2. Mapping Causal Loci within Families
For rare Mendelian disorders, familial information can be crucial. The availability of an extended pedigree can be very informative in mapping which variant(s) fits the mode of inheritance in the case(s) and not in the unaffected members of the family (e.g., for autosomal recessive mutations, confirming heterozygosity in the parents is a must). This will provide linkage data where its importance is best displayed by Sobreira et al. where WES data from a single proband was sufficient in discovering the causal variants in two different families . Where available, previously published linkage data (i.e., associating a chromosomal region to a Mendelian disorder) should also made use of.
Traditionally a LOD score of 3 (Prob. = 1/1000) is required for a variant/region to be accepted as causal. Reaching this threshold requires many large families with many affected individuals. However this is not feasible for most disease causal variants (which are very rare by nature) and other lines of evidence such as animal knockouts, molecular studies and local sequence alignments (by using UniProt as mentioned above) are required to make a case for the causality of variants, especially mutations which are not stop gains (e.g., missense).
As mentioned previously, understanding the characteristics of a Mendelian disorder is important. If the disorder is categorised as “familial” (i.e., occurs more in families than by chance alone), which are usually very rare by nature, then availability of familial data becomes crucial, as unaffected members of the family are going to be the main source of information when determining neutral alleles. Any homozygous (and rare) stop gains, splice-site acceptor/donor variants and start losses in previously identified genes would be prime candidates.
Approach taken in families is different from the approaches taken when analysing common Mendelian disorders using unrelated individuals. For common Mendelian disorders (e.g., Finnish Heritage disorders [63–65]), fitting the dataset into a recessive inheritance model requires most (if not all) affected individuals to have two copies of the disease allele, enabling the identification of founder mutations as they will be overrepresented in the cases. These variants will be homozygous through endogamy and not consanguinity.
4.3. Autozygosity Mapping
For consanguineous subjects, the causal mutation usually lies within an autozygous region (characterised by long regions of homozygosity, LRoH, which are generally >5 Mb, see ), thus checking whether any candidate genes overlaps with an LRoH can narrow region(s) of interest. There are several tools which can identify LRoHs such as Plink, AutoSNPa and AgilentVariantMapper. We have made available a user-friendly python script (AutoZplotter) to plot heterozygosity/homozygosity status of variants in VCF files to allow for manual screening of short autozygous regions as well as LRoHs.
There are several software which can detect long runs of homozygosity reliably (>5 Mb), however they struggle to identify regions that are shorter. Therefore we developed AutoZplotter which plots homozygosity/heterozygosity state and enables quick visualisation of suspected autozygous regions (requires Xming or other X11 display server). These regions can then be followed up in more detail if any overlaps with a candidate gene/region. The input format of AutoZplotter is VCF thus it suits any type of genetic data (e.g., SNP array, WES, WGS). AutoZplotter was used for this purpose in a previous study by Alsaadi et al. .
4.5. Exceptional Cases
There can always be exceptional cases (in consanguineous families also) such as compound heterozygotes (i.e., individuals carrying different variants in the two copies of the same gene). This would require haplotype phasing and the confirmation of variant status (i.e., heterozygosity for one allele and absence of the other) in the parents and the proband(s) by sequencing of PCR amplicons containing variant or genotyping the variant directly. Beagle and HAPI-UR are two widely used haplotype phasing tools for their efficiency and speed [67, 68].
4.6. Identifying Highly Penetrant Variants for Common-Complex Disorders
For common complex disorders, identifying causal variants in outbred populations has proven to be a difficult and costly process (Supp. Figure S2); and these disorders can have many unknowns such as the significance of environmental factors on the disorder (see two examples of differential environmental influence on disease/traits in references [69, 70]) and epistasis . Many of the causal variants may be relatively rare (and almost always in heterozygous state) in the population introducing issues with statistical power. Traditional GWAS do not attempt to analyse them thus they are largely ignored, leaving a lot of heritability of common complex disorders unexplained. Analysing individuals with extreme phenotypes where the segregation of disease mimics autosomal recessive disorders (e.g., in consanguineous families) can be useful in identifying highly penetrant causal genes/mutations for complex disorders (e.g., obesity and leptin gene mutations ). The genetic influence in these individuals is predicted to be higher and is expected to have a single highly penetrant variant in homozygous state. These highly penetrant mutations can mimic Mendelian disorders causal variants. Therefore similar study designs can be used as stated above (e.g., Autozygosity/homozygosity mapping).
The NGS era has brought data management problems to traditional geneticists. Many data formats and bioinformatics tools have been developed to tackle this problem. One can easily be lost in the plethora of databases, data formats and tools. “Which tools are out there? How do I use it? What do I do next with the data I have?” are continually asked questions. This review aims to guide the reader in the rapidly changing and ever expanding world of bioinformatics. Figure 4 depicts a summary of the analysis process from DNA extraction to finding the causal variant, putting into perspective which file formats are expected at each step and which bioinformatics tools we prefer due to reasons mentioned before. Researchers can then appreciate the stage that they are at and how many other steps are required for completion as well as knowing what to do at each step.
Whole exome sequencing is the current gold standard in the discovery of highly penetrant disease causal mutations. As knowledge on the noncoding parts of the genome can still be considered to be in its early days, the human exome is still a pragmatic target for many. As approx. 1600 known Mendelian disorders (and ~3500 when suspected ones are included) and most common-complex disorders are still waiting for their molecular basis to be figured out (from http://omim.org/statistics/entry, true as of 15/07/14), future genetic studies have much to discover. However for these projects to be fruitful, careful planning is needed to make full use of available tools and databases (see Table 4).
|The most important factors when carrying out a genetic association study are (i) the availability of reliable data (ii) bioinformatics and biological expertise, and (iii) careful planning.|
Finally, with this paper we have also made AutoZplotter available (input format: VCF), which plots homozygosity/heterozygosity state and enables quick visualisation of suspected autozygous regions. This can be important for shorter autozygous regions where other autozygosity mappers struggle.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
A. Mesut Erzurumluoglu, Denis Baird, and Tom G. Richardson are PhD students funded by the Medical Research Council (MRC). This work was also supported by the MRC (MC_UU_12013/8 and G1000427).
S1: List of useful bioinformatics command, parameter and files for variant calling, annotation and analysis.
Figure S1: Pedigree trees demonstrating examples of relevant consanguineous unions.
Figure S2: Example of the routes taken when analysing complex disorders.
- P. Danecek, A. Auton, G. Abecasis et al., “The variant call format and VCFtools,” Bioinformatics, vol. 27, no. 15, pp. 2156–2158, 2011.
- A. McKenna, M. Hanna, E. Banks et al., “The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data,” Genome Research, vol. 20, no. 9, pp. 1297–1303, 2010.
- M. L. Metzker, “Sequencing technologies—the next generation,” Nature Reviews Genetics, vol. 11, no. 1, pp. 31–46, 2010.
- L. Bonetta, “Whole-genome sequencing breaks the cost barrier,” Cell, vol. 141, no. 6, pp. 917–919, 2010.
- E. Pettersson, J. Lundeberg, and A. Ahmadian, “Generations of sequencing technologies,” Genomics, vol. 93, no. 2, pp. 105–111, 2009.
- D. J. Hedges, “Comparison of three targeted enrichment strategies on the SOLiD sequencing platform,” PLoS ONE, vol. 6, Article ID e18595, 2011.
- J. K. Teer and J. C. Mullikin, “Exome sequencing: the sweet spot before whole genomes,” Human Molecular Genetics, vol. 19, no. 2, pp. R145–R151, 2010.
- D. Bick and D. Dimmock, “Whole exome and whole genome sequencing,” Current Opinion in Pediatrics, vol. 23, no. 6, pp. 594–600, 2011.
- M. Choi, U. I. Scholl, W. Ji et al., “Genetic diagnosis by whole exome capture and massively parallel DNA sequencing,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 45, pp. 19096–19101, 2009.
- The ENCODE Project Consortium, “The ENCODE (ENCyclopedia Of DNA Elements) project,” Science, vol. 306, no. 5696, pp. 636–640, 2004.
- M. Kircher, D. M. Witten, P. Jain, B. J. O'Roak, G. M. Cooper, and J. Shendure, “A general framework for estimating the relative pathogenicity of human genetic variants,” Nature Genetics, vol. 46, no. 3, pp. 310–315, 2014.
- H. A. Shihab, M. F. Rogers, J. Gough et al., “An integrative approach to predicting the functional effects of non-coding and coding sequence variation,” Bioinformatics, 2015.
- G. R. Ritchie, I. Dunham, E. Zeggini, and P. Flicek, “Functional annotation of noncoding sequence variants,” Nature Methods, vol. 11, no. 3, pp. 294–296, 2014.
- A. Kiezun, K. Garimella, R. Do et al., “Exome sequencing and the genetic basis of complex traits,” Nature Genetics, vol. 44, no. 6, pp. 623–630, 2012.
- C.-S. Ku, N. Naidoo, and Y. Pawitan, “Revisiting Mendelian disorders through exome sequencing,” Human Genetics, vol. 129, no. 4, pp. 351–370, 2011.
- C. Gilissen, A. Hoischen, H. G. Brunner, and J. A. Veltman, “Disease gene identification strategies for exome sequencing,” European Journal of Human Genetics, vol. 20, no. 5, pp. 490–497, 2012.
- The 1000 Genomes Project Consortium, “A map of human genome variation from population-scale sequencing,” Nature, vol. 467, no. 7319, pp. 1061–1073, 2010.
- A. Brandstätter, T. Sänger, S. Lutz-Bonengel et al., “Phantom mutation hotspots in human mitochondrial DNA,” Electrophoresis, vol. 26, no. 18, pp. 3414–3429, 2005.
- M. M. Alsaadi, T. R. Gaunt, C. R. Boustred et al., “From a single whole exome read to notions of clinical screening: primary ciliary dyskinesia and RSPH9 p.Lys268del in the Arabian Peninsula,” Annals of Human Genetics, vol. 76, no. 3, pp. 211–220, 2012.
- V. H. Castleman, L. Romio, R. Chodhari et al., “Mutations in radial spoke head protein genes RSPH9 and RSPH4A cause primary ciliary dyskinesia with central-microtubular-pair abnormalities,” American Journal of Human Genetics, vol. 84, no. 2, pp. 197–209, 2008.
- A. Gnirke, A. Melnikov, J. Maguire et al., “Solution hybrid selection with ultra-long oligonucleotides for massively parallel targeted sequencing,” Nature Biotechnology, vol. 27, no. 2, pp. 182–189, 2009.
- M. N. Bainbridge, M. Wang, D. L. Burgess et al., “Whole exome capture in solution with 3 Gbp of data,” Genome Biology, vol. 11, article R62, 2010.
- A.-M. Sulonen, P. Ellonen, H. Almusa et al., “Comparison of solution-based exome capture methods for next generation sequencing,” Genome biology, vol. 12, no. 9, article R94, 2011.
- E. Y. Chan, “Next-generation sequencing methods: impact of sequencing accuracy on SNP discovery,” in Single Nucleotide Polymorphisms, A. A. Komar, Ed., pp. 95–111, Humana Press, 2009.
- K. Ye, M. H. Schulz, Q. Long, R. Apweiler, and Z. Ning, “Pindel: a pattern growth approach to detect break points of large deletions and medium sized insertions from paired-end short reads,” Bioinformatics, vol. 25, no. 21, pp. 2865–2871, 2009.
- H. Thorvaldsdóttir, J. T. Robinson, and J. P. Mesirov, “Integrative genomics viewer (IGV): high-performance genomics data visualization and exploration,” Briefings in Bioinformatics, vol. 14, no. 2, pp. 178–192, 2013.
- H. Li, B. Handsaker, A. Wysoker et al., “The Sequence Alignment/Map format and SAMtools,” Bioinformatics, vol. 25, no. 16, pp. 2078–2079, 2009.
- X. Liu, S. Han, Z. Wang, J. Gelernter, and B.-Z. Yang, “Variant callers for next-generation sequencing data: a comparison study,” PLoS ONE, vol. 8, no. 9, Article ID e75619, 2013.
- S. Purcell, B. Neale, K. Todd-Brown et al., “PLINK: a tool set for whole-genome association and population-based linkage analyses,” The American Journal of Human Genetics, vol. 81, no. 3, pp. 559–575, 2007.
- C. G. Woods, E. M. Valente, J. Bond, and E. Roberts, “A new method for autozygosity mapping using single nucleotide polymorphisms (SNPs) and EXCLUDEAR,” Journal of Medical Genetics, vol. 41, no. 8, article e101, 2004.
- I. M. Carr, S. Bhaskar, J. O' Sullivan et al., “Autozygosity Mapping with Exome Sequence Data,” Human Mutation, vol. 34, no. 1, pp. 50–56, 2013.
- I. M. Carr, K. J. Flintoff, G. R. Taylor, A. F. Markham, and D. T. Bonthron, “Interactive visual analysis of SNP data for rapid autozygosity mapping in consanguineous families,” Human Mutation, vol. 27, no. 10, pp. 1041–1046, 2006.
- P. D. Stenson, E. V. Ball, M. Mort, A. D. Phillips, K. Shaw, and D. N. Cooper, “UNIT 1.13 The Human Gene Mutation Database (HGMD) and its exploitation in the fields of personalized genomics and molecular evolution,” in Current Protocols in Bioinformatics, chapter 1, John Wiley & Sons, 2012.
- M. J. Landrum, J. M. Lee, G. R. Riley et al., “ClinVar: public archive of relationships among sequence variation and human phenotype,” Nucleic Acids Research, vol. 42, no. 1, pp. D980–D985, 2014.
- I. F. A. C. Fokkema, P. E. M. Taschner, G. C. P. Schaafsma, J. Celli, J. F. J. Laros, and J. T. den Dunnen, “LOVD v.2.0: the next generation in gene variant databases,” Human Mutation, vol. 32, no. 5, pp. 557–563, 2011.
- Online Mendelian Inheritance in Man OM-NIoGM, Johns Hopkins University, Baltimore, Md, USA, World Wide Web, OMIM, 2013, http://omim.org/.
- D. J. McCarthy, P. Humburg, A. Kanapin et al., “Choice of transcripts and software has a large effect on variant annotation,” Genome Medicine, vol. 6, article 26, 2014.
- S. B. Ng, E. H. Turner, P. D. Robertson et al., “Targeted capture and massively parallel sequencing of 12 human exomes,” Nature, vol. 461, no. 7261, pp. 272–276, 2009.
- M. M. Alsaadi, A. M. Erzurumluoglu, S. Rodriguez et al., “Nonsense mutation in coiled-coil domain containing 151 gene (CCDC151) causes primary ciliary dyskinesia,” Human Mutation, vol. 35, no. 12, pp. 1446–1448, 2014.
- P. C. Ng and S. Henikoff, “Predicting the effects of amino acid substitutions on protein function,” Annual Review of Genomics and Human Genetics, vol. 7, pp. 61–80, 2006.
- P. C. Ng and S. Henikoff, “SIFT: predicting amino acid changes that affect protein function,” Nucleic Acids Research, vol. 31, no. 13, pp. 3812–3814, 2003.
- I. A. Adzhubei, S. Schmidt, L. Peshkin et al., “A method and server for predicting damaging missense mutations,” Nature Methods, vol. 7, no. 4, pp. 248–249, 2010.
- H. A. Shihab, J. Gough, D. N. Cooper et al., “Predicting the functional, molecular, and phenotypic consequences of amino acid substitutions using hidden Markov models,” Human Mutation, vol. 34, no. 1, pp. 57–65, 2013.
- B. Li, V. G. Krishnan, M. E. Mort et al., “Automated inference of molecular mechanisms of disease from amino acid substitutions,” Bioinformatics, vol. 25, no. 21, pp. 2744–2750, 2009.
- A. González-Pérez and N. López-Bigas, “Improving the assessment of the outcome of nonsynonymous SNVs with a consensus deleteriousness score, Condel,” The American Journal of Human Genetics, vol. 88, no. 4, pp. 440–449, 2011.
- J. Thusberg, A. Olatubosun, and M. Vihinen, “Performance of mutation pathogenicity prediction methods on missense variants,” Human Mutation, vol. 32, no. 4, pp. 358–368, 2011.
- F.-O. Desmet, D. Hamroun, M. Lalande, G. Collod-Bëroud, M. Claustres, and C. Béroud, “Human Splicing Finder: an online bioinformatics tool to predict splicing signals,” Nucleic Acids Research, vol. 37, no. 9, article e67, 2009.
- W. McLaren, B. Pritchard, D. Rios, Y. Chen, P. Flicek, and F. Cunningham, “Deriving the consequences of genomic variants with the ensembl API and SNP effect predictor,” Bioinformatics, vol. 26, no. 16, pp. 2069–2070, 2010.
- H. A. Shihab, J. Gough, D. N. Cooper, I. N. M. Day, and T. R. Gaunt, “Predicting the functional consequences of cancer-associated amino acid substitutions,” Bioinformatics, vol. 29, no. 12, pp. 1504–1510, 2013.
- E. Capriotti and R. B. Altman, “A new disease-specific machine learning approach for the prediction of cancer-causing missense variants,” Genomics, vol. 98, no. 4, pp. 310–317, 2011.
- Z. E. Sauna and C. Kimchi-Sarfaty, “Understanding the contribution of synonymous mutations to human disease,” Nature Reviews Genetics, vol. 12, no. 10, pp. 683–691, 2011.
- O. J. Buske, A. Manickaraj, S. Mital, P. N. Ray, and M. Brudno, “Identification of deleterious synonymous variants in human genomes,” Bioinformatics, vol. 29, no. 15, pp. 1843–1850, 2013.
- I. N. M. Day, “dbSNP in the detail and copy number complexities,” Human Mutation, vol. 31, no. 1, pp. 2–4, 2010.
- L. Musumeci, J. W. Arthur, F. S. G. Cheung, A. Hoque, S. Lippman, and J. K. V. Reichardt, “Single Nucleotide Differences (SNDs) in the dbSNP database may lead to errors in genotyping and haplotyping studies,” Human Mutation, vol. 31, no. 1, pp. 67–73, 2010.
- P. N. Inglis, K. A. Boroevich, and M. R. Leroux, “Piecing together a ciliome,” Trends in Genetics, vol. 22, no. 9, pp. 491–500, 2006.
- A. Franceschini, D. Szklarczyk, S. Frankild et al., “STRING v9.1: protein-protein interaction networks, with increased coverage and integration,” Nucleic Acids Research, vol. 41, no. 1, pp. D808–D815, 2013.
- M. Kanehisa and S. Goto, “KEGG: kyoto encyclopedia of genes and genomes,” Nucleic Acids Research, vol. 28, no. 1, pp. 27–30, 2000.
- P. Yue, E. Melamud, and J. Moult, “SNPs3D: candidate gene and SNP selection for association studies,” BMC Bioinformatics, vol. 7, article 166, 2006.
- A. Marques-Pinheiro, M. Marduel, J.-P. Rabès et al., “A fourth locus for autosomal dominant hypercholesterolemia maps at 16q22.1,” European Journal of Human Genetics, vol. 18, no. 11, pp. 1236–1242, 2010.
- M.-P. Audrézet, J.-M. Chen, O. Raguénès et al., “Genomic rearrangements in the CFTR gene: extensive allelic heterogeneity and diverse mutational mechanisms,” Human Mutation, vol. 23, no. 4, pp. 343–357, 2004.
- X. L. Zheng and J. E. Sadler, “Pathogenesis of thrombotic microangiopathies,” Annual Review of Pathology: Mechanisms of Disease, vol. 3, pp. 249–277, 2008.
- N. L. M. Sobreira, E. T. Cirulli, D. Avramopoulos et al., “Whole-genome sequencing of a single proband together with linkage analysis identifies a Mendelian disease gene,” PLoS Genetics, vol. 6, no. 6, Article ID e1000991, 2010.
- R. Norio, “Finnish disease heritage I: characteristics, causes, background,” Human Genetics, vol. 112, no. 5-6, pp. 441–456, 2003.
- R. Norio, “Finnish Disease Heritage II: population prehistory and genetic roots of Finns,” Human Genetics, vol. 112, no. 5-6, pp. 457–469, 2003.
- R. Norio, “The finnish disease heritage III: the individual diseases,” Human Genetics, vol. 112, no. 5-6, pp. 470–526, 2003.
- C. G. Woods, J. Cox, K. Springell et al., “Quantification of homozygosity in consanguineous individuals with autosomal recessive disease,” American Journal of Human Genetics, vol. 78, no. 5, pp. 889–896, 2006.
- A. L. Williams, N. Patterson, J. Glessner, H. Hakonarson, and D. Reich, “Phasing of many thousands of genotyped samples,” The American Journal of Human Genetics, vol. 91, no. 2, pp. 238–251, 2012.
- S. R. Browning and B. L. Browning, “Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering,” The American Journal of Human Genetics, vol. 81, no. 5, pp. 1084–1097, 2007.
- M. N. Blumenthal, “Genetic, epigenetic, and environmental factors in asthma and allergy,” Annals of Allergy, Asthma & Immunology, vol. 108, no. 2, pp. 69–73, 2012.
- J. Kettunen, T. Tukiainen, A.-P. Sarin et al., “Genome-wide association study identifies multiple loci influencing human serum metabolite levels,” Nature Genetics, vol. 44, no. 3, pp. 269–276, 2012.
- O. Combarros, M. Cortina-Borja, A. D. Smith, and D. J. Lehmann, “Epistasis in sporadic Alzheimer's disease,” Neurobiology of Aging, vol. 30, no. 9, pp. 1333–1349, 2009.
- S. Farooqi, H. Rau, J. Whitehead, and S. O'Rahilly, “ob gene mutations and human obesity,” Proceedings of the Nutrition Society, vol. 57, no. 3, pp. 471–475, 1998.
- N. Homer, B. Merriman, and S. F. Nelson, “Local alignment of two-base encoded DNA sequence,” BMC Bioinformatics, vol. 10, article 175, 2009.
- B. Langmead and S. L. Salzberg, “Fast gapped-read alignment with Bowtie 2,” Nature Methods, vol. 9, no. 4, pp. 357–359, 2012.
- N. O. Stitziel, A. Kiezun, and S. R. Sunyaev, “Computational and statistical approaches to analyzing variants identified by exome sequencing,” Genome Biology, vol. 12, article 227, 2011.
- H. Li and R. Durbin, “Fast and accurate short read alignment with Burrows-Wheeler transform,” Bioinformatics, vol. 25, no. 14, pp. 1754–1760, 2009.
- H. Li, J. Ruan, and R. Durbin, “Mapping short DNA sequencing reads and calling variants using mapping quality scores,” Genome Research, vol. 18, no. 11, pp. 1851–1858, 2008.
- R. Li, C. Yu, Y. Li et al., “SOAP2: an improved ultrafast tool for short read alignment,” Bioinformatics, vol. 25, no. 15, pp. 1966–1967, 2009.
- S. Q. Le and R. Durbin, “SNP detection and genotyping from low-coverage sequencing data on multiple diploid samples,” Genome Research, vol. 21, no. 6, pp. 952–960, 2011.
- A. R. Quinlan, D. A. Stewart, M. P. Strömberg, and G. T. Marth, “Pyrobayes: an improved base caller for SNP discovery in pyrosequences,” Nature Methods, vol. 5, no. 2, pp. 179–181, 2008.
- R. Li, Y. Li, X. Fang et al., “SNP detection for massively parallel whole-genome resequencing,” Genome Research, vol. 19, no. 6, pp. 1124–1132, 2009.
- D. Li, Y. Guo, H. Shao et al., “Genetic diversity, molecular phylogeny and selection evidence of the silkworm mitochondria implicated by complete resequencing of 41 genomes,” BMC Evolutionary Biology, vol. 10, no. 1, article 81, 2010.
- S. Li, S. Wang, Q. Deng et al., “Identification of genome-wide variations among three elite restorer lines for hybrid-rice,” PLoS ONE, vol. 7, no. 2, Article ID e30952, 2012.
- V. Boeva, T. Popova, K. Bleakley et al., “Control-FREEC: a tool for assessing copy number and allelic content using next-generation sequencing data,” Bioinformatics, vol. 28, no. 3, pp. 423–425, 2012.
- D. Challis, J. Yu, U. S. Evani et al., “An integrative variant analysis suite for whole exome next-generation sequencing data,” BMC Bioinformatics, vol. 13, article 8, 2012.
- B. Howie, J. Marchini, and M. Stephens, “Genotype imputation with thousands of genomes,” G3: Genes, Genomes, Genetics, vol. 1, no. 6, pp. 457–470, 2011.
- Y. Li, C. J. Willer, J. Ding, P. Scheet, and G. R. Abecasis, “MaCH: using sequence and genotype data to estimate haplotypes and unobserved genotypes,” Genetic Epidemiology, vol. 34, no. 8, pp. 816–834, 2010.
- N.-L. Sim, P. Kumar, J. Hu, S. Henikoff, G. Schneider, and P. C. Ng, “SIFT web server: predicting effects of amino acid substitutions on proteins,” Nucleic Acids Research, vol. 40, no. 1, pp. W452–W457, 2012.
- P. Kumar, S. Henikoff, and P. C. Ng, “Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm,” Nature Protocols, vol. 4, no. 7, pp. 1073–1081, 2009.
- E. V. Davydov, D. L. Goode, M. Sirota, G. M. Cooper, A. Sidow, and S. Batzoglou, “Identifying a high fraction of the human genome to be under selective constraint using GERP++,” PLoS Computational Biology, vol. 6, no. 12, Article ID e1001025, 2010.
- G. M. Cooper, D. L. Goode, S. B. Ng et al., “Single-nucleotide evolutionary constraint scores highlight disease-causing mutations,” Nature Methods, vol. 7, no. 4, pp. 250–251, 2010.
- G. M. Cooper, E. A. Stone, G. Asimenos, E. D. Green, S. Batzoglou, and A. Sidow, “Distribution and intensity of constraint in mammalian genomic sequence,” Genome Research, vol. 15, no. 7, pp. 901–913, 2005.
- K. S. Pollard, M. J. Hubisz, K. R. Rosenbloom, and A. Siepel, “Detection of nonneutral substitution rates on mammalian phylogenies,” Genome Research, vol. 20, no. 1, pp. 110–121, 2010.
- Y. Bromberg and B. Rost, “SNAP: predict effect of non-synonymous polymorphisms on function,” Nucleic Acids Research, vol. 35, no. 11, pp. 3823–3835, 2007.
- L. Conde, J. M. Vaquerizas, H. Dopazo et al., “PupaSuite: finding functional single nucleotide polymorphisms for large-scale genotyping purposes,” Nucleic Acids Research, vol. 34, pp. W621–W625, 2006.
- J. Reumers, J. Schymkowitz, J. Ferkinghoff-Borg, F. Stricher, L. Serrano, and F. Rousseau, “SNPeffect: a database mapping molecular phenotypic effects of human non-synonymous coding SNPs,” Nucleic Acids Research, vol. 33, pp. D527–D532, 2005.
- B. Reva, Y. Antipin, and C. Sander, “Predicting the functional impact of protein mutations: application to cancer genomics,” Nucleic Acids Research, vol. 39, no. 17, article e118, 2011.
- P. D. Thomas, A. Kejariwal, M. J. Campbell et al., “PANTHER: a browsable database of gene products organized by biological function, using curated protein family and subfamily classification,” Nucleic Acids Research, vol. 31, pp. 334–341, 2003.
- H. Mi, Q. Dong, A. Muruganujan, P. Gaudet, S. Lewis, and P. D. Thomas, “PANTHER version 7: improved phylogenetic trees, orthologs and collaboration with the Gene Ontology Consortium,” Nucleic Acids Research, vol. 38, no. 1, Article ID gkp1019, pp. D204–D210, 2009.
- R. Calabrese, E. Capriotti, P. Fariselli, P. L. Martelli, and R. Casadio, “Functional annotations improve the predictive score of human disease-related mutations in proteins,” Human Mutation, vol. 30, no. 8, pp. 1237–1244, 2009.
- L. Bao, M. Zhou, and Y. Cui, “nsSNPAnalyzer: identifying disease-associated nonsynonymous single nucleotide polymorphisms,” Nucleic Acids Research, vol. 33, no. 2, pp. W480–W482, 2005.
- E. Capriotti, R. Calabrese, and R. Casadio, “Predicting the insurgence of human genetic diseases associated to single point protein mutations with support vector machines and evolutionary information,” Bioinformatics, vol. 22, no. 22, pp. 2729–2734, 2006.
- V. Ramensky, P. Bork, and S. Sunyaev, “Human non-synonymous SNPs: server and survey,” Nucleic Acids Research, vol. 30, no. 17, pp. 3894–3900, 2002.
- C. Ferrer-Costa, J. L. Gelpí, L. Zamakola, I. Parraga, X. de la Cruz, and M. Orozco, “PMUT: a web-based tool for the annotation of pathological mutations on proteins,” Bioinformatics, vol. 21, no. 14, pp. 3176–3178, 2005.
Copyright © 2015 A. Mesut Erzurumluoglu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.