The splicing code

Alternative splicing of the Snap23 microexon is regulated by MBNL, QKI, and RBFOX2 in a tissue-specific manner and is altered in striated muscle diseases

The reprogramming of alternative splicing networks during development is a hallmark of tissue maturation and identity. Alternative splicing of microexons (small, genomic regions ≤ 51 nucleotides) functionally regulate protein-protein interactions in the brain and is altered in several neuronal diseases. However, little is known about the regulation and function of alternatively spliced microexons in striated muscle. Here, we investigated alternative splicing of a microexon in the synaptosome-associated protein 23 (Snap23) encoded gene. We found that inclusion of this microexon is developmentally regulated and tissue-specific, as it occurs exclusively in adult heart and skeletal muscle. The alternative region is highly conserved in mammalian species and encodes an in-frame sequence of 11 amino acids. Furthermore, we showed that alternative splicing of this microexon is mis-regulated in mouse models of heart and skeletal muscle diseases. We identified the RNA-binding proteins (RBPs) quaking (QKI) and RNA binding fox-1 homolog 2 (RBFOX2) as the primary splicing regulators of the Snap23 microexon. We found that QKI and RBFOX2 bind downstream of the Snap23 microexon to promote its inclusion, and this regulation can be escaped when the weak splice donor is mutated to the consensus 5' splice site. Finally, we uncovered the interplay between QKI and muscleblind-like splicing regulator (MBNL) as an additional, but minor layer of Snap23 microexon splicing control. Our results are one of the few reports detailing microexon alternative splicing regulation during mammalian striated muscle development.

Data-driven insights to inform splice-altering variant assessment

Disease-causing genetic variants often disrupt mRNA splicing, an intricate process that is incompletely understood. Thus, accurate inference of which genetic variants will affect splicing and what their functional consequences will be is challenging, particularly for variants outside of the essential splice sites. Here, we describe a set of data-driven heuristics that inform the interpretation of human splice-altering variants (SAVs) based on the analysis of annotated exons, experimentally validated SAVs, and the currently understood principles of splicing biology. We defined requisite splicing criteria by examining around 202,000 canonical protein-coding exons and 19,000 experimentally validated splicing branchpoints. This analysis defined the sequence, spacing, and motif strength required for splicing, with 95.9% of the exons examined meeting these criteria. By considering over 12,000 experimentally validated variants from the SpliceVarDB, we defined a set of heuristics that inform the evaluation of putative SAVs. To ensure the applicability of each heuristic, only those supported by at least 10 experimentally validated variants were considered. This allowed us to establish a measure of spliceogenicity: the proportion of variants at a location (or motif site) that affected splicing in a given context. This study makes considerable advances toward bridging the gap between computational predictions and the biological process of splicing, offering an evidence-based approach to identifying SAVs and evaluating their impact. Our splicing heuristics enhance the current framework for genetic variant evaluation with a robust, detailed, and comprehensible analysis by adding valuable context over traditional binary prediction tools.

Tailored antisense oligonucleotides designed to correct aberrant splicing reveal actionable groups of mutations for rare genetic disorders

Effective translation of rare disease diagnosis knowledge into therapeutic applications is achievable within a reasonable timeframe; where mutations are amenable to current antisense oligonucleotide technology. In our study, we identified five distinct types of abnormal splice-causing mutations in patients with rare genetic disorders and developed a tailored antisense oligonucleotide for each mutation type using phosphorodiamidate morpholino oligomers with or without octa-guanidine dendrimers and 2'-O-methoxyethyl phosphorothioate. We observed variations in treatment effects and efficiencies, influenced by both the chosen chemistry and the specific nature of the aberrant splicing patterns targeted for correction. Our study demonstrated the successful correction of all five different types of aberrant splicing. Our findings reveal that effective correction of aberrant splicing can depend on altering the chemical composition of oligonucleotides and suggest a fast, efficient, and feasible approach for developing personalized therapeutic interventions for genetic disorders within short time frames.

Natural languages and RNA virus evolution

Information concepts from physics, mathematics and computer science support many areas of research in biology. Their focus is on objective information, which provides correlations and patterns related to objects, processes, marks and signals. In these approaches only the quantitative aspects of the meaning of the information is relevant. In other areas of biology, 'meaningful information', which is subjective in nature, relies on the physiology of the organism's sensory organs and on the interpretation of the perceived signals, which is then translated into action, even if this is only mental (in brained animals). Information is involved, in terms of both amount and quality. Here we contextualize and review the main theories that deal with 'meaningful-information' at a molecular level from different areas of natural language research, namely biosemiotics, code-biology, biocommunication and biohermeneutics. As this information mediates between the organism and its environment, we emphasize how such theories compare with the neo-Darwinian treatment of genetic information, and how they project onto the rapid evolution of RNA viruses.

Regulating PCCA gene expression by modulation of pseudoexon splicing patterns to rescue enzyme activity in propionic acidemia

Pseudoexons are nonfunctional intronic sequences that can be activated by deep-intronic sequence variation. Activation increases pseudoexon inclusion in mRNA and interferes with normal gene expression. The PCCA c.1285-1416A>G variation activates a pseudoexon and causes the severe metabolic disorder propionic acidemia by deficiency of the propionyl-CoA carboxylase enzyme encoded by PCCA and PCCB. We characterized this pathogenic pseudoexon activation event in detail and identified hnRNP A1 to be important for normal repression. The PCCA c.1285-1416A>G variation disrupts an hnRNP A1-binding splicing silencer and simultaneously creates a splicing enhancer. We demonstrate that blocking this region of regulation with splice-switching antisense oligonucleotides restores normal splicing and rescues enzyme activity in patient fibroblasts and in a cellular model created by CRISPR gene editing. Interestingly, the PCCA pseudoexon offers an unexploited potential to upregulate gene expression because healthy tissues show relatively high inclusion levels. By blocking inclusion of the nonactivated wild-type pseudoexon, we can increase both PCCA and PCCB protein levels, which increases the activity of the heterododecameric enzyme. Surprisingly, we can increase enzyme activity from residual levels in not only patient fibroblasts harboring PCCA missense variants but also those harboring PCCB missense variants. This is a potential treatment strategy for propionic acidemia.

The tri-flow adaptiveness of codes in major evolutionary transitions

Life codes increase in both number and variety with biological complexity. Although our knowledge of codes is constantly expanding, the evolutionary progression of organic, neural, and cultural codes in response to selection pressure remains poorly understood. Greater clarification of the selective mechanisms is achieved by investigating how major evolutionary transitions reduce spatiotemporal and energetic constraints on transmitting heritable code to offspring. Evolution toward less constrained flows is integral to enduring flow architecture everywhere, in both engineered and natural flow systems. Beginning approximately 4 billion years ago, the most basic level for transmitting genetic material to offspring was initiated by protocell division. Evidence from ribosomes suggests that protocells transmitted comma-free or circular codes, preceding the evolution of standard genetic code. This rudimentary information flow within protocells is likely to have first emerged within the geo-energetic and geospatial constraints of hydrothermal vents. A broad-gauged hypothesis is that major evolutionary transitions overcame such constraints with tri-flow adaptations. The interconnected triple flows incorporated energy-converting, spatiotemporal, and code-based informational dynamics. Such tri-flow adaptations stacked sequence splicing code on top of protein-DNA recognition code in eukaryotes, prefiguring the transition to sexual reproduction. Sex overcame the spatiotemporal-energetic constraints of binary fission with further code stacking. Examples are tubulin code and transcription initiation code in vertebrates. In a later evolutionary transition, language reduced metabolic-spatiotemporal constraints on inheritance by stacking phonetic, phonological, and orthographic codes. In organisms that reproduce sexually, each major evolutionary transition is shown to be a tri-flow adaptation that adds new levels of code-based informational exchange. Evolving biological complexity is also shown to increase the nongenetic transmissibility of code.

