The purpose of this exercise is to introduce tools for analyzing differential gene expression in RNA-seq data. You will analyze RNA-seq data from human reference (a mix of tissues) and brain tissue to identify genes for which expression is enriched in the brain. Due to time constraints, the data we will analyze is a subset (chr22) of a larger RNA-seq dataset.
We will use a suite of tools called the Tuxedo pipeline. For an additional tutorial, see the following paper:
Trapnell et al. 2014, Nature Protocols. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks.
Note: this pipeline is no longer updated and has been replaced with a more efficient and accurate pipeline, however, we will use the original Tuxedo pipeline because there is far more support currently available for it and it has fewer bugs to contend with. The new and improved pipeline consists of a similar suite of tools: HISAT, StringTie, and Ballgown. See https://www.nature.com/articles/nprot.2016.095
1. Assess data quality using FastQC.
2. Quality filter datasets using Trimmomatic.
3. Align the RNA-seq reads to the human genome using TopHat2.
4. Assemble transcripts based on RNA-seq data using cufflinks and cuffmerge.
5. Compare expression differences using cuffdiff.
6. Visualize data using the genome viewing software IGV.
7. Plot data with R and cummeRbund.
1. Open a terminal window and change into the
There are 12 RNA-seq datasets corresponding to paired-end data for 3 replicates from two sample sets (brain and ref). Examine a few lines of one of the files using
zless. What information is contained in each line?
2. Assess the quality of the data using
See FastQC tutorial for additional details: https://www.youtube.com/watch?v=bz93ReOv87Y
3. Trim adapter sequences and quality filter the RNA-seq data (fastq files) using
Trim adapter sequences and quality filter each dataset using Trimmomatic (you will run trimmomatic 6 times in total).
$ trimmomatic PE 'input_fastq_1' 'input_fastq_2' 'trimmed_P1' 'trimmed_P2' 'trimmed_U1' 'trimmed_U2' ILLUMINACLIP:/Users/bz360/Documents/TruSeq3-PE.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36
For output file names, use:
See the Trimmomatic manual for a detailed description of options: http://www.usadellab.org/cms/uploads/supplementary/Trimmomatic/TrimmomaticManual_V0.32.pdf
Submit an answer to the following question on Canvas:
What proportion of the reads in each library was retained?
4. Assess the quality of one of the datasets after quality filtering using
5. Create a
bowtie index for the human chromosome 22 sequence:
$ bowtie2-build 'sequence.fa' 'prefix'
The chr22 sequence should in the
You will need to decompress the file if it has the
If you don't have the chromosome sequence, you can download it here:
For the bowtie prefix, use