Position-dependent splicing activation and repression by SR and hnRNP proteins rely on common mechanisms

The HIV 5' Gag Region Displays a Specific Nucleotide Bias Regulating Viral Splicing and Infectivity

Alternative splicing and the expression of intron-containing mRNAs is one hallmark of HIV gene expression. To facilitate the otherwise hampered nuclear export of non-fully processed mRNAs, HIV encodes the Rev protein, which recognizes its intronic response element and fuels the HIV RNAs into the CRM-1-dependent nuclear protein export pathway. Both alternative splicing and Rev-dependency are regulated by the primary HIV RNA sequence. Here, we show that these processes are extremely sensitive to sequence alterations in the 5’coding region of the HIV genomic RNA. Increasing the GC content by insertion of either GFP or silent mutations activates a cryptic splice donor site in gag, entirely deregulates the viral splicing pattern, and lowers infectivity. Interestingly, an adaptation of the inserted GFP sequence toward an HIV-like nucleotide bias reversed these phenotypes completely. Of note, the adaptation yielded completely different primary sequences although encoding the same amino acids. Thus, the phenotypes solely depend on the nucleotide composition of the two GFP versions. This is a strong indication of an HIV-specific mRNP code in the 5′ gag region wherein the primary RNA sequence bias creates motifs for RNA-binding proteins and controls the fate of the HIV-RNA in terms of viral gene expression and infectivity.

Modifying splice site usage with ModCon: Maintaining the genetic code while changing the underlying mRNP code

Codon degeneracy of amino acid sequences permits an additional “mRNP code” layer underlying the genetic code that is related to RNA processing. In pre-mRNA splicing, splice site usage is determined by both intrinsic strength and sequence context providing RNA binding sites for splicing regulatory proteins. In this study, we systematically examined modification of splicing regulatory properties in the neighborhood of a GT site, i.e. potential splice site, without altering the encoded amino acids.

We quantified the splicing regulatory properties of the neighborhood around a potential splice site by its Splice Site HEXplorer Weight (SSHW) based on the HEXplorer score algorithm. To systematically modify GT site neighborhoods, either minimizing or maximizing their SSHW, we designed the novel stochastic optimization algorithm ModCon that applies a genetic algorithm with stochastic crossover, insertion and random mutation elements supplemented by a heuristic sliding window approach.

To assess the achievable range in SSHW in human splice donors without altering the encoded amino acids, we applied ModCon to a set of 1000 randomly selected Ensembl annotated human splice donor sites, achieving substantial and accurate changes in SSHW. Using ModCon optimization, we successfully switched splice donor usage in a splice site competition reporter containing coding sequences from FANCA, FANCB or BRCA2, while retaining their amino acid coding information.

The ModCon algorithm and its R package implementation can assist in reporter design by either introducing novel splice sites, silencing accidental, undesired splice sites, and by generally modifying the entire mRNP code while maintaining the genetic code.

Splicing factor SRSF1 promotes breast cancer progression via oncogenic splice switching of PTPMT1

Background

Intensive evidence has highlighted the effect of aberrant alternative splicing (AS) events on cancer progression when triggered by dysregulation of the SR protein family. Nonetheless, the underlying mechanism in breast cancer (BRCA) remains elusive. Here we sought to explore the molecular function of SRSF1 and identify the key AS events regulated by SRSF1 in BRCA.

Methods

We conducted a comprehensive analysis of the expression and clinical correlation of SRSF1 in BRCA based on the TCGA dataset, Metabric database and clinical tissue samples. Functional analysis of SRSF1 in BRCA was conducted in vitro and in vivo. SRSF1-mediated AS events and their binding motifs were identified by RNA-seq, RNA immunoprecipitation-PCR (RIP-PCR) and in vivo crosslinking followed by immunoprecipitation (CLIP), which was further validated by the minigene reporter assay. PTPMT1 exon 3 (E3) AS was identified to partially mediate the oncogenic role of SRSF1 by the P-AKT/C-MYC axis. Finally, the expression and clinical significance of these AS events were validated in clinical samples and using the TCGA database.

Results

SRSF1 expression was consistently upregulated in BRCA samples, positively associated with tumor grade and the Ki-67 index, and correlated with poor prognosis in a hormone receptor-positive (HR+) cohort, which facilitated proliferation, cell migration and inhibited apoptosis in vitro and in vivo. We identified SRSF1-mediated AS events and discovered the SRSF1 binding motif in the regulation of splice switching of PTPMT1. Furthermore, PTPMT1 splice switching was regulated by SRSF1 by binding directly to its motif in E3 which partially mediated the oncogenic role of SRSF1 by the AKT/C-MYC axis. Additionally, PTPMT1 splice switching was validated in tissue samples of BRCA patients and using the TCGA database. The high-risk group, identified by AS of PTPMT1 and expression of SRSF1, possessed poorer prognosis in the stage I/II TCGA BRCA cohort.

Conclusions

SRSF1 exerts oncogenic roles in BRCA partially by regulating the AS of PTPMT1, which could be a therapeutic target candidate in BRCA and a prognostic factor in HR+ BRCA patient.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-021-01978-8.

Transcription Regulation Through Nascent RNA Folding

Folding of RNA into secondary structures through intramolecular base pairing determines an RNA's three-dimensional architecture and associated function. Simple RNA structures like stem loops can provide specialized functions independent of coding capacity, such as protein binding, regulation of RNA processing and stability, stimulation or inhibition of translation. RNA catalysis is dependent on tertiary structures found in the ribosome, tRNAs and group I and II introns. While the extent to which non-coding RNAs contribute to cellular maintenance is generally appreciated, the fact that both non-coding and coding RNA can assume relevant structural states has only recently gained attention. In particular, the co-transcriptional folding of nascent RNA of all classes has the potential to regulate co-transcriptional processing, RNP (ribonucleoprotein particle) formation, and transcription itself. Riboswitches are established examples of co-transcriptionally folded coding RNAs that directly regulate transcription, mainly in prokaryotes. Here we discuss recent studies in both prokaryotes and eukaryotes showing that structure formation may carry a more widespread regulatory logic during RNA synthesis. Local structures forming close to the catalytic center of RNA polymerases have the potential to regulate transcription by reducing backtracking. In addition, stem loops or more complex structures may alter co-transcriptional RNA processing or its efficiency. Several examples of functional structures have been identified to date, and this review provides an overview of physiologically distinct processes where co-transcriptionally folded RNA plays a role. Experimental approaches such as single-molecule FRET and in vivo structural probing to further advance our insight into the significance of co-transcriptional structure formation are discussed.

Splicing at the phase-separated nuclear speckle interface: a model

Phase-separated membraneless bodies play important roles in nucleic acid biology. While current models for the roles of phase separation largely focus on the compartmentalization of constituent proteins, we reason that other properties of phase separation may play functional roles. Specifically, we propose that interfaces of phase-separated membraneless bodies could have functional roles in spatially organizing biochemical reactions. Here we propose such a model for the nuclear speckle, a membraneless body implicated in RNA splicing. In our model, sequence-dependent RNA positioning along the nuclear speckle interface coordinates RNA splicing. Our model asserts that exons are preferentially sequestered into nuclear speckles through binding by SR proteins, while introns are excluded through binding by nucleoplasmic hnRNP proteins. As a result, splice sites at exon-intron boundaries are preferentially positioned at nuclear speckle interfaces. This positioning exposes splice sites to interface-localized spliceosomes, enabling the subsequent splicing reaction. Our model provides a simple mechanism that seamlessly explains much of the complex logic of splicing. This logic includes experimental results such as the antagonistic duality between splicing factors, the position dependence of splicing sequence motifs, and the collective contribution of many motifs to splicing decisions. Similar functional roles for phase-separated interfaces may exist for other membraneless bodies.