Parvovirus B19 and Human Parvovirus 4 Encode Similar Proteins in a Reading Frame Overlapping the VP1 Capsid Gene

Viruses frequently contain overlapping genes, which encode functionally unrelated proteins from the same DNA or RNA region but in different reading frames. Yet, overlapping genes are often overlooked during genome annotation, in particular in DNA viruses. Here we looked for the presence of overlapping genes likely to encode a functional protein in human parvovirus B19 (genus Erythroparvovirus), using an experimentally validated software, Synplot2. Synplot2 detected an open reading frame, X, conserved in all erythroparvoviruses, which overlaps the VP1 capsid gene and is under highly significant selection pressure. In a related virus, human parvovirus 4 (genus Tetraparvovirus), Synplot2 also detected an open reading frame under highly significant selection pressure, ARF1, which overlaps the VP1 gene and is conserved in all tetraparvoviruses. These findings provide compelling evidence that the X and ARF1 proteins must be expressed and functional. X and ARF1 have the exact same location (they overlap the region of the VP1 gene encoding the phospholipase A2 domain), are both in the same frame (+1) with respect to the VP1 frame, and encode proteins with similar predicted properties, including a central transmembrane region. Further studies will be needed to determine whether they have a common origin and similar function. X and ARF1 are probably translated either from a polycistronic mRNA by a non-canonical mechanism, or from an unmapped monocistronic mRNA. Finally, we also discovered proteins predicted to be expressed from a frame overlapping VP1 in other species related to parvovirus B19: porcine parvovirus 2 (Z protein) and bovine parvovirus 3 (X-like protein).

A novel HPV16 splicing enhancer critical for viral oncogene expression and cell immortalization

High-risk carcinogenic human papillomaviruses (HPVs), e.g. HPV16, express the E6 and E7 oncogenes from two mRNAs that are generated in a mutually exclusive manner by splicing. The HPV16 E7 mRNA, also known as the E6*I/E7 mRNA, is produced by splicing between splice sites SD226 and SA409, while E6 mRNAs retain the intron between these splice sites. We show that splicing between HPV16 splice sites SD226 and SA409 is controlled by a splicing enhancer consisting of a perfect repeat of an adenosine-rich, 11 nucleotide sequence: AAAAGCAAAGA. Two nucleotide substitutions in both 11 nucleotide sequences specifically inhibited production of the spliced E6*I/E7 mRNA. As a result, production of E7 protein was reduced and the ability of HPV16 to immortalize human primary keratinocytes was abolished. The splicing-enhancing effect was mediated by the cellular TRAP150/THRAP3 protein that also enhanced splicing of other high-risk HPV E6*I/E7 mRNAs, but had no effect on low-risk HPV mRNAs. In summary, we have identified a novel splicing enhancer in the E6 coding region that is specific for high-risk HPVs and that is critically linked to HPV16 carcinogenic properties.

Hold out the genome: a roadmap to solving the cis-regulatory code

Gene expression is regulated by transcription factors that work together to read cis-regulatory DNA sequences. The 'cis-regulatory code' - how cells interpret DNA sequences to determine when, where and how much genes should be expressed - has proven to be exceedingly complex. Recently, advances in the scale and resolution of functional genomics assays and machine learning have enabled substantial progress towards deciphering this code. However, the cis-regulatory code will probably never be solved if models are trained only on genomic sequences; regions of homology can easily lead to overestimation of predictive performance, and our genome is too short and has insufficient sequence diversity to learn all relevant parameters. Fortunately, randomly synthesized DNA sequences enable testing a far larger sequence space than exists in our genomes, and designed DNA sequences enable targeted queries to maximally improve the models. As the same biochemical principles are used to interpret DNA regardless of its source, models trained on these synthetic data can predict genomic activity, often better than genome-trained models. Here we provide an outlook on the field, and propose a roadmap towards solving the cis-regulatory code by a combination of machine learning and massively parallel assays using synthetic DNA.

How ambiguity codes specify molecular descriptors and information flow in Code Biology

The article presents IUPAC ambiguity codes for incomplete nucleic acid specification, and their use in Code Biology. It is shown how to use this nomenclature in order to extract accurate information on different properties of the biological systems. We investigated the use of ambiguity codes, as mathematical and logical operators and truth table elements, for the encoding of amino acids by means of the Standard Genetic Code. It is explained how to use ambiguity codes and truth functions in order to obtain accurate information on different properties of the biological systems. Nucleotide ambiguity codes could be applied to: 1. encoding descriptive information of nucleotides, amino acids and proteins (e.g., of polarity, relative solvent accessibility, atom depth, etc.), and 2. system modelling ranging from standard bioinformatics tools to classic evolutionary models (i.e. from Miyazawa-Jernigan statistical potential to Kimura three-substitution-type model, respectively). It is shown that the algorithms based on IUPAC ambiguity codes, Boolean functions and truth table, Probabilistic Square of Opposition/Semiotic Square and Klein 4-groups-could be used for the bioinformatics analyses and Relational data modelling in natural science. Underlying mathematical, logical and semiotic concepts of interest are presented and addressed.

How biological codes break causal chains to enable autonomy for organisms

Autonomy, meaning freedom from exogenous control, requires independence of both constitution and cybernetic regulation. Here, the necessity of biological codes to achieve both is explained, assuming that Aristotelian efficient cause is 'formal cause empowered by physical force'. Constitutive independence requires closure to efficient causation (in the Rosen sense); cybernetic independence requires transformation of cause-effect into signal-response relations at the organism boundary; the combination of both kinds of independence enables adaptation and evolution. Codes and cyphers translate information from one form of physical embodiment (domain) to another. Because information can only contribute as formal cause to efficient cause within the domain of its embodiment, translation can extend or restrict the range over which information is effective. Closure to efficient causation requires internalised information to be isolated from the cycle of efficient causes that it informs: e.g. Von Neumann self-replicator requires a (template) source of information that is causally isolated from the physical replication system. Life operationalises this isolation with the genetic code translating from the (isolated) domain of codons to that of protein interactions. Separately, cybernetic freedom is achieved at the cell boundary because transducers, which embody molecular coding, translate exogenous information into a domain where it no longer has the power of efficient cause. Information, not efficient cause, passes through the boundary to serve as stimulus for an internally generated response. Coding further extends freedom by enabling historically accumulated information to be selectively transformed into efficient cause under internal control, leaving it otherwise stored inactive. Code-based translation thus enables selective causal isolation, controlling the flow from cause to effect. Genetic code, cell-signalling codes and, in eukaryotes, the histone code, signal sequence based protein sorting and other code-dependent processes all regulate and separate causal chains. The existence of life can be seen as an expression of the power of molecular codes to selectively isolate and thereby organise causal relations among molecular interactions to form an organism.

