Automated splicing mutation analysis by information theory

Variant Impact Predictor database (VIPdb), version 2: trends from three decades of genetic variant impact predictors

BACKGROUND

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

RESULTS

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past three decades, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 190 VIPs, resulting in a total of 407 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

CONCLUSIONS

VIPdb version 2 summarizes 407 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications. VIPdb is available at  https://genomeinterpretation.org/vipdb.

Pharmacogenomics: A Genetic Approach to Drug Development and Therapy

The majority of the well-known pharmacogenomics research used in the medical sciences contributes to our understanding of medication interactions. It has a significant impact on treatment and drug development. The broad use of pharmacogenomics is required for the progress of therapy. The main focus is on how genes and an intricate gene system affect the body's reaction to medications. Novel biomarkers that help identify a patient group that is more or less likely to respond to a certain medication have been discovered as a result of recent developments in the field of clinical therapeutics. It aims to improve customized therapy by giving the appropriate drug at the right dose at the right time and making sure that the right prescriptions are issued. A combination of genetic, environmental, and patient variables that impact the pharmacokinetics and/or pharmacodynamics of medications results in interindividual variance in drug response. Drug development, illness susceptibility, and treatment efficacy are all impacted by pharmacogenomics. The purpose of this work is to give a review that might serve as a foundation for the creation of new pharmacogenomics applications, techniques, or strategies.

CYP46A1 variants influence Alzheimer’s disease risk and brain cholesterol metabolism

Background

Cholesterol 24S-hydroxylase (CYP46) catalyzes the conversion of cholesterol to 24S-hydroxycholesterol, the primary cerebral cholesterol elimination product. Only few gene variations in CYP46 gene (CYP46A1) have been investigated for their relevance as genetic risk factors of Alzheimer’s disease (AD) and results are contradictory.

Methods

We performed a gene variability screening in CYP46A1 and investigated the effect of gene variants on the risk of AD and on CSF levels of cholesterol and 24S-hydroxycholesterol.

Results

Two of the identified 16 SNPs in CYP46A1 influenced AD risk in our study (rs7157609: p = 0.016; rs4900442: p = 0.019). The interaction term of both SNPs was also associated with an increased risk of AD (p = 0.006). Haplotypes including both SNPs were calculated and haplotype G–C was identified to influence the risk of AD (p = 0.005). AD patients and non-demented controls, who were carriers of the G–C haplotype, presented with reduced CSF levels of 24S-hydroxycholesterol (p = 0.001) and cholesterol (p < 0.001).

Conclusion

Our results suggest that CYP46A1 gene variations might act as risk factor for AD via an influence on brain cholesterol metabolism.

Expression Changes Confirm Genomic Variants Predicted to Result in Allele-Specific, Alternative mRNA Splicing

