CORNAS: coverage-dependent RNA-Seq analysis of gene expression data without biological replicates.

Area of Science:

• Bioinformatics
• Computational Biology
• Genomics

Background:

• Current RNA-Seq analysis methods assume normalized counts represent true gene counts, ignoring potential variations.
• Existing statistical models for RNA sequencing (RNA-Seq) do not explicitly incorporate sequencing coverage information.
• This can lead to inaccuracies in identifying differentially expressed genes, especially with limited replicates or low coverage.

Purpose of the Study:

• To develop a novel statistical method for RNA-Seq analysis that incorporates sequencing coverage information.
• To improve the accuracy of estimating true gene counts from observed RNA-Seq data.
• To enhance differential gene expression analysis in challenging experimental conditions.

Main Methods:

• Developed a fast Bayesian approach to estimate the posterior distribution of true gene counts.
• Utilized sequencing coverage data, derived from RNA sample concentration, as a key parameter.
• Integrated this coverage parameter into a differential gene expression analysis pipeline.

Main Results:

• The new Bayesian method demonstrates superior or comparable performance to existing tools like NOISeq and GFOLD.
• Simulations and experiments with real unreplicated data validate the method's effectiveness.
• Successfully incorporated sequencing coverage into RNA-Seq differential expression analysis.

Conclusions:

• The developed method effectively addresses analytical bottlenecks in RNA-Seq experiments with few replicates and low sequencing coverage.
• The CORNAS (Coverage-dependent RNA-Seq) software implements this novel approach.
• This method offers a valuable tool for more accurate gene expression analysis in RNA-Seq studies.