Exon recognition in vertebrate splicing

Timing is everything: advances in quantifying splicing kinetics

Splicing is a highly regulated process critical for proper pre-mRNA maturation and the maintenance of a healthy cellular environment. Splicing events are impacted by ongoing transcription, neighboring splicing events, and cis and trans regulatory factors on the respective pre-mRNA transcript. Within this complex regulatory environment, splicing kinetics have the potential to influence splicing outcomes but have historically been challenging to study in vivo. In this review, we highlight recent technological advancements that have enabled measurements of global splicing kinetics and of the variability of splicing kinetics at single introns. We demonstrate how identifying features that are correlated with splicing kinetics has increased our ability to form potential models for how splicing kinetics may be regulated in vivo.

Branch site recognition by the spliceosome

The spliceosome is a eukaryotic multimegadalton RNA-protein complex that removes introns from transcripts. The spliceosome ensures the selection of each exon-intron boundary through multiple recognition events. Initially, the 5' splice site (5' SS) and branch site (BS) are bound by the U1 small nuclear ribonucleoprotein (snRNP) and the U2 snRNP, respectively, while the 3' SS is mostly determined by proximity to the branch site. A large number of splicing factors recognize the splice sites and recruit the snRNPs before the stable binding of the snRNPs occurs by base-pairing the snRNA to the transcript. Fidelity of this process is crucial, as mutations in splicing factors and U2 snRNP components are associated with many diseases. In recent years, major advances have been made in understanding how splice sites are selected in Saccharomyces cerevisiae and humans. Here, I review and discuss the current understanding of the recognition of splice sites by the spliceosome with a focus on recognition and binding of the branch site by the U2 snRNP in humans.

OpenASO: RNA Rescue - designing splice-modulating antisense oligonucleotides through community science

Splice-modulating antisense oligonucleotides (ASOs) are precision RNA-based drugs that are becoming an established modality to treat human disease. Previously, we reported the discovery of ASOs that target a novel, putative intronic RNA structure to rescue splicing of multiple pathogenic variants of F8 exon 16 that cause hemophilia A. However, the conventional approach to discovering splice-modulating ASOs is both laborious and expensive. Here, we describe an alternative paradigm that integrates data-driven RNA structure prediction and community science to discover splice-modulating ASOs. Using a splicing-deficient pathogenic variant of F8 exon 16 as a model, we show that 25% of the top-scoring molecules designed in the Eterna OpenASO challenge have a statistically significant impact on enhancing exon 16 splicing. Additionally, we show that a distinct combination of ASOs designed by Eterna players can additively enhance the inclusion of the splicing-deficient exon 16 variant. Together, our data suggests that crowdsourcing designs from a community of citizen scientists may accelerate the discovery of splice-modulating ASOs with potential to treat human disease.

Conserved role for spliceosomal component PRPF40A in microexon splicing

Microexons (exons ≤30 nts) are important features of neuronal transcriptomes, but pose mechanistic challenges to the splicing machinery. We previously showed that PRP-40, a component of the U1 spliceosome, is globally required for microexon splicing in C. elegans. Here we show that the homologous PRPF40A is also globally required for microexon splicing in mouse neuroblastoma cells. We find that PRPF40A co-regulates microexons along with SRRM4, a neuron-specific regulator of microexon splicing. The relationship between exon size and dependence on PRPF40A/SRRM4 is distinct, with SRRM4-dependence exhibiting a size threshold (~30 nts) and PRPF40A-dependence exhibiting a graded decrease as exon size increases. Finally, we show that PRPF40A knockdown causes an increase in productive splicing of its spliceosomal binding partner Luc7l by skipping of a small "poison exon." Similar homeostatic cross-regulation is often observed across paralogous RNA binding proteins. Here we find this concept likewise applies across evolutionarily unrelated but functionally and physically coupled spliceosomal components.

Study of the RNA splicing kinetics via in vivo 5-EU labeling

Splicing is an important step of gene expression in all eukaryotes. Splice sites might be used with different efficiency, giving rise to alternative splicing products. At the same time, splice sites might be used at a variable rate. We used 5-ethynyl uridine labeling to sequence a nascent transcriptome of HeLa cells and deduced the rate of splicing for each donor and acceptor splice site. The following correlation analysis showed a correspondence of primary transcript features with the rate of splicing. Some dependencies we revealed were anticipated, such as a splicing rate decrease with a decreased complementarity of the donor splice site to U1 and acceptor sites to U2 snRNAs. Other dependencies were more surprising, like a negative influence of a distance to the 5' end on the rate of the acceptor splicing site utilization, or the differences in splicing rate between long, short, and RBM17-dependent introns. We also observed a deceleration of last intron splicing with an increase of the distance to the poly(A) site, which might be explained by the cooperativity of the splicing and polyadenylation. Additional analysis of splicing kinetics of SF3B4 knockdown cells suggested the impairment of a U2 snRNA recognition step. As a result, we deconvoluted the effects of several examined features on the splicing rate into a single regression model. The data obtained here are useful for further studies in the field, as they provide general splicing rate dependencies as well as help to justify the existence of slowly removed splice sites.

Single-Base Substitution Causing Dual-Exon Skipping Event in PKD2 Gene: Unusual Molecular Finding from Exome Sequencing in a Patient with Autosomal Dominant Polycystic Kidney Disease

Background: Pathogenic variants in the Polycystic Kidney Disease 2 (PKD2) gene are associated with Autosomal Dominant Polycystic Kidney Disease (ADPKD) in approximately 30% of cases. In recent years, the high-throughput sequencing techniques have significantly increased the number of variants identified in affected patients. Here, we described the peculiar effect of a PKD2 splicing variant, the c.1717-2A>G, identified in an Italian male patient with ADPKD. This variant led to the unusual and rare skipping of two consecutive exons, causing a large in-frame deletion. Methods: The genetic evaluation of the patient was performed using the Next-Generation Sequencing (NGS) assay Clinical Exome Solution® (SOPHiA Genetics). Bioinformatics analysis was performed using the SOPHiA DDM platform (SOPHiA Genetics). Prediction of pathogenicity was carried out by integrating several in silico tools. RNA evaluation was performed to test the effect of the variant on the PKD2 splicing using a Reverse-Transcription PCR coupled with cDNA sequencing. Results: NGS revealed the presence of the PKD2 c.1717-2A>G variant that lies in the canonical splice site of intron 7. This rare variant was predicted to have a significant impact on the splicing, proved by the RNA-based analysis. We identified the presence of a transcript characterised by the simultaneous skipping of exons 8 and 9, with a retained reading frame and the merging of exons 7-10. Conclusions: We described for the first time a dual-exon skip event related to the presence of a single-base substitution in the PKD2 gene in an ADPKD-affected patient. We assumed that the molecular basis of such a rare mechanism lies in the specific order of intron removal. The finding represents novel evidence of an alternative and unusual splicing mechanism in the PKD2 gene, adding insights to the pathogenesis of the ADPKD.

