Expression of 24,426 human alternative splicing events and predicted cis regulation in 48 tissues and cell lines

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.

RBFOX1 Regulates Calcium Signaling and Enhances SERCA2 Translation

RBFOX1 is an RNA-binding protein that regulates alternative splicing and RNA processing in the neurons, skeletal muscle, and heart. We intended to define the role of RBFOX1 in regulating calcium homeostasis to maintain normal cardiac function. We generated cardiomyocyte-specific Rbfox1 gene-deletion mice (cKO). The cardiomyocyte-specific deletion of RBFOX1 was confirmed by Western blotting and immunohistochemistry. The cKO mice showed mild hypertrophy and depressed cardiac function under homeostatic conditions, which did not deteriorate with age. Pressure overload by trans-aortic constriction (TAC) caused exaggerated cardiac hypertrophy and accelerated heart failure in cKO compared with wild-type mice. We performed Western blotting to assess the expression of important Ca2+-handling proteins, which showed alterations in the phosphorylation of PLN and CAMKII and decreased expression of SERCA2. We measured the Ca2+ dynamics and noted significantly delayed Ca2+ reuptake into the sarcoplasmic reticulum. Importantly, the decrease in SERCA2 expression was not due to reduced mRNA expression or altered splicing. To assess the possibility of the post-transcriptional regulation of SERCA2 expression by RBFOX1, we performed RNA immunoprecipitation (RIP), which showed the binding of RBFOX1 protein to Serca2 mRNA, which was confirmed in luciferase assays with the Serca2a 3'-untranslated region fused to luciferase. Finally, we performed a puromycin incorporation experiment, which showed that RBFOX1 enhances SERCA2 protein translation. Our results show that RBFOX1 plays a crucial role in regulating the expression of Ca2+-handling genes to maintain normal cardiac function. We show an important post-transcriptional role of RBFOX1 in regulating SERCA2 expression.

Alternative splicing landscape in mouse skeletal muscle and adipose tissue: Effects of intermittent fasting and exercise

Alternative splicing contributes to diversify the cellular protein landscape, but aberrant splicing is implicated in many diseases. To which extent mis-splicing contributes to insulin resistance as the causal defect of type 2 diabetes and whether this can be reversed by lifestyle interventions is largely unknown. Therefore, RNA sequencing data from skeletal muscle and adipose tissue of diabetes-susceptible NZO mice treated with or without intermittent fasting and of healthy C57BL/6J mice subjected to exercise were analyzed for alternative splicing differences using Whippet and rMATS. Diet and exercise interventions triggered comparable levels of splicing changes, although the splicing profile of skeletal muscle appeared to be more flexible than that of adipose tissue, with 72-114 differential splicing events in muscle and less than 25 in adipose tissue. Splicing changes induced by time-restricted feeding, alternate-day fasting and exercise were generally mild, with a maximal percent spliced in (PSI) difference of 67%, indicating that alternative splicing plays a rather minor role in lifestyle-induced adaptations of muscle and adipose tissue in mice. However, intron retention contributed to the regulation of gene expression, influencing genes whose expression was directly linked to phenotypic parameters (e.g. Eno2 and Pan2). Alternate-day fasting promoted skipping of exon 7 in Mlxipl (coding for ChREBP), thereby affecting the glucose sensing module of this carbohydrate-responsive transcription factor. Both intermittent fasting and exercise training led to alternative splicing of known diabetes-related GWAS genes (e.g. Abcc8, Ifnar2, Smarcad1), highlighting the potential metabolic relevance of these changes.

The Spatial-Temporal Alternative Splicing Profile Reveals the Functional Diversity of FXR1 Isoforms in Myogenesis

Alternative splicing (AS) is a fundamental mechanism contributing to proteome diversity, yet its comprehensive landscape and regulatory dynamics during skeletal muscle development remain largely unexplored. Here, the temporal AS profiles are investigated during myogenesis in five vertebrates, conducting comprehensive profiling across 27 developmental stages in skeletal muscle and encompassing ten tissues in adult pigs. The analysis reveals a pervasive and evolutionarily conserved pattern of alternative exon usage throughout myogenic differentiation, with hundreds of skipped exons (SEs) showing developmental regulation, particularly within skeletal muscle. Notably, this study identifies a muscle-specific SE (exon 15) within the Fxr1 gene, whose AS generates two dynamically expressed isoforms with distinct functions: the isoform without exon 15 (Fxr1E15 -) regulates myoblasts proliferation, while the isoform incorporating exon 15 (Fxr1E15+) promotes myogenic differentiation and fusion. Transcriptome analysis suggests that specifically knocking-down Fxr1E15+ isoform in myoblasts modulates differentiation by influencing gene expression and splicing of specific targets. The increased inclusion of exon 15 during differentiation is mediated by the binding of Rbm24 to the intron. Furthermore, in vivo experiments indicate that the Fxr1E15+ isoform facilitates muscle regeneration. Collectively, these findings provide a comprehensive resource for AS studies in skeletal muscle development, underscoring the diverse functions and regulatory mechanisms governing distinct Fxr1 isoforms in myogenesis.

Long-read subcellular fractionation and sequencing reveals the translational fate of full-length mRNA isoforms during neuronal differentiation

Alternative splicing (AS) alters the cis-regulatory landscape of mRNA isoforms, leading to transcripts with distinct localization, stability, and translational efficiency. To rigorously investigate mRNA isoform-specific ribosome association, we generated subcellular fractionation and sequencing (Frac-seq) libraries using both conventional short reads and long reads from human embryonic stem cells (ESCs) and neural progenitor cells (NPCs) derived from the same ESCs. We performed de novo transcriptome assembly from high-confidence long reads from cytosolic, monosomal, light, and heavy polyribosomal fractions and quantified their abundance using short reads from their respective subcellular fractions. Thousands of transcripts in each cell type exhibited association with particular subcellular fractions relative to the cytosol. Of the multi-isoform genes, 27% and 19% exhibited significant differential isoform sedimentation in ESCs and NPCs, respectively. Alternative promoter usage and internal exon skipping accounted for the majority of differences between isoforms from the same gene. Random forest classifiers implicated coding sequence (CDS) and untranslated region (UTR) lengths as important determinants of isoform-specific sedimentation profiles, and motif analyses reveal potential cell type-specific and subcellular fraction-associated RNA-binding protein signatures. Taken together, our data demonstrate that alternative mRNA processing within the CDS and UTRs impacts the translational control of mRNA isoforms during stem cell differentiation, and highlight the utility of using a novel long-read sequencing-based method to study translational control.

Insights into the epitranscriptomic role of N6-methyladenosine on aging skeletal muscle

The modification of RNA through the N6-methyladenosine (m6A) has emerged as a growing area of research due to its regulatory role in gene expression and various biological processes regulating the expression of genes. m6A RNA methylation is a post-transcriptional modification that is dynamic and reversible and found in mRNA, tRNA, rRNA, and other non-coding RNA of most eukaryotic cells. It is executed by special proteins known as "writers," which initiate methylation; "erasers," which remove methylation; and "readers," which recognize it and regulate the expression of the gene. Modification by m6A regulates gene expression by affecting the splicing, translation, stability, and localization of mRNA. Aging causes molecular and cellular damage, which forms the basis of most age-related diseases. The decline in skeletal muscle mass and functionality because of aging leads to metabolic disorders and morbidities. The inability of aged muscles to regenerate and repair after injury poses a great challenge to the geriatric populace. This review seeks to explore the m6A epigenetic regulation in the myogenesis and regeneration processes in skeletal muscle as well as the progress made on the m6A epigenetic regulation of aging skeletal muscles.

Human-specific paralogs of SRGAP2 induce neotenic features of microglia structural and functional maturation

