An atlas of alternative splicing profiles and functional associations reveals new regulatory programs and genes that simultaneously express multiple major isoforms

Plourde M, Joly DL, Germain H. Comprehensive mapping of Arabidopsis alternative splicing landscape reveals key insights into plant development and immunity

The different steps of alternative splicing (AS) in plants and its regulatory mechanisms have already been studied extensively. Its broader impact on cell identity, plant immunity-related genes, and their study as a whole remains to be investigated. Using transgenic plants, we sorted 11 different Arabidopsis thaliana cell types ranging from root to aerial organs using fluorescence-activated cell sorting. RNA-seq data were analyzed with vast-tools and enabled us to generate a high-resolution AS landscape across multiple cell types, all collected through the same experimental procedure. The analysis of cell type-specific gene expression and alternative splicing events highlights the importance of AS on transcription and AS regulation itself. AS is also shown to be tightly linked to cell identity. By using closely related cell types, we captured alternative splicing events involved in specific stages of plant development. The columella cells, among others, show intensified AS regulation and an interesting splicing profile, especially regarding immunity-related genes. Overall, our analysis brings a valuable tool in the study of cell type identity, plant immunity, and AS.

Mutations in the microexon splicing regulator srrm4 have minor phenotypic effects on zebrafish neural development

Achieving a diversity of neuronal cell types and circuits during brain development requires alternative splicing of developmentally regulated mRNA transcripts. Microexons are a type of alternatively spliced exon that are 3-27 nucleotides in length and are predominantly expressed in neuronal tissues. A key regulator of microexon splicing is the RNA-binding protein Serine/arginine repetitive matrix 4 (Srrm4). Srrm4 is a highly conserved, vertebrate splicing factor that is part of an ancient family of splicing proteins. To better understand the function of Srrm4 during brain development, we examined the neural expression of zebrafish srrm4 from 1 to 5 days of development using fluorescence in situ hybridization. We found that srrm4 has a dynamically changing expression pattern, with expression in diverse cell types and stages during development. We then used CRISPR-based mutagenesis to generate zebrafish srrm4 mutants. Unlike previously described morphant phenotypes, srrm4 mutants did not show overt morphological defects. Whole-brain morphometric analysis revealed a reduction in optic tectum neuropil in G0 crispants that, unexpectedly, was also not replicated in stable mutants. Sequencing of wild-type and mutant transcriptomes revealed only minor changes in splicing and did not support a hypothesis of transcriptional adaptation, suggesting that another, as yet, unidentified mechanism of compensation is occurring. srrm4 thus appears to have a limited role in zebrafish neural development.

A highly conserved neuronal microexon in DAAM1 controls actin dynamics, RHOA/ROCK signaling, and memory formation

Actin cytoskeleton dynamics is essential for proper nervous system development and function. A conserved set of neuronal-specific microexons influences multiple aspects of neurobiology; however, their roles in regulating the actin cytoskeleton are unknown. Here, we study a microexon in DAAM1, a formin-homology-2 (FH2) domain protein involved in actin reorganization. Microexon inclusion extends the linker region of the DAAM1 FH2 domain, altering actin polymerization. Genomic deletion of the microexon leads to neuritogenesis defects and increased calcium influx in differentiated neurons. Mice with this deletion exhibit postsynaptic defects, fewer immature dendritic spines, impaired long-term potentiation, and deficits in memory formation. These phenotypes are associated with increased RHOA/ROCK signaling, which regulates actin-cytoskeleton dynamics, and are partially rescued by treatment with a ROCK inhibitor. This study highlights the role of a conserved neuronal microexon in regulating actin dynamics and cognitive functioning.

Gene Duplication and Alternative Splicing as Evolutionary Drivers of Proteome Specialization

Animals comprise hundreds of cell types, each with specialized biological functions. However, many genes expressed in each cell type belong to widely conserved gene families with ancestrally ubiquitous expression. This raises a paradox: how have these genes evolved to shape cell type-specific traits without compromising their ancestral function in all other cells? This can be achieved through gene duplication and the origin of regulated, alternatively spliced exons, which generate new related proteins in the form of paralogous genes and alternative isoforms, respectively. Here, we explore how such new related proteins can contribute to the evolution of specific cell types while preserving broader ancestral roles. Specifically, we separately classify possible expression and functional fates for new related proteins and discuss their interplays and evolutionary likelihood. Our primary hypothesis is that expression specialization, mostly coupled with functional specialization, is the predominant fate for both paralogous genes and alternative isoforms throughout animal evolution.

TAF2 condensation in nuclear speckles links basal transcription factor TFIID to RNA splicing factors

TFIID is an essential basal transcription factor, crucial for RNA polymerase II (pol II) promoter recognition and transcription initiation. The TFIID complex consists of the TATA binding protein (TBP) and 13 TBP-associated factors (TAFs) that contain intrinsically disordered regions (IDRs) with currently unknown functions. Here, we show that a conserved IDR drives TAF2 to nuclear speckle condensates independently of other TFIID subunits. Quantitative mass spectrometry analyses reveal TAF2 proximity to RNA splicing factors including specific interactions of the TAF2 IDR with SRRM2 in nuclear speckles. Deleting the IDR from TAF2 does not majorly impact global gene expression but results in changes of alternative splicing events. Further, genome-wide binding analyses suggest that the TAF2 IDR impedes TAF2 promoter association by guiding TAF2 to nuclear speckles. This study demonstrates that an IDR within the large multiprotein complex TFIID controls nuclear compartmentalization and thus links distinct molecular processes, namely transcription initiation and RNA splicing.

Functional impact of splicing variants in the elaboration of complex traits in cattle

GWAS conducted directly on imputed whole genome sequence have led to the identification of numerous genetic variants associated with agronomic traits in cattle. However, such variants are often simply markers in linkage disequilibrium with the actual causal variants, which is a limiting factor for the development of accurate genomic predictions. It is possible to identify causal variants by integrating information on how variants impact gene expression into GWAS output. RNA splicing plays a major role in regulating gene expression. Thus, assessing the effect of variants on RNA splicing may explain their function. Here, we use a high-throughput strategy to functionally analyse putative splice-disrupting variants in the bovine genome. Using GWAS, massively parallel reporter assay and deep learning algorithms designed to predict splice-disrupting variants, we identify 38 splice-disrupting variants associated with complex traits in cattle, three of which could be classified as causal. Our results indicate that splice-disrupting variants are widely found in the quantitative trait loci related to these phenotypes. Using our combined approach, we also assess the validity of splicing predictors originally developed to analyse human variants in the context of the bovine genome.

The splicing factor Acin1 is essential for embryonic development but has limited effects on muscle structure and homeostasis

Apoptotic chromatin condensation inducer 1 (Acin1) is an RNA-binding protein involved in the regulation of alternative splicing, but its physiological function remains unclear. Global deletion of Acin1 causes embryonic lethality around E11.5, with mutants exhibiting developmental delays and increased apoptosis. Using conditional knockout mice, we found that skeletal muscle myofiber-specific Acin1 knockout mice (Acin1 MKO) are viable and fertile and that Acin1 MKO mice show enlarged myofibers and ongoing muscle damage and regeneration, characterized by increased central nuclei and embryonic myosin heavy chain expression. RNA-seq analysis revealed that Acin1 deletion altered the expression and splicing patterns of genes crucial for muscle function. Notable changes included modified splicing of genes associated with muscle disease and mitochondrial function, often resulting in the expression of gene variants typical of immature or diseased muscle. These findings suggest that Acin1 is essential for embryonic development and has limited effects on muscle structure and homeostasis via its regulation of gene expression and alternative splicing.