Tnpo3 controls splicing of the pre-mRNA encoding the canonical TCR α chain of iNKT cells

Unconventional T cells, such as innate natural killer T cells (iNKT) cells, are an important part of vertebrate immune defences. iNKT recognise glycolipids through a T cell receptor (TCR) that is composed of a semi-invariant TCR α chain, paired with a restricted set of TCR β chains. Here, we show that splicing of the cognate Trav11-Traj18-Trac pre-mRNA encoding the characteristic Vα14Jα18 variable region of this semi-invariant TCR depends on the presence of Tnpo3. The Tnpo3 gene encodes a nuclear transporter of the β-karyopherin family whose cargo includes various splice regulators. The block of iNKT cell development in the absence of Tnpo3 can be overcome by transgenic provision of a rearranged Trav11-Traj18-Trac cDNA, indicating that Tnpo3 deficiency does not interfere with the development of iNKT cells per se. Our study thus identifies a role for Tnpo3 in regulating the splicing of the pre-mRNA encoding the cognate TCRα chain of iNKT cells.

Alternative splicing outcomes across an RNA-binding protein concentration gradient

Alternative splicing (AS) is a dynamic RNA processing step that produces multiple RNA isoforms from a single pre-mRNA transcript and contributes to the complexity of the cellular transcriptome and proteome. This process is regulated through a network of cis-regulatory sequence elements and trans-acting factors, most-notably RNA binding proteins (RBPs). The muscleblind-like (MBNL) and RNA binding fox-1 homolog (RBFOX) are two well characterized families of RBPs that regulate fetal to adult AS transitions critical for proper muscle, heart, and central nervous system development. To better understand how the concentration of these RBPs influences AS transcriptome wide, we engineered a MBNL1 and RBFOX1 inducible HEK-293 cell line. Modest induction of exogenous RBFOX1 in this cell line modulated MBNL1-dependent AS outcomes in 3 skipped exon events, despite significant levels of endogenous RBFOX1 and RBFOX2. Due to background RBFOX levels, we conducted a focused analysis of dose-dependent MBNL1 skipped exon AS outcomes and generated transcriptome wide dose-response curves. Analysis of this data demonstrates that MBNL1-regulated exclusion events may require higher concentrations of MBNL1 protein to properly regulate AS outcomes compared to inclusion events and that multiple arrangements of YGCY motifs can produce similar splicing outcomes. These results suggest that rather than a simple relationship between the organization of RBP binding sites and a specific splicing outcome, that complex interaction networks govern both AS inclusion and exclusion events across a RBP gradient.

Alternative Splicing Changes Promoted by NOVA2 Upregulation in Endothelial Cells and Relevance for Gastric Cancer

Angiogenesis is crucial for cancer progression. While several anti-angiogenic drugs are in use for cancer treatment, their clinical benefits are unsatisfactory. Thus, a deeper understanding of the mechanisms sustaining cancer vessel growth is fundamental to identify novel biomarkers and therapeutic targets. Alternative splicing (AS) is an essential modifier of human proteome diversity. Nevertheless, AS contribution to tumor vasculature development is poorly known. The Neuro-Oncological Ventral Antigen 2 (NOVA2) is a critical AS regulator of angiogenesis and vascular development. NOVA2 is upregulated in tumor endothelial cells (ECs) of different cancers, thus representing a potential driver of tumor blood vessel aberrancies. Here, we identified novel AS transcripts generated upon NOVA2 upregulation in ECs, suggesting a pervasive role of NOVA2 in vascular biology. In addition, we report that NOVA2 is also upregulated in ECs of gastric cancer (GC), and its expression correlates with poor overall survival of GC patients. Finally, we found that the AS of the Rap Guanine Nucleotide Exchange Factor 6 (RapGEF6), a newly identified NOVA2 target, is altered in GC patients and associated with NOVA2 expression, tumor angiogenesis, and poor patient outcome. Our findings provide a better understanding of GC biology and suggest that AS might be exploited to identify novel biomarkers and therapeutics for anti-angiogenic GC treatments.

Minigene splicing assays reveal new insights into exonic variants of the SLC12A3 gene in Gitelman syndrome

BACKGROUND

Gitelman syndrome (GS) is a type of salt-losing tubular disease, most of which is caused by SLC12A3 gene variants, and missense variants account for the majority. Recently, the phenomenon of exon skipping, in which variants disrupt normal pre-mRNA splicing, has been related to a variety of diseases. Therefore, we hypothesize that a certain proportion of SLC12A3 variants can result in disease via interfering with the normal splicing process.

METHODS

We analyzed 342 previously presumed SLC12A3 missense variants using bioinformatics programs and identified candidate variants that may alter the splicing of pre-mRNA through minigene assays.

RESULTS

Our study revealed that, among ten candidate variants, six variants (c.602G>A, c.602G>T, c.1667C>T, c.1925G>A, c.2548G>C, and c.2549G>C) led to complete or incomplete exon skipping by affecting exonic splicing regulatory elements and/or disturbing canonical splice sites.

CONCLUSION

It is worth mentioning that this is the largest study on pre-mRNA splicing of SLC12A3 exonic variants. In addition, our study emphasizes the importance of detecting splicing function at the mRNA level in GS and indicates that minigene analysis is a valuable tool for splicing functional assays of variants in vitro.

A history of the genetic and molecular identification of genes and their functions controlling insect sex determination

The genetics of the sex determination regulatory cascade in Drosophila melanogaster has a fascinating history, interlinked with the foundation of the Genetics discipline itself. The discovery that alternative splicing rather than differential transcription is the molecular mechanism underlying the upstream control of sex differences in the Drosophila model system was surprising. This notion is now fully integrated into the scientific canon, appearing in many genetics textbooks and online education resources. In the last three decades, it was a key reference point for starting evolutionary studies in other insect species by using homology-based approaches. This review will introduce a very brief history of Drosophila genetics. It will describe the genetic and molecular approaches applied for the identifying and cloning key genes involved in sex determination in Drosophila and in many other insect species. These comparative analyses led to supporting the idea that sex-determining pathways have evolved mainly by recruiting different upstream signals/genes while maintaining widely conserved intermediate and downstream regulatory genes. The review also provides examples of the link between technological advances and research achievements, to stimulate reflections on how science is produced. It aims to hopefully strengthen the related historical and conceptual knowledge of general readers of other disciplines and of younger geneticists, often focused on the latest technical-molecular approaches.

SRSF1 acts as an IFN-I-regulated cellular dependency factor decisively affecting HIV-1 post-integration steps