Microglia play key roles in shaping synaptic connectivity during neural circuits development. Whether microglia display human-specific features of structural and functional maturation is currently unknown. We show that the ancestral gene SRGAP2A and its human-specific (HS) paralogs SRGAP2B/C are not only expressed in cortical neurons but are the only HS gene duplications expressed in human microglia. Here, using combination of xenotransplantation of human induced pluripotent stem cell (hiPSC)-derived microglia and mouse genetic models, we demonstrate that (1) HS SRGAP2B/C are necessary and sufficient to induce neotenic features of microglia structural and functional maturation in a cell-autonomous manner, and (2) induction of SRGAP2-dependent neotenic features of microglia maturation non-cell autonomously impacts synaptic development in cortical pyramidal neurons. Our results reveal that, during human brain evolution, human-specific genes SRGAP2B/C coordinated the emergence of neotenic features of synaptic development by acting as genetic modifiers of both neurons and microglia.

Reverse engineering neuron type-specific and type-orthogonal splicing-regulatory networks using single-cell transcriptomes

Cell type-specific alternative splicing (AS) enables differential gene isoform expression between diverse neuron types with distinct identities and functions. Current studies linking individual RNA-binding proteins (RBPs) to AS in a few neuron types underscore the need for holistic modeling. Here, we use network reverse engineering to derive a map of the neuron type-specific AS regulatory landscape from 133 mouse neocortical cell types defined by single-cell transcriptomes. This approach reliably inferred the regulons of 350 RBPs and their cell type-specific activities. Our analysis revealed driving factors delineating neuronal identities, among which we validated Elavl2 as a key RBP for MGE-specific splicing in GABAergic interneurons using an in vitro ESC differentiation system. We also identified a module of exons and candidate regulators specific for long- and short-projection neurons across multiple neuronal classes. This study provides a resource for elucidating splicing regulatory programs that drive neuronal molecular diversity, including those that do not align with gene expression-based classifications.

Long-read sequencing for 29 immune cell subsets reveals disease-linked isoforms

Alternative splicing events are a major causal mechanism for complex traits, but they have been understudied due to the limitation of short-read sequencing. Here, we generate a full-length isoform annotation of human immune cells from an individual by long-read sequencing for 29 cell subsets. This contains a number of unannotated transcripts and isoforms such as a read-through transcript of TOMM40-APOE in the Alzheimer's disease locus. We profile characteristics of isoforms and show that repetitive elements significantly explain the diversity of unannotated isoforms, providing insight into the human genome evolution. In addition, some of the isoforms are expressed in a cell-type specific manner, whose alternative 3'-UTRs usage contributes to their specificity. Further, we identify disease-associated isoforms by isoform switch analysis and by integration of several quantitative trait loci analyses with genome-wide association study data. Our findings will promote the elucidation of the mechanism of complex diseases via alternative splicing.

RNA Panel Sequencing Is an Effective Tool to Help Classify Splice Variants for Clinical Oncogenetic Diagnosis

Routine gene panel analysis identifies pathogenic variants in clinically relevant genes. However, variants of unknown significance (VUSs) are commonly observed, many of which potentially have an impact on mRNA transcription and splicing. Several software programs attempt to predict the impact of variants on splicing and thus make it possible to select the variants for which it is important to study the effect on the transcripts. Transcript analysis is also necessary to show the tandem character of large duplications, and it can be useful for the search for deep intronic variants that are difficult to identify in a DNA panel. We analyzed 53 variants of unknown significance by targeted sequencing of 48 genes using RNA extracted from patient blood samples. RT-PCR and Sanger sequencing of patient mRNA or minigene monoallelic analysis was also carried out when necessary. For the 53 VUSs, 21 could be classified as likely neutral and 10 as pathogenic or likely pathogenic. Data are comprehensively presented for four variants: PTEN c.206+6T>G, MLH1 c.791-489_791-20del, BRCA2 c.68-8_68-7delinsAA, and MSH2 c.(1076+1_1077-1)_(1276+1_1277-1)dup. These four examples illustrate the usefulness of blood RNA panel sequencing in clinical oncogenetics to help classify VUSs with predicted splice effects. It could also be useful for characterizing large duplications and for detecting deep intronic variants with an impact on expressed transcripts.

scaRNA20 promotes pseudouridylatory modification of small nuclear snRNA U12 and improves cardiomyogenesis

Non-coding RNAs, particularly small Cajal-body associated RNAs (scaRNAs), play a significant role in spliceosomal RNA modifications. While their involvement in ischemic myocardium regeneration is known, their role in cardiac development is unexplored. We investigated scaRNA20's role in iPSC differentiation into cardiomyocytes (iCMCs) via overexpression and knockdown assays. We measured scaRNA20-OE-iCMCs and scaRNA20-KD-iCMCs contractility using Particle Image Velocimetry (PIV), comparing them to control iCMCs. We explored scaRNA20's impact on alternative splicing via pseudouridylation (Ψ) of snRNA U12, analyzing its functional consequences in cardiac differentiation. scaRNA20-OE-iPSC differentiation increased beating colonies, upregulated cardiac-specific genes, activated TP53 and STAT3, and preserved contractility under hypoxia. Conversely, scaRNA20-KD-iCMCs exhibited poor differentiation and contractility. STAT3 inhibition in scaRNA20-OE-iPSCs hindered cardiac differentiation. RNA immunoprecipitation revealed increased Ψ at the 28th uridine of U12 RNA in scaRNA20-OE iCMCs. U12-KD iCMCs had reduced cardiac differentiation, which improved upon U12 RNA introduction. In summary, scaRNA20-OE in iPSCs enhances cardiomyogenesis, preserves iCMC function under hypoxia, and may have implications for ischemic myocardium regeneration.

STIGMA: Single-cell tissue-specific gene prioritization using machine learning

Clinical exome and genome sequencing have revolutionized the understanding of human disease genetics. Yet many genes remain functionally uncharacterized, complicating the establishment of causal disease links for genetic variants. While several scoring methods have been devised to prioritize these candidate genes, these methods fall short of capturing the expression heterogeneity across cell subpopulations within tissues. Here, we introduce single-cell tissue-specific gene prioritization using machine learning (STIGMA), an approach that leverages single-cell RNA-seq (scRNA-seq) data to prioritize candidate genes associated with rare congenital diseases. STIGMA prioritizes genes by learning the temporal dynamics of gene expression across cell types during healthy organogenesis. To assess the efficacy of our framework, we applied STIGMA to mouse limb and human fetal heart scRNA-seq datasets. In a cohort of individuals with congenital limb malformation, STIGMA prioritized 469 variants in 345 genes, with UBA2 as a notable example. For congenital heart defects, we detected 34 genes harboring nonsynonymous de novo variants (nsDNVs) in two or more individuals from a set of 7,958 individuals, including the ortholog of Prdm1, which is associated with hypoplastic left ventricle and hypoplastic aortic arch. Overall, our findings demonstrate that STIGMA effectively prioritizes tissue-specific candidate genes by utilizing single-cell transcriptome data. The ability to capture the heterogeneity of gene expression across cell populations makes STIGMA a powerful tool for the discovery of disease-associated genes and facilitates the identification of causal variants underlying human genetic disorders.

Comparative analyses of dynamic transcriptome profiles highlight key response genes and dominant isoforms for muscle development and growth in chicken

BACKGROUND

Modern breeding strategies have resulted in significant differences in muscle mass between indigenous chicken and specialized broiler. However, the molecular regulatory mechanisms that underlie these differences remain elusive. The aim of this study was to identify key genes and regulatory mechanisms underlying differences in breast muscle development between indigenous chicken and specialized broiler.

RESULTS