PRMT5 promotes full-length HTT expression by repressing multiple proximal intronic polyadenylation sites

Expansion of the CAG trinucleotide repeat tract in exon 1 of the Huntingtin (HTT) gene causes Huntington's disease (HD) through the expression of a polyglutamine-expanded form of the HTT protein. This mutation triggers cellular and biochemical pathologies, leading to cognitive, motor, and psychiatric symptoms in HD patients. Targeting HTT splicing with small molecule drugs is a compelling approach to lowering HTT protein levels to treat HD, and splice modulators are currently being tested in the clinic. Here, we identify PRMT5 as a novel regulator of HTT messenger RNA (mRNA) splicing and alternative polyadenylation. PRMT5 inhibition disrupts the splicing of HTT introns 9 and 10, leading to the activation of multiple proximal intronic polyadenylation sites within these introns and promoting premature termination, cleavage, and polyadenylation of the HTT mRNA. This suggests that HTT protein levels may be lowered due to this mechanism. We also detected increasing levels of these truncated HTT transcripts across a series of neuronal differentiation samples, which correlated with lower PRMT5 expression. Notably, PRMT5 inhibition in glioblastoma stem cells potently induced neuronal differentiation. We posit that PRMT5-mediated regulation of intronic polyadenylation, premature termination, and cleavage of the HTT mRNA modulates HTT expression and plays an important role during neuronal differentiation.

Circular RNAs exhibit exceptional stability in the aging brain and serve as reliable age and experience indicators

Circular RNAs (circRNAs) increase in the brain with age across various animal systems. To elucidate the reasons behind this phenomenon, we profile circRNAs from fly heads at six time points throughout their lifespan. Our results reveal a linear increase in circRNA levels with age, independent of changes in mRNA levels, overall transcription, intron retention, or host gene splicing, demonstrating that the age-related accumulation is due to high stability rather than increased biogenesis. This remarkable stability suggests that circRNAs can serve as markers of environmental experience. Indeed, flies exposed to a 10-day regimen at 29°C exhibit higher levels of specific circRNAs even 6 weeks after returning to standard conditions, indicating that circRNAs can reveal past environmental stimuli. Moreover, half-life measurements show circRNA stability exceeding 20 days, with some displaying virtually no degradation. These findings underscore the remarkable stability of circRNAs in vivo and their potential as markers for stress and life experiences.

Biosurfer for systematic tracking of regulatory mechanisms leading to protein isoform diversity

Long-read RNA-seq has shed light on transcriptomic complexity, but questions remain about the functionality of downstream protein products. We introduce Biosurfer, a computational approach for comparing protein isoforms, while systematically tracking the transcriptional, splicing, and translational variations that underlie differences in the sequences of the protein products. Using Biosurfer, we analyzed the differences in 35,082 pairs of GENCODE annotated protein isoforms, finding a majority (70%) of variable N-termini are due to the alternative transcription start sites, while only 9% arise from 5' UTR alternative splicing (AS). Biosurfer's detailed tracking of nucleotide-to-residue relationships helps reveal an uncommonly tracked source of single amino acid residue changes arising from the codon splits at junctions. For 17% of internal sequence changes, such split codon patterns lead to single residue differences, termed "ragged codons." Of variable C-termini, 72% involve splice- or intron retention-induced reading frameshifts. We systematically characterize an unusual pattern of reading frame changes, in which the first frameshift is closely followed by a distinct second frameshift that restores the original frame, which we term a "snapback" frameshift. We analyze the long-read RNA-seq-predicted proteome of a human cell line and find similar trends as compared to our GENCODE analysis, with the exception of a higher proportion of transcripts predicted to undergo nonsense-mediated decay. Biosurfer's comprehensive characterization of long-read RNA-seq data sets should accelerate insights of the functional role of protein isoforms, providing mechanistic explanation of the origins of the proteomic diversity driven by the AS. Biosurfer is available as a Python package.

Drosophila ovarian stem cell niche ageing involves coordinated changes in transcription and alternative splicing

Gene expression (GE) and alternative splicing (AS) contribute to the formation of new interaction networks with potentially significant cellular functions. Here, we investigate ageing in the Drosophila female germline stem cell (GSC) niche and describe functional changes in both GE and AS. The GSC niche comprises three types of support cells, whose ageing transcriptomes reveal differential GE and AS variations related to cell adhesion, cytoskeleton and neural signalling. Because each population show distinctive GE and AS changes, niche cell types possess unique ageing signatures. Depending on the cell population, groups of genes display changes in both GE and AS, revealing a coordinated regulation of transcription and splicing during niche ageing. One such gene is Fasciclin 2, a neural adhesion molecule that we find is essential for niche functioning. Furthermore, genes involved in AS undergo changes in GE and/or AS themselves, providing a mechanistic explanation for the coordination of these two processes during niche ageing. This is the case of the splicing factor Smu1, described here as a key element necessary for ovarian niche homeostasis.

The main sources of molecular organization in the cell. Atlas of self-organized and self-regulated dynamic biostructures

One of the most important goals of contemporary biology is to understand the principles of the molecular order underlying the complex dynamic architecture of cells. Here, we present an overview of the main driving forces involved in the cellular molecular complexity and in the emergent functional dynamic structures, spanning from the most basic molecular organization levels to the complex emergent integrative systemic behaviors. First, we address the molecular information processing which is essential in many complex fundamental mechanisms such as the epigenetic memory, alternative splicing, regulation of transcriptional system, and the adequate self-regulatory adaptation to the extracellular environment. Next, we approach the biochemical self-organization, which is central to understand the emergency of metabolic rhythms, circadian oscillations, and spatial traveling waves. Such a complex behavior is also fundamental to understand the temporal compartmentalization of the cellular metabolism and the dynamic regulation of many physiological activities. Numerous examples of biochemical self-organization are considered here, which show that practically all the main physiological processes in the cell exhibit this type of dynamic molecular organization. Finally, we focus on the biochemical self-assembly which, at a primary level of organization, is a basic but important mechanism for the order in the cell allowing biomolecules in a disorganized state to form complex aggregates necessary for a plethora of essential structures and physiological functions. In total, more than 500 references have been compiled in this review. Due to these main sources of order, systemic functional structures emerge in the cell, driving the metabolic functionality towards the biological complexity.

From computational models of the splicing code to regulatory mechanisms and therapeutic implications

Since the discovery of RNA splicing and its role in gene expression, researchers have sought a set of rules, an algorithm or a computational model that could predict the splice isoforms, and their frequencies, produced from any transcribed gene in a specific cellular context. Over the past 30 years, these models have evolved from simple position weight matrices to deep-learning models capable of integrating sequence data across vast genomic distances. Most recently, new model architectures are moving the field closer to context-specific alternative splicing predictions, and advances in sequencing technologies are expanding the type of data that can be used to inform and interpret such models. Together, these developments are driving improved understanding of splicing regulatory mechanisms and emerging applications of the splicing code to the rational design of RNA- and splicing-based therapeutics.

