SCNM1, a putative RNA splicing factor that modifies disease severity in mice

Escaping from CRISPR-Cas-mediated knockout: the facts, mechanisms, and applications

Clustered regularly interspaced short palindromic repeats and associated Cas protein (CRISPR-Cas), a powerful genome editing tool, has revolutionized gene function investigation and exhibits huge potential for clinical applications. CRISPR-Cas-mediated gene knockout has already become a routine method in research laboratories. However, in the last few years, accumulating evidences have demonstrated that genes knocked out by CRISPR-Cas may not be truly silenced. Functional residual proteins could be generated in such knockout organisms to compensate the putative loss of function, termed herein knockout escaping. In line with this, several CRISPR-Cas-mediated knockout screenings have discovered much less abnormal phenotypes than expected. How does knockout escaping happen and how often does it happen have not been systematically reviewed yet. Without knowing this, knockout results could easily be misinterpreted. In this review, we summarize these evidences and propose two main mechanisms allowing knockout escaping. To avoid the confusion caused by knockout escaping, several strategies are discussed as well as their advantages and disadvantages. On the other hand, knockout escaping also provides convenient tools for studying essential genes and treating monogenic disorders such as Duchenne muscular dystrophy, which are discussed in the end.

Human voltage-gated Na+ and K+ channel properties underlie sustained fast AP signaling

Human cortical pyramidal neurons are large, have extensive dendritic trees, and yet have unexpectedly fast input-output properties: Rapid subthreshold synaptic membrane potential changes are reliably encoded in timing of action potentials (APs). Here, we tested whether biophysical properties of voltage-gated sodium (Na+) and potassium (K+) currents in human pyramidal neurons can explain their fast input-output properties. Human Na+ and K+ currents exhibited more depolarized voltage dependence, slower inactivation, and faster recovery from inactivation compared with their mouse counterparts. Computational modeling showed that despite lower Na+ channel densities in human neurons, the biophysical properties of Na+ channels resulted in higher channel availability and contributed to fast AP kinetics stability. Last, human Na+ channel properties also resulted in a larger dynamic range for encoding of subthreshold membrane potential changes. Thus, biophysical adaptations of voltage-gated Na+ and K+ channels enable fast input-output properties of large human pyramidal neurons.

A spinal muscular atrophy modifier implicates the SMN protein in SNARE complex assembly at neuromuscular synapses

Reduced survival motor neuron (SMN) protein triggers the motor neuron disease, spinal muscular atrophy (SMA). Restoring SMN prevents disease, but it is not known how neuromuscular function is preserved. We used model mice to map and identify an Hspa8G470R synaptic chaperone variant, which suppressed SMA. Expression of the variant in the severely affected mutant mice increased lifespan >10-fold, improved motor performance, and mitigated neuromuscular pathology. Mechanistically, Hspa8G470R altered SMN2 splicing and simultaneously stimulated formation of a tripartite chaperone complex, critical for synaptic homeostasis, by augmenting its interaction with other complex members. Concomitantly, synaptic vesicular SNARE complex formation, which relies on chaperone activity for sustained neuromuscular synaptic transmission, was found perturbed in SMA mice and patient-derived motor neurons and was restored in modified mutants. Identification of the Hspa8G470R SMA modifier implicates SMN in SNARE complex assembly and casts new light on how deficiency of the ubiquitous protein causes motor neuron disease.

Mutations in SCNM1 cause orofaciodigital syndrome due to minor intron splicing defects affecting primary cilia

Orofaciodigital syndrome (OFD) is a genetically heterogeneous ciliopathy characterized by anomalies of the oral cavity, face, and digits. We describe individuals with OFD from three unrelated families having bi-allelic loss-of-function variants in SCNM1 as the cause of their condition. SCNM1 encodes a protein recently shown to be a component of the human minor spliceosome. However, so far the effect of loss of SCNM1 function on human cells had not been assessed. Using a comparative transcriptome analysis between fibroblasts derived from an OFD-affected individual harboring SCNM1 mutations and control fibroblasts, we identified a set of genes with defective minor intron (U12) processing in the fibroblasts of the affected subject. These results were reproduced in SCNM1 knockout hTERT RPE-1 (RPE-1) cells engineered by CRISPR-Cas9-mediated editing and in SCNM1 siRNA-treated RPE-1 cultures. Notably, expression of TMEM107 and FAM92A encoding primary cilia and basal body proteins, respectively, and that of DERL2, ZC3H8, and C17orf75, were severely reduced in SCNM1-deficient cells. Primary fibroblasts containing SCNM1 mutations, as well as SCNM1 knockout and SCNM1 knockdown RPE-1 cells, were also found with abnormally elongated cilia. Conversely, cilia length and expression of SCNM1-regulated genes were restored in SCNM1-deficient fibroblasts following reintroduction of SCNM1 via retroviral delivery. Additionally, functional analysis in SCNM1-retrotransduced fibroblasts showed that SCNM1 is a positive mediator of Hedgehog (Hh) signaling. Our findings demonstrate that defective U12 intron splicing can lead to a typical ciliopathy such as OFD and reveal that primary cilia length and Hh signaling are regulated by the minor spliceosome through SCNM1 activity.

Integrative Analysis of Nanopore and Illumina Sequencing Reveals Alternative Splicing Complexity in Pig Longissimus Dorsi Muscle

Alternative splicing (AS) is a key step in the post-transcriptional regulation of gene expression that can affect intramuscular fat (IMF). In this study, longissimus dorsi muscles from 30 pigs in high- and low- IMF groups were used to perform Oxford Nanopore Technologies (ONT) full-length sequencing and Illumina strand-specific RNA-seq. A total of 43,688 full-length transcripts were identified, with 4,322 novel genes and 30,795 novel transcripts. Using AStalavista, a total of 14,728 AS events were detected in the longissimus dorsi muscle. About 17.79% of the genes produced splicing isoforms, in which exon skipping was the most frequent AS event. By analyzing the expression differences of mRNAs and splicing isoforms, we found that differentially expressed mRNAs with splicing isoforms could participate in skeletal muscle development and fatty acid metabolism, which might determine muscle-related traits. SERBP1, MYL1, TNNT3, and TNNT1 were identified with multiple splicing isoforms, with significant differences in expression. AS events occurring in IFI6 and GADD45G may cause significant differences in gene expression. Other AS events, such as ONT.15153.3, may regulate the function of ART1 by regulating the expression of different transcripts. Moreover, co-expression and protein-protein interaction (PPI) analysis indicated that several genes (MRPL27, AAR2, PYGM, PSMD4, SCNM1, and HNRNPDL) may be related to intramuscular fat. The splicing isoforms investigated in our research provide a reference for the study of alternative splicing regulation of intramuscular fat deposition.

Challenges of CYP21A2 genotyping in children with 21-hydroxylase deficiency: determination of genotype-phenotype correlation using next generation sequencing in Southeastern Anatolia

