A novel approach to describe a U1 snRNA binding site

Expanding the genetic and phenotypic spectrum of congenital myasthenic syndrome: new homozygous VAMP1 splicing variants in 2 novel individuals

We report the cases of two Spanish pediatric patients with hypotonia, muscle weakness and feeding difficulties at birth. Whole-exome sequencing (WES) uncovered two new homozygous VAMP1 (Vesicle Associated Membrane Protein 1) splicing variants, NM_014231.5:c.129+5 G > A in the boy patient (P1) and c.341-24_341-16delinsAGAAAA in the girl patient (P2). This gene encodes the vesicle-associated membrane protein 1 (VAMP1) that is a component of a protein complex involved in the fusion of synaptic vesicles with the presynaptic membrane. VAMP1 has a highly variable C-terminus generated by alternative splicing that gives rise to three main isoforms (A, B and D), being VAMP1A the only isoform expressed in the nervous system. In order to assess the pathogenicity of these variants, expression experiments of RNA for VAMP1 were carried out. The c.129+5 G > A and c.341-24_341-16delinsAGAAAA variants induced aberrant splicing events resulting in the deletion of exon 2 (r.5_131del; p.Ser2TrpfsTer7) in the three isoforms in the first case, and the retention of the last 14 nucleotides of the 3' of intron 4 (r.340_341ins341-14_341-1; p.Ile114AsnfsTer77) in the VAMP1A isoform in the second case. Pathogenic VAMP1 variants have been associated with autosomal dominant spastic ataxia 1 (SPAX1) and with autosomal recessive presynaptic congenital myasthenic syndrome (CMS). Our patients share the clinical manifestations of CMS patients with two important differences: they do not show the typical electrophysiological pattern that suggests pathology of pre-synaptic neuromuscular junction, and their muscular biopsies present hypertrophic fibers type 1. In conclusion, our data expand both genetic and phenotypic spectrum associated with VAMP1 variants.

Srsf1 and Elavl1 act antagonistically on neuronal fate choice in the developing neocortex by controlling TrkC receptor isoform expression

The seat of higher-order cognitive abilities in mammals, the neocortex, is a complex structure, organized in several layers. The different subtypes of principal neurons are distributed in precise ratios and at specific positions in these layers and are generated by the same neural progenitor cells (NPCs), steered by a spatially and temporally specified combination of molecular cues that are incompletely understood. Recently, we discovered that an alternatively spliced isoform of the TrkC receptor lacking the kinase domain, TrkC-T1, is a determinant of the corticofugal projection neuron (CFuPN) fate. Here, we show that the finely tuned balance between TrkC-T1 and the better known, kinase domain-containing isoform, TrkC-TK+, is cell type-specific in the developing cortex and established through the antagonistic actions of two RNA-binding proteins, Srsf1 and Elavl1. Moreover, our data show that Srsf1 promotes the CFuPN fate and Elavl1 promotes the callosal projection neuron (CPN) fate in vivo via regulating the distinct ratios of TrkC-T1 to TrkC-TK+. Taken together, we connect spatio-temporal expression of Srsf1 and Elavl1 in the developing neocortex with the regulation of TrkC alternative splicing and transcript stability and neuronal fate choice, thus adding to the mechanistic and functional understanding of alternative splicing in vivo.

Alu RNA fold links splicing with signal recognition particle proteins

Transcriptomic diversity in primates was considerably expanded by exonizations of intronic Alu elements. To better understand their cellular mechanisms we have used structure-based mutagenesis coupled with functional and proteomic assays to study the impact of successive primate mutations and their combinations on inclusion of a sense-oriented AluJ exon in the human F8 gene. We show that the splicing outcome was better predicted by consecutive RNA conformation changes than by computationally derived splicing regulatory motifs. We also demonstrate an involvement of SRP9/14 (signal recognition particle) heterodimer in splicing regulation of Alu-derived exons. Nucleotide substitutions that accumulated during primate evolution relaxed the conserved left-arm AluJ structure including helix H1 and reduced the capacity of SRP9/14 to stabilize the closed Alu conformation. RNA secondary structure-constrained mutations that promoted open Y-shaped conformations of the Alu made the Alu exon inclusion reliant on DHX9. Finally, we identified additional SRP9/14 sensitive Alu exons and predicted their functional roles in the cell. Together, these results provide unique insights into architectural elements required for sense Alu exonization, identify conserved pre-mRNA structures involved in exon selection and point to a possible chaperone activity of SRP9/14 outside the mammalian signal recognition particle.

Splicing accuracy varies across human introns, tissues and age

Alternative splicing impacts most multi-exonic human genes. Inaccuracies during this process may have an important role in ageing and disease. Here, we investigated mis-splicing using RNA-sequencing data from ~14K control samples and 42 human body sites, focusing on split reads partially mapping to known transcripts in annotation. We show that mis-splicing occurs at different rates across introns and tissues and that these splicing inaccuracies are primarily affected by the abundance of core components of the spliceosome assembly and its regulators. Using publicly available data on short-hairpin RNA-knockdowns of numerous spliceosomal components and related regulators, we found support for the importance of RNA-binding proteins in mis-splicing. We also demonstrated that age is positively correlated with mis-splicing, and it affects genes implicated in neurodegenerative diseases. This in-depth characterisation of mis-splicing can have important implications for our understanding of the role of splicing inaccuracies in human disease and the interpretation of long-read RNA-sequencing data.

Modeling splicing outcome by combining 5'ss strength and splicing regulatory elements

Correct pre-mRNA processing in higher eukaryotes vastly depends on splice site recognition. Beyond conserved 5'ss and 3'ss motifs, splicing regulatory elements (SREs) play a pivotal role in this recognition process. Here, we present in silico designed sequences with arbitrary a priori prescribed splicing regulatory HEXplorer properties that can be concatenated to arbitrary length without changing their regulatory properties. We experimentally validated in silico predictions in a massively parallel splicing reporter assay on more than 3000 sequences and exemplarily identified some SRE binding proteins. Aiming at a unified 'functional splice site strength' encompassing both U1 snRNA complementarity and impact from neighboring SREs, we developed a novel RNA-seq based 5'ss usage landscape, mapping the competition of pairs of high confidence 5'ss and neighboring exonic GT sites along HBond and HEXplorer score coordinate axes on human fibroblast and endothelium transcriptome datasets. These RNA-seq data served as basis for a logistic 5'ss usage prediction model, which greatly improved discrimination between strong but unused exonic GT sites and annotated highly used 5'ss. Our 5'ss usage landscape offers a unified view on 5'ss and SRE neighborhood impact on splice site recognition, and may contribute to improved mutation assessment in human genetics.

