How to handle replicates and treatment conditions in TissueEnrich package
1
0
Entering edit mode
@andrewoler-19454
Last seen 5.9 years ago

Hi,

I'm trying to use the TissueEnrich package to determine which genes are differentially responsive to a drug treatment and I'm not quite sure how to get the data into a format that will work for TissueEnrich. It appears that TissueEnrich requires a table with TPM or FPKM for each tissue. However what I have is 4 biological replicates for each sample, so do you know how I would indicate that all of these samples are the same tissue? The other question (which is maybe more important) is can TissueEnrich work with differential expression results? For example, in our experiment we have treated each sample with either drug or vehicle and we're interested in knowing whether tissues behave differently in response to the drug. So our metric rather than being TPM or FPKM, which are absolute expression levels, would be a log fold ratio of expression in drug treatment / vehicle treatment, which will include some negative values. Will that work for TissueEnrich, or if not then do you have other suggestions for applying our data to TissueEnrich?

I read the vignette here and the manual but didn't see any reference to dealing with replicates or differential expression. In the vignette, it does say that for the example dataset, "In order to make the tissue-specific gene calculations more robust, we only used tissues that had ≥2 biological replicates." However it doesn't say how these replicates were handled. Maybe the data for the sample replicates was merged?

Thanks,

Andrew

rna-seq TissueEnrich tissue-specific expression differential expression • 1.7k views
ADD COMMENT
0
Entering edit mode

Hi Andrew,

As from your description, it seems like you are trying to find the differentially expressed genes based on the drug treatment. I recommend that you should use tools like DESeq2 or edgeR for that. Also, we took care of replicates by taking the average of the gene expression across replicates. TissueEnrich in its implementation only works with positive values. One way you can use TissueEnrich is to take the absolute fold-change values.

Regards, Ashish Jain

ADD REPLY
0
Entering edit mode
Ashish Jain ▴ 10
@ashish-jain-14944
Last seen 3.1 years ago
Iowa State University

Hi Andrew,

As from your description, it seems like you are trying to find the differentially expressed genes based on the drug treatment. I recommend that you should use tools like DESeq2 or edgeR for that. Also, we took care of replicates by taking the average of the gene expression across replicates. TissueEnrich in its implementation only works with positive values. One way you can use TissueEnrich is to take the absolute fold-change values.

Regards, Ashish Jain

ADD COMMENT
0
Entering edit mode

Thanks for the reply. Yes, I have already used DESeq2 to determine differentially expressed genes. I am looking for a tool to organize genes based on whether they are differentially expressed in one or more tissues. I can try your idea for dealing with replicates (average fold-change) and if I analyze up-regulated and down-regulated genes separately then I can use positive values for fold-change.

ADD REPLY

Login before adding your answer.

Traffic: 685 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