Mechanisms for U2AF to define 3' splice sites and regulate alternative splicing in the human genome

Representative cancer-associated U2AF2 mutations alter RNA interactions and splicing

High-throughput sequencing of hematologic malignancies and other cancers has revealed recurrent mis-sense mutations of genes encoding pre-mRNA splicing factors. The essential splicing factor U2AF2 recognizes a polypyrimidine-tract splice-site signal and initiates spliceosome assembly. Here, we investigate representative, acquired U2AF2 mutations, namely N196K or G301D amino acid substitutions associated with leukemia or solid tumors, respectively. We determined crystal structures of the wild-type (WT) compared with N196K- or G301D-substituted U2AF2 proteins, each bound to a prototypical AdML polypyrimidine tract, at 1.5, 1.4, or 1.7 Å resolutions. The N196K residue appears to stabilize the open conformation of U2AF2 with an inter-RNA recognition motif hydrogen bond, in agreement with an increased apparent RNA-binding affinity of the N196K-substituted protein. The G301D residue remains in a similar position as the WT residue, where unfavorable proximity to the RNA phosphodiester could explain the decreased RNA-binding affinity of the G301D-substituted protein. We found that expression of the G301D-substituted U2AF2 protein reduces splicing of a minigene transcript carrying prototypical splice sites. We further show that expression of either N196K- or G301D-substituted U2AF2 can subtly alter splicing of representative endogenous transcripts, despite the presence of endogenous, WT U2AF2 such as would be present in cancer cells. Altogether, our results demonstrate that acquired U2AF2 mutations such as N196K and G301D are capable of dysregulating gene expression for neoplastic transformation.

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

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

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.

Complex landscape of alternative splicing in myeloid neoplasms

Myeloid neoplasms are characterized by frequent mutations in at least seven components of the spliceosome that have distinct roles in the process of pre-mRNA splicing. Hotspot mutations in SF3B1, SRSF2, U2AF1 and loss of function mutations in ZRSR2 have revealed widely different aberrant splicing signatures with little overlap. However, previous studies lacked the power necessary to identify commonly mis-spliced transcripts in heterogeneous patient cohorts. By performing RNA-Seq on bone marrow samples from 1,258 myeloid neoplasm patients and 63 healthy bone marrow donors, we identified transcripts frequently mis-spliced by mutated splicing factors (SF), rare SF mutations with common alternative splicing (AS) signatures, and SF-dependent neojunctions. We characterized 17,300 dysregulated AS events using a pipeline designed to predict the impact of mis-splicing on protein function. Meta-splicing analysis revealed a pattern of reduced levels of retained introns among disease samples that was exacerbated in patients with splicing factor mutations. These introns share characteristics with “detained introns,” a class of introns that have been shown to promote differentiation by detaining pro-proliferative transcripts in the nucleus. In this study, we have functionally characterized 17,300 targets of mis-splicing by the SF mutations, identifying a common pathway by which AS may promote maintenance of a proliferative state.

Spliceosomal factor mutations and mis-splicing in MDS

Somatic, heterozygous missense and nonsense mutations in at least seven proteins that function in the spliceosome are found at high frequency in MDS patients. These proteins act at various steps in the process of splicing by the spliceosome and lead to characteristic alterations in the alternative splicing of a subset of genes. Several studies have investigated the effects of these mutations and have attempted to identify a commonly affected gene or pathway. Here, we summarize what is known about the normal function of these proteins and how the mutations alter the splicing landscape of the genome. We also summarize the commonly mis-spliced gene targets and discuss the state of mechanistic unification that has been achieved. Finally, we discuss alternative mechanisms by which these mutations may lead to disease.

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.

An autoinhibitory intramolecular interaction proof-reads RNA recognition by the essential splicing factor U2AF2

The recognition of cis-regulatory RNA motifs in human transcripts by RNA binding proteins (RBPs) is essential for gene regulation. The molecular features that determine RBP specificity are often poorly understood. Here, we combined NMR structural biology with high-throughput iCLIP approaches to identify a regulatory mechanism for U2AF2 RNA recognition. We found that the intrinsically disordered linker region connecting the two RNA recognition motif (RRM) domains of U2AF2 mediates autoinhibitory intramolecular interactions to reduce nonproductive binding to weak Py-tract RNAs. This proofreading favors binding of U2AF2 at stronger Py-tracts, as required to define 3' splice sites at early stages of spliceosome assembly. Mutations that impair the linker autoinhibition enhance the affinity for weak Py-tracts result in promiscuous binding of U2AF2 along mRNAs and impact on splicing fidelity. Our findings highlight an important role of intrinsically disordered linkers to modulate RNA interactions of multidomain RBPs.

The zinc finger domains in U2AF26 and U2AF35 have diverse functionalities including a role in controlling translation

Recent work has associated point mutations in both zinc fingers (ZnF) of the spliceosome component U2AF35 with malignant transformation. However, surprisingly little is known about the functionality of the U2AF35 ZnF domains in general. Here we have analysed key functionalities of the ZnF domains of mammalian U2AF35 and its paralog U2AF26. Both ZnFs are required for splicing regulation, whereas only ZnF2 controls protein stability and contributes to the interaction with U2AF65. These features are confirmed in a naturally occurring splice variant of U2AF26 lacking ZnF2, that is strongly induced upon activation of primary mouse T cells and localized in the cytoplasm. Using Ribo-Seq in a model T cell line we provide evidence for a role of U2AF26 in activating cytoplasmic steps in gene expression, notably translation. Consistently, an MS2 tethering assay shows that cytoplasmic U2AF26/35 increase translation when localized to the 5'UTR of a model mRNA. This regulation is partially dependent on ZnF1 thus providing a connection between a core splicing factor, the ZnF domains and the regulation of translation. Altogether, our work reveals unexpected functions of U2AF26/35 and their ZnF domains, thereby contributing to a better understanding of their role and regulation in mammalian cells.

Alternative Splicing Regulatory Networks: Functions, Mechanisms, and Evolution

High-throughput sequencing-based methods and their applications in the study of transcriptomes have revolutionized our understanding of alternative splicing. Networks of functionally coordinated and biologically important alternative splicing events continue to be discovered in an ever-increasing diversity of cell types in the context of physiologically normal and disease states. These studies have been complemented by efforts directed at defining sequence codes governing splicing and their cognate trans-acting factors, which have illuminated important combinatorial principles of regulation. Additional studies have revealed critical roles of position-dependent, multivalent protein-RNA interactions that direct splicing outcomes. Investigations of evolutionary changes in RNA binding proteins, splice variants, and associated cis elements have further shed light on the emergence, mechanisms, and functions of splicing networks. Progress in these areas has emphasized the need for a coordinated, community-based effort to systematically address the functions of individual splice variants associated with normal and disease biology.

