From words to complete phrases: insight into single-cell isoforms using short and long reads

Single-cell Rapid Capture Hybridization sequencing reliably detects isoform usage and coding mutations in targeted genes

Single-cell long-read sequencing has transformed our understanding of isoform usage and the mutation heterogeneity between cells. Despite unbiased in-depth analysis, the low sequencing throughput often results in insufficient read coverage, thereby limiting our ability to perform mutation calling for specific genes. Here, we developed a single-cell Rapid Capture Hybridization sequencing (scRaCH-seq) method that demonstrates high specificity and efficiency in capturing targeted transcripts using long-read sequencing, allowing an in-depth analysis of mutation status and transcript usage for genes of interest. The method includes creating a probe panel for transcript capture, using barcoded primers for pooling and efficient sequencing via Oxford Nanopore Technologies platforms. scRaCH-seq is applicable to stored and indexed single-cell cDNA, which allows analysis to be combined with existing short-read RNA-seq data sets. In our investigation of BTK and SF3B1 genes in samples from patients with chronic lymphocytic leukemia (CLL), we detect SF3B1 isoforms and mutations with high sensitivity. Integration with short-read single-cell RNA sequencing (scRNA-seq) data reveals significant gene expression differences in SF3B1-mutated CLL cells, although it does not impact the sensitivity of the anticancer drug venetoclax. scRaCH-seq's capability to study long-read transcripts of multiple genes makes it a powerful tool for single-cell genomics.

De novo detection of somatic variants in high-quality long-read single-cell RNA sequencing data

In cancer, genetic and transcriptomic variations generate clonal heterogeneity, leading to treatment resistance. Long-read single-cell RNA sequencing (LR scRNA-seq) has the potential to detect genetic and transcriptomic variations simultaneously. Here, we present LongSom, a computational workflow leveraging high-quality LR scRNA-seq data to call de novo somatic single-nucleotide variants (SNVs), including in mitochondria (mtSNVs), copy number alterations (CNAs), and gene fusions, to reconstruct the tumor clonal heterogeneity. Before somatic variant calling, LongSom reannotates marker gene-based cell types using cell mutational profiles. LongSom distinguishes somatic SNVs from noise and germline polymorphisms by applying an extensive set of hard filters and statistical tests. Applying LongSom to human ovarian cancer samples, we detected clinically relevant somatic SNVs that were validated against matched DNA samples. Leveraging somatic SNVs and fusions, LongSom found subclones with different predicted treatment outcomes. In summary, LongSom enables de novo variant detection without the need for normal samples, facilitating the study of cancer evolution, clonal heterogeneity, and treatment resistance.

Transcriptomics in the era of long-read sequencing

Transcriptome sequencing revolutionized the analysis of gene expression, providing an unbiased approach to gene detection and quantification that enabled the discovery of novel isoforms, alternative splicing events and fusion transcripts. However, although short-read sequencing technologies have surpassed the limited dynamic range of previous technologies such as microarrays, they have limitations, for example, in resolving full-length transcripts and complex isoforms. Over the past 5 years, long-read sequencing technologies have matured considerably, with improvements in instrumentation and analytical methods, enabling their application to RNA sequencing (RNA-seq). Benchmarking studies are beginning to identify the strengths and limitations of long-read RNA-seq, although there remains a need for comprehensive resources to guide newcomers through the intricacies of this approach. In this Review, we provide a comprehensive overview of the long-read RNA-seq workflow, from library preparation and sequencing challenges to core data processing, downstream analyses and emerging developments. We present an extensive inventory of experimental and analytical methods and discuss current challenges and prospects.

AsaruSim: a single-cell and spatial RNA-Seq Nanopore long-reads simulation workflow

MOTIVATION

The combination of long-read sequencing technologies like Oxford Nanopore with single-cell RNA sequencing (scRNAseq) assays enables the detailed exploration of transcriptomic complexity, including isoform detection and quantification, by capturing full-length cDNAs. However, challenges remain, including the lack of advanced simulation tools that can effectively mimic the unique complexities of scRNAseq long-read datasets. Such tools are essential for the evaluation and optimization of isoform detection methods dedicated to single-cell long-read studies.

RESULTS

We developed AsaruSim, a workflow that simulates synthetic single-cell long-read Nanopore datasets, closely mimicking real experimental data. AsaruSim employs a multi-step process that includes the creation of a synthetic count matrix, generation of perfect reads, optional PCR amplification, introduction of sequencing errors, and comprehensive quality control reporting. Applied to a dataset of human peripheral blood mononuclear cells, AsaruSim accurately reproduced experimental read characteristics.

AVAILABILITY AND IMPLEMENTATION

The source code and full documentation are available at https://github.com/GenomiqueENS/AsaruSim.

Interneuron FGF13 regulates seizure susceptibility via a sodium channel-independent mechanism

