Question: Robust RNA-seq normalization
gravatar for rasmus.agren
9 months 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 • 356 views
ADD COMMENTlink modified 9 months ago by kristoffer.vittingseerup3.0k • written 9 months ago by rasmus.agren0

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

ADD REPLYlink written 9 months ago by ATpoint28k
gravatar for kristoffer.vittingseerup
9 months ago by
European Union
kristoffer.vittingseerup3.0k 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 9 months ago by kristoffer.vittingseerup3.0k
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: 709 users visited in the last hour