The plant epitranscriptome: revisiting pseudouridine and 2'-O-methyl RNA modifications

There is growing evidence that post-transcriptional RNA modifications are highly dynamic and can be used to improve crop production. Although more than 172 unique types of RNA modifications have been identified throughout the kingdom of life, we are yet to leverage upon the understanding to optimize RNA modifications in crops to improve productivity. The contributions of internal mRNA modifications such as N6-methyladenosine (m⁶ A) and 5-methylcytosine (m⁵ C) methylations to embryonic development, root development, leaf morphogenesis, flowering, fruit ripening and stress response are sufficiently known, but the roles of the two most abundant RNA modifications, pseudouridine (Ψ) and 2'-O-methylation (Nm), in the cell remain unclear due to insufficient advances in high-throughput technologies in plant development. Therefore, in this review, we discuss the latest methods and insights gained in mapping internal Ψ and Nm and their unique properties in plants and other organisms. In addition, we discuss the limitations that remain in high-throughput technologies for qualitative and quantitative mapping of these RNA modifications and highlight future challenges in regulating the plant epitranscriptome.

Splice site proximity influences alternative exon definition

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

Principles and correction of 5'-splice site selection

In Eukarya, immature mRNA transcripts (pre-mRNA) often contain coding sequences, or exons, interleaved by non-coding sequences, or introns. Introns are removed upon splicing, and further regulation of the retained exons leads to alternatively spliced mRNA. The splicing reaction requires the stepwise assembly of the spliceosome, a macromolecular machine composed of small nuclear ribonucleoproteins (snRNPs). This review focuses on the early stage of spliceosome assembly, when U1 snRNP defines each intron 5’-splice site (5ʹss) in the pre-mRNA. We first introduce the splicing reaction and the impact of alternative splicing on gene expression regulation. Thereafter, we extensively discuss splicing descriptors that influence the 5ʹss selection by U1 snRNP, such as sequence determinants, and interactions mediated by U1-specific proteins or U1 small nuclear RNA (U1 snRNA). We also include examples of diseases that affect the 5ʹss selection by U1 snRNP, and discuss recent therapeutic advances that manipulate U1 snRNP 5ʹss selectivity with antisense oligonucleotides and small-molecule splicing switches.

Altered HIV-1 mRNA Splicing Due to Drug-Resistance-Associated Mutations in Exon 2/2b

The underlying molecular mechanism and their general effect on the replication capacity of HIV 1 drug-resistance-associated mutations is often poorly understood. To elucidate the effect of two such mutations located in a region with a high density of spicing regulatory elements on the HIV-1-splicing outcome, bioinformatic predictions were combined with transfection and infection experiments. Results show that the previously described R263K drug-resistance-associated integrase mutation has additionally a severe effect on the ESE2b splicing regulatory element (SRE) in exon 2b, which causes loss of SD2b recognition. This was confirmed by an R263R silent mutation with a similar predicted effect on the exon 2b SRE. In contrast, a V260I mutation and its silent counterpart with a lower effect on ESS2b did not exhibit any differences in the splicing pattern. Since HIV-1 highly relies on a balanced splicing reaction, changes in the splicing outcome can contribute to changes in viral replication and might add to the effect of escape mutations toward antiviral drugs. Thus, a classification of mutations purely addressing proteins is insufficient.

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.

Functional analyses of human LUC7-like proteins involved in splicing regulation and myeloid neoplasms

Vertebrates have evolved three paralogs, termed LUC7L, LUC7L2, and LUC7L3, of the essential yeast U1 small nuclear RNA (snRNA)-associated splicing factor Luc7p. We investigated the mechanistic and regulatory functions of these putative splicing factors, of which one (LUC7L2) is mutated or deleted in myeloid neoplasms. Protein interaction data show that all three proteins bind similar core but distinct regulatory splicing factors, probably mediated through their divergent arginine-serine-rich domains, which are not present in Luc7p. Knockdown of each factor reveals mostly unique sets of significantly dysregulated alternative splicing events dependent on their binding locations, which are largely non-overlapping. Notably, knockdown of LUC7L2 alone significantly upregulates the expression of multiple spliceosomal factors and downregulates glycolysis genes, possibly contributing to disease pathogenesis. RNA binding studies reveal that LUC7L2 and LUC7L3 crosslink to weak 5' splice sites and to the 5' end of U1 snRNA, establishing an evolutionarily conserved role in 5' splice site selection.

Spliceosomal snRNA Epitranscriptomics

Small nuclear RNAs (snRNAs) are critical components of the spliceosome that catalyze the splicing of pre-mRNA. snRNAs are each complexed with many proteins to form RNA-protein complexes, termed as small nuclear ribonucleoproteins (snRNPs), in the cell nucleus. snRNPs participate in pre-mRNA splicing by recognizing the critical sequence elements present in the introns, thereby forming active spliceosomes. The recognition is achieved primarily by base-pairing interactions (or nucleotide-nucleotide contact) between snRNAs and pre-mRNA. Notably, snRNAs are extensively modified with different RNA modifications, which confer unique properties to the RNAs. Here, we review the current knowledge of the mechanisms and functions of snRNA modifications and their biological relevance in the splicing process.

Mechanisms behind polyphagia in a pest insect: Responses of Spodoptera frugiperda (J.E. Smith) strains to preferential and alternative larval host plants assessed with gene regulatory networks

