Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Variance modelling in DESeq2

    Hi there,

    My apologies if this is a stupid question. I've just moved into RNA-seq analysis, and am trying to figure out some issues.

    I'm working with a sub-optimal design (the sequencing was done prior to me coming on to this project). I have bone marrow that have been differentiated in two different medias. I have sequence data for three passages of each, and want to compare how close the differentiated cells look between the culture conditions. Effectively I think this means I no biological replicates (but with three technical replicates assuming the effect of passaging is negligible), which I know means I should be weary of any stats I get. I also have some controls samples with biological replicates. All of the samples were sequenced to a depth of ~30 million reads, except for one which had >200 million reads (assuming our sequencing vendor mis-quantified this sample).

    When I perform a PCA using DESeq2, all samples cluster pretty closely together, and away from the controls. Our collaborator performed the same analyses using the limma/Voom pipeline and found that the sample with >200 million reads did not cluster with everything else. This only seemed to occur after TMM normalization (without TMM, our data matches).

    I'm now not sure if my analysis is correct. Does DESeq2 normalize the data in a fashion similar to TMM? If so is this 'baked in', or should I be calling a particular function. I tried using estimateSizeFactors(dds), which adds a 'normalizationFactors' assay, but it doesn't seem to affect the rlogged data I'm using to generate the PCA.

    I also tried to recapitulate the variance plots (Fig 4. from the vignette), but sizeFactors(dds) returns NULL, and I'm not sure what to use to generate values for this attribute.

    These data are Salmon-aligned, and imported in to DESeq2 with:

    Code:
    txi <- tximport(files, type="salmon", tx2gene=t2g)
    dds <- DESeqDataSetFromTximport(txi, sampleTable, ~condition)
    dds <- estimateSizeFactors(dds)
    rld <- rlog(dds, blind=FALSE)
    plotPCA(rld, intgroup=c("condition", "passage" ))
    Any suggestions or help would be greatly appreciated!

    Thanks!

Latest Articles

Collapse

  • seqadmin
    Strategies for Sequencing Challenging Samples
    by seqadmin


    Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
    03-22-2024, 06:39 AM
  • seqadmin
    Techniques and Challenges in Conservation Genomics
    by seqadmin



    The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

    Avian Conservation
    Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
    03-08-2024, 10:41 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Yesterday, 06:37 PM
0 responses
7 views
0 likes
Last Post seqadmin  
Started by seqadmin, Yesterday, 06:07 PM
0 responses
7 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-22-2024, 10:03 AM
0 responses
49 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-21-2024, 07:32 AM
0 responses
66 views
0 likes
Last Post seqadmin  
Working...
X