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.

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.

Improved genome assembly of American alligator genome reveals conserved architecture of estrogen signaling.

The American alligator, Alligator mississippiensis, like all crocodilians, has temperature-dependent sex determination, in which the sex of an embryo is determined by the incubation temperature of the egg during a critical period of development. The lack of genetic differences between male and female alligators leaves open the question of how the genes responsible for sex determination and differentiation are regulated. Insight into this question comes from the fact that exposing an embryo incubated at male-producing temperature to estrogen causes it to develop ovaries. Because estrogen response elements are known to regulate genes over long distances, a contiguous genome assembly is crucial for predicting and understanding their impact. We present an improved assembly of the American alligator genome, scaffolded with in vitro proximity ligation (Chicago) data. We use this assembly to scaffold two other crocodilian genomes based on synteny. We perform RNA sequencing of tissues from American alligator embryos to find genes that are differentially expressed between embryos incubated at male- versus female-producing temperature. Finally, we use the improved contiguity of our assembly along with the current model of CTCF-mediated chromatin looping to predict regions of the genome likely to contain estrogen-responsive genes. We find that these regions are significantly enriched for genes with female-biased expression in developing gonads after the critical period during which sex is determined by incubation temperature. We thus conclude that estrogen signaling is a major driver of female-biased gene expression in the post-temperature sensitive period gonads.

Modulation of nonsense mediated decay by rapamycin.

Rapamycin is a naturally occurring macrolide whose target is at the core of nutrient and stress regulation in a wide range of species. Despite well-established roles as an inhibitor of cap-dependent mRNA translation, relatively little is known about its effects on other modes of RNA processing. Here, we characterize the landscape of rapamycin-induced post-transcriptional gene regulation. Transcriptome analysis of rapamycin-treated cells reveals genome-wide changes in alternative mRNA splicing and pronounced changes in NMD-sensitive isoforms. We demonstrate that despite well-documented attenuation of cap-dependent mRNA translation, rapamycin can augment NMD of certain transcripts. Rapamycin-treatment significantly reduces the levels of both endogenous and exogenous Premature Termination Codon (PTC)-containing mRNA isoforms and its effects are dose-, UPF1- and 4EBP-dependent. The PTC-containing SRSF6 transcript exhibits a shorter half-life upon rapamycin-treatment as compared to the non-PTC isoform. Rapamycin-treatment also causes depletion of PTC-containing mRNA isoforms from polyribosomes, underscoring the functional relationship between translation and NMD. Enhanced NMD activity also correlates with an enrichment of the nuclear Cap Binding Complex (CBC) in rapamycin-treated cells. Our data demonstrate that rapamycin modulates global RNA homeostasis by NMD.

Frac-seq reveals isoform-specific recruitment to polyribosomes.

Pre-mRNA splicing is required for the accurate expression of virtually all human protein coding genes. However, splicing also plays important roles in coordinating subsequent steps of pre-mRNA processing such as polyadenylation and mRNA export. Here, we test the hypothesis that nuclear pre-mRNA processing influences the polyribosome association of alternative mRNA isoforms. By comparing isoform ratios in cytoplasmic and polyribosomal extracts, we determined that the alternative products of ∼30% (597/1954) of mRNA processing events are differentially partitioned between these subcellular fractions. Many of the events exhibiting isoform-specific polyribosome association are highly conserved across mammalian genomes, underscoring their possible biological importance. We find that differences in polyribosome association may be explained, at least in part by the observation that alternative splicing alters the cis-regulatory landscape of mRNAs isoforms. For example, inclusion or exclusion of upstream open reading frames (uORFs) in the 5'UTR as well as Alu-elements and microRNA target sites in the 3'UTR have a strong influence on polyribosome association of alternative mRNA isoforms. Taken together, our data demonstrate for the first time the potential link between alternative splicing and translational control of the resultant mRNA isoforms.

The splicing factor SF2/ASF regulates translation initiation by enhancing phosphorylation of 4E-BP1.