Comprehensive characterisation of intronic mis-splicing mutations in human cancers

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

Intronic tRNAs of mitochondrial origin regulate constitutive and alternative splicing

BACKGROUND

The presence of nuclear mitochondrial DNA (numtDNA) has been reported within several nuclear genomes. Next to mitochondrial protein-coding genes, numtDNA sequences also encode for mitochondrial tRNA genes. However, the biological roles of numtDNA remain elusive.

RESULTS

Employing in silico analysis, we identify 281 mitochondrial tRNA homologs in the human genome, which we term nimtRNAs (nuclear intronic mitochondrial-derived tRNAs), being contained within introns of 76 nuclear host genes. Despite base changes in nimtRNAs when compared to their mtRNA homologs, a canonical tRNA cloverleaf structure is maintained. To address potential functions of intronic nimtRNAs, we insert them into introns of constitutive and alternative splicing reporters and demonstrate that nimtRNAs promote pre-mRNA splicing, dependent on the number and positioning of nimtRNA genes and splice site recognition efficiency. A mutational analysis reveals that the nimtRNA cloverleaf structure is required for the observed splicing increase. Utilizing a CRISPR/Cas9 approach, we show that a partial deletion of a single endogenous nimtRNALys within intron 28 of the PPFIBP1 gene decreases inclusion of the downstream-located exon 29 of the PPFIBP1 mRNA. By employing a pull-down approach followed by mass spectrometry, a 3'-splice site-associated protein network is identified, including KHDRBS1, which we show directly interacts with nimtRNATyr by an electrophoretic mobility shift assay.

CONCLUSIONS

We propose that nimtRNAs, along with associated protein factors, can act as a novel class of intronic splicing regulatory elements in the human genome by participating in the regulation of splicing.

The role of hnRNPs in frontotemporal dementia and amyotrophic lateral sclerosis

Dysregulated RNA metabolism is emerging as a crucially important mechanism underpinning the pathogenesis of frontotemporal dementia (FTD) and the clinically, genetically and pathologically overlapping disorder of amyotrophic lateral sclerosis (ALS). Heterogeneous nuclear ribonucleoproteins (hnRNPs) comprise a family of RNA-binding proteins with diverse, multi-functional roles across all aspects of mRNA processing. The role of these proteins in neurodegeneration is far from understood. Here, we review some of the unifying mechanisms by which hnRNPs have been directly or indirectly linked with FTD/ALS pathogenesis, including their incorporation into pathological inclusions and their best-known roles in pre-mRNA splicing regulation. We also discuss the broader functionalities of hnRNPs including their roles in cryptic exon repression, stress granule assembly and in co-ordinating the DNA damage response, which are all emerging pathogenic themes in both diseases. We then present an integrated model that depicts how a broad-ranging network of pathogenic events can arise from declining levels of functional hnRNPs that are inadequately compensated for by autoregulatory means. Finally, we provide a comprehensive overview of the most functionally relevant cellular roles, in the context of FTD/ALS pathogenesis, for hnRNPs A1-U.

Viral Perturbation of Alternative Splicing of a Host Transcript Benefits Infection

Pathogens disturb alternative splicing patterns of infected eukaryotic hosts. However, in plants it is unknown if this is incidental to infection or represents a pathogen-induced remodeling of host gene expression needed to support infection. Here, we compared changes in transcription and protein accumulation with changes in transcript splicing patterns in maize (Zea mays) infected with the globally important pathogen sugarcane mosaic virus (SCMV). Our results suggested that changes in alternative splicing play a major role in determining virus-induced proteomic changes. Focusing on maize phytoene synthase1 (ZmPSY1), which encodes the key regulatory enzyme in carotenoid biosynthesis, we found that although SCMV infection decreases total ZmPSY1 transcript accumulation, the proportion of splice variant T001 increases by later infection stages so that ZmPSY1 protein levels are maintained. We determined that ZmPSY1 has two leaf-specific transcripts, T001 and T003, distinguished by differences between the respective 3'-untranslated regions (UTRs). The shorter 3'-UTR of T001 makes it the more efficient mRNA. Nonsense ZmPSY1 mutants or virus-induced silencing of ZmPSY1 expression suppressed SCMV accumulation, attenuated symptoms, and decreased chloroplast damage. Thus, ZmPSY1 acts as a proviral host factor that is required for virus accumulation and pathogenesis. Taken together, our findings reveal that SCMV infection-modulated alternative splicing ensures that ZmPSY1 synthesis is sustained during infection, which supports efficient virus infection.

Allosteric regulation of U1 snRNP by splicing regulatory proteins controls spliceosomal assembly

Alternative splicing is responsible for much of the transcriptomic and proteomic diversity observed in eukaryotes and involves combinatorial regulation by many cis-acting elements and trans-acting factors. SR and hnRNP splicing regulatory proteins often have opposing effects on splicing efficiency depending on where they bind the pre-mRNA relative to the splice site. Position-dependent splicing repression occurs at spliceosomal E-complex, suggesting that U1 snRNP binds but cannot facilitate higher order spliceosomal assembly. To test the hypothesis that the structure of U1 snRNA changes during activation or repression, we developed a method to structure-probe native U1 snRNP in enriched conformations that mimic activated or repressed spliceosomal E-complexes. While the core of U1 snRNA is highly structured, the 5' end of U1 snRNA shows different SHAPE reactivities and psoralen crosslinking efficiencies depending on where splicing regulatory elements are located relative to the 5' splice site. A motif within the 5' splice site binding region of U1 snRNA is more reactive toward SHAPE electrophiles when repressors are bound, suggesting U1 snRNA is bound, but less base-paired. These observations demonstrate that splicing regulators modulate splice site selection allosterically.

The role of alternative splicing in cancer: From oncogenesis to drug resistance

Alternative splicing is a tightly regulated process whereby non-coding sequences of pre-mRNA are removed and protein-coding segments are assembled in diverse combinations, ultimately giving rise to proteins with distinct or even opposing functions. In the past decade, whole genome/transcriptome sequencing studies revealed the high complexity of splicing regulation, which occurs co-transcriptionally and is influenced by chromatin status and mRNA modifications. Consequently, splicing profiles of both healthy and malignant cells display high diversity and alternative splicing was shown to be widely deregulated in multiple cancer types. In particular, mutations in pre-mRNA regulatory sequences, splicing regulators and chromatin modifiers, as well as differential expression of splicing factors are important contributors to cancer pathogenesis. It has become clear that these aberrations contribute to many facets of cancer, including oncogenic transformation, cancer progression, response to anticancer drug treatment as well as resistance to therapy. In this respect, alternative splicing was shown to perturb the expression a broad spectrum of relevant genes involved in drug uptake/metabolism (i.e. SLC29A1, dCK, FPGS, and TP), activation of nuclear receptor pathways (i.e. GR, AR), regulation of apoptosis (i.e. MCL1, BCL-X, and FAS) and modulation of response to immunotherapy (CD19). Furthermore, aberrant splicing constitutes an important source of novel cancer biomarkers and the spliceosome machinery represents an attractive target for a novel and rapidly expanding class of therapeutic agents. Small molecule inhibitors targeting SF3B1 or splice factor kinases were highly cytotoxic against a wide range of cancer models, including drug-resistant cells. Importantly, these effects are enhanced in specific cancer subsets, such as splicing factor-mutated and c-MYC-driven tumors. Furthermore, pre-clinical studies report synergistic effects of spliceosome modulators in combination with conventional antitumor agents. These strategies based on the use of low dose splicing modulators could shift the therapeutic window towards decreased toxicity in healthy tissues. Here we provide an extensive overview of the latest findings in the field of regulation of splicing in cancer, including molecular mechanisms by which cancer cells harness alternative splicing to drive oncogenesis and evade anticancer drug treatment as well as splicing-based vulnerabilities that can provide novel treatment opportunities. Furthermore, we discuss current challenges arising from genome-wide detection and prediction methods of aberrant splicing, as well as unravelling functional relevance of the plethora of cancer-related splicing alterations.

