Question: Robust RNA-seq normalization
gravatar for rasmus.agren
12 days ago by
rasmus.agren0 wrote:

I have an experiment with significant technical issues. There are large differences in RIN values, as well as in QC metrics from the sequencing. Coverage is bad and varying, with samples ranging from 2 - 6 M uniquely mapping reads (human samples). Unfortunately the samples are few and expensive, and the experimental design requires that all samples have to be used. The purpose of the analysis is global differential expression analysis.

Which is the most robust normalization method to deal with this? I can assume that the amount of RNA per cell should be constant and I expect a relatively small number of DE genes. I don't expect any correlation between the biological question and the quality metrics, but I only have 24 samples and multiple conditions so that could occur by chance. I would typically use TMM normalization and then include RIN as a covariate in the DE testing. Previously I've used SVA for correction, but I'm worried about distorting the biological signal. Are there any better approaches to this?

rna-seq • 154 views
ADD COMMENTlink modified 10 days ago by kristoffer.vittingseerup1.7k • written 12 days ago by rasmus.agren0

Can you show some diagnostic plots? Do the samples really separate on a PCA based on RIN?

ADD REPLYlink written 12 days ago by ATpoint15k
gravatar for kristoffer.vittingseerup
10 days ago by
European Union
kristoffer.vittingseerup1.7k wrote:

The best way of dealing with it is to add those factors to the design matrix when you make your differential expression analysis model

ADD COMMENTlink written 10 days ago by kristoffer.vittingseerup1.7k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1645 users visited in the last hour