Structural basis for polypyrimidine tract recognition by the essential pre-mRNA splicing factor U2AF65

From benign to pathogenic variants and vice versa: pyrimidine transitions at position -3 of TAG and CAG 3' splice sites

In the human genome, CAG 3' splice sites (3'ss) are more than twice as frequent as TAG 3'ss. The greater abundance of the former has been attributed to a higher probability of exon skipping upon cytosine-to-thymine transitions at intron position -3 (-3C > T) than thymine-to-cytosine variants (-3T > C). However, molecular mechanisms underlying this bias and its clinical impact are poorly understood. In this study, base-pairing probabilities (BPPs) and RNA secondary structures were compared between CAG 3'ss that produced more skipping of downstream exons than their mutated UAG versions (termed "laggard" CAG 3'ss) and UAG 3'ss that resulted in more skipping than their mutated CAG counterparts (canonical 3'ss). The laggard CAG 3'ss showed significantly higher BPPs across intron-exon boundaries than canonical 3'ss. The difference was centered on positions -5 to -1 relative to the intron-exon junction, the region previously shown to exhibit the strongest high-resolution ultraviolet crosslinking to the small subunit of auxiliary factor of U2 snRNP (U2AF1). RNA secondary structure predictions suggested that laggard CAG 3'ss were more often sequestered in paired conformations and in longer stem structures while canonical 3'ss were more frequently unpaired. Taken together, the excess of base-pairing at 3'ss has a potential to alter the hierarchy in intrinsic splicing efficiency of human YAG 3'ss from canonical CAG > UAG to non-canonical UAG > CAG, to modify the clinical impact of transitions at this position and to change their classification from pathogenic to benign or vice versa.

RNAelem: an algorithm for discovering sequence-structure motifs in RNA bound by RNA-binding proteins

Motivation

RNA-binding proteins (RBPs) play a crucial role in the post-transcriptional regulation of RNA. Given their importance, analyzing the specific RNA patterns recognized by RBPs has become a significant research focus in bioinformatics. Deep Neural Networks have enhanced the accuracy of prediction for RBP-binding sites, yet understanding the structural basis of RBP-binding specificity from these models is challenging due to their limited interpretability. To address this, we developed RNAelem, which combines profile context-free grammar and the Turner energy model for RNA secondary structure to predict sequence-structure motifs in RBP-binding regions.

Results

RNAelem exhibited superior detection accuracy compared to existing tools for RNA sequences with structural motifs. Upon applying RNAelem to the eCLIP database, we were not only able to reproduce many known primary sequence motifs in the absence of secondary structures, but also discovered many secondary structural motifs that contained sequence-nonspecific insertion regions. Furthermore, the high interpretability of RNAelem yielded insightful findings such as long-range base-pairing interactions in the binding region of the U2AF protein.

Availability and implementation

The code is available at https://github.com/iyak/RNAelem.

Inefficient recruitment of DDX39B impedes pre-spliceosome assembly on FOXP3 introns

Forkhead box P3 (FOXP3) is the master fate-determining transcription factor in regulatory T (Treg) cells and is essential for their development, function, and homeostasis. Mutations in FOXP3 cause immunodysregulation polyendocrinopathy enteropathy X-linked (IPEX) syndrome, and aberrant expression of FOXP3 has been implicated in other diseases such as multiple sclerosis and cancer. We previously demonstrated that pre-mRNA splicing of FOXP3 RNAs is highly sensitive to levels of DExD-box polypeptide 39B (DDX39B), and here we investigate the mechanism of this sensitivity. FOXP3 introns have cytidine (C)-rich/uridine (U)-poor polypyrimidine (py) tracts that are responsible for their inefficient splicing and confer sensitivity to DDX39B. We show that there is a deficiency in the assembly of commitment complexes (CCs) on FOXP3 introns, which is consistent with the lower affinity of U2AF2 for C-rich/U-poor py tracts. Our data indicate an even stronger effect on the conversion of CCs to pre-spliceosomes. We propose that this is due to an altered conformation that U2AF2 adopts when it binds to C-rich/U-poor py tracts and that this conformation has a lower affinity for DDX39B. As a consequence, CCs assembled on FOXP3 introns are defective in recruiting DDX39B, and this leads to the inefficient assembly of pre-spliceosome complexes.

Beyond reader proteins: RNA binding proteins and RNA modifications in conversation to regulate gene expression

Post-transcriptional mRNA modifications play diverse roles in gene expression and RNA function. In many cases, RNA modifications function by altering how cellular machinery such as RNA binding proteins (RBPs) interact with RNA substrates. For instance, N6-methyladenosine (m6A) is recognized by the well-characterized YTH domain-containing family of "reader" proteins. For other mRNA modifications, similar global readers of modification status have not been clearly defined. Rather, most interactions between RBPs and RNA modifications have a more complicated dependence on sequence context and binding modality. The current handful of studies that demonstrate modifications impacting protein binding likely represent only a fraction of the full landscape. In this review, we dissect the known instances of RNA modifications altering RBP binding, specifically m6A, N1-methyladenosine (m1A), 5-methylcytosine (m5C), pseudouridine (Ψ), and internal N7-methylguanosine. We then review the biochemical properties of these and other identified mRNA modifications including dihydrouridine (D), N4-acetylcytosine (ac4C), and 2'-O-Methylation (Nme). We focus on how these properties would be likely to impact RNA:RBP interactions, including by changes to hydrogen bond potential, base-stacking efficiency, and RNA conformational preferences. The effects of RNA modifications on secondary structure have been well-studied, and we briefly discuss how structural effects imparted by modifications can lead to protein binding changes. Finally, we discuss strategies for uncovering as-yet-to-be identified modification-sensitive RBP:RNA Interactions. Coordinating future efforts to intersect the epitranscriptome and the RNA-protein interactome will illuminate the rules governing RNA modification recognition and the mechanisms responsible for the biological consequences of mRNA modification. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Processing > RNA Editing and Modification.

Simultaneous Localization of Two High Affinity Divalent Metal Ion Binding Sites in the Tetracycline RNA Aptamer with Mn2+-Based Pulsed Dipolar EPR Spectroscopy

Mg2+ ions play an essential part in stabilizing the tertiary structure of nucleic acids. While the importance of these ions is well documented, their localization and elucidation of their role in the structure and dynamics of nucleic acids are often challenging. In this work, pulsed electron-electron double resonance spectroscopy (PELDOR, also known as DEER) was used to localize two high affinity divalent metal ion binding sites in the tetracycline RNA aptamer with high accuracy. For this purpose, the aptamer was labeled at different positions with a semirigid nitroxide spin label and diamagnetic Mg2+ was replaced with paramagnetic Mn2+, which did not alter the folding process or ligand binding. Out of the several divalent metal ion binding sites that are known from the crystal structure, two binding sites with high affinity were detected: one that is located at the ligand binding center and another at the J1/2 junction of the RNA.

Functional analysis of 3'-UTR hairpins supports a two-tiered model for posttranscriptional regulation of MAT2A by METTL16

S-adenosylmethionine (SAM) is the methyl donor for nearly all cellular methylation events, so cells need to carefully control SAM levels. MAT2A encodes the only SAM synthetase expressed in the majority of human cells, and its 3'-UTR has six conserved regulatory hairpins (hp1-6) that can be methylated by the N6-methyladenosine methyltransferase METTL16. Hp1 begins 8 nt from the stop codon, whereas hp2-6 are clustered further downstream (∼800 nt). These hairpins have been proposed to regulate MAT2A mRNA levels in response to intracellular SAM levels by regulating intron detention of the last intron of MAT2A and by modulating the stability of the fully spliced mRNA. However, a dissection of these two posttranscriptional mechanisms has not been previously reported. Using a modular reporter system, we show that hp1 functions primarily when the detained intron is included in the reporter and when that intron has a suboptimal polypyrimidine tract. In contrast, the hp2-6 cluster modulates mRNA stability independent of the detained intron, although hp1 may make a minor contribution to the regulation of decay as well. Taken with previously published reports, these data support a two-tiered model for MAT2A posttranscriptional regulation by METTL16 through its interactions with hp1 and hp2-6. In the upstream tier, hp1 and METTL16 control MAT2A intron detention, whereas the second tier involves METTL16-dependent methylation of hp2-6 to control MAT2A mRNA stability. Thus, cells use a similar set of molecular factors to achieve considerable complexity in the posttranscriptional regulation of SAM homeostasis.

Molecular basis of RNA-binding and autoregulation by the cancer-associated splicing factor RBM39

Pharmacologic depletion of RNA-binding motif 39 (RBM39) using aryl sulfonamides represents a promising anti-cancer therapy but requires high levels of the adaptor protein DCAF15. Consequently, novel approaches to deplete RBM39 in an DCAF15-independent manner are required. Here, we uncover that RBM39 autoregulates via the inclusion of a poison exon into its own pre-mRNA and identify the cis-acting elements that govern this regulation. We also determine the NMR solution structures of RBM39's tandem RNA recognition motifs (RRM1 and RRM2) bound to their respective RNA targets, revealing how RRM1 recognises RNA stem loops whereas RRM2 binds specifically to single-stranded N(G/U)NUUUG. Our results support a model where RRM2 selects the 3'-splice site of a poison exon and the RRM3 and RS domain stabilise the U2 snRNP at the branchpoint. Our work provides molecular insights into RBM39-dependent 3'-splice site selection and constitutes a solid basis to design alternative anti-cancer therapies.

Purifying selection against spurious splicing signals contributes to the base composition evolution of the polypyrimidine tract

Among eukaryotes, the major spliceosomal pathway is highly conserved. While long introns may contain additional regulatory sequences, the ones in short introns seem to be nearly exclusively related to splicing. Although these regulatory sequences involved in splicing are well-characterized, little is known about their evolution. At the 3' end of introns, the splice signal nearly universally contains the dimer AG, which consists of purines, and the polypyrimidine tract upstream of this 3' splice signal is characterized by over-representation of pyrimidines. If the over-representation of pyrimidines in the polypyrimidine tract is also due to avoidance of a premature splicing signal, we hypothesize that AG should be the most under-represented dimer. Through the use of DNA-strand asymmetry patterns, we confirm this prediction in fruit flies of the genus Drosophila and by comparing the asymmetry patterns to a presumably neutrally evolving region, we quantify the selection strength acting on each motif. Moreover, our inference and simulation method revealed that the best explanation for the base composition evolution of the polypyrimidine tract is the joint action of purifying selection against a spurious 3' splice signal and the selection for pyrimidines. Patterns of asymmetry in other eukaryotes indicate that avoidance of premature splicing similarly affects the nucleotide composition in their polypyrimidine tracts.

Early Splicing Complexes and Human Disease

Over the last decade, our understanding of spliceosome structure and function has significantly improved, refining the study of the impact of dysregulated splicing on human disease. As a result, targeted splicing therapeutics have been developed, treating various diseases including spinal muscular atrophy and Duchenne muscular dystrophy. These advancements are very promising and emphasize the critical role of proper splicing in maintaining human health. Herein, we provide an overview of the current information on the composition and assembly of early splicing complexes-commitment complex and pre-spliceosome-and their association with human disease.

Monovalent Ionic Atmosphere Modulates the Selection of Suboptimal RNA Sequences by Splicing Factors' RNA Recognition Motifs

The U2AF2 splicing factor is involved in the RNA recognition of the pre-mRNA poly-pyrimidine signaling sequence. This protein contains two RRM domains connected by a flexible linker, which ensure the preferential selection of a poly-uridine sequence over a poly-cytosine one. In this work, all-atom simulations provide insights into the U2AF2 recognition mechanism and on the features underlying its selectivity. Our outcomes show that U2AF2's RNA recognition is driven by cooperative events modulated by RNA-protein and RNA-ion interactions. Stunningly, monovalent ions contribute to mediating the binding of the weakly binding polyC strand, thus contributing to the selection of suboptimal poly-pyrimidine tracts. This finding broadens our understanding of the diverse traits tuning splicing factors' selectivity and adaptability to precisely handle and process diverse pre-mRNA sequences.

Reliability and accuracy of single-molecule FRET studies for characterization of structural dynamics and distances in proteins

Single-molecule Förster-resonance energy transfer (smFRET) experiments allow the study of biomolecular structure and dynamics in vitro and in vivo. We performed an international blind study involving 19 laboratories to assess the uncertainty of FRET experiments for proteins with respect to the measured FRET efficiency histograms, determination of distances, and the detection and quantification of structural dynamics. Using two protein systems with distinct conformational changes and dynamics, we obtained an uncertainty of the FRET efficiency ≤0.06, corresponding to an interdye distance precision of ≤2 Å and accuracy of ≤5 Å. We further discuss the limits for detecting fluctuations in this distance range and how to identify dye perturbations. Our work demonstrates the ability of smFRET experiments to simultaneously measure distances and avoid the averaging of conformational dynamics for realistic protein systems, highlighting its importance in the expanding toolbox of integrative structural biology.

A splicing silencer in SMN2 intron 6 is critical in spinal muscular atrophy

Spinal muscular atrophy (SMA) is a fatal neuromuscular disease caused by homozygous deletions or mutations of the SMN1 gene. SMN2 is a paralogous gene of SMN1 and a modifying gene of SMA. A better understanding of how SMN2 exon 7 splicing is regulated helps discover new therapeutic targets for SMA therapy. Based on an antisense walk method to map exonic and intronic splicing silencers (ESSs and ISSs) in SMN2 exon 7 and the proximal regions of its flanking introns, we identified one ISS (ISS6-KH) at upstream of the branch point site in intron 6. By using mutagenesis-coupled RT-PCR with SMN1/2 minigenes, immunochromatography, overexpression and siRNA-knockdown, we found this ISS consists of a bipartite hnRNP A1 binding cis-element and a poly-U sequence located between the proximal hnRNP A1 binding site (UAGCUA) and the branch site. Both HuR and hnRNP C1 proteins promote exon 7 skipping through the poly-U stretch. Mutations or deletions of these motifs lead to efficient SMN2 exon 7 inclusion comparable to SMN1 gene. Furthermore, we identified an optimal antisense oligonucleotide that binds the intron six ISS and causes striking exon 7 inclusion in the SMN2 gene in patient fibroblasts and SMA mouse model. Our findings demonstrate that this novel ISS plays an important role in SMN2 exon 7 skipping and highlight a new therapeutic target for SMA therapy.

All-Atom Simulations Elucidate the Impact of U2AF2 Cancer-Associated Mutations on Pre-mRNA Recognition

The U2AF2 splicing factor, made of two tandem RNA recognition motifs (RRMs) joined by a flexible linker, selects the intronic polypyrimidine sequence of premature mRNA, thus ensuring splicing fidelity. Increasing evidence links mutations of key splicing factors, including U2AF2, to a variety of cancers. Nevertheless, the impact of U2AF2 cancer-associated mutations on polypyrimidine recognition remains unclear. Here, we combined extensive (18 μs-long) all-atom molecular dynamics simulations and dynamical network theory analysis (NWA) of U2AF2, in its wild-type form and in the presence of the six most frequent cancer-associated mutations, bound to a poly-U strand. Our results reveal that the selected mutations affect the pre-mRNA binding at two hot spot regions, irrespectively of where these mutants are placed on the distinct U2AF2 domains. Complementarily, NWA traced the existence of cross-communication pathways, connecting each mutation site to these recognition hot spots, whose strength is altered by the mutations. Our outcomes suggest the existence of a structural/dynamical interplay of the two U2AF2's RRMs underlying the recognition of the polypyrimidine tract and reveal that the cancer-associated mutations affect the polypyrimidine selection by altering the RRMs' cooperativity. This mechanism may be shared by other RNA binding proteins hallmarked, like U2AF2, by multidomain architecture and high plasticity.

Long non-coding RNA and RNA-binding protein interactions in cancer: Experimental and machine learning approaches

Understanding the complex and specific roles played by non-coding RNAs (ncRNAs), which comprise the bulk of the genome, is important for understanding virtually every hallmark of cancer. This large group of molecules plays pivotal roles in key regulatory mechanisms in various cellular processes. Regulatory mechanisms, mediated by long non-coding RNA (lncRNA) and RNA-binding protein (RBP) interactions, are well documented in several types of cancer. Their effects are enabled through networks affecting lncRNA and RBP stability, RNA metabolism including N⁶-methyladenosine (m⁶A) and alternative splicing, subcellular localization, and numerous other mechanisms involved in cancer. In this review, we discuss the reciprocal interplay between lncRNAs and RBPs and their involvement in epigenetic regulation via histone modifications, as well as their key role in resistance to cancer therapy. Other aspects of RBPs including their structural domains, provide a deeper knowledge on how lncRNAs and RBPs interact and exert their biological functions. In addition, current state-of-the-art knowledge, facilitated by machine and deep learning approaches, unravels such interactions in better details to further enhance our understanding of the field, and the potential to harness RNA-based therapeutics as an alternative treatment modality for cancer are discussed.

The nucleolus is the site for inflammatory RNA decay during infection

Inflammatory cytokines are key signaling molecules that can promote an immune response, thus their RNA turnover must be tightly controlled during infection. Most studies investigate the RNA decay pathways in the cytosol or nucleoplasm but never focused on the nucleolus. Although this organelle has well-studied roles in ribosome biogenesis and cellular stress sensing, the mechanism of RNA decay within the nucleolus is not completely understood. Here, we report that the nucleolus is an essential site of inflammatory pre-mRNA instability during infection. RNA-sequencing analysis reveals that not only do inflammatory genes have higher intronic read densities compared with non-inflammatory genes, but their pre-mRNAs are highly enriched in nucleoli during infection. Notably, nucleolin (NCL) acts as a guide factor for recruiting cytosine or uracil (C/U)-rich sequence-containing inflammatory pre-mRNAs and the Rrp6-exosome complex to the nucleolus through a physical interaction, thereby enabling targeted RNA delivery to Rrp6-exosomes and subsequent degradation. Consequently, Ncl depletion causes aberrant hyperinflammation, resulting in a severe lethality in response to LPS. Importantly, the dynamics of NCL post-translational modifications determine its functional activity in phases of LPS. This process represents a nucleolus-dependent pathway for maintaining inflammatory gene expression integrity and immunological homeostasis during infection.

The nucleolus is the traditional site for ribosomal RNA biogenesis. Here, the authors find that the nucleolus is a site of inflammatory pre-mRNA turnover and elucidated how immune homeostasis can be maintained by controlling inflammatory gene expression.

Insight into the Tumor Suppression Mechanism from the Structure of Human Polypyrimidine Splicing Factor (PSF/SFPQ) Complexed with a 30mer RNA from Murine Virus-like 30S Transcript-1

Human polypyrimidine-binding splicing factor (PSF/SFPQ) is a tumor suppressor protein that regulates the gene expression of several proto-oncogenes and binds to the 5'-polyuridine negative-sense template (5'-PUN) of some RNA viruses. The activity of PSF is negatively regulated by long-noncoding RNAs, human metastasis associated in lung adenocarcinoma transcript-1 and murine virus-like 30S transcript-1 (VL30-1). PSF is a 707-amino acid protein that has a DNA-binding domain and two RNA recognition motifs (RRMs). Although the structure of the apo-truncated PSF is known, how PSF recognizes RNA remains elusive. Here, we report the 2.8 Å and 3.5 Å resolution crystal structures of a biologically active truncated construct of PSF (sPSF, consisting of residues 214-598) alone and in a complex with a 30mer fragment of VL30-1 RNA, respectively. The structure of the complex reveals how the 30mer RNA is recognized at two U-specific induced-fit binding pockets, located at the previously unrecognized domain-swapped, inter-subunit RRM1 (of the first subunit)-RRM2 (of the second subunit) interfaces that do not exist in the apo structure. Thus, the sPSF dimer appears to have two conformations in solution: one in a low-affinity state for RNA binding, as seen in the apo-structure, and the other in a high-affinity state for RNA binding, as seen in the sPSF-RNA complex. PSF undergoes an all or nothing transition between having two or no RNA-binding pockets. We predict that the RNA binds with a high degree of positive cooperativity. These structures provide an insight into a new regulatory mechanism that is likely involved in promoting malignancies and other human diseases.

Pre-mRNA splicing factor U2AF2 recognizes distinct conformations of nucleotide variants at the center of the pre-mRNA splice site signal

The essential pre-mRNA splicing factor U2AF2 (also called U2AF65) identifies polypyrimidine (Py) tract signals of nascent transcripts, despite length and sequence variations. Previous studies have shown that the U2AF2 RNA recognition motifs (RRM1 and RRM2) preferentially bind uridine-rich RNAs. Nonetheless, the specificity of the RRM1/RRM2 interface for the central Py tract nucleotide has yet to be investigated. We addressed this question by determining crystal structures of U2AF2 bound to a cytidine, guanosine, or adenosine at the central position of the Py tract, and compared U2AF2-bound uridine structures. Local movements of the RNA site accommodated the different nucleotides, whereas the polypeptide backbone remained similar among the structures. Accordingly, molecular dynamics simulations revealed flexible conformations of the central, U2AF2-bound nucleotide. The RNA binding affinities and splicing efficiencies of structure-guided mutants demonstrated that U2AF2 tolerates nucleotide substitutions at the central position of the Py tract. Moreover, enhanced UV-crosslinking and immunoprecipitation of endogenous U2AF2 in human erythroleukemia cells showed uridine-sensitive binding sites, with lower sequence conservation at the central nucleotide positions of otherwise uridine-rich, U2AF2-bound splice sites. Altogether, these results highlight the importance of RNA flexibility for protein recognition and take a step towards relating splice site motifs to pre-mRNA splicing efficiencies.

Somatic Mutations in Core Spliceosome Components Promote Tumorigenesis and Generate an Exploitable Vulnerability in Human Cancer

Simple Summary

High throughput exome sequencing approaches have disclosed recurrent cancer-associated mutations in spliceosomal components, which drive aberrant pre-mRNA processing events and support the tumor phenotype. At the same time, mutations in spliceosome genes and aberrant splicing regulation establish a selective vulnerability of cancer cells to splicing-targeting approaches, which could be exploited therapeutically. It is conceivable that a better understanding of the mechanisms and roles of abnormal splicing in tumor metabolism will facilitate the development of a novel generation of tumor-targeting drugs. In this review, we describe recent advances in the elucidation of the biological impact and biochemical effects of somatic mutations in core spliceosome components on splicing choices and their associated targetable vulnerabilities.

Abstract

Alternative pre-mRNA processing enables the production of distinct mRNA and protein isoforms from a single gene, thus greatly expanding the coding potential of eukaryotic genomes and fine-tuning gene expression programs. Splicing is carried out by the spliceosome, a complex molecular machinery which assembles step-wise on mRNA precursors in the nucleus of eukaryotic cells. In the last decade, exome sequencing technologies have allowed the identification of point mutations in genes encoding splicing factors as a recurrent hallmark of human cancers, with higher incidence in hematological malignancies. These mutations lead to production of splicing factors that reduce the fidelity of the splicing process and yield splicing variants that are often advantageous for cancer cells. However, at the same time, these mutations increase the sensitivity of transformed cells to splicing inhibitors, thus offering a therapeutic opportunity for novel targeted strategies. Herein, we review the recent literature documenting cancer-associated mutations in components of the early spliceosome complex and discuss novel therapeutic strategies based on small-molecule spliceosome inhibitors that exhibit strong anti-tumor effects, particularly against cancer cells harboring mutations in spliceosomal components.

RBM22, a Key Player of Pre-mRNA Splicing and Gene Expression Regulation, Is Altered in Cancer

RNA-Binding Proteins (RBP) are very diverse and cover a large number of functions in the cells. This review focuses on RBM22, a gene encoding an RBP and belonging to the RNA-Binding Motif (RBM) family of genes. RBM22 presents a Zinc Finger like and a Zinc Finger domain, an RNA-Recognition Motif (RRM), and a Proline-Rich domain with a general structure suggesting a fusion of two yeast genes during evolution: Cwc2 and Ecm2. RBM22 is mainly involved in pre-mRNA splicing, playing the essential role of maintaining the conformation of the catalytic core of the spliceosome and acting as a bridge between the catalytic core and other essential protein components of the spliceosome. RBM22 is also involved in gene regulation, and is able to bind DNA, acting as a bona fide transcription factor on a large number of target genes. Undoubtedly due to its wide scope in the regulation of gene expression, RBM22 has been associated with several pathologies and, notably, with the aggressiveness of cancer cells and with the phenotype of a myelodysplastic syndrome. Mutations, enforced expression level, and haploinsufficiency of RBM22 gene are observed in those diseases. RBM22 could represent a potential therapeutic target in specific diseases, and, notably, in cancer.

Phylogenetic Analysis of the Plant U2 snRNP Auxiliary Factor Large Subunit A Gene Family in Response to Developmental Cues and Environmental Stimuli

In all organisms, splicing occurs through the formation of spliceosome complexes, and splicing auxiliary factors are essential during splicing. U2AF65 is a crucial splicing cofactor, and the two typical RNA-recognition motifs at its center recognize and bind the polypyrimidine sequence located between the intron branch site and the 3'-splice site. U2AF65A is a member of the U2AF65 gene family, with pivotal roles in diseases in mammals, specifically humans; however, few studies have investigated plant U2AF65A, and its specific functions are poorly understood. Therefore, in the present study, we systematically identified U2AF65A in plant species from algae to angiosperms. Based on 113 putative U2AF65A sequences from 33 plant species, phylogenetic analyses were performed, followed by basic bioinformatics, including the comparisons of gene structure, protein domains, promoter motifs, and gene expression levels. In addition, using rice as the model crop, we demonstrated that the OsU2AF65A protein is localized to the nucleus and cytoplasm, and it is involved in responses to various stresses, such as drought, high salinity, low temperature, and heavy metal exposure (e.g., cadmium). Using Arabidopsis thaliana and rice mutants, we demonstrated that U2AF65A is involved in the accumulation of plant biomass, growth of hypocotyl upon thermal stimulation, and reduction of tolerance of high temperature stress. These findings offer an overview of the U2AF65 gene family and its stress response functions, serving as the reference for further comprehensive functional studies of the essential specific splicing cofactor U2AF65A in the plant kingdom.

A synthetic small molecule stalls pre-mRNA splicing by promoting an early-stage U2AF2-RNA complex

Dysregulated pre-mRNA splicing is an emerging Achilles heel of cancers and myelodysplasias. To expand the currently limited portfolio of small-molecule drug leads, we screened for chemical modulators of the U2AF complex, which nucleates spliceosome assembly and is mutated in myelodysplasias. A hit compound specifically enhances RNA binding by a U2AF2 subunit. Remarkably, the compound inhibits splicing of representative substrates and stalls spliceosome assembly at the stage of U2AF function. Computational docking, together with structure-guided mutagenesis, indicates that the compound bridges the tandem U2AF2 RNA recognition motifs via hydrophobic and electrostatic moieties. Cells expressing a cancer-associated U2AF1 mutant are preferentially killed by treatment with the compound. Altogether, our results highlight the potential of trapping early spliceosome assembly as an effective pharmacological means to manipulate pre-mRNA splicing. By extension, we suggest that stabilizing assembly intermediates may offer a useful approach for small-molecule inhibition of macromolecular machines.

Computational and Mass Spectrometry-Based Approach Identify Deleterious Non-Synonymous Single Nucleotide Polymorphisms (nsSNPs) in JMJD6

The jumonji domain-containing protein 6 (JMJD6) gene catalyzes the arginine demethylation and lysine hydroxylation of histone and a growing list of its known substrate molecules, including p53 and U2AF65, suggesting a possible role in mRNA splicing and transcription in cancer progression. Mass spectrometry-based technology offers the opportunity to detect SNP variants accurately and effectively. In our study, we conducted a combined computational and filtration workflow to predict the nonsynonymous single nucleotide polymorphisms (nsSNPs) present in JMJD6, followed by a liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis and validation. The computational approaches SIFT, PolyPhen-2, SNAP, I-Mutant 2.0, PhD-SNP, PANTHER, and SNPS&GO were integrated to screen out the predicted damaging/deleterious nsSNPs. Through the three-dimensional structure of JMJD6, H187R (rs1159480887) was selected as a candidate for validation. The validation experiments showed that the mutation of this nsSNP in JMJD6 obviously affected mRNA splicing or the transcription of downstream genes through the reduced lysyl-hydroxylase activity of its substrates, U2AF65 and p53, further indicating the accuracy of this prediction method. This research provides an effective computational workflow for researchers with an opportunity to select prominent deleterious nsSNPs and, thus, remains promising for examining the dysfunction of proteins.

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.

Global developmental delay, systemic dysmorphism and epilepsy in a patient with a de novo U2AF2 variant

U2 small nuclear RNA auxiliary factor 2 (U2AF2) is an essential pre-mRNA splicing factor in an early step of splicing. Alternative splicing plays an important role in neuronal development, and disorders of RNA processing steps are implicated in neurological disorders. Recently, the large trio whole-exome sequencing study reported U2AF2 as a novel gene significantly associated with developmental disorders: however, the clinical details of patients with U2AF2 variants were not available. Here, we report an individual with a de novo U2AF2 variant (c.445C>T, p.(Arg149Trp)) using trio-based whole-exome sequencing. This residue was positioned in the RNA recognition motif 1 which recognizes a polypyrimidine-tract splice site signal. The patient showed global developmental delay, intellectual disability, epilepsy, short stature, microcephaly, facial dysmorphism, intermittent exotropia, bilateral ptosis, muscle hypotonia and thin corpus callosum, indicating that U2AF2-related disorder could include systemic dysmorphisms, epilepsy and brain malformation along with global developmental delay.

Splice site m6A methylation prevents binding of U2AF35 to inhibit RNA splicing

The N⁶-methyladenosine (m⁶A) RNA modification is used widely to alter the fate of mRNAs. Here we demonstrate that the C. elegans writer METT-10 (the ortholog of mouse METTL16) deposits an m⁶A mark on the 3′ splice site (AG) of the S-adenosylmethionine (SAM) synthetase pre-mRNA, which inhibits its proper splicing and protein production. The mechanism is triggered by a rich diet and acts as an m⁶A-mediated switch to stop SAM production and regulate its homeostasis. Although the mammalian SAM synthetase pre-mRNA is not regulated via this mechanism, we show that splicing inhibition by 3′ splice site m⁶A is conserved in mammals. The modification functions by physically preventing the essential splicing factor U2AF35 from recognizing the 3′ splice site. We propose that use of splice-site m⁶A is an ancient mechanism for splicing regulation.

•

m⁶A deposited at 3′ splice site by worm METT-10 inhibits splicing

•

Methylation blocks 3′ splice site recognition by splicing factor U2AF35

•

Methylation and splicing inhibition is a response to change in worm diet

•

Splicing inhibition by 3′ splice site m⁶A is conserved in mammals

The Alter Retina: Alternative Splicing of Retinal Genes in Health and Disease

Alternative splicing of mRNA is an essential mechanism to regulate and increase the diversity of the transcriptome and proteome. Alternative splicing frequently occurs in a tissue- or time-specific manner, contributing to differential gene expression between cell types during development. Neural tissues present extremely complex splicing programs and display the highest number of alternative splicing events. As an extension of the central nervous system, the retina constitutes an excellent system to illustrate the high diversity of neural transcripts. The retina expresses retinal specific splicing factors and produces a large number of alternative transcripts, including exclusive tissue-specific exons, which require an exquisite regulation. In fact, a current challenge in the genetic diagnosis of inherited retinal diseases stems from the lack of information regarding alternative splicing of retinal genes, as a considerable percentage of mutations alter splicing or the relative production of alternative transcripts. Modulation of alternative splicing in the retina is also instrumental in the design of novel therapeutic approaches for retinal dystrophies, since it enables precision medicine for specific mutations.

Comprehensive characterisation of intronic mis-splicing mutations in human cancers

Previous studies studying mis-splicing mutations were based on exome data and thus our current knowledge is largely limited to exons and the canonical splice sites. To comprehensively characterise intronic mis-splicing mutations, we analysed 1134 pan-cancer whole genomes and transcriptomes together with 3022 normal control samples. The ratio-based splicing analysis resulted in 678 somatic intronic mutations, with 46% residing in deep introns. Among the 309 deep intronic single nucleotide variants, 245 altered core splicing codes, with 38% activating cryptic splice sites, 12% activating cryptic polypyrimidine tracts, and 36% and 12% disrupting authentic polypyrimidine tracts and branchpoints, respectively. All the intronic cryptic splice sites were created at pre-existing GT/AG dinucleotides or by GC-to-GT conversion. Notably, 85 deep intronic mutations indicated gain of splicing enhancers or loss of splicing silencers. We found that 64 tumour suppressors were affected by intronic mutations and blood cancers showed higher proportion of deep intronic mutations. In particular, a telomere maintenance gene, POT1, was recurrently mis-spliced by deep intronic mutations in blood cancers. We validated a pseudoexon activation involving a splicing silencer in POT1 by CRISPR/Cas9. Our results shed light on previously unappreciated mechanisms by which noncoding mutations acting on splicing codes in deep introns contribute to tumourigenesis.

Unraveling the mechanism of recognition of the 3' splice site of the adenovirus major late promoter intron by the alternative splicing factor PUF60

Pre-mRNA splicing is critical for achieving required amounts of a transcript at a given time and for regulating production of encoded protein. A given pre-mRNA may be spliced in many ways, or not at all, giving rise to multiple gene products. Numerous splicing factors are recruited to pre-mRNA splice sites to ensure proper splicing. One such factor, the 60 kDa poly(U)-binding splicing factor (PUF60), is recruited to sites that are not always spliced, but rather function as alternative splice sites. In this study, we characterized the interaction of PUF60 with a splice site from the adenovirus major late promoter (the AdML 3' splice site, AdML3'). We found that the PUF60-AdML3' dissociation constants are in the micromolar range, with the binding affinity predominantly provided by PUF60's two central RNA recognition motifs (RRMs). A 1.95 Å crystal structure of the two PUF60 RRMs in complex with AdML3' revealed a dimeric organization placing two stretches of nucleic acid tracts in opposing directionalities, which can cause looping of nucleic acid and explain how PUF60 affects pre-mRNA geometry to effect splicing. Solution characterization of this complex by light-scattering and UV/Vis spectroscopy suggested a potential 2:1 (PUF602:AdML3') stoichiometry, consistent with the crystal structure. This work defines the sequence specificity of the alternative splicing factor PUF60 at the pre-mRNA 3' splice site. Our observations suggest that control of pre-mRNA directionality is important in the early stage of spliceosome assembly, and advance our understanding of the molecular mechanism by which alternative and constitutive splicing factors differentiate among 3' splice sites.

