An analysis of tissue-specific alternative splicing at the protein level

Systematic analysis of the effects of splicing on the diversity of post-translational modifications in protein isoforms using PTM-POSE

Post-translational modifications (PTMs) and splicing are important regulatory processes for controlling protein function and activity. Despite examples of interplay between alternative splicing and cell signaling in literature, there have been few detailed analyses of the impacts of alternative splicing on PTMs, partly due to difficulties in extracting PTM information from splicing measurements. We developed a computational pipeline, PTM Projection Onto Splice Events (PTM-POSE), to identify "prospective" PTM sites in alternative isoforms and splice events recorded in databases using only the genomic coordinates of a splice event or isoform of interest. Importantly, PTM-POSE integrates various PTM-specific databases and tools to allow for deeper analysis of the individual and global impact of spliced PTMs on isoform function, protein interactions, and regulation by enzymes like kinases. Using PTM-POSE, we performed a systematic analysis of PTM diversification across isoforms annotated in the Ensembl database. We found that 32% of PTMs are excluded from at least one Ensembl isoform, with palmitoylation being most likely to be excluded (49%) and glycosylation and crotonylation exhibiting the highest constitutive rates (75% and 94%, respectively). Further, approximately 2% of prospective PTM sites exhibited altered regulatory sequences surrounding the modification site, suggesting that regulatory or binding interactions might be different in these proteoforms. When comparing splicing of phosphorylation sites to measured phosphorylation abundance in KRAS-expressing lung cells, differential inclusion of phosphorylation sites correlated with phosphorylation levels, particularly for larger changes in inclusion (> 20%). To better understand how splicing diversification of PTMs may alter protein function and regulatory networks in specific biological contexts, we applied PTM-POSE to exon utilization measurements from TCGASpliceSeq of prostate tumor samples from The Cancer Genome Atlas (TCGA) and identified 1,489 PTMs impacted by ESRP1-correlated splicing, a splicing factor associated with worsened prognosis. We identified protein interaction and regulatory networks that may be rewired as a result of differential inclusion of PTM sites in ribosomal and cytoskeletal proteins. We also found instances in which ESRP1-mediated splicing impacted PTMs by altering flanking residues surrounding specific phosphorylation sites that may be targets of 14-3-3 proteins and SH2 domains. In addition, SGK1 signaling was found to be influenced by ESRP1 expression through increased inclusion of SGK1 substrates in ESRP1-expressing patients. Based on validation in a separate prostate cancer cohort from the Chinese Prostate Cancer Genome and EpiGenome Atlas (CPGEA), this correlated with increased phosphorylation of SGK1 substrates, particularly when SGK1 was predicted to be active. From this work, we highlighted the extensive splicing-control of PTM sites across the transcriptome and the novel information that can be gained through inclusion of PTMs in the analysis of alternative splicing. Importantly, we have provided a publicly available python package (PTM-POSE: https://github.com/NaegleLab/PTM-POSE) and all associated data for use by the broader scientific community to allow for continued exploration of the relationship between splicing and PTMs.

LincRNA-ASAO promotes dental pulp repair through interacting with PTBP1 to increase ALPL alternative splicing

BACKGROUND

Alternative splicing not only expands the genetic encoding of genes but also determines cellular activities. This study aimed to elucidate the regulation mechanism and biological functions of lincRNA-ASAO in the process of odontogenesis-related genes alternative splicing mediated odontogenic differentiation of hDPSCs.

METHODS

RACE, RNA-seq, FISH and bioinformatics techniques were used to identify novel lincRNA-ASAO. ALP staining, alizarin red staining, qRT-PCR and western blot were used to identify the role of lincRNA-ASAO in regulating the odontoblast differentiation of hDPSCs. The binding protein PTBP1 of lincRNA-ASAO was screened by RNA-Pulldown, protein profiling and bioinformatics. The target gene ALPL of lincRNA-ASAO/PTBP1 was identified by RNA-seq, bioinformatics technology and DNA agarose gel electrophoresis. FISH, IF, PAR-CLIP and bioinformatics techniques were used to determine the roles of lincRNA-ASAO, PTBP1 and ALPL pre-mRNA in the odontoblast differentiation of hDPSCs.

RESULTS

We identified a novel lincRNA-ASAO that could promote the odontogenic differentiation of human Dental Pulp Stem Cells (hDPSCs). And, the interaction between lincRNA-ASAO and alternative splicing factor PTBP1 promoted the odontoblast differentiation of hDPSCs. In addition, lincRNA-ASAO forms duplexes with ALPL pre-mRNA, targeting PTBP1 to exonic splicing silencer (ESS) of ALPL and regulating exon 2 skipping. Notably, lincRNA-ASAO/PTBP1 regulated ALPL production to increase the type 2 splice variant, which promoted the odontoblast differentiation of hDPSCs.

CONCLUSIONS

We have identified the novel lincRNA-ASAO, which can promote the odontoblast differentiation of hDPSCs. The mechanism study found that lincRNA-ASAO/PTBP1 mediated the exon 2 skipping of ALPL pre-mRNA, resulting in the type 2 splice variant of ALPL. Our results enrich the understanding of lncRNAs and alternative splicing in regulating the odontoblast differentiation of hDPSCs, and provide clues to improve the clinical therapeutic potential of hDPSCs for dental pulp restoration.

Alternative splicing of Clock transcript mediates the response of circadian clocks to temperature changes

Circadian clocks respond to temperature changes over the calendar year, allowing organisms to adjust their daily biological rhythms to optimize health and fitness. In Drosophila, seasonal adaptations are regulated by temperature-sensitive alternative splicing (AS) of period (per) and timeless (tim) genes that encode key transcriptional repressors of clock gene expression. Although Clock (Clk) gene encodes the critical activator of circadian gene expression, AS of its transcripts and its potential role in temperature regulation of clock function have not been explored. Here, we observed that Clk transcripts undergo temperature-sensitive AS. Specifically, cold temperature leads to the production of an alternative Clk transcript, hereinafter termed Clk-cold, which encodes a CLK isoform with an in-frame deletion of four amino acids proximal to the DNA binding domain. Notably, serine 13 (S13), which we found to be a CK1α-dependent phosphorylation site, is deleted in CLK-cold protein. We demonstrated that upon phosphorylation at CLK(S13), CLK-DNA interaction is reduced, thus decreasing transcriptional activity of CLK. This is in agreement with our findings that CLK occupancy at clock genes and transcriptional output are elevated at cold temperature likely due to higher amounts of CLK-cold isoforms that lack S13 residue. Finally, we showed that PER promotes CK1α-dependent phosphorylation of CLK(S13), supporting kinase-scaffolding role of repressor proteins as a conserved feature in the regulation of eukaryotic circadian clocks. This study provides insights into the complex collaboration between AS and phospho-regulation in shaping temperature responses of the circadian clock.

In silico and in cellulo approaches for functional annotation of human protein splice variants

The elegance of pre-mRNA splicing mechanisms continues to interest scientists even after over a half century, since the discovery of the fact that coding regions in genes are interrupted by non-coding sequences. The vast majority of human genes have several mRNA variants, coding structurally and functionally different protein isoforms in a tissue-specific manner and with a linkage to specific developmental stages of the organism. Alteration of splicing patterns shifts the balance of functionally distinct proteins in living systems, distorts normal molecular pathways, and may trigger the onset and progression of various pathologies. Over the past two decades, numerous studies have been conducted in various life sciences disciplines to deepen our understanding of splicing mechanisms and the extent of their impact on the functioning of living systems. This review aims to summarize experimental and computational approaches used to elucidate the functions of splice variants of a single gene based on our experience accumulated in the laboratory of interactomics of proteoforms at the Institute of Biomedical Chemistry (IBMC) and best global practices.

The Role of Sirtuin-1 Isoforms in Regulating Mitochondrial Function

The sirtuin-1 (SIRT1) gene contains multiple exons that usually undergo alternative splicing. The exclusion of one or more exons causes domain loss in the alternatively spliced isoforms and may change their functions. However, it is not completely established to what extent the loss of a non-catalytic domain could affect its regulatory function. Using muscle cells and SIRT1-knockout cells, we examined the function of the constitutively spliced isoform (SIRT1-v1) versus the alternatively spliced isoforms SIRT1-v2 and SIRT1-v3 that had lost part of the N-terminal region. Our data indicate that partial loss of the N-terminal domains in SIRT1-v2 and SIRT1-v3 attenuated their function. The full-length SIRT1-v1 significantly increased the oxidative phosphorylation and ATP production rate. Furthermore, SIRT1-v1 specifically upregulated the mitochondrial respiratory complex I without affecting the activity of complexes II, III, and IV. Additionally, domain loss affected the regulation of site-specific lysine acetylation in the histone H4 protein, the gene expression of respiratory complex I subunits, and the metabolic balance of oxidative phosphorylation versus glycolysis. Since alternatively spliced isoforms tend to increase with advancing age, the impact of SIRT1 isoforms on mitochondrial respiratory complexes warrants further investigation.

Tissue Usage Preference and Intrinsically Disordered Region Remodeling of Alternative Splicing Derived Proteoforms in the Heart

A computational analysis of mass spectrometry data was performed to uncover alternative splicing derived protein variants across chambers of the human heart. Evidence for 216 non-canonical isoforms was apparent in the atrium and the ventricle, including 52 isoforms not documented on SwissProt and recovered using an RNA sequencing derived database. Among non-canonical isoforms, 29 show signs of regulation based on statistically significant preferences in tissue usage, including a ventricular enriched protein isoform of tensin-1 (TNS1) and an atrium-enriched PDZ and LIM Domain 3 (PDLIM3) isoform 2 (PDLIM3-2/ALP-H). Examined variant regions that differ between alternative and canonical isoforms are highly enriched with intrinsically disordered regions. Moreover, over two-thirds of such regions are predicted to function in protein binding and RNA binding. The analysis here lends further credence to the notion that alternative splicing diversifies the proteome by rewiring intrinsically disordered regions, which are increasingly recognized to play important roles in the generation of biological function from protein sequences.

Mining alternative splicing patterns in scRNA-seq data using scASfind

Single-cell RNA-seq (scRNA-seq) is widely used for transcriptome profiling, but most analyses focus on gene-level events, with less attention devoted to alternative splicing. Here, we present scASfind, a novel computational method to allow for quantitative analysis of cell type-specific splicing events using full-length scRNA-seq data. ScASfind utilizes an efficient data structure to store the percent spliced-in value for each splicing event. This makes it possible to exhaustively search for patterns among all differential splicing events, allowing us to identify marker events, mutually exclusive events, and events involving large blocks of exons that are specific to one or more cell types.

Characterization of alternative splicing during mammalian brain development reveals the extent of isoform diversity and potential effects on protein structural changes

Regulation of gene expression is critical for fate commitment of stem and progenitor cells during tissue formation. In the context of mammalian brain development, a plethora of studies have described how changes in the expression of individual genes characterize cell types across ontogeny and phylogeny. However, little attention has been paid to the fact that different transcripts can arise from any given gene through alternative splicing (AS). Considered a key mechanism expanding transcriptome diversity during evolution, assessing the full potential of AS on isoform diversity and protein function has been notoriously difficult. Here, we capitalize on the use of a validated reporter mouse line to isolate neural stem cells, neurogenic progenitors and neurons during corticogenesis and combine the use of short- and long-read sequencing to reconstruct the full transcriptome diversity characterizing neurogenic commitment. Extending available transcriptional profiles of the mammalian brain by nearly 50,000 new isoforms, we found that neurogenic commitment is characterized by a progressive increase in exon inclusion resulting in the profound remodeling of the transcriptional profile of specific cortical cell types. Most importantly, we computationally infer the biological significance of AS on protein structure by using AlphaFold2, revealing how radical protein conformational changes can arise from subtle changes in isoforms sequence. Together, our study reveals that AS has a greater potential to impact protein diversity and function than previously thought, independently from changes in gene expression.

Alternative splicing of clock transcript mediates the response of circadian clocks to temperature changes

Circadian clocks respond to temperature changes over the calendar year, allowing organisms to adjust their daily biological rhythms to optimize health and fitness. In Drosophila, seasonal adaptations and temperature compensation are regulated by temperature-sensitive alternative splicing (AS) of period (per) and timeless (tim) genes that encode key transcriptional repressors of clock gene expression. Although clock (clk) gene encodes the critical activator of clock gene expression, AS of its transcripts and its potential role in temperature regulation of clock function have not been explored. We therefore sought to investigate whether clk exhibits AS in response to temperature and the functional changes of the differentially spliced transcripts. We observed that clk transcripts indeed undergo temperature-sensitive AS. Specifically, cold temperature leads to the production of an alternative clk transcript, hereinafter termed clk-cold, which encodes a CLK isoform with an in-frame deletion of four amino acids proximal to the DNA binding domain. Notably, serine 13 (S13), which we found to be a CK1α-dependent phosphorylation site, is among the four amino acids deleted in CLK-cold protein. Using a combination of transgenic fly, tissue culture, and in vitro experiments, we demonstrated that upon phosphorylation at CLK(S13), CLK-DNA interaction is reduced, thus decreasing CLK occupancy at clock gene promoters. This is in agreement with our findings that CLK occupancy at clock genes and transcriptional output are elevated at cold temperature, which can be explained by the higher amounts of CLK-cold isoforms that lack S13 residue. This study provides new insights into the complex collaboration between AS and phospho-regulation in shaping temperature responses of the circadian clock.

Small molecules modulating RNA splicing: a review of targets and future perspectives

In eukaryotic cells, RNA splicing is crucial for gene expression. Dysregulation of this process can result in incorrect mRNA processing, leading to aberrant gene expression patterns. Such abnormalities are implicated in many inherited diseases and cancers. Historically, antisense oligonucleotides, which bind to specific RNA targets, have been used to correct these splicing abnormalities. Despite their high specificity of action, these oligonucleotides have drawbacks, such as lack of oral bioavailability and the need for chemical modifications to enhance cellular uptake and stability. As a result, recent efforts focused on the development of small organic molecules that can correct abnormal RNA splicing event under disease conditions. This review discusses known and potential targets of these molecules, including RNA structures, trans-acting splicing factors, and the spliceosome - the macromolecular complex responsible for RNA splicing. We also rely on recent advances to discuss therapeutic applications of RNA-targeting small molecules in splicing correction. Overall, this review presents an update on strategies for RNA splicing modulation, emphasizing the therapeutic promise of small molecules.

GFAP-isoforms in the nervous system: Understanding the need for diversity

Glial fibrillary acidic protein (GFAP) is an intermediate filament (IF) protein expressed in specific types of glial cells in the nervous system. The expression of GFAP is highly regulated during brain development and in neurological diseases. The presence of distinct GFAP-isoforms in various cell types, developmental stages, and diseases indicates that GFAP (post-)transcriptional regulation has a role in glial cell physiology and pathology. GFAP-isoforms differ in sub-cellular localisation, IF-network assembly properties, and IF-dynamics which results in distinct molecular interactions and mechanical properties of the IF-network. Therefore, GFAP (post-)transcriptional regulation is likely a mechanism by which radial glia, astrocytes, and glioma cells can modulate cellular function.

Environmentally relevant concentrations of 2,3,7,8-TCDD induced inhibition of multicellular alternative splicing and transcriptional dysregulation

Alternative splicing (AS), found in approximately 95 % of human genes, significantly amplifies protein diversity and is implicated in disease pathogenesis when dysregulated. However, the precise involvement of AS in the toxic mechanisms induced by TCDD (2,3,7,8-tetrachlorodibenzo-p-dioxin) remains incompletely elucidated. This study conducted a thorough global AS analysis in six human cell lines following TCDD exposure. Our findings revealed that environmentally relevant concentration (0.1 nM) of TCDD significantly suppressed AS events in all cell types, notably inhibiting diverse splicing events and reducing transcript diversity, potentially attributed to modifications in the splicing patterns of the inhibitory factor family, particularly hnRNP. And we identified 151 genes with substantial AS alterations shared among these cell types, particularly enriched in immune and metabolic pathways. Moreover, TCDD induced cell-specific changes in splicing patterns and transcript levels, with increased sensitivity notably in THP-1 monocyte, potentially linked to aberrant expression of pivotal genes within the spliceosome pathway (DDX5, EFTUD2, PUF60, RBM25, SRSF1, and CRNKL1). This study extends our understanding of disrupted alternative splicing and its relation to the multisystem toxicity of TCDD. It sheds light on how environmental toxins affect post-transcriptional regulatory processes, offering a fresh perspective for toxicology and disease etiology investigations.

Alternative Transcripts Diversify Genome Function for Phenome Relevance to Health and Diseases

Manipulation using alternative exon splicing (AES), alternative transcription start (ATS), and alternative polyadenylation (APA) sites are key to transcript diversity underlying health and disease. All three are pervasive in organisms, present in at least 50% of human protein-coding genes. In fact, ATS and APA site use has the highest impact on protein identity, with their ability to alter which first and last exons are utilized as well as impacting stability and translation efficiency. These RNA variants have been shown to be highly specific, both in tissue type and stage, with demonstrated importance to cell proliferation, differentiation and the transition from fetal to adult cells. While alternative exon splicing has a limited effect on protein identity, its ubiquity highlights the importance of these minor alterations, which can alter other features such as localization. The three processes are also highly interwoven, with overlapping, complementary, and competing factors, RNA polymerase II and its CTD (C-terminal domain) chief among them. Their role in development means dysregulation leads to a wide variety of disorders and cancers, with some forms of disease disproportionately affected by specific mechanisms (AES, ATS, or APA). Challenges associated with the genome-wide profiling of RNA variants and their potential solutions are also discussed in this review.

Tissue Usage Preference and Intrinsically Disordered Region Remodeling of Alternative Splicing Derived Proteoforms in the Heart

A computational analysis of mass spectrometry data was performed to uncover alternative splicing derived protein variants across chambers of the human heart. Evidence for 216 non-canonical isoforms was apparent in the atrium and the ventricle, including 52 isoforms not documented on SwissProt and recovered using an RNA sequencing derived database. Among non-canonical isoforms, 29 show signs of regulation based on statistically significant preferences in tissue usage, including a ventricular enriched protein isoform of tensin-1 (TNS1) and an atrium-enriched PDZ and LIM Domain 3 (PDLIM3) isoform 2 (PDLIM3-2/ALP-H). Examined variant regions that differ between alternative and canonical isoforms are highly enriched in intrinsically disordered regions, and over two-thirds of such regions are predicted to function in protein binding and/or RNA binding. The analysis here lends further credence to the notion that alternative splicing diversifies the proteome by rewiring intrinsically disordered regions, which are increasingly recognized to play important roles in the generation of biological function from protein sequences.

Nanopore sequencing unveils the complexity of the cold-activated murine brown adipose tissue transcriptome

Alternative transcription increases transcriptome complexity by expression of multiple transcripts per gene. Annotation and quantification of transcripts using short-read sequencing is non-trivial. Long-read sequencing aims at overcoming these problems by sequencing full-length transcripts. Activation of brown adipose tissue (BAT) thermogenesis involves major transcriptomic remodeling and positively affects metabolism via increased energy expenditure. We benchmark Oxford Nanopore Technology (ONT) long-read sequencing protocols to Illumina short-read sequencing assessing alignment characteristics, gene and transcript detection and quantification, differential gene and transcript expression, transcriptome reannotation, and differential transcript usage (DTU). We find ONT sequencing is superior to Illumina for transcriptome reassembly, reducing the risk of false-positive events by unambiguously mapping reads to transcripts. We identified novel isoforms of genes undergoing DTU in cold-activated BAT including Cars2, Adtrp, Acsl5, Scp2, Aldoa, and Pde4d, validated by real-time PCR. The reannotated murine BAT transcriptome established here provides a framework for future investigations into the regulation of BAT.

Protein diversification through post-translational modifications, alternative splicing, and gene duplication

Proteins provide the basis for cellular function. Having multiple versions of the same protein within a single organism provides a way of regulating its activity or developing novel functions. Post-translational modifications of proteins, by means of adding/removing chemical groups to amino acids, allow for a well-regulated and controlled way of generating functionally distinct protein species. Alternative splicing is another method with which organisms possibly generate new isoforms. Additionally, gene duplication events throughout evolution generate multiple paralogs of the same genes, resulting in multiple versions of the same protein within an organism. In this review, we discuss recent advancements in the study of these three methods of protein diversification and provide illustrative examples of how they affect protein structure and function.

New insights into the centrosome-associated spliceosome components as regulators of ciliogenesis and tissue identity

Biomolecular condensates are membrane-less assemblies of proteins and nucleic acids. Centrosomes are biomolecular condensates that play a crucial role in nuclear division, cytoskeletal remodeling, and cilia formation in animal cells. Spatial omics technology is providing new insights into the dynamic exchange of spliceosome components between the nucleus and the centrosome/cilium. Intriguingly, centrosomes are emerging as cytoplasmic sites for information storage, enriched with RNA molecules and RNA-processing proteins. Furthermore, growing evidence supports the view that nuclear spliceosome components assembled at the centrosome function as potential coordinators of splicing subprograms, pluripotency, and cell differentiation. In this article, we first discuss the current understanding of the centrosome/cilium complex, which controls both stem cell differentiation and pluripotency. We next explore the molecular mechanisms that govern splicing factor assembly and disassembly around the centrosome and examine how RNA processing pathways contribute to ciliogenesis. Finally, we discuss numerous unresolved compelling questions regarding the centrosome-associated spliceosome components and transcript variants within the cytoplasm as sources of RNA-based secondary messages in the regulation of cell identity and cell fate determination. This article is categorized under: RNA-Based Catalysis > RNA Catalysis in Splicing and Translation RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes RNA Processing > Splicing Regulation/Alternative Splicing RNA Processing > RNA Processing.

The identification of protein and RNA interactors of the splicing factor Caper in the adult Drosophila nervous system

Post-transcriptional gene regulation is a fundamental mechanism that helps regulate the development and healthy aging of the nervous system. Mutations that disrupt the function of RNA-binding proteins (RBPs), which regulate post-transcriptional gene regulation, have increasingly been implicated in neurological disorders including amyotrophic lateral sclerosis, Fragile X Syndrome, and spinal muscular atrophy. Interestingly, although the majority of RBPs are expressed widely within diverse tissue types, the nervous system is often particularly sensitive to their dysfunction. It is therefore critical to elucidate how aberrant RNA regulation that results from the dysfunction of ubiquitously expressed RBPs leads to tissue specific pathologies that underlie neurological diseases. The highly conserved RBP and alternative splicing factor Caper is widely expressed throughout development and is required for the development of Drosophila sensory and motor neurons. Furthermore, caper dysfunction results in larval and adult locomotor deficits. Nonetheless, little is known about which proteins interact with Caper, and which RNAs are regulated by Caper. Here we identify proteins that interact with Caper in both neural and muscle tissue, along with neural specific Caper target RNAs. Furthermore, we show that a subset of these Caper-interacting proteins and RNAs genetically interact with caper to regulate Drosophila gravitaxis behavior.

SUsPECT: a pipeline for variant effect prediction based on custom long-read transcriptomes for improved clinical variant annotation

Our incomplete knowledge of the human transcriptome impairs the detection of disease-causing variants, in particular if they affect transcripts only expressed under certain conditions. These transcripts are often lacking from reference transcript sets, such as Ensembl/GENCODE and RefSeq, and could be relevant for establishing genetic diagnoses. We present SUsPECT (Solving Unsolved Patient Exomes/gEnomes using Custom Transcriptomes), a pipeline based on the Ensembl Variant Effect Predictor (VEP) to predict variant impact on custom transcript sets, such as those generated by long-read RNA-sequencing, for downstream prioritization. Our pipeline predicts the functional consequence and likely deleteriousness scores for missense variants in the context of novel open reading frames predicted from any transcriptome. We demonstrate the utility of SUsPECT by uncovering potential mutational mechanisms of pathogenic variants in ClinVar that are not predicted to be pathogenic using the reference transcript annotation. In further support of SUsPECT's utility, we identified an enrichment of immune-related variants predicted to have a more severe molecular consequence when annotating with a newly generated transcriptome from stimulated immune cells instead of the reference transcriptome. Our pipeline outputs crucial information for further prioritization of potentially disease-causing variants for any disease and will become increasingly useful as more long-read RNA sequencing datasets become available.

Alternative splicing shapes the transcriptome complexity in blackgram [Vigna mungo (L.) Hepper]

Vigna mungo, a highly consumed crop in the pan-Asian countries, is vulnerable to several biotic and abiotic stresses. Understanding the post-transcriptional gene regulatory cascades, especially alternative splicing (AS), may underpin large-scale genetic improvements to develop stress-resilient varieties. Herein, a transcriptome based approach was undertaken to decipher the genome-wide AS landscape and splicing dynamics in order to establish the intricacies of their functional interactions in various tissues and stresses. RNA sequencing followed by high-throughput computational analyses identified 54,526 AS events involving 15,506 AS genes that generated 57,405 transcripts isoforms. Enrichment analysis revealed their involvement in diverse regulatory functions and demonstrated that transcription factors are splicing-intensive, splice variants of which are expressed differentially across tissues and environmental cues. Increased expression of a splicing regulator NHP2L1/SNU13 was found to co-occur with lower intron retention events. The host transcriptome is significantly impacted by differential isoform expression of 1172 and 765 AS genes that resulted in 1227 (46.8% up and 53.2% downregulated) and 831 (47.5% up and 52.5% downregulated) transcript isoforms under viral pathogenesis and Fe²⁺ stressed condition, respectively. However, genes experiencing AS operate differently from the differentially expressed genes, suggesting AS is a unique and independent mode of regulatory mechanism. Therefore, it can be inferred that AS mediates a crucial regulatory role across tissues and stressful situations and the results would provide an invaluable resource for future endeavours in V. mungo genomics.

Definition of the transcriptional units of inherited retinal disease genes by meta-analysis of human retinal transcriptome data

BACKGROUND

Inherited retinal diseases (IRD) are genetically heterogeneous disorders that cause the dysfunction or loss of photoreceptor cells and ultimately lead to blindness. To date, next-generation sequencing procedures fail to detect pathogenic sequence variants in coding regions of known IRD disease genes in about 30-40% of patients. One of the possible explanations for this missing heritability is the presence of yet unidentified transcripts of known IRD genes. Here, we aimed to define the transcript composition of IRD genes in the human retina by a meta-analysis of publicly available RNA-seq datasets using an ad-hoc designed pipeline.

RESULTS

We analysed 218 IRD genes and identified 5,054 transcripts, 3,367 of which were not previously reported. We assessed their putative expression levels and focused our attention on 435 transcripts predicted to account for at least 5% of the expression of the corresponding gene. We looked at the possible impact of the newly identified transcripts at the protein level and experimentally validated a subset of them.

CONCLUSIONS

This study provides an unprecedented, detailed overview of the complexity of the human retinal transcriptome that can be instrumental in contributing to the resolution of some cases of missing heritability in IRD patients.

Regulation of pre-mRNA splicing: roles in physiology and disease, and therapeutic prospects

The removal of introns from mRNA precursors and its regulation by alternative splicing are key for eukaryotic gene expression and cellular function, as evidenced by the numerous pathologies induced or modified by splicing alterations. Major recent advances have been made in understanding the structures and functions of the splicing machinery, in the description and classification of physiological and pathological isoforms and in the development of the first therapies for genetic diseases based on modulation of splicing. Here, we review this progress and discuss important remaining challenges, including predicting splice sites from genomic sequences, understanding the variety of molecular mechanisms and logic of splicing regulation, and harnessing this knowledge for probing gene function and disease aetiology and for the design of novel therapeutic approaches.

The physiology of alternative splicing

Alternative splicing is a substantial contributor to the high complexity of transcriptomes of multicellular eukaryotes. In this Review, we discuss the accumulated evidence that most of this complexity is reflected at the protein level and fundamentally shapes the physiology and pathology of organisms. This notion is supported not only by genome-wide analyses but, mainly, by detailed studies showing that global and gene-specific modulations of alternative splicing regulate highly diverse processes such as tissue-specific and species-specific cell differentiation, thermal regulation, neuron self-avoidance, infrared sensing, the Warburg effect, maintenance of telomere length, cancer and autism spectrum disorders (ASD). We also discuss how mastering the control of alternative splicing paved the way to clinically approved therapies for hereditary diseases.

Identification of Alternative Splicing in Proteomes of Human Melanoma Cell Lines without RNA Sequencing Data

Alternative splicing is one of the main regulation pathways in living cells beyond simple changes in the level of protein expression. Most of the approaches proposed in proteomics for the identification of specific splicing isoforms require a preliminary deep transcriptomic analysis of the sample under study, which is not always available, especially in the case of the re-analysis of previously acquired data. Herein, we developed new algorithms for the identification and validation of protein splice isoforms in proteomic data in the absence of RNA sequencing of the samples under study. The bioinformatic approaches were tested on the results of proteome analysis of human melanoma cell lines, obtained earlier by high-resolution liquid chromatography and mass spectrometry (LC-MS). A search for alternative splicing events for each of the cell lines studied was performed against the database generated from all known transcripts (RefSeq) and the one composed of peptide sequences, which included all biologically possible combinations of exons. The identifications were filtered using the prediction of both retention times and relative intensities of fragment ions in the corresponding mass spectra. The fragmentation mass spectra corresponding to the discovered alternative splicing events were additionally examined for artifacts. Selected splicing events were further validated at the mRNA level by quantitative PCR.

GENCODE: reference annotation for the human and mouse genomes in 2023

GENCODE produces high quality gene and transcript annotation for the human and mouse genomes. All GENCODE annotation is supported by experimental data and serves as a reference for genome biology and clinical genomics. The GENCODE consortium generates targeted experimental data, develops bioinformatic tools and carries out analyses that, along with externally produced data and methods, support the identification and annotation of transcript structures and the determination of their function. Here, we present an update on the annotation of human and mouse genes, including developments in the tools, data, analyses and major collaborations which underpin this progress. For example, we report the creation of a set of non-canonical ORFs identified in GENCODE transcripts, the LRGASP collaboration to assess the use of long transcriptomic data to build transcript models, the progress in collaborations with RefSeq and UniProt to increase convergence in the annotation of human and mouse protein-coding genes, the propagation of GENCODE across the human pan-genome and the development of new tools to support annotation of regulatory features by GENCODE. Our annotation is accessible via Ensembl, the UCSC Genome Browser and https://www.gencodegenes.org.

Reversal of splicing infidelity is a pre-activation step in B cell differentiation

Introduction

B cell activation and differentiation is central to the adaptive immune response. Changes in exon usage can have major impacts on cellular signaling and differentiation but have not been systematically explored in differentiating B cells.

Methods

We analyzed exon usage and intron retention in RNA-Seq data from subsets of human B cells at various stages of differentiation, and in an in vitro laboratory model of B cell activation and differentiation (Epstein Barr virus infection).

Results

Blood naïve B cells were found to have an unusual splicing profile, with unannotated splicing events in over 30% of expressed genes. Splicing changed substantially upon naïve B cell entry into secondary lymphoid tissue and before activation, involving significant increases in exon commitment and reductions in intron retention. These changes preferentially involved short introns with weak splice sites and were likely mediated by an overall increase in splicing efficiency induced by the lymphoid environment. The majority of transcripts affected by splicing changes showed restoration of encoded conserved protein domains and/or reduced targeting to the nonsense-mediated decay pathway. Affected genes were enriched in functionally important immune cell activation pathways such as antigen-mediated signaling, cell cycle control and mRNA processing and splicing.

Discussion

Functional observations from donor B cell subsets in progressive states of differentiation and from timecourse experiments using the in vitro model suggest that these widespread changes in mRNA splicing play a role in preparing naïve B cells for the decisive step of antigen-mediated activation and differentiation.

Where protein structure and cell diversity meet

Protein-protein interaction networks - interactomes - are charted with the hope to understand how phenotypes emerge and how they are altered in disease states. Early efforts to map interactomes have focused on the assembly of context agnostic, reference networks. However, recent studies have mapped interactomes across different cell lines and tissues, finding highly variable interactomes due to the rewiring of protein-protein interactions in different contexts. Increasing evidence points to significant links between protein structure and interactome diversity seen across cell types and tissues. We discuss how recent findings support the key role of alternative splicing and phosphorylation, two well-established regulators of protein structural and functional diversity, in defining cell type- and tissue-specific interactomes. Moreover, we show that intrinsically disordered protein regions are most favorably equipped to support interactome rewiring by acting as hubs of protein structure and function regulation.

Origins and Evolution of Human Tandem Duplicated Exon Substitution Events

The mutually exclusive splicing of tandem duplicated exons produces protein isoforms that are identical save for a homologous region that allows for the fine tuning of protein function. Tandem duplicated exon substitution events are rare, yet highly important alternative splicing events. Most events are ancient, their isoforms are highly expressed, and they have significantly more pathogenic mutations than other splice events. Here, we analyzed the physicochemical properties and functional roles of the homologous polypeptide regions produced by the 236 tandem duplicated exon substitutions annotated in the human gene set. We find that the most important structural and functional residues in these homologous regions are maintained, and that most changes are conservative rather than drastic. Three quarters of the isoforms produced from tandem duplicated exon substitution events are tissue-specific, particularly in nervous and cardiac tissues, and tandem duplicated exon substitution events are enriched in functional terms related to structures in the brain and skeletal muscle. We find considerable evidence for the convergent evolution of tandem duplicated exon substitution events in vertebrates, arthropods, and nematodes. Twelve human gene families have orthologues with tandem duplicated exon substitution events in both Drosophila melanogaster and Caenorhabditis elegans. Six of these gene families are ion transporters, suggesting that tandem exon duplication in genes that control the flow of ions into the cell has an adaptive benefit. The ancient origins, the strong indications of tissue-specific functions, and the evidence of convergent evolution suggest that these events may have played important roles in the evolution of animal tissues and organs.

Clinical variant interpretation and biologically relevant reference transcripts

Clinical variant interpretation is highly dependent on the choice of reference transcript. Although the longest transcript has traditionally been chosen as the reference, APPRIS principal and MANE Select transcripts, biologically supported reference sequences, are now available. In this study, we show that MANE Select and APPRIS principal transcripts are the best reference transcripts for clinical variation. APPRIS principal and MANE Select transcripts capture almost all ClinVar pathogenic variants, and they are particularly powerful over the 94% of coding genes in which they agree. We find that a vanishingly small number of ClinVar pathogenic variants affect alternative protein products. Alternative isoforms that are likely to be clinically relevant can be predicted using TRIFID scores, the highest scoring alternative transcripts are almost 700 times more likely to house pathogenic variants. We believe that APPRIS, MANE and TRIFID are essential tools for clinical variant interpretation.

Differential Expression Pattern of linc-ROR Spliced Variants in Pluripotent and Non-Pluripotent Cell Lines

OBJECTIVE

The human large intergenic non-coding RNA-regulator of reprogramming program (linc-ROR) is known as a stem cell specific linc-RNA. linc-ROR counteracts differentiation via sequestering microRNA-145 (miR-145) that targets OCT4 transcript. Despite the research on the expression and function, the exact structure of linc-ROR transcripts is not clear. Considering the contribution of alternative splicing in transcripts structures and function, identifying different spliced variants of linc-ROR is necessary for further functional analyses. We aimed to find the alternatively spliced transcripts of linc-ROR and investigate their expression pattern in stem and cancer cell lines and during neural differentiation of NT2 cells as a model for understanding linc-ROR role in stem cell and differentiation.

MATERIALS AND METHODS

In this experimental study, linc-ROR locus was scanned for identifying novel exons. Different primer sets were used to detect new spliced variants by reverse transcription polymerase chain reaction (RT-PCR) and direct sequencing. Quantitative PCR (qPCR) and RT-PCR were employed to profile expression of linc-ROR transcripts in different cell lines and during neural differentiation of stem cells.

RESULTS

We could discover 13 novel spliced variants of linc-ROR harboring unique array of exons. Our work uncovered six novel exons, some of which were the product of exonized transposable elements. Monitoring expression profile of the linc-ROR spliced variants in a panel of pluripotent and non-pluripotent cells exhibited that all transcripts were primarily expressed in pluripotent cells. Moreover, the examined linc-ROR spliced variants showed a similar downregulation during neural differentiation of NT2 cells.

CONCLUSION

Altogether, our data showed despite the difference in the structure and composition of exons, various spliced variants of linc-ROR showed similar expression pattern in stem cells and through differentiation.

TEx-MST: tissue expression profiles of MANE select transcripts

Recently, a new reference transcript dataset [Matched Annotation from the NCBI and EMBL-EBI (MANE) select] was released by NCBI and EMBL-EBI to make available a new unified representative transcript for human protein-coding genes. While the main purpose of MANE project is to provide a harmonized gene and transcript information standard, there is no explicit tissue expression information about these MANE select transcripts. In this report, we tried to provide useful expression profiles of MANE select transcripts in various normal human tissues to allow further interrogation of their molecular modulations and functional significance. We obtained the new V9 transcript expression dataset from the Genotype-Tissue Expression (GTEx) web portal. This new GTEx dataset, based on a long-read sequencing platform, affords better assessment of the expression of alternative spliced transcripts. This tissue expression profiles of MANE select transcripts (TEx-MST) database not only provides the basic information of MANE select transcripts but also tissue expression profiles on alternative transcripts in protein-coding genes. Users can initiate the interrogation by gene symbol searches or by browsing the MANE genes with various criteria (such as genome locations or expression rankings). We further utilized the GENCODE biotype feature to identify the top-ranked protein-coding transcripts by choosing the most expressed protein-coding transcripts from GTEx datasets (both V8 and V9 datasets). In summary, there are 18 083 genes matched between MANE and GTEx. Among them, 13 245 MANE select transcripts matched with the top-ranked protein-coding transcripts in GTEx V9 dataset, which underlined the dominate expression of MANE select transcripts. This TEx-MST web bioinformatic database provides a visualized user interface for the normal tissue expression patterns of MANE select transcripts using the newly released GTEx dataset.

Database URL: TEx-MST is available at https://texmst.ibms.sinica.edu.tw/

APPRIS principal isoforms and MANE Select transcripts define reference splice variants

MOTIVATION

Selecting the splice variant that best represents a coding gene is a crucial first step in many experimental analyses, and vital for mapping clinically relevant variants. This study compares the longest isoforms, MANE Select transcripts, APPRIS principal isoforms, and expression data, and aims to determine which method is best for selecting biological important reference splice variants for large-scale analyses.

RESULTS

Proteomics analyses and human genetic variation data suggest that most coding genes have a single main protein isoform. We show that APPRIS principal isoforms and MANE Select transcripts best describe these main cellular isoforms, and find that using the longest splice variant as the representative is a poor strategy. Exons unique to the longest splice isoforms are not under selective pressure, and so are unlikely to be functionally relevant. Expression data are also a poor means of selecting the main splice variant. APPRIS principal and MANE Select exons are under purifying selection, while exons specific to alternative transcripts are not. There are MANE and APPRIS representatives for almost 95% of genes, and where they agree they are particularly effective, coinciding with the main proteomics isoform for over 98.2% of genes.

AVAILABILITY AND IMPLEMENTATION

APPRIS principal isoforms for human, mouse and other model species can be downloaded from the APPRIS database (https://appris.bioinfo.cnio.es), GENCODE genes (https://www.gencodegenes.org/) and the Ensembl website (https://www.ensembl.org). MANE Select transcripts for the human reference set are available from the Ensembl, GENCODE and RefSeq databases (https://www.ncbi.nlm.nih.gov/refseq/). Lists of splice variants where MANE and APPRIS coincide are available from the APPRIS database.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The whole blood transcriptional regulation landscape in 465 COVID-19 infected samples from Japan COVID-19 Task Force

Coronavirus disease 2019 (COVID-19) is a recently-emerged infectious disease that has caused millions of deaths, where comprehensive understanding of disease mechanisms is still unestablished. In particular, studies of gene expression dynamics and regulation landscape in COVID-19 infected individuals are limited. Here, we report on a thorough analysis of whole blood RNA-seq data from 465 genotyped samples from the Japan COVID-19 Task Force, including 359 severe and 106 non-severe COVID-19 cases. We discover 1169 putative causal expression quantitative trait loci (eQTLs) including 34 possible colocalizations with biobank fine-mapping results of hematopoietic traits in a Japanese population, 1549 putative causal splice QTLs (sQTLs; e.g. two independent sQTLs at TOR1AIP1), as well as biologically interpretable trans-eQTL examples (e.g., REST and STING1), all fine-mapped at single variant resolution. We perform differential gene expression analysis to elucidate 198 genes with increased expression in severe COVID-19 cases and enriched for innate immune-related functions. Finally, we evaluate the limited but non-zero effect of COVID-19 phenotype on eQTL discovery, and highlight the presence of COVID-19 severity-interaction eQTLs (ieQTLs; e.g., CLEC4C and MYBL2). Our study provides a comprehensive catalog of whole blood regulatory variants in Japanese, as well as a reference for transcriptional landscapes in response to COVID-19 infection.

Proteotranscriptomics - A facilitator in omics research

Applications in omics research, such as comparative transcriptomics and proteomics, require the knowledge of the species-specific gene sequence and benefit from a comprehensive high-quality annotation of the coding genes to achieve high coverage. While protein-coding genes can in simple cases be detected by scanning the genome for open reading frames, in more complex genomes exonic sequences are separated by introns. Despite advances in sequencing technologies that allow for ever-growing numbers of genomes, the quality of many of the provided genome assemblies do not reach reference quality. These non-contiguous assemblies with gaps and the necessity to predict splice sites limit accurate gene annotation from solely genomic data. In contrast, the transcriptome only contains transcribed gene regions, is devoid of introns and thus provides the optimal basis for the identification of open reading frames. The additional integration of proteomics data to validate predicted protein-coding genes further enriches for accurate gene models. This review outlines the principles of the proteotranscriptomics approach, discusses common challenges and suggests methods for improvement.

Alternative splicing as a source of phenotypic diversity

A major goal of evolutionary genetics is to understand the genetic processes that give rise to phenotypic diversity in multicellular organisms. Alternative splicing generates multiple transcripts from a single gene, enriching the diversity of proteins and phenotypic traits. It is well established that alternative splicing contributes to key innovations over long evolutionary timescales, such as brain development in bilaterians. However, recent developments in long-read sequencing and the generation of high-quality genome assemblies for diverse organisms has facilitated comparisons of splicing profiles between closely related species, providing insights into how alternative splicing evolves over shorter timescales. Although most splicing variants are probably non-functional, alternative splicing is nonetheless emerging as a dynamic, evolutionarily labile process that can facilitate adaptation and contribute to species divergence.

Dominant transcript expression profiles of human protein-coding genes interrogated with GTEx dataset

The discovery and quantification of mRNA transcripts using short-read next-generation sequencing (NGS) data is a complicated task. There are far more alternative mRNA transcripts expressed by human genes than can be identified from NGS transcriptome data and various bioinformatic pipelines, while the numbers of annotated human protein-coding genes has gradually declined in recent years. It is essential to learn more about the thorough tissue expression profiles of alternative transcripts in order to obtain their molecular modulations and actual functional significance. In this report, we present a bioinformatic database for interrogating the representative tissue of human protein-coding transcripts. The database allows researchers to visually explore the top-ranked transcript expression profiles in particular tissue types. Most transcripts of protein-coding genes were found to have certain tissue expression patterns. This observation demonstrated that many alternative transcripts were particularly modulated in different cell types. This user-friendly tool visually represents transcript expression profiles in a tissue-specific manner. Identification of tissue specific protein-coding genes and transcripts is a substantial advance towards interpreting their biological functions and further functional genomics studies.

NeuroSCORE is a genome-wide omics-based model that identifies candidate disease genes of the central nervous system

To identify candidate disease genes of central nervous system (CNS) phenotypes, we created the Neurogenetic Systematic Correlation of Omics-Related Evidence (NeuroSCORE). We identified five genome-wide metrics highly associated with CNS phenotypes to score 19,601 protein-coding genes. Genes scored one point per metric (range: 0-5), identifying 8298 scored genes (scores ≥ 1) and 1601 "high scoring" genes (scores ≥ 3). Using logistic regression, we determined the odds ratio that genes with a NeuroSCORE from 1 to 5 would be associated with known CNS-related phenotypes compared to genes that scored zero. We tested NeuroSCORE using microarray copy number variants (CNVs) in case-control cohorts and aggregate mouse model data. High scoring genes are associated with CNS phenotypes (OR = 5.5, p < 2E-16), enriched in case CNVs, and mouse ortholog genes that cause behavioral and nervous system abnormalities. We identified 1058 high scoring genes with no disease association in OMIM. Transforming the logistic regression results indicates high scoring genes have an 84-92% chance of being associated with a CNS phenotype. Top scoring genes include GRIA1, MAP4K4, SF1, TNPO2, and ZSWIM8. Finally, we interrogated CNVs in the Clinical Genome Resource, finding the majority of clinically significant CNVs contain high scoring genes. These findings can direct future research and improve molecular diagnostics.

ErbB4 alternative splicing mediates fetal mouse alveolar type II cell differentiation in vitro

BACKGROUND

Alternative splicing (AS) creates different protein isoforms, an important mechanism regulating cell-specific function. Little is known about AS in lung development, particularly in alveolar type II (ATII) cells. ErbB4 receptor isoforms Jma and Jmb have significant and opposing functions in the brain, heart, and lung development and/or disease. However, the regulators of ErbB4 AS are unknown. ErbB4 AS regulators in fetal mouse ATII cells control its function in ATII cell maturation.

METHODS

Candidate ErbB4 AS regulators were found using in silico analysis. Their developmental expression was studied in fetal mouse ATII cells. The effects of splice factor downregulation and upregulation on ATII cell maturation were analyzed.

RESULTS

ErbB4-Jma increased significantly in ATII cells after gestation E16.5. In silico analysis found four candidate splice factors: FOX2, CUG/CELF1, TIAR, and HUB. Fetal ATII cells expressed these factors in distinct developmental profiles. HUB downregulation in E17.5 ATII cells increased Jma isoform levels and Sftpb gene expression and decreased Jmb. HUB overexpression decreased Jma and Sftpb.

CONCLUSIONS

ErbB4 AS is developmentally controlled by HUB in fetal ATII cells, promoting ATII differentiation. Regulated AS expression during ATII cell differentiation suggests novel therapeutic strategies to approach human disease.

IMPACT

Alternative splicing (AS) of the ErbB4 receptor, involving mutually exclusive exon inclusion, creates Jma and Jmb isoforms with distinct differences in receptor processing and function. The Jma isoform of ErbB4 promotes differentiation of fetal lung alveolar type II cells. The AS is mediated in part by the RNA-binding protein HUB. The molecular mechanism of AS for ErbB4 has not been previously described. The regulation of ErbB4 AS has important implications in the development of organs, such as the lung, brain, and heart, and for disease, including cancer.

Characterization of Alternative Splicing Events in Porcine Skeletal Muscles with Different Intramuscular Fat Contents

Meat quality is one of the most important economic traits in pig breeding and production. Intramuscular fat (IMF) is a major factor that improves meat quality. To better understand the alternative splicing (AS) events underlying meat quality, long-read isoform sequencing (Iso-seq) was used to identify differential (D)AS events between the longissimus thoracis (LT) and semitendinosus (ST), which differ in IMF content, together with short-read RNA-seq. Through Iso-seq analysis, we identified a total of 56,789 novel transcripts covering protein-coding genes, lncRNA, and fusion transcripts that were not previously annotated in pigs. We also identified 456,965 AS events, among which 3930 were DAS events, corresponding to 2364 unique genes. Through integrative analysis of Iso-seq and RNA-seq, we identified 1174 differentially expressed genes (DEGs), among which 122 were DAS genes, i.e., DE-DAS genes. There are 12 overlapped pathways between the top 20 DEGs and DE-DAS genes, as revealed by KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis, indicating that DE-DAS genes play important roles in the differential phenotype of LT and ST. Further analysis showed that upregulated DE-DAS genes are more important than downregulated ones in IMF deposition. Fatty acid degradation and the PPAR (peroxisome proliferator-activated receptor) signaling pathway were found to be the most important pathways regulating the differential fat deposition of the two muscles. The results update the existing porcine genome annotations and provide data for the in-depth exploration of the mechanisms underlying meat quality and IMF deposition.

APPRIS: selecting functionally important isoforms

APPRIS (https://appris.bioinfo.cnio.es) is a well-established database housing annotations for protein isoforms for a range of species. APPRIS selects principal isoforms based on protein structure and function features and on cross-species conservation. Most coding genes produce a single main protein isoform and the principal isoforms chosen by the APPRIS database best represent this main cellular isoform. Human genetic data, experimental protein evidence and the distribution of clinical variants all support the relevance of APPRIS principal isoforms. APPRIS annotations and principal isoforms have now been expanded to 10 model organisms. In this paper we highlight the most recent updates to the database. APPRIS annotations have been generated for two new species, cow and chicken, the protein structural information has been augmented with reliable models from the EMBL-EBI AlphaFold database, and we have substantially expanded the confirmatory proteomics evidence available for the human genome. The most significant change in APPRIS has been the implementation of TRIFID functional isoform scores. TRIFID functional scores are assigned to all splice isoforms, and APPRIS uses the TRIFID functional scores and proteomics evidence to determine principal isoforms when core methods cannot.

BaseScope™ Approach to Visualize Alternative Splice Variants in Tissue

Defining the cell-specific alternative splicing landscape in complex tissues is an important goal to gain functional insights. Deep-sequencing techniques coupled to genetic strategies for cell identification has provided important cues on cell-specific exon usage in complex tissues like the nervous system. BaseScope™ has emerged as a powerful and highly sensitive alternative to in situ hybridization to determine exon composition in tissue with spatial and morphological context. In this protocol, we will review how BaseScope was utilized to detect the e37a-Cacna1b splice variant of the presynaptic calcium channel CaV2.2 or N-type. This splice variant arises from a pair of mutually exclusive exons (e37a and e37b). E37a-Cacna1b is heavily underrepresented relative to e37b-Cacna1b and both exons share 60% of their sequence. By using BaseScope™, we were able to discover that e37a-Cacna1b is expressed in excitatory pyramidal neurons of hippocampus and cortex, as well as motor neurons of the ventral horn of the spinal cord.

How alternative splicing changes the properties of plant proteins

Most plant primary transcripts undergo alternative splicing (AS), and its impact on protein diversity is a subject of intensive investigation. Several studies have uncovered various mechanisms of how particular protein splice isoforms operate. However, the common principles behind the AS effects on protein function in plants have rarely been surveyed. Here, on the selected examples, we highlight diverse tissue expression patterns, subcellular localization, enzymatic activities, abilities to bind other molecules and other relevant features. We describe how the protein isoforms mutually interact to underline their intriguing roles in altering the functionality of protein complexes. Moreover, we also discuss the known cases when these interactions have been placed inside the autoregulatory loops. This review is particularly intended for plant cell and developmental biologists who would like to gain inspiration on how the splice variants encoded by their genes of interest may coordinately work.

Transcriptome-Wide Analyses Identify Dominant as the Predominantly Non-Conservative Alternative Splicing Inheritance Patterns in F1 Chickens

Transcriptome analysis has been used to investigate many economically traits in chickens; however, alternative splicing still lacks a systematic method of study that is able to promote proteome diversity, and fine-tune expression dynamics. Hybridization has been widely utilized in chicken breeding due to the resulting heterosis, but the dynamic changes in alternative splicing during this process are significant yet unclear. In this study, we performed a reciprocal crossing experiment involving the White Leghorn and Cornish Game chicken breeds which exhibit major differences in body size and reproductive traits, and conducted RNA sequencing of the brain, muscle, and liver tissues to identify the inheritance patterns. A total of 40 515 and 42 612 events were respectively detected in the brain and muscle tissues, with 39 843 observed in the liver; 2807, 4242, and 4538 events significantly different between two breeds were identified in the brain, muscle, and liver tissues, respectively. The hierarchical cluster of tissues from different tissues from all crosses, based on the alternative splicing profiles, suggests high tissue and strain specificity. Furthermore, a comparison between parental strains and hybrid crosses indicated that over one third of alternative splicing genes showed conserved patterns in all three tissues, while the second prevalent pattern was non-additive, which included both dominant and transgressive patterns; this meant that the dominant pattern plays a more important role than suppression. Our study provides an overview of the inheritance patterns of alternative splicing in layer and broiler chickens, to better understand post-transcriptional regulation during hybridization.

Exploring the Diverse Functional and Regulatory Consequences of Alternative Splicing in Development and Disease

Alternative splicing is a fundamental mechanism of eukaryotic RNA regulation that increases the transcriptomic and proteomic complexity within an organism. Moreover, alternative splicing provides a framework for generating unique yet complex tissue- and cell type-specific gene expression profiles, despite using a limited number of genes. Recent efforts to understand the negative consequences of aberrant splicing have increased our understanding of developmental and neurodegenerative diseases such as spinal muscular atrophy, frontotemporal dementia and Parkinsonism linked to chromosome 17, myotonic dystrophy, and amyotrophic lateral sclerosis. Moreover, these studies have led to the development of innovative therapeutic treatments for diseases caused by aberrant splicing, also known as spliceopathies. Despite this, a paucity of information exists on the physiological roles and specific functions of distinct transcript spliceforms for a given gene. Here, we will highlight work that has specifically explored the distinct functions of protein-coding spliceforms during development. Moreover, we will discuss the use of alternative splicing of noncoding exons to regulate the stability and localization of RNA transcripts.

The Full-Length Transcriptome Provides New Insights Into the Transcript Complexity of Abdominal Adipose and Subcutaneous Adipose in Pekin Ducks

Adipose tissues have a central role in organisms, and adipose content is a crucial economic trait of poultry. Pekin duck is an ideal model to study the mechanism of abdominal and subcutaneous adipose deposition for its high ability of adipose synthesis and deposition. Alternative splicing contributes to functional diversity in abdominal and subcutaneous adipose. However, there has been no systematic analysis of the dynamics of differential alternative splicing of abdominal and subcutaneous adipose in Pekin duck. In our study, the Pacific Biosciences (PacBio) Iso-Seq technology was applied to explore the transcriptional complexity of abdominal and subcutaneous adipose in Pekin ducks. In total, 143,931 and 111,337 full-length non-chimeric transcriptome sequences of abdominal and subcutaneous adipocytes were obtained from 41.78 GB raw data, respectively. These data led us to identify 19,212 long non-coding RNAs (lncRNAs) and 74,571 alternative splicing events. In addition, combined with the next-generation sequencing technology, we correlated the structure and function annotation with the differential expression profiles of abdominal and subcutaneous adipose transcripts. This study identified lots of novel alternative splicing events and major transcripts of transcription factors related to adipose synthesis. STAT3 was reported as a vital gene for adipogenesis, and we found that its major transcript is STAT3-1, which may play a considerable role in the process of adipose synthesis in Pekin duck. This study greatly increases our understanding of the gene models, genome annotations, genome structures, and the complexity and diversity of abdominal and subcutaneous adipose in Pekin duck. These data provide insights into the regulation of alternative splicing events, which form an essential part of transcript diversity during adipogenesis in poultry. The results of this study provide an invaluable resource for studying alternative splicing and tissue-specific expression.

Functional enrichment of alternative splicing events with NEASE reveals insights into tissue identity and diseases

Alternative splicing (AS) is an important aspect of gene regulation. Nevertheless, its role in molecular processes and pathobiology is far from understood. A roadblock is that tools for the functional analysis of AS-set events are lacking. To mitigate this, we developed NEASE, a tool integrating pathways with structural annotations of protein-protein interactions to functionally characterize AS events. We show in four application cases how NEASE can identify pathways contributing to tissue identity and cell type development, and how it highlights splicing-related biomarkers. With a unique view on AS, NEASE generates unique and meaningful biological insights complementary to classical pathways analysis.

Molecular Profiling of DNA Methylation and Alternative Splicing of Genes in Skeletal Muscle of Obese Rabbits

DNA methylation and the alternative splicing of precursor messenger RNAs (pre-mRNAs) are two important genetic modification mechanisms. However, both are currently uncharacterized in the muscle metabolism of rabbits. Thus, we constructed the Tianfu black rabbit obesity model (obese rabbits fed with a 10% high-fat diet and control rabbits from 35 days to 70 days) and collected the skeletal muscle samples from the two groups for Genome methylation sequencing and RNA sequencing. DNA methylation data showed that the promoter regions of 599 genes and gene body region of 2522 genes had significantly differential methylation rates between the two groups, of which 288 genes had differential methylation rates in promoter and gene body regions. Analysis of alternative splicing showed 555 genes involved in exon skipping (ES) patterns, and 15 genes existed in differential methylation regions. Network analysis showed that 20 hub genes were associated with ubiquitinated protein degradation, muscle development pathways, and skeletal muscle energy metabolism. Our findings suggest that the two types of genetic modification have potential regulatory effects on skeletal muscle development and provide a basis for further mechanistic studies in the rabbit.

The clinical importance of tandem exon duplication-derived substitutions

Most coding genes in the human genome are annotated with multiple alternative transcripts. However, clear evidence for the functional relevance of the protein isoforms produced by these alternative transcripts is often hard to find. Alternative isoforms generated from tandem exon duplication-derived substitutions are an exception. These splice events are rare, but have important functional consequences. Here, we have catalogued the 236 tandem exon duplication-derived substitutions annotated in the GENCODE human reference set. We find that more than 90% of the events have a last common ancestor in teleost fish, so are at least 425 million years old, and twenty-one can be traced back to the Bilateria clade. Alternative isoforms generated from tandem exon duplication-derived substitutions also have significantly more clinical impact than other alternative isoforms. Tandem exon duplication-derived substitutions have >25 times as many pathogenic and likely pathogenic mutations as other alternative events. Tandem exon duplication-derived substitutions appear to have vital functional roles in the cell and may have played a prominent part in metazoan evolution.