Solution structures of the SURP domains and the subunit-assembly mechanism within the splicing factor SF3a complex in 17S U2 snRNP

Biallelic variants in SNUPN cause a limb girdle muscular dystrophy with myofibrillar-like features

Alterations in RNA-splicing are a molecular hallmark of several neurological diseases, including muscular dystrophies, where mutations in genes involved in RNA metabolism or characterized by alterations in RNA splicing have been described. Here, we present five patients from two unrelated families with a limb-girdle muscular dystrophy (LGMD) phenotype carrying a biallelic variant in SNUPN gene. Snurportin-1, the protein encoded by SNUPN, plays an important role in the nuclear transport of small nuclear ribonucleoproteins (snRNPs), essential components of the spliceosome. We combine deep phenotyping, including clinical features, histopathology and muscle MRI, with functional studies in patient-derived cells and muscle biopsies to demonstrate that variants in SNUPN are the cause of a new type of LGMD according to current definition. Moreover, an in vivo model in Drosophila melanogaster further supports the relevance of Snurportin-1 in muscle. SNUPN patients show a similar phenotype characterized by proximal weakness starting in childhood, restrictive respiratory dysfunction and prominent contractures, although inter-individual variability in terms of severity even in individuals from the same family was found. Muscle biopsy showed myofibrillar-like features consisting of myotilin deposits and Z-disc disorganization. MRI showed predominant impairment of paravertebral, vasti, sartorius, gracilis, peroneal and medial gastrocnemius muscles. Conservation and structural analyses of Snurportin-1 p.Ile309Ser variant suggest an effect in nuclear-cytosol snRNP trafficking. In patient-derived fibroblasts and muscle, cytoplasmic accumulation of snRNP components is observed, while total expression of Snurportin-1 and snRNPs remains unchanged, which demonstrates a functional impact of SNUPN variant in snRNP metabolism. Furthermore, RNA-splicing analysis in patients' muscle showed widespread splicing deregulation, in particular in genes relevant for muscle development and splicing factors that participate in the early steps of spliceosome assembly. In conclusion, we report that SNUPN variants are a new cause of limb girdle muscular dystrophy with specific clinical, histopathological and imaging features, supporting SNUPN as a new gene to be included in genetic testing of myopathies. These results further support the relevance of splicing-related proteins in muscle disorders.

In silico identification, characterization, and expression analysis of RNA recognition motif (RRM) containing RNA-binding proteins in Aedes aegypti

RNA-binding proteins (RBPs) are the proteins that bind RNAs and regulate their functioning. RBPs in mosquitoes are gaining attention due to their ability to bind flaviviruses and regulate their replication and transmission. Despite their relevance, RBPs in mosquitoes are not explored much. In this study, we screened the whole genome of Aedes aegypti, the primary vector of several pathogenic viruses, and identified the proteins containing RNA recognition motif (RRM), the most abundant protein domain in eukaryotes. Using several in silico strategies, a total of 135 RRM-containing RBPs were identified in Ae. aegypti. The proteins were characterized based on their available annotations and the sequence similarity with Drosophila melanogaster. Ae. aegypti RRM-containing RBPs included serine/arginine-rich (SR) proteins, polyadenylate-binding proteins (PABP), heteronuclear ribonucleoproteins (hnRNP), small nuclear ribonucleoproteins (snRNP), splicing factors, eukaryotic initiation factors, transformers, and nucleolysins. Phylogenetic analysis revealed that the proteins and the domain organization are conserved among Ae. aegypti, Bombyx mori, and Drosophila melanogaster. However, the gene length and the intron-exon organization varied across the insect species. Expression analysis of the genes encoding RBPs using publicly available RNA sequencing data for different developmental time points of the mosquito life cycle starting from the ovary and eggs up to the adults revealed stage-specific expression with several genes preferentially expressed in early embryonic stages and blood-fed female ovaries. This is the first database for the Ae. aegypti RBPs that can serve as the reference base for future investigations. Stage-specific genes can be further explored to determine their role in mosquito growth and development with a focus on developing novel mosquito control strategies.

Structural basis for the interaction between the first SURP domain of the SF3A1 subunit in U2 snRNP and the human splicing factor SF1

