Genome-wide detection of tissue-specific alternative splicing in the human transcriptome

A New Phase for WNK Kinase Signaling Complexes as Biomolecular Condensates

The purpose of this review is to highlight transformative advances that have been made in the field of biomolecular condensates, with special emphasis on condensate material properties, physiology, and kinases, using the With-No-Lysine (WNK) kinases as a prototypical example. To convey how WNK kinases illustrate important concepts for biomolecular condensates, we start with a brief history, focus on defining features of biomolecular condensates, and delve into some examples of how condensates are implicated in cellular physiology (and pathophysiology). We then highlight how WNK kinases, through the action of "WNK droplets" that ubiquitously regulate intracellular volume and kidney-specific "WNK bodies" that are implicated in distal tubule salt reabsorption and potassium homeostasis, exemplify many of the defining features of condensates. Finally, this review addresses the controversies within this emerging field and questions to address.

Regulation of DNA damage response by RNA/DNA-binding proteins: Implications for neurological disorders and aging

RNA-binding proteins (RBPs) are evolutionarily conserved across most forms of life, with an estimated 1500 RBPs in humans. Traditionally associated with post-transcriptional gene regulation, RBPs contribute to nearly every known aspect of RNA biology, including RNA splicing, transport, and decay. In recent years, an increasing subset of RBPs have been recognized for their DNA binding properties and involvement in DNA transactions. We refer to these RBPs with well-characterized DNA binding activity as RNA/DNA binding proteins (RDBPs), many of which are linked to neurological diseases. RDBPs are associated with both nuclear and mitochondrial DNA repair. Furthermore, the presence of intrinsically disordered domains in RDBPs appears to be critical for regulating their diverse interactions and plays a key role in controlling protein aggregation, which is implicated in neurodegeneration. In this review, we discuss the emerging roles of common RDBPs from the heterogeneous nuclear ribonucleoprotein (hnRNP) family, such as TAR DNA binding protein-43 (TDP43) and fused in sarcoma (FUS) in controlling DNA damage response (DDR). We also explore the implications of RDBP pathology in aging and neurodegenerative diseases and provide a prospective on the therapeutic potential of targeting RDBP pathology mediated DDR defects for motor neuron diseases and aging.

Single-cell long-read targeted sequencing reveals transcriptional variation in ovarian cancer

Single-cell RNA sequencing predominantly employs short-read sequencing to characterize cell types, states and dynamics; however, it is inadequate for comprehensive characterization of RNA isoforms. Long-read sequencing technologies enable single-cell RNA isoform detection but are hampered by lower throughput and unintended sequencing of artifacts. Here we develop Single-cell Targeted Isoform Long-Read Sequencing (scTaILoR-seq), a hybridization capture method which targets over a thousand genes of interest, improving the median number of on-target transcripts per cell by 29-fold. We use scTaILoR-seq to identify and quantify RNA isoforms from ovarian cancer cell lines and primary tumors, yielding 10,796 single-cell transcriptomes. Using long-read variant calling we reveal associations of expressed single nucleotide variants (SNVs) with alternative transcript structures. Phasing of SNVs across transcripts enables the measurement of allelic imbalance within distinct cell populations. Overall, scTaILoR-seq is a long-read targeted RNA sequencing method and analytical framework for exploring transcriptional variation at single-cell resolution.

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.

Long-read RNA sequencing identifies region- and sex-specific C57BL/6J mouse brain mRNA isoform expression and usage

Alternative splicing (AS) contributes to the biological heterogeneity between species, sexes, tissues, and cell types. Many diseases are either caused by alterations in AS or by alterations to AS. Therefore, measuring AS accurately and efficiently is critical for assessing molecular phenotypes, including those associated with disease. Long-read sequencing enables more accurate quantification of differentially spliced isoform expression than short-read sequencing approaches, and third-generation platforms facilitate high-throughput experiments. To assess differences in AS across the cerebellum, cortex, hippocampus, and striatum by sex, we generated and analyzed Oxford Nanopore Technologies (ONT) long-read RNA sequencing (lrRNA-Seq) C57BL/6J mouse brain cDNA libraries. From > 85 million reads that passed quality control metrics, we calculated differential gene expression (DGE), differential transcript expression (DTE), and differential transcript usage (DTU) across brain regions and by sex. We found significant DGE, DTE, and DTU across brain regions and that the cerebellum had the most differences compared to the other three regions. Additionally, we found region-specific differential splicing between sexes, with the most sex differences in DTU in the cortex and no DTU in the hippocampus. We also report on two distinct patterns of sex DTU we observed, sex-divergent and sex-specific, that could potentially help explain sex differences in the prevalence and prognosis of various neurological and psychiatric disorders in future studies. Finally, we built a Shiny web application for researchers to explore the data further. Our study provides a resource for the community; it underscores the importance of AS in biological heterogeneity and the utility of long-read sequencing to better understand AS in the brain.

The Roles of hnRNP Family in the Brain and Brain-Related Disorders

Heterogeneous nuclear ribonucleoproteins (hnRNPs) belong to a complex family of RNA-binding proteins that are essential to control alternative splicing, mRNA trafficking, synaptic plasticity, stress granule formation, cell cycle regulation, and axonal transport. Over the past decade, hnRNPs have been associated with different brain disorders such as Alzheimer's disease, multiple sclerosis, and schizophrenia. Given their essential role in maintaining cell function and integrity, it is not surprising that dysregulated hnRNP levels lead to neurological implications. This review aims to explore the primary functions of hnRNPs in neurons, oligodendrocytes, microglia, and astrocytes, and their roles in brain disorders. We also discuss proteomics and other technologies and their potential for studying and evaluating hnRNPs in brain disorders, including the discovery of new therapeutic targets and possible pharmacological interventions.

Quantifying transcriptome diversity: a review

Following the central dogma of molecular biology, gene expression heterogeneity can aid in predicting and explaining the wide variety of protein products, functions and, ultimately, heterogeneity in phenotypes. There is currently overlapping terminology used to describe the types of diversity in gene expression profiles, and overlooking these nuances can misrepresent important biological information. Here, we describe transcriptome diversity as a measure of the heterogeneity in (1) the expression of all genes within a sample or a single gene across samples in a population (gene-level diversity) or (2) the isoform-specific expression of a given gene (isoform-level diversity). We first overview modulators and quantification of transcriptome diversity at the gene level. Then, we discuss the role alternative splicing plays in driving transcript isoform-level diversity and how it can be quantified. Additionally, we overview computational resources for calculating gene-level and isoform-level diversity for high-throughput sequencing data. Finally, we discuss future applications of transcriptome diversity. This review provides a comprehensive overview of how gene expression diversity arises, and how measuring it determines a more complete picture of heterogeneity across proteins, cells, tissues, organisms and species.

Long-read RNA sequencing identifies region- and sex-specific C57BL/6J mouse brain mRNA isoform expression and usage