Shiba: a versatile computational method for systematic identification of differential RNA splicing across platforms

Alternative pre-mRNA splicing (AS) is a fundamental regulatory process that generates transcript diversity and cell type variation. We developed Shiba, a comprehensive method that integrates transcript assembly, splicing event identification, read counting, and differential splicing analysis across RNA-seq platforms. Shiba excels in capturing annotated and unannotated AS events with superior accuracy, sensitivity, and reproducibility. It addresses the often-overlooked issue of junction read imbalance, significantly reducing false positives to aid target prioritization and downstream analyses. Unlike other tools that require large numbers of biological replicates or resulting in low sensitivity and high false positives, Shiba's statistics framework is agnostic to sample size, as demonstrated by simulated data and its effective application to real n= 1 RNA-seq datasets. To extend its utility to single-cell RNA-seq, we developed scShiba, which applies Shiba's pseudobulk approach to analyze splicing at the cluster level. scShiba successfully revealed AS regulation in developmental dopaminergic neurons and differences between excitatory and inhibitory neurons. Both Shiba and scShiba are available in Docker/Singularity containers and Snakemake pipelines, ensuring reproducibility. With their comprehensive capabilities, Shiba and scShiba enable systematic quantification of alternative splicing events across various platforms, laying a solid foundation for mechanistic exploration of the functional complexity in RNA splicing.

Toward a comprehensive profiling of alternative splicing proteoform structures, interactions and functions

The mRNA splicing machinery has been estimated to generate 100,000 known protein-coding transcripts for 20,000 human genes (Ensembl, Sept. 2024). However, this set is expanding with the massive and rapidly growing data coming from high-throughput technologies, particularly single-cell and long-read sequencing. Yet, the implications of splicing complexity at the protein level remain largely uncharted. In this review, we describe the current advances toward systematically assessing the contribution of alternative splicing to proteome function diversification. We discuss the potential and challenges of using artificial intelligence-based techniques in identifying alternative splicing proteoforms and characterising their structures, interactions, and functions.

Loss of cped1 does not affect bone and lean tissue in zebrafish

Human genetic studies have nominated cadherin-like and PC-esterase domain-containing 1 (CPED1) as a candidate target gene mediating bone mineral density (BMD) and fracture risk heritability. Recent efforts to define the role of CPED1 in bone in mouse and human models have revealed complex alternative splicing and inconsistent results arising from gene targeting, making its function in bone difficult to interpret. To better understand the role of CPED1 in adult bone mass and morphology, we conducted a comprehensive genetic and phenotypic analysis of cped1 in zebrafish, an emerging model for bone and mineral research. We analyzed two different cped1 mutant lines and performed deep phenotyping to characterize more than 200 measures of adult vertebral, craniofacial, and lean tissue morphology. We also examined alternative splicing of zebrafish cped1 and gene expression in various cell/tissue types. Our studies fail to support an essential role of cped1 in adult zebrafish bone. Specifically, homozygous mutants for both cped1 mutant alleles, which are expected to result in loss-of-function and impact all cped1 isoforms, exhibited no significant differences in the measures examined when compared to their respective wildtype controls, suggesting that cped1 does not significantly contribute to these traits. We identified sequence differences in critical residues of the catalytic triad between the zebrafish and mouse orthologs of CPED1, suggesting that differences in key residues, as well as distinct alternative splicing, could underlie different functions of CPED1 orthologs in the two species. Our studies fail to support a requirement of cped1 in zebrafish bone and lean tissue, adding to evidence that variants at 7q31.31 can act independently of CPED1 to influence BMD and fracture risk.

Exon inclusion signatures enable accurate estimation of splicing factor activity

Splicing factors control exon inclusion in messenger RNA, shaping transcriptome and proteome diversity. Their catalytic activity is regulated by multiple layers, making single-omic measurements on their own fall short in identifying which splicing factors underlie a phenotype. Here, we propose splicing factor activity can be estimated by interpreting changes in exon inclusion. We benchmark methods to construct splicing factor→exon networks and calculate activity. Combining RNA-seq perturbation-based networks with VIPER (virtual inference of protein activity by enriched regulon analysis) accurately captures splicing factor activation modulated by different regulatory layers. This approach consolidates splicing factor regulation into a single score derived solely from exon inclusion signatures, allowing functional interpretation of heterogeneous conditions. As a proof of concept, we identify recurrent cancer splicing programs, revealing oncogenic- and tumor suppressor-like splicing factors missed by conventional methods. These programs correlate with patient survival and key cancer hallmarks: initiation, proliferation, and immune evasion. Altogether, we show splicing factor activity can be accurately estimated from exon inclusion changes, enabling comprehensive analyses of splicing regulation with minimal data requirements.

Coordinated neuron-specific splicing events restrict nucleosome engagement of the LSD1 histone demethylase complex

Chromatin regulatory proteins are expressed broadly and assumed to exert the same intrinsic function across cell types. Here, we report that 14 chromatin regulators undergo evolutionary-conserved neuron-specific splicing events involving microexons. Among them are two components of a histone demethylase complex: LSD1 H3K4 demethylase and the H3K4me0-reader PHF21A. We found that neuronal LSD1 splicing reduces the enzymes' affinity to the nucleosome. Meanwhile, neuronal PHF21A splicing significantly attenuates histone H3 binding and further ablates the DNA-binding function exerted by an AT-hook motif. Furthermore, in vitro reconstitution of the canonical and neuronal PHF21A-LSD1 complexes, combined with in vivo methylation mapping, identified the neuronal complex as a hypomorphic H3K4 demethylating machinery. The neuronal PHF21A, albeit with its weaker nucleosome binding, is necessary for normal gene expression and the H3K4 landscape in the developing brain. Thus, ubiquitously expressed chromatin regulatory complexes can exert neuron-specific functions via alternative splicing of their subunits.

Shiba: A versatile computational method for systematic identification of differential RNA splicing across platforms

Alternative pre-mRNA splicing (AS) is a fundamental regulatory process that generates transcript diversity and cell type variation. We developed Shiba, a comprehensive method that integrates transcript assembly, splicing event identification, read counting, and differential splicing analysis across RNA-seq platforms. Shiba excels in capturing annotated and unannotated AS events with superior accuracy, sensitivity, and reproducibility. It addresses the often-overlooked issue of junction read imbalance, significantly reducing false positives to aid target prioritization and downstream analyses. Unlike other tools that require large numbers of biological replicates or resulting in low sensitivity and high false positives, Shiba's statistics framework is agnostic to sample size, as demonstrated by simulated data and its effective application to real n=1 RNA-seq datasets. To extend its utility to single-cell RNA-seq, we developed scShiba, which applies Shiba's pseudobulk approach to analyze splicing at the cluster level. scShiba successfully revealed AS regulation in developmental dopaminergic neurons and differences between excitatory and inhibitory neurons. Both Shiba and scShiba are available in Docker/Singularity containers and Snakemake pipelines, ensuring reproducibility. With their comprehensive capabilities, Shiba and scShiba enable systematic quantification of alternative splicing events across various platforms, laying a solid foundation for mechanistic exploration of the functional complexity in RNA splicing.

