Exon inclusion is dependent on predictable exonic splicing enhancers

Data-driven insights to inform splice-altering variant assessment

Disease-causing genetic variants often disrupt mRNA splicing, an intricate process that is incompletely understood. Thus, accurate inference of which genetic variants will affect splicing and what their functional consequences will be is challenging, particularly for variants outside of the essential splice sites. Here, we describe a set of data-driven heuristics that inform the interpretation of human splice-altering variants (SAVs) based on the analysis of annotated exons, experimentally validated SAVs, and the currently understood principles of splicing biology. We defined requisite splicing criteria by examining around 202,000 canonical protein-coding exons and 19,000 experimentally validated splicing branchpoints. This analysis defined the sequence, spacing, and motif strength required for splicing, with 95.9% of the exons examined meeting these criteria. By considering over 12,000 experimentally validated variants from the SpliceVarDB, we defined a set of heuristics that inform the evaluation of putative SAVs. To ensure the applicability of each heuristic, only those supported by at least 10 experimentally validated variants were considered. This allowed us to establish a measure of spliceogenicity: the proportion of variants at a location (or motif site) that affected splicing in a given context. This study makes considerable advances toward bridging the gap between computational predictions and the biological process of splicing, offering an evidence-based approach to identifying SAVs and evaluating their impact. Our splicing heuristics enhance the current framework for genetic variant evaluation with a robust, detailed, and comprehensible analysis by adding valuable context over traditional binary prediction tools.

All exons are not created equal-exon vulnerability determines the effect of exonic mutations on splicing

It is now widely accepted that aberrant splicing of constitutive exons is often caused by mutations affecting cis-acting splicing regulatory elements (SREs), but there is a misconception that all exons have an equal dependency on SREs and thus a similar vulnerability to aberrant splicing. We demonstrate that some exons are more likely to be affected by exonic splicing mutations (ESMs) due to an inherent vulnerability, which is context dependent and influenced by the strength of exon definition. We have developed VulExMap, a tool which is based on empirical data that can designate whether a constitutive exon is vulnerable. Using VulExMap, we find that only 25% of all exons can be categorized as vulnerable, whereas two-thirds of 359 previously reported ESMs in 75 disease genes are located in vulnerable exons. Because VulExMap analysis is based on empirical data on splicing of exons in their endogenous context, it includes all features important in determining the vulnerability. We believe that VulExMap will be an important tool when assessing the effect of exonic mutations by pinpointing whether they are located in exons vulnerable to ESMs.

The human genome contains over a million autonomous exons

Mammalian mRNA and lncRNA exons are often small compared to introns. The exon definition model predicts that exons splice autonomously, dependent on proximal exon sequence features, explaining their delineation within large introns. This model has not been examined on a genome-wide scale, however, leaving open the question of how often mRNA and lncRNA exons are autonomous. It is also unknown how frequently such exons can arise by chance. Here, we directly assayed large fragments (500-1000 bp) of the human genome by exon trapping, which detects exons spliced into a heterologous transgene, here designed with a large intron context. We define the trapped exons as "autonomous." We obtained ∼1.25 million trapped exons, including most known mRNA and well-annotated lncRNA internal exons, demonstrating that human exons are predominantly autonomous. mRNA exons are trapped with the highest efficiency. Nearly a million of the trapped exons are unannotated, most located in intergenic regions and antisense to mRNA, with depletion from the forward strand of introns. These exons are not conserved, suggesting they are nonfunctional and arose from random mutations. They are nonetheless highly enriched with known splicing promoting sequence features that delineate known exons. Novel autonomous exons are more numerous than annotated lncRNA exons, and computational models also indicate they will occur with similar frequency in any randomly generated sequence. These results show that most human coding exons splice autonomously, and provide an explanation for the existence of many unconserved lncRNAs, as well as a new annotation and inclusion levels of spliceable loci in the human genome.

Identification of Novel MET Exon 14 Skipping Variants in Non-Small Cell Lung Cancer Patients: A Prototype Workflow Involving in Silico Prediction and RT-PCR

Background and aims: The MET exon 14 skipping (METex14) is an oncogenic driver mutation that provides a therapeutic opportunity in non-small cell lung cancer (NSCLCs) patients. This event often results from sequence changes at the MET canonical splicing sites. We characterize two novel non-canonical splicing site variants of MET that produce METex14. Materials and Methods: Two variants were identified in three advanced-stage NSCLC patients in a next-generation sequencing panel. The potential impact on splicing was predicted using in silico tools. METex14 mutation was confirmed using reverse transcription (RT)-PCR and a Sanger sequencing analysis on RNA extracted from stained cytology smears. Results: The interrogated MET (RefSeq ID NM_000245.3) variants include a single nucleotide substitution, c.3028+3A>T, in intron 14 and a deletion mutation, c.3012_3028del, in exon 14. The in silico prediction analysis exhibited reduced splicing strength in both variants compared with the MET normal transcript. The RT-PCR and subsequent Sanger sequencing analyses confirmed METex14 skipping in all three patients carrying these variants. Conclusion: This study reveals two non-canonical MET splice variants that cause exon 14 skipping, concurrently also proposes a clinical workflow for the classification of such non-canonical splicing site variants detected by routine DNA-based NGS test. It shows the usefulness of in silico prediction to identify potential METex14 driver mutation and exemplifies the opportunity of routine cytology slides for RNA-based testing.

Stretching muscle cells induces transcriptional and splicing transitions and changes in SR proteins

Alternative splicing is an RNA processing mechanism involved in skeletal muscle development and pathology. Muscular diseases exhibit splicing alterations and changes in mechanobiology leading us to investigate the interconnection between mechanical forces and RNA processing. We performed deep RNA-sequencing after stretching muscle cells. First, we uncovered transcriptional changes in genes encoding proteins involved in muscle function and transcription. Second, we observed that numerous mechanosensitive genes were part of the MAPK pathway which was activated in response to stretching. Third, we revealed that stretching skeletal muscle cells increased the proportion of alternatively spliced cassette exons and their inclusion. Fourth, we demonstrated that the serine and arginine-rich proteins exhibited stronger transcriptional changes than other RNA-binding proteins and that SRSF4 phosphorylation is mechanosensitive. Identifying SRSF4 as a mechanosensitive RNA-binding protein that might contribute to crosstalk between mechanotransduction, transcription, and splicing could potentially reveal novel insights into muscular diseases, particularly those with unknown etiologies.

Ghandhi S, Shimada N, Hubbard K. Identification of RNA Species That Bind to the hnRNP A1 in Normal and Senescent Human Fibroblasts