A model of k-mer surprisal to quantify local sequence information content surrounding splice regions

Molecular sequences carry information. Analysis of sequence conservation between homologous loci is a proven approach with which to explore the information content of molecular sequences. This is often done using multiple sequence alignments to support comparisons between homologous loci. These methods therefore rely on sufficient underlying sequence similarity with which to construct a representative alignment. Here we describe a method using a formal metric of information, surprisal, to analyse biological sub-sequences without alignment constraints. We applied our model to the genomes of five different species to reveal similar patterns across a panel of eukaryotes. As the surprisal of a sub-sequence is inversely proportional to its occurrence within the genome, the optimal size of the sub-sequences was selected for each species under consideration. With the model optimized, we found a strong correlation between surprisal and CG dinucleotide usage. The utility of our model was tested by examining the sequences of genes known to undergo splicing. We demonstrate that our model can identify biological features of interest such as known donor and acceptor sites. Analysis across all annotated coding exon junctions in Homo sapiens reveals the information content of coding exons to be greater than the surrounding intron regions, a consequence of increased suppression of the CG dinucleotide in intronic space. Sequences within coding regions proximal to exon junctions exhibited novel patterns within DNA and coding mRNA that are not a function of the encoded amino acid sequence. Our findings are consistent with the presence of secondary information encoding features such as DNA and RNA binding sites, multiplexed through the coding sequence and independent of the information required to define the corresponding amino-acid sequence. We conclude that surprisal provides a complementary methodology with which to locate regions of interest in the genome, particularly in situations that lack an appropriate multiple sequence alignment.