BACKGROUND/PURPOSE

Although it is known that there is generally a good correlation between genotypes and phenotypes, the number of studies reporting discrepancies has recently increased, exclusively between milder genotypes and their phenotypes due to the complex nature of the CYP21A2 gene and methodological pitfalls. This study aimed to assess CYP21A2 genotyping in children with 21-hydroxylase deficiency (21-OHD) and establish their predictive genotype-phenotype correlation features using a large cohort in Southeastern Anatolia's ethnically diverse population.

METHODS

The patients were classified into three groups: salt-wasting (SW), simple virilizing (SV) and non-classical (NC). The genotypes were categorized into six groups due to residual enzyme activity: null-A-B-C-D-E. CYP21A2 genotyping was performed by sequence-specific primer and sequenced with next generation sequencing (NGS), and the expected phenotypes were compared to the observed phenotypes.

RESULTS

A total of 118 unrelated children with 21-OHD were included in this study (61% SW, 24.5% SV and 14.5% NC). The pathogenic variants were found in 79.5% of 171 mutated alleles (60.2%, 22.2%, and 17.6% in SW, SV and NC, respectively). Patient distribution based on genotype groups was as follows: null-16.1%, A-41.4%, B-6.0%, C-14.4%, E-22%). In2G was the most common pathogenic variant (33.9% of all alleles) and the most common variant in the three phenotype groups (SW-38.8%, SV-22.2% and NC-23.3%). The total genotype-phenotype correlation was 81.5%. The correlations of the null and A groups were 100% and 76.1%, respectively, while it was lower in group B and poor in group C (71.4% and 23.5%, respectively).

CONCLUSION

This study revealed that the concordance rates of the severe genotypes with their phenotypes were good, while those of the milder genotypes were poor. The discrepancies could have resulted from the complex characteristics of 21-OHD genotyping and the limitations of using NGS alone without integrating with other comprehensive methods.

Evidence in disease and non-disease contexts that nonsense mutations cause altered splicing via motif disruption

Transcripts containing premature termination codons (PTCs) can be subject to nonsense-associated alternative splicing (NAS). Two models have been evoked to explain this, scanning and splice motif disruption. The latter postulates that exonic cis motifs, such as exonic splice enhancers (ESEs), are disrupted by nonsense mutations. We employ genome-wide transcriptomic and k-mer enrichment methods to scrutinize this model. First, we show that ESEs are prone to disruptive nonsense mutations owing to their purine richness and paucity of TGA, TAA and TAG. The motif model correctly predicts that NAS rates should be low (we estimate 5–30%) and approximately in line with estimates for the rate at which random point mutations disrupt splicing (8–20%). Further, we find that, as expected, NAS-associated PTCs are predictable from nucleotide-based machine learning approaches to predict splice disruption and, at least for pathogenic variants, are enriched in ESEs. Finally, we find that both in and out of frame mutations to TAA, TGA or TAG are associated with exon skipping. While a higher relative frequency of such skip-inducing mutations in-frame than out of frame lends some credence to the scanning model, these results reinforce the importance of considering splice motif modulation to understand the etiology of PTC-associated disease.

Sodium channelopathies in neurodevelopmental disorders

The voltage-gated sodium channel α-subunit genes comprise a highly conserved gene family. Mutations of three of these genes, SCN1A, SCN2A and SCN8A, are responsible for a significant burden of neurological disease. Recent progress in identification and functional characterization of patient variants is generating new insights and novel approaches to therapy for these devastating disorders. Here we review the basic elements of sodium channel function that are used to characterize patient variants. We summarize a large body of work using global and conditional mouse mutants to characterize the in vivo roles of these channels. We provide an overview of the neurological disorders associated with mutations of the human genes and examples of the effects of patient mutations on channel function. Finally, we highlight therapeutic interventions that are emerging from new insights into mechanisms of sodium channelopathies.

Antisense oligonucleotide development for the selective modulation of CYP3A5 in renal disease