Splice isoform structure and abundance can be affected by either noncoding or masquerading coding variants that alter the structure or abundance of transcripts. When these variants are common in the population, these nonconstitutive transcripts are sufficiently frequent so as to resemble naturally occurring, alternative mRNA splicing. Prediction of the effects of such variants has been shown to be accurate using information theory-based methods. Single nucleotide polymorphisms (SNPs) predicted to significantly alter natural and/or cryptic splice site strength were shown to affect gene expression. Splicing changes for known SNP genotypes were confirmed in HapMap lymphoblastoid cell lines with gene expression microarrays and custom designed q-RT-PCR or TaqMan assays. The majority of these SNPs (15 of 22) as well as an independent set of 24 variants were then subjected to RNAseq analysis using the ValidSpliceMut web beacon (http://validsplicemut.cytognomix.com), which is based on data from the Cancer Genome Atlas and International Cancer Genome Consortium. SNPs from different genes analyzed with gene expression microarray and q-RT-PCR exhibited significant changes in affected splice site use. Thirteen SNPs directly affected exon inclusion and 10 altered cryptic site use. Homozygous SNP genotypes resulting in stronger splice sites exhibited higher levels of processed mRNA than alleles associated with weaker sites. Four SNPs exhibited variable expression among individuals with the same genotypes, masking statistically significant expression differences between alleles. Genome-wide information theory and expression analyses (RNAseq) in tumor exomes and genomes confirmed splicing effects for 7 of the HapMap SNP and 14 SNPs identified from tumor genomes. q-RT-PCR resolved rare splice isoforms with read abundance too low for statistical significance in ValidSpliceMut. Nevertheless, the web-beacon provides evidence of unanticipated splicing outcomes, for example, intron retention due to compromised recognition of constitutive splice sites. Thus, ValidSpliceMut and q-RT-PCR represent complementary resources for identification of allele-specific, alternative splicing.

VIPdb, a genetic Variant Impact Predictor Database

Genome sequencing identifies vast number of genetic variants. Predicting these variants' molecular and clinical effects is one of the preeminent challenges in human genetics. Accurate prediction of the impact of genetic variants improves our understanding of how genetic information is conveyed to molecular and cellular functions, and is an essential step towards precision medicine. Over one hundred tools/resources have been developed specifically for this purpose. We summarize these tools as well as their characteristics, in the genetic Variant Impact Predictor Database (VIPdb). This database will help researchers and clinicians explore appropriate tools, and inform the development of improved methods. VIPdb can be browsed and downloaded at https://genomeinterpretation.org/vipdb.

RNA splicing in human disease and in the clinic

Defects at the level of the pre-mRNA splicing process represent a major cause of human disease. Approximately 15-50% of all human disease mutations have been shown to alter functioning of basic and auxiliary splicing elements. These elements are required to ensure proper processing of pre-mRNA splicing molecules, with their disruption leading to misprocessing of the pre-mRNA molecule and disease. The splicing process is a complex process, with much still to be uncovered before we are able to accurately predict whether a reported genomic sequence variant (GV) represents a splicing-associated disease mutation or a harmless polymorphism. Furthermore, even when a mutation is correctly identified as affecting the splicing process, there still remains the difficulty of providing an exact evaluation of the potential impact on disease onset, severity and duration. In this review, we provide a brief overview of splicing diagnostic methodologies, from in silico bioinformatics approaches to wet lab in vitro and in vivo systems to evaluate splicing efficiencies. In particular, we provide an overview of how the latest developments in high-throughput sequencing can be applied to the clinic, and are already changing clinical approaches.

Analysis of mutations in 7 genes associated with neuronal excitability and synaptic transmission in a cohort of children with non-syndromic infantile epileptic encephalopathy

Epileptic Encephalopathy (EE) is a heterogeneous condition in which cognitive, sensory and/or motor functions deteriorate as a consequence of epileptic activity, which consists of frequent seizures and/or major interictal paroxysmal activity. There are various causes of EE and they may occur at any age in early childhood. Genetic mutations have been identified to contribute to an increasing number of children with early onset EE which had been previously considered as cryptogenic. We identified 26 patients with Infantile Epileptic Encephalopathy (IEE) of unknown etiology despite extensive workup and without any specific epilepsy syndromic phenotypes. We performed genetic analysis on a panel of 7 genes (ARX, CDKL5, KCNQ2, PCDH19, SCN1A, SCN2A, STXBP1) and identified 10 point mutations [ARX (1), CDKL5 (3), KCNQ2 (2), PCDH19 (1), SCN1A (1), STXBP1 (2)] as well as one microdeletion involving both SCN1A and SCN2A. The high rate (42%) of mutations suggested that genetic testing of this IEE panel of genes is recommended for cryptogenic IEE with no etiology identified. These 7 genes are associated with channelopathies or synaptic transmission and we recommend early genetic testing if possible to guide the treatment strategy.

Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis

The interpretation of genomic variants has become one of the paramount challenges in the post-genome sequencing era. In this review we summarize nearly 20 years of research on the applications of information theory (IT) to interpret coding and non-coding mutations that alter mRNA splicing in rare and common diseases. We compile and summarize the spectrum of published variants analyzed by IT, to provide a broad perspective of the distribution of deleterious natural and cryptic splice site variants detected, as well as those affecting splicing regulatory sequences. Results for natural splice site mutations can be interrogated dynamically with Splicing Mutation Calculator, a companion software program that computes changes in information content for any splice site substitution, linked to corresponding publications containing these mutations. The accuracy of IT-based analysis was assessed in the context of experimentally validated mutations. Because splice site information quantifies binding affinity, IT-based analyses can discern the differences between variants that account for the observed reduced (leaky) versus abolished mRNA splicing. We extend this principle by comparing predicted mutations in natural, cryptic, and regulatory splice sites with observed deleterious phenotypic and benign effects. Our analysis of 1727 variants revealed a number of general principles useful for ensuring portability of these analyses and accurate input and interpretation of mutations. We offer guidelines for optimal use of IT software for interpretation of mRNA splicing mutations.

Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis

The interpretation of genomic variants has become one of the paramount challenges in the post-genome sequencing era. In this review we summarize nearly 20 years of research on the applications of information theory (IT) to interpret coding and non-coding mutations that alter mRNA splicing in rare and common diseases. We compile and summarize the spectrum of published variants analyzed by IT, to provide a broad perspective of the distribution of deleterious natural and cryptic splice site variants detected, as well as those affecting splicing regulatory sequences. Results for natural splice site mutations can be interrogated dynamically with Splicing Mutation Calculator, a companion software program that computes changes in information content for any splice site substitution, linked to corresponding publications containing these mutations. The accuracy of IT-based analysis was assessed in the context of experimentally validated mutations. Because splice site information quantifies binding affinity, IT-based analyses can discern the differences between variants that account for the observed reduced (leaky) versus abolished mRNA splicing. We extend this principle by comparing predicted mutations in natural, cryptic, and regulatory splice sites with observed deleterious phenotypic and benign effects. Our analysis of 1727 variants revealed a number of general principles useful for ensuring portability of these analyses and accurate input and interpretation of mutations. We offer guidelines for optimal use of IT software for interpretation of mRNA splicing mutations.

Three new PLP1 splicing mutations demonstrate pathogenic and phenotypic diversity of Pelizaeus-Merzbacher disease

Pelizaeus-Merzbacher disease is a severe X-linked disorder of central myelination caused by mutations affecting the proteolipid protein gene. We describe 3 new PLP1 splicing mutations, their effect on splicing and associated phenotypes. Mutation c.453_453+6del7insA affects the exon 3B donor splice site and disrupts the PLP1-transcript without affecting the DM20, was found in a patient with severe Pelizaeus-Merzbacher disease and in his female cousin with early-onset spastic paraparesis. Mutation c.191+1G>A causes exon 2 skipping with a frame shift, is expected to result in a functionally null allele, and was found in a patient with mild Pelizaeus-Merzbacher disease and in his aunt with late-onset spastic paraparesis. Mutation c.696+1G>A utilizes a cryptic splice site in exon 5, causes partial exon 5 skipping and in-frame deletion, and was found in an isolated patient with a severe classical Pelizaeus-Merzbacher. PLP1 splice-site mutations express a variety of disease phenotypes mediated by different molecular pathogenic mechanisms.

PMD patient mutations reveal a long-distance intronic interaction that regulates PLP1/DM20 alternative splicing

Alternative splicing of the proteolipid protein 1 gene (PLP1) produces two forms, PLP1 and DM20, due to alternative use of 5' splice sites with the same acceptor site in intron 3. The PLP1 form predominates in central nervous system RNA. Mutations that reduce the ratio of PLP1 to DM20, whether mutant or normal protein is formed, result in the X-linked leukodystrophy Pelizaeus-Merzbacher disease (PMD). We investigated the ability of sequences throughout PLP1 intron 3 to regulate alternative splicing using a splicing minigene construct transfected into the oligodendrocyte cell line, Oli-neu. Our data reveal that the alternative splice of PLP1 is regulated by a long-distance interaction between two highly conserved elements that are separated by 581 bases within the 1071-base intron 3. Further, our data suggest that a base-pairing secondary structure forms between these two elements, and we demonstrate that mutations of either element designed to destabilize the secondary structure decreased the PLP1/DM20 ratio, while swap mutations designed to restore the structure brought the PLP1/DM20 ratio to near normal levels. Sequence analysis of intron 3 in families with clinical symptoms of PMD who did not have coding-region mutations revealed mutations that segregated with disease in three families. We showed that these patient mutations, which potentially destabilize the secondary structure, also reduced the PLP1/DM20 ratio. This is the first report of patient mutations causing disease by disruption of a long-distance intronic interaction controlling alternative splicing. This finding has important implications for molecular diagnostics of PMD.

Fast Entropic Profiler: An Information Theoretic Approach for the Discovery of Patterns in Genomes

Information theory has been used for quite some time in the area of computational biology. In this paper we present a pattern discovery method, named Fast Entropic Profiler, that is based on a local entropy function that captures the importance of a region with respect to the whole genome. The local entropy function has been introduced by Vinga and Almeida in , here we discuss and improve the original formulation. We provide a linear time and linear space algorithm called Fast Entropic Profiler ( FastEP), as opposed to the original quadratic implementation. Moreover we propose an alternative normalization that can be also efficiently implemented. We show that FastEP is suitable for large genomes and for the discovery of patterns with unbounded length. FastEP is available at http://www.dei.unipd.it/~ciompin/main/FastEP.html.

Information theory applications for biological sequence analysis

Information theory (IT) addresses the analysis of communication systems and has been widely applied in molecular biology. In particular, alignment-free sequence analysis and comparison greatly benefited from concepts derived from IT, such as entropy and mutual information. This review covers several aspects of IT applications, ranging from genome global analysis and comparison, including block-entropy estimation and resolution-free metrics based on iterative maps, to local analysis, comprising the classification of motifs, prediction of transcription factor binding sites and sequence characterization based on linguistic complexity and entropic profiles. IT has also been applied to high-level correlations that combine DNA, RNA or protein features with sequence-independent properties, such as gene mapping and phenotype analysis, and has also provided models based on communication systems theory to describe information transmission channels at the cell level and also during evolutionary processes. While not exhaustive, this review attempts to categorize existing methods and to indicate their relation with broader transversal topics such as genomic signatures, data compression and complexity, time series analysis and phylogenetic classification, providing a resource for future developments in this promising area.

In silico tools for splicing defect prediction: a survey from the viewpoint of end users

RNA splicing is the process during which introns are excised and exons are spliced. The precise recognition of splicing signals is critical to this process and mutations affecting splicing comprise a considerable proportion of genetic disease etiology. Analysis of RNA samples from the patient is the most straightforward and reliable method to detect splicing defects. However, currently the technical limitation prohibits its use in routine clinical practice. In silico tools that predict potential consequences of splicing mutations may be useful in daily diagnostic activities. In this review, we provide medical geneticists with some basic insights into some of the most popular in silico tools for splicing defect prediction, from the viewpoint of end-users. Bioinformaticians in relevant areas who are working on huge datasets may also benefit from this review. Specifically, we focus on those tools whose primary goal is to predict the impact of mutations within the 5′ and 3′ splicing consensus regions: the algorithms used by different tools as well as their major advantages and disadvantages are briefly introduced; the formats of their input and output are summarized; and the interpretation, evaluation, and prospection are also discussed.

Low-frequency intermediate penetrance variants in the ROCK1 gene predispose to Tetralogy of Fallot

BACKGROUND

Epidemiological studies indicate a substantial excess familial recurrence of non-syndromic Tetralogy of Fallot (TOF), implicating genetic factors that remain largely unknown. The Rho induced kinase 1 gene (ROCK1) is a key component of the planar cell polarity signalling pathway, which plays an important role in normal cardiac development. The aim of this study was to investigate the role of genetic variation in ROCK1 on the risk of TOF.

RESULTS

ROCK1 was sequenced in a discovery cohort of 93 non-syndromic TOF probands to identify rare variants. TagSNPs were selected to capture commoner variation in ROCK1. Novel variants and TagSNPs were genotyped in a discovery cohort of 458 TOF cases and 1331 healthy controls, and positive findings were replicated in a further 209 TOF cases and 1290 healthy controls. Association between genotypes and TOF was assessed using LAMP.A rare SNP (c.807C > T; rs56085230) discovered by sequencing was associated with TOF risk (p = 0.006) in the discovery cohort. The variant was also significantly associated with the risk of TOF in the replication cohort (p = 0.018). In the combined cohorts the odds ratio for TOF was 2.61 (95% CI 1.58-4.30); p < 0.0001. The minor allele frequency of rs56085230 in the cases was 0.02, and in the controls it was 0.007. The variant accounted for 1% of the population attributable risk (PAR) of TOF. We also found significant association with TOF for an uncommon TagSNP in ROCK1, rs288979 (OR 1.64 [95% CI 1.15-2.30]; p = 1.5x10⁻⁵). The minor allele frequency of rs288979 in the controls was 0.043, and the variant accounted for 11% of the PAR of TOF. These association signals were independent of each other, providing additional internal validation of our result.

CONCLUSIONS

Low frequency intermediate penetrance (LFIP) variants in the ROCK1 gene predispose to the risk of TOF.

Interpretation, stratification and evidence for sequence variants affecting mRNA splicing in complete human genome sequences

Information theory-based methods have been shown to be sensitive and specific for predicting and quantifying the effects of non-coding mutations in Mendelian diseases. We present the Shannon pipeline software for genome-scale mutation analysis and provide evidence that the software predicts variants affecting mRNA splicing. Individual information contents (in bits) of reference and variant splice sites are compared and significant differences are annotated and prioritized. The software has been implemented for CLC-Bio Genomics platform. Annotation indicates the context of novel mutations as well as common and rare SNPs with splicing effects. Potential natural and cryptic mRNA splicing variants are identified, and null mutations are distinguished from leaky mutations. Mutations and rare SNPs were predicted in genomes of three cancer cell lines (U2OS, U251 and A431), which were supported by expression analyses. After filtering, tractable numbers of potentially deleterious variants are predicted by the software, suitable for further laboratory investigation. In these cell lines, novel functional variants comprised 6–17 inactivating mutations, 1–5 leaky mutations and 6–13 cryptic splicing mutations. Predicted effects were validated by RNA-seq analysis of the three aforementioned cancer cell lines, and expression microarray analysis of SNPs in HapMap cell lines.

Comparative in vitro and in silico analyses of variants in splicing regions of BRCA1 and BRCA2 genes and characterization of novel pathogenic mutations

Several unclassified variants (UVs) have been identified in splicing regions of disease-associated genes and their characterization as pathogenic mutations or benign polymorphisms is crucial for the understanding of their role in disease development. In this study, 24 UVs located at BRCA1 and BRCA2 splice sites were characterized by transcripts analysis. These results were used to evaluate the ability of nine bioinformatics programs in predicting genetic variants causing aberrant splicing (spliceogenic variants) and the nature of aberrant transcripts. Eleven variants in BRCA1 and 8 in BRCA2, including 8 not previously characterized at transcript level, were ascertained to affect mRNA splicing. Of these, 16 led to the synthesis of aberrant transcripts containing premature termination codons (PTCs), 2 to the up-regulation of naturally occurring alternative transcripts containing PTCs, and one to an in-frame deletion within the region coding for the DNA binding domain of BRCA2, causing the loss of the ability to bind the partner protein DSS1 and ssDNA. For each computational program, we evaluated the rate of non-informative analyses, i.e. those that did not recognize the natural splice sites in the wild-type sequence, and the rate of false positive predictions, i.e., variants incorrectly classified as spliceogenic, as a measure of their specificity, under conditions setting sensitivity of predictions to 100%. The programs that performed better were Human Splicing Finder and Automated Splice Site Analyses, both exhibiting 100% informativeness and specificity. For 10 mutations the activation of cryptic splice sites was observed, but we were unable to derive simple criteria to select, among the different cryptic sites predicted by the bioinformatics analyses, those actually used. Consistent with previous reports, our study provides evidences that in silico tools can be used for selecting splice site variants for in vitro analyses. However, the latter remain mandatory for the characterization of the nature of aberrant transcripts.

Prediction of mutant mRNA splice isoforms by information theory-based exon definition

Mutations that affect mRNA splicing often produce multiple mRNA isoforms, resulting in complex molecular phenotypes. Definition of an exon and its inclusion in mature mRNA relies on joint recognition of both acceptor and donor splice sites. This study predicts cryptic and exon-skipping isoforms in mRNA produced by splicing mutations from the combined information contents (R(i), which measures binding-site strength, in bits) and distribution of the splice sites defining these exons. The total information content of an exon (R(i),total) is the sum of the R(i) values of its acceptor and donor splice sites, adjusted for the self-information of the distance separating these sites, that is, the gap surprisal. Differences between total information contents of an exon (ΔR(i,total)) are predictive of the relative abundance of these exons in distinct processed mRNAs. Constraints on splice site and exon selection are used to eliminate nonconforming and poorly expressed isoforms. Molecular phenotypes are computed by the Automated Splice Site and Exon Definition Analysis (http://splice.uwo.ca) server. Predictions of splicing mutations were highly concordant (85.2%; n = 61) with published expression data. In silico exon definition analysis will contribute to streamlining assessment of abnormal and normal splice isoforms resulting from mutations.

A cryptic BAP1 splice mutation in a family with uveal and cutaneous melanoma, and paraganglioma

Inactivating germ line BRCA1-associated protein-1 (BAP1) mutations have recently been reported in families with uveal or cutaneous malignant melanoma (UMM, CMM), mesothelioma, and meningioma. Although apparently predisposing to a wide range of tumors, the exact tumor spectrum associated with germ line BAP1 mutations has yet to be established. Here, we report a novel germ line BAP1 splice mutation, c.1708C>G (p.Leu570fs*40), in a multiple-case Danish UMM family with a spectrum of other tumors. Whole-exome sequencing identified an apparent missense mutation of BAP1 in UMM, CMM, as well as paraganglioma, breast cancer, and suspected mesothelioma cases in the family. Bioinformatic analysis and splicing assays demonstrated that this mutation creates a strong cryptic splice donor, resulting in aberrant splicing and a truncating frameshift of the BAP1 transcript. Somatic loss of the wild-type allele was also confirmed in the UMM and paraganglioma tumors. Our findings further support BAP1 as a melanoma susceptibility gene and extend the potential predisposition spectrum to paraganglioma.

Identification of SCN1A and PCDH19 mutations in Chinese children with Dravet syndrome

BACKGROUND

Dravet syndrome is a severe form of epilepsy. Majority of patients have a mutation in SCN1A gene, which encodes a voltage-gated sodium channel. A recent study has demonstrated that 16% of SCN1A-negative patients have a mutation in PCDH19, the gene encoding protocadherin-19. Mutations in other genes account for only a very small proportion of families. TSPYL4 is a novel candidate gene within the locus 6q16.3-q22.31 identified by linkage study.

OBJECTIVE

The present study examined the mutations in epileptic Chinese children with emphasis on Dravet syndrome.

METHODS

A hundred children with severe epilepsy were divided into Dravet syndrome and non-Dravet syndrome groups and screened for SCN1A mutations by direct sequencing. SCN1A-negative Dravet syndrome patients and patients with phenotypes resembling Dravet syndrome were checked for PCDH19 and TSPYL4 mutations.

RESULTS

Eighteen patients (9 males, 9 females) were diagnosed to have Dravet syndrome. Among them, 83% (15/18) had SCN1A mutations including truncating (7), splice site (2) and missense mutations (6). The truncating/splice site mutations were associated with moderate to severe degree of intellectual disability (p<0.05). During the progression of disease, 73% (11/15) had features fitting into the diagnostic criteria of autism spectrum disorder and 53% (8/15) had history of vaccination-induced seizures. A novel PCDH19 p.D377N mutation was identified in one SCN1A-negative female patient with Dravet syndrome and a known PCDH19 p.N340S mutation in a female non-Dravet syndrome patient. The former also inherited a TSPYL4 p.G60R variant.

CONCLUSION

A high percentage of SCN1A mutations was identified in our Chinese cohort of Dravet syndrome patients but none in the rest of patients. We demonstrated that truncating/splice site mutations were linked to moderate to severe intellectual disability in these patients. A de novo PCDH19 missense mutation together with an inherited TSPYL4 missense variant were identified in a patient with Dravet syndrome.

Categorization of 77 dystrophin exons into 5 groups by a decision tree using indexes of splicing regulatory factors as decision markers

Background

Duchenne muscular dystrophy, a fatal muscle-wasting disease, is characterized by dystrophin deficiency caused by mutations in the dystrophin gene. Skipping of a target dystrophin exon during splicing with antisense oligonucleotides is attracting much attention as the most plausible way to express dystrophin in DMD. Antisense oligonucleotides have been designed against splicing regulatory sequences such as splicing enhancer sequences of target exons. Recently, we reported that a chemical kinase inhibitor specifically enhances the skipping of mutated dystrophin exon 31, indicating the existence of exon-specific splicing regulatory systems. However, the basis for such individual regulatory systems is largely unknown. Here, we categorized the dystrophin exons in terms of their splicing regulatory factors.

Results

Using a computer-based machine learning system, we first constructed a decision tree separating 77 authentic from 14 known cryptic exons using 25 indexes of splicing regulatory factors as decision markers. We evaluated the classification accuracy of a novel cryptic exon (exon 11a) identified in this study. However, the tree mislabeled exon 11a as a true exon. Therefore, we re-constructed the decision tree to separate all 15 cryptic exons. The revised decision tree categorized the 77 authentic exons into five groups. Furthermore, all nine disease-associated novel exons were successfully categorized as exons, validating the decision tree. One group, consisting of 30 exons, was characterized by a high density of exonic splicing enhancer sequences. This suggests that AOs targeting splicing enhancer sequences would efficiently induce skipping of exons belonging to this group.

Conclusions

The decision tree categorized the 77 authentic exons into five groups. Our classification may help to establish the strategy for exon skipping therapy for Duchenne muscular dystrophy.

Identification of activated cryptic 5' splice sites using structure profiles and odds measure

The activation of cryptic 5′ splice sites (5′ SSs) is often related to human hereditary diseases. The DNA-based mutation screening strategies are commonly used to recognize the cryptic 5′ SSs, because features of the local DNA sequence can influence the choice of cryptic 5′ SSs. To improve the identification of the cryptic 5′ SSs, we developed a structure-based method, named SPO (structure profiles and odds measure), which combines two parameters, the structural feature derived from hydroxyl radical cleavage pattern and odds measure, to assess the likelihood of a cryptic 5′ SS activation in competing with its paired authentic 5′ SS. Compared to the current tools for identifying activated cryptic 5′ SSs, the SPO algorithm achieves higher prediction accuracy than the other methods, including MaxEnt, MDD, Markov model, weight matrix model, Shapiro and Senapathy matrix, R_i and ΔG. In addition, the predicted ΔSPO scores from the SPO algorithm exhibited a greater degree of correlation with the strength of cryptic 5′ SS activation than that measured from the other seven methods. In conclusion, the SPO algorithm provides an optimal identification of cryptic 5′ SSs, can be applied in designing mutagenesis experiments for various splicing events and may be helpful to investigate the relationship between structural variants and human hereditary diseases.

Screening of the SOD1, FUS, TARDBP, ANG, and OPTN mutations in Korean patients with familial and sporadic ALS

About 5% of amyotrophic lateral sclerosis (ALS) cases are known to be familial (fALS) and mutations in SOD1 and other genes are found in more than 20% of fALS patients and in 2%-4% of apparently sporadic ALS (sALS) cases. However, there are few reports on the proportion of fALS and the frequency of mutations in Korean patients with ALS. We screened mutations in the SOD1, FUS, TARDBP, ANG, and OPTN genes in 258 consecutively enrolled Korean patients with ALS from October 2006 to November 2010. The frequency of fALS was estimated to be 3.5% (9/258), and mutations were identified in 88.9% (8/9) of fALS patients but only in 2.8% (7/249) of sALS patients. Seven fALS and 3 sALS patients had mutations in SOD1 gene while all the others had FUS gene. The proportion of fALS was lower than that reported in Caucasian populations but the frequency of SOD1 gene mutations in Korean fALS patients (77.8%, 7/9) was much higher than that reported in other ethnic groups. These findings might suggest that there is an ethnic difference in the proportion of fALS and the genetic background of ALS.

VAX1 mutation associated with microphthalmia, corpus callosum agenesis, and orofacial clefting: the first description of a VAX1 phenotype in humans

Vax1 and Vax2 have been implicated in eye development and the closure of the choroid fissure in mice and zebrafish. We sequenced the coding exons of VAX1 and VAX2 in 70 patients with anophthalmia/microphthalmia (A/M). In VAX1, we observed homozygosity for two successive nucleotide substitutions c.453G>A and c.454C>A, predicting p.Arg152Ser, in a proband of Egyptian origin with microphthalmia, small optic nerves, cleft lip/palate, and corpus callosum agenesis. This mutation affects an invariant residue in the homeodomain of VAX1 and was absent from 96 Egyptian controls. It is likely that the mutation results in a loss of function, as the mutation results in a phenotype similar to the Vax1 homozygous null mouse. We did not identify any mutations in VAX2. This is the first description of a phenotype associated with a VAX1 mutation in humans and establishes VAX1 as a new causative gene for A/M.

Molecular genetic characterization of SMAD signaling molecules in pulmonary arterial hypertension

Heterozygous germline mutations of BMPR2 contribute to familial clustering of pulmonary arterial hypertension (PAH). To further explore the genetic basis of PAH in isolated cases, we undertook a candidate gene analysis to identify potentially deleterious variation. Members of the bone morphogenetic protein (BMP) pathway, namely SMAD1, SMAD4, SMAD5, and SMAD9, were screened by direct sequencing for gene defects. Four variants were identified in SMADs 1, 4, and 9 among a cohort of 324 PAH cases, each not detected in a substantial control population. Of three amino acid substitutions identified, two demonstrated reduced signaling activity in vitro. A putative splice site mutation in SMAD4 resulted in moderate transcript loss due to compromised splicing efficiency. These results demonstrate the role of BMPR2 mutation in the pathogenesis of PAH and indicate that variation within the SMAD family represents an infrequent cause of the disease.

Clinical expression and new SPINK5 splicing defects in Netherton syndrome: unmasking a frequent founder synonymous mutation and unconventional intronic mutations

Netherton syndrome (NS) is a severe skin disease caused by loss-of-function mutations in SPINK5 (serine protease inhibitor Kazal-type 5) encoding the serine protease inhibitor LEKTI (lympho-epithelial Kazal type-related inhibitor). Here, we disclose new SPINK5 defects in 12 patients, who presented a clinical triad suggestive of NS with variations in inter- and intra-familial disease expression. We identified a new and frequent synonymous mutation c.891C>T (p.Cys297Cys) in exon 11 of the 12 NS patients. This mutation disrupts an exonic splicing enhancer sequence and causes out-of-frame skipping of exon 11. Haplotype analysis indicates that this mutation is a founder mutation in Greece. Two other new deep intronic mutations, c.283-12T>A in intron 4 and c.1820+53G>A in intron 19, induced partial intronic sequence retention. A new nonsense c.2557C>T (p.Arg853X) mutation was also identified. All mutations led to a premature termination codon resulting in no detectable LEKTI on skin sections. Two patients with deep intronic mutations showed residual LEKTI fragments in cultured keratinocytes. These fragments retained some functional activity, and could therefore, together with other determinants, contribute to modulate the disease phenotype. This new founder mutation, the most frequent mutation described in European populations so far, and these unusual intronic mutations, widen the clinical and molecular spectrum of NS and offer new diagnostic perspectives for NS patients.

The association of LRP5 gene polymorphisms with ankylosing spondylitis in a Chinese Han population

OBJECTIVE

To clarify the association between polymorphisms of low-density lipoprotein receptor-related protein 5 (LRP5) with ankylosing spondylitis (AS) in a Chinese Han population.

METHODS

Sixteen patients with AS were recruited for preliminary screening through gene sequencing. Then 14 single-nucleotide polymorphisms (SNP) of LRP5 were followed up in 296 patients and 170 controls.

RESULTS

Sequencing the LRP5 showed 24 SNP including 3 novel SNP [LRP5SNP1 (c.-1596T > C), LRP5SNP2 (c.3764-30G > A), and LRP5SNP3 (c.4488+74G > A)]. Validation of SNP showed that the LRP5SNP3 were associated with AS after multiple testing correction (allele P(c) = 0.0087, genotype P(c) = 0.0316, haplotype AGA, P(c) = 0.0051, OR = 2.54 and haplotype AGG, P(c) = 0.048, OR = 0.63, respectively). The SNP rs686921 was associated with male predominance in both patients with AS (p = 0.032, OR 1.54) and controls (p = 0.014, OR 1.94).

CONCLUSION

LRP5 may be involved in the pathogenesis of AS. Further study will be required to clarify the effect of LRP5 on the pathogenic mechanism of AS.

SCN1A mutational analysis in Korean patients with Dravet syndrome

OBJECTIVE

The aim of this study was to characterize the SCN1A mutation spectrum in Korean patients with Dravet syndrome.

METHODS

Twenty-nine patients diagnosed with Dravet syndrome at the Seoul National University Children's Hospital were included in the study. Direct sequencing and multiplex ligation-dependent probe amplification (MLPA) were used to identify SCN1A mutations. Mutations were classified as either truncation (nonsense and frameshift) or missense mutations.

RESULTS

Nineteen pathogenic mutations (19/29; 66%) and three unclassified variants were identified. One large deletion mutation spanning exons 1-20 was detected using MLPA. Fifteen of these 19 SCN1A mutations were novel. Eleven mutations were classified as truncations (seven frameshift and four nonsense mutations) and seven were classified as missense mutations. Truncating mutations spanned the whole span of subunits of the SCN1A protein, whereas all missense mutations were localized at either the voltage sensor (S4) or the ion pore (S5-S6) regions. Analysis according to clinical phenotype revealed that SCN1A mutations were more frequent in the classic group than in the borderline group (78% vs. 45%).

CONCLUSIONS

SCN1A mutational analysis of Korean Dravet syndrome patients resulted in the identification of 15 novel mutations, which could expand the spectrum of SCN1A mutations and confirms the current understanding of genotype-phenotype correlations.

Identification of germline alterations of the mad homology 2 domain of SMAD3 and SMAD4 from the Ontario site of the breast cancer family registry (CFR)

INTRODUCTION

A common feature of neoplastic cells is that mutations in SMADs can contribute to the loss of sensitivity to the anti-tumor effects of transforming growth factor-β (TGF-β). However, germline mutation analysis of SMAD3 and SMAD4, the principle substrates of the TGF-β signaling pathway, has not yet been conducted in breast cancer. Thus, it is currently unknown whether germline SMAD3 and SMAD4 mutations are involved in breast cancer predisposition.

METHODS

We performed mutation analysis of the highly conserved mad-homology 2 (MH2) domains for both genes in genomic DNA from 408 non-BRCA1/BRCA2 breast cancer cases and 710 population controls recruited by the Ontario site of the breast cancer family registry (CFR) using denaturing high-performance liquid chromatography (DHPLC) and direct DNA sequencing. The results were interpreted in several ways. First, we adapted nucleotide diversity analysis to quantitatively assess whether the frequency of alterations differ between the two genes. Next, in silico tools were used to predict variants' effect on domain function and mRNA splicing. Finally, 37 cases or controls harboring alterations were tested for aberrant splicing using reverse-transcription polymerase chain reaction (PCR) and real-time PCR statistical comparison of germline expressions by non-parametric Mann-Whitney test of independent samples.

RESULTS

We identified 27 variants including 2 novel SMAD4 coding variants c.1350G > A (p.Gln450Gln), and c.1701A > G (p.Ile525Val). There were no inactivating mutations even though c.1350G > A was predicted to affect exonic splicing enhancers. However, several additional findings were of note: 1) nucleotide diversity estimate for SMAD3 but not SMAD4 indicated that coding variants of the MH2 domain were more infrequent than expected; 2) in breast cancer cases SMAD3 was significantly over-expressed relative to controls (P < 0.05) while the case harboring SMAD4 c.1350G > A was associated with elevated germline expression (> 5-fold); 3) separate analysis using tissue expression data showed statistically significant over-expression of SMAD3 and SMAD4 in breast carcinomas.