Autosomal dominant retinitis pigmentosa with toxic gain of function: Mechanisms and therapeutics

Autosomal dominant retinitis pigmentosa is a form of retinitis pigmentosa, an inherited retinal degenerative disorder characterized by progressive loss of photoreceptors eventually leading to irreversible loss of vision. Mutations in genes involved in the basic functions of the visual system give rise to this condition. These mutations can either lead to loss of function or toxic gain of function phenotypes. While autosomal dominant retinitis pigmentosa caused by loss of function can be ideally treated by gene supplementation with a single vector to address a different spectrum of mutations in a gene, the same strategy cannot be applied to toxic gain of function phenotypes. In toxic gain of function phenotypes, the mutation in the gene results in the acquisition of a new function that can interrupt the functioning of the wildtype protein by various mechanisms leading to cell toxicity, thus making a single approach impractical. This review focuses on the genes and mechanisms that cause toxic gain of function phenotypes associated with autosomal dominant retinitis pigmentosa and provide a bird's eye view on current therapeutic strategies and ongoing clinical trials.

Systematic characterization of the branch point binding protein, splicing factor 1, gene family in plant development and stress responses

BACKGROUND

Among eukaryotic organisms, alternative splicing is an important process that can generate multiple transcripts from one same precursor messenger RNA, which greatly increase transcriptome and proteome diversity. This process is carried out by a super-protein complex defined as the spliceosome. Specifically, splicing factor 1/branchpoint binding protein (SF1/BBP) is a single protein that can bind to the intronic branchpoint sequence (BPS), connecting the 5' and 3' splice site binding complexes during early spliceosome assembly. The molecular function of this protein has been extensively investigated in yeast, metazoa and mammals. However, its counterpart in plants has been seldomly reported.