hnRNP A1 is a member of the hnRNPs (heterogeneous nuclear ribonucleoproteins) family of proteins that play a central role in regulating genes responsible for cell proliferation, DNA repair, apoptosis, and telomere biogenesis. Previous studies have shown that hnRNPA1 had reduced protein levels and increased cytoplasmic accumulation in senescent human diploid fibroblasts. The consequence of reduced protein expression and altered cellular localization may account for the alterations in gene expression observed during senescence. There is limited information for gene targets of hnRNP A1 as well as its in vivo function. In these studies, we performed RNA co-immunoprecipitation experiments using hnRNP A1 as the target protein to identify potential mRNA species in ribonucleoprotein (RNP) complexes. Using this approach, we identified the human double minute 2 (HDM2) mRNA as a binding target for hnRNP A1 in young and senescent human diploid fibroblasts cells. It was also observed that alterations of hnRNP A1 expression modulate HDM2 mRNA levels in young IMR-90 cells. We also demonstrated that the levels of HDM2 mRNA increased with the downregulation of hnRNP A1 and decrease with the overexpression of hnRNP A1. Although we did not observe a significant decrease in HDM2 protein level, a concomitant increase in p53 protein level was detected with the overexpression of hnRNP A1. Our studies also show that hnRNP A1 directly interacts with HDM2 mRNA at a region corresponding to its 3′ UTR (untranslated region of a gene). The results from this study demonstrate that hnRNP A1 has a novel role in participating in the regulation of HDM2 gene expression.

HDAC3 inhibition improves urinary-concentrating defect in hypokalaemia by promoting AQP2 transcription

AIM

This study investigated whether enhanced histone acetylation, achieved by inhibiting histone deacetylases (HDACs), could prevent decreased aquaporin-2 (AQP2) expression during hypokalaemia.

METHODS

Male Wistar rats were fed a potassium-free diet with or without 4-phenylbutyric acid (4-PBA) or the selective HDAC3 inhibitor RGFP966 for 4 days. Primary renal inner medullary collecting duct (IMCD) cells and immortalized mouse cortical collecting duct (mpkCCD) cells were cultured in potassium-deprivation medium with or without HDAC inhibitors.

RESULTS

4-PBA increased the levels of AQP2 mRNA and protein in the kidney inner medullae in hypokalaemic (HK) rats, which was associated with decreased urine output and increased urinary osmolality. The level of acetylated H3K27 (H3K27ac) protein was decreased in the inner medullae of HK rat kidneys; this decrease was mitigated by 4-PBA. The H3K27ac levels were decreased in IMCD and mpkCCD cells cultured in potassium-deprivation medium. Decreased H3K27ac in the Aqp2 promoter region was associated with reduced Aqp2 mRNA levels. HDAC3 protein expression was upregulated in mpkCCD and IMCD cells in response to potassium deprivation, and the binding of HDAC3 to the Aqp2 promoter was also increased. RGFP966 increased the levels of H3K27ac and AQP2 proteins and enhanced binding between H3K27ac and AQP2 in mpkCCD cells. Furthermore, RGFP966 reversed the hypokalaemia-induced downregulation of AQP2 and H3K27ac and alleviated polyuria in rats. RGFP966 increased interstitial osmolality in the kidney inner medullae of HK rats but did not affect urinary cAMP levels.

CONCLUSION

HDAC inhibitors prevented the downregulation of AQP2 induced by potassium deprivation, probably by enhancing H3K27 acetylation.

Genetically programmed alternative splicing of NEMO mediates an autoinflammatory disease phenotype

Host defense and inflammation are regulated by the NF-κB essential modulator (NEMO), a scaffolding protein with a broad immune cell and tissue expression profile. Hypomorphic mutations in inhibitor of NF-κB kinase regulatory subunit gamma (IKBKG) encoding NEMO typically present with immunodeficiency. Here, we characterized a pediatric autoinflammatory syndrome in 3 unrelated male patients with distinct X-linked IKBKG germline mutations that led to overexpression of a NEMO protein isoform lacking the domain encoded by exon 5 (NEMO-Δex5). This isoform failed to associate with TANK binding kinase 1 (TBK1), and dermal fibroblasts from affected patients activated NF-κB in response to TNF but not TLR3 or RIG-I–like receptor (RLR) stimulation when isoform levels were high. By contrast, T cells, monocytes, and macrophages that expressed NEMO-Δex5 exhibited increased NF-κB activation and IFN production, and blood cells from these patients expressed a strong IFN and NF-κB transcriptional signature. Immune cells and TNF-stimulated dermal fibroblasts upregulated the inducible IKK protein (IKKi) that was stabilized by NEMO-Δex5, promoting type I IFN induction and antiviral responses. These data revealed how IKBKG mutations that lead to alternative splicing of skipping exon 5 cause a clinical phenotype we have named NEMO deleted exon 5 autoinflammatory syndrome (NDAS), distinct from the immune deficiency syndrome resulting from loss-of-function IKBKG mutations.

IntSplice2: Prediction of the Splicing Effects of Intronic Single-Nucleotide Variants Using LightGBM Modeling

Prediction of the effect of a single-nucleotide variant (SNV) in an intronic region on aberrant pre-mRNA splicing is challenging except for an SNV affecting the canonical GU/AG splice sites (ss). To predict pathogenicity of SNVs at intronic positions −50 (Int-50) to −3 (Int-3) close to the 3’ ss, we developed light gradient boosting machine (LightGBM)-based IntSplice2 models using pathogenic SNVs in the human gene mutation database (HGMD) and ClinVar and common SNVs in dbSNP with 0.01 ≤ minor allelic frequency (MAF) < 0.50. The LightGBM models were generated using features representing splicing cis-elements. The average recall/sensitivity and specificity of IntSplice2 by fivefold cross-validation (CV) of the training dataset were 0.764 and 0.884, respectively. The recall/sensitivity of IntSplice2 was lower than the average recall/sensitivity of 0.800 of IntSplice that we previously made with support vector machine (SVM) modeling for the same intronic positions. In contrast, the specificity of IntSplice2 was higher than the average specificity of 0.849 of IntSplice. For benchmarking (BM) of IntSplice2 with IntSplice, we made a test dataset that was not used to train IntSplice. After excluding the test dataset from the training dataset, we generated IntSplice2-BM and compared it with IntSplice using the test dataset. IntSplice2-BM was superior to IntSplice in all of the seven statistical measures of accuracy, precision, recall/sensitivity, specificity, F1 score, negative predictive value (NPV), and matthews correlation coefficient (MCC). We made the IntSplice2 web service at https://www.med.nagoya-u.ac.jp/neurogenetics/IntSplice2.

Mutations primarily alter the inclusion of alternatively spliced exons