CONCLUSIONS

This study shows that inactivating germline alterations in SMAD3 and SMAD4 are rare, suggesting a limited role in driving tumorigenesis. Nevertheless, aberrant germline expressions of SMAD3 and SMAD4 may be more common in breast cancer than previously suspected and offer novel insight into their roles in predisposition and/or progression of breast cancer.

Disruption of PTPRO causes childhood-onset nephrotic syndrome

Idiopathic nephrotic syndrome (INS) is a genetically heterogeneous group of disorders characterized by proteinuria, hypoalbuminemia, and edema. Because it typically results in end-stage kidney disease, the steroid-resistant subtype (SRNS) of INS is especially important when it occurs in children. The present study included 29 affected and 22 normal individuals from 17 SRNS families; genome-wide analysis was performed with Affymetrix 250K SNP arrays followed by homozygosity mapping. A large homozygous stretch on chromosomal region 12p12 was identified in one consanguineous family with two affected siblings. Direct sequencing of protein tyrosine phosphatase receptor type O (PTPRO; also known as glomerular epithelial protein-1 [GLEPP1]) showed homozygous c.2627+1G>T donor splice-site mutation. This mutation causes skipping of the evolutionarily conserved exon 16 (p.Glu854_Trp876del) at the RNA level. Immunohistochemistry with GLEPP1 antibody showed a similar staining pattern in the podocytes of the diseased and control kidney tissues. We used a highly polymorphic intragenic DNA marker-D12S1303-to search for homozygosity in 120 Turkish and 13 non-Turkish individuals in the PodoNet registry. This analysis yielded 17 candidate families, and a distinct homozygous c.2745+1G>A donor splice-site mutation in PTPRO was further identified via DNA sequencing in a second Turkish family. This mutation causes skipping of exon 19, and this introduces a premature stop codon at the very beginning of exon 20 (p.Asn888Lysfs*3) and causes degradation of mRNA via nonsense-mediated decay. Immunohistochemical analysis showed complete absence of immunoreactive PTPRO. Ultrastructural alterations, such as diffuse foot process fusion and extensive microvillus transformation of podocytes, were observed via electron microscopy in both families. The present study introduces mutations in PTPRO as another cause of autosomal-recessive nephrotic syndrome.