Regulation of flowering transition by alternative splicing: the role of the U2 auxiliary factor

Flowering transition is regulated by complex genetic networks in response to endogenous and environmental signals. Pre-mRNA splicing is an essential step for the post-transcriptional regulation of gene expression. Alternative splicing of key flowering genes has been investigated in detail over the past decade. However, few splicing factors have been identified as being involved in flowering transition. Human heterodimeric splicing factor U2 snRNP auxiliary factor (U2AF) consists of two subunits, U2AF35 and U2AF65, and functions in 3' splice site recognition in mRNA splicing. Recent studies reveal that Arabidopsis U2AF65a/b and U2AF35a/b play important roles in the splicing of key flowering genes. We summarize recent advances in research on splicing-regulated flowering transition by focusing on the role of Arabidopsis U2AF in the splicing of key flowering-related genes at ambient temperature and in the abscisic acid signaling pathways.

Genetic basis of myelodysplastic syndromes

During the past decade, substantial progress has been made in the field of the genetics of myelodysplastic syndromes (MDS). These comprise a group of chronic myeloid neoplasms with abnormal cell morphology and progression to acute myeloid leukemia (AML), where revolutionary sequencing technologies have played a major role. Through extensive sequencing of a large number of MDS genomes, a comprehensive registry of driver mutations involved in the pathogenesis of MDS has been revealed, along with their impacts on clinical phenotype and prognosis. The most frequently affected molecules are involved in DNA methylations, chromatin modification, RNA splicing, transcription, signal transduction, cohesin regulation, and DNA repair. These mutations show strong positive and negative correlations with each other, suggesting the presence of functional interactions between mutations, which dictate disease progression. Because these mutations are associated with disease phenotype, drug response, and clinical outcomes, it is essential to be familiar with MDS genetics not only for better understanding of MDS pathogenesis but also for management of patients.

Functional significance of U2AF1 S34F mutations in lung adenocarcinomas

The functional role of U2AF1 mutations in lung adenocarcinomas (LUADs) remains incompletely understood. Here, we report a significant co-occurrence of U2AF1 S34F mutations with ROS1 translocations in LUADs. To characterize this interaction, we profiled effects of S34F on the transcriptome-wide distribution of RNA binding and alternative splicing in cells harboring the ROS1 translocation. Compared to its wild-type counterpart, U2AF1 S34F preferentially binds and modulates splicing of introns containing CAG trinucleotides at their 3' splice junctions. The presence of S34F caused a shift in cross-linking at 3' splice sites, which was significantly associated with alternative splicing of skipped exons. U2AF1 S34F induced expression of genes involved in the epithelial-mesenchymal transition (EMT) and increased tumor cell invasion. Finally, S34F increased splicing of the long over the short SLC34A2-ROS1 isoform, which was also associated with enhanced invasiveness. Taken together, our results suggest a mechanistic interaction between mutant U2AF1 and ROS1 in LUAD.

mTOR-regulated U2af1 tandem exon splicing specifies transcriptome features for translational control

U2 auxiliary factor 1 (U2AF1) functions in 3′-splice site selection during pre-mRNA processing. Alternative usage of duplicated tandem exons in U2AF1 produces two isoforms, U2AF1a and U2AF1b, but their functional differences are unappreciated due to their homology. Through integrative approaches of genome editing, customized-transcriptome profiling and crosslinking-mediated interactome analyses, we discovered that the expression of U2AF1 isoforms is controlled by mTOR and they exhibit a distinctive molecular profile for the splice site and protein interactomes. Mechanistic dissection of mutually exclusive alternative splicing events revealed that U2AF1 isoforms’ inherent differential preferences of nucleotide sequences and their stoichiometry determine the 3′-splice site. Importantly, U2AF1a-driven transcriptomes feature alternative splicing events in the 5′-untranslated region (5′-UTR) that are favorable for translation. These findings unveil distinct roles of duplicated tandem exon-derived U2AF1 isoforms in the regulation of the transcriptome and suggest U2AF1a-driven 5′-UTR alternative splicing as a molecular mechanism of mTOR-regulated translational control.

Characterizing the interplay between gene nucleotide composition bias and splicing

BACKGROUND

Nucleotide composition bias plays an important role in the 1D and 3D organization of the human genome. Here, we investigate the potential interplay between nucleotide composition bias and the regulation of exon recognition during splicing.

RESULTS

By analyzing dozens of RNA-seq datasets, we identify two groups of splicing factors that activate either about 3200 GC-rich exons or about 4000 AT-rich exons. We show that splicing factor-dependent GC-rich exons have predicted RNA secondary structures at 5' ss and are dependent on U1 snRNP-associated proteins. In contrast, splicing factor-dependent AT-rich exons have a large number of decoy branch points, SF1- or U2AF2-binding sites and are dependent on U2 snRNP-associated proteins. Nucleotide composition bias also influences local chromatin organization, with consequences for exon recognition during splicing. Interestingly, the GC content of exons correlates with that of their hosting genes, isochores, and topologically associated domains.

CONCLUSIONS

We propose that regional nucleotide composition bias over several dozens of kilobase pairs leaves a local footprint at the exon level and induces constraints during splicing that can be alleviated by local chromatin organization at the DNA level and recruitment of specific splicing factors at the RNA level. Therefore, nucleotide composition bias establishes a direct link between genome organization and local regulatory processes, like alternative splicing.

Cooperative evolution of two different TEs results in lineage-specific novel transcripts in the BLOC1S2 gene

BACKGROUND

The BLOC1S2 gene encodes the multifunctional protein BLOS2, a shared subunit of two lysosomal trafficking complexes: i) biogenesis of lysosome-related organelles complex-1 and i) BLOC-1-related complex. In our previous study, we identified an intriguing unreported transcript of the BLOC1S2 gene that has a novel exon derived from two transposable elements (TEs), MIR and AluSp. To investigate the evolutionary footprint and molecular mechanism of action of this transcript, we performed PCR and RT-PCR experiments and sequencing analyses using genomic DNA and RNA samples from humans and various non-human primates.

RESULTS

The results showed that the MIR element had integrated into the genome of our common ancestor, specifically in the BLOC1S2 gene region, before the radiation of all primate lineages and that the AluSp element had integrated into the genome of our common ancestor, fortunately in the middle of the MIR sequences, after the divergence of Old World monkeys and New World monkeys. The combined MIR and AluSp sequences provide a 3' splice site (AG) and 5' splice site (GT), respectively, and generate the Old World monkey-specific transcripts. Moreover, branch point sequences for the intron removal process are provided by the MIR and AluSp combination.

CONCLUSIONS

We show for the first time that sequential integration into the same location and sequence divergence events of two different TEs generated lineage-specific transcripts through sequence collaboration during primate evolution.