Long non-coding RNA AGAP2-AS1 accelerates cell proliferation, migration, invasion and the EMT process in colorectal cancer via regulating the miR-4,668-3p/SRSF1 axis

BACKGROUND

Colorectal cancer (CRC) is a frequently occurring tumor. Although a number of long noncoding RNAs have been researched in CRC, the expression, function and mechanism of AGAP2-AS1 remains poorly investigated.

METHODS

Gene expression was analyzed by a quantitative reverse transcriptase-polymerase chain rreaction and western blot analyses. Cell counting kit-8, colony formation and Transwell assays were utilized to explore the functional role of AGAP2-AS1 in CRC. Luciferase reporter, RNA pull down and RNA immunoprecipitation assays were implemented to verify relationships between RNA molecules.

RESULTS

In the present study, AGAP2-AS1 was unveiled as highly expressed in CRC cell lines compared to normal cells. AGAP2-AS1 knockdown suppressed cell proliferation, migration, invasion and the epithelial-to-mesenchymal transition process. Interestingly, AGAP2-AS1 sponges miR-4,668-3p to release SRSF1 in CRC. Furthermore, in the rescue functional assay, miR-4,668-3p down-regulation exacerbated the malignant behaviors of AGAP2-AS1-depleted CRC cells, whereas such effects were further offset by SRSF1 knockdown.

CONCLUSIONS

AGAP2-AS1 facilitates cell proliferation, motility and EMT in CRC via targeting the miR-4,668-3p/SRSF1 axis. AGAP2-AS1 or SRSF1 may have potential as underlying therapeutic targets for CRC patients.

Influenza virus natural sequence heterogeneity in segment 8 affects interactions with cellular RNA-binding proteins and splicing efficiency

Segment 8 mRNAs of influenza virus A/Brevig Misson/1918/1 (H1N1) are poorly spliced compared to segment 8 mRNAs of influenza virus A/Netherlands/178/95 (H3N2). Using oligonucleotide-mediated protein pull down with oligos spanning the entire length of segment 8 of either influenza virus H1N1 or influenza virus H3N2 we identified cellular RNA binding proteins that interacted with oligonucleotides derived from either H1N1 or H3N2 sequences. When the identified hot spots for RNA binding proteins in H1N1 segment 8 mRNAs were replaced by H3N2 sequences, splicing efficiency increased significantly. Replacing as few as three nucleotides of the H1N1 mRNA with sequences from H3N2 mRNA, enhanced splicing of the H1N1 mRNAs. Cellular proteins U2AF65 and HuR interacted preferentially with the 3'-splice site of H3N2 and overexpression of HuR reduced the levels of unspliced H1N1 mRNAs, suggesting that U2AF65 and HuR contribute to control of influenza virus mRNA splicing.

Profiling of cis- and trans-acting factors supporting noncanonical splice site activation

Recently, by combining transcriptomics with functional splicing reporter assays we were able to identify GT > GC > TT as the three highest ranked dinucleotides of human 5' splice sites (5'ss). Here, we have extended our investigations to the proteomic characterization of nuclear proteins that bind to canonical and noncanonical 5'ss. Surprisingly, we found that U1 snRNP binding to functional 5'ss sequences prevented components of the DNA damage response (DDR) from binding to the RNA, suggesting a close link between spliceosome arrangement and genome stability. We demonstrate that all tested noncanonical 5'ss sequences are bona-fide targets of the U2-type spliceosome and are bound by U1 snRNP, including U1-C, in the presence of splicing enhancers. The quantity of precipitated U1-C protein was similar for all noncanonical 5'ss dinucleotides, so that the highly different 5'ss usage was likely due to a later step after early U1 snRNP binding. In addition, we show that an internal GT at positions +5/+6 can be advantageous for splicing at position +1 of noncanonical splice sites. Likewise, and in agreement with previous observations, splicing inactive U1 snRNP binding sites could serve as splicing enhancers, which may also explain the higher abundance of U1 snRNPs compared to other U snRNPs. Finally, we observe that an arginine-serine (RS)-rich domain recruitment to stem loop I of the U1 snRNA is functionally sufficient to promote exon-definition and upstream 3'ss activation.

Serine/arginine-rich splicing factors: the bridge linking alternative splicing and cancer

The serine/arginine-rich splicing factors (SRs) belong to the serine arginine-rich protein family, which plays an extremely important role in the splicing process of precursor RNA. The SRs recognize the splicing elements on precursor RNA, then recruit and assemble spliceosome to promote or inhibit the occurrence of splicing events. In tumors, aberrant expression of SRs causes abnormal splicing of RNA, contributing to proliferation, migration and apoptosis resistance of tumor cells. Here, we reviewed the vital role of SRs in various tumors and discussed the promise of analyzing mRNA alternative splicing events in tumor. Further, we highlight the challenges and discussed the perspectives for the identification of new potential targets for cancer therapy via SRs family members.

A large-scale binding and functional map of human RNA-binding proteins

Many proteins regulate the expression of genes by binding to specific regions encoded in the genome1. Here we introduce a new data set of RNA elements in the human genome that are recognized by RNA-binding proteins (RBPs), generated as part of the Encyclopedia of DNA Elements (ENCODE) project phase III. This class of regulatory elements functions only when transcribed into RNA, as they serve as the binding sites for RBPs that control post-transcriptional processes such as splicing, cleavage and polyadenylation, and the editing, localization, stability and translation of mRNAs. We describe the mapping and characterization of RNA elements recognized by a large collection of human RBPs in K562 and HepG2 cells. Integrative analyses using five assays identify RBP binding sites on RNA and chromatin in vivo, the in vitro binding preferences of RBPs, the function of RBP binding sites and the subcellular localization of RBPs, producing 1,223 replicated data sets for 356 RBPs. We describe the spectrum of RBP binding throughout the transcriptome and the connections between these interactions and various aspects of RNA biology, including RNA stability, splicing regulation and RNA localization. These data expand the catalogue of functional elements encoded in the human genome by the addition of a large set of elements that function at the RNA level by interacting with RBPs.

TIMP1 intron 3 retention is a marker of colon cancer progression controlled by hnRNPA1

We previously reported a 40-transcripts signature marking the normal mucosa to colorectal adenocarcinoma transition. Eight of these mRNAs also showed splicing alterations, including a specific intron 3 retention in tissue metalloprotease inhibitor I (TIMP1), which decreased during the early steps of colorectal cancer progression. To decipher the mechanism of intron 3 retention/splicing, we first searched for putative RNA binding protein binding sites onto the TIMP1 sequence. We identified potential serine arginine rich splicing factor 1 (SRSF1) and heterogeneous nuclear RiboNucleoProtein A1 (hnRNPA1) binding sites at the end of intron 3 and the beginning of exon 4, respectively. RNA immunoprecipitation showed that hnRNPA1, but not SRSF1 could bind to the corresponding region in TIMP1 pre-mRNA in live cells. Furthermore, using a TIMP1-based ex vivo minigene approach, together with a plasmon resonance in vitro RNA binding assay, we confirmed that hnRNPA1 could indeed bind to wild type TIMP1 exon 4 pre-mRNA and control TMP1 intron 3 splicing, the interaction being abolished in presence of a mutant sequence that disrupted this site. These results indicated that hnRNPA1, upon binding to TIMP1 exon 4, was a positive regulator of intron 3 splicing. We propose that this TIMP1-intron 3 + transcript belongs to the class of nuclear transcripts with "detained" introns, an abundant molecular class, including in cancer.

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.