RESULTS

To this end, we conducted a systematic characterization of the SF1 gene family across plant lineages. In this work, a total of 92 sequences from 59 plant species were identified. Phylogenetic relationships of these sequences were constructed, and subsequent bioinformatic analysis suggested that this family likely originated from an ancient gene transposition duplication event. Most plant species were shown to maintain a single copy of this gene. Furthermore, an additional RNA binding motif (RRM) existed in most members of this gene family in comparison to their animal and yeast counterparts, indicating that their potential role was preserved in the plant lineage.

CONCLUSION

Our analysis presents general features of the gene and protein structure of this splicing factor family and will provide fundamental information for further functional studies in plants.

A splice site-sensing conformational switch in U2AF2 is modulated by U2AF1 and its recurrent myelodysplasia-associated mutation

An essential heterodimer of the U2AF1 and U2AF2 pre-mRNA splicing factors nucleates spliceosome assembly at polypyrimidine (Py) signals preceding the major class of 3′ splice sites. U2AF1 frequently acquires an S34F-encoding mutation among patients with myelodysplastic syndromes (MDS). The influence of the U2AF1 subunit and its S34F mutation on the U2AF2 conformations remains unknown. Here, we employ single molecule Förster resonance energy transfer (FRET) to determine the influence of wild-type or S34F-substituted U2AF1 on the conformational dynamics of U2AF2 and its splice site RNA complexes. In the absence of RNA, the U2AF1 subunit stabilizes a high FRET value, which by structure-guided mutagenesis corresponds to a closed conformation of the tandem U2AF2 RNA recognition motifs (RRMs). When the U2AF heterodimer is bound to a strong, uridine-rich splice site, U2AF2 switches to a lower FRET value characteristic of an open, side-by-side arrangement of the RRMs. Remarkably, the U2AF heterodimer binds weak, uridine-poor Py tracts as a mixture of closed and open U2AF2 conformations, which are modulated by the S34F mutation. Shifts between open and closed U2AF2 may underlie U2AF1-dependent splicing of degenerate Py tracts and contribute to a subset of S34F-dysregulated splicing events in MDS patients.

