I am working with a set of (bulk) RNA-Seq data collected across multiple runs, run at different times of the year. I have normalized my data using library size / quantile / RUV normalization, and would like to check (quantitatively and/or qualitatively) whether or not normalization has succeeded in removing the batch effects. It is important to note that by normalization has succeeded, I simply mean that unwanted variation has been removed - further analysis is required to ensure that biological variation has not been removed. What are some plots / statistical tests / software packages which provide a first-step QC for normalization?
Confirm success or failure RNA-Seq normalization
rna seq;normalization
You should use box plots and PCA plot. Let's take a look at the RUV paper:https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4404308/Before normalization and after UQ normalization:Libraries do not cluster as expected according to treatment. ... for UQ-normalized counts. UQ normalization does not lead to better clustering of the samples...Before normalization, the medians in the box-plot obviously look very different among replicates.After UQ normalization, the medians look closer but Trt.11 look like an outlier. Furthermore, the treatments aren't clustered on the PCA plot. Since they are replicates, you'd like them be close on the plot.After RUV normalization... RUVg shrinks the expression measures for Library 11 towards the median across libraries, suggesting robustness against outliers. ... Libraries cluster as expected by treatment. ...The RUV has made the distribution more robust and the samples closer on the PCA plot. However, it's still not perfect as one of the treatments is not close to the other two on the first PC.The vignettes for Bioconductor RUVSeq describes the two functions: plotRLE and plotPCA.