Welcome to the blog

Posts

My thoughts and ideas

De Novo Transcript Assembly | Griffith Lab

RNA-seq Bioinformatics

Introduction to bioinformatics for RNA sequence analysis

De Novo Transcript Assembly


RNA-seq_Flowchart5


Running Stringtie in De Novo Mode

Note, to discover novel transcripts with Stringtie using the alignments we generated in the previous modules we will now run Stringtie in de novo mode. To use de novo mode do NOT specify either of the -G OR -e options.

Extra options specified below

  • –rf tells StringTie that our data is stranded and to use the correct strand specific mode (i.e. assume a stranded library fr-firststrand).
  • -p 4 tells Stringtie to use eight CPUs
  • -l name prefix for output transcripts (default: STRG)
  • -o output path/file name for the assembled transcripts GTF (default: stdout)
cd $RNA_HOME/
mkdir -p expression/stringtie/de_novo/
cd expression/stringtie/de_novo/

stringtie --rf -p 4 -l HBR_Rep1 -o HBR_Rep1/transcripts.gtf $RNA_ALIGN_DIR/HBR_Rep1.bam
stringtie --rf -p 4 -l HBR_Rep2 -o HBR_Rep2/transcripts.gtf $RNA_ALIGN_DIR/HBR_Rep2.bam
stringtie --rf -p 4 -l HBR_Rep3 -o HBR_Rep3/transcripts.gtf $RNA_ALIGN_DIR/HBR_Rep3.bam

stringtie --rf -p 4 -l UHR_Rep1 -o UHR_Rep1/transcripts.gtf $RNA_ALIGN_DIR/UHR_Rep1.bam
stringtie --rf -p 4 -l UHR_Rep2 -o UHR_Rep2/transcripts.gtf $RNA_ALIGN_DIR/UHR_Rep2.bam
stringtie --rf -p 4 -l UHR_Rep3 -o UHR_Rep3/transcripts.gtf $RNA_ALIGN_DIR/UHR_Rep3.bam

Reference Guided Transcript Assembly | Griffith Lab

RNA-seq Bioinformatics

Introduction to bioinformatics for RNA sequence analysis

Reference Guided Transcript Assembly


RNA-seq_Flowchart4


De novo transcript discovery and differential expression using Stringtie and Ballgown

In the previous module, we ran Stringtie in reference only mode using the -G and -e Stringtie options.

In this module we will run Stringtie in two additional modes: (1) reference guided mode and (2) de novo mode. Stringtie can predict the transcripts present in each library with or without help from knowledge of known transcripts. Stringtie will then assign arbitrary transcript IDs to each transcript assembled from the data and estimate expression for those transcripts. One complication with this method is that in each library, a different set of transcripts is likely to be predicted. There may be a lot of similarities but the number of transcripts and their exact structure will differ in the output files for each library. Before you can compare across libraries you therefore need to determine which transcripts correspond to each other across the libraries. Stringtie provides a merge command to combine predicted transcript GTF files from across different libraries.

Once you have a merged GTF file you can run Stringtie with this instead of the known transcripts GTF file we used before. The merged GTF is used to recalculate expression estimates in prepartion for running Ballgown using the merged, novel transcripts.

To run Stringtie in ‘reference guided’ mode: use the -G option WITHOUT -e.

To run Stringtie in ‘de novo’ mode do NOT specify either of the -G OR -e options.

Refer to the Stringtie manual for a more detailed explanation: https://ccb.jhu.edu/software/stringtie/index.shtml?t=manual

Running Stringtie in Reference Guided Mode

Using the alignments we generated in the previous modules we will now run Stringtie in reference guided mode using the -G option ONLY.

Extra options specified below:

First, create an output directory and then run stringtie in reference-guided mode.

cd $RNA_HOME/
mkdir -p expression/stringtie/ref_guided/
cd expression/stringtie/ref_guided/

stringtie --rf -p 4 -G $RNA_REF_GTF -l HBR_Rep1 -o HBR_Rep1/transcripts.gtf $RNA_ALIGN_DIR/HBR_Rep1.bam
stringtie --rf -p 4 -G $RNA_REF_GTF -l HBR_Rep2 -o HBR_Rep2/transcripts.gtf $RNA_ALIGN_DIR/HBR_Rep2.bam
stringtie --rf -p 4 -G $RNA_REF_GTF -l HBR_Rep3 -o HBR_Rep3/transcripts.gtf $RNA_ALIGN_DIR/HBR_Rep3.bam

stringtie --rf -p 4 -G $RNA_REF_GTF -l UHR_Rep1 -o UHR_Rep1/transcripts.gtf $RNA_ALIGN_DIR/UHR_Rep1.bam
stringtie --rf -p 4 -G $RNA_REF_GTF -l UHR_Rep2 -o UHR_Rep2/transcripts.gtf $RNA_ALIGN_DIR/UHR_Rep2.bam
stringtie --rf -p 4 -G $RNA_REF_GTF -l UHR_Rep3 -o UHR_Rep3/transcripts.gtf $RNA_ALIGN_DIR/UHR_Rep3.bam