Efficient HIV-1 replication depends on balanced levels of host cell components including cellular splicing factors as the family of serine/arginine-rich splicing factors (SRSF, 1-10). Type I interferons (IFN-I) play a crucial role in the innate immunity against HIV-1 by inducing the expression of IFN-stimulated genes (ISGs) including potent host restriction factors. The less well known IFN-repressed genes (IRepGs) might additionally affect viral replication by downregulating host dependency factors that are essential for the viral life cycle; however, so far, the knowledge about IRepGs involved in HIV-1 infection is very limited. In this work, we could demonstrate that HIV-1 infection and the associated ISG induction correlated with low SRSF1 levels in intestinal lamina propria mononuclear cells (LPMCs) and peripheral blood mononuclear cells (PBMCs) during acute and chronic HIV-1 infection. In HIV-1-susceptible cell lines as well as primary monocyte-derived macrophages (MDMs), expression levels of SRSF1 were transiently repressed upon treatment with specific IFNα subtypes in vitro. Mechanically, 4sU labeling of newly transcribed mRNAs revealed that IFN-mediated SRSF1 repression is regulated on early RNA level. SRSF1 knockdown led to an increase in total viral RNA levels, but the relative proportion of the HIV-1 viral infectivity factor (Vif) coding transcripts, which is essential to counteract APOBEC3G-mediated host restriction, was significantly reduced. In the presence of high APOBEC3G levels, however, increased LTR activity upon SRSF1 knockdown facilitated the overall replication, despite decreased vif mRNA levels. In contrast, SRSF1 overexpression significantly impaired HIV-1 post-integration steps including LTR transcription, alternative splice site usage, and virus particle production. Since balanced SRSF1 levels are crucial for efficient viral replication, our data highlight the so far undescribed role of SRSF1 acting as an IFN-modulated cellular dependency factor decisively regulating HIV-1 post-integration steps.

Cooperative engagement and subsequent selective displacement of SR proteins define the pre-mRNA 3D structural scaffold for early spliceosome assembly

We recently reported that serine–arginine-rich (SR) protein-mediated pre-mRNA structural remodeling generates a pre-mRNA 3D structural scaffold that is stably recognized by the early spliceosomal components. However, the intermediate steps between the free pre-mRNA and the assembled early spliceosome are not yet characterized. By probing the early spliceosomal complexes in vitro and RNA-protein interactions in vivo, we show that the SR proteins bind the pre-mRNAs cooperatively generating a substrate that recruits U1 snRNP and U2AF65 in a splice signal-independent manner. Excess U1 snRNP selectively displaces some of the SR protein molecules from the pre-mRNA generating the substrate for splice signal-specific, sequential recognition by U1 snRNP, U2AF65 and U2AF35. Our work thus identifies a novel function of U1 snRNP in mammalian splicing substrate definition, explains the need for excess U1 snRNP compared to other U snRNPs in vivo, demonstrates how excess SR proteins could inhibit splicing, and provides a conceptual basis to examine if this mechanism of splicing substrate definition is employed by other splicing regulatory proteins.

Uncovering the impacts of alternative splicing on the proteome with current omics techniques

The high-throughput sequencing of cellular RNAs has underscored a broad effect of isoform diversification through alternative splicing on the transcriptome. Moreover, the differential production of transcript isoforms from gene loci has been recognized as a critical mechanism in cell differentiation, organismal development, and disease. Yet, the extent of the impact of alternative splicing on protein production and cellular function remains a matter of debate. Multiple experimental and computational approaches have been developed in recent years to address this question. These studies have unveiled how molecular changes at different steps in the RNA processing pathway can lead to differences in protein production and have functional effects. New and emerging experimental technologies open exciting new opportunities to develop new methods to fully establish the connection between messenger RNA expression and protein production and to further investigate how RNA variation impacts the proteome and cell function. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing Translation > Regulation RNA Evolution and Genomics > Computational Analyses of RNA.

Identifying genes with conserved splicing structure and orthologous isoforms in human, mouse and dog

BACKGROUND

In eukaryote transcriptomes, a significant amount of transcript diversity comes from genes' capacity to generate different transcripts through alternative splicing. Identifying orthologous alternative transcripts across multiple species is of particular interest for genome annotators. However, there is no formal definition of transcript orthology based on the splicing structure conservation. Likewise there is no public dataset benchmark providing groups of orthologous transcripts sharing a conserved splicing structure.

RESULTS

We introduced a formal definition of splicing structure orthology and we predicted transcript orthologs in human, mouse and dog. Applying a selective strategy, we analyzed 2,167 genes and their 18,109 known transcripts and identified a set of 253 gene orthologs that shared a conserved splicing structure in all three species. We predicted 6,861 transcript CDSs (coding sequence), mainly for dog, an emergent model species. Each predicted transcript was an ortholog of a known transcript: both share the same CDS splicing structure. Evidence for the existence of the predicted CDSs was found in external data.

CONCLUSIONS

We generated a dataset of 253 gene triplets, structurally conserved and sharing all their CDSs in human, mouse and dog, which correspond to 879 triplets of spliced CDS orthologs. We have released the dataset both as an SQL database and as tabulated files. The data consists of the 879 CDS orthology groups with their detailed splicing structures, and the predicted CDSs, associated with their experimental evidence. The 6,861 predicted CDSs are provided in GTF files. Our data may contribute to compare highly conserved genes across three species, for comparative transcriptomics at the isoform level, or for benchmarking splice aligners and methods focusing on the identification of splicing orthologs. The data is available at https://data-access.cesgo.org/index.php/s/V97GXxOS66NqTkZ .

OTC intron 4 variations mediate pathogenic splicing patterns caused by the c.386G>A mutation in humans and spfash mice, and govern susceptibility to RNA-based therapies

Background

Aberrant splicing is a common outcome in the presence of exonic or intronic variants that might hamper the intricate network of interactions defining an exon in a specific gene context. Therefore, the evaluation of the functional, and potentially pathological, role of nucleotide changes remains one of the major challenges in the modern genomic era. This aspect has also to be taken into account during the pre-clinical evaluation of innovative therapeutic approaches in animal models of human diseases. This is of particular relevance when developing therapeutics acting on splicing, an intriguing and expanding research area for several disorders. Here, we addressed species-specific splicing mechanisms triggered by the OTC c.386G>A mutation, relatively frequent in humans, leading to Ornithine TransCarbamylase Deficiency (OTCD) in patients and spf^ash mice, and its differential susceptibility to RNA therapeutics based on engineered U1snRNA.

Methods

Creation and co-expression of engineered U1snRNAs with human and mouse minigenes, either wild-type or harbouring different nucleotide changes, in human (HepG2) and mouse (Hepa1-6) hepatoma cells followed by analysis of splicing pattern. RNA pulldown studies to evaluate binding of specific splicing factors.

Results

Comparative nucleotide analysis suggested a role for the intronic +10-11 nucleotides, and pull-down assays showed that they confer preferential binding to the TIA1 splicing factor in the mouse context, where TIA1 overexpression further increases correct splicing. Consistently, the splicing profile of the human minigene with mouse +10-11 nucleotides overlapped that of mouse minigene, and restored responsiveness to TIA1 overexpression and to compensatory U1snRNA. Swapping the human +10-11 nucleotides into the mouse context had opposite effects.