Preimplantation genetic diagnosis for hereditary endocrine disease

OBJECTIVE

To describe the process of preimplantation genetic diagnosis (PGD), which allows the selection of embryos without mutations for implantation, with specific application for mutations in GNAS.

METHODS

We identified a GNAS mutation in a patient with a severe form of Albright hereditary osteodystrophy and pseudohypoparathyroidism type 1a with phocomelia and performed PGD on embryos derived by in vitro fertilization in order to deliver an unaffected infant.

RESULTS

After in vitro fertilization, embryos that were homozygous normal for GNAS were identified and implanted into the mother. Ultrasonography 34 days after embryo transfer showed a viable singleton intrauterine pregnancy. A normal-appearing male infant was born at 36.5 weeks of gestation. Newborn screening revealed normal results of thyroid function tests, and a buccal swab obtained when the child was 1 year old verified normal GNAS gene sequences.

CONCLUSION

PGD is an alternative that can be offered for many genetic diseases and represents a method to decrease and potentially eliminate the transmission of severe genetic diseases. Patients with multiple endocrine neoplasia (MEN) type 2 with known RET gene mutations as well as those with other heritable disorders are candidates for PGD. Successful PGD in patients with MEN has not yet been reported and has met with some early difficulties, but it is believed that this technique will eventually be successful for MEN and other hereditary endocrine disorders.

