Boquila: NGS Read Simulator to Eliminate Read Nucleotide Bias in Sequence Analysis

Overview

Journal Turk J Biol

Specialty Biology

Date 2023 Aug 2

PMID 37529166

Authors

Umit Akkose

Ogun Adebali

Affiliations

Soon will be listed here.

Abstract

Sequence content is heterogeneous throughout genomes. Therefore, genome-wide next-generation sequencing (NGS) reads biased towards specific nucleotide profiles are affected by the genome-wide heterogeneous nucleotide distribution. Boquila generates sequences that mimic the nucleotide profile of true reads, which can be used to correct the nucleotide-based bias of genome-wide distribution of NGS reads. Boquila can be configured to generate reads from only specified regions of the reference genome. It also allows the use of input DNA sequencing to correct the bias due to the copy number variations in the genome. Boquila uses standard file formats for input and output data, and it can be easily integrated into any workflow for high-throughput sequencing applications.

Citing Articles

UV-induced reorganization of 3D genome mediates DNA damage response.

Kaya V, Adebali O Nat Commun. 2025; 16(1):1376.

PMID: 39910043 PMC: 11799157. DOI: 10.1038/s41467-024-55724-7.

Global repair is the primary nucleotide excision repair subpathway for the removal of pyrimidine-pyrimidone (6-4) damage from the Arabidopsis genome.

Kaya S, Erdogan D, Sancar A, Adebali O, Oztas O Sci Rep. 2024; 14(1):3308.

PMID: 38332020 PMC: 10853524. DOI: 10.1038/s41598-024-53472-8.

The interplay of 3D genome organization with UV-induced DNA damage and repair.

Akkose U, Adebali O J Biol Chem. 2023; 299(5):104679.

PMID: 37028766 PMC: 10192929. DOI: 10.1016/j.jbc.2023.104679.

References

Xia Y, Liu Y, Deng M, Xi R . Pysim-sv: a package for simulating structural variation data with GC-biases. BMC Bioinformatics. 2017; 18(Suppl 3):53. PMC: 5374556. DOI: 10.1186/s12859-017-1464-8. View

Adebali O, Chiou Y, Hu J, Sancar A, Selby C . Genome-wide transcription-coupled repair in is mediated by the Mfd translocase. Proc Natl Acad Sci U S A. 2017; 114(11):E2116-E2125. PMC: 5358382. DOI: 10.1073/pnas.1700230114. View

Hu J, Adebali O, Adar S, Sancar A . Dynamic maps of UV damage formation and repair for the human genome. Proc Natl Acad Sci U S A. 2017; 114(26):6758-6763. PMC: 5495279. DOI: 10.1073/pnas.1706522114. View

Polz M, Cavanaugh C . Bias in template-to-product ratios in multitemplate PCR. Appl Environ Microbiol. 1998; 64(10):3724-30. PMC: 106531. DOI: 10.1128/AEM.64.10.3724-3730.1998. View

Yuan X, Zhang J, Yang L . IntSIM: An Integrated Simulator of Next-Generation Sequencing Data. IEEE Trans Biomed Eng. 2016; 64(2):441-451. DOI: 10.1109/TBME.2016.2560939. View

Yu Z, Du F, Ban R, Zhang Y . SimuSCoP: reliably simulate Illumina sequencing data based on position and context dependent profiles. BMC Bioinformatics. 2020; 21(1):331. PMC: 7379788. DOI: 10.1186/s12859-020-03665-5. View

Yang C, Chu J, Warren R, Birol I . NanoSim: nanopore sequence read simulator based on statistical characterization. Gigascience. 2017; 6(4):1-6. PMC: 5530317. DOI: 10.1093/gigascience/gix010. View

Hu X, Yuan J, Shi Y, Lu J, Liu B, Li Z . pIRS: Profile-based Illumina pair-end reads simulator. Bioinformatics. 2012; 28(11):1533-5. DOI: 10.1093/bioinformatics/bts187. View

Shcherbina A . FASTQSim: platform-independent data characterization and in silico read generation for NGS datasets. BMC Res Notes. 2014; 7:533. PMC: 4246604. DOI: 10.1186/1756-0500-7-533. View

10.

Ross M, Russ C, Costello M, Hollinger A, Lennon N, Hegarty R . Characterizing and measuring bias in sequence data. Genome Biol. 2013; 14(5):R51. PMC: 4053816. DOI: 10.1186/gb-2013-14-5-r51. View

11.

Ivakhno S, Colombo C, Tanner S, Tedder P, Berri S, Cox A . tHapMix: simulating tumour samples through haplotype mixtures. Bioinformatics. 2016; 33(2):280-282. DOI: 10.1093/bioinformatics/btw589. View

12.

Wang J, Quach B, Furey T . Correcting nucleotide-specific biases in high-throughput sequencing data. BMC Bioinformatics. 2017; 18(1):357. PMC: 5540620. DOI: 10.1186/s12859-017-1766-x. View

13.

Pattnaik S, Gupta S, Rao A, Panda B . SInC: an accurate and fast error-model based simulator for SNPs, Indels and CNVs coupled with a read generator for short-read sequence data. BMC Bioinformatics. 2014; 15:40. PMC: 3926339. DOI: 10.1186/1471-2105-15-40. View

14.

Hu J, Lieb J, Sancar A, Adar S . Cisplatin DNA damage and repair maps of the human genome at single-nucleotide resolution. Proc Natl Acad Sci U S A. 2016; 113(41):11507-11512. PMC: 5068337. DOI: 10.1073/pnas.1614430113. View

15.

Grant G, Farkas M, Pizarro A, Lahens N, Schug J, Brunk B . Comparative analysis of RNA-Seq alignment algorithms and the RNA-Seq unified mapper (RUM). Bioinformatics. 2011; 27(18):2518-28. PMC: 3167048. DOI: 10.1093/bioinformatics/btr427. View

16.

Pratas D, Pinho A, Rodrigues J . XS: a FASTQ read simulator. BMC Res Notes. 2014; 7:40. PMC: 3927261. DOI: 10.1186/1756-0500-7-40. View

17.

Mu J, Mohiyuddin M, Li J, Asadi N, Gerstein M, Abyzov A . VarSim: a high-fidelity simulation and validation framework for high-throughput genome sequencing with cancer applications. Bioinformatics. 2014; 31(9):1469-71. PMC: 4410653. DOI: 10.1093/bioinformatics/btu828. View

18.

Caboche S, Audebert C, Lemoine Y, Hot D . Comparison of mapping algorithms used in high-throughput sequencing: application to Ion Torrent data. BMC Genomics. 2014; 15:264. PMC: 4051166. DOI: 10.1186/1471-2164-15-264. View

19.

Johnson S, Trost B, Long J, Pittet V, Kusalik A . A better sequence-read simulator program for metagenomics. BMC Bioinformatics. 2014; 15 Suppl 9:S14. PMC: 4168713. DOI: 10.1186/1471-2105-15-S9-S14. View

20.

Jia B, Xuan L, Cai K, Hu Z, Ma L, Wei C . NeSSM: a Next-generation Sequencing Simulator for Metagenomics. PLoS One. 2013; 8(10):e75448. PMC: 3790878. DOI: 10.1371/journal.pone.0075448. View