In my lab we have a question about this (we are not experts in bioinformatics by the way, we still learning):
We are doing rna-seq of resistant vs susceptible peppers. I did a DEG analysis using DESEQ2, I got up and down regulated genes using a FDR:0.01 and LFC:0.5. Now I have my respective list of DE genes, but a questions springed to us: with our information is possible to do a "contrast" of genes expression levels using as "control" some house keeping genes ? Like a RT-PCR.
1 answer
Yes, it can be easily done. If you have a list of genes that you assume have the same expression between the two conditions you can use them as controlGenes in the estimateSizeFactors function call. See https://www.rdocumentation.org/packages/DESeq2/versions/1.12.3/topics/estimateSizeFactors:
estimateSizeFactors"(object, type = c("ratio", "iterate"), locfunc = stats::median, geoMeans, **controlGenes**, normMatrix)
If you have a list of genes that you assume have the same expression between the two conditions
I am really interested in this. Can you do such assumptions on non model organisms (peppers in this case) without an experimental validation on the same samples used for the DEG analysis?
You have to make an assumption eventually, either implicit or explicit. The implicit assumption is that most of the genes have the same expression level. Well, if it doesn't hold then the conditions you are comparing might be too distinct. An explicit assumption using a list of genes might be beneficial especially if you have a lot of experience and know what to expect for.
Thank you for the answer!
Log in to answer this question.
-
Limma for DEG analysis
written by rahma.khoualdia •
-
How to approach DESeq2 analysis on 2 plant lines infected with virus, with 2 time points (days post…
written by duplessisantzel •
-
Why we are using filtering >0 for up and <0 for down after TopTags() to extract de genes ids?
written by Ann •
-
Detect RNAi-derived transcripts in RNA-seq
written by boscariolfp •
-
limma - voom with TPM data
written by AlexStar •
-
best practice for applying FDR to highly similar group DEG analysis
written by hs960201 •
-
Which is the best approach to get DEG across one condition, in a three-factor design experiment?
written by antonioggsousa
-
TCGA RNA-Seq data for gene correlation test
written by JackieMe •
-
What Sample Purity Is Needed For An Accurate Gene Expression Assay
written by Obi Griffith
-
Too Many Genes Under Positive Selection
written by Plantae