Alternative splicing (AS) contributes to the biological heterogeneity between species, sexes, tissues, and cell types. Many diseases are either caused by alterations in AS or by alterations to AS. Therefore, measuring AS accurately and efficiently is critical for assessing molecular phenotypes, including those associated with disease. Long-read sequencing enables more accurate quantification of differentially spliced isoform expression than short-read sequencing approaches, and third-generation platforms facilitate high-throughput experiments. To assess differences in AS across the cerebellum, cortex, hippocampus, and striatum by sex, we generated and analyzed Oxford Nanopore Technologies (ONT) long-read RNA sequencing (lrRNA-Seq) C57BL/6J mouse brain cDNA libraries. From >85 million reads that passed quality control metrics, we calculated differential gene expression (DGE), differential transcript expression (DTE), and differential transcript usage (DTU) across brain regions and by sex. We found significant DGE, DTE, and DTU across brain regions and that the cerebellum had the most differences compared to the other three regions. Additionally, we found region-specific differential splicing between sexes, with the most sex differences in DTU in the cortex and no DTU in the hippocampus. We also report on two distinct patterns of sex DTU we observed, sex-divergent and sex-specific, that could potentially help explain sex differences in the prevalence and prognosis of various neurological and psychiatric disorders in future studies. Finally, we built a Shiny web application for researchers to explore the data further. Our study provides a resource for the community; it underscores the importance of AS in biological heterogeneity and the utility of long-read sequencing to better understand AS in the brain.

Specific recognition and sensitive quantification of mRNA splice variants via one-pot ligation-dependent loop-mediated isothermal amplification

Specific recognition and sensitive quantification of mRNA alternative splice variants have been a necessity for exploring the regulatory mechanism of RNA splicing and revealing the association between pre-mRNA splicing and transcriptome function, as well as disease diagnosis. However, their wide abundance range and high sequence homology pose enormous challenges for high sensitivity and selectivity quantification of splice variants. Herein, taking advantage of the excellent specificity of ligation and the powerful nucleic acid replication feature of loop-mediated isothermal amplification (LAMP), we developed a one-pot method (termed one-pot ligation-LAMP) for specific recognition and sensitive quantification of mRNA splicing variants based on two splicing junction-specific stem-loop DNA probe ligation and the subsequently initiating LAMP. The one-pot ligation-LAMP can specifically detect as low as 100 aM mRNA splice variants without any nonspecific signals and quantify them with a wide dynamics range spanning at least six orders of magnitude. We have demonstrated that the one-pot ligation-LAMP is a versatile and practical strategy for accurately quantifying different splicing variants in complex biological samples with high sensitivity all in one tube within 90 min, thereby providing an attractive tool for mRNA splice variant-related studies.

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

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

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.

A splicing transcriptome-wide association study identifies novel altered splicing for Alzheimer's disease susceptibility

Alzheimer's disease (AD) is a common neurodegenerative disease in aging individuals. Alternative splicing is reported to be relevant to AD development while their roles in etiology of AD remain largely elusive. We performed a comprehensive splicing transcriptome-wide association study (spTWAS) using intronic excision expression genetic prediction models of 12 brain tissues developed through three modelling strategies, to identify candidate susceptibility splicing introns for AD risk. A total of 111,326 (46,828 proxy) cases and 677,663 controls of European ancestry were studied. We identified 343 associations of 233 splicing introns (143 genes) with AD risk after Bonferroni correction (0.05/136,884 = 3.65 × 10^-7). Fine-mapping analyses supported 155 likely causal associations corresponding to 83 splicing introns of 55 genes. Eighteen causal splicing introns of 15 novel genes (EIF2D, WDR33, SAP130, BYSL, EPHB6, MRPL43, VEGFB, PPP1R13B, TLN2, CLUHP3, LRRC37A4P, CRHR1, LINC02210, ZNF45-AS1, and XPNPEP3) were identified for the first time to be related to AD susceptibility. Our study identified novel genes and splicing introns associated with AD risk, which can improve our understanding of the etiology of AD.

Hidden Genetic Regulation of Human Complex Traits via Brain Isoforms

Alternative splicing exists in most multi-exonic genes, and exploring these complex alternative splicing events and their resultant isoform expressions is essential. However, it has become conventional that RNA sequencing results have often been summarized into gene-level expression counts mainly due to the multiple ambiguous mapping of reads at highly similar regions. Transcript-level quantification and interpretation are often overlooked, and biological interpretations are often deduced based on combined transcript information at the gene level. Here, for the most variable tissue of alternative splicing, the brain, we estimate isoform expressions in 1,191 samples collected by the Genotype-Tissue Expression (GTEx) Consortium using a powerful method that we previously developed. We perform genome-wide association scans on the isoform ratios per gene and identify isoform-ratio quantitative trait loci (irQTL), which could not be detected by studying gene-level expressions alone. By analyzing the genetic architecture of the irQTL, we show that isoform ratios regulate educational attainment via multiple tissues including the frontal cortex (BA9), cortex, cervical spinal cord, and hippocampus. These tissues are also associated with different neuro-related traits, including Alzheimer's or dementia, mood swings, sleep duration, alcohol intake, intelligence, anxiety or depression, etc. Mendelian randomization (MR) analysis revealed 1,139 pairs of isoforms and neuro-related traits with plausible causal relationships, showing much stronger causal effects than on general diseases measured in the UK Biobank (UKB). Our results highlight essential transcript-level biomarkers in the human brain for neuro-related complex traits and diseases, which could be missed by merely investigating overall gene expressions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s43657-023-00100-6.

Cell-type-specific alternative splicing in the Arabidopsis germline

During sexual reproduction in flowering plants, the two haploid sperm cells (SCs) embedded within the cytoplasm of a growing pollen tube are carried to the embryo sac for double fertilization. Pollen development in flowering plants is a dynamic process that encompasses changes at transcriptome and epigenome levels. While the transcriptome of pollen and SCs in Arabidopsis (Arabidopsis thaliana) is well documented, previous analyses have mostly been based on gene-level expression. In-depth transcriptome analysis, particularly the extent of alternative splicing (AS) at the resolution of SC and vegetative nucleus (VN), is still lacking. Therefore, we performed RNA-seq analysis to generate a spliceome map of Arabidopsis SCs and VN isolated from mature pollen grains. Based on our de novo transcriptome assembly, we identified 58,039 transcripts, including 9,681 novel transcripts, of which 2,091 were expressed in SCs and 3,600 in VN. Four hundred and sixty-eight genes were regulated both at gene and splicing levels, with many having functions in mRNA splicing, chromatin modification, and protein localization. Moreover, a comparison with egg cell RNA-seq data uncovered sex-specific regulation of transcription and splicing factors. Our study provides insights into a gamete-specific AS landscape at unprecedented resolution.

Characterization of an RNA binding protein interactome reveals a context-specific post-transcriptional landscape of MYC-amplified medulloblastoma

Pediatric medulloblastoma (MB) is the most common solid malignant brain neoplasm, with Group 3 (G3) MB representing the most aggressive subgroup. MYC amplification is an independent poor prognostic factor in G3 MB, however, therapeutic targeting of the MYC pathway remains limited and alternative therapies for G3 MB are urgently needed. Here we show that the RNA-binding protein, Musashi-1 (MSI1) is an essential mediator of G3 MB in both MYC-overexpressing mouse models and patient-derived xenografts. MSI1 inhibition abrogates tumor initiation and significantly prolongs survival in both models. We identify binding targets of MSI1 in normal neural and G3 MB stem cells and then cross referenced these data with unbiased large-scale screens at the transcriptomic, translatomic and proteomic levels to systematically dissect its functional role. Comparative integrative multi-omic analyses of these large datasets reveal cancer-selective MSI1-bound targets sharing multiple MYC associated pathways, providing a valuable resource for context-specific therapeutic targeting of G3 MB.

SRPK1 regulates RNA binding in a pre-spliceosomal complex using a catalytic bypass mechanism

