Ularcirc: Visualisation and enhanced analysis of circular RNAs via back and canonical forward splicing

Circular RNAs (circRNAs) are a unique class of transcripts that can only be identified from sequence alignments spanning discordant junctions, commonly referred to as backsplice junctions (BSJ). The challenges of detecting a BSJ from short read high throughput sequencing (HTS) data has steered software development to focus primarily on algorithmic methods to accurately capture BSJs. Here we present Ularcirc, the first software tool that provides a complete circRNA workflow from detection, integrated visualization, quality filtering of BSJ and forward splicing junctions (FSJ), through to sequence retrieval and downstream functional analysis. More importantly, Ularcirc uses an innovative method to filter out false positive circRNAs coined read alignment distribution (RAD) score which allows detection of circRNAs independent of gene annotations. We used Ularcirc to characterise circRNAs from public and in-house generated data sets and demonstrate how to discover (i) novel splicing patterns of parental transcripts, (ii) internal splicing patterns of circRNA, and (iii) the complexity of BSJ formation. Furthermore, we identify circRNAs that have potential open reading frames longer than their linear sequence. Finally, we have identified and validated the presence of a novel class of circRNA generated from ApoA4 transcripts whose BSJ derive from multiple sites within coding exons. Ularcirc can be accessed via https://github.com/VCCRI/Ularcirc.