Aryal, Sandeep

Abstract Recent advances in long-read transcriptome sequencing enable high-throughput profiling of full-length RNA isoforms in bulk, single-cell, and single-nucleus samples. However, long-read datasets typically contain a mixture of complete and partial transcripts, leading to pervasive ambiguity in read-to-isoform assignment and complicating accurate isoform identification and quantification, particularly in the absence of reliable reference annotations. These challenges are further amplified in single-cell and single-nucleus samples, where coverage is sparse and transcriptional heterogeneity is high. Here, we present the Long Read Alignment Assembler (LRAA), a unified and versatile computational framework for isoform identification and quantification from long-read RNA sequencing data across bulk, single-cell, and single-nucleus transcriptomic samples. LRAA combines splice-graph based structural modeling with expectation maximization based optimization to probabilistically resolve ambiguous read assignments and improve isoform abundance estimation. The framework supports quantification-only, reference-guided, and fully reference-free (de novo) modes of analysis within a single methodological paradigm. We benchmarked LRAA using both simulated and genuine long-read datasets spanning sequencing standards and whole transcriptomes. Central to this evaluation is a novel benchmarking strategy based on Multiplexed Overexpression of Regulatory Factors (MORFs), which provides biologically expressed, barcoded isoforms with unambiguous read-level ground truth. Across all benchmarks, including MORFs, synthetic spike-ins, and whole-transcriptome datasets, LRAA consistently outperformed state-of-the-art methods in isoform identification accuracy, sensitivity, and expression quantification. Finally, we demonstrate the biological utility of LRAA by resolving cell-type-specific isoform usage across peripheral blood immune cell populations and by detecting a pathogenic cryptic isoform of STMN2 with associated transcriptional changes in single-nucleus RNA-seq data from frontal cortex tissue of an individual with frontotemporal dementia (FTD). Together, these results establish LRAA as a robust and general solution for resolving transcript diversity in complex biological systems, from development to disease.

GGGGCC (G 4 C 2 ) repeat expansion in C9ORF72 is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Toxicity is thought to result from the accumulation of either repeat RNAs and/or dipeptide repeat proteins (DPRs) translated from repeat-containing transcripts through repeat-associated non-AUG (RAN) translation. To disentangle RNA from DPR toxicity, we mutated a CUG codon predominantly used to initiate DPR translation from all three reading frames. This mutation disrupted DPR synthesis while preserving the expression of repeat-containing RNAs. Despite the accumulation of RNA foci, behavioral deficits and pathological abnormalities, including p-TDP-43 inclusions, STING activation, motor neuron loss, neuroinflammation, and increased plasma neurofilament concentration, were alleviated in C9ORF72 mice. Base editing of the CUG codon also improved molecular phenotypes and survival in patient induced pluripotent stem cell–derived neurons, which highlights the potential of therapeutically targeting DPR production rather than repeat RNAs.