The SR protein SF2/ASF has been initially characterized as a splicing factor but has also been shown to mediate postsplicing activities such as mRNA export and translation. Here we demonstrate that SF2/ASF promotes translation initiation of bound mRNAs and that this activity requires the presence of the cytoplasmic cap-binding protein eIF4E. SF2/ASF promotes translation initiation by suppressing the activity of 4E-BP, a competitive inhibitor of cap-dependent translation. This activity is mediated by interactions of SF2/ASF with both mTOR and the phosphatase PP2A, two key regulators of 4E-BP phosphorylation. These findings suggest the model whereby SF2/ASF functions as an adaptor protein to recruit the signaling molecules responsible for regulation of cap-dependent translation of specific mRNAs. Taken together, these data suggest a novel mechanism for the activation of translation initiation of a subset of mRNAs bound by the shuttling protein SF2/ASF.

Loss of exon identity is a common mechanism of human inherited disease.

It is widely accepted that at least 10% of all mutations causing human inherited disease disrupt splice-site consensus sequences. In contrast to splice-site mutations, the role of auxiliary cis-acting elements such as exonic splicing enhancers (ESE) and exonic splicing silencers (ESS) in human inherited disease is still poorly understood. Here we use a top-down approach to determine rates of loss or gain of known human exonic splicing regulatory (ESR) sequences associated with either disease-causing mutations or putatively neutral single nucleotide polymorphisms (SNPs). We observe significant enrichment toward loss of ESEs and gain of ESSs among inherited disease-causing variants relative to neutral polymorphisms, indicating that exon skipping may play a prominent role in aberrant gene regulation. Both computational and biochemical approaches underscore the relevance of exonic splicing enhancer loss and silencer gain in inherited disease. Additionally, we provide direct evidence that both SRp20 (SRSF3) and possibly PTB (PTBP1) are involved in the function of a splicing silencer that is created de novo by a total of 83 different inherited disease mutations in 67 different disease genes. Taken together, we find that ~25% (7154/27,681) of known mis-sense and nonsense disease-causing mutations alter functional splicing signals within exons, suggesting a much more widespread role for aberrant mRNA processing in causing human inherited disease than has hitherto been appreciated.

junctionCounts: comprehensive alternative splicing analysis and prediction of isoform-level impacts to the coding sequence.

Alternative splicing (AS) is emerging as an important regulatory process for complex biological processes. Transcriptomic studies therefore commonly involve the identification and quantification of alternative processing events, but the need for predicting the functional consequences of changes to the relative inclusion of alternative events remains largely unaddressed. Many tools exist for the former task, albeit each constrained to its own event type definitions. Few tools exist for the latter task; each with significant limitations. To address these issues we developed junctionCounts, which captures both simple and complex pairwise AS events and quantifies them with straightforward exon-exon and exon-intron junction reads in RNA-seq data, performing competitively among similar tools in terms of sensitivity, false discovery rate and quantification accuracy. Its partner utility, cdsInsertion, identifies transcript coding sequence (CDS) information via in silico translation from annotated start codons, including the presence of premature termination codons. Finally, findSwitchEvents connects AS events with CDS information to predict the impact of individual events to the isoform-level CDS. We used junctionCounts to characterize splicing dynamics and NMD regulation during neuronal differentiation across four primates, demonstrating junctionCounts' capacity to robustly characterize AS in a variety of organisms and to predict its effect on mRNA isoform fate.

Targeting IGF2BP3 enhances antileukemic effects of menin-MLL inhibition in MLL-AF4 leukemia.