Cancer-Associated Substitutions in RNA Recognition Motifs of PUF60 and U2AF65 Reveal Residues Required for Correct Folding and 3' Splice-Site Selection

U2AF65 (U2AF2) and PUF60 (PUF60) are splicing factors important for recruitment of the U2 small nuclear ribonucleoprotein to lariat branch points and selection of 3′ splice sites (3′ss). Both proteins preferentially bind uridine-rich sequences upstream of 3′ss via their RNA recognition motifs (RRMs). Here, we examined 36 RRM substitutions reported in cancer patients to identify variants that alter 3′ss selection, RNA binding and protein properties. Employing PUF60- and U2AF65-dependent 3′ss previously identified by RNA-seq of depleted cells, we found that 43% (10/23) and 15% (2/13) of independent RRM mutations in U2AF65 and PUF60, respectively, conferred splicing defects. At least three RRM mutations increased skipping of internal U2AF2 (~9%, 2/23) or PUF60 (~8%, 1/13) exons, indicating that cancer-associated RRM mutations can have both cis- and trans-acting effects on splicing. We also report residues required for correct folding/stability of each protein and map functional RRM substitutions on to existing high-resolution structures of U2AF65 and PUF60. These results identify new RRM residues critical for 3′ss selection and provide relatively simple tools to detect clonal RRM mutations that enhance the mRNA isoform diversity.