Manitoba-oculo-tricho-anal (MOTA) syndrome is caused by mutations in FREM1

BACKGROUND

Manitoba-oculo-tricho-anal (MOTA) syndrome is a rare condition defined by eyelid colobomas, cryptophthalmos and anophthalmia/microphthalmia, an aberrant hairline, a bifid or broad nasal tip, and gastrointestinal anomalies such as omphalocele and anal stenosis. Autosomal recessive inheritance had been assumed because of consanguinity in the Oji-Cre population of Manitoba and reports of affected siblings, but no locus or cytogenetic aberration had previously been described.

METHODS AND RESULTS

This study shows that MOTA syndrome is caused by mutations in FREM1, a gene previously mutated in bifid nose, renal agenesis, and anorectal malformations (BNAR) syndrome. MOTA syndrome and BNAR syndrome can therefore be considered as part of a phenotypic spectrum that is similar to, but distinct from and less severe than, Fraser syndrome. Re-examination of Frem1(bat/bat) mutant mice found new evidence that Frem1 is involved in anal and craniofacial development, with anal prolapse, eyelid colobomas, telecanthus, a shortened snout and reduced philtral height present in the mutant mice, similar to the human phenotype in MOTA syndrome.

CONCLUSIONS

The milder phenotypes associated with FREM1 deficiency in humans (MOTA syndrome and BNAR syndrome) compared to that resulting from FRAS1 and FREM2 loss of function (Fraser syndrome) are also consistent with the less severe phenotypes resulting from Frem1 loss of function in mice. Together, Fraser, BNAR and MOTA syndromes constitute a clinically overlapping group of FRAS-FREM complex diseases.