Srsf7 Establishes the Juvenile Transcriptome through Age-Dependent Alternative Splicing in Mice

The juvenile phase is characterized by continuously progressing physiological processes such as growth and maturation, which are accompanied by transitions in gene expression. The contribution of transcriptome dynamics to the establishment of juvenile properties remains unclear. Here, we investigated alternative splicing (AS) events in postnatal growth and elucidated the landscape of age-dependent alternative splicing (ADAS) in C57BL/6 mice. Our analysis of ADAS in the cerebral cortex, cardiomyocytes, and hepatocytes revealed numerous juvenile-specific splicing isoforms that shape the juvenile transcriptome, which in turn functions as a basis for the highly anabolic status of juvenile cells. Mechanistically, the juvenile-expressed splicing factor Srsf7 mediates ADAS, as exemplified by switching from juvenile to adult forms of anabolism-associated genes Eif4a2 and Rbm7. Suppression of Srsf7 results in “fast-forwarding” of this transcriptome transition, causing impaired anabolism and growth in mice. Thus, juvenile-specific AS is indispensable for the anabolic state of juveniles and differentiates juveniles from adults.

•

Age-dependent alternative splicing (ADAS) was determined in mice

•

Srsf7 depletion causes loss of cellular juvenescence

•

Srsf7 mutation causes a shift from juvenile to adult-type transcriptome

•

Srsf7 promotes juvenile growth and anabolism through ADAS

A Conserved Kinase-Based Body-Temperature Sensor Globally Controls Alternative Splicing and Gene Expression

Homeothermic organisms maintain their core body temperature in a narrow, tightly controlled range. Whether and how subtle circadian oscillations or disease-associated changes in core body temperature are sensed and integrated in gene expression programs remain elusive. Furthermore, a thermo-sensor capable of sensing the small temperature differentials leading to temperature-dependent sex determination (TSD) in poikilothermic reptiles has not been identified. Here, we show that the activity of CDC-like kinases (CLKs) is highly responsive to physiological temperature changes, which is conferred by structural rearrangements within the kinase activation segment. Lower body temperature activates CLKs resulting in strongly increased phosphorylation of SR proteins in vitro and in vivo. This globally controls temperature-dependent alternative splicing and gene expression, with wide implications in circadian, tissue-specific, and disease-associated settings. This temperature sensor is conserved across evolution and adapted to growth temperatures of diverse poikilotherms. The dynamic temperature range of reptilian CLK homologs suggests a role in TSD.

The SR-protein FgSrp2 regulates vegetative growth, sexual reproduction and pre-mRNA processing by interacting with FgSrp1 in Fusarium graminearum

Serine/arginine (SR) proteins play significant roles in pre-mRNA splicing in eukaryotes. To investigate how gene expression influences fungal development and pathogenicity in Fusarium graminearum, a causal agent of Fusarium head blight (FHB) of wheat and barley, our previous study identified a SR protein FgSrp1 in F. graminearum, and showed that it is important for conidiation, plant infection and pre-mRNA processing. In this study, we identified another SR protein FgSrp2 in F. graminearum, which is orthologous to Schizosaccharomyces pombe Srp2. Our data showed that, whereas yeast Srp2 is essential for growth, deletion of FgSRP2 resulted in only slight defects in vegetative growth and perithecia melanization. FgSrp2 localized to the nucleus and both its N- and C-terminal regions were important for the localization to the nucleus. FgSrp2 interacted with FgSrp1 to form a complex in vivo. Double deletion of FgSRP1 and FgSRP2 revealed that they had overlapping functions in vegetative growth and sexual reproduction. RNA-seq analysis revealed that, although deletion of FgSRP2 alone had minimal effects, deletion of both FgSRP1 and FgSRP2 caused significant changes in gene transcription and RNA splicing. Overall, our results indicated that FgSrp2 regulates vegetative growth, sexual reproduction and pre-mRNA processing by interacting with FgSrp1.

Alternative splicing in aging and longevity

Alternative pre-mRNA splicing increases the complexity of the proteome that can be generated from the available genomic coding sequences. Dysregulation of the splicing process has been implicated in a vast repertoire of diseases. However, splicing has recently been linked to both the aging process itself and pro-longevity interventions. This review focuses on recent research towards defining RNA splicing as a new hallmark of aging. We highlight dysfunctional alternative splicing events that contribute to the aging phenotype across multiple species, along with recent efforts toward deciphering mechanistic roles for RNA splicing in the regulation of aging and longevity. Further, we discuss recent research demonstrating a direct requirement for specific splicing factors in pro-longevity interventions, and specifically how nutrient signaling pathways interface to splicing factor regulation and downstream splicing targets. Finally, we review the emerging potential of using splicing profiles as a predictor of biological age and life expectancy. Understanding the role of RNA splicing components and downstream targets altered in aging may provide opportunities to develop therapeutics and ultimately extend healthy lifespan in humans.

A perspective on crassulacean acid metabolism photosynthesis evolution of orchids on different continents: Dendrobium as a case study

Members of the Orchidaceae, one of the largest families of flowering plants, evolved the crassulacean acid metabolism (CAM) photosynthesis strategy. It is thought that CAM triggers adaptive radiation into new niche spaces, yet very little is known about its origin and diversification on different continents. Here, we assess the prevalence of CAM in Dendrobium, which is one of the largest genera of flowering plants and found in a wide range of environments, from the high altitudes of the Himalayas to relatively arid habitats in Australia. Based on phylogenetic time trees, we estimated that CAM, as determined by δ 13C values less negative than -20.0‰, evolved independently at least eight times in Dendrobium. The oldest lineage appeared in the Asian clade during the middle Miocene, indicating the origin of CAM was associated with a pronounced climatic cooling that followed a period of aridity. Divergence of the four CAM lineages in the Asian clade appeared to be earlier than divergence of those in the Australasian clade. However, CAM species in the Asian clade are much less diverse (25.6%) than those in the Australasian clade (57.9%). These findings shed new light on CAM evolutionary history and the aridity levels of the paleoclimate on different continents.

Mechanisms and Functions of Long Non-Coding RNAs at Multiple Regulatory Levels

Long non-coding (lnc) RNAs are non-coding RNAs longer than 200 nt. lncRNAs primarily interact with mRNA, DNA, protein, and miRNA and consequently regulate gene expression at the epigenetic, transcriptional, post-transcriptional, translational, and post-translational levels in a variety of ways. They play important roles in biological processes such as chromatin remodeling, transcriptional activation, transcriptional interference, RNA processing, and mRNA translation. lncRNAs have important functions in plant growth and development; biotic and abiotic stress responses; and in regulation of cell differentiation, the cell cycle, and the occurrence of many diseases in humans and animals. In this review, we summarize the functions and mechanisms of lncRNAs in plants, humans, and animals at different regulatory levels.

Targeting mRNA processing as an anticancer strategy

Discoveries in the past decade have highlighted the potential of mRNA as a therapeutic target for cancer. Specifically, RNA sequencing revealed that, in addition to gene mutations, alterations in mRNA can contribute to the initiation and progression of cancer. Indeed, precursor mRNA processing, which includes the removal of introns by splicing and the formation of 3' ends by cleavage and polyadenylation, is frequently altered in tumours. These alterations result in numerous cancer-specific mRNAs that generate altered levels of normal proteins or proteins with new functions, leading to the activation of oncogenes or the inactivation of tumour-suppressor genes. Abnormally spliced and polyadenylated mRNAs are also associated with resistance to cancer treatment and, unexpectedly, certain cancers are highly sensitive to the pharmacological inhibition of splicing. This Review summarizes recent progress in our understanding of how splicing and polyadenylation are altered in cancer and highlights how this knowledge has been translated for drug discovery, resulting in the production of small molecules and oligonucleotides that modulate the spliceosome and are in clinical trials for the treatment of cancer.