CYP3A5 is the primary CYP3A subfamily enzyme expressed in the human kidney and its aberrant expression may contribute to a broad spectrum of renal disorders. Pharmacogenetic studies have reported inconsistent linkages between CYP3A5 expression and hypertension, however, most investigators have considered CYP3A5*1 as active and CYP3A5*3 as an inactive allele. Observations of gender specific differences in CYP3A5*3/*3 protein expression suggest additional complexity in gene regulation that may underpin an environmentally responsive role for CYP3A5 in renal function. Reconciliation of the molecular mechanism driving conditional restoration of functional CYP3A5*3 expression from alternatively spliced transcripts, and validation of a morpholino-based approach for selectively suppressing renal CYP3A5 expression, is the focus of this work. Morpholinos targeting a cryptic splice acceptor created by the CYP3A5*3 mutation in intron 3 rescued functional CYP3A5 expression in vitro, and salt-sensitive cellular mechanisms regulating splicing and conditional expression of CYP3A5*3 transcripts are reported. The potential for a G-quadruplex (G4) in intron 3 to mediate restored splicing to exon 4 in CYP3A5*3 transcripts was also investigated. Finally, a proximal tubule microphysiological system (PT-MPS) was used to evaluate the safety profile of morpholinos in proximal tubule epithelial cells, highlighting their potential as a therapeutic platform for the treatment of renal disease.

Structure of the activated human minor spliceosome

The minor spliceosome mediates splicing of the rare but essential U12-type precursor messenger RNA. Here, we report the atomic features of the activated human minor spliceosome determined by cryo-electron microscopy at 2.9-angstrom resolution. The 5' splice site and branch point sequence of the U12-type intron are recognized by the U6atac and U12 small nuclear RNAs (snRNAs), respectively. Five newly identified proteins stabilize the conformation of the catalytic center: The zinc finger protein SCNM1 functionally mimics the SF3a complex of the major spliceosome, the RBM48-ARMC7 complex binds the γ-monomethyl phosphate cap at the 5' end of U6atac snRNA, the U-box protein PPIL2 coordinates loop I of U5 snRNA and stabilizes U5 small nuclear ribonucleoprotein (snRNP), and CRIPT stabilizes U12 snRNP. Our study provides a framework for the mechanistic understanding of the function of the human minor spliceosome.

Plant Cyclophilins: Multifaceted Proteins With Versatile Roles

Cyclophilins constitute a family of ubiquitous proteins that bind cyclosporin A (CsA), an immunosuppressant drug. Several of these proteins possess peptidyl-prolyl cis-trans isomerase (PPIase) activity that catalyzes the cis-trans isomerization of the peptide bond preceding a proline residue, essential for correct folding of the proteins. Compared to prokaryotes and other eukaryotes studied until now, the cyclophilin gene families in plants exhibit considerable expansion. With few exceptions, the role of the majority of these proteins in plants is still a matter of conjecture. However, recent studies suggest that cyclophilins are highly versatile proteins with multiple functionalities, and regulate a plethora of growth and development processes in plants, ranging from hormone signaling to the stress response. The present review discusses the implications of cyclophilins in different facets of cellular processes, particularly in the context of plants, and provides a glimpse into the molecular mechanisms by which these proteins fine-tune the diverse physiological pathways.

Normal audiogram but poor sensitivity to brief sounds in mice with compromised voltage-gated sodium channels (Scn8amedJ)

The Scn8a^medJ mutation of the gene for sodium channels at the nodes of Ranvier slows nerve conduction, resulting in motor abnormalities. This mutation is also associated with loss of spontaneous bursting activity in the dorsal cochlear nucleus. However initial tests of auditory sensitivity in mice homozygous for this mutation, using standard 400-ms tones, demonstrated normal hearing sensitivity. Further testing, reported here, revealed a severely compromised sensitivity to short-duration tones of 10 and 2 ms durations. Such a deficit might be expected to interfere with auditory functions that depend on rapid processing of auditory signals.

Predicting the likelihood of postoperative seizure status based on mRNA sequencing in low-grade gliomas

AIM

No comprehensive and objective methods yet exist for predicting postoperative seizure.

PATIENTS & METHODS

mRNA-seq data and corresponding postoperative seizure status of 109 low-grade glioma samples were obtained from Chinese Glioma Genome Atlas database and divided into two sets randomly. Logistic regression and receiver operating characteristic analysis with risk score method were used to develop a ten-gene prediction model.

RESULTS

Considering gene number and area under the curve of receiver operating characteristic, a ten-gene model was generated which showed an area under the curve of 0.9965 in training set. Patients with high-risk scores had higher probability of postoperative seizure compared with those with low-risk scores.

CONCLUSION

This is the first prediction model for postoperative seizures in gliomas, integrating multiple genes. Clinical application may help patients with postoperative seizure control.

mRNA processing in mutant zebrafish lines generated by chemical and CRISPR-mediated mutagenesis produces unexpected transcripts that escape nonsense-mediated decay

As model organism-based research shifts from forward to reverse genetics approaches, largely due to the ease of genome editing technology, a low frequency of abnormal phenotypes is being observed in lines with mutations predicted to lead to deleterious effects on the encoded protein. In zebrafish, this low frequency is in part explained by compensation by genes of redundant or similar function, often resulting from the additional round of teleost-specific whole genome duplication within vertebrates. Here we offer additional explanations for the low frequency of mutant phenotypes. We analyzed mRNA processing in seven zebrafish lines with mutations expected to disrupt gene function, generated by CRISPR/Cas9 or ENU mutagenesis methods. Five of the seven lines showed evidence of altered mRNA processing: one through a skipped exon that did not lead to a frame shift, one through nonsense-associated splicing that did not lead to a frame shift, and three through the use of cryptic splice sites. These results highlight the need for a methodical analysis of the mRNA produced in mutant lines before making conclusions or embarking on studies that assume loss of function as a result of a given genomic change. Furthermore, recognition of the types of adaptations that can occur may inform the strategies of mutant generation.

The recent rise of reverse genetic, gene targeting methods has allowed researchers to readily generate mutations in any gene of interest with relative ease. Should these mutations have the predicted effect on the mRNA and encoded protein, we would expect many more abnormal phenotypes than are typically being seen in reverse genetic screens. Here we set out to explore some of the reasons for this discrepancy by studying seven separate mutations in zebrafish. We present evidence that thorough cDNA sequence analysis is a key step in assessing the likelihood that a given mutation will produce hypomorphic or null alleles. This study reveals that mRNA processing in the mutant background often produces transcripts that escape nonsense-mediated decay, thereby potentially preserving gene function. By understanding the ways that cells avoid the deleterious consequences of mutations, researchers can better design reverse genetic strategies to increase the likelihood of gene disruption.

Genetics of Congenital Adrenal Hyperplasia

Congenital adrenal hyperplasia (CAH) refers to a group of autosomal recessive disorders due to single-gene defects in the various enzymes required for cortisol biosynthesis. CAH represents a continuous phenotypic spectrum with more than 95% of all cases caused by 21-hydroxylase deficiency. Genotyping is an important tool in confirming the diagnosis or carrier state, provides prognostic information on disease severity, and is essential for genetic counseling. In this article, the authors provide an in-depth discussion on the genetics of CAH, including genetic diagnosis, molecular analysis, genotype-phenotype relationships, and counseling of patients and their families.

The neurogenetics of alternative splicing

Alternative precursor-mRNA splicing is a key mechanism for regulating gene expression in mammals and is controlled by specialized RNA-binding proteins. The misregulation of splicing is implicated in multiple neurological disorders. We describe recent mouse genetic studies of alternative splicing that reveal its critical role in both neuronal development and the function of mature neurons. We discuss the challenges in understanding the extensive genetic programmes controlled by proteins that regulate splicing, both during development and in the adult brain.

Bringing genetic background into focus

Researchers should embrace differences in genetic background to build richer disease models that more accurately reflect the level of variation in the human population, posits Clement Chow.

An update on transcriptional and post-translational regulation of brain voltage-gated sodium channels

Voltage-gated sodium channels are essential proteins in brain physiology, as they generate the sodium currents that initiate neuronal action potentials. Voltage-gated sodium channels expression, localisation and function are regulated by a range of transcriptional and post-translational mechanisms. Here, we review our understanding of regulation of brain voltage-gated sodium channels, in particular SCN1A (Na_V1.1), SCN2A (Na_V1.2), SCN3A (Na_V1.3) and SCN8A (Na_V1.6), by transcription factors, by alternative splicing, and by post-translational modifications. Our focus is strongly centred on recent research lines, and newly generated knowledge.

Genetic dissection of horizontal cell inhibitory signaling in mice in complete darkness in vivo

PURPOSE

To test the hypothesis that horizontal cell (HC) inhibitory signaling controls the degree to which rod cell membranes are depolarized as measured by the extent to which L-type calcium channels (LTCCs) are open in complete darkness in the mouse retina in vivo.

METHODS

Dark-adapted wild-type (wt), CACNA1F (Ca(v)1.4(-/-)), arrestin-1 (Arr1(-/-)), and CACNA1D (Ca(v)1.3(-/-)) C57Bl/6 mice were studied. Manganese-enhanced MRI (MEMRI) evaluated the extent that rod LTCCs are open as an index of loss of HC inhibitory signaling. Subgroups were pretreated with D-cis-diltiazem (DIL) at a dose that specifically antagonizes Ca(v)1.2 channels in vivo.

RESULTS

Knockout mice predicted to have impaired HC inhibitory signaling (Ca(v)1.4(-/-) or Arr1(-/-)) exhibited greater than normal rod manganese uptake; inner retinal uptake was also supernormal. Genetically knocking out a closely associated gene not expected to impact HC inhibitory signaling (CACNA1D) did not generate this phenotype. The Arr1(-/-) mice exhibited the largest rod uptake of manganese. Manganese-enhanced MRI of DIL-treated Arr1(-/-) mice suggested a greater number of operant LTCC subtypes (i.e., Ca(v)1.2, 1.3, and 1.4) in rods and inner retina than that in DIL-treated Ca(v)1.4(-/-) mice (i.e., Ca(v)1.3). The Ca(v)1.3(-/-) + DIL-treated mice exhibited evidence for a compensatory contribution from Ca(v)1.2 LTCCs.

CONCLUSIONS

The data suggest that loss of HC inhibitory signaling is the proximate cause leading to maximally open LTCCs in rods, and possibly inner retinal cells, in mice in total darkness in vivo, regardless of compensatory changes in LTCC subtype manifested in the mutant mice.

Reduced Nav1.6 Sodium Channel Activity in Mice Increases In Vivo Sensitivity to Volatile Anesthetics

Na_v1.6 is a major voltage-gated sodium channel in the central and peripheral nervous systems. Within neurons, the channel protein is concentrated at the axon initial segment and nodes of Ranvier, where it functions in initiation and propagation of action potentials. We examined the role of Na_v1.6 in general anesthesia using two mouse mutants with reduced activity of Na_v1.6, Scn8a^medJ/medJ and Scn8a^9J/9J. The mice were exposed to the general anesthetics isoflurane and sevoflurane in step-wise increments; the concentration required to produce loss of righting reflex, a surrogate for anesthetic-induced unconsciousness in rodents, was determined. Mice homozygous for these mutations exhibited increased sensitivity to both isoflurane and sevoflurane. The increased sensitivity was observed during induction of unconsciousness but not during the recovery phase, suggesting that the effect is not attributable to compromised systemic physiology. Electroencephalographic theta power during baseline waking was lower in mutants, suggesting decreased arousal and reduced neuronal excitability. This is the first report linking reduced activity of a specific voltage-gated sodium channel to increased sensitivity to general anesthetics in vivo.

Contrasting genetic architectures in different mouse reference populations used for studying complex traits

Quantitative trait loci (QTLs) are being used to study genetic networks, protein functions, and systems properties that underlie phenotypic variation and disease risk in humans, model organisms, agricultural species, and natural populations. The challenges are many, beginning with the seemingly simple tasks of mapping QTLs and identifying their underlying genetic determinants. Various specialized resources have been developed to study complex traits in many model organisms. In the mouse, remarkably different pictures of genetic architectures are emerging. Chromosome Substitution Strains (CSSs) reveal many QTLs, large phenotypic effects, pervasive epistasis, and readily identified genetic variants. In contrast, other resources as well as genome-wide association studies (GWAS) in humans and other species reveal genetic architectures dominated with a relatively modest number of QTLs that have small individual and combined phenotypic effects. These contrasting architectures are the result of intrinsic differences in the study designs underlying different resources. The CSSs examine context-dependent phenotypic effects independently among individual genotypes, whereas with GWAS and other mouse resources, the average effect of each QTL is assessed among many individuals with heterogeneous genetic backgrounds. We argue that variation of genetic architectures among individuals is as important as population averages. Each of these important resources has particular merits and specific applications for these individual and population perspectives. Collectively, these resources together with high-throughput genotyping, sequencing and genetic engineering technologies, and information repositories highlight the power of the mouse for genetic, functional, and systems studies of complex traits and disease models.

Glycine receptor mouse mutants: model systems for human hyperekplexia

Human hyperekplexia is a neuromotor disorder caused by disturbances in inhibitory glycine-mediated neurotransmission. Mutations in genes encoding for glycine receptor subunits or associated proteins, such as GLRA1, GLRB, GPHN and ARHGEF9, have been detected in patients suffering from hyperekplexia. Classical symptoms are exaggerated startle attacks upon unexpected acoustic or tactile stimuli, massive tremor, loss of postural control during startle and apnoea. Usually patients are treated with clonazepam, this helps to dampen the severe symptoms most probably by up-regulating GABAergic responses. However, the mechanism is not completely understood. Similar neuromotor phenotypes have been observed in mouse models that carry glycine receptor mutations. These mouse models serve as excellent tools for analysing the underlying pathomechanisms. Yet, studies in mutant mice looking for postsynaptic compensation of glycinergic dysfunction via an up-regulation in GABAA receptor numbers have failed, as expression levels were similar to those in wild-type mice. However, presynaptic adaptation mechanisms with an unusual switch from mixed GABA/glycinergic to GABAergic presynaptic terminals have been observed. Whether this presynaptic adaptation explains the improvement in symptoms or other compensation mechanisms exist is still under investigation. With the help of spontaneous glycine receptor mouse mutants, knock-in and knock-out studies, it is possible to associate behavioural changes with pharmacological differences in glycinergic inhibition. This review focuses on the structural and functional characteristics of the various mouse models used to elucidate the underlying signal transduction pathways and adaptation processes and describes a novel route that uses gene-therapeutic modulation of mutated receptors to overcome loss of function mutations.

A mutation in a splicing factor that causes retinitis pigmentosa has a transcriptome-wide effect on mRNA splicing

BACKGROUND

Substantial progress has been made in the identification of sequence elements that control mRNA splicing and the genetic variants in these elements that alter mRNA splicing (referred to as splicing quantitative trait loci - sQTLs). Genetic variants that affect mRNA splicing in trans are harder to identify because their effects can be more subtle and diffuse, and the variants are not co-located with their targets. We carried out a transcriptome-wide analysis of the effects of a mutation in a ubiquitous splicing factor that causes retinitis pigmentosa (RP) on mRNA splicing, using exon microarrays.

RESULTS

Exon microarray data was generated from whole blood samples obtained from four individuals with a mutation in the splicing factor PRPF8 and four sibling controls. Although the mutation has no known phenotype in blood, there was evidence of widespread differences in splicing between cases and controls (affecting approximately 20% of exons). Most probesets with significantly different inclusion (defined as the expression intensity of the exon divided by the expression of the corresponding transcript) between cases and controls had higher inclusion in cases and corresponded to exons that were shorter than average, AT rich, located towards the 5' end of the gene and flanked by long introns. Introns flanking affected probesets were particularly depleted for the shortest category of introns, associated with splicing via intron definition.

CONCLUSIONS

Our results show that a mutation in a splicing factor, with a phenotype that is restricted to retinal tissue, acts as a trans-sQTL cluster in whole blood samples. Characteristics of the affected exons suggest that they are spliced co-transcriptionally and via exon definition. However, due to the small sample size available for this study, further studies are required to confirm the widespread impact of this PRPF8 mutation on mRNA splicing outside the retina.

Cerebellar ataxia and functional genomics: Identifying the routes to cerebellar neurodegeneration

Cerebellar ataxias are progressive neurodegenerative disorders characterized by atrophy of the cerebellum leading to motor dysfunction, balance problems, and limb and gait ataxia. These include among others, the dominantly inherited spinocerebellar ataxias, recessive cerebellar ataxias such as Friedreich's ataxia, and X-linked cerebellar ataxias. Since all cerebellar ataxias display considerable overlap in their disease phenotypes, common pathological pathways must underlie the selective cerebellar neurodegeneration. Therefore, it is important to identify the molecular mechanisms and routes to neurodegeneration that cause cerebellar ataxia. In this review, we discuss the use of functional genomic approaches including whole-exome sequencing, genome-wide gene expression profiling, miRNA profiling, epigenetic profiling, and genetic modifier screens to reveal the underlying pathogenesis of various cerebellar ataxias. These approaches have resulted in the identification of many disease genes, modifier genes, and biomarkers correlating with specific stages of the disease. This article is part of a Special Issue entitled: From Genome to Function.

Mouse genetic background influences the dental phenotype

Dental enamel covers the crown of the vertebrate tooth and is considered to be the hardest tissue in the body. Enamel develops during secretion of an extracellular matrix by ameloblast cells in the tooth germ, prior to eruption of the tooth into the oral cavity. Secreted enamel proteins direct mineralization patterns during the maturation stage of amelogenesis as the tooth prepares to erupt. The amelogenins are the most abundant enamel proteins and are required for normal enamel development. Phenotypic differences were observed between incisors from individual Amelx (amelogenin) null mice that had a mixed 129xC57BL/6J genetic background and between inbred wild-type (WT) mice with different genetic backgrounds (C57BL/6J, C3H/HeJ, FVB/NJ). We hypothesized that this could be due to modifier genes, as human patients with a mutation in an enamel protein gene causing the enamel defect amelogenesis imperfecta (AI) can also have varied appearance of dentitions within a kindred. Enamel density measurements varied for all WT inbred strains midway during incisor development. Enamel thickness varied between some WT strains, and, unexpectedly, dentin density varied extensively between incisors and molars of all WT and Amelx null strains studied. WTFVB/NJ incisors were more similar to those of Amelx null mice than to those of the other WT strains in terms of incisor height/width ratio and pattern of enamel mineralization. Strain-specific differences led to the conclusion that modifier genes may be implicated in determining both normal development and severity of enamel appearance in AI mouse models and may in future studies be related to phenotypic heterogeneity within human AI kindreds reported in the literature.

Mouse models of PI(3,5)P2 deficiency with impaired lysosome function

The endolysosomal system and autophagy are essential components of macromolecular turnover in eukaryotic cells. The low-abundance signaling lipid PI(3,5)P2 is a key regulator of this pathway. Analysis of mouse models with defects in PI(3,5)P2 biosynthesis has revealed the unique dependence of the mammalian nervous system on this signaling pathway. This insight led to the discovery of the molecular basis for several human neurological disorders, including Charcot-Marie-Tooth disease and Yunis-Varon syndrome. Spontaneous mutants, conditional knockouts, transgenic lines, and gene-trap alleles of Fig4, Vac14, and Pikfyve (Fab1) in the mouse have provided novel information regarding the role of PI(3,5)P2in vivo. This review summarizes what has been learned from mouse models and highlights the utility of manipulating complex signaling pathways in vivo.

Sodium channel SCN8A (Nav1.6): properties and de novo mutations in epileptic encephalopathy and intellectual disability

The sodium channel Nav1.6, encoded by the gene SCN8A, is one of the major voltage-gated channels in human brain. The sequences of sodium channels have been highly conserved during evolution, and minor changes in biophysical properties can have a major impact in vivo. Insight into the role of Nav1.6 has come from analysis of spontaneous and induced mutations of mouse Scn8a during the past 18 years. Only within the past year has the role of SCN8A in human disease become apparent from whole exome and genome sequences of patients with sporadic disease. Unique features of Nav1.6 include its contribution to persistent current, resurgent current, repetitive neuronal firing, and subcellular localization at the axon initial segment (AIS) and nodes of Ranvier. Loss of Nav1.6 activity results in reduced neuronal excitability, while gain-of-function mutations can increase neuronal excitability. Mouse Scn8a (med) mutants exhibit movement disorders including ataxia, tremor and dystonia. Thus far, more than ten human de novo mutations have been identified in patients with two types of disorders, epileptic encephalopathy and intellectual disability. We review these human mutations as well as the unique features of Nav1.6 that contribute to its role in determining neuronal excitability in vivo. A supplemental figure illustrating the positions of amino acid residues within the four domains and 24 transmembrane segments of Nav1.6 is provided to facilitate the location of novel mutations within the channel protein.

Aberrant splicing in neurological diseases

Splicing of precursor messenger RNA (pre-mRNA) removes the intervening sequences (introns) and joins the expressed regions (exons) in the nucleus, before an intron-containing eukaryotic mRNA transcript can be exported and translated into proteins in the cytoplasm. While some sequences are always included or excluded (constitutive splicing), others can be selectively used (alternative splicing) in this process. Particularly by alternative splicing, up to tens of thousands of variant transcripts can be produced from a single gene, which contributes greatly to the proteomic diversity for such complex cellular functions as 'wiring' neurons in the nervous system. Disruption of this process leads to aberrant splicing, which accounts for the defects of up to 50% of mutations that cause certain human genetic diseases. In this review, we describe the different mechanisms of aberrant splicing that cause or have been associated with neurological diseases.

Rapid screening of yeast mutants with reporters identifies new splicing phenotypes

Nuclear precursor mRNA splicing requires the stepwise assembly of a large complex, the spliceosome. Recent large-scale analyses, including purification of splicing complexes, high-throughput genetic screens and interactomic studies, have linked numerous factors to this dynamic process, including a well-defined core conserved from yeast to human. Intriguingly, despite extensive studies, no splicing defects were reported for some of the corresponding yeast mutants. To resolve this paradox, we screened a collection of viable yeast strains carrying mutations in splicing-related factors with a set of reporters including artificial constructs carrying competing splice sites. Previous analyses have indeed demonstrated that this strategy identifies yeast factors able to regulate alternative splicing and whose properties are conserved in human cells. The method, sensitive to subtle defects, revealed new splicing phenotypes for most analyzed factors such as the Urn1 protein. Interestingly, a mutant of PRP8 specifically lacking an N-terminal proline-rich region stimulated the splicing of a reporter containing competing branchpoint/3' splice site regions. Thus, using appropriate reporters, yeast can be used to quickly delineate the effect of various factors on splicing and identify those with the propensity to regulate alternative splicing events.

Modifier genes and the plasticity of genetic networks in mice

Modifier genes are an integral part of the genetic landscape in both humans and experimental organisms, but have been less well explored in mammals than other systems. A growing number of modifier genes in mouse models of disease nonetheless illustrate the potential for novel findings, while new technical advances promise many more to come. Modifier genes in mouse models include induced mutations and spontaneous or wild-derived variations captured in inbred strains. Identification of modifiers among wild-derived variants in particular should detect disease modifiers that have been shaped by selection and might therefore be compatible with high fitness and function. Here we review selected examples and argue that modifier genes derived from natural variation may provide a bias for nodes in genetic networks that have greater intrinsic plasticity and whose therapeutic manipulation may therefore be more resilient to side effects than conventional targets.

MATS: a Bayesian framework for flexible detection of differential alternative splicing from RNA-Seq data

Ultra-deep RNA sequencing has become a powerful approach for genome-wide analysis of pre-mRNA alternative splicing. We develop MATS (multivariate analysis of transcript splicing), a Bayesian statistical framework for flexible hypothesis testing of differential alternative splicing patterns on RNA-Seq data. MATS uses a multivariate uniform prior to model the between-sample correlation in exon splicing patterns, and a Markov chain Monte Carlo (MCMC) method coupled with a simulation-based adaptive sampling procedure to calculate the P-value and false discovery rate (FDR) of differential alternative splicing. Importantly, the MATS approach is applicable to almost any type of null hypotheses of interest, providing the flexibility to identify differential alternative splicing events that match a given user-defined pattern. We evaluated the performance of MATS using simulated and real RNA-Seq data sets. In the RNA-Seq analysis of alternative splicing events regulated by the epithelial-specific splicing factor ESRP1, we obtained a high RT–PCR validation rate of 86% for differential exon skipping events with a MATS FDR of <10%. Additionally, over the full list of RT–PCR tested exons, the MATS FDR estimates matched well with the experimental validation rate. Our results demonstrate that MATS is an effective and flexible approach for detecting differential alternative splicing from RNA-Seq data.

Rbfox proteins regulate alternative splicing of neuronal sodium channel SCN8A

The SCN8A gene encodes the voltage-gated sodium channel Na(v)1.6, a major channel in neurons of the CNS and PNS. SCN8A contains two alternative exons,18N and 18A, that exhibit tissue specific splicing. In brain, the major SCN8A transcript contains exon 18A and encodes the full-length sodium channel. In other tissues, the major transcript contains exon 18N and encodes a truncated protein, due to the presence of an in-frame stop codon. Selection of exon 18A is therefore essential for generation of a functional channel protein, but the proteins involved in this selection have not been identified. Using a 2.6 kb Scn8a minigene containing exons 18N and 18A, we demonstrate that co-transfection with Fox-1 or Fox-2 initiates inclusion of exon 18A. This effect is dependent on the consensus Fox binding site located 28 bp downstream of exon 18A. We examined the alternative splicing of human SCN8A and found that the postnatal switch to exon 18A is completed later than 10 months of age. In purified cell populations, transcripts containing exon 18A predominate in neurons but are not present in oligodendrocytes or astrocytes. Transcripts containing exon 18N appear to be degraded by nonsense-mediated decay in HEK cells. Our data indicate that RBFOX proteins contribute to the cell-specific expression of Na(v)1.6 channels in mature neurons.

A modifier gene alleviates hypothyroidism-induced hearing impairment in Pou1f1dw dwarf mice

Thyroid hormone has pleiotropic effects on cochlear development, and genomic variation influences the severity of associated hearing deficits. DW/J-Pou1f1dw/dw mutant mice lack pituitary thyrotropin, which causes severe thyroid hormone deficiency and profound hearing impairment. To assess the genetic complexity of protective effects on hypothyroidism-induced hearing impairment, an F1 intercross was generated between DW/J-Pou1f1dw/+ carriers and an inbred strain with excellent hearing derived from Mus castaneus, CAST/EiJ. Approximately 24% of the (DW/J×CAST/EiJ) Pou1f1dw/dw F2 progeny had normal hearing. A genome scan revealed a locus on chromosome 2, named modifier of dw hearing, or Mdwh, that rescues hearing despite persistent hypothyroidism. This chromosomal region contains the modifier of tubby hearing 1 (Moth1) locus that encodes a protective allele of the microtubule-associated protein MTAP1A. DW/J-Pou1f1dw/+ carriers were crossed with the AKR strain, which also carries a protective allele of Mtap1a, and we found that AKR is not protective for hearing in the (DW/J×AKR) Pou1f1dw/dw F2 progeny. Thus, protective alleles of Mtap1a are not sufficient to rescue DW/J-Pou1f1dw/dw hearing. We expect that identification of protective modifiers will enhance our understanding of the mechanisms of hypothyroidism-induced hearing impairment.

Modifier genes and non-genetic factors reshape anatomical deficits in Zfp423-deficient mice

Development of neural circuitry depends on the integration of signaling pathways to coordinate specification, proliferation and differentiation of cell types in the right number, in the right place, at the right time. Zinc finger protein 423 (Zfp423), a 30-zinc finger transcription factor, forms alternate complexes with components of several developmental signaling pathways, suggesting it as a point of signal integration during brain development. We previously showed that mice lacking Zfp423 have reduced proliferation of cerebellar precursor cells, resulting in complete loss of vermis and variable hypoplasia of cerebellar hemispheres. Here, we show that Zfp423(-/-) hemisphere malformations are shaped by both genetic and non-genetic factors, producing distinct phenotype distributions in different inbred genetic backgrounds. In genetic mapping studies, we identify four additive modifier loci (Amzn1-4) and seven synthetically interacting loci (Smzn1.1-3.1) that together explain approximately one-third of the phenotypic variance. Strain-specific sequence polymorphism and expression data provide a reduced list of functional variant candidate genes at each modifier locus. Environmental covariates add only modest explanatory power, suggesting an additional stochastic component. These results provide a comprehensive analysis of sources of phenotype variation in a model of hindbrain malformation.

Pathogenic mechanism of the FIG4 mutation responsible for Charcot-Marie-Tooth disease CMT4J

CMT4J is a severe form of Charcot-Marie-Tooth neuropathy caused by mutation of the phosphoinositide phosphatase FIG4/SAC3. Affected individuals are compound heterozygotes carrying the missense allele FIG4-I41T in combination with a null allele. Analysis using the yeast two-hybrid system demonstrated that the I41T mutation impairs interaction of FIG4 with the scaffold protein VAC14. The critical role of this interaction was confirmed by the demonstration of loss of FIG4 protein in VAC14 null mice. We developed a mouse model of CMT4J by expressing a Fig4-I41T cDNA transgene on the Fig4 null background. Expression of the mutant transcript at a level 5× higher than endogenous Fig4 completely rescued lethality, whereas 2× expression gave only partial rescue, providing a model of the human disease. The level of FIG4-I41T protein in transgenic tissues is only 2% of that predicted by the transcript level, as a consequence of the protein instability caused by impaired interaction of the mutant protein with VAC14. Analysis of patient fibroblasts demonstrated a comparably low level of mutant I41T protein. The abundance of FIG4-I41T protein in cultured cells is increased by treatment with the proteasome inhibitor MG-132. The data demonstrate that FIG4-I41T is a hypomorphic allele encoding a protein that is unstable in vivo. Expression of FIG4-I41T protein at 10% of normal level is sufficient for long-term survival, suggesting that patients with CMT4J could be treated by increased production or stabilization of the mutant protein. The transgenic model will be useful for testing in vivo interventions to increase the abundance of the mutant protein.

Charcot-Marie-Tooth disease type 4J is a severe neurological disorder with childhood or adult onset and progression to loss of mobility and death. Patients inherit a mutation that changes amino acid residue 41 of the FIG4 protein from isoleucine to threonine. We report that this mutation destabilizes the FIG4 protein by blocking its interaction with a stabilizing protein partner. We developed a mouse model of CMT4J and found that a low level of expression of the mutant protein, 10% of wildtype level, is sufficient to prevent lethality. This work provides the scientific basis for development of a directed treatment for this rare, lethal disorder.

Genetic modifiers of neurological disease

Genetic modifiers make an important contribution to neurological disease phenotypes. Significant progress has been made by studying genetic modifiers in model organisms. The ability to study complex genetic interactions in model systems contributes to our understanding of the genetic factors that influence neurological disease. This will lead to the development of novel therapeutic strategies and personalized treatment based on genetic risk.

Developments in RNA splicing and disease

Pre-mRNA processing, including 5'-end capping, splicing, editing, and polyadenylation, consists of a series of orchestrated and primarily cotranscriptional steps that ensure both the high fidelity and extreme diversity characteristic of eukaryotic gene expression. Alternative splicing and editing allow relatively small genomes to encode vast proteomic arrays while alternative 3'-end formation enables variations in mRNA localization, translation, and stability. Of course, this mechanistic complexity comes at a high price. Mutations in the myriad of RNA sequence elements that regulate mRNA biogenesis, as well as the trans-acting factors that act upon these sequences, underlie a number of human diseases. In this review, we focus on one of these key RNA processing steps, splicing, to highlight recent studies that describe both conventional and novel pathogenic mechanisms that underlie muscle and neurological diseases.

Epileptogenic ion channel mutations: from bedside to bench and, hopefully, back again

Mutations of genes coding for ion channels cause several genetically determined human epileptic syndromes. The identification of a gene variant linked to a particular disease gives important information, but it is usually necessary to perform functional studies in order to completely disclose the pathogenic mechanisms. The functional consequences of epileptogenic mutations have been studied both in vitro and in vivo with several experimental systems, studies that have provided significant knowledge on the pathogenic mechanisms that leads to inherited human epilepsies, and possibly also on the pathogenic mechanisms of non-genetic human epilepsies due to "acquired channelopathies". However, several open issues remain and difficulties in the interpretation of the experimental data have arisen that limit translational applications. We will highlight the value and the limits of different approaches to the study of epileptogenic channelopathies, focussing on the importance of the experimental systems in the assessment of the functional effects of the mutations and on the possible applications of the obtained results to the clinical practice.

Evolution of alternative splicing in primate brain transcriptomes

Alternative splicing is a predominant form of gene regulation in higher eukaryotes. The evolution of alternative splicing provides an important mechanism for the acquisition of novel gene functions. In this work, we carried out a genome-wide phylogenetic survey of lineage-specific splicing patterns in the primate brain, via high-density exon junction array profiling of brain transcriptomes of humans, chimpanzees and rhesus macaques. We identified 509 genes showing splicing differences among these species. RT-PCR analysis of 40 exons confirmed the predicted splicing evolution of 33 exons. Of these 33 exons, outgroup analysis using rhesus macaques confirmed 13 exons with human-specific increase or decrease in transcript inclusion levels after humans diverged from chimpanzees. Some of the human-specific brain splicing patterns disrupt domains critical for protein-protein interactions, and some modulate translational efficiency of their host genes. Strikingly, for exons showing splicing differences across species, we observed a significant increase in the rate of silent substitutions within exons, coupled with accelerated sequence divergence in flanking introns. This indicates that evolution of cis-regulatory signals is a major contributor to the emergence of human-specific splicing patterns. In one gene (MAGOH), using minigene reporter assays, we demonstrated that the combination of two human-specific cis-sequence changes created its human-specific splicing pattern. Together, our data reveal widespread human-specific changes of alternative splicing in the brain and suggest an important role of splicing in the evolution of neuronal gene regulation and functions.

Modifiers of von Willebrand factor identified by natural variation in inbred strains of mice

Type 1 von Willebrand disease (VWD) is the most common inherited human bleeding disorder. However, diagnosis is complicated by incomplete penetrance and variable expressivity, as well as wide variation in von Willebrand factor (VWF) levels among the normal population. Previous work has exploited the highly variable plasma VWF levels among inbred strains of mice to identify 2 major regulators, Mvwf1 and Mvwf2 (modifier of VWF). Mvwf1 is a glycosyltransferase and Mvwf2 is a natural variant in Vwf that alters biosynthesis. We report the identification of an additional alteration at the Vwf locus (Mvwf5), as well as 2 loci unlinked to Vwf (Mvwf6-7) using a backcross approach with the inbred mouse strains WSB/EiJ and C57BL/6J. Through positional cloning, we show that Mvwf5 is a cis-regulatory variant that alters Vwf mRNA expression. A similar mechanism could potentially explain a significant percentage of human VWD cases, especially those with no detectable mutation in the VWF coding sequence. Mvwf6 displays conservation of synteny with potential VWF modifier loci identified in human pedigrees, suggesting that its ortholog may modify VWF in human populations.

Fine mapping of an epilepsy modifier gene on mouse Chromosome 19

Mutations in voltage-gated sodium channels are associated with several types of human epilepsy. Variable expressivity and penetrance are common features of inherited epilepsy caused by sodium channel mutations, suggesting that genetic modifiers may influence clinical severity. The mouse model Scn2a(Q54) has an epilepsy phenotype due to a mutation in Scn2a that results in elevated persistent sodium current. Phenotype severity in Scn2a(Q54) mice is dependent on the genetic background. Congenic C57BL/6J.Q54 mice have delayed onset and low seizure frequency compared to (C57BL/6J x SJL/J)F1.Q54 mice. Previously, we identified two modifier loci that influence the Scn2a(Q54) epilepsy phenotype: Moe1 (modifier of epilepsy 1) on Chromosome 11 and Moe2 on Chromosome 19. We have constructed interval-specific congenic strains to further refine the position of Moe2 on Chromosome 19 to a 5-Mb region. Sequencing and expression analyses of genes in the critical interval suggested two potential modifier candidates: (1) voltage-gated potassium channel subunit subfamily V, member 2 (Kcnv2), and (2) SWI/SNF-related, matrix-associated, actin-dependent regulator of chromatin, subfamily a, member 2 (Smarca2). Based on its biological role in regulating membrane excitability and the association between ion channel variants and seizures, Kcnv2 is a strong functional candidate for Moe2. Modifier genes affecting the epilepsy phenotype of Scn2a(Q54) mice may contribute to variable expressivity and penetrance in human epilepsy patients with sodium channel mutations.

The ataxia3 mutation in the N-terminal cytoplasmic domain of sodium channel Na(v)1.6 disrupts intracellular trafficking

The ENU-induced neurological mutant ataxia3 was mapped to distal mouse chromosome 15. Sequencing of the positional candidate gene Scn8a encoding the sodium channel Na(v)1.6 identified a T>C transition in exon 1 resulting in the amino acid substitution p.S21P near the N terminus of the channel. The cytoplasmic N-terminal region is evolutionarily conserved but its function has not been well characterized. ataxia3 homozygotes exhibit a severe disorder that includes ataxia, tremor, and juvenile lethality. Unlike Scn8a null mice, they retain partial hindlimb function. The mutant transcript is stable but protein abundance is reduced and the mutant channel is not detected in its usual site of concentration at nodes of Ranvier. In whole-cell patch-clamp studies of transfected ND7/23 cells that were maintained at 37 degrees C, the mutant channel did not produce sodium current, and function was not restored by coexpression of beta1 and beta2 subunits. However, when transfected cells were maintained at 30 degrees C, the mutant channel generated voltage-dependent inward sodium currents with an average peak current density comparable with wild type, demonstrating recovery of channel activity. Immunohistochemistry of primary cerebellar granule cells from ataxia3 mice demonstrated that the mutant protein is retained in the cis-Golgi. This trafficking defect can account for the low level of Na(v)1.6-S21P at nodes of Ranvier in vivo and at the surface of transfected cells. The data demonstrate that the cytoplasmic N-terminal domain of the sodium channel is required for anterograde transport from the Golgi complex to the plasma membrane.

Genetics of congenital adrenal hyperplasia

Congenital adrenal hyperplasia (CAH) is one of the most common inherited metabolic disorders. It comprises a group of autosomal recessive disorders caused by the deficiency of one of four steroidogenic enzymes involved in cortisol biosynthesis or in the electron donor enzyme P450 oxidoreductase (POR) that serves as electron donor to steroidogenic cytochrome P450 (CYP) type II enzymes. The biochemical and clinical phenotype depends on the specific enzymatic defect and the impairment of specific enzyme activity. Defects of steroid 21-hydroxylase (CYP21A2) and 11beta-hydroxylase (CYP11B1) only affect adrenal steroidogenesis, whereas 17alpha-hydroxylase (CYP17A1) and 3beta-hydroxysteroid dehydrogenase type 2 (HSD3B2) deficiency also impact on gonadal steroid biosynthesis. Inactivating POR gene mutations are the cause of CAH manifesting with apparent combined CYP17A1-CYP21A2 deficiency. P450 oxidoreductase deficiency (ORD) has a complex phenotype including two unique features not observed in any other CAH variant: skeletal malformations and severe genital ambiguity in both sexes.

A targeted deleterious allele of the splicing factor SCNM1 in the mouse

The auxiliary spliceosomal protein SCNM1 contributes to recognition of nonconsensus splice donor sites. SCNM1 was first identified as a modifier of the severity of a sodium channelopathy in the mouse. The most severely affected strain, C57BL/6J, carries the variant allele SCNM1R187X, which is defective in splicing the mutated donor site in the Scn8a(medJ) transcript. To further probe the in vivo function of SCNM1, we constructed a floxed allele and generated a mouse with constitutive deletion of exons 3-5. The SCNM1Delta3-5 protein is produced and correctly localized to the nucleus, but is more functionally impaired than the C57BL/6J allele. Deficiency of SCNM1 did not significantly alter other brain transcripts. We characterized an ENU-induced allele of Scnm1 and evaluated the ability of wild-type SCNM1 to rescue lethal mutations of I-mfa and Brunol4. The phenotypes of the Scnm1Delta3-5 mutant confirm the role of this splice factor in processing the Scn8a(medJ) transcript and provide a new allele of greater severity for future studies.

Genome-wide survey of allele-specific splicing in humans

Background

Accurate mRNA splicing depends on multiple regulatory signals encoded in the transcribed RNA sequence. Many examples of mutations within human splice regulatory regions that alter splicing qualitatively or quantitatively have been reported and allelic differences in mRNA splicing are likely to be a common and important source of phenotypic diversity at the molecular level, in addition to their contribution to genetic disease susceptibility. However, because the effect of a mutation on the efficiency of mRNA splicing is often difficult to predict, many mutations that cause disease through an effect on splicing are likely to remain undiscovered.

Results

We have combined a genome-wide scan for sequence polymorphisms likely to affect mRNA splicing with analysis of publicly available Expressed Sequence Tag (EST) and exon array data. The genome-wide scan uses published tools and identified 30,977 SNPs located within donor and acceptor splice sites, branch points and exonic splicing enhancer elements. For 1,185 candidate splicing polymorphisms the difference in splicing between alternative alleles was corroborated by publicly available exon array data from 166 lymphoblastoid cell lines. We developed a novel probabilistic method to infer allele-specific splicing from EST data. The method uses SNPs and alternative mRNA isoforms mapped to EST sequences and models both regulated alternative splicing as well as allele-specific splicing. We have also estimated heritability of splicing and report that a greater proportion of genes show evidence of splicing heritability than show heritability of overall gene expression level. Our results provide an extensive resource that can be used to assess the possible effect on splicing of human polymorphisms in putative splice-regulatory sites.

Conclusion

We report a set of genes showing evidence of allele-specific splicing from an integrated analysis of genomic polymorphisms, EST data and exon array data, including several examples for which there is experimental evidence of polymorphisms affecting splicing in the literature. We also present a set of novel allele-specific splicing candidates and discuss the strengths and weaknesses of alternative technologies for inferring the effect of sequence variants on mRNA splicing.

Mouse forward genetics in the study of the peripheral nervous system and human peripheral neuropathy

Forward genetics, the phenotype-driven approach to investigating gene identity and function, has a long history in mouse genetics. Random mutations in the mouse transcend bias about gene function and provide avenues towards unique discoveries. The study of the peripheral nervous system is no exception; from historical strains such as the trembler mouse, which led to the identification of PMP22 as a human disease gene causing multiple forms of peripheral neuropathy, to the more recent identification of the claw paw and sprawling mutations, forward genetics has long been a tool for probing the physiology, pathogenesis, and genetics of the PNS. Even as spontaneous and mutagenized mice continue to enable the identification of novel genes, provide allelic series for detailed functional studies, and generate models useful for clinical research, new methods, such as the piggyBac transposon, are being developed to further harness the power of forward genetics.