Altered mRNA splicing in lipoprotein disorders

PURPOSE OF REVIEW

To review recent publications concerning the functional assessment on pre-mRNA splicing of genomic variants found in some monogenic dyslipidemias. Examples are derived from familial hypercholesterolemia, familial HDL deficiency/Tangier disease and familial hypobetalipoproteinemia.

RECENT FINDINGS

About 5-10% of genomic variants found in familial hypercholesterolemia, FHD/Tangier disease and familial hypobetalipoproteinemia are located in the introns of the candidate genes and are classified as splicing mutations. Although variants affecting highly conserved GT/AG dinucleotides at the splice sites are likely to be pathogenic, it is difficult to predict the effects of variants located deep in the introns. Algorithms were developed to predict the effect of these variants and to provide the rationale for functional studies. Combined in-silico and wet bench analysis revealed that some intronic variants classified as pathogenic have no effect, whereas others generated abnormal transcripts. Nucleotide substitutions at the 5' and the 3' of exons might change the splice site consensus sequence, causing splicing defects. Rare silent mutations were identified which create new splice sites within exons, with the consequent production of abnormal transcripts.

SUMMARY

Intronic variants, even if located deep in introns, as well as exonic variants could affect splicing with the formation of abnormal transcripts encoding structurally abnormal proteins.

