AG-dependent 3'-splice sites are predisposed to aberrant splicing due to a mutation at the first nucleotide of an exon

FexSplice: A LightGBM-Based Model for Predicting the Splicing Effect of a Single Nucleotide Variant Affecting the First Nucleotide G of an Exon

Single nucleotide variants (SNVs) affecting the first nucleotide G of an exon (Fex-SNVs) identified in various diseases are mostly recognized as missense or nonsense variants. Their effect on pre-mRNA splicing has been seldom analyzed, and no curated database is available. We previously reported that Fex-SNVs affect splicing when the length of the polypyrimidine tract is short or degenerate. However, we cannot readily predict the splicing effects of Fex-SNVs. We here scrutinized the available literature and identified 106 splicing-affecting Fex-SNVs based on experimental evidence. We similarly identified 106 neutral Fex-SNVs in the dbSNP database with a global minor allele frequency (MAF) of more than 0.01 and less than 0.50. We extracted 115 features representing the strength of splicing cis-elements and developed machine-learning models with support vector machine, random forest, and gradient boosting to discriminate splicing-affecting and neutral Fex-SNVs. Gradient boosting-based LightGBM outperformed the other two models, and the length and nucleotide compositions of the polypyrimidine tract played critical roles in the discrimination. Recursive feature elimination showed that the LightGBM model using 15 features achieved the best performance with an accuracy of 0.80 ± 0.12 (mean and SD), a Matthews Correlation Coefficient (MCC) of 0.57 ± 0.15, an area under the curve of the receiver operating characteristics curve (AUROC) of 0.86 ± 0.08, and an area under the curve of the precision-recall curve (AUPRC) of 0.87 ± 0.09 using a 10-fold cross-validation. We developed a web service program, named FexSplice that accepts a genomic coordinate either on GRCh37/hg19 or GRCh38/hg38 and returns a predicted probability of aberrant splicing of A, C, and T variants.

Multiomics analysis of the NAD+-PARP1 axis reveals a role for site-specific ADP-ribosylation in splicing in embryonic stem cells

The differentiation of embryonic stem cells (ESCs) into a lineage-committed state is a dynamic process involving changes in cellular metabolism, epigenetic modifications, post-translational modifications, gene expression, and RNA processing. Here we integrated data from metabolomic, proteomic, and transcriptomic assays to characterize how alterations in NAD⁺ metabolism during the differentiation of mouse ESCs lead to alteration of the PARP1-mediated ADP-ribosylated (ADPRylated) proteome and mRNA isoform specialization. Our metabolomic analyses indicate that mESCs use distinct NAD⁺ biosynthetic pathways in different cell states: the de novo pathway in the pluripotent state, and the salvage and Preiss-Handler pathways as differentiation progresses. We observed a dramatic induction of PARP1 catalytic activity driven by enhanced nuclear NAD⁺ biosynthesis during the early stages of mESC differentiation (e.g., within 12 h of LIF removal). PARP1-modified proteins in mESCs are enriched for biological processes related to stem cell maintenance, transcriptional regulation, and RNA processing. The PARP1 substrates include core spliceosome components, such as U2AF35 and U2AF65, whose splicing functions are modulated by PARP1-mediated site-specific ADP-ribosylation. Finally, we observed that splicing is dysregulated genome-wide in Parp1 knockout mESCs. Together, these results demonstrate a role for the NAD⁺-PARP1 axis in the maintenance of mESC state, specifically in the splicing program during differentiation.

Study of novel NARS2 variants in patient of combined oxidative phosphorylation deficiency 24

BACKGROUND

NARS2 catalyzes the attachment of asparagine amino acids to mitochondrial tRNA^Asn and is critical for efficient mitochondrial protein synthesis. Biallelic variants in NARS2 are associated with combined oxidative phosphorylation deficiency 24 (COXPD24) and autosomal recessive deafness-94.

METHODS