Moreover, the interplay between the authentic and the adjacent cryptic 5′ss in the human OTC dictates pathogenic mechanisms of several OTCD-causing 5′ss mutations, and only the c.386+5G>A change, abrogating the cryptic 5′ss, was rescuable by engineered U1snRNA.

Conclusions

Subtle intronic variations explain species-specific OTC splicing patterns driven by the c.386G>A mutation, and the responsiveness to engineered U1snRNAs, which suggests careful elucidation of molecular mechanisms before proposing translation of tailored therapeutics from animal models to humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s10020-021-00418-9.

Alternative splicing and liver disease

Alternative splicing produces complex and dynamic changes in the protein isoforms that are necessary for the proper biological functioning of the metabolic pathways involved in liver development and hepatocyte homeostasis. Changes in the physiological state of alternatively spliced forms are increasingly linked to liver pathologies. This may occur when the expression or function of the set of proteins controlling the alternative splicing processes are altered by external effectors such as oxidative stress and other environmental variations. Studies addressing these modifications reveal a complex interplay between the expression levels of different proteins that regulate the alternative splicing process as well as the changes in alternative splicing. This interplay results in a cascade of different protein isoforms that correlate with the progression of non-alcoholic fatty liver disease, hepatocellular carcinoma, and alcoholic liver disease. However, research on the detailed molecular mechanism underlying the production of these isoforms is needed. It is imperative to identify the physiological processes affected by the differentially spliced isoforms and confirm their role on the onset and maintenance of the pathology. This is required to design potential therapeutic approaches targeting the key splicing changes to revert the pathological condition as well as identify prognostic markers. In this review, we describe the complexity of the splicing process through an example to encourage researchers to go down this path. Subsequently, rather than a catalog of splicing events we have hand-picked and discuss a few selected studies of specific liver pathologies and suggested ways to focus research on these areas.

Pseudoexon activation in disease by non-splice site deep intronic sequence variation - wild type pseudoexons constitute high-risk sites in the human genome

Accuracy of pre-messenger RNA (pre-mRNA) splicing is crucial for normal gene expression. Complex regulation supports the spliceosomal distinction between authentic exons and the many seemingly functional splice sites delimiting pseudoexons. Pseudoexons are nonfunctional intronic sequences that can be activated for aberrant inclusion in mRNA, which may cause disease. Pseudoexon activation is very challenging to predict, in particular when activation occurs by sequence variants that alter the splicing regulatory environment without directly affecting splice sites. As pseudoexon inclusion often evades detection due to activation of nonsense-mediated mRNA decay, and because conventional diagnostic procedures miss deep intronic sequence variation, pseudoexon activation is a heavily underreported disease mechanism. Pseudoexon characteristics have mainly been studied based on in silico predicted sequences. Moreover, because recognition of sequence variants that create or strengthen splice sites is possible by comparison with well-established consensus sequences, this type of pseudoexon activation is by far the most frequently reported. Here we review all known human disease-associated pseudoexons that carry functional splice sites and are activated by deep intronic sequence variants located outside splice site sequences. We delineate common characteristics that make this type of wild type pseudoexons distinct high-risk sites in the human genome.

The major evolutionary transitions and codes of life

Major evolutionary transitions as well as the evolution of codes of life are key elements in macroevolution which are characterized by increase in complexity Major evolutionary transitions ensues by a transition in individuality and by the evolution of a novel mode of using, transmitting or storing information. Here is where codes of life enter the picture: they are arbitrary mappings between different (mostly) molecular species. This flexibility allows information to be employed in a variety of ways, which can fuel evolutionary innovation. The collation of the list of major evolutionary transitions and the list of codes of life show a clear pattern: codes evolved prior to a major evolutionary transition and then played roles in the transition and/or in the transformation of the new individual. The evolution of a new code of life is in itself not a major evolutionary transition but allow major evolutionary transitions to happen. This could help us to identify new organic codes.

Discovery of a pre-mRNA structural scaffold as a contributor to the mammalian splicing code

The specific recognition of splice signals at or near exon-intron junctions is not explained by their weak conservation and instead is postulated to require a multitude of features embedded in the pre-mRNA strand. We explored the possibility of 3D structural scaffold of AdML-a model pre-mRNA substrate-guiding early spliceosomal components to the splice signal sequences. We find that mutations in the non-cognate splice signal sequences impede recruitment of early spliceosomal components due to disruption of the global structure of the pre-mRNA. We further find that the pre-mRNA segments potentially interacting with the early spliceosomal component U1 snRNP are distributed across the intron, that there is a spatial proximity of 5' and 3' splice sites within the pre-mRNA scaffold, and that an interplay exists between the structural scaffold and splicing regulatory elements in recruiting early spliceosomal components. These results suggest that early spliceosomal components can recognize a 3D structural scaffold beyond the short splice signal sequences, and that in our model pre-mRNA, this scaffold is formed across the intron involving the major splice signals. This provides a conceptual basis to analyze the contribution of recognizable 3D structural scaffolds to the splicing code across the mammalian transcriptome.

Removing quote marks from the RNA polymerase II CTD 'code'

In eukaryotes, RNA polymerase II (Pol II) is responsible for the synthesis of all mRNAs and myriads of short and long untranslated RNAs, whose fabrication involves close spatiotemporal coordination between transcription, RNA processing and chromatin modification. Crucial for such a coordination is an unusual C-terminal domain (CTD) of the Pol II largest subunit, made of tandem repetitions (26 in yeast, 52 in chordates) of the heptapeptide with the consensus sequence YSPTSPS. Although largely unstructured and with poor sequence content, the Pol II CTD derives its extraordinary functional versatility from the fact that each amino acid in the heptapeptide can be posttranslationally modified, and that different combinations of CTD covalent marks are specifically recognized by different protein binding partners. These features have led to propose the existence of a Pol II CTD code, but this expression is generally used by authors with some caution, revealed by the frequent use of quote marks for the word 'code'. Based on the theoretical framework of code biology, it is argued here that the Pol II CTD modification system meets the requirements of a true organic code, where different CTD modification states represent organic signs whose organic meanings are biological reactions contributing to the many facets of RNA biogenesis in coordination with RNA synthesis by Pol II. Importantly, the Pol II CTD code is instantiated by adaptor proteins possessing at least two distinct domains, one of which devoted to specific recognition of CTD modification profiles. Furthermore, code rules can be altered by experimental interchange of CTD recognition domains of different adaptor proteins, a fact arguing in favor of the arbitrariness, and thus bona fide character, of the Pol II CTD code. Since the growing family of CTD adaptors includes RNA binding proteins and histone modification complexes, the Pol II CTD code is by its nature integrated with other organic codes, in particular the splicing code and the histone code. These issues will be discussed taking into account fascinating developments in Pol II CTD research, like the discovery of novel modifications at non-consensus sites, the recently recognized CTD physicochemical properties favoring liquid-liquid phase separation, and the discovery that the Pol II CTD, originated before the divergence of most extant eukaryotic taxa, has expanded and diversified with developmental complexity in animals and plants.