A dataset of gene expression from Spodoptera frugiperda, a highly generalist pest moth, was used to understand how gene regulation is related to larval host plant preference. Transcriptomic data of corn and rice strains of S. frugiperda larvae, reared on different diets, were analysed with three different approaches of gene network inference, namely co-expression, weighted co-expression and Bayesian networks, since each methodology provides a different visualization of the data. Using these approaches, it was possible to identify two loosely interconnected co-expression networks, one of them responsible for fast response to herbivory and anti-herbivory mechanisms and the other related to housekeeping genes, which present slower response to environmental variations. Integrating different levels of information such as gene expression patterns, gene assembly, transcriptomics, relationship among genes and phenotypes, functional relationships, among other information, enabled a wider visualization of S. frugiperda response to diet stimuli. The biological properties in the proposed networks are here described and discussed, as well as patterns of gene expression related to larval performance attributes.

VarCon: An R Package for Retrieving Neighboring Nucleotides of an SNV

Reporting of a single nucleotide variant (SNV) follows the Sequence Variant Nomenclature (http://varnomen.hgvs.org/), using an unambiguous numbering scheme specific for coding and noncoding DNA. However, the corresponding sequence neighborhood of a given SNV, which is required to assess its impact on splicing regulation, is not easily accessible from this nomenclature. Providing fast and easy access to this neighborhood just from a given SNV reference, the novel tool VarCon combines information of the Ensembl human reference genome and the corresponding transcript table for accurate retrieval. VarCon also displays splice site scores (HBond and MaxEnt scores) and HEXplorer profiles of an SNV neighborhood, reflecting position-dependent splice enhancing and silencing properties.

Regulation and Function of RNA Pseudouridylation in Human Cells

Recent advances in pseudouridine detection reveal a complex pseudouridine landscape that includes messenger RNA and diverse classes of noncoding RNA in human cells. The known molecular functions of pseudouridine, which include stabilizing RNA conformations and destabilizing interactions with varied RNA-binding proteins, suggest that RNA pseudouridylation could have widespread effects on RNA metabolism and gene expression. Here, we emphasize how much remains to be learned about the RNA targets of human pseudouridine synthases, their basis for recognizing distinct RNA sequences, and the mechanisms responsible for regulated RNA pseudouridylation. We also examine the roles of noncoding RNA pseudouridylation in splicing and translation and point out the potential effects of mRNA pseudouridylation on protein production, including in the context of therapeutic mRNAs.

A scale-free analysis of the HIV-1 genome demonstrates multiple conserved regions of structural and functional importance

HIV-1 replicates via a low-fidelity polymerase with a high mutation rate; strong conservation of individual nucleotides is highly indicative of the presence of critical structural or functional properties. Identifying such conservation can reveal novel insights into viral behaviour. We analysed 3651 publicly available sequences for the presence of nucleic acid conservation beyond that required by amino acid constraints, using a novel scale-free method that identifies regions of outlying score together with a codon scoring algorithm. Sequences with outlying score were further analysed using an algorithm for producing local RNA folds whilst accounting for alignment properties. 11 different conserved regions were identified, some corresponding to well-known cis-acting functions of the HIV-1 genome but also others whose conservation has not previously been noted. We identify rational causes for many of these, including cis functions, possible additional reading frame usage, a plausible mechanism by which the central polypurine tract primes second-strand DNA synthesis and a conformational stabilising function of a region at the 5′ end of env.

HIV-1 is a very rapidly mutating organism, however some parts of its genetic material change more than others. We looked for coding regions of HIV-1 that change relatively little, by turning the problem of finding such regions into a problem in signal processing, and solving this using a novel analytical approach that we recently described. We investigated why the regions we identified change less, including using the genetic code in the regions we found to prime an algorithm to predict their structures. In some cases there are already known functions for the features we found, in others they provide new insights into the properties of known regions, and in some cases we identify new regions that vary less for as yet unknown functional reasons.

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.

Modifications in small nuclear RNAs and their roles in spliceosome assembly and function

Modifications in cellular RNAs have emerged as key regulators of all aspects of gene expression, including pre-mRNA splicing. During spliceosome assembly and function, the small nuclear RNAs (snRNAs) form numerous dynamic RNA-RNA and RNA-protein interactions, which are required for spliceosome assembly, correct positioning of the spliceosome on substrate pre-mRNAs and catalysis. The human snRNAs contain several base methylations as well as a myriad of pseudouridines and 2'-O-methylated nucleotides, which are largely introduced by small Cajal body-specific ribonucleoproteins (scaRNPs). Modified nucleotides typically cluster in functionally important regions of the snRNAs, suggesting that their presence could optimise the interactions of snRNAs with each other or with pre-mRNAs, or may affect the binding of spliceosomal proteins. snRNA modifications appear to play important roles in snRNP biogenesis and spliceosome assembly, and have also been proposed to influence the efficiency and fidelity of pre-mRNA splicing. Interestingly, alterations in the modification status of snRNAs have recently been observed in different cellular conditions, implying that some snRNA modifications are dynamic and raising the possibility that these modifications may fine-tune the spliceosome for particular functions. Here, we review the current knowledge on the snRNA modification machinery and discuss the timing, functions and dynamics of modifications in snRNAs.

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.

Exon Junction Complex Shapes the Transcriptome by Repressing Recursive Splicing

Recursive splicing (RS) starts by defining an "RS-exon," which is then spliced to the preceding exon, thus creating a recursive 5' splice site (RS-5ss). Previous studies focused on cryptic RS-exons, and now we find that the exon junction complex (EJC) represses RS of hundreds of annotated, mainly constitutive RS-exons. The core EJC factors, and the peripheral factors PNN and RNPS1, maintain RS-exon inclusion by repressing spliceosomal assembly on RS-5ss. The EJC also blocks 5ss located near exon-exon junctions, thus repressing inclusion of cryptic microexons. The prevalence of annotated RS-exons is high in deuterostomes, while the cryptic RS-exons are more prevalent in Drosophila, where EJC appears less capable of repressing RS. Notably, incomplete repression of RS also contributes to physiological alternative splicing of several human RS-exons. Finally, haploinsufficiency of the EJC factor Magoh in mice is associated with skipping of RS-exons in the brain, with relevance to the microcephaly phenotype and human diseases.

Old gene, new phenotype: splice-altering variants in CEACAM16 cause recessive non-syndromic hearing impairment

BACKGROUND

Hearing loss is a genetically and phenotypically heterogeneous disorder.

OBJECTIVES

The purpose of this study was to determine the genetic cause underlying the postlingual progressive hearing loss in two Iranian families.

METHODS

We used OtoSCOPE, a next-generation sequencing platform targeting >150 genes causally linked to deafness, to screen two deaf probands. Data analysis was completed using a custom bioinformatics pipeline, and variants were functionally assessed using minigene splicing assays.

RESULTS

We identified two homozygous splice-altering variants (c.37G>T and c.662-1G>C) in the CEACAM16 gene, segregating with the deafness in each family. The minigene splicing results revealed the c.37G>T results in complete skipping of exon 2 and loss of the AUG start site. The c.662-1G>C activates a cryptic splice site inside exon 5 resulting in a shift in the mRNA reading frame.

CONCLUSIONS

These results suggest that loss-of-function mutations in CEACAM16 result in postlingual progressive hearing impairment and further support the role of CEACAM16 in auditory function.

A comprehensive characterization of cis-acting splicing-associated variants in human cancer

Although many driver mutations are thought to promote carcinogenesis via abnormal splicing, the landscape of splicing-associated variants (SAVs) remains unknown due to the complexity of splicing abnormalities. Here, we developed a statistical framework to systematically identify SAVs disrupting or newly creating splice site motifs and applied it to matched whole-exome and transcriptome sequencing data from 8976 samples across 31 cancer types, generating a catalog of 14,438 SAVs. Such a large collection of SAVs enabled us to characterize their genomic features, underlying mutational processes, and influence on cancer driver genes. In fact, ∼50% of SAVs identified were those disrupting noncanonical splice sites (non-GT-AG dinucleotides), including the third and fifth intronic bases of donor sites, or newly creating splice sites. Mutation signature analysis revealed that tobacco smoking is more strongly associated with SAVs, whereas ultraviolet exposure has less impact. SAVs showed remarkable enrichment of cancer-related genes, and as many as 14.7% of samples harbored at least one SAVs affecting them, particularly in tumor suppressors. In addition to intron retention, whose association with tumor suppressor inactivation has been previously reported, exon skipping and alternative splice site usage caused by SAVs frequently affected tumor suppressors. Finally, we described high-resolution distributions of SAVs along the gene and their splicing outcomes in commonly disrupted genes, including TP53, PIK3R1, GATA3, and CDKN2A, which offers genetic clues for understanding their functional properties. Collectively, our findings delineate a comprehensive portrait of SAVs, novel insights into transcriptional de-regulation in cancer.

Functional Interplay between Small Non-Coding RNAs and RNA Modification in the Brain

Small non-coding RNAs are essential for transcription, translation and gene regulation in all cell types, but are particularly important in neurons, with known roles in neurodevelopment, neuroplasticity and neurological disease. Many small non-coding RNAs are directly involved in the post-transcriptional modification of other RNA species, while others are themselves substrates for modification, or are functionally modulated by modification of their target RNAs. In this review, we explore the known and potential functions of several distinct classes of small non-coding RNAs in the mammalian brain, focusing on the newly recognised interplay between the epitranscriptome and the activity of small RNAs. We discuss the potential for this relationship to influence the spatial and temporal dynamics of gene activation in the brain, and predict that further research in the field of epitranscriptomics will identify interactions between small RNAs and RNA modifications which are essential for higher order brain functions such as learning and memory.

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.

The Role of Noncoding RNA Pseudouridylation in Nuclear Gene Expression Events

Pseudouridine is the most abundant internal RNA modification in stable noncoding RNAs (ncRNAs). It can be catalyzed by both RNA-dependent and RNA-independent mechanisms. Pseudouridylation impacts both the biochemical and biophysical properties of RNAs and thus influences RNA-mediated cellular processes. The investigation of nuclear-ncRNA pseudouridylation has demonstrated that it is critical for the proper control of multiple stages of gene expression regulation. Here, we review how nuclear-ncRNA pseudouridylation contributes to transcriptional regulation and pre-mRNA splicing.

Succession of splicing regulatory elements determines cryptic 5΄ss functionality

A critical step in exon definition is the recognition of a proper splice donor (5΄ss) by the 5’ end of U1 snRNA. In the selection of appropriate 5΄ss, cis-acting splicing regulatory elements (SREs) are indispensable. As a model for 5΄ss recognition, we investigated cryptic 5΄ss selection within the human fibrinogen Bβ-chain gene (FGB) exon 7, where we identified several exonic SREs that simultaneously acted on up- and downstream cryptic 5΄ss. In the FGB exon 7 model system, 5΄ss selection iteratively proceeded along an alternating sequence of U1 snRNA binding sites and interleaved SREs which in principle supported different 3’ exon ends. Like in a relay race, SREs either suppressed a potential 5΄ss and passed the splicing baton on or splicing actually occurred. From RNA-Seq data, we systematically selected 19 genes containing exons with silent U1 snRNA binding sites competing with nearby highly used 5΄ss. Extensive SRE analysis by different algorithms found authentic 5΄ss significantly more supported by SREs than silent U1 snRNA binding sites, indicating that our concept may permit generalization to a model for 5΄ss selection and 3’ exon end definition.

RNA splicing in human disease and in the clinic

Defects at the level of the pre-mRNA splicing process represent a major cause of human disease. Approximately 15-50% of all human disease mutations have been shown to alter functioning of basic and auxiliary splicing elements. These elements are required to ensure proper processing of pre-mRNA splicing molecules, with their disruption leading to misprocessing of the pre-mRNA molecule and disease. The splicing process is a complex process, with much still to be uncovered before we are able to accurately predict whether a reported genomic sequence variant (GV) represents a splicing-associated disease mutation or a harmless polymorphism. Furthermore, even when a mutation is correctly identified as affecting the splicing process, there still remains the difficulty of providing an exact evaluation of the potential impact on disease onset, severity and duration. In this review, we provide a brief overview of splicing diagnostic methodologies, from in silico bioinformatics approaches to wet lab in vitro and in vivo systems to evaluate splicing efficiencies. In particular, we provide an overview of how the latest developments in high-throughput sequencing can be applied to the clinic, and are already changing clinical approaches.

Analysis of Competing HIV-1 Splice Donor Sites Uncovers a Tight Cluster of Splicing Regulatory Elements within Exon 2/2b

The HIV-1 accessory protein Vif is essential for viral replication by counteracting the host restriction factor APOBEC3G (A3G), and balanced levels of both proteins are required for efficient viral replication. Noncoding exons 2/2b contain the Vif start codon between their alternatively used splice donors 2 and 2b (D2 and D2b). For vif mRNA, intron 1 must be removed while intron 2 must be retained. Thus, splice acceptor 1 (A1) must be activated by U1 snRNP binding to either D2 or D2b, while splicing at D2 or D2b must be prevented. Here, we unravel the complex interactions between previously known and novel components of the splicing regulatory network regulating HIV-1 exon 2/2b inclusion in viral mRNAs. In particular, using RNA pulldown experiments and mass spectrometry analysis, we found members of the heterogeneous nuclear ribonucleoparticle (hnRNP) A/B family binding to a novel splicing regulatory element (SRE), the exonic splicing silencer ESS2b, and the splicing regulatory proteins Tra2/SRSF10 binding to the nearby exonic splicing enhancer ESE2b. Using a minigene reporter, we performed bioinformatics HEXplorer-guided mutational analysis to narrow down SRE motifs affecting splice site selection between D2 and D2b. Eventually, the impacts of these SREs on the viral splicing pattern and protein expression were exhaustively analyzed in viral particle production and replication experiments. Masking of these protein binding sites by use of locked nucleic acids (LNAs) impaired Vif expression and viral replication.IMPORTANCE Based on our results, we propose a model in which a dense network of SREs regulates vif mRNA and protein expression, crucial to maintain viral replication within host cells with varying A3G levels and at different stages of infection. This regulation is maintained by several serine/arginine-rich splicing factors (SRSF) and hnRNPs binding to those elements. Targeting this cluster of SREs with LNAs may lead to the development of novel effective therapeutic strategies.

A Novel Alzheimer-Associated SNP in Tmp21 Increases Amyloidogenesis

Recent studies suggest that TMP21 is a selective modulator of γ-secretase and its dysregulation affects APP processing, leading to increased Aβ generation. However, the genetic association between Tmp21 and Alzheimer's disease (AD) remains elusive. In this study, we identified that a novel single-nucleotide polymorphism (SNP) rs12435391 (IVS4-28T>C) in intron 4 of Tmp21 was genetically associated with AD. We found that allele C of the SNP rs12435391 did not affect splicing site recognition, but it significantly increased TMP21 gene expression. The stability of Tmp21 pre-mRNA and the transcription of Tmp21 were not affected by allele C of the SNP rs12435391. However, allele C of the SNP rs12435391 significantly increased the splicing efficiency of Tmp21 pre-mRNA, leading to the elevation of mature mRNA. Furthermore, allele C of the SNP rs12435391 significantly reduced C83 level and increased Aβ generation. Taken together, our study suggests that TMP21 is genetically associated with Alzheimer's disease, with the novel Tmp21 SNP as a risk factor for Alzheimer's pathogenesis.

Differential hnRNP D isoform incorporation may confer plasticity to the ESSV-mediated repressive state across HIV-1 exon 3

Even though splicing repression by hnRNP complexes bound to exonic sequences is well-documented, the responsible effector domains of hnRNP proteins have been described for only a select number of hnRNP constituents. Thus, there is only limited information available for possible varying silencer activities amongst different hnRNP proteins and composition changes within possible hnRNP complex assemblies. In this study, we identified the glycine-rich domain (GRD) of hnRNP proteins as a unifying feature in splice site repression. We also show that all four hnRNP D isoforms can act as genuine splicing repressors when bound to exonic positions. The presence of an extended GRD, however, seemed to potentiate the hnRNP D silencer activity of isoforms p42 and p45. Moreover, we demonstrate that hnRNP D proteins associate with the HIV-1 ESSV silencer complex, probably through direct recognition of "UUAG" sequences overlapping with the previously described "UAGG" motifs bound by hnRNP A1. Consequently, this spatial proximity seems to cause mutual interference between hnRNP A1 and hnRNP D. This interplay between hnRNP A1 and D facilitates a dynamic regulation of the repressive state of HIV-1 exon 3 which manifests as fluctuating relative levels of spliced vpr- and unspliced gag/pol-mRNAs.

Human immunodeficiency virus type 1 splicing at the major splice donor site is controlled by highly conserved RNA sequence and structural elements

Human immunodeficiency virus type 1 (HIV-1) splicing has to be strictly controlled to ensure the balanced production of the unspliced and all differently spliced viral RNAs. Splicing at the major 59 splice site (59ss) that is used for the synthesis of all spliced RNAs is modulated by the local RNA structure and binding of regulatory SR proteins. Here, we demonstrate that the suboptimal sequence complementarity between this 59ss and U1 small nuclear RNA (snRNA) also contributes to prevent excessive splicing. Analysis of a large set of HIV-1 sequences revealed that all three regulatory features of the 59ss region (RNA structure, SR protein binding and sequence complementarity with U1 snRNA) are highly conserved amongst virus isolates, which supports their importance. Combined mutations that destabilize the local RNA structure, remove binding sites for inhibitory SR proteins and optimize the U1 snRNA complementarity resulted in almost complete splicing and accordingly reduced virus replication.

Exon-skipping and mRNA decay in human liver tissue: molecular consequences of pathogenic bile salt export pump mutations

The bile salt export pump BSEP mediates bile formation. Over 150 BSEP mutations are associated with progressive familial intrahepatic cholestasis type 2 (PFIC-2), with few characterised specifically. We examined liver tissues from two PFIC-2 patients compound heterozygous for the splice-site mutation c.150 + 3A > C and either c.2783_2787dup5 resulting in a frameshift with a premature termination codon (child 1) or p.R832C (child 2). Splicing was analysed with a minigene system and mRNA sequencing from patients’ livers. Protein expression was shown by immunofluorescence. Using the minigene, c.150 + 3A > C causes complete skipping of exon 3. In liver tissue of child 1, c.2783_2787dup5 was found on DNA but not on mRNA level, implying nonsense-mediated mRNA decay (NMD) when c.2783_2787dup5 is present. Still, BSEP protein as well as mRNA with and without exon 3 were detectable and can be assigned to the c.150 + 3A > C allele. Correctly spliced transcripts despite c.150 + 3A > C were also confirmed in liver of child 2. In conclusion, we provide evidence (1) for effective NMD due to a BSEP frameshift mutation and (2) partial exon-skipping due to c.150 + 3A > C. The results illustrate that the extent of exon-skipping depends on the genomic and cellular context and that regulation of splicing may have therapeutic potential.

Functional Dissection of an Alternatively Spliced Herpesvirus Gene by Splice Site Mutagenesis

Herpesviruses have large and complex DNA genomes. The largest among the herpesviruses, those of the cytomegaloviruses, include over 170 genes. Although most herpesvirus gene products are expressed from unspliced transcripts, a substantial number of viral transcripts are spliced. Some viral transcripts are subject to alternative splicing, which leads to the expression of several proteins from a single gene. Functional analysis of individual proteins derived from an alternatively spliced gene is difficult, as deletion and nonsense mutagenesis, both common methods used in the generation of viral gene knockout mutants, affect several or all gene products at the same time. Here, we show that individual gene products of an alternatively spliced herpesvirus gene can be inactivated selectively by mutagenesis of the splice donor or acceptor site and by intron deletion or substitution mutagenesis. We used this strategy to dissect the essential M112/113 gene of murine cytomegalovirus (MCMV), which encodes the MCMV Early 1 (E1) proteins. The expression of each of the four E1 protein isoforms was inactivated individually, and the requirement for each isoform in MCMV replication was analyzed in fibroblasts, endothelial cells, and macrophages. We show that the E1 p87 isoform, but not the p33, p36, and p38 isoforms, is essential for viral replication in cell culture. Moreover, the presence of one of the two medium-size isoforms (p36 or p38) and the presence of intron 1, but not its specific sequence, are required for viral replication. This study demonstrates the usefulness of splice site mutagenesis for the functional analysis of alternatively spliced herpesvirus genes.

IMPORTANCE

Herpesviruses include up to 170 genes in their DNA genomes. The functions of most viral gene products remain poorly defined. The construction of viral gene knockout mutants has thus been an important tool for functional analysis of viral proteins. However, this strategy is of limited use when viral gene transcripts are alternatively spliced, leading to the expression of several proteins from a single gene. In this study, we showed, as a proof of principle, that each protein product of an alternatively spliced gene can be eliminated individually by splice site mutagenesis. Mutant viruses lacking individual protein products displayed different phenotypes, demonstrating that the products of alternatively spliced genes have nonredundant functions.

Human endogenous retrovirus HERV-K(HML-2) activity in prostate cancer is dominated by a few loci

BACKGROUND

Increased expression of human endogenous retroviruses, especially HERV-K(HML-2) proviruses, has recently been associated with prostate carcinoma progression. In particular, a HML-2 locus in chromosome 22q11.23 (H22q) is upregulated in many cases. We therefore aimed at delineating the extent and repertoire of HML-2 transcription in prostate cancer tissues and cell lines and to define the transcription pattern and biological effects of H22q.

METHODS

Sanger and high throughput amplicon sequencing was used to define the repertoire of expressed HML-2 in a selected set of samples. qRT-PCR was used to quantify expression of selected proviruses in an extended set of prostate cancer tissues. Transcription factor binding sites (TFBS) were compared bioinformatically using the Transfac database. Expression of H22q was further characterized by siRNA-mediated knockdown, 5' RACE mapping of transcriptional start sites (TSS) and identification of splice sites. Functional effects of H22q knockdown were investigated by viability and apoptosis assays.

RESULTS

In addition to H22q, a limited number of other proviruses were found expressed by sequencing. Of these, provirus ERVK-5 and to a lesser degree ERVK-15 were frequently upregulated in prostate cancer. In contrast, expression of ERVK-24, predominant in germ cell tumors, was not detectable in prostatic tissues. While HML-2 LTRs contain binding sites for the androgen receptor and cofactors, no consistent differences in transcription factor binding sites were found between expressed and non-expressed proviruses. The H22q locus contains two 5'-LTRs of which the upstream LTR is predominantly used in prostatic cells, with an imprecise TSS. Splicing of H22q transcripts is complex, generating, among others, a transcript with an Np9-like ORF. Knockdown of H22q did not significantly affect proliferation or apoptosis of prostate cancer cells.

CONCLUSIONS

Our findings further underline that HML-2 expression is commonly highly tissue-specific. In prostate cancer, a limited number of loci become activated, especially H22q and ERVK-5. As expressed and non-expressed proviruses do not differ significantly in TFBS, tissue- and tumor-specific expression may be governed primarily by chromatin context. Overexpression of HML-2 H22q is more likely consequence than cause of prostate cancer progression.

Intron retention is a widespread mechanism of tumor-suppressor inactivation

A substantial fraction of disease-causing mutations are pathogenic through aberrant splicing. Although genome profiling studies have identified somatic single-nucleotide variants (SNVs) in cancer, the extent to which these variants trigger abnormal splicing has not been systematically examined. Here we analyzed RNA sequencing and exome data from 1,812 patients with cancer and identified ∼900 somatic exonic SNVs that disrupt splicing. At least 163 SNVs, including 31 synonymous ones, were shown to cause intron retention or exon skipping in an allele-specific manner, with ∼70% of the SNVs occurring on the last base of exons. Notably, SNVs causing intron retention were enriched in tumor suppressors, and 97% of these SNVs generated a premature termination codon, leading to loss of function through nonsense-mediated decay or truncated protein. We also characterized the genomic features predictive of such splicing defects. Overall, this work demonstrates that intron retention is a common mechanism of tumor-suppressor inactivation.

HIV-1 splicing is controlled by local RNA structure and binding of splicing regulatory proteins at the major 5' splice site

The 5' leader region of the human immunodeficiency virus 1 (HIV-1) RNA genome contains the major 5' splice site (ss) that is used in the production of the many spliced viral RNAs. This splice-donor (SD) region can fold into a stable stem-loop structure and the thermodynamic stability of this RNA hairpin influences splicing efficiency. In addition, splicing may be modulated by binding of splicing regulatory (SR) proteins, in particular SF2/ASF (SRSF1), SC35 (SRSF2), SRp40 (SRSF5) and SRp55 (SRSF6), to sequence elements in the SD region. The role of RNA structure and SR protein binding in splicing control was previously studied by functional analysis of mutant SD sequences. The interpretation of these studies was complicated by the fact that most mutations simultaneously affect both structure and sequence elements. We therefore tried to disentangle the contribution of these two variables by designing more precise SD region mutants with a single effect on either the sequence or the structure. The current analysis indicates that HIV-1 splicing at the major 5'ss is modulated by both the stability of the local RNA structure and the binding of splicing regulatory proteins.

rbamtools: an R interface to samtools enabling fast accumulative tabulation of splicing events over multiple RNA-seq samples

The open source environment R isf the most widely used software to statistically explore biological data sets including sequence alignments. BAM is the de facto standard file format for sequence alignment. With rbamtools, we provide now a full spectrum of accessibility to BAM for R users such as reading, writing, extraction of subsets and plotting of alignment depth where the script syntax closely follows the SAM/BAM format. Additionally, rbamtools enables fast accumulative tabulation of splicing events over multiple BAM files.

HIV-1 splicing at the major splice donor site is restricted by RNA structure

The 5' leader region of the HIV-1 RNA contains the major 5' splice site (ss) that is used in the production of all spliced viral RNAs. This splice-donor (SD) region can fold a stem-loop structure. We demonstrate that whereas stabilization of this SD hairpin reduces splicing efficiency, destabilization increases splicing. Both stabilization and destabilization reduce viral fitness. These results demonstrate that the stability of the SD hairpin can modulate the level of splicing, most likely by controlling the accessibility of the 5'ss for the splicing machinery. The natural stability of the SD hairpin restricts splicing and this stability seems to be fine-tuned to reach the optimal balance between unspliced and spliced RNAs for efficient virus replication. The 5'ss region of different HIV-1 isolates and the related SIVmac239 can fold a similar structure. This evolutionary conservation supports the importance of this structure in viral replication.

Genomic HEXploring allows landscaping of novel potential splicing regulatory elements

Effective splice site selection is critically controlled by flanking splicing regulatory elements (SREs) that can enhance or repress splice site use. Although several computational algorithms currently identify a multitude of potential SRE motifs, their predictive power with respect to mutation effects is limited. Following a RESCUE-type approach, we defined a hexamer-based ‘HEXplorer score’ as average Z-score of all six hexamers overlapping with a given nucleotide in an arbitrary genomic sequence. Plotted along genomic regions, HEXplorer score profiles varied slowly in the vicinity of splice sites. They reflected the respective splice enhancing and silencing properties of splice site neighborhoods beyond the identification of single dedicated SRE motifs. In particular, HEXplorer score differences between mutant and reference sequences faithfully represented exonic mutation effects on splice site usage. Using the HIV-1 pre-mRNA as a model system highly dependent on SREs, we found an excellent correlation in 29 mutations between splicing activity and HEXplorer score. We successfully predicted and confirmed five novel SREs and optimized mutations inactivating a known silencer. The HEXplorer score allowed landscaping of splicing regulatory regions, provided a quantitative measure of mutation effects on splice enhancing and silencing properties and permitted calculation of the mutationally most effective nucleotide.

Molecular events accompanying rous sarcoma virus rescue from rodent cells and the role of viral gene complementation

Transformation of rodent cells with avian Rous sarcoma virus (RSV) opened new ways to studying virus integration and expression in nonpermissive cells. We were interested in (i) the molecular changes accompanying fusion of RSV-transformed mammalian cells with avian cells leading to virus rescue and (ii) enhancement of this process by retroviral gene products. The RSV-transformed hamster RSCh cell line was characterized as producing only a marginal amount of env mRNA, no envelope glycoprotein, and a small amount of unprocessed Gag protein. Egress of viral unspliced genomic RNA from the nucleus was hampered, and its stability decreased. Cell fusion of the chicken DF-1 cell line with RSCh cells led to production of env mRNA, envelope glycoprotein, and processed Gag and virus-like particle formation. Proteosynthesis inhibition in DF-1 cells suppressed steps leading to virus rescue. Furthermore, new aberrantly spliced env mRNA species were found in the RSCh cells. Finally, we demonstrated that virus rescue efficiency can be significantly increased by complementation with the env gene and the highly expressed gag gene and can be increased the most by a helper virus infection. In summary, Env and Gag synthesis is increased after RSV-transformed hamster cell fusion with chicken fibroblasts, and both proteins provided in trans enhance RSV rescue. We conclude that the chicken fibroblast yields some factor(s) needed for RSV replication, particularly Env and Gag synthesis, in nonpermissive rodent cells.

IMPORTANCE

One of the important issues in retrovirus heterotransmission is related to cellular factors that prevent virus replication. Rous sarcoma virus (RSV), a member of the avian sarcoma and leukosis family of retroviruses, is able to infect and transform mammalian cells; however, such transformed cells do not produce infectious virus particles. Using the well-defined model of RSV-transformed rodent cells, we established that the lack of virus replication is due to the absence of chicken factor(s), which can be supplemented by cell fusion. Cell fusion with permissive chicken cells led to an increase in RNA splicing and nuclear export of specific viral mRNAs, as well as synthesis of respective viral proteins and production of virus-like particles. RSV rescue by cell fusion can be potentiated by in trans expression of viral genes in chicken cells. We conclude that rodent cells lack some chicken factor(s) required for proper viral RNA processing and viral protein synthesis.

In silico tools for splicing defect prediction: a survey from the viewpoint of end users

RNA splicing is the process during which introns are excised and exons are spliced. The precise recognition of splicing signals is critical to this process and mutations affecting splicing comprise a considerable proportion of genetic disease etiology. Analysis of RNA samples from the patient is the most straightforward and reliable method to detect splicing defects. However, currently the technical limitation prohibits its use in routine clinical practice. In silico tools that predict potential consequences of splicing mutations may be useful in daily diagnostic activities. In this review, we provide medical geneticists with some basic insights into some of the most popular in silico tools for splicing defect prediction, from the viewpoint of end-users. Bioinformaticians in relevant areas who are working on huge datasets may also benefit from this review. Specifically, we focus on those tools whose primary goal is to predict the impact of mutations within the 5′ and 3′ splicing consensus regions: the algorithms used by different tools as well as their major advantages and disadvantages are briefly introduced; the formats of their input and output are summarized; and the interpretation, evaluation, and prospection are also discussed.

Regulation of constitutive and alternative splicing by PRMT5 reveals a role for Mdm4 pre-mRNA in sensing defects in the spliceosomal machinery

The tight control of gene expression at the level of both transcription and post-transcriptional RNA processing is essential for mammalian development. We here investigate the role of protein arginine methyltransferase 5 (PRMT5), a putative splicing regulator and transcriptional cofactor, in mammalian development. We demonstrate that selective deletion of PRMT5 in neural stem/progenitor cells (NPCs) leads to postnatal death in mice. At the molecular level, the absence of PRMT5 results in reduced methylation of Sm proteins, aberrant constitutive splicing, and the alternative splicing of specific mRNAs with weak 5' donor sites. Intriguingly, the products of these mRNAs are, among others, several proteins regulating cell cycle progression. We identify Mdm4 as one of these key mRNAs that senses the defects in the spliceosomal machinery and transduces the signal to activate the p53 response, providing a mechanistic explanation of the phenotype observed in vivo. Our data demonstrate that PRMT5 is a master regulator of splicing in mammals and uncover a new role for the Mdm4 pre-mRNA, which could be exploited for anti-cancer therapy.

Pick one, but be quick: 5' splice sites and the problems of too many choices

Splice site selection is fundamental to pre-mRNA splicing and the expansion of genomic coding potential. 5' Splice sites (5'ss) are the critical elements at the 5' end of introns and are extremely diverse, as thousands of different sequences act as bona fide 5'ss in the human transcriptome. Most 5'ss are recognized by base-pairing with the 5' end of the U1 small nuclear RNA (snRNA). Here we review the history of research on 5'ss selection, highlighting the difficulties of establishing how base-pairing strength determines splicing outcomes. We also discuss recent work demonstrating that U1 snRNA:5'ss helices can accommodate noncanonical registers such as bulged duplexes. In addition, we describe the mechanisms by which other snRNAs, regulatory proteins, splicing enhancers, and the relative positions of alternative 5'ss contribute to selection. Moreover, we discuss mechanisms by which the recognition of numerous candidate 5'ss might lead to selection of a single 5'ss and propose that protein complexes propagate along the exon, thereby changing its physical behavior so as to affect 5'ss selection.