Patient information was obtained after recruitment. Genetic tests were performed using whole exome sequencing (WES) and Sanger sequencing. Structure prediction was based on the RaptorX and SWISS-MODEL platforms. The mRNA analysis of paternal variant was performed. Expression levels and dimerization of wild-type (WT) and mutant NARS2 were detected in human embryonic kidney (HEK) 293T cells. Mitochondrial localization of NARS2 variants was determined using immunofluorescence staining.

RESULTS

The patient presented early onset generalized epilepsy, myoclonic seizures, severe bilateral hearing impairment and affected liver and heart. WES identified two compound heterozygous variants in NARS2: c.1141A>G/p.Asn381Asp and c.1290G>C/p.Trp430Cys. In silico analysis predicted that both variants would cause significant and pathogenic changes in secondary structure. NARS2 c.1290G>C is a variant at the first nucleotide of an exon, a location thought to affect mRNA splicing. Although transcriptional experiments did not identify aberrant splicing, we observed a lower proportion of transcripts from the NARS2 c.1290G>C variant. An in vitro experiment showed that both variants impaired NARS2 expression, while mitochondrial localization and dimerization remained unaffected.

CONCLUSIONS

The patient was diagnosed with COXPD24 caused by novel NARS2 variations. The cardiac dysfunction is identified for the first time. In vitro study revealed impairment of variants on NARS2 expression. These data enrich our knowledge regarding the phenotypic and genotypic spectra of NARS2.

AKR1A1 Variant Associated With Schizophrenia Causes Exon Skipping, Leading to Loss of Enzymatic Activity

Schizophrenia is a heterogeneous psychiatric disorder characterized by positive symptoms such as hallucinations and delusions, negative symptoms such as anhedonia and flat affect, and cognitive impairment. Recently, glucuronate (GlucA) levels were reported to be significantly higher in serum of patients with schizophrenia than those in healthy controls. The accumulation of GlucA is known to be related to treatment-resistant schizophrenia, since GlucA is known to promote drug excretion by forming conjugates with drugs. However, the cause of GlucA accumulation remains unclear. Aldo-keto reductase family one member A1 (AKR1A1) is an oxidoreductase that catalyzes the reduction of GlucA. Genetic loss of AKR1A1 function is known to result in the accumulation of GlucA in rodents. Here, we aimed to explore genetic defects in AKR1A1 in patients with schizophrenia, which may result in the accumulation of GlucA. We identified 28 variants of AKR1A1 in patients with schizophrenia and control subjects. In particular, we identified a silent c.753G > A (rs745484618, p. Arg251Arg) variant located at the first position of exon 8 to be associated with schizophrenia. Using a minigene assay, we found that the c.753G > A variant induced exon 8 skipping in AKR1A1, resulting in a frameshift mutation, which in turn led to truncation of the AKR1A1 protein. Using the recombinant protein, we demonstrated that the truncated AKR1A1 completely lost its activity. Furthermore, we showed that AKR1A1 mRNA expression in the whole blood cells of individuals with the c.753G > A variant tended to be lower than that in those without the variants, leading to lower AKR activity. Our findings suggest that AKR1A1 carrying the c.753G > A variant induces exon skipping, leading to a loss of gene expression and enzymatic activity. Thus, GlucA patients with schizophrenia with the c.753G > A variant may show higher GlucA levels, leading to drug-resistant schizophrenia, since drug excretion by GlucA is enhanced.

U5 snRNA Interactions With Exons Ensure Splicing Precision

