Single-nuclei isoform RNA sequencing unlocks barcoded exon connectivity in frozen brain tissue

Exon-variant interplay and multi-modal evidence identify endocrine dysregulation in severe psychiatric disorders impacting excitatory neurons

Bipolar disorder (BD), major depressive disorder (MDD), and schizophrenia share genetic architecture, yet their molecular mechanisms remain elusive. Both common and rare genetic variants contribute to neural dysfunction, impacting cognition and behavior. This study investigates the molecular effects of genetic variants on human cortical single-cell types using a single-exon analysis approach. Integrating exon-level eQTLs (common variants influencing exon expression) and joint exon eQT-Scores (combining polygenic risk scores with exon-level gene expression) from a postmortem psychiatric cohort (BD = 15, MDD = 24, schizophrenia = 68, controls = 62) with schizophrenia-focused rare variant data from the SCHEMA consortium, we identified 110 core genes enriched in pathways including circadian entrainment (FDR = 0.02), cortisol synthesis and secretion (FDR = 0.026), and dopaminergic synapse (FDR = 0.038). Additional enriched pathways included hormone signaling (FDRs < 0.0298, including insulin, GnRH, aldosterone, and growth hormone pathways) and, notably, adrenergic signaling in cardiomyocytes (FDR = 0.0028). These pathways highlight shared molecular mechanisms in the three disorders. Single-nuclei RNA sequencing data from three cortical regions revealed that these core set genes are predominantly expressed in excitatory neuron layers 2-6 of the dorsolateral prefrontal cortex, linking molecular changes to cell types involved in cognitive dysfunction. Our results demonstrate the power of integrating multimodal genetic and transcriptomic data at the exon level. This approach moves beyond symptom-based diagnoses toward molecular classifications, identifying potential therapeutic targets for psychiatric disorders.

A molecular brain atlas reveals cellular shifts during the repair phase of stroke

Ischemic stroke triggers a cascade of pathological events that affect multiple cell types and often lead to incomplete functional recovery. Despite advances in single-cell technologies, the molecular and cellular responses that contribute to long-term post-stroke impairment remain poorly understood. To gain better insight into the underlying mechanisms, we generated a single-cell transcriptomic atlas from distinct brain regions using a mouse model of permanent focal ischemia at one month post-injury. Our findings reveal cell- and region-specific changes within the stroke-injured and peri-infarct brain tissue. For instance, GABAergic and glutamatergic neurons exhibited upregulated genes in signaling pathways involved in axon guidance and synaptic plasticity, and downregulated pathways associated with aerobic metabolism. Using cell-cell communication analysis, we identified increased strength in predicted interactions within stroke tissue among both neural and non-neural cells via signaling pathways such as those involving collagen, protein tyrosine phosphatase receptor, neuronal growth regulator, laminin, and several cell adhesion molecules. Furthermore, we found a strong correlation between mouse transcriptome responses after stroke and those observed in human nonfatal brain stroke lesions. Common molecular features were linked to inflammatory responses, extracellular matrix organization, and angiogenesis. Our findings provide a detailed resource for advancing our molecular understanding of stroke pathology and for discovering therapeutic targets in the repair phase of stroke recovery.

Genetic regulation of nascent RNA maturation revealed by direct RNA nanopore sequencing

Quantitative trait loci analyses have revealed an important role for genetic variants in regulating alternative splicing (AS) and alternative cleavage and polyadenylation (APA) in humans. Yet, these studies are generally performed with mature mRNA, so they report on the outcome rather than the processes of RNA maturation and thus may overlook how variants directly modulate pre-mRNA processing. The order in which the many introns of a human gene are removed can substantially influence AS, while nascent RNA polyadenylation can affect RNA stability and decay. However, how splicing order and poly(A) tail length are regulated by genetic variation has never been explored. Here, we used direct RNA nanopore sequencing to investigate allele-specific pre-mRNA maturation in 12 human lymphoblastoid cell lines. We find frequent splicing order differences between alleles and uncover significant single-nucleotide polymorphism (SNP)-splicing order associations in 17 genes. This includes SNPs located in or near splice sites as well as more distal intronic and exonic SNPs. Moreover, several genes showed allele-specific poly(A) tail lengths, many of which also have a skewed allelic abundance ratio. HLA class I transcripts, which encode proteins that play an essential role in antigen presentation, show the most allele-specific splicing orders, which frequently co-occur with allele-specific AS, APA, or poly(A) tail length differences. Together, our results expose new layers of genetic regulation of pre-mRNA maturation and highlight the power of long-read RNA sequencing for allele-specific analyses.

Alternative splicing in autism spectrum disorder: Recent insights from mechanisms to therapy

Alternative splicing (AS) is a vital and highly dynamic RNA regulatory mechanism that allows a single gene to generate multiple mRNA and protein isoforms. Dysregulation of AS has been identified as a key contributor to the pathogenesis of autism spectrum disorders (ASD). A comprehensive understanding of aberrant splicing mechanisms and their functional consequences in ASD can help uncover the molecular basis of the disorder and facilitate the development of therapeutic strategies. This review focuses on the major aberrant splicing events and key splicing regulators associated with ASD, highlighting their roles in linking defective splicing to ASD pathogenesis. In addition, a discussion of how emerging technologies, such as long-read sequencing, single-cell sequencing, spatial transcriptomics and CRISPR-Cas systems are offering novel insights into the role and mechanisms of AS in ASD is presented. Finally, the RNA splicing-based therapeutic strategies are evaluated, emphasizing their potential to address unmet clinical needs in ASD treatment.

Evolutionary divergence in CTCF-mediated chromatin topology drives transcriptional innovation in humans

Chromatin topology can impact gene regulation, but how evolutionary divergence in chromatin topology has shaped gene regulatory landscapes for distinctive human traits remains poorly understood. CTCF sites determine chromatin topology by forming domains and loops. Here, we show evolutionary divergence in CTCF-mediated chromatin topology at the domain and loop scales during primate evolution, elucidating distinct mechanisms for shaping regulatory landscapes. Human-specific divergent domains lead to a broad rewiring of transcriptional landscapes. Divergent CTCF loops concord with species-specific enhancer activity, influencing enhancer connectivity to target genes in a concordant yet constrained manner. Under this concordant mechanism, we establish the role of human-specific CTCF loops in shaping transcriptional isoform diversity, with functional implications for disease susceptibility. Furthermore, we validate the function of these human-specific CTCF loops using human forebrain organoids. This study advances our understanding of genetic evolution from the perspective of genome architecture.

Poison exons: tuning RNA splicing for targeted gene regulation