Promotion of adipogenesis by JMJD6 requires the AT hook-like domain and is independent of its catalytic function

JMJD6 is a member of the Jumonji C domain containing enzymes that demethylate and/or hydroxylate substrate proteins. It is a multi-functional protein that has been implicated in disparate aspects of transcriptional and post-transcriptional control of gene expression, including but not limited to enhancer and promoter binding, release of paused RNA polymerase II, control of splicing, and interaction with the translation machinery. JMJD6 contributes to multiple aspects of animal development, including adipogenesis modeled in culture. We mutated proposed or characterized domains in the JMJD6 protein to better understand the requirement for JMJD6 in adipogenic differentiation. Mutation of JMJD6 amino acids that mediate binding of iron and 2-oxogluterate, which are required cofactors for enzymatic activity, had no impact on JMJD6 function, showing that catalytic activity is not required for JMJD6 contributions to adipogenic differentiation. In addition, we documented the formation of JMJD6 oligomers and showed that catalytic activity is not required for oligomerization, as has been reported previously. We also observed no effect of mutations in the sumoylation site and in the poly-serine stretch. In contrast, mutation of the AT hook-like structure, which mediates interaction with DNA and/or RNA, compromised JMJD6 function by blocking its ability to interact with chromatin at genes that express regulators of adipogenesis. The ability of JMJD6 to interact with nucleic acids may be a critical requirement for its function in adipogenic differentiation. The requirement for the AT hook-like domain and the lack of requirement for catalytic activity giving rise to the idea that co-activation of transcription by JMJD6 may be functioning as a scaffold protein that supports the interactions of other critical regulators.

SF3B2-Mediated RNA Splicing Drives Human Prostate Cancer Progression

Androgen receptor splice variant-7 (AR-V7) is a constitutively active AR variant implicated in castration-resistant prostate cancers. Here, we show that the RNA splicing factor SF3B2, identified by in silico and CRISPR/Cas9 analyses, is a critical determinant of AR-V7 expression and is correlated with aggressive cancer phenotypes. Transcriptome and PAR-CLIP analyses revealed that SF3B2 controls the splicing of target genes, including AR, to drive aggressive phenotypes. SF3B2-mediated aggressive phenotypes in vivo were reversed by AR-V7 knockout. Pladienolide B, an inhibitor of a splicing modulator of the SF3b complex, suppressed the growth of tumors addicted to high SF3B2 expression. These findings support the idea that alteration of the splicing pattern by high SF3B2 expression is one mechanism underlying prostate cancer progression and therapeutic resistance. This study also provides evidence supporting SF3B2 as a candidate therapeutic target for treating patients with cancer. SIGNIFICANCE: RNA splicing factor SF3B2 is essential for the generation of an androgen receptor (AR) variant that renders prostate cancer cells resistant to AR-targeting therapy.Graphical Abstract: http://cancerres.aacrjournals.org/content/canres/79/20/5204/F1.large.jpg.

U2AF65 assemblies drive sequence-specific splice site recognition

The essential splicing factor U2AF65 is known to help anchoring U2 snRNP at the branch site. Its C-terminal UHM domain interacts with ULM motifs of SF3b155, an U2 snRNP protein. Here, we report a cooperative binding of U2AF65 and the related protein CAPERα to the multi-ULM domain of SF3b155. In addition, we show that the RS domain of U2AF65 drives a liquid-liquid phase separation that is amplified by intronic RNA with repeated pyrimidine tracts. In cells, knockdown of either U2AF65 or CAPERα improves the inclusion of cassette exons that are preceded by such repeated pyrimidine-rich motifs. These results support a model in which liquid-like assemblies of U2AF65 and CAPERα on repetitive pyrimidine-rich RNA sequences are driven by their RS domains, and facilitate the recruitment of the multi-ULM domain of SF3b155. We anticipate that posttranslational modifications and proteins recruited in dynamical U2AF65 and CAPERα condensates may further contribute to the complex mechanisms leading to specific splice site choice that occurs in cells.

PUF60-activated exons uncover altered 3' splice-site selection by germline missense mutations in a single RRM

PUF60 is a splicing factor that binds uridine (U)-rich tracts and facilitates association of the U2 small nuclear ribonucleoprotein with primary transcripts. PUF60 deficiency (PD) causes a developmental delay coupled with intellectual disability and spinal, cardiac, ocular and renal defects, but PD pathogenesis is not understood. Using RNA-Seq, we identify human PUF60-regulated exons and show that PUF60 preferentially acts as their activator. PUF60-activated internal exons are enriched for Us upstream of their 3′ splice sites (3′ss), are preceded by longer AG dinucleotide exclusion zones and more distant branch sites, with a higher probability of unpaired interactions across a typical branch site location as compared to control exons. In contrast, PUF60-repressed exons show U-depletion with lower estimates of RNA single-strandedness. We also describe PUF60-regulated, alternatively spliced isoforms encoding other U-bound splicing factors, including PUF60 partners, suggesting that they are co-regulated in the cell, and identify PUF60-regulated exons derived from transposed elements. PD-associated amino-acid substitutions, even within a single RNA recognition motif (RRM), altered selection of competing 3′ss and branch points of a PUF60-dependent exon and the 3′ss choice was also influenced by alternative splicing of PUF60. Finally, we propose that differential distribution of RNA processing steps detected in cells lacking PUF60 and the PUF60-paralog RBM39 is due to the RBM39 RS domain interactions. Together, these results provide new insights into regulation of exon usage by the 3′ss organization and reveal that germline mutation heterogeneity in RRMs can enhance phenotypic variability at the level of splice-site and branch-site selection.

AtU2AF65b functions in abscisic acid mediated flowering via regulating the precursor messenger RNA splicing of ABI5 and FLC in Arabidopsis

In mammalians and yeast, the splicing factor U2AF65/Mud2p functions in precursor messenger RNA (pre-mRNA) processing. Arabidopsis AtU2AF65b encodes a putative U2AF65 but its specific functions in plants are unknown. This paper examines the function of AtU2AF65b as a negative regulator of flowering time in Arabidopsis. We investigated the expression and function of AtU2AF65b in abscisic acid (ABA)-regulated flowering as well as the transcript abundance and pre-mRNA splicing of flowering-related genes in the knock-out mutants of AtU2AF65b. The atu2af65b mutants show early-flowering phenotype under both long-day and short-day conditions. The transcript accumulation of the flowering repressor gene FLOWERING LOCUS C (FLC) is reduced in the shoot apex of atu2af65b, due to both increased intron retention and reduced transcription activation. Reduced transcription of FLC results, at least partially, from the abnormal splicing and reduced transcript abundance of ABSCISIC ACID-INSENSITIVE 5 (ABI5), which encodes an activator of FLC in ABA-regulated flowering signaling. Additionally, the expression of AtU2AF65b is promoted by ABA. Transition to flowering and splicing of FLC and ABI5 in the atu2af65b mutants are compromised during ABA-induced flowering. ABA-responsive AtU2AF65b functions in the pre-mRNA splicing of FLC and ABI5 in shoot apex, whereby AtU2AF65b is involved in ABA-mediated flowering transition in Arabidopsis.

