Pathogenicity Of Missense Variants Research Articles

REVEL: An Ensemble Method for Predicting the Pathogenicity of Rare Missense Variants

The vast majority of coding variants are rare, and assessment of the contribution of rare variants to complex traits is hampered by low statistical power and limited functional data. Improved methods for predicting the pathogenicity of rare coding variants are needed to facilitate the discovery of disease variants from exome sequencing studies. We developed REVEL (rare exome variant ensemble learner), an ensemble method for predicting the pathogenicity of missense variants on the basis of individual tools: MutPred, FATHMM, VEST, PolyPhen, SIFT, PROVEAN, MutationAssessor, MutationTaster, LRT, GERP, SiPhy, phyloP, and phastCons. REVEL was trained with recently discovered pathogenic and rare neutral missense variants, excluding those previously used to train its constituent tools. When applied to two independent test sets, REVEL had the best overall performance (p < 10-12) as compared to any individual tool and seven ensemble methods: MetaSVM, MetaLR, KGGSeq, Condel, CADD, DANN, and Eigen. Importantly, REVEL also had the best performance for distinguishing pathogenic from rare neutral variants with allele frequencies <0.5%. The area under the receiver operating characteristic curve (AUC) for REVEL was 0.046-0.182 higher in an independent test set of 935 recent SwissVar disease variants and 123,935 putatively neutral exome sequencing variants and 0.027-0.143 higher in an independent test set of 1,953 pathogenic and 2,406 benign variants recently reported in ClinVar than the AUCs for other ensemble methods. We provide pre-computed REVEL scores for all possible human missense variants to facilitate the identification of pathogenic variants in the sea of rare variants discovered as sequencing studies expand in scale.

Arginine-Glycine Amidinotransferase Deficiency and Functional Characterization of Missense Variants in GATM.

Arginine-glycine amidinotransferase (GATM) deficiency is an autosomal-recessive disorder caused by pathogenic variants in GATM. Clinical features include intellectual disability, hypotonia, and myopathy. Due to normal neurodevelopment in asymptomatic individuals on creatine monotherapy, GATM deficiency is a good candidate for newborn screening. To determine the carrier frequency of GATM deficiency, we performed functional characterization of rare missense variants in GATM reported as heterozygous in the Exome Variant Server database. To assess phenotype and genotype correlation, we developed a clinical severity scoring system. Two patients with mild phenotype had a nonsense missense variant. Severe phenotype was present in patients with missense as well as truncating variants. There seems to be no phenotype and genotype correlation. We cloned a novel GATM transcript. We found seven missense variants retaining 0% of wild-type GATM activity indicating putative pathogenicity. Based on our study results, high Genomic Evolutionary Rate Profiling conservation score, conserved amino acid substitution in species, and low allele frequency in exome databases would be the most sensitive in silico analysis tools to predict pathogenicity of missense variants. We present first study of the functional characterization of missense variants in GATM as well as clinical severity score of patients with GATM deficiency.

Common Pitfalls and Novel Opportunities for Predicting Variant Pathogenicity

The prediction of missense variant pathogenicity is normally performed using analyses of multiple sequence alignments optionally augmented with analyses of the (predicted) protein structure. The most straightforward way, though, is to search the literature to see whether this variant has already been described. Variant data from homologous proteins are also valuable because mutations in a homologous protein often have similar effects as mutations at the equivalent residues of the protein of interest. Transferring variant data seems trivial but is seriously hampered by the fact that homologous residue positions have different numbers in different species. This problem is even bigger when to proteins have such low sequence identities that they can no longer be aligned based on their sequences only and their structures need to be compared to align them accurately. The protein superfamily analysis software suite 3DM solves these problems, because 3DM is a system that combines high quality structure based multiple sequence alignments in which aligned residues have the same number, with all published mutant and variant data for human and all other species. We have used 3DM to analyze nine human proteins for which many disease-related variants are known. This study reveals that mutation data can be transferred even between very distant homologous proteins. Thus, protein superfamily information systems, such as 3DM, offer a wealth of unused information that can be used in the analysis of human variants.

A pilot study to estimate incidence of guanidinoacetate methyltransferase deficiency in newborns by direct sequencing of the GAMT gene