Proteomic Profiling Towards a Better Understanding of Genetic Based Muscular Diseases: The Current Picture and a Look to the Future

Proteomics accelerates diagnosis and research of muscular diseases by enabling the robust analysis of proteins relevant for the manifestation of neuromuscular diseases in the following aspects: (i) evaluation of the effect of genetic variants on the corresponding protein, (ii) prediction of the underlying genetic defect based on the proteomic signature of muscle biopsies, (iii) analysis of pathophysiologies underlying different entities of muscular diseases, key for the definition of new intervention concepts, and (iv) patient stratification according to biochemical fingerprints as well as (v) monitoring the success of therapeutic interventions. This review presents-also through exemplary case studies-the various advantages of mass proteomics in the investigation of genetic muscle diseases, discusses technical limitations, and provides an outlook on possible future application concepts. Hence, proteomics is an excellent large-scale analytical tool for the diagnostic workup of (hereditary) muscle diseases and warrants systematic profiling of underlying pathophysiological processes. The steady development may allow to overcome existing limitations including a quenched dynamic range and quantification of different protein isoforms. Future directions may include targeted proteomics in diagnostic settings using not only muscle biopsies but also liquid biopsies to address the need for minimally invasive procedures.

Pol II degradation activates cell death independently from the loss of transcription

Pol II-mediated transcription is essential for eukaryotic life. While loss of transcription is thought to be universally lethal, the associated mechanisms promoting cell death are not yet known. Here, we show that death following loss of Pol II is not caused by dysregulated gene expression. Instead, death occurs in response to the loss of Pol II protein itself, specifically loss of the enzymatic subunit, Rbp1. Loss of Pol II exclusively activates apoptosis, and expression of a transcriptionally inactive version of Rpb1 rescues cell viability. Using functional genomics, we identify a previously uncharacterized mechanism that regulates lethality following loss of Pol II, which we call the Pol II Degradation-dependent Apoptotic Response (PDAR). Using the genetic dependencies of PDAR, we identify clinically used drugs that owe their efficacy to a PDAR-dependent mechanism. Our findings unveil a novel apoptotic signaling response that contributes to the efficacy of a wide array of anti-cancer therapies.

Nucleo-cytoplasmic distribution of SAP18 reveals its dual function in splicing regulation and heat-stress response in Arabidopsis

Dynamic shuttling of proteins between the nucleus and cytoplasm orchestrates vital functions in eukaryotes. Here, we reveal the multifaceted functions of Arabidopsis Sin3-associated protein 18 kDa (SAP18) in the regulation of development and heat-stress tolerance. Proteomic analysis demonstrated that SAP18 is a core component of the nuclear apoptosis- and splicing-associated protein (ASAP) complex in Arabidopsis, contributing to the precise splicing of genes associated with leaf development. Genetic analysis further confirmed the critical role of SAP18 in different developmental processes as part of the ASAP complex, including leaf morphogenesis and flowering time. Interestingly, upon heat shock, SAP18 translocates from the nucleus to cytoplasmic stress granules and processing bodies. The heat-sensitive phenotype of a SAP18 loss-of-function mutant revealed a novel role for SAP18 in plant thermoprotection. These findings significantly expand our understanding of the relevance of SAP18 for plant growth, linking nuclear splicing with cytoplasmic stress responses and providing new perspectives for future exploration of plant thermotolerance mechanisms.

Phosphorylation of a nuclear condensate regulates cohesion and mRNA retention

Nuclear speckles are membraneless organelles that associate with active transcription sites and participate in post-transcriptional mRNA processing. During the cell cycle, nuclear speckles dissolve following phosphorylation of their protein components. Here, we identify the PP1 family as the phosphatases that counteract kinase-mediated dissolution. PP1 overexpression increases speckle cohesion and leads to retention of mRNA within speckles and the nucleus. Using APEX2 proximity labeling combined with RNA-sequencing, we characterize the recruitment of specific RNAs. We find that many transcripts are preferentially enriched within nuclear speckles compared to the nucleoplasm, particularly chromatin- and nucleus-associated transcripts. While total polyadenylated RNA retention increases with nuclear speckle cohesion, the ratios of most mRNA species to each other are constant, indicating non-selective retention. We further find that cellular responses to heat shock, oxidative stress, and hypoxia include changes to the phosphorylation and cohesion of nuclear speckles and to mRNA retention. Our results demonstrate that tuning the material properties of nuclear speckles provides a mechanism for the acute control of mRNA localization.

CDC40 suppression induces CDCA5 splicing defects and anti-proliferative effects in lung cancer cells

High mortality and low response rates in lung cancer patients call for novel therapeutic targets. Data mining of whole-genome genetic dependency screens suggest Cell Division Cycle 40 (CDC40) to be an essential protein for lung cancer cell survival. We characterized CDC40 knockdown effects in multiple lung cancer cell lines, revealing induced cell cycle defects that resulted in strong growth inhibition and activation of apoptosis. Global transcriptional and splicing changes were also investigated, where CDC40 knockdown resulted in perturbation of splicing- and translation-related genes as well as more transcripts with intron retention. In the transcript of the cell cycle regulatory protein CDCA5, CDC40 knockdown was shown to induce retention of the first intron, leading to an increase in the unspliced CDCA5 transcript and subsequent decrease in CDCA5 protein expression. Additionally, protein-protein interactions of CDC40 were explored and spliceosome components were found to be its main binding partners, further highlighting the role of CDC40 in splicing. CDC40 mutation analysis suggests that it may be difficult to disrupt key interactions using small molecules within a large complex. Our results demonstrate that CDC40 is essential for lung cancer cell growth, and that its inhibition may represent a viable therapeutic strategy for lung cancer.

Proteomic study identifies Aurora-A-mediated regulation of alternative splicing through multiple splicing factors

The cell cycle regulator Aurora-A kinase presents an attractive target for cancer therapies, though its inhibition is also associated with toxic side effects. To gain a more nuanced understanding of Aurora-A function, we applied shotgun proteomics to identify 407 specific protein partners, including several splicing factors. Supporting a role in alternative splicing, we found that Aurora-A localizes to nuclear speckles, the storehouse of splicing proteins. Aurora-A interacts with and phosphorylates splicing factors both in vitro and in vivo, suggesting that it regulates alternative splicing by modulating the activity of these splicing factors. Consistently, Aurora-A inhibition significantly impacts the alternative splicing of 505 genes, with RNA motif analysis revealing an enrichment for Aurora-A interacting splicing factors. Additionally, we observed a significant positive correlation between the splicing events regulated by Aurora-A and those modulated by its interacting splicing factors. An interesting example is represented by CLK1 exon 4, which appears to be regulated by Aurora-A through SRSF3. Collectively, our findings highlight a broad role of Aurora-A in the regulation of alternative splicing.

Human introns contain conserved tissue-specific cryptic poison exons

