This protocol is successfully used to quantitatively detect levels and spatial patterns of mRNA expression in multiple tissue types across vertebrate species. The method can detect low abundance transcripts and allows processing of hundreds of slides simultaneously. We present this protocol using expression profiling of avian embryonic brain formation as an example.