ABSTRACT: RNA-binding proteins (RBPs) are emerging as a novel class of therapeutic targets in cancer, including in leukemia, given their important role in posttranscriptional gene regulation, and have the unexplored potential to be combined with existing therapies. The RBP insulin-like growth factor 2 messenger RNA-binding protein 3 (IGF2BP3) has been found to be a critical regulator of MLL-AF4 leukemogenesis and represents a promising therapeutic target. Here, we study the combined effects of targeting IGF2BP3 and menin-MLL interaction in MLL-AF4-driven leukemia in vitro and in vivo, using genetic inhibition with CRISPR-Cas9-mediated deletion of Igf2bp3 and pharmacologic inhibition of the menin-MLL interaction with multiple commercially available inhibitors. Depletion of Igf2bp3 sensitized MLL-AF4 leukemia to the effects of menin-MLL inhibition on cell growth and leukemic initiating cells in vitro. Mechanistically, we found that both Igf2bp3 depletion and menin-MLL inhibition led to increased differentiation in vitro and in vivo, seen in functional readouts and by gene expression analyses. IGF2BP3 knockdown had a greater effect on increasing survival and attenuating disease than pharmacologic menin-MLL inhibition with small molecule MI-503 alone and showed enhanced antileukemic effects in combination. Our work shows that IGF2BP3 is an oncogenic amplifier of MLL-AF4-mediated leukemogenesis and a potent therapeutic target, providing a paradigm for targeting leukemia at both the transcriptional and posttranscriptional level.

Site identification in high-throughput RNA-protein interaction data.

MOTIVATION: Post-transcriptional and co-transcriptional regulation is a crucial link between genotype and phenotype. The central players are the RNA-binding proteins, and experimental technologies [such as cross-linking with immunoprecipitation- (CLIP-) and RIP-seq] for probing their activities have advanced rapidly over the course of the past decade. Statistically robust, flexible computational methods for binding site identification from high-throughput immunoprecipitation assays are largely lacking however. RESULTS: We introduce a method for site identification which provides four key advantages over previous methods: (i) it can be applied on all variations of CLIP and RIP-seq technologies, (ii) it accurately models the underlying read-count distributions, (iii) it allows external covariates, such as transcript abundance (which we demonstrate is highly correlated with read count) to inform the site identification process and (iv) it allows for direct comparison of site usage across cell types or conditions. AVAILABILITY AND IMPLEMENTATION: We have implemented our method in a software tool called Piranha. Source code and binaries, licensed under the GNU General Public License (version 3) are freely available for download from http://smithlab.usc.edu. CONTACT: andrewds@usc.edu SUPPLEMENTARY INFORMATION: Supplementary data available at Bioinformatics online.

Distinct oncogenic phenotypes in hematopoietic specific deletions of Trp53.

Loss of function in the tumor suppressor gene TP53 is the most common alteration seen in human cancer. In mice, P53 deletion in all cells leads predominantly to the development of T-cell lymphomas, followed by B-cell lymphomas, sarcomas and teratomas. In order to dissect the role of P53 in the hematopoietic system, we generated and analyzed two different mouse models deficient for P53. A pan-hematopoietic P53 deletion mouse was created using Vav1-Cre based deletion; and a B-cell-specific deletion mouse was created using a CD19-Cre based deletion. The Vav1-P53CKO mice predominantly developed T-cell malignancies in younger mice, and myeloid malignancies in older mice. In T-cell malignancies, there was accelerated thymic cell maturation with overexpression of Notch1 and its downstream effectors. CD19-P53CKO mice developed marginal zone expansion in the spleen, followed by marginal zone lymphoma, some of which progressed to diffuse large B-cell lymphomas. Interestingly, marginal zone and diffuse large B-cell lymphomas had a unique gene expression signature characterized by activation of the PI3K pathway, compared with wild type marginal zone or follicular cells of the spleen. This study demonstrates lineage specific P53 deletion leading to distinct phenotypes secondary to unique gene expression programs set in motion.

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.

A potential role for SARS-CoV-2 small viral RNAs in targeting host microRNAs and modulating gene expression.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes coronavirus disease (COVID-19) in humans, with symptoms ranging from mild to severe, including fatality. The molecular mechanisms surrounding the effects of viral infection on the host RNA machinery remain poorly characterized. We used a comparative transcriptomics approach to investigate the effects of SARS-CoV-2 infection on the host mRNA and sRNA expression machinery in a human lung epithelial cell line (Calu-3) and an African green monkey kidney cell line (Vero-E6). Upon infection, we observed global changes in host gene expression and differential expression of dozens of host miRNAs, many with known links to viral infection and immune response. Additionally, we discovered an expanded landscape of more than a hundred SARS-CoV-2-derived small viral RNAs (svRNAs) predicted to interact with differentially expressed host mRNAs and miRNAs. svRNAs are derived from distinct regions of the viral genome and sequence signatures suggest they are produced by a non-canonical biogenesis pathway. 52 of the 67 svRNAs identified in Calu-3 cells are predicted to interact with differentially expressed miRNAs, with many svRNAs having multiple targets. Accordingly, we speculate that these svRNAs may play a role in SARS-CoV-2 propagation by modulating post-transcriptional gene regulation, and that methods for antagonizing them may have therapeutic value.