Mechanisms and regulation of spliceosome-mediated pre-mRNA splicing in Saccharomyces cerevisiae

Pre-mRNA splicing, the removal of introns and ligation of flanking exons, is a crucial step in eukaryotic gene expression. The spliceosome, a macromolecular complex made up of five small nuclear RNAs (snRNAs) and dozens of proteins, assembles on introns via a complex pathway before catalyzing the two transesterification reactions necessary for splicing. All of these steps have the potential to be highly regulated to ensure correct mRNA isoform production for proper cellular function. While Saccharomyces cerevisiae (yeast) has a limited set of intron-containing genes, many of these genes are highly expressed, resulting in a large number of transcripts in a cell being spliced. As a result, splicing regulation is of critical importance for yeast. Just as in humans, yeast splicing can be influenced by protein components of the splicing machinery, structures and properties of the pre-mRNA itself, or by the action of trans-acting factors. It is likely that further analysis of the mechanisms and pathways of splicing regulation in yeast can reveal general principles applicable to other eukaryotes. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

Co-transcriptional gene regulation in eukaryotes and prokaryotes

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascent RNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

Structural insights into the cross-exon to cross-intron spliceosome switch

Early spliceosome assembly can occur through an intron-defined pathway, whereby U1 and U2 small nuclear ribonucleoprotein particles (snRNPs) assemble across the intron1. Alternatively, it can occur through an exon-defined pathway2-5, whereby U2 binds the branch site located upstream of the defined exon and U1 snRNP interacts with the 5' splice site located directly downstream of it. The U4/U6.U5 tri-snRNP subsequently binds to produce a cross-intron (CI) or cross-exon (CE) pre-B complex, which is then converted to the spliceosomal B complex6,7. Exon definition promotes the splicing of upstream introns2,8,9 and plays a key part in alternative splicing regulation10-16. However, the three-dimensional structure of exon-defined spliceosomal complexes and the molecular mechanism of the conversion from a CE-organized to a CI-organized spliceosome, a pre-requisite for splicing catalysis, remain poorly understood. Here cryo-electron microscopy analyses of human CE pre-B complex and B-like complexes reveal extensive structural similarities with their CI counterparts. The results indicate that the CE and CI spliceosome assembly pathways converge already at the pre-B stage. Add-back experiments using purified CE pre-B complexes, coupled with cryo-electron microscopy, elucidate the order of the extensive remodelling events that accompany the formation of B complexes and B-like complexes. The molecular triggers and roles of B-specific proteins in these rearrangements are also identified. We show that CE pre-B complexes can productively bind in trans to a U1 snRNP-bound 5' splice site. Together, our studies provide new mechanistic insights into the CE to CI switch during spliceosome assembly and its effect on pre-mRNA splice site pairing at this stage.

Structural insights into human exon-defined spliceosome prior to activation

Spliceosome is often assembled across an exon and undergoes rearrangement to span a neighboring intron. Most states of the intron-defined spliceosome have been structurally characterized. However, the structure of a fully assembled exon-defined spliceosome remains at large. During spliceosome assembly, the pre-catalytic state (B complex) is converted from its precursor (pre-B complex). Here we report atomic structures of the exon-defined human spliceosome in four sequential states: mature pre-B, late pre-B, early B, and mature B. In the previously unknown late pre-B state, U1 snRNP is already released but the remaining proteins are still in the pre-B state; unexpectedly, the RNAs are in the B state, with U6 snRNA forming a duplex with 5'-splice site and U5 snRNA recognizing the 3'-end of the exon. In the early and mature B complexes, the B-specific factors are stepwise recruited and specifically recognize the exon 3'-region. Our study reveals key insights into the assembly of the exon-defined spliceosomes and identifies mechanistic steps of the pre-B-to-B transition.

Structural insights into branch site proofreading by human spliceosome

Selection of the pre-mRNA branch site (BS) by the U2 small nuclear ribonucleoprotein (snRNP) is crucial to prespliceosome (A complex) assembly. The RNA helicase PRP5 proofreads BS selection but the underlying mechanism remains unclear. Here we report the atomic structures of two sequential complexes leading to prespliceosome assembly: human 17S U2 snRNP and a cross-exon pre-A complex. PRP5 is anchored on 17S U2 snRNP mainly through occupation of the RNA path of SF3B1 by an acidic loop of PRP5; the helicase domain of PRP5 associates with U2 snRNA; the BS-interacting stem-loop (BSL) of U2 snRNA is shielded by TAT-SF1, unable to engage the BS. In the pre-A complex, an initial U2-BS duplex is formed; the translocated helicase domain of PRP5 stays with U2 snRNA and the acidic loop still occupies the RNA path. The pre-A conformation is specifically stabilized by the splicing factors SF1, DNAJC8 and SF3A2. Cancer-derived mutations in SF3B1 damage its association with PRP5, compromising BS proofreading. Together, these findings reveal key insights into prespliceosome assembly and BS selection or proofreading by PRP5.

Tumor Antigens beyond the Human Exome

With the advent of immunotherapeutics, a new era in the combat against cancer has begun. Particularly promising are neo-epitope-targeted therapies as the expression of neo-antigens is tumor-specific. In turn, this allows the selective targeting and killing of cancer cells whilst healthy cells remain largely unaffected. So far, many advances have been made in the development of treatment options which are tailored to the individual neo-epitope repertoire. The next big step is the achievement of efficacious "off-the-shelf" immunotherapies. For this, shared neo-epitopes propose an optimal target. Given the tremendous potential, a thorough understanding of the underlying mechanisms which lead to the formation of neo-antigens is of fundamental importance. Here, we review the various processes which result in the formation of neo-epitopes. Broadly, the origin of neo-epitopes can be categorized into three groups: canonical, noncanonical, and viral neo-epitopes. For the canonical neo-antigens that arise in direct consequence of somatic mutations, we summarize past and recent findings. Beyond that, our main focus is put on the discussion of noncanonical and viral neo-epitopes as we believe that targeting those provides an encouraging perspective to shape the future of cancer immunotherapeutics.