mRNA Editing, Processing and Quality Control in Caenorhabditis elegans

While DNA serves as the blueprint of life, the distinct functions of each cell are determined by the dynamic expression of genes from the static genome. The amount and specific sequences of RNAs expressed in a given cell involves a number of regulated processes including RNA synthesis (transcription), processing, splicing, modification, polyadenylation, stability, translation, and degradation. As errors during mRNA production can create gene products that are deleterious to the organism, quality control mechanisms exist to survey and remove errors in mRNA expression and processing. Here, we will provide an overview of mRNA processing and quality control mechanisms that occur in Caenorhabditis elegans, with a focus on those that occur on protein-coding genes after transcription initiation. In addition, we will describe the genetic and technical approaches that have allowed studies in C. elegans to reveal important mechanistic insight into these processes.

Cis- and trans-regulations of pre-mRNA splicing by RNA editing enzymes influence cancer development

RNA editing and splicing are the two major processes that dynamically regulate human transcriptome diversity. Despite growing evidence of crosstalk between RNA editing enzymes (mainly ADAR1) and splicing machineries, detailed mechanistic explanations and their biological importance in diseases, such as cancer are still lacking. Herein, we identify approximately a hundred high-confidence splicing events altered by ADAR1 and/or ADAR2, and ADAR1 or ADAR2 protein can regulate cassette exons in both directions. We unravel a binding tendency of ADARs to dsRNAs that involves GA-rich sequences for editing and splicing regulation. ADAR1 edits an intronic splicing silencer, leading to recruitment of SRSF7 and repression of exon inclusion. We also present a mechanism through which ADAR2 binds to dsRNA formed between GA-rich sequences and polypyrimidine (Py)-tract and precludes access of U2AF65 to 3' splice site. Furthermore, we find these ADARs-regulated splicing changes per se influence tumorigenesis, not merely byproducts of ADARs editing and binding.