The RNA binding protein EWS is broadly involved in the regulation of pri-miRNA processing in mammalian cells

The Ewing Sarcoma protein (EWS) is a multifaceted RNA binding protein (RBP) with established roles in transcription, pre-mRNA processing and DNA damage response. By generating high quality EWS-RNA interactome, we uncovered its specific and prevalent interaction with a large subset of primary microRNAs (pri-miRNAs) in mammalian cells. Knockdown of EWS reduced, whereas overexpression enhanced, the expression of its target miRNAs. Biochemical analysis revealed that multiple elements in target pri-miRNAs, including the sequences flanking the stem-loop region, contributed to high affinity EWS binding and sequence swap experiments between target and non-target demonstrated that the flanking sequences provided the specificity for enhanced pri-miRNA processing by the Microprocessor Drosha/DGCR8. Interestingly, while repressing Drosha expression, as reported earlier, we found that EWS was able to enhance the recruitment of Drosha to chromatin. Together, these findings suggest that EWS may positively and negatively regulate miRNA biogenesis via distinct mechanisms, thus providing a new foundation to understand the function of EWS in development and disease.

Myelodysplastic Syndrome-Associated SRSF2 Mutations Cause Splicing Changes by Altering Binding Motif Sequences

Serine/arginine-rich splicing factor 2 (SRSF2) is a member of the SR protein family that is involved in both constitutive and alternative mRNA splicing. Mutations in SRSF2 gene are frequently reported in myelodysplastic syndromes (MDS) and acute myeloid leukemia (AML). It is imperative to understand how these mutations affect SRSF2-mediated splicing and cause MDS. In this study, we characterized MDS-associated SRSF2 mutants (P95H, P95L, and P95R). We found that those mutants and wild-type SRSF2 proteins showed nuclear localization in HeLa cells. In vitro splicing reaction also revealed that mutant proteins associated with both precursor and spliced mRNAs, suggesting that the mutants directly participate in splicing. We established the human myeloid leukemia K562 cell lines that stably expressed myc-tagged wild-type or mutant SRSF2 proteins, and then performed RNA-sequence to analyze the splicing pattern of each cell line. The results revealed that both wild-type and mutants affected splicing of approximately 3,000 genes. Although splice site sequences adjacent to the affected exons showed no significant difference compared to the total exons, exonic motif analyses with both inclusion- and exclusion-enhanced exons demonstrated that wild-type and mutants have different binding sequences in exons. These results indicate that mutations of SRSF2 in MDS change binding properties of SRSF2 to exonic motifs and this causes aberrant splicing.

Targeting an RNA-Binding Protein Network in Acute Myeloid Leukemia

RNA-binding proteins (RBPs) are essential modulators of transcription and translation frequently dysregulated in cancer. We systematically interrogated RBP dependencies in human cancers using a comprehensive CRISPR/Cas9 domain-focused screen targeting RNA-binding domains of 490 classical RBPs. This uncovered a network of physically interacting RBPs upregulated in acute myeloid leukemia (AML) and crucial for maintaining RNA splicing and AML survival. Genetic or pharmacologic targeting of one key member of this network, RBM39, repressed cassette exon inclusion and promoted intron retention within mRNAs encoding HOXA9 targets as well as in other RBPs preferentially required in AML. The effects of RBM39 loss on splicing further resulted in preferential lethality of spliceosomal mutant AML, providing a strategy for treatment of AML bearing RBP splicing mutations.

U2AF65 enhances milk synthesis and growth of bovine mammary epithelial cells by positively regulating the mTOR-SREBP-1c signalling pathway

U2 snRNP auxiliary factor 65 kDa (U2AF65) is a splicing factor that promotes prespliceosome assembly. The function of U2AF65 in alternative splicing has been identified; however, the essential physiological role of U2AF65 remains poorly understood. In this study, we investigated the regulatory role of U2AF65 in milk synthesis and growth of bovine mammary epithelial cells (BMECs). Our results showed that U2AF65 localizes in the nucleus. Treatment with amino acids (Met and Leu) and hormones (prolactin and β-estradiol) upregulated the expression of U2AF65 in these cells. U2AF65 overexpression increased the synthesis of β-casein, triglycerides, and lactose; increased cell viability; and promoted proliferation of BMECs. Furthermore, our results showed that U2AF65 positively regulated mTOR phosphorylation and expression of mature mRNA of mTOR and SREBP-1c. Collectively, our findings demonstrate that U2AF65 regulates the mRNA expression of signalling molecules (mTOR and SREBP-1c) involved in milk synthesis and growth of BMECs, possibly via controlling the splicing and maturation of these mRNAs. U2 snRNP auxiliary factor 65 kDa (U2AF65) is a splicing factor that promotes prespliceosome assembly. The essential physiological role of U2AF65 remains poorly understood. In the present study, we confirmed that U2AF65 functions as a positive regulator of milk synthesis in and proliferation of bovine mammary epithelial cells via the mTOR-SREBP-1c signalling pathway. Therefore, our study uncovers the regulatory role of U2AF65 in milk synthesis and cell proliferation.

PTBP1 acts as a dominant repressor of the aberrant tissue-specific splicing of ISCU in hereditary myopathy with lactic acidosis

BACKGROUND

Hereditary myopathy with lactic acidosis (HML) is an autosomal recessive disease caused by an intron mutation in the iron-sulfur cluster assembly (ISCU) gene. The mutation results in aberrant splicing, where part of the intron is retained in the final mRNA transcript, giving rise to a truncated nonfunctional ISCU protein. Using an ISCU mini-gene system, we have previously shown that PTBP1 can act as a repressor of the mis-splicing of ISCU, where overexpression of PTBP1 resulted in a decrease of the incorrect splicing. In this study, we wanted to, in more detail, analyze the role of PTBP1 in the regulation of endogenous ISCU mis-splicing.

METHODS

Overexpression and knockdown of PTBP1 was performed in myoblasts from two HML patients and a healthy control. Quantification of ISCU mis-splicing was done by qRTPCR. Biotinylated ISCU RNA, representing wildtype and mutant intron sequence, was used in a pull-down assay with nuclear extracts from myoblasts. Levels of PTBP1 in human cell lines and mice tissues were analyzed by qRTPCR and western blot.

RESULTS

