Mutation prediction by PolyPhen or functional assay, a detailed comparison of CYP27B1 missense mutations

Mutation of c.244G>T in NR5A1 gene causing 46, XY DSD by affecting RNA splicing

Objective

To identify the pathogenic mechanism of the c.244G>T mutation in NR5A1 gene found in a Chinese patient with 46, XY disorders of sex development (DSD). Subjects and methods: Genomic DNA was extracted from a Chinese 46, XY DSD patient. Targeted next-generation and Sanger sequencing were performed to investigate and validate the gene mutation causing 46, XY DSD, respectively. In silico tools were used to predict the pathogenicity of the variant. Dual luciferase reporter gene assay and minigene splicing reporter assay were used to identify the pathogenicity of the variant.

Results

A novel heterozygous variant, c.244G>T (p.Ala82Ser), in NR5A1 gene was detected in the 46, XY DSD patient. Four of five silico tools predicting pathogenicity of missense variants indicated that the variant was pathogenic. However, in vitro functional study showed that p.Ala82Ser did not affect the transcriptional activity of NR5A1. In silico tools predicting the potential splicing loci revealed that c.244G>T led to aberrant splicing of NR5A1 RNA. Minigene splicing reporter assay confirmed that c.244G>T resulted in the deletion of exon2 or deletion of 19 nucleotides in 3′ end of exon2.

Conclusions

Mutation of c.244G>T in NR5A1 results in 46, XY DSD by inducing abnormal splicing of NR5A1 RNA instead of amino acid substitution of NR5A1.

Identification of genetic variants associated with tacrolimus metabolism in kidney transplant recipients by extreme phenotype sampling and next generation sequencing

An extreme phenotype sampling (EPS) model with targeted next-generation sequencing (NGS) identified genetic variants associated with tacrolimus (Tac) metabolism in subjects from the Deterioration of Kidney Allograft Function (DeKAF) Genomics cohort which included 1,442 European Americans (EA) and 345 African Americans (AA). This study included 48 subjects separated into 4 groups of 12 (AA high, AA low, EA high, EA low). Groups were selected by the extreme phenotype of dose-normalized Tac trough concentrations after adjusting for common genetic variants and clinical factors. NGS spanned >3 Mb of 28 genes and identified 18,661 genetic variants (3,961 previously unknown). A group of 125 deleterious variants, by SIFT analysis, were associated with Tac troughs in EAs (burden test, p=0.008), CYB5R2 was associated with Tac troughs in AAs (SKAT, p=0.00079). In CYB5R2, rs61733057 (increased allele frequency in AAs) was predicted to disrupt protein function by SIFT and PolyPhen2 analysis. The variants merit further validation.

Elucidation of the genetic and epigenetic landscape alterations in RNA binding proteins in glioblastoma

RNA binding proteins (RBPs) have been implicated in cancer development. An integrated bioinformatics analysis of RBPs (n = 1756) in various datasets (n = 11) revealed several genetic and epigenetically altered events among RBPs in glioblastoma (GBM). We identified 13 mutated and 472 differentially regulated RBPs in GBM samples. Mutations in AHNAK predicted poor prognosis. Copy number variation (CNV), DNA methylation and miRNA targeting contributed to RBP differential regulation. Two sets of differentially regulated RBPs that may be implicated in initial astrocytic transformation and glioma progression were identified. We have also identified a four RBP (NOL3, SUCLG1, HERC5 and AFF3) signature, having a unique expression pattern in glioma stem-like cells (GSCs), to be an independent poor prognostic indicator in GBM. RBP risk score derived from the signature also stratified GBM into low-risk and high-risk groups with significant survival difference. Silencing NOL3, SUCLG1 and HERC5 inhibited GSC maintenance. Gene set enrichment analysis of differentially regulated genes between high-risk and low-risk underscored the importance of inflammation, EMT and hypoxia in high-risk GBM. Thus, we provide a comprehensive overview of genetic and epigenetic regulation of RBPs in glioma development and progression.