Serine-arginine protein kinase 1 (SRPK1) phosphorylates serine-arginine (SR) proteins in the cytoplasm, directing them to the nucleus for splicing function. SRPK1 has also been detected in the nucleus but its function here is still not fully understood. We now demonstrate that nuclear SRPK1 can regulate U1-70K, a protein component of the uridine-rich 1 small nuclear ribonucleoprotein (U1 snRNP) that binds SR proteins and facilitates 5' splice-site selection in precursor mRNA. We found that SRPK1 uses a large, disordered domain to bind U1-70K, regulating the interaction of an exonic splicing enhancer (ESE) to the associated SR protein. Surprisingly, the catalytic activity of SRPK1 is not required for this phenomenon. Instead, SRPK1 associates directly with the N-terminus of U1-70K and alters the regulatory function of the distal C-terminus, modifying interactions between the U1-70K:SR protein complex and the ESE. Disruption of SRPK1 binding to this complex affects the alternative splicing of genes modulated by the C-terminus of U1-70K. Such findings suggest that, in addition to operating as a traditional serine-modifying catalyst, SRPK1 can also bypass this intrinsic activity to regulate RNA contacts in an early pre-spliceosomal complex.

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.

Comprehensive Transcriptome Analysis Reveals Sex-Specific Alternative Splicing Events in Zebrafish Gonads

Alternative splicing is an important way of regulating gene functions in eukaryotes. Several key genes involved in sex determination and gonadal differentiation, such as nr5a1 and ddx4, have sex-biased transcripts between males and females, suggesting a potential regulatory role of alternative splicing in gonads. Currently, the sex-specific alternative splicing events and genes have not been comprehensively studied at the genome-wide level in zebrafish. In this study, through global splicing analysis on three independent sets of RNA-seq data from matched zebrafish testes and ovaries, we identified 120 differentially spliced genes shared by the three datasets, most of which haven’t been reported before. Functional enrichment analysis showed that the GO terms of mRNA processing, mRNA metabolism and microtubule-based process were strongly enriched. The testis- and ovary-biased alternative splicing genes were identified, and part of them (tp53bp1, tpx2, mapre1a, kif2c, and ncoa5) were further validated by RT-PCR. Sequence characteristics analysis suggested that the lengths, GC contents, and splice site strengths of the alternative exons or introns may have different influences in different types of alternative splicing events. Interestingly, we identified an unexpected high proportion (over 70%) of non-frameshift exon-skipping events, suggesting that in these cases the two protein isoforms derived from alternative splicing may both have functions. Furthermore, as a representative example, we found that the alternative splicing of ncoa5 causes the loss of a conserved RRM domain in the short transcript predominantly produced in testes. Our study discovers novel sex-specific alternative splicing events and genes with high reliabilities in zebrafish testes and ovaries, which would provide attractive targets for follow-up studies to reveal the biological significances of alternative splicing events and genes in sex determination and gonadal differentiation.

The Post-Translational Modification Networking in WNK-Centric Hypertension Regulation and Electrolyte Homeostasis

The with-no-lysine (WNK) kinase family, comprising four serine-threonine protein kinases (WNK1-4), were first linked to hypertension due to their mutations in association with pseudohypoaldosteronism type II (PHAII). WNK kinases regulate crucial blood pressure regulators, SPAK/OSR1, to mediate the post-translational modifications (PTMs) of their downstream ion channel substrates, such as sodium chloride co-transporter (NCC), epithelial sodium chloride (ENaC), renal outer medullary potassium channel (ROMK), and Na/K/2Cl co-transporters (NKCCs). In this review, we summarize the molecular pathways dysregulating the WNKs and their downstream target renal ion transporters. We summarize each of the genetic variants of WNK kinases and the small molecule inhibitors that have been discovered to regulate blood pressure via WNK-triggered PTM cascades.

Screening and Identification of Survival-Associated Splicing Factors in Lung Squamous Cell Carcinoma

Lung squamous cell carcinoma (LUSC) is a disease with high morbidity and mortality. Many studies have shown that aberrant alternative splicing (AS) can lead to tumorigenesis, and splicing factors (SFs) serve as an important function during AS. In this research, we propose an analysis method based on synergy to screen key factors that regulate the initiation and progression of LUSC. We first screened alternative splicing events (ASEs) associated with survival in LUSC patients by bivariate Cox regression analysis. Then an association network consisting of OS-ASEs, SFs, and their targeting relationship was constructed to identify key SFs. Finally, 10 key SFs were selected in terms of degree centrality. The validation on TCGA and cross-platform GEO datasets showed that some SFs were significantly differentially expressed in cancer and paracancer tissues, and some of them were associated with prognosis, indicating that our method is valid and accurate. It is expected that our method would be applied to a wide range of research fields and provide new insights in the future.

OUP accepted manuscript

Alternative splicing is frequently involved in the diversification of protein function and can also be modulated for therapeutic purposes. Here we develop a predictive model, called Exon ByPASS (predicting Exon skipping Based on Protein amino acid SequenceS), to assess the criticality of exon inclusion based solely on information contained in the amino acid sequence upstream and downstream of the exon junctions. By focusing on protein sequence, Exon ByPASS predicts exon skipping independent of tissue and species in the absence of any intronic information. We validate model predictions using transcriptomic and proteomic data and show that the model can capture exon skipping in different tissues and species. Additionally, we reveal potential therapeutic opportunities by predicting synthetically skippable exons and neo-junctions arising in cancer cells.

Genome annotation with long RNA reads reveals new patterns of gene expression and improves single-cell analyses in an ant brain

BACKGROUND

Functional genomic analyses rely on high-quality genome assemblies and annotations. Highly contiguous genome assemblies have become available for a variety of species, but accurate and complete annotation of gene models, inclusive of alternative splice isoforms and transcription start and termination sites, remains difficult with traditional approaches.

RESULTS

Here, we utilized full-length isoform sequencing (Iso-Seq), a long-read RNA sequencing technology, to obtain a comprehensive annotation of the transcriptome of the ant Harpegnathos saltator. The improved genome annotations include additional splice isoforms and extended 3' untranslated regions for more than 4000 genes. Reanalysis of RNA-seq experiments using these annotations revealed several genes with caste-specific differential expression and tissue- or caste-specific splicing patterns that were missed in previous analyses. The extended 3' untranslated regions afforded great improvements in the analysis of existing single-cell RNA-seq data, resulting in the recovery of the transcriptomes of 18% more cells. The deeper single-cell transcriptomes obtained with these new annotations allowed us to identify additional markers for several cell types in the ant brain, as well as genes differentially expressed across castes in specific cell types.

CONCLUSIONS

Our results demonstrate that Iso-Seq is an efficient and effective approach to improve genome annotations and maximize the amount of information that can be obtained from existing and future genomic datasets in Harpegnathos and other organisms.

RNA Targeting in Inherited Neuromuscular Disorders: Novel Therapeutic Strategies to Counteract Mis-Splicing

Neuromuscular disorders represent multifaceted abnormal conditions, with little or no cure, leading to patient deaths from complete muscle wasting and atrophy. Despite strong efforts in the past decades, development of effective treatments is still urgently needed. Advent of next-generation sequencing technologies has allowed identification of novel genes and mutations associated with neuromuscular pathologies, highlighting splicing defects as essential players. Deciphering the significance and relative contributions of defective RNA metabolism will be instrumental to address and counteract these malignancies. We review here recent progress on the role played by alternative splicing in ensuring functional neuromuscular junctions (NMJs), and its involvement in the pathogenesis of NMJ-related neuromuscular disorders, with particular emphasis on congenital myasthenic syndromes and muscular dystrophies. We will also discuss novel strategies based on oligonucleotides designed to bind their cognate sequences in the RNA or targeting intermediary of mRNA metabolism. These efforts resulted in several chemical classes of RNA molecules that have recently proven to be clinically effective, more potent and better tolerated than previous strategies.

