Lerna: transformer architectures for configuring error correction tools for short- and long-read genome sequencing
Autor: | Atul Sharma, Pranjal Jain, Ashraf Mahgoub, Zihan Zhou, Kanak Mahadik, Somali Chaterji |
---|---|
Rok vydání: | 2021 |
Předmět: |
FOS: Computer and information sciences
Computer Science - Machine Learning QH301-705.5 Computer applications to medicine. Medical informatics R858-859.7 Biochemistry Machine Learning (cs.LG) Natural language processing (NLP) Nanopore reads Structural Biology Quantitative Biology - Genomics Perplexity metric Biology (General) Error correction Molecular Biology Genomics (q-bio.GN) Base Sequence Applied Mathematics Methodology Article High-Throughput Nucleotide Sequencing Genomics Sequence Analysis DNA Automated configuration tuning Computer Science Applications Transformer networks FOS: Biological sciences Parameter search space PacBio reads Algorithms Software |
Zdroj: | BMC Bioinformatics BMC Bioinformatics, Vol 23, Iss 1, Pp 1-26 (2022) |
ISSN: | 1471-2105 |
Popis: | Sequencing technologies are prone to errors, making error correction (EC) necessary for downstream applications. EC tools need to be manually configured for optimal performance. We find that the optimal parameters (e.g., k-mer size) are both tool- and dataset-dependent. Moreover, evaluating the performance (i.e., Alignment-rate or Gain) of a given tool usually relies on a reference genome, but quality reference genomes are not always available. We introduce Lerna for the automated configuration of k-mer-based EC tools. Lerna first creates a language model (LM) of the uncorrected genomic reads; then, calculates the perplexity metric to evaluate the corrected reads for different parameter choices. Next, it finds the one that produces the highest alignment rate without using a reference genome. The fundamental intuition of our approach is that the perplexity metric is inversely correlated with the quality of the assembly after error correction. Results: First, we show that the best k-mer value can vary for different datasets, even for the same EC tool. Second, we show the gains of our LM using its component attention-based transformers. We show the model's estimation of the perplexity metric before and after error correction. The lower the perplexity after correction, the better the k-mer size. We also show that the alignment rate and assembly quality computed for the corrected reads are strongly negatively correlated with the perplexity, enabling the automated selection of k-mer values for better error correction, and hence, improved assembly quality. Additionally, we show that our attention-based models have significant runtime improvement for the entire pipeline -- 18X faster than previous works, due to parallelizing the attention mechanism and the use of JIT compilation for GPU inferencing. 26 pages, 5 figures, 10 tables. Accepted to BMC Bioinformatics |
Databáze: | OpenAIRE |
Externí odkaz: |