Neurogenetics in Argentina: diagnostic yield in a personalized research based clinic

As a whole neurogenetic diseases are a common group of neurological disorders. However, the recognitionand molecular diagnosis of these disorders is not always straightforward. Besides, there is a paucity of informationregarding the diagnostic yield that specialized neurogenetic clinics could obtain. We performed a prospective,observational, analytical study of the patients seen in a neurogenetic clinic at a tertiary medicalcentre to assess the diagnostic yield of a comprehensive diagnostic evaluation that included a personalizedclinical assessment along with traditional and next-generation sequencing diagnostic tests. We included a cohortof 387 patients from May 2008 to June 2014. For sub-group analysis we selected a sample of patientswhose main complaint was the presence of progressive ataxia, to whom we applied a systematic moleculardiagnostic algorithm. Overall, a diagnostic mutation was identified in 27·4% of our cohort. However, if weonly considered those patients where a molecular test could be performed, the success rate rises to 45%. Weobtained diagnostic yields of 23·5 and 57·5% in the global group of ataxic patients and in the subset of ataxicpatients with a positive family history, respectively. Thus, about a third of patients evaluated in a neurogeneticclinic could be successfully diagnosed.

Computational Pathology: A Path Ahead

CONTEXT

We define the scope and needs within the new discipline of computational pathology, a discipline critical to the future of both the practice of pathology and, more broadly, medical practice in general.

OBJECTIVE

To define the scope and needs of computational pathology.

DATA SOURCES

A meeting was convened in Boston, Massachusetts, in July 2014 prior to the annual Association of Pathology Chairs meeting, and it was attended by a variety of pathologists, including individuals highly invested in pathology informatics as well as chairs of pathology departments.

CONCLUSIONS

The meeting made recommendations to promote computational pathology, including clearly defining the field and articulating its value propositions; asserting that the value propositions for health care systems must include means to incorporate robust computational approaches to implement data-driven methods that aid in guiding individual and population health care; leveraging computational pathology as a center for data interpretation in modern health care systems; stating that realizing the value proposition will require working with institutional administrations, other departments, and pathology colleagues; declaring that a robust pipeline should be fostered that trains and develops future computational pathologists, for those with both pathology and nonpathology backgrounds; and deciding that computational pathology should serve as a hub for data-related research in health care systems. The dissemination of these recommendations to pathology and bioinformatics departments should help facilitate the development of computational pathology.

Whole exome sequencing identifies mutation of EDNRA involved in ACTH-independent macronodular adrenal hyperplasia

ACTH independent macronodular adrenal hyperplasia (AIMAH) is a rare disorder characterized by bilateral macronodular hyperplasia of the adrenal glands and increased cortisol production with subclinical or overt Cushing's syndrome. Although the family clustering of AIMAH is infrequent, we have tried our best to find such a familial affected pedigree with complete clinical information and successfully collect adrenalectomy tissue samples from two members of this family. Using whole exome sequencing and several variant prioritization strategies based on disease network analysis, we identified Endothelin receptor type A (EDNRA) Ser420Thr mutation as a causative mutation of AIMAH. EDNRA is a member of G protein coupled receptor family and is involved in cardiovascular or polycystic kidney disease. Our findings indicate that the mutation of EDNRA at S420T site should be regard as a potential AIMAH causative variation in familial and sporadic affected patients.

Single nucleotide polymorphism and its dynamics for pharmacogenomics

