GaussF: An Unbiased Gaussian Framework for High-Precision Transcript Quantification with a Programmable Self-Benchmark Model
- 1Shenzhen university
Abstract
Short-read RNA sequencing has revolutionized transcriptome analysis, but accurate isoform-level quantification faces substantial challenges from overlapping sequences among complex isoforms and inherent sequencing biases, leading to uneven coverage. Current alignment-based and alignment-free methods, while valuable, largely rely on empirical, data-derived factors for bias correction and isoform resolution. This dependency limits their capacity to fully capture the intricate spectrum of sequencing data. We introduce GaussF, a novel alignment-free methodology specifically designed for precise isoform quantification by directly addressing these biases and the isoform recognition problem without relying on empirical data. GaussF is predicated on the observation that k-mer counts categorized by GC content across a natural transcript exhibit an intrinsic Gaussian distribution. It establishes a transcript-specific self-benchmarking workflow utilizing a theoretical Gaussian cumulative distribution function (CDF). Crucially, the parameters defining this CDF are derived independently of any empirical sequencing data. This parameter-fixed theoretical CDF serves as an unbiased benchmark, capturing only the genuine contribution of abundance while systematically excluding all bias-related data. To improve isoform resolution despite shared regions, GaussF characterizes isoforms through identifying minimal-shared regions and restructures the k-mer reference accordingly. This approach simplifies complexity and enhances accuracy. Rigorous validation using artificial and real RNA samples, alongside comparative analysis with leading quantification tools like Kallisto, Salmon, Cufflinks, StringTie, and RSEM, demonstrates GaussF's superior performance in accurate isoform detection and quantification. GaussF is available on GitHub and PyPI, representing a significant advancement in RNA-seq analysis.
Related
- Genetically programmable molecular framework for constructing precision recombinant human ECM mimics
- SEME enables programmable transcript reconstruction through synthetic microexon incorporation
- Bramble: projection of spliced genomic alignments into transcriptomic space for improved transcript quantification
- Transcript quantification using context-aware MAPQ modeling
- Systematic evaluation of transcript normalization methods and quantification pipelines across long-read sequencing platforms
BibTeX
@misc{su2026gaussf,
title = {{GaussF: An Unbiased Gaussian Framework for High-Precision Transcript Quantification with a Programmable Self-Benchmark Model}},
author = {Q. Su},
howpublished = {bioRxiv},
year = {2026},
month = oct,
doi = {10.64898/2026.10.28.762697},
url = {https://www.biorxiv.org/content/10.64898/2026.10.28.762697v1},
note = {CC BY-NC-ND},
}
RIS
TY - UNPB
TI - GaussF: An Unbiased Gaussian Framework for High-Precision Transcript Quantification with a Programmable Self-Benchmark Model
AU - Q. Su
PY - 2026
DA - 2026/10/28
PB - bioRxiv
DO - 10.64898/2026.10.28.762697
UR - https://www.biorxiv.org/content/10.64898/2026.10.28.762697v1
AB - Short-read RNA sequencing has revolutionized transcriptome analysis, but accurate isoform-level quantification faces substantial challenges from overlapping sequences among complex isoforms and inherent sequencing biases, leading to uneven coverage. Current alignment-based and alignment-free methods, while valuable, largely rely on empirical, data-derived factors for bias correction and isoform resolution. This dependency limits their capacity to fully capture the intricate spectrum of sequencing data. We introduce GaussF, a novel alignment-free methodology specifically designed for precise isoform quantification by directly addressing these biases and the isoform recognition problem without relying on empirical data. GaussF is predicated on the observation that k-mer counts categorized by GC content across a natural transcript exhibit an intrinsic Gaussian distribution. It establishes a transcript-specific self-benchmarking workflow utilizing a theoretical Gaussian cumulative distribution function (CDF). Crucially, the parameters defining this CDF are derived independently of any empirical sequencing data. This parameter-fixed theoretical CDF serves as an unbiased benchmark, capturing only the genuine contribution of abundance while systematically excluding all bias-related data. To improve isoform resolution despite shared regions, GaussF characterizes isoforms through identifying minimal-shared regions and restructures the k-mer reference accordingly. This approach simplifies complexity and enhances accuracy. Rigorous validation using artificial and real RNA samples, alongside comparative analysis with leading quantification tools like Kallisto, Salmon, Cufflinks, StringTie, and RSEM, demonstrates GaussF's superior performance in accurate isoform detection and quantification. GaussF is available on GitHub and PyPI, representing a significant advancement in RNA-seq analysis.
ER -