Installation and Configuration of FusionSeq
From GersteinInfo
User documentation main
Contents |
Installation and configuration
In order to install FusionSeq the external packages need to be installed first. Please, follow the instruction provided by the single packages. After they are installed, the first step for FusionSeq is the installation and configuration of BIOS. BIOS is a C library of useful general definitions for manipulating strings, arrays, and parser and more related to bioinformatic analysis. It requires the GSL library, which, in most systems, can be installed with the following commands (for details, please refer to the specific instructions at the GNU Scientific Library website):
$ cd path2gslSource/ $ ./configure --prefix=path2installation $ make $ make install
If a 64bit system is used, add CFLAGS=-m64 in the ./configure command. Similarly, the GD libraries can be installed in most systems with:
$ cd path2gdSource $ ./configure --prefix=path2installation --with-jpeg=/path2jpegLib/ $ make $ make install
although, they are NOT required for the core analysis. If you want to use them, please make sure that png, jpeg, zlib, freetype 2.x, and xpm are properly installed and linked. These are required by gfr2images in order to create a figure depicting wich regions of two genes are connected by PE-reads. See Installation and configuration of FusionSeq for setting the appropriate environmental variables.
Installing and configuring BIOS
To install BIOS a few variables need to be set before compiling the library. Here is an example of the procedure on a bash shell:
$ export BIOINFOCONFDIR=/pathToBios/bios/conf $ export BIOINFOGSLDIR=/pathTo/gsl-1.12 $ cd /pathToBios/bios $ make $ make prod
Please refer to BIOS documentation for additional information.
Installing and configuring ROOT
To install ROOT, please follow the instructions on the website. You may also include GSL library when compiling, but it is not a requirement. Once ROOT is installed, a few variables need to be defined in order to properly use this library with FusionSeq.
$ export ROOTSYS=/path2ROOT/ $ export PATH=$ROOTSYS/bin:$PATH
Installing and configuring FusionSeq
FusionSeq is composed of several programs divided into a set of core modules (to identify the candidate fusion transcripts), and a set of additional modules (to create images, BED, GFF and other auxiliary files for visualization and analysis) and CGIs for visualization of the results. To run the analysis, only the core modules are required. For the auxiliary modules, one needs to specify the location of the Drawing tool libraries by editing the specific section in the Makefile. For installing the visualization tools, please read Installing CGIs.
To install FusionSeq, the file geneFusionsConfig.h needs to be edited to specify to locations of the annotation files and other required tools:
geneFusionConfig.h: // --------------------------------- This section is required --------------------------------- // Location of the bowtie indexes #define BOWTIE_INDEXES "/path2bowtieIndexes" // Location of the bowtie indexes // #define BLAT_TWO_BIT_TO_FA "path2blat/twoBit2fa" // Location and filename for the reference genome in 2bit format (to be used by blat) #define BLAT_DATA_DIR "path2blat/DATA/blat" #define BLAT_TWO_BIT_DATA_FILENAME "hg18.2bit" // Location of the hg18 genome and name of the composite transcript model sequence file. Note that the file should be in the same directory #define GENOME_DIR "path2genome/" #define KNOWN_GENE_TRANSCRIPT_COMPOSITE_MODEL_FA_FILENAME "knownGeneAnnotationTranscriptCompositeModel.fa" // location of the annotation files #define ANNOTATION_DIR "/path2annotationFiles" // knownGeneAnnotation filename #define KNOWN_GENE_ANNOTATION_FILENAME "knownGeneAnnotation.txt" // conversion of knownGenes to gene symbols, description, etc. #define KNOWN_GENE_XREF_FILENAME "kgXref.txt" // conversion of knownGenes to TreeFam #define KNOWN_GENE_TREE_FAM_FILENAME "knownToTreefam.txt" // file name of the composite model file #define KNOWN_GENE_TRANSCRIPT_COMPOSITE_MODEL_FILENAME "knownGeneAnnotationTranscriptCompositeModel.txt" // ----------------------- This section is optional: visualization tools ------------------------- // URL of the cgi directory on the web server #define WEB_URL_CGI "http://cgiURL" // location of the data directory on the web server, as seen from the web server #define WEB_DATA_DIR "/path2data" // URL of the data directory on the web server #define WEB_DATA_LINK "http://dataURL" // Number of nucleotides flanking the region (for UCSC Genome Browser) #define UCSC_GENOME_BROWSER_FLANKING_REGION 500 // URL of the public website (non cgi) #define WEB_PUB_DIR "http://publicURL" // Location of the structural data for Circos #define WEB_SDATA_DIR "/path2structuralDataCircos" // Location of Circos installation #define WEB_CIRCOS_DIR "/path2circos"
Once the configuration file is ready, the core modules can be compiled and installed. For compiling the auxilliary modules, the Makefile needs to be properly edited (see Auxiliary modules). For the visualization tools, i.e. the CGIs, some additional variables need to be defined (see Installing CGIs). Once the configuration file is set up, the compilation just requires:
$ make // for the core analysis elements $ make all // for the core analysis elements as well as the additional programs and the visualization tools $ make cgi // for compiling the visualization/summary tools (see Installing CGIs) $ make deploy// for installing the visualization/summary tools to the web server
Auxiliary modules
These modules generate a set of useful data files for interpreting and visualizing the results. For example, gfr2gff generates the GFF files that can be displayed with the UCSC Genome Browser to show the location and the connection between the paired reads; gfr2fasta generates two fasta files containing the sequences of the reads, one for each end. Most of these modules do not require additional configuration, but gfr2images. This tool creates a schematic for each candidate showing what are the exons connected by paired-end reads. It uses graphic libraries whose locations need to be specified in the Makefile (section "optional parameters"). Here is an example on how to edit the Makefile.
GDDIR = -L/path2gd/gd-2.0.35/ GDINC = -I$(GDDIR)/include GDLIB = -L$(GDDIR)/lib PNGLIB = -L/usr/lib64 JPEGLIB = -L/usr/X11/lib ZLIB = -L/usr/lib FREETYPELIB = -L/usr/lib64
Note that GDINC and GDLIB are automatically defined once GDDIR is set. Once the Makefile is properly defined, the installation goes on as usual. However, to fully exploit the auxiliary modules, the CGIs should also be installed.
Installing CGIs
The visualization tools are rather useful to display the results of the analysis, although they are completely independent from the analysis itself. These tools require a web server able to interpret CGI programs. We tested our tools on an Apache Web server. First, a set of variables need to be specfied, describing the locations of the different tools:
$ export FUSIONSEQWEBSERVER=web_server_name $ export FUSIONSEQWEBUSER=webuserID $ export FUSIONSEQCGIDIR=/path2cgiDir $ export FUSIONSEQCGITARGET="-B target_architecture" // optional: to be specified only if using the CGIs on a machine with a different architecture than the core programs. On the webserver with the CGIs, excecute gcc -dumpmachine to get the target_architecture
The CGI programs assumes a certain directory structure in your data directory (WEB_DATA_DIR):
./ALIGNMENTS ./BED ./FASTA ./GFF ./IMAGES ./WIGS
BED, FASTA, GFF, and IMAGES contain the data generated by gfr2bed, gfr2fasta, gfr2gff and gfr2images, respectively. ALIGNMENTS and WIGS contain the results of the junction-sequence identification analysis, namely bp2alignment and bp2wigs. The user is required to ensure that these directories have the expected files.
One of the CGI applications is SeqViz, which is used to visualize Paired-End RNA-Seq reads. SeqViz requires the software package Circos in order to perform visualization. Please visit the Circos website to download the latest version of Circos and for detailed information on installing Circos and the required Perl modules.
There is a set of CSS style sheets, JavaScript scripts, and annotation files required for the CGIs. They may be downloaded here. In the Web Files tarball, the following files are included: The /web folder contains the CSS style sheets, images, and a distribution of the JQuery and JQuery UI Javascript libraries that are used by the CGI applications. Copy the contents of this folder to a non-CGI directory such as public_html. Set WEB_PUB_DIR to this directory. The /IMAGES folder contains images required by showDetails_cgi. Copy this folder to the directory specified by WEB_DATA_LINK. The /structdata folder contains genomic structure and annotation data files for Circos. Copy the contents of this folder to a directory for which Circos has sufficient permissions. Set WEB_SDATA_DIR to this directory.