SURP domains are exclusively found in splicing-related proteins in all eukaryotes. SF3A1, a component of the U2 snRNP, has two tandem SURP domains, SURP1, and SURP2. SURP2 is permanently associated with a specific short region of SF3A3 within the SF3A protein complex whereas, SURP1 binds to the splicing factor SF1 for recruitment of U2 snRNP to the early spliceosomal complex, from which SF1 is dissociated during complex conversion. Here, we determined the solution structure of the complex of SURP1 and the human SF1 fragment using nuclear magnetic resonance (NMR) methods. SURP1 adopts the canonical topology of α1-α2-310 -α3, in which α1 and α2 are connected by a single glycine residue in a particular backbone conformation, allowing the two α-helices to be fixed at an acute angle. A hydrophobic patch, which is part of the characteristic surface formed by α1 and α2, specifically contacts a hydrophobic cluster on a 16-residue α-helix of the SF1 fragment. Furthermore, whereas only hydrophobic interactions occurred between SURP2 and the SF3A3 fragment, several salt bridges and hydrogen bonds were found between the residues of SURP1 and the SF1 fragment. This finding was confirmed through mutational studies using bio-layer interferometry. The study also revealed that the dissociation constant between SURP1 and the SF1 fragment peptide was approximately 20 μM, indicating a weak or transient interaction. Collectively, these results indicate that the interplay between U2 snRNP and SF1 involves a transient interaction of SURP1, and this transient interaction appears to be common to most SURP domains, except for SURP2.

Chromatin-Associated Protein Sugp2 Involved in mRNA Alternative Splicing During Mouse Spermatogenesis

Mammalian spermatogenesis is a highly ordered process that is determined by chromatin-associated moderators which still remain poorly understood. Through a multi-control group proteomics strategy, we confirmed that Sugp2 was a chromatin-associated candidate protein, and its signal arose along spermatogenesis. The expression results showed that Sugp2, which is mainly expressed in the testis, had two transcripts, encoding one protein. During spermatogenesis, Sugp2 was enriched in the nucleus of male germ cells. With the depletion of Sugp2 by CRISPER-Cas9 technology, we found that Sugp2 controlled a network of genes on metal ion and ATP binding, suggesting that alternative splicing regulation by Sugp2 is involved in cellular ion and energy metabolism during spermatogenesis, while it had a little effect on meiotic progression and male fertility. Collectively, these data demonstrated that, as a chromatin-associated protein, Sugp2 mediated the alternative splicing regulatory network during spermatogenesis.

Disease-Causing Mutations in SF3B1 Alter Splicing by Disrupting Interaction with SUGP1

SF3B1, which encodes an essential spliceosomal protein, is frequently mutated in myelodysplastic syndromes (MDS) and many cancers. However, the defect of mutant SF3B1 is unknown. Here, we analyzed RNA sequencing data from MDS patients and confirmed that SF3B1 mutants use aberrant 3' splice sites. To elucidate the underlying mechanism, we purified complexes containing either wild-type or the hotspot K700E mutant SF3B1 and found that levels of a poorly studied spliceosomal protein, SUGP1, were reduced in mutant spliceosomes. Strikingly, SUGP1 knockdown completely recapitulated the splicing errors, whereas SUGP1 overexpression drove the protein, which our data suggest plays an important role in branchsite recognition, into the mutant spliceosome and partially rescued splicing. Other hotspot SF3B1 mutants showed similar altered splicing and diminished interaction with SUGP1. Our study demonstrates that SUGP1 loss is a common defect of spliceosomes with disease-causing SF3B1 mutations and, because this defect can be rescued, suggests possibilities for therapeutic intervention.

Structure and Conformational Dynamics of the Human Spliceosomal Bact Complex

The spliceosome is a highly dynamic macromolecular complex that precisely excises introns from pre-mRNA. Here we report the cryo-EM 3D structure of the human B^act spliceosome at 3.4 Å resolution. In the B^act state, the spliceosome is activated but not catalytically primed, so that it is functionally blocked prior to the first catalytic step of splicing. The spliceosomal core is similar to the yeast B^act spliceosome; important differences include the presence of the RNA helicase aquarius and peptidyl prolyl isomerases. To examine the overall dynamic behavior of the purified spliceosome, we developed a principal component analysis-based approach. Calculating the energy landscape revealed eight major conformational states, which we refined to higher resolution. Conformational differences of the highly flexible structural components between these eight states reveal how spliceosomal components contribute to the assembly of the spliceosome, allowing it to generate a dynamic interaction network required for its subsequent catalytic activation.

Transcriptome and long noncoding RNA sequencing of three extracellular vesicle subtypes released from the human colon cancer LIM1863 cell line