PTBP1 overexpression in HML patient myoblasts resulted in a substantial decrease of ISCU mis-splicing while knockdown of PTBP1 resulted in a drastic increase. The effect could be observed in both patient and control myoblasts. We could also show that PTBP1 interacts with both the mutant and wild-type ISCU intron sequence, but with a higher affinity to the mutant sequence. Furthermore, low levels of PTBP1 among examined mouse tissues correlated with high levels of incorrect splicing of ISCU.

CONCLUSION

Our results show that PTBP1 acts as a dominant repressor of ISCU mis-splicing. We also show an inverse correlation between the levels of PTBP1 and ISCU mis-splicing, suggesting that the high level of mis-splicing in the skeletal muscle is primarily due to the low levels of PTBP1.

PRMT1 Deficiency in Mouse Juvenile Heart Induces Dilated Cardiomyopathy and Reveals Cryptic Alternative Splicing Products

Protein arginine methyltransferase 1 (PRMT1) catalyzes the asymmetric dimethylation of arginine residues in proteins and methylation of various RNA-binding proteins and is associated with alternative splicing in vitro. Although PRMT1 has essential in vivo roles in embryonic development, CNS development, and skeletal muscle regeneration, the functional importance of PRMT1 in the heart remains to be elucidated. Here, we report that juvenile cardiomyocyte-specific PRMT1-deficient mice develop severe dilated cardiomyopathy and exhibit aberrant cardiac alternative splicing. Furthermore, we identified previously undefined cardiac alternative splicing isoforms of four genes (Asb2, Fbxo40, Nrap, and Eif4a2) in PRMT1-cKO mice and revealed that eIF4A2 protein isoforms translated from alternatively spliced mRNA were differentially ubiquitinated and degraded by the ubiquitin-proteasome system. These findings highlight the essential roles of PRMT1 in cardiac homeostasis and alternative splicing regulation.

•

PRMT1 deficiency in cardiomyocytes causes dilated cardiomyopathy in juvenile mice

•

PRMT1-deficient heart shows abnormal alternative splicing patterns

•

Previously undefined cardiac splicing events are revealed by transcriptome analysis

•

eIF4A2 isoforms are differentially ubiquitinated and degraded

Aberrant splicing and defective mRNA production induced by somatic spliceosome mutations in myelodysplasia

Spliceosome mutations are frequently found in myelodysplasia. Splicing alterations induced by these mutations, their precise targets, and the effect at the transcript level have not been fully elucidated. Here we report transcriptomic analyses of 265 bone marrow samples from myelodysplasia patients, followed by a validation using CRISPR/Cas9-mediated gene editing and an assessment of nonsense-mediated decay susceptibility. Small but widespread reduction of intron-retaining isoforms is the most frequent splicing alteration in SF3B1-mutated samples. SF3B1 mutation is also associated with 3′ splice site alterations, leading to the most pronounced reduction of canonical transcripts. Target genes include tumor suppressors and genes of mitochondrial iron metabolism or heme biosynthesis. Alternative exon usage is predominant in SRSF2- and U2AF1-mutated samples. Usage of an EZH2 cryptic exon harboring a premature termination codon is increased in both SRSF2- and U2AF1-mutated samples. Our study reveals a landscape of splicing alterations and precise targets of various spliceosome mutations.

Mutations to the splicing machinery may have an important role in myelodysplasia. Here, the authors describe splicing factor gene mutations in myelodysplasia and report tumor suppressor, epigenetic, iron metabolism and heme biosynthesis genes as their targets.

An important class of intron retention events in human erythroblasts is regulated by cryptic exons proposed to function as splicing decoys

During terminal erythropoiesis, the splicing machinery in differentiating erythroblasts executes a robust intron retention (IR) program that impacts expression of hundreds of genes. We studied IR mechanisms in the SF3B1 splicing factor gene, which expresses ∼50% of its transcripts in late erythroblasts as a nuclear isoform that retains intron 4. RNA-seq analysis of nonsense-mediated decay (NMD)-inhibited cells revealed previously undescribed splice junctions, rare or not detected in normal cells, that connect constitutive exons 4 and 5 to highly conserved cryptic cassette exons within the intron. Minigene splicing reporter assays showed that these cassettes promote IR. Genome-wide analysis of splice junction reads demonstrated that cryptic noncoding cassettes are much more common in large (>1 kb) retained introns than they are in small retained introns or in nonretained introns. Functional assays showed that heterologous cassettes can promote retention of intron 4 in the SF3B1 splicing reporter. Although many of these cryptic exons were spliced inefficiently, they exhibited substantial binding of U2AF1 and U2AF2 adjacent to their splice acceptor sites. We propose that these exons function as decoys that engage the intron-terminal splice sites, thereby blocking cross-intron interactions required for excision. Developmental regulation of decoy function underlies a major component of the erythroblast IR program.

Spliceosome Mutations Induce R Loop-Associated Sensitivity to ATR Inhibition in Myelodysplastic Syndromes

Heterozygous somatic mutations in spliceosome genes (U2AF1, SF3B1, ZRSR2, or SRSF2) occur in >50% of patients with myelodysplastic syndrome (MDS). These mutations occur early in disease development, suggesting that they contribute to MDS pathogenesis and may represent a unique genetic vulnerability for targeted therapy. Here, we show that RNA splicing perturbation by expression of the U2AF1(S34F) mutant causes accumulation of R loops, a transcription intermediate containing RNA:DNA hybrids and displaced single-stranded DNA, and elicits an ATR response. ATR inhibitors (ATRi) induced DNA damage and cell death in U2AF1(S34F)-expressing cells, and these effects of ATRi were enhanced by splicing modulating compounds. Moreover, ATRi-induced DNA damage was suppressed by overexpression of RNaseH1, an enzyme that specifically removes the RNA in RNA:DNA hybrids, suggesting that the ATRi sensitivity of U2AF1(S34F)-expressing cells arises from R loops. Taken together, our results demonstrate that ATR may represent a novel therapeutic target in patients with MDS carrying the U2AF1(S34F) mutation and potentially other malignancies harboring spliceosome mutations.Significance: This study provides preclinical evidence that patients with MDS or other myeloid malignancies driven by spliceosome mutations may benefit from ATR inhibition to exploit the R loop-associated vulnerability induced by perturbations in splicing. Cancer Res; 78(18); 5363-74. ©2018 AACR.

Nuclear PTEN safeguards pre-mRNA splicing to link Golgi apparatus for its tumor suppressive role