Pharmacogenomics is the study of how the genetic makeup determines the response to a therapeutic intervention. It has the capability to revolutionize the practice of medicine by personalized approach for treatment through the use of novel diagnostic tools. Pharmacogenomic based approaches reduce the trial-and-error approach and restrict the exposure of patients to those drugs which are not effective or are toxic for them. Single Nucleotide Polymorphisms (SNPs) hold the key in defining the risk of an individual's susceptibility to various illnesses and response to drugs. There is an ongoing process of identifying the common, biologically relevant SNPs, in particular those that are associated with the risk of disease and adverse drug reaction. The identification and characterization of these SNPs are necessary before their use as genetic tools. Most of the ongoing SNP related studies are biased deliberately towards coding regions and the data generated from them are therefore unlikely to reflect genome wide distribution of SNPs. To avoid this biasing towards the coding regions SNP, SNP consortium protocol was designed. Though, projects like the HapMap increase credibility and use of SNPs, still there are some concern like the required sample (patient) sizes, the number of SNPs required for mapping, number of association studies, the cost of SNP genotyping, and the interpretation and explanation of results are some of the challenges that surround this field.

A novel compound mutation of CYP27B1 in a Chinese family with vitamin D-dependent rickets type 1A

OBJECTIVES

Mutations in the CYP27B1 gene, which encodes vitamin D 1α-hydroxylase, are the genetic basis of vitamin D-dependent rickets type 1A (VDDR1A, MIM 264700). The aim of this study was to investigate a novel CYP27B1 mutation and its clinical manifestations.

METHODS

VDDR1A was diagnosed based on clinical presentation, a physical examination, bone characteristics on an X-ray, and laboratory results. A molecular model of the CYP27B1 protein was constructed using the SWISS-MODEL server and Swiss-PdbViewer.

RESULTS

We sequenced the CYP27B1 gene in a 5-year-old male child who presented with growth retardation and a history of frequent hand, leg, and perioral twitching since the age of 12 months. We identified a compound heterozygous mutation consisting of two missense mutations: one in exon 7 (R389C [c.1165C>T]) and one in exon 8 (R459C [c.1375C>T]). We used the wild-type CYP27B1 as a receptor and calcidiol as a ligand to predict the interaction between the R459 site and calcidiol. According to the predicted structure, the wild-type R459 residue localizes to the pocket where CYP27B1 binds to its ligand.

CONCLUSIONS

According to the Human Gene Mutation Database, the compound heterozygous mutation identified in our patient is novel and has not yet been reported in the literature. This mutation provides a new basis for further research on VDDR1A and for the development of clinical diagnostics.

Targeted exome sequencing for mitochondrial disorders reveals high genetic heterogeneity

Background

Mitochondrial disorders are difficult to diagnose due to extreme genetic and phenotypic heterogeneities.

Methods

We explored the utility of targeted next-generation sequencing for the diagnosis of mitochondrial disorders in 148 patients submitted for clinical testing. A panel of 447 nuclear genes encoding mitochondrial respiratory chain complexes, and other genes inducing secondary mitochondrial dysfunction or that cause diseases which mimic mitochondrial disorders were tested.

Results

We identified variants considered to be possibly disease-causing based on family segregation data and/or variants already known to cause disease in twelve genes in thirteen patients. Rare or novel variants of unknown significance were identified in 45 additional genes for various metabolic, genetic or neurogenetic disorders.

Conclusions

Primary mitochondrial defects were confirmed only in four patients indicating that majority of patients with suspected mitochondrial disorders are presumably not the result of direct impairment of energy production. Our results support that clinical and routine laboratory ascertainment for mitochondrial disorders are challenging due to significant overlapping non-specific clinical symptoms and lack of specific biomarkers. While next-generation sequencing shows promise for diagnosing suspected mitochondrial disorders, the challenges remain as the underlying genetic heterogeneity may be greater than suspected and it is further confounded by the similarity of symptoms with other conditions as we report here.

Novel CYP2B6 enzyme variants in a Rwandese population: functional characterization and assessment of in silico prediction tools