In silico prediction of splice-affecting nucleotide variants

It appears that all types of genomic nucleotide variations can be deleterious by affecting normal pre-mRNA splicing via disruption/creation of splice site consensus sequences. As it is neither pertinent nor realistic to perform functional testing for all of these variants, it is important to identify those that could lead to a splice defect in order to restrict experimental transcript analyses to the most appropriate cases. In silico tools designed to provide this type of prediction are available. In this chapter, we present in silico splice tools integrated in the Alamut (Interactive Biosoftware) application and detail their use in routine diagnostic applications. At this time, in silico predictions are useful for variants that decrease the strength of wild-type splice sites or create a cryptic splice site. Importantly, in silico predictions are not sufficient to classify variants as neutral or deleterious: they should be used as part of the decision-making process to detect potential candidates for splicing anomalies, prompting molecular geneticists to carry out transcript analyses in a limited and pertinent number of cases which could be managed in routine settings.

Splicing mutations in glycogen-storage disease type II: evaluation of the full spectrum of mutations and their relation to patients' phenotypes

Glycogen-storage disease type II is an autosomal recessive-inherited disorder due to the deficiency of acid α-glucosidase. A large number of mutations in the acid α-glucosidase gene have been described to date. Among them, ~15% are variations that may affect mRNA splicing process. In this study, we have for the first time comprehensively reviewed the available information on splicing mutations of the acid α-glucosidase gene and we have evaluated their possible impact on the splicing process using different in silico approaches. Out of the 39 different GAA-sequence variations described, an in silico analysis using seven different programs showed that 97% of them are predicted to have an impact on the splicing process. Moreover, this analysis showed a quite good correlation between the impact of the mutation on the splicing process and the clinical phenotype. In addition, we have performed the functional characterization of three novel sequence variants found in Italian patients and still uncharacterized. Using a minigene system, we have confirmed their pathogenic nature. In conclusion, this study has shown that in silico analysis represents a useful tool to select mutations that affect the splicing process of the acid α-glucosidase gene and provides an updated picture of all this kind of mutations reported till now.