Genome-wide evolution of wobble base-pairing nucleotides of branchpoint motifs with increasing organismal complexity

How have the branchpoint motifs evolved in organisms of different complexity? Here we identified and examined the consensus motifs (R₁C₂T₃R₄A₅Y₆, R: A or G, Y: C or T) of 898 fungal genomes. In Ascomycota unicellular yeasts, the G₄/A₄ ratio is mostly (98%) below 0.125 but increases sharply in multicellular species by about 40 times on average, and in the more complex Basidiomycota, it increases further by about 7 times. The global G₄ increase is consistent with A₄ to G₄ transitions in evolution. Of the G₄/A₄-interacting amino acids of the branchpoint binding protein MSL5 (SF1) and the HSH155 (SF3B1), as well as the 5' splice sites (SS) and U2 snRNA genes, the 5' SS G³/A³ co-vary with the G₄ to some extent. However, corresponding increase of the G₄-complementary GCAGTA-U2 gene is rare, suggesting wobble-base pairing between the G₄-containing branchpoint motif and GTAGTA-U2 in most of these species. Interestingly, the G₄/A₄ ratio correlates well with the abundance of alternative splicing in the two phyla, and G₄ enriched significantly at the alternative 3' SS of genes in RNA metabolism, kinases and membrane proteins. Similar wobble nucleotides also enriched at the 3' SS of multicellular fungi with only thousands of protein-coding genes. Thus, branchpoint motifs have evolved U2-complementarity in unicellular Ascomycota yeasts, but have gradually gained more wobble base-pairing nucleotides in fungi of higher complexity, likely to destabilize branchpoint motif-U2 interaction and/or branchpoint A protrusion for alternative splicing. This implies an important role of relaxing the branchpoint signals in the multicellularity and further complexity of fungi.

Cancer the'RBP'eutics-RNA-binding proteins as therapeutic targets for cancer

RNA-binding proteins (RBPs) play a critical role in the regulation of various RNA processes, including splicing, cleavage and polyadenylation, transport, translation and degradation of coding RNAs, non-coding RNAs and microRNAs. Recent studies indicate that RBPs not only play an instrumental role in normal cellular processes but have also emerged as major players in the development and spread of cancer. Herein, we review the current knowledge about RNA binding proteins and their role in tumorigenesis as well as the potential to target RBPs for cancer therapeutics.

HPV18 Utilizes Two Alternative Branch Sites for E6*I Splicing to Produce E7 Protein

Human papillomavirus 18 (HPV18) E6 and E7 oncogenes are transcribed as a single bicistronic E6E7 pre-mRNA. The E6 ORF region in the bicistronic E6E7 pre-mRNA contains an intron. Splicing of this intron disrupts the E6 ORF integrity and produces a spliced E6*I RNA for efficient E7 translation. Here we report that the E6 intron has two overlapped branch point sequences (BPS) upstream of its 3′ splice site, with an identical heptamer AACUAAC, for E6*I splicing. One heptamer has a branch site adenosine (underlined) at nt 384 and the other at nt 388. E6*I splicing efficiency correlates to the expression level of E6 and E7 proteins and depends on the selection of which branch site. In general, E6*I splicing prefers the 3′ss-proximal branch site at nt 388 over the distal branch site at nt 384. Inactivation of the nt 388 branch site was found to activate a cryptic acceptor site at nt 636 for aberrant RNA splicing. Together, these data suggest that HPV18 modulates its production ratio of E6 and E7 proteins by alternative selection of the two mapped branch sites for the E6*I splicing, which could be beneficial in its productive or oncogenic infection according to the host cell environment.

Structural identification of conserved RNA binding sites in herpesvirus ORF57 homologs: implications for PAN RNA recognition

Kaposi's sarcoma-associated herpesvirus (KSHV) transcribes a long noncoding polyadenylated nuclear (PAN) RNA, which promotes the latent to lytic transition by repressing host genes involved in antiviral responses as well as viral proteins that support the latent state. KSHV also expresses several early proteins including ORF57 (Mta), a member of the conserved multifunctional ICP27 protein family, which is essential for productive replication. ORF57/Mta interacts with PAN RNA via a region termed the Mta responsive element (MRE), stabilizing the transcript and supporting nuclear accumulation. Here, using a close homolog of KSHV ORF57 from herpesvirus saimiri (HVS), we determined the crystal structure of the globular domain in complex with a PAN RNA MRE, revealing a uracil specific binding site that is also conserved in KSHV. Using solution NMR, RNA binding was also mapped within the disordered N-terminal domain of KSHV ORF57, and showed specificity for an RNA fragment containing a GAAGRG motif previously known to bind a homologous region in HVS ORF57. Together these data located novel differential RNA recognition sites within neighboring domains of herpesvirus ORF57 homologs, and revealed high-resolution details of their interactions with PAN RNA, thus providing insight into interactions crucial to viral function.

Molecular Mechanisms of pre-mRNA Splicing through Structural Biology of the Spliceosome

Precursor messenger RNA (pre-mRNA) splicing is executed by the spliceosome. In the past 3 years, cryoelectron microscopy (cryo-EM) structures have been elucidated for a majority of the yeast spliceosomal complexes and for a few human spliceosomes. During the splicing reaction, the dynamic spliceosome has an immobile core of about 20 protein and RNA components, which are organized around a conserved splicing active site. The divalent metal ions, coordinated by U6 small nuclear RNA (snRNA), catalyze the branching reaction and exon ligation. The spliceosome also contains a mobile but compositionally stable group of about 13 proteins and a portion of U2 snRNA, which facilitate substrate delivery into the splicing active site. The spliceosomal transitions are driven by the RNA-dependent ATPase/helicases, resulting in the recruitment and dissociation of specific splicing factors that enable the reaction. In summary, the spliceosome is a protein-directed metalloribozyme.

Extensive Alternative Splicing of KIR Transcripts