Cytochrome P450 CYP2B6 is a highly polymorphic enzyme that metabolizes numerous drugs, pesticides, and environmental toxins. Sequence analysis of a Rwandese population identified eight functionally uncharacterized nonsynonymous variants c.329G>T (p.G110V), c.341T>C (p.I114T), c.444G>T (p.E148D), c.548T>G (p.V183G), c.637T>C (p.F213L), c.758G>A (p.R253H), c.835G>C (p.A279P), and c.1459C>A (p.R487S), and five novel alleles termed CYP2B6*33 to CYP2B6*37 were assigned. Recombinant expression in COS-1 cells and functional characterization using the antidepressant bupropion and the antiretroviral efavirenz (EFV) as substrates demonstrated complete or almost complete loss-of-function for variants p.G110V, p.I114T, p.V183G, and p.F213L, whereas p.E148D, p.R253H, p.A279P, and p.R487S variants were functional. The data were used to assess the predictive power of eight online available functional prediction programs for amino-acid changes. Although none of the programs correctly predicted the functionality of all variants, substrate docking simulation analyses indicated similar conformational changes by all four deleterious mutations within the enzyme's active site, thus explaining lack of enzymatic function for both substrates. Because low-activity alleles of CYP2B6 are associated with impaired EFV metabolism and adverse drug response, these results are of potential utility for personalized treatment strategies in HIV/AIDS therapy.

Novel variations in the adiponectin gene (ADIPOQ) may affect distribution of oligomeric complexes

Adiponectin is an obesity related protein that mediates the risk of type 2 diabetes in obese individuals with its anti-inflammatory and insulin-sensitizing properties. To date, five functional variations have been identified in the adiponectin gene. However, these variations are rare, and fail to fully explain adiponectin variability, suggesting unidentified causal variations exist. Thus, our objective was to identify novel, potentially functional amino acid-changing variations in ADIPOQ exonic regions and relate them to oligomeric forms of adiponectin in serum. We sequenced ADIPOQ exons in 30 adolescents chosen from a school-based cohort based on serum adiponectin and insulin levels. Four coding region changes were identified: a methionine initiation skip (MIS), P32L, R55C, and Y111H, of which R55C and Y111H have been previously identified. Individuals with the novel variations and R55C had low levels of adiponectin and decreased adiponectin oligomerization compared to adolescents with similar body mass index and insulin levels. Further, bioinformatic analysis predicted putative functionality of these variations. In our study, Y111H was unrelated to total circulating adiponectin or adiponectin oligomerization. Given the disruption of adiponectin oligomerization in the individuals with MIS, P32L, and R55C coding changes, these variations may lead to increased metabolic disease risk and warrant further examination in larger cohorts.

THAP1 mutations and dystonia phenotypes: genotype phenotype correlations

THAP1 mutations have been shown to be the cause of DYT6. A number of different mutation types and locations in the THAP1 gene have been associated with a range of severity and dystonia phenotypes, but, as yet, it has been difficult to identify clear genotype phenotype patterns. Here, we screened the THAP1 gene in a further series of dystonia cases and evaluated the mutation pathogenicity in this series as well as previously reported mutations to investigate possible phenotype-genotype correlations. THAP1 mutations have been identified throughout the coding region of the gene, with the greatest concentration of variants localized to the THAP1 domain. In the additional cases analyzed here, a further two mutations were found. No obvious, indisputable genotype-phenotype correlation emerged from these data. However, we managed to find a correlation between the pathogenicity of mutations, distribution, and age of onset of dystonia. THAP1 mutations are an important cause of dystonia, but, as yet, no clear genotype-phenotype correlations have been identified. Greater mutation numbers in different populations will be important and mutation-specific functional studies will be essential to identify the pathogenicity of the various THAP1 mutations. © 2012 Movement Disorder Society

Identifying disease mutations in genomic medicine settings: current challenges and how to accelerate progress

