Differential Expression of Single Gene
1
0
Entering edit mode
5.8 years ago
gtasource ▴ 60

Hi all,

I have 6 aligned BAM files for a species I am working with. We have a single gene we want to see if it's differential. I'd rather not have to run an entire DE analysis when we just want to see if one specific gene is differential. My plan was to extract the gene coordinates and run BedTools coverage to get raw counts. From there, I'd do some type of statistical test (not sure which one) to see if there is a significant difference in the counts.

expression gene expression de-seq • 1.4k views
ADD COMMENT
2
Entering edit mode

This does not really sound like a good approach to me.

For any DEG analysis, the mapping step is the largest part of the analysis and since you apparently already have that I don't think you will gain much be filtering before running DEG. It might even be so that your statistics will turn out unreliable (eg. the library normalisation).

Personally I would run the whole DEG and only then select your gene of interest.

ADD REPLY
1
Entering edit mode

How were the BAMs aligned? I would not run any statistical test on raw counts. Why not take the opportunity to do a 'robust' (yep, that magical word again...) differential expression analysis for all genes and have that resource for future purposes? What if, in 1 month, you want to go back and check another gene?

With an aligned BAM from, e.g., STAR, you could run featureCounts or HTseq for the purposes of obtaining raw counts over transcripts listed in an input GTF. You could then normalise these raw counts in EdgeR or DESeq2.

ADD REPLY
0
Entering edit mode

this could have been a (bad though) solution for the problem, but you have to normalize by sequencing depth. Otherwise the counts for 1 gene don't tell you anything.

ADD REPLY
0
Entering edit mode
5.8 years ago

Best is to run differential analysis on all the data and check for your gene of interest in the DE results.

ADD COMMENT

Login before adding your answer.

Traffic: 2234 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6