Coordinated Regulation of Pre-mRNA Splicing by the SFPS-RRC1 Complex to Promote Photomorphogenesis

Light signals perceived by the phytochrome (phy) family of photoreceptors control gene expression at both transcriptional and posttranscriptional levels to promote photomorphogenesis. Recently, we identified a factor called SPLICING FACTOR FOR PHYTOCHROME SIGNALING (SFPS) that directly interacts with the photoreceptor phyB and regulates pre-mRNA splicing in Arabidopsis (Arabidopsis thaliana). To identify SFPS-interacting proteins, we performed an immunoprecipitation followed by a mass spectrometry and identified the Ser/Arg-like protein REDUCED RED-LIGHT RESPONSES IN CRY1CRY2 BACKGROUND1 (RRC1). Genetic analyses revealed that the sfps-2 rrc1-3 phenotypes are similar to those of the single mutants, suggesting that RRC1 and SFPS might function together. RNA sequence analyses of rrc1-3 identified a large number of genes whose pre-mRNA splicing is altered under dark and light conditions. Comparison of the sequence data revealed a subset of common genes coregulated by SFPS and RRC1 under dark and light conditions. Similar to SFPS, RRC1 also interacts with phyB, colocalizes in nuclear photobodies, and regulates light-dependent pre-mRNA splicing of a subset of genes. Taken together, these data suggest that although SFPS and RRC1 can regulate distinct subsets of genes, they also form a complex and coordinately control pre-mRNA splicing of a subset of genes involved in light signaling and circadian clock pathways to promote photomorphogenesis.

CAGI 5 splicing challenge: Improved exon skipping and intron retention predictions with MMSplice

Pathogenic genetic variants often primarily affect splicing. However, it remains difficult to quantitatively predict whether and how genetic variants affect splicing. In 2018, the fifth edition of the Critical Assessment of Genome Interpretation proposed two splicing prediction challenges based on experimental perturbation assays: Vex-seq, assessing exon skipping, and MaPSy, assessing splicing efficiency. We developed a modular modeling framework, MMSplice, the performance of which was among the best on both challenges. Here we provide insights into the modeling assumptions of MMSplice and its individual modules. We furthermore illustrate how MMSplice can be applied in practice for individual genome interpretation, using the MMSplice VEP plugin and the Kipoi variant interpretation plugin, which are directly applicable to VCF files.

How RNA structure dictates the usage of a critical exon of spinal muscular atrophy gene

Role of RNA structure in pre-mRNA splicing has been implicated for several critical exons associated with genetic disorders. However, much of the structural studies linked to pre-mRNA splicing regulation are limited to terminal stem-loop structures (hairpins) sequestering splice sites. In few instances, role of long-distance interactions is implicated as the major determinant of splicing regulation. With the recent surge of reports of circular RNA (circRNAs) generated by backsplicing, role of Alu-associated RNA structures formed by long-range interactions are taking central stage. Humans contain two nearly identical copies of Survival Motor Neuron (SMN) genes, SMN1 and SMN2. Deletion or mutation of SMN1 coupled with the inability of SMN2 to compensate for the loss of SMN1 due to exon 7 skipping causes spinal muscular atrophy (SMA), one of the leading genetic diseases of children. In this review, we describe how structural elements formed by both local and long-distance interactions are being exploited to modulate SMN2 exon 7 splicing as a potential therapy for SMA. We also discuss how Alu-associated secondary structure modulates generation of a vast repertoire of SMN circRNAs. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Combinatorial regulation of alternative splicing

The generation of protein coding mRNAs from pre-mRNA is a fundamental biological process that is required for gene expression. Alternative pre-mRNA splicing is responsible for much of the transcriptomic and proteomic diversity observed in higher order eukaryotes. Aberrations that disrupt regular alternative splicing patterns are known to cause human diseases, including various cancers. Alternative splicing is a combinatorial process, meaning many factors affect which two splice sites are ligated together. The features that dictate exon inclusion are comprised of splice site strength, intron-exon architecture, RNA secondary structure, splicing regulatory elements, promoter use and transcription speed by RNA polymerase and the presence of post-transcriptional nucleotide modifications. A comprehensive view of all of the factors that influence alternative splicing decisions is necessary to predict splicing outcomes and to understand the molecular basis of disease. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

A mechanism underlying position-specific regulation of alternative splicing

Many RNA-binding proteins including a master regulator of splicing in developing brain and muscle, polypyrimidine tract-binding protein 1 (PTBP1), can either activate or repress alternative exons depending on the pre-mRNA recruitment position. When bound upstream or within regulated exons PTBP1 tends to promote their skipping, whereas binding to downstream sites often stimulates inclusion. How this switch is orchestrated at the molecular level is poorly understood. Using bioinformatics and biochemical approaches we show that interaction of PTBP1 with downstream intronic sequences can activate natural cassette exons by promoting productive docking of the spliceosomal U1 snRNP to a suboptimal 5' splice site. Strikingly, introducing upstream PTBP1 sites to this circuitry leads to a potent splicing repression accompanied by the assembly of an exonic ribonucleoprotein complex with a tightly bound U1 but not U2 snRNP. Our data suggest a molecular mechanism underlying the transition between a better-known repressive function of PTBP1 and its role as a bona fide splicing activator. More generally, we argue that the functional outcome of individual RNA contacts made by an RNA-binding protein is subject to extensive context-specific modulation.

Long noncoding RNA HOXA11-AS promotes gastric cancer cell proliferation and invasion via SRSF1 and functions as a biomarker in gastric cancer

BACKGROUND

Gastric cancer (GC) is the fourth most frequent malignancy all over the world. The diagnosis of GC is challenging and the prognosis of GC is very unfavorable. Accumulating evidence reveals that serum long noncoding RNAs (lncRNAs) can function as biomarkers in various types of cancers, including GC.

AIM

To explore the level and molecular mechanism of the lncRNA HOXA11-AS in GC and the diagnostic and prognostic significance of serum HOXA11-AS in GC.

METHODS

HOXA11-AS levels in GC tissue, cell lines, and serum samples were measured. The correlation between HOXA11-AS expression and clinicopathological characteristics was analyzed. The role of HOXA11-AS in the diagnosis and prognosis of GC was evaluated. Cell function assays were performed for exploration of the roles of HOXA11-AS in GC cells. Moreover, Western blot was performed to explore the target regulated by HOXA11-AS in GC cells.

RESULTS

Up-regulation of HOXA11-AS was found in GC tissues, cell lines, and serum samples. In GC patients, decreased serum HOXA11-AS levels were negatively related with tumor size, TNM stage, and lymph node metastasis. The area under the receiver operating characteristic curve of serum HOXA11-AS in the diagnosis of GC was 0.924 (95%CI: 0.881-0.967; sensitivity, 0.787; specificity 0.978). Results of the Kaplan-Meier survival curves suggested the GC patients with a lower HOXA11-AS level having a better overall survival rate. HOXA11-AS promoted GC cell proliferation and invasion. SRSF1 may be the target regulated by HOXA11-AS in GC cells.

CONCLUSION

HOXA11-AS promotes GC cell proliferation and invasion via SRSF1 and may function as a promising marker in GC.

Context matters: Regulation of splice donor usage