Two time-series RNA-sequencing profiles of breast muscles were generated from commercial Arbor Acres (AA) broiler (fast-growing) and Chinese indigenous Lushi blue-shelled-egg (LS) chicken (slow-growing) at embryonic days 10, 14, and 18, and post-hatching day 1 and weeks 1, 3, and 5. Principal component analysis of the transcriptome profiles showed that the top four principal components accounted for more than 80% of the total variance in each breed. The developmental axes between the AA and LS chicken overlapped at the embryonic stages but gradually separated at the adult stages. Integrative investigation of differentially-expressed transcripts contained in the top four principal components identified 44 genes that formed a molecular network associated with differences in breast muscle mass between the two breeds. In addition, alternative splicing analysis revealed that genes with multiple isoforms always had one dominant transcript that exhibited a significantly higher expression level than the others. Among the 44 genes, the TNFRSF6B gene, a mediator of signal transduction pathways and cell proliferation, harbored two alternative splicing isoforms, TNFRSF6B-X1 and TNFRSF6B-X2. TNFRSF6B-X1 was the dominant isoform in both breeds before the age of one week. A switching event of the dominant isoform occurred at one week of age, resulting in TNFRSF6B-X2 being the dominant isoform in AA broiler, whereas TNFRSF6B-X1 remained the dominant isoform in LS chicken. Gain-of-function assays demonstrated that both isoforms promoted the proliferation of chicken primary myoblasts, but only TNFRSF6B-X2 augmented the differentiation and intracellular protein content of chicken primary myoblasts.

CONCLUSIONS

For the first time, we identified several key genes and dominant isoforms that may be responsible for differences in muscle mass between slow-growing indigenous chicken and fast-growing commercial broiler. These findings provide new insights into the regulatory mechanisms underlying breast muscle development in chicken.

Cell-type specific regulator RBPMS switches alternative splicing via higher-order oligomerization and heterotypic interactions with other splicing regulators

Alternative pre-mRNA splicing decisions are regulated by RNA binding proteins (RBPs) that can activate or repress regulated splice sites. Repressive RBPs typically harness multivalent interactions to bind stably to target RNAs. Multivalency can be achieved by homomeric oligomerization and heteromeric interactions with other RBPs, often mediated by intrinsically disordered regions (IDRs), and by possessing multiple RNA binding domains. Cell-specific splicing decisions often involve the action of widely expressed RBPs, which are able to bind multivalently around target exons, but without effect in the absence of a cell-specific regulator. To address how cell-specific regulators can collaborate with constitutive RBPs in alternative splicing regulation, we used the smooth-muscle specific regulator RBPMS. Recombinant RBPMS is sufficient to confer smooth muscle cell specific alternative splicing of Tpm1 exon 3 in cell-free assays by preventing assembly of ATP-dependent splicing complexes. This activity depends upon a C-terminal IDR that facilitates dynamic higher-order self-assembly, cooperative binding to multivalent RNA and interactions with widely expressed splicing co-regulators, including MBNL1 and RBFOX2, allowing cooperative assembly of stable cell-specific regulatory complexes.

To RNA-binding and beyond: Emerging facets of the role of Rbfox proteins in development and disease

The RNA-binding Fox-1 homologue (Rbfox) proteins represent an ancient family of splicing factors, conserved through evolution. All members share an RNA recognition motif (RRM), and a particular affinity for the GCAUG signature in target RNA molecules. The role of Rbfox, as a splice factor, deciding the tissue-specific inclusion/exclusion of an exon, depending on its binding position on the flanking introns, is well known. Rbfox often acts in concert with other splicing factors, and forms splicing regulatory networks. Apart from this canonical role, recent studies show that Rbfox can also function as a transcription co-factor, and affects mRNA stability and translation. The repertoire of Rbfox targets is vast, including genes involved in the development of tissue lineages, such as neurogenesis, myogenesis, and erythropoeiesis, and molecular processes, including cytoskeletal dynamics, and calcium handling. A second layer of complexity is added by the fact that Rbfox expression itself is regulated by multiple mechanisms, and, in vertebrates, exhibits tissue-specific expression. The optimum dosage of Rbfox is critical, and its misexpression is etiological to various disease conditions. In this review, we discuss the contextual roles played by Rbfox as a tissue-specific regulator for the expression of many important genes with diverse functions, through the lens of the emerging data which highlights its involvement in many human diseases. Furthermore, we explore the mechanistic details provided by studies in model organisms, with emphasis on the work with Drosophila. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Turnover and Surveillance > Regulation of RNA Stability RNA Processing > Splicing Regulation/Alternative Splicing.

The regulation of expression and splicing of transcription factors are related to the muscle damage caused by radiation in tree shrews

Radiotherapy-induced muscle injury (RIMI) is a major complication of radiotherapy for nasopharyngeal carcinoma. Transcription factor (TF) expression and alternative splicing are crucial events in transcriptional and posttranscriptional regulation, respectively, and are known to be involved in key signaling pathways contributing to a variety of human disorders, including radiation injury. To investigate the TFs and alternative splicing events involved in RIMI, we constructed a tree shrew model as described previously in which the RIMI group received 20 Gy of irradiation on the tensor veli palatini (TVP) muscles. The irradiated muscles were evaluated by RNA sequencing (RNA-seq) 6 months later, and the results compared with those for normal TVP muscles. The alt5p and alt3p events were the two main types of differentially regulated alternative splicing events (RASEs) identified via the Splice sites Usage Variation Analysis (SUVA) software, and these RASEs were highly conserved in RIMI. According to functional enrichment analysis, the differentially RASEs were primarily enriched in pathways related to transcriptional regulation. Furthermore, we identified 16 alternative splicing TFs (ASTFs) in ASTF-differentially expressed gene (DEG) networks based on co-expression analysis, and the regulatory networks were chiefly enriched in pathways linked to cell proliferation and differentiation. This study revealed that RASEs and ASTF-DEG networks may both play important regulatory roles in gene expression network alteration in RIMI. Future studies on the targeting mechanisms and early interventions directed at RASEs and ASTF-DEG networks may aid in the treatment of RIMI.

Multi-omics data analysis reveals the biological implications of alternative splicing events in lung adenocarcinoma

Cancer is characterized by the dysregulation of alternative splicing (AS). However, the comprehensive regulatory mechanisms of AS in lung adenocarcinoma (LUAD) are poorly understood. Here, we displayed the AS landscape in LUAD based on the integrated analyses of LUAD's multi-omics data. We identified 13,995 AS events in 6309 genes as differentially expressed alternative splicing events (DEASEs) mainly covering protein-coding genes. These DEASEs were strongly linked to "cancer hallmarks", such as apoptosis, DNA repair, cell cycle, cell proliferation, angiogenesis, immune response, generation of precursor metabolites and energy, p53 signaling pathway and PI3K-AKT signaling pathway. We further built a regulatory network connecting splicing factors (SFs) and DEASEs. In addition, RNA-binding protein (RBP) mutations that can affect DEASEs were investigated to find some potential cancer drivers. Further association analysis demonstrated that DNA methylation levels were highly correlated with DEASEs. In summary, our results can bring new insight into understanding the mechanism of AS and provide novel biomarkers for personalized medicine of LUAD.

Using the ACMG/AMP framework to capture evidence related to predicted and observed impact on splicing: Recommendations from the ClinGen SVI Splicing Subgroup

The American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology (AMP) framework for classifying variants uses six evidence categories related to the splicing potential of variants: PVS1, PS3, PP3, BS3, BP4, and BP7. However, the lack of guidance on how to apply such codes has contributed to variation in the specifications developed by different Clinical Genome Resource (ClinGen) Variant Curation Expert Panels. The ClinGen Sequence Variant Interpretation Splicing Subgroup was established to refine recommendations for applying ACMG/AMP codes relating to splicing data and computational predictions. We utilized empirically derived splicing evidence to (1) determine the evidence weighting of splicing-related data and appropriate criteria code selection for general use, (2) outline a process for integrating splicing-related considerations when developing a gene-specific PVS1 decision tree, and (3) exemplify methodology to calibrate splice prediction tools. We propose repurposing the PVS1_Strength code to capture splicing assay data that provide experimental evidence for variants resulting in RNA transcript(s) with loss of function. Conversely, BP7 may be used to capture RNA results demonstrating no splicing impact for intronic and synonymous variants. We propose that the PS3/BS3 codes are applied only for well-established assays that measure functional impact not directly captured by RNA-splicing assays. We recommend the application of PS1 based on similarity of predicted RNA-splicing effects for a variant under assessment in comparison with a known pathogenic variant. The recommendations and approaches for consideration and evaluation of RNA-assay evidence described aim to help standardize variant pathogenicity classification processes when interpreting splicing-based evidence.

