BEERS2: RNA-Seq simulation through high fidelity in silico modeling

Area of Science:

• Bioinformatics
• Computational Biology
• Genomics

Background:

• Accurate simulation of RNA-sequencing (RNA-seq) reads is essential for developing and benchmarking bioinformatics tools.
• Existing RNA-seq simulators have seen limited advancement in the past decade, necessitating improved solutions.

Approach:

• Developed BEERS2, a novel RNA-seq simulator with a flexible, configurable design.
• BEERS2 simulates the entire RNA library preparation and sequencing pipeline in detail.
• Incorporates effects of polyA selection, ribosomal depletion, hexamer priming biases, PCR amplification biases (GC-content), barcode errors, and PCR errors.

Key Points:

• BEERS2 generates realistic RNA-seq reads in FASTQ, SAM, or BAM formats.
• Outputs true alignment information and transcript-level quantification values.
• Simulates various library preparation steps and technical biases for enhanced realism.

Conclusions:

• BEERS2 represents a comprehensive RNA-seq simulation tool to date.
• Demonstrates the utility of BEERS2 by analyzing the impact of simulation parameters on the Salmon pseudoalignment algorithm.