Dysregulation of pre-mRNA alternative splicing (AS) is closely associated with cancers. However, the relationships between the AS and classic oncogenes/tumor suppressors are largely unknown. Here we show that the deletion of tumor suppressor PTEN alters pre-mRNA splicing in a phosphatase-independent manner, and identify 262 PTEN-regulated AS events in 293T cells by RNA sequencing, which are associated with significant worse outcome of cancer patients. Based on these findings, we report that nuclear PTEN interacts with the splicing machinery, spliceosome, to regulate its assembly and pre-mRNA splicing. We also identify a new exon 2b in GOLGA2 transcript and the exon exclusion contributes to PTEN knockdown-induced tumorigenesis by promoting dramatic Golgi extension and secretion, and PTEN depletion significantly sensitizes cancer cells to secretion inhibitors brefeldin A and golgicide A. Our results suggest that Golgi secretion inhibitors alone or in combination with PI3K/Akt kinase inhibitors may be therapeutically useful for PTEN-deficient cancers.

HNRNPA1 promotes recognition of splice site decoys by U2AF2 in vivo

Alternative pre-mRNA splicing plays a major role in expanding the transcript output of human genes. This process is regulated, in part, by the interplay of trans-acting RNA binding proteins (RBPs) with myriad cis-regulatory elements scattered throughout pre-mRNAs. These molecular recognition events are critical for defining the protein-coding sequences (exons) within pre-mRNAs and directing spliceosome assembly on noncoding regions (introns). One of the earliest events in this process is recognition of the 3′ splice site (3′ss) by U2 small nuclear RNA auxiliary factor 2 (U2AF2). Splicing regulators, such as the heterogeneous nuclear ribonucleoprotein A1 (HNRNPA1), influence spliceosome assembly both in vitro and in vivo, but their mechanisms of action remain poorly described on a global scale. HNRNPA1 also promotes proofreading of 3′ss sequences though a direct interaction with the U2AF heterodimer. To determine how HNRNPA1 regulates U2AF–RNA interactions in vivo, we analyzed U2AF2 RNA binding specificity using individual-nucleotide resolution crosslinking immunoprecipitation (iCLIP) in control and HNRNPA1 overexpression cells. We observed changes in the distribution of U2AF2 crosslinking sites relative to the 3′ss of alternative cassette exons but not constitutive exons upon HNRNPA1 overexpression. A subset of these events shows a concomitant increase of U2AF2 crosslinking at distal intronic regions, suggesting a shift of U2AF2 to “decoy” binding sites. Of the many noncanonical U2AF2 binding sites, Alu-derived RNA sequences represented one of the most abundant classes of HNRNPA1-dependent decoys. We propose that one way HNRNPA1 regulates exon definition is to modulate the interaction of U2AF2 with decoy or bona fide 3′ss.

In vitro iCLIP-based modeling uncovers how the splicing factor U2AF2 relies on regulation by cofactors

Alternative splicing generates distinct mRNA isoforms and is crucial for proteome diversity in eukaryotes. The RNA-binding protein (RBP) U2AF2 is central to splicing decisions, as it recognizes 3′ splice sites and recruits the spliceosome. We establish “in vitro iCLIP” experiments, in which recombinant RBPs are incubated with long transcripts, to study how U2AF2 recognizes RNA sequences and how this is modulated by trans-acting RBPs. We measure U2AF2 affinities at hundreds of binding sites and compare in vitro and in vivo binding landscapes by mathematical modeling. We find that trans-acting RBPs extensively regulate U2AF2 binding in vivo, including enhanced recruitment to 3′ splice sites and clearance of introns. Using machine learning, we identify and experimentally validate novel trans-acting RBPs (including FUBP1, CELF6, and PCBP1) that modulate U2AF2 binding and affect splicing outcomes. Our study offers a blueprint for the high-throughput characterization of in vitro mRNP assembly and in vivo splicing regulation.

Novel approach to functional SNPs discovery from genome-wide data reveals promising variants for colon cancer risk

In the majority of colorectal cancer (CRC) cases, the genetic basis of predisposition remains unexplained. The goal of the study was to assess the regulatory SNPs (rSNPs) in the human genome and to reveal СRC drivers based on the available chromatin immunoprecipitation sequencing (ChIP-Seq, ChIA-PET) and transcriptional profiling (RNA-Seq) data. We combined positional (locations within genome regulatory elements) and functional (associated with allele-specific binding and expression) criteria followed by an analysis using genome-wide association studies (GWAS) and minor allele frequency (MAF) datasets. DeSeq2 analysis through 70 CRC patients reinforced the regulatory potential. rSNPs (1,476) that were associated with significant (P < 0.01) allele-specific events resulting in thirty that exhibited a link with CRC according to the MAF and 27, with a risk of malignancy in general according to GWAS. Selected rSNPs may modify the expression of genes for tumor suppressors and the regulators of signaling pathways, including noncoding RNAs. However, the rSNPs from the most represented group affect the expression of genes related to splicing. Our findings strongly suggest that the identified variants might contribute to CRC susceptibility, which indicates that aberrant splicing is one of the key mechanisms for unraveling disease etiopathogenesis and provides useful inputs for interpreting how genotypic variation corresponds to phenotypic outcome.

PolyC-binding proteins enhance expression of the CDK2 cell cycle regulatory protein via alternative splicing

The PolyC binding proteins (PCBPs) impact alternative splicing of a subset of mammalian genes that are enriched in basic cellular functions. Here, we focus our analysis on PCBP-controlled cassette exon-splicing within the cell cycle control regulator cyclin-dependent kinase-2 (CDK2) transcript. We demonstrate that PCBP binding to a C-rich polypyrimidine tract (PPT) preceding exon 5 of the CDK2 transcript enhances cassette exon inclusion. This splice enhancement is U2AF65-independent and predominantly reflects actions of the PCBP1 isoform. Remarkably, PCBPs' control of CDK2 ex5 splicing has evolved subsequent to mammalian divergence via conversion of constitutive exon 5 inclusion in the mouse CDK2 transcript to PCBP-responsive exon 5 alternative splicing in humans. Importantly, exclusion of exon 5 from the hCDK2 transcript dramatically represses the expression of CDK2 protein with a corresponding perturbation in cell cycle kinetics. These data highlight a recently evolved post-transcriptional pathway in primate species with the potential to modulate cell cycle control.

The cancer-associated U2AF35 470A>G (Q157R) mutation creates an in-frame alternative 5' splice site that impacts splicing regulation in Q157R patients

Recent work has identified cancer-associated U2AF35 missense mutations in two zinc-finger (ZnF) domains, but little is known about Q157R/P substitutions within the second ZnF. Surprisingly, we find that the c.470A>G mutation not only leads to the Q157R substitution, but also creates an alternative 5' splice site (ss) resulting in the deletion of four amino acids (Q157Rdel). Q157P, Q157R, and Q157Rdel control alternative splicing of distinct groups of exons in cell culture and in human patients, suggesting that missplicing of different targets may contribute to cellular aberrations. Our data emphasize the importance to explore missense mutations beyond altered protein sequence.

Developing new targeting strategy for androgen receptor variants in castration resistant prostate cancer

