Here is the scenario: I have RNAseq data from seven tumor samples. I've put them through the cufflinks pipeline including both cuffdiff and cuffnorm. For cuffdiff, I treated each sample as its own condition with no replicates. I'm aware this cannot produce statistically significant results, but beyond that differential expression isn't really what I want anyway as all these samples derive from the same disease.
Broad question: Is there anything meaningful that can be done in the given scenario?
More specific:
I'm grasping at straws here because of the lack of planning in this experiment (there is normal and tumor DNA but not matched RNA!). Thanks for your help.
Broad question: Is there anything meaningful that can be done in the given scenario?
More specific:
- If I have a list of a few genes that I'm interested in. As I understand, I cannot simply compare (within each sample separately) the FPKM of a gene of interest with other random genes (or e.g. an average across all genes). Can I do anything with this list?
- I have access to other RNAseq data from public databases. Does it make sense (for exploratory purposes anyway) to do differential expression between my tumors and e.g. breast cancer tumors? Given tissue- and individual-dependency of expression, I'm guessing not, but I'm wondering if any of you had more insight.
I'm grasping at straws here because of the lack of planning in this experiment (there is normal and tumor DNA but not matched RNA!). Thanks for your help.
Comment