Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • Axolotl
    Junior Member
    • Nov 2013
    • 3

    Help understanding my DESeq data

    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!
    Last edited by Axolotl; 11-15-2013, 01:40 AM.
  • Simon Anders
    Senior Member
    • Feb 2010
    • 995

    #2
    Are you doing the qPCR validation on the same samples? If so, this is pointless -- you will just find what you already know. If not, you will soon learn why using raw instead of adjusted p values is not permitted,

    Comment

    • Axolotl
      Junior Member
      • Nov 2013
      • 3

      #3
      No, I'm validating on new samples, treated the same way though... So you are saying I could have picked the putative target genes blindly, by chance? Oh God...!
      What about the base mean values in my control sample? I can't figure out why these should vary...

      Comment

      • Simon Anders
        Senior Member
        • Feb 2010
        • 995

        #4
        You do know what a p value is, don't you? You wouldn't decide that a 10% threshold on a p value is a good idea without knowing what a p value is, or would you?

        Ok, seriously: Imagine you perform an RNA-Seq experiment, comparing samples treated with some drug with control samples, but some evil colleague has swapped your drug with water. How many of your organism's 20,000 genes will nevertheless have a p value below 10%?

        Answer: 10%, i.e., 2000 genes.

        This is the definition of "p value": Even if the treatment had no effect, x% of the genes will have a p value below x%.

        And this is why you must never use raw p values in high-throughput experiments.

        (And sorry that I sound annoyed but I explain this basic fact about once every week, and I cannot fathom why so many people don't know this, because -- it's somewhat important to know this, isn't it?)
        Last edited by Simon Anders; 11-15-2013, 06:18 AM.

        Comment

        • Axolotl
          Junior Member
          • Nov 2013
          • 3

          #5
          I am very sorry to annoy you. Yes, I new what the p Value is, although I don't know how the adjusted one is estimated. But shouldn't the adjusted pValue be somehow dependent on the raw value? And what about the fact, that these genes appear in DESeq and edgeR? Is this likely? And why is "my" bioinformatician not screeming?
          Again, I'm really sorry about my "stupid" questions, I already realized that my feeling of probabilities goes systematically against the laws.
          Thank you for being patient.

          Comment

          Latest Articles

          Collapse

          • SEQadmin2
            Nine Things a Sample Prep Scientist Thinks About Before Sequencing
            by SEQadmin2


            I’m not a sequencing expert. I’m a purification scientist who uses NGS to evaluate workflows my group develops. With this perspective, we think about the sample first and the NGS workflow second. The sequencer is an exceptionally honest reporter, but it can only report on what you give it, so whether you get clean, interpretable data from an NGS workflow is largely determined before you begin.

            Here are nine questions we think about, in roughly the order they matter, before...
            06-18-2026, 07:11 AM
          • SEQadmin2
            From Collection to Sequencing: Why Sample Preparation and Preservation Define Sequencing Data
            by SEQadmin2


            Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.


            The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
            ...
            06-02-2026, 10:05 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by SEQadmin2, Today, 11:10 AM
          0 responses
          5 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 06-17-2026, 06:09 AM
          0 responses
          41 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 06-09-2026, 11:58 AM
          0 responses
          102 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 06-05-2026, 10:09 AM
          0 responses
          123 views
          0 reactions
          Last Post SEQadmin2  
          Working...