The presence of androgen receptor variant 7 (AR-V7) variants becomes a significant hallmark of castration-resistant prostate cancer (CRPC) relapsed from hormonal therapy and is associated with poor survival of CRPC patients because of lacking a ligand-binding domain. Currently, it still lacks an effective agent to target AR-V7 or AR-Vs in general. Here, we showed that a novel class of agents (thailanstatins, TSTs and spliceostatin A analogs) can significantly suppress the expression of AR-V7 mRNA and protein but in a less extent on the full-length AR expression. Mechanistically, TST-D is able to inhibit AR-V7 gene splicing by interfering the interaction between U2AF65 and SAP155 and preventing them from binding to polypyrimidine tract located between the branch point and the 3' splice site. In vivo, TST-D exhibits a potent tumor inhibitory effect on human CRPC xenografts leading to cell apoptosis. The machinery associated with AR gene splicing in CRPC is a potential target for drugs. Based on their potency in the suppression of AR-V7 responsible for the growth/survival of CRPC, TSTs representing a new class of anti-AR-V agents warrant further development into clinical application.

CLIP-seq analysis of multi-mapped reads discovers novel functional RNA regulatory sites in the human transcriptome

Crosslinking or RNA immunoprecipitation followed by sequencing (CLIP-seq or RIP-seq) allows transcriptome-wide discovery of RNA regulatory sites. As CLIP-seq/RIP-seq reads are short, existing computational tools focus on uniquely mapped reads, while reads mapped to multiple loci are discarded. We present CLAM (CLIP-seq Analysis of Multi-mapped reads). CLAM uses an expectation–maximization algorithm to assign multi-mapped reads and calls peaks combining uniquely and multi-mapped reads. To demonstrate the utility of CLAM, we applied it to a wide range of public CLIP-seq/RIP-seq datasets involving numerous splicing factors, microRNAs and m⁶A RNA methylation. CLAM recovered a large number of novel RNA regulatory sites inaccessible by uniquely mapped reads. The functional significance of these sites was demonstrated by consensus motif patterns and association with alternative splicing (splicing factors), transcript abundance (AGO2) and mRNA half-life (m⁶A). CLAM provides a useful tool to discover novel protein–RNA interactions and RNA modification sites from CLIP-seq and RIP-seq data, and reveals the significant contribution of repetitive elements to the RNA regulatory landscape of the human transcriptome.

Stoichiometries of U2AF35, U2AF65 and U2 snRNP reveal new early spliceosome assembly pathways

The selection of 3΄ splice sites (3΄ss) is an essential early step in mammalian RNA splicing reactions, but the processes involved are unknown. We have used single molecule methods to test whether the major components implicated in selection, the proteins U2AF35 and U2AF65 and the U2 snRNP, are able to recognize alternative candidate sites or are restricted to one pre-specified site. In the presence of adenosine triphosphate (ATP), all three components bind in a 1:1 stoichiometry with a 3΄ss. Pre-mRNA molecules with two alternative 3΄ss can be bound concurrently by two molecules of U2AF or two U2 snRNPs, so none of the components are restricted. However, concurrent occupancy inhibits splicing. Stoichiometric binding requires conditions consistent with coalescence of the 5΄ and 3΄ sites in a complex (I, initial), but if this cannot form the components show unrestricted and stochastic association. In the absence of ATP, when complex E forms, U2 snRNP association is unrestricted. However, if protein dephosphorylation is prevented, an I-like complex forms with stoichiometric association of U2 snRNPs and the U2 snRNA is base-paired to the pre-mRNA. Complex I differs from complex A in that the formation of complex A is associated with the loss of U2AF65 and 35.

Identification of high-confidence RNA regulatory elements by combinatorial classification of RNA-protein binding sites

Crosslinking immunoprecipitation sequencing (CLIP-seq) technologies have enabled researchers to characterize transcriptome-wide binding sites of RNA-binding protein (RBP) with high resolution. We apply a soft-clustering method, RBPgroup, to various CLIP-seq datasets to group together RBPs that specifically bind the same RNA sites. Such combinatorial clustering of RBPs helps interpret CLIP-seq data and suggests functional RNA regulatory elements. Furthermore, we validate two RBP–RBP interactions in cell lines. Our approach links proteins and RNA motifs known to possess similar biochemical and cellular properties and can, when used in conjunction with additional experimental data, identify high-confidence RBP groups and their associated RNA regulatory elements.

NEAT1 scaffolds RNA-binding proteins and the Microprocessor to globally enhance pri-miRNA processing

MicroRNA biogenesis is known to be modulated by a variety of RNA binding proteins (RBPs), but in most cases, individual RBPs appear to influence the processing of a small subset of target miRNAs. We herein report that the RNA binding NONO/PSF heterodimer binds a large number of expressed pri-miRNAs in HeLa cells to globally enhance pri-miRNA processing by the Drosha/DGCR8 Microprocessor. Because NONO/PSF are key components of paraspeckles organized by the lncRNA NEAT1, we further demonstrate that NEAT1 also has a profound effect on global pri-miRNA processing. Mechanistic dissection reveals that NEAT1 broadly interacts with NONO/PSF as well as many other RBPs, and that multiple RNA segments in NEAT1, including a “pseudo pri-miRNA” near its 3′ end, help attract the Microprocessor. These findings suggest a bird nest model for a large non-coding RNA to orchestrate efficient processing of almost an entire class of small non-coding RNAs in the nucleus.

JMJD6 and U2AF65 co-regulate alternative splicing in both JMJD6 enzymatic activity dependent and independent manner

JMJD6, a jumonji C (Jmj C) domain-containing protein demethylase and hydroxylase, has been implicated in an array of biological processes. It has been shown that JMJD6 interacts with and hydroxylates multiple serine/arginine-rich (SR) proteins and SR related proteins, including U2AF65, all of which are known to function in alternative splicing regulation. However, whether JMJD6 is widely involved in alternative splicing and the molecular mechanism underlying JMJD6-regulated alternative splicing have remained incompletely understood. Here, by using RASL-Seq, we investigated the functional impact of RNA-dependent interaction between JMJD6 and U2AF65, revealing that JMJD6 and U2AF65 co-regulated a large number of alternative splicing events. We further demonstrated the JMJD6 function in alternative splicing in jmjd6 knockout mice. Mechanistically, we showed that the enzymatic activity of JMJD6 was required for a subset of JMJD6-regulated splicing, and JMJD6-mediated lysine hydroxylation of U2AF65 could account for, at least partially, their co-regulated alternative splicing events, suggesting both JMJD6 enzymatic activity-dependent and independent control of alternative splicing. These findings reveal an intimate link between JMJD6 and U2AF65 in alternative splicing regulation, which has important implications in development and disease processes.

DDX54 regulates transcriptome dynamics during DNA damage response

