Question: suggestions to compare gene counts
gravatar for b10hazard
5 weeks ago by
United States
b10hazard20 wrote:

I'm new to RNA-seq and I'm playing around with STAR Aligner. I was changing some of the command line options on STAR and trying to determine their impact on my gene counts (counted using featureCounts command line tool). My question is, what is the best way to compare these data sets?

For example, in one trial I ran STAR with default settings and then I reran with alignIntronMin set to 100 (default is 21). So now I have two count text files...



I came up with a metric where I calculate the fold-difference in gene counts between sets and calculate the median fold-difference of the 500 most changed genes. This can be called something like Median-500-Most-Fold-Different-Genes metric. In this case the value is really low; 0.0000001, probably because my command line change had little effect.

However, is there a better way to assess these changes? Am I on the right track? Keep in mind I only have one RNA-seq dataset. I can't do comparisons beyond how my command line changes alter the gene counts. Any suggestions would be welcome. Thanks!

ADD COMMENTlink modified 5 weeks ago by Charles Plessy2.4k • written 5 weeks ago by b10hazard20
gravatar for Charles Plessy
5 weeks ago by
Charles Plessy2.4k
Charles Plessy2.4k wrote:

The field of differential expression analysis is quite mature now, so I recommend to first look at existing solutions, that are based on one or two decades of experience. In case you do not know where to start and you are fine using R and Bioconductor, I would suggest to have a look at their RNA-seq workflow.

ADD COMMENTlink written 5 weeks ago by Charles Plessy2.4k

As I mentioned in my post; I only have one file. So I cannot look at differential gene expression.

ADD REPLYlink written 5 weeks ago by b10hazard20

You have only one file, but with it, you generated multiple expression tables. Have you tried to compare the expression tables using software for differential gene expression ?

ADD REPLYlink written 5 weeks ago by Charles Plessy2.4k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1333 users visited in the last hour