Previously we reported that LIM1863 colorectal cancer (CRC) cells secrete three distinct extracellular vesicle subtypes – two subpopulations of exosomes (apical EpCAM-Exos and basolateral A33-Exos) and shed microvesicles (sMVs) – with distinct protein and miRNA signatures. Here, we extend our omics approach to understand the fundamental role of LIM1863-derived EVs by performing a comprehensive analysis of their mRNAs and long non-coding RNAs (lncRNAs) using RNA-Seq. We show that 2,389 mRNAs, 317 pseudogene transcripts, 1,028 lncRNAs and 206 short non-coding RNAs selectively distributed to (i.e., are enriched in) LIM1863 EVs, relative to the parent cell. An Ensembl/UniProtKB analysis revealed 1,937 mRNAs encode canonical proteins, 348 isoforms (including splice-variant proteins), and 119 ‘missing proteins’ (i.e., annotated in Ensembl but not UniProtKB). Further dissection of our protein/RNA data revealed that 6/151 observed RNA binding proteins have the potential to interact with ~75% of EV-enriched RNAs. Intriguingly, the co-existence of U1 and U2 ribonucleoproteins and their cognate snRNAs in LIM1863 EVs suggests a possible association of CRC EVs with recipient cell splicing events. Our data reveal several potential lncRNA CRC biomarkers and novel splicing/fusion genes that, collectively, will advance our understanding of EV biology in CRC and accelerate the development of EV-based diagnostics and therapeutics.

Solution structure of the first RNA recognition motif domain of human spliceosomal protein SF3b49 and its mode of interaction with a SF3b145 fragment

The spliceosomal protein SF3b49, a component of the splicing factor 3b (SF3b) protein complex in the U2 small nuclear ribonucleoprotein, contains two RNA recognition motif (RRM) domains. In yeast, the first RRM domain (RRM1) of Hsh49 protein (yeast orthologue of human SF3b49) reportedly interacts with another component, Cus1 protein (orthologue of human SF3b145). Here, we solved the solution structure of the RRM1 of human SF3b49 and examined its mode of interaction with a fragment of human SF3b145 using NMR methods. Chemical shift mapping showed that the SF3b145 fragment spanning residues 598–631 interacts with SF3b49 RRM1, which adopts a canonical RRM fold with a topology of β1‐α1‐β2‐β3‐α2‐β4. Furthermore, a docking model based on NOESY measurements suggests that residues 607–616 of the SF3b145 fragment adopt a helical structure that binds to RRM1 predominantly via α1, consequently exhibiting a helix–helix interaction in almost antiparallel. This mode of interaction was confirmed by a mutational analysis using GST pull‐down assays. Comparison with structures of all RRM domains when complexed with a peptide found that this helix–helix interaction is unique to SF3b49 RRM1. Additionally, all amino acid residues involved in the interaction are well conserved among eukaryotes, suggesting evolutionary conservation of this interaction mode between SF3b49 RRM1 and SF3b145.

PDB Code(s): 5GVQ

Cellular stress response 1 down-regulates the expression of epidermal growth factor receptor and platelet-derived growth factor receptor through inactivation of splicing factor 3A3

Cellular stress response 1 (CSR1) is a tumor suppressor gene that plays an important role in regulating cell death. In this report, we show that the N-terminus of CSR1 interacts with splicing factor 3A, subunit 3 (SF3A3). The SF3A3 binding motif was identified in the region of amino acids 62-91 of CSR1 through cell-free binding analyses. The interaction between CSR1 and SF3A3 led to migration of SF3A3 from nucleus to cytoplasm. The cytoplasmic redistribution of SF3A3 significantly reduced the splicing efficiency of epidermal growth factor receptor and platelet-derived growth factor receptor. Induction of CSR1 or down-regulation of SF3A3 also significantly reduced the splicing activity of oxytocin reporter gene both in vivo and in vitro. Mutant CSR1 that lacks the SF3A3 binding motif contained no RNA splicing regulatory activity, while the peptide corresponding to the SF3A3 binding motif in CSR1 interfered with the wild-type CSR1 mediated inhibition of RNA splicing. Interaction of CSR1 and SF3A3 is essential for CSR1 mediated cell death. To our knowledge, this is the first report demonstrating that RNA splicing is negatively regulated by redistribution of a splicing factor. © 2016 Wiley Periodicals, Inc.

Mammalian splicing factor SF1 interacts with SURP domains of U2 snRNP-associated proteins