Deciphering the impact of coding and non-coding SCN1A gene variants on RNA splicing

Variants that disrupt normal pre-mRNA splicing are increasingly being recognized as a major cause of monogenic disorders. The SCN1A gene, a key epilepsy gene that is linked to various epilepsy phenotypes, is no exception. Approximately 10% of all reported variants in the SCN1A gene are designated as splicing variants, with many located outside of the canonical donor and acceptor splice sites, and most have not been functionally investigated. However, given its restricted expression pattern, functional analysis of splicing variants in the SCN1A gene could not be routinely performed. In this study, we conducted a comprehensive analysis of all reported SCN1A variants and their potential to impact SCN1A splicing and conclude that splicing variants are substantially misannotated and under-represented. We created a splicing reporter system consisting of 18 splicing vectors covering all 26 protein-coding exons with different genomic contexts and several promoters of varying strengths in order to reproduce the wild-type splicing pattern of the SCN1A gene, revealing cis-regulatory elements essential for proper recognition of SCN1A exons. Functional analysis of 95 SCN1A variants was carried out, including all 68 intronic variants reported in the literature, located outside of the splice sites canonical dinucleotides; 21 exonic variants of different classes (synonymous, missense, nonsense and in-frame deletion) and six variants observed in patients with epilepsy. Interestingly, almost 20% of tested intronic variants had no influence on SCN1A splicing, despite being reported as causative in the literature. Moreover, we confirmed that the majority of predicted exonic variants affect splicing unravelling their true molecular mechanism. We used functional data to perform genotype-phenotype correlation, revealing distinct distribution patterns for missense and splice-affecting 'missense' variants and observed no difference in the phenotype severity of variants leading to in-frame and out-of-frame isoforms, indicating that the Nav1.1 protein is highly intolerant to structural variations. Our work demonstrates the importance of functional analysis in proper variant annotation and provides a tool for high-throughput delineation of splice-affecting variants in SCN1A in a whole-gene manner.

Exonic splicing code and coordination of divalent metals in proteins

Exonic sequences contain both protein-coding and RNA splicing information but the interplay of the protein and splicing code is complex and poorly understood. Here, we have studied traditional and auxiliary splicing codes of human exons that encode residues coordinating two essential divalent metals at the opposite ends of the Irving-Williams series, a universal order of relative stabilities of metal-organic complexes. We show that exons encoding Zn2+-coordinating amino acids are supported much less by the auxiliary splicing motifs than exons coordinating Ca2+. The handicap of the former is compensated by stronger splice sites and uridine-richer polypyrimidine tracts, except for position -3 relative to 3' splice junctions. However, both Ca2+ and Zn2+ exons exhibit close-to-constitutive splicing in multiple tissues, consistent with their critical importance for metalloprotein function and a relatively small fraction of expendable, alternatively spliced exons. These results indicate that constraints imposed by metal coordination spheres on RNA splicing have been efficiently overcome by the plasticity of exon-intron architecture to ensure adequate metalloprotein expression.

RNA elements that control human papillomavirus mRNA splicing-targets for therapy?

Human papillomaviruses (HPVs) cause more than 4.5% of all cancer in the world and more than half of these cases are attributed to human papillomavirus type 16 (HPV16). Prophylactic vaccines are available but antiviral drugs are not. Novel targets for therapy are urgently needed. Alternative RNA splicing is extensively used by HPVs to express all their genes and HPV16 is no exception. This process must function to perfection since mis-splicing could perturb the HPV gene expression program by altering mRNA levels or by generating dysfunctional mRNAs. Cis-acting RNA elements on the viral mRNAs and their cognate cellular trans-acting factors control papillomavirus RNA splicing. The precise but delicate nature of the splicing process renders splicing sensitive to interference. As such, papillomavirus RNA splicing is a potential target for therapy. Here we summarize our current understanding of cis-acting HPV16 RNA elements that control HPV16 mRNA splicing via cellular proteins and discuss how they may be exploited as targets for therapy to papillomavirus infections and cancer.

The human papillomavirus late life cycle and links to keratinocyte differentiation

Regulation of human papillomavirus (HPV) gene expression is tightly linked to differentiation of the keratinocytes the virus infects. HPV late gene expression is confined to the cells in the upper layers of the epithelium where the virus capsid proteins are synthesized. As these proteins are highly immunogenic, and the upper epithelium is an immune-privileged site, this spatial restriction aids immune evasion. Many decades of work have contributed to the current understanding of how this restriction occurs at a molecular level. This review will examine what is known about late gene expression in HPV-infected lesions and will dissect the intricacies of late gene regulation. Future directions for novel antiviral approaches will be highlighted.

An intronic RNA element modulates Factor VIII exon-16 splicing

Pathogenic variants in the human Factor VIII (F8) gene cause Hemophilia A (HA). Here, we investigated the impact of 97 HA-causing single-nucleotide variants on the splicing of 11 exons from F8. For the majority of F8 exons, splicing was insensitive to the presence of HA-causing variants. However, splicing of several exons, including exon-16, was impacted by variants predicted to alter exonic splicing regulatory sequences. Using exon-16 as a model, we investigated the structure-function relationship of HA-causing variants on splicing. Intriguingly, RNA chemical probing analyses revealed a three-way junction structure at the 3'-end of intron-15 (TWJ-3-15) capable of sequestering the polypyrimidine tract. We discovered antisense oligonucleotides (ASOs) targeting TWJ-3-15 partially rescue splicing-deficient exon-16 variants by increasing accessibility of the polypyrimidine tract. The apical stem loop region of TWJ-3-15 also contains two hnRNPA1-dependent intronic splicing silencers (ISSs). ASOs blocking these ISSs also partially rescued splicing. When used in combination, ASOs targeting both the ISSs and the region sequestering the polypyrimidine tract, fully rescue pre-mRNA splicing of multiple HA-linked variants of exon-16. Together, our data reveal a putative RNA structure that sensitizes F8 exon-16 to aberrant splicing.

The Nuclear Cap-Binding Complex, a multitasking binding partner of RNA polymerase II transcripts

In eukaryotic cells, RNAs transcribed by RNA polymerase-II receive the modification at the 5' end. This structure is called the cap structure. The cap structure has a fundamental role for translation initiation by recruiting eukaryotic translation initiation factor 4F (eIF4F). The other important mediator of the cap structure is a nuclear cap-binding protein complex (CBC). CBC consists of two proteins, which are renamed as NCBP1 and NCBP2 (previously called as CBP80/NCBP and CBP20/NIP1, respectively). This review article discusses the multiple roles CBC mediates and co-ordinates in several gene expression steps in eukaryotes.