The RNA-binding protein IGF2BP3 is critical for MLL-AF4-mediated leukemogenesis.

Despite recent advances in therapeutic approaches, patients with MLL-rearranged leukemia still have poor outcomes. Here, we find that the RNA-binding protein IGF2BP3, which is overexpressed in MLL-translocated leukemia, strongly amplifies MLL-Af4-mediated leukemogenesis. Deletion of Igf2bp3 significantly increases the survival of mice with MLL-Af4-driven leukemia and greatly attenuates disease, with a minimal impact on baseline hematopoiesis. At the cellular level, MLL-Af4 leukemia-initiating cells require Igf2bp3 for their function in leukemogenesis. At the molecular level, IGF2BP3 regulates a complex posttranscriptional operon governing leukemia cell survival and proliferation. IGF2BP3-targeted mRNA transcripts include important MLL-Af4-induced genes, such as those in the Hoxa locus, and the Ras signaling pathway. Targeting of transcripts by IGF2BP3 regulates both steady-state mRNA levels and, unexpectedly, pre-mRNA splicing. Together, our findings show that IGF2BP3 represents an attractive therapeutic target in this disease, providing important insights into mechanisms of posttranscriptional regulation in leukemia.

Predicting sequence and structural specificities of RNA binding regions recognized by splicing factor SRSF1.

BACKGROUND: RNA-binding proteins (RBPs) play diverse roles in eukaryotic RNA processing. Despite their pervasive functions in coding and noncoding RNA biogenesis and regulation, elucidating the sequence specificities that define protein-RNA interactions remains a major challenge. Recently, CLIP-seq (Cross-linking immunoprecipitation followed by high-throughput sequencing) has been successfully implemented to study the transcriptome-wide binding patterns of SRSF1, PTBP1, NOVA and fox2 proteins. These studies either adopted traditional methods like Multiple EM for Motif Elicitation (MEME) to discover the sequence consensus of RBP's binding sites or used Z-score statistics to search for the overrepresented nucleotides of a certain size. We argue that most of these methods are not well-suited for RNA motif identification, as they are unable to incorporate the RNA structural context of protein-RNA interactions, which may affect to binding specificity. Here, we describe a novel model-based approach--RNAMotifModeler to identify the consensus of protein-RNA binding regions by integrating sequence features and RNA secondary structures. RESULTS: As an example, we implemented RNAMotifModeler on SRSF1 (SF2/ASF) CLIP-seq data. The sequence-structural consensus we identified is a purine-rich octamer 'AGAAGAAG' in a highly single-stranded RNA context. The unpaired probabilities, the probabilities of not forming pairs, are significantly higher than negative controls and the flanking sequence surrounding the binding site, indicating that SRSF1 proteins tend to bind on single-stranded RNA. Further statistical evaluations revealed that the second and fifth bases of SRSF1octamer motif have much stronger sequence specificities, but weaker single-strandedness, while the third, fourth, sixth and seventh bases are far more likely to be single-stranded, but have more degenerate sequence specificities. Therefore, we hypothesize that nucleotide specificity and secondary structure play complementary roles during binding site recognition by SRSF1. CONCLUSION: In this study, we presented a computational model to predict the sequence consensus and optimal RNA secondary structure for protein-RNA binding regions. The successful implementation on SRSF1 CLIP-seq data demonstrates great potential to improve our understanding on the binding specificity of RNA binding proteins.

A method for campus-wide SARS-CoV-2 surveillance at a large public university.