APPLICATION OF THE ACMG/AMP FRAMEWORK TO CAPTURE EVIDENCE RELEVANT TO PREDICTED AND OBSERVED IMPACT ON SPLICING: RECOMMENDATIONS FROM THE CLINGEN SVI SPLICING SUBGROUP

The American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) framework for classifying variants uses six evidence categories related to the splicing potential of variants: PVS1 (null variant in a gene where loss-of-function is the mechanism of disease), PS3 (functional assays show damaging effect on splicing), PP3 (computational evidence supports a splicing effect), BS3 (functional assays show no damaging effect on splicing), BP4 (computational evidence suggests no splicing impact), and BP7 (silent change with no predicted impact on splicing). However, the lack of guidance on how to apply such codes has contributed to variation in the specifications developed by different Clinical Genome Resource (ClinGen) Variant Curation Expert Panels. The ClinGen Sequence Variant Interpretation (SVI) Splicing Subgroup was established to refine recommendations for applying ACMG/AMP codes relating to splicing data and computational predictions. Our study utilised empirically derived splicing evidence to: 1) determine the evidence weighting of splicing-related data and appropriate criteria code selection for general use, 2) outline a process for integrating splicing-related considerations when developing a gene-specific PVS1 decision tree, and 3) exemplify methodology to calibrate bioinformatic splice prediction tools. We propose repurposing of the PVS1_Strength code to capture splicing assay data that provide experimental evidence for variants resulting in RNA transcript(s) with loss of function. Conversely BP7 may be used to capture RNA results demonstrating no impact on splicing for both intronic and synonymous variants, and for missense variants if protein functional impact has been excluded. Furthermore, we propose that the PS3 and BS3 codes are applied only for well-established assays that measure functional impact that is not directly captured by RNA splicing assays. We recommend the application of PS1 based on similarity of predicted RNA splicing effects for a variant under assessment in comparison to a known Pathogenic variant. The recommendations and approaches for consideration and evaluation of RNA assay evidence described aim to help standardise variant pathogenicity classification processes and result in greater consistency when interpreting splicing-based evidence.

Re-evaluating the impact of alternative RNA splicing on proteomic diversity

Alternative splicing (AS) constitutes a mechanism by which protein-coding genes and long non-coding RNA (lncRNA) genes produce more than a single mature transcript. From plants to humans, AS is a powerful process that increases transcriptome complexity. Importantly, splice variants produced from AS can potentially encode for distinct protein isoforms which can lose or gain specific domains and, hence, differ in their functional properties. Advances in proteomics have shown that the proteome is indeed diverse due to the presence of numerous protein isoforms. For the past decades, with the help of advanced high-throughput technologies, numerous alternatively spliced transcripts have been identified. However, the low detection rate of protein isoforms in proteomic studies raised debatable questions on whether AS contributes to proteomic diversity and on how many AS events are really functional. We propose here to assess and discuss the impact of AS on proteomic complexity in the light of the technological progress, updated genome annotation, and current scientific knowledge.

Comprehensive analysis of pre-mRNA alternative splicing regulated by m6A methylation in pig oxidative and glycolytic skeletal muscles

BACKGROUND

Different types of skeletal myofibers exhibit distinct physiological and metabolic properties that are associated with meat quality traits in livestock. Alternative splicing (AS) of pre-mRNA can generate multiple transcripts from an individual gene by differential selection of splice sites. N6-methyladenosine (m6A) is the most abundant modification in mRNAs, but its regulation for AS in different muscles remains unknown.  RESULTS: We characterized AS events and m6A methylation pattern in pig oxidative and glycolytic muscles. A tota1 of 1294 differential AS events were identified, and differentially spliced genes were significantly enriched in processes related to different phenotypes between oxidative and glycolytic muscles. We constructed the regulatory network between splicing factors and corresponding differential AS events and identified NOVA1 and KHDRBS2 as key splicing factors. AS event was enriched in m6A-modified genes, and the methylation level was positively correlated with the number of AS events in genes. The dynamic change in m6A enrichment was associated with 115 differentially skipping exon (SE-DAS) events within 92 genes involving in various processes, including muscle contraction and myofibril assembly. We obtained 23.4% SE-DAS events (27/115) regulated by METTL3-meditaed m6A and experimentally validated the aberrant splicing of ZNF280D, PHE4DIP, and NEB. The inhibition of m6A methyltransferase METTL3 could induce the conversion of oxidative fiber to glycolytic fiber in PSCs.

CONCLUSION

Our study suggested that m6A modification could contribute to significant difference in phenotypes between oxidative and glycolytic muscles by mediating the regulation of AS. These findings would provide novel insights into mechanisms underlying muscle fiber conversion.

The Role of Alternative Splicing Factors hnRNP G and Fox-2 in the Progression and Prognosis of Esophageal Cancer

Aim. Alternative splicing (AS) has been widely demonstrated in the occurrence and progression of many cancers. Nevertheless, the involvement of cancer-associated splicing factors in the development of esophageal carcinoma (ESCA) remains to be explored. Method. RNA-Seq data and the corresponding clinical information of the ESCA cohort were downloaded from The Cancer Genome Atlas database. Bioinformatics methods were used to further analyzed the differently expressed AS (DEAS) events and their splicing network. Kaplan–Meier, Cox regression, and unsupervised cluster analyses were used to assess the association between AS events and clinical characteristics of ESCA patients. The splicing factors screened out were verified in vitro at the cellular level. Results. A total of 50,342 AS events were identified, of which 3,988 were DEAS events and 46 of these were associated with overall survival (OS) of ESCA patients, with a 5-year OS rate of 0.941. By constructing a network of AS events with survival-related splicing factors, the AS factors related to prognosis can be further identified. In vitro experiments and database analysis confirmed that the high expression of hnRNP G in ESCA is related to the high invasion ability of ESCA cells and the poor prognosis of ESCA patients. In contrast, the low expression of fox-2 in esophageal cancer is related to a better prognosis. Conclusion. ESCA-associated AS factors hnRNP G and Fox-2 are of great value in deciphering the underlying mechanisms of AS in ESCA and providing clues for therapeutic goals for further validation.

Multi-omics analysis reveals RNA splicing alterations and their biological and clinical implications in lung adenocarcinoma

Alternative RNA splicing is one of the most important mechanisms of posttranscriptional gene regulation, which contributes to protein diversity in eukaryotes. It is well known that RNA splicing dysregulation is a critical mechanism in tumor pathogenesis and the rationale for the promising splice-switching therapeutics for cancer treatment. Although we have a comprehensive understanding of DNA mutations, abnormal gene expression profiles, epigenomics, and proteomics in lung adenocarcinoma (LUAD), little is known about its aberrant alternative splicing profiles. Here, based on the multi-omics data generated from over 1000 samples, we systematically studied the RNA splicing alterations in LUAD and revealed their biological and clinical implications. We identified 3688 aberrant alternative splicing events (AASEs) in LUAD, most of which were alternative promoter and exon skip. The specific regulatory roles of RNA binding proteins, somatic mutations, and DNA methylations on AASEs were comprehensively interrogated. We dissected the functional implications of AASEs and concluded that AASEs mainly affected biological processes related to tumor proliferation and metastasis. We also found that one subtype of LUAD with a particular AASEs pattern was immunogenic and had a better prognosis and response rate to immunotherapy. These findings revealed novel events related to tumorigenesis and tumor immune microenvironment and laid the foundation for the development of splice-switching therapies for LUAD.

Maternal Microbiota Modulate a Fragile X-like Syndrome in Offspring Mice