Splicing factor 1 (SF1) recognizes the branch point sequence (BPS) at the 3′ splice site during the formation of early complex E, thereby pre-bulging the BPS adenosine, thought to facilitate subsequent base-pairing of the U2 snRNA with the BPS. The 65-kDa subunit of U2 snRNP auxiliary factor (U2AF65) interacts with SF1 and was shown to recruit the U2 snRNP to the spliceosome. Co-immunoprecipitation experiments of SF1-interacting proteins from HeLa cell extracts shown here are consistent with the presence of SF1 in early splicing complexes. Surprisingly almost all U2 snRNP proteins were found associated with SF1. Yeast two-hybrid screens identified two SURP domain-containing U2 snRNP proteins as partners of SF1. A short, evolutionarily conserved region of SF1 interacts with the SURP domains, stressing their role in protein–protein interactions. A reduction of A complex formation in SF1-depleted extracts could be rescued with recombinant SF1 containing the SURP-interaction domain, but only partial rescue was observed with SF1 lacking this sequence. Thus, SF1 can initially recruit the U2 snRNP to the spliceosome during E complex formation, whereas U2AF65 may stabilize the association of the U2 snRNP with the spliceosome at later times. In addition, these findings may have implications for alternative splicing decisions.

Proteins associated with SF3a60 in T. brucei

Trypanosoma brucei relies on Spliced leader trans splicing to generate functional messenger RNAs. Trans splicing joins the specialized SL exon from the SL RNA to pre-mRNAs and is mediated by the trans-spliceosome, which is made up of small nuclear ribonucleoprotein particles and non-snRNP factors. Although the trans spliceosome is essential for trypanosomatid gene expression, not all spliceosomal protein factors are known and of these, only a few are completely characterized. In this study, we have characterized the trypanosome Splicing Factor, SF3a60, the only currently annotated SF3a component. As expected, epitope-tagged SF3a60 localizes in the trypanosome nucleus. SF3a60 is essential for cell viability but its depletion seem to have no detectable effect on trans-splicing. In addition, we used SF3a60 as bait in a Yeast-2-hybrid system screen and identified its interacting protein factors. The interactions with SF3a120, SF3a66 and SAP130 were confirmed by tandem affinity purification and mass spectrometry.

Re-evaluation of the role of calcium homeostasis endoplasmic reticulum protein (CHERP) in cellular calcium signaling

Changes in cytoplasmic Ca(2+) concentration, resulting from activation of intracellular Ca(2+) channels within the endoplasmic reticulum, regulate several aspects of cellular growth and differentiation. Ca(2+) homeostasis endoplasmic reticulum protein (CHERP) is a ubiquitously expressed protein that has been proposed as a regulator of both major families of endoplasmic reticulum Ca(2+) channels, inositol 1,4,5-trisphosphate receptors (IP(3)Rs) and ryanodine receptors (RyRs), with resulting effects on mitotic cycling. However, the manner by which CHERP regulates intracellular Ca(2+) channels to impact cellular growth is unknown. Here, we challenge previous findings that CHERP acts as a direct cytoplasmic regulator of IP(3)Rs and RyRs and propose that CHERP acts in the nucleus to impact cellular proliferation by regulating the function of the U2 snRNA spliceosomal complex. The previously reported effects of CHERP on cellular growth therefore are likely indirect effects of altered spliceosomal function, consistent with prior data showing that loss of function of U2 snRNP components can interfere with cell growth and induce cell cycle arrest.

Intrinsic disorder in the human spliceosomal proteome

The spliceosome is a molecular machine that performs the excision of introns from eukaryotic pre-mRNAs. This macromolecular complex comprises in human cells five RNAs and over one hundred proteins. In recent years, many spliceosomal proteins have been found to exhibit intrinsic disorder, that is to lack stable native three-dimensional structure in solution. Building on the previous body of proteomic, structural and functional data, we have carried out a systematic bioinformatics analysis of intrinsic disorder in the proteome of the human spliceosome. We discovered that almost a half of the combined sequence of proteins abundant in the spliceosome is predicted to be intrinsically disordered, at least when the individual proteins are considered in isolation. The distribution of intrinsic order and disorder throughout the spliceosome is uneven, and is related to the various functions performed by the intrinsic disorder of the spliceosomal proteins in the complex. In particular, proteins involved in the secondary functions of the spliceosome, such as mRNA recognition, intron/exon definition and spliceosomal assembly and dynamics, are more disordered than proteins directly involved in assisting splicing catalysis. Conserved disordered regions in spliceosomal proteins are evolutionarily younger and less widespread than ordered domains of essential spliceosomal proteins at the core of the spliceosome, suggesting that disordered regions were added to a preexistent ordered functional core. Finally, the spliceosomal proteome contains a much higher amount of intrinsic disorder predicted to lack secondary structure than the proteome of the ribosome, another large RNP machine. This result agrees with the currently recognized different functions of proteins in these two complexes.