Genetic analyses and systematic mutagenesis have revealed that synonymous, non-synonymous and intronic mutations frequently alter the inclusion levels of alternatively spliced exons, consistent with the concept that altered splicing might be a common mechanism by which mutations cause disease. However, most exons expressed in any cell are highly-included in mature mRNAs. Here, by performing deep mutagenesis of highly-included exons and by analysing the association between genome sequence variation and exon inclusion across the transcriptome, we report that mutations only very rarely alter the inclusion of highly-included exons. This is true for both exonic and intronic mutations as well as for perturbations in trans. Therefore, mutations that affect splicing are not evenly distributed across primary transcripts but are focussed in and around alternatively spliced exons with intermediate inclusion levels. These results provide a resource for prioritising synonymous and other variants as disease-causing mutations.

A Depletion of Stop Codons in lincRNA is Owing to Transfer of Selective Constraint from Coding Sequences

Although the constraints on a gene’s sequence are often assumed to reflect the functioning of that gene, here we propose transfer selection, a constraint operating on one class of genes transferred to another, mediated by shared binding factors. We show that such transfer can explain an otherwise paradoxical depletion of stop codons in long intergenic noncoding RNAs (lincRNAs). Serine/arginine-rich proteins direct the splicing machinery by binding exonic splice enhancers (ESEs) in immature mRNA. As coding exons cannot contain stop codons in one reading frame, stop codons should be rare within ESEs. We confirm that the stop codon density (SCD) in ESE motifs is low, even accounting for nucleotide biases. Given that serine/arginine-rich proteins binding ESEs also facilitate lincRNA splicing, a low SCD could transfer to lincRNAs. As predicted, multiexon lincRNA exons are depleted in stop codons, a result not explained by open reading frame (ORF) contamination. Consistent with transfer selection, stop codon depletion in lincRNAs is most acute in exonic regions with the highest ESE density, disappears when ESEs are masked, is consistent with stop codon usage skews in ESEs, and is diminished in both single-exon lincRNAs and introns. Owing to low SCD, the maximum lengths of pseudo-ORFs frequently exceed null expectations. This has implications for ORF annotation and the evolution of de novo protein-coding genes from lincRNAs. We conclude that not all constraints operating on genes need be explained by the functioning of the gene but may instead be transferred owing to shared binding factors.

CAGI 5 splicing challenge: Improved exon skipping and intron retention predictions with MMSplice

Pathogenic genetic variants often primarily affect splicing. However, it remains difficult to quantitatively predict whether and how genetic variants affect splicing. In 2018, the fifth edition of the Critical Assessment of Genome Interpretation proposed two splicing prediction challenges based on experimental perturbation assays: Vex-seq, assessing exon skipping, and MaPSy, assessing splicing efficiency. We developed a modular modeling framework, MMSplice, the performance of which was among the best on both challenges. Here we provide insights into the modeling assumptions of MMSplice and its individual modules. We furthermore illustrate how MMSplice can be applied in practice for individual genome interpretation, using the MMSplice VEP plugin and the Kipoi variant interpretation plugin, which are directly applicable to VCF files.

Predicting the change of exon splicing caused by genetic variant using support vector regression

Alternative splicing can be disrupted by genetic variants that are related to diseases like cancers. Discovering the influence of genetic variations on the alternative splicing will improve the understanding of the pathogenesis of variants. Here, we developed a new approach, PredPSI-SVR to predict the impact of variants on exon skipping events by using the support vector regression. From the sequence of a particular exon and its flanking regions, 42 comprehensive features related to splicing events were extracted. By using a greedy feature selection algorithm, we found eight features contributing most to the prediction. The trained model achieved a Pearson correlation coefficient (PCC) of 0.570 in the 10-fold cross-validation based on the training data set provided by the "vex-seq" challenge of the 5th Critical Assessment of Genome Interpretation. In the blind test also held by the challenge, our prediction ranked the 2nd with a PCC of 0.566 that demonstrates the robustness of our method. A further test indicated that the PredPSI-SVR is helpful in prioritizing deleterious synonymous mutations. The method is available on https://github.com/chenkenbio/PredPSI-SVR.

MMSplice: modular modeling improves the predictions of genetic variant effects on splicing

Predicting the effects of genetic variants on splicing is highly relevant for human genetics. We describe the framework MMSplice (modular modeling of splicing) with which we built the winning model of the CAGI5 exon skipping prediction challenge. The MMSplice modules are neural networks scoring exon, intron, and splice sites, trained on distinct large-scale genomics datasets. These modules are combined to predict effects of variants on exon skipping, splice site choice, splicing efficiency, and pathogenicity, with matched or higher performance than state-of-the-art. Our models, available in the repository Kipoi, apply to variants including indels directly from VCF files.

Predicting gene structure changes resulting from genetic variants via exon definition features

Motivation

Genetic variation that disrupts gene function by altering gene splicing between individuals can substantially influence traits and disease. In those cases, accurately predicting the effects of genetic variation on splicing can be highly valuable for investigating the mechanisms underlying those traits and diseases. While methods have been developed to generate high quality computational predictions of gene structures in reference genomes, the same methods perform poorly when used to predict the potentially deleterious effects of genetic changes that alter gene splicing between individuals. Underlying that discrepancy in predictive ability are the common assumptions by reference gene finding algorithms that genes are conserved, well-formed and produce functional proteins.

Results

We describe a probabilistic approach for predicting recent changes to gene structure that may or may not conserve function. The model is applicable to both coding and non-coding genes, and can be trained on existing gene annotations without requiring curated examples of aberrant splicing. We apply this model to the problem of predicting altered splicing patterns in the genomes of individual humans, and we demonstrate that performing gene-structure prediction without relying on conserved coding features is feasible. The model predicts an unexpected abundance of variants that create de novo splice sites, an observation supported by both simulations and empirical data from RNA-seq experiments. While these de novo splice variants are commonly misinterpreted by other tools as coding or non-coding variants of little or no effect, we find that in some cases they can have large effects on splicing activity and protein products and we propose that they may commonly act as cryptic factors in disease.

Availability and implementation

The software is available from geneprediction.org/SGRF.

Supplementary information

Supplementary information is available at Bioinformatics online.

Saturation mutagenesis reveals manifold determinants of exon definition

To illuminate the extent and roles of exonic sequences in the splicing of human RNA transcripts, we conducted saturation mutagenesis of a 51-nt internal exon in a three-exon minigene. All possible single and tandem dinucleotide substitutions were surveyed. Using high-throughput genetics, 5560 minigene molecules were assayed for splicing in human HEK293 cells. Up to 70% of mutations produced substantial (greater than twofold) phenotypes of either increased or decreased splicing. Of all predicted secondary structural elements, only a single 15-nt stem–loop showed a strong correlation with splicing, acting negatively. The in vitro formation of exon-protein complexes between the mutant molecules and proteins associated with spliceosome formation (U2AF35, U2AF65, U1A, and U1-70K) correlated with splicing efficiencies, suggesting exon definition as the step affected by most mutations. The measured relative binding affinities of dozens of human RNA binding protein domains as reported in the CISBP-RNA database were found to correlate either positively or negatively with splicing efficiency, more than could fit on the 51-nt test exon simultaneously. The large number of these functional protein binding correlations point to a dynamic and heterogeneous population of pre-mRNA molecules, each responding to a particular collection of binding proteins.