Maternal microbial dysbiosis has been implicated in adverse postnatal health conditions in offspring, such as obesity, cancer, and neurological disorders. We observed that the progeny of mice fed a Westernized diet (WD) with low fiber and extra fat exhibited higher frequencies of stereotypy, hyperactivity, cranial features and lower FMRP protein expression, similar to what is typically observed in Fragile X Syndrome (FXS) in humans. We hypothesized that gut dysbiosis and inflammation during pregnancy influenced the prenatal uterine environment, leading to abnormal phenotypes in offspring. We found that oral in utero supplementation with a beneficial anti-inflammatory probiotic microbe, Lactobacillus reuteri, was sufficient to inhibit FXS-like phenotypes in offspring mice. Cytokine profiles in the pregnant WD females showed that their circulating levels of pro-inflammatory cytokine interleukin (Il)-17 were increased relative to matched gravid mice and to those given supplementary L. reuteri probiotic. To test our hypothesis of prenatal contributions to this neurodevelopmental phenotype, we performed Caesarian (C-section) births using dissimilar foster mothers to eliminate effects of maternal microbiota transferred during vaginal delivery or nursing after birth. We found that foster-reared offspring still displayed a high frequency of these FXS-like features, indicating significant in utero contributions. In contrast, matched foster-reared progeny of L. reuteri-treated mothers did not exhibit the FXS-like typical features, supporting a key role for microbiota during pregnancy. Our findings suggest that diet-induced dysbiosis in the prenatal uterine environment is strongly associated with the incidence of this neurological phenotype in progeny but can be alleviated by addressing gut dysbiosis through probiotic supplementation.

The Splicing of the Mitochondrial Calcium Uniporter Genuine Activator MICU1 Is Driven by RBFOX2 Splicing Factor during Myogenic Differentiation

Alternative splicing, the process by which exons within a pre-mRNA transcript are differentially joined or skipped, is crucial in skeletal muscle since it is required both during myogenesis and in post-natal life to reprogram the transcripts of contractile proteins, metabolic enzymes, and transcription factors in functionally distinct muscle fiber types. The importance of such events is underlined by the numerosity of pathological conditions caused by alternative splicing aberrations. Importantly, many skeletal muscle Ca²⁺ homeostasis genes are also regulated by alternative splicing mechanisms, among which is the Mitochondrial Ca²⁺ Uniporter (MCU) genuine activator MICU1 which regulates MCU opening upon cell stimulation. We have previously shown that murine skeletal muscle MICU1 is subjected to alternative splicing, thereby generating a splice variant-which was named MICU1.1-that confers unique properties to the mitochondrial Ca²⁺ uptake and ensuring sufficient ATP production for muscle contraction. Here we extended the analysis of MICU1 alternative splicing to human tissues, finding two additional splicing variants that were characterized by their ability to regulate mitochondrial Ca²⁺ uptake. Furthermore, we found that MICU1 alternative splicing is induced during myogenesis by the splicing factor RBFOX2. These results highlight the complexity of the alternative splicing mechanisms in skeletal muscle and the regulation of mitochondrial Ca²⁺ among tissues.

Rbfox1 is required for myofibril development and maintaining fiber type-specific isoform expression in Drosophila muscles

Protein isoform transitions confer muscle fibers with distinct properties and are regulated by differential transcription and alternative splicing. RNA-binding Fox protein 1 (Rbfox1) can affect both transcript levels and splicing, and is known to contribute to normal muscle development and physiology in vertebrates, although the detailed mechanisms remain obscure. In this study, we report that Rbfox1 contributes to the generation of adult muscle diversity in Drosophila Rbfox1 is differentially expressed among muscle fiber types, and RNAi knockdown causes a hypercontraction phenotype that leads to behavioral and eclosion defects. Misregulation of fiber type-specific gene and splice isoform expression, notably loss of an indirect flight muscle-specific isoform of Troponin-I that is critical for regulating myosin activity, leads to structural defects. We further show that Rbfox1 directly binds the 3'-UTR of target transcripts, regulates the expression level of myogenic transcription factors myocyte enhancer factor 2 and Salm, and both modulates expression of and genetically interacts with the CELF family RNA-binding protein Bruno1 (Bru1). Rbfox1 and Bru1 co-regulate fiber type-specific alternative splicing of structural genes, indicating that regulatory interactions between FOX and CELF family RNA-binding proteins are conserved in fly muscle. Rbfox1 thus affects muscle development by regulating fiber type-specific splicing and expression dynamics of identity genes and structural proteins.

Computational Analysis of Alternative Splicing Using VAST-TOOLS and the VastDB Framework

Alternative splicing (AS) can vastly expand animal transcriptomes and proteomes. Two main open questions in the field are how AS is regulated across cell/tissue types and disease, and what roles different AS events play. To facilitate AS research, we have created the computational VastDB framework, which comprises a series of complementary software and resources that we describe in this chapter. The VastDB framework is especially designed to aid biomedical researchers without a strong computational background. It offers tools and resources to: (a) quantify AS and identify differentially spliced AS events using RNA-seq data (vast-tools), (b) perform multiple genomic and sequence analyses for investigating AS events (Matt), (c) identify AS events with genomic and regulatory conservation among species (ExOrthist), and (d) help with the biological interpretation of the results, and, ultimately, with the identification of interesting AS events to design wet-lab experiments (VastDB and PastDB).

RNA self-splicing by engineered hairpin ribozyme variants

Small RNAs capable of self-cleavage and ligation might have been the precursors for the much more complex self-splicing group I and II introns in an early RNA world. Here, we demonstrate the activity of engineered hairpin ribozyme variants, which as self-splicing introns are removed from their parent RNA. In the process, two cleavage reactions are supported at the two intron-exon junctions, followed by ligation of the two generated exon fragments. As a result, the hairpin ribozyme, here acting as the self-splicing intron, is cut out. Two self-splicing hairpin ribozyme variants were investigated, one designed by hand, the other by a computer-aided approach. Both variants perform self-splicing, generating a cut-out intron and ligated exons.

Complexity and graded regulation of neuronal cell-type-specific alternative splicing revealed by single-cell RNA sequencing

The enormous cellular diversity in the mammalian brain, which is highly prototypical and organized in a hierarchical manner, is dictated by cell-type-specific gene-regulatory programs at the molecular level. Although prevalent in the brain, the contribution of alternative splicing (AS) to the molecular diversity across neuronal cell types is just starting to emerge. Here, we systematically investigated AS regulation across over 100 transcriptomically defined neuronal types of the adult mouse cortex using deep single-cell RNA-sequencing data. We found distinct splicing programs between glutamatergic and GABAergic neurons and between subclasses within each neuronal class. These programs consist of overlapping sets of alternative exons showing differential splicing at multiple hierarchical levels. Using an integrative approach, our analysis suggests that RNA-binding proteins (RBPs) Celf1/2, Mbnl2, and Khdrbs3 are preferentially expressed and more active in glutamatergic neurons, while Elavl2 and Qk are preferentially expressed and more active in GABAergic neurons. Importantly, these and additional RBPs also contribute to differential splicing between neuronal subclasses at multiple hierarchical levels, and some RBPs contribute to splicing dynamics that do not conform to the hierarchical structure defined by the transcriptional profiles. Thus, our results suggest graded regulation of AS across neuronal cell types, which may provide a molecular mechanism to specify neuronal identity and function that are orthogonal to established classifications based on transcriptional regulation.

Transcriptome Analysis in a Primary Human Muscle Cell Differentiation Model for Myotonic Dystrophy Type 1

Myotonic dystrophy type 1 (DM1) is caused by CTG-repeat expansions leading to a complex pathology with a multisystemic phenotype that primarily affects the muscles and brain. Despite a multitude of information, especially on the alternative splicing of several genes involved in the pathology, information about additional factors contributing to the disease development is still lacking. We performed RNAseq and gene expression analyses on proliferating primary human myoblasts and differentiated myotubes. GO-term analysis indicates that in myoblasts and myotubes, different molecular pathologies are involved in the development of the muscular phenotype. Gene set enrichment for splicing reveals the likelihood of whole, differentiation stage specific, splicing complexes that are misregulated in DM1. These data add complexity to the alternative splicing phenotype and we predict that it will be of high importance for therapeutic interventions to target not only mature muscle, but also satellite cells.