Deep learning of human polyadenylation sites at nucleotide resolution reveals molecular determinants of site usage and relevance in disease

The genomic distribution of cleavage and polyadenylation (polyA) sites should be co-evolutionally optimized with the local gene structure. Otherwise, spurious polyadenylation can cause premature transcription termination and generate aberrant proteins. To obtain mechanistic insights into polyA site optimization across the human genome, we develop deep/machine learning models to identify genome-wide putative polyA sites at unprecedented nucleotide-level resolution and calculate their strength and usage in the genomic context. Our models quantitatively measure position-specific motif importance and their crosstalk in polyA site formation and cleavage heterogeneity. The intronic site expression is governed by the surrounding splicing landscape. The usage of alternative polyA sites in terminal exons is modulated by their relative locations and distance to downstream genes. Finally, we apply our models to reveal thousands of disease- and trait-associated genetic variants altering polyadenylation activity. Altogether, our models represent a valuable resource to dissect molecular mechanisms mediating genome-wide polyA site expression and characterize their functional roles in human diseases.

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.

Pre-mRNA splicing and its cotranscriptional connections

Transcription of eukaryotic genes by RNA polymerase II (Pol II) yields RNA precursors containing introns that must be spliced out and the flanking exons ligated together. Splicing is catalyzed by a dynamic ribonucleoprotein complex called the spliceosome. Recent evidence has shown that a large fraction of splicing occurs cotranscriptionally as the RNA chain is extruded from Pol II at speeds of up to 5 kb/minute. Splicing is more efficient when it is tethered to the transcription elongation complex, and this linkage permits functional coupling of splicing with transcription. We discuss recent progress that has uncovered a network of connections that link splicing to transcript elongation and other cotranscriptional RNA processing events.

Pre-mRNA splicing order is predetermined and maintains splicing fidelity across multi-intronic transcripts

Combinatorially, intron excision within a given nascent transcript could proceed down any of thousands of paths, each of which would expose different dynamic landscapes of cis-elements and contribute to alternative splicing. In this study, we found that post-transcriptional multi-intron splicing order in human cells is largely predetermined, with most genes spliced in one or a few predominant orders. Strikingly, these orders were conserved across cell types and stages of motor neuron differentiation. Introns flanking alternatively spliced exons were frequently excised last, after their neighboring introns. Perturbations to the spliceosomal U2 snRNA altered the preferred splicing order of many genes, and these alterations were associated with the retention of other introns in the same transcript. In one gene, early removal of specific introns was sufficient to induce delayed excision of three proximal introns, and this delay was caused by two distinct cis-regulatory mechanisms. Together, our results demonstrate that multi-intron splicing order in human cells is predetermined, is influenced by a component of the spliceosome and ensures splicing fidelity across long pre-mRNAs.

Identification of neoantigens and immunological subtypes in clear cell renal cell carcinoma for mRNA vaccine development and patient selection

Clear cell renal cell carcinoma (ccRCC) is a common urological malignancy with diverse histological types. This study aimed to detect neoantigens in ccRCC to develop mRNA vaccines and distinguish between ccRCC immunological subtypes for construction of an immune landscape to select patients suitable for vaccination. Using The Cancer Genome Atlas SpliceSeq database, The Cancer Genome Atlas, and the International Cancer Genome Consortium cohorts, we comprehensively analysed potential tumour antigens of ccRCC associated with aberrant alternative splicing, somatic mutation, nonsense-mediated mRNA decay factors, antigen-presenting cells, and overall survival. Immune subtypes (C1/C2) and nine immune gene modules of ccRCC were identified by consistency clustering and weighted correlation network analysis. The immune landscape as well as molecular and cellular characteristics of immunotypes were assessed. Rho-guanine nucleotide exchange factor 3 (ARHGEF3) was identified as a new ccRCC antigen for development of an mRNA vaccine. A higher tumour mutation burden, differential expression of immune checkpoints, and immunogenic cell death were observed in cases with the C2 immunotype. Cellular characteristics increased the complexity of the immune environment, and worse outcomes were observed in ccRCC cases with the C2 immunotype. We constructed the immune landscape for selecting patients with the C2 immunotype suitable for vaccination.

Potential role of N6-adenine DNA methylation in alternative splicing and endosymbiosis in Paramecium bursaria

N⁶-adenine DNA methylation (6mA), a rediscovered epigenetic mark in eukaryotic organisms, diversifies in abundance, distribution, and function across species, necessitating its study in more taxa. Paramecium bursaria is a typical model organism with endosymbiotic algae of the species Chlorella variabilis. This consortium therefore serves as a valuable system to investigate the functional role of 6mA in endosymbiosis, as well as the evolutionary importance of 6mA among eukaryotes. In this study, we report the first genome-wide, base pair-resolution map of 6mA in P. bursaria and identify its methyltransferase PbAMT1. Functionally, 6mA exhibits a bimodal distribution at the 5' end of RNA polymerase II-transcribed genes and possibly participates in transcription by facilitating alternative splicing. Evolutionarily, 6mA co-evolves with gene age and likely serves as a reverse mark of endosymbiosis-related genes. Our results offer new insights for the functional diversification of 6mA in eukaryotes as an important epigenetic mark.

Introns: the "dark matter" of the eukaryotic genome

The emergence of introns was a significant evolutionary leap that is a major distinguishing feature between prokaryotic and eukaryotic genomes. While historically introns were regarded merely as the sequences that are removed to produce spliced transcripts encoding functional products, increasingly data suggests that introns play important roles in the regulation of gene expression. Here, we use an intron-centric lens to review the role of introns in eukaryotic gene expression. First, we focus on intron architecture and how it may influence mechanisms of splicing. Second, we focus on the implications of spliceosomal snRNAs and their variants on intron splicing. Finally, we discuss how the presence of introns and the need to splice them influences transcription regulation. Despite the abundance of introns in the eukaryotic genome and their emerging role regulating gene expression, a lot remains unexplored. Therefore, here we refer to introns as the "dark matter" of the eukaryotic genome and discuss some of the outstanding questions in the field.

Structural basis of catalytic activation in human splicing