Imperfect conservation of human pre-mRNA splice sites is necessary to produce alternative isoforms. This flexibility is combined with the precision of the message reading frame. Apart from intron-termini GU_AG and the branchpoint A, the most conserved are the exon-end guanine and +5G of the intron start. Association between these guanines cannot be explained solely by base-pairing with U1 snRNA in the early spliceosome complex. U6 succeeds U1 and pairs +5G in the pre-catalytic spliceosome, while U5 binds the exon end. Current U5 snRNA reconstructions by CryoEM cannot explain the conservation of the exon-end G. Conversely, human mutation analyses show that guanines of both exon termini can suppress splicing mutations. Our U5 hypothesis explains the mechanism of splicing precision and the role of these conserved guanines in the pre-catalytic spliceosome. We propose: (1) optimal binding register for human exons and U5-the exon junction positioned at U5Loop1 C₃₉|C₃₈; (2) common mechanism for base-pairing of human U5 snRNA with diverse exons and bacterial Ll.LtrB intron with new loci in retrotransposition-guided by base pair geometry; and (3) U5 plays a significant role in specific exon recognition in the pre-catalytic spliceosome. Statistical analyses showed increased U5 Watson-Crick pairs with the 5'exon in the absence of +5G at the intron start. In 5'exon positions -3 and -5, this effect is specific to U5 snRNA rather than U1 snRNA of the early spliceosome. Increased U5 Watson-Crick pairs with 3'exon position +1 coincide with substitutions of the conserved -3C at the intron 3'end. Based on mutation and X-ray evidence, we propose that -3C pairs with U2 G₃₁ juxtaposing the branchpoint and the 3'intron end. The intron-termini pair, formed in the pre-catalytic spliceosome to be ready for transition after branching, and the early involvement of the 3'intron end ensure that the 3'exon contacts U5 in the pre-catalytic complex. We suggest that splicing precision is safeguarded cooperatively by U5, U6, and U2 snRNAs that stabilize the pre-catalytic complex by Watson-Crick base pairing. In addition, our new U5 model explains the splicing effect of exon-start +1G mutations: U5 Watson-Crick pairs with exon +2C/+3G strongly promote exon inclusion. We discuss potential applications for snRNA therapeutics and gene repair by reverse splicing.

Molecular analysis of GALT gene in Argentinian population: Correlation with enzyme activity and characterization of a novel Duarte-like allele

Background

Classical galactosemia is an autosomal recessive inherited metabolic disorder caused by mutations in the galactose-1-phosphate uridyltransferase (GALT) gene. GALT enzyme deficiency leads to the accumulation of galactose-1-phosphate in various organs, causing hepatic, renal and cerebral impairment. Over 300 mutations have been reported in the GALT gene. The aim of this study was to describe molecular characterization of GALT gene in Argentinian patients with decreased GALT activity, and to correlate molecular results with enzyme activity.

Methods

37 patients with enzyme activity below 6.3 μmol/h/g Hb (35% of normal value) were included. GALT activity was measured on red blood cells. DNA was extracted from peripheral blood. p.Gln188Arg mutation was studied by PCR-RFLP and, on samples negative or heterozygous, GALT gene was sequenced. In vivo splicing analysis of the GALT gene was performed on RNA extracted from leukocytes of one patient.

Results

14 different sequence variations were identified among 72 unrelated alleles. The two most common disease-causing mutations were p.Gln188Arg (24/72) and p.Lys285Asn (9/72). Three novel mutations were detected. One of them, c.688G>A, caused partial skipping of exon 9 of the GALT gene. Enzyme activity correlated with GALT genotype in 36 of the 37 patients.

Conclusion

This is the first report of sequence variations in the GALT gene in the Argentinian population. This study highlights the contribution of the molecular analysis to the diagnosis of Galactosemia and reveals c.688G>A as a novel Duarte-like variant, with a high prevalence in our population.

Prevalence of germline variants in consensus moderate-to-high-risk predisposition genes to hereditary breast and ovarian cancer in BRCA1/2-negative Brazilian patients

PURPOSE

This study aimed to identify and classify genetic variants in consensus moderate-to-high-risk predisposition genes associated with Hereditary Breast and Ovarian Cancer Syndrome (HBOC), in BRCA1/2-negative patients from Brazil.

METHODS

The study comprised 126 index patients who met NCCN clinical criteria and tested negative for all coding exons and intronic flanking regions of BRCA1/2 genes. Multiplex PCR-based assays were designed to cover the complete coding regions and flanking splicing sites of six genes implicated in HBOC. Sequencing was performed on HiSeq2500 Genome Analyzer.

RESULTS

