Improving spliced alignment by modeling splice sites with deep learning
Journal:
arXiv
Published Date:
Jun 15, 2025
Abstract
Motivation: Spliced alignment refers to the alignment of messenger RNA (mRNA)
or protein sequences to eukaryotic genomes. It plays a critical role in gene
annotation and the study of gene functions. Accurate spliced alignment demands
sophisticated modeling of splice sites, but current aligners use simple models,
which may affect their accuracy given dissimilar sequences.
Results: We implemented minisplice to learn splice signals with a
one-dimensional convolutional neural network (1D-CNN) and trained a model with
7,026 parameters for vertebrate and insect genomes. It captures conserved
splice signals across phyla and reveals GC-rich introns specific to mammals and
birds. We used this model to estimate the empirical splicing probability for
every GT and AG in genomes, and modified minimap2 and miniprot to leverage
pre-computed splicing probability during alignment. Evaluation on human
long-read RNA-seq data and cross-species protein datasets showed our method
greatly improves the junction accuracy especially for noisy long RNA-seq reads
and proteins of distant homology.
Availability and implementation: https://github.com/lh3/minisplice