Eukaryotic cells express a large number of transcripts from a single gene due to alternative splicing. Despite hundreds of thousands of splice isoforms being annotated in databases, it has been reported that the current exon catalogs remain incomplete. At the same time, introns of human protein-coding (PC) genes contain a large number of evolutionarily conserved elements with unknown function. Here, we explore the possibility that some of them represent cryptic exons that are expressed in rare conditions. We identified a group of cryptic exons that are similar to the annotated exons in terms of evolutionary conservation and RNA-seq read coverage in the Genotype-Tissue Expression dataset. Most of them were poison, i.e. generated an nonsense-mediated decay (NMD) isoform upon inclusion, and many showed signs of tissue-specific and cancer-specific expression and regulation. We performed RNA-seq in A549 cell line treated with cycloheximide to inactivate NMD and confirmed using quantitative polymerase chain reaction that seven of eight exons tested are, indeed, expressed. This study shows that introns of human PC genes contain cryptic poison exons, which reside in conserved intronic regions and remain not fully annotated due to insufficient representation in RNA-seq libraries.

A chemical screen identifies PRMT5 as a therapeutic vulnerability for paclitaxel-resistant triple-negative breast cancer

Paclitaxel-resistant triple negative breast cancer (TNBC) remains one of the most challenging breast cancers to treat. Here, using an epigenetic chemical probe screen, we uncover an acquired vulnerability of paclitaxel-resistant TNBC cells to protein arginine methyltransferases (PRMTs) inhibition. Analysis of cell lines and in-house clinical samples demonstrates that resistant cells evade paclitaxel killing through stabilizing mitotic chromatin assembly. Genetic or pharmacologic inhibition of PRMT5 alters RNA splicing, particularly intron retention of aurora kinases B (AURKB), leading to a decrease in protein expression, and finally results in selective mitosis catastrophe in paclitaxel-resistant cells. In addition, type I PRMT inhibition synergies with PRMT5 inhibition in suppressing tumor growth of drug-resistant cells through augmenting perturbation of AURKB-mediated mitotic signaling pathway. These findings are fully recapitulated in a patient-derived xenograft (PDX) model generated from a paclitaxel-resistant TNBC patient, providing the rationale for targeting PRMTs in paclitaxel-resistant TNBC.

Functional analysis of the GPAT4 gene mutation predicted to affect splicing

The GPAT4 gene is considered as a potential functional candidate for single nucleotide polymorphism (SNP) studies in dairy cattle breeding due to its association with dairy performance in cattle by encoding an enzyme responsible for the presence of diacylglycerols and triacylglycerols in milk. Using the example of the GPAT4 gene, we applied the minigene splicing assay to analyze the functional consequences of its variant that was predicted to affect normal splicing. The results of functional analysis revealed the sequence variations (rs442541537), transfection experiments in a wild type and mutant cell line model system demonstrated that the investigated mutation in the second intron of the GPAT4 gene was responsible for the presence of a second exon in mature messenger RNA (mRNA). The cases of its absence in the spliced mature mRNA transcript resulted in a truncated dysfunctional protein due to the appearance of a stop codon. Thus, the discovered SNP led to alternative splicing in pre-mRNA by the 'cassette exon' ('exon skipping') mechanism. The studied mutation can potentially be a molecular genetic marker for alternative splicing for the GPAT4 gene and, therefore contributes to economic benefits in cattle breeding programs.

Transcriptome-wide splicing network reveals specialized regulatory functions of the core spliceosome

The spliceosome is the complex molecular machinery that sequentially assembles on eukaryotic messenger RNA precursors to remove introns (pre-mRNA splicing), a physiologically regulated process altered in numerous pathologies. We report transcriptome-wide analyses upon systematic knock down of 305 spliceosome components and regulators in human cancer cells and the reconstruction of functional splicing factor networks that govern different classes of alternative splicing decisions. The results disentangle intricate circuits of splicing factor cross-regulation, reveal that the precise architecture of late-assembling U4/U6.U5 tri-small nuclear ribonucleoprotein (snRNP) complexes regulates splice site pairing, and discover an unprecedented division of labor among protein components of U1 snRNP for regulating exon definition and alternative 5' splice site selection. Thus, we provide a resource to explore physiological and pathological mechanisms of splicing regulation.

Exon nomenclature and classification of transcripts database (ENACTdb): a resource for analyzing alternative splicing mediated proteome diversity

Motivation

Gene transcripts are distinguished by the composition of their exons, and this different exon composition may contribute to advancing proteome complexity. Despite the availability of alternative splicing information documented in various databases, a ready association of exonic variations to the protein sequence remains a mammoth task.

Results

To associate exonic variation(s) with the protein systematically, we designed the Exon Nomenclature and Classification of Transcripts (ENACT) framework for uniquely annotating exons that tracks their loci in gene architecture context with encapsulating variations in splice site(s) and amino acid coding status. After ENACT annotation, predicted protein features (secondary structure/disorder/Pfam domains) are mapped to exon attributes. Thus, ENACTdb provides trackable exonic variation(s) association to isoform(s) and protein features, enabling the assessment of functional variation due to changes in exon composition. Such analyses can be readily performed through multiple views supported by the server. The exon-centric visualizations of ENACT annotated isoforms could provide insights on the functional repertoire of genes due to alternative splicing and its related processes and can serve as an important resource for the research community.

Availability and implementation

The database is publicly available at https://www.iscbglab.in/enactdb/. It contains protein-coding genes and isoforms for Caenorhabditis elegans, Drosophila melanogaster, Danio rerio, Mus musculus, and Homo sapiens.

ZMYND11 Functions in Bimodal Regulation of Latent Genes and Brain-like Splicing to Safeguard Corticogenesis

Despite the litany of pathogenic variants linked to neurodevelopmental disorders (NDD) including autism (ASD) and intellectual disability 1,2 , our understanding of the underlying mechanisms caused by risk genes remain unclear. Here, we leveraged a human pluripotent stem cell model to uncover the neurodevelopmental consequences of mutations in ZMYND11 , a newly implicated risk gene 3,4 . ZMYND11, known for its tumor suppressor function, encodes a histone-reader that recognizes sites of transcriptional elongation and acts as a co-repressor 5,6 . Our findings reveal that ZMYND11-deficient cortical neural stem cells showed upregulation of latent developmental pathways, impairing progenitor and neuron production. In addition to its role on histones, ZMYND11 controls a brain-specific isoform switch involving the splicing regulator RBFOX2. Extending our findings to other chromatin-related ASD risk factors revealed similar developmental pathway activation and splicing dysregulation, partially rescuable through ZMYND11's regulatory functions.

An ancient competition for the conserved branchpoint sequence influences physiological and evolutionary outcomes in splicing

