Hello everybody!
I hope somone can help me... I am really bad in this bioinformatics stuff.
I've made an RNAseq experiment, among other goals, to get target genes regulated by a chromatin protein complex:
- 1 control, 3 different mutants, 3 biological replicates.
- Illumina sequencing > reads mapped with TopHat (by a statistician)
- Differential gene expression analysis with DESeq and edgeR (by a bioinformatician). I got 2 lists (DESeq + edgeR) of differentially expressed genes for each mutant compared to the control.
- Because of the small number of significantly deregulated genes (according to adjusted pValue or FDR), I've set a rather loose threshold for the pValue (10%) and the log2 fold change (-1.5 > log2 FC > 1.5) in the two lists. Then I merged them, ending up with lists of about 100-400 deregulated genes per list.
- I am now validating some interesting Targets by qPCR. Further biological experiments are planned.
My questions:
1. DESeq data: Is it normal that the "base mean" and the "base mean A" (corresponding to the control sample, which was the same for the 3 treatments) are not the same in the 3 lists (for the 3 mutants) of differentially expressed genes?
2. qPCR validation: What is the relationship between the normalized expression ratio of the "2^-deltadeltaCt" method and the fold change of the DESeq/edgeR analysis. Should there be a quantitative correlation?
Thank you for any help!
I hope somone can help me... I am really bad in this bioinformatics stuff.
I've made an RNAseq experiment, among other goals, to get target genes regulated by a chromatin protein complex:
- 1 control, 3 different mutants, 3 biological replicates.
- Illumina sequencing > reads mapped with TopHat (by a statistician)
- Differential gene expression analysis with DESeq and edgeR (by a bioinformatician). I got 2 lists (DESeq + edgeR) of differentially expressed genes for each mutant compared to the control.
- Because of the small number of significantly deregulated genes (according to adjusted pValue or FDR), I've set a rather loose threshold for the pValue (10%) and the log2 fold change (-1.5 > log2 FC > 1.5) in the two lists. Then I merged them, ending up with lists of about 100-400 deregulated genes per list.
- I am now validating some interesting Targets by qPCR. Further biological experiments are planned.
My questions:
1. DESeq data: Is it normal that the "base mean" and the "base mean A" (corresponding to the control sample, which was the same for the 3 treatments) are not the same in the 3 lists (for the 3 mutants) of differentially expressed genes?
2. qPCR validation: What is the relationship between the normalized expression ratio of the "2^-deltadeltaCt" method and the fold change of the DESeq/edgeR analysis. Should there be a quantitative correlation?
Thank you for any help!
Comment