2025 Hyper Recent •CC0 1.0 Universal

This work is dedicated to the public domain. No rights reserved.

Access Preprint From Server
April 13th, 2025
Version: 1
Walter and Eliza Hall Institute of Medical Research
bioinformatics
biorxiv

Dividing out quantification uncertainty enables assessment of differential transcript usage with diffSplice

Baldoni, P. L.Open in Google Scholar•Chen, L.Open in Google Scholar•Li, M.Open in Google Scholar•Chen, Y.Open in Google Scholar•Smyth, G. K.Open in Google Scholar

Differential transcript usage (DTU) refers to changes in the relative abundance of transcript isoforms of the same gene between experimental conditions, even when the total expression of the gene doesn\'t change. DTU analysis requires the quantification of individual isoforms from RNA-seq data, which has a high level of uncertainty due to transcript overlap and read-to-transcript ambiguity (RTA). Popular DTU analysis methods do not directly account for the RTA overdispersion within their statistical frameworks, leading to reduced statistical power or poor error rate control, particularly in scenarios with small sample sizes. This article presents limma and edgeR analysis pipelines that account for RTA during DTU assessment. Leveraging recent advancements in the limma and edgeR Bioconductor packages, we propose DTU analysis pipelines optimized for small and large datasets with a unified interface via the diffSplice function. The pipelines make use of divided counts to remove RTA-induced dispersion from transcript isoform counts and account for the sparsity in transcript-level counts. Simulations and analysis of real data from mouse mammary epithelial cells demonstrate that the diffSplice pipelines provide greater power, improved efficiency, and improved FDR control compared to existing specialized DTU methods.

Similar Papers

biorxiv
Tue Apr 15 2025
Hepatic isomiR landscaping reveals new biological insights into metabolic dysfunction in steatotic liver disease
Post-transcriptionally modified microRNA (miRNA), called isomiRs, expand the repertoire of transcripts that can leveraged for therapeutic targets and biological insights. However, the expression of isomiRs has not been characterized in metabolic dysfunction-associated steatotic liver disease (MASLD). Therefore, we assessed the isomiR expression profile in liver biopsies from 79 patients with MASLD...
Brion, C.
•
Hoang, S. A.
•
Wang, G.
•
Mirshahi, F.
...•
Sanyal, A. J.
biorxiv
Tue Apr 15 2025
RNAchat: Integrating machine learning algorithms to identify metapathways based on clinical and multi-omics data
Background: Understanding the global interactions among cellular pathways and types is essential for unraveling the complex biological mechanisms that underlie various diseases. While many existing studies have primarily focused on omics data, they often fail to integrate clinical heterogeneity, which is crucial for a comprehensive understanding of disease biology. Although significant progress ha...
Tang, M.
biorxiv
Tue Apr 15 2025
ANOMALY: A Snakemake pipeline for identifying NuMTs from Long-Read Sequencing Data
Motivation Nuclear mitochondrial DNA segments (NuMT) can significantly affect cellular processes, including cancer development and disease progression. Current methods to call NuMTs rely on short-read sequencing data but struggle to resolve complex NuMTs. These limitations can be overcome by employing long-read sequencing data. However, no such workflow exists to capture NuMTs from long-read seque...
Mahar, N. S.
•
Singh, R.
•
Gupta, I.
•
Ramdas, S.
biorxiv
Tue Apr 15 2025
CellLoop: Identifying single-cell 3D genome chromatin loops
Single-cell 3D genome mapping offers insights into chromatin architecture, yet the inherent sparsity and high dimensionality of contact maps hinder reliable identification of chromatin loops at the single-cell level. Here, we present CellLoop, a single-cell chromatin loop detection algorithm based on a density-based center detection framework that integrates intra-cellular and neighboring inter-ce...
Ye, Y.
•
Hu, Y.
•
Zhang, S.
•
Chen, H.
•
Gao, L.
biorxiv
Tue Apr 15 2025
BINSEQ: A Family of High-Performance Binary Formats for Nucleotide Sequences
Modern genomics produces billions of sequencing records per run, which are typically stored as gzip-compressed FASTQ files. While this format is widely used, it is not optimal for high-throughput processing due to its reliance on single-threaded decompression and sequential parsing of irregularly sized records. This limitation is particularly problematic for applications that would benefit from pa...
Teyssier, N.
•
Dobin, A.
biorxiv
Tue Apr 15 2025
ganon2: up-to-date and scalable metagenomics analysis
The fast growth of public genomic sequence repositories greatly contributes to the success of metagenomics. However, they are growing at a faster pace than the computational resources to use them. This challenges current methods, which struggle to take full advantage of massive and fast data generation. We propose a generational leap in performance and usability with ganon2, a sequence classificat...
Piro, V. C.
•
Reinert, K.
biorxiv
Tue Apr 15 2025
Circuit Mining in Transcriptomics Data
A central goal in neuropharmacological research is to alter brain function by targeting genes whose expression is specific to the corresponding brain circuit. Identifying such genes in large spatially resolved transcriptomics data requires the expertise of bioinformaticians for handling data complexity and to perform statistical tests. This time-consuming process is often decoupled from the routin...
Peherstorfer, T.
•
Ulonska, S.
•
Burger, B.
•
Lucato, S.
...•
Bühler, K.
biorxiv
Tue Apr 15 2025
Comparison and Optimization of Cellular Neighbor Preference Methods for Quantitative Tissue Analysis
Studying the spatial distribution of cell types in tissues is essential for understanding their function in health and disease. A widely used spatial feature for quantifying tissue organization is the pairwise neighbor preference (NEP) of cell types, commonly referred to as co-occurrence or colocalization. Various methods to infer NEP have proved their utility in spatial omics studies, but despite...
Schiller, C.
•
Ibarra-Arellano, M. A.
•
Bestak, K.
•
Tanevski, J.
•
Schapiro, D.
biorxiv
Tue Apr 15 2025
FLASH-MM: fast and scalable single-cell differential expression analysis using linear mixed-effects models
Single-cell RNA sequencing (scRNA-seq) enables detailed comparisons of gene expression across cells and conditions. Single-cell differential expression analysis faces challenges like sample correlation, individual variation, and scalability. We developed a fast and scalable linear mixed-effects model (LMM) estimation algorithm, FLASH-MM, to address these issues. We reformulate aspects of the model...
Xu, C.
•
Pouyabahar, D.
•
Voisin, V.
•
Heydari, H.
•
Bader, G.
biorxiv
Tue Apr 15 2025
Male XIST expression in cardiac pseudo-glia does not induce X chromosome inactivation
XIST is a long non-coding RNA that mediates the process of X chromosome inactivation in females, and has not been previously observed in males outside of developmental, pathological, and ectopic contexts. We report the unprecedented endogenous expression of XIST in 204 male heart samples across 14 single-cell RNA sequencing studies. The male expression of XIST is specific to a rare cell type typic...
Gorin, G.
•
Goodman, L.