Estimation of alternative splicing isoform frequencies from RNA-Seq data

Area of Science:

• Genomics
• Bioinformatics
• Molecular Biology

Background:

• RNA sequencing (RNA-Seq) is a leading technology for gene expression profiling.
• Estimating expression levels of alternative splicing gene isoforms is challenging due to short read lengths in current RNA-Seq technologies.

Purpose of the Study:

• To develop a novel algorithm for accurate inference of isoform- and gene-specific expression levels from RNA-Seq data.
• To address the challenges in gene expression profiling posed by short sequencing reads.

Main Methods:

• Developed IsoEM, an expectation-maximization algorithm for RNA-Seq data analysis.
• IsoEM utilizes insert size distribution, base quality scores, strand, and read pairing information.
• An open-source Java implementation of IsoEM is available.

Main Results:

• IsoEM demonstrates scalable running time and superior performance compared to existing methods for isoform and gene expression level estimation.
• Empirical experiments on synthetic and real RNA-Seq datasets validate IsoEM's effectiveness.
• Simulation results indicate that longer reads (>25-36 bases) do not always improve accuracy for isoform and gene expression estimation at a fixed sequencing cost.

Conclusions:

• IsoEM provides a robust and efficient solution for gene and isoform expression level estimation from RNA-Seq data.
• The findings suggest optimal read length strategies for cost-effective and accurate gene expression profiling.