Introduction to Bulk RNA-seq data analysis

15th, 22nd and 29th November 2021

Taught remotely

Bioinformatics Training, Craik-Marshall Building, Downing Site, University of Cambridge

Instructors

Outline

In this workshop, you will be learning how to analyse RNA-seq data. This will include read alignment, quality control, quantification against a reference, reading the count data into R, performing differential expression analysis, and gene set testing, with a focus on the DESeq2 analysis workflow. You will learn how to generate common plots for analysis and visualisation of gene expression data, such as boxplots and heatmaps.

This workshop is aimed at biologists interested in learning how to perform differential expression analysis of RNA-seq data.

Whilst we have run this course for several years, we are still learning how to teach it remotely and this is the first time split over three weeks.
Please bear with us if there are any technical hitches, and be aware that timings for different sections laid out in the schedule below may not be adhered to. There may be some necessity to make adjusments to the course as we go.

Prerequisites

**Some basic experience of using a UNIX/LINUX command line is assumed**

**Some R knowledge is assumed and essential. Without it, you will struggle on this course.** If you are not familiar with the R statistical programming language we strongly encourage you to work through an introductory R course before attempting these materials. We recommend our Introduction to R course

Shared Google Document

This Google Document contains useful information and links..

Please use it to post any questions you have during the course.

The trainers will be monitoring the document and will answer questions as quickly as they can.

Introduce Yourself

There is another Google Doc here. Please write a couple sentences here to introduce yourself to the class, tell us a bit about your background and what you hope to get out of this course. If you are a student or staff at the University of Cambridge, tell us which Department you are in.

Course etiquette

As this course is being taught online and there are a large number of participants, we will all need to follow a few simple rules to ensure things run as smoothly as possible:

  1. Please mute your microphone

  2. To get help from a tutor, please click the “Raise Hand” button in Zoom:

    This can be found by clicking on the ‘Participants’ button. A tutor will then contact you in the chat. If necessary, you and the tutor can be moved to a breakout room where you can discuss your issue in more detail.

  3. Please ask any general question by typing it into the Google Doc mentioned above

  4. During practicals, when you are done, please press the green “Yes” button:

    This way we will know when we can move on.

Timetable

We are still learning how to teach this course remotely, all times here should be regarded as aspirations

Day 1

9:30 - 9:45 - Welcome!

9:45 - 10:15 - Introduction to RNAseq Methods - Ashley D Sawle

10:15 - 11:15 Raw read file format and QC - Chandra Chilamakuri - Practical (pdf)
- Practical solutions (pdf)

11:15 - 12:45 Short read alignment with HISAT2 - Jon Price
- Practical (pdf)
- Practical solutions (pdf)

12:45 - 13:45 Lunch

13:45 - 15:30 QC of alignment - Jon Price
- Practical (pdf)
<- Practical solutions (pdf)

15:30 - 17:00 Quantification of Gene Expression with Salmon - Abbi Edwards
- Practical (pdf)

Day 2

9:30 - 10:15 Introduction to RNAseq Analysis in R - Jon Price

10:15 - 12:15 - RNA-seq Data Exploration (pdf) - Ashley Sawle
Practical solutions (pdf)
Live script

12:15 - 13:15 Lunch

13:15 - 15:45 Statistical Analysis of Bulk RNAseq Data

15:45 - 17:00 - Differential Expression for RNA-seq - Part 1 (pdf) - Abbi Edwards Live script

Day 3

9.30 - 10.00 - Recap of Day 1 and 2 - Hugo Tavares

10.00 - 12:00 - Differential Expression for RNA-seq - Part 2 (pdf) - Abbi Edwards
- practical solutions (pdf)
- live script - extra models plots

12.00 - 12.45 Annotation and Visualisation of RNA-seq results - Chandra Chilamakuri

12.45 - 13.45 Lunch

13.45 - 15.00 Annotation and Visualisation of RNA-seq results (pdf) - Chandra Chilamakuri

15:00 - 17:00 Gene-set testing - Stephane Ballereau

Source Materials for Practicals

The lecture slides and other source materials, including R code and practical solutions, can be found in the course’s Github repository

Extended materials

The Extended Materials contain extensions to some of the sessions and additional materials, including instruction on downloading and processing the raw data for this course, a link to an excellent R course, and where to get further help after the course.

Additional Resources

Acknowledgements

This course is based on the course RNAseq analysis in R prepared by Combine Australia and delivered on May 11/12th 2016 in Carlton. We are extremely grateful to the authors for making their materials available; Maria Doyle, Belinda Phipson, Matt Ritchie, Anna Trigos, Harriet Dashnow, Charity Law.

The materials have been rewritten/modified/corrected/updated by various contributors over the past 5 years including:

Abigail Edwards Ashley D Sawle Chandra Chilamakuri Dominique-Laurent Couturier Guillermo Parada González Hugo Tavares Jon Price Mark Dunning Mark Fernandes Oscar Rueda Sankari Nagarajan Stephane Ballereau Zeynep Kalender Atak

Apologies if we have missed anyone!