Unconfigured Ad

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

    #1

    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
            Beyond CRISPR/Cas9: Understand, Choose, and Use the Right Genome Editing Tool
            by SEQadmin2



            CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).

            Despite this, “CRISPR helped turn genome editing from a specialized technique into
            ...
            07-31-2026, 11:01 AM
          • SEQadmin2
            Proteomic Platforms: How to Choose the Right Analytical Strategy to Improve Detection and Clinical Applications
            by SEQadmin2


            Proteomics platforms are evolving rapidly, with advances in mass spectrometry and affinity-based approaches expanding what researchers can detect and at what scale. As the field moves toward deeper proteome coverage and clinical applications, scientists face an increasingly complex landscape of tools. This article will explore how researchers are navigating these choices to find the right platform for their work.

            The systematic characterization of the human proteome has
            ...
            07-20-2026, 11:48 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by SEQadmin2, Today, 10:05 AM
          0 responses
          8 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 08-13-2026, 12:22 PM
          0 responses
          32 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 08-11-2026, 10:35 AM
          0 responses
          24 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 08-06-2026, 07:41 AM
          0 responses
          38 views
          0 reactions
          Last Post SEQadmin2  
          Working...