The pace of exome and genome sequencing is accelerating, with the identification of many new disease-causing mutations in research settings, and it is likely that whole exome or genome sequencing could have a major impact in the clinical arena in the relatively near future. However, the human genomics community is currently facing several challenges, including phenotyping, sample collection, sequencing strategies, bioinformatics analysis, biological validation of variant function, clinical interpretation and validity of variant data, and delivery of genomic information to various constituents. Here we review these challenges and summarize the bottlenecks for the clinical application of exome and genome sequencing, and we discuss ways for moving the field forward. In particular, we urge the need for clinical-grade sample collection, high-quality sequencing data acquisition, digitalized phenotyping, rigorous generation of variant calls, and comprehensive functional annotation of variants. Additionally, we suggest that a 'networking of science' model that encourages much more collaboration and online sharing of medical history, genomic data and biological knowledge, including among research participants and consumers/patients, will help establish causation and penetrance for disease causal variants and genes. As we enter this new era of genomic medicine, we envision that consumer-driven and consumer-oriented efforts will take center stage, thus allowing insights from the human genome project to translate directly back into individualized medicine.

Glucokinase (GCK) mutations and their characterization in MODY2 children of southern Italy

Type 2 Maturity Onset Diabetes of the Young (MODY2) is a monogenic autosomal disease characterized by a primary defect in insulin secretion and hyperglycemia. It results from GCK gene mutations that impair enzyme activity. Between 2006 and 2010, we investigated GCK mutations in 66 diabetic children from southern Italy with suspected MODY2. Denaturing High Performance Liquid Chromatography (DHPLC) and sequence analysis revealed 19 GCK mutations in 28 children, six of which were novel: p.Glu40Asp, p.Val154Leu, p.Arg447Glyfs, p.Lys458_Cys461del, p.Glu395_Arg397del and c.580-2A>T. We evaluated the effect of these 19 mutations using bioinformatic tools such as Polymorphism Phenotyping (Polyphen), Sorting Intolerant From Tolerant (SIFT) and in silico modelling. We also conducted a functional study to evaluate the pathogenic significance of seven mutations that are among the most severe mutations found in our population, and have never been characterized: p.Glu70Asp, p.His137Asp, p.Phe150Tyr, p.Val154Leu, p.Gly162Asp, p.Arg303Trp and p.Arg392Ser. These seven mutations, by altering one or more kinetic parameters, reduced enzyme catalytic activity by >40%. All mutations except p.Glu70Asp displayed thermal-instability, indeed >50% of enzyme activity was lost at 50°C/30 min. Thus, these seven mutations play a pathogenic role in MODY2 insurgence. In conclusion, this report revealed six novel GCK mutations and sheds some light on the structure-function relationship of human GCK mutations and MODY2.

Performance of computational tools in evaluating the functional impact of laboratory-induced amino acid mutations

Summary: Site-directed mutagenesis is frequently used by scientists to investigate the functional impact of amino acid mutations in the laboratory. Over 10 000 such laboratory-induced mutations have been reported in the UniProt database along with the outcomes of functional assays. Here, we explore the performance of state-of-the-art computational tools (Condel, PolyPhen-2 and SIFT) in correctly annotating the function-altering potential of 10 913 laboratory-induced mutations from 2372 proteins. We find that computational tools are very successful in diagnosing laboratory-induced mutations that elicit significant functional change in the laboratory (up to 92% accuracy). But, these tools consistently fail in correctly annotating laboratory-induced mutations that show no functional impact in the laboratory assays. Therefore, the overall accuracy of computational tools for laboratory-induced mutations is much lower than that observed for the naturally occurring human variants. We tested and rejected the possibilities that the preponderance of changes to alanine and the presence of multiple base-pair mutations in the laboratory were the reasons for the observed discordance between the performance of computational tools for natural and laboratory mutations. Instead, we discover that the laboratory-induced mutations occur predominately at the highly conserved positions in proteins, where the computational tools have the lowest accuracy of correct prediction for variants that do not impact function (neutral). Therefore, the comparisons of experimental-profiling results with those from computational predictions need to be sensitive to the evolutionary conservation of the positions harboring the amino acid change.

Contact:s.kumar@asu.edu