Poison exons (PEs) are a class of alternatively spliced exons whose inclusion targets mRNA transcripts for degradation via the nonsense-mediated decay (NMD) pathway. Although a role for NMD as an essential mRNA quality control pathway has long been appreciated, recent advances in RNA sequencing (RNA-seq) strategies and analyses have revealed that its coupling to RNA splicing is broadly used to regulate mRNA stability and abundance. Regulation of PE splicing affects patterns of targeted degradation across the transcriptome and influences gene expression in both healthy and disease states. Importantly, PEs represent a novel therapeutic opportunity to modulate the expression of disease-relevant genes with sequence-specific resolution. We review the emergence of PE splicing in endogenous gene regulation, its misregulation in disease, and the ways in which it can be leveraged for therapeutic benefit.

Are we there yet? Exploring astrocyte heterogeneity one cell at a time

Astrocytes are a highly abundant cell type in the brain and spinal cord. Like neurons, astrocytes can be molecularly and functionally distinct to fulfill specialized roles. Recent technical advances in sequencing-based single cell assays have driven an explosion of omics data characterizing astrocytes in the healthy, aged, injured, and diseased central nervous system. In this review, we will discuss recent studies which have furthered our understanding of astrocyte biology and heterogeneity, as well as discuss the limitations and challenges of sequencing-based single cell and spatial genomics methods and their potential future utility.

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.

A single-cell atlas to map sex-specific gene-expression changes in blood upon neurodegeneration