RNA splicing in human disease and in the clinic

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

Investigating polymorphisms by bioinformatics is a potential cost-effective method to screen for germline mutations in Chinese familial adenomatous polyposis patients

The aim of this study was to investigate germline mutations of the APC, MUTYH and AXIN2 genes in Chinese patients with familial adenomatous polyposis (FAP), and further assess the value of bioinformatics in screening the pathogenic changes predisposing to FAP. APC genes from 11 unrelated FAP patients in Yunnan province in China were firstly examined by exon-specific DNA sequencing. For samples without already known pathogenic changes predisposing to FAP in the APC gene, whole-gene sequencing of MUTYH and AXIN2 was performed. Mutational analysis of each gene was performed by bioinformatics. Eleven different types of APC polymorphisms were observed in the cohort of families analyzed. Of these polymorphisms, four were missense substitutions (V1822D, V1173G, P1760H and K2057), one was a nonsense substitution (S1196X), and six were silent substitutions (Y486Y, T449T, T1493T, G1678G, S1756S and P1960P). One missense mutation (Q335H) and two intronic substitutions (c.264+11G>A and c.420+35A>G) were detected in the MUTYH gene, and four synonymous mutations (I144I, P455P, P462P and L688L) and three intonic mutations (c.1060–77G>T, c.1060–287A>G and c.1060–282 A>G) of the AXIN2 gene were observed. In addition to the already reported pathogenic mutations, by using function assessment tools and databases, the synonymous substitutions observed in the APC gene of our samples were predicted to affect splicing regulation in the translation of mRNA, while the missense mutations observed in the APC gene and MUTYH gene were predicted to be disease-related polymorphisms; however, no functional effect of the mutations was observed in the AXIN2 gene. Comprehensive screening for germline mutations in APC, MUTYH and AXIN2 genes followed by prediction of pathogenicity using bioinformatic tools contributes to a cost-effective way of screening germline mutations in Chinese familial adenomatous polyposis patients.

IntSplice: prediction of the splicing consequences of intronic single-nucleotide variations in the human genome

