Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Differential gene expression analysis with bioreplicates using EdgeR/DESeq

    Hi everyone,

    I have read and search a lot about this topic but can not find any solution to my problem. May you will be able to help me.

    I am doing an intern-ship in bioinformatics for my master and I have to deal with RNA-seq data. I have 2 sets of experiments (A and B), both having 2 illumina runs of two stages (1 and 2) of a plant. A and B has not been done at the same time and the technology is a bit different, coming up with:
    runs about 30M reads for A,
    runs about 80M reads for B.

    For a given stage the log(RPKM) of the replicates are very well correlated.

    When I use EdgeR to obtain a common dispersion from the counts of each runs searching for differential expressed genes between each stage I obtain 0.86. Which seems far too big regarding the correlation of the RPKM. Moreover the number of differentially expressed genes is not consistent with our affymetrix knowledge (about 250 genes when we expected about 1000 genes).

    I first think about filtering the list of genes from the one having a count per million below 1 in all conditions. I then obtain a dispersion of 0.76 : still to high...

    I also think about getting variance stabilized data (with DESeq) to use with limma but it does not make sense if the samples are not paired, does it?

    I am wondering if I am doing something wrong here and if there are any filtration/computation that I should have done to obtain a more consistent common dispersion.

    Any idea would be really appreciate,

    François

  • #2
    A few points:

    edgeR is a Bioconductor package, so more detailed help is available on the Bioconductor mailing list than on SEQanswers.

    If you want to get your RNA-seq data into limma, the way to do this is use the voom() function of the limma package. See the limma User's Guide.

    There are any number of things that might be causing problems with your analysis, but there's no to way know from the information that you give. Your dispersion values are very high indeed. Have you used an MDS plot to look at your data?

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Exploring the Dynamics of the Tumor Microenvironment
      by seqadmin




      The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
      07-08-2024, 03:19 PM
    • seqadmin
      Exploring Human Diversity Through Large-Scale Omics
      by seqadmin


      In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
      06-25-2024, 06:43 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 07-19-2024, 07:20 AM
    0 responses
    38 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 07-16-2024, 05:49 AM
    0 responses
    49 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 07-15-2024, 06:53 AM
    0 responses
    61 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 07-10-2024, 07:30 AM
    0 responses
    43 views
    0 likes
    Last Post seqadmin  
    Working...
    X