Upgrade Trinity + RSEM

* Upgrade Trinity

cd /root/
wget  http://garr.dl.sourceforge.net/project/trinityrnaseq/trinityrnaseq_r20131110.tar.gz 
tar xzf trinityrnaseq_r20131110.tar.gz 
cd trinityrnaseq_r20131110/ 
make

cd ..
mv /usr/share/trinity /usr/share/trinityOLD
mv trinityrnaseq_r20131110 /usr/share/trinity 
chown -R root:root /usr/share/trinity
chmod -R 755 /usr/share/trinity

* Install RSEM

cd /root/
wget  http://deweylab.biostat.wisc.edu/rsem/src/rsem-1.2.7.tar.gz  
tar -xvzf rsem-1.2.7.tar.gz 
cd rsem-1.2.7/  
make

cd ..
mv rsem-1.2.7 /usr/local/bin/  
chown -R root:root /usr/local/bin/rsem-1.2.7 
cd /usr/local/bin/
ln -s rsem-1.2.7 rsem 

* User script to activate Trinity and RSEM

export PATH=/usr/share/trinity:/usr/share/trinity/util:/usr/share/trinity/util/RSEM_util:/usr/local/bin/rsem:$PATH
export TRINITY_HOME=/usr/share/trinity
ulimit -s unlimited 

* Test installation

Logon and activate Trinity

user@rachaelx:~$ source trinityactivate.sh

Copy sample directory

user@rachaelx:~$ cp -R /usr/share/trinity/sample_data ~/

Edit test script and remove script path

user@rachaelx:~$ cd sample_data/test_Trinity_Assembly/
user@rachaelx:~/sample_data/test_Trinity_Assembly$ nano -w runMe.sh
#!/bin/bash -ve
 
if [ -e reads.left.fq.gz ] && ! [ -e reads.left.fq ]
then
    gunzip -c reads.left.fq.gz > reads.left.fq
fi
if [ -e reads.right.fq.gz ] && ! [ -e reads.right.fq ]
then
    gunzip -c reads.right.fq.gz > reads.right.fq
fi
 
#######################################################
##  Run Trinity to Generate Transcriptome Assemblies ##
#######################################################
 
## use jellyfish
Trinity.pl --seqType fq --JM 2G --left reads.left.fq.gz --right reads.right.fq.gz --SS_lib_type RF --CPU 4 --no_cleanup
##### Done Running Trinity #####
 
if [ ! $* ]; then
    exit 0
fi
sleep 2
 
######################################################
## align reads back to the transcripts using Bowtie ##
######################################################
 
sleep 2
alignReads.pl --left reads.left.fq --right reads.right.fq --target trinity_out_dir/Trinity.fasta --aligner bowtie --seqType fq --SS_lib_type RF
##### Done aligning reads #######
sleep 2
 
###########################################
# use RSEM to estimate read abundance  ####
###########################################
 
sleep 2
run_RSEM.pl --transcripts trinity_out_dir/Trinity.fasta --name_sorted_bam bowtie_out/bowtie_out.nameSorted.sam.+.sam.PropMapPairsForRSEM.bam --paired
###### Done running RSEM ########

Exec test

user@rachaelx:~/sample_data/test_Trinity_Assembly$ ./runMe.sh

Result

...
##### Done Running Trinity #####

Clean test dirs

user@rachaelx:~/sample_data/test_Trinity_Assembly$ ./cleanme.pl
 
 
rachaelx/uptrinityrsem.txt ยท Last modified: 2013/11/13 18:22 by giancarlo

Developers: CNR Ceris IT Office and Library
Giancarlo Birello (g.birello _@_ ceris.cnr.it) and Anna Perin (a.perin _@_ ceris.cnr.it)
BioInfo@TO.CNR is licensed under: Creative Commons License
Recent changes RSS feed Creative Commons License Valid XHTML 1.0 Valid CSS Driven by DokuWiki
Drupal Garland Theme for Dokuwiki