Precise spatiotemporal regulation of splicing is mediated by splicing cis-elements on pre-mRNA. Single-nucleotide variations (SNVs) affecting intronic cis-elements possibly compromise splicing, but no efficient tool has been available to identify them. Following an effect-size analysis of each intronic nucleotide on annotated alternative splicing, we extracted 105 parameters that could affect the strength of the splicing signals. However, we could not generate reliable support vector regression models to predict the percent-splice-in (PSI) scores for normal human tissues. Next, we generated support vector machine (SVM) models using 110 parameters to directly differentiate pathogenic SNVs in the Human Gene Mutation Database and normal SNVs in the dbSNP database, and we obtained models with a sensitivity of 0.800±0.041 (mean and s.d.) and a specificity of 0.849±0.021. Our IntSplice models were more discriminating than SVM models that we generated with Shapiro-Senapathy score and MaxEntScan::score3ss. We applied IntSplice to a naturally occurring and nine artificial intronic mutations in RAPSN causing congenital myasthenic syndrome. IntSplice correctly predicted the splicing consequences for nine of the ten mutants. We created a web service program, IntSplice (http://www.med.nagoya-u.ac.jp/neurogenetics/IntSplice) to predict splicing-affecting SNVs at intronic positions from -50 to -3.

CoSREM: a graph mining algorithm for the discovery of combinatorial splicing regulatory elements

BACKGROUND

Alternative splicing (AS) is a post-transcriptional regulatory mechanism for gene expression regulation. Splicing decisions are affected by the combinatorial behavior of different splicing factors that bind to multiple binding sites in exons and introns. These binding sites are called splicing regulatory elements (SREs). Here we develop CoSREM (Combinatorial SRE Miner), a graph mining algorithm to discover combinatorial SREs in human exons. Our model does not assume a fixed length of SREs and incorporates experimental evidence as well to increase accuracy. CoSREM is able to identify sets of SREs and is not limited to SRE pairs as are current approaches.

RESULTS

We identified 37 SRE sets that include both enhancer and silencer elements. We show that our results intersect with previous results, including some that are experimental. We also show that the SRE set GGGAGG and GAGGAC identified by CoSREM may play a role in exon skipping events in several tumor samples. We applied CoSREM to RNA-Seq data for multiple tissues to identify combinatorial SREs which may be responsible for exon inclusion or exclusion across tissues.

CONCLUSION

The new algorithm can identify different combinations of splicing enhancers and silencers without assuming a predefined size or limiting the algorithm to find only pairs of SREs. Our approach opens new directions to study SREs and the roles that AS may play in diseases and tissue specificity.

Identifying splicing regulatory elements with de Bruijn graphs

Splicing regulatory elements (SREs) are short, degenerate sequences on pre-mRNA molecules that enhance or inhibit the splicing process via the binding of splicing factors, proteins that regulate the functioning of the spliceosome. Existing methods for identifying SREs in a genome are either experimental or computational. Here, we propose a formalism based on de Bruijn graphs that combines genomic structure, word count enrichment analysis, and experimental evidence to identify SREs found in exons. In our approach, SREs are not restricted to a fixed length (i.e., k-mers, for a fixed k). As a result, we identify 2001 putative exonic enhancers and 3080 putative exonic silencers for human genes, with lengths varying from 6 to 15 nucleotides. Many of the predicted SREs overlap with experimentally verified binding sites. Our model provides a novel method to predict variable length putative regulatory elements computationally for further experimental investigation.

129-Derived Mouse Strains Express an Unstable but Catalytically Active DNA Polymerase Iota Variant

Mice derived from the 129 strain have a nonsense codon mutation in exon 2 of the polymerase iota (Polι) gene and are therefore considered Polι deficient. When we amplified Polι mRNA from 129/SvJ or 129/Ola testes, only a small fraction of the full-length cDNA contained the nonsense mutation; the major fraction corresponded to a variant Polι isoform lacking exon 2. Polι mRNA lacking exon 2 contains an open reading frame, and the corresponding protein was detected using a polyclonal antibody raised against the C terminus of the murine Polι protein. The identity of the corresponding protein was further confirmed by mass spectrometry. Although the variant protein was expressed at only 5 to 10% of the level of wild-type Polι, it retained de novo DNA synthesis activity, the capacity to form replication foci following UV irradiation, and the ability to rescue UV light sensitivity in Polι(-/-) embryonic fibroblasts derived from a new, fully deficient Polι knockout (KO) mouse line. Furthermore, in vivo treatment of 129-derived male mice with Velcade, a drug that inhibits proteasome function, stabilized and restored a substantial amount of the variant Polι in these animals, indicating that its turnover is controlled by the proteasome. An analysis of two xeroderma pigmentosum-variant (XPV) cases corresponding to missense mutants of Polη, a related translesion synthesis (TLS) polymerase in the same family, similarly showed a destabilization of the catalytically active mutant protein by the proteasome. Collectively, these data challenge the prevailing hypothesis that 129-derived strains of mice are completely deficient in Polι activity. The data also document, both for 129-derived mouse strains and for some XPV patients, new cases of genetic defects corresponding to the destabilization of an otherwise functional protein, the phenotype of which is reversible by proteasome inhibition.

Dissecting alternative splicing in the formation of Miltenberger glycophorin subtype III (GYP.Mur)

BACKGROUND AND OBJECTIVES

Miltenberger subtype III (Mi.III, GP.Mur) is one of the most important red cell phenotypes in the fields of transfusion in South-East Asia. GP.Mur is believed to evolve from homologous gene recombination events between glycophorin A (GYPA) and glycophorin B (GYPB). GYP.Mur differs from GYPB in only seven nucleotides dispersed near the region of 3' exon 3 of GYP.Mur. The goal of this study was to dissect how these nucleotide variants affected splicing of exon 3.

MATERIALS AND METHODS

We first designed two minigene constructs: one containing GYP.Mur from exon 2 to exon 4 and the other containing GYPB in the same region. To test how these nucleotide variations between GYP.Mur and GYPB affected the splicing, a repertoire of the GYP.Mur-like minigene constructs with different point mutations were created. These minigene variants were evaluated for their abilities to induce splicing of exon 3 using a heterologous expression system.

RESULTS

(1) GYP.Mur minigene expressed exons 2, 3 and 4, whereas GYPB minigene expressed only exon 2 and exon 4. (2) The single nucleotide alteration at the position of the 5' splice site of glycophorin intron 3 reversed the splicing decision. (3) The nucleotide variations between GYP.Mur and GYPB other than that at the 5' splice site showed very little or no effect on splicing of exon 3.

CONCLUSION

Splicing of the glycophorin B-A-B hybrids (GYP.Mur and GYP.BUN) and unsplicing of GYPB follow the GU-AG rule strictly.

Functional annotation of putative regulatory elements at cancer susceptibility Loci

Most cancer-associated genetic variants identified from genome-wide association studies (GWAS) do not obviously change protein structure, leading to the hypothesis that the associations are attributable to regulatory polymorphisms. Translating genetic associations into mechanistic insights can be facilitated by knowledge of the causal regulatory variant (or variants) responsible for the statistical signal. Experimental validation of candidate functional variants is onerous, making bioinformatic approaches necessary to prioritize candidates for laboratory analysis. Thus, a systematic approach for recognizing functional (and, therefore, likely causal) variants in noncoding regions is an important step toward interpreting cancer risk loci. This review provides a detailed introduction to current regulatory variant annotations, followed by an overview of how to leverage these resources to prioritize candidate functional polymorphisms in regulatory regions.

Absence of an intron splicing silencer in porcine Smn1 intron 7 confers immunity to the exon skipping mutation in human SMN2

Spinal Muscular Atrophy is caused by homozygous loss of SMN1. All patients retain at least one copy of SMN2 which produces an identical protein but at lower levels due to a silent mutation in exon 7 which results in predominant exclusion of the exon. Therapies targeting the splicing of SMN2 exon 7 have been in development for several years, and their efficacy has been measured using either in vitro cellular assays or in vivo small animal models such as mice. In this study we evaluated the potential for constructing a mini-pig animal model by introducing minimal changes in the endogenous porcine Smn1 gene to maintain the native genomic structure and regulation. We found that while a Smn2-like mutation can be introduced in the porcine Smn1 gene and can diminish the function of the ESE, it would not recapitulate the splicing pattern seen in human SMN2 due to absence of a functional ISS immediately downstream of exon 7. We investigated the ISS region and show here that the porcine ISS is inactive due to disruption of a proximal hnRNP A1 binding site, while a distal hnRNP A1 binding site remains functional but is unable to maintain the functionality of the ISS as a whole.

Exon identity crisis: disease-causing mutations that disrupt the splicing code

Cis-acting RNA elements control the accurate expression of human multi-exon protein coding genes. Single nucleotide variants altering the fidelity of this regulatory code and, consequently, pre-mRNA splicing are expected to contribute to the etiology of numerous human diseases.

The evolution, impact and properties of exonic splice enhancers

BACKGROUND

In humans, much of the information specifying splice sites is not at the splice site. Exonic splice enhancers are one of the principle non-splice site motifs. Four high-throughput studies have provided a compendium of motifs that function as exonic splice enhancers, but only one, RESCUE-ESE, has been generally employed to examine the properties of enhancers. Here we consider these four datasets to ask whether there is any consensus on the properties and impacts of exonic splice enhancers.

RESULTS

While only about 1% of all the identified hexamer motifs are common to all analyses we can define reasonably sized sets that are found in most datasets. These consensus intersection datasets we presume reflect the true properties of exonic splice enhancers. Given prior evidence for the properties of enhancers and splice-associated mutations, we ask for all datasets whether the exonic splice enhancers considered are purine enriched; enriched near exon boundaries; able to predict trends in relative codon usage; slow evolving at synonymous sites; rare in SNPs; associated with weak splice sites; and enriched near longer introns. While the intersect datasets match expectations, only one original dataset, RESCUE-ESE, does. Unexpectedly, a fully experimental dataset identifies motifs that commonly behave opposite to the consensus, for example, being enriched in exon cores where splice-associated mutations are rare.

CONCLUSIONS

Prior analyses that used the RESCUE-ESE set of hexamers captured the properties of consensus exonic splice enhancers. We estimate that at least 4% of synonymous mutations are deleterious owing to an effect on enhancer functioning.

In silico to in vivo splicing analysis using splicing code models

With the growing appreciation of RNA splicing's role in gene regulation, development, and disease, researchers from diverse fields find themselves investigating exons of interest. Commonly, researchers are interested in knowing if an exon is alternatively spliced, if it is differentially included in specific tissues or in developmental stages, and what regulatory elements control its inclusion. An important step towards the ability to perform such analysis in silico was made with the development of computational splicing code models. Aimed as a practical how-to guide, we demonstrate how researchers can now use these code models to analyze a gene of interest, focusing on Bin1 as a case study. Bridging integrator 1 (BIN1) is a nucleocytoplasmic adaptor protein known to be functionally regulated through alternative splicing in a tissue-specific manner. Specific Bin1 isoforms have been associated with muscular diseases and cancers, making the study of its splicing regulation of wide interest. Using AVISPA, a recently released web tool based on splicing code models, we show that many Bin1 tissue-dependent isoforms are correctly predicted, along with many of its known regulators. We review the best practices and constraints of using the tool, demonstrate how AVISPA is used to generate high confidence novel regulatory hypotheses, and experimentally validate predicted regulators of Bin1 alternative splicing.

Novel GUCA1A mutations suggesting possible mechanisms of pathogenesis in cone, cone-rod, and macular dystrophy patients

Here, we report two novel GUCA1A (the gene for guanylate cyclase activating protein 1) mutations identified in unrelated Spanish families affected by autosomal dominant retinal degeneration (adRD) with cone and rod involvement. All patients from a three-generation adRD pedigree underwent detailed ophthalmic evaluation. Total genome scan using single-nucleotide polymorphisms and then the linkage analysis were undertaken on the pedigree. Haplotype analysis revealed a 55.37 Mb genomic interval cosegregating with the disease phenotype on chromosome 6p21.31-q15. Mutation screening of positional candidate genes found a heterozygous transition c.250C>T in exon 4 of GUCA1A, corresponding to a novel mutation p.L84F. A second missense mutation, c.320T>C (p.I107T), was detected by screening of the gene in a Spanish patients cohort. Using bioinformatics approach, we predicted that either haploinsufficiency or dominant-negative effect accompanied by creation of a novel function for the mutant protein is a possible mechanism of the disease due to c.250C>T and c.320T>C. Although additional functional studies are required, our data in relation to the c.250C>T mutation open the possibility that transacting factors binding to de novo created recognition site resulting in formation of aberrant splicing variant is a disease model which may be more widespread than previously recognized as a mechanism causing inherited RD.

Leucine to proline substitution by SNP at position 197 in Caspase-9 gene expression leads to neuroblastoma: a bioinformatics analysis

To understand the role of CASP9 (Caspase-9) gene products in relation to neuroblastoma disease, we have analyzed the single nucleotide polymorphisms (SNPs) associated with this gene. This can help us understand the genetic variations that can alter the function of the gene products. A total of 941 SNPs are investigated for CASP9 gene. To determine whether a non-synonymous SNP (nsSNP) in this gene affects its protein product, we used certain computational tools which predicted one nsSNP, rs1052574, to have deleterious phenotypic effect. This polymorphic variant results in amino acid substitution from leucine to proline at 197 position, i.e., from acyclic amino acid to a 5-membered amino acid which resides in the buried area of the protein with a high level of conservation. This amino acid substitution shows a transition from helix to coil in the mutant protein. Hence, due to the complete alteration in the structural property of the amino acid side chain, the stability of the protein is reduced which may affect the function of CASP9 protein, leading to deregulation of apoptosis and neuroblastoma development.

Refining the role of PMS2 in Lynch syndrome: germline mutational analysis improved by comprehensive assessment of variants

BACKGROUND AND AIM

The majority of mismatch repair (MMR) gene mutations causing Lynch syndrome (LS) occur either in MLH1 or MSH2. However, the relative contribution of PMS2 is less well defined. The aim of this study was to evaluate the role of PMS2 in LS by assessing the pathogenicity of variants of unknown significance (VUS) detected in the mutational analysis of PMS2 in a series of Spanish patients.

METHODS

From a cohort of 202 LS suspected patients, 13 patients showing loss of PMS2 expression in tumours were screened for germline mutations in PMS2, using a long range PCR based strategy and multiplex ligation dependent probe amplification (MLPA). Pathogenicity assessment of PMS2 VUS was performed evaluating clinicopathological data, frequency in control population and in silico and in vitro analyses at the RNA and protein level.

RESULTS

Overall 25 different PMS2 DNA variants were detected. Fourteen were classified as polymorphisms. Nine variants were classified as pathogenic: seven alterations based on their molecular nature and two after demonstrating a functional defect (c.538-3C>G affected mRNA processing and c.137G>T impaired MMR activity). The c.1569C>G variant was classified as likely neutral while the c.384G>A remained as a VUS. We have also shown that the polymorphic variant c.59G>A is MMR proficient.

CONCLUSIONS

Pathogenic PMS2 mutations were detected in 69% of patients harbouring LS associated tumours with loss of PMS2 expression. In all, PMS2 mutations account for 6% of the LS cases identified. The comprehensive functional analysis shown here has been useful in the classification of PMS2 VUS and contributes to refining the role of PMS2 in LS.

BRCA1 polymorphisms and breast cancer epidemiology in the Western New York exposures and breast cancer (WEB) study

Results of studies for the association of BRCA1 genotypes and haplotypes with sporadic breast cancer have been inconsistent. Therefore, a candidate single nucleotide polymorphism (SNP) approach was used in a breast cancer case-control study to explore genotypes and haplotypes that have the potential to affect protein functions or levels. In a breast cancer case-control study, genotyping of BRCA1 polymorphisms Q356R, D693N, and E1038G was performed on 1,005 cases and 1,765 controls. Unconditional, polytomous logistic regression and χ(2) -tests were used to examine the associations of breast cancer with genotypes and haplotypes. In addition, interactions between genotype and smoking, benign breast disease, family history of breast cancer, body mass index (BMI), alcohol consumption, and hormonal risk factors, hormone receptor status, and breast cancer pathology were calculated also using logistic regression and χ(2) . Although sporadic breast cancer was not associated with BRCA1 genotypes or haplotypes overall or by menopausal status, there was evidence of an interaction between the E1038G BRCA1 genotype, smoking, and BMI among premenopausal women (P for interaction = 0.01 and 0.045, respectively) and between E1038G and D693N BRCA1 genotypes and hormone therapy use among postmenopausal women (P for interaction = 0.01 and 0.02, respectively). There were no other associations found between BRCA1 genotypes and stage, histological grade, or nuclear grade. However, the D693N SNP was associated with the risk of triple negative breast cancer (odds ratio = 2.31 95% confidence interval 1.08-4.93). The BRCA1 variants studied may play a role in the etiology of triple negative breast cancer and may interact with environmental factors such as hormone therapy or smoking and increase sporadic breast cancer risk.

The mechanism of alternative splicing of the X-linked NDUFB11 gene of the respiratory chain complex I, impact of rotenone treatment in neuroblastoma cells

A study is presented on the regulation of alternative splicing (AS) of the Ndufb11 gene of complex I of the mitochondrial respiratory chain and the impact on this process of rotenone treatment in neuroblastoma cells. In physiological conditions the Ndufb11 gene produces at high level a short transcript isoform encoding for a 153 aa protein. This subunit is essential for the assembly of a functional and stable mammalian complex I. The gene produces also, at low level, a longer transcript isoform encoding for a 163 aa protein whose role is unknown. Evidence is presented here showing that the level of the two isoforms is regulated by three DGGGD ESS elements located in exon 2 which can bind the hnRNPH1 protein. In neuronal cells rotenone treatment affects the Ndufb11 alternative splicing pathway, with the increase of the 163/153 mRNAs ratio. This effect appears to be due to the down-regulation of the hnRNPH1 protein. Since rotenone induces apoptosis in neuronal cells, the post-transcriptional regulation of the Ndufb11 gene can be involved in the programmed cell death process.

An exon splice enhancer primes IGF2:IGF2R binding site structure and function evolution

Placental development and genomic imprinting coevolved with parental conflict over resource distribution to mammalian offspring. The imprinted genes IGF2 and IGF2R code for the growth promoter insulin-like growth factor 2 (IGF2) and its inhibitor, mannose 6-phosphate (M6P)/IGF2 receptor (IGF2R), respectively. M6P/IGF2R of birds and fish do not recognize IGF2. In monotremes, which lack imprinting, IGF2 specifically bound M6P/IGF2R via a hydrophobic CD loop. We show that the DNA coding the CD loop in monotremes functions as an exon splice enhancer (ESE) and that structural evolution of binding site loops (AB, HI, FG) improved therian IGF2 affinity. We propose that ESE evolution led to the fortuitous acquisition of IGF2 binding by M6P/IGF2R that drew IGF2R into parental conflict; subsequent imprinting may then have accelerated affinity maturation.

Analysis of the effects of rare variants on splicing identifies alterations in GABAA receptor genes in autism spectrum disorder individuals

A large-scale sequencing screen of X-linked synaptic genes in individuals with autism spectrum disorder (ASD) or schizophrenia (SCZ), two common neurodevelopmental disorders, identified many variants most of which have no easily predictable effect on gene function. In this report, we evaluated the impact of these rare missense and silent variants on gene splicing. For this purpose, we used complementary in silico analyses, in vitro minigene-based assays and RNA prepared from lymphoblastoid cells derived from patients with these mutations. Our goal was to identify the variants which might either create or disrupt an acceptor splice site, a donor splice site or an exonic splicing enhancer, thus leading to aberrant splicing that could be involved in the pathogenesis of ASD or SCZ. We identified truncating mutations in distinct X-linked gamma-aminobutyric acid A (GABAA) receptor subunit-encoding genes, GABRQ and GABRA3, in two different families. Furthermore, missense and silent variants in nuclear RNA export factor 5 and histone deacetylase 6 were shown to partially disrupt the protein. While genes from the GABAergic pathway have previously been thought to be involved in the pathophysiology of ASD, this is the first report of ASD patients with truncating mutations in GABA receptors genes.

Comprehensive functional assessment of MLH1 variants of unknown significance

Lynch syndrome is associated with germline mutations in DNA mismatch repair (MMR) genes. Up to 30% of DNA changes found are variants of unknown significance (VUS). Our aim was to assess the pathogenicity of eight MLH1 VUS identified in patients suspected of Lynch syndrome. All of them are novel or not previously characterized. For their classification, we followed a strategy that integrates family history, tumor pathology, and control frequency data with a variety of in silico and in vitro analyses at RNA and protein level, such as MMR assay, MLH1 and PMS2 expression, and subcellular localization. Five MLH1 VUS were classified as pathogenic: c.[248G>T(;)306G>C], c.[780C>G;788A>C], and c.791-7T>A affected mRNA processing, whereas c.218T>C (p.L73P) and c.244A>G [corrected] (p.T82A) impaired MMR activity. Two other VUS were considered likely neutral: the silent c.702G>A variant did not affect mRNA processing or stability, and c.974G>A (p.R325Q) did not influence MMR function. In contrast, variant c.25C>T (p.R9W) could not be classified, as it associated with intermediate levels of MMR activity. Comprehensive functional assessment of MLH1 variants was useful in their classification and became relevant in the diagnosis and genetic counseling of carrier families.

Identification of evolutionarily conserved exons as regulated targets for the splicing activator tra2β in development

Alternative splicing amplifies the information content of the genome, creating multiple mRNA isoforms from single genes. The evolutionarily conserved splicing activator Tra2β (Sfrs10) is essential for mouse embryogenesis and implicated in spermatogenesis. Here we find that Tra2β is up-regulated as the mitotic stem cell containing population of male germ cells differentiate into meiotic and post-meiotic cells. Using CLIP coupled to deep sequencing, we found that Tra2β binds a high frequency of exons and identified specific G/A rich motifs as frequent targets. Significantly, for the first time we have analysed the splicing effect of Sfrs10 depletion in vivo by generating a conditional neuronal-specific Sfrs10 knock-out mouse (Sfrs10(fl/fl); Nestin-Cre(tg/+)). This mouse has defects in brain development and allowed correlation of genuine physiologically Tra2β regulated exons. These belonged to a novel class which were longer than average size and importantly needed multiple cooperative Tra2β binding sites for efficient splicing activation, thus explaining the observed splicing defects in the knockout mice. Regulated exons included a cassette exon which produces a meiotic isoform of the Nasp histone chaperone that helps monitor DNA double-strand breaks. We also found a previously uncharacterised poison exon identifying a new pathway of feedback control between vertebrate Tra2 proteins. Both Nasp-T and the Tra2a poison exon are evolutionarily conserved, suggesting they might control fundamental developmental processes. Tra2β protein isoforms lacking the RRM were able to activate specific target exons indicating an additional functional role as a splicing co-activator. Significantly the N-terminal RS1 domain conserved between flies and humans was essential for the splicing activator function of Tra2β. Versions of Tra2β lacking this N-terminal RS1 domain potently repressed the same target exons activated by full-length Tra2β protein.

Mutations in SDCCAG8/NPHP10 Cause Bardet-Biedl Syndrome and Are Associated with Penetrant Renal Disease and Absent Polydactyly

The ciliopathies are an expanding group of disorders caused by mutations in genes implicated in the biogenesis and function of primary cilia. Bardet-Biedl syndrome (BBS) is a model ciliopathy characterized by progressive retinal degeneration, obesity, polydactyly, cognitive impairment, kidney anomalies and hypogonadism. Mutations in SDCCAG8(NPHP10) were described recently in patients with nephronophthisis and retinal degeneration (Senior-Loken syndrome; SLS). Given the phenotypic and genetic overlap between known ciliopathy genes, we hypothesized that mutations in SDCCAG8 might also contribute alleles to more severe, multisystemic ciliopathies. We performed genetic and phenotypic analyses of 2 independent BBS cohorts. Subsequent to mutation screening, we made a detailed phenotypic analysis of 5 families mutated for SDCCAG8 (3 homozygous and 2 compound heterozygous mutations) and conducted statistical analyses across both cohorts to examine possible phenotype-genotype correlations with mutations at this locus. All patients with mutations in SDCCAG8 fulfilled the diagnostic criteria for BBS (retinal degeneration, obesity, cognitive defects, renal failure, hypogonadism). Interestingly, none of the patients with primary SDCCAG8 mutations had polydactyly, a frequent but not obligatory BBS feature. In contrast, the same patients displayed early-onset renal failure, obesity, as well as recurrent pulmonary and ENT infections. Comparison of the phenotypes of these families with our entire BBS cohort indicated that renal impairment and absent polydactyly correlated significantly with causal SDCCAG8 mutations. Thus, SDCCAG8 mutations are sufficient to cause BBS in 1-2% of our combined cohorts, and define this gene as the sixteenth BBS locus (BBS16). The absence of polydactyly and the concomitant, apparently fully penetrant association with early kidney failure represents the first significant genotype-phenotype correlation in BBS that potentially represents an indicator for phenotype-driven priority screening and informs specific patient management.

Evolution of SR protein and hnRNP splicing regulatory factors

The splicing of pre-mRNAs is an essential step of gene expression in eukaryotes. Introns are removed from split genes through the activities of the spliceosome, a large ribonuclear machine that is conserved throughout the eukaryotic lineage. While unicellular eukaryotes are characterized by less complex splicing, pre-mRNA splicing of multicellular organisms is often associated with extensive alternative splicing that significantly enriches their proteome. The alternative selection of splice sites and exons permits multicellular organisms to modulate gene expression patterns in a cell type-specific fashion, thus contributing to their functional diversification. Alternative splicing is a regulated process that is mainly influenced by the activities of splicing regulators, such as SR proteins or hnRNPs. These modular factors have evolved from a common ancestor through gene duplication events to a diverse group of splicing regulators that mediate exon recognition through their sequence-specific binding to pre-mRNAs. Given the strong correlations between intron expansion, the complexity of pre-mRNA splicing, and the emergence of splicing regulators, it is argued that the increased presence of SR and hnRNP proteins promoted the evolution of alternative splicing through relaxation of the sequence requirements of splice junctions.

Novel DMRT1 3'UTR+11insT mutation associated to XY partial gonadal dysgenesis

The Y-chromosome-located SRY gene encodes a small testis-specific protein containing a DNA-binding motif known as the HMG (high mobility group) box. However, mutations in SRY are not frequent especially in cases of 46,XY partial gonadal dysgenesis. Several sex-determining genes direct the fate of the bipotential gonad to either testis or ovary. In addition, heterozygous small deletions in 9p can cause complete and partial XY gonadal dysgenesis without other symptoms. Human DMRT1 gene, which is located at 9p24.3, is expressed in testis and ovary and has been considered, among others, a candidate autosomal gene responsible for gonadal dysgenesis. In this report we describe a nucleotide insertion in DMRT1 3'UTR in a patient of XY partial gonadal dygenesis. The 3'UTR+11insT is located within a conserved motif important for mRNA stabilization.

The evolution of alternative splicing exons in vascular endothelial growth factor A

The C-terminus alternative splicing in VEGFA (vascular endothelial growth factor A) is known for its impact on physiological and pathological angiogenesis. Based on our prediction and RT-PCR verification, we identified anti-angiogenic VEGFA165b isoforms in mouse and rabbit for the first time. We also found that the relative expression level of VEGFA165b isoform had been increasing from rodents to human, and exon8b may have experienced a minor-to-major form exon conversion, possibly correlated with its gain-of-function. It is suggested that introduction of alternative splicing exons (esp. exon6 and exon8b) made important contributions to the transcriptional diversity of VEGFA and played a crucial role in the evolution of its regulatory mechanism.

Molecular design of a splicing switch responsive to the RNA binding protein Tra2β

Tra2β regulates a number of splicing switches including activation of the human testis-specific exon HIPK3-T in the Homeodomain Interacting Protein Kinase 3 gene. By testing HIPK3-T exons of different intrinsic strengths, we found Tra2β most efficiently activated splicing inclusion of intrinsically weak exons, although these were spliced at a lower overall level. Both the RRM and N-terminal RS-rich region of Tra2β were required for splicing activation. Bioinformatic searches for splicing enhancers and repressors mapped four physically distinct exonic splicing enhancers (ESEs) within HIPK3-T, each containing the known Tra2β AGAA-rich binding site. Surprisingly disruption of each single ESE prevented Tra2β-mediated activation, although single mutated exons could still bind Tra2β protein by gel shifts and functional splicing analyses. Titration experiments indicate an additive model of HIPK3-T splicing activation, requiring availability of an array of four distinct ESEs to enable splicing activation. To enable this efficient Tra2β-mediated splicing switch to operate, a closely adjacent downstream and potentially competitive stronger 5'-splice site is actively repressed. Our data indicate that a novel arrangement of multiple mono-specific AGAA-rich ESEs coupled to a weak 5'-splice site functions as a responsive gauge. This gauge monitors changes in the specific nuclear concentration of the RNA binding protein Tra2β, and co-ordinately regulates HIPK3-T exon splicing inclusion.

Using positional distribution to identify splicing elements and predict pre-mRNA processing defects in human genes

We present an intuitive strategy for predicting the effect of sequence variation on splicing. In contrast to transcriptional elements, splicing elements appear to be strongly position dependent. We demonstrated that exonic binding of the normally intronic splicing factor, U2AF65, inhibits splicing. Reasoning that the positional distribution of a splicing element is a signature of its function, we developed a method for organizing all possible sequence motifs into clusters based on the genomic profile of their positional distribution around splice sites. Binding sites for serine/arginine rich (SR) proteins tended to be exonic whereas heterogeneous ribonucleoprotein (hnRNP) recognition elements were mostly intronic. In addition to the known elements, novel motifs were returned and validated. This method was also predictive of splicing mutations. A mutation in a motif creates a new motif that sometimes has a similar distribution shape to the original motif and sometimes has a different distribution. We created an intraallelic distance measure to capture this property and found that mutations that created large intraallelic distances disrupted splicing in vivo whereas mutations with small distances did not alter splicing. Analyzing the dataset of human disease alleles revealed known splicing mutants to have high intraallelic distances and suggested that 22% of disease alleles that were originally classified as missense mutations may also affect splicing. This category together with mutations in the canonical splicing signals suggest that approximately one third of all disease-causing mutations alter pre-mRNA splicing.