Is there any straightforward way to identify housekeeping genes based on a series of genome-wide gene expression profiles? Literature searches suggest that meta-analysis of large-scale expression data can do the task, but what if the organism doesn't have enough expression data to mine? Is it safe to grab the genes with fold changes bouncing around one (0.9~1.1 range) if comparing all samples pairwise?
I think this answer depends on why you need to know the housekeeping genes. What is the plan for using them once you have them?
The main purpose is to identify a set of housekeeping genes to help normalize mRNA-Seq expression profiles across tissues and experiments. Much like what TMM normalization does, we'd like to use a set of housekeeping genes to calibrate or accommodate the RNA composition and/or sequencing depth difference across samples.