Global alternative splicing landscape of skeletal muscle atrophy induced by hindlimb unloading

Background

Long-term exposure to microgravity will cause skeletal muscle atrophy, which can cause serious harm to astronauts in space travel. Therefore, it is important to explore skeletal muscle atrophy’s molecular mechanism for its prevention and treatment. However, as an important regulatory approach of skeletal muscle physiology, the role of alternative splicing in skeletal muscle atrophy, especially skeletal muscle atrophy caused by disuse, is unclear.

Methods

We established a rat hindlimb unloading model and performed RNA sequencing on soleus muscle, which was seriously atrophied during unloading. Several bioinformatics methods were used to identify alternative splicing events and determine their gene functions.

Results

Many alternative splicing events were found to occur at different time points (12 h, 24 h, 36 h, 3 days, and 7 days) after hindlimb unloading. These differential alternative splicing events mainly occurred in the gene's coding domain sequence region, and 59% of the alternative splicing events caused open reading frameshift. Bioinformatics analysis results showed that genes with different alternative splicing events were enriched in multiple pathways related to muscle atrophy, including the insulin signaling pathway, endocytosis, mitophagy, and ubiquitin-proteasome pathway. Moreover, alternative splicing of several deubiquitinase genes persisted during skeletal muscle atrophy induced by unloading. Additionally, we identified 10 differentially expressed RNA binding proteins during skeletal muscle atrophy induced by unloading, mainly containing Xpo4, Eif4e2, P4ha1, Lrrfip1, Zc3h14, Emg1, Hnrnp h1, Mbnl2, RBfox1, and Mbnl1. Hnrnp h1 and Mbnl2 were significantly downregulated, and RBfox1 and Mbnl1 were significantly upregulated during skeletal muscle atrophy caused by unloading.

Conclusions

To the best of our knowledge, the present study is the first to propose alternative splicing alterations related to disuse-induced muscle atrophy, emphasizing that alternative splicing is a new focus of attention in the occurrence of muscle atrophy.

MTSplice predicts effects of genetic variants on tissue-specific splicing

We develop the free and open-source model Multi-tissue Splicing (MTSplice) to predict the effects of genetic variants on splicing of cassette exons in 56 human tissues. MTSplice combines MMSplice, which models constitutive regulatory sequences, with a new neural network that models tissue-specific regulatory sequences. MTSplice outperforms MMSplice on predicting tissue-specific variations associated with genetic variants in most tissues of the GTEx dataset, with largest improvements on brain tissues. Furthermore, MTSplice predicts that autism-associated de novo mutations are enriched for variants affecting splicing specifically in the brain. We foresee that MTSplice will aid interpreting variants associated with tissue-specific disorders.

hnRNP L is essential for myogenic differentiation and modulates myotonic dystrophy pathologies

INTRODUCTION

RNA-binding proteins (RBPs) play an important role in skeletal muscle development and disease by regulating RNA splicing. In myotonic dystrophy type 1 (DM1), the RBP MBNL1 (muscleblind-like) is sequestered by toxic CUG repeats, leading to missplicing of MBNL1 targets. Mounting evidence from the literature has implicated other factors in the pathogenesis of DM1. Herein we sought to evaluate the functional role of the splicing factor hnRNP L in normal and DM1 muscle cells.

METHODS

Co-immunoprecipitation assays using hnRNPL and MBNL1 expression constructs and splicing profiling in normal and DM1 muscle cell lines were performed. Zebrafish morpholinos targeting hnrpl and hnrnpl2 were injected into one-cell zebrafish for developmental and muscle analysis. In human myoblasts downregulation of hnRNP L was achieved with shRNAi. Ascochlorin administration to DM1 myoblasts was performed and expression of the CUG repeats, DM1 splicing biomarkers, and hnRNP L expression levels were evaluated.

RESULTS

Using DM1 patient myoblast cell lines we observed the formation of abnormal hnRNP L nuclear foci within and outside the expanded CUG repeats, suggesting a role for this factor in DM1 pathology. We showed that the antiviral and antitumorigenic isoprenoid compound ascochlorin increased MBNL1 and hnRNP L expression levels. Drug treatment of DM1 muscle cells with ascochlorin partially rescued missplicing of established early biomarkers of DM1 and improved the defective myotube formation displayed by DM1 muscle cells.

DISCUSSION

Together, these studies revealed that hnRNP L can modulate DM1 pathologies and is a potential therapeutic target.

Restriction of an intron size en route to endothermy

Ca2+-insensitive and -sensitive E1 subunits of the 2-oxoglutarate dehydrogenase complex (OGDHC) regulate tissue-specific NADH and ATP supply by mutually exclusive OGDH exons 4a and 4b. Here we show that their splicing is enforced by distant lariat branch points (dBPs) located near the 5' splice site of the intervening intron. dBPs restrict the intron length and prevent transposon insertions, which can introduce or eliminate dBP competitors. The size restriction was imposed by a single dominant dBP in anamniotes that expanded into a conserved constellation of four dBP adenines in amniotes. The amniote clusters exhibit taxon-specific usage of individual dBPs, reflecting accessibility of their extended motifs within a stable RNA hairpin rather than U2 snRNA:dBP base-pairing. The dBP expansion took place in early terrestrial species and was followed by a uridine enrichment of large downstream polypyrimidine tracts in mammals. The dBP-protected megatracts permit reciprocal regulation of exon 4a and 4b by uridine-binding proteins, including TIA-1/TIAR and PUF60, which promote U1 and U2 snRNP recruitment to the 5' splice site and BP, respectively, but do not significantly alter the relative dBP usage. We further show that codons for residues critically contributing to protein binding sites for Ca2+ and other divalent metals confer the exon inclusion order that mirrors the Irving-Williams affinity series, linking the evolution of auxiliary splicing motifs in exons to metallome constraints. Finally, we hypothesize that the dBP-driven selection for Ca2+-dependent ATP provision by E1 facilitated evolution of endothermy by optimizing the aerobic scope in target tissues.

Post-Transcriptional Regulation in Skeletal Muscle Development, Repair, and Disease

Skeletal muscle formation is a complex process that requires tight spatiotemporal control of key myogenic factors. Emerging evidence suggests that RNA processing is crucial for the regulation of these factors, and that multiple post-transcriptional regulatory pathways work dependently and independently of one another to enable precise control of transcripts throughout muscle development and repair. Moreover, disruption of these pathways is implicated in neuromuscular disease, and the recent development of RNA-mediated therapies shows enormous promise in the treatment of these disorders. We discuss the overlapping post-transcriptional regulatory pathways that mediate muscle development, how these pathways are disrupted in neuromuscular disorders, and advances in RNA-mediated therapies that present a novel approach to the treatment of these diseases.

Intronic Determinants Coordinate Charme lncRNA Nuclear Activity through the Interaction with MATR3 and PTBP1

Chromatin architect of muscle expression (Charme) is a muscle-restricted long noncoding RNA (lncRNA) that plays an important role in myogenesis. Earlier evidence indicates that the nuclear Charme isoform, named pCharme, acts on the chromatin by assisting the formation of chromatin domains where myogenic transcription occurs. By combining RNA antisense purification (RAP) with mass spectrometry and loss-of-function analyses, we have now identified the proteins that assist these chromatin activities. These proteins—which include a sub-set of splicing regulators, principally PTBP1 and the multifunctional RNA/DNA binding protein MATR3—bind to sequences located within the alternatively spliced intron-1 to form nuclear aggregates. Consistent with the functional importance of pCharme interactome in vivo, a targeted deletion of the intron-1 by a CRISPR-Cas9 approach in mouse causes the release of pCharme from the chromatin and results in cardiac defects similar to what was observed upon knockout of the full-length transcript.

•

pCharme is the chromatin-retained isoform of the muscle-specific Charme lncRNA

•

Intronic signals coordinate the association of pCharme with MATR3 and PTBP1

•

The particle assembly prompts pCharme intron-1 chromatin retention

•