The systematic screening of asymptomatic and pre-symptomatic individuals is a powerful tool for controlling community transmission of infectious disease on college campuses. Faced with a paucity of testing in the beginning of the COVID-19 pandemic, many universities developed molecular diagnostic laboratories focused on SARS-CoV-2 diagnostic testing on campus and in their broader communities. We established the UC Santa Cruz Molecular Diagnostic Lab in early April 2020 and began testing clinical samples just five weeks later. Using a clinically-validated laboratory developed test (LDT) that avoided supply chain constraints, an automated sample pooling and processing workflow, and a custom laboratory information management system (LIMS), we expanded testing from a handful of clinical samples per day to thousands per day with the testing capacity to screen our entire campus population twice per week. In this report we describe the technical, logistical, and regulatory processes that enabled our pop-up lab to scale testing and reporting capacity to thousands of tests per day.

Concerted regulation of nuclear and cytoplasmic activities of SR proteins by AKT.

Serine/arginine-rich (SR) proteins are important regulators of mRNA splicing. Several postsplicing activities have been described for a subset of shuttling SR proteins, including regulation of mRNA export and translation. Using the fibronectin gene to study the links between signal-transduction pathways and SR protein activity, we show that growth factors not only modify the alternative splicing pattern of the fibronectin gene but also alter translation of reporter messenger RNAs in an SR protein-dependent fashion, providing two coregulated levels of isoform-specific amplification. These effects are inhibited by specific small interfering RNAs against SR proteins and are mediated by the AKT kinase, which elicits opposite effects to those evoked by overexpressing SR protein kinases Clk and SRPK. These results show how SR protein activity is modified in response to extracellular stimulation, leading to a concerted regulation of splicing and translation.

Genome-wide prediction of cis-acting RNA elements regulating tissue-specific pre-mRNA alternative splicing.

BACKGROUND: Human genes undergo various patterns of pre-mRNA splicing across different tissues. Such variation is primarily regulated by trans-acting factors that bind on exonic and intronic cis-acting RNA elements (CAEs). Here we report a computational method to mechanistically identify cis-acting RNA elements that contribute to the tissue-specific alternative splicing pattern. This method is an extension of our previous model, SplicingModeler, which predicts the significant CAEs that contribute to the splicing differences between two tissues. In this study, we introduce tissue-specific functional levels estimation step, which allows evaluating regulatory functions of predicted CAEs that are involved in more than two tissues. RESULTS: Using a publicly available Affymetrix Genechip Human Exon Array dataset, our method identifies 652 cis-acting RNA elements (CAEs) across 11 human tissues. About one third of predicted CAEs can be mapped to the known RBP (RNA binding protein) binding sites or match with other predicted exonic splicing regulator databases. Interestingly, the vast majority of predicted CAEs are in intronic regulatory regions. A noticeable exception is that many exonic elements are found to regulate the alternative splicing between cerebellum and testes. Most identified elements are found to contribute to the alternative splicing between two tissues, while some are important in multiple tissues. This suggests that genome-wide alternative splicing patterns are regulated by a combination of tissue-specific cis-acting elements and "general elements" whose functional activities are important but differ across multiple tissues. CONCLUSION: In this study, we present a model-based computational approach to identify potential cis-acting RNA elements by considering the exon splicing variation as the combinatorial effects of multiple cis-acting regulators. This methodology provides a novel evaluation on the functional levels of cis-acting RNA elements by estimating their tissue-specific functions on various tissues.

Using RNase sequence specificity to refine the identification of RNA-protein binding regions.

Massively parallel pyrosequencing is a high-throughput technology that can sequence hundreds of thousands of DNA/RNA fragments in a single experiment. Combining it with immunoprecipitation-based biochemical assays, such as cross-linking immunoprecipitation (CLIP), provides a genome-wide method to detect the sites at which proteins bind DNA or RNA. In a CLIP-pyrosequencing experiment, the resolutions of the detected protein binding regions are partially determined by the length of the detected RNA fragments (CLIP amplicons) after trimming by RNase digestion. The lengths of these fragments usually range from 50-70 nucleotides. Many genomic regions are marked by multiple RNA fragments. In this paper, we report an empirical approach to refine the localization of protein binding regions by using the distribution pattern of the detected RNA fragments and the sequence specificity of RNase digestion. We present two regions to which multiple amplicons map as examples to demonstrate this approach.