Molecular profiling of RNA tumors using high-throughput RNA sequencing: From raw data to systems level analyses

Abstract

RNAseq is a powerful technique enabling global profiles of transcriptomes in healthy and diseased states. In this chapter we review pipelines to analyze the data generated by sequencing RNA, from raw data to a system level analysis. We first give an overview of workflow to generate mapped reads from FASTQ files, including quality control of FASTQ, filtering and trimming of reads, and alignment of reads to a genome. Then, we compare and contrast three popular options to determine differentially expressed (DE) transcripts (The Tuxedo Pipeline, DESeq2, and Limma/voom). Finally, we examine four tool sets to extrapolate biological meaning from the list of DE genes (Genecards, The Human Protein Atlas, GSEA, and ToppGene). We emphasize the need to ask a concise scientific question and to clearly under stand the strengths and limitations of the methods.