Question: DGE analysis suggestion : linearity between TF and target gene expression level
0
gravatar for bharata1803
4.0 years ago by
bharata1803420
Japan
bharata1803420 wrote:

Hello,

So, I have done gene expression analysis using DESeq2. Now, I have several gene which is differentially expressed between normal and cancer sample. I tried to analyze what is the cause of the differentially expressed. I also have data about variation from exome seq  with same sample (I got it from NCBI GEO). To find the cause of DE gene, I tried to list all of TF for that gene that I can found from gene card website. My method is really simple, use Pearson correlation too check the correlation coeff and plot it in scatter plot to find the linearity. From this method, I found quite interesting result which is I found 1 TF correlates strongly with downregulation of my target gene (the target gene and TF gene both downregulated). Other TF didn't show this result. I use assay function from DESeq2 object to get the normalized readcount, not the raw count or logarithmic expression level.

My question is, is my method acceptable, both biologically and mathematically?

Also, I found most of the TF and target gene didn't give any strong linear correlation, do I need to calculate it with non-linear correlation method?

 

Thank you for your suggestion.

rna-seq • 1.3k views
ADD COMMENTlink modified 4.0 years ago by Devon Ryan90k • written 4.0 years ago by bharata1803420

what did you calculate the Pearson correlation for? i.e. which values are you comparing with each other?

ADD REPLYlink written 4.0 years ago by Friederike4.2k

I use the value from DESeq2 object I get it using assay function from DESeq2. I use Pearson correlation coefficient just to check which TF has strong correlation in linearity befor I plot it in scatter plot.

ADD REPLYlink written 4.0 years ago by bharata1803420

and what value did you use for the TF?
 

ADD REPLYlink written 4.0 years ago by Friederike4.2k

Same, I compare expression level from DESeq2 for target gene and TF

ADD REPLYlink written 4.0 years ago by bharata1803420
0
gravatar for Devon Ryan
4.0 years ago by
Devon Ryan90k
Freiburg, Germany
Devon Ryan90k wrote:

While a non-linear relationship is probably more biologically realistic, if you get a good candidate using a simple linear relationship then that's a good start. You'll obviously need to do more experiments (e.g., is it expressed in your cells of interest, what happens if you inhibit/activate it in your system, ...), but you seem to have a good start on things. BTW, keep in mind that when you test many candidates, you'll randomly get some significant findings simply due to testing many things (i.e., always be wary of multiple comparisons).

ADD COMMENTlink written 4.0 years ago by Devon Ryan90k

Thank you for your reply. I've read from another post that Spearman correlation also god  for checking because it can interpret a bit non linearity. 

ADD REPLYlink written 4.0 years ago by bharata1803420

Yeah, I would agree that a rank-based method like Spearman's correlation would be preferred.

ADD REPLYlink written 4.0 years ago by Devon Ryan90k
Please log in to add an answer.

Help
Access

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