  Although RNA-Seq data provide unprecedented isoform-level expression
information, detection of alternative isoform regulation (AIR) remains
difficult, particularly when working with an incomplete transcript annotation.
We introduce JunctionSeq, a new method that builds on the statistical
techniques used by the well-established DEXSeq package to detect differential
usage of both exonic regions and splice junctions. In particular, JunctionSeq
is capable of detecting differentials in novel splice junctions without the
need for an additional isoform assembly step, greatly improving performance
when the available transcript annotation is flawed or incomplete. JunctionSeq
also provides a powerful and streamlined visualization toolset that allows
bioinformaticians to quickly and intuitively interpret their results. We tested
our method on publicly available data from several experiments performed on the
rat pineal gland and Toxoplasma gondii, successfully detecting known and
previously validated AIR genes in 19 out of 19 gene-level hypothesis tests. Due
to its ability to query novel splice sites, JunctionSeq is still able to detect
these differentials even when all alternative isoforms for these genes were not
included in the transcript annotation. JunctionSeq thus provides a powerful
method for detecting alternative isoform regulation even with low-quality
annotations. An implementation of JunctionSeq is available as an R/Bioconductor
package.