Pre-mRNA splicing follows a pathway driven by ATP-dependent RNA helicases. A crucial event of the splicing pathway is the catalytic activation, which takes place at the transition between the activated Bact and the branching-competent B* spliceosomes. Catalytic activation occurs through an ATP-dependent remodelling mediated by the helicase PRP2 (also known as DHX16)1-3. However, because PRP2 is observed only at the periphery of spliceosomes3-5, its function has remained elusive. Here we show that catalytic activation occurs in two ATP-dependent stages driven by two helicases: PRP2 and Aquarius. The role of Aquarius in splicing has been enigmatic6,7. Here the inactivation of Aquarius leads to the stalling of a spliceosome intermediate-the BAQR complex-found halfway through the catalytic activation process. The cryogenic electron microscopy structure of BAQR reveals how PRP2 and Aquarius remodel Bact and BAQR, respectively. Notably, PRP2 translocates along the intron while it strips away the RES complex, opens the SF3B1 clamp and unfastens the branch helix. Translocation terminates six nucleotides downstream of the branch site through an assembly of PPIL4, SKIP and the amino-terminal domain of PRP2. Finally, Aquarius enables the dissociation of PRP2, plus the SF3A and SF3B complexes, which promotes the relocation of the branch duplex for catalysis. This work elucidates catalytic activation in human splicing, reveals how a DEAH helicase operates and provides a paradigm for how helicases can coordinate their activities.

Rescue of blood coagulation Factor VIII exon-16 mis-splicing by antisense oligonucleotides

The human Factor VIII ( F8 ) protein is essential for the blood coagulation cascade and specific F8 mutations cause the rare bleeding disorder Hemophilia A (HA). Here, we investigated the impact of HA-causing single-nucleotide mutations on F8 pre-mRNA splicing. We found that 14/97 (∼14.4%) coding sequence mutations tested in our study induced exon skipping. Splicing patterns of 4/11 (∼36.4%) F8 exons tested were especially sensitive to the presence of common disease-causing mutations. RNA-chemical probing analyses revealed a three-way junction structure at the 3' end of intron 15 (TWJ-3-15). TWJ-3-15 sequesters the polypyrimidine tract, a key determinant of 3' splice site strength. Using exon-16 of the F8 gene as a model, we designed specific antisense oligonucleotides (ASOs) that target TWJ-3-15 and identified three that promote the splicing of F8 exon-16. Interaction of TWJ-3-15 with ASOs increases accessibility of the polypyrimidine tract and inhibits the binding of hnRNPA1-dependent splicing silencing factors. Moreover, ASOs targeting TWJ-3-15 rescue diverse splicing-sensitive HA-causing mutations, most of which are distal to the 3' splice site being impacted. The TWJ-3-15 structure and its effect on mRNA splicing provide a model for HA etiology in patients harboring specific F8 mutations and provide a framework for precision RNA-based HA therapies.

Regulation of pre-mRNA splicing: roles in physiology and disease, and therapeutic prospects

The removal of introns from mRNA precursors and its regulation by alternative splicing are key for eukaryotic gene expression and cellular function, as evidenced by the numerous pathologies induced or modified by splicing alterations. Major recent advances have been made in understanding the structures and functions of the splicing machinery, in the description and classification of physiological and pathological isoforms and in the development of the first therapies for genetic diseases based on modulation of splicing. Here, we review this progress and discuss important remaining challenges, including predicting splice sites from genomic sequences, understanding the variety of molecular mechanisms and logic of splicing regulation, and harnessing this knowledge for probing gene function and disease aetiology and for the design of novel therapeutic approaches.

Exon definitive regions for MPC1 microexon splicing and its usage for splicing modulation

Alternative splicing of microexons (3-30 base pairs [bp]) is involved in important biological processes in brain development and human cancers. However, understanding a splicing process of non-3x bp microexons is scarce. We showed that 4 bp microexon of mitochondrial pyruvate carrier1 (MPC1) is constitutively included in mRNA. Based on our studies with minigene and exon island constructs, we found the strong exon definition region in the proximal introns bordering MPC1 microexon. Ultimately, we defined a nucleotide fragment from the 3'ss 67 bp of MPC1 microexon to the 5'ss consensus sequence, as a core exon island, which can concatenate its microexon and neighboring exons by splicing. Furthermore, we showed that insertion of the core exon island into a target exon or intron induced skip the target exon or enhance the splicing of an adjacent exon, respectively. Collectively, we suggest that the exon island derived from MPC1 microexon modifies genuine splicing patterns depending on its position, thereby providing insights on strategies for splicing-mediated gene correction.

Position-dependent effects of RNA-binding proteins in the context of co-transcriptional splicing

Alternative splicing is an important step in eukaryotic mRNA pre-processing which increases the complexity of gene expression programs, but is frequently altered in disease. Previous work on the regulation of alternative splicing has demonstrated that splicing is controlled by RNA-binding proteins (RBPs) and by epigenetic DNA/histone modifications which affect splicing by changing the speed of polymerase-mediated pre-mRNA transcription. The interplay of these different layers of splicing regulation is poorly understood. In this paper, we derived mathematical models describing how splicing decisions in a three-exon gene are made by combinatorial spliceosome binding to splice sites during ongoing transcription. We additionally take into account the effect of a regulatory RBP and find that the RBP binding position within the sequence is a key determinant of how RNA polymerase velocity affects splicing. Based on these results, we explain paradoxical observations in the experimental literature and further derive rules explaining why the same RBP can act as inhibitor or activator of cassette exon inclusion depending on its binding position. Finally, we derive a stochastic description of co-transcriptional splicing regulation at the single-cell level and show that splicing outcomes show little noise and follow a binomial distribution despite complex regulation by a multitude of factors. Taken together, our simulations demonstrate the robustness of splicing outcomes and reveal that quantitative insights into kinetic competition of co-transcriptional events are required to fully understand this important mechanism of gene expression diversity.

Roles of minor spliceosome in intron recognition and the convergence with the better understood major spliceosome

