Whole-body gene expression pattern registration in Platynereis larvae
© Asadulina et al.; licensee BioMed Central Ltd. 2012
Received: 1 August 2012
Accepted: 23 October 2012
Published: 3 December 2012
Skip to main content
© Asadulina et al.; licensee BioMed Central Ltd. 2012
Received: 1 August 2012
Accepted: 23 October 2012
Published: 3 December 2012
Digital anatomical atlases are increasingly used in order to depict different gene expression patterns and neuronal morphologies within a standardized reference template. In evo-devo, a discipline in which the comparison of gene expression patterns is a widely used approach, such standardized anatomical atlases would allow a more rigorous assessment of the conservation of and changes in gene expression patterns during micro- and macroevolutionary time scales. Due to its small size and invariant early development, the annelid Platynereis dumerilii is particularly well suited for such studies. Recently a reference template with registered gene expression patterns has been generated for the anterior part (episphere) of the Platynereis trochophore larva and used for the detailed study of neuronal development.
Here we introduce and evaluate a method for whole-body gene expression pattern registration for Platynereis trochophore and nectochaete larvae based on whole-mount in situ hybridization, confocal microscopy, and image registration. We achieved high-resolution whole-body scanning using the mounting medium 2,2’-thiodiethanol (TDE), which allows the matching of the refractive index of the sample to that of glass and immersion oil thereby reducing spherical aberration and improving depth penetration. This approach allowed us to scan entire whole-mount larvae stained with nitroblue tetrazolium/5-bromo-4-chloro-3-indolyl phosphate (NBT/BCIP) in situ hybridization and counterstained fluorescently with an acetylated-tubulin antibody and the nuclear stain 4’6-diamidino-2-phenylindole (DAPI). Due to the submicron isotropic voxel size whole-mount larvae could be scanned in any orientation. Based on the whole-body scans, we generated four different reference templates by the iterative registration and averaging of 40 individual image stacks using either the acetylated-tubulin or the nuclear-stain signal for each developmental stage. We then registered to these templates the expression patterns of cell-type specific genes. In order to evaluate the gene expression pattern registration, we analyzed the absolute deviation of cell-center positions. Both the acetylated-tubulin- and the nuclear-stain-based templates allowed near-cellular-resolution gene expression registration. Nuclear-stain-based templates often performed significantly better than acetylated-tubulin-based templates. We provide detailed guidelines and scripts for the use and further expansion of the Platynereis gene expression atlas.
We established whole-body reference templates for the generation of gene expression atlases for Platynereis trochophore and nectochaete larvae. We anticipate that nuclear-staining-based image registration will be applicable for whole-body alignment of the embryonic and larval stages of other organisms in a similar size range.
Three-dimensional digital anatomical atlases represent biological structures within a reference space, and can provide morphological, molecular and functional information for the components of these structures. Atlases are used to represent specimens that possess the same set of morphological features, such as embryos of the same developmental stage or adult brains with stereotypic morphology [1, 2]. Such atlases can visualize complex anatomical relations or provide a reference for inter-individual and inter-species comparisons (for example, sexual dimorphism) [3, 4]. The generation of atlases requires volume imaging using a marker that visualizes the entire structure of a specimen (for example, DNA stains, immunostaining of all synapses).
Atlases show their full potential when data can be compiled across several individuals labeled or manipulated in various ways. The mapping of additional features to the reference template requires co-labeling of the specimens with the reference marker and the markers of interest. Markers can include neuron-filling dyes labeling individual neurons [5–7], fluorescent proteins expressed in a subset of cells by transgenic techniques , or gene expression patterns visualized by in situ hybridization [9–11], immunostaining , or transgenic techniques . Digital atlases overcome the limitations of traditional single- or multicolor labeling light microscopy methods that are often restricted in the number of channels (genes, cells) that can be simultaneously detected (an exception is the Brainbow technology ).
When the morphology of individual neurons is projected onto a reference template, such a representation can give insights into the organization of neural circuits and their influence on behavior [3, 5, 6, 8]. The representation of multiple gene expression patterns in a common reference can elucidate gene regulatory interactions , inter-individual variation in gene expression , or reveal the ‘molecular fingerprint’ of neuron types, allowing evolutionary comparisons . Such gene expression atlases have been generated for Drosophila, mouse , chicken , zebrafish  and the marine annelid Platynereis.
The integration of imaging data acquired across multiple individuals and experiments requires sophisticated image registration techniques. Simple superimposition would be insufficient due to variations in morphology and experimental conditions. The image registration methods can be classified as ‘voxel value-based’ [9, 11, 15] or ‘segmentation-based’ [1–3, 7, 10]. In the voxel value-based methods, the registration procedure relies on a metric calculated from pixel intensities (for example, mutual information ). For the registration procedure, the stained axonal scaffold or neuropil [8, 11], nuclear stain , or differential interference contrast (DIC) images  can be used as a reference signal. Segmentation-based registration methods require the prior segmentation and annotation of objects . The corresponding segmented objects from the individual images then determine the transformations that should be applied. For segmentation, either cells expressing a marker gene , stained nuclei combined with a transgenic muscle label , or other anatomical structures  can serve as a reference.
Image registration protocols usually start with coarse, rigid registration that is further refined by non-rigid registration. For example, the virtual insect brain (VIB) protocol  begins with global and local rigid registration, refined by non-rigid registration. This protocol was used to generate brain atlases for the fruit fly (Drosophila melanogaster) , the sphinx moth (Manduca sexta) , the flour beetle (Tribolium castaneum) , and the desert locust (Schistocerca gregaria) . Another image registration protocol, the iterative shape-averaging method (ISO) , starts with affine registration followed by iterative non-rigid registration. This protocol was used to generate atlases for the honeybee (Apis mellifera) , the desert locust (S. gregaria)  and the tobacco budworm (Heliothis virescens) .
A similar protocol was applied for the registration of the gene expression patterns in Platynereis[11, 15]. This voxel value-based method applied rigid and affine transformations followed by B-Spline deformable transformation, as implemented in the Insight Toolkit (ITK), an open-source tool for image analysis . The reference signal used in this protocol was the larval axonal scaffold and ciliary bands stained with an acetylated-tubulin antibody.
In Platynereis, image registration has to date been restricted to the trochophore (48 hours post fertilization, hpf) larval episphere, due to the limited depth-penetration achieved by using imaging setups with refractive index mismatches (for example, glycerol mounting with an oil immersion objective). Here we push confocal microscopy closer to its depth-limit by using the mounting medium 2,2′-thiodiethanol (TDE) with a refractive index matched to that of glass and immersion oil . The use of TDE enabled us to perform high-quality whole-body scans on Platynereis trochophore (48 hpf) and nectochaete (72 hpf) larvae using isotropic voxel size. Isotropic voxel size also alleviated the need for precise body orientation, allowing us to scan and register larvae in any orientation. These whole-body scans were first oriented along their antero-posterior (AP) and dorso-ventral (DV) axes based on anatomical landmarks using a problem-specific algorithm developed in ImageJ . The oriented larval images were then registered with ITK using either the nuclear stain DAPI, or acetylated tubulin, as the reference signal. We developed reference templates for 48 and 72 hpf larvae and demonstrated, using several marker genes, that the registration method allows near-cellular-resolution registration. We also compared the accuracy of the different templates, finding that the nuclear-stain template often outperforms the tubulin template. We provide detailed instructions and scripts to allow the use of our registration protocol in Platynereis, and to facilitate its transfer to other organisms. This work provides the foundation for the development of a near-cellular-resolution whole-body gene expression atlas for Platynereis larvae.
In situ hybridization using nitroblue tetrazolium (NBT)/5-bromo-4-chloro-3-indolyl phosphate (BCIP) staining combined with anti-acetylated-tubulin and DAPI staining was performed as previously described , with slight modifications ([see Additional file 1 and ). Following in situ hybridization, larvae were transferred into 97% TDE diluted with PBS plus 0.2% Tween. Larvae were placed on a glass slide with several layers (three to five) of adhesive tape on both sides, forming a chamber with the coverslip.
Imaging was performed with an Olympus FV1000 confocal system on an IX81 inverted microscope using an UPlanSApo 60× Oil objective with 1.35 N.A. and 0.15 mm working distance (Olympus Deutschland GmbH, Hamburg, Germany). For excitation we used 405, 559 and 635 nm diode lasers. Three-channel excitation was performed simultaneously and emission light passed through a 405/488/559/635 main dichroic beam splitter and through a pinhole of diameter 120 μm before going to the detection module. Fluorescent light was split into three beams directed to three PMTs by the secondary dichroic filters 490 and 640. Channel 1 for DAPI detection had a further 425 to 475 band-pass filter, channel 2 for TRITC detection had a 570 to 625 band-pass filter and channel 3 for the detection of NBT/BCIP had a 780 long-pass filter. Images were scanned with 1.0× zoom and recorded as 512 × 512 pixels corresponding to a pixel size of 0.414 μm × 0.414 μm. The entire volume of the larvae was scanned with a z-step size of 0.41 μm, in order to have an isotropic voxel size. This pixel size is above the X-Y optical resolution (approximately 0.23 μm for our system), but is small enough to provide good resolution, tolerable scanning times and minimal bleaching. We chose to use an isotropic voxel size to allow the scanning of larvae in any orientation. Using a Z-step of 0.41 μm (the Z-resolution is approximately 1 μm), all larvae in all orientations have the same pixel size for all three axes.
All average image stacks are available for download . Raw image stacks are available upon request.
For image registration we used a voxel value-based protocol that sequentially applies rigid, affine and deformable transformations. The method extends the published protocol used for the registration of the genes expressed in the episphere of 48 hours post fertilization (hpf) larvae of Platynereis. The protocol presented here utilizes either staining of the larval axonal scaffold with acetylated tubulin, or DAPI-stained nuclei, as a reference marker for registration. The protocol was implemented using ImageJ  and the Insight Toolkit (ITK) .
Since larvae were scanned in random orientations, we developed two ImageJ algorithms for the rapid coarse orientation of whole-body scans. The scans were oriented based on the features provided by the acetylated-tubulin signal and DAPI stained nuclei. For 48 hpf larvae, we defined the AP orientation based on the prominent acetylated-tubulin signal of the prototroch ciliary band. We subsequently found the position of the ventral nerve cord in the acetylated-tubulin signal to define the DV axis. For 72 hpf larvae, we defined the AP axis by finding the body major axis based on the DAPI signal. For DV orientation we then used the position of the nerve cord. The algorithms were implemented in ImageJ using the Particle Analyzer, the TransformJ , and the Orientation plugins (Additional file 2). The rigid orientation of one image stack takes approximately 50 sec on a PC with 32 Gb RAM and an Intel® Core™ i5-2500 CPU 3.3 GHz x 4 processor.
The initial coarse orientation was refined using affine and deformable transformations (for detailed instructions and scripts [see Additional file 3) as implemented in ITK. The following are the components of the image registration procedure in ITK: two input images (template image and sample image), the optimizer, that searches for the optimal transformation for the registration (we used the ITK class RegularStepGradientDescentOptimizer for both affine and deformable image registration steps), the interpolator, that estimates the intensity of the pixels at non-grid positions after transformation (we used the ITK class LinearInterpolateImageFunction for both affine and deformable image registration steps), and the metric, that evaluates the alignment at each optimizer step (we used the ITK class MattesMutualInformationImageToImageMetric for both affine and deformable image registration steps). The Mattes mutual information metric  is a form of mutual information metric that is evaluated from a subset of pixels uniformly sampled from the image. The main difference to the other mutual information metric is that in the Mattes metric the subset of pixels sampled at the first iteration is reused in the subsequent optimizer iterations.
We used the multi-resolution approach (MultiResolutionImageRegistrationMethod class) for both affine and deformable image registration steps. This method first performs registration at low resolution and then stepwise registrations at higher resolution levels. This approach helps to avoid local minima. Fine registration takes approximately 50 min on a PC with 32 Gb RAM and an Intel® Core™ i5-2500 CPU 3.3 GHz × 4 processor.
To generate the templates, we used an iterative procedure similar to those described in [5, 11]. We selected 50 scans for both developmental stages and for both reference markers. The images were first registered to a manually selected high quality image. The registered images were sorted according to the metric (Mattes mutual information metric), and the top 40 scans were then used to generate an average image that served as the reference for the next iteration step.
After each iteration step we calculated the metric between each registered image and the average to which it was registered to determine the number of iterations necessary to generate a template. The value of the similarity metric reduces considerably after the second iteration, whereas the reduction after the third iteration is less substantial. Therefore we performed three iterations to generate a representative template for Platynereis.
To register gene expression patterns to a template, we first registered the corresponding reference marker (nuclear stain or acetylated tubulin) to the template and then applied the obtained transformations to the gene expression channel. We scanned up to ten samples for each gene in both developmental stages and averaged them to avoid bias towards one sample and compensate for natural and technical variation.
To evaluate the precision of the registration, we introduced a metric that provides information about the spatial precision of the registration, which cannot be determined based on the Mattes mutual information metric. First we determined the center positions of corresponding cells marked by the same gene expressed in single cells using ImageJ. We thresholded the signal in the image stacks and determined the X and Y coordinates in a Z-projection of the expressing cell. The Z coordinate was defined as the middle of the Z-span of the cell. We first calculated the average cell center coordinates for each group of corresponding cells from the individual registered samples and then calculated the distances of the individual cells to this average coordinate (absolute deviation).
To visualize registered gene expression patterns, we developed a ChannelMerger plugin for ImageJ that allows merging multiple channels in one RGB stack. The number of channels that can be merged is not restricted by the plugin. The plugin provides basic options such as color change and visibility setting. The merged image is compatible with ImageJ functionality (that is, saving as tif stack, adjusting brightness/contrast, viewing in 3DViewer, etcetera). The ChannelMerger plugin can be downloaded from .
The sequence of the TrpC cDNA for Platynereis has been submitted to GenBank [GenBank: JX916288].
In confocal microscopy, the absorption and scatter of light, as well as spherical and chromatic aberration, limit the attainable scanning depth. These aberrations are particularly damaging for signal intensity and resolution if the stimulus and emitted light has to pass through media with different refractive indices. The use of 97% TDE as a mounting medium with a refractive index of 1.515, closely matched to that of glass and oil, can greatly increase depth penetration when using an oil objective. TDE is miscible with water, easily penetrates across cell membranes, and is compatible with most fluorophores .
We tested TDE for use as a mounting medium for whole-mount in situ hybridization samples of Platynereis larvae. We counterstained the samples with the nuclear marker DAPI and an antibody against acetylated-tubulin labeled with a TRITC secondary antibody, to label neurites and cilia. We then performed three-channel confocal imaging of DAPI and TRITC fluorescence and the NBT/BCIP in situ hybridization signal. To scan the NBT/BCIP signal, we used a 635 nm excitation laser and a 780 nm long-pass filter, taking advantage of the far-red fluorescent emission of NBT/BCIP [21, 26]. We found that the far red signal is less prone to artifacts than reflection imaging, a method that takes advantage of reflection from the NBT/BCIP precipitate . Using reflection imaging, the cover glass, dirt particles, and the chaetae emit a strong background signal. Using far-red fluorescence, the glass and dirt particles have no signal, and only the spinning glands and the chaetae contribute to a background signal.
Spherical aberration also leads to the deterioration of optical resolution and image sharpness as a function of depth. To quantify the difference in image sharpness of glycerol- versus TDE-mounted samples, we used an edge detection filter (FeatureJ Edges implemented in Fiji ) on the DAPI channel (Figure 1D,H). We plotted the mean edge-filtered signal as a function of imaging depth. This analysis showed a more pronounced loss in signal quality with depth for glycerol- than for TDE-mounted specimens (Figure 1D,H and [see Additional file 4). TDE mounting thus allows the high-resolution confocal microscopic imaging of the entire volume of a Platynereis larva.
For image registration, we extended the protocol used in previous studies [11, 15]. Given the random orientation of our scans, we introduced an initial automatic rigid orientation step using a custom ImageJ script [see Additional file 2. This procedure performed the rough AP and DV orientation of the larvae based on landmarks in the DAPI and acetylated-tubulin channels. The oriented stacks were further processed using affine and deformable transformations as implemented in ITK [see Additional file 3.
Using this procedure we progressively aligned and averaged 40 individual scans from both 48 and 72 hpf larvae (Figure 3A,B,C,F,G,H). Each channel was aligned and averaged separately. Using an iterative procedure for template refinement we obtained high quality templates after three iterations (Figure 3B,C,G,H and [see Additional files 5, 6, 7 and 8]).
We next registered the expression patterns of test genes to the nuclear-stain and tubulin templates. We used five neuropeptide precursor genes (FMRFamide, SPY, L11, RGWamide and WLD) [15, 28], rhabdomeric opsin 1 (r-opsin-1), and a C-type transient receptor potential channel (TrpC). All of these genes show expression in a restricted number of neurons in the episphere and the trunk, allowing quantitative analyses [see Additional file 9.
The NBT/BCIP precipitate is known to block the path of light , therefore it was also important to test the method for a gene with a broader coherent expression domain. We chose prohormone convertase 2 (phc2), a gene with broad expression in the Platynereis nervous system ([see Additional file 10 and ). We scanned up to ten whole-mount in situ hybridization samples for each gene for both 48 and 72 hpf larvae and registered them independently to the nuclear-stain and the tubulin templates. We ran both the affine and the deformable registration steps until convergence of the metric or maximum 100 iterations [see Additional file 11. Visual inspection of the aligned image stacks and the deformation field showed that the individual samples showed good overlap with the template [see Additional file 11, even for the broadly expressed phc2 gene [see Additional file 10. We then averaged each gene’s expression and projected the averages onto the respective templates (Figure 3D,E,I,J and [see Additional files 121314 and 15).
The average absolute deviation for all cells was lower than the average cell body diameter (10.2 μm, s.d. = 2.3 μm, n = 100 along the X-axis, defined from the thresholded in situ signal). We also found that the nuclear-stain templates performed equally well or significantly better than the tubulin templates. We therefore used the nuclear-stain templates for further analyses. These results show that it is possible to register gene expression patterns in both the head and the trunk of Platynereis larvae with near-cellular resolution using the whole-body templates.
Taking advantage of such whole-body near-cellular resolution image registration, we next developed image analysis and statistical tools to quantitatively assess the overlap of gene expression pattern averages.
To facilitate the colocalization analysis of a large number of genes, we first developed a multichannel visualization plugin (ChannelMerger) for the open source image analysis platform ImageJ [20, 31]. The plugin allows merging any number of channels in one RGB stack, where each original channel is displayed in a different color [see Additional files 12 and 13. Alternatively, one can use commercial image analysis tools to display multiple channels (for example, Imaris) [see Additional files 14 and 15. Such multichannel views allow the visual inspection of a large number of genes, and the identification of potentially coexpressed genes.
To further enable fast and unbiased identification of colocalizing gene expression signals, we also developed a pair-wise channel merging macro [see Additional file 16 for the Fiji image processing package (Fiji is just ImageJ) [32, 33]. The macro first thresholds each average gene expression stack (Yen multilevel thresholding ), and displays the overlapping region together with the original unthresholded patterns on maximum projections [see Additional files 17 and 18.
The above methods rely on visual inspection and image thresholding to detect potentially coexpressing genes. Given the accuracy of the registration method, such approaches are reliable for broadly expressed genes that show considerable overlap. The analysis of genes expressed in single cells may require further statistical tests, or eventual experimental validation (for example, by double in situ hybridization).
One possibility to further test the potential coexpression of two genes showing single-cell expression in the same area in the atlas is to analyze their combined cell-coordinate statistics. If two genes coexpress in the same cell, the average absolute deviation of the combined cell body coordinates should be similar to that of the individual genes. If two genes are expressed in adjacent cells, then the average absolute deviation of the combined cell body coordinates should be larger than that of the individual genes.
Another approach is to calculate the absolute deviation of a ‘reference cell’ expressing one gene and compare this value to the average absolute deviation of a ‘test cell’ expressing another gene relative to the ‘reference cell’ average coordinate. If two genes coexpress (that is, the ‘reference’ and the ‘test’ cells are the same), then the two measures should be similar. If two genes are expressed in adjacent cells, then the average absolute deviation of the ‘test cell’ relative to the ‘reference cell’ average should be larger than the average absolute deviation of the ‘reference cell’.
Here we extended the gene expression registration protocol that was originally applicable to the episphere of Platynereis trochophores  to the whole-body of 48 and 72 hpf larvae. Cellular resolution whole-body registration will allow the development of a Platynereis gene expression atlas containing expression information for a large number of genes at different stages. Such an atlas will provide an overview of the spatial distribution of gene expression patterns, giving insights into the dynamics of gene interactions during development, and will provide a catalog of cell types, providing the foundation for understanding physiology.
A particularly important aspect of the Platynereis atlas is to identify genes that coexpress in particular cells, allowing the characterization of the molecular fingerprint and function of cells. Our quantitative analyses showed that the average absolute deviation for any registered cell is approximately half a cell diameter, meaning that gene expression patterns could be registered at or close to cellular resolution. The convincing demonstration of coexpression may nevertheless require experimental validation in some cases. For example, the atlas revealed a strong colocalization in the average signals of r opsin 1 and TrpC in the eye photoreceptors, which we subsequently verified experimentally. The coexpression of r opsin 1 and TrpC indicates that Platynereis rhabdomeric photoreceptors may also employ a phototransduction cascade similar to insect rhabdomeric photoreceptors  and intrinsically photosensitive retinal ganglion cells , strengthening the idea that these photoreceptors represent a homologous cell type .
Several different reference markers have been used for the generation of anatomical templates [2, 3, 8–11]. In Platynereis, the sharp and invariant acetylated-tubulin signal was used for the generation of the first gene expression map . Our results show that besides acetylated tubulin, nuclear stain can also be used as a reliable reference marker for Platynereis larvae. Nuclear stains (TOTO™-3 and Sytox green) have also been used in the recently established registration protocol for zebrafish larvae . Given that acetylated tubulin is not universally applicable (for example, it does not work in brachiopod larvae following in situ hybridization; A Hejnol, personal communication), nuclear stains may also be more suitable as reference markers in other animals. Nuclear stains have several other advantages over antibody markers. By labeling every cell, such stains give the broadest possible label of the anatomy. In contrast, tubulin is absent from large parts of the body, as in Platynereis larvae, and consequently provides less information for the registration. As our templates clearly show, nuclear stain gives a fuller representation of the anatomy, yet it is not a homogeneous label and reveals detailed internal structure (for example, is absent from neuropil and muscle fibers). Additionally, nuclear stain is cheap and easy to perform, and is compatible with standard in situ hybridization protocols (although a protocol with a sodium dodecyl sulfate (SDS)-containing hybridization solution yielded poor DAPI signal in Platynereis in our hands). We therefore propose nuclear stains as the first markers of choice for future image registration projects.
Here we introduced whole-body gene expression registration for Platynereis larvae, imaged by confocal microscopy. We hope that the reference templates and detailed instructions and scripts we provide will facilitate the building of a community resource for Platynereis. The possibility of near-cellular-resolution coexpression analysis makes Platynereis a powerful experimental model for the detailed characterization of cell types. We anticipate that by using optimized confocal imaging and whole-body scans it will be possible to develop reference templates for other developmental stages in Platynereis as well as other animals of a similar size range.
Differential interference contrast
Hours post fertilization
Iterative shape-averaging method
Nitroblue tetrazolium/5-bromo-4-chloro-3-indolyl phosphate
Sodium dodecyl sulfate
Virtual insect brain.
We thank Steffen Schmidt, Thomas Holder and Andre Noll for IT support. We thank Elizabeth Williams for her comments on the manuscript and Markus Conzelmann and Timea Keller for the initial characterization of TrpC. The research leading to these results received funding from the European Research Council under the European Union's Seventh Framework Programme (FP7/2007-2013)/European Research Council Grant Agreement 260821.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.