Recognition of the intron branchpoint during spliceosome assembly is a multistep process that defines both mRNA structure and amount. A branchpoint sequence motif UACUAAC is variably conserved in eukaryotic genomes, but in some organisms more than one protein can recognize it. Here we show that SF1 and Quaking (QKI) compete for a subset of intron branchpoints with the sequence ACUAA. SF1 activates exon inclusion through this sequence, but QKI represses the inclusion of alternatively spliced exons with this intron branchpoint sequence. Using mutant reporters derived from a natural intron with two branchpoint-like sequences, we find that when either branchpoint sequence is mutated, the other is used as a branchpoint, but when both are present, neither is used due to high affinity binding and strong splicing repression by QKI. QKI occupancy at the dual branchpoint site directly prevents SF1 binding and subsequent recruitment of spliceosome-associated factors. Finally, the ectopic expression of QKI in budding yeast (which lacks QKI) is lethal, due at least in part to widespread splicing repression. In conclusion, QKI can function as a splicing repressor by directly competing with SF1/BBP for a subset of branchpoint sequences that closely mirror its high affinity binding site. This suggests that QKI and degenerate branchpoint sequences may have co-evolved as a means through which specific gene expression patterns could be maintained in QKI-expressing or non-expressing cells in metazoans, plants, and animals.

C2H2-zinc-finger transcription factors bind RNA and function in diverse post-transcriptional regulatory processes

Cys2-His2 zinc-finger proteins (C2H2-ZNFs) constitute the largest class of DNA-binding transcription factors (TFs) yet remain largely uncharacterized. Although certain family members, e.g., GTF3A, have been shown to bind both DNA and RNA, the extent to which C2H2-ZNFs interact with-and regulate-RNA-associated processes is not known. Using UV crosslinking and immunoprecipitation (CLIP), we observe that 148 of 150 analyzed C2H2-ZNFs bind directly to RNA in human cells. By integrating CLIP sequencing (CLIP-seq) RNA-binding maps for 50 of these C2H2-ZNFs with data from chromatin immunoprecipitation sequencing (ChIP-seq), protein-protein interaction assays, and transcriptome profiling experiments, we observe that the RNA-binding profiles of C2H2-ZNFs are generally distinct from their DNA-binding preferences and that they regulate a variety of post-transcriptional processes, including pre-mRNA splicing, cleavage and polyadenylation, and m6A modification of mRNA. Our results thus define a substantially expanded repertoire of C2H2-ZNFs that bind RNA and provide an important resource for elucidating post-transcriptional regulatory programs.

Differences in alternative splicing and their potential underlying factors between animals and plants

BACKGROUND

Alternative splicing (AS), a posttranscriptional process, contributes to the complexity of transcripts from a limited number of genes in a genome, and AS is considered a great source of genetic and phenotypic diversity in eukaryotes. In animals, AS is tightly regulated during the processes of cell growth and differentiation, and its dysregulation is involved in many diseases, including cancers. Likewise, in plants, AS occurs in all stages of plant growth and development, and it seems to play important roles in the rapid reprogramming of genes in response to environmental stressors. To date, the prevalence and functional roles of AS have been extensively reviewed in animals and plants. However, AS differences between animals and plants, especially their underlying molecular mechanisms and impact factors, are anecdotal and rarely reviewed.

AIM OF REVIEW

This review aims to broaden our understanding of AS roles in a variety of biological processes and provide insights into the underlying mechanisms and impact factors likely leading to AS differences between animals and plants.

KEY SCIENTIFIC CONCEPTS OF REVIEW

We briefly summarize the roles of AS regulation in physiological and biochemical activities in animals and plants. Then, we underline the differences in the process of AS between plants and animals and especially analyze the potential impact factors, such as gene exon/intron architecture, 5'/3' untranslated regions (UTRs), spliceosome components, chromatin dynamics and transcription speeds, splicing factors [serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs)], noncoding RNAs, and environmental stimuli, which might lead to the differences. Moreover, we compare the nonsense-mediated mRNA decay (NMD)-mediated turnover of the transcripts with a premature termination codon (PTC) in animals and plants. Finally, we summarize the current AS knowledge published in animals versus plants and discuss the potential development of disease therapies and superior crops in the future.

Protein isoform-centric therapeutics: expanding targets and increasing specificity

Most protein-coding genes produce multiple protein isoforms; however, these isoforms are commonly neglected in drug discovery. The expression of protein isoforms can be specific to a disease, tissue and/or developmental stage, and this specific expression can be harnessed to achieve greater drug specificity than pan-targeting of all gene products and to enable improved treatments for diseases caused by aberrant protein isoform production. In recent years, several protein isoform-centric therapeutics have been developed. Here, we collate these studies and clinical trials to highlight three distinct but overlapping modes of action for protein isoform-centric drugs: isoform switching, isoform introduction or depletion, and modulation of isoform activity. In addition, we discuss how protein isoforms can be used clinically as targets for cell type-specific drug delivery and immunotherapy, diagnostic biomarkers and sources of cancer neoantigens. Collectively, we emphasize the value of a focus on isoforms as a route to discovering drugs with greater specificity and fewer adverse effects. This approach could enable the targeting of proteins for which pan-inhibition of all isoforms is toxic and poorly tolerated.

Autoregulated splicing of TRA2β programs T cell fate in response to antigen-receptor stimulation

T cell receptor (TCR) sensitivity to peptide-major histocompatibility complex (MHC) dictates T cell fate. Canonical models of TCR sensitivity cannot be fully explained by transcriptional regulation. In this work, we identify a posttranscriptional regulatory mechanism of TCR sensitivity that guides alternative splicing of TCR signaling transcripts through an evolutionarily ultraconserved poison exon (PE) in the RNA-binding protein (RBP) TRA2β in mouse and human. TRA2β-PE splicing, seen during cancer and infection, was required for TCR-induced effector T cell expansion and function. Tra2β-PE skipping enhanced T cell response to antigen by increasing TCR sensitivity. As antigen levels decreased, Tra2β-PE reinclusion allowed T cell survival. Finally, we found that TRA2β-PE was first included in the genome of jawed vertebrates that were capable of TCR gene rearrangements. We propose that TRA2β-PE splicing acts as a gatekeeper of TCR sensitivity to shape T cell fate.

junctionCounts: comprehensive alternative splicing analysis and prediction of isoform-level impacts to the coding sequence

Alternative splicing (AS) is emerging as an important regulatory process for complex biological processes. Transcriptomic studies therefore commonly involve the identification and quantification of alternative processing events, but the need for predicting the functional consequences of changes to the relative inclusion of alternative events remains largely unaddressed. Many tools exist for the former task, albeit each constrained to its own event type definitions. Few tools exist for the latter task; each with significant limitations. To address these issues we developed junctionCounts, which captures both simple and complex pairwise AS events and quantifies them with straightforward exon-exon and exon-intron junction reads in RNA-seq data, performing competitively among similar tools in terms of sensitivity, false discovery rate and quantification accuracy. Its partner utility, cdsInsertion, identifies transcript coding sequence (CDS) information via in silico translation from annotated start codons, including the presence of premature termination codons. Finally, findSwitchEvents connects AS events with CDS information to predict the impact of individual events to the isoform-level CDS. We used junctionCounts to characterize splicing dynamics and NMD regulation during neuronal differentiation across four primates, demonstrating junctionCounts' capacity to robustly characterize AS in a variety of organisms and to predict its effect on mRNA isoform fate.

Accurate long-read transcript discovery and quantification at single-cell, pseudo-bulk and bulk resolution with Isosceles