The drawbridge of nature: Evolutionary complexification as a generation and novel interpretation of coding systems

The phenomenon of evolutionary complexification corresponds to the generation of new coding systems (defined as а codepoiesis by Marcello Barbieri). The whole process of generating novel coding statements that substantiate organizational complexification leads to an expansion of the system that incorporates externality to support newly generated complex structures. During complexifying evolution, the values are assigned to the previously unproven statements via their encoding by using new codes or rearranging the old ones. In this perspective, living systems during evolution continuously realize the proof of Gödel's theorem. In the real physical world, this realization is grounded in the irreversible reduction of the fundamental uncertainty appearing in the self-referential process of internal measurement performed by living systems. It leads to the formation of reflexive loops that establish novel interrelations between the biosystem and the external world and provide a possibility of active anticipatory transformation of externality. We propose a metamathematical framework that can account for the underlying logic of codepoiesis, outline the basic principles of the generation of new coding systems, and describe main codepoietic events in the course of progressive biological evolution. The evolutionary complexification is viewed as a metasystem transition that results in the increase of external work by the system based on the division of labor between its components.

Insights into the role of alternative splicing in plant temperature response

Alternative splicing occurs in all eukaryotic organisms. Since the first description of multiexon genes and the splicing machinery, the field has expanded rapidly, especially in animals and yeast. However, our knowledge about splicing in plants is still quite fragmented. Though eukaryotes show some similarity in the composition and dynamics of the splicing machinery, observations of unique plant traits are only starting to emerge. For instance, plant alternative splicing is closely linked to their ability to perceive various environmental stimuli. Due to their sessile lifestyle, temperature is a central source of information allowing plants to adjust their development to match current growth conditions. Hence, seasonal temperature fluctuations and day-night cycles can strongly influence plant morphology across developmental stages. Here we discuss the available data about temperature-dependent alternative splicing in plants. Given its fragmented state it is not always possible to fit specific observations into a coherent picture, yet it is sufficient to estimate the complexity of this field and the need of further research. Better understanding of alternative splicing as a part of plant temperature response and adaptation may also prove to be a powerful tool for both, fundamental and applied sciences.

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.

Potential role of the X circular code in the regulation of gene expression

The X circular code is a set of 20 trinucleotides (codons) that has been identified in the protein-coding genes of most organisms (bacteria, archaea, eukaryotes, plasmids, viruses). It has been shown previously that the X circular code has the important mathematical property of being an error-correcting code. Thus, motifs of the X circular code, i.e. a series of codons belonging to X and called X motifs, allow identification and maintenance of the reading frame in genes. X motifs are significantly enriched in protein-coding genes, but have also been identified in many transfer RNA (tRNA) genes and in important functional regions of the ribosomal RNA (rRNA), notably in the peptidyl transferase center and the decoding center. Here, we investigate the potential role of X motifs as functional elements of protein-coding genes. First, we identify the codons of the X circular code which are frequent or rare in each domain of life (archaea, bacteria, eukaryota) and show that, for the amino acids with the highest codon bias, the preferred codon is often an X codon. We also observe a correlation between the 20 X codons and the optimal codons/dicodons that have been shown to influence translation efficiency. Then, we examined recently published experimental results concerning gene expression levels in diverse organisms. The approach used is the analysis of X motifs according to their density d_s(X), i.e. the number of X motifs per kilobase in a gene sequence s. Surprisingly, this simple parameter identifies several unexpected relations between the X circular code and gene expression. For example, the X motifs are significantly enriched in the minimal gene set belonging to the three domains of life, and in codon-optimized genes. Furthermore, the density of X motifs generally correlates with experimental measures of translation efficiency and mRNA stability. Taken together, these results lead us to propose that the X motifs may represent a genetic signal contributing to the maintenance of the correct reading frame and the optimization and regulation of gene expression.

Six Exonic Variants in the SLC5A2 Gene Cause Exon Skipping in a Minigene Assay

Background

Familial renal glucosuria is a rare renal tubular disorder caused by SLC5A2 gene variants. Most of them are exonic variants and have been classified as missense variants. However, there is growing evidence that some of these variants can be detrimental by affecting the pre-mRNA splicing process. Therefore, we hypothesize that a certain proportion of SLC5A2 exonic variants can result in disease via interfering with the normal splicing process of the pre-mRNA.

Methods

We used bioinformatics programs to analyze 77 previously described presumed SLC5A2 missense variants and identified candidate variants that may alter the splicing of pre-mRNA through minigene assays.

Results

Our study indicated six of 7 candidate variants induced splicing alterations. Variants c.216C > A, c.294C > A, c.886G > C, c.932A > G and c.962A > G may disrupt splicing enhancer motifs and generate splicing silencer sequences resulting in the skipping of exon 3. Variants c.305C > T and c.1129G > A probably disturb splice sites leading to exon skipping.

Conclusion

To our knowledge, we report, for the first time, SLC5A2 exonic variants that produce alterations in pre-mRNA. Our research reinforces the importance of assessing the consequences for putative point variants at the mRNA level. Additionally, we propose that minigenes function analysis may be valuable to evaluate the impact of SLC5A2 exonic variants on pre-mRNA splicing without patients’ RNA samples.

A global survey of the transcriptome of allopolyploid Brassica napus based on single-molecule long-read isoform sequencing and Illumina-based RNA sequencing data

Brassica napus is a recent allopolyploid derived from the hybridization of Brassica rapa (A_r A_r ) and Brassica oleracea (C_o C_o ). Because of the high sequence similarity between the A_n and C_n subgenomes, it is difficult to provide an accurate landscape of the whole transcriptome of B. napus. To overcome this problem, we applied a single-molecule long-read isoform sequencing (Iso-Seq) technique that can produce long reads to explore the complex transcriptome of B. napus at the isoform level. From the Iso-Seq data, we obtained 147 698 non-redundant isoforms, capturing 37 403 annotated genes. A total of 18.1% (14 934/82 367) of the multi-exonic genes showed alternative splicing (AS). In addition, we identified 549 long non-coding RNAs, the majority of which displayed tissue-specific expression profiles, and detected 7742 annotated genes that possessed isoforms containing alternative polyadenylation sites. Moreover, 31 591 AS events located in open reading frames (ORFs) lead to potential protein isoforms by in-frame or frameshift changes in the ORF. Illumina RNA sequencing of five tissues that were pooled for Iso-Seq was also performed and showed that 69% of the AS events were tissue-specific. Our data provide abundant transcriptome resources for a transcript isoform catalog of B. napus, which will facilitate genome reannotation, strengthen our understanding of the B. napus transcriptome and be applied for further functional genomic research.

Blastocladiella emersonii spliceosome is regulated in response to the splicing inhibition caused by the metals cadmium, cobalt and manganese