Overall, we identified 488 unique variants. We identified five patients (3.97%) that harbored pathogenic or likely pathogenic variants in four genes: ATM (1), CHEK2 (2), PALB2 (1), and TP53 (1). One hundred and thirty variants were classified as variants of uncertain significance (VUS), 10 of which were predicted to disrupt mRNA splicing (seven non-coding variants and three coding variants), while other six missense VUS were classified as probably damaging by prediction algorithms.

CONCLUSION

A detailed mutational profile of non-BRCA genes is still being described in Brazil. In this study, we contributed to filling this gap, by providing important data on the diversity of genetic variants in a Brazilian high-risk patient cohort. ATM, CHEK2, PALB2 and TP53 are well established as HBOC predisposition genes, and the identification of deleterious variants in such actionable genes contributes to clinical management of probands and relatives.

Splicing Enhancers at Intron-Exon Borders Participate in Acceptor Splice Sites Recognition

Acceptor splice site recognition (3′ splice site: 3′ss) is a fundamental step in precursor messenger RNA (pre-mRNA) splicing. Generally, the U2 small nuclear ribonucleoprotein (snRNP) auxiliary factor (U2AF) heterodimer recognizes the 3′ss, of which U2AF35 has a dual function: (i) It binds to the intron–exon border of some 3′ss and (ii) mediates enhancer-binding splicing activators’ interactions with the spliceosome. Alternative mechanisms for 3′ss recognition have been suggested, yet they are still not thoroughly understood. Here, we analyzed 3′ss recognition where the intron–exon border is bound by a ubiquitous splicing regulator SRSF1. Using the minigene analysis of two model exons and their mutants, BRCA2 exon 12 and VARS2 exon 17, we showed that the exon inclusion correlated much better with the predicted SRSF1 affinity than 3′ss quality, which were assessed using the Catalog of Inferred Sequence Binding Preferences of RNA binding proteins (CISBP-RNA) database and maximum entropy algorithm (MaxEnt) predictor and the U2AF35 consensus matrix, respectively. RNA affinity purification proved SRSF1 binding to the model 3′ss. On the other hand, knockdown experiments revealed that U2AF35 also plays a role in these exons’ inclusion. Most probably, both factors stochastically bind the 3′ss, supporting exon recognition, more apparently in VARS2 exon 17. Identifying splicing activators as 3′ss recognition factors is crucial for both a basic understanding of splicing regulation and human genetic diagnostics when assessing variants’ effects on splicing.

Impact of acceptor splice site NAGTAG motif on exon recognition