Developmental and epileptic encephalopathies (DEEs), a class of devastating neurological disorders characterized by recurrent seizures and exacerbated by disruptions to excitatory/inhibitory balance in the brain, are commonly caused by mutations in ion channels. Disruption of, or variants in, FGF13 were implicated as causal for a set of DEEs, but the underlying mechanisms were clouded because FGF13 is expressed in both excitatory and inhibitory neurons, FGF13 undergoes extensive alternative splicing producing multiple isoforms with distinct functions, and the overall roles of FGF13 in neurons are incompletely cataloged. To overcome these challenges, we generated a set of novel cell-type-specific conditional knockout mice. Interneuron-targeted deletion of Fgf13 led to perinatal mortality associated with extensive seizures and impaired the hippocampal inhibitory/excitatory balance while excitatory neuron-targeted deletion of Fgf13 caused no detectable seizures and no survival deficits. While best studied as a voltage-gated sodium channel (Nav) regulator, we observed no effect of Fgf13 ablation in interneurons on Navs but rather a marked reduction in K+ channel currents. Re-expressing different Fgf13 splice isoforms could partially rescue deficits in interneuron excitability and restore K+ channel current amplitude. These results enhance our understanding of the molecular mechanisms that drive the pathogenesis of Fgf13-related seizures and expand our understanding of FGF13 functions in different neuron subsets.

Profiling the epigenome using long-read sequencing

The advent of single-molecule, long-read sequencing (LRS) technologies by Oxford Nanopore Technologies and Pacific Biosciences has revolutionized genomics, transcriptomics and, more recently, epigenomics research. These technologies offer distinct advantages, including the direct detection of methylated DNA and simultaneous assessment of DNA sequences spanning multiple kilobases along with their modifications at the single-molecule level. This has enabled the development of new assays for analyzing chromatin states and made it possible to integrate data for DNA methylation, chromatin accessibility, transcription factor binding and histone modifications, thereby facilitating comprehensive epigenomic profiling. Owing to recent advancements, alternative, nascent and translating transcripts can be detected using LRS approaches. This Review discusses LRS-based experimental and computational strategies for characterizing chromatin states and highlights their advantages over short-read sequencing methods. Furthermore, we demonstrate how various long-read methods can be integrated to design multi-omics studies to investigate the relationship between chromatin states and transcriptional dynamics.

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.

De novo detection of somatic variants in high-quality long-read single-cell RNA sequencing data

In cancer, genetic and transcriptomic variations generate clonal heterogeneity, leading to treatment resistance. Long-read single-cell RNA sequencing (LR scRNA-seq) has the potential to detect genetic and transcriptomic variations simultaneously. Here, we present LongSom, a computational workflow leveraging high-quality LR scRNA-seq data to call de novo somatic single-nucleotide variants (SNVs), including in mitochondria (mtSNVs), copy-number alterations (CNAs), and gene fusions, to reconstruct the tumor clonal heterogeneity. Before somatic variants calling, LongSom re-annotates marker gene based cell types using cell mutational profiles. LongSom distinguishes somatic SNVs from noise and germline polymorphisms by applying an extensive set of hard filters and statistical tests. Applying LongSom to human ovarian cancer samples, we detected clinically relevant somatic SNVs that were validated against matched DNA samples. Leveraging somatic SNVs and fusions, LongSom found subclones with different predicted treatment outcomes. In summary, LongSom enables de novo variants detection without the need for normal samples, facilitating the study of cancer evolution, clonal heterogeneity, and treatment resistance.

Interneuron FGF13 regulates seizure susceptibility via a sodium channel-independent mechanism

Developmental and Epileptic Encephalopathies (DEEs), a class of devastating neurological disorders characterized by recurrent seizures and exacerbated by disruptions to excitatory/inhibitory balance in the brain, are commonly caused by mutations in ion channels. Disruption of, or variants in, FGF13 were implicated as causal for a set of DEEs, but the underlying mechanisms were clouded because FGF13 is expressed in both excitatory and inhibitory neurons, FGF13 undergoes extensive alternative splicing producing multiple isoforms with distinct functions, and the overall roles of FGF13 in neurons are incompletely cataloged. To overcome these challenges, we generated a set of novel cell type-specific conditional knockout mice. Interneuron-targeted deletion of Fgf13 led to perinatal mortality associated with extensive seizures and impaired the hippocampal inhibitory/excitatory balance while excitatory neuron-targeted deletion of Fgf13 caused no detectable seizures and no survival deficits. While best studied as a voltage-gated sodium channel (Nav) regulator, we observed no effect of Fgf13 ablation in interneurons on Navs but rather a marked reduction in K+ channel currents. Re-expressing different Fgf13 splice isoforms could partially rescue deficits in interneuron excitability and restore K+ channel current amplitude. These results enhance our understanding of the molecular mechanisms that drive the pathogenesis of Fgf13-related seizures and expand our understanding of FGF13 functions in different neuron subsets.

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.

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.