The clinical course and treatment of neurodegenerative disease are complicated by immune-system interference and chronic inflammatory processes, which remain incompletely understood. Mapping immune signatures in larger human cohorts through single-cell gene expression profiling supports our understanding of observed peripheral changes in neurodegeneration. Here, we employ single-cell gene expression profiling of over 909k peripheral blood mononuclear cells (PBMCs) from 121 healthy individuals, 48 patients with mild cognitive impairment (MCI), 46 with Parkinson's disease (PD), 27 with Alzheimer's disease (AD), and 15 with both PD and MCI. The dataset is interactively accessible through a freely available website ( https://www.ccb.uni-saarland.de/adrcsc ). In this work, we identify disease-associated changes in blood cell type composition and the gene expression in a sex-specific manner, offering insights into peripheral and solid tissue signatures in AD and PD.

Hybrid exons evolved by coupling transcription initiation and splicing at the nucleotide level

Exons within transcripts are traditionally classified as first, internal or last exons, each governed by different regulatory mechanisms. We recently described the widespread usage of 'hybrid' exons that serve as terminal or internal exons in different transcripts. Here, we employ an interpretable deep learning pipeline to dissect the sequence features governing the co-regulation of transcription initiation and splicing in hybrid exons. Using ENCODE data from human tissues, we identified 80 000 hybrid first-internal exons. These exons often possess a relaxed chromatin state, allowing transcription initiation within the gene body. Interestingly, transcription start sites of hybrid exons are typically centered at the 3' splice site, suggesting tight coupling between splicing and transcription initiation. We identified two subcategories of hybrid exons: the majority resemble internal exons, maintaining strong 3' splice sites, while a minority show enrichment in promoter elements, resembling first exons. Diving into the evolution of their sequences, we found that human hybrid exons with orthologous first exons in other species usually gained 3' splice sites or whole exons upstream, while those with orthologous internal exons often gained promoter elements. Overall, our findings unveil the intricate regulatory landscape of hybrid exons and reveal stronger connections between transcription initiation and RNA splicing than previously acknowledged.

Neuronal cell adhesion molecule (NRCAM) variant defined by microexon skipping is an essential, antigenically distinct, and targetable proteoform in high-grade glioma

To overcome the paucity of known tumor-specific surface antigens in pediatric high-grade glioma (pHGG), we contrasted splicing patterns in pHGGs and normal brain samples. Among alternative splicing events affecting extracellular protein domains, the most pervasive alteration was the skipping of ≤30 nucleotide-long microexons. Several of these skipped microexons mapped to L1-IgCAM family members, such as NRCAM . Bulk and single-nuclei short- and long-read RNA-seq revealed uniform skipping of NRCAM microexons 5 and 19 in virtually every pHGG sample. Importantly, the Δex5Δex19 (but not the full-length) NRCAM proteoform was essential for pHGG cell migration and invasion in vitro and tumor growth in vivo. We developed a monoclonal antibody selective for Δex5Δex19 NRCAM and demonstrated that "painting" of pHGG cells with this antibody enables killing by T cells armed with an FcRI-based universal immune receptor. Thus, pHGG-specific NRCAM and possibly other L1-IgCAM proteoforms are promising and highly selective targets for adoptive immunotherapies.

Statement of significance

Existing targets for chimeric antigen receptors (CAR)-armed T cells are often shared by CNS tumors and normal tissues, creating the potential for on-target/off-tumor toxicities. Here we demonstrate that in CNS tumors of glial origin, cell adhesion molecules have alternatively spliced proteoforms, which could be targeted by highly selective therapeutic antibodies.

Detecting gene expression in Caenorhabditis elegans

Reliable methods for detecting and analyzing gene expression are necessary tools for understanding development and investigating biological responses to genetic and environmental perturbation. With its fully sequenced genome, invariant cell lineage, transparent body, wiring diagram, detailed anatomy, and wide array of genetic tools, Caenorhabditis elegans is an exceptionally useful model organism for linking gene expression to cellular phenotypes. The development of new techniques in recent years has greatly expanded our ability to detect gene expression at high resolution. Here, we provide an overview of gene expression methods for C. elegans, including techniques for detecting transcripts and proteins in situ, bulk RNA sequencing of whole worms and specific tissues and cells, single-cell RNA sequencing, and high-throughput proteomics. We discuss important considerations for choosing among these techniques and provide an overview of publicly available online resources for gene expression data.

RNA Isoform Diversity in Human Neurodegenerative Diseases

Single-nucleus RNA-sequencing (snRNA-seq) has revealed new levels of cellular organization and diversity within the human brain. However, full-length mRNA isoforms are not resolved in typical snRNA-seq analyses using short-read sequencing that cannot capture full-length transcripts. Here we combine standard 10x Genomics short-read snRNA-seq with targeted PacBio long-read snRNA-seq to examine isoforms of genes associated with neurological diseases at the single-cell level from prefrontal cortex samples of diseased and nondiseased human brain, assessing over 165,000 cells. Samples from 25 postmortem donors with Alzheimer's disease (AD), dementia with Lewy bodies (DLB), or Parkinson's disease (PD), along with age-matched controls, were compared. Analysis of the short-read libraries identified shared and distinct gene expression changes across the diseases. The same libraries were then assayed using enrichment probes to target 50 disease-related genes followed by long-read PacBio sequencing, enabling linkage between cell type and isoform expression. Vast mRNA isoform diversity was observed in all 50 targeted genes, even those that were not differentially expressed in the short-read data. We also developed an informatics method for detection of isoform structural differences in novel isoforms versus the reference annotation. These data expand available single-cell datasets of the human prefrontal cortical transcriptome with combined short- and long-read sequencing across AD, DLB, and PD, revealing increased mRNA isoform diversity that may contribute to disease features and could potentially represent therapeutic targets for neurodegenerative diseases.

(Alternative) transcription start sites as regulators of RNA processing

Alternative transcription start site usage (ATSS) is a widespread regulatory strategy that enables genes to choose between multiple genomic loci for initiating transcription. This mechanism is tightly controlled during development and is often altered in disease states. In this review, we examine the growing evidence highlighting a role for transcription start sites (TSSs) in the regulation of mRNA isoform selection during and after transcription. We discuss how the choice of transcription initiation sites influences RNA processing and the importance of this crosstalk for cell identity and organism function. We also speculate on possible mechanisms underlying the integration of transcriptional and post-transcriptional processes.

Steering research on mRNA splicing in cancer towards clinical translation

Splicing factors are affected by recurrent somatic mutations and copy number variations in several types of haematologic and solid malignancies, which is often seen as prima facie evidence that splicing aberrations can drive cancer initiation and progression. However, numerous spliceosome components also 'moonlight' in DNA repair and other cellular processes, making their precise role in cancer difficult to pinpoint. Still, few would deny that dysregulated mRNA splicing is a pervasive feature of most cancers. Correctly interpreting these molecular fingerprints can reveal novel tumour vulnerabilities and untapped therapeutic opportunities. Yet multiple technological challenges, lingering misconceptions, and outstanding questions hinder clinical translation. To start with, the general landscape of splicing aberrations in cancer is not well defined, due to limitations of short-read RNA sequencing not adept at resolving complete mRNA isoforms, as well as the shallow read depth inherent in long-read RNA-sequencing, especially at single-cell level. Although individual cancer-associated isoforms are known to contribute to cancer progression, widespread splicing alterations could be an equally important and, perhaps, more readily actionable feature of human cancers. This is to say that in addition to 'repairing' mis-spliced transcripts, possible therapeutic avenues include exacerbating splicing aberration with small-molecule spliceosome inhibitors, targeting recurrent splicing aberrations with synthetic lethal approaches, and training the immune system to recognize splicing-derived neoantigens.

Understanding isoform expression by pairing long-read sequencing with single-cell and spatial transcriptomics

RNA isoform diversity, produced via alternative splicing, and alternative usage of transcription start and poly(A) sites, results in varied transcripts being derived from the same gene. Distinct isoforms can play important biological roles, including by changing the sequences or expression levels of protein products. The first single-cell approaches to RNA sequencing-and later, spatial approaches-which are now widely used for the identification of differentially expressed genes, rely on short reads and offer the ability to transcriptomically compare different cell types but are limited in their ability to measure differential isoform expression. More recently, long-read sequencing methods have been combined with single-cell and spatial technologies in order to characterize isoform expression. In this review, we provide an overview of the emergence of single-cell and spatial long-read sequencing and discuss the challenges associated with the implementation of these technologies and interpretation of these data. We discuss the opportunities they offer for understanding the relationships between the distinct variable elements of transcript molecules and highlight some of the ways in which they have been used to characterize isoforms' roles in development and pathology. Single-nucleus long-read sequencing, a special case of the single-cell approach, is also discussed. We attempt to cover both the limitations of these technologies and their significant potential for expanding our still-limited understanding of the biological roles of RNA isoforms.

Diffusion imaging genomics provides novel insight into early mechanisms of cerebral small vessel disease

Cerebral small vessel disease (cSVD) is a leading cause of stroke and dementia. Genetic risk loci for white matter hyperintensities (WMH), the most common MRI-marker of cSVD in older age, were recently shown to be significantly associated with white matter (WM) microstructure on diffusion tensor imaging (signal-based) in young adults. To provide new insights into these early changes in WM microstructure and their relation with cSVD, we sought to explore the genetic underpinnings of cutting-edge tissue-based diffusion imaging markers across the adult lifespan. We conducted a genome-wide association study of neurite orientation dispersion and density imaging (NODDI) markers in young adults (i-Share study: N = 1 758, (mean[range]) 22.1[18-35] years), with follow-up in young middle-aged (Rhineland Study: N = 714, 35.2[30-40] years) and late middle-aged to older individuals (UK Biobank: N = 33 224, 64.3[45-82] years). We identified 21 loci associated with NODDI markers across brain regions in young adults. The most robust association, replicated in both follow-up cohorts, was with Neurite Density Index (NDI) at chr5q14.3, a known WMH locus in VCAN. Two additional loci were replicated in UK Biobank, at chr17q21.2 with NDI, and chr19q13.12 with Orientation Dispersion Index (ODI). Transcriptome-wide association studies showed associations of STAT3 expression in arterial and adipose tissue (chr17q21.2) with NDI, and of several genes at chr19q13.12 with ODI. Genetic susceptibility to larger WMH volume, but not to vascular risk factors, was significantly associated with decreased NDI in young adults, especially in regions known to harbor WMH in older age. Individually, seven of 25 known WMH risk loci were associated with NDI in young adults. In conclusion, we identified multiple novel genetic risk loci associated with NODDI markers, particularly NDI, in early adulthood. These point to possible early-life mechanisms underlying cSVD and to processes involving remyelination, neurodevelopment and neurodegeneration, with a potential for novel approaches to prevention.

The multicellular compartmentation of plant specialized metabolism

The phenomenon of multicellular compartmentation in biosynthetic pathways has been documented for only a limited subset of specialized metabolites, despite its hypothesized significance in facilitating plant survival and adaptation to environmental stress. Transporters that shuttle metabolic intermediates between cells are hypothesized to be integral components enabling compartmentalized biosynthesis. Nevertheless, our understanding of the multicellular compartmentation of plant specialized metabolism and the associated intermediate transporters remains incomplete. The emergence of single-cell and spatial multiomics techniques holds promise for shedding light on unresolved questions in this field, such as the prevalence of multicellular compartmentation across the plant kingdom and the specific types of specialized metabolites whose biosynthetic pathways are prone to compartmentation. Advancing our understanding of the mechanisms underlying multicellular compartmentation will contribute to improving the production of specialized target metabolites through metabolic engineering or synthetic biology.

An integrative single-cell atlas for exploring the cellular and temporal specificity of genes related to neurological disorders during human brain development

Single-cell technologies have enhanced comprehensive knowledge regarding the human brain by facilitating an extensive transcriptomic census across diverse brain regions. Nevertheless, understanding the cellular and temporal specificity of neurological disorders remains ambiguous due to developmental variations. To address this gap, we illustrated the dynamics of disorder risk gene expression under development by integrating multiple single-cell RNA sequencing datasets. We constructed a comprehensive single-cell atlas of the developing human brain, encompassing 393,060 single cells across diverse developmental stages. Temporal analysis revealed the distinct expression patterns of disorder risk genes, including those associated with autism, highlighting their temporal regulation in different neuronal and glial lineages. We identified distinct neuronal lineages that diverged across developmental stages, each exhibiting temporal-specific expression patterns of disorder-related genes. Lineages of nonneuronal cells determined by molecular profiles also showed temporal-specific expression, indicating a link between cellular maturation and the risk of disorder. Furthermore, we explored the regulatory mechanisms involved in early brain development, revealing enriched patterns of fetal cell types associated with neuronal disorders indicative of the prenatal stage's influence on disease determination. Our findings facilitate unbiased comparisons of cell type‒disorder associations and provide insight into dynamic alterations in risk genes during development, paving the way for a deeper understanding of neurological disorders.

Fine mapping of RNA isoform diversity using an innovative targeted long-read RNA sequencing protocol with novel dedicated bioinformatics pipeline

BACKGROUND

Solving the structure of mRNA transcripts is a major challenge for both research and molecular diagnostic purposes. Current approaches based on short-read RNA sequencing and RT-PCR techniques cannot fully explore the complexity of transcript structure. The emergence of third-generation long-read sequencing addresses this problem by solving this sequence directly. However, genes with low expression levels are difficult to study with the whole transcriptome sequencing approach. To fix this technical limitation, we propose a novel method to capture transcripts of a gene panel using a targeted enrichment approach suitable for Pacific Biosciences and Oxford Nanopore Technologies platforms.

RESULTS

We designed a set of probes to capture transcripts of a panel of genes involved in hereditary breast and ovarian cancer syndrome. We present SOSTAR (iSofOrmS annoTAtoR), a versatile pipeline to assemble, quantify and annotate isoforms from long read sequencing using a new tool specially designed for this application. The significant enrichment of transcripts by our capture protocol, together with the SOSTAR annotation, allowed the identification of 1,231 unique transcripts within the gene panel from the eight patients sequenced. The structure of these transcripts was annotated with a resolution of one base relative to a reference transcript. All major alternative splicing events of the BRCA1 and BRCA2 genes described in the literature were found. Complex splicing events such as pseudoexons were correctly annotated. SOSTAR enabled the identification of abnormal transcripts in the positive controls. In addition, a case of unexplained inheritance in a family with a history of breast and ovarian cancer was solved by identifying an SVA retrotransposon in intron 13 of the BRCA1 gene.

CONCLUSIONS

We have validated a new protocol for the enrichment of transcripts of interest using probes adapted to the ONT and PacBio platforms. This protocol allows a complete description of the alternative structures of transcripts, the estimation of their expression and the identification of aberrant transcripts in a single experiment. This proof-of-concept opens new possibilities for RNA structure exploration in both research and molecular diagnostics.

Genetic regulation of nascent RNA maturation revealed by direct RNA nanopore sequencing

Quantitative trait loci analyses have revealed an important role for genetic variants in regulating alternative splicing (AS) and alternative cleavage and polyadenylation (APA) in humans. Yet, these studies are generally performed with mature mRNA, so they report on the outcome rather than the processes of RNA maturation and thus may overlook how variants directly modulate pre-mRNA processing. The order in which the many introns of a human gene are removed can substantially influence AS, while nascent RNA polyadenylation can affect RNA stability and decay. However, how splicing order and poly(A) tail length are regulated by genetic variation has never been explored. Here, we used direct RNA nanopore sequencing to investigate allele-specific pre-mRNA maturation in 12 human lymphoblastoid cell lines. We found frequent splicing order differences between alleles and uncovered significant single nucleotide polymorphism (SNP)-splicing order associations in 17 genes. This included SNPs located in or near splice sites as well as more distal intronic and exonic SNPs. Moreover, several genes showed allele-specific poly(A) tail lengths, many of which also had a skewed allelic abundance ratio. HLA class I transcripts, which encode proteins that play an essential role in antigen presentation, showed the most allele-specific splicing orders, which frequently co-occurred with allele-specific AS, APA or poly(A) tail length differences. Together, our results expose new layers of genetic regulation of pre-mRNA maturation and highlight the power of long-read RNA sequencing for allele-specific analyses.

A molecular brain atlas reveals cellular shifts during the repair phase of stroke

Ischemic stroke triggers a cascade of pathological events that affect multiple cell types and often lead to incomplete functional recovery. Despite advances in single-cell technologies, the molecular and cellular responses that contribute to long-term post-stroke impairment remain poorly understood. To gain better insight into the underlying mechanisms, we generated a single-cell transcriptomic atlas from distinct brain regions using a mouse model of permanent focal ischemia at one month post-injury. Our findings reveal cell- and region-specific changes within the stroke-injured and peri-infarct brain tissue. For instance, GABAergic and glutamatergic neurons exhibited upregulated genes in signaling pathways involved in axon guidance and synaptic plasticity, and downregulated pathways associated with aerobic metabolism. Using cell-cell communication analysis, we identified increased strength in predicted interactions within stroke tissue among both neural and non-neural cells via signaling pathways such as those involving collagen, protein tyrosine phosphatase receptor, neuronal growth regulator, laminin, and several cell adhesion molecules. Furthermore, we found a strong correlation between mouse transcriptome responses after stroke and those observed in human nonfatal brain stroke lesions. Common molecular features were linked to inflammatory responses, extracellular matrix organization, and angiogenesis. Our findings provide a detailed resource for advancing our molecular understanding of stroke pathology and for discovering therapeutic targets in the repair phase of stroke recovery.

Multi-omic profiling of pathogen-stimulated primary immune cells

We performed long-read transcriptome and proteome profiling of pathogen-stimulated peripheral blood mononuclear cells (PBMCs) from healthy donors to discover new transcript and protein isoforms expressed during immune responses to diverse pathogens. Long-read transcriptome profiling reveals novel sequences and isoform switching induced upon pathogen stimulation, including transcripts that are difficult to detect using traditional short-read sequencing. Widespread loss of intron retention occurs as a common result of all pathogen stimulations. We highlight novel transcripts of NFKB1 and CASP1 that may indicate novel immunological mechanisms. RNA expression differences did not result in differences in the amounts of secreted proteins. Clustering analysis of secreted proteins revealed a correlation between chemokine (receptor) expression on the RNA and protein levels in C. albicans- and poly(I:C)-stimulated PBMCs. Isoform aware long-read sequencing of pathogen-stimulated immune cells highlights the potential of these methods to identify novel transcripts, revealing a more complex transcriptome landscape than previously appreciated.

Targeted DNA-seq and RNA-seq of Reference Samples with Short-read and Long-read Sequencing

Next-generation sequencing (NGS) has revolutionized genomic research by enabling high-throughput, cost-effective genome and transcriptome sequencing accelerating personalized medicine for complex diseases, including cancer. Whole genome/transcriptome sequencing (WGS/WTS) provides comprehensive insights, while targeted sequencing is more cost-effective and sensitive. In comparison to short-read sequencing, which still dominates the field due to high speed and cost-effectiveness, long-read sequencing can overcome alignment limitations and better discriminate similar sequences from alternative transcripts or repetitive regions. Hybrid sequencing combines the best strengths of different technologies for a more comprehensive view of genomic/transcriptomic variations. Understanding each technology's strengths and limitations is critical for translating cutting-edge technologies into clinical applications. In this study, we sequenced DNA and RNA libraries of reference samples using various targeted DNA and RNA panels and the whole transcriptome on both short-read and long-read platforms. This study design enables a comprehensive analysis of sequencing technologies, targeting protocols, and library preparation methods. Our expanded profiling landscape establishes a reference point for assessing current sequencing technologies, facilitating informed decision-making in genomic research and precision medicine.

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.

Neuronal RNA processing: cross-talk between transcriptional regulation and RNA-binding proteins

In the nervous system, alternative RNA processing is particularly prevalent, which results in the expression of thousands of transcript variants found in no other tissue. Neuron-specific RNA-binding proteins co-transcriptionally regulate alternative splicing, alternative polyadenylation, and RNA editing, thereby shaping the RNA identity of nervous system cells. Recent evidence suggests that interactions between RNA-binding proteins and cis-regulatory elements such as promoters and enhancers play a role in the determination of neuron-specific expression profiles. Here, we discuss possible mechanisms through which transcription and RNA processing cross-talk to generate the uniquely complex neuronal transcriptome, with a focus on alternative 3'-end formation.

Neuron type-specific proteomics reveals distinct Shank3 proteoforms in iSPNs and dSPNs lead to striatal synaptopathy in Shank3B-/- mice

Combinatorial expression of postsynaptic proteins underlies synapse diversity within and between neuron types. Thus, characterization of neuron-type-specific postsynaptic proteomes is key to obtaining a deeper understanding of discrete synaptic properties and how selective dysfunction manifests in synaptopathies. To overcome the limitations associated with bulk measures of synaptic protein abundance, we developed a biotin proximity protein tagging probe to characterize neuron-type-specific postsynaptic proteomes in vivo. We found Shank3 protein isoforms are differentially expressed by direct and indirect pathway spiny projection neurons (dSPNs and iSPNs). Investigation of Shank3B-/- mice lacking exons 13-16 within the Shank3 gene, reveal distinct Shank3 protein isoform expression in iSPNs and dSPNs. In Shank3B-/- striatum, Shank3E and Shank3NT are expressed by dSPNs but are undetectable in iSPNs. Proteomic analysis indicates significant and selective alterations in the postsynaptic proteome of Shank3B-/- iSPNs. Correspondingly, the deletion of exons 13-16 diminishes dendritic spine density, reduces spine head diameter, and hampers corticostriatal synaptic transmission in iSPNs. Remarkably, reintroducing Shank3E in adult Shank3B-/- iSPNs significantly rectifies the observed dendritic spine morphological and corticostriatal synaptic transmission deficits. We report unexpected cell-type specific synaptic protein isoform expression which could play a key causal role in specifying synapse diversity and selective synapse dysfunction in synaptopathies.

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

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

Advances in single-cell long-read sequencing technologies

With an increase in accuracy and throughput of long-read sequencing technologies, they are rapidly being assimilated into the single-cell sequencing pipelines. For transcriptome sequencing, these techniques provide RNA isoform-level information in addition to the gene expression profiles. Long-read sequencing technologies not only help in uncovering complex patterns of cell-type specific splicing, but also offer unprecedented insights into the origin of cellular complexity and thus potentially new avenues for drug development. Additionally, single-cell long-read DNA sequencing enables high-quality assemblies, structural variant detection, haplotype phasing, resolving high-complexity regions, and characterization of epigenetic modifications. Given that significant progress has primarily occurred in single-cell RNA isoform sequencing (scRiso-seq), this review will delve into these advancements in depth and highlight the practical considerations and operational challenges, particularly pertaining to downstream analysis. We also aim to offer a concise introduction to complementary technologies for single-cell sequencing of the genome, epigenome and epitranscriptome. We conclude by identifying certain key areas of innovation that may drive these technologies further and foster more widespread application in biomedical science.

Single-cell long-read sequencing-based mapping reveals specialized splicing patterns in developing and adult mouse and human brain

RNA isoforms influence cell identity and function. However, a comprehensive brain isoform map was lacking. We analyze single-cell RNA isoforms across brain regions, cell subtypes, developmental time points and species. For 72% of genes, full-length isoform expression varies along one or more axes. Splicing, transcription start and polyadenylation sites vary strongly between cell types, influence protein architecture and associate with disease-linked variation. Additionally, neurotransmitter transport and synapse turnover genes harbor cell-type variability across anatomical regions. Regulation of cell-type-specific splicing is pronounced in the postnatal day 21-to-postnatal day 28 adolescent transition. Developmental isoform regulation is stronger than regional regulation for the same cell type. Cell-type-specific isoform regulation in mice is mostly maintained in the human hippocampus, allowing extrapolation to the human brain. Conversely, the human brain harbors additional cell-type specificity, suggesting gain-of-function isoforms. Together, this detailed single-cell atlas of full-length isoform regulation across development, anatomical regions and species reveals an unappreciated degree of isoform variability across multiple axes.

Developmental isoform diversity in the human neocortex informs neuropsychiatric risk mechanisms

RNA splicing is highly prevalent in the brain and has strong links to neuropsychiatric disorders; yet, the role of cell type-specific splicing and transcript-isoform diversity during human brain development has not been systematically investigated. In this work, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone and cortical plate regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 distinct isoforms, of which 72.6% were novel (not previously annotated in Gencode version 33), and uncovered a substantial contribution of transcript-isoform diversity-regulated by RNA binding proteins-in defining cellular identity in the developing neocortex. We leveraged this comprehensive isoform-centric gene annotation to reprioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders.

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

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

Single-cell and spatial transcriptomics: Bridging current technologies with long-read sequencing

Single-cell technologies have transformed biomedical research over the last decade, opening up new possibilities for understanding cellular heterogeneity, both at the genomic and transcriptomic level. In addition, more recent developments of spatial transcriptomics technologies have made it possible to profile cells in their tissue context. In parallel, there have been substantial advances in sequencing technologies, and the third generation of methods are able to produce reads that are tens of kilobases long, with error rates matching the second generation short reads. Long reads technologies make it possible to better map large genome rearrangements and quantify isoform specific abundances. This further improves our ability to characterize functionally relevant heterogeneity. Here, we show how researchers have begun to combine single-cell, spatial transcriptomics, and long-read technologies, and how this is resulting in powerful new approaches to profiling both the genome and the transcriptome. We discuss the achievements so far, and we highlight remaining challenges and opportunities.

Relapse to cocaine seeking is regulated by medial habenula NR4A2/NURR1 in mice

Drugs of abuse can persistently change the reward circuit in ways that contribute to relapse behavior, partly via mechanisms that regulate chromatin structure and function. Nuclear orphan receptor subfamily4 groupA member2 (NR4A2, also known as NURR1) is an important effector of histone deacetylase 3 (HDAC3)-dependent mechanisms in persistent memory processes and is highly expressed in the medial habenula (MHb), a region that regulates nicotine-associated behaviors. Here, expressing the Nr4a2 dominant negative (Nurr2c) in the MHb blocks reinstatement of cocaine seeking in mice. We use single-nucleus transcriptomics to characterize the molecular cascade following Nr4a2 manipulation, revealing changes in transcriptional networks related to addiction, neuroplasticity, and GABAergic and glutamatergic signaling. The network controlled by NR4A2 is characterized using a transcription factor regulatory network inference algorithm. These results identify the MHb as a pivotal regulator of relapse behavior and demonstrate the importance of NR4A2 as a key mechanism driving the MHb component of relapse.

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.

Advances in single-cell omics and multiomics for high-resolution molecular profiling

Single-cell omics technologies have revolutionized molecular profiling by providing high-resolution insights into cellular heterogeneity and complexity. Traditional bulk omics approaches average signals from heterogeneous cell populations, thereby obscuring important cellular nuances. Single-cell omics studies enable the analysis of individual cells and reveal diverse cell types, dynamic cellular states, and rare cell populations. These techniques offer unprecedented resolution and sensitivity, enabling researchers to unravel the molecular landscape of individual cells. Furthermore, the integration of multimodal omics data within a single cell provides a comprehensive and holistic view of cellular processes. By combining multiple omics dimensions, multimodal omics approaches can facilitate the elucidation of complex cellular interactions, regulatory networks, and molecular mechanisms. This integrative approach enhances our understanding of cellular systems, from development to disease. This review provides an overview of the recent advances in single-cell and multimodal omics for high-resolution molecular profiling. We discuss the principles and methodologies for representatives of each omics method, highlighting the strengths and limitations of the different techniques. In addition, we present case studies demonstrating the applications of single-cell and multimodal omics in various fields, including developmental biology, neurobiology, cancer research, immunology, and precision medicine.

ScISOr-ATAC reveals convergent and divergent splicing and chromatin specificities between matched cell types across cortical regions, evolution, and in Alzheimer's Disease

Multimodal measurements have become widespread in genomics, however measuring open chromatin accessibility and splicing simultaneously in frozen brain tissues remains unconquered. Hence, we devised Single-Cell-ISOform-RNA sequencing coupled with the Assay-for-Transposase-Accessible-Chromatin (ScISOr-ATAC). We utilized ScISOr-ATAC to assess whether chromatin and splicing alterations in the brain convergently affect the same cell types or divergently different ones. We applied ScISOr-ATAC to three major conditions: comparing (i) the Rhesus macaque (Macaca mulatta) prefrontal cortex (PFC) and visual cortex (VIS), (ii) cross species divergence of Rhesus macaque versus human PFC, as well as (iii) dysregulation in Alzheimer's disease in human PFC. We found that among cortical-layer biased excitatory neuron subtypes, splicing is highly brain-region specific for L3-5/L6 IT_RORB neurons, moderately specific in L2-3 IT_CUX2.RORB neurons and unspecific in L2-3 IT_CUX2 neurons. In contrast, at the chromatin level, L2-3 IT_CUX2.RORB neurons show the highest brain-region specificity compared to other subtypes. Likewise, when comparing human and macaque PFC, strong evolutionary divergence on one molecular modality does not necessarily imply strong such divergence on another molecular level in the same cell type. Finally, in Alzheimer's disease, oligodendrocytes show convergently high dysregulation in both chromatin and splicing. However, chromatin and splicing dysregulation most strongly affect distinct oligodendrocyte subtypes. Overall, these results indicate that chromatin and splicing can show convergent or divergent results depending on the performed comparison, justifying the need for their concurrent measurement to investigate complex systems. Taken together, ScISOr-ATAC allows for the characterization of single-cell splicing and chromatin patterns and the comparison of sample groups in frozen brain samples.

Multi-ancestry genetic analysis of gene regulation in coronary arteries prioritizes disease risk loci

Genome-wide association studies (GWASs) have identified hundreds of risk loci for coronary artery disease (CAD). However, non-European populations are underrepresented in GWASs, and the causal gene-regulatory mechanisms of these risk loci during atherosclerosis remain unclear. We incorporated local ancestry and haplotypes to identify quantitative trait loci for expression (eQTLs) and splicing (sQTLs) in coronary arteries from 138 ancestrally diverse Americans. Of 2,132 eQTL-associated genes (eGenes), 47% were previously unreported in coronary artery; 19% exhibited cell-type-specific expression. Colocalization revealed subgroups of eGenes unique to CAD and blood pressure GWAS. Fine-mapping highlighted additional eGenes, including TBX20 and IL5. We also identified sQTLs for 1,690 genes, among which TOR1AIP1 and ULK3 sQTLs demonstrated the importance of evaluating splicing to accurately identify disease-relevant isoform expression. Our work provides a patient-derived coronary artery eQTL resource and exemplifies the need for diverse study populations and multifaceted approaches to characterize gene regulation in disease processes.

mRNA initiation and termination are spatially coordinated

The expression of a precise mRNA transcriptome is crucial for establishing cell identity and function, with dozens of alternative isoforms produced for a single gene sequence. The regulation of mRNA isoform usage occurs by the coordination of co-transcriptional mRNA processing mechanisms across a gene. Decisions involved in mRNA initiation and termination underlie the largest extent of mRNA isoform diversity, but little is known about any relationships between decisions at both ends of mRNA molecules. Here, we systematically profile the joint usage of mRNA transcription start sites (TSSs) and polyadenylation sites (PASs) across tissues and species. Using both short and long read RNA-seq data, we observe that mRNAs preferentially using upstream TSSs also tend to use upstream PASs, and congruently, the usage of downstream sites is similarly paired. This observation suggests that mRNA 5' end choice may directly influence mRNA 3' ends. Our results suggest a novel "Positional Initiation-Termination Axis" (PITA), in which the usage of alternative terminal sites are coupled based on the order in which they appear in the genome. PITA isoforms are more likely to encode alternative protein domains and use conserved sites. PITA is strongly associated with the length of genomic features, such that PITA is enriched in longer genes with more area devoted to regions that regulate alternative 5' or 3' ends. Strikingly, we found that PITA genes are more likely than non-PITA genes to have multiple, overlapping chromatin structural domains related to pairing of ordinally coupled start and end sites. In turn, PITA coupling is also associated with fast RNA Polymerase II (RNAPII) trafficking across these long gene regions. Our findings indicate that a combination of spatial and kinetic mechanisms couple transcription initiation and mRNA 3' end decisions based on ordinal position to define the expression mRNA isoforms.

ZEBRA: a hierarchically integrated gene expression atlas of the murine and human brain at single-cell resolution

The molecular causes and mechanisms of neurodegenerative diseases remain poorly understood. A growing number of single-cell studies have implicated various neural, glial, and immune cell subtypes to affect the mammalian central nervous system in many age-related disorders. Integrating this body of transcriptomic evidence into a comprehensive and reproducible framework poses several computational challenges. Here, we introduce ZEBRA, a large single-cell and single-nucleus RNA-seq database. ZEBRA integrates and normalizes gene expression and metadata from 33 studies, encompassing 4.2 million human and mouse brain cells sampled from 39 brain regions. It incorporates samples from patients with neurodegenerative diseases like Alzheimer's disease, Parkinson's disease, and Multiple sclerosis, as well as samples from relevant mouse models. We employed scVI, a deep probabilistic auto-encoder model, to integrate the samples and curated both cell and sample metadata for downstream analysis. ZEBRA allows for cell-type and disease-specific markers to be explored and compared between sample conditions and brain regions, a cell composition analysis, and gene-wise feature mappings. Our comprehensive molecular database facilitates the generation of data-driven hypotheses, enhancing our understanding of mammalian brain function during aging and disease. The data sets, along with an interactive database are freely available at https://www.ccb.uni-saarland.de/zebra.

Regulation of alternative splicing and polyadenylation in neurons

Cell-type-specific gene expression is a fundamental feature of multicellular organisms and is achieved by combinations of regulatory strategies. Although cell-restricted transcription is perhaps the most widely studied mechanism, co-transcriptional and post-transcriptional processes are also central to the spatiotemporal control of gene functions. One general category of expression control involves the generation of multiple transcript isoforms from an individual gene, whose balance and cell specificity are frequently tightly regulated via diverse strategies. The nervous system makes particularly extensive use of cell-specific isoforms, specializing the neural function of genes that are expressed more broadly. Here, we review regulatory strategies and RNA-binding proteins that direct neural-specific isoform processing. These include various classes of alternative splicing and alternative polyadenylation events, both of which broadly diversify the neural transcriptome. Importantly, global alterations of splicing and alternative polyadenylation are characteristic of many neural pathologies, and recent genetic studies demonstrate how misregulation of individual neural isoforms can directly cause mutant phenotypes.

Global detection of human variants and isoforms by deep proteome sequencing

An average shotgun proteomics experiment detects approximately 10,000 human proteins from a single sample. However, individual proteins are typically identified by peptide sequences representing a small fraction of their total amino acids. Hence, an average shotgun experiment fails to distinguish different protein variants and isoforms. Deeper proteome sequencing is therefore required for the global discovery of protein isoforms. Using six different human cell lines, six proteases, deep fractionation and three tandem mass spectrometry fragmentation methods, we identify a million unique peptides from 17,717 protein groups, with a median sequence coverage of approximately 80%. Direct comparison with RNA expression data provides evidence for the translation of most nonsynonymous variants. We have also hypothesized that undetected variants likely arise from mutation-induced protein instability. We further observe comparable detection rates for exon-exon junction peptides representing constitutive and alternative splicing events. Our dataset represents a resource for proteoform discovery and provides direct evidence that most frame-preserving alternatively spliced isoforms are translated.

Single-cell long-read sequencing in human cerebral organoids uncovers cell-type-specific and autism-associated exons

Dysregulation of alternative splicing has been repeatedly associated with neurodevelopmental disorders, but the extent of cell-type-specific splicing in human neural development remains largely uncharted. Here, single-cell long-read sequencing in induced pluripotent stem cell (iPSC)-derived cerebral organoids identifies over 31,000 uncatalogued isoforms and 4,531 cell-type-specific splicing events. Long reads uncover coordinated splicing and cell-type-specific intron retention events, which are challenging to study with short reads. Retained neuronal introns are enriched in RNA splicing regulators, showing shorter lengths, higher GC contents, and weaker 5' splice sites. We use this dataset to explore the biological processes underlying neurological disorders, focusing on autism. In comparison with prior transcriptomic data, we find that the splicing program in autistic brains is closer to the progenitor state than differentiated neurons. Furthermore, cell-type-specific exons harbor significantly more de novo mutations in autism probands than in siblings. Overall, these results highlight the importance of cell-type-specific splicing in autism and neuronal gene regulation.

Detection of isoforms and genomic alterations by high-throughput full-length single-cell RNA sequencing in ovarian cancer

Understanding the complex background of cancer requires genotype-phenotype information in single-cell resolution. Here, we perform long-read single-cell RNA sequencing (scRNA-seq) on clinical samples from three ovarian cancer patients presenting with omental metastasis and increase the PacBio sequencing depth to 12,000 reads per cell. Our approach captures 152,000 isoforms, of which over 52,000 were not previously reported. Isoform-level analysis accounting for non-coding isoforms reveals 20% overestimation of protein-coding gene expression on average. We also detect cell type-specific isoform and poly-adenylation site usage in tumor and mesothelial cells, and find that mesothelial cells transition into cancer-associated fibroblasts in the metastasis, partly through the TGF-β/miR-29/Collagen axis. Furthermore, we identify gene fusions, including an experimentally validated IGF2BP2::TESPA1 fusion, which is misclassified as high TESPA1 expression in matched short-read data, and call mutations confirmed by targeted NGS cancer gene panel results. With these findings, we envision long-read scRNA-seq to become increasingly relevant in oncology and personalized medicine.

Expression and splicing mediate distinct biological signals

BACKGROUND

Through alternative splicing, most human genes produce multiple isoforms in a cell-, tissue-, and disease-specific manner. Numerous studies show that alternative splicing is essential for development, diseases, and their treatments. Despite these important examples, the extent and biological relevance of splicing are currently unknown.

RESULTS

To solve this problem, we developed pairedGSEA and used it to profile transcriptional changes in 100 representative RNA-seq datasets. Our systematic analysis demonstrates that changes in splicing, on average, contribute to 48.1% of the biological signal in expression analyses. Gene-set enrichment analysis furthermore indicates that expression and splicing both convey shared and distinct biological signals.

CONCLUSIONS

These findings establish alternative splicing as a major regulator of the human condition and suggest that most contemporary RNA-seq studies likely miss out on critical biological insights. We anticipate our results will contribute to the transition from a gene-centric to an isoform-centric research paradigm.

Developmental isoform diversity in the human neocortex informs neuropsychiatric risk mechanisms

RNA splicing is highly prevalent in the brain and has strong links to neuropsychiatric disorders, yet the role of cell-type-specific splicing or transcript-isoform diversity during human brain development has not been systematically investigated. Here, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone (GZ) and cortical plate (CP) regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 unique isoforms, of which 72.6% are novel (unannotated in Gencode-v33), and uncovered a substantial contribution of transcript-isoform diversity, regulated by RNA binding proteins, in defining cellular identity in the developing neocortex. We leveraged this comprehensive isoform-centric gene annotation to re-prioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders.

One-Sentence Summary

A cell-specific atlas of gene isoform expression helps shape our understanding of brain development and disease.

Structured Abstract

INTRODUCTION: The development of the human brain is regulated by precise molecular and genetic mechanisms driving spatio-temporal and cell-type-specific transcript expression programs. Alternative splicing, a major mechanism increasing transcript diversity, is highly prevalent in the human brain, influences many aspects of brain development, and has strong links to neuropsychiatric disorders. Despite this, the cell-type-specific transcript-isoform diversity of the developing human brain has not been systematically investigated.RATIONALE: Understanding splicing patterns and isoform diversity across the developing neocortex has translational relevance and can elucidate genetic risk mechanisms in neurodevelopmental disorders. However, short-read sequencing, the prevalent technology for transcriptome profiling, is not well suited to capturing alternative splicing and isoform diversity. To address this, we employed third-generation long-read sequencing, which enables capture and sequencing of complete individual RNA molecules, to deeply profile the full-length transcriptome of the germinal zone (GZ) and cortical plate (CP) regions of the developing human neocortex at tissue and single-cell resolution.RESULTS: We profiled microdissected GZ and CP regions of post-conception week (PCW) 15-17 human neocortex in bulk and at single-cell resolution across six subjects using high-fidelity long-read sequencing (PacBio IsoSeq). We identified 214,516 unique isoforms, of which 72.6% were novel (unannotated in Gencode), and >7,000 novel exons, expanding the proteome by 92,422 putative proteoforms. We uncovered thousands of isoform switches during cortical neurogenesis predicted to impact RNA regulatory domains or protein structure and implicating previously uncharacterized RNA-binding proteins in cellular identity and neuropsychiatric disease. At the single-cell level, early-stage excitatory neurons exhibited the greatest isoform diversity, and isoform-centric single-cell clustering led to the identification of previously uncharacterized cell states. We systematically assessed the contribution of transcriptomic features, and localized cell and spatio-temporal transcript expression signatures across neuropsychiatric disorders, revealing predominant enrichments in dynamic isoform expression and utilization patterns and that the number and complexity of isoforms per gene is strongly predictive of disease. Leveraging this resource, we re-prioritized thousands of rare de novo risk variants associated with autism spectrum disorders (ASD), intellectual disability (ID), and neurodevelopmental disorders (NDDs), more broadly, to potentially more severe consequences and revealed a larger proportion of cryptic splice variants with the expanded transcriptome annotation provided in this study.CONCLUSION: Our study offers a comprehensive landscape of isoform diversity in the human neocortex during development. This extensive cataloging of novel isoforms and splicing events sheds light on the underlying mechanisms of neurodevelopmental disorders and presents an opportunity to explore rare genetic variants linked to these conditions. The implications of our findings extend beyond fundamental neuroscience, as they provide crucial insights into the molecular basis of developmental brain disorders and pave the way for targeted therapeutic interventions. To facilitate exploration of this dataset we developed an online portal ( https://sciso.gandallab.org/ ).

Analyzing alternative splicing in Alzheimer's disease postmortem brain: a cell-level perspective

Alzheimer's disease (AD) is a neurodegenerative disease with no effective cure that attacks the brain's cells resulting in memory loss and changes in behavior and language skills. Alternative splicing is a highly regulated process influenced by specific cell types and has been implicated in age-related disorders such as neurodegenerative diseases. A comprehensive detection of alternative splicing events (ASEs) at the cellular level in postmortem brain tissue can provide valuable insights into AD pathology. Here, we provided cell-level ASEs in postmortem brain tissue by employing bioinformatics pipelines on a bulk RNA sequencing study sorted by cell types and two single-cell RNA sequencing studies from the prefrontal cortex. This comprehensive analysis revealed previously overlooked splicing and expression changes in AD patient brains. Among the observed alterations were changed in the splicing and expression of transcripts associated with chaperones, including CLU in astrocytes and excitatory neurons, PTGDS in astrocytes and endothelial cells, and HSP90AA1 in microglia and tauopathy-afflicted neurons, which were associated with differential expression of the splicing factor DDX5. In addition, novel, unknown transcripts were altered, and structural changes were observed in lncRNAs such as MEG3 in neurons. This work provides a novel strategy to identify the notable ASEs at the cell level in neurodegeneration, which revealed cell type-specific splicing changes in AD. This finding may contribute to interpreting associations between splicing and neurodegenerative disease outcomes.

Coordination of alternative splicing and alternative polyadenylation revealed by targeted long read sequencing

Nervous system development is associated with extensive regulation of alternative splicing (AS) and alternative polyadenylation (APA). AS and APA have been extensively studied in isolation, but little is known about how these processes are coordinated. Here, the coordination of cassette exon (CE) splicing and APA in Drosophila was investigated using a targeted long-read sequencing approach we call Pull-a-Long-Seq (PL-Seq). This cost-effective method uses cDNA pulldown and Nanopore sequencing combined with an analysis pipeline to quantify inclusion of alternative exons in connection with alternative 3' ends. Using PL-Seq, we identified genes that exhibit significant differences in CE splicing depending on connectivity to short versus long 3'UTRs. Genomic long 3'UTR deletion was found to alter upstream CE splicing in short 3'UTR isoforms and ELAV loss differentially affected CE splicing depending on connectivity to alternative 3'UTRs. This work highlights the importance of considering connectivity to alternative 3'UTRs when monitoring AS events.