PyDESeq2: a python package for bulk RNA-seq differential expression analysis.

Area of Science:

• Bioinformatics
• Computational Biology
• Genomics

Background:

• Differential gene expression analysis is crucial for understanding biological processes and disease mechanisms.
• RNA-sequencing (RNA-seq) is a widely used technology for quantifying gene expression levels.
• Existing tools like DESeq2 are powerful but may lack seamless integration with Python-based data science ecosystems.

Purpose of the Study:

• To present PyDESeq2, a novel Python implementation of the DESeq2 workflow for differential expression analysis.
• To evaluate PyDESeq2's performance, accuracy, and compatibility with modern Python data science tools.
• To provide an open-source, accessible tool for researchers analyzing bulk RNA-seq data.

Main Methods:

• Re-implementation of the DESeq2 differential expression analysis workflow in Python.
• Performance evaluation using large-scale datasets, including The Cancer Genome Atlas (TCGA) data.
• Assessment of model likelihood, computational speed, and interoperability with Python libraries.

Main Results:

• PyDESeq2 produces results similar to, though not identical to, the original DESeq2.
• Achieves higher model likelihood compared to the original DESeq2.
• Demonstrates significant speed improvements on large datasets and offers easier integration with Python tools.

Conclusions:

• PyDESeq2 provides a viable and efficient Python-based alternative for differential gene expression analysis.
• The tool enhances computational speed and facilitates integration within the Python data science ecosystem.
• PyDESeq2 is released as open-source software, promoting accessibility and further development within the scverse ecosystem.