Elaborate research on splicing, starting in the late seventies, evolved from the discovery that 5' splice sites are recognized by their complementarity to U1 snRNA towards the realization that RNA duplex formation cannot be the sole basis for 5'ss selection. Rather, their recognition is highly influenced by a number of context factors including transcript architecture as well as splicing regulatory elements (SREs) in the splice site neighborhood. In particular, proximal binding of splicing regulatory proteins highly influences splicing outcome. The importance of SRE integrity especially becomes evident in the light of human pathogenic mutations where single nucleotide changes in SREs can severely affect the resulting transcripts. Bioinformatics tools nowadays greatly assist in the computational evaluation of 5'ss, their neighborhood and the impact of pathogenic mutations. Although predictions are already quite robust, computational evaluation of the splicing regulatory landscape still faces challenges to increase future reliability. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Alternative splicing coupled to nonsense-mediated mRNA decay contributes to the high-altitude adaptation of maca (Lepidium meyenii)

Alpine plants remain the least studied plant communities in terrestrial ecosystems. However, how they adapt to high-altitude environments is far from clear. Here, we used RNA-seq to investigate a typical alpine plant maca (Lepidium meyenii) to understand its high-altitude adaptation at transcriptional and post-transcriptional level. At transcriptional level, we found that maca root significantly up-regulated plant immunity genes in day-time comparing to night-time, and up-regulated abiotic (cold/osmotic) stress response genes in Nov and Dec comparing to Oct. In addition, 17 positively selected genes were identified, which could be involved in mitochondrion. At post-transcriptional level, we found that maca had species-specific characterized alternative splicing (AS) profile which could be influenced by stress environments. For example, the alternative 3' splice site events (A3SS, 39.62%) were predominate AS events in maca, rather than intron retention (IR, 23.17%). Interestingly, besides serine/arginine-rich (SR) proteins and long non-coding RNAs (lncRNAs), a lot of components in nonsense-mediated mRNA decay (NMD) were identified under differential alternative splicing (DAS), supporting AS coupled to NMD as essential mechanisms for maca's stress responses and high-altitude adaptation. Taken together, we first attempted to unveil maca's high-altitude adaptation mechanisms based on transcriptome and post-transcriptome evidence. Our data provided valuable insights to understand the high-altitude adaptation of alpine plants.

A novel role of U1 snRNP: Splice site selection from a distance

Removal of introns by pre-mRNA splicing is fundamental to gene function in eukaryotes. However, understanding the mechanism by which exon-intron boundaries are defined remains a challenging endeavor. Published reports support that the recruitment of U1 snRNP at the 5'ss marked by GU dinucleotides defines the 5'ss as well as facilitates 3'ss recognition through cross-exon interactions. However, exceptions to this rule exist as U1 snRNP recruited away from the 5'ss retains the capability to define the splice site, where the cleavage takes place. Independent reports employing exon 7 of Survival Motor Neuron (SMN) genes suggest a long-distance effect of U1 snRNP on splice site selection upon U1 snRNP recruitment at target sequences with or without GU dinucleotides. These findings underscore that sequences distinct from the 5'ss may also impact exon definition if U1 snRNP is recruited to them through partial complementarity with the U1 snRNA. In this review we discuss the expanded role of U1 snRNP in splice-site selection due to U1 ability to be recruited at more sites than predicted solely based on GU dinucleotides.

Interplay between coding and exonic splicing regulatory sequences

The inclusion of exons during the splicing process depends on the binding of splicing factors to short low-complexity regulatory sequences. The relationship between exonic splicing regulatory sequences and coding sequences is still poorly understood. We demonstrate that exons that are coregulated by any given splicing factor share a similar nucleotide composition bias and preferentially code for amino acids with similar physicochemical properties because of the nonrandomness of the genetic code. Indeed, amino acids sharing similar physicochemical properties correspond to codons that have the same nucleotide composition bias. In particular, we uncover that the TRA2A and TRA2B splicing factors that bind to adenine-rich motifs promote the inclusion of adenine-rich exons coding preferentially for hydrophilic amino acids that correspond to adenine-rich codons. SRSF2 that binds guanine/cytosine-rich motifs promotes the inclusion of GC-rich exons coding preferentially for small amino acids, whereas SRSF3 that binds cytosine-rich motifs promotes the inclusion of exons coding preferentially for uncharged amino acids, like serine and threonine that can be phosphorylated. Finally, coregulated exons encoding amino acids with similar physicochemical properties correspond to specific protein features. In conclusion, the regulation of an exon by a splicing factor that relies on the affinity of this factor for specific nucleotide(s) is tightly interconnected with the exon-encoded physicochemical properties. We therefore uncover an unanticipated bidirectional interplay between the splicing regulatory process and its biological functional outcome.

Gymnotic Delivery of LNA Mixmers Targeting Viral SREs Induces HIV-1 mRNA Degradation

Transcription of the HIV-1 provirus generates a viral pre-mRNA, which is alternatively spliced into more than 50 HIV-1 mRNAs encoding all viral proteins. Regulation of viral alternative splice site usage includes the presence of splicing regulatory elements (SREs) which can dramatically impact RNA expression and HIV-1 replication when mutated. Recently, we were able to show that two viral SREs, G_I3-2 and ESE_tat, are important players in the generation of viral vif, vpr and tat mRNAs. Furthermore, we demonstrated that masking these SREs by transfected locked nucleic acid (LNA) mixmers affect the viral splicing pattern and viral particle production. With regard to the development of future therapeutic LNA mixmer-based antiretroviral approaches, we delivered the G_I3-2 and the ESE_tat LNA mixmers “nakedly”, without the use of transfection reagents (gymnosis) into HIV-1 infected cells. Surprisingly, we observed that gymnotically-delivered LNA mixmers accumulated in the cytoplasm, and seemed to co-localize with GW bodies and induced degradation of mRNAs containing their LNA target sequence. The G_I3-2 and the ESE_tat LNA-mediated RNA degradation resulted in abrogation of viral replication in HIV-1 infected Jurkat and PM1 cells as well as in PBMCs.

Splicing of long non-coding RNAs primarily depends on polypyrimidine tract and 5' splice-site sequences due to weak interactions with SR proteins

Many nascent long non-coding RNAs (lncRNAs) undergo the same maturation steps as pre-mRNAs of protein-coding genes (PCGs), but they are often poorly spliced. To identify the underlying mechanisms for this phenomenon, we searched for putative splicing inhibitory sequences using the ncRNA-a2 as a model. Genome-wide analyses of intergenic lncRNAs (lincRNAs) revealed that lincRNA splicing efficiency positively correlates with 5'ss strength while no such correlation was identified for PCGs. In addition, efficiently spliced lincRNAs have higher thymidine content in the polypyrimidine tract (PPT) compared to efficiently spliced PCGs. Using model lincRNAs, we provide experimental evidence that strengthening the 5'ss and increasing the T content in PPT significantly enhances lincRNA splicing. We further showed that lincRNA exons contain less putative binding sites for SR proteins. To map binding of SR proteins to lincRNAs, we performed iCLIP with SRSF2, SRSF5 and SRSF6 and analyzed eCLIP data for SRSF1, SRSF7 and SRSF9. All examined SR proteins bind lincRNA exons to a much lower extent than expression-matched PCGs. We propose that lincRNAs lack the cooperative interaction network that enhances splicing, which renders their splicing outcome more dependent on the optimality of splice sites.

View from an mRNP: The Roles of SR Proteins in Assembly, Maturation and Turnover