Pre-mRNA splicing is an essential step in gene expression, when introns are removed and exons joined by the complex of proteins called spliceosome. Correct splicing requires a precise exon/intron junction definition, which is determined by a consensual donor and acceptor splice site at the 5' and 3' end, respectively. An acceptor splice site (3'ss) consists of highly conserved AG nucleotides in positions E^-2 and E^-1. These nucleotides can appear in tandem, located 3 bp from each other. Then they are referred to as NAGNAG or tandem 3'ss, which can be alternatively spliced. NAG/TAG 3'ss motif abundance is extremely low and cannot be easily explained by just a nucleotide preference in this position. We tested artificial NAG/TAG motif's potential negative effect on exon recognition using a minigene assay. Introducing the NAG/TAG motif into seven different exons revealed no general negative effect on exon recognition. The only observed effect was the partial use of the newly formed distal 3'ss. We can conclude that this motif's extremely low preference in a natural 3'ss is not a consequence of the NAG/TAG motif's negative effect on exon recognition, but more likely the result of other RNA processing aspects, such as an alternative 3'ss choice, decreased 3'ss strength, or incorporating an amber stop codon.

Integrative Profiling of Alternative Splicing Induced by U2AF1 S34F Mutation in Lung Adenocarcinoma Reveals a Mechanistic Link to Mitotic Stress

Mutations in spliceosome components have been implicated in carcinogenesis of various types of cancer. One of the most frequently found is U2AF1 S34F missense mutation. Functional analyses of this mutation have been largely limited to hematological malignancies although the mutation is also frequently seen in other cancer types including lung adenocarcinoma (LUAD). We examined the impact of knockdown (KD) of wild type (wt) U2AF1 and ectopic expression of two splice variant S34F mutant proteins in terms of alternative splicing (AS) pattern and cell cycle progression in A549 lung cancer cells. We demonstrate that induction of distinct AS events and disruption of mitosis at distinct sub-stages result from KD and ectopic expression of the mutant proteins. Importantly, when compared with the splicing pattern seen in LUAD patients with U2AF1 S34F mutation, ectopic expression of S34F mutants but not KD was shown to result in common AS events in several genes involved in cell cycle progression. Our study thus points to an active role of U2AF1 S34F mutant protein in inducing cell cycle dysregulation and mitotic stress. In addition, alternatively spliced genes which we describe here may represent novel potential markers of lung cancer development.

Linkage Study Revealed Complex Haplotypes in a Multifamily due to Different Mutations in CAPN3 Gene in an Iranian Ethnic Group

Calpainopathy is an autosomal recessive form of limb girdle muscular dystrophies which is caused by mutation in CAPN3 gene. In the present study, co-segregation of this disorder was analyzed with four short tandem repeat markers linked to the CAPN3 gene. Three apparently unrelated Iranian families with same ethnicity were investigated. Haplotype analysis and sequencing of the CAPN3 gene were performed. DNA sample from one of the patients was simultaneously sent for next-generation sequencing. DNA sequencing identified two mutations. It was seen as a homozygous c.2105C>T in exon 19 in one family, a homozygous novel mutation c.380G>A in exon 3 in another family, and a compound heterozygote form of these two mutations in the third family. Next-generation sequencing also confirmed our results. It was expected that, due to the rare nature of limb girdle muscular dystrophies, affected individuals from the same ethnic group share similar mutations. Haplotype analysis showed two different homozygote patterns in two families, yet a compound heterozygote pattern in the third family as seen in the mutation analysis. This study shows that haplotype analysis would help in determining presence of different founders.

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.

The Activation-Induced Assembly of an RNA/Protein Interactome Centered on the Splicing Factor U2AF2 Regulates Gene Expression in Human CD4 T Cells

Activation of CD4 T cells is a reaction to challenges such as microbial pathogens, cancer and toxins that defines adaptive immune responses. The roles of T cell receptor crosslinking, intracellular signaling, and transcription factor activation are well described, but the importance of post-transcriptional regulation by RNA-binding proteins (RBPs) has not been considered in depth. We describe a new model expanding and activating primary human CD4 T cells and applied this to characterizing activation-induced assembly of splicing factors centered on U2AF2. We immunoprecipitated U2AF2 to identify what mRNA transcripts were bound as a function of activation by TCR crosslinking and costimulation. In parallel, mass spectrometry revealed the proteins incorporated into the U2AF2-centered RNA/protein interactome. Molecules that retained interaction with the U2AF2 complex after RNAse treatment were designated as "central" interactome members (CIMs). Mass spectrometry also identified a second class of activation-induced proteins, "peripheral" interactome members (PIMs), that bound to the same transcripts but were not in physical association with U2AF2 or its partners. siRNA knockdown of two CIMs and two PIMs caused changes in activation marker expression, cytokine secretion, and gene expression that were unique to each protein and mapped to pathways associated with key aspects of T cell activation. While knocking down the PIM, SYNCRIP, impacts a limited but immunologically important set of U2AF2-bound transcripts, knockdown of U2AF1 significantly impairs assembly of the majority of protein and mRNA components in the activation-induced interactome. These results demonstrated that CIMs and PIMs, either directly or indirectly through RNA, assembled into activation-induced U2AF2 complexes and play roles in post-transcriptional regulation of genes related to cytokine secretion. These data suggest an additional layer of regulation mediated by the activation-induced assembly of RNA splicing interactomes that is important for understanding T cell activation.

Aberrant splicing of U12-type introns is the hallmark of ZRSR2 mutant myelodysplastic syndrome

Somatic mutations in the spliceosome gene ZRSR2 — located on the X chromosome — are associated with myelodysplastic syndrome (MDS). ZRSR2 is involved in the recognition of 3΄ splice site during the early stages of spliceosome assembly; however, its precise role in RNA splicing has remained unclear. Here, we characterize ZRSR2 as an essential component of the minor spliceosome (U12-dependent) assembly. shRNA mediated knockdown of ZRSR2 leads to impaired splicing of the U12-type introns, and RNA-Sequencing of MDS bone marrow reveals that loss of ZRSR2 activity causes increased mis-splicing. These splicing defects involve retention of the U12-type introns while splicing of the U2-type introns remain mostly unaffected. ZRSR2 deficient cells also exhibit reduced proliferation potential and distinct alterations in myeloid and erythroid differentiation in vitro. These data identify a specific role for ZRSR2 in RNA splicing and highlight dysregulated splicing of U12-type introns as a characteristic feature of ZRSR2 mutations in MDS.

U2AF1 mutations alter sequence specificity of pre-mRNA binding and splicing

We previously identified missense mutations in the U2AF1 splicing factor affecting codons S34 (S34F and S34Y) or Q157 (Q157R and Q157P) in 11% of patients with de novo myelodysplastic syndromes (MDS). Although the role of U2AF1 as an accessory factor in the U2 snRNP is well established, it is not yet clear how mutations affect splicing or contribute to MDS pathophysiology. We analyzed splice junctions in RNA-seq data generated from transfected CD34+ hematopoietic cells and found significant differences in the abundance of known and novel junctions in samples expressing mutant U2AF1 (S34F). For selected transcripts, splicing alterations detected by RNA-seq were confirmed by analysis of primary de novo MDS patient samples. These effects were not due to impaired U2AF1 (S34F) localization as it co-localized normally with U2AF2 within nuclear speckles. We further found evidence in the RNA-seq data for decreased affinity of U2AF1 (S34F) for uridine (relative to cytidine) at the e-3 position immediately upstream of the splice acceptor site and corroborated this finding using affinity binding assays. These data suggest that the S34F mutation alters U2AF1 function in the context of specific RNA sequences, leading to aberrant alternative splicing of target genes, some of which may be relevant for MDS pathogenesis.

A novel functional role for MMSET in RNA processing based on the link between the REIIBP isoform and its interaction with the SMN complex

The chromosomal translocation t(4;14) deregulates MMSET (WHSC1/NSD2) expression and is a poor prognostic factor in multiple myeloma (MM). MMSET encodes two major protein isoforms. We have characterized the role of the shorter isoform (REIIBP) in myeloma cells and identified a clear and novel interaction of REIIBP with members of the SMN (survival of motor neuron) complex that directly affects the assembly of the spliceosomal ribonucleic particles. Using RNA-seq we show that REIIBP influences the RNA splicing pattern of the cell. This new discovery provides novel insights into the understanding of MM pathology, and potential new leads for therapeutic targeting.

Exon first nucleotide mutations in splicing: evaluation of in silico prediction tools

Mutations in the first nucleotide of exons (E⁺¹) mostly affect pre-mRNA splicing when found in AG-dependent 3′ splice sites, whereas AG-independent splice sites are more resistant. The AG-dependency, however, may be difficult to assess just from primary sequence data as it depends on the quality of the polypyrimidine tract. For this reason, in silico prediction tools are commonly used to score 3′ splice sites. In this study, we have assessed the ability of sequence features and in silico prediction tools to discriminate between the splicing-affecting and non-affecting E⁺¹ variants. For this purpose, we newly tested 16 substitutions in vitro and derived other variants from literature. Surprisingly, we found that in the presence of the substituting nucleotide, the quality of the polypyrimidine tract alone was not conclusive about its splicing fate. Rather, it was the identity of the substituting nucleotide that markedly influenced it. Among the computational tools tested, the best performance was achieved using the Maximum Entropy Model and Position-Specific Scoring Matrix. As a result of this study, we have now established preliminary discriminative cut-off values showing sensitivity up to 95% and specificity up to 90%. This is expected to improve our ability to detect splicing-affecting variants in a clinical genetic setting.

Identification of six novel PTH1R mutations in families with a history of primary failure of tooth eruption

Primary Failure of tooth Eruption (PFE) is a non-syndromic disorder which can be caused by mutations in the parathyroid hormone receptor 1 gene (PTH1R). Traditionally, the disorder has been identified clinically based on post-emergent failure of eruption of permanent molars. However, patients with PTH1R mutations will not benefit from surgical and/or orthodontic treatment and it is therefore clinically important to establish whether a given failure of tooth eruption is caused by a PTH1R defect or not. We analyzed the PTH1R gene in six patients clinically diagnosed with PFE, all of which had undergone surgical and/or orthodontic interventions, and identified novel PTH1R mutations in all. Four of the six mutations were predicted to abolish correct mRNA maturation either through introduction of premature stop codons (c.947C>A and c.1082G>A), or by altering correct mRNA splicing (c.544-26_544-23del and c.989G>T). The latter was validated by transfection of minigenes. The six novel mutations expand the mutation spectrum for PFE from eight to 14 pathogenic mutations. Loss-of-function mutations in PTH1R are also associated with recessively inherited Blomstrand chondrodysplasia. We compiled all published PTH1R mutations and identified a mutational overlap between Blomstrand chondrodysplasia and PFE. The results suggest that a genetic approach to preclinical diagnosis will have important implication for surgical and orthodontic treatment of patients with failure of tooth eruption.

Patterns of missplicing due to somatic U2AF1 mutations in myeloid neoplasms

Recently, recurrent mutations of spliceosomal genes were frequently identified in myeloid malignancies, as well as other types of cancers. One of these spliceosomal genes, U2AF1, was affected by canonical somatic mutations in aggressive type of myeloid malignancies. We hypothesized that U2AF1 mutations causes defects of splicing (missplicing) in specific genes and that such misspliced genes might be important in leukemogenesis. We analyzed RNA deep sequencing to compare splicing patterns of 201 837 exons between the cases with U2AF1 mutations (n = 6) and wild type (n = 14). We identified different alternative splicing patterns in 35 genes comparing cells with mutant and wild-type U2AF1. U2AF1 mutations are associated with abnormal splicing of genes involved in functionally important pathways, such as cell cycle progression and RNA processing. In addition, many of these genes are somatically mutated or deleted in various cancers. Of note is that the alternative splicing patterns associated with U2AF1 mutations were associated with specific sequence signals at the affected splice sites. These novel observations support the hypothesis that U2AF1 mutations play a significant role in myeloid leukemogenesis due to selective missplicing of tumor-associated genes.

Defects in spliceosomal machinery: a new pathway of leukaemogenesis

Proper splicing of pre-mRNA is required for protein synthesis and therefore is a fundamental cellular function. The discovery of a variety of somatic spliceosomal mutations in haematological malignancies, including myeloid neoplasms and chronic lymphocytic leukaemia has pointed to a new leukaemogenic pathway involving spliceosomal dysfunction. Theoretically, spliceosomal mutations can lead to activation of incorrect splice sites, intron retention or aberrant alternative splicing occurring in patterns generated by mutations of individual spliceosomal proteins. Such events can produce a defective balance between protein isoforms leading to functional consequences including defective regulation of proliferation and differentiation. The observed pattern of occurrence of highly specific missense mutations, coupled with the lack of nonsense mutations and deletions, implies a gain-of-function or better gain-of-dysfunction mechanism. Incorrect splicing of downstream genes, such as tumour suppressor genes, may result in haploinsufficient expression through nonsense-mediated mRNA decay. Thus, spliceosomal mutations may, depending on the pattern of affected proteins, lead to similar functional effects on tumour suppressor genes as chromosomal deletions, epigenetic silencing or inactivating/hypomorphic mutations. The prognostic value of the most common mutations and their phenotypic association in the clinical setting is currently under investigation. It is likely that spliceosomal mutations may indicate sensitivity to spliceosome inhibitors applied in the form of a synthetic lethal approach. This review discusses the most current aspects of spliceosomal research in the context of haematological malignancies.

Mutations in the spliceosome machinery, a novel and ubiquitous pathway in leukemogenesis

Myelodysplastic syndromes (MDSs) are chronic and often progressive myeloid neoplasms associated with remarkable heterogeneity in the histomorphology and clinical course. Various somatic mutations are involved in the pathogenesis of MDS. Recently, mutations in a gene encoding a spliceosomal protein, SF3B1, were discovered in a distinct form of MDS with ring sideroblasts. Whole exome sequencing of 15 patients with myeloid neoplasms was performed, and somatic mutations in spliceosomal genes were identified. Sanger sequencing of 310 patients was performed to assess phenotype/genotype associations. To determine the functional effect of spliceosomal mutations, we evaluated pre-mRNA splicing profiles by RNA deep sequencing. We identified additional somatic mutations in spliceosomal genes, including SF3B1, U2AF1, and SRSF2. These mutations alter pre-mRNA splicing patterns. SF3B1 mutations are prevalent in low-risk MDS with ring sideroblasts, whereas U2AF1 and SRSF2 mutations are frequent in chronic myelomonocytic leukemia and advanced forms of MDS. SF3B1 mutations are associated with a favorable prognosis, whereas U2AF1 and SRSF2 mutations are predictive for shorter survival. Mutations affecting spliceosomal genes that result in defective splicing are a new leukemogenic pathway. Spliceosomal genes are probably tumor suppressors, and their mutations may constitute diagnostic biomarkers that could potentially serve as therapeutic targets.

Recurrent mutations in the U2AF1 splicing factor in myelodysplastic syndromes

Myelodysplastic syndromes (MDS) are hematopoietic stem cell disorders that often progress to chemotherapy-resistant secondary acute myeloid leukemia (sAML). We used whole-genome sequencing to perform an unbiased comprehensive screen to discover the somatic mutations in a sample from an individual with sAML and genotyped the loci containing these mutations in the matched MDS sample. Here we show that a missense mutation affecting the serine at codon 34 (Ser34) in U2AF1 was recurrently present in 13 out of 150 (8.7%) subjects with de novo MDS, and we found suggestive evidence of an increased risk of progression to sAML associated with this mutation. U2AF1 is a U2 auxiliary factor protein that recognizes the AG splice acceptor dinucleotide at the 3' end of introns, and the alterations in U2AF1 are located in highly conserved zinc fingers of this protein. Mutant U2AF1 promotes enhanced splicing and exon skipping in reporter assays in vitro. This previously unidentified, recurrent mutation in U2AF1 implicates altered pre-mRNA splicing as a potential mechanism for MDS pathogenesis.

Genetic variation of pre-mRNA alternative splicing in human populations

The precise splicing outcome of a transcribed gene is controlled by complex interactions between cis regulatory splicing signals and trans-acting regulators. In higher eukaryotes, alternative splicing is a prevalent mechanism for generating transcriptome and proteome diversity. Alternative splicing can modulate gene function, affect organismal phenotype and cause disease. Common genetic variation that affects splicing regulation can lead to differences in alternative splicing between human individuals and consequently impact expression level or protein function. In several well-documented examples, such natural variation of alternative splicing has indeed been shown to influence disease susceptibility and drug response. With new microarray and sequencing-based genomic technologies that can analyze eukaryotic transcriptomes at the exon or nucleotide level, it has become possible to globally compare the alternative splicing profiles across human individuals in any tissue or cell type of interest. Recent large-scale transcriptome studies using high-density splicing-sensitive microarray and deep RNA sequencing (RNA-Seq) have revealed widespread genetic variation of alternative splicing in humans. In the future, an extensive catalog of alternative splicing variation in human populations will help elucidate the molecular underpinnings of complex traits and human diseases, and shed light on the mechanisms of splicing regulation in human cells.