Catalyzed by spliceosomes in the nucleus, RNA splicing removes intronic sequences from precursor RNAs in eukaryotes to generate mature RNA, which also significantly increases proteome complexity and fine-tunes gene expression. Most metazoans have two coexisting spliceosomes; the major spliceosome, which removes >99.5% of introns, and the minor spliceosome, which removes far fewer introns (only 770 at present have been predicted in the human genome). Both spliceosomes are large and dynamic machineries, each consisting of five small nuclear RNAs (snRNAs) and more than 100 proteins. However, the dynamic assembly, catalysis, and protein composition of the minor spliceosome are still poorly understood. With different splicing signals, minor introns are rare and usually distributed alone and flanked by major introns in genes, raising questions of how they are recognized by the minor spliceosome and how their processing deals with the splicing of neighboring major introns. Due to large numbers of introns and close similarities between the two machinery, cooperative, and competitive recognition by the two spliceosomes has been investigated. Functionally, many minor-intron-containing genes are evolutionarily conserved and essential. Mutations in the minor spliceosome exhibit a variety of developmental defects in plants and animals and are linked to numerous human diseases. Here, we review recent progress in the understanding of minor splicing, compare currently known components of the two spliceosomes, survey minor introns in a wide range of organisms, discuss cooperation and competition of the two spliceosomes in splicing of minor-intron-containing genes, and contributions of minor splicing mutations in development and diseases. This article is categorized under: RNA Processing > Processing of Small RNAs RNA Processing > Splicing Mechanisms RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry.

Profiling lariat intermediates reveals genetic determinants of early and late co-transcriptional splicing

Long introns with short exons in vertebrate genes are thought to require spliceosome assembly across exons (exon definition), rather than introns, thereby requiring transcription of an exon to splice an upstream intron. Here, we developed CoLa-seq (co-transcriptional lariat sequencing) to investigate the timing and determinants of co-transcriptional splicing genome wide. Unexpectedly, 90% of all introns, including long introns, can splice before transcription of a downstream exon, indicating that exon definition is not obligatory for most human introns. Still, splicing timing varies dramatically across introns, and various genetic elements determine this variation. Strong U2AF2 binding to the polypyrimidine tract predicts early splicing, explaining exon definition-independent splicing. Together, our findings question the essentiality of exon definition and reveal features beyond intron and exon length that are determinative for splicing timing.

Structural studies of the spliceosome: Bridging the gaps

The spliceosome is a multi-megadalton RNA-protein complex responsible for the removal of non-coding introns from pre-mRNAs. Due to its complexity and dynamic nature, it has proven to be a very challenging target for structural studies. Developments in single particle cryo-EM have overcome these previous limitations and paved the way towards a structural characterisation of the splicing machinery. Despite tremendous progress, many aspects of spliceosome structure and function remain elusive. In particular, the events leading to the definition of exon-intron boundaries, alternative and non-canonical splicing events, and cross-talk with other cellular machineries. Efforts are being made to address these knowledge gaps and further our mechanistic understanding of the spliceosome. Here, we summarise recent progress in the structural and functional analysis of the spliceosome.

Human PRPF39 is an alternative splicing factor recruiting U1 snRNP to weak 5' splice sites

Human PRPF39 is a homolog of the yeast Prp39 and Prp42 paralogs. We have previously shown that human PRPF39 forms a homodimer that interacts with the CTD of U1C, mirroring the yeast Prp39/Prp42 heterodimer. We demonstrate here that PRPF39 knockdown in HEK293 cells affects many alternative splicing events primarily by reducing the usage of weak 5'ss. Additionally, PRPF39 preferentially binds to a GC-rich RNA, likely at the interface between its NTD and CTD. These data indicate that PRPF39 potentially recruits U1 snRNP to a weak 5' ss, serving as a previously unrecognized alternative splicing factor. We further demonstrate that human TIA1 binds to U1C through its RRM1 and RRM3+Q domains but has no significant binding to PRPF39. Finally, all three human LUC7L isoforms directly interact with U1C. These results reveal significant parallels to the yeast U1 snRNP structure and support the use of yeast U1 snRNP as a model for understanding the mechanism of human alternative splicing.

Transcriptome-Wide Detection of Intron/Exon Definition in the Endogenous Pre-mRNA Transcripts of Mammalian Cells and Its Regulation by Depolarization

Pairing of splice sites across an intron or exon is the central point of intron or exon definition in pre-mRNA splicing with the latter mode proposed for most mammalian exons. However, transcriptome-wide pairing within endogenous transcripts has not been examined for the prevalence of each mode in mammalian cells. Here we report such pairings in rat GH3 pituitary cells by measuring the relative abundance of nuclear RNA-Seq reads at the intron start or end (RISE). Interestingly, RISE indexes are positively correlated between 5′ and 3′ splice sites specifically across introns or exons but inversely correlated with the usage of adjacent exons. Moreover, the ratios between the paired indexes were globally modulated by depolarization, which was disruptible by 5-aza-Cytidine. The nucleotide matrices of the RISE-positive splice sites deviate significantly from the rat consensus, and short introns or exons are enriched with the cross-intron or -exon RISE pairs, respectively. Functionally, the RISE-positive genes cluster for basic cellular processes including RNA binding/splicing, or more specifically, hormone production if regulated by depolarization. Together, the RISE analysis identified the transcriptome-wide regulation of either intron or exon definition between weak splice sites of short introns/exons in mammalian cells. The analysis also provides a way to further track the splicing intermediates and intron/exon definition during the dynamic regulation of alternative splicing by extracellular factors.

Structural advances in transcription elongation

Transcription is the first step of gene expression and involves RNA polymerases. After transcription initiation, RNA polymerase enters elongation followed by transcription termination at the end of the gene. Only recently, structures of transcription elongation complexes bound to key transcription elongation factors have been determined in bacterial and eukaryotic systems. These structures have revealed numerous insights including the basis for transcriptional pausing, RNA polymerase interaction with large complexes such as the ribosome and the spliceosome, and the transition into productive elongation. Here, we review these structures and describe areas for future research.

Alternative splicing as a source of phenotypic diversity

A major goal of evolutionary genetics is to understand the genetic processes that give rise to phenotypic diversity in multicellular organisms. Alternative splicing generates multiple transcripts from a single gene, enriching the diversity of proteins and phenotypic traits. It is well established that alternative splicing contributes to key innovations over long evolutionary timescales, such as brain development in bilaterians. However, recent developments in long-read sequencing and the generation of high-quality genome assemblies for diverse organisms has facilitated comparisons of splicing profiles between closely related species, providing insights into how alternative splicing evolves over shorter timescales. Although most splicing variants are probably non-functional, alternative splicing is nonetheless emerging as a dynamic, evolutionarily labile process that can facilitate adaptation and contribute to species divergence.

Exonic splicing code and protein binding sites for calcium

