We present a method and protocol for fluorescent in situ hybridization (FISH) in zebrafish embryos to enable three-dimensional imaging of patterns of gene expression using confocal laser scanning microscopy. We describe the development of our protocol and the processing workflow of the three-dimensional images from the confocal microscope. We refer to this protocol as zebraFISH. FISH is based on the use of tyramide signal amplification (TSA), which results in highly sensitive and very localized fluorescent staining. The zebraFISH protocol was extensively tested and here we present a panel of five probes for genes expressed in different tissues or single cells. FISH in combination with confocal laser scanning microscopy provides an excellent tool to generate three-dimensional images of patterns of gene expression. We propose that such three-dimensional images are suitable for building a repository of gene expression patterns, complementary to our previously published three-dimensional anatomical atlas of zebrafish development (bio-imaging.liacs.nl/). Our methodology for image processing of three-dimensional confocal images allows an analytical approach to the definition of gene expression domains based on the three-dimensional anatomical atlas.