Question

Programs used in RNA-seq workflow

0

Entering edit mode

cirisarri95 • 0

@cirisarri95-10027

Last seen 8.6 years ago

Good morning everyone,

I was interested in knowing which programs is the bioconductor package RNA-seq workflow: gene exploratory analysis and differential expression for:

-Quality control

-Mapping

-Counting

-Read Analysis

Thanks everyone!

Regards

deseq2 bioconductor workflows differential gene expression rnaseqgene • 3.3k views

ADD COMMENT • link 8.7 years ago cirisarri95 • 0

score 1 · Answer 1 · 2016-04-03

1

Entering edit mode

Michael Love 43k

@mikelove

Last seen 2 days ago

United States

hi,

I'm not sure I follow your question.

There are a number of R packages used in the workflow.

Every time you see code, such as library("Rsamtools"), that is loading an R package called, "Rsamtools". We describe each package used in each respective section, and when there are multiple options, we describe the one we use as well as provide links to alternatives (counting, differential expression analysis, distance calculation, etc.).

You can just search the page here for "library":

http://www.bioconductor.org/help/workflows/rnaseqGene/

And at the end of the page you can see all the libraries which were loaded throughout the R session, under "other attached packages".

Additionally we describe the use of the STAR read aligner in the workflow, with citation and a link to the STAR aligner software. This software is not an R package, but a standalone tool which is run on the command line.

ADD COMMENT • link 8.7 years ago Michael Love 43k

0

Entering edit mode

Hi,

Then, if I want to run the whole package workflow, do I need to first run STAR read aligner? Or will it be runned with the rest of the package?

Thanks.

ADD REPLY • link 8.6 years ago cirisarri95 • 0

0

Entering edit mode

STAR is not actually part of the workflow which is run "live". The live workflow contains only commands in R. We just show you example commands that were used for STAR to make the BAM files used in the workflow. You should prepare the alignment yourself using whatever alignment software you prefer to generate the BAM files.

ADD REPLY • link 8.6 years ago Michael Love 43k

score 0 · Answer 2 · 2016-04-05

0

Entering edit mode

cirisarri95 • 0

@cirisarri95-10027

Last seen 8.6 years ago

Hi,

Thanks for your help.

Also, for running this workflow, is it set up so that I input data once and it runs the whole process with no more help? Or should I be inputting for each step the data again and again?

ADD COMMENT • link 8.7 years ago cirisarri95 • 0

1

Entering edit mode

You only build the DESeqDataSet (dds) object once, and then the different functions do something and either return the object with new information stored inside of it, or they return a different object.

You'll have to just look at the workflow and follow along. You'll see by the names of the R variables when things are created, when they are modified and when you get a new object as an output of a function.

If you are brand new to R, you may benefit from taking a beginner's R course. I have some R resources collected here:

http://genomicsclass.github.io/book/pages/resources.html

ADD REPLY • link 8.7 years ago Michael Love 43k