Deletion of the intron-1 by CRISPR-Cas9 leads to heart defects in mouse

RNA-binding motif protein 39 (RBM39): An emerging cancer target

RNA-binding motif protein 39 (RBM39) is an RNA-binding protein involved in transcriptional co-regulation and alternative RNA splicing. Recent studies have revealed that RBM39 is the unexpected target of aryl sulphonamides, which act as molecular glues between RBM39 and the DCAF15-associated E3 ubiquitin ligase complex leading to selective degradation of the target. Loss of RBM39 leads to aberrant splicing events and differential gene expression, thereby inhibiting cell cycle progression and causing tumour regression in a number of preclinical models. Many clinical studies have shown that aryl sulphonamides were well tolerated, but their clinical performance was limited due to an insufficient understanding of the target, RBM39 biology and a lack of predictive biomarkers. This review summarises the current knowledge of RBM39 function and discusses the therapeutic potential of this spliceosome target in cancer therapy.

A model of k-mer surprisal to quantify local sequence information content surrounding splice regions

Molecular sequences carry information. Analysis of sequence conservation between homologous loci is a proven approach with which to explore the information content of molecular sequences. This is often done using multiple sequence alignments to support comparisons between homologous loci. These methods therefore rely on sufficient underlying sequence similarity with which to construct a representative alignment. Here we describe a method using a formal metric of information, surprisal, to analyse biological sub-sequences without alignment constraints. We applied our model to the genomes of five different species to reveal similar patterns across a panel of eukaryotes. As the surprisal of a sub-sequence is inversely proportional to its occurrence within the genome, the optimal size of the sub-sequences was selected for each species under consideration. With the model optimized, we found a strong correlation between surprisal and CG dinucleotide usage. The utility of our model was tested by examining the sequences of genes known to undergo splicing. We demonstrate that our model can identify biological features of interest such as known donor and acceptor sites. Analysis across all annotated coding exon junctions in Homo sapiens reveals the information content of coding exons to be greater than the surrounding intron regions, a consequence of increased suppression of the CG dinucleotide in intronic space. Sequences within coding regions proximal to exon junctions exhibited novel patterns within DNA and coding mRNA that are not a function of the encoded amino acid sequence. Our findings are consistent with the presence of secondary information encoding features such as DNA and RNA binding sites, multiplexed through the coding sequence and independent of the information required to define the corresponding amino-acid sequence. We conclude that surprisal provides a complementary methodology with which to locate regions of interest in the genome, particularly in situations that lack an appropriate multiple sequence alignment.

First Identification of RNA-Binding Proteins That Regulate Alternative Exons in the Dystrophin Gene

The Duchenne muscular dystrophy (DMD) gene has a complex expression pattern regulated by multiple tissue-specific promoters and by alternative splicing (AS) of the resulting transcripts. Here, we used an RNAi-based approach coupled with DMD-targeted RNA-seq to identify RNA-binding proteins (RBPs) that regulate splicing of its skeletal muscle isoform (Dp427m) in a human muscular cell line. A total of 16 RBPs comprising the major regulators of muscle-specific splicing events were tested. We show that distinct combinations of RBPs maintain the correct inclusion in the Dp427m of exons that undergo spatio-temporal AS in other dystrophin isoforms. In particular, our findings revealed the complex networks of RBPs contributing to the splicing of the two short DMD exons 71 and 78, the inclusion of exon 78 in the adult Dp427m isoform being crucial for muscle function. Among the RBPs tested, QKI and DDX5/DDX17 proteins are important determinants of DMD exon inclusion. This is the first large-scale study to determine which RBP proteins act on the physiological splicing of the DMD gene. Our data shed light on molecular mechanisms contributing to the expression of the different dystrophin isoforms, which could be influenced by a change in the function or expression level of the identified RBPs.

Systematics for types and effects of RNA variations

Systematics is described for annotation of variations in RNA molecules. The conceptual framework is part of Variation Ontology (VariO) and facilitates depiction of types of variations, their functional and structural effects and other consequences in any RNA molecule in any organism. There are more than 150 RNA related VariO terms in seven levels, which can be further combined to generate even more complicated and detailed annotations. The terms are described together with examples, usually for variations and effects in human and in diseases. RNA variation type has two subcategories: variation classification and origin with subterms. Altogether six terms are available for function description. Several terms are available for affected RNA properties. The ontology contains also terms for structural description for affected RNA type, post-transcriptional RNA modifications, secondary and tertiary structure effects and RNA sugar variations. Together with the DNA and protein concepts and annotations, RNA terms allow comprehensive description of variations of genetic and non-genetic origin at all possible levels. The VariO annotations are readable both for humans and computer programs for advanced data integration and mining.

A long-read RNA-seq approach to identify novel transcripts of very large genes

RNA-seq is widely used for studying gene expression, but commonly used sequencing platforms produce short reads that only span up to two exon junctions per read. This makes it difficult to accurately determine the composition and phasing of exons within transcripts. Although long-read sequencing improves this issue, it is not amenable to precise quantitation, which limits its utility for differential expression studies. We used long-read isoform sequencing combined with a novel analysis approach to compare alternative splicing of large, repetitive structural genes in muscles. Analysis of muscle structural genes that produce medium (Nrap: 5 kb), large (Neb: 22 kb), and very large (Ttn: 106 kb) transcripts in cardiac muscle, and fast and slow skeletal muscles identified unannotated exons for each of these ubiquitous muscle genes. This also identified differential exon usage and phasing for these genes between the different muscle types. By mapping the in-phase transcript structures to known annotations, we also identified and quantified previously unannotated transcripts. Results were confirmed by endpoint PCR and Sanger sequencing, which revealed muscle-type-specific differential expression of these novel transcripts. The improved transcript identification and quantification shown by our approach removes previous impediments to studies aimed at quantitative differential expression of ultralong transcripts.

Identification of Prognostic Signatures of Alternative Splicing in Glioma

Alternative splicing (AS) is a ubiquitous mechanism in which pre-mRNA can be spliced into divergent variants and involved in carcinogenesis and progression in several cancers. In the present study, we systematically profiled prognostic AS signatures involving both low grade glioma (LGG) and glioblastoma (GBM) and investigated the association of AS signatures with tumor grade and IDH1 status in glioma. Percent spliced in (PSI) values and corresponding clinical data were obtained from TCGA SpliceSeq and TCGA data portal, respectively. Prognostic AS signatures were identified using univariate and stepwise multivariate Cox regression. Heatmap analysis was performed based on prognostic AS signatures. A prognostic signature was established with 69 and 88 AS events, including specific splicing events of MUTYH, STEAP3, and CTNNB1, in LGG and GBM cohorts, respectively. The area under the curve (AUC) of the prediction model was 0.968 at 2000 days of overall survival (OS) in the LGG cohort and 0.966 at 450 days of OS in the GBM cohort. In addition, these prognostic AS signatures could complement current molecular classification, such as IDH1 mutation, 1p/19q codeletion, and ATRX loss, of glioma and further identify potential subgroups of glioma with the same molecular features. In conclusion, our study systematically profiled prognostic AS events involving both low grade glioma and glioblastoma for the first time, which also shed light on the crosstalk between AS signatures and molecular features of glioma.

Alternative Splicing of the SLCO1B1 Gene: An Exploratory Analysis of Isoform Diversity in Pediatric Liver

The hepatic influx transporter OATP1B1 (SLCO1B1) plays an important role in the disposition of endogenous substrates and drugs prescribed to children. Alternative splicing increases the diversity of protein products from > 90% of human genes and may be triggered by developmental signals. As concentrations of several endogenous OATP1B1 substrates change during growth and development, with this exploratory study we investigated age-dependent alternative splicing of SLCO1B1 mRNA in 97 postmortem livers (fetus-adolescents). Twenty-seven splice variants were detected; 10 were confirmed by additional bioinformatic analyses and verified by quantitative polymerase chain reaction, and selected for detailed analysis based on relative abundance, association with age, and overlap with an adjacent gene. Two splice variants code for reference OATP1B1 protein, and eight code for truncated proteins. The expression of eight isoforms was associated with age. We conclude that alternative splicing of SLCO1B1 occurs frequently in children; although the functional consequences remain unknown, the data raise the possibility of a regulatory role for alternative splicing in mediating developmental changes in drug disposition.