An intronic G run within HIV-1 intron 2 is critical for splicing regulation of vif mRNA

Within target T lymphocytes, human immunodeficiency virus type I (HIV-1) encounters the retroviral restriction factor APOBEC3G (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3G; A3G), which is counteracted by the HIV-1 accessory protein Vif. Vif is encoded by intron-containing viral RNAs that are generated by splicing at 3' splice site (3'ss) A1 but lack splicing at 5'ss D2, which results in the retention of a large downstream intron. Hence, the extents of activation of 3'ss A1 and repression of D2, respectively, determine the levels of vif mRNA and thus the ability to evade A3G-mediated antiviral effects. The use of 3'ss A1 can be enhanced or repressed by splicing regulatory elements that control the recognition of downstream 5'ss D2. Here we show that an intronic G run (G(I2)-1) represses the use of a second 5'ss, termed D2b, that is embedded within intron 2 and, as determined by RNA deep-sequencing analysis, is normally inefficiently used. Mutations of G(I2)-1 and activation of D2b led to the generation of transcripts coding for Gp41 and Rev protein isoforms but primarily led to considerable upregulation of vif mRNA expression. We further demonstrate, however, that higher levels of Vif protein are actually detrimental to viral replication in A3G-expressing T cell lines but not in A3G-deficient cells. These observations suggest that an appropriate ratio of Vif-to-A3G protein levels is required for optimal virus replication and that part of Vif level regulation is effected by the novel G run identified here.

