The Digital Bee Brain: Integrating and Managing Neurons in a Common 3D Reference System

The honeybee standard brain (HSB) serves as an interactive tool for relating morphologies of bee brain neurons and provides a reference system for functional and bibliographical properties (http://www.neurobiologie.fu-berlin.de/beebrain/). The ultimate goal is to document not only the morphological network properties of neurons collected from separate brains, but also to establish a graphical user interface for a neuron-related data base. Here, we review the current methods and protocols used to incorporate neuronal reconstructions into the HSB. Our registration protocol consists of two separate steps applied to imaging data from two-channel confocal microscopy scans: (1) The reconstruction of the neuron, facilitated by an automatic extraction of the neuron's skeleton based on threshold segmentation, and (2) the semi-automatic 3D segmentation of the neuropils and their registration with the HSB. The integration of neurons in the HSB is performed by applying the transformation computed in step (2) to the reconstructed neurons of step (1). The most critical issue of this protocol in terms of user interaction time – the segmentation process – is drastically improved by the use of a model-based segmentation process. Furthermore, the underlying statistical shape models (SSM) allow the visualization and analysis of characteristic variations in large sets of bee brain data. The anatomy of neural networks composed of multiple neurons that are registered into the HSB are visualized by depicting the 3D reconstructions together with semantic information with the objective to integrate data from multiple sources (electrophysiology, imaging, immunocytochemistry, molecular biology). Ultimately, this will allow the user to specify cell types and retrieve their morphologies along with physiological characterizations.

The honeybee standard brain (HSB) serves as an interactive tool for relating morphologies of bee brain neurons and provides a reference system for functional and bibliographical properties (http://www.neurobiologie.fu-berlin.de/beebrain/). The ultimate goal is to document not only the morphological network properties of neurons collected from separate brains, but also to establish a graphical user interface for a neuron-related data base. Here, we review the current methods and protocols used to incorporate neuronal reconstructions into the HSB. Our registration protocol consists of two separate steps applied to imaging data from two-channel confocal microscopy scans: (1) The reconstruction of the neuron, facilitated by an automatic extraction of the neuron's skeleton based on threshold segmentation, and (2) the semi-automatic 3D segmentation of the neuropils and their registration with the HSB. The integration of neurons in the HSB is performed by applying the transformation computed in step (2) to the reconstructed neurons of step (1). The most critical issue of this protocol in terms of user interaction timethe segmentation process -is drastically improved by the use of a model-based segmentation process. Furthermore, the underlying statistical shape models (SSM) allow the visualization and analysis of characteristic variations in large sets of bee brain data. The anatomy of neural networks composed of multiple neurons that are registered into the HSB are visualized by depicting the 3D reconstructions together with semantic information with the objective to integrate data from multiple sources (electrophysiology, imaging, immunocytochemistry, molecular biology). Ultimately, this will allow the user to specify cell types and retrieve their morphologies along with physiological characterizations.
The averaging method applied in the HSB is based on ideas of Ashburner (2000) and Guimond et al. (2000), who derive an average-shape image through an iteration of one affine registration, followed by multiple elastic registrations (Rohlfing et al., 2001(Rohlfing et al., , 2004. To create an initial average image, the method first registers all images to an (arbitrarily chosen) initial reference using affine registration. The method then registers all images non-rigidly to this average, generating a new average, and so forth. The underlying idea is that after several such iterations the average converges to the shape centroid of the population, which is, up to affine components (position, orientation, scaling, and shearing), independent of the choice of the initial reference image (Guimond et al., 2000, for more details on applying these methods for creating the honeybee standard brain: see Brandt et al., 2005).
The first step towards filling the atlas with structural information about neurons requires the semi-automatic segmentation of the neurons of interest and their related neuropils. The neuropil segmentations are spatially registered onto structures of the atlas. Then, the transformation coordinates produced by registration are used to fit neurons from different experiments into the atlas. Protocols for integrating genetically labeled neuron populations and single cell reconstructions into brain atlases have been described (e.g., Jenett et al., 2006;Kuß et al., 2007;Rybak et al., 2009). These protocols provide the electronic resources and tools for reconstructing neurons (Schmitt et al., 2004;Evers et al., 2005) as well as registration techniques that enable spatial normalization of structures using geometric warping algorithms (Rohlfing et al., 2001;Toga and Mazziotta, 2002;Westerhoff, 2003;Maye et al., 2006).
To facilitate the segmentation process, a statistical shape model (SSM) was developed. In order to develop the SSM a method was used that is based on successful procedures for automatic segmentation of medical imaging data (Lamecker et al., 2004;Kainmüller et al., 2007Kainmüller et al., , 2009Seim et al., 2008). The method integrates a priori information about variable neuropil shapes as imaged from confocal microscope imaging data (Neubert, 2007;Singer et al., 2008). Features in the imaging data are compared to intensity profiles of the confocal gray-value data, which have been learned by the model from a training data set. These comparisons are used to adapt all neuropil boundaries contained in the SSM to individual neuropil boundaries in the present imaging data.
In this study, emphasis was placed on the strategy of standardizing and optimizing the registration process and the subsequent fitting procedures of neuronal data. These issues are important to make the HSB usable for researchers at other labs. The applicability of statistical shape atlases that contain information about brain structures and their variability are discussed.
Furthermore, we introduce an ontology-based approach integrates vast amounts of data from various experimental sources in a structured way into a single coherent database.
The HSB atlas and examples of registered neurons can be downloaded and visualized at http://www.neurobiologie.fu-berlin.de/ beebrain/

Materials and Methods
All animals (workerbee foragers, Apis mellifera carnica) were taken from the hives at the Institute for Neurobiology, Free University, Berlin, Germany.

Lucifer yellow histology
For preparations that were utilized for construction of the SSM Lucifer yellow was used as a neuropil background stain. Brains were fixed in 4% para-formaldehyde (pFA) (Sigma) for either 2 h at room temperature or overnight at 4ºC. After washing, brains were dehydrated in ascending ethanol series, and cleared in methylsalicylate. 4% Lucifer yellow was added either to the fixative at a dilution of 1:500 or to PBS-TritonX treated overnight.
Ethyl gallate histology (Wigglesworth, 1957) Brains were fixed for 4 h in 2.5% glutaraldehyde in cacodylate buffer. After several washes in buffer, brains were osmicated in 2% OsO 4 in cacodylate buffer for 1 h in the dark. Tissue was then transferred to 0.5% ethyl gallate (Merck) in distilled water for 1-4 h. The solution was changed until the blue-gray color disappeared. After thorough washing in distilled water, the specimens were dehydrated and embedded in Durcupan (Fluca). Ethyl gallate preparations were sectioned at 10-25 μm.

Confocal microscopy
Whole-mount brains used for creating the SSM were counterstained with Lucifer yellow and imaged sequentially with the Leica TCS-SP2 confocal microscope using a 10× dry or 10× oil Leica objective (HC PL APO 10 × /0.4, Leica, Bensheim, Germany). For Lucifer yellow-stained tissue, the Ar-Kr 488-nm laser line was used at a voxel resolution of approximately 1.5 × 1.5 × 3 μm. The dyefilled neurons were excited using the 543-nm line (detected with an emission spectrum of 550-620 nm) or the 633-nm line (detected with an emission spectrum of 650-750 nm) of the HeNe laser. For high-resolution scans of intracellularly stained neurons, we used the 20× oil (HC PL APO 20 × /0.70) 40× oil HCX PL APO CS 40.0×/1.25) and 63× oil objectives (HCX PL APO 63×/1.32-0.60) (Leica, Bensheim, Germany). Depending on the zoom factor (1-4) the voxel resolution was approximately 0.1-0.4 × 0.1-0.4 × 1 μm. In all confocal scans we used a pixel resolution of 1024 × 1024 in xy axes and an 8 bit intensity resolution. Because of the refractive index mismatch in the optical path, dry lenses usually introduce a shortening of distances in the z axis. According to Bucher et al. (2000), shortening can be considered as a linear scaling in the z direction. Therefore, the scaling factor from preparations that were scanned with dry lenses is estimated to be 1.6. (refraction index: methylsalicylate = 1.51, oil = 1.54, air = 1).
A description of the registration process and parameter settings can be found in the Amira User Guide (Visage Imaging, Berlin; San Diego, CA, USA). A detailed protocol is found in the section "Registration protocol" of supplementary material (see also Kuß et al., 2007).

Registration using landmarks
Histological sections stained with ethyl gallate were registered into the HSB using the LandmarkWarp module of Amira. The corresponding anatomical locations in the histological sections and the gray-value dataset of the HSB were defined by two sets of landmarks.

the statIstIcal shape Model
A SSM captures the mean shape and the geometric variability of a given set of input geometries by a limited number of parameters and can therefore be used for automated and robust image segmentation. The SSM of the bee brain applied in this work was presented in Lienhard (2008). A detailed description of the procedure and references to other approaches can be found in Lamecker (2008). Our strategy for the creation of an SSM from three-dimensional image data includes several steps: (1) Triangular surfaces were reconstructed from 16 manually labeled image stacks of central neuropils of the worker bee. In contrast to the previously presented HSB, the neuropil counter stain was achieved using 5% Lucifer yellow. Most labels for neuropils were chosen as defined in the HSB (Brandt et al., 2005). The median and lateral calyces of each mushroom body were segmented without subdivisions. The optic and antennal lobes were not included in this model. A subdivision of the calyces into lip, collar, and basal ring was not performed.
(2) Point-to-point correspondences were established between all training surfaces, i.e., vertices with the same index share the same anatomical position on each training surface. To achieve this, all surfaces were equally partitioned into regions (patches), which represent homologous biological compartments. A reference surface triangulation was then mapped onto each training shape using surface parameterization techniques (Lamecker et al., 2003). (3) After alignment of all training shapes to one reference, principal component analysis was applied on the set of shape vectors representing the surfaces' vertex coordinates. This process generated the SSM representing the average shape plus a linear combination of the most characteristic modes of shape variation (shape modes) contained in the training set.
The SSM allows a highly compact representation of shape variations among a large number of individuals. The most challenging step in the SSM generation is the identification of corresponding points (step 2). Our approach is interactive as it involves manual specification of patch boundaries. Yet, this allows the production of accurate SSMs even for very complex geometries with large deformations and arbitrary topologies, as in the case of the bee brain.

SSM-based automated image segmentation
Image segmentation can be automated by using a priori knowledge, particularly about geometrical shapes and intensity profiles. The general idea is to roughly position an SSM in the imaging data and subsequently vary the shape parameters (weights of the shape modes) and the spatial location until the SSM matches the object in the imaging data as closely as possible (Lamecker, 2008). The intensity distributions of the underlying imaging data are evaluated around tracIng of neurons All confocal scans were digitized as double channels after which each channel was analyzed separately using the three-dimensional visualization and segmentation modules in Amira (version 4.1; Visage Imaging, Berlin; San Diego, CA, USA). Tracing and reconstruction of the neurons, including topology, lengths, and diameters, were performed using a module that integrates methods presented in (Schmitt et al., 2004;Evers et al., 2005). Traced single neurons, which were Amira data SkeletonTree format, were converted to the LineSet format and then triangulated to meshed surfaces. These surface files were exported as wavefront (obj) files. Wavefront files of neurons and neuropil surfaces were imported with the Adobe 3D Reviewer to Adobe Acrobat Pro Extended (Adobe Systems, Inc.). The images in the PDF version of this manuscript can be viewed by using the 3D viewer mode of the Acrobat Reader (version 9 and higher, which is freely available at http://get.adobe.com/de/reader/.

Semi-automatic segmentation
Image segmentation was performed semi-automatically using Amira 4.1. The segmentation results were image stacks of type LabelField. LabelFields assign a label, which represents a distinct (brain) structure, to each voxel.
In most cases, no image preprocessing on the raw images was necessary, except for an adjustment of the gray-scale window. In some cases, Gaussian smoothing and unsharp masking from Amira's DigitalFilters tool were applied to enhance faint contours. Image stacks were then loaded into the interactive segmentation editor. Using the segmentation editor's BrushTool, the neuropil areas of interest were traced manually, slice by slice. To facilitate and speed up this manual segmentation process, a method that automatically interpolates segmentations between image slices was applied. In a post-processing step, connected areas of voxels containing only a small number of voxels were eliminated with the RemoveIslands tool. Finally, the LabelFields were smoothed (SmoothLabels).
An image series showing the brain neuropils of the HSB with superimposed LabelFields using the ColorWash module is provided in the Movie 1 in supplementary material. The final segmentations were supervised by a segmentation expert (curator of the HSB, JR). Depending on the staining quality of the tissue, an experienced user needs a minimum of 8-10 h to reconstruct all neuropils defined for the HSB.

Registration and transformation of neurons into the HSB
To fit a neuron into the HSB, two steps are applied. First, the neuron's related neuropils are registered to corresponding parts of the HSB. This requires an affine and a subsequent elastic registration, which respectively result in a 9-degree of freedom transformation matrix and a deformation field (VectorField). Second, the registration resulted are used to transform the neuron's geometric representation (SkeletonTree or LineSet).
The affine and the elastic registration procedures use a metric that takes the spatial correspondence of two label fields into account. Similarity measures were used as label consistency. The affine registration procedure further uses a hierarchical optimization algorithm going from coarser to finer resolution. such as Lucifer yellow. Lucifer yellow staining, though reduced in contrast compared to the synapsin antibody staining, allows more rapid and easier histological processing. It further detects neuronal structures in much more detail than autofluorescence, and results in homogenous staining throughout the brain.
The challenge for transforming neurons into the common frame of the HSB (Figure 1) is to segment the complete Gestalt of each neuron at low resolution together with its spatial relations to adjacent structures (often gained by high-resolution microscopy) such that a warping algorithm allows for precise matching between samples from different individual brains.
In Figure 2 several examples of neurons are shown, each stemming from a different preparation, and each having been separately transformed into the HSB. To capture all parts of a neuron at different resolutions, several confocal scans were necessary. This poses the problem of registering several label fields leading to subsequent errors during realignment of all parts of the neuron that accumulate during the registration processes, which may induce artifacts at the edges of the label fields (Maye et al., 2006). The strategy used for each of the examples of Figure 2 was to first scan an overview of the brain containing the whole neuron, and then reconstruct the main parts of the neuron. Afterwards, the neuron's regions of interest, for example, the dendritic tree of the L3 neuron, were reconstructed from highresolution scans ( Figure 2B), downscaled, aligned, and merged to the main neuron reconstructed from low resolution scans. Deviations that are caused by using several independent regions of the same neuron in several steps of registration are kept to a minimum. As a result, the spatial accuracy of the warping process is enhanced. This iterative procedure allows us to compose compounds of registrations at different levels of resolution, for example, the target areas of olfactory and mechanosensory interneurons (Ai et al., 2009) in subregions of the protocerebral lobe (arrows in Figure 2A) together with registrations at high resolution of their fine structures for the analysis of their local topological features (Figures 2C,D). The movie the current SSM, and the surface of the SSM is displaced according to a set of given rules leading to the displacement model, which will be discussed in the next section. From the computed displacement, new shape weights or new locations are computed. This procedure guarantees that the segmentation indeed represents a plausible shape (robustness). In order to overcome possible mismatches due to individual variations not captured by the SSM, a post-processing step usually provides some fine tuning for accuracy (Kainmüller et al., 2007(Kainmüller et al., , 2009Seim et al., 2008). Apart from the SSM itself, the main ingredient is a rule for displacing the SSM according to the underlying imaging data to be segmented. This is provided by the displacement model.

The displacement model
A simple method for computing displacements of the surface model in the imaging data is to determine a normal displacement for each vertex of the model such that the new vertex position coincides with a strong gradient in the imaging data. Here, the only assumption made is that object boundaries in the imaging data are reflected by significant local variations in the image intensity. In the case of confocal imaging of bee brains more information about the imaging process can be included in order to refine this model. Such extensions have been proposed by Neubert (2007) and Singer (2008) for different imaging protocols.

Segmentation performance
Accuracy. Cross-validation tests are used to estimate how reliably the SSM-based segmentation of the bee brain performs in practice.
In the leave-one-out test one image is removed from the training set, and a calculation is performed to determine how accurately a reduced SSM, which is constructed from the remaining images, can be adapted to that removed image. The ability of the model to describe arbitrary shapes is described as completeness or generality. In order to estimate the quality of the displacement models (intensity profile analysis) leave-all-in tests were performed. In contrast to the leave-one-out test, the known image is not removed from the SSM. This way the performance of the displacement strategy can be measured, separated from the quality of the SSM itself.

Performance measures.
Measures for the mean and maximal surface distance were calculated using the Amira SurfaceDistance module. This module computes several different distance measures between two surfaces. The following measures were computed from the histogram of these values: mean distance and standard deviation; root mean square distance; maximum distance (Hausdorff distance); medial distance; area deviation (percentage of area that deviates more than a given threshold).

ConfoCal miCRosCopy
Because the registration process is based on label fields, the neuropils have to be stained in such a way that neuropil borders can be identified and segmented. We first applied an antibody against synapsin which nicely stains neuropils dense in synapses, and thus contrasts its border to surrounding tissue (Brandt et al., 2005). A simpler method involves imaging the autofluorescence of the tissue induced by glutaraldehyde. One could also enhance the autofluorescence using dyes  of supplementary material S3 and the interactive viewing mode in the PDF file (Figure 2A) allows visualization of the neurons' spatial relationships.
Estimating the accuracy of the registration process is a difficult task. Certainly, the process involves numerous steps, which are prone to induce inaccuracies (see above). For example, the histological procedure may induce variable distortions due to local shrinkage differences that are not fully compensated by the affine and elastic registration. The experimenter may not segment correctly; the number of label fields (and, thus, the number of registrations steps) and thus distortions of the reconstructed surfaces of neuropils and neurons will lead to incorrect locations of the transformed neuron in the HSB. In addition, the neuropils and neurons themselves will differ from animal to animal, and it is this variability that determines the fundamental limit regarding the reliability of any brain atlas -besides the methodological problems. Therefore, it is not possible to derive a measure of reliability in the composition of neurons registered sequentially into the brain atlas. The best way of checking the spatial accuracy of such neurons is by comparing the relative positions of the neuron to the neuropil border lines in the original preparation to the situation after the registration process. We provide circumstantial evidence for the reliability of the segmentation and registration process by describing an example in which we compared the locations of intracellularly stained neurons with cross sections in high-resolution ethyl gallate-stained paraffin sections using two different registration methods.

IntegratIng data collected by dIfferent hIstologIcal methods
In Figure 3 a registered olfactory projection neuron (L5) is visualized together with a registered ethyl gallate section showing midbrain regions (mushroom bodies, central body, and protocerebral lobe). The ethyl gallate method (Wigglesworth, 1957) provides detailed information about the neural architecture revealing the composition of neuropils, somata, and tracts, thus capturing the spatial context information. A whole series of ethyl gallate-stained sections was first used to identify the median and lateral antennocerebralis tracts (ACTs) in a correlative light and electron microscopy study (Rybak, 1994). Our future goal is to integrate data from these histological procedures into the HSB. Here, we demonstrate the spatial accuracy of the registration process. A horizontal ethyl gallate section was warped into the HSB using a landmark-based registration by finding corresponding points or landmarks in the HSB and the histology section (Figure 3). Separately, a single stained L3 axon, which was transformed to the HSB using a July 2010 | Volume 4 | Article 30 | 6 Rybak et al.
Honeybee standard brain To achieve the accuracy required for neuronal connectivity estimates, the development of higher resolution atlases that define subregions of the brain and allow their integration into a common coordinate system is needed. An example is given in this special issue of Frontiers in Neuroscience by el Jundi et al. (2010) who constructed a high-resolution atlas of the central complex of Locusta and determined the spatial relations of two central complex neurons.

Advantage of the SSM
The accurate and reliable localization of region boundaries during the registration process is a prerequisite for fitting neurons into any Standard Atlas. This is true for manual (or semi-automatic) and fully automatic segmentation techniques, though the advantage of the latter is that the level of human expert interaction is reduced. Labeling performed by different individuals often leads to variable results. A model-based auto-segmentation of neuropil boundaries utilizes a priori knowledge about the 3D shape of an object, in our case the bee brain, and characteristic features of the imaging data. Such a model provides a measure of the variability of the object set and can therefore be used to analyze and quantify morphological volumetric changes in neuropiles of the adult animal.
An SSM of the central bee brain (excluding the antennal and optic lobes) was calculated from 16 training shapes, resulting in 17 shape modes. These were extended to 32 shapes (31 shape modes) by mirroring the right and left brain hemisphere along the neuraxis in each preparation. Each training shape was manually segmented by labeling the neuropil boundaries of the confocal imaging data stained with Lucifer yellow. Triangulated polygonal surfaces were reconstructed from the labeled images and simplified to 150000 triangles. The surfaces were then affinely aligned using label field registration, runs through the corresponding ascending and descending parts of the median ACT (m-ACT) as seen in the ethyl gallate section (Figures 3A-C). The spatial accuracy of the registration process is indeed very high, and allows identification of m-ACT neurons in the median and lateral antenno-cerebralis tract (l-ACT) even at the single neuron level (see black arrows in Figure 3B).

analyzIng putatIve synaptIc connectIons
Fitting neurons into the HSB can be achieved with a certain degree of accuracy with regard to spatial relationships, but thus far cannot replace studies on synaptic connectivity. This must be achieved by electron microscopy (e.g., Ganeshina and Menzel, 2001) or by some approximation in confocal co-localization studies on the light microscopy level. Combining high-resolution confocal laser scanning microscopy with precise three-dimensional dendritic surface reconstruction (Schmitt et al., 2004) allows for automated co-localization analysis in order to map the distribution of potential synaptic contacts onto dendritic trees or axon terminals (Evers et al., 2005;Meseke et al., 2009).
This technique was used to estimate the distribution of putative GABAergic synaptic contacts on the dendrites of the Pe1 neuron, a single identified mushroom body extrinsic neuron (Mauelshagen, 1993;Rybak and Menzel, 1998 , Figure 4). GABA-like immunoreactivity has been shown for the A3 feedback neurons (Schäfer and Bicker, 1986;Grünewald, 1999). These neurons also innervate the mushroom body lobes and peduncle and may provide local inhibitory input to the Pe1 neuron ( Figure 4B, green label). In Figure 4D the distribution of putative inhibitory input synapses onto the Pe1 dendritic tree is highlighted by red dots, indicating a distance of GABA-like immunoreactivity profiles of up to about 300 nm (Okada et al., 2007). The SSM-based segmentation algorithm, as described in the section "Methods", was applied until no further improvement could be achieved, meaning that the change of the shape between consecutive iterations fell under a defined threshold. The displacement vectors were computed via analyzing 1D profiles of image intensities along surface normals at each vertex ( Figure 5D). Based on these displacements the SSM was iteratively adapted (see Figures 5C,D and the movie in section "Automatic segmentation" of supplementary material). The process of initial positioning and the adjustment of the model in an exemplary training dataset image (LY12) are shown in Figure 6.
Application of the segmentation algorithm on the training set images results in an average surface distance of 4.06 ± 0.95 μm of the fitted model to the manual segmented shapes. A leave-one-out test simulating new image data yields to a distance of 8.82 ± 1.02 μm. In comparison, two manual segmentations of the same imaging data lead to a mean surface distance of 4.79 μm.
The SSM and displacement algorithm was also tested with confocal data not contained in the training set that was used for creation of the shape model. We used extracted parts of the SSM in combination with the displacement model in order to automatically segment the mushroom body neuropil in high-resolution confocal scans ( Figure 7A). Figure 7B shows the results for a single section using either the whole mushroom body calyces (red intersects) or the calyces (yellow intersects) to auto-segment the structures. Slightly better results were achieved using the reduced model of the calyx (arrows in Figure 7B). Figures 7C,D shows manually segmented neuropil borders in comparison to the automatic segmentation. Figure 8 provides a direct comparison of the HSB and SSM. The mean surface distance between HSB and SSM amounts to 8.5 μm. Particularly large distances were found at the median calyx (MC) and subesophageal ganglion (SOG), possibly indicating stronger shrinkage dependencies induced by the different histological procedures. It takes an experienced segmentor around 8 h to manually segment those central brain structures used for the SSM.
In contrast, a further evaluation of the quality of automatic segmented bee brains and the estimated post-processing times amounts currently to approximately 3-4 h by an experienced segmenter.

HierarcHical structure labeling and browsing
The HSB created so far contains only geometric and topological information about neurons and neuropils. For many applications, semantic information also needs to be integrated. This semantic information includes information about the hierarchical organization of brain structures and information about relations between structures. An example would be the description of the anatomical proximity of neuron A and neuron B and the possible communication between them. Often neurobiologists have this knowledge, but the information needs to be made explicit.
In recent years, the development of ontologies has been an appropriate choice for capturing and representing semantic information in many fields, including biology. In the information sciences, an ontology is a formal representation of concepts or structures and relations among those structures in a defined application domain. Visually, ontologies can be described as graphs in which structures are represented by nodes and relationships between structures by edges. In ontology modeling, there are two different types of nodes: the geometrical center as a reference point, and then transferred to a common coordinate system. In order to map the surfaces of the training shape in a proper way, certain conditions are required for creating the surfaces (for details see Lienhard, 2008). In order to achieve a correct model of a biological structure, one needs to find the corresponding anatomical points on all training shapes. To determine correspondences, surfaces were divided into 89 regions according to shape features and anatomical landmarks ( Figure 5A). A principal component analysis provided a linear model of the shape variability of the training set ( Figure 5B, see also the Movie in section "The statistical shape model" of supplementary material).
In order to place the SSM of the central bee brain roughly into the confocal images an affine registration was used applying a nondeformable model that contained only brain regions that represent borders to exterior structures. A positioning algorithm recomputed rotation, scaling, and translation parameters using characteristic image features of brain tissues. In order to reduce noise in the lowcontrast Lucifer yellow stain a non-linear isotropic filter was applied to the data (Weickert, 1997;Lamecker et al., 2004). classes and instances. Classes describe common concepts, such as brain. An instance of such a class could be, for example, the brain of animal A. Relations appear among classes, among instances, and among classes and instances. Relations can have different types, as for example the type isA. With these basic tools we can create the statement brain of animal A isA brain.
To support the understanding and analysis of structural and functional characteristics of brain structures, ontologies need to fulfill several requirements. An ideal brain ontology would include a complete set of structural parts and neuron types. It would also contain axonal projections between regions and neuron types, and it would include morphological, connectional, and functional properties of these particular neurons. According to Bota and Swanson (2008), an ideal ontology would be species specific. These authors also state that the development of such an ontology is a long-term goal for a community project. Indeed, we consider our attempts as an early step only.
Our ontology uses predefined classes of the foundational model of anatomy (FMA) (Rosse and Mejino, 2003). The most important of these are Cell, Cell_Part, Organ, and Organ_Part where Cell and

dIscussIon
The motivation for creating a digital atlas of the bee brain originated from the experience that a large amount of information is lost when single, intracellularly marked neurons are drawn on paper traced in camera lucida projections or just photographed. Neurons are three-dimensional entities embedded in a network of other neurons, and it is this information that is required in the future to interpret functional properties of neurons and neural networks in relation to their structure and connectivity (Abel et al., 2001;Müller et al., 2002;Krofczik et al., 2008). Insect brains are small enough to be scanned fully with confocal microscopy at a reasonable resolution. Therefore, no border problem appears, or at least it is reduced to the spatially limited connections with the ventral chord. Furthermore, many neurons in the insect brain are individually identifiable, and quite a number of them have already been identified (e.g. Hammer, 1993;Mauelshagen, 1993;Menzel, 2001;Heinze and Homberg, 2008;Homberg, 2008; see also this issue). Often neural tracts or compositions of local neurons consist of a few hundred neurons allowing for the possibility that in the not too distant future all neurons of a particular neuropil or part of the brain will be described in their morphology. In that case one would need this "description" in a digital 3D format so that the full power of mass data computation can be applied to visualize zoom Cell_Part only consider neurons and Organ and Organ_Part only consider neuropils of the bee brain. We further restrict our relations to describe spatial, morphological and, if available, connectional properties. Table 1 lists the most important relations used in our bee brain ontology. Figure 9 shows a scheme of how classes, instances and relations are connected using the Pe1 neuron as an example (see also Figures 2 and 4). The editor Protégé was used to create the ontology. Currently, the ontology contains 100 classes attached to 600 instances, and 1300 relations of 17 types. Integrated are several neuron types including the location representation of their somata, axons, and dendrites. This ontology has been linked to the reconstructions of the HSB by assigning the reconstructions' ID and file name to appropriate instances of the ontology. This step enables ontology-based browsing of the atlas (Kuß et al., , 2009. In a first usage approach of the ontology-linked HSB, we addressed the automatic creation of meaningful visualizations. Good visualizations transport a large amount of information and form an important communication medium. This information can be used to present research results, to communicate with research partners or to teach neurobiology. Often the process of creating such meaningful and expressive visualizations is time-consuming and requires sound knowledge of the visualization software used. In our approach, the user only selects a structure to be visualized and a predefined query, such as "Show overview". Then, an algorithm automatically creates a visualization that contains the selected structure highlighted as a focus object and further structures forming the context. This works as follows: Each predefined query owns a set of relation types considered to be relevant. Starting at the from creating the brain atlas and filling it with useful information is limited to understanding the spatial relationship. But this will soon change, because the information stored in the framework of the atlas will allow us to pose new questions, to discover novel patterns of neural connections, to assemble and organize large amounts of information, and to relate function to structure as proposed by us as well as other authors (e.g. Namiki and Kanzaki, 2008;Staudacher et al., 2009).
We began our project by creating an average atlas from 20 bee brains whose 22 neuropils were segmented manually and then used for the averaging process. The composition of these neuropils made it possible to calculate rigid and elastic transformations that provided enough information for faithful registration of neurons. The average-shape property ensures that the deformation applied determine potential connectivity patterns, to derive quantitative measures of distances, diameters, branching patterns, potential synaptic sites, and relate structures to functional components such as distribution of transmitters, receptors, channels, and intracellular molecules. At the moment we are far from getting even close to these goals, but important groundwork has been done, and one can hope that the time-consuming (and tedious) steps towards reaching these goals, such as manual segmentation, correction of errors, complicated procedures during registration, neuron tracing, (see Maye et al., 2006) will soon be overcome or become less cumbersome (for an review on automated registration and neuron tracing methods, see Peng, 2008). Apart from the aesthetic pleasure one experiences in visualizing single neurons and their compositions within the 3D atlas, right now the reward gained Figure 7 | evaluation of the auto-segmentation process using a reduced shape model for confocal imaging data. (A) Confocal Lucifer yellow image of the median and left brain hemisphere scanned with a 20× oil objective. (B) Result of the automatic segmentation process. Arrows indicate refinement of automatic segmentation using isolated neuropils: yellow lines when only the calyx was used, red line when the whole mushroom bodies were used. (C, D) Surface representations of shape differences between manually (transparent) and automatically segmented (solid) labels. The mean surface distance of the manually segmented and the automatically segmented surface is 11 μm using the mushroom body model of the SSM as compared to 10.5 μm using the calyx model. Scale bar in (C): surface distance in μm.
inaccuracies are due to surface reconstruction (smoothing) from segmented label fields, and to cropping of areas of interest. For the latter the edges of cut regions are not well defined for the individual brain and the standard reference and are therefore difficult to operate for the registration algorithm (Maye et al., 2006).
The question of whether two close neurons are potentially in synaptic contact is much more difficult to answer and may well be beyond the scope of the registration process. Double markings in the same brain combined with electron microscopy will also be necessary in the future to prove such contacts, but the registration process already provides conclusive evidence that will either provide the motivation to start such a demanding project or forsake it altogether (e.g., Hohensee et al., 2008;Meseke et al., 2009). We provided one example to document that registered neurons are so precisely embedded in the histology of high-resolution light microscopy cross sections (Figure 3) that one may well conclude that their location relative to surrounding neurons is down to the precision of a few microns. In a similar approach single-cell labeling was used in combination with non-rigid registration techniques to estimate the synaptic density and spatial relationship of Drosophila olfactory interneurons (projection neurons) in central brain areas (Jefferis, 2007). Using intensity-based image registration for averaging the brains, they estimated the accuracy of registration up to a few microns. Greater accuracy might not be possible because the protocols always depend on neurons from different brains, and the variance of neuron Gestalt from brain to brain will limit the resolution.
A general problem with neuron reconstructions relates to the fact that high-resolution imaging microscopy requires lenses whose working distance is often not large enough to cover the whole dendritic tree of the respective neuron, whereas low resolution images are necessary to connect parts of the same neuron or different neurons. Physical sections, e.g., vibrotome offer a solution to this problem, but consecutive sections need to be aligned such that the neuron can be fully reconstructed and registered into the atlas (el Jundi et al., 2009(el Jundi et al., , 2010. So far we have registered 50 neurons into the HSB, a small proportion indeed of the approximately 950000 neurons of the bee brain (Witthöft, 1967). However, even this small number calls for more sophisticated means of visualization, selection of combinations of neurons and ways of highlighting particular properties of the network arising from these neurons (potential contacts, estimated information flow, combination with data from, for example, immunocytochemistry and electron microscopy, electrophysiology, and Ca-imaging).

the statIstIcal shape atlas
An atlas derived from an averaging process (as the HSB, Figure 1) contains a large amount of information about spatial relations of structures. It is highly suggestive to use this information for one of the most time-consuming, difficult and tedious steps, the segmentation process of the structures (neuropils) necessary for the registration of any individual brain. We took up this argument and implemented a procedure, a model-based auto-segmentation, originally developed for the analysis of shape variability and modeling of structures in medical imaging (Lamecker, 2008). This method was adapted and applied to to the individuals remains small. Furthermore, manual segmentation as applied in our first approach is subject to noise, i.e., contouring between slices varies according to criterion variability of the experimenter. Averaging several such noisy label images reduces random parts of the contours, thus increasing the reliability of the standard. We have shown (Rohlfing et al., 2001;Brandt et al., 2005) that the non-rigid registration is able to increase the distinctness of inner structures such as tracts and strata even though the algorithm does not "know" about those structures, because it is applied on the label images without interior structures. We deduced from our observations that registration fidelity is sufficient for the spatial scale level of the standard brain. This result also makes us optimistic that a non-rigid registration of neuropil boundaries to the standard yields a reliable and reasonably accurate estimate of the "true" position of a co-stained neuron within the standard.
As pointed out above, it is not easy to evaluate how accurate the registrations of neurons are using the average neuropil borders as guiding posts. When a neuron runs close to the border of a neuropil used for registration a small deviation from its relative position becomes very important, that is the neuron lies either inside or outside the particular neuropil. We have observed these inaccuracies, and they can be corrected by repeating the segmentation and registration processes (see section "Registration protocol" of supplementary material). Preparation artifacts due to dissecting of the specimen and histological processing can lead to distortion effects that are only partly corrected by the registration algorithm. Further Figure 8 | A comparison of the iterative average shape brain (HSB) and the statistical shape model (SSM). The mean surface distance measured after rigid registration of HSB onto the SSM amounts to 8.5 μm. Only the surface model of the HSB is shown here, and the surface distances between HSB and SSM are indicated by colored vectors (arrow and false-color scale bar). Note the particularly large distance at the median calyx (MC) and subesophageal ganglion (SOG), which might be due to the stronger shrinkage process caused by the different histological procedures employed for the two models (see section "Methods" and text). the central bee brain to generate a SSM (Figure 5) that reflects the shape variability of a restricted set of 16 bee brains. In combination with a displacement model (Kainmüller et al., 2007(Kainmüller et al., , 2009Lamecker, 2008), based on evaluation of intensity gradients within the confocal images, the SSM allows automatic neuropil segmentation (Figures 5-7). Our approach was combined with alterations in the histological procedure. Lucifer yellow was used instead of synaptic antibodies as a neuropil background stain. Furthermore, only a selected part of the brain was used (the central brain excluding the antennal lobes). Lucifer yellow treatment during the histological procedure provides us with the same information as the antibody, but is faster and the results in more homogenously stained neuropils. The focus on the central brain allowed us to test the power of the approach for the most important and most complex structures of the bee brain, the mushroom bodies (Figures 4 and 7). An automatic procedure for segmenting neuropils provides the following additional advantages: (1). A priori knowledge about the 3D shape of the objects in question, in our case, the bee brain with its characteristic features of the image data. (2). Measures about the variability of the object set, and thus provides us with information which can be used to analyze and quantify any changes induced during development or on the basis of different genetic backgrounds (as described by Kurylas et al.,, 2008;el Jundi et al., 2010 for the VIB). (3). Information useful for across-species investigations in order to analyze evolutionary changes by comparative analysis of brain structures at least in closely related groups (e.g., hymenoptera, Gronenberg, 2001;beewolf, Rybak et al., 2003).
We found that SSM can well be used to detect neuropil borders in these preparations. Nevertheless there are deviations in the automatic segmentation process (see Figure 7). An analysis of the post-processing time by a segmentation expert shows that one still saves considerable time and gains accuracy (Rybak, personal observation).
A quantitative comparison of HSB and SSM (Figure 8) reveals that the volume of the HSB is smaller relative to the SSM. Stronger shrinkage of the HSB might be due to the prolonged incubation time required for the synaptic antibody procedure (1 week). Nevertheless, since shape differences in segmented brains by either the synaptic antibodies or the Lucifer yellow method seem to be small, it is reasonable to take segmented data used for the HSB and include it into the enlarged model of the SSM. Moreover, all imaging data used for the HSB were segmented by experts, thus providing very reliable definitions of the neuropil boundaries, and the quality of the SSM will be enhanced by adding the HSB dataset, since shape variability represented by the enlarged SSM will capture a higher number of histological procedures used for insect preparations (i.e., fixation, incubation of antibodies, use of fluorescent dyes).
One disadvantage of the current HSB is that it is closed to improvements and adaptations, which will certainly result from more appropriate histological procedures such as the shift from antibody staining of neuropil borders to easy-to-use fluorescent dyes. Working with the SSM will allow us to create a novel form of HSB that grows with each brain, and which adapts stepwise any morphological changes with histological procedures. Thus, enlarging the set of brains included in the SSM by already segmented brains used for the HSB and by developing more elaborated displacement algorithms based on intensity profile analysis will allow us to create such a new atlas based on many more brains, particularly in the context of different experiments (electrophysiology, immunocytochemistry, etc). Additionally, combining a shape atlas with proper registration techniques will allow us to use such an atlas, initially created to replace manual segmentation (for the label field registration of neurons). Once the deformation field that fits the SSM model to a confocal image is calculated it can be used to integrate the neuron into the average-shape atlas. Such an approach can forego and eventually replace the label field registration as in the current HSB. http://life.bio.sunysb.edu/morph/glossary/gloss1.html. http://life.bio.sunysb.edu/morph/glossary/gloss2.html.

Supplement 1: ColorwaSh hSB movies in wmv format
These movies show confocal images of the honeybee brain in frontal, horizontal and sagittal directions. Neuropils were stained with synaptic antibodies. The colored label fields of brain neuropils as defined for the Honey Bee Standard Brain (HSB, Brandt et al., 2005) are superimposed. Note that some neuropil areas that are not defined in the HSB are labeled. Abbreviations: AL: antennal lobe, a: alpha-lobe, b: beta-lobe, CB: central body, li: lip, co: collar, br: basal ring, Me: medulla, Lo: lobula, PL: protocerebral lobe, SOG: subesophageal ganglion. MC: median calyx, LC: lateral calyx. ot: optic tubercle,, pb: protocerebral bridge, DL: dorsal lobe, lac: lateral accessory lobe, pe: peduncle, lh: lateral horn.

Supplement 2 pdf
Registration Protocol that describes the incorporation of neuronal morphologies into the Honeybee Standard brain (HSB).

Supplement 3 movies in mpg format
a The spatial relationship of olfactory L3 neuron (blue), mechanosensory (DL and 2: green and red, respectively) and central interneuron Pe1 (yellow) after transformation to the Honeybee Standard Brain (HSB). b the same movie in stereo mode

Supplement 4 movie in wmv format
The Statistical Shape Model: Visualization of the central brain and morphological variations of the mushroom bodies and protocerebral neuropils (for more details: see text and Lienhard, 2008).

Supplement 5 movie in wmv format
Displacement model: Positioning and Displacement of the central brain model in the image data. Reddish-white flickering indicates the deformation of the model during the adjustment of the SSM to the image data.

aCknowledgmentS
We thank Alvar Prönneke for his help in evaluating the SSM Vincent Dercksen for landmark registration and Dirk Drenske for help with the 3D-PDF. We are particularly grateful to Gisela Manz and Astrid Klawitter for help with the reconstructions. We would also like to thank Daniel Münch and Hiro Ai for providing data for Figures 2 and 4. note An Interactive three-dimensional view of the Honeybee Standard Brain and integrated neurons for Figures 1,2A and 4A can be found in the 3D PDF file in the supplementary section. In order to utilises the 3D tools requires viewing with Adobe Acrobat Reader 8.0 or greater. gloSSary Labelfield There are two ways to represent segmentations of images: either boundaries are added to an image that enclose sets of pixels that are considered to belong together, or to pixels/voxels a label (a specific ID) is assigned that represents a particular object. In Amira the label field representation is used. Ontology In the information sciences, an ontology is a formal representation of concepts or structures and relations among those structures in a defined application domain. They represent semantic information using a controlled vocabulary. Visually, ontologies can be described as graphs in which structures are presented by nodes and relationships between structures are represented by edges. Registration The process of computing a coordinate transformation that maps the coordinates of one image onto the anatomically equivalent point in another image. Segmentation Classification of regions (intensity values) within the image data and partition of homogenous regions. The use of a single threshold means binarization of the image, i.e., separation of background and structure of interest. Segmentation may separate many structures within the image (connected components or segments of the image). Transformation Mathematical operation that moves, rotates, scales, and/or even deforms an object in such a way, that it will be aligned to another similar one. Warping Reformatting of an image under a deformation given by a non-rigid coordinate transformation. This is also known as also elastic or free-form deformation. Glossary on the Internet: http://en.wikipedia.org/wiki/ Image_ registration.