Hi everyone,

I am a PhD student in biochemistry, and I am learning about gene expression signature. My lab generated a 36-gene mouse signature. These genes are all highly expressed. I am interested in identifying "mouse-like" human samples from a large set of primary breast tumors.

I was wondering if someone could please give me the general guidelines on how to apply a gene signature. I can write code, but don't understand the principles (I am reading tough). Is it based on the gene names and their fold-change or just the names? I am sorry for asking such a basic question, but I am learning this aspect of bioinformatics. I read that a naive Bayes classifier is a good idea? Alternatively, ranking the samples (based on how well they express the signature) and using bootstrap resampling?

I would also greatly appreciate to be redirected to a former post or tutorial.

Thank you!
One possibility (that wouldn't even require writing much code) is to use this signature as a gene set and use GSEA on the human samples to look for samples in which that set is more highly expressed than expected. The general idea is to perform GSEA on a large number of samples, many of which you expect to not show enrichment, and then look at the resulting enrichment score (or p-value) distribution. From that, you should be able to get an idea of whether the expression of this set generally follows a normal distribution or whether there's a bimodal distribution...meaning that there's a subset of samples that you're going to be very interested in. You could alternatively use resampling there, though I think it'll be quicker and easier to just have a look at the distributions first (nothing is preventing you from doing both).

That's one fairly straight forward possibility, though there are others.