hnRNP A1B, a Splice Variant of HNRNPA1, Is Spatially and Temporally Regulated

RNA binding proteins (RBPs) play a key role in cellular growth, homoeostasis and survival and are tightly regulated. A deep understanding of their spatiotemporal regulation is needed to understand their contribution to physiology and pathology. Here, we have characterized the spatiotemporal expression pattern of hnRNP A1 and its splice variant hnRNP A1B in mice. We have found that hnRNP A1B expression is more restricted to the CNS compared to hnRNP A1, and that it can form an SDS-resistant dimer in the CNS. Also, hnRNP A1B expression becomes progressively restricted to motor neurons in the ventral horn of the spinal cord, compared to hnRNP A1 which is more broadly expressed. We also demonstrate that hnRNP A1B is present in neuronal processes, while hnRNP A1 is absent. This finding supports a hypothesis that hnRNP A1B may have a cytosolic function in neurons that is not shared with hnRNP A1. Our results demonstrate that both isoforms are differentially expressed across tissues and have distinct localization profiles, suggesting that the two isoforms may have specific subcellular functions that can uniquely contribute to disease progression.

Alternative Splicing in Cardiovascular Disease-A Survey of Recent Findings

Alternative splicing, a driver of posttranscriptional variance, differs from canonical splicing by arranging the introns and exons of an immature pre-mRNA transcript in a multitude of different ways. Although alternative splicing was discovered almost half a century ago, estimates of the proportion of genes that undergo alternative splicing have risen drastically over the last two decades. Deep sequencing methods and novel bioinformatic algorithms have led to new insights into the prevalence of spliced variants, tissue-specific splicing patterns and the significance of alternative splicing in development and disease. Thus far, the role of alternative splicing has been uncovered in areas ranging from heart development, the response to myocardial infarction to cardiac structural disease. Circular RNAs, a product of alternative back-splicing, were initially discovered in 1976, but landmark publications have only recently identified their regulatory role, tissue-specific expression, and transcriptomic abundance, spurring a renewed interest in the topic. The aim of this review is to provide a brief insight into some of the available findings on the role of alternative splicing in cardiovascular disease, with a focus on atherosclerosis, myocardial infarction, heart failure, dilated cardiomyopathy and circular RNAs in myocardial infarction.

Long-read RNA sequencing reveals widespread sex-specific alternative splicing in threespine stickleback fish

Alternate isoforms are important contributors to phenotypic diversity across eukaryotes. Although short-read RNA-sequencing has increased our understanding of isoform diversity, it is challenging to accurately detect full-length transcripts, preventing the identification of many alternate isoforms. Long-read sequencing technologies have made it possible to sequence full-length alternative transcripts, accurately characterizing alternative splicing events, alternate transcription start and end sites, and differences in UTR regions. Here, we use Pacific Biosciences (PacBio) long-read RNA-sequencing (Iso-Seq) to examine the transcriptomes of five organs in threespine stickleback fish (Gasterosteus aculeatus), a widely used genetic model species. The threespine stickleback fish has a refined genome assembly in which gene annotations are based on short-read RNA sequencing and predictions from coding sequence of other species. This suggests some of the existing annotations may be inaccurate or alternative transcripts may not be fully characterized. Using Iso-Seq we detected thousands of novel isoforms, indicating many isoforms are absent in the current Ensembl gene annotations. In addition, we refined many of the existing annotations within the genome. We noted many improperly positioned transcription start sites that were refined with long-read sequencing. The Iso-Seq-predicted transcription start sites were more accurate and verified through ATAC-seq. We also detected many alternative splicing events between sexes and across organs. We found a substantial number of genes in both somatic and gonadal samples that had sex-specific isoforms. Our study highlights the power of long-read sequencing to study the complexity of transcriptomes, greatly improving genomic resources for the threespine stickleback fish.

The RNA-Binding Protein HuD Regulates Alternative Splicing and Alternative Polyadenylation in the Mouse Neocortex

The neuronal Hu/ELAV-like proteins HuB, HuC and HuD are a class of RNA-binding proteins that are crucial for proper development and maintenance of the nervous system. These proteins bind to AU-rich elements (AREs) in the untranslated regions (3'-UTRs) of target mRNAs regulating mRNA stability, transport and translation. In addition to these cytoplasmic functions, Hu proteins have been implicated in alternative splicing and alternative polyadenylation in the nucleus. The purpose of this study was to identify transcriptome-wide effects of HuD deletion on both of these nuclear events using RNA sequencing data obtained from the neocortex of Elavl4-/- (HuD KO) mice. HuD KO affected alternative splicing of 310 genes, including 17 validated HuD targets such as Cbx3, Cspp1, Snap25 and Gria2. In addition, deletion of HuD affected polyadenylation of 53 genes, with the majority of significantly altered mRNAs shifting towards usage of proximal polyadenylation signals (PAS), resulting in shorter 3'-UTRs. None of these genes overlapped with those showing alternative splicing events. Overall, HuD KO had a greater effect on alternative splicing than polyadenylation, with many of the affected genes implicated in several neuronal functions and neuropsychiatric disorders.

Mutant PRPF8 Causes Widespread Splicing Changes in Spliceosome Components in Retinitis Pigmentosa Patient iPSC-Derived RPE Cells

Retinitis pigmentosa (RP) is a rare, progressive disease that affects photoreceptors and retinal pigment epithelial (RPE) cells with blindness as a final outcome. Despite high medical and social impact, there is currently no therapeutic options to slow down the progression of or cure the disease. The development of effective therapies was largely hindered by high genetic heterogeneity, inaccessible disease tissue, and unfaithful model organisms. The fact that components of ubiquitously expressed splicing factors lead to the retina-specific disease is an additional intriguing question. Herein, we sought to correlate the retinal cell-type-specific disease phenotype with the splicing profile shown by a patient with autosomal recessive RP, caused by a mutation in pre-mRNA splicing factor 8 (PRPF8). In order to get insight into the role of PRPF8 in homeostasis and disease, we capitalize on the ability to generate patient-specific RPE cells and reveal differentially expressed genes unique to RPE cells. We found that spliceosomal complex and ribosomal functions are crucial in determining cell-type specificity through differential expression and alternative splicing (AS) and that PRPF8 mutation causes global changes in splice site selection and exon inclusion that particularly affect genes involved in these cellular functions. This finding corroborates the hypothesis that retinal tissue identity is conferred by a specific splicing program and identifies retinal AS events as a framework toward the design of novel therapeutic opportunities.

Integration of Transformative Platforms for the Discovery of Causative Genes in Cardiovascular Diseases

Cardiovascular diseases are the leading cause of morbidity and mortality worldwide. Genome-wide association studies (GWAS) are powerful epidemiological tools to find genes and variants associated with cardiovascular diseases while follow-up biological studies allow to better understand the etiology and mechanisms of disease and assign causality. Improved methodologies and reduced costs have allowed wider use of bulk and single-cell RNA sequencing, human-induced pluripotent stem cells, organoids, metabolomics, epigenomics, and novel animal models in conjunction with GWAS. In this review, we feature recent advancements relevant to cardiovascular diseases arising from the integration of genetic findings with multiple enabling technologies within multidisciplinary teams to highlight the solidifying transformative potential of this approach. Well-designed workflows integrating different platforms are greatly improving and accelerating the unraveling and understanding of complex disease processes while promoting an effective way to find better drug targets, improve drug design and repurposing, and provide insight towards a more personalized clinical practice.

An extended catalogue of tandem alternative splice sites in human tissue transcriptomes

Tandem alternative splice sites (TASS) is a special class of alternative splicing events that are characterized by a close tandem arrangement of splice sites. Most TASS lack functional characterization and are believed to arise from splicing noise. Based on the RNA-seq data from the Genotype Tissue Expression project, we present an extended catalogue of TASS in healthy human tissues and analyze their tissue-specific expression. The expression of TASS is usually dominated by one major splice site (maSS), while the expression of minor splice sites (miSS) is at least an order of magnitude lower. Among 46k miSS with sufficient read support, 9k (20%) are significantly expressed above the expected noise level, and among them 2.5k are expressed tissue-specifically. We found significant correlations between tissue-specific expression of RNA-binding proteins (RBP), tissue-specific expression of miSS, and miSS response to RBP inactivation by shRNA. In combination with RBP profiling by eCLIP, this allowed prediction of novel cases of tissue-specific splicing regulation including a miSS in QKI mRNA that is likely regulated by PTBP1. The analysis of human primary cell transcriptomes suggested that both tissue-specific and cell-type-specific factors contribute to the regulation of miSS expression. More than 20% of tissue-specific miSS affect structured protein regions and may adjust protein-protein interactions or modify the stability of the protein core. The significantly expressed miSS evolve under the same selection pressure as maSS, while other miSS lack signatures of evolutionary selection and conservation. Using mixture models, we estimated that not more than 15% of maSS and not more than 54% of tissue-specific miSS are noisy, while the proportion of noisy splice sites among non-significantly expressed miSS is above 63%.

Pre-mRNA splicing is an important step in the processing of the genomic information during gene expression. During splicing, introns are excised from a gene transcript, and the remaining exons are ligated. Our work concerns one its particular subtype, which involves the so-called tandem alternative splice sites, a group of closely located exon borders that are used alternatively. We analyzed RNA-seq measurements of gene expression provided by the Genotype-Tissue Expression (GTEx) project, the largest to-date collection of such measurements in healthy human tissues, and constructed a detailed catalogue of tandem alternative splice sites. Within this catalogue, we characterized patterns of tissue-specific expression, regulation, impact on protein structure, and evolutionary selection acting on tandem alternative splice sites. In a number of genes, we predicted regulatory mechanisms that could be responsible for choosing one of many tandem alternative splice sites. The results of this study provide an invaluable resource for molecular biologists studying alternative splicing.

Transcriptomic analysis of the trade-off between endurance and burst-performance in the frog Xenopus allofraseri

BACKGROUND

Variation in locomotor capacity among animals often reflects adaptations to different environments. Despite evidence that physical performance is heritable, the molecular basis of locomotor performance and performance trade-offs remains poorly understood. In this study we identify the genes, signaling pathways, and regulatory processes possibly responsible for the trade-off between burst performance and endurance observed in Xenopus allofraseri, using a transcriptomic approach.

RESULTS

We obtained a total of about 121 million paired-end reads from Illumina RNA sequencing and analyzed 218,541 transcripts obtained from a de novo assembly. We identified 109 transcripts with a significant differential expression between endurant and burst performant individuals (FDR ≤ 0.05 and logFC ≥2), and blast searches resulted in 103 protein-coding genes. We found major differences between endurant and burst-performant individuals in the expression of genes involved in the polymerization and ATPase activity of actin filaments, cellular trafficking, proteoglycans and extracellular proteins secreted, lipid metabolism, mitochondrial activity and regulators of signaling cascades. Remarkably, we revealed transcript isoforms of key genes with functions in metabolism, apoptosis, nuclear export and as a transcriptional corepressor, expressed in either burst-performant or endurant individuals. Lastly, we find two up-regulated transcripts in burst-performant individuals that correspond to the expression of myosin-binding protein C fast-type (mybpc2). This suggests the presence of mybpc2 homoeologs and may have been favored by selection to permit fast and powerful locomotion.

CONCLUSION

These results suggest that the differential expression of genes belonging to the pathways of calcium signaling, endoplasmic reticulum stress responses and striated muscle contraction, in addition to the use of alternative splicing and effectors of cellular activity underlie locomotor performance trade-offs. Ultimately, our transcriptomic analysis offers new perspectives for future analyses of the role of single nucleotide variants, homoeology and alternative splicing in the evolution of locomotor performance trade-offs.

Prpf31 is essential for the survival and differentiation of retinal progenitor cells by modulating alternative splicing

Dysfunction of splicing factors often result in abnormal cell differentiation and apoptosis, especially in neural tissues. Mutations in pre-mRNAs processing factor 31 (PRPF31) cause autosomal dominant retinitis pigmentosa, a progressive retinal degeneration disease. The transcriptome-wide splicing events specifically regulated by PRPF31 and their biological roles in the development and maintenance of retina are still unclear. Here, we showed that the differentiation and viability of retinal progenitor cells (RPCs) are severely perturbed in prpf31 knockout zebrafish when compared with other tissues at an early embryonic stage. At the cellular level, significant mitotic arrest and DNA damage were observed. These defects could be rescued by the wild-type human PRPF31 rather than the disease-associated mutants. Further bioinformatic analysis and experimental verification uncovered that Prpf31 deletion predominantly causes the skipping of exons with a weak 5′ splicing site. Moreover, genes necessary for DNA repair and mitotic progression are most enriched among the differentially spliced events, which may explain the cellular and tissular defects in prpf31 mutant retinas. This is the first time that Prpf31 is demonstrated to be essential for the survival and differentiation of RPCs during retinal neurogenesis by specifically modulating the alternative splicing of genes involved in DNA repair and mitosis.

Consequences for Pancreatic β-Cell Identity and Function of Unregulated Transcript Processing

Mounting evidence suggests a role for alternative splicing (AS) of transcripts in the normal physiology and pathophysiology of the pancreatic β-cell. In the apparent absence of RNA repair systems, RNA decay pathways are likely to play an important role in controlling the stability, distribution and diversity of transcript isoforms in these cells. Around 35% of alternatively spliced transcripts in human cells contain premature termination codons (PTCs) and are targeted for degradation via nonsense-mediated decay (NMD), a vital quality control process. Inflammatory cytokines, whose levels are increased in both type 1 (T1D) and type 2 (T2D) diabetes, stimulate alternative splicing events and the expression of NMD components, and may or may not be associated with the activation of the NMD pathway. It is, however, now possible to infer that NMD plays a crucial role in regulating transcript processing in normal and stress conditions in pancreatic β-cells. In this review, we describe the possible role of Regulated Unproductive Splicing and Translation (RUST), a molecular mechanism embracing NMD activity in relationship to AS and translation of damaged transcript isoforms in these cells. This process substantially reduces the abundance of non-functional transcript isoforms, and its dysregulation may be involved in pancreatic β-cell failure in diabetes.

Splicing alterations in healthy aging and disease

Alternative RNA splicing is a key step in gene expression that allows generation of numerous messenger RNA transcripts encoding proteins of varied functions from the same gene. It is thus a rich source of proteomic and functional diversity. Alterations in alternative RNA splicing are observed both during healthy aging and in a number of human diseases, several of which display premature aging phenotypes or increased incidence with age. Age-associated splicing alterations include differential splicing of genes associated with hallmarks of aging, as well as changes in the levels of core spliceosomal genes and regulatory splicing factors. Here, we review the current known links between alternative RNA splicing, its regulators, healthy biological aging, and diseases associated with aging or aging-like phenotypes. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Processing > Splicing Regulation/Alternative Splicing.

Splicing alterations in healthy aging and disease

Alternative RNA splicing is a key step in gene expression that allows generation of numerous messenger RNA transcripts encoding proteins of varied functions from the same gene. It is thus a rich source of proteomic and functional diversity. Alterations in alternative RNA splicing are observed both during healthy aging and in a number of human diseases, several of which display premature aging phenotypes or increased incidence with age. Age-associated splicing alterations include differential splicing of genes associated with hallmarks of aging, as well as changes in the levels of core spliceosomal genes and regulatory splicing factors. Here, we review the current known links between alternative RNA splicing, its regulators, healthy biological aging, and diseases associated with aging or aging-like phenotypes. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Processing > Splicing Regulation/Alternative Splicing.

MeDAS: a Metazoan Developmental Alternative Splicing database

Alternative splicing is widespread throughout eukaryotic genomes and greatly increases transcriptomic diversity. Many alternative isoforms have functional roles in developmental processes and are precisely temporally regulated. To facilitate the study of alternative splicing in a developmental context, we created MeDAS, a Metazoan Developmental Alternative Splicing database. MeDAS is an added-value resource that re-analyses publicly archived RNA-seq libraries to provide quantitative data on alternative splicing events as they vary across the time course of development. It has broad temporal and taxonomic scope and is intended to assist the user in identifying trends in alternative splicing throughout development. To create MeDAS, we re-analysed a curated set of 2232 Illumina polyA+ RNA-seq libraries that chart detailed time courses of embryonic and post-natal development across 18 species with a taxonomic range spanning the major metazoan lineages from Caenorhabditis elegans to human. MeDAS is freely available at https://das.chenlulab.com both as raw data tables and as an interactive browser allowing searches by species, tissue, or genomic feature (gene, transcript or exon ID and sequence). Results will provide details on alternative splicing events identified for the queried feature and can be visualised at the gene-, transcript- and exon-level as time courses of expression and inclusion levels, respectively.

Global regulatory features of alternative splicing across tissues and within the nervous system of C. elegans

Alternative splicing plays a major role in shaping tissue-specific transcriptomes. Among the broad tissue types present in metazoans, the central nervous system contains some of the highest levels of alternative splicing. Although many documented examples of splicing differences between broad tissue types exist, there remains much to be understood about the splicing factors and the cis sequence elements controlling tissue and neuron subtype-specific splicing patterns. By using translating ribosome affinity purification coupled with deep-sequencing (TRAP-seq) in Caenorhabditis elegans, we have obtained high coverage profiles of ribosome-associated mRNA for three broad tissue classes (nervous system, muscle, and intestine) and two neuronal subtypes (dopaminergic and serotonergic neurons). We have identified hundreds of splice junctions that exhibit distinct splicing patterns between tissue types or within the nervous system. Alternative splicing events differentially regulated between tissues are more often frame-preserving, are more highly conserved across Caenorhabditis species, and are enriched in specific cis regulatory motifs, when compared with other types of exons. By using this information, we have identified a likely mechanism of splicing repression by the RNA-binding protein UNC-75/CELF via interactions with cis elements that overlap a 5′ splice site. Alternatively spliced exons also overlap more frequently with intrinsically disordered peptide regions than constitutive exons. Moreover, regulated exons are often shorter than constitutive exons but are flanked by longer intron sequences. Among these tissue-regulated exons are several highly conserved microexons <27 nt in length. Collectively, our results indicate a rich layer of tissue-specific gene regulation at the level of alternative splicing in C. elegans that parallels the evolutionary forces and constraints observed across metazoa.

Transcriptome analysis of alternative splicing in the pathogen life cycle in human foreskin fibroblasts infected with Trypanosoma cruzi

Trypanosoma cruzi is an intracellular protozoan parasite that causes Chagas disease as a zoonotic pathogen. The parasite has been shown to remodel expression in the host transcriptome under different conditions. Although alternative splicing (AS) is involved in virtually every biological function in eukaryotes, including cellular differentiation and responses to immune reactions, host AS events that occur as a result of T. cruzi infection have yet to be explored. In this study, we bioinformatically investigated the transcriptome AS dynamics of T. cruzi (Y strain) infected human foreskin fibroblasts using RNA-Seq data captured over four timepoints (4, 24, 48, and 72 h post infection (hpi)). We identified 1768, 399, 250, and 299 differentially expressed exons (AS exons) at 4, 24, 48, and 72 hpi, respectively, showing that host AS mechanism may have a significant role in the intracellular life cycle of the parasite. We present an exon skipping event in HDAC7, which is a candidate gene that is important in the parasite’s cell cycle. To sum up, this bioinformatics analysis of transcriptome may provide new potential insight into AS regulation in human foreskin fibroblast (HFF) cells infected by T. cruzi and into its implication to the parasite life cycle. Moreover, identified AS genes may provide new potential molecular candidates for improving treatment.

Real-time detection of mRNA splicing variants with specifically designed reverse-transcription loop-mediated isothermal amplification

Alternative splicing is a ubiquitous and crucial process in cellular processes and has a specific linkage with diseases. To date, developing cost-effective methods with high sensitivity and specificity for detection of splicing variants has been needed. Herein, we report a novel splicing variant assay based on specifically designed reverse-transcription loop-mediated isothermal amplification. After reverse transcribing the splicing variant into cDNA, four DNA primers are specifically designed to recognize six distinct regions. The four DNA primers can hybridize with corresponding sequences for extension and strand displacement DNA synthesis to form stem-loop DNA and then LAMP amplification is started. The proposed method can detect as low as 100 aM splicing variants in real-time fashion with high specificity, showing great potential in biological function and clinical studies.

A novel splicing variant assay is developed based on specifically designed reverse-transcription (RT) loop-mediated isothermal amplification.

Global Analysis of Alternative Splicing Difference in Peripheral Immune Organs between Tongcheng Pigs and Large White Pigs Artificially Infected with PRRSV In Vivo

Alternative splicing (AS) plays a significant role in regulating gene expression at the transcriptional level in eukaryotes. Flexibility and diversity of transcriptome and proteome can be significantly increased through alternative splicing of genes. In the present study, transcriptome data of peripheral immune organs including spleen and inguinal lymph nodes (ILN) were used to identify AS difference between PRRSV-resistant Tongcheng (TC) pigs and PRRSV-susceptible Large White (LW) pigs artificially infected with porcine reproductive and respiratory syndrome virus (PRRSV) in vivo. The results showed that PRRSV infection induced global alternative splicing events (ASEs) with different modes. Among them, 373 genes and 595 genes in the spleen and ILN of TC pigs, while 458 genes and 560 genes in the spleen and ILN of LW pigs had significantly differential ASEs. Alternative splicing was subject to tissue-specific and lineage-specific regulation in response to PRRSV infection. Enriched GO terms and pathways showed that genes with differential ASEs played important roles in transcriptional regulation, immune response, metabolism, and apoptosis. Furthermore, a splicing factor associated with apoptosis, SRSF4, was significantly upregulated in LW pigs. Functional analysis on apoptosis associated genes was validated by RT-PCR and DNA sequencing. These findings revealed different response to PRRSV between PRRSV-resistant TC pigs and PRRSV-susceptible LW pigs at the level of alternative splicing, suggesting the potential relationship between AS and disease resistance to PRRSV.

Kidney Is Essential for Blood Pressure Modulation by Dietary Potassium

Eating more potassium may reduce blood pressure and the occurrence of other cardiovascular diseases by actions on various systems, including the vasculature, the sympathetic nervous system, systemic metabolism, and body fluid volume. Among these, the kidney plays a major role in the potassium-rich diet-mediated blood pressure reduction. PURPOSE OF REVIEW: To provide an overview of recent discoveries about the mechanisms by which a potassium-rich diet leads to natriuresis. RECENT FINDINGS: Although the distal convoluted tubule (DCT) is a short part of the nephron that reabsorbs salt, via the sodium-chloride cotransporter (NCC), it is highly sensitive to changes in plasma potassium concentration. Activation or inhibition of NCC raises or lowers blood pressure. Recent work suggests that extracellular potassium concentration is sensed by the DCT via intracellular chloride concentration which regulates WNK kinases in the DCT. High-potassium diet targets NCC in the DCT, resulting in natriuresis and fluid volume reduction, which are protective from hypertension and other cardiovascular problems.

Alternative Splicing Enhances the Transcriptome Complexity of Liriodendron chinense

Alternative splicing (AS) plays pivotal roles in regulating plant growth and development, flowering, biological rhythms, signal transduction, and stress responses. However, no studies on AS have been performed in Liriodendron chinense, a deciduous tree species that has high economic and ecological value. In this study, we used multiple tools and algorithms to analyze transcriptome data derived from seven tissues via hybrid sequencing. Although only 17.56% (8,503/48,408) of genes in L. chinense were alternatively spliced, these AS genes occurred in 37,844 AS events. Among these events, intron retention was the most frequent AS event, producing 1,656 PTC-containing and 3,310 non-PTC-containing transcripts. Moreover, 183 long noncoding RNAs (lncRNAs) also underwent AS events. Furthermore, weighted gene coexpression network analysis (WGCNA) revealed that there were great differences in the activities of transcription and post-transcriptional regulation between pistils and leaves, and AS had an impact on many physiological and biochemical processes in L. chinense, such as photosynthesis, sphingolipid metabolism, fatty acid biosynthesis and metabolism. Moreover, our analysis showed that the features of genes may affect AS, as AS genes and non-AS genes had differences in the exon/intron length, transcript length, and number of exons/introns. In addition, the structure of AS genes may impact the frequencies and types of AS because AS genes with more exons or introns tended to exhibit more AS events, and shorter introns tended to be retained, whereas shorter exons tended to be skipped. Furthermore, eight AS genes were verified, and the results were consistent with our analysis. Overall, this study reveals that AS and gene interaction are mutual-on one hand, AS can affect gene expression and translation, while on the other hand, the structural characteristics of the gene can also affect AS. This work is the first to comprehensively report on AS in L. chinense, and it can provide a reference for further research on AS in L. chinense.

Tissue-specific mouse mRNA isoform networks

Alternative Splicing produces multiple mRNA isoforms of genes which have important diverse roles such as regulation of gene expression, human heritable diseases, and response to environmental stresses. However, little has been done to assign functions at the mRNA isoform level. Functional networks, where the interactions are quantified by their probability of being involved in the same biological process are typically generated at the gene level. We use a diverse array of tissue-specific RNA-seq datasets and sequence information to train random forest models that predict the functional networks. Since there is no mRNA isoform-level gold standard, we use single isoform genes co-annotated to Gene Ontology biological process annotations, Kyoto Encyclopedia of Genes and Genomes pathways, BioCyc pathways and protein-protein interactions as functionally related (positive pair). To generate the non-functional pairs (negative pair), we use the Gene Ontology annotations tagged with "NOT" qualifier. We describe 17 Tissue-spEcific mrNa iSoform functIOnal Networks (TENSION) following a leave-one-tissue-out strategy in addition to an organism level reference functional network for mouse. We validate our predictions by comparing its performance with previous methods, randomized positive and negative class labels, updated Gene Ontology annotations, and by literature evidence. We demonstrate the ability of our networks to reveal tissue-specific functional differences of the isoforms of the same genes. All scripts and data from TENSION are available at: https://doi.org/10.25380/iastate.c.4275191 .

Minor intron splicing revisited: identification of new minor intron-containing genes and tissue-dependent retention and alternative splicing of minor introns

BACKGROUND

Mutations in minor spliceosome components such as U12 snRNA (cerebellar ataxia) and U4atac snRNA (microcephalic osteodysplastic primordial dwarfism type 1 (MOPD1)) result in tissue-specific symptoms. Given that the minor spliceosome is ubiquitously expressed, we hypothesized that these restricted phenotypes might be caused by the tissue-specific regulation of the minor spliceosome targets, i.e. minor intron-containing genes (MIGs). The current model of inefficient splicing is thought to apply to the regulation of the ~ 500 MIGs identified in the U12DB. However this database was created more than 10 years ago. Therefore, we first wanted to revisit the classification of minor introns in light of the most recent reference genome. We then sought to address specificity of MIG expression, minor intron retention, and alternative splicing (AS) across mouse and human tissues.

RESULTS

We employed position-weight matrices to obtain a comprehensive updated list of minor introns, consisting of 722 mouse and 770 human minor introns. These can be found in the Minor Intron DataBase (MIDB). Besides identification of 99% of the minor introns found in the U12DB, we also discovered ~ 150 new MIGs. We then analyzed the RNAseq data from eleven different mouse tissues, which revealed tissue-specific MIG expression and minor intron retention. Additionally, many minor introns were efficiently spliced compared to their flanking major introns. Finally, we identified several novel AS events across minor introns in both mouse and human, which were also tissue-dependent. Bioinformatics analysis revealed that several of the AS events could result in the production of novel tissue-specific proteins. Moreover, like the major introns, we found that these AS events were more prevalent in long minor introns, while retention was favoured in shorter introns.

CONCLUSION

Here we show that minor intron splicing and AS across minor introns is a highly organised process that might be regulated in coordination with the major spliceosome in a tissue-specific manner. We have provided a framework to further study the impact of the minor spliceosome and the regulation of MIG expression. These findings may shed light on the mechanism underlying tissue-specific phenotypes in diseases associated with minor spliceosome inactivation. MIDB can be accessed at https://midb.pnb.uconn.edu .

Evolution of Alternative Splicing in Eudicots

Alternative pre-mRNA splicing (AS) is prevalent in plants and is involved in many interactions between plants and environmental stresses. However, the patterns and underlying mechanisms of AS evolution in plants remain unclear. By analyzing the transcriptomes of four eudicot species, we revealed that the divergence of AS is largely due to the gains and losses of AS events among orthologous genes. Furthermore, based on a subset of AS, in which AS can be directly associated with specific transcripts, we found that AS that generates transcripts containing premature termination codons (PTC), are likely more conserved than those that generate non-PTC containing transcripts. This suggests that AS coupled with nonsense-mediated decay (NMD) might play an important role in affecting mRNA levels post-transcriptionally. To understand the mechanisms underlying the divergence of AS, we analyzed the key determinants of AS using a machine learning approach. We found that the presence/absence of alternative splice site (SS) within the junction, the distance between the authentic SS and the nearest alternative SS, the size of exon-exon junctions were the major determinants for both alternative 5' donor site and 3' acceptor site among the studied species, suggesting a relatively conserved AS mechanism. The comparative analysis further demonstrated that variations of the identified AS determinants significantly contributed to the AS divergence among closely related species in both Solanaceae and Brassicaceae taxa. Together, these results provide detailed insights into the evolution of AS in plants.

Expression of NOTCH3 exon 16 differentiates Diffuse Large B-cell Lymphoma into molecular subtypes and is associated with prognosis

Diffuse large B-cell lymphoma (DLBCL) is a heterogeneous disease with diverse clinical presentation and outcome. Bio-clinical prognostic models including oncogene expression and cell-of-origin phenotyping has been developed, however, approximately 30% of all patients still die from their disease, illustrating the need for additional prognostic biomarkers associating oncogenesis and phenotypic subclasses. Hence, we tested if alternative splice variations have biomarker potential. Initial alternative splicing analysis of human exon array from clinical DLBCL samples identified candidate genes. Experimental validation by ddPCR was performed in a DLBCL cohort classified into ABC/GCB subclasses, B-cell associated gene signatures (BAGS: naive, centroblast, centrocyte, memory, and plasmablast), and vincristine resistant gene signatures. Prognostic potential was assessed for aberrantly spliced transcripts. Thus, NOTCH3 was identified as alternatively spliced, with differential exon 16 depletion (−exon 16) between differentiation associated BAGS subtypes. Predicted vincristine resistant patients of the GCB subclass had significantly downregulated NOTCH3 −exon 16 transcript expression and tended to display adverse overall survival for R-CHOP treated patients. In conclusion, we have identified a specific alternatively spliced NOTCH3 event that differentiate molecular subtypes of DLBCL and display prognostic and predictive biomarker potential in GCB DLBCL.

Prader-Willi locus Snord116 RNA processing requires an active endogenous allele and neuron-specific splicing by Rbfox3/NeuN

Prader-Willi syndrome (PWS), an imprinted neurodevelopmental disorder characterized by metabolic, sleep and neuropsychiatric features, is caused by the loss of paternal SNORD116, containing only non-coding RNAs (ncRNAs). The primary SNORD116 transcript is processed into small nucleolar RNAs (snoRNAs), which localize to nucleoli, and their spliced host gene 116HG, which is retained at its site of transcription. While functional complementation of the SNORD116 ncRNAs is a desirable goal for treating PWS, the mechanistic requirements of SNORD116 RNA processing are poorly understood. Here we developed and tested a novel transgenic mouse which ubiquitously expresses Snord116 on both a wild-type and a Snord116 paternal deletion (Snord116+/-) background. Interestingly, while the Snord116 transgene was ubiquitously expressed in multiple tissues, splicing of the transgene and production of snoRNAs was limited to brain tissues. Knockdown of Rbfox3, encoding neuron-specific splicing factor neuronal nuclei (NeuN) in Snord116+/--derived neurons, reduced splicing of the transgene in neurons. RNA fluorescence in situ hybridization for 116HG revealed a single significantly larger signal in transgenic mice, demonstrating colocalization of transgenic and endogenous 116HG RNAs. Similarly, significantly increased snoRNA levels were detected in transgenic neuronal nucleoli, indicating that transgenic Snord116 snoRNAs were effectively processed and localized. In contrast, neither transgenic 116HG nor snoRNAs were detectable in either non-neuronal tissues or Snord116+/- neurons. Together, these results demonstrate that exogenous expression and neuron-specific splicing of the Snord116 locus are insufficient to rescue the genetic deficiency of Snord116 paternal deletion. Elucidating the mechanisms regulating Snord116 processing and localization is essential to develop effective gene replacement therapies for PWS.

Phylotranscriptomics suggests the jawed vertebrate ancestor could generate diverse helper and regulatory T cell subsets

BACKGROUND

The cartilaginous fishes diverged from other jawed vertebrates ~ 450 million years ago (mya). Despite this key evolutionary position, the only high-quality cartilaginous fish genome available is for the elephant shark (Callorhinchus milii), a chimaera whose ancestors split from the elasmobranch lineage ~ 420 mya. Initial analysis of this resource led to proposals that key components of the cartilaginous fish adaptive immune system, most notably their array of T cell subsets, was primitive compared to mammals. This proposal is at odds with the robust, antigen-specific antibody responses reported in elasmobranchs following immunization. To explore this discrepancy, we generated a multi-tissue transcriptome for small-spotted catshark (Scyliorhinus canicula), a tractable elasmobranch model for functional studies. We searched this, and other newly available sequence datasets, for CD4+ T cell subset-defining genes, aiming to confirm the presence or absence of each subset in cartilaginous fishes.

RESULTS

We generated a new transcriptome based on a normalised, multi-tissue RNA pool, aiming to maximise representation of tissue-specific and lowly expressed genes. We utilized multiple transcriptomic datasets and assembly variants in phylogenetic reconstructions to unambiguously identify several T cell subset-specific molecules in cartilaginous fishes for the first time, including interleukins, interleukin receptors, and key transcription factors. Our results reveal the inability of standard phylogenetic reconstruction approaches to capture the site-specific evolutionary processes of fast-evolving immune genes but show that site-heterogeneous mixture models can adequately do so.

CONCLUSIONS

Our analyses reveal that cartilaginous fishes are capable of producing a range of CD4+ T cell subsets comparable to that of mammals. Further, that the key molecules required for the differentiation and functioning of these subsets existed in the jawed vertebrate ancestor. Additionally, we highlight the importance of considering phylogenetic diversity and, where possible, utilizing multiple datasets for individual species, to accurately infer gene presence or absence at higher taxonomic levels.

Identification of novel transcripts and peptides in developing murine lens

We previously investigated the transcriptome and proteome profiles of the murine ocular lens at six developmental time points including two embryonic (E15 and E18) and four postnatal time points (P0, P3, P6, and P9). Here, we extend our analyses to identify novel transcripts and peptides in developing  mouse lens. We identified a total of 9,707 novel transcripts and 325 novel fusion genes in developing mouse lens. Additionally, we identified 13,281 novel alternative splicing (AS) events in mouse lens including 6,990 exon skipping (ES), 2,447 alternative 3' splice site (A3SS), 1,900 alternative 5' splice site (A5SS), 1,771 mutually exclusive exons (MXE), and 173 intron retention (IR). Finally, we integrated our OMIC (Transcriptome and Proteome) datasets identifying 20 novel peptides in mouse lens. All 20 peptides were validated through matching MS/MS spectra of synthetic peptides. To the best of our knowledge, this is the first report integrating OMIC datasets to identify novel peptides in developing murine lens.

Comprehensive map of age-associated splicing changes across human tissues and their contributions to age-associated diseases

Alternative splicing contributes to phenotypic diversity at multiple biological scales, and its dysregulation is implicated in both ageing and age-associated diseases in human. Cross-tissue variability in splicing further complicates its links to age-associated phenotypes and elucidating these links requires a comprehensive map of age-associated splicing changes across multiple tissues. Here, we generate such a map by analyzing ~8500 RNA-seq samples across 48 tissues in 544 individuals. Employing a stringent model controlling for multiple confounders, we identify 49,869 tissue-specific age-associated splicing events of 7 distinct types. We find that genome-wide splicing profile is a better predictor of biological age than the gene and transcript expression profiles, and furthermore, age-associated splicing provides additional independent contribution to age-associated complex diseases. We show that the age-associated splicing changes may be explained, in part, by concomitant age-associated changes of the upstream splicing factors. Finally, we show that our splicing-based model of age can successfully predict the relative ages of cells in 8 of the 10 paired longitudinal data as well as in 2 sets of cell passage data. Our study presents the first systematic investigation of age-associated splicing changes across tissues, and further strengthening the links between age-associated splicing and age-associated diseases.