Serine- and arginine-rich proteins (SR proteins) are a family of multitasking RNA-binding proteins (RBPs) that are key determinants of messenger ribonucleoprotein (mRNP) formation, identity and fate. Apart from their essential functions in pre-mRNA splicing, SR proteins display additional pre- and post-splicing activities and connect nuclear and cytoplasmic gene expression machineries. Through changes in their post-translational modifications (PTMs) and their subcellular localization, they provide functional specificity and adjustability to mRNPs. Transcriptome-wide UV crosslinking and immunoprecipitation (CLIP-Seq) studies revealed that individual SR proteins are present in distinct mRNPs and act in specific pairs to regulate different gene expression programmes. Adopting an mRNP-centric viewpoint, we discuss the roles of SR proteins in the assembly, maturation, quality control and turnover of mRNPs and describe the mechanisms by which they integrate external signals, coordinate their multiple tasks and couple subsequent mRNA processing steps.

SRSF1 modulates PTPMT1 alternative splicing to regulate lung cancer cell radioresistance

BACKGROUND

Radioresistance is the major cause of cancer treatment failure. Additionally, splicing dysregulation plays critical roles in tumorigenesis. However, the involvement of alternative splicing in resistance of cancer cells to radiotherapy remains elusive. We sought to investigate the key role of the splicing factor SRSF1 in the radioresistance in lung cancer.

METHODS

Lung cancer cell lines, xenograft mice models, and RNA-seq were employed to study the detailed mechanisms of SRSF1 in lung cancer radioresistance. Clinical tumor tissues and TCGA dataset were utilized to determine the expression levels of distinct SRSF1-regulated splicing isoforms. KM-plotter was applied to analyze the survival of cancer patients with various levels of SRSF1-regulated splicing isoforms.

FINDINGS

Splicing factors were screened to identify their roles in radioresistance, and SRSF1 was found to be involved in radioresistance in cancer cells. The level of SRSF1 is elevated in irradiation treated lung cancer cells, whereas knockdown of SRSF1 sensitizes cancer cells to irradiation. Mechanistically, SRSF1 modulates various cancer-related splicing events, particularly the splicing of PTPMT1, a PTEN-like mitochondrial phosphatase. Reduced SRSF1 favors the production of short isoforms of PTPMT1 upon irradiation, which in turn promotes phosphorylation of AMPK, thereby inducing DNA double-strand break to sensitize cancer cells to irradiation. Additionally, the level of the short isoform of PTPMT1 is decreased in cancer samples, which is correlated to cancer patients' survival.

CONCLUSIONS

Our study provides mechanistic analyses of aberrant splicing in radioresistance in lung cancer cells, and establishes SRSF1 as a potential therapeutic target for sensitization of patients to radiotherapy.

Ranking noncanonical 5' splice site usage by genome-wide RNA-seq analysis and splicing reporter assays

Most human pathogenic mutations in 5' splice sites affect the canonical GT in positions +1 and +2, leading to noncanonical dinucleotides. On the other hand, noncanonical dinucleotides are observed under physiological conditions in ∼1% of all human 5'ss. It is therefore a challenging task to understand the pathogenic mutation mechanisms underlying the conditions under which noncanonical 5'ss are used. In this work, we systematically examined noncanonical 5' splice site selection, both experimentally using splicing competition reporters and by analyzing a large RNA-seq data set of 54 fibroblast samples from 27 subjects containing a total of 2.4 billion gapped reads covering 269,375 exon junctions. From both approaches, we consistently derived a noncanonical 5'ss usage ranking GC > TT > AT > GA > GG > CT. In our competition splicing reporter assay, noncanonical splicing was strictly dependent on the presence of upstream or downstream splicing regulatory elements (SREs), and changes in SREs could be compensated by variation of U1 snRNA complementarity in the competing 5'ss. In particular, we could confirm splicing at different positions (i.e., -1, +1, +5) of a splice site for all noncanonical dinucleotides "weaker" than GC. In our comprehensive RNA-seq data set analysis, noncanonical 5'ss were preferentially detected in weakly used exon junctions of highly expressed genes. Among high-confidence splice sites, they were 10-fold overrepresented in clusters with a neighboring, more frequently used 5'ss. Conversely, these more frequently used neighbors contained only the dinucleotides GT, GC, and TT, in accordance with the above ranking.

Characterization of cis-acting elements that control oscillating alternative splicing

Alternative splicing (AS) in response to changing external conditions often requires alterations in the ability of sequence-specific RNA-binding proteins to bind to cis-acting sequences in their target pre-mRNA. While daily oscillations in AS events have been described in several organisms, cis-acting sequences that control time of the day-dependent AS remain largely elusive. Here we define cis-regulatory RNA elements that control body-temperature driven rhythmic AS using the mouse U2af26 gene as a model system. We identify a complex network of cis-regulatory sequences that regulate AS of U2af26, and show that the activity of two enhancer elements is necessary for oscillating AS. A minigene comprising these U2af26 regions recapitulates rhythmic splicing of the endogenous gene, which is controlled through temperature-regulated SR protein phosphorylation. Mutagenesis of the minigene delineates the cis-acting enhancer element for SRSF2 within exon 6 to single nucleotide resolution and reveals that the combined activity of SRSF2 and SRSF7 is required for oscillating U2af26 AS. By combining RNA-Seq with an siRNA screen and individual-nucleotide resolution cross-linking and immunoprecipitation (iCLIP), we identify a complex network of SR proteins that globally controls temperature-dependent rhythmic AS, with the direction of splicing depending on the position of the cis-acting elements. Together, we provide detailed insights into the sequence requirements that allow trans-acting factors to generate daily rhythms in AS.

HMGA1a induces alternative splicing of estrogen receptor alpha in MCF-7 human breast cancer cells

The high-mobility group A protein 1a (HMGA1a) protein is known as an oncogene whose expression level in cancer tissue correlates with the malignant potential, and known as a component of senescence-related structures connecting it to tumor suppressor networks in fibroblasts. HMGA1 protein binds to DNA, but recent studies have shown it exerts novel functions through RNA-binding. Our previous studies have shown that sequence-specific RNA-binding of HMGA1a induces exon-skipping of Presenilin-2 exon 5 in sporadic Alzheimer disease. Here we show that HMGA1a induced exon-skipping of the estrogen receptor alpha (ERα) gene and increased ERα46 mRNA expression in MCF-7 breast cancer cells. An RNA-decoy of HMGA1a efficiently blocked this event and reduced ERα46 protein expression. Blockage of HMGA1a RNA-binding property consequently induced cell growth through reduced ERα46 expression in MCF-7 cells and increased sensitivity to tamoxifen in the tamoxifen-resistant cell line, MCF-7/TAMR1. Stable expression of an HMGA1a RNA-decoy in MCF-7 cells exhibited decreased ERα46 protein expression and increased estrogen-dependent tumor growth when these cells were implanted in nude mice. These results show HMGA1a is involved in alternative splicing of the ERα gene and related to estrogen-related growth as well as tamoxifen sensitivity in MCF-7 breast cancer cells.

Genome-Wide Association Analyses Reveal the Importance of Alternative Splicing in Diversifying Gene Function and Regulating Phenotypic Variation in Maize

Alternative splicing (AS) enhances transcriptome diversity and plays important roles in regulating plant processes. Although widespread natural variation in AS has been observed in plants, how AS is regulated and contribute to phenotypic variation is poorly understood. Here, we report a population-level transcriptome assembly and genome-wide association study to identify splicing quantitative trait loci (sQTLs) in developing maize (Zea mays) kernels from 368 inbred lines. We detected 19,554 unique sQTLs for 6570 genes. Most sQTLs showed small isoform usage changes without involving major isoform switching between genotypes. The sQTL-affected isoforms tend to display distinct protein functions. We demonstrate that nonsense-mediated mRNA decay, microRNA-mediated regulation, and small interfering peptide-mediated peptide interference are frequently involved in sQTL regulation. The natural variation in AS and overall mRNA level appears to be independently regulated with different cis-sequences preferentially used. We identified 214 putative trans-acting splicing regulators, among which ZmGRP1, encoding an hnRNP-like glycine-rich RNA binding protein, regulates the largest trans-cluster. Knockout of ZmGRP1 by CRISPR/Cas9 altered splicing of numerous downstream genes. We found that 739 sQTLs colocalized with previous marker-trait associations, most of which occurred without changes in overall mRNA level. Our findings uncover the importance of AS in diversifying gene function and regulating phenotypic variation.

