Accurate spliced alignment of long RNA sequencing reads

Sahlin K, Mäkinen V

Bioinformatics 37 (24) 4643-4651 [2021-12-11; online 2021-07-24]

Long-read RNA sequencing technologies are establishing themselves as the primary techniques to detect novel isoforms, and many such analyses are dependent on read alignments. However, the error rate and sequencing length of the reads create new challenges for accurately aligning them, particularly around small exons. We present an alignment method uLTRA for long RNA sequencing reads based on a novel two-pass collinear chaining algorithm. We show that uLTRA produces higher accuracy over state-of-the-art aligners with substantially higher accuracy for small exons on simulated and synthetic data. On simulated data, uLTRA achieves an accuracy of about 60% for exons of length 10 nucleotides or smaller and close to 90% accuracy for exons of length between 11 to 20 nucleotides. On biological data where true read location is unknown, we show several examples where uLTRA aligns to known and novel isoforms containing small exons that are not detected with other aligners. While uLTRA obtains its accuracy using annotations, it can also be used as a wrapper around minimap2 to align reads outside annotated regions. uLTRA is available at https://github.com/ksahlin/ultra. Supplementary data are available at Bioinformatics online.

Kristoffer Sahlin

SciLifeLab Fellow

PubMed 34302453

DOI 10.1093/bioinformatics/btab540

Crossref 10.1093/bioinformatics/btab540

pii: 6327681
pmc: PMC8665758


Publications 9.5.0