Structure and assembly of the SF3a splicing factor complex of U2 snRNP

SF3a is an evolutionarily conserved heterotrimeric complex essential for pre-mRNA splicing. It functions in spliceosome assembly within the mature U2 snRNP (small nuclear ribonucleoprotein particle), and its displacement from the spliceosome initiates the first step of the splicing reaction. We have identified a core domain of the yeast SF3a complex required for complex assembly and determined its crystal structure. The structure shows a bifurcated assembly of three subunits, Prp9, Prp11 and Prp21, with Prp9 interacting with Prp21 via a bidentate-binding mode, and Prp21 wrapping around Prp11. Structure-guided biochemical analysis also shows that Prp9 harbours a major binding site for stem-loop IIa of U2 snRNA. These findings provide mechanistic insights into the assembly of U2 snRNP.

Interaction domains and nuclear targeting signals in subunits of the U2 small nuclear ribonucleoprotein particle-associated splicing factor SF3a

Human splicing factor SF3a is a component of the mature U2 small nuclear ribonucleoprotein particle (snRNP) and its three subunits of 60, 66, and 120 kDa are essential for splicing in vitro and in vivo. The SF3a heterotrimer forms in the cytoplasm and enters the nucleus independently of the U2 snRNP. Here, we have analyzed domains required for in vitro interactions between the SF3a subunits. Our results indicate that the SF3a66-SF3a120 interaction is mediated by a 27-amino acid region in SF3a120 C-terminal to the second suppressor-of-white-apricot and prp21/spp91 domain and amino acids 108-210 of SF3a66. Neither of these sequences contains known structural motifs, suggesting that the interaction domains are novel. Moreover, an ∼100-amino acid region, including the SURP2 domain of SF3a120 but extending into neighboring regions, is sufficient for binding to SF3a60. Analysis of determinants for nuclear import of SF3a demonstrates that SF3a120 provides the major nuclear localization signal and SF3a60 contributes to nuclear import.

Detailed close-ups and the big picture of spliceosomes

The spliceosome is the huge macromolecular assembly responsible for the removal of introns from pre-mRNA transcripts. The size and complexity of this dynamic cellular machine dictate that structural analysis of the spliceosome is best served by a combination of techniques. Electron microscopy is providing a more global, albeit less detailed, view of spliceosome assemblies. X-ray crystallographers and NMR spectroscopists are steadily reporting more atomic resolution structures of individual spliceosome components and fragments. Increasingly, structures of these individual pieces in complex with binding partners are yielding insights into the interfaces that hold the entire spliceosome assembly together. Although the information arising from the various structural studies of splicing machinery has not yet fully converged into a complete model, we can expect that a detailed understanding of spliceosome structure will arise at the juncture of structural and computational modeling methods.

Solution structure of the second RNA recognition motif (RRM) domain of murine T cell intracellular antigen-1 (TIA-1) and its RNA recognition mode

T cell intracellular antigen-1 (TIA-1), an apoptosis promoting factor, functions as a splicing regulator for the Fas pre-mRNA. TIA-1 possesses three RNA recognition motifs (RRMs) and a glutamine-rich domain. The second RRM (RRM2) is necessary and sufficient for tight, sequence-specific binding to the uridine-rich sequences buried around the 5' splice sites. In the present study, we solved the solution structure of the murine TIA-1 RRM2 by heteronuclear-nuclear magnetic resonance spectroscopy. The TIA-1 RRM2 adopts the RRM fold (betaalphabetabetaalphabeta) and possesses an extra beta-strand between beta2 and beta3, which forms an additional beta-sheet with the C-terminal part of beta2. We refer to this structure as the beta2-beta2' beta-loop. Interestingly, this characteristic beta-loop structure is conserved among a number of RRMs, including the U2AF65 RRM2 and the Sex-lethal RRM1 and RRM2, which also bind to uridine-rich RNAs. Furthermore, we identified a new sequence motif in the beta2-beta2' beta-loop, the DxxT motif. Chemical shift perturbation analyses of both the main and side chains upon binding to the uridine pentamer RNA revealed that most of the beta-sheet surface, including the beta2-beta2' beta-loop, is involved in the RNA binding. An investigation of the chemical shift perturbation revealed similarity in the RNA recognition modes between the TIA-1 and U2AF65 RRMs.