Growth Performance, Serum Biochemical Profile, Oxidative Status, and Fertility Traits in Male Japanese Quail Fed on Ginger (Zingiber officinale, Roscoe) Essential Oil

This study was designed to evaluate the effect of ginger (Zingiber officinale, Rosc.) essential oil on growth performance, serum biochemical profile, oxidative stress, and histological structure of testes and fertility traits in Japanese quail (Coturnix coturnix japonica). 96 three-week-old male Japanese quail weighing between 120 and 130 g were randomly assigned to 4 dietary treatment groups in a completely randomized design. Each group was divided into 4 replicates of 6 quails. Quails in control group received orally 100 μl/kg bw of distilled water, while the three test groups received, respectively, by gastric intubation 50, 100, and 150 μl/kg bw of ginger essential oil. At 12 weeks old, twelve birds per treatment were randomly selected and fasted for 24 hours, weighed, and slaughtered to assess organ and biochemical parameters. At the same period, 4 mature male quails per treatment were chosen at random and individually housed in cages, each with four untreated females for fertility and hatchability traits. The main results revealed that growth characteristics were not markedly (P > 0.05) affected by essential oil whatever the dose. The left testis weight increased significantly (P < 0.05) with 100 and 150 μl/kg bw of essential oil compared to the control. The serum content in total cholesterol and triglycerides, the liver weight, the serum content in transaminases, and malondialdehyde decreased in treated quails. The serum content in total protein and globulin and the antioxidant enzymes activities increased in treated birds compared to the control. The histological changes in the testis were less visible in treated Japanese quails. At the doses of 100 and 150 μl/kg bw, this essential oil induced a significant increase (P < 0.05) in fertility rate compared to the control. Under the conditions of this study, the ginger rhizomes essential oil can be used in poultry to reduce the lipid peroxidation in reproductive tissues and improve the fertility traits.

Mutually exclusive acetylation and ubiquitylation of the splicing factor SRSF5 control tumor growth

Most tumor cells take up more glucose than normal cells. Splicing dysregulation is one of the molecular hallmarks of cancer. However, the role of splicing factor in glucose metabolism and tumor development remains poorly defined. Here, we show that upon glucose intake, the splicing factor SRSF5 is specifically induced through Tip60-mediated acetylation on K125, which antagonizes Smurf1-mediated ubiquitylation. SRSF5 promotes the alternative splicing of CCAR1 to produce CCAR1S proteins, which promote tumor growth by enhancing glucose consumption and acetyl-CoA production. Conversely, upon glucose starvation, SRSF5 is deacetylated by HDAC1, and ubiquitylated by Smurf1 on the same lysine, resulting in proteasomal degradation of SRSF5. The CCAR1L proteins accumulate to promote apoptosis. Importantly, SRSF5 is hyperacetylated and upregulated in human lung cancers, which correlates with increased CCAR1S expression and tumor progression. Thus, SRSF5 responds to high glucose to promote cancer development, and SRSF5-CCAR1 axis may be valuable targets for cancer therapeutics.

Genetics of alternative splicing evolution during sunflower domestication

Alternative splicing enables organisms to produce the diversity of proteins necessary for multicellular life by using relatively few protein-coding genes. Although differences in splicing have been identified among divergent taxa, the shorter-term evolution of splicing is understudied. The origins of novel splice forms, and the contributions of alternative splicing to major evolutionary transitions, are largely unknown. This study used transcriptomes of wild and domesticated sunflowers to examine splice differentiation and regulation during domestication. We identified substantial splicing divergence between wild and domesticated sunflowers, mainly in the form of intron retention. Transcripts with divergent splicing were enriched for seed-development functions, suggesting that artificial selection impacted splicing patterns. Mapping of quantitative trait loci (QTLs) associated with 144 differential splicing cases revealed primarily trans-acting variation affecting splicing patterns. A large proportion of identified QTLs contain known spliceosome proteins and are associated with splicing variation in multiple genes. Examining a broader set of wild and domesticated sunflower genotypes revealed that most differential splicing patterns in domesticated sunflowers likely arose from standing variation in wild Helianthus annuus and gained frequency during the domestication process. However, several domesticate-associated splicing patterns appear to be introgressed from other Helianthus species. These results suggest that sunflower domestication involved selection on pleiotropic regulatory alleles. More generally, our findings indicate that substantial differences in isoform abundances arose rapidly during a recent evolutionary transition and appear to contribute to adaptation and population divergence.

Blocking the 5' splice site of exon 4 by a morpholino oligomer triggers APOL1 protein isoform switch

APOL1 risk alleles G1 or G2 are associated with a kidney disease phenotype exclusively in people of recent African ancestry. Here we show that exon 4 encoding a part of the APOL1 signal peptide is constitutively spliced in major APOL1 transcripts expressed in kidney glomerular and tubular cells. We demonstrate that constitutive splicing of exon 4 results from a suboptimal hnRNP A1 binding motif found in exon 4. Accordingly, a robust binding of hnRNP A1 protein to a consensus hnRNP A1 cis-acting element in exon 4 results in almost complete exclusion of exon 4 from the APOL1 minigene transcripts. Blocking the 5' splice site at the exon 4/intron boundary with a specific antisense morpholino oligonucleotide excludes exon 4 from the splicing pattern of endogenous APOL1 transcripts. These transcripts are fully functional and produce APOL1 protein isoform that is not normally detectable in podocytes. Together with our previous data showing no cytotoxicity of overexpressed APOL1 isoform lacking exon 4, we propose that morpholino-induced APOL1 isoform switch may provide a new tool to identify in vivo molecular mechanism(s) by which risk alleles promote or mediate the kidney disease phenotype.

Behind the scenes of HIV-1 replication: Alternative splicing as the dependency factor on the quiet

Alternative splicing plays a key role in the HIV-1 life cycle and is essential to maintain an equilibrium of mRNAs that encode viral proteins and polyprotein-isoforms. In particular, since all early HIV-1 proteins are expressed from spliced intronless and late enzymatic and structural proteins from intron containing, i.e. splicing repressed viral mRNAs, cellular splicing factors and splicing regulatory proteins are crucial for the replication capacity. In this review, we will describe the complex network of cis-acting splicing regulatory elements (SREs), which are mainly localized in the neighbourhoods of all HIV-1 splice sites and warrant the proper ratio of individual transcript isoforms. Since SREs represent binding sites for trans-acting cellular splicing factors interacting with the cellular spliceosomal apparatus we will review the current knowledge of interactions between viral RNA and cellular proteins as well as their impact on viral replication. Finally, we will discuss potential therapeutic approaches targeting HIV-1 alternative splicing.

SRSF2 directly inhibits intron splicing to suppresses cassette exon inclusion

SRSF2, a Serine-Arginine rich (SR) protein, is a splicing activator that mediates exon inclusion and exclusion events equally well. Here we show SRSF2 directly suppresses intron splicing to suppress cassette exon inclusion in SMN pre-mRNA. Through a serial mutagenesis, we demonstrate that a 10 nt RNA sequence surrounding the branch-point (BP), is important for SRSF2-mediated inhibition of cassette exon inclusion through directly interacting with SRSF2. We conclude that SRSF2 inhibits intron splicing to promote exon exclusion.