The killer-cell Ig-like receptors (KIR) form a multigene entity involved in modulating immune responses through interactions with MHC class I molecules. The complexity of the KIR cluster is reflected by, for instance, abundant levels of allelic polymorphism, gene copy number variation, and stochastic expression profiles. The current transcriptome study involving human and macaque families demonstrates that KIR family members are also subjected to differential levels of alternative splicing, and this seems to be gene dependent. Alternative splicing may result in the partial or complete skipping of exons, or the partial inclusion of introns, as documented at the transcription level. This post-transcriptional process can generate multiple isoforms from a single KIR gene, which diversifies the characteristics of the encoded proteins. For example, alternative splicing could modify ligand interactions, cellular localization, signaling properties, and the number of extracellular domains of the receptor. In humans, we observed abundant splicing for KIR2DL4, and to a lesser extent in the lineage III KIR genes. All experimentally documented splice events are substantiated by in silico splicing strength predictions. To a similar extent, alternative splicing is observed in rhesus macaques, a species that shares a close evolutionary relationship with humans. Splicing profiles of Mamu-KIR1D and Mamu-KIR2DL04 displayed a great diversity, whereas Mamu-KIR3DL20 (lineage V) is consistently spliced to generate a homolog of human KIR2DL5 (lineage I). The latter case represents an example of convergent evolution. Although just a single KIR splice event is shared between humans and macaques, the splicing mechanisms are similar, and the predicted consequences are comparable. In conclusion, alternative splicing adds an additional layer of complexity to the KIR gene system in primates, and results in a wide structural and functional variety of KIR receptors and its isoforms, which may play a role in health and disease.

RBM-5 modulates U2AF large subunit-dependent alternative splicing in C. elegans

A key step in pre-mRNA splicing is the recognition of 3' splicing sites by the U2AF large and small subunits, a process regulated by numerous trans-acting splicing factors. How these trans-acting factors interact with U2AF in vivo is unclear. From a screen for suppressors of the temperature-sensitive (ts) lethality of the C. elegans U2AF large subunit gene uaf-1(n4588) mutants, we identified mutations in the RNA binding motif gene rbm-5, a homolog of the tumor suppressor gene RBM5. rbm-5 mutations can suppress uaf-1(n4588) ts-lethality by loss of function and neuronal expression of rbm-5 was sufficient to rescue the suppression. Transcriptome analyses indicate that uaf-1(n4588) affected the expression of numerous genes and rbm-5 mutations can partially reverse the abnormal gene expression to levels similar to that of wild type. Though rbm-5 mutations did not obviously affect alternative splicing per se, they can suppress or enhance, in a gene-specific manner, the altered splicing of genes in uaf-1(n4588) mutants. Specifically, the recognition of a weak 3' splice site was more susceptible to the effect of rbm-5. Our findings provide novel in vivo evidence that RBM-5 can modulate UAF-1-dependent RNA splicing and suggest that RBM5 might interact with U2AF large subunit to affect tumor formation.

BPP: a sequence-based algorithm for branch point prediction

Motivation

Although high-throughput sequencing methods have been proposed to identify splicing branch points in the human genome, these methods can only detect a small fraction of the branch points subject to the sequencing depth, experimental cost and the expression level of the mRNA. An accurate computational model for branch point prediction is therefore an ongoing objective in human genome research.

Results

We here propose a novel branch point prediction algorithm that utilizes information on the branch point sequence and the polypyrimidine tract. Using experimentally validated data, we demonstrate that our proposed method outperforms existing methods. Availability and implementation: https://github.com/zhqingit/BPP.

Contact

djguo@cuhk.edu.hk.

Supplementary information

Supplementary data are available at Bioinformatics online.

Capturing dynamic conformational shifts in protein–ligand recognition using integrative structural biology in solution

In recent years, a dynamic view of the structure and function of biological macromolecules is emerging, highlighting an essential role of dynamic conformational equilibria to understand molecular mechanisms of biological functions. The structure of a biomolecule, i.e. protein or nucleic acid in solution, is often best described as a dynamic ensemble of conformations, rather than a single structural state. Strikingly, the molecular interactions and functions of the biological macromolecule can then involve a shift between conformations that pre-exist in such an ensemble. Upon external cues, such population shifts of pre-existing conformations allow gradually relaying the signal to the downstream biological events. An inherent feature of this principle is conformational dynamics, where intrinsically disordered regions often play important roles to modulate the conformational ensemble. Unequivocally, solution-state NMR spectroscopy is a powerful technique to study the structure and dynamics of such biomolecules in solution. NMR is increasingly combined with complementary techniques, including fluorescence spectroscopy and small angle scattering. The combination of these techniques provides complementary information about the conformation and dynamics in solution and thus affords a comprehensive description of biomolecular functions and regulations. Here, we illustrate how an integrated approach combining complementary techniques can assess the structure and dynamics of proteins and protein complexes in solution.

Universal correction of enzymatic sequence bias reveals molecular signatures of protein/DNA interactions

Coupling molecular biology to high-throughput sequencing has revolutionized the study of biology. Molecular genomics techniques are continually refined to provide higher resolution mapping of nucleic acid interactions and structure. Sequence preferences of enzymes can interfere with the accurate interpretation of these data. We developed seqOutBias to characterize enzymatic sequence bias from experimental data and scale individual sequence reads to correct intrinsic enzymatic sequence biases. SeqOutBias efficiently corrects DNase-seq, TACh-seq, ATAC-seq, MNase-seq and PRO-seq data. We show that seqOutBias correction facilitates identification of true molecular signatures resulting from transcription factors and RNA polymerase interacting with DNA.

snRNP proteins in health and disease

Split gene architecture of most human genes requires removal of intervening sequences by mRNA splicing that occurs on large multiprotein complexes called spliceosomes. Mutations compromising several spliceosomal components have been recorded in degenerative syndromes and haematological neoplasia, thereby highlighting the importance of accurate splicing execution in homeostasis of assorted adult tissues. Moreover, insufficient splicing underlies defective development of craniofacial skeleton and upper extremities. This review summarizes recent advances in the understanding of splicing factor function deduced from cryo-EM structures. We combine these data with the characterization of splicing factors implicated in hereditary or somatic disorders, with a focus on potential functional consequences the mutations may elicit in spliceosome assembly and/or performance. Given aberrant splicing or perturbations in splicing efficiency substantially underpin disease pathogenesis, profound understanding of the mis-splicing principles may open new therapeutic vistas. In three major sections dedicated to retinal dystrophies, hereditary acrofacial syndromes, and haematological malignancies, we delineate the noticeable variety of conditions associated with dysfunctional splicing and accentuate recurrent patterns in splicing defects.

Cancer-Associated Mutations Mapped on High-Resolution Structures of the U2AF2 RNA Recognition Motifs

Acquired point mutations of pre-mRNA splicing factors recur among cancers, leukemias, and related neoplasms. Several studies have established that somatic mutations of a U2AF1 subunit, which normally recognizes 3' splice site junctions, recur among myelodysplastic syndromes. The U2AF2 splicing factor recognizes polypyrimidine signals that precede most 3' splice sites as a heterodimer with U2AF1. In contrast with those of the well-studied U2AF1 subunit, descriptions of cancer-relevant U2AF2 mutations and their structural relationships are lacking. Here, we survey databases of cancer-associated mutations and identify recurring missense mutations in the U2AF2 gene. We determine ultra-high-resolution structures of the U2AF2 RNA recognition motifs (RRM1 and RRM2) at 1.1 Å resolution and map the structural locations of the mutated U2AF2 residues. Comparison with prior, lower-resolution structures of the tandem U2AF2 RRMs in the RNA-bound and apo states reveals clusters of cancer-associated mutations at the U2AF2 RRM-RNA or apo-RRM1-RRM2 interfaces. Although the role of U2AF2 mutations in malignant transformation remains uncertain, our results show that cancer-associated mutations correlate with functionally important surfaces of the U2AF2 splicing factor.