Blastocladiella emersonii is an aquatic fungus of the phylum Blastocladiomycota, localized near the base of the fungal tree. Previous studies have shown that B. emersonii responds to heat shock and cadmium exposure inducing the transcription of a high number of genes. EST sequencing from heat shocked and cadmium exposed B. emersonii cells has shown that exposure to cadmium causes strong splicing inhibition. Despite the knowledge about splicing inhibition by cadmium, it is still unclear if other metal contaminants can cause the same response. In the present study, we have demonstrated that the effect of cadmium exposure on splicing inhibition is much stronger than that of other divalent metals such as cobalt and manganese. Data presented here also indicate that intron retention occurs randomly among the fungal transcripts, as verified by analyzing differently affected transcripts. In addition, we identified in the genome of B. emersonii the genes encoding the snRNA splicing components U1, U2, U4, U5 and U6 and observed that spliceosome snRNAs are upregulated in the presence of metals, in particular snRNA U1 in cells under cadmium exposure. This observation suggests that snRNA upregulation might be a defense of the fungal cell against the metal stress condition.

Croes C, Petersen N, Bosma P, Riccardi F, Pagani F, Pinotti M, van de Graaf SFJ. A Compensatory U1snRNA Partially Rescues FAH Splicing and Protein Expression in a Splicing-Defective Mouse Model of Tyrosinemia Type I

The elucidation of aberrant splicing mechanisms, frequently associated with disease has led to the development of RNA therapeutics based on the U1snRNA, which is involved in 5' splice site (5'ss) recognition. Studies in cellular models have demonstrated that engineered U1snRNAs can rescue different splicing mutation types. However, the assessment of their correction potential in vivo is limited by the scarcity of animal models with the targetable splicing defects. Here, we challenged the U1snRNA in the FAH5961SB mouse model of hepatic fumarylacetoacetate hydrolase (FAH) deficiency (Hereditary Tyrosinemia type I, HT1) due to the FAH c.706G>A splicing mutation. Through minigene expression studies we selected a compensatory U1snRNA (U1F) that was able to rescue this mutation. Intriguingly, adeno-associated virus-mediated delivery of U1F (AAV8-U1F), but not of U1wt, partially rescued FAH splicing in mouse hepatocytes. Consistently, FAH protein was detectable only in the liver of AAV8-U1F treated mice, which displayed a slightly prolonged survival. Moreover, RNA sequencing revealed the negligible impact of the U1F on the splicing profile and overall gene expression, thus pointing toward gene specificity. These data provide early in vivo proof-of-principle of the correction potential of compensatory U1snRNAs in HTI and encourage further optimization on a therapeutic perspective, and translation to other splicing-defective forms of metabolic diseases.

5' and 3' splicing signals evolution in vertebrates: Analysis in a conserved gene family

The mitochondrial solute carrier genes (SLC25) are highly conserved during vertebrate evolution. In most SLC25 genes of zebrafish, chicken, mouse, and human, the introns are located at exactly superimposable positions. In these topographically corresponding introns we studied the composition of the initial and terminal hexanucleotides (5'ss and 3'ss) which are instrumental in splicing signaling, focusing on the evolutionary conservation/mutation dynamics of these genetically related sequences. At each position, the per cent conservation of zebrafish individual nucleotides in chicken, mouse and human is proportional to their percent frequency in zebrafish; furthermore, nucleotide mutations are biased in favor of the more represented nucleotides, thus compensating for those highly represented zebrafish nucleotides which have not been conserved. As a result of these evolutionary dynamics, the general nucleotide composition at each position has remained relatively conserved throughout vertebrates. At 5'ss, following the canonical GT, A and G are largely prevailing at position +3, A at +4 and G at +5 (GT[A/G]AGx). At 3'ss, T and C are largely prevailing at positions -6, -5 and -3, preceding the canonical intron terminal AG ([C/T] [C/T]x[C/T]AG). However, the actual composition of the tetranucleotides at 5' and 3' often does not conform to the above scheme. At 5'ss the more canonical sequence is completely expressed in 63% of cases and partially (2 or 1 matches) in 37 % of cases. At 3'ss the more canonical sequence is completely expressed in 71 % of cases and partially (2 or 1 matches) in 29 % of cases. The nucleotide conservation loss (nucleotide mutation) is higher in the evolution from fish to the last common ancestor of birds and mammals (58 %), then diminishes in the successive evolution steps up to the mammalian common ancestor (10 %), and becomes still lower at the divergence of rodents and primates (5 %).

Alternative Splicing and DNA Damage Response in Plants

Plants are exposed to a variety of abiotic and biotic stresses that may result in DNA damage. Endogenous processes - such as DNA replication, DNA recombination, respiration, or photosynthesis - are also a threat to DNA integrity. It is therefore essential to understand the strategies plants have developed for DNA damage detection, signaling, and repair. Alternative splicing (AS) is a key post-transcriptional process with a role in regulation of gene expression. Recent studies demonstrate that the majority of intron-containing genes in plants are alternatively spliced, highlighting the importance of AS in plant development and stress response. Not only does AS ensure a versatile proteome and influence the abundance and availability of proteins greatly, it has also emerged as an important player in the DNA damage response (DDR) in animals. Despite extensive studies of DDR carried out in plants, its regulation at the level of AS has not been comprehensively addressed. Here, we provide some insights into the interplay between AS and DDR in plants.

Assessment of the functional impact on the pre-mRNA splicing process of 28 nucleotide variants associated with Pompe disease in GAA exon 2 and their recovery using antisense technology

Glycogen storage disease II (GSDII), also called Pompe disease, is an autosomal recessive inherited disease caused by a defect in glycogen metabolism due to the deficiency of the enzyme acid alpha-glucosidase (GAA) responsible for its degradation. So far, more than 500 sequence variants of the GAA gene have been reported but their possible involvement on the pre-messenger RNA splicing mechanism has not been extensively studied. In this work, we have investigated, by an in vitro functional assay, all putative splicing variants within GAA exon 2 and flanking introns. Our results show that many variants falling in the canonical splice site or the exon can induce GAA exon 2 skipping. In these cases, therefore, therapeutic strategies aimed at restoring protein folding of partially active mutated GAA proteins might not be sufficient. Regarding this issue, we have tested the effect of antisense oligonucleotides (AMOs) that were previously shown capable of rescuing splicing misregulation caused by the common c.-32-13T>G variant associated with the childhood/adult phenotype of GSDII. Interestingly, our results show that these AMOs are also quite effective in rescuing the splicing impairment of several exonic splicing variants, thus widening the potential use of these effectors for GSDII treatment.

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.

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.

Exonic CLDN16 mutations associated with familial hypomagnesemia with hypercalciuria and nephrocalcinosis can induce deleterious mRNA alterations

BACKGROUND

Familial hypomagnesaemia with hypercalciuria and nephrocalcinosis type 1 is an autosomal recessive disease characterized by excessive renal magnesium and calcium excretion, bilateral nephrocalcinosis, and progressive chronic renal failure. This rare disease is caused by mutations in CLDN16 that encodes claudin-16, a tight-junction protein involved in paracellular reabsorption of magnesium and calcium in the renal tubule. Most of these variants are located in exons and have been classified as missense mutations. The functional consequences of some of these claudin-16 mutant proteins have been analysed after heterologous expression showing indeed a significant loss of function compared to the wild-type claudin-16. We hypothesize that a number of CLDN16 exonic mutations can be responsible for the disease phenotype by disrupting the pre-mRNA splicing process.