Contributions of alternative splicing to muscle type development and function

Animals possess a wide variety of muscle types that support different kinds of movements. Different muscles have distinct locations, morphologies and contractile properties, raising the question of how muscle diversity is generated during development. Normal aging processes and muscle disorders differentially affect particular muscle types, thus understanding how muscles normally develop and are maintained provides insight into alterations in disease and senescence. As muscle structure and basic developmental mechanisms are highly conserved, many important insights into disease mechanisms in humans as well as into basic principles of muscle development have come from model organisms such as Drosophila, zebrafish and mouse. While transcriptional regulation has been characterized to play an important role in myogenesis, there is a growing recognition of the contributions of alternative splicing to myogenesis and the refinement of muscle function. Here we review our current understanding of muscle type specific alternative splicing, using examples of isoforms with distinct functions from both vertebrates and Drosophila. Future exploration of the vast potential of alternative splicing to fine-tune muscle development and function will likely uncover novel mechanisms of isoform-specific regulation and a more holistic understanding of muscle development, disease and aging.

Genome-wide analysis reveals the association between alternative splicing and DNA methylation across human solid tumors

Background

Dysregulation of alternative splicing (AS) is a critical signature of cancer. However, the regulatory mechanisms of cancer-specific AS events, especially the impact of DNA methylation, are poorly understood.

Methods

By using The Cancer Genome Atlas (TCGA) SpliceSeq and TCGA data for ten solid tumor types, association analysis was performed to characterize the potential link between cancer-specific AS and DNA methylation. Functional and pathway enrichment analyses were performed, and the protein-protein interaction (PPI) network was constructed with the String website. The prognostic analysis was carried out with multivariate Cox regressions models.

Results

15,818 AS events in 3955 annotated genes were identified across ten solid tumor types. The different DNA methylation patterns between tumor and normal tissues at the corresponding alternative spliced exon boundaries were shown, and 51.3% of CpG sites (CpGs) revealed hypomethylated in tumors. Notably, 607 CpGs were found to be highly correlated with 369 cancer-specific AS events after permutation tests. Among them, the hypomethylated CpGs account for 52.7%, and the number of down-regulated exons was 173. Furthermore, we found 38 AS events in 35 genes could serve as new molecular biomarkers to predict patient survival.

Conclusions

Our study described the relationship between DNA methylation and AS events across ten human solid tumor types and provided new insights into intragenic DNA methylation and exon usage during the AS process.

miR-1/206 downregulates splicing factor Srsf9 to promote C2C12 differentiation

Background

Myogenesis is driven by specific changes in the transcriptome that occur during the different stages of muscle differentiation. In addition to controlled transcriptional transitions, several other post-transcriptional mechanisms direct muscle differentiation. Both alternative splicing and miRNA activity regulate gene expression and production of specialized protein isoforms. Importantly, disruption of either process often results in severe phenotypes as reported for several muscle diseases. Thus, broadening our understanding of the post-transcriptional pathways that operate in muscles will lay the foundation for future therapeutic interventions.

Methods

We employed bioinformatics analysis in concert with the well-established C2C12 cell system for predicting and validating novel miR-1 and miR-206 targets engaged in muscle differentiation. We used reporter gene assays to test direct miRNA targeting and studied C2C12 cells stably expressing one of the cDNA candidates fused to a heterologous, miRNA-resistant 3′ UTR. We monitored effects on differentiation by measuring fusion index, myotube area, and myogenic gene expression during time course differentiation experiments.

Results

Gene ontology analysis revealed a strongly enriched set of putative miR-1 and miR-206 targets associated with RNA metabolism. Notably, the expression levels of several candidates decreased during C2C12 differentiation. We discovered that the splicing factor Srsf9 is a direct target of both miRNAs during myogenesis. Persistent Srsf9 expression during differentiation impaired myotube formation and blunted induction of the early pro-differentiation factor myogenin as well as the late differentiation marker sarcomeric myosin, Myh8.

Conclusions

Our data uncover novel miR-1 and miR-206 cellular targets and establish a functional link between the splicing factor Srsf9 and myoblast differentiation. The finding that miRNA-mediated clearance of Srsf9 is a key myogenic event illustrates the coordinated and sophisticated interplay between the diverse components of the gene regulatory network.

Rbfox3/NeuN Regulates Alternative Splicing of Tau Exon 10

Alternative splicing of tau exon 10 generates tau isoforms with three or four microtubule-binding repeats, 3R-tau or 4R-tau, which are under developmental regulation. Dysregulation of tau exon 10 splicing is sufficient to cause neurodegenerative disorders. The RNA-binding Fox3 (Rbfox3), identified as NeuN, regulates RNA processing. However, whether Rbfox3/NeuN regulates tau exon 10 splicing is unknown. In the present study, we found that the developmental expression of 4R-tau coincided with the expression of Rbfox3 in rat brains. Rbfox3 enhanced tau exon 10 inclusion. Tau intron 10 contains UGCAUG, the conservative binding sequence of Rbfox3. Intron 10 of tau pre-mRNA was co-immunoprecipitated by Rbfox3/NeuN. Deletion mutants of the RNA recognition motif (RRM) or three RNA-binding sites of the RRM in Rbfox3/NeuN failed to enhance tau exon 10 inclusion. Rbfox3, specifically expressed in the fetal brain, did not affect tau exon 10 splicing. The level of Rbfox3/NeuN was reduced and was associated with the ratio of 4R-tau/3R-tau in the excitotoxic mouse brains induced by kainic acid. These findings suggest that Rbfox3/NeuN regulates the alternative splicing of tau exon 10 and that decreased Rbfox3/NeuN may lower the ratio of 4R-tau/3R-tau.

Sequence and Evolutionary Features for the Alternatively Spliced Exons of Eukaryotic Genes

Alternative splicing of pre-mRNAs is a crucial mechanism for maintaining protein diversity in eukaryotes without requiring a considerable increase of genes in the number. Due to rapid advances in high-throughput sequencing technologies and computational algorithms, it is anticipated that alternative splicing events will be more intensively studied to address different kinds of biological questions. The occurrences of alternative splicing mean that all exons could be classified to be either constitutively or alternatively spliced depending on whether they are virtually included into all mature mRNAs. From an evolutionary point of view, therefore, the alternatively spliced exons would have been associated with distinctive biological characteristics in comparison with constitutively spliced exons. In this paper, we first outline the representative types of alternative splicing events and exon classification, and then review sequence and evolutionary features for the alternatively spliced exons. The main purpose is to facilitate understanding of the biological implications of alternative splicing in eukaryotes. This knowledge is also helpful to establish computational approaches for predicting the splicing pattern of exons.

Deep Splicing Code: Classifying Alternative Splicing Events Using Deep Learning

Alternative splicing (AS) is the process of combining different parts of the pre-mRNA to produce diverse transcripts and eventually different protein products from a single gene. In computational biology field, researchers try to understand AS behavior and regulation using computational models known as “Splicing Codes”. The final goal of these algorithms is to make an in-silico prediction of AS outcome from genomic sequence. Here, we develop a deep learning approach, called Deep Splicing Code (DSC), for categorizing the well-studied classes of AS namely alternatively skipped exons, alternative 5’ss, alternative 3’ss, and constitutively spliced exons based only on the sequence of the exon junctions. The proposed approach significantly improves the prediction and the obtained results reveal that constitutive exons have distinguishable local characteristics from alternatively spliced exons. Using the motif visualization technique, we show that the trained models learned to search for competitive alternative splice sites as well as motifs of important splicing factors with high precision. Thus, the proposed approach greatly expands the opportunities to improve alternative splicing modeling. In addition, a web-server for AS events prediction has been developed based on the proposed method.