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
            Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
            by SEQadmin2



            Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

            Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
            ...
            10-14-2026, 09:55 AM
          • SEQadmin2
            New Genomics Technologies Take Aim at Long-Standing Limits
            by SEQadmin2


            Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

            We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
            ...
            09-28-2026, 10:25 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by SEQadmin2, 10-22-2026, 11:13 AM
          0 responses
          28 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-22-2026, 10:48 AM
          0 responses
          21 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-15-2026, 10:52 AM
          0 responses
          38 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-14-2026, 05:39 AM
          0 responses
          44 views
          0 reactions
          Last Post SEQadmin2  
          Working...