Differential Usage of Exons and Splice Junctions
> v0.6.16 (Revised Wed Nov 25 14:52:43 EST 2015)
JunctionSeq is an R package designed to detect and assess differential usage of exons and splice junctions in High-Throughput, Next-Generation RNA-Seq datasets. The methodology is based on the methods used by the DEXSeq bioconductor package, originally proposed by Anders, Reyes, and Huber.
Help, documentation, and the most recent release of JunctionSeq is available on the JunctionSeq github pages.
If you do not find an answer to your question here, you can email the developer at JunctionSeq-Contact (at) list.nih.gov.
Note: the current version of JunctionSeq is ONLY compatible with Bioconductor 3.2 or higher. For older versions of JunctionSeq compatible with Bioconductor 3.0 and 3.1, see JunctionSeq release v0.5.1.
For more information see the JunctionSeq vignette or the online reference documentation.
There is also a comprehensive walkthrough of the entire analysis pipeline, along with a full example dataset with example bam files.
For help with individual R functions in the R utility, use the R command:
> help(functionname);
For a full listing of all help topics for the R utility, use the R command:
> help(package="JunctionSeq");
JunctionSeq can be installed automatically from source using the R commands:
> source("http://hartleys.github.io/JunctionSeq/install/JS.install.R");
> JS.install();
In order to successfully install from source on windows, you must have Rtools installed. Alternatively you can install the pre-compiled binaries, see next section. Installation on OSX is more complex, and requires Xcode and gfortran 4.8.2 (see here). For installation on OSX, you may want to install the pre-compiled binaries (see below).
If you encounter problems with installation, you can install all prerequisite packages manually using the R commands:
#Install CRAN packages:
install.packages("statmod")
install.packages("plotrix")
install.packages("stringr")
install.packages("Rcpp")
install.packages("RcppArmadillo")
install.packages("locfit")
#Install Bioconductor packages:
source("http://bioconductor.org/biocLite.R");
biocLite();
biocLite("Biobase");
biocLite("BiocGenerics");
biocLite("BiocParallel");
biocLite("GenomicRanges");
biocLite("IRanges");
biocLite("S4Vectors");
biocLite("genefilter");
biocLite("geneplotter");
biocLite("SummarizedExperiment");
You can then install JunctionSeq itself with the command:
#Install JunctionSeq (from source):
install.packages("http://hartleys.github.io/JunctionSeq/install/JunctionSeq_LATEST.tar.gz",
repos = NULL,
type="source")
If you are installing to windows, you will also require Rtools which allows advanced packages to be installed from source-code. Installing from source to OSX is a little more complex, and requires Xcode and gfortran 4.8.2. See here for a description of some of these issues.
Alternatively, you can also install the pre-compiled binaries. This only works for windows and recent versions of OSX. It has been tested on OSX v14.5 (Darwin).
#Install Windows binary:
install.packages("http://hartleys.github.io/JunctionSeq/install/JunctionSeq_0.6.16.zip",
repos=NULL);
#OR
#Install OSX binary:
install.packages("http://hartleys.github.io/JunctionSeq/install/JunctionSeq_0.6.16.tgz",
repos=NULL);
JunctionSeq may use large amounts of RAM. You can reduce the memory usage considerably by reducing the number of cores used with the nCores parameter. Unfortunetely, BiocParallel duplicates the entire environment whenever it runs in multicore mode, so amount of RAM required is multiplied by the number of cores in use.
The exact memory requirements will vary depending on a large number of different factors, such as genome size/complexity, number of replicates, and the number of novel splice junctions.
This software package is licensed under the GNU-GPL v3:
This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.
This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
You should have received a copy of the GNU General Public License along with this program. If not, see http://www.gnu.org/licenses/.
Portions of this software are "United States Government Work"
under the terms of the United States Copyright Act.
It was written as part of the authors' official duties for the
United States Government and thus those portions cannot be
copyrighted. Those portions of this software are freely
available to the public for use without a copyright notice.
Restrictions cannot be placed on its present or future use.
Although all reasonable efforts have been taken to ensure the accuracy and reliability of the software and data, the National Human Genome Research Institute (NHGRI) and the U.S. Government does not and cannot warrant the performance or results that may be obtained by using this software or data. NHGRI and the U.S. Government disclaims all warranties as to performance, merchantability or fitness for any particular purpose.