Question: Co-expression analysis with paired samples
2
gravatar for agatawesol
5.4 years ago by
agatawesol30
United States
agatawesol30 wrote:

Dear all,

I have whole transcriptome RNA-Seq data for 20 paired samples (10 controls and 10 stimulated with an agent). I have ~2,500 differentially expressed genes between the controls and stimulated samples, and I would like to identify some smaller modules of genes acting together within this big number of genes. Ideally those would be then matched with a common upstream transcription factor, or other regulator. I would like to perform a co-expression analysis on my differentially expressed genes, but I am not sure what is the best way to do this with paired samples. 

1) What is the best input data for calculating correlations and co-expression? I have considered: the read counts mapped to my genes (that would not take into account pairing between samples), the difference in read counts between my stimulated sample and its control (one value per sample pair) or log2 fold change between the pair. Maybe there is an even better solution? 

2) Which correlation metric would be better to use with this small sample size? Pearson or Spearman?

3) When clustering my genes - is it better to consider genes correlated in one direction only, or both (positively and negatively correlated)?

Any help would be greatly appreciated!

ADD COMMENTlink modified 5.4 years ago • written 5.4 years ago by agatawesol30

I agree with Ashutosh's suggestion of WGCNA, you can get normalized counts (CPM) from EdgeR and DESeq(2), that you can use for WGCNA.

You could also do hierarchical clustering to see which genes cluster together. It will handle positive/negative.

ADD REPLYlink written 5.4 years ago by brentp23k
1
gravatar for Ashutosh Pandey
5.4 years ago by
Philadelphia
Ashutosh Pandey11k wrote:

Have you ever tried WGCNA (http://labs.genetics.ucla.edu/horvath/CoexpressionNetwork/Rpackages/WGCNA/)

It is a well established tool for generating coexpression network. You should definitely give it a try. 

ADD COMMENTlink written 5.4 years ago by Ashutosh Pandey11k
0
gravatar for agatawesol
5.4 years ago by
agatawesol30
United States
agatawesol30 wrote:

Yes, WGCNA is definitely on my list of things to try with this data. I didn't see if it can somehow deal with paired samples?

My question is more fundamental than which software to use. I would like to know: how to best utilize the power of paired design for co-expression analysis? And what should I be correlating here - the counts or some measure of change of expression in a pair?

ADD COMMENTlink written 5.4 years ago by agatawesol30

You should make this a comment rather than a new answer see the [ADD COMMENT] box below any existing comment or reply.

ADD REPLYlink written 5.4 years ago by brentp23k

Sorry, will do that next time. 

ADD REPLYlink written 5.4 years ago by agatawesol30
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: 1122 users visited in the last hour