METHODS

We selected 12 previously described presumed CLDN16 missense mutations and analysed their potential effect on pre-mRNA splicing using a minigene assay.

RESULTS

Our results indicate that five of these mutations induce significant splicing alterations. Mutations c.453G > T and c.446G > T seem to inactivate exonic splicing enhancers and promote the use of an internal cryptic acceptor splice site resulting in inclusion of a truncated exon 3 in the mature mRNA. Mutation c.571G > A affects an exonic splicing enhancer resulting in partial skipping of exon 3. Mutations c.593G > C and c.593G > A disturb the acceptor splice site of intron 3 and cause complete exon 4 skipping.

CONCLUSIONS

To our knowledge, this is the first report of CLDN16 exonic mutations producing alterations in splicing. We suggest that in the absence of patients RNA samples, splicing functional assays with minigenes could be valuable for evaluating the effect of exonic CLDN16 mutations on pre-mRNA splicing.

Networks of mRNA Processing and Alternative Splicing Regulation in Health and Disease

mRNA processing events introduce an intricate layer of complexity into gene expression processes, supporting a tremendous level of diversification of the genome's coding and regulatory potential, particularly in vertebrate species. The recent development of massive parallel sequencing methods and their adaptation to the identification and quantification of different RNA species and the dynamics of mRNA metabolism and processing has generated an unprecedented view over the regulatory networks that are established at this level, which contribute to sustain developmental, tissue specific or disease specific gene expression programs. In this chapter, we provide an overview of the recent evolution of transcriptome profiling methods and the surprising insights that have emerged in recent years regarding distinct mRNA processing events - from the 5' end to the 3' end of the molecule.

RNA splicing analysis in genomic medicine

High-throughput next-generation sequencing technologies have led to a rapid increase in the number of sequence variants identified in clinical practice via diagnostic genetic tests. Current bioinformatic analysis pipelines fail to take adequate account of the possible splicing effects of such variants, particularly where variants fall outwith canonical splice site sequences, and consequently the pathogenicity of such variants may often be missed. The regulation of splicing is highly complex and as a result, in silico prediction tools lack sufficient sensitivity and specificity for reliable use. Variants of all kinds can be linked to aberrant splicing in disease and the need for correct identification and diagnosis grows ever more crucial as novel splice-switching antisense oligonucleotide therapies start to enter clinical usage. RT-PCR provides a useful targeted assay of the splicing effects of identified variants, while minigene assays, massive parallel reporter assays and animal models can also be used for more detailed study of a particular splicing system, given enough time and resources. However, RNA-sequencing (RNA-seq) has the potential to be used as a rapid diagnostic tool in genomic medicine. By utilising data science approaches and machine learning, it may prove possible to finally understand and interpret the 'splicing code' and apply this knowledge in human disease diagnostics.

PRP4KA, a Putative Spliceosomal Protein Kinase, Is Important for Alternative Splicing and Development in Arabidopsis thaliana

Prp4 kinase (Prp4k) is the first spliceosome-associated kinase shown to regulate splicing in fungi and metazoans, but nothing is yet known about its functions in plants. Here, Kanno and Venhuizen et al. report...

Splicing of precursor messenger RNAs (pre-mRNAs) is an essential step in the expression of most eukaryotic genes. Both constitutive splicing and alternative splicing, which produces multiple messenger RNA (mRNA) isoforms from a single primary transcript, are modulated by reversible protein phosphorylation. Although the plant splicing machinery is known to be a target for phosphorylation, the protein kinases involved remain to be fully defined. We report here the identification of pre-mRNA processing 4 (PRP4) KINASE A (PRP4KA) in a forward genetic screen based on an alternatively spliced GFP reporter gene in Arabidopsis thaliana (Arabidopsis). Prp4 kinase is the first spliceosome-associated kinase shown to regulate splicing in fungi and mammals but it has not yet been studied in plants. In the same screen we identified mutants defective in SAC3A, a putative mRNA export factor that is highly coexpressed with PRP4KA in Arabidopsis. Whereas the sac3a mutants appear normal, the prp4ka mutants display a pleiotropic phenotype featuring atypical rosettes, late flowering, tall final stature, reduced branching, and lowered seed set. Analysis of RNA-sequencing data from prp4ka and sac3a mutants identified widespread and partially overlapping perturbations in alternative splicing in the two mutants. Quantitative phosphoproteomic profiling of a prp4ka mutant detected phosphorylation changes in several serine/arginine-rich proteins, which regulate constitutive and alternative splicing, and other splicing-related factors. Tests of PRP4KB, the paralog of PRP4KA, indicated that the two genes are not functionally redundant. The results demonstrate the importance of PRP4KA for alternative splicing and plant phenotype, and suggest that PRP4KA may influence alternative splicing patterns by phosphorylating a subset of splicing regulators.

Alternative Splicing of Differentiated Myeloid Cell Transcripts after Infection by Anaplasma phagocytophilum Impacts a Selective Group of Cellular Programs

Eukaryotic proteome diversity exceeds that encoded within individual genes, and results in part from alternative splicing events of pre-messenger RNA. The diversity of these splicing events can shape the outcome in development and differentiation of normal tissues, and is important in pathogenic circumstances such as cancer and some heritable conditions. A role for alternative splicing of eukaryotic genes in response to viral and intracellular bacterial infections has only recently been recognized, and plays an important role in providing fitness for microbial survival, while potentially enhancing pathogenicity. Anaplasma phagocytophilum survives within mammalian neutrophils by reshaping transcriptional programs that govern cellular functions. We applied next generation RNAseq to ATRA-differentiated HL-60 cells established to possess transcriptional and functional responses similar to A. phagocytophilum-infected human neutrophils. This demonstrated an increase in transcripts with infection and high proportion of alternatively spliced transcript events (ASEs) for which predicted gene ontology processes were in part distinct from those identified by evaluation of single transcripts or gene-level analyses alone. The alternative isoforms are not on average shorter, and no alternative splicing in genes encoding spliceosome components is noted. Although not evident at gene-level analyses, individual spliceosome transcripts that impact nearly all spliceosome components were significantly upregulated. How the distinct GO processes predicted by ASEs are regulated by infection and whether they are relevant to fitness or pathogenicity of A. phagocytophilum should be addressed in more detailed studies.

The Role of E6 Spliced Isoforms (E6*) in Human Papillomavirus-Induced Carcinogenesis

Persistent infections with High Risk Human Papillomaviruses (HR-HPVs) are the main cause of cervical cancer development. The E6 and E7 oncoproteins of HR-HPVs are derived from a polycistronic pre-mRNA transcribed from an HPV early promoter. Through alternative splicing, this pre-mRNA produces a variety of E6 spliced transcripts termed E6*. In pre-malignant lesions and HPV-related cancers, different E6/E6* transcriptional patterns have been found, although they have not been clearly associated to cancer development. Moreover, there is a controversy about the participation of E6* proteins in cancer progression. This review addresses the regulation of E6 splicing and the different functions that have been found for E6* proteins, as well as their possible role in HPV-induced carcinogenesis.