Skip to content

Differential Splicing

obigriffith edited this page Nov 18, 2014 · 27 revisions

RNA-seq Flowchart - Module 5

#5-iv. Differential Splicing

Use Cuffdiff to compare the UHR and HBR conditions.

Refer to the Cufflinks manual for a more detailed explanation: http://cufflinks.cbcb.umd.edu/manual.html#cuffdiff

Cuffdiff basic usage:

cuffdiff [options] <transcripts.gtf> <sample1_hits.sam> <sample2_hits.sam> [... sampleN_hits.sam]
  • Supply replicate SAMs as comma separated lists for each condition: sample1_rep1.sam,sample1_rep2.sam,...sample1_repM.sam
  • '-p 8' tells cuffdiff to use eight CPUs
  • '-L' tells cuffdiff the labels to use for samples

For reference guided mode:

cd $RNA_HOME/
mkdir -p de/tophat_cufflinks/ref_guided/
cd $RNA_HOME/alignments/tophat/

Perform UHR vs. HBR comparison, for known/novel (reference guided mode) transcripts:

cuffdiff -p 8 -L UHR,HBR -o $RNA_HOME/de/tophat_cufflinks/ref_guided/ --frag-len-mean 262 --frag-len-std-dev 80 --no-update-check $RNA_HOME/expression/tophat_cufflinks/ref_guided/merged/merged.gtf UHR_Rep1_ERCC-Mix1/accepted_hits.bam,UHR_Rep2_ERCC-Mix1/accepted_hits.bam,UHR_Rep3_ERCC-Mix1/accepted_hits.bam HBR_Rep1_ERCC-Mix2/accepted_hits.bam,HBR_Rep2_ERCC-Mix2/accepted_hits.bam,HBR_Rep3_ERCC-Mix2/accepted_hits.bam

For de novo mode: cd $RNA_HOME/ mkdir de/tophat_cufflinks/de_novo/ cd $RNA_HOME/alignments/tophat/

Perform Tumor vs. Normal comparison, using just the capture libraries, for known/novel (de novo mode) transcripts cuffdiff -p 8 -L UHR,HBR -o $RNA_HOME/de/tophat_cufflinks/de_novo/ --frag-len-mean 262 --frag-len-std-dev 80 --no-update-check $RNA_HOME/expression/tophat_cufflinks/de_novo/merged/merged.gtf UHR_Rep1_ERCC-Mix1/accepted_hits.bam,UHR_Rep2_ERCC-Mix1/accepted_hits.bam,UHR_Rep3_ERCC-Mix1/accepted_hits.bam HBR_Rep1_ERCC-Mix2/accepted_hits.bam,HBR_Rep2_ERCC-Mix2/accepted_hits.bam,HBR_Rep3_ERCC-Mix2/accepted_hits.bam

| Previous Section | This Section | Next Section | |:----------------------------------------------:|:-----------------------------------------------:|:-----------------------------------------------------------------:| | Merging | Differential Splicing | Splicing Visualization |