Explore the vast range of titles available.

RNA sequencing studies with complex designs and transcript-resolution analyses involve multiple hypotheses per gene; however, conventional approaches fail to control the false discovery rate (FDR) at gene level. We propose stageR, a two-stage testing paradigm that leverages the increased power of aggregated gene-level tests and allows post hoc assessment for significant genes. This method provides gene-level FDR control and boosts power for testing interaction effects. In transcript-level analysis, it provides a framework that performs powerful gene-level tests while maintaining biological interpretation at transcript-level resolution. The procedure is applicable whenever individual hypotheses can be aggregated, providing a unified framework for complex high-throughput experiments.

Alternative splicing produces multiple functional transcripts from a single gene. Dysregulation of splicing is known to be associated with disease and as a hallmark of cancer. Existing tools for differential transcript usage (DTU) analysis either lack in performance, cannot account for complex experimental designs or do not scale to massive scRNA-seq data. We introduce satuRn, a fast and flexible quasi-binomial generalized linear modelling framework that is on par with the best performing DTU methods from the bulk RNA-seq realm, while providing good false discovery rate control, addressing complex experimental designs and scaling to scRNA-seq applications.

Alternative splicing (AS) is a fundamental mechanism governing transcriptomic diversity and cellular identity. Although 293T (human embryonic kidney) and A549 (human lung adenocarcinoma) cell lines are widely used, cell-type-specific splicing dynamics—including responses to receptor overexpression—remain incompletely characterized. To address this, we integrated Oxford Nanopore long-read sequencing with BGI short-read data to profile transcriptomes under both basal and GPCR-overexpressing conditions (ADORA3 in 293T; P2RY12 in A549). Full-length isoform analysis using FLAIR and SQANTI3 revealed extensive transcriptomic complexity, including 18.02% novel isoforms in 293T and 19.52% in A549 cells. We found that 293T cells exhibited a stable transcriptome architecture enriched in splicing-related pathways, whereas A549 cells underwent broader transcriptional remodeling linked to tumorigenic processes. These findings suggest that 293T cells may be a suitable model for investigating splicing regulation, while A549 cells could serve as a relevant system for exploring tumor-related transcriptome dynamics. Our work elucidates context-dependent AS regulation and underscores the value of integrating long-read sequencing with FLAIR/SQANTI3 for dissecting cell-state-specific transcriptome dynamics.

Over the past 8 years, multiple studies examined the phenomenon of isoform switching in human cancers and discovered that isoform switching is widespread, with hundreds to thousands of such events per cancer type. Although all of these studies used slightly different definitions of isoform switching, which in part led to a rather poor overlap of their results, they all leveraged transcript usage, a proportion of the transcript’s expression in the total expression level of the parent gene, to detect isoform switching. However, how changes in transcript usage correlate with changes in transcript expression is not sufficiently explored. In this article, we adopt the most common definition of isoform switching and use a state-of-the-art tool for the analysis of differential transcript usage, SatuRn, to detect isoform switching events in 12 cancer types. We analyze the detected events in terms of changes in transcript usage and the relationship between transcript usage and transcript expression on a global scale. The results of our analysis suggest that the relationship between changes in transcript usage and changes in transcript expression is far from straightforward, and that such quantitative information can be effectively used for prioritizing isoform switching events for downstream analyses.

Detection of differential transcript usage (DTU) from RNA-seq data is an important bioinformatic analysis that complements differential gene expression analysis. Here we present a simple workflow using a set of existing R/Bioconductor packages for analysis of DTU. We show how these packages can be used downstream of RNA-seq quantification using the Salmon software package. The entire pipeline is fast, benefiting from inference steps by Salmon to quantify expression at the transcript level. The workflow includes live, runnable code chunks for analysis using DRIMSeq and DEXSeq, as well as for performing two-stage testing of DTU using the stageR package, a statistical framework to screen at the gene level and then confirm which transcripts within the significant genes show evidence of DTU. We evaluate these packages and other related packages on a simulated dataset with parameters estimated from real data.