The HIV-1 major splice donor D1 is activated by splicing enhancer elements within the leader region and the p17-inhibitory sequence

Usage of the HIV-1 major 5' splice site D1 is a prerequisite for generation of all spliced viral mRNAs encoding essential regulatory and structural proteins. We set out to determine whether flanking sequences ensure D1-activation. We found that an exonic splicing enhancer function is exerted by the region upstream of D1, which is crucially required for its activation. Additionally, we identified an intronic splicing regulatory element within the p17-instability element of the Gag-ORF enhancing D1-activation. Furthermore, our experimental data demonstrated that sequence motifs displaying high similarity to consensus binding sites for SR protein SC35 (SRSF2) overlapping with D1 fine-tune its activation. Our results reveal that D1-activation is safe-guarded by the interplay of upstream and downstream located splicing enhancer elements ensuring usage of D1 even if its strength is decreased upon mutation. The identification of sequence elements activating D1-usage sheds further light on the balanced expression of alternatively spliced HIV-1 mRNAs.

DNA diagnostics and exon skipping

The predominate form of DNA diagnostics remains nucleic acid sequencing in the research and clinical setting. While DNA sequencing allows a mutation to be correctly identified, only RNA sequencing can confirm the effect of that mutation on the resulting mRNA transcript. In the absence of RNA sequencing, predictions are reliant on either experimental studies or bioinformatic modelling. While each of these approaches provides insights into cellular splicing choices, of which exon skipping is but one, both possess inherent weaknesses. A method which is able to integrate and appropriately weigh the various factors influencing cellular splicing choices into an accurate, comprehensive modelling tool still remains elusive.In this overview chapter, the current methods utilised for DNA diagnostics and the impact of the emerging next-generation sequencing techniques are considered. We explore why RNA remains a problematic medium with which to work. To understand how exon skipping can be predicted from a DNA sequence, the key cis-acting elements influencing splicing are reviewed. Finally, the current methods used to predict exon skipping including RNA-based studies, experimental studies, and bioinformatic modelling approaches are outlined.