Analysis of mismatch repair gene mutations in Turkish HNPCC patients

Hereditary non-polyposis colorectal cancer (HNPCC or Lynch syndrome) is caused by the inheritance of a mutant allele of a DNA mismatch repair gene. We aimed to investigate types and frequencies of mismatch repair (MMR) gene mutations in Turkish patients with HNPCC and to identify specific biomarkers for early diagnosis of their non-symptomatic kindred's. The molecular characteristics of 28 Turkish colorectal cancer patients at high-risk for HNPCC were investigated by analysis of microsatellite instability (MSI), immunohistochemistry and methylation-specific PCR in order to select tumors for mutation analysis. Ten cases (35.7%) were classified as MSI (+). Lack of expression of the main MMR proteins was observed in MSI (+) tumors. Hypermethylation of the MLH1 promoter region was observed in one tumor. Nine Lynch syndrome cases showed novel germ-line alterations of the MMR gene: two frame-shifts (MLH1 c.1843dupC and MLH1 c.1743delG) and three missense mutations (MLH1 c.293G>C, MLH1 c.954_955delinsTA and MSH2 c.2210G>A). Unclassified variants were evaluated as likely to be pathogenic by using the in-silico analyses. In addition, the MSH2 c.2210G>A alteration could be considered as a founder mutation for the Turkish population due to its identification in five different Lynch syndrome families and absence in control group. The present study adds new information about MMR gene mutation types and their role in Lynch syndrome. This is the first detailed research on Turkish Lynch syndrome families.

An improved in silico selection of phenotype affecting polymorphisms in SLC6A4, HTR1A and HTR2A genes

OBJECTIVE

Among the experimentally assessed DNA variations in serotonin related genes, some influence physiological expression of personality and mental disorders, others alter the responses to pharmacological and/or psychotherapeutic treatments. Because of the huge number of polymorphisms lying in genes and of the great length of time necessary to perform association studies, a selection of the variations being studied is a necessary and crucial step.

METHODS

In this work we used the most updated and assessed bioinformatic tools to predict the phenotype affecting polymorphisms of the human HTR1A, HTR2A and SLC6A4 serotonin related genes. Moreover, we carried out a literature search to collect information about the recent association studies to compare it versus our prediction data.

RESULTS

Gene polymorphism analysis indicated the variations that are worth considering in the association studies in the field of psychiatry, psychology and pharmacogenomics. The literature revision allowed to show both the few well and the most not enough investigated polymorphisms.

CONCLUSIONS

Our data can be useful to select polymorphisms for new association studies, especially those not yet investigated that can be related to behaviour, mental disorders and individual treatment response.

Clinical and genetic analysis of Korean patients with congenital insensitivity to pain with anhidrosis

Congenital insensitivity to pain with anhidrosis (CIPA) is a rare autosomal recessive disease characterized by anhidrosis, insensitivity to noxious stimuli, and mental retardation. Mutations in the NTRK1 gene are associated with the pathogenesis of CIPA. In this study, we performed a clinical and genetic analysis on the NTRK1 gene in four Korean patients with CIPA. All patients had typical clinical manifestations of CIPA, including anhidrosis, recurrent fever, absent pain perception, and developmental delay. Sequencing analysis revealed one predominant mutation, c.851-33T>A, in four affected alleles and three novel mutations, including c.287+2dupT, c.2155G>A (p.Glu719Lys), and c.1218delC (p.Pro407ArgfsX), in each affected allele. For one patient, who was heterozygous for c.851-33T>A, another mutation could not be identified, suggesting that a possible hidden intronic or large genomic mutation may have been present. This study extends the spectrum of mutations in the NTRK1 gene and confirms that Korean patients with CIPA have the same genetic background as other ethnicities.