MntJULiP and Jutils: Differential splicing analysis of RNA-seq data with covariates

Differences in alternative splicing patterns can reveal important markers of phenotypic differentiation, including biomarkers of disease. Emerging large and complex RNA-seq datasets from disease and population studies include multiple confounders such as sex, age, ethnicity and clinical attributes, which demand highly specialized data analysis tools. However, few methods are equipped to handle the new challenges. We describe an implementation of our programs MntJULiP and Jutils for differential splicing detection and visualization from RNA-seq data that takes into account covariates. MntJULiP detects intron-level differences in alternative splicing from RNA-seq data using a Bayesian mixture model. Jutils visualizes alternative splicing variation with heatmaps, PCA and sashimi plots, and Venn diagrams. Our tools are scalable and can process thousands of samples within hours. We applied our methods to the collection of GTEx brain RNA-seq samples to deconvolute the effects of sex and age at death on the splicing patterns. In particular, clustering of covariate adjusted data identifies a subgroup of individuals undergoing a distinct splicing program during aging. MntJULiP and Jutils are implemented in Python and are available from https://github.com/splicebox/.