Background RNA sequencing of whole blood has been increasingly employed to find transcriptomic signatures of disease states. These studies traditionally utilize short-read sequencing of cDNA, missing important aspects of RNA expression such as differential isoform abundance and poly(A) tail length variation. Methods We used Oxford Nanopore Technologies sequencing to sequence native mRNA extracted from whole blood from 12 patients with definite bacterial and viral sepsis and compared with results from matching Illumina short-read cDNA sequencing data. Additionally, we explored poly(A) tail length variation, novel transcript identification, and differential transcript usage. Results The correlation of gene count data between Illumina cDNA- and Nanopore RNA-sequencing strongly depended on the choice of analysis pipeline; NanoCount for Nanopore and Kallisto for Illumina data yielded the highest mean Pearson’s correlation of 0.927 at the gene level and 0.736 at the transcript isoform level. We identified 2 genes with differential polyadenylation, 9 genes with differential expression and 4 genes with differential transcript usage between bacterial and viral infection. Gene ontology gene set enrichment analysis of poly(A) tail length revealed enrichment of long tails in mRNA of genes involved in signaling and short tails in oxidoreductase molecular functions. Additionally, we detected 240 non-artifactual novel transcript isoforms. Conclusions Nanopore RNA- and Illumina cDNA-gene counts are strongly correlated, indicating that both platforms are suitable for discovery and validation of gene count biomarkers. Nanopore direct RNA-seq provides additional advantages by uncovering additional post- and co-transcriptional biomarkers, such as poly(A) tail length variation and transcript isoform usage.

Alternative splicing is a biological process during gene expression that allows a single gene to code for multiple proteins. However, splicing patterns can be altered in some conditions or diseases. Here, we present BANDITS, a R/Bioconductor package to perform differential splicing, at both gene and transcript level, based on RNA-seq data. BANDITS uses a Bayesian hierarchical structure to explicitly model the variability between samples and treats the transcript allocation of reads as latent variables. We perform an extensive benchmark across both simulated and experimental RNA-seq datasets, where BANDITS has extremely favourable performance with respect to the competitors considered.

For steady transgenic expression, lentiviral vector-mediated gene delivery is a commonly used technique. One question that needs to be explored is how external lentiviral vectors and overexpressed genes perturb cellular homeostasis, potentially altering transcriptional networks. In this study, two Human Embryonic Kidney 293T (HEK293T)-derived cell lines were established via lentiviral transduction, one overexpressing green fluorescent protein (GFP) and the other co-overexpressing GFP and ADORA3 following puromycin selection to ensure stable genomic integration. Genes with differentially transcript utilization (gDTUs) and differentially expressed genes (DEGs) across cell lines were identified after short-read and long-read RNA-seq. Only 31 genes were discovered to have changed in expression when GFP was expressed, although hundreds of genes showed variations in transcript use. In contrast, even when co-overexpression of GFP and ADORA3 alters the expression of more than 1000 genes, there are still less than 1000 gDTUs. Moreover, DEGs linked to ADORA3 overexpression play a major role in RNA splicing, whereas gDTUs are highly linked to a number of malignancies and the molecular mechanisms that underlie them. For the analysis of gene expression data from stable cell lines derived from HEK293T, our findings provide important insights into changes in gene expression and alternative splicing.

Dysregulated RNA splicing is a well-recognized characteristic of colorectal cancer (CRC); however, its intricacies remain obscure, partly due to challenges in profiling full-length transcript variants at the single-cell level. Here, we employ high-depth long-read scRNA-seq to define the full-length transcriptome of colorectal epithelial cells in 12 CRC patients, revealing extensive isoform diversities and splicing alterations. Cancer cells exhibited increased transcript complexity, with widespread 3′-UTR shortening and reduced intron retention. Distinct splicing regulation patterns were observed between intrinsic-consensus molecular subtypes (iCMS), with iCMS3 displaying even higher splicing factor activities and more pronounced 3′-UTR shortening. Furthermore, we revealed substantial shifts in isoform usage that result in alterations of protein sequences from the same gene with distinct carcinogenic effects during tumorigenesis of CRC. Allele-specific expression analysis revealed dominant mutant allele expression in key oncogenes and tumor suppressors. Moreover, mutated PPIG was linked to widespread splicing dysregulation, and functional validation experiments confirmed its critical role in modulating RNA splicing and tumor-associated processes. Our findings highlight the transcriptomic plasticity in CRC and suggest novel candidate targets for splicing-based therapeutic strategies.

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.