The cellular response to genotoxic stress is mediated by a well-characterized network of DNA surveillance pathways. The contribution of post-transcriptional gene regulatory networks to the DNA damage response (DDR) has not been extensively studied. Here, we systematically identified RNA-binding proteins differentially interacting with polyadenylated transcripts upon exposure of human breast carcinoma cells to ionizing radiation (IR). Interestingly, more than 260 proteins, including many nucleolar proteins, showed increased binding to poly(A)+ RNA in IR-exposed cells. The functional analysis of DDX54, a candidate genotoxic stress responsive RNA helicase, revealed that this protein is an immediate-to-early DDR regulator required for the splicing efficacy of its target IR-induced pre-mRNAs. Upon IR exposure, DDX54 acts by increased interaction with a well-defined class of pre-mRNAs that harbor introns with weak acceptor splice sites, as well as by protein-protein contacts within components of U2 snRNP and spliceosomal B complex, resulting in lower intron retention and higher processing rates of its target transcripts. Because DDX54 promotes survival after exposure to IR, its expression and/or mutation rate may impact DDR-related pathologies. Our work indicates the relevance of many uncharacterized RBPs potentially involved in the DDR.

Splicing Factor Mutations in Cancer

Many cancers demonstrate aberrant splicing patterns that contribute to their development and progression. Recently, recurrent somatic mutations of genes encoding core subunits of the spliceosome have been identified in several different cancer types. These mutations are most common in hematologic malignancies like the myelodysplastic syndromes (MDS), acute myeloid leukemia, and chronic lymphocytic leukemia, but also in occur in several solid tumors at lower frequency. The most frequent mutations occur in SF3B1, U2AF1, SRSF2, and ZRSR2 and are largely exclusive of each other. Mutations in SF3B1, U2AF1, and SRSF2 acquire heterozygous missense mutations in specific codons, resembling oncogenes. ZRSR2 mutations include clear loss-of-function variants, a pattern more common to tumor suppressor genes. These splicing factors are associated with distinct clinical phenotypes and patterns of mutation in different malignancies. Mutations have both diagnostic and prognostic relevance. Splicing factor mutations appear to affect only a minority of transcripts which show little overlap by mutation type. How differences in splicing caused by somatic mutations of spliceosome subunits lead to oncogenesis is not clear and may involve different targets in each disease type. However, cells with mutated splicing machinery may be particularly vulnerable to further disruption of the spliceosome suggesting a novel strategy for the targeted therapy of cancers.

Splicing Factor Mutations in Myelodysplasias: Insights from Spliceosome Structures

Somatic mutations of pre-mRNA splicing factors recur among patients with myelodysplastic syndrome (MDS) and related malignancies. Although these MDS-relevant mutations alter the splicing of a subset of transcripts, the mechanisms by which these single amino acid substitutions change gene expression remain controversial. New structures of spliceosome intermediates and associated protein complexes shed light on the molecular interactions mediated by 'hotspots' of the SF3B1 and U2AF1 pre-mRNA splicing factors. The frequently mutated SF3B1 residues contact the pre-mRNA splice site. Based on structural homology with other spliceosome subunits, and recent findings of altered RNA binding by mutant U2AF1 proteins, we suggest that affected U2AF1 residues also contact pre-mRNA. Altered pre-mRNA recognition emerges as a molecular theme among MDS-relevant mutations of pre-mRNA splicing factors.

RNA and Proteins: Mutual Respect

Proteins and RNA are often found in ribonucleoprotein particles (RNPs), where they function in cellular processes to synthesize proteins (the ribosome), chemically modify RNAs (small nucleolar RNPs), splice pre-mRNAs (the spliceosome), and, on a larger scale, sequester RNAs, degrade them, or process them (P bodies, Cajal bodies, and nucleoli). Each RNA–protein interaction is a story in itself, as both molecules can change conformation, compete for binding sites, and regulate cellular functions. Recent studies of Xist long non-coding RNP, the U4/5/6 tri-small nuclear RNP complex, and an activated state of a spliceosome reveal new features of RNA interactions with proteins, and, although their stories are incomplete, they are already fascinating.

Large-scale analysis of branchpoint usage across species and cell lines

The coding sequence of each human pre-mRNA is interrupted, on average, by 11 introns that must be spliced out for proper gene expression. Each intron contains three obligate signals: a 5′ splice site, a branch site, and a 3′ splice site. Splice site usage has been mapped exhaustively across different species, cell types, and cellular states. In contrast, only a small fraction of branch sites have been identified even once. The few reported annotations of branch site are imprecise as reverse transcriptase skips several nucleotides while traversing a 2–5 linkage. Here, we report large-scale mapping of the branchpoints from deep sequencing data in three different species and in the SF3B1 K700E oncogenic mutant background. We have developed a novel method whereby raw lariat reads are refined by U2snRNP/pre-mRNA base-pairing models to return the largest current data set of branchpoint sequences with quality metrics. This analysis discovers novel modes of U2snRNA:pre-mRNA base-pairing conserved in yeast and provides insight into the biogenesis of intron circles. Finally, matching branch site usage with isoform selection across the extensive panel of ENCODE RNA-seq data sets offers insight into the mechanisms by which branchpoint usage drives alternative splicing.

Insights into the design and interpretation of iCLIP experiments

Background

Ultraviolet (UV) crosslinking and immunoprecipitation (CLIP) identifies the sites on RNAs that are in direct contact with RNA-binding proteins (RBPs). Several variants of CLIP exist, which require different computational approaches for analysis. This variety of approaches can create challenges for a novice user and can hamper insights from multi-study comparisons. Here, we produce data with multiple variants of CLIP and evaluate the data with various computational methods to better understand their suitability.

Results

We perform experiments for PTBP1 and eIF4A3 using individual-nucleotide resolution CLIP (iCLIP), employing either UV-C or photoactivatable 4-thiouridine (4SU) combined with UV-A crosslinking and compare the results with published data. As previously noted, the positions of complementary DNA (cDNA)-starts depend on cDNA length in several iCLIP experiments and we now find that this is caused by constrained cDNA-ends, which can result from the sequence and structure constraints of RNA fragmentation. These constraints are overcome when fragmentation by RNase I is efficient and when a broad cDNA size range is obtained. Our study also shows that if RNase does not efficiently cut within the binding sites, the original CLIP method is less capable of identifying the longer binding sites of RBPs. In contrast, we show that a broad size range of cDNAs in iCLIP allows the cDNA-starts to efficiently delineate the complete RNA-binding sites.

Conclusions

We demonstrate the advantage of iCLIP and related methods that can amplify cDNAs that truncate at crosslink sites and we show that computational analyses based on cDNAs-starts are appropriate for such methods.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-1130-x) contains supplementary material, which is available to authorized users.