genomics for systems biology / ISB2020 RNA sequencing (RNA-seq)

Transcription

1 RNA sequencing (RNA-seq)

2 Module Outline MO 13-Mar-2017 RNA sequencing: Introduction 1 WE 15-Mar-2017 RNA sequencing: Introduction 2 MO 20-Mar-2017 Paper: PMID : Human genomics. The human transcriptome across tissues and individuals. WE 22-Mar-2017 RNA sequencing: Models for RNA-sequencing Data MO 27-Mar-2017 Paper: PMID : A survey of best practices for RNA-seq data analysis WE 29-Mar-2017 Paper: PMID : Understanding development and stem cells using single cell-based analyses of gene expression

3 High Throughput Technologies Measuring things in parallel, for example: Northern blot: Expression of one gene. Microarray: Expression of all known genes. NGS (e.g. RNASeq): The transcriptome. NGS: Next Generation Sequencing: Generates thousands of sequences in parallel Underlies many applications, not only RNASeq Dramatically reduced cost

4 Next Generation Sequencing - Sequencing target (e.g. a genome) - Library preparation (includes fragmentation) - Single vs. paired end - Different technologies / platforms Differ in: - Read length - Accuracy - Cost Each fragment gets sequenced in parallel Much more efficient than sequential sequencing

5 Next generation sequencing [NHGRI: ]

6 What can we measure with NGS?

7 What can we measure with NGS? [Soon et al. 2013]

8 HighThroughput Technologies High-throughput technologies measure things in parallel Many based on Next Generation Sequencing We will focus on measuring transcription, BUT: NGS much more versatile. Experimental design, data processing and data analysis are increasingly demanding.

9 RNA sequencing (RNA-seq) Gene expression

10 Gene expression [Wikipedia] Transcription: pre-mrna RNA processing: mrna - 5 capping - 3 cleavage and polyadenylation - RNA splicing Translation: protein - Ribosome contains rrna. - rrna depletion for RNASeq - Poly(A) enrichment for RNASeq

11 Gene expression [Licatalosi & Darnell, NRG 2010] - Exons, Introns and UTRs - DNA -> pre-mrna -> mrna - Poly-A tails and sites

12 RNA sequencing RNA sequencing

13 RNA sequencing Sequence cdna to get information about RNA content of the sample. Rough idea: Total RNA of sample Poly(A) selection to target mrna (not always) rrna removal (not always) Reverse transcription (mrna to cdna) Fragmentation Sequencing

14 RNA sequencing steps: Quality control and then: Computational pipeline on the left - Many different methods, but steps essentially similar. - Alignment - Disambiguation / filtering - Abundance quantification [Mortazavi et al., Nature Methods 2008] RPKM: reads (mrna abundance) kilo base (transcript length) million reads mapped (library size) FPKM: single -> paired end reads -> fragments more reliable alignment TPM: Transcripts per million

15 RNA sequencing steps: Assign reads (fragments) to annotated exons

16 RNA sequencing steps: Ambiguous reads: Alternative: constitutively expressed exons [gene-level only].

17 RNA sequencing steps: Make new gene(models):

18 RNA sequencing steps: Calculate gene abundances (here: RPKM)

19 RNA sequencing (differential) RNA abundance and what RNA-seq can do

20 RNA expression / abundance RPKM FPKM TPM [Pachter 2013]

21 Differential gene expression 2.0 First things first: Differential gene expression: Using constitutively expressed exons. Far from trivial. Easier than -> [Trapnell et al., Nature Protocols 2012]

22 Comparison to microarray => RNA similar to microarray; better correlation for higher expression levels. [Guo et al. 2013]

23 Differential isoforms with RNA-seq RNAseq can distinguish isoforms. [Katz et al. 2010]

24 SNPs from gene expression data RNAseq can detect polymorphisms. [Zhao et al. 2014]

25 RNA sequencing Characterization of entire transcriptome Can be used for: Discovery (existence / expression of) genes, exons, and transcripts (alternative splicing) Change detection between conditions in Gene expression Isoform expression, and more Different types of RNAs (e.g., mirnas, ncrnas)