Auxilliary splicing sequences in exons, known as enhancers (ESEs) and silencers (ESSs), have been subject to strong selection pressures at the RNA and protein level. The protein component of this splicing code is substantial, recently estimated at ∼50% of the total information within ESEs, but remains poorly understood. The ESE/ESS profiles were previously associated with the Irving-Williams (I-W) stability series for divalent metals, suggesting that the ESE/ESS evolution was shaped by metal binding sites. Here, we have examined splicing activities of exonic sequences that encode protein binding sites for Ca2+, a weak binder in the I-W affinity order. We found that predicted exon inclusion levels for the EF-hand motifs and for Ca2+-binding residues in nonEF-hand proteins were higher than for average exons. For canonical EF-hands, the increase was centred on the EF-hand chelation loop and, in particular, on Ca2+-coordinating residues, with a 1>12>3∼5>9 hierarchy in the 12-codon loop consensus and usage bias at codons 1 and 12. The same hierarchy but a lower increase was observed for noncanonical EF-hands, except for S100 proteins. EF-hand loops preferentially accumulated exon splits in two clusters, one located in their N-terminal halves and the other around codon 12. Using splicing assays and published crosslinking and immunoprecipitation data, we identify candidate trans-acting factors that preferentially bind conserved GA-rich motifs encoding negatively charged amino acids in the loops. Together, these data provide evidence for the high capacity of codons for Ca2+-coordinating residues to be retained in mature transcripts, facilitating their exon-level expansion during eukaryotic evolution.

Gene architecture directs splicing outcome in separate nuclear spatial regions

How the splicing machinery defines exons or introns as the spliced unit has remained a puzzle for 30 years. Here, we demonstrate that peripheral and central regions of the nucleus harbor genes with two distinct exon-intron GC content architectures that differ in the splicing outcome. Genes with low GC content exons, flanked by long introns with lower GC content, are localized in the periphery, and the exons are defined as the spliced unit. Alternative splicing of these genes results in exon skipping. In contrast, the nuclear center contains genes with a high GC content in the exons and short flanking introns. Most splicing of these genes occurs via intron definition, and aberrant splicing leads to intron retention. We demonstrate that the nuclear periphery and center generate different environments for the regulation of alternative splicing and that two sets of splicing factors form discrete regulatory subnetworks for the two gene architectures. Our study connects 3D genome organization and splicing, thus demonstrating that exon and intron definition modes of splicing occur in different nuclear regions.

Transcription and splicing dynamics during early Drosophila development

Widespread cotranscriptional splicing has been demonstrated from yeast to human. However, most studies to date addressing the kinetics of splicing relative to transcription used either Saccharomyces cerevisiae or metazoan cultured cell lines. Here, we adapted native elongating transcript sequencing technology (NET-seq) to measure cotranscriptional splicing dynamics during the early developmental stages of Drosophila melanogaster embryos. Our results reveal the position of RNA polymerase II (Pol II) when both canonical and recursive splicing occur. We found heterogeneity in splicing dynamics, with some RNAs spliced immediately after intron transcription, whereas for other transcripts no splicing was observed over the first 100 nt of the downstream exon. Introns that show splicing completion before Pol II has reached the end of the downstream exon are necessarily intron-defined. We studied the splicing dynamics of both nascent pre-mRNAs transcribed in the early embryo, which have few and short introns, as well as pre-mRNAs transcribed later in embryonic development, which contain multiple long introns. As expected, we found a relationship between the proportion of spliced reads and intron size. However, intron definition was observed at all intron sizes. We further observed that genes transcribed in the early embryo tend to be isolated in the genome whereas genes transcribed later are often overlapped by a neighboring convergent gene. In isolated genes, transcription termination occurred soon after the polyadenylation site, while in overlapped genes, Pol II persisted associated with the DNA template after cleavage and polyadenylation of the nascent transcript. Taken together, our data unravel novel dynamic features of Pol II transcription and splicing in the developing Drosophila embryo.

Genetic analysis of 76 Spanish Pompe disease patients: Identification of 12 novel pathogenic GAA variants and functional characterization of splicing variants

Glycogenosis type II (GSDII), or Pompe disease (MIM 232300), is an inherited autosomal recessive disorder caused by deficiency of the lysosomal acid-α-glucosidase. Mutations in the GAA gene alter normal enzyme production and lead to progressive buildup of intralysosomal glycogen, which plays an essential role in the severity and progression of the disease. We report here the study of 76 patients from Spain with either infantile or late onset form of Pompe disease. The analysis consisted in the molecular study of exons and intron flanking fragments of GAA gene. We have identified 55 different molecular pathogenic variants, 12 of them not previously described. In addition, we have determined a frequency of 84.37% for the c.-32-13T>G mutation in patients with the late-onset form of the disease. Functional characterization of some splice mutations showed deleterious mechanisms on the processing of mRNA.

Splice site proximity influences alternative exon definition

Alternative splicing enables higher eukaryotes to expand mRNA diversity from a finite number of genes through highly combinatorial splice site selection mechanisms that are influenced by the sequence of competing splice sites, cis-regulatory elements binding trans-acting factors, the length of exons and introns harbouring alternative splice sites and RNA secondary structures at putative splice junctions. To test the hypothesis that the intron definition or exon definition modes of splice site recognition direct the selection of alternative splice patterns, we created a database of alternative splice site usage (ALTssDB). When alternative splice sites are embedded within short introns (intron definition), the 5' and 3' splice sites closest to each other across the intron preferentially pair, consistent with previous observations. However, when alternative splice sites are embedded within large flanking introns (exon definition), the 5' and 3' splice sites closest to each other across the exon are preferentially selected. Thus, alternative splicing decisions are influenced by the intron and exon definition modes of splice site recognition. The results demonstrate that the spliceosome pairs splice sites that are closest in proximity within the unit of initial splice site selection.

Alternative Splicing in Human Biology and Disease

Alternative pre-mRNA splicing allows for the production of multiple mRNAs from an individual gene, which not only expands the protein-coding potential of the genome but also enables complex mechanisms for the post-transcriptional control of gene expression. Regulation of alternative splicing entails a combinatorial interplay between an abundance of trans-acting splicing factors, cis-acting regulatory sequence elements and their concerted effects on the core splicing machinery. Given the extent and biological significance of alternative splicing in humans, it is not surprising that aberrant splicing patterns can cause or contribute to a wide range of diseases. In this introductory chapter, we outline the mechanisms that govern alternative pre-mRNA splicing and its regulation and discuss how dysregulated splicing contributes to human diseases affecting the motor system and the brain.

