Two silent substitutions in the PDHA1 gene cause exon 5 skipping by disruption of a putative exonic splicing enhancer

Long-term use of investigational β-Hydroxybutyrate salts in children with multiple acyl-CoA dehydrogenase or pyruvate dehydrogenase deficiency

Several disorders of energy metabolism have been treated with exogenous ketone bodies. The benefit of this treatment is best documented in multiple acyl-CoA dehydrogenase deficiency (MADD) (MIM#231680). One might also expect ketone bodies to help in other disorders with impaired ketogenesis or in conditions that profit from a ketogenic diet. Here, we report the use of a novel preparation of dextro-β-hydroxybutyrate (D-βHB) salts in two cases of MADD and one case of pyruvate dehydrogenase (PDH) deficiency (MIM#312170). The two patients with MADD had previously been on a racemic mixture of D- and L‑sodium hydroxybutyrate. Patient #1 found D-βHB more palatable, and the change in formulation corrected hypernatraemia in patient #2. The patient with PDH deficiency was on a ketogenic diet but had not previously been given hydroxybutyrate. In this case, the addition of D-βHB improved ketosis. We conclude that NHS101 is a good candidate for further clinical studies in this group of diseases of inborn errors of metabolism.

Phenotype-driven approaches to enhance variant prioritization and diagnosis of rare disease

Rare disease diagnostics and disease gene discovery have been revolutionized by whole-exome and genome sequencing but identifying the causative variant(s) from the millions in each individual remains challenging. The use of deep phenotyping of patients and reference genotype-phenotype knowledge, alongside variant data such as allele frequency, segregation, and predicted pathogenicity, has proved an effective strategy to tackle this issue. Here we review the numerous tools that have been developed to automate this approach and demonstrate the power of such an approach on several thousand diagnosed cases from the 100,000 Genomes Project. Finally, we discuss the challenges that need to be overcome if we are going to improve detection rates and help the majority of patients that still remain without a molecular diagnosis after state-of-the-art genomic interpretation.

Code inside the codon: The role of synonymous mutations in regulating splicing machinery and its impact on disease

In eukaryotes, precise pre-mRNA processing, including alternative splicing, is essential to carry out the intricate protein translation process. Both point mutations (that alter the translated protein sequence) and synonymous mutations (that do not alter the translated protein sequence) are capable of affecting the splicing process. Synonymous mutations are known to affect gene expression via altering mRNA stability, mRNA secondary structure, splicing processes, and translational kinetics. In higher eukaryotes, precise splicing is regulated by three weakly conserved cis-elements, 5' and 3' splice sites and the branch site. Many other cis-acting elements (exonic/intronic splicing enhancers and silencers) and trans-acting splicing factors (serine and arginine-rich proteins and heterogeneous nuclear ribonucleoproteins) have also been found to enhance or suppress the splicing process. The appearance of synonymous mutations in cis-acting elements can alter the splicing process by changing the binding pattern of splicing factors to exonic splicing enhancers or silencer motifs. This results in exon skipping, intron retention, and various other forms of alternative splicing, eventually leading to the emergence of a wide range of diseases. The focus of this review is to elucidate the role of synonymous mutations and their impact on abnormal splicing mechanisms. Further, this study highlights the function of synonymous mutation in mediating abnormal splicing in cancer and development of X-linked, and autosomal inherited diseases.

ECHS1 disease in two unrelated families of Samoan descent: Common variant - rare disorder

Mutations in the short-chain enoyl-CoA hydratase (SCEH) gene, ECHS1, cause a rare autosomal recessive disorder of valine catabolism. Patients usually present with developmental delay, regression, dystonia, feeding difficulties, and abnormal MRI with bilateral basal ganglia involvement. We present clinical, biochemical, molecular, and functional data for four affected patients from two unrelated families of Samoan descent with identical novel compound heterozygous mutations. Family 1 has three affected boys while Family 2 has an affected daughter, all with clinical and MRI findings of Leigh syndrome and intermittent episodes of acidosis and ketosis. WES identified a single heterozygous variant in ECHS1 at position c.832G > A (p.Ala278Thr). However, western blot revealed significantly reduced ECHS1 protein for all affected family members. Decreased SCEH activity in fibroblasts and a mild increase in marker metabolites in urine further supported ECHS1 as the underlying gene defect. Additional investigations at the DNA (aCGH, WGS) and RNA (qPCR, RT-PCR, RNA-Seq, RNA-Array) level identified a silent, common variant at position c.489G > A (p.Pro163=) as the second mutation. This substitution, present at high frequency in the Samoan population, is associated with decreased levels of normally spliced mRNA. To our understanding, this is the first report of a novel, hypomorphic allele c.489G > A (p.Pro163=), associated with SCEH deficiency.

Potential role of stress-induced gluconeogenesis in disease aggravation and mortality in pyruvate dehydrogenase deficiency: A case-based hypothesis

Pyruvate dehydrogenase (PDH) deficiency is an inherited metabolic disorder caused by a defect in any subunit of the pyruvate dehydrogenase complex (PDHC), which has an essential role in glucose metabolism. The causes of disease progression in PDH deficiency are not fully understood yet. Based on repeated observations of a patient with PDH deficiency at our center, we hypothesized that stress-induced gluconeogenesis contributes to rapid exacerbation of the disease. This link has not been established previously.

Massive parallel sequencing identifies RAPSN and PDHA1 mutations causing fetal akinesia deformation sequence

INTRODUCTION

Fetal akinesia deformation sequence (FADS) or arthrogryposis multiplex congenita (AMC) is characterized by clinical ambiguity and genetic heterogeneity, hampering genetic diagnosis via traditional sequencing methods. Next generation sequencing (NGS) of all known disease-causing genes offers an elegant solution to identify the genetic etiology of AMC/FADS in a diagnostic setting.

METHODS

An in-house developed disease-associated gene panel was conducted in two unrelated fetuses with FADS. First, a de novo analysis was performed on the entire disease-associated gene panel. If no pathogenic mutation was identified, analysis of variants retained in a specific subpanel with arthrogryposis/fetal akinesia-causing genes was performed.

RESULTS

In the first family, FADS relates to a homozygous c.484G > A (p.Glu162Lys) mutation in the gene RAPSN. The second case concerns a sporadic patient with brain anomalies and arthrogryposis due to a de novo hemizygous c.498C > T splice-site mutation in the pyruvate dehydrogenase-alpha 1 (PDHA1) gene.

DISCUSSION

NGS facilitated genetic diagnosis, and hence genetic counseling, for both families with AMC/FADS. Biallelic RAPSN mutations typically result in congenital myasthenia syndrome, or occasionally in FADS. This is the first report attributing the RAPSN mutation c.484G > A, identified in a homozygous state in patient 1, to FADS. The second patient represents the first case of AMC due to a PDHA1 mutation, advocating that pyruvate dehydrogenase deficiency should be considered in the differential diagnosis of fetal akinesia. This study illustrates the relevance of a disease-associated-gene panel as a diagnostic tool in pregnancies complicated by this genetically heterogeneous condition.

Synonymous Somatic Variants in Human Cancer Are Not Infamous: A Plea for Full Disclosure in Databases and Publications

Single-nucleotide variants (SNVs) are the most frequent genetic changes found in human cancer. Most driver alterations are missense and nonsense variants localized in the coding region of cancer genes. Unbiased cancer genome sequencing shows that synonymous SNVs (sSNVs) can be found clustered in the coding regions of several cancer oncogenes or tumor suppressor genes suggesting purifying selection. sSNVs are currently underestimated, as they are usually discarded during analysis. Furthermore, several public databases do not display sSNVs, which can lead to analytical bias and the false assumption that this mutational event is uncommon. Recent progress in our understanding of the deleterious consequences of these sSNVs for RNA stability and protein translation shows that they can act as strong drivers of cancer, as demonstrated for several cancer genes such as TP53 or BCL2L12. It is therefore essential that sSNVs be properly reported and analyzed in order to provide an accurate picture of the genetic landscape of the cancer genome.

Decoding mechanisms by which silent codon changes influence protein biogenesis and function

SCOPE

Synonymous codon usage has been a focus of investigation since the discovery of the genetic code and its redundancy. The occurrences of synonymous codons vary between species and within genes of the same genome, known as codon usage bias. Today, bioinformatics and experimental data allow us to compose a global view of the mechanisms by which the redundancy of the genetic code contributes to the complexity of biological systems from affecting survival in prokaryotes, to fine tuning the structure and function of proteins in higher eukaryotes. Studies analyzing the consequences of synonymous codon changes in different organisms have revealed that they impact nucleic acid stability, protein levels, structure and function without altering amino acid sequence. As such, synonymous mutations inevitably contribute to the pathogenesis of complex human diseases. Yet, fundamental questions remain unresolved regarding the impact of silent mutations in human disorders. In the present review we describe developments in this area concentrating on mechanisms by which synonymous mutations may affect protein function and human health.

PURPOSE

This synopsis illustrates the significance of synonymous mutations in disease pathogenesis. We review the different steps of gene expression affected by silent mutations, and assess the benefits and possible harmful effects of codon optimization applied in the development of therapeutic biologics.

PHYSIOLOGICAL AND MEDICAL RELEVANCE

Understanding mechanisms by which synonymous mutations contribute to complex diseases such as cancer, neurodegeneration and genetic disorders, including the limitations of codon-optimized biologics, provides insight concerning interpretation of silent variants and future molecular therapies.

Somatic mosaicism for a novel PDHA1 mutation in a male with severe pyruvate dehydrogenase complex deficiency

Pyruvate dehydrogenase complex (PDC) deficiencies are mostly due to mutations in the X-linked PDHA1 gene. Males with hemizygous PDHA1 mutations are clinically more severely affected, while those with mosaic PDHA1 mutations may manifest milder phenotypes. We report a patient harboring a novel, mosaic missense PDHA1 mutation, c.523G > A (p.A175T), with a severe clinical presentation of congenital microcephaly, significant brain abnormalities, persistent seizures, profound developmental delay, and failure to thrive. We review published cases of PDHA1 mosaicism.

Phenylbutyrate increases pyruvate dehydrogenase complex activity in cells harboring a variety of defects

Objective

Deficiency of pyruvate dehydrogenase complex (PDHC) is the most common genetic disorder leading to lactic acidosis. PDHC deficiency is genetically heterogenous and most patients have defects in the X-linked E1-α gene but defects in the other components of the complex encoded by PDHB, PDHX, DLAT, DLD genes or in the regulatory enzyme encoded by PDP1 have also been found. Phenylbutyrate enhances PDHC enzymatic activity in vitro and in vivo by increasing the proportion of unphosphorylated enzyme through inhibition of pyruvate dehydrogenase kinases and thus, has potential for therapy of patients with PDHC deficiency. In the present study, we investigated response to phenylbutyrate of multiple cell lines harboring all known gene defects resulting in PDHC deficiency.

Methods

Fibroblasts of patients with PDHC deficiency were studied for their enzyme activity at baseline and following phenylbutyrate incubation. Drug responses were correlated with genotypes and protein levels by Western blotting.

Results

Large deletions affecting PDHA1 that result in lack of detectable protein were unresponsive to phenylbutyrate, whereas increased PDHC activity was detected in most fibroblasts harboring PDHA1 missense mutations. Mutations affecting the R349-α residue were directed to proteasome degradation and were consistently unresponsive to short-time drug incubation but longer incubation resulted in increased levels of enzyme activity and protein that may be due to an additional effect of phenylbutyrate as a molecular chaperone.

Interpretation

PDHC enzyme activity was enhanced by phenylbutyrate in cells harboring missense mutations in PDHB, PDHX, DLAT, DLD, and PDP1 genes. In the prospect of a clinical trial, the results of this study may allow prediction of in vivo response in patients with PDHC deficiency harboring a wide spectrum of molecular defects.

Cauli: a mouse strain with an Ift140 mutation that results in a skeletal ciliopathy modelling Jeune syndrome

Cilia are architecturally complex organelles that protrude from the cell membrane and have signalling, sensory and motility functions that are central to normal tissue development and homeostasis. There are two broad categories of cilia; motile and non-motile, or primary, cilia. The central role of primary cilia in health and disease has become prominent in the past decade with the recognition of a number of human syndromes that result from defects in the formation or function of primary cilia. This rapidly growing class of conditions, now known as ciliopathies, impact the development of a diverse range of tissues including the neural axis, craniofacial structures, skeleton, kidneys, eyes and lungs. The broad impact of cilia dysfunction on development reflects the pivotal position of the primary cilia within a signalling nexus involving a growing number of growth factor systems including Hedgehog, Pdgf, Fgf, Hippo, Notch and both canonical Wnt and planar cell polarity. We have identified a novel ENU mutant allele of Ift140, which causes a mid-gestation embryonic lethal phenotype in homozygous mutant mice. Mutant embryos exhibit a range of phenotypes including exencephaly and spina bifida, craniofacial dysmorphism, digit anomalies, cardiac anomalies and somite patterning defects. A number of these phenotypes can be attributed to alterations in Hedgehog signalling, although additional signalling systems are also likely to be involved. We also report the identification of a homozygous recessive mutation in IFT140 in a Jeune syndrome patient. This ENU-induced Jeune syndrome model will be useful in delineating the origins of dysmorphology in human ciliopathies.

Skeletal ciliopathies are an emerging field of human disease in which skeletal birth defects arise due to abnormal communication between cells. This failure in communication arises following mutation in components of the primary cilia, a hair-like structure present on every cell. The skeletal ciliopathies are debilitating and in severe cases lead to death in early infancy. However, the mechanisms by which these malformations come about remains unclear. Mouse models are often used to delineate the causes of human birth defects and we have identified a model that mimics one of these conditions known as Jeune syndrome. It is the first mouse model with a mutation in the Ift140 gene, and these mice exhibit phenotypes that are often seen in this set of human syndromes. We have complimented this study with the discovery of a patient that presents with Jeune Syndrome resulting from mutation of human IFT140. This model will allow us to explore the role of IFT140 and the primary cilia in normal human development and provide insight into the field of human skeletal ciliopathies.

The spectrum of pyruvate dehydrogenase complex deficiency: clinical, biochemical and genetic features in 371 patients

CONTEXT

Pyruvate dehydrogenase complex (PDC) deficiency is a genetic mitochondrial disorder commonly associated with lactic acidosis, progressive neurological and neuromuscular degeneration and, usually, death during childhood. There has been no recent comprehensive analysis of the natural history and clinical course of this disease.

OBJECTIVE

We reviewed 371 cases of PDC deficiency, published between 1970 and 2010, that involved defects in subunits E1α and E1β and components E1, E2, E3 and the E3 binding protein of the complex.

DATA SOURCES AND EXTRACTION

English language peer-reviewed publications were identified, primarily by using PubMed and Google Scholar search engines.

RESULTS

Neurodevelopmental delay and hypotonia were the commonest clinical signs of PDC deficiency. Structural brain abnormalities frequently included ventriculomegaly, dysgenesis of the corpus callosum and neuroimaging findings typical of Leigh syndrome. Neither gender nor any clinical or neuroimaging feature differentiated the various biochemical etiologies of the disease. Patients who died were younger, presented clinically earlier and had higher blood lactate levels and lower residual enzyme activities than subjects who were still alive at the time of reporting. Survival bore no relationship to the underlying biochemical or genetic abnormality or to gender.

CONCLUSIONS

Although the clinical spectrum of PDC deficiency is broad, the dominant clinical phenotype includes presentation during the first year of life; neurological and neuromuscular degeneration; structural lesions revealed by neuroimaging; lactic acidosis and a blood lactate:pyruvate ratio ≤ 20.

The spectrum of pyruvate dehydrogenase complex deficiency: clinical, biochemical and genetic features in 371 patients

CONTEXT

Pyruvate dehydrogenase complex (PDC) deficiency is a genetic mitochondrial disorder commonly associated with lactic acidosis, progressive neurological and neuromuscular degeneration and, usually, death during childhood. There has been no recent comprehensive analysis of the natural history and clinical course of this disease.

OBJECTIVE

We reviewed 371 cases of PDC deficiency, published between 1970 and 2010, that involved defects in subunits E1α and E1β and components E1, E2, E3 and the E3 binding protein of the complex.

DATA SOURCES AND EXTRACTION

English language peer-reviewed publications were identified, primarily by using PubMed and Google Scholar search engines.

RESULTS

Neurodevelopmental delay and hypotonia were the commonest clinical signs of PDC deficiency. Structural brain abnormalities frequently included ventriculomegaly, dysgenesis of the corpus callosum and neuroimaging findings typical of Leigh syndrome. Neither gender nor any clinical or neuroimaging feature differentiated the various biochemical etiologies of the disease. Patients who died were younger, presented clinically earlier and had higher blood lactate levels and lower residual enzyme activities than subjects who were still alive at the time of reporting. Survival bore no relationship to the underlying biochemical or genetic abnormality or to gender.

CONCLUSIONS

Although the clinical spectrum of PDC deficiency is broad, the dominant clinical phenotype includes presentation during the first year of life; neurological and neuromuscular degeneration; structural lesions revealed by neuroimaging; lactic acidosis and a blood lactate:pyruvate ratio ≤20.

Molecular characterization of 82 patients with pyruvate dehydrogenase complex deficiency. Structural implications of novel amino acid substitutions in E1 protein

BACKGROUND

Pyruvate dehydrogenase complex (PDHc) deficiencies are an important cause of primary lactic acidosis. Most cases result from mutations in the X-linked gene for the pyruvate dehydrogenase E1α subunit (PDHA1) while a few cases result from mutations in genes for E1β (PDHB), E2 (DLAT), E3 (DLD) and E3BP (PDHX) subunits or PDH-phosphatase (PDP1).

AIM

To report molecular characterization of 82 PDHc-deficient patients and analyze structural effects of novel missense mutations in PDHA1.

METHODS

PDHA1 variations were investigated first, by exon sequencing using a long range PCR product, gene dosage assay and cDNA analysis. Mutation scanning in PDHX, PDHB, DLAT and DLD cDNAs was further performed in unsolved cases. Novel missense mutations in PDHA1 were located on the tridimensional model of human E1 protein to predict their possible functional consequences.

RESULTS

PDHA1 mutations were found in 30 girls and 35 boys. Three large rearrangements, including two contiguous gene deletion syndrome were identified. Novel missense, frameshift and splicing mutations were also delineated and a nonsense mutation in a mosaic male. Mutations p.Glu75Ala, p.Arg88Ser, p.Arg119Trp, p.Gly144Asp, p.Pro217Arg, p.Arg235Gly, p.Tyr243Cys, p.Tyr243Ser, p.Arg245Gly, p.Pro250Leu, p.Gly278Arg, p.Met282Val, p.Gly298Glu in PDHA1 were predicted to impair active site channel conformation or subunit interactions. Six out of the seven patients with PDHB mutations displayed the recurrent p.Met101Val mutation; 9 patients harbored PDHX mutations and one patient DLD mutations.

CONCLUSION

We provide an efficient stepwise strategy for mutation screening in PDHc genes and expand the growing list of PDHA1 mutations analyzed at the structural level.

A novel point mutation within the EDA gene causes an exon dropping in mature RNA in Holstein Friesian cattle breed affected by X-linked anhidrotic ectodermal dysplasia

Background

X-linked anhidrotic ectodermal dysplasia is a disorder characterized by abnormal development of tissues and organs of ectodermal origin caused by mutations in the EDA gene. The bovine EDA gene encodes the ectodysplasin A, a membrane protein expressed in keratinocytes, hair follicles and sweat glands, which is involved in the interactions between cell and cell and/or cell and matrix. Four mutations causing ectodermal dysplasia in cattle have been described so far.

Results

We identified a new single nucleotide polymorphism (SNP) at the 9^th base of exon 8 in the EDA gene in two calves of Holstein Friesian cattle breed affected by ectodermal dysplasia. This SNP is located in the exonic splicing enhancer (ESEs) recognized by SRp40 protein. As a consequence, the spliceosome machinery is no longer able to recognize the sequence as exonic and causes exon skipping. The mutation determines the deletion of the entire exon (131 bp) in the RNA processing, causing a severe alteration of the protein structure and thus the disease.

Conclusion

We identified a mutation, never described before, that changes the regulation of alternative splicing in the EDA gene and causes ectodermal dysplasia in cattle. The analysis of the SNP allows the identification of carriers that can transmit the disease to the offspring. This mutation can thus be exploited for a rational and efficient selection of unequivocally healthy cows for breeding.

MRI features of 4 female patients with pyruvate dehydrogenase E1 alpha deficiency

Pyruvate dehydrogenase complex is a key intramitochondrial multienzyme complex required for the conversion of pyruvate to acetyl-CoA. Most patients with pyruvate dehydrogenase deficiency have a defect in the E1 alpha subunit, associated with mutations in the PDHA1 gene. In this report, we submit detailed magnetic resonance images in 4 affected female patients with PDHA1 mutations who had with severe cortical atrophy, dilated ventricles, and an incomplete corpus callosum. In one of these patients, the magnetic resonance imaging pattern prompted molecular diagnostic testing when enzymatic testing was normal. We underscore that this constellation of features, which may be misdiagnosed as periventricular leukomalacia, illustrates a pattern highly suggestive of a deficiency of pyruvate dehydrogenase E1 alpha in female patients and should trigger appropriate diagnostic investigations.

A novel mutation (c.951C>T) in an exonic splicing enhancer results in exon 10 skipping in the human mitochondrial acetoacetyl-CoA thiolase gene

Mitochondrial acetoacetyl-CoA thiolase (T2) deficiency is an inherited disorder affecting isoleucine catabolism and ketone body metabolism. A Japanese female developed a severe ketoacidotic attack at the age of 7 months. Urinary organic acid analysis showed elevated excretion of 2-methyl-3-hydroxybutyrate but not tiglylglycine. She was diagnosed as having T2 deficiency by enzyme assay using fibroblasts. Mutation analysis revealed a compound heterozygote of c.556G>T(D186Y) and c.951C>T(D317D). Since c.951C>T does not cause amino acid change, we performed cDNA analysis and found that exon 10 skipping had occurred in the c.951C>T allele. A computer search using an ESE finder showed that an exonic splicing enhancer sequence, SF2/ASF, was located in CTGA(951)CGC. We hypothesized that the exonic splicing enhancer is necessary for accurate splicing since the first nucleotide of exon 10 is C, which weakens the splice acceptor site of intron 9. We made a mini gene construct including exon 9-truncated intron 9-exon 10-truncated intron 10-exon 11 for a splicing experiment. We also made three mutant constructs which alter the SF2/ASF site (947C>T, 951C>T, 952G>A). An min-gene splicing experiment clearly showed that exon 10 skipping was induced in all three mutant constructs. Moreover, additional substitution of G for C at the first nucleotide of exon 10 resulted in normal splicing in these three mutants. These results confirmed that c.951C>T diminished the effect of the exonic splicing enhancer and caused exon 10 skipping.

Independent contribution of common CFTR variants to chronic pancreatitis

OBJECTIVE

We have assessed whether CFTR gene has a major impact on chronic pancreatitis (CP) pathogenesis than that provided by the CFTR mutations. For this aim, we have evaluated clinical parameters, CFTR mutations, and 3 potential regulatory CFTR variants (coding single-nucleotide polymorphisms): c.1540A>G, c.2694T>G, and c.4521G>A.

METHODS

CFTR gene analysis was performed in a cohort of 136 CP patients and 93 controls from Spanish population using current scanning techniques (single-strand conformation polymorphism/heteroduplex, denaturing gradient gel electrophoresis, and denaturing high-performance liquid chromatography) and direct sequencing.

RESULTS

A higher frequency of CFTR mutations were observed in patients (39%) than in controls (15%; P < or = 0.001), differences being mostly attributable to the prevalence of the cystic fibrosis (CF)-causing mutations (P = 0.009). The analysis of variants has shown statistically significant differences between patients and controls for c.4521G>A (Pcorrected = 0.036). Furthermore, the multi-marker analysis revealed that the 1540A;2694G;4521A (AGA) haplotype was more prevalent in CP than controls (Pcorrected = 0.042). Remarkably, this association was unrelated to CF-causing mutations (P = 0.006).

CONCLUSIONS

Our results corroborate the higher susceptibility of CF carriers to CP and, furthermore, suggest that the AGA haplotype could contribute to an increased risk in the development of CP irrespective of other CF-causing mutations.

Genomic features defining exonic variants that modulate splicing

A comparative analysis of SNPs and their exonic and intronic environments identifies the features predictive of splice affecting variants.

Background

Single point mutations at both synonymous and non-synonymous positions within exons can have severe effects on gene function through disruption of splicing. Predicting these mutations in silico purely from the genomic sequence is difficult due to an incomplete understanding of the multiple factors that may be responsible. In addition, little is known about which computational prediction approaches, such as those involving exonic splicing enhancers and exonic splicing silencers, are most informative.

Results

We assessed the features of single-nucleotide genomic variants verified to cause exon skipping and compared them to a large set of coding SNPs common in the human population, which are likely to have no effect on splicing. Our findings implicate a number of features important for their ability to discriminate splice-affecting variants, including the naturally occurring density of exonic splicing enhancers and exonic splicing silencers of the exon and intronic environment, extensive changes in the number of predicted exonic splicing enhancers and exonic splicing silencers, proximity to the splice junctions and evolutionary constraint of the region surrounding the variant. By extending this approach to additional datasets, we also identified relevant features of variants that cause increased exon inclusion and ectopic splice site activation.

Conclusions

We identified a number of features that have statistically significant representation among exonic variants that modulate splicing. These analyses highlight putative mechanisms responsible for splicing outcome and emphasize the role of features important for exon definition. We developed a web-tool, Skippy, to score coding variants for these relevant splice-modulating features.

Identification of novel mutations and sequence variation in the Zellweger syndrome spectrum of peroxisome biogenesis disorders

Peroxisome biogenesis disorders (PBD) are a heterogeneous group of autosomal recessive neurodegenerative disorders that affect multiple organ systems. Approximately 80% of PBD patients are classified in the Zellweger syndrome spectrum (PBD-ZSS). Mutations in the PEX1, PEX6, PEX10, PEX12, or PEX26 genes are found in approximately 90% of PBD-ZSS patients. Here, we sequenced the coding regions and splice junctions of these five genes in 58 PBD-ZSS cases previously subjected to targeted sequencing of a limited number of PEX gene exons. In our cohort, 71 unique sequence variants were identified, including 18 novel mutations predicted to disrupt protein function and 2 novel silent variants. We identified 4 patients who had two deleterious mutations in one PEX gene and a third deleterious mutation in a second PEX gene. For two such patients, we conducted cell fusion complementation analyses to identify the defective gene responsible for aberrant peroxisome assembly. Overall, we provide empirical data to estimate the relative fraction of disease-causing alleles that occur in the coding and splice junction sequences of these five PEX genes and the frequency of cases where mutations occur in multiple PEX genes. This information is beneficial for efforts aimed at establishing rapid and sensitive clinical diagnostics for PBD-ZSS patients and interpreting the results from these genetic tests.

Characterization of a novel splicing variant in the RAPTOR gene

The mammalian target of rapamycin (mTOR) plays an essential role in the regulation of cell growth, proliferation and apoptosis. Raptor, the regulatory associated protein of mTOR, is an important member in this signaling pathway. In the present report,we identified and characterized a novel splicing variant of this gene, RAPTOR v2, in which exons 14-17, 474 bp in total, are omitted from the mRNA. This deletion does not change the open reading frame, but causes a nearly complete absence of HEAT repeats, which were shown to be involved in the binding of mTOR substrates. Real time PCR performed on 48 different human tissues demonstrated the ubiquitous presence of this splice variant. Quantification of mRNA levels in lymphoblastoid cell lines (LCL) from 56 unrelated HapMap individuals revealed that the expression of this splicing form is quite variable. One synonymous SNP, rs2289759 in exon 14, was predicted by ESEfinder to cause a significant gain/loss of SRp55 and/or SF2/ASF binding sites, and thus potentially influence splicing. This prediction was confirmed by linear regression analysis between the ratio of RAPTOR v2 to total RAPTOR mRNA levels and the SNP genotype in the above 56 individuals (r=0.281 and P=0.036). Moreover, the functional evaluation indicated that this splicing isoform is expected to retain the ability to bind mTOR, but is unlikely to bind mTOR substrates, hence affecting signal transduction and further cell proliferation.