Accurate detection and quantification of mRNA isoforms from nanopore long-read sequencing remains challenged by technical noise, particularly in single cells. To address this, we introduce Isosceles, a computational toolkit that outperforms other methods in isoform detection sensitivity and quantification accuracy across single-cell, pseudo-bulk and bulk resolution levels, as demonstrated using synthetic and biologically-derived datasets. Here we show Isosceles improves the fidelity of single-cell transcriptome quantification at the isoform-level, and enables flexible downstream analysis. As a case study, we apply Isosceles, uncovering coordinated splicing within and between neuronal differentiation lineages. Isosceles is suitable to be applied in diverse biological systems, facilitating studies of cellular heterogeneity across biomedical research applications.

Loss of cped1 does not affect bone and lean tissue in zebrafish

Human genetic studies have nominated Cadherin-like and PC-esterase Domain-containing 1 (CPED1) as a candidate target gene mediating bone mineral density (BMD) and fracture risk heritability. Recent efforts to define the role of CPED1 in bone in mouse and human models have revealed complex alternative splicing and inconsistent results arising from gene targeting, making its function in bone difficult to interpret. To better understand the role of CPED1 in adult bone mass and morphology, we conducted a comprehensive genetic and phenotypic analysis of cped1 in zebrafish, an emerging model for bone and mineral research. We analyzed two different cped1 mutant lines and performed deep phenotyping to characterize more than 200 measures of adult vertebral, craniofacial, and lean tissue morphology. We also examined alternative splicing of zebrafish cped1 and gene expression in various cell/tissue types. Our studies fail to support an essential role of cped1 in adult zebrafish bone. Specifically, homozygous mutants for both cped1 mutant alleles, which are expected to result in loss-of-function and impact all cped1 isoforms, exhibited no significant differences in the measures examined when compared to their respective wildtype controls, suggesting that cped1 does not significantly contribute to these traits. We identified sequence differences in critical residues of the catalytic triad between the zebrafish and mouse orthologs of CPED1, suggesting that differences in key residues, as well as distinct alternative splicing, could underlie different functions of CPED1 orthologs in the two species. Our studies fail to support a requirement of cped1 in zebrafish bone and lean tissue, adding to evidence that variants at 7q31.31 can act independently of CPED1 to influence BMD and fracture risk.

In silico RNA isoform screening to identify potential cancer driver exons with therapeutic applications

Alternative splicing is crucial for cancer progression and can be targeted pharmacologically, yet identifying driver exons genome-wide remains challenging. We propose identifying such exons by associating statistically gene-level cancer dependencies from knockdown viability screens with splicing profiles and gene expression. Our models predict the effects of splicing perturbations on cell proliferation from transcriptomic data, enabling in silico RNA screening and prioritizing targets for splicing-based therapies. We identified 1,073 exons impacting cell proliferation, many from genes not previously linked to cancer. Experimental validation confirms their influence on proliferation, especially in highly proliferative cancer cell lines. Integrating pharmacological screens with splicing dependencies highlights the potential driver exons affecting drug sensitivity. Our models also allow predicting treatment outcomes from tumor transcriptomes, suggesting applications in precision oncology. This study presents an approach to identifying cancer driver exon and their therapeutic potential, emphasizing alternative splicing as a cancer target.

The V-ATPase/ATG16L1 axis is controlled by the V1H subunit

Defects in organellar acidification indicate compromised or infected compartments. Recruitment of the autophagy-related ATG16L1 complex to pathologically neutralized organelles targets ubiquitin-like ATG8 molecules to perturbed membranes. How this process is coupled to proton gradient disruption is unclear. Here, we reveal that the V1H subunit of the vacuolar ATPase (V-ATPase) proton pump binds directly to ATG16L1. The V1H/ATG16L1 interaction only occurs within fully assembled V-ATPases, allowing ATG16L1 recruitment to be coupled to increased V-ATPase assembly following organelle neutralization. Cells lacking V1H fail to target ATG8s during influenza infection or after activation of the immune receptor stimulator of interferon genes (STING). We identify a loop within V1H that mediates ATG16L1 binding. A neuronal V1H isoform lacks this loop and is associated with attenuated ATG8 targeting in response to ionophores in primary murine and human iPSC-derived neurons. Thus, V1H controls ATG16L1 recruitment following proton gradient dissipation, suggesting that the V-ATPase acts as a cell-intrinsic damage sensor.

Integrative multi-omic analysis reveals conserved cell-projection deficits in human Down syndrome brains

Down syndrome (DS) is the most common genetic cause of cognitive disability. However, it is largely unclear how triplication of a small gene subset may impinge on diverse aspects of DS brain physiopathology. Here, we took a multi-omic approach and simultaneously analyzed by RNA-seq and proteomics the expression signatures of two diverse regions of human postmortem DS brains. We found that the overexpression of triplicated genes triggered global expression dysregulation, differentially affecting transcripts, miRNAs, and proteins involved in both known and novel biological candidate pathways. Among the latter, we observed an alteration in RNA splicing, specifically modulating the expression of genes involved in cytoskeleton and axonal dynamics in DS brains. Accordingly, we found an alteration in axonal polarization in neurons from DS human iPSCs and mice. Thus, our study provides an integrated multilayer expression database capable of identifying new potential targets to aid in designing future clinical interventions for DS.

Matrix stiffness-dependent regulation of immunomodulatory genes in human MSCs is associated with the lncRNA CYTOR

Cell-matrix interactions in 3D environments significantly differ from those in 2D cultures. As such, mechanisms of mechanotransduction in 2D cultures are not necessarily applicable to cell-encapsulating hydrogels that resemble features of tissue architecture. Accordingly, the characterization of molecular pathways in 3D matrices is expected to uncover insights into how cells respond to their mechanical environment in physiological contexts, and potentially also inform hydrogel-based strategies in cell therapies. In this study, a bone marrow-mimetic hydrogel was employed to systematically investigate the stiffness-responsive transcriptome of mesenchymal stromal cells. High matrix rigidity impeded integrin-collagen adhesion, resulting in changes in cell morphology characterized by a contractile network of actin proximal to the cell membrane. This resulted in a suppression of extracellular matrix-regulatory genes involved in the remodeling of collagen fibrils, as well as the upregulation of secreted immunomodulatory factors. Moreover, an investigation of long noncoding RNAs revealed that the cytoskeleton regulator RNA (CYTOR) contributes to these 3D stiffness-driven changes in gene expression. Knockdown of CYTOR using antisense oligonucleotides enhanced the expression of numerous mechanoresponsive cytokines and chemokines to levels exceeding those achievable by modulating matrix stiffness alone. Taken together, our findings further our understanding of mechanisms of mechanotransduction that are distinct from canonical mechanotransductive pathways observed in 2D cultures.

R2Dtool: integration and visualization of isoform-resolved RNA features

MOTIVATION

Long-read RNA sequencing enables the mapping of RNA modifications, structures, and protein-interaction sites at the resolution of individual transcript isoforms. To understand the functions of these RNA features, it is critical to analyze them in the context of transcriptomic and genomic annotations, such as open reading frames and splice junctions.

RESULTS