Novel compound heterozygous mutations in the CYP4F22 gene in a patient with autosomal recessive congenital ichthyosis

Autosomal recessive congenital ichthyosis (ARCI) is a rare form of keratinization disorder of the skin, which can be caused by mutations in 14 ARCI genes. We present a rare case of ARCI that carried a novel null mutation and a novel splice site mutation in the CYP4F22 gene.

Detecting circRNA in purified spliceosomal P complex

Circular RNAs (circRNAs) generated from back-splicing of exons have been found in a wide range of eukaryotic species and exert a variety of biological functions. Unlike canonical splicing, the mechanism of back-splicing has long remained elusive. We recently determined the cryo-EM structure of the yeast spliceosomal E complex assembled on introns, leading us to hypothesize that the same E complex can assemble across an exon forming the exon-definition complex. This complex, when assembled on long exons, goes through the splicing cycle and catalyzes back-splicing to generate circRNAs. Supporting this hypothesis, we purified the yeast post-catalytic spliceosomal P complex (the best complex in the splicing cycle to trap splicing products and intermediates) and detected canonical and back-splicing products as well as splicing intermediates. Here we describe in detail this procedure, which may be applied to other organisms to facilitate research on the biogenesis and regulation of circRNA.

Evolution of the Early Spliceosomal Complex-From Constitutive to Regulated Splicing

Pre-mRNA splicing is a major process in the regulated expression of genes in eukaryotes, and alternative splicing is used to generate different proteins from the same coding gene. Splicing is a catalytic process that removes introns and ligates exons to create the RNA sequence that codifies the final protein. While this is achieved in an autocatalytic process in ancestral group II introns in prokaryotes, the spliceosome has evolved during eukaryogenesis to assist in this process and to finally provide the opportunity for intron-specific splicing. In the early stage of splicing, the RNA 5' and 3' splice sites must be brought within proximity to correctly assemble the active spliceosome and perform the excision and ligation reactions. The assembly of this first complex, termed E-complex, is currently the least understood process. We focused in this review on the formation of the E-complex and compared its composition and function in three different organisms. We highlight the common ancestral mechanisms in S. cerevisiae, S. pombe, and mammals and conclude with a unifying model for intron definition in constitutive and regulated co-transcriptional splicing.

An Exon Skipping in CRS1 Is Associated with Perturbed Chloroplast Development in Maize

Chloroplasts of higher plants are semi-autonomous organelles that perform photosynthesis and produce hormones and metabolites. They play crucial roles in plant growth and development. Although many seedling-lethal nuclear genes or regulators required for chloroplast development have been characterized, the understanding of chloroplast development is still limited. Using a genetic screen, we isolated a mutant named ell1, with etiolated leaves and a seedling-lethal phenotype. Analysis by BN-PAGE and transmission electron microscopy revealed drastic morphological defects of chloroplasts in ell1 mutants. Genetic mapping of the mutant gene revealed a single mutation (G-to-A) at the 5′ splice site of intron 5 in CRS1, resulting in an exon skipping in CRS1, indicating that this mutation in CRS1 is responsible for the observed phenotype, which was further confirmed by genetic analysis. The incorrectly spliced CRS1 failed to mediate the splicing of atpF intron. Moreover, the quantitative analysis suggested that ZmCRS1 may participate in chloroplast transcription to regulate the development of chloroplast. Taken together, these findings improve our understanding of the ZmCRS1 protein and shed new light on the regulation of chloroplast development in maize.

Wilms' Tumor 1-Associating Protein complex regulates alternative splicing and polyadenylation at potential G-quadruplex-forming splice site sequences

Wilms’ tumor 1-associating protein (WTAP) is a core component of the N6-methyladenosine (m6A)-methyltransferase complex, along with VIRMA, CBLL1, ZC3H13 (KIAA0853), RBM15/15B, and METTL3/14, which generate m6A, a key RNA modification that affects various processes of RNA metabolism. WTAP also interacts with splicing factors; however, despite strong evidence suggesting a role of Drosophila WTAP homolog fl(2)d in alternative splicing (AS), its role in splicing regulation in mammalian cells remains elusive. Here we demonstrate using RNAi coupled with RNA-seq that WTAP, VIRMA, CBLL1, and ZC3H13 modulate AS, promoting exon skipping and intron retention in AS events that involve short introns/exons with higher GC content and introns with weaker polypyrimidine-tract and branch points. Further analysis of GC-rich sequences involved in AS events regulated by WTAP, together with minigene assay analysis, revealed potential G-quadruplex formation at splice sites where WTAP has an inhibitory effect. We also found that several AS events occur in the last exon of one isoform of MSL1 and WTAP, leading to competition for polyadenylation. Proteomic analysis also suggested that WTAP/CBLL1 interaction promotes recruitment of the 3′-end processing complex. Taken together, our results indicate that the WTAP complex regulates AS and alternative polyadenylation via inhibitory mechanisms in GC-rich sequences.

Molecular Screening of PAX2 Gene Polymorphism in Primary Vesicoureteral Reflux Patients in Taif Governorate, KSA

<b>Background and Objective:</b> Primary Nonsyndromic Vesicoureteral Reflux (PVUR) is a widespread genetic malformation and considered a prevalent Congenital Abnormality of the Kidney and Urinary Tract (CAKUT). Mutations in the <i>PAX2 </i>gene have been associated with abnormalities in the kidney extending from CAKUT to oncogenic processes. The present study analyzes the <i>PAX2</i> polymorphisms and their association with primary VUR in Saudi children patients from the Taif governorate. <b>Materials and Methods:</b> Fifteen children with primary VUR were identified and screened for gene mutations in the <i>PAX2</i> gene by direct sequencing method of purified Polymerase Chain Reaction (PCR) products of all exons to elucidate the correlation between <i>PAX2</i> gene and VUR. <b>Results:</b> Seven new variants have been defined. Three polymorphic missense variants in homozygous genotype form were found in intron 8 and detected in eight patients, One missense mutation was found in exon 10 in the site of transactivation domain and detected in ten patients and <i>in-silico</i> analysis predicted it as a pathogenic one, Three mutations were found in exon 11 and detected in all patients as a compound homozygous. <b>Conclusion:</b> <i>PAX2</i>is important for normal kidney development and mutations in the gene possibly lead to disturbance in the protein structure and could be non-functional thus mutations in <i>PAX2</i> may be one of the causes of PVUR in Saudi Arabia. Further investigation is necessary to understand the aetiology of disease and maybe other genes implicated in VUR.