Hi everyone,
I have a RNA-seq expression count matrix of 2 contrast conditions (~10 biological samples per condition), but these conditions are affected by (severe) batch effect from different sequencing experiments.
I looked up for some batch-effect removal tools, but they could only fix batch-effect for samples of same condition group (different conditions may contain large true biological variations that account for most of batch-effect difference).
I plan to choose a group of housekeeping genes to adjust for this group difference, but I am still confusing about practical steps to do that. Could you please give me some suggestions? Here are some thoughts I am still questioning:
- Should I perform TPM then TMM cross-sample normalization before considering expression value of these housekeeping gene?
- In these housekeeping genes, there is probably a large difference in expression value between them, how could I straighten all of them down to one scaling factor for each sample, and then scale expression level of all other genes by this factor?
Thank you very much.
rna-seq
normalization
What separates these "batch groups"? Different library preps, sequencing machine?
Yes, because I obtained them from different public databases, so 2 experiment protocols were totally different (however both data are Hiseq Illumina reads)
Then I see little chance of using them in the same analysis, especially because you have absolutely no way of validating the results by an independent experimental approach.