BackgroundGAMT deficiency is an autosomal recessive disorder of creatine biosynthesis causing developmental delays or intellectual disability in untreated patients as a result of irreversible brain damage occurring prior to diagnosis. Normal neurodevelopmental outcome has been reported in patients treated from neonatal period highlighting the importance of early treatment. MethodsFive hundred anonymized newborns from the National Newborn Screening Program of The Netherlands were included into this pilot study. Direct sequencing of the coding region of the GAMT gene was applied following DNA extraction. The disease causing nature of novel missense variants in the GAMT gene was studied by overexpression studies. GAA and creatine was measured in blood dot spots. ResultsWe detected two carriers, one with a known common (c.327G>A) and one with a novel mutation (c.297_309dup (p.Arg105Glyfs*) in the GAMT gene. The estimated incidence of GAMT deficiency was 1:250,000. We also detected five novel missense variants. Overexpression of these variants in GAMT deficient fibroblasts did restore GAMT activity and thus all were considered rare, but not disease causing variants including the c.131G>T (p.Arg44Leu) variant. Interestingly, this variant was predicted to be pathogenic by in silico analysis. The variants were included in the Leiden Open Variation Database (LOVD) database (www.LOVD.nl/GAMT). The average GAA level was 1.14μmol/L±0.45 standard deviations. The average creatine level was 408μmol/L±106. The average GAA/creatine ratio was 2.94±0.136. ConclusionThe estimated incidence of GAMT deficiency is 1:250,000 newborns based on our pilot study. The newborn screening for GAMT deficiency should be implemented to identify patients at the asymptomatic stage to achieve normal neurodevelopmental outcome for this treatable neurometabolic disease. Biochemical investigations including GAA, creatine and GAMT enzyme activity measurements are essential to confirm the diagnosis of GAMT deficiency. According to availability, all missense variants can be assessed functionally, as in silico prediction analysis of missense variants is not sufficient to confirm the pathogenicity of missense variants.

QDPR gene mutation and clinical follow-up in Chinese patients with dihydropteridine reductase deficiency.

This study aimed to investigate the mutation spectrum of the QDPR gene, to determine the effect of mutations on dihydropteridine reductase (DHPR) structure/function, to discuss the potential genotypephenotype correlation, and to evaluate the clinical outcome of Chinese patients after treatment. Nine DHPR-deficient patients were enrolled in this study and seven of them underwent neonatal screening. QDPR gene mutations were analyzed and confirmed by routine methods. The potential pathogenicity of missense variants was analyzed using Clustal X, PolyPhen program and Swiss-PDB Viewer 4.04_OSX software, respectively. The clinical outcomes of the patients were evaluated after long-term treatment. In 10 mutations of the 9 patients, 4 were novel mutations (G20V, V86D, G130S and A175R), 4 were reported by us previously, and 2 known mutations were identified. R221X was a hotspot mutation (27.7%) in our patients. Eight missense mutations probably had damage to protein. Six patients in this series were treated with a good control of phenylalanine level. The height and weight of the patients were normal at the age of 4 months to 7.5 years. Four patients, who underwent a neonatal screening and were treated early, showed a normal mental development. In 2 patients diagnosed late, neurological symptoms were significantly improved. The mutation spectrum of the QDPR gene is different in the Chinese population. Most mutations are related to severe phenotype. The determination of DHPR activity should be performed in patients with hyperphenylalaninemia. DHPR-deficient patients who were treated below the age of 2 months may have a near normal mental development.

Genetic and bioinformatics analysis of four novel GCK missense variants detected in Caucasian families with GCK-MODY phenotype.

Heterozygous loss-of-function mutations in the glucokinase (GCK) gene cause maturity-onset diabetes of the young (MODY) subtype GCK (GCK-MODY/MODY2). GCK sequencing revealed 16 distinct mutations (13 missense, 1 nonsense, 1 splice site, and 1 frameshift-deletion) co-segregating with hyperglycaemia in 23 GCK-MODY families. Four missense substitutions (c.718A>G/p.Asn240Asp, c.757G>T/p.Val253Phe, c.872A>C/p.Lys291Thr, and c.1151C>T/p.Ala384Val) were novel and a founder effect for the nonsense mutation (c.76C>T/p.Gln26*) was supposed. We tested whether an accurate bioinformatics approach could strengthen family-genetic evidence for missense variant pathogenicity in routine diagnostics, where wet-lab functional assays are generally unviable. In silico analyses of the novel missense variants, including orthologous sequence conservation, amino acid substitution (AAS)-pathogenicity predictors, structural modeling and splicing predictors, suggested that the AASs and/or the underlying nucleotide changes are likely to be pathogenic. This study shows how a careful bioinformatics analysis could provide effective suggestions to help molecular-genetic diagnosis in absence of wet-lab validations.

Classifying MMR Variants: Time for Revised Nomenclature in Lynch Syndrome

Inactivating germline mutations in DNA mismatch repair (MMR) genes are diagnostic for Lynch syndrome. However, the clinical significance of missense variants is uncertain. A threshold level of compromised MLH1 expression, correlating with greater protein instability and MMR functional defect, has been identified to help classify the pathogenicity of missense variants.

Prediction of Pathogenicity of Missense Variants in Arrhythmogenic Right Ventricular Cardiomyopathy

Background Pathogenicity of truncating or splice site mutations in arrhythmogenic right ventricular cardiomyopathy (ARVC) patients has been widely accepted. However, pathogenicity of missense variants remains controversial. We assessed pathogenicity of ARVC-related missense variants and associated phenotype. Methods Pathogenic mutations were found in 571 patients (239 families, age 43 ± 18 years) from US/Dutch ARVC cohorts by screening of desmosomal (PKP2, DSC2, DSG2, DSP, JUP) and/or nondesmosomal (PLN, TMEM43) genes. Patients with >1 mutation or variant were excluded. Missense variants were considered pathogenic when (1) algorithms SIFT and PolyPhen-2 predicted pathogenicity (SIFT: tolerance index score ≤0.05, PolyPhen-2: pathogenicity score >0.90) and (2) minor allele frequency in 1000 genomes was Results A pathogenic missense mutation was found in 86/571 patients (15%). Frequently occurring mutation was PKP2 c.2386T>C (p.Cys796Arg; 34/86). Survival curves (P = .065) and, more importantly, sustained arrhythmia (ventricular tachycardia/fibrillation or appropriate ICD therapy) free survival curves (P = .810) of patients with a missense mutation were similar to those with pathogenic non-missense (330/571 truncating, 129/571 splice site) mutations (Figure). Conclusions ARVC-related missense mutations have a similar pathogenic character as non-missense mutations regarding survival and substantial arrhythmic burden. Pathogenicity of missense mutations can be correctly predicted by combined criteria from SIFT, PolyPhen-2, and 1000 genomes. Download : Download high-res image (120KB) Download : Download full-size image

PON-P: Integrated predictor for pathogenicity of missense variants

High-throughput sequencing data generation demands the development of methods for interpreting the effects of genomic variants. Numerous computational methods have been developed to assess the impact of variations because experimental methods are unable to cope with both the speed and volume of data generation. To harness the strength of currently available predictors, the Pathogenic-or-Not-Pipeline (PON-P) integrates five predictors to predict the probability that nonsynonymous variations affect protein function and may consequently be disease related. Random forest methodology-based PON-P shows consistently improved performance in cross-validation tests and on independent test sets, providing ternary classification and statistical reliability estimate of results. Applied to missense variants in a melanoma cancer cell line, PON-P predicts variants in 17 genes to affect protein function. Previous studies implicate nine of these genes in the pathogenesis of various forms of cancer. PON-P may thus be used as a first step in screening and prioritizing variants to determine deleterious ones for further experimentation.

Using SIFT and PolyPhen to Predict Loss-of-Function and Gain-of-Function Mutations

The interpretation of novel missense variants is a challenge with increasing numbers of such variants being identified and a responsibility to report the findings in the context of all available scientific evidence. Various in silico bioinformatic tools have been developed that predict the likely pathogenicity of missense variants; however, their utility within the diagnostic setting requires further investigation. The aim of our study was to test the predictive value of two of these tools, sorting intolerant from tolerant (SIFT) and polymorphism phenotyping (PolyPhen), in a set of 141 missense variants (131 pathogenic, 8 benign) identified in the ABCC8, GCK, and KCNJ11 genes. Sixty-six of the mutations caused a gain of protein function, while 67 were loss-of-function mutations. The evolutionary conservation at each residue was also investigated using multiple sequence alignments from the UCSC genome browser. The sensitivity of SIFT and PolyPhen was reasonably high (69% and 68%, respectively), but their specificity was low (13% and 16%). Both programs were significantly better at predicting loss-of-function mutations than gain-of-function mutations (SIFT, p = 0.001; PolyPhen, p < or = 0.0001). The most reliable method for assessing the likely pathogenicity of a missense variant was to investigate the degree of conservation at the affected residue. Eighty-eight percent of the mutations affected highly conserved residues, while all of the benign variants occurred at residues that were polymorphic across multiple species. Although SIFT and PolyPhen may be useful in prioritizing changes that are likely to cause a loss of protein function, their low specificity means that their predictions should be interpreted with caution and further evidence to support/refute pathogenicity should be sought before reporting novel missense changes.

PKHD1 sequence variations in 78 children and adults with autosomal recessive polycystic kidney disease and congenital hepatic fibrosis

PKHD1, the gene mutated in autosomal recessive polycystic kidney disease (ARPKD)/congenital hepatic fibrosis (CHF), is an exceptionally large and complicated gene that consists of 86 exons and has a number of alternatively spliced transcripts. Its longest open reading frame contains 67 exons that encode a 4074 amino acid protein called fibrocystin or polyductin. The phenotypes caused by PKHD1 mutations are similarly complicated, ranging from perinatally-fatal PKD to CHF presenting in adulthood with mild kidney disease. To date, more than 300 mutations have been described throughout PKHD1. Most reported cohorts include a large proportion of perinatal-onset ARPKD patients; mutation detection rates vary between 42% and 87%. Here we report PKHD1 sequencing results on 78 ARPKD/CHF patients from 68 families. Differing from previous investigations, our study required survival beyond 6 months and included many adults with a CHF-predominant phenotype. We identified 77 PKHD1 variants (41 novel) including 19 truncating, 55 missense, 2 splice, and 1 small in-frame deletion. Using computer-based prediction tools (GVGD, PolyPhen, SNAP), we achieved a mutation detection rate of 79%, ranging from 63% in the CHF-predominant group to 82% in the remaining families. Prediction of the pathogenicity of missense variants will remain challenging until a functional assay is available. In the meantime, use of PKHD1 sequencing data for clinical decisions requires caution, especially when only novel or rare missense variants are identified.

No evidence of the APC D1822V missense variant's pathogenicity in Tunisian patients with sporadic colorectal cancer

Objectives Sporadic colorectal cancer is influenced by numerous single nucleotide polymorphisms (SNPs), each with minor effects on the cancer risk. This study seeks to determine whether there is any association of the I1307K, E1317Q and D1822V variants within the Adenomatous polyposis coli gene ( APC) and risk to develop colorectal cancer in Tunisian population. Methods Direct sequencing was used to investigate three SNPs in the APC in 48 Tunisian sporadic colorectal cancer cases and 63 controls. Results There was no statistically significant association between the I1307K, E1317Q and D1822V variants investigated and colorectal cancer risk. Conclusion The lack of association may show that these variants selected for this study are not involved in the colorectal carcinogenic process. Otherwise, the eventual biological effect is so little to go undetected, unless increasing the sample size.

Evaluating the clinical utility of a molecular genetic test for polycystic kidney disease

Autosomal dominant polycystic kidney disease (ADPKD) is estimated to affect 1/600–1/1000 individuals worldwide. The disease is characterized by age dependent renal cyst formation that results in kidney failure during adulthood. Although ultrasound imaging may be an adequate diagnostic tool in at risk individuals older than 30, this modality may not be sufficiently sensitive in younger individuals or for those from PKD2 families who have milder disease. DNA based assays may be indicated in certain clinical situations where imaging cannot provide a definitive clinical diagnosis. The goal of this study was to evaluate the utility of direct DNA analysis in a test sample of 82 individuals who were judged to have polycystic kidney disease by standard clinical criteria. The samples were analyzed using a commercially available assay that employs sequencing of both genes responsible for the disorder. Definite disease causing mutations were identified in 34 (∼42%) study participants. An additional 30 (∼37%) subjects had either in frame insertions/deletions, non-canonical splice site alterations or a combination of missense changes that were also judged likely to be pathogenic. We noted striking sequence variability in the PKD1 gene, with a mean of 13.1 variants per participant (range 0–60). Our results and analysis highlight the complexity of assessing the pathogenicity of missense variants particularly when individuals have multiple amino acid substitutions. We conclude that a significant fraction of ADPKD mutations are caused by amino acid substitutions that need to be interpreted carefully when utilized in clinical decision-making.