Question: Where to find a list of stress response genes?
gravatar for goodez
4 weeks ago by
United States
goodez420 wrote:

I need a list of heat shock genes to exclude from my analysis. I'm working with hg38 genome. I'm not sure the best approach to find all heat shock genes, or if I can download it from somewhere.


EDIT -- Not just heat shock, perhaps more generally I need "stress response genes".

annotation gene • 138 views
ADD COMMENTlink modified 4 weeks ago • written 4 weeks ago by goodez420

Stress is such a general term that it cannot be answered without more details. Stress can be the presence of heat or the presence of a chemotherapeutic agent or just the wrong pH in the incubator. I suggest you define what stress is in your context and then search NCBI for datasets where cells have been exposed to this stress and RNA-seq with proper experimental design (e.g. three biological replicates per stress and control condtion) was performed. Run it through a standard RNA-seq pipeline and extract the genes that come out differentially expressed.

ADD REPLYlink written 29 days ago by ATpoint11k

Thanks for the idea. I do agree it is general. However I really do not have more details. I was simply asked to "remove stress response genes" because the biology expert thinks one of our replicates was subjected to some kind of stress during the experiment.

Let's assume I just want to focus on heat shock genes. Is there really no existing source of annotated heat shock genes?

ADD REPLYlink written 29 days ago by goodez420

Don't allow that the wet-lab guys fool you. If they want things removed, ask them to clearly specify what they exactly want. If you eventually remove the wrong genes, you'll get 100% of the blame so ask them to be specific ;-)

ADD REPLYlink written 29 days ago by ATpoint11k

If one sample is suspected to have been subjected to uncontrolled factors, remove the whole sample, or remove nothing.

ADD REPLYlink written 29 days ago by h.mon21k

You can also do some diagnostics like PCA on the transformed counts (rlog or vst in DESeq2) and see if on the global level you see evidence for a stress exposure (that would be the respective samples clustering away from the unaffected ones of the same treatment group).

ADD REPLYlink written 29 days ago by ATpoint11k

I don't have counts. This is ChIP-seq data of Pol II. For the most part the replicates agree with each other, but they think the heat shock genes are activated in one replicate, which affects our average gene metaplot, for example. I was just hoping there was some source of heat shock genes existing... Doesn't need to be a perfect list of genes or anything.

ADD REPLYlink written 29 days ago by goodez420

Overall they correlate very well. It's not a huge concern for us if this is just limited to a small subset of genes. We aren't doing differential expression. We're looking at the general profile of Pol II transcription with and without a treatment. If I remove these genes, and still see the differences, then we'll throw out the replicate.

ADD REPLYlink written 29 days ago by goodez420
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: 1495 users visited in the last hour