We have developed R2Dtool, a bioinformatics tool that integrates transcript-mapped information with transcript and genome annotations, allowing for the isoform-resolved analytics and graphical representation of RNA features in their genomic context. We illustrate R2Dtool's capability to integrate and expedite RNA feature analysis using epitranscriptomics data. R2Dtool facilitates the comprehensive analysis and interpretation of alternative transcript isoforms.

AVAILABILITY AND IMPLEMENTATION

R2Dtool is freely available under the MIT license at github.com/comprna/R2Dtool.

Long-read transcript sequencing identifies differential isoform expression in the entorhinal cortex in a transgenic model of tau pathology

Increasing evidence suggests that alternative splicing plays an important role in Alzheimer's disease (AD) pathology. We used long-read sequencing in combination with a novel bioinformatics tool (FICLE) to profile transcript diversity in the entorhinal cortex of female transgenic (TG) mice harboring a mutant form of human tau. Our analyses revealed hundreds of novel isoforms and identified differentially expressed transcripts - including specific isoforms of Apoe, App, Cd33, Clu, Fyn and Trem2 - associated with the development of tau pathology in TG mice. Subsequent profiling of the human cortex from AD individuals and controls revealed similar patterns of transcript diversity, including the upregulation of the dominant TREM2 isoform in AD paralleling the increased expression of the homologous transcript in TG mice. Our results highlight the importance of differential transcript usage, even in the absence of gene-level expression alterations, as a mechanism underpinning gene regulation in the development of AD neuropathology.

Biallelic variants in SNUPN cause a limb girdle muscular dystrophy with myofibrillar-like features

Alterations in RNA-splicing are a molecular hallmark of several neurological diseases, including muscular dystrophies, where mutations in genes involved in RNA metabolism or characterized by alterations in RNA splicing have been described. Here, we present five patients from two unrelated families with a limb-girdle muscular dystrophy (LGMD) phenotype carrying a biallelic variant in SNUPN gene. Snurportin-1, the protein encoded by SNUPN, plays an important role in the nuclear transport of small nuclear ribonucleoproteins (snRNPs), essential components of the spliceosome. We combine deep phenotyping, including clinical features, histopathology and muscle MRI, with functional studies in patient-derived cells and muscle biopsies to demonstrate that variants in SNUPN are the cause of a new type of LGMD according to current definition. Moreover, an in vivo model in Drosophila melanogaster further supports the relevance of Snurportin-1 in muscle. SNUPN patients show a similar phenotype characterized by proximal weakness starting in childhood, restrictive respiratory dysfunction and prominent contractures, although inter-individual variability in terms of severity even in individuals from the same family was found. Muscle biopsy showed myofibrillar-like features consisting of myotilin deposits and Z-disc disorganization. MRI showed predominant impairment of paravertebral, vasti, sartorius, gracilis, peroneal and medial gastrocnemius muscles. Conservation and structural analyses of Snurportin-1 p.Ile309Ser variant suggest an effect in nuclear-cytosol snRNP trafficking. In patient-derived fibroblasts and muscle, cytoplasmic accumulation of snRNP components is observed, while total expression of Snurportin-1 and snRNPs remains unchanged, which demonstrates a functional impact of SNUPN variant in snRNP metabolism. Furthermore, RNA-splicing analysis in patients' muscle showed widespread splicing deregulation, in particular in genes relevant for muscle development and splicing factors that participate in the early steps of spliceosome assembly. In conclusion, we report that SNUPN variants are a new cause of limb girdle muscular dystrophy with specific clinical, histopathological and imaging features, supporting SNUPN as a new gene to be included in genetic testing of myopathies. These results further support the relevance of splicing-related proteins in muscle disorders.

Development of a genome atlas for discriminating benign, preinvasive, and invasive lung nodules

To tackle misdiagnosis in lung cancer screening with low-dose computed tomography (LDCT), we aimed to compile a genome atlas for differentiating benign, preinvasive, and invasive lung nodules and characterize their molecular pathogenesis. We collected 432 lung nodule tissue samples from Chinese patients, spanning benign, atypical adenomatous hyperplasia (AAH), adenocarcinoma in situ (AIS), minimally invasive adenocarcinoma (MIA), and invasive adenocarcinoma (IA). We performed comprehensive sequencing, examining somatic variants, gene expressions, and methylation levels. Our findings uncovered EGFR and TP53 mutations as key drivers in - early lung cancer development, with EGFR mutation frequency increasing with disease progression. Both EGFR mutations and EGF/EGFR hypo-methylation activated the EGFR pathway, fueling cancer growth. Transcriptome analysis identified four lung nodule subtypes (G1-4) with distinct molecular features and immune cell infiltrations: EGFR-driven G1, EGFR/TP53 co-mutation G2, inflamed G3, stem-like G4. Estrogen/androgen response was associated with the EGFR pathway, proposing a new therapy combining tyrosine kinase inhibitors with antiestrogens. Preinvasive nodules exhibited stem cell pathway enrichment, potentially hindering invasion. Epigenetic regulation of various genes was essential for lung cancer initiation and development. This study provides insights into the molecular mechanism of neoplastic progression and identifies potential diagnostic biomarkers and therapeutic targets for lung cancer.

YTHDC1 m6A-dependent and m6A-independent functions converge to preserve the DNA damage response

Cells have evolved a robust and highly regulated DNA damage response to preserve their genomic integrity. Although increasing evidence highlights the relevance of RNA regulation, our understanding of its impact on a fully efficient DNA damage response remains limited. Here, through a targeted CRISPR-knockout screen, we identify RNA-binding proteins and modifiers that participate in the p53 response. Among the top hits, we find the m6A reader YTHDC1 as a master regulator of p53 expression. YTHDC1 binds to the transcription start sites of TP53 and other genes involved in the DNA damage response, promoting their transcriptional elongation. YTHDC1 deficiency also causes the retention of introns and therefore aberrant protein production of key DNA damage factors. While YTHDC1-mediated intron retention requires m6A, TP53 transcriptional pause-release is promoted by YTHDC1 independently of m6A. Depletion of YTHDC1 causes genomic instability and aberrant cancer cell proliferation mediated by genes regulated by YTHDC1. Our results uncover YTHDC1 as an orchestrator of the DNA damage response through distinct mechanisms of co-transcriptional mRNA regulation.

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.

High-throughput sensitive screening of small molecule modulators of microexon alternative splicing using dual Nano and Firefly luciferase reporters

Disruption of alternative splicing frequently causes or contributes to human diseases and disorders. Consequently, there is a need for efficient and sensitive reporter assays capable of screening chemical libraries for compounds with efficacy in modulating important splicing events. Here, we describe a screening workflow employing dual Nano and Firefly luciferase alternative splicing reporters that affords efficient, sensitive, and linear detection of small molecule responses. Applying this system to a screen of ~95,000 small molecules identified compounds that stimulate or repress the splicing of neuronal microexons, a class of alternative exons often disrupted in autism and activated in neuroendocrine cancers. One of these compounds rescues the splicing of several analyzed microexons in the cerebral cortex of an autism mouse model haploinsufficient for Srrm4, a major activator of brain microexons. We thus describe a broadly applicable high-throughput screening system for identifying candidate splicing therapeutics, and a resource of small molecule modulators of